key: cord-0030092-e249diao
authors: Ding, Pengpeng; Guo, Yutong; Wang, Canghai; Chen, Jianhong; Guo, Chunmei; Liu, Hong; Shi, Qi
title: A Network Pharmacology Approach for Uncovering the Antitumor Effects and Potential Mechanisms of the Sijunzi Decoction for the Treatment of Gastric Cancer
date: 2022-04-12
journal: Evid Based Complement Alternat Med
DOI: 10.1155/2022/9364313
sha: 32b3c82d1b1bc813dd94e13834d939e8152f2601
doc_id: 30092
cord_uid: e249diao

BACKGROUND: Sijunzi decoction (SJZD), a classic Chinese formula, has been clinically used for the treatment of gastrointestinal disorders. However, few studies have uncovered its antitumor effects and its potential mechanisms against gastric cancer (GC). Therefore, this work aimed to identify the active compounds and putative targets of the SJZD and to further explore the potential mechanisms involved in the treatment of GC. MATERIALS AND METHODS: The active compounds and potential targets of the SJZD and related genes for GC treatment were collected from a public database. Traditional Chinese medicine (TCM)-compound-target-disease networks, Venn diagrams, protein–protein interactions (PPIs), gene ontology (GO), and the Kyoto Encyclopedia of Genes and Genomes (KEGG) were used to obtain the bioactive compounds, key targets, and potential pathways. Next, the human gastric adenocarcinoma cell line NUGC-4 was inoculated subcutaneously into the right flank of NCG mice to build a tumor-bearing mouse model to further verify the findings. RESULTS: There were 117 compounds in the SJZD in total. The SJZD and GC had 161 and 3288 potential targets, respectively, among which 123 targets overlapped. The network analysis showed that quercetin, kaempferol formononetin, ginsenoside, atractylenolide III, etc., were bioactive molecules. The tumor necrosis factor (TNF), interleukin-6 (IL-6), cellular tumor antigen p53 (TP53), transcription factor AP-1 (JUN), and vascular endothelial growth factor A (VEGFA) were potential targets. A KEGG pathway enrichment analysis revealed 110 pathways involved in the pathways for cancer, including the PI3K-AKT signaling pathway. Validation experiments showed that the SJZD inhibited tumor growth and induced apoptosis in tumor cells. In addition, the SJZD downregulated expressions of VEGFA, iNOS, COX-2, and Bax/Bcl2 and inhibited the expressions of p-PI3K and p-AKT. CONCLUSION: The SJZD treats GC by inhibiting blood vessel hyperplasia and inducing cell apoptosis by regulating the PI3K/AKT pathway.

With its heterogenous nature, gastric cancer (GC) is a highly aggressive malignancy, and it is the fourth most common cancer in the world, with more than 70% of the cases occurring in the developing world [1, 2] . Due to increases in the standards of hygiene, in food conservation, and in the intake of fresh fruits, incidences of GC have decreased in most parts of the world [3, 4] . Smoking, drinking, Helicobacter pylori infections, etc., can still significantly increase the risk of developing GC [5, 6] . GC patients mainly receive treatments that include surgical resection, chemotherapy, and targeted therapy [7] [8] [9] . However, these therapy approaches impose physical, psychological, and financial burdens on GC patients [10] . Traditional Chinese medicine (TCM), which is low-cost and has low toxicity, has been utilized in China for a long time to prevent and treat GC and warrants further investigations [11, 12] . e Sijunzi decoction (SJZD) is a classic Chinese formula for "reinforcing asthenia qi" and nourishing the spleen; it consists of Renshen (Codonopsis pilosula), Baizhu (Atractylodes macrocephala), Fulin (Poria cocos), and Gancao (Glycyrrhiza uralensis) [13] . In clinical applications, researchers have recently observed that the SJZD has some therapeutic effects in the treatment of malignant digestive tumors [14] . Zhao Aiguang et al. also demonstrated that SJZD inhibited GC cell growth in vivo [15] . However, due to the multiple pathways by which TCM acts, the mechanism of the SJZD action for the treatment of GC remains to be explored.

