Metadynamics modelling of the solvent effect on primary hydroxyl rotamer equilibria in hexopyranosides

Abstract

Accurate modelling of rotamer equilibria for the primary hydroxyl groups of monosaccharides continues to be a great challenge of computational glycochemistry. The metadynamics technique was applied to study the conformational free energy surfaces of methyl α- d-glucopyranoside and methyl α- d-galactopyranoside, employing the glycam06 force field. For both molecules, seven to eight conformational free-energy minima, differing in the ω (O-5–C-5–C-6–O-6) and χ (C-3–C-4–O-4–HO-4) dihedral angles, were identified in vacuum or in a water environment. The calculated rotamer equilibrium of the primary hydroxyl group is significantly different in vacuum than in water. The major effect of a water environment is the destabilisation of a hydrogen bond between O-4–HO-4 and O-6–HO-6 groups. It was possible to calculate the free-energy differences of individual rotamers with an accuracy of better than 2 kJ/mol. The calculated gg, gt and tg rotamer populations in water are in close agreement with experimental measurements, and therefore support the theoretical background of metadynamics.