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.     

Folding subdomains of thioredoxin characterized by native-state hydrogen exchange

Native-state hydrogen exchange (HX) studies, used in conjunction with NMR spectroscopy, have been carried out on Escherichia coli thioredoxin (Trx) for characterizing two folding subdomains of the protein. The backbone amide protons of only the slowest-exchanging 24 amino acid residues, of a total of 108 amino acid residues, could be followed at pH 7. The free energy of the opening event that results in an amide hydrogen exchanging with solvent (DeltaG(op)) was determined at each of the 24 amide hydrogen sites. The values of DeltaG(op) for the amide hydrogens belonging to residues in the helices alpha(1), alpha(2), and alpha(4) are consistent with them exchanging with the solvent only when the fully unfolded state is sampled transiently under native conditions. The denaturant-dependences of the values of DeltaG(op) provide very little evidence that the protein samples partially unfolded forms, lower in energy than the unfolded state. The amide hydrogens belonging to the residues in the beta strands, which form the core of the protein, appear to have higher values of DeltaG(op) than amide hydrogens belonging to residues in the helices, suggesting that they might be more stable to exchange. This apparently higher stability to HX of the beta strands might be either because they exchange out their amide hydrogens in a high energy intermediate preceding the globally unfolded state, or, more likely, because they form residual structure in the globally unfolded state. In either case, the central beta strands-beta(3,) beta(2), and beta(4)-would appear to form a cooperatively folding subunit of the protein. The native-state HX methodology has made it possible to characterize the free energy landscape that Trx can sample under equilibrium native conditions.     

Deuterium exchange and mass spectrometry reveal the interaction differences of two synthetic modulators of RXRalpha LBD

Protein amide hydrogen/deuterium (H/D) exchange was used to compare the interactions of two antagonists, UVI 2112 and UVI 3003, with that of the agonist, 9-cis-retinoic acid, upon binding to the human retinoid X receptor alpha ligand-binding domain (hRXRalpha LBD) homodimer. Analysis of the H/D content by mass spectrometry showed that in comparison to 9-cis-retinoic acid, the antagonists provide much greater protection toward deuterium exchange-in throughout the protein, suggesting that the protein-antagonist complex adopts a more restricted conformation or ensemble of conformations in which solvent accesses to amide protons are reduced. A comparison between the two antagonists shows that UVI 3003 is more protective in the C-terminal region due to the extra hydrophobic interactions derived from the atoms in the benzene ring of the carboxylic acid chain. It was less protective within regions comprising peptides 271-278 and 326-330 due to differences in conformational orientation, and/or shorter carboxylic acid chain length. Decreased deuterium exchange-in in the segment 234-239 where the residues do not involve interactions with the ligand was observed with the two antagonists, but not with 9-cis-RA. The amide protons of helix 12 of the agonist- or antagonist-occupied protein in solution have the same deuterium exchange rates as the unliganded protein, supporting a suggestion made previously that helix 12 can cover the occupied binding cavity only with the cofactor present to adjust its location.     

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.     

Photoactivation of rhodopsin causes an increased hydrogen-deuterium exchange of buried peptide groups.

A key step in visual transduction is the light-induced conformational changes of rhodopsin that lead to binding and activation of the G-protein transducin. In order to explore the nature of these conformational changes, time-resolved Fourier transform infrared spectroscopy was used to measure the kinetics of hydrogen/deuterium exchange in rhodopsin upon photoexcitation. The extent of hydrogen/deuterium exchange of backbone peptide groups can be monitored by measuring the integrated intensity of the amide II and amide II' bands. When rhodopsin films are exposed to D2O in the dark for long periods, the amide II band retains at least 60% of its integrated intensity, reflecting a core of backbone peptide groups that are resistant to H/D exchange. Upon photoactivation, rhodopsin in the presence of D2O exhibits a new phase of H/D exchange which at 10 degrees C consists of fast (time constant approximately 30 min) and slow (approximately 11 h) components. These results indicate that photoactivation causes buried portions of the rhodopsin backbone structure to become more accessible.     

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).     

Hydrogen/deuterium exchange and aggregation of a polyvaline and a polyleucine alpha-helix investigated by matrix-assisted laser desorption ionization mass spectrometry

