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.

Acetonitrile adduct analysis of underivatised amino acids offers improved sensitivity for hydrophilic interaction liquid chromatography tandem mass-spectrometry

LC-MS/MS method development for native amino acid detection can be problematic due to low ionisation efficiencies, in source fragmentation, potential for cluster ion formation and incorrect application of chromatography techniques. This has led to the majority of the scientific community derivatising amino acids for more sensitive analysis. Derivatisation has several benefits including reduced signal-to-noise ratios, more efficient ionisation, and a change in polarity, allowing the use of reverse phase chromatography. However, derivatisation of amino acids can be expensive, requires additional sample preparation steps, is more time consuming and increases sample instability, due to the most derivatised amino acids only be stable for finite amount of time. While showing initial promise, development of reliable hydrophilic interaction liquid chromatography (HILIC) separation methods has presented difficulties for the analyst including irreproducible separation and poor sensitivity. This study aimed to find a means to improve the detection sensitivity of the 20 protein amino acids by HILIC-MS/MS. We describe the use of previously undescribed amino acid-acetonitrile (ACN) adducts to improve detection of 16 out of the 20 amino acids. While all amino acids examined did form an ACN adduct, 4 had low intensity adduct formation compared to their protonated state, 3 of which are classified as basic amino acids. For 15 of the 20 amino acids tested, we used the ACN adduct for both quantification and qualification ions and demonstrated a significant enhancement in signal-to-noise ratio, ranging from 23 to 1762% improvement. Lower LODs, LOQs and lower ranges of linearity were also achieved for these amino acids. The optimised method was applied to a human neuroblastoma cell line (SH-SY5Y) with the potential to be applied to other complex sample types. The improved sensitivity this method offers simplifies sample preparation and reduces the costs of amino acid analysis compared to those methods that rely on derivatisation for sensitivity.

Sodium Binding Interactions with Aliphatic Amino Acids: A Guided Ion Beam and Computational Study

Metal binding affinities play a vital role in medicinal, biological, and industrial applications. In particular, metal cation-amino acid (AA) interactions contribute to protein stability such that analyzing analogous prototypical interactions is important. Here, we present a full description of the interactions of sodium cations (Na⁺) and six aliphatic amino acids (AA), where AA = glycine (Gly), alanine (Ala), homoalanine (hAla), valine (Val), leucine (Leu), and isoleucine (Ile). Experimentally, these interactions are evaluated using threshold collision-induced dissociation carried out in a guided ion beam tandem mass spectrometer, allowing for the determination of the kinetic-energy-dependent behavior of Na⁺-AA dissociation. Analysis of these dissociation cross sections, after accounting for multiple ion-molecule collisions, internal energy of reactant ions, and unimolecular decay rates, allows the determination of absolute Na⁺-AA bond dissociation energies (BDEs) in kJ/mol of Gly (164.0), Ala (166.9), hAla (167.9), Val (172.7), Leu (173.7), and Ile (174.6). These are favorably compared to quantum chemical calculations conducted at the B3LYP, B3P86, MP2(full), B3LYP-GD3BJ, and M06-2X levels of theory. Our combination of structural and energetic analyses provides information regarding the specific factors responsible for Na⁺ interactions with amino acids. Specifically, we find that the BDEs increase linearly with increasing polarizability of the amino acid.

Experimental and theoretical investigations of infrared multiple photon dissociation spectra of lysine complexes with Zn2+ and Cd2

The gas-phase structures of zinc and cadmium complexes of lysine (Lys) are investigated via a combination of infrared multiple photon dissociation action spectroscopy and ab initio quantum chemical calculations. In order to unambiguously identify the experimentally observed species, [Zn(Lys-H)]⁺ and CdCl⁺(Lys), the action spectra were compared to linear absorption spectra calculated at the B3LYP level of theory, using 6-311+G(d,p) and def2-TVZP basis sets for the zinc and cadmium systems, respectively. Single point energies were also calculated at the B3LYP, B3P86, MP2, and B3LYP-GD3BJ (accounting for empirical dispersion) levels of theory using larger basis sets. Identification of the experimentally formed isomers is possible through good agreement between infrared multiple photon dissociation action spectra and the theoretically predicted spectra. The [Zn(Lys-H)]⁺ complex adopts a tridentate orientation involving the amino acid backbone amine and deprotonated carboxylic acid groups as well as the side-chain amine group, [N_α,CO^-,N_ɛ]. The CdCl⁺(Lys) complex similarly adopts a tridentate chelation involving the amino acid backbone amine and carbonyl groups, as well as the side-chain amine group, [N_α,CO,N_ɛ]. In both cases, the identified complexes are the lowest energy gas-phase structures at all levels of theory.

Cationic Noncovalent Interactions: Energetics and Periodic Trends

