Application of restrained minimization, simulated annealing and molecular dynamics simulations for the conformational analysis of oligosaccharides

Abstract

The purpose of the present study was to determine the confidence with which the small number of ¹H NMR nuclear Overhauser effect (NOE) distance constraints measurable across glycosidic linkages in oligosaccharides could be used for solution conformational analysis. This was assessed by use of these constraints in restrained molecular mechanical minimization of the tetrasaccharide Galβ1–4(Fucα1-3)Glc-NAcβ1–3Gal, a model compound of the Lewis-X antigenic determinant. This presents a particularly severe test case in view of extreme resonance overlap and a dearth of inter-residue distance constraints. It is concluded that these constraints, when used in conventional restrained minimization, result in the generation of ‘virtual conformations’ and local minima about glycosidic linkages. However, these restraints are nevertheless found to be useful in the initial stages of a conformational analysis strategy involving restrained minimization combined with dynamical simulated annealing to define more accurately the global minimum energy configuration, together with molecular dynamics simulation to explore conformational mobility about this minimum. Theoretical ROE values calculated over the time course of the MD simulation, using a formalism appropriate for the time scale of the internal motion, are compared with those obtained experimentally in the oligosaccharide.