Accessing protein conformational ensembles using room-temperature X-ray crystallography

Unveiling Putative Excited State and Transmission of Binding Information in the Fluoride Riboswitch

Riboswitches regulate downstream gene expression by binding to specific small molecules or ions with multiple mechanisms to transfer the binding information. In the case of the fluoride riboswitch, the transcription termination signal is conveyed through a transient excited state (ES). In this work, we performed conventional molecular dynamics (MD) simulations, totaling 180 μs, to obtain the ES structure and investigate the mechanism underlying information transmission in Mg2+/F- binding within the fluoride riboswitch aptamer. The Mg2+/F- binding pocket exhibits various conformations in its apo form. A series of ES structures were extracted from the MD trajectories of the apo form. The dynamics of the Mg2+/F- binding pocket influenced key pair A40-U48 in ES structures. The pathway connecting the binding pocket to the pair involves interactions between the phosphate groups of U7 and G8 and the nucleobases of G8-C47-U48. Our work presents a structural ensemble of the ES and elucidates a pathway for transferring Mg2+/F- binding information, thereby facilitating the understanding of how the holo-like apo state achieves transcriptional repression.

Structural and functional effects of phosphopriming and scaffolding in the kinase GSK-3β

Glycogen synthase kinase 3β (GSK-3β) targets specific signaling pathways in response to distinct upstream signals. We used structural and functional studies to dissect how an upstream phosphorylation step primes the Wnt signaling component β-catenin for phosphorylation by GSK-3β and how scaffolding interactions contribute to this reaction. Our crystal structure of GSK-3β bound to a phosphoprimed β-catenin peptide confirmed the expected binding mode of the phosphoprimed residue adjacent to the catalytic site. An aspartate phosphomimic in the priming site of β-catenin adopted an indistinguishable structure but reacted approximately 1000-fold slower than the native phosphoprimed substrate. This result suggests that substrate positioning alone is not sufficient for catalysis and that native phosphopriming interactions are necessary. We also obtained a structure of GSK-3β with an extended peptide from the scaffold protein Axin that bound with greater affinity than that of previously crystallized Axin fragments. This structure neither revealed additional contacts that produce the higher affinity nor explained how substrate interactions in the GSK-3β active site are modulated by remote Axin binding. Together, our findings suggest that phosphopriming and scaffolding produce small conformational changes or allosteric effects, not captured in the crystal structures, that activate GSK-3β and facilitate β-catenin phosphorylation. These results highlight limitations in our ability to predict catalytic activity from structure and have potential implications for the role of natural phosphomimic mutations in kinase regulation and phosphosite evolution.

In situ serial crystallography facilitates 96-well plate structural analysis at low symmetry

The advent of serial crystallography has rejuvenated and popularized room-temperature X-ray crystal structure determination. Structures determined at physiological temperature reveal protein flexibility and dynamics. In addition, challenging samples (e.g. large complexes, membrane proteins and viruses) form fragile crystals that are often difficult to harvest for cryo-crystallography. Moreover, a typical serial crystallography experiment requires a large number of microcrystals, mainly achievable through batch crystallization. Many medically relevant samples are expressed in mammalian cell lines, producing a meager quantity of protein that is incompatible with batch crystallization. This can limit the scope of serial crystallography approaches. Direct in situ data collection from a 96-well crystallization plate enables not only the identification of the best diffracting crystallization condition but also the possibility for structure determination under ambient conditions. Here, we describe an in situ serial crystallography (iSX) approach, facilitating direct measurement from crystallization plates mounted on a rapidly exchangeable universal plate holder deployed at a microfocus beamline, ID23-2, at the European Synchrotron Radiation Facility. We applied our iSX approach on a challenging project, autotaxin, a therapeutic target expressed in a stable human cell line, to determine the structure in the lowest-symmetry P1 space group at 3.0 Å resolution. Our in situ data collection strategy provided a complete dataset for structure determination while screening various crystallization conditions. Our data analysis reveals that the iSX approach is highly efficient at a microfocus beamline, improving throughput and demonstrating how crystallization plates can be routinely used as an alternative method of presenting samples for serial crystallography experiments at synchrotrons.

The probability of cryptic pocket opening controls functional tradeoffs in filovirus immune evasion

Cryptic pockets are of growing interest as potential drug targets. However, it remains unclear whether they contribute to protein function or if they are merely happenstantial features that can easily be evolved away to achieve drug resistance. Here, we explore whether a cryptic pocket in the Interferon Inhibitory Domain (IID) of viral protein 35 (VP35) of Zaire ebolavirus has a functional role. We use simulations and experiments to study the relationship between cryptic pocket opening and double-stranded RNA (dsRNA) binding of the IIDs of two other filoviruses, Reston and Marburg. These homologs have structures nearly identical to Zaire but block different interferon pathways because Marburg IID only binds the backbone, while Zaire and Reston IIDs bind both backbone and blunt ends of dsRNA. We hypothesized that these differences arise from alterations in their pocket dynamics. Simulations and thiol-labeling experiments demonstrate that Reston has a lower probability of opening the cryptic pocket while Marburg has a higher probability than Zaire. Subsequent dsRNA-binding assays with different length substrates suggest that closed conformations preferentially bind dsRNA blunt ends while open conformations prefer binding the backbone. Further, a point mutation that increases the probability of cryptic pocket opening in Reston shows that the open states prefer backbone binding, while a substitution that lowers the probability of pocket opening in Zaire improves binding to blunt ends. These results demonstrate the cryptic pocket controls a functional tradeoff, suggesting cryptic pockets are under selective pressure and may be difficult to evolve away to achieve drug resistance.

A structural perspective on the temperature-dependent activity of enzymes

Enzymes are biomolecular catalysts whose activity varies with temperature. Unlike for small-molecule catalysts, the structural ensembles of enzymes can vary substantially with temperature, and it is in general unclear how this modulates the temperature dependence of activity. Here multi-temperature X-ray crystallography was used to record structural changes from -20°C to 40°C for a mesophilic enzyme in complex with inhibitors mimicking substrate-, intermediate-, and product-bound states, representative of major complexes underlying the kinetic constantk c a t . Both inhibitors, substrates and catalytically relevant loop motifs increasingly populate catalytically competent conformations as temperature increases. These changes occur even in temperature ranges where kinetic measurements show roughly linear Arrhenius/Eyring behavior where parameters characterizing the system are assumed to be temperature independent. Simple analysis shows that linear Arrhenius/Eyring behavior can still be observed when the underlying activation energy / enthalpy values vary with temperature, e.g., due to structural changes, and that the underlying thermodynamic parameters can be far from values derived from Arrhenius/Eyring model fits. Our results indicate a critical role for temperature-dependent atomic-resolution structural data in interpreting temperature-dependent kinetic data from enzymatic systems.

