Thermodynamics and mechanisms of protonated diglycine decomposition: a guided ion beam study

Perspective: intrinsic interactions of metal ions with biological molecules as studied by threshold collision-induced dissociation and infrared multiple photon dissociation

In this perspective, gas-phase studies of group 1 monocations and group 12 dications with amino acids and small peptides are highlighted. Although the focus is on two experimental techniques, threshold collision-induced dissociation and infrared multiple photon dissociation action spectroscopy, these methods as well as complementary approaches are summarized. The synergistic interplay with theory, made particularly powerful by the small sizes of the systems explored and the absence of solvent and support, is also elucidated. Importantly, these gas-phase methods permit quantitative insight into the structures and thermodynamics of metal cations interacting with biological molecules. Periodic trends in how these interactions vary as the metal cations get heavier are discussed as are quantitative trends with changes in the amino acid side chain and effects of hydration. Such trends allow these results to transcend the limitations associated with the biomimetic model systems.

Energetics and mechanisms for decomposition of cationized amino acids and peptides explored using guided ion beam tandem mass spectrometry

Fragmentation studies of cationized amino acids and small peptides as studied using guided ion beam tandem mass spectrometry (GIBMS) are reviewed. After a brief examination of the key attributes of the GIBMS approach, results for a variety of systems are examined, compared, and contrasted. Cationization of amino acids, diglycine, and triglycine with alkali cations generally leads to dissociations in which the intact biomolecule is lost. Exceptions include most lithiated species as well as a few examples for sodiated and one example for potassiated species. Like the lithiated species, cationization by protons leads to numerous dissociation channels. Results for protonated glycine, cysteine, asparagine, diglycine, and a series of tripeptides are reviewed, along with the thermodynamic consequences that can be gleaned. Finally, the important physiological process of the deamidation of asparagine (Asn) residues is explored by the comparison of five dipeptides in which the C-terminal partner (AsnXxx) is altered. The GIBMS thermochemistry is shown to correlate well with kinetic results from solution phase studies.

Potassium Binding Interactions with Aliphatic Amino Acids: Thermodynamic and Entropic Effects Analyzed via a Guided Ion Beam and Computational Study

Noncovalent interactions between alkali metals and amino acids are critical for many biological processes, especially for proper function of protein ion channels; however, many precise binding affinities between alkali metals and amino acids still need to be measured. This study addresses this need by using threshold collision-induced dissociation with a guided ion beam tandem mass spectrometer to measure binding affinities between potassium cations and the aliphatic amino acids: Gly, Ala, hAla, Val, Leu, and Ile. These measurements are supplemented by theoretical calculations and include commentary on effects of enthalpy, entropy, and structural preference. Notably, all levels of theory indicate that the lowest-lying isomers at 298 K have K⁺ binding to the carbonyl oxygen in either a monodentate ([CO]) or bidentate ([CO,OH]) fashion, isomers that are linked in a double-well potential. This complicates the analysis of the data, although does not greatly influence the final results. Analysis of the resulting cross sections includes accounting for multiple ion-molecule collisions, internal energy of reactant ions, and unimolecular decay rates. The resulting experimental bond dissociation energies generally increase as the polarizability of the amino acid increases, results that agree well with quantum chemical calculations done at the B3LYP, B3P86, and MP2(full) levels of theory, with B3LYP-GD3BJ predicting systematically larger values.

Thermodynamics and Reaction Mechanisms for Decomposition of a Simple Protonated Tripeptide, H+GGA: From H+GGG to H+GAG to H+GGA

