Supercooled liquid-like solvent in trypsin crystals: implications for crystal annealing and temperature-controlled X-ray radiation damage studies

Refinement of multiconformer ensemble models from multi-temperature X-ray diffraction data

Conformational ensembles underlie all protein functions. Thus, acquiring atomic-level ensemble models that accurately represent conformational heterogeneity is vital to deepen our understanding of how proteins work. Modeling ensemble information from X-ray diffraction data has been challenging, as traditional cryo-crystallography restricts conformational variability while minimizing radiation damage. Recent advances have enabled the collection of high quality diffraction data at ambient temperatures, revealing innate conformational heterogeneity and temperature-driven changes. Here, we used diffraction datasets for Proteinase K collected at temperatures ranging from 313 to 363 K to provide a tutorial for the refinement of multiconformer ensemble models. Integrating automated sampling and refinement tools with manual adjustments, we obtained multiconformer models that describe alternative backbone and sidechain conformations, their relative occupancies, and interconnections between conformers. Our models revealed extensive and diverse conformational changes across temperature, including increased bound peptide ligand occupancies, different Ca2+ binding site configurations and altered rotameric distributions. These insights emphasize the value and need for multiconformer model refinement to extract ensemble information from diffraction data and to understand ensemble-function relationships.

Refinement of Multiconformer Ensemble Models from Multi-temperature X-ray Diffraction Data

Conformational ensembles underlie all protein functions. Thus, acquiring atomic-level ensemble models that accurately represent conformational heterogeneity is vital to deepen our understanding of how proteins work. Modeling ensemble information from X-ray diffraction data has been challenging, as traditional cryo-crystallography restricts conformational variability while minimizing radiation damage. Recent advances have enabled the collection of high quality diffraction data at ambient temperatures, revealing innate conformational heterogeneity and temperature-driven changes. Here, we used diffraction datasets for Proteinase K collected at temperatures ranging from 313 to 363K to provide a tutorial for the refinement of multiconformer ensemble models. Integrating automated sampling and refinement tools with manual adjustments, we obtained multiconformer models that describe alternative backbone and sidechain conformations, their relative occupancies, and interconnections between conformers. Our models revealed extensive and diverse conformational changes across temperature, including increased bound peptide ligand occupancies, different Ca2+ binding site configurations and altered rotameric distributions. These insights emphasize the value and need for multiconformer model refinement to extract ensemble information from diffraction data and to understand ensemble-function relationships.

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.

The impact of cryosolution thermal contraction on proteins and protein crystals: volumes, conformation and order

Cryocooling of macromolecular crystals is commonly employed to limit radiation damage during X-ray diffraction data collection. However, cooling itself affects macromolecular conformation and often damages crystals via poorly understood processes. Here, the effects of cryosolution thermal contraction on macromolecular conformation and crystal order in crystals ranging from 32 to 67% solvent content are systematically investigated. It is found that the solution thermal contraction affects macromolecule configurations and volumes, unit-cell volumes, crystal packing and crystal order. The effects occur through not only thermal contraction, but also pressure caused by the mismatched contraction of cryosolvent and pores. Higher solvent-content crystals are more affected. In some cases the solvent contraction can be adjusted to reduce mosaicity and increase the strength of diffraction. Ice formation in some crystals is found to cause damage via a reduction in unit-cell volume, which is interpreted through solvent transport out of unit cells during cooling. The results point to more deductive approaches to cryoprotection optimization by adjusting the cryosolution composition to reduce thermal contraction-induced stresses in the crystal with cooling.

Crystallization of recombinant green mamba ρ-Da1a toxin during a lyophilization procedure and its structure determination

ρ-Da1a toxin from eastern green mamba (Dendroaspis angusticeps) venom is a polypeptide of 65 amino acids with a strong affinity for the G-protein-coupled α(1A)-adrenoceptor. This neurotoxin has been crystallized from resolubilized lyophilized powder, but the best crystals grew spontaneously during lyophilization. The crystals belonged to the trigonal space group P3(1)21, with unit-cell parameters a = b = 37.37, c = 66.05 Å, and diffracted to 1.95 Å resolution. The structure solved by molecular replacement showed strong similarities to green mamba muscarinic toxins.

Global radiation damage: temperature dependence, time dependence and how to outrun it

