- phenixbb - phenix-online.org

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 23 Jun '10

23 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 22 Jun '10

22 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

LAST CHANCE TO REGISTER FOR DIFFRACTION METHODS GRC
by A Leslie 21 Jun '10

21 Jun '10

The biannual Gordon Conference on Diffraction Methods in Structural Biology will be held at Bates College in Lewiston, Maine from 18-23rd July. The program covers all aspects of macromolecular structure solution from crystallization to structure solution and refinement and includes sessions on complementary techniques and the potential of free electron lasers. Many of the leading methods developers will be present and the format of the meeting, with free afternoons, provides ample opportunities for informal discussions. The registration deadline is *** THIS COMING SUNDAY, 27th June *** = ======================================================================== Limited funds are still available to support new student and postdoctoral applications, contact andrew(a)mrc-lmb.cam.ac.uk for details. = ======================================================================== Further details, including a full program and details on the application procedure can be obtained from: http://www.grc.org/programs.aspx?year=2010&program=diffrac We hope to see you there ! Andrew Leslie (Chair) & Ana Gonzalez (Vice Chair) Outline Program Diffraction Methods in Structural Biology, 18-23rd July 2010, Bates College, Lewiston, Maine Sunday pm: Macromolecular Structures, Pushing the boundaries Discussion leader: Michael Rossmann 1. ABC Transporters: Doug Rees 2. The Spliceosome : Kyoshi Nagai Monday am: Crystallisation and the diffraction experiment Discussion leader: Bob Sweet 1. Nanolitre crystallisation: Seth Harris 2. Twinning and disorder: Todd Yeates 3. Simulating the diffraction experiment: James Holton 4. Microcrystallography: Gwyndaf Evans Monday pm: Radiation damage and data collection strategies Discussion leader: Vivian Stojanoff 1. Radiation damage, experimental: Martin Weik 2. Microbeams: Yanhui Zou 3. Testing data collection strategies with the Pilatus: Marcus Mueller Tues am: Synchrotron developments and automated data processing Discussion leader: Thomas Schneider 1. Synchrotrons, latest developments: Sean Mc Sweeney 2. Synchrotrons, latest developments: Janet Smith 3. Data processing with Xia2/CHEF: Graeme Winter 4. Data processing with AUTOPROC: Clemens Vonrhein Tues pm: Complementary techniques Discussion leader: Ana Gonzales 1. Electron tomography: Sriram Subramaniam 2. SAXS: Hiro Tsuruta 3. Spectroscopy: Carrie Wilmot Weds am: Structure solution and refinement Discussion leader: Paul Adams 1. SAD Phasing: Randy Read 2. Molecular replacement with BALBES: Garib Murshudov 3. PHENIX: Tom Terwilliger 4. TLS: Ethan Merritt Weds pm: Challenging problems/Membrane Proteins Discussion leader: Michael Wiener 1. Using cubic lipidic phase: Martin Caffrey 2. Membrane Proteins: Stephen White 3. Heterogeneity, low solubility and low resolution: structural pursuit on the ternary Myddosome complex in Toll-like receptor signaling: Hao Wu Thurs am: Selected Posters and Map improvement and Model building Discussion leader: Tassos Perrakis 1. COOT: Paul Emsley 2. Modelling nucleic acids: Johan Hattne Thursday pm: Future Methods, FELs, Diffraction imaging Discussion leader: Ed Lattman 1. Opportunities fof membrane protein analysis with FELs: Petra Fromme 2. Single particle imaging with pulsed hard Xray lasers: Anton Barty 3. Femtosecond nanocrystallography of membrane proteins using LCSL: John Spence

1 0

"Open in Coot" not working on 1.6.1-357
by Ben Eisenbraun 19 Jun '10

19 Jun '10

Howdy PHENIXians, I can't seem to get the "Open in Coot" option to work in PHENIX. I tested version 1.6.1-357 on 32 and 64-bit linux as well as the OS X Intel version, and while Coot will open, it won't display the working data. 'phenix.find_coot_command' works fine from the command line. On the linux machines, when Coot opens, it doesn't have the "Connected to PHENIX" button. On the OS X Intel machine, it does have the button, but it still won't display the working data. Any ideas? -ben -- | Ben Eisenbraun | Software Sysadmin | | Structural Biology Grid | http://sbgrid.org | | Harvard Medical School | http://hms.harvard.edu |

6 25

is derivative data useful for looking for the missingpart of the model?
by Jacob Keller 16 Jun '10

16 Jun '10

> Hi Fengyun, > > How have you determined whether you in fact have the full-length protein > in your crystals? It sounds like, from your statistics, that everything is > already modelled--it seems hard to believe that 80% of the whole protein > would give such good R values. Do others agree? > > You could run some of your crystals on a gel next to a sample which is > known to be full-length, and you could get the answer. > > Jacob > > ******************************************* > Jacob Pearson Keller > Northwestern University > Medical Scientist Training Program > Dallos Laboratory > F. Searle 1-240 > 2240 Campus Drive > Evanston IL 60208 > lab: 847.491.2438 > cel: 773.608.9185 > email: j-keller2(a)northwestern.edu > ******************************************* > > ----- Original Message ----- > From: <fn1(a)rice.edu> > To: <phenixbb(a)phenix-online.org> > Sent: Tuesday, June 15, 2010 5:59 PM > Subject: [phenixbb] is derivative data useful for looking for the > missingpart of the model? > > >> Hi everyone, >> >> I got a partial model, consisting of 80% of the whole molecule, from the >> molecular replacement. The native data reflects to 2.3 A, and now the >> R/Rfree= 0.22/0.28. My problem is that the missing part of the model is >> important for us, but there is not any density for this region now. >> >> Meantime, we have several heavy atom derivative data. My question is >> that if refine the model against the derivative data, will the heavy >> atom appear in the difference map if existing? Any suggestions on how to >> use the derivative data to locate the missing part of model? >> >> Thanks in advance! >> Regards, >> Fengyun >> >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://phenix-online.org/mailman/listinfo/phenixbb >

2 1