Is N-acetyl-d-glucosamine a rigid 4C1 chair?

Understanding microsecond-timescale dynamics is crucial to establish three-dimensional (3D) structure–activity relationships in sugars but has been intractable to experiments and simulations. As a consequence, whether arguably the most important chemical scaffold in glycobiology, N-acetyl-d-glucosamine (GlcNAc), deviates from a rigid 4C1 chair is unknown. Here, conformer populations and exchange kinetics were quantified from the longest aqueous carbohydrate simulations to date (0.2 ms total) of GlcNAc, four derivatives from heparan sulfate and their methylglycosides. Unmodified GlcNAc took 3–5 μs to reach a conformational equilibrium, which comprised a metastable 4C1 chair that underwent 4C1 ↔ 1C4 transitions at a predicted forward rate of 0.8 μs−1 with an average 1C4-chair lifetime of 3 ns. These predictions agree with high-resolution crystallography and nuclear magnetic resonance but not with the hypothesis that GlcNAc is a rigid 4C1 chair, concluded from previous experimental analyses and non-aqueous modeling. The methylglycoside was calculated to have a slower forward rate (0.3 μs−1) and a more stable 4C1 conformer (0.2 kcal mol−1), suggesting that pivotal 3D intermediates (particularly 2SO, 1S5 and B2,5) increased in energy, and water was implicated as a major cause. Sulfonation (N-, 3-O and 6-O) significantly augmented this effect by blocking pseudorotation, but did not alter the rotational preferences of hydroyxl or hydroxymethyl groups. We therefore propose that GlcNAc undergoes puckering exchange that is dependent on polymerization and sulfo substituents. Our analyses, and 3D model of the equilibrium GlcNAc conformer in water, can be used as dictionary data and present new opportunities to rationally modify puckering and carbohydrate bioactivity, with diverse applications from improving crop yields to disease amelioration.

[1]  Jacopo Tomasi,et al.  A new integral equation formalism for the polarizable continuum model: Theoretical background and applications to isotropic and anisotropic dielectrics , 1997 .

[2]  M. Yamabhai,et al.  Chitin Research Revisited , 2010, Marine drugs.

[3]  Matthias Rief,et al.  Single-molecule force spectroscopy on polysaccharides by AFM – nanomechanical fingerprint of α-(1,4)-linked polysaccharides , 1999 .

[4]  M J Harvey,et al.  ACEMD: Accelerating Biomolecular Dynamics in the Microsecond Time Scale. , 2009, Journal of chemical theory and computation.

[5]  H. Nothaft,et al.  Protein glycosylation in bacteria: sweeter than ever , 2010, Nature Reviews Microbiology.

[6]  T. Maruyama,et al.  Conformational changes and anticoagulant activity of chondroitin sulfate following its O-sulfonation. , 1998, Carbohydrate research.

[7]  W. B. Gleason,et al.  Carbohydrates with relevance to the structure of glycosaminoglycans: The crystal structures of 2-deoxy-2-(sulfoamino)-α-d-glucopyranose sodium salt dihydrate, 2-amino-2-deoxy-α,β-d-glucopyranose 3-(hydrogen sulfate) monohydrate, and 2-amino-2-deoxy-α-d-glucopyranose 6-(hydrogen sulfate) monohydrate , 1995 .

[8]  Frank A. Momany,et al.  B3LYP/6-311++G** study of α- and β-d-glucopyranose and 1,5-anhydro-d-glucitol: 4C1 and 1C4 chairs, 3,OB and B3,O boats, and skew-boat conformations , 2004 .

[9]  Jeffrey D. Esko,et al.  Heparan sulphate proteoglycans fine-tune mammalian physiology , 2007, Nature.

[10]  P. Kollman,et al.  An approach to computing electrostatic charges for molecules , 1984 .

[11]  T. Ghosh,et al.  Focus on antivirally active sulfated polysaccharides: from structure-activity analysis to clinical evaluation. , 2009, Glycobiology.

[12]  M. Hricovíni,et al.  Relationship between structure and three-bond proton-proton coupling constants in glycosaminoglycans. , 2007, Carbohydrate research.

[13]  Haiying Liu,et al.  Lessons Learned from the Contamination of Heparin , 2009 .

[14]  Thomas J Boltje,et al.  Opportunities and challenges in synthetic oligosaccharide and glycoconjugate research. , 2009, Nature chemistry.

[15]  F. Mo,et al.  A refined model for N-acetyl-α-d-glucosamine , 1975 .

[16]  L. Poppe Modeling carbohydrate conformations from NMR data : maximum entropy rotameric distribution about the C5-C6 bond in gentiobiose , 1993 .

[17]  A. Almond,et al.  Hyaluronan: the absence of amide-carboxylate hydrogen bonds and the chain conformation in aqueous solution are incompatible with stable secondary and tertiary structure models. , 2006, The Biochemical journal.

[18]  Crystal structure of the CD2-binding domain of CD58 (lymphocyte function-associated antigen 3) at 1.8-A resolution. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[19]  Giovanni Scalmani,et al.  New developments in the polarizable continuum model for quantum mechanical and classical calculations on molecules in solution , 2002 .

