- Open Access
QUDeX-MS: hydrogen/deuterium exchange calculation for mass spectra with resolved isotopic fine structure
BMC Bioinformatics volume 15, Article number: 403 (2014)
Hydrogen/deuterium exchange (HDX) coupled to mass spectrometry permits analysis of structure, dynamics, and molecular interactions of proteins. HDX mass spectrometry is confounded by deuterium exchange-associated peaks overlapping with peaks of heavy, natural abundance isotopes, such as carbon-13. Recent studies demonstrated that high-performance mass spectrometers could resolve isotopic fine structure and eliminate this peak overlap, allowing direct detection and quantification of deuterium incorporation.
Here, we present a graphical tool that allows for a rapid and automated estimation of deuterium incorporation from a spectrum with isotopic fine structure. Given a peptide sequence (or elemental formula) and charge state, the mass-to-charge ratios of deuterium-associated peaks of the specified ion is determined. Intensities of peaks in an experimental mass spectrum within bins corresponding to these values are used to determine the distribution of deuterium incorporated. A theoretical spectrum can then be calculated based on the estimated distribution of deuterium exchange to confirm interpretation of the spectrum. Deuterium incorporation can also be detected for ion signals without a priori specification of an elemental formula, permitting detection of exchange in complex samples of unidentified material such as natural organic matter. A tool is also incorporated into QUDeX-MS to help in assigning ion signals from peptides arising from enzymatic digestion of proteins. MATLAB-deployable and standalone versions are available for academic use at qudex-ms.sourceforge.net and agarlabs.com.
Isotopic fine structure HDX-MS offers the potential to increase sequence coverage of proteins being analyzed through mass accuracy and deconvolution of overlapping ion signals. As previously demonstrated, however, the data analysis workflow for HDX-MS data with resolved isotopic fine structure is distinct. QUDeX-MS we hope will aid in the adoption of isotopic fine structure HDX-MS by providing an intuitive workflow and interface for data analysis.
Incubation of molecules in deuterium oxide and measurement of exchange between hydrogen and deuterium (i.e. hydrogen/deuterium exchange, abbreviated as HDX) is a versatile technique for characterizing molecular structure and reactivity. Coupled with mass spectrometry (HDX MS), the development and application of HDX MS to the study of proteins has led to a diverse range of methodologies for the characterization of protein structural dynamics, stability, and molecular interactions. The versatility of HDX MS has been demonstrated through a range of analytical applications including studying folding/unfolding - and oligomerization -; how mutation -, posttranslational modification -, and interaction with other molecules affect protein structures -; and the quality of biopharmaceuticals -. In contrast to other techniques capable of mapping protein dynamics in solution such as nuclear magnetic resonance, HDX MS consumes far less sample, requires lower concentrations, and can analyze larger protein complexes. In a typical HDX MS workflow examining protein structure, proteins are labeled with deuterium oxide around physiological pH for a range of periods of time (from seconds to days), followed by a reaction quenching including a drop to acidic pH (~2.5) and lowered temperature (0°C) . To map deuterium incorporation, and in turn, map extent of solvent exposure of residues, a bottom-up strategy is often adopted whereby proteases (capable of digesting proteins under the acidic quench conditions) break down proteins prior to MS analysis. MS analysis is often performed by electrospray ionization (ESI), although matrix-assisted laser desorption/ionization (MALDI) can also be beneficial ,. Liquid chromatography (LC) is often utilized prior to MS to separate digest peptides and increase protein coverage. A number of computational utilities have been introduced to streamline the analysis of HDX MS data, including HDX-Analyzer , HX-Express , Deuterator , HD Desktop , DEX , Hydra , TOF2H , and Hexicon ,.
One caveat of HDX MS data analysis is that at typical resolving powers (including relatively high resolving powers routinely obtained on FT-Orbitrap instruments), peptide deuteration is detected as a convolution of deuterated peaks with peaks corresponding to the natural abundance distribution of heavy isotopes (i.e. 13C, 15 N, 18O, and 34S/36S). Thus, in order to estimate the distribution of deuterium incorporation, typical HDX analysis determines the mass centroid shift between the deuterated and the native isotope pattern of a peptide to measure average deuterium incorporation. Methods for estimating the full distribution of deuterium exchange, to detect potential bimodal exchange behaviors, have also be devised ,, although these methods are dependent upon sequence information and/or the undeuterated isotope distribution. To minimize back-exchange during chromatography, fast gradients are necessary, increasing the potential for overlapping peak distributions in spectra. High mass resolution can help to alleviate this problem, and algorithms have been developed to estimate accurate deuterium content in spectra with overlapping distributions at low signal-to-noise . At sufficiently high resolving powers (m/∆m 50% ≥ 100,000), elemental compositions of peptides can be used to reliably assign overlapping distributions in a rapid and automated procedure . However, as resolving power increases, the mass defect associated with different atomic nuclei binding energies leads to the detection to the isotopic fine structure. Indeed, at ultrahigh resolving powers (Figure 1A) deuterium exchange-associated peaks become distinct from natural abundance heavy isotopic species, enabling a more direct detection of the distribution of deuterium incorporation .
In the case of Fourier transform ion cyclotron resonance (FTICR) mass spectrometers operating at low-pressure, resolving power m/∆m 50% (where ∆m 50% is the full width at half-maximum peak height) at a particular m/z is limited in magnitude mode by the strength of the magnetic field (B 0) and the time-domain ICR signal acquisition period (T acq) ,:
While long acquisition periods can be used to produce spectra of sufficient quality to resolve deuterium exchange-associated peaks, limiting the acquisition period is necessary in order to accurately detect these peaks on the fast LC-time scale used during HDX experiments. Increasingly powerful magnetic fields provide one means of lowering acquisition time, with fields now up to 21 T , although this involves substantial investment. Solving the phase of the Fourier transform spectrum in absorption mode can further increase resolving power up to 2-fold beyond the magnitude mode limit -. As deuterium-associated peaks at a given charge state are separated by a constant m HDX/z, where m HDX is the mass difference between 2H and 1H minus the mass difference between 13C and 12C, the acquisition period necessary to resolve deuterium-associated peaks of a given absolute mass M can be estimated using:
where is the effective field strength in absorption mode, m adduct is the ionization adduct mass (typically proton mass = ~1.00728 Da) and
is the peak full-width at half height necessary to resolve Lorentzian-shaped peaks of equal height separated by m HDX/z with valley between peaks of 50% maximum height. Thus, in absorption mode at 21 T, <1 sec acquisition periods are required to resolve deuterium exchange-associated peaks on digest peptides less than 1 kDa (Figure 1B). Indeed, for peptides less than 1 kDa, depending on charge state, absorption mode can also yield improvements significant enough to make deuterium-exchange associated peaks resolvable on the more common 9.4 T FTICR with acquisition periods of <1 sec. On high performing Orbitrap instruments, resolving powers in excess of 1 M can be achieved under certain conditions with acquisition times on the time scale of a few seconds , with resolving powers ~240,000 at m/z 400 sufficient to resolve small enough peptides with less than one second acquisition times . Thus, with small enough digest fragments, HDX can be monitored based on resolved isotopic fine structure on a UPLC-time scale. However, in monitoring deuterium incorporation in this fashion, a novel data analysis workflow is necessary to take maximum advantage of this ultrahigh mass resolution data. As yet, no computational utility for this workflow exists.
Here, we present QUDeX-MS (quick ultrahigh-resolution hydrogen/deuterium exchange mass spectrometry), which was developed in order to facilitate the analysis of isotopic fine structure-resolved HDX MS data. The graphical interface, launched either as a stand-alone program or within MATLAB, guides users through the process of interpreting spectra obtained from peptides, proteins, and other molecules, providing the user with the distribution of deuterium incorporated as output.
For a given charge state of an ion, the distance between deuterium exchange-associated peaks in an isotopically resolved distribution is constant (the mass difference between 2H and 1H divided by charge). Thus, starting from the monoisotopic m/z, the m/z’s associated with deuterium-incorporated forms of the natural abundance monoisotopic molecule, the “pseudomonoisotopic peaks” , can be calculated. For example, in the case of human substance P (C-terminally amidated RPKPQQFFGLM), the monoisotopic [M + 2H]2+ is found at m/z 674.3713 and the pseudomonoisotopic m/z’s corresponding to incorporation of 1, 2, or 3 deuterium are m/z 674.8745, m/z 675.3776, and m/z 675.8808, respectively (Figure 2A). In isotopic distributions where fine structure is resolved, these deuterium-associated peaks are distinct from the peaks corresponding to the presence of natural abundance isotopes such as the [M + 2H]2+ of substance P with exactly two 13C at m/z 675.3747, substance P with exactly one 13C and one deuterium at m/z 675.3762, and substance P with exact two deuterium at m/z 675.3776 (Figure 2B-C).
QUDeX-MS determines deuterium incorporation by first finding the percent contribution each pseudomonoisotopic peak has relative to the sum of all the pseudomonoisotopic peaks plus the monoisotopic peak (Figure 3). The average number of deuterium is then determined by multiplying these percentages by the corresponding number of deuterium incorporated for each pseudomonoisotopic m/z. For each of these peaks, the relative distribution of deuterium-associated peak intensities is found, and from that the mean number of deuteriums incorporated is calculated . In addition to estimating the number of deuterium incorporated into the natural abundance monoisotopic species, this process is repeated for each natural abundance isotopomer that consists of all the most abundant naturally occurring stable isotopes (i.e. 14 N, 16O, and 32S) except for an increasing number of 13C (i.e. the isotopomers corresponding to one, two, etc. 13C and no other heavy atoms). A final estimated mean number of deuterium incorporated is calculated based on the most abundant of the 13C-associated isotopomer distributions, which is often the monoisotopic (13C = 0) species for smaller peptides.
In order to calculate the distribution of deuterium incorporated in a given molecule such as a peptide or protein, a spectrum to be analyzed can be loaded in QUDeX-MS. To demonstrate usage of QUDeX-MS, a spectrum of deuterated substance P is provided as sample data with the installation package. Spectra to be analyzed should be formatted as a two-column (m/z and corresponding intensity) ASCII file with no headers and standard column delimiters (i.e. whitespace or commas). Spectra can be denoised using wavelets - and baseline subtracted within QUDeX-MS if necessary using the included Preprocessing Tool. Adjustments can be made to slide the m/z alignment within QUDeX-MS, but if proper assignments are to be made to a reference sequence or formula, it is necessary the spectra be well-calibrated beforehand. If necessary, a table of deuterium-associated m/z’s (based on a specified monoisotopic mass and charge state) can be exported from QUDeX and used to calibrate spectra in another program prior to loading into QUDeX.
After loading an experimental spectrum to analyze, several routes can be chosen for estimating deuterium incorporation of ion signals in the spectrum (Figure 4). If the exact elemental formula (or amino acid sequence) of the ion signal to be analyzed is known, this can be provided to QUDeX-MS. If the elemental composition of an ion signal being analyzed is not known, or if the user does not want to specify it a priori, deuterium incorporation can be estimated by detecting isotopically resolved ion signals in the spectrum within QUDeX-MS, or by simply specifying a monoisotopic mass from which to calculate a deuterium-associated m/z distribution.
The precise location of the deuterium-associated peaks for a given natural abundance isotopomer are then projected on the spectrum. The m/z alignment (mass calibration) of the spectrum can be adjusted, if necessary, to better fit with the calculated deuterium-associated peaks. Once the user is satisfied with the alignment, the distribution of deuterium incorporation can be determined from the spectrum by QUDeX-MS. If the elemental formula has been specified, the distribution of deuterium incorporation can then be used to calculate a theoretical line spectrum, which is projected onto the loaded spectrum. Calculated theoretical spectra based on estimated deuterium incorporation can also be exported as two-column ASCII files for additional analyses.
For automated processing, spectra can be batched loaded using a simple configuration file (sample file and instructions included with QUDeX-MS releases). If the sequences or formulas are known for spectra being loaded in batch, this can also be specified in the configuration file, and spectra will be processed automatically. A digest tool is provided with QUDeX-MS that can be used to determine potential cleavage fragments from a protein sequence (with or without specifying a protease-specific cleavage rule) that can be used to assign ions in a spectrum or set of spectra loaded into QUDeX-MS. Digest results from QUDeX can be exported and formatted (including adding assignments from tandem MS) to assess sequence coverage using an external tool such as MSTools . Likewise, the estimated deuterium incorporation for analyzed ions in any of the proposed workflows can be exported in a report that includes the estimated mean and standard deviation of incorporation for each ion and the root-mean-squared error determined by comparing the experimental mass of detected peaks in the spectrum with their theoretical values.
For LC-MS experiments, QUDeX-MS will only calculate the deuterium uptake within a given spectrum. Given the acquisition times necessary to achieve isotopic fine structure-quality spectra (and the fast peptide elution times necessary in a bottom-up HDX experiment to prevent back-exchange), the number of scans a given peptide will be in is limited. The user can average scans of ions detected in an LC-MS run and load that as a single spectrum to calculate deuterium uptake. Given the considerable file sizes associated with LC-MS runs at this resolving power, we recommend limiting the m/z range of the resulting averaged spectrum to only include a particular ion’s distribution before loading into QUDeX-MS. Alternatively, the user can process individual (i.e. not averaged) scans, including full scans (although we recommend only loading individual isotopic distributions, if possible), across an LC-MS run and interpret scan-to-scan variability in deuterium uptake calculated for a given ion appropriately.
Results and discussion
Calculation of deuterium exchange of a molecule of known composition
An early step in a typical workflow for HDX MS experiments includes performing tandem mass spectrometry (MS/MS) on a (non-deuterated) digested protein sample in order to identify peptides that will be monitored in later HDX experiments. In this case, where sequence and modification state of peptides has been determined, spectra imported into QUDeX-MS can be analyzed by inputting these known characteristics to determine deuterium exchange (Figure 5A). In example spectra provided with QUDeX-MS, substance P spectra have been acquired both with and without deuterium exchange. The default sequence when loading ions by sequence or formula is set to substance P to help when first learning how to use QUDeX-MS. The substance P analyzed for the example spectra was amidated and the appropriate modification (addition of hydrogen and nitrogen atoms and removal of one oxygen) is also provided as default. A specific charge state must also be chosen, which should correspond to a charge state observed in the mass spectrum to be analyzed. By default, QUDeX-MS assumes the user would like to analyze the protonated form of the input molecule, and thus will adjust the elemental formula by adding the appropriate number of protons (depending on charge state). However, this automatic addition of protons can be disabled, allowing the user to add alternative adducts (such as Na+ or K+) directly into the elemental formula. Additional modification to a peptide sequence can be performed this way as well.
Once specification of the desired ion to be analyzed is completed, the distribution of m/z values corresponding to the deuterium-incorporated forms of the natural abundance monoisotopic molecule are calculated and projected onto the spectrum (Figure 5B). Now that the ion to be analyzed is loaded, the deuterium incorporation can be calculated automatically. In the example of substance P, we see the mean number of incorporated deuterium is estimated from the spectrum to be 5.21 based on the distribution of intensities observed from the deuterium peaks associated with the natural abundance monoisotopic substance P. This was consistent with the value calculated previously using manual annotation of the spectrum .
Modeling theoretical spectrum from calculated deuterium exchange
Included within QUDeX-MS is a method for simulating an isotopic fine structure line spectrum based on the estimated deuterium incorporation distribution and the specified elemental formula. Using the percent contributions of the pseudomonoisotopic peaks, a set of isotopic fine structure spectra are simulated  using a set of isotope abundance parameters with a fixed number of deuterium corresponding to the particular pseudomonoisotopic peak. That is, given that the elemental formula for amidated substance P is C63H98N18O13S, the [M + 2H]2+ distribution based on the population of substance P with no deuterium exchange is modeled using the formula C63H100N18O13S and accepted natural isotope abundances. In turn, the [M + 2H]2+ distributions for substance P with exactly 1, 2, and 3 deuterium exchanged are modeled using the formulae C63DH99N18O13S, C63D2H98N18O13S, and C63D3H97N18O13S, and so forth. These distributions are then added together, weighted relative to one another based on their percent contribution estimated from the experimental spectrum. This composite spectrum is projected onto the experimental spectrum (Figure 5C). In the analysis of the example spectrum of deuterated substance P, the three major peaks in fine structure distributions starting from the right moving leftwards are generally the deuterium-associated peaks corresponding to substance P with exactly 0, 1, or 2 13C (and no other natural abundance heavy isotopes).
Calculation of deuterium exchange from spectra without specifying chemical formula a priori
Assuming i) a charge state can be determined from an ion signal and ii) the monoisotopic peak in a particular isotopic distribution can be accurately assigned, the m/z’s of the deuterium-associated pseudomonoisotopic peaks can be calculated as these are offset by a constant distance determined by the charge and the mass of deuterium. Thus, deuterium incorporation can be estimated for many molecules without a priori knowledge of formula or the spectrum of the undeuterated species. The user can directly specify a monoisotopic m/z and charge state to examine, or all isotopically resolved distributions detectable in the spectra can be analyzed (Figure 6). This could be applied to estimating exchange of ions in complex samples of unidentified material such as natural organic matter, as well as studying proteins with unspecified posttranslational modifications or sequence variations. However, it is important to make sure the monoisotopic m/z, even if it is not detected, is properly specified in order to correctly estimate deuterium exchange.
Identification of digest peptides in a deuterated sample from a reference protein sequence
In bottom-up HDX experiments studying protein conformation, a proteolytic enzyme is used to digest the protein after the deuterium exchange reaction has been quenched by low pH. To help facilitate use of QUDeX-MS in the context of a typical HDX MS experiment, we have incorporated a “Digest Tool” within the program that allows the user to specify a protein and digest it in silico. The protein can be digested using a specified enzyme definition. However, HDX experiments require enzymes such as pepsin that are active at low pH but often cleave with broad specificity, leading to a wide variety of peptides with difficult to predict cleavage sites. High mass accuracy can be particularly useful in this case to help in the accurate assignment of peptide fragments ,, although some ambiguities will inevitably remain (such as peptides flanked by the same amino acid on either side). Thus, the QUDeX-MS Digest Tool also allows determining all possible digest fragments of a certain length. Digest peptides can be imported directly into the Command Window, edited further to add modifications, or the digest fragments can be matched to peaks found in one or all loaded spectra (Figure 7). We also specifically incorporate the ability to match digest fragments against a deuterated spectrum in the case where the monoisotopic peak is not detected. Digest sequences can be matched against detected isotopic envelopes in a spectrum (even when they are missing a detectable monoisotopic peak) using a combination of (high mass accuracy) mass defect and information regarding the number of exchange sites. This results in the monoisotopic peak being “assigned” (based on the theoretical mass of the peptide) even though it is not detectable in the spectrum. Included with QUDeX-MS is a simulated set of spectra (with resolving power = 1.5 M) from a pepsin digest of deuterated SOD1 based on results we previously observed  in order to illustrate how to use the Digest Tool within QUDeX-MS to process fine structure HDX MS data.
Directions for future development
Currently, QUDeX-MS is specialized for the estimation of deuterium incorporation in the limit where isotopic fine structure is well-resolved, which can be challenging to achieve experimentally. Without well-resolved fine structure, the ultrahigh mass accuracy of FT-ICR instruments can be useful in increasing the number of unambiguous peptide assignments . Indeed, automated and robust procedures have been developed to leverage the ability to resolve closely overlapping distributions (on the order of a few milliDaltons) that can be resolved on FT-ICR instruments . QUDeX-MS offers the opportunity to analyze deuterium-exchange peaks directly, but requires these peaks be resolved in order to function accurately. While resolution of deuterium-associated fine structure is not yet routine, being able to work interchangeably between the resolving power regimes where these peaks are resolved and when they are not would be useful in order to fully utilize datasets.
Applications already available for analyzing and managing data from HDX MS experiments already have many features that enable the visualization and statistical analysis of HDX MS data. Further incorporation of features from these programs would help further bridging the gap between the already existing capabilities of traditional HDX MS and isotopic fine structure HDX MS. Likewise, high mass accuracy afforded on this instrumentation platform can also aid in the unambiguous assignment of peptides, including in the case where only the deuterium-exchanged sample is analyzed, as described here. Both the high mass accuracy and ability to resolve isotopic composition could also be used, combined with HDX data, for the unambiguous identification of the elemental composition of small molecules, as well as their structural elucidation, particularly in complex mixtures such as in natural organic matter or during metabolite profiling. Combining tools used for elemental composition prediction - with the ability for QUDeX-MS to estimate deuterium incorporation would provide an integrated data analysis environment for data of this type. Indeed, if elemental composition could be determined from accurate mass and/or isotopic composition, small molecule structure prediction combined with HDX data could be used to potentially assign structure as well.
Isotopic fine structure HDX-MS offers the potential to further increase sequence coverage of proteins being analyzed through high mass accuracy and deconvolution of overlapping ion signals. Many recent advancements in mass spectrometers, including the commercialization of the dynamically harmonized Fourier transform ion cyclotron resonance cell by Ivan Boldin and Eugene Nikolaev , have further helped researchers to routinely resolve isotopic fine structure of peptides. Isotopic fine structure HDX-MS has also been demonstrated to be useful in the characterization of molecules in complex mixtures of natural organic matter ,. As previously demonstrated, however, the data analysis workflow for HDX-MS data with resolved isotopic fine structure is distinct. QUDeX-MS we hope will help facilitate the adoption of this technique via streamlining the data analysis process for researchers.
Availability and requirements
Project name: QUDeX-MS
Project home page: qudex-ms.sourceforge.net, www.agarlabs.com
Operating system(s): Precompiled versions available for PC (32- and 64-bit) and OSX
Programming language: MATLAB
Other requirements: None
Any restrictions to use by non-academics: None
Fourier transform ion cyclotron resonance
Yi Q, Baker D: Direct evidence for a two-state protein unfolding transition from hydrogen-deuterium exchange, mass spectrometry, and NMR. Protein Sci. 1996, 5 (6): 1060-1066. 10.1002/pro.5560050608.
Engen JR, Wales TE, Chen SG, Marzluff EM, Hassell KM, Weis DD, Smithgall TE: Partial cooperative unfolding in proteins as observed by hydrogen exchange mass spectrometry. Int Rev Phys Chem. 2013, 32 (1): 96-127. 10.1080/0144235X.2012.751175.
Khanal A, Pan Y, Brown LS, Konermann L: Pulsed hydrogen/deuterium exchange mass spectrometry for time-resolved membrane protein folding studies. J Mass Spectrom. 2012, 47 (12): 1620-1626. 10.1002/jms.3127.
Eghiaian F, Daubenfeld T, Quenet Y, van Audenhaege M, Bouin AP, van der Rest G, Grosclaude J, Rezaei H: Diversity in prion protein oligomerization pathways results from domain expansion as revealed by hydrogen/deuterium exchange and disulfide linkage. Proc Natl Acad Sci U S A. 2007, 104 (18): 7414-7419. 10.1073/pnas.0607745104.
Hodkinson JP, Radford SE, Ashcroft AE: The role of conformational flexibility in beta(2)-microglobulin amyloid fibril formation at neutral pH. Rapid Commun Mass Spectrom. 2012, 26 (16): 1783-1792. 10.1002/rcm.6282.
Sitkiewicz E, Tarnowski K, Poznanski J, Kulma M, Dadlez M: Oligomerization interface of RAGE receptor revealed by MS-monitored hydrogen deuterium exchange. PLoS One 2013, 8(10):e76353.,
Smirnovas V, Kim JI, Lu XJ, Atarashi R, Caughey B, Surewicz WK: Distinct structures of scrapie prion protein (PrPSc)-seeded versus spontaneous recombinant prion protein fibrils revealed by hydrogen/deuterium exchange. J Biol Chem. 2009, 284 (36): 24233-24241. 10.1074/jbc.M109.036558.
Pan J, Han J, Borchers CH, Konermann L: Structure and dynamics of small soluble abeta(1–40) oligomers studied by top-down hydrogen exchange mass spectrometry. Biochemistry. 2012, 51 (17): 3694-3703. 10.1021/bi3002049.
Molnar KS, Karabacak NM, Johnson JL, Wang Q, Tiwari A, Hayward LJ, Coales SJ, Hamuro Y, Agar JN: A common property of amyotrophic lateral sclerosis-associated variants: destabilization of the copper/zinc superoxide dismutase electrostatic loop. J Biol Chem. 2009, 284 (45): 30965-30973. 10.1074/jbc.M109.023945.
Tang L, Coales SJ, Morrow JA, Edmunds T, Hamuro Y: Characterization of the N370S mutant of glucocerebrosidase by hydrogen/deuterium exchange mass spectrometry. Chem BioChem. 2012, 13 (15): 2243-2250.
Rose RJ, van Berkel PHC, van den Bremer ETJ, Labrijn AF, Vink T, Schuurman J, Heck AJR, Parren PWHI: Mutation of Y407 in the CH3 domain dramatically alters glycosylation and structure of human IgG. MAbs. 2013, 5 (2): 219-228. 10.4161/mabs.23532.
Chetty PS, Ohshiro M, Saito H, Dhanasekaran P, Lund-Katz S, Mayne L, Englander W, Phillips MC: Effects of the Iowa and Milano mutations on apolipoprotein a-I structure and dynamics determined by hydrogen exchange and mass spectrometry. Biochemistry. 2012, 51 (44): 8993-9001. 10.1021/bi300926j.
Hamuro Y, Wong L, Shaffer J, Kim JS, Stranz DD, Jennings PA, Woods VL, Adams JA: Phosphorylation driven motions in the COOH-terminal Src kinase, Csk, revealed through enhanced hydrogen-deuterium exchange and mass spectrometry (DXMS). J Mol Biol. 2002, 323 (5): 871-881. 10.1016/S0022-2836(02)01003-3.
Bobst CE, Thomas JJ, Salinas PA, Savickas P, Kaltashov IA: Impact of oxidation on protein therapeutics: Conformational dynamics of intact and oxidized acid-beta-glucocerebrosidase at near-physiological pH. Protein Sci. 2010, 19 (12): 2366-2378. 10.1002/pro.517.
Mitra G, Muralidharan M, Narayanan S, Pinto J, Srinivasan K, Mandal AK: Glutathionylation induced structural changes in oxy human hemoglobin analyzed by backbone amide hydrogen/deuterium exchange and MALDI-mass spectrometry. Bioconjug Chem. 2012, 23 (12): 2344-2353. 10.1021/bc300291u.
Balchin D, Stoychev SH, Dirr HW: S-nitrosation destabilizes glutathione transferase p 1–1. Biochemistry. 2013, 52 (51): 9394-9402. 10.1021/bi401414c.
Pan JX, Borchers CH: Top-down structural analysis of posttranslationally modified proteins by Fourier transform ion cyclotron resonance-MS with hydrogen/deuterium exchange and electron capture dissociation. Proteomics. 2013, 13 (6): 974-981. 10.1002/pmic.201200246.
Chalmers MJ, Busby SA, Pascal BD, West GM, Griffin PR: Differential hydrogen/deuterium exchange mass spectrometry analysis of protein-ligand interactions. Expert Rev Proteomics. 2011, 8 (1): 43-59. 10.1586/epr.10.109.
Lanman J, Prevelige PE: High-sensitivity mass spectrometry for imaging subunit interactions: hydrogen/deuterium exchange. Curr Opin Struct Biol. 2004, 14 (2): 181-188. 10.1016/j.sbi.2004.03.006.
Xiao H, Kaltashov IA, Eyles SJ: Indirect assessment of small hydrophobic ligand binding to a model protein using a combination of ESI MS and HDX/ESI MS. J Am Soc Mass Spectrom. 2003, 14 (5): 506-515. 10.1016/S1044-0305(03)00135-1.
Tropak MB, Kornhaber GJ, Rigat BA, Maegawa GH, Buttner JD, Blanchard JE, Murphy C, Tuske SJ, Coales SJ, Hamuro Y, Brown ED, Mahuran DJ: Identification of pharmacological chaperones for Gaucher disease and characterization of their effects on beta-glucocerebrosidase by hydrogen/deuterium exchange mass spectrometry. Chem BioChem. 2008, 9 (16): 2650-2662.
Chalmers MJ, Busby SA, Pascal BD, He Y, Hendrickson CL, Marshall AG, Griffin PR: Probing protein ligand interactions by automated hydrogen/deuterium exchange mass spectrometry. Anal Chem. 2006, 78 (4): 1005-1014. 10.1021/ac051294f.
Weinreb PH, Li S, Gao SX, Liu T, Pepinsky RB, Caravella JA, Lee JH, Woods VL: Dynamic structural changes are observed upon collagen and metal ion binding to the integrin alpha 1 I domain. J Biol Chem. 2012, 287 (39): 32897-32912. 10.1074/jbc.M112.354365.
Trelle MB, Dupont DM, Madsen JB, Andreasen PA, Jorgensen TJD: Dissecting the effect of RNA aptamer binding on the dynamics of plasminogen activator inhibitor 1 using hydrogen/deuterium exchange mass spectrometry. ACS Chem Biol. 2014, 9 (1): 174-182. 10.1021/cb400619v.
Frego L, Gautschi E, Martin L, Davidson W: The determination of high-affinity protein/inhibitor binding constants by electrospray ionization hydrogen/deuterium exchange mass spectrometry. Rapid Commun Mass Spectrom. 2006, 20 (16): 2478-2482. 10.1002/rcm.2614.
Kaltashov IA, Bobst CE, Abzalimov RR, Berkowitz SA, Houde D: Conformation and dynamics of biopharmaceuticals: transition of mass spectrometry-based tools from academe to industry. J Am Soc Mass Spectrom. 2010, 21 (3): 323-337. 10.1016/j.jasms.2009.10.013.
Kaltashov IA, Bobst CE, Abzalimov RR, Wang G, Baykal B, Wang S: Advances and challenges in analytical characterization of biotechnology products: mass spectrometry-based approaches to study properties and behavior of protein therapeutics. Biotechnol Adv. 2012, 30 (1): 210-222. 10.1016/j.biotechadv.2011.05.006.
Bobst CE, Abzalimov RR, Houde D, Kloczewiak M, Mhatre R, Berkowitz SA, Kaltashov IA: Detection and characterization of altered conformations of protein pharmaceuticals using complementary mass spectrometry-based approaches. Anal Chem. 2008, 80 (19): 7473-7481. 10.1021/ac801214x.
Engen JR: Analysis of protein conformation and dynamics by hydrogen/deuterium exchange MS. Anal Chem. 2009, 81 (19): 7870-7875. 10.1021/ac901154s.
Mandell JG, Falick AM, Komives EA: Measurement of amide hydrogen exchange by MALDI-TOF mass spectrometry. Anal Chem. 1998, 70 (19): 3987-3995. 10.1021/ac980553g.
Nazabal A, Laguerre M, Schmitter JM: Hydrogen/deuterium exchange on yeast ATPase supramolecular protein complex analyzed at high sensitivity by MALDI mass spectrometry. J Am Soc Mass Spectrom. 2003, 14 (5): 471-481. 10.1016/S1044-0305(03)00124-7.
Liu S, Liu L, Uzuner U, Zhou X, Gu M, Shi W, Zhang Y, Dai SY, Yuan JS: HDX-analyzer: a novel package for statistical analysis of protein structure dynamics. BMC Bioinformatics 2011, 12(Suppl 1):S43.,
Weis DD, Engen JR, Kass IJ: Semi-automated data processing of hydrogen exchange mass spectra using HX-Express. J Am Soc Mass Spectrom. 2006, 17 (12): 1700-1703. 10.1016/j.jasms.2006.07.025.
Pascal BD, Chalmers MJ, Busby SA, Mader CC, Southern MR, Tsinoremas NF, Griffin PR: The Deuterator: software for the determination of backbone amide deuterium levels from H/D exchange MS data. BMC Bioinformatics 2007, 8:156.,
Pascal BD, Chalmers MJ, Busby SA, Griffin PR: HD desktop: an integrated platform for the analysis and visualization of H/D exchange data. J Am Soc Mass Spectrom. 2009, 20 (4): 601-610. 10.1016/j.jasms.2008.11.019.
Hotchko M, Anand GS, Komives EA, Ten Eyck LF: Automated extraction of backbone deuteration levels from amide H/2H mass spectrometry experiments. Protein Sci. 2006, 15 (3): 583-601. 10.1110/ps.051774906.
Slysz GW, Baker CA, Bozsa BM, Dang A, Percy AJ, Bennett M, Schriemer DC: Hydra: software for tailored processing of H/D exchange data from MS or tandem MS analyses. BMC Bioinformatics 2009, 10:162.,
Nikamanon P, Pun E, Chou W, Koter MD, Gershon PD: “TOF2H”: a precision toolbox for rapid, high density/high coverage hydrogen-deuterium exchange mass spectrometry via an LC-MALDI approach, covering the data pipeline from spectral acquisition to HDX rate analysis. BMC Bioinformatics 2008, 9:387.,
Kreshuk A, Stankiewicz M, Lou XH, Kirchner M, Hamprecht FA, Mayer MP: Automated detection and analysis of bimodal isotope peak distributions in H/D exchange mass spectrometry using HeXicon. Int J Mass Spectrom. 2011, 302 (1–3): 125-131. 10.1016/j.ijms.2010.08.025.
Lindner R, Lou X, Reinstein J, Shoeman RL, Hamprecht FA, Winkler A: Hexicon 2: automated processing of hydrogen-deuterium exchange mass spectrometry data with improved deuteration distribution estimation. J Am Soc Mass Spectrom. 2014, 25 (6): 1018-1028. 10.1007/s13361-014-0850-y.
Guttman M, Weis DD, Engen JR, Lee KK: Analysis of overlapped and noisy hydrogen/deuterium exchange mass spectra. J Am Soc Mass Spectrom. 2013, 24 (12): 1906-1912. 10.1007/s13361-013-0727-5.
Kazazic S, Zhang HM, Schaub TM, Emmett MR, Hendrickson CL, Blakney GT, Marshall AG: Automated data reduction for hydrogen/deuterium exchange experiments, enabled by high-resolution Fourier transform ion cyclotron resonance mass spectrometry. J Am Soc Mass Spectrom. 2010, 21 (4): 550-558. 10.1016/j.jasms.2009.12.016.
Liu Q, Easterling ML, Agar JN: Resolving isotopic fine structure to detect and quantify natural abundance- and hydrogen/deuterium exchange-derived isotopomers. Anal Chem. 2014, 86 (1): 820-825. 10.1021/ac403365g.
Marshall AG, Comisarow MB, Parisod G: Relaxation and spectral-line shape in Fourier-transform ion resonance spectroscopy. J Chem Phys. 1979, 71 (11): 4434-4444. 10.1063/1.438196.
Marshall AG, Hendrickson CL, Jackson GS: Fourier transform ion cyclotron resonance mass spectrometry: A primer. Mass Spectrom Rev. 1998, 17 (1): 1-35. 10.1002/(SICI)1098-2787(1998)17:1<1::AID-MAS1>3.0.CO;2-K.
Hsu CS, Hendrickson CL, Rodgers RP, McKenna AM, Marshall AG: Petroleomics: advanced molecular probe for petroleum heavy ends. J Mass Spectrom. 2011, 46 (4): 337-343. 10.1002/jms.1893.
Xian F, Hendrickson CL, Blakney GT, Beu SC, Marshall AG: Automated broadband phase correction of Fourier transform ion cyclotron resonance mass spectra. Anal Chem. 2010, 82 (21): 8807-8812. 10.1021/ac101091w.
Qi YL, Thompson CJ, Van Orden SL, O’Connor PB: Phase correction of Fourier transform ion cyclotron resonance mass spectra using MatLab. J Am Soc Mass Spectrom. 2011, 22 (1): 138-147. 10.1007/s13361-010-0006-7.
Qi YL, Barrow MP, Li HL, Meier JE, Van Orden SL, Thompson CJ, O’Connor PB: Absorption-mode: the next generation of Fourier transform mass spectra. Anal Chem. 2012, 84 (6): 2923-2929. 10.1021/ac3000122.
Kilgour DPA, Wills R, Qi YL, O’Connor PB: Autophaser: an algorithm for automated generation of absorption mode spectra for FT-ICR MS. Anal Chem. 2013, 85 (8): 3903-3911. 10.1021/ac303289c.
Denisov E, Damoc E, Lange O, Makarov A: Orbitrap mass spectrometry with resolving powers above 1,000,000. Int J Mass Spectrom. 2012, 325: 80-85. 10.1016/j.ijms.2012.06.009.
Michalski A, Damoc E, Lange O, Denisov E, Nolting D, Muller M, Viner R, Schwartz J, Remes P, Belford M, Dunyach JJ, Cox J, Horning S, Mann M, Makarov A: Ultra high resolution linear ion trap Orbitrap mass spectrometer (Orbitrap Elite) facilitates top down LC MS/MS and versatile peptide fragmentation modes. Mol Cell Proteomics 2012, 11(3):O111.01369.,
Barclay VJ, Bonner RF, Hamilton IP: Application of wavelet transforms to experimental spectra: smoothing, denoising, and data set compression. Anal Chem. 1997, 69 (1): 78-90. 10.1021/ac960638m.
Donoho DL, Johnstone IM: Adapting to unknown smoothness via wavelet shrinkage. J Am Stat Assoc. 1995, 90 (432): 1200-1224. 10.1080/01621459.1995.10476626.
Coombes KR, Tsavachidis S, Morris JS, Baggerly KA, Hung MC, Kuerer HM: Improved peak detection and quantification of mass spectrometry data acquired from surface-enhanced laser desorption and ionization by denoising spectra with the undecimated discrete wavelet transform. Proteomics. 2005, 5 (16): 4107-4117. 10.1002/pmic.200401261.
Kavan D, Man P: MSTools – Web based application for visualization and presentation of HXMS data. Int J Mass Spectrom. 2011, 302 (1–3): 53-58. 10.1016/j.ijms.2010.07.030.
Rockwood AL, VanOrden SL, Smith RD: Ultrahigh resolution isotope distribution calculations. Rapid Commun Mass Spectrom. 1996, 10 (1): 54-59. 10.1002/(SICI)1097-0231(19960115)10:1<54::AID-RCM444>3.0.CO;2-Z.
Lam TT, Lanman JK, Emmett MR, Hendrickson CL, Marshall AG, Prevelige PE: Mapping of protein:protein contact surfaces by hydrogen/deuterium exchange, followed by on-line high-performance liquid chromatography-electrospray ionization Fourier-transform ion-cyclotron-resonance mass analysis. J Chromatogr A. 2002, 982 (1): 85-95. 10.1016/S0021-9673(02)01357-2.
Zhang HM, Kazazic S, Schaub TM, Tipton JD, Emmett MR, Marshall AG: Enhanced digestion efficiency, peptide ionization efficiency, and sequence resolution for protein hydrogen/deuterium exchange monitored by Fourier transform ion cyclotron resonance mass spectrometry. Anal Chem. 2008, 80 (23): 9034-9041. 10.1021/ac801417d.
Okawa S, Fischer B, Krijgsveld J: Properties of isotope patterns and their utility for peptide identification in large-scale proteomic experiments. Rapid Commun Mass Spectrom. 2013, 27 (9): 1067-1075. 10.1002/rcm.6551.
Pluskal T, Uehara T, Yanagida M: Highly accurate chemical formula prediction tool utilizing high-resolution mass spectra, MS/MS fragmentation, heuristic rules, and isotope pattern matching. Anal Chem. 2012, 84 (10): 4396-4403. 10.1021/ac3000418.
Patiny L, Borel A: ChemCalc: a building block for tomorrow’s chemical infrastructure. J Chem Inf Model. 2013, 53 (5): 1223-1228. 10.1021/ci300563h.
Kind T, Fiehn O: Seven golden rules for heuristic filtering of molecular formulas obtained by accurate mass spectrometry. BMC Bioinformatics 2007, 8:105.,
Kind T, Fiehn O: Metabolomic database annotations via query of elemental compositions: mass accuracy is insufficient even at less than 1 ppm. BMC Bioinformatics 2006, 7:234.,
Boldin IA, Nikolaev EN: Fourier transform ion cyclotron resonance cell with dynamic harmonization of the electric field in the whole volume by shaping of the excitation and detection electrode assembly. Rapid Commun Mass Spectrom. 2011, 25 (1): 122-126. 10.1002/rcm.4838.
Cho Y, Ahmed A, Kim S: Application of atmospheric pressure photo ionization hydrogen/deuterium exchange high-resolution mass spectrometry for the molecular level speciation of nitrogen compounds in heavy crude oils. Anal Chem. 2013, 85 (20): 9758-9763. 10.1021/ac402157r.
Kostyukevich Y, Kononikhin A, Popov I, Kharybin O, Perminova I, Konstantinov A, Nikolaev E: Enumeration of labile hydrogens in natural organic matter by use of hydrogen/deuterium exchange Fourier transform ion cyclotron resonance mass spectrometry. Anal Chem. 2013, 85 (22): 11007-11013. 10.1021/ac402609x.
The authors would like to express great appreciation for the insightful comments provided by the reviewers, which helped expand the functionality of QUDeX-MS immensely. The authors would like to acknowledge Jeniffer Quijada for initial testing of QUDeX-MS and her suggestions regarding the QUDeX-MS interface. Research reported in this publication was supported by the National Institute of Neurological Disorders and Stroke of the National Institutes of Health under award number R01NS065263.
The authors declare that they have no competing interests.
QL obtained sample spectra and performed initial analysis. JPS implemented the complete data analysis workflow and interface. QL, JPS, and JNA developed the data analysis framework. JPS and JNA wrote the manuscript. All authors have given approval to the final version of the manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Salisbury, J.P., Liu, Q. & Agar, J.N. QUDeX-MS: hydrogen/deuterium exchange calculation for mass spectra with resolved isotopic fine structure. BMC Bioinformatics 15, 403 (2014). https://doi.org/10.1186/s12859-014-0403-1
- Hydrogen/deuterium exchange
- Isotopic fine structure
- Fourier transform ion cyclotron resonance
- Mass spectrometry
- Protein structure
- Protein analysis