Room-temperature X-ray Crystallography Research Articles

Going around the Kok cycle of the water oxidation reaction with femtosecond X-ray crystallography.

The water oxidation reaction in photosystem II (PS II) produces most of the molecular oxygen in the atmosphere, which sustains life on Earth, and in this process releases four electrons and four protons that drive the downstream process of CO2 fixation in the photosynthetic apparatus. The catalytic center of PS II is an oxygen-bridged Mn4Ca complex (Mn4CaO5) which is progressively oxidized upon the absorption of light by the chlorophyll of the PS II reaction center, and the accumulation of four oxidative equivalents in the catalytic center results in the oxidation of two waters to dioxygen in the last step. The recent emergence of X-ray free-electron lasers (XFELs) with intense femtosecond X-ray pulses has opened up opportunities to visualize this reaction in PS II as it proceeds through the catalytic cycle. In this review, we summarize our recent studies of the catalytic reaction in PS II by following the structural changes along the reaction pathway via room-temperature X-ray crystallography using XFELs. The evolution of the electron density changes at the Mn complex reveals notable structural changes, including the insertion of OX from a new water molecule, which disappears on completion of the reaction, implicating it in the O-O bond formation reaction. We were also able to follow the structural dynamics of the protein coordinating with the catalytic complex and of channels within the protein that are important for substrate and product transport, revealing well orchestrated conformational changes in response to the electronic changes at the Mn4Ca cluster.

Computational remodeling of an enzyme conformational landscape for altered substrate selectivity

Structural plasticity of enzymes dictates their function. Yet, our ability to rationally remodel enzyme conformational landscapes to tailor catalytic properties remains limited. Here, we report a computational procedure for tuning conformational landscapes that is based on multistate design of hinge-mediated domain motions. Using this method, we redesign the conformational landscape of a natural aminotransferase to preferentially stabilize a less populated but reactive conformation and thereby increase catalytic efficiency with a non-native substrate, resulting in altered substrate selectivity. Steady-state kinetics of designed variants reveals activity increases with the non-native substrate of approximately 100-fold and selectivity switches of up to 1900-fold. Structural analyses by room-temperature X-ray crystallography and multitemperature nuclear magnetic resonance spectroscopy confirm that conformational equilibria favor the target conformation. Our computational approach opens the door to targeted alterations of conformational states and equilibria, which should facilitate the design of biocatalysts with customized activity and selectivity.

Revealing protonation states and tracking substrate in serine hydroxymethyltransferase with room-temperature X-ray and neutron crystallography

Pyridoxal 5’-phosphate (PLP)-dependent enzymes utilize a vitamin B6-derived cofactor to perform a myriad of chemical transformations on amino acids and other small molecules. Some PLP-dependent enzymes, such as serine hydroxymethyltransferase (SHMT), are promising drug targets for the design of small-molecule antimicrobials and anticancer therapeutics, while others have been used to synthesize pharmaceutical building blocks. Understanding PLP-dependent catalysis and the reaction specificity is crucial to advance structure-assisted drug design and enzyme engineering. Here we report the direct determination of the protonation states in the active site of Thermus thermophilus SHMT (TthSHMT) in the internal aldimine state using room-temperature joint X-ray/neutron crystallography. Conserved active site architecture of the model enzyme TthSHMT and of human mitochondrial SHMT (hSHMT2) were compared by obtaining a room-temperature X-ray structure of hSHMT2, suggesting identical protonation states in the human enzyme. The amino acid substrate serine pathway through the TthSHMT active site cavity was tracked, revealing the peripheral and cationic binding sites that correspond to the pre-Michaelis and pseudo-Michaelis complexes, respectively. At the peripheral binding site, the substrate is bound in the zwitterionic form. By analyzing the observed protonation states, Glu53, but not His residues, is proposed as the general base catalyst, orchestrating the retro-aldol transformation of L-serine into glycine.

Contribution of the catalytic dyad of SARS-CoV-2 main protease to binding covalent and noncovalent inhibitors