An expanded trove of fragment-bound structures for the allosteric enzyme PTP1B from computational reanalysis of large-scale crystallographic data

Due to their low binding affinities, detecting small-molecule fragments bound to protein structures from crystallographic datasets has been a challenge. Here, we report a trove of 65 new fragment hits for PTP1B, an "undruggable" therapeutic target enzyme for diabetes and cancer. These structures were obtained from computational analysis of data from a large crystallographic screen, demonstrating the power of this approach to elucidate many (∼50% more) "hidden" ligand-bound states of proteins. Our new structures include a fragment hit found in a novel binding site in PTP1B with a unique location relative to the active site, one that links adjacent allosteric sites, and, perhaps most strikingly, a fragment that induces long-range allosteric protein conformational responses. Altogether, our research highlights the utility of computational analysis of crystallographic data, makes publicly available dozens of new ligand-bound structures of a high-value drug target, and identifies novel aspects of ligandability and allostery in PTP1B.

Cryo2RT: a high-throughput method for room-temperature macromolecular crystallography from cryo-cooled crystals

Advances in structural biology have relied heavily on synchrotron cryo-crystallography and cryogenic electron microscopy to elucidate biological processes and for drug discovery. However, disparities between cryogenic and room-temperature (RT) crystal structures pose challenges. Here, Cryo2RT, a high-throughput RT data-collection method from cryo-cooled crystals that leverages the cryo-crystallography workflow, is introduced. Tested on endothiapepsin crystals with four soaked fragments, thaumatin and SARS-CoV-2 3CLpro, Cryo2RT reveals unique ligand-binding poses, offers a comparable throughput to cryo-crystallography and eases the exploration of structural dynamics at various temperatures.

A snapshot love story: what serial crystallography has done and will do for us

Serial crystallography, born from groundbreaking experiments at the Linac Coherent Light Source in 2009, has evolved into a pivotal technique in structural biology. Initially pioneered at X-ray free-electron laser facilities, it has now expanded to synchrotron-radiation facilities globally, with dedicated experimental stations enhancing its accessibility. This review gives an overview of current developments in serial crystallography, emphasizing recent results in time-resolved crystallography, and discussing challenges and shortcomings.

Mapping protein conformational landscapes from crystallographic drug fragment screens

Proteins are dynamic macromolecules. Knowledge of a protein's thermally accessible conformations is critical to determining important transitions and designing therapeutics. Accessible conformations are highly constrained by a protein's structure such that concerted structural changes due to external perturbations likely track intrinsic conformational transitions. These transitions can be thought of as paths through a conformational landscape. Crystallographic drug fragment screens are high-throughput perturbation experiments, in which thousands of crystals of a drug target are soaked with small-molecule drug precursors (fragments) and examined for fragment binding, mapping potential drug binding sites on the target protein. Here, we describe an open-source Python package, COLAV (COnformational LAndscape Visualization), to infer conformational landscapes from such large-scale crystallographic perturbation studies. We apply COLAV to drug fragment screens of two medically important systems: protein tyrosine phosphatase 1B (PTP-1B), which regulates insulin signaling, and the SARS CoV-2 Main Protease (MPro). With enough fragment-bound structures, we find that such drug screens also enable detailed mapping of proteins' conformational landscapes.

A curated rotamer library for common post-translational modifications of proteins

MOTIVATION

Sidechain rotamer libraries of the common amino acids of a protein are useful for folded protein structure determination and for generating ensembles of intrinsically disordered proteins (IDPs). However, much of protein function is modulated beyond the translated sequence through the introduction of post-translational modifications (PTMs).

RESULTS

In this work, we have provided a curated set of side chain rotamers for the most common PTMs derived from the RCSB PDB database, including phosphorylated, methylated, and acetylated sidechains. Our rotamer libraries improve upon existing methods such as SIDEpro, Rosetta, and AlphaFold3 in predicting the experimental structures for PTMs in folded proteins. In addition, we showcase our PTM libraries in full use by generating ensembles with the Monte Carlo Side Chain Entropy (MCSCE) for folded proteins, and combining MCSCE with the Local Disordered Region Sampling algorithms within IDPConformerGenerator for proteins with intrinsically disordered regions.

AVAILABILITY AND IMPLEMENTATION

The codes for dihedral angle computations and library creation are available at https://github.com/THGLab/ptm_sc.git.

Fast Quantitative Validation of 3D Models of Low-Affinity Protein-Ligand Complexes by STD NMR Spectroscopy

Low-affinity protein-ligand interactions are important for many biological processes, including cell communication, signal transduction, and immune responses. Structural characterization of these complexes is also critical for the development of new drugs through fragment-based drug discovery (FBDD), but it is challenging due to the low affinity of fragments for the binding site. Saturation transfer difference (STD) NMR spectroscopy has revolutionized the study of low-affinity receptor-ligand interactions enabling binding detection and structural characterization. Comparison of relaxation and exchange matrix calculations with 1H STD NMR experimental data is essential for the validation of 3D structures of protein-ligand complexes. In this work, we present a new approach based on the calculation of a reduced relaxation matrix, in combination with funnel metadynamics MD simulations, that allows a very fast generation of experimentally STD-NMR-validated 3D structures of low-affinity protein-ligand complexes.

Automated multiconformer model building for X-ray crystallography and cryo-EM

