Factors Controlling the Relative Stability of Anomers and Hydroxymethyl Conformers of Glucopyranose

Barrows, S. E.; Storer, J. W.; Cramer, C. J.; French, A. D.; Truhlar, D. G.
J. Comput. Chem. 1998, 19, 1111.

The relative energies of eleven different conformers of D-glucose, including different exo-anomeric orientations and at least one of each hydroxymethyl conformer (G^-, G⁺, and T) for each of the two anomeric forms (a and b), were calculated at much more complete levels of quantum mechanical (QM) electronic structure theory than previously available, and relative free energies in solution were calculated by adding rotational, vibrational, and aqueous solvation effects. The gas-phase results are based on very large basis sets (up to 624 contracted basis functions) and the coupled cluster method for electron correlation. Solvation Model 5.4/AM1 was used to calculate the effects of aqueous solvation. Factors contributing to the relative energies of these conformers have been analyzed. Relative energies vary considerably (up to 4.5 kcal/mol) depending on the theoretical level, and different levels of theory disagree as to which anomer is the lower in energy. The highest-level gas-phase calculations predict the a anomer to be lower in free energy by 0.4 kcal/mol (Boltzmann average). Gas-phase energies from several different force fields are compared to the QM results. The QM structures optimized at the MP2/cc-pVDZ level of theory compare well with experiment for three different crystal structures. In water, the b anomers are better solvated than the a anomers by 0.6 kcal/mol (Boltzmann average). Contributions of individual hydrophilic groups to the solvation free energies are analyzed.

To request a copy of this article, send e-mail to the Research Reports Coordinator at the Minnesota Supercomputer Institute (requests@msi.umn.edu). Please provide a mailing address and specify that you would like UMSI report 97/153.