Conformational Properties of Glycosaminoglycan Disaccharides: A Molecular Dynamics Study
- 22 September 2021
- journal article
- research article
- Published by American Chemical Society (ACS) in The Journal of Physical Chemistry B
- Vol. 125 (39) , 10900-10916
Abstract
The structure and conformation of glycosaminoglycans (GAGs) are of central importance to understand the mechanisms behind their functions in biological systems. Due to the inherent chemical and structural heterogeneity of GAGs, focusing on longer, naturally existing GAG chains hinders drawing conclusions on the influence of the chemical functionalization on the basic conformational degree of freedom, that is, the dynamic shape of glycosidic linkage present in the particular disaccharide repeating unit. In the present study, we have considered the complete set of 106 GAG-related disaccharides, being potential building blocks for longer GAG chains (including hyaluronan, chondroitin, keratan, dermatan, and heparan). Both the unfunctionalized units and all possible combinations of either partially or fully sulfated derivatives contribute to this number. The unbiased and enhanced sampling molecular dynamics simulations provide a link to understand the influence of sulfation on the conformational properties of GAG glycosidic linkages. Residue–residue hydrogen bonding is not significant for either the glycosidic linkage conformation or its flexibility. It was found that in the majority of cases, the dominating conformation of the linkage is weakly affected by sulfation and the main role is played by the steric and stereoelectronic effects. However, there exist numerous cases where sulfation increases the contribution of alternative conformations to a nonnegligible extent and, in some rare cases (restricted to disaccharides building heparan), leads to the reorientation of the glycosidic linkage. The identified sulfation sites, being the most important in this context, are C6 and C3 at the GlcNAc residue. Finally, the full set of free energy maps relying on the glycosidic dihedral angle values for diverse GAG disaccharides are provided; they may be used for further studies, focused on longer GAG chains.Keywords
Funding Information
- Narodowe Centrum Nauki (2019/35/B/ST4/01149 OPUS 18)
This publication has 42 references indexed in Scilit:
- Intramolecular hydrogen-bonding in aqueous carbohydrates as a cause or consequence of conformational preferences: a molecular dynamics study of cellobiose stereoisomersEuropean Biophysics Journal, 2013
- CHARMM Additive All-Atom Force Field for Carbohydrate Derivatives and Its Utility in Polysaccharide and Carbohydrate–Protein ModelingJournal of Chemical Theory and Computation, 2011
- Free Energy Landscapes of Iduronic Acid and Related MonosaccharidesJournal of the American Chemical Society, 2010
- CHARMM Additive All-Atom Force Field for Glycosidic Linkages between HexopyranosesJournal of Chemical Theory and Computation, 2009
- Canonical sampling through velocity rescalingThe Journal of Chemical Physics, 2007
- Free-Energy Landscape for β Hairpin Folding from Combined Parallel Tempering and MetadynamicsJournal of the American Chemical Society, 2006
- Hyaluronan: The Local Solution Conformation Determined by NMR and Computer Modeling is Close to a Contracted Left-handed 4-Fold HelixJournal of Molecular Biology, 2006
- Replica-exchange molecular dynamics method for protein foldingChemical Physics Letters, 1999
- D-Idose: a one- and two-dimensional NMR investigation of solution composition and conformationThe Journal of Organic Chemistry, 1986
- Comparison of simple potential functions for simulating liquid waterThe Journal of Chemical Physics, 1983