An improved OPLS–AA force field for carbohydrates

This work describes an improved version of the original OPLS–all atom (OPLS–AA) force field for carbohydrates (Damm et al., J Comp Chem 1997, 18, 1955). The improvement is achieved by applying additional scaling factors for the electrostatic interactions between 1,5‐ and 1,6‐interactions. This new model is tested first for improving the conformational energetics of 1,2‐ethanediol, the smallest polyol. With a 1,5‐scaling factor of 1.25 the force field calculated relative energies are in excellent agreement with the ab initio‐derived data. Applying the new 1,5‐scaling makes it also necessary to use a 1,6‐scaling factor for the interactions between the C4 and C6 atoms in hexopyranoses. After torsional parameter fitting, this improves the conformational energetics in comparison to the OPLS–AA force field. The set of hexopyranoses included in the torsional parameter derivation consists of the two anomers of D‐glucose, D‐mannose, and D‐galactose, as well as of the methyl‐pyranosides of D‐glucose, D‐mannose. Rotational profiles for the rotation of the exocyclic group and of different hydroxyl groups are also compared for the two force fields and at the ab initio level of theory. The new force field reduces the overly high barriers calculated using the OPLS–AA force field. This leads to better sampling, which was shown to produce more realistic conformational behavior for hexopyranoses in liquid simulation. From 10‐ns molecular dynamics (MD) simulations of α‐D‐glucose and α‐D‐galactose the ratios for the three different conformations of the hydroxymethylene group and the average 3JH,H coupling constants are derived and compared to experimental values. The results obtained for OPLS–AA–SEI force field are in good agreement with experiment whereas the properties derived for the OPLS–AA force field suffer from sampling problems. The undertaken investigations show that the newly derived OPLS–AA–SEI force field will allow simulating larger carbohydrates or polysaccharides with improved sampling of the hydroxyl groups. © 2002 Wiley Periodicals, Inc. J Comput Chem 23: 1416–1429, 2002

[1]  Jan Kroon,et al.  An extension of the GROMOS force field for carbohydrates, resulting in improvement of the crystal structure determination of α-D-galactose , 1995 .

[2]  Robert J. Woods,et al.  Molecular Mechanical and Molecular Dynamic Simulations of Glycoproteins and Oligosaccharides. 1. GLYCAM_93 Parameter Development , 1995 .

[3]  F. D. Leeuw,et al.  The relationship between proton-proton NMR coupling constants and substituent electronegativities—I : An empirical generalization of the karplus equation , 1980 .

[4]  W. L. Jorgensen,et al.  Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids , 1996 .

[5]  Donald G. Truhlar,et al.  Exo-anomeric effects on energies and geometries of different conformations of glucose and related systems in the gas phase and aqueous solution , 1997 .

[6]  Dan N. Bernardo,et al.  An Anisotropic Polarizable Water Model: Incorporation of All-Atom Polarizabilities into Molecular Mechanics Force Fields , 1994 .

[7]  L Wang,et al.  The early stage of folding of villin headpiece subdomain observed in a 200-nanosecond fully solvated molecular dynamics simulation. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[8]  P. Kollman,et al.  Encyclopedia of computational chemistry , 1998 .

[9]  Jenn-Huei Lii,et al.  Importance of selecting proper basis set in quantum mechanical studies of potential energy surfaces of carbohydrates , 1999 .

[10]  Michael Schlenkrich,et al.  Force field parameters for carbohydrates , 1996 .

[11]  X. Daura,et al.  Reversible peptide folding in solution by molecular dynamics simulation. , 1998, Journal of molecular biology.

[12]  James W. Brown,et al.  AB INITIO STUDIES OF THE EXOCYCLIC HYDROXYMETHYL ROTATIONAL SURFACE IN ALPHA -D-GLUCOPYRANOSE , 1996 .

[13]  James W. Brown,et al.  Exocyclic Hydroxymethyl Rotational Conformers of β- and α-d-Glucopyranose in the Gas Phase and Aqueous Solution , 1998 .

