Lessons learned from comparing molecular dynamics engines on the SAMPL5 dataset

详细信息    查看全文

• 作者：Michael R. Shirts ; Christoph Klein…
• 关键词：Molecular dynamics ; Simulation validation ; Molecular simulation ; SAMPL5
• 刊名：Journal of Computer-Aided Molecular Design
• 出版年：2017
• 出版时间：January 2017
• 年：2017
• 卷：31
• 期：1
• 页码：147-161
• 全文大小：
• 刊物类别：Chemistry and Materials Science
• 刊物主题：Physical Chemistry; Computer Applications in Chemistry; Animal Anatomy / Morphology / Histology;
• 出版者：Springer International Publishing
• ISSN：1573-4951
• 卷排序：31

文摘

We describe our efforts to prepare common starting structures and models for the SAMPL5 blind prediction challenge. We generated the starting input files and single configuration potential energies for the host-guest in the SAMPL5 blind prediction challenge for the GROMACS, AMBER, LAMMPS, DESMOND and CHARMM molecular simulation programs. All conversions were fully automated from the originally prepared AMBER input files using a combination of the ParmEd and InterMol conversion programs. We find that the energy calculations for all molecular dynamics engines for this molecular set agree to better than 0.1 % relative absolute energy for all energy components, and in most cases an order of magnitude better, when reasonable choices are made for different cutoff parameters. However, there are some surprising sources of statistically significant differences. Most importantly, different choices of Coulomb’s constant between programs are one of the largest sources of discrepancies in energies. We discuss the measures required to get good agreement in the energies for equivalent starting configurations between the simulation programs, and the energy differences that occur when simulations are run with program-specific default simulation parameter values. Finally, we discuss what was required to automate this conversion and comparison.