- phenixbb - phenix-online.org

SAD dataset
by sbiswas2＠ncsu.edu 28 Jun '10

28 Jun '10

Hi all, I was trying to run autosol with a sca file. I get the following error message: GUESS FILE TYPE MERGE TYPE sca premerged TARGET LABELS ['IPLUS', 'SIGIPLUS', 'IMINU', 'SIGIMINU'] Unit cell: (45.947, 88.031, 129.436, 90, 90, 90) Space group: F 2 2 2 (No. 22) CONTENTS: ['AutoSol_run_26_/TEMP0/F222peak_PHX_I_476(+).sca', 'sca', 'premerged', 'F 2 2 2', [45.947000000000003, 88.031000000000006, 129.43600000000001, 90.0, 90.0, 90.0], 2.000299364734313, ['IPLUS', 'SIGIPLUS', 'IMINU', 'SIGIMINU']] File F222peak_PHX_I_476(+).sca RES: 2.0 Guess of most complete dataset: 1 with resolution of 2.0 for dataset 1 FILE IS: AutoSol_run_26_/LAST.LOG ******************************************************************************** Failed to carry out AutoSol_scale_and_analyze_mad: None of the solve versions worked I guess this has something to do with the column label of sca file. How do I fix this? Any help is appreciated. thanks, Shya

1 0

Molecule will not load into coot from Phenix 1.6.2-432
by Subhangi Ghosh 25 Jun '10

25 Jun '10

Hi, I'm using Phenix version 1.6.2-432. When I run phenix.refine, coot opens up, but my molecule does not load. If I open in coot after refinement, molecule still does not load in coot. Please help. Thank you, Subhangi

2 6

NAG clashes: C2 with ND2
by Sam Stampfer 25 Jun '10

25 Jun '10

Hi, I'm in the process of preparing 4 structures for submission to rcsb, and the validation server has pointed out numerous issues with the N-acetyl-glucosamine moieties on my structures. While I have fixed the torsion angles accordingly (using CARP from http://www.glycosciences.de/tools/carp/ as a validation checker) and they stay fixed after phenix refinement, phenix keeps generating 2-2.2 Angstrom clashes between ND2 on asparagine with C2 on the NAG (the linkage is ND2-C1 and the usual distances are at least 2.28 Angstroms). Occasionally, the ND2 to C2 distance is as low as 1.8 Angstroms. This is happening in all my structures and happens independently of the setting of wxc_scale (even for wxc_scale=0.05!) My NAG groups fit the density well, both pre- and post- refinement. Any suggestions? Will rcsb care? These clashes are present in 12 out of 39 NAG moieties. Thanks! -Sam

3 3

PHENIX refine: detached sidechain atoms
by Oliver King 25 Jun '10

25 Jun '10