[14]  W. F. V. Gunsteren,et al.  Moleküldynamik‐Computersimulationen; Methodik, Anwendungen und Perspektiven in der Chemie , 1990 .

[15]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[16]  P. Mark Rodger,et al.  Lifetime of a Hydrogen Bond in Aqueous Solutions of Carbohydrates , 1999 .

[17]  Hiroshi Ohrui,et al.  1H-NMR studies of (6r)- and (6s)-deuterated d-hexoses: assignment of the preferred rotamers about C5C6 bond of D-glucose and D-galactose derivatives in solutions , 1984 .

[18]  Charles W. Bauschlicher,et al.  A comparison of the accuracy of different functionals , 1995 .

[19]  P. Schleyer.  Encyclopedia of computational chemistry , 1998 .

[20]  J. Brady,et al.  A revised potential-energy surface for molecular mechanics studies of carbohydrates. , 1988, Carbohydrate research.

[21]  Alan E. Mark,et al.  The GROMOS96 Manual and User Guide , 1996 .

[22]  Imre G. Csizmadia,et al.  Density functional conformational analysis of 1,2-ethanediol , 1995 .

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

[24]  Hiroshi Ohrui,et al.  The preferred rotamer about the C5—C6 bond of D-galactopyranoses and the stereochemistry of dehydrogenation by D-galactose oxidase , 1987 .

[25]  Hanoch Senderowitz,et al.  Carbohydrates: United Atom AMBER* Parameterization of Pyranoses and Simulations Yielding Anomeric Free Energies , 1996 .

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

[27]  Hiroshi Ohrui,et al.  1H NMR Analyses of Rotameric Distribution of C5-C6 bonds of D-Glucopyranoses in Solution , 1988 .

[28]  P. Kollman,et al.  Application of RESP charges to calculate conformational energies, hydrogen bond energies, and free energies of solvation , 1993 .

[29]  X. Daura,et al.  Peptide Folding: When Simulation Meets Experiment , 1999 .

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

[31]  Donald G. Truhlar,et al.  Quantum Chemical Conformational Analysis of 1,2-Ethanediol: Correlation and Solvation Effects on the Tendency To Form Internal Hydrogen Bonds in the Gas Phase and in Aqueous Solution , 1994 .

[32]  Karl-Heinz Ott,et al.  Parametrization of GROMOS force field for oligosaccharides and assessment of efficiency of molecular dynamics simulations , 1996, J. Comput. Chem..

[33]  Aatto Laaksonen,et al.  Molecular Dynamics Simulation of the α-d-Manp-(1 → 3)- β-d-Glcp-OMe Disaccharide in Water and Water/DMSO Solution , 1999 .

[34]  Heather A. Carlson,et al.  Accuracy of free energies of hydration for organic molecules from 6‐31g*‐derived partial charges , 1993, J. Comput. Chem..

[35]  Kenneth M. Merz,et al.  A carbohydrate force field for amber and its application to the study of saccharide to surface adsorption , 1997 .

[36]  Kenneth M. Merz,et al.  A force field for monosaccharides and (1 → 4) linked polysaccharides , 1994, J. Comput. Chem..

[37]  P. Kollman,et al.  Pathways to a protein folding intermediate observed in a 1-microsecond simulation in aqueous solution. , 1998, Science.

[38]  Wolfgang Damm,et al.  Reversible peptide folding: Dependence on molecular force field used , 2000 .

[39]  William L. Jorgensen,et al.  A comprehensive study of the rotational energy profiles of organic systems by ab initio MO theory, forming a basis for peptide torsional parameters , 1995, J. Comput. Chem..

[40]  Donald G. Truhlar,et al.  Factors controlling relative stability of anomers and hydroxymethyl conformers of glucopyranose , 1998, J. Comput. Chem..

[41]  Jan Kroon,et al.  Upack program package for crystal structure prediction: Force fields and crystal structure generation for small carbohydrate molecules , 1999, J. Comput. Chem..

[42]  Peter A. Kollman,et al.  Use of Locally Enhanced Sampling in Free Energy Calculations: Testing and Application to the α → β Anomerization of Glucose , 1998 .