Abstract:It is challenging to parameterize the force field for calcium ions (Ca2+) in calcium-binding proteins because of their unique coordination chemistry that involves the surrounding atoms required for stability. In this work, we observed wide variation in Ca2+ binding loop conformations of the Ca2+-binding protein calmodulin (CaM), which adopts the most populated ternary structures determined from the MD simulations, followed by ab initio quantum mechanical (QM) calculations on all twelve amino acids in the loop that coordinate Ca2+ in aqueous solution. Ca2+ charges were derived by fitting to the electrostatic potential (ESP) in the context of a classical or polarizable force field (PFF). We discovered that the atomic radius of Ca2+ in conventional force fields is too large for the QM calculation to capture the variation in the coordination geometry of Ca2+ in its ionic form, leading to unphysical charges. Specifically, we found that the fitted atomic charges of Ca2+ in the context of PFF depend on the coordinating geometry of electronegative atoms from the amino acids in the loop. Although nearby water molecules do not influence the atomic charge of Ca2+, they are crucial for compensating for the coordination of Ca2+ due to the conformational flexibility in the EF-hand loop. Our method advances the development of force fields for metal ions and protein binding sites in dynamic environments.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: how to accurately parameterize the force field of calcium ions (Ca²⁺) in calcium - binding proteins, especially considering their unique coordination chemistry and the influence of surrounding atoms on their stability. Specifically, the study focuses on the conformational changes of the Ca²⁺ - binding loop in calmodulin (CaM) and their impact on the Ca²⁺ charge. ### Main problems: 1. **Changes in the conformation of the Ca²⁺ - binding loop**: - The binding loop (EF - hand motif) of Ca²⁺ in calcium - binding proteins has various conformational changes, which affect the coordination geometry of Ca²⁺. - Through molecular dynamics (MD) simulations, the authors found that the conformation of the Ca²⁺ - binding loop shows significant differences in different environments (such as Ca²⁺ - retaining and - releasing environments). 2. **Limitations in traditional force fields**: - The atomic radius of Ca²⁺ used in traditional force fields is too large, resulting in the inability of quantum mechanics (QM) calculations to accurately capture the coordination geometry changes of Ca²⁺, thus producing unreasonable charge values. - These unreasonable charge values will affect the dynamic study of Ca²⁺ - protein interactions. 3. **Importance of polarization effects**: - In classical or polarizable force fields (PFF), the charge of Ca²⁺ depends on the coordination geometry of electronegative atoms from amino acid rings. - Although nearby water molecules do not affect the atomic charge of Ca²⁺, they are crucial for compensating the coordination of Ca²⁺, especially in terms of the conformational flexibility of the EF - hand loop. ### Solutions: - **Quantum mechanics calculations**: - The authors used the ab initio quantum mechanics (ab initio QM) method to calculate all amino acids involved in Ca²⁺ coordination to determine a more accurate charge distribution. - By adjusting the vdW radius of Ca²⁺ (from 1.7 Å to 1.0 Å), the problem of unreasonable charge values was solved, and a more realistic charge distribution was obtained. - **Molecular dynamics simulations**: - MD simulations were used to generate a large number of different conformations of the Ca²⁺ - binding loop, and representative structures were selected from them for QM calculations. - Through neural network clustering analysis, the most representative conformations were further screened out for subsequent calculations. - **Polarizable force field development**: - The i_RESP method was used to fit partial atomic charges, taking into account the polarization effect, making the charge distribution more reasonable. - It was verified that the theoretical calculations at the B3LYP/SVP level are highly consistent with the results at the MP2/aug - cc - pVTZ level, ensuring the accuracy of the calculations. ### Conclusion: This study, by combining MD simulations and QM calculations, has solved the problem of Ca²⁺ charge parameterization in traditional force fields and provided a new method for developing more accurate force fields for metal ions and protein - binding sites. This method is not only applicable to Ca²⁺, but can also be extended to the study of other divalent ions (such as Mg²⁺).

Determining the atomic charge of calcium ion requires the information of its coordination geometry in an EF-hand motif

Determining the atomic charge of calcium ion in a calcium-binding protein requires the information of its dynamic coordination geometry in aqueous solution

Binding Energy and Free Energy of Calcium Ion to Calmodulin EF-Hands with the Drude Polarizable Force Field

A New Quantum Calibrated Force Field for Zinc-Protein Complex

Study of interactions between metal ions and protein model compounds by energy decomposition analyses and the AMOEBA force field.

Coarse-Grained Modeling and Molecular Dynamics Simulations of Ca2+-Calmodulin

Quantum effects in cation interactions with first and second coordination shell ligands in metalloproteins

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

Identification of the critical structural determinants of the EF-hand domain arrangements in calcium binding proteins.

Calcium ions in aqueous solutions: Accurate force field description aided by ab initio molecular dynamics and neutron scattering.

New force field for calcium binding sites in annexin–membrane complexes

Calcium Ion Binding to the Mutants of Calmodulin: A Structure-Based Computational Predictive Model of Binding Affinity Using a Charge Scaling Approach in Molecular Dynamics Simulation

Many-body effect determines the selectivity for Ca2+ and Mg2+ in proteins

Extended magnesium and calcium force field parameters for accurate ion-nucleic acid interactions in biomolecular simulations.

The Effect of Macromolecular Crowding, Ionic Strength and Calcium Binding on Calmodulin Dynamics

Coarse-Grained Parameters for Divalent Cations within the SIRAH Force Field

Coarse-grained modeling of the calcium, sodium, magnesium and potassium cations interacting with proteins

Chemistry-informed Machine Learning Explains Calcium-binding Proteins Fuzzy Shape for Communicating Changes in the Atomic States of Calcium Ions

Improved model of hydrated calcium ion for molecular dynamics simulations using classical biomolecular force fields

Interaction of divalent cations and amino acids in bulk water: Molecular simulations with neural network potentials

Force fields for monovalent and divalent metal cations in TIP3P water based on thermodynamic and kinetic properties