Hydrogen-deuterium exchange reveals long-range dynamical allostery in soybean lipoxygenase

Transitioning enzyme catalysis towards photocatalysis

Enzyme biocatalysis is being industrialized at a phenomenal rate. Biocatalysis offers routes to chemical transformations that avoid the use of expensive metal catalysts, high temperatures and pressures, while providing impressive enantio-, regio- and chemo-selectivities. Working individually or as cascades, in live cells or cell-free preparations, to manufacture everyday chemicals, materials, healthcare products, fuels and pharmaceuticals and in diagnostic and industrial sensing applications, enzymes are key enablers in a circular bioeconomy. An ability to exploit and tailor biocatalysts rapidly and predictably requires knowledge of structure-mechanism relationships and the physical chemistry of enzyme action. This knowledge has advanced since our millennium article on this topic (Sutcliffe and Scrutton Phil Trans R. Soc. Lond. A. 2000. 358, 367-386). Here, we discuss an emerging frontier-enzyme photobiocatalysis. Photoenzymes are rarely found in nature. This limits 'difficult-to-achieve' reactions in biology that are generally accessible to chemical photocatalysts. We discuss here the emergence of photobiocatalysis as a new frontier. We review knowledge of natural photoenzymes and identify challenges and limitations in their use as photobiocatalysts. We consider emerging reports on repurposing natural enzymes as photobiocatalysts. We also discuss prospects for de novo design of photobiocatalysts which as a general concept would transform catalysis science.This article is part of the theme issue 'Science into the next millennium: 25 years on'.

A Foundational Shift in Models for Enzyme Function

This Perspective addresses the unresolved, and still hotly contested, question of how enzymes transition from stable enzyme-substrate (ES) complexes to successful, femtosecond barrier crossings. By extending Marcus theory to enzyme-catalyzed reactions, we argue that environmental reorganization of the protein scaffold, together with associated water molecules, achieves the intersection of reactant and product potential energy surfaces. After discussing the experimentally demonstrated importance of reduced activation enthalpy in enzyme-catalyzed transformations, we describe new methodologies that measure the temperature dependence of (i) time-averaged hydrogen/deuterium exchange into backbone amides and (ii) time-dependent Stokes shifts to longer emission wavelengths in appended chromophores at the protein/water interface. These methods not only identify specific pathways for the transfer of thermal energy from solvent to the reacting bonds of bound substrates but also suggest that collective thermally activated protein restructuring must occur very rapidly (on the ns-ps time scale) over long distances. Based on these findings, we introduce a comprehensive model for how barrier crossing takes place from the ES complex. This exploits the structural preorganization inherent in protein folding and subsequent conformational sampling, which optimally positions essential catalytic components within ES ground states and correctly places reactive bonds in the substrate(s) relative to embedded energy transfer networks connecting the protein surface to the active site. The existence of these anisotropic energy distribution pathways introduces a new dimension into the ongoing quest for improved de novo enzyme design.

Identification of Scaffold Specific Energy Transfer Networks in the Enthalpic Activation of Orotidine 5'-Monophosphate Decarboxylase

