Home > Research > Publications & Outputs > Evaluation and optimization of a force field fo...


Text available via DOI:

View graph of relations

Evaluation and optimization of a force field for crystalline forms of mannitol and sorbitol

Research output: Contribution to Journal/MagazineJournal articlepeer-review

  • H. de Waard
  • A. Amani
  • J. Kendrick
  • W. L. J. Hinrichs
  • H. W. Frijlink
  • J. Anwar
<mark>Journal publication date</mark>14/01/2010
<mark>Journal</mark>Journal of Physical Chemistry B
Issue number1
Number of pages8
Pages (from-to)429-436
Publication StatusPublished
<mark>Original language</mark>English


Two force fields, the GROMOS53A5/53A6 (United atom) and the AMBER95 (all atom) parameter sets, Coupled with partial atomic charges derived from quantum mechanical calculations were evaluated for their ability to reproduce the known crystalline forms of the polyols mannitol and sorbitol. The force fields were evaluated using molecular dynamics Simulations at 10 K (which is akin to potential energy minimization) with the simulation cell lengths and angles free to evolve. Both force fields performed relatively poorly, not being able to Simultaneously reproduce all of the crystal structures within a 5% deviation level. The parameter sets were then systematically optimized using sensitivity analysis, and a revised AMBER95 set was found to reproduce the crystal Structures with less than 5% deviation from experiment. The stability of the various crystalline forms for each of the parameter sets (original and revised) was then assessed in extended MD simulations at 298 K and I bar covering 1 ns simulation time. The AMBER95 parameter sets (original and revised) were found to be effective ill reproducing the crystal Structures in these more stringent tests. Remarkably, the performance of the original AMBER95 parameter set was found to be slightly better than that of the revised set in these simulations at 298 K. The results of this study Suggest that, whenever feasible, one should include Molecular simulations at elevated temperatures when optimizing parameters.