Solution structure of hyaluronic acid oligomers by experimental and theoretical NMR, and molecular dynamics simulation.

The conformational properties of hyaluronic acid (HA) oligomers in aqueous solution were investigated by combining high-resolution NMR experimental results, theoretical simulation of NMR two-dimensional (2D) spectra by Complete Relaxation Matrix Analysis (CORMA), and molecular dynamics calculations. New experimental findings recorded for the tetra- and hexasaccharides enabled the stiffness of the HA and its viscoelastic properties to be interpreted. In particular, rotating frame nuclear Overhauser effect spectroscopy spectra provided new information about the arrangement of the glycosidic linkage. From (13)C NMR relaxation the rotational correlation time (tau(c)) were determined. The tau(c) were employed in the calculation of geometrical constraints, by using the MARDIGRAS algorithm. Restrained simulated annealing and 1 ns of unrestrained molecular dynamic simulations were performed on the hexasaccharide in a box of 1215 water molecules. The beta(1 --> 3) and beta(1 --> 4) glycosidic links were found to be rigid. The lack of rotational degree of freedom is due to direct and/or water-mediated interresidue hydrogen bonding. Both single or tandem water bridges were found between carboxylate group and N-acetil group. The carboxylate group of glucuronic acid is not involved in a direct link with the amide group of N-acetyl glucosamine and this facilitated bonding between the residue and the water molecules.