Orotidine 5'-monophosphate decarboxylase (OMPDC) is one of the most efficient enzyme systems studied, enhancing the decarboxylation of OMP to uridine 5'-monophosphate (UMP) by ca. 17 orders of magnitude, primarily by reducing the enthalpy of activation by ca. 28 kcal/mol. Despite a substantial reduction in activation enthalpy, OMPDC requires 15 kcal/mol of activation energy post-ES complex formation. This study investigates the physical basis of how thermal energy from solvent collisions is directed into the active site of enzyme to enable efficient thermal activation of the reaction. Comparative study of temperature-dependent hydrogen-deuterium exchange mass spectrometry (TDHDX) for WT and mutant forms of enzymes has recently been shown to uncover site specific protein networks for thermal energy transfer from solvent to enzyme active sites. In this study, we interrogate region-specific changes in the enthalpic barrier for local protein flexibility using a native OMPDC from Methanothermobacter thermautotrophicus (Mt-OMPDC) and a single site variant (Leu123Ala) that alters the activation enthalpy for catalytic turnover. The data obtained implicate four spatially resolved, thermally sensitive networks that originate at different protein/solvent interfaces and terminate at sites surrounding the substrate near the substrate phosphate-binding region (R203), the substrate- ribose binding region (K42), and a reaction enhancing loop5 (S127). These are proposed to act synergistically, transiently optimizing the position and electrostatics of the reactive carboxylate of the substrate to facilitate activated complex formation. The uncovered complexity of thermal activation networks in Mt-OMPDC distinguishes this enzyme from other members of the TIM barrel family previously investigated by TDHDX. The new findings extend the essential role of protein scaffold dynamics in orchestrating enzyme activity, with broad implications for the design of highly efficient biocatalysts.

3D variability analysis reveals a hidden conformational change controlling ammonia transport in human asparagine synthetase

Advances in X-ray crystallography and cryogenic electron microscopy (cryo-EM) offer the promise of elucidating functionally relevant conformational changes that are not easily studied by other biophysical methods. Here we show that 3D variability analysis (3DVA) of the cryo-EM map for wild-type (WT) human asparagine synthetase (ASNS) identifies a functional role for the Arg-142 side chain and test this hypothesis experimentally by characterizing the R142I variant in which Arg-142 is replaced by isoleucine. Support for Arg-142 playing a role in the intramolecular translocation of ammonia between the active site of the enzyme is provided by the glutamine-dependent synthetase activity of the R142 variant relative to WT ASNS, and MD simulations provide a possible molecular mechanism for these findings. Combining 3DVA with MD simulations is a generally applicable approach to generate testable hypotheses of how conformational changes in buried side chains might regulate function in enzymes.

3D Variability Analysis Reveals a Hidden Conformational Change Controlling Ammonia Transport in Human Asparagine Synthetase

How motions in enzymes might be linked to catalytic function is of considerable general interest. Advances in X-ray crystallography and cryogenic electron microscopy (cryo-EM) offer the promise of elucidating functionally relevant conformational changes that are not easily studied by other biophysical methods. Here we use 3D variability analysis (3DVA) of the cryo-EM map for wild-type (WT) human asparagine synthetase (ASNS) to identify a functional role for the Arg-142 side chain as a gate that mediates ammonia access to a catalytically relevant intramolecular tunnel. Our 3DVA-derived hypothesis is assessed experimentally, using the R142I variant in which Arg-142 is replaced by isoleucine, and by molecular dynamics (MD) simulations on independent, computational models of the WT human ASNS monomer and its catalytically relevant, ternary complex with β-aspartyl-AMP and MgPPi. Residue fluctuations in the MD trajectories for the human ASNS monomer are consistent with those determined for 3DVA-derived structures. These MD simulations also indicate that the gating function of Arg-142 is separate from the molecular events that form a continuous tunnel linking the two active sites. Experimental support for Arg-142 playing a role in intramolecular ammonia translocation is provided by the glutamine-dependent synthetase activity of the R142 variant relative to WT ASNS. MD simulations of computational models for the R142I variant and the R142I/β-aspartyl-AMP/MgPPi ternary complex provide a possible molecular basis for this observation. Overall, the combination of 3DVA with MD simulations is a generally applicable approach to generate testable hypotheses of how conformational changes in buried side chains might regulate function in enzymes.

Impact of N-Glycosylation on Protein Structure and Dynamics Linked to Enzymatic C-H Activation in the M. oryzae Lipoxygenase

