These are preliminary reports that have not been peer-reviewed. They should not be regarded as conclusive, guide clinical practice/health-related behavior, or be reported in news media as established information. For more information, please see our FAQs.
2 files

Reproducibility of Free Energy Calculations Across Different Molecular Simulation Software

submitted on 01.06.2018, 16:23 and posted on 04.06.2018, 22:16 by Hannes H. Loeffler, Stefano Bosisio, Guilherme Duarte Ramos Matos, Donghyuk Suh, Benoît Roux, David L. Mobley, Julien Michel

Alchemical free energy calculations are an increasingly important modern simulation technique. Contemporary molecular simulation software such as AMBER, CHARMM, GROMACS and SOMD include support for the method. Implementation details vary among those codes but users expect reliability and reproducibility, i.e. for a given molec- ular model and set of forcefield parameters, comparable free energy should be obtained within statistical bounds regardless of the code used. Relative alchemical free energy (RAFE) simulation is increasingly used to support molecule discovery projects, yet the reproducibility of the methodology has been less well tested than its absolute counter- part. Here we present RAFE calculations of hydration free energies for a set of small organic molecules and demonstrate that free energies can be reproduced to within about 0.2 kcal/mol with aforementioned codes. Achieving this level of reproducibility requires considerable attention to detail and package–specific simulation protocols, and no uni- versally applicable protocol emerges. The benchmarks and protocols reported here should be useful for the community to validate new and future versions of software for free energy calculations.


Email Address of Submitting Author


University of Edinburgh


United Kingdom

ORCID For Submitting Author


Declaration of Conflict of Interest

no conflict of interest


Read the published paper

in Journal of Chemical Theory and Computation