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

MatchMaps: non-isomorphous difference maps for X-ray crystallography

Conformational change mediates the biological functions of macromolecules. Crystallographic measurements can map these changes with extraordinary sensitivity as a function of mutations, ligands and time. A popular method for detecting structural differences between crystallographic data sets is the isomorphous difference map. These maps combine the phases of a chosen reference state with the observed changes in structure factor amplitudes to yield a map of changes in electron density. Such maps are much more sensitive to conformational change than structure refinement is, and are unbiased in the sense that observed differences do not depend on refinement of the perturbed state. However, even modest changes in unit-cell properties can render isomorphous difference maps useless. This is unnecessary. Described here is a generalized procedure for calculating observed difference maps that retains the high sensitivity to conformational change and avoids structure refinement of the perturbed state. This procedure is implemented in an open-source Python package, MatchMaps, that can be run in any software environment supporting PHENIX [Liebschner et al. (2019). Acta Cryst. D75, 861-877] and CCP4 [Agirre et al. (2023). Acta Cryst. D79, 449-461]. Worked examples show that MatchMaps 'rescues' observed difference electron-density maps for poorly isomorphous crystals, corrects artifacts in nominally isomorphous difference maps, and extends to detecting differences across copies within the asymmetric unit or across altogether different crystal forms.

Combining temperature perturbations with X-ray crystallography to study dynamic macromolecules: A thorough discussion of experimental methods

Temperature is an important state variable that governs the behavior of microscopic systems, yet crystallographers rarely exploit temperature changes to study the structure and dynamics of biological macromolecules. In fact, approximately 90% of crystal structures in the Protein Data Bank were determined under cryogenic conditions, because sample cryocooling makes crystals robust to X-ray radiation damage and facilitates data collection. On the other hand, cryocooling can introduce artifacts into macromolecular structures, and can suppress conformational dynamics that are critical for function. Fortunately, recent advances in X-ray detector technology, X-ray sources, and computational data processing algorithms make non-cryogenic X-ray crystallography easier and more broadly applicable than ever before. Without the reliance on cryocooling, high-resolution crystallography can be combined with various temperature perturbations to gain deep insight into the conformational landscapes of macromolecules. This Chapter reviews the historical reasons for the prevalence of cryocooling in macromolecular crystallography, and discusses its potential drawbacks. Next, the Chapter summarizes technological developments and methodologies that facilitate non-cryogenic crystallography experiments. Finally, the chapter discusses the theoretical underpinnings and practical aspects of multi-temperature and temperature-jump crystallography experiments, which are powerful tools for understanding the relationship between the structure, dynamics, and function of proteins and other biological macromolecules.

MD simulations of macromolecular crystals: Implications for the analysis of Bragg and diffuse scattering

Some of our most detailed information about structure and dynamics of macromolecules comes from X-ray-diffraction studies in crystalline environments. More than 170,000 atomic models have been deposited in the Protein Data Bank, and the number of observations (typically of intensities of Bragg diffraction peaks) is generally quite large, when compared to other experimental methods. Nevertheless, the general agreement between calculated and observed intensities is far outside the experimental precision, and the majority of scattered photons fall between the sharp Bragg peaks, and are rarely taken into account. This chapter considers how molecular dynamics simulations can be used to explore the connections between microscopic behavior in a crystalline lattice and observed scattering intensities, and point the way to new atomic models that could more faithfully recapitulate Bragg intensities and extract useful information from the diffuse scattering that lies between those peaks.

Emerging Time-Resolved X-Ray Diffraction Approaches for Protein Dynamics

Proteins guide the flows of information, energy, and matter that make life possible by accelerating transport and chemical reactions, by allosterically modulating these reactions, and by forming dynamic supramolecular assemblies. In these roles, conformational change underlies functional transitions. Time-resolved X-ray diffraction methods characterize these transitions either by directly triggering sequences of functionally important motions or, more broadly, by capturing the motions of which proteins are capable. To date, most successful have been experiments in which conformational change is triggered in light-dependent proteins. In this review, I emphasize emerging techniques that probe the dynamic basis of function in proteins lacking natively light-dependent transitions and speculate about extensions and further possibilities. In addition, I review how the weaker and more distributed signals in these data push the limits of the capabilities of analytical methods. Taken together, these new methods are beginning to establish a powerful paradigm for the study of the physics of protein function.

Determining biomolecular structures near room temperature using X-ray crystallography: concepts, methods and future optimization

