Simulated Isotope Exchange Patterns Enable Protein Structure Determination

Integration of Hydrogen-Deuterium Exchange Mass Spectrometry with Molecular Dynamics Simulations and Ensemble Reweighting Enables High Resolution Protein-Ligand Modeling

Hydrogen-Deuterium exchange mass spectrometry's (HDX-MS) utility in identifying and characterizing protein-small molecule interaction sites has been established. The regions that are seen to be protected from exchange upon ligand binding indicate regions that may be interacting with the ligand, giving a qualitative understanding of the ligand binding pocket. However, quantitatively deriving an accurate high-resolution structure of the protein-ligand complex from the HDX-MS data remains a challenge, often limiting its use in applications such as small molecule drug design. Recent efforts have focused on the development of methods to quantitatively model Hydrogen-Deuterium exchange (HDX) data from computationally modeled structures to garner atomic level insights from peptide-level resolution HDX-MS. One such method, HDX ensemble reweighting (HDXer), employs maximum entropy reweighting of simulated HDX data to experimental HDX-MS to model structural ensembles. In this study, we implement and validate a workflow which quantitatively leverages HDX-MS data to accurately model protein-small molecule ligand interactions. To that end, we employ a strategy combining computational protein-ligand docking, molecular dynamics simulations, HDXer, and dimensional reduction and clustering approaches to extract high-resolution drug binding poses that most accurately conform with HDX-MS data. We apply this workflow to model the interaction of ERK2 and FosA with small molecule compounds and inhibitors they are known to bind. In five out of six of the protein-ligand pairs tested, the HDX derived protein-ligand complexes result in a ligand root-mean-square deviation (RMSD) within 2.5 Å of the known crystal structure ligand.

Protein complex prediction using Rosetta, AlphaFold, and mass spectrometry covalent labeling

Covalent labeling (CL) in combination with mass spectrometry can be used as an analytical tool to study and determine structural properties of protein-protein complexes. However, data from these experiments is sparse and does not unambiguously elucidate protein structure. Thus, computational algorithms are needed to deduce structure from the CL data. In this work, we present a hybrid method that combines models of protein complex subunits generated with AlphaFold with differential CL data via a CL-guided protein-protein docking in Rosetta. In a benchmark set, the RMSD (root-mean-square deviation) of the best-scoring models was below 3.6 Å for 5/5 complexes with inclusion of CL data, whereas the same quality was only achieved for 1/5 complexes without CL data. This study suggests that our integrated approach can successfully use data obtained from CL experiments to distinguish between nativelike and non-nativelike models.

Computational Structure Prediction for Antibody-Antigen Complexes From Hydrogen-Deuterium Exchange Mass Spectrometry: Challenges and Outlook

Although computational structure prediction has had great successes in recent years, it regularly fails to predict the interactions of large protein complexes with residue-level accuracy, or even the correct orientation of the protein partners. The performance of computational docking can be notably enhanced by incorporating experimental data from structural biology techniques. A rapid method to probe protein-protein interactions is hydrogen-deuterium exchange mass spectrometry (HDX-MS). HDX-MS has been increasingly used for epitope-mapping of antibodies (Abs) to their respective antigens (Ags) in the past few years. In this paper, we review the current state of HDX-MS in studying protein interactions, specifically Ab-Ag interactions, and how it has been used to inform computational structure prediction calculations. Particularly, we address the limitations of HDX-MS in epitope mapping and techniques and protocols applied to overcome these barriers. Furthermore, we explore computational methods that leverage HDX-MS to aid structure prediction, including the computational simulation of HDX-MS data and the combination of HDX-MS and protein docking. We point out challenges in interpreting and incorporating HDX-MS data into Ab-Ag complex docking and highlight the opportunities they provide to build towards a more optimized hybrid method, allowing for more reliable, high throughput epitope identification.

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.

Deuterium Isotope Effects on Acid-Base Equilibrium of Organic Compounds

