ISC, CAS, Google Scholar     h-index: 20

Document Type : Original Research Article

Authors

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

Abstract

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

Highlights

  • 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.

Keywords

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. https://doi.org/10.1371/journal.pone.0259008.
[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. https://doi.org/10.1186/s42269-022-00874-1.
[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. https://doi.org/10.1007/s40203-022-00123-3.
[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. https://doi.org/10.3390/molecules23112878.
[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. https://doi.org/10.1186/s42269-022-00779-z.
[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). https://doi.org/10.1080/07391102.2020.1852963.
[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. https://doi.org/10.1016/j.cdc.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). https://doi.org/10.1186/s42269-021-00541-x.
[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). https://doi.org/10.1186/s43088-020-00073-9.
[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. https://doi.org/10.1186/s42269-020-00297-w.
[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. https://doi.org/10.1155/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. https://doi.org/10.1186/s43094-021-00279-3.