Quantum chemical study of ligand-receptor electrostatic interactions in molecular recognition of the glutamate receptor

To understand the mechanism of molecular recognition of glutamate receptors, we calculated the dipole moments, net charges, and electrostatic potentials of the agonists and antagonists of the glutamate receptors under aqueous conditions by the ab initio molecular orbital method at the HF/6-311++G(3df,2pd) level. All of the ligands had negative net charges at both termini and, consequently, formed negative electrostatic potentials at these termini. We further calculated the net charges and electrostatic potential of the S1S2 lobes of the glutamate receptor subunit GluR2 under aqueous conditions by the semiempirical PM3 molecular orbital method, and found that the ligand-binding cleft of the S1S2 lobes formed a primarily positive electrostatic potential. A strongly positive electrostatic potential was formed particularly around Arg485 in the S1 lobe. A negative electrostatic potential was observed only in a small region around Glu657 in the S2 lobe in the ligand-binding cleft. When a ligand approaches the ligand-binding cleft, it may proceed to Arg485 in the S1 lobe, due to both repulsion by the negative electrostatic potential of Glu657 as well as the attraction of the strongly positive electrostatic potential of Arg485 itself.