We present a thorough characterization of fragmentations observed in threshold collision-induced dissociation (TCID) experiments of protonated glycylglycylalanine (H⁺GGA) with Xe using a guided ion beam tandem mass spectrometer. Kinetic energy dependent cross sections for nine ionic products were obtained and analyzed to provide 0 K barriers for the five primary products, [b₂]⁺, [y₁ + 2H]⁺, [b₃]⁺, [y₂ + 2H]⁺, and [a₁]⁺; and four secondary products, [a₂]⁺, [a₃]⁺, high-energy [y₁ + 2H]⁺, and CH₃CHNH₂⁺, after accounting for multiple ion-molecule collisions, the internal energy of reactant ions, unimolecular decay rates, competition between channels, and sequential dissociations. Relaxed potential energy surface scans performed at the B3LYP-GD3BJ/6-311+G(d,p) level of theory are used to identify transition states (TSs) and intermediates of the five primary and three secondary products (with the mechanism of the other secondary product previously established). Geometry optimizations and single point energy calculations of reactants, products, intermediates, and TSs were performed at several levels of theory. These theoretical energies are compared with experimental threshold energies and found to give reasonable agreement, with B3LYP-GD3BJ and M06-2X levels of theory performing slightly better than MP2 and better than B3LYP. The results obtained here are compared with previous results for decomposition of H⁺GGG and H⁺GAG to probe the effect of changing the amino acid sequence. Methylation in H⁺GGA has a significant effect on the competition between the primary sequence products, [b₂]⁺ and [y₁ + 2H]⁺, suppressing the [b₂]⁺ cross section by raising its threshold energy, while enhancing that of [y₁ + 2H]⁺ by lowering its threshold energy.

Backbone Cleavages of Protonated Peptoids upon Collision-Induced Dissociation: Competitive and Consecutive B-Y and A1-YX Reactions

Mass spectrometric techniques and more particularly collision-induced dissociation (CID) experiments represent a powerful method for the determination of the primary sequence of (bio)molecules. However, the knowledge of the ion fragmentation patterns say the dissociation reaction mechanisms is a prerequisite to reconstitute the sequence based on fragment ions. Previous papers proposed that protonated peptoids dissociate following an oxazolone-ring mechanism starting from the O-protonation species and leading to high mass Y sequence ions. Here we revisit this backbone cleavage mechanism by performing CID and ion mobility experiments, together with computational chemistry, on tailor-made peptoids. We demonstrated that the B/Y cleavages of collisionally activated O-protonated peptoids must involve the amide nitrogen protonated structures as the dissociating species, mimicking the CID behavior of protonated peptides. Upon the nucleophilic attack of the oxygen atom of the N-terminal adjacent carbonyl group on the carbonyl carbon atom of the protonated amide, the peptoid ions directly dissociate to form an ion-neutral complex associating an oxazolone ion to the neutral truncated peptoid residue. Dissociation of the ion/neutral complex predominantly produces Y ions due to the high proton affinity of the secondary amide function characteristic of truncated peptoids. Whereas the production of Y_x ions from acetylated peptoids also involves the B/Y pathway, the observation of abundant Y_x ions from non-acetylated peptoid ions is shown in the present study to arise from an A₁-Y_x mechanism. The consecutive and competitive characters of the A₁-Y_x and the B/Y mechanisms are also investigated by drift time-aligned CID experiments.

Thermodynamics and Reaction Mechanisms for Decomposition of a Simple Protonated Tripeptide, H+GAG: a Guided Ion Beam and Computational Study

We present a thorough characterization of fragmentations observed in threshold collision-induced dissociation (TCID) experiments of protonated glycylalanylglycine (H⁺GAG) with Xe using a guided ion beam tandem mass spectrometer. Kinetic energy dependent cross sections for nine ionic products were observed and analyzed to provide 0 K barriers for the six primary products: [b₂]⁺, [y₁ + 2H]⁺, [b₃]⁺, CO loss, [y₂ + 2H]⁺, and [a₁]⁺; and three secondary products: [a₂]⁺, [a₃]⁺, and CH₃CHNH₂⁺, after accounting for multiple ion-molecule collisions, internal energy of reactant ions, unimolecular decay rates, competition between channels, and sequential dissociations. Relaxed potential energy surface scans performed at the B3LYP-GD3BJ/6-311+G(d,p) level of theory are used to identify transition states (TSs) and intermediates of the six primary and one secondary products (where the other two secondary products have mechanisms previously established). Geometry optimizations and single-point energy calculations were performed at several levels of theory. These theoretical energies are compared with experimental threshold energies and are found to give reasonably good agreement, with B3LYP-GD3BJ and M06-2X levels of theory performing better than other levels. The results obtained here are also compared with previous results for decomposition of H⁺GGG. The primary difference observed is a lowering of the threshold for the [b₂]⁺ product ion and a concomitant suppression of the directly competing [y₁ + 2H]⁺ product, the result of specific methylation of the [b₂]⁺ product ion.