Lipoxygenases (LOXs) from pathogenic fungi are potential therapeutic targets for defense against plant and select human diseases. In contrast to the canonical LOXs in plants and animals, fungal LOXs are unique in having appended N-linked glycans. Such important post-translational modifications (PTMs) endow proteins with altered structure, stability, and/or function. In this study, we present the structural and functional outcomes of removing or altering these surface carbohydrates on the LOX from the devastating rice blast fungus, M. oryzae, MoLOX. Alteration of the PTMs did notinfluence the active site enzyme-substrate ground state structures as visualized by electron-nuclear double resonance (ENDOR) spectroscopy. However, removal of the eight N-linked glycans by asparagine-to-glutamine mutagenesis nonetheless led to a change in substrate selectivity and an elevated activation energy for the reaction with substrate linoleic acid, as determined by kinetic measurements. Comparative hydrogen-deuterium exchange mass spectrometry (HDX-MS) analysis of wild-type and Asn-to-Gln MoLOX variants revealed a regionally defined impact on the dynamics of the arched helix that covers the active site. Guided by these HDX results, a single glycan sequon knockout was generated at position 72, and its comparative substrate selectivity from kinetics nearly matched that of the Asn-to-Gln variant. The cumulative data from model glyco-enzyme MoLOX showcase how the presence, alteration, or removal of even a single N-linked glycan can influence the structural integrity and dynamics of the protein that are linked to an enzyme's catalytic proficiency, while indicating that extensive glycosylation protects the enzyme during pathogenesis by protecting it from protease degradation.

Identification of the Thermal Activation Network in Human 15-Lipoxygenase-2: Divergence from Plant Orthologs and Its Relationship to Hydrogen Tunneling Activation Barriers

The oxidation of polyunsaturated fatty acids by lipoxygenases (LOXs) is initiated by a C-H cleavage step in which the hydrogen atom is transferred quantum mechanically (i.e., via tunneling). In these reactions, protein thermal motions facilitate the conversion of ground-state enzyme-substrate complexes to tunneling-ready configurations and are thus important for transferring energy from the solvent to the active site for the activation of catalysis. In this report, we employed temperature-dependent hydrogen-deuterium exchange mass spectrometry (TDHDX-MS) to identify catalytically linked, thermally activated peptides in a representative animal LOX, human epithelial 15-LOX-2. TDHDX-MS of wild-type 15-LOX-2 was compared to two active site mutations that retain structural stability but have increased activation energies (Ea) of catalysis. The Ea value of one variant, V427L, is implicated to arise from suboptimal substrate positioning by increased active-site side chain rotamer dynamics, as determined by X-ray crystallography and ensemble refinement. The resolved thermal network from the comparative Eas of TDHDX-MS between wild-type and V426A is localized along the front face of the 15-LOX-2 catalytic domain. The network contains a clustering of isoleucine, leucine, and valine side chains within the helical peptides. This thermal network of 15-LOX-2 is different in location, area, and backbone structure compared to a model plant lipoxygenase from soybean that exhibits a low Ea value of catalysis compared to the human ortholog. The presented data provide insights into the divergence of thermally activated protein motions in plant and animal LOXs and their relationships to the enthalpic barriers for facilitating hydrogen tunneling.

Biochemical and hydrogen-deuterium exchange studies of the single nucleotide polymorphism Y649C in human platelet 12-lipoxygenase linked to a bleeding disorder

Human platelet 12-lipoxygenase (h12-LOX) is responsible for the formation of oxylipin products that play an important role in platelet aggregation. Single nucleotide polymorphisms (SNPs) of h12-LOX have been implicated in several diseases. In this study, we investigate the structural, dynamical, and functional impact of a h12-LOX SNP that generates a tyrosine-to-cysteine mutation at a buried site (Y649C h12-LOX) and was previously ascribed with reduced levels of 12(S)-hydroxyeicosatetraenoic acid (12S-HETE) production in isolated platelets. Herein, in vitro Michaelis-Menten kinetics show reduced catalytic rates for Y649C compared to WT h12-LOX at physiological or lower temperatures. Both proteins exhibited similar melting temperatures, metal content, and oligomerization state. Liposome binding for both proteins was also dependent upon the presence of calcium, temperature, and liposome composition; however, the Y649C variant was found to have lowered binding capacity to liposomes compared to WT at physiological temperatures. Further, hydrogen-deuterium exchange mass spectrometry (HDX-MS) experiments revealed a regional defined enhancement in the peptide mobility caused by the mutation. This increased instability for the mutation stemmed from a change in an interaction with an arched helix that lines the substrate binding site, located ≥15 Å from the mutation site. Finally, differential scanning calorimetry demonstrated a reduced protein (un)folding enthalpy, consistent with the HDX results. Taken together, these results demonstrate remarkable similarity between the mutant and WT h12-LOX, and yet, subtle changes in activity, membrane affinity and protein stability may be responsible for the significant physiological changes that the Y649C SNP manifests in platelet biology.