In their folded state, biomolecules exchange between multiple conformational states that are crucial for their function. Traditional structural biology methods, such as X-ray crystallography and cryogenic electron microscopy (cryo-EM), produce density maps that are ensemble averages, reflecting molecules in various conformations. Yet, most models derived from these maps explicitly represent only a single conformation, overlooking the complexity of biomolecular structures. To accurately reflect the diversity of biomolecular forms, there is a pressing need to shift toward modeling structural ensembles that mirror the experimental data. However, the challenge of distinguishing signal from noise complicates manual efforts to create these models. In response, we introduce the latest enhancements to qFit, an automated computational strategy designed to incorporate protein conformational heterogeneity into models built into density maps. These algorithmic improvements in qFit are substantiated by superior Rfree and geometry metrics across a wide range of proteins. Importantly, unlike more complex multicopy ensemble models, the multiconformer models produced by qFit can be manually modified in most major model building software (e.g., Coot) and fit can be further improved by refinement using standard pipelines (e.g., Phenix, Refmac, Buster). By reducing the barrier of creating multiconformer models, qFit can foster the development of new hypotheses about the relationship between macromolecular conformational dynamics and function.

High Pressure Promotes Binding of the Allosteric Inhibitor Zn2+-Cyclen in Crystals of Activated H-Ras

In this work, we experimentally investigate the potency of high pressure to drive a protein toward an excited state where an inhibitor targeted for this state can bind. Ras proteins are small GTPases cycling between active GTP-bound and inactive GDP-bound states. Various states of GTP-bound Ras in active conformation coexist in solution, amongst them, state 2 which binds to effectors, and state 1, weakly populated at ambient conditions, which has a low affinity for effectors. Zn2+-cyclen is an allosteric inhibitor of Ras protein, designed to bind specifically to the state 1. In H-Ras(wt).Mg2+.GppNHp crystals soaked with Zn2+-cyclen, no binding could be observed, as expected in the state 2 conformation which is the dominant state at ambient pressure. Interestingly, Zn2+-cyclen binding is observed at 500 MPa pressure, close to the nucleotide, in Ras protein that is driven by pressure to a state 1 conformer. The unknown binding mode of Zn2+-cyclen to H-Ras can thus be fully characterized in atomic details. As a more general conjunction from our study, high pressure x-ray crystallography turns out to be a powerful method to induce transitions allowing drug binding in proteins that are in low-populated conformations at ambient conditions, enabling the design of specific inhibitors.

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.

Analysis of the effect of cations on protein conformational stability using solid-state nanopores

The conformation of proteins is closely related to their biological functions, and it is affected by many factors, including the type of cations in solution. However, it is difficult to detect the conformational changes of a protein in situ. As a single-molecule sensing technology, nanopores can convert molecular structural information into analyzable current signals within a reasonable time range. Herein, we detect and analyze the effects of two different types of monovalent cations (Na+ and Li+) on a model protein bovine serum albumin (BSA) conformation using SiNx nanopores with different diameters. The quantitative analysis results show that the excluded volume of BSA in LiCl salt solutions is larger than the value in NaCl solution, indicating that Li+ is more prone to unfolding the proteins and making them unstable. This study demonstrated that nanopores enable the in situ detection of the structure of proteins at the single-molecule level and provide a new approach for the quantitative analysis of proteins.

A Curated Rotamer Library for Common Post-Translational Modifications of Proteins

Sidechain rotamer libraries of the common amino acids of a protein are useful for folded protein structure determination and for generating ensembles of intrinsically disordered proteins (IDPs). However much of protein function is modulated beyond the translated sequence through thFiguree introduction of post-translational modifications (PTMs). In this work we have provided a curated set of side chain rotamers for the most common PTMs derived from the RCSB PDB database, including phosphorylated, methylated, and acetylated sidechains. Our rotamer libraries improve upon existing methods such as SIDEpro and Rosetta in predicting the experimental structures for PTMs in folded proteins. In addition, we showcase our PTM libraries in full use by generating ensembles with the Monte Carlo Side Chain Entropy (MCSCE) for folded proteins, and combining MCSCE with the Local Disordered Region Sampling algorithms within IDPConformerGenerator for proteins with intrinsically disordered regions.

Identifying and avoiding radiation damage in macromolecular crystallography

Radiation damage remains one of the major impediments to accurate structure solution in macromolecular crystallography. The artefacts of radiation damage can manifest as structural changes that result in incorrect biological interpretations being drawn from a model, they can reduce the resolution to which data can be collected and they can even prevent structure solution entirely. In this article, we discuss how to identify and mitigate against the effects of radiation damage at each stage in the macromolecular crystal structure-solution pipeline.

A delayed decoupling methyl-TROSY pulse sequence for atomic resolution studies of folded proteins and RNAs in condensates

Solution NMR spectroscopy has tremendous potential for providing atomic resolution insights into the interactions between proteins and nucleic acids partitioned into condensed phases of phase-separated systems. However, the highly viscous nature of the condensed phase challenges applications, and in particular, the extraction of quantitative, site-specific information. Here, we present a delayed decoupling-based HMQC pulse sequence for methyl-TROSY studies of 'client' proteins and nucleic acids partitioned into 'scaffold' proteinaceous phase-separated solvents. High sensitivity and excellent quality spectra are recorded of a nascent form of superoxide dismutase and of a small RNA fragment partitioned into CAPRIN1 condensates.

Uncovering Protein Ensembles: Automated Multiconformer Model Building for X-ray Crystallography and Cryo-EM

In their folded state, biomolecules exchange between multiple conformational states that are crucial for their function. Traditional structural biology methods, such as X-ray crystallography and cryogenic electron microscopy (cryo-EM), produce density maps that are ensemble averages, reflecting molecules in various conformations. Yet, most models derived from these maps explicitly represent only a single conformation, overlooking the complexity of biomolecular structures. To accurately reflect the diversity of biomolecular forms, there is a pressing need to shift towards modeling structural ensembles that mirror the experimental data. However, the challenge of distinguishing signal from noise complicates manual efforts to create these models. In response, we introduce the latest enhancements to qFit, an automated computational strategy designed to incorporate protein conformational heterogeneity into models built into density maps. These algorithmic improvements in qFit are substantiated by superior R f r e e and geometry metrics across a wide range of proteins. Importantly, unlike more complex multicopy ensemble models, the multiconformer models produced by qFit can be manually modified in most major model building software (e.g. Coot) and fit can be further improved by refinement using standard pipelines (e.g. Phenix, Refmac, Buster). By reducing the barrier of creating multiconformer models, qFit can foster the development of new hypotheses about the relationship between macromolecular conformational dynamics and function.

