Gas-phase acidities of cysteine-polyglycine peptides: the effect of the cysteine position

Gas-Phase Dissociation Chemistry of Deprotonated RGD

We investigate the structure and dissociation pathways of the deprotonated amphoteric peptide arginylglycylasparic acid, [RGD-H]^-. We model the pertinent gas-phase structures and fragmentation chemistry of the precursor anions and predominant sequence-informative bond cleavages (b₂+H₂O, c₂, and z₁ peaks) and compare these predictions to our tandem mass spectra and infrared spectroscopy experiments. Formation of the b₂+H₂O anions requires rate-limiting intramolecular back biting to cleave the second amide bond and generate an anhydride structure. Facile cleavage of the newly formed ester bond with concerted expulsion of a cyclic anhydride neutral generates the product structure. IR spectroscopy supports this b₂+H₂O anion having structures that are essentially identical to C-terminally deprotonated arginylglycine, [RG-H]^-. Formation of the c₂ anion is predicted to require concerted expulsion of CO₂ from the aspartyl side chain carboxylate and cleavage of the N-C_alpha bond to produce a proton-bound dimer of arginylglycinamide and acrylate. Proton transfers within the dimer then enable predominant detection of a c₂ anion with the negative charge nominally on the central, glycine nitrogen (amidate structure) as the proton affinity of this structure is predicted to be lower than acrylate by ∼27 kJ mol^-1. Alternate means of cleaving the same N-C_alpha bond produce deprotonated cis-1,4-dibut-2-enoic acid z₁ anion structures. These lowest energy processes involve C-H proton mobilization from the aspartyl side chain prior to N-C_alpha bond cleavage consistent with proposals from the literature.

Comparison of the pH-dependent formation of His and Cys heptapeptide complexes of nickel(II), copper(II), and zinc(II) as determined by ion mobility-mass spectrometry

The analog methanobactin (amb) peptide with the sequence ac-His₁ -Cys₂ -Gly₃ -Pro₄ -Tyr₅ -His₆ -Cys₇ (amb_5A ) will bind the metal ions of zinc, nickel, and copper. To further understand how amb_5A binds these metals, we have undertaken a series of studies of structurally related heptapeptides where one or two of the potential His or Cys binding sites have been replaced by Gly, or the C-terminus has been blocked by amidation. The studies were designed to compare how these metals bind to these sequences in different pH solutions of pH 4.2 to 10 and utilized native electrospray ionization (ESI) with ion mobility-mass spectrometry (IM-MS) which allows for the quantitative analysis of the charged species produced during the reactions. The native ESI conditions were chosen to conserve as much of the solution-phase behavior of the amb peptides as possible and an analysis of how the IM-MS results compare with the expected solution-phase behavior is discussed. The oligopeptides studied here have applications for tag-based protein purification methods, as therapeutics for diseases caused by elevated metal ion levels or as inhibitors for metal-protein enzymes such as matrix metalloproteinases.

An Experimental and Computational Study of the Gas-Phase Acidities of the Common Amino Acid Amides

Using proton-transfer reactions in a Fourier transform ion cyclotron resonance mass spectrometer and correlated molecular orbital theory at the G3(MP2) level, gas-phase acidities (GAs) and the associated structures for amides corresponding to the common amino acids have been determined for the first time. These values are important because amino acid amides are models for residues in peptides and proteins. For compounds whose most acidic site is the C-terminal amide nitrogen, two ions populations were observed experimentally with GAs that differ by 4-7 kcal/mol. The lower energy, more acidic structure accounts for the majority of the ions formed by electrospray ionization. G3(MP2) calculations predict that the lowest energy anionic conformer has a cis-like orientation of the [-C(═O)NH](-) group whereas the higher energy, less acidic conformer has a trans-like orientation of this group. These two distinct conformers were predicted for compounds with aliphatic, amide, basic, hydroxyl, and thioether side chains. For the most acidic amino acid amides (tyrosine, cysteine, tryptophan, histidine, aspartic acid, and glutamic acid amides) only one conformer was observed experimentally, and its experimental GA correlates with the theoretical GA related to side chain deprotonation.

Kinetic Method Analysis of the Effect of cis- and trans-Hydroxylation on the Proton Affinity of Proline

Proton affinities of proline and hydroxyproline were measured using the Cooks’ kinetic method. The measurements show that hydroxylation increases the proton affinity, which is consistent with X3LYP computation results. This work supports findings from a previous study (S. Mezzache et al., Rapid Commun. Mass Spectrom. 2005, 19, 2279) that modification of proline increases its proton affinity, but it does not provide compelling evidence for the prediction in the same study that proton affinity of the molecule is substantially influenced by intramolecular interactions involving the proton. PBE1PBE calculations suggest that isotropic polarizability rather than intramolecular hydrogen-bonding interactions provides a more suitable diagnosis for trends in proton affinity changes associated with modifications.

Determination of the gas-phase acidities of oligopeptides

