In silico modeling studies of N9 -substituted harmine derivatives as potential anticancer agents: combination of ligand-based and structure-based approaches

A computational study was carried out to develop quantitative-structure activity relationship (QSAR), pharmacophore, molecular docking and molecular dynamics simulations of a series of N9-substituted harmine derivatives in order to investigate the structural factors involved in the cytotoxic activity and thus design new active derivatives. A valid 3 D-QSAR (R-2= 0.89, q(2)=0.67, R(2)pred = 0.72) and 2 D-QSAR (R-2= 0.81, q(2)=0.69, R(2)pred = 0.76) models were obtained correlating the cytotoxic activity with hydrophobic and hydrogen bond acceptor (HBA) features for 3 D-QSAR and SlogP and a_acc descriptors for 2 D-QSAR. Analysis of the selected descriptors for both models highlighted that lipophilicity and hydrogen bonding acceptor atoms remain the crucial properties and those on which cytotoxic activity depends. Also, these findings are in agreement with the characteristics of the generated pharmacophore. Furthermore, molecular docking revealed that the binding energy (-9.74 kcal/mol) and inhibition constant (0.071 mu mol) correlate with the activity of the most active compound that forms hydrophobic interactions and two hydrogen bonds with the the dual specificity tyrosine phosphorylation regulated kinase 1 A (DYRK1A). The molecular dynamics simulations revealed that the protein-ligand equilibrium is stable after 100000 fs of trajectories. Based on these results, we designed new N-9 -substituted harmine derivatives with improved properties: predicted cytotoxic activities, estimated binding energies, estimated inhibition constants and interaction modes with amino acid residues of DYRK1A, compared to the best compound in the studied dataset. Additionally, these newly designed inhibitors showed promising results in the preliminary in silico Absorption, Distribution, Metabolism, Excretion and Toxicity (ADMET) evaluations. Communicated by Ramaswamy H. Sarma