Structures of protonated arginine dimer and bradykinin investigated by density functional theory: further support for stable gas-phase salt bridges

Exploiting Charge State Distribution To Probe Intramolecular Interactions in Gas-Phase Phosphopeptides and Enhance Proteomics Analyses

Charging of analytes is a prerequisite for performing mass spectrometry analysis. In proteomics, electrospray ionization is the dominant technique for this process. Although the observation of differences in the peptide charge state distribution (CSD) is well-known among experimentalists, its analytical value remains underexplored. To investigate the utility of this dimension, we analyzed several public data sets, comprising over 250,000 peptide CSD profiles from the human proteome. We found that the dimensions of the CSD demonstrate high reproducibility across multiple laboratories, mass analyzers, and extensive time intervals. The general observation was that the CSD enabled effective partitioning of the peptide property space, resulting in enhanced discrimination between sequence and constitutional peptide isomers. Next, by evaluating the CSD values of phosphorylated peptides, we were able to differentiate between phosphopeptides that indicate the formation of intramolecular structures in the gas phase and those that do not. The reproducibility of the CSD values (mean cosine similarity above 0.97 for most of the experiments) qualified CSD data suitable to train a deep-learning model capable of accurately predicting CSD values (mean cosine similarity - 0.98). When we applied the CSD dimension to MS1- and MS2-based proteomics experiments, we consistently observed around a 5% increase in protein and peptide identification rate. Even though the CSD dimension is not as effective a discriminator as the widely used retention time dimension, it still holds the potential for application in direct infusion proteomics.

THE IMS PARADOX: A PERSPECTIVE ON STRUCTURAL ION MOBILITY-MASS SPECTROMETRY

Studies of large proteins, protein complexes, and membrane protein complexes pose new challenges, most notably the need for increased ion mobility (IM) and mass spectrometry (MS) resolution. This review covers evolutionary developments in IM-MS in the authors' and key collaborators' laboratories with specific focus on developments that enhance the utility of IM-MS for structural analysis. IM-MS measurements are performed on gas phase ions, thus "structural IM-MS" appears paradoxical-do gas phase ions retain their solution phase structure? There is growing evidence to support the notion that solution phase structure(s) can be retained by the gas phase ions. It should not go unnoticed that we use "structures" in this statement because an important feature of IM-MS is the ability to deal with conformationally heterogeneous systems, thus providing a direct measure of conformational entropy. The extension of this work to large proteins and protein complexes has motivated our development of Fourier-transform IM-MS instruments, a strategy first described by Hill and coworkers in 1985 (Anal Chem, 1985, 57, pp. 402-406) that has proved to be a game-changer in our quest to merge drift tube (DT) and ion mobility and the high mass resolution orbitrap MS instruments. DT-IMS is the only method that allows first-principles determinations of rotationally averaged collision cross sections (CSS), which is essential for studies of biomolecules where the conformational diversities of the molecule precludes the use of CCS calibration approaches. The Fourier transform-IM-orbitrap instrument described here also incorporates the full suite of native MS/IM-MS capabilities that are currently employed in the most advanced native MS/IM-MS instruments. © 2020 John Wiley & Sons Ltd. Mass Spec Rev.

Theoretical studies of infrared signatures of proton-bound amino acid dimers with homochiral and heterochiral moieties

Proton-bound homochiral and heterochiral dimers, X-H⁺ -X, of five amino acids (X = Ser, Ala, Thr, Phe, and Arg) are investigated theoretically using quantum chemical density functional theory (DFT) calculations and molecular dynamics simulations with the aim to unveil diastereomer-specific mid-infrared (mid-IR) absorption bands in the spectral range of 1000 to 1800 cm^-1 . The theoretical calculations performed in this work imply that all systems, except Ala₂ H⁺ , have distinct mid-IR absorption bands in homochiral and heterochiral configurations, which make them appropriate systems to be studied experimentally with mid-IR spectroscopy. We show that intermolecular interaction with the side chain, in the form of hydrogen bonding or cation-π interaction, is necessary for chiral effects to be present in the mid-IR spectra of proton-bound dimers of amino acids. We also report new conformers for Ala₂ H⁺ , Thr₂ H⁺ , Phe₂ H⁺ , and Arg₂ H⁺ , which were not found in earlier studies of these dimers.

Effect of chemical interaction between oleic acid and L-Arginine on oral perception, as a function of polymorphisms of CD36 and OBPIIa and genetic ability to taste 6-n-propylthiouracil

