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.     

Competition between pi and non-pi cation-binding sites in aromatic amino acids: a theoretical study of alkali metal cation (Li+, Na+, K+)-phenylalanine complexes

To understand the cation-pi interaction in aromatic amino acids and peptides, the binding of M(+) (where M(+) = Li(+), Na(+), and K(+)) to phenylalanine (Phe) is studied at the best level of density functional theory reported so far. The different modes of M(+) binding show the same order of binding affinity (Li(+)>Na(+)>K(+)), in the approximate ratio of 2.2:1.5:1.0. The most stable binding mode is one in which the M(+) is stabilized by a tridentate interaction between the cation and the carbonyl oxygen (O[double bond]C), amino nitrogen (--NH(2)), and aromatic pi ring; the absolute Li(+), Na(+), and K(+) affinities are estimated theoretically to be 275, 201, and 141 kJ mol(-1), respectively. Factors affecting the relative stabilities of various M(+)-Phe binding modes and conformers have been identified, with ion-dipole interaction playing an important role. We found that the trend of pi and non-pi cation bonding distances (Na(+)-pi>Na(+)-N>Na(+)-O and K(+)-pi>K(+)-N>K(+)-O) in our theoretical Na(+)/K(+)-Phe structures are in agreement with the reported X-ray crystal structures of model synthetic receptors (sodium and potassium bound lariat ether complexes), even though the average alkali metal cation-pi distance found in the crystal structures is longer. This difference between the solid and the gas-phase structures can be reconciled by taking the higher coordination number of the cations in the lariat ether complexes into account.     

Thermodynamic and ab initio analysis of the controversial enthalpy of formation of formaldehyde.

There are two values, -26.0 and -27.7 kcal mol(-1), that are routinely reported in literature evaluations for the standard enthalpy of formation, Delta(f) H(o)(298), of formaldehyde (CH(2)=O), where error limits are less than the difference in values. In this study, we summarize the reported literature for formaldehyde enthalpy values based on evaluated measurements and on computational studies. Using experimental reaction enthalpies for a series of reactions involving formaldehyde, in conjunction with known enthalpies of formation, its enthalpy is determined to be -26.05+/-0.42 kcal mol(-1), which we believe is the most accurate enthalpy currently available. For the same reaction series, the reaction enthalpies are evaluated using six computational methods: CBS-Q, CBS-Q//B3, CBS-APNO, G2, G3, and G3B3 yield Delta(f) H(o)(298)=-25.90+/-1.17 kcal mol(-1), which is in good agreement to our experimentally derived result. Furthermore, the computational chemistry methods G3, G3MP2B3, CCSD/6-311+G(2df,p)//B3LYP/6-31G(d), CCSD(T)/6-311+G(2df,p)//B3LYP/6-31G(d), and CBS-APNO in conjunction with isodesmic and homodesmic reactions are used to determine Delta(f) H(o)(298). Results from a series of five work reactions at the higher levels of calculation are -26.30+/-0.39 kcal mol(-1) with G3, -26.45+/-0.38 kcal mol(-1) with G3MP2B3, -26.09+/-0.37 kcal mol(-1) with CBS-APNO, -26.19+/-0.48 kcal mol(-1) with CCSD, and -26.16+/-0.58 kcal mol(-1) with CCSD(T). Results from heat of atomization calculations using seven accurate ab initio methods yields an enthalpy value of -26.82+/-0.99 kcal mol(-1). The results using isodesmic reactions are found to give enthalpies more accurate than both other computational approaches and are of similar accuracy to atomization enthalpy calculations derived from computationally intensive W1 and CBS-APNO methods. Overall, our most accurate calculations provide an enthalpy of formation in the range of -26.2 to -26.7 kcal mol(-1), which is within computational error of the suggested experimental value. The relative merits of each of the three computational methods are discussed and depend upon the accuracy of experimental enthalpies of formation required in the calculations and the importance of systematic computational errors in the work reaction. Our results also calculate Delta(f) H(o)(298) for the formyl anion (HCO(-)) as 1.28+/-0.43 kcal mol(-1).     

