Scope and utility of hydrogen exchange as a tool for mapping landscapes

Lessons from a quarter century of being human in protein science

Over the past quarter century, my engagement with the protein society has allowed me to witness first-hand the evolution of our deepening understanding of the complexity of protein folding landscapes. During my own evolution as a protein scientist, my passion for protein folding has deepened into an obsession with mapping and decoding the thermodynamic and kinetic secrets of protein landscapes - especially those of rebel proteins, whose "non-traditional" behavior has challenged our paradigms and inspired the expansion of our models and methods. It is perhaps not surprising that I see parallels in the evolution of the landscape framework and in the development of our own trajectories as humans in STEM. Just as with proteins however, we need to recognize that our individual human landscapes are not isolated from our local departmental and institutional communities, and are integrated into the larger networks of our STEM disciplines, academia, industry and/or government, not to mention society. My experience with hundreds of participants in the Being Human in STEM initiative that Amherst College undergraduates and I co-founded in 2016 has helped me find hope for STEM and humanity. If we commit to reconciling our identities as scientists with our responsibilities as human beings, together we can accelerate the evolution of individual, community and societal landscapes to contribute to addressing the dire challenges facing our planet. This article is protected by copyright. All rights reserved.

Prediction and Validation of a Protein's Free Energy Surface Using Hydrogen Exchange and (Importantly) Its Denaturant Dependence

The denaturant dependence of hydrogen-deuterium exchange (HDX) is a powerful measurement to identify the breaking of individual H-bonds and map the free energy surface (FES) of a protein including the very rare states. Molecular dynamics (MD) can identify each partial unfolding event with atomic-level resolution. Hence, their combination provides a great opportunity to test the accuracy of simulations and to verify the interpretation of HDX data. For this comparison, we use Upside, our new and extremely fast MD package that is capable of folding proteins with an accuracy comparable to that of all-atom methods. The FESs of two naturally occurring and two designed proteins are so generated and compared to our NMR/HDX data. We find that Upside's accuracy is considerably improved upon modifying the energy function using a new machine-learning procedure that trains for proper protein behavior including realistic denatured states in addition to stable native states. The resulting increase in cooperativity is critical for replicating the HDX data and protein stability, indicating that we have properly encoded the underlying physiochemical interactions into an MD package. We did observe some mismatch, however, underscoring the ongoing challenges faced by simulations in calculating accurate FESs. Nevertheless, our ensembles can identify the properties of the fluctuations that lead to HDX, whether they be small-, medium-, or large-scale openings, and can speak to the breadth of the native ensemble that has been a matter of debate.

Exploiting the Propagation of Constrained Variables for Enhanced HDX-MS Data Optimization

Nonlinear programming has found useful applications in protein biophysics to help understand the microscopic exchange kinetics of data obtained using hydrogen-deuterium exchange mass spectrometry (HDX-MS). Finding a microscopic kinetic solution for HDX-MS data provides a window into local protein stability and energetics allowing them to be quantified and understood. Optimization of HDX-MS data is a significant challenge, however, due to the requirement to solve a large number of variables simultaneously with exceptionally large variable bounds. Modeled rates are frequently uncertain with an explicate dependency on the initial guess values. In order to enhance the search for a minimum solution in HDX-MS optimization, the ability of selected constrained variables to propagate throughout the data is considered. We reveal that locally bound constrained optimization induces a global effect on all variables. The global response to local constraints is large and surprisingly long-range, but the outcome is unpredictable, unexpectedly decreasing the overall accuracy of certain data sets depending on the stringency of the constraints. Utilizing previously described in-house validation criteria based on covariance matrices, a method is described that is able to accurately determine whether constraints benefit or impair the optimization of HDX-MS data. From this, we establish a new two-stage method for our online optimizer HDXmodeller that can effectively leverage locally bound variables to enhance HDX-MS data modeling.

How cooperative are protein folding and unfolding transitions?

A thermodynamically and kinetically simple picture of protein folding envisages only two states, native (N) and unfolded (U), separated by a single activation free energy barrier, and interconverting by cooperative two-state transitions. The folding/unfolding transitions of many proteins occur, however, in multiple discrete steps associated with the formation of intermediates, which is indicative of reduced cooperativity. Furthermore, much advancement in experimental and computational approaches has demonstrated entirely non-cooperative (gradual) transitions via a continuum of states and a multitude of small energetic barriers between the N and U states of some proteins. These findings have been instrumental towards providing a structural rationale for cooperative versus noncooperative transitions, based on the coupling between interaction networks in proteins. The cooperativity inherent in a folding/unfolding reaction appears to be context dependent, and can be tuned via experimental conditions which change the stabilities of N and U. The evolution of cooperativity in protein folding transitions is linked closely to the evolution of function as well as the aggregation propensity of the protein. A large activation energy barrier in a fully cooperative transition can provide the kinetic control required to prevent the accumulation of partially unfolded forms, which may promote aggregation. Nevertheless, increasing evidence for barrier-less "downhill" folding, as well as for continuous "uphill" unfolding transitions, indicate that gradual non-cooperative processes may be ubiquitous features on the free energy landscape of protein folding.