Helical remodeling augments 5-lipoxygenase activity in the synthesis of pro-inflammatory mediators

The synthesis of pro-inflammatory leukotrienes implicated in asthma, allergic rhinitis, and atherosclerosis is initiated by the enzyme 5-lipoxygenase (5-LOX). The crystal structure of human Stable-5-LOX revealed a conformation where the catalytic iron was inaccessible to bulk solvent as two aromatic residues on a conserved helix-α2 (Hα2) plugged the substrate access portal. Whether 5-LOX can also adopt a more open conformation has not been resolved. Here, we present a new conformation of 5-LOX where Hα2 adopts an elongated conformation equivalent to that described in other animal lipoxygenase structures. Our observation of the sigmoidal kinetic behavior of 5-LOX, which is indicative of positive cooperativity, is consistent with a substrate-induced conformational change that shifts the ensemble of enzyme populations to favor the catalytically competent state. Strategic point mutations along Hα2 designed to unlock the closed conformation and elongate Hα2 resulted in improved kinetic parameters, altered limited-proteolysis data, and a drastic reduction in the length of the lag phase yielding the most active Stable-5-LOX to date. Structural predictions by AlphaFold2 of these variants statistically favor an elongated Hα2 and reinforce a model in which improved kinetic parameters correlate with a more readily adopted, open conformation. Taken together, these data provide valuable insights into the synthesis of leukotrienes.

Integration of machine learning with computational structural biology of plants

Computational structural biology of proteins has developed rapidly in recent decades with the development of new computational tools and the advancement of computing hardware. However, while these techniques have widely been used to make advancements in human medicine, these methods have seen less utilization in the plant sciences. In the last several years, machine learning methods have gained popularity in computational structural biology. These methods have enabled the development of new tools which are able to address the major challenges that have hampered the wide adoption of the computational structural biology of plants. This perspective examines the remaining challenges in computational structural biology and how the development of machine learning techniques enables more in-depth computational structural biology of plants.

Thermodynamic and biophysical study of fatty acid effector binding to soybean lipoxygenase: implications for allostery driven by helix α2 dynamics

Previous comparative kinetic isotope effects have inferred an allosteric site for fatty acids and their derivatives that modulates substrate selectivity in 15-lipoxygenases. Hydrogen-deuterium exchange also previously revealed regionally defined enhanced protein flexibility, centred at helix α2 - a gate to the substrate entrance. Direct evidence for allosteric binding and a complete understanding of its mechanism remains elusive. In this study, we examine the binding thermodynamics of the fatty acid mimic, oleyl sulfate (OS), with the monomeric model plant 15-LOX, soybean lipoxygenase (SLO), using isothermal titration calorimetry. Dynamic light scattering and differential scanning calorimetry rule out OS-induced oligomerization or structural changes. These data provide evidence that the fatty acid allosteric regulation of SLO is controlled by the dynamics of helix α2.

Kinetic and structural investigations of novel inhibitors of human epithelial 15-lipoxygenase-2