Hi All, I've noticed that when refining a model in phenix.refine using a PDB file from Refmac, the sidechain atoms of certain residues become detached from the rest of the molecule and appear to float on their own, at least when viewing in Coot. I think this is down to the format of the PDB file. For example the atoms of a Leu residue from a phenix.refine file which displays well is of the form: ATOM 1 N LEU A 8 50.022 -34.247 -5.817 1.00 58.12 N ATOM 2 CA LEU A 8 49.788 -34.905 -4.539 1.00 55.84 C ATOM 3 C LEU A 8 48.339 -35.358 -4.348 1.00 60.39 C ATOM 4 O LEU A 8 48.008 -36.016 -3.360 1.00 68.40 O ATOM 5 CB LEU A 8 50.219 -34.011 -3.373 1.00 57.98 C ATOM 6 CG LEU A 8 51.690 -34.126 -2.975 1.00 65.66 C ATOM 7 CD1 LEU A 8 52.014 -33.170 -1.836 1.00 59.99 C ATOM 8 CD2 LEU A 8 52.020 -35.566 -2.594 1.00 68.57 C where as from a Refmac PDB which also behaves well, it is of the form: ATOM 1 N LEU A 8 50.453 -35.722 -5.617 1.00 20.00 N ATOM 2 CA LEU A 8 49.649 -35.131 -4.482 1.00 20.00 C ATOM 3 CB LEU A 8 50.190 -33.735 -4.147 1.00 20.00 C ATOM 4 CG LEU A 8 51.461 -33.755 -3.275 1.00 20.00 C ATOM 5 CD1 LEU A 8 52.556 -32.768 -3.741 1.00 20.00 C ATOM 6 CD2 LEU A 8 51.082 -33.546 -1.799 1.00 20.00 C ATOM 7 C LEU A 8 48.166 -35.063 -4.824 1.00 20.00 C ATOM 8 O LEU A 8 47.822 -34.257 -5.568 1.00 20.00 O but after refinement in phenix it becomes ATOM 1 N LEU A 8 50.734 -35.936 -5.935 1.00 54.45 N ATOM 2 CA LEU A 8 49.846 -35.185 -4.892 1.00 51.89 C ATOM 3 CB LEU A 8 50.377 -33.761 -4.516 1.00 51.22 C ATOM 4 CG LEU A 8 51.550 -33.730 -3.431 1.00 59.27 C ATOM 5 CD1 LEU A 8 52.457 -32.407 -3.635 1.00 60.88 C ATOM 6 CD2 LEU A 8 51.120 -34.090 -1.879 1.00 55.50 C ATOM 7 C LEU A 8 48.268 -35.042 -5.221 1.00 47.92 C ATOM 8 O LEU A 8 47.989 -34.282 -6.195 1.00 47.32 O and the CD2 atom is now too far away to be part of the residue in Coot Is there an easy way to convert a PDB from Refmac into one which will behave itself when put through phenix.refine? I've tried using phenix.pdb_tools and also putting the model through Molprobity and hoping that that the output would be corrected. Thanks, Olly King

6 8

Defining planarity for Phenix.refine
by Peter Grey 22 Jun '10

22 Jun '10

Dear all, I have an RNA/Protein big complex at low resolution (roughly 4A). I would like to have the base pairs in the RNA as close as possible to ideal base pairs. I added restraints for distances between base pair hydrogen-bonding atoms but this is not enough to ensure that the bases will be in the same plane. Could you suggest how to define for Phenix.refine this planarity ? I am grateful for your advice, Peter

7 16

Re: [phenixbb] AutoMR:Multi-subunit model
by S. Saif Hasan 22 Jun '10

22 Jun '10

