The potential mechanism of Longsheyangquan Decoction on the treatment of bladder cancer: Systemic network pharmacology and molecular docking

Our goal was to explore the bioactive constituents of Longsheyangquan (LSYQ) Decoction and elucidate its mechanisms on the treatment of bladder cancer (BCa). A total of 38 compounds were selected based on their pharmacokinetic properties in three large traditional Chinese medicine (TCM) databases. 654 putative targets of LSYQ Decoction were predicted using a structure-based, reverse-docking algorithm online, of which 343 overlapped with BCa-related protein-coding genes. The protein-protein interaction (PPI) network was constructed to perform module analysis for further Gene Ontology (GO) annotations and Kyoto Encyclopedia Genes and Genomes (KEGG) pathway enrichment analysis, which identified CDK2, EGFR, MMP9 and PTGS2 as hub targets. The TCM-compound-target network and compound-target-pathway network together revealed that quercetin, diosmetin, enhydrin and luteolin were the main components of LSYQ Decoction. Finally, molecular docking showed the affinity between the key compounds and the hub target proteins to verify the accuracy of drug target prediction in the first place. The present study deciphered the core components and targets of LSYQ Decoction on the treatment of BCa in a comprehensive systemic pharmacological manner.

Introduction

It is estimated that there were 573,278 new cases and 212,536 disease-related deaths of bladder cancer worldwide in 2020 alone, making it the 10th most commonly diagnosed cancer globally and the most common malignancy of the urinary system (Sung et al., 2021). The past 2 decades have witnessed a transition of the cancer spectrum in China that is shifting from a less-developed country to a developed one (Feng et al., 2019). At the same time, the incidence of BCa has increased to a level similar to that of Western Europe and the United States, due to an increasingly westernized lifestyle and an upward trend of smoking among women (Feng et al., 2019). Despite being confined to the mucosa (stage Ta, carcinoma in situ) or submucosa (stage T1) at the early stage, BCa often recurs and requires frequent imaging examinations, such as urinary ultrasonography, contrast-enhanced CT, or sometimes even transurethral cystoscopy. However, approximately 30% of non-muscle invasive bladder cancers would gradually progress into muscle-invasive ones during the course of the disease (van der Heijden and Witjes, 2009). Postoperative immediate and regular intravesical administration of Bacillus Calmette-Guerin, pirarubicin, mitomycin or gemcitabine could effectively kill the residual tumor cells after transurethral resection of bladder tumors for superficial BCa (Babjuk et al., 2022). Radical cystectomy or bladder-sparing surgery followed by platinum-based chemotherapy remains the first-line treatment for advanced-stage patients. Although neoadjuvant chemotherapy and especially immunotherapy with checkpoint inhibitors have improved the prognosis of muscle-invasive and metastatic BCa in recent years (Singh and Black, 2016; Burdett et al., 2022), the 5-year overall survival and response rate to the immunotherapy remain low (Gakis, 2020). In addition to the uncertain therapeutic outcome, the local and systemic side effects associated with chemo- or immunotherapy have become a major reason for discontinued treatment and have limited their further clinical applications (Singh and Black, 2016). Therefore, it is of critical urgency to develop novel cures with high efficacy while reducing toxicity.

Traditional Chinese medicine usually consists of various herbal medicines derived from natural plants, and has been used for thousands of years to treat different types of diseases, with the goal of establishing the dynamic balance of the whole body. Longsheyangquan (LSYQ) Decoction whose composition is primarily Solanum Nigrum Linn., Duchesnea Indica Focke, Herba Solani Lyrati, Spora Lygodii, Smilacis Glabrae Rhizoma and Juncus effusus L., is an empirical prescription for preventing postoperative BCa recurrence on the basis of TCM theory (Zhou, 2006). While some of these ingredients have been proven to exert anti-cancer effect, the underlying mechanisms of LSYQ Decoction as a whole and its promising therapeutic effect on BCa are still largely unclear.

Network pharmacology, first proposed by Dr. Hopkins, is an efficient tool for systemic pharmacological research in which the interactions among target genes, pathways, diseases, and drugs are studied in detail using a network-based method to unravel the mechanisms of complex preparations at multiple levels (Hopkins, 2007). Due to the in-depth analysis and holistic view shared with classical Chinese medicine, it is quite suitable for the research of TCM herbs and formulas. Molecular docking is a frequently-used method in the field of structure-based drug design that predicts the preferred conformation when a small molecule ligand binds to a target to form a stable complex, and the strength of binding, i.e., the affinity between them, can be calculated (Lengauer and Rarey, 1996).

Herein, we aimed to comprehensively investigate LSYQ Decoction as well as its influence on BCa using an in-silico approach, as neither had been reported before. Initially, we screened out and selected the putative bioactive components of LSYQ Decoction. Then, we acquired the overlap between the drug targets and BCa-related target genes for further GO and KEGG enrichment analyses. Meanwhile, we visualized the TCM-compound-target network as well as the compound-target-pathway network to provide a general overview of the potential anti-tumor properties and molecular mechanisms of LSYQ Decoction. Finally, we employed molecular docking to confirm the affinity between the core components and the hub targets. The current study not only elucidated the possible mechanisms of LSYQ Decoction on the treatment of BCa, but also serve to promote more research into TCM in order to discover highly effective herbal medications with mild side effects and relatively low research and development investment.

A systemic scheme of the pharmacological network for the present study is illustrated in Figure 1.

FIGURE 1

FIGURE 1. The scheme for investigation of mechanisms of LSYQ Decoction on the treatment of BCa.

Methods

Chemical compounds of LSYQ Decoction and targets prediction