The effect of mutations of the catalytic dyad residues of SARS-CoV-2 main protease (MProWT) on the thermodynamics of binding of covalent inhibitors comprising nitrile [nirmatrelvir (NMV), NBH2], aldehyde (GC373), and ketone (BBH1) warheads to MPro is examined together with room temperature X-ray crystallography. When lacking the nucleophilic C145, NMV binding is ∼400-fold weaker corresponding to 3.5kcal/mol and 13.3 °C decrease in free energy (ΔG) and thermal stability (Tm), respectively, relative to MProWT. The H41A mutation results in a 20-fold increase in the dissociation constant (Kd), and 1.7kcal/mol and 1.4 °C decreases in ΔG and Tm, respectively. Increasing the pH from 7.2 to 8.2 enhances NMV binding to MProH41A, whereas no significant change is observed in binding to MProWT. Structures of the four inhibitor complexes with MPro1-304/C145A show that the active site geometries of the complexes are nearly identical to that of MProWT with the nucleophilic sulfur of C145 positioned to react with the nitrile or the carbonyl carbon. These results support a two-step mechanism for the formation of the covalent complex involving an initial non-covalent binding followed by a nucleophilic attack by the thiolate anion of C145 on the warhead carbon. Noncovalent inhibitor ensitrelvir (ESV) exhibits a binding affinity to MProWT that is similar to NMV but differs in its thermodynamic signature from NMV. The binding of ESV to MProC145A also results in a significant, but smaller, increase in Kd and decrease in ΔG and Tm, relative to NMV.

AI-Accelerated Design of Targeted Covalent Inhibitors for SARS-CoV-2.

Direct-acting antivirals for the treatment of the COVID-19 pandemic caused by the SARS-CoV-2 virus are needed to complement vaccination efforts. Given the ongoing emergence of new variants, automated experimentation, and active learning based fast workflows for antiviral lead discovery remain critical to our ability to address the pandemic's evolution in a timely manner. While several such pipelines have been introduced to discover candidates with noncovalent interactions with the main protease (Mpro), here we developed a closed-loop artificial intelligence pipeline to design electrophilic warhead-based covalent candidates. This work introduces a deep learning-assisted automated computational workflow to introduce linkers and an electrophilic "warhead" to design covalent candidates and incorporates cutting-edge experimental techniques for validation. Using this process, promising candidates in the library were screened, and several potential hits were identified and tested experimentally using native mass spectrometry and fluorescence resonance energy transfer (FRET)-based screening assays. We identified four chloroacetamide-based covalent inhibitors of Mpro with micromolar affinities (KI of 5.27 μM) using our pipeline. Experimentally resolved binding modes for each compound were determined using room-temperature X-ray crystallography, which is consistent with the predicted poses. The induced conformational changes based on molecular dynamics simulations further suggest that the dynamics may be an important factor to further improve selectivity, thereby effectively lowering KI and reducing toxicity. These results demonstrate the utility of our modular and data-driven approach for potent and selective covalent inhibitor discovery and provide a platform to apply it to other emerging targets.

Room-temperature serial synchrotron crystallography of the human phosphatase PTP1B.

Room-temperature X-ray crystallography provides unique insights into protein conformational heterogeneity, but obtaining sufficiently large protein crystals is a common hurdle. Serial synchrotron crystallography (SSX) helps to address this hurdle by allowing the use of many medium- to small-sized crystals. Here, a recently introduced serial sample-support chip system has been used to obtain the first SSX structure of a human phosphatase, specifically protein tyrosine phosphatase 1B (PTP1B) in the unliganded (apo) state. In previous apo room-temperature structures, the active site and allosteric sites adopted alternate conformations, including open and closed conformations of the active-site WPD loop and of a distal allosteric site. By contrast, in our SSX structure the active site is best fitted with a single conformation, but the distal allosteric site is best fitted with alternate conformations. This observation argues for additional nuance in interpreting the nature of allosteric coupling in this protein. Overall, our results illustrate the promise of serial methods for room-temperature crystallography, as well as future avant-garde crystallography experiments, for PTP1B and other proteins.

Unraveling Complex Local Protein Environments with 4-Cyano-l-phenylalanine.

We present a multifaceted approach to effectively probe complex local protein environments utilizing the vibrational reporter unnatural amino acid (UAA) 4-cyano-l-phenylalanine (pCNPhe) in the model system superfolder green fluorescent protein (sfGFP). This approach combines temperature-dependent infrared (IR) spectroscopy, X-ray crystallography, and molecular dynamics (MD) simulations to provide a molecular interpretation of the local environment of the nitrile group in the protein. Specifically, a two-step enantioselective synthesis was developed that provided an 87% overall yield of pCNPhe in high purity without the need for chromatography. It was then genetically incorporated individually at three unique sites (74, 133, and 149) in sfGFP to probe these local protein environments. The incorporation of the UAA site-specifically in sfGFP utilized an engineered, orthogonal tRNA synthetase in E. coli using the Amber codon suppression protocol, and the resulting UAA-containing sfGFP constructs were then explored with this approach. This methodology was effectively utilized to further probe the local environments of two surface sites (sites 133 and 149) that we previously explored with room temperature IR spectroscopy and X-ray crystallography and a new interior site (site 74) featuring a complex local environment around the nitrile group of pCNPhe. Site 133 was found to be solvent-exposed, while site 149 was partially buried. Site 74 was found to consist of three distinct local environments around the nitrile group including nonspecific van der Waals interactions, hydrogen-bonding to a structural water, and hydrogen-bonding to a histidine side chain.