In this review, noncovalent interactions of ions with neutral molecules are discussed. After defining the scope of the article, which excludes anionic and most protonated systems, methods associated with measuring thermodynamic information for such systems are briefly recounted. An extensive set of tables detailing available thermodynamic information for the noncovalent interactions of metal cations with a host of ligands is provided. Ligands include small molecules (H2, NH3, CO, CS, H2O, CH3CN, and others), organic ligands (O- and N-donors, crown ethers and related molecules, MALDI matrix molecules), π-ligands (alkenes, alkynes, benzene, and substituted benzenes), miscellaneous inorganic ligands, and biological systems (amino acids, peptides, sugars, nucleobases, nucleosides, and nucleotides). Hydration of metalated biological systems is also included along with selected proton-based systems: 18-crown-6 polyether with protonated peptides and base-pairing energies of nucleobases. In all cases, the literature thermochemistry is evaluated and, in many cases, reanchored or adjusted to 0 K bond dissociation energies. Trends in these values are discussed and related to a variety of simple molecular concepts.

Gas-phase lithium cation affinity of glycine

The gas-phase lithium cation binding thermochemistry of glycine has been determined theoretically by quantum chemical calculations at the G4 level and experimentally by the extended kinetic method using electrospray ionization quadrupole time-of-flight tandem mass spectrometry. The lithium cation affinity of glycine, ∆(Li)H°(298)(GLY), i.e. the∆(Li)H°(298) of the reaction GlyLi(+)→ Gly + Li(+)) given by the G4 method is equal to 241.4 kJ.mol(-1) if only the most stable conformer of glycine is considered or to 242.3 kJ.mol(-1) if the 298K equilibrium mixture of neutral conformers is included in the calculation. The ∆(Li)H°(298)(GLY) deduced from the extended kinetic method is obviously dependent on the choice of the Li(+) affinity scale, thus∆(Li)H°(298)(GLY) is equal to 228.7±0.9(2.0) kJ.mol(- 1) if anchored to the recently re-evaluated lithium cation affinity scale but shifted to 235.4±1.0 kJ.mol(-1) if G4 computed lithium cation affinities of the reference molecules is used. This difference of 6.3 kJ.mol(-1) may originate from a compression of the experimental lithium affinity scale in the high ∆(Li)H°(298) region. The entropy change associated with the reaction GlyLi(+)→Gly + Li(+) reveals a gain of approximately 15 J.mol(-) 1.K(-1) with respect to monodentate Li(+) acceptors. The origin of this excess entropy is attributed to the bidentate interaction between the Li(+) cation and both the carbonyl oxygen and the nitrogen atoms of glycine. The computed G4 Gibbs free energy,∆(Li)G°(298)(GLY) is equal to 205.3 kJ.mol(-1), a similar result, 201.0±3.4 kJ.mol(-1), is obtained from the experiment if the∆(Li)G°(298) of the reference molecules is anchored on the G4 results.

Desalting large protein complexes during native electrospray mass spectrometry by addition of amino acids to the working solution

A simple method for mitigating the adverse effects of salt adduction during native protein mass spectrometry by addition of amino-acids.

Extension of metal oxide laser ionization mass spectrometry to analytes with varied chemical functionalities

RATIONALE

Applications of metal oxide laser ionization mass spectrometry (MOLI MS) have been limited to compounds that provide protons for ionization upon adsorption to a metal oxide surface. The addition of a small molecule, which can act as a proton source and extend MOLI MS applications to a wider variety of analytes, was investigated.

METHODS

Mass spectrometric measurements were made with a PerSeptive Biosystems Voyager DE-STR MALDI mass spectrometer, using NiO-based MOLI to generate ions. NiO was pretreated in this study with a number of small molecules to act as proton sources during ionization. Peptides, monosaccharides, and a crude oil sample were analyzed with pretreated and untreated NiO.

RESULTS

Candidate small molecule proton donors were evaluated based on their effect on signal-to-noise ratio and peak intensity for tristearin. The best performing proton source was used for subsequent analyses. The analysis of monosaccharides, peptides and compounds in heavy crude oil was enhanced by the use of proton source treated NiO for MOLI MS. Proton/cation exchange following ionization was observed.

CONCLUSIONS

Addition of a proton source to NiO enhanced the MOLI MS response for peptides, carbohydrates, and compounds in a crude oil sample. Efficient ionization of these samples effectively extended applications of MOLI MS beyond lipid analytes containing esters.

DFT STUDY ON GAS-PHASE INTERACTION BETWEEN HISTIDINE AND ALKALI METAL IONS (Li+, Na+, K+); AND INFLUENCE OF THESE IONS ON HISTIDINE ACIDITY

The gas-phase metal affinities of histidine Li +, Na + and K + ions have been determined theoretically employing the hybrid B3LYP exchange–correlation functional and using 6-311++G** basis sets. All computations indicate that the metal ion affinity decreases on going from Li + to Na + and K + for the considered amino acid. Different types of M + coordinations on several histidine conformers/tautomers were considered in detail.

