Applying a new algorithm for obtaining site specific rate constants for H/D exchange of the gas phase proton-bound arginine dimer

Resolution of protein hydrogen/deuterium exchange by fitting amide exchange probabilities to the peptide isotopic envelopes

RATIONALE

Mass spectra processing in protein hydrogen/deuterium (H/D) exchange has been remarkably improved by the introduction of fitting of the amide exchange probabilities to peptide isotopic envelope intensities (Kan et al., 2013), in contrast to methods in which only the peptide deuterium uptakes (centroid shifts of isotopic envelopes) are used. However, the known implementations are based on the general fitting routines that use only the objective function values. Besides, applicability of more than one fitting method makes necessary their comparative evaluation.

METHODS

Two fitting methods were considered: the common least squares and the fitting of the multinomial distribution representing the number of deuterium atoms exchanged in the individual peptides. Both methods were applied either directly to the isotopic envelope data or to the deuterium distributions obtained by envelope deconvolution (i.e. de-isotoping).

RESULTS

An autonomous Matlab script was prepared, based on the exact expressions for the gradient and Hessian of the objective function, with the trust-region algorithm implemented in the compact analytical form recently made available. The least-squares fitting to the envelope data produced the best results, with the greatest precision and good coverage of exact values by the confidence intervals. The deuterium distributions were sensitive to the (simulated) experimental error whose progression by envelope deconvolution caused degradation in accuracy. The multinomial distribution fitting exhibited poor performance due to inadequate representation of the experimental error and missing of the appropriate weight parameters. Some specific peptide arrangement details were discussed as potential sources of ambiguity in the fitting results.

CONCLUSIONS

The method of fitting to peptide isotopic envelopes has been improved by using the exact gradient and Hessian of the objective function. The fitting should be repeated with different initial guesses in order to find not only the global minimum, but also the local minima with similar depths which may exist due to eventual ambiguity of the fitting results.

Computational methods and challenges in hydrogen/deuterium exchange mass spectrometry

Hydrogen/Deuterium exchange (HDX) has been applied, since the 1930s, as an analytical tool to study the structure and dynamics of (small) biomolecules. The popularity of using HDX to study proteins increased drastically in the last two decades due to the successful combination with mass spectrometry (MS). Together with this growth in popularity, several technological advances have been made, such as improved quenching and fragmentation. As a consequence of these experimental improvements and the increased use of protein-HDXMS, large amounts of complex data are generated, which require appropriate analysis. Computational analysis of HDXMS requires several steps. A typical workflow for proteins consists of identification of (non-)deuterated peptides or fragments of the protein under study (local analysis), or identification of the deuterated protein as a whole (global analysis); determination of the deuteration level; estimation of the protection extent or exchange rates of the labile backbone amide hydrogen atoms; and a statistically sound interpretation of the estimated protection extent or exchange rates. Several algorithms, specifically designed for HDX analysis, have been proposed. They range from procedures that focus on one specific step in the analysis of HDX data to complete HDX workflow analysis tools. In this review, we provide an overview of the computational methods and discuss outstanding challenges. © 2016 Wiley Periodicals, Inc. Mass Spec Rev 36:649-667, 2017.

Time-resolved pulsed hydrogen/deuterium exchange mass spectrometry probes gaseous proteins structural kinetics

A pulsed hydrogen/deuterium exchange (HDX) method has been developed for rapid monitoring of the exchange kinetics of protein ions with D2O a few milliseconds after electrospray ionization (ESI). The stepwise gradual evolution of HDX of multiply charged protein ions was monitored using the pulsed HDX mass spectrometry technique. Upon introducing a very short pulse of D2O (in the μs to ms time scale) into the linear ion trap (LIT) of a time-of-flight (TOF) mass spectrometer, bimodal distributions were detected for the ions of cytochrome c and ubiquitin. Mechanistic details of HDX reactions for ubiquitin and cytochrome c in the gas phase were uncovered and the structural transitions were followed by analyzing the kinetics of HDX.

