Structural and energetic effects in the molecular recognition of protonated peptidomimetic bases by 18-crown-6

Absolute 18-crown-6 (18C6) affinities of nine protonated peptidomimetic bases are determined using guided ion beam tandem mass spectrometry techniques. The bases (B) included in this work are mimics for the n-terminal amino group and the side chains of the basic amino acids, i.e., the favorable sites for binding of 18C6 to peptides and proteins. Isopropylamine is chosen as a mimic for the n-terminal amino group, imidazole and 4-methylimidazole are chosen as mimics for the side chain of histidine (His), 1-methylguanidine is chosen as a mimic for the side chain of arginine (Arg), and several primary amines including methylamine, ethylamine, n-propylamine, n-butylamine, and 1,5-diamino pentane as mimics for the side chain of lysine (Lys). Theoretical electronic structure calculations are performed to determine stable geometries and energetics for neutral and protonated 18C6 and the peptidomimetic bases, as well as the proton bound complexes comprised of these species, (B)H(+)(18C6). The measured 18C6 binding affinities of the Lys side chain mimics are larger than the measured binding affinities of the mimics for Arg and His. These results suggest that the Lys side chains should be the preferred binding sites for 18C6 complexation to peptides and proteins. Present results also suggest that competition between Arg or His and Lys for 18C6 is not significant. The mimic for the n-terminal amino group exhibits a measured binding affinity for 18C6 that is similar to or greater than that of the Lys side chain mimics. However, theory suggests that binding to n-terminal amino group mimic is weaker than that to all of the Lys mimics. These results suggest that the n-terminal amino group may compete with the Lys side chains for 18C6 complexation.     

Metal cation dependence of interactions with amino acids: bond energies of Cs+ to Gly, Pro, Ser, Thr, and Cys

The interactions of cesium cations with five amino acids (AA) including glycine (Gly), proline (Pro), serine (Ser), threonine (Thr), and cysteine (Cys) are examined in detail. Experimentally, the bond dissociation energies (BDEs) are determined using threshold collision-induced dissociation of the Cs(+)(AA) complexes with xenon in a guided ion beam tandem mass spectrometer. Analyses of the energy-dependent cross sections include consideration of unimolecular decay rates, internal energy of the reactant ions, and multiple ion-neutral collisions. Bond dissociation energies (0 K) of 93.3 ± 2.5, 107.9 ± 4.6, 102.3 ± 4.1, 105.4 ± 4.3, and 96.8 ± 4.2 kJ/mol are determined for complexes of Cs(+) with Gly, Pro, Ser, Thr, and Cys, respectively. Quantum chemical calculations are conducted at the B3LYP, B3P86, MP2(full), and M06 levels of theory with geometries and zero-point energies calculated at the B3LYP level using both HW*/6-311+G(2d,2p) and def2-TZVPPD basis sets. Results obtained using the former basis sets are systematically low compared to the experimental bond energies, whereas the latter basis sets show good agreement. For Cs(+)(Gly), theory predicts the ground-state conformer has the cesium cation binding to the carbonyl group of the carboxylic acid. For Cs(+)(Pro), the secondary nitrogen accepts the carboxylic acid hydrogen to form the zwitterionic structure, and the metal cation binds to both oxygens. Cs(+)(Ser), Cs(+)(Thr), and Cs(+)(Cys) are found to have tridentate binding at the MP2(full) level, whereas the density functional approaches slightly prefer bidentate binding of Cs(+) at the carboxylic acid moiety. Comparison of these results to those for the smaller alkali cations provides insight into the trends in binding affinities and structures associated with metal cation variations.     

Structural and Energetic Effects in the Molecular Recognition of Acetylated Amino Acids by 18-Crown-6