[20]  Benedict M. Sattelle,et al.  Free energy landscapes of iduronic acid and related monosaccharides. , 2010, Journal of the American Chemical Society.

[21]  Karl Nicholas Kirschner,et al.  GLYCAM06: A generalizable biomolecular force field. Carbohydrates , 2008, J. Comput. Chem..

[22]  T. Toida,et al.  Effect of 6-O-sulfonate hexosamine residue on anticoagulant activity of fully O-sulfonated glycosaminoglycans , 2000, Glycoconjugate Journal.

[23]  Cornelis Altona,et al.  Prediction of anti and gauche vicinal proton‐proton coupling constants in carbohydrates: A simple additivity rule for pyranose rings , 1980 .

[24]  Tamar Schlick,et al.  Molecular dynamics-based approaches for enhanced sampling of long-time, large-scale conformational changes in biomolecules , 2009, F1000 biology reports.

[25]  B. Mulloy,et al.  Assignment of the 1H-n.m.r. spectra of heparin and heparan sulphate. , 1987, Carbohydrate research.

[26]  Lijiang Yang,et al.  Thermodynamics and kinetics simulations of multi-time-scale processes for complex systems , 2008 .

[27]  A. Almond,et al.  Use of 15N-NMR to resolve molecular details in isotopically-enriched carbohydrates: sequence-specific observations in hyaluronan oligomers up to decasaccharides. , 2004, Glycobiology.

[28]  M. A. Meledeo,et al.  Hexosamine analogs: from metabolic glycoengineering to drug discovery. , 2009, Current opinion in chemical biology.

[29]  Zachary Shriver,et al.  Roles of heparan-sulphate glycosaminoglycans in cancer , 2002, Nature Reviews Cancer.

[30]  H. Wessel,et al.  Conformational flexibility in highly sulfated β-d-glucopyranoside derivatives , 1995 .

[31]  E. Yates,et al.  The crystal and molecular structure of 2-sulfamino-2-deoxy-alpha-D-glucopyranose sodium salt.2H2O (glucosamine 2-sulfate). , 1995, International journal of biological macromolecules.

[32]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[33]  Karl N. Kirschner,et al.  Solvent interactions determine carbohydrate conformation , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[34]  Pauline M Rudd,et al.  A systematic approach to protein glycosylation analysis: a path through the maze. , 2010, Nature chemical biology.

[35]  Christopher B. Barnett,et al.  Ring puckering: a metric for evaluating the accuracy of AM1, PM3, PM3CARB-1, and SCC-DFTB carbohydrate QM/MM simulations. , 2010, The journal of physical chemistry. B.

[36]  D. Cremer,et al.  General definition of ring puckering coordinates , 1975 .

[37]  A. Minton,et al.  Quantitative Characterization of Heparin Binding to Tau Protein , 2009, The Journal of Biological Chemistry.

[38]  R. Linhardt,et al.  Synthetic oligosaccharides as heparin-mimetics displaying anticoagulant properties. , 2003, Current pharmaceutical design.

[39]  Vojtech Spiwok,et al.  Modelling of beta-D-glucopyranose ring distortion in different force fields: a metadynamics study. , 2010, Carbohydrate research.

[40]  J. Herbert,et al.  Synthesis of Conformationally Locked L-Iduronic Acid Derivatives. Direct Evidence for a Critical Role of the Skew-Boat 2S0 Conformer in the Activation of Antithrombin by Heparin. , 2002 .

[41]  M. Forster,et al.  Conformation and dynamics of heparin and heparan sulfate. , 2000, Glycobiology.

[42]  J. Herbert,et al.  Synthesis of conformationally locked L-iduronic acid derivatives: direct evidence for a critical role of the skew-boat 2S0 conformer in the activation of antithrombin by heparin. , 2001, Chemistry.

[43]  A. Imberty,et al.  NMR and molecular modeling reveal key structural features of synthetic nodulation factors. , 2011, Glycobiology.

[44]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[45]  Wayne Boucher,et al.  The CCPN data model for NMR spectroscopy: Development of a software pipeline , 2005, Proteins.

[46]  A. Wlodawer,et al.  Crystal Structure of Interleukin-19 Defines a New Subfamily of Helical Cytokines* , 2003, The Journal of Biological Chemistry.

[47]  G. Colombo,et al.  Conformation of heparin pentasaccharide bound to antithrombin III , 2001 .

[48]  W. V. van Gunsteren,et al.  A fast SHAKE algorithm to solve distance constraint equations for small molecules in molecular dynamics simulations , 2001 .

[49]  Christopher B. Barnett,et al.  Pyranose ring transition state is derived from cellobiohydrolase I induced conformational stability and glycosidic bond polarization. , 2010, Journal of the American Chemical Society.

[50]  M. Nilsson,et al.  The structural plasticity of heparan sulfate NA-domains and hence their role in mediating multivalent interactions is confirmed by high-accuracy 15N-NMR relaxation studies , 2007, Glycoconjugate Journal.

