Identification of Potential Leishmania N-Myristoyltransferase Inhibitors from Withania somnifera (L.) Dunal: A Molecular Docking and Molecular Dynamics Investigation

Leishmaniasis is a group of infectious diseases caused by Leishmania protozoa. The ineffectiveness, high toxicity, and/or parasite resistance of the currently available antileishmanial drugs has created an urgent need for safe and effective leishmaniasis treatment. Currently, the molecular-docking technique is used to predict the proper conformations of small-molecule ligands and the strength of the contact between a protein and a ligand, and the majority of research for the development of new drugs is centered on this type of prediction. Leishmania N-myristoyltransferase (NMT) has been shown to be a reliable therapeutic target for investigating new anti-leishmanial molecules through this kind of virtual screening. Natural products provide an incredible source of affordable chemical scaffolds that serve in the development of effective drugs. Withania somnifera leaves, roots, and fruits have been shown to contain withanolide and other phytomolecules that are efficient anti-protozoal agents against Malaria, Trypanosoma, and Leishmania spp. Through a review of previously reported compounds from W. somnifera-afforded 35 alkaloid, phenolic, and steroid compounds and 132 withanolides/derivatives, typical of the Withania genus. These compounds were subjected to molecular docking screening and molecular dynamics against L. major NMT. Calycopteretin-3-rutinoside and withanoside IX showed the highest affinity and binding stability to L. major NMT, implying that these compounds could be used as antileishmanial drugs and/or as a scaffold for the design of related parasite NMT inhibitors with markedly enhanced binding affinity.


Introduction
Leishmaniasis is a serious neglected tropical disease that causes several illnesses connected to immune-system dysfunction and poverty, which are sometimes associated with high fatality rates [1]. The disease is manifested in three different ways: Cutaneous leishmaniasis (CL), which is characterized by lesions at the site of infection; mucocutaneous leishmaniasis, which is characterized by invasion and destruction of the mucosa; and visceral leishmaniasis, which is the most severe form because the infection spreads to other organs like the spleen and liver [2]. Each year, around two million new cases are reported worldwide, with cutaneous leishmaniasis accounting for 75% of these infections. The primary cause of cutaneous leishmaniasis is Leishmania major, and out of the 89 countries where it is present, Afghanistan, Brazil, Iran, Peru, Saudi Arabia, and Syria account for 90% Table 1. Standard antileishmanial chemotherapeutic drugs and their suggested mode of antileishmanial effect.

Drug
Mode of Action Sodium stibogluconate

Amphotericin B
Binds the membrane sterols of the parasite and alters its permeability to K + and Mg 2+ selectively

Pentamidine
Interferes with DNA synthesis and modifies the morphology of the kinetoplast

Sodium stibogluconate
Metabolites 2023, 13, x FOR PEER REVIEW 2 of 23 cases are reported worldwide, with cutaneous leishmaniasis accounting for 75% of these infections. The primary cause of cutaneous leishmaniasis is Leishmania major, and out of the 89 countries where it is present, Afghanistan, Brazil, Iran, Peru, Saudi Arabia, and Syria account for 90% of CL cases [1,3]. Despite the long history of the disease, no vaccine is available, and all of the current medications are either ineffective, have unfavorable side effects, or are losing efficacy as resistance emerges [4].
To produce a new anti-Leishmania drug, long isolation, purification, or synthesis processes, as well as excessively expensive in vitro and in vivo biological evaluations, are required [5]. Thus, the present challenge is to identify novel molecules with potential antileishmanial activity while conserving time and money. Utilizing in silico techniques can shorten the time and lower the cost associated with developing new medications [6].
Molecular docking is an in silico technique used to identify correct conformations of small-molecule ligands and estimate the strength of the protein-ligand interaction [7]. Presently, the majority of research for the development of novel drugs focuses on this kind of analysis. Molecular docking has been employed in multiple studies to create new antiparasitic medications while considering a variety of targets to find new Leishmania treatments [8]. Discovering new leishmaniasis treatments is made more enticing by molecules implicated in parasite-specific metabolic processes [9]. The affection of this kind of target implies the death of the parasite and the control of the infection [5]. Worth mentioning is that standard chemotherapeutic antileishmanial drugs operate via related mechanisms-namely, interference in parasite metabolic processes (Table 1) [10]. Targeting specific molecular pathways is a common approach in rational drug design and discovery for developing such leishmaniasis-treating compounds. Pteridine reductase, trypanothione reductase, N-myristoyltransferase (NMT), trypanothione synthetase, inosine-uridine nucleoside hydrolase, and topoisomerases are just a few of the more than 21 potential therapeutic targets of antileishmanial drug discovery that have been reported in the literature and gathered in a book chapter [9].

