- phenixbb - phenix-online.org

on TLS for Phenix refine
by Smith Lee 22 Feb '15

22 Feb '15

Dear All, For TLS, the document says, "Can I use both TLS and anisotropic ADPs? Yes, but not for the same atoms - since TLS is essentially constrained anisotropic refinement, the two methods are mutually exclusive". In the "Refinment settings" of the Phenix refine graphical interface, will you please tell me which button is for the anisotropic ADP? The individual B-factors, or the Group B-factors, or only the "Optimize X-ray/ADP weight"? I am looking forward to getting your reply. Smith

3 2

Odd behaviour of real-space refinement
by Tristan Croll 22 Feb '15

22 Feb '15

Hi all, I'm just finalising the re-refinement of a 3.6 Angstrom structure. The stats in general are excellent (Rwork/Rfree = 0.219/0.260; Ramachandran outliers 0.12%; Clashscore < 2; rotamer outliers 2.1%) and I've gone through the whole structure to check and double-check fit to the density. I'm very confident that it's effectively correct. When I run through a fairly standard refinement protocol (XYZ coords, Individual B-factors, optimizing weights, secondary structure restraints), all seems fine - the Rfree rattles around within 0.2% or so. But when I include the real-space refinement option, very strange things start to happen. Each rsrl step causes a 3-4% jump in both Rwork and Rfree, and jumps the Clashscore up to 40 or so. Rotamer outliers briefly go to zero but then increase substantially (to 6+%) in subsequent steps, and the number of Ramachandran outliers increases. I thought I'd missed some dreadful error in my structure until I looked more closely at the results. Some very nonsensical changes have been made - for example, an arginine in the middle of a beta-sheet (admittedly with no observable density associated with its sidechain, but register is well determined by surrounding residues and backbone density is clear) was flipped from its original orientation to point directly into the hydrophobic core of the domain. Somewhat surprising, to say the least. Is this a bug in the algorithm, or is it simply that 3.6A is too low resolution for it to work reliably? My version of Phenix is a nightly build from about 2 weeks ago. Happy to share coordinates and log files if necessary. Thanks, Tristan

1 0

Odd behaviour of real-space refinement
by Tristan Croll 21 Feb '15

21 Feb '15

Hi all, I'm just finalising the re-refinement of a 3.6 Angstrom structure. The stats in general are excellent (Rwork/Rfree = 0.219/0.260; Ramachandran outliers 0.12%; Clashscore < 2) and I've gone through the whole structure to check and double-check fit to the density. The only stat that bugs me a little is the rotamer outliers (2.3%). When I run through a fairly standard refinement protocol (XYZ coords, Individual B-factors, optimizing weights, secondary structure restraints), all seems fine - the Rfree rattles around within 0.2% or so. But when I include the real-space option, very strange things start to happen. Each rsrl step causes a 3-4% jump in both Rwork and Rfree, and jumps the Clashscore up to 40 or so. Rotamer outliers briefly go to zero but then increase substantially in subsequent steps, and the number of Ramachandran outliers increases. I thought I'd missed some dreadful error in my structure until I looked more closely at the results. The attached image is an example of what's happened: an exposed arginine in the middle of a beta-sheet (admittedly with no observable density associated with its sidechain, but register is well fixed by surrounding residues) has been flipped from its original orientation (yellow) to point directly into the hydrophobic core of the domain. Somewhat surprising, to say the least. Is this a bug in the algorithm, or is it simply that 3.6A is too low resolution for it to work reliably? My version of Phenix is a nightly build from about 2 weeks ago. Thanks, Tristan

1 0

Re: [phenixbb] Autobuild iterative omit map error
by Terwilliger, Thomas Charles 20 Feb '15

20 Feb '15

