Modelling of beta-D-glucopyranose ring distortion in different force fields: a metadynamics study.

Modelling of carbohydrate conformations is a challenging task for force field developers. Three carbohydrate force fields, namely GLYCAM06, GROMOS 45a4 and OPLS were evaluated. Free energies of different ring conformations of beta-D-glucopyranose were calculated using metadynamics in vacuum as well as in explicitly modelled water. All three force fields model the (4)C(1) conformation as the most stable by at least 6kJ/mol, as compared to other conformations. Interconversion from the (4)C(1) to any other conformation is associated with a barrier of no lower than 26kJ/mol. The free energy surface calculated in the GLYCAM06 force field is in remarkably good agreement with the recent Car-Parrinello metadynamics study. The effect of a water environment is relatively low and analogous in all tested force fields. Namely, the presence of water stabilizes the upper-left ((3,O)B) versus bottom-right (B(3,O)) area of Stoddard's plot, relative to the situation in vacuum. Comparison of free and potential surfaces is also provided for vacuum calculations.

[1]  A. Laio,et al.  Metadynamics: a method to simulate rare events and reconstruct the free energy in biophysics, chemistry and material science , 2008 .

[2]  H. Gilbert,et al.  Substrate Distortion by a Lichenase Highlights the Different Conformational Itineraries Harnessed by Related Glycoside Hydrolases , 2006 .

[3]  H. Verli,et al.  Depiction of the forces participating in the 2-O-sulfo-alpha-L-iduronic acid conformational preference in heparin sequences in aqueous solutions. , 2008, Carbohydrate research.

[4]  G. P. Johnson,et al.  Comparison of different force fields for the study of disaccharides. , 2009, Carbohydrate research.

[5]  Teppei Suzuki.  The hydration of glucose: the local configurations in sugar-water hydrogen bonds. , 2008, Physical chemistry chemical physics : PCCP.

[6]  Eric J. Sorin,et al.  Exploring the helix-coil transition via all-atom equilibrium ensemble simulations. , 2005, Biophysical journal.

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

[8]  J. Kastrup,et al.  Structural Rearrangements of Sucrose Phosphorylase from Bifidobacterium adolescentis during Sucrose Conversion* , 2006, Journal of Biological Chemistry.

[9]  Gerrit Groenhof,et al.  GROMACS: Fast, flexible, and free , 2005, J. Comput. Chem..

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

[11]  William L. Jorgensen,et al.  OPLS all‐atom force field for carbohydrates , 1997 .

[12]  Piotr E Marszalek,et al.  Identification of sugar isomers by single-molecule force spectroscopy. , 2006, Journal of the American Chemical Society.

[13]  G. A. Jeffrey,et al.  Stereographic representation of the cremer-pople ring-puckering parameters for pyranoid rings☆ , 1979 .

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

[15]  P. Grootenhuis,et al.  Parametrization and application of CHEAT95, and extended atom force field for hydrated oligosaccharides , 1995 .

[16]  T. A. Jones,et al.  Crystallographic evidence for substrate ring distortion and protein conformational changes during catalysis in cellobiohydrolase Ce16A from trichoderma reesei. , 1999, Structure.

[17]  Vojtech Spiwok,et al.  Conformational free energy surface of alpha-N-acetylneuraminic acid: an interplay between hydrogen bonding and solvation. , 2009, The journal of physical chemistry. B.

[18]  Martin Frank,et al.  Carbohydrate Structure Suite (CSS): analysis of carbohydrate 3D structures derived from the PDB , 2004, Nucleic Acids Res..

[19]  Vojtech Spiwok,et al.  Metadynamics in essential coordinates: free energy simulation of conformational changes. , 2007, The journal of physical chemistry. B.

[20]  Philippe H. Hünenberger,et al.  Using the local elevation method to construct optimized umbrella sampling potentials: Calculation of the relative free energies and interconversion barriers of glucopyranose ring conformers in water , 2010, J. Comput. Chem..

[21]  Alessandra Magistrato,et al.  Dissociation of minor groove binders from DNA: insights from metadynamics simulations , 2008, Nucleic acids research.

[22]  Paolo Carloni,et al.  Mechanism of Action of Cyclophilin A Explored by Metadynamics Simulations , 2009, PLoS Comput. Biol..

[23]  Kjeld Rasmussen,et al.  A comparison and chemometric analysis of several molecular mechanics force fields and parameter sets applied to carbohydrates , 1998 .

[24]  P. Hünenberger,et al.  Conformation, dynamics, solvation and relative stabilities of selected beta-hexopyranoses in water: a molecular dynamics study with the GROMOS 45A4 force field. , 2007, Carbohydrate research.

[25]  Jim Pfaendtner,et al.  Nucleotide-dependent conformational states of actin , 2009, Proceedings of the National Academy of Sciences.