Oral sensitivity to fats varies in individuals influencing nutritional status and health. Variations in oleic acid perception are associated with CD36 and odorant binding protein (OBPIIa) polymorphisms, and 6-n-propylthiouracil (PROP) sensitivity, which is mediated by TAS2R38 receptor. L-Arginine (L-Arg) supplementation was shown to modify the perception of the five taste qualities. Here we analyzed the effect of three concentrations (5, 10, 15 mmol/L) of L-Arg on oral perception of oleic acid in forty-six subjects classified for PROP taster status and genotyped for TAS2R38, CD36 and OBPIIa polymorphisms. L-Arg supplementation was effective in increasing the perceived intensity of oleic acid in most subjects. The lowest concentration was the most effective, especially in PROP non-tasters or medium tasters, and in subjects with at least an allele A in CD36 and OBPIIa loci. Density Functional Theory (DFT) calculations were exploited to characterize the chemical interaction between L-Arg and oleic acid, showing that a stable 1:1 oleate·ArgH⁺ adduct can be formed, stabilized by a pair of hydrogen bonds. Results indicate that L-Arg, acting as a ‘carrier’ of fatty acids in saliva, can selectively modify taste response, and suggest that it may to be used in personalized dietetic strategies to optimize eating behaviors and health.

Competition between salt bridge and non-zwitterionic structures in deprotonated amino acid dimers

The effect of side chain functional groups on salt bridge structures in deprotonated amino acid homodimers is investigated using both infrared multiple photon dissociation spectroscopy between 650 and 1850 cm⁻¹ and theory.

Structure of protonated tryptophan dimer in the gas phase investigated by IRPD spectroscopy and theoretical calculations

RATIONAL

In addition to proton affinity (PA), side chains may also have an effect on the structures of amino acid complex ions in the gas phase. A previous study showed that the most stable isomer of Pro2 H(+) favored a salt-bridged structure. Tryptophan has a PA close to that of proline, but a quite different side chain. Thus, it will be interesting to discover whether the salt-bridged or charge-solvated conformation is energetically more favorable for Trp2 H(+) in the gas phase.

METHODS

The infrared photodissociation (IRPD) spectrum of Trp2 H(+) was obtained using a Fourier transform ion cyclotron resonance mass spectrometer equipped with a tunable OPO laser. The non-covalent cluster ions were generated by electrospray ionization. Structural optimization and frequency calculation of the selected isomers were performed at the M062X/6-311++G(d,p) level.

RESULTS

The experimental IRPD spectrum of Trp2 H(+) was reported in the region of 2700-3750 cm(-1) . Theoretical calculations show that the most stable isomer has a charge-solvated structure. Its energy was found to be 9 kcal/mol lower than that of the most stable salt-bridged isomer. The experimental spectrum is consistent with the predicted spectra of the most stable charge-solvated structures. Temperature effect on the stability of isomers was also evaluated and it was revealed that the contribution from salt-bridged isomers can be neglected at a temperature of 300 K.

CONCLUSIONS

Combining the method of IRPD spectroscopy with theoretical calculations, the structures of Trp2 H(+) were investigated. It is shown that the structures of Trp2 H(+) are dominated by charge-solvated forms. The results also indicate that the side chain may considerably affect the stability of the zwitterionic forms. Copyright © 2016 John Wiley & Sons, Ltd.

Stereoselective recognition of the Ac-Glu-Tyr-OH dipeptide by pseudopeptidic cages

Pseudopeptidic molecular cages are appealing receptors since they can display different polar and non-polar interaction sites in a modular framework and a controlled disposition. Inspired by previous host-guest knowledge, two pseudopeptidic molecular cages based on serine and threonine (CySer and CyThr, respectively) were designed and synthesized as hosts for the binding of the four possible stereoisomers of the Ac-Glu-Tyr-OH dipeptide, a target sequence of tyrosine kinases. The careful NMR titration experiments in aqueous acetonitrile allowed the determination of the binding constants and reflected a difference in the stability of the corresponding diastereomeric host-guest complexes. The CySer cage proved to be slightly more efficient than the CyThr counterpart, although both showed similar stereoselectivity trends: LL > DD ≥ LD > DL. This stereoselective binding was retained in the gas phase, as shown by ESI-MS competition experiments using the enantiomer-labelled method (EL), as well as CID experiments. Thus, the MS-determined discriminations follow the same trends observed by NMR, suggesting that the stereoselectivity observed for these systems must be mainly dictated by the polar host-guest interactions. Despite the stereoselective binding of short peptide sequences in competitive media being a challenging issue in supramolecular chemistry, our results demonstrate the power of pseudopeptidic cages in molecular recognition with foreseen implications in chemical biology.

Suppression of peptide ion dissociation under electron capture: role of backbone amide hydrogen

RATIONALE

The electron capture dissociation (ECD) of proteins/peptides is affected by the nature and sequence of amino acid residues. Electron capture/transfer with no dissociation is an intriguing phenomenon that has occasionally been observed. We have previously identified that diarginated peptides enriched with glutamic acid residues were found to show suppression of backbone fragmentation. In this paper, we report the effect of geometrical parameters of a peptide, including chain length, conformation and amide hydrogen, on the suppression of ECD fragmentation using synthetic model peptides.

METHODS

Glycine containing model polypeptides were used to probe the mechanism. Molecular-mechanics was used to obtain the conformation of the precursor ions. The ECD experiments were performed on a Bruker APEX III 4.7 T Fourier transform ion cyclotron resonance (FTICR) mass spectrometer.

RESULTS