Protein structural dynamics at the gas/water interface examined by hydrogen exchange mass spectrometry

Gas/water interfaces (such as air bubbles or foam) are detrimental to the stability of proteins, often causing aggregation. This represents a potential problem for industrial processes, for example, the production and handling of protein drugs. Proteins possess surfactant-like properties, resulting in a high affinity for gas/water interfaces. The tendency of previously buried nonpolar residues to maximize contact with the gas phase can cause significant structural distortion. Most earlier studies in this area employed spectroscopic tools that could only provide limited information. Here we use hydrogen/deuterium exchange (HDX) mass spectrometry (MS) for probing the conformational dynamics of the model protein myoglobin (Mb) in the presence of N(2) bubbles. HDX/MS relies on the principle that unfolded and/or highly dynamic regions undergo faster deuteration than tightly folded segments. In bubble-free solution Mb displays EX2 behavior, reflecting the occurrence of short-lived excursions to partially unfolded conformers. A dramatically different behavior is seen in the presence of N(2) bubbles; EX2 dynamics still take place, but in addition the protein shows EX1 behavior. The latter results from interconversion of the native state with conformers that are globally unfolded and long-lived. These unfolded species likely correspond to Mb that is adsorbed to the surface of gas bubbles. N(2) sparging also induces aggregation. To explain the observed behavior we propose a simple model, that is, "semi-unfolded" ↔ "native" ↔ "globally unfolded" → "aggregated". This model quantitatively reproduces the experimentally observed kinetics. To the best of our knowledge, the current study marks the first exploration of surface denaturation phenomena by HDX/MS.

Conformer-specific characterization of nonnative protein states using hydrogen exchange and top-down mass spectrometry

Characterization of structure and dynamics of nonnative protein states is important for understanding molecular mechanisms of processes as diverse as folding, binding, aggregation, and enzyme catalysis to name just a few; however, selectively probing local minima within rugged energy landscapes remains a problem. Mass spectrometry (MS) coupled with hydrogen/deuterium exchange (HDX) offers a unique advantage of being able to make a distinction among multiple protein conformers that coexist in solution; however, detailed structural interrogation of such states previously remained out of reach of HDX MS. In this work, we exploited the aforementioned unique feature of HDX MS in combination with the ability of MS to isolate narrow populations of protein ions to characterize individual protein conformers coexisting in solution in equilibrium. Subsequent fragmentation of the protein ions using electron-capture dissociation allowed us to allocate the deuterium distribution along the protein backbone, yielding a backbone-amide protection map for the selected conformer unaffected by contributions from other protein states present in solution. The method was tested with the small regulatory protein ubiquitin (Ub), which is known to form nonnative intermediate states under a variety of mildly denaturing conditions. Protection maps of these intermediate states obtained at residue-level resolution provide clear evidence that they are very similar to the so-called A-state of Ub that is formed in solutions with low pH and high alcohol. Method validation was carried out by comparing the backbone-amide protection map of native Ub with those deduced from high-resolution NMR measurements.

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.

Biological insights from hydrogen exchange mass spectrometry

Over the past two decades, hydrogen exchange mass spectrometry (HXMS) has achieved the status of a widespread and routine approach in the structural biology toolbox. The ability of hydrogen exchange to detect a range of protein dynamics coupled with the accessibility of mass spectrometry to mixtures and large complexes at low concentrations result in an unmatched tool for investigating proteins challenging to many other structural techniques. Recent advances in methodology and data analysis are helping HXMS deliver on its potential to uncover the connection between conformation, dynamics and the biological function of proteins and complexes. This review provides a brief overview of the HXMS method and focuses on four recent reports to highlight applications that monitor structure and dynamics of proteins and complexes, track protein folding, and map the thermodynamics and kinetics of protein unfolding at equilibrium. These case studies illustrate typical data, analysis and results for each application and demonstrate a range of biological systems for which the interpretation of HXMS in terms of structure and conformational parameters provides unique insights into function. This article is part of a Special Issue entitled: Mass spectrometry in structural biology.

Cutting off functional loops from homodimeric enzyme superoxide dismutase 1 (SOD1) leaves monomeric β-barrels

