<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" id="owaParaStyle"></style>
</head>
<body fpstyle="1" ocsi="0">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">
<div>Hi Mohamed,</div>
<div><br>
</div>
<div>This looks generally fine to me. &nbsp;One thing that stands out though is the line:</div>
<div><br>
</div>
<div>4.0- 3.5 &nbsp;0.79 &nbsp; 3.23 &nbsp;0.24 &nbsp; &nbsp;4001 &nbsp;2577</div>
<div><br>
</div>
<div>where the half-dataset CC is 0.24, much greater than for lower resolutions. &nbsp;This suggests that something is systematically wrong in this resolution shell (or elsewhere). &nbsp;It may be useful to cut the data at 4 A just for this reason. &nbsp;As Diana points out
 it may be useful for other reasons as well to cut the resolution at 4 A for finding the sites.</div>
<div><br>
</div>
<div>I would suggest taking all of the datasets you have that are more or less isomorphous and putting them in a directory and then using scale_and_merge to put them all together (this is what it is designed for). &nbsp;It should be able to down-weight the ones
 that are most different from the average.</div>
<div><br>
</div>
<div>You could also try the brute-force option in phenix.hyss to try very hard to find the sites.</div>
<div><br>
</div>
<div>All the best,</div>
<div>Tom T</div>
<div><font face="Times" size="3"><br>
</font></div>
<div>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div id="divRpF906049" style="direction: ltr;"><font face="Tahoma" size="2" color="#000000"><b>From:</b> mohamed noor [mohamed.noor34@gmail.com]<br>
<b>Sent:</b> Friday, August 07, 2015 10:55 AM<br>
<b>To:</b> Terwilliger, Thomas Charles<br>
<b>Cc:</b> PHENIX user mailing list<br>
<b>Subject:</b> Re: [phenixbb] Strong anomalous signal but AutoSol fails<br>
</font><br>
</div>
<div></div>
<div>
<div dir="ltr">
<div>
<div>Hi Tom<br>
<br>
</div>
Just to check that I am on the right track, I have attached below the output from phenix.anomalous_signal using a single XDS_ASCII.HKL file from a single crystal. It is one of a few datasets, but based on Xtriage seems to be the best one.<br>
<br>
</div>
The space group is P 3(1/2) 2 1. All datasets were collected at the peak wavelength (fluorescence scan). My initial feeling that there is a strong signal comes from CORRECT.LP. Aimless and Xtriage. IIRC, the resolution difference between the optimistic and
 pessimistic measurability is about 0.3-0.4 A.<br>
