Document Type : Original Research Article


Department of Chemistry, Faculty of Physical Sciences, Ahmadu Bello University, P.M.B. 1044, Zaria, Kaduna State, Nigeria


The current research was conducted as part of the anti-leishmanial drug discovery effort towards new drug molecules with attributes that overcome the limitations of existing therapies. This work utilizes a combined approach of Quantitative Structure-Activity Relationship (QSAR), virtual docking screening, and pharmacokinetics analysis to design some novel 2,6-diarylidene cyclohexanone analogs using ligand-based drug design methods, while also performing docking investigation, drug-likeness analysis, and Molecular Dynamic (MD) simulation to evaluate their anti-leishmanial potential. Some crucial parameters were calculated for the built QSAR model, including R2 = 0.7827, R2adj = 0.7206, Q2cv = 0.6414, and R2test = 0.8539, which indicate an acceptable QSAR model. The combined results of QSAR, docking, and pharmacokinetics analysis suggested compound 1 as the template. The Six (6) newly designed analogs possessed higher binding scores than the reference drug Pentamidine in the order; 1a (-10.2 kcal/mol) > 1e (-9.6) > 1d (-9.4) > 1c (-9.2) > Template (-9.1) > 1f (-9) > 1b (-8.5) > Pentamidine (-6.9 kcal/mol), while their predicted pIC50 followed the order; 1e (8.7321) > 1c (7.6772) > 1f (7.1602) > 1a (6.8289) > 1d (6.7738) > 1b (6.5772) > Template (5.3824). The results of the drug-likeness testing suggest 1 and the new analogs (especially 1a) as being orally bioavailable with excellent pharmacokinetic profiles. These molecules equally showed good pharmacological interactions with the receptor, Pyridoxal kinase (PDB: 6K91).  In addition, the MD simulation results confirmed the stability and rigidity of 1_6K91 and 1a_6K91. Therefore, the new analogs could be considered as potent anti-leishmanial inhibitors.

Graphical Abstract

Combined QSAR Modeling, Molecular Docking Screening, and Pharmacokinetics Analyses for the Design of Novel 2, 6-Diarylidene Cyclohexanone Analogs as Potent Anti-Leishmanial Agents


  • The high limitations of existing therapies resulting in pending cases of leishmanial diseases necessitated the study.


  • A computer-aided design of novel 2, 6-diarylidene cyclohexanone analogs as anti-leishmanial agents; 2-D QSAR, molecular docking, pharmacokinetics, and MD simulation were conducted to investigate the anti-leishmanial activities, binding interaction pattern, drug-likeness properties, and the stability of protein-ligand complexes of the newly designed analogs.


  • The resulting docking scores were used to identify the best protein-ligand interaction pairs, while the docked poses were post-screened and analyzed for helpful information on binding/pharmacological interactions.


  • The generated 2-D QSAR data were subjected to some relevant validation tests in order to qualify the model.


  • Also, the newly designed analogs were subjected to pharmacokinetic analysis in order to ascertain their drug-likeness properties, while the protein-ligand complexes of 1 and 1a were analyzed for binding stability using the result of MD simulation.


  • The combined results of the study lead to the identification of compound 1a as a potential drug candidate for the treatment of leishmanial infections.


Main Subjects