Thermodynamics and Reaction Mechanisms of Decomposition of the Simplest Protonated Tripeptide, Triglycine: A Guided Ion Beam and Computational Study

We present a thorough characterization of fragmentations observed in threshold collision-induced dissociation (TCID) experiments of protonated triglycine (H⁺GGG) with Xe using a guided ion beam tandem mass spectrometer (GIBMS). Kinetic energy-dependent cross-sections for 10 ionic products are observed and analyzed to provide 0 K barriers for six primary products: [b₂]⁺, [y₁ + 2H]⁺, [b₃]⁺, CO loss, [y₂ + 2H]⁺, and [a₁]⁺; three secondary products: [a₂]⁺, [a₃]⁺, and [y₂ + 2H - CO]⁺; and two tertiary products: high energy [y₁ + 2H]⁺ and [a₂ - CO]⁺ after accounting for multiple ion-molecule collisions, internal energy of reactant ions, unimolecular decay rates, competition between channels, and sequential dissociations. Relaxed potential energy surface scans performed at the B3LYP-D3/6-311+G(d,p) level of theory are used to identify transition states (TSs) and intermediates of the six primary and one secondary products. Geometry optimizations and single point energy calculations were performed at several levels of theory. These theoretical energies are compared with experimental energies and are found to give reasonably good agreement, in particular for the M06-2X level of theory. This good agreement between experiment and theory validates the reaction mechanisms explored computationally here and elsewhere and allows identification of the product structures formed at threshold energies. The present work presents the first measurement of absolute experimental threshold energies of important sequence ions and non-sequence ions: [y₁ + 2H]⁺, [b₃]⁺, CO loss, [a₁]⁺, and [a₃]⁺, and refines those for [b₂]⁺ and [y₂ + 2H]⁺ previously measured. Graphical Abstract ᅟ.

Thermodynamics and Mechanisms of Protonated Asparaginyl-Glycine Decomposition

Deamidation at asparagine residues, a spontaneous post-translational modification in proteins, plays a significant role in various biological processes and degenerative diseases. In the current work, we present a full description of the deamidation process as well as other key fragmentations (dehydration, peptide bond cleavage, and loss of 2 NH3) from protonated asparaginyl-glycine, H(+)(AsnGly), by studying its kinetic energy dependent collision-induced dissociation (CID) with Xe using a guided ion beam tandem mass spectrometer. These results are compared with those for sustained off-resonance irradiation (SORI)-CID of H(+)(AsnGly) with Ar in a Fourier transform ion cyclotron resonance mass spectrometer. Computationally, simulating annealing methodology and a series of relaxed potential energy scans at the B3LYP/6-31G(d) level were performed to identify all intermediate and transition state (TS) structures for each key reaction. All species were further optimized at the B3LYP and B3LYP-GD3BJ/6-311+G(d,p) levels of theory. Single point energies of all major reaction species were calculated at the B3LYP, B3P86, MP2(full), and B3LYP-GD3BJ levels of theory and using M06-2X for rate-limiting species. Relative energies of intermediates, TSs, and products allow characterization of the elementary and rate limiting steps in H(+)(AsnGly) decomposition. By combining experimental and computational results, the complete mechanistic nature of H(+)(AsnGly) deamidation and other fragmentations is explored and compared to the previously studied H(+)(Asn) complex. The influence of water solvation on key TSs is also explored. On a fundamental level, this analysis will aid in understanding the thermodynamic and kinetic characteristics of the key intramolecular interactions involved in deamidation, dehydration, and other important fragmentations of peptides.

Unimolecular dissociation of peptides: statistical vs. non-statistical fragmentation mechanisms and time scales

