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

Trimethyllysine Reader Proteins Exhibit Widespread Charge-Agnostic Binding via Different Mechanisms to Cationic and Neutral Ligands

Chemical Basis For The Recognition Of Trimethyllysine By Epigenetic Reader Proteins

Manipulating Cation- Interactions with Genetically Encoded Tryptophan Derivatives (vol 144, pg 6742, 2022)

Manipulating Cation- Interactions of Reader Proteins in Living Cells with Genetic Code Expansion

Binding Methylarginines and Methyllysines as Free Amino Acids: A Comparative Study of Multiple Host Classes

Mimicking the Nucleosomal Context in Peptide-Based Binders of a H3K36me Reader Increases Binding Affinity While Altering the Binding Mode

Contribution of Electrostatic CH 3 -π Interactions to Recognition of Histone Asymmetric Dimethylarginine by the SPIN1 Triple Tudor Domain

Mechanistic Insight on General Protein-Binding Ability of ATP and the Impacts of Arginine Residues

Negative Charges, Not Necessary Phosphorylation, are Required for Ligand Recognition by 14-3-3 Proteins

Manipulating Cation-π Interactions of Reader Proteins in Living Cells with Genetic Code Expansion

Polarizability Plays a Decisive Role in Modulating Association Between Molecular Cations and Anions

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

How Cations Change Peptide Structure

Insights into Ligand-Mediated Activation of an Oligomeric Ring-Shaped Gene-Regulatory Protein from Solution- and Solid-State NMR

Positive charges promote the recognition of proteins by the chaperone SlyD from Escherichia coli

Anion Binding to Ammonium and Guanidinium Hosts: Implications for the Reverse Hofmeister Effects Induced by Lysine and Arginine Residues

Consistent Picture of Phosphate-Divalent Cation Binding from Models with Implicit and Explicit Electronic Polarization.

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

Distinct Mode of Methylated Lysine-4 of Histone H3 Recognition by Tandem Tudor-Like Domains of Spindlin1.

Quantitative Interaction Proteomics and Genome-wide Profiling of Epigenetic Histone Marks and Their Readers

A consistent picture of phosphate–divalent cation binding from models with implicit and explicit electronic polarization