Biological insights from hydrogen exchange mass spectrometry

A helical fulcrum in eIF2B coordinates allosteric regulation of stress signaling

The integrated stress response (ISR) enables cells to survive a variety of acute stresses, but chronic activation of the ISR underlies age-related diseases. ISR signaling downregulates translation and activates expression of stress-responsive factors that promote return to homeostasis and is initiated by inhibition of the decameric guanine nucleotide exchange factor eIF2B. Conformational and assembly transitions regulate eIF2B activity, but the allosteric mechanisms controlling these dynamic transitions and mediating the therapeutic effects of the small-molecule ISR inhibitor ISRIB are unknown. Using hydrogen-deuterium exchange-mass spectrometry and cryo-electron microscopy, we identified a central α-helix whose orientation allosterically coordinates eIF2B conformation and assembly. Biochemical and cellular signaling assays show that this 'switch-helix' controls eIF2B activity and signaling. In sum, the switch-helix acts as a fulcrum of eIF2B conformational regulation and is a highly conserved actuator of ISR signal transduction. This work uncovers a conserved allosteric mechanism and unlocks new therapeutic possibilities for ISR-linked diseases.

The conformational landscape of a serpin N-terminal subdomain facilitates folding and in-cell quality control

Many multi-domain proteins including the serpin family of serine protease inhibitors contain non-sequential domains composed of regions that are far apart in sequence. Because proteins are translated vectorially from N- to C-terminus, such domains pose a particular challenge: how to balance the conformational lability necessary to form productive interactions between early and late translated regions while avoiding aggregation. This balance is mediated by the protein sequence properties and the interactions of the folding protein with the cellular quality control machinery. For serpins, particularly α 1 -antitrypsin (AAT), mutations often lead to polymer accumulation in cells and consequent disease suggesting that the lability/aggregation balance is especially precarious. Therefore, we investigated the properties of progressively longer AAT N-terminal fragments in solution and in cells. The N-terminal subdomain, residues 1-190 (AAT190), is monomeric in solution and efficiently degraded in cells. More β -rich fragments, 1-290 and 1-323, form small oligomers in solution, but are still efficiently degraded, and even the polymerization promoting Siiyama (S53F) mutation did not significantly affect fragment degradation. In vitro, the AAT190 region is among the last regions incorporated into the final structure. Hydrogen-deuterium exchange mass spectrometry and enhanced sampling molecular dynamics simulations show that AAT190 has a broad, dynamic conformational ensemble that helps protect one particularly aggregation prone β -strand from solvent. These AAT190 dynamics result in transient exposure of sequences that are buried in folded, full-length AAT, which may provide important recognition sites for the cellular quality control machinery and facilitate degradation and, under favorable conditions, reduce the likelihood of polymerization.

Integrating Hydrogen Deuterium Exchange-Mass Spectrometry with Molecular Simulations Enables Quantification of the Conformational Populations of the Sugar Transporter XylE

A yet unresolved challenge in structural biology is to quantify the conformational states of proteins underpinning function. This challenge is particularly acute for membrane proteins owing to the difficulties in stabilizing them for in vitro studies. To address this challenge, we present an integrative strategy that combines hydrogen deuterium exchange-mass spectrometry (HDX-MS) with ensemble modeling. We benchmark our strategy on wild-type and mutant conformers of XylE, a prototypical member of the ubiquitous Major Facilitator Superfamily (MFS) of transporters. Next, we apply our strategy to quantify conformational ensembles of XylE embedded in different lipid environments. Further application of our integrative strategy to substrate-bound and inhibitor-bound ensembles allowed us to unravel protein-ligand interactions contributing to the alternating access mechanism of secondary transport in atomistic detail. Overall, our study highlights the potential of integrative HDX-MS modeling to capture, accurately quantify, and subsequently visualize co-populated states of membrane proteins in association with mutations and diverse substrates and inhibitors.

Anti-CRISPR proteins function through thermodynamic tuning and allosteric regulation of CRISPR RNA-guided surveillance complex

CRISPR RNA-guided detection and degradation of foreign DNA is a dynamic process. Viruses can interfere with this cellular defense by expressing small proteins called anti-CRISPRs. While structural models of anti-CRISPRs bound to their target complex provide static snapshots that inform mechanism, the dynamics and thermodynamics of these interactions are often overlooked. Here, we use hydrogen deuterium exchange-mass spectrometry (HDX-MS) and differential scanning fluorimetry (DSF) experiments to determine how anti-CRISPR binding impacts the conformational landscape of the type IF CRISPR RNA guided surveillance complex (Csy) upon binding of two different anti-CRISPR proteins (AcrIF9 and AcrIF2). The results demonstrate that AcrIF2 binding relies on enthalpic stabilization, whereas AcrIF9 uses an entropy driven reaction to bind the CRISPR RNA-guided surveillance complex. Collectively, this work reveals the thermodynamic basis and mechanistic versatility of anti-CRISPR-mediated immune suppression. More broadly, this work presents a striking example of how allosteric effectors are employed to regulate nucleoprotein complexes.

Computational Modeling of Molecular Structures Guided by Hydrogen-Exchange Data

