Hindawi Evidence-Based Complementary and Alternative Medicine Volume 2020, Article ID 3670309, 9 pages https://doi.org/10.1155/2020/3670309

Research Article Potential Molecular Mechanisms of Chaihu-Shugan-San in Treatment of Breast Cancer Based on Network Pharmacology

Kunmin Xiao,1,2 Kexin Li,1 Sidan Long,1 Chenfan Kong,1 and Shijie Zhu 1,2

1Graduate School, Beijing University of Chinese Medicine, Beijing 100029, China 2Department of Oncology, Wangjing Hospital, China Academy of Chinese Medical Sciences, Beijing 100102, China

Correspondence should be addressed to Shijie Zhu; [email protected]

Received 16 May 2020; Accepted 5 August 2020; Published 25 September 2020

Guest Editor: Azis Saifudin

Copyright © 2020 Kunmin Xiao et al. /is is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Breast cancer is one of the most common cancers endangering women’s health all over the world. Traditional Chinese medicine is increasingly recognized as a possible complementary and alternative therapy for breast cancer. Chaihu-Shugan-San is a traditional Chinese medicine prescription, which is extensively used in clinical practice. Its therapeutic effect on breast cancer has attracted extensive attention, but its mechanism of action is still unclear. In this study, we explored the molecular mechanism of Chaihu- Shugan-San in the treatment of breast cancer by network pharmacology. /e results showed that 157 active ingredients and 8074 potential drug targets were obtained in the TCMSP database according to the screening conditions. 2384 disease targets were collected in the TTD, OMIM, DrugBank, GeneCards disease database. We applied the Bisogenet plug-in in Cytoscape 3.7.1 to obtain 451 core targets. /e biological process of ontology (GO) involves the mRNA catabolic process, RNA catabolic process, telomere organization, nucleobase-containing compound catabolic process, heterocycle catabolic process, and so on. In cellular component, cytosolic part, focal adhesion, cell-substrate adherens junction, and cell-substrate junction are highly correlated with breast cancer. In the molecular function category, most were addressed to ubiquitin-like ligase binding, protein domain specific binding, and Nop56p-associated pre-rRNA complex. Besides, the results of the KEGG pathway analysis showed that the pathways mainly involved in apoptosis, cell cycle, transcriptional dysregulation, endocrine resistance, and viral infection. In conclusion, the treatment of breast cancer by Chaihu-Shugan-San is the result of multicomponent, multitarget, and multipathway interaction. /is study provides a certain theoretical basis for the treatment of breast cancer by Chaihu-Shugan-San and has certain reference value for the development and application of new drugs.

1. Introduction Traditional Chinese medicine (TCM) has a long history in the etiology, pathogenesis, prevention, and treatment of Breast cancer is one of the most common cancers en- breast cancer. According to the principle of TCM syndrome dangering women’s health all over the world. /e GLOBO- differentiation and treatment, the clinical syndrome of CAN 2018 statistics show alarming results that there are 8.6 breast cancer is mainly “Liver-Qi” stagnation. Chaihu- million new cases of female cancer and 4.2 million female Shugan-San is one of the classical prescriptions for the cancer deaths worldwide. /e proportion of breast cancer is treatment of “Liver-Qi” stagnation. It has the effect of 24.2% and 15.0%, respectively, ranking first in female soothing “Liver-Qi.” It has a history of 485 years and is cancer incidence and death [1]. It is predicted that, by the widely used in clinical practice [3–5]. Chaihu-Shugan-San 2050s, the global incidence of breast cancer will reach includes seven kinds of traditional Chinese medicine such as nearly 3,200,000 new cases of breast cancer each year. /ese Bupleurum chinense DC (Chinese name: Chaihu), Radix datasets reflect the high incidence of breast cancer and the Paeoniae Alba (Chinese name: Baishao), Citrus reticulata urgent global need for breast cancer prevention and Blanco (Chinese name Chenpi), Cyperus rotundus L treatment measures [2]. (Chinese name: Xiangfu), Glycyrrhiza uralensis Fisch 2 Evidence-Based Complementary and Alternative Medicine

