Indian Journal of Animal Research

  • Chief EditorK.M.L. Pathak

  • Print ISSN 0367-6722

  • Online ISSN 0976-0555

  • NAAS Rating 6.50

  • SJR 0.263

  • Impact Factor 0.4 (2024)

Frequency :
Monthly (January, February, March, April, May, June, July, August, September, October, November and December)
Indexing Services :
Science Citation Index Expanded, BIOSIS Preview, ISI Citation Index, Biological Abstracts, Scopus, AGRICOLA, Google Scholar, CrossRef, CAB Abstracting Journals, Chemical Abstracts, Indian Science Abstracts, EBSCO Indexing Services, Index Copernicus
Indian Journal of Animal Research, volume 58 issue 8 (august 2024) : 1340-1345

In silico Discovery of Potential Natural Inhibitors against Trypanothione Synthetase in Canine leishmaniasis

Pravas Ranjan Sahoo1,*, M. Pashupathi1, Ritun Patra2, Anika1, Komal1, Ajay Kumar1
1Division of Biochemistry, ICAR-Indian Veterinary Research Institute, Bareilly-243 001, Uttar Pradesh, India.
2Department of Veterinary Anatomy, Odisha University of Agriculture and Technology, Bhubaneswar-751 001, Odisha, India.
Cite article:- Sahoo Ranjan Pravas, Pashupathi M., Patra Ritun, Anika, Komal, Kumar Ajay (2024). In silico Discovery of Potential Natural Inhibitors against Trypanothione Synthetase in Canine leishmaniasis . Indian Journal of Animal Research. 58(8): 1340-1345. doi: 10.18805/IJAR.B-4562.
Background: Leishmania infantum, one important intracellular parasite causes most potentially lethal diseases such as leishmaniasis i.e. Visceral leishmaniasis in dogs. Although certain chemical drugs such as pentostam, amphotericin B, miltefosine have been trailed against this disease, but all these drugs induced antibiotic resistance and toxicity in the host. Further, the trypanothione synthetase, a key enzyme of this parasite which catalyzes a reaction, maintaining thiol redox within the cell. 

Methods: The binding study was carried out with selected natural/synthetic phytocompounds/drugs against the modeled trypanothione synthetase though molecular docking. 

Result: The generated protein model with lowest discrete potential energy (DOPE) -19960.97 was found good in quality with z score of -5.19 and quality factor of 61.83% and it was found that the natural inhibitors glycyrrhetic acid (GRA) and Theaflavin (TFN) showed highest binding energies of -7.34 and -6.95 Kcal/mol. This study may be concluded that the natural inhibitors glycyrrhetic acid (GRA) and Theaflavin (TFN) would be potential regimen in treatment of canine leishmaniasis.
Leishmania infantum, one important vector-borne protozoan parasite, belong to the genus Leishmania (Kinetoplastida: Trypanosomatidae) causes deadly disease i.e. visceral leishmaniasis in dogs which causes 1.3 million new cases and 20,000 deaths every year according to WHO 2015 (Singh et al., 2012). This protozoa is transmitted by species of sandfly belonging to the genus Phlebotomus and prevalent throughout the tropical and temperate regions including Africa, China, India, Nepal, Southern Europe, Russia and South America (Johan et al. 2012). The infected dogs showed the clinical signs of fever, low RBC count, skin ulcer and an enlarged liver (Barrett and Croft, 2012). Further, the trypanothione synthase (TS), which catalyses glutathione to spermidine, a key intermediate in maintaining thiol redox within the cell and defending against harmful oxidative effects in such protozoa and can be utilized as drug target against these protozoa (Tetaud et al., 1998). This enzyme has important role in generation of free energy from ATP hydrolysis which conjugates glutathione and spermidine to form the intermediate of glutathionyl spermidine and the final product of trypanothione, helps in the survival of these protozoa (Oza et al., 2006). The lack of crystal structure of Leishmania infantum TS prompts to perform homology modelling of Leishmania infantum TS by utilizing the crystal structure of Leishmania major (Fyf et al., 2008). Although certain drugs such as pentavalent antimonials, Amphotericin B, Miltefsine, Paromomycin, Pentamidine and Sitamaquine have been used against this disease since long time, but all these drugs have limitations of increased drug resistance, cytotoxic side effects, cost and avaibility (Singh et al., 2012). So now this traditional paradigm has changed to search for natural herbal compounds which would be trailed against leishmania infantum. For this, different herbal as well as synthetic compounds were thought to be useful for treating heart diseases and an effect on the permeability of capillaries (Hooper et al., 2012), can be utilized as identifiable natural inhibitor against Leishmania infantum. So this prompts to undertake this current study with objective of identification of most potent herbal compound against trypanothione synthetase of Leishmania infantum through molecular docking experiment.
The study was conducted to discover natural inhibitors against trypanothione synthetase (TS) of Leishmania infantum at the Biochemistry Division, Indian Veterinary Research Institute, Bareilly, U.P. Bhubaneswar from December 2020 to July 2021.
 
