Ice in biomolecular cryocrystallography.

Diffraction data acquired from cryocooled protein crystals often include diffraction from ice. Analysis of ice diffraction from crystals of three proteins shows that the ice formed within solvent cavities during rapid cooling is comprised of a stacking-disordered mixture of hexagonal and cubic planes, with the cubic plane fraction increasing with increasing cryoprotectant concentration and increasing cooling rate. Building on the work of Thorn and coworkers [Thorn et al. (2017), Acta Cryst. D73, 729-727], a revised metric is defined for detecting ice from deposited protein structure-factor data, and this metric is validated using full-frame diffraction data from the Integrated Resource for Reproducibility in Macromolecular Crystallography. Using this revised metric and improved algorithms, an analysis of structure-factor data from a random sample of 89 827 PDB entries collected at cryogenic temperatures indicates that roughly 16% show evidence of ice contamination, and that this fraction increases with increasing solvent content and maximum solvent-cavity size. By examining the ice diffraction-peak positions at which structure-factor perturbations are observed, it is found that roughly 25% of crystals exhibit ice with primarily hexagonal character, indicating that inadequate cooling rates and/or cryoprotectant concentrations were used, while the remaining 75% show ice with a stacking-disordered or cubic character.

Ice formation and solvent nanoconfinement in protein crystals.

Ice formation within protein crystals is a major obstacle to the cryocrystallographic study of protein structure, and has limited studies of how the structural ensemble of a protein evolves with temperature in the biophysically interesting range from ∼260 K to the protein-solvent glass transition near 200 K. Using protein crystals with solvent cavities as large as ∼70 Å, time-resolved X-ray diffraction was used to study the response of protein and internal solvent during rapid cooling. Solvent nanoconfinement suppresses freezing temperatures and ice-nucleation rates so that ice-free, low-mosaicity diffraction data can be reliably collected down to 200 K without the use of cryoprotectants. Hexagonal ice (Ih) forms in external solvent, but internal crystal solvent forms stacking-disordered ice (Isd) with a near-random stacking of cubic and hexagonal planes. Analysis of powder diffraction from internal ice and single-crystal diffraction from the host protein structure shows that the maximum crystallizable solvent fraction decreases with decreasing crystal solvent-cavity size, and that an ∼6 Šthick layer of solvent adjacent to the protein surface cannot crystallize. These results establish protein crystals as excellent model systems for the study of nanoconfined solvent. By combining fast cooling, intense X-ray beams and fast X-ray detectors, complete structural data sets for high-value targets, including membrane proteins and large complexes, may be collected at ∼220-240 K that have much lower mosaicities and comparable B factors, and that may allow more confident identification of ligand binding than in current cryocrystallographic practice.

Solvent flows, conformation changes and lattice reordering in a cold protein crystal.

When protein crystals are abruptly cooled, the unit-cell, protein and solvent-cavity volumes all contract, but the volume of bulk-like internal solvent may expand. Outflow of this solvent from the unit cell and its accumulation in defective interior crystal regions has been suggested as one cause of the large increase in crystal mosaicity on cooling. It is shown that when apoferritin crystals are abruptly cooled to temperatures between 220 and 260 K, the unit cell contracts, solvent is pushed out and the mosaicity grows. On temperature-dependent timescales of 10 to 200 s, the unit-cell and solvent-cavity volume then expand, solvent flows back in, and the mosaicity and B factor both drop. Expansion and reordering at fixed low temperature are associated with small-amplitude but large-scale changes in the conformation and packing of apoferritin. These results demonstrate that increases in mosaicity on cooling arise due to solvent flows out of or into the unit cell and to incomplete, arrested relaxation of protein conformation. They indicate a critical role for time in variable-temperature crystallographic studies, and the feasibility of probing interactions and cooperative conformational changes that underlie cold denaturation in the presence of liquid solvent at temperatures down to ∼200 K.

Resolution and dose dependence of radiation damage in biomolecular systems.

The local Fourier-space relation between diffracted intensity I, diffraction wavevector q and dose D, , is key to probing and understanding radiation damage by X-rays and energetic particles in both diffraction and imaging experiments. The models used in protein crystallography for the last 50 years provide good fits to experimental I(q) versus nominal dose data, but have unclear physical significance. More recently, a fit to diffraction and imaging experiments suggested that the maximum tolerable dose varies as q -1 or linearly with resolution. Here, it is shown that crystallographic data have been strongly perturbed by the effects of spatially nonuniform crystal irradiation and diffraction during data collection. Reanalysis shows that these data are consistent with a purely exponential local dose dependence, = I 0(q)exp[-D/D e(q)], where D e(q) ∝ q α with α ≃ 1.7. A physics-based model for radiation damage, in which damage events occurring at random locations within a sample each cause energy deposition and blurring of the electron density within a small volume, predicts this exponential variation with dose for all q values and a decay exponent α ≃ 2 in two and three dimensions, roughly consistent with both diffraction and imaging experiments over more than two orders of magnitude in resolution. The B-factor model used to account for radiation damage in crystallographic scaling programs is consistent with α = 2, but may not accurately capture the dose dependencies of structure factors under typical nonuniform illumination conditions. The strong q dependence of radiation-induced diffraction decays implies that the previously proposed 20-30 MGy dose limit for protein crystallography should be replaced by a resolution-dependent dose limit that, for atomic resolution data sets, will be much smaller. The results suggest that the physics underlying basic experimental trends in radiation damage at T ≃ 100 K is straightforward and universal. Deviations of the local I(q, D) from strictly exponential behavior may provide mechanistic insights, especially into the radiation-damage processes responsible for the greatly increased radiation sensitivity observed at T ≃ 300 K.

