An improved OPLS–AA force field for carbohydrates
- 27 September 2002
- journal article
- Published by Wiley in Journal of Computational Chemistry
- Vol. 23 (15) , 1416-1429
- https://doi.org/10.1002/jcc.10139
Abstract
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 alpha-D-glucose and alpha-D-galactose the ratios for the three different conformations of the hydroxymethylene group and the average (3)J(H,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.Keywords
This publication has 39 references indexed in Scilit:
- Reversible peptide folding: Dependence on molecular force field usedJournal of Computational Chemistry, 2000
- Importance of selecting proper basis set in quantum mechanical studies of potential energy surfaces of carbohydratesJournal of Computational Chemistry, 1999
- Upack program package for crystal structure prediction: Force fields and crystal structure generation for small carbohydrate moleculesJournal of Computational Chemistry, 1999
- OPLS all-atom force field for carbohydratesJournal of Computational Chemistry, 1997
- A carbohydrate force field for amber and its application to the study of saccharide to surface adsorptionJournal of Molecular Structure: THEOCHEM, 1997
- Parametrization of GROMOS force field for oligosaccharides and assessment of efficiency of molecular dynamics simulationsJournal of Computational Chemistry, 1996
- An extension of the GROMOS force field for carbohydrates, resulting in improvement of the crystal structure determination of α-D-galactoseActa crystallographica Section B, Structural science, crystal engineering and materials, 1995
- A force field for monosaccharides and (1 → 4) linked polysaccharidesJournal of Computational Chemistry, 1994
- Accuracy of free energies of hydration for organic molecules from 6‐31g*‐derived partial chargesJournal of Computational Chemistry, 1993
- Moleküldynamik‐Computersimulationen; Methodik, Anwendungen und Perspektiven in der ChemieAngewandte Chemie, 1990