Hydrogen Exchange Data Research Articles

Effects of ligand binding on dynamics of fatty acid binding protein and interactions with membranes

Intracellular transport of fatty acids involves binding of ligands to their carrier fatty acid binding proteins (FABPs) and interactions of ligand-free and -bound FABPs with membranes. Previous studies focused on ligand-free FABPs. Here, our amide hydrogen exchange data showed that oleic acid binding to human intestinal FABP (hIFABP) stabilizes the protein, most likely through enhancing the hydrogen-bonding network, and induces rearrangement of sidechains even far away from the ligand binding site. Using NMR relaxation techniques, we found that the ligand binding affects not only conformational exchanges between major and minor states but also the affinity of hIFABP to nanodiscs. Analyses of the relaxation and amide exchange data suggested that two minor native-like states existing in both ligand-free and -bound hIFABPs originate from global “breathing” motions, while one minor native-like state comes from local motions. The amide hydrogen exchange data also indicated that helix αII undergoes local unfolding through which ligands can exit from the binding cavity.

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.

Toward a proper interpretation of hydrogen exchange data in disordered proteins

Amide protons readily exchange with solvent water. Above pH ~4 the exchange occurs as a direct result of ionic attack by a hydroxide ion, OH–. Measurements of solvent exchange, often involving deuterium-for-hydrogen substitution, is one of the essential NMR and mass spectrometry techniques. Of note, these experiments allow one to “see the invisible,” i.e., obtain information about the effects that are otherwise difficult or impossible to observe. For example, these experiments alerted protein science community to the fact that fully structured globular proteins occasionally experience major meltdowns (1).

Modeling hydrogen exchange of proteins by a multiscale method* *Project supported by the National Natural Science Foundation of China (Grant Nos. 11974173 and 11934008) and the HPC Center of Nanjing University.

We proposed a practical way for mapping the results of coarse-grained molecular simulations to the observables in hydrogen change experiments. By combining an atomic-interaction based coarse-grained model with an all-atom structure reconstruction algorithm, we reproduced the experimental hydrogen exchange data with reasonable accuracy using molecular dynamics simulations. We also showed that the coarse-grained model can be further improved by imposing experimental restraints from hydrogen exchange data via an iterative optimization strategy. These results suggest that it is feasible to develop an integrative molecular simulation scheme by incorporating the hydrogen exchange data into the coarse-grained molecular dynamics simulations and therefore help to overcome the accuracy bottleneck of coarse-grained models.

SWOTein: a structure-based approach to predict stability Strengths and Weaknesses of prOTEINs.

Although structured proteins adopt their lowest free energy conformation in physiological conditions, the individual residues are generally not in their lowest free energy conformation. Residues that are stability weaknesses are often involved in functional regions, whereas stability strengths ensure local structural stability. The detection of strengths and weaknesses provides key information to guide protein engineering experiments aiming to modulate folding and various functional processes. We developed the SWOTein predictor which identifies strong and weak residues in proteins on the basis of three types of statistical energy functions describing local interactions along the chain, hydrophobic forces and tertiary interactions. The large-scale analysis of the different types of strengths and weaknesses demonstrated their complementarity and the enhancement of the information they provide. Moreover, a good average correlation was observed between predicted and experimental strengths and weaknesses obtained from native hydrogen exchange data. SWOTein application to three test cases further showed its suitability to predict and interpret strong and weak residues in the context of folding, conformational changes and protein-protein binding. In summary, SWOTein is both fast and accurate and can be applied at small and large scale to analyze and modulate folding and molecular recognition processes. The SWOTein webserver provides the list of predicted strengths and weaknesses and a protein structure visualization tool that facilitates the interpretation of the predictions. It is freely available for academic use at http://babylone.ulb.ac.be/SWOTein/.

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.

POPPeT: a New Method to Predict the Protection Factor of Backbone Amide Hydrogens

Hydrogen exchange (HX) has become an important tool to monitor protein structure and dynamics. The interpretation of HX data with respect to protein structure requires understanding of the factors that influence exchange. Simulated protein structures can be validated by comparing experimental deuteration profiles with the profiles derived from the modeled protein structure. To do this, we propose here a new method, POPPeT, for protection factor prediction based on protein motions that enable HX. By comparing POPPeT with two existing methods, the phenomenological approximation and COREX, we show enhanced predictability measured at both protection factor and deuteration level. This method can be subsequently used by modeling strategies for protein structure prediction.Graphical ᅟ

