Automated extraction of backbone deuteration levels from amide H/2H mass spectrometry experiments

A Fourier deconvolution method has been developed to explicitly determine the amount of backbone amide deuterium incorporated into protein regions or segments by hydrogen/deuterium (H/D) exchange with high-resolution mass spectrometry. Determination and analysis of the level and number of backbone amide exchanging in solution provide more information about the solvent accessibility of the protein than do previous centroid methods, which only calculate the average deuterons exchanged. After exchange, a protein is digested into peptides as a way of determining the exchange within a local area of the protein. The mass of a peptide upon deuteration is a sum of the natural isotope abundance, fast exchanging side-chain hydrogens (present in MALDI-TOF H/2H data) and backbone amide exchange. Removal of the components of the isotopic distribution due to the natural isotope abundances and the fast exchanging side-chains allows for a precise quantification of the levels of backbone amide exchange, as is shown by an example from protein kinase A. The deconvoluted results are affected by overlapping peptides or inconsistent mass envelopes, and evaluation procedures for these cases are discussed. Finally, a method for determining the back exchange corrected populations is presented, and its effect on the data is discussed under various circumstances.     

Determinants of protein hydrogen exchange studied in equine cytochrome c.

The exchange of a large number of amide hydrogens in oxidized equine cytochrome c was measured by NMR and compared with structural parameters. Hydrogens known to exchange through local structural fluctuations and through larger unfolding reactions were separately considered. All hydrogens protected from exchange by factors greater than 10(3) are in defined H-bonds, and almost all H-bonded hydrogens including those at the protein surface were measured to exchange slowly. H-exchange rates do not correlate with H-bond strength (length) or crystallographic B factors. It appears that the transient structural fluctuation necessary to bring an exchangeable hydrogen into H-bonding contact with the H-exchange catalyst (OH(-)-ion) involves a fairly large separation of the H-bond donor and acceptor, several angstroms at least, and therefore depends on the relative resistance to distortion of immediately neighboring structure. Accordingly, H-exchange by way of local fluctuational pathways tends to be very slow for hydrogens that are neighbored by tightly anchored structure and for hydrogens that are well buried. The slowing of buried hydrogens may also reflect the need for additional motions that allow solvent access once the protecting H-bond is separated, although it is noteworthy that burial in a protein like cytochrome c does not exceed 4 angstroms. When local fluctuational pathways are very slow, exchange can become dominated by a different category of larger, cooperative, segmental unfolding reactions reaching up to global unfolding.     

Hydrogen exchange shows peptide binding stabilizes motions in Hck SH2.

Src-homology-2 domains are small, 100 amino acid protein modules that are present in a number of signal transduction proteins. Previous NMR studies of SH2 domain dynamics indicate that peptide binding decreases protein motions in the pico- to nanosecond, and perhaps slower, time range. We suggest that amide hydrogen exchange and mass spectrometry may be useful for detecting changes in protein dynamics because hydrogen exchange rates are relatively insensitive to the time domains of the dynamics. In the present study, hydrogen exchange and mass spectrometry were used to probe hematopoietic cell kinase SH2 that was either free or bound to a 12-residue high-affinity peptide. Hydrogen exchange rates were determined by exposing free and bound SH2 to D(2)O, fragmenting the SH2 with pepsin, and determining the deuterium level in the peptic fragments. Binding generally decreased hydrogen exchange along much of the SH2 backbone, indicating a widespread reduction in dynamics. Alterations in the exchange of the most rapidly exchanging amide hydrogens, which was detected following acid quench and analysis by mass spectrometry, were used to locate differences in low-amplitude motion when SH2 was bound to the peptide. In addition, the results indicate that hydrogen exchange from the folded form of SH2 is an important process along the entire SH2 backbone.     

Re-examination of rhodopsin structure by hydrogen exchange

