Predicting single mutation effects on binding affinity at protein-protein interfaces based on MMPBSA calculations
Abstract
Accurate prediction of mutation-induced changes in protein-protein affinity remains a central challenge in computational biophysics and protein engineering. Here we present a physics-based scoring method for predicting the effects of single amino acid substitutions on protein-protein and protein-peptide binding affinity. The method is based on MMPBSA energy calculations and implemented in a user-friendly, publicly available workflow (https://github.com/Hoecker-Lab/mmpbsa_scoring). To derive the scoring function, MMPBSA energy terms were fitted using a large and diverse dataset of experimentally measured binding affinity changes upon point mutation. The terms were then evaluated on independent systems. Specifically, we assessed the ability to predict mutational effects at the interface of SARS-CoV-2 spike protein with its cellular receptor, as well as computationally predicted Armadillo-repeat protein-peptide complexes. The approach performs well across diverse targets and shows promising predictive performance on both experimentally determined and computationally predicted structures. It matches or exceeds the accuracy of established physics-based protein design software and state-of-the-art deep learning-based scoring functions when applied to data outside their training domain. The method provides a computationally efficient and interpretable framework for prioritizing interface mutations in protein engineering and binding affinity optimization.
Related
- MPFIL: Multi-Level Protein Feature Integration for Predicting Mutation-Induced Protein-Protein Binding Affinity Changes
- A hyper-IgM2 AID mutation limits IgH break formation but preserves affinity maturation
- Benchmarking Large Language Models for Assay-Aware Drug-Target Binding Affinity Prediction
- Leakage-controlled benchmarking reveals generalization limits of deep learning for protein-ligand binding affinity prediction
- Expanding gmx_MMPBSA capabilities for end-state binding free energy calculations with GROMACS and AMBER
BibTeX
@misc{noske2026predicting,
title = {{Predicting single mutation effects on binding affinity at protein-protein interfaces based on MMPBSA calculations}},
author = {J. Noske and M. Janzen and T. Lepoivre and B. Höcker},
howpublished = {bioRxiv},
year = {2026},
month = sep,
doi = {10.64898/2026.09.22.753561},
url = {https://www.biorxiv.org/content/10.64898/2026.09.22.753561v1},
note = {CC BY},
}
RIS
TY - UNPB
TI - Predicting single mutation effects on binding affinity at protein-protein interfaces based on MMPBSA calculations
AU - J. Noske
AU - M. Janzen
AU - T. Lepoivre
AU - B. Höcker
PY - 2026
DA - 2026/09/22
PB - bioRxiv
DO - 10.64898/2026.09.22.753561
UR - https://www.biorxiv.org/content/10.64898/2026.09.22.753561v1
AB - Accurate prediction of mutation-induced changes in protein-protein affinity remains a central challenge in computational biophysics and protein engineering. Here we present a physics-based scoring method for predicting the effects of single amino acid substitutions on protein-protein and protein-peptide binding affinity. The method is based on MMPBSA energy calculations and implemented in a user-friendly, publicly available workflow (https://github.com/Hoecker-Lab/mmpbsa_scoring). To derive the scoring function, MMPBSA energy terms were fitted using a large and diverse dataset of experimentally measured binding affinity changes upon point mutation. The terms were then evaluated on independent systems. Specifically, we assessed the ability to predict mutational effects at the interface of SARS-CoV-2 spike protein with its cellular receptor, as well as computationally predicted Armadillo-repeat protein-peptide complexes. The approach performs well across diverse targets and shows promising predictive performance on both experimentally determined and computationally predicted structures. It matches or exceeds the accuracy of established physics-based protein design software and state-of-the-art deep learning-based scoring functions when applied to data outside their training domain. The method provides a computationally efficient and interpretable framework for prioritizing interface mutations in protein engineering and binding affinity optimization.
ER -