Relative Stability of Alternative Chair Forms and Hydroxymethyl Conformations of D-Glucopyranose

Barrows, S. E.; Dulles, F. J.; Cramer, C. J.; Truhlar, D. G.; French, A. D.
Carbohydr. Res. 1995, 276, 219.

The relative energies of two hydroxymethyl conformers for each of the two chair forms (⁴C1 and ¹C4) of [[beta]]-D-glucose 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 vibrational, rotational, and solvent effects. The gas-phase results are based on very large basis sets (up to 624 contracted basis functions), the coupled cluster method for electron correlation, and Solvation Model 4 for the effects of hydration or nonpolar solvation. Molecular mechanics (MM) and quantum mechanical (QM) electronic structure theory have been applied to analyze the factors contributing to the relative energies of these conformers. Relative energies varied widely (up to 35 kcal/mol) depending on theoretical level, and several levels of theory predict the experimentally unobserved ¹C4 ring conformation to be the lower in energy. The highest-level calculations predict the ⁴C1 chair to be lower in free energy by about 8 kcal/mol, and we also find that the gauche⁺, i.e., gt, conformer of ⁴C1 is lower than the trans, i.e., tg. Low-energy structures optimized by either quantum mechanical or molecular mechanical methods were commonly characterized by multiple intramolecular hydrogen bonds. Superior hydrogen bonding geometries are available in the ¹C4 chair, but are counteracted by increased steric repulsions between axial substituents; MM calculations also indicate increased torsional strain in the ¹C4 chair. Manifestations of greater steric strain in the calculated ¹C4 structures compared to the ⁴C1 structures include longer ring bonds, a larger bond angle at the ring oxygen atom, and smaller puckering amplitudes. The MM and QM ⁴C1 structures compare well with each other and with available X-ray diffraction data. The largest discrepancies between the two kinds of models occur for geometric parameters associated with the anomeric center--the QM structure agrees better with experiment. Greater differences between QM and MM structures are observed for ¹C4 structures, especially in the relative orientations of hydroxyl groups serving as hydrogen bond acceptors. In water, the ⁴C1 chairs are better solvated than the ¹C4 chairs by about 5 to 9 kcal/mol because of both larger polarization free energies and improved hydrogen bonding interactions with the first solvation shell. In n-hexadecane, the ⁴C1 chairs are better solvated by about 2 to 4 kcal/mol both because of larger polarization free energies and because the larger solvent accessible surface areas of the ⁴C1 conformers allow increased favorable dispersion interactions. The differential polarization free energies are associated primarily with the hydroxyl groups; the greater steric congestion in the ¹C4 chairs reduces opportunities for favorable dielectric screening.

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 95/39.