Abstract:Despite decades of investigations, the principal mechanisms responsible for the high affinity and specificity of proteins for key physiological cations K(+), Na(+), and Ca(2+) remain a hotly debated topic. At the core of the debate is an apparent need (or lack thereof) for an accurate description of the electrostatic response of the charge distribution in a protein to the binding of an ion. These effects range from partial electronic polarization of the directly ligating atoms to long-range effects related to partial charge transfer and electronic delocalization effects. While accurate modeling of cation recognition by metalloproteins warrants the use of quantum-mechanics (QM) calculations, the most popular approximations used in major biomolecular simulation packages rely on the implicit modeling of electronic polarization effects. That is, high-level QM computations for ion binding to proteins are desirable, but they are often unfeasible, because of the large size of the reactive-site models and the need to sample conformational space exhaustively at finite temperature. Several solutions to this challenge have been proposed in the field, ranging from the recently developed Drude polarizable force-field for simulations of metalloproteins to approximate tight-binding density functional theory (DFTB). To delineate the usefulness of different approximations, we examined the accuracy of three recent and commonly used theoretical models and numerical algorithms, namely, CHARMM C36, the latest developed Drude polarizable force fields, and DFTB3 with the latest 3OB parameters. We performed MD simulations for 30 cation-selective proteins with high-resolution X-ray structures to create ensembles of structures for analysis with different levels of theory, e.g., additive and polarizable force fields, DFTB3, and DFT. The results from DFT computations were used to benchmark CHARMM C36, Drude, and DFTB3 performance. The explicit modeling of quantum effects unveils the key electrostatic properties of the protein sites and the importance of specific ion-protein interactions. One of the most interesting findings is that secondary coordination shells of proteins are noticeably perturbed in a cation-dependent manner, showing significant delocalization and long-range effects of charge transfer and polarization upon binding Ca(2+).

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

Parameterization of Divalent Cations for Coarse-Grained Simulations

Chelator-Based Parameterization of the 12-6-4 Lennard-Jones Molecular Mechanics Potential for More Realistic Metal Ion–Protein Interactions

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

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

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

The SIRAH force field 2.0: Altius, Fortius, Citius

Modelling the interaction between silver(I) ion and proteins with the 12-6 Lennard-Jones potential: a bottom-up parameterization approach

Optimized Magnesium Force Field Parameters for Biomolecular Simulations with Accurate Solvation, Ion-Binding, and Water-Exchange Properties

Upgrading and Validation of the AMBER Force Field for Histidine and Cysteine Zinc(II)-Binding Residues in Sites with Four Protein Ligands

A New Quantum Calibrated Force Field for Zinc-Protein Complex

Taking into Account the Ion-Induced Dipole Interaction in the Nonbonded Model of Ions

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

Parameterization of a Fluctuating Charge Model for Complexes Containing 3d Transition Metals

Structural Survey of Zinc Containing Proteins and the Development of the Zinc AMBER Force Field (ZAFF).

Solution NMR refinement of a metal ion bound protein using metal ion inclusive restrained molecular dynamics methods

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

Metallicious : Automated Force-Field Parameterization of Covalently Bound Metals for Supramolecular Structures

Systematic Derivation of AMBER Force Field Parameters Applicable to Zinc-Containing Systems.

metallicious: Automated force-field parametrization of covalently bound metals for supramolecular structures

Development of Nonbonded Models for Metal Cations Using the Electronic Continuum Correction