- phenixbb - phenix-online.org

message after refinement
by Almudena Ponce Salvatierra 14 Jul '15

14 Jul '15

Dear all, after refinement of my nucleic acid structure with phenix.refine I get: "warning: rms deviation from ideal geometry (bonds and/or angles) is suspisciously high. You may need to optimize X-ray stereochemistry weighting" However when I move to the molprobity tab I see no outliers no anything... looking weird. I even more wonder where this comes from because the fact is that I had one refinement with no warnings or anything and today I repeated the same, exactly the same, except for the fact that I renumbered one of the chains. Now all out of a sudden (I know it doesn't sound scientific at all), I get this warning regarding the rmsd when I was actually considering that the refinement stage was over. Does anybody have a theory regarding where the problem might be here? Thanks a lot in advance. Best wishes, Almudena -- Almudena Ponce-Salvatierra Macromolecular crystallography and Nucleic acid chemistry Max Planck Institute for Biophysical Chemistry Am Fassberg 11 37077 Göttingen Germany

2 1

Refinement with new mtz.file
by Aleksandar Bijelic 14 Jul '15

14 Jul '15

Dear CCP4 user and experts, I refined (with PHENIX) a 3.0 A dataset obtaining Rfree of about 0.24 (with good geometry according to Ramachandran, Beta Outliers, etc.) .... Everything seems to be ok (especially in relation to the resolution) .... because the mtz.file I used was quite old and I cannot find my xscale.hkl file, I processed the data set again (this time with optimizng and polishing) and received a "better" file according to almost everything (resolution limit, I/sigma, CC(1/2), Rmeas), thus I used this new mtz.file and put it in my last refinement step (the refinement which led to the above mentioned Rfree = 0.24). Suprisingly, the refinement starts at Rwork = 0.18 and Rfree = 0.19 but ending up with 0.20 and 0.22, respectively. So I wanted to know if this is usual? I was expecting my data to become slightlly better but what is irritating me is the starting R-values of the refinement and that it get worse during refinement. Maybe I did something wrong? Is it reasonable to replace the mtz.file with a new one in the last refinement step or should I start the refinement from the scratch? In the sake of completeness, I deleted the header of the used pdb.file because of the R-flag error which occurs, since Phenix reconizes that the pdb.file was already used with other Rflags. Sorry, but I am still a beginner in this field, so I would be very grateful if somebody could explain me this situation and my mistake and if I need to start refinement from the beginning. Thank you in advance! Best Regards. Aleksandar -- ------------------------------------------- Aleksandar Bijelic, MSc. Institut für Biophysikalische Chemie Universität Wien Althanstrasse 14 A-1090 Wien Tel: +43 1 4277 52536 e-Mail: aleksandar.bijelic(a)univie.ac.at --------------------------------------------

7 10

Set up Simulated Annealing starting tempreture
by Xiao Lei 14 Jul '15

14 Jul '15

Dear Phenixbb members, I'm learning phenix this days. I use phenix refinement and found there is a simulated annealing option (SA Cartesian or SA torsion option). However, I found the SA cartesian starts at 5000K and SA torsion starts at 2500K by default, and I'd like to change the starting temperature by myself. I search phenix manual and found: To change starting SA temperature: % phenix.refine data.hkl model.pdb simulated_annealing=true \ strategy=individual_sites simulated_annealing.start_temperature=10000 I do not understand this, does this mean I have to type the above command in somewhere in phenix? is there a place I can simply put SA start temperature numbers in? Xiao

2 2

NCS and R-Values Hung up at 0.42/0.46 for 2.1Å resolution
by John Bruning 11 Jul '15

11 Jul '15

I collected two datasets of the same ligand bound structure at 2.4 and 2.1 Å resolution. The apo structure has already been solved to 2.2Å resolution in P31 (and P3121 for SeMet). After processing each dataset as both P31 and P3121 in MOSFLM, I proceeded with MR with Phaser using sequential partial search models to place either 3 copies in the ASU for P3121, or 6 copies for P31. (The previously solved apo structure had 6 copies of the same protein in the ASU). *** in summary two datasets, processed as both P31 or P3121 After running xtriage for either .mtz, I get this warning about NCS: Distance to origin : 56.728 Height relative to origin : 26.388 % Translational pseudo-symmetry is very likely present in these data. Be aware that this will change the intensity statistics and may impact subsequent analyses, and in practice may lead to higher R-factors in refinement. Regardless, I can proceed with the MR using sequential partial search models to place either 3 copies (P3121) or 6 copies (P31) with decent solutions (i.e. packing looks good upon inspection, scoring looks appropriate -> SOLU SET RFZ=2.8 TFZ=8.3 PAK=13*LLG*=-37 *LLG*=2199 as an example). However, I did get the following somewhat expected warnings: Number of known components (1) not divisible by 2: Translational *NCS* correction not applied Large non-origin Patterson peak indicates that translational *NCS* is present Upon any attempt at initial refinement (with a number of different strategies), the R-factors are hung up > 0.42 (for example in the P3121, 2.1Å case): Start R-work = 0.4966, R-free = 0.5133; Final R-work = 0.4231, R-free = 0.4477 My concern is that there is some pathology that is not being recognized (merohedral twinning) or that the NCS/pseudosymmetry is not being addressed appropriately. -- Cheers, John

1 0

Re: [phenixbb] Fwd: [ccp4bb] Phenix; converting a map to structure factors
by Pavel Afonine 11 Jul '15

11 Jul '15

Hi Mark, Tom copied me your email. > As I'm sure you are aware there is tremendous interest in building > into and refining against low resolution (by X-ray standards) EM maps. > Could you comment to the phenix BB about whether its appropriate to > convert maps to reciprocal space and use ML refinement? While it is a possible way of doing this (and this is how it's been done in the past) I don't think it is a good route to follow. In case of cryo-EM your data is map, so why not just refine against the map? Why take a detour through reciprocal-space by converting map into structure factors (=manipulating your original data) and doing reciprocal space using refinement programs that were designed for crystallography? With this in mind we have developed a new refinement tool in Phenix (it is still in active development!) called phenix.real_space_refine that is designed to refine atomic models directly into maps (any map: X-ray, neutron, cryo-EM, low or high resolution), for quick summary see http://phenix-online.org/presentations/real_space_refine.pdf as well as documentation page (at the bottom): http://phenix-online.org/documentation/ and newsletter article (page 43): http://phenix-online.org/newsletter/CCN_2013_07.pdf Let me know if you have any questions or comments! Pavel

