<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
</head>
<body dir="auto">
<div>Hi Pavel,</div>
<div><br>
</div>
<div>I think the scaling issue only occurs with pseudo-centering where a significant number of reflections would be systematically weak, which should not be the case with rotational NCS I suppose.&nbsp;</div>
<div><br>
</div>
<div>Best,</div>
<div>Esko<br>
<br>
Sent from my iPhone</div>
<div><br>
On 3 apr 2015, at 20:58, Pavel Afonine &lt;<a href="mailto:pafonine@lbl.gov">pafonine@lbl.gov</a>&gt; wrote:<br>
<br>
</div>
<blockquote type="cite">
<div>Hi Esko,<br>
<br>
thanks for the paper - very interesting!<br>
<br>
Perhaps accounting for ptNCS matters a lot for structure solution or/and refinement/model building at early stages, however I did not see any problem with refinement in this particular example.
<br>
<br>
Here is what I did:<br>
<br>
1) Get model and data from PDB:<br>
<br>
phenix.fetch_pdb 1yup --mtz<br>
<br>
2) Get initial statistics: R-factors from PDB file header, geometry using command:<br>
<br>
phenix.pdbtools model_stat=true 1yup.pdb<br>
<br>
Rwork=0.235 Rfree=0.311<br>
&nbsp;all-atom clashscore : 8.21<br>
&nbsp;ramachandran plot:<br>
&nbsp;&nbsp; outliers : 1.52&nbsp; %<br>
&nbsp;&nbsp; allowed&nbsp; : 7.43&nbsp; %<br>
&nbsp;&nbsp; favored&nbsp; : 91.05 %<br>
&nbsp;rotamer outliers : 14.11 %<br>
&nbsp;cbeta deviations : 3<br>
<br>
3) Then I perturbed the model and refined it with phenix.refine (restrained coordinates and B-factors, water update, weights optimization, local real-space refinement). Here is the statistics for refined model:<br>
<br>
Rwork=0.2380 Rfree=0.2709<br>
&nbsp;all-atom clashscore : 10.61<br>
&nbsp;ramachandran plot:<br>
&nbsp;&nbsp; outliers : 0.80&nbsp; %<br>
&nbsp;&nbsp; allowed&nbsp; : 5.20&nbsp; %<br>
&nbsp;&nbsp; favored&nbsp; : 94.00 %<br>
&nbsp;rotamer outliers : 2.38 %<br>
&nbsp;cbeta deviations : 0<br>
<br>
Just to make it clear: I'm *not* saying we shouldn't account for ptNCS in refinement; in fact, we really should especially given that Randy Read et al worked out all the methods for this!<br>
<br>
I will re-visit this once we implement ptNCS in phenix.refine.<br>
<br>
All the best,<br>
Pavel<br>
<br>
<br>
<div class="moz-cite-prefix">On 4/3/15 9:01 AM, Esko Oksanen wrote:<br>
</div>
<blockquote cite="mid:21782C8D-0497-4913-94B4-81545C8AFCFE@esss.se" type="cite">&nbsp;Hi Mark,
<div><br>
</div>
<div>I believe at least part of the problem arises already in the scaling stage, where the programs assume a unimodal distribution which is not the case in the presence of pseudo-translations. We encountered a similar issue with refinement and pseudo-translation,
 where I simply scaled the systematically weak data separately from the strong and did rigid body refinement against the weak data and restrained against the strong. The details can be found in&nbsp;<a moz-do-not-send="true" href="http://scripts.iucr.org/cgi-bin/paper?ea5052">http://scripts.iucr.org/cgi-bin/paper?ea5052</a>.</div>
<div><br>
</div>
<div>&nbsp; Best,</div>
<div>&nbsp; Esko<br>
<div apple-content-edited="true">
<div style="color: rgb(0, 0, 0); letter-spacing: normal;
            orphans: auto; text-align: start; text-indent: 0px;
            text-transform: none; white-space: normal; widows: auto;
            word-spacing: 0px; -webkit-text-stroke-width: 0px;
            word-wrap: break-word; -webkit-nbsp-mode: space;
            -webkit-line-break: after-white-space;">
<span class="x_Apple-style-span" style="border-collapse:
              separate; orphans: 2; widows: 2;">
<div style="word-wrap: break-word;"><span class="x_Apple-style-span" style="border-collapse:
                  separate;">
<div style="word-wrap: break-word;"><span class="x_Apple-style-span" style="border-collapse:
                      separate;">
