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

A modified phase-retrieval algorithm to facilitate automatic de novo macromolecular structure determination in single-wavelength anomalous diffraction

The success of experimental phasing in macromolecular crystallography relies primarily on the accurate locations of heavy atoms bound to the target crystal. To improve the process of substructure determination, a modified phase-retrieval algorithm built on the framework of the relaxed alternating averaged reflection (RAAR) algorithm has been developed. Importantly, the proposed algorithm features a combination of the π-half phase perturbation for weak reflections and enforces the direct-method-based tangent formula for strong reflections in reciprocal space. The proposed algorithm is extensively demonstrated on a total of 100 single-wavelength anomalous diffraction (SAD) experimental datasets, comprising both protein and nucleic acid structures of different qualities. Compared with the standard RAAR algorithm, the modified phase-retrieval algorithm exhibits significantly improved effectiveness and accuracy in SAD substructure determination, highlighting the importance of additional constraints for algorithmic performance. Furthermore, the proposed algorithm can be performed without human intervention under most conditions owing to the self-adaptive property of the input parameters, thus making it convenient to be integrated into the structural determination pipeline. In conjunction with the IPCAS software suite, we demonstrated experimentally that automatic de novo structure determination is possible on the basis of our proposed algorithm.

Experimental phasing opportunities for macromolecular crystallography at very long wavelengths

Despite recent advances in cryo-electron microscopy and artificial intelligence-based model predictions, a significant fraction of structure determinations by macromolecular crystallography still requires experimental phasing, usually by means of single-wavelength anomalous diffraction (SAD) techniques. Most synchrotron beamlines provide highly brilliant beams of X-rays of between 0.7 and 2 Å wavelength. Use of longer wavelengths to access the absorption edges of biologically important lighter atoms such as calcium, potassium, chlorine, sulfur and phosphorus for native-SAD phasing is attractive but technically highly challenging. The long-wavelength beamline I23 at Diamond Light Source overcomes these limitations and extends the accessible wavelength range to λ = 5.9 Å. Here we report 22 macromolecular structures solved in this extended wavelength range, using anomalous scattering from a range of elements which demonstrate the routine feasibility of lighter atom phasing. We suggest that, in light of its advantages, long-wavelength crystallography is a compelling option for experimental phasing.

Obtaining anomalous and ensemble information from protein crystals from 220 K up to physiological temperatures

X-ray crystallography has been invaluable in delivering structural information about proteins. Previously, an approach has been developed that allows high-quality X-ray diffraction data to be obtained from protein crystals at and above room temperature. Here, this previous work is built on and extended by showing that high-quality anomalous signal can be obtained from single protein crystals using diffraction data collected at 220 K up to physiological temperatures. The anomalous signal can be used to directly determine the structure of a protein, i.e. to phase the data, as is routinely performed under cryoconditions. This ability is demonstrated by obtaining diffraction data from model lysozyme, thaumatin and proteinase K crystals, the anomalous signal from which allowed their structures to be solved experimentally at 7.1 keV X-ray energy and at room temperature with relatively low data redundancy. It is also demonstrated that the anomalous signal from diffraction data obtained at 310 K (37°C) can be used to solve the structure of proteinase K and to identify ordered ions. The method provides useful anomalous signal at temperatures down to 220 K, resulting in an extended crystal lifetime and increased data redundancy. Finally, we show that useful anomalous signal can be obtained at room temperature using X-rays of 12 keV energy as typically used for routine data collection, allowing this type of experiment to be carried out at widely accessible synchrotron beamline energies and enabling the simultaneous extraction of high-resolution data and anomalous signal. With the recent emphasis on obtaining conformational ensemble information for proteins, the high resolution of the data allows such ensembles to be built, while the anomalous signal allows the structure to be experimentally solved, ions to be identified, and water molecules and ions to be differentiated. Because bound metal-, phosphorus- and sulfur-containing ions all have anomalous signal, obtaining anomalous signal across temperatures and up to physiological temperatures will provide a more complete description of protein conformational ensembles, function and energetics.

Native SAD phasing at room temperature

