Hydra: software for tailored processing of H/D exchange data from MS or tandem MS analyses

Background  

Hydrogen/deuterium exchange mass spectrometry (H/DX-MS) experiments implemented to characterize protein interaction and protein folding generate large quantities of data. Organizing, processing and visualizing data requires an automated solution, particularly when accommodating new tandem mass spectrometry modes for H/DX measurement. We sought to develop software that offers flexibility in defining workflows so as to support exploratory treatments of H/DX-MS data, with a particular focus on the analysis of very large protein systems and the mining of tandem mass spectrometry data.     

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.     

Anion modulation of the 1H/2H exchange rates in backbone amide protons monitored by NMR spectroscopy

The ability of three anionic cosolutes (sulfate, thiocyanate, and chloride) in modulating the (1)H/(2)H exchange rates for backbone amide protons has been investigated using nuclear magnetic resonance (NMR) for two different proteins: the IGg-binding domain of protein L (ProtL) and the glucose-galactose-binding protein (GGBP). Our results show that moderate anion concentrations (0.2 M-1 M) regulate the exchange rate following the Hofmeister series: Addition of thiocyanate increases the exchange rates for both proteins, while sulfate and chloride (to a less extent) slow down the exchange reaction. In the presence of the salt, no alteration of the protein structure and minimal variations in the number of measurable peaks are observed. Experiments with model compounds revealed that the unfolded state is modulated in an equivalent way by these cosolutes. For ProtL, the estimated values for the local free energy change upon salt addition (m (3,DeltaG )) are consistent with the previously reported free energy contribution from the cosolute's preferential interaction/exclusion term indicating that nonspecific weak interactions between the anion and the amide groups constitute the dominant mechanism for the exchange-rate modulation. The same trend is also found for GGBP in the presence of thiocyanate, underlining the generality of the exchange-rate modulation mechanism, complementary to more investigated effects like the electrostatic interactions or specific anion binding to protein sites.     

Polymorphic triple beta-sheet structures contribute to amide hydrogen/deuterium (H/D) exchange protection in the Alzheimer amyloid beta42 peptide

Characterization of the polymorphic structural range of Aβ oligomers is important to the understanding of the mechanisms of toxicity. Yet for highly polymorphic ensembles, experimental structural elucidation is difficult. Here, we use a combination of NMR solvent protection experiments and computational structural screening to identify major species in the amyloid conformational ensemble. We examined the polymorphic pentamer and fibril seeds of Aβ42 and its mutants and compared the theoretical backbone amide protection obtained from simulations with experimental hydrogen/deuterium (H/D) exchange protection ratio. We observed that highly flexible pentamers do not share structural similarities with fibril seed oligomers, except the turn regions. We found that a novel amyloid structural motif of a triple β-sheet, with the N-terminal residues interacting with the core (Lys(17)-Glu(22)) β-sheet region, correlates with H/D exchange protection. The triple β-sheet Aβ42 oligomer has a minimal exposure of hydrophobic residues and is further stabilized by the E22Q (Dutch) mutation in Alzheimer disease. The experimental H/D exchange solvent protection ratio implies that triple β-sheet fibrils and globulomers could coexist in the Aβ42 ensemble, pointing to a broad heterogeneous aggregate population. Our results suggest that an approach that combines computational modeling with NMR protection data can be a useful strategy for obtaining clues to the preferred conformational species of the assemblies in solution and help in alleviating experimental difficulties and consequently possible errors in the exchange data for Aβ42 fibrils.     

Conformational analysis of Epac activation using amide hydrogen/deuterium exchange mass spectrometry

Exchange proteins directly activated by cAMP (Epac) play important roles in mediating the effects of cAMP through the activation of downstream small GTPases, Rap. To delineate the mechanism of Epac activation, we probed the conformation and structural dynamics of Epac using amide hydrogen/deuterium exchange and structural modeling. Our studies show that cAMP induces significant conformational changes that lead to a spatial rearrangement of the regulatory components of Epac and allows the exposure of the catalytic core for effector binding without imposing significant conformational change on the catalytic core. Homology modeling and comparative structural analyses of the cAMP binding domains of Epac and cAMP-dependent protein kinase (PKA) lead to a model of Epac activation, in which Epac and PKA activation by cAMP employs the same underlying principle, although the detailed structural and conformational changes associated with Epac and PKA activation are significantly different.     

Construct optimization for protein NMR structure analysis using amide hydrogen/deuterium exchange mass spectrometry

