feature articles
What happens when the signs of anomalous differences or the handedness of
are inverted?aSAIC-Frederick Inc., Basic Research Program, Argonne National Laboratory, Argonne, IL 60439, USA, bProtein Structure Section, MCL, National Cancer Institute, NCI-Frederick, Frederick, MD 21702, USA, and cSynchrotron Radiation Research Section, MCL, National Cancer Institute, Argonne National Laboratory, Argonne, IL 60439, USA
*Correspondence e-mail: dauter@anl.gov
Proper solution of a macromolecular
based on and/or isomorphous differences requires that the anomalous differences in reflection amplitudes be measured properly and that the correct of the be selected. If this information is wrong then the resulting electron-density maps will not show the correct structural features, but the reflection phases and map features will be related to the correct ones in a specific way. This text aims to explain how misinterpretation of the Bijvoet differences or of the affects the resulting phases and electron-density maps.Keywords: phasing; anomalous scattering; SAD; indexing errors.
1. Introduction
Recently, Matthews (2007) published an excellent and thorough discussion of the principles of phasing macromolecular crystal structures based on the anomalous diffraction signal, showing that change of the signs of anomalous differences or misinterpretation of the handedness of the constellation of anomalous scatterers will never lead to correctly interpretable electron-density maps. His article is a reaction to the recent retraction of erroneous structures resulting from the wrong interpretation of experimental diffraction data, leading to inversion of the signs of anomalous differences (Chang et al., 2006; Miller, 2007).
As Matthews correctly stated, phasing methods based on the anomalous diffraction signal will only lead to the correct structure if the reflections are indexed properly (in the right-handed coordinate system) and if the handedness of the constellation of anomalous scatterers, one of two possibilities, is correctly chosen. A special case occurs if the
has a centrosymmetric arrangement, where the problem of its handedness vanishes. However, the phases and electron-density maps obtained with wrongly assigned anomalous differences or handedness of the are related to the correct solution in a specific way. These relationships will be elaborated with emphasis on the single-wavelength anomalous diffraction (SAD) method of phasing.2. SAD phasing
2.1. Correct interpretation
Fig. 1 illustrates the relationships between various contributions to a if the crystal contains atoms displaying an effect (for simplicity, it is assumed that all anomalous scatterers are of the same kind). The contribution of the normally scattering part is represented by the brown vector FN and the normal scattering of the anomalous atoms by the red vector FA. The real part of the is represented by the vector FA′ and the imaginary part by the purple vectors FA′′. They result from the anomalous corrections to the atomic form factor: fA = fAo + fA′ + ifA′′. The vector FA′ is collinear and antiparallel to FA, since fA′ is real and usually negative. The purple vector +FA′′ is perpendicular and rotated anticlockwise with respect to FA, since the correction fA′′ is multiplied by the imaginary factor i = (−1)1/2. The purple vector −FA′′ represents the anomalous contribution of the negative Friedel mate or, more strictly, it represents the complex conjugate of the Friedel-related reflection, which in fact has an opposite phase and lies on the other side of the horizontal axis of the Argand diagram. If it were drawn as the original reflection, not its mirror image, the vector FA′′(−h) would also be rotated anticlockwise with respect to FA(−h). The black vectors FT+ and FT− therefore represent the total of a pair of Friedel-related reflections. Only the lengths of these two vectors (reflection amplitudes) are available in a single-wavelength diffraction data set. The anomalous difference ΔF± is equal to the difference between the two amplitudes within a Friedel (or Bijvoet) pair. A more detailed discussion of the relevant vector and phase relations can be found elsewhere (Blundell & Johnson, 1976; Drenth, 1999; Dauter, 2002).
The measurable anomalous difference ΔF± is approximately related to the total contribution 2FA′′ through the sinusoidal function, depending on the difference in the directions of the total average FT and of the factor representing the anomalous (blue and red vectors in Fig. 1): ΔF± = |F+| − |F−| ≃ 2FA′′sin(φT − φA) = 2δAFAsin(φT − φA), where δA = FA′′/FA is characteristic for a particular scatterer. In principle, to locate the positions of anomalous scatterers one should use the amplitudes |FA| representing the total contribution of the anomalous These amplitudes may be estimated in the multi-wavelength anomalous diffraction (MAD) experiment, although not from a single SAD data set. However, because the measurable anomalous differences are related to the FA values, it is possible to use them for the location of anomalous sites using Patterson or Since sin2(α) = [1 − cos(2α)]/2, the Patterson synthesis calculated with |ΔF±|2 ≃ 4δA2FA2sin2(φT − φA) = 2δA2FA2 − 2δA2FA2cos[2(φT − φA)] will reproduce the proper vectors representing the anomalous scatterers (first term) accompanied by (more-or-less random) noise resulting from the second cosine term. Location of the anomalous scatterer sites is the first step in the SAD phasing protocol.
If the anomalous FA′′ and −FA′′ must lie on the two circles with radii corresponding to |FT+| and |FT−|. As illustrated in Fig. 2(a), since sin(φT − φA) = sin(π − φT + φA), there are two solutions to this geometrical problem, giving rise to the well known phase ambiguity inherent to SAD (as well as SIR) phasing (Drenth, 1999).
is located, it is possible to calculate its contribution; therefore, the length and orientation of the red and purple vectors are known. In addition to this, the lengths, but not the directions, of the total amplitudes of both Friedel mates are available. It is therefore required that the ends of the (correctly oriented) purple vectors +The two solutions in the SAD approach are not absolutely symmetrical, as is evidenced by the different lengths of the brown vectors FN corresponding to the normal scatterers. In effect, the anomalous and normal substructures represent different fractions of the total crystal in the two possible phase solutions, which increases the probability that the solution with smaller FN and a smaller angle between FA and FT is correct (Ramachandran & Raman, 1956). If the anomalous is small, this influence of the known (anomalous) described by Sim (1959, 1964), is very small. In the simplest approach, one can accept the average phase φSAD = φA − π/2 for calculation of the initial electron-density map. As pointed out by Wang (1985), such a map will correspond to a superposition of the correct image of the structure (contributed by the set of correctly chosen phases) and of the featureless noisy map (resulting from the contribution of the wrong phases). This map, although not necessarily directly interpretable, may at least indicate which parts of the cell correspond to the macromolecule and which to the much flatter solvent regions. The subsequent iterative density-modification procedures (such as solvent flattening) should enhance the correct features and indicate the correct choice of the alternative SAD phase. In conclusion, the properly indexed data set and correctly chosen handedness of the lead to the SAD map showing the positively contoured features of the macromolecule of correct with L-amino acids in proteins or D-sugars in nucleic acids.
2.2. of opposite handedness
In the SAD approach the anomalous scatterer sites can be located either by analysis of the Patterson syntheses or by the use of
with both methods utilizing anomalous differences. Neither of these methods differentiates between the two possible enantiomorphic constellations of atoms in the Both outcomes are equally probable and both have to be checked for correctness. If the Patterson or direct-methods solution correctly represents the anomalous in the crystal, the initial SAD map should disclose correct (although not very accurate) features of the However, the opposite handedness of the results in a different set of phases and a map with different features.This phenomenon is illustrated in Fig. 2(b). All vectors representing various scattering contributors have the same length as previously, but the phase φA of the anomalous is negated, corresponding to the opposite As a consequence, the two alternative phase solutions are different and the resulting vector FcSAD lies on the other side of the vertical axis of the Argand diagram, with the corresponding phase φcSAD = π − φoSAD. Since cos(π − φ) = −cos(φ), this has the following consequence for the calculated Fourier map:
The Fourier synthesis will therefore be related to the correct map by centrosymmetric inversion and will have negative values. It will show the structural features of inverted
in the negative contours of the map.Since the positive contours of such a map do not show any significant features except noise and the solvent-flattening procedure removes all negative features, such a map will not lead to solution of the crystal structure.
2.3. Data with inverted signs of anomalous differences
If the indices of Friedel-related reflections are inadvertently exchanged before SAD phasing, the resulting SAD phase φiSAD will be increased by π in comparison with the original, correct phase φoSAD, as illustrated in Fig. 2(c). Since cos(π + φ) = −cos(φ), the map calculated with such phases will have the following relation to the original map:
It will therefore represent the correct structure, but in negative contours. Again, the solvent flattening will fail to show any interpretable structural features.
If, in addition, the wrong handedness of the anomalous φicSAD will have a negative value, as shown in Fig. 2(d). Taking into account that cos(−φ) = cos(φ), the resulting map will have the following property:
is chosen, the SAD phaseThis map will show the positive features, but inverted by the center of symmetry. The solvent-flattening procedure will enhance these features, but all chiral centers will have wrong handedness, helices in proteins will be left-handed etc.
2.4. Centrosymmetric constellation of anomalous sites
A special case may occur if all anomalous scatterers are centrosymmetically arranged in the crystal P1 cell or one unique site in the P21 cell are examples of this case. Another well known example is provided by the rhombohedral R3 insulin crystal with two Zn atoms, both of which lie at the threefold axis (Dauter et al., 2002). The relevant phase diagram is shown in Fig. 2(e). As a consequence of the centrosymmetric arrangement of the anomalous sites, the red vector FA must be horizontal, with its phase φA equal either to 0 or π (if the center lies at the cell origin). The average phase φSAD, which is always orthogonal to φA, must therefore be equal to ±π/2. Since the in this case is centrosymmetric, its inversion does not change anything; indeed, if φSAD = ±π/2 then π − φSAD = ±π/2 as well and both conditions described above in §2.1 and §2.2 are fulfilled; therefore, ρo(r) = −ρo(−r). If the anomalous signal is measured correctly, the Fourier map will be anti-centrosymmetric and contain both the correct positive structural features and the inverted features in its negative contours. Such a map is therefore amenable to solvent flattening in the usual way. However, if the signs of the anomalous differences are exchanged [both cases in §2.3 fulfilled, with ρo(−r) = −ρo(r)], the Fourier map will also be anti-centrosymmetric and contain correct features in the negative contours, but the positive features will have inverted chirality.
One or two anomalous scatterer sites in the3. MAD phasing
Conceptually, multiwavelength anomalous diffraction phasing can be analyzed as the result of several SAD cases with different amounts of anomalous signal and this is illustrated in Fig. 3. Each wavelength indicates a pair of possible phase solutions (as in the SAD case), but only one solution in each pair coincides for all wavelengths. If more than one data set is available, there is no need to use the average `SAD' phase; it is immediately possible to select proper values of the initial phases. In principle, the MAD map should show the correct structure with a low level of noise even before it is density-modified.
The behavior of all individual phase estimations for all wavelengths after inversion of the anomalous signal or of the handedness of a 2 for SAD phasing, with the same consequences for the map features.
is therefore exactly the same as discussed in §4. phasing
The
method utilizes differences in the amplitudes measured from a native crystal and from one or more derivatives containing heavy atoms attached to the crystallized macromolecule. If only one derivative is used the technique is called single (SIR); in the case of more than one derivative it is called multiple (MIR). These phasing approaches do not make use of anomalous diffraction effects, even if the heavy atoms can potentially display properties.4.1. SIR phasing
The principle of SIR phasing is illustrated in Fig. 4(a). From the diffraction experiment, only the amplitudes of the native data and of one derivative data set are available. For a particular reflection it is initially only known that the end of the brown vector (representing the native protein amplitude FP) must lie on the Argand diagram at the circle of radius |FP| and, analogously, that the blue vector corresponding to the derivative amplitude FPH must end at the circle of radius |FPH|. The location of the heavy atoms can be obtained from Patterson or direct-methods calculations based on the observed isomorphous differences (|FPH| − |FP|), since they are related to the heavy-atom amplitudes FH through a trigonometric relation depending on the mutual disposition of FP and FH, i.e. on the difference of the phases (φH − φP). This is analogous to the relationship between ΔF± and FA in the SAD case. If the heavy-atom positions are known, it is possible to calculate their contribution FH and locate on the Argand diagram two possible solutions of the protein phase. Since it is a priori not known which of them is correct, the `average' phase φSIR can be used for map calculation, which should show the correct image of the structure blurred by a degree of noise. A density-modification procedure may lead to an interpretable electron-density map and therefore to solution of the structure. This procedure is closely analogous to SAD phasing, except that in SIR the two alternative phase solutions are fully symmetric with respect to the lengths of all contributing vectors and the Sim probability contribution is equal for both possibilities. This results from the fact that in SIR the average phase φSIR is the same as the phase φA, whereas in SAD φSAD is shifted by π/2 from φA.
Since in the SIR approach the φH, is negated, which is shown in Fig. 4(b). The resulting phases will also have opposite signs and the corresponding Fourier map will show the image of the structure inverted by the center of symmetry.
effect is not used at all, it does not depend on the proper assignment of the signs of Bijvoet differences. However, the location of heavy atoms by Patterson or may equally probably lead to the correct or centrosymmetrically inverted constellation. If the of heavy atoms has inverted handedness, the phase of its contribution,If the heavy atoms are arranged in the π. In this case the Fourier map will also be centrosymmetric and will contain two superimposed images of the structure. In this case density-modification procedures will not lead to the structure solution, in contrast to the situation with SAD/MAD, where a centrosymmetric results in the anti-centrosymmetric Fourier map.
centrosymmetrically, the resulting SIR phases will have values of 0 or4.2. MIR phasing
Similarly to the second-wavelength data in MAD, the second heavy-atom derivative data in MIR indicate which of the alternative SIR phases is correct, as illustrated in Fig. 4(c). If there are more derivatives, such an indication becomes stronger. It is assumed that the contributions of the various heavy atoms are not equivalent, since obviously for reliable phase indication the heavy-atom vectors (purple and red in Fig. 4c) must differ. Moreover, the contribution of the macromolecule itself (represented by the brown vectors) must be the same in all the derivatives, which means that all the crystals should be isomorphous, hence the name of this method. If the introduction of heavy atoms affects the orientation of the protein molecule in the or causes changes in the crystal unit-cell parameters, the isomorphicity is no longer preserved and all MIR/SIR phase indications will be aberrated by errors. This effect is not so detrimental for the MAD/SAD techniques, where all data are usually collected from a single crystal. Obviously, all these phasing methods are negatively affected by inaccuracies in the measured diffraction intensities and by the radiation damage inflicted on crystals by exposure to X-rays.
The effect of wrong handedness of the heavy-atom
is the same as in the SIR case and results in positive but centrosymetrically inverted map features.5. Combination of and anomalous diffraction signals
Fig. 5 illustrates the simplest case of SIRAS, with one derivative displaying a significant effect. Individually, each SAD and SIR contribution leads to two possible phase solutions, but only one of them is common to both approaches, so that the phase indication in this case is in principle unequivocal. In general, in the MIRAS approach the consequences of incorrect assignment of anomalous differences or may be more complicated than explained above. A simple summary is presented in Table 1.
|
Obviously, with correct signs of Bijvoet differences and proper handedness of the
the Fourier map will represent the true However, another clear case results if the anomalous signs as well as the are inverted, when the phase indications contributed by both isomorphous and anomalous differences are wrong but in agreement. The resulting Fourier map will show the image of the structure inverted by the center of symmetry.In the two remaining combinations, the contributions of isomorphous and anomalous effects contradict each other. The resulting phases will depend on the relative strength of the two effects. If one of the sources of phasing significantly prevails, the Fourier map may resemble to some extent the correct structural features in its positive or negative contours. However, it is certainly best to make sure that the interpretation of the anomalous differences and of the
is correct.6. Conclusions
Inversion of handedness of the anomalous
and change of sign of the anomalous differences result in a change of the phases of all reflections. The corresponding features of the Fourier maps calculated with such phases depend on the property of the trigonometric cosine function appearing in the Fourier synthesis in the following way:As clearly stated by Matthews (2007), in each of these cases it should be easily possible to identify the type of error at the early stage of map interpretation.
In the special case of a centrosymmetric constellation of the anomalous
the problem of its handedness obviously does not exist and the phasing procedures should lead to the correct However, if the signs of the anomalous differences are changed, the map will show features with wrong handedness.The latter situation caused the unfortunate need to retract the wrongly interpreted structures of ABC and EmrE transporters (Chang et al., 2006). Two Os atoms in the P1 cell form a centrosymmetric constellation, assuming comparable occupancies. The inadvertent exchange of the signs of all anomalous differences at the data-processing stage must have led to a nonrandom map that, however, exhibited inverted Thus, the only explanation of the origin of the incorrect model would be that the right-handed helices were actually built into the left-handed map features, something that would be possible in a map at low resolution.
All detectors and data-processing programs in current routine use produce data that are indexed in the correct right-handed coordinate system. Inversion of the anomalous signal may occur only as a result of using programs or procedures that have not been properly validated.
Acknowledgements
This work was supported in part by the Intramural Research Program of the NIH, National Cancer Institute, Center for Cancer Research and with Federal funds from the National Cancer Institute, National Institutes of Health under Contract No. NO1-CO-12400. The content of this publication does not necessarily reflect the views or policies of the Department of Health and Human Services, nor does the mention of trade names, commercial products or organizations imply endorsement by the US Government.
References
Blundell, T. L. & Johnson, L. N. (1976). Protein Crystallography. New York: Academic Press. Google Scholar
Chang, G., Roth, C. B., Reyes, C. L., Pornillos, O., Chen, Y.-J. & Chen, A. P. (2006). Science, 314, 1875. CrossRef PubMed Google Scholar
Dauter, Z. (2002). Acta Cryst. D58, 1958–1967. 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
Drenth, J. (1999). Principles of Protein X-ray Crystallography, 2nd ed. Heidelberg: Springer. Google Scholar
Matthews, B. W. (2007). Protein Sci. 16, 1013–1016. Web of Science CrossRef PubMed CAS Google Scholar
Miller, G. (2007). Science, 314, 1856–1857. CrossRef Google Scholar
Ramachandran, G. N. & Raman, S. (1956). Curr. Sci. 25, 348–351. CAS Google Scholar
Sim, G. A. (1959). Acta Cryst. 12, 813–815. CrossRef IUCr Journals Web of Science Google Scholar
Sim, G. A. (1964). Acta Cryst. 17, 1072–1073. CrossRef IUCr Journals Web of Science Google Scholar
Wang, B.-C. (1985). Methods Enzymol. 115, 90–112. CrossRef CAS PubMed Google Scholar
© International Union of Crystallography. Prior permission is not required to reproduce short quotations, tables and figures from this article, provided the original authors and source are cited. For more information, click here.