A series of studies that provide a consistent and illuminating picture of global radiation damage to protein crystals, especially at temperatures above ∼200 K, are described. The radiation sensitivity shows a transition near 200 K, above which it appears to be limited by solvent-coupled diffusive processes. Consistent with this interpretation, a component of global damage proceeds on timescales of several minutes at 180 K, decreasing to seconds near room temperature. As a result, data collection times of order 1 s allow up to half of global damage to be outrun at 260 K. Much larger damage reductions near room temperature should be feasible using larger dose rates delivered using microfocused beams, enabling a significant expansion of structural studies of proteins under more nearly native conditions.

Spatial distribution of radiation damage to crystalline proteins at 25-300 K

The spatial distribution of radiation damage (assayed by increases in atomic B factors) to thaumatin and urease crystals at temperatures ranging from 25 to 300 K is reported. The nature of the damage changes dramatically at approximately 180 K. Above this temperature the role of solvent diffusion is apparent in thaumatin crystals, as solvent-exposed turns and loops are especially sensitive. In urease, a flap covering the active site is the most sensitive part of the molecule and nearby loops show enhanced sensitivity. Below 180 K sensitivity is correlated with poor local packing, especially in thaumatin. At all temperatures, the component of the damage that is spatially uniform within the unit cell accounts for more than half of the total increase in the atomic B factors and correlates with changes in mosaicity. This component may arise from lattice-level, rather than local, disorder. The effects of primary structure on radiation sensitivity are small compared with those of tertiary structure, local packing, solvent accessibility and crystal contacts.

Dark progression reveals slow timescales for radiation damage between T = 180 and 240 K

Can radiation damage to protein crystals be `outrun' by collecting a structural data set before damage is manifested? Recent experiments using ultra-intense pulses from a free-electron laser show that the answer is yes. Here, evidence is presented that significant reductions in global damage at temperatures above 200 K may be possible using conventional X-ray sources and current or soon-to-be available detectors. Specifically, `dark progression' (an increase in damage with time after the X-rays have been turned off) was observed at temperatures between 180 and 240 K and on timescales from 200 to 1200 s. This allowed estimation of the temperature-dependent timescale for damage. The rate of dark progression is consistent with an Arrhenius law with an activation energy of 14 kJ mol(-1). This is comparable to the activation energy for the solvent-coupled diffusive damage processes responsible for the rapid increase in radiation sensitivity as crystals are warmed above the glass transition near 200 K. Analysis suggests that at T = 300 K data-collection times of the order of 1 s (and longer at lower temperatures) may allow significant reductions in global radiation damage, facilitating structure solution on crystals with liquid solvent. No dark progression was observed below T = 180 K, indicating that no important damage process is slowed through this timescale window in this temperature range.

Glass transition in thaumatin crystals revealed through temperature-dependent radiation-sensitivity measurements

The temperature-dependence of radiation damage to thaumatin crystals between T = 300 and 100 K is reported. The amount of damage for a given dose decreases sharply as the temperature decreases from 300 to 220 K and then decreases more gradually on further cooling below the protein-solvent glass transition. Two regimes of temperature-activated behavior were observed. At temperatures above ∼200 K the activation energy of 18.0 kJ mol(-1) indicates that radiation damage is dominated by diffusive motions in the protein and solvent. At temperatures below ∼200 K the activation energy is only 1.00 kJ mol(-1), which is of the order of the thermal energy. Similar activation energies describe the temperature-dependence of radiation damage to a variety of solvent-free small-molecule organic crystals over the temperature range T = 300-80 K. It is suggested that radiation damage in this regime is vibrationally assisted and that the freezing-out of amino-acid scale vibrations contributes to the very weak temperature-dependence of radiation damage below ∼80 K. Analysis using the radiation-damage model of Blake and Phillips [Blake & Phillips (1962), Biological Effects of Ionizing Radiation at the Molecular Level, pp. 183-191] indicates that large-scale conformational and molecular motions are frozen out below T = 200 K but become increasingly prevalent and make an increasing contribution to damage at higher temperatures. Possible alternative mechanisms for radiation damage involving the formation of hydrogen-gas bubbles are discussed and discounted. These results have implications for mechanistic studies of proteins and for studies of the protein glass transition. They also suggest that data collection at T ≃ 220 K may provide a viable alternative for structure determination when cooling-induced disorder at T = 100 is excessive.

