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

Quantum Mechanics of Proteins in Explicit Water: The Role of Plasmon-Like Solute-Solvent Interactions

Including Explicit Water Molecules as Part of the Protein Structure in MM/PBSA Calculations.

Water at interface with proteins

Solvation and Desolvation Effects in Protein Folding: Native Flexibility, Kinetic Cooperativity and Enthalpic Barriers under Isostability Conditions

A quantum mechanical computational method for modeling electrostatic and solvation effects of protein

Dynamic Entity Formed by Protein and Its Hydration Water

Anomalous Protein-Protein Interactions in Multivalent Salt Solution

Decoupling Between The Translation And Rotation Of Water In The Proximity Of A Protein Molecule

Protein Sequence, Structure, Stability and Functionality

Combining Neutron Scattering, Deuteration Technique, and Molecular Dynamics Simulations to Study Dynamics of Protein and Its Surface Water Molecules

High-Precision Megahertz-to-Terahertz Dielectric Spectroscopy of Protein Collective Motions and Hydration Dynamics

Dynamical coupling between protein conformational fluctuation and hydration water: Heterogeneous dynamics of biological water

Hydrophobic and ionic-interactions in bulk and confined water with implications for collapse and folding of proteins

Quantum transport and utilization of free energy in protein $α$-helices

Quantum dynamics of electronic excitations in biomolecular chromophores: role of the protein environment and solvent

Equations for Stochastic Macromolecular Mechanics of Single Proteins: Equilibrium Fluctuations, Transient Kinetics and Nonequilibrium Steady-State

How Does Water Boost the Protein Dynamics?

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

Energy, water, and protein folding: A molecular dynamics‐based quantitative inventory of molecular interactions and forces that make proteins stable

Hydration Water Dynamics and Instigation of Protein Structural Relaxation

Protein and hydration-water dynamics are decoupled: A new model connecting dynamics and biochemical function is required