Relative stability of alternative chair forms and hydroxymethyl conformations of β- d-glucopyranose

Abstract

The relative energies of two hydroxymethyl conformers for each of the two chair forms ( 4C 1 and 1C 4) of β- d-glucose were calculated at much more complete levels of quantum mechanical (QM) electronic structure theory than previously, 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), and the coupled cluster method for electron correlation. Solvation Model 4 was used to calculate the effects of hydration or nonpolar solvation. Molecular mechanics (MM) and 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 1 C 4 ring conformation to be the lower in energy. The highest level calculations predict the 4 C 1 chair to be lower in free energy by about 8 kcal/mol, and we also find that the gauche + ( gt) conformer of 4 C 1 is lower than the trans (tg) conformer. 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 13 C 4 chair, but are counteracted by increased steric repulsions between axial substituents; MM calculations also indicate increased torsional strain in the 1 C 4 chair. Manifestations of greater steric strain in the calculated 1 C 4 structures compared to the 4 C 1 structures include longer ring bonds, a larger bond angle at the ring oxygen atom, and smaller puckering amplitudes. The MM and QM 4 C 1 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 1 C 4 structures, especially in the relative orientations of hydroxyl groups serving as hydrogen bond acceptors. In water, the 4 C 1 chairs are better solvated than the 4 C 4 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 (a hypothetical) n-hexadecane solution, the 4 C 1 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 4 C 1, conformers allow increased favorable dispersion interactions. The differential polarization free energies are associated primarily with the hydroxyl groups; the greater steric congestion in the 1 C 4 chairs reduces opportunities for favorable dielectric screening.