<div style="word-wrap: break-word;">
<div><b><font class="x_Apple-style-span" color="#00426C"><font class="x_Apple-style-span" face="Tahoma"><font class="x_Apple-style-span" size="2"><span class="x_Apple-style-span" style="font-size: 10px;"><span class="x_Apple-style-span" style="color: rgb(0, 0, 0);
                                      font-weight: normal;"><font class="x_Apple-style-span" color="#0B466C">_</font><font class="x_Apple-style-span" color="#164A6C">_</font><font class="x_Apple-style-span" color="#204E6C">_</font><font class="x_Apple-style-span" color="#2B536C">_</font><font class="x_Apple-style-span" color="#36576C">_</font><font class="x_Apple-style-span" color="#415B6C">_</font><font class="x_Apple-style-span" color="#4B5F6C">_</font><font class="x_Apple-style-span" color="#56646C">_</font><font class="x_Apple-style-span" color="#61686C">_</font><font class="x_Apple-style-span" color="#6C6C6C">_</font><font class="x_Apple-style-span" color="#858585">_</font><font class="x_Apple-style-span" color="#929292">_</font><font class="x_Apple-style-span" color="#9F9F9F">_</font><font class="x_Apple-style-span" color="#AAAAAA">_</font><font class="x_Apple-style-span" color="#B7B7B7">_</font><font class="x_Apple-style-span" color="#C4C4C4">_</font><font class="x_Apple-style-span" color="#D1D1D1">_</font><font class="x_Apple-style-span" color="#DEDEDE">_</font><font class="x_Apple-style-span" color="#EBEBEB">_</font><font class="x_Apple-style-span" color="#F8F8F8">_</font><font class="x_Apple-style-span" color="#FFFFFF">_</font></span></span></font></font></font></b></div>
<div style=""><b><font class="x_Apple-style-span" color="#00426C"><font class="x_Apple-style-span" face="Tahoma"><font class="x_Apple-style-span" size="2"><span class="x_Apple-style-span" style="font-size: 10px;"><span class="x_Apple-style-span" style="color: rgb(0, 0, 0);
                                      font-weight: normal;"><font class="x_Apple-style-span" color="#FFFFFF"><br>
</font></span></span></font></font></font></b></div>
<div style=""><b><font class="x_Apple-style-span" color="#00426C"><font class="x_Apple-style-span" face="Tahoma"><font class="x_Apple-style-span" size="2"><span class="x_Apple-style-span" style="font-size: 10px;">Esko Oksanen, PhD</span></font></font></font></b></div>
<div><font class="x_Apple-style-span" face="Tahoma"><font class="x_Apple-style-span" color="#00426C"><font class="x_Apple-style-span" size="2"><span class="x_Apple-style-span" style="font-size: 10px;">Instrument Scientist - Neutron Macromolecular Crystallography</span></font></font></font></div>
<div><font class="x_Apple-style-span" face="Tahoma"><font class="x_Apple-style-span" color="#00426C"><font class="x_Apple-style-span" size="2"><span class="x_Apple-style-span" style="font-size: 10px;"><br>
</span></font></font></font></div>
<div><font class="x_Apple-style-span" face="Tahoma"><font class="x_Apple-style-span" color="#00426C"><font class="x_Apple-style-span" size="2"><span class="x_Apple-style-span" style="font-size: 10px;"></span></font></font></font></div>
<div><font class="x_Apple-style-span" face="Tahoma"><font class="x_Apple-style-span" color="#00426C"><font class="x_Apple-style-span" size="2"><span class="x_Apple-style-span" style="font-size: 10px;">European Spallation Source ESS AB</span></font></font></font></div>
<div><font class="x_Apple-style-span" face="Tahoma"><font class="x_Apple-style-span" color="#00426C"><font class="x_Apple-style-span" size="2"><span class="x_Apple-style-span" style="font-size: 10px;">P.O Box 176, SE-221 00 Lund, Sweden</span></font></font></font></div>
<div><font class="x_Apple-style-span" face="Tahoma"><font class="x_Apple-style-span" color="#00426C"><font class="x_Apple-style-span" size="2"><span class="x_Apple-style-span" style="font-size: 10px;">Visiting address: Tunav�gen 24, Lund</span></font></font></font></div>
<div><font class="x_Apple-style-span" color="#00426c" face="Tahoma" size="2"><br>
</font></div>
</div>
</span></div>
</span></div>
</span><span class="x_Apple-style-span" style="border-collapse: separate; orphans: 2; widows: 2;">
<div style="word-wrap: break-word;"><span class="x_Apple-style-span" style="border-collapse:
                  separate;">
<div style="word-wrap: break-word;"><span class="x_Apple-style-span" style="border-collapse:
                      separate;">