Drug
Mode of Action Sodium stibogluconate

Amphotericin B
Binds the membrane sterols of the parasite and alters its permeability to K + and Mg 2+ selectively

Pentamidine
Interferes with DNA synthesis and modifies the morphology of the kinetoplast

Amphotericin B
Metabolites 2023, 13, x FOR PEER REVIEW 2 of 23 cases are reported worldwide, with cutaneous leishmaniasis accounting for 75% of these infections. The primary cause of cutaneous leishmaniasis is Leishmania major, and out of the 89 countries where it is present, Afghanistan, Brazil, Iran, Peru, Saudi Arabia, and Syria account for 90% of CL cases [1,3]. Despite the long history of the disease, no vaccine is available, and all of the current medications are either ineffective, have unfavorable side effects, or are losing efficacy as resistance emerges [4].
To produce a new anti-Leishmania drug, long isolation, purification, or synthesis processes, as well as excessively expensive in vitro and in vivo biological evaluations, are required [5]. Thus, the present challenge is to identify novel molecules with potential antileishmanial activity while conserving time and money. Utilizing in silico techniques can shorten the time and lower the cost associated with developing new medications [6].
Molecular docking is an in silico technique used to identify correct conformations of small-molecule ligands and estimate the strength of the protein-ligand interaction [7]. Presently, the majority of research for the development of novel drugs focuses on this kind of analysis. Molecular docking has been employed in multiple studies to create new antiparasitic medications while considering a variety of targets to find new Leishmania treatments [8]. Discovering new leishmaniasis treatments is made more enticing by molecules implicated in parasite-specific metabolic processes [9]. The affection of this kind of target implies the death of the parasite and the control of the infection [5]. Worth mentioning is that standard chemotherapeutic antileishmanial drugs operate via related mechanisms-namely, interference in parasite metabolic processes (Table 1) [10]. Targeting specific molecular pathways is a common approach in rational drug design and discovery for developing such leishmaniasis-treating compounds. Pteridine reductase, trypanothione reductase, N-myristoyltransferase (NMT), trypanothione synthetase, inosine-uridine nucleoside hydrolase, and topoisomerases are just a few of the more than 21 potential therapeutic targets of antileishmanial drug discovery that have been reported in the literature and gathered in a book chapter [9].

Drug
Mode of Action Sodium stibogluconate

Amphotericin B
Binds the membrane sterols of the parasite and alters its permeability to K + and Mg 2+ selectively

Pentamidine
Interferes with DNA synthesis and modifies the morphology of the kinetoplast Binds the membrane sterols of the parasite and alters its permeability to K + and Mg 2+ selectively Pentamidine cases are reported worldwide, with cutaneous leishmaniasis accounting for 75% of these infections. The primary cause of cutaneous leishmaniasis is Leishmania major, and out of the 89 countries where it is present, Afghanistan, Brazil, Iran, Peru, Saudi Arabia, and Syria account for 90% of CL cases [1,3]. Despite the long history of the disease, no vaccine is available, and all of the current medications are either ineffective, have unfavorable side effects, or are losing efficacy as resistance emerges [4].
To produce a new anti-Leishmania drug, long isolation, purification, or synthesis processes, as well as excessively expensive in vitro and in vivo biological evaluations, are required [5]. Thus, the present challenge is to identify novel molecules with potential antileishmanial activity while conserving time and money. Utilizing in silico techniques can shorten the time and lower the cost associated with developing new medications [6].
Molecular docking is an in silico technique used to identify correct conformations of small-molecule ligands and estimate the strength of the protein-ligand interaction [7]. Presently, the majority of research for the development of novel drugs focuses on this kind of analysis. Molecular docking has been employed in multiple studies to create new antiparasitic medications while considering a variety of targets to find new Leishmania treatments [8]. Discovering new leishmaniasis treatments is made more enticing by molecules implicated in parasite-specific metabolic processes [9]. The affection of this kind of target implies the death of the parasite and the control of the infection [5]. Worth mentioning is that standard chemotherapeutic antileishmanial drugs operate via related mechanisms-namely, interference in parasite metabolic processes (Table 1) [10]. Targeting specific molecular pathways is a common approach in rational drug design and discovery for developing such leishmaniasis-treating compounds. Pteridine reductase, trypanothione reductase, N-myristoyltransferase (NMT), trypanothione synthetase, inosine-uridine nucleoside hydrolase, and topoisomerases are just a few of the more than 21 potential therapeutic targets of antileishmanial drug discovery that have been reported in the literature and gathered in a book chapter [9].