The hydrogen exchange behavior of rhodopsin was re-examined by studies of the protein in the disc membrane and after solubilization in octyl glucoside. The methods used measure either the peptide hydrogens alone (hydrogen-deuterium exchange by infrared spectroscopy) or all slowly exchanging hydrogens (hydrogen-tritium exchange by hel filtration). Under mild exchange conditions, disc membranes and solubilized lipid-free proteins show very similar exchange behavior, indicating the absence of slowly exchanging lipid protons. At high temperature, exchange of an additional large group of very slow peptide NH can be detected. The total number of slow hydrogens significantly exceeds the amide content, and apparently includes slowly exchanging protons from perhaps 40% of the protein's non-amide side chains. This is thought to require the involvement of many polar side chains in internal H-bonding. The exchange rates of the non-amide side chains sites have not been determined. However, to the extent that these contribute to the fast time region of the measured kinetic H-exchange curve, previously identified with exposed, non-H-bonded peptides, the estimate of freely exposed rhodopsin peptides must be reduced. The fraction of free peptides could range from a remarkably high value of 70% down to about 45%.     

Stopped-flow NMR measurement of hydrogen exchange rates in reduced horse cytochrome c under strongly destabilizing conditions

A procedure to measure exchange rates of fast exchanging protein amide hydrogens by time-resolved NMR spectroscopy following in situ initiation of the reaction by diluting a native protein solution into an exchanging deuterated buffer is described. The method has been used to measure exchange rates of a small set of amide hydrogens of reduced cytochrome c, maintained in a strictly anaerobic atmosphere, in the presence of an otherwise inaccessible range of guanidinium deuterochloride concentrations. The results for the measured protons indicate that hydrogen exchange in the unfolding transition region of cytochrome c reach the EX2 limit, but emphasize the difficulty in interpretation of the exchange mechanism in protein hydrogen exchange studies. Comparison of free energies of structure opening for the measured hydrogens with the global unfolding free energy monitored by far-UV CD measurements has indicated the presence of at least one partially unfolded equilibrium species of reduced cytochrome c. The results provide the first report of measurement of free energy of opening of structure to exchange in the 0–2-kcal/mol range. Proteins 32:241–247, 1998. © 1998 Wiley-Liss, Inc.     

Hydrogen/deuterium exchange and mass spectrometric analysis of a protein containing multiple disulfide bonds: Solution structure of recombinant macrophage colony stimulating factor-beta (rhM-CSFbeta)

Studies with the homodimeric recombinant human macrophage colony-stimulating factor beta (rhM-CSFbeta), show for the first time that a large number (9) of disulfide linkages can be reduced after amide hydrogen/deuterium (H/D) exchange, and the protein digested and analyzed successfully for the isotopic composition by electrospray mass spectrometry. Analysis of amide H/D after exchange-in shows that in solution the conserved four-helix bundle of (rhM-CSFbeta) has fast and moderately fast exchangeable sections of amide hydrogens in the alphaA helix, and mostly slow exchanging sections of amide hydrogens in the alphaB, alphaC, and alphaD helices. Most of the amide hydrogens in the loop between the beta1 and beta4 sheets exhibited fast or moderately fast exchange, whereas in the amino acid 63-67 loop, located at the interface of the two subunits, the exchange was slow. Solvent accessibility as measured by H/D exchange showed a better correlation with the average depth of amide residues calculated from reported X-ray crystallographic data for rhM-CSFalpha than with the average B-factor. The rates of H/D exchange in rhM-CSFbeta appear to correlate well with the exposed surface calculated for each amino acid residue in the crystal structure except for the alphaD helix. Fast hydrogen isotope exchange throughout the segment amino acids 150-221 present in rhM-CSFbeta, but not rhM-CSFalpha, provides evidence that the carboxy-terminal region is unstructured. It is, therefore, proposed that the anomalous behavior of the alphaD helix is due to interaction of the carboxy-terminal tail with this helical segment.     

Hydrogen exchange-mass spectrometry: optimization of digestion conditions

The direct linkage between folded structures of proteins and their function has increased the need for high resolution structures. In addition, there is a need for analytical methods for detecting and locating changes in the folded structures of proteins under a wide variety of conditions. The rates at which hydrogens located at peptide amide linkages undergo isotopic exchange has become the basis for an important method for detecting such structural changes. When detected by mass spectrometry, hydrogen exchange can be used to study dilute solutions of large proteins and protein complexes with very high sensitivity. To locate structural changes, labeled proteins are often digested with acid proteases to form peptides whose hydrogen/deuterium levels are determined by mass spectrometry. This approach is successful only when the protein can be digested rapidly under conditions where isotope exchange is slow. This study describes how columns packed with immobilized pepsin can be used to reduce the digestion time and to provide an effective means for separating the pepsin from the isotopically labeled fragments. These columns are part of an on-line system that facilitates both rapid digestion of low concentrations of protein and concentration of the peptides.     

