Hydrogen-deuterium exchange mass spectrometry reveals folding and allostery in protein-protein interactions

HDX-MS Analysis of Catalytic Activation of IKK2 in the IκB Kinase Complex

The IκB Kinase (IKK) complex, containing catalytic IKK2 and noncatalytic NEMO subunits, plays essential roles in the induction of transcription factors of the NF-κB family. Catalytic activation of IKK2 via phosphorylation of its activation loop is promoted upon noncovalent association of linear or K63-linked polyubiquitin chains to NEMO within the IKK complex. The mechanisms of this activation remain speculative. To investigate interaction dynamics within the IKK complex during activation of IKK2, we conducted hydrogen-deuterium exchange coupled with mass spectrometry (HDX-MS) on NEMO and IKK2 proteins in their free and complex-bound states. Altered proton exchange profiles were observed in both IKK2 and NEMO upon complex formation, and changes were consistent with the involvement of distinct regions throughout the entire length of both proteins, including previously uncharacterized segments, in direct or allosteric interactions. Association with linear tetraubiquitin (Ub4) affected multiple regions of the IKK2:NEMO complex, in addition to previously identified interaction sites on NEMO. Intriguingly, observed enhanced solvent accessibility of the IKK2 activation loop within the IKK2:NEMO:Ub4 complex, coupled with contrasting protection of surrounding segments of the catalytic subunit, suggests an allosteric role for NEMO:Ub4 in priming IKK2 for phosphorylation-dependent catalytic activation.

Phosphorylation-induced flexibility of proto-oncogenic Bcl3 regulates transcriptional activation by NF-κB p52 homodimer

The overexpression of proto-oncogene Bcl3 is observed in various cancers. Bcl3 is extensively phosphorylated and associates with homodimers of NF-κB p50 and p52 to regulate transcription. Through cellular and biochemical assays, we observed that phospho-mimetic Glu substitution at Ser366 in addition to previously studied Ser33, 114 and 446 is necessary to switch Bcl3 from an IκB-like inhibitor to a transcriptional activator. To study interactive features of p52 and Bcl3, and phosphorylation- mediated changes in Bcl3 that regulate DNA-binding by p52, we performed HDX-MS of both Bcl3 and p52 within various complexes. Nature of interactions within Bcl3:(p52:p52) complex in presence and absence of DNA, differential flexibility of Bcl3, and allosteric changes in Bcl3 upon phospho-modifications revealed why a facile accommodation of DNA requires phosphorylation. The inhibitory nature of unphosphorylated Bcl3 on DNA binding by p52:p52 also relieved by a C-terminal deletion of Bcl3. Overall, this study revealed mechanistic bases of how Bcl3 phosphorylation regulates transcriptional potential of NF-κB and intricate cell physiology, a dysregulation of which can lead to cancers.

Exploring the structural acrobatics of fold-switching proteins using simplified structure-based models

Metamorphic proteins are a paradigm of the protein folding process, by encoding two or more native states, highly dissimilar in terms of their secondary, tertiary, and even quaternary structure, on a single amino acid sequence. Moreover, these proteins structurally interconvert between these native states in a reversible manner at biologically relevant timescales as a result of different environmental cues. The large-scale rearrangements experienced by these proteins, and their sometimes high mass interacting partners that trigger their metamorphosis, makes the computational and experimental study of their structural interconversion challenging. Here, we present our efforts in studying the refolding landscapes of two quintessential metamorphic proteins, RfaH and KaiB, using simplified dual-basin structure-based models (SBMs), rigorously footed on the energy landscape theory of protein folding and the principle of minimal frustration. By using coarse-grained models in which the native contacts and bonded interactions extracted from the available experimental structures of the two native states of RfaH and KaiB are merged into a single Hamiltonian, dual-basin SBM models can be generated and savvily calibrated to explore their fold-switch in a reversible manner in molecular dynamics simulations. We also describe how some of the insights offered by these simulations have driven the design of experiments and the validation of the conformational ensembles and refolding routes observed using this simple and computationally efficient models.

