In Silico Comparison of Bioactive Compounds Characterized from Azadirachta indica with an FDA-Approved Drug against Schistosomal Agents: New Insight into Schistosomiasis Treatment

The burden of human schistosomiasis, a known but neglected tropical disease in Sub-Saharan Africa, has been worrisome in recent years. It is becoming increasingly difficult to tackle schistosomiasis with praziquantel, a drug known to be effective against all Schistosoma species, due to reports of reduced efficacy and resistance. Therefore, this study seeks to investigate the antischistosomal potential of phytochemicals from Azadirachta indica against proteins that have been implicated as druggable targets for the treatment of schistosomiasis using computational techniques. In this study, sixty-three (63) previously isolated and characterized phytochemicals from A. indica were identified from the literature and retrieved from the PubChem database. In silico screening was conducted to assess the inhibitory potential of these phytochemicals against three receptors (Schistosoma mansoni Thioredoxin glutathione reductase, dihydroorotate dehydrogenase, and Arginase) that may serve as therapeutic targets for schistosomiasis treatment. Molecular docking, ADMET prediction, ligand interaction, MMGBSA, and molecular dynamics simulation of the hit compounds were conducted using the Schrodinger molecular drug discovery suite. The results show that Andrographolide possesses a satisfactory pharmacokinetic profile, does not violate the Lipinski rule of five, binds with favourable affinity with the receptors, and interacts with key amino acids at the active site. Importantly, its interaction with dihydroorotate dehydrogenase, an enzyme responsible for the catalysis of the de novo pyrimidine nucleotide biosynthetic pathway rate-limiting step, shows a glide score and MMGBSA of −10.19 and −45.75 Kcal/mol, respectively. In addition, the MD simulation shows its stability at the active site of the receptor. Overall, this study revealed that Andrographolide from Azadirachta indica could serve as a potential lead compound for the development of an anti-schistosomal drug.


Introduction
Human schistosomiasis, commonly known as bilharzia, is a major cause of morbidity in sub-Saharan Africa.This disease is triggered by blood-dwelling parasitic flatworms of the genus Schistosoma.The complex life cycles of Schistosoma species include the infection of both a mammalian definitive host, such as humans, and a freshwater snail as an intermediate host.Restricted availability and access to healthcare facilities, growing poverty, poor sanitary infrastructure, and contact with water activities, especially swimming and irrigated farming, among other factors, predispose people to this debilitating disease [1,2].In terms of the economic impact of parasite infections that warrant considerable concern, schistosomiasis is the second most significant human parasitic disease after malaria.It accounts for over 280,000 deaths each year, while approximately 250 million people are infected with schistosomiasis annually [3,4].
Research from the World Health Organization (WHO) 2013 indicates that more than 40% of cases of schistosomiasis and other terrestrial parasite diseases, excluding malaria, are caused by skin-penetrating cercariae that are present in contaminated water [5,6].
Additionally, studies reveal that schistosomiasis affected over 229 million people in 2015, thereby accounting for a high percentage of those cases in sub-Saharan Africa between the ages of 5 and 14 years [7,8].Ameliorating the schistosomiasis challenge was a target for 2020 for the WHO, with the main target, point of action, and milestones focused on scaling up the work to overcome the global effect posed by schistosomiasis.Moreover, with few successes recorded in eradicating schistosomiasis in most countries such as Japan, Tunisia, and some Caribbean countries, more countries like China, Egypt, and Nigeria are taking steps to eradicate the disease as well [9].Among the vulnerable to schistosomiasis are poor people from rural areas, especially in regions where farming and fishing are the main industries.Household practices like clothes washing in this area also predispose one to contact with contaminated water and risk of infection [10].Children are also more susceptible to schistosomiasis because of poor hygiene and recreational activities like swimming [11,12].
Common species of schistosomes are S. mansoni, S. haematobium, S. japonicum, S. mekongi, and S. intercalatum.Intestinal schistosomiasis is typically more prevalent in sub-Saharan Africa due to the presence of S. mansoni and S. haematobium, among other parasites [7].S. mansoni, a cause of liver and intestine schistosomiasis in South America and Africa, is spread by biomphalaria snails.While biomphalaria snails include many species, such as B. alexandrina, B. sudanica, B. pfeifferi, and B. hoanomphala, the genus Bulinus contains the following species: B. tropicus, B. globosus, B. truncatus, B. forskalli, and B. africanus.S. haematobium is a major cause of urinary schistosomiasis [11,13,14].The microscopic examination of faeces and urine remains the most reliable method for the diagnosis of schistosomal infection.The distinctive size, shape, and presence of lateral or terminal spines make schistosome eggs easily identifiable when observed under a microscope [15,16].
Praziquantel has been one of the treatment choices for schistosomiasis since the mid-1980s.However, compared to praziquantel, which targets adult parasites, the new schistosomiasis medication Artemether targets larvae parasites more efficiently.As a result, Artemether's use for treating schistosomiasis should be limited to prevent the development of drug resistance in the malaria parasite [17][18][19].Major concerns are raised by frequent reports of resistance to oxamniquine, another authorized drug that is used to treat S. mansoni [20,21].The difficulty of dealing with the parasite and the pharmaceutical industry's habitual indifference to tropical diseases prevent the development of a novel strategy for antischistosomal medication research.
Due to its therapeutic and pesticidal properties, the traditional medicinal plant Azadirachta indica possesses the therapeutic potency to treat several human illnesses and act against organisms that are drug-resistant and capable of forming biofilms, such as Schistosomiasis.It also has pharmacological significance and is a highly rich source of bioactive ligands.A. indica has been studied in both experimental and clinical settings with an active chemical isolated from the A. indica tree and has been recognized for some years to have antiviral, antifungal, antibacterial, and insecticidal properties among natural products [22][23][24].The plant has also traditionally been used for parasite infection treatment and gastrointestinal problems [25,26]; however, an effort to drastically control schistosome infections is still based largely on the administration of a chemotherapeutic agent such as praziquantel (PZQ), an oral therapy that is both safe and effective against all human schistosome species.Unfortunately, fast re-infection following treatment with PZQ, as well as certain cases of PZQ resistance in the treatment of specific strains of S. haematobium and S. mansoni as a result of drug pressure, remain a concern.Due to the aforementioned, it has become vital to look for a new antischistosomal drug to replace praziquantel therapy.In this study, we validated the anti-schistosomal activity of compounds identified from A. indica using various in silico tools.