Single-wavelength anomalous diffraction (SAD) is a routine method for overcoming the phase problem when solving macromolecular structures. This technique requires the accurate measurement of intensities to determine differences between Bijvoet pairs. Although SAD experiments are commonly conducted at cryogenic temperatures to mitigate the effects of radiation damage, such temperatures can alter the conformational ensemble of the protein and may impede the merging of data from multiple crystals due to non-uniform freezing. Here, a strategy is presented to obtain high-quality data from room-temperature, single-crystal experiments. To illustrate the strengths of this approach, native SAD phasing at 6.55 keV was used to solve four structures of three model systems at 295 K. The resulting data sets allow automatic phasing and model building, and reveal alternate conformations that reflect the structure of proteins at room temperature.

Likelihood-based estimation of substructure content from single-wavelength anomalous diffraction (SAD) intensity data

SAD phasing can be challenging when the signal-to-noise ratio is low. In such cases, having an accurate estimate of the substructure content can determine whether or not the substructure of anomalous scatterer positions can successfully be determined. Here, a likelihood-based target function is proposed to accurately estimate the strength of the anomalous scattering contribution directly from the measured intensities, determining a complex correlation parameter relating the Bijvoet mates as a function of resolution. This gives a novel measure of the intrinsic anomalous signal. The SAD likelihood target function also accounts for correlated errors in the measurement of intensities from Bijvoet mates, which can arise from the effects of radiation damage. When the anomalous signal is assumed to come primarily from a substructure comprising one anomalous scatterer with a known value of f'' and when the protein composition of the crystal is estimated correctly, the refined complex correlation parameters can be interpreted in terms of the atomic content of the primary anomalous scatterer before the substructure is known. The maximum-likelihood estimation of substructure content was tested on a curated database of 357 SAD cases with useful anomalous signal. The prior estimates of substructure content are highly correlated to the content determined by phasing calculations, with a correlation coefficient (on a log-log basis) of 0.72.

Phosphorus and sulfur SAD phasing of the nucleic acid-bound DNA-binding domain of interferon regulatory factor 4

Solution of the structure of the DNA-binding domain of interferon regulatory factor 4 bound to its interferon-stimulated response element by native intrinsic phosphorus and sulfur single-wavelength anomalous dispersion methods (native SAD) is described.

Pivotal to the regulation of key cellular processes such as the transcription, replication and repair of DNA, DNA-binding proteins play vital roles in all aspects of genetic activity. The determination of high-quality structures of DNA-binding proteins, particularly those in complexes with DNA, provides crucial insights into the understanding of these processes. The presence in such complexes of phosphate-rich oligonucleotides offers the choice of a rapid method for the routine solution of DNA-binding proteins through the use of long-wavelength beamlines such as I23 at Diamond Light Source. This article reports the use of native intrinsic phosphorus and sulfur single-wavelength anomalous dispersion methods to solve the complex of the DNA-binding domain (DBD) of interferon regulatory factor 4 (IRF4) bound to its interferon-stimulated response element (ISRE). The structure unexpectedly shows three molecules of the IRF4 DBD bound to one ISRE. The sole reliance on native intrinsic anomalous scattering elements that belong to DNA–protein complexes renders the method of general applicability to a large number of such protein complexes that cannot be solved by molecular replacement or by other phasing methods.

Native de novo structural determinations of non-canonical nucleic acid motifs by X-ray crystallography at long wavelengths

Obtaining phase information remains a formidable challenge for nucleic acid structure determination. The introduction of an X-ray synchrotron beamline designed to be tunable to long wavelengths at Diamond Light Source has opened the possibility to native de novo structure determinations by the use of intrinsic scattering elements. This provides opportunities to overcome the limitations of introducing modifying nucleotides, often required to derive phasing information. In this paper, we build on established methods to generate new tools for nucleic acid structure determinations. We report on the use of (i) native intrinsic potassium single-wavelength anomalous dispersion methods (K-SAD), (ii) use of anomalous scattering elements integral to the crystallization buffer (extrinsic cobalt and intrinsic potassium ions), (iii) extrinsic bromine and intrinsic phosphorus SAD to solve complex nucleic acid structures. Using the reported methods we solved the structures of (i) Pseudorabies virus (PRV) RNA G-quadruplex and ligand complex, (ii) PRV DNA G-quadruplex, and (iii) an i-motif of human telomeric sequence. Our results highlight the utility of using intrinsic scattering as a pathway to solve and determine non-canonical nucleic acid motifs and reveal the variability of topology, influence of ligand binding, and glycosidic angle rearrangements seen between RNA and DNA G-quadruplexes of the same sequence.