Human epithelial 15-lipoxygenase-2 (h15-LOX-2, ALOX15B) is expressed in many tissues and has been implicated in atherosclerosis, cystic fibrosis and ferroptosis. However, there are few reported potent/selective inhibitors that are active ex vivo. In the current work, we report newly discovered molecules that are more potent and structurally distinct from our previous inhibitors, MLS000545091 and MLS000536924 (Jameson et al, PLoS One, 2014, 9, e104094), in that they contain a central imidazole ring, which is substituted at the 1-position with a phenyl moiety and with a benzylthio moiety at the 2-position. The initial three molecules were mixed-type, non-reductive inhibitors, with IC₅₀ values of 0.34 ± 0.05 μM for MLS000327069, 0.53 ± 0.04 μM for MLS000327186 and 0.87 ± 0.06 μM for MLS000327206 and greater than 50-fold selectivity versus h5-LOX, h12-LOX, h15-LOX-1, COX-1 and COX-2. A small set of focused analogs was synthesized to demonstrate the validity of the hits. In addition, a binding model was developed for the three imidazole inhibitors based on computational docking and a co-structure of h15-LOX-2 with MLS000536924. Hydrogen/deuterium exchange (HDX) results indicate a similar binding mode between MLS000536924 and MLS000327069, however, the latter restricts protein motion of helix-α2 more, consistent with its greater potency. Given these results, we designed, docked, and synthesized novel inhibitors of the imidazole scaffold and confirmed our binding mode hypothesis. Importantly, four of the five inhibitors mentioned above are active in an h15-LOX-2/HEK293 cell assay and thus they could be important tool compounds in gaining a better understanding of h15-LOX-2's role in human biology. As such, a suite of similar pharmacophores that target h15-LOX-2 both in vitro and ex vivo are presented in the hope of developing them as therapeutic agents.

Advances in Hydrogen/Deuterium Exchange Mass Spectrometry and the Pursuit of Challenging Biological Systems

Solution-phase hydrogen/deuterium exchange (HDX) coupled to mass spectrometry (MS) is a widespread tool for structural analysis across academia and the biopharmaceutical industry. By monitoring the exchangeability of backbone amide protons, HDX-MS can reveal information about higher-order structure and dynamics throughout a protein, can track protein folding pathways, map interaction sites, and assess conformational states of protein samples. The combination of the versatility of the hydrogen/deuterium exchange reaction with the sensitivity of mass spectrometry has enabled the study of extremely challenging protein systems, some of which cannot be suitably studied using other techniques. Improvements over the past three decades have continually increased throughput, robustness, and expanded the limits of what is feasible for HDX-MS investigations. To provide an overview for researchers seeking to utilize and derive the most from HDX-MS for protein structural analysis, we summarize the fundamental principles, basic methodology, strengths and weaknesses, and the established applications of HDX-MS while highlighting new developments and applications.

Fatty Acid Allosteric Regulation of C-H Activation in Plant and Animal Lipoxygenases

Lipoxygenases (LOXs) catalyze the (per) oxidation of fatty acids that serve as important mediators for cell signaling and inflammation. These reactions are initiated by a C-H activation step that is allosterically regulated in plant and animal enzymes. LOXs from higher eukaryotes are equipped with an N-terminal PLAT (Polycystin-1, Lipoxygenase, Alpha-Toxin) domain that has been implicated to bind to small molecule allosteric effectors, which in turn modulate substrate specificity and the rate-limiting steps of catalysis. Herein, the kinetic and structural evidence that describes the allosteric regulation of plant and animal lipoxygenase chemistry by fatty acids and their derivatives are summarized.

In Vitro Analytical Approaches to Study Plant Ligand-Receptor Interactions

State-of-the-art in vitro methods characterize receptor-ligand interactions, highlighting experiment strategies, advantages and limitations.

Comparative kinetic isotope effects on first- and second-order rate constants of soybean lipoxygenase variants uncover a substrate-binding network