Empirical Correction for Differences in Chemical Exchange Rates in Hydrogen Exchange-Mass Spectrometry Measurements.

A barrier to the use of hydrogen exchange-mass spectrometry (HX-MS) in many contexts, especially analytical characterization of various protein therapeutic candidates, is that differences in temperature, pH, ionic strength, buffering agent, or other additives can alter chemical exchange rates, making HX data gathered under differing solution conditions difficult to compare. Here, we present data demonstrating that HX chemical exchange rates can be substantially altered not only by the well-established variables of temperature and pH but also by additives including arginine, guanidine, methionine, and thiocyanate. To compensate for these additive effects, we have developed an empirical method to correct the hydrogen-exchange data for these differences. First, differences in chemical exchange rates are measured by use of an unstructured reporter peptide, YPI. An empirical chemical exchange correction factor, determined by use of the HX data from the reporter peptide, is then applied to the HX measurements obtained from a protein of interest under different solution conditions. We demonstrate that the correction is experimentally sound through simulation and in a proof-of-concept experiment using unstructured peptides under slow-exchange conditions (pD 4.5 at ambient temperature). To illustrate its utility, we applied the correction to HX-MS excipient screening data collected for a pharmaceutically relevant IgG4 mAb being characterized to determine the effects of different formulations on backbone dynamics.

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.

NMR Investigation of Structures of G-protein Coupled Receptor Folding Intermediates

Folding of G-protein coupled receptors (GPCRs) according to the two-stage model (Popot, J. L., and Engelman, D. M. (1990) Biochemistry 29, 4031-4037) is postulated to proceed in 2 steps: partitioning of the polypeptide into the membrane followed by diffusion until native contacts are formed. Herein we investigate conformational preferences of fragments of the yeast Ste2p receptor using NMR. Constructs comprising the first, the first two, and the first three transmembrane (TM) segments, as well as a construct comprising TM1-TM2 covalently linked to TM7 were examined. We observed that the isolated TM1 does not form a stable helix nor does it integrate well into the micelle. TM1 is significantly stabilized upon interaction with TM2, forming a helical hairpin reported previously (Neumoin, A., Cohen, L. S., Arshava, B., Tantry, S., Becker, J. M., Zerbe, O., and Naider, F. (2009) Biophys. J. 96, 3187-3196), and in this case the protein integrates into the hydrophobic interior of the micelle. TM123 displays a strong tendency to oligomerize, but hydrogen exchange data reveal that the center of TM3 is solvent exposed. In all GPCRs so-far structurally characterized TM7 forms many contacts with TM1 and TM2. In our study TM127 integrates well into the hydrophobic environment, but TM7 does not stably pack against the remaining helices. Topology mapping in microsomal membranes also indicates that TM1 does not integrate in a membrane-spanning fashion, but that TM12, TM123, and TM127 adopt predominantly native-like topologies. The data from our study would be consistent with the retention of individual helices of incompletely synthesized GPCRs in the vicinity of the translocon until the complete receptor is released into the membrane interior.

Mechanism of Protein Denaturation: Partial Unfolding of the P22 Coat Protein I-Domain by Urea Binding

The I-domain is an insertion domain of the bacteriophage P22 coat protein that drives rapid folding and accounts for over half of the stability of the full-length protein. We sought to determine the role of hydrogen bonds (H-bonds) in the unfolding of the I-domain by examining 3JNC’ couplings transmitted through H-bonds, the temperature and urea-concentration dependence of 1HN and 15N chemical shifts, and native-state hydrogen exchange at urea concentrations where the domain is predominantly folded. The native-state hydrogen-exchange data suggest that the six-stranded β-barrel core of the I-domain is more stable against unfolding than a smaller subdomain comprised of a short α-helix and three-stranded β-sheet. H-bonds, separately determined from solvent protection and 3JNC’ H-bond couplings, are identified with an accuracy of 90% by 1HN temperature coefficients. The accuracy is improved to 95% when 15N temperature coefficients are also included. In contrast, the urea dependence of 1HN and 15N chemical shifts is unrelated to H-bonding. The protein segments with the largest chemical-shift changes in the presence of urea show curved or sigmoidal titration curves suggestive of direct urea binding. Nuclear Overhauser effects to urea for these segments are also consistent with specific urea-binding sites in the I-domain. Taken together, the results support a mechanism of urea unfolding in which denaturant binds to distinct sites in the I-domain. Disordered segments bind urea more readily than regions in stable secondary structure. The locations of the putative urea-binding sites correlate with the lower stability of the structure against solvent exchange, suggesting that partial unfolding of the structure is related to urea accessibility.