Data produced by hydrogen-exchange monitoring experiments have been used in structural studies of molecules for several decades. Despite uncertainties about the structural determinants of hydrogen exchange itself, such data have successfully helped guide the structural modeling of challenging molecular systems, such as membrane proteins or large macromolecular complexes. As hydrogen-exchange monitoring provides information on the dynamics of molecules in solution, it can complement other experimental techniques in so-called integrative modeling approaches. However, hydrogen-exchange data have often only been used to qualitatively assess molecular structures produced by computational modeling tools. In this paper, we look beyond qualitative approaches and survey the various paradigms under which hydrogen-exchange data have been used to quantitatively guide the computational modeling of molecular structures. Although numerous prediction models have been proposed to link molecular structure and hydrogen exchange, none of them has been widely accepted by the structural biology community. Here, we present as many hydrogen-exchange prediction models as we could find in the literature, with the aim of providing the first exhaustive list of its kind. From purely structure-based models to so-called fractional-population models or knowledge-based models, the field is quite vast. We aspire for this paper to become a resource for practitioners to gain a broader perspective on the field and guide research toward the definition of better prediction models. This will eventually improve synergies between hydrogen-exchange monitoring and molecular modeling.

A Review on Computational Approaches for Analyzing Hydrogen- Deuterium (H/D) Exchange of Proteins

Native state Hydrogen-Deuterium (H/D) exchange method has been used to study the structures and the unfolding pathways for quite a number of proteins. The H/D exchange method is generally monitored using nuclear magnetic resonance (NMR) spectroscopy and mass spectrometry (MS) techniques. NMR-assisted H/D exchange methods primarily monitor the residue level fluctuation of proteins, whereas MS-assisted H/D exchange methods analyze multifold ensemble conformations of proteins. In this connection, quite a large number of computational tools and algorithms have been developed for processing and analyzing huge amount of the H/D exchange data generated from these techniques. In this review, most of the freely available computational tools associated with the H/D exchange of proteins have been comprehensively reviewed and scopes to improve/ develop novel computational approaches for analyzing the H/D exchange data of proteins have also been brought into fore.

Characterizing ligand-induced conformational changes in clinically relevant galectin-1 by HN/H2O (D2O) exchange

Glycans of cellular glycoconjugates serve as biochemical signals for a multitude of (patho)physiological processes via binding to their receptors (e.g. lectins). In the case of human adhesion/growth-regulatory galectin-1 (Gal-1), small angle neutron scattering and fluorescence correlation spectroscopy have revealed a significant decrease of its gyration radius and increase of its diffusion coefficient upon binding lactose, posing the pertinent question on the nature and region(s) involved in the underlying structural alterations. Requiring neither a neutron source nor labeling, diffusion measurements by ¹H NMR spectroscopy are shown here to be sufficiently sensitive to detect this ligand-induced change. In order to figure out which region(s) of Gal-1 is (are) affected at the level of peptides, we first explored the use of H/D exchange mass spectrometry (HDX MS). Hereby, we found a reduction in proton exchange kinetics beyond the lactose-binding site. The measurement of fast H^N/H₂O exchange by phase-modulated NMR clean chemical exchange (CLEANEX) NMR on ¹⁵N-labeled Gal-1 then increased the spatial resolution to the level of individual amino acids. The mapped regions with increased protection from H^N/H₂O (D₂O) exchange that include the reduction of solvent exposure around the interface can underlie the protein's compaction. These structural changes have potential to modulate this galectin's role in lattice formation on the cell surface and its interaction(s) with protein(s) at the F-face.

Bovine Serum Albumin Catalysed Hydrogen and Deuterium Evolution at Mercury Electrodes

The hydrogen evolution reaction (HER), catalysed by proteins at mercury electrodes and reflected in chronopotentiometric stripping peak H, provides a label-free and reagentless analytical technique that is sensitive to protein structure. Here we show how the kinetic isotope effect affected the HER catalysed by the protein bovine serum albumin (BSA). We found that the deuteron bond, which is stronger than that of a proton, contributed to less effective transport of deuterons mediated by BSA at the Hg|D₂ O interface, and enhanced structural stability of the surface-attached native BSA in D₂ O solution. A structural transition was also observed in the surface-attached urea-denatured BSA, and is probably due to the destabilisation of some secondary structural remnants retained by the 17 SS-bonds. Because the catalytically active groups involved in proton or deuteron transfer in native proteins are often exposed towards solutions and their protons exchange almost instantly, no signs of H/D exchange were observed in native BSA using peak H under the given conditions.

Assessing Site-Specific Enhancements Imparted by Hyperpolarized Water in Folded and Unfolded Proteins by 2D HMQC NMR