Lipoxygenases are widespread enzymes found in virtually all eukaryotes, including fungi, and, more recently, in prokaryotes. These enzymes act on long-chain polyunsaturated fatty acid substrates (C18 to C20), raising questions regarding how the substrate threads its way from solvent to the active site. Herein, we report a comparison of the temperature dependence of isotope effects on first- and second-order rate constants among single-site variants of the prototypic plant enzyme soybean lipoxygenase-1 substituted at amino acid residues inferred to impact substrate binding. We created 10 protein variants including four amino acid positions, Val-750, Ile-552, Ile-839, and Trp-500, located within a previously proposed substrate portal. The conversion of these bulky hydrophobic side chains to smaller side chains is concluded to increase the mobility of flanking helices, giving rise to increased off rates for substrate dissociation from the enzyme. In this manner, we identified a specific "binding network" that can regulate movement of the substrate from the solvent to the active site. Taken together with our previous findings on C-H and O2 activation of soybean lipoxygenase-1, these results support the emergence of multiple complementary networks within a single protein scaffold that modulate different steps along the enzymatic reaction coordinate.

Kinetic Characterization of the C-H Activation Step for the Lipoxygenase from the Pathogenic Fungus Magnaporthe oryzae: Impact of N-Linked Glycosylation

Lipoxygenases from pathogenic fungi belong to the lipoxygenase family of enzymes, which catalyze C-H activation of polyunsaturated fatty acids to form a diverse set of cell-signaling hydroperoxides. While the lipoxygenase catalytic domains are structurally and functionally similar, these fungal enzymes are decorated with N-linked glycans. The impact of N-linked glycans on the structure and function of these enzymes remains largely unknown. One exemplary system is MoLOX, a lipoxygenase from the fungus Magnaporthe oryzae, that is emerging as an important target for the devastating rice blast disease. Herein, we demonstrate that hydrogen transfer, associated with C-H cleavage of the substrate linoleic acid by MoLOX, is rate-determining and occurs by a hydrogen tunneling mechanism. Using the differential enthalpic barrier for hydrogen and deuterium transfer, ΔE_a, as a kinetic reporter of tunneling efficiency, a disproportionate increase in the activation energy for deuterium transfer is observed upon treatment of MoLOX with a peptide:N-glycosidase that cleaves N-linked carbohydrates from the protein. This increased ΔE_a is consistent with an impairment of substrate positioning in the enzyme-substrate complex for both the tunneling ready state and the ground state. These results provide new insight into the functional consequences of N-linked glycosylation on lipoxygenase C-H activation and have important implications for MoLOX inhibitor design.

HDX-MS reveals orthosteric and allosteric changes in apolipoprotein-D structural dynamics upon binding of progesterone

Apolipoprotein-D is a glycosylated tetrameric lipocalin that binds and transports small hydrophobic molecules such as progesterone and arachidonic acid. Like other lipocalins, apolipoprotein-D adopts an eight-stranded β-barrel fold stabilized by two intramolecular disulphide bonds, with an adjacent α-helix. Crystallography studies of recombinant apolipoprotein-D demonstrated no major conformational changes upon progesterone binding. Amide hydrogen-deuterium exchange mass spectrometry (HDX-MS) reports structural changes of proteins in solution by monitoring exchange of amide hydrogens in the protein backbone with deuterium. HDX-MS detects changes in conformation and structural dynamics in response to protein function such as ligand binding that may go undetected in X-ray crystallography, making HDX-MS an invaluable orthogonal technique. Here, we report an HDX-MS protocol for apolipoprotein-D that solved challenges of high protein rigidity and low pepsin cleavage using rigorous quenching conditions and longer deuteration times, yielding 85% sequence coverage and 50% deuterium exchange. The relative fractional deuterium exchange of ligand-free apolipoprotein-D revealed apolipoprotein-D to be a highly structured protein. Progesterone binding was detected by significant reduction in deuterium exchange in eight peptides. Stabilization of apolipoprotein-D dynamics can be interpreted as a combined orthosteric effect in the ligand binding pocket and allosteric effect at the N-terminus and C-terminus. Together, our experiments provide insight into apolipoprotein-D structural dynamics and map the effects of progesterone binding that are relayed to distal parts of the protein. The observed stabilization of apolipoprotein-D dynamics upon progesterone binding demonstrates a common behaviour in the lipocalin family and may have implications for interactions of apolipoprotein-D with receptors or lipoprotein particles. Statement: We reveal for the first time how apolipoprotein-D, which is protective in Alzheimer's disease, becomes more ordered when bound to a molecule of steroid hormone. These results significantly extend the understanding of apolipoprotein-D structure from X-ray crystallography studies by incorporating information on how protein motion changes over time. To achieve these results an improved protocol was developed, suitable for proteins similar to apolipoprotein-D, to elucidate how proteins change flexibility when binding to small molecules.