For roughly two decades, cryocrystallography has been the overwhelmingly dominant method for determining high-resolution biomolecular structures. Competition from single-particle cryo-electron microscopy and micro-electron diffraction, increased interest in functionally relevant information that may be missing or corrupted in structures determined at cryogenic temperature, and interest in time-resolved studies of the biomolecular response to chemical and optical stimuli have driven renewed interest in data collection at room temperature and, more generally, at temperatures from the protein-solvent glass transition near 200 K to ∼350 K. Fischer has recently reviewed practical methods for room-temperature data collection and analysis [Fischer (2021), Q. Rev. Biophys. 54, e1]. Here, the key advantages and physical principles of, and methods for, crystallographic data collection at noncryogenic temperatures and some factors relevant to interpreting the resulting data are discussed. For room-temperature data collection to realize its potential within the structural biology toolkit, streamlined and standardized methods for delivering crystals prepared in the home laboratory to the synchrotron and for automated handling and data collection, similar to those for cryocrystallography, should be implemented.

The temperature-dependent conformational ensemble of SARS-CoV-2 main protease (Mpro)

The COVID-19 pandemic, instigated by the SARS-CoV-2 coronavirus, continues to plague the globe. The SARS-CoV-2 main protease, or Mpro, is a promising target for the development of novel antiviral therapeutics. Previous X-ray crystal structures of Mpro were obtained at cryogenic tem-per-ature or room tem-per-ature only. Here we report a series of high-resolution crystal structures of unliganded Mpro across multiple tem-per-atures from cryogenic to physiological, and another at high humidity. We inter-rogate these data sets with parsimonious multiconformer models, multi-copy ensemble models, and isomorphous difference density maps. Our analysis reveals a perturbation-dependent conformational landscape for Mpro, including a mobile zinc ion inter-leaved between the catalytic dyad, mercurial conformational heterogeneity at various sites including a key substrate-binding loop, and a far-reaching intra-molecular network bridging the active site and dimer inter-face. Our results may inspire new strategies for antiviral drug development to aid preparation for future coronavirus pandemics.

Complementarity of neutron, XFEL and synchrotron crystallography for defining the structures of metalloenzymes at room temperature

Room-temperature macromolecular crystallography allows protein structures to be determined under close-to-physiological conditions, permits dynamic freedom in protein motions and enables time-resolved studies. In the case of metalloenzymes that are highly sensitive to radiation damage, such room-temperature experiments can present challenges, including increased rates of X-ray reduction of metal centres and site-specific radiation-damage artefacts, as well as in devising appropriate sample-delivery and data-collection methods. It can also be problematic to compare structures measured using different crystal sizes and light sources. In this study, structures of a multifunctional globin, dehaloperoxidase B (DHP-B), obtained using several methods of room-temperature crystallographic structure determination are described and compared. Here, data were measured from large single crystals and multiple microcrystals using neutrons, X-ray free-electron laser pulses, monochromatic synchrotron radiation and polychromatic (Laue) radiation light sources. These approaches span a range of 18 orders of magnitude in measurement time per diffraction pattern and four orders of magnitude in crystal volume. The first room-temperature neutron structures of DHP-B are also presented, allowing the explicit identification of the hydrogen positions. The neutron data proved to be complementary to the serial femtosecond crystallography data, with both methods providing structures free of the effects of X-ray radiation damage when compared with standard cryo-crystallography. Comparison of these room-temperature methods demonstrated the large differences in sample requirements, data-collection time and the potential for radiation damage between them. With regard to the structure and function of DHP-B, despite the results being partly limited by differences in the underlying structures, new information was gained on the protonation states of active-site residues which may guide future studies of DHP-B.

