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.     

Insights into enzyme structure and dynamics elucidated by amide H/D exchange mass spectrometry

Conformational changes and protein dynamics play an important role in the catalytic efficiency of enzymes. Amide H/D exchange mass spectrometry (H/D exchange MS) is emerging as an efficient technique to study the local and global changes in protein structure and dynamics due to ligand binding, protein activation-inactivation by modification, and protein-protein interactions. By monitoring the selective exchange of hydrogen for deuterium along a peptide backbone, this sensitive technique probes protein motions and structural elements that may be relevant to allostery and function. In this report, several applications of H/D exchange MS are presented which demonstrate the unique capability of amide hydrogen/deuterium exchange mass spectrometry for examining dynamic and structural changes associated with enzyme catalysis.     

Amide H/2H exchange reveals communication between the cAMP and catalytic subunit-binding sites in the R(I)alpha subunit of protein kinase A

The changes in backbone hydrogen/deuterium (H/2H) exchange in the regulatory subunit (R(I)alpha(94-244)) of cyclic AMP-dependent protein kinase A (PKA) were probed by MALDI-TOF mass spectrometry. The three naturally occurring states of the regulatory subunit were studied: (1) free R(I)alpha(94-244), which likely represents newly synthesized protein, (2) R(I)alpha(94-244) bound to the catalytic (C) subunit, or holoenzyme, and (3) R(I)alpha(94-244) bound to cAMP. Protection from amide exchange upon C-subunit binding was observed for the helical subdomain, including the A-helix and B-helix, pointing to regions adjacent to those shown to be important by mutagenesis. In addition, C-subunit binding caused changes in observed amide exchange in the distal cAMP-binding pocket. Conversely, cAMP binding caused protection in the cAMP-binding pocket and increased exchange in the helical subdomain. These results suggest that the mutually exclusive binding of either cAMP or C-subunit is controlled by binding at one site transmitting long distance changes to the other site.     

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.     

Dynamics of cAPK type IIbeta activation revealed by enhanced amide H/2H exchange mass spectrometry (DXMS)

cAMP-dependent protein kinase (cAPK) is a key component in numerous cell signaling pathways. The cAPK regulatory (R) subunit maintains the kinase in an inactive state until cAMP saturation of the R-subunit leads to activation of the enzyme. To delineate the conformational changes associated with cAPK activation, the amide hydrogen/deuterium exchange in the cAPK type IIbeta R-subunit was probed by electrospray mass spectrometry. Three states of the R-subunit, cAMP-bound, catalytic (C)-subunit bound, and apo, were incubated in deuterated water for various lengths of time and then, prior to mass spectrometry analysis, subjected to digestion by pepsin to localize the deuterium incorporation. High sequence coverage (>99%) by the pepsin-digested fragments enables us to monitor the dynamics of the whole protein. The effects of cAMP binding on RIIbeta amide hydrogen exchange are restricted to the cAMP-binding pockets, while the effects of C-subunit binding are evident across both cAMP-binding domains and the linker region. The decreased amide hydrogen exchange for residues 253-268 within cAMP binding domain A and for residues 102-115, which include the pseudosubstrate inhibitory site, support the prediction that these two regions represent the conserved primary and peripheral C-subunit binding sites. An increase in amide hydrogen exchange for a broad area within cAMP-binding domain B and a narrow area within cAMP-binding domain A (residues 222-232) suggest that C-subunit binding transmits long-distance conformational changes throughout the protein.     

Mass spectrometric measurement of protein amide hydrogen exchange rates of apo- and holo-myoglobin.