Template selection
 
The lack of crystal structures of Leishmania infantum trypanothione synthetase (LiTS) caused to opt homology modeling to determine the structure of LiTS. The sequence of LiTS was retrieved from the NCBI protein database with accession number A0A3S7WXK4 and PSI-BLAST (Altschul et al., 1997) was performed against Protein Data Bank. Trypanothione synthetase of Leishmania major (PDB ID: 2VPM) was considered for template for the reason that the resolution was 2.80 Å and the identity and similarity with Leishmania infantum were found 31% and 67% respectively.
 
Homology modelling and structure validation
 
Homology modeling was performed for the LiTS protein sequence upon Modeller v9.21 (Sali and Blundell, 1993) using Trypanothione synthetase of Leishmania major as template. Ten 3D homology models were generated and the validation for the best model was done using GA341 (Melo et al., 2002) and DOPE (Shen and Sali, 2006) scores. Further the structures of best model were analysed by SAVES validation package (Eisenberg et al., 1997) and RAMPAGE server (http://mordred.bioc.cam.ac.uk/~rapper/rampage.php). ProSA program was also used to analyse the quality of the best model along with the template (Wiederstein and Sippl, 2007). ProQ analysis , ERRAT and Verify 3D analysis was done to check the quality and stability of the best model. Finally, the Model-template superimposition was performed using Pymol v2.1 molecular visualization program (De Lano, 2002).
 
Secondary structure prediction
 
 The secondary structure of LiTS protein was predicted from its complete amino acid sequence (Accession: A0A3S7WXK4) using PSIPRED 4.0 algorithm (Buchan and Jones, 2019).
 
Selection and retrieval of ligand molecules
 
In this study, six natural compounds and four synthetic compounds and with antiprotozoal, anti-inflammatatory and antileishmania activities and one FDA approved antiprotozoal drug were selected for docking studies.
 
Molecular docking
 
In this study, molecular docking between above selected ligands and the drug target Leishmania infantum trypanothione synthetase (LiTS) was performed under AutoDock 4.2 (http://autodock.scripps.edu/) platform using Auto-Dock Tools 4 (Rizvi et al., 2013). For this, prior to docking, the macromolecule and the ligands were prepared with addition of Kollman charges and polar hydrogen atoms using ADT tool. The grid map was assigned around the drug binding cavity of the target protein with size of 126, 126, 126 in x, y, z direction with grid spacing of 0.375Å. In this, docking was performed with lamarckian genetic algorithm (LGA) of 10 runs (25000000 energy evaluation steps for each run), keeping the LiTS protein rigid and  ligand molecules flexible with in drug binding pocket. Binding interactions between atoms and molecules and visualization were carried out using UCSF Chimera v1.14 and PyMol molecular graphics tool (www.pymol.org).
In the current scenario, the canine visceral Leishmaniosis that causes approximately 20,000 deaths, annual incidence of 300,000 new cases and 1 billion of people at risk of infection drives the major concern around the glove for control of this disease (Hailu et al., 2016). Although various chemotherapeutics have been trailed against this disease, but all compounds cause serious side effects, such as renal, pancreatic and hepatic toxicity, teratogenicity, cardiac as well as gastrointestinal problems (Copeland and Aronson, 2015). Further, a single FDI approved drug miltefosine which has been used since 2014 against this disease but it has been facing lot of difficulties in its implementation due to problems of affordability and accessibility (Sunyoto et al., 2018). It has also been reported that the trypanothione synthetase has vital role in the survival and growth of protozoa as this enzyme mediates the final step of tryptophan synthesis. Moreover, it was reported that humans do not have tryptophan synthase, so this enzyme was explored as a potential drug target against leishmania infantum (Chaudhary and Roos, 2005). Although the binding interaction studies, targeting different enzymes of leishmania infantum such as dihydrofolate reductase, thymidylate synthase (Vadloori et al. 2018) Glutamyl cysteine Synthetase (Agnihotri et al., 2017), Inositol phosphoryl ceramide synthase (Norcliffe et al., 2018) with natural inhibitors have been reported in previous finding, but here authors focused on the discovery of potential natural compounds against thetryptophan synthase enzyme. As the crystal structure of leishmania infantum trypanothione synthetase (LiTs) was not avaible, the primary sequence analysis resulted a conserved domain i.e CHAP domain present between 109th to 257th amino acid which mediates the utilization of catalytic cysteine residue in a nucleophilic-attack mechanism during tryptophan synthesis (Bateman et al., 2003). To its support, the avaible crystal structure of leishmania major was retrieved and utilized for generation of protein models.

The model structure of Leishmania infantum trypanothione synthetase (LiTS) was generated using trypanothione synthetase of Leishmania major (PDB ID: 2VPM chain: A) as template. Out of ten models, Model-9 with lowest DOPE and GA341 score of -19960.97 and 0.02942 was considered for best model shown in Table 1 and performed for further structural validation studies. Further the structural superimposition of LdTS model with template 2VPM-A chain both before and after energy minimization (Fig 1) revealed a Root Mean Square Deviation (RMSD) score of 0.112 Å. The RAMPAGE analysis of this model showed (88.5%, 7.1%, 4.3%) amino acid residues in the favored, allowed in the outlier region (Fig 2A). The ProSA-web analysis revealed a Z score of -5.19 (Fig 2B). The ProQ analysis resulted Levitt-Gerstein (LG) and Max sub score of 2.86 and 0.73.Further the ERRAT programme showed the overall quality value for LiTS model of Leishmania infantum was 61.83% (Fig 2C). In addition, Verify3D plot of the modelled protein (Fig 2D) showed PASS and the 3D environment profile resulted 81.93% of the residues have averaged 3D-1D score~*0.2. The structural validation study suggested our model is best in terms of quality and stability.

Table 1: Ten LiTs protein models generated by homology modeling.



Fig 1: The structural super imposition between identified drug target trypanothione synthetase (blue color) of L. infantum and template (PDB ID: 2VPM, Chain A) (red color) of L. major represented.



Fig 2: Ramachandran plot by RAMPAGE presented in (A), z plot which describes the overall quality of model evaluated and deciphered (B), Quality verification plot of the energy minimized model of the LiTs performed using ERRAT shown in (C), Verify 3D plot (D).



Due to lack of crystal structure, the secondary structure of Leishmania infantum trypanothione synthetase (LiTS) protein was predicted from its primary sequence using PSIPRED web server which showed 1 long, 3 medium, 2 short helical regions and 3 medium, 11 short β-sheets within the structure of LiTS protein (Fig 3).

Fig 3: Predicted secondary structural elements for full length trypanothione synthetase of Leishmania infantum. Helix: Pink cylinder; Sheet: Yellow cylinder.



Few drugs such as pentavalent antimonial derivatives, sodium stibogluconate, paromomycin, pentamidine, miltefosine and amphotericin-B showed promising effect against this protozoa, but all these drugs showed potent toxic with serious side effects and also exhibit drug resistance (Oh et al., 2014). These drugs also cause severe adverse reactions such as reversible peripheral neuropathy; pancreatitis, nephrotoxicity cardiotoxicity, pancytopenia, myalgia and bone pain (Croft et al., 2006).

This study has been able to find the binding efficiency of selected phytochemicals (10-hydroxycamptothecin, Theaflavin, Hecogenin acetate, glycyrrhizic acid, convallatoxin, tubocurarine, cafestol, mundulone, pomiferin, catechin) against tryptophan synthase of the protozoa. The separate binding interaction study between selected natural inhibitors and antiprotozoal drug with LiTS protein was done to assay the better therapeutic agent against L. infantum based upon their free binding energy and the inhibition constant of each binding complex which was reported in Table 2.

Table 2: Docking scores of eleven ligands against Leishmania infantum trypanothione synthetase (LiTS) protein are reported.



The result showed that among all selected inhibitor, Glycyrrhetic acid (-7.34 kcal/mol; KI: 4.18uM), Theaflavin (-6.95 kcal/mol; KI: 8.04uM) showed best binding efficiency with LiTS protein. The binding interaction between LiTS with glycyrrhetic acid and Theaflavin are shown in surface and ribbon structure presented in Fig 4 A, B, C and D respectively.

Fig 4: Binding interaction between Leishmania infantum trypanothione synthetase with glycyrrhetic acid showing by surface (A), ribbon (B) and with Theaflavin showing by surface (C), ribbon (D).



In this study it was showed that among natural inhibitor, glycyrrhetic acid and Theaflavin showed highest binding affinity against LiTs which is a good agreement with the findings of Venkatesan et al., 2011 and It has been reported that Glycyrrhetic acid which is a pentacyclic triterpenoid aglycone, a product derived from the plant Glycyrrhiza glabra showed tremendous antiparasitic activity by activation through nitric oxide (NO) upregulation, proinflammatory cytokine expression and NF-κB activation through p38 kinase (Gupta et al., 2015). Similarly theaflavin is a class of natural flavonoids derived from the dried leaves of the plant Camellia sinensis (tea) and related plants with potent antioxidant properties. The antiprozoal activity of theaflavin might be due to inhibition of 1-deoxy-D-xylulose 5-phosphate reductoisomerase, the key enzyme of the MEP terpenoid biosynthetic pathway previously reported by Hui et al., 2016.
This study concludes that the enzyme trypanothione synthase plays a major role in the tryptophan biosynthesis for the survival of L. infantum. Out of all the compounds in this study, Glycyrrhizic acid and Theaflavin which showed higher binding affinity, have been identified likely to be potential inhibitors of LiTs. This docking study hypothesized that Glycyrrhizic acid and Theaflavin could act as potent inhibitors of LiTs enzyme and would provide cost effective natural therapeutics against canine leishmaniasis in nearest future.
The authors are very much thankful to Biochemistry Division, Indian Veterinary Research Institute for providing necessary facilities.
All authors declare that they have no conflict of interest.

  1. Agnihotri, P., Mishra, A.K., Mishra, S., Sirohi, V. K, Sahasrabuddhe, A.A. and Pratap, J.V. (2017). Identification of novel inhibitors of Leishmania donovani γ-Glutamyl cysteine Synthetase using structure based virtual screening, docking, molecular dynamics simulation and in vitro studies. Journal of Chemical Information and Modeling. 57(4): 815-825.

  2. Altschul, S.F., Madden, T.L., Schaffer, A.A., Zhang, J., Zhang, Z. and Miller, W. (1997). Gapped Blast and PSI Blast: A new generation of protein database search programs. Nucleic Acids Research. 25(17): 3389-402.

  3. Barrett, M.P. and Croft, S.L. (2012). Management of trypanosomiasis and leishmaniasis. British Medical Bulletin. 104(1):175-196.

  4. Bateman, A. and Rawlings, N.D. (2003). The CHAP domain: A large family of amidases including GSP amidase and Peptidoglycan hydrolases. Trends in Biochemical Sciences. 28(5): 234-237.

  5. Buchan, D.W.A. and Jones, D.T. (2019). The PSIPRED protein analysis workbench: 20 years on. Nucleic Acids Research. 47(1): 402-407.

  6. Chaudhary, K. and Roos, D.S. (2005). Protozoan genomics for drug discovery. Nature Biotechnology. 23(9): 1089-1091.

  7. Copeland, N.K. and Aronson, N.E. (2015). Leishmaniasis: Treatment updates and clinical practice guidelines review. Current Opinion in Infectious Disease. 28(5): 426-437. 

  8. Croft, S.L., Sundar, S., Fairland, A.H. (2006). Drug resistance in Leishmaniasis. Clinical Microbiology Reviews. 19(1): 111-126.

  9. De Lano, W.L. (2002). Pymol: An open-source molecular graphics tool. CCP4 Newsletter on Protein Crystallography. 40: 82-92.

  10. Eisenberg, D., Luthy, R., Bowie, J.U. (1997). Verify 3D: Assessment of protein models with three-dimensional profiles. Methods in Enzymology. 277: 396-404.

  11. Fyf, P.K., Oza, S.L., Fairlamb, A.H., Hunter, W.N. (2008). Leishmania Trypanothione Synthetase-Amidase structure reveals a basis for regulation of conflicting synthetic and hydrolytic activities. Journal of Biological Chemistry. 285(25): 17672-80.

  12. Gupta, P., Das, P.K., Ukil, A. (2015). Anti-leishmanial effect of 18 β-glycyrrhetinic acid is mediated by toll-like receptor-dependent canonical and non-canonical p38 activation. Antimicrobial Agents and Chemotherapy. 59(5): 2531-9.

  13. Hailu, A., Dagne, D.A., Boelaert, M. (2016). Leishmaniasis in Neglected Tropical Diseases-sub-saharan Africa, [(eds) Gyapong, J. and Oatin, B.] (Berlin: Springer). 87-112. 

  14. Hooper, L., Kay, C., Abdelhamid, A., Kroon, P.A., Cohn, J.S. and Rimm, E.B, et al. (2012). Effect of chocolate, cocoa and flavan-3-ols on cardiovascular health: Asystematic review and meta-analysis of randomized trials. American Journal of Clinical Nutrition. 95(3): 740-751.

  15. Hui, X., Yue, Q., Zhang, D.D., Li, H., Yang, S.Q., Gao, W.Y. (2016). Antimicrobial mechanism of the flavins: They target 1-deoxy-D-xylulose 5-phosphate reductoisomerase, the key enzyme of the MEP terpenoid biosynthetic pathway. Scientific Report. 6: 38945.

  16. Johan, V.G. and Ermias, D. (2012). Visceral Leishmaniasis. Infectious Disease Clinics of North America. 26(2): 309-322.  

  17. Melo, F., Sanchez, R., Sali, A. (2002). Statistical potentials for fold assessment. Protein Science. 11(2): 430-448.

  18. Norcliffe, J.L., Mina, J.G, Alvarez, E., Cantizani, J., Anton, F.D., Colmenarejo, G, et al. (2018). Identifying inhibitors of the Leishmania Inositol Phosphorylceramide Synthase with antiprotozoal activity using a yeast-based assay and ultra- high throughput screening platform. Scientific Report. 8(3938): 1-10.

  19. Oh, S., Kim, S., Kong, S., Yang, G., Lee, N. and Han, D, et al. (2014). Synthesis and biological evaluation of 2-3-dihydroimidazo [1,2-a] benzimidazole derivatives against Leishmania donovani and Trypanosoma cruzi. European Journal of Medicinal Chemistry. 84: 395-403.

  20. Oza, S.L., Wyllie, S., Fairlamb, A.H. (2006). Mapping the functional synthetase domain of Trypanothione synthetase from Leishmania major. Molecular and Biochemical Parasitology. 149: 117-120.

  21. Rizvi, S.M.D., Shakil, S. Haneef, M. (2013). A simple click by click protocol to perform docking: Auto Dock 4.2 made easy for non-bioinformaticians. Excli Journal. 12: 831-857.

  22. Sali, A. and Blundell, T.L. (1993). Comparative protein modeling by satisfaction of spatial restraints. Journal of Molecular Biology. 234(3): 779-815.

  23. Shen, M.Y. and Sali, A. (2006). Statistical potential for assessment and prediction of protein structures. Protein Science. 15(11): 2507-2524.

  24. Singh, N., Kumar, M. and Singh, R.K. (2012). Leishmaniasis/: Current status of available drugs and new potential drug targets. Asian Pacific Journal of Tropical Medicine. 5(6): 485-497.

  25. Sunyoto, T., Potet, J., Boelaert, M. (2018). Why Miltefosine-a life- saving drug for Leishmaniasisis unavailable to people who need it the most. BMJ Global Health. 3: 1-10.

  26. Tetaud, E., Manai, F., Barrett, M.P., Nadeau, K., Walsh, C.T. Fairlamb, A.H. (1998). Cloning and characterization of the two enzymes responsible for trypanothione biosynthesis in Crithidiafa sciculata. Journal Biological Chemistry. 31: 19383-90. 

  27. Vadloori, B., Sharath, A.K., Prabhu, N.P., Maurya, R. (2018). Homology modelling, molecular docking and molecular dynamics simulations reveal the inhibition of Leishmania donovani Dihydrofolate reductase - Thymidylate synthase enzyme by Withaferin-A. BMC Research Notes. 246: 1-7.

  28. Venkatesan, S.K., Saudagar, P., Dubey, V.K. (2011). Identification of novel inhibitor of Trypanothionine synthase from two Leishmania species: Comparative in silico analysis. Journal of Proteins and Proteomics. 2: 41-48.

  29. Wiederstein, M. and Sippl, M.J. (2007). ProSA-web: Interactive web service for the recognition of errors in three dimensional structures of proteins. Nucleic Acids Research. 35(2): W407-W10. 

  30. Yan, L., Sun, Y., Li, X., Liang, Y., Zhu, G., Wang, J, et al. (2015). The effect of hydroxyl camptothecin on wound healing following reduction of the knee intra-articular adhesion in rabbits. Cell Biochemistry and Biophysics. 73(1): 221-7.

Editorial Board

View all (0)