Probing ligand binding of endothiapepsin by `temperature-resolved' macromolecular crystallography

A room-temperature X-ray crystallographic method using temperature as a trigger to record movie-like structural snapshots has been developed and applied to study ligand binding and protein plasticity.

Continuous developments in cryogenic X-ray crystallography have provided most of our knowledge of 3D protein structures, which has recently been further augmented by revolutionary advances in cryoEM. However, a single structural conformation identified at cryogenic temperatures may introduce a fictitious structure as a result of cryogenic cooling artefacts, limiting the overview of inherent protein physiological dynamics, which play a critical role in the biological functions of proteins. Here, a room-temperature X-ray crystallo­graphic method using temperature as a trigger to record movie-like structural snapshots has been developed. The method has been used to show how TL00150, a 175.15 Da fragment, undergoes binding-mode changes in endothiapepsin. A surprising fragment-binding discrepancy was observed between the cryo-cooled and physiological temperature structures, and multiple binding poses and their interplay with DMSO were captured. The observations here open up new promising prospects for structure determination and interpretation at physiological temperatures with implications for structure-based drug discovery.

Evaluating the impact of X-ray damage on conformational heterogeneity in room-temperature (277 K) and cryo-cooled protein crystals

Cryo-cooling has been nearly universally adopted to mitigate X-ray damage and facilitate crystal handling in protein X-ray crystallography. However, cryo X-ray crystallographic data provide an incomplete window into the ensemble of conformations that is at the heart of protein function and energetics. Room-temperature (RT) X-ray crystallography provides accurate ensemble information, and recent developments allow conformational heterogeneity (the experimental manifestation of ensembles) to be extracted from single-crystal data. Nevertheless, high sensitivity to X-ray damage at RT raises concerns about data reliability. To systematically address this critical issue, increasingly X-ray-damaged high-resolution data sets (1.02–1.52 Å resolution) were obtained from single proteinase K, thaumatin and lysozyme crystals at RT (277 K). In each case a modest increase in conformational heterogeneity with X-ray damage was observed. Merging data with different extents of damage (as is typically carried out) had negligible effects on conformational heterogeneity until the overall diffraction intensity decayed to ∼70% of its initial value. These effects were compared with X-ray damage effects in cryo-cooled crystals by carrying out an analogous analysis of increasingly damaged proteinase K cryo data sets (0.9–1.16 Å resolution). X-ray damage-associated heterogeneity changes were found that were not observed at RT. This property renders it difficult to distinguish real from artefactual conformations and to determine the conformational response to changes in temperature. The ability to acquire reliable heterogeneity information from single crystals at RT, together with recent advances in RT data collection at accessible synchrotron beamlines, provides a strong motivation for the widespread adoption of RT X-ray crystallography to obtain conformational ensemble information.

Mechanism of Side Chain-Controlled Proton Conductivity in Bioinspired Peptidic Nanostructures

Bioinspired peptide assemblies are promising candidates for use as proton-conducting materials in electrochemical devices and other advanced technologies. Progress toward applications requires establishing foundational structure-function relationships for transport in these materials. This experimental-theoretical study sheds light on how the molecular structure and proton conduction are linked in three synthetic cyclic peptide nanotube assemblies that comprise the three canonical basic amino acids (lysine, arginine, and histidine). Experiments find an order of magnitude higher proton conductivity for lysine-containing peptide assemblies compared to histidine and arginine containing assemblies. The simulations indicate that, upon peptide assembly, the basic amino acid side chains are close enough to enable direct proton transfer. The proton transfer kinetics is determined in the simulations to be governed by the structure and flexibility of the side chains. Together, experiments and theory indicate that the proton mobility is the main determinant of proton conductivity, critical for the performance of peptide-based devices.

Native or Non-Native Protein-Protein Docking Models? Molecular Dynamics to the Rescue

Molecular docking excels at creating a plethora of potential models of protein-protein complexes. To correctly distinguish the favorable, native-like models from the remaining ones remains, however, a challenge. We assessed here if a protocol based on molecular dynamics (MD) simulations would allow distinguishing native from non-native models to complement scoring functions used in docking. To this end, the first models for 25 protein-protein complexes were generated using HADDOCK. Next, MD simulations complemented with machine learning were used to discriminate between native and non-native complexes based on a combination of metrics reporting on the stability of the initial models. Native models showed higher stability in almost all measured properties, including the key ones used for scoring in the Critical Assessment of PRedicted Interaction (CAPRI) competition, namely the positional root mean square deviations and fraction of native contacts from the initial docked model. A random forest classifier was trained, reaching a 0.85 accuracy in correctly distinguishing native from non-native complexes. Reasonably modest simulation lengths of the order of 50-100 ns are sufficient to reach this accuracy, which makes this approach applicable in practice.

MD simulations reveal the basis for dynamic assembly of Hfq-RNA complexes

The conserved protein Hfq is a key factor in the RNA-mediated control of gene expression in most known bacteria. The transient intermediates Hfq forms with RNA support intricate and robust regulatory networks. In Pseudomonas, Hfq recognizes repeats of adenine–purine–any nucleotide (ARN) in target mRNAs via its distal binding side, and together with the catabolite repression control (Crc) protein, assembles into a translation–repression complex. Earlier experiments yielded static, ensemble-averaged structures of the complex, but details of its interface dynamics and assembly pathway remained elusive. Using explicit solvent atomistic molecular dynamics simulations, we modeled the extensive dynamics of the Hfq–RNA interface and found implications for the assembly of the complex. We predict that syn/anti flips of the adenine nucleotides in each ARN repeat contribute to a dynamic recognition mechanism between the Hfq distal side and mRNA targets. We identify a previously unknown binding pocket that can accept any nucleotide and propose that it may serve as a ‘status quo’ staging point, providing nonspecific binding affinity, until Crc engages the Hfq–RNA binary complex. The dynamical components of the Hfq–RNA recognition can speed up screening of the pool of the surrounding RNAs, participate in rapid accommodation of the RNA on the protein surface, and facilitate competition among different RNAs. The register of Crc in the ternary assembly could be defined by the recognition of a guanine-specific base–phosphate interaction between the first and last ARN repeats of the bound RNA. This dynamic substrate recognition provides structural rationale for the stepwise assembly of multicomponent ribonucleoprotein complexes nucleated by Hfq–RNA binding.

Macromolecular room temperature crystallography

X-ray crystallography enables detailed structural studies of proteins to understand and modulate their function. Conducting crystallographic experiments at cryogenic temperatures has practical benefits but potentially limits the identification of functionally important alternative protein conformations that can be revealed only at room temperature (RT). This review discusses practical aspects of preparing, acquiring, and analyzing X-ray crystallography data at RT to demystify preconceived impracticalities that freeze progress of routine RT data collection at synchrotron sources. Examples are presented as conceptual and experimental templates to enable the design of RT-inspired studies; they illustrate the diversity and utility of gaining novel insights into protein conformational landscapes. An integrative view of protein conformational dynamics enables opportunities to advance basic and biomedical research.

qFit 3: Protein and ligand multiconformer modeling for X-ray crystallographic and single-particle cryo-EM density maps

New X-ray crystallography and cryo-electron microscopy (cryo-EM) approaches yield vast amounts of structural data from dynamic proteins and their complexes. Modeling the full conformational ensemble can provide important biological insights, but identifying and modeling an internally consistent set of alternate conformations remains a formidable challenge. qFit efficiently automates this process by generating a parsimonious multiconformer model. We refactored qFit from a distributed application into software that runs efficiently on a small server, desktop, or laptop. We describe the new qFit 3 software and provide some examples. qFit 3 is open-source under the MIT license, and is available at https://github.com/ExcitedStates/qfit-3.0.

Instrumentation and experimental procedures for robust collection of X-ray diffraction data from protein crystals across physiological temperatures

Traditional X-ray diffraction data collected at cryo-temperatures have delivered invaluable insights into the three-dimensional structures of proteins, providing the backbone of structure-function studies. While cryo-cooling mitigates radiation damage, cryo-temperatures can alter protein conformational ensembles and solvent structure. Furthermore, conformational ensembles underlie protein function and energetics, and recent advances in room-temperature X-ray crystallography have delivered conformational heterogeneity information that can be directly related to biological function. Given this capability, the next challenge is to develop a robust and broadly applicable method to collect single-crystal X-ray diffraction data at and above room temperature. This challenge is addressed herein. The approach described provides complete diffraction data sets with total collection times as short as ∼5 s from single protein crystals, dramatically increasing the quantity of data that can be collected within allocated synchrotron beam time. Its applicability was demonstrated by collecting 1.09-1.54 Å resolution data over a temperature range of 293-363 K for proteinase K, thaumatin and lysozyme crystals at BL14-1 at the Stanford Synchrotron Radiation Lightsource. The analyses presented here indicate that the diffraction data are of high quality and do not suffer from excessive dehydration or radiation damage.

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.

Major conformational changes in the structure of lysozyme obtained from a crystal with a very low solvent content

A new crystal form of lysozyme with a very low solvent content (26.35%) has been obtained in the orthorhombic space group P212121 (with unit-cell parameters a = 30.04, b = 51.68, c = 61.53 Å). The lysozyme structure obtained from these crystals does not show the typical overall fold. Instead, major conformational changes take place in some elements of the secondary structure and in the hydrophobic core of the protein. At the end of the central α-helix (α2), Glu35 is usually buried in the catalytic site and shows an abnormally high pKa value, which is key to the activity of the enzyme. The high pKa value of this glutamate residue is favoured by the hydrophobic environment, particularly by its neighbour Trp108, which is important for structural stability and saccharide binding. In this new structure, Trp108 shows a 90° rotation of its side chain, which results in the rearrangement of the hydrophobic core. Conformational changes also result in the exposure of Glu35 to the solvent, which impairs the catalytic site by increasing the distance between Glu35 and Asp52 and lowering the pKa value of the glutamate. Altogether, this new lysozyme structure reveals major conformational changes in the hydrophobic core and catalytic site that might play a role in the folding and bactericidal function of the protein.

Structural Evidence of Photoisomerization Pathways in Fluorescent Proteins

Double-bond photoisomerization in molecules such as the green fluorescent protein (GFP) chromophore can occur either via a volume-demanding one-bond-flip pathway or via a volume-conserving hula-twist pathway. Understanding the factors that determine the pathway of photoisomerization would inform the rational design of photoswitchable GFPs as improved tools for super-resolution microscopy. In this communication, we reveal the photoisomerization pathway of a photoswitchable GFP, rsEGFP2, by solving crystal structures of cis and trans rsEGFP2 containing a monochlorinated chromophore. The position of the chlorine substituent in the trans state breaks the symmetry of the phenolate ring of the chromophore and allows us to distinguish the two pathways. Surprisingly, we find that the pathway depends on the arrangement of protein monomers within the crystal lattice: in a looser packing, the one-bond-flip occurs, whereas, in a tighter packing (7% smaller unit cell size), the hula-twist occurs.

A comparison of gas stream cooling and plunge cooling of macromolecular crystals

Cryocooling for macromolecular crystallography is usually performed via plunging the crystal into a liquid cryogen or placing the crystal in a cold gas stream. These two approaches are compared here for the case of nitro-gen cooling. The results show that gas stream cooling, which typically cools the crystal more slowly, yields lower mosaicity and, in some cases, a stronger anomalous signal relative to rapid plunge cooling. During plunging, moving the crystal slowly through the cold gas layer above the liquid surface can produce mosaicity similar to gas stream cooling. Annealing plunge cooled crystals by warming and recooling in the gas stream allows the mosaicity and anomalous signal to recover. For tetragonal thermolysin, the observed effects are less pronounced when the cryosolvent has smaller thermal contraction, under which conditions the protein structures from plunge cooled and gas stream cooled crystals are very similar. Finally, this work also demonstrates that the resolution dependence of the reflecting range is correlated with the cooling method, suggesting it may be a useful tool for discerning whether crystals are cooled too rapidly. The results support previous studies suggesting that slower cooling methods are less deleterious to crystal order, as long as ice formation is prevented and dehydration is limited.

A temperature-controlled cold-gas humidifier and its application to protein crystals with the humid-air and glue-coating method

A new temperature-controllable humidifier for X-ray diffraction has been developed. It is shown that the humidifier can successfully maintain protein crystal growth at a temperature lower than room temperature.

The room-temperature experiment has been revisited for macromolecular crystallography. Despite being limited by radiation damage, such experiments reveal structural differences depending on temperature, and it is expected that they will be able to probe structures that are physiologically alive. For such experiments, the humid-air and glue-coating (HAG) method for humidity-controlled experiments is proposed. The HAG method improves the stability of most crystals in capillary-free experiments and is applicable at both cryogenic and ambient temperatures. To expand the thermal versatility of the HAG method, a new humidifier and a protein-crystal-handling workbench have been developed. The devices provide temperatures down to 4°C and successfully maintain growth at that temperature of bovine cytochrome c oxidase crystals, which are highly sensitive to temperature variation. Hence, the humidifier and protein-crystal-handling workbench have proved useful for temperature-sensitive samples and will help reveal temperature-dependent variations in protein structures.

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.

In situ detection of a novel lysozyme monoclinic crystal form upon controlled relative humidity variation

The effect of relative humidity (rH) on protein crystal structures, an area that has attracted high scientific interest during the past decade, is investigated in this study on hen egg-white lysozyme (HEWL) polycrystalline precipitates via in situ laboratory X-ray powder diffraction (XRPD) measurements. For this purpose, HEWL was crystallized at room temperature and pH 4.5, leading to a novel monoclinic HEWL phase which, to our knowledge, has not been reported before. Analysis of XRPD data collected upon rH variation revealed several structural modifications. These observations, on a well-studied molecule like HEWL, underline not only the high impact of humidity levels on biological crystal structures, but also the significance of in-house XRPD as an analytical tool in industrial drug development and its potential to provide information for enhancing manufacturing of pharmaceuticals.