Absolute 18-crown-6 (18C6) binding affinities of four protonated acetylated amino acids (AcAAs) are determined using guided ion beam tandem mass spectrometry techniques. The AcAAs examined in this work include: N-terminal acetylated lysine (N_???CAcLys), histidine (N_???CAcHis), and arginine (N_???CAcArg) as well as side chain acetylated lysine (N_???CAcLys). The kinetic-energy-dependent cross sections for collision-induced dissociation (CID) of the (AcAA)H⁺(18C6) complexes are analyzed using an empirical threshold law to extract absolute 0 and 298?K (AcAA)H⁺?18C6 bond dissociation energies (BDEs) after accounting for the effects of multiple collisions, kinetic and internal energy distributions of the reactants, and unimolecular dissociation lifetimes. Theoretical electronic structure calculations are performed to determine stable geometries and energetics for neutral and protonated 18C6 and the AcAAs as well as the proton bound complexes of these species, (AcAA)H⁺(18C6), at the B3LYP/6-311+G(2d,2p)//B3LYP/6-31?G* and M06/6-311+G(2d,2p)//B3LYP/6-31G* levels of theory. For all four (AcAA)H⁺(18C6) complexes, loss of neutral 18C6 corresponds to the most favorable dissociation pathway. At elevated energies, products arising from sequential dissociation of the primary CID product, H⁺(AcAA), are also observed. Protonated N_???CAcLys exhibits a greater 18C6 binding affinity than other protonated N_???CAcAAs, suggesting that the side chains of Lys residues are the preferred binding sites for 18C6 complexation to peptides and proteins. N_???CAcLys exhibits a greater 18C6 binding affinity than N_???CAcLys, suggesting that binding of 18C6 to the side chain of Lys residues is more favorable than to the N-terminal amino group of Lys.     

Tautomerization in the formation and collision-induced dissociation of alkali metal cation-cytosine complexes

Noncovalent interactions between alkali metal cations and the various low-energy tautomeric forms of cytosine are investigated both experimentally and theoretically. Threshold collision-induced dissociation (CID) of M(+)(cytosine) complexes with Xe is studied using guided ion beam tandem mass spectrometry, where M(+) = Li(+), Na(+), and K(+). In all cases, the only dissociation pathway observed corresponds to endothermic loss of the intact cytosine molecule. The cross-section thresholds are interpreted to yield 0 and 298 K bond dissociation energies (BDEs) for the M(+)(cytosine) complexes after accounting for the effects of multiple ion-neutral collisions, the kinetic and internal energy distributions of the reactants, and dissociation lifetimes. Ab initio calculations are performed at the MP2(full)/6-31G* level of theory to determine the structures of the neutral cytosine tautomers, the M(+)(cytosine) complexes, and the TSs for unimolecular tautomerization. The molecular parameters derived from these structures are employed for the calculation of the unimolecular rates for tautomerization and the thermochemical analysis of the experimental data. Theoretical BDEs of the various M(+)(cytosine) complexes and the energy barriers for the unimolecular tautomerization of these complexes are determined at MP2(full)/6-311+G(2d,2p) level of theory using the MP2(full)/6-31G* optimized geometries. In addition, BDEs for the Li(+)(cytosine) complexes are also determined at the G3 level of theory. Based upon the tautomeric mixture generated upon thermal vaporization of cytosine, calculated M(+)-cytosine BDEs and barriers to tautomerization for the low-energy tautomeric forms of M(+)(cytosine), and measured thresholds for CID of M(+)(cytosine) complexes, we conclude that tautomerization occurs during both complex formation and CID.     

Noncovalent interactions of Zn+ with N-donor ligands (pyridine, 4,4'-dipyridyl, 2,2'-dipyridyl, and 1,10-phenanthroline): collision-induced dissociation and theoretical studies

The binding interactions in complexes of Zn(+) with nitrogen donor ligands, (N-L) = pyridine (x = 1-4), 4,4'-dipyridyl (x = 1-3), 2,2'-dipyridyl (x = 1-2), and 1,10-phenanthroline (x = 1-2), are examined in detail. The bond dissociation energies (BDEs) for loss of an intact ligand from the Zn(+)(N-L)(x) complexes are reported. Experimental BDEs are obtained from thermochemical analyses of the threshold regions of the collision-induced dissociation cross sections of Zn(+)(N-L)(x) complexes. Density functional theory calculations at the B3LYP/6-31G* level of theory are performed to determine stable structures of these species and to provide molecular parameters needed for the thermochemical analysis of experimental data. Relative stabilities of the various conformations of these N-donor ligands and their complexes to Zn(+) as well as theoretical BDEs are determined from single point energy calculations at the B3LYP/6-311+G(2d,2p) and M06/6-311+G(2d,2p) levels of theory using the B3LYP/6-31G* optimized geometries. The experimental BDEs for the Zn(+)(N-L)(x) complexes are in reasonably good agreement with values derived from density functional theory calculations. BDEs derived from M06 calculations provide better agreement with the measured values than those based on B3LYP calculations. Trends in the sequential BDEs are explained in terms of sp polarization of Zn(+) and repulsive ligand-ligand interactions. Comparisons are made to the analogous Cu(+)(N-L)(x) and Ni(+)(N-L)(x) complexes previously studied.     

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.     