Hi Emily, (Back on-list so others can see the answer.) It seems that your PDB file has a "TER" record between residues ASN 446 and ALA 447 of chain B. Our parsing system interprets that as the end of chain B and thinks that ALA 447 is another chain. You can see this by running: phenix.print_sequence xxxxx.pdb I tried deleting the TER record and it looks as though autobuild will get past the interpretation stage with it taken out. Let me know if that doesn't do it! (Also...do you know why the TER record might be there? It is a little surprising!) All the best, Tom T On Feb 20, 2015, at 9:10 AM, Emilia C. Arturo (Emily) wrote: Thank you Tom! Attached are the reflections file, sequence file, input model pdb file, and the log file from phenix. Emily. On Fri, Feb 20, 2015 at 11:01 AM, Terwilliger, Thomas Charles <terwilliger(a)lanl.gov<mailto:[email protected]>> wrote: Hi Emily, I'm sorry for the trouble! I'll be happy to try and help with it. If you send me (off-list) the input files you are using and the full log file I'll have a look. All the best, Tom T On Feb 20, 2015, at 7:53 AM, Emilia C. Arturo (Emily) wrote: > Hello. > > I'd like to use phenix's Autobuild iterative omit map tool, but have run into an error I could use some help with. When the process aborts, the message that phenix outputs is that the input file has duplicate chains. I looked at the pdb file and I do not see duplicate chains; could the ANISOU records for each atom generated by TLS refinement be hanging up the process? I had no error with Autobuild when I started the iterative map process using an input file that had been generated by phenix.refine using the same configuration as the TLS refinement had used (except, of course, for the TLS parameters). > > Thanks for your attention. > Emily. > > _______________________________________________ > phenixbb mailing list > phenixbb(a)phenix-online.org<mailto:[email protected]> > http://phenix-online.org/mailman/listinfo/phenixbb > Unsubscribe: phenixbb-leave(a)phenix-online.org<mailto:[email protected]> <AutoBuild_run_85_1.log><III-40-C2_justscan3_xscale_cutresto29.mtz><PAH.fa><PAH_C2_ref22plus_refine_84.pdb>

1 0

Autobuild iterative omit map error
by Emilia C. Arturo (Emily) 20 Feb '15

20 Feb '15

Hello. I'd like to use phenix's Autobuild iterative omit map tool, but have run into an error I could use some help with. When the process aborts, the message that phenix outputs is that the input file has duplicate chains. I looked at the pdb file and I do not see duplicate chains; could the ANISOU records for each atom generated by TLS refinement be hanging up the process? I had no error with Autobuild when I started the iterative map process using an input file that had been generated by phenix.refine using the same configuration as the TLS refinement had used (except, of course, for the TLS parameters). Thanks for your attention. Emily.

2 1

on Bulk-solvent and scaling
by Smith Lee 20 Feb '15

20 Feb '15

Dear All, Will you please introduce to me briefly on the function of the step of Bulk-solvent and scaling in the Phenix refine and what is the corresponding mechanism? Smith

2 1

On Phenix refinement
by Smith Liu 20 Feb '15

20 Feb '15

Dear All, The default number of cycles for the Phenix refinement graphical interface is 3 (should be macro-cycles). I have done a refinement with 3 default cycles, and I get the acceptable refinement results. Then I kept everything same, except I changed the cycles of refinement from 3 to 5, however the refinement result was ugly which was unacceptable (Rfree increased). The Phenix document regards 5 macrocycles refine is fine. For my situation, do you think the 3 macrocycles refine result is acceptable or not? If not acceptable, what do you suggest for improvement? I am looking forward to getting your reply. Smith

2 3

erraser for DNA?
by Phoebe A. Rice 20 Feb '15

20 Feb '15

I know phenix has lovely DNA restraints, but I think the backbone of my lowish-resolution DNA structure could use more help. Is there anything like erraser for DNA that will make my DNA backbone as B-like as possible? ++++++++++++++++++++++++++++++++++++++++++ Phoebe A. Rice Dept. of Biochemistry & Molecular Biology The University of Chicago PRice(a)uchicago.edu<mailto:[email protected]>

2 1

The discrepancy for determination of Ramachandran outliers by Coot and MolProbity
by Smith Lee 19 Feb '15

19 Feb '15

Dear All, It often finds for the Ramachandran favored determined by Coot, MolProbity regards as Ramachandran outliers. There are earlier posts regards Coot and MolProbity has different database for the determination of the Ramachandran plots. Then will you please let me know the correct way to correct the Ramachandran outliers by Coot in order to meet the MolProbity Standards? I am looking forward to getting your reply. Smith

1 0

SA (torsion angle) along with B factor reset to low value.
by Tarek DawoD 19 Feb '15

19 Feb '15

to remove memory of the old original R-free flags (which is not available) and using a new test set, i reset the B factor and use SA (torsion angle) i got this error: not implemented: support of different sites. individual and sites.torsion_angles selection. what does that mean? thank you in advance TK_DawoD

2 1

on HIS ** both conformations clash in the flipping
by Smith Lee 17 Feb '15

17 Feb '15

Dear All, In the log output of the Phenix refine, for the flipping of the N/Q/H residues before XYZ refinement, it often occurs something like "HIS ** both conformations clash, **PLEASE CHECK MANUALLY**". Will you please tell me what will clash with 1 of the 2 conformations of the specified HIS? I have checked the specified HIS in Coot, and I still do not know what does it mean for "both conformations clash". I am looking forward to getting a reply from you for clarification. Smith

2 1

Computational Crystallography Newsletter - Volume 6, Number 1
by Nigel Moriarty 17 Feb '15

17 Feb '15