The optimized structures indicate that Li + and Na + prefer a tri-dentate coordination, bonding with a nitrogen atom of imidazole ring ( N τ), – NH 2, and an oxygen atom of a carbonyl, while in the K +-histidine lowest-energy conformer, the cation appears to be bi-coordinated to both oxygen atoms of the zwitterionic form by the energy values not too far from that of tri-coordination. We also performed the DFT calculations for proton dissociation energy of histidine both in the presence and absence of alkali metal ions. Our results also reveal that the acidity of histidine dramatically increases upon metal ion complexation.

Encapsulation of metal cations by the PhePhe ligand: a cation-π ion cage

Structures and binding thermochemistry are investigated for protonated PhePhe and for complexes of PhePhe with the alkaline-earth ions Ba(2+) and Ca(2+), the alkali-metal ions Li(+), Na(+), K(+), and Cs(+), and the transition-metal ion Ag(+). The two neighboring aromatic side chains open the possibility of a novel encapsulation motif of the metal ion in a double cation-π configuration, which is found to be realized for the alkaline-earth complexes and, in a variant form, for the Ag(+) complex. Experimentally, complexes are formed by electrospray ionization, trapped in an FT-ICR mass spectrometer, and characterized by infrared multiple photon dissociation (IRMPD) spectroscopy using the free electron laser FELIX. Interpretation is assisted by thermochemical and IR spectral calculations using density functional theory (DFT). The IRMPD spectrum of protonated PhePhe is reproduced with good fidelity by the calculated spectrum of the most stable conformation, although the additional presence of the secondmost stable conformation is not excluded. All metal-ion complexes have charge-solvated binding modes, with zwitterion (salt bridge) forms being much less stable. The amide oxygen always coordinates to the metal ion, as well as at least one phenyl ring (cation-π interaction). At least one additional chelation site is always occupied, which may be either the amino nitrogen or the carboxy carbonyl oxygen. The alkaline-earth complexes prefer a highly compact caged structure with both phenyl rings providing cation-π stabilization in a "sandwich" configuration (OORR chelation). The alkali-metal complexes prefer open-cage structures with only one cation-π interaction, except perhaps Cs(+). The Ag(+) complex shows a unique preference for the closed-cage amino-bound NORR structure. Ligand-driven perturbations of normal-mode frequencies are generally found to correlate linearly with metal-ion binding energy.

The Curtin-Hammett principle in mass spectrometry

The Curtin-Hammett principle (CHP) is an important concept in physical organic chemistry and is often utilized in the investigation of reaction mechanisms. Two reactants, A and B, in rapid equilibrium, react to form products P(A) and P(B) with rates k(A) and k(B), respectively. If the reaction is under kinetic control and the rate of equilibration between the two reactants is much faster than the reactions to form products, then the branching ratio of products P(A) and P(B) depends solely on the difference in barrier heights for the two product channels. The CHP is based on the fact that the ratio of products formed is not determined by the reactant population ratio. However, the CHP also applies to studies in other areas of chemistry, including mass spectrometry. This Account describes work from our groups in which the results must be interpreted in light of the CHP. These studies illustrate two important implications of the CHP. First, they demonstrate how product distributions cannot be used to assess reactant structure in mechanistic studies in Curtin-Hammett systems. A recent investigation of the structure of hydroxysiliconate anions demonstrated that it was not possible to distinguish between the possible reactant ion structures. A second important implication of the CHP is that the structure of the reactant does not affect the product branching ratio and therefore does not need to be a consideration if the CHP applies. We address this aspect of the discussion through kinetic method studies of the acidities of amino acids and proton affinities of bifunctional compounds. Recently reported mass spectrometric studies illustrate how the CHP puts limitations on what conclusions can be drawn from product distribution studies but also allows experimental methods, such as the kinetic method, to be carried out for complicated systems without having to know all the details of the reactant ion structures. These studies show that although the CHP is most commonly applied in mechanistic studies in physical organic chemistry, it also applies to other areas of chemistry, including mass spectrometry. Although the CHP in some cases limits the conclusions that can be drawn from an experimental study, its proper application can often be used to greatly simplify very complicated chemical systems. Therefore, it is important in mass spectrometry, and indeed, in all areas of chemistry, to recognize those systems in which the CHP should and should not apply.

The sodium ion affinities of cytosine and its methylated derivatives