Hydrogen Deuterium Exchange and other Mass Spectrometry-based Approaches for Epitope Mapping

Antigen-antibody interactions are a fundamental subset of protein-protein interactions responsible for the "survival of the fittest". Determining the interacting interface of the antigen, called an epitope, and that on the antibody, called a paratope, is crucial to antibody development. Because each antigen presents multiple epitopes (unique footprints), sophisticated approaches are required to determine the target region for a given antibody. Although X-ray crystallography, Cryo-EM, and nuclear magnetic resonance can provide atomic details of an epitope, they are often laborious, poor in throughput, and insensitive. Mass spectrometry-based approaches offer rapid turnaround, intermediate structural resolution, and virtually no size limit for the antigen, making them a vital approach for epitope mapping. In this review, we describe in detail the principles of hydrogen deuterium exchange mass spectrometry in application to epitope mapping. We also show that a combination of MS-based approaches can assist or complement epitope mapping and push the limit of structural resolution to the residue level. We describe in detail the MS methods used in epitope mapping, provide our perspective about the approaches, and focus on elucidating the role that HDX-MS is playing now and in the future by organizing a discussion centered around several improvements in prototype instrument/applications used for epitope mapping. At the end, we provide a tabular summary of the current literature on HDX-MS-based epitope mapping.

SARS-CoV-2 spike protein as a bacterial lipopolysaccharide delivery system in an overzealous inflammatory cascade

Accumulating evidence indicates a potential role for bacterial lipopolysaccharide (LPS) in the overactivation of the immune response during SARS-CoV-2 infection. LPS is recognized by Toll-like receptor 4, mediating proinflammatory effects. We previously reported that LPS directly interacts with SARS-CoV-2 spike (S) protein and enhances proinflammatory activities. Using native gel electrophoresis and hydrogen-deuterium exchange mass spectrometry, we showed that LPS binds to multiple hydrophobic pockets spanning both the S1 and S2 subunits of the S protein. Molecular simulations validated by a microscale thermophoresis binding assay revealed that LPS binds to the S2 pocket with a lower affinity compared to S1, suggesting a role as an intermediate in LPS transfer. Congruently, nuclear factor-kappa B (NF-κB) activation in monocytic THP-1 cells is strongly boosted by S2. Using NF-κB reporter mice followed by bioimaging, a boosting effect was observed for both S1 and S2, with the former potentially facilitated by proteolysis. The Omicron S variant binds to LPS, but with reduced affinity and LPS boosting in vitro and in vivo. Taken together, the data provide a molecular mechanism by which S protein augments LPS-mediated hyperinflammation.

EGCG has Dual and Opposing Effects on the N-terminal Region of Self-associating α-synuclein Oligomers

Oligomers of the protein α-synuclein (α-syn) are thought to be a major toxic species in Parkinson's disease, particularly through their ability to permeabilize cell membranes. The green tea polyphenol epigallocatechin gallate (EGCG) has been found to reduce this ability. We have analyzed α-syn oligomer dynamics and interconversion by H/D exchange monitored by mass spectrometry (HDX-MS). Our results show that the two oligomers OI and OII co-exist in equilibrium; OI is a multimer of OII and its dissociation can be followed by HDX-MS by virtue of the correlated exchange of the N-terminal region. Urea destabilizes the α-syn oligomers, dissociating OI to OII and monomers. Oligomers exposed to EGCG undergo Met oxidation. Intriguingly, EGCG induces an oxidation-dependent effect on the structure of the N-terminal region. For the non-oxidized N-terminal region, EGCG increases the stability of the folded structure as measured by a higher level of protection against H/D exchange. In contrast, protection is clearly abrogated in the Met oxidized N-terminal region. Having a non-oxidized and disordered N-terminal region is known to be essential for efficient membrane binding. Therefore, our results suggest that the combined effect of a structural stabilization of the non-oxidized N-terminal region and the presence of a disordered oxidized N-terminal region renders the oligomers less cytotoxic by decreasing the ability of the N-terminal region to bind to cell membranes and facilitate their permeabilization.

