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 alpha-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(-1). On the other hand, 4-aminoquinoleine seems to have a smaller binding energy.

No abstract

The cannabinoid receptors (CB1/CB2) and the T-type calcium channels are involved in disorders associated with both physiological pain and depressive behaviors. Valuable pharmacological species carbazole derivatives such as the NMP-4, NMP-7, and NMP-181 (Neuro Molecular Production) regulate both biological entities. In this work, DFT calculations were performed to characterize theoretically their structural and chemical reactivity properties using the BP86/cc-pVTZ level of theory. The molecular orbital contributions and the chemical reactivity analysis reveal that a major participation of the carbazole group is in the donor-acceptor interactions of the NMP compounds. The DFT analysis on the NMP compounds provides insights into the relevant functional groups involved during the ligand-receptor interactions. Molecular docking analysis is used to reveal possible sites of interaction of the NMP compounds with the Cav3.2 calcium channel. The interaction energy values and reported experimental evidence indicate that the site denominated as “Pore-blocking”, which is formed mainly by hydrophobic residues and the T586 residue, is a probable binding site for the NMP compounds.

Regulation of cellular excitability and oscillatory behavior of resting membrane potential in nerve cells are largely mediated by the low-voltage activated T-type calcium channels. This calcium channel family is constituted by three isoforms, namely, Ca V 3.1, Ca V 3.2, and Ca V 3.3, that are largely distributed in the nervous system and other parts of the body. Dysfunction of T-type calcium channels is associated with a wide range of pathophysiologies including epilepsy, neuropathic pain, cardiac problems, and major depressive disorders. Due to their pharmacological relevance, finding molecular agents able to modulate the channel's function may provide therapeutic means to ameliorate their related disorders. Here we used electrophysiological experiments to show that genistein, a canonical tyrosine kinase inhibitor, reduces the activity of the human Ca V 3.3 channel in a concentration-dependent manner. The inhibitory effect of genistein is independent of tyrosine kinase modulation and does not affect the voltage-dependent gating of the channel. Subsequently, we used computational methods to identify plausible molecular poses for the interaction of genistein and the Ca V 3.3 channel. Starting from different molecular poses, we carried out all-atom molecular dynamics (MD) simulations to identify the interacting determinants for the Ca V 3.3/genistein complex formation. Our extensive (microsecond-length) simulations suggest specific binding interactions that seem to stabilize the protein/inhibitor complex. Furthermore, our results from the unbiased MD simulations are in good agreement with the recently solved cryoelectron microscopy structure of the Ca V 3.1/Z944 complex in terms of both the location of the ligand binding site and the role of several equivalent amino acid residues. Proposed interacting complex loci were subsequently tested and corroborated by electrophysiological experiments using another naturally occurring isoflavone derivative, daidzein. Thus, by using a combination of in vitro and in silico techniques, we have identified interacting determinants relevant to the Ca V 3.3/genistein complex formation and propose that genistein directly blocks the function of the human Ca V 3.3 channel as a result of such interaction. Specifically, we proposed that a combination of polar interactions involving the three hydroxyl groups of genistein and an aromatic interaction with the fused rings are the main binding interactions in the complex formation. Our results pave the way for the rational development of improved and novel low-voltage activated T-type calcium channel inhibitors.

No abstract