Effects of protein-crystal hydration and temperature on side-chain conformational heterogeneity in monoclinic lysozyme crystals.

The modulation of main-chain and side-chain conformational heterogeneity and solvent structure in monoclinic lysozyme crystals by dehydration (related to water activity) and temperature is examined. Decreasing the relative humidity (from 99 to 11%) and decreasing the temperature both lead to contraction of the unit cell, to an increased area of crystal contacts and to remodeling of primarily contact and solvent-exposed residues. Both lead to the depopulation of some minor side-chain conformers and to the generation of new conformations. Side-chain modifications and main-chain r.m.s.d.s associated with cooling from 298 to 100 K depend on relative humidity and are minimized at 85% relative humidity (r.h.). Dehydration from 99 to 93% r.h. and cooling from 298 to 100 K result in a comparable number of remodeled residues, with dehydration-induced remodeling somewhat more likely to arise from contact interactions. When scaled to equivalent temperatures based on unit-cell contraction, the evolution of side-chain order parameters with dehydration shows generally similar features to those observed on cooling to T = 100 K. These results illuminate the qualitative and quantitative similarities between structural perturbations induced by modest dehydration, which routinely occurs in samples prepared for 298 and 100 K data collection, and cryocooling. Differences between these perturbations in terms of energy landscapes and occupancies, and implications for variable-temperature crystallography between 180 and 298 K, are discussed. It is also noted that remodeling of a key lysozyme active-site residue by dehydration, which is associated with a radical decrease in the enzymatic activity of lysozyme powder, arises due to a steric clash with the residue of a symmetry mate.

Lifetimes and spatio-temporal response of protein crystals in intense X-ray microbeams.

Serial synchrotron-based crystallography using intense microfocused X-ray beams, fast-framing detectors and protein microcrystals held at 300 K promises to expand the range of accessible structural targets and to increase overall structure-pipeline throughputs. To explore the nature and consequences of X-ray radiation damage under microbeam illumination, the time-, dose- and temperature-dependent evolution of crystal diffraction have been measured with maximum dose rates of 50 MGy s-1. At all temperatures and dose rates, the integrated diffraction intensity for a fixed crystal orientation shows non-exponential decays with dose. Non-exponential decays are a consequence of non-uniform illumination and the resulting spatial evolution of diffracted intensity within the illuminated crystal volume. To quantify radiation-damage lifetimes and the damage state of diffracting crystal regions, a revised diffraction-weighted dose (DWD) is defined and it is shown that for Gaussian beams the DWD becomes nearly independent of actual dose at large doses. An apparent delayed onset of radiation damage seen in some intensity-dose curves is in fact a consequence of damage. Intensity fluctuations at high dose rates may arise from the impulsive release of gaseous damage products. Accounting for these effects, data collection at the highest dose rates increases crystal radiation lifetimes near 300 K (but not at 100 K) by a factor of ∼1.5-2 compared with those observed at conventional dose rates. Improved quantification and modeling of the complex spatio-temporal evolution of protein microcrystal diffraction in intense microbeams will enable more efficient data collection, and will be essential in improving the accuracy of structure factors and structural models.

Active-site protein dynamics and solvent accessibility in native Achromobacter cycloclastes copper nitrite reductase.

Microbial nitrite reductases are denitrifying enzymes that are a major component of the global nitrogen cycle. Multiple structures measured from one crystal (MSOX data) of copper nitrite reductase at 240 K, together with molecular-dynamics simulations, have revealed protein dynamics at the type 2 copper site that are significant for its catalytic properties and for the entry and exit of solvent or ligands to and from the active site. Molecular-dynamics simulations were performed using different protonation states of the key catalytic residues (AspCAT and HisCAT) involved in the nitrite-reduction mechanism of this enzyme. Taken together, the crystal structures and simulations show that the AspCAT protonation state strongly influences the active-site solvent accessibility, while the dynamics of the active-site 'capping residue' (IleCAT), a determinant of ligand binding, are influenced both by temperature and by the protonation state of AspCAT. A previously unobserved conformation of IleCAT is seen in the elevated temperature series compared with 100 K structures. DFT calculations also show that the loss of a bound water ligand at the active site during the MSOX series is consistent with reduction of the type 2 Cu atom.

Thermal contraction of aqueous glycerol and ethylene glycol solutions for optimized protein-crystal cryoprotection.

The thermal contraction of aqueous cryoprotectant solutions on cooling to cryogenic temperatures is of practical importance in protein cryocrystallography and in biological cryopreservation. In the former case, differential contraction on cooling of protein molecules and their lattice relative to that of the internal and surrounding solvent may lead to crystal damage and the degradation of crystal diffraction properties. Here, the amorphous phase densities of aqueous solutions of glycerol and ethylene glycol at T = 77 K have been determined. Densities with accuracies of <0.5% to concentrations as low as 30%(w/v) were determined by rapidly cooling drops with volumes as small as 70 pl, assessing their optical clarity and measuring their buoyancy in liquid nitrogen-argon solutions. The use of these densities in contraction matching of internal solvent to the available solvent spaces is complicated by several factors, most notably the exclusion of cryoprotectants from protein hydration shells and the expected deviation of the contraction behavior of hydration water from bulk water. The present methods and results will assist in developing rational approaches to cryoprotection and an understanding of solvent behavior in protein crystals.