Drude polarizable force field for cation–π interactions of alkali and quaternary ammonium ions with aromatic amino acid side chains

Optimizing cation-π force fields for molecular dynamics studies of competitive solvation in conjugated organosulfur polymers for lithium-sulfur batteries

Lithium-sulfur (Li/S) batteries are emerging as a next-generation energy storage technology due to their high theoretical energy density and cost-effectiveness. π-Conjugated organosulfur polymers, such as poly(4-(thiophene-3-yl)benzenethiol) (PTBT), have shown promise in overcoming challenges such as the polysulfide shuttle effect by providing a conductive framework and enabling sulfur copolymerization. In these cathodes, cation-π interactions significantly influence Li+ diffusion and storage properties in π-conjugated cathodes, but classical OPLS-AA force fields fail to capture these effects. This study employs a bottom-up approach based on density functional theory (DFT) to optimize the nonbonded interaction parameters (OPLS-AA/corr.), particularly for the Li+-π interactions with the PTBT polymer. Following prior work, we used an ion-induced dipole potential to model the cation-π interactions. The impact of the solvent on the PTBT monomers was examined by computing the potential of mean force (PMF) between PTBT monomers and Li+ ions in both explicit and implicit solvents using the Boltzmann inversion of probability distributions close to room temperature. In the implicit solvent case, the magnitude of the binding free energy decreased with increasing dielectric constant, as the dominant electrostatics scaled with the dielectric constant. In contrast, in the explicit solvent case, considering the mixtures of organic solvent DME and DOL, the binding free energy shows minimal dependence on solvent composition due to the competing interaction of TBT and Li+ with the solvent molecules. However, increasing salt concentration decreases the binding free energy due to Debye-Hückel screening effects. In general, this work suggests that the optimized parameters can be widely used in the simulation of polymers in electrolytes for the Li/S battery to enhance the representation of cation-π interactions for a fixed charge force field.

Complexes of Alkaline and Ammonium Cations with Dopamine and Eumelanin Precursors: Dissecting the Role of Noncovalent Cation-π and Cation-Lone Pair (σ-Type) Interactions

Cation-π interactions and their possible competition with other noncovalent interactions (NCI) might play a key role in both dopamine- and eumelanin-based bioinspired materials. In this contribution, to unravel the delicate interplay between cation-π interactions and other possible competing forces, the configurational space of noncovalent complexes formed by dopamine or eumelanin precursors (o-benzoquinone, DHI and a semiquinone dimer) and three different cations (Na⁺, K⁺, and NH₄⁺) is sampled by means of accurate ab initio calculations. To this end, we resort to the mp2^mod method, recently validated by us for benzene-, phenol-, and catechol-cation complexes, whose computational convenience allows for an extensive exploration of the cation-molecule interaction energy surface, by sampling a total of more than 10⁴ arrangements. The mp2^mod interaction energy landscapes reveal that, besides the expected cation-π driven arrangements, for all considered molecule-cation pairs the most stable complexes are found when the cation lies within the plane containing the six-membered ring, thus maximizing the σ-type interaction with the oxygen's lone pairs. Due to the loss of aromaticity, the σ-type/cation-π strength ratio is remarkably large in o-benzoquinone, where cation-π complexes seem unlikely to be formed. The above features are shared among all considered cations but are significantly larger when considering the smaller Na⁺. Besides delivering a deeper insight onto the NCI network established by the considered precursors in the presence of ions, the present results can serve as a reference database to validate or refine lower level methods, as, for instance, the force fields employed in classical simulations.

Corrections in the CHARMM36 Parametrization of Chloride Interactions with Proteins, Lipids, and Alkali Cations, and Extension to Other Halide Anions

The nonpolarizable CHARMM force field is one of the most widely used energy functions for all-atom biomolecular simulations. Chloride is the only halide ion included in the latest version, CHARMM36m, and is used widely in simulation studies, often as an electrolyte ion but also as the biological substrate of transport proteins and enzymes. Here, we find that existing parameters systematically underestimate the interaction of Cl^- with proteins and lipids. Accordingly, when examined in solution, little to no Cl^-association can be observed with most components of the protein, including backbone, polar side chains and aromatic rings. The strength of the interaction with cationic side chains and with alkali ions is also incongruent with experimental measurements, specifically osmotic coefficients of concentrated solutions. Consistent with these findings, a 4-μs trajectory of the Cl^--specific transport protein CLC-ec1 shows irreversible Cl^- dissociation from the so-called S_cen binding site, even in a 150 mM NaCl buffer. To correct for these deficiencies, we formulate a series of pair-specific Lennard-Jones parameters that override those resulting from the conventional Lorentz-Berthelot combination rules. These parameters, referred to as NBFIX, are systematically calibrated against available experimental data as well as ab initio geometry optimizations and energy evaluations, for a wide set of binary and ternary Cl^- complexes with protein and lipid analogs and alkali cations. Analogously, we also formulate parameter sets for the other three biological halide ions, namely, fluoride, bromide, and iodide. The resulting parameters are used to calculate the potential of mean force defining the interaction of each anion and each of the protein and lipid analogues in bulk water, revealing association free energies in the range of -0.3 to -3.3 kcal/mol, with the F^- complexes being the least stable. The NBFIX corrections also preserve the Cl^- occupancy of CLC-ec1 in a second 4-μs trajectory. We posit that these optimized molecular-mechanics models provide a more realistic foundation for all-atom simulation studies of processes entailing changes in hydration, recognition, or transport of halide anions.

