Journal of Physical Chemistry A, Vol.102, No.26, 5086-5092, 1998
Exocyclic hydroxymethyl rotational conformers of beta- and alpha-D-glucopyranose in the gas phase and aqueous solution
The intrinsic exocyclic hydroxymethyl rotational surface for beta-D-glucoppranose as well as the beta-alpha anomer energy difference for D-glucopyranose has been studied using ab initio quantum mechanical methods including continuum solvation, Relevant stationary points, including rotational transition states, have been characterized by full geometry optimization using the 6-31G(d) basis set for the most stable counterclockwise (cc) overall conformation. Effects of dynamic electron correlation on both the geometric structures and the relative energetics of this system are also explored using Moller-Plesset perturbation theory (MP2 through MP4-(SDTQ)) and density functional methods (BLYP). A total of six stationary points, including three minima and three transition states, were identified along the exocyclic rotational surface. All three minima were found to be very close in energy with a final order of GG (0.0) < CT (2.84) < TG (3.05) based on the relative free energy, Delta G(298)degrees determined at the MP4(SDTQ)/6-31G(d)//MP2/6-31G(d) level of theory. The rotational transition state free energy differences varied from 18.8 to 28.9 kJ mol-l at the same level of theory with the transition state connecting the TG and GG minima being the lowest. The intrinsic gas-phase beta-alpha anomer free energy difference for the cc-TG conformer of D-glucopyranose was also determined at various levels of theory. On the basis of the convergence of the MP series, this energetic quantity has been estimated at 8 +/- 2 kJ mol(-1) favoring the or-anomer, and is insensitive to exocyclic hydroxymethyl rotation. Two different solvation models were used to explore the effects of aqueous solvation on the energetic parameters mentioned above. The Onsager continuum solvation model and the self-consistent isodensity polarized continuum model (SCIPCM) showed significant, yet predictable, effects on the exocyclic hydroxymethyl rotational surface for beta-D-glucopyranose. Shifts in the relative energetics compared to those for the gas phase ranged from -0.8 to + 2.8 kJ mol(-1) for Onsager dipole model and -1.6 to as much as +4.7 kJ mol(-1) for the SCIPCM model at the MP2 6-31G(d) level, resulting in a qualitative change in the ordering of the relative stability of the three stable minima. The effects of the solvation models on the beta-alpha anomer energy difference were also significant, showing a relative decrease in the beta-alpha anomer energy difference from the intrinsic gas-phase result. However, it is clear that these specific continuum solvation models alone cannot account for the experimentally observed preference of the beta-anomer in aqueous solution.
Keywords:MOLECULAR-DYNAMICS SIMULATIONS;FORCE-FIELD;AB-INITIO;CONFORMATIONAL-ANALYSIS;ANOMERIC EQUILIBRIUM;D-GLUCOSE;STEREOCHEMICAL ASPECTS;ORBITAL CALCULATIONS;STRUCTURAL MOIETIES;FREE-ENERGY