2 4

atom vdW dictionary for Python
by David Shin 11 Jul '15

11 Jul '15

Hi - I'm writing a Python script to measure the volume of pockets within protein structures - plus extra bells and whistles... thus a custom job. I was wondering just for speed, if someone could point me to values for van der Waals radii that are somewhat agreed upon for protein atoms, or if there is some dictionary already made that someone could point to that resides in phenix or some other structural biology package that would be ok for me to copy/use. Thanks, Dave

4 3

Re: [phenixbb] [ccp4bb] Phenix; converting a map to structure factors
by Terwilliger, Thomas Charles 11 Jul '15

11 Jul '15

Hi Reza, Just use: phenix.map_to_structure_factors map.ccp4 d_min=2.5 or something like that, All the best, Tom T On Jul 10, 2015, at 11:27 AM, Reza Khayat wrote: Hi, Is it possible to convert a map to structure factors using Phenix? Thanks. Best wishes, Reza Reza Khayat, PhD Assistant Professor City College of New York Department of Chemistry New York, NY 10031

3 4

NCS and R-Values Hung up at ~0.42/0.46
by John Bruning 10 Jul '15

10 Jul '15

Hi Everyone, I collected two datasets of the same ligand bound structure at 2.4 and 2.1 Å resolution. The apo structure has already been solved to 2.2Å resolution in P31 (and P3121 for SeMet). After processing each dataset as both P31 and P3121 in MOSFLM, I proceeded with MR with Phaser using sequential partial search models to place either 3 copies in the ASU for P3121, or 6 copies for P31. (The previously solved apo structure had 6 copies of the same protein in the ASU). *** in summary two datasets, processed as both P31 or P3121 After running xtraige for either .mtz, I get this warning about NCS: Largest Patterson peak with length larger than 15 Angstrom: Frac. coord. : 0.333 -0.333 -0.153 Distance to origin : 56.728 Height relative to origin : 26.388 % p_value(height) : 2.659e-03symmetry. Translational pseudo-symmetry is very likely present in these data. Be aware that this will change the intensity statistics and may impact subsequent analyses, and in practice may lead to higher R-factors in refinement. Regardless, I can proceed with the MR using sequential partial search models to place either 3 copies (P3121) or 6 copies (P31) with decent solutions (i.e. packing looks good upon inspection, scoring looks appropriate -> SOLU SET RFZ=2.8 TFZ=8.3 PAK=13*LLG*=-37 *LLG*=2199 as an example). However, I did get the following somewhat expected warnings: Number of known components (1) not divisible by 2: Translational *NCS* correction not applied Large non-origin Patterson peak indicates that translational *NCS* is present Upon any attempt at initial refinement (with a number of different strategies), the R-factors are hung up > 0.42 (for example in the P3121, 2.1Å case): Start R-work = 0.4966, R-free = 0.5133 Final R-work = 0.4231, R-free = 0.4477 *** best R-factors from ~ 15 different strategies My concern is that there is some pathology that is not being recognized (merohedral twinning) or that the NCS/pseudosymmetry is not being addressed appropriately. Any suggestions in troubleshooting this would be much appreciated, as I have pretty much exhausted everything I could think of. -- Cheers, John

2 1

(no subject)
by Ashok Nayak 10 Jul '15

10 Jul '15

Dear Phenixites, I have a helical bundle structure, where in I wish to understand by what rotation and about which axis the dimers are inter-related. What I get from the superposition with Superpose is the transformation matrix and not the axis, and the angles are also elusive. Ashok Nayak CSIR-CDRI Lucknow, India --

2 1

Axis of Rotation
by Ashok Nayak 10 Jul '15

10 Jul '15

Dear Phenixites, I have a helical bundle structure, where in I wish to understand by what rotation and about which axis the dimers are inter-related. What I get from the superposition with Superpose is the transformation matrix and not the axis, and the angles are also elusive. Thanks in advance -- Ashok Senior Research Fellow - Dr JV Pratap, Lab No-LSN 008 Molecular and Structural Biology Division Central Drug Research Institute, Janakipuram Extension Lucknow-226031 India

2 1