I am pleased to announce the publication of the latest issue of the Computational Crystallography Newsletter: http://www.phenix-online.org/newsletter/ A listing of the articles and short communications is given below. Please note that the newsletter accepts articles of a general nature of interest to all crystallographers. Please send any articles to me at NWMoriarty(a)lbl.gov noting that there is a Word Template on the website to streamline production. January 2015 Articles Improved Probabilistic Estimates of Biocrystal Solvent Content Rapid evaluation of non-bonded overlaps in atomic models Short Communications Plan a SAD experiment, scale SAD data, and analyze your anomalous signal Validation of carbohydrate structures in CCP4 6.5 Disulfide bond restraints Fitting Tip #9: Avoiding excess cis peptides at low resolution or high B FAQ - Tips for coordinated metal refinement Cheers Nigel --- Nigel W. Moriarty Building 64R0246B, Physical Biosciences Division Lawrence Berkeley National Laboratory Berkeley, CA 94720-8235 Phone : 510-486-5709 Email : NWMoriarty(a)LBL.gov Fax : 510-486-5909 Web : CCI.LBL.gov

1 0

Quick request for addition to the line at the end of all these mails
by Morten Grøftehauge 16 Feb '15

16 Feb '15

Dear administrators of this list, Could you please add a line at the very end of these mails with: To unsubscribe click: www.somelink.com As it is it seems that some people have trouble getting off the list after they sign up. Because my spam filter keeps sending phenixbb mails into the spam folder. Which is likely because someone out there keeps flagging them as spam. If there's a link like that then gmail will even offer to attempt to unsubscribe for you. Just a thought. Cheers, Morten

2 1

Please comment on the twinning analysis of my data by Xtriage
by Smith Liu 15 Feb '15

15 Feb '15

Dear All, I have done a Xtriage twinning analysis on my data. I summarize the result as following, Twinning analysis summary: The results of the L-test indicate that the intensity statistics below as expected. No twinning is suspected. Multivariate Z score L-test: 1.625. It further shows, Possible twin operators o merohedral twin operator(s) found, 1 pseudo-merohedral twin operator(s) found. Type Axis R metric(%) delta(le Pa) delta (lebed) Twin law PM 2-fold 4.934 2.547 0.036 l -k h Based on the above information, will you please advise me whether I need to input "l -k h" in the "Twin law" in "Other options" in "Refinement settings" in the Phenix refine graphical interface? I am looking forward to getting your reply. Smith

2 1

On the cycles of the Phenix refinement
by Smith Lee 15 Feb '15

15 Feb '15

Dear All, The default number of cycles for the Phenix refinement graphical interface is 3 (should be macro-cycles). I have done a refinement with 3 default cycles, and I get the acceptable refinement results. Then I kept everything same, except I changed the cycles of refinement from 3 to 5, however the refinement result was ugly which was unacceptable. Someone told me once I have corrected the outliers by Coot, I can have a Phenix refinement with macro-cycles as 1, in order to avoid the Ramachandran corrected ones changing backing to Ramachandran outliers (Will you please further comment whether this strategy to correct Ramachandran outliers is acceptable or not). Thus can we regard the macro-cycles of Phenix refinement as one parameter for the refinement? And why the macro-cycles of the refinement can affect the refinement result? I am looking forward to getting your reply. Smith

1 0

Re: [phenixbb] Fw: Fw: Question related to Phenix refinement
by Pavel Afonine 15 Feb '15

15 Feb '15

> From the attached the image, on the x-axis you can see the 4' rsrg and > 5'rsrg. > > I am looking forward to getting an explanation from you. Dear Smith, hope this will lead you to the answer: chevy:/net/chevy/raid1/afonine/sources % which fifig fifig: aliased to libtbx.find_files -g chevy:/net/chevy/raid1/afonine/sources % fifig rsrg ./cctbx_project/mmtbx/refinement/print_statistics.py: "realsrg" : "rsrg", chevy:/net/chevy/raid1/afonine/sources % All the best, Pavel

1 0

Question related to Phenix refinement
by Smith Lee 15 Feb '15

15 Feb '15

Dear All. In the Phenix refine "run status" window, with the progression of the refinement, it occurs 3' bss, 3' rsrl, 3' xyz, 3' adp (for 3rd round of refinement). I guess xyz means the strategy is for x, y, z coordinates? Then what do bss, rsrl and adp mean? I am looking forward to getting your reply. Smith

9 14

Test
by Paul Adams 14 Feb '15

14 Feb '15