Estimation of Hydrogen-Exchange Protection Factors from MD Simulation Based on Amide Hydrogen Bonding Analysis.

Hydrogen exchange (HX) studies have provided critical insight into our understanding of protein folding, structure, and dynamics. More recently, hydrogen exchange mass spectrometry (HX-MS) has become a widely applicable tool for HX studies. The interpretation of the wealth of data generated by HX-MS experiments as well as other HX methods would greatly benefit from the availability of exchange predictions derived from structures or models for comparison with experiment. Most reported computational HX modeling studies have employed solvent-accessible-surface-area based metrics in attempts to interpret HX data on the basis of structures or models. In this study, a computational HX-MS prediction method based on classification of the amide hydrogen bonding modes mimicking the local unfolding model is demonstrated. Analysis of the NH bonding configurations from molecular dynamics (MD) simulation snapshots is used to determine partitioning over bonded and nonbonded NH states and is directly mapped into a protection factor (PF) using a logistics growth function. Predicted PFs are then used for calculating deuteration values of peptides and compared with experimental data. Hydrogen exchange MS data for fatty acid synthase thioesterase (FAS-TE) collected for a range of pHs and temperatures was used for detailed evaluation of the approach. High correlation between prediction and experiment for observable fragment peptides is observed in the FAS-TE and additional benchmarking systems that included various apo/holo proteins for which literature data were available. In addition, it is shown that HX modeling can improve experimental resolution through decomposition of in-exchange curves into rate classes, which correlate with prediction from MD. Successful rate class decompositions provide further evidence that the presented approach captures the underlying physical processes correctly at the single residue level. This assessment is further strengthened in a comparison of residue resolved protection factor predictions for staphylococcal nuclease with NMR data, which was also used to compare prediction performance with other algorithms described in the literature. The demonstrated transferable and scalable MD based HX prediction approach adds significantly to the available tools for HX-MS data interpretation based on available structures and models.

HX‐MS2 for high performance conformational analysis of complex protein states

Water-mediated hydrogen exchange (HX) processes involving the protein main chain are sensitive to structural dynamics and molecular interactions. Measuring deuterium uptake in amide bonds provides information on conformational states, structural transitions and binding events. Increasingly, deuterium levels are measured by mass spectrometry (MS) from proteolytically generated peptide fragments of large molecular systems. However, this bottom-up method has limited spectral capacity and requires a burdensome manual validation exercise, both of which restrict analysis of protein systems to generally less than 150 kDa. In this study, we present a bottom-up HX-MS(2) method that improves peptide identification rates, localizes high-quality HX data and simplifies validation. The method combines a new peptide scoring algorithm (WUF, weighted unique fragment) with data-independent acquisition of peptide fragmentation data. Scoring incorporates the validation process and emphasizes identification accuracy. The HX-MS(2) method is illustrated using data from a conformational analysis of microtubules treated with dimeric kinesin MCAK. When compared to a conventional Mascot-driven HX-MS method, HX-MS(2) produces two-fold higher α/β-tubulin sequence depth at a peptide utilization rate of 74%. A Mascot approach delivers a utilization rate of 44%. The WUF score can be constrained by false utilization rate (FUR) calculations to return utilization values exceeding 90% without serious data loss, indicating that automated validation should be possible. The HX-MS(2) data confirm that N-terminal MCAK domains anchor kinesin force generation in kinesin-mediated depolymerization, while the C-terminal tails regulate MCAK-tubulin interactions.

Protein destabilisation in ionic liquids: the role of preferential interactions in denaturation

