Abstract:Cation-pi interactions play a significant role in a host of processes eminently relevant to biology. However, polarization effects arising from the interaction of cations with aromatic moieties have long been recognized to be inadequately described by pairwise additive force fields. In the present work, we address this longstanding shortcoming through the nonbonded FIX (NBFIX) feature of the CHARMM36 force field, modifying pair-specific Lennard-Jones (LJ) parameters, while circumventing the limitations of the Lorentz-Berthelot combination rules. The potentials of mean force (PMFs) characterizing prototypical Cation-pi interactions in aqueous solutions are first determined using a hybrid quantum mechanical/molecular mechanics (QM/MM) strategy in conjunction with an importance-sampling algorithm. The LJ parameters describing the Cation-pi pairs are then optimized to match the QM/MM PMFs. The standard binding free energies of nine Cation-pi complexes, i.e., toluene, para-cresol, and 3-methyl-indole interacting with either ammonium, guanidinium, or tetramethylammonium, determined with this new set of parameters agree well with the experimental measurements. Additional simulations were carried out on three different classes of biological objects featuring Cation-pi interactions, including five individual proteins, three protein-ligand complexes, and two protein-protein complexes. Our results indicate that the description of Cation-pi interactions is overall improved using NBFIX corrections, compared with the standard pairwise additive force field. Moreover, an accurate binding free energy calculation for a protein-ligand complex containing Cation-pi interactions (2BOK) shows that using the new parameters, the experimental binding affinity can be reproduced quantitatively. Put together, the present work suggests that the NBFIX parameters optimized here can be broadly utilized in the simulation of proteins in an aqueous solution to enhance the representation of Cation-pi interactions, at no additional computational cost.

Accurate Description of Solvent-Exposed Salt Bridges with a Non-polarizable Force Field Incorporating Solvent Effects

Simulation Study of Ion Pairing in Concentrated Aqueous Salt Solutions with a Polarizable Force Field

Free Energy of Solvated Salt Bridges: a Simulation and Experimental Study.

Salt Bridge in Aqueous Solution: Strong Structural Motifs but Weak Enthalpic Effect

Electrostatic solvation energy for two oppositely charged ions in a solvated protein system: salt bridges can stabilize proteins.

Ion Solvation Thermodynamics from Simulation with a Polarizable Force Field

The multi-scale polarizable pseudo-particle solvent coarse-grained approach: From NaCl salt solutions to polyelectrolyte hydration

Small Molecule Solvation Changes Due to the Presence of Salt Are Governed by the Cost of Solvent Cavity Formation and Dispersion

Solvation-Enhanced Salt Bridges

Simulating Monovalent and Divalent Ions in Aqueous Solution Using a Drude Polarizable Force Field.

Structural Analysis of Confined Monovalent Salts: Combined Effects of Steric Hindrance, Surface Charge Representation, and Dielectric Response

Accurate Description of Cation-π Interactions in Proteins with a Nonpolarizable Force Field at No Additional Cost.

Ionic force field optimization based on single-ion and ion-pair solvation properties

Scaling Solute-Solvent Distances to Improve Solubility and Ion Paring Predictions in Rigid Ion Models

Solvent Effect on Dative and Ionic Bond Strengths: A Unified Theory from Potential Analysis and Valence-Bond Computations

Efficient Molecular Approach to Quantifying Solvent-Mediated Interactions

Force field for halide and alkali ions in water based on single-ion and ion-pair thermodynamic properties for a wide range of concentrations

Solvent effect on the UV/Vis absorption spectra in aqueous solution: The nonequilibrium polarization with an explicit representation of the solvent environment

Balanced polarizable Drude force field parameters for molecular anions: phosphates, sulfates, sulfamates, and oxides

Self-Consistent Implementation of a Solvation Free Energy Framework to Predict the Salt Solubilities of Six Alkali Halides