Measurement of backbone amide hydrogen exchange rates can provide detailed information concerning protein structure, dynamics, and interactions. Although nuclear magnetic resonance is typically used to provide these data, its use is restricted to lower molecular weight proteins that are soluble at millimolar concentrations. Not subject to these limitations is a mass spectrometric approach for measuring deuterium incorporation into proteins that are subsequently proteolyzed by pepsin; the resulting peptide masses are measured using a flowing-fast atom bombardment ionization source (Zhang Z, Smith DL, 1993, Protein Sci 2:522-531). In the current study, amide deuterium incorporation for intact apo- and holo-myoglobin was measured using liquid chromatography coupled directly to an electrospray ionization (LC/MS) source. Electrospray ionization provided a more complete coverage of the protein sequence and permitted the measurement of deuterium incorporation into intact proteins. Tandem mass spectrometry was used to rapidly identify the peptic peptides. It was found that within 30 s, the amides in apo-myoglobin were 47% deuterated, whereas holo-myoglobin was 12% deuterated. Peptic digestion and LC/MS demonstrated that regions represented by peptic peptides encompassing positions 1-7, 12-29, and 110-134 were not significantly altered by removal of the heme. Likewise, destabilized regions were identified within positions 33-106 and 138-153.     

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.     

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.     

Hydrogen bonding and equilibrium protium-deuterium fractionation factors in the immunoglobulin G binding domain of protein G

Protium-deuterium fractionation factors (phi) were determined for more than 85% of the backbone amide protons in the IgG binding domains of protein G, GB1 and GB2, from NMR spectra recorded over a range of H2O/D2O solvent ratios. Previous studies suggest a correlation between phi and hydrogen bond strength; amide and hydroxyl groups in strong hydrogen bonds accumulate protium (phi < 1), while weak hydrogen bonds accumulate deuterium (phi > 1). Our results show that the alpha-helical residues have slightly lower phi values (1.03 +/- 0.05) than beta-sheet residues (1.12 +/- 0.07), on average. The lowest phi value obtained (0.65) does not involve a backbone amide but rather is for the interaction between two side chains, Y45 and D47. Fractionation factors for solvent-exposed residues are between the alpha-helix and beta-sheet values, on average, and are close to those for random coil peptides. Further, the difference in phiav between alpha-helix and solvent-exposed residues is small, suggesting that differences in hydrogen bond strength for intrachain hydrogen bonds and amide...water hydrogen bonds are also small. Overall, the enrichment for deuterium suggests that most backbone...backbone hydrogen bonds are weak.     

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.     

Epitope mapping of a monoclonal antibody against human thrombin by H/D-exchange mass spectrometry reveals selection of a diverse sequence in a highly conserved protein

The epitope of a monoclonal antibody raised against human thrombin has been determined by hydrogen/deuterium exchange coupled to MALDI mass spectrometry. The antibody epitope was identified as the surface of thrombin that retained deuterium in the presence of the monoclonal antibody compared to control experiments in its absence. Covalent attachment of the antibody to protein G beads and efficient elution of the antigen after deuterium exchange afforded the analysis of all possible epitopes in a single MALDI mass spectrum. The epitope, which was discontinuous, consisting of two peptides close to anion-binding exosite I, was readily identified. The epitope overlapped with, but was not identical to, the thrombomodulin binding site, consistent with inhibition studies. The antibody bound specifically to human thrombin and not to murine or bovine thrombin, although these proteins share 86% identity with the human protein. Interestingly, the epitope turned out to be the more structured of two surface regions in which higher sequence variation between the three species is seen.     

Topographic study of arrestin using differential chemical modifications and hydrogen/deuterium exchange.

Arrestin is involved in the quenching of phototransduction by binding to photoactivated and phosphorylated rhodopsin (P-Rho*). To study its conformational changes and regions interacting with P-Rho*, arrestin was subjected to (1) differential acetylation at lysine residues in the presence and absence of P-Rho*, and (2) amide hydrogen/deuterium exchange. Labeled protein was proteolysed and analyzed by mass spectrometry. Three Lys residues, 28, 176, and 211, were protected from acetylation in native arrestin, although they were not located in regions exhibiting slow amide hydrogen exchange rates. The presence of P-Rho* protected lysine 201 from acetylation and partially protected 14 other lysyl residues, including (2, 5), (163, 166, 167), (232, 235, 236, 238), (267, 276), (298, 300), and 367, where parentheses indicate lysine residues found within the same peptide. In contrast, in the C-terminal region of arrestin, lysyl residues (386, 392, 395) were more exposed upon binding to P-Rho*. These data allowed us to identify functional regions in the arrestin molecule.     

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.     

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.     

Hydrogen/deuterium exchange studies of native rabbit MM-CK dynamics

