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

The binding of Na+ 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 sodiated amino acid 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 the complexes after accounting for unimolecular decay rates, internal energy of the reactant ions, and multiple ion-molecule collisions. Quantum chemical calculations for a number of geometric conformations of each Na+(L) complex 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 the energetic contributions of individual functionalities as well as steric influences of relative chain lengths to be thoroughly explored. Na+ binding affinities for the amide complexes are systematically stronger than those for the acid complexes by 14 +/- 1 kJ/mol, which is attributed to an inductive effect of the OH group in the carboxylic acid side chain. Additionally, the Na+ binding affinity for the longer-chain amino acids (Glx) is enhanced by 4 +/- 1 kJ/mol compared to the shorter-chain Asx because steric effects are reduced.     

Experimental and theoretical studies of sodium cation complexes of the deamidation and dehydration products of asparagine, glutamine, aspartic acid, and glutamic acid

The deamidation and dehydration products of Na+(L), where L = asparagine (Asn), glutamine (Gln), aspartic acid (Asp), and glutamic acid (Glu), are examined in detail utilizing collision-induced dissociation (CID) with Xe in a guided ion beam tandem mass spectrometer (GIBMS). Results establish that the Na+(L) complexes decompose upon formation in our dc discharge/flow tube ion source to form a bis-ligand complex, Na+(L-HX)(HX), composed of a sodium cation, the (L-HX) decomposition product, and HX, where HX = NH3 for the amides and H2O for the acids. Analysis of the energy-dependent CID cross sections for the Na+(L-HX)(HX) complexes provides unambiguous identification of the (L-HX) fragmentation products as 3-amino succinic anhydride (a-SA) for Asx and oxo-proline (O-Pro) for Glx. Furthermore, these experiments establish the 0 K sodium cation affinities for these five-membered ring decomposition products and the H2O and NH3 binding affinities of the Na+(a-SA) and Na+(O-Pro) complexes after accounting for unimolecular decay rates, the internal energy of reactant ions, and multiple ion-molecule collisions. Quantum chemical calculations are determined for a number of geometric conformations of all reaction species as well as a number of candidate species for (L-HX) 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 the experimental work and quantum chemical calculations allows for a complete characterization of the products of deamidation and dehydration of Asx and Glx, as well as the details of Na+, H2O, and NH3 binding to the decomposition species.     

Cation-pi interactions: structures and energetics of complexation of Na+ and K+ with the aromatic amino acids, phenylalanine, tyrosine, and tryptophan

Threshold collision-induced dissociation of M(+)(AAA) with Xe is studied using guided ion beam tandem mass spectrometry. M(+) include the alkali metal ions Na(+) and K(+). The three aromatic amino acids are examined, AAA = phenylalanine, tyrosine, or tryptophan. In all cases, endothermic loss of the intact aromatic amino acid is the dominant reaction pathway. The threshold regions of the cross sections are interpreted to extract 0 and 298 K bond dissociation energies for the M(+)-AAA complexes after accounting for the effects of multiple ion-neutral collisions, internal energy of the reactant ions, and dissociation lifetimes. Density functional theory calculations at the B3LYP/6-31G level of theory are used to determine the structures of the neutral aromatic amino acids and their complexes to Na(+) and K(+) and to provide molecular constants required for the thermochemical analysis of the experimental data. Theoretical bond dissociation energies are determined from single-point energy calculations at the B3LYP/6-311++G(3df,3pd) level using the B3LYP/6-31G geometries. Good agreement between theory and experiment is found for all systems. The present results are compared to earlier studies of these systems performed via kinetic and equilibrium methods. The present results are also compared to the analogous Na(+) and K(+) complexes to glycine, benzene, phenol, and indole to elucidate the relative contributions that each of the functional components of these aromatic amino acids make to the overall binding in these complexes.     

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.     

The special five-membered ring of proline: An experimental and theoretical investigation of alkali metal cation interactions with proline and its four- and six-membered ring analogues

The interaction of the alkali metal cations, Li+, Na+, and K+, with the amino acid proline (Pro) and its four- and six-membered ring analogues, azetidine-2-carboxylic acid (Aze) and pipecolic acid (Pip), are examined in detail. Experimentally, threshold collision-induced dissociation of the M+(L) complexes, where M = Li, Na, and K and L = Pro, Aze, and Pip, with Xe are studied using a guided ion beam tandem mass spectrometer. From analysis of the kinetic energy dependent cross sections, M(+)-L bond dissociation energies are measured. These analyses account for unimolecular decay rates, internal energy of reactant ions, and multiple ion-molecule collisions. Ab initio calculations for a number of geometric conformations of the M+(L) complexes were determined at the B3LYP/6-311G(d,p) level with single-point energies calculated at MP2(full), B3LYP, and B3P86 levels using a 6-311+G(2d,2p) basis set. Theoretical bond energies show good agreement with the experimental bond energies, which establishes that the zwitterionic form of the alkali metal cation/amino acid, the lowest energy conformation, is formed in all cases. Despite the increased conformational mobility in the Pip systems, the Li+, Na+, and K+ complexes of Pro show higher binding energies. A meticulous examination of the zwitterionic structures of these complexes provides an explanation for the stability of the five-membered ring complexes.     