Defining a Canonical Ligand-Binding Pocket in the Orphan Nuclear Receptor Nurr1

Nuclear receptor-related 1 protein (Nurr1/NR4A2) is an orphan nuclear receptor (NR) that is considered to function without a canonical ligand-binding pocket (LBP). A crystal structure of the Nurr1 ligand-binding domain (LBD) revealed no physical space in the conserved region where other NRs with solvent accessible apo-protein LBPs bind synthetic and natural ligands. Using solution nuclear magnetic resonance spectroscopy, hydrogen/deuterium exchange mass spectrometry, and molecular dynamics simulations, we show that the putative canonical Nurr1 LBP is dynamic with high solvent accessibility, exchanges between two or more conformations on the microsecond-to-millisecond timescale, and can expand from the collapsed crystallized conformation to allow binding of unsaturated fatty acids. These findings should stimulate future studies to probe the ligandability and druggability of Nurr1 for both endogenous and synthetic ligands, which could lead to new therapeutics for Nurr1-related diseases, including Parkinson's disease and schizophrenia.

Role of Computational Methods in Going beyond X-ray Crystallography to Explore Protein Structure and Dynamics

Protein structural biology came a long way since the determination of the first three-dimensional structure of myoglobin about six decades ago. Across this period, X-ray crystallography was the most important experimental method for gaining atomic-resolution insight into protein structures. However, as the role of dynamics gained importance in the function of proteins, the limitations of X-ray crystallography in not being able to capture dynamics came to the forefront. Computational methods proved to be immensely successful in understanding protein dynamics in solution, and they continue to improve in terms of both the scale and the types of systems that can be studied. In this review, we briefly discuss the limitations of X-ray crystallography in studying protein dynamics, and then provide an overview of different computational methods that are instrumental in understanding the dynamics of proteins and biomacromolecular complexes.

Mass Spectral Detection of Forward- and Reverse-Hydrogen/Deuterium Exchange Resulting from Residual Solvent Vapors in Electrospray Sources

Characterizing glycans is analytically challenging since glycans are heterogeneous, branched polymers with different three-dimensional conformations. Hydrogen/deuterium exchange-mass spectrometry (HDX-MS) has been used to analyze native conformations and dynamics of biomolecules by measuring the mass increase of analytes as labile protons are replaced with deuterium following exposure to deuterated solvents. The rate of exchange is dependent on the chemical functional group, the presence of hydrogen bonds, pH, temperature, charge, and solvent accessibility. HDX-MS of carbohydrates is challenging due to the rapid exchange rate of hydroxyls. Here, we describe an observed HDX reaction associated with residual solvent vapors saturating electrospray sources. When undeuterated melezitose was infused after infusing D₂O, samples with up to 73% deuterium exchange were detected. This residual solvent HDX was observed for both carbohydrates and peptides in multiple instruments and dependent on sample infusion rate, infusion time, and deuterium content of the solvent. This residual solvent HDX was observed over several minutes of sample analysis and persisted long enough to alter the measured deuterium labeling and possibly change the interpretation of the results. This work illustrates that residual solvent HDX competes with in-solution HDX for rapidly exchanging functional groups. Thus, we propose conditions to minimize this effect, specifically for top-down, in-electrospray ionization, and quench-flow HDX experiments. Graphical Abstract ᅟ.