Low-energy tautomers and conformers of neutral and protonated arginine.

The relative stabilities of zwitterionic and canonical forms of neutral arginine and of its protonated derivative were studied by using ab initio electronic structure methods. Trial structures were first identified at the PM3 level of theory with use of a genetic algorithm to systematically vary geometrical parameters. Further geometry optimizations of these structures were performed at the MP2 and B3LYP levels of theory with basis sets of the 6-31++G** quality. The final energies were determined at the CCSD/6-31++G** level and corrected for thermal effects determined at the B3LYP level. Two new nonzwitterionic structures of the neutral were identified, and one of them is the lowest energy structure found so far. The five lowest energy structures of neutral arginine are all nonzwitterionic in nature and are clustered within a narrow energy range of 2.3 kcal/mol. The lowest energy zwitterion structure is less stable than the lowest nonzwitterion structure by 4.0 kcal/mol. For no level of theory is a zwitterion structure suggested to be the global minimum. The calculated proton affinity of 256.3 kcal/mol and gas-phase basicity of 247.8 kcal/mol of arginine are in reasonable agreement with the measured values of 251.2 and 240.6 kcal/mol, respectively. The calculated vibrational characteristics of the low-energy structures of neutral arginine provide an alternative interpretation of the IR-CRLAS spectrum (Chapo et al. J. Am. Chem. Soc. 1998, 120, 12956-12957).     

Structure of the boronic acid dimer and the relative stabilities of its conformers

Despite the widespread use of boronic acids in materials science and as pharmaceutical agents, many aspects of their structure and reactivity are not well understood. In this research the boronic acid dimer, [HB(OH)(2)](2), was studied by second-order M?ller-Plesset (MP2) perturbation theory and coupled cluster methodology with single and double excitations (CCSD). Pople split-valence 6-31+G*, 6-311G**, and 6-311++G** and Dunning-Woon correlation-consistent cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, and aug-cc-pVTZ basis sets were employed for the calculations. A doubly hydrogen-bonded conformer (1) of the dimer was consistently found to be lowest in energy; the structure of 1 was planar (C(2h)) at most computational levels employed but was significantly nonplanar (C(2)) at the MP2/6-311++G** and CCSD/6-311++G** levels, the result of an intrinsic problem with Pople-type sp-diffuse basis functions on heavy atoms. The dimerization energy, enthalpy, and free energy for the formation of (1) from the exo-endo conformer of the monomer were -10.8, -9.2, and +1.2 kcal/mol, respectively, at the MP2/aug-cc-pVTZ level. Several other hydrogen-bonded conformers of the dimer were local minima on the potential energy surface (PES) and ranged from 2 to 5 kcal/mol higher in energy than 1. Nine doubly OH-bridged conformers, in which the boron atoms were tetracoordinated, were also local minima on the PES, but they were all greater than 13 kcal/mol higher in energy than 1; doubly H-bridged structures proved to be transition states. MP2 and CCSD results were compared to those from the BLYP, B3LYP, OLYP, O3LYP, PBE1PBE, and TPSS functionals with the 6-311++G** and aug-cc-pVTZ basis sets; the PBE1PBE functional performed best relative to the MP2 and CCSD results. Self-consistent reaction field (SCRF) calculations predict that boronic acid dimerization is less favorable in solution than in vacuo.     

Structure of cationized arginine (arg.m, m = h, li, na, k, rb, and cs) in the gas phase: further evidence for zwitterionic arginine