Miltefosine
Associated with leishmanial alkyl-lipid metabolism and phospholipid biosynthesis

Paromomycin
Inhibits the protein biosynthesis in sensitive Leishmania parasites N-myristoyltransferase is an example of the molecules implicated in parasite metabolic processes. It is an enzyme that catalyzes the co-translational transfer of a C14 fatty acid from myristoyl-CoA onto the N-terminal glycine residue of a significant subset of proteins [11]. N-myristoylation is believed to play a crucial role in the correct cellular localization and biological function of such modified proteins. Since 1997, it has been established that Leishmania major contains NMT [12]. Later, it was determined that Leishmania sp. needed the NMT enzyme to survive. It has been extensively studied in some organisms, including Trypanosoma brucei, T. cruzi, and Leishmania spp. [13][14][15][16][17]. In a highthroughput screening campaign run by Pfizer, four different scaffolds, aminoacylpyrrolidine (PF-03402623), piperidinylindole (PF-03393842), and thienopyrimidine (PF-003494M), together with a subsequently created fused hybrid compound, 43, were discovered to be highly effective inhibitors of the Labrus donovani NMT (Figure 1) [18]. A different experiment on a collection of 1600 pyrazolyl sulfonamide compounds led to the discovery of compound 2, a potent inhibitor of Leishmania major NMT with negligible activity against Leishmania donovani intracellular amastigotes [5]. Furthermore, several publications and patents indicated the existence of additional NMT inhibitor scaffolds, including pyrrolidines, piperidinylindoles, azetidinopyrimidines, aminomethylindazoles, benzimidazoles, thienopyrimidines, biphenyl derivatives, benzofuranes, benzothiophenes, oxadiazoles, (pyrazolomethyl)-1,3,4-oxadiazoles, and thienopyrimidine [19][20][21][22].  [9].
Inhibits the protein biosynthesis in sensitive Leishmania parasites N-myristoyltransferase is an example of the molecules implicated in parasite metabolic processes. It is an enzyme that catalyzes the co-translational transfer of a C 14 fatty acid from myristoyl-CoA onto the N-terminal glycine residue of a significant subset of proteins [11]. N-myristoylation is believed to play a crucial role in the correct cellular localization and biological function of such modified proteins. Since 1997, it has been established that Leishmania major contains NMT [12]. Later, it was determined that Leishmania sp. needed the NMT enzyme to survive. It has been extensively studied in some organisms, including Trypanosoma brucei, T. cruzi, and Leishmania spp. [13][14][15][16][17]. In a high-throughput screening campaign run by Pfizer, four different scaffolds, aminoacylpyrrolidine (PF-03402623), piperidinylindole (PF-03393842), and thienopyrimidine (PF-003494M), together with a subsequently created fused hybrid compound, 43, were discovered to be highly effective inhibitors of the Labrus donovani NMT (Figure 1) [18]. A different experiment on a collection of 1600 pyrazolyl sulfonamide compounds led to the discovery of compound 2, a potent inhibitor of Leishmania major NMT with negligible activity against Leishmania donovani intracellular amastigotes [5]. Furthermore, several publications and patents indicated the existence of additional NMT inhibitor scaffolds, including pyrrolidines, piperidinylindoles, azetidinopyrimidines, aminomethylindazoles, benzimidazoles, thienopyrimidines, biphenyl derivatives, benzofuranes, benzothiophenes, oxadiazoles, (pyrazolomethyl)-1,3,4oxadiazoles, and thienopyrimidine [19][20][21][22].

Miltefosine
Associated with leishmanial alkyl-lipid metabolism and phospholipid biosynthesis

