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

Implementation of the QUBE Force Field in SOMD for High-Throughput Alchemical Free Energy Calculations

revised on 16.04.2021, 15:40 and posted on 19.04.2021, 05:59 by Lauren Nelson, Sofia Bariami, Chris Ringrose, Joshua Horton, Vadiraj Kurdekar, Antonia Mey, Julien Michel, Daniel Cole

The quantum mechanical bespoke (QUBE) force field approach has been developed to facilitate the automated derivation of potential energy function parameters for modelling protein-ligand binding. To date the approach has been validated in the context of Monte Carlo simulations of protein-ligand complexes. We describe here the implementation of the QUBE force field in the alchemical free energy calculation molecular dynamics simulation package SOMD. The implementation is validated by demonstrating the reproducibility of absolute hydration free energies computed with the QUBE force field across the SOMD and GROMACS software packages. We further demonstrate, by way of a case study involving two series of non-nucleoside inhibitors of HIV-1 reverse transcriptase, that the availability of QUBE in a modern simulation package that makes efficient use of GPU acceleration will facilitate high-throughput alchemical free energy calculations.


Peta-5: A National Facility for Petascale Data Intensive Computation and Analytics

Engineering and Physical Sciences Research Council

Find out more...

Application of large-scale quantum mechanical simulation to the development of future drug therapies

Engineering and Physical Sciences Research Council

Find out more...



Email Address of Submitting Author


Newcastle University


United Kingdom

ORCID For Submitting Author


Declaration of Conflict of Interest

No conflict of interest