[26]  Alessandro Laio,et al.  Efficient exploration of reactive potential energy surfaces using Car-Parrinello molecular dynamics. , 2003, Physical review letters.

[27]  D. Cremer,et al.  Molecular orbital theory of the electronic structure of organic compounds. XXIII. Pseudorotation in saturated five-membered ring compounds , 1975 .

[28]  H. Berendsen,et al.  Interaction Models for Water in Relation to Protein Hydration , 1981 .

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

[30]  Francesco L Gervasio,et al.  Metadynamics simulation of prion protein: beta-structure stability and the early stages of misfolding. , 2006, Journal of the American Chemical Society.

[31]  M. Lascombe,et al.  Atomic (0.94 A) resolution structure of an inverting glycosidase in complex with substrate. , 2002, Journal of molecular biology.

[32]  A. Laio,et al.  Escaping free-energy minima , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[33]  Alessandro Laio,et al.  The conformational free energy landscape of beta-D-glucopyranose. Implications for substrate preactivation in beta-glucoside hydrolases. , 2007, Journal of the American Chemical Society.

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

[35]  W. Hol,et al.  Crystal structure of heat‐labile enterotoxin from Escherichia coli with increased thermostability introduced by an engineered disulfide bond in the A subunit , 1997, Protein science : a publication of the Protein Society.

[36]  D. Klug,et al.  Structural transformations in carbon under extreme pressure: beyond diamond. , 2009, The Journal of chemical physics.

[37]  Régis Pomès,et al.  Enhancing the accuracy, the efficiency and the scope of free energy simulations. , 2005, Current opinion in structural biology.

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

[39]  F. Pederiva,et al.  On the calculation of puckering free energy surfaces. , 2009, The Journal of chemical physics.

[40]  R. Broglia,et al.  Exploring the protein G helix free‐energy surface by solute tempering metadynamics , 2007, Proteins.

[41]  V. Spiwok,et al.  Metadynamics modelling of the solvent effect on primary hydroxyl rotamer equilibria in hexopyranosides. , 2009, Carbohydrate research.

[42]  Lars Olsen,et al.  Evaluation of carbohydrate molecular mechanical force fields by quantum mechanical calculations. , 2004, Carbohydrate research.

[43]  W. L. Jorgensen.  Free energy calculations: a breakthrough for modeling organic chemistry in solution , 1989 .

[44]  K. Sakka,et al.  Crystal Structure of Cel44A, a Glycoside Hydrolase Family 44 Endoglucanase from Clostridium thermocellum* , 2007, Journal of Biological Chemistry.

[45]  V. Spiwok,et al.  Potential-energy and free-energy surfaces of glycyl-phenylalanyl-alanine (GFA) tripeptide: experiment and theory. , 2008, Chemistry.

[46]  G. Volckaert,et al.  Structural analysis of a glycoside hydrolase family 43 arabinoxylan arabinofuranohydrolase in complex with xylotetraose reveals a different binding mechanism compared with other members of the same family. , 2009, The Biochemical journal.

[47]  K. Ueda,et al.  Several transition states from (4)C(1) to skew conformations of beta-D-glucopyranose. , 2009, Carbohydrate Research.

[48]  Roberto D. Lins,et al.  A new GROMOS force field for hexopyranose‐based carbohydrates , 2005, J. Comput. Chem..

[49]  Robert J Woods,et al.  Structural glycobiology: A game of snakes and ladders , 2008, Glycobiology.

[50]  M. Hricovíni.  B3LYP/6-311++G** study of structure and spin-spin coupling constant in methyl 2-O-sulfo-alpha-L-iduronate. , 2006, Carbohydrate research.

[51]  B. Henrissat,et al.  Structural determinants of substrate specificity in family 1 beta-glucosidases: novel insights from the crystal structure of sorghum dhurrinase-1, a plant beta-glucosidase with strict specificity, in complex with its natural substrate. , 2004, The Journal of biological chemistry.

[52]  G. Davies,et al.  Mapping the conformational itinerary of β-glycosidases by X-ray crystallography , 2003 .

[53]  D. Quigley,et al.  Metadynamics simulations of calcite crystallization on self-assembled monolayers. , 2009, The Journal of chemical physics.

[54]  Clemens Lakner,et al.  An Enzymatic Atavist Revealed in Dual Pathways for Water Activation , 2008, PLoS biology.

[55]  G. Davies,et al.  Structure of the Humicola insolens cellobiohydrolase Cel6A D416A mutant in complex with a non-hydrolysable substrate analogue, methyl cellobiosyl-4-thio-beta-cellobioside, at 1.9 A. , 2002, Acta crystallographica. Section D, Biological crystallography.

[56]  Peter A. Kollman,et al.  FREE ENERGY CALCULATIONS : APPLICATIONS TO CHEMICAL AND BIOCHEMICAL PHENOMENA , 1993 .