Protein hydrogen exchange: testing current models

To investigate the determinants of protein hydrogen exchange (HX), HX rates of most of the backbone amide hydrogens of Staphylococcal nuclease were measured by NMR methods. A modified analysis was used to improve accuracy for the faster hydrogens. HX rates of both near surface and well buried hydrogens are spread over more than 7 orders of magnitude. These results were compared with previous hypotheses for HX rate determination. Contrary to a common assumption, proximity to the surface of the native protein does not usually produce fast exchange. The slow HX rates for unprotected surface hydrogens are not well explained by local electrostatic field. The ability of buried hydrogens to exchange is not explained by a solvent penetration mechanism. The exchange rates of structurally protected hydrogens are not well predicted by algorithms that depend only on local interactions or only on transient unfolding reactions. These observations identify some of the present difficulties of HX rate prediction and suggest the need for returning to a detailed hydrogen by hydrogen analysis to examine the bases of structure-rate relationships, as described in the companion paper (Skinner et al., Protein Sci 2012;21:996-1005).     

Human recombinant [C22A] FK506-binding protein amide hydrogen exchange rates from mass spectrometry match and extend those from NMR.

Hydrogen/deuterium exchange behavior of human recombinant [C22A] FK506 binding protein (C22A FKBP) has been determined by protein fragmentation, combined with electrospray Fourier transform ion cyclotron resonance mass spectrometry (MS). After a specified period of H/D exchange in solution, C22A FKBP was digested by pepsin under slow exchange conditions (pH 2.4, 0 degree C), and then subjected to on-line HPLC/MS for deuterium analysis of each proteolytic peptide. The hydrogen exchange rate of each individual amide hydrogen was then determined independently by heteronuclear two-dimensional NMR on 15N-enriched C22A FKBP. A maximum entropy method (MEM) algorithm makes it possible to derive the distributions of hydrogen exchange rate constants from the MS-determined deuterium exchange-in curves in either the holoprotein or its proteolytic segments. The MEM-derived rate constant distributions of C22A FKBP and different segments of C22A FKBP are compared to the rate constants determined by NMR for individual amide protons. The rate constant distributions determined by both methods are consistent and complementary, thereby validating protein fragmentation/mass spectrometry as a reliable measure of hydrogen exchange in proteins.     

Structural basis of human erythrocyte glucose transporter function in reconstituted system. Hydrogen exchange

Hydrogen exchange kinetic behavior of human erythrocyte glucose transporter protein in vesicles was studied in the absence and in the presence of D-glucose or a well known inhibitor, cytochalasin B. This is to detect a proposed channel of water penetrating into the protein through which the sugar molecule passes and to monitor any conformational changes induced by the substrate or inhibitor. Analyses of the kinetic data revealed several classes of hydrogens which exchange with readily distinguishable rates. Of 660 hydrogens detected per transporter, approximately 30% exchanged with rates generally characterized as those of free amide hydrogens indicating they are interfaced to solvent water. Since the transporter is known to be embedded deep in the hydrophobic area of the membrane with minimum exposure to the outside of the membrane lipid bilayer, a significant portion of these free amide hydrogens must be at the purported channel rather than outside of the membrane. D-Glucose and cytochalasin B affected the exchange kinetics of these presumably channel-associated free amide hydrogens rather differently. D-Glucose reduced the apparent rate constants, but not the total number. Cytochalasin B on the other hand reduced the total number to one-half without significantly changing the apparent rate constants. The remaining 70% of the labeled hydrogens exchanged with much slower rates which vary 10-10,000-fold, indicating that they are internally structured peptide amide and side chain hydrogens. Both D-glucose and cytochalasin B further reduced the rates of these hydrogens, indicating a global stabilization of the protein structure.     

Estimating Constraints for Protection Factors from HDX-MS Data