Deuterium isotope effects on acid-base equilibrium have been investigated using a combined path integral and free-energy perturbation simulation method. To understand the origin of the linear free-energy relationship of ΔpKa=pKaD2O-pKaH2O versus pKaH2O, we examined two theoretical models for computing the deuterium isotope effects. In Model 1, only the intrinsic isotope exchange effect of the acid itself in water was included by replacing the titratable protons with deuterons. Here, the dominant contribution is due to the difference in zero-point energy between the two isotopologues. In Model 2, the medium isotope effects are considered, in which the free energy change as a result of replacing H₂O by D₂O in solute-solvent hydrogen-bonding complexes is determined. Although the average ΔpKa change from Model 1 was found to be in reasonable agreement with the experimental average result, the pKaH2O dependence of the solvent isotope effects is absent. A linear free-energy relationship is obtained by including the medium effect in Model 2, and the main factor is due to solvent isotope effects in the anion-water complexes. The present study highlights the significant roles of both the intrinsic isotope exchange effect and the medium solvent isotope effect.

Protein structure prediction using residue-resolved protection factors from hydrogen-deuterium exchange NMR

Hydrogen-deuterium exchange (HDX) measured by nuclear magnetic resonance (NMR) provides structural information for proteins relating to solvent accessibility and flexibility. While this structural information is beneficial, the data cannot be used exclusively to elucidate structures. However, the structural information provided by the HDX-NMR data can be supplemented by computational methods. In previous work, we developed an algorithm in Rosetta to predict structures using qualitative HDX-NMR data (categories of exchange rate). Here we expand on the effort, and utilize quantitative protection factors (PFs) from HDX-NMR for structure prediction. From observed correlations between PFs and solvent accessibility/flexibility measures, we present a scoring function to quantify the agreement with HDX data. Using a benchmark set of 10 proteins, an average improvement of 5.13 Å in root-mean-square deviation (RMSD) is observed for cases of inaccurate Rosetta predictions. Ultimately, seven out of 10 predictions are accurate without including HDX data, and nine out of 10 are accurate when using our PF-based HDX score.

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.

Protein Structure Prediction from NMR Hydrogen-Deuterium Exchange Data

Amide hydrogen-deuterium exchange (HDX) has long been used to determine regional flexibility and binding sites in proteins; however, the data are too sparse for full structural characterization. Experiments that measure HDX rates, such as HDX-NMR, have far higher throughput compared to structure determination via X-ray crystallography, cryo-EM, or a full suite of NMR experiments. Data from HDX-NMR experiments encode information on the protein structure, making HDX a prime candidate to be supplemented by computational algorithms for protein structure prediction. We have developed a methodology to incorporate HDX-NMR data into ab initio protein structure prediction using the Rosetta software framework to predict structures based on experimental agreement. To demonstrate the efficacy of our algorithm, we examined 38 proteins with HDX-NMR data available, comparing the predicted model with and without the incorporation of HDX data into scoring. The root-mean-square deviation (rmsd, a measure of the average atomic distance between superimposed models) of the predicted model improved by 1.42 Å on average after incorporating the HDX-NMR data into scoring. The average rmsd improvement for the proteins where the selected model rmsd changed after incorporating HDX data was 3.63 Å, including one improvement of more than 11 Å and seven proteins improving by greater than 4 Å, with 12/15 proteins improving overall. Additionally, for independent verification, two proteins that were not part of the original benchmark were scored including HDX data, with a dramatic improvement of the selected model rmsd of nearly 9 Å for one of the proteins. Moreover, we have developed a confidence metric allowing us to successfully identify near-native models in the absence of a native structure. Improvement in model selection with a strong confidence measure demonstrates that protein structure prediction with HDX-NMR is a powerful tool which can be performed with minimal additional computational strain and expense.

HDXmodeller: an online webserver for high-resolution HDX-MS with auto-validation

The extent to which proteins are protected from hydrogen deuterium exchange (HDX) provides valuable insight into their folding, dynamics and interactions. Characterised by mass spectrometry (MS), HDX benefits from negligible mass restrictions and exceptional throughput and sensitivity but at the expense of resolution. Exchange mechanisms which naturally transpire for individual residues cannot be accurately located or understood because amino acids are characterised in differently sized groups depending on the extent of proteolytic digestion. Here we report HDXmodeller, the world's first online webserver for high-resolution HDX-MS. HDXmodeller accepts low-resolution HDX-MS input data and returns high-resolution exchange rates quantified for each residue. Crucially, HDXmodeller also returns a set of unique statistics that can correctly validate exchange rate models to an accuracy of 99%. Remarkably, these statistics are derived without any prior knowledge of the individual exchange rates and facilitate unparallel user confidence and the capacity to evaluate different data optimisation strategies.