Hyperpolarized water can be a valuable aid in protein NMR, leading to amide group ¹H polarizations that are orders of magnitude larger than their thermal counterparts. Suitable procedures can exploit this to deliver 2D ¹H–¹⁵N correlations with good resolution and enhanced sensitivity. These enhancements depend on the exchange rates between the amides and the water, thereby yielding diagnostic information about solvent accessibility. This study applied this “HyperW” method to four proteins exhibiting a gamut of exchange behaviors: PhoA^(350–471), an unfolded 122-residue fragment; barstar, a fully folded ribonuclease inhibitor; R17, a 13.3 kDa system possessing folded and unfolded forms under slow interconversion; and drkN SH3, a protein domain whose folded and unfolded forms interchange rapidly and with temperature-dependent population ratios. For PhoA4^(350–471) HyperW sensitivity enhancements were ≥300×, as expected for an unfolded protein sequence. Though fully folded, barstar also exhibited substantial enhancements; these, however, were not uniform and, according to CLEANEX experiments, reflected the solvent-exposed residues. R17 showed the expected superposition of ≥100-fold enhancements for its unfolded form, coexisting with more modest enhancements for their folded counterparts. Unexpected, however, was the behavior of drkN SH3, for which HyperW enhanced the unfolded but, surprisingly, enhanced even more certain folded protein sites. These preferential enhancements were repeatedly and reproducibly observed. A number of explanations—including three-site exchange magnetization transfers between water and the unfolded and folded states; cross-correlated relaxation processes from hyperpolarized “structural” waters and labile side-chain protons; and the possibility that faster solvent exchange rates characterize certain folded sites over their unfolded counterparts—are considered to account for them.

Hydrogen/Deuterium Exchange Mass Spectrometry for the Structural Analysis of Detergent-Solubilized Membrane Proteins

Integral membrane proteins are involved in numerous biological functions and represent important drug targets. Despite their abundance in the human proteome, the number of integral membrane protein structures is largely underrepresented in the Protein Data Bank. The challenges associated with the biophysical characterization of such biological systems are well known. Most structural approaches, including X-ray crystallography, SAXS, or mass spectrometry (MS), require the complete solubilization of membrane proteins in aqueous solutions. Detergents are frequently used for this task, but may interfere with the analysis, as is the case with MS. The use of "MS-friendly" detergents, such as non-ionic alkyl glycoside detergents, has greatly facilitated the analysis of detergent-solubilized membrane proteins. Here, we describe a protocol, which we have successfully implemented in our laboratory to study the structure and dynamics of detergent-solubilized integral membrane proteins by Hydrogen/Deuterium eXchange and Mass Spectrometry (HDX-MS). The procedure does not require detergent removal prior to MS analysis, instead taking advantage of the ultra-high pressure chromatographic system to separate deuterated peptides from "MS-friendly" detergents.

Thermodynamics of Helix formation in small peptides of varying lengthin vacuo, implicit solvent and explicit solvent: Comparison between AMBER force fields

Helix formation is of great significance in protein folding. The helix-forming tendencies of amino acids are accumulated along the sequence to determine the helix-forming tendency of peptides. Computer simulation can be used to model this process in atomic details and give structural insights. In the current work, we employ equilibrate-state free energy simulation to systematically study the folding/unfolding thermodynamics of a series of mutated peptides. Two AMBER force fields including AMBER99SB and AMBER14SB are compared. The new 14SB force field uses refitted torsion parameters compared with 99SB and they share the same atomic charge scheme. We find that in vacuo the helix formation is mutation dependent, which reflects the different helix propensities of different amino acids. In general, there are helix formers, helix indifferent groups and helix breakers. The helical structure becomes more favored when the length of the sequence becomes longer, which arises from the formation of additional backbone hydrogen bonds in the lengthened sequence. Therefore, the helix indifferent groups and helix breakers will become helix formers in long sequences. Also, protonation-dependent helix formation is observed for ionizable groups. In 14SB, the helical structures are more stable than in 99SB and differences can be observed in their grouping schemes, especially in the helix indifferent group. In solvents, all mutations are helix indifferent due to protein–solvent interactions. The decrease in the number of backbone hydrogen bonds is the same with the increase in the number of protein–water hydrogen bonds. The 14SB in explicit solvent is able to capture the free energy minima in the helical state while 14SB in implicit solvent, 99SB in explicit solvent and 99SB in implicit solvent cannot. The helix propensities calculated under 14SB agree with the corresponding experimental values, while the 99SB results obviously deviate from the references. Hence, implicit solvent models are unable to correctly describe the thermodynamics even for the simple helix formation in isolated peptides. Well-developed force fields and explicit solvents are needed to correctly describe the protein dynamics. Aside from the free energy, differences in conformational ensemble under different force fields in different solvent models are observed. The numbers of hydrogen bonds formed under different force fields agree and they are mostly determined by the solvent model.

Probing Cascade complex composition and stability using native mass spectrometry techniques

