Conformational properties of α- or β-(1→6)-linked oligosaccharides: Hamiltonian replica exchange MD simulations and...

Conformational sampling for a set of 10 α- or β-(1→6)-linked oligosaccharides has been studied using explicit solvent Hamiltonian replica exchange (HREX) simulations and NMR spectroscopy techniques. Validation of the force field and simulation methodology is done by comparing calculated transglycosidic J coupling constants and proton-proton distances with the corresponding NMR data. Initial calculations showed poor agreement, for example, with >3 Hz deviation of the calculated (3)J(H5,H6R) values from the experimental data, prompting optimization of the ω torsion angle parameters associated with (1→6)-linkages. The resulting force field is in overall good agreement (i.e., within ∼0.5 Hz deviation) from experimental (3)J(H5,H6R) values, although some small limitations are evident. Detailed hydrogen bonding analysis indicates that most of the compounds lack direct intramolecular H-bonds between the two monosaccharides; however, minor sampling of the O6···HO2' hydrogen bond is present in three compounds. The results verify the role of the gauche effect between O5 and O6 atoms in gluco- and manno-configured pyranosides causing the ω torsion angle to sample an equilibrium between the gt and gg rotamers. Conversely, galacto-configured pyranosides sample a population distribution in equilibrium between gt and tg rotamers, while the gg rotamer populations are minor. Water radial distribution functions suggest decreased accessibility to the O6 atom in the (1→6)-linkage as compared to the O6' atom in the nonreducing sugar. The role of bridging water molecules between two sugar moieties on the distributions of ω torsion angles in oligosaccharides is also explored.