Beyond integration: modeling every pixel to obtain better structure factors from stills

Most crystallographic data processing methods use pixel integration. In serial femtosecond crystallography (SFX), the intricate interaction between the reciprocal lattice point and the Ewald sphere is integrated out by averaging symmetrically equivalent observations recorded across a large number (104-106) of exposures. Although sufficient for generating biological insights, this approach converges slowly, and using it to accurately measure anomalous differences has proved difficult. This report presents a novel approach for increasing the accuracy of structure factors obtained from SFX data. A physical model describing all observed pixels is defined to a degree of complexity such that it can decouple the various contributions to the pixel intensities. Model dependencies include lattice orientation, unit-cell dimensions, mosaic structure, incident photon spectra and structure factor amplitudes. Maximum likelihood estimation is used to optimize all model parameters. The application of prior knowledge that structure factor amplitudes are positive quantities is included in the form of a reparameterization. The method is tested using a synthesized SFX dataset of ytterbium(III) lysozyme, where each X-ray laser pulse energy is centered at 9034 eV. This energy is 100 eV above the Yb3+ L-III absorption edge, so the anomalous difference signal is stable at 10 electrons despite the inherent energy jitter of each femtosecond X-ray laser pulse. This work demonstrates that this approach allows the determination of anomalous structure factors with very high accuracy while requiring an order-of-magnitude fewer shots than conventional integration-based methods would require to achieve similar results.

Advances in long-wavelength native phasing at X-ray free-electron lasers

Long-wavelength pulses from the Swiss X-ray free-electron laser (XFEL) have been used for de novo protein structure determination by native single-wavelength anomalous diffraction (native-SAD) phasing of serial femtosecond crystallography (SFX) data. In this work, sensitive anomalous data-quality indicators and model proteins were used to quantify improvements in native-SAD at XFELs such as utilization of longer wavelengths, careful experimental geometry optimization, and better post-refinement and partiality correction. Compared with studies using shorter wavelengths at other XFELs and older software versions, up to one order of magnitude reduction in the required number of indexed images for native-SAD was achieved, hence lowering sample consumption and beam-time requirements significantly. Improved data quality and higher anomalous signal facilitate so-far underutilized de novo structure determination of challenging proteins at XFELs. Improvements presented in this work can be used in other types of SFX experiments that require accurate measurements of weak signals, for example time-resolved studies.

Experimental phasing with vanadium and application to nucleotide-binding membrane proteins

The structure determination of soluble and membrane proteins can be hindered by the crystallographic phase problem, especially in the absence of a suitable homologous structure. Experimental phasing is the method of choice for novel structures; however, it often requires heavy-atom derivatization, which can be difficult and time-consuming. Here, a novel and rapid method to obtain experimental phases for protein structure determination by vanadium phasing is reported. Vanadate is a transition-state mimic of phosphoryl-transfer reactions and it has the advantage of binding specifically to the active site of numerous enzymes catalyzing this reaction. The applicability of vanadium phasing has been validated by determining the structures of three different protein-vanadium complexes, two of which are integral membrane proteins: the rabbit sarcoplasmic reticulum Ca²⁺-ATPase, the antibacterial peptide ATP-binding cassette transporter McjD from Escherichia coli and the soluble enzyme RNAse A from Bos taurus. Vanadium phasing was successful even at low resolution and despite severe anisotropy in the data. This method is principally applicable to a large number of proteins, representing six of the seven Enzyme Commission classes. It relies exclusively on the specific chemistry of the protein and it does not require any modifications, making it a very powerful addition to the phasing toolkit. In addition to the phasing power of this technique, the protein-vanadium complexes also provide detailed insights into the reaction mechanisms of the studied proteins.

Towards the automation of in situ experimental phasing

