<div dir="ltr"><div>Thanks for that suggestion, Phil. After several weeks, I've finally found solid answers to all my questions, which I'll write here in case anyone reading this message is curious:</div><div><br></div><div>Q1: Is there any way to apply different map weights to different pieces of the structure during refinement?</div>A1: Pavel tells me the answer is no—this feature is not available in Phenix yet. The way I was able to accomplish something to a similar effect was by performing full-model refinements at different weights. Then, in a final refinement job, use torsion restraints from each of those output models (as reference models) across the desired stretches of residues. Clunky, but it solved a few problems for me.<div><br></div><div>Q2: How is the generic "weight" parameter at the bottom of the Phenix GUI applied to the refinement?</div><div>A2: The target function is <span style="color:rgb(0,0,0)">T = Tmap * weight + Trestraints. At the time of writing this message, there is a typo in the Phenix Real-Space Refine documentation, which incorrectly reports the target function as "</span><span style="color:rgb(0,0,0)">T = Tmap + weight * Trestraints" </span><span style="color:rgb(0,0,0)">(Pavel confirmed).</span></div><div><span style="color:rgb(0,0,0)"><br></span></div><div>Q3: Why is the PDB validation telling me I have geometric issues that were not acknowledged by Phenix?</div><div>A3: Phil is right. The ESD for several nucleic acid bond lengths and angles is wide as compared to the sigma tolerated by wwPDB. To fix this, I went into the CIF files within the geostd folder and made the following changes to the troublesome RNA and DNA nucleotide files:</div><div><div><u>data_AD.cif, data_TD.cif, data_CD.cif, data_GD.cif:</u></div><div>length:</div><div>C1'-N1, ESD=0.020</div><div>C3'-O3', ESD=0.020</div><div>angles:</div><div>O4'-C1'-N9/N1, ESD=0.30</div><div><br></div><div><u>data_U.cif</u></div><div>angles:</div><div>C5-C6-N1, ESD=0.5</div><div>C6-N1-C2, ESD=0.5</div><div><br></div><div><u>data_C.cif</u></div><div>angles:</div><div>C6-N1-C2, ESD=0.5</div><div><br></div><div><u>data_A.cif</u></div><div>angles:</div><div>N7-C8-N9, ESD=0.5</div></div><div><br></div><div><br></div><div>Best,</div><div>Josh</div><div><div><br></div><div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Sep 13, 2021 at 7:23 AM Philip D. Jeffrey <<a href="mailto:pjeffrey@princeton.edu">pjeffrey@princeton.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div dir="ltr">
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
:: <span style="font-size:15px;background-color:rgb(255,255,255);display:inline">most of the outliers are not acknowledged
as such (by Phenix), and I only become </span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="font-size:15px;background-color:rgb(255,255,255);display:inline">:: aware of them when the wwPDB validation
tells me they are there.</span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="font-size:15px;background-color:rgb(255,255,255);display:inline"><br>
</span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="font-size:15px;background-color:rgb(255,255,255);display:inline">Check that the equilibrium values that
PDB are using and the ones that Phenix are using are the same. If they're not (or if the sigmas are way different) you're not going to be able to change the outcome just by changing weights.</span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="font-size:15px;background-color:rgb(255,255,255);display:inline"><br>
</span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="font-size:15px;background-color:rgb(255,255,255);display:inline">Phil Jeffrey</span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="font-size:15px;background-color:rgb(255,255,255);display:inline">Princeton</span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="font-size:15px;background-color:rgb(255,255,255);display:inline"><br>
</span></div>
<div id="gmail-m_4782705509461365504appendonsend"></div>
<hr style="display:inline-block;width:98%">
<div id="gmail-m_4782705509461365504divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> <a href="mailto:phenixbb-bounces@phenix-online.org" target="_blank">phenixbb-bounces@phenix-online.org</a> <<a href="mailto:phenixbb-bounces@phenix-online.org" target="_blank">phenixbb-bounces@phenix-online.org</a>> on behalf of Josh Cofsky <<a href="mailto:josh.cofsky@gmail.com" target="_blank">josh.cofsky@gmail.com</a>><br>
<b>Sent:</b> Sunday, September 12, 2021 9:42 PM<br>
<b>To:</b> <a href="mailto:phenixbb@phenix-online.org" target="_blank">phenixbb@phenix-online.org</a> <<a href="mailto:phenixbb@phenix-online.org" target="_blank">phenixbb@phenix-online.org</a>><br>
<b>Subject:</b> Re: [phenixbb] Restricting real-space refinement and weights to specific atom selections</font>
<div> </div>
</div>
<div>
<div dir="ltr">
<div>I now realize that I am confused as to what the mentioned "weight" parameter refers to. The Phenix RSR documentation contains a line that reads "<span style="color:rgb(0,0,0)">T = Tmap + weight * Trestraints"—if I'm understanding correctly (?), this means
a larger "weight" value should yield greater conformity to geometric restraints at the expense of model/map fit. In reality, I observe the opposite trend in my refinements (which is why I assumed it was a "map weight"). Can someone clarify what</span><span style="color:rgb(0,0,0)"> "weight"
(Refinement Settings > "Other Options") refers to?</span></div>
<div><br>
</div>
<div><span style="color:rgb(0,0,0)">The reason I am exploring these questions is that Phenix RSR is introducing bond angle outliers into a low-resolution area of a DNA molecule in my model—most of the outliers are not acknowledged as such (by Phenix), and I
only become aware of them when the wwPDB validation tells me they are there. I assumed this was the result of overfitting to the map because setting "weight" to 0.0 prevented these outliers from appearing, but at the expense of map-model fit in the remainder
of the structure (which is why I wanted to apply different weights to different segments of the model).</span><span style="color:rgb(0,0,0)"><br>
</span></div>
<div><span style="color:rgb(0,0,0)"><br>
</span></div>
<div><span style="color:rgb(0,0,0)">Thanks,</span></div>
<div>Josh<br>
</div>
</div>
<br>
<div>
<div dir="ltr">On Sun, Sep 12, 2021 at 12:21 PM Josh Cofsky <<a href="mailto:josh.cofsky@gmail.com" target="_blank">josh.cofsky@gmail.com</a>> wrote:<br>
</div>
<blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div dir="ltr">Hi all,
<div><br>
</div>
<div>I'm using phenix.real_space_refine to refine a model into a cryo-EM map. My local map resolution is highly variable across the molecule, and using a single map weight (the parameter shown in the Refinement Settings > "Other Options" in the GUI) on the
whole model leads to poor map-to-model fit in some locations and overfitting in others. Potential solutions could be to:</div>
<div>1. assign different map weights to different segments of the model (based on a user-specified atom selection)</div>
<div>OR</div>
<div>2. perform iterative refinements with different weight settings, holding different (user-specified) segments of the model fixed during each iteration (I suspect this isn't possible based on Pavel's message here: <a href="http://phenix-online.org/pipermail/phenixbb/2021-May/025037.html" target="_blank">http://phenix-online.org/pipermail/phenixbb/2021-May/025037.html</a>)</div>
<div><br>
</div>
<div>Is either of these approaches (or another one that accomplishes the same goal) possible? I am using the Phenix GUI, but if this is only possible in command line, and you're willing to share an example command that accomplishes either of these tasks, please
let me know.</div>
<div><br>
</div>
<div>Thanks!</div>
<font color="#888888">
<div>-Josh</div>
<div><br>
</div>
<div>(also, apologies if this is a duplicate message...still trying to figure out the phenixbb listserv)</div>
</font></div>
</blockquote>
</div>
</div>
</div>
</blockquote></div></div></div></div>