Results and Discussion
Bioinformatics has provided a rapid and versatile platform for drug discovery and the investigation of diseases at the molecular level via the use of computational techniques.This platform has played a key role in the understanding of protein-ligand interaction, which is critical for translational drug development [27][28][29][30][31].There has recently been an upsurge of interest in the exploration of small molecules as an alternative therapy for schistosomiasis treatment, with current therapies showing inadequacy in the treatment and ameliorating the resurgence and spread of this disease [32,33].Therefore, the search for a novel antischistosomal drug to use as an alternative to praziquantel therapy has become necessary.This work, therefore, explored the therapeutic potential of compounds identified from Azadirachta indica as novel anti-schistosomal lead compounds using a computational screening approach [34][35][36].Three clinical targets, namely Schistosoma mansoni Thioredoxin glutathione reductase, dihydroorotate dehydrogenase, and Arginase, were chosen based on their key roles in the maintenance of redox homeostasis, the pyrimidine nucleotide biosynthetic pathway, and metabolism, respectively.The intermolecular interaction of the top three hit compounds (5318517: Andrographolide; 5318767: Nicotiflorin; and 5280804: Isoquercitrin) as well as standard drugs 4891: Praziquantel (a conventional schistosomal medication), 4612: Oxamniquine, and 68827: Artemisinin were studied.
Figure 1 shows the structure of the hit compounds and that of the standard drugs on the market.The pharmacokinetic profile and drug-likeness of a compound have been identified as key factors in identifying potential lead compounds.High/good gastrointestinal absorption is desirable in oral formulations and therefore plays an important role in the delivery of the drug to the target site [37,38].These properties are presented in Tables 1-3 for the hit compounds and the standards.The results of the ADMET analysis (Tables 1-3) show that Andrographolide has a favourable profile, with high gastrointestinal absorption, and conforms with Lipinski's rule (zero violation), like the standard drugs, in comparison with other hit compounds (Nicotiflorin and Isoquercitrin have low gastrointestinal absorption profile), as shown in Table 3.In general, the drug-likeness and ADMET profile of Andrographolide and the standard drug, Pranziquantel, are similar.Moreover, Andrographolide is not a BBB permeant, and does not inhibit CYP enzymes.This makes Andrographolide a promising compound for further analysis.
Figure 1 shows the structure of the hit compounds and that of the standard drugs on the market.The pharmacokinetic profile and drug-likeness of a compound have been identified as key factors in identifying potential lead compounds.High/good gastrointestinal absorption is desirable in oral formulations and therefore plays an important role in the delivery of the drug to the target site [37,38].These properties are presented in Tables 1-3 for the hit compounds and the standards.The results of the ADMET analysis (Tables 1-3) show that Andrographolide has a favourable profile, with high gastrointestinal absorption, and conforms with Lipinski's rule (zero violation), like the standard drugs, in comparison with other hit compounds (Nicotiflorin and Isoquercitrin have low gastrointestinal absorption profile), as shown in Table 3.In general, the drug-likeness and ADMET profile of Andrographolide and the standard drug, Pranziquantel, are similar.Moreover, Andrographolide is not a BBB permeant, and does not inhibit CYP enzymes.This makes Andrographolide a promising compound for further analysis.The docking score could be used to predict the inhibitory potentials of ligands at the receptor site [39][40][41].Thus, to predict the therapeutic efficacy of the hit compounds, molecular docking with the selected proteins was carried out using the Schrodinger's Maestro suite (Glide), and the observed interactions are shown in Table 4.As presented in Table 4, these results suggest a good binding affinity of the hit compounds to the active site of the target proteins.Andrographolide emerged as the best candidate, not just because of its superior glide score, but because of its favourable pharmacokinetic profiles, as discussed and shown in Tables 1-3.The intermolecular interaction of this compound was compared to Pranziquantel.The SmTGR-Andrographolide complex showed a glide score of −4.65 kcal/mol and MMGBSA of −30.49kcal/mol, with the formation of three hydrogen bonds with SER85, ASP84, and GLN86 residues.However, Praziquantel has a glide score of −2.58 kcal/mol, MMGBSA of −28.18 kcal/mol, and one hydrogen bond with SER85 residue.In addition, the SmDHODG-Andrographolide complex had −10.19 kcal/mol docking energy and −45.75 kcal/mol free binding energy.It was also able to produce a similar pose to the co-crystalized inhibitor and interact with similar amino acid residues at the active site, forming hydrogen bonds with PHE61, ARG130, ALA49, and SER53.These, in addition to its hydrophobic interactions, were projected to provide complex stability when compared to Praziquantel.Schistosoma mansoni's arginase-andrographolide complex, on the other hand, had a glide score of −4.34 kcal/mol and hydrogen bonds with ASP211, SER165, and ASN160 when compared with Praziquantel complex with a glide score of −3.41 kcal/mol and a hydrogen bond with SER167, as shown in Table 4.The docking scores observed for Andrographolide across these receptors are comparable with known inhibitors from the CHEMBL database (Tables 5 and 6).The docking validation and ligand interactions for these complexes are shown in Figures 2-12.Following a promising result from the molecular docking, Molecular Dynamic (MD) Simulation analysis was performed for 100 ns to evaluate the atomistic level changes in receptors' binding sites and the stability of the ligands at the binding site of the receptor with respect to the time scale.The root mean square deviation (RMSD), root mean square fluctuation (RMSF), radius of gyration (rGyr), pressure swing adsorption (PSA), and solvent-accessibility surface area (SASA) were analysed, as shown in Table 5 (mean ± SEM) and Figure 13.
The RMSD data give valuable information on variability when applied to very similar proteins.They reveal the conformational variation of a protein-ligand complex over time.The lower the value, the better the stability of the protein-ligand complex.The MD simulation data, as shown in Figure 13, reveal Andrographolide's stable binding conformation with the receptors, with average RMSD values of 1.151, 1.301, and 1.137, respectively, as shown through the 100 ns MD simulation.The RMSF indicates the crucial residues involved in the strongest interactions with a ligand, while the rGyr was analysed to deduce their changes in compactness and to ascertain the protein's stability.Figure 13B reveals the gyration radius spectrum of the Andrographolide complexes with 4Q3T, 6UY4, and 2X8H, with an average rGyr of 4.022, 3.728, and 3.399, respectively (Table 7).It is important to note that the Andrographolide complexes are the most compact of all the ligands, having shown the lowest average rGyr value.The plot showed stability throughout the simulation period for Andrographolide-6UY4 and Andrographolide-2X8H, and therefore it could be argued that it does not cause any distortion in the structure of the proteins.This supports our RMSD and molecular docking findings that Andrographolide may be an effective inhibitor of the Schistosoma mansoni thioredoxin glutathione reductase and dihydroorotate dehydrogenase.