The Traditional Chinese Medicine System Pharmacology (TCMSP) online database (https://tcmsp-e.com/tcmsp.php) was utilized in search of chemical constituents of LSYQ Decoction, which is composed of Solanum Nigrum Linn., Duchesnea Indica Focke, Herba Solani Lyrati, Spora Lygodii, Smilacis Glabrae Rhizoma and Juncus effusus L. The database integrates TCMs, chemical components and their pharmacokinetic properties, including oral bioavailability (OB), gastrointestinal absorption, blood-brain barrier permeability, drug-likeness (DL) etc. The active compounds were selected by criteria as follows: OB ≥ 30%, DL ≥ 0.18, molecular weight ≤500 Dal, hydrogen bond donors ≤5, hydrogen bond acceptors ≤10 and AlogP ≤5, with the last four corresponding to the “Rule of Five” first raised by Lipinski in 1997 (Lipinski, 2016). We then searched HERB and the Traditional Chinese Medicine Information Database (TCM-ID) for additional chemical constituents as supplementary, the former launched by the Institute of Computing Technology of Chinese Academy of Sciences as well as Beijing University of Chinese Medicine, and the latter by the National University of Singapore. Literature mining was also conducted on the China National Knowledge Internet (CNKI, https://www.cnki.net/) for more possible constituents, especially in the case of Duchesnea Indica Focke whose components were re-input into the TCMSP to obtain pharmacokinetic properties for next-step selection. Additional compounds were input into SwissADME (http://www.swissadme.ch/) (Daina et al., 2017)for pharmacokinetic and drug-likeness parameters. Only compounds with high gastrointestinal absorption that met at least three of the five drug-likeness principles were selected.

Targets for a total of 38 potentially-bioactive compounds were predicted by SwissTargetPrediction (http://www.swisstargetprediction.ch/) (Daina et al., 2019), which determines the most likely protein targets of a small molecule based on a combination of 2D and 3D similarity. Homo Sapiens was chosen as the only species.

Construction of TCM-compound-target network

Cytoscape software (version 3.8.0) was used to vividly display the relationship between TCMs, potentially-bioactive compounds and targets associated with the compounds.

Identification of BCa-related targets

Human BCa-related targets were searched using the following terms: “bladder cancer” AND “bladder carcinoma” AND “malignant neoplasm of bladder” AND “bladder carcinoma in situ”. Targets were retrieved from GeneCards (https://www.genecards.org/), Therapeutic Target Database (TTD, http://db.idrblab.net/ttd/) and GenCLiP 3 (http://ci.smu.edu.cn/genclip3/analysis.php). For GeneCards, only targets with relevance scores greater than the median were selected. The universal Protein Resources (UniProt, https://www.uniprot.org/) was subsequently used to distinguish protein-coding target genes from non-coding ones.

Establishment of protein-protein interaction network and module analysis

Overlapping targets of LSYQ Decoction and BCa were acquired and considered as candidate targets by means of a Venn diagram via an online plotting resource (http://jvenn.toulouse.inra.fr/app/index.html). An original PPI network involving candidate targets was created by the Search Tool for the Retrieval of Interacting Genes (STRING, https://string-db.org/) in order to explore the relationship between the imported protein-coding genes. Afterwards, the original PPI network was analyzed by Cytoscape and its CytoNCA plugin. Nodes whose Betweenness, Closeness and Degree were less than the corresponding medians were filtered out. Furthermore, the Molecular Complex Detection (MCODE) plugin helped us to select and visualize significant modules with criteria as follows: degree cut off = 2, node score cut off = 0.2, K-core = 2 and maximum depth = 100. Modules with score ≥10 were included for further analyses.

Functional annotation and enrichment analyses

Gene ontology analysis is a widely-accepted approach to define and describe genes across species from three aspects, namely molecular function (MF), biological process (BP), and cellular component (CC). KEGG is in fact an enormous collection of databases involving genomes, drugs, enzymes, biological pathways etc. Herein, the Metascape (https://metascape.org/gp/index.html#/main/step1) was applied to perform GO and KEGG enrichment analyses for two clusters with adjusted p values <0.01 as the cut-off.

Construction of compound-target-pathway network

A network was plotted by Cytoscape for the top 25 KEGG pathways in the top-ranked module with the corresponding targets and chemical compounds so as to highlight the hub target genes and the core compounds.

Molecular docking analysis

Molecular docking could intuitively display the interaction between the compounds and their potential target proteins with the help of Autodock (version 4.2.6). The 3D structures of the core compounds were downloaded from PubChem (https://pubchem.ncbi.nlm.nih.gov/), while the crystal structures of human proteins encoded by the hub genes were obtained from the RCSB Protein Data Bank (PDB, https://www.rcsb.org/). The protein structures were then pre-processed by removing the irrelevant water molecules and ligands, adding hydrogen atoms, and merging non-polar hydrogen atoms. After adding hydrogen atoms for the compounds, protein-compound docking was conducted with the method of Lamarckian genetic algorithm in Autodock Tools (version 1.5.6). The number of runs for each docking was set to 50 as officially recommended, in order to gain the best affinity (i.e., the lowest binding energy), the conformation of which was then visualized in PyMOL software (version 2.2.0).

Results

Chemical compounds and potential targets of LSYQ Decoction

38 of the 386 compounds in LSYQ Decoction were identified as drug-like molecules (Tables 1, 2) and had their putative targets successfully predicted by SwissTargetPrediction, which yielded 654 human proteins after the removal of duplicates. There were six components in Solanum Nigrum Linn., four components in Duchesnea Indica Focke, three components in Herba Solani Lyrati, four components in Spora Lygodi, five components in Smilacis Glabrae Rhizoma and 19 components in Juncus effusus L. It is worth noting that quercetin and diosgenin were present in two or three herbal medicines.

TABLE 1

TABLE 1. Bioactive compounds retrieved from TCMSP.

TABLE 2

TABLE 2. Bioactive compounds retrieved from other two databases and through literature mining.

TCM-compound-target network

The intricate relationship between six TCMs, 38 compounds, and 654 compound-related targets was clearly illustrated by Cytoscape in Figure 2. TCMs, active compounds and targets were represented by green round rectangular-, orange circular- and blue diamond-shaped nodes, respectively. The greater the Degree value of a node was, the larger the size and lower the transparency were in the network. This network emphasized the essence of Chinese medicine formulas, that is, multi-components of a formula exerting multi-target effects.

FIGURE 2

FIGURE 2. The TCM-compound-target network with 698 nodes and 2,861 edges. TCMs, active compounds and targets were represented by green round rectangular-, orange circular- and blue diamond-shaped nodes, respectively. The greater the Degree value of a node was, the larger the size and lower the transparency were in the network.

Identification of candidate targets

A total of 7,078 BCa-related genes were extracted from the GeneCards, TTD and GenCLiP three databases. After filtering based on the median of relevance scores and removing duplicates, 2,987 human protein-coding target genes were finally enrolled. Integration of potential TCM target genes and BCa-associated genes presented by a Venn diagram revealed that there were 343 genes in common, indicating their potential therapeutic role in the treatment of BCa by LSYQ Decoction (Figure 3 and Supplementary Table S1).

FIGURE 3

FIGURE 3. The 343 overlapping targets of LSYQ Decoction and BCa were identified by a Venn diagram.

Construction and analysis of PPI network

An original PPI network containing 339 nodes and 6,471 edges was constructed after importing 343 candidate target genes into the STRING website while discarding four disconnected genes. The network owned an average node degree of 37.7 and an average local clustering coefficient of 0.535. The CytoNCA plugin in Cytoscape software was then employed to generate three fundamental topological measures of the centrality of a node in the network: Betweenness, Closeness and Degree. Nodes whose three aforementioned parameters were greater than their corresponding medians formed a core network with 133 nodes and 3,323 edges, in which the MCODE plugin detected two crucial modules with scores ≥10 (Figure 4). The highest scoring module (Cluster One) had a score of 40.85 with 48 nodes and 960 edges, while the other (Cluster Two) had 43 nodes and 261 edges with a score of 12.43 (Figure 6A and Supplementary Table S2).

FIGURE 4

FIGURE 4. The process of screening the hub targets. The left panel showed the initial PPI network (339 nodes and 6,471 edges) of overlapping candidate targets. The middle panel showed the median-filtered PPI network with 133 nodes and 3,323 edges. The right panel (Cluster One) was generated as a result of MCODE module analysis from the middle panel. Cluster One had 48 nodes and 960 edges with a MCODE score of 40.85.

Functional annotation and enrichment analyses

The target genes in two modules were each analyzed on the Metascape website with the log(adjusted p-value) ≤ −2 as the cut-off. For Cluster One, the BP analysis revealed that the hub genes were enriched in cellular response to organonitrogen compound, response to reactive oxygen species, positive regulation of cell migration and so on (Figure 5A). The CC analysis showed that hub genes were mainly related with cyclin-dependent protein kinase holoenzyme complex, transcription regulator complex, protein kinase complex and the like (Figure 5B). The MF analysis found that the hub genes were concentrated in protein kinase binding, protein serine/threonine/tyrosine kinase activity, transcription factor binding etc. (Figure 5C). KEGG pathway enrichment analysis indicated that hub genes were remarkably enriched in bladder cancer, endocrine resistance, proteoglycans in cancer and pathways in cancers, such as PI3K-Akt signaling pathway, ErbB signaling pathway, VEGF signaling pathway and MAPK signaling pathway (Figure 5D and Supplementary Table S3). As for Cluster Two, the genes were significantly enriched in pathways in cancer, PI3K-Akt signaling pathway, EGFR tyrosine kinase inhibitor resistance. In addition, a total of 457, eight and 25 terms were enriched in BP, CC and MF, respectively (Figure 6B, full data not shown).

FIGURE 5

FIGURE 5. Enrichment analyses of Cluster One for the top 20 GO annotations and top 25 KEGG pathways (A) GO analysis of biological process terms (B) GO analysis of cellular component terms (C) GO analysis of molecular function terms (D) KEGG pathway analysis.

FIGURE 6

FIGURE 6. Composition and enrichment analyses of Cluster Two (A) Cluster Two had 43 nodes and 261 edges with a MCODE score of 12.43 (B) Significantly enriched GO BP/CC/MF terms and KEGG pathways of Cluster Two were displayed in the Nightingale rose diagram, with the number above each petal representing the gene count enriched in this term. The larger the petal, the lower the adjusted p-value (all the adjusted p-value were under 0.01).

Meanwhile, a network concerning the top 25 KEGG pathways in Cluster One was drawn along with the corresponding targets and drug compounds to accurately illustrate the underlying regulatory mechanisms of LSYQ Decoction for BCa (Figure 7). According to the two PPI networks, the Degree values of Cyclin-dependent kinase 2 (CDK2), Epidermal growth factor receptor (EGFR), Matrix metalloproteinase 9 (MMP9) and Prostaglandin-endoperoxide synthase 2 (PTGS2) were among the highest, suggesting their roles as hub targets.

FIGURE 7

FIGURE 7. The compound-target-pathway network for the top 25 KEGG pathways in Cluster One. Pathways, compounds and targets were represented by green triangular-, orange circular- and blue diamond-shaped nodes, respectively. The greater the Degree value of a node was, the larger the size and lower the transparency were in the network.

Molecular docking

To verify the accuracy of SwissTargetPrediction and to ascertain whether key molecules in LSYQ Decoction could directly bind to their protein targets, a widely-accepted molecular docking approach was performed via Autodock. Four typical critical compounds (quercetin, diosmetin, enhydrin and luteolin) attempted to dock with four hub target proteins (CDK2, EGFR, MMP9 and PTGS2). The lower the binding energy, the stronger the bond. The lowest binding energy, in other words, the best affinity, of 16 dockings were shown in a heatmap (Figure 8A). Binding energies lower than −5 kcal/mol, especially < −7 kcal/mol, are considered to form more stable structures than those with higher binding energies. The docking results suggested that all molecules except enhydrin could more or less form stable conformations with four protein receptors, which is consistent with the SwissTargetPrediction results (Figure 8B). Since both quercetin and diosmetin are flavonoids, the sites of the amino acid residues forming hydrogen bonds with either of them in CDK2 were completely identical. The same happened to the arginine residue at the 424th location of MMP9 when docking with quercetin or diosmetin.

FIGURE 8

FIGURE 8. Molecular docking of four typical core compounds with four hub target proteins (A) The binding energy heatmap with a deeper color suggesting a more stable binding (B) The molecular docking between CDK2, MMP9 and quercetin or diosmetin. The yellow and orange ligands represented quercetin and diosmetin, respectively. The enlarged view on the right showed the interactions of ligands and amino acid residues of CDK2 and MMP9.

Discussion

Bladder cancer is the most common malignancy of the urinary system with painless gross hematuria as the cardinal symptom; however, patients may be asymptomatic in the early stage, leading to delayed diagnosis and eventual progression of the disease (van der Heijden and Witjes, 2009). While early detection depends primarily on imaging technologies, some patients are diagnosed as advanced BCa in the first place (Gakis, 2020). Furthermore, men are generally thought to be at higher risk for BCa than women. This gender pattern of incidence persists worldwide, but divided trends have been observed, particularly in developed countries, with consistent rates in men and increasing tendencies among women (Sung et al., 2021). The rising prevalence of smoking in women may partly explain this. The focus of management for BCa should hence be on the prevention of relapse, progression into the advanced stage, and metastasis, as non-muscle invasive BCa is always associated with a promising prognosis, whereas muscle invasive BCa remains refractory, regardless of whether the bladder has been resected or not (van der Heijden and Witjes, 2009). Although great progress has been made in treatment, for instance the emergence of brand-new adjuvant immunotherapy, only a small number of patients benefit from it, not to mention the toxic side effects it brings about (Rey-Cárdenas et al., 2021). With its precious clinical treatment experience and abundant theoretical resources dating back to ancient China, TCM could be a perfect complementary remedy and a key entry point for the development of natural medicines with multi-target potential.

From the perspective of TCM, the bladder is an organ responsible for eliminating the “damp evil” that arises daily in the human body. Should the bladder fail to eliminate it, it would accumulate and transform into a combination of damp and heat evils, eventually causing the cancer together, during which time the two evils could probably block neighboring vessels. This is how hematuria comes into being, according to TCM theory. The botanical medicine in LSYQ Decoction serve precisely the TCM principle of clearing the heat and prompting diuresis. Solanum Nigrum Linn., a herbal plant rich in flavonoids, is widely adopted in TCM for its antineoplastic effect (Butt et al., 2018). Previous studies shown that alpha-solanine, an extract of Solanum Nigrum Linn., could inhibit invasion of human prostate cancer cells by suppressing epithelial-mesenchymal transition and matrix metalloproteinase (MMP) expression (Shen et al., 2014), while Solanum Nigrum Linn. Polysaccharide-1, another active ingredient therein, exhibited immunomodulatory activity by promoting the release of NO, tumor necrosis factor-α and interleukin (IL)-6 in macrophages and inducing Th1 responses mediated by IL -2 and interferon-γ (Pu et al., 2020). Extracts of Duchesnea Indica Focke, also known as mock strawberry, have been reported to have anti-metastatic potential by prohibiting MMP2 activity in human oral squamous cells (Yang et al., 2019). As one of the most important elements in Baiying Qinghou Decoction, Herba Solani Lyrati (Baiying in Chinese) could target TP53, EGFR, IL1B and NOS3 proteins for the treatment of laryngeal squamous cell carcinoma (Gao et al., 2021). Smilacis Glabrae Rhizoma, a member of the Smilacaceae family native to Southeast Asia, has long been utilized to treat sore throat and syphilis. It also exhibited an anti-tumor effect in a dose-dependent manner in human breast cancer cell line MCF7, colon carcinoma cell line HT -29, and gastric cancer cell line BGC-823 by inducing mitochondrial apoptosis (Gao et al., 2011). Juncus effusus L. has been attributed with a variety of pharmacological functions, including anti-inflammatory, anti-allergic, anti-cancer and others (Ma et al., 2016; Park et al., 2016). For example, dehydroeffusol, a low molecular weight phenanthrene from Juncus effusus L., effectively induced endoplasmic reticulum stress, resulting in decreased proliferation of gastric cancer cells (Zhang et al., 2016). Nevertheless, the effect of LSYQ Decoction on bladder cancer has not been clearly elucidated.

Considering the composition diversity of this formula, we made use of network pharmacology to systemically explore it at herbal, molecular, and genetic levels in the present study. Upon strict screening, 38 active compounds were finally identified, among which 654 human proteins were predicted as potential targets. It again proved that TCM formulas possess the characteristics of multi-drug and multi-target. The TCM targets are intersected with BCa-related protein-coding genes via a Venn diagram. The hub target genes, namely CDK2, EGFR, MMP9 and PTGS2, were then determined through the visualization of the PPI network and module analysis. Meanwhile, quercetin, diosmetin, enhydrin, as well as luteolin were regarded as core constituents in LSYQ Decoction through topological mining. The target genes in both of the two modules, whose cluster scores were greater than 10, have accumulated in significant pathways related to cancers, such as the PI3K-Akt signaling pathway, the ErbB signaling pathway, the VEGF signaling pathway, the MAPK signaling pathway, and the HIF-1 signaling pathway.

The anti-tumor activity of flavonoids has been extensively described and widely accepted (Teillet et al., 2008). Quercetin, a flavonoid found in abundance in herbal plants and daily consumed foods, can not only restrain the proliferation and induce the apoptosis of various tumor cells via different pathways (Su et al., 2016; Mirazimi et al., 2022; Reyes-Avendaño et al., 2022), but can also act as a free-radical scavenger to exert antioxidant and anti-inflammatory effects, therefore, inhibiting oncogenesis (Hsieh et al., 2022)and inflammatory diseases (Yang et al., 2022). Citrus fruits and extracts of many herbs are also rich in diosmetin (Roowi and Crozier, 2011), another commonly encountered flavone with numerous pharmacological functions, such as anti-inflammation, anti-oxidative stress, and anti-neoplasia (Roma et al., 2018; Chen et al., 2020; Ma and Zhang, 2020). Luteolin, also a natural flavonoid, has been proven to play an important role in anti-inflammatory, antioxidant, antiallergic, and anticancer processes (Juszczak et al., 2022). It could suppress the growth of human bladder cancer cell line T24 by upregulating p21 and inhibiting mTOR signaling (Iida et al., 2020). According to the two PPI networks, the degrees of CDK2, EGFR, MMP9 and PTGS2 were among the highest, suggesting their roles as hub targets. EGFR overexpression is found in various human malignancies (Levantini et al., 2022), including BCa (Seidl, 2020), and a few monoclonal antibody drugs have been invented to target its protein in breast cancer and non-small cell lung cancer, either approved by the Food and Drug Administration of the US or under clinical trials. CDK2 is a member of the cyclin-dependent kinase family that combines with cyclin E and A in turn during the cell cycle to prompt the transition from the G1 phase to the S phase (Hume et al., 2020), making it a controller that regulates mitosis and cell proliferation. MMP9, a protein from the calcium-dependent zinc-containing endopeptidase family, is engaged in the degradation of extracellular matrix, especially type-IV collagen, resulting in multiple physiological and pathological processes, including tissue remodeling, angiogenesis, cell migration, rheumatoid arthritis, and tumor invasion (Dufour et al., 2010; Xue et al., 2014; Li et al., 2021). Higher MMP9 expression is associated with a poorer prognosis for patients with breast cancer or colorectal cancer (Dufour et al., 2011; Pezeshkian et al., 2021). PTGS2, also known as cyclooxygenase-2, converts arachidonic acid to prostaglandin endoperoxide H2 and engages in downstream inflammatory responses. Not only is the expression of PTGS2 remarkably upregulated during inflammation, but also in various tumors such as gallbladder, and breast cancers (Legan, 2010; Glover et al., 2011). What’s more, prostaglandin E2, the product of PTGS2, was reported to stimulate tumor progression (Menter et al., 2010). In addition to the four targets mentioned above, there are many other proteins participating in the occurrence and progression of BCa that also scored high in Degree values in the two PPI networks, such as AKT1, SRC, CDK1, ESR1, CCNB1, and so on. Ultimately, molecular docking was employed to examine the binding affinity of typical core compounds and hub targets. With the exception of enhydrin, the results suggested that all of the typical compounds had a relatively good binding capacity to four target proteins (affinity less than -5 kcal/mol), and the four molecules even shared stronger bonds with CDK2 and MMP9 with all the binding energies less than -7 kcal/mol (Alam MS. et al., 2022). Molecular docking analysis has increasingly become an important in-silico tool to preliminarily check the structural fitness between target proteins and small molecules regardless whether the latter were newly synthesized or repurposed for a different disease (Alam M. S. et al., 2022; Majid et al., 2022).

Our study still leaves some room for improvement. First, both the components of LSYQ Decoction and BCa-related genes were retrieved from databases. We tried our best to search for all of them and set a rational threshold for filtration. Second, the targets of LSYQ Decoction were predicted by a reverse-docking algorithm at SwissPredictionTarget. We cross-validated them by molecular docking. Last but not least, the mechanisms discovered by systemic network pharmacology are still not complete and need further experimental exploration and verification.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Author contributions

ZC involved in design of the research, literature mining, data collection and drafting of the manuscript. Others contributed to the data collection. HJ revised the manuscript and supervised the whole research. All the authors have read and approved the final manuscript.

Funding

This work was sponsored by Shanghai ShenKang Hospital Development Center (SHDC 12021104).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

The first author pays great tribute to all that care and support him along the way.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2022.932039/full#supplementary-material

Abbreviations

BCa, Bladder cancer; BP Biological Process; CC, Cellular Component; DL, Drug-likeness; GO, Gene Ontology; HIF, Hypoxia-inducible factor; IL, Interleukin; KEGG, Kyoto Encyclopedia Genes and Genomes; LSYQ, Longsheyangquan; MF, Molecular Function; mTOR, Mechanistic target of rapamycin kinase; OB, Oral bioavailability; PPI, Protein-protein interaction; TCM, Traditional Chinese Medicine; VEGF, Vascular endothelial growth factor; CDK1, Cyclin-dependent kinase one; CDK2, Cyclin-dependent kinase two; EGFR, Epidermal growth factor receptor; MMP9, Matrix metalloproteinase nine; PTGS2, Prostaglandin-endoperoxide synthase two; ESR, Estrogen receptor one.; CCNB1, Cyclin B1.

References

Alam, M. S., Rahaman, M. M., Sultana, A., Wang, G., and Mollah, M. N. H. (2022a). Statistics and network-based approaches to identify molecular mechanisms that drive the progression of breast cancer. Comput. Biol. Med. 145, 105508. doi:10.1016/j.compbiomed.2022.105508

PubMed Abstract | CrossRef Full Text | Google Scholar

Alam, M. S., Sultana, A., Reza, M. S., Amanullah, M., Kabir, S. R., and Mollah, M. N. H. (2022b). Integrated bioinformatics and statistical approaches to explore molecular biomarkers for breast cancer diagnosis, prognosis and therapies. PLOS ONE 17, e0268967. doi:10.1371/journal.pone.0268967

PubMed Abstract | CrossRef Full Text | Google Scholar

Babjuk, M., Burger, M., Capoun, O., Cohen, D., Compérat, E. M., Dominguez Escrig, J. L., et al. (2022). European association of urology guidelines on non-muscle-invasive bladder cancer (Ta, T1, and carcinoma in situ). Eur. Urol. 81, 75–94. doi:10.1016/j.eururo.2021.08.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Burdett, S., Fisher, D. J., Vale, C. L., Sternberg, C. N., Clarke, N. W., Parmar, M. K. B., et al. (2022). Adjuvant chemotherapy for muscle-invasive bladder cancer: A systematic review and meta-analysis of individual participant data from randomised controlled trials. Eur. Urol. 81, 50–61. doi:10.1016/j.eururo.2021.09.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Butt, G., Romero, M. A., Tahir, F., and Farooqi, A. A. (2018). Emerging themes of regulation of oncogenic proteins by Solanum nigrum and its bioactive molecules in different cancers. J. Cell Biochem. 119, 9640–9644. doi:10.1002/jcb.27258

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y., Wang, Y., Liu, M., Zhou, B., and Yang, G. (2020). Diosmetin exhibits anti-proliferative and anti-inflammatory effects on TNF-α-stimulated human rheumatoid arthritis fibroblast-like synoviocytes through regulating the Akt and NF-κB signaling pathways. Phytother. Res. 34, 1310–1319. doi:10.1002/ptr.6596

PubMed Abstract | CrossRef Full Text | Google Scholar

Daina, A., Michielin, O., and Zoete, V. (2017). SwissADME: A free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules. Sci. Rep. 7, 42717. doi:10.1038/srep42717

PubMed Abstract | CrossRef Full Text | Google Scholar

Daina, A., Michielin, O., and Zoete, V. (2019). SwissTargetPrediction: Updated data and new features for efficient prediction of protein targets of small molecules. Nucleic Acids Res. 47, W357. doi:10.1093/nar/gkz382

PubMed Abstract | CrossRef Full Text | Google Scholar

Dufour, A., Sampson, N. S., Li, J., Kuscu, C., Rizzo, R. C., DeLeon, J. L., et al. (2011). Small-molecule anticancer compounds selectively target the hemopexin domain of matrix metalloproteinase-9. Cancer Res. 71, 4977–4988. doi:10.1158/0008-5472.CAN-10-4552

PubMed Abstract | CrossRef Full Text | Google Scholar

Dufour, A., Zucker, S., Sampson, N. S., Kuscu, C., and Cao, J. (2010). Role of matrix metalloproteinase-9 dimers in cell migration: Design of inhibitory peptides. J. Biol. Chem. 285, 35944–35956. doi:10.1074/jbc.M109.091769

PubMed Abstract | CrossRef Full Text | Google Scholar

Feng, R. M., Zong, Y. N., Cao, S. M., and Xu, R. H. (2019). Current cancer situation in China: Good or bad news from the 2018 global cancer statistics? Cancer Commun. (Lond) 39, 22. doi:10.1186/s40880-019-0368-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Gakis, G. (2020). Management of muscle-invasive bladder cancer in the 2020s: Challenges and perspectives. Eur. Urol. Focus 6, 632–638. doi:10.1016/j.euf.2020.01.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, K., Zhu, Y., Wang, H., Gong, X., Yue, Z., Lv, A., et al. (2021). Network pharmacology reveals the potential mechanism of Baiying Qinghou decoction in treating laryngeal squamous cell carcinoma. Aging (Albany NY) 13 (24), 26003–26021. doi:10.18632/aging.203786

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, Y., Su, Y., Qu, L., Xu, S., Meng, L., Cai, S. Q., et al. (2011). Mitochondrial apoptosis contributes to the anti-cancer effect of Smilax glabra Roxb. Toxicol. Lett. 207, 112–120. doi:10.1016/j.toxlet.2011.08.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Glover, J. A., Hughes, C. M., Cantwell, M. M., and Murray, L. J. (2011). A systematic review to establish the frequency of cyclooxygenase-2 expression in normal breast epithelium, ductal carcinoma in situ, microinvasive carcinoma of the breast and invasive breast cancer. Br. J. Cancer 105, 13–17. doi:10.1038/bjc.2011.204

PubMed Abstract | CrossRef Full Text | Google Scholar

Hopkins, A. L. (2007). Network pharmacology. Nat. Biotechnol. 25, 1110–1111. doi:10.1038/nbt1007-1110

PubMed Abstract | CrossRef Full Text | Google Scholar

Hsieh, H. L., Yu, M. C., Cheng, L. C., Chu, M. Y., Huang, T. H., Yeh, T. S., et al. (2022). Quercetin exerts anti-inflammatory effects via inhibiting tumor necrosis factor-α-induced matrix metalloproteinase-9 expression in normal human gastric epithelial cells. World J. Gastroenterol. 28, 1139–1158. doi:10.3748/wjg.v28.i11.1139

PubMed Abstract | CrossRef Full Text | Google Scholar

Hume, S., Dianov, G. L., and Ramadan, K. (2020). A unified model for the G1/S cell cycle transition. Nucleic Acids Res. 48, 12483–12501. doi:10.1093/nar/gkaa1002

PubMed Abstract | CrossRef Full Text | Google Scholar

Iida, K., Naiki, T., Naiki-Ito, A., Suzuki, S., Kato, H., Nozaki, S., et al. (2020). Luteolin suppresses bladder cancer growth via regulation of mechanistic target of rapamycin pathway. Cancer Sci. 111, 1165–1179. doi:10.1111/cas.14334

PubMed Abstract | CrossRef Full Text | Google Scholar

Juszczak, A. M., Wöelfle, U., Končić, M. Z., and Tomczyk, M. (2022). Skin cancer, including related pathways and therapy and the role of luteolin derivatives as potential therapeutics. Med. Res. Rev. 42, 1423–1462. doi:10.1002/med.21880

PubMed Abstract | CrossRef Full Text | Google Scholar

Legan, M. (2010). Cyclooxygenase-2, p53 and glucose transporter-1 as predictors of malignancy in the development of gallbladder carcinomas. Bosn. J. Basic Med. Sci. 10, 192–196. doi:10.17305/bjbms.2010.2684

PubMed Abstract | CrossRef Full Text | Google Scholar

Lengauer, T., and Rarey, M. (1996). Computational methods for biomolecular docking. Curr. Opin. Struct. Biol. 6, 402–406. doi:10.1016/S0959-440X(96)80061-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Levantini, E., Maroni, G., del Re, M., and Tenen, D. G. (2022). EGFR signaling pathway as therapeutic target in human cancers. Seminars Cancer Biol. 1, 1. doi:10.1016/j.semcancer.2022.04.002

CrossRef Full Text | Google Scholar

Li, T., Tan, K. S., Tu, Y. Y., Zhao, L., Liu, J., Ong, H. H., et al. (2021). Overexpression of neutrophil MMP-9 and HIF-1α may contribute to the finger-like projections formation and histo-pathogenesis in nasal inverted papilloma. J. Inflamm. Res. 14, 2979–2991. doi:10.2147/JIR.S312605

PubMed Abstract | CrossRef Full Text | Google Scholar

Lipinski, C. A. (2016). Rule of five in 2015 and beyond: Target and ligand structural limitations, ligand chemistry structure and drug discovery project decisions. Adv. Drug Deliv. Rev. 101, 34–41. doi:10.1016/j.addr.2016.04.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, A., and Zhang, R. (2020). Diosmetin inhibits cell proliferation, induces cell apoptosis and cell cycle arrest in liver cancer. Cancer Manag. Res. 12, 3537–3546. doi:10.2147/CMAR.S240064

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, W., Zhang, Y., Ding, Y. Y., Liu, F., and Li, N. (2016). Cytotoxic and anti-inflammatory activities of phenanthrenes from the medullae of Juncus effusus L. Arch. Pharm. Res. 39, 154–160. doi:10.1007/s12272-015-0680-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Majid, A., Lakshmikanth, M., Lokanath, N. K., and Poornima Priyadarshini, C. G. (2022). Generation, characterization and molecular binding mechanism of novel dipeptidyl peptidase-4 inhibitory peptides from sorghum bicolor seed protein. Food Chem. 369, 130888. doi:10.1016/j.foodchem.2021.130888

PubMed Abstract | CrossRef Full Text | Google Scholar

Menter, D. G., Schilsky, R. L., and DuBois, R. N. (2010). Cyclooxygenase-2 and cancer treatment: Understanding the risk should be worth the reward. Clin. Cancer Res. 16, 1384–1390. doi:10.1158/1078-0432.CCR-09-0788

PubMed Abstract | CrossRef Full Text | Google Scholar

Mirazimi, S. M. A., Dashti, F., Tobeiha, M., Shahini, A., Jafari, R., Khoddami, M., et al. (2022). Application of quercetin in the treatment of gastrointestinal cancers. Front. Pharmacol. 13, 860209. doi:10.3389/fphar.2022.860209

PubMed Abstract | CrossRef Full Text | Google Scholar

Park, N. Y., Kim, S. G., Park, H. H., Jeong, K. T., Lee, Y. J., and Lee, E. (2016). Anti-inflammatory effects of Juncus effusus extract (JEE) on LPS-stimulated RAW 264.7 cells and edema models. Pharm. Biol. 54, 243–250. doi:10.3109/13880209.2015.1029053

PubMed Abstract | CrossRef Full Text | Google Scholar

Pezeshkian, Z., Nobili, S., Peyravian, N., Shojaee, B., Nazari, H., Soleimani, H., et al. (2021). Insights into the role of matrix metalloproteinases in precancerous conditions and in colorectal cancer. Cancers (Basel) 13, 6226. doi:10.3390/cancers13246226

PubMed Abstract | CrossRef Full Text | Google Scholar

Pu, Y., Liu, Z., Zhong, C., Zhang, X., and Bao, Y. (2020). Immunomodulatory effects of a polysaccharide from Solanum nigrum Linne through TLR4-MyD88 signaling pathway. Int. Immunopharmacol. 88, 106973. doi:10.1016/j.intimp.2020.106973

PubMed Abstract | CrossRef Full Text | Google Scholar

Rey-Cárdenas, M., Guerrero-Ramos, F., Gómez de Liaño Lista, A., Carretero-González, A., Bote, H., Herrera-Juárez, M., et al. (2021). Recent advances in neoadjuvant immunotherapy for urothelial bladder cancer: What to expect in the near future. Cancer Treat. Rev. 93, 102142. doi:10.1016/j.ctrv.2020.102142

PubMed Abstract | CrossRef Full Text | Google Scholar

Reyes-Avendaño, I., Reyes-Jiménez, E., González-García, K., Pérez-Figueroa, D. C., Baltiérrez-Hoyos, R., Tapia-Pastrana, G., et al. (2022). Quercetin regulates key components of the cellular microenvironment during early hepatocarcinogenesis. Antioxidants 11, 358. doi:10.3390/antiox11020358

PubMed Abstract | CrossRef Full Text | Google Scholar

Roma, A., Rota, S. G., and Spagnuolo, P. A. (2018). Diosmetin induces apoptosis of acute myeloid leukemia cells. Mol. Pharm. 15, 1353–1360. doi:10.1021/acs.molpharmaceut.7b01151

PubMed Abstract | CrossRef Full Text | Google Scholar

Roowi, S., and Crozier, A. (2011). Flavonoids in tropical Citrus species. J. Agric. Food Chem. 59, 12217–12225. doi:10.1021/jf203022f

PubMed Abstract | CrossRef Full Text | Google Scholar

Seidl, C. (2020). Targets for therapy of bladder cancer. Semin. Nucl. Med. 50, 162–170. doi:10.1053/j.semnuclmed.2020.02.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, K. H., Liao, A. C., Hung, J. H., Lee, W. J., Hu, K. C., Lin, P. T., et al. (2014). α-Solanine inhibits invasion of human prostate cancer cell by suppressing epithelial-mesenchymal transition and MMPs expression. Molecules 19, 11896–11914. doi:10.3390/molecules190811896

PubMed Abstract | CrossRef Full Text | Google Scholar

Singh, P., and Black, P. (2016). Emerging role of checkpoint inhibition in localized bladder cancer. Urol. Oncol. 34, 548–555. doi:10.1016/j.urolonc.2016.09.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Su, Q., Peng, M., Zhang, Y., Xu, W., Darko, K. O., Tao, T., et al. (2016). Quercetin induces bladder cancer cells apoptosis by activation of AMPK signaling pathway. Am. J. Cancer Res. 6 (2), 498–508.

PubMed Abstract | Google Scholar

Sung, H., Ferlay, J., Siegel, R. L., Laversanne, M., Soerjomataram, I., Jemal, A., et al. (2021). Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA A Cancer J. Clin. 71, 209–249. doi:10.3322/caac.21660

CrossRef Full Text | Google Scholar

Teillet, F., Boumendjel, A., Boutonnat, J., and Ronot, X. (2008). Flavonoids as RTK inhibitors and potential anticancer agents. Med. Res. Rev. 28, 715–745. doi:10.1002/med.20122

PubMed Abstract | CrossRef Full Text | Google Scholar

van der Heijden, A. G., and Witjes, J. A. (2009). Recurrence, progression, and follow-up in non-muscle-invasive bladder cancer. Eur. Urol. Suppl. 8, 556–562. doi:10.1016/j.eursup.2009.06.010

CrossRef Full Text | Google Scholar

Xue, M., McKelvey, K., Shen, K., Minhas, N., March, L., Park, S. Y., et al. (2014). Endogenous MMP-9 and not MMP-2 promotes rheumatoid synovial fibroblast survival, inflammation and cartilage degradation. Rheumatol. Oxf. 53, 2270–2279. doi:10.1093/rheumatology/keu254

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, W. E., Ho, Y. C., Tang, C. M., Hsieh, Y. S., Chen, P. N., Lai, C. T., et al. (2019). Duchesnea indica extract attenuates oral cancer cells metastatic potential through the inhibition of the matrix metalloproteinase-2 activity by down-regulating the MEK/ERK pathway. Phytomedicine 63, 152960. doi:10.1016/j.phymed.2019.152960

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Y., Shi, G. n., Wu, X., Xu, M., Chen, C. j., Zhou, Y., et al. (2022). Quercetin impedes Th17 cell differentiation to mitigate arthritis involving pparγ‐driven transactivation of SOCS3 and redistribution corepressor SMRT from PPARγ to STAT3. Mol. Nutr. Food Res. 66, 2100826. doi:10.1002/mnfr.202100826

CrossRef Full Text | Google Scholar

Zhang, B., Han, H., Fu, S., Yang, P., Gu, Z., Zhou, Q., et al. (2016). Dehydroeffusol inhibits gastric cancer cell growth and tumorigenicity by selectively inducing tumor-suppressive endoplasmic reticulum stress and a moderate apoptosis. Biochem. Pharmacol. 104, 8–18. doi:10.1016/j.bcp.2016.01.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, Y. Q. (2006). Practical oncology of traditional Chinese medicine. Beijing: Traditional Chinese Medicine Classics Publishing House.

Google Scholar