2,684
edits
Line 588: | Line 588: | ||
[[File:1y13-raddam-contrast-raddam.png]] | [[File:1y13-raddam-contrast-raddam.png]] | ||
== Automatically building | == Automatically building the main chain of 452 out of 519 residues == | ||
Based on the sites obtained by SHELXD, we run | Based on the sites obtained by SHELXD, we run | ||
shelxe.beta -a -q -h -b -s0.585 -m40 raddam raddam_fa | shelxe.beta -a -q -h -b -s0.585 -m40 raddam raddam_fa | ||
This already builds a significant number of residues, but also gives an improved list of heavy atom sites - there are actually 6 sites instead of the 5 that SHELXD wrote out (yes, we had asked SHELXD for 3 sites since there are 3 Met residues, but SHELXD as always was smarter than we are). We "mv raddam.hat raddam_fa.res" for another run of SHELXE: | This already builds a significant number of residues, but also gives an improved list of heavy atom sites - there are actually 6 sites instead of the 5 that SHELXD wrote out (yes, we had asked SHELXD for 3 sites since there are 3 Met residues, but SHELXD as always was smarter than we are). We "mv raddam.hat raddam_fa.res" for another run of SHELXE: | ||
shelxe.beta -a -q -h6 -b -s0.585 -m40 raddam raddam_fa | shelxe.beta -a -q -h6 -b -s0.585 -m40 -n3 raddam raddam_fa | ||
and get | and get | ||
<pre> | <pre> | ||
452 residues left after pruning, divided into chains as follows: | |||
A: | A: 15 B: 5 C: 22 D: 22 E: 27 F: 62 G: 263 H: 36 | ||
CC for partial structure against native data = | CC for partial structure against native data = 39.83 % | ||
------------------------------------------------------------------------------ | ------------------------------------------------------------------------------ | ||
Line 607: | Line 605: | ||
Global autotracing cycle 4 | Global autotracing cycle 4 | ||
<wt> = 0.300, Contrast = 0. | <wt> = 0.300, Contrast = 0.447, Connect. = 0.705 for dens.mod. cycle 1 | ||
<wt> = 0.300, Contrast = 0. | <wt> = 0.300, Contrast = 0.660, Connect. = 0.781 for dens.mod. cycle 2 | ||
<wt> = 0.300, Contrast = 0. | <wt> = 0.300, Contrast = 0.723, Connect. = 0.801 for dens.mod. cycle 3 | ||
<wt> = 0.300, Contrast = 0. | <wt> = 0.300, Contrast = 0.762, Connect. = 0.807 for dens.mod. cycle 4 | ||
Pseudo-free CC = | Pseudo-free CC = 64.88 % | ||
<wt> = 0.300, Contrast = 0. | <wt> = 0.300, Contrast = 0.785, Connect. = 0.810 for dens.mod. cycle 5 | ||
<wt> = 0.300, Contrast = 0. | <wt> = 0.300, Contrast = 0.806, Connect. = 0.813 for dens.mod. cycle 6 | ||
<wt> = 0.300, Contrast = 0. | <wt> = 0.300, Contrast = 0.820, Connect. = 0.815 for dens.mod. cycle 7 | ||
<wt> = 0.300, Contrast = 0. | <wt> = 0.300, Contrast = 0.831, Connect. = 0.817 for dens.mod. cycle 8 | ||
<wt> = 0.300, Contrast = 0. | <wt> = 0.300, Contrast = 0.839, Connect. = 0.819 for dens.mod. cycle 9 | ||
Pseudo-free CC = | Pseudo-free CC = 69.74 % | ||
<wt> = 0.300, Contrast = 0. | <wt> = 0.300, Contrast = 0.845, Connect. = 0.820 for dens.mod. cycle 10 | ||
<wt> = 0.300, Contrast = 0. | <wt> = 0.300, Contrast = 0.849, Connect. = 0.821 for dens.mod. cycle 11 | ||
<wt> = 0.300, Contrast = 0. | <wt> = 0.300, Contrast = 0.851, Connect. = 0.822 for dens.mod. cycle 12 | ||
<wt> = 0.300, Contrast = 0. | <wt> = 0.300, Contrast = 0.853, Connect. = 0.823 for dens.mod. cycle 13 | ||
<wt> = 0.300, Contrast = 0. | <wt> = 0.300, Contrast = 0.854, Connect. = 0.823 for dens.mod. cycle 14 | ||
Pseudo-free CC = 70. | Pseudo-free CC = 70.80 % | ||
<wt> = 0.300, Contrast = 0. | <wt> = 0.300, Contrast = 0.854, Connect. = 0.824 for dens.mod. cycle 15 | ||
<wt> = 0.300, Contrast = 0. | <wt> = 0.300, Contrast = 0.855, Connect. = 0.824 for dens.mod. cycle 16 | ||
<wt> = 0.300, Contrast = 0. | <wt> = 0.300, Contrast = 0.855, Connect. = 0.824 for dens.mod. cycle 17 | ||
<wt> = 0.300, Contrast = 0. | <wt> = 0.300, Contrast = 0.854, Connect. = 0.824 for dens.mod. cycle 18 | ||
<wt> = 0.300, Contrast = 0. | <wt> = 0.300, Contrast = 0.854, Connect. = 0.824 for dens.mod. cycle 19 | ||
Pseudo-free CC = | Pseudo-free CC = 71.03 % | ||
<wt> = 0.300, Contrast = 0. | <wt> = 0.300, Contrast = 0.854, Connect. = 0.824 for dens.mod. cycle 20 | ||
Estimated mean FOM and mapCC as a function of resolution | Estimated mean FOM and mapCC as a function of resolution | ||
d inf - 4.62 - 3.64 - 3.17 - 2.88 - 2.67 - 2.51 - 2.38 - 2.27 - 2.18 - 2.11 | d inf - 4.62 - 3.64 - 3.17 - 2.88 - 2.67 - 2.51 - 2.38 - 2.27 - 2.18 - 2.11 | ||
<FOM> 0. | <FOM> 0.736 0.786 0.768 0.721 0.701 0.681 0.618 0.595 0.587 0.540 | ||
<mapCC> 0. | <mapCC> 0.862 0.932 0.946 0.934 0.924 0.924 0.922 0.913 0.882 0.858 | ||
N 4206 4227 4214 4135 4185 4207 4292 4406 4320 3702 | N 4206 4227 4214 4135 4185 4207 4292 4406 4320 3702 | ||
Estimated mean FOM = 0. | Estimated mean FOM = 0.674 Pseudo-free CC = 71.18 % | ||
Density (in map sigma units) at input heavy atom sites | Density (in map sigma units) at input heavy atom sites | ||
Site x y z occ*Z density | Site x y z occ*Z density | ||
1 0. | 1 0.2276 0.7578 0.1189 34.0000 29.98 | ||
2 0. | 2 0.1568 0.6345 0.3049 32.2898 30.44 | ||
3 0. | 3 0.1767 0.5344 0.2160 32.2388 29.67 | ||
4 0.3059 0. | 4 0.3059 0.4535 0.1297 26.0746 23.51 | ||
5 0. | 5 0.0280 0.8243 0.1410 22.7324 21.02 | ||
6 0. | 6 0.0383 0.9748 0.0492 21.5050 21.18 | ||
Site x y z h(sig) near old near new | Site x y z h(sig) near old near new | ||
1 0. | 1 0.1569 0.6345 0.3048 30.4 2/0.02 9/13.36 3/15.73 2/19.52 7/22.13 | ||
2 0. | 2 0.2278 0.7578 0.1188 30.0 1/0.02 1/19.52 6/21.97 7/22.48 9/25.02 | ||
3 0.1767 0. | 3 0.1767 0.5345 0.2158 29.7 3/0.03 9/2.90 1/15.73 4/19.45 2/26.88 | ||
4 0. | 4 0.3060 0.4536 0.1292 23.5 4/0.07 3/19.45 9/21.16 8/26.49 5/26.83 | ||
5 0. | 5 0.0382 0.9748 0.0490 21.2 6/0.02 8/2.63 8/15.66 5/15.88 6/19.80 | ||
6 0. | 6 0.0278 0.8240 0.1416 21.1 5/0.08 5/19.80 8/21.59 7/21.87 2/21.97 | ||
7 0. | 7 0.1854 0.9571 0.1787 -5.0 5/21.86 6/21.87 1/22.13 2/22.48 8/22.57 | ||
8 0. | 8 0.0427 0.9993 0.0530 -5.0 6/2.62 5/2.63 8/15.31 5/15.66 6/21.59 | ||
9 0. | 9 0.1787 0.5611 0.2228 -4.7 3/2.91 3/2.90 1/13.36 4/21.16 2/25.02 | ||
</pre> | </pre> | ||
At this point the structure is obviously solved, and we could use buccanneer or Arp/wArp to add side chains and the rest of the model. 3-fold NCS surely helps! | At this point the structure is obviously solved, and we could use buccanneer or Arp/wArp to add side chains and the rest of the model. 3-fold NCS surely helps! | ||
== Could we do better? == | == Could we do better? == |