To consider whether existing molecular force fields can adequately reproduce cation−π interactions without adding special interaction terms, theoretical calculations with geometry optimization were performed on three configurations of tetramethylammonium (TMA) interacting via one, two, or three N-methyl groups with a benzene ring, by use of density-functional theory (DFT) methods B3LYP/6-31G* and B3LYP/6-311G**, ab initio method MP2/6-31G*, and molecular mechanic methods EFF, Tinker's Amber and MM3. Only the first configuration was found to be stable from the DFT and MP2 results, and its geometry was found to be highly flexible. ESP CHELPG charges estimated from the DFT and MP2 calculations were used to modify the atomic charges of the force fields employed in the molecular mechanics calculations to improve agreement with the BSSE-corrected binding energies deduced from the DFT and MP2 results. After this modification, the molecular mechanics results were found to be in good agreement with those obtained by DFT and MP2, without a requirement to add any additional terms to the force fields. This was confirmed by comparing the energy profiles of the complex as benzene was moved away from TMA in 0.2 Å intervals. Hence it is possible to use existing force fields to represent cation−π interactions by a simple adjustment of certain partial atomic charge parameters. In this context, we discuss the high flexibility of the cation−π interactions in the framework of molecular recognition in biological systems.
Read full abstract