Hello, Thanks. When I use the PDB structure with 8 polypeptides as the search model and specify the 8 sequences together, I get this warning from AutoMR- "ASYMMETRIC UNIT TOO FULL (scattering from ensembles =1.82 of total scattering)- TOTAL SCATTERING INCREASED: To get rid of this warning message do one of the following: (1) Increasing the scattering in the ASU (keyword COMPOSITION) because you have more scattering than you thought (2) Decrease the number of ensembles you are packing into the ASU (keywords SOLUTION and SEARCH)" Could this be a problem in missing some additional prosthetic groups ? I did not specify any non-protein groups but the search model has them. -- S. Saif Hasan Graduate Student Department of Biological Sciences Hockmeyer Hall of Structural Biology Purdue University 240 S. Martin Jischke Drive West Lafayette, IN 47907 Quoting "Thomas C. Terwilliger" <terwilliger(a)lanl.gov>: > Hi Saif, > > If I understand you correctly, you can specify the search model with 8 > polypeptides as one PDB file and all 8 sequences (separated by blank lines > or lines starting with ">" in one sequence file. Also component 1 can be > all 8 polypeptides and copies=1. > > All the best, > Tom T > > >> Hello, > >> > >> I have a membrane protein complex made up of 8 different polypeptides > >> present in > >> two copies each (16mer). The asymmetric unit is the octamer, as > determined > >> from > >> previous crystallographic studies. To carry out molecular replacement > with > >> the > >> phases from the published structure in AutoMR, should I specify the > >> coordinates > >> of the search model and its overall sequence identity as one file or > >> should I > >> break it down for each subunit ? Also, for the asymmetric unit contents, > >> should > >> I specify the sequence or mass of each subunit separately (as component > 1, > >> 2 and > >> so on) keeping the number of copies fixed at 1 ? There is no NCS. > >> > >> Thanks. > >> -- > >> S. Saif Hasan > >> > >> Graduate Student > >> Department of Biological Sciences > >> Hockmeyer Hall of Structural Biology > >> Purdue University > >> 240 S. Martin Jischke Drive > >> West Lafayette, IN 47907 > >> > >> > >> > >> > >> > >> _______________________________________________ > >> phenixbb mailing list > >> phenixbb(a)phenix-online.org > >> http://phenix-online.org/mailman/listinfo/phenixbb > >> > > _______________________________________________ > phenixbb mailing list > phenixbb(a)phenix-online.org > http://phenix-online.org/mailman/listinfo/phenixbb >

1 0

AutoMR:Multi-subunit model
by S. Saif Hasan 21 Jun '10

21 Jun '10

Hello, I have a membrane protein complex made up of 8 different polypeptides present in two copies each (16mer). The asymmetric unit is the octamer, as determined from previous crystallographic studies. To carry out molecular replacement with the phases from the published structure in AutoMR, should I specify the coordinates of the search model and its overall sequence identity as one file or should I break it down for each subunit ? Also, for the asymmetric unit contents, should I specify the sequence or mass of each subunit separately (as component 1, 2 and so on) keeping the number of copies fixed at 1 ? There is no NCS. Thanks. -- S. Saif Hasan Graduate Student Department of Biological Sciences Hockmeyer Hall of Structural Biology Purdue University 240 S. Martin Jischke Drive West Lafayette, IN 47907

3 2

secondary restraints
by Jianghai Zhu 21 Jun '10

21 Jun '10

Hi, I am using secondary restraints in the phenix.refine. I got the following error message. hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'A' and resseq 141:143) and (altloc 'A' or altloc ' ') and name H" => 0 donors "(chain 'A' and resseq 141:143) and (altloc 'A' or altloc ' ') and name O" => 3 acceptors hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'A' and resseq 160:162) and (altloc 'A' or altloc ' ') and name H" => 0 donors "(chain 'A' and resseq 160:162) and (altloc 'A' or altloc ' ') and name O" => 3 acceptors hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'A' and resseq 242:244) and (altloc 'A' or altloc ' ') and name H" => 0 donors "(chain 'A' and resseq 242:244) and (altloc 'A' or altloc ' ') and name O" => 3 acceptors hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'A' and resseq 360:363) and (altloc 'A' or altloc ' ') and name H" => 0 donors "(chain 'A' and resseq 360:363) and (altloc 'A' or altloc ' ') and name O" => 4 acceptors hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'A' and resseq 387:389) and (altloc 'A' or altloc ' ') and name H" => 0 donors "(chain 'A' and resseq 387:389) and (altloc 'A' or altloc ' ') and name O" => 3 acceptors hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'A' and resseq 455:457) and (altloc 'A' or altloc ' ') and name H" => 0 donors "(chain 'A' and resseq 455:457) and (altloc 'A' or altloc ' ') and name O" => 3 acceptors hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'B' and resseq 143:152) and (altloc 'A' or altloc ' ') and name H" => 0 donors "(chain 'B' and resseq 143:152) and (altloc 'A' or altloc ' ') and name O" => 10 acceptors hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'B' and resseq 156:165) and (altloc 'A' or altloc ' ') and name H" => 0 donors "(chain 'B' and resseq 156:165) and (altloc 'A' or altloc ' ') and name O" => 10 acceptors hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'B' and resseq 190:193) and (altloc 'A' or altloc ' ') and name H" => 1 donors "(chain 'B' and resseq 190:193) and (altloc 'A' or altloc ' ') and name O" => 4 acceptors hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'B' and resseq 222:227) and (altloc 'A' or altloc ' ') and name H" => 0 donors "(chain 'B' and resseq 222:227) and (altloc 'A' or altloc ' ') and name O" => 6 acceptors hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'B' and resseq 243:251) and (altloc 'A' or altloc ' ') and name H" => 0 donors "(chain 'B' and resseq 243:251) and (altloc 'A' or altloc ' ') and name O" => 9 acceptors hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'B' and resseq 253:256) and (altloc 'A' or altloc ' ') and name H" => 0 donors "(chain 'B' and resseq 253:256) and (altloc 'A' or altloc ' ') and name O" => 4 acceptors hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'B' and resseq 278:282) and (altloc 'A' or altloc ' ') and name H" => 0 donors "(chain 'B' and resseq 278:282) and (altloc 'A' or altloc ' ') and name O" => 5 acceptors hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'B' and resseq 301:305) and (altloc 'A' or altloc ' ') and name H" => 0 donors "(chain 'B' and resseq 301:305) and (altloc 'A' or altloc ' ') and name O" => 5 acceptors hydrogen_bonds_from_selections: incomplete non-PRO residues in helix. "(chain 'B' and resseq 311:321) and (altloc 'A' or altloc ' ') and name H" => 0 donors "(chain 'B' and resseq 311:321) and (altloc 'A' or altloc ' ') and name O" => 11 acceptors Traceback (most recent call last): File "/nfs/home/jzhu/phenix/phenix-1.6.2-432/phenix/phenix/command_line/refine.py", line 11, in <module> command_line.run(command_name="phenix.refine", args=sys.argv[1:]) File "/nfs/home/jzhu/phenix/phenix-1.6.2-432/phenix/phenix/refinement/command_line.py", line 71, in run overwrite = command_line_interpreter.command_line.options.overwrite) File "/nfs/home/jzhu/phenix/phenix-1.6.2-432/phenix/phenix/refinement/driver.py", line 121, in __init__ h_bond_table = self.sec_str.get_bonds_table(log=log) File "/nfs/home/jzhu/phenix/phenix-1.6.2-432/cctbx_project/mmtbx/secondary_structure.py", line 1023, in get_bonds_table log=log) File "/nfs/home/jzhu/phenix/phenix-1.6.2-432/cctbx_project/mmtbx/secondary_structure.py", line 444, in hydrogen_bonds_from_selections ss_type="helix") File "/nfs/home/jzhu/phenix/phenix-1.6.2-432/cctbx_project/mmtbx/secondary_structure.py", line 652, in _donors_and_acceptors donor_isel.insert(k, n_atoms) IndexError: Index out of range. How do I fix it? Thanks. -- Jianghai