Metamorphic proteins under a computational microscope: Lessons from a fold-switching RfaH protein

Metamorphic proteins constitute unexpected paradigms of the protein folding problem, as their sequences encode two alternative folds, which reversibly interconvert within biologically relevant timescales to trigger different cellular responses. Once considered a rare aberration, metamorphism may be common among proteins that must respond to rapidly changing environments, exemplified by NusG-like proteins, the only transcription factors present in every domain of life. RfaH, a specialized paralog of bacterial NusG, undergoes an all-α to all-β domain switch to activate expression of virulence and conjugation genes in many animal and plant pathogens and is the quintessential example of a metamorphic protein. The dramatic nature of RfaH structural transformation and the richness of its evolutionary history makes for an excellent model for studying how metamorphic proteins switch folds. Here, we summarize the structural and functional evidence that sparked the discovery of RfaH as a metamorphic protein, the experimental and computational approaches that enabled the description of the molecular mechanism and refolding pathways of its structural interconversion, and the ongoing efforts to find signatures and general properties to ultimately describe the protein metamorphome.

Cullin-independent recognition of HHARI substrates by a dynamic RBR catalytic domain

RING-between-RING (RBR) E3 ligases mediate ubiquitin transfer through an obligate E3-ubiquitin thioester intermediate prior to substrate ubiquitination. Although RBRs share a conserved catalytic module, substrate recruitment mechanisms remain enigmatic, and the relevant domains have yet to be identified for any member of the class. Here we characterize the interaction between the auto-inhibited RBR, HHARI (AriH1), and its target protein, 4EHP, using a combination of XL-MS, HDX-MS, NMR, and biochemical studies. The results show that (1) a di-aromatic surface on the catalytic HHARI Rcat domain forms a binding platform for substrates and (2) a phosphomimetic mutation on the auto-inhibitory Ariadne domain of HHARI promotes release and reorientation of Rcat for transthiolation and substrate modification. The findings identify a direct binding interaction between a RING-between-RING ligase and its substrate and suggest a general model for RBR substrate recognition.

Functional roles of enzyme dynamics in accelerating active site chemistry: Emerging techniques and changing concepts

With the growing acceptance of the contribution of protein conformational ensembles to enzyme catalysis and regulation, research in the field of protein dynamics has shifted toward an understanding of the atomistic properties of protein dynamical networks and the mechanisms and time scales that control such behavior. A full description of an enzymatic reaction coordinate is expected to extend beyond the active site and include site-specific networks that communicate with the protein/water interface. Advances in experimental tools for the spatial resolution of thermal activation pathways are being complemented by biophysical methods for visualizing dynamics in real time. An emerging multidimensional model integrates the impacts of bound substrate/effector on the distribution of protein substates that are in rapid equilibration near room temperature with reaction-specific protein embedded heat transfer conduits.

Allosteric couplings upon binding of RfaH to transcription elongation complexes

In every domain of life, NusG-like proteins bind to the elongating RNA polymerase (RNAP) to support processive RNA synthesis and to couple transcription to ongoing cellular processes. Structures of factor-bound transcription elongation complexes (TECs) reveal similar contacts to RNAP, consistent with a shared mechanism of action. However, NusG homologs differ in their regulatory roles, modes of recruitment, and effects on RNA synthesis. Some of these differences could be due to conformational changes in RNAP and NusG-like proteins, which cannot be captured in static structures. Here, we employed hydrogen-deuterium exchange mass spectrometry to investigate changes in local and non-local structural dynamics of Escherichia coli NusG and its paralog RfaH, which have opposite effects on expression of xenogenes, upon binding to TEC. We found that NusG and RfaH regions that bind RNAP became solvent-protected in factor-bound TECs, whereas RNAP regions that interact with both factors showed opposite deuterium uptake changes when bound to NusG or RfaH. Additional changes far from the factor-binding site were observed only with RfaH. Our results provide insights into differences in structural dynamics exerted by NusG and RfaH during binding to TEC, which may explain their different functional outcomes and allosteric regulation of transcriptional pausing by RfaH.

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.