Disordered or unstructured regions of proteins, while often very important biologically, can pose significant challenges for resonance assignment and three‐dimensional structure determination of the ordered regions of proteins by NMR methods. In this article, we demonstrate the application of ¹H/²H exchange mass spectrometry (DXMS) for the rapid identification of disordered segments of proteins and design of protein constructs that are more suitable for structural analysis by NMR. In this benchmark study, DXMS is applied to five NMR protein targets chosen from the Northeast Structural Genomics project. These data were then used to design optimized constructs for three partially disordered proteins. Truncated proteins obtained by deletion of disordered N‐ and C‐terminal tails were evaluated using ¹H‐¹⁵N HSQC and ¹H‐¹⁵N heteronuclear NOE NMR experiments to assess their structural integrity. These constructs provide significantly improved NMR spectra, with minimal structural perturbations to the ordered regions of the protein structure. As a representative example, we compare the solution structures of the full length and DXMS‐based truncated construct for a 77‐residue partially disordered DUF896 family protein YnzC from Bacillus subtilis, where deletion of the disordered residues (ca. 40% of the protein) does not affect the native structure. In addition, we demonstrate that throughput of the DXMS process can be increased by analyzing mixtures of up to four proteins without reducing the sequence coverage for each protein. Our results demonstrate that DXMS can serve as a central component of a process for optimizing protein constructs for NMR structure determination. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

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.     

AMASS: software for automatically validating the quality of MS/MS spectrum from SEQUEST results

Time-consuming and experience-dependent manual validations of tandem mass spectra are usually applied to SEQUEST results. This inefficient method has become a significant bottleneck for MS/MS data processing. Here we introduce a program AMASS (advanced mass spectrum screener), which can filter the tandem mass spectra of SEQUEST results by measuring the match percentage of high-abundant ions and the continuity of matched fragment ions in b, y series. Compared with Xcorr and DeltaCn filter, AMASS can increase the number of positives and reduce the number of negatives in 22 datasets generated from 18 known protein mixtures. It effectively removed most noisy spectra, false interpretations, and about half of poor fragmentation spectra, and AMASS can work synergistically with Rscore filter. We believe the use of AMASS and Rscore can result in a more accurate identification of peptide MS/MS spectra and reduce the time and energy for manual validation.     

Ca2+ binding to calbindin D9k strongly affects backbone dynamics: measurements of exchange rates of individual amide protons using 1H NMR

