Help Privacy Policy Disclaimer
  Advanced SearchBrowse




Journal Article

Benchmarking polarizable and non-polarizable force fields for Ca2+-peptides against a comprehensive QM dataset


Hu,  Xiaojuan
NOMAD, Fritz Haber Institute, Max Planck Society;


Baldauf,  Carsten
Molecular Physics, Fritz Haber Institute, Max Planck Society;

External Resource
No external resources are shared
Fulltext (restricted access)
There are currently no full texts shared for your IP range.
Fulltext (public)

(Publisher version), 7MB

Supplementary Material (public)
There is no public supplementary material available

Amin, K. S., Hu, X., Salahub, D. R., Baldauf, C., Lim, C., & Noskov, S. (2020). Benchmarking polarizable and non-polarizable force fields for Ca2+-peptides against a comprehensive QM dataset. The Journal of Chemical Physics, 153(14): 144102. doi:10.1063/5.0020768.

Cite as: https://hdl.handle.net/21.11116/0000-0007-4ECF-A
Explicit description of atomic polarizability is critical for the accurate treatment of inter-molecular interactions by force fields (FFs) in molecular dynamics (MD) simulations aiming to investigate complex electrostatic environments such as metal-binding sites of metalloproteins. Several models exist to describe key monovalent and divalent cations interacting with proteins. Many of these models have been developed from ion–amino-acid interactions and/or aqueous-phase data on cation solvation. The transferability of these models to cation–protein interactions remains uncertain. Herein, we assess the accuracy of existing FFs by their abilities to reproduce hierarchies of thousands of Ca2+–dipeptide interaction energies based on density-functional theory calculations. We find that the Drude polarizable FF, prior to any parameterization, better approximates the QM interaction energies than any of the non-polarizable FFs. Nevertheless, it required improvement in order to address polarization catastrophes where, at short Ca2+–carboxylate distances, the Drude particle of oxygen overlaps with the divalent cation. To ameliorate this, we identified those conformational properties that produced the poorest prediction of interaction energies to reduce the parameter space for optimization. We then optimized the selected cation–peptide parameters using Boltzmann-weighted fitting and evaluated the resulting parameters in MD simulations of the N-lobe of calmodulin. We also parameterized and evaluated the CTPOL FF, which incorporates charge-transfer and polarization effects in additive FFs. This work shows how QM-driven parameter development, followed by testing in condensed-phase simulations, may yield FFs that can accurately capture the structure and dynamics of ion–protein interactions.