<html><head><meta http-equiv="Content-Type" content="text/html charset=iso-8859-1"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Hi,<div><br></div><div>The syntax and behaviour have changed since 2008. Most of the time, people run LLG completion until Phaser finds no more sites, in which case the log-likelihood-gradient map is by definition close to flat at that point. So the default is not to write out the map coefficients. A special flag (LLGMAPS ON) has to be set to get the map coefficients, and then they go into a separate file.</div><div><br></div><div>Here's a sample script to make an initial log-likelihood-gradient map for our MR-SAD tutorial:</div><div><br></div><div><div>phenix.phaser << eof > LLGmap.log</div><div>ROOT LLGmap</div><div>TITLE MR_SAD starting from MR sol with goat alpha-lactalbumin model</div><div>MODE EP_AUTO</div><div>HKLIN lyso2001_scala1.mtz</div><div>LLGCOMPLETE COMPLETE ON</div><div>LLGCOMPLETE NCYC 0</div><div>LLGMAPS ON</div><div>LLGCOMPLETE SCATTERER S</div><div>PARTIAL PDB eptute_1.1.pdb IDENT 0.4</div><div>CRYSTAL four DATASET CuKa LABIN F+ = F_CuKa(+) SIG+ = SIGF_CuKa(+) F- = F_CuKa(-) SIG- = SIGF_CuKa(-)</div><div>COMPOSITION PROTEIN SEQ hewl.pir NUMBER 1</div><div>WAVELENGTH 1.5418</div><div>eof</div><div><br></div><div>You have to say "LLGCOMPLETE COMPLETE ON" to get the LLG map calculated, but "LLGCOMPLETE NCYC 0" to avoid any sites being added, and "LLGMAPS ON" to get the output LLG map, which has coefficients FLLG_S and PHLLG_S for this case (with "S" for the type of scatterer).</div><div><br></div><div>However, we've tested various options as part of the work on ion identification (led by Nat Echols), and if you want to identify anomalous scatterers from a final model it's better to do iterative completion adding purely imaginary scatterers to a model of real scatterers. The refined occupancy gives you a measure of the f" (or f" times occupancy for a partially occupied anomalous scatterer) for the site identified by the imaginary scatterer location, and by doing iterative completion the inclusion of sites found in the initial cycles improves the signal-to-noise in later cycles. Here's a script that could be used for a final model:</div><div><br></div><div><div>phenix.phaser << eof > MRSAD_imag.log</div><div>ROOT MRSAD_imag</div><div>TITLE MR_SAD adding purely imaginary scatterers to final model</div><div>MODE EP_AUTO</div><div>HKLIN lyso2001_scala1.mtz</div><div>LLGCOMPLETE COMPLETE ON</div><div>LLGCOMPLETE SCATTERER AX</div><div>PARTIAL PDB HEWL.pdb IDENT 1.0</div><div>CRYSTAL four DATASET CuKa LABIN F+ = F_CuKa(+) SIG+ = SIGF_CuKa(+) F- = F_CuKa(-) SIG- = SIGF_CuKa(-)</div><div>COMPOSITION PROTEIN SEQ hewl.pir NUMBER 1</div><div>WAVELENGTH 1.5418</div><div>eof</div></div><div><br></div><div>Here what you want to look at is the set of sites (MRSAD_imag.1.pdb) with their occupancies, and not an LLG map (which isn't produced by this script). There are occasionally some spurious sites, especially if the anomalous signal is very weak, but it usually provides good evidence when a site falls on top of what you expect to be a metal ion.</div><div><br></div><div>Best wishes,</div><div><br></div><div>Randy</div><div><br></div><div><div>On 3 Apr 2014, at 17:28, Alastair Fyfe <<a href="mailto:afyfe@SOE.UCSC.EDU">afyfe@SOE.UCSC.EDU</a>> wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite">
<meta content="text/html; charset=ISO-8859-1" http-equiv="Content-Type">
<div bgcolor="#FFFFFF" text="#000000">
<div class="moz-cite-prefix">A related question: is there a current
example of using phaser to calculate SAD LLG maps from a refined
model not refined with phenix.refine? In the past I have used the
description posted here:<br>
<a class="moz-txt-link-freetext" href="http://www.phenix-online.org/pipermail/phenixbb/2008-July/012399.html">http://www.phenix-online.org/pipermail/phenixbb/2008-July/012399.html</a><br>
<br>
However phaser syntax and keywords seem to have changed a bit as
this no longer works. Also, have the derived FLLG and PHLLG
columns been renamed ?<br>
thanks,<br>
Alastair Fyfe<br>
<br>
On 04/03/2014 09:13 AM, Nathaniel Echols wrote:<br>
</div>
<blockquote cite="mid:CALeAa1MSXet8_non0B5-zY+Phe_AE2j8h2Jxm=ApgSyaFXG7Lg@mail.gmail.com" type="cite">
<pre wrap="">PS. Tom pointed out that the anomalous measurability in Xtriage depends on
the sigmas, which is obviously a problem for synthetic data - you can
generate fake sigmas with "add_sigmas=True", but the resulting statistics
will be meaningless.
-Nat
On Thu, Apr 3, 2014 at 8:27 AM, Nathaniel Echols <a class="moz-txt-link-rfc2396E" href="mailto:nechols@lbl.gov"><nechols@lbl.gov></a> wrote:
</pre>
<blockquote type="cite">
<pre wrap="">I guess it depends on what you're looking for as the final output. It's
easy to generate an MTZ file with anomalous Fcalc (this is in the GUI too,
of course):
phenix.fmodel model.pdb high_resolution=2.0 type=real wavelength=0.9792
Extracting some kind of useful summary from the data might require a
little extra scripting - although this may be the kind of thing we should
just add to Xtriage (which only reports "anomalous measurability" right
now).
-Nat
On Thu, Apr 3, 2014 at 7:20 AM, Jonathan Grimes <a class="moz-txt-link-rfc2396E" href="mailto:jonathan@strubi.ox.ac.uk"><jonathan@strubi.ox.ac.uk></a>wrote:
</pre>
<blockquote type="cite">
<pre wrap=""> Given a refined protein structure, is there an straightforward way to
calculate the anomalous
differences as a function of resolution, at wavelength X.
many thanks
jon
Dr. Jonathan M. Grimes,
NDM Senior Reseach Fellow
University Research Lecturer
DIAMOND Research Fellow
Division of Structural Biology
Wellcome Trust Centre for Human Genetics
University of Oxford
Roosevelt Drive,
Oxford OX3 7BN, UK
Email: <a class="moz-txt-link-abbreviated" href="mailto:Jonathan@strubi.ox.ac.uk">Jonathan@strubi.ox.ac.uk</a>, Web: <a class="moz-txt-link-abbreviated" href="http://www.strubi.ox.ac.uk/">www.strubi.ox.ac.uk</a>
Tel: (+44) - 1865 - 287561, FAX: (+44) - 1865 - 287547
_______________________________________________
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>
</pre>
</blockquote>
<pre wrap="">
</pre>
</blockquote>
<pre wrap=""></pre>
<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>
</pre>
</blockquote>
<br>
</div>
_______________________________________________<br>phenixbb mailing list<br><a href="mailto:phenixbb@phenix-online.org">phenixbb@phenix-online.org</a><br>http://phenix-online.org/mailman/listinfo/phenixbb<br></blockquote></div><br><div apple-content-edited="true">
<span class="Apple-style-span" style="border-collapse: separate; font-size: 12px; border-spacing: 0px; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div>------</div><div>Randy J. Read</div><div>Department of Haematology, University of Cambridge</div><div>Cambridge Institute for Medical Research Tel: + 44 1223 336500</div><div>Wellcome Trust/MRC Building Fax: + 44 1223 336827</div><div>Hills Road E-mail: <a href="mailto:rjr27@cam.ac.uk">rjr27@cam.ac.uk</a></div><div>Cambridge CB2 0XY, U.K. www-structmed.cimr.cam.ac.uk</div></span></span>
</div>
<br></div></body></html>