Human FoxP Transcription Factors as Tractable Models of the Evolution and Functional Outcomes of Three-Dimensional Domain Swapping

The association of two or more proteins to adopt a quaternary complex is one of the most widespread mechanisms by which protein function is modulated. In this scenario, three-dimensional domain swapping (3D-DS) constitutes one plausible pathway for the evolution of protein oligomerization that exploits readily available intramolecular contacts to be established in an intermolecular fashion. However, analysis of the oligomerization kinetics and thermodynamics of most extant 3D-DS proteins shows its dependence on protein unfolding, obscuring the elucidation of the emergence of 3D-DS during evolution, its occurrence under physiological conditions, and its biological relevance. Here, we describe the human FoxP subfamily of transcription factors as a feasible model to study the evolution of 3D-DS, due to their significantly faster dissociation and dimerization kinetics and lower dissociation constants in comparison to most 3D-DS models. Through the biophysical and functional characterization of FoxP proteins, relevant structural aspects highlighting the evolutionary adaptations of these proteins to enable efficient 3D-DS have been ascertained. Most biophysical studies on FoxP suggest that the dynamics of the polypeptide chain are crucial to decrease the energy barrier of 3D-DS, enabling its fast oligomerization under physiological conditions. Moreover, comparison of biophysical parameters between human FoxP proteins in the context of their minute sequence differences suggests differential evolutionary strategies to favor homoassociation and presages the possibility of heteroassociations, with direct impacts in their gene regulation function.

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 Footprinting, Conformational Dynamics, and Core Interface-Adjacent Neutralization "Hotspots" in the SARS-CoV-2 Spike Protein Receptor Binding Domain/Human ACE2 Interaction

The novel severe respiratory syndrome-like coronavirus (SARS-CoV-2) causes COVID-19 in humans and is responsible for one of the most destructive pandemics of the last century. At the root of SARS-CoV infection is the interaction between the viral spike protein and the human angiotensin converting enzyme 2 protein, which allows the virus to gain entry into host cells through endocytosis. In this work, we apply hydrogen-deuterium exchange mass spectrometry (HDX-MS) to provide a detailed view of the functional footprint and conformational dynamics associated with this interaction. Our results broadly agree with the binding interface derived from high resolution X-ray crystal structure data but also provide insights into shifts in structure and dynamics that accompany complexation, including some that occur immediately outside of the core binding interface. We propose that dampening of these "binding-site adjacent" dynamic shifts could represent a mechanism for neutralizing activity in a multitude of spike protein-targeted mAbs that have been found to specifically bind these "peripheral" sites. Our results highlight the unique capacity of HDX-MS to detect potential neutralization "hotspots" outside of the core binding interfaces defined by high resolution structural data.

A nexus of intrinsic dynamics underlies translocase priming

The cytoplasmic ATPase SecA and the membrane-embedded SecYEG channel assemble to form the Sec translocase. How this interaction primes and catalytically activates the translocase remains unclear. We show that priming exploits a nexus of intrinsic dynamics in SecA. Using atomistic simulations, smFRET, and HDX-MS, we reveal multiple dynamic islands that cross-talk with domain and quaternary motions. These dynamic elements are functionally important and conserved. Central to the nexus is a slender stem through which rotation of the preprotein clamp of SecA is biased by ATPase domain motions between open and closed clamping states. An H-bonded framework covering most of SecA enables multi-tier dynamics and conformational alterations with minimal energy input. As a result, cognate ligands select preexisting conformations and alter local dynamics to regulate catalytic activity and clamp motions. These events prime the translocase for high-affinity reception of non-folded preprotein clients. Dynamics nexuses are likely universal and essential in multi-liganded proteins.

Global Profiling of Lysine Accessibility to Evaluate Protein Structure Changes in Alzheimer's Disease

