Automated Adaptive Absolute Binding Free Energy Calculations

Julien Michel,Finlay Clark,Daniel Cole,Graeme Robb
DOI: https://doi.org/10.26434/chemrxiv-2024-3ft7f
2024-06-18
Abstract:Alchemical absolute binding free energy (ABFE) calculations have substantial potential in drug discovery, but are often prohibitively computationally expensive. To unlock their potential, efficient automated ABFE workflows are required to reduce both computational cost and human intervention. We present a fully-automated ABFE workflow based on the automated selection of λ windows, the ensemble-based detection of equilibration, and the adaptive allocation of sampling time based on inter-replicate statistics. We find that the automated selection of intermediate states with consistent overlap is rapid, robust, and simple to implement. Robust detection of equilibration is achieved with a paired t-test between the free energy estimates at initial and final portions of a an ensemble of runs. We determine reasonable default parameters for all algorithms and show that the full workflow produces equivalent results to a non-adaptive scheme over a variety of test systems, while often accelerating equilibration. Our complete workflow is implemented in the open-source package A3FE (https://github.com/michellab/a3fe).
Chemistry
What problem does this paper attempt to address?
This paper mainly discusses the problem of automating the adaptation of Absolute Binding Free Energy (ABFE) calculations. In drug discovery, ABFE calculations can predict the binding affinity between molecules and targets, but they are often computationally expensive. To overcome this challenge, a fully automated ABFE workflow is proposed in this study, which is based on the selection of intermediate states, cluster-based equilibrium detection, and adaptive sampling time allocation based on replica statistics. The study found that selecting intermediate states with consistent overlap is fast, robust, and easy to implement. Balanced and robust detection is achieved by performing paired t-tests between the initial and final parts of a series of runs. The study determined the default parameters for all algorithms and demonstrated that the complete workflow produces comparable results to non-adaptive schemes on various test systems while frequently accelerating equilibration. This comprehensive workflow is implemented in the open-source package A3FE. In summary, the paper attempts to address how to improve the efficiency of ABFE calculations, reduce computational costs and human intervention, in order to predict the binding affinity between drugs and targets more effectively.