Significant decreases in the intensities of the fragment ions were observed for the 23-mer polypeptide with only one E residue. This implied that the E:R ratio was no longer the sole determining factor for the occurrence of suppression effects. Results of conformational searches showed that there was a hydrogen-bonding 'ladder' formed in the 23-mer polypeptide, which was not found in the 15-mer peptide. Substituting the normal amino acid residues by the corresponding N-methylated amino acid residues in the model peptide, the suppression effect disappeared.

CONCLUSIONS

Our results indicate that survival of the intact reduced peptide ion after electron capture depends also on the length of the peptide. The amide hydrogen was critical in forming the resonance structure that suppressed the ECD fragmentation.

Structure of Protonated Threonine Dimers in the Gas Phase: Salt-Bridged or Charge-Solvated?

For homodimers of amino acids, their salt-bridged structures are gradually stabilized as the proton affinity of the component amino acid increases. Threonine has a proton affinity value located in the middle of the list of 20 natural amino acids. Thus, identifying whether the most stable isomer of protonated threonine dimer (Thr(2)H(+)) has a charge-solvated or salt-bridged structure is important and helpful for understanding the structures of other homodimers. By combining infrared photodissociation (IRPD) spectroscopy and theoretical calculations, the structures of Thr(2)H(+) were investigated. Based on calculations at the M062X/6-311++G(d,p)//M062X/6-311++G(d,p) level, the most stable isomer of Thr(2)H(+) was computed to be a charge-solvated structure, with an energy 3.87 kcal/mol lower than the most stable salt-bridged isomer. The predicted infrared spectrum is in good agreement with the experimental spectrum. To evaluate the temperature effect on the distribution of different isomers, the relative concentrations of the six isomers of Thr(2)H(+) were calculated at different temperatures, according to their partition functions and enthalpies. The results show that the isomers are dominated by charge-solvated structures at a temperature of 300 K.

Dose-Dependent Effects of L-Arginine on PROP Bitterness Intensity and Latency and Characteristics of the Chemical Interaction between PROP and L-Arginine

Genetic variation in the ability to taste the bitterness of 6-n-propylthiouracil (PROP) is a complex trait that has been used to predict food preferences and eating habits. PROP tasting is primarily controlled by polymorphisms in the TAS2R38 gene. However, a variety of factors are known to modify the phenotype. Principle among them is the salivary protein Ps-1 belonging to the basic proline-rich protein family (bPRP). Recently, we showed that oral supplementation with Ps-1 as well as its related free amino acids (L-Arg and L-Lys) enhances PROP bitterness perception, especially for PROP non-tasters who have low salivary levels of Ps-1. Here, we show that salivary L-Arg levels are higher in PROP super-tasters compared to medium tasters and non-tasters, and that oral supplementation with free L-Arg enhances PROP bitterness intensity as well as reduces bitterness latency in a dose-dependent manner, particularly in individuals with low salivary levels of both free L-Arg and Ps-1 protein. Supplementation with L-Arg also enhanced the bitterness of caffeine. We also used 1H-NMR spectroscopy and quantum-mechanical calculations carried out by Density Functional Theory (DFT) to characterize the chemical interaction between free L-Arg and the PROP molecule. Results showed that the -NH2 terminal group of the L-ArgH+ side chain interacts with the carbonyl or thiocarbonyl groups of PROP by forming two hydrogen bonds with the resulting charged adduct. The formation of this PROP•ArgH+ hydrogen-bonded adduct could enhance bitterness intensity by increasing the solubility of PROP in saliva and its availability to receptor sites. Our data suggest that L-Arg could act as a 'carrier' of various bitter molecules in saliva.

The use of chromium(III) to supercharge peptides by protonation at low basicity sites

The addition of chromium(III) nitrate to solutions of peptides with seven or more residues greatly increases the formation of doubly protonated peptides, [M + 2H](2+), by electrospray ionization. The test compound heptaalanine has only one highly basic site (the N-terminal amino group) and undergoes almost exclusive single protonation using standard solvents. When Cr(III) is added to the solution, abundant [M + 2H](2+) forms, which involves protonation of the peptide backbone or the C-terminus. Salts of Al(III), Mn(II), Fe(III), Fe(II), Cu(II), Zn (II), Rh(III), La(III), Ce(IV), and Eu(III) were also studied. Although several metal ions slightly enhance protonation, Cr(III) has by far the greatest ability to generate [M + 2H](2+). Cr(III) does not supercharge peptide methyl esters, which suggests that the mechanism involves interaction of Cr(III) with a carboxylic acid group. Other factors may include the high acidity of hexa-aquochromium(III) and the resistance of Cr(III) to reduction. Nitrate salts enhance protonation more than chloride salts and a molar ratio of 10:1 Cr(III):peptide produces the most intense [M + 2H](2+). Cr(III) also supercharges numerous other small peptides, including highly acidic species. For basic peptides, Cr(III) increases the charge state (2+ versus 1+) and causes the number of peptide molecules being protonated to double or triple. Chromium(III) does not supercharge the proteins cytochrome c and myoglobin. The ability of Cr(III) to enhance [M + 2H](2+) intensity may prove useful in tandem mass spectrometry because of the resulting overall increase in signal-to-noise ratio, the fact that [M + 2H](2+) generally dissociate more readily than [M + H](+), and the ability to produce [M + 2H](2+) precursors for electron-based dissociation techniques.