Alternate conformations found in protein structures implies biological functions: A case study using cyclophilin A

Protein dynamics linked to numerous biomolecular functions, such as ligand binding, allosteric regulation, and catalysis, must be better understood at the atomic level. Reactive atoms of key residues drive a repertoire of biomolecular functions by flipping between alternate conformations or conformational substates, seldom found in protein structures. Probing such sparsely sampled alternate conformations would provide mechanistic insight into many biological functions. We are therefore interested in evaluating the instance of amino acids adopted alternate conformations, either in backbone or side-chain atoms or in both. Accordingly, over 70000 protein structures appear to contain alternate conformations only 'A' and 'B' for any atom, particularly the instance of amino acids that adopted alternate conformations are more for Arg, Cys, Met, and Ser than others. The resulting protein structure analysis depicts that amino acids with alternate conformations are mainly found in the helical and β-regions and are often seen in high-resolution X-ray crystal structures. Furthermore, a case study on human cyclophilin A (CypA) was performed to explain the pre-existing intrinsic dynamics of catalytically critical residues from the CypA and how such intrinsic dynamics perturbed upon Ser99Thr mutation using molecular dynamics simulations on the ns-μs timescale. Simulation results demonstrated that the Ser99Thr mutation had impaired the alternate conformations or the catalytically productive micro-environment of Phe113, mimicking the experimentally observed perturbation captured by X-ray crystallography. In brief, a deeper comprehension of alternate conformations adopted by the amino acids may shed light on the interplay between protein structure, dynamics, and function.

Functional protein dynamics in a crystal

Proteins are molecular machines and to understand how they work, we need to understand how they move. New pump-probe time-resolved X-ray diffraction methods open up ways to initiate and observe protein motions with atomistic detail in crystals on biologically relevant timescales. However, practical limitations of these experiments demands parallel development of effective molecular dynamics approaches to accelerate progress and extract meaning. Here, we establish robust and accurate methods for simulating dynamics in protein crystals, a nontrivial process requiring careful attention to equilibration, environmental composition, and choice of force fields. With more than seven milliseconds of sampling of a single chain, we identify critical factors controlling agreement between simulation and experiments and show that simulated motions recapitulate ligand-induced conformational changes. This work enables a virtuous cycle between simulation and experiments for visualizing and understanding the basic functional motions of proteins.

The switch states of the GDP-bound HRAS affected by point mutations: a study from Gaussian accelerated molecular dynamics simulations and free energy landscapes

Point mutations play a vital role in the conformational transformation of HRAS. In this work, Gaussian accelerated molecular dynamics (GaMD) simulations followed by constructions of free energy landscapes (FELs) were adopted to explore the effect of mutations D33K, A59T and L120A on conformation states of the GDP-bound HRAS. The results from the post-processing analyses on GaMD trajectories suggest that mutations alter the flexibility and motion modes of the switch domains from HRAS. The analyses from FELs show that mutations induce more disordered states of the switch domains and affect interactions of GDP with HRAS, implying that mutations yield a vital effect on the binding of HRAS to effectors. The GDP-residue interaction network revealed by our current work indicates that salt bridges and hydrogen bonding interactions (HBIs) play key roles in the binding of GDP to HRAS. Furthermore, instability in the interactions of magnesium ions and GDP with the switch SI leads to the extreme disorder of the switch domains. This study is expected to provide the energetic basis and molecular mechanism for further understanding the function of HRAS.Communicated by Ramaswamy H. Sarma.

Structural dynamics and functional cooperativity of human NQO1 by ambient temperature serial crystallography and simulations

The human NQO1 (hNQO1) is a flavin adenine nucleotide (FAD)-dependent oxidoreductase that catalyzes the two-electron reduction of quinones to hydroquinones, being essential for the antioxidant defense system, stabilization of tumor suppressors, and activation of quinone-based chemotherapeutics. Moreover, it is overexpressed in several tumors, which makes it an attractive cancer drug target. To decipher new structural insights into the flavin reductive half-reaction of the catalytic mechanism of hNQO1, we have carried serial crystallography experiments at new ID29 beamline of the ESRF to determine, to the best of our knowledge, the first structure of the hNQO1 in complex with NADH. We have also performed molecular dynamics simulations of free hNQO1 and in complex with NADH. This is the first structural evidence that the hNQO1 functional cooperativity is driven by structural communication between the active sites through long-range propagation of cooperative effects across the hNQO1 structure. Both structural results and MD simulations have supported that the binding of NADH significantly decreases protein dynamics and stabilizes hNQO1 especially at the dimer core and interface. Altogether, these results pave the way for future time-resolved studies, both at x-ray free-electron lasers and synchrotrons, of the dynamics of hNQO1 upon binding to NADH as well as during the FAD cofactor reductive half-reaction. This knowledge will allow us to reveal unprecedented structural information of the relevance of the dynamics during the catalytic function of hNQO1.

FLEXR GUI: a graphical user interface for multi-conformer modeling of proteins

Proteins are well known 'shapeshifters' which change conformation to function. In crystallography, multiple conformational states are often present within the crystal and the resulting electron-density map. Yet, explicitly incorporating alternative states into models to disentangle multi-conformer ensembles is challenging. We previously reported the tool FLEXR, which, within a few minutes, automatically separates conformational signal from noise and builds the corresponding, often missing, structural features into a multi-conformer model. To make the method widely accessible for routine multi-conformer building as part of the computational toolkit for macromolecular crystallography, we present a graphical user interface (GUI) for FLEXR, designed as a plugin for Coot 1. The GUI implementation seamlessly connects FLEXR models with the existing suite of validation and modeling tools available in Coot. We envision that FLEXR will aid crystallographers by increasing access to a multi-conformer modeling method that will ultimately lead to a better representation of protein conformational heterogeneity in the Protein Data Bank. In turn, deeper insights into the protein conformational landscape may inform biology or provide new opportunities for ligand design. The code is open source and freely available on GitHub at https://github.com/TheFischerLab/FLEXR-GUI.

Changes in an enzyme ensemble during catalysis observed by high-resolution XFEL crystallography