(Chinese name: Gancao), Citrus aurantium L (Chinese Chinese Medicine Systems Pharmacology Database and name: Zhiqiao), and Ligusticum chuanxiong Hort (Chinese Analysis Platform (TCMSP, http://tcmspw.com/tcmsp.php, name: Chuanxiong) [6]. Traditional Chinese medicine updated on May 31, 2014) [17]. Search keywords are Chaihu, meridian tropism is one of the core components of the Baishao, Chenpi, Xiangfu, Gancao, Zhiqiao, and Chuan- theory of TCM. According to the theory of TCM, the xiong, and only oral bioavailability (OB) ≥30% and drug- characteristics of the selective distribution of the effective likeness (DL) ≥0.18 were considered in this study. /e components of traditional Chinese medicine in the body are Canonical SMILES sequence of the compound was searched basically consistent with the relationship between the viscera in the PubChem database (https://pubchem.ncbi.nlm.nih. and viscera of the corresponding meridian tropism. Chaihu, gov/) [18], and this sequence was used to predict the target of Chenpi, Xiangfu, Chuanxiong, and Baishao in Chaihu- the compound in the Swiss Target Prediction online data- Shugan-San belong to the liver meridian, which can soothe base (http://www.swisstargetprediction.ch/) [19] and collect the liver and regulate Qi well. In traditional Chinese target protein gene names and UniProt ID in the prediction medicine theory, the main location of breast cancer is in the results. liver and often related to the spleen and kidney in the process of disease. Chaihu-Shugan-San meridian attribution is consistent with breast cancer meridian attribution and can 2.2. Network Construction of Components and Targets. better play the therapeutic effect. /e chemical composition and potential targets of the above As a classical prescription, Chaihu-Shugan-San has been Chaihu-Shugan-San were uploaded to Cytoscape 3.7.1 [20] widely studied in pharmacology. Chaihu-Shugan-San has software to build up the component-target network. In the the pharmacological effects of antidepressant [7], regulation network, the degree centrality (DC) represents the number of neuro-endocrine-immune network [8, 9], anti- of nodes in the network that directly interacts with the node. inflammation, antioxidative stress [10], lipid-lowering and /e greater the degree, the more the biological functions it glucose-lowering, and antifibrosis [11]. In terms of antitu- participates in; the betweenness centrality (BC) refers to the mor, in a study of 86 patients with stage III breast cancer, the proportion of the number of nodes passing through the control group was given a standardized CAF regimen shortest path in the network, and the larger the BC is, the (cyclophosphamide + doxorubicin + 5-fluorouracil), and the more influential the node is. Closeness centrality (CC) re- experimental group was treated with Chaihu-Shugan-San on flects the degree of proximity between nodes, and the re- the basis of CAF. /e short-term effective rate of the ob- ciprocal of the shortest path distance from one node to other servation group was significantly higher than that of the nodes is CC. /e closer the nodes are, the larger the CC is; control group (81.4% vs. 58.14%); the Karnofsky improve- the average shortest path length (ASPL) is the average of the ment rate of the observation group was significantly higher shortest path length between all points in the network. /e than that of the control group (48.84% vs. 34.88%) [12]. As a smaller the average path of a node, the more crucial this safe complementary alternative therapy, Chaihu-Shugan- node is in the network. San combined with other chemotherapy regimens can im- prove the therapeutic effect, alleviate the myelosuppression 2.3. Prediction of Breast Cancer Targets. With “breast cancer” caused by chemotherapeutic drugs, and improve the or “malignant breast tumors” as keywords, we searched in prognosis of breast cancer patients [13, 14]. However, the Online Mendelian Inheritance in Man (OMIM, http:// mechanism of Chaihu-Shugan-San in the treatment of www.omim.org/, updated on May 4, 2018) [21], DrugBank breast cancer remains to be further explored. (https://www.drugbank.ca/, version 5.1.6, updated on Apr Classical prescriptions are currently the preferred way to treat 22, 2020) [22], /erapeutic Target Database (TTD, http://db. diseases in TCM clinic, but they lack a scientific basis to reasonably idrblab.net/ttd/, updated on Nov 11, 2019) [23], and Gen- explain the mechanism of TCM prescriptions from the whole to the eCards (https://www.genecards.org/, version 4.14.0) [24] to local level or from the cellular to the molecular level [15]. Network collect breast cancer-related targets. In the GeneCards da- pharmacology is an emerging discipline based on the integration of tabase, the higher the score value is, the closer the rela- systems biology, molecular biology, pharmacology, and a variety of tionship between the target and disease is, and the score network computing platforms in the context of the era of big data, value greater than the median is used as the screening which can more directly explain the association between TCM condition to extract the key target. /e above retrieval results prescriptions and diseases [16]. /erefore, this study constructed a were combined to remove duplicates and serve as the multidimensional network of “ingredient-target-pathway” through prediction target library of breast cancer. network pharmacology to explore the potential molecular mecha- nism of Chaihu-Shugan-San in the treatment of breast cancer and to provide a certain theoretical basis for Chaihu-Shugan-San in the 2.4. Protein-Protein Interaction (PPI) Network Construction treatment of breast cancer. and Selection of Core Targets. /e BisoGenet plug-in in Cytoscape 3.7.1 draws the PPI network and maps the 2. Materials and Methods Chaihu-Shugan decoction component targets and breast cancer-related disease targets into the protein interaction 2.1. Screening of Active Components and Target Prediction in relationship network, using Cytoscape 3.7.1 merge two Chaihu-Shugan-San. In this study, the chemical compo- protein interaction networks, to extract the intersection of nents of the seven herbs were searched on Traditional the network. Based on the intersection network, the Evidence-Based Complementary and Alternative Medicine 3

CytoNCA plug-in [25] is used to screen out the nodes whose disease targets were screened from Drugbank database, and degree centrality (DC) is greater than 2 times the median of 1286 disease targets with score >13.96 were obtained from all nodes. After multiple screening, the core PPI network is GeneCards database. /e duplicates were deleted after finally obtained. merging, and 2,384 breast cancer-related targets were finally obtained. 2.5. GO Functional Enrichment Analysis and KEGG Pathway Enrichment Analysis. /e core target of PPI network se- 3.4. Construction of PPI Network of Chaihu-Shugan-San and lected above was imported into Metascape (https:// Disease Targets. To further explore the pharmacological metascape.org/gp/index.html, updated on March 20, mechanism of Chaihu-Shugan-San on breast cancer, 2020) [26] database for KEGG (Kyoto Encyclopedia of Chaihu-Shugan-San and breast cancer protein were input and Genomes) pathway analysis and GO (Gene into the BisoGenet plug-in of Cytoscape 7.2.1 software for Ontology) biological process enrichment analysis. Param- merging. CytoNCA plug-in performs topological analysis eter is set to min overlap >3, p value cutoff <0.01, and min and takes 2 times of the average degree value as the screening enrichment >1.5. Taking p value as parameter and sorting condition. In the first screening, a network composed of from small to large as screening condition, KEGG pathway 2,728 nodes and 109,005 edges was obtained by the median and GO biological process of the top 20 eligible were selected DC >46. Finally, a PPI network with 451 nodes and 17,140 and uploaded to OmicShare (http://www.omicshare.com/ edges was constructed by further screening with the median tools) platform for data visualization. DC >156. /e process is shown in Figure 2. Topological parameters of the top 10 targets with high DC are shown in 2.6. Constructing PPI Network of Ingredient-Disease-KEGG Table 2, and other detailed results are shown in Table S2. Pathway. /e top 20 KEGG pathways, Chaihu-Shugan-San active ingredients, and disease common targets were 3.5. GO Biological Process and Enrichment Analysis of KEGG uploaded to Cytoscape 3.7.1 software to obtain the multi- Pathway. GO biological process consists of molecular dimensional network diagram of component-disease function (MF), biological process (BP), and cellular com- -KEGG pathway. ponent (CC) to interpret antitumor biological processes at key targets. Kyoto Encyclopedia of Genes and Genomes 3. Results (KEGG) enrichment analysis studies the key target of an- titumor signaling pathways. /e results of GO enrichment 3.1. Active Ingredient and Target of Chaihu-Shugan-San. analysis showed that there were 2000 biological processes, OB ≥30%, DL ≥0.18 as the screening conditions, after 274 cell components, and 564 molecular functions. KEGG searching TCMSP database, ChaiHu-ShuGan-San obtained enrichment has 154 signaling pathways. According to the a total of 157 chemical components, 13 compounds from ranking of p values, the top 20 were selected to plot the Baishao, 17 compounds from Chaihu, 5 compounds from bubble chart (Figure 3). /e left side of each chart is the top Chenpi, 7 compounds from Chuanxiong, 5 compounds enriched name. /e color of bubbles from blue to red from Zhiqiao, 92 qualified compounds from Gancao, and 18 represents the p value from large to small. /e larger the compounds from Xiangfu (as shown in Table S1 in Sup- bubbles, the larger the gene count of the pathway. /e plementary Materials). 8074 targets were obtained by in- horizontal axis represents the ratio of the pathway genes to putting 158 chemical components into Swiss Target the total input genes. /e top 20 signal pathways of KEGG Prediction online database. enrichment are shown in Table S3.

3.2. Compound-Target Network Construction. /e com- 3.6. 7e Multidimensional Network of “Component-Target pound-target network consists of 945 nodes and 8200 edges. Disease-KEGG Signaling Pathway” Was Constructed. 29 of 157 compounds were not found in the database and not Combining the component-target network and the first 20 involved in the network construction (Figure 1). In this KEGG signaling pathway targets, a multidimensional net- network, the average degree value is 15.647, and most of the work of “component-disease target-KEGG signaling path- proteins share common ligands with other proteins, which way” was obtained by Cytoscape 7.2.1 software (as shown in reflects the mechanism of the joint action between multi- Figure 4). /e results showed that the effective components components and multitargets of Chaihu-Shugan-San, and of Chaihu-Shugan-San could treat breast cancer by multi- conform to the characteristics of the traditional Chinese target and multisignal pathways. formula. Table 1 shows the detailed topological parameters of the top 10 compounds with high DC. 4. Discussion 3.3. Screening of Breast Cancer Targets. Breast cancer or Traditional Chinese medicine compound acts on diseases malignant breast tumors were used as keywords to search in through multimolecule, multitarget, and multipathway and TTD, OMIM, DrugBank, and GeneCards databases. 37 plays a certain therapeutic effect. Network pharmacology is a disease targets were obtained from TTD database, 1163 science based on the macroconnection under the back- disease targets were screened from OMIM database, 202 ground of the big data era. It systematically analyzes the 4 Evidence-Based Complementary and Alternative Medicine

MOL000359 MOL004328 MOL000354

MOL000358 MOL005828 MOL000422

MOL000098 MOL000211 MOL003542

MOL004891 MOL004903MOL004898 MOL004885MOL004884 MOL004904 MOL004883 MOL0000098 MOL004905 MOL004882 MOL004907 MOL004879 MOL001923 MOL004908 MOL004866 MOL004910 MOL004864 MOL004911 MOL004863 KIT NEK1 IGFBP3 GBA TACR1 CDK8 ATP12A KCNN3 ABCC9 PPARA PPARG FNTA MMP1 PIM1 AHR CYP3A4 CA1 WDR5 CHRNB4 MYLK MARK1 KCNN2 RORC CXCL8 KCNN1HSD17B14 ADRB3 HDAC2 BS MOL004912 MOL004857 MOL001930 MOL001924 PRCP LYPLA2 ALDH2 HRH4 LGALS3 CETP TTK DRD2 TLR4 DAPK1 TOP2A DHODH KLK1 PRKCZ FUCA1 GLRA1 CAMK2B PDE5A SCN2A ALOX15 HCRTR2 PRKCG PLA2G7 EPHA8 SCD ACVRL1 TYK2 HTT MOL004913 MOL004856 PRKACA ADRA2B ESRRG ECE1 UGT2B7 CTSV PREP WEE1 TBK1 PPP1CC HTR6 CXCR3 NR3C2 EPHA7 HPSE HSP90AB1 PTGER1 GSK3B MME FPR2 ADRA1A PRKCQ MIF ALPG PRKCH GRK6 CALCRL PTPRF

BDKRB1 MAPK9 AGPAT2 RXRA LNPEP PIM2 PTGDR PDE4A HDAC1 CYP2C9 SIRT3 CD38 ADRA2A EZH2 RBBP9 STS SREBF2 SERPINE1SMARCA2 CRHR1 KCNK9 PTPRS FAAH PLA2G2A SRD5A2 MAPK14 YES1 CCKBR MOL004914 MOL004855

TRPV1 MTOR RARA PPM1A PGR CASP6 PTK2 GIPR ATP4B PIK3CA BACE2 NTRK1 GRK7 MAPT OPRK1 CA2 PDE3B ESRRA CYP1B1 APEX1 GLI2 SLC28A3 RARB KCNQ3 ATM PIK3CD CDK6 IGFBP4 MOL004915 MOL004849 MOL001919 CYP2D6 CASP1 PYGL HPGDS DRD4 CDK2 QPCT NCOR2 PIP4K2C ACP1 TRPA1 PRKCA CASP8 FYN SOAT1 SSTR2 MMP14 HIPK4 RPS6KA4 PKM CALCA EPHA5 HDAC11 NPY2R NR1D1 HSD17B1 PDE2A FPGS MOL004917 MOL004848 PTGER2 AKR1C1 IMPDH1 ADRA1B PITRM1 CNOT7 TNNC1 GRIA2 MALT1 HTR7 NDUFA4 MDM2 HSD17B2 TNK2 CACNA2D1 DGAT1 S1PR4 TERT ATIC SLC5A2 TACR3 ROCK1 CYP51A1 CPT1A CA13 CCND1 CPT1B FABP4 MOL004924 MOL004841 GPER1 P2RY1 PDE7A PLK1 EP300 JAK1 RORA LTB4R GRK2 PDE4B TGFB1 AURKA MAP3K9 NOX1 ADAM33 NR1H4 TNKS MGLL TUBB3 DNTT RARG CHRM4 PIP5K1C S1PR1 PLA2G5 EPHA6 PAK4 ADORA2B

NAMPT AURKB BDKRB2 MMP13 DCTPP1 MPG DUSP3 CFTR CA12 AXL CXCR4 MAOB PRSS1 TAS2R31 BRD4 RAF1 PI4KB MAPK1 SSTR1 CASR SLC22A12 EPHA2 HSPA8 JAK3 SLC5A4 TYR MMP26 CYP17A1 MOL004935 MOL004838

TSPO SSTR3 CA14 LCK ADK CTSK SLC16A1 CXCR1 EPHA4 SLC22A6 FCER2 SAE1 ABCC1 BCL2 SLC6A3 CXCR2 MAP2K2 TKT CDK5 ADA ELANE CCNE1 MCL1 RIPK2 CA6 PDE4D TYMS GSR MOL004609 MOL013187 MOL004941 MOL004835 TRPM8 CBFB GART MPO DAGLA NR4A1 CCND3 BCHE EWS-Fli1 STAT6 FUT4 DAPK2 ABCG2 CCNE2 NR1I3 ASAH1 HDAC5 F10 SQLE CHRNA4 PARP2 EDNRA KIF11 HMGCR ANPEP NFKB1 FBP1 FGF2 MOL004945 MOL004833 MOL004598 MOL004624 SLC5A7 FEN1 ADRB1 CHEK2 THRA TYRO3 CDK9 CLK1 YWHAG GSTP1 PIM3 HSP90AA1 PTK2B SLC47A1 PGGT1B MMP12 EPHB2 PIK3C2G SGK1 FKBP1A SIGMAR1 BRS3 PIK3R1 IL2 AKT2 PLA2G10 PARP1 EPHX1 INSR MDM4 NOS3 SLC6A4 GABBR1 USP1 BLK AMY1A XIAP NAAA KISS1R KLK2 RPS6KA2 FOLR1 SF3B3 EPHB3 KCNA3 MAPK8 CYP27A1 HTR3A CHEK1 CSF1R SSTR4 NOX4 DYRK1B SELP SIRT2 RXRG MOL004948 MOL004828 GC GNRHR GABRG2RASGRP3 CMA1 GSTM2 IGFBP2 PTGER3 CCKAR HAO1 BAD GPBAR1 DRD1 HDAC4 NOS2 PRKCE LDLR GRIA1 CA3 AHCY LIMK1 PDPK1 BCL2L1 CCNB3 MMP25 PLAT GCGR PPP2R5A ITK MOL000049 MOL004644 MOL004949 MOL004827 CH SLC5A1 IMPDH2 KCNK3 PRKG2 MAPK10 THRB ACVR1 HCK PTGIR ACACB HNF4A HTR1B FDFT1 MAP4K4 ABAT TXK S1PR5 DHFR ADRB2 EPHA1 FLT4 IGFBP1 DAO ADRA2C VCP HDAC8 CTRC ST3GAL3

PTPN6 BTK ENPP1 FGFR2 MPI RET CA4 DRD3 TNKS2 TBXA2R PGK1 ROS1 PIK3CB ALPL AURKC GRIN2B IKBKE ARG1 STAT1 MMP8 HPGD MAOA PRKDC FASN IDO1 HK1 MAPKAPK2 PLK2 MOL004957 MOL004824 CBR1 MAP2K4 AOC3 MELK MMP2 CES1 EGLN1 MAP3K7 KDR NPC1L1 ACHE CSNK1G2 CREBBP CYP2C19 LIMK2 POLB CHRM1 CA7 COQ8B SCARB1 SHBG KDM4E MGMT PTPN1 SRD5A1 HTR2B KCNMA1 C5AR1 MOL002776 MOL004653 MOL004959 MOL004820 CFD SLC2A1 SLC1A3 PSEN2 IDH1 BRAF PDE4C CYP11B1 CDC7 TUBB1 MAST3 RPS6KB1 NR1H3 DPP4 SLC6A2 PLG ADCY10 CCR4 ABL1 P2RX3 AKR1C2 ADRA1D PTP4A3 MMP10 HDAC6 MERTK MGAM PDE6D

PRF1 PHLPP2 SCN9A ADAM17 GSK3A HTR1A PPP2CA HRH3 MMP16 ERN1 KDM4D SLC9A1 JUN MMP15 TOP1 EDNRB EIF2AK3 F3 IGF1R BMP4 DPYD HDAC9 PLAU ADAMTS5 CA9 IKBKB EPHX2 AKT1 MOL004961 MOL004815 MOL001645 MOL0004718 GRM4 GABRB3 PLK3 FABP5 FUT7 WNT3A NQO2 FGR CNR1 MKNK1 PDE1B GRM2 LGALS9 GAPDH CDC25A PTAFR VDR PTGS1 EIF4H LIPE F2R MAP3K5 IRAK4 SCN5A CTSB AVPR2 CDC25C SLC29A1 MOL000490 MOL004966 MOL004814 CALM1 ERBB2 G6PD OPRM1 ERCC5 CDK4 SSTR5 FOLR2 PLAA CAPN1 CASP3 STK17B PROKR1 DNM1 REN IGFBP5 RPS6KA5 KCNH2 PNP GPR35 TDP2 PTGES PDGFRA MAPK3 FGFR1 HK2 APP MAP3K12 MOL004974 MOL004811 ALOX12 SLC28A2 VEGFA CES2 ABHD6 SHH ODC1 BACE1 CLK2 TLR9 UTS2R HSD11B1 CLK3 LDHA SNCA SLC46A1 CTSL CDK5R1 MMP9 CASK TAAR1 CCNT1 GCK GUSB RBP4 HSD17B3 PLCG1 PPM1B

ATR PDF NEK6 PPIA KDM4C KAT2B DBF4 DHCR7 RELA CA5B DSTYK AKR1A1 PRKD1 AGTR1 NCOR1 SLC6A9 CYP1A2 TTR STAT3 LTA4H EGFR TBXAS1 ADCY5 AVPR1A CCNA2 HDAC7 ALOX5 BMX MOL004978 MOL004808

SERPINA6 MAP2K3 IGFBP6 CYP11B2 SIRT1 SLC18A2 CTSD PDE6C TGM2 HSD11B2 JAK2 DYRK2 PABPC1 KDM5A PTPN2 EPHB4 S1PR3 PLA2G1B FLT1 CYP19A1 CNR2 MPEG1MAPKAPK5 ITGAL DRD5 DYRK3 NR1H2 PIK3CG MOL004980 MOL004805 MOL004058 MGAT2 LYPLA1 PDE1A GRM1 PGF TNF PDGFRB GRIN2A PRKCD ESR1 HSPA5 NOS1 CD81 ALK ACPP HDAC10 ILK LYN PTGS2 CASP7 PSEN1 RPS6KA3 TACR2 SRC MMP7 ME1 FABP2 SOAT2 MOL004985 MOL003896 MOL004068 MOL004053 ESRRB EPHA3 EPHB1 KDM5C LRRK2 AKR1C4 ABCB1 EPAS1 AKR1C3 GLO1 AKR1B1 P2RX7 ADAMTS4 XDH KCNA5 PLEC CHRNA7 CHRM2 FABP3 CSNK1D PLK4 FGFR3 PDE10A PTPRCAP COMT MET KDM4A HRAS MOL004990 MOL003656 TGFBR1 GLI1 ADORA2AHSP90B1 PRKCB KDM1A SPHK2 ESR2 OXTR DNMT1 NEK2 MT-CO2 METAP2 PTGER4 DUT DMPK MBD2 FOLH1 HDAC3 RPS6KA1 CYP1A1 ADAM10CDC42BPA PKN1 MARS HCRTR1 ADORA1 MTNR1B MOL004991 MOL002565 XF PTPN11 NAE1 LGALS7 FFAR1 FGF1 PDE3A TNNI3K BMP1 OPRD1 HIF1A ALOX15B TAOK1 GABRA2 EZR CA5A ROCK2 CSNK1A1AKR1B10 TAOK3 NUAK1 MTNR1A PDK1 CDC25B F2 SYK XPO1 CLK4 GPR84 MOL004996 MOL002311 MOL004071 MOL003044 SPHK1 NR3C1 CSNK1G1 SLC19A1 PFKFB3 ADAM9 CAPN2 GSTA1 GRM5 ADORA3 CHRM3 CDK7 MAP2K1 PPARD ALOX5AP CTSS AR EIF4A1 HTR2C DAPK3 FABP1 CHRNA3 PHKG2 PTK6 TYMP HTR2A TAOK2 TEK MOL005000 MOL001792 MOL005001 MOL001484 FLT3 DYRK1A MAP3K8 CCR2 MMP3 CSK CDK1 NPY5R OGA CCR1 ST6GAL1 PON1 PRKD2 CSNK2A1 MOL005003 MOL000500 MOL005012 MOL000497 MOL004074 MOL000449 MOL005016 MOL000417 MOL005017MOL000239MOL000392 MOL000006

MOL002135 MOL001494 CX ZQ CP

MOL002341 MOL005100 MOL000433

Figure 1: Compound-target network: green diamond nodes represent compounds, red diamond nodes represent common compounds, yellow v nodes represent herb names, and purple circular nodes represent corresponding potential targets of compounds. molecular mechanism of action from all levels, which is consistent with the holistic view of TCM and the thought of Table 1: Topological parameter of top 10 compounds. syndrome differentiation and treatment. Chaihu-Shugan- San prescription, in which Chaihu is the monarch drug, is ID Molecule name DC BC CC ASPL good at soothing the “Liver-Qi.” Xiangfu and Chuanxiong MOL000422 Kaempferol 317 0.011 0.377 2.651 are the minister drugs, which can relieve “Liver-Qi.” Chenpi MOL000354 Isorhamnetin 303 0.008 0.373 2.681 and Zhiqiao, regulating “Qi” stagnation and Baishao, MOL000359 Sitosterol 233 0.005 0.360 2.774 nourishing “Blood” and softening the “Liver,” are adjuvants. MOL004328 Naringenin 199 0.058 0.397 2.520 Gancao is used as a guide medicine to reconcile various MOL000358 Beta-sitosterol 139 0.004 0.357 2.802 drugs. /e combination of various drugs can regulate “Liver- MOL000098 Quercetin 116 0.009 0.377 2.654 Qi” and smooth “Qi.” Chaihu-Shugan-San has a history of MOL004609 Areapillin 101 0.008 0.372 2.690 MOL003044 Chrysoeriol 101 0.008 0.371 2.690 more than 480 years. /e basic compatibility of clinical MOL000006 Luteolin 101 0.008 0.371 2.690 medication is Chenpi 6 g, Chaihu 6 g, Chuanxiong 4.5 g, MOL004071 Hyndarin 101 0.063 0.374 2.677 Xiangfu 4.5 g, Zhiqiao 4.5 g, Shaoyao 4.5 g, and Gancao 1.5 g, which is added or subtracted according to the actual situ- ation of patients. /e official preparation method is to add breast cancer cells by activating the Fas signaling pathway 250 ml of water to the herbs and boil them for 30 min. Su and caspase-8 activity [32] and is expected to be an orphan [27] and his team identified 33 chemical constituents in nutrition drug against cancer [33]. Kaempferol has shown a Chaihu-Shugan-San. Among them, gallic acid (source: good affinity for PAK4 in molecular docking and is con- Shaoyao), oxidized paeoniflorin (source: Shaoyao), paeo- sidered to be a potential inhibitor in triple-negative breast niflorin (source: Shaoyao), paeoniflorin (source: Shaoyao), cancer [34], and kaempferol can prevent G2/M phase of the glycyrrhizin (source: Gancao), naringin (source: Chenpi, cell cycle by downregulating CDK1 in human breast cancer Zhiqiao), hesperidin (source: Chenpi, Zhiqiao), and ferulic MDA-MB-453 cells [35], and blocking RhoA and Rac1 acid (source: Chuanxiong) had higher contents, all above signaling pathways to inhibit breast cancer cell migration 1000 mg/g [28]. Although some studies have comprehen- and invasion [36] is a powerful antioxidant inducer and can sively elucidated the treatment of depression [6, 29], non- inhibit oncogene transformation and induce cancer cell alcoholic fatty liver disease [30, 31] and functional dyspepsia apoptosis and DNA damage. Quercetin, naringenin, and [5] by Chaihu-Shugan-San, no studies have comprehen- isorhamnetin, such as flavonoids, can prevent breast cancer sively elucidated the mechanism of Chaihu-Shugan-San in cell migration through inflammatory and apoptotic cell the treatment of breast cancer. /erefore, with the aid of signaling [37, 38], and quercetin can induce autophagy by network pharmacology, this study analyzed the specific inhibiting the Akt-mTOR pathway [39]. molecular mechanism of Chaihu-Shugan-San in the treat- /e PPI network showed that the active components in ment of breast cancer from a microscopic perspective. Chaihu-Shugan-San may function through the core targets /e results of network analysis showed that the active such as histone deacetylase 1 (HDAC1), huntingtin (HTT), ingredients in Chaihu-Shugan-San mainly included beta- RAC-alpha serine/threonine-protein kinase (AKT1), hepa- sitosterol, kaempferol, quercetin, naringenin, isorhamnetin, toma-derived growth factor (HDGF), roquin-1 (RC3H1), and nobiletin. Beta-sitosterol can promote the apoptosis of chromobox protein homolog 8 (CBX8), histone deacetylase Evidence-Based Complementary and Alternative Medicine 5 DC > 46 DC > 156

9787 nodes and 215324 edges 2728 nodes and 109005 edges 451 nodes and 17140 edges Figure2: Network topology analysis of PPI.

treatment of breast cancer by Chaihu-Shugan-San is the Table 2: Topological parameter of top 10 core targets. result of multimolecular, multitarget, and multipathway interaction. Target DC BC CC ASPL /e biological process of (GO) involves HDAC1 1976 0.075 0.534 1.874 the mRNA catabolic process, RNA catabolic process, telo- HTT 1695 0.094 0.526 1.900 mere organization, nucleobase-containing compound cat- AKT1 1136 0.031 0.495 2.019 abolic process, heterocycle catabolic process, and so on. In HDGF 1095 0.021 0.499 2.002 cellular component, cytosolic part, focal adhesion, cell- RC3H1 930 0.026 0.501 1.996 substrate adherens junction, and cell-substrate junction are CBX8 836 0.013 0.488 2.050 highly correlated with breast cancer. In the molecular HDAC2 813 0.015 0.493 2.029 SUMO1 810 0.021 0.500 1.999 function category, most proteins were addressed to ubiq- RPSA 761 0.013 0.489 2.047 uitin-like protein ligase binding, protein domain specific RPLP0 736 0.011 0.481 2.081 binding, and Nop56p-associated pre-rRNA complex. In addition, the results of KEGG pathway analysis showed that the pathways mainly involved in apoptosis, cell cycle, 2 (HDAC2), small ubiquitin-related modifier 1 (SUMO1), transcriptional dysregulation, endocrine resistance, and viral 40 S ribosomal protein SA (RPSA), and 60 S acidic ribosomal infection. Estrogen receptor (ER) signal transduction protein P0(RPLP0). HDAC1 plays an important role in pathway plays a central role in the development of breast transcriptional regulation and cell cycle progression [40]. cancer. ER can not only regulate the expression of certain HDAC1 can promote the proliferation and migration of genes through its mediated signal transduction pathway but breast cancer cells by activating the Snail/IL-8 signaling also has extensive connections with many other signal pathway [41]. Downregulation of HTT transcription and transduction pathways, forming a signal transduction reg- protein levels is a key factor in poor prognosis and metastasis ulatory network [45]. ER binds to receptor proteins in the development of breast cancer [42]. AKT1 is involved in the nucleus, and the receptor is activated. Activated ER-α and regulation of many tumor processes, including tumor ER-β form homodimers or heterodimers. Some coregulators proliferation, cell survival, metabolism, growth, and an- form complexes with dimers. /e complexes bind to ER giogenesis. /e mutation frequency of AKT1 in Chinese response elements to initiate transcription, thereby regu- breast cancer patients is 3.2%, and it is considered to be a lating the function of target genes, leading to abnormal cell sensitive target for the treatment of breast cancer. A study proliferation and differentiation, and ultimately leading to involving 313 Chinese breast cancer patients found that the tumorigenesis [46]. /e increased mutation rate of ER-α in mutation frequency of AKT1 in Chinese breast cancer pa- precancerous lesions of breast cancer affects the junction tients was 3.2%, and it is considered a sensitive target for the between ER-α zinc finger region and ligand binding domain, treatment of breast cancer [43]. HDAC2 is a poor prognostic resulting in high sensitivity of the body to estrogen. Under factor in patients receiving anthracycline therapy and is the action of low levels of hormones, ER-α is highly bound to positively correlated with breast cancer metastasis, pro- TNF-2 co-activator, which leads to the occurrence of tu- gression, increased Ki-67, multidrug resistance protein, and mors. For ER receptor-positive breast cancer patients, negatively correlated with overall survival of patients [44]. quantitative expression of ER receptor is an independent /e occurrence and development of breast cancer are closely imaging factor to evaluate their prognosis, recurrence, and related to the core proteins, which fully prove that the metastasis [47]. Abnormal activation of MAPK signal 6 Evidence-Based Complementary and Alternative Medicine

Top 20 of GO enrichment Top 20 of GO enrichment

GO:0071824:protein−DNA complex subunit organization GO:0015935:small ribosomal subunit

GO:0019221:cytokine−mediated signaling pathway GO:0044815:DNA packaging complex GO:1903320:regulation of protein modification GO:0000786:nucleosome by small protein conjugation or removal GO:0019941:modification−dependent protein catabolic process GO:0022627:cytosolic small ribosomal subunit

GO:0044257:cellular protein catabolic process GO:1990234:transferase complex

GO:0019083:viral transcription GO:0044391:ribosomal subunit GO:0000184:nuclear−transcribed mRNA catabolic GO:0016604:nuclear body process, nonsense−mediated decay GO:0006281:DNA repair GO:0005840:ribosome

GO:0065004:protein−DNA complex assembly GO:0000788:nuclear nucleosome

GO:0097190:apoptotic signaling pathway GO:0000781:, telomeric region

GO:0010564:regulation of cell cycle process GO:0000784:nuclear chromosome, telomeric region GO term GO term

GO:1901361:organic cyclic compound catabolic process p value GO:0098687:chromosomal region 7.947150e − 43 p value 5.12632e − 21 GO:0019439:aromatic compound catabolic process 5.960362e − 43 GO:0032993:protein−DNA complex 3.84474e − 21 3.973575e − 43 GO:0044270:cellular nitrogen compound catabolic process GO:0022626:cytosolic ribosome 2.56316e − 21 1.986787e − 43 1.28158e − 21 GO:0046700:heterocycle catabolic process 3.457230e − 63 GO:0070161:anchoring junction 4.98804e − 50 GO:0034655:nucleobase−containing compound catabolic process GO:0005912:adherens junction

GO:0032200:telomere organization GO:0030055:cell−substrate junction

GO:0006401:RNA catabolic process GO:0005924:cell−substrate adherens junction

GO:0006402:mRNA catabolic process GO:0005925:focal adhesion GO:0072594:establishment of protein GO:0044445:cytosolic part localization to organelle

0.1 0.2 0.3 0.4 0.2 0.4 0.6 Rich factor Rich factor

Gene number Gene number 50 80 30 60 60 90 40 70 70 50 (a) (b)

Top 20 of GO enrichment Top 20 of KEGG Enrichment

GO:0001085:RNA polymerase II transcription factor binding hsa01522:Endocrine resistance

GO:0005198:structural molecule activity hsa05202:transcriptional misregulation in cancer

CORUM:5266:TNF−alpha/NF−kappa B signaling complex 6 hsa04068:FoxO signaling pathway

GO:0051427:hormone receptor binding hsa04010:MAPK signaling pathway

GO:0042826:histone deacetylase binding hsa04919:thyroid hormone signaling pathway

GO:0003682:chromatin binding hsa04210:apoptosis

CORUM:305:40S ribosomal subunit, cytoplasmic hsa05215:prostate cancer

CORUM:338:40S ribosomal subunit, cytoplasmic hsa05166:HTLV−I infection

GO:0046982:protein heterodimerization activity hsa05212:pancreatic cancer

GO:0003735:structural constituent of ribosome p value hsa04120:ubiquitin mediated proteolysis 4.266670e − 19 p value 5.849230e − 22

GO:0050839:cell adhesion molecule binding e Pathway hsa05168:herpes simplex infection GO term 3.200003 − 19 4.386922e − 22 2.133335e − 19 GO:0045296:cadherin binding hsa05220:chronic myeloid leukemia 2.924615e − 22 1.066667e − 19 1.462308e − 22 CORUM:306:ribosome, cytoplasmic 2.554400e − 78 hsa03010:ribosome 4.782950e − 67 GO:0008134:transcription factor binding hsa04110:cell cycle

GO:0019901:protein kinase binding hsa05205:proteoglycans in cancer

GO:0019900:kinase binding hsa05034:alcoholism

CORUM:3055:Nop56p−associated pre−rRNA complex hsa05200:pathways in cancer

GO:0019904:protein domain specific binding hsa05161:hepatitis B

GO:0031625:ubiquitin protein ligase binding hsa05169:Epstein−Barr virus infection

GO:0044389:ubiquitin−like protein ligase binding hsa05203:viral carcinogenesis

0.00 0.25 0.50 0.75 0.1 0.2 0.3 0.4 Rich factor Rich factor

Gene number Gene number 25 75 30 50 50 100 40 60 (c) (d)

Figure 3: GO function enrichment analysis and enrichment analysis of KEGG signaling pathway (top 20). (a) BP. (b) CC. (c) MF. (d) KEGG. transduction pathway can lead to cell loss of apoptosis and migration, and apoptosis of breast cancer cells through differentiation ability, promote malignant transformation, MAPK and are expected to become new targets for the abnormal proliferation, produce tumors, and further pro- treatment of breast cancer. Studies found that the activation mote the proliferation of tumor cells. /erefore, inhibitors of or loss of FOXO function can inhibit the growth and me- some key kinases in the MAPK signaling pathway have tastasis of breast tumors [50], and the dysregulation of become a hotspot in the treatment of breast cancer in recent FOXO transcription factors has also become a key molecule years. Studies have found that Kruppel-like factor 4 [48] and in the endocrine resistance mechanism [51]. More and pre-mRNA processing factor 4 [49] affect the growth, more attention has been paid to the relationship between Evidence-Based Complementary and Alternative Medicine 7

AHR EPHB2TOP1 VEGFA ALDH2 TYMPBRAF DUSP3 HPSE MAP2K2 EPHA2 DHFR SREBF2 TLR4 COMT KDM4C NR1H2 IGFBP2 NCOR2 CYP2D6 FGFR3 PGR LIMK1 HSD11B2 FABP4TYMSMPO CDK7SERPINE1 CALCACTSD F2 CYP11B1 CDC25C CASP7 HDAC9 ESR2 CBFB H4C11H4C12H4C4H2BC10 MAPK9 TTR ABCC1 H4C2H4C8 H2BC8 H4C5 USP7 JAK3 TLR9 VCP H4C13 YWHAZ PRSS1 YWHAG DNMT1 IDH1 RPS6KA3 H4C14 YWHAE YWHAQ MTOR CTSB SLC6A3 YWHAB H4C15 MMP2 KISS1RDAPK1ROS1PTPN1HMGCR G6PD KIT TRAF2 PLAU CDC25A CTSL CHUK HDAC3IKBKGH4C3 PLK3 WEE1 KAT2B H4C6 RB1 VDR IGFBP5 TGFB1 CCKBR HSPA5 SNW1 H4C1 CSNK2A2 CCNE2 AURKB NOS3 H4-16 H4C9 RAC1 BAD PIK3CD BCL2 HSPA1B TOP2A AURKA BMP1 CSNK2B YWHAH LYN CCND1 SLC6A4 HSPB1 SLC6A2 CA2 EPHA3 HLA-B HSPA1A TRAF1 F3 HSD17B1 FGFR1 ROCK2 CHEK1 ICAM1 TYR H2BC5 MYC HSPA8 REL PSEN2 SLC19A1 TKT ESRRA PARP1 RET GRB2 CYP19A1 RARB NEDD4 PSMC5 PXN RPS6KB1 DPYD ALOX5 DDB1 SRD5A2 SMARCA2 MOL004898MOL004904MOL004974MOL004833MOL001792MOL004905MOL003896FASN POLR2A Chronic Regulation STAT6 LGALS3 PTGS2 PTEN Prostate Focal NFKBIA ATF2 myeloid of actin HIF1A ATM STAT1 HDAC6 PRKCA ESR1 PRKACA SLC29A1 MAPK10 FLT1 PSMC3 cancer adhesion MOL000098MOL004058MOL004978MOL004848MOL002341MOL004910MOL004985 leukemiaSignaling cytoskeleton MMP9 PPARG NR3C1 PML pathways HNRNPK CDKN2A FGF2 TGM2 PSMD2 SIRT1 PLCG1 PTPN11 MAPK14MAPK1 IKBKB PKM MMP1 regulating Colorectal Pancreatic PRKCQ MOL000422MOL000211MOL004071MOL001923MOL001919MOL000392MOL004911 ODC1 Alcoholism CDKN1A TNF SPHK2 PTK6 CUL2 pluripotency cancer cancer DDX3X EGFR HDAC4 CDK1 HDAC2 PIK3R1 EZR JUN PSMD4 of stem CFTR CDC42 PIK3CA MOL004053MOL004841MOL004838MOL004849MOL003656MOL004996MOL004924 FGF1 cells Epstein-Barr EP300 CASP3 IKBKE MAPK3 HSP90AA1AR CDK4 MAP2K1 JAK1 GAPDH TRAF6 Proteoglycans HTLV-I Viral PTK2B SHC1 virus IQGAP1 MOL005003MOL005000MOL000359MOL000500MOL004857MOL004814MOL004991 in cancer infection carcinogenesis MYH9 NTRK1 RELA ILK ERBB2 SRC PLK1 HDAC5 MYLK CCND3 LDHA CXCR3 CSF1R VIM infection Bacterial CCNA2 PI3K-Akt Systemic CLTC RARA MOL005012MOL004966MOL004864MOL001484MOL004856MOL000433MOL001494 STS CANX invasion of Adherens FLNA CDK2 RAF1 MDM2 GSK3B CREBBP NFKB1 STAT3 IGFBP4 EDNRB MMP3 XPO1 signaling lupus epithelial junction EDNRA pathway erythematosus SSX2IP HDAC1 PRKCD MAPK8 TGFBR1ABL1 AKT1 CASP8 CHRNA3 NOS2 MOL000497MOL000239MOL004935MOL004828MOL013187MOL002311MOL004879FABP3 INSR FOS MAP3K3 cells AGE-RAGE DDX5 CHEK2 FoxO signaling Pathways in VCL BTK PRKD1 HRAS PDPK1 CDK6 signaling pathway in Hepatitis B MOL004908MOL004882MOL004959MOL004913MOL002565MOL004885 SMAD3 cancer PTK2 VCAM1 pathway diabetic H3C7 FYN JAK2 EZH2 IL2 MMP13 SHBG SMAD4 complications AKR1C3 KAT5 H3C2 PLA2G2A WNT3A ERBB3 HPGD EPAS1 CYP1A2 MAP3K1 IGFBP1 BRS3 PIK3CB PRMT1 H3C12 CAV1 AXL ALK NFATC2 CXCR4 MMP14 PKN1 ACTG1 IGF1R PCNA PRKAB1 H3C8 CXCR2 CSNK1D GSTP1 XIAP KDM1A BRCA1 H3C6 ACTB PDGFRB DRD2 TERT MDM4 AKT2 APC IGFBP3 PLA2G10 SYK PPP2CB H3C4 H3C10 SHH MAP2K4 ANPEP GPER1 ESRRG CBL MMP7 FLT3 CYP3A4 PPP2R1A SOX2 H3C11 GLI1MCL1 CDC25BTEK CRK RPS6 POU5F1 BCL2L1 FOLH1 ADA ERCC5 CDC37BMI1 SMAD1 H3C3 SLC2A1 CTNNB1 SSTR2 SRD5A1 H3C1 PTGS1 FN1 CXCL8 KDR GNRHR H3-4 FLT4 HSP90AB1 BMP4 NCOR1 PPP1CC KCNH2 S1PR3 ABCG2 SMAD2 MGMT ADAM17 MAP3K8 VHL PPP1CB CA9 PTPRF SPHK1 RBX1 PPP1CA ATR CYP1B1EPHB4CCNE1PON1ABCB1 TNFRSF1A NANOGCOMMD3-BMI1 RARG FBXO25ITGA4PPP2CARPS6KB2POU5F1B APEX1 EPHA4 MIF FGFR2 PIK3CG GLI2 PRKCB GRM1 NQO2 HSD17B3 PLG F2R ROCK1CYP1A1 CYP17A1 PDGFRACSK NR3C2MET MME Figure 4: Component-disease target-KEGG signaling pathway. Green represents active ingredients of drugs; purple represents targets; yellow represents common targets; and red represents signaling pathways. viral infection and breast cancer. In particular, human Conflicts of Interest papilloma virus (HPV) has a strong cause-and-effect re- lationship with breast cancer. Many studies have found that All authors state that they have no conflicts of interest re- different HPV genotypes are associated with the prevalence garding the publication of this paper. of breast cancer and the nuclear prognosis. /e relationship between viral infection and breast cancer has been paid Authors’ Contributions more and more attention [52–54]. /e relationship be- tween Epstein–Barr Virus (EBV) and breast cancer has also Shijie Zhu conceived and designed the experiments. Kunmin been extensively studied, but the current evidence is less Xiao and Kexin Li performed the experiments and wrote the and more controversial [55]. It is proved again that Chaihu- manuscript. Sidan Long and Chenfan Kong contributed to Shugan-San treatment of breast cancer is through a analysis tools. combination of multiple biological pathways and multiple signaling pathways, but this multifeature is not only found Acknowledgments in a single disease. Chaihu-Shugan-San is mainly involved in the regulation of neurotransmitters, regulation of in- /is work was supported by grants from the China National flammatory mediators of TRP channels, calcium signaling Natural Science Foundation (Grant no. 81973640). pathways, cyclic adenosine monophosphate signaling pathways, and neuroactive ligand-receptor interactions to Supplementary Materials play an antidepressant role [56]. Chaihu-Shugan-San can improve cognitive impairment in Alzheimer’s disease Table S1: active ingredients of Chaihu-Shugan-San. Table S2: through multitarget action, and its effect is verified by topological parameters of Chaihu-Shugan-San targets. Table biological experiments [57]. /ese all embody the principle S3: the top 20 signal pathways of KEGG enrichment. of “treating different diseases with the same treatment” in (Supplementary Materials) TCM. References Data Availability [1] F. Bray, J. Ferlay, I. Soerjomataram et al., “Global cancer statistics 2018: GLOBOCAN estimates of incidence and All data generated or analyzed during this study are included mortality worldwide for 36 cancers in 185 countries,” CA: A in this paper. Cancer Journal for Clinicians, vol. 68, no. 6, pp. 394–424, 2018. [2] J. R. Benson and I. Jatoi, “/e global breast cancer burden,” Future Oncology, vol. 8, no. 6, pp. 697–702, 2012. Disclosure [3] N. Yang, X. Jiang, X. Qiu, Z. Hu, L. Wang, and M. Song, “Modified Chaihu shugan powder for functional dyspepsia: Kunmin Xiao and Kexin Li are the first authors. meta-analysis for randomized controlled trial,” Evidence- 8 Evidence-Based Complementary and Alternative Medicine

Based Complementary and Alternative Medicine, [17] J. Ru, P. Li, J. Wang et al., “TCMSP: a database of systems vol. 201310 pages, Article ID 791724, 2013. pharmacology for drug discovery from herbal medicines,” [4] Y. Wang, R. Fan, and X. Huang, “Meta-analysis of the clinical Journal of Cheminformatics, vol. 6, no. 1, 2014. effectiveness of traditional Chinese medicine formula Chaihu- [18] W.-D Ihlenfeldt, “A virtual file system for the PubChem Shugan-San in depression,” Journal of Ethnopharmacology, chemical structure and bioassay database,” Chemistry Central vol. 141, no. 2, pp. 571–577, 2012. Journal, vol. 2, no. 1, 2008. [5] F. Qin, J. Y. Liu, and J. H. Yuan, “Chaihu-Shugan-San, an [19] D. Gfeller, O. Michielin, and V. Zoete, “Shaping the inter- oriental herbal preparation, for the treatment of chronic action landscape of bioactive molecules,” Bioinformatics, gastritis: a meta-analysis of randomized controlled trials,” vol. 29, no. 23, pp. 3073–3079, 2013. Journal of Ethnopharmacology, vol. 146, no. 2, pp. 433–439, [20] P. Shannon, A. Markiel, O. Ozier et al., “Cytoscape: a software 2013. environment for integrated models of biomolecular interac- [6] Z. He, R. Fan, C. Zhang et al., “Chaihu-Shugan-San reinforces tion networks,” Genome Research, vol. 13, no. 11, CYP3A4 expression via pregnane X receptor in depressive pp. 2498–2504, 2003. treatment of liver-Qi stagnation syndrome,” Evidence-Based [21] A. Hamosh, A. F. Scott, J. S. Amberger, C. A. Bocchini, and Complementary and Alternative Medicine, vol. 201911 pages, V. A. McKusick, “Online Mendelian Inheritance in Man Article ID 9781675, 2019. (OMIM), a knowledgebase of human genes and genetic [7] K. K. Jia, S. M. Pan, H. Ding et al., “Chaihu-shugan san disorders,” Nucleic Acids Research, vol. 33, pp. D514–D517, inhibits inflammatory response to improve insulin signaling 2004. in liver and prefrontal cortex of CUMS rats with glucose [22] D. S. Wishart, Y. D. Feunang, A. C. Guo et al., “DrugBank 5.0: intolerance,” Biomedicine & Pharmacotherapy, vol. 103, a major update to the DrugBank database for 2018,” Nucleic pp. 1415–1428, 2018. Acids Research, vol. 46, no. 1, pp. D1074–D1082, 2018. [8] X. Huang, Y. Guo, W. H. Huang et al., “Searching the cy- [23] Y. Wang, S. Zhang, F. Li et al., “/erapeutic target database tochrome P450 enzymes for the metabolism of meranzin 2020: enriched resource for facilitating research and early hydrate: a prospective antidepressant originating from development of targeted therapeutics,” Nucleic Acids Re- Chaihu-Shugan-San,” PLoS One, vol. 9, no. 11, Article ID search, vol. 48, no. 1, pp. D1031–D1041, 2020. [24] G. Stelzer, N. Rosen, I. Plaschkes et al., “/e GeneCards suite: e113819, 2014. [9] Z. H. Su, H. M. Jia, H. W. Zhang, Y. F. Feng, L. An, and from gene data mining to disease genome sequence analyses,” Current Protocols in Bioinformatics, vol. 54, no. 1, 2016. Z. M. Zou, “Hippocampus and serum metabolomic studies to [25] Y. Tang, M. Li, J. Wang, Y. Pan, and F. X. Wu, “CytoNCA: a explore the regulation of Chaihu-Shu-Gan-San on metabolic cytoscape plugin for centrality analysis and evaluation of network disturbances of rats exposed to chronic variable protein interaction networks,” Biosystems, vol. 127, pp. 67–72, stress,” Molecular BioSystems, vol. 10, no. 3, pp. 549–561, 2014. 2015. [10] Y. Liang, Y. Zhang, Y. Deng et al., “Chaihu-Shugan-San [26] Y. Zhou, B. Zhou, L. Pache et al., “Metascape provides a decoction modulates intestinal microbe dysbiosis and alle- biologist-oriented resource for the analysis of systems-level viates chronic metabolic inflammation in NAFLD rats via the datasets,” Nature Communications, vol. 10, no. 1, p. 1523, NLRP3 inflammasome pathway,” Evidence-Based Comple- 2019. mentary and Alternative Medicine, vol. 201811 pages, Article [27] Z. H. Su, G. A. Zou, A. Preiss, H. W. Zhang, and Z. M. Zou, ID 9390786, 2018. “Online identification of the antioxidant constituents of [11] Q. H. Yang, Y. J. Xu, G. H. Feng et al., “p38 MAPK signal traditional Chinese medicine formula Chaihu-Shu-Gan-San pathway involved in anti-inflammatory effect of Chaihu- by LC-LTQ-Orbitrap mass spectrometry and microplate Shugan-San and Shen-ling-Bai-zhu-San on hepatocyte in spectrophotometer,” Journal of Pharmaceutical and Bio- non-alcoholic steatohepatitis rats,” African Journal of Tradi- medical Analysis, vol. 53, no. 3, pp. 454–461, 2010. tional, Complementary and Alternative Medicines, vol. 11, [28] X. Ni, M. Cao, Z. Wu, and L. Li, “Research progress on no. 1, pp. 213–221, 2013. chemical components and pharmacological effects of Chaihu [12] L. Liu, Y. Ni, and F. Yin, “Effect of Chaihu Shugan powder Shugan powder,” Shanghai Journal of Traditional Chinese combined with CAF in the treatment of stage III breast cancer Medicine, vol. 51, no. 9, pp. 109–113, 2017. and its effect on serum MMPs, estrogen level and tumor [29] Q. Liu, N. N. Sun, Z. Z. Wu, D. H. Fan, and M. Q. Cao, markers,” Journal of Chinese Medicinal Materials, vol. 43, “Chaihu-Shugan-San exerts an antidepressive effect by no. 1, pp. 220-221, 2020. downregulating miR-124 and releasing inhibition of the [13] F. Quan, B. Shi, and M. Yu, “Clinical effects of Chaihu Shugan MAPK14 and Gria3 signaling pathways,” Neural Regeneration powder combined with Paclitaxel and cyclophosphamide in Research, vol. 13, no. 5, pp. 837–845, 2018. treatment of breast cancer,” Grassroots Medical Forum, [30] H. Nie, Y. Deng, C. Zheng et al., “A network pharmacology- vol. 23, no. 26, pp. 3703-3704, 2019. based approach to explore the effects of Chaihu Shugan [14] W. He, “Advanced breast cancer treated with Bupleurum powder on a non-alcoholic fatty liver rat model through powder for soothing the liver combined with chemotherapy nuclear receptors,” Journal of Cellular and Molecular Medi- randomized parallel group study,” Journal of Practical Tradi- cine, vol. 24, no. 9, pp. 5168–5184, 2020. tional Chinese Internal Medicine, vol. 28, no. 8, pp. 83–86, 2014. [31] W. N. Jiang, D. Li, T. Jiang et al., “Protective effects of chaihu [15] S. LI and B. Zhang, “Traditional Chinese medicine network shugan san () on nonalcoholic fatty liver disease in rats with pharmacology: theory, methodology and application,” Chi- insulin resistance,” Chinese Journal of Integrative Medicine, nese Journal of Natural Medicines, vol. 11, no. 2, pp. 110–120, vol. 24, no. 2, pp. 125–132, 2018. 2014. [32] A. B. Awad, M. Chinnam, C. S. Fink, and P. G. Bradford, [16] A. L. Hopkins, “Network pharmacology: the next paradigm in “Beta-Sitosterol activates Fas signaling in human breast drug discovery,” Nature Chemical Biology, vol. 4, no. 11, cancer cells,” Phytomedicine, vol. 14, no. 11, pp. 747–754, pp. 682–690, 2008. 2007. Evidence-Based Complementary and Alternative Medicine 9

[33] M. S. B. Sayeed and S. S. Ameen, “Beta-sitosterol: a promising [48] Y. Jia, J. Zhou, X. Luo et al., “KLF4 overcomes tamoxifen but orphan nutraceutical to fight against cancer,” Nutrition resistance by suppressing MAPK signaling pathway and and Cancer, vol. 67, no. 8, pp. 1216–1222, 2015. predicts good prognosis in breast cancer,” Cellular Signalling, [34] M. A. Arowosegbe, O. T. Amusan, S. A. Adeola et al., vol. 42, pp. 165–175, 2018. “Kaempferol as a potential PAK4 inhibitor in triple negative [49] S. Park, S. H. Han, H. G. Kim et al., “PRPF4 is a novel breast cancer: extra precision glide docking and free energy therapeutic target for the treatment of breast cancer by calculation,” Current Drug Discovery Technologies, vol. 16, influencing growth, migration, invasion, and apoptosis of 2019. breast cancer cells via p38 MAPK signaling pathway,” Mo- [35] L. Zhu and L. Xue, “Kaempferol suppresses proliferation and lecular and Cellular Probes, vol. 47, Article ID 101440, 2019. induces cell cycle arrest, apoptosis, and DNA damage in breast [50] M. Hornsveld, L. M. Smits, M. Meerlo et al., “FOXO tran- cancer cells,” Oncology Research Featuring Preclinical and scription factors both suppress and support breast cancer Clinical Cancer 7erapeutics, vol. 27, no. 6, pp. 629–634, 2019. progression,” Progression. Cancer Research, vol. 78, no. 9, [36] S. Li, T. Yan, R. Deng et al., “Low dose of kaempferol sup- pp. 2356–2369, 2018. presses the migration and invasion of triple-negative breast [51] M. Bullock, “FOXO factors and breast cancer: outfoxing cancer cells by downregulating the activities of RhoA and endocrine resistance,” Endocrine-Related Cancer, vol. 23, Rac1,” Onco Targets and 7erapy, vol. 10, pp. 4809–4819, 2017. no. 2, pp. R113–R130, 2016. [37] Z. Zhao, G. Jin, Y. Ge, and Z. Guo, “Naringenin inhibits [52] K. U. Petry, “HPV and cervical cancer,” Scandinavian Journal migration of breast cancer cells via inflammatory and apo- of Clinical and Laboratory Investigation, vol. 74, pp. 59–62, 2014. ptosis cell signaling pathways,” Inflammopharmacology, [53] C. Y. Kan, B. J. Iacopetta, J. S. Lawson, and N. J. Whitaker, vol. 27, no. 5, pp. 1021–1036, 2019. “Identification of human papillomavirus DNA gene se- [38] J. Hu, Y. Zhang, X. Jiang et al., “ROS-mediated activation and quences in human breast cancer,” British Journal of Cancer, mitochondrial translocation of CaMKII contributes to Drp1- vol. 93, no. 8, pp. 946–948, 2005. dependent mitochondrial fission and apoptosis in triple- [54] N. Khodabandehlou, S. Mostafaei, A. Etemadi et al., “Human negative breast cancer cells by isorhamnetin and chloro- papilloma virus and breast cancer: the role of inflammation quine,” Journal of Experimental & Clinical Cancer Research, and viral expressed proteins,” BMC Cancer, vol. 19, no. 1, vol. 38, no. 1, p. 225, 2019. 61 pages, 2019. [39] L. Jia, S. Huang, X. Yin, Y. Zan, Y. Guo, and L. Han, [55] S. L. Glaser, J. L. Hsu, and M. L. Gulley, “Epstein-barr virus “Quercetin suppresses the mobility of breast cancer by sup- and breast cancer: state of the evidence for viral carcino- pressing glycolysis through Akt-mTOR pathway mediated genesis,” Cancer Epidemiology, Biomarkers & Prevention, autophagy induction,” Life Sciences, vol. 208, pp. 123–130, vol. 13, no. 5, pp. 688–697, 2004. 2018. [56] Y. Y. Liu, D. Hu, Q. Q. Fan et al., “Study on mechanism of [40] N. Zhang, H. Zhang, Y. Liu et al., “SREBP1, targeted by miR- chaihu shugan powder () for treating depression based on 18a-5p, modulates epithelial-mesenchymal transition in network pharmacology,” Chinese Journal of Integrative breast cancer via forming a co-repressor complex with Snail Medicine, 2019. and HDAC1/2,” Cell Death & Differentiation, vol. 26, no. 5, [57] Q. Zeng, L. Li, W. Siu et al., “A combined molecular biology pp. 843–859, 2019. and network pharmacology approach to investigate the multi- [41] Z. Tang, S. Ding, H. Huang et al., “HDAC1 triggers the target mechanisms of Chaihu Shugan San on Alzheimer’s proliferation and migration of breast cancer cells via upre- disease,” Biomedicine & Pharmacotherapy, vol. 120, Article ID gulation of interleukin-8,” Biological Chemistry, vol. 398, 109370, 2019. no. 12, pp. 1347–1356, 2017. [42] M. S. /ion, J. R. McGuire, C. M. Sousa et al., “Unraveling the role of Huntingtin in breast cancer metastasis,” Journal of the National Cancer Institute, vol. 107, no. 10, djv208 pages, 2015. [43] G. Li, X. Guo, M. Chen et al., “Prevalence and spectrum of AKT1, PIK3CA, PTEN and TP53 somatic mutations in Chinese breast cancer patients,” PLoS One, vol. 13, no. 9, Article ID e0203495, 2018. [44] H. Zhao, Z. Yu, L. Zhao et al., “HDAC2 overexpression is a poor prognostic factor of breast cancer patients with increased multidrug resistance-associated protein expression who re- ceived anthracyclines therapy,” Japanese Journal of Clinical Oncology, vol. 46, no. 10, pp. 893–902, 2016. [45] S. I. Hayashi and Y. Yamaguchi, “Estrogen signaling and prediction of endocrine therapy,” Cancer Chemotherapy & Pharmacology, vol. 56, no. 1, pp. 27–31, 2005. [46] A. J. Knox, M. J. Meegan, and D. G. Lloyd, “Estrogen re- ceptors: molecular interactions, virtual screening and future prospects,” Current Topics in Medicinal Chemistry, vol. 6, no. 3, pp. 217–243, 2006. [47] M. V. Dieci, F. Piacentini, M. Dominici et al., “Quantitative expression of estrogen receptor on relapse biopsy for ER- positive breast cancer: prognostic impact,” Anticancer re- search, vol. 34, no. 7, pp. 3657–3662, 2014.