Predictable Accuracy of Martini 3 Absolute Protein-Ligand Binding Free Energies across 15 Targets
Abstract
Alchemical free-energy calculations predict protein-ligand binding affinities with close-to-experimental accuracy, but their cost restricts them to small sets of ligands. Coarse-grained (CG) models such as Martini 3 retain a rigorous statistical-mechanical treatment at a fraction of this cost, yet their accuracy for protein-ligand binding has not been assessed systematically. Here, we compute absolute binding free energies for 172 ligands across 15 protein targets with a Martini 3 double-decoupling protocol and compare them against atomistic FEP+ and OpenFE references. Martini underbinds for every target, by 2.7 to 8.7 kcal/mol. This offset is systematic: it comprises a uniform contribution that grows with the number of charged pocket residues. A linear model in this number, which requires only the protein structure, reduces the mean absolute error from 5.7 to 1.2 kcal/mol in leave-one-out cross-validation, close to the 1.0 kcal/mol obtained with exact per-target offsets. The ranking of ligands within a target varies more widely (median Spearman ρ of 0.32, from −0.41 to 0.89) and is rationalized by the structure-activity relationship of each series: Martini ranks ligands well when their affinity is driven by hydrophobicity, π-stacking, or size, and poorly when it is driven by hydrogen bonding or charge. Both trends follow from Martini's calibration against bulk partitioning, which captures neither the dielectric environment of a binding pocket nor the directionality of hydrogen bonds. While MMGBSA retains the atomistic detail needed to rank congeneric ligands, Martini yields physically meaningful free energies, and its coarse-graining maps many atomistically distinct molecules onto one representation, making it a complementary tool for rapidly spanning chemical space.
Related
- Water-Model Dependence of Absolute Binding Free Energies across the Bromodomain Family
- Ranking versus Absolute Accuracy in Matter–Light Simulation Cost Prediction
- Finite-Field Second Hyperpolarizability via Adaptive Bayesian Quantum Phase Estimation for Early Fault-Tolerant Quantum Computers: Ground-State Energies Beyond Chemical Accuracy
- Coarse-Graining in Drug Discovery: Opportunities and Challenges for Martini 3 in Small-Molecule Binding
- A Structure-Programmable Palette of Compact Cationic Dyes Enables Predictable Organelle Addressing and Multicolor Mitochondrial Imaging.
BibTeX
@misc{walter2026predictable,
title = {{Predictable Accuracy of Martini 3 Absolute Protein-Ligand Binding Free Energies across 15 Targets}},
author = {Luis J. Walter and Tristan Bereau},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.11714},
eprint = {2610.11714},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.11714},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Predictable Accuracy of Martini 3 Absolute Protein-Ligand Binding Free Energies across 15 Targets
AU - Luis J. Walter
AU - Tristan Bereau
PY - 2026
DA - 2026/10/08
PB - arXiv
DO - 10.48550/arxiv.2610.11714
UR - https://arxiv.org/abs/2610.11714
AB - Alchemical free-energy calculations predict protein-ligand binding affinities with close-to-experimental accuracy, but their cost restricts them to small sets of ligands. Coarse-grained (CG) models such as Martini 3 retain a rigorous statistical-mechanical treatment at a fraction of this cost, yet their accuracy for protein-ligand binding has not been assessed systematically. Here, we compute absolute binding free energies for 172 ligands across 15 protein targets with a Martini 3 double-decoupling protocol and compare them against atomistic FEP+ and OpenFE references. Martini underbinds for every target, by 2.7 to 8.7 kcal/mol. This offset is systematic: it comprises a uniform contribution that grows with the number of charged pocket residues. A linear model in this number, which requires only the protein structure, reduces the mean absolute error from 5.7 to 1.2 kcal/mol in leave-one-out cross-validation, close to the 1.0 kcal/mol obtained with exact per-target offsets. The ranking of ligands within a target varies more widely (median Spearman ρ of 0.32, from −0.41 to 0.89) and is rationalized by the structure-activity relationship of each series: Martini ranks ligands well when their affinity is driven by hydrophobicity, π-stacking, or size, and poorly when it is driven by hydrogen bonding or charge. Both trends follow from Martini's calibration against bulk partitioning, which captures neither the dielectric environment of a binding pocket nor the directionality of hydrogen bonds. While MMGBSA retains the atomistic detail needed to rank congeneric ligands, Martini yields physically meaningful free energies, and its coarse-graining maps many atomistically distinct molecules onto one representation, making it a complementary tool for rapidly spanning chemical space.
ER -