Hydrogen/deuterium exchange monitored by mass spectrometry is a promising technique for rapidly fingerprinting structural and dynamical properties of proteins. The time-dependent change in the mass of any fragment of the polypeptide chain depends uniquely on the rate of exchange of its amide hydrogens, but determining the latter from the former is generally not possible. Here, we show that, if time-resolved measurements are available for a number of overlapping peptides that cover the whole sequence, rate constants for each amide hydrogen exchange (or equivalently, their protection factors) may be extracted and the uniqueness of the solutions obtained depending on the degree of peptide overlap. However, in most cases, the solution is not unique, and multiple alternatives must be considered. We provide a statistical method that clusters the solutions to further reduce their number. Such analysis always provides meaningful constraints on protection factors and can be used in situations in which obtaining more refined experimental data is impractical. It also provides a systematic way to improve data collection strategies to obtain unambiguous information at single-residue level (e.g., for assessing protein structure predictions at atomistic level).     

Hydrogen exchange behavior of [U-15N]-labeled oxidized and reduced iso-1-cytochrome c

Heteronuclear NMR spectroscopy was used to measure the hydrogen-deuterium exchange rates of backbone amide hydrogens in both oxidized and reduced [U-15N]iso-1-cytochrome c from the yeast Saccharomyces cerevisiae. The exchange data confirm previously reported data [Marmorino et al. (1993) Protein Sci. 2, 1966-1974], resolve several inconsistencies, and provide more thorough coverage of exchange rates throughout the cytochrome c protein in both oxidation states. Combining the data previously collected on unlabeled C102T with the current data collected on [U-15N]C102T, exchange rates for 53 protons in the oxidized state and 52 protons in the reduced state can now be reported. Most significantly, hydrogen exchange measurements on [U-15N]iso-1-cytochrome c allowed the observation of exchange behavior of the secondary structures, such as large loops, that are not extensively hydrogen-bonded. For the helices, the most slowly exchanging protons are found in the middle of the helix, with more rapidly exchanging protons at the helix ends. The observation for the Omega-loops in cytochrome c is just the opposite. In the loops, the ends contain the most slowly exchanging protons and the loop middles allow more rapid exchange. This is found to be true in cytochrome c loops, even though the loop ends are not attached to any regular secondary structures. Some of the exchange data are strikingly inconsistent with data collected on the C102S variant at a different pH, which suggests pH-dependent dynamic differences in the protein structure. This new hydrogen exchange data for loop residues could have implications for the substructure model of eukaryotic cytochrome c folding. Isotopic labeling of variant forms of cytochrome c can now be used to answer many questions about the structure and folding of this model protein.     

Capsid structure and dynamics of a human rhinovirus probed by hydrogen exchange mass spectrometry

Viral capsids are dynamic protein assemblies surrounding viral genomes. Despite the high-resolution structures determined by X-ray crystallography and cryo-electron microscopy, their in-solution structure and dynamics can be probed by hydrogen exchange. We report here using hydrogen exchange combined with protein enzymatic fragmentation and mass spectrometry to determine the capsid structure and dynamics of a human rhinovirus, HRV14. Capsid proteins (VP1-4) were labeled with deuterium by incubating intact virus in D(2)O buffer at neutral pH. The labeled proteins were digested by immobilized pepsin to give peptides analyzed by capillary reverse-phase HPLC coupled with nano-electrospray mass spectrometry. Deuterium levels incorporated at amide linkages in peptic fragments were measured for different exchange times from 12 sec to 30 h to assess the amide hydrogen exchange rates along each of the four protein backbones. Exchange results generally agree with the crystal structure of VP1-4,with extended, flexible terminal and surface-loop regions in fast exchange and folded helical and sheet structures in slow exchange. In addition, three alpha-helices, one from each of VP1-3, exhibited very slow exchange, indicating high stability of the protomeric interface. The beta-strands at VP3 N terminus also had very slow exchange, suggesting stable pentamer contacts. It was noted, however, that the interface around the fivefold axis had fast and intermediate exchange, indicating relatively more flexibility. Even faster exchange rates were found in the N terminus of VP1 and most segments of VP4, suggesting high flexibilities, which may correspond to their potential roles in virus uncoating.     

Hydrogen-isotope exchange of oxidized and reduced cytochrome c. A comparison of mass spectrometry and infrared methods.