Paromomycin
Inhibits the protein biosynthesis in sensitive Leishmania parasites N-myristoyltransferase is an example of the molecules implicated in parasite metabolic processes. It is an enzyme that catalyzes the co-translational transfer of a C14 fatty acid from myristoyl-CoA onto the N-terminal glycine residue of a significant subset of proteins [11]. N-myristoylation is believed to play a crucial role in the correct cellular localization and biological function of such modified proteins. Since 1997, it has been established that Leishmania major contains NMT [12]. Later, it was determined that Leishmania sp. needed the NMT enzyme to survive. It has been extensively studied in some organisms, including Trypanosoma brucei, T. cruzi, and Leishmania spp. [13][14][15][16][17]. In a highthroughput screening campaign run by Pfizer, four different scaffolds, aminoacylpyrrolidine (PF-03402623), piperidinylindole (PF-03393842), and thienopyrimidine (PF-003494M), together with a subsequently created fused hybrid compound, 43, were discovered to be highly effective inhibitors of the Labrus donovani NMT ( Figure 1) [18]. A different experiment on a collection of 1600 pyrazolyl sulfonamide compounds led to the discovery of compound 2, a potent inhibitor of Leishmania major NMT with negligible activity against Leishmania donovani intracellular amastigotes [5]. Furthermore, several publications and patents indicated the existence of additional NMT inhibitor scaffolds, including pyrrolidines, piperidinylindoles, azetidinopyrimidines, aminomethylindazoles, benzimidazoles, thienopyrimidines, biphenyl derivatives, benzofuranes, benzothiophenes, oxadiazoles, (pyrazolomethyl)-1,3,4-oxadiazoles, and thienopyrimidine [19][20][21][22].  N-myristoyltransferase inhibitor compounds from medicinal plants may be useful in the quest for alternatives to current treatments for CL. Withania somnifera (L.) Dunal (Solanaceae), also known as Ashwagandha or Indian ginseng, has been used as a traditional herb against a plethora of human medical conditions [23]. It is one of the most extensively used plants in the Unani and Ayurvedic systems of medicine. It has been reported as an important source of withanolides, alkaloids, steroids, flavonoids, nitrogen-containing compounds, and others [23]. Among them, withanolides are a class of highly oxygenated steroids generated from a C28 ergostane skeleton ( Figure 2). They are marker compounds characteristic of Solanaceae plants, particularly those of the genus Withania [23]. They are thought to be responsible for the majority of W. somnifera bioactivity [16]. W. somnifera happens to be one of the prime examples of Rasayana, a branch of Ayurvedic science, a medicinal plant that possesses immunomodulation, anti-cancer, anti-depressant, and neuroprotective properties; promotes the body's resistance to diseases; increases strength and intellect; and delays aging; as well as other biological properties [23]. Closer examination of pertinent in vitro and in vivo studies revealed that W. somnifera extracts molecules' significant bioactivity against several metabolic, reproductive, cardiovascular, neurological, and psychological conditions. Additionally, it has been shown to have antibacterial and antiparasitic properties against Trypanosoma, Leishmania species, and Malaria [24]. The alleged value of W. somnifera in treating leishmaniasis, as proposed by earlier research, is highlighted in the following themes.
Metabolites 2023, 13, x FOR PEER REVIEW 4 of 23 N-myristoyltransferase inhibitor compounds from medicinal plants may be useful in the quest for alternatives to current treatments for CL. Withania somnifera (L.) Dunal (Solanaceae), also known as Ashwagandha or Indian ginseng, has been used as a traditional herb against a plethora of human medical conditions [23]. It is one of the most extensively used plants in the Unani and Ayurvedic systems of medicine. It has been reported as an important source of withanolides, alkaloids, steroids, flavonoids, nitrogencontaining compounds, and others [23]. Among them, withanolides are a class of highly oxygenated steroids generated from a C28 ergostane skeleton ( Figure 2). They are marker compounds characteristic of Solanaceae plants, particularly those of the genus Withania [23]. They are thought to be responsible for the majority of W. somnifera bioactivity [16]. W. somnifera happens to be one of the prime examples of Rasayana, a branch of Ayurvedic science, a medicinal plant that possesses immunomodulation, anti-cancer, antidepressant, and neuroprotective properties; promotes the body's resistance to diseases; increases strength and intellect; and delays aging; as well as other biological properties [23]. Closer examination of pertinent in vitro and in vivo studies revealed that W. somnifera extracts molecules' significant bioactivity against several metabolic, reproductive, cardiovascular, neurological, and psychological conditions. Additionally, it has been shown to have antibacterial and antiparasitic properties against Trypanosoma, Leishmania species, and Malaria [24]. The alleged value of W. somnifera in treating leishmaniasis, as proposed by earlier research, is highlighted in the following themes.