In the present work we have investigated mechanisms of gas phase unimolecular dissociation of a relatively simple dipeptide, the di-proline anion, by means of chemical dynamics simulations, using the PM3 semi-empirical Hamiltonian. In particular, we have considered two activation processes that are representative limits of what occurs in collision induced dissociation experiments: (i) thermal activation, corresponding to several low energy collisions, in which the system is prepared with a microcanonical distribution of energy; (ii) collisional activation where a single shock of hundreds of kcal mol⁻¹ (300 kcal mol⁻¹ in the present case) can transfer sufficient energy to allow dissociation. From these two activation processes we obtained different product abundances, and for one particular fragmentation pathway a clear mechanistic difference for the two activation processes. This mechanism corresponds to the leaving of an OH⁻ group and subsequent formation of water by taking a proton from the remaining molecule. This last reaction is always observed in thermal activation while in collisional activation it is less favoured and the formation of OH⁻ as a final product is observed. More importantly, we show that while in thermal activation unimolecular dissociation follows exponential decay, in collision activation the initial population decays with non-exponential behaviour. Finally, from the thermal activation simulations it was possible to obtain rate constants as a function of temperature that show Arrhenius behaviour. Thus activation energies have also been extracted from these simulations.

The protonated and sodiated dimers of proline studied by IRMPD spectroscopy in the N-H and O-H stretching region and computational methods

IRMPD spectroscopy and computational chemistry techniques have been used to determine that the proton- and sodium-bound dimers of proline exist as a mixture of a number of different structures. Simulated annealing computations were found to be helpful in determining the unique structures of the protonated and sodiated dimers, augmenting chemical intuition. The experimental and computational results are consistent with the proton-bound dimer of N-protonated proline bound to zwitterionic proline. There was no spectroscopic evidence in the 3200-3800 cm(-1) region for a canonical structure which is predicted to have a weak N-H stretch at about 3440 cm(-1). A well resolved band at 1733 cm(-1) from a previous spectroscopic study (DOI: 10.1021/ja068715a ) was reassigned from a high energy canonical isomer to the C=O stretch of a lower energy zwitterionic structure. This band is a free carboxylate C=O stretch where protonated proline is hydrogen bonded to the other carboxylate oxygen which is also involved in an intramolecular hydrogen bond. Fifteen structures of the sodium bound proline dimer were computed to be within 10 kJ mol(-1) of Gibbs energy and eight structures were within 5 kJ mol(-1). None of these structures can be ruled out based on the experimental IRMPD spectrum. They all have an N-H stretching band predicted in a position that agrees with the experimental spectrum. However, only structures where one of the proline monomers is in the canonical form and having a free O-H bond can produce the band at ∼3600 cm(-1).

IR action spectroscopy shows competitive oxazolone and diketopiperazine formation in peptides depends on peptide length and identity of terminal residue in the departing fragment

The interplay between the entropically and enthalpically favored products of peptide fragmentation is probed using a combined experimental and theoretical approach. These b2 ion products can take either an oxazolone or diketopiperazine structure. Cleavage after the second amide bond is often a favorable process because the products are small ring structures that are particularly stable. These structures are structurally characterized by action IRMPD spectroscopy and semi-quantified using gas-phase hydrogen-deuterium exchange. The formation of the oxazolone and diketopiperazine has been thought to be largely governed by the identity of the first two residues at the N-terminus of the peptide. We show here that the length of the precursor peptide and identity of the third residue play a significant role in the formation of the diketopiperazine structure in peptides containing an N-terminal asparagine residue. This is additionally the first instance showing an N-terminal residue with an amide side chain can promote formation of the diketopiperazine b2 ion structure.

Thermodynamics and mechanism of protonated cysteine decomposition: a guided ion beam and computational study

A quantitative molecular description of the decomposition of protonated cysteine, H(+)Cys, is provided by studying the kinetic energy dependence of threshold collision-induced dissociation (CID) with Xe using a guided ion beam tandem mass spectrometer (GIBMS). Primary dissociation channels are deamidation (yielding both NH3 loss and NH4(+) formation) and (H2O + CO) loss reactions, followed by an additional six subsequent decompositions. Analysis of the kinetic energy-dependent CID cross sections provides the 0 K barriers for six different reactions after accounting for unimolecular decay rates, internal energy of reactant ions, multiple ion-molecule collisions, and competition among the decay channels. To identify the mechanisms associated with these reactions, quantum chemical calculations performed at the B3LYP/6-311 + G(d,p) level were used to locate the transition states (TSs) and intermediates for these processes. Single point energies of the reactants, products, and key optimized TSs and intermediates are calculated at B3LYP, B3P86, and MP2(full) levels using a 6-311 + G(2d,2p) basis set. The computational characterization of the elementary steps of these reactions, including the structures of the final products, is validated by quantitative agreement with the experimental energetics. In agreement with previous work, deamidation is facilitated by anchimeric assistance of the thio group, which also leads to an interesting rearrangement of the intact amino acid identified computationally.