Statistical analysis of ion mobility spectrometry. II. Adaptively biased methods and shape correlations

Following a recent effort [J. Am. Soc. Mass Spectrom. 23, 386-396 (2012)], we continue to explore computational methodologies for generating molecular conformations to support collisional cross sections suggested by ion mobility measurements. Here, adaptively biased molecular dynamics (ABMD) simulations are used to sample the configuration space and to achieve flat-histogram sampling along the reaction coordinates of the first two moments of the gyration tensor. The method is tested and compared with replica-exchange simulations on triply-protonated bradykinin and on a larger 25-residue peptide. It is found to have a much higher efficiency for producing large sets of conformations in a broad range of diffusion cross-sections, whereas it does not compete with conventional replica-exchange molecular dynamics in locating the lowest-energy structure. Nevertheless, the broad sampling obtained from the ABMD method allows to quantitatively correlate the diffusion cross-section Ω with other geometric order parameters that have simpler interpretation. The strong correlations found between the diffusion cross-section and the radius of gyration, the surface area and the volume of the convex hull suggest an optimal template for accurately mimicking the variations of Ω in a broad range of conformations, using only geometrical information and doing so at a very moderate computational cost. The existence of such a correlation is confirmed on the much larger protein α-lactalbumin.

Statistical analysis of ion mobility spectrometry. I. Unbiased and guided replica-exchange molecular dynamics

Achieving (bio)macromolecular structural assignment from the interpretation of ion mobility spectrometry (IMS) experiments requires successful comparison with computer modeling. Replica-exchange molecular dynamics simulations with suitable force fields not only offer a convenient framework to locate relevant conformations, especially in the case of multiple-funnel energy landscapes, but they are also well suited to statistical analyses. In the present paper, we discuss two extensions of the method used to improve its efficiency in the context of IMS. Two doubly-protonated polyalanines [RA(4)XA(4)K + 2H](2+) with X=V and D appear as favorable cases for which the calculated collision cross-section distributions naturally agree with the measurements, providing reliable candidate structures. For these compounds, a careful consideration of other order parameters based on the weighted histogram method resolves several otherwise hidden underlying conformational families. In the case of a much larger peptide exhibiting bistability, assignment is more difficult but could be achieved by guiding the sampling with an umbrella potential using the square gyration radius as the biasing coordinate. Applied to triply protonated bradykinine, the two presented methods indicate that different conformations compatible with the measurements are very close in energy.

Ions in size-selected aqueous nanodrops: sequential water molecule binding energies and effects of water on ion fluorescence

The effects of water on ion fluorescence were investigated, and average sequential water molecule binding energies to hydrated ions, M(z)(H(2)O)(n), at large cluster size were measured using ion nanocalorimetry. Upon 248-nm excitation, nanodrops with ~25 or more water molecules that contain either rhodamine 590(+), rhodamine 640(+), or Ce(3+) emit a photon with average energies of approximately 548, 590, and 348 nm, respectively. These values are very close to the emission maxima of the corresponding ions in solution, indicating that the photophysical properties of these ions in the nanodrops approach those of the fully hydrated ions at relatively small cluster size. As occurs in solution, these ions in nanodrops with 8 or more water molecules fluoresce with a quantum yield of ~1. Ce(3+) containing nanodrops that also contain OH(-) fluoresce, whereas those with NO(3)(-) do not. This indirect fluorescence detection method has the advantages of high sensitivity, and both the size of the nanodrops as well as their constituents can be carefully controlled. For ions that do not fluoresce in solution, such as protonated tryptophan, full internal conversion of the absorbed 248-nm photon occurs, and the average sequential water molecule binding energies to the hydrated ions can be accurately obtained at large cluster sizes. The average sequential water molecule binding energies for TrpH(+)(H(2)O)(n) and a doubly protonated tripeptide, [KYK + 2H](2+)(H(2)O)(n), approach asymptotic values of ~9.3 (n ≥ 11) and ~10.0 kcal/mol (n ≥ 25), respectively, consistent with a liquidlike structure of water in these nanodrops.

On the zwitterionic nature of gas-phase peptides and protein ions

Determining the total number of charged residues corresponding to a given value of net charge for peptides and proteins in gas phase is crucial for the interpretation of mass-spectrometry data, yet it is far from being understood. Here we show that a novel computational protocol based on force field and massive density functional calculations is able to reproduce the experimental facets of well investigated systems, such as angiotensin II, bradykinin, and tryptophan-cage. The protocol takes into account all of the possible protomers compatible with a given charge state. Our calculations predict that the low charge states are zwitterions, because the stabilization due to intramolecular hydrogen bonding and salt-bridges can compensate for the thermodynamic penalty deriving from deprotonation of acid residues. In contrast, high charge states may or may not be zwitterions because internal solvation might not compensate for the energy cost of charge separation.

