FEP-SPell-ABFE: An Open-Source Automated Alchemical Absolute Binding Free Energy Calculation Workflow for Drug Discovery

Pengfei Li,Tingting Pu,Ye Mei

DOI: https://doi.org/10.26434/chemrxiv-2024-tkvrh

2024-10-23

Abstract:The binding affinity between a drug molecule and its target, measured by absolute binding free energy (ABFE), is a crucial factor in the lead discovery phase of drug development. Recent research has highlighted the potential of in silico ABFE predictions to directly aid drug development by allowing for the ranking and prioritization of promising candidates. This paper introduces an open-source Python workflow called FEP-SPell-ABFE, designed to automate ABFE calculations with minimal user involvement. The workflow requires only three key inputs: a receptor protein structure in PDB format, candidate ligands in SDF format, and a configuration file (config.yaml) that governs both workflow and molecular dynamics simulation parameters. It produces a ranked list of ligands along with their binding free energies in Comma-Separated Values (CSV) format. The workflow leverages SLURM (Simple Linux Utility for Resource Management) for automating task execution and resource allocation across modules. A usage example and several benchmark systems for validation are provided. The FEP-SPell-ABFE workflow, along with a practical example, is publicly accessible on GitHub at https://github.com/freeenergylab/FEP-SPell-ABFE, distributed under the MIT License.

Chemistry

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: **How to accurately predict the binding affinity between drug molecules and their targets, especially in the lead compound discovery stage of drug discovery**. Specifically, the paper introduces an open - source automated workflow named FEP - SPell - ABFE, which aims to simplify the calculation process of absolute binding free energy (ABFE), thereby helping researchers screen and prioritize potential drug candidates more efficiently.

### Problem Background

In drug design, accurately predicting protein - ligand binding affinity is a fundamental goal. Binding free energy (especially absolute binding free energy, ABFE) is an important indicator for measuring the strength of the interaction between drug molecules and targets. In recent years, computer - simulation - based ABFE prediction methods have received extensive attention because they can directly assist drug development and accelerate the research process by ranking and prioritizing potential candidates.

### Existing Challenges

Although the relative binding free energy (RBFE) method is relatively mature, its application scope is limited because it depends on the chemical similarity between ligands and is usually only applicable to homologous series. In contrast, ABFE methods do not require chemical similarity between ligands and can be applied to a wider range of molecules, but are more computationally complex, involving additional steps such as introducing or removing constraints between proteins and ligands, which increase the complexity of the simulation protocol and computational requirements.

### Solution

To address these challenges, this paper proposes FEP - SPell - ABFE, an open - source Python workflow for automating ABFE calculations. This workflow only requires three key input files:
1. Receptor protein structure file (PDB format)
2. Candidate ligand structure file (SDF format)
3. Configuration file (config.yaml)

The workflow simplifies ABFE calculations in the following ways:
- **Automated task execution**: Utilize SLURM for task scheduling and resource management.
- **Reduce user intervention**: Minimize the input that users need to provide and simplify the submission process.
- **Modular design**: Include multiple modules such as topology construction, system equilibration, alchemical transformation, FEP simulation, and analysis, and each module can be run independently or in parallel.

Finally, the workflow generates a CSV file containing ligand binding free energy, helping researchers quickly evaluate and compare the binding abilities of different ligands.

### Formula Representation

The formulas involved in the paper are presented in Markdown format as follows:

\[
\Delta G_{\text{bind}}=\Delta G_{\text{int,sol}}+\Delta G_{\text{restr,dum}}-\Delta G_{\text{int,com}}+\Delta G_{\text{restr,com}}
\]

where:
- \(\Delta G_{\text{int,sol}}\): Free energy change of decoupling van der Waals and electrostatic interactions in the solvent
- \(\Delta G_{\text{restr,dum}}\): Free energy change of applying Boresch - style constraints in aqueous solution
- \(\Delta G_{\text{int,com}}\): Free energy change of decoupling van der Waals and electrostatic interactions in the binding pocket
- \(\Delta G_{\text{restr,com}}\): Free energy change of removing constraints in the binding pocket

Through these steps, the FEP - SPell - ABFE workflow can calculate the binding free energy between drug molecules and targets efficiently and accurately, thus providing strong support for drug discovery.