Enzymes populate ensembles of structures necessary for catalysis that are difficult to experimentally characterize. We use time-resolved mix-and-inject serial crystallography at an x-ray free electron laser to observe catalysis in a designed mutant isocyanide hydratase (ICH) enzyme that enhances sampling of important minor conformations. The active site exists in a mixture of conformations, and formation of the thioimidate intermediate selects for catalytically competent substates. The influence of cysteine ionization on the ICH ensemble is validated by determining structures of the enzyme at multiple pH values. Large molecular dynamics simulations in crystallo and time-resolved electron density maps show that Asp17 ionizes during catalysis and causes conformational changes that propagate across the dimer, permitting water to enter the active site for intermediate hydrolysis. ICH exhibits a tight coupling between ionization of active site residues and catalysis-activated protein motions, exemplifying a mechanism of electrostatic control of enzyme dynamics.

Functional Protein Dynamics in a Crystal

Proteins are molecular machines and to understand how they work, we need to understand how they move. New pump-probe time-resolved X-ray diffraction methods open up ways to initiate and observe protein motions with atomistic detail in crystals on biologically relevant timescales. However, practical limitations of these experiments demands parallel development of effective molecular dynamics approaches to accelerate progress and extract meaning. Here, we establish robust and accurate methods for simulating dynamics in protein crystals, a nontrivial process requiring careful attention to equilibration, environmental composition, and choice of force fields. With more than seven milliseconds of sampling of a single chain, we identify critical factors controlling agreement between simulation and experiments and show that simulated motions recapitulate ligand-induced conformational changes. This work enables a virtuous cycle between simulation and experiments for visualizing and understanding the basic functional motions of proteins.

A redox switch allows binding of Fe(II) and Fe(III) ions in the cyanobacterial iron-binding protein FutA from Prochlorococcus

The marine cyanobacterium Prochlorococcus is a main contributor to global photosynthesis, whilst being limited by iron availability. Cyanobacterial genomes generally encode two different types of FutA iron-binding proteins: periplasmic FutA2 ABC transporter subunits bind Fe(III), while cytosolic FutA1 binds Fe(II). Owing to their small size and their economized genome Prochlorococcus ecotypes typically possess a single futA gene. How the encoded FutA protein might bind different Fe oxidation states was previously unknown. Here, we use structural biology techniques at room temperature to probe the dynamic behavior of FutA. Neutron diffraction confirmed four negatively charged tyrosinates, that together with a neutral water molecule coordinate iron in trigonal bipyramidal geometry. Positioning of the positively charged Arg103 side chain in the second coordination shell yields an overall charge-neutral Fe(III) binding state in structures determined by neutron diffraction and serial femtosecond crystallography. Conventional rotation X-ray crystallography using a home source revealed X-ray-induced photoreduction of the iron center with observation of the Fe(II) binding state; here, an additional positioning of the Arg203 side chain in the second coordination shell maintained an overall charge neutral Fe(II) binding site. Dose series using serial synchrotron crystallography and an XFEL X-ray pump-probe approach capture the transition between Fe(III) and Fe(II) states, revealing how Arg203 operates as a switch to accommodate the different iron oxidation states. This switching ability of the Prochlorococcus FutA protein may reflect ecological adaptation by genome streamlining and loss of specialized FutA proteins.

Droplet microfluidics for time-resolved serial crystallography

Serial crystallography requires large numbers of microcrystals and robust strategies to rapidly apply substrates to initiate reactions in time-resolved studies. Here, we report the use of droplet miniaturization for the controlled production of uniform crystals, providing an avenue for controlled substrate addition and synchronous reaction initiation. The approach was evaluated using two enzymatic systems, yielding 3 µm crystals of lysozyme and 2 µm crystals of Pdx1, an Arabidopsis enzyme involved in vitamin B6 biosynthesis. A seeding strategy was used to overcome the improbability of Pdx1 nucleation occurring with diminishing droplet volumes. Convection within droplets was exploited for rapid crystal mixing with ligands. Mixing times of <2 ms were achieved. Droplet microfluidics for crystal size engineering and rapid micromixing can be utilized to advance time-resolved serial crystallography.

Charge neutralization and β-elimination cleavage mechanism of family 42 L-rhamnose-α-1,4-D-glucuronate lyase revealed using neutron crystallography

Gum arabic (GA) is widely used as an emulsion stabilizer and edible coating and consists of a complex carbohydrate moiety with a rhamnosyl-glucuronate group capping the non-reducing ends. Enzymes that can specifically cleave the glycosidic chains of GA and modify their properties are valuable for structural analysis and industrial application. Cryogenic X-ray crystal structure of GA-specific L-rhamnose-α-1,4-D-glucuronate lyase from Fusarium oxysporum (FoRham1), belonging to the polysaccharide lyase (PL) family 42, has been previously reported. To determine the specific reaction mechanism based on its hydrogen-containing enzyme structure, we performed joint X-ray/neutron crystallography of FoRham1. Large crystals were grown in the presence of L-rhamnose (a reaction product), and neutron and X-ray diffraction datasets were collected at room temperature at 1.80 and 1.25 Å resolutions, respectively. The active site contained L-rhamnose and acetate, the latter being a partial analog of glucuronate. Incomplete H/D exchange between Arg166 and acetate suggested that a strong salt-bridge interaction was maintained. Doubly deuterated His105 and deuterated Tyr150 supported the interaction between Arg166 and the acetate. The unique hydrogen-rich environment functions as a charge neutralizer for glucuronate and stabilizes the oxyanion intermediate. The NE2 atom of His85 was deprotonated and formed a hydrogen bond with the deuterated O1 hydroxy of L-rhamnose, indicating the function of His85 as the base/acid catalyst for bond cleavage via β-elimination. Asp83 functions as a pivot between the two catalytic histidine residues by bridging them. This His-His-Asp structural motif is conserved in the PL 24, 25, and 42 families.

ATP-dependent conformational dynamics in a photoactivated adenylate cyclase revealed by fluorescence spectroscopy and small-angle X-ray scattering