Network pharmacology is an integrated discipline that combines systems biology, multidirectional pharmacology, computational biology, and network analysis to investigate the relationships between drugs and diseases to find drug targets [16] . Chinese herbs usually possess multiple components and targets, which are consistent with the characteristics of network pharmacology [17] . In recent years, substantial advances have been made in exploring the mechanisms of TCM action and searching for the active compounds and therapeutic targets involved by utilizing network pharmacology [18, 19] . erefore, in this study, we used network pharmacology to identify the active compounds and putative targets of SJZD and further explored the potential mechanisms involved in the treatment of GC. Furthermore, we validated the results of network pharmacology in vivo by building a GC murine model ( Figure 1 ).

e keywords "Renshen" (Renshen is the root of the acanthaceae plant), "Baizhu" (Baizhu is the rhizome of Atractylodes macrocephala Koidz.), "Fulin" (Fulin is the dried sclerotia of Poria cocos (Schw.)Wolf ), and "Gancao" (Gancao is the dried root and rhizome of Glycyrrhiza uralensis Fisch of legumes) were selected, and the herbs' compounds and targets were retrieved from the Traditional Chinese Medicine System Pharmacology Database (TCMSP, https://tcmspw.com/ index.php). e compounds were screened according to the pharmacokinetic parameters "oral bioavailability (OB) ≥ 30%" and "drug likeness (DL) ≥ 0.18". en, the targets of GC were obtained from the GeneCards (https://www. genecards.org/), Online Mendelian Inheritance in Man (OMIM, https://omim.org/), and erapeutic Target Database (TTD, https://db.idrblab.org/ttd/) databases with "Gastric Cancer" as the keyword. Next, the corresponding gene names were found in the UniProt Knowledgebase (UniProtKB) (https://www.UniProt.org/) by searching for "Homo sapiens" [20] .

e intersecting targets of the SJZD and GC were searched to make a Venn diagram. e above targets and corresponding compounds as well as the Chinese herbs were imported into a Cytoscape 3.6.1 software, which is a useful tool for analyzing and visualizing biological networks, to generate a TCM-compound-target-disease network.

