Accurate modeling of ion-protein binding in solution: molecular mechanics force field innovation and validation on Mg2+
5 affiliations
- 1Department of Molecular Biosciences, University of South Florida, 4202 E. Fowler Ave., Tampa, USA, FL-33620
- 2Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Műegyetem rkp. 3., Budapest, Hungary, H-1111
- 3HUN-REN-BME Quantum Chemistry Research Group, Műegyetem rkp. 3., Budapest, Hungary, H-1111
- 4MTA-BME Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., Budapest, Hungary, H-1111
- 5Department of physics, University of South Florida, 4202 E. Fowler Ave., Tampa, Florida, United States, 33620
Abstract
Multivalent ions bind to proteins to control many biological processes. Molecular simulations can effectively supplement experiments to resolve underlying mechanisms, but suffer from reliability and feasibility issues in describing ion-protein binding in the condensed phase. Here we present a molecular mechanics based generalizable solution to this long-standing challenge. Using the divalent Mg 2+ as a model, we demonstrate that our approach yields unprecedented accuracy in describing its binding to proteins in solution – compared to high level quantum mechanics, the mean error in local interaction energies is 0.5%, and the binding modes and binding free energies match experiment. Our study also shows that key to achieving this accuracy requires a delicate balance between the large repulsive and attractive forces at short distances near ions for which explicit modeling of polarization is essential. Overall, our work lays foundation for future development of molecular simulation methods for ions and opens the way to quantitatively study Mg-dependent mechanisms in hundreds of enzymes and proteins.
Related
- Optical Characterization Measurements of Commercial Polyvinyltoluene Scintillators for Accurate Geant4 Detector Modeling
- DNNsolver: Accurate and Efficient Deep Learning Modeling of Polarization-Sensitive Diffractive Neural Networks
- Neutrino Direct Simulation Monte Carlo: Accurate modeling of out-of-equilibrium decaying cosmological relics
- Modeling Transport and Mechanics in Anion Exchange Membranes to Optimize Continuous CO2 Separations
- Accurate Molecular Size Determination by Diffusion Ordered NMR Spectroscopy Based on an Improved Diffusion Model
BibTeX
@misc{schillaci2026accurate,
title = {{Accurate modeling of ion-protein binding in solution: molecular mechanics force field innovation and validation on Mg2+}},
author = {Giorgio Schillaci and Julian Mélendez Delgado and Matthew Saunders and Irene Arias Gutierrez and Péter R. Nagy and Sameer Varma},
howpublished = {ChemRxiv},
year = {2026},
month = sep,
doi = {10.26434/chemrxiv.15009072/v1},
url = {https://chemrxiv.org/doi/full/10.26434/chemrxiv.15009072/v1},
note = {CC BY},
}
RIS
TY - UNPB
TI - Accurate modeling of ion-protein binding in solution: molecular mechanics force field innovation and validation on Mg2+
AU - Giorgio Schillaci
AU - Julian Mélendez Delgado
AU - Matthew Saunders
AU - Irene Arias Gutierrez
AU - Péter R. Nagy
AU - Sameer Varma
PY - 2026
DA - 2026/09/18
PB - ChemRxiv
DO - 10.26434/chemrxiv.15009072/v1
UR - https://chemrxiv.org/doi/full/10.26434/chemrxiv.15009072/v1
AB - Multivalent ions bind to proteins to control many biological processes. Molecular simulations can effectively supplement experiments to resolve underlying mechanisms, but suffer from reliability and feasibility issues in describing ion-protein binding in the condensed phase. Here we present a molecular mechanics based generalizable solution to this long-standing challenge. Using the divalent Mg 2+ as a model, we demonstrate that our approach yields unprecedented accuracy in describing its binding to proteins in solution – compared to high level quantum mechanics, the mean error in local interaction energies is 0.5%, and the binding modes and binding free energies match experiment. Our study also shows that key to achieving this accuracy requires a delicate balance between the large repulsive and attractive forces at short distances near ions for which explicit modeling of polarization is essential. Overall, our work lays foundation for future development of molecular simulation methods for ions and opens the way to quantitatively study Mg-dependent mechanisms in hundreds of enzymes and proteins.
ER -