A combined nuclear magnetic resonance and molecular dynamics study of the two structural motifs for mixed-linkage β-glucans: methyl β-cellobioside and methyl β-laminarabioside Original Research Article

The conformational and hydration properties of the two disaccharides methyl β-cellobioside and methyl β-laminarabioside were investigated by NMR spectroscopy and explicit solvation molecular dynamics simulations using the carbohydrate solution force field (CSFF). Adiabatic maps produced with this force field displayed 4 minima A: (Φ = 300°, Ψ = 280°), B: (Φ = 280°, Ψ = 210°), C: (Φ = 260°, Ψ = 60°), and D: (Φ = 60°, Ψ = 260°) for methyl β-cellobioside and 3 minima A: (Φ = 290°, Ψ = 130°), B: (Φ = 270°, Ψ = 290°), and C: (Φ = 60°, Ψ = 120°) for methyl β-laminarabioside. Molecular dynamics simulations were initiated from all minima. For each disaccharide, the simulation started from the A minimum was conducted for 50 ns, while the other minima were explored for 10 ns. The simulations revealed two stable minima for both compounds. For methyl β-cellobioside, the simulation minima in aqueous solution were shifted from their adiabatic map counterparts, while the simulation minima for methyl β-laminarabioside coincided with the corresponding adiabatic map minima. To validate the simulation results, NMR-derived NOEs and coupling constants across the glycoside linkage, 3JHC and 3JCH, were compared with values calculated from the MD trajectories. For each disaccharide, the best agreement was obtained for the simulations started at the A minimum. For both compounds, inter-ring water bridges in combination with the direct hydrogen bonds between the same groups were found to be determining factors for the overall solution structure of the disaccharides which differed from solid-state structures. Comparison with helical parameters showed that the preferred glycosidic dihedral configurations in the methyl β-cellobioside simulation were not highly compatible with the structure of cellulose, but that curdlan helix structures agreed relatively well with the methyl β-laminarabioside simulation. Polymers generated using glycosidic dihedral angles from the simulations revealed secondary structure motifs that that may help to elucidate polymer associations and small-molecule binding.