The gas-phase structures of cationized arginine, Arg.M(+), M = Li, Na, K, Rb, and Cs, were studied both by hybrid method density functional theory calculations and experimentally using low-energy collisionally activated and thermal radiative dissociation. Calculations at the B3LYP/LACVP++** level of theory show that the salt-bridge structures in which the arginine is a zwitterion (protonated side chain, deprotonated C-terminus) become more stable than the charge-solvated structures with increasing metal ion size. The difference in energy between the most stable charge-solvated structure and salt-bridge structure of Arg.M(+) increases from -0.7 kcal/mol for Arg.Li(+) to +3.3 kcal/mol for Arg.Cs(+). The stabilities of the salt-bridge and charge-solvated structures reverse between M = Li and Na. These calculations are in good agreement with the results of dissociation experiments. The low-energy dissociation pathways depend on the cation size. Arginine complexed with small cations (Li and Na) loses H(2)O, while arginine complexed with larger cations (K, Rb, and Cs) loses NH(3). Loss of H(2)O must come from a charge-solvated ion, whereas the loss of NH(3) can come from the protonated side chain of a salt-bridge structure. The results of dissociation experiments using several cationized arginine derivatives are consistent with the existence of these two distinct structures. In particular, arginine methyl esters, which cannot form salt bridges, dissociate by loss of methanol, analogous to loss of H(2)O from Arg.M(+); no loss of NH(3) is observed. Although dissociation experiments probe gas-phase structure indirectly, the observed fragmentation pathways are in good agreement with the calculated lowest energy isomers. The combination of the results from experiment and theory provides strong evidence that the structure of arginine-alkali metal ion complexes in the gas phase changes from a charge-solvated structure to a salt-bridge structure as the size of the metal ion increases.     

A quantum-chemical study of the LiH + LiH⋅− ⇄ LiH2− + Li⋅ reaction