Adaptive prokaryotic immune systems rely on clustered regularly interspaced short palindromic repeats and their associated genes to provide the components necessary to clear infection by foreign genetic elements. These immune systems are based on highly specific nucleases that bind DNA or RNA and, upon sequence recognition, degrade the bound nucleic acid. Because of their specificity, CRISPR-Cas systems are being co-opted to edit genes in eukaryotic cells. While the general function of these systems is well understood, an understanding of mechanistic details to facilitate engineering and application to this new arena remains a topic of intense study. Here, we present two methods that have been successfully used to study the structure and mechanism of the Type IE CRISPR system, Cascade, from Escherichia coli. We provide the protocol for a typical native mass spectrometry experiment which, because it allows for analysis of a protein complex without disruption of the noncovalent interactions within the complex, can be used to determine complex composition, architecture, and relative affinity between subunits. We, also, provide the protocol for intact protein hydrogen-deuterium exchange mass spectrometry, which provides insight into the overall conformational stability of the complex and changes in complex stability based on conditions such as substrate binding. Investigating the solution-phase structure, stability, and dynamics of these complexes improves the overall understanding of the mechanism facilitating engineered adjustments to function or utility.

Fine-tuning of the stability of β-strands by Y181 in perfringolysin O directs the prepore to pore transition

Perfringolysin O (PFO) is a toxic protein that forms β-barrel transmembrane pores upon binding to cholesterol-containing membranes. The formation of lytic pores requires conformational changes in PFO that lead to the conversion of water-soluble monomers into membrane-bound oligomers. Although the general outline of stepwise pore formation has been established, the underlying mechanistic details await clarification. To extend our understanding of the molecular mechanisms that control the pore formation, we compared the hydrogen-deuterium exchange patterns of PFO with its derivatives bearing mutations in the D3 domain. In the case of two of these mutations F318A, Y181A, known from previous work to lead to a decreased lytic activity, global destabilization of all protein domains was observed in their water-soluble forms. This was accompanied by local changes in D3 β-sheet, including unexpected stabilization of functionally important β1 strand in Y181A. In case of the double mutation (F318A/Y181A) that completely abolished the lytic activity, several local changes were retained, but the global destabilization effects of single mutations were reverted and hydrogen-deuterium exchange (HDX) pattern returned to PFO level. Strong structural perturbations were not observed in case of remaining variants in which other residues of the hydrophobic core of D3 domain were substituted by alanine. Our results indicate the existence in PFO of a well-tuned H-bonding network that maintains the stability of the D3 β-strands at appropriate level at each transformation step. F318 and Y181 moieties participate in this network and their role extends beyond their direct intermolecular interaction during oligomerization that was identified previously.

Revealing Unknown Protein Structures Using Computational Conformational Sampling Guided by Experimental Hydrogen-Exchange Data

Both experimental and computational methods are available to gather information about a protein's conformational space and interpret changes in protein structure. However, experimentally observing and computationally modeling large proteins remain critical challenges for structural biology. Our work aims at addressing these challenges by combining computational and experimental techniques relying on each other to overcome their respective limitations. Indeed, despite its advantages, an experimental technique such as hydrogen-exchange monitoring cannot produce structural models because of its low resolution. Additionally, the computational methods that can generate such models suffer from the curse of dimensionality when applied to large proteins. Adopting a common solution to this issue, we have recently proposed a framework in which our computational method for protein conformational sampling is biased by experimental hydrogen-exchange data. In this paper, we present our latest application of this computational framework: generating an atomic-resolution structural model for an unknown protein state. For that, starting from an available protein structure, we explore the conformational space of this protein, using hydrogen-exchange data on this unknown state as a guide. We have successfully used our computational framework to generate models for three proteins of increasing size, the biggest one undergoing large-scale conformational changes.

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.

Reference Parameters for Protein Hydrogen Exchange Rates

The analysis of many hydrogen exchange (HX) experiments depends on knowledge of exchange rates expected for the unstructured protein under the same conditions. We present here some minor adjustments to previously calibrated values and a stringent test of their accuracy. Graphical Abstract ᅟ.

Native State of Complement Protein C3d Analysed via Hydrogen Exchange and Conformational Sampling

Hydrogen/deuterium exchange detected by mass spectrometry (HDXMS) provides valuable information on protein structure and dynamics. Although HDX-MS data is often interpreted using crystal structures, it was suggested that conformational ensembles produced by molecular dynamics simulations yield more accurate interpretations. In this paper, we analyse the complement protein C3d by performing an HDX-MS experiment, and evaluate several interpretation methodologies using an existing prediction model to derive HDX-MS data from protein structure. To interpret and refine C3d's HDX-MS data, we look for a conformation (or conformational ensemble) of C3d that allows computationally replicating this data. We confirm that crystal structures are not a good choice and suggest that conformational ensembles produced by molecular dynamics simulations might not always be satisfactory either. Finally, we show that coarse-grained conformational sampling of C3d produces a conformation from which its HDX-MS data can be replicated and refined.

Hydrogen-Deuterium Exchange Mass Spectrometry to Study Protein Complexes

Hydrogen-deuterium exchange (HDX) mass spectrometry (MS) can provide valuable information about binding, allostery, and other conformational effects of interaction in protein complexes. For protein-ligand complexes, where the ligand may be a small molecule, peptide, nucleotide, or another protein(s), a typical experiment measures HDX in the protein alone and then compares that with HDX for the protein when part of the complex. Multiple factors are critical in the design and implementation of such experiments, including thoughtful consideration of the percent protein bound, the effects of the labeling protocol on the protein complex, and the dynamic range of the analysis method. With careful planning and techniques, HDX MS analysis of protein complexes can be very informative.