The sodium ion affinities of cytosine (Cyt), 5-methylcytosine (5MeCyt) and 1-methylcytosine (1MeCyt) have been determined by experimental and quantum chemical methods. Na(+)-bound heterodimers were produced carrying one cytosine or methylated cytosine ligand (designated as C) and one peptide or amino acid reference base (designated as Pep); the Pep molecules included the peptides GlyLeu, GlyPhe, SerGly, and PheGly, and the amino acid His. The dissociation kinetics of these C--Na(+)--Pep ions were determined by collisionally activated dissociation (CAD) and converted to relative and absolute Na(+) affinities via kinetic method approaches. Relative Na(+) affinities increase in the order (kJ/mol): GlyLeu (0) < Cyt (3) < GlyPhe (4) < SerGly (6) < 5MeCyt (8) < PheGly (11) < 1MeCyt (13) < His (17). Anchoring the relative values of the nucleobases to the absolute affinities of the reference bases leads to absolute Na(+) affinities of 214 +/- 8, 219 +/- 8, and 224 +/- 8 kJ/mol for Cyt, 5MeCyt, and 1MeCyt, respectively. Ab initio calculations were used to confirm these results. The computed affinities of Cyt (213 kJ/mol) and 1MeCyt (217 kJ/mol) are in very good agreement with the experiments. These values unambiguously correspond to Na(+) complexes with the keto form of cytosine and its methyl derivatives. Ab initio calculations on tautomerization mechanisms in the gas versus condensed phase are used to discuss why the sodiated keto isomers were formed in the present electrospray ionization (ESI) experiments, but the enol isomers in previous fast atom bombardment (FAB) experiments.

Hydration Energies of Sodiated Amino Acids from Gas-Phase Equilibria Determinations

The sequential hydration of a number of sodiated amino acids is investigated using a high-pressure mass spectrometer. Ions produced continuously by electrospray are injected into the reaction chamber in the pulsed mode where the hydration equilibria, AANa+(H2O)n-1+H2O=AANa+(H2O)n (AA=Val, Pro, Met, Phe, and Gln), and the temperature dependence of the equilibrium constants are measured in the gas phase at 10 mbar (N2 bath gas and known pressure of H2O). The thermochemical properties, DeltaH degrees n, DeltaS degrees n, and DeltaG degrees n, for the hydrated systems are determined and discussed in conjunction with the structural forms. The results show that the binding energies of water to the AANa+ complexes decrease with the increasing number of water molecules. The present results from equilibrium measurements are compared to those from earlier studies obtained by other techniques. A correlation between the free energy changes for the addition of the first and second water molecules to AANa+, and the corresponding sodium ion affinities, is observed. Generally, the hydration free energy becomes weaker as the AA-Na+ bond strength increases.

Effect of metal ions (Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+) and water coordination on the structure of glycine and zwitterionic glycine

Interactions between metal ions and amino acids are common both in solution and in the gas phase. Here, the effect of metal ions and water on the structure of glycine is examined. The effect of metal ions (Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+) and water on structures of Gly.Mn+(H2O)m and GlyZwitt.Mn+(H2O)m (m = 0, 2, 5) complexes have been determined theoretically by employing the hybrid B3LYP exchange-correlation functional and using extended basis sets. Selected calculations were carried out also by means of CBS-QB3 model chemistry. The interaction enthalpies, entropies, and Gibbs energies of eight complexes Gly.Mn+ (Mn+ = Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+) were determined at the B3LYP density functional level of theory. The computed Gibbs energies DeltaG degrees are negative and span a rather broad energy interval (from -90 to -1100 kJ mol(-1)), meaning that the ions studied form strong complexes. The largest interaction Gibbs energy (-1076 kJ mol(-1)) was computed for the NiGly2+ complex. Calculations of the molecular structure and relative stability of the Gly.Mn+(H2O)m and GlyZwitt.Mn+(H2O)m (Mn+ = Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+; m = 0, 2, and 5) systems indicate that in the complexes with monovalent metal cations the most stable species are the NO coordinated metal cations in non-zwitterionic glycine. Divalent cations Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+ prefer coordination via the OO bifurcated bonds of the zwitterionic glycine. Stepwise addition of two and five water molecules leads to considerable changes in the relative stability of the hydrated species. Addition of two water molecules at the metal ion in both Gly.Mn+ and GlyZwitt.Mn+ complexes reduces the relative stability of metallic complexes of glycine. For Mn+ = Li+ or Na+, the addition of five water molecules does not change the relative order of stability. In the Gly.K+ complex, the solvation shell of water molecules around K+ ion has, because of the larger size of the potassium cation, a different structure with a reduced number of hydrogen-bonded contacts. This results in a net preference (by 10.3 kJ mol(-1)) of the GlyZwitt.K+H2O5 system. Addition of five water molecules to the glycine complexes containing divalent cations Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+ results in a net preference for non-zwitterionic glycine species. The computed relative Gibbs energies are quite high (-10 to -38 kJ mol(-1)), and the NO coordination is preferred in the Gly.Mn+(H2O)5 (Mn+ = Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+) complexes over the OO coordination.

The sodium ion affinities of simple di-, tri-, and tetrapeptides