Radiation damage in macromolecular crystallography: what is it and why should we care?

Radiation damage inflicted during diffraction data collection in macromolecular crystallography has re-emerged in the last decade as a major experimental and computational challenge, as even for crystals held at 100 K it can result in severe data-quality degradation and the appearance in solved structures of artefacts which affect biological interpretations. Here, the observable symptoms and basic physical processes involved in radiation damage are described and the concept of absorbed dose as the basic metric against which to monitor the experimentally observed changes is outlined. Investigations into radiation damage in macromolecular crystallography are ongoing and the number of studies is rapidly increasing. The current literature on the subject is compiled as a resource for the interested researcher.

Progress in rational methods of cryoprotection in macromolecular crystallography

Cryogenic cooling of macromolecular crystals is commonly used for X-ray data collection both to reduce crystal damage from radiation and to gather functional information by cryogenically trapping intermediates. However, the cooling process can damage the crystals. Limiting cooling-induced crystal damage often requires cryoprotection strategies, which can involve substantial screening of solution conditions and cooling protocols. Here, recent developments directed towards rational methods for cryoprotection are described. Crystal damage is described in the context of the temperature response of the crystal as a thermodynamic system. As such, the internal and external parts of the crystal typically have different cryoprotection requirements. A key physical parameter, the thermal contraction, of 26 different cryoprotective solutions was measured between 294 and 72 K. The range of contractions was 2-13%, with the more polar cryosolutions contracting less. The potential uses of these results in the development of cryocooling conditions, as well as recent developments in determining minimum cryosolution soaking times, are discussed.

Temperature-dependent macromolecular X-ray crystallography

X-ray crystallography provides structural details of biological macromolecules. Whereas routine data are collected close to 100 K in order to mitigate radiation damage, more exotic temperature-controlled experiments in a broader temperature range from 15 K to room temperature can provide both dynamical and structural insights. Here, the dynamical behaviour of crystalline macromolecules and their surrounding solvent as a function of cryo-temperature is reviewed. Experimental strategies of kinetic crystallography are discussed that have allowed the generation and trapping of macromolecular intermediate states by combining reaction initiation in the crystalline state with appropriate temperature profiles. A particular focus is on recruiting X-ray-induced changes for reaction initiation, thus unveiling useful aspects of radiation damage, which otherwise has to be minimized in macromolecular crystallography.

Slow cooling of protein crystals

Cryoprotectant-free thaumatin crystals have been cooled from 300 to 100 K at a rate of 0.1 K s(-1) - 10(3)-10(4) times slower than in conventional flash cooling - while continuously collecting X-ray diffraction data, so as to follow the evolution of protein lattice and solvent properties during cooling. Diffraction patterns show no evidence of crystalline ice at any temperature. This indicates that the lattice of protein molecules is itself an excellent cryoprotectant, and with sodium potassium tartrate incorporated from the 1.5 M mother liquor ice nucleation rates are at least as low as in a 70% glycerol solution. Crystal quality during slow cooling remains high, with an average mosaicity at 100 K of 0.2 degrees . Most of the mosaicity increase occurs above approximately 200 K, where the solvent is still liquid, and is concurrent with an anisotropic contraction of the unit cell. Near 180 K a crossover to solid-like solvent behavior occurs, and on further cooling there is no additional degradation of crystal order. The variation of B factor with temperature shows clear evidence of a protein dynamical transition near 210 K, and at lower temperatures the slope dB/dT is a factor of 3-6 smaller than has been reported for any other protein. These results establish the feasibility of fully temperature controlled studies of protein structure and dynamics between 300 and 100 K.

Evidence for liquid water during the high-density to low-density amorphous ice transition

Polymorphism of water has been extensively studied, but controversy still exists over the phase transition between high-density amorphous (HDA) and low-density amorphous (LDA) ice. We report the phase behavior of HDA ice inside high-pressure cryocooled protein crystals. Using X-ray diffraction, we demonstrate that the intermediate states in the temperature range from 80 to 170 K can be reconstructed as a linear combination of HDA and LDA ice, suggesting a first-order transition. We found evidence for a liquid state of water during the ice transition based on the protein crystallographic data. These observations open the possibility that the HDA ice induced by high-pressure cryocooling is a genuine glassy form of high-density liquid.