Interpretation of HDX Data by Maximum-Entropy Reweighting of Simulated Structural Ensembles

Hydrogen-deuterium exchange combined with mass spectrometry (HDX-MS) is a widely applied biophysical technique that probes the structure and dynamics of biomolecules without the need for site-directed modifications or bio-orthogonal labels. The mechanistic interpretation of HDX data, however, is often qualitative and subjective, owing to a lack of quantitative methods to rigorously translate observed deuteration levels into atomistic structural information. To help address this problem, we have developed a methodology to generate structural ensembles that faithfully reproduce HDX-MS measurements. In this approach, an ensemble of protein conformations is first generated, typically using molecular dynamics simulations. A maximum-entropy bias is then applied post hoc to the resulting ensemble such that averaged peptide-deuteration levels, as predicted by an empirical model, agree with target values within a given level of uncertainty. We evaluate this approach, referred to as HDX ensemble reweighting (HDXer), for artificial target data reflecting the two major conformational states of a binding protein. We demonstrate that the information provided by HDX-MS experiments and by the model of exchange are sufficient to recover correctly weighted structural ensembles from simulations, even when the relevant conformations are rarely observed. Degrading the information content of the target data—e.g., by reducing sequence coverage, by averaging exchange levels over longer peptide segments, or by incorporating different sources of uncertainty—reduces the structural accuracy of the reweighted ensemble but still allows for useful insights into the distinctive structural features reflected by the target data. Finally, we describe a quantitative metric to rank candidate structural ensembles according to their correspondence with target data and illustrate the use of HDXer to describe changes in the conformational ensemble of the membrane protein LeuT. In summary, HDXer is designed to facilitate objective structural interpretations of HDX-MS data and to inform experimental approaches and further developments of theoretical exchange models.

Structural mass spectrometry goes viral

Over the last 20 years, mass spectrometry (MS), with its ability to analyze small sample amounts with high speed and sensitivity, has more and more entered the field of structural virology, aiming to investigate the structure and dynamics of viral proteins as close to their native environment as possible. The use of non-perturbing labels in hydrogen-deuterium exchange MS allows for the analysis of interactions between viral proteins and host cell factors as well as their dynamic responses to the environment. Cross-linking MS, on the other hand, can analyze interactions in viral protein complexes and identify virus-host interactions in cells. Native MS allows transferring viral proteins, complexes and capsids into the gas phase and has broken boundaries to overcome size limitations, so that now even the analysis of intact virions is possible. Different MS approaches not only inform about size, stability, interactions and dynamics of virus assemblies, but also bridge the gap to other biophysical techniques, providing valuable constraints for integrative structural modeling of viral complex assemblies that are often inaccessible by single technique approaches. In this review, recent advances are highlighted, clearly showing that structural MS approaches in virology are moving towards systems biology and ever more experiments are performed on cellular level.

Interlaboratory Comparison of Hydrogen-Deuterium Exchange Mass Spectrometry Measurements of the Fab Fragment of NISTmAb

Hydrogen-deuterium exchange mass spectrometry (HDX-MS) is an established, powerful tool for investigating protein-ligand interactions, protein folding, and protein dynamics. However, HDX-MS is still an emergent tool for quality control of biopharmaceuticals and for establishing dynamic similarity between a biosimilar and an innovator therapeutic. Because industry will conduct quality control and similarity measurements over a product lifetime and in multiple locations, an understanding of HDX-MS reproducibility is critical. To determine the reproducibility of continuous-labeling, bottom-up HDX-MS measurements, the present interlaboratory comparison project evaluated deuterium uptake data from the Fab fragment of NISTmAb reference material (PDB: 5K8A ) from 15 laboratories. Laboratories reported ∼89 800 centroid measurements for 430 proteolytic peptide sequences of the Fab fragment (∼78 900 centroids), giving ∼100% coverage, and ∼10 900 centroid measurements for 77 peptide sequences of the Fc fragment. Nearly half of peptide sequences are unique to the reporting laboratory, and only two sequences are reported by all laboratories. The majority of the laboratories (87%) exhibited centroid mass laboratory repeatability precisions of ⟨ sLab⟩ ≤ (0.15 ± 0.01) Da (1σx̅). All laboratories achieved ⟨sLab⟩ ≤ 0.4 Da. For immersions of protein at THDX = (3.6 to 25) °C and for D2O exchange times of tHDX = (30 s to 4 h) the reproducibility of back-exchange corrected, deuterium uptake measurements for the 15 laboratories is σreproducibility15 Laboratories( tHDX) = (9.0 ± 0.9) % (1σ). A nine laboratory cohort that immersed samples at THDX = 25 °C exhibited reproducibility of σreproducibility25C cohort( tHDX) = (6.5 ± 0.6) % for back-exchange corrected, deuterium uptake measurements.