Ensemble-function relationships to dissect mechanisms of enzyme catalysis.

Decades of structure-function studies have established our current extensive understanding of enzymes. However, traditional structural models are snapshots of broader conformational ensembles of interchanging states. We demonstrate the need for conformational ensembles to understand function, using the enzyme ketosteroid isomerase (KSI) as an example. Comparison of prior KSI cryogenic x-ray structures suggested deleterious mutational effects from a misaligned oxyanion hole catalytic residue. However, ensemble information from room-temperature x-ray crystallography, combined with functional studies, excluded this model. Ensemble-function analyses can deconvolute effects from altering the probability of occupying a state (P-effects) and changing the reactivity of each state (k-effects); our ensemble-function analyses revealed functional effects arising from weakened oxyanion hole hydrogen bonding and substrate repositioning within the active site. Ensemble-function studies will have an integral role in understanding enzymes and in meeting the future goals of a predictive understanding of enzyme catalysis and engineering new enzymes.

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.

High-Throughput Virtual Screening and Validation of a SARS-CoV-2 Main Protease Noncovalent Inhibitor.

Despite the recent availability of vaccines against the acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the search for inhibitory therapeutic agents has assumed importance especially in the context of emerging new viral variants. In this paper, we describe the discovery of a novel noncovalent small-molecule inhibitor, MCULE-5948770040, that binds to and inhibits the SARS-Cov-2 main protease (Mpro) by employing a scalable high-throughput virtual screening (HTVS) framework and a targeted compound library of over 6.5 million molecules that could be readily ordered and purchased. Our HTVS framework leverages the U.S. supercomputing infrastructure achieving nearly 91% resource utilization and nearly 126 million docking calculations per hour. Downstream biochemical assays validate this Mpro inhibitor with an inhibition constant (Ki) of 2.9 μM (95% CI 2.2, 4.0). Furthermore, using room-temperature X-ray crystallography, we show that MCULE-5948770040 binds to a cleft in the primary binding site of Mpro forming stable hydrogen bond and hydrophobic interactions. We then used multiple μs-time scale molecular dynamics (MD) simulations and machine learning (ML) techniques to elucidate how the bound ligand alters the conformational states accessed by Mpro, involving motions both proximal and distal to the binding site. Together, our results demonstrate how MCULE-5948770040 inhibits Mpro and offers a springboard for further therapeutic design.

Michaelis-like complex of SARS-CoV-2 main protease visualized by room-temperature X-ray crystallography.

SARS-CoV-2 emerged at the end of 2019 to cause an unprecedented pandemic of the deadly respiratory disease COVID-19 that continues to date. The viral main protease (Mpro) is essential for SARS-CoV-2 replication and is therefore an important drug target. Understanding the catalytic mechanism of Mpro, a cysteine protease with a catalytic site comprising the noncanonical Cys145-His41 dyad, can help in guiding drug design. Here, a 2.0 Å resolution room-temperature X-ray crystal structure is reported of a Michaelis-like complex of Mpro harboring a single inactivating mutation C145A bound to the octapeptide Ac-SAVLQSGF-CONH2 corresponding to the nsp4/nsp5 autocleavage site. The peptide substrate is unambiguously defined in subsites S5 to S3' by strong electron density. Superposition of the Michaelis-like complex with the neutron structure of substrate-free Mpro demonstrates that the catalytic site is inherently pre-organized for catalysis prior to substrate binding. Induced fit to the substrate is driven by P1 Gln binding in the predetermined subsite S1 and rearrangement of subsite S2 to accommodate P2 Leu. The Michaelis-like complex structure is ideal for in silico modeling of the SARS-CoV-2 Mpro catalytic mechanism.

A capillary-based microfluidic device enables primary high-throughput room-temperature crystallographic screening.