The sodium ion affinities (binding energies) of nineteen peptides containing 2-4 residues have been determined by experimental and computational approaches. Na(+)-bound heterodimers with amino acid and peptide ligands (Pep(1), Pep(2)) were produced by electrospray ionization. The dissociations of these Pep(1)-Na(+)-Pep(2) ions to Pep(1)-Na(+) and Pep(2)-Na(+) were examined by collisionally activated dissociation to construct a ladder of relative affinities via the kinetic method. The accuracy of this ladder was subsequently ascertained by experiments using several excitation energies for four peptide pairs. The relative scale was converted to absolute affinities by anchoring the relative values to the known Na(+) affinity of GlyGly. The Na(+) affinities of AlaAla, HisGly, GlyHis, GlyGlyGly, AlaAlaAla, GlyGlyGlyGly, and AlaAlaAlaAla were also calculated at the MP2(full)/6-311 + G(2d,2p) level of ab initio theory using geometries that were optimized at the MP2(full)/6-31G(d) level for AlaAla or HF/6-31G(d) level for the other peptides; the resulting values agree well with experimental Na(+) affinities. Increasing the peptide size is found to dramatically augment the Na(+) binding energy. The calculations show that in nearly all cases, all available carbonyl oxygens are sodium binding sites in the most stable structures. Whenever side chains are available, as in HisGly and GlyHis, specific additional binding sites are provided to the cation. Oligoglycines and oligoalanines have similar binding modes for the di- and tripeptides, but differ significantly for the tetrapeptides: while the lowest energy structure of GlyGlyGlyGly-Na(+) has the peptide folded around the ion with all four carbonyl oxygens in close contact with Na(+), that of AlaAlaAlaAla-Na(+) involves a pseudo-cyclic peptide in which the C and N termini interact via hydrogen bonding, while Na(+) sits on top of the oxygens of three nearly parallel C=O bonds.

Relative calcium-binding strengths of amino acids determined using the kinetic method

We have measured the relative calcium-binding energies of amino acids using tandem mass spectrometry of Ca(2+)-bound trimeric amino acids. Although calcium-bound dimeric amino acid complexes coordinated too strongly to allow observation of the two competing dissociation products (calcium-bound monomeric ions) required for analysis of their metal binding affinities using the conventional kinetic method, the Ca(2+)-bound trimeric cluster ions dissociated readily to form dimeric cluster ions through simple ligand losses. The calcium-binding energies were obtained by comparing the ratio of the [Ca(2+)(A(1))(2) - H(+)](+) and [Ca(2+) (A(1))(A(2)) - H(+)](+) ions that dissociated from the [Ca(2+) (A(1))(2)(A(2)) - H(+)](+) ion and the ratio of the [Ca(2+)(A(2))(2) - H(+)](+) and [Ca(2+)(A(1)) (A(2)) - H(+)](+) ions that dissociated from the [Ca(2+)(A(1))(A(2))(2) - H(+)](+) ion, where A(1) and A(2) represent two amino acids. The energies deduced from this analysis represent the relative average binding energies of complexes having the form [Ca(2+)(A(1))(2) - H(+)](+). The relative Ca(2+)-binding strengths of the alpha-amino acid complexes follow the order Cys < Ser < Thr < Ile < Leu < Val < Gly < Ala < Pro < Phe < Met < Tyr < Asn < His < Gln < Trp < Lys < Arg. To our knowledge, this report provides the first example of using kinetic methods to determine the relative binding strengths of divalent metal-amino acid complexes.

Proton Affinity of Canavanine and Canaline, Oxyanalogues of Arginine and Ornithine, from the Extended Kinetic Method

The absolute proton affinities of the nonprotein amino acids canavanine and canaline have been determined using the extended kinetic method in an electrospray ionization quadrupole ion trap instrument. Canavanine results from the substitution of an oxygen atom for the delta-CH2 group in the side chain of the protein amino acid arginine, whereas canaline results from a similar substitution at the delta-CH2 group in the side chain of ornithine. Absolute proton affinities of 1001+/-9 and 950+/-7 kJ/mol are obtained for canavanine and canaline, respectively. For canaline, this proton affinity is in excellent agreement with theoretical predictions obtained using the hybrid density functional theory method B3LYP/6-311++G**//B3LYP/6-31+G*. For canavanine, theory predicts a somewhat larger proton affinity of 1015 kJ/mol. Oxygen atom substitution in these nonprotein amino acids results in a decrease in their proton affinities of 40-50 kJ/mol compared to arginine and ornithine.

Structures of Lithiated Lysine and Structural Analogues in the Gas Phase: Effects of Water and Proton Affinity on Zwitterionic Stability

The structures of lithiated lysine, ornithine, and related molecules, both with and without a water molecule, are investigated using both density functional theory and blackbody infrared radiative dissociation experiments. The lowest-energy structure of lithiated lysine without a water molecule is nonzwitterionic; the metal ion interacts with both nitrogen atoms and the carbonyl oxygen. Structures in which lysine is zwitterionic are higher in energy by more than 29 kJ/mol. In contrast, the singly hydrated clusters with the zwitterionic and nonzwitterionic forms of lysine are more similar in energy, with the nonzwitterionic form more stable by only approximately 7 kJ/mol. Thus, a single water molecule can substantially stabilize the zwitterionic form of an amino acid. Analogous molecules that have methyl groups attached to either the N-terminus (NMeLys) or the side-chain amine (Lys(Me)) have proton affinities greater than that of lysine. In the lithiated clusters with a water molecule attached, the zwitterionic forms of NMeLys and Lys(Me) are calculated to be approximately 4 and approximately 11 kJ/mol more stable than the nonzwitterionic forms, respectively. Calculations of the potential-energy pathway for interconversion between the different forms of lysine in the lithiated complex indicate multiple stable intermediates with an overall barrier height of approximately 83 kJ/mol between the lowest-energy nonzwitterionic form and the most accessible zwitterionic form. Experimentally determined binding energies of water are similar for all these complexes and range from 57 to 64 kJ/mol. These results suggest that loss of a water molecule from the lysine complexes is both energetically and entropically favored compared to interconversion between the nonzwitterionic and zwitterionic structures. Comparisons to calculated binding energies of water to the various structures show that the experimental results are most consistent with the nonzwitterionic forms.