Hydration of potassiated amino acids in the gas phase

The thermochemistry of stepwise hydration of several potassiated amino acids was studied by measuring the gas-phase equilibria, AAK(+)(H(2)O)(n-1) + H(2)O = AAK(+)(H(2)O)(n) (AA = Gly, AL, Val, Met, Pro, and Phe), using a high-pressure mass spectrometer. The AAK(+) ions were obtained by electrospray and the equilibrium constants K(n-1,n) were measured in a pulsed reaction chamber at 10 mbar bath gas, N(2), containing a known partial pressure of water vapor. Determination of the equilibrium constants at different temperatures was used to obtain the DeltaH(n)(o), DeltaS(n)(o), and DeltaG(n)(o) values. The results indicate that the water binding energy in AAK(+)(H(2)O) decreases as the K(+) affinity to AA increases. This trend in binding energies is explained in terms of changes in the side-chain substituent, which delocalize the positive charge from K(+) to AA in AAK(+) complexes, varying the AAK(+)-H(2)O electrostatic interaction.     

Cation-pi interactions with a model for the side chain of tryptophan: structures and absolute binding energies of alkali metal cation-indole complexes

Threshold collision-induced dissociation techniques are employed to determine bond dissociation energies (BDEs) of mono- and bis-complexes of alkali metal cations, Li+, Na+, K+, Rb+, and Cs+, with indole, C8H7N. The primary and lowest energy dissociation pathway in all cases is endothermic loss of an intact indole ligand. Sequential loss of a second indole ligand is observed at elevated energies for the bis-complexes. Density functional theory calculations at the B3LYP/6-31G level of theory are used to determine the structures, vibrational frequencies, and rotational constants of these complexes. Theoretical BDEs are determined from single point energy calculations at the MP2(full)/6-311+G(2d,2p) level using the B3LYP/6-31G* geometries. The agreement between theory and experiment is very good for all complexes except Li+ (C8H7N), where theory underestimates the strength of the binding. The trends in the BDEs of these alkali metal cation-indole complexes are compared with the analogous benzene and naphthalene complexes to examine the influence of the extended pi network and heteroatom on the strength of the cation-pi interaction. The Na+ and K+ binding affinities of benzene, phenol, and indole are also compared to those of the aromatic amino acids, phenylalanine, tyrosine, and tryptophan to elucidate the factors that contribute to the binding in complexes to the aromatic amino acids. The nature of the binding and trends in the BDEs of cation-pi complexes between alkali metal cations and benzene, phenol, and indole are examined to help understand nature's preference for engaging tryptophan over phenylalanine and tyrosine in cation-pi interactions in biological systems.     

Binding Motifs in the Naked Complexes of Target Amino Acids with an Excerpt of Antitumor Active Biomolecule: An Ion Vibrational Spectroscopy Assay