A new approach to in situ experimental phasing introduced by Lawrence et al. (2020, Acta Cryst. D76, 790–801) will be helpful for the macromolecular crystallography community.

The predictive power of data-processing statistics

This study describes a method to estimate the likelihood of success in determining a macromolecular structure by X-ray crystallography and experimental single-wavelength anomalous dispersion (SAD) or multiple-wavelength anomalous dispersion (MAD) phasing based on initial data-processing statistics and sample crystal properties. Such a predictive tool can rapidly assess the usefulness of data and guide the collection of an optimal data set. The increase in data rates from modern macromolecular crystallography beamlines, together with a demand from users for real-time feedback, has led to pressure on computational resources and a need for smarter data handling. Statistical and machine-learning methods have been applied to construct a classifier that displays 95% accuracy for training and testing data sets compiled from 440 solved structures. Applying this classifier to new data achieved 79% accuracy. These scores already provide clear guidance as to the effective use of computing resources and offer a starting point for a personalized data-collection assistant.

SAD phasing of XFEL data depends critically on the error model

A nonlinear least-squares method for refining a parametric expression describing the estimated errors of reflection intensities in serial crystallographic (SX) data is presented. This approach, which is similar to that used in the rotation method of crystallographic data collection at synchrotrons, propagates error estimates from photon-counting statistics to the merged data. Here, it is demonstrated that the application of this approach to SX data provides better SAD phasing ability, enabling the autobuilding of a protein structure that had previously failed to be built. Estimating the error in the merged reflection intensities requires the understanding and propagation of all of the sources of error arising from the measurements. One type of error, which is well understood, is the counting error introduced when the detector counts X-ray photons. Thus, if other types of random errors (such as readout noise) as well as uncertainties in systematic corrections (such as from X-ray attenuation) are completely understood, they can be propagated along with the counting error, as appropriate. In practice, most software packages propagate as much error as they know how to model and then include error-adjustment terms that scale the error estimates until they explain the variance among the measurements. If this is performed carefully, then during SAD phasing likelihood-based approaches can make optimal use of these error estimates, increasing the chance of a successful structure solution. In serial crystallography, SAD phasing has remained challenging, with the few examples of de novo protein structure solution each requiring many thousands of diffraction patterns. Here, the effects of different methods of treating the error estimates are estimated and it is shown that using a parametric approach that includes terms proportional to the known experimental uncertainty, the reflection intensity and the squared reflection intensity to improve the error estimates can allow SAD phasing even from weak zinc anomalous signal.

Structure determination of the nucleosome core particle by selenium SAD phasing

The eukaryotic genome is compacted inside the nucleus of the cell in the form called chromatin. The fundamental unit of chromatin is the nucleosome, which contains four types of histones (H3, H4, H2A and H2B) and approximately 150 base pairs of DNA wrapped around the histone complex. The structure of the nucleosome is highly conserved across several eukaryotic species, and molecular replacement has been the primary phasing method used to solve nucleosome structures by X-ray crystallography. However, there is currently no simple, widely applicable experimental phasing method for the nucleosome. In the present study, it is demonstrated that selenomethionine-incorporated histones H3, H2A and H2B can be reconstituted into nucleosomes and crystallized for structural determination. Unexpectedly, it was found that the nucleosome can be phased with a relatively small number of Se atoms. The structures of nucleosome core particles containing 12 and 16 Se atoms were solved by SAD phasing at 2.5 and 2.4 Å resolution, respectively. The present study demonstrates a simple method for determining nucleosome structures by experimental phasing, which may be particularly useful for noncanonical structures that cannot be solved by molecular replacement.

Macromolecular structure determination using X-rays, neutrons and electrons: recent developments in Phenix

Diffraction (X-ray, neutron and electron) and electron cryo-microscopy are powerful methods to determine three-dimensional macromolecular structures, which are required to understand biological processes and to develop new therapeutics against diseases. The overall structure-solution workflow is similar for these techniques, but nuances exist because the properties of the reduced experimental data are different. Software tools for structure determination should therefore be tailored for each method. Phenix is a comprehensive software package for macromolecular structure determination that handles data from any of these techniques. Tasks performed with Phenix include data-quality assessment, map improvement, model building, the validation/rebuilding/refinement cycle and deposition. Each tool caters to the type of experimental data. The design of Phenix emphasizes the automation of procedures, where possible, to minimize repetitive and time-consuming manual tasks, while default parameters are chosen to encourage best practice. A graphical user interface provides access to many command-line features of Phenix and streamlines the transition between programs, project tracking and re-running of previous tasks.