Factors that influence helical preferences for singly charged gas-phase peptide ions: the effects of multiple potential charge-carrying sites

Ion mobility-mass spectrometry is used to investigate the structure(s) of a series of model peptide [M + H](+) ions to better understand how intrinsic properties affect structure in low dielectric environments. The influence of peptide length, amino acid sequence, and composition on gas-phase structure is examined for a series of model peptides that have been previously studied in solution. Collision cross sections for the [M + H](+) ions of Ac-(AAKAA)(n)Y-NH(2) (n = 3-6) and Ac-Y(AEAAKA)(n)F-NH(2) (n = 2-5) are reported and correlated with candidate structures generated using molecular modeling techniques. The [M + H](+) ions of the AAKAA peptide series each exhibit a single, dominant ion mobility arrival time distribution (ATD) which correlates to partial helical structures, whereas the [M + H](+) ions of the AEAAKA ion series are composed of ATDs which correlate to charge-solvated globules (i.e., the charge is coordinated or solvated by polar peptide functional groups). These data raise numerous questions concerning intrinsic properties (amino acid sequence and composition as well as charge location) that dictate gas-phase peptide ion structure, which may reflect trends for peptide ion structure in low dielectric environments, such as transmembrane segments.

A mechanistic study of the H/D exchange reactions of protonated arginine and arginine-containing di- and tripeptides

The gas-phase H/D exchange reactions of arginine (R) and arginine-containing di- and tri-peptide (gly-arg (GR), arg-gly (RG), gly-gly-arg (GGR), gly-arg-gly (GRG) and arg-gly-gly (RGG)) [M + H]+ ions with deuterated ammonia (ND3) were investigated by using Fourier-transform ion cyclotron resonance mass spectrometry (FT-ICR), ion mobility-mass spectrometry (IM-MS), ab initio and density functional theory-based molecular orbital calculations and molecular modeling. Three exchanges are observed for arginine and arginine-containing tri-peptide [M + H]+ ions, whereas the di-peptide [M + H]+ ions undergo a single H/D exchange. In addition, C-terminal methylation blocks H/D exchange of arginine and the arginine-containing peptide [M + H]+ ions, and a single H/D exchange is observed for N-terminal acetylated arginine [M + H]+ ions. A general mechanism for H/D exchange involving a collision complex that is best described as a "solvated salt-bridge" structure is proposed.

Modifying the charge state distribution of proteins in electrospray ionization mass spectrometry by chemical derivatization

Electrospray ionization (ESI) of denatured proteins produces a broad distribution of multiply-charged ions leading to multiple peaks in the mass spectrum. We investigated changes in the positive-mode ESI charge state distribution produced by several chemical modifications of denatured proteins. Capping carboxylic acid groups with neutral functional groups yields little change in charge state distribution compared with unmodified proteins. The results indicate that carboxyl groups do not play a significant role in the positive charging of denatured proteins in ESI. The modification of proteins with additional basic sites or fixed positive charges generates substantially higher charge states, providing evidence that the number of ionizable sites, rather than molecular size and shape, determines ESI charging for denatured proteins. Fixed charge modification also significantly reduces the number of protons acquired by a protein, in that the charge state envelope is not increased by the full number of fixed charges appended. This result demonstrates that Coulombic repulsion between positive charges plays a significant role in determining charge state distribution by affecting the gas-phase basicity of ionizable sites. Addition of fixed-charge moieties to a protein is a useful approach for shifting protein charge state distributions to higher charge states, and with further work, it may help limit the distribution of protein ions to fewer charge states.

A density functional theory (DFT) study on gas-phase proton transfer reactions of derivatized and underivatized peptide ions generated by matrix-assisted laser desorption ionization

In this study, classic molecular dynamics (MD) simulations followed by density functional theory (DFT) calculations are employed to calculate the proton transfer reaction enthalpy shifts for native and derivatized peptide ions in the MALDI plume. First, absolute protonation and deprotonation enthalpies are calculated for native peptides (RPPGF and AFLDASR), the corresponding hexyl esters and three common matrices alpha-cyano-4-hydroxycinnamic acid (4HCCA), 2,5-dihydroxybenzoic acid (DHB), and 6 aza-2-thiothymine (ATT). From the proton exchange reaction calculations, protonation and deprotonation of the neutral peptides are thermodynamically favorable in the gas phase as long as the corresponding protonated/deprotonated matrix ions are present in the plume. Moreover, the gain in proton affinity shown by the ester ions suggests that the increase in ion yield is likely to be related to an easier proton transfer from the matrix to the peptide.

Application of FT-ICR-MS for the study of proton-transfer reactions involving biomolecules