The preferential binding of anions and cations in aqueous solutions of the ionic liquids (ILs) 1-butyl-3-methylimidazolium ([C4mim](+)) and 1-ethyl-3-methylimidazolium ([C2mim](+)) chloride and dicyanamide (dca(-)) with the small alpha-helical protein Im7 was investigated using a combination of differential scanning calorimetry, NMR spectroscopy and molecular dynamics (MD) simulations. Our results show that direct ion interactions are crucial to understand the effects of ILs on the stability of proteins and that an anion effect is dominant. We show that the binding of weakly hydrated anions to positively charged or polar residues leads to the partial dehydration of the backbone groups, and is critical to control stability, explaining why dca(-) is more denaturing than Cl(-). Direct cation-protein interactions also mediate stability; cation size and hydrophobicity are relevant to account for destabilisation as shown by the effect of [C4mim](+) compared to [C2mim](+). The specificity in the interaction of IL ions with protein residues established by weak favourable interactions is confirmed by NMR chemical shift perturbation, amide hydrogen exchange data and MD simulations. Differences in specificity are due to the balance of interaction established between ion pairs and ion-solvent that determine the type of residues affected. When the interaction of both cation and anion with the protein is strong the net result is similar to a non-specific interaction, leading ultimately to unfolding. Since the nature of the ions is a determinant of the level of interaction with the protein towards denaturation or stabilisation, ILs offer a unique possibility to modulate protein stabilisation or even folding events.

Structure and Dynamics of the ATP-Bound Open Conformation of Hsp70 Chaperones

Central to the chaperone function of Hsp70s is the transition between open and closed conformations of their polypeptide substrate binding domain (SBD), which is regulated through an allosteric mechanism via ATP binding and hydrolysis in their nucleotide binding domain (NBD). Although the structure of the closed conformation of Hsp70s is well studied, the open conformation has remained elusive. Here, we report on the 2.4Å crystal structure of the ATP-bound open conformation of the Escherichia coli Hsp70 homolog DnaK. In the open DnaK structure, the β sheet and α-helical lid subdomains of the SBD are detached from one another and docked to different faces of the NBD. The contacts between the β sheet subdomain and the NBD reveal the mechanism of allosteric regulation. In addition, we demonstrate that docking of the β sheet and α-helical lid subdomains to the NBD is a sequential process influenced by peptide and protein substrates.

The N-terminal Capping Propensities of the D-helix Modulate the Allosteric Activation of the Escherichia coli cAMP Receptor Protein

Transduction of biological signals at the molecular level involves the activation and/or inhibition of allosteric proteins. In the transcription factor cAMP receptor protein (CRP) from Escherichia coli, the allosteric activation, or apo-holo transition, involves rigid body motions of domains and structural rearrangements within the hinge region connecting the cAMP- and DNA-binding domains. During this apo-holo transition, residue 138 is converted as part of the elongated D-helix to the position of the N-terminal capping residue of a shorter D-helix. The goal of the current study is to elucidate the role of residue 138 in modulating the allostery between cAMP and DNA binding. By systematically mutating residue 138, we found that mutants with higher N-terminal capping propensities lead to increased cooperativity of cAMP binding and a concomitant increase in affinity for lac-DNA. Furthermore, mutants with higher N-terminal capping propensity correlate with properties characteristic of holo-CRP, particularly, increase in protein structural dynamics. Overall, our results provide a quantitative characterization of the role of residue 138 in the isomerization equilibrium between the apo and holo forms of CRP, and in turn the thermodynamic underpin to the molecular model of allostery revealed by the high resolution structural studies.

Protein dynamics viewed by hydrogen exchange

To examine the relationship between protein structural dynamics and measurable hydrogen exchange (HX) data, the detailed exchange behavior of most of the backbone amide hydrogens of Staphylococcal nuclease was compared with that of their neighbors, with their structural environment, and with other information. Results show that H-bonded hydrogens are protected from exchange, with HX rate effectively zero, even when they are directly adjacent to solvent. The transition to exchange competence requires a dynamic structural excursion that removes H-bond protection and allows exposure to solvent HX catalyst. The detailed data often make clear the nature of the dynamic excursion required. These range from whole molecule unfolding, through smaller cooperative unfolding reactions of secondary structural elements, and down to local fluctuations that involve as little as a single peptide group or side chain or water molecule. The particular motion that dominates the exchange of any hydrogen is the one that allows the fastest HX rate. The motion and the rate it produces are determined by surrounding structure and not by nearness to solvent or the strength of the protecting H-bond itself or its acceptor type (main chain, side chain, structurally bound water). Many of these motions occur over time scales that are appropriate for biochemical function.

