<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=us-ascii"><meta name=Generator content="Microsoft Word 12 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman","serif";
color:black;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p.MsoAcetate, li.MsoAcetate, div.MsoAcetate
{mso-style-priority:99;
mso-style-link:"Balloon Text Char";
margin:0in;
margin-bottom:.0001pt;
font-size:8.0pt;
font-family:"Tahoma","sans-serif";
color:black;}
span.BalloonTextChar
{mso-style-name:"Balloon Text Char";
mso-style-priority:99;
mso-style-link:"Balloon Text";
font-family:"Tahoma","sans-serif";
color:black;}
span.EmailStyle19
{mso-style-type:personal;
font-family:"Calibri","sans-serif";
color:#1F497D;}
span.EmailStyle20
{mso-style-type:personal-reply;
font-family:"Calibri","sans-serif";
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body bgcolor=white lang=EN-US link=blue vlink=purple><div class=WordSection1><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>Pavel,<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'> Your answers have been helpful, thank you. Regarding why I’m not happy with the default settings for the ordered solvent parameters, I am actually expecting them to be appropriate in most cases. However, I think it’s good to “trust but verify” ;) I’m more interested in seeing firsthand how map density and CC correlate for waters in general for a given model/dataset. <o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'> It’s low priority (and probably won’t make too much difference), but if you run across the answer to which 2mFo-DFc map is used for water picking, I’d be interested. I would have guessed it was the original Fobs that was used.<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>Best Regards,<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>-Andy Torelli<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><div><div style='border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in'><p class=MsoNormal><b><span style='font-size:10.0pt;font-family:"Tahoma","sans-serif";color:windowtext'>From:</span></b><span style='font-size:10.0pt;font-family:"Tahoma","sans-serif";color:windowtext'> phenixbb-bounces@phenix-online.org [mailto:phenixbb-bounces@phenix-online.org] <b>On Behalf Of </b>Pavel Afonine<br><b>Sent:</b> Wednesday, January 05, 2011 6:08 PM<br><b>To:</b> phenixbb@phenix-online.org<br><b>Subject:</b> Re: [phenixbb] Filtering ordered solvent molecules based on secondary map<o:p></o:p></span></p></div></div><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Hi Andrew,<br><br><br><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif"'> Thanks for the helpful comments. Do both phenix.model_vs_data and phenix.real_space_correlation produce the same values for CC that phenix.refine uses to keep/remove waters when the ordered solvent routine is used?</span><o:p></o:p></p><p class=MsoNormal><br>yes, they should. They use the same code - that's the joy of having libraries: multiple applications can re-use the same core routines. <br><span style='font-size:11.0pt;font-family:"Calibri","sans-serif"'> </span> <o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif"'> Also, are the default 2mFo-DFc and mFo-DFc map coefficients output by phenix.refine (in the “_maps_coeffs.mtz” file) the same coefficients used by phenix.refine during the ordered solvent checking?</span><o:p></o:p></p><p class=MsoNormal><br>This I would need to check since I don't remember. As I hope all phenix.refine users know, by default phenix.refine writes out two 2mFo-DFc maps (one is using original Fobs set, and one using Fobs set where missing Fobs are filled with DFc), and one mFo-DFc map. Most likely, for water picking phenix.refine uses the one 2mFo-DFc map that uses original set of Fobs.<br><br><br><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif"'> Specifically, does phenix.refine use “filled” maps for ordered solvent checking?</span><o:p></o:p></p><p class=MsoNormal><br>Unlikely.<br><span style='font-size:11.0pt;font-family:"Calibri","sans-serif"'> </span> <o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif"'> If you haven’t guessed, I am thinking of a way to intelligently pick values for the ordered solvent parameters such as primary_map_cutoff, poor_cc_threshold and poor_map_threshold. </span><o:p></o:p></p><p class=MsoNormal><br>Yes, I kind of figured this out... -:) Although I'm wondering why you are not happy with the default settings which are supposed to be good enough most of the time?<br><br><br><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif"'>My idea is to first run a quick round of refinement with generous ordered solvent parameters so that I get a model that is somewhat overpopulated with automatically-picked waters. Then I can manually inspect the mFobs-DFmodel and 2mFobs-DFmodel maps (i.e. primary and secondary maps) and the CC values in order to decide where to select appropriate cutoffs to limit addition of spurious waters for a given model/dataset, which will differ based on map quality, resolution, etc. So obviously I want to calculate CC values and maps in the same way phenix.refine would to judge waters so that I can “see what phenix sees” in this regard.</span><o:p></o:p></p><p class=MsoNormal style='margin-bottom:12.0pt'><br>I see. I think you are on the right track to achieve this. Depending how generous you are setting the water selection criteria you can get as many waters as you want. This is not always bad (as many typically think) but just emulates some earlier versions of ARP idea which is in the end is pretty successful. That is: modeling some density peaks that you cant interpret in terms o your model right now with "dummy waters" may improve the overall map quality which might be helpful.<br><br>Pavel.<br><br><o:p></o:p></p></div></body></html>