Hydration at glycosidic linkages of malto- and cello-oligosaccharides in aqueous solution from molecular dynamics simulation: Effect of conformational flexibility

Molecular simulations of dodecyl-β-maltoside micelles in water: influence of the headgroup conformation and force field parameters

This paper deals with the development and validation of new potential parameter sets, based on the CHARMM36 and GLYCAM06 force fields, to simulate micelles of the two anomeric forms (α and β) of N-dodecyl-β-maltoside (C(12)G(2)), a surfactant widely used in the extraction and purification of membrane proteins. In this context, properties such as size, shape, internal structure, and hydration of the C(12)G(2) anomer micelles were thoroughly investigated by molecular dynamics simulations and the results compared with experiments. Additional simulations were also performed with the older CHARMM22 force field for carbohydrates (Kuttel, M.; et al. J. Comput. Chem. 2002, 23, 1236-1243). We find that our CHARMM and GLYCAM parameter sets yield similar results in the case of properties related to the micelle structure but differ for other properties such as the headgroup conformation or the micelle hydration. In agreement with experiments, our results show that for all model potentials the β-C(12)G(2) micelles have a more pronounced ellipsoidal shape than those containing α anomers. The computed radius of gyration is 20.2 and 25.4 Å for the α- and β-anomer micelles, respectively. Finally, we show that depending on the potential the water translational diffusion of the interfacial water is 7-11.5 times slower than that of bulk water due to the entrapment of the water in the micelle crevices. This retardation is independent of the headgroup in α- or β-anomers.

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

Explicit-solvent molecular dynamics (MD) simulations of the 11 glucose-based disaccharides in water at 300K and 1bar are reported. The simulations were carried out with the GROMOS 45A4 force-field and the sampling along the glycosidic dihedral angles phi and psi was artificially enhanced using the local elevation umbrella sampling (LEUS) method. The trajectories are analyzed in terms of free-energy maps, stable and metastable conformational states (relative free energies and estimated transition timescales), intramolecular H-bonds, single molecule configurational entropies, and agreement with experimental data. All disaccharides considered are found to be characterized either by a single stable (overwhelmingly populated) state ((1-->n)-linked disaccharides with n=1, 2, 3, or 4) or by two stable (comparably populated and differing in the third glycosidic dihedral angle omega ; gg or gt) states with a low interconversion barrier ((1-->6)-linked disaccharides). Metastable (anti-phi or anti-psi) states are also identified with relative free energies in the range of 8-22 kJ mol(-1). The 11 compounds can be classified into four families: (i) the alpha(1-->1)alpha-linked disaccharide trehalose (axial-axial linkage) presents no metastable state, the lowest configurational entropy, and no intramolecular H-bonds; (ii) the four alpha(1-->n)-linked disaccharides (n=1, 2, 3, or 4; axial-equatorial linkage) present one metastable (anti-psi) state, an intermediate configurational entropy, and two alternative intramolecular H-bonds; (iii) the four beta(1-->n)-linked disaccharides (n=1, 2, 3, or 4; equatorial-equatorial linkage) present two metastable (anti-phi and anti-psi) states, an intermediate configurational entropy, and one intramolecular H-bond; (iv) the two (1-->6)-linked disaccharides (additional glycosidic dihedral angle) present no (isomaltose) or a pair of (gentiobiose) metastable (anti-phi) states, the highest configurational entropy, and no intramolecular H-bonds. The observed conformational preferences appear to be dictated by four main driving forces (ring conformational preferences, exo-anomeric effect, steric constraints, and possible presence of a third glycosidic dihedral angle), leaving a secondary role to intramolecular H-bonding and specific solvation effects. In spite of the weak conformational driving force attributed to solvent-exposed H-bonds in water (highly polar protic solvent), intramolecular H-bonds may still have a significant influence on the physico-chemical properties of the disaccharide by decreasing its hydrophilicity. Along with previous work, the results also complete the suggestion of a spectrum of approximate transition timescales for carbohydrates up to the disaccharide level, namely: approximately 30 ps (hydroxyl groups), approximately 1 ns (free lactol group, free hydroxymethyl groups, glycosidic dihedral angleomega in (1-->6)-linked disaccharides), approximately 10 ns to 2 micros (ring conformation, glycosidic dihedral angles phi and psi). The calculated average values of the glycosidic torsional angles agree well with the available experimental data, providing validation for the force-field and simulation methodology employed.

GLYCAM06: a generalizable biomolecular force field. Carbohydrates

A new derivation of the GLYCAM06 force field, which removes its previous specificity for carbohydrates, and its dependency on the AMBER force field and parameters, is presented. All pertinent force field terms have been explicitly specified and so no default or generic parameters are employed. The new GLYCAM is no longer limited to any particular class of biomolecules, but is extendible to all molecular classes in the spirit of a small-molecule force field. The torsion terms in the present work were all derived from quantum mechanical data from a collection of minimal molecular fragments and related small molecules. For carbohydrates, there is now a single parameter set applicable to both alpha- and beta-anomers and to all monosaccharide ring sizes and conformations. We demonstrate that deriving dihedral parameters by fitting to QM data for internal rotational energy curves for representative small molecules generally leads to correct rotamer populations in molecular dynamics simulations, and that this approach removes the need for phase corrections in the dihedral terms. However, we note that there are cases where this approach is inadequate. Reported here are the basic components of the new force field as well as an illustration of its extension to carbohydrates. In addition to reproducing the gas-phase properties of an array of small test molecules, condensed-phase simulations employing GLYCAM06 are shown to reproduce rotamer populations for key small molecules and representative biopolymer building blocks in explicit water, as well as crystalline lattice properties, such as unit cell dimensions, and vibrational frequencies.