Extending, and Repositioning, a Thermochemical Ladder: High-Level Quantum Chemical Calculations on the Sodium Cation Affinity Scale

High-level ab initio quantum chemical calculations, at the CP-dG2thaw level of theory, are reported for coordination of Na+ to a wide assortment of small organic and inorganic ligands. The ligands range in size from H to C6H6, and include 22 of the ligands for which precise relative sodium ion binding free energies have been determined by recent Fourier transform ion cyclotron resonance and guided ion beam studies. Agreement with the relative experimental values is excellent (+/-1.1 kJ mol(-1)), and agreement with the absolute scale (obtained when these relative values are pegged to the CH3NH2 "anchor" value measured in a high-pressure mass spectrometric study) is only marginally poorer, with CP-dG2thaw values exceeding the absolute experimental DeltaG(298) values by an average of 2.1 kJ mol(-1). The excellent agreement between experiment and the CP-dG2thaw technique also suggests that the additional 97 ligands surveyed here (which, in many cases, are not readily susceptible to laboratory investigation) can also be reliably fitted to the existing experimental scale. However, while CP-dG2thaw and the experimental ladder are in close accord, a small set of higher level ab initio calculations on sodium ion/ligand complexes (including several values obtained here using the W1 protocol) suggests that the CP-dG2thaw values are themselves too low by approximately 2.5 kJ mol(-1), thereby implying that the accepted laboratory values are typically 4.6 kJ mol(-1) too low. The present work also highlights the importance of Na+/ligand binding energy determinations (whether by experimental or theoretical approaches) on a case-by-case basis: trends in increasing binding energy along homologous series of compounds are not reliably predictable, nor are binding site preferences or chelating tendencies in polyfunctional compounds.

Carbonate, carbamate, urea, and guanidine as model species for functional groups in biological molecules — A combined density functional theory and mass spectrometry examination of polysodiation and gas-phase dissociation

Electrospray mass spectrometric analyses of carbonate, carbamate, urea, and guanidine in the presence of sodium ions results in polysodiated complexes of the type, [M – (n – 1)H + nNa]+, where M is the species analyzed and n = 2 and (or) 3. These complexes are unusual in their high sodium contents and their attractive structures. The [M – H + 2Na]+ complexes fragment to produce [Na2NH2]+, Na+, and a product ion resulting from the loss of NH3. The [M – 2H + 3Na]+ complexes dissociate to give a number of fragment ions, including [Na3O]+, [Na2OH]+, [NaOCNNa]+, [Na2NH2]+, [Na3CN2]+, and Na+, depending on the elemental composition of M. These dissociation reactions were examined both experimentally by mass spectrometry and theoretically by density functional theory. Reaction profiles as calculated at B3LYP/6-311G** for disodiated carabamate, urea, and guanidine, and for all the bases in trisodiated form are given. Both experimental and theory show [Na3CO3]+ to have high stability towards collisional activation and to preferentially dissociate to give [Na3O]+ and CO2, rather than Na+ and Na2CO3.Key words: carbonate, carbamate, urea, guanidine, electrospray, collision-induced dissociation, density functional theory, sodium ion binding energies.

Structures and Hydration Enthalpies of Cationized Glutamine and Structural Analogues in the Gas Phase

The structures of lithiated and sodiated glutamine, both with and without a water molecule, are investigated using experiment and theory. Loss of water from these complexes and from lithiated and sodiated complexes of asparagine methyl ester, asparagine ethyl ester, and glutamine methyl ester is probed with blackbody infrared radiative dissociation experiments performed over a wide temperature range. Threshold dissociation energies, E(o), for loss of a water molecule from these complexes are obtained from master equation modeling of these data. The values of E(o) are 63 +/- 1 and 53 +/- 1 kJ/mol for the lithiated and sodiated glutamine complexes, respectively. These values are similar to those for the nonzwitterionic model complexes and are in excellent agreement with calculated values. In contrast, water binding to the zwitterionic form is calculated to be significantly higher. These results indicate that glutamine in these lithiated and sodiated complexes with a water molecule are nonzwitterionic. Complexes with the asparagine side chain have slightly higher E(o) values than those with the glutamine side chain, a result consistent with more effective solvation of the metal ion due to the slightly longer side chain of glutamine. Calculations indicate that lithiated and sodiated glutamine are nonzwitterionic, with the metal ion interacting with the amine nitrogen and carbonyl oxygen from the amino acid backbone and the amide oxygen of the side chain. Addition of a water molecule does not affect the lowest-energy structure of lithiated glutamine, whereas, for sodiated glutamine, the lowest-energy zwitterionic and nonzwitterionic structures are essentially isoenergetic.