A novel capillary-based microfluidic strategy to accelerate the process of small-molecule-compound screening by room-temperature X-ray crystallography using protein crystals is reported. The ultra-thin microfluidic devices are composed of a UV-curable polymer, patterned by cleanroom photolithography, and have nine capillary channels per chip. The chip was designed for ease of sample manipulation, sample stability and minimal X-ray background. 3D-printed frames and cassettes conforming to SBS standards are used to house the capillary chips, providing additional mechanical stability and compatibility with automated liquid- and sample-handling robotics. These devices enable an innovative in situ crystal-soaking screening workflow, akin to high-throughput compound screening, such that quantitative electron density maps sufficient to determine weak binding events are efficiently obtained. This work paves the way for adopting a room-temperature microfluidics-based sample delivery method at synchrotron sources to facilitate high-throughput protein-crystallography-based screening of compounds at high concentration with the aim of discovering novel binding events in an automated manner.

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.

Malleability of the SARS-CoV-2 3CL Mpro Active-Site Cavity Facilitates Binding of Clinical Antivirals

The COVID-19 pandemic caused by SARS-CoV-2 requires rapid development of specific therapeutics and vaccines. The main protease of SARS-CoV-2, 3CL Mpro, is an established drug target for the design of inhibitors to stop the virus replication. Repurposing existing clinical drugs can offer a faster route to treatments. Here, we report on the binding mode and inhibition properties of several inhibitors using room temperature X-ray crystallography and in vitro enzyme kinetics. The enzyme active-site cavity reveals a high degree of malleability, allowing aldehyde leupeptin and hepatitis C clinical protease inhibitors (telaprevir, narlaprevir, and boceprevir) to bind and inhibit SARS-CoV-2 3CL Mpro. Narlaprevir, boceprevir, and telaprevir are low-micromolar inhibitors, whereas the binding affinity of leupeptin is substantially weaker. Repurposing hepatitis C clinical drugs as COVID-19 treatments may be a useful option to pursue. The observed malleability of the enzyme active-site cavity should be considered for the successful design of specific protease inhibitors.

Ensemble-based enzyme design can recapitulate the effects of laboratory directed evolution in silico

The creation of artificial enzymes is a key objective of computational protein design. Although de novo enzymes have been successfully designed, these exhibit low catalytic efficiencies, requiring directed evolution to improve activity. Here, we use room-temperature X-ray crystallography to study changes in the conformational ensemble during evolution of the designed Kemp eliminase HG3 (kcat/KM 146 M−1s−1). We observe that catalytic residues are increasingly rigidified, the active site becomes better pre-organized, and its entrance is widened. Based on these observations, we engineer HG4, an efficient biocatalyst (kcat/KM 103,000 M−1s−1) containing key first and second-shell mutations found during evolution. HG4 structures reveal that its active site is pre-organized and rigidified for efficient catalysis. Our results show how directed evolution circumvents challenges inherent to enzyme design by shifting conformational ensembles to favor catalytically-productive sub-states, and suggest improvements to the design methodology that incorporate ensemble modeling of crystallographic data.

Room-temperature X-ray crystallography reveals the oxidation and reactivity of cysteine residues in SARS-CoV-2 3CL Mpro: insights into enzyme mechanism and drug design.

The emergence of the novel coronavirus SARS-CoV-2 has resulted in a worldwide pandemic not seen in generations. Creating treatments and vaccines to battle COVID-19, the disease caused by the virus, is of paramount importance in order to stop its spread and save lives. The viral main protease, 3CL Mpro, is indispensable for the replication of SARS-CoV-2 and is therefore an important target for the design of specific protease inhibitors. Detailed knowledge of the structure and function of 3CL Mpro is crucial to guide structure-aided and computational drug-design efforts. Here, the oxidation and reactivity of the cysteine residues of the protease are reported using room-temperature X-ray crystallography, revealing that the catalytic Cys145 can be trapped in the peroxysulfenic acid oxidation state at physiological pH, while the other surface cysteines remain reduced. Only Cys145 and Cys156 react with the alkylating agent N-ethylmaleimide. It is suggested that the zwitterionic Cys145-His45 catalytic dyad is the reactive species that initiates catalysis, rather than Cys145-to-His41 proton transfer via the general acid-base mechanism upon substrate binding. The structures also provide insight into the design of improved 3CL Mpro inhibitors.

Structural plasticity of SARS-CoV-2 3CL Mpro active site cavity revealed by room temperature X-ray crystallography