Demetallation of the homodimeric enzyme Cu/Zn-superoxide dismutase (SOD1) is known to unleash pronounced dynamic motions in the long active-site loops that comprise almost a third of the folded structure. The resulting apo species, which shows increased propensity to aggregate, stands out as the prime disease precursor in amyotrophic lateral sclerosis (ALS). Even so, the detailed structural properties of the apoSOD1 framework have remained elusive and controversial. In this study, we examine the structural interplay between the central apoSOD1 barrel and the active-site loops by simply cutting them off; loops IV and VII were substituted with short Gly-Ala-Gly linkers. The results show that loop removal breaks the dimer interface and leads to soluble, monomeric β-barrels with high structural integrity. NMR-detected nuclear Overhauser effects are found between all of the constituent β-strands, confirming ordered interactions across the whole barrel. Moreover, the breathing motions of the SOD1 barrel are overall insensitive to loop removal and yield hydrogen/deuterium protection factors typical for cooperatively folded proteins (i.e. the active-site loops act as a "bolt-on" domain with little dynamic influence on its structural foundation). The sole exceptions are the relatively low protection factors in β-strand 5 and the turn around Gly-93, a hot spot for ALS-provoking mutations, which decrease even further upon loop removal. Taken together, these data suggest that the cytotoxic function of apoSOD1 does not emerge from its folded ground state but from a high energy intermediate or even from the denatured ensemble.

Native state dynamics drive the unfolding of the SH3 domain of PI3 kinase at high denaturant concentration

Little is known about the role of protein dynamics in directing protein unfolding along a specific pathway and about the role played by chemical denaturants in modulating the dynamics and the initiation of unfolding. In this study, deuterium-hydrogen exchange (HX) detected by electrospray ionization mass spectrometry (ESI-MS) was used to study the unfolding of the SH3 domain of the PI3 kinase. Unfolding on the principal unfolding pathway occurs in 2 steps, both in the absence and in the presence of 1.8 M guanidine hydrochloride (GdnHCl). In both cases, the first step leads to the formation of an intermediate, I(N), with 5 fewer protected amide hydrogen sites than in N. In the second step, I(N) loses the structure protecting the remaining 14 amide hydrogen sites from HX as it unfolds completely. ESI-MS analysis of fragments of the protein created by proteolytic digestion, after completion of the HX reaction, shows that I(N) has lost protection against HX in the same segments of native structure during unfolding in the absence and presence of 1.8 M GdnHCl. Hence, GdnHCl does not appear to play a direct active role in the initiation of unfolding. However, at higher GdnHCl concentrations, a second unfolding pathway is shown to compete effectively with the N <--> I(N) <--> U pathway. In this way, the denaturant modulates the energy landscape of unfolding.

The ALS-associated mutation G93A in human copper-zinc superoxide dismutase selectively destabilizes the remote metal binding region

More than 100 distinct mutations in the gene (SOD1) for human copper-zinc superoxide dismutase (CuZnSOD) have been associated with familial amyotrophic lateral sclerosis (fALS). Studies of these mutant proteins, which often have been performed under far from physiological conditions, have indicated effects on protein stabilities, catalytic activity, and metal binding affinities but with no common pattern. Also, with the knowledge that ALS is a late onset disease it is apparent that protein interactions which contribute to the disorder might, in the natural cellular milieu, depend on a delicate balance between intrinsic protein properties. In this study, we have used experimental conditions as near as possible to the in vivo conditions to reduce artifacts emanating from the experimental setup. Using 1H-15N HSQC NMR spectroscopy, we have analyzed hydrogen exchange at the amide groups of wild-type (wt) CuZnSOD and the fALS-associated G93A SOD variant in their fully metalated states. From analyses of the exchange pattern, we have characterized the local dynamics at 64% of all positions in detail in both the wt and G93A protein. The results show that the G93A mutation had no effect on the dynamics at a majority of the investigated positions. However, the mutation results in local destabilization at the site of the mutation and also in stabilization at a few positions that were apparently scattered over the entire protein surface. Most remarkably, the mutation selectively destabilized the remote metal binding region. The results indicate that the metal binding region may affect the intermolecular protein-protein interactions which cause formation of protein aggregates.

Protein structure and dynamics studied by mass spectrometry: H/D exchange, hydroxyl radical labeling, and related approaches

Mass spectrometry (MS) plays a central role in studies on protein structure and dynamics. This review highlights some of the recent developments in this area, with focus on applications involving the use of electrospray ionization (ESI) MS. Although this technique involves the transformation of analytes into highly nonphysiological species (desolvated gas-phase ions in the vacuum), ESI-MS can provide detailed insights into the solution-phase behavior of proteins. Notably, the ionization process itself occurs in a structurally sensitive manner. An increased degree of solution-phase unfolding is correlated with a higher level of protonation. Also, ESI allows the transfer of intact noncovalent complexes into the gas phase, thereby yielding information on binding partners, stoichiometries, and even affinities. A particular focus of this article is the use of hydrogen/deuterium exchange (HDX) methods and hydroxyl radical (.OH) labeling for monitoring dynamic and structural aspect of solution-phase proteins. Conceptual similarities and differences between the two methods are discussed. We describe a simple method for the computational simulation of protein HDX patterns, a tool that can be helpful for the interpretation of isotope exchange data recorded under mixed EX1/EX2 conditions. Important aspects of .OH labeling include a striking dependence on protein concentration, and the tendency of commonly used solvent additives to act as highly effective radical scavengers. If not properly controlled, both of these factors may lead to experimental artifacts.