The structures of proton-bound complexes of 5,7-dimethoxy-4H-chromen-4-one ( 1 ) and basic amino acids (AAs), namely, histidine (His) and lysine (Lys), have been examined by means of mass spectrometry coupled with IR ion spectroscopy and quantum chemical calculations. This selection of systems is based on the fact that 1 represents a portion of glabrescione B, a natural small molecule of promising antitumor activity, while His and Lys are protein residues lining the cavity of the alleged receptor binding site. These species are thus a model of the bioactive adduct, although clearly the isolated state of the present study bears little resemblance to the complex biological environment. A common feature of [ 1 +AA+H]⁺ complexes is the presence of a protonated AA bound to neutral 1 , in spite of the fact that the gas-phase basicity of 1 is comparable to those of Lys and His. The carbonyl group of 1 acts as a powerful hydrogen-bond acceptor. Within [ 1 +AA+H]⁺ the side-chain substituents (imidazole group for His and terminal amino group for Lys) present comparable basic properties to those of the α-amino group, taking part to a cooperative hydrogen-bond network. Structural assignment, relying on the comparative analysis of the infrared multiple photon dissociation (IRMPD) spectrum and calculated IR spectra for the candidate geometries, derives from an examination over two frequency ranges: 900–1800 and 2900–3700 cm⁻¹. Information gained from the latter one proved especially valuable, for example, pointing to the contribution of species characterized by an unperturbed carboxylic OH or imidazole NH stretching mode.     

Putative One‐Pot Prebiotic Polypeptides with Ribonucleolytic Activity

KIA7, a peptide with a highly restricted set of amino acids (Lys, Ile, Ala, Gly and Tyr), adopts a specifically folded structure. Some amino acids, including Lys, Ile, Ala, Gly and His, form under the same putative prebiotic conditions, whereas different conditions are needed for producing Tyr, Phe and Trp. Herein, we report the 3D structure and conformational stability of the peptide KIA7H, which is composed of only Lys, Ile, Ala, Gly and His. When the imidazole group is neutral, this 20‐mer peptide adopts a four‐helix bundle with a specifically packed hydrophobic core. Therefore, one‐pot prebiotic proteins with well‐defined structures might have arisen early in chemical evolution. The Trp variant, KIA7W, was also studied. It adopts a 3D structure similar to that of KIA7H and its previously studied Tyr and Phe variants, but is remarkably more stable. When tested for ribonucleolytic activity, KIA7H, KIA7W and even short, unstructured peptides rich in His and Lys, in combination with Mg⁺⁺, Mn⁺⁺ or Ni⁺⁺ (but not Cu⁺⁺, Zn⁺⁺ or EDTA) specifically cleave the single‐stranded region in an RNA stem–loop. This suggests that prebiotic peptide–divalent cation complexes with ribonucleolytic activity might have co‐inhabited the RNA world.     

Experimental and theoretical studies of potassium cation interactions with the acidic amino acids and their amide derivatives

The binding of K(+) to aspartic acid (Asp), glutamic acid (Glu), asparagine (Asn), and glutamine (Gln) is examined in detail by studying the collision-induced dissociation (CID) of the four potassium cation-bound amino acid complexes with Xe using a guided ion beam tandem mass spectrometer (GIBMS). Formed by electrospray ionization, these complexes have energy-dependent CID cross sections that are analyzed to provide 0 K bond energies after accounting for unimolecular decay rates, internal energy of reactant ions, and multiple ion-molecule collisions. Quantum chemical calculations for a number of geometric conformations of each K(+)(L) complex are determined at the B3LYP/6-311+G(d,p) level with single-point energies calculated at B3LYP, B3P86, and MP2(full) levels using a 6-311+G(2d,2p) basis set. Theoretical bond dissociation energies are in good agreement with the experimental values. This coordinated examination of both experimental work and quantum chemical calculations allows for a comprehensive understanding of the molecular interactions of K(+) with the Asx and Glx amino acids. K(+) binding affinities for the amide complexes are systematically stronger than those for the acid complexes by 9+/-1 kJ/mol, which is attributed to an inductive effect of the OH group in the carboxylic acid side chain. Additionally, the K(+) binding affinity for the longer-chain amino acids (Glx) is enhanced by 5+/-1 kJ/mol compared to the shorter-chain Asx because steric effects are reduced. Further, a detailed comparison between experimental and theoretical results reveals interesting differences in the binding of K(+) and Na(+) to these amino acids.     

Reaction of Cu+ with dimethoxyethane: competition between association and multiple dissociation channels