Making routine native SAD a reality: lessons from beamline X06DA at the Swiss Light Source

Native single-wavelength anomalous dispersion (SAD) is the most attractive de novo phasing method in macromolecular crystallography, as it directly utilizes intrinsic anomalous scattering from native crystals. However, the success of such an experiment depends on accurate measurements of the reflection intensities and therefore on careful data-collection protocols. Here, the low-dose, multiple-orientation data-collection protocol for native SAD phasing developed at beamline X06DA (PXIII) at the Swiss Light Source is reviewed, and its usage over the last four years on conventional crystals (>50 µm) is reported. Being experimentally very simple and fast, this method has gained popularity and has delivered 45 de novo structures to date (13 of which have been published). Native SAD is currently the primary choice for experimental phasing among X06DA users. The method can address challenging cases: here, native SAD phasing performed on a streptavidin-biotin crystal with P21 symmetry and a low Bijvoet ratio of 0.6% is highlighted. The use of intrinsic anomalous signals as sequence markers for model building and the assignment of ions is also briefly described.

Choosing your (Friedel) mates wisely: grouping data sets to improve anomalous signal

Single-wavelength anomalous diffraction (SAD) phasing from multiple crystals can be especially challenging in samples with weak anomalous signals and/or strong non-isomorphism. Here, advantage is taken of the combinatorial diversity possible in such experiments to study the relationship between merging statistics and downstream metrics of phasing signals. It is furthermore shown that a genetic algorithm (GA) can be used to optimize the grouping of data sets to enhance weak anomalous signals based on these merging statistics.

Cadmium SAD phasing at CuKα wavelength

Single-wavelength anomalous diffraction (SAD) is the most common method for de novo elucidation of macromolecular structures by X-ray crystallography. It requires an anomalous scatterer in a crystal to calculate phases. A recent study by Panneerselvam et al. emphasized the utility of cadmium ions for SAD phasing at the standard synchrotron wavelength of 1 Å. Here we show that cadmium is also useful for phasing of crystals collected in-house with CuKα radiation. Using a crystal of single-domain antibody as an experimental model, we demonstrate how cadmium SAD can be conveniently employed to solve a CuKα dataset. We then discuss the factors which make this method generally applicable.

Fast and accurate data collection for macromolecular crystallography using the JUNGFRAU detector

The accuracy of X-ray diffraction data is directly related to how the X-ray detector records photons. Here we describe the application of a direct-detection charge-integrating pixel-array detector (JUNGFRAU) in macromolecular crystallography (MX). JUNGFRAU features a uniform response on the subpixel level, linear behavior toward high photon rates, and low-noise performance across the whole dynamic range. We demonstrate that these features allow accurate MX data to be recorded at unprecedented speed. We also demonstrate improvements over previous-generation detectors in terms of data quality, using native single-wavelength anomalous diffraction (SAD) phasing, for thaumatin, lysozyme, and aminopeptidase N. Our results suggest that the JUNGFRAU detector will substantially improve the performance of synchrotron MX beamlines and equip them for future synchrotron light sources.

In situ serial crystallography for rapid de novo membrane protein structure determination

De novo membrane protein structure determination is often limited by the availability of large crystals and the difficulties in obtaining accurate diffraction data for experimental phasing. Here we present a method that combines in situ serial crystallography with de novo phasing for fast, efficient membrane protein structure determination. The method enables systematic diffraction screening and rapid data collection from hundreds of microcrystals in in meso crystallization wells without the need for direct crystal harvesting. The requisite data quality for experimental phasing is achieved by accumulating diffraction signals from isomorphous crystals identified post-data collection. The method works in all experimental phasing scenarios and is particularly attractive with fragile, weakly diffracting microcrystals. The automated serial data collection approach can be readily adopted at most microfocus macromolecular crystallography beamlines.