Fourier transform ion cyclotron resonance mass spectrometry, combined with modern ionization (fast atom bombardment , electrospray ionization, matrix-assisted laser desorption-ionization), fragmentation (collision-induced dissociation, surface-induced dissociation, one-photon ultraviolet photodissociation, infrared multiphoton dissociation, blackbody infrared radiative dissociation, electron-capture dissociation), and separation (high-performance liquid chromatography, liquid chromatography, capillary electrophoresis) techniques is now becoming one of the most attractive and frequently used instrumental platforms for gas-phase studies of biomolecules such as amino acids, bioamines, peptides, polypeptides, proteins, nucleobases, nucleosides, nucleotides, polynucleotides, nucleic acids, saccharides, polysaccharides, etc. Since it gives the possibilities to trap the ions from a few seconds up to thousands of seconds, it is often applied to study ion/molecule reactions in the gas phase, particularly proton-transfer reactions which provide important information on acid-base properties. These properties determine in part the three-dimensional structure of biomolecules, most of their intramolecular and intermolecular interactions, and consequently their biological activity. They also indicate the form (unionized, zwitterionic, protonated, or deprotonated) which the biomolecule may take in a nonpolar environment.

Infrared spectroscopy of cationized arginine in the gas phase: direct evidence for the transition from nonzwitterionic to zwitterionic structure

The gas-phase structures of protonated and alkali metal cationized arginine (Arg) and arginine methyl ester (ArgOMe) are investigated with infrared spectroscopy and ab initio calculations. Infrared spectra, measured in the hydrogen-stretch region, provide compelling evidence that arginine changes from its nonzwitterionic to zwitterionic form with increasing metal ion size, with the transition in structure occurring between lithium and sodium. For sodiated arginine, evidence for both forms is obtained from spectral deconvolution, although the zwitterionic form is predominant. Comparisons of the photodissociation spectra with spectra calculated for low-energy candidate structures provide additional insights into the detailed structures of these ions. Arg*Li+, ArgOMe*Li+, and ArgOMe*Na+ exist in nonzwitterionic forms in which the metal ion is tricoordinated with the amino acid, whereas Arg*Na+ and Arg*K+ predominately exist in a zwitterionic form where the protonated side chain donates one hydrogen bond to the N terminus of the amino acid and the metal ion is bicoordinated with the carboxylate group. Arg*H+ and ArgOMe*H+ have protonated side chains that form the same interaction with the N terminus as zwitterionic, alkali metal cationized arginine, yet both are unambiguously determined to be nonzwitterionic. Calculations indicate that for clusters with protonated side chains, structures with two strong hydrogen bonds are lowest in energy, in disagreement with these experimental results. This study provides new detailed structural assignments and interpretations of previously observed fragmentation patterns for these ions.

Conformational analysis of arginine in gas phase—A strategy for scanning the potential energy surface effectively

The determination of all possible low-lying energy conformers of flexible molecules is of fundamental interest for various applications. It necessitates a reliable conformational search that is able to detect all important minimum structures and calculates the energies on an adequate level of theory. This work presents a strategy to identify low-energy conformers using arginine as an example by means of a force-field based conformational search in combination with high-level geometry optimizations (RI-MP2/TZVPP+). The methods used for various stages in the conformational search strategy are shown and various pitfalls are discussed. We can show that electronic energies calculated on a DFT level of theory with standard exchange-correlation functionals strongly underestimate the intramolecular stabilization resulting from stacked orientations of the guanidine and carbonyl moiety of arginine due to the deficiency of DFT to describe dispersion effects. In this case by usage of electron correlation methods, low energy conformers comprising stacked arrangements that are counterintuitive become favorable.

Zwitterionic States in Gas-Phase Polypeptide Ions Revealed by 157-nm Ultra-Violet Photodissociation

A new method of detecting the presence of deprotonation and determining its position in gas-phase polypeptide cations is described. The method involves 157-nm ultra-violet photodissociation (UVPD) and is based on monitoring the losses of CO2 (44 Da) from electronically excited deprotonated carboxylic groups relative to competing COOH losses (45 Da) from neutral carboxylic groups. Loss of CO2 is a strong indication of the presence of a zwitterionic [(+)...(-)...(+)] salt bridge in the gas-phase polypeptide cation. This method provides a tool for studying, for example, the nature of binding within polypeptide clusters. Collision-activated dissociation (CAD) of decarboxylated cations localizes the position of deprotonation. Fragment abundances can be used for the semiquantitative assessment of the branching ratio of deprotonation among different acidic sites, however, the mechanism of the fragment formation should be taken into account. Cations of Trp-cage proteins exist preferentially as zwitterions, with the deprotonation position divided between the Asp9 residue and the C terminus in the ratio 3:2. The majority of dications of the same molecule are not zwitterions. Furthermore, 157-nm UVPD produces abundant radical cations M*+ from protonated molecules through the loss of a hydrogen atom. This method of producing M*+ ions is general and can be applied to any gas-phase peptide cation. The abundance of the molecular radical cations M*+ produced is sufficient for further tandem mass spectrometry (MS/MS), which, in the cases studied, yielded side-chain loss of a basic amino acid as the most abundant fragmentation channel together with some backbone cleavages.

Gas-phase ion/ion reactions of multiply protonated polypeptides with metal containing anions