Amino acid residues located at different positions in folded proteins often exhibit different degrees of acidities. For example, a cysteine residue located at or near the N-terminus of a helix is often more acidic than that at or near the C-terminus (1-6). Although extensive experimental studies on the acid-base properties of peptides have been carried out in the condensed phase, in particular in aqueous solutions (6-8), the results are often complicated by solvent effects (7). In fact, most of the active sites in proteins are located near the interior region where solvent effects have been minimized (9,10). In order to understand intrinsic acid-base properties of peptides and proteins, it is important to perform the studies in a solvent-free environment. We present a method to measure the acidities of oligopeptides in the gas-phase. We use a cysteine-containing oligopeptide, Ala3CysNH2 (A3CH), as the model compound. The measurements are based on the well-established extended Cooks kinetic method (Figure 1) (11-16). The experiments are carried out using a triple-quadrupole mass spectrometer interfaced with an electrospray ionization (ESI) ion source (Figure 2). For each peptide sample, several reference acids are selected. The reference acids are structurally similar organic compounds with known gas-phase acidities. A solution of the mixture of the peptide and a reference acid is introduced into the mass spectrometer, and a gas-phase proton-bound anionic cluster of peptide-reference acid is formed. The proton-bound cluster is mass isolated and subsequently fragmented via collision-induced dissociation (CID) experiments. The resulting fragment ion abundances are analyzed using a relationship between the acidities and the cluster ion dissociation kinetics. The gas-phase acidity of the peptide is then obtained by linear regression of the thermo-kinetic plots (17,18). The method can be applied to a variety of molecular systems, including organic compounds, amino acids and their derivatives, oligonucleotides, and oligopeptides. By comparing the gas-phase acidities measured experimentally with those values calculated for different conformers, conformational effects on the acidities can be evaluated.

The applicability of the kinetic method for measuring relative affinities of macromolecules for polyatomic substrates

This paper is a review of the kinetic method for the determination of thermochemical values for gas-phase molecules. In addition, we have explored the utility of the kinetic method to obtain meaningful relative binding energies of macromolecules for polyatomic substrates using a system comprising poly(methylmethacrylate) (PMMA) oligomers and doubly protonated diaminoalkanes. The major factors which determined the suitability of the kinetic method for this system were identified as (i) the structural arrangement of the parent ion complex, (ii) possible reverse activation barriers, and (iii) the evaluations of Δ(ΔS‡). Molecular mechanics/molecular dynamics (MM/MD) simulations, together with ion mobility spectrometry, suggests the parent ion complexes represent a relatively equal sharing of the substrate between two the PMMA oligomers within the complex and that the two PMMA oligomers interact almost exclusively with the substrate, and not with each other. MS/MS of the trimeric parent complexes resulted in one PMMA unit leaving as a neutral which suggests very limited coulombic repulsion (that would contribute to a reverse activation barrier). The drift times of PMMA-diaminoalkane complexes that were generated directly by ESI-MS or by dissociation of a trimeric PMMA-diaminoalkane-PMMA complex were found to be identical, and when combined with MM/MD simulations suggested that the product PMMA-diaminoalkane dication has the same conformation as it does when part of a trimeric complex. This is evidence for Δ(ΔS‡) ≃ Δ(ΔS) and using a statistical mechanics approach, Δ(ΔS) ≃ 0. The effective temperature variable in the kinetic method expression was found to decrease as a function of the size of the trimeric complex, suggesting that the population distribution of the dissociating ensemble of complexes narrows as size increases.

Matrix-assisted laser desorption/ionisation mass spectrometric response factors of peptides generated using different proteolytic enzymes

Matrix-assisted laser desorption/ionisation (MALDI) mechanisms and the factors that influence the intensity of the ion signal in the mass spectrum remain imperfectly understood. In proteomics, it is often necessary to maximise the peptide response in the mass spectrum, especially for low abundant proteins or for proteolytic peptides of particular significance. We set out to determine which of the common proteolytic enzymes give rise to peptides with the best response factors under MALDI conditions. Standard proteins were enzymatically digested using four common proteases. We assessed relative response factors by coanalyzing the resulting digests. Thus, when tryptic peptides were added in equimolar quantities to their corresponding Asp-N, chymotrypsin and Glu-C digests, tryptic peptide signals were always predominant in the resulting MALDI mass spectra. Observable peaks attributable to non-tryptic peptides generally contained a terminal basic residue. It was proposed that a terminal basic residue has a disproportionate influence upon gas-phase basicity, and this hypothesis was supported by experiments with model isotopically labelled peptides. Experiments applying Cook's kinetic method showed that the peptide with a C-terminal arginine residue was more basic than the equivalent peptide with an N-terminal arginine, which was more basic than the peptide in which the arginine was mid-chain. Thus, the observation of the higher MALDI mass spectrometry response factors of tryptic peptides in comparison with peptides derived using other proteolytic enzymes corresponds with higher gas-phase basicities and may, along with other factors such as the complexity of the digest, influence the choice of enzyme in "bottom-up" proteomic experiments.

Conformational preferences and pK(a) value of selenocysteine residue

The conformational preferences of the L-selenocysteine (Sec) dipeptides with selenol and selenolate groups (Ac-Sec-NHMe and Ac-Sec(-) -NHMe, respectively) and the apparent (i.e., macroscopic) pK(a) value of the Sec residue have been studied using the dispersion-corrected density functionals M06-2X and B2PLYP-D with the implicit solvation method in the gas phase and in water. In the gas phase, the backbone-to-backbone and/or side chain-to-backbone hydrogen bonds are found to contribute in stabilizing the most preferred conformations for the Sec and Sec(-) residues, as seen for the Cys and Cys(-) residues. However, the polyproline II-like conformations prevail over the conformations with the backbone-to-backbone hydrogen bonds in water because of the weakened hydrogen bonds by the favorable direct interactions between the backbone CO and HN groups and water molecules. The Sec and Sec(-) residues are found to adopt more various conformations than the Cys and Cys(-) residues in water, although the most preferred conformations of the neutral and/or anionic forms of the two residues are similar each other in the gas phase and in water. Using the statistically weighted free energies of the Sec and Sec(-) dipeptides in the gas phase and their solvation free energies, the pK(a) value of the Sec residue is estimated to be 5.47 at 25°C, which is in good agreement with the experimental value of 5.43 ± 0.02. It is found that the lower pK(a) value of the selenol side chain for the Sec residue by ∼3 units than the thiol side chain for the Cys residue is ascribed to the higher gas-phase acidity of the Sec residue.