Home | About | Journals | Submit | Contact Us | Français |

**|**IUCr Open Access Articles**|**PMC2852309

Formats

Article sections

- Abstract
- Introduction
- Basics of experimental phasing methods
- Phasing from unmerged data
- ‘Broken symmetries’ in macromolecular crystallography
- Implementation
- Discussion and conclusion
- References

Authors

Related links

Acta Crystallogr D Biol Crystallogr. 2010 April 1; 66(Pt 4): 447–457.

Published online 2010 March 24. doi: 10.1107/S0907444909053578

PMCID: PMC2852309

Correspondence e-mail: marc.schiltz/at/epfl.ch

Experimental phasing and radiation damage

Received 2009 April 14; Accepted 2009 December 11.

Copyright © Schiltz & Bricogne 2010

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

This article has been cited by other articles in PMC.

The space-group symmetry of a crystal structure imposes a point-group symmetry on its diffraction pattern, giving rise to so-called symmetry-equivalent reflections. Instances in macromolecular crystallography are discussed in which the symmetry in reciprocal space is broken, *i.e.* where symmetry-related reflections are no longer equivalent. Such a situation occurs when the sample suffers from site-specific radiation damage during the X-ray measurements. Another example of broken symmetry arises from the polarization anisotropy of anomalous scattering. In these cases, the genuine intensity differences between symmetry-related reflections can be exploited to yield phase information in the structure-solution process. In this approach, the usual separation of the data merging and phasing steps is abandoned. The data are kept unmerged down to the Harker construction, where the symmetry-breaking effects are explicitly modelled and refined and become a source of supplementary phase information.

The space-group symmetry of a crystal structure generally imposes a point-group symmetry in reciprocal space, giving rise to so-called symmetry-equivalent reflections. This usually introduces a certain redundancy in diffraction data recorded with an area detector and is exploited in X-ray crystallography to increase the accuracy of the data by averaging over the symmetry-equivalent measurements (data merging). However, genuine intensity differences between symmetry-related reflections can occur when the sample suffers from site-specific radiation-induced changes during X-ray data collection (Schiltz *et al.*, 2004 ; Schiltz & Bricogne, 2007 ) or in the presence of the rather obscure phenomenon of anisotropy of anomalous scattering (AAS; Bricogne *et al.*, 2005 ; Sanishvili *et al.*, 2007 ; Schiltz & Bricogne, 2008 ). In such cases, the use of unmerged data can generate phasing power through the intensity differences of symmetry-related reflections.

In the present communication, we will review the symmetry-breaking effects in reciprocal space which arise from site-specific radiation damage and AAS. We start with a brief introduction to experimental phasing methods and show how the effects of site-specific radiation damage and AAS can be incorporated into the general scheme of experimental phasing methods using an extended Harker construction. In line with the general philosophy of contributions to the CCP4 Study Weekend, this is intended to be a ‘tutorial’ paper, in which formal developments are kept to a minimum. The presentation given here is intended to serve as an introduction and guide to the more comprehensive papers that we have published on this subject.

Experimental phasing methods are, in their bare essence, based on interference effects between the wave scattered from the unknown structure and the wave scattered from a set of atoms whose positions are, at least approximately, known. This is called the *substructure* (Fig. 1 ). The X-rays scattered by the substructure provide a reference wave of known phase and amplitude which interferes with the wave scattered from the unknown structure. Depending on the relative phase difference between these two waves, the superposition will yield partially constructive or destructive interference effects which produce measurable intensity modulations. In this way, information about the unknown (and not directly observable) phases is transferred to measurable intensity variations.

Experimental phasing viewed as an interference experiment. The X-rays scattered by a subset of atoms (the substructure) provide a reference wave of known phase and amplitude which interferes with the wave scattered from the unknown structure (in red). **...**

Formally, for a given reflection **h**, let us denote the structure amplitude of the wave scattered from the unknown part by **P**(**h**) = *P*(**h**)exp[*i*ϕ_{P}(**h**)] and the structure amplitude of the reference wave, scattered from the substructure, by **H**(**h**) = *H*(**h**)exp[*i*ϕ_{H}(**h**)]. Their superposition then yields for the structure amplitude of the diffracted wave

Its square modulus, which is proportional to the diffracted intensity, yields an expression that contains an interference term which depends on the relative phase difference ϕ_{P}(**h**) − ϕ_{H}(**h**),