Hydrogen-deuterium exchange in 2H20 solutions of the two redox states of horse heart cytochrome c was investigated at 20 degrees C, pH 7, by mass spectrometry and infrared spectroscopy. Mass spectrometry indicates that ferricytochrome has 20 hydrogens unexchanged after 24 h, 28 hydrogens exchanging between 10 min and 24 h, and 156 hydrogens exchanging within 10 min; comparative values for ferrocytochrome are 45, 19 and 140. The displacement of the exchange curves obtained by infrared corresponds to 8 to 9 peptide hydrogens. These combined methods show many non-peptide hydrogens exchanging rapidly (87 and 79 for ferricytochrome c and ferrocytochrome c respectively), whereas others, probably buried inside the molecule and involved in hydrogen bonds, are not exchanged, even after 24 h (14 and 30 hydrogens respectively, which is relatively large for a small protein). Infrared results are given in terms of changes of standard free energy for the transconformational reaction which exposes the peptide hydrogens to solvent: in ferricytochrome c and ferrycoytochrome c, 30% and 40% respectively of the peptide hydrogens are protected by conformational transitions stabilized by more than 5 kcal/mol (21 kJ/mol), which implies a large increase in rigidity for the reduced form.     

A water-lipid interface induces a highly dynamic folded state in apocytochrome c and cytochrome c, which may represent a common folding intermediate.

In this study, we have used CD and NMR techniques to investigate the secondary structure of (apo-) cytochrome c both in solution and when associated with micelles. In aqueous solution, the holoprotein cytochrome c is tightly folded at secondary and tertiary levels and differs strongly from its random-coiled precursor. However, in the presence of 12-PN/12-Pglycol (9:1) micelles, we observed a remarkable resemblance between the CD spectra of these partially helical proteins. The water-lipid interface induces a secondary folding of apocytochrome c, whereas cytochrome c is suggested to partially lose its tertiary structure. The exchange of all amide protons and, using deuterium-labeled proteins, of all amide deuterons with the solvent was monitored by NMR. A rapid exchange rate was observed, indicating that these folding states are highly dynamic. Saturation-transfer NMR of micelle-associated apocytochrome c showed that the exchange takes place at the (sub-) second time scale. The holoprotein in the presence of micelles was found to have two distinct exchange rates: (1) a fast rate, comparable to that found for the micelle-associated precursor and 4.5 times slower than that of the random-coiled apocytochrome c, and (2) a slow rate which is 75 times slower than the precursor in solution. Urea denaturation studies showed the micelle-bound proteins to have a low helix stability, which explains the inability of the lipid-induced secondary structure to prevent its labile protons from rapid exchange. The uniqueness of this lipid-induced highly dynamic folding state of (apo-) cytochrome c is demonstrated by comparison with amphiphilic polypeptides like melittin, and its implications for membrane translocation and functioning are discussed.     

Carbon-13 NMR method for the detection of correlated hydrogen exchange at adjacent backbone peptide amides and its application to hydrogen exchange in five antiparallel beta strands within the hydrophobic core of Streptomyces subtilisin inhibitor (SSI)

A novel method for monitoring proton-deuteron (H/D) exchange at backbone amides is based on the observation of H/D isotope effects on the (13)C NMR signals from peptide carbonyls. The line shape of the carbonyl (13)C(i) signal is influenced by differential H/D occupancy at the two adjacent amides: the H(N)(i)(+1) (beta site) and the H(N)(i) (gamma site). At a carbon frequency of 75.4 MHz, the H --> D isotope shifts on the (13)C signal are about 5-7 Hz for exchange at the beta site and 2 Hz or less for exchange at the gamma site. Because the effects at the two sites are additive, the time dependence of the line shape of a particular carbonyl resonance can report not only the exchange rates at the individual sites but also the level of dual exchange. Therefore, the data can be analyzed to determine the rate (k(c)) and degree of correlated exchange (X(betagamma)) at the two sites. We have applied this approach to the investigation of the pH dependence of hydrogen exchange at several adjacent residues in Streptomyces subtilisin inhibitor (SSI). Two selectively labeled SSI proteins were produced: one with selective (13)C' labeling at all valyl residues and one with selective (13)C' labeling at all leucyl residues. This permitted the direct observation by one-dimensional (13)C NMR of selected carbonyl signals from residues with slowly exchanging amides at the i and i + 1 positions. The residues investigated were located in an alpha helix and in a five-stranded antiparallel beta sheet. Samples of the two labeled proteins were prepared at various pH values, and (13)C NMR spectra were collected at 50 degrees C prior to and at various times after transferring the sample from H(2)O to (2)H(2)O. Most of the slowly exchanging amides studied were intramolecular hydrogen-bond donors. In agreement with prior studies, the results indicated that the exchange rates of the amide hydrogens in proteins are governed not only by hydrogen bonding but also by other factors. For example, the amide hydrogen of Thr34 exchanges rapidly even though it is an intramolecular hydrogen-bond donor. Over nearly the whole pH range studied, the apparent rates of uncorrelated exchange (k(beta) and k(gamma)) were proportional to [OH(-)] and the apparent rates of correlated exchange at two adjacent sites (k(c)) were roughly proportional to [OH(-)](2). This enabled us to extract the pH-independent exchange rates (k(beta) degrees , k(gamma) degrees , and k(c) degrees ). In all cases in which correlated exchange could be measured, the observed sigmoidal pH dependence of X(betagamma) could be replicated roughly from the derived pH-independent rates.     