Experimental and theoretical investigation of alkali metal cation interactions with hydroxyl side-chain amino acids

Absolute bond dissociation energies of serine (Ser) and threonine (Thr) to alkali metal cations are determined experimentally by threshold collision-induced dissociation of M+AA complexes, where M+=Li+, Na+, and K+ and AA=Ser and Thr, with xenon in a guided ion beam tandem mass spectrometer. Experimental results show that the binding energies of both amino acids to the alkali metal cations are very similar to one another and follow the order of Li+>Na+>K+. Quantum chemical calculations at three different levels, B3LYP, B3P86, and MP2(full), using the 6-311+G(2d,2p) basis set with geometries and zero-point energies calculated at the B3LYP/6-311+G(d,p) level show good agreement with the experimental bond energies. Theoretical calculations show that all M+AA complexes have charge-solvated structures (nonzwitterionic) with [CO, N, O] tridentate coordination.     

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.     

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.     

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).     

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.     

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.     

Solvation of copper ions by acetone. structures and sequential binding energies of Cu+(acetone)x,x = 1-4 from collision-induced dissociation and theoretical studies

Collision-induced dissociation of Cu+(acetone)(x), x = 1-4, with Xe is studied as a function of kinetic energy using guided ion beam mass spectrometry. In all cases, the primary and lowest energy dissociation channel observed is endothermic loss of one acetone molecule. The primary cross section thresholds are interpreted to yield 0 and 298 K bond energies after accounting for the effects of multiple ion-neutral collisions, internal energy of the complexes, and dissociation lifetimes. Density functional calculations at the B3LYP/6-31G* level of theory are used to determine the structures of these complexes and provide molecular constants necessary for the thermodynamic analysis of the experimental data. Theoretical bond dissociation energies are determined from single point calculations 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 experimental bond energies determined here are in good agreement with previous experimental measurements made in a high-pressure mass spectrometer for the sum of the first and second bond energy (i.e., Cu+(acetone)2 --> Cu+ + 2 acetone) when these results are properly anchored. The agreement between theory and experiment is reasonable in all cases, but varies both with the size of the cluster and the level of theory employed. B3LYP does an excellent job for the x = 1 and 3 clusters, but is systematically low for the x = 2 and 4 clusters such that the overall trends in sequential binding energies are not parallel. In contrast, all MP2 values are somewhat low, but the overall trends parallel the measured values for all clusters. The trends in the measured Cu+(acetone), binding energies are explained in terms of 4s-3d sigma hybridization effects and ligand-ligand repulsion in the clusters.     

Dipole effects on cation-pi interactions: absolute bond dissociation energies of complexes of alkali metal cations to N-methylaniline and N,N-dimethylaniline

Threshold collision-induced dissociation of M (+)( nMA) x with Xe is studied using guided ion beam mass spectrometry, where nMA = N-methylaniline and N, N-dimethylaniline and x = 1 and 2. M (+) includes the following alkali metal cations: Li (+), Na (+), K (+), Rb (+), and Cs (+). In all cases, the primary dissociation pathway corresponds to the endothermic loss of an intact nMA ligand. The primary cross section thresholds are interpreted to yield 0 and 298 K bond dissociation energies (BDEs) for ( nMA) x-1 M (+)-( nMA) after accounting for the effects of multiple ion-neutral collisions, the internal and kinetic energy distributions of the reactants, and the dissociation lifetimes. Density functional theory calculations at the B3LYP/6-31G* level of theory are used to determine the structures of these complexes, which are also used in single-point calculations at the MP2(full)/6-311+G(2d,2p) level to determine theoretical BDEs. The results of these studies are compared to previous studies of the analogous M (+)(aniline) x complexes to examine the effects of methylation of the amino group on the binding interactions. Comparisons are also made to a wide variety of cation-pi complexes previously studied to elucidate the contributions that ion-dipole, ion-induced-dipole, and ion-quadrupole interactions make to the overall binding.     

Absolute thermodynamic measurements of alkali metal cation interactions with a simple dipeptide and tripeptide