Overview of the Antileishmanial Properties of Withania somnifera
Screening of W. somnifera from different geographical regions has shown that methanolic extract exhibits in vitro antileishmanial action against free-living promastigotes [25,26] and intracellular amastigotes of L. major [26]. Investigation of a solvent-soluble fraction of root and fruit hydromethanolic extracts demonstrated significant antileishmanial promastigote properties of the butanol-soluble fraction from roots and fruits, whereas in vitro growth-inhibitory assessment on axenic amastigotes

Overview of the Antileishmanial Properties of Withania somnifera
Screening of W. somnifera from different geographical regions has shown that methanolic extract exhibits in vitro antileishmanial action against free-living promastigotes [25,26] and intracellular amastigotes of L. major [26]. Investigation of a solvent-soluble fraction of root and fruit hydromethanolic extracts demonstrated significant antileishmanial promastigote properties of the butanol-soluble fraction from roots and fruits, whereas in vitro growth-inhibitory assessment on axenic amastigotes revealed promising activity of the root ethyl acetate-soluble and butanol-soluble fractions [1]. Further bio-guided fractionation revealed that the withanolide-enriched fraction from W. somnifera ethanolic extract was effective [17].
Studies have shown that pure withaferin-A, a prominent withanolide in W. somnifera, possesses antileishmanial activity. It inhibits protein kinase C, which allows the apoptotic topoisomerase I-DNA complex to induce apoptosis [25,27]. Additionally, withanolide Z, a chlorinated withanolide, has been reported to exert an inhibitory effect against L. donovani topoisomerase-1 [28].
The anti-leishmanial activity of the withanolide-enriched extract and pure withanolide of W. somnifera was found to be mediated by the induction of morphological alterations from a spindle to round shape and the loss of flagella and cell integrity in promastigotes. They also induce apoptosis-like cellular death in L. donovani by inducing DNA nicks, apoptosis, and cell-cycle arrest in a dose-and time-dependent manner, actions that were mediated by increasing reactive oxygen species (ROS) production and decreasing mitochondrial potential [29]. Additionally, these withanolide-enriched fractions and pure withanolide of W. somnifera, alone or in combination with other herbal products or standard anti-Leishmanial drugs, were found to modulate hamsters' immunological response to infection with L. donovani.

Immunomodulatory Effects of W. somnifera in Leishmaniasis Infections
One of the immunopathological consequences of active visceral leishmaniasis is the suppression of protective T-helper (Th)-1 cells and the induction of disease-promoting Th-2 cells [30]. Therefore, host immunomodulation is crucial for the treatment of visceral leishmaniasis. In animal-model research, W. somnifera chemotype NMITLI-101 R extract, withaferin A, and the chemotype NMITLI-101 R extract in combination with an ED 50 dose of miltefosine were investigated for their immunoprotective and therapeutic effects against L. donovani infection. It has been observed that the efficacy of W. somnifera was linked to the compelling Th1 immune responses driven by interferon-gamma (IFN-γ) and interleukin-12 (IL-12), as well as dramatically reduced levels of Th2 cytokines (IL-4, IL-10) and transforming growth factor beta (TGF-β). Meanwhile, they significantly increased the levels of NO production, ROS creation, and Leishmania-specific IgG2 antibodies, along with profoundly delayed-type hypersensitivity (DTH) and strong T-cell responses [31][32][33][34].
In another study, when Asparagus racemosus and W. somnifera were used to treat infected mice, the parasite burden was successfully reduced and protective Th1-type immune responses were elicited, resulting in the normalization of biochemical and hematological parameters [30]. In the treatment of visceral leishmaniasis, the combination of W. somnifera extract with cisplatin led to a notable selective elevation of the Th1 type of immunity, verified immunomodulatory action, and a protective impact against the adverse effects of cisplatin on multiple bodily organs. The percentage of CD4 and CD8 T-lymphocytes, as well as the natural killer (NK) cell-associated marker NK1, increased significantly [33].

