# BindFlow: A Free, User-Friendly Pipeline for Absolute Binding Free Energy Calculations Using Free Energy Perturbation or MM(PB/GB)SA

**Authors:** Alejandro Martínez León, Lucas Andersen, Jochen S. Hub

PMC · DOI: 10.1021/acs.jctc.5c02026 · 2026-01-07

## TL;DR

BindFlow is a free, user-friendly software for calculating binding free energies of molecules with high accuracy and lower computational cost.

## Contribution

BindFlow introduces an accessible pipeline for ABFE calculations and demonstrates that MM(PB/GB)SA can rival FEP in accuracy at lower cost.

## Key findings

- BindFlow predictions match experimental results as well as gold-standard methods.
- MM(PB/GB)SA achieves correlations comparable to FEP for some systems and force fields.
- BindFlow was validated on 139 receptor–ligand pairs across diverse targets.

## Abstract

We present BindFlow,
a Python-based software for automated
absolute
binding free energy (ABFE) calculations at the free energy perturbation
(FEP) or at the molecular mechanics Poisson–Boltzmann/generalized
Born surface area [MM­(PB/GB)­SA] level of theory. BindFlow is free,
open-source, user-friendly, and easily customizable, runs on workstations
or distributed computing platforms, and provides extensive documentation
and tutorials. BindFlow uses GROMACS as a molecular dynamics engine
and provides built-in support for the small-molecule force fields
GAFF, OpenFF, and Espaloma, as well as support for user-provided custom
force fields. We test BindFlow by computing affinities for 139 receptor–ligand
pairs, involving eight different targets, including six soluble proteins,
one membrane protein, and one nonprotein host–guest system.
We find that the agreement of BindFlow predictions with experiments
is overall similar to gold standards in the field. Interestingly,
we find that MM­(PB/GB)­SA achieves correlations that, for some systems
and force fields, approach those obtained with FEP while requiring
only a fraction of the computational cost. This study establishes
BindFlow as a validated and accessible tool for ABFE calculations.

## Full-text entities

- **Genes:** TYK2 (tyrosine kinase 2) [NCBI Gene 7297] {aka IMD35, JTK1}, GPBAR1 (G protein-coupled bile acid receptor 1) [NCBI Gene 151306] {aka BG37, GPCR19, GPR131, M-BAR, TGR5}, PPIF (peptidylprolyl isomerase F) [NCBI Gene 10105] {aka CYP3, CyP-M, Cyp-D, CypD}, IGKV2D-29 (immunoglobulin kappa variable 2D-29) [NCBI Gene 28882] {aka A2a, A2c, IGKV2D29}, PTPN1 (protein tyrosine phosphatase non-receptor type 1) [NCBI Gene 5770] {aka PTP1B}, LIG4 (DNA ligase 4) [NCBI Gene 3981] {aka LIG4S}, MCL1 (MCL1 apoptosis regulator, BCL2 family member) [NCBI Gene 4170] {aka BCL2L3, EAT, MCL1-ES, MCL1L, MCL1S, Mcl-1}, MAPK14 (mitogen-activated protein kinase 14) [NCBI Gene 1432] {aka CSBP, CSBP1, CSBP2, CSPB1, EXIP, Mxi2}, F2 (coagulation factor II, thrombin) [NCBI Gene 2147] {aka PT, RPRGL2, THPH1}
- **Diseases:** FEP (MESH:D011502)
- **Chemicals:** tecovirimat (MESH:C505045), lipids (MESH:D008055), Na+ (MESH:D012964), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (MESH:C028694), Amber99sb-ildn (-), Asp (MESH:D001224), C (MESH:D002244), SA (MESH:D000077145), Water (MESH:D014867), NaCl (MESH:D012965), Hydrogen (MESH:D006859)
- **Species:** Monkeypox virus (no rank) [taxon 10244], Homo sapiens (human, species) [taxon 9606]
- **Mutations:** A2A

## Figures

17 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12854747/full.md

---
Source: https://tomesphere.com/paper/PMC12854747