Infrared Fingerprint Spectroscopy and Theoretical Studies of Potassium Ion Tagged Amino Acids and Peptides in the Gas Phase

Infrared multiple-photon dissociation spectroscopy is effected on the K(+) tagged aromatic amino acids tyrosine and phenylalanine, as well as the K(+) tagged peptides bradykinin fragment 1-5 and [Leu]-enkephalin. The fingerprint (800-1800 cm(-1)) infrared spectra of these species are compared to density-functional theory (DFT) calculated spectra to determine whether the complex is in the charge solvation (CS) or salt bridge (SB) (i.e. zwitterionic) configuration. For the aromatic amino acids the CS structure is favored and the tridentate N/O/ring structure is found to be the preferred binding geometry for K(+). The experimental and theoretical evidence for bradykinin fragment 1-5 tagged with K(+) suggests that the SB structure is favored; the calculations indicate a head-to-tail looped structure stabilized by a salt bridge between the protonated guanidine group and the deprotonated C-terminus, which allows K(+) to sit in a binding pocket with five C=O electrostatic interactions. For K(+) tagged [Leu]-enkephalin the spectroscopic evidence is not as clear. While the calculations clearly favor a CS structure and the observation of a weak carboxylic acid C=O stretching band in the infrared spectrum matches this finding, the prominence of a band at 1600 cm(-1) renders the analysis more ambiguous, and hence the presence of some salt bridge ions cannot be excluded. Another striking feature in the [Leu]-enkephalin spectrum is the high infrared activity of the tyrosine side-chain modes, which can be clearly identified from comparison to the [Tyr + K](+) experimental spectrum, but which is not reproduced by the DFT calculations.

Effects of charge state and cationizing agent on the electron capture dissociation of a peptide

Electron capture dissociation (ECD) is a promising method for de novo sequencing proteins and peptides and for locating the positions of labile posttranslational modifications and binding sites of noncovalently bound species. We report the ECD of a synthetic peptide containing 10 alanine residues and 6 lysine residues uniformly distributed across the sequence. ECD of the (M + 2H)(2+) produces a limited range of c (c(7)-c(15)) and z (z(9)-z(15)) fragment ions, but ECD of higher charge states produces a wider range of c (c(2)-c(15)) and z (z(2)-z(6), z(9)-z(15)) ions. Fragmentation efficiency increases with increasing precursor charge state, and efficiencies up to 88% are achieved. Heating the (M + 2H)(2+) to 150 degrees C does not increase the observed range of ECD fragment ions, indicating that the limited products are due to backbone cleavages occurring near charges and not due to effects of tertiary structure. ECD of the (M + 2Li)(2+) and (M + 2Cs)(2+) produces di- and monometalated analogues of the same c and z ions observed from the (M + 2H)(2+), with the abundance of dimetalated fragment ions increasing with fragment ion mass, a result consistent with the metal cations being located near the peptide termini to minimize Coulombic repulsion. In stark contrast to the ECD results, collisional activation of cesiated dications overwhelmingly results in ejection of Cs(+). The abundance of cesiated fragment ions formed from ECD of the (M + Cs + Li)(2+) exceeds that of lithiated fragment ions by 10:1. ECD of the (M + H + Li)(2+) results in exclusively lithiated c and z ions, indicating an overwhelming preference for neutralization and cleavage at protonated sites over metalated sites. These results are consistent with preferential neutralization of the cation with the highest recombination energy.

Experimental validation of Gaussian-3 lithium cation affinities of amides: implications for the gas-phase lithium cation basicity scale

Using a refined Gaussian-3 (G3) protocol, the highest level of ab initio calculations reported so far, we have established the Li+ cation binding enthalpy (affinity) at 0 K (in kJ mol-1) for formamide (195.7), N-methylformamide (209.2), N,N'-dimethylformamide (220.0), acetamide (211.7), N-methylacetamide (222.5), and N,N'-dimethylacetamide (230.1), with an estimated maximum uncertainty of +/-8 kJ mol-1. With these six theoretical lithium cation binding affinities as reference values, the absolute Li+ affinities of imidazole and dimethoxyethane were determined by the extended kinetic method, and by adopting the statistical data treatment protocol recently proposed by Armentrout. The Li+ affinities obtained for these two ligands are in good agreement (within 6 kJ mol-1) with recent values determined by the threshold collision-induced dissociation method, and consistent with the Li+ basicity values first reported by Taft and co-workers in 1990. Our study confirms that the previously suggested, and recently implemented, downward revision of Taft's original basicity scale by 10.9 kJ mol-1 is justified for ligands with revised basicities less than 151 kJ mol-1. However, for selected ligands with Li+ basicities greater than 151 kJ mol-1, including some of the six amides studied in this work, the reported discrepancy between theoretical and experimental estimates in the revised Li+ basicity scale of Burk et al. is likely to arise from experimental uncertainties.