Weighing-up protein dynamics: the combination of native mass spectrometry and molecular dynamics simulations

Structural dynamics underpin biological function at the molecular level, yet many biophysical and structural biology approaches give only a static or averaged view of proteins. Native mass spectrometry yields spectra of the many states and interactions in the structural ensemble, but its spatial resolution is limited. Conversely, molecular dynamics simulations are innately high-resolution, but have a limited capacity for exploring all structural possibilities. The two techniques hence differ fundamentally in the information they provide, returning data that reflect different length scales and time scales, making them natural bedfellows. Here we discuss how the combination of native mass spectrometry with molecular dynamics simulations is enabling unprecedented insights into a range of biological questions by interrogating the motions of proteins, their assemblies, and interactions.

Quantitative Evaluation of Native Protein Folds and Assemblies by Hydrogen Deuterium Exchange Mass Spectrometry (HDX-MS)

Hydrogen deuterium exchange mass spectrometry (HDX-MS) has significant potential for protein structure initiatives but its relationship with protein conformations is unclear. We report on the efficacy of HDX-MS to distinguish between native and non-native proteins using a popular approach to calculate HDX protection factors (PFs) from protein structures. The ability of HDX-MS to identify native protein conformations is quantified by binary structural classification such that merits of the approach for protein modelling can be quantified and better understood. We show that highly accurate PF calculations are not a prerequisite for HDX-MS simulations that are capable of effectively discriminating between native and non-native protein folds. The simulations can also be performed directly on unique structures facilitating high-throughput evaluation of many alternate conformations. The ability of HDX-MS to classify the conformations of homo-protein assemblies is also investigated. In contrast to protein monomers, we show a significant lack of correspondence between the simulated and experimental HDX-MS data for these systems with a subsequent decrease in the ability of HDX-MS to identify native states. However, we demonstrate surprisingly high diagnostic ability of the simulated data for assemblies in which a significant proportion of the individual chains occupy protein-protein interfaces. We relate this to the number of peptides that can sample alternate subunit orientations and discuss these observations within the larger context of applying HDX-MS to evaluate protein structures. Graphical Abstract.

Structural model of human dUTPase in complex with a novel proteinaceous inhibitor

Human deoxyuridine 5'-triphosphate nucleotidohydrolase (dUTPase), essential for DNA integrity, acts as a survival factor for tumor cells and is a target for cancer chemotherapy. Here we report that the Staphylococcal repressor protein Stl_SaPIBov1 (Stl) forms strong complex with human dUTPase. Functional analysis reveals that this interaction results in significant reduction of both dUTPase enzymatic activity and DNA binding capability of Stl. We conducted structural studies to understand the mechanism of this mutual inhibition. Small-angle X-ray scattering (SAXS) complemented with hydrogen-deuterium exchange mass spectrometry (HDX-MS) data allowed us to obtain 3D structural models comprising a trimeric dUTPase complexed with separate Stl monomers. These models thus reveal that upon dUTPase-Stl complex formation the functional homodimer of Stl repressor dissociates, which abolishes the DNA binding ability of the protein. Active site forming dUTPase segments were directly identified to be involved in the dUTPase-Stl interaction by HDX-MS, explaining the loss of dUTPase activity upon complexation. Our results provide key novel structural insights that pave the way for further applications of the first potent proteinaceous inhibitor of human dUTPase.