One- and two-dimensional 1H NMR have been used to study the backbone dynamics in Ca2(+)-free (apo) and Ca2(+)-loaded (Ca2) calbindin D9k at pH 7.5 and 25 degrees C. Hydrogen exchange rates of all 71 backbone amide protons (NH's) have been measured for the Ca2 form by both a direct exchange-out experiment and another experiment that measures the transfer of saturation from water protons to amide protons. A large number of NH's are found to be highly protected against exchange with solvent protons. The results for the Ca2 form are related to solvent accessibility and hydrogen bonding obtained in molecular dynamics simulations of calcium-loaded calbindin. The correlation with these parameters is strong within the N-terminal half of calbindin, which is found to be more stable than the C-terminal half. The amide proton exchange in the apo form is much faster than in the Ca2 form and was studied in a series of experiments in which the exchange was quenched after different times by Ca2+ addition. This experiment is applicable to all amide hydrogens that exchange slowly in the Ca2 form. For these NH's the effects of Ca2+ removal span from a 10(2)-fold decrease to a 10(5)-fold increase of the exchange rate, and the average is a 220-fold increase. The effects on individual NH exchange rates show that the four alpha-helices are almost intact after calcium removal and that the changes in dynamics involve not only the Ca2(+)-binding region. Hydrogen bonds involving backbone NH's in the Ca2+ loops appear to be broken or weakened when calbindin releases Ca2+, whereas the beta-sheet between the Ca2+ loops is found to be present in both the Ca2 and apo forms. Large Ca2(+)-induced effects on NH exchange rates were measured for a few residues at alpha-helix ends far from the two Ca2(+)-binding sites. This may be the result of a change in interhelix angles (or the rate of interhelix angle fluctuations) on calcium binding.     

Mapping of discontinuous conformational epitopes by amide hydrogen/deuterium exchange mass spectrometry and computational docking

Understanding antigen-antibody interactions at the sub-molecular level is of particular interest for scientific, regulatory, and intellectual property reasons, especially with increasing demand for monoclonal antibody therapeutic agents. Although various techniques are available for the determination of an epitope, there is no widely applicable, high-resolution, and reliable method available. Here, a combination approach using amide hydrogen/deuterium exchange coupled with proteolysis and mass spectrometry (HDX-MS) and computational docking was applied to investigate antigen-antibody interactions. HDX-MS is a widely applicable, medium-resolution, medium-throughput technology that can be applied to epitope identification. First, the epitopes of cytochrome c-E8, IL-13-CNTO607, and IL-17A-CAT-2200 interactions identified using the HDX-MS method were compared with those identified by X-ray co-crystal structures. The identified epitopes are in good agreement with those identified using high-resolution X-ray crystallography. Second, the HDX-MS data were used as constraints for computational docking. More specifically, the non-epitope residues of an antigen identified using HDX-MS were designated as binding ineligible during computational docking. This approach, termed HDX-DOCK, gave more tightly clustered docking poses than stand-alone docking for all antigen-antibody interactions examined and improved docking results significantly for the cytochrome c-E8 interaction.     

Crystal structure and amide H/D exchange of binary complexes of alcohol dehydrogenase from Bacillus stearothermophilus: insight into thermostability and cofactor binding

The crystal structure of NAD(+)-dependent alcohol dehydrogenase from Bacillus stearothermophilus strain LLD-R (htADH) was determined using X-ray diffraction data at a resolution of 2.35 A. The structure of homotetrameric htADH is highly homologous to those of bacterial and archaeal homotetrameric alcohol dehydrogenases (ADHs) and also to the mammalian dimeric ADHs. There is one catalytic zinc atom and one structural zinc atom per enzyme subunit. The enzyme was crystallized as a binary complex lacking the nicotinamide adenine dinucleotide (NAD(+)) cofactor but including a zinc-coordinated substrate analogue trifluoroethanol. The binary complex structure is in an open conformation similar to ADH structures without the bound cofactor. Features important for the thermostability of htADH are suggested by a comparison with a homologous mesophilic enzyme (55% identity), NAD(+)-dependent alcohol dehydrogenase from Escherichia coli. To gain insight into the conformational change triggered by NAD(+) binding, amide hydrogen-deuterium exchange of htADH, in the presence and absence of NAD(+), was studied by HPLC-coupled electrospray mass spectrometry. When the deuteron incorporation of the protein-derived peptides was analyzed, it was found that 9 of 21 peptides show some decrease in the level of deuteron incorporation upon NAD(+) binding, and another 4 peptides display slower exchange rates. With one exception (peptide number 8), none of the peptides that are altered by bound NAD(+) are in contact with the alcohol-substrate-binding pocket. Furthermore, peptides 5 and 8, which are located outside the NAD(+)-binding pocket, are notable by displaying changes upon NAD(+) binding. This suggests that the transition from the open to the closed conformation caused by cofactor binding has some long-range effects on the protein structure and dynamics.     

Analysis of conformational changes during activation of protein kinase Pak2 by amide hydrogen/deuterium exchange

During apoptotic stress, protein kinase Pak2 is cleaved by caspase 3 to form a heterotetramer that is constitutively activated following autophosphorylation. The active protein kinase migrates slightly slower than the inactive holoenzyme when analyzed by gel filtration, suggesting an expanded conformation. Activation of Pak2 comprises a series of structural changes resulting from caspase cleavage, ATP binding, and autophosphorylation of Pak2. Changes at each step were individually analyzed by amide hydrogen/deuterium exchange coupled with mass spectrometry and compared with inactive Pak2. The auto-inhibited form was shown to bind ATP in the active site, with minor changes in the glycine loop and the autoinhibitory domain (AID). Caspase cleavage produced significant changes in solvent accessibility in the AID and upper lobe of the catalytic domain. Cleavage of ATP-bound Pak2 relaxes the allosteric inhibition, as shown by increased solvent accessibility in the upper and lower lobes, including the G-helix, facilitating the autophosphorylation of two sites required for activation, Ser-141 in the regulatory domain and Thr-402 in the catalytic domain. Autophosphorylation increased the amide hydrogen/deuterium exchange solvent accessibility of the contact region between the AID and the G-helix, the E-F loop, and the N terminus. Thus, activation of Pak2 via caspase cleavage is associated with structural relaxation of Pak2 that allows for complete auto-phosphorylation, resulting in a more comprehensive solvent-exposed and conformationally dynamic enzyme.     

Influence of the dimer interface on glutathione transferase structure and dynamics revealed by amide H/D exchange mass spectrometry

Mammalian glutathione (GSH) transferases are dimeric proteins, many of which share a common hydrophobic interaction motif that is important for dimer stability. In the rGSTM1-1 enzyme this motif involves the side chain of F56, located on the 56 loop of the N-terminal domain, which is intercalated between the alpha4- and alpha5-helices of the C-terminal domain of the opposing subnuit. Disruption of the complementary interactions in this motif by mutation of F56 to serine, arginine, or glutamate is known to have deleterious effects on catalytic efficiency but remarkably different effects on the stability of the dimer [Hornby et al. (2002) Biochemistry 41, 14238-14247]. The structural basis for the behavior of the mutants by amide H/D exchange mass spectrometry is described. A substantial decrease in H/D exchange is observed in the GSH binding domain and in parts of the dimer interface upon substrate binding. The F56S and F56R mutants exhibit enhanced H/D exchange kinetics in the GSH binding domain and at the dimer interface. In contrast, the F56E mutant shows a decrease in the rate and extent of amide H/D exchange at the dimer interface and enhanced exchange kinetics in the GSH binding domain. The results suggest that the F56E mutant has a restructured dimer interface with decreased solvent accessibility and dynamics. Although all of the F56 mutations disrupt the GSH binding site, the effects of the mutations on the structure of the subunit interface and dimer stability are quite distinct.     

Gaining knowledge from previously unexplained spectra-application of the PTM-Explorer software to detect PTM in HUPO BPP MS/MS data

A novel software tool named PTM-Explorer has been applied to LC-MS/MS datasets acquired within the Human Proteome Organisation (HUPO) Brain Proteome Project (BPP). PTM-Explorer enables automatic identification of peptide MS/MS spectra that were not explained in typical sequence database searches. The main focus was detection of PTMs, but PTM-Explorer detects also unspecific peptide cleavage, mass measurement errors, experimental modifications, amino acid substitutions, transpeptidation products and unknown mass shifts. To avoid a combinatorial problem the search is restricted to a set of selected protein sequences, which stem from previous protein identifications using a common sequence database search. Prior to application to the HUPO BPP data, PTM-Explorer was evaluated on excellently manually characterized and evaluated LC-MS/MS data sets from Alpha-A-Crystallin gel spots obtained from mouse eye lens. Besides various PTMs including phosphorylation, a wealth of experimental modifications and unspecific cleavage products were successfully detected, completing the primary structure information of the measured proteins. Our results indicate that a large amount of MS/MS spectra that currently remain unidentified in standard database searches contain valuable information that can only be elucidated using suitable software tools.     

MassSieve: Panning MS/MS peptide data for proteins

We present MassSieve, a Java‐based platform for visualization and parsimony analysis of single and comparative LC‐MS/MS database search engine results. The success of mass spectrometric peptide sequence assignment algorithms has led to the need for a tool to merge and evaluate the increasing data set sizes that result from LC‐MS/MS‐based shotgun proteomic experiments. MassSieve supports reports from multiple search engines with differing search characteristics, which can increase peptide sequence coverage and/or identify conflicting or ambiguous spectral assignments.     

A strategy to obtain backbone resonance assignments of deuterated proteins in the presence of incomplete amide 2H/1H back-exchange

Replacement of non-exchangeable protons by deuterons has become a standard tool in structural studies of proteins on the order of 30–40 kDa to overcome problems arising from rapid ¹H and ¹³C transverse relaxation. However, ¹H nuclei are required at exchangeable sites to maintain the benefits of proton detection. Protein expression in D₂O-based media containing deuterated carbon sources yields protein deuterated in all positions. Subsequent D/H-exchange is commonly used to reintroduce protons in labile positions. Since this strategy may fail for large proteins with strongly inhibited exchange we propose to express the protein in fully deuterated algal lysate medium in 100% H₂O. As a side-effect partial C protonation occurs in a residue-type dependent manner. Samples obtained by this protocol are suitable for complementary ¹H^N- and ¹H-based triple resonance experiments allowing complete backbone resonance assignments in cases where back-exchange of amide protons is very slow after expression in D₂O and refolding of chemically denatured protein is not feasible. This approach is explored using a 35-kDa protein as a test case. The degree of C protonation of individual amino acids is determined quantitatively and transverse relaxation properties of ¹H^N and ¹⁵N nuclei of the partially deuterated protein are investigated and compared to the fully protonated and perdeuterated species. Based on the deviations of assigned chemical shifts from random coil values its solution secondary structure can be established.     

High-sensitivity mass spectrometry for imaging subunit interactions: hydrogen/deuterium exchange

In recent years, advances in mass spectrometry have provided unprecedented knowledge of protein expression within cells. It has become apparent that many proteins function as macromolecular complexes. Structural genomics programs are determining the fold of these proteins at an increasing rate and electron microscopic tomography potentially provides a means to determine the location of these complexes within the cell. A complete understanding of the molecular mechanism of these proteins requires detailed information on the interactions and dynamics within the complex. Recent advances in mass spectrometry now make it possible to use hydrogen/deuterium exchange to detect intersubunit interfaces and dynamics within supramolecular complexes.