The membrane-associated pulmonary surfactant protein C (SP-C), containing a polyvaline alpha-helix, and a synthetic SP-C analogue with a polyleucine helix (SP-C(Leu)) were studied by hydrogen/deuterium exchange matrix-assisted laser desorption ionization (MALDI) mass spectrometry. SP-C, but not SP-C(Leu), formed abundant amyloid fibrils under experimental conditions. In CD(3)OD/D(2)O, 91:9 (v/v), containing 2 mM ammonium acetate, SP-C(Leu) and SP-C exchanged 40% of their exchangeable hydrogens within 1 min. This corresponds to exchange of labile side-chain hydrogen atoms, hydrogens on the N- and C-terminal heteroatoms, and amide hydrogen atoms in the unstructured N-terminal regions. After approximately 300 h, four exchangeable hydrogen atoms in SP-C(Leu) and 10 in SP-C remained unexchanged. During this time period the ion current corresponding to singly charged SP-C decreased to <10% of the initial value due to the formation of insoluble aggregates that are not detected by MALDI mass spectrometry. In contrast, the ion current for SP-C(Leu) was maintained over this time period, although the peptides were incubated together. In combination, hydrogen/deuterium exchange and aggregation data indicate that the polyleucine peptide refolds into a helix after opening, while the unfolded polyvaline peptide forms insoluble beta-sheet aggregates rather than refolding into a helix. The SP-C helix, but not the SP-C(Leu) helix, is thus in a metastable state, which may contribute to the recently observed tendency of SP-C and its precursor to misfold and aggregate in vivo.     

Comparison of hydrogen exchange rates for bovine pancreatic trypsin inhibitor in crystals and in solution.

Hydrogen exchange rate constants for the 17 slowest exchanging amide NH groups in bovine pancreatic trypsin inhibitor (BPTI) were measured in solution and in form II and form III crystals. All 17 amide hydrogens are buried and intramolecularly hydrogen bonded in the crystal structure, except Lys 41 which is buried and hydrogen bonded to a buried water. Large-scale crystallization procedures were developed for these experiments, and rate constants for both crystal and solution exchange were measured by 1H NMR spectroscopy of exchange-quenched samples in solution. Two conditions of pH and temperature, pH 9.8 and 35 degrees C, and pH 9.4 and 25 degrees C, bring two groups of hydrogens into the experimental time window (minutes to weeks). One consists of the 10 slowest exchanging hydrogens, all of which are associated with the central beta-sheet of BPTI. The second group consists of seven more rapidly exchanging hydrogens, which are distributed throughout the molecule, primarily in a loop or turn. In both groups, most hydrogens exchange more slowly in crystals, but there is considerable variation in the degree to which the exchange is depressed in crystals. Many differences observed for the more rapidly exchanging hydrogens can be attributed to local surface effects arising from intermolecular contacts in the crystal lattice. Within the slower group, however, a very large effect on exchange of Ile 18 and Tyr 35 appears to be selectively transmitted through the matrix of the molecule.(ABSTRACT TRUNCATED AT 250 WORDS)     

The secondary structure of the inhibited mitochondrial ADP/ATP transporter from yeast analyzed by FTIR spectroscopy

Fourier transform infrared spectroscopy has been applied to the study of the carboxyatractyloside-inhibited mitochondrial ADP/ATP transporter from the yeast Saccharomyces cerevisiae, either solubilized in dodecyl maltoside or reconstituted in phosphatidylcholine liposomes. Its secondary structure has been estimated by means of Fourier self-deconvolution followed by curve fit. A Voigt function was used to fit the components of the deconvoluted spectrum, aiming to account for any distortions introduced by deconvolution. For any of the states analyzed, reconstituted or solubilized, in solution or in dry films, 60-70% of the amino acids are found to adopt alpha-helix plus unordered structures, coherent with the six transmembrane spanning helix model. Moreover, the problem of structure preservation on drying was addressed, and several observations pointed to a maintenance of the protein structure in dry films. Comparison of reconstituted and solubilized samples indicated the presence of both lipid-induced changes in the protein (decrease of the beta-sheets and increase of unordered structures) and protein-induced changes in the lipids (strong hydrogen bonding of lipid C=O groups). To obtain a better discrimination of alpha-helix and unordered structure contributions for the reconstituted form, H/D exchange experiments were performed. Between 35% and 45% of the amino acids were finally assigned to alpha-helix structures, compatible with the existence of five or six transmembrane spanning helices in the transporter. The level of H/D exchange was determined after 15 h of exposure to D(2)O vapor to be 85%, reflecting a high accessibility of the amide hydrogens even for the carboxyatractyloside-inhibited state.     

