Stacking interaction and its role in kynurenic acid binding to glutamate ionotropic receptors

Stacking interaction is known to play an important role in protein folding, enzyme-substrate and ligand-receptor complex formation. It has been shown to make a contribution into the aromatic antagonists binding with glutamate ionotropic receptors (iGluRs), in particular, the complex of NMDA receptor NR1 subunit with the kynurenic acid (KYNA) derivatives. The specificity of KYNA binding to the glutamate receptors subtypes might partially result from the differences in stacking interaction. We have calculated the optimal geometry and binding energy of KYNA dimers with the four types of aromatic amino acid residues in Rattus and Drosophila ionotropic iGluR subunits. All ab initio quantum chemical calculations were performed taking into account electron correlations at MP2 and MP4 perturbation theory levels. We have also investigated the potential energy surfaces (PES) of stacking and hydrogen bonds (HBs) within the receptor binding site and calculated the free energy of the ligand-receptor complex formation. The energy of stacking interaction depends both on the size of aromatic moieties and the electrostatic effects. The distribution of charges was shown to determine the geometry of polar aromatic ring dimers. Presumably, stacking interaction is important at the first stage of ligand binding when HBs are weak. The freedom of ligand movements and rotation within receptor site provides the precise tuning of the HBs pattern, while the incorrect stacking binding prohibits the ligand-receptor complex formation.