Disulfide driven folding for a conditionally disordered protein

Conditionally disordered proteins are either ordered or disordered depending on the environmental context. The substrates of the mitochondrial intermembrane space (IMS) oxidoreductase Mia40 are synthesized on cytosolic ribosomes and diffuse as intrinsically disordered proteins to the IMS, where they fold into their functional conformations; behaving thus as conditionally disordered proteins. It is not clear how the sequences of these polypeptides encode at the same time for their ability to adopt a folded structure and to remain unfolded. Here we characterize the disorder-to-order transition of a Mia40 substrate, the human small copper chaperone Cox17. Using an integrated real-time approach, including chromatography, fluorescence, CD, FTIR, SAXS, NMR, and MS analysis, we demonstrate that in this mitochondrial protein, the conformational switch between disordered and folded states is controlled by the formation of a single disulfide bond, both in the presence and in the absence of Mia40. We provide molecular details on how the folding of a conditionally disordered protein is tightly regulated in time and space, in such a way that the same sequence is competent for protein translocation and activity.

Disulfide driven folding for a conditionally disordered protein

Conditionally disordered proteins are either ordered or disordered depending on the environmental context. The substrates of the mitochondrial intermembrane space (IMS) oxidoreductase Mia40 are synthesized on cytosolic ribosomes and diffuse as intrinsically disordered proteins to the IMS, where they fold into their functional conformations; behaving thus as conditionally disordered proteins. It is not clear how the sequences of these polypeptides encode at the same time for their ability to adopt a folded structure and to remain unfolded. Here we characterize the disorder-to-order transition of a Mia40 substrate, the human small copper chaperone Cox17. Using an integrated real-time approach, including chromatography, fluorescence, CD, FTIR, SAXS, NMR, and MS analysis, we demonstrate that in this mitochondrial protein, the conformational switch between disordered and folded states is controlled by the formation of a single disulfide bond, both in the presence and in the absence of Mia40. We provide molecular details on how the folding of a conditionally disordered protein is tightly regulated in time and space, in such a way that the same sequence is competent for protein translocation and activity.

H/D exchange mass spectrometry and statistical coupling analysis reveal a role for allostery in a ferredoxin-dependent bifurcating transhydrogenase catalytic cycle

Recent investigations into ferredoxin-dependent transhydrogenases, a class of enzymes responsible for electron transport, have highlighted the biological importance of flavin-based electron bifurcation (FBEB). FBEB generates biomolecules with very low reduction potential by coupling the oxidation of an electron donor with intermediate potential to the reduction of high and low potential molecules. Bifurcating systems can generate biomolecules with very low reduction potentials, such as reduced ferredoxin (Fd), from species such as NADPH. Metabolic systems that use bifurcation are more efficient and confer a competitive advantage for the organisms that harbor them. Structural models are now available for two NADH-dependent ferredoxin-NADP⁺ oxidoreductase (Nfn) complexes. These models, together with spectroscopic studies, have provided considerable insight into the catalytic process of FBEB. However, much about the mechanism and regulation of these multi-subunit proteins remains unclear. Using hydrogen/deuterium exchange mass spectrometry (HDX-MS) and statistical coupling analysis (SCA), we identified specific pathways of communication within the model FBEB system, Nfn from Pyrococus furiosus, under conditions at each step of the catalytic cycle. HDX-MS revealed evidence for allosteric coupling across protein subunits upon nucleotide and ferredoxin binding. SCA uncovered a network of co-evolving residues that can provide connectivity across the complex. Together, the HDX-MS and SCA data show that protein allostery occurs across the ensemble of iron‑sulfur cofactors and ligand binding sites using specific pathways that connect domains allowing them to function as dynamically coordinated units.

Computational methods and challenges in hydrogen/deuterium exchange mass spectrometry

Hydrogen/Deuterium exchange (HDX) has been applied, since the 1930s, as an analytical tool to study the structure and dynamics of (small) biomolecules. The popularity of using HDX to study proteins increased drastically in the last two decades due to the successful combination with mass spectrometry (MS). Together with this growth in popularity, several technological advances have been made, such as improved quenching and fragmentation. As a consequence of these experimental improvements and the increased use of protein-HDXMS, large amounts of complex data are generated, which require appropriate analysis. Computational analysis of HDXMS requires several steps. A typical workflow for proteins consists of identification of (non-)deuterated peptides or fragments of the protein under study (local analysis), or identification of the deuterated protein as a whole (global analysis); determination of the deuteration level; estimation of the protection extent or exchange rates of the labile backbone amide hydrogen atoms; and a statistically sound interpretation of the estimated protection extent or exchange rates. Several algorithms, specifically designed for HDX analysis, have been proposed. They range from procedures that focus on one specific step in the analysis of HDX data to complete HDX workflow analysis tools. In this review, we provide an overview of the computational methods and discuss outstanding challenges. © 2016 Wiley Periodicals, Inc. Mass Spec Rev 36:649-667, 2017.

