Abstract
Cation–π interactions involving the tetramethylammonium motif are common in biological systems, where they play essential roles in membrane protein function, DNA expression regulation, and protein folding. However, accurately modeling cation–π interactions, where electronic polarization is crucial, presents a significant computational challenge, especially in large biomolecular systems. This study applies a physically justified electronic continuum correction (ECC) to the CHARMM36 force field, scaling ionic charges by a factor of 0.75 to effectively account for electronic polarization without adding computational overhead. While not specifically designed for cation–π interactions, this approach significantly improves predictions of the structure of an aqueous tetramethylammonium–pyridine complex compared to neutron diffraction data. These results, along with computational predictions for the structure of the aqueous tetramethylammonium–phenol complex, highlight the potential of ECC as a versatile method for improving the description of cation–π interactions in biomolecular simulations.
Supplementary materials
Title
Supporting Information
Description
Additional metrics with TIP4P water model and discussion about the roles of both the cation and the anion.
Actions
Supplementary weblinks
Title
Molecular Dynamics Simulations Set Up and Neutron Scattering Data
Description
Additional files for simulations and neutron scattering data along with its analysis.
Actions
View