[51]  Alexander Rich,et al.  The double helix: a tale of two puckers , 2003, Nature Structural Biology.

[52]  David P. Anderson,et al.  High-Throughput All-Atom Molecular Dynamics Simulations Using Distributed Computing , 2010, J. Chem. Inf. Model..

[53]  R. Dror,et al.  Long-timescale molecular dynamics simulations of protein structure and function. , 2009, Current opinion in structural biology.

[54]  J. Turnbull,et al.  Interactions of heparin/heparan sulfate with proteins: appraisal of structural factors and experimental approaches. , 2004, Glycobiology.

[55]  T. Maruyama,et al.  Preparation and anticoagulant activity of fully O-sulphonated glycosaminoglycans. , 1999, International journal of biological macromolecules.

[56]  D. Uhrín,et al.  Residual dipolar coupling investigation of a heparin tetrasaccharide confirms the limited effect of flexibility of the iduronic acid on the molecular shape of heparin , 2009, Glycobiology.

[57]  Thomas Dierks,et al.  The heparanome--the enigma of encoding and decoding heparan sulfate sulfation. , 2007, Journal of biotechnology.

[58]  Holger Gohlke,et al.  The Amber biomolecular simulation programs , 2005, J. Comput. Chem..

[59]  J. Turnbull Heparan sulfate glycomics: towards systems biology strategies. , 2010, Biochemical Society transactions.

[60]  Roy Jefferis,et al.  Glycosylation as a strategy to improve antibody-based therapeutics , 2009, Nature Reviews Drug Discovery.

[61]  Nagarajan Vaidehi,et al.  Sulfation patterns of glycosaminoglycans encode molecular recognition and activity , 2006, Nature chemical biology.

[62]  L. Johnson The crystal structure of N-acetyl-a-D-glycosamine , 1966 .

[63]  Yuan-Ping Pang,et al.  Polysaccharide elasticity governed by chair–boat transitions of the glucopyranose ring , 1998, Nature.

[64]  X. Biarnés,et al.  The conformational free-energy landscape of β-D-mannopyranose: evidence for a (1)S(5) → B(2,5) → (O)S(2) catalytic itinerary in β-mannosidases. , 2010, Journal of the American Chemical Society.

[65]  Benedict M. Sattelle,et al.  Less is more when simulating unsulfated glycosaminoglycan 3D‐structure: Comparison of GLYCAM06/TIP3P, PM3‐CARB1/TIP3P, and SCC‐DFTB‐D/TIP3P predictions with experiment , 2010, J. Comput. Chem..

[66]  S. Kennel,et al.  The Heparan Sulfate Motif (GlcNS6S-IdoA2S)3, Common in Heparin, Has a Strict Topography and Is Involved in Cell Behavior and Disease* , 2010, The Journal of Biological Chemistry.

[67]  H. Grubmüller,et al.  ’Chair-boat‘ transitions and side groups affect the stiffness of polysaccharides. , 1999 .

[68]  I. Roberts,et al.  Investigating the Molecular Basis for the Virulence of Escherichia coli K5 by Nuclear Magnetic Resonance Analysis of the Capsule Polysaccharide , 2009, Journal of Molecular Microbiology and Biotechnology.

[69]  S. Grzesiek,et al.  NMRPipe: A multidimensional spectral processing system based on UNIX pipes , 1995, Journal of biomolecular NMR.

[70]  M. Reed,et al.  NMR spectra of oligosaccharides at ultra-high field (900 MHz) have better resolution than expected due to favourable molecular tumbling , 2006, Carbohydrate research.

[71]  W. Goddard,et al.  UFF, a full periodic table force field for molecular mechanics and molecular dynamics simulations , 1992 .

[72]  L. Pedersen,et al.  Anticoagulant heparan sulfate: structural specificity and biosynthesis , 2007, Applied Microbiology and Biotechnology.

[73]  Modesto Orozco,et al.  Recent advances in the study of nucleic acid flexibility by molecular dynamics. , 2008, Current opinion in structural biology.

[74]  Peter J. Reilly,et al.  Modeling of aldopyranosyl ring puckering with MM3 (92) , 1994 .

[75]  L. Hsieh‐Wilson,et al.  The Chemical Neurobiology of Carbohydrates , 2008 .

[76]  M. Wada,et al.  X-ray Structure of Anhydrous β-Chitin at 1 Å Resolution , 2011 .

[77]  D. Rabenstein Heparin and heparan sulfate: structure and function. , 2002, Natural product reports.

[78]  D. Svergun,et al.  Structural basis for sugar recognition, including the Tn carcinoma antigen, by the lectin SNA‐II from Sambucus nigra , 2009, Proteins.

[79]  I. Wilson,et al.  Structure, Receptor Binding, and Antigenicity of Influenza Virus Hemagglutinins from the 1957 H2N2 Pandemic , 2009, Journal of Virology.

[80]  Donald G. Truhlar,et al.  Relative stability of alternative chair forms and hydroxymethyl conformations of β-d-glucopyranose , 1995 .