Abstract

Background. Glioma is the most common primary intracranial tumor in adult patients. Among them, glioblastoma is a highly malignant one with a poor prognosis. Flavonoids are a class of phenolic compounds widely distributed in plants and have many biological functions, such as anti-inflammatory, antioxidant, antiaging, and anticancer. Nowadays, flavonoids have been applied to the therapy of glioma; however, the molecular mechanism underlying the therapeutic effects has not been fully elaborated. This study was carried out to explore the mechanism of selected active flavonoid compounds in treating glioma using network pharmacology and molecular docking approaches. Methods. Active ingredients and associated targets of flavonoids were acquired by using the Traditional Chinese Medicine Database and Analysis Platform (TCMSP) and Swiss TargetPrediction platform. Genes related to glioma were obtained from the GeneCards and DisGeNET databases. The intersection targets between flavonoid targets and glioma-related genes were used to construct protein-protein interaction (PPI) network via the STRING database, and the results were analyzed by Cytoscape software. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed and displayed by utilizing the Metascape portal and clusterProfiler R package. Molecular docking was carried out by iGEMDOCK and SwissDock, and the results were visually displayed by UCSF Chimera software. Results. Eighty-four active flavonoid compounds and 258 targets overlapped between flavonoid targets and glioma-related genes were achieved. PPI network revealed potential therapeutic targets, such as AKT1, EGFR, VEGFA, MAPK3, and CASP3, based on their node degree. GO and KEGG analyses showed that core targets were mainly enriched in the PI3K-Akt signaling pathway. Molecular docking simulation indicated that potential glioma-related targets-MAPK1 and HSP90AA1 were bounded more firmly with epigallocatechin-3-gallate (EGCG) than with quercetin. Conclusions. The findings of this study indicated that selected active flavonoid compounds might play therapeutic roles in glioma mainly through the PI3K-Akt signaling pathway. Moreover, EGCG had the potential antiglioma activity by targeting MAPK1 and HSP90AA1.

1. Introduction

Glioma, originated from the neuroepithelium, accounts for 40%∼50% of brain tumors and is one of the most common primary intracranial tumors; among them, glioblastoma is a highly malignant one with poor clinical outcome [1]. Extensive studies have shown that flavonoids have a good therapeutic effect on glioma [24]; however, the underlying therapeutic molecular mechanisms of flavonoids on glioma are not stated clearly. Hence, a systematic exploration of the molecular mechanisms of flavonoids on glioma is critical. Network pharmacology is an emerging interdisciplinary discipline and has been applied to comprehensively analyze the functional mechanisms of traditional Chinese medicine [5]. It is also used to reveal the active ingredients of natural medicine treating glioma [6].

In this study, we tried to systematically identify the molecular mechanisms of flavonoids’ antiglioma effects based on findings from network pharmacology and molecular docking. The flowchart of this study is shown in Figure 1. Our work portrays the ground view of antiglioma molecular mechanisms of flavonoids, which provided active compounds and therapy targets curing glioma.

2. Materials and Methods

2.1. Screening Active Ingredients and Predicting Related Targets