Experimentally assessing molecular dynamics sampling of the protein native state conformational distribution

The acute sensitivity to conformation exhibited by amide hydrogen exchange reactivity provides a valuable test for the physical accuracy of model ensembles developed to represent the Boltzmann distribution of the protein native state. A number of molecular dynamics studies of ubiquitin have predicted a well-populated transition in the tight turn immediately preceding the primary site of proteasome-directed polyubiquitylation Lys 48. Amide exchange reactivity analysis demonstrates that this transition is 103-fold rarer than these predictions. More strikingly, for the most populated novel conformational basin predicted from a recent 1ms MD simulation of bovine pancreatic trypsin inhibitor (at 13% of total), experimental hydrogen exchange data indicates a population below 10−6. The most sophisticated efforts to directly incorporate experimental constraints into the derivation of model protein ensembles have been applied to ubiquitin, as illustrated by three recently deposited studies (PDB codes 2NR2, 2K39 and 2KOX2K392KOX). Utilizing the extensive set of experimental NOE constraints, each of these three ensembles yields a modestly more accurate prediction of the exchange rates for the highly exposed amides than does a standard unconstrained molecular simulation. However, for the less frequently exposed amide hydrogens, the 2NR2 ensemble offers no improvement in rate predictions as compared to the unconstrained MD ensemble. The other two NMR-constrained ensembles performed markedly worse, either underestimating (2KOX) or overestimating (2K39) the extent of conformational diversity.

Quantitative Assessment of Protein Structural Models by Comparison of H/D Exchange MS Data with Exchange Behavior Accurately Predicted by DXCOREX

Peptide amide hydrogen/deuterium exchange mass spectrometry (DXMS) data are often used to qualitatively support models for protein structure. We have developed and validated a method (DXCOREX) by which exchange data can be used to quantitatively assess the accuracy of three-dimensional (3-D) models of protein structure. The method utilizes the COREX algorithm to predict a protein's amide hydrogen exchange rates by reference to a hypothesized structure, and these values are used to generate a virtual data set (deuteron incorporation per peptide) that can be quantitatively compared with the deuteration level of the peptide probes measured by hydrogen exchange experimentation. The accuracy of DXCOREX was established in studies performed with 13 proteins for which both high-resolution structures and experimental data were available. The DXCOREX-calculated and experimental data for each protein was highly correlated. We then employed correlation analysis of DXCOREX-calculated versus DXMS experimental data to assess the accuracy of a recently proposed structural model for the catalytic domain of a Ca(2+)-independent phospholipase A(2). The model's calculated exchange behavior was highly correlated with the experimental exchange results available for the protein, supporting the accuracy of the proposed model. This method of analysis will substantially increase the precision with which experimental hydrogen exchange data can help decipher challenging questions regarding protein structure and dynamics.

Conformer-Specific Hydrogen Exchange Analysis of Aβ(1–42) Oligomers by Top-Down Electron Capture Dissociation Mass Spectrometry

Protein structural studies are particularly challenging under conditions in which several conformational species (e.g., monomers and aggregated forms) coexist in solution. Most spectroscopic techniques provide population-averaged data. Hence, it is usually not possible to obtain detailed structural information on individual protein species in heterogeneous samples. The current work employs an experimental strategy that addresses this issue. Solution-phase hydrogen exchange (HX) is used in combination with tandem mass spectrometry. Electrosprayed intact ions exhibiting specific HX mass shifts are selected in the gas phase, followed by electron capture dissociation. The resulting fragment ion deuteration pattern provides amide hydrogen bonding information in a conformer-specific and spatially resolved fashion. The feasibility of this approach is demonstrated by applying it to neurotoxic Aβ(1-42) oligomers that coexist with disordered monomers in solution. The findings of this study point to similarities between oligomers and mature amyloid fibrils with regard to the Aβ(1-42) backbone organization. Specifically, fibrils and oligomers appear to share a β-loop-β secondary structure motif. The spatial resolution obtained with the "top-down" approach used here exceeds that of earlier proteolysis-based HX data on Aβ.