Structural insights into the photoactivated adenylate cyclases can be used to develop new ways of controlling cellular cyclic adenosine monophosphate (cAMP) levels for optogenetic and other applications. In this work, we use an integrative approach that combines biophysical and structural biology methods to provide insight on the interaction of adenosine triphosphate (ATP) with the dark-adapted state of the photoactivated adenylate cyclase from the cyanobacterium Oscillatoria acuminata (OaPAC). A moderate affinity of the nucleotide for the enzyme was calculated and the thermodynamic parameters of the interaction have been obtained. Stopped-flow fluorescence spectroscopy and small-angle solution scattering have revealed significant conformational changes in the enzyme, presumably in the adenylate cyclase (AC) domain during the allosteric mechanism of ATP binding to OaPAC with small and large-scale movements observed to the best of our knowledge for the first time in the enzyme in solution upon ATP binding. These results are in line with previously reported drastic conformational changes taking place in several class III AC domains upon nucleotide binding.

Insights into Modeling Approaches in Chemistry: Assessing Ligand-Protein Binding Thermodynamics Based on Rigid-Flexible Model Molecules

In the field of chemistry, model compounds find extensive use for investigating complex objects. One prime example of such object is the protein-ligand supramolecular interaction. Prediction the enthalpic and entropic contribution to the free energy associated with this process, as well as the structural and dynamic characteristics of protein-ligand complexes poses considerable challenges. This review exemplifies modeling approaches used to study protein-ligand binding (PLB) thermodynamics by employing pairs of conformationally constrained/flexible model molecules. Strategically designing the model molecules can reduce the number of variables that influence thermodynamic parameters. This enables scientists to gain deeper insights into the enthalpy and entropy of PLB, which is relevant for medicinal chemistry and drug design. The model studies reviewed here demonstrate that rigidifying ligands may induce compensating changes in the enthalpy and entropy of binding. Some "rules of thumb" have started to emerge on how to minimize entropy-enthalpy compensation and design efficient rigidified or flexible ligands.

From femtoseconds to minutes: time-resolved macromolecular crystallography at XFELs and synchrotrons

Over the last decade, the development of time-resolved serial crystallography (TR-SX) at X-ray free-electron lasers (XFELs) and synchrotrons has allowed researchers to study phenomena occurring in proteins on the femtosecond-to-minute timescale, taking advantage of many technical and methodological breakthroughs. Protein crystals of various sizes are presented to the X-ray beam in either a static or a moving medium. Photoactive proteins were naturally the initial systems to be studied in TR-SX experiments using pump-probe schemes, where the pump is a pulse of visible light. Other reaction initiations through small-molecule diffusion are gaining momentum. Here, selected examples of XFEL and synchrotron time-resolved crystallography studies will be used to highlight the specificities of the various instruments and methods with respect to time resolution, and are compared with cryo-trapping studies.

Pushed to extremes: distinct effects of high temperature versus pressure on the structure of STEP

Protein function hinges on small shifts of three-dimensional structure. Elevating temperature or pressure may provide experimentally accessible insights into such shifts, but the effects of these distinct perturbations on protein structures have not been compared in atomic detail. To quantitatively explore these two axes, we report the first pair of structures at physiological temperature versus. high pressure for the same protein, STEP (PTPN5). We show that these perturbations have distinct and surprising effects on protein volume, patterns of ordered solvent, and local backbone and side-chain conformations. This includes interactions between key catalytic loops only at physiological temperature, and a distinct conformational ensemble for another active-site loop only at high pressure. Strikingly, in torsional space, physiological temperature shifts STEP toward previously reported active-like states, while high pressure shifts it toward a previously uncharted region. Altogether, our work indicates that temperature and pressure are complementary, powerful, fundamental macromolecular perturbations.

X-ray Crystallography Module in MD Simulation Program Amber 2023. Refining the Models of Protein Crystals

The MD simulation package Amber offers an attractive platform to refine crystallographic structures of proteins: (i) state-of-the-art force fields help to regularize protein coordinates and reconstruct the poorly diffracting elements of the structure, such as flexible loops; (ii) MD simulations restrained by the experimental diffraction data provide an effective strategy to optimize structural models of protein crystals, including explicitly modeled interstitial solvent as well as crystal contacts. Here, we present the new crystallography module xray, released as a part of the Amber 2023 package. This module contains functions to calculate and scale structure factors (including the contributions from bulk solvent), evaluate the maximum-likelihood-type crystallographic potential, and compute its derivative forces. The X-ray functionality of Amber no longer relies on external dependencies so that the full advantage of GPU acceleration can be taken. This makes it possible to refine in a short time hundreds of crystal models, including supercell models comprised of multiple unit cells. The new automated Amber-based refinement procedure leads to an appreciable improvement in Rfree (in some cases, by as much as 0.067) as well as MolProbity scores.

An expanded view of ligandability in the allosteric enzyme PTP1B from computational reanalysis of large-scale crystallographic data

The recent advent of crystallographic small-molecule fragment screening presents the opportunity to obtain unprecedented numbers of ligand-bound protein crystal structures from a single high-throughput experiment, mapping ligandability across protein surfaces and identifying useful chemical footholds for structure-based drug design. However, due to the low binding affinities of most fragments, detecting bound fragments from crystallographic datasets has been a challenge. Here we report a trove of 65 new fragment hits across 59 new liganded crystal structures for PTP1B, an "undruggable" therapeutic target enzyme for diabetes and cancer. These structures were obtained from computational analysis of data from a large crystallographic screen, demonstrating the power of this approach to elucidate many (~50% more) "hidden" ligand-bound states of proteins. Our new structures include a fragment hit found in a novel binding site in PTP1B with a unique location relative to the active site, one that validates another new binding site recently identified by simulations, one that links adjacent allosteric sites, and, perhaps most strikingly, a fragment that induces long-range allosteric protein conformational responses via a previously unreported intramolecular conduit. Altogether, our research highlights the utility of computational analysis of crystallographic data, makes publicly available dozens of new ligand-bound structures of a high-value drug target, and identifies novel aspects of ligandability and allostery in PTP1B.

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. The isomorphous difference map remains the gold standard for detecting structural differences between datasets. Isomorphous difference 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 minute changes in unit cell properties can render isomorphous difference maps useless. This is unnecessary. Here we describe a generalized procedure for calculating observed difference maps that retains the high sensitivity to conformational change and avoids structure refinement of the perturbed state. We have implemented this procedure in an open-source python package, MatchMaps, that can be run in any software environment supporting PHENIX and CCP4. Through examples, we 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.