Even if the complex amplitude of the reference wave, **H**(**h**), is known, a single measurement *F*
^{2}(**h**) will not yield all the information about the unknown structure since there are two unknowns [*P*(**h**) and ϕ* _{P}*(

each with a different reference wave **H**
_{j}(**h**). Putting it differently, the complex structure factor, **F**(**h**), is split up into an unknown part, **P**(**h**), which is constant, and a known part, **H**(**h**), which is varied between different measurements. Each equation (3) defines a circle on the complex plane with radius *F*
_{j}(**h**) and centre −**H**
_{j}(**h**). The set of all such *N*
_{h} circles intersect at the point **P**(**h**) and represents what is known as the *Harker construction* (Harker, 1956 ; Fig. 2 ).

In the method of isomorphous replacement (Green *et al.*, 1954 ; Harker, 1956 ), the variations in the reference wave **H**(**h**) are achieved by collecting diffraction data from several crystals (different *derivatives*) containing different substructures (usually consisting of heavier atoms). In modern treatments of this method a native crystal has no special status and is simply considered as a derivative for which **H**(**h**) = 0 (de La Fortelle & Bricogne, 1997 ). The standard radiation-induced phasing (RIP) method (Ravelli *et al.*, 2003 ; Nanao & Ravelli, 2006 ), in which a diffraction data set is first collected from a fresh crystal followed by exposure to X-ray or UV radiation and subsequent collection of a second data set containing site-specific radiation damage, can be subsumed under the category of isomorphous replacement.

In these cases, the index *j* labels the various derivatives. For success of the method, it is important that the constant part, whose structure factor is **P**(**h**), remains essentially unaltered in the different derivatives: a condition known as isomorphism. In other words, the variations in the diffracted intensities must solely arise from modulations of the substructure.

In the multi-wavelength anomalous diffraction (MAD) method (Hendrickson, 1991 ), the reference wave is modulated by exploiting the wavelength-dependence of the atomic scattering factors *f*′ and *f*′′ in the vicinity of an absorption edge of those atoms that constitute the substructure (Fig. 3 ). In this case the index *j* labels the various data sets recorded at different wavelengths. Again, it is important that the constant part **P**(**h**) remains essentially unaffected by the wavelength changes.

In the above examples the label *j* refers to a *data set*, *e.g.* a set of reflections recorded on a particular isomorphous derivative or a set of reflections recorded at a particular wavelength. It is often assumed that the set of all structure-factor amplitudes *F*
_{j}(**h**) which share the same label *j* would need to form a coherent data set related to a particular crystal structure. However, this is not a necessary requirement since the Harker construction (3) is set up for each individual reflection **h**. The index *j* is really just a generic label which encodes book-keeping information about the experimental parameters that selectively modulate the reference wave **H**
* _{j}*(

In many circumstances, variations in the intensities of symmetry-related reflections or of repeated measurements of the same reflection simply arise from geometric or experimental factors such as differences in absorption, irradiated crystal volume, incident beam flux or decay owing to overall (nonspecific) radiation damage. These variations do not specifically affect the substructure and therefore do not produce differential modulations which are exploitable for phase determination *via* the Harker construction. They can be corrected for by multiplicative (scale) factors which are usually determined empirically and applied to the measured intensities. Since such intensity variations cannot generate phase information, the various symmetry-related intensity measurements and repeated measurements of the same reflection are usually merged into a single structure-factor amplitude after the correction factors have been applied.

However, there are several instances in macromolecular crystallography where symmetry-related reflections and/or repeated measurements of the same reflection display specific variations in the substructure amplitude and therefore give rise to intensity differences for which adequate correction cannot be made by multiplicative (scale) factors. Provided that such symmetry-breaking effects in the substructure can be modelled and refined by a set of parameters in real space (*e.g.* coordinates of atomic positions, occupancy factors, atomic scattering factors *etc*.) they can become a source of phase information (Fig. 4 ). This requires a paradigm shift in the data-processing strategy, since the usual separation of the data-merging and phasing steps is abandoned. The data are kept unmerged down to the Harker construction, where the symmetry-breaking is explicitly modelled and refined and becomes a source of supplementary phase information. Data merging is effectively carried out on the complex plane, *i.e.* through the generalized Harker construction: from all the data items *F _{j}*(

Anomalous scattering is the best-known example of a symmetry-breaking effect (Bijvoet, 1954 ) that can be exploited for phase determination (Okaya *et al.*, 1955 ; Blow & Rossmann, 1961 ; North, 1965 ; Matthews, 1966 ). The structure factor of a single atom is given by

where *T*(**h**) denotes the Debye–Waller factor and all other symbols have their usual meanings. If the anomalous scattering factor *f*′′ is negligibly small, we can write Friedel’s law as

where * denotes complex conjugation and the pair of reflections **h** and −**h** are called Friedel opposites. Let us denote *f*
^{+}(**h**) = *f*(**h**) and *f*
^{−}(**h**) = *f**(−**h**), so that Friedel’s law is written as

Anomalous scattering breaks the Friedel (or Laue group) equivalence between the reflections **h** and −**h** (Fig. 5 ) and can therefore give rise to measurable intensity differences [except for centric reflections, where *f*(**h**) = *f*(−**h**) even in the presence of anomalous scattering]. Again, it is of importance that only the atoms of the substructure are significantly affected by anomalous scattering. If all atoms display similar anomalous scattering, then it is impossible to generate phase information through this effect^{1}. In the case of anomalous scattering, we can therefore write

where **P**
^{+}(**h**) = **P**(**h**), **P**
^{−}(**h**) = **P***(−**h**), **H**
^{+}(**h**) = **H**(**h**) and **H**
^{−}(**h**) = **H***(−**h**). Thus, the Harker construction is actually set up with **H**
^{+}(**h**) and **H**
^{−}(**h**) as the centres of the circles whose radii correspond to the experimentally observed *F*(**h**) and *F*(−**h**) data (Blow & Rossmann, 1961 ; North, 1965 ).

Friedel symmetry-breaking arising from anomalous scattering. (*a*) Structure factor of an atom which has negligible *f*′′: Friedel’s law is satisfied. (*b*) Structure factor of an atom exhibiting anomalous scattering: *f*
**...**

It is important to note that the Friedel symmetry-breaking effect of anomalous scattering is distinct from its wavelength-dependence, which was mentioned earlier, although in the MAD method both effects are exploited for phase determination. In contrast, in a single-wavelength anomalous dispersion (SAD) experiment only one data set is recorded from the same sample at a single wavelength. The data are left partially unmerged in the sense that Friedel opposites are kept separate. More specifically, the reflections are merged according to the crystal point group, not its Laue group.

A number of case studies have established that X-ray-induced damage to crystalline protein samples initially occurs at discrete well localized sites (Burmeister, 2000 ; Ravelli & McSweeney, 2000 ; Weik *et al.*, 2000 ), leading to the breakage of disulfide bonds, the decarboxylation of acidic residues, the loss of hydroxyl groups from tyrosines and the loss of methylthio groups from methionines. Heavier atoms such as selenium in selenosubstituted proteins (Rice *et al.*, 2000 ; Ravelli *et al.*, 2005 ), bromine in brominated nucleic acids (Ennifar *et al.*, 2002 ; Ravelli *et al.*, 2003 ; Schiltz *et al.*, 2004 ), metals in metalloproteins (Penner-Hahn *et al.*, 1989 ; Schlichting *et al.*, 2000 ; Berglund *et al.*, 2002 ; Yano *et al.*, 2005 ) and iodine (Evans *et al.*, 2003 ; Zwart *et al.*, 2004 ) and mercury (Ramagopal *et al.*, 2005 ) in isomorphous derivatives often exhibit a particularly pronounced sensitivity to X-ray damage.

The structure factor of a radiation-sensitive atom changes continuously with time or, more precisely, as a function of the X-ray dose *d*. In most cases, site-specific radiation damage simply leads to an increased disorder of the atom so that we can write its structure factor as

where is the zero-dose occupancy of the atom and ξ(*d*) is a continuous function of *d* that varies between the values ξ(0) = 0 and ξ(∞) = 1.

In the presence of site-specific radiation damage, symmetry-related reflections or repeated measurements of the same reflection that are recorded at different X-ray doses will no longer be equivalent since they pertain to different stages of the radiation-induced changes in the substructure. When these data are kept unmerged, and provided that the site-specific modulations of the substructure [*i.e.* the function ξ(*d*) for each atom of the substructure] can be modelled, it is possible to generate phase information from the observed intensity differences *via* the generalized Harker construction.

In the case of the crystal structure determination of Br-DIS, a brominated 23-nucleotide RNA fragment, standard three-wavelength MAD phasing was unsuccessful because of rapid X-ray-induced debromination (Ennifar *et al.*, 2002 ). We demonstrated that a substantial enhancement of the phasing power was achieved by modelling the site-specific changes in a continuous dose-dependent fashion (Fig. 6 ) and by keeping the diffraction data unmerged (Schiltz *et al.*, 2004 ). The evolution of site-specific radiation damage was explicitly modelled in real space through a very simple exponential decay model with only one refineable parameter, β, for each of the two independent Br atoms,

In this example, exploiting the symmetry-breaking effects of site-specific radiation damage enabled the successful phasing of a previously resistant problem (Fig. 7 ).

Other situations can arise in which the positions and/or the scattering factors of the substructure atoms also vary as a function of X-ray dose. More elaborate models that are appropriate for such cases have been discussed in Schiltz & Bricogne (2007 ).

The anomalous scattering terms for isolated atoms are scalars but, generally, chemical bonding and the symmetry of the atom’s environment induce a directional dependence of *f*′ and *f*′′ on the direction of linear polarization of the X-ray beam. This anisotropy of anomalous scattering (AAS) is significant in the near-edge region of absorption maxima, as has been experimentally demonstrated in numerous investigations on inorganic and small-molecule compounds (Templeton & Templeton, 1982 , 1988 , 1995 ; Kirfel *et al.*, 1991 ; Dmitrienko, 1983 ; Dmitrienko *et al.*, 2005 ). AAS has also been observed in selenated proteins (Hendrickson *et al.*, 1990 ; Fanchon & Hendrickson, 1990 ; Bricogne *et al.*, 2005 ; Schiltz & Bricogne, 2008 ), in metalloproteins (Hendrickson *et al.*, 1988 ) and in brominated nucleotides (Bricogne *et al.*, 2005 ; Sanishvili *et al.*, 2007 ; Oliéric *et al.*, 2007 ). AAS arises from resonant transitions between the core electrons and antibonding valence molecular orbitals that are rendered nonspherically symmetric by the chemical bonding of the absorbing atom. The anomalous scattering thus depends on the relative orientation of the electric field vector of the incident X-ray beam (the polarization direction) with respect to these molecular orbitals. This is illustrated in Figs. 8 and 9 , respectively, for Se in selenomethionine and Br in brominated nucleotides, which are the two most important anomalous scatterers used for phasing purposes in macromolecular crystallography. The variations in *f*′ and *f*′′ as a function of molecular orientation are very substantial in the near-edge region and are of the same order of magnitude as those that can be obtained by changing the wavelength of the incident beam in a MAD experiment.

Anomalous scattering factors *f*′ and *f*′′ for Se in selenomethionine residues. The curves represent the anomalous scattering factors when the polarization direction of the incident X-ray beam is aligned with one of the molecular **...**

Anomalous scattering factors *f*′ and *f*′′ for Br in brominated nucleotides. The curves represent the anomalous scattering factors when the polarization direction of the incident X-ray beam is aligned with one of the molecular principal **...**

In brominated nucleotides, the white line arises through a resonant transition of a 1*s* core electron to an empty antibonding σ* molecular orbital, which is a linear combination of mainly the C 2*sp*
_{z}
^{2} and Br 4*p*
_{z} orbitals (*z* being oriented along the C—Br bond) and therefore has a pronounced *p*
_{z} symmetry (Sanishvili *et al.*, 2007 ). In accordance, in the experimental absorption spectra the white line is observed to be most pronounced along the direction parallel to the C—Br bond, whereas it completely disappears when the polarization vector of the X-ray beam is perpendicular to the C—Br bond. Concomitantly, a large shift of more than 7 eV in the energy position of the minimum (the so-called *inflection point*) is observed in the *f*′ spectra when the direction of polarization is changed from parallel to perpendicular to the C—Br bond.

Synchrotron X-rays are linearly polarized in the plane defined by the orbit of the electron beam, *i.e.* in the horizontal plane. In principle, the effects of AAS can be revealed by variations in the near-edge absorption spectra and by variations in the diffracted intensities that occur upon changing the orientation of the crystal with respect to the direction of polarization of the X-ray beam. However, at protein crystallography beamlines a single rotation (spindle) axis is usually employed and this axis is almost universally horizontally oriented, *i.e.* exactly parallel to the direction of polarization of the X-ray beam. Thus, as the crystal is rotated during a data collection the direction of polarization of the incident beam does not change with respect to the crystal. Nevertheless, the polarization of the incident beam can break the crystal symmetry. Since symmetry-related anomalously scattering atoms may experience the incident electric field under different polarization orientations, they are not necessarily equivalent any longer as far as their scattering amplitudes are concerned.

The symmetry-breaking effects of AAS were cogently demonstrated in the case of a brominated Z-DNA duplex d(CGCG[BrU]G) (Sanishvili *et al.*, 2007 ; Schiltz & Bricogne, 2008 ; Fig. 10 ). X-ray diffraction data were recorded from a single cryocooled crystal at a wavelength corresponding to the position of the Br *K*-edge white line in brominated nucleotides. The data were recorded using a single scan axis oriented parallel to the direction of polarization of the X-ray beam. Hence, during the data collection the direction of polarization remained constant with respect to the crystal and thus also with respect to the C—Br bonds. However, the various symmetry-related C—Br bonds ‘saw’ the polarization direction of the incident beam from different relative orientations and thus gave rise to different effective anomalous scattering factors. The Br atoms with C—Br bonds nearly aligned with the polarization direction exhibited a strong white line (large *f*′′), but this was not the case for symmetry-related Br atoms, where the C—Br bonds were oriented more closely perpendicular to the polarization direction. In a sense, the symmetry operations acted as an ‘internal’ goniometer allowing the AAS properties of symmetry-related sites to be sampled at different orientations, even though the orientation of the crystal with respect to the polarization direction remained fixed. In an anomalous difference Fourier map computed with the data merged in point group 1 (*i.e. not* imposing any symmetry), symmetry-related Br sites displayed widely differing peak heights. A clear correlation can be established between the height of each peak and the angle between the corresponding C—Br bond direction and the direction of X-ray polarization.

Similar AAS-induced symmetry-breaking effects were also observed in selenated protein crystals (Schiltz & Bricogne, 2008 ).

When reflection data are merged in a certain point or Laue group one actually imposes a symmetry on the crystal structure, so that any genuine intensity differences between symmetry-related reflections are averaged out. Therefore, the widespread practice of merging data prior to phasing completely scrambles the effects of AAS and this may explain why conventional SAD or MAD phasing strategies have not substantially suffered from ignoring AAS altogether. On the other hand, if the data are kept unmerged the intensity differences in symmetry-related reflections can be exploited to model the AAS of anomalously scattering atoms and to generate phase information. Since data sets are usually recorded with a certain degree of redundancy, this additional phase information essentially comes ‘for free’. This can be of particular use in resolving the phase ambiguity in SAD experiments. In Schiltz & Bricogne (2008 ), examples were presented of improvements in the quality of phases which were typically of the same order of magnitude as those obtained in a conventional approach by adding a second-wavelength data set to a SAD experiment. Thus, the exploitation of AAS can give access to two-wavelength map quality with single-wavelength measurements. Such a gain is particularly significant, since radiation damage can frequently preclude the collection of a second-wavelength data set.

While AAS has been extensively studied in physical, inorganic and small-molecule crystallography, the existence of this phenomenon has until recently remained relatively obscure to protein crystallographers (with the notable exception of Fanchon & Hendrickson, 1990 ). Because the details of AAS have not always been well understood, a widely held misconception states that AAS in protein crystals is only observable if the molecular groups which surround the anomalously scattering atoms are all aligned in the crystal. A related fallacy is the belief that AAS in protein crystals of high symmetry and/or containing a large number of anomalously scattering atoms in the asymmetric unit would ‘average out’ to isotropy. This is indeed the case for linear X-ray dichroism (as, for example, observed in polarized absorption spectra) and birefringence, which are global (macroscopic) consequences of AAS and which follow the point-group symmetry of the crystal (Bricogne *et al.*, 2005 ; Sanishvili *et al.*, 2007 ). There can thus be no dichroism in cubic crystals (except for higher order effects, which are usually weak). However, this is clearly not the case for AAS effects in *diffraction*, which are microscopic (local) effects to which each individual atom contributes with its own phase shift. In Schiltz & Bricogne (2008 ), we reported substantial AAS effects in the diffracted intensities from crystals of selenomethionine phosphopantetheine adenylyltransferase (PPAT; Izard & Geerlof, 1999 ), which crystallizes in the cubic space group *I*23 and contains 384 ordered selenomethionine groups in the unit cell.

The simplest model for AAS uses ‘symmetry-unrolled’ anomalous scattering factors (Schiltz & Bricogne, 2008 ) where for each anomalously scattering atom *in the unit cell* individual anomalous scattering factors *f*′ and *f*′′ are defined and can be freely refined. This model is only valid if all reflections have been recorded with the same polarization direction and therefore applies to the particular, but not uncommon, case in which a data set has been collected by rotating about a horizontal spindle axis. It should be noted that this model is not equivalent to refining all the substructure atoms in space group *P*1. Only the anomalous scattering factors of symmetry-related atoms are refined individually; the positional and thermal parameters are constrained to obey the space-group symmetry, *i.e.* the symmetry-breaking effects are assumed to only originate from the anomalous scattering properties.

In a more general model for AAS the anomalous scattering properties of an atom are described by a second-rank tensor **f**, represented by a symmetric 3 × 3 matrix with complex-valued entries (Templeton & Templeton, 1982 ; Dmitrienko, 1983 ),

The AAS tensors of two sites *s* and *k* which are symmetry-related through the space-group operation (**R**
_{g}, **t**
* _{g}*) (such that

where the left superscript *t* stands for matrix transposition.

We have shown in Schiltz & Bricogne (2008 ) that in the context of macromolecular crystallography the anomalous scattering factor of an atom that displays AAS can be approximated by

where it is assumed that the incident beam is completely linearly polarized along a direction given by the unit vector **p**. The unit vector **p**′ is obtained by projecting **p** onto a plane perpendicular to the scattered-beam direction and corresponds to the direction of linear polarization of the diffracted beam in the absence of AAS (Schiltz & Bricogne, 2008 , 2009 ). Unless the polarization direction **p** is oriented along a symmetry axis, the similarity transformation (11) does not in general give rise to identical scattering factors *f* for symmetry-related sites. The symmetry-breaking effects of AAS can thus be properly modelled by adopting a tensor description for the anomalous scattering factors.

The models for site-specific radiation damage and AAS outlined in the previous sections have been implemented in the heavy-atom refinement and phasing program *SHARP* (de La Fortelle & Bricogne, 1997 ; Bricogne *et al.*, 2003 ). The program has been extended for the use of unmerged data and is able to read and processes X-ray dose information for each reflection measurement as well as goniometric information in various forms (Schiltz *et al.*, 200; Schiltz & Bricogne, 2009 ). By applying the generalized Harker construction, symmetry-related reflections and repeated measurements of the same reflection can be used, together with data recorded at another wavelength and/or from another heavy-atom derivative or native.

Since the intensity differences between symmetry-related reflections are usually not very large (of the same order of magnitude as Friedel differences), data scaling can be carried out along conventional lines, *i.e.* by minimizing the disagreement between symmetry-related reflections in the crystal Laue group. The parameters of the substructure atoms are usually first refined in conventional mode, *i.e.* without modelling site-specific radiation damage or AAS. Once the refinement of the positional parameters of the substructure atoms has converged, it is possible to switch on the refinement of dose-dependent or AAS parameters.

The modelling and parametrization of non-isomorphism in the case of data affected by site-specific radiation damage or AAS is significantly more complex than for standard cases. The error model that is currently implemented in *SHARP* assumes that the effects of all sources of non-isomorphism are uncorrelated between different observations of a given reflection (de La Fortelle & Bricogne, 1997 ; Bricogne *et al.*, 2003 ). In essence, a diagonal approximation is used for the non-isomorphism covariance matrix. Such an approximation may not always be entirely justified since non-isomorphism can be correlated across observations that have been recorded under similar geometric conditions (Bricogne *et al.*, 2003 ). For a more general treatment it will be necessary to resort to multivariate likelihood functions which are capable of accommodating adequate patterns of covariances between the various observations (Bricogne, 2000 ; Pannu *et al.*, 2003 ). The implementation of these functions in *SHARP* is currently under way.

Although anomalous scattering, site-specific radiation damage and AAS in the substructure are different phenomena with different physical origins, they have the common property of inducing intensity differences between symmetry-related reflections in a diffraction experiment. Provided that these effects are included in a parametrized model of the substructure, they can become a source of phase information. With the current practice of recording diffraction data using the single-axis rotation method, it is almost always the case that data sets with a certain degree of redundancy (multiplicity) are collected. Redundancy is often achieved before completeness, *i.e.* several symmetry-equivalent observations of certain reflections are recorded while for other reflections no observations have yet been measured. Thus, redundancy is usually a byproduct of striving to collect a complete data set. In this sense, the additional phase information that may be contained in the intensity differences between symmetry-related reflections essentially comes ‘for free’. Since overall radiation damage is in many cases the main limiting factor in the amount of data that can be collected from a single sample for the purpose of experimental phasing, it is of the utmost importance to be able to derive the maximum amount of phase information from this limited amount of data. A current limitation in the implementation of these methods is the approximate treatment of correlated sources of non-isomorphism.

Although in many cases the standard anomalous signal generated through Bijvoet differences is likely to be the main source of phase information, this can be supplemented by exploiting the symmetry-breaking effects in unmerged data as described above. In particular cases, such as that of the brominated RNA fragment described earlier, site-specific and overall radiation damage evolve on significantly different timescales. The former can then become a very significant source of phase information to complement the anomalous phasing signal. However, such favourable cases are rather atypical. In many ‘real-life’ situations overall radiation damage unfolds at a rate that is not significantly different in comparison to the evolution of site-specific radiation damage or in comparison to the total time that is required to record a complete data set. In such cases the quality of the phases will ultimately be limited by the effects of overall radiation damage, although the proper modelling and exploitation of site-specific radiation damage can still yield a noticeable improvement of phases, as was for instance the case in the structure determination of the PP2A phosphatase activator Ypa2 (Leulliot *et al.*, 2006 ; Schiltz & Bricogne, 2007 ).

The general question then arises of how to design a data-collection strategy that enables the optimal exploitation of the various possible sources of phase information when the lifetime of the crystal is limited. Crystals can be intentionally misaligned in order to maximize the AAS-induced inequivalence between symmetry-related reflections. However, since the standard anomalous signal is the most important source of phase information, the reduction of systematic errors in Bijvoet intensity differences is of prime importance. In crystals belonging to high-symmetry space groups, Bijvoet pairs are usually recorded in close temporal proximity, even if the crystal is not specially aligned. In such cases there will always be some symmetries that are broken by AAS and these effects can then be exploited for additional phase information. For crystals of lower symmetry the deliberate alignment of a symmetry axis along the spindle allows the simultaneous recording of Bijvoet pairs. However, if the spindle is oriented horizontally (*i.e.* along the direction of linear polarization of the synchrotron beam) such a geometry will partly or fully neutralize the symmetry-breaking effects of AAS. A more ideal geometry, which would minimize systematic errors in Bijvoet intensity differences by *aligning* a symmetry axis *with respect to the spindle axis*, while at the same time maximizing the symmetry-breaking effects of AAS by *misaligning* the symmetry axis *with respect to the direction of X-ray polarization*, can be achieved with multi-axis goniometers, where the scan axis is not constrained to be aligned with the X-ray polarization direction. Thus, with the future use of goniometers with a vertical scan axis designed for the purposes of gaining mechanical stability in the handling of microcrystals, the effects of AAS will become truly ubiquitous in all data sets collected at an absorption edge of a covalently bonded anomalous scatterer and their proper treatment in experimental phasing will be imperative if a major waste of phase information is to be avoided.

We thank P. Pattison, R. Sanishvili and the members of Global Phasing Ltd for their many contributions at various stages of this project. We acknowledge financial support from the Swiss National Science Foundation Grant No. 200021-101519/1.

Even in the absence of symmetry-breaking effects, the complex structure factors of symmetry-equivalent reflections are, in general, not identical; only their moduli are. Some care therefore has to be exercised when setting up the Harker construction with symmetry-related reflections. All complex structure factors have to be transformed to a common unique representative. The situation is similar to that encountered in building up the Harker construction with Friedel pairs in the case of anomalous scattering. In that case, the Harker construction is actually set up with **F**(**h**) and **F***(−**h**), *i.e.*
**F**(−**h**) is reflected through the real axis on the complex plane (Blow & Rossmann, 1961 ; North, 1965 ).

Let denote the space group of the crystal. The operation of an element *g* of will be written as

Reflections **h**
* _{u}* and

The set of all reflections **h**
* _{v}* which satisfy (14) is called the

In the absence of symmetry-breaking effects, the structure factors of symmetry-equivalent reflections are related by (Waser, 1955 )

If the symmetry-related reflections **h** and **h**
_{g} are to be used on the same Harker construction, it is thus necessary to rotate **F**(**h**
_{g}) back to the phase angle of **F**(**h**), *i.e.* undo the phase shift exp(−2π*i*
**h**·**t**
* _{g}*) and therefore define

In addition, we also want to include the Friedel opposite −**h** and its symmetry-related reflections on the Harker construction. For acentric reflections, **h** and −**h** form distinct orbits, which we label by σ = ±, respectively. We therefore extend our definition to

In the absence of symmetry-breaking effects, all defined in this way will be identical for all and for any sign σ = {+, −} (strict symmetry equivalence, including Friedel equivalence). If we assume, as before, that the unknown part of the structure factor **P**(**h**) is unaffected by any symmetry-breaking effects, we can write

Isotropic anomalous scattering in the substructure breaks the Friedel symmetry of acentric reflections, so that

but otherwise the symmetry equivalence is preserved,

In the presence of AAS or site-specific radiation damage, all are potentially different for the various and signs σ = {+, −}.

In conclusion then, the generalized Harker construction is set up with circles whose radii correspond to observed amplitudes *F*
_{σ,g}(**h**) of all the reflections which are symmetry-related (including Friedel symmetry) to **h**. Each circle is centred at a position given by (18). In an ideal (error-free) case, all the circles will then intersect at the point **P**(**h**).

- Berglund, G. I., Carlsson, G. H., Smith, A. T., Szöke, H., Henriksen, A. & Hajdu, J. (2002).
*Nature (London)*,**417**, 463–468. [PubMed] - Bijvoet, J. (1954).
*Nature (London)*,**173**, 888–891. - Blow, D. M. & Rossmann, M. G. (1961).
*Acta Cryst.***14**, 1195–1202. - Bricogne, G. (2000).
*Proceedings of the Workshop on Advanced Special Functions and Applications, Melfi (PZ), Italy, 9–12 May 1999*, edited by D. Cocolicchio, G. Dattoli & H. M. Srivastava, pp. 315–321. Rome: Arcane Editrice. - Bricogne, G., Capelli, S. C., Evans, G., Mitschler, A., Pattison, P., Roversi, P. & Schiltz, M. (2005).
*J. Appl. Cryst.***38**, 168–182. - Bricogne, G., Vonrhein, C., Flensburg, C., Schiltz, M. & Paciorek, W. (2003).
*Acta Cryst.*D**59**, 2023–2030. [PubMed] - Burmeister, W. P. (2000).
*Acta Cryst.*D**56**, 328–341. [PubMed] - Dmitrienko, V. E. (1983).
*Acta Cryst.*A**39**, 29–35. - Dmitrienko, V. E., Ishida, K., Kirfel, A. & Ovchinnikova, E. N. (2005).
*Acta Cryst.*A**61**, 481–493. [PubMed] - Ennifar, E., Carpentier, P., Ferrer, J.-L., Walter, P. & Dumas, P. (2002).
*Acta Cryst.*D**58**, 1262–1268. [PubMed] - Evans, G., Polentarutti, M., Djinovic Carugo, K. & Bricogne, G. (2003).
*Acta Cryst.*D**59**, 1429–1434. [PubMed] - Fanchon, E. & Hendrickson, W. A. (1990).
*Acta Cryst.*A**46**, 809–820. [PubMed] - Flack, H. D. & Shmueli, U. (2007).
*Acta Cryst.*A**63**, 257–265. [PubMed] - Green, D., Ingram, V. & Perutz, M. (1954).
*Proc. R. Soc. Lond. Ser. A*,**225**, 287–307. - Harker, D. (1956).
*Acta Cryst.***9**, 1–9. - Hendrickson, W. A. (1991).
*Science*,**254**, 51–58. [PubMed] - Hendrickson, W. A., Horton, J. R. & LeMaster, D. M. (1990).
*EMBO J.***9**, 1665–1672. [PubMed] - Hendrickson, W. A., Smith, J. L., Phizackerley, R. P. & Merritt, E. A. (1988).
*Proteins*,**4**, 77–88. [PubMed] - Izard, T. & Geerlof, A. (1999).
*EMBO J.***18**, 2021–2030. [PubMed] - Kirfel, A., Petcov, A. & Eichhorn, K. (1991).
*Acta Cryst.*A**47**, 180–195. - La Fortelle, E. de & Bricogne, G. (1997).
*Methods Enzymol.***276**, 472–494. - Leulliot, N., Vicentini, G., Jordene, J., Cheruel, S., Liger, D., Schiltz, M., Van Tilbeurgh, H., Barford, D. & Goris, J. (2006).
*Mol. Cell*,**23**, 413–424. [PubMed] - Matthews, B. W. (1966).
*Acta Cryst.***20**, 82–86. - Nanao, M. & Ravelli, R. B. G. (2006).
*Structure*,**14**, 791–800. [PubMed] - North, A. C. T. (1965).
*Acta Cryst.***18**, 212–216. - Okaya, Y., Saito, Y. & Pepinsky, R. (1955).
*Phys. Rev.***98**, 1857–1858. - Oliéric, V., Ennifar, E., Meents, A., Fleurant, M., Besnard, C., Pattison, P., Schiltz, M., Schulze-Briese, C. & Dumas, P. (2007).
*Acta Cryst.*D**63**, 759–768. [PubMed] - Pannu, N. S., McCoy, A. J. & Read, R. J. (2003).
*Acta Cryst.*D**59**, 1801–1808. [PubMed] - Penner-Hahn, J. E., Murata, M., Hodgson, K. O. & Freeman, H. C. (1989).
*Inorg. Chem.***28**, 1826–1832. - Ramagopal, U. A., Dauter, Z., Thirumuruhan, R., Fedorov, E. & Almo, S. C. (2005).
*Acta Cryst.*D**61**, 1289–1298. [PubMed] - Ravelli, R. B. G. & McSweeney, S. M. (2000).
*Structure*,**8**, 315–328. [PubMed] - Ravelli, R. B. G., Nanao, M. H., Lovering, A., White, S. & McSweeney, S. (2005).
*J. Synchrotron Rad.***12**, 276–284. [PubMed] - Ravelli, R. B. G., Schrøder Leiros, H.-K., Pan, B., Caffrey, M. & McSweeney, S. M. (2003).
*Structure*,**11**, 217–224. [PubMed] - Rice, L. M., Earnest, T. N. & Brunger, A. T. (2000).
*Acta Cryst.*D**56**, 1413–1420. [PubMed] - Sanishvili, R., Besnard, C., Camus, F., Fleurant, M., Pattison, P., Bricogne, G. & Schiltz, M. (2007).
*J. Appl. Cryst.***40**, 552–558. - Schiltz, M. & Bricogne, G. (2007).
*J. Synchrotron Rad.***14**, 34–42. [PubMed] - Schiltz, M. & Bricogne, G. (2008).
*Acta Cryst.*D**64**, 711–729. [PubMed] - Schiltz, M. & Bricogne, G. (2009).
*J. Appl. Cryst.***42**, 101–108. - Schiltz, M., Dumas, P., Ennifar, E., Flensburg, C., Paciorek, W., Vonrhein, C. & Bricogne, G. (2004).
*Acta Cryst.*D**60**, 1024–1031. [PubMed] - Schlichting, I., Berendzen, J., Chu, K., Stock, A. M., Maves, S. A., Benson, D. E., Sweet, R. M., Ringe, D., Petsko, G. A. & Sligar, S. G. (2000).
*Science*,**287**, 1615–1622. [PubMed] - Templeton, D. H. & Templeton, L. K. (1982).
*Acta Cryst.*A**38**, 62–67. - Templeton, L. K. & Templeton, D. H. (1988).
*Acta Cryst.*A**44**, 1045–1051. [PubMed] - Templeton, L. K. & Templeton, D. H. (1995).
*J. Synchrotron Rad.***2**, 31–35. [PubMed] - Waser, J. (1955).
*Acta Cryst.***8**, 595. - Weik, M., Ravelli, R. B. G., Kryger, G., McSweeney, S., Raves, M. L., Harel, M., Gros, P., Silman, I., Kroon, J. & Sussman, J. L. (2000).
*Proc. Natl Acad. Sci. USA*,**97**, 623–628. [PubMed] - Yano, J., Kern, J., Irrgang, K.-D., Latimer, M. J., Bergmann, U., Glatzel, P., Pushkar, Y., Biesiadka, J., Loll, B., Sauer, K., Messinger, J., Zouni, A. & Yachandra, V. (2005).
*Proc. Natl Acad. Sci. USA*,**102**, 12047–12052. [PubMed] - Zwart, P. H., Banumathi, S., Dauter, M. & Dauter, Z. (2004).
*Acta Cryst.*D**60**, 1958–1963. [PubMed]

Articles from Acta Crystallographica Section D: Biological Crystallography are provided here courtesy of **International Union of Crystallography**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's Canada Institute for Scientific and Technical Information in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |