Can I solve my structure by SAD phasing? Anomalous signal in SAD phasing

Terwilliger, T.C.; Bunkóczi, G.; Hung, L.-W.; Zwart, P.H.; Smith, J.L.; Akey, D.L.; Adams, P.D.

doi:10.1107/S2059798315019269

research papers

STRUCTURAL
BIOLOGY

ISSN: 2059-7983

Volume 72| Part 3| March 2016| Pages 346-358

doi:10.1107/S2059798315019269

Open

access

Can I solve my structure by SAD phasing? Anomalous signal in SAD phasing

Thomas C. Terwilliger,^a ^* Gábor Bunkóczi,^b Li-Wei Hung,^c Peter H. Zwart,^d Janet L. Smith,^e,^f David L. Akey ^e and Paul D. Adams ^d

^aBioscience Division, Los Alamos National Laboratory, Mail Stop M888, Los Alamos, NM 87545, USA, ^bDepartment of Haematology, University of Cambridge, Cambridge Institute for Medical Research, Wellcome Trust/MRC Building, Cambridge CB2 0XY, England, ^cPhysics Division, Los Alamos National Laboratory, Mail Stop D454, Los Alamos, NM 87545, USA, ^dPhysical Biosciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA, ^eLife Sciences Institute, University of Michigan, Ann Arbor, MI 48109, USA, and ^fDepartment of Biological Chemistry, University of Michigan, Ann Arbor, MI 48109, USA
^*Correspondence e-mail: terwilliger@lanl.gov

(Received 6 April 2015; accepted 12 October 2015; online 1 March 2016)

A key challenge in the SAD phasing method is solving a structure when the anomalous signal-to-noise ratio is low. A simple theoretical framework for describing measurements of anomalous differences and the resulting useful anomalous correlation and anomalous signal in a SAD experiment is presented. Here, the useful anomalous correlation is defined as the correlation of anomalous differences with ideal anomalous differences from the anomalous substructure. The useful anomalous correlation reflects the accuracy of the data and the absence of minor sites. The useful anomalous correlation also reflects the information available for estimating crystallographic phases once the substructure has been determined. In contrast, the anomalous signal (the peak height in a model-phased anomalous difference Fourier at the coordinates of atoms in the anomalous substructure) reflects the information available about each site in the substructure and is related to the ability to find the substructure. A theoretical analysis shows that the expected value of the anomalous signal is the product of the useful anomalous correlation, the square root of the ratio of the number of unique reflections in the data set to the number of sites in the substructure, and a function that decreases with increasing values of the atomic displacement factor for the atoms in the substructure. This means that the ability to find the substructure in a SAD experiment is increased by high data quality and by a high ratio of reflections to sites in the substructure, and is decreased by high atomic displacement factors for the substructure.

Keywords: SAD phasing; anomalous signal; anomalous phasing; solving structures.

1. Introduction

1.1. Single-wavelength anomalous diffraction

The single-wavelength anomalous diffraction (SAD) method is a remarkably powerful approach to macromolecular structure determination (Hendrickson & Teeter, 1981 ; Wang, 1985 ; reviewed in Dauter et al., 2002 ; Hendrickson, 2014 ). It has become the dominant method for the determination by X-ray diffraction of structures that are not closely related to any structure already present in the Protein Data Bank (PDB; Berman et al., 2000 ), accounting for 73% of such new structures by 2013. In the SAD approach, the X-ray wavelength is tuned to be at or near an absorption edge of an element that is present at a limited number of sites in the macromolecule. Depending on the element and the wavelength, part of the scattering from the atoms in this substructure will then be shifted in phase by π/2 relative to the `normal' scattering from other atoms in the structure. This results in a difference in intensities for reflections related by inversion which otherwise would have identical intensities (Bijvoet, 1954 ). The anomalous differences, which correspond to the atoms in the substructure (Kartha & Parthasarathy, 1965 ; North, 1965 ; Strahs & Kraut, 1968 ), are then used to find the locations of these substructure atoms (Weeks et al., 1993 ; Terwilliger & Berendzen, 1999 ; Schneider & Sheldrick, 2002 ; Grosse-Kunstleve & Adams, 2003 ). The substructure and the anomalous differences are then used together to estimate crystallographic phases for the entire structure (Hendrickson & Teeter, 1981; Wang, 1985; Otwinowski, 1991 ; de La Fortelle & Bricogne, 1997 ; Furey & Swaminathan, 1997 ; McCoy et al., 2004 ; Pannu & Read, 2004 ). An electron-density map can then be calculated with these phases and the averages of structure factors for Bijvoet pairs. This electron-density map is normally then improved by density-modification techniques (Wang, 1985; Cowtan & Main, 1996 ; Abrahams & Leslie, 1996 ; Terwilliger, 2000 ; Cowtan, 2010 ) and interpreted in terms of an atomic model (see, for example, Jones et al., 1991 ; Perrakis et al., 1999 ; Emsley et al., 2010 ; Terwilliger et al., 2008 ; Cowtan, 2006 ; Langer et al., 2008 ).

There are several crucial steps in carrying out a SAD experiment (Dauter et al., 2002; Liu et al., 2011 ; Hendrickson, 2014). Some of these are experimental steps, such as having a crystal with well ordered anomalous scatterers and making accurate measurements of the intensities of Bijvoet pairs of reflections without substantial radiation damage (see, for example, Debreczeni et al., 2003 ; González, 2003 ; Garman, 2003 ; Krojer et al., 2013 ). Subsequent crucial steps involve analysis of the experimental data, decisions about which data to include, finding the locations of the atoms in the substructure, choosing the correct hand of the substructure and obtaining sufficiently accurate phase information to allow density-modification procedures to improve the phases and yield an interpretable electron-density map (Hendrickson, 2014). Here, we describe a formulation for the useful anomalous signal and test it with several solved structures. In the accompanying manuscript (Terwilliger et al., 2016 ), the application of the formulation to data sets for unsolved structures is described.

1.2. Measures of signal and noise in anomalous data

The anomalous differences between Bijvoet pairs of reflections are generally small (typically 1–5%), so a key overall consideration in a SAD experiment is obtaining sufficient signal. Additionally, anomalous differences, even if perfectly measured, are not the same as the structure factors for the anomalously scattering substructure, so these differences can be thought of as having a high level of intrinsic noise. Although these factors have been recognized for some time, it has not been fully clear which metrics best describe the signal in a SAD experiment and which values of these metrics indicate that the substructure will be solved and a sufficiently accurate electron-density map obtained.

1.2.1. Metrics predictable from experimental setup

A measure of signal in anomalous data that can be estimated in advance of carrying out an experiment is the Bijvoet ratio 〈|F⁺ − F⁻|〉/〈F〉 (Hendrickson & Teeter, 1981; Wang, 1985). As noted by Zwart (2005 ) and Dauter (2006 ), this measure is useful for obtaining a general idea of how large the anomalous signal might be, but the experiment and errors in measurement substantially affect the actual anomalous signal. Additionally, this ratio is strongly affected by the atomic displacement factors of atoms in the anomalous substructure (Shen et al., 2003 ; Zwart, 2005).

1.2.2. Metrics calculated from measured anomalous differences

An important set of metrics for signal in anomalous data are based on estimates of the accuracy of measured anomalous differences. One of these is the ratio of the mean anomalous difference to the mean estimated uncertainty in the difference 〈|Δ_ano|〉/〈σ_ano〉 (Schneider & Sheldrick, 2002; Wang, 1985). This can be used to identify the resolution to which the anomalous differences are useful (Schneider & Sheldrick, 2002). A related measure is based on a normal probability plot of normalized anomalous differences (Howell & Smith, 1992 ). Another related metric used to identify which anomalous differences are useful is the `measurability' of an anomalous data set (Parthasarathy & Parthasarathi, 1974 ; Zwart, 2005), which describes the fraction of anomalous differences that are very accurately measured (those with anomalous differences at least three times the magnitude of their uncertainties). A different metric used to identify whether anomalous differences are useful is a comparison of the merging χ² considering Bijvoet pairs as being equivalent with the χ² considering them separately (Otwinowski & Minor, 1997 ).

An estimate of signal that is based on experimental measurements but that does not require the use of experimental estimates of uncertainty is the anomalous scattering ratio, R_as (Fu et al., 2004 ). This measure is the ratio of differences among measurements of equivalent acentric reflections compared with measurements of equivalent centric reflections. As centric reflections have no anomalous differences but the errors in measurement are likely to be similar to those of the acentric reflections, this comparison can be a good indicator of the signal in the anomalous data. Its value for a series of data sets for zinc-free insulin crystals (Fu et al., 2004) was closely related to whether the substructure could be determined using SHELXD (Schneider & Sheldrick, 2002).

A very powerful measure of signal in an anomalous data set that also does not require estimates of experimental uncertainties is the correlation of anomalous differences measured at different wavelengths, obtained from different crystals, measured in different regions of reciprocal space, or simply measured more than once. The correlation of anomalous differences between data collected in different X-ray diffraction images obtained from crystals derivatized with heavy atoms was used some time ago to confirm the presence of anomalous differences using data collected on film (Colman et al., 1972 ; Buehner et al., 1974 ). The correlation of anomalous differences measured at different X-ray wavelengths of a MAD experiment were used in SOLVE (Terwilliger & Berendzen, 1999) and SHELXD (Schneider & Sheldrick, 2002) to identify the resolution to which the anomalous differences were likely to be useful (Dauter, 2006). More recently, the half-data-set anomalous correlation (Evans, 2006 ), obtained by dividing an unmerged data set into two parts and calculating the correlation of anomalous differences in the two parts, has been widely used to evaluate the utility of anomalous differences in SAD experiments.