Rapid Analysis of Protein Structure and Dynamics by Hydrogen/Deuterium Exchange Mass Spectrometry

An automated approach for the rapid analysis of protein structure has been developed and used to study acid-induced conformational changes in human growth hormone. The labeling approach involves hydrogen/deuterium exchange (H/D-Ex) of protein backbone amide hydrogens with rapid and sensitive detection by mass spectrometry (MS). Briefly, the protein is incubated for defined intervals in a deuterated environment. After rapid quenching of the exchange reaction, the partially deuterated protein is enzymatically digested and the resulting peptide fragments are analyzed by liquid chromatography mass spectrometry (LC-MS). The deuterium buildup curve measured for each fragment yields an average amide exchange rate that reflects the environment of the peptide in the intact protein. Additional analyses allow mapping of the free energy of folding on localized segments along the protein sequence affording unique dynamic and structural information. While amide H/D-Ex coupled with MS is recognized as a powerful technique for studying protein structure and protein–ligand interactions, it has remained a labor-intensive task. The improvements in the amide H/D-Ex methodology described here include solid phase proteolysis, automated liquid handling and sample preparation, and integrated data reduction software that together improve sequence coverage and resolution, while achieving a sample throughput nearly 10-fold higher than the commonly used manual methods.     

Secondary structure of filamentous bacteriophage coat protein is preserved in lipid environments

1H nuclear magnetic resonance experiments have shown that the amide hydrogens of residues 30 to 40 of bacteriophage Pf1 coat protein in micelles undergo very slow exchange with solvent deuterons. The amide 1H resonances from these residues were used to monitor the structural stability of the membrane-spanning helix of the coat protein during the transition of the coat protein from its structural form, in the virus particle, to the membrane-bound form, in micelles. The helix was found to remain folded on the 10(-3) second time-scale of the experiment, which indicates that no major disruption or rearrangement of the central part of the protein structure occurs during the process of coat protein solubilization by detergent. The results also suggest that a helical peptide can associate with lipids without reorganization of its secondary structure. However, a general model for the insertion of proteins into membranes cannot be established from these results, because the mechanism of the detergent solubilization process may differ somewhat from that of the membrane insertion process.     

Hydrogen exchange in BPTI variants that do not share a common disulfide bond.

Bovine pancreatic trypsin inhibitor (BPTI) is stabilized by 3 disulfide bonds, between cysteines 30-51, 5-55, and 14-38. To better understand the influence of disulfide bonds on local protein structure and dynamics, we have measured amide proton exchange rates in 2 folded variants of BPTI, [5-55]Ala and [30-51; 14-38]V5A55, which share no common disulfide bonds. These proteins resemble disulfide-bonded intermediates that accumulate in the BPTI folding pathway. Essentially the same amide hydrogens are protected from exchange in both of the BPTI variants studied here as in native BPTI, demonstrating that the variants adopt fully folded, native-like structures in solution. However, the most highly protected amide protons in each variant differ, and are contained within the sequences of previously studied peptide models of related BPTI folding intermediates containing either the 5-55 or the 30-51 disulfide bond.