The linear sequence of amino acids in a protein folds into a 3D structure to execute protein activity and function, but it is still challenging to profile the 3D structure at the proteome scale. Here, we present a method of native protein tandem mass tag (TMT) profiling of Lys accessibility and its application to investigate structural alterations in human brain specimens of Alzheimer's disease (AD). In this method, proteins are extracted under a native condition, labeled by TMT reagents, followed by trypsin digestion and peptide analysis using two-dimensional liquid chromatography and tandem mass spectrometry (LC/LC-MS/MS). The method quantifies Lys labeling efficiency to evaluate its accessibility on the protein surface, which may be affected by protein conformations, protein modifications, and/or other molecular interactions. We systematically optimized the amount of TMT reagents, reaction time, and temperature and then analyzed protein samples under multiple conditions, including different labeling time (5 and 30 min), heat treatment, AD and normal human cases. The experiment profiled 15370 TMT-labeled peptides in 4475 proteins. As expected, the heat treatment led to extensive changes in protein conformations, with 17% of the detected proteome displaying differential labeling. Compared to the normal controls, AD brain showed different Lys accessibility of tau and RNA splicing complexes, which are the hallmarks of AD pathology, as well as proteins involved in transcription, mitochondrial, and synaptic functions. To eliminate the possibility that the observed differential Lys labeling was caused by protein level change, the whole proteome was quantified with standard TMT-LC/LC-MS/MS for normalization. Thus, this native protein TMT method enables the proteome-wide measurement of Lys accessibility, representing a straightforward strategy to explore protein structure in any biological system.

Diverse Native Ensembles Dictate the Differential Functional Responses of Nuclear Receptor Ligand-Binding Domains

Native states of folded proteins are characterized by a large ensemble of conformations whose relative populations and interconversion dynamics determine the functional output. This is more apparent in transcription factors that have evolved to be inherently sensitive to small perturbations, thus fine-tuning gene expression. To explore the extent to which such functional features are imprinted on the folding landscape of transcription factor ligand-binding domains (LBDs), we characterize paralogous LBDs of the nuclear receptor (NR) family employing an energetically detailed and ensemble-based Ising-like statistical mechanical model. We find that the native ensembles of the LBDs from glucocorticoid receptor, PPAγ, and thyroid hormone receptor display a remarkable diversity in the width of the native wells, the number and nature of partially structured states, and hence the degree of conformational order. Monte Carlo simulations employing the full state representation of the ensemble highlight that many of the functional conformations coexist in equilibrium, whose relative populations are sensitive to both temperature and the strength of ligand binding. Allosteric modulation of the degree of structure at a coregulator binding site on ligand binding is shown to arise via a redistribution of populations in the native ensembles of glucocorticoid and PPAγ LBDs. Our results illustrate how functional requirements can drive the evolution of conformationally diverse native ensembles in paralogs.

Intrinsically Disordered Regions of the DNA-Binding Domain of Human FoxP1 Facilitate Domain Swapping

Forkhead box P (FoxP) proteins are unique transcription factors that spatiotemporally regulate gene expression by tethering two chromosome loci together via functional domain-swapped dimers formed through their DNA-binding domains. Further, the differential kinetics on this dimerization mechanism underlie an intricate gene regulation network at physiological conditions. Nonetheless, poor understanding of the structural dynamics and steps of the association process impedes to link the functional domain swapping to human-associated diseases. Here, we have characterized the DNA-binding domain of human FoxP1 by integrating single-molecule Förster resonance energy transfer and hydrogen-deuterium exchange mass spectrometry data with molecular dynamics simulations. Our results confirm the formation of a previously postulated domain-swapped (DS) FoxP1 dimer in solution and reveal the presence of highly populated, heterogeneous, and locally disordered dimeric intermediates along the dimer dissociation pathway. The unique features of FoxP1 provide a glimpse of how intrinsically disordered regions can facilitate domain swapping oligomerization and other tightly regulated association mechanisms relevant in biological processes.

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.

Differential Local Stability Governs the Metamorphic Fold Switch of Bacterial Virulence Factor RfaH