1.2.3. Metrics requiring known structure

Even in the absence of measurement errors, the anomalous substructure does not fully account for the observed anomalous differences, since this also contains a contribution from atoms that are not detectable from the signal owing to their individual contributions being very weak (for example having a low anomalous scattering factor, low occupancy or high atomic displacement parameters) or from anomalously scattering atoms in the solvent continuum (Fourme et al., 1995 ). However, the cumulative effect of such atoms can be substantial in experiments such as sulfur SAD carried out at longer wavelengths (for example 1.8 Å), where the anomalous signal from C, N and O atoms becomes significant (Hendrickson, 2014), or in SAD experiments with heavy-atom soaks, where there are many minor sites, and even in selenomethionine SAD experiments, where the selenomethionine side chains might have multiple conformations.

In this work, we will define the `useful anomalous correlation' as the correlation between measured anomalous differences and the ideal anomalous differences that correspond to the final refined structure considering anomalous scattering from the detectable anomalous substructure in the crystal only. This correlation of course cannot be calculated directly unless the structure has been solved. Nevertheless, we will show here that it is a useful way of quantifying the information that is present in individual anomalous differences. The useful anomalous correlation can be thought of as describing the fraction of the measured anomalous differences that correspond to the ideal anomalous differences coming just from the substructure. It is also affected (typically decreased) by measurement errors or radiation damage. As the useful anomalous correlation reflects the information available for phase calculation, it is anticipated to be related to the quality of the electron-density map that can be obtained from a SAD experiment (cf. Zwart, 2005, where the correlation of the anomalous differences to the values of the true anomalous structure factor F_A is related to the ease of SAD structure determination).

A second important metric is the `anomalous signal', defined here, as in the work of others, as the peak height in an anomalous difference Fourier at the coordinates of the atoms in the anomalous substructure (Yang et al., 2003 ). The anomalous signal can also only be calculated after the structure has been solved. We will show here that it is a good measure of the total information present per site in the substructure , Jane S. Richardsonin all the anomalous differences in a data set. It can be calculated from the observed anomalous differences and an atomic model for the structure (without the need to consider anomalous scattering from the structure or substructure). The anomalous signal has been found to be closely related to whether the anomalous substructure can be determined (Yang et al., 2003; Liu et al., 2011, 2013 ; Akey et al., 2014 ; Bunkóczi et al., 2014 ; Weinert et al., 2015 ).

1.3. Objectives

In this work, we develop a simple framework for describing the relationships between measured anomalous differences, useful anomalous correlation and the anomalous signal. We show that the anomalous signal can be used to estimate the probability of determining the anomalous substructure and that the useful anomalous correlation can be used to estimate the expected quality of the electron-density map that can be calculated once the substructure is known.

2. Methods

2.1. Structure-factor relationships and anomalous differences in anomalous scattering

We represent the scattering factor (form factor) for the anomalously scattering atoms in the structure as