AlphaFold predictions are valuable hypotheses and accelerate but do not replace experimental structure determination

Artificial intelligence-based protein structure prediction methods such as AlphaFold have revolutionized structural biology. The accuracies of these predictions vary, however, and they do not take into account ligands, covalent modifications or other environmental factors. Here, we evaluate how well AlphaFold predictions can be expected to describe the structure of a protein by comparing predictions directly with experimental crystallographic maps. In many cases, AlphaFold predictions matched experimental maps remarkably closely. In other cases, even very high-confidence predictions differed from experimental maps on a global scale through distortion and domain orientation, and on a local scale in backbone and side-chain conformation. We suggest considering AlphaFold predictions as exceptionally useful hypotheses. We further suggest that it is important to consider the confidence in prediction when interpreting AlphaFold predictions and to carry out experimental structure determination to verify structural details, particularly those that involve interactions not included in the prediction.

Advanced manufacturing provides tailor-made solutions for crystallography with x-ray free-electron lasers

Serial crystallography at large facilities, such as x-ray free-electron lasers and synchrotrons, evolved as a powerful method for the high-resolution structural investigation of proteins that are critical for human health, thus advancing drug discovery and novel therapies. However, a critical barrier to successful serial crystallography experiments lies in the efficient handling of the protein microcrystals and solutions at microscales. Microfluidics are the obvious approach for any high-throughput, nano-to-microliter sample handling, that also requires design flexibility and rapid prototyping to deal with the variable shapes, sizes, and density of crystals. Here, we discuss recent advances in polymer 3D printing for microfluidics-based serial crystallography research and present a demonstration of emerging, large-scale, nano-3D printing approaches leading into the future of 3D sample environment and delivery device fabrication from liquid jet gas-dynamic virtual nozzles devices to fixed-target sample environment technology.

High-resolution double vision of the allosteric phosphatase PTP1B

Protein tyrosine phosphatase 1B (PTP1B) plays important roles in cellular homeostasis and is a highly validated therapeutic target for multiple human ailments, including diabetes, obesity and breast cancer. However, much remains to be learned about how conformational changes may convey information through the structure of PTP1B to enable allosteric regulation by ligands or functional responses to mutations. High-resolution X-ray crystallography can offer unique windows into protein conformational ensembles, but comparison of even high-resolution structures is often complicated by differences between data sets, including non-isomorphism. Here, the highest resolution crystal structure of apo wild-type (WT) PTP1B to date is presented out of a total of ∼350 PTP1B structures in the PDB. This structure is in a crystal form that is rare for PTP1B, with two unique copies of the protein that exhibit distinct patterns of conformational heterogeneity, allowing a controlled comparison of local disorder across the two chains within the same asymmetric unit. The conformational differences between these chains are interrogated in the apo structure and between several recently reported high-resolution ligand-bound structures. Electron-density maps in a high-resolution structure of a recently reported activating double mutant are also examined, and unmodeled alternate conformations in the mutant structure are discovered that coincide with regions of enhanced conformational heterogeneity in the new WT structure. These results validate the notion that these mutations operate by enhancing local dynamics, and suggest a latent susceptibility to such changes in the WT enzyme. Together, these new data and analysis provide a detailed view of the conformational ensemble of PTP1B and highlight the utility of high-resolution crystallography for elucidating conformational heterogeneity with potential relevance for function.

An RNA excited conformational state at atomic resolution

Sparse and short-lived excited RNA conformational states are essential players in cell physiology, disease, and therapeutic development, yet determining their 3D structures remains challenging. Combining mutagenesis, NMR spectroscopy, and computational modeling, we determined the 3D structural ensemble formed by a short-lived (lifetime ~2.1 ms) lowly-populated (~0.4%) conformational state in HIV-1 TAR RNA. Through a strand register shift, the excited conformational state completely remodels the 3D structure of the ground state (RMSD from the ground state = 7.2 ± 0.9 Å), forming a surprisingly more ordered conformational ensemble rich in non-canonical mismatches. The structure impedes the formation of the motifs recognized by Tat and the super elongation complex, explaining why this alternative TAR conformation cannot activate HIV-1 transcription. The ability to determine the 3D structures of fleeting RNA states using the presented methodology holds great promise for our understanding of RNA biology, disease mechanisms, and the development of RNA-targeting therapeutics.

Representing structures of the multiple conformational states of proteins

Biomolecules exhibit dynamic behavior that single-state models of their structures cannot fully capture. We review some recent advances for investigating multiple conformations of biomolecules, including experimental methods, molecular dynamics simulations, and machine learning. We also address the challenges associated with representing single- and multiple-state models in data archives, with a particular focus on NMR structures. Establishing standardized representations and annotations will facilitate effective communication and understanding of these complex models to the broader scientific community.

Structural analysis of wild-type and Val120Thr mutant Candida boidinii formate dehydrogenase by X-ray crystallography

