Calculations of relative intensities of fragment ions in the MSMS spectra of a doubly charged penta-peptide

Size Dependent Fragmentation Chemistry of Short Doubly Protonated Tryptic Peptides

Tandem mass spectrometry of electrospray ionized multiply charged peptide ions is commonly used to identify the sequence of peptide(s) and infer the identity of source protein(s). Doubly protonated peptide ions are consistently the most efficiently sequenced ions following collision-induced dissociation of peptides generated by tryptic digestion. While the broad characteristics of longer (N ≥ 8 residue) doubly protonated peptides have been investigated, there is comparatively little data on shorter systems where charge repulsion should exhibit the greatest influence on the dissociation chemistry. To address this gap and further understand the chemistry underlying collisional-dissociation of doubly charged tryptic peptides, two series of analytes ([G_xR+2H]²⁺ and [A_xR+2H]²⁺, x = 2-5) were investigated experimentally and with theory. We find distinct differences in the preference of bond cleavage sites for these peptides as a function of size and to a lesser extent composition. Density functional calculations at two levels of theory predict that the threshold relative energies required for bond cleavages at the same site for peptides of different size are quite similar (for example, b₂-y_N-2). In isolation, this finding is inconsistent with experiment. However, the predicted extent of entropy change of these reactions is size dependent. Subsequent RRKM rate constant calculations provide a far clearer picture of the kinetics of the competing bond cleavage reactions enabling rationalization of experimental findings. The M06-2X data were substantially more consistent with experiment than were the B3LYP data.

N-Protonated Isomers and Coulombic Barriers to Dissociation of Doubly Protonated Ala8Arg

Collision-induced dissociation (or tandem mass spectrometry, MS/MS) of a protonated peptide results in a spectrum of fragment ions that is useful for inferring amino acid sequence. This is now commonplace and a foundation of proteomics. The underlying chemical and physical processes are believed to be those familiar from physical organic chemistry and chemical kinetics. However, first-principles predictions remain intractable because of the conflicting necessities for high accuracy (to achieve qualitatively correct kinetics) and computational speed (to compensate for the high cost of reliable calculations on such large molecules). To make progress, shortcuts are needed. Inspired by the popular mobile proton model, we have previously proposed a simplified theoretical model in which the gas-phase fragmentation pattern of protonated peptides reflects the relative stabilities of N-protonated isomers, thus avoiding the need for transition-state information. For singly protonated Ala _n (n = 3-11), the resulting predictions were in qualitative agreement with the results from low-energy MS/MS experiments. Here, the comparison is extended to a model tryptic peptide, doubly protonated Ala₈Arg. This is of interest because doubly protonated tryptic peptides are the most important in proteomics. In comparison with experimental results, our model seriously overpredicts the degree of backbone fragmentation at N₉. We offer an improved model that corrects this deficiency. The principal change is to include Coulombic barriers, which hinder the separation of the product cations from each other. Coulombic barriers may be equally important in MS/MS of all multiply charged peptide ions. Graphical Abstract ᅟ.

Predicting collision-induced dissociation spectra: semi-empirical calculations as a rapid and effective tool in software-aided mass spectral interpretation

RATIONALE

Fifteen molecules were modelled using quantum chemistry, prior to interpreting their collision-induced dissociation (CID) product ion spectra, in a 'blind trial' to establish if calculated protonation-induced bond elongation could be used to predict which bonds cleaved during CID. Bond elongation has the potential to be used as a descriptor predicting bond cleavage.

METHODS

The 15 molecules were modelled with respect to protonation-induced bond length changes using Density Functional Theory (DFT). Significant bond elongations were highlighted to flag potential bond cleavages. CID product ion spectra, obtained using positive ion electrospray ionisation (Waters Synapt G1), were interpreted to establish if observed bond cleavages correlated with calculated bond elongations. Calculations were also undertaken using AM1 (Austin Model 1) to see if this rapid approach gave similar results to the computationally demanding DFT.

RESULTS

The AM1-calculated bond elongations were found to be similar to those generated by DFT. All the polarised bonds observed to cleave (n = 82) had been calculated to elongate significantly. Protonation, possibly via proton migration, on the most electronegative atom in the bond appeared to initiate cleavage, leading to a 100% success rate in predicting the bonds that broke as a result of protonation on a heteroatom. Cleavage of carbon-carbon bonds was not predicted.

CONCLUSIONS

Cleavage of the polarised bonds appears to result from protonation on the more electronegative atom of the bond, inducing conformational changes leading to bond weakening. AM1-calculated bond length changes act as a descriptor for predicting bond cleavage. However, the impetus for cleavage of the unpolarised bonds may be product ion stability rather than bond weakening.