<div><br>
<br>
Estimation of anomalous signal in a dataset<br>
<br>
Estimating B-value for anomalous substructure as&nbsp;&nbsp;&nbsp; 91.2&nbsp; based on <br>
overall B-value of&nbsp;&nbsp;&nbsp; 75.4 (Note: you can set this with b_value_anomalous=xx)<br>
<br>
Getting scaled data and half-datasets with scale_and_merge<br>
Log file will be: scale.log<br>
<br>
Files for half_dataset CC: ['TEMP4/XDS_ASCII_SG150_1800frames_reprocFriedel_0_1.sca']
<br>
Files for half_dataset CC: ['TEMP4/XDS_ASCII_SG150_1800frames_reprocFriedel_0_1.sca']
<br>
Files for half_dataset CC: ['TEMP4/XDS_ASCII_SG150_1800frames_reprocFriedel_0_1.sca']
<br>
Scaled data are in: scaled_data.mtz<br>
Half-dataset A is in: half_dataset_a.mtz<br>
Half-dataset B is in: half_dataset_b.mtz<br>
Using scaled data in analysis<br>
<br>
<br>
Setting up estimator for CC*<br>
<br>
-------------------Summary of signal in this dataset ------------------------<br>
<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Shell<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; CCano&nbsp;&nbsp; Nrefl Nrefl<br>
Resolution Esqr I/sigI&nbsp; half&nbsp;&nbsp; anom&nbsp;&nbsp; half<br>
47.8- 7.0&nbsp; 0.50&nbsp; 19.55&nbsp; 0.32&nbsp;&nbsp;&nbsp; 2389&nbsp; 2360<br>
&nbsp;7.0- 6.5&nbsp; 0.69&nbsp; 10.63&nbsp; 0.11&nbsp;&nbsp;&nbsp;&nbsp; 621&nbsp;&nbsp; 608<br>
&nbsp;6.5- 6.0&nbsp; 0.87&nbsp;&nbsp; 8.47&nbsp; 0.08&nbsp;&nbsp;&nbsp;&nbsp; 850&nbsp;&nbsp; 818<br>
&nbsp;6.0- 5.5&nbsp; 0.88&nbsp;&nbsp; 7.52&nbsp; 0.08&nbsp;&nbsp;&nbsp; 1189&nbsp; 1141<br>
&nbsp;5.5- 5.0&nbsp; 0.72&nbsp;&nbsp; 5.88&nbsp; 0.09&nbsp;&nbsp;&nbsp; 1709&nbsp; 1588<br>
&nbsp;5.0- 4.5&nbsp; 0.76&nbsp;&nbsp; 4.59&nbsp; 0.07&nbsp;&nbsp;&nbsp; 2521&nbsp; 2247<br>
&nbsp;4.5- 4.0&nbsp; 0.98&nbsp;&nbsp; 3.25&nbsp; 0.05&nbsp;&nbsp;&nbsp; 3881&nbsp; 3268<br>
&nbsp;4.0- 3.5&nbsp; 0.79&nbsp;&nbsp; 3.23&nbsp; 0.24&nbsp;&nbsp;&nbsp; 4001&nbsp; 2577<br>
&nbsp;3.5- 3.3&nbsp; 1.72&nbsp;&nbsp; 1.99 -0.02&nbsp;&nbsp;&nbsp; 3420&nbsp; 2118<br>
<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Cumulative<br>
<br>
----------------------Data quality-----------------&nbsp;&nbsp;&nbsp; Best guess of expected<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; results of finding sites<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ------ and phasing--------<br>
<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; CCano&nbsp;&nbsp; Nrefl&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; P(Substr)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>
Resolution Skew Esqr&nbsp; half&nbsp;&nbsp; anom&nbsp;&nbsp;&nbsp; CC* Signal&nbsp; &#43;/-&nbsp;&nbsp;&nbsp;&nbsp; (%)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; FOM*&nbsp; &#43;/-<br>
47.8- 7.0&nbsp; 0.02 0.47&nbsp; 0.32&nbsp;&nbsp;&nbsp; 2389&nbsp; 0.51&nbsp; 10.1&nbsp;&nbsp; 1.3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 68&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.2&nbsp;&nbsp; 0.0<br>
47.8- 6.5&nbsp; 0.02 0.50&nbsp; 0.28&nbsp;&nbsp;&nbsp; 3010&nbsp; 0.48&nbsp; 10.7&nbsp;&nbsp; 1.7&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 72&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.2&nbsp;&nbsp; 0.1<br>
47.8- 6.0&nbsp; 0.02 0.57&nbsp; 0.24&nbsp;&nbsp;&nbsp; 3860&nbsp; 0.43&nbsp; 10.8&nbsp;&nbsp; 2.8&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 72&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.2&nbsp;&nbsp; 0.1<br>
47.8- 5.5&nbsp; 0.05 0.64&nbsp; 0.20&nbsp;&nbsp;&nbsp; 5049&nbsp; 0.44&nbsp; 12.6&nbsp;&nbsp; 2.4&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 78&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.2&nbsp;&nbsp; 0.1<br>
47.8- 5.0&nbsp; 0.01 0.66&nbsp; 0.17&nbsp;&nbsp;&nbsp; 6758&nbsp; 0.39&nbsp; 12.8&nbsp;&nbsp; 3.3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 78&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.2&nbsp;&nbsp; 0.1<br>
47.8- 4.5&nbsp; 0.01 0.69&nbsp; 0.14&nbsp;&nbsp;&nbsp; 9279&nbsp; 0.36&nbsp; 13.4&nbsp;&nbsp; 3.4&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 82&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.2&nbsp;&nbsp; 0.1<br>
47.8- 4.0&nbsp; 0.00 0.77&nbsp; 0.12&nbsp;&nbsp; 13160&nbsp; 0.33&nbsp; 14.3&nbsp;&nbsp; 3.6&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 88&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.2&nbsp;&nbsp; 0.1<br>
47.8- 3.5&nbsp; 0.00 0.77&nbsp; 0.15&nbsp;&nbsp; 17161&nbsp; 0.39&nbsp; 18.4&nbsp;&nbsp; 3.8&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 98&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.3&nbsp;&nbsp; 0.1<br>
47.8- 3.3&nbsp; 0.00 0.88&nbsp; 0.14&nbsp;&nbsp; 20581&nbsp; 0.31&nbsp; 15.7&nbsp;&nbsp; 2.1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 95&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.2&nbsp;&nbsp; 0.0<br>
<br>
<div class="gmail_extra"><br>
<div class="gmail_quote">On Fri, Aug 7, 2015 at 2:30 PM, Terwilliger, Thomas Charles
<span dir="ltr">&lt;<a href="mailto:terwilliger@lanl.gov" target="_blank">terwilliger@lanl.gov</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex; border-left:1px #ccc solid; padding-left:1ex">
Hi Mohamed,<br>
<br>
You might try running phenix.anomalous_signal on your data (may require finding your unmerged data and running phenix.scale_and_merge first).&nbsp; This will give you an idea if you should be able to solve your SAD dataset.<br>
<br>
See:&nbsp; <a href="http://www.phenix-online.org/version_docs/1.10pre-2124/reference/anomalous_signal.html" rel="noreferrer" target="_blank">
http://www.phenix-online.org/version_docs/1.10pre-2124/reference/anomalous_signal.html</a><br>
<br>
All the best,<br>
Tom T<br>
<br>
<br>
<br>
<br>
<br>
From: <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 mohamed noor [<a href="mailto:mohamed.noor34@gmail.com" target="_blank">mohamed.noor34@gmail.com</a>]<br>
<br>
Sent: Thursday, August 06, 2015 3:16 PM<br>
<br>
To: PHENIX user mailing list<br>
<br>
Subject: [phenixbb] Strong anomalous signal but AutoSol fails<br>
<div class="HOEnZb">
<div class="h5"><br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
Dear developers<br>
<br>
<br>
<br>
<br>
I have a low resolution anomalous dataset which Aimless suggests has an effective resolution to 3.3 A and anomalous signal to 3.5 A. However, SAD phasing with AutoSol is not successful with the final R factor around 50 %.<br>
<br>
<br>
<br>
<br>
<br>
I also have another dataset collected at a remote wavelength without anomalous signal to 3 A but they are not isomorphous (&gt; 2 A difference in c axis).<br>
<br>
<br>
<br>
<br>
<br>
<br>
The anomalous signal comes from the ligand heme c, which is bound covalently to the protein, so its occupancy should be 1. The protein is quite small with about 120 residues. Xtriage suggests an NCS of 6 to 20 with most likely number to be 13.<br>
<br>
<br>
<br>
<br>
Is there any reason why a reasonable solution cannot be found? There is no twinning.<br>
<br>
<br>
<br>
<br>
I am using the latest nightly 1.10 pre2124.<br>
<br>
<br>
<br>
<br>
Thanks.<br>
<br>
<br>
<br>
<br>
<br>
</div>
</div>
</blockquote>
</div>
<br>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>