$[f^{\rm tot} = f^{\rm o} + f' + if'', \eqno(1a)]$

where f^o + f′ is the real part of the form factor and if′′ is the imaginary part. Fig. 1 illustrates the structure-factor relationships contributing to anomalous scattering for a particular reflection in the simple case where there is a single type of anomalous scatterer in the structure (for example sulfur or selenium). The structure factor for all of the nonscattering or weakly anomalously scattering atoms (carbon, nitrogen, oxygen etc.) is written here as F_P. The structure factor arising from the real part of the scattering factor f^o + f′ for the anomalously scattering atoms is F_H, and the part arising from the imaginary part of the scattering factor if′′ for these atoms is F_A. The structure factors (F⁺, F⁻) for this Bijvoet pair of reflections are then given (after reflection of the F⁻ member across the x axis; c.f. Kartha & Parthasarathy, 1965; Dauter et al., 2002) by

$[\eqalignno { {\bf F}^+ &= {\bf F}_{\rm P} + {\bf F}_{\rm H} + {\bf F}_{\rm A}, & (1b) \cr {\bf F}^{-*} &= {\bf F}_{\rm P} + {\bf F}_{\rm H} - {\bf F}_{\rm A}. & (1c)}]$

We will assume here that the structure factor arising from the imaginary part of the scattering factor for these atoms (F_A) is small relative to the other components. In this case, we can write an approximation for the difference in magnitudes between the members of this Bijvoet pair of reflections (F⁺ − F⁻). This `anomalous difference' (Δ_ano) is approximately given by

$[\Delta_{\rm ano} \equiv F^+ - F^- \simeq - 2F_{\rm A}\sin(\alpha), \eqno (1d)]$

where α is the angle between the structure factors F_P and F_H (Fig. 1).

Figure 1
Diagram of the relationships between structure factors and anomalous differences. Structure factors corresponding to an acentric reflection with indices (h, k, l) and to its Bijvoet mate with indices (−h, −k, −l) are given. The structure factors for the Bijvoet mate are reflected across the x axis for clarity in presentation. The structure factor for the non-anomalous atoms in the structure is designated as F_P. The part of the structure factor originating from the real part of the form factor for the atoms in the anomalous substructure (f^o + f′) is shown as F_H and the part of this structure factor coming from the imaginary part of the form factor (if′′) is shown as F_A.

2.2. Contributions to the anomalous signal

The anomalous signal (S_ano; Yang et al., 2003) is the mean peak height in an anomalous difference Fourier map at the coordinates of the anomalously scattering atoms, normalized to the r.m.s. value of the anomalous difference Fourier map. The anomalous difference Fourier map is calculated with coefficients based on the anomalous differences in (1d),

$[\rho_{\rm ano}(x) = {1 \over V}{\textstyle\sum\limits_h} \Delta_{{\rm ano},h}\exp\left[i(\varphi_h^{\rm c}) - {{\pi} \over 2}\right]\exp[-2\pi i(h \cdot x)], \eqno (2a)]$

where Δ_ano,h is the anomalous difference for this reflection for reflection h and φ_h^c is the phase of the structure factor for the non-anomalous part of the structure (F_P). The anomalous signal S_ano is then

$[S_{\rm ano} \equiv {{\langle \rho_{\rm ano}(x_j)\rangle } \over {\langle\rho_{\rm ano}^2\rangle^{1/2}}}, \eqno (2b)]$

where ρ_ano(x_j) is the value of the anomalous difference Fourier at the position of the jth anomalously scattering atom and 〈ρ²_ano〉^1/2 is the r.m.s. of the map.

An anomalous difference Fourier is typically calculated in order to show the positions of the atoms in the anomalous substructure. It is related to the ideal Fourier map for these atoms, except that the coefficients F_H,hexp[i(φ_h^c + α_h)] in the ideal Fourier map are replaced by Δ_ano,hexp[i(φ_h^c − π/2)] in the anomalous difference Fourier. We can see that this is reasonable by rearranging these coefficients slightly and by considering (1d). The coefficients in the ideal Fourier map are

$[F_{{\rm H},h}\exp[i(\varphi_h^{\rm c} + \alpha _h)] = F_{{\rm H},h}[\cos(\alpha_h) + i\sin(\alpha _h)]\exp(i\varphi_h^{\rm c}) \eqno (2c)]$

and the coefficients in the anomalous difference Fourier map (using 1d) are

$[\Delta_{{\rm ano},h}\exp\left[i\left(\varphi_h^{\rm c} - {\pi \over 2}\right)\right] = 2F_{{\rm A},h}i\sin(\alpha _h)\exp(i\varphi_h^{\rm c}),\eqno (2d)]$

where h = (h, k, l) are the indices of a reflection, Δ_ano,h is the anomalous difference for this reflection, φ_h^c is the phase of the structure factor for the non-anomalous part of the structure (F_P) calculated from the known structure, α_h is the angle between the structure factor arising from the real part of the form factor for the anomalous substructure (F_H) and the structure factor for the non-anomalous part of the structure factors (F_P; cf. Fig. 1). Inspecting (2c) and (2d), it can be seen that the anomalous difference Fourier represents just the sine term in (2c), and with a factor of 2F_A,h (twice the anomalous structure-factor amplitude) instead of F_H,h (the real structure-factor amplitude). It is therefore reasonable to expect that this map will have peaks at the positions of atoms in the anomalous substructure, but that the map will have a high level of intrinsic noise owing to missing the cosine term from (2c).

2.2.1. Anomalous signal for an ideal anomalous difference Fourier

We can now calculate how high the anomalous signal would be expected to be for an ideal anomalous difference Fourier (i.e. one with no measurement errors). Assuming that there is only one type of anomalous scatterer present, the structure factor F_A is perpendicular to F_H (Fig. 1) and their magnitudes are related for a particular reflection h by a factor a that corresponds to the ratio of anomalous to real scattering for the anomalous substructure and that depends on the resolution of the reflection,

$[\eqalignno {F_{\rm A} &= aF_{\rm H}, & (3a) \cr a &= {{f''}\over {f^{\rm o} + f'}}, & (3b)}]$

where f^o + f′ and f′′ are the real and imaginary parts of the scattering factor for the atoms in the anomalous substructure at the resolution of reflection h (1a). As the imaginary part of the scattering is normally owing to core electrons close to the nucleus, its real-space image can be adequately represented by a delta function whose Fourier transform is constant as a function of resolution. Consequently, the values of f′′ and f′ are typically assumed to be constant with resolution but the value of the form factor f^o falls off with resolution (Chantler, 1995 ; Hendrickson, 2014).

We can use these relationships and (2c) and (2d) to calculate the expected value of an anomalous difference Fourier map at the coordinates of the atoms in the substructure as well as the expected r.m.s. value of the map. Without loss of generality, in this analysis we will assume that the space-group setting chosen is a primitive setting (without centering). The value of the density in an anomalous difference Fourier map at coordinates represented by x = (x, y, z) can be written using (1d) and (2b) as

$[\eqalignno {\rho_{\rm ano}&(x) = {1 \over V}{\textstyle\sum\limits_h} \Delta_{{\rm ano},h} \exp \left[ i\left(\varphi_h^{\rm c} - {\pi \over 2}\right)\right] \exp[-2\pi i(h \cdot x)] \cr &\simeq {1 \over V}{\textstyle \sum\limits_h} -2F_{{\rm A},h}\sin(\alpha_h)\exp\left[i\left(\varphi_h^{\rm c} - {\pi \over 2} \right)\right] \exp[-2\pi i(h \cdot x)], \cr &&(4)}]$

where h, Δ_ano,h, φ_h^c and α_h are as in (2c) and (2d) and V is the volume of the unit cell. After rearrangement and using the scale factor a given above between the anomalous structure-factor amplitude (F_A) and the real part of the structure-factor amplitude for the substructure (F_H), this expression for the density in an anomalous difference Fourier becomes

$[\eqalignno {\rho_{\rm ano}(x) &\simeq {1 \over V}{\textstyle\sum\limits_h} {{f''_{h}}\over{f^{\rm o}_h + f'_{h}}} F_{{\rm H},h}\exp[i(\varphi_{h}^{\rm c} + \alpha_{h})]\exp[-2\pi i(h\cdot x)] \cr &\ \quad {\times}\ [1- \exp (-2i \alpha_{h})].&(5)}]$

The structure factor arising from the real part of the form factor for the anomalous substructure (F_H,h; Fig. 1) can be calculated from the scattering factors (f_h^o + f_h′), coordinates (x_j) and atomic displacement factors (B_j) of the substructure as

$[\eqalignno {{\bf F}_{{\rm H},h} & \equiv F_{{\rm H},h}\exp[i(\varphi_h^{\rm c} + \alpha_h)] \cr & = {\textstyle\sum\limits_j} (f_h^{\rm o} + f'_h)\exp[-B_j(\sin^2\theta_h/\lambda^2)]\exp[2\pi i(h \cdot x_j)], &(6)}]$

where sin(θ_h) is the scattering angle for reflection h and λ is the wavelength of the X-rays used in the experiment. Substituting this expression for F_H,h into the approximation for the density in an anomalous difference Fourier gives

$[\eqalignno {\rho_{\rm ano}(x) & \simeq {1 \over V} {\textstyle\sum\limits_h} f''_{h}\exp[-2\pi i (h \cdot x)][1-\exp(-2i\alpha_{h})]\cr &\ \quad {\times}\ {\textstyle\sum\limits_j}\exp[-B_{j}(\sin^{2}\theta_{h}/\lambda^{2})]\exp[2\pi i (h \cdot x_{j})]. & (7)}]$

If all the atoms in the anomalous substructure are assumed to have identical atomic displacement factors (B_j = B), then this can be simplified slightly and rearranged to read

$[\rho_{\rm ano}(x) \simeq {1 \over V}{\textstyle\sum\limits_h} f_{h,B}[1 - \exp(-2i\alpha_h)] {\textstyle \sum\limits_j} \exp\{-2\pi i[h \cdot (x - x_j)]\}, \eqno(8)]$

where the factors f_h,B are the anomalous scattering factors adjusted for the effects of the atomic displacement factor B at the resolution of reflection h and are given by

$[f_{h,B} = f''_{h}\exp[-B(\sin^{2}\theta_{h}/\lambda^{2})].\eqno (9)]$

Noting that the expected value of the sum over all reflections $[\textstyle\sum_h \exp\{-2\pi i[h\cdot(x - x_j)]\}]$ is zero unless x is approximately equal to one of the coordinates x_j of an atom in the anomalous substructure, and noting that the phase angle α_h is the phase difference between the structure factor for the non-anomalous part of the structure and the anomalous substructure and is therefore essentially independent of the anomalous substructure, the expected value of the map 〈ρ_ano(x_j)〉 at atomic positions corresponding to the substructure is given by

$[\langle\rho_{\rm ano}(x_j)\rangle \simeq {N \over V}\langle f_{h,B}\rangle, \eqno(10)]$

where N is the number of reflections (including the entire sphere, not just the asymmetric unit of the structure factors) used to calculate the map, V is the volume of the unit cell and 〈f_h,B〉 is the mean value of the anomalous scattering factors (9) over all of the reflections included in the analysis.

A similar calculation can be made to estimate the mean-square value of the anomalous difference Fourier map and from it the r.m.s. of the map. The mean-square value of the density in this map is given by

$[\eqalignno {\langle\rho_{\rm ano}^2\rangle &= {1 \over V}{\textstyle\int\limits_V} {{{\rm d}V} \over V}{\textstyle\sum\limits_h \sum\limits_j} f_j^h\exp[2\pi i(h \cdot x_j)] \exp[-2\pi i(h \cdot x)]\cr &\ \quad {\times}\ [1 - \exp(-2i\alpha_h)] \cr &\ \quad {\times}\ {1 \over V}\textstyle \sum\limits_{h'} \sum\limits_{j'} f_{j'}^{h'}\exp[-2\pi i(h' \cdot x_{j'})] \exp[2\pi i(h' \cdot x)] \cr &\ \quad {\times}\ [1 - \exp(2i\alpha_{h'})]. & (11)}]$

The expected value of the mean-square density in the anomalous difference Fourier map is then given by

$[\langle \rho_{\rm ano}^2\rangle \simeq {{2Nn} \over {V^2}}\langle f_{h,B}^2\rangle, \eqno(12)]$

where, as before, N is the number of reflections and V is the volume of the unit cell. The number of atoms in the anomalous substructure in the entire unit cell is n, and 〈f²_h,B〉 is the mean-square value of the imaginary component of the anomalous scattering factors, including atomic displacement factors (9), over all of the reflections included in the analysis.

The anomalous signal is the ratio of the expected value of the density at coordinates of atoms in the anomalous substructure to the r.m.s. value of the map for the model-phased anomalous difference Fourier map. This can now be estimated,

$[\langle S_{\rm ano}^{\rm ideal} \rangle \simeq \left ({{N} \over {2nf_B}}\right)^{1/2}, \eqno(13)]$

where N and n refer to the number of reflections including the entire sphere and the number of sites in the entire unit cell, respectively, and where the factor f_B is the second moment of the values of the scattering factors,

$[f_B = {{\langle f_{h,B}^2\rangle } \over {\langle f_{h,B}\rangle^2}}. \eqno(14)]$

Writing the total number of reflections as the number of symmetry operators N_sym times the number of unique acentric reflections N_refl times two (for Bijvoet pairs of reflections),

$[N = 2N_{\rm refl}N_{\rm sym} \eqno(15)]$

and writing the number of atoms in the anomalous substructure for the entire crystal as the number of atoms in the substructure in the asymmetric unit times the number of symmetry operators,

$[n = n_{\rm sites}N_{\rm sym}, \eqno(16)]$

we can rewrite the expected anomalous signal in an ideal anomalous difference Fourier map as

$[\langle S_{\rm ano}^{\rm ideal}\rangle \simeq \left ({{N_{\rm refl}} \over {n_{\rm sites}f_B}}\right)^{1/2}, \eqno(17)]$

where N_refl is the number of unique noncentrosymmetric reflections and n_sites is the number of unique sites in the substructure. As noted at the beginning of the section, this analysis assumes that the space-group setting chosen is chosen to be primitive. As the number of unique reflections and the number of unique sites do not depend on the setting chosen, (17) can be applied equally well to centered space groups. Sites that are at special positions do require special treatment in (16), however (a site on a twofold will count as half the amount of a site on a general position).

2.2.2. Anomalous signal for a realistic anomalous difference Fourier

In this section, we will expand the calculation of expected anomalous signal to realistic cases where errors can be present in the measurement and where there may be significant anomalous scattering from minor sites. We can write a simple expression for contributions to an observed anomalous difference (Δ_ano^obs),

$[\Delta^{\rm obs}_{\rm ano} = \Delta_{\rm ano} + \Delta^{\rm other}_{\rm ano} + \varepsilon. \eqno(18)]$

Here, Δ_ano is the `useful' anomalous difference owing to the anomalous substructure, Δ_ano^other is the (true, but not useful) anomalous difference owing to all other anomalously scattering atoms and minor sites for the anomalous substructure and ∊ is the error in measurement (including the effects of radiation damage). From this expression it can be seen that the anomalous contributions (Δ_ano^obs) from sites that are not part of the anomalous substructure and errors in measurement have similar overall effects on the observed anomalous differences.

We will assume that the observed anomalous difference and each term contributing to it have expected values of zero. Further, we will assume that each term contributing to the anomalous difference is uncorrelated with the other terms. We define a normalized variance, E², as the ratio of the sum of the variances of the anomalous contribution from sites not part of the substructure Δ_ano^other and the errors in measurement ∊ to the mean-square value of the useful anomalous difference Δ_ano,

$[E^2 \equiv {{\langle(\Delta_{\rm ano}^{\rm other})^2\rangle + \langle(\varepsilon_{\rm ano})^2\rangle } \over {\langle \Delta_{\rm ano}^2\rangle}}. \eqno(19)]$

With this definition, we are in a position to calculate the expected value of the anomalous signal in the case where minor sites and measurement errors are present. The density in such a realistic anomalous difference Fourier map can be written (see equation 4) as

$[\eqalignno {\rho_{\rm ano}^{\rm obs}(x) &= {1 \over V}{\textstyle\sum\limits_h} \Delta^{\rm obs}_{{\rm ano},h} \exp\left[i\left(\varphi_h^{\rm c} - {\pi \over 2}\right)\right] \exp[-2\pi i(h \cdot x)] \cr & \simeq {1 \over V}{\textstyle \sum\limits_h} [-2F_{{\rm A},h}\sin(\alpha_h) + \Delta^{\rm other}_{{\rm ano},h} + \varepsilon _h] \cr &\ \quad{\times}\ \exp\left[i\left(\varphi_h^{\rm c} - {\pi \over 2}\right)\right] \exp[-2\pi i(h \cdot x)]. &(20)}]$

As in the previous section, this can be simplified and rearranged, yielding an expression for the density in an anomalous difference Fourier that is the sum of the density for a perfect anomalous difference Fourier and a term that contains the contributions from minor sites and errors in measurement,

$[\eqalignno {\rho_{\rm ano}^{\rm obs}(x) &\simeq {1 \over V}{\textstyle\sum\limits_h} f_h [1 - \exp(-2i\alpha_h)] {\textstyle \sum\limits_j} \exp\{-2\pi i[h \cdot (x - x_j)]\} \cr &\ +\ {1 \over V}{\textstyle \sum\limits_h} (\Delta^{\rm other}_{{\rm ano},h} + \varepsilon_h) \exp\left[i\left(\varphi_h^{\rm c} - {\pi \over 2}\right)\right]\exp[-2\pi i(h \cdot x)]. \cr &&(21)}]$

As the expected value of each of the two error terms is zero, the expected value of this density at the coordinates of atoms in the anomalous substructure is the same as in the case without errors (8), with the expected value given in (10). The expected value of the mean-square density in this map, however, is now higher owing to the error terms. It is given (compare with equation 12) by

$[\langle \rho_{\rm ano}^2\rangle \simeq {{2Nn} \over {V^2}} \langle f_{h,B}^2\rangle + {N \over {V^2}} [\langle (\Delta_{\rm ano}^{\rm other})^2\rangle + \langle\varepsilon_{\rm ano}^2\rangle]. \eqno(22)]$

We can simplify this in several steps. From (1d), it may be seen that the expected mean-square value of the anomalous difference corresponding to the substructure atoms (〈Δ²_ano〉) is given by

$[\langle \Delta^2_{\rm ano}\rangle = 4\langle F_{\rm A}^2\rangle\langle\sin^2(\alpha)\rangle = 2\langle F_{\rm A}^2 \rangle. \eqno(23)]$

The value of 〈F_A²〉 on the right-hand side of (23) can in turn be calculated in two steps. Firstly, substituting (3a) and (3b) into (6), assuming again that the atoms in the anomalous substructure are assumed to have identical atomic displacement factors (B_j = B), and then substituting in (9) yields an expression relating the magnitude F_A to the anomalous scattering factor f_h,B from (9) and the coordinates of the anomalously scattering atoms x_j,

$[\eqalignno {F_{{\rm A},h}\exp[i(\varphi_h^{\rm c} + \alpha_h)] & = \textstyle\sum\limits_j f''_h \exp[-B(\sin^2\theta_h/\lambda ^2)]\exp[2\pi i(h \cdot x_j)] \cr & = \textstyle\sum\limits_j f_{h,B}\exp[2\pi i(h \cdot x_j)]. & (24)}]$

The mean-square value 〈F_A²〉 can then be calculated as

$[\eqalignno{\langle F_{\rm A}^2\rangle &= \textstyle\sum\limits_j f_{h,B}\exp[2\pi i(h \cdot x_j)] \textstyle\sum\limits_k f_{h,B}\exp[-2\pi i(h \cdot x_k)] \cr &= n\langle f_{h,B}^2\rangle, &(25)}]$

where n is the number of sites in the asymmetric unit of the crystal. Then, substituting (25) into (23), we have

$[\langle \Delta^2_{\rm ano}\rangle = 2n\langle f_{h,B}^2\rangle. \eqno(26)]$

Finally, using (26) along with the definition of E² in (19) and rearranging, we can simplify the estimate of the mean-square density in the anomalous difference Fourier map with errors in (22) to read

$[\langle \rho_{\rm ano}^2\rangle \simeq {{2Nn} \over {V^2}}\langle f_{h,B}^2\rangle (1 + {E^2}). \eqno(27)]$

Comparing (12) and (27), we see that the mean-square value of the density overall in the anomalous difference Fourier map is larger than that in a perfect anomalous difference Fourier map by the factor 1 + E². As noted above (see equations 8, 10 and 21), the mean value of the density in the anomalous difference Fourier map with errors at coordinates of atoms in the anomalous substructure is the same as that in a perfect map. This leads to an expression for the anomalous signal in the presence of errors in measurement and minor sites,

$[\langle S_{\rm ano}^{\rm obs}\rangle \simeq \left[{{N_{\rm refl}} \over {(1 + E^2)n_{\rm sites}f_B}}\right]^{1/2}, \eqno(28)]$

where the second moment of the values of the scattering factors f_B is given in (14). We can relate the error term 1 + E² to the correlation (CC_ano) between the measured anomalous differences and the anomalous differences owing to the substructure atoms. The useful anomalous correlation CC_ano is given by

$[{\rm CC}_{\rm ano} \equiv {{\langle \Delta_{\rm ano}\Delta_{\rm ano}^{\rm obs}\rangle}\over {\langle\Delta_{\rm ano}^2\rangle^{1/2}\langle(\Delta_{\rm ano}^{\rm obs})^2\rangle^{1/2}}}. \eqno(29)]$

Using (18) and (19), it may be shown that the expected value of the useful anomalous correlation is given by

$[{\rm CC}_{\rm ano} \simeq {1 \over {(1 + E^2)^{1/2}}}. \eqno(30)]$

This yields a simple formula for the expected anomalous signal in the presence of minor sites and errors in measurement,

$[\langle S_{\rm ano}^{\rm obs}\rangle \simeq {\rm CC}_{\rm ano}\left({{N_{\rm refl}} \over {n_{\rm sites}f_B}}\right)^{1/2}. \eqno(31)]$

Comparison with (17) shows that the expected anomalous signal in a realistic case is simply equal to its expected value in an ideal case multiplied by the useful anomalous correlation CC_ano.

2.3. Test data from the PDB

We downloaded data sets from the PDB to serve as test cases for our analyses. The data consisted of 218 MAD and SAD data sets from 113 different PDB entries with diffraction data extending to resolutions from 1.2 to 4.5 Å. The MAD PDB entries were split into individual data sets for each wavelength of X-rays used to measure diffraction data. The PHENIX (Adams et al., 2010 ) tool phenix.sad_data_from_pdb was used to extract the individual data sets from PDB entries. The PDB entries used were 1vjn , 1vjr , 1vjz , 1vk4 , 1vkm (Levin et al., 2005 ), 1vlm , 1vqr (Xu et al., 2006 ), 1xri (Aceti et al., 2008 ), 1y7e , 1z82 , 1zy9 , 1zyb , 2a2o , 2a3n , 2a6b , 2aj2 , 2aml , 2avn , 2b8m , 2etd , 2etj , 2ets (Kozbial et al., 2008 ), 2etv , 2evr (Xu, Sudek et al., 2009 ), 2f4p , 2fea (Xu et al., 2007 ), 2ffj , 2fg0 (Xu, Sudek et al., 2009), 2fg9 , 2fna (Xu, Rife et al., 2009 ), 2fqp , 2fur , 2fzt , 2g0t , 2g42 , 2gc9 , 2nlv (Hwang et al., 2014 ), 2nuj , 2nwv , 2o08 , 2o1q , 2o2x , 2o2z , 2o3l , 2o62 , 2o7t , 2o8q , 2obp , 2oc5 , 2od5 , 2od6 , 2oh3 , 2okc , 2okf (Hwang et al., 2014), 2ooj , 2opk , 2osd , 2otm , 2ozg , 2ozj , 2p10 , 2p4o , 2p7h , 2p7i , 2p97 , 2pg3 , 2pg4 , 2pgc , 2pim , 2pn1 , 2pnk , 2ppv , 2pr1 , 2pr7 , 2prv , 2pv4 , 2pw4 , 2wcd (Mueller et al., 2009 ), 2xdd (Fineran et al., 2009 ), 2zxh (Osawa et al., 2009 ), 3caz , 3din (Zimmer et al., 2008 ), 3dto , 3fx0 (Lo et al., 2009 ), 3guw , 3gw7 , 3hxk , 3hxp , 3lml , 3mv3 (Hsia & Hoelz, 2010 ), 3ov0 (Pokkuluri et al., 2011 ), 3pg5 , 3zgx (Bürmann et al., 2013 ), 3zxu (Schmitzberger & Harrison, 2012 ), 4acb (Leibundgut et al., 2004 ), 4asn (Aylett & Lowe, 2012 ), 4bql (Lindås et al., 2014 ), 4cb0 (Mechaly et al., 2014 ), 4cbv (Boudes et al., 2014 ), 4fsx (Du et al., 2012 ), 4g9i (Tominaga et al., 2012 ), 4gkw (Qiao et al., 2012 ), 4h6y (He et al., 2013 ), 4hkr (Hou et al., 2012 ), 4hnd (Zhou et al., 2014 ), 4ifq (Sampathkumar et al., 2013 ), 4lck (Zhang & Ferré-D'Amaré, 2013 ), 4nsc (Wang et al., 2014 ), 4nt5 (Zhou & Springer, 2014 ), 4px7 (Fan et al., 2014 ), 4q8j (Schäfer et al., 2014 ), 4qka (Gao & Serganov, 2014 ) and 4tq5 (Huang et al., 2014 ).

2.4. Software availability

We have developed software as part of the PHENIX software suite (Adams et al., 2010) that calculates the anomalous signal and anomalous correlation in a data set (phenix.anomalous_signal).

3. Results and discussion

3.1. Dependence of the anomalous signal S_ano^obs on CC_ano, N_refl, n_sites and f_B

The key theoretical result from this work (equation 31) is that the expected anomalous signal S_ano^obs for a SAD data set is related in a very simple way to the useful anomalous correlation CC_ano, the number of unique reflections N_refl, the number of sites in the anomalous substructure n_sites and the second moment of the scattering factors for the anomalous substructure f_B,

$[\langle S_{\rm ano}^{\rm obs}\rangle \simeq {\rm CC}_{\rm ano}\left({{N_{\rm refl}} \over {n_{\rm sites}f_B}}\right)^{1/2}. \eqno(31)]$

In (31) the anomalous signal S_ano^obs is the mean value of a model-phased anomalous difference Fourier at the coordinates of atoms in the anomalous substructure. The anomalous signal is the principal metric in this work for the useful signal in a SAD data set. The useful anomalous correlation CC_ano is the correlation between measured anomalous differences and those expected of an ideal structure where the only anomalous scatterers are those in the anomalous substructure and where there are no errors in measurement (29). The useful anomalous correlation CC_ano is the principal metric in this work for similarity between measured and ideal anomalous differences. The number of reflections N_refl includes all acentric reflections that are unique under space-group symmetry, and the number of sites n_sites is the number contained in the asymmetric unit of the crystal. The factor f_B (14) is the second moment of the scattering factors corresponding to the anomalous substructure. It is related to the fall-off with resolution of the scattering from the anomalous substructure because it includes the atomic displacement factors (assumed to all be equal; B_j = B; equation 9).

3.1.1. Anomalous signal for an idealized crystal

The significance of (31) is that it shows how the anomalous signal depends on CC_ano, N_refl, n_sites and f_B. If the data were measured perfectly and if there were no anomalous scatterers other than those in the substructure, then the useful anomalous correlation CC_ano would be unity. Further, if all of the anomalous scatterers had atomic displacement factors of zero (B = 0), then the second moment of the scattering factors f_B would also be unity. The expected value 〈S_ano^B=0〉 of the anomalous signal in this ideal case with atomic displacement factors of zero is then simply equal to the square root of the ratio of unique reflections to unique sites,

$[\langle S_{\rm ano}^{B = 0}\rangle \simeq \left({{N_{\rm refl}} \over {n_{\rm sites}}}\right)^{1/2}. \eqno(32)]$

(32) indicates that anomalous signal increases with the number of unique reflections and decreases with the number of unique sites in the anomalous substructure in a simple way. In particular, if there are more sites in the substructure then correspondingly more reflections must be measured to obtain the same anomalous signal.

In this ideal case with atomic displacement factors of zero, the maximum possible expected anomalous signal 〈S_ano^max〉 for a crystal is obtained if there is a single site in the anomalous substructure (32). The value of the maximum possible expected anomalous signal is just the square root of the number of reflections,

$[\langle S_{\rm ano}^{\rm max}\rangle \simeq N_{\rm refl}^{1/2}. \eqno(33)]$

3.1.2. Dependence of the anomalous signal for an idealized crystal on resolution

Up to this point, the influence of the resolution of the data on the anomalous signal has not been considered. If the anomalous scatterers in an otherwise idealized crystal have nonzero atomic displacement factors, then the contributions from these anomalous scatterers will be resolution-dependent (simply owing to the fall-off of intensities with resolution). Correspondingly, the anomalous signal (the peak heights at coordinates of anomalously scattering atoms in the anomalous difference Fourier) will not increase as the square root of the number of reflections as in (33), but rather by some smaller factor. This relationship is given in (17), where the effect of the fall-off of intensities with resolution is captured by the factor f_B, the value of the second moment of the scattering factors including atomic displacement factors. If the anomalously scattering atoms have nonzero atomic displacement factors, the value of the second moment of the scattering factors, f_B (14), will be greater than one. In such an idealized case, the anomalous signal 〈S_ano^ideal〉 will be depend on the resolution of the data through the factor f_B in (17),

$[\langle S_{\rm ano}^{\rm ideal}\rangle \simeq \left ({{N_{\rm refl}} \over {n_{\rm sites}f_B}}\right)^{1/2}. \eqno(17)]$

The second moment of the scattering factors f_B depends on how much the scattering factors vary over the resolution of the measured reflections (cf. equation 14, noting that a distribution with many large and many small values will have a large value of the second moment while a narrow distribution will have a small second moment).

3.1.3. Anomalous signal for anomalous data measured with errors from a real crystal

For a real crystal, some of the true anomalous differences will come from minor sites of the principal anomalously scattering atom and from the weak anomalous scattering from atoms not considered to be part of the anomalous substructure. Furthermore, the anomalous differences will be measured with experimental errors. In this case the measured anomalous differences will have a correlation CC_ano to those measured perfectly from an ideal crystal. The consequence of a useful anomalous correlation of less than unity is that the anomalous signal will be reduced based on this correlation,

$[\langle S_{\rm ano}^{\rm obs}\rangle \simeq {\rm CC}_{\rm ano}\left({{N_{\rm refl}} \over {n_{\rm sites}f_B}}\right)^{1/2}. \eqno(31)]$

In this realistic case, the overall expected anomalous signal is just the product of the useful anomalous correlation and the square root of the number of reflections (N_refl) divided by the number of sites in the anomalous substructure and the second moment f_B of the scattering factors for the anomalous substructure.

3.1.4. Anomalous signal and resolution of the data included in the calculation

In this analysis (31), the anomalous correlation CC_ano is the overall average for the entire data set. As the accuracies of anomalous differences typically decrease at high resolution, if additional high-resolution, low-accuracy anomalous differences are included in a data set, then the overall value of CC_ano will decrease. Inspecting (31), it can be seen that as the resolution of data included in the calculation increases, the net effect on the anomalous signal is a combination of three factors. These are (i) an increase owing to the number of reflections N_refl, (ii) a decrease owing to the decrease in average anomalous correlation CC_ano and (iii) a decrease owing to the increase in the second moment f_B of the scattering factors for the anomalous substructure, averaged over all reflections included. The combined effects of increasing the number of reflections and increasing the average second moment of the scattering factors for the anomalous substructure, as discussed above, is that adding reflections that have negligible contributions from the anomalous substructure (e.g. at resolutions where the structure factors are essentially zero) has no effect on the overall anomalous signal.

The net effect of increasing the number of reflections and decreasing the anomalous correlation, however, can be either an increase or a decrease in the anomalous signal. If high-resolution data were measured with accuracy similar to that of lower resolution data, then the anomalous correlation CC_ano would be approximately constant as additional data are included and the anomalous signal would increase as the square root of the number of reflections. In contrast, if high-resolution data were measured with poor accuracy and the anomalous differences at these resolutions were essentially random, then the additional data would contribute only noise to the anomalous difference Fourier map (2b). This would neither increase nor decrease the average peak height at positions of anomalously scattering atoms (the numerator in equation 2b), but it would increase the overall r.m.s. of the map (the denominator in equation 2b). Therefore, as expected, the inclusion of essentially random anomalous differences will decrease the overall anomalous signal S_ano. This behavior can also be seen if (29) and (31) are examined. The anomalous correlation is defined in (29),

$[{\rm CC}_{\rm ano} \equiv {{\langle \Delta_{\rm ano}\Delta_{\rm ano}^{\rm obs}\rangle}\over {\langle\Delta_{\rm ano}^2\rangle^{1/2}\langle(\Delta_{\rm ano}^{\rm obs})^2\rangle^{1/2}}}. \eqno(29)]$

Imagine adding random anomalous differences, and suppose that they are about the same in magnitude as those that are already in the data set. In this case, the value of the r.m.s. observed anomalous difference 〈(Δ_ano^obs)²〉^1/2 will not change. In contrast, the mean value of the product of the true and observed anomalous difference, 〈Δ_anoΔ_ano^obs〉, will decrease. If n random anomalous differences are added to N well measured anomalous differences its value will decrease by a factor of N/(N + n). This means again that adding in anomalous differences that are random will reduce the overall anomalous signal S_ano. We emphasize again that all of these effects are captured in (31), where the values of the anomalous correlation, number of reflections and second moment of scattering factors of the anomalously scattering atoms are all the overall values for all reflections considered.

3.1.5. Comparison of expected and actual anomalous signal

We carried out a comparison of the anomalous signal expected from (31) with the actual anomalous signal in data sets downloaded from the PDB. We used 113 MAD and SAD data sets from the PDB to serve as test cases for evaluating our theoretical analysis. The MAD data sets were split into separate `data sets', each containing the data measured at one X-ray wavelength, yielding a total of 218 data sets for analysis. The high-resolution limits of these data sets ranged from 1.2 to 4.5 Å and the anomalously scattering atoms included selenium, sulfur, cobalt, mercury, zinc, nickel, iron, calcium, barium and iridium. Each data set was used to the full resolution limit except as noted below.

As all of the structures are known, we could calculate ideal anomalous differences Δ_ano based on the anomalous substructure in the model. We used these ideal differences along with the measured anomalous differences Δ_ano^obs to calculate the useful anomalous correlation (CC_ano; equation 17). From the X-ray data and model, we could identify the number of sites in the anomalous substructure (n_sites) and the number of unique measured reflections (N_refl). We calculated the second moment of the scattering factors of atoms in the substructure (f_B) using (9) and (14) and using the mean value of the atomic displacement factors for the atoms in the substructure in (9). These calculations allowed us to estimate the value of the anomalous signal S_ano^obs using (31). Note that all of the data are included in the calculation of the anomalous signal S_ano^obs. The resolution dependence of the anomalous signal is captured in the resolution dependence of the number of reflections, the anomalous correlation and the second moment of scattering factors of atoms in the substructure. We then obtained the actual value of the anomalous signal from the mean peak height at coordinates of atoms in the substructure in an anomalous difference Fourier map (2b). Fig. 2 shows that the values of the anomalous signal obtained from (31) are very similar to the actual values obtained from an anomalous difference Fourier map using (2b), indicating that our simple theoretical analysis gives a good description of the overall contributions to the anomalous signal.

Figure 2
Comparison of anomalous signal estimated with (31)

and the actual anomalous signal. Each point in the figure corresponds to data measured at one X-ray wavelength taken from one entry in the PDB. The x coordinate of each point is the anomalous signal estimated as described in the text using (31)

and the y coordinate is the anomalous signal calculated directly from a model-phased anomalous difference Fourier map using (2b)

3.2. Relationship between the anomalous signal and the solution of the anomalous substructure

As noted above, the anomalous signal has been found to be a good indicator of whether the anomalous substructure can be obtained in a SAD experiment (Yang et al., 2003; Liu et al., 2011, 2013; Akey et al., 2014; Bunkóczi et al., 2014; Weinert et al., 2015). Fig. 3 illustrates this relationship for 1874 complete and truncated data sets extracted from the 218 SAD data sets used above. These data sets were constructed by truncating the data at resolutions ranging from 1.5 to 6 Å. For each complete or truncated data set, the anomalous signal was calculated using an anomalous difference Fourier map and the coordinates of the atoms in the anomalous substructure. The same data were then used as input to the HySS substructure-search tool using likelihood-based substructure completion (Bunkóczi et al., 2014), and the fraction of the true sites in the substructure is plotted in Fig. 3 as a function of the anomalous signal of the corresponding data set. Fig. 3 also illustrates the mean fraction of complete and truncated SAD data sets that are solved as a function of the anomalous signal. Those cases where at least 50% of the sites were found (sites placed within 3 Å of an atom or a symmetry-equivalent atom in the known substructure) were considered to be `solved' for this analysis.

Figure 3
Success in substructure determination as a function of the anomalous signal in SAD data sets. Each point in the plot represents the fraction of the anomalous substructure found using likelihood-based methods for a complete or a resolution-truncated SAD data set as described in the text. The line represents the fraction of data sets where at least 50% of the sites in the substructure were found, as calculated in bins of resolution.

It can be seen from Fig. 3 that for data sets where the anomalous signal is less than about 7–10 few of the SAD data sets could be solved with likelihood-based methods, while for those where the anomalous signal is greater than about 10–15 most of the data sets could be solved, supporting the earlier observations that the anomalous signal is a good indicator of whether a SAD data set can be solved. The fraction of data sets that could be solved (the solid line in Fig. 3) reaches about 50% at an anomalous signal of about 9.

3.3. Relationship between useful anomalous correlation and the quality of phase calculation

The useful anomalous correlation (CC_ano^obs; equation 29) is a measure of how well the measured anomalous differences reflect those expected from the anomalous substructure in the crystal. As discussed by Zwart (2005), the accuracy of crystallographic phases that can be calculated from a set of anomalous differences using knowledge of the anomalous substructure is related to the useful anomalous correlation. Fig. 4 illustrates the relationship between useful anomalous correlation and the map correlation to a model-phased map obtained after phase calculations are carried out using Phaser (McCoy et al., 2004) with the known anomalous substructure and the measured anomalous differences for the data sets shown in Fig. 2. Fig. 4 shows that for these data sets the map quality increases substantially with useful anomalous correlation.

Figure 4
Phase accuracy as a function of useful anomalous correlation. Each point represents one SAD data set. The useful anomalous correlation is calculated from the correlation of model-based and measured anomalous differences. The measured anomalous differences are used along with the known anomalous substructure to calculate crystallographic phases. The phase accuracy is represented as the correlation between a map calculated using these phases and a map calculated using model phases.

3.4. Relationship between number of sites in the anomalous substructure and the anomalous signal

(31) indicates that the anomalous signal is proportional to the inverse of the square root of the number of sites in the substructure. It is appropriate to consider how structure determination would be affected in a realistic situation by varying the number of sites. This is not quite as straightforward as it might seem, as increasing the number of sites in the substructure would normally also increase the anomalous differences. Consequently, if all other crystal-size and data-collection conditions were the same, the crystal with more sites in the substructure would have a higher useful anomalous correlation (CC_ano^obs; equation 29), partially offsetting the decrease in anomalous signal expected from (31). Two limiting situations can be considered. In the case where anomalous differences are measured very accurately and the useful anomalous correlation is already very near unity, increasing the number of sites cannot increase the useful anomalous correlation (CC_ano^obs) further, so the increase in the number of sites will decrease the anomalous signal according to the square root of the number of sites as in (31).

In the more common case where the useful anomalous correlation (CC_ano^obs) is lower, increasing the number of sites will have a smaller effect on the anomalous signal than would be found with a constant value of the anomalous correlation. The extent of this effect can be seen by rewriting (31) to explicitly account for the contributions of the useful anomalous differences to the useful anomalous correlation by including the relationship between the normalized variance (E²; equation 19) and the number of sites. Assuming independence of the terms in (18), we can write an expression for the expected mean-square measured anomalous difference 〈(Δ_ano^obs)²〉,

$[\langle(\Delta_{\rm ano}^{\rm obs})^2\rangle = \langle \Delta_{\rm ano}^2\rangle + \langle(\Delta_{\rm ano}^{\rm other})^2\rangle + \langle(\varepsilon_{\rm ano})^2\rangle. \eqno(34)]$

Then, considering (28), we can rewrite this in terms of the number of sites in the anomalous substructure n,

$[\langle(\Delta_{\rm ano}^{\rm obs})^2\rangle = 2n_{\rm sites}N_{\rm sym}\langle f_{h,B}^2\rangle +\langle(\Delta _{\rm ano}^{\rm other})^2\rangle + \langle(\varepsilon_{\rm ano})\rangle. \eqno(35)]$

Substituting into (19), we obtain an expression for the normalized variance E² that depends on the number of sites, the mean-square anomalous difference Δ_ano^other from sources other than the substructure, the errors in measurement ∊_ano and the scattering from an individual site in the substructure f_h,B,

$[E^2 \simeq {{\langle (\Delta_{\rm ano}^{\rm other})^2\rangle + \langle \varepsilon_{\rm ano})^2\rangle} \over {2n_{\rm sites}N_{\rm sym}\langle f_{h,B}^2\rangle}}. \eqno(36)]$

Factoring out the dependence of the normalized variance on the number of sites n_sites, we can write that

$[E^2 \simeq {{e^2} \over {n_{\rm sites}}}, \eqno(37)]$

where e² is the ratio of the total mean-square errors in the anomalous differences to the expected mean-square useful anomalous differences for a single site in the substructure,

$[e^2 \simeq {{\langle (\Delta_{\rm ano}^{\rm other})^2\rangle + \langle (\varepsilon_{\rm ano})^2\rangle } \over {2N_{\rm sym}\langle f_{h,B}^2\rangle}}. \eqno(38)]$

Using (37) and (38), we can rewrite (31) (or equation 28) in a way that explicitly includes the number of sites and the ratio of the total mean-square errors in the anomalous differences to the expected mean-square useful anomalous differences for a single site in the substructure,

$[\langle S_{\rm ano}^{\rm obs} \rangle \simeq \left [ {{N_{\rm refl}} \over {(n_{\rm sites} + e^2)f_B}}\right]^{1/2}. \eqno(39)]$

It can be seen from (39) that if the ratio e² is much smaller than the number of sites n_sites then increasing the number of sites will decrease the anomalous signal S_ano^obs approximately according to the square root of the number of sites. This corresponds to the situation where the data are very accurately measured and there are only small anomalous differences arising from any atoms other than those in the substructure, as in §3.1.2. In contrast, if the ratio e² is comparable to or larger than the number of sites, then changing the number of sites will have a smaller effect, but the anomalous signal will still always decrease with increasing numbers of sites in the substructure. This corresponds to the situation in which the total mean-square error in the anomalous differences is comparable to or larger than the total mean-square useful anomalous difference.

4. Conclusions

Our simple theory (31) shows how the anomalous signal in a SAD data set depends on the correlation of anomalous differences with those corresponding only to the anomalous substructure, the number of unique reflections measured, the number of sites in the substructure and the atomic displacement factors for the atoms in the substructure. Combining this with the empirical observation that the anomalous signal is a predictor of solving the anomalous substructure (Fig. 2) gives a clear idea of the features of the crystal and the experiment that determine whether the substructure can be obtained. (31) shows that even if the data are measured precisely, the anomalous signal is limited by several factors. These are the number of reflections measured, the number of sites in the substructure, the atomic displacement factors of the atoms in the substructure and the presence of minor sites. The limits on obtainable values of the anomalous signal can be used to ensure that substructure solution is at least possible when designing an experiment. (31) further shows that if errors in measurement are present, or if not all the anomalous scattering comes from the anomalous substructure, the anomalous signal is reduced by the value of the useful anomalous correlation CC_ano^obs (the correlation of anomalous differences with those expected from a crystal where anomalous differences come only from the substructure and where there are no experimental errors). We emphasize that throughout this analysis the value of each parameter in (31) is the value corresponding to all reflections in the data set, and the parameters in (31) may change as higher resolution, lower accuracy data are included. Although including high-resolution, low-accuracy anomalous differences will increase the total number of reflections, which would seem to increase the anomalous signal through (31), this effect could be offset by the resulting lower value of the overall anomalous correlation CC_ano^obs. Consequently, in order to use (31) effectively to decide which data to include in an analysis it is important to consider how each parameter in (31) would be expected to change as additional data are included.

There are two principal bottlenecks in structure determination using the SAD method. One is finding the locations of atoms in the substructure, as discussed above, and the other is the calculation of crystallographic phases (Liu et al., 2013). Fig. 4 indicates that the accuracy of experimental phases in the SAD method are closely related to the useful anomalous correlation CC_ano^obs. This is consistent with the observations of Zwart (2005) and provides a basis for predicting experimental map quality before and after the collection of SAD data.

Taken as a whole, our theoretical treatment of the contributions to the anomalous signal and useful anomalous correlation provide a foundation for evaluating whether a SAD experiment is likely to lead to successful substructure and phase determination.

Acknowledgements

The authors appreciate the support received from the US National Institutes of Health (grant P01GM063210 to PDA, TCT, Jane S. Richardson and Randy J. Read, and grant P01AI055672 to JLS). This work was partially supported by the US Department of Energy under contract DE-AC02-05CH11231.

References

Abrahams, J. P. & Leslie, A. G. W. (1996). Acta Cryst. D52, 30–42. CrossRef CAS Web of Science IUCr Journals Google Scholar
Aceti, D. J., Bitto, E., Yakunin, A. F., Proudfoot, M., Bingman, C. A., Frederick, R. O., Sreenath, H. K., Vojtik, F. C., Wrobel, R. L., Fox, B. G., Markley, J. L. & Phillips, G. N. Jr (2008). Proteins, 73, 241–253. Web of Science CrossRef PubMed CAS Google Scholar
Adams, P. D. et al. (2010). Acta Cryst. D66, 213–221. Web of Science CrossRef CAS IUCr Journals Google Scholar
Akey, D. L., Brown, W. C., Konwerski, J. R., Ogata, C. M. & Smith, J. L. (2014). Acta Cryst. D70, 2719–2729. Web of Science CrossRef IUCr Journals Google Scholar
Aylett, C. H. S. & Lowe, J. (2012). Proc. Natl Acad. Sci. USA, 109, 16522–16527. CrossRef Google Scholar
Berman, H. M., Westbrook, J., Feng, Z., Gilliland, G., Bhat, T. N., Weissig, H., Shindyalov, I. N. & Bourne, P. E. (2000). Nucleic Acids Res. 28, 235–242. Web of Science CrossRef PubMed CAS Google Scholar
Bijvoet, J. M. (1954). Nature (London), 173, 888–891. CrossRef Web of Science Google Scholar
Boudes, M., Sanchez, D., Graille, M., van Tilbeurgh, H., Durand, D. & Quevillon-Cheruel, S. (2014). Nucleic Acids Res. 42, 5302–5313. CrossRef Google Scholar
Buehner, M., Ford, G. C., Moras, D., Olsen, K. W. & Rossmann, M. G. (1974). J. Mol. Biol. 82, 563–585. CrossRef CAS PubMed Web of Science Google Scholar
Bunkóczi, G., McCoy, A. J., Echols, N., Grosse-Kunstleve, R. W., Adams, P. D., Holton, J. M., Read, R. J. & Terwilliger, T. C. (2014). Nature Methods, 12, 127–130. Web of Science PubMed Google Scholar
Bürmann, F., Shin, H., Basquin, J., Soh, Y., Giménez-Oya, V., Kim, Y., Oh, B. & Gruber, S. (2013). Nature Struct. Mol. Biol. 20, 371–379. Google Scholar
Chantler, C. T. (1995). J. Phys. Chem. Ref. Data, 24, 71. CrossRef Google Scholar
Colman, P. M., Jansonius, J. N. & Matthews, B. M. (1972). J. Mol. Biol. 70, 701–724. CrossRef Google Scholar
Cowtan, K. (2006). Acta Cryst. D62, 1002–1011. Web of Science CrossRef CAS IUCr Journals Google Scholar
Cowtan, K. (2010). Acta Cryst. D66, 470–478. Web of Science CrossRef CAS IUCr Journals Google Scholar
Cowtan, K. D. & Main, P. (1996). Acta Cryst. D52, 43–48. CrossRef CAS Web of Science IUCr Journals Google Scholar
Dauter, Z. (2006). Acta Cryst. D62, 867–876. Web of Science CrossRef CAS IUCr Journals Google Scholar
Dauter, Z., Dauter, M. & Dodson, E. J. (2002). Acta Cryst. D58, 494–506. Web of Science CrossRef CAS IUCr Journals Google Scholar
Debreczeni, J. É., Bunkóczi, G., Ma, Q., Blaser, H. & Sheldrick, G. M. (2003). Acta Cryst. D59, 688–696. Web of Science CrossRef CAS IUCr Journals Google Scholar
Du, J., Zhong, X., Bernatavichute, Y. V., Stroud, H., Feng, S., Caro, E., Vashisht, A. A., Terragni, J., Chin, H. G., Tu, A., Hetzel, J., Wohlschlegel, J. A., Pradhan, S., Patel, D. J. & Jacobsen, S. E. (2012). Cell, 151, 167–180. CrossRef Google Scholar
Emsley, P., Lohkamp, B., Scott, W. G. & Cowtan, K. (2010). Acta Cryst. D66, 486–501. Web of Science CrossRef CAS IUCr Journals Google Scholar
Evans, P. (2006). Acta Cryst. D62, 72–82. Web of Science CrossRef CAS IUCr Journals Google Scholar
Fan, J., Jiang, D., Zhao, Y., Liu, J. & Zhang, X. C. (2014). Proc. Natl Acad. Sci. USA, 111, 7636–7640. CrossRef Google Scholar
Fineran, P. C., Blower, T. R., Foulds, I. J., Humphreys, D. P., Lilley, K. S. & Salmond, G. P. C. (2009). Proc. Natl Acad. Sci. USA, 106, 894–899. CrossRef Google Scholar
Fourme, R., Shepard, W., Kahn, R., l'Hermite, G. & Li de La Sierra, I. (1995). J. Synchrotron Rad. 2, 36–48. CrossRef CAS Web of Science IUCr Journals Google Scholar
Fu, Z.-Q., Rose, J. P. & Wang, B.-C. (2004). Acta Cryst. D60, 499–506. Web of Science CrossRef CAS IUCr Journals Google Scholar
Furey, W. & Swaminathan, S. (1997). Methods Enzymol. 276, 590–620. CrossRef Web of Science Google Scholar
Gao, A. & Serganov, A. (2014). Nature Chem. Biol. 10, 787–792. CrossRef Google Scholar
Garman, E. (2003). Curr. Opin. Struct. Biol. 13, 545–551. Web of Science CrossRef PubMed CAS Google Scholar
González, A. (2003). Acta Cryst. D59, 1935–1942. Web of Science CrossRef IUCr Journals Google Scholar
Grosse-Kunstleve, R. W. & Adams, P. D. (2003). Acta Cryst. D59, 1966–1973. Web of Science CrossRef CAS IUCr Journals Google Scholar
He, X., Kuo, Y.-C., Rosche, T. J. & Zhang, X. (2013). Structure, 21, 355–364. CrossRef Google Scholar
Hendrickson, W. A. (2014). Q. Rev. Biophys. 47, 49–93. Web of Science CrossRef PubMed Google Scholar
Hendrickson, W. A. & Teeter, M. M. (1981). Nature (London), 290, 107–113. CrossRef CAS Web of Science Google Scholar
Hou, X., Pedi, L., Diver, M. M. & Long, S. (2012). Science, 338, 1308–1313. CrossRef Google Scholar
Howell, P. L. & Smith, G. D. (1992). J. Appl. Cryst. 25, 81–86. CrossRef Web of Science IUCr Journals Google Scholar
Hsia, K.-C. & Hoelz, A. (2010). Proc. Natl Acad. Sci. USA, 107, 11271–11276. Web of Science CrossRef CAS PubMed Google Scholar
Huang, H., Levin, E. J., Liu, S., Bai, Y., Lockless, S. W. & Zhou, M. (2014). PLoS Biol. 12, e1001911. CrossRef Google Scholar
Hwang, W. C., Golden, J. W., Pascual, J., Xu, D., Cheltsov, A. & Godzik, A. (2014). Proteins, doi:10.1002/prot.24679. Google Scholar
Jones, T. A., Zou, J.-Y., Cowan, S. W. & Kjeldgaard, M. (1991). Acta Cryst. A47, 110–119. CrossRef CAS Web of Science IUCr Journals Google Scholar
Kartha, G. & Parthasarathy, R. (1965). Acta Cryst. 18, 745–749. CrossRef CAS IUCr Journals Web of Science Google Scholar
Kozbial, P. et al. (2008). Proteins, 71, 1589–1596. CrossRef Google Scholar
Krojer, T., Pike, A. C. W. & von Delft, F. (2013). Acta Cryst. D69, 1303–1313. CrossRef IUCr Journals Google Scholar
La Fortelle, E. de & Bricogne, G. (1997). Methods Enzymol. 276, 472–494. Google Scholar
Langer, G., Cohen, S. X., Lamzin, V. S. & Perrakis, A. (2008). Nature Protoc. 3, 1171–1179. Web of Science CrossRef CAS Google Scholar
Leibundgut, M., Frick, C., Thanbichler, M., Böck, A. & Ban, N. (2004). EMBO J. 24, 11–22. CrossRef Google Scholar
Levin, I. et al. (2005). Proteins, 59, 864–868. CrossRef Google Scholar
Lindås, A.-C., Chruszcz, M., Bernander, R. & Valegård, K. (2014). Acta Cryst. D70, 492–500. CrossRef IUCr Journals Google Scholar
Liu, Q., Liu, Q. & Hendrickson, W. A. (2013). Acta Cryst. D69, 1314–1332. Web of Science CrossRef CAS IUCr Journals Google Scholar
Liu, Z.-J., Chen, L., Wu, D., Ding, W., Zhang, H., Zhou, W., Fu, Z.-Q. & Wang, B.-C. (2011). Acta Cryst. A67, 544–549. Web of Science CrossRef CAS IUCr Journals Google Scholar
Lo, Y.-C., Lin, S.-C., Rospigliosi, C. C., Conze, D. B., Wu, C.-J., Ashwell, J. D., Eliezer, D. & Wu, H. (2009). Mol. Cell, 33, 602–615. CrossRef Google Scholar
McCoy, A. J., Storoni, L. C. & Read, R. J. (2004). Acta Cryst. D60, 1220–1228. Web of Science CrossRef CAS IUCr Journals Google Scholar
Mechaly, A. E., Sassoon, N., Betton, J.-M. & Alzari, P. M. (2014). PLoS Biol. 12, e1001776. CrossRef Google Scholar
Mueller, M., Grauschopf, U., Maier, T., Glockshuber, R. & Ban, N. (2009). Nature (London), 459, 726–730. Web of Science CrossRef PubMed CAS Google Scholar
North, A. C. T. (1965). Acta Cryst. 18, 212–216. CrossRef IUCr Journals Web of Science Google Scholar
Osawa, T., Ito, K., Inanaga, H., Nureki, O., Tomita, K. & Numata, T. (2009). Structure, 17, 713–724. CrossRef Google Scholar
Otwinowski, Z. (1991). Proceedings of the CCP4 Study Weekend. Isomorphous Replacement and Anomalous Scattering, edited by W. Wolf, P. R. Evans & A. G. W. Leslie, pp. 80–86. Warrington: Daresbury Laboratory. Google Scholar
Otwinowski, Z. & Minor, W. (1997). Methods Enzymol. 276, 307–326. CrossRef CAS Web of Science Google Scholar
Pannu, N. S. & Read, R. J. (2004). Acta Cryst. D60, 22–27. Web of Science CrossRef CAS IUCr Journals Google Scholar
Parthasarathy, S. & Parthasarathi, V. (1974). Acta Cryst. A30, 649–654. CrossRef IUCr Journals Web of Science Google Scholar
Perrakis, A., Morris, R. & Lamzin, V. S. (1999). Nature Struct. Biol. 6, 458–463. Web of Science CrossRef PubMed CAS Google Scholar
Pokkuluri, P. R., Londer, Y. Y., Duke, N. E., Pessanha, M., Yang, X., Orshonsky, V., Orshonsky, L., Erickson, J., Zagyanskiy, Y., Salgueiro, C. A. & Schiffer, M. (2011). J. Struct. Biol. 174, 223–233. CrossRef Google Scholar
Qiao, R., Cabral, G., Lettman, M. M., Dammermann, A. & Dong, G. (2012). EMBO J. 31, 4334–4347. CrossRef Google Scholar
Sampathkumar, P. et al. (2013). Structure, 21, 560–571. CrossRef Google Scholar
Schäfer, I. B., Rode, M., Bonneau, F., Schüssler, S. & Conti, E. (2014). Nature Struct. Mol. Biol. 21, 591–598. Google Scholar
Schmitzberger, F. & Harrison, S. C. (2012). EMBO Rep. 13, 216–222. CrossRef Google Scholar
Schneider, T. R. & Sheldrick, G. M. (2002). Acta Cryst. D58, 1772–1779. Web of Science CrossRef CAS IUCr Journals Google Scholar
Shen, Q., Wang, J. & Ealick, S. E. (2003). Acta Cryst. A59, 371–373. Web of Science CrossRef CAS IUCr Journals Google Scholar
Strahs, G. & Kraut, J. (1968). J. Mol. Biol. 35, 503–512. CrossRef CAS PubMed Web of Science Google Scholar
Terwilliger, T. C. (2000). Acta Cryst. D56, 965–972. Web of Science CrossRef CAS IUCr Journals Google Scholar
Terwilliger, T. C. & Berendzen, J. (1999). Acta Cryst. D55, 849–861. Web of Science CrossRef CAS IUCr Journals Google Scholar
Terwilliger, T. C., Bunkóczi, G., Hung, L.-W., Zwart, P. H., Smith, J., Akey, D. & Adams, P. D. (2016). Acta Cryst. D72, 359–374. CrossRef IUCr Journals Google Scholar
Terwilliger, T. C., Grosse-Kunstleve, R. W., Afonine, P. V., Moriarty, N. W., Zwart, P. H., Hung, L.-W., Read, R. J. & Adams, P. D. (2008). Acta Cryst. D64, 61–69. Web of Science CrossRef CAS IUCr Journals Google Scholar
Tominaga, T., Watanabe, S., Matsumi, R., Atomi, H., Imanaka, T. & Miki, K. (2012). Acta Cryst. F68, 1153–1157. CrossRef IUCr Journals Google Scholar
Wang, B.-C. (1985). Methods Enzymol. 115, 90–112. CrossRef CAS PubMed Google Scholar
Wang, L., Yang, X., Li, S., Wang, Z., Liu, Y., Feng, J., Zhu, Y. & Shen, Y. (2014). EMBO J. 33, 594–604. CrossRef Google Scholar
Weeks, C. M., DeTitta, G. T., Miller, R. & Hauptman, H. A. (1993). Acta Cryst. D49, 179–181. CrossRef CAS Web of Science IUCr Journals Google Scholar
Weinert, T. et al. (2015). Nature Methods, 12, 131–133. Web of Science CrossRef CAS PubMed Google Scholar
Xu, Q. et al. (2006). Proteins, 62, 292–296. CrossRef Google Scholar
Xu, Q. et al. (2007). Proteins, 69, 433–439. CrossRef Google Scholar
Xu, Q., Rife, C. L. et al. (2009). Proteins, 74, 1041–1049. CrossRef Google Scholar
Xu, Q., Sudek, S. et al. (2009). Structure, 17, 303–313. Web of Science CrossRef PubMed CAS Google Scholar
Yang, C., Pflugrath, J. W., Courville, D. A., Stence, C. N. & Ferrara, J. D. (2003). Acta Cryst. D59, 1943–1957. Web of Science CrossRef CAS IUCr Journals Google Scholar
Zhang, J. & Ferré-D'Amaré, A. R. (2013). Nature (London), 500, 363–366. CrossRef Google Scholar
Zhou, Q., Li, J., Yu, H., Zhai, Y., Gao, Z., Liu, Y., Pang, X., Zhang, L., Schulten, K., Sun, F. & Chen, C. (2014). Nature Commun. 5, 3552. Google Scholar
Zhou, Y. F. & Springer, T. A. (2014). Blood, 123, 1785–1793. CrossRef Google Scholar
Zimmer, J., Nam, Y. & Rapoport, T. A. (2008). Nature (London), 455, 936–943. CrossRef Google Scholar
Zwart, P. H. (2005). Acta Cryst. D61, 1437–1448. Web of Science CrossRef CAS IUCr Journals Google Scholar

This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.

STRUCTURAL
BIOLOGY

ISSN: 2059-7983

Volume 72| Part 3| March 2016| Pages 346-358

doi:10.1107/S2059798315019269

Open

access

Format		BIBTeX
		EndNote
		RefMan
		Refer
		Medline
		CIF
		SGML
		Plain Text
		Text

Format		BIBTeX
		EndNote
		RefMan
		Refer
		Medline
		CIF
		SGML
		Plain Text
		Text

Search IUCr Journals		doi		Advanced search
Author		volume	page

research papers\(\def\hfill{\hskip 5em}\def\hfil{\hskip 3em}\def\eqno#1{\hfil {#1}}\)

Can I solve my structure by SAD phasing? Anomalous signal in SAD phasing

1. Introduction

1.1. Single-wavelength anomalous diffraction

1.2. Measures of signal and noise in anomalous data

1.2.1. Metrics predictable from experimental setup

1.2.2. Metrics calculated from measured anomalous differences

1.2.3. Metrics requiring known structure

1.3. Objectives

2. Methods

2.1. Structure-factor relationships and anomalous differences in anomalous scattering

2.2. Contributions to the anomalous signal

2.2.1. Anomalous signal for an ideal anomalous difference Fourier

2.2.2. Anomalous signal for a realistic anomalous difference Fourier

2.3. Test data from the PDB

2.4. Software availability

3. Results and discussion

3.1. Dependence of the anomalous signal Sanoobs on CCano, Nrefl, nsites and fB

3.1.1. Anomalous signal for an idealized crystal

3.1.2. Dependence of the anomalous signal for an idealized crystal on resolution

3.1.3. Anomalous signal for anomalous data measured with errors from a real crystal

3.1.4. Anomalous signal and resolution of the data included in the calculation

3.1.5. Comparison of expected and actual anomalous signal

3.2. Relationship between the anomalous signal and the solution of the anomalous substructure

3.3. Relationship between useful anomalous correlation and the quality of phase calculation

3.4. Relationship between number of sites in the anomalous substructure and the anomalous signal

4. Conclusions

Acknowledgements

References

research papers

3.1. Dependence of the anomalous signal S_ano^obs on CC_ano, N_refl, n_sites and f_B