Gas-phase hydrogen/deuterium exchange in a traveling wave ion guide for the examination of protein conformations

Accumulating evidence suggests that solution-phase conformations of small globular proteins and large molecular protein assemblies can be preserved for milliseconds after electrospray ionization. Thus, the study of proteins in the gas phase on this time scale is highly desirable. Here we demonstrate that a traveling wave ion guide (TWIG) of a Synapt mass spectrometer offers a highly suitable environment for rapid and efficient gas-phase hydrogen/deuterium exchange (HDX). Gaseous ND(3) was introduced into either the source TWIG or the TWIG located just after the ion mobility cell, such that ions underwent HDX as they passed through the ND(3) on the way to the time-of-flight analyzer. The extent of deuterium labeling could be controlled by varying the quantity of ND(3) or the speed of the traveling wave. The gas-phase HDX of model peptides corresponded to labeling of primarily fast exchanging sites due to the short labeling times (ranging from 0.1 to 10 ms). In addition to peptides, gas-phase HDX of ubiquitin, cytochrome c, lysozyme, and apomyoglobin were examined. We conclude that HDX of protein ions in a TWIG is highly sensitive to protein conformation, enables the detection of conformers present on submilliseconds time scales, and can readily be combined with ion mobility spectrometry.

Fragmentation of doubly-protonated peptide ion populations labeled by H/D exchange with CD(3)OD

Doubly-protonated bradykinin (RPPGFSPFR) and an angiotensin III analogue (RVYIFPF) were subjected to hydrogen/deuterium (H/D) exchange with CD(3)OD in a Fourier transform ion cyclotron resonance (FT-ICR) mass spectrometer. A bimodal distribution of deuterium incorporation was present for bradykinin after H/D exchange for 90 s at a CD(3)OD pressure of 4 x 10(-7) Torr, indicating the existence of at least two distinct populations. Bradykinin ion populations corresponding to 0-2 and 5-11 deuteriums (i.e., D(0), D(1), D(2), D(5), D(6), D(7), D(8), D(9), D(10), and D(11)) were each monoisotopically selected and fragmented via sustained off-resonance irradiation (SORI) collision-induced dissociation (CID). The D(0)-D(2) ion populations, which correspond to the slower exchanging population, consistently require lower SORI amplitude to achieve a similar precursor ion survival yield as the faster-reacting (D(5)-D(11)) populations. These results demonstrate that conformation/protonation motif has an effect on fragmentation efficiency for bradykinin. Also, the partitioning of the deuterium atoms into fragment ions suggests that the C-terminal arginine residue exchanges more rapidly than the N-terminal arginine. Total deuterium incorporation in the b(1)/y(8) and b(2)/y(7) ion pairs matches very closely the theoretical values for all ion populations studied, indicating that the ions of a complementary pair are likely formed during the same fragmentation event, or that no scrambling occurs upon SORI. Deuterium incorporation into the y(1)/a(8) pseudo-ion pair does not closely match the expected theoretical values. The other peptide, doubly-protonated RVYIFPF, has a trimodal distribution of deuterium incorporation upon H/D exchange with CD(3)OD at a pressure of 1 x 10(-7) Torr for 600 s, indicating at least three distinct ion populations. After 90 s of H/D exchange where at least two distinct populations are detected, the D(0)-D(7) ion populations were monoisotopically selected and fragmented via SORI-CID over a range of SORI amplitudes. The precursor ion survival yield as a function of SORI amplitude falls into two distinct behaviors corresponding to slower- and faster-reacting ion populations. The slower-reacting population requires larger SORI amplitudes to achieve the same precursor ion survival yield as the faster exchanging population. Total deuterium incorporation into the y(2)/b(5) ion pairs matches closely the theoretical values over all ion populations and SORI amplitudes studied. This result indicates the y(2) and b(5) ions are likely formed by the same mechanism over the SORI amplitudes studied.

The gas-phase H/D exchange mechanism of protonated amino acids