Gas-phase reactions of multiply protonated polypeptides and metal containing anions represent a new methodology for manipulating the cationizing agent composition of polypeptides. This approach affords greater flexibility in forming metal containing ions than commonly used methods, such as electrospray ionization of a metal salt/peptide mixture and matrix-assisted laser desorption. Here, the effects of properties of the polypeptide and anionic reactant on the nature of the reaction products are investigated. For a given metal, the identity of the ligand in the metal containing anion is the dominant factor in determining product distributions. For a given polypeptide ion, the difference between the metal ion affinity and the proton affinity of the negatively charged ligand in the anionic reactant is of predictive value in anticipating the relative contributions of proton transfer and metal ion transfer. Furthermore, the binding strength of the ligand anion to charge sites in the polypeptide correlates with the extent of observed cluster ion formation. Polypeptide composition, sequence, and charge state can also play a notable role in determining the distribution of products. In addition to their usefulness in gas-phase ion synthesis strategies, the reactions of protonated polypeptides and metal containing anions represent an example of a gas-phase ion/ion reaction that is sensitive to polypeptide structure. These observations are noteworthy in that they allude to the possibility of obtaining information, without requiring fragmentation of the peptide backbone, about ion structure as well as the relative ion affinities associated with the reactants.

Gas-phase ion/ion reactions of multiply protonated polypeptides with metal containing anions

Gas-phase reactions of multiply protonated polypeptides and metal containing anions represent a new methodology for manipulating the cationizing agent composition of polypeptides. This approach affords greater flexibility in forming metal containing ions than commonly used methods, such as electrospray ionization of a metal salt/peptide mixture and matrix-assisted laser desorption. Here, the effects of properties of the polypeptide and anionic reactant on the nature of the reaction products are investigated. For a given metal, the identity of the ligand in the metal containing anion is the dominant factor in determining product distributions. For a given polypeptide ion, the difference between the metal ion affinity and the proton affinity of the negatively charged ligand in the anionic reactant is of predictive value in anticipating the relative contributions of proton transfer and metal ion transfer. Furthermore, the binding strength of the ligand anion to charge sites in the polypeptide correlates with the extent of observed cluster ion formation. Polypeptide composition, sequence, and charge state can also play a notable role in determining the distribution of products. In addition to their usefulness in gas-phase ion synthesis strategies, the reactions of protonated polypeptides and metal containing anions represent an example of a gas-phase ion/ion reaction that is sensitive to polypeptide structure. These observations are noteworthy in that they allude to the possibility of obtaining information, without requiring fragmentation of the peptide backbone, about ion structure as well as the relative ion affinities associated with the reactants.

Binding energies of water to lithiated valine: formation of solution-phase structure in vacuo

Dissociation kinetics for loss of a water molecule from hydrated ions of lithiated valine, alanine ethyl ester and betaine are determined using blackbody infrared radiative dissociation at temperatures between -60 and 110 degrees C. From master equation modeling of these data, values of the threshold dissociation energy are obtained for clusters containing one through three water molecules. By comparing the values for valine with its two isomers, one a model for the nonzwitterion structure, the other a model for the zwitterion structure, information about the structure of valine in these hydrated clusters is inferred. Structures, relative energies, and water binding energies for these ions are also calculated at the B3LYP/6-31++G** level of theory. With one water molecule, both experiment and theory indicate that valine is not a zwitterion and that the lithium ion coordinates with the amino nitrogen and the carbonyl oxygen (NO coordinated) and the water molecule interacts directly with the lithium ion. With two water molecules, the zwitterion and nonzwitterion structures are nearly isoenergetic, but the experiment clearly indicates a NO-coordinated nonzwitterion structure. With three water molecules, both the experimental data and theory indicate that the lithium ion binds to the carboxylate group of valine, i.e., valine is zwitterionic with three water molecules. The agreement between the experimentally determined and calculated binding energies is good for all the clusters, with deviations of <== 0.12 eV.

H-bonding cooperativity and energetics of alpha-helix formation of five 17-amino acid peptides

Five peptides, each containing 17 amino acids, have been completely geometrically optimized in their alpha-helical and beta-strand forms using a mixed DFT/AM1 procedure. B3LYP/D95** was used for the entire helical structures, while AM1 was initially used to optimize the side chains, followed by reoptimization at the DFT level. The energetic and structural results show (1) that the helices are favored over the strands by 29.5 to 37.4 kcal/mol; (2) that alkyl groups on the amino acid side chains favor helix formation even in the absence of solvent; (3) that C-H...O hydrogen bonds contribute to the relative stability of the helices that contain amino acids (val, leu and ile) with beta-hydrogens in their alkyl side chains; (4) that formation of these helices entails approximately 6.6 kcal/mol of strain within the backbone per hydrogen bond; and (5) that H-bond cooperativity is essential for the alpha-helix to become more stable than a corresponding beta-strand. This last observation strongly suggests that pairwise potentials are inadequate for modeling of peptides and proteins.

Correlation of secondary structures of bradykinin B1 receptor antagonists with their activity