Mass spectrometry on hydrogen/deuterium exchange of dihydrofolate reductase: effects of ligand binding

To address the effects of ligand binding on the structural fluctuations of Escherichia coli dihydrofolate reductase (DHFR), the hydrogen/deuterium (H/D) exchange kinetics of its binary and ternary complexes formed with various ligands (folate, dihydrofolate, tetrahydrofolate, NADPH, NADP(+), and methotrexate) were examined using electrospray ionization mass spectrometry. The kinetic parameters of H/D exchange reactions, which consisted of two phases with fast and slow rates, were sensitively influenced by ligand binding, indicating that changes in the structural fluctuation of the DHFR molecule are associated with the alternating binding and release of the cofactor and substrate. No additivity was observed in the kinetic parameters between a ternary complex and its constitutive binary complexes, indicating that ligand binding cooperatively affects the structural fluctuation of the DHFR molecule via long-range interactions. The local H/D exchange profile of pepsin digestion fragments was determined by matrix-assisted laser desorption/ionization mass spectrometry, and the helix and loop regions that appear to participate in substrate binding, largely fluctuating in the apo-form, are dominantly influenced by ligand binding. These results demonstrate that the structural fluctuation of kinetic intermediates plays an important role in enzyme function, and that mass spectrometry on H/D exchange coupled with ligand binding and protease digestion provide new insight into the structure-fluctuation-function relationship of enzymes.     

NMR characterization of three-disulfide variants of lysozyme, C64A/C80A, C76A/C94A, and C30A/C115A--a marginally stable state in folded proteins

Our earlier NMR study showed that a two-disulfide variant of hen lysozyme containing intra-alpha-domain disulfide bridges, C6-C127 and C30-C115, is partially folded, with the alpha domain tightly folded to the nativelike conformation and the beta domain flexible or unfolded. With a view that the formation of a third disulfide bridge is a key for the accomplishment of the overall chain fold, three-dimensional structures of three-disulfide variants of hen lysozyme lacking one disulfide bridge (C64A/C80A, C76A/C94A, and C30A/C115A) were studied in detail using NMR spectroscopy. Amide hydrogen exchange rates were measured to estimate the degree of conformational fluctuation in a residue-specific manner. The structure of C76A/C94A was found to be quite similar to that of the wild type, except for the peptide segment of residues 74-78. The structure of C64A/C80A was considerably disordered in the entire region of the loop (residues 62-79). Further, it was found that a network of hydrogen bonds within the beta sheet and the 3(10) helix in the beta domain were disrupted and fluctuating. In C30A/C115A, the D helix was unstructured and the interface of the B helix with the D helix was significantly perturbed. However, the structural disorder generated in the hydrophobic core of the alpha domain was prevented by the C helix from propagating toward the beta domain. A marginally stable state in folded proteins is discussed based on the structures remaining in each three-disulfide variant.     

The kinetic and equilibrium molten globule intermediates of apoleghemoglobin differ in structure

An important question in protein folding is whether molten globule states formed under equilibrium conditions are good structural models for kinetic folding intermediates. The structures of the kinetic and equilibrium intermediates in the folding of the plant globin apoleghemoglobin have been compared at high resolution by quench-flow pH-pulse labeling and interrupted hydrogen/deuterium exchange analyzed in dimethyl sulfoxide. Unlike its well studied homolog apomyoglobin, where the equilibrium and kinetic intermediates are quite similar, there are striking structural differences between the intermediates formed by apoleghemoglobin. In the kinetic intermediate, formed during the burst phase of the quench-flow experiment, protected amides and helical structure are found mainly in the regions corresponding to the G and H helices of the folded protein, and in parts of the E helix and CE loop regions, whereas in the equilibrium intermediate, amide protection and helical structure are seen in parts of the A and B helix regions, as well as in the G and H regions, and the E helix remains largely unfolded. These results suggest that the structure of the molten globule intermediate of apoleghemoglobin is more plastic than that of apomyoglobin, so that it is readily transformed depending on the solution conditions, particularly pH. Thus, in the case of apoleghemoglobin at least, the equilibrium molten globule formed under destabilizing conditions at acid pH is not a good model for the compact intermediate formed during kinetic refolding experiments. Our high-precision kinetic analysis also reveals an additional slow phase during the folding of apoleghemoglobin, which is not observed for apomyoglobin. Hydrogen exchange pulse-labeling experiments show that the slow-folding phase is associated with residues in the CE loop, which probably forms non-native structure in the intermediate that must be resolved before folding can proceed to completion.     

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%.     