Molecular Docking Studies in the Development of Antileishmanial Drugs
According to molecular modeling and dynamic investigations of the leishmanial protein kinase C structure, withaferin A and withanone were suggested to disrupt the protein kinase C (PKC) pathway [23,35]. Molecular-docking studies of the binding mode of withaferin-A with pteridine reductase 1 (PTR1), which is involved in pteridine salvage, a crucial enzyme for parasite proliferation, revealed that withaferin-A inhibits PTR-1 enzymes through the uncompetitive mode of inhibition in the parasites.
In summary, in rational antileishmanial drug design and discovery, targeting particular biochemical pathways was found to be a typical strategy for creating leishmaniasis-treating molecules. In order to screen molecules from both natural and synthetic sources, more than 21 molecular targets were used [9]. The structure-guided creation of new lead compounds discovered in high-throughput screening efforts aimed at L. major and L. donovani NMT has led to the identification of effective inhibitors [9,36]. These newly discovered Leishmania NMT inhibitors did not possess the same cellular activity as the enzyme against Leishmania donovani axenic amastigotes, a fact that has been explained by the restricted cellular absorption related to the basic nature of the compounds [19].
W. somnifera's wide variety of the non-basic withanolide derivative and phenolics would be an important source of new hits of Leishmania NMT inhibitors with higher cellular activity against Leishmania parasites. Thus, the goals of this study are to highlight W. somnifera anti-leishmanial properties as well as to discover new L. major NMT inhibitors using in silico molecular-docking and molecular-dynamics analysis of W. somnifera metabolites as potential anti-L. major medicine.

Phytochemical Review and Data Collection
A thorough review of the previously reported compounds from W. somnifera (L.) Dunal in peer-reviewed research papers and international databases, including Science Direct, Pubchem, Google Scholar, SciFinder, etc., afforded a total of 167 different phytoconstituents, which were considered for the present study.