RfaH, a two-domain protein from a universally conserved NusG/Spt5 family of regulators, is required for the transcription and translation of long virulence and conjugation operons in many Gram-negative bacterial pathogens. Escherichia coli RfaH action is controlled by a unique large-scale structural rearrangement triggered by recruitment to transcription elongation complexes through a specific DNA element. Upon recruitment, the C-terminal domain of RfaH refolds from an α-hairpin, which is bound to RNA polymerase binding site within the N-terminal domain, into an unbound β-barrel that interacts with the ribosome. Although structures of the autoinhibited (α-hairpin) and active (β-barrel) states and plausible refolding pathways have been reported, how this reversible switch is encoded within RfaH sequence and structure is poorly understood. Here, we combined hydrogen-deuterium exchange measurements by mass spectrometry and nuclear magnetic resonance with molecular dynamics to evaluate the differential local stability between both RfaH folds. Deuteron incorporation reveals that the tip of the C-terminal hairpin (residues 125-145) is stably folded in the autoinhibited state (∼20% deuteron incorporation), whereas the rest of this domain is highly flexible (>40% deuteron incorporation), and its flexibility only decreases in the β-folded state. Computationally predicted ΔG agree with these results by displaying similar anisotropic stability within the tip of the α-hairpin and on neighboring N-terminal domain residues. Remarkably, the β-folded state shows comparable structural flexibility than nonmetamorphic homologs. Our findings provide information critical for understanding the metamorphic behavior of RfaH and other chameleon proteins and for devising targeted strategies to combat bacterial infections.

Advances in MS Based Strategies for Probing Ligand-Target Interactions: Focus on Soft Ionization Mass Spectrometric Techniques

The non-covalent interactions between small drug molecules and disease-related proteins (ligand-target interactions) mediate various pharmacological processes in the treatment of different diseases. The development of the analytical methods to assess those interactions, including binding sites, binding energies, stoichiometry and association-dissociation constants, could assist in clarifying the mechanisms of action, precise treatment of targeted diseases as well as the targeted drug discovery. For the last decades, mass spectrometry (MS) has been recognized as a powerful tool to study the non-covalent interactions of the ligand-target complexes with the characteristics of high sensitivity, high-resolution, and high-throughput. Soft ionization mass spectrometry, especially the electrospray mass spectrometry (ESI-MS) and matrix assisted laser desorption ionization mass spectrometry (MALDI-MS), could achieve the complete transformation of the target analytes into the gas phase, and subsequent detection of the small drug molecules and disease-related protein complexes, and has exerted great advantages for studying the drug ligands-protein targets interactions, even in case of identifying active components as drug ligands from crude extracts of medicinal plants. Despite of other analytical techniques for this purpose, such as the NMR and X-ray crystallography, this review highlights the principles, research hotspots and recent applications of the soft ionization mass spectrometry and its hyphenated techniques, including hydrogen-deuterium exchange mass spectrometry (HDX-MS), chemical cross-linking mass spectrometry (CX-MS), and ion mobility spectrometry mass spectrometry (IMS-MS), in the study of the non-covalent interactions between small drug molecules and disease-related proteins.

Accurate Prediction of Amide Exchange in the Fast Limit Reveals Thrombin Allostery

Amide hydrogen/deuterium exchange mass spectrometry (HDXMS) of proteins has become extremely popular for identifying ligand-binding sites, protein-protein interactions, intrinsic disorder, and allosteric changes upon protein modification. Such phenomena are revealed when amide exchange is measured in the fast limit, that is, within a few minutes of exchange in deuterated buffer. The HDXMS data have a resolution of the length of peptides and are difficult to interpret because many different phenomena lead to changes in hydrogen/deuterium exchange. We present a quantitative analysis of accelerated molecular dynamics simulations that provides impressive agreement with peptide-length HDXMS data. Comparative analysis of thrombin and a single-point mutant reveals that the simulation analysis can distinguish the subtle differences in exchange due to mutation. In addition, the results provide a deeper understanding of the underlying changes in dynamics revealed by the HDXMS that extend far from the site of mutation.