2 1

how to extend R-free set
by Robert Immormino 21 Jun '10

21 Jun '10

Hi, Hopefully I'm just missing something obvious.... I'm trying to copy the Rfree set from an old set of data, to a new mutant structure....and extend to higher resolution. I've tried: phenix.reflection_file_converter old.mtz --label="R-free-flags" new.sca --label="i_obs,sigma" --mtz new.mtz and several other combination, but I am still failing to get phenix to read multiple data files and output the desired mtz. I can get just the Rfree flags or the I and SigI, but not both. I also tried with phenix.refine phenix.refine --dry-run old.pdb old.mtz new.sca ligands.cif refinement.input.xray_data.labels="i_obs,sigma" refinement.input.xray_data.r_free_flags.label="R-free-flags" This seems closer, but I get the error: Sorry: R-free flags not compatible with F-obs array: missing flag for 492 F-obs selected for refinement: Thanks, -bob

7 11

find_all_ligands vs. ligand_identification
by Mark A Saper 21 Jun '10

21 Jun '10

Could someone explain the difference between the GUI "ligand_identification" routine and the documented command line "find_all_ligands". The latter appears to have more flexible input parameters. I'm not sure the GUI for ligand_identification is working correctly. Thanks, Mark _________________________________ Mark A. Saper, Ph.D. Associate Professor of Biological Chemistry University of Michigan Biophysics, 3040 Chemistry Building 930 N University Ave Ann Arbor MI 48109-1055 U.S.A. saper(a)umich.edu phone (734) 764-3353 fax (734) 764-3323 http://www.biochem.med.umich.edu/?q=saper http://www.strucbio.umich.edu/

2 1