<HTML>
<HEAD>
<TITLE>Re: [phenixbb] fixing geometry for deposition</TITLE>
</HEAD>
<BODY>
<FONT FACE="Geneva, Verdana, Helvetica, Arial"><SPAN STYLE='font-size:11pt'>Hi Raja, <BR>
I would like to add to Pavel&#8217;s excellent advice and ask if you have tried TLS refinement? It always improves the model in our hands, and in most structures ( <a href="http://journals.iucr.org/j/issues/2009/03/00/he5437/index.html">http://journals.iucr.org/j/issues/2009/03/00/he5437/index.html</a> ). You might also consider excluding the regions with poor geometry from the NCS groups, to see if they are in fact different. You are also forcing NCS regardless of the RMS distance. I would refine a few rounds without NCS and with the default refinement parameters, and then then turn the NCS back on. Look in the .geo file for excessive distances, and examine these regions in Coot using the NCS ghosts and maps to decide if they should be excluded from the NCS groups. <BR>
<BR>
Another suggestion would be to try the kicked maps, or generate &nbsp;composite omit maps to look at the regions with poor geometry. <BR>
<BR>
You might need to subtly change half a dozen residues surrounding the geometry outlier to fix the problem. This might be done with simulated annealing, or try &nbsp;autobuild in place. <BR>
<BR>
Lastly, if its only a few residues, you can&#8217;t fix it, and the maps support your model, don&#8217;t worry about it. <BR>
<BR>
Best Regards,<BR>
Kendall Nettles<BR>
</SPAN></FONT><SPAN STYLE='font-size:11pt'><FONT FACE="Calibri, Verdana, Helvetica, Arial"><BR>
<BR>
On 11/17/09 3:02 PM, &quot;Raja Dey&quot; &lt;<a href="rdey@usc.edu">rdey@usc.edu</a>&gt; wrote:<BR>
<BR>
</FONT></SPAN><BLOCKQUOTE><SPAN STYLE='font-size:11pt'><FONT FACE="Calibri, Verdana, Helvetica, Arial">Hi Pavel,<BR>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;I used the following commands at the last stage of refinement:<BR>
<BR>
phenix.refine 1205A-p21212.mtz nov138_001.pdb simulated_annealing=false ncs_groups.params main.ncs=true ncs.find_automatically=false refinement.ncs.excessive_distance_limit=None main.bulk_solvent_and_scale=false strategy=rigid_body+individual_sites xray_data.high_resolution=2.9 refinement.main.number_of_macro_cycles=1 output.write_maps=true output.prefix=nov139 --overwrite<BR>
<BR>
Regarding amino acids phy, psy..... I tried real space refinement/regularization in coot, but getting the same list of amino acids again and again in the pdb validation summary letter. I also tried mutation on its own to get the ideal geometry, still got the same result.<BR>
<BR>
Did I do anything wrong?<BR>
<BR>
<BR>
Raja<BR>
<BR>
----- Original Message -----<BR>
From: Pavel Afonine &lt;<a href="PAfonine@lbl.gov">PAfonine@lbl.gov</a>&gt;<BR>
Date: Tuesday, November 17, 2009 11:50 am<BR>
Subject: Re: [phenixbb] fixing geometry for deposition<BR>
To: PHENIX user mailing list &lt;<a href="phenixbb@phenix-online.org">phenixbb@phenix-online.org</a>&gt;<BR>
<BR>
&gt; Hi Raja,<BR>
&gt;<BR>
&gt; &gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;pdb validation summary letter gives a list<BR>
&gt; of covalent bond angles greater than 6 times standard deviation for<BR>
&gt; dna part of my molecule. I am wondering if someone can tell how to<BR>
&gt; fix that. I used PHENIX to refine my molecule. I did not use any<BR>
&gt; restraints for the dna part.<BR>
&gt; &gt; &nbsp;<BR>
&gt;<BR>
&gt; did you completely turn off all the geometry restraints for that<BR>
&gt; selected part of your structure? If so, then it is not surprising<BR>
&gt; that<BR>
&gt; it got distorted during refinement at 2.9A resolution, and<BR>
&gt; therefore I<BR>
&gt; would suggest re-running phenix.refine using all the restraints<BR>
&gt; that<BR>
&gt; phenix.refine normally uses by default.<BR>
&gt;<BR>
&gt; &gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Also a few amino acids are listed for which<BR>
&gt; phy and psy fall outside Ramachandran plot. I like to know what is<BR>
&gt; the best way to fix these phy, psy angles.<BR>
&gt; &gt; &nbsp;<BR>
&gt;<BR>
&gt; 1) Can't you do it in Coot? (especially when there are only a few<BR>
&gt; such<BR>
&gt; outliers, as you say).<BR>
&gt;<BR>
&gt; 2) Try adding riding H atoms and re-run phenix.refine (with and w/o<BR>
&gt; weights optimization to see which option gives better results).<BR>
&gt;<BR>
&gt; 3) Try a quick geometry regularization going into trying &quot;2)&quot;.<BR>
&gt; Repeat<BR>
&gt; the same with and w/o H.<BR>
&gt;<BR>
&gt; This is a pretty frequent question, so I hope someone someone can<BR>
&gt; tell a<BR>
&gt; success story about it...<BR>
&gt;<BR>
&gt; &gt; Final refined parameters are<BR>
&gt; &gt;<BR>
&gt; &gt; REMARK Start: r_work = 0.3576 r_free = 0.3345 bonds = 0.009<BR>
&gt; angles = 1.344 &nbsp;&nbsp;&nbsp;&nbsp;<BR>
&gt; &gt; REMARK Final: r_work = 0.2517 r_free = 0.2845 bonds = 0.008<BR>
&gt; angles = 1.209 &nbsp;&nbsp;&nbsp;&nbsp;<BR>
&gt; &gt;<BR>
&gt; &gt; at 2.9A resolution<BR>
&gt;<BR>
&gt; Looks reasonable given the resolution...<BR>
&gt;<BR>
&gt; Pavel.<BR>
&gt;<BR>
&gt; _______________________________________________<BR>
&gt; phenixbb mailing list<BR>
&gt; <a href="phenixbb@phenix-online.org">phenixbb@phenix-online.org</a><BR>
&gt; <a href="http://www.phenix-online.org/mailman/listinfo/phenixbb">http://www.phenix-online.org/mailman/listinfo/phenixbb</a><BR>
&gt;<BR>
_______________________________________________<BR>
phenixbb mailing list<BR>
<a href="phenixbb@phenix-online.org">phenixbb@phenix-online.org</a><BR>
<a href="http://www.phenix-online.org/mailman/listinfo/phenixbb">http://www.phenix-online.org/mailman/listinfo/phenixbb</a><BR>
<BR>
</FONT></SPAN></BLOCKQUOTE>
</BODY>
</HTML>