-- Paul Adams Deputy Division Director, Physical Biosciences Division, Lawrence Berkeley Lab Division Deputy for Biosciences, Advanced Light Source, Lawrence Berkeley Lab Adjunct Professor, Department of Bioengineering, U.C. Berkeley Vice President for Technology, the Joint BioEnergy Institute Laboratory Research Manager, ENIGMA Science Focus Area Building 64, Room 248 Building 80, Room 247 Building 978, Room 4126 Tel: 1-510-486-4225, Fax: 1-510-486-5909 http://cci.lbl.gov/paul Lawrence Berkeley Laboratory 1 Cyclotron Road BLDG 64R0121 Berkeley, CA 94720, USA. Executive Assistant: Louise Benvenue [ LBenvenue(a)lbl.gov ][ 1-510-495-2506 ] --

1 0

Hydrogen adding and real space refinement
by Smith Lee 13 Feb '15

13 Feb '15

Dear All, I find when I add H in the PDB file, in a refinement by Phenix, there is no real space refinement step (rsrl). If I do not add H in the PDB file, there is the rsrl step in the refinement. Somebody told me it was because there was a bug in the Phenix I downloaded, somebody told me if H added, there would be no real space refinement. I myself cannot clarify it by reading documents. Will you please clarify whether no rsrl step after H addition was because of bug or because of H addition leading to no rsrl step? And why H addition can lead to no rsrl step? I am looking forward to getting your reply. Best regards. Smith

1 0

How should we estimate the "uncertainty" of the occupancy of an atom?
by Masaki UNNO 12 Feb '15

12 Feb '15

Hello We determined the structure of an enzyme-substrate complex. The substrate of the enzyme contains some partially occupied atoms. We were able to calculate the occupancies of these atoms using PHENIX_refine. By the way, I would like to know the uncertainties of the occupancies of these atoms. How can I estimate the uncertainties? For example, the occupancy of an atom was calculated as 0.65; can we write it as 0.65 (+/- 0.05) or 0.65 (+/- 0.1)? I know that it is related to B-factor and depends on the resolution (number of diffraction data per parameters), but I would like to know how to estimate (calculate?) the uncertainty. Could someone tell me about this issue? Best regards ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Masaki UNNO Graduate School of Science and Engineering Ibaraki University, JAPAN 4-12-1 Nakanarusawa, Hitachi, Ibaraki 316-8511, Japan Tel&Fax +81-294-38-5041 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

6 17

Re: [phenixbb] How should we estimate the "uncertainty" of the occupancy of an atom?
by Dale Tronrud 11 Feb '15

11 Feb '15

On 2/11/2015 12:34 PM, Keller, Jacob wrote: >> If your interpretation of your model depends on whether the occupancy of this atom is 0.65 or 0.75 crystallography is the tool you need. > > Did you mean *not* the tool...? > > JPK > You are right. I meant crystallography is NOT your tool. Sorry. Dale Tronrud

1 0

Specifying Multiple Twin Operators
by Keller, Jacob 11 Feb '15

11 Feb '15

Dear List, Is it possible to specify multiple twin operators in Phenix? I have tetartohedrally-twinned data. If so, what's the syntax, and an example? And can it be done through the GUI? Thanks, Jacob ******************************************* Jacob Pearson Keller, PhD Looger Lab/HHMI Janelia Research Campus 19700 Helix Dr, Ashburn, VA 20147 email: kellerj(a)janelia.hhmi.org *******************************************

2 1

Question on Torsion-angle NCS
by Smith Liu 11 Feb '15

11 Feb '15

Dear All, For the graphical Phenix refine, if we click the "NCS options" button, a small window will occur. In the small window, will you please tell me how much should I input for "Weigh on NCS-related B-factor restraints"? I am looking forward to getting your reply. Smith

2 1

Calculating map with electron scattering factors
by Omar Davulcu 10 Feb '15

10 Feb '15

Hi, everyone. Probably a silly question, but here goes... I'm trying to calculate a simulated map around a set of coordinates using the electron scattering table. Phenix.maps doesn't appear to allow selection of that particular scattering table. I noticed that phenix.refine allows me to select that table and I thought I would set the number of macro cycles to 0 and turn off all the refinement options and produce a map that way. I thought I would generate the structure factors needed as input with phenix.fmodel, but that also doesn't allow selection of the electron scattering table. So, my (naïve) question is: what am I doing wrong? Is there a way I can calculate this map in phenix? Any help would be greatly appreciated. Thanks much! Omar Davulcu

2 3

on XYZ individual (reciprocal space) refinement
by Dialing Pretty 10 Feb '15

10 Feb '15

Dear All, Will you please tell me the XYZ individual (reciprocal space) refinement is executed by which strategy of the Phnix refinement? By XYZ coordinates, or by Optimize X-ray/ADP weight,or by Optimize X-ray/sterochemistry weight? Dialing

2 1