title: Lessons learned from comparing molecular dynamics engines on the SAMPL5 dataset.
abstract: 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.
authors: Department of Chemical and Biological Engineering, University of Colorado Boulder, Boulder, CO, USA. michael.shirts@colorado.edu.
authors: Michael R
authors: MR
authors: Shirts
authors: Department of Chemical Engineering, Vanderbilt University, Nashville, TN, USA.
authors: Christoph
authors: C
authors: Klein
authors: Department of Chemistry and Chemical Biology, Rutgers University, New Brunswick, NJ, USA.
authors: Jason M
authors: JM
authors: Swails
authors: Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California San Diego, La Jolla, CA, USA.
authors: Jian
authors: J
authors: Yin
authors: Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California San Diego, La Jolla, CA, USA.
authors: Michael K
authors: MK
authors: Gilson
authors: Departments of Pharmaceutical Sciences and Chemistry, University of California, Irvine, Irvine, CA, USA.
authors: David L
authors: DL
authors: Mobley
authors: Department of Chemistry and Chemical Biology, Rutgers University, New Brunswick, NJ, USA.
authors: David A
authors: DA
authors: Case
authors: Department of Chemical Engineering, University of Virginia, Charlottesville, VA, USA.
authors: Ellen D
authors: ED
authors: Zhong
mesh_terms: D001665:Binding Sites; D008024:Ligands; D008968:Molecular Conformation; D056004:Molecular Dynamics Simulation; D015394:Molecular Structure; D011485:Protein Binding; D011506:Proteins / Q000737:chemistry*; D012984:Software; D012997:Solvents / Q000737:chemistry; D013816:Thermodynamics; D014867:Water / Q000737:chemistry
chemical_list: D008024:Ligands; D011506:Proteins; D012997:Solvents; D014867:Water
keywords: Molecular dynamics; Molecular simulation; SAMPL5; Simulation validation