Structural changes in the C-terminus of Ca2+-bound rat S100B (beta beta) upon binding to a peptide derived from the C-terminal regulatory domain of p53.

S100B(beta beta) is a dimeric Ca2+-binding protein that interacts with p53, inhibits its phosphorylation by protein kinase C (PKC) and promotes disassembly of the p53 tetramer. Likewise, a 22 residue peptide derived from the C-terminal regulatory domain of p53 has been shown to interact with S100B(beta beta) in a Ca2+-dependent manner and inhibits its phosphorylation by PKC. Hence, structural studies of Ca2+-loaded S100B(beta beta) bound to the p53 peptide were initiated to characterize this interaction. Analysis of nuclear Overhauser effect (NOE) correlations, amide proton exchange rates, 3J(NH-H alpha) coupling constants, and chemical shift index data show that, like apo- and Ca2+-bound S100B(beta beta), S100B remains a dimer in the p53 peptide complex, and each subunit has four helices (helix 1, Glu2-Arg20; helix 2, Lys29-Asn38; helix 3, Gln50-Asp61; helix 4, Phe70-Phe87), four loops (loop 1, Glu21-His25; loop 2, Glu39-Glu49; loop 3, Glu62-Gly66; loop 4, Phe88-Glu91), and two beta-strands (beta-strand 1, Lys26-Lys28; beta-strand 2, Glu67-Asp69), which forms a short antiparallel beta-sheet. However, in the presence of the p53 peptide helix 4 is longer by five residues than in apo- or Ca2+-bound S100B(beta beta). Furthermore, the amide proton exchange rates in helix 3 (K55, V56, E58, T59, L60, D61) are significantly slower than those of Ca2+-bound S100B(beta beta). Together, these observations plus intermolecular NOE correlations between the p53 peptide and S100B(beta beta) support the notion that the p53 peptide binds in a region of S100B(beta beta), which includes residues in helix 2, helix 3, loop 2, and the C-terminal loop, and that binding of the p53 peptide interacts with and induces the extension of helix 4.     

Determination of amide hydrogen exchange by mass spectrometry: a new tool for protein structure elucidation.

A new method based on protein fragmentation and directly coupled microbore high-performance liquid chromatography-fast atom bombardment mass spectrometry (HPLC-FABMS) is described for determining the rates at which peptide amide hydrogens in proteins undergo isotopic exchange. Horse heart cytochrome c was incubated in D2O as a function of time and temperature to effect isotopic exchange, transferred into slow exchange conditions (pH 2-3, 0 degrees C), and fragmented with pepsin. The number of peptide amide deuterons present in the proteolytic peptides was deduced from their molecular weights, which were determined following analysis of the digest by HPLC-FABMS. The present results demonstrate that the exchange rates of amide hydrogens in cytochrome c range from very rapid (k > 140 h-1) to very slow (k < 0.002 h-1). The deuterium content of specific segments of the protein was determined as a function of incubation temperature and used to indicate participation of these segments in conformational changes associated with heating of cytochrome c. For the present HPLC-FABMS system, approximately 5 nmol of protein were used for each determination. Results of this investigation indicate that the combination of protein fragmentation and HPLC-FABMS is relatively free of constraints associated with other analytical methods used for this purpose and may be a general method for determining hydrogen exchange rates in specific segments of proteins.     

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.     

Fourier transform infrared studies of ribonuclease in H2O and 2H2O solutions

Fourier transform infrared transmission spectra have been obtained of the enzyme ribonuclease in both H₂O and ²H₂O. The resolution of the spectra have been enhanced by Fourier self-deconvolution procedures. The infrared spectrum of ribonuclease changes during exchange of the enzyme's amide hydrogens for deuterium and the exchange has been followed in the amide I and amide II spectral regions. The amide I band shifts towards lower wavenumbers during both the fast and slow phases of hydrogen exchange and the interpretation of these shifts has aided the band assignments. In particular these studies have allowed an assignment to be made for the high frequency component of the β-strand absorption that differs from that proposed previously. This paper represents the first example of the use of deconvoluted Fourier transform infrared spectra in conjunction with hydrogen-deuterium exchange in order to aid in the assignment of a proteins's infrared bands.     

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.