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+).

Benchmarking biomolecular force field‐based Zn2+ for mono‐ and bimetallic ligand binding sites

A Comparison of Bonded and Nonbonded Zinc(II) Force Fields with NMR Data

A Transferable Non-bonded Pairwise Force Field to Model Zinc Interactions in Metalloproteins.

Comparing Non-Bonded Metal Ion Models in the Divalent Cation Binding Protein PsaA.

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

Development of a Fluctuating Charge Model for Zinc-Containing Metalloproteins

Comparison of Empirical Zn2+ Models in Protein–DNA Complexes

Protocols for metallo- and serine-β-lactamase free energy predictions: insights from cross-class inhibitors

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

Simulating Metal-Imidazole Complexes

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

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

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

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

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

Upgraded AMBER Force Field for Zinc-Binding Residues and Ligands for Predicting Structural Properties and Binding Affinities in Zinc-Proteins

A New Quantum Calibrated Force Field for Zinc-Protein Complex

Benchmark Data for Interactions in Zeolite Model Complexes and Their Use for Assessment and Validation of Electronic Structure Methods

Prediction of structures of zinc-binding proteins through explicit modeling of metal coordination geometry.

Rational Design of Nonbonded Point Charge Models for Divalent Metal Cations with Lennard-Jones 12-6 Potential

Interactions of Zn(II) with Single and Multiple Amino Acids. Insights from Density Functional and Ab Initio Calculations