The reaction of Cu+ with dimethoxyethane (DXE) is studied using kinetic-energy dependent guided ion beam mass spectrometry. The bimolecular reaction forms an associative Cu(+)(DXE) complex that is long-lived and dissociates into several competitive channels: C4H9O2(+)+CuH, Cu(+)(C3H6O)+CH3OH, back to reactants, and other minor channels. The kinetic-energy dependences of the cross sections for the three largest product channels are interpreted with several different models (including rigorous phase space theory) to yield 0 K bond energies after accounting for the effects of multiple ion-molecule collisions, internal energy of the reactant ions, Doppler broadening, and dissociation lifetimes. These values are compared with bond energies obtained from collision-induced dissociation (CID) studies of the Cu(+)(DXE) complex and found to be self-consistent. Although all models provide reasonable thermochemistry, phase space theory reproduces the details of the cross sections most accurately. We also examine the dynamics of this reaction using time-of-flight methods and a retarding potential analysis. This provides additional insight into the unimolecular decay of the long-lived Cu(+)(DXE) association complex. Comparison of results from this study with those from the complementary CID study, thus forming the same energized Cu(+)(DXE) complex in two distinct ways, allows an assessment of the models used to interpret CID thresholds.     

Experimental and theoretical studies of sodium cation interactions with D-arabinose, xylose, glucose, and galactose

The binding of Na (+) to arabinose (Ara), xylose (Xyl), glucose (Glc), and galactose (Gal) is examined in detail by studying the collision-induced dissociation (CID) of the four sodiated monosaccharide complexes with Xe using a guided ion beam tandem mass spectrometer (GIBMS). Analysis of the energy-dependent CID cross-sections provides 0 K sodium cation affinities for experimental complexes after accounting for unimolecular decay rates, internal energy of reactant ions, and multiple ion-neutral collisions. Quantum chemical calculations for a number of geometric conformations of each Na (+)(L) complex with a comprehensive analysis of the alpha and beta anomeric forms are determined at the B3LYP/6-311+G(d,p) level with single-point energies calculated at MP2(full), B3LYP, and B3P86 levels using a 6-311+G(2d,2p) basis set. This coordinated examination of both experimental work and quantum chemical calculations allows for determination of the bond energy for both the alpha and beta forms of each monosaccharide studied here. An understanding of the energetic contributions of individual structural characteristics as well as the energetic trends in binding among the monosaccharides is developed. Structural characteristics that affect the energetics of binding involve multidentate sodium cation coordination, ring sterics, and hydrogen bonding schemes. The overall trend in sodium binding affinities for the eight ligands follows beta-Ara < alpha-Ara < beta-Xyl < beta-Glc < alpha-Glc < alpha;-Xyl < alpha-Gal < beta-Gal.     

Amino acids of ricin and its polypeptides

Ricin and its corresponding polypeptides (A & B chain) were purified from castor seed. The molecular weight of ricin subunits were 29,000 and 28,000 daltons. The amino acids in ricin determined were Asp45 The22 Ser40 Glu53 Cys4 Gly96 His5 Ile21 Leu33 Lys20 Met4 Phe13 Pro37 Tyr11 Ala45 Val23 Arg20 indicating that ricin contains approximately 516 amino acid residues. The amino acids of the two subunits of ricin A and B chains were Asp23 The12 Ser21 Glu29 Cys2 Gly48 His3 Ile12, Leu17 Lys10 Met2 Phe6 Pro17 Tyr7 Ala35 Val13 Arg13 while in B chain the amino acids were Asp22 The10 Ser19 Glu25 Cys2 Gly47 His1 Ile10, Leu15 Lys11 Met1 Phe7 Pro6 Tyr5 Ala32Val11 Arg10. The total helical content of ricin came around 53.6% which is a new observation.     

Hydration energies of deprotonated amino acids from gas phase equilibria measurements

Singly hydrated clusters of deprotonated amino acids were studied using an electrospray high-pressure mass spectrometer equipped with a pulsed ion-beam reaction chamber. Thermochemical data, DeltaH(o), DeltaS(o), and DeltaG(o), for the hydration reaction [AA - H](-) + H(2)O = [AA - H](-).(H(2)O) were obtained from gas-phase equilibria determinations for AA = Gly, Ala, Val, Pro, Phe, Lys, Met, Trp, Gln, Arg, and Asp. The hydration free-energy changes are found to depend significantly on the side-chain substituents. The water binding energy in [AA - H](-).(H(2)O) increases with the gas-phase acidity of AA. The anionic hydrogen bond strengths in [AA - H](-).(H(2)O) are compared with those of the cationic bonds in the corresponding AAH(+).(H(2)O) systems.     