[1]          L.D.C. Clementino, G.F.S. Fernandes, I.M. Prokopczyk, W.C. Laurindo, D. Toyama, B.P. Motta, Design, synthesis and biological evaluation of N-oxide derivatives with potent in vivo antileishmanial activity. PLoS ONE, 16 (2021) e0259008.
[2] A. Upadhyay, P. Chandrakar, S. Gupta, N. Parmar, S.K. Singh, M. Rashid et al., Synthesis, biological evaluation, structure−activity relationship, and mechanism of action studies of quinoline−metronidazole derivatives against experimental visceral Leishmaniasis. J. Med. Chem., 62 (2019) 5655−5671.
[3] M.D. Yousuf, D. Mukherjee, A. Pal, S. Dey, S. Mandal, C. Pal et al., Synthesis and biological evaluation of Ferrocenylquinoline as a potential antileishmanial agent. Chem. Med. Chem., 10 (2015) 546-554. 
[4] F.A. Ugbe, G.A. Shallangwa, A. Uzairu, I. Abdulkadir, A combined 2-D and 3-D QSAR modeling, molecular docking study, design, and pharmacokinetic profiling of some arylimidamide-azole hybrids as superior L. donovani inhibitors. Bull. Natl. Res. Cent., 46 (2022) 189.
[5] C.C.B. Brito, H.V.C. da Silva, D.J. Brondani, A.R. de Faria, R.M. Ximenes, I.M. da Silva et al., Synthesis and biological evaluation of thiazole derivatives as LbSOD inhibitors. Journal of Enzyme Inhibition and Medicinal Chemistry, 34 (2019) 333-342.
[6] F.A. Ugbe, G.A. Shallangwa, A. Uzairu, I. Abdulkadir, Theoretical modeling and design of some pyrazolopyrimidine derivatives as Wolbachia inhibitors, targeting lymphatic flariasis and onchocerciasis. In Silico Pharmacol., 10 (2022) 8.
[7] Y. Fan, Y. Lu, X. Chen, B. Tekwani, X. Li, Y. Shen, Anti-Leishmanial and cytotoxic activities of a series of maleimides: synthesis, biological evaluation and structure-activity relationship.  Molecules, 23 (2018) 2878.
[8] F.A. Ugbe, G.A. Shallangwa, A. Uzairu, I. Abdulkadir, Theoretical activity prediction, structure-based design, molecular docking and pharmacokinetic studies of some maleimides against Leishmania donovani for the treatment of leishmaniasis. Bull. Natl. Res. Cent., 46 (2022) 92.
[9] S.E. Adeniji, D.E. Arthur, M. Abdullahi, A. Abdullahi, F.A. Ugbe, Computer-aided modeling of triazole analogues, docking studies of the compounds on DNA gyrase enzyme and design of new hypothetical compounds with efficient activities. J. Biomol. Struct. Dyn., 2020 (2020).
[10] F.A. Ugbe, G.A. Shallangwa, A. Uzairu, I. Abdulkadir, Activity modeling, molecular docking and pharmacokinetic studies of some boron-pleuromutilins as anti-wolbachia agents with potential for treatment of filarial diseases. Chemical Data Collections, 36 (2021) 100783.
[11] L. Abdel-Ilah, E. Veljovic, L. Gurbeta, A. Badnjevic, Applications of QSAR study in drug design. International Journal of Engineering Research and Technology, 6 (6) (2017).
[12] M.T. Ibrahim, A. Uzairu, G.A. Shallangwa, S. Uba, Lead identification of some anti-cancer agents with prominent activity against Non-small Cell Lung Cancer (NSCLC) and structure-based design. Chemistry Africa, 3 (2020) 1023–1044.
[13] H.A. Lawal, A. Uzairu, S. Uba, QSAR, molecular docking studies, ligand-based design and pharmacokinetic analysis on Maternal Embryonic Leucine Zipper Kinase (MELK) inhibitors as potential anti-triple-negative breast cancer (MDA-MB-231cell line) drug compounds. Bulletin of the National Research Centre, 45 (90) (2021).
[14] C.A. Lipinski, F. Lombardo, B.W. Dominy, P.J. Feeney, Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Advanced Drug Delivery Reviews, 46 (2001) 3-26.
[15] E.I. Edache, A. Uzairu, P.A. Mamza, G.A. Shallangwa, Theoretical investigation of the cooperation of           iminoguanidine with the enzymes-binding domain of Covid-19 and bacterial lysozyme inhibitors  and their pharmacokinetic properties. J. Mex. Chem. Soc., 66(4) (2022) 513-542.
[16] S. Are, S. Gatreddi, P. Jakkula, I.A. Qureshi, Structural attributes and substrate specificity of pyridoxal kinase from Leishmania donovani. Int. J. Biol. Macromol., 152 (2020) 812-827.
[17] V. Kumar, M. Sharma, B.R. Rakesh, C.K. Malik, S. Neelagiri, K.B. Neerupudi et al., Pyridoxal kinase: a vitamin B6 salvage pathway enzyme from Leishmania donovani. Int. J. Biol. Macromol., 119 (2018) 320-334.
[18] Z.U. Din, M.A. Trapp, L. Soman de Medeiros, D. Lazarin-Bidóia, F.P. Garcia, F. Peron et al., Symmetrical and unsymmetrical substituted 2,5-diarylidene cyclohexanones as anti-parasitic compounds.  European Journal of Medicinal Chemistry, 155 (2018) 596–608.
[19] H. Chandru, A.C. Sharada, B.K. Bettadaiah, C.S.A. Kumar, K.S. Rangappa, J.K. Sunila, In vivo growth inhibitory and anti-angiogenic effects of synthetic novel dienone cyclopropoxy curcumin analogs on mouse Ehrlich ascites tumor. Bioorg Med Chem., 15 (2007) 7696–703.
[20] X. Wang, H. Dong, Q. Qin, QSAR models on aminopyrazole-substituted resorcylate compounds as Hsp90 inhibitors. J. Comp. Sci. & Eng., 48 (2020) 1146-1156.
[21] R.W. Kennard, L.A. Stone, Computer aided design of experiments. Technometrics, 11 (1969) 137–148.
[22] D.E. Arthur, A. Uzairu, P. Mamza, S.E. Abechi, G.A. Shallangwa, Activity and toxicity modeling of some NCI selected compounds against leukemia P388ADR cell line using genetic algorithm-multiple linear regressions.  Journal of King Saud University – Science, 32 (2020) 324-331.
[23] S.N. Adawara, G.A. Shallangwa, P. Mamza, A. Ibrahim, Molecular docking and QSAR theoretical model for prediction of phthalazinone derivatives as new class of potent dengue virus inhibitors. Beni-Suef University Journal of Basic and Applied Sciences, 9 (2020).
[24] Y. Isyaku, A. Uzairu, S. Uba, M.T. Ibrahim, A.B. Umar, QSAR, molecular docking, and design of novel 4-(N, N-diarylmethyl amines) Furan-2(5H)-one derivatives as insecticides against Aphis craccivora. Bull. Natl. Res. Cent., 44 (2020) 44.
[25] K. Roy, P. Chakraborty, I. Mitra, P.K. Ojha, S. Kar, R.N. Das, Some case studies on application of “rm2” metrics for judging quality of quantitative structure–activity relationship predictions: Emphasis on scaling of response data. J. Comp. Chem., 34 (2013) 1071-1082.
[26] E.I. Edache, H. Samuel, Y.I. Sulyman, O. Arinze, O.I. Ayine, QSAR and molecular docking analysis of substituted tetraketone and benzyl-benzoate analogs as anti-tyrosine: A novel approach to anti-tyrosine kinase drug design and discovery.Chemistry Research Journal, 5 (2020) 79-100.
[27] S.E. Adeniji, S. Uba, A. Uzairu, Activity modeling of some potent inhibitors against mycobacterium tuberculosis using genetic function approximation approach. Adıyaman University Journal of Science, 9 (2019) 77-98.
[28]  A. Abdullahi, G.A. Shallangwa, M.T. Ibrahim, A.U. Bello, D.E. Arthur, A. Uzairu et al., QSAR studies on some C14-urea tetrandrine compounds as potent anti-cancer agents against Leukemia cell line (K562). JOTCSA, 5 (2019) 1391-402.
[29] A. Tropsha, P. Gramatica, V.K. Gombar, The importance of being earnest: validation is the absolute essential for successful application and interpretation of QSPR models. Mol. Inform., 22 (2003) 69-77.
[30] R. Veerasamy, H. Rajak, A. Jain, S. Sivadasan, C.P. Varghese, R.K. Agrawal, Validation of     QSAR models-strategies and importance.  Int. J. Drug. Des. Discov., 3 (2011) 511–519.
[31] Y. Sun, A.W. Yang,  A.  Hung, G.B. Lenon, Screening for a potential therapeutic agent from the herbal formula in the 4th edition of the Chinese national guidelines for the initial-stage management of COVID-19 via molecular docking.  Evid Based Complement Alternat Med, (2020) 3219840.
[32] J. Lee, X. Cheng, J.M. Swails, M.S. Yeom, P.K. Eastman, J.A. Lemkul et al., CHARMM-GUI input generator for NAMD, GROMACS, AMBER, OpenMM, and    CHARMM/OpenMM simulations using the CHARMM36 additive force field. J. Chem. Theory Comput., 12 (2016) 405-413.
[33]       E.I. Edache, A. Uzairu, P.A. Mamza, G.A. Shallangwa, Computational modeling and analysis of the theoretical structure of thiazolino 2‑pyridone amide inhibitors for Yersinia pseudo-tuberculosis and Chlamydia trachomatis Infectivity. Bull. Sci. Res., 4 (2022) 14-39.
[34] F. Muniba, Molecular dynamics (MD) simulation using Gromacs. Bioinformatics Review, 5 (2019) 12.
[35] Q. Bai, S. Tan, T. Xu, H. Liu, J. Huang, X. Yao, MolAICal: a soft tool for 3D drug design of protein targets by artificial intelligence and classical algorithm. Briefings in Bioinformatics, 00 (2020) 1-12. 
[36] S.E. Adeniji, S. Uba, A. Uzairu, QSAR modeling and molecular docking analysis of some     active compounds against mycobacterium tuberculosis receptor (Mtb CYP121).  Journal of Pathogens, (2018) 24-64.
[37]       F.A. Ugbe, G.A. Shallangwa, A. Uzairu, I. Abdulkadir, A 2-D QSAR modeling, molecular docking study and design of 2-Arylbenzimidazole derivatives as novel leishmania inhibitors: A molecular dynamics study. Adv. J. Chem. A, 6(1) (2023) 50-64.
[38] M.T. Ibrahim, A. Uzairu, S. Uba, G.A. Shallangwa, Design of more potent quinazoline derivative EGFRWT inhibitors for the treatment of NSCLC: a computational approach. Future Journal of Pharmaceutical Sciences, 7 (2021) 140.