The secondary structure of a bradykinin B(1)receptor antagonist B-10324 (F5C-Lys-(1)- Lys(0)-Arg(1)-Pro(2)- Hyp(3)-Gly(4)-CpG(5)- Ser(6)-DTic(7)-CpG(8)) was determined by NMR at 800MHz. The conformational data are compared with those obtained previously for two bradykinin B(1) receptor antagonists, namely B-9858 (Lys-(1)- Lys(0)-Arg(1)-Pro(2)- Hyp(3)-Gly(4)-Igl(5)- Ser(6)-DIgl(7)-Oic(8)) and B-10148 (Lys-(1)-Lys(0)-Arg(1)- Pro(2)-Hyp(3)-Gly(4)- Igl(5)-Ser(6)-DF5F(7)- Oic(8)). The abnormal amino acids are: Hyp, trans-4- hydroxyproline; Tic, 1, 2, 3, 4-tetrahydroisoquinoline-3-carboxylic acid; Oic, (2S, 3aS, 7aS)-octahydroindole-2-carboxylic acid; Igl, alpha(2- indanyl)glycine; F5F, 2,3,4,5,6-pentafluorophenylalanine; CpG, alpha- cyclopentylglycine. F5C, pentafluorocinnamoyl, is the N-terminal protecting group and is not involved in the peptide secondary structure. B-10324 contains an N-terminal Pro(2)- CpG(5) distorted type II beta-turn whereas the rest of the peptide is random. A salt bridge is not observed between the carboxylate group at the C-terminal end and the Arg(1) side chain, in contrast to that previously observed for B-9858 and B- 10148. The conformations are correlated with the measured B(1) receptor antagonist activities (J.-F. Larrivée, L. Gera, S. Houle, J. Bouthillier, D. R. Bachvarov, J. M. Stewart and F. Marc au, Br. J. Pharmacol. 131, 885-892 (2000)). The importance of the N-terminal beta-turn is highlighted.

Gas-phase basicities for ions from bradykinin and its des-arginine analogues

Apparent gas-phase basicities (GB(app)s) for [M + H]+ of bradykinin, des-Arg1-bradykinin and des-Arg9-bradykinin have been assigned by deprotonation reactions of [M + 2H]2+ in a Fourier transform ion cyclotron resonance mass spectrometer. With a GB(app) of 225.8 +/- 4.2 kcal x mol(-1), bradykinin [M + H]+ is the most basic of the ions studied. Ions from des-Arg1-bradykinin and des-Arg9-bradykinin have GB(app) values of 222.8 +/- 4.3 kcal x mol(-1) and 214.9 +/- 2.3 kcal x mol(-1), respectively. One purpose of this work was to determine a suitable reaction efficiency 'break point' for assigning GB(app) values to peptide ions using the bracketing method. An efficiency value of 0.1 (i.e. approximately 10% of all collisions resulting in a deprotonation reaction) was used to assign GB(app)s. Support for this criterion is provided by the fact that our GB(app) values for des-Arg1-bradykinin and des-Arg9-bradykinin are identical, within experimental error, to literature values obtained using a modified kinetic method. However, the GB(app)s for bradykinin ions from the two studies differ by 10.3 kcal x mol(-1). The reason for this is not clear, but may involve conformation differences produced by experimental conditions. The results may be influenced by salt-bridge conformers and/or by conformational changes caused by the use of a proton-bound heterodimer in the kinetic method. Factors affecting the basicities of these peptide ions are also discussed, and molecular modeling is used to provide information on protonation sites and conformations. The presence of two highly basic arginine residues on bradykinin results in its high GB(app), while the basicity of des-Arg1-bradykinin ions is increased by the presence of two proline residues at the N-terminus. The proline residue in the second position folds the peptide chain in a manner that increases intramolecular hydrogen bonding to the protonated N-terminal amino group of the proline at the first position.

Structure of cationized glycine, gly.m (m = be, mg, ca, sr, ba), in the gas phase: intrinsic effect of cation size on zwitterion stability

Interactions between divalent metal ions and biomolecules are common both in solution and in the gas phase. Here, the intrinsic effect of divalent alkaline earth metal ions (Be, Mg, Ca, Sr, Ba) on the structure of glycine in the absence of solvent is examined. Results from both density functional and Moller-Plesset theories indicate that for all metal ions except beryllium, the salt-bridge form of the ion, in which glycine is a zwitterion, is between 5 and 12 kcal/mol more stable than the charge-solvated structure in which glycine is in its neutral form. For beryllium, the charge-solvated structure is 5-8 kcal/mol more stable than the salt-bridge structure. Thus, there is a dramatic change in the structure of glycine with increased metal cation size. Using a Hartree-Fock-based partitioning method, the interaction between the metal ion and glycine is separated into electrostatic, charge transfer and deformation components. The charge transfer interactions are more important for stabilizing the charge-solvated structure of glycine with beryllium relative to magnesium. In contrast, the difference in stability between the charge-solvated and salt-bridge structure for magnesium is mostly due to electrostatic interactions that favor formation of the salt-bridge structure. These results indicate that divalent metal ions dramatically influence the structure of this simplest amino acid in the gas phase.