Revising the proton affinity scale of the naturally occurring α-amino acids

The proton affinities (PA) of the 20 naturally occurring alpha-amino acids (AA) have been determined computationally by means of density functional theory (DFT) and high-level G2(MP2) calculations. These theoretical PAs, together with data that have appeared since 1997 in the literature, are used to validate the most reasonable currently available PA scale for AAs (Harrison, A. G. Mass Spectrom. Rev. 1997, 16, 201-217.). Significant scatter is observed for the PAs of Ser, Asp, Phe, Asn, Met, Pro, Gln, Glu, Trp, His, Lys, and Arg, many of which have a basic side-chain functionality. Critical review of the available data leads to new consensus PAs for Asn, Gln, Met, and Arg of 222.4, 230.5, 223.7, and 250.2 kcal/mol, respectively.     

Gas‐phase doubly charged complexes of cyclic peptides with copper in +1, +2 and +3 formal oxidation states: formation,structures and electron capture dissociation

Copper complexes with a cyclic D‐His‐β‐Ala‐L‐His‐L‐Lys and all‐L‐His‐β‐Ala‐His‐Lys peptides were generated by electrospray which were doubly charged ions that had different formal oxidation states of Cu(I), Cu(II) and Cu(III) and different protonation states of the peptide ligands. Electron capture dissociation showed no substantial differences between the D‐His and L‐His complexes. All complexes underwent peptide cross‐ring cleavages upon electron capture. The modes of ring cleavage depended on the formal oxidation state of the Cu ion and peptide protonation. Density functional theory (DFT) calculations, using the B3LYP with an effective core potential at Cu and M06‐2X functionals, identified several precursor ion structures in which the Cu ion was threecoordinated to pentacoordinated by the His and Lys side‐chain groups and the peptide amide or enolimine groups. The electronic structure of the formally Cu(III) complexes pointed to an effective Cu(I) oxidation state with the other charge residing in the peptide ligand. The relative energies of isomeric complexes of the [Cu(c‐HAHK + H)]²⁺ and [Cu(c‐HAHK ? H)]²⁺ type with closed electronic shells followed similar orders when treated by the B3LYP and M06‐2X functionals. Large differences between relative energies calculated by these methods were obtained for open‐shell complexes of the [Cu(c‐HAHK)]²⁺ type. Charge reduction resulted in lowering the coordination numbers for some Cu complexes that depended on the singlet or triplet spin state being formed. For [Cu(c‐HAHK ? H)]²⁺ complexes, solution H/D exchange involved only the N–H protons, resulting in the exchange of up to seven protons, as established by ultra‐high mass resolution measurements. Contrasting the experiments, DFT calculations found the lowest energy structures for the gas‐phase ions that were deprotonated at the peptide C_α positions. Copyright © 2012 John Wiley & Sons, Ltd.     

Sodium cation affinities of MALDI matrices determined by guided ion beam tandem mass spectrometry: application to benzoic acid derivatives