Selective binding of monovalent cations to the stacking G-quartet structure formed by guanosine 5'-monophosphate: a solid-state NMR study

We report a solid-state multinuclear ((23)Na, (15)N, (13)C, and (31)P) NMR study on the relative affinity of monovalent cations for a stacking G-quartet structure formed by guanosine 5'-monophosphate (5'-GMP) self-association at pH 8. Two major types of cations are bound to the 5'-GMP structure: one at the surface and the other within the channel cavity between two G-quartets. The channel cation is coordinated to eight carbonyl oxygen atoms from the guanine bases, whereas the surface cation is close to the phosphate group and likely to be only partially hydrated. On the basis of solid-state (23)Na NMR results from a series of ion titration experiments, we have obtained quantitative thermodynamic parameters concerning the relative cation binding affinity for each of the two major binding sites. For the channel cavity site, the values of the free energy difference (Delta G degrees at 25 degrees C) for ion competition between M(+) and Na(+) ions are K(+) (-1.9 kcal mol(-1)), NH(4)(+) (-1.8 kcal mol(-1)), Rb(+) (-0.3 kcal mol(-1)), and Cs(+) (1.8 kcal mol(-1)). For the surface site, the values Delta G degrees are K(+) (2.5 kcal mol(-1)), NH(4)(+) (-1.3 kcal mol(-1)), Rb(+) (1.1 kcal mol(-1)), and Cs(+) (0.9 kcal mol(-1)). Solid-state NMR data suggest that the affinity of monovalent cations for the 5'-GMP structure follows the order NH(4)(+) > Na(+) > Cs(+) > Rb(+) > K(+) at the surface site and K(+) > NH(4)(+) > Rb(+) > Na(+) > Cs(+) > Li(+) at the channel cavity site. We have found that the cation-induced stability of a 5'-GMP structure is determined only by the affinity of monovalent cations for the channel site and that the binding of monovalent cations to phosphate groups plays no role in 5'-GMP self-ordered structure. We have demonstrated that solid-state (23)Na and (15)N NMR can be used simultaneously to provide mutually complementary information about competitive binding between Na(+) and NH(4)(+) ions.     

Enthalpies of formation of gas-phase N3, N3-, N5+, and N5- from Ab initio molecular orbital theory, stability predictions for N5(+)N3(-) and N5(+)N5(-), and experimental evidence for the instability of N5(+)N3(-)