TCMSP (https://tcmsp-e.com/) is a unique systematic pharmacology platform for Chinese herbal medicines and is characterized by exploring relationships between compounds, targets, and diseases. TCMSP has been utilized to screen active flavonoid ingredients [7]. Flavonoids are categorized according to their molecular structures into flavones, flavonols, isoflavones, chalcones, flavones, and anthocyanidins [4], as shown in Table 1. Key parameters were taken into account, such as oral bioavailability (OB), drug-likeness (DL), and blood-brain barrier (BBB). OB is a measurement of the proportion of drugs entering the blood circulation. DL indicates the potential of a compound to be developed into a therapeutic drug with respect to its physical and chemical properties. BBB impedes drug distribution between blood and brain [8], which is a very important parameter in the treatment of glioma. In general, one compound is considered permeable across the BBB when its BBB permeability is larger than −0.30. Active flavonoids were acquired by the following criteria: OB ≥ 30%, DL ≥ 0.18, and BBB ≥ −0.3 [9, 10]. Active compounds of flavonoids reported publicly in recent five years were also taken into analyses from 1st January 2016 to 30th May 2021. The active flavonoids were confirmed by PubChem (https://pubchem.ncbi.nlm.nih.gov/). The structures of compounds saved in SDF format were used to predict potential targets of active molecules by using Swiss TargetPrediction platform (http://www.swisstargetprediction.ch/), and the species were set as “Homo sapiens” and the probability was set larger than 0. All the targets from UniProt database (https://www.uniprot.org/) were in standardized format.

2.2. Determination of Glioma-Related Targets and Common Targets with Active Components

Glioma-related targets were retrieved by using the GeneCards (https://www.genecards.org/) and DisGeNET (https://www.disgenet.org/home/) databases with the keyword “glioma.” Flavonoid targets and glioma-related genes were imported into Venny 2.1 (https://bioinfogp.cnb.csic.es/tools/venny/) to acquire common targets as the potential targets for further analyses.

2.3. Network Construction of Common Targets

A protein-protein interaction (PPI) network was constructed using the STRING platform (version 11.5, https://string-db.org/), and “Homo sapiens” and Medium Confidence (0.4) were set. All the information derived from the STRING database was then imported into Cytoscape software (version 3.8.2) for visual display.

2.4. GO Function and KEGG Pathway Enrichment Analyses

GO is wildly used to study gene functions, including the biological process (BP), molecular function (MF), and cell component (CC) [11]. KEGG (http://www.kegg.jp/) is an integrated database of genomic, chemical, and system functional information and is extensively used to capture significantly enriched biological pathways [12]. The common targets of flavonoid-glioma were imported into the Metascape portal (http://metascape.org/gp/index.html/) for enrichment analysis. The top 20 GO and KEGG pathway enrichment analysis results were visualized by clusterProfiler R package as the histogram or bubble graph with [13].

2.5. Molecular Docking

The most potential pathway and its related genes of flavonoids treating glioma were obtained. Then, these targets with the most promising ingredients were reconfirmed by using molecular docking. Crystal structures of related proteins were obtained from the RCSB Protein Data Bank (PDB, https://www.rcsb.org/) with high resolution and score, water was removed, while hydrogens were added by MGLTools software (version 1.5.6). Verified compounds in.mol2 format were acquired from the TCM@Taiwan database (https://tcm.cmu.edu.tw/). Molecular docking was carried out using iGEMDOCK software (version 2.1) with default parameters. We selected the most potential proteins which had the lowest energy and determined their docking ligands using the SwissDock platform (http://www.swissdock.ch/docking/). The results were visually displayed by UCSF Chimera software (version 1.15).

3. Results

3.1. Active Ingredients of Flavonoid

According to the TCMSP database, 55 compounds were screened out with the thresholds of OB ≥ 30%, DL ≥ 0.18, and BBB ≥ −0.3 (Table 2). Twenty-nine compounds were ruled out due to the aforementioned screening conditions, but they have been reported to have antiglioma properties in previous studies [2, 3, 1465]. We added them into our study to decipher the whole view of flavonoids' antiglioma molecular mechanism. As a result, a total of 84 active compounds were selected for further analysis. Quercetin, epigallocatechin-3-gallate, isoliquiritigenin, genistein, apigenin, kaempferol, and luteolin had 154, 140, 124, 97, 80, 63, and 57 targets, respectively. It revealed that these seven flavonoids probably played significant roles in curing glioma (Table 3).

3.2. Overlapping Common Targets of Flavonoid-Glioma

After exclusion of duplicated data, 5086 and 3097 glioma-related targets were identified from GeneCards and DisGeNET databases, respectively, and 569 candidate targets of active flavonoids were integrated from Swiss TargetPrediction. Two hundred and fifty-eight intersection targets were obtained among these three gene sets and were used for further analysis (Figure 2).

3.3. Common Targets Network Construction

These 258 putative gene targets correlated with glioma were analyzed using the STRING database. A total of 258 nodes and 4407 edges were embodied with the average node degree 34.2. These results were imported into Cytoscape software (version 3.8.2) for further analysis. The network is shown in Figure 3. The node color reflected the number of interacted nodes, and the more nodes to one node linked with, the deeper colored it became, as shown in Figure 4.

The potential targets were AKT1, EGFR, VEGFA, MAPK3, CASP3, SRC, HRAS, TNF, MAPK1, CCND1, ESR1, HSP90AA1, and MTOR as their degrees were above 100 (Table 4). Node degrees were counted by Cytoscape. The greater a node degree is, the more important biological functions the node has in the PPI network (Table 5).

3.4. GO and KEGG Pathway Enrichment Analyses

GO terms were enriched by the Metascape platform. The results showed that BP terms enriched in glioma-flavonoids overlapping targets mainly included peptidyl-tyrosine phosphorylation and modification, response to oxidative stress and oxygen levels. The top five enriched CC terms were membrane raft, membrane microdomain, membrane region, neuronal cell body, and transferase complex, transferring phosphorus-containing groups. MF terms displayed the intersection genes that were mainly enriched in protein tyrosine kinase activity, protein serine/threonine kinase activity, transmembrane receptor protein kinase activity, transmembrane receptor protein tyrosine kinase activity, and phosphatase binding (Figure 5).

KEGG pathway enrichment analysis of the 258 intersection gene targets was carried out by Metascape. The main pathways among these genes included PI3K-Akt, Ras, HIF-1, and Neurotrophin signaling pathways (Figure 6).

3.5. Molecular Docking

The result of the KEGG pathway enrichment analysis indicated that the PI3K-Akt signaling pathway was the main pathway through which flavonoids affected the glioma. Among the potential targets, AKT1, EGFR, MAPK1, MAPK3, CCND1, MTOR, VEGFA, HRAS, and HSP90AA1 were enriched in the PI3K-Akt signaling pathway. We selected seven potential active molecules, including quercetin, epigallocatechin-3-gallate, isoliquiritigenin, genistein, apigenin, kaempferol, and luteolin, to dock with nine target proteins, and chose temozolomide as the control. Lower binding energy indicates a stabler conformation. We used the quantitative value of fitness to evaluate the binding level. Fitness is the total energy of a predicted pose in the binding site. The empirical scoring function of iGEMDOCK is estimated as follows: Fitness = vdW + Hbond + Elec. The vdW term is van der Waal energy; Hbond and Elect terms are hydrogen bonding energy and electrostatic energy, respectively [66]. The results were visually displayed with a heatmap (Figure 7). It was interesting to note that epigallocatechin-3-gallate (EGCG) had a good bonding ability to most target proteins enriched in the PI3K-Akt signaling pathway, while the opposite pattern was observed for isoliquiritigenin. Experimental studies show that quercetin could induce autophagy and apoptosis in human neuroglioma cells through the PI3K-Akt signaling pathway [67]. AKT1, MTOR, CCND1, and EGFR are closely associated with autophagy and apoptosis in glioma [6870]. Our findings obviously showed that EGCG had a better docking score to these proteins than quercetin. For further research, quercetin and EGCG were selected in this study to dock with MAPK1 and HSP90AA1 target proteins individually due to their relatively lower energy value in the molecular docking. The results of SwissDock revealed that the estimated Gibbs free energies (ΔG) of best binding modes of EGCG with two targeting proteins were −9.27 kcal/mol and −8.53 kcal/mol, respectively, while the binding energies of quercetin with two targeting proteins were −8.23 kcal/mol and −7.95 kcal/mol, separately. In addition, EGCG had one backbone hydrogen bond (HB) interacting with Glu33 of MAPK1, and the distance was 2.031 Å. One backbone HB bounded with Gly97 of HSP90AA1 with a distance of 2.117 Å. The results displayed by UCSF Chimera software were shown in Figure 8.

4. Discussion

In recent years, flavonoids are widely used for antiglioma treatment. The mechanisms of flavonoids are very complex because they have multiple potential targets and active components. Network pharmacology together with bioinformatics has superiority in the systematic elucidation of the mechanism of TCM at the molecular level and representation of interactions between active compounds, potential targets, and various pathways.

In our study, potential targets of active components analysis revealed that quercetin, EGCG, isoliquiritigenin, genistein, apigenin, kaempferol, and luteolin interacted with multiple targets in the network. These findings showed that they may play important roles in the treatment of glioma. As reported, these seven active flavonoid compounds have distinct ways of treating glioma. Quercetin, a flavonol, had the most potential targets in this study. It plays antiglioma effects by inducing cell apoptosis [17], inhibiting proliferation and migration [14], and modulating the inflammatory process [71]. Moreover, quercetin could affect human glioma cells through the PI3K-Akt signaling pathway [72]. EGCG is a polyphenol flavonoid, which is generally distributed in green tea and has shown great properties in cancer prevention due to its safety, low cost, and excellent bioavailability [73]. EGCG in high doses (>40 μmol/L) could suppress cancer cells by inducing apoptosis and by inhibiting autophagic processes [74] and regulate apoptosis-related and autophagy-related proteins (caspase3, caspase 9, Bax, LC3B II, and Beclin) [75]. MTOR is a key regulator of autophagy, EGCG may enhance the phosphorylation of eNOS and mTOR via the activation of the PI3K-Akt pathway [76]. Furthermore, it has the effects of antiglioma through inhibiting proliferation and decreasing invasion of glioma cells [77]. Isoliquiritigenin, isolated from licorice, has been found to be a potent stimulator of cell differentiation and has potential application for treating human brain glioma by inhibiting proliferation and blocking angiogenic through Notch1 and Akt signaling pathway, respectively [22, 23]. Genistein, an isoflavone in legumes and some herbal medicines, suppresses the expression of matrix metalloproteinase 2 (MMP-2) and vascular endothelial growth factor (VEGF) to serve antigiloma role [78]. Genistein sensitizes glioblastoma cells to carbon ions through inhibiting DNA-PKcs phosphorylation and subsequently repressing the nonhomologous end-joining and delaying the homologous recombination repair pathways [24]. Apigenin, a flavone, has been shown to take part in restoring the immune system and weakening the self-renewal and invasiveness capacity of glioblastoma stem-like cells (GSCs) [25, 29]. It was reported to inhibit the expression of STAT3, AKT, and MAPK in the GSCs [26]. Kaempferol has also been demonstrated to possess good antiglioma effects by inducing reactive oxygen species (ROS) and subsequently leads to autophagy and cell death [30, 79]. Luteolin is a flavone and has an inhibitory effect on downstream signal molecules activated by EGFR, particularly the Akt and MAPK signal pathways [33, 80]. It induces a lethal endoplasmic reticulum stress response and mitochondrial dysfunction in glioblastoma cells by increasing intracellular ROS levels [31].

Immune factors have been considered as a significant factor contributing to the development and progression of glioma [81]. In the PPI network, most potential targets were closely related to immunity, including AKT1, TNF, EGFR, VEGFA, MAPK1, MAPK3, CASP3, SRC, HRAS, CCND1, ESR1, HSP90AA1, and MTOR [8287]. And these proteins were regarded as core proteins in our study and might play important roles in the therapeutic effect of flavonoids on glioma. Recent studies have shown that luteolin decreased the expression of immune-related genes including MMP9, MAPK1, HSP90AA1, CASP3, ALB, EGFR, SRC, HRAS, and ESR1. And among these genes, MMP9, MAPK1, HSP90AA1, EGFR, SRC, and HRAS are confirmed in vivo at the protein and mRNA levels [88].

To further indicate the potential mechanism of flavonoids in treating glioma, KEGG analysis discovered that PI3K-Akt was the main signaling pathway. It is a classic signal transduction pathway involved in cell proliferation, apoptosis, migration, invasion, and angiogenesis in glioma and plays an important role in the occurrence and development of glioma [89]. The result of molecular docking showed that EGCG had good bonding with MAPK1 and HSP90AA1 in the PI3K-Akt signaling pathway. Relevant studies confirmed that EGCG induces apoptosis, inhibits proliferation, and decreases invasion of glioma cells via the MAPK pathway in vivo [77]. Kim et al. also found that EGCG induced the expression of MAPK1 in glioma cells [90]. Heat Shock Protein 90 can promote oncogenesis since it interacts and supports numerous proteins and is essential for malignant transformation and progression. However, the HSP90AA1 gene is not altered in a major of tumors according to the Cancer Genome Atlas (TCGA) [91]. To evaluate its role in the treatment of glioma, downregulation of HSP90AA1-IT1 (HSP90AA1 intronic transcript (1) was done, which could significantly suppress cell viability, proliferation, EMT, invasion, and migration of glioma [92]. Thus, there might be a correlation between HSP90AA1 and glioma; however, there is no report about EGCG curing glioma via targeting HSP90AA1.

Although there is an abundance of information and the analysis process is complex, some useful and credible conclusions have been drawn. Due to limitations of compounds screening and accuracy of target prediction, the results obtained in this study are general, and in vitro and in vivo experiments are needed for verification. In short, our study portrayed the ground view of flavonoids in the treatment of glioma.

5. Conclusions

This study elaborated the mechanisms of active flavonoids on antiglioma using network pharmacology and molecular docking by constructing a compound‒target‒pathway network. Active components have particular advantages in curing glioma by targeting MAPK1, MAPK3, EGFR, MTOR, AKT1, VEGFA, CCND1, HSP90AA1, and HRAS. In addition, EGCG can target HSP90AA1 and MAPK1 via the PI3K-Akt signaling pathway. These findings offered a research foundation for further investigation of flavonoids on antiglioma.

Abbreviations

BBB:Blood-brain barrier
BP:Biological process
CC:Cell component
DL:Drug-likeness
EGCG:Epigallocatechin-3-gallate
EMT:Epithelial-mesenchymal transition
eNOS:Endothelial nitric oxide synthase
GO:Gene Ontology
GSCs:Glioblastoma stem-like cells
HB:Hydrogen bond
HIF-1:Hypoxia-inducible factor-1
KEGG:Kyoto Encyclopedia of Genes and Genomes
MF:Molecular function
MMP-2:Matrix metalloproteinase 2
OB:Oral bioavailability
PDB:Protein data bank
PI3K-Akt:Phosphatidylinositol 3-kinase-protein kinase B
PPI:Protein-protein interaction
ROS:Reactive oxygen species
SDF:Structure-data file
STAT3:Signal transducer and activator of transcription 3
TCM:Traditional Chinese medicine
TCMSP:Traditional Chinese medicine database and analysis platform
ΔG:Gibbs free energy.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

Hong Xiao and Hongyi Liu designed the study. Taiping Li and Yong Xiao participated in the study design and wrote the article. Zhen Wang conducted data analysis. All authors read and approved the final manuscript. Taiping Li and Yong Xiao contributed equally to this work.

Acknowledgments

This study was supported by the Medical Research Foundation of Jiangsu Health Commission (H2019059), Science and Technology Development Foundation of Nanjing Medical University (NMUB2018095), Medical Science and Technology Development Foundation of Nanjing (YKK19088), and Postgraduate Research & Practice Innovation Program of Jiangsu Province (SJCX19-0331 and KYCX20-1419). In addition, the authors are very grateful to Dr. Wang (Lijun Wang) for her careful revision of our manuscript.