Preparation and Generation of the Protein Grid
The protein preparation and grid generation were carried out using the protein preparation wizard and Glide Grid tools in Maestro's Schrödinger Suite, as previously described [42,43].

Ligand Preparation
Ligand preparation was conducted by OPLS3 optimization at a physiologically relevant PH range of 7.2 ± 0.2.For each ligand structure, potential ionization states were generated, and stereoisomers were computed by modifying specific chiralities while keeping others constant [44].

Computational Approach in Extra Precision (XP) Docking Analysis and Docking Validation
The library of 63 ligands downloaded from PubChem was used for molecular docking studies to discover potential interactions with the active site of the receptors.First, the 63 compounds were screened using glide high throughput virtual screening (HTVS), and the top 10% of the ranking ligands were chosen for (XP) docking as a scoring algorithm to carry out more upscale docking simulation with the receptors.To validate the molecular docking protocol, the co-crystallised ligands were extracted and redocked to the active site of the receptors [45].Also, the FASTA sequence of the target proteins was blasted with the CHEMBL database server (https://www.ebi.ac.uk/chembl/ accessed on 3 April 2024).A set of compounds known to have inhibitory activities on these targets was retrieved and docked to the active site to validate the docking protocol.

Absorption, Distribution, Metabolism, Excretion, and Toxicological Prediction
The Qikprop option of the Maestro tools was used to estimate the (ADME) properties and toxicological potency of ligands together with molecular properties of the lead ligands [46].

Free Binding Energy Determination
The Prime MM-GBSA software (https://www.schrodinger.com/platform/) of the Schrödinger suite was used to calculate the potential free binding energy of the receptorligand complexes to ascertain the balance of their complexes.The OPLS3 and VSGB were chosen as the force field and solvent model, respectively, with other options set at default settings for the free binding energy calculation.This was calculated using the equation ∆Gbind = ∆E + ∆G solv + ∆GSA (1) ∆E = E complex − E protein − E ligand , where E complex , E protein , and E ligand are the minimized energies of the protein-inhibitor complex, protein, and inhibitor, respectively, and ∆G solv = G solv (complex) − G solv (protein) − G solv (ligand), where G solv (complex), G solv (protein), and G solv (ligand) are the solvation free energies of the complex, protein, and inhibitor.∆GSA = GSA (complex) − GSA (protein) − GSA (ligand), where GSA (complex), GSA (protein), and GSA (ligand) represent the surface area energies for the complex, protein, and inhibitor, respectively [47].

Conclusions
Based on the molecular docking analysis, this study revealed that some phytochemicals from Azadirachta indica possess a similar or even higher affinity for SmTGR, SmDHODG, and Schistosoma mansoni Arginase than Praziquantel, the drug available on the market for the treatment of schistosomiasis.The docking analysis shows that Andrographolide, Isoquercitrin, and Nicotiflorin possess a high affinity for these receptors.The MD simulation confirmed this observation by showing that these phytochemicals generated good binding affinity towards the receptors.Andrographolide stands out from all these phytochemicals because of its favourable pharmacokinetic profile and stable binding conformation at the receptor site.Our results therefore suggest that Andrographolide can act as a multi-target inhibitor of Schistosoma mansoni key metabolic enzymes.Therefore, this study revealed that Andrographolide from Azadirachta indica could serve as a potential lead compound for the development of an anti-schistosomal drug, owing to its satisfactory pharmacokinetic profile, the non-violation of the Lipinski rule of five, its good binding affinity, and stability at the receptor site of the Schistosoma mansoni key metabolic enzymes.These make it a viable compound worthy of potential optimization and the experimental evaluation of its biological activity for the possible development of anti-schistosomal drugs.Therefore, in vitro and/or in vivo experimental validation of Andrographolide's anti-schistosomal activity is recommended.

Figure 3 . 23 Figure 4 .
Figure 3.The 3D (top) and 2D (bottom) ligand interactions of the Artemisinin-thioredoxin glutathione reductase complex.The binding site residues' charges are indicated in the 2D image by the colours red for negative, blue for positive, and white for neutral.

2X8H_OxamFigure 4 .Figure 5 .Figure 5 .
Figure 4.The 3D (top) and 2D (bottom) ligand interactions of the Oxamniquine-thioredoxin glutathione reductase complex.The binding site residues' charges are indicated in the 2D image by the colours red for negative, blue for positive, and white for neutral.

Figure 6 .Figure 6 .
Figure 6.The 3D (top) and 2D (bottom) ligand interactions of the Andrographolide-thiore glutathione reductase complex.The binding site residues' charges are indicated in the 2D ima the colours red for negative, blue for positive, and white for neutral.

Figure 7 .Figure 7 .
Figure 7.The 3D (top) and 2D (bottom) ligand interactions of the Nicotiflorin-thiored glutathione reductase complex.The binding site residues' charges are indicated in the 2D imag the colours red for negative, blue for positive, and white for neutral.

Figure 8 .
Figure 8.The 3D (top) and 2D (bottom) ligand interactions of the Oxamniquine-Schistosoma man Arginase complex.The binding site residues' charges are indicated in the 2D image by the colo red for negative, blue for positive, and white for neutral.

Figure 8 .
Figure 8.The 3D (top) and 2D (bottom) ligand interactions of the Oxamniquine-Schistosoma mansoni Arginase complex.The binding site residues' charges are indicated in the 2D image by the colours red for negative, blue for positive, and white for neutral.

Figure 9 .
Figure 9.The 3D (top) and 2D (bottom) ligand interactions of the Andrographolide-Schistoso mansoni Arginase complex.The binding site residues' charges are indicated in the 2D image by t colours red for negative, blue for positive, and white for neutral.

Figure 9 .
Figure 9.The 3D (top) and 2D (bottom) ligand interactions of the Andrographolide-Schistosoma mansoni Arginase complex.The binding site residues' charges are indicated in the 2D image by the colours red for negative, blue for positive, and white for neutral.

Figure 10 .Figure 10 .
Figure 10.The 3D (top) and 2D (bottom) ligand interactions of the Andrographolide-Schistoso mansoni dihydroorotate dehydrogenase complex.The binding site residues' charges are indica in the 2D image by the colours red for negative, blue for positive, and white for neutral.

Figure 11 .Figure 11 .
Figure 11.The 3D (top) and 2D (bottom) ligand interactions of the Artemisinin-Schistosoma man dihydroorotate dehydrogenase complex.The binding site residues' charges are indicated in th image by the colours red for negative, blue for positive, and white for neutral.

Figure 12 .Figure 12 .
Figure 12.The 3D (top) and 2D (bottom) ligand interactions of the Isoquercitrin-Schistosoma ma dihydroorotate dehydrogenase complex.The binding site residues' charges are indicated in th image by the colours red for negative, blue for positive, and white for neutral.Following a promising result from the molecular docking, Molecular Dynamic (M Simulation analysis was performed for 100 ns to evaluate the atomistic level change receptors' binding sites and the stability of the ligands at the binding site of the rece

Table 1 .
Hit compounds from Azadirachta indica with their physicochemical properties.

Table 1 .
Hit compounds from Azadirachta indica with their physicochemical properties.

Table 2 .
Medicinal chemistry profile and drug-likeness of reference drugs and hit compounds from Azadirachta indica.

Table 3 .
Absorptivity, distribution, and metabolism properties of identified ligands compared with reference ligands.

Table 4 .
Molecular interaction of the docking output of the complexes.

Table 5 .
Binding energies of known inhibitors for Schistosoma mansoni thioredoxin glutathione reductase (2X8H) retrieved from the CHEMBL database.

Table 6 .
Binding energies of known inhibitors for Schistosoma mansoni dihydroorotate dehydrogenase (6UY4) retrieved from the CHEMBL database.

Table 7 .
Interactive properties of MDs of the native receptors and protein-ligand interactions.