Bifunctional derivatives of the alkaloid galanthamine, designed to interact with both the active site of the enzyme acetylcholinesterase (AChE) and its peripheral cation binding site, have been assayed with Torpedo californica AChE (TcAChE), and the three-dimensional structures of their complexes with the enzyme have been solved by X-ray crystallography. Differences were noted between the IC(50) values obtained for TcAChE and those for Electrophorus electricus AChE. These differences are ascribed to sequence differences in one or two residues lining the active-site gorge of the enzyme. The binding of one of the inhibitors disrupts the native conformation of one wall of the gorge, formed by the loop Trp279-Phe290. It is proposed that flexibility of this loop may permit the binding of inhibitors such as galanthamine, which are too bulky to penetrate the narrow neck of the gorge formed by Tyr121 and Phe330 as seen in the crystal structure.

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.