From the mobile proton to wandering hydride ion: mechanistic aspects of gas-phase ion chemistry

Structural characterization of molecular species by mass spectrometry supposes the knowledge of the type of ions generated and the mechanism by which they dissociate. In this context, a need for a rationalization of electrospray ionization(+)(-) mass spectra of small molecules has been recently expressed. Similarly, at the other end of the mass scale, efforts are currently made to interpret the major fragmentation processes of protonated and deprotonated peptides and their reduced forms produced in electron capture or electron transfer experiments. Most fragmentation processes of molecular and pseudo-molecular ions produced in the ion source of a mass spectrometer may be described by a combination of several key mechanistic steps: simple bond dissociation, formation of ion-neutral complex intermediates, hydrogen atom, hydride ion or proton migrations and nucleophilic attack. Selected crucial aspects of these elementary reactions, occurring inside positively charged ions, will be recalled and illustrated by examples taken in recent mass spectrometry literature. Emphasis will be given on the protonation process and its consequence in terms of structure and energetic.

The power of accurate energetics (or thermochemistry: what is it good for?)

The utility of measuring the energetics of ion-molecule reactions is discussed. After distinguishing between the terms of thermodynamics (macroscopic, equilibrium quantities) and energetics (microscopic and kinetically relevant quantities), the potential energy surfaces for ion-molecule reactions are reviewed and their implications discussed. Equations describing the kinetic energy dependence of ion-molecule reactions are introduced and the effects of entropy on reaction rates and branching ratios are discussed. Several case histories allow an exploration of the utility of accurate thermochemical information and probe how accurate such energetic information must be to be predictive. These case studies include decomposition of hydrated metal dications, the reaction of FeO(+) with H(2), and fragmentation of a small protonated peptide (GG). These illustrate a range of interesting systems for which accurate energetic information has been influential in understanding the observed reactivity. Comparisons with theory demonstrate that experimental information is still required for truly predictive capability.

Using dissociation energies to predict observability of b- and y-peaks in mass spectra of short peptides. II. Results for hexapeptides with non-polar side chains

RATIONALE

The hypothesis that dissociation energies can serve as a predictor of observability of b- and y-peaks is tested for seven hexapeptides. If the hypothesis holds true for large classes of peptides, one would be able to improve the scoring accuracy of peptide identification tools by excluding theoretical peaks that cannot be observed in practical product ion spectra due to various physical, chemical or thermodynamic considerations.

METHODS

Product ion m/z spectra of hexapeptides AAAAAA, AAAFAA, AAAVAA, AAFAAA, AAVAAA, AAFFAA and AAVVAA have been acquired on a Finnigan LTQ XL mass spectrometer in the collision-induced dissociation (CID) activation mode on a grid of activation times 0.05 to 100 ms and normalized collision energy 10 to 35%. Dissociation energies were calculated for all fragmentation channels leading to b- and y-fragments at the TPSS/6-31G(d,p) level of the density functional theory.

RESULTS

It was demonstrated that the m/z peaks observed in the product ion spectra correspond to the fragmentation channels with dissociation energies below a certain threshold value. However, there is no direct correlation between the most intense m/z peaks and the lowest dissociation energies. Using the dissociation energies, it was predicted that out of 63 theoretically possible peaks in the b- and y-series of the seven hexapeptides, 19 should not be observable in practical spectra. In the experiments, 24 peaks were not observed, including all 19 predicted.

CONCLUSIONS

Dissociation energies alone are not sufficient for predicting ion intensity relationships in product ion m/z spectra. Nevertheless, the present data suggest that dissociation energies appear to be good predictors of observability of b- and y-peaks and potentially very useful for filtering theoretical peaks of each candidate peptide in peptide identification tools. Published 2012. This article is a US Government work and is in the public domain in the USA.