Parts of the potential energy surface of the title process and related processes have been investigated at the SCF /6-31G **, SCF /6-31++G **, and MP 2/6-31++G ** levels. The investigated reaction is exothermic (?6.23 kcal/mol, MP 4/6-31++G **//MP 2/6–31++G** level, ZPE included): A linear intermediate radical anion, Li? H? Li? H??, is significantly stabilized with respect to LiH + LiH?? (?38.74 kcal/mol, the same level as above). The BSSE at MP 2/6–31++G **//MP 2/6–31++G ** amounts to 1.8 kcal/mol. The title process seems to be suitable for experimental study in molecular beams.     

Conformational energies for 2-substituted butanes

The conformational free energies for some 2-substituted butanes where X = F, Cl, CN, and CCH were calculated using G3-B3, CBS-QB3, and CCSD(T)/6-311++G(2d,p) as well as other theoretical levels. The above methods gave consistent results with free energies relative to the trans conformers as follows: X = CCH, g+ = 0.77 +/- 0.05 kcal/mol. g- = 0.88 +/- 0.05 kcal/mol; X = CN, g+ = 0.85 +/- 0.05 kcal/mol, g- = 0.75 +/- 0.05 kcal/mol; X = Cl, g+ = 0.70 +/- 0.05 kcal/ml, g- = 0.80 +/- 0.05 kcal/mol; and X = F, g+ = 0.53 +/- 0.05 kcal/mol, g- = 0.83 +/- 0.05 kcal/mol. The conformational free energies also were estimated using the observed liquid phase IR spectra and intensities calculated using B3LYP/6-311++G** and MP2/6-311++G**. The rotational free energy profiles for all of the compounds were estimated at the G3-B3 level.     

Interaction of anions with perfluoro aromatic compounds

The complexes formed by a variety of anions with perfluoro derivatives of benzene, naphthalene, pyridine, thiophene, and furan have been calculated using DFT (B3LYP/6-31++G**) and MP2 (MP2/6-31++G** and MP2/6-311++G**) ab initio methods. The minimum structures show the anion interacting with the pi-cloud of the aromatic compounds. The interaction energies obtained range between -8 and -19 kcal mol(-1). The results obtained at the MP2/6-31++G** and MP2/6-311++G** levels are similar. However, the B3LYP/6-31++G** results provide longer interaction distances and smaller interaction energies than do the MP2 results. The interaction energies have been partitioned using an electrostatic, polarization, and van der Waals scheme. The AIM analysis of the electron density shows a variety of topologies depending on the aromatic system considered.     

The C7-C10 cycloalkanes revisited

The conformers of cycloheptane through cyclodecane have been examined at the B3LYP/6-311+G* and MP2/6-311+G* theoretical levels, with some additional calculations at the CCD/6-311+G* and CCSD(T)/6-311++G** levels. With cyclooctane, B3LYP predicts that the boat-chair and crown conformers have similar energies, whereas MP2 and CCSD(T) predict that the crown conformer is 2 kcal/mol higher in energy. The latter is in agreement with the electron diffraction data. With cyclononane, B3LYP predicts that two of the higher-energy conformers found in molecular mechanics calculations should convert to one of the lower-energy conformers. However, MP2/6-311+G* optimizations find them to be true minima on the potential energy surface. B3LYP systematically predicts larger C-C-C bond angles for these compounds than either MP2 or CCD. The results of molecular mechanics MM4 calculations are generally in good agreement with those obtained using MP2.     

A semi-empirical and ab initio combined approach for the full conformational searches of gaseous lysine and lysine–H2O complex

A full structural search of the canonical, zwitterionic, protonated and deprotonated lysine conformers in gas phase is presented. A total of 17,496 canonical, 972 zwitterionic, 11,664 protonated and 1458 trial deprotonated structures were generated by allowing for all combinations of internal single-bond rotamers. All the trial structures were initially optimized at the AM1 level, and the resulting structures were determined at the B3LYP/6-311G^* level. A total of 927 canonical, 730 protonated and 193 deprotonated conformers were found, but there were no stable zwitterionic structures in the gas phase. The most stable conformers of the canonical, protonated and deprotonated lysine were further optimized at the B3LYP/6-311++G^** level. The energies of the most stable structures were determined at the MP2/6-311G(2df,p) level and the vibrational frequencies were calculated at the B3LYP/6-311++G^** level. The rotational constants, dipole moments, zero-point vibrational energies, harmonic frequencies, vertical ionization energies, enthalpies, Gibbs free energies and conformational distributions of gaseous lysine were presented. Numerous new structures are found and the lowest-energy lysine conformer is more stable than the existing one by 1.1 kcal/mol. Hydrogen bonds are classified and may cause significant red-shifts to the associated vibrational frequencies. The calculated proton affinity/dissociation energy and gas-phase basicity/acidity are in good agreement with the experiments. Calculations are also presented for the canonical lysine–H₂O and zwitterionic lysine–H₂O clusters. Interaction between lysine and H₂O significantly affects the relative conformational stabilities. Only one water molecule is sufficient to produce the stable zwitterionic structures in gas phase. The lowest-energy structure is found to be zwitterions when applying the conductor-like polarized continuum solvent model (CPCM) to the lysine–H₂O complexes.     

Dihydrogen trioxide clusters, (HOOOH)n (n = 2-4), and the hydrogen-bonded complexes of HOOOH with acetone and dimethyl ether: implications for the decomposition of HOOOH

Hydrogen-bonded gas-phase molecular clusters of dihydrogen trioxide (HOOOH) have been investigated using DFT (B3LYP/6-311++G(3df,3pd)) and MP2/6-311++G(3df,3pd) methods. The binding energies, vibrational frequencies, and dipole moments for the various dimer, trimer, and tetramer structures, in which HOOOH acts as a proton donor as well as an acceptor, are reported. The stronger binding interaction in the HOOOH dimer, as compared to that in the analogous cyclic structure of the HOOH dimer, indicates that dihydrogen trioxide is a stronger acid than hydrogen peroxide. A new decomposition pathway for HOOOH was explored. Decomposition occurs via an eight-membered ring transition state for the intermolecular (slightly asynchronous) transfer of two protons between the HOOOH molecules, which form a cyclic dimer, to produce water and singlet oxygen (Delta (1)O 2). This autocatalytic decomposition appears to explain a relatively fast decomposition (Delta H a(298K) = 19.9 kcal/mol, B3LYP/6-311+G(d,p)) of HOOOH in nonpolar (inert) solvents, which might even compete with the water-assisted decomposition of this simplest of polyoxides (Delta H a(298K) = 18.8 kcal/mol for (H 2O) 2-assisted decomposition) in more polar solvents. The formation of relatively strongly hydrogen-bonded complexes between HOOOH and organic oxygen bases, HOOOH-B (B = acetone and dimethyl ether), strongly retards the decomposition in these bases as solvents, most likely by preventing such a proton transfer.     

Quantum-Chemical Study of Alkyl Carbenium Ions in 100% Sulfuric Acid

A quantum-chemical study of neutral and protonated monoalkyl sulfates RHSO₄and [RH₂SO₄]⁺(where R = CH₃, C₂H₅, iso-C₃H₇, and tert-C₄H₉) is carried out. Calculations are performed using the Hartree–Fock method in the 6-31G** and 6-31++G** basis sets taking into account electron correlation according to the Müller–Plesset perturbation theory MP2/6-31+G*//6-31+G*. Protonated tert-butyl sulfate was also calculated by the DFT B3LYP/6-31++G** method. It was found that monoalkyl sulfates are covalent compounds, and the complete abstraction of alkyl carbenium ions from them has substantial energy cost: 196.4, 161.7, 150.8 and 136.0 kcal/mol, respectively. Protonated methyl and ethyl sulfates are also covalent compounds according to the calculation. They have lower but still high energies of heterolytic dissociation (65.0 and 33.5 kcal/mol, respectively). The energy of R⁺abstraction from protonated isopropyl sulfate is much lower: 23.6 kcal/mol. The main covalent state and the ion–molecular pair, which is a carbenium ion [C(CH₃)₂H]⁺solvated by the H₂SO₄molecule, have about the same energy. The ground state of protonated tert-butyl sulfate corresponds to the ion–molecular complex [C(CH₃)⁺ ₃H₂SO₄] with still lower energy of carbenium ion [C(CH₃)₃]⁺abstraction, which is equal to 10.0 kcal/mol. Calculation according to the DFT B3LYP/6-31++G** method shows the absence of a minimum for the protonated tert-butyl sulfate with a covalent structure on the potential energy surface.     

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.     

气相中O3与HSO自由基间的氢键复合物

气相中O3与HSO自由基之间的相互作用及其反应在大气化学中非常重要.在DFT-B3LYP/6-311++G**和MP2/6-311++G**水平上求得O3+HSO复合物势能面上的稳定构型,B3LYP方法得到了三种构型(复合物Ⅰ,Ⅱ和Ⅲ),而MP2方法只能得到一种构犁(复合物Ⅱ).在复合物Ⅰ和Ⅲ中,HSO单元中的1H原子作为质子供体.与O3分子中的端基O原子作为质子受体相互作用,形成红移氢键复合物;而在复合物Ⅱ中,虽与复合物Ⅰ和Ⅲ中具有相间的质子供体和质子受体,却形成了蓝移氢键复合物.B3LYP/6-311++G**水平上计算的单体间相互作用能的计算考虑了基组重甍误差(BSSE)和零点振动能(ZPVE)校正,其值在-3.37到-4.55 kJ·mol-1之间.采用自然键轨道理论(NBO)对单体间相互作用的本质进行了考查,并通过分子中原子理论(AIM)分析了三种复合物中氢键的电子密度拓扑性质.     

An ab initio study of intermolecular interactions of nitromethane dimer and nitromethane trimer

Different geometries of nitromethane dimer and nitromethane trimer have been fully optimized employing the density functional theory B3LYP method and the 6-31++G** basis set. Three-body interaction energy has been obtained with the ab initio supermolecular approach at the levels of MP2/6-31++G**//B3LYP/6-31++G** and MP2/aug-cc-pVDZ//B3LYP/6-31++G**. The internal rotation of methyl group induced by intermolecular interaction has been observed theoretically. For the optimized structures of nitromethane dimer, the strength of C--H...O--N H-bond ranges from -9.0 to -12.4 kJ mol(-1) at the MP2/aug-cc-pVDZ//B3LYP/6-31++G** level, and the B3LYP method underestimates the interaction strength compared with the MP2 method, while MP2/6-31++G**//B3LYP/6-31++G** calculated DeltaE(C) is within 2.5 kJ mol(-1) of the corresponding value at the MP4(SDTQ)/6-31G**//B3LYP/6-31++G** level. The analytic atom-atom intermolecular potential has been successfully regressed by using the MP2/6-31++G**//B3LYP/6-31++G** calculated interaction energies of nitromethane dimer. For the optimized structures of nitromethane trimer the three-body interaction energies occupy small percentage of corresponding total binding energies, but become important for the compressed nitromethane explosive. In addition, it has been discovered that the three-body interaction energy in the cyclic nitromethane trimer is more and more negative as intermolecular distances decrease from 2.2 to 1.7 A.     

CH···O interaction lowers hydrogen transfer barrier to keto-enol tautomerization of β-cyclohexanedione: combined infrared spectroscopic and electronic structure calculation study

Molecular association and keto-enol tautomerization of β-cyclohexanedione (β-CHD) have been investigated in argon matrix and also in a thin solid film prepared by depositing pure β-CHD vapor on a cold (8 K) KBr window. Infrared spectra reveal that, in low-pressure vapor and argon matrix, the molecules are exclusively in diketo tautomeric form. The CH···O hydrogen bonded dimers of the diketo tautomer are produced by annealing the matrix at 28 K. No indication is found for keto-enol tautomerization of β-CHD in dimeric complexes in argon matrix within the temperature range of 8-28 K. On the other hand, in thin film of pure diketo tautomer, the conversion initiates only when the film is heated at temperatures above 165 K. The observed threshold appears to be associated with excitation of the intermolecular modes, and the IR spectra recorded at high temperatures display narrowing of vibrational bandwidths, which has been associated with reorientations of the molecules in the film. The nonoccurrence of tautomerization of the matrix isolated dimer is consistent with the barrier predicted by electronic structure calculations at B3LYP/6-311++G** and MP2/6-311++G** levels of theory. The transition state calculation predicts that CH···O interaction has a dramatic effect on lowering of the tautomerization barrier, from more than 60 kcal/mol for the bare molecule to ~35-45 kcal/mol for dimers.     

Pi-pi interaction in pyridine

pi-pi Interaction in pyridine dimer and trimer has been investigated in different geometries and orientations at the ab initio (HF, MP2) and DFT (B3LYP) levels of theory using various basis sets (6-31G, 6-31G, 6-311++G) and corrected for basis set superposition error (BSSE). While the HF and DFT calculations show the pyridine dimer and the trimer to be unstable with respect to the monomer, the MP2 calculations show them to be clearly stable, thus emphasizing the need to include electron correlation while determining stacking interaction in such systems. The calculated MP2/6-311++G binding energy (100% BSSE corrected) of the parallel-sandwich, antiparallel-sandwich, parallel-displaced, antiparallel-displaced, T-up and T-down geometries for pyridine dimer are 1.53, 3.05, 2.39, 3.97, 1.91, 1.47 kcal/mol, respectively. The results show the antiparallel-displaced geometry to be the most stable. The binding energies for the trimer in parallel-sandwich, antiparallel-sandwich, and antiparallel-displaced geometry are found to be 3.18, 6.14, and 8.04 kcal/mol, respectively.     

One water molecule stabilizes the cationized arginine zwitterion

Singly hydrated clusters of lithiated arginine, sodiated arginine, and lithiated arginine methyl ester are investigated using infrared action spectroscopy and computational chemistry. Whereas unsolvated lithiated arginine is nonzwitterionic, these results provide compelling evidence that attachment of a single water molecule to this ion makes the zwitterionic form of arginine, in which the side chain is protonated, more stable. The experimental spectra of lithiated and sodiated arginine with one water molecule are very similar and contain spectral signatures for protonated side chains, whereas those of lithiated arginine and singly hydrated lithiated arginine methyl ester are different and contain spectral signatures for neutral side chains. Calculations at the B3LYP/6-31++G** level of theory indicate that solvating lithiated arginine with a single water molecule preferentially stabilizes the zwitterionic forms of this ion by 25-32 kJ/mol and two essentially isoenergetic zwitterionic structure are most stable. In these structures, the metal ion either coordinates with the N-terminal amino group and an oxygen atom of the carboxylate group (NO coordinated) or with both oxygen atoms of the carboxylate group (OO coordinated). In contrast, the OO-coordinated zwitterionic structure of sodiated arginine, both with and without a water molecule, is clearly lowest in energy for both ions. Hydration of the metal ion in these clusters weakens the interactions between the metal ion and the amino acid, whereas hydrogen-bond strengths are largely unaffected. Thus, hydration preferentially stabilizes the zwitterionic structures, all of which contain strong hydrogen bonds. Metal ion size strongly affects the relative propensity for these ions to form NO or OO coordinated structures and results in different zwitterionic structures for lithiated and sodiated arginine clusters containing one water molecule.     

Structural studies on ethyl isovalerate by microwave spectroscopy and quantum chemical calculations

We observed the microwave spectrum of ethyl isovalerate by molecular beam Fourier transform microwave spectroscopy. The rotational and centrifugal distortion constants of the most abundant conformer were determined. Its structure was investigated by comparison of the experimental rotational constants with those obtained by ab initio methods. In a first step, the rotational constants of various conformers were calculated at the MP2/6-311++G** level of theory. Surprisingly, no agreement with the experimental results was found. Therefore, we concluded that in the case of ethyl isovalerate more advanced quantum chemical methods are required to obtain a reliable molecular geometry. Ab initio calculations carried out at MP3/6-311++G**, MP4/6-311++G**, and CCSD/6-311++G** levels and also density functional theory calculations using the B3LYP/6-311++G** method gave similar results for the rotational constants, but they were clearly distinct from those obtained at the MP2/6-311++G** level. With use of these more advanced methods, the rotational constants of the lowest energy conformer were in good agreement with those obtained from the microwave spectrum.     

Theoretical Studies on the Intermolecular Interactions of Aza-calix[2]arene[2]-triazines with RDX

Six fully optimized structures of the aza-calix[2]arene[2]-triazines/RDX supramo-lecular complexes have been obtained at the DFT-B3LYP/6-311++G** level,and the corresponding intermolecular interactions have been investigated using the B3LYP,mPWPW91 and MP2 methods at the 6-311++G** level,respectively.The natural bond orbital(NBO) and atoms in molecules(AIM) analyses have been performed to reveal the origin of interactions.To our interest,the result indicates that the strongest interaction is up to -22.34 kJ/mol after basis set superposition error(BSSE) and zero point energy(ZPE) correction at the MP2/6-311++G** level.Furthermore,the intermolecular interactions between aza-calix[2]arene[2]-triazines with the substituted amidos and RDX are stronger than those of other complexes.Thus,the complexes with amidos can be used as the candidates to increase the stability of explosive and eliminate the explosive wastewater.     

On the question of salt bridges of cationized amino acids in the gas phase: glycine and arginine

The geometrical shapes of the sodiated and cesiated amino acids glycine and arginine were probed in the gas phase by using the ion mobility based ion chromatography method. The data were compared to those obtained for alkali cationized methyl esters and for all the protonated species. Molecular mechanics, semiempirical, and ab initio/density functional theory (DFT) calculations were carried out to generate model structures for comparison with experiment and to determine the relative energies of different structures. For alkali cationized glycine the experimental cross sections agreed with charge solvation structures which were found by calculation to be the most stable forms as well. Both experiment and theory indicated that sodium is solvated by both the amino and the carbonyl groups, while cesium is solvated by one or both oxygen(s) of the carboxyl group. Alkali cationized arginine was found to form a salt bridge structure. The carboxylate group is stabilized by both the charged guanidinium group and the alkali ion. High level (6-311++G7 and DZVP) ab initio/DFT calculations carried out on sodiated and rubidiated N amidinoglycine, which contains a guanidino group and which was used as a model for the larger arginine molecule, indicated that the salt bridge structures are ∼10 kcal/mol more stable than the charge solvation forms for both alkali ions. The structure of protonated arginine, i.e. salt bridge or charge solvation, could not be unambiguously determined.     

Computational studies on carbohydrates: I. Density functional ab initio geometry optimization on maltose conformations

Ab initio geometry optimization was carried out on 10 selected conformations of maltose and two 2‐methoxytetrahydropyran conformations using the density functional denoted B3LYP combined with two basis sets. The 6‐31G* and 6‐311++G** basis sets make up the B3LYP/6‐31G* and B3LYP/6‐311++G** procedures. Internal coordinates were fully relaxed, and structures were gradient optimized at both levels of theory. Ten conformations were studied at the B3LYP/6‐31G* level, and five of these were continued with full gradient optimization at the B3LYP/6‐311++G** level of theory. The details of the ab initio optimized geometries are presented here, with particular attention given to the positions of the atoms around the anomeric center and the effect of the particular anomer and hydrogen bonding pattern on the maltose ring structures and relative conformational energies. The size and complexity of the hydrogen‐bonding network prevented a rigorous search of conformational space by ab initio calculations. However, using empirical force fields, low‐energy conformers of maltose were found that were subsequently gradient optimized at the two ab initio levels of theory. Three classes of conformations were studied, as defined by the clockwise or counterclockwise direction of the hydroxyl groups, or a flipped conformer in which the ψ‐dihedral is rotated by ∼180°. Different combinations of ω side‐chain rotations gave energy differences of more than 6 kcal/mol above the lowest energy structure found. The lowest energy structures bear remarkably close resemblance to the neutron and X‐ray diffraction crystal structures. © 2000 John Wiley & Sons, Inc. * J Comput Chem 21: 1204–1219, 2000