Maximum-likelihood determination of anomalous substructures

A fast Fourier transform (FFT) method is described for determining the substructure of anomalously scattering atoms in macromolecular crystals that allows successful structure determination by X-ray single-wavelength anomalous diffraction (SAD). This method is based on the maximum-likelihood SAD phasing function, which accounts for measurement errors and for correlations between the observed and calculated Bijvoet mates. Proof of principle is shown that this method can improve determination of the anomalously scattering substructure in challenging cases where the anomalous scattering from the substructure is weak but the substructure also constitutes a significant fraction of the real scattering. The method is deterministic and can be fast compared with existing multi-trial dual-space methods for SAD substructure determination.

Substructure determination using phase-retrieval techniques

Thus far, the application of phase-retrieval methods in crystallography has mainly been aimed at variants of charge flipping or structure-factor flipping. In this work, the relaxed averaged alternating reflections (RAAR) algorithm is applied to determine anomalously scattering substructures from single-wavelength anomalous diffraction (SAD) data of macromolecules. The algorithm has been implemented in a new program, PRASA, and has been shown to significantly outperform charge flipping in determining anomalously scattering substructures on a test sample of 169 SAD data sets with resolutions up to 3.88 Å.

An introduction to experimental phasing of macromolecules illustrated by SHELX; new autotracing features

For the purpose of this article, experimental phasing is understood to mean the determination of macromolecular structures by exploiting small intensity differences of Friedel opposites and possibly of reflections measured at different wavelengths or for heavy-atom derivatives, without the use of specific structural models. The SHELX programs provide a robust and efficient route for routine structure solution by the SAD, MAD and related methods, but involve a number of simplifying assumptions that may limit their applicability in borderline cases. The substructure atoms (i.e. those with significant anomalous scattering) are first located by direct methods, and the experimental data are then used to estimate phase shifts that are added to the substructure phases to obtain starting phases for the native reflections. These are then improved by density modification and, if the resolution of the data and the type of structure permit, polyalanine tracing. A number of extensions to the tracing algorithm are discussed; these are designed to improve its performance at low resolution. Given native data to 2.5 Å resolution or better, a correlation coefficient greater than 25% between the structure factors calculated from such a trace and the native data is usually a good indication that the structure has been solved.

Experimental phasing of serial femtosecond crystallography data

A synopsis of and prospects for de novo phasing using diffraction data collected at X-ray free-electron lasers are given.

The hidden treasure in your data: phasing with unexpected weak anomalous scatterers from routine data sets

Single-wavelength anomalous dispersion (SAD) utilizing anomalous signal from native S atoms, or other atoms with Z ≤ 20, generally requires highly redundant data collected using relatively long-wavelength X-rays. Here, the results from two proteins are presented where the anomalous signal from serendipitously acquired surface-bound Ca atoms with an anomalous data multiplicity of around 10 was utilized to drive de novo structure determination. In both cases, the Ca atoms were acquired from the crystallization solution, and the data-collection strategy was not optimized to exploit the anomalous signal from these scatterers. The X-ray data were collected at 0.98 Å wavelength in one case and at 1.74 Å in the other (the wavelength was optimized for sulfur, but the anomalous signal from calcium was exploited for structure solution). Similarly, using a test case, it is shown that data collected at ∼1.0 Å wavelength, where the f'' value for sulfur is 0.28 e, are sufficient for structure determination using intrinsic S atoms from a strongly diffracting crystal. Interestingly, it was also observed that SHELXD was capable of generating a substructure solution from high-exposure data with a completeness of 70% for low-resolution reflections extending to 3.5 Å resolution with relatively low anomalous multiplicity. Considering the fact that many crystallization conditions contain anomalous scatterers such as Cl, Ca, Mn etc., checking for the presence of fortuitous anomalous signal in data from well diffracting crystals could prove useful in either determining the structure de novo or in accurately assigning surface-bound atoms.

Experimental Phasing: Substructure Solution and Density Modification as Implemented in SHELX

