Conformational properties of glucose-based disaccharides investigated using molecular dynamics simulations with local elevation umbrella sampling

Conformational properties of methyl β-maltoside and methyl α- and β-cellobioside disaccharides

An investigation of the conformational properties of methyl β-maltoside, methyl α-cellobioside, and methyl β-cellobioside disaccharides using NMR spectroscopy and molecular dynamics (MD) techniques, is presented. Emphasis is placed on validation of a recently presented force field for hexopyranose disaccharides followed by elucidation of the conformational properties of two different types of glycosidic linkages, α-(1 → 4) and β-(1 → 4). Both gas-phase and aqueous-phase simulations are performed to gain insight into the effect of solvent on the conformational properties. A number of transglycosidic J-coupling constants and proton-proton distances are calculated from the simulations and are used to identify the percent sampling of the three glycosidic conformations (syn, anti-φ, and anti-ψ) and, in turn, describe the flexibility around the glycosidic linkage. The results show the force field to be in overall good agreement with experiment, although some very small limitations are evident. Subsequently, a thorough hydrogen bonding analysis is performed to obtain insights into the conformational properties of the disaccharides. In methyl β-maltoside, competition between HO2'-O3 intramolecular hydrogen bonding and intermolecular hydrogen bonding of those groups with solvent leads to increased sampling of syn, anti-φ, and anti-ψ conformations and better agreement with NMR J-coupling constants. In methyl α- and β-cellobioside, O5'-HO6 and HO2'-O3 hydrogen bonding interactions are in competition with intermolecular hydrogen bonding involving the solvent molecules. This competition leads to retention of the O5'-HO3 hydrogen bond and increased sampling of the syn region of the φ/ψ map. Moreover, glycosidic torsions are correlated to the intramolecular hydrogen bonding occurring in the molecules. The present results verify that in the β-(1 → 4)-linkage intramolecular hydrogen bonding in the aqueous phase is due to the decreased ability of water to successfully compete for the O5' and HO3 hydrogen bonding moieties, in contrast to that occurring between the O5' and HO6 atoms in this α-(1 → 4)-linkage.

A reoptimized GROMOS force field for hexopyranose-based carbohydrates accounting for the relative free energies of ring conformers, anomers, epimers, hydroxymethyl rotamers, and glycosidic linkage conformers

This article presents a reoptimization of the GROMOS 53A6 force field for hexopyranose-based carbohydrates (nearly equivalent to 45A4 for pure carbohydrate systems) into a new version 56A(CARBO) (nearly equivalent to 53A6 for non-carbohydrate systems). This reoptimization was found necessary to repair a number of shortcomings of the 53A6 (45A4) parameter set and to extend the scope of the force field to properties that had not been included previously into the parameterization procedure. The new 56A(CARBO) force field is characterized by: (i) the formulation of systematic build-up rules for the automatic generation of force-field topologies over a large class of compounds including (but not restricted to) unfunctionalized polyhexopyranoses with arbritrary connectivities; (ii) the systematic use of enhanced sampling methods for inclusion of experimental thermodynamic data concerning slow or unphysical processes into the parameterization procedure; and (iii) an extensive validation against available experimental data in solution and, to a limited extent, theoretical (quantum-mechanical) data in the gas phase. At present, the 56A(CARBO) force field is restricted to compounds of the elements C, O, and H presenting single bonds only, no oxygen functions other than alcohol, ether, hemiacetal, or acetal, and no cyclic segments other than six-membered rings (separated by at least one intermediate atom). After calibration, this force field is shown to reproduce well the relative free energies of ring conformers, anomers, epimers, hydroxymethyl rotamers, and glycosidic linkage conformers. As a result, the 56A(CARBO) force field should be suitable for: (i) the characterization of the dynamics of pyranose ring conformational transitions (in simulations on the microsecond timescale); (ii) the investigation of systems where alternative ring conformations become significantly populated; (iii) the investigation of anomerization or epimerization in terms of free-energy differences; and (iv) the design of simulation approaches accelerating the anomerization process along an unphysical pathway.