- phenixbb - phenix-online.org

ramachandran_restraints
by José Trincão 28 Dec '10

28 Dec '10

Dear all, I am refining a structure to 2.9 A, with pretty bad experimental data (P1, ~5400 residues in the A.U. (4 mol), anisotropic, fairly incomplete...). I have been running phenix.refine with NCS, 1 ADP per residue, lower target weights, TLS. My R/Rfree are stuck at about 2.4/3.05. I was getting around 2% ramachandran outliers at every cycle so I have been giving the nightly builds a shot. Now I get 0.2% outliers without changing any parameters for the ramachandra_restraints but it is forcing some of the peptide bonds out of the density. What can I tweak to improve the refinement? Thanks all in advance! Jose José Trincão, PhD CQFB@FCT-UNL 2829-516 Caparica, Portugal "It's very hard to make predictions... especially about the future" - Niels Bohr

3 2

Re: [phenixbb] Simple Phaser Question (Ed Pozharski
by Joseph Noel 27 Dec '10

27 Dec '10

Hi Ed, Yes, that is what I generally do but it sometimes is a slightly different so I wanted a more general approach. Joe ___________________________________________________________ Joseph P. Noel, Ph.D. Investigator, Howard Hughes Medical Institute Professor, The Jack H. Skirball Center for Chemical Biology and Proteomics The Salk Institute for Biological Studies 10010 North Torrey Pines Road La Jolla, CA 92037 USA Phone: (858) 453-4100 extension 1442 Cell: (858) 349-4700 Fax: (858) 597-0855 E-mail: noel(a)salk.edu Web Site (Salk): http://www.salk.edu/faculty/faculty_details.php?id=37 Web Site (HHMI): http://hhmi.org/research/investigators/noel.html ___________________________________________________________ On Dec 25, 2010, at 12:00 PM, phenixbb-request(a)phenix-online.org wrote: > Send phenixbb mailing list submissions to > phenixbb(a)phenix-online.org > > To subscribe or unsubscribe via the World Wide Web, visit > http://phenix-online.org/mailman/listinfo/phenixbb > or, via email, send a message with subject or body 'help' to > phenixbb-request(a)phenix-online.org > > You can reach the person managing the list at > phenixbb-owner(a)phenix-online.org > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of phenixbb digest..." > > > Today's Topics: > > 1. Simple Phaser Question (Joseph Noel) > 2. Re: Simple Phaser Question (Ed Pozharski) > 3. Questions about phenix.refine with twin_law (Keitaro Yamashita) > 4. Re: Questions about phenix.refine with twin_law (Peter Zwart) > 5. Re: Questions about phenix.refine with twin_law (Pavel Afonine) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Fri, 24 Dec 2010 17:11:37 -0800 > From: Joseph Noel <noel(a)salk.edu> > To: phenixbb(a)phenix-online.org > Subject: [phenixbb] Simple Phaser Question > Message-ID: <5E3B81ED-4C13-488E-8C28-33DB29263AEF(a)salk.edu> > Content-Type: text/plain; charset="us-ascii" > > Hi All, > > I am working with a protein that is a physiological dimer. On occasion I obtain a unit cell with one monomer in the asymmetric unit but most often I obtain an enantiomorphic space group with an expanded C axis that contains two molecules in the asymmetric unit. Anyway, that is all moot to some extent. I am wondering if for the case of my NCS dimer, is there a way using Phaser in Phenix to ensure that when both monomers are found, they are "close together" representative of the physiological dimer? When I search with individual monomers I often obtain a solution with two monomers but in different asymmetric units. Of course, I could easily do the transformation to form the physiological NCS dimer for later refinement but just wondering if I could avoid this from the get go with the appropriate Phaser keyword input (GUI if possible). > > Thanks and Happy Holidays! > > Joe > ___________________________________________________________ > Joseph P. Noel, Ph.D. > Investigator, Howard Hughes Medical Institute > Professor, The Jack H. Skirball Center for Chemical Biology and Proteomics > The Salk Institute for Biological Studies > 10010 North Torrey Pines Road > La Jolla, CA 92037 USA > > Phone: (858) 453-4100 extension 1442 > Cell: (858) 349-4700 > Fax: (858) 597-0855 > E-mail: noel(a)salk.edu > > Web Site (Salk): http://www.salk.edu/faculty/faculty_details.php?id=37 > Web Site (HHMI): http://hhmi.org/research/investigators/noel.html > ___________________________________________________________ > >

3 2

Simple Phaser Question
by Joseph Noel 25 Dec '10

25 Dec '10

Hi All, I am working with a protein that is a physiological dimer. On occasion I obtain a unit cell with one monomer in the asymmetric unit but most often I obtain an enantiomorphic space group with an expanded C axis that contains two molecules in the asymmetric unit. Anyway, that is all moot to some extent. I am wondering if for the case of my NCS dimer, is there a way using Phaser in Phenix to ensure that when both monomers are found, they are "close together" representative of the physiological dimer? When I search with individual monomers I often obtain a solution with two monomers but in different asymmetric units. Of course, I could easily do the transformation to form the physiological NCS dimer for later refinement but just wondering if I could avoid this from the get go with the appropriate Phaser keyword input (GUI if possible). Thanks and Happy Holidays! Joe ___________________________________________________________ Joseph P. Noel, Ph.D. Investigator, Howard Hughes Medical Institute Professor, The Jack H. Skirball Center for Chemical Biology and Proteomics The Salk Institute for Biological Studies 10010 North Torrey Pines Road La Jolla, CA 92037 USA Phone: (858) 453-4100 extension 1442 Cell: (858) 349-4700 Fax: (858) 597-0855 E-mail: noel(a)salk.edu Web Site (Salk): http://www.salk.edu/faculty/faculty_details.php?id=37 Web Site (HHMI): http://hhmi.org/research/investigators/noel.html ___________________________________________________________

2 1

Non-standard residue numbers in NCS restraints
by Sampson, Jared 23 Dec '10

23 Dec '10

Greetings all, My lab is working on a 3.2Å resolution structure with 6 NCS copies of an antibody Fab fragment with insertions in its CDR loops. The residue numbers of CDR loop H3, for example, include 100A, 100B, etc. in each of the heavy chains (Kabat numbering). There are observable differences between NCS copies in parts of the model due to crystal contacts and so forth, so we want to exclude those residues from the NCS. Phenix.refine, when given NCS restraints like the example below, produces a result with much higher R/Rfree values than when the same structure was run with "normal" continuous numbering (R/Rfree of .26/.30 compared to .21/.26 before changing to Kabat numbering). Reference selection: chain H and (resseq 1:100A or resseq 100H:113) Restrained selection: chain I and (resseq 1:100A or resseq 100H:113) Phenix version is 1.6.4-486. Any thoughts on what might be causing this or how to account for the non-standard numbers in NCS restraints? Thanks, Jared Sampson -- Jared Sampson Xiangpeng Kong Lab NYU Langone Medical Center New York, NY 10016 212-263-7898 ------------------------------------------------------------ This email message, including any attachments, is for the sole use of the intended recipient(s) and may contain information that is proprietary, confidential, and exempt from disclosure under applicable law. Any unauthorized review, use, disclosure, or distribution is prohibited. If you have received this email in error please notify the sender by return email and delete the original message. Please note, the recipient should check this email and any attachments for the presence of viruses. The organization accepts no liability for any damage caused by any virus transmitted by this email. =================================

2 1

Problem with Molprobity analysis in phenix 1.6.4
by Mark A Saper 23 Dec '10

23 Dec '10

I am running phenix 1.6.4-486 with GUI and obtained high clash scores in the Molprobity calculations done by phenix.refine. I think I know why after doing some experiments with the current Molprobity running at http://molprobity.biochem.duke.edu/. If I take my well-refined structure as output by phenix (.pdb) and use it directly on the Molprobity web site, the clash score is much lower: 2.75 for 39 total bumps. Many of these involve residues or waters with high B factors. If I first remove all hydrogens from the phenix-refined .pdb file with a text editor, and upload it to the Molprobity web site, then add hydrogens, and then calculate clash score, the score is now 8.69 or 123 bumps for 14150 atoms. This is *identical* to what phenix calculates. Therefore, why doesn't phenix feed the refined coordinates directly into the Molprobity (probe) calculations? Has this change been made in one of the nightly builds? 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/

6 7

Fwd: Fourth annual CCP4 summer school in USA, at APS, June 7-15
by Nathaniel Echols 22 Dec '10

22 Dec '10

---------- Forwarded message ---------- From: "Sanishvili, Ruslan" <rsanishvili(a)anl.gov> To: <phenixbb(a)phenix-online.org> Date: Wed, 22 Dec 2010 11:02:15 -0600 Subject: Fourth annual CCP4 summer school in USA, at APS, June 7-15 Dear Colleagues, Building on the success of the past 3 years, we are pleased to announce the fourth annual CCP4 summer school at Advanced Photon Source (APS), Argonne National Laboratory (ANL). All details can be found at http://www.ccp4.ac.uk/schools/APS-2011/index.php Title: "CCP4 school: From data collection to structure refinement and beyond" Dates: June 7 to 15. Site: Advanced Photon Source, Argonne National Laboratory, Argonne, Illinois (Near Chicago), USA The school content: Data collection workshop the first two days: beamline training and data collection on GM/CA-CAT beamlines 23ID-B and 23ID-D. For data collection, only the participants' crystals will be used. Software workshop: The rest of the time after data collection will feature many modern crystallographic software packages taught by authors and other experts. It will be organized in three Sections – lectures, tutorials and hands-on trouble-shooting. There will be model data sets available for tutorials but data, provided by participants, will have higher priority for the hands-on sessions. Applicants: Graduate students, postdoctoral researchers and young scientists at the assistant professor level are encouraged to apply. Only 20 applicants will be selected for participation. Participants of the workshop are strongly encouraged to bring their own problem data sets or crystals so the problems can be addressed during data collection workshop and/or hands-on sessions. Application: Application deadline is April 8. Application form, the program, contact info and other details can be found at http://www.ccp4.ac.uk/schools/APS-2011/index.php Fees: There is no fee for the workshop. The students will be responsible for their transportation and lodging. The workshop organizers will arrange economical lodging at the Argonne Guest House. The workshop will also cover the expenses for all meals and refreshments. Garib, Ronan and Nukri Ruslan Sanishvili (Nukri), Ph.D. GM/CA-CAT Biosciences Division, ANL 9700 S. Cass Ave. Argonne, IL 60439 Tel: (630)252-0665 Fax: (630)252-0667 rsanishvili(a)anl.gov

1 0

ADP restraints - distance power, average power, etc
by Joseph Noel 21 Dec '10

21 Dec '10

Hi All, I am refining a very well ordered structure at 1.45 A and find that after doing just about everything including optimization of weights, etc that there are still areas of positive density residing over S atoms of Met, backbone atoms, etc. I have added hydrogens and used them during refinement as well. The maps are of a high quality and Free R factors are quite good, ~ 17%. Is there anything that can be played with more such as values used in the ADP restraints window to try and achieve a difference map with far fewer areas of significant positive density (all greater then or equal to 4 sigma)? I am not sure exactly what effect the values of distance power, average power, etc will have on refinement. The other option I was thinking of was the scattering table options. I have other structures of this protein that extend to about 1.2. Thanks! Joe ___________________________________________________________ Joseph P. Noel, Ph.D. Investigator, Howard Hughes Medical Institute Professor, The Jack H. Skirball Center for Chemical Biology and Proteomics The Salk Institute for Biological Studies 10010 North Torrey Pines Road La Jolla, CA 92037 USA Phone: (858) 453-4100 extension 1442 Cell: (858) 349-4700 Fax: (858) 597-0855 E-mail: noel(a)salk.edu Web Site (Salk): http://www.salk.edu/faculty/faculty_details.php?id=37 Web Site (HHMI): http://hhmi.org/research/investigators/noel.html ___________________________________________________________

7 16

Re: [phenixbb] ADP restraints - distance power, average power, etc
by Joseph Noel 20 Dec '10

20 Dec '10

Dear Ed, Thanks so much for the suggestion about REFMAC but I seemed to have found out the issue. The positive density on Met sulfurs and backbone carbonyls and amide NHs does not appear when I calculate a "normal" sigmaa weighted mFo-dFc map. The positive density was appearing only when I checked the kicked map option. So the "extra" density seems to be due to generating kicked maps. Sorry about the confusion related to the sigma level. Just meant that the positive density was not insignificant in the kicked mFo-dFc maps. Density was visible on these specific areas starting around 3.5-4.0 sigma and visible in several places contoured up to 8 sigma. Joe ___________________________________________________________ Joseph P. Noel, Ph.D. Investigator, Howard Hughes Medical Institute Professor, The Jack H. Skirball Center for Chemical Biology and Proteomics The Salk Institute for Biological Studies 10010 North Torrey Pines Road La Jolla, CA 92037 USA Phone: (858) 453-4100 extension 1442 Cell: (858) 349-4700 Fax: (858) 597-0855 E-mail: noel(a)salk.edu Web Site (Salk): http://www.salk.edu/faculty/faculty_details.php?id=37 Web Site (HHMI): http://hhmi.org/research/investigators/noel.html ___________________________________________________________ On Dec 20, 2010, at 12:00 PM, phenixbb-request(a)phenix-online.org wrote: > Send phenixbb mailing list submissions to > phenixbb(a)phenix-online.org > > To subscribe or unsubscribe via the World Wide Web, visit > http://phenix-online.org/mailman/listinfo/phenixbb > or, via email, send a message with subject or body 'help' to > phenixbb-request(a)phenix-online.org > > You can reach the person managing the list at > phenixbb-owner(a)phenix-online.org > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of phenixbb digest..." > > > Today's Topics: > > 1. ADP restraints - distance power, average power, etc (Joseph Noel) > 2. Re: ADP restraints - distance power, average power, etc > (Ralf W. Grosse-Kunstleve) > 3. Re: ADP restraints - distance power, average power, etc > (Pavel Afonine) > 4. Secondary structure restraints - where to specify file using > GUI? (Luke Rice) > 5. Re: Secondary structure restraints - where to specify file > using GUI? (Nathaniel Echols) > 6. Re: ADP restraints - distance power, average power, etc > (Joseph Noel) > 7. Re: ADP restraints - distance power, average power, etc > (Joseph Noel) > 8. Re: ADP restraints - distance power, average power, etc > (Ed Pozharski) > 9. Re: Secondary structure restraints - where to specify file > using GUI? (Luke Rice) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Sun, 19 Dec 2010 12:58:43 -0800 > From: Joseph Noel <noel(a)salk.edu> > To: phenixbb(a)phenix-online.org > Subject: [phenixbb] ADP restraints - distance power, average power, > etc > Message-ID: <2E132727-1B34-4ED5-98DA-2210B7D1FF48(a)salk.edu> > Content-Type: text/plain; charset="us-ascii" > > Hi All, > > I am refining a very well ordered structure at 1.45 A and find that after doing just about everything including optimization of weights, etc that there are still areas of positive density residing over S atoms of Met, backbone atoms, etc. I have added hydrogens and used them during refinement as well. The maps are of a high quality and Free R factors are quite good, ~ 17%. Is there anything that can be played with more such as values used in the ADP restraints window to try and achieve a difference map with far fewer areas of significant positive density (all greater then or equal to 4 sigma)? > > I am not sure exactly what effect the values of distance power, average power, etc will have on refinement. The other option I was thinking of was the scattering table options. I have other structures of this protein that extend to about 1.2. > > Thanks! > > Joe > ___________________________________________________________ > Joseph P. Noel, Ph.D. > Investigator, Howard Hughes Medical Institute > Professor, The Jack H. Skirball Center for Chemical Biology and Proteomics > The Salk Institute for Biological Studies > 10010 North Torrey Pines Road > La Jolla, CA 92037 USA > > Phone: (858) 453-4100 extension 1442 > Cell: (858) 349-4700 > Fax: (858) 597-0855 > E-mail: noel(a)salk.edu > > Web Site (Salk): http://www.salk.edu/faculty/faculty_details.php?id=37 > Web Site (HHMI): http://hhmi.org/research/investigators/noel.html > ___________________________________________________________ > >

1 0

phenix.refine with simulated annealing
by zhangh1＠umbc.edu 18 Dec '10

18 Dec '10

Hi there, I was trying to perform a SA refinement for PDB# 2VZ8 using the following phenix script: random_seed=2679941 phenix.refine 2VZ8.pdb 2vz8.mtz simulated_annealing=true main.random_seed=${random_seed} output.prefix=${random_seed} >${PDB}-phx-sa.log where 2VZ8.pdb was downloaded from rcsb, and 2vz8.mtz was converted from the cif file by ccp4-cif2mtz. But it ended up with: ******************* =============================== refinement start ============================== Traceback (most recent call last): File "/users/hz9/TOOLS/PHENIX/phenix-1.6.4-486/phenix/phenix/command_line/refine.py", line 11, in <module> command_line.run(command_name="phenix.refine", args=sys.argv[1:]) File "/users/hz9/TOOLS/PHENIX/phenix-1.6.4-486/phenix/phenix/refinement/command_line.py", line 89, in run call_back_handler=call_back_handler) File "/users/hz9/TOOLS/PHENIX/phenix-1.6.4-486/phenix/phenix/refinement/driver.py", line 1182, in run call_back_handler = call_back_handler) File "/users/hz9/TOOLS/PHENIX/phenix-1.6.4-486/phenix/phenix/refinement/strategies.py", line 229, in refinement_machine xray_structure = model.xray_structure,) File "/users/hz9/TOOLS/PHENIX/phenix-1.6.4-486/phenix/phenix/refinement/strategies.py", line 133, in extract_from_pdb_header xray_structure = xray_structure) File "/users/hz9/TOOLS/PHENIX/phenix-1.6.4-486/cctbx_project/mmtbx/f_model.py", line 302, in __init__ twin_law_xyz = sgtbx.rt_mx(symbol=self.twin_law, r_den=12, t_den=144) RuntimeError: cctbx Error: Parse error: unexpected character: NULL ************************ I tried another smaller system and it works well. Not sure what happens, and thanks for any suggestion! Hailiang

4 3

phenix.map_to_object for HA sites
by Michael Hothorn 17 Dec '10

17 Dec '10

Hi, I have a .pdb with heavy atoms that are dispersed all over the unit cell. I would like to move all of them (about 40) as close as possible to my refined structure. Could phenix.map_to_object do that for all heavy atoms at once, or would I need to do this in an atom by atom manner? Thanks Michael

2 1