Potential energy surfaces for the gas-phase interaction between alpha-alanine and alkali metal Ions (Li+, Na+, K+). A density functional study

The conformations of alpha-alanine and those of the eight most stable adducts that it forms with lithium, sodium, and potassium cations were studied by means of the density functional theory using the hybrid B3LYP exchange correlation potential and the 6-311++G basis set. Minima and transition states characterizing the energetic paths for the interaction of the three metal ions with the free amino acid were explored in detail, and the results show that they are almost the same in the cases of lithium and sodium ions. At the absolute minimum, the Li+ and Na+ cations appear to be contemporaneously linked to the carbonyl oxygen and nitrogen atoms of alpha-alanine. Slight differences were found for potassium; the most stable adduct has the cation coordinated to both oxygen atoms of carboxyl function. The influence of the low-energy conformers in the determination of the gas-phase absolute affinities is demonstrated.

Gas phase hydrogen deuterium exchange reactions of a model peptide: FT-ICR and computational analyses of metal induced conformational mutations

We utilized gas phase hydrogen/deuterium (H/D) exchange reactions and ab initio calculations to investigate the complexation between a model peptide (Arg-Gly-Asp[triple bond]RGD) with various alkali metal ions. The peptide conformation is drastically altered upon alkali metal ion complexation. The associated conformational changes depend on both the number and type of complexing alkali metal ions. Sodium has a smaller ionic diameter and prefers a multidentate interaction that involves all three amino acids of the peptide. Conversely, potassium and cesium form different types of complexes with the RGD. The [RGD + 2Cs - H]+ species exhibit the slowest H/D exchange reactivity (reaction rate constant of approximately 6 x 10(-13) cm3molecule(-1)s(-1) for the fastest exchanging labile hydrogen with ND3). The reaction rate constant of the protonated RGD is two orders of magnitude faster than that of the [RGD + 2Cs - H]+. Addition of the first cesium to the RGD reduces the H/D exchange reaction rate constant (i.e., D0) by a factor of seven whereas sodium reduces this value by a factor of thirty. Conversely, addition of the second alkali metal ions has the opposite effect; the rate of D0 disappearance for all [RGD + 2Met - H]+ species (Met[triple bond]Na, K, and Cs) decreases with the alkali metal ion size.

Cation-pi interactions and the gas-phase thermochemistry of the Na(+)/phenylalanine complex

The complex of Na(+) with phenylalanine (Phe) is a prototype for the participation of cation-pi interactions in metal-ion binding to biological molecules. A recent comparison of this complex with the Na(+)/alanine (Na(+)/Ala) counterpart suggested only a small contribution of the phenyl ring interaction to binding, casting doubt on the extent of the cation-pi effect. The present work reexamines this thermochemistry using ligand-exchange equilibrium measurements in the Fourier transform ion cyclotron resonance (FT-ICR) ion trapping mass spectrometer. An increment of 7 +/- 2 kcal mol(-1) was found in the Ala/Phe comparison of binding enthalpies, confirming the importance of cation-pi binding enhancement in the Phe case. Absolute Na(+) binding enthalpies of 38 +/- 2 and 45 +/- 2 kcal mol(-1) were assigned for Ala and Phe, respectively, using pyridine as the thermochemical reference ligand. All of these results were supported by quantum calculations using both density functional and Hartree-Fock/MP2 methods, improved in several respects over previous calculations. Alanine methyl ester (AlaMe) was also observed, and found to have an Na(+) ion affinity larger by 2.3 kcal mol(-1) than Ala. New, lower energy conformations of neutral Phe were discovered in the computations.

Cation-pi effects in the complexation of Na+ and K+ with Phe, Tyr, and Trp in the gas phase

Na+ and K+ gas-phase affinities of the three aromatic amino acids Phe, Tyr, and Trp were measured by the kinetic method. Na+ binds these amino acids much more strongly than K+, and for both metal ions the binding strength was found to follow the order Phe < or = Tyr < Trp. Quantum chemical calculations by density functional theory (DFT) gave the same qualitative ordering, but suggested a somewhat larger Phe/Trp increment. These results are in acceptable agreement with predictions based on the binding of Na+ and K+ to the side chain model molecules benzene, phenol, and indole, and are also in reasonable agreement with the predictions from purely electrostatic calculations of the side-chain binding effects. The binding energies were compared with those to the aliphatic amino acids glycine and alanine. Binding to the aromatic amino acids was found to be stronger both experimentally and computationally, but the DFT calculations indicate substantially larger increments relative to alanine than shown by the experiments. Possible reasons for this difference are discussed. The metal ion binding energies show the same trends as the proton affinities.