Kinetics of conformer formation of glucose and maltose in aqueous solutions

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.

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.

Ultrasonic spectrometry of aqueous solutions of alkyl maltosides: kinetics of micelle formation and head-group isomerization

At frequencies between 100 kHz and 400 MHz, ultrasonic attenuation spectra are measured at 25 degrees C for aqueous solutions of hexyl-, heptyl-, octyl-, nonyl-, and decyl-beta-D-maltopyranoside as well as of decyl-alpha-D-maltopyranoside. The spectra with surfactant concentration c above the relevant critical micelle concentration (cmc) display three relaxation terms with discrete relaxation times. That with a relaxation time between 0.1 and 1.2 micros is due to exchange of monomers between micelles and the suspending phase. It is discussed in the light of the Teubner-Kahlweit-Aniansson-Wall model of the formation/decay kinetics of systems with Gaussian size distribution of micelles. The relaxation parameters are compared to those for solutions of other non-ionic surfactants, such as alkyl monoglycosides and poly(ethylene glycol) monoalkyl ethers. At c < cmc this low-frequency relaxation term is missing and at c approximately = cmc it is broadened, as is characteristic of solutions of oligomeric molecular structures rather than proper micelles. The relaxation terms with relaxation times between 6 and 15 ns and 0.7 and 2.3 ns reveal head-group rotations around glycosidic angles and isomerization of the exocyclic hydroxymethyl group, respectively. These unimolecular reactions are also examined with a view to solutions of alkyl monoglycosides as well as of glucose and maltose.

Octylglucopyranoside and cyclodextrin in water. Self-aggregation and complex formation

At 25 degrees C ultrasonic attenuation spectra between 100 kHz and 400 MHz as well as sound velocities and densities of aqueous solutions of the surfactant n-octyl-beta-d-glucopyranoside and of the cage compound alpha-cyclodextrin have been measured. The liquid system reveals a critical association concentration (cac) exceeding the cmc of the surfactant almost by the cyclodextrin concentration and thus indicating a significant formation of cyclodextrin-surfactant inclusion complexes. The ultrasonic spectra display altogether four relaxation regions. The one with largest relaxation time (0.27 mus < or = tau(1) < or = 1.6 mus) exhibits a noticeable amplitude at surfactant concentrations larger than the cac only. It is assigned to the monomer exchange between the micelles and the suspending phase. A term with relaxation time tau(2) (33 ns < or = tau(2) < or = 135 ns) is characteristic for solutions containing both solutes. It is assumed to reflect the inclusion complex formation. Complexes with 1:1, 2:1, and 1:2 stoichiometry appear to exist. The terms at higher frequencies (4.8 ns < or = tau(3) < or = 9.8 ns; 0.8 ns <or = tau(4) < or = 2 ns) are due to fluctuations of the carbohydrate residues around the glucosidic bond angles and to the rotational isomerization of the exocyclic hydroxymethyl groups, respectively.