A mass spectrometry and Density Functional Theory study of gas-phase H/D exchange in protonated Ala, Cys, Ile, Leu, Met, and Val is reported. Site-specific rate constants were determined and results identify the alpha-amino group as the protonation site. Lack of exchange on the Cys thiol group is explained by the absence of strong intramolecular hydrogen bonding within the reaction complex. In aliphatic amino acids the presence of a methyl group at the beta-C atom was found to lower the site-specific H/D exchange rate for amino hydrogens. Study of the exchange mechanism showed that isotopic exchange occurs in two independent reactions: in one, only the carboxylic hydrogen is exchanged and in the other, both carboxylic and amino group hydrogens exchange. The proposed reaction mechanisms, calculated structures of various species, and a number of structural findings are consistent with experimental data.

A Fast Flow Tube Study of Gas Phase H/D Exchange of Multiply Protonated Ubiquitin

An electrospray ionization (ESI)/fast-flow technique has been applied to the study of gas phase hydrogen/deuterium (H/D) exchange kinetics. Multiply charged ubiquitin ions [ubiquitin + nH](n)(+), in charge states n = 7-13, were reacted with ND(3). The behavior of ND(3) as exchange reagent is different from that of the previously studied reagents, D(2)O and CH(3)OD. Contrary to those, the maximum number of exchanged hydrogen atoms and the overall exchange rate were observed to increase with increasing charge state of the ubiquitin ions. The results are reagent-dependent because the exchange mechanisms are different for the different reagents. This observation is in agreement with a recent conclusion by Beauchamp and co-workers that contrary to the assumption often expressed in earlier studies, H/D exchange kinetics may not directly reflect ion structures. The results for all three reagents are, however, consistent with observations of previous ion mobility experiments that with increasing charge state the conformers change from more compact, partially folded structures to elongated nearly linear ones. H/D exchange of (ubiquitin + 13H)(13+) with ND(3) leads to two separated ion populations reflecting the possible existence of two conformers with different exchange rates. The ions (ubiquitin + 8H)(8+) and (ubiquitin + 11H)(11+) represent a partially folded structure and an unfolded structure, respectively, and were studied in greater detail. The relative abundances of ions were measured in steps of 0.5 m/z (mass-to-charge ratio), as a function of the ND(3) flow rate. The experimental results were simulated by computer fitted curves based on a recently developed algorithm. The algorithm allows the extraction of sets of grouped rate constants. Eight rate constant groups were deduced for each of the two ions. These rate constants correspond to 32 and 44 H/D exchanges for the 8+ and 11+ charged ions, respectively. The results indicate higher individual rates for most of the exchanged atoms in the 11+ ion compared to the 8+ ion.

Protonated Serine Octamer Cluster: Structure Elucidation by Gas-Phase H/D Exchange Reactions

The H/D exchange kinetics of the protonated serine octamer was investigated by both flow-tube and FT-ICR experiments. The exchange was observed to be bimodal in agreement with previous observations. Quantitative analysis of the experimental results led to site-specific H/D exchange rate constants on the basis of which the structures of both ion populations were deduced. We observe the two separate conformers exchanging 33 hydrogens each-in an independent manner and at different rates. This result was achieved through a probabilistic algorithm that groups together equivalent hydrogen atoms having equal rate constants. The slower exchanging population A is assigned an all-zwitterionic structure. Its faster exchanging counterpart B is assigned an all-neutral structure. Population A was found to be more stable toward collision-induced activation than population B. All of these findings are consistent with previous experimental results, thus comprising a self-consistent picture of the protonated serine octamer and its gas-phase properties.

Direct quantitative determination of surface Brønsted acidity of solids by H/D exchange using D2O

A highly flexible method for direct and quantitative determination of surface Brønsted acidity of solids in terms of number, type and strength of the acid sites based on quantitative H/D exchange kinetics between the acid solid and gaseous D2O has been developed and applied to materials covering the whole range of acidity.