Shoot-and-Trap: use of specific x-ray damage to study structural protein dynamics by temperature-controlled cryo-crystallography

Although x-ray crystallography is the most widely used method for macromolecular structure determination, it does not provide dynamical information, and either experimental tricks or complementary experiments must be used to overcome the inherently static nature of crystallographic structures. Here we used specific x-ray damage during temperature-controlled crystallographic experiments at a third-generation synchrotron source to trigger and monitor (Shoot-and-Trap) structural changes putatively involved in an enzymatic reaction. In particular, a nonhydrolyzable substrate analogue of acetylcholinesterase, the "off-switch" at cholinergic synapses, was radiocleaved within the buried enzymatic active site. Subsequent product clearance, observed at 150 K but not at 100 K, indicated exit from the active site possibly via a "backdoor." The simple strategy described here is, in principle, applicable to any enzyme whose structure in complex with a substrate analogue is available and, therefore, could serve as a standard procedure in kinetic crystallography studies.

Structural characterization of water and ice in mesoporous SBA-15 silicas: II. The 'almost-filled' case for 86 Å pore diameter

Neutron diffraction measurements for D(2)O in SBA-15 silica of pore diameter 86 Å have been made in a temperature range from 300 to 100 K. The pore-filling factor for the liquid phase is 0.95, resulting in an 'almost-filled' sample. The nucleation and transformation of the ice phase were determined for cooling and warming cycles at two different rates. The primary nucleation event at 258 K leads to a defective form of ice-I with predominantly cubic ice features. For temperatures below the main nucleation event, the data indicate the formation of an interfacial layer of disordered water/ice that varies with temperature and is reversible. The main diffraction peak for the water phase shows similar features to those observed in earlier studies, indicating enhanced hydrogen bonding and network correlations for the confined phase as the temperature is decreased. A detailed profile analysis of the triplet peak is presented in the accompanying paper (Seyed-Yazdi et al 2008 J. Phys.: Condens. Matter 20 205108).

Hyperquenching for protein cryocrystallography

When samples having volumes characteristic of protein crystals are plunge cooled in liquid nitrogen or propane, most cooling occurs in the cold gas layer above the liquid. By removing this cold gas layer, cooling rates for small samples and modest plunge velocities are increased to 1.5 × 10(4) K s(-1), with increases of a factor of 100 over current best practice possible with 10 μm samples. Glycerol concentrations required to eliminate water crystallization in protein-free aqueous mixtures drop from ∼28% w/v to as low as 6% w/v. These results will allow many crystals to go from crystallization tray to liquid cryogen to X-ray beam without cryoprotectants. By reducing or eliminating the need for cryoprotectants in growth solutions, they may also simplify the search for crystallization conditions and for optimal screens. The results presented here resolve many puzzles, such as why plunge cooling in liquid nitrogen or propane has, until now, not yielded significantly better diffraction quality than gas-stream cooling.

Liquid-like water confined in stacks of biological membranes at 200 k and its relation to protein dynamics

Confined water is of considerable current interest owing to its biophysical importance and relevance to cryopreservation. It can be studied in its amorphous or supercooled state in the "no-man's land", i.e., in the temperature range between 150 and 235 K, in which bulk water is always crystalline. Amorphous deuterium oxide (D(2)O) was obtained in the intermembrane spaces of a stack of purple membranes from Halobacterium salinarum by flash cooling to 77 K. Neutron diffraction showed that upon heating to 200 K the intermembrane water space decreased sharply with an associated strengthening of ice diffraction, indicating that water beyond the first membrane hydration layer flowed out of the intermembrane space to form crystalline ice. It was concluded that the confined water undergoes a glass transition at or below 200 K to adopt an ultraviscous liquid state from which it crystallizes to form ice as soon as it finds itself in an unconfined, bulk-water environment. Our results provide model-free evidence for translational diffusion of confined water in the no-man's land. Potential effects of the confined-water glass transition on nanosecond membrane dynamics were investigated by incoherent elastic neutron scattering experiments. These revealed no differences between flash-cooled and slow-cooled samples (in the latter, the intermembrane space at temperatures <250 K is occupied only by the first membrane hydration layers), with dynamical transitions at 150 and 260 K, but not at 200 K, suggesting that nanosecond membrane dynamics are not sensitive to the state of the water beyond the first hydration shell at cryotemperatures.