This chapter describes experimental phasing methods as implemented in SHELX. After introducing fundamental concepts underlying all experimental phasing approaches, the methods used by SHELXC/D/E are described in greater detail, such as dual-space direct methods, Patterson seeding and density modification with the sphere of influence algorithm. Intensity differences from data for experimental phasing can also be used for the generation and usage of difference maps with ANODE for validation and phasing purposes. A short section describes how molecular replacement can be combined with experimental phasing methods. The second half covers practical challenges, such as prerequisites for successful experimental phasing, evaluation of potential solutions, and what to do if substructure search or density modification fails. It is also shown how auto-tracing in SHELXE can improve automation and how it ties in with automatic model building after phasing.

Selenourea: a convenient phasing vehicle for macromolecular X-ray crystal structures

Majority of novel X-ray crystal structures of proteins are currently solved using the anomalous diffraction signal provided by selenium after incorporation of selenomethionine instead of natural methionine by genetic engineering methods. However, selenium can be inserted into protein crystals in the form of selenourea (SeC(NH₂)₂), by adding the crystalline powder of selenourea into mother liquor or cryo-solution with native crystals, in analogy to the classic procedure of heavy-atom derivatization. Selenourea is able to bind to reactive groups at the surface of macromolecules primarily through hydrogen bonds, where the selenium atom may serve as acceptor and amide groups as donors. Selenourea has different chemical properties than heavy-atom reagents and halide ions and provides a convenient way of phasing crystal structures of macromolecules.

Efficient merging of data from multiple samples for determination of anomalous substructure

The benefits of using local scaling and optimization of anomalous signal (as implemented in PHENIX) for merging data sets from many crystals for determination of the substructure for weak anomalous scatterers are examined.

Merging of data from multiple crystals has proven to be useful for determination of the anomalously scattering atomic substructure for crystals with weak anomalous scatterers (e.g. S and P) and/or poor diffraction. Strategies for merging data from many samples, which require assessment of sample isomorphism, rely on metrics of variability in unit-cell parameters, anomalous signal correlation and overall data similarity. Local scaling, anomalous signal optimization and data-set weighting, implemented in phenix.scale_and_merge, provide an efficient protocol for merging data from many samples. The protein NS1 was used in a series of trials with data collected from 28 samples for phasing by single-wavelength anomalous diffraction of the native S atoms. The local-scaling, anomalous-optimization protocol produced merged data sets with higher anomalous signal quality indicators than did standard global-scaling protocols. The local-scaled data were also more successful in substructure determination. Merged data quality was assessed for data sets where the multiplicity was reduced in either of two ways: by excluding data from individual crystals (to reduce errors owing to non-isomorphism) or by excluding the last-recorded segments of data from each crystal (to minimize the effects of radiation damage). The anomalous signal was equivalent at equivalent multiplicity for the two procedures, and structure-determination success correlated with anomalous signal metrics. The quality of the anomalous signal was strongly correlated with data multiplicity over a range of 12-fold to 150-fold multiplicity. For the NS1 data, the local-scaling and anomalous-optimization protocol handled sample non-isomorphism and radiation-induced decay equally well.

Can I solve my structure by SAD phasing? Planning an experiment, scaling data and evaluating the useful anomalous correlation and anomalous signal

A key challenge in the SAD phasing method is solving a structure when the anomalous signal-to-noise ratio is low. Here, algorithms and tools for evaluating and optimizing the useful anomalous correlation and the anomalous signal in a SAD experiment are described. A simple theoretical framework [Terwilliger et al. (2016), Acta Cryst. D72, 346-358] is used to develop methods for planning a SAD experiment, scaling SAD data sets and estimating the useful anomalous correlation and anomalous signal in a SAD data set. The phenix.plan_sad_experiment tool uses a database of solved and unsolved SAD data sets and the expected characteristics of a SAD data set to estimate the probability that the anomalous substructure will be found in the SAD experiment and the expected map quality that would be obtained if the substructure were found. The phenix.scale_and_merge tool scales unmerged SAD data from one or more crystals using local scaling and optimizes the anomalous signal by identifying the systematic differences among data sets, and the phenix.anomalous_signal tool estimates the useful anomalous correlation and anomalous signal after collecting SAD data and estimates the probability that the data set can be solved and the likely figure of merit of phasing.