Empirical formulation and parameterization of cation–π interactions for protein modeling

Cation–π interaction is comparable and as important as other main molecular interaction types, such as hydrogen bond, electrostatic interaction, van der Waals interaction, and hydrophobic interaction. Cation–π interactions frequently occur in protein structures, because six (Phe, Tyr, Trp, Arg, Lys, and His) of 20 natural amino acids and all metallic cations could be involved in cation–π interaction. Cation–π interactions arise from complex physicochemical nature and possess unique interaction behaviors, which cannot be modeled and evaluated by existing empirical equations and force field parameters that are widely used in the molecular dynamics. In this study, the authors present an empirical approach for cation–π interaction energy calculations in protein interactions. The accurate cation–π interaction energies of aromatic amino acids (Phe, Tyr, and Try) with protonated amino acids (Arg and Lys) and metallic cations (Li+, Na+, K+, and Ca2+) are calculated using B3LYP/6‐311+G(d,p) method as the benchmark for the empirical formulization and parameterization. Then, the empirical equations are built and the parameters are optimized based on the benchmark calculations. The cation–π interactions are distance and orientation dependent. Correspondingly, the empirical equations of cation–π interactions are functions of two variables, the distance r and the orientation angle θ. Two types of empirical equations of cation–π interactions are proposed. One is a modified distance and orientation dependent Lennard–Jones equation. The second is a polynomial function of two variables r and θ. The amino acid‐based empirical equations and parameters provide simple and useful tools for evaluations of cation–π interaction energies in protein interactions. © 2011 Wiley Periodicals, Inc. J Comput Chem , 2011