Candida boidinii NAD+-dependent formate dehydrogenase (CbFDH) has gained significant attention for its potential application in the production of biofuels and various industrial chemicals from inorganic carbon dioxide. The present study reports the atomic X-ray crystal structures of wild-type CbFDH at cryogenic and ambient temperatures, as well as that of the Val120Thr mutant at cryogenic temperature, determined at the Turkish Light Source `Turkish DeLight'. The structures reveal new hydrogen bonds between Thr120 and water molecules in the active site of the mutant CbFDH, suggesting increased stability of the active site and more efficient electron transfer during the reaction. Further experimental data is needed to test these hypotheses. Collectively, these findings provide invaluable insights into future protein-engineering efforts that could potentially enhance the efficiency and effectiveness of CbFDH.

Kinetic Resolution of the Atomic 3D Structures Formed by Ground and Excited Conformational States in an RNA Dynamic Ensemble

Knowing the 3D structures formed by the various conformations populating the RNA free-energy landscape, their relative abundance, and kinetic interconversion rates is required to obtain a quantitative and predictive understanding of how RNAs fold and function at the atomic level. While methods integrating ensemble-averaged experimental data with computational modeling are helping define the most abundant conformations in RNA ensembles, elucidating their kinetic rates of interconversion and determining the 3D structures of sparsely populated short-lived RNA excited conformational states (ESs) remains challenging. Here, we developed an approach integrating Rosetta-FARFAR RNA structure prediction with NMR residual dipolar couplings and relaxation dispersion that simultaneously determines the 3D structures formed by the ground-state (GS) and ES subensembles, their relative abundance, and kinetic rates of interconversion. The approach is demonstrated on HIV-1 TAR, whose six-nucleotide apical loop was previously shown to form a sparsely populated (∼13%) short-lived (lifetime ∼ 45 μs) ES. In the GS, the apical loop forms a broad distribution of open conformations interconverting on the pico-to-nanosecond time scale. Most residues are unpaired and preorganized to bind the Tat-superelongation protein complex. The apical loop zips up in the ES, forming a narrow distribution of closed conformations, which sequester critical residues required for protein recognition. Our work introduces an approach for determining the 3D ensemble models formed by sparsely populated RNA conformational states, provides a rare atomic view of an RNA ES, and kinetically resolves the atomic 3D structures of RNA conformational substates, interchanging on time scales spanning 6 orders of magnitude, from picoseconds to microseconds.

Radiation damage to biological macromolecules∗

In this review, we describe recent research developments into radiation damage effects in macromolecular X-ray crystallography observed at synchrotrons and X-ray free electron lasers. Radiation damage in small molecule X-ray crystallography, small angle X-ray scattering experiments, microelectron diffraction, and single particle cryo-electron microscopy is briefly covered.

Grasping cryptic binding sites to neutralize drug resistance in the field of anticancer

Drug resistance is a significant obstacle to successful cancer treatment. The utilization and development of cryptic binding sites (CBSs) in proteins involved in cancer-related drug-resistance (CRDR) could help to overcome that drug resistance. However, there is no comprehensive review of the successful use of CBSs in addressing CRDR. Here, we have systematically summarized and analyzed the opportunities and challenges of using CBSs in addressing CRDR and revealed the key role that CBSs have in targeting CRDR. First, we have identified the CRDR targets and the corresponding CBSs. Second, we discuss the mechanisms by which CBSs can overcome CRDR. Finally, we have provided examples of successful CBS applications in addressing CRDR. We hope that this approach will provide guidance to biologists and chemists in effectively utilizing CBSs for the development of new drugs to alleviate CRDR.

Introduction to diffuse scattering and data collection

A long-standing goal in X-ray crystallography has been to extract information about the collective motions of proteins from diffuse scattering: the weak, textured signal that is found in the background of diffraction images. In the past few years, the field of macromolecular diffuse scattering has seen dramatic progress, and many of the past challenges in measurement and interpretation are now considered tractable. However, the concept of diffuse scattering is still new to many researchers, and a general set of procedures needed to collect a high-quality dataset has never been described in detail. Here, we provide the first guidelines for performing diffuse scattering experiments, which can be performed at any macromolecular crystallography beamline that supports room-temperature studies with a direct detector. We begin with a brief introduction to the theory of diffuse scattering and then walk the reader through the decision-making processes involved in preparing for and conducting a successful diffuse scattering experiment. Finally, we define quality metrics and describe ways to assess data quality both at the beamline and at home. Data obtained in this way can be processed independently by crystallographic software and diffuse scattering software to produce both a crystal structure, which represents the average atomic coordinates, and a three-dimensional diffuse scattering map that can then be interpreted in terms of models for protein motions.

Molecular-dynamics simulations of macromolecular diffraction, part II: Analysis of protein crystal simulations

Molecular-dynamics (MD) simulations have contributed substantially to our understanding of protein structure and dynamics, yielding insights into many biological processes including protein folding, drug binding, and mechanisms of protein-protein interactions. Much of what is known about protein structure comes from macromolecular crystallography (MX) experiments. MD simulations of protein crystals are useful in the study of MX because the simulations can be analyzed to calculate almost any crystallographic observable of interest, from atomic coordinates to structure factors and densities, B-factors, multiple conformations and their populations/occupancies, and diffuse scattering intensities. Computing resources and software to support crystalline MD simulations are now readily available to many researchers studying protein structure and dynamics and who may be interested in advanced interpretation of MX data, including diffuse scattering. In this work, we outline methods of analyzing MD simulations of protein crystals and provide accompanying Jupyter notebooks as practical resources for researchers wishing to perform similar analyses on their own systems of interest.

Molecular-dynamics simulations of macromolecular diffraction, part I: Preparation of protein crystal simulations

Molecular-dynamics (MD) simulations of protein crystals enable the prediction of structural and dynamical features of both the protein and the solvent components of macromolecular crystals, which can be validated against diffraction data from X-ray crystallographic experiments. The simulations have been useful for studying and predicting both Bragg and diffuse scattering in protein crystallography; however, the preparation is not yet automated and includes choices and tradeoffs that can impact the results. Here we examine some of the intricacies and consequences of the choices involved in setting up MD simulations of protein crystals for the study of diffraction data, and provide a recipe for preparing the simulations, packaged in an accompanying Jupyter notebook. This article and the accompanying notebook are intended to serve as practical resources for researchers wishing to put these models to work.

Changes in an Enzyme Ensemble During Catalysis Observed by High Resolution XFEL Crystallography

Enzymes populate ensembles of structures with intrinsically different catalytic proficiencies that are difficult to experimentally characterize. We use time-resolved mix-and-inject serial crystallography (MISC) at an X-ray free electron laser (XFEL) to observe catalysis in a designed mutant (G150T) isocyanide hydratase (ICH) enzyme that enhances sampling of important minor conformations. The active site exists in a mixture of conformations and formation of the thioimidate catalytic intermediate selects for catalytically competent substates. A prior proposal for active site cysteine charge-coupled conformational changes in ICH is validated by determining structures of the enzyme over a range of pH values. A combination of large molecular dynamics simulations of the enzyme in crystallo and time-resolved electron density maps shows that ionization of the general acid Asp17 during catalysis causes additional conformational changes that propagate across the dimer interface, connecting the two active sites. These ionization-linked changes in the ICH conformational ensemble permit water to enter the active site in a location that is poised for intermediate hydrolysis. ICH exhibits a tight coupling between ionization of active site residues and catalysis-activated protein motions, exemplifying a mechanism of electrostatic control of enzyme dynamics.

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.