Creatine kinase (CK) isoenzymes catalyse the reversible transfer of a phosphoryl group from ATP onto creatine. This reaction plays a very important role in the regulation of intracellular ATP concentrations in excitable tissues. CK isoenzymes are highly resistant to proteases in native conditions. To appreciate localized backbone dynamics, kinetics of amide hydrogen exchange with deuterium was measured by pulse-labeling the dimeric cytosolic muscle CK isoenzyme. Upon exchange, the protein was digested with pepsin, and the deuterium content of the resulting peptides was determined by liquid chromatography coupled to mass spectrometry (MS). The deuteration kinetics of 47 peptides identified by MS/MS and covering 96% of the CK backbone were analyzed. Four deuteration patterns have been recognized: The less deuterated peptides are located in the saddle-shaped core of CK, whereas most of the highly deuterated peptides are close to the surface and located around the entrance to the active site. Their exchange kinetics are discussed by comparison with the known secondary and tertiary structures of CK with the goal to reveal the conformational dynamics of the protein. Some of the observed dynamic motions may be linked to the conformational changes associated with substrate binding and catalytic mechanism.     

Distinct interaction modes of an AKAP bound to two regulatory subunit isoforms of protein kinase A revealed by amide hydrogen/deuterium exchange

The structure of an AKAP docked to the dimerization/docking (D/D) domain of the type II (RIIalpha) isoform of protein kinase A (PKA) has been well characterized, but there currently is no detailed structural information of an AKAP docked to the type I (RIalpha) isoform. Dual-specific AKAP2 (D-AKAP2) binds in the nanomolar range to both isoforms and provided us with an opportunity to characterize the isoform-selective nature of AKAP binding using a common docked ligand. Hydrogen/deuterium (H/D) exchange combined with mass spectrometry (DXMS) was used to probe backbone structural changes of an alpha-helical A-kinase binding (AKB) motif from D-AKAP2 docked to both RIalpha and RIIalpha D/D domains. The region of protection upon complex formation and the magnitude of protection from H/D exchange were determined for both interacting partners in each complex. The backbone of the AKB ligand was more protected when bound to RIalpha compared to RIIalpha, suggesting an increased helical stabilization of the docked AKB ligand. This combined with a broader region of backbone protection induced by the AKAP on the docking surface of RIalpha indicated that there were more binding constraints for the AKB ligand when bound to RIalpha. This was in contrast to RIIalpha, which has a preformed, localized binding surface. These distinct modes of AKAP binding may contribute to the more discriminating nature of the RIalpha AKAP-docking surface. DXMS provides valuable structural information for understanding binding specificity in the absence of a high-resolution structure, and can readily be applied to other protein-ligand and protein-protein interactions.     

The Deuterator: software for the determination of backbone amide deuterium levels from H/D exchange MS data

Background  

The combination of mass spectrometry and solution phase amide hydrogen/deuterium exchange (H/D exchange) experiments is an effective method for characterizing protein dynamics, and protein-protein or protein-ligand interactions. Despite methodological advancements and improvements in instrumentation and automation, data analysis and display remains a tedious process. The factors that contribute to this bottleneck are the large number of data points produced in a typical experiment, each requiring manual curation and validation, and then calculation of the level of backbone amide exchange. Tools have become available that address some of these issues, but lack sufficient integration, functionality, and accessibility required to address the needs of the H/D exchange community. To date there is no software for the analysis of H/D exchange data that comprehensively addresses these issues.     

Solution structure and dynamics of a heat shock protein assembly probed by hydrogen exchange and mass spectrometry

