Molecular docking study of the binding of aminopyridines within the K+ channel

We present a molecular docking study aimed to identify the binding site of protonated aminopyridines for the blocking of voltage dependent K⁺ channels. Several active aminopyridines are considered: 2-aminopyridine, 3-aminopyridine, 4-aminopyridine, 3,4-diaminopyridine, and 4-aminoquinoleine. We apply the AutoDock force field with a lamarckian genetic algorithm, using atomic charges for the ligands derived from the electrostatic potential obtained at the B3LYP/cc-pVDZ level. We find a zone in the α-subunit of the K⁺ channel bearing common binding sites. This zone corresponds to five amino acids comprised between residuals Thr107 and Ala111, in the KcsA K⁺ channel (1J95 pdb structure). The 2-aminopyridine, 3-aminopyridine, 4-aminopyridine, and 3,4-diaminopyridine bind to the carboxylic oxygens of Thr107 and Ala111. In all cases aminopyridines are perpendicular to the axis of the pore. 4-aminoquinoleine binds to the carboxylic oxygen of Ala111. Due to its large size, the molecular plane is parallel to the axis of the pore. The charge distributions and the structures of the binding complexes suggest that the interaction is driven by formation of several hydrogen bonds. We find 2-aminopyridine, 3-aminopyridine, 4-aminopyridine, and 3,4-diaminopyridine with similar binding energy. Considering the standard error of the estimate of the AutoDock force field, this energy should lie, as a rough estimation, in the interval 3–7 kcal mol⁻¹. On the other hand, 4-aminoquinoleine seems to have a smaller binding energy. Figure Three-dimensional structure of the complex between 4-AQH⁺ and the binding sites of the K⁺ pore. Only the amino acid sequence from Thr107 to Ala111 is considered. Two different representations are included. In the left, the Thr107 position is marked. The right representation shows the CO oxygens of the peptide bond as spacefilled structures.