Absolute bond dissociation energies (BDEs) of glycylglycine (GG) and glycylglycylglycine (GGG) to sodium and potassium cations and sequential bond energies of glycine (G) in Na+G2 were determined experimentally by threshold collision-induced dissociation (TCID) in a guided ion beam tandem mass spectrometer. Experimental results showed that the binding energies follow the order of Na+ > K+ and M+GGG > M+GG > M+G. Theoretical calculations at the B3LYP/6-311+G(d) level show that all complexes had charge-solvated structures (nonzwitterionic) with either [CO,CO] bidentate or [N,CO,CO] tridentate coordination for M+GG complexes, [CO,CO,CO] tridentate or [N,CO,CO,CO] tetradentate coordination for M+GGG complexes, and [N,CO,N,CO] tetradentate coordination for Na+G2. Ab initio calculations at three different levels of theory (B3LYP, B3P86, and MP2(full) using the 6-311+G(2d,2p) basis set with geometries and zero-point energies calculated at the B3LYP/6-311+G(d) level) show good agreement with the experimental bond energies. This study demonstrates for the first time that TCID measurements of absolute BDEs can be successfully extended to biological molecules as complex as a tripeptide.     

Potassium cation affinities of matrix assisted laser desorption ionization matrices determined by threshold collision-induced dissociation: application to benzoic acid derivatives

Threshold collision-induced dissociation of K+(xBA) complexes with xenon 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 K+-xBA after accounting for the effects of multiple ion-neutral collisions, the kinetic and internal 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 the xBA ligands and their complexes with K+. Theoretical BDEs are determined from single-point energy calculations at the B3LYP/6-311+G(2d,2p) and MP2(full)/6-311+G(2d,2p) levels using B3LYP/6-31G* optimized geometries. Four favorable binding modes for the K+(xBA) complexes are found. In all complexes to an xBA ligand that does not have a 2-hydroxyl substituent, the most favorable binding mode corresponds to a single interaction with the carbonyl oxygen atom. Formation of a 4-membered ring via chelation interactions with both oxygen atoms of the carboxylic acid group is found to be the most favorable binding mode for all of the 2-hydroxy-substituted systems except K+(2,3-dihydroxybenzoic acid). In these complexes, a hydrogen-bonding interaction between the hydrogen atom of the carboxylic acid moiety and the oxygen atom of the 2-hydroxy substituent provides additional stabilization. Formation of a 5-membered chelation ring via interaction of K+ with the oxygen atoms of adjacent hydroxyl substituents is also favorable and corresponds to the ground-state geometry for the K+(23DHBA) complex. Formation of a 6-membered chelation ring via interaction of K+ with the carbonyl and 2-hydroxyl oxygen atoms is also quite favorable but does not correspond to the ground-state geometry for any of the systems examined here. The experimental BDEs determined here are in very good agreement with the calculated values.     

Structural and energetic effects in the molecular recognition of amino acids by 18-crown-6

Absolute 18-crown-6 (18C6) affinities of five amino acids (AAs) are determined using guided ion beam tandem mass spectrometry techniques. The AAs examined in this work include glycine (Gly), alanine (Ala), lysine (Lys), histidine (His), and arginine (Arg). Theoretical electronic structure calculations are performed to determine stable geometries and energetics for neutral and protonated 18C6 and the AAs as well as the proton bound complexes comprised of these species, (AA)H(+)(18C6). The proton affinities (PAs) of Gly and Ala are lower than the PA of 18C6, whereas the PAs of Lys, His, and Arg exceed that of 18C6. Therefore, the collision-induced dissociation (CID) behavior of the (AA)H(+)(18C6) complexes differs markedly across these systems. CID of the complexes to Gly and Ala produces H(+)(18C6) as the dominant and lowest energy pathway. At elevated energies, H(+)(AA) was produced in competition with H(+)(18C6) as a result of the relatively favorable entropy change in the formation of H(+)(AA). In contrast, CID of the complexes to the protonated basic AAs results in the formation of H(+)(AA) as the only direct CID product. H(+)(18C6) was not observed, even at elevated energies, as a result of unfavorable enthalpy and entropy change associated with its formation. Excellent agreement between the measured and calculated (AA)H(+)-18C6 bond dissociation energies (BDEs) is found with M06 theory for all complexes except (His)H(+)(18C6), where theory overestimates the strength of binding. In contrast, B3LYP theory significantly underestimates the (AA)H(+)-18C6 BDEs in all cases. Among the basic AAs, Lys exhibits the highest binding affinity for 18C6, suggesting that the side chains of Lys residues are the preferred binding site for 18C6 complexation in peptides and proteins. Gly and Ala exhibit greater 18C6 binding affinities than Lys, suggesting that the N-terminal amino group provides another favorable binding site for 18C6. Trends in the 18C6 binding affinities among the five AAs examined here exhibit an inverse correlation with the polarizability and proton affinity of the AA. Therefore, the ability of the N-terminal amino group to compete for 18C6 complexation is best for Gly and should become increasing less favorable as the size of the side chain substituent increases.     

Complexes between lithium cation and diphenylalkanes in the gas phase: the pincer effect