The solution conformation and dynamics of the 16.9 kDa small heat shock protein from wheat have been studied using a combination of hydrogen/deuterium exchange, proteolytic digestion, and mass spectrometry. At room temperature, HSP16.9 exists as a dodecameric assembly. Regions of HSP16.9 that form extensive and essential intersubunit contacts in the assembly, including residues 1-40 and 131-151, show little or no protection against hydrogen/deuterium exchange after incubation in D(2)O for 5 s. The high levels of hydrogen/deuterium exchange indicate that these regions have experienced large conformational fluctuations in solution, breaking intersubunit contacts and exposing buried amide hydrogens to solvent. When HSP16.9 is pulse labeled for 10 ms, residues 1-40 and 131-151 are substantially more protected than they are after 5 s. Thus, the breaking of intersubunit contacts occurs on a time scale between 10 milliseconds and 5 s. At 42 degrees C, HSP16.9 exists in a suboligomeric form. When the intrinsic temperature dependence of hydrogen/deuterium exchange is taken into account, exchange patterns at 25 and 42 degrees C are identical within experimental error, suggesting that the conformation of individual HSP16.9 subunits is the same in both the dodecameric and subdodecameric forms. Significant protection is seen in regions that form the dimeric interface, suggesting that the stable suboligomeric form is a dimer. Taken together, these results suggest that heat activation of HSP16.9 occurs by shifting the dodecamer <--> dimer equilibrium in favor of free dimers. The conformation of the dimers themselves does not appear to be altered with an increase in temperature.     

Solvent accessibility of the thrombin-thrombomodulin interface

The kinetics of solvent accessibility at the protein-protein interface between thrombin and a fragment of thrombomodulin, TMEGF45, have been monitored by amide hydrogen/deuterium (H/2H) exchange detected by MALDI-TOF mass spectrometry. The interaction is rapid and reversible, requiring development of theory and experimental methods to distinguish H/2H exchange due to solvent accessibility at the interface from H/2H exchange due to complex dissociation. Association and dissociation rate constants were measured by surface plasmon resonance and amide H/2H exchange rates were measured at different pH values and concentrations of TMEGF45. When essentially 100% of the thrombin was bound to TMEGF45, two segments of thrombin became completely solvent-inaccessible, as evidenced by the pH insensitivity of the amide H/2H exchange rates. These segments form part of anion-binding exosite I and contain the residues for which alanine substitution abolishes TM binding. Several other regions of thrombin showed slowing of amide exchange upon TMEGF45 binding, but the exchange remained pH-dependent, suggesting that these regions of thrombin were rendered only partially solvent-inaccessible by TMEGF45 binding. These partially inaccessible regions of thrombin form both surface and buried contacts into the active site of thrombin and contain residues implicated in allosteric changes in thrombin upon TM binding.     

Solution structure and backbone dynamics of long-[Arg(3)]insulin-like growth factor-I

Long-[Arg(3)]insulin-like growth factor-I (IGF-I) is a potent analog of insulin-like growth factor-I that has been modified by a Glu(3) --> Arg mutation and a 13-amino acid extension appended to the N terminus. We have determined the solution structure of (15)N-labeled Long-[Arg(3)]-IGF-I using high resolution NMR and restrained molecular dynamics techniques to a precision of 0.82 +/- 0.28 A root mean square deviation for the backbone heavy atoms in the three alpha-helices and 3.5 +/- 0.9 A root mean square deviation for all backbone heavy atoms excluding the 8 N-terminal residues and the 8 C-terminal eight residues. Overall, the structure of the IGF-I domain is consistent with earlier studies of IGF-I with some minor changes remote from the N terminus. The major variations in the structure, compared with IGF-I, occur at the N terminus with a substantial reorientation of the N-terminal three residues of the IGF-I domain. These results are interpreted in terms of the lower binding affinity for insulin-like growth factor-binding proteins. The backbone dynamics of Long-[Arg(3)]IGF-I were investigated using (15)N nuclear spin relaxation and the heteronuclear nuclear Overhauser enhancement (NOE). There is a considerable degree of flexibility in Long-[Arg(3)]IGF-I, even in the alpha-helices, as indicated by an average ((1)H)(15)N NOE of 0.55 for the regions. The largest heteronuclear NOEs are observed in the helical regions, lower heteronuclear NOEs are observed in the C-domain loop separating helix 1 from helix 2, and negative heteronuclear NOEs are observed in the N-terminal extension and at the C terminus. Despite these data indicating conformational flexibility for the N-terminal extension, slow amide proton exchange was observed for some residues in this region, suggesting some transitory structure does exist, possibly a molten helix. A certain degree of flexibility may be necessary in all insulin-like growth factors to enable association with various receptors and binding proteins.