SHORT COMMUNICATION   Open Access    

In silico exploration of Elaeocarpus ganitrus extract phytochemicals on STAT3, to assess their anticancer potential

More Information
  • Elaeocarpus ganitrus Rox of the Elaeocarpaceae family is a broad-leaved medicinal plant and exhaustively used in orthodox systems of treating diseases. However, its anticancer impact and propensity to STAT3 has not yet been analyzed. The plant's extracts were in vitro assayed on the HeLa cell line and subsequently, GC-MS chromatogram of the methanolic, and chloroform extracts of the plant revealed that 106 compounds were present in the extracts. Subsequent filtration using Lipinski rules resulted in 81 phytochemicals being selected for the docking process with pre-selected receptor STAT3 (6NJS). Twenty-six out of 81 phyto-ligands showed high binding energy. Many drugs have weak pharmacokinetic properties and cellular toxicity and consequently, cannot pass through clinical trials. Hence, it is essential to determine the pharmacokinetic parameters of the phytoligands showing preferred binding with receptor 6NJS to consider the apparent bioavailability. The data for pharmacokinetics behavior, bioavailability extent, drug-likeness properties, medicinal chemistry friendliness, and toxicity of 26 phytochemicals with referenced inhibitors was explored. These 26 compounds were further checked for their ADMET properties by using the swissADME and PROTOX-II web server with the known inhibitors plumbagin and sanguinarine to determine the lead phytocompounds. The predictions of ADMET properties obtained six suitable phytocompounds (EG-9, EG-12, EG-13, EG-15, EG-16 and EG-26) of E. ganitrus, and found to be a perfect fit in the bioavailability radar. 2D and 3D interaction of phytoligands with the STAT3 show that the binding is through lys97, suggesting NH2-terminal domain binding of STAT3 with ligands which is the main mono-ubiquitin conjugation spot. Most of the phytoligands interactions exist in the Linker domain and Transactivation domain of the STAT3.
  • Sugars are predominant carbon and energy sources and support plant vegetative and reproductive growth[1]. Transport of sugars across the plant bio-membrane needs the assistance of specific transporters[1,2]. These transporters act as bridges that mediate the distribution of sugars between source–sink organs, which is critical for sugar homeostasis and the cellular exchange of sugar efflux in multicellular organisms[1,3,4].

    SWEETs (Sugars Will Eventually be Exported Transporter) and SUTs (sucrose transporters), MSTs (monosaccharide transporters) are the main known sugar transporters in eukaryotes[5] . Unlike SUTs and MSTs, the relatively newly reported sugar transporter SWEETs are pH-independent transporters. SWEETs play important roles in phloem transport and act as bidirectional transmembrane transporters of sugars along the concentration gradient[3,4]. AtSWEET1 was first identified as a glucose transporter with clear functional characterization[2]. In addition, the SWEET multi-gene family was identified and classified into four clades and the functional divergence of these paralogs were also revealed in Arabidopsis at the same time[2]. For example, Clade II AtSWEET8 contributes to pollen viability and Clade III AtSWEET15 is involved in leaf senescence. Thereafter, another SWEET in Arabidopsis Clade III AtSWEET9 was characterized as an important transporter involved in nectar secretion[6]. Besides the SWEET family in Arabidopsis, the SWEET gene family has been identified in many plants including Tea (Camellia sinensis, Cs), tomato (Solanum lycopersicum, Sl), wheat (Triticum aestivum, Ta) barrel medic (Medicago truncatula, Mt), cabbage (Brassica rapa, Br), daylily (Hemerocallis fulva, Hf), grapevine (Vitis vinifera, Vv), rice (Oryza sativa, Os) and poplar (Populus trichocarpa, Pt and P. alba × P. glandulosa, Pag)[713]. These SWEET homologs belong to the MtN3/saliva family and consist of seven α-helical transmembrane domains (TMs): a tandem repeat of three transmembrane domains (TMs) connected with a linker-inversion TM[2,14].

    SWEETs participate in various biological processes including development, flowering, stress responses and plant-pathogen interaction in plants[4,15]. In addition, different SWEET gene family members show functional divergence or redundancy. In Arabidopsis, AtSWEET8 and 13 support pollen development; AtSWEET11 and 12 provide sucrose to the SUTs for phloem loading and play distinct roles in seed filling; and AtSWEET9 is essential for nectar secretion[24,6,16]. BrSWEET9 in Brassica rapa was also reported to be involved in nectar secretion[17]. Overexpression of PagSWEET7 promotes secondary growth and xylem sugar content[12]. OsSWEET11 and 15 have functions affecting pollen development and are key players in seed filling in rice[18]. SWEET homologs also play important roles in abiotic stress responses. Overexpression of AtSWEET16 promotes freezing tolerance in Arabidopsis[19] while AtSWEET11 and 12 mutants exhibit greater freezing tolerance[20]. AtSWEET15 could be induced by various abiotic stresses including osmotic, drought, salinity, and cold stresses and overexpression of AtSWEET15 results in transgenic plants with hypersensitiveness to cold and salinity stresses[21]. CsSWEET1a and CsSWEET16 were also reported to mediate freezing tolerance[22,23]. Recently, more studies have shown that SWEETs are involved in plant-pathogen interaction and are known as susceptibility (S) genes, acting as targets of effector proteins during host–microbe interactions in many plant species[15]. GhSWEET10, is induced by Avrb6, a transcription activator-like (TAL) effectors from Xanthomonas citri subsp. Malvacearum (Xcm) and is responsible for maintaining virulence of Xcm avrb6 and the cotton susceptibility to infections[24]. OsSWEET11–15 belonging to clade III in rice have been shown to be induced by TAL effector from Xanthomonas oryzae and support pathogen growth[25]. In contrast, some SWEETs could also function as resistance genes. Overexpression of IbSWEET10 can promote resistance to F. oxysporum in sweet potato[26]. Mutation of AtSWEET2 resulted in increased susceptibility to the root necrotrophic pathogen Pythium irregulare[27].

    Mulberry (Morus spp., Moraceae) is a traditional economic crop plant and a new beverage plant. In addition, its fruits are rich in nutrient and bioactive components and the ripening process of mulberry fruits along with sugar accumulation and distribution. Mulberry suffers various abiotic stresses and the disasterous fungal disease sclerotiniose which bursts at the early stage of mulberry fruit development[2830]. Mulberry fruits with sclerotiniose lose their color and flavor and turn pale instead of ripening. C. shiraiana is the dominant causal agent of mulberry sclerotiniose in China, and it results in hypertrophy sorosis sclerotiniose. SWEETs as the important transporters involved in sugar homeostasis are expected to be involved in mulberry fruit development and interaction with sclerotiniose pathogens. However, to date, few studies on SWEETs have been reported in mulberry, although the SWEET gene family may play important roles in mulberry fruit development and responses to abiotic and biotic stresses. Mulberry genome information has been released successively since the Morus notabilis genome was reported in 2013[31]. The chromosome-level genome of M.alba (Ma) was released by Jiao et al. and the genome of M. yunnanensis was recently released by Xia et al.[32,33]. Released genome information makes it possible to perform genome-wide characterization of the SWEET gene family in mulberry. In the present study, a total of 24 SWEET genes were identified in the Morus alba genome and their phylogenetic classification, conserved motifs, gene structures, distribution on chromosomes, cis-elements in promoter regions and tissue expression profile were revealed. In addition, the responses of MaSWEETs to various abiotic stresses and sclerotiniose pathogen infection were also detected. MaSWEET1a was functionally characterized as a negative regulator which increased the mulberry susceptibility to C. carunculoides infection.

    The xylem, phloem, fruits at four different developmental stages (S0, inflorescence; S1, green fruits; S2, reddish fruits; S3, purple fruits)) and diseased fruit infected with C. shiraiana of Morus atropurpurea variety Zhongshen 1 (Mazs) were collected from the National Mulberry Genebank (NMGB) in Zhenjiang, China for expression profiling. Seedlings of the M. alba var. Fengchi and tobacco (Nicotiana Benthamiana) were grown in a chamber at 22 °C with a 16/8 day/night cycle and 40%–60% humidity. C. shiraiana was provided by Professor Zhao and was cultured in potato dextrose agar (PDA) medium.

    Tobacco at the four-leaf stage was used for transient overexpression. M. alba var. Fengchi seedlings at the four-euphylla stage were used for virus-induced gene silencing (VIGS). Four-week-old seedlings with similar growth conditions (~12–15 cm high) were used for treatments under different abiotic stresses. Detailed information for abiotic stress treatments were reported in our previous study[34]. All the above samples were immediately frozen in liquid nitrogen after collection and then stored at −80 °C until use. Three biological replicates were performed for each experiment.

    The M. alba genome sequences (.fasta) and annotation file (.gff) were generously provided by Professor Jiao, who released this genome information. The Hidden Markov Model (HMM) profiles of the SWEET domain (PF03083) were downloaded from the Pfam database (http://pfam.xfam.org/) and used to search the candidate SWEET proteins in the M. alba proteome with HMMER software. In addition, the protein sequences of AtSWEETs, OsSWEETs and PtSWEETs were downloaded from TAIR (www.arabidopsis.org/), TIGR (http://rice.plantbiology.msu.edu/) and phyto-zome (https://phytozome-next.jgi.doe.gov/) respectively, and used as queries to search against the M. alba proteome. The Toolbox for Biologists v1.098774[35] was used to analyze the sequence length, molecular weight and theoretical isoelectric point (pI) values of each MaSWEET protein. The distributions of TM helices were predicted by the TMHMM Server v. 2.0 (www.cbs.dtu.dk/services/TMHMM). Prediction of subcellular localization of MaSWEET proteins using the online Tool WoLF PSORT (www.genscript.com/wolf-psort.html)[36].

    Chromosome location information of MaSWEETs was extracted based on the Morus alba genome annotation file. Tbtools v1.098774 were used to identify syntenic blocks and tandem duplication events using default parameters[37,38]. The results were visualized using Tbtools v1.098774 and both the tandem duplication and block duplication gene pairs were marked.

    MaSWEETs were aligned using clustal W assembled in MEGA11.0. The alignment result was exported and manually speculated for scanning the MtN3 repeats. The online MEME Suite version 5.5.0 was used to identify 7 conserved motifs from 24 amino acid sequences of SWEET genes in Morus alba. The Hidden Markov Model (HMM) profiles of the SWEET domain (PF03083) were downloaded from the Pfam database.

    The gene structure of each MaSWEET was displayed based on the genome sequence and its annotation file using Gene Structure View assembled in Tbtools v1.098774. The upstream 2000 bp sequences were extracted for in silico promoter region analysis. Cis-acting elements were predicted using PlantCARE (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/).

    A neighbor-joining (NJ) phylogenetic tree was constructed using full-length SWEETs protein sequences from A. thaliana, P. trichocarpa, O. sativa and M. alba using MEGA11.0[26] with JTT + G model and bootstrap test with 1000 replicates.

    RNA extraction and cDNA synthesis were performed as in our previous report using Plant RN52 Kit (Aidlab, Beijing, China) and PC54-TRUEscript RT kit (Aidlab, Beijing, China) according to the manual[39]. RT-qPCR (quantitative real-time PCR) was performed to validate the expression patterns of MaSWEETs in different tissues, fruit development stages and stresses using ABI StepOnePlus™ Real-Time PCR System (USA). The primers are available in Supplemental Table S1. Actin was used as a reference gene[40]. Graphpad Prism8.0 was used to visualize the RT-qPCR results and to perform T-test and ANOVA. P value < 0.05 was marked as significant. At least three individuals were used and three technical replicates respectively were performed for RT-qPCR.

    The recombinant plasmids pNC-1304-35S:MaSWEET1a were constructed using nimble cloning[41]. Both recombinant plasmids pNC-1304-35S: SWEET1a and empty vector pNC-Cam1304-35SMCS, as the negative control, were transformed into Agrobacterium tumefaciens GV3101 and then transferred into N. benthamiana leaves via Agrobacterium-mediated transient transformation, as previously reported[41]. Overexpression of MaSWEET1a was determined using RT-qPCR by comparing the expression levels of target genes in transgenic plants with those in the negative controls.

    Virus-induced gene silencing (VIGS) was used to obtain MaSWEET1a/b down-regulated mulberry, in accordance with our previous report[42]. Agrobacterium tumefaciens containing recombinant plasmids pTRV2-MaSWEET1a/b, pTRV1 and pTRV2 (negative control) were cultured in transient transformation buffer and then transferred into mulberry leaves by means of pressure injection. Ten independent mulberry plants were injected. The knock-down efficiency for the target genes was determined by RT-qPCR 15 d after injection by comparing the transgenic plants with the negative controls and wild types.

    Cell death symptoms and the growth condition of C. shiraiana were recorded to estimate the resistance of transgenic plants to C. shiraiana infection[43,44]. C. shiraiana was inoculated at 2 d after infiltration in tobacco and at 10 d after infiltration in mulberry. The cell death symptoms were photographed after inoculation until the sclerotia appeared. The results are representative of at least three biological replicates.

    A total of 24 SWEET homologs were identified based on the genome information of M. alba and named according to their orthologs in A. thaliana, P. trichocarpa or V. vinifera[45]. These MaSWEETs encode proteins ranging from 197 aa to 304 aa with molecular weight from 21.45 to 34.12 kDa and theoretical isoelectric points from 7.16 to 9.61 (Table 1). Subcellular localization prediction of these MaSWEETs showed that most of them (18/24) distributed on membrane structures such as plasma membrane (PM), tonoplast membrane (TM) and chloroplast thylakoid membrane (CTM). Phylogenetic analysis of MaSWEETs and SWEETs from model plants such as A. thaliana, P. trichocarpa and O. sativa showed that four clades were formed by these SWEETs (Fig. 1). According to previous studies, SWEETs in plants were generally classed as four phylogenetic clades which is in agreement with our results[8]. Major MaSWEETs (10/24) together with AtSWEET1, 2 and 3 belong to clade I. Clade II and IV contain five MaSWEETs each and Clade III contains four MaSWEETs (Fig. 1 and Table 1).

    Table 1.  SWEET gene family in Morus alba.
    CladeGene nameAccession no.Gene IDCDS Size
    (bp)
    Protein physicochemical characteristicsTMHsSubcellularMtN3/Saliva (PQ-Loop Repeat)
    Length (aa)MW (kDa)pIAliphatic indexLocalization*Domain Position
    IMaSWEET1aXM_024170697.1-0M.alba_G001204972924226.559.61113.517CTM6-94, 131-209
    IMaSWEET1bXM_024170698.1-0M.alba_G001204967822521.459.51115.746CTM1-49, 86-164
    IMaSWEET2aXM_024163709.1-0M.alba_G001924470823525.948.39129.797CTM18-104, 137-221
    IMaSWEET2bXM_024164777.1-0M.alba_G001086370523425.928.98122.867TM17-101, 137-218
    IMaSWEET2cXM_024163707.1-0M.alba_G001924477425728.588.58132.337PM54-126, 159-243
    IMaSWEET2dXM_024163712.1-0M.alba_G001924468122625.278.22128.896EX23-95, 128-212
    IMaSWEET2eXM_024163708.1-0M.alba_G001924472924226.728.49128.067CTM39-111, 144-228
    IMaSWEET2fXM_024163703.1-0M.alba_G001924477725828.88.49132.568PM41-127, 160-244
    IMaSWEET2gXM_024163711.1-0M.alba_G001924468422725.497.62129.167PM10-96, 129-213
    IMaSWEET3XM_010099554.2-0M.alba_G000306378326029.018.89115.737PM9-98, 132-216
    IIMaSWEET4aXM_010091939.1-0M.alba_G000927673524427.459.28109.397CTM10-98, 134-218
    IIMaSWEET4bXM_010113461.2-0M.alba_G000153673824527.48.98122.087EX11-97, 134-216
    IIMaSWEET5XM_024168739.1-0M.alba_G001829571123626.647.63120.937CY10-93, 131-127
    IIMaSWEET7aXM_010108966.2-0M.alba_G000511077425728.329.57128.567CTM11-95, 134-218
    IIMaSWEET7bXM_010108964.1-0M.alba_G000510978926228.999124.967CTM10-97, 134-218
    IIIMaSWEET10XM_010095631.2-0M.alba_G001801688829533.128.86120.317CTM11-96, 132-216
    IIIMaSWEET11aXM_010114440.2-0M.alba_G001690180426729.639.47124.087CTM11-99, 135-220
    IIIMaSWEET11bXM_010095633.2-0M.alba_G001801591530434.127.57112.897CTM12-99, 134-219
    IIIMaSWEET15XM_010092381.2-0M.alba_G000676788529433.427.16109.017CTM12-99, 133-219
    IVMaSWEET16XM_024167733.1-0M.alba_G001461790930233.29.08114.277CTM20-92, 129-212
    IVMaSWEET17aXM_024171451.1-0M.alba_G001461470823526.468.71119.875CY5-78, 116-198
    IVMaSWEET17bXM_024167902.1-0M.alba_G001461375325028.058.71120.886PM8-93, 131-213
    IVMaSWEET17cXM_024171286.1-0M.alba_G000819572023926.728.94111.727CY6-92, 129-212
    IVMaSWEET17dXM_024171287.1-1M.alba_G0008196723240279.43122.677CTM6-92, 127-213
    * The subcellular localizations were predicted by WoLFPSORT. PM, plasma membrane; EX, extracellular; CY, cytoplasmic; TM, tonoplast membrane; CTM, chloroplast thylakoid membrane.
     | Show Table
    DownLoad: CSV
    Figure 1.  Phylogenetic relationships of the SWEET family genes in Arabidopsis, Oryza sativa, Populus, Vitis vinifera, and Morus alba. The sequences of the 104 SWEET proteins from the above four plant species were aligned by Clustal Omega, and the phylogenetic tree was constructed by the MEGA 11.0 using the NJ method with 1000 bootstrap replicates. The proteins from Arabidopsis, Oryza sativa, Populus, Vitis vinifera, and Morus alba are indicated with the prefixes of At, Os, Pt, Vv, and Ma, respectively.

    MaSWEETs distributed on 12 chromosomes except chromosome 1 and 3. Chromosome 5 occupied five MaSEETs which formed a gene cluster. Chromosome 5 is the chromosome that had the most SWEETs and the following is chromosome 6 and 8 which had four MaWSEETs each. There was only one MaSWEET locating on chromosome 4 (Fig. 2). In addition, there were three MaSWEETs on chromosome 2 and 7 respectively and two MaSWEETs on chromosome 1 and 5 respectively. Gene duplication including block duplication and tandem duplication is the main cause for gene family expansion. Tandem duplications were found on chromosome 2, 6 and 12 (linked by red lines in Fig. 2). It is also interesting to find several possible gene clusters such as MaSWEET1a-b on chromosome 5, MaSWEET2a-g on chromosome 9, MaSWEET7 a-b on chromosome 13, MaSWEET17a-b on chromosome 12 and MaSWEET17c-d on chromosome 6. Two gene pairs (MaSWEET4b/MaSWEET5 and MaSWEET11b/ MaSWEET15) resulting from block duplications were also marked (linked by black lines in Fig. 2).

    Figure 2.  Distribution of MaSWEET genes in Morus alba chromosomes. The tandem gene pairs are linked by red lines. The block duplications gene pairs are marked by black lines. The scale is provided in megabase (Mb).

    MaSWEETs always located on membrane structures with transmembrane domains. The prediction results of MaSWEETs using DeepTMHMM showed that most MaSWEETs posed seven types of transmembrane helices (TMH) named TMH1-7 (Table 1, Supplemental Fig. S1). Alignment and conserved motif analysis showed that almost all MaSWEETs kept the conserved TMH and active sites Tyr and Asp indicating by red full triangles (Fig. 3). The active residues Tyr and Asp were reported to be involved in forming hydrogen bonds to ensure sugar transport activity[14]. In addition, all MaSWEETs except MaSWEET4a had a conserved Ser in each triple helix bundle (THB) which can be phosphorylated and is important for SWEET activity (Fig. 3). MaSWEET4a replaced Ser with Thr at the first Ser phosphorylation site between TMH1 and TMH2 which may also retain similar activity as Thr is also a common phosphorylation site. All MaSWEETs retained the conserved second Ser phosphorylation site between TMH5 and TMH6 (Fig. 3).

    Figure 3.  Multiple sequence alignment of MaSWEET proteins. The positions of the TMHs are underlined. The positions of the active sites of tyrosine (Y) and aspartic acid (D) are indicated by red triangles. The conserved serine (S) phosphorylation sites are indicated by blue triangles.

    MaSWEETs gene structures were identified based on the annotation information of the M. alba genome. In summary, there were six MaSWEETs with six introns, 12 MaSWEETs with five introns and five MaSWEETs with four introns (Fig. 4). In addition, genes clustered together based on phylogenetic analysis are likely to show similar gene structures and length. For example, MaSWEET2a, c, d, e, f and g with six introns, MaSWEET7a and MaSWEET7b with four introns, and MaSWEET17a and MaSWEET17b with four introns.

    Figure 4.  Gene organization of MaSWEETs and cis-elements in promoter regions of MaSWEETs. (a) Phylogenetic tree using 24 MaSWEETs. (b) Exon/intron structures of Morus alba L. SWEETs. (c) Cis-element distribution in the promoter regions of MaSWEETs.

    Promoter region analysis of MaSWEETs indicated the possible function of MaSWEETs in response to hormones and environment stimulus. Among all the 22 types of cis-elements identified, most of them are light response elements accounting for 44% of the total elements (Supplemental Table S2). In addition, hormone response elements were also widely identified in the promoters of MaSWEETs (Fig. 4c). Most MaSWEETs had abscisic acid (ABA), salicylic acid (SA) or methyl jasmonate (MeJA) related response elements in their promoter regions. Especially, MaSWEET1a-b, MaSWEET16, MaSWEET17a-d had cis-elements involved in response to five types of hormones (ABA, SA, MeJA, auxin and gibberellins). Several Myb binding cis-elements were also identified in promoter regions of MaSWEET2a and MaSWEET10 (Fig. 4c, Supplemental Table S2).

    The tissue or organ expression profiles of MaSWEETs were revealed. The MaSWEETs with high sequence identity (> 91%) are hard to distinguished by RT-qPCR and were determined by common primers to reveal their total transcription levels. MaSWEETs in phylogenetic clade I showed quite similar expression patterns with highest expression levels in leaf and relatively higher expression levels in early stages (S0 and S1) of fruit development except MaSWEET1a/b (Fig. 5ad). MaSWEET1a/b showed higher expression levels in fruits during whole fruit development with highest expression level in fruit at S0 stage (Fig. 5a). However, MaSWEETs (MaSWEET10, 11a-b, and 15) in phylogenetic clade III showed preferential expression in fruits especially at the late stages (S2 and S3) (Fig. 5jm). MaSWEETs in phylogenetic clade II had different expression patterns in different tissues or organs. MaSWEET4a and b showed highest expression level in fruit at the S1 stage while MaSWEET5 showed obviously preferential expression in xylem (Fig. 5eg). MaSWEET7a and b showed similar expression pattern with MaSWEET2 cluster and MaSWEET3 from phylogenetic clade I. MaSWEETs in phylogenetic clade IV showed similar expression pattern with highest expression levels in leaf (Fig. 5b, d, h, i). MaSWEET16 also had higher expression in fruit with similar expression level at four different development stages (Fig. 5n).

    Figure 5.  Transcript levels of MaSWEETs in leaves, xylem, phloem, and different development stages of fruit. Three technical replicates were analyzed. Error bars represent SE. Different letters indicate statistically significant differences (Duncan's test, p < 0.05).

    Mulberry sclerotiniose is a fungal disease resulting from fungal pathogen infection. Most (20/24) MaSWEETs showed positive or negative responses to the fungal infection. MaSWEET1a/b, MaSWEET2 cluster, MaSWEET4b, and MaSWEET17 a-d showed a significant decrease of expression levels in diseased fruits with sclerotiniose compared with the expression levels in healthy fruits (Fig. 6). In contrast, MaSWEET2b, MaSWEET3, MaSWEET7b, MaSWEET10 and MaSWEET11a-b showed significant increases of expression levels in diseased fruits. MaSWEETs also played roles in response to various abiotic stresses including drought, water logging, cold and high temperature. MaSWEET1a/b showed a positive response to drought with significant increasing expression levels while other clade I MaSWEETs, MaSWEET2b and 3 significantly decreased their expression level under detected abiotic stresses (Fig. 7ad). In contrast, MaSWEET16 significantly increased its expression level under detected abiotic stresses. MaSWEET4a-b, MaSWEET5 in phylogenetic clade II and MaSWEET11a-b in clade III showed similar response patterns with a significant increase of expression levels in response to low temperature (4 °C), high temperature (40 °C) or drought (Fig. 7eg). MaSWEET15 showed high sensitivity for drought and significant increase of expression levels under drought stress. MaSWEET17a-d showed a negative response to temperature change with a significant decrease of expression levels under low temperature and high temperature treatments (Fig. 7oq).

    Figure 6.  Expression levels of 24 selected MaSWEET genes in response to fungi stress conditions. Three technical replicates were analyzed. Error bars represent SE. Asterisks indicate significant difference as determined by Student’s t-test (* p < 0.05; ** p < 0.01; *** p < 0.001; **** p < 0.0001).
    Figure 7.  Expression levels of 24 selected MaSWEET genes in response to 4 °C, 40 °C , drought, and flood stress conditions. Three technical replicates were analyzed. Error bars represent SE. Asterisks indicate significant difference as determined by Student’s t-test (* p < 0.05 ; ** p < 0.01 ; *** p < 0.001 ; **** p < 0.0001 ).

    MaSWEET1a/b is quite different from other MaSWEETs in clade I based on expression profile analysis, which showed preferential expression in fruits and a negative response to sclerotiniose pathogen (Ciboria shiraiana) infection. Our unpublished data indicated MaSWEET1a as key genes involved in the pathogen infection process based on comparative transcriptome analysis. Transient overexpression of MaSWEET1a in tobacco and VIGS knock-down of MaSWEET1a/b in mulberry were performed. RT-qPCR results validated the successful overexpression of MaSWEET1a in tobacco and knock-down of MaSWEET1a/b in mulberry (Fig. 8b, d). The expression level of MaSWEET1a affected the resistance to C. shiraiana infection in both tobacco and mulberry (Fig. 8a, d). Overexpression of MaSWEET1a decreased the resistance to C. shiraiana infection with more severe cell death symptoms observed in OE-line tobacco (Fig. 8a). Knock-down of MaSWEET1a/b in mulberry could increase the resistance to C. shiraiana infection in mulberry (Fig. 8d). These results proved that MaSWEET1a is an important negative regulator of resistance to C. shiraiana infection in mulberry.

    Figure 8.  Functional characterization of MaSWEET1a/b in tobacco and mulberry. (a) Damage of tobacco overexpressing MaSWEET1a after infection by C. shiraiana. (b) Expression levels of MaSWEET1a in tobacco. MaSWEET1a-T1, T2, T3 are three independent treatments with transient overexpressed MaSWEET1a. (c) Damage of mulberry with knock-down MaSWEET1a/b after infection by C. shiraiana. (d) The expression level of MaSWEET1a/b in mulberry. MaSWEET1a/b-T1, T3, T4 are three independent treatments with down-regulation of MaSWEET1a/b using VIGS.

    Functional studies on SWEETs have revealed that SWEET homologs not only act as loading and unloading transporters of sugars but also play critical roles in various biological processes. Typically, angiosperm genomes contain about 15–25 SWEET genes[4]. In the present study, a total of 24 SWEET genes were identified and clustered into four clades corresponding with the knowledge of SWEETs in angiosperm. Several SWEETs including MaSWEET1a/b, MaSWEET4a-b, MaSWEET10, MaSWEET11a-b and MaSWEET15 showed preferential expression in fruits indicating their possible roles in fruit development. Similar expression preference of AtSWEETs was also reported in Arabidopsis[2]. It is noted that fruit-preferential expressed MaSWEETs still showed temporal expression difference during fruit development indicating time-course regulation of MaSWEETs for fruit ripening in mulberry. Early-stage expressed MaSWEET1a/b and MaSWEET4a/b further differed from late-stage expressed MaSWEET10, MaSWEET11a/b and MaSWEET15 in terms of detailed expression patterns during fruit ripening. MaSWEET2b and 2 cluster genes (MaSWEET2a, c-g) showed preferential expression in leaves which is similar with the ortholog ZjSWEET2.2 in Ziziphus jujuba. ZjSWEET2.2 was reported to be involved in mediating sugar loading in leaves[46]. MaSWEET3, 7a-b,16 and 17a-d also showed highest expression levels in leaves indicating their possible roles in sugar source loading or unloading. MaSWEET11b showed higher expression levels in phloem and its ortholog AtSWEET11 in Arabidopsis was reported to be involved in sugar phloem loading[3].

    Sugar signal is critical for plants in response to various stresses. Previous studies have shown that SWEETs participated in abiotic and biotic responses in many plant species including arabidopsis and rice[21, 25]. AtSWEET11, 12, 15 and 16 were reported to be involved in affecting cold tolerance in Arabidopsis[1921]. HfSWEET17 was also reported as a positive regulator of resistance to cold stress in daylily[11]. Cold environment (4 °C) induced expression of MaSWEET4a, 4b, 5, 11a, 11b and 16 in mulberry. Interestingly, these cold-induced MaSWEETs can also be induced by high temperature. MaSWEET15 which is the ortholog of AtSWEET15 can be induced by drought as well as low or high temperature. MaSWEET15, MaSWEET1a/b. 4a, 4b, 5, 7a, 11a, 11b, and 16 also showed positive responses to drought. It is obvious that some SWEET genes can be induced by different stresses. AtSWEET15 was also reported to be induced by osmotic, drought and salinity[21]. MaSWEET4a, 4b, 11a, 11b and 16 can be induced by low or high temperature and drought indicating their important roles in response to various abiotic stresses in mulberry.

    SWEETs were generally thought to 'support the enemy' during infection. SWEETs especially those that function as exporters generally facilitate the export of sugars out of host cells, which support pathogen growth in the apoplasm[15,47,48]. Clade III SWEETs including AtSWEET11, 12, OsSWEET11 were characterized as negative regulators of resistance to fungal infection and Clade III SWEETs including OsSWEET11, 13, 14 and GhSWEET10 were characterized as negative regulators of resistance to bacterial pathogen infection[2,24,49,50]. In contrast, clade I AtSWEET2, a glucose importer and clade III IbSWEET10 were reported as positive regulators of resistance to fungal infection[48,51]. Therefore, roles of SWEETs in response to pathogen infection may be quite different. Most MaSWEETs were disturbed in diseased fruits that resulted from sclerotiniose pathogen infection. MaSWEET2b, MaSWEET3 in clade I, MaSWEET7b in clade II, MaSWEET10 and MaSWEET11a-b in clade III showed significant increase in expression levels in diseased fruits while MaSWEET1a/b, MaSWEET2 cluster (MaSWEET2a, c-g) in clade I showed significant decrease of expression levels in diseased fruits. MaSWEET1a was further validated as a negative regulator of resistance to C. shiraiana infection. Given the fact that MaSWEET1a/b was repressed in diseased fruits, it is likely that a possible pathway through repression of MaSWEET1a exists in mulberry to defense pathogen infection.

    In conclusion, we have performed a genome-wide investigation of SWEET genes in Morus and a comprehensive analysis including phylogenetic analysis, promoter analysis and expression profile analysis was also carried out. Their possible roles in development and response to abiotic and biotic stresses were addressed. In particular, the functional role of MaSWEET1a in regulation of tolerance to C. shiraiana infection was validated using both VIGS knock-down and transient overexpression in tobacco combined with inoculation of C. shiraiana. The results in this study provides a foundation for studying the function of the SWEET family in mulberry plants and provides a negative regulator of resistance to C. shiraiana infection for further genetic modification.

    This work was jointly supported by the National Natural Science Foundation of China (32201526), Crop Germplasm Resources Protection Project of the Ministry of Agriculture and Rural Affairs of the People's Republic of China (19200382), National Infrastructure for Crop Germplasm Resources (NCGRC-2020-041), and China Agriculture Research System of MOF and MARA (CARS-18). We thank Professor Aichun Zhao at Southwest University for providing us C. shiraiana strains and Professor Feng Jiao at Northwest University of Agriculture and Forestry who provided us the genome annotation file of M. alba.

  • The authors declare that they have no conflict of interest.

  • [1]

    Yu H, Lee H, Herrmann A, Buettner R, Jove R. 2014. Revisiting STAT3 signalling in cancer: new and unexpected biological functions. Nature Reviews Cancer 14:736−46

    doi: 10.1038/nrc3818

    CrossRef   Google Scholar

    [2]

    Zou S, Tong Q, Liu B, Huang W, Tian Y, et al. 2020. Targeting STAT3 in cancer immunotherapy. Molecular Cancer 19:145

    doi: 10.1186/s12943-020-01258-7

    CrossRef   Google Scholar

    [3]

    Fan Y, Mao R, Yang J. 2013. NF-κB and STAT3 signaling pathways collaboratively link inflammation to cancer. Protein & Cell 4:176−85

    doi: 10.1007/s13238-013-2084-3

    CrossRef   Google Scholar

    [4]

    Yang X, Xu L, Yang L, Xu S. 2023. Research progress of STAT3-based dual inhibitors for cancer therapy. Bioorganic & Medicinal Chemistry 91:117382

    doi: 10.1016/j.bmc.2023.117382

    CrossRef   Google Scholar

    [5]

    Dong J, Cheng XD, Zhang WD, Qin JJ. 2021. Recent update on development of small-molecule STAT3 inhibitors for cancer therapy: from phosphorylation inhibition to protein degradation. Journal of Medicinal Chemistry 64:8884−915

    doi: 10.1021/acs.jmedchem.1c00629

    CrossRef   Google Scholar

    [6]

    Siveen KS, Sikka S, Surana R, Dai X, Zhang J, et al. 2014. Targeting the STAT3 signaling pathway in cancer: role of synthetic and natural inhibitors. Biochimica et Biophysica Acta 1845:136−54

    doi: 10.1016/j.bbcan.2013.12.005

    CrossRef   Google Scholar

    [7]

    Ijaz S, Akhtar N, Khan MS, Hameed A, Irfan M, et al. 2018. Plant derived anticancer agents: a green approach towards skin cancers. Biomedicine & Pharmacotherapy 103:1643−51

    doi: 10.1016/j.biopha.2018.04.113

    CrossRef   Google Scholar

    [8]

    Salehi B, Machin L, Monzote L, Sharifi-Rad J, Ezzat SM, et al. 2020. Therapeutic potential of quercetin: new insights and perspectives for human health. ACS Omega 5:11849−72

    doi: 10.1021/acsomega.0c01818

    CrossRef   Google Scholar

    [9]

    Majolo F, de Oliveira Becker Delwing LK, Marmitt DJ, Bustamante-Filho IC, Goettert MI. 2019. Medicinal plants and bioactive natural compounds for cancer treatment: important advances for drug discovery. Phytochemistry Letters 31:196−207

    doi: 10.1016/j.phytol.2019.04.003

    CrossRef   Google Scholar

    [10]

    Patel B, Das S, Prakash R, Yasir M. 2010. Natural bioactive compound with anticancer potential. International Journal of Advances in Pharmaceutical Sciences 1:32−41

    doi: 10.5138/ijaps.2010.0976.1055.01003

    CrossRef   Google Scholar

    [11]

    Mousavi SM, Hashemi SA, Behbudi G, Mazraedoost S, Omidifar N, et al. 2021. A review on health benefits of Malva sylvestris L. nutritional compounds for metabolites, antioxidants, and anti-inflammatory, anticancer, and antimicrobial applications. Evidence-Based Complementary and Alternative Medicine 2021:5548404

    doi: 10.1155/2021/5548404

    CrossRef   Google Scholar

    [12]

    Krishna P, Kumari NR, Manisree V, Rani KS, Deepthi BVP, Sharma JVC. 2019. Medicinal benefits of Elaeocarpus Ganitrus (Rudraksha) - A divine herb. World Journal of Pharmaceutical Research 8:552−65

    Google Scholar

    [13]

    Mahajanakatti AB, Deepak TS, Achar RR, Pradeep S, Prasad SK, et al. 2022. Nanoconjugate synthesis of Elaeocarpus ganitrus and the assessment of its antimicrobial and antiproliferative properties. Molecules 27:2442

    doi: 10.3390/molecules27082442

    CrossRef   Google Scholar

    [14]

    Das PK. 2015. Phytochemical screening of methanolic extracts of different parts of rudraksh plant (Elaeocarpus ganitrus). Journal of Biological Sciences 15:111−12

    doi: 10.3844/OJBSCI.2015.111.112

    CrossRef   Google Scholar

    [15]

    Motallebi M, Bhia M, Rajani HF, Bhia I, Tabarraei H, et al. 2022. Naringenin: a potential flavonoid phytochemical for cancer therapy. Life Sciences 305:120752

    doi: 10.1016/j.lfs.2022.120752

    CrossRef   Google Scholar

    [16]

    Zhang Y, Liu X, Ruan J, Zhuang X, Zhang X, et al. 2020. Phytochemicals of garlic: promising candidates for cancer therapy. Biomedicine & Pharmacotherapy 123:109730

    doi: 10.1016/j.biopha.2019.109730

    CrossRef   Google Scholar

    [17]

    Fulda S, Debatin KM. 2006. Resveratrol modulation of signal transduction in apoptosis and cell survival: a mini-review. Cancer Detection and Prevention 30:217−23

    doi: 10.1016/j.cdp.2006.03.007

    CrossRef   Google Scholar

    [18]

    Ahmad K, Bhat AR, Athar F. 2017. Pharmacokinetic evaluation of Callistemon viminalis derived natural compounds as targeted inhibitors against δ-opioid receptor and farnesyl transferase. Letters in Drug Design & Discovery 14:488−99

    doi: 10.2174/1570180814666161214114322

    CrossRef   Google Scholar

    [19]

    Sudradjat SE, Timotius KH. 2022. Pharmacological properties and phytochemical components of Elaeocarpus: a comparative study. Phytomedicine Plus 2:100365

    doi: 10.1016/j.phyplu.2022.100365

    CrossRef   Google Scholar

    [20]

    Kumar TS, Shanmugam S, Palvannan T, Bharathi Kumar VM. 2008. Evaluation of antioxidant properties of Elaeocarpus ganitrus roxb. leaves. Iranian Journal of Pharmaceutical Research 7(3):211−15

    Google Scholar

    [21]

    Sultana B, Anwar F, Ashraf M. 2009. Effect of extraction solvent/technique on the antioxidant activity of selected medicinal plant extracts. Molecules 14:2167−80

    doi: 10.3390/molecules14062167

    CrossRef   Google Scholar

    [22]

    Ozigis HO, Olaifa KA, Agbeja AO, Asabia LO, Akindolu DR, et al. 2023. Qualitative phytochemical analysis of leave and stem bark of Zanthoxylum zanthoxyloides and Zanthoxylum gilletti. Journal of Chemical Society of Nigeria 48(3):891

    doi: 10.46602/jcsn.v48i3.891

    CrossRef   Google Scholar

    [23]

    Dhivya R, Jaividhya P, Riyasdeen A, Palaniandavar M, Mathan G, et al. 2015. In vitro antiproliferative and apoptosis-inducing properties of a mononuclear copper(II) complex with dppz ligand, in two genotypically different breast cancer cell lines. BioMetals 28:929−43

    doi: 10.1007/s10534-015-9877-1

    CrossRef   Google Scholar

    [24]

    Daina A, Michielin O, Zoete V. 2017. SwissADME: a free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules. Scientific Reports 7:42717

    doi: 10.1038/srep42717

    CrossRef   Google Scholar

    [25]

    Lipinski CA. 2004. Lead- and drug-like compounds: the rule-of-five revolution. Drug Discovery Today Technologies 1:337−41

    doi: 10.1016/j.DDTec.2004.11.007

    CrossRef   Google Scholar

    [26]

    Khan A, Mohammad T, Shamsi A, Hussain A, Alajmi MF, et al. 2022. Identification of plant-based hexokinase 2 inhibitors: combined molecular docking and dynamics simulation studies. Journal of Biomolecular Structure & Dynamics 40:10319−31

    doi: 10.1080/07391102.2021.1942217

    CrossRef   Google Scholar

    [27]

    Banerjee P, Eckert AO, Schrey AK, Preissner R. 2018. ProTox-II: a webserver for the prediction of toxicity of chemicals. Nucleic Acids Research 46:W257−W263

    doi: 10.1093/nar/gky318

    CrossRef   Google Scholar

    [28]

    Beg A, Khan FI, Lobb KA, Islam A, Ahmad F, et al. 2019. High throughput screening, docking, and molecular dynamics studies to identify potential inhibitors of human calcium/calmodulin-dependent protein kinase IV. Journal of Biomolecular Structure & Dynamics 37:2179−92

    doi: 10.1080/07391102.2018.1479310

    CrossRef   Google Scholar

    [29]

    Tsaioun Katya, Kates SA. (Eds) 2011. ADMET for medicinal chemists: a practical guide. Hoboken, New Jersey (simultaneously in Canada): John Wiley & Sons. https://doi.org/10.1002/9780470915110

    [30]

    Onawole AT, Sulaiman KO, Adegoke RO, Kolapo TU. 2017. Identification of potential inhibitors against the Zika virus using consensus scoring. Journal of Molecular Graphics & Modelling 73:54−61

    doi: 10.1016/j.jmgm.2017.01.018

    CrossRef   Google Scholar

    [31]

    Siramshetty VB, Nickel J, Omieczynski C, Gohlke BO, Drwal MN, et al. 2016. WITHDRAWN—a resource for withdrawn and discontinued drugs. Nucleic Acids Research 44:D1080−D1086

    doi: 10.1093/nar/gkv1192

    CrossRef   Google Scholar

    [32]

    Ray S, Zhao Y, Jamaluddin M, Edeh CB, Lee C, et al. 2014. Inducible STAT3 NH2 terminal mono-ubiquitination promotes BRD4 complex formation to regulate apoptosis. Cellular Signalling 26:1445−55

    doi: 10.1016/j.cellsig.2014.03.007

    CrossRef   Google Scholar

  • Cite this article

    Mehnaj, Bhat AR, Athar F. 2024. In silico exploration of Elaeocarpus ganitrus extract phytochemicals on STAT3, to assess their anticancer potential. Medicinal Plant Biology 3: e009 doi: 10.48130/mpb-0024-0010
    Mehnaj, Bhat AR, Athar F. 2024. In silico exploration of Elaeocarpus ganitrus extract phytochemicals on STAT3, to assess their anticancer potential. Medicinal Plant Biology 3: e009 doi: 10.48130/mpb-0024-0010

Figures(4)  /  Tables(11)

Article Metrics

Article views(2472) PDF downloads(461)

Other Articles By Authors

SHORT COMMUNICATION   Open Access    

In silico exploration of Elaeocarpus ganitrus extract phytochemicals on STAT3, to assess their anticancer potential

Medicinal Plant Biology  3 Article number: e009  (2024)  |  Cite this article

Abstract: Elaeocarpus ganitrus Rox of the Elaeocarpaceae family is a broad-leaved medicinal plant and exhaustively used in orthodox systems of treating diseases. However, its anticancer impact and propensity to STAT3 has not yet been analyzed. The plant's extracts were in vitro assayed on the HeLa cell line and subsequently, GC-MS chromatogram of the methanolic, and chloroform extracts of the plant revealed that 106 compounds were present in the extracts. Subsequent filtration using Lipinski rules resulted in 81 phytochemicals being selected for the docking process with pre-selected receptor STAT3 (6NJS). Twenty-six out of 81 phyto-ligands showed high binding energy. Many drugs have weak pharmacokinetic properties and cellular toxicity and consequently, cannot pass through clinical trials. Hence, it is essential to determine the pharmacokinetic parameters of the phytoligands showing preferred binding with receptor 6NJS to consider the apparent bioavailability. The data for pharmacokinetics behavior, bioavailability extent, drug-likeness properties, medicinal chemistry friendliness, and toxicity of 26 phytochemicals with referenced inhibitors was explored. These 26 compounds were further checked for their ADMET properties by using the swissADME and PROTOX-II web server with the known inhibitors plumbagin and sanguinarine to determine the lead phytocompounds. The predictions of ADMET properties obtained six suitable phytocompounds (EG-9, EG-12, EG-13, EG-15, EG-16 and EG-26) of E. ganitrus, and found to be a perfect fit in the bioavailability radar. 2D and 3D interaction of phytoligands with the STAT3 show that the binding is through lys97, suggesting NH2-terminal domain binding of STAT3 with ligands which is the main mono-ubiquitin conjugation spot. Most of the phytoligands interactions exist in the Linker domain and Transactivation domain of the STAT3.

    • Cancer is a cluster of diseases and STAT3 protein has significant roles in all types of cancer. Signal transducer and activator of transcription (STAT), belong to the family of cytoplasmic transcription factors, activate and transduce extracellular growth factor, and also affect cytokine signals and affect gene transcriptional events. STAT3 mutant intrinsically alone is enough to instigate oncogenic transformation, and tumorigenesis[13]. A survey of the current literature reveals that STATs have transactivated domains and play a significant role in cancer migration and invasion. Hampering of c-Src kinase activity or downregulation of STAT3 signaling stimulates apoptosis[4]. The study of chemical interactions between STAT3 receptor and phytochemicals assist in drug designing and hence in cancer therapy[5]. There are a variety of phytochemicals that have a high propensity to modulate directly or indirectly the STAT3 signaling pathway. Triterpenoids like betulinic acid, polyphenols curcuminoids, plumbagin a naphthoquinone, diosgenin a steroid, hydroxycinnamic acid, and thymoquinone are the phytochemicals that suppress STAT3 expression[6]. Many plant-derived phytochemicals manifest high anticancer activity and lead researchers to adopt integrated multifaceted research techniques[711]. Though, Elaeocarpus ganitrus Roxb. (also known as Rudraksha) constitutively placed in Ayurvedic system, also has anticancer potential[12]. Recently its silver nanoparticle has been assessed for anticancer and antiproliferative activities[13]. Its impact on STAT is yet to be explored. In the last few decades, phytochemical composition of Elaeocarpus genus has been extensively investigated. Phytochemicals of various extracts of different parts of the plant showed the presence of alkaloids, flavonoids, carbohydrates, glycosides, proteins, quinine, coumarins, tannins, minerals, vitamins, saponins, phenolic compounds, and fixed oils in a high concentration, thus adding to its medicinal value[14]. The pharmacological screening of metabolites like polyphenols, alkaloids, terpenoids and flavonoids have been explored to demonstrate cancer pathways to ascertain possible mechanism[1518]. As stated in the literature, the beads and the bark of the plants have been extensively studied while the leaves of the E. ganitrus have not been studied for their anticancer efficiency. Besides, leaves of the plants were shown to have good antioxidant potential[19,20]. The emphasis of the study is to identify phytochemicals retrieved from Elaeocarpus ganitrus leaf research data and GC-MS profiling. To accentuate, the binding role of Elaeocarpus ganitrus phytochemicals with STAT3 receptor, their ADME properties and pharmacokinetic studies were investigated.

    • The chemicals and solvents used in the extraction and phytochemical analysis were of analytical grade, sourced from Sigma-Aldrich. MTT (3-[4,5-dimethylthiazol-2-1yl]-2,5 diphenyl tetrazolium bromide) was also procured from Sigma-Aldrich. HeLa cells were obtained from the National Centre for Cell Sciences (NCCS), Pune, India. Fetal bovine serum and Dulbecco's Modified Eagle's Media were acquired from Gibco-life technologies.

    • Fresh leaves of E. ganitrus were purchased from Patanjali Herbal Garden Nursery in Panchayanpur, Uttarakhand, India. Authentication of the Elaeocarpus ganitrus was conducted by the Department of Botany, Jamia Hamdard, New Delhi, India, and the voucher specimen was deposited at the University.

    • Leaves of E. ganitrus were carefully washed, air-dried for ten days, and ground to a fine powder. A sample of 1,000 grams of powder was exhaustively extracted three times with 100% methanol (10 times weight/volume) at room temperature for 72 h using a soxhlet apparatus. The resulting crude methanol extract was fractionated successively with solvents in increasing polarity order: heptane, chloroform, ethyl acetate, methanol, and water. The residue was air-dried and utilized for the subsequent solvents. The fractions obtained from each solvent were filtered, dried under vacuum using a rotary evaporator, and stored at 40 °C until use[21].

    • The presence or absence of phytochemicals such as terpenoids, steroids, saponins, flavonoids, glycosides, tannins, and phenols in the chloroform and methanol extracts of E. ganitrus leaves was determined following the standard methodology[22].

    • The HeLa cell line was stored in Dulbecco's Modified Eagle's Medium which is rich in 10% Fetal Bovine Serum, 1% antibiotic solution, 25 mM sodium bicarbonate, and 10 mM HEPES in a 5% CO2 humidified atmosphere at 37 °C in an air jacketed incubator. The stock culture was perpetuated in the exponentially growing phase by passaging as, monolayer culture with 0.02% EDTA. Dislodged cells suspended in complete medium were routinely reseeded.

    • The cytotoxic effects of the various fractions of E. ganitrus leaf on the HeLa cancer cell line were evaluated using the MTT assay. Cells were seeded overnight, and exposed to different concentrations of the prepared fractions (ranging from 50 to 200 μg/ml), and incubated for 48 h. After treatment, cells were incubated with MTT solution and the formazan crystals were solubilized and the absorbance was read at 570 nm[23].

    • Binding energies of phytochemicals retrieved from plant leave extract with STAT3 were calculated by using software InstaDock for molecular docking. Discovery Studio Visualizer, and PyMOL, were used to visualize the chemical interactions of ligands and proteins. SWISS-ADME tool and ProTox-II were used for pharmacokinetic profiling studies. The X-ray crystal structure of STAT3 (PDB ID: 6NJS) was downloaded from Protein Data Bank (PDB). All co-crystallized hetero atoms and attached water molecules and co-crystallized ligands, were eliminated from the original coordinates. The Polar hydrogen atoms were inculcated, the residue structures having lower occupancy were removed, and the incomplete side chains were then substituted by using ADT. Three-dimensional structures of phytocompounds were sketched using Chem3D.

    • Determination of the analogous behavior to the drug of phytoligands with the help of cheminformatics was done using online tool SwissADME developed by the Molecular Modelling Group, Swiss Institute of Bioinformatics[24]. The computation of pharmacokinetics and physicochemical molecular properties help medicinal chemists in their routine drug discovery processes. Significant basic molecular information can be excavated from the chemical structure. The methods were preferred over other methods because of the speed, but also for the ease of interpretating results by fingerprinting method to enable researchers move through translation to medicinal chemistry and in molecular designing[25].

    • The rationale behind molecular docking is to steer medicinal chemists for translational research. The affinity of a molecule to the receptor changes with small structural changes in the molecule[26]. For molecular docking, STAT3 core complex PDB id (PMID: 31715132) was remodeled to ascertain binding energies with the best conformational poses of Elaeocarpus ganitrus leaves phytoligands. The InstaDock software is used to dock phytoligands with blind search space having a grid size of 110, 70, and 108 Å for X, Y, Z coordinates, correspondingly. The center of the grid was confined to X: 63.09, Y: 14.98, and Z: −76.91 axis, which covers all the heavy atoms embedded in the protein. The conformational site selected was so that the movement of the ligands was free to probe their best binding coordinates. Default docking specifications were employed to calculate various parameters. All the docking conformational poses were generated using PyMOL, a molecular visualization system and Discovery Studio Predictor.

    • Physicochemical parameters, water solubility, lipophilicity, pharmacokinetics, and drug-likeness were elicited from SwissADME. To retrieve the toxicological profile of the phytoligands ProTox-II servers were employed[27]. Early estimation of the Absorption, Distribution, Metabolism, Excretion and Toxicity abbreviated as ADMET imperative to ascribe the pharmacodynamics success of the lead phytoligands. (SMILES) strings to encode chemical structures were imported from PubChem, open chemistry database and implemented in SWISS-ADME tool[24] to auspicate lipophilicity to show hydrophobic and hydrophilic nature, water solubility, necessary for absorption across membranes, and drug-likeness rules to assess metabolic profiles. Toxicology prediction of phytoligands is a crucial and fundamental aspect in the drug discovery process. ProTox-II is used to estimate computational toxicity, to accelerate the course to drug discovery, compute animal toxicity, and also help to attenuate animal experiments. In the PROTOX-II web server, toxicity classes are designated into four segments. Category I comprised of chemical entities with LD50 (LD = lethal Dose) values (LD50 ≤ 5) mg/kg, Category II comprised of compounds with LD50 values (5 < LD50 ≤ 50) mg/kg, Category III comprised of chemical entities having LD50 values (50 < LD50 ≤ 300) mg/kg, Category IV comprised of compounds which have LD50 values (300 < LD50 ≤ 2,000) mg/kg, Category V comprised of compounds with LD50 values (2,000 < LD50 ≤ 5,000) mg/kg and Category VI comprised of compounds showing LD50 values (LD50 > 5,000) mg/kg[28]. Category I and II manifested high toxicity, Category III and IV are comparatively less toxic and Category V and VI are considered to be non-toxic.

    • The solvent extraction technique is usually employed to prepare extracts from plant materials attributable to its convenience to operate. The importance lies in that a large amount of plant material can be extracted with minimal solvent[26]. Fresh leaves of Elaeocarpus ganitrus were purchased from Patanjali Herbal Garden Site Nursery located in Panchayanpur, Uttarakhand 249405, India. The confirmation of the authenticity of the Elaeocarpus ganitrus was done by the Department of Botany, Jamia Hamdard, New Delhi, India, and the leaf specimens deposited in the University. The crude methanol extract was unintermittedly fractionated in the solvents heptane, chloroform, ethyl acetate, methanol, and water according to their increasing polarity[16]. The anticancer activity of extracts was analyzed on the basis of their IC50 values. Cancerous HeLa cell line when treated with E. ganitrus leaf extracts exhibited a substantial inhibition of cells. The half maximal inhibitory concentration of chloroform and methanol extracts of E. ganitrus was (IC50 = 304.39 μg/ml) and (IC50 = 308.59 μg/ml) respectively followed by water (IC50 = 340.14 μg/ml), ethyl acetate (IC50 = 350.72 μg/ml) and heptane (IC50 = 381.76 μg/ml) extracts (Fig. 1ae & Fig. 2). The qualitative investigation using standard methodology[22] of chloroform and methanol fractions of E. ganitrus leaves disinterred the presence of major phytochemicals namely steroids, saponins, terpenoids, tannins, phenols, glycosides and flavonoids Table 1. GC-MS analysis of the chloroform and methanolic fractions was done based on their lowest half maximal inhibitory concentration to get a complete profiling of the plant compounds. The peaks in the total ion current (TIC) chromatogram of GC-MS profile of the phytoligands commensurate with the spectrum of known chemical databases stockpiled in the GC-MS library. The gas chromatogram depicts the relative concentrations of different phytoligands getting eluted according to the retention time. The heights of the peak represent the comparative concentrations of the compounds present in the plant appear as peaks at different m/z ratios. The components present with their retention time, molecular formula, molecular weight and concentration (peak area %) are provided in Tables 2 & 3 showing the presence of 56 and 50 bioactive phytochemicals in the chloroform and methanol extracts respectively. Of 106 phytoligands obtained from chloroform and methanol extracts of E. ganitrus leaves, 81 phytoligands were identified has having the best drug-like properties following Lipinski's rule of five. Lipinski's rule states that molecular properties, physical or chemical of a compound are significant for a drug's pharmacokinetics behavior inside a biological system. The drug molecules that go along with the RO5 have fewer attrition rates when undergoing clinical trials. The cheminformatics study to identify potential chemical entities having propensity for predefined biological targets is called virtual screening[28]. To endeavor in vitro experiments time diminution, molecular docking-based virtual screening of 81 selected compounds with two reference inhibitors having substantial binding energies with 6NJS were preferred for further analysis. The STAT3 has dual nature as an oncogene or as a tumor suppressor during cancer progression. It has a SH2 domain, linker domain, DNA binding domain, and all-alpha domain. The total energy of binding, Vander Waals forces, hydrogen bonding, electrostatic attraction, desolvation, and also a number of rotatable bonds present in the phytoligand, contribute to observe the free energy of binding of phytoligands with the receptor. Twenty-six (EG-1 to EG-26) compounds were selected as having appreciable binding affinities towards the 6NJS receptor (Table 4).

      Figure 1. 

      Effects of (a) heptane, (b) chloroform, (c) methanol, (d) ethyl acetate and (e) water fractions of E. ganitrus leaves on the human cancer cell lines HeLa using MTT assay.

      Figure 2. 

      IC50 values of different extracts of E. ganitrus leaves against human cancer cell lines HeLa.

      Table 1.  Qualitative analysis of phytochemicals in E. ganitrus leaf extracts.

      Tested compounds Chloroform extractMethanol extract
      Steroids++
      Terpenoids++
      Saponins++
      Glycosides++
      Tannins++
      Flavonoids++
      Phenols++
      + → Present; − → Absent.

      Table 2.  GC–MS analysis of chloroform fraction of E. ganitrus leaves.

      Peak no.R. TimeAreaArea %Name
      17.32833924513.60Phenol, 2-methoxy-4-(2-propenyl)-
      27.4944325420.46Cyclododecane
      37.92550188745.32Bicyclo[7.2.0]undec-4-ene,4,11,11-trimethyl-8-methylene-
      48.3922645730.281,4,8-Cycloundecatriene, 2,6,6,9-tetramethyl-,(e,e,e)-
      59.13715685461.66Phenol, 3,5-bis(1,1-dimethylethyl)-
      610.01627271302.891-Heptadecene
      712.22129599333.141-Octadecene
      812.6701789630.19Neophytadiene
      912.7821478930.162-Pentadecanone, 6,10,14-trimethyl-
      1013.4961275110.147,9-Di-tert-butyl-1-oxaspiro(4,5)deca-6,9-diene-2,8-dione
      1113.59421053742.23Hexadecanoic acid, methyl ester
      1213.8011860960.20Isophytol
      1314.00325829872.74Dibutyl phthalate
      1414.2331264020.141-Nonadecene
      1515.1432629120.281-Octadecanol
      1615.1964093250.439,12-Octadecadienoic acid (z,z)-, methyl ester
      1715.25614188331.509,12,15-Octadecatrienoic acid, methyl ester, (z,z,z)-
      1815.39682671938.76P-menth-1-ene-3,3-d2
      1915.77642684864.53Cholest-24-ene, (5.alpha.,20.xi.)-
      2016.08118351251.95Behenic alcohol
      2117.0155746820.61Glycidyl palmitate
      2217.5023567620.384,8,12,16-Tetramethylheptadecan-4-olide
      2317.79215482661.64N-tetracosanol-1
      2418.5078929660.95Glycidyl oleate
      2518.6333870750.41Pentacosane
      2618.88511608121.23Hexadecanoic acid, 2-hydroxy-1-(hydroxymethyl)ethyl ester
      2718.94925201262.671,2-Benzenedicarboxylic acid
      2819.38317228801.83Hexacosyl pentafluoropropionate
      2919.9972957890.31Carbonic acid, propyl 3,5-difluophenyl ester
      3020.15216559811.76Tetracontane
      3120.2877183500.769-Otadecenoic acid (z)-, 2,3-dihydroxypropyl ester
      3220.43716334241.73Octadecanoic acid, 2,3-dihydroxypropyl ester
      3320.87585309029.04Carbonic acid, eicosyl prop-1-en-2-yl ester
      3421.22511697911.24 .alpha.-tocospiro b
      3521.37114757101.56 .alpha.-tocospiro b
      3621.56627907932.96Tetracosane
      3721.6195665540.601-Heptacosanol
      3821.9782082850.22Tetracontane
      3922.24431468313.34Tetracontane
      4022.3391988570.21Triacontyl acetate
      4122.7586025580.64 .gamma.-tocopherol
      4222.98737696614.00Tetracontane
      4323.0833075630.33Octacosanol
      4423.3688043680.852,5,7,8-Tetramethyl-2-(4,8,12-trimethyltridecyl)-3,4-dihydro-2h-chromen-6-yl hexofuranoside
      4523.83228257813.00Hexatriacontane
      4624.4422247490.24Ergost-5-en-3-ol
      4724.6971463260.162,6,10,15,19,23-Hexamethyl-tetracosa-2,10,14,18,22-pentaene-6,7-diol
      4824.81625792062.73Tetracontane
      4925.34444036834.67 .gamma.-sitosterol
      5025.8972474510.26Phenol, 2,4-bis(1,1-dimethylethyl)-, phosphite (3:1)
      5125.97110509901.11Tetracontane
      5227.36910622631.13Tetracontane
      5329.04944683384.74Benzenepropanoic acid, 3,5-bis(1,1-dimethylethyl)-4-hydroxy-,octadecyl ester
      5431.0618597260.91Tetrapentacontane
      5533.5056890170.73Tetrapentacontane
      5636.5155702640.60Tetrapentacontane

      Table 3.  GC–MS analysis of methanol fraction of E. ganitrus leaves.

      Peak no.R. timeAreaArea %Name
      14.56217160622.754h-pyran-4-one, 2,3-dihydro-3,5-dihydroxy-6-methyl-
      25.5451328910.211,5-Dimethyl-1-vinyl-4-hexenyl 2-aminobenzoate
      36.130664200.11E-6-octadecen-1-ol acetate
      46.658982694115.754-Hydroxy-3-methylacetophenone
      57.421592600.091-Undecanol
      67.7462900200.46Methyl2,3,6,7-tetra-o-acetyl-4-o-methyl-.beta.-glycero-d-glucoheptopyranoside
      78.93519327693.10Guanosine
      89.1785358220.861,3:2,5-Dimethylene-l-rhamnitol
      99.9494638700.74Octadecanoic acid
      1010.1442790410.451,2-Benzenedicarboxylic acid, diethyl este
      1110.37016839662.70 .alpha.-methyl-l-sorboside
      1210.60613007272.08 .alpha.-d-galactopyranoside, methyl
      1310.9202805730.45Butanoic acid, 3-methyl-, hexahydro-4- methylspiro[cyclopenta[c]pyran-7(1h),2'-oxirane]-1,6-diyl ester
      1411.090803800.13Tricyclo[7.2.0.0(2,6)]undecan-5-ol, 2,6,10,10-tetramethyl- (isomer 3)
      1511.2241616060.26 .alpha.-d-galactopyranoside, methyl
      1611.4921894850.30Octadecanoic acid, methyl ester
      1712.4372606340.422(4h)-benzofuranone, 5,6,7,7a-tetrahydro-6- hydroxy-4,4,7a-trimethyl-, (6s-cis)-
      1812.6431237680.20Neophytadiene
      1913.56528660454.59Hexadecanoic acid, methyl ester
      2013.780340440.051-hexadecen-3-ol, 3,5,11,15-tetramethyl-
      2113.910475350.08Silane, ethenylethyldimethyl-
      2214.555735350.12Pentadecanoic acid, methyl ester
      2315.18817688992.839,12-Octadecadienoic acid (z,z)-, methyl ester
      2415.24959714079.57(9e,12e)-9,12-octadecadienoyl chloride #
      2515.385793321212.711,1'-Bicyclohexyl, 2-methyl-, cis-
      2615.4819100731.46Methyl stearate
      2715.758823962913.21Cholest-24-ene, (5.alpha.,20.xi.)-
      2816.0753964310.64Methyl octadeca-9,12-dienoate
      2916.444772460.12Methyl 4-(dimethylamino)bicyclo[2.2.2]oct- 5-ene-2-carboxylate
      3016.612548090.09Hexadecanoic acid, 2-hydroxy-1-(hydroxymethyl)ethyl ester
      3116.9992156190.3517-octadecynoic acid
      3217.2492190610.35Eicosanoic acid, methyl ester
      3318.1272324800.37Oleoyl chloride
      3418.5093011250.48Undec-10-ynoic acid, undec-2-en-1-yl ester
      3518.7051353320.22Hexadecanoic acid, 1-(hydroxymethyl)-1,2-ethanediyl ester
      3618.89946078647.38Hexadecanoic acid, 2-hydroxy-1-(hydroxymethyl)ethyl ester
      3719.6551063120.17Hexadecanoic acid, methyl ester
      3820.3005806710.93Oleoyl chloride
      3920.4629943041.59Octadecanoic acid, 2,3-dihydroxypropyl ester
      4020.91235840275.749-octadecenamide
      4121.2252214500.35 .alpha.-tocospiro b
      4221.3783841650.62 .alpha.-tocospiro b
      4321.6261759730.28Eicosyl heptafluorobutyrate
      4421.803899870.14Hexacosanoic acid, methyl ester
      4522.7691609600.26 .gamma.-tocopherol
      4622.989867570.14Tetracontane
      4723.174924420.15Stigmast-5-en-3-ol, (3.beta.)-
      4823.38010301431.65Vitamin e
      4925.37212372631.98 .gamma.-sitosterol
      5027.0841838260.29Di-o-acetyltetrahydrostapelogenin

      Table 4.  Docking results of 81 phytoligands.

      S. no.Name of the ligandBinding free energy
      (kcal/mol)
      pKiLigand efficieny (kcal/mo/non-H atom)Torsional energy
      1Phenol, 2-methoxy-4-(2-propenyl)-–5.64.110.46671.2452
      2Cyclododecane–5.94.330.49170
      3Bicyclo[7.2.0]undec-4-ene,4,11,11-trimethyl-8-methylene-–6.64.840.440
      41,4,8-Cycloundecatriene, 2,6,6,9-tetramethyl-,(e,e,e)-–6.54.770.43330
      5Phenol, 3,5-bis(1,1-dimethylethyl)-–6.54.770.43330.9339
      61-Heptadecene–4.63.370.27064.3582
      71-Octadecene–42.930.24.0469
      8Neophytadiene–5.94.330.2954.0469
      92-Pentadecanone, 6,10,14-trimethyl-–5.43.960.28423.7356
      107,9-Di-tert-butyl-1-oxaspiro(4,5)deca-6,9-diene-2,8-dione–6.54.770.3250.6226
      11Hexadecanoic acid, methyl ester–4.93.590.25794.6695
      12Isophytol–4.93.590.23334.3582
      13Dibutyl phthalate–5.23.810.263.113
      141-Nonadecene–4.93.590.25794.9808
      151-octadecanol–53.670.26325.2921
      169,12-Octadecadienoic acid (z,z)-, methyl ester–53.670.23814.6695
      179,12,15-Octadecatrienoic acid, methyl ester, (z,z,z)-–5.43.960.25714.3582
      18P-Menth-1-ene-3,3-d2–4.93.590.490.3113
      19Behenic alcohol–4.73.450.20436.5373
      20Glycidyl palmitate–5.43.960.28423.7356
      214,8,12,16-Tetramethylheptadecan-4-olide–6.34.620.27393.7356
      22N-tetracosanol-1–4.43.230.1767.1599
      23Glycidyl oleate–4.63.370.19175.6034
      24Pentacosane–4.73.450.1886.8486
      25Hexadecanoic acid, 2-hydroxy-1-(hydroxymethyl)ethyl ester–4.63.370.26.226
      261,2-Benzenedicarboxylic acid–5.74.180.4751.2452
      27Carbonic acid, propyl 3,5-difluophenyl ester–6.14.470.40671.5565
      289-Octadecenoic acid (z)-, 2,3-dihydroxypropyl ester–53.670.26.5373
      29Octadecanoic acid, 2,3-dihydroxypropyl ester–4.93.590.1966.8486
      30Carbonic acid, eicosyl prop-1-en-2-yl ester–5.13.740.18896.8486
      31Tetracosane–4.93.590.20426.5373
      321-Heptacosanol–5.13.740.18218.0938
      33Triacontyl acetate–3.92.860.11479.339
      34Gamma.-tocopherol–6.84.990.22674.0469
      35Octacosanol–4.23.080.14488.4051
      362,5,7,8-Tetramethyl-2-(4,8,12-trimethyltridecyl)-3,4-dihydro-2h-chromen-6-yl hexofuranoside–7.25.280.17146.226
      37Ergost-5-en-3-ol–7.35.350.25171.8678
      382,6,10,15,19,23-Hexamethyl-tetracosa-2,10,14,18,22-pentaene-6,7-diol–64.40.18755.6034
      39Gamma.-sitosterol–96.60.32.1791
      404h-pyran-4-one,2,3-dihydro-3,5-dihydroxy-6-methyl-–53.670.50.6226
      411,5-Dimethyl-1-vinyl-4-hexenyl 2-aminobenzoate–6.44.690.322.4904
      42E-6-octadecen-1-ol acetate–4.73.450.21365.2921
      434-Hydroxy-3-methylacetophenone–5.74.180.51820.6226
      441-undecanol–4.53.30.3753.113
      45Methyl2,3,6,7-tetra-o-acetyl-4-o-methyl-.beta.-glycero-d-glucoheptopyranoside–5.54.030.19643.7356
      46Guanosine–6.84.990.26151.5565
      471,3:2,5-Dimethylene-l-rhamnitol–5.43.960.41540.3113
      48Octadecanoic acid–5.33.890.2655.2921
      491,2-benzenedicarboxylic acid, diethyl este–5.43.960.33751.8678
      50 .alpha.-methyl-l-sorboside–4.73.450.36151.8678
      51 .alpha.-d-galactopyranoside, methyl–5.23.810.41.8678
      52Butanoic acid, 3-methyl-, hexahydro-4- methylspiro[cyclopenta[c]pyran-7(1h),2'-oxirane]-1,6-diyl ester–6.84.990.22673.4243
      53Tricyclo[7.2.0.0(2,6)]undecan-5-ol, 2,6,10,10-tetramethyl- (isomer 3)–6.54.770.40620.3113
      54 .alpha.-d-galactopyranoside, methyl–5.33.890.40771.8678
      55Octadecanoic acid, methyl ester–4.13.010.19525.2921
      562(4h)-benzofuranone, 5,6,7,7a-tetrahydro-6- hydroxy-4,4,7a-trimethyl-, (6s-cis)-–6.54.770.46430.3113
      57Neophytadiene–53.670.254.0469
      58Hexadecanoic acid, methyl ester–4.93.590.25794.6695
      591-Hexadecen-3-ol, 3,5,11,15-tetramethyl-–5.74.180.27144.3582
      60Pentadecanoic acid, methyl ester–4.43.230.24444.3582
      619,12-Octadecadienoic acid (z,z)-, methyl ester–5.43.960.25714.6695
      62(9e,12e)-9,12-octadecadienoyl chloride #–4.73.450.2354.3582
      631,1'-bicyclohexyl, 2-methyl-, cis-–5.64.110.43080.3113
      64Methyl stearate–53.670.23815.2921
      65Cholest-24-ene, (5.alpha.,20.xi.)-–9.26.750.34071.2452
      66Methyl octadeca-9,12-dienoate–4.53.30.21434.6695
      67Methyl 4-(dimethylamino)bicyclo[2.2.2]oct- 5-ene-2-carboxylate–5.74.180.380.9339
      68Hexadecanoic acid, 2-hydroxy-1-(hydroxymethyl)ethyl ester–4.83.520.20876.226
      6917-octadecynoic acid–5.13.740.2555.2921
      70Eicosanoic acid, methyl ester–4.83.520.20875.9147
      71Undec-10-ynoic acid, undec-2-en-1-yl ester–5.13.740.21255.9147
      72Hexadecanoic acid, 2-hydroxy-1-(hydroxymethyl)ethyl ester–4.73.450.20436.226
      73Hexadecanoic acid, methyl ester–4.53.30.23684.6695
      74Oleoyl chloride–4.63.370.234.6695
      75Octadecanoic acid, 2,3-dihydroxypropyl ester–4.63.370.1846.8486
      769-octadecenamide–4.73.450.2354.6695
      77 .alpha.-tocospiro b–6.44.690.19394.3582
      78Eicosyl heptafluorobutyrate–5.64.110.16977.1599
      79Hexacosanoic acid, methyl ester–4.73.450.16217.7825
      80Stigmast-5-en-3-ol, (3.beta.)-–7.65.570.25332.1791
      81Vitamin e–7.15.210.2294.0469
      82Plumbagin–5.94.330.42140.3113
      83Sanguinarine–8.96.530.3560

      The absorption of drugs by the body is related to their pharmacokinetic properties and also cellular toxicity. The potency of the drug depends mostly on the pharmacokinetic parameters because ADME processes command the rate and extent of absorption when an administered dose of a drug approaches to its action site. Hence, in silico pharmacokinetic profile of filtered compounds was surveyed to gather the putative bioavailability data for receptor 6NJS. The cumulative findings for pharmacokinetics profiling, bioavailability data, drug-likeness properties and drug friendliness and toxicity effects of selected 26 phytoligands with known inhibitors (Plumbagin and Sanguinarine) are given in Tables 510. The prediction revealed that the six molecules (EG-9, EG-12, EG-13, EG-16, and EG-26) can be lead compounds for new drug candidates for anti-cancer phytomedicine. The Half maximal Inhibitory concentration of EG-13 was (IC50 = 254.29 µg/ml) further support our results (Fig. 3).

      Table 5.  Pharmacokinetics prediction of phytoligands established in E. ganitrus.

      S. no.Phytochemical
      Gastro- intestinal absorptionBlood-brain permeantP-glycoprotein substrateCYP450 1A2 inhibitorCYP450 2C19 inhibitorCYP450 2C9 inhibitorCYP450 2D6 inhibitorCYP450 3A4 inhibitorSkin permeation as log Kp (cm/s)
      EG-1Cholest-24-ene, (5.alpha.,20.xi.)-LowNoNoNoNoYesNoNo–1.02
      EG-2gamma.-sitosterolLowNoNoNoNoNoNoNo–2.65
      EG-3Stigmast-5-en-3-ol, (3.beta.)-LowNoNoNoNoNoNoNo–2.20
      EG-4Ergost-5-en-3-olLowNoNoNoNoNoNoNo–2.50
      EG-52,5,7,8-Tetramethyl-2-(4,8,12-trimethyltridecyl)-3,4-dihydro-2h-chromen-6-yl hexofuranosideLowNoNoNoNoNoNoYes–3.60
      EG-6Vitamin eLowNoYesNoNoNoNoNo–1.33
      EG-7GuanosineLowNoNoNoNoNoNoNo–9.37
      EG-8gamma.-tocopherolLowNoYesNoNoNoNoNo–1.51
      EG-9Butanoic acid, 3-methyl-, hexahydro-4- methylspiro[cyclopenta[c]pyran-7(1h),2'-oxirane]-1,6-diyl esterHighYesNoNoNoNoYesYes–6.18
      EG-10Bicyclo[7.2.0]undec-4-ene,4,11,11-trimethyl-8-methylene-LowNoNoNoYesYesNoNo–4.44
      EG-111,4,8-Cycloundecatriene, 2,6,6,9-tetramethyl-,(e,e,e)-LowNoNoNoNoYesNoNo–4.32
      EG-12Phenol, 3,5-bis(1,1-dimethylethyl)-HighYesNoNoNoNoYesNo–4.07
      EG-137,9-Di-tert-butyl-1-oxaspiro(4,5)deca-6,9-diene-2,8-dioneHighYesNoNoYesYesNoNo–5.28
      EG-142(4h)-benzofuranone, 5,6,7,7a-tetrahydro-6- hydroxy-4,4,7a-trimethyl-, (6s-cis)-HighYesNoNoNoNoNoNo–6.79
      EG-15Tricyclo[7.2.0.0(2,6)]undecan-5-ol, 2,6,10,10-tetramethyl- (isomer 3)HighYesNoNoYesYesNoNo–4.75
      EG-161,5-Dimethyl-1-vinyl-4-hexenyl 2-aminobenzoateHighYesNoNoYesYesNoNo–4.54
      EG-17alpha.-tocospiro bHighNoNoNoNoNoNoNo–3.90
      EG-184,8,12,16-Tetramethylheptadecan-4-olideLowNoNoYesNoYesNoNo–2.70
      EG-19Carbonic acid, propyl 3,5-difluophenyl esterHighYesNoYesYesNoNoNo–5.37
      EG-202,6,10,15,19,23-Hexamethyl-tetracosa-2,10,14,18,22-pentaene-6,7-diolLowNoNoYesNoYesNoNo–2.37
      EG-21CyclododecaneLowNoNoNoNoNoNoNo–4.42
      EG-22NeophytadieneLowNoYesNoNoYesNoNo–1.17
      EG-231,2-Benzenedicarboxylic acidHighNoNoNoNoNoNoNo–6.80
      EG-244-Hydroxy-3-methylacetophenoneHighYesNoYesNoNoNoNo–6.54
      EG-251-Hexadecen-3-ol, 3,5,11,15-tetramethyl-LowNoYesNoNoYesNoNo–2.41
      EG-26Methyl 4-(dimethylamino)bicyclo[2.2.2]oct- 5-ene-2-carboxylateHighYesNoNoNoNoNoNo–6.65
      PlumbaginHighYesNoYesNoNoNoNo–5.82
      SanguinarineHighYesYesYesYesNoNoNo–5.17

      Table 6.  Bioavailability prediction of phytoligands established in E. ganitrus.

      Phyto-ligands Bioavailability scoreWater solubility as logSiLOGPXLOGP3WLOGPMLOGPSILICOS-IT
      EG-10.55Poorly soluble as –6.255.1210.628.428.327.14
      EG-20.55Poorly soluble as –6.194.758.867.965.807.04
      EG-30.55Poorly soluble as –6.194.799.348.026.737.04
      EG-40.55Moderately soluble as –5.794.928.807.636.546.63
      EG-50.55Poorly soluble as –7.376.148.896.313.498.12
      EG-60.55Poorly soluble as –9.165.9210.708.846.149.75
      EG-70.55Very Soluble as 0.51–0.23–1.89–3.00–2.76–2.22
      EG-80.55Poorly soluble as –8.795.7610.338.535.949.20
      EG-90.55Soluble as –2.863.873.342.932.073.34
      EG-100.55Soluble as –3.773.294.384.734.634.19
      EG-110.55Soluble as –3.523.274.555.044.533.91
      EG-120.55Soluble as –4.252.864.913.993.873.81
      EG-130.55Soluble as –3.812.913.813.592.873.82
      EG-140.55Very Soluble as –1.821.881.001.411.491.86
      EG-150.55Soluble as –3.183.014.093.613.813.40
      EG-160.55Moderately soluble as –4.283.374.834.123.633.75
      EG-170.55Poorly soluble as –7.194.947.246.583.677.85
      EG-180.55Poorly soluble as –6.314.157.866.524.966.99
      EG-190.55Soluble as –3.592.843.173.732.912.81
      EG-200.55Poorly soluble as –6.306.119.388.776.019.10
      EG-210.55Soluble as –3.213.014.104.685.004.00
      EG-220.55Poorly soluble as –6.115.059.627.176.217.30
      EG-230.85Soluble as –1.140.600.731.081.200.61
      EG-240.55Very Soluble as –2.531.540.951.901.442.14
      EG-250.55Moderately soluble as –5.514.978.026.365.256.57
      EG-260.55Very Soluble as –1.352.701.311.451.771.11
      Plumbagin0.55Soluble as –2.851.792.291.720.592.22
      Sanguinarine0.55Poorly soluble as –6.09–0.044.453.432.723.85

      Table 7.  Drug-likeness prediction of phytoligands established in E. ganitrus.

      Phyto-ligandsLipinski
      rule
      Ghose
      filter
      Veber
      filter
      Egan
      filter
      Muegge
      filter
      EG-1YesNoYesNoNo
      EG-2YesNoYesNoNo
      EG-3YesNoYesNoNo
      EG-4YesNoYesNoNo
      EG-5YesNoNoNoNo
      EG-6YesNoNoNoNo
      EG-7YesNoNoNoNo
      EG-8YesNoNoNoNo
      EG-9YesYesYesYesYes
      EG-10YesYesYesYesYes
      EG-11YesYesYesYesYes
      EG-12YesYesYesYesNo
      EG-13YesYesYesYesYes
      EG-14YesYesYesYesNo
      EG-15YesYesYesYesNo
      EG-16YesYesYesYesYes
      EG-17YesNoNoNoNo
      EG-18YesNoNoNoNo
      EG-19YesYesYesYesYes
      EG-20YesNoNoNoNo
      EG-21YesYesYesYesNo
      EG-22YesNoNoNoNo
      EG-23YesNoYesYesNo
      EG-24YesNoYesYesNo
      EG-25YesNoNoNoNo
      EG-26YesYesYesYesYes
      PlumbaginYesYesYesYesNo
      SanguinarineYesYesYesYesYes

      Table 8.  Medicinal chemistry prediction of phytoligands established in E. ganitrus.

      SI. No.PAINS structural alertBrenk structural alertLead-
      likeness
      Synthetic accessibility score
      EG-10125.61
      EG-20126.42
      EG-30126.30
      EG-40126.17
      EG-50037.10
      EG-60035.17
      EG-70003.86
      EG-80035.00
      EG-90225.59
      EG-100124.51
      EG-110123.66
      EG-120021.37
      EG-130014.35
      EG-140013.63
      EG-150023.77
      EG-160212.91
      EG-170036.76
      EG-180024.12
      EG-190112.23
      EG-200135.52
      EG-210022.21
      EG-220124.08
      EG-230011.00
      EG-240011.00
      EG-250123.89
      EG-260114.38
      Plumbagin2012.41
      Sanguinarine0212.59

      Table 9.  Toxicity prediction of phytoligands established in E. ganitrus.

      Phyto-ligandsLD50 (mg/kg)Toxicity classHepatotoxicityCarcinogenicityImmunotoxicityMutagenicityCytotoxicity
      EG-150005InactiveInactiveActiveInactiveInactive
      EG-28904InactiveInactiveActiveInactiveInactive
      EG-38904InactiveInactiveActiveInactiveInactive
      EG-48904InactiveInactiveActiveInactiveInactive
      EG-530005InactiveInactiveActiveInactiveInactive
      EG-650005InactiveInactiveInactiveInactiveInactive
      EG-7132InactiveInactiveInactiveInactiveInactive
      EG-850005InactiveInactiveInactiveInactiveInactive
      EG-980006InactiveActiveInactiveActiveInactive
      EG-1053005InactiveInactiveActiveInactiveInactive
      EG-1136505InactiveInactiveInactiveInactiveInactive
      EG-128004InactiveInactiveInactiveInactiveInactive
      EG-139004InactiveInactiveInactiveInactiveInactive
      EG-14342InactiveActiveInactiveInactiveInactive
      EG-1520505InactiveInactiveInactiveInactiveInactive
      EG-1642505InactiveInactiveInactiveInactiveInactive
      EG-173003InactiveInactiveInactiveInactiveActive
      EG-1844005InactiveInactiveInactiveInactiveInactive
      EG-1915004InactiveInactiveInactiveInactiveInactive
      EG-2043005InactiveInactiveInactiveInactiveInactive
      EG-217503InactiveActiveInactiveInactiveInactive
      EG-2250506InactiveInactiveInactiveInactiveInactive
      EG-2325305InactiveInactiveInactiveInactiveInactive
      EG-2428305InactiveInactiveInactiveInactiveInactive
      EG-253404InactiveInactiveInactiveInactiveInactive
      EG-2620004InactiveInactiveInactiveInactiveInactive
      Plumbagin162InactiveActiveInactiveActiveInactive
      Sanguinarine7784InactiveActiveActiveActiveInactive

      Table 10.  Bioavailability prediction of phytoligands established in E. ganitrus.

      Phyto-ligandLipophilicity
      (XLOGP3)
      Size
      (MW g/mol)
      Polarity
      (TPSA)
      Insolubility
      [Log S (ESOL)]
      Insaturation
      (Fraction Csp3)
      Flexibility
      (Num. rotatable bonds)
      EG-93.34368.4674.36–3.700.908
      EG-124.91206.3220.23–4.380.572
      EG-133.81276.3743.37–3.820.652
      EG-154.09222.3720.23–3.801.000
      EG-164.83273.3752.32–4.340.357
      EG-261.31209.2829.54–1.760.753
      Plumbagin2.29188.1854.37–2.770.090
      Sanguinarine4.45332.3340.80–5.240.150

      Figure 3. 

      IC50 values of EG-13 phytochemical of E. ganitrus leaves against human cancer cell lines HeLa.

      In Table 5, for pharmacokinetics prognostication, the gastrointestinal (GI) absorption rate was fetched for all preferred six phytoligands and both reference drugs. The blood-brain permeability was seen as positive for all the six phytoligands and both reference drugs. The prediction of bioavailability (Table 6) demonstrated that similar bioavailability scores were observed for all the filtered six phytoligands (0.55) like reference drugs. The water solubility data showed all the six compounds and plumbagin are soluble while Sanguinarine is poorly soluble. For drug-likeness prediction (Table 7), all the six compounds and both known inhibitors were obtained suitable for the Lipinski rule as zero violation. For Ghose, Veber, and Egan filter 0 violation was obtained for all the six phytoligands and both inhibitors. In the case of medicinal chemistry friendliness prediction (Table 8), the PAINS structural alert obtained 0 violations for all the six phytoligands and sanguinarine while two alerts for plumbagin. Table 9 shows EG-9 belongs to the non-toxic class VI, EG-15, and EG-16 also belong to the non-toxic class V, EG-12, EG-13, EG-26 and Sanguinarine belongs to the less toxic class IV while plumbagin belongs to the high-toxic class II. The bioavailability radar (Fig. 3) for phytoligands depicting bioavailability prognostic showed that all six phytoligands were found within the data range of lipophility nature (−0.7 < XLOGP3 < +5.0), molecule size (150 g/mol < MW < 500g/mol), polarity (20 Ų < TPSA < 130Ų), insolubility [−6 < LogS (ESOL) < 0], insaturation (0.25 < Fraction Csp3 < 1) and flexible bonds (0 < Num. rotatable bonds < 9) and colored part of radar while known inhibitors plumbagin and sanguinarine does not fit the bioavailability radar (Table 10). As mentioned in Table 5, all the phytoligands and reference compounds have higher gastrointestinal (GI) absorption rates, therefore they can instantly be absorbed by the human intestine. All phytoligands have the ability to pass the blood brain barrier (BBB permeant) and values for the aqueous solubility (log S) of the phytochemicals fall in the recommended range that is −1 to −5[29], thus, have improved absorption and distribution properties. The bioavailability scores were identical for all six molecules, standing at 0.55, similar to the reference drugs. In drug-likeness prediction, none of the six compounds and both known inhibitors violated the Lipinski rule, Ghose, Veber, and Egan filters. Regarding medicinal chemistry friendliness, the PAINS structural alert identified zero violations for all six phytoligands and Sanguinarine, whereas Plumbagin had two alerts. Table 9 revealed that EG-9 belonged to the non-toxic class VI, while EG-15 and EG-16 were in harmless class V. Other compounds EG-12, EG-13, EG-26 and sanguinarine was from less harmful class IV which could be modified to a non-toxic class during the lead optimization stage of drug discovery[30] while selected standard plumbagin showed high toxic class II. Drug-induced hepatotoxicity often lead to abrupt liver failure and drug rejections[31]. Drug-induced liver injury might be long-term or occur only once. Obviously, the selected compounds and standards are non-hepatotoxic. The bioavailability radar (Fig. 4) depicted that all six phytoligands were within the data range for oral bioavailability prediction. Conversely, standards plumbagin and sanguinarine did not fit within the bioavailability radar. The pink area shown in the radar corresponds to the most promising zone for all the bioavailability properties. In Table 10, all the phytochemicals satisfied 150 g/mol and 500 g/mol criteria for (SIZE) of good drug candidates. The polarity (POLAR) was observed with the Total Polarity Surface Area (TPSA) and all the phytochemicals show good TPSA values. Besides, the flexibility (FLEX) property evaluated by the number of rotatable bonds falls within the recommended range. Lipophilicity (LIPO) and insolubility (INSOLU) were evaluated and come in the range The Unsaturation (INSATU) was calculated using Fraction Csp3 falls within a recommended range of 0.25 < Fraction Csp3 < 1) for all phytoligands. However, Plumbagin and Sanguinarine exhibit lower values (0.09 and 0.15, respectively).

      Figure 4. 

      Bioavailability radar (pink area exhibits optimal range of particular property) for leading phytocompounds molecules. LIPO = lipophilicity as XLOGP3, SIZE = size as molecular weight, POLAR = polarity as TPSA (topological polar surface area), INSOLU = insolubility in water by log S scale, INSATU = insaturation as per fraction of carbons in the sp3 hybridization, and FLEX = flexibility as per rotatable bonds.

      2D and 3D interactions of the five phytoligands (EG-9, EG-12, EG-13, EG-15, EG-16 and EG-26) with 6njs are shown in Table 11. EG-9 divulged two assenting hydrogen bond interactions at the active site having amino acids of Glu96 and Lys97. In additon to that a non-classical C-H bond Vander Waals interaction was also noticed at the active site involving Arg93 residue and alkyl and pi-alkyl interactions were observed at Leu525 and Trp501 respectively. In EG-12 a conventional hydrogen bond interaction was observed at Asn538, a pi-pi T-shaped, two alkyl and a pi-alkyl interactions were observed at Tyr539, Ile522, Trp501 and Leu525 respectively. EG-13 showed one favorable hydrogen bond interaction and two hydrophobic alkyl interactions at the active site with the residue of Glu96, Leu95 and Lys97 respectively. EG-15 showed two alkyl and two pi-alkyl interactions at the active site of the residues of Leu95, Ile522, Trp501 and Tyr539 respectively. In EG-16 two conventional hydrogen bonds were observed at Leu731 and Thr716. EG-26 formed three favorable hydrogen bonds with Asp369, Asp370 and Asp371 at the active site of the receptor. Plumbagin showed a conventional hydrogen bond interaction, a pi-pi T-shaped and a carbon-hydrogen bond interaction at Tyr539, Trp501 and Ser540 respectively. Sanguinarine showed a carbon hydrogen bond, a pi-sigma, a alkyl, and a pi-alkyl interaction at the site of Glu696, Leu731, Pro769 and Pro695 respectively (Table 11). Previously it has been shown that residue at 97 could have amprospective ubiquitin acceptor position in STAT3 NH2 terminal domain, suggesting lysine amino acid may have a significant role and location in a sumolation/ubiquitination consensus sequence[32]. The majority of phytoligand interactions exist in the Linker domain and Transactivation domain of the STAT3.

      Table 11.  2D and 3D binding interactions between the receptor 6NJS and molecules.

      Phyto-ligands2D- Binding interaction3D- Binding interaction
      EG-9 (-6.8)
      EG-12 (-6.5)
      EG-13 (-6.5)
      EG-15 (-6.5)
      EG-16 (-6.4)
      EG-26
      (-5.7)
      Plumbagin
      Sanguinarine
    • All the six compounds (EG-9, EG-12, EG-13, EG-15, EG-16 and EG-26) significantly bind with STAT3. The phytochemicals epitomized good in silico results as reflected by their promising binding affinity, considerable inhibitory constant with optimum protein-ligand stabilization energy. Consecutively, binding signifies that phytoligands interact with STAT3 by the NH2 terminal and boosts its transcriptional activity and interferes with the cellular proliferation process and apoptosis[32]. Bioavailability radar and toxicological profiles of the preferred phytoligands revealed that these compounds compel to have ample drug likeliness properties. Moreover, EG-9, EG-13, EG-15, EG-16 and EG-26 have not been explored for their anticancer potential and can be derivatized or have the probability of being used as lead compounds.

    • The authors confirm contribution to the paper as follows: study design and draft manuscript preparation (equal): Mehnaj, Bhat AR, Athar F; supervision: Athar F; experimentation and writing of manuscript: Mehnaj; characterization and editing: Bhat AR. All authors reviewed the results and approved the final version of the manuscript.

    • This study involved the use of established human cell lines. The cell lines used in this research were obtained from the National Centre for Cell Sciences (NCCS), Pune, India and were used in accordance with institutional and national ethical standards. The cell lines have been previously published or validated, and no new human tissues were used in this study.

    • The supplementary data will be made available by the authors to all upon reasonable request.

    • Miss Mehnaj is grateful to UGC for obtaining the non-NET fellowship allowing completion of this work.

      • The authors declare that they have no conflict of interest.

      • Copyright: © 2024 by the author(s). Published by Maximum Academic Press, Fayetteville, GA. This article is an open access article distributed under Creative Commons Attribution License (CC BY 4.0), visit https://creativecommons.org/licenses/by/4.0/.
    Figure (4)  Table (11) References (32)
  • About this article
    Cite this article
    Mehnaj, Bhat AR, Athar F. 2024. In silico exploration of Elaeocarpus ganitrus extract phytochemicals on STAT3, to assess their anticancer potential. Medicinal Plant Biology 3: e009 doi: 10.48130/mpb-0024-0010
    Mehnaj, Bhat AR, Athar F. 2024. In silico exploration of Elaeocarpus ganitrus extract phytochemicals on STAT3, to assess their anticancer potential. Medicinal Plant Biology 3: e009 doi: 10.48130/mpb-0024-0010

Catalog

  • About this article

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return