Threshold collision-induced dissociation of Na(+)(xBA) complexes with Xe is studied using guided ion beam mass spectrometry. The xBA ligands studied include benzoic acid and all of the mono- and dihydroxy-substituted benzoic acids: 2-, 3-, and 4-hydroxybenzoic acid and 2,3-, 2,4-, 2,5-, 2,6-, 3,4-, and 3,5-dihydroxybenzoic acid. In all cases, the primary product corresponds to endothermic loss of the intact xBA ligand. The cross section thresholds are interpreted to yield 0 and 298 K bond dissociation energies (BDEs) for Na(+)-xBA after accounting for the effects of multiple ion-neutral collisions, internal and kinetic energy distributions of the reactants, and dissociation lifetimes. Density functional theory calculations at the B3LYP/6-31G* level of theory are used to determine the structures of these complexes and provide the molecular constants necessary for the thermodynamic analysis of the experimental data. Theoretical BDEs are determined at the B3LYP/6-311+G(2d,2p) and MP2(full)/6-311+G(2d,2p) levels using the B3LYP/6-31G* optimized geometries. The trends in the measured BDEs suggest two very different binding modes for the Na(+)(xBA) complexes, while theory finds four. In general, the most stable binding conformation involves the formation of a six-membered chelation ring via interaction with the carbonyl and 2-hydroxyl oxygen atoms. The ground state geometries of the Na(+)(xBA) complexes in which the ligand does not possess a 2-hydroxyl group generally involve binding of Na(+) to either the carbonyl oxygen atom or to both oxygen atoms of the carboxylic acid group. These binding modes tend to be competitive because the enhancement in binding associated with the chelation interactions in the latter is mediated by steric repulsion between the hydroxyl and ortho hydrogen atoms. When possible, hydrogen bonding interactions with the ring hydroxyl group(s) enhance the stability of these complexes. The agreement between the theoretical and experimental BDEs is quite good for B3LYP and somewhat less satisfactory for MP2(full).     

Influence of thioketo substitution on the properties of uracil and its noncovalent interactions with alkali metal ions: threshold collision-induced dissociation and theoretical studies

Experimental and theoretical studies are carried out to determine the influence of thioketo substitution on the properties of uracil and its noncovalent interactions with alkali metal ions. Bond dissociation energies of alkali metal ion-thiouracil complexes, M(+)(SU), are determined using threshold collision-induced dissociation techniques in a guided ion beam mass spectrometer, where M(+) = Li(+), Na(+), and K(+) and SU = 2-thiouracil, 4-thiouracil, 2,4-dithiouracil, 5-methyl-2-thiouracil, and 6-methyl-2-thiouracil. Ab initio electronic structure calculations are performed to determine the structures and theoretical bond dissociation energies of these complexes and provide molecular constants necessary for thermodynamic analysis of the experimental data. Theoretical calculations are also performed to examine the influence of thioketo substitution on the acidities, proton affinities, and A::SU Watson-Crick base pairing energies. In general, thioketo substitution leads to an increase in both the proton affinity and the acidity of uracil. 2-Thio substitution generally results in an increase in the alkali metal ion binding affinities but has almost no affect on the stability of the A::SU base pair. In contrast, 4-thio substitution results in a decrease in the alkali metal ion binding affinities and a significant decrease in the stability of the A::SU base pair. In addition, alkali metal ion binding is expected to lead to an increase in the stability of both single-stranded and double-stranded nucleic acids by reducing the charge on the nucleic acid in a zwitterion effect as well as through additional noncovalent interactions between the alkali metal ion and the nucleobases.     

Water molecule adsorption on protonated dipeptides

Equilibrium constants for the adsorption of the first water molecule on six protonated dipeptides (Gly-Gly+H(+), Gly-Ala+H(+), Ala-Gly+H(+), Ala-Ala+H(+), Pro-Gly+H(+), and Gly-Trp+H(+)) have been measured as a function of temperature, and DeltaH(o) and DeltaS(o) determined. Density functional theory calculations were performed for both the unsolvated peptides and the peptide water complexes at the B3LYP/6-311++G level. MP2/6-311++G** calculations were also carried out for Gly/Ala peptides. The calculations suggest that adsorption of a water molecule by these simple dipeptides is a complex process, both the unsolvated peptide and the peptide-water complexes have multiple conformations with similar free energies. Average DeltaH(o) and DeltaS(o) values derived from the calculations are in reasonable agreement with the experimental results. According to the calculations, the dominant water adsorption process involves a significant conformational change to accommodate a bridging water molecule. DeltaH(o) is diminished for Pro-Gly+H(+) mainly because the water interacts with a secondary amine, whereas for Gly-Trp+H(+), DeltaH(o) is significantly decreased by the loss of cation-pi interactions upon water adsorption. For unsolvated peptides the proton affinities of the N-terminus and the backbone carbonyl groups are known to be similar. Addition of a single water molecule causes a significant stabilization of the N-terminus protonation site.