Accuracy of Alternate Nonpolarizable Force Fields for the Determination of Protein-Ligand Binding Affinities Dominated by Cation-π Interactions

Modifying pair-specific Lennard-Jones parameters through the nonbonded FIX (NBFIX) feature of the CHARMM36 force field has proven cost-effective for improving the description of cation-π interactions in biological objects by means of pairwise additive potential energy functions. Here, two sets of newly optimized CHARMM36 force-field parameters including NBFIX corrections, coined CHARMM36m-NBF and CHARMM36-WYF, and the original force fields, namely CHARMM36m and Amber ff14SB, are used to determine the standard binding free energies of seven protein-ligand complexes containing cation-π interactions. Compared with precise experimental measurements, our results indicate that the uncorrected, original force fields significantly underestimate the binding free energies, with a mean error of 5.3 kcal/mol, while the mean errors of CHARMM36m-NBF and CHARMM36-WYF amount to 0.8 and 2.1 kcal/mol, respectively. The present study cogently demonstrates that the use of modified parameters jointly with NBFIX corrections dramatically increases the accuracy of the standard binding free energy of protein-ligand complexes dominated by cation-π interactions, most notably with CHARMM36m-NBF.

The phenoxyl group-modulated interplay of cation-π and σ-type interactions in the alkali metal series

The interaction potential energy surfaces (IPESs) of four alkaline metal cations (Na+, K+, Rb+ and Cs+) complexed with phenol and catechol were explored by accurate ab initio calculations to investigate the interplay of different noncovalent interactions and their behavior along the alkali metal series and upon -OH substitution. Selected one-dimensional interaction energy curves revealed two different minimum energy configurations for all phenol- and catechol-metal complexes, characterized either by cation-π or σ-type interactions. For each investigated complex several two-dimensional IPES maps were also computed, exploiting the computational advantages of the MP2mod approach. The size of the alkali cation was found to play a similar role in modulating both kinds of complexes, as the interaction strength always decreases along the metal series, from Na+ to Cs+. Conversely, the number of hydroxyl substituents markedly affected cation-π complexes vs. σ-type ones. As a most relevant finding, in catechol-metal complexes the strength of cation-π interactions is around half that of the σ-type ones. It is argued that the combined effect of cation dimensions and hydroxyl substitution in catechol-Na+ complexes makes σ-type configurations remarkably more stable and easily accessible than cation-π ones. Besides shedding new light on the origin of biological phenomena connected with underwater adhesion, the quantum mechanical interaction energy database provided herein may offer a useful reference for tuning accurate force fields, suitable for molecular dynamics simulations, where environmental effects might be also taken into account.

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

Cation-π 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-π 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-π pairs are then optimized to match the QM/MM PMFs. The standard binding free energies of nine cation-π 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-π interactions, including five individual proteins, three protein-ligand complexes, and two protein-protein complexes. Our results indicate that the description of cation-π 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-π 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-π interactions, at no additional computational cost.

New Molecular-Mechanics Model for Simulations of Hydrogen Fluoride in Chemistry and Biology

Hydrogen fluoride (HF) is the most polar diatomic molecule and one of the simplest molecules capable of hydrogen-bonding. HF deviates from ideality both in the gas phase and in solution and is thus of great interest from a fundamental standpoint. Pure and aqueous HF solutions are broadly used in chemical and industrial processes, despite their high toxicity. HF is a stable species also in some biological conditions, because it does not readily dissociate in water unlike other hydrogen halides; yet, little is known about how HF interacts with biomolecules. Here, we set out to develop a molecular-mechanics model to enable computer simulations of HF in chemical and biological applications. This model is based on a comprehensive high-level ab initio quantum chemical investigation of the structure and energetics of the HF monomer and dimer; (HF)_n clusters, for n = 3-7; various clusters of HF and H₂O; and complexes of HF with analogs of all 20 amino acids and of several commonly occurring lipids, both neutral and ionized. This systematic analysis explains the unique properties of this molecule: for example, that interacting HF molecules favor nonlinear geometries despite being diatomic and that HF is a strong H-bond donor but a poor acceptor. The ab initio data also enables us to calibrate a three-site molecular-mechanics model, with which we investigate the structure and thermodynamic properties of gaseous, liquid, and supercritical HF in a wide range of temperatures and pressures; the solvation structure of HF in water and of H₂O in liquid HF; and the free diffusion of HF across a lipid bilayer, a key process underlying the high cytotoxicity of HF. Despite its inherent simplifications, the model presented significantly improves upon previous efforts to capture the properties of pure and aqueous HF fluids by molecular-mechanics methods and to our knowledge constitutes the first parameter set calibrated for biomolecular simulations.