Molecular-Docking Simulation
PyRx software was used for the docking experiments [37]. The RCSB Protein Data Bank (https://www.rcsb.org/ accessed on 20 March 2022) was used to retrieve the threedimensional (3D) structure of L. major N-myristoyltransferase (NMT) in complex with the thienopyrimidine inhibitor IMP-0000096 (PDB ID: 6QDF) ( Figure 3) determined by X-RAY diffraction (Resolution: 1.49 Å) as a target for molecular-docking studies [21]. The 3D chemical structures of the selected molecules were retrieved from PubChem, with polar hydrogen added, partial charge corrected, and energy minimized using the Merck molecular force field (MMFF94x). Molecular-docking analysis was carried out via flexible ligand-fixed receptor-docking parameters using an active complexed ligand active site. The most stable affinity-binding interactions were selected based on the best pose scores. The docking scores and 2D and 3D interactions were recorded [38,39]. BIOVIA Discovery Studio (v21.1.0.20298) was used for 2D and 3D interaction visualization [40]. synthetic sources, more than 21 molecular targets were used [9]. The structure-guided creation of new lead compounds discovered in high-throughput screening efforts aimed at L. major and L. donovani NMT has led to the identification of effective inhibitors [9,36]. These newly discovered Leishmania NMT inhibitors did not possess the same cellular activity as the enzyme against Leishmania donovani axenic amastigotes, a fact that has been explained by the restricted cellular absorption related to the basic nature of the compounds [19].
W. somnifera's wide variety of the non-basic withanolide derivative and phenolics would be an important source of new hits of Leishmania NMT inhibitors with higher cellular activity against Leishmania parasites. Thus, the goals of this study are to highlight W. somnifera anti-leishmanial properties as well as to discover new L. major NMT inhibitors using in silico molecular-docking and molecular-dynamics analysis of W. somnifera metabolites as potential anti-L. major medicine.

Phytochemical Review and Data Collection
A thorough review of the previously reported compounds from W. somnifera (L.) Dunal in peer-reviewed research papers and international databases, including Science Direct, Pubchem, Google Scholar, SciFinder, etc., afforded a total of 167 different phytoconstituents, which were considered for the present study.

Molecular-Docking Simulation
PyRx software was used for the docking experiments [37]. The RCSB Protein Data Bank (https://www.rcsb.org/ accessed on 20 March 2022) was used to retrieve the threedimensional (3D) structure of L. major N-myristoyltransferase (NMT) in complex with the thienopyrimidine inhibitor IMP-0000096 (PDB ID: 6QDF) ( Figure 3) determined by X-RAY diffraction (Resolution: 1.49 Å ) as a target for molecular-docking studies [21]. The 3D chemical structures of the selected molecules were retrieved from PubChem, with polar hydrogen added, partial charge corrected, and energy minimized using the Merck molecular force field (MMFF94x). Molecular-docking analysis was carried out via flexible ligand-fixed receptor-docking parameters using an active complexed ligand active site. The most stable affinity-binding interactions were selected based on the best pose scores. The docking scores and 2D and 3D interactions were recorded [38,39]. BIOVIA Discovery Studio (v21.1.0.20298) was used for 2D and 3D interaction visualization [40].

Molecular-Dynamics Simulation
Molecular-dynamics simulations (MDS) for the best three generated ligand-enzyme complexes were performed using the Nanoscale Molecular Dynamics (NAMD) 3.0 software, applying the CHARMM27 force field, and the MDS was continued for 100 ns following the previously described protocol [41]. The trajectory was stored every 0.1 ns and further analyzed with the VMD 1.9 software. The Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) embedded in the MMPBSA.py module of AMBER18 was utilized to calculate the binding free energy of the docked complex [42]. One hundred frames were processed from the trajectories in total, and the binding free energy was estimated using the following equation: Each of the aforementioned terms requires the calculation of multiple energy components, including van der Waals energy, electrostatic energy, internal energy from molecular mechanics, and a polar contribution to solvation energy.

Results
Earlier phytochemical research on W. somnifera led to the isolation of many phytochemical compounds. According to their chemical structures and the number of isolates, to the best of our knowledge, it could be distinguished as 10 alkaloids, 15 phenolic compounds, 10 sterols, 6 withanones, 6 chloro-containing withanolides, 5 sulfur-containing withanolides, 9 withanamide, 87 withanolides, and 19 withanosides. To discover potential hits that could be used as scaffolds for developing antileishmanial drug candidates, we evaluated these compounds against the Leishmania major NMT utilizing in silico molecular docking and molecular dynamics. The compounds are sorted according to their phytochemical groups parallel to their docking results in descending order according to pose score in Table 2.     Sitoindoside X −12.9 [121] * Compounds with top-scoring docking results.

Discussion
Cutaneous leishmaniasis, caused by more than 20 species of Leishmania parasites, is a derelict tropical disease endemic in most world nations with high incidence rates. In the Kingdom of Saudi Arabia, CL is a significant public-health issue because of several risk factors, including rapid population growth and migration. The disease is endemic in many parts of the kingdom, and L. major and L. tropica are the most commonly detected species [123]. Current chemotherapeutic drugs used for leishmaniasis treatment are discouraging due to associated toxicity, a rase of drug resistance, and high cost. The development of new medications with enough safety and affordability is immediately needed to tackle this disease. W. somnifera is a perennial shrub found in open fields and deserts from the Mediterranean region to Southeast Asia [1]. A reasonable level of safety has been supported by up to 30 clinical studies; there were no obvious side effects or changes in hematological, biochemical, or vital indicators. Pre-clinical studies on chronic toxicity that lasted up to 8 months [16] further supported the safety. As stated above, both in vitro and in vivo tests were used to determine the antileishmanial activity of various organ extracts and sub-extracts of W. somnifera. However, only a small subset of metabolites was tested for antileishmanial activity. Together, these factors make it more important than ever to thoroughly research the phytoconstituents of the plant as possible sustainable and affordable antileishmanial remedies. Due to the economic issue and the hazards of handling a live Leishmania parasite in laboratory work, we adopted computational analysis, molecular docking, and molecular dynamics as a quick and cost-effective scheme to identify potential anti-L. major compounds from W. somnifera.
N-Myristoyltransferase is a prevalent essential enzyme in all Leishmania species [21]. NMT catalyzes the binding of myristate to the amino-terminal glycine residue of a subset of eukaryotic proteins involved in various cellular processes, such as vesicular protein trafficking and signal transduction. NMT has been demonstrated to be essential for viability by classical gene knockout and RNA interference, suggesting that this enzyme has the potential as a target for drug development. It has been proven that NMT is a valid therapeutic target for the treatment of fungus and parasite infections [15,124]. The discovery of the NMT crystal structure establishes a method for structural analysis of inhibitor complexes and structure-assisted drug discovery [21]. Amazingly, it has been demonstrated that the inhibitors only operate on the host enzyme despite the high degree of conservation between the active regions of the parasite and human NMTs [125]. According to the determined structure of β-sheet in the C-terminal domain and protein loops, the residues that are predicted to be the targeted inhibitor site in the NMT are Tyr-80, Val-81, Glu-82, Phe-88, Phe-90, Tyr-92, Asn1-67, Thr-203, Tyr-217, His-219, Phe-232, Tyr-326, Ile-328, Ser-330, Leu-341, Ala-343, Tyr-345, Val-374, Asn-376, Asp-396, Leu-399, Met-420, and Leu421 [126].
As a first step toward utilizing one or more of them as an effective CL medication, the wide range of previously reported compounds from W. somnifera species was chosen for virtual screening as L. major NMT inhibitors. In total, 167 compounds were studied for their binding affinity as NMT enzyme inhibitors using molecular-docking simulation ( Table 2). The results revealed a wide range of compounds' affinity to the NMT enzyme (−9.0 to −24.0 kcal/mol).
Compared with the pose score of the complexed thienopyrimidine inhibitor IMP-0000096 (−14.2 kcal/mol), these compounds could be promising antileishmanial scaffolds.
The introduction of promastigote into the body, the blood circulation of amastigote, phagocytosis and macrophage cell growth, lysis, and blood flow as amastigote once more are among the steps that occur in Leishmania infection [1]. However, NMT inhibitors obtained by synthesis from various research groups were unable to kill parasites inside macrophages at reasonable therapeutic doses, as shown by in vivo assay [2]. On the other hand, withanolide-enriched extracts from W. somnifera reduced the intracellular parasite load by~50% compared to the infected control [34]. Withaferin-A produced a dose-dependent decrease in parasite number inside the macrophages at concentrations of 0.5-1.5 µM [34]. These literature results, together with our findings from docking studies, suggest that the compounds with the top docking scores could be leishmanicidal agents that operate by inhibiting L. major NMT in the free promastigote and the intercellular amastigote forms. Another exciting outcome of our research is the potential to use these molecules as a promising framework for the development of novel leishmanial NMT inhibitors and possibly NMT inhibitors of other protozoa.     Compared with the pose score of the complexed thienopyrimidine inhibitor IMP-0000096 (−14.2 kcal/mol), these compounds could be promising antileishmanial scaffolds.
The introduction of promastigote into the body, the blood circulation of amastigote, phagocytosis and macrophage cell growth, lysis, and blood flow as amastigote once more are among the steps that occur in Leishmania infection [1]. However, NMT inhibitors obtained by synthesis from various research groups were unable to kill parasites inside macrophages at reasonable therapeutic doses, as shown by in vivo assay [2]. On the other hand, withanolide-enriched extracts from W. somnifera reduced the intracellular parasite load by ~50% compared to the infected control [34]. Withaferin-A produced a dosedependent decrease in parasite number inside the macrophages at concentrations of 0.5-1.5 µ M [34]. These literature results, together with our findings from docking studies, suggest that the compounds with the top docking scores could be leishmanicidal agents that operate by inhibiting L. major NMT in the free promastigote and the intercellular

Conclusions
The urgent need to discover a new anti-L. major drug in a reasonable time and at a low cost encouraged us to consider in silico techniques to examine W. somnifera metabolites as an inhibitor of L. major NMT. Among the 167 virtually screened compounds, the phenolic glycoside calycopteretin-3-rutinoside, withanoside IX, and 4,16dihydroxy-5β,6β-epoxyphysagulin D showed promising binding affinity towards Leishmania NMT. In light of the fact that the standard chemotherapeutic antileishmanial medications, such as antimonials, amphotericin B, miltefosine, paromomycin, and pentamidine, work by interfering with parasite metabolic processes, exploring new natural substances that are connected to parasite-specific metabolic processes, such as our discovered NMT inhibitors from W. somnifera, is extremely promising. The results of our investigation will help scientists design in vitro and preclinical animal studies employing Figure 10. Interaction energies of calycopteretin-3-rutinoside, withanoside IX, and 4,16-dihydroxy-5β,6β-epoxyphysagulin D along with the co-crystallized inhibitor inside the active site of NMT throughout 100 ns MDS runs (A-D), respectively). The total interaction energy is the sum of both electrostatic and van der Waals energies.

Conclusions
The urgent need to discover a new anti-L. major drug in a reasonable time and at a low cost encouraged us to consider in silico techniques to examine W. somnifera metabolites as an inhibitor of L. major NMT. Among the 167 virtually screened compounds, the phenolic glycoside calycopteretin-3-rutinoside, withanoside IX, and 4,16-dihydroxy-5β,6βepoxyphysagulin D showed promising binding affinity towards Leishmania NMT. In light of the fact that the standard chemotherapeutic antileishmanial medications, such as antimonials, amphotericin B, miltefosine, paromomycin, and pentamidine, work by interfering with parasite metabolic processes, exploring new natural substances that are connected to parasite-specific metabolic processes, such as our discovered NMT inhibitors from W. somnifera, is extremely promising. The results of our investigation will help scientists design in vitro and preclinical animal studies employing W. somnifera metabolites, with an emphasis on the promising NMT inhibitors we have discovered.