Ab initio molecular orbital theory has been used to calculate accurate enthalpies of formation and adiabatic electron affinities or ionization potentials for N3, N3-, N5+, and N5- from total atomization energies. The calculated heats of formation of the gas-phase molecules/ions at 0 K are DeltaHf(N3(2Pi)) = 109.2, DeltaHf(N3-(1sigma+)) = 47.4, DeltaHf(N5-(1A1')) = 62.3, and DeltaHf(N5+(1A1)) = 353.3 kcal/mol with an estimated error bar of +/-1 kcal/mol. For comparison purposes, the error in the calculated bond energy for N2 is 0.72 kcal/mol. Born-Haber cycle calculations, using estimated lattice energies and the adiabatic ionization potentials of the anions and electron affinities of the cations, enable reliable stability predictions for the hypothetical N5(+)N3(-) and N5(+)N5(-) salts. The calculations show that neither salt can be stabilized and that both should decompose spontaneously into N3 radicals and N2. This conclusion was experimentally confirmed for the N5(+)N3(-) salt by low-temperature metathetical reactions between N5SbF6 and alkali metal azides in different solvents, resulting in violent reactions with spontaneous nitrogen evolution. It is emphasized that one needs to use adiabatic ionization potentials and electron affinities instead of vertical potentials and affinities for salt stability predictions when the formed radicals are not vibrationally stable. This is the case for the N5 radicals where the energy difference between vertical and adiabatic potentials amounts to about 100 kcal/mol per N5.     

Insight into the cation-π interaction at the metal binding site of the copper metallochaperone CusF

The periplasmic Cu(+)/Ag(+) chaperone CusF features a novel cation-π interaction between a Cu(+)/Ag(+) ion and Trp44 at the metal binding site. The nature and strength of the Cu(+)/Ag(+)-Trp44 interactions were investigated using computational methodologies. Quantum-mechanical (QM) calculations showed that the Cu(+) and Ag(+) interactions with Trp44 are of similar strength (~14 kcal/mol) and bond order. Quantum-mechanical/molecular-mechanical (QM/MM) calculations showed that Cu(+) binds in a distorted tetrahedral coordination environment in the Trp44Met mutant, which lacks the cation-π interaction. Molecular dynamics (MD) simulations of CusF in the apo and Cu(+)-bound states emphasized the importance of the Cu(+)-Trp44 interaction in protecting Cu(+) from water oxidation. The protein structure does not change over the time scale of hundreds of nanoseconds in the metal-bound state. The metal recognition site exhibits small motions in the apo state but remains largely preorganized toward metal binding. Trp44 remains oriented to form the cation-π interaction in the apo state and faces an energetic penalty to move away from the metal ion. Cu(+) binding quenches the protein's internal motions in regions linked to binding CusB, suggesting that protein motions play an essential role in Cu(+) transfer to CusB.     

Car-Parrinello MD simulations for the Na+ -phenylalanine complex in aqueous solution

Car-Parrinello molecular dynamics (CPMD) calculations are presented for a Na (+)(Phe) complex in aqueous solution and for various stable Na (+)(Phe) complexes and Na (+)(H 2O) n clusters in the gas phase (with up to six water molecules). The CPMD results are compared to available experimental and ab initio reference data, to DFT results obtained with various combinations of density functionals and basis sets, and to previous classical mechanics MD simulations. The agreement with the reference data in the gas phase validates the CPMD method, showing that it is a valid approach for studying these systems and that it describes correctly the competing Na (+)-Phe and Na (+)-H 2O interactions. Analysis of MD trajectories reveals that the Na (+)(Phe) complex in aqueous solution maintains a stable configuration in which the Na (+) cation hovers above the phenyl ring, at an average distance of 3.85 A from the ring center, while remaining strongly bound to one of the carboxylic oxygens of Phe. Constrained MD simulations indicate that the free energy barrier opposing dissociation of the complex exceeds 5.5 kcal/mol. We thus confirm that "cation- pi" interactions between alcali cations and the pi ring, combined with other kinds of interactions, may allow aromatic amino acids to overcome the competition with water in binding a cation.     

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.     

Design of an adenosine analogue that selectively improves the affinity of a mutant U1A protein for RNA

The RNA recognition motif (RRM), one of the most common RNA binding domains, contains three highly conserved aromatic amino acids that participate in stacking interactions with RNA bases. We have investigated the contribution of these highly conserved aromatic amino acids to the affinity of the complex formed between the N-terminal RRM of the U1A protein and stem loop 2 of U1 snRNA. Previously, we found that substitution of one of these conserved aromatic amino acids, Phe56, with Ala resulted in a large destabilization of the complex. Here, we have modified A6, the base in stem loop 2 RNA that stacks with Phe56, to compensate for a portion of the destabilization caused by the Phe56Ala mutation. We have designed two modified adenosines, A-3CPh and A-4CPh, in which a phenyl group is linked to the adenosine such that it may replace the phenyl group that is eliminated by the Phe56Ala mutation in the complex. We have found that incorporation of A-3CPh into stem loop 2 RNA stabilizes the complex formed with Phe56Ala by 0.6 kcal/mol, while incorporation of A-4CPh into stem loop 2 RNA stabilizes this complex by 1.8 kcal/mol. Either base modification destabilizes the wild-type complex by 0.8-0.9 kcal/mol. Experiments with other U1A mutant proteins suggest that the stabilization of the complex between the Phe56Ala U1A protein and stem loop 2 RNA is due to a specific interaction between the Phe56Ala U1A protein and A6-4CPh stem loop 2 RNA.     

The heat of formation of the uranyl dication: theoretical evaluation based on relativistic density functional calculations

By using a set of model reactions, we estimated the heat of formation of gaseous UO2(2+) from quantum-chemical reaction enthalpies and experimental heats of formation of reference species. For this purpose, we performed relativistic density functional calculations for the molecules UO2(2+), UO2, UF6, and UF5. We used two gradient-corrected exchange-correlation functionals (revised Perdew-Burke-Ernzerhof (PBEN) and Becke-Perdew (BP)) and we accounted for spin-orbit interaction in a self-consistent fashion. Indeed, spin-orbit interaction notably affects the energies of the model reactions, especially if compounds of U(IV) are involved. Our resulting theoretical estimates for delta fH(o)0 (UO2(2+)), 365+/-10 kcal mol(-1) (PBEN) and 370+/-12 kcal mol(-1) (BP), are in quantitative agreement with a recent experimental result, 364+/-15 kcal mol(-1). Agreement between the results of the two different exchange-correlation functionals PBEN and BP supports the reliability of our approach. The procedure applied offers a general means to derive unknown enthalpies of formation of actinide species based on the available well-established data for other compounds of the element in question.     

Accurate computational determination of the binding energy of the SO3 x H2O complex

Reliable thermochemical data for the reaction SO3 + H2O<-->SO3 x H2O (1a) are of crucial importance for an adequate modeling of the homogeneous H2SO4 formation in the atmosphere. We report on high-level quantum chemical calculations to predict the binding energy of the SO3 x H2O complex. The electronic binding energy is accurately computed to De = 40.9+/-1.0 kJ/mol = 9.8+/-0.2 kcal/mol. By using harmonic frequencies from density functional theory calculations (B3LYP/cc-pVTZ and TPSS/def2-TZVP), zero-point and thermal energies were calculated. From these data, we estimate D0 = -Delta H(1a)0(0 K) = 7.7+/-0.5 kcal/mol and Delta H(1a)0(298 K) = -8.3+/-1.0 kcal/mol.     

Negative-ion photoelectron spectroscopy, gas-phase acidity, and thermochemistry of the peroxyl radicals CH(3)OO and CH(3)CH(2)OO

Methyl, methyl-d(3), and ethyl hydroperoxide anions (CH(3)OO(-), CD(3)OO(-), and CH(3)CH(2)OO(-)) have been prepared by deprotonation of their respective hydroperoxides in a stream of helium buffer gas. Photodetachment with 364 nm (3.408 eV) radiation was used to measure the adiabatic electron affinities: EA[CH(3)OO, X(2)A' '] = 1.161 +/- 0.005 eV, EA[CD(3)OO, X(2)A' '] = 1.154 +/- 0.004 eV, and EA[CH(3)CH(2)OO, X(2)A' '] = 1.186 +/- 0.004 eV. The photoelectron spectra yield values for the term energies: Delta E(X(2)A' '-A (2)A')[CH(3)OO] = 0.914 +/- 0.005 eV, Delta E(X(2)A' '-A (2)A')[CD(3)OO] = 0.913 +/- 0.004 eV, and Delta E(X(2)A' '-A (2)A')[CH(3)CH(2)OO] = 0.938 +/- 0.004 eV. A localized RO-O stretching mode was observed near 1100 cm(-1) for the ground state of all three radicals, and low-frequency R-O-O bending modes are also reported. Proton-transfer kinetics of the hydroperoxides have been measured in a tandem flowing afterglow-selected ion flow tube (FA-SIFT) to determine the gas-phase acidity of the parent hydroperoxides: Delta(acid)G(298)(CH(3)OOH) = 367.6 +/- 0.7 kcal mol(-1), Delta(acid)G(298)(CD(3)OOH) = 367.9 +/- 0.9 kcal mol(-1), and Delta(acid)G(298)(CH(3)CH(2)OOH) = 363.9 +/- 2.0 kcal mol(-1). From these acidities we have derived the enthalpies of deprotonation: Delta(acid)H(298)(CH(3)OOH) = 374.6 +/- 1.0 kcal mol(-1), Delta(acid)H(298)(CD(3)OOH) = 374.9 +/- 1.1 kcal mol(-1), and Delta(acid)H(298)(CH(3)CH(2)OOH) = 371.0 +/- 2.2 kcal mol(-1). Use of the negative-ion acidity/EA cycle provides the ROO-H bond enthalpies: DH(298)(CH(3)OO-H) = 87.8 +/- 1.0 kcal mol(-1), DH(298)(CD(3)OO-H) = 87.9 +/- 1.1 kcal mol(-1), and DH(298)(CH(3)CH(2)OO-H) = 84.8 +/- 2.2 kcal mol(-1). We review the thermochemistry of the peroxyl radicals, CH(3)OO and CH(3)CH(2)OO. Using experimental bond enthalpies, DH(298)(ROO-H), and CBS/APNO ab initio electronic structure calculations for the energies of the corresponding hydroperoxides, we derive the heats of formation of the peroxyl radicals. The "electron affinity/acidity/CBS" cycle yields Delta(f)H(298)[CH(3)OO] = 4.8 +/- 1.2 kcal mol(-1) and Delta(f)H(298)[CH(3)CH(2)OO] = -6.8 +/- 2.3 kcal mol(-1).     

Accurate measurement of the relative bond energies of CO and H2O ligands in Fe+ mono- and bis-ligated complexes

The systems Fe(H(2)O)(n) (+)/CO[bond]H(2)O and Fe(CO)(n) (+)/CO[bond]H(2)O (n = 1 and 2) were investigated in a triple cell Fourier transform ion cyclotron resonance (FT-ICR) mass spectrometer. Using mixtures of CO with a very small amount of water, the ligand exchange equilibrium was reached, allowing experimental determination of the relevant equilibrium constants and free energies of reaction. Quantum chemical calculations at the B3LYP level of theory on the reactant and product species allowed us to determine the entropic terms and to derive the relative bond energies of CO and H(2)O in the mono- and bis-ligated complexes. For n = 1, H(2)O is more strongly bound to Fe(+) than CO by 4.1 +/- 1.6 kJ x mol(-1) at 298 K. For n = 2, at the same temperature, H(2)O is more strongly bound than CO to (H(2)O)Fe(+) by 7.6 +/- 1.6 kJ x mol(-1), and to (CO)Fe(+) by more than 20.1 kJ x mol(-1).     

A solid-state 23Na NMR study of monovalent cation binding to double-stranded DNA at low relative humidity

We report a solid-state (23)Na NMR study of monovalent cation (Li(+), Na(+), K(+), Rb(+), Cs(+) and NH(4) (+)) binding to double-stranded calf thymus DNA (CT DNA) at low relative humidity, ca 0-10%. Results from (23)Na--(31)P rotational echo double resonance (REDOR) NMR experiments firmly establish that, at low relative humidity, monovalent cations are directly bound to the phosphate group of CT DNA and are partially dehydrated. On the basis of solid-state (23)Na NMR titration experiments, we obtain quantitative thermodynamic parameters concerning the cation-binding affinity for the phosphate group of CT DNA. The free energy difference (DeltaG degrees ) between M(+) and Na(+) ions is as follows: Li(+) (-1.0 kcal mol(-1)), K(+) (7.2 kcal mol(-1)), NH(4) (+) (1.0 kcal mol(-1)), Rb(+) (4.5 kcal mol(-1)) and Cs(+) (1.5 kcal mol(-1)). These results suggest that, at low relative humidity, the binding affinity of monovalent cations for the phosphate group of CT DNA follows the order: Li(+) > Na(+) > NH(4) (+) > Cs(+) > Rb(+) > K(+). This sequence is drastically different from that observed for CT DNA in solution. This discrepancy is attributed to the different modes of cation binding in dry and wet states of DNA. In the wet state of DNA, cations are fully hydrated. Our results suggest that the free energy balance between direct cation-phosphate contact and dehydration interactions is important. The reported experimental results on relative ion-binding affinity for the DNA backbone may be used for testing theoretical treatment of cation-phosphate interactions in DNA.     

Composition and thermochemistry of silver bromide vapor

The mole fractions of AgBr and Ag3Br3 in the saturated vapor at 840 K have been evaluated from the vapor mass spectrum, by comparison with the corresponding spectrum of AgCl vapor, where the monomer/trimer ratio is known accurately from vapor molecular weight measurements. Combination of these results with new measurements of the vapor pressure of molten AgBr by the torsion-effusion method in the range 805-936 K yielded the third law enthalpies of vaporization and the standard enthalpies of formation DeltafH degrees 298(AgBr, g) = 27.8 +/- 0.3 kcal mol(-1) and DeltafH degrees 298(Ag3Br3, g) = -19.0 +/- 1 kcal mol(-1). The dissociation energy, D degrees 0(AgBr), is found to be 66.4 +/- 0.3 kcal mol(-1), or 2.88 +/- 0.01 eV, some 3.5-5 kcal mol(-1) lower than previous literature values. Approximate thermochemical stabilities of the dimer species Ag2Cl2 and Ag2Br2 have also been evaluated.     

The role of cation-pi interactions in biomolecular association. Design of peptides favoring interactions between cationic and aromatic amino acid side chains.

Cation-pi interactions between amino acid side chains are increasingly being recognized as important structural and functional features of proteins and other biomolecules. Although these interactions have been found in static protein structures, they have not yet been detected in dynamic biomolecular systems. We determined, by (1)H NMR spectroscopic titrations, the energies of cation-pi interactions of the amino acid derivative AcLysOMe (1) with AcPheOEt (2) and with AcTyrOEt (3) in aqueous and three organic solvents. The interaction energy is substantial; it ranges from -2.1 to -3.4 kcal/mol and depends only slightly on the dielectric constant of the solvent. To assess the effects of auxiliary interactions and structural preorganization on formation of cation-pi interactions, we studied these interactions in the association of pentapeptides. Upon binding of the positively-charged peptide AcLysLysLysLysLysNH(2) (5) to the negatively-charged partner AcAspAspXAspAspNH(2) (6), in which X is Leu (6a), Tyr (6b), and Phe (6c), multiple interactions occur. Association of the two pentapeptides is dynamic. Free peptides and their complex are in fast exchange on the NMR time-scale, and 2D (1)H ROESY spectra of the complex of the two pentapeptides do not show intermolecular ROESY peaks. Perturbations of the chemical shifts indicated that the aromatic groups in peptides 6b and 6c were affected by the association with 5. The association constants K(A) for 5 with 6a and with 6b are nearly equal, (4.0 +/- 0.7) x 10(3) and (5.0 +/- 1.0) x 10(3) M(-)(1), respectively, while K(A) for 5 with 6c is larger, (8.3 +/- 1.3) x 10(3) M(-)(1). Molecular-dynamics (MD) simulations of the pentapeptide pairs confirmed that their association is dynamic and showed that cation-pi contacts between the two peptides are stereochemically possible. A transient complex between 5 and 6 with a prominent cation-pi interaction, obtained from MD simulations, was used as a template to design cyclic peptides C(X) featuring persistent cation-pi interactions. The cyclic peptide C(X) had a sequence in which X is Tyr, Phe, and Leu. The first two peptides do, but the third does not, contain the aromatic residue capable of interacting with a cationic Lys residue. This covalent construct offered conformational stability over the noncovalent complexes and allowed thorough studies by 2D NMR spectroscopy. Multiple conformations of the cyclic peptides C(Tyr) and C(Phe) are in slow exchange on the NMR time-scale. In one of these conformations, cation-pi interaction between Lys3 and Tyr9/Phe9 is clearly evident. Multiple NOEs between the side chains of residues 3 and 9 are observed; chemical-shift changes are consistent with the placement of the side chain of Lys3 over the aromatic ring. In contrast, the cyclic peptide C(Leu) showed no evidence for close approach of the side chains of Lys3 and Leu9. The cation-pi interaction persists in both DMSO and aqueous solvents. When the disulfide bond in the cyclic peptide C(Phe) was removed, the cation-pi interaction in the acyclic peptide AC(Phe) remained. To test the reliability of the pK(a) criterion for the existence of cation-pi interactions, we determined residue-specific pK(a) values of all four Lys side chains in all three cyclic peptides C(X). While NOE cross-peaks and perturbations of the chemical shifts clearly show the existence of the cation-pi interaction, pK(a) values of Lys3 in C(Tyr) and in C(Phe) differ only marginally from those values of other lysines in these dynamic peptides. Our experimental results with dynamic peptide systems highlight the role of cation-pi interactions in both intermolecular recognition at the protein-protein interface and intramolecular processes such as protein folding.     

Cation mobility and kinetics of ion exchange in zirconium hydrogen monothiophosphate hydrate, Zr(HPO(3)S)(2)x1.5H(2)O

The ion conductivity of zirconium hydrogen monothiophosphate (Zr(HPO(3)S)(2)x1.5H(2)O) has been measured by impedance spectroscopy. The measured value of proton conductivity is 3 x 10(-5) S/cm at 298 K. Conductivity was shown to decrease with increasing temperature due to a dehydration process. Above 450 K, the conductivity is likely governed by proton transport in the anhydrous phase Zr(HPO(3)S)(2). The activation energies of proton conductivity were measured to be 18 +/- 2 kJ/mol for Zr(HPO(3)S)(2)x1.5H(2)O and 60 +/- 3 kJ/mol for the anhydrous compound. The kinetics of ion exchange was studied with the use of potentiometric titration for several ion pairs, H(+)/Na(+), H(+)/Zn(2+), and Na(+)/Zn(2+) in Zr(HPO(3)S)(2)x1.5H(2)O. The diffusion coefficient values for H(+)/Na(+) ion exchange in Zr(HPO(3)S)(2)x1.5H(2)O are lower than those reported in alpha-zirconium phosphate. At the same time, the mobility of zinc ions in Zr(HPO(3)S)(2)x1.5H(2)O is higher than sodium ion mobility. The ion exchange H(+)/Zn(2+) is accompanied by the slow hydrolysis of the initial compound. In all cases, the powdered solids were evaluated by powder X-ray diffraction, and particle sizes were controlled by grinding and sieving the powders.     

Photoelectron spectroscopy and density functional theory of puckered ring structures of Group 13 metal-ethylenediamine

The ethylenediamine (en) complexes of Al, Ga, and In atoms were prepared in laser-vaporization supersonic molecular beams and studied with pulsed field ionization zero electron kinetic energy photoelectron spectroscopy and density functional theory. Several conformers of each metal complex are obtained by B3LYP calculations, and a five-membered cyclic structure is identified by combining the experimental measurements and theoretical calculations. Adiabatic ionization potentials, vibrational frequencies, and bond dissociation energies are determined for the ring structure. The ionization potentials of the Al, Ga, and In species are measured to be 32 784 (5), 33 324 (5), and 33 637 (7) cm(-1), respectively, and metal-ligand dissociation energies of the ionic and neutral complexes are calculated to be 60.2/16.2 (Al(+)/Al), 55.5/13.0 (Ga(+)/Ga), and 50.0/11.4 (In(+)/In) kcal mol(-1). Metal-ligand stretch and bend as well as a number of ligand-based vibrations are measured. Harmonic frequencies and anharmonicities of the M(+)-N (M=Al,Ga,In) stretch are determined for all three M(+)-en ions and the C-C-N bend of Ga(+)-en and In(+)-en. In comparison to monodentate methylamine, the bidentate binding of ethylenediamine leads to a significantly lower ionization potential and higher metal-ligand bond strength of the metal complexes.     

A joint experimental and theoretical study of cation-pi interactions: multiple-decker sandwich complexes of ferrocene with alkali metal ions (Li+, Na+, K+, Rb+, Cs+)

A systematic study of cation-pi interactions between alkali metal ions and the cyclopentadienyl ring of ferrocene is presented. The alkali metal (Li+, Na+, K+, Rb+, Cs+) salts of the ditopic mono(pyrazol-1-yl)borate ligand [1,1'-fc(BMe2pz)2]2- crystallize from dimethoxyethane as multiple-decker sandwich complexes with the M+ ions bound to the pi faces of the ferrocene cyclopentadienyl rings in an eta5 manner (fc = (C5H4)2Fe; pz = pyrazolyl). X-ray crystallography of the lithium complex reveals discrete trimetallic entities with each lithium ion being coordinated by only one cyclopentadienyl ring. The sodium salt forms polyanionic zigzag chains where each Na+ ion bridges the cyclopentadienyl rings of two ferrocene moieties. Linear columns [-CpR-Fe-CpR-M+-CpR-Fe-CpR-M+-](infinity) (R = [-BMe2pz]-) are established by the K+, Rb+, and Cs+ derivatives in the solid state. According to DFT calculations, the binding enthalpies of M+-eta5(ferrocene) model complexes are about 20% higher as compared to the corresponding M+-eta6(benzene) aggregates when M+ = Li+ or Na+. For K+ and Rb+, the degree of cation-pi interaction with both aromatics is about the same. The binding sequence along the M+-eta5(ferrocene) series follows a classical electrostatic trend with the smaller ions being more tightly bound.     

Arene-cation interactions of positive quadrupole moment aromatics and arene-anion interactions of negative quadrupole moment aromatics

Intermolecular interactions involving aromatic pi-electron density are widely believed to be governed by the aromatic molecular quadrupole moment, Theta(zz). Arene-cation binding is believed to occur primarily with negative Theta(zz) aromatics, and arene-anion binding is believed to occur largely with positive Theta(zz) aromatics. We have performed quantum mechanical computations that show the cation binding of positive Theta(zz) aromatics and the anion binding of negative Theta(zz) aromatics is quite common in the gas phase. The pi-electron density of hexafluorobenzene, the prototypical positive Theta(zz ) aromatic (experimental Theta(zz) = 9.5 +/- 0.5 DA), has a Li+ binding enthalpy of -4.37 kcal/mol at the MP2(full)/6-311G**level of theory. The RHF/6-311G** calculated Theta(zz) value of 1,4-dicyanobenzene is +11.81 DA, yet it has an MP2(full)/6-311G** Li+ binding enthalpy of -12.65 kcal/mol and a Na+ binding enthalpy of -3.72 kcal/mol. The pi-electron density of benzene, the prototypical negative Theta(zz) aromatic (experimental Theta(zz) = -8.7 +/- 0.5 DA), has a F- binding enthalpy of -5.51 kcal/mol. The RHF/6-311G** calculated Theta(zz) of C6H2I4 is -10.45 DA, yet it has an MP2(full)/6-311++G** calculated F- binding enthalpy of -20.13 kcal/mol. Our results show that as the aromatic Theta(zz) value increases the cation binding enthalpy decreases; a plot of cation binding enthalpies versus aromatic Theta(zz) gives a line of best of fit with R2 = 0.778. No such correlation exists between the aromatic Theta(zz) value and the anion binding enthalpy; the line of best fit has R2 = 0.297. These results are discussed in terms of electrostatic and polarizability contributions to the overall binding enthalpies.