e intersection targets of the SJZD and GC were uploaded to the STRING database (https://cn.string-db.org/cgi/) with the following limiting condition: species (Homo sapiens). en, the PPI network was analyzed by using Cytoscape 3.6.1, which obtained "Nodes" representing the target protein molecules and "Edges" representing the mutual relationships.

A gene ontology (GO) analysis including biological process (BP), cellular component (CC), and molecular function (MF) was performed in the DAVID database (https://david.ncifcrf. gov/home.jsp) by inputting the above targets. en, the top 10 terms were selected to make graphical structures in GraphPad Prism 8.0. e results of the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were obtained from the DAVID database, and a bubble chart was made on imageGP (http://www.ehbio.com) to explore the correlated functional pathways of the selected genes.

e human gastric adenocarcinoma cell line NUGC-4 was maintained in a RPMI-1640 medium supplemented with 10% fetal bovine serum at 37°C in a humidified atmosphere containing 5% CO 2 . e NUGC-4 cell suspension for ultimate collection had a total of 1 × 108 NUGC-4 cells in 1 mL of PBS [21] . e cells were inoculated subcutaneously into the right flanks of NCG (NOD/ShiLtJGpt-Prkdcem26Il2rgem26/Gpt, male, 5 weeks of age) mice (100 μL/mouse). e study was approved by the Medical Research Ethics Review Committee of the Beijing Shijitan Hospital of the Capital Medical University (Beijing, China). When the tumor size reached 100 mm 3 , the tumor-bearing mice were randomly divided into 4 groups: model group (n � 6), SJZD-L group (SJZD-low dose, 10 mg/ kg, n � 6), SJZD-H group (SJZD high dose, 40 mg/kg, n � 6), and Urelumab group (20 mg/kg, n � 6). Mice in the SJZD groups were administered SZJD once a day, and mice in the Model group were given an equal volume of distilled water. e mice in the Urelumab group were injected with Urelumab twice a week. e tumor size was measured every week. After the administration of the drugs for 21 days, the mice were sacrificed, and the tumors were isolated, photographed, and weighed. e tumor growth inhibition weight value (TGIW) and tumor growth inhibition volume value (TGIV) were statistically analyzed. TGIW% � (1-Treatment group weight/Model group weight) × 100%; TGIV% � (1-Treatment group RTV/Model group RTV) × 100%; Relative tumor volume (RTV) � V 21 day /V 0 day ; Tumor volume � 0.5 × major axis × minor axis 2 .

2.6. H&E Staining. Tissues were fixed with 10% neutral buffered formalin overnight, embedded in paraffin, and cut into 4.5 μm sections using a microtome. After dewaxing and rehydration, the sections were stained with hematoxylin & eosin solution (R20570, Yuanye, Shanghai, China). Tumors were identified on sections. ese sections were sealed using neutral gum and photographed under a microscope (Nikon, Tokyo, Japan) at × 200.

Tumors were lysed using a radioimmunoprecipitation assay (RIPA) buffer to obtain the proteins. e protein concentrations were detected by a BCA protein assay kit. e proteins were separated by SDS/PAGE and blotted onto polyvinylidene fluoride (PVDF) membranes. Next, the membranes were blocked in skim milk for 2 hours and then incubated with primary antibodies against Bax (1 : 1000, 50599-2-ig, Proteintech, Wuhan, China), Bcl2 (1 : 1000, 12789-1-AP, Proteintech, Wuhan, China), p-PI3K (1 : 1000, 67071-1-Ig, Proteintech, Wuhan, China), p-Akt (1 : 1000, 10716-2-AP, Proteintech, Wuhan, China), and VEGFA (1 : 1000, 19003-1-AP, Proteintech, Wuhan, China) at 4°C overnight, followed by HRP-conjugated secondary antibodies. Protein bands were detected using an ECL chemiluminescence system (Beyotime, Shanghai, China). β-actin (1 : 1000, 66009-1-Ig, Proteintech, Wuhan, China) was used as an internal control.

After dehydration, hydration, and antigen recovery, tumor biopsies were incubated with the primary antibodies at the following dilution ratios: iNOS (1 : 800, 18985-1-AP, Proteintech, Wuhan, China) and COX-2 (1 : 800, 66351-1-Ig, Proteintech, Wuhan, China). Subsequently, the sections were incubated with a secondary antibody for 2 hours at room temperature. Finally, these sections were sealed using neutral gum and photographed under a microscope at a magnification × 200.

Analysis. Data are presented as the mean-± standard deviation. Differences between the two groups were investigated using Student's t-test. P < 0.05 was considered to indicate a statistically significant difference. Table 1 , there are 117 compounds in the SJZD in total (5 compounds in Baizhu, 6 compounds in Fulin, 90 compounds in Gancao, and 17 compounds in Renshen; Renshen and Gancao have one common compound). In addition, the SJZD and GC had 161 and 3288 potential targets, respectively. A total of 123 targets overlapped, which were the main related targets in the SJZD for the treatment of GC ( Figure 2 ).

To perform a comprehensive analysis, a TCM-compound-target-disease network was constructed (Supplement Table 2 and Figure 3 ). e network had 245 nodes, including 1 disease, 4 Chinese herbs, 117 compounds, and 123 targets. Each compound had multiple effective targets, such as quercetin (MOL000098), kaempferol (MOL000422), 7methoxy-2-methyl isoflavone (MOL003896), formononetin (MOL000392), and isorhamnetin (MOL 000354), with 60, 36, 34, 27, and 25 targets, respectively. On the other hand, each target also interacted with many compounds, such as the estrogen receptor (ESR1), calmodulin (CAMKK2), androgen receptor (AR), glycogen synthase kinase-3 beta (GSK3B), nuclear receptor coactivator 2 (NCOA2), coagulation factor Xa (F10), mitogen-activated protein kinase 14 (MAPK 14) , and beta-2 adrenergic receptor (ADRB2), acting on 78, 73, 69, 56, 50, 52, 47, and 38 compounds, respectively. e results demonstrated that the SJZD treats GC with multiple compounds and multiple targets. Evidence-Based Complementary and Alternative Medicine

e overlapping targets were input into Cytoscape software to achieve a visualized analysis. e network had 123 nodes and 1275 edges, and the average node degree was 20.7 (Supplement Table 3 and Figure 4 ). en, after topology analysis, it was found that the tumor necrosis factor (TNF, degree value � 67), interleukin-6 (IL-6, degree value � 66), cellular tumor antigen p53 (TP53, degree value � 63), transcription factor AP-1 (JUN, degree value � 62), and vascular endothelial growth factor A (VEGFA, degree value � 60) had higher degree values than the other targets. ese targets were mainly involved in the inflammatory response, malignant tumors, cell apoptosis, and proliferation.

To further explore the mechanisms of SJZD action for the treatment of GC, the overlapping targets were analyzed for GO functional and KEGG pathway enrichment. e enrichment results showed that there were 440 BP terms involved in the positive regulation of transcription from the RNA polymerase II promoter, signal transduction, and so on; 59 CC terms covered the extracellular space, plasma membrane, etc.; and 100 items of MF including protein binding, protein homodimerization activity, and so on were present. e top enrichment results are displayed in Figures 5(a)-5(c) . e related pathways for SJZD action for the treatment of GC were obtained by KEGG pathway enrichment analysis. A total of 110 pathways were enriched (P < 0.05), including pathways in cancer and the PI3K-AKT signaling pathway ( Figure 5(d) ).

After treatment with the SJZD for 21 days, the treated mice were sacrificed to isolate the tumor. e mice in the SJZD-H group had smaller and lighter tumors than the model group mice (P < 0.01, Figures 6(a)-6(c) ). In addition, the TGIW and TGIV results indicated that the SJZD-H treatment substantially inhibited tumor growth (Table 1 , P < 0.01). e H&E results showed that the tumor cells in the model group were closely arranged and grew well, with hyperchromatic nuclei and normal morphologies. After treatment with a high dose of the SJZD, the cells showed Evidence-Based Complementary and Alternative Medicine necrosis and exhibited diffuse and ambiguous distributions (Figures 6(b)-6(d) ). e results suggested that the SJZD had an inhibitory effect on gastric tumor proliferation.

According to the results of the PPI network and KEGG pathway enrichment analyses, specific targets, including VEGFA, Bax, Bcl2, and the PI3K/AKT pathway, were involved in the SJZD treatment of GC. erefore, we used WB to measure the expression of these proteins. As shown in Figures 7(a)-7(e), compared with the model group, the expression of VEGFA, p-PI3K, and p-AKT in the SJZD-H group was distinctly reduced along with an elevation in the ratio of Bax/Bcl2 (P < 0.01). e results suggested that the SJZD regulated the PI3K/AKT/VEGF pathway to treat GC.

Furthermore, the expression levels of iNOS and COX-2 were evaluated by immunohistochemistry. Figures 8(a) and 8(b) show that the expression levels of iNOS and COX-2 in the model group were distinctly higher than those in the SJZD-H group and Urelumab group (P < 0.01). ese results demonstrated that SJZD treatment downregulated the iNOS and COX-2 levels in the gastric tumor models.

In this study, network pharmacology and in vivo experiments were performed to reveal the bioactivity and mechanism of SJZD action for treating GC. e TCM-compoundtarget-disease network showed connections between multiple compounds and targets, which is consistent with the holism of TCM [22] . PPI results showed that TNF, IL-6, VEGFA, etc., played an important role in the SJZD treatment of GC. e GO functional and KEGG pathway enrichment results revealed the function of targets and potential mechanisms of the SJZD for the treatment of GC involving with pathways in cancer, such as the PI3K-AKT signaling pathway and so on. Furthermore, the experiment verified that the SJZD had a therapeutic effect on the gastric cancer model mice involving the VEGFA/PI3K/AKT pathways. e SJZD, a traditional Chinese prescription, has been used for the clinical treatment of gastrointestinal disorders for over 2000 years [23] . Due to its high therapeutic value, in recent years, the bioactive compounds in the SJZD have attracted a considerable amount of attention. Zhibo Guan et al. demonstrated that formononetin, ginsenoside, atractylenolide III, and so on were the characteristic components in the SJZD [24] . As shown in the analysis of the TCM compound-target-disease network, formononetin (MOL000392), ginsenoside (MOL005344), and atractylenolide III (MOL000072) were also the major active compounds in the SJZD treatment of GC. Furthermore, it has been reported that formononetin, ginsenoside, and atractylenolide III have antineoplastic efficacy in GC cells [25] [26] [27] . Our results also showed that the SJZD had an marked inhibitory effect on tumor cell proliferation, as evidenced by the reduced tumor weight and volume. In addition, pathology results also confirmed the SJZD-induced apoptosis of tumor cells.

e PPI network and KEGG results indicated that VEGFA was the main target in the treatment with SZJD against GC. VEGF is a well-known angiogenic factor and plays an important role in the proliferation, migration, and invasion of vascular endothelial cells [28] . e angiogenic pathway is a sound target for obstructing the excessive proliferation of cells because nutrients and growth factors are supplied through blood vessels to tumor cells [29] . Previous studies have indicated that during the development of GC, COX-2 and increased proangiogenic growth factors, particularly VEGF, have close relationships in which the downregulation of COX-2 inhibits VEGF generation [30, 31] . In addition, it has been reported that COX-2 overexpression in tumors significantly correlates with iNOS overexpression [32] . e nitric oxide (NO) produced through iNOS might increase vascular permeability by increasing the enzymatic activity of COX-2, accelerating the nutrient supply to tumor tissues, and ultimately promoting tumor growth [33] . In our animal experiments, we found that the SJZD suppressed the expression of VEGFA, COX-2, and iNOS to prevent tumor cell proliferation during GC development, which was consistent with the results of Yang et al. [34] .

To further explore the mechanisms of SJZD action for the treatment of GC, KEGG enrichment analysis was performed, which showed that the PI3K-AKT signaling pathway was the most enriched gene involved in the SJZD treatment of GC. e PI3K-AKT signaling pathway participates in many cellular activities, including cell growth, differentiation, apoptosis, and cytoskeletal rearrangement [35] . PI3K is an intracellular phosphatidylinositol kinase that can activate AKT, its downstream target, to regulate cell apoptosis during the development of human tumors [36] . LiRong et al. demonstrated that by targeting the PI3K/AKT pathway, the proliferation of GC in vivo and vitro can be inhibited, and cell apoptosis is induced by increasing the Bax/Bcl2 ratio [37] . Our results also emphasized that SZJD inhibited PI3K and AKT expression and increased the Bax/ Bcl2 ratio in a GC model.

In summary, our study combined network pharmacologybased predictions and in vivo experiments to verify the antineoplastic efficacy and the potential mechanisms of SJZD action for the treatment of GC. We found that the SJZD exerted its anti-GC effects by downregulating the expression of VEGFA, iNOS, and COX-2 to inhibit blood vessel hyperplasia and by upregulating the Bax/Bcl2 ratio to induce cell apoptosis by regulating the PI3K/AKT pathway. However, there were still deficiencies in this study. e potential compounds involved in the SJZD treatment of GC Evidence-Based Complementary and Alternative Medicine 9 need to be identified and validated in vivo or in vitro. In addition, the related targets should be further explored to provide a theory for the SJZD treatment of GC.

All data generated or analyzed during this study are included in this published article.

e authors declare no Conflicts of Interest.

Characteristics of gastric cancer in Asia

Tumor heterogeneity of gastric cancer: from the perspective of tumor-initiating cell

Worldwide trends in gastric cancer mortality

Gastric cancer: epidemiology, prevention, classification, and treatment

Alcohol and tobacco use in relation to gastric cancer: a prospective study of men in Shanghai, China

Helicobacter pylori and gastric cancer: a state of the art review

Current role of surgical therapy in gastric cancer

Neoadjuvant chemotherapy: survival benefit in gastric cancer

Ramucirumab and paclitaxel in patients with gastric cancer and prior trastuzumab: subgroup analysis from RAINBOW study

Cancer burden of major cancers in China: a need for sustainable actions

Weipiling ameliorates gastric precancerous lesions in Atp4a (−/−) mice

e traditional Chinese medicine baicalein potently inhibits gastric cancer cells

Sijunzi decoction-treated rat serum induces apoptosis of side population cells in gastric carcinoma

Progress in research on applying Sijunzi decoction in treating digestive malignant tumor

Effects of Chinese Jianpi herbs on cell apoptosis and related gene expression in human gastric cancer grafted onto nude mice

Network-based methods for prediction of drug-target interactions

e mechanisms of Sijunzi decoction in the treatment of chronic gastritis revealed by network pharmacology

A systems biology-based approach to uncovering molecular mechanisms underlying effects of traditional Chinese medicine Qingdai in chronic myelogenous leukemia, involving integration of network pharmacology and molecular docking technology

e important herbal pair for the treatment of COVID-19 and its possible mechanisms

UniProt: the universal protein knowledgebase

Effects of combined simultaneous and sequential endostar and cisplatin treatment in a mice model of gastric cancer peritoneal metastases

Progress in research on the effects of traditional Chinese medicine on the tumor microenvironment

Sijunzi Decoction attenuates 2,4,6-trinitrobenzene sulfonic acid (TNBS)-induced colitis in rats and ameliorates TNBS-induced claudin-2 damage via NF-κB pathway in Caco 2 cells

Rapid characterization of the chemical constituents of Sijunzi decoction by UHPLC coupled with Fourier transform ion cyclotron resonance mass spectrometry

Formononetin exhibits anticancer activity in gastric carcinoma cell and regulating miR-542-5p

Ginsenoside-Rh2 inhibits proliferation and induces apoptosis of human gastric cancer SGC-7901 side population cells

Atractylenolide III enhances the anti-neoplastic efficacy of docetaxel in gastric cancer cell by inhibiting fibroblast growth factor receptors 1, -2, and -4 expression

Lentivirus-mediated VEGF knockdown suppresses gastric cancer cell proliferation and tumor growth in vitro and in vivo

Revisiting tumor angiogenesis: vessel co-option, vessel remodeling, and cancer cell-derived vasculature formation

Molecular basis of therapeutic approaches to gastric cancer

e function and mechanism of COX-2 in angiogenesis of gastric cancer cells

Upregulation of cyclooxygenase-2 expression in lymphocytic thyroiditis and thyroid tumors

Relationship between the expression of iNOS, VEGF, tumor angiogenesis and gastric cancer

Phenethyl isothiocyanate inhibits migration and invasion of human gastric cancer AGS cells through suppressing MAPK and NF-kappaB signal pathways

Proliferation, survival and metabolism: the role of PI3K/AKT/mTOR signalling in pluripotency and cell fate determination

LncRNA AK023391 promotes tumorigenesis and invasion of gastric cancer through activation of the PI3K/Akt signaling pathway

Salidroside induces apoptosis and protective autophagy in human gastric cancer AGS cells through the PI3K/Akt/mTOR pathway