Skip to content
DraconiticResearch
Search Sources API
206,690 studies31 sourcesUpdated 18 min ago
PreprintCC BYChemRxiv

Water-Model Dependence of Absolute Binding Free Energies across the Bromodomain Family

  1. 1Faculty of Engineering
  2. 2Ege University

Abstract

The choice of explicit water model may affect absolute binding free energy (ABFE) predictions for bromodomains, whose conserved acetyl-lysine pockets contain structured waters. We assembled a benchmark of ten distinct bromodomain–ligand complexes with experimental isothermal titration calorimetry thermodynamics and compared four TIP-family water models, TIP3P, TIP4P, TIP5P, and TIP7P, using a common alchemical protocol. TIP3P gave the lowest point-estimate mean unsigned error (MUE), 1.56 kcal mol−1, against 1.89–2.11 kcal mol−1 for the other models, at the lowest cost. TIP7P correlated best with experiment (r = 0.73), but the experimental affinities span only 2.1 kcal mol−1, about the size of the errors of the models. The 0.55 kcal mol−1 MUE difference between TIP3P and TIP7P had a nominal 95% paired bootstrap interval that excluded zero when only the benchmark composition was resampled, but not when the calculated affinities were perturbed by their sampling-sensitivity measures or the four charged ligands were left out. No other paired difference was resolved, including after finite-size corrections for the charged ligands. The leg decomposition shows that solvent and complex contributions can compensate, so water-model effects must be judged over the complete binding cycle. With one calculation per complex and model, the higher-site TIP models gave no demonstrated accuracy benefit with this force field for 1.4–2.7 times the cost.

Related

  1. Predictable Accuracy of Martini 3 Absolute Protein-Ligand Binding Free Energies across 15 Targets
    Luis J. Walter, Tristan Bereau
    arXivPreprintCC BY 4.0Full textPhysical chemistryarXiv
  2. Molecular Mechanics-Based Symmetric Two-State Energy-Gap Estimation of Binding Free Energies
    Eiko Koizumi, Jonathan Creasey, Mohammad Alam +1
    ChemRxivPreprintLicence not statedGeneral chemistryChemRxiv
  3. Ionization Energies, Electron Affinities, Bandgaps, Exciton Binding Energies, and Polarization Energies of Orientation-Controlled Picene, [6]-Phenacene, and [7]-Phenacene Thin Films
    Rintaro Makino, Mihiro Kubo, Keiichirou Yonezawa +2
    arXivPreprintCC BY 4.0Full textPhysical chemistryCondensed matterarXiv
  4. Enhanced sampling enables binding site water reorganization in polarizable absolute binding free energy calculations
    Marharyta Blazhynska, Narjes Ansari, Louis Lagardère +1
    Communications ChemistryPeer-reviewedCC BY 4.0Full text5 citedTheoretical & computational chemistryBiological & medical physicsNatureCrossrefHALOpenAlex
  5. Site-Selective Functional Group Classification of Auger-Electron Spectra and Core-Electron Binding Energies with Convolutional Neural Networks
    Adam E. A. Fouda
    arXivPreprintCC BY 4.0Full textPhysical chemistryarXiv
BibTeX
@misc{altun2026watermodel,
  title = {{Water-Model Dependence of Absolute Binding Free Energies across the Bromodomain Family}},
  author = {Ege Altun and İpek Sarıbal and Asude Sezer and Suleyman Selim Cinaroglu},
  howpublished = {ChemRxiv},
  year = {2026},
  month = oct,
  doi = {10.26434/chemrxiv.15010882/v1},
  url = {https://chemrxiv.org/doi/full/10.26434/chemrxiv.15010882/v1},
  note = {CC BY},
}
RIS
TY  - UNPB
TI  - Water-Model Dependence of Absolute Binding Free Energies across the Bromodomain Family
AU  - Ege Altun
AU  - İpek Sarıbal
AU  - Asude Sezer
AU  - Suleyman Selim Cinaroglu
PY  - 2026
DA  - 2026/10/22
PB  - ChemRxiv
DO  - 10.26434/chemrxiv.15010882/v1
UR  - https://chemrxiv.org/doi/full/10.26434/chemrxiv.15010882/v1
AB  - The choice of explicit water model may affect absolute binding free energy (ABFE) predictions for bromodomains, whose conserved acetyl-lysine pockets contain structured waters. We assembled a benchmark of ten distinct bromodomain–ligand complexes with experimental isothermal titration calorimetry thermodynamics and compared four TIP-family water models, TIP3P, TIP4P, TIP5P, and TIP7P, using a common alchemical protocol. TIP3P gave the lowest point-estimate mean unsigned error (MUE), 1.56 kcal mol−1, against 1.89–2.11 kcal mol−1 for the other models, at the lowest cost. TIP7P correlated best with experiment (r = 0.73), but the experimental affinities span only 2.1 kcal mol−1, about the size of the errors of the models. The 0.55 kcal mol−1 MUE difference between TIP3P and TIP7P had a nominal 95% paired bootstrap interval that excluded zero when only the benchmark composition was resampled, but not when the calculated affinities were perturbed by their sampling-sensitivity measures or the four charged ligands were left out. No other paired difference was resolved, including after finite-size corrections for the charged ligands. The leg decomposition shows that solvent and complex contributions can compensate, so water-model effects must be judged over the complete binding cycle. With one calculation per complex and model, the higher-site TIP models gave no demonstrated accuracy benefit with this force field for 1.4–2.7 times the cost.
ER  -