Hydrogen/Deuterium Exchange Mass Spectrometry of Human Green Opsin Reveals a Conserved Pro-Pro Motif in Extracellular Loop 2 of Monostable Visual G Protein-Coupled Receptors

Opsins comprise the protein component of light sensitive G protein-coupled receptors (GPCRs) in the retina of the eye that are responsible for the transduction of light into a biochemical signal. Here, we used hydrogen/deuterium (H/D) exchange coupled with mass spectrometry to map conformational changes in green cone opsin upon light activation. We then compared these findings with those reported for rhodopsin. The extent of H/D exchange in green cone opsin was greater than in rhodopsin in the dark and bleached states, suggesting a higher structural heterogeneity for green cone opsin. Further analysis revealed that green cone opsin exists as a dimer in both dark (inactive) and bleached (active) states, and that the predicted glycosylation sites at N32 and N34 are indeed glycosylated. Comparison of deuterium uptake between inactive and active states of green cone opsin also disclosed a reduced solvent accessibility of the extracellular N-terminal region and an increased accessibility of the chromophore binding site. Increased H/D exchange at the extracellular side of transmembrane helix four (TM4) combined with an analysis of sequence alignments revealed a conserved Pro-Pro motif in extracellular loop 2 (EL2) of monostable visual GPCRs. These data present new insights into the locus of chromophore release at the extracellular side of TM4 and TM5 and provide a foundation for future functional evaluation.

Coarse-Grained Conformational Sampling of Protein Structure Improves the Fit to Experimental Hydrogen-Exchange Data

Monitoring hydrogen/deuterium exchange (HDX) undergone by a protein in solution produces experimental data that translates into valuable information about the protein's structure. Data produced by HDX experiments is often interpreted using a crystal structure of the protein, when available. However, it has been shown that the correspondence between experimental HDX data and crystal structures is often not satisfactory. This creates difficulties when trying to perform a structural analysis of the HDX data. In this paper, we evaluate several strategies to obtain a conformation providing a good fit to the experimental HDX data, which is a premise of an accurate structural analysis. We show that performing molecular dynamics simulations can be inadequate to obtain such conformations, and we propose a novel methodology involving a coarse-grained conformational sampling approach instead. By extensively exploring the intrinsic flexibility of a protein with this approach, we produce a conformational ensemble from which we extract a single conformation providing a good fit to the experimental HDX data. We successfully demonstrate the applicability of our method to four small and medium-sized proteins.

R468A mutation in perfringolysin O destabilizes toxin structure and induces membrane fusion

Perfringolysin O (PFO) belongs to the family of cholesterol-dependent cytolysins. Upon binding to a cholesterol-containing membrane, PFO undergoes a series of structural changes that result in the formation of a β-barrel pore and cell lysis. Recognition and binding to cholesterol are mediated by the D4 domain, one of four domains of PFO. The D4 domain contains a conserved tryptophan-rich loop named undecapeptide (E₄₅₈CTGLAWEWWR₄₆₈) in which arginine 468 is essential for retaining allosteric coupling between D4 and other domains during interaction of PFO with the membrane. In this report we studied the impact of R468A mutation on the whole protein structure using hydrogen-deuterium exchange coupled with mass spectrometry. We found that in aqueous solution, compared to wild type (PFO), PFO^R468A showed increased deuterium uptake due to exposure of internal toxin regions to the solvent. This change reflected an overall structural destabilization of PFO^R468A in solution. Conversely, upon binding to cholesterol-containing membranes, PFO^R468A revealed a profound decrease of hydrogen-deuterium exchange when compared to PFO. This block of deuterium uptake resulted from PFO^R468A-induced aggregation and fusion of liposomes, as found by dynamic light scattering, microscopic observations and FRET measurements. In the result of liposome aggregation and fusion, the entire PFO^R468A molecule became shielded from aqueous solution and thereby was protected against proteolytic digestion and deuteration. We have established that structural changes induced by the R468A mutation lead to exposure of an additional cholesterol-independent liposome-binding site in PFO that confers its fusogenic property, altering the mode of the toxin action.

Conformational dynamics of a neurotransmitter:sodium symporter in a lipid bilayer

Neurotransmitter:sodium symporters (NSSs) are integral membrane proteins responsible for the sodium-dependent reuptake of small-molecule neurotransmitters from the synaptic cleft. The symporters for the biogenic amines serotonin (SERT), dopamine (DAT), and norepinephrine (NET) are targets of multiple psychoactive agents, and their dysfunction has been implicated in numerous neuropsychiatric ailments. LeuT, a thermostable eubacterial NSS homolog, has been exploited as a model protein for NSS members to canvass the conformational mechanism of transport with a combination of X-ray crystallography, cysteine accessibility, and solution spectroscopy. Despite yielding remarkable insights, these studies have primarily been conducted with protein in the detergent-solubilized state rather than embedded in a membrane mimic. In addition, solution spectroscopy has required site-specific labeling of nonnative cysteines, a labor-intensive process occasionally resulting in diminished transport and/or binding activity. Here, we overcome these limitations by reconstituting unlabeled LeuT in phospholipid bilayer nanodiscs, subjecting them to hydrogen-deuterium exchange coupled with mass spectrometry (HDX-MS), and facilitating interpretation of the data with molecular dynamics simulations. The data point to changes of accessibility and dynamics of structural elements previously implicated in the transport mechanism, in particular transmembrane helices (TMs) 1a and 7 as well as extracellular loops (ELs) 2 and 4. The results therefore illuminate the value of this strategy for interrogating the conformational mechanism of the more clinically significant mammalian membrane proteins including SERT and DAT, neither of which tolerates complete removal of endogenous cysteines, and whose activity is heavily influenced by neighboring lipids.