<div style="word-wrap: break-word;">
<div><font class="x_Apple-style-span" face="Tahoma"><font class="x_Apple-style-span" color="#00426C"><font class="x_Apple-style-span" size="2"><span class="x_Apple-style-span" style="font-size: 10px;">Switchboard: &#43;46 46 888 00 00</span></font></font></font></div>
<div><font class="x_Apple-style-span" face="Tahoma"><font class="x_Apple-style-span" color="#00426C"><font class="x_Apple-style-span" size="2"><span class="x_Apple-style-span" style="font-size: 10px;">Phone: &#43;46 46 888 30 36</span></font></font></font></div>
<div><font class="x_Apple-style-span" face="Tahoma"><font class="x_Apple-style-span" color="#00426C"><font class="x_Apple-style-span" size="2"><span class="x_Apple-style-span" style="font-size: 10px;">Mobile: &#43;46 721 79 20 36</span></font></font></font></div>
<div><font class="x_Apple-style-span" face="Tahoma"><font class="x_Apple-style-span" color="#00426C"><font class="x_Apple-style-span" size="2"><span class="x_Apple-style-span" style="font-size: 10px;">E-mail:&nbsp;<a moz-do-not-send="true" href="mailto:esko.oksanen@esss.se">esko.oksanen@esss.se</a></span></font></font></font></div>
<div><font class="x_Apple-style-span" face="Tahoma"><font class="x_Apple-style-span" color="#00426C"><font class="x_Apple-style-span" size="2"><span class="x_Apple-style-span" style="font-size: 10px;"><br>
</span></font></font></font></div>
<div><font class="x_Apple-style-span" face="Tahoma"><font class="x_Apple-style-span" color="#00426C"><font class="x_Apple-style-span" size="2"><span class="x_Apple-style-span" style="font-size: 10px;"><a moz-do-not-send="true" href="http://www.esss.se/">www.esss.se</a></span></font></font></font></div>
<div><font class="x_Apple-style-span" face="Tahoma"><font class="x_Apple-style-span" color="#00426C"><font class="x_Apple-style-span" size="2"><br>
</font></font></font></div>
</div>
</span></div>
</span></div>
</span></div>
<br class="Apple-interchange-newline">
<br class="Apple-interchange-newline">
</div>
<br>
<div>
<div>On 03 Apr 2015, at 15:59, Mark A Saper &lt;<a moz-do-not-send="true" href="mailto:saper@umich.edu">saper@umich.edu</a>&gt; wrote:</div>
<br class="Apple-interchange-newline">
<blockquote type="cite">I know that Phaser takes account of NCS translation in the likelihood statistics ( Read, Adams, McCoy (2013) Intensity statistics in the presence of translational non crystallographic symmetry�). Has this also been incorporated into
 phenix.refine? &nbsp;I have a structure with pseudo NCS translations because of NCS rotational symmetry. One suggestion that I found in the archives, was to explicitly LS target refinement rather than ML. &nbsp;Any comments?<br>
<br>
Mark<br>
_________________________________<br>
Mark A. Saper, Ph.D.<span class="Apple-tab-span" style="white-space:pre"> </span>
<span class="Apple-tab-span" style="white-space:pre"></span><span class="Apple-tab-span" style="white-space:pre"></span><span class="Apple-tab-span" style="white-space:pre"></span><br>
Associate Professor of Biological Chemistry, University of Michigan<br>
Ann Arbor MI &nbsp;48109-1055 U.S.A.<br>
<a moz-do-not-send="true" href="mailto:saper@umich.edu">saper@umich.edu</a> &nbsp;&nbsp;&nbsp;&nbsp;phone (734) 764-3353 &nbsp;&nbsp;&nbsp;&nbsp;fax (734) 764-3323<br>
<br>
<br>
<br>
<br>
_______________________________________________<br>
phenixbb mailing list<br>
<a moz-do-not-send="true" href="mailto:phenixbb@phenix-online.org">phenixbb@phenix-online.org</a><br>
<a class="moz-txt-link-freetext" href="http://phenix-online.org/mailman/listinfo/phenixbb">http://phenix-online.org/mailman/listinfo/phenixbb</a><br>
Unsubscribe: <a class="moz-txt-link-abbreviated" href="mailto:phenixbb-leave@phenix-online.org">
phenixbb-leave@phenix-online.org</a><br>
</blockquote>
</div>
<br>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset> <br>
<pre wrap="">_______________________________________________
phenixbb mailing list
<a class="moz-txt-link-abbreviated" href="mailto:phenixbb@phenix-online.org">phenixbb@phenix-online.org</a>
<a class="moz-txt-link-freetext" href="http://phenix-online.org/mailman/listinfo/phenixbb">http://phenix-online.org/mailman/listinfo/phenixbb</a>
Unsubscribe: <a class="moz-txt-link-abbreviated" href="mailto:phenixbb-leave@phenix-online.org">phenixbb-leave@phenix-online.org</a></pre>
</blockquote>
<br>
</div>
</blockquote>
</body>
</html>