The COVID-19 disease caused by the SARS-CoV-2 coronavirus has become a pandemic health crisis. An attractive target for antiviral inhibitors is the main protease 3CL Mpro due to its essential role in processing the polyproteins translated from viral RNA. Here we report the room temperature X-ray structure of unliganded SARS-CoV-2 3CL Mpro, revealing the ligand-free structure of the active site and the conformation of the catalytic site cavity at near-physiological temperature. Comparison with previously reported low-temperature ligand-free and inhibitor-bound structures suggest that the room temperature structure may provide more relevant information at physiological temperatures for aiding in molecular docking studies.

Room-temperature photo-induced martensitic transformation in a protein crystal.

Martensitic transformations are the first-order crystal-to-crystal phase transitions that occur mostly in materials such as steel, alloys and ceramics, thus having many technological applications. These phase transitions are rarely observed in molecular crystals and have not been detected in protein crystals. Reversibly switchable fluorescent proteins are widely used in biotechnology, including super-resolution molecular imaging, and hold promise as candidate biomaterials for future high-tech applications. Here, we report on a reversibly switchable fluorescent protein, Tetdron, whose crystals undergo a photo-induced martensitic transformation at room temperature. Room-temperature X-ray crystallography demonstrates that at equilibrium Tetdron chromophores are all in the trans configuration, with an ∼1:1 mixture of their protonated and deprotonated forms. Irradiation of a Tetdron crystal with 400 nm light induces a martensitic transformation, which results in Tetdron tetramerization at room temperature revealed by X-ray photocrystallography. Crystal and solution spectroscopic measurements provide evidence that the photo-induced martensitic phase transition is coupled with the chromophore deprotonation, but no trans-cis isomerization is detected in the structure of an irradiated crystal. It is hypothesized that protein dynamics assists in the light-induced proton transfer from the chromophore to the bulk solvent and in the ensuing martensitic phase transition. The unique properties of Tetdron may be useful in developing novel biomaterials for optogenetics, data storage and nanotechnology.

Biophysical Characterization of a Disabled Double Mutant of Soybean Lipoxygenase: The "Undoing" of Precise Substrate Positioning Relative to Metal Cofactor and an Identified Dynamical Network.

Soybean lipoxygenase (SLO) has served as a prototype for understanding the molecular origin of enzymatic rate accelerations. The double mutant (DM) L546A/L754A is considered a dramatic outlier, due to the unprecedented size and near temperature-independence of its primary kinetic isotope effect, low catalytic efficiency, and elevated enthalpy of activation. To uncover the physical basis of these features, we herein apply three structural probes: hydrogen-deuterium exchange mass spectrometry, room-temperature X-ray crystallography and EPR spectroscopy on four SLO variants (wild-type (WT) enzyme, DM, and the two parental single mutants, L546A and L754A). DM is found to incorporate features of each parent, with the perturbation at position 546 predominantly influencing thermally activated motions that connect the active site to a protein-solvent interface, while mutation at position 754 disrupts the ligand field and solvation near the cofactor iron. However, the expanded active site in DM leads to more active site water molecules and their associated hydrogen bond network, and the individual features from L546A and L754A alone cannot explain the aggregate kinetic properties for DM. Using recently published QM/MM-derived ground-state SLO-substrate complexes for WT and DM, together with the thorough structural analyses presented herein, we propose that the impairment of DM is the combined result of a repositioning of the reactive carbon of linoleic acid substrate with regard to both the iron cofactor and a catalytically linked dynamic region of protein.

Rescue of conformational dynamics in enzyme catalysis by directed evolution

Rational design and directed evolution have proved to be successful approaches to increase catalytic efficiencies of both natural and artificial enzymes. Protein dynamics is recognized as important, but due to the inherent flexibility of biological macromolecules it is often difficult to distinguish which conformational changes are directly related to function. Here, we use directed evolution on an impaired mutant of the proline isomerase CypA and identify two second-shell mutations that partially restore its catalytic activity. We show both kinetically, using NMR spectroscopy, and structurally, by room-temperature X-ray crystallography, how local perturbations propagate through a large allosteric network to facilitate conformational dynamics. The increased catalysis selected for in the evolutionary screen is correlated with an accelerated interconversion between the two catalytically essential conformational sub-states, which are both captured in the high-resolution X-ray ensembles. Our data provide a glimpse of an evolutionary trajectory and show how subtle changes can fine-tune enzyme function.