Defining Low-Dimensional Projections to Guide Protein Conformational Sampling

Exploring the conformational space of proteins is critical to characterize their functions. Numerous methods have been proposed to sample a protein's conformational space, including techniques developed in the field of robotics and known as sampling-based motion-planning algorithms (or sampling-based planners). However, these algorithms suffer from the curse of dimensionality when applied to large proteins. Many sampling-based planners attempt to mitigate this issue by keeping track of sampling density to guide conformational sampling toward unexplored regions of the conformational space. This is often done using low-dimensional projections as an indirect way to reduce the dimensionality of the exploration problem. However, how to choose an appropriate projection and how much it influences the planner's performance are still poorly understood issues. In this article, we introduce two methodologies defining low-dimensional projections that can be used by sampling-based planners for protein conformational sampling. The first method leverages information about a protein's flexibility to construct projections that can efficiently guide conformational sampling, when expert knowledge is available. The second method builds similar projections automatically, without expert intervention. We evaluate the projections produced by both methodologies on two conformational search problems involving three middle-size proteins. Our experiments demonstrate that (i) defining projections based on expert knowledge can benefit conformational sampling and (ii) automatically constructing such projections is a reasonable alternative.

MEMHDX: an interactive tool to expedite the statistical validation and visualization of large HDX-MS datasets

Motivation: With the continued improvement of requisite mass spectrometers and UHPLC systems, Hydrogen/Deuterium eXchange Mass Spectrometry (HDX-MS) workflows are rapidly evolving towards the investigation of more challenging biological systems, including large protein complexes and membrane proteins. The analysis of such extensive systems results in very large HDX-MS datasets for which specific analysis tools are required to speed up data validation and interpretation.

Results: We introduce a web application and a new R-package named ‘MEMHDX’ to help users analyze, validate and visualize large HDX-MS datasets. MEMHDX is composed of two elements. A statistical tool aids in the validation of the results by applying a mixed-effects model for each peptide, in each experimental condition, and at each time point, taking into account the time dependency of the HDX reaction and number of independent replicates. Two adjusted P-values are generated per peptide, one for the ‘Change in dynamics’ and one for the ‘Magnitude of ΔD’, and are used to classify the data by means of a ‘Logit’ representation. A user-friendly interface developed with Shiny by RStudio facilitates the use of the package. This interactive tool allows the user to easily and rapidly validate, visualize and compare the relative deuterium incorporation on the amino acid sequence and 3D structure, providing both spatial and temporal information.

Availability and Implementation: MEMHDX is freely available as a web tool at the project home page http://memhdx.c3bi.pasteur.fr

Contact:marie-agnes.dillies@pasteur.fr or sebastien.brier@pasteur.fr

Supplementary information:Supplementary data is available at Bioinformatics online.

Conformational effects in protein electrospray-ionization mass spectrometry

Electrospray-ionization mass spectrometry (ESI-MS) is a key tool of structural biology, complementing the information delivered by conventional biochemical and biophysical methods. Yet, the mechanism behind the conformational effects in protein ESI-MS is an object of debate. Two parameters-solvent-accessible surface area (As) and apparent gas-phase basicity (GBapp)-are thought to play a role in controlling the extent of protein ionization during ESI-MS experiments. This review focuses on recent experimental and theoretical investigations concerning the influence of these parameters on ESI-MS results and the structural information that can be derived. The available evidence supports a unified model for the ionization mechanism of folded and unfolded proteins. These data indicate that charge-state distribution (CSD) analysis can provide valuable structural information on normally folded, as well as disordered structures.

Analytical Aspects of Hydrogen Exchange Mass Spectrometry

This article reviews the analytical aspects of measuring hydrogen exchange by mass spectrometry (HX MS). We describe the nature of analytical selectivity in hydrogen exchange, then review the analytical tools required to accomplish fragmentation, separation, and the mass spectrometry measurements under restrictive exchange quench conditions. In contrast to analytical quantitation that relies on measurements of peak intensity or area, quantitation in HX MS depends on measuring a mass change with respect to an undeuterated or deuterated control, resulting in a value between zero and the maximum amount of deuterium that can be incorporated. Reliable quantitation is a function of experimental fidelity and to achieve high measurement reproducibility, a large number of experimental variables must be controlled during sample preparation and analysis. The method also reports on important qualitative aspects of the sample, including conformational heterogeneity and population dynamics.

The future of drug discovery: enabling technologies for enhancing lead characterization and profiling therapeutic potential

