Ab initio modeling of glycosyl torsions and anomeric effects in a model carbohydrate: 2-ethoxy tetrahydropyran.

A range of ab initio calculations were carried out on the axial and equatorial anomers of the model carbohydrate 2-ethoxy tetrahydropyran to evaluate the level of theory required to accurately evaluate the glycosyl dihedral angle and the anomeric ratio. Vacuum CCSD(T)/CBS extrapolations at the global minimum yield DeltaE = E(equatorial) - E(axial) = 1.42 kcal/mol. When corrected for solvent (by the IEFPCM model), zero-point vibrations and entropy, DeltaG(298) = 0.49 kcal/mol, in excellent agreement with the experimental value of 0.47 +/- 0.3 kcal/mol. A new additivity scheme, the layered composite method (LCM), yields DeltaE to within 0.1 kcal/mol of the CCSD(T)/CBS result at a fraction of the computer requirements. Anomeric ratios and one-dimensional torsional surfaces generated by LCM and the even more efficient MP2/cc-pVTZ level of theory are in excellent agreement, indicating that the latter is suitable for force-field parameterization of carbohydrates. Hartree-Fock and density functional theory differ from CCSD(T)/CBS for DeltaE by approximately 1 kcal/mol; they show similar deviations in torsional surfaces evaluated from LCM. A comparison of vacuum and solvent-corrected one- and two-dimensional torsional surfaces indicates the equatorial form of 2-ethoxy tetrahydropyran is more sensitive to solvent than the axial.

[1]  Donald G. Truhlar,et al.  Exo-anomeric effects on energies and geometries of different conformations of glucose and related systems in the gas phase and aqueous solution , 1997 .

[2]  J. Praly,et al.  Influence of solvent on the magnitude of the anomeric effect , 1987 .

[3]  Feliu Maseras,et al.  IMOMM: A new integrated ab initio + molecular mechanics geometry optimization scheme of equilibrium structures and transition states , 1995, J. Comput. Chem..

[4]  L. Radom,et al.  G3-RAD and G3X-RAD: Modified Gaussian-3 (G3) and Gaussian-3X (G3X) procedures for radical thermochemistry , 2003 .

[5]  R. Lemieux,et al.  Effects of unshared pairs of electrons and their solvation on conformational equilibria , 1971 .

[6]  Donald Voet,et al.  Biochemistry, 2nd ed. , 1995 .

[7]  Carlos A. Stortz,et al.  Comparative performance of MM3(92) and two TINKER™ MM3 versions for the modeling of carbohydrates , 2005, J. Comput. Chem..

[8]  Shawn T. Brown,et al.  Advances in methods and algorithms in a modern quantum chemistry program package. , 2006, Physical chemistry chemical physics : PCCP.

[9]  O. Runquist,et al.  Conformational analysis of some 2-alkoxytetrahydropyrans , 1968 .

[10]  Jacopo Tomasi,et al.  Continuum solvation models: A new approach to the problem of solute’s charge distribution and cavity boundaries , 1997 .

[11]  A. J. Kirby,et al.  The Anomeric Effect and Related Stereoelectronic Effects at Oxygen , 1982 .

[12]  Jacopo Tomasi,et al.  A new integral equation formalism for the polarizable continuum model: Theoretical background and applications to isotropic and anisotropic dielectrics , 1997 .

[13]  David Feller,et al.  The use of systematic sequences of wave functions for estimating the complete basis set, full configuration interaction limit in water , 1993 .

[14]  Frank A. Momany,et al.  Computational studies on carbohydrates: I. Density functional ab initio geometry optimization on maltose conformations , 2000 .

[15]  Rengaswami Chandrasekaran,et al.  Conformation of Carbohydrates , 1998 .

[16]  Nathan J DeYonker,et al.  The correlation consistent composite approach (ccCA): an alternative to the Gaussian-n methods. , 2006, The Journal of chemical physics.

[17]  Kjeld Rasmussen,et al.  A comparison and chemometric analysis of several molecular mechanics force fields and parameter sets applied to carbohydrates , 1998 .

[18]  J. Brady,et al.  The role of hydrogen bonding in carbohydrates: molecular dynamics simulations of maltose in aqueous solution , 1993 .

[19]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[20]  U. Salzner,et al.  Ab Initio Examination of Anomeric Effects in Tetrahydropyrans, 1,3-Dioxanes, and Glucose , 1994 .

[21]  I. Tvaroška,et al.  The anomeric and exo-anomeric effects of a hydroxyl group and the stereochemistry of the hemiacetal linkage , 1998 .

[22]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[23]  L. Curtiss,et al.  Gaussian-3 (G3) theory for molecules containing first and second-row atoms , 1998 .

[24]  T. H. Dunning.  Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[25]  Hanoch Senderowitz,et al.  Carbohydrates: United Atom AMBER* Parameterization of Pyranoses and Simulations Yielding Anomeric Free Energies , 1996 .

[26]  K. Wiberg,et al.  Solvent Effects. 5. Influence of Cavity Shape, Truncation of Electrostatics, and Electron Correlation on ab Initio Reaction Field Calculations , 1996 .

[27]  Trygve Helgaker,et al.  Basis-set convergence of correlated calculations on water , 1997 .

[28]  Robert J. Woods,et al.  Molecular Mechanical and Molecular Dynamic Simulations of Glycoproteins and Oligosaccharides. 1. GLYCAM_93 Parameter Development , 1995 .

[29]  J. Tomasi,et al.  The IEF version of the PCM solvation method: an overview of a new method addressed to study molecular solutes at the QM ab initio level , 1999 .

[30]  I. Tvaroška,et al.  Ab initio molecular orbital calculation of carbohydrate model compounds. 2. Conformational analysis of axial and equatorial 2-methoxytetrahydropyrans , 1994 .

[31]  W. Szarek,et al.  Anomeric Effect: Origin and Consequences , 1979 .

[32]  J. C. Martin,et al.  Solvation effects on conformational equilibria. Studies related to the conformational properties of 2-methoxytetrahydropyran and related methyl glycopyranosides , 1969 .

[33]  Gregory S. Tschumper,et al.  Intrinsic conformational preferences of substituted cyclohexanes and tetrahydropyrans evaluated at the CCSD(T) complete basis set limit: implications for the anomeric effect. , 2005, The journal of physical chemistry. A.

[34]  Jacopo Tomasi,et al.  Evaluation of Solvent Effects in Isotropic and Anisotropic Dielectrics and in Ionic Solutions with a Unified Integral Equation Method: Theoretical Bases, Computational Implementation, and Numerical Applications , 1997 .

[35]  E. Eliel.  Conformational Analysis in Heterocyclic Systems: Recent Results and Applications , 1972 .

[36]  H. Booth,et al.  Experimental studies of the anomeric effect , 1987 .

[37]  David Feller,et al.  Application of systematic sequences of wave functions to the water dimer , 1992 .

[38]  L. Radom,et al.  An Assessment of the Performance of High-Level Theoretical Procedures in the Computation of the Heats of Formation of Small Open-Shell Molecules , 2002 .