The gas-phase lithium cation basicities (LCB values, Gibbs free energies of binding) of alpha,omega-diphenylalkanes Ph-(CH(2))(n)-Ph (n=2, 3, or 7) and 1,1-diphenylethane Ph-CH(Me)-Ph were investigated by means of Fourier-transform ion cyclotron resonance (FTICR) mass spectrometry. Their structures, and those of the corresponding Li(+) complexes were optimized at the B3LYP/6-31G(d) level and their relative stabilities calculated at the B3LYP/6-311+G(3df,2p)//B3LYP/6-31G(d) level. Whereas the most stable conformers of the free diphenylalkanes were found to adopt a completely stretched aliphatic chain connecting the two benzene rings, the most stable Li(+) complexes correspond to conformers in which the alkali metal cation interacts simultaneously with both benzene rings through the folding of the aliphatic chain ("pincer effect"). This chelation brings about a significant enhancement of the Li(+) binding enthalpies (LBE values), which were calculated to be approximately 75 kJ mol(-1) higher than those evaluated for conventional (singly coordinated) pi complexes in which the metal cation interacts with only one of the benzene rings. The increase of the corresponding lithium cation basicities, however, (Gibbs free energies of Li(+) binding, LCB values) was calculated to be smaller by approximately 15 kJ mol(-1) as the pincer effect is entropically disfavored. The good agreement between the calculated LCB values, assuming a statistical distribution of the different conformers present in the gas phase, and the experimental LCB values measured by means of FTICR mass spectrometry are considered indirect evidence of the existence of the pincer effect.     

Influence of the d orbital occupation on the nature and strength of copper cation-pi interactions: threshold collision-induced dissociation and theoretical studies

Threshold collision-induced dissociation techniques are employed to determine the bond dissociation energies of a wide variety of copper cation-pi complexes, Cu(+)(pi-ligand), where pi-ligand = benzene, flurobenzene, chlorobenzene, bromobenzene, iodobenzene, phenol, toluene, anisole, pyrrole, N-methylpyrrole, indole, naphthalene, aniline, N-methylaniline, and N,N-dimethylaniline. The primary and lowest energy dissociation pathway corresponds to the endothermic loss of the intact neutral pi-ligand for all complexes except those to N-methylpyrrole, indole, aniline, N-methylaniline, and N,N-dimethylaniline. In the latter complexes, the primary dissociation pathway corresponds to loss of the intact ligand accompanied by charge transfer, thereby producing a neutral copper atom and ionized pi-ligand. Fragmentation of the pi-ligands is also observed at elevated energies in several cases. Theoretical calculations at the B3LYP/6-311G(d,p) level of theory are used to determine the structures, vibrational frequencies, and rotational constants of these complexes. Multiple low-energy conformers are found for all of the copper cation-pi complexes. Theoretical bond dissociation energies are determined from single point energy calculations at the B3LYP/6-311+G(3df,2p) level of theory using the B3LYP/6-311G(d,p) optimized geometries. The agreement between theory and experiment is very good for most complexes. The nature and strength of the binding in these copper cation-pi complexes are studied and compared with the corresponding cation-pi complexes to Na(+). Natural bond orbital analyses are carried out to examine the influence of the d orbital occupation on copper cation-pi interactions.     

Experimental validation of theoretical potassium and sodium cation affinities of amides by mass spectrometric kinetic method measurements

In this study the theoretical Gaussian-2 K(+)/Na(+) binding affinities (enthalpies) at 0 K (in kJ mol(-1)) for six amides in the order: formamide (109.2/138.5) < N-methylformamide (117.7/148.6) < acetamide (118.7/149.5) < N,N-dimethylformamide (123.9/156.4) < N-methylacetamide (125.6/157.7) < N,N-dimethylacetamide (129.2/162.6), reported previously (Siu et al., J. Chem. Phys. 2001; 114: 7045-7051), were validated experimentally by mass spectrometric kinetic method measurements. By monitoring the collision-induced dissociation (CID) of K(+)/Na(+)-bound heterodimers of the amides, the relative affinities were shown to be accurate to within +/-2 kJ mol(-1). With these six theoretical K(+)/Na(+) binding affinities as reference values, the absolute K(+)/Na(+) affinities of imidazole, 1-methylimidazole, pyridazine and 1,2-dimethoxyethane were determined by the extended kinetic method, and found to be consistent (to within +/-9 kJ mol(-1)) with literature experimental values obtained by threshold-CID, equilibrium high-pressure mass spectrometry, and Fourier transform ion cyclotron resonance/ligand-exchange equilibrium methods. A self-consistent resolution is proposed for the inconsistencies in the relative order of K(+)/Na(+) affinities of amides reported in the literature. These two sets of validated K(+) and Na(+) affinity values are useful as reference values in kinetic method measurements of K(+)/Na(+) affinity of model biological ligands, such as the K(+) affinities of aliphatic amino acids.