Technology often serves as a handmaiden and catalyst of invention. The discovery of safe, effective medications depends critically upon experimental approaches capable of providing high-impact information on the biological effects of drug candidates early in the discovery pipeline. This information can enable reliable lead identification, pharmacological compound differentiation and successful translation of research output into clinically useful therapeutics. The shallow preclinical profiling of candidate compounds promulgates a minimalistic understanding of their biological effects and undermines the level of value creation necessary for finding quality leads worth moving forward within the development pipeline with efficiency and prognostic reliability sufficient to help remediate the current pharma-industry productivity drought. Three specific technologies discussed herein, in addition to experimental areas intimately associated with contemporary drug discovery, appear to hold particular promise for strengthening the preclinical valuation of drug candidates by deepening lead characterization. These are: i) hydrogen-deuterium exchange mass spectrometry for characterizing structural and ligand-interaction dynamics of disease-relevant proteins; ii) activity-based chemoproteomics for profiling the functional diversity of mammalian proteomes; and iii) nuclease-mediated precision gene editing for developing more translatable cellular and in vivo models of human diseases. When applied in an informed manner congruent with the clinical understanding of disease processes, technologies such as these that span levels of biological organization can serve as valuable enablers of drug discovery and potentially contribute to reducing the current, unacceptably high rates of compound clinical failure.

Patterns of structural dynamics in RACK1 protein retained throughout evolution: a hydrogen-deuterium exchange study of three orthologs

RACK1 is a member of the WD repeat family of proteins and is involved in multiple fundamental cellular processes. An intriguing feature of RACK1 is its ability to interact with at least 80 different protein partners. Thus, the structural features enabling such interactomic flexibility are of great interest. Several previous studies of the crystal structures of RACK1 orthologs described its detailed architecture and confirmed predictions that RACK1 adopts a seven-bladed β-propeller fold. However, this did not explain its ability to bind to multiple partners. We performed hydrogen-deuterium (H-D) exchange mass spectrometry on three orthologs of RACK1 (human, yeast, and plant) to obtain insights into the dynamic properties of RACK1 in solution. All three variants retained similar patterns of deuterium uptake, with some pronounced differences that can be attributed to RACK1's divergent biological functions. In all cases, the most rigid structural elements were confined to B-C turns and, to some extent, strands B and C, while the remaining regions retained much flexibility. We also compared the average rate constants for H-D exchange in different regions of RACK1 and found that amide protons in some regions exchanged at least 1000-fold faster than in others. We conclude that its evolutionarily retained structural architecture might have allowed RACK1 to accommodate multiple molecular partners. This was exemplified by our additional analysis of yeast RACK1 dimer, which showed stabilization, as well as destabilization, of several interface regions upon dimer formation.

A new approach to measuring protein backbone protection with high spatial resolution using H/D exchange and electron capture dissociation

Inadequate spatial resolution remains one of the most serious limitations of hydrogen/deuterium exchange-mass spectrometry (HDX-MS), especially when applied to larger proteins (over 30 kDa). Supplementing proteolytic fragmentation of the protein in solution with ion dissociation in the gas phase has been used successfully by several groups to obtain near-residue level resolution. However, the restrictions imposed by the LC-MS/MS mode of operation on the data acquisition time frame makes it difficult in many cases to obtain a signal-to-noise ratio adequate for reliable assignment of the backbone amide protection levels at individual residues. This restriction is lifted in the present work by eliminating the LC separation step from the workflow and taking advantage of the high resolving power and dynamic range of a Fourier transform ion cyclotron resonance-mass spectrometer (FTICR-MS). A residue-level resolution is demonstrated for a peptic fragment of a 37 kDa recombinant protein (N-lobe of human serum transferrin), using electron-capture dissociation as an ion fragmentation tool. The absence of hydrogen scrambling in the gas phase prior to ion dissociation is verified using redundant HDX-MS data generated by FTICR-MS. The backbone protection pattern generated by direct HDX-MS/MS is in excellent agreement with the known crystal structure of the protein but also provides information on conformational dynamics, which is not available from the static X-ray structure.

Analysis of disordered proteins using a simple apparatus for millisecond quench-flow H/D exchange

Measurement of amide H/D exchange on the ms time scale can provide valuable information about the dynamic behavior of the most flexible regions of proteins. We describe here a simple mixing apparatus, assembled solely from off-the-shelf components, that can be used for H/D exchange mass spectrometry to measure exchange on the 50-5000 ms time scale. Our apparatus utilizes flow-injection to minimize sample consumption. Although the mixer operates at low Reynolds numbers (less than 10(2)) where laminar flow is expected, H/D exchange kinetics were well-approximated using the assumption of plug-flow. We validated this approximation using fluorescence imaging of fluorescein-conjugated bovine serum albumin in the delay line and by demonstrating agreement between measured and calculated H/D exchange kinetics for a mixture of peptides. The performance of the apparatus was further validated by measuring rapid H/D exchange kinetics by an intrinsically disordered protein, murine CBP(2059-2117) (UniProt CBP_MOUSE). H/D exchange data from CBP, both free and in complex with human ACTR(1018-1088) (UniProt NCOA3_HUMAN), were consistent with previous biophysical studies of this protein.