# Exploration of the Molecular Mechanism of Danzhi Xiaoyao Powder in Endometrial Cancer through Network Pharmacology

Authors: Lanyu Li; Lukai Yang; Fang Liu; Jinfeng Qu
Journal: Evidence-Based Complementary and Alternative Medicine
 (2022)

---

## Abstract

Endometrial cancer (EC) is a common malignant tumor of the female reproductive system. Current treatments such as surgery and long-term hormone therapy are ineffective and have side effects. Danzhi Xiaoyao powder (DXP) can inhibit the growth of EC cells and induce apoptosis, but the pharmacological and molecular mechanisms of anticancer effects are still unclear. In this study, active components and potential targets of DXP were obtained from public databases. Protein effects and regulatory pathways of common targets were analyzed by protein-protein interaction (PPI), GO and KEGG. The results of network pharmacology showed that there are 87 common targets between EC and DXP. GO enrichment analysis showed that these targets were associated with response to oxidative stress, response to nutrient levels, hormone receptor binding and nuclear hormone receptor binding, etc. The results of KEGG analysis indicated that IL-17, TNF, PI3K/AKT, and RAS/RAF/MEK/ERK (ERK) signaling pathway were enriched in the anti-EC of DXP. Additionally, we cultured HEC-1B and KLE cells for validate experiments. DXP showed an inhibition of proliferation, migration, and cell cycle of both cells. Moreover, the expression of RAS, p-RAF, p-MEK, ERK, and p-ERK related proteins were downregulated. In conclusion, DXP might inhibit the proliferation of EC cells via apoptosis. Furthermore, DXP-induced inhibition of EC development might involve RAS/RAF/MEK/ERK pathway.


---

## Body

## 1. Introduction

Endometrial cancer (EC) is one of the malignant tumors of the female reproductive system, and its prevalence and mortality have rapidly risen worldwide [1, 2]. Globally, there were an estimated 417,000 incident cases and 97,000 deaths in 2020 [3]. The occurrence of EC is usually associated with overexposure to endogenous or exogenous estrogens [4]. Currently, the main treatment methods for EC are surgery and chemotherapy [5]. Chemotherapy drugs have certain toxicity, and due to the evolutionary nature of tumors, tumors cells may differentiate into subclonal cells with certain drug resistance, resulting in reduced antitumor efficacy of these chemotherapy drugs and possible adverse reactions [6, 7]. Therefore, it is important to find drugs with low side effects and better efficacy to improve the condition of EC.With the development of cellular biotechnology, breakthroughs have been made in the research of anticancer herbal medicines. In the past 20 years, natural antitumor drugs have accounted for 66.29% of the anticancer drugs on the market [8]. Traditional Chinese medicine (TCM) has unique advantages in terms of efficacy and safety. Its reasonable application can help control the development of the disease, reduce the toxic effect of radiotherapy and avoid the recurrence of cancer metastasis to a large extent [9]. Danzhi Xiaoyao powder (DXP) is a formula of TCM including Cortex Moutan, Gardeniae Fructus, Radix Bupleuri, Angelicae Sinensis Radix, Paeoniae Radix Alba, Atractylodes macrocephala Koidz., Poria Cocos (Schw.) Wolf., and Licorice [10]. Studies have found that DXP has antidepressant [11] effects, and some studies have found that DXP has some efficacy against breast cancer and precancerous lesions [12, 13], especially ER-positive breast cancer [13], suggesting that DXP may also have some anticancer effects. Clinical studies have shown that DXP has good therapeutic effects on patients with endometriosis liver depression and blood heat type [14]. Therefore, DXP has the potential to be a drug for the treatment of EC. However, its active components, possible targets, and molecular mechanism have not been elucidated.TCM plays a therapeutic role by restoring homeostasis in the biological process during disease development through multiple targets [15]. Network pharmacology is a new method to analyze the relationship between drugs, ingredients, targets, and diseases in the network based on biological networks and is gaining importance [16–18]. Therefore, the combining of TCM and network pharmacology is a new approach for the preferential selection of disease-related genes, revealing the association between diseases and drugs and further elucidating the regulatory role of TCM [19, 20].In this study, network pharmacology combined within vitro cellular experiment validation was performed to explore the active components, potential targets, and molecular mechanism of the anti-EC effect of DXP.


## 2. Materials and Methods

### 2.1. Collection of EC-Related Targets

We searched the Online Mendelian Inheritance in Man (OMIM,https://www.omim.org/) [21], DisGeNET (https://www.disgenet.org/) and GeneCards (https://www.genecards.org/) with the key word “endometrial cancer” to obtain EC-related targets. The duplicate disease targets were merged and deleted.


### 2.2. Collection of Active Components and Potential Targets of DXP

To collect all the relative components of DXP, we searched the traditional Chinese medicine systems pharmacology database and analysis platform (TCMSP,https://lsp.nwu.edu.cn/tcmsp.php) [22], Traditional Chinese Medicine Integrated Database (TCMID, https://tcm.cz3.nus.edu.sg/group/tcm-id/tcmid) [23], and Herbal Ingredients’ Targets Database (HIT, https://lifecenter.sgst.cn/hit) [24], which includes a great deal of information about Chinese herbal medicine. All predicted targets of active components of DXP were searched from TCM-ID, TCMSP, HIT, and Search Tool for Interacting Chemicals (STITCH, https://stitch.embl.de) databases and selected according to the component-target correlation degree above 400 from STITCH database.


#### 2.2.1. Drug-Likeness Screening

In the early stages of drug discovery, drug-likeness is a key consideration in the selection of compounds. A quantitative estimate of drug-likeness (QED), reported by Bickerton et al. in 2012, was used to evaluate a component’s drug-likeness [25]. In this work, DXP components with QED ≥ 0.3 were included for further study.


#### 2.2.2. Binomial Statistical Mode Analysis

According to previous studies, a binomial statistical mode was used for secondary screening for active components and potential targets of DXP [26, 27]. In this study, we set targets with P<0.05 that could be used as main targets for the formula.


### 2.3. Common Target Screening and Protein-Protein Interaction (PPI) Construction

The gene targets of components of DXP and EC were analyzed by Venny (version 2.1) to obtain the common genes, which were the targets of action of DXP for the treatment of EC.The PPI network was constructed by importing the common targets to the String database (https://string-db.org/cgi/input.pl), limiting the species to “Homo sapiens.” The results were imported into Cytoscape (version 3.8.0) for further analysis and optimization of the protein interaction network.


### 2.4. Construction of Herb-Components-Target (H-C-T) Network

An H-C-T network was constructed based on the active components of DXP and common targets. The above information was imported into Cytoscape for further analysis.


### 2.5. Enrichment Analysis

The clusterProfiler package of R (version 3.15.4) was applied for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. The GO analysis enriching for 3 aspects of potential targets is as follows: biological process (BP), cellular component (CC), and molecular function (MF).P<0.01 was set as a screening threshold and the functional terms and pathways with this value were considered statistically significant and retained.


### 2.6. Cell Culture

The EC cell line HEC-1B and KLE cells were purchased from the Shanghai Institute for Biological Sciences, Chinese Academy of Sciences. HEC-1B cells were cultured in Dulbecco’s Modified Eagle’s Medium (DMEM) medium (Gibco, MD, USA) with 10% fetal bovine serum (Gibco); the KLE cells were cultured in DMEM/F12 medium (Gibco) with 10% fetal bovine serum (Gibco), both at 37°C with 95% air and 5% CO2.


### 2.7. Detection of Cell Viability

HEC-1B and KLE cells (2 × 103/well) were seeded into 96-well plates and cultured overnight at 37°C with 95% air and 5% CO2. Then, cells were treated with different concentrations of DXP (0, 0.05, 0.1, 0.2, 0.4, 0.8, and 1.6 mg/mL) for 24 h and 48 h. Next, cells were incubated with CCK-8 solution (Biosharp, China) for 2 h at 37°C and the absorbance at 450 nm was measured by an automatic microplate reader (BioTek, VT, USA).


### 2.8. Flow Cytometry

HEC-1B and KLE cells were cultured in 6-well plates at 3 × 105 cells/well and treated with 0.4 or 0.8 mg/mL DXP for 48 h, respectively.For cell cycle analysis, cells were collected and washed with phosphate-buffered saline (PBS) followed by prechilled 70% ethanol at 4°C for more than 4 h. After fixation, cells were incubated with 100μg/mL RNase (Nanjing KeyGen Biotech Co., Ltd., China) for 30 min, followed by 50 μg/mL propidium iodide (PI) (Becton-Dickinson, NJ, USA) for 30 min in the dark. Finally, a flow cytometer (Becton-Dickinson) was used to detect the cell cycle of HEC-1B and KLE cells.For apoptosis analysis, cells were digested with trypsin without EDTA, washed twice with prechilled PBS, and then resuspended in the binding buffer. Next, cells were incubated with Annexin V-FITC (Becton-Dickinson) and PI for 15 min at room temperature under light-proof conditions. Finally, a flow cytometer (Becton-Dickinson) was used to detect the cell apoptosis of HEC-1B and KLE cells.


### 2.9. Cell Migration

The ability of migration in HEC-1B and KLE cells was detected by transwell assays. HEC-1B and KLE cells were resuspended in serum-free medium. 200μL·cell suspensions were added to the upper layer of the transwell chamber (Corning, NY, USA). The lower chamber contained 500 μL of serum-containing medium. After 48 h incubation, the cells in the lower chamber were stained with 0.1% crystal violet for 20 min and counted under a 400 × microscope (Olympus, Japan) in 5 random fields.


### 2.10. Western Blot

The cells were collected and lysed with cell lysis buffer (Beyotime, China) for 30 minutes on ice, then centrifuged at 12000 rpm and 4°C for 10 minutes. The supernatants were collected and the protein concentrations were quantified by the BCA kit (Pierce, IL, USA). The protein samples were separated by the 10% SDS-PAGE gradient gel. The separated protein was transferred onto PVDF membranes. The membranes were blocked with 5% skim milk for 90 min, and incubated with primary antibodies (1 : 1000, Abcam, UK) overnight at 4°C. Following this, they were incubated with secondary antibody (1 : 500,Hangzhou Multi Sciences Co., Ltd., China, GAR007) for 1 h. Finally, the protein bands were detected by an enhanced chemiluminescent (ECL) system (Beyotime). The primary antibodies were listed as below: RAS (ab52939), p-RAF (ab200653), MEK1 (ab32091), p-MEK1 (ab96379), ERK (ab32537), p-ERK (ab131438), and GAPDH (ab9485).


### 2.11. Statistical Analysis

The results were presented as means ± SD and data comparison of multiple groups was adopted by one-way ANOVA followed by Tukey’s test. All statistical analyses were obtained by GraghPad 7.0. The results were considered statistically significant if theP value was less than 0.05.


## 2.1. Collection of EC-Related Targets

We searched the Online Mendelian Inheritance in Man (OMIM,https://www.omim.org/) [21], DisGeNET (https://www.disgenet.org/) and GeneCards (https://www.genecards.org/) with the key word “endometrial cancer” to obtain EC-related targets. The duplicate disease targets were merged and deleted.


## 2.2. Collection of Active Components and Potential Targets of DXP

To collect all the relative components of DXP, we searched the traditional Chinese medicine systems pharmacology database and analysis platform (TCMSP,https://lsp.nwu.edu.cn/tcmsp.php) [22], Traditional Chinese Medicine Integrated Database (TCMID, https://tcm.cz3.nus.edu.sg/group/tcm-id/tcmid) [23], and Herbal Ingredients’ Targets Database (HIT, https://lifecenter.sgst.cn/hit) [24], which includes a great deal of information about Chinese herbal medicine. All predicted targets of active components of DXP were searched from TCM-ID, TCMSP, HIT, and Search Tool for Interacting Chemicals (STITCH, https://stitch.embl.de) databases and selected according to the component-target correlation degree above 400 from STITCH database.


### 2.2.1. Drug-Likeness Screening

In the early stages of drug discovery, drug-likeness is a key consideration in the selection of compounds. A quantitative estimate of drug-likeness (QED), reported by Bickerton et al. in 2012, was used to evaluate a component’s drug-likeness [25]. In this work, DXP components with QED ≥ 0.3 were included for further study.


### 2.2.2. Binomial Statistical Mode Analysis

According to previous studies, a binomial statistical mode was used for secondary screening for active components and potential targets of DXP [26, 27]. In this study, we set targets with P<0.05 that could be used as main targets for the formula.


## 2.2.1. Drug-Likeness Screening

In the early stages of drug discovery, drug-likeness is a key consideration in the selection of compounds. A quantitative estimate of drug-likeness (QED), reported by Bickerton et al. in 2012, was used to evaluate a component’s drug-likeness [25]. In this work, DXP components with QED ≥ 0.3 were included for further study.


## 2.2.2. Binomial Statistical Mode Analysis

According to previous studies, a binomial statistical mode was used for secondary screening for active components and potential targets of DXP [26, 27]. In this study, we set targets with P<0.05 that could be used as main targets for the formula.


## 2.3. Common Target Screening and Protein-Protein Interaction (PPI) Construction

The gene targets of components of DXP and EC were analyzed by Venny (version 2.1) to obtain the common genes, which were the targets of action of DXP for the treatment of EC.The PPI network was constructed by importing the common targets to the String database (https://string-db.org/cgi/input.pl), limiting the species to “Homo sapiens.” The results were imported into Cytoscape (version 3.8.0) for further analysis and optimization of the protein interaction network.


## 2.4. Construction of Herb-Components-Target (H-C-T) Network

An H-C-T network was constructed based on the active components of DXP and common targets. The above information was imported into Cytoscape for further analysis.


## 2.5. Enrichment Analysis

The clusterProfiler package of R (version 3.15.4) was applied for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. The GO analysis enriching for 3 aspects of potential targets is as follows: biological process (BP), cellular component (CC), and molecular function (MF).P<0.01 was set as a screening threshold and the functional terms and pathways with this value were considered statistically significant and retained.


## 2.6. Cell Culture

The EC cell line HEC-1B and KLE cells were purchased from the Shanghai Institute for Biological Sciences, Chinese Academy of Sciences. HEC-1B cells were cultured in Dulbecco’s Modified Eagle’s Medium (DMEM) medium (Gibco, MD, USA) with 10% fetal bovine serum (Gibco); the KLE cells were cultured in DMEM/F12 medium (Gibco) with 10% fetal bovine serum (Gibco), both at 37°C with 95% air and 5% CO2.


## 2.7. Detection of Cell Viability

HEC-1B and KLE cells (2 × 103/well) were seeded into 96-well plates and cultured overnight at 37°C with 95% air and 5% CO2. Then, cells were treated with different concentrations of DXP (0, 0.05, 0.1, 0.2, 0.4, 0.8, and 1.6 mg/mL) for 24 h and 48 h. Next, cells were incubated with CCK-8 solution (Biosharp, China) for 2 h at 37°C and the absorbance at 450 nm was measured by an automatic microplate reader (BioTek, VT, USA).


## 2.8. Flow Cytometry

HEC-1B and KLE cells were cultured in 6-well plates at 3 × 105 cells/well and treated with 0.4 or 0.8 mg/mL DXP for 48 h, respectively.For cell cycle analysis, cells were collected and washed with phosphate-buffered saline (PBS) followed by prechilled 70% ethanol at 4°C for more than 4 h. After fixation, cells were incubated with 100μg/mL RNase (Nanjing KeyGen Biotech Co., Ltd., China) for 30 min, followed by 50 μg/mL propidium iodide (PI) (Becton-Dickinson, NJ, USA) for 30 min in the dark. Finally, a flow cytometer (Becton-Dickinson) was used to detect the cell cycle of HEC-1B and KLE cells.For apoptosis analysis, cells were digested with trypsin without EDTA, washed twice with prechilled PBS, and then resuspended in the binding buffer. Next, cells were incubated with Annexin V-FITC (Becton-Dickinson) and PI for 15 min at room temperature under light-proof conditions. Finally, a flow cytometer (Becton-Dickinson) was used to detect the cell apoptosis of HEC-1B and KLE cells.


## 2.9. Cell Migration

The ability of migration in HEC-1B and KLE cells was detected by transwell assays. HEC-1B and KLE cells were resuspended in serum-free medium. 200μL·cell suspensions were added to the upper layer of the transwell chamber (Corning, NY, USA). The lower chamber contained 500 μL of serum-containing medium. After 48 h incubation, the cells in the lower chamber were stained with 0.1% crystal violet for 20 min and counted under a 400 × microscope (Olympus, Japan) in 5 random fields.


## 2.10. Western Blot

The cells were collected and lysed with cell lysis buffer (Beyotime, China) for 30 minutes on ice, then centrifuged at 12000 rpm and 4°C for 10 minutes. The supernatants were collected and the protein concentrations were quantified by the BCA kit (Pierce, IL, USA). The protein samples were separated by the 10% SDS-PAGE gradient gel. The separated protein was transferred onto PVDF membranes. The membranes were blocked with 5% skim milk for 90 min, and incubated with primary antibodies (1 : 1000, Abcam, UK) overnight at 4°C. Following this, they were incubated with secondary antibody (1 : 500,Hangzhou Multi Sciences Co., Ltd., China, GAR007) for 1 h. Finally, the protein bands were detected by an enhanced chemiluminescent (ECL) system (Beyotime). The primary antibodies were listed as below: RAS (ab52939), p-RAF (ab200653), MEK1 (ab32091), p-MEK1 (ab96379), ERK (ab32537), p-ERK (ab131438), and GAPDH (ab9485).


## 2.11. Statistical Analysis

The results were presented as means ± SD and data comparison of multiple groups was adopted by one-way ANOVA followed by Tukey’s test. All statistical analyses were obtained by GraghPad 7.0. The results were considered statistically significant if theP value was less than 0.05.


## 3. Results

### 3.1. EC-Related Targets

By retrieving three public databases, 836 EC-related targets were obtained, which included 237 in GeneCards, 197 in OMIM, and 402 in DisGeNET with the key word “endometrial cancer.” A total of 679 targets were screened by combining and deleting the repetitive values.


### 3.2. Active Components and Potential Targets of DXP

A total of 321 DXP-related components were searched from public databases including HIT, TCMSP, and TCMID databases. The potential targets of DXP were obtained from four public databases (HIT, TCMSP, TCMID, and STITCH), and there were 5,558 targets related to active components of DXP. Then, 274 active components were obtained according to drug-likeness screening (QED ≥ 0.3). Moreover, 589 main targets were obtained from binomial statistical mode analysis. Following it, we retrieved 214 active components of DXP. Among 589 screened components-related targets and 679 EC-related targets, 87 common targets were identified and considered as the key targets for further study (Figure1). Accordingly, 168 active components were screened and identified as the main active components of DXP that act on EC (Table 1).Figure 1
Venn diagram of common targets of danzhi xiaoyao powder (DXP) and endometrial cancer (EC).Table 1
Active components of DXP.
Chemical_nameQEDChemical_nameQEDBeta-elemene0.5799Phosphatidic acid0.4043Phenylalanine0.5664Glabrone0.8364Scopoletin0.5425Echinatin0.6610Guercetol0.5064Glabrol0.5836Quercetin0.5064Nonanoic acid0.5775Ursolic acid0.4433NON0.5672Calycosin0.8850(-)-Catechin0.5139Formononetin0.8529Gancaonin L0.4938Glabridin0.8323Hemo-sol0.4838Chrysin0.8206d-limonene0.4838Eugenol0.6955Alpha-limonene0.4838Thymol0.6510Stearic acid0.3017Kaempferol0.6372Catechin0.5139Farnesol0.6157WLN: QV40.5826isoliquiritigenin0.5824hexanoic acid0.5687Scopoletol0.5425Pentadecylic acid0.4059Capsaicin0.5370Myristicin0.6887IPH0.5172Scatole0.4871Puerarin0.4049Coumarin0.4124Glucuronic acid0.3250PAC0.6684Z-ligustilide0.6667Caprylic acid0.5818Linalool0.6172Octanoic acid0.5818Salicylic acid0.6129Oleanolic acid0.4460BUA0.5445Anethole0.6262L-Carvone0.5247Vanillin0.5173BuOH0.5171Gallocatechin0.4426Adenine0.5125WLN: VHR0.4967Crocetin0.5030Isorhamnetin0.6678Ethanol0.4049Geraniol0.6172IFP0.3920Mairin0.4635glutamate0.3835GLB0.3046Palmitic acid0.3653Histidine0.4184Esculetin0.3579EUG0.6993LPG0.3562HEX0.4628Alanine0.3562Vitexin0.3497Deoxygomisin a0.7145Methyleugenol0.6599Butylated hydroxytoluene0.7531EB0.5161Apocynin0.67362-Heptanone0.5103Paeonol0.5478Alpha-glycyrrhetinic acid0.5188Alpha-humulene0.4851Naphthalene0.5114Stigmasterol0.4599OYA0.3958Hyacinthin0.4290DFV0.8226Acetic acid0.4199Liquiritigenin0.8226Baicalin0.3617WLN: QVR BVQ0.6891Isovitexin0.34971,2-Benzenedicarboxylic acid0.6891Choline0.3405Licocoumarone0.5799Tetradecane0.3217Glabranin0.8174Caffeic acid0.4750PHPH0.5905Pinocembrin0.8226Genkwanin0.8862PHB0.61295,4′-Dihydroxy-7-methoxyflavone0.88625-Indolol0.5651PLO0.75023,4,5-Trihydroxybenzoic acid0.4656O-xylene0.4767Gallic acid0.4656Angelicin0.4670Umbelliferone0.3628Dodekan0.3809Morusin0.6040Dodecane0.3809DTY0.5110Epigallocatechin0.4426Genipin0.5093Beta-glycyrrhetinic acid0.5188Licochalcone A0.4563Liquiritin0.4838Beta-sitosterol0.4354Guasol0.5771Angelicin0.4354Guaiacol0.5771Zoomaric acid0.3256BU30.4701Uralenic acid0.6249PYG0.456018Beta-glycyrrhetinic acid0.5188Isoeugenol0.7273Glycyrrhetinic acid0.5188Ethyl protocatechuate0.549418Alpha-glycyrrhetic acid0.5188Safrol0.6221DBP0.4752Malic acid0.4307Myristic acid0.4490Fraxetin0.4141Fructose0.3101PCR0.5390DEP0.6925P-cresol0.5390Py0.4453O-thymol0.6510FER0.7180Carvacrol0.6510Isoimperatorin0.4856Licochalcone B0.4567Lauric acid0.3925Methylglyoxal0.3144Alpha-linolenic acid0.3326ICO0.6690Succinic acid0.5303Hexenal0.3885PHA0.5664Adonitol0.3082L-Ile0.4718Butal0.4435h-Met-h0.4506WLN: VH60.3957L-valine0.4266Nonanal0.3938Istidina0.4207Decanal0.3868L-valin0.4120PTL0.3745Prolinum0.3867Trans-2-nonenal0.3144Glutamine0.38352,4-Heptdienal0.3076Protocatechuic acid0.5265Prunetin0.8850


### 3.3. Construction of PPI Network

The 87 common targets of DXP and EC were imported to the String database for analysis, and Cytoscape (version 3.8.0) was used for PPI construction (Figure2). It contained 87 nodes and 1,587 edges, and the average degree value was 36.5.Figure 2
Common target protein-protein interaction (PPI) network between DXP and EC. Each node represents a target. The lines between the nodes represent the association among the different targets. The size and color of the node and its color shade change represent the size of the degree value.


### 3.4. Construction of Herb-Components-Target (H-C-T) Network

The H-C-T network was obtained through the import of 168 active components and 87 potential common targets of DXP into Cytoscape (version 3.8.0) and is shown in Figure3. There were 262 nodes and 1159 edges. The node size was proportional to its degree value. The red rhombus stood for different herbs which were contained in DXP. The cyan hexagon represented key targets of DXP acting on EC. Other different colored rhombuses represented the compounds contained in different herbs. Among them, the orange rhombus stood for the components that were contained in many kinds of herbs. The result showed that the top 5 active components nodes with the largest size of degree were quercetin, ethanol, kaempferol, IFP, and palmitic acid, with degree values of 64, 52, 43, 42, and 41, respectively. CASP3, AKT1, CYP1A1, JUN, and PPARG exhibited the highest degree values of 48, 25, 24, 22, and 22, respectively. These active compounds and targets may play a vital role in the effect of DXP against EC.Figure 3
Herb-components-target (H-C-T) network. The red rhombus represents the different herbs contained in DXP. The cyan hexagon represents the key targets of DXP acting on EC. Other different colors rhombuses represent the compounds contained in different herbs. Among them, the orange rhombus represents the components that are contained in many kinds of herbs.


### 3.5. Enrichment Analysis

In order to illuminate the biological functions and signaling pathways during the process of anti-EC of DXP, we employed the clusterProfiler package of the R language to perform GO and KEGG enrichment analysis on the 87 targets. A total of significant GO 1123 terms were obtained, including 1066 terms of BP, 18 terms of CC, and 39 terms of MF (P<0.01). The top 15 enriched terms of BP, CC, and MF were selected and visualized. The results showed that the targets were related to response to oxidative stress, response to nutrient levels, response to steroid hormone, and other biological processes; enriched in the transcription factor complex, membrane raft, membrane microdomain, and other cellular components; and related to RNA polymerase II transcription factor binding, hormone receptor binding, ubiquitin-like protein ligase binding, and other molecular functions (Figures 4(a)–4(c)).Figure 4
Enrichment analysis of GO and KEGG pathways of key target genes. The size of the circle represents the number of genes, and the color of the circle represents the correctedP value.
(a)(b)(c)(d)In KEGG enrichment analysis, a total of 143 pathways of anti-EC targets of DXP were identified according to theP value less than 0.01. The top 15 significantly enriched pathways with the highest gene counts were shown in Figure 4(d), and the specific information was shown in Table 2. The results showed that the targets may be associated with the AGE-RAGE signaling pathway in diabetic complications, the IL-17 signaling pathway, the TNF signaling pathway, and other signaling pathways. Among them, the endometrial cancer, apoptosis, EGFR tyrosine kinase inhibitor resistance, and proteoglycans in the cancer signaling pathway map showed that the PI3K-AKT and RAS-RAF-MEK-ERK (ERK) signaling pathways may play a crucial role in the treatment of EC with DXP (Figure 5).Table 2
Enrichment signaling pathways of KEGG.
IDDescriptionP.adjustCounthsa05215Prostate cancer6.24E − 3027hsa04933AGE-RAGE signaling pathway in diabetic complications3.62E − 2524hsa01522Endocrine resistance5.68E − 2423hsa05205Proteoglycans in cancer7.11E − 1925hsa05219Bladder cancer9.54E − 1915hsa01521EGFR tyrosine kinase inhibitor resistance1.86E − 1818hsa05212Pancreatic cancer2.61E − 1717hsa05225Hepatocellular carcinoma2.61E − 1722hsa05224Breast cancer2.61E − 1721hsa04657IL-17 signaling pathway3.59E − 1718hsa04917Prolactin signaling pathway1.39E − 1616hsa05213Endometrial cancer1.94E − 1615hsa05214Glioma4.10E − 1616hsa04668TNF signaling pathway7.33E − 1624hsa04210Apoptosis1.32E − 1518Figure 5
The map of endometrial cancer (a), apoptosis (b), EGFR tyrosine kinase inhibitor resistance (c), and proteoglycans in cancer (d) signaling pathway. The red boxes represent key targets.
(a)(b)(c)(d)The network of target-BP-pathways was obtained by Cytoscape (version 3.8.0) and shown in Figure6, which contained 105 nodes and 713 edges and included 75 common targets, the top 15 BP terms, and the top 15 signaling pathways. The square stood for disease targets, the triangle stood for signaling pathway, and the arrow represented biological processes.Figure 6
Core target-biological process-signal pathway network. The square represents DXP anti-EC targets, the triangle represents the top 15 signaling pathway, and the arrow represents the top 15 biological processes.


### 3.6. Correlation Analysis of Targets with BP and Signaling Pathways

To explore the correlation of the hub targets with the signaling pathway and BP more clearly, a hub target network was constituted with the top 10 targets, including TP53 (n = 71), AKT1 (n = 71), INS (n = 70), IL6 (n = 67), EGFR (n = 66), CASP3 (n = 65), ESR1 (n = 65), VEGFA (n = 65), MAPK3 (n = 64), and JUN (n = 63) (Figure 7(a)). The correlations were analyzed by combining them with the top 15 pathways and BP terms. The results showed that AKT1 had a higher correlation with the BP, and its response to oxidative stress was more relevant to the targets (Figure 7(b)). In the analysis of correlation with signaling pathways, AKT1 and MAPK3 were more important to the pathways (Figure 7(c)).Figure 7
Hub network of top 10 high-degree targets (a), the correlation between the core network and the first 15 BP terms (b), and the correlation between the core network and the first pathway (c).
(a)(b)(c)


### 3.7. DXP Suppresses the Viability, Cell Cycle, and Metastasis While Triggers the Apoptosis of EC Cells

Different concentrations of DXP were used to treat HEC-1B and KLE cells and assay their survival rates. As shown in Figure8(a), the cell viabilities of both types were decreased with the increasing concentration of the DXP at 24 h and 48 h. The LC50 values were calculated to be 0.717 mg/mL for HEC-1B and 0.603 mg/mL for KLE when DXP treated for 48 h. So, DXP (0.4 mg/mL and 0.8 mg/mL) treated HEC-1B and KLE cells for 48 h were selected in the subsequent experimental.Figure 8
HEC-1B and KLE cells were treated with different concentrations of DXP for cell proliferation (a), cycle (b), migration (c), and apoptosis (d) assays.∗P<0.05, ∗∗P<0.01 vs. Control group; ##P<0.01 vs. 0.4 mg/mL group.
(a)(b)(c)(d)Cell cycles results showed DXP significantly increased the proportion of G0/G1 phase cells (P<0.05), and DXP with a high-dose have a better inhibitory effect than low-dose (P<0.01) (Figure 8(b)). Meanwhile, DXP also inhibited the migration of EC cells. The apoptosis rate in high-dose was significantly higher than low-dose (P<0.01) (Figure 8(c)). Additionally, the apoptosis rate of EC cells was promoted by DXP addition (P<0.05) (Figure 8(d)). These results indicated that DXP suppressed the malignant progress of EC cells in vitro.


### 3.8. DXP Inhibited the Protein Expression of ERK Signaling Pathway

Network pharmacology analysis indicated that the RAS pathway is the most important pathway in the treatment of DXP on EC. Thus, the expression levels of RAS, p-RAF, MEK, p-MEK, ERK, and p-ERK were examined by Western blot. As shown in Figure9, DXP decreased the expression of RAS, p-RAF, p-MEK, and p-ERK proteins, but had no effect on MEK and ERK protein expression. It suggested that DXP against EC through ERK signaling pathway.Figure 9
The expression of RAS/RAF/MEK/ERK (ERK) pathway related proteins was detected using Western blot.∗P<0.05, ∗∗P<0.01 vs. Control group; #P<0.05, ##P<0.01 vs. 0.4 mg/mL group.


## 3.1. EC-Related Targets

By retrieving three public databases, 836 EC-related targets were obtained, which included 237 in GeneCards, 197 in OMIM, and 402 in DisGeNET with the key word “endometrial cancer.” A total of 679 targets were screened by combining and deleting the repetitive values.


## 3.2. Active Components and Potential Targets of DXP

A total of 321 DXP-related components were searched from public databases including HIT, TCMSP, and TCMID databases. The potential targets of DXP were obtained from four public databases (HIT, TCMSP, TCMID, and STITCH), and there were 5,558 targets related to active components of DXP. Then, 274 active components were obtained according to drug-likeness screening (QED ≥ 0.3). Moreover, 589 main targets were obtained from binomial statistical mode analysis. Following it, we retrieved 214 active components of DXP. Among 589 screened components-related targets and 679 EC-related targets, 87 common targets were identified and considered as the key targets for further study (Figure1). Accordingly, 168 active components were screened and identified as the main active components of DXP that act on EC (Table 1).Figure 1
Venn diagram of common targets of danzhi xiaoyao powder (DXP) and endometrial cancer (EC).Table 1
Active components of DXP.
Chemical_nameQEDChemical_nameQEDBeta-elemene0.5799Phosphatidic acid0.4043Phenylalanine0.5664Glabrone0.8364Scopoletin0.5425Echinatin0.6610Guercetol0.5064Glabrol0.5836Quercetin0.5064Nonanoic acid0.5775Ursolic acid0.4433NON0.5672Calycosin0.8850(-)-Catechin0.5139Formononetin0.8529Gancaonin L0.4938Glabridin0.8323Hemo-sol0.4838Chrysin0.8206d-limonene0.4838Eugenol0.6955Alpha-limonene0.4838Thymol0.6510Stearic acid0.3017Kaempferol0.6372Catechin0.5139Farnesol0.6157WLN: QV40.5826isoliquiritigenin0.5824hexanoic acid0.5687Scopoletol0.5425Pentadecylic acid0.4059Capsaicin0.5370Myristicin0.6887IPH0.5172Scatole0.4871Puerarin0.4049Coumarin0.4124Glucuronic acid0.3250PAC0.6684Z-ligustilide0.6667Caprylic acid0.5818Linalool0.6172Octanoic acid0.5818Salicylic acid0.6129Oleanolic acid0.4460BUA0.5445Anethole0.6262L-Carvone0.5247Vanillin0.5173BuOH0.5171Gallocatechin0.4426Adenine0.5125WLN: VHR0.4967Crocetin0.5030Isorhamnetin0.6678Ethanol0.4049Geraniol0.6172IFP0.3920Mairin0.4635glutamate0.3835GLB0.3046Palmitic acid0.3653Histidine0.4184Esculetin0.3579EUG0.6993LPG0.3562HEX0.4628Alanine0.3562Vitexin0.3497Deoxygomisin a0.7145Methyleugenol0.6599Butylated hydroxytoluene0.7531EB0.5161Apocynin0.67362-Heptanone0.5103Paeonol0.5478Alpha-glycyrrhetinic acid0.5188Alpha-humulene0.4851Naphthalene0.5114Stigmasterol0.4599OYA0.3958Hyacinthin0.4290DFV0.8226Acetic acid0.4199Liquiritigenin0.8226Baicalin0.3617WLN: QVR BVQ0.6891Isovitexin0.34971,2-Benzenedicarboxylic acid0.6891Choline0.3405Licocoumarone0.5799Tetradecane0.3217Glabranin0.8174Caffeic acid0.4750PHPH0.5905Pinocembrin0.8226Genkwanin0.8862PHB0.61295,4′-Dihydroxy-7-methoxyflavone0.88625-Indolol0.5651PLO0.75023,4,5-Trihydroxybenzoic acid0.4656O-xylene0.4767Gallic acid0.4656Angelicin0.4670Umbelliferone0.3628Dodekan0.3809Morusin0.6040Dodecane0.3809DTY0.5110Epigallocatechin0.4426Genipin0.5093Beta-glycyrrhetinic acid0.5188Licochalcone A0.4563Liquiritin0.4838Beta-sitosterol0.4354Guasol0.5771Angelicin0.4354Guaiacol0.5771Zoomaric acid0.3256BU30.4701Uralenic acid0.6249PYG0.456018Beta-glycyrrhetinic acid0.5188Isoeugenol0.7273Glycyrrhetinic acid0.5188Ethyl protocatechuate0.549418Alpha-glycyrrhetic acid0.5188Safrol0.6221DBP0.4752Malic acid0.4307Myristic acid0.4490Fraxetin0.4141Fructose0.3101PCR0.5390DEP0.6925P-cresol0.5390Py0.4453O-thymol0.6510FER0.7180Carvacrol0.6510Isoimperatorin0.4856Licochalcone B0.4567Lauric acid0.3925Methylglyoxal0.3144Alpha-linolenic acid0.3326ICO0.6690Succinic acid0.5303Hexenal0.3885PHA0.5664Adonitol0.3082L-Ile0.4718Butal0.4435h-Met-h0.4506WLN: VH60.3957L-valine0.4266Nonanal0.3938Istidina0.4207Decanal0.3868L-valin0.4120PTL0.3745Prolinum0.3867Trans-2-nonenal0.3144Glutamine0.38352,4-Heptdienal0.3076Protocatechuic acid0.5265Prunetin0.8850


## 3.3. Construction of PPI Network

The 87 common targets of DXP and EC were imported to the String database for analysis, and Cytoscape (version 3.8.0) was used for PPI construction (Figure2). It contained 87 nodes and 1,587 edges, and the average degree value was 36.5.Figure 2
Common target protein-protein interaction (PPI) network between DXP and EC. Each node represents a target. The lines between the nodes represent the association among the different targets. The size and color of the node and its color shade change represent the size of the degree value.


## 3.4. Construction of Herb-Components-Target (H-C-T) Network

The H-C-T network was obtained through the import of 168 active components and 87 potential common targets of DXP into Cytoscape (version 3.8.0) and is shown in Figure3. There were 262 nodes and 1159 edges. The node size was proportional to its degree value. The red rhombus stood for different herbs which were contained in DXP. The cyan hexagon represented key targets of DXP acting on EC. Other different colored rhombuses represented the compounds contained in different herbs. Among them, the orange rhombus stood for the components that were contained in many kinds of herbs. The result showed that the top 5 active components nodes with the largest size of degree were quercetin, ethanol, kaempferol, IFP, and palmitic acid, with degree values of 64, 52, 43, 42, and 41, respectively. CASP3, AKT1, CYP1A1, JUN, and PPARG exhibited the highest degree values of 48, 25, 24, 22, and 22, respectively. These active compounds and targets may play a vital role in the effect of DXP against EC.Figure 3
Herb-components-target (H-C-T) network. The red rhombus represents the different herbs contained in DXP. The cyan hexagon represents the key targets of DXP acting on EC. Other different colors rhombuses represent the compounds contained in different herbs. Among them, the orange rhombus represents the components that are contained in many kinds of herbs.


## 3.5. Enrichment Analysis

In order to illuminate the biological functions and signaling pathways during the process of anti-EC of DXP, we employed the clusterProfiler package of the R language to perform GO and KEGG enrichment analysis on the 87 targets. A total of significant GO 1123 terms were obtained, including 1066 terms of BP, 18 terms of CC, and 39 terms of MF (P<0.01). The top 15 enriched terms of BP, CC, and MF were selected and visualized. The results showed that the targets were related to response to oxidative stress, response to nutrient levels, response to steroid hormone, and other biological processes; enriched in the transcription factor complex, membrane raft, membrane microdomain, and other cellular components; and related to RNA polymerase II transcription factor binding, hormone receptor binding, ubiquitin-like protein ligase binding, and other molecular functions (Figures 4(a)–4(c)).Figure 4
Enrichment analysis of GO and KEGG pathways of key target genes. The size of the circle represents the number of genes, and the color of the circle represents the correctedP value.
(a)(b)(c)(d)In KEGG enrichment analysis, a total of 143 pathways of anti-EC targets of DXP were identified according to theP value less than 0.01. The top 15 significantly enriched pathways with the highest gene counts were shown in Figure 4(d), and the specific information was shown in Table 2. The results showed that the targets may be associated with the AGE-RAGE signaling pathway in diabetic complications, the IL-17 signaling pathway, the TNF signaling pathway, and other signaling pathways. Among them, the endometrial cancer, apoptosis, EGFR tyrosine kinase inhibitor resistance, and proteoglycans in the cancer signaling pathway map showed that the PI3K-AKT and RAS-RAF-MEK-ERK (ERK) signaling pathways may play a crucial role in the treatment of EC with DXP (Figure 5).Table 2
Enrichment signaling pathways of KEGG.
IDDescriptionP.adjustCounthsa05215Prostate cancer6.24E − 3027hsa04933AGE-RAGE signaling pathway in diabetic complications3.62E − 2524hsa01522Endocrine resistance5.68E − 2423hsa05205Proteoglycans in cancer7.11E − 1925hsa05219Bladder cancer9.54E − 1915hsa01521EGFR tyrosine kinase inhibitor resistance1.86E − 1818hsa05212Pancreatic cancer2.61E − 1717hsa05225Hepatocellular carcinoma2.61E − 1722hsa05224Breast cancer2.61E − 1721hsa04657IL-17 signaling pathway3.59E − 1718hsa04917Prolactin signaling pathway1.39E − 1616hsa05213Endometrial cancer1.94E − 1615hsa05214Glioma4.10E − 1616hsa04668TNF signaling pathway7.33E − 1624hsa04210Apoptosis1.32E − 1518Figure 5
The map of endometrial cancer (a), apoptosis (b), EGFR tyrosine kinase inhibitor resistance (c), and proteoglycans in cancer (d) signaling pathway. The red boxes represent key targets.
(a)(b)(c)(d)The network of target-BP-pathways was obtained by Cytoscape (version 3.8.0) and shown in Figure6, which contained 105 nodes and 713 edges and included 75 common targets, the top 15 BP terms, and the top 15 signaling pathways. The square stood for disease targets, the triangle stood for signaling pathway, and the arrow represented biological processes.Figure 6
Core target-biological process-signal pathway network. The square represents DXP anti-EC targets, the triangle represents the top 15 signaling pathway, and the arrow represents the top 15 biological processes.


## 3.6. Correlation Analysis of Targets with BP and Signaling Pathways

To explore the correlation of the hub targets with the signaling pathway and BP more clearly, a hub target network was constituted with the top 10 targets, including TP53 (n = 71), AKT1 (n = 71), INS (n = 70), IL6 (n = 67), EGFR (n = 66), CASP3 (n = 65), ESR1 (n = 65), VEGFA (n = 65), MAPK3 (n = 64), and JUN (n = 63) (Figure 7(a)). The correlations were analyzed by combining them with the top 15 pathways and BP terms. The results showed that AKT1 had a higher correlation with the BP, and its response to oxidative stress was more relevant to the targets (Figure 7(b)). In the analysis of correlation with signaling pathways, AKT1 and MAPK3 were more important to the pathways (Figure 7(c)).Figure 7
Hub network of top 10 high-degree targets (a), the correlation between the core network and the first 15 BP terms (b), and the correlation between the core network and the first pathway (c).
(a)(b)(c)


## 3.7. DXP Suppresses the Viability, Cell Cycle, and Metastasis While Triggers the Apoptosis of EC Cells

Different concentrations of DXP were used to treat HEC-1B and KLE cells and assay their survival rates. As shown in Figure8(a), the cell viabilities of both types were decreased with the increasing concentration of the DXP at 24 h and 48 h. The LC50 values were calculated to be 0.717 mg/mL for HEC-1B and 0.603 mg/mL for KLE when DXP treated for 48 h. So, DXP (0.4 mg/mL and 0.8 mg/mL) treated HEC-1B and KLE cells for 48 h were selected in the subsequent experimental.Figure 8
HEC-1B and KLE cells were treated with different concentrations of DXP for cell proliferation (a), cycle (b), migration (c), and apoptosis (d) assays.∗P<0.05, ∗∗P<0.01 vs. Control group; ##P<0.01 vs. 0.4 mg/mL group.
(a)(b)(c)(d)Cell cycles results showed DXP significantly increased the proportion of G0/G1 phase cells (P<0.05), and DXP with a high-dose have a better inhibitory effect than low-dose (P<0.01) (Figure 8(b)). Meanwhile, DXP also inhibited the migration of EC cells. The apoptosis rate in high-dose was significantly higher than low-dose (P<0.01) (Figure 8(c)). Additionally, the apoptosis rate of EC cells was promoted by DXP addition (P<0.05) (Figure 8(d)). These results indicated that DXP suppressed the malignant progress of EC cells in vitro.


## 3.8. DXP Inhibited the Protein Expression of ERK Signaling Pathway

Network pharmacology analysis indicated that the RAS pathway is the most important pathway in the treatment of DXP on EC. Thus, the expression levels of RAS, p-RAF, MEK, p-MEK, ERK, and p-ERK were examined by Western blot. As shown in Figure9, DXP decreased the expression of RAS, p-RAF, p-MEK, and p-ERK proteins, but had no effect on MEK and ERK protein expression. It suggested that DXP against EC through ERK signaling pathway.Figure 9
The expression of RAS/RAF/MEK/ERK (ERK) pathway related proteins was detected using Western blot.∗P<0.05, ∗∗P<0.01 vs. Control group; #P<0.05, ##P<0.01 vs. 0.4 mg/mL group.


## 4. Discussion

As an epithelial malignant tumor originating from the endometrium, EC ranks first in the incidence of malignant tumors of the female reproductive system, with a trend of younger incidence [28–30]. In recent years, the extraction of effective antitumor components with low toxicity from TCM has become a new research orientation [31, 32]. DXP is attracting attention for its anticancer and immune-enhancing activities [12]. In this work, we conducted network pharmacology and cellular experiments to identify the active components, core targets, and potential mechanisms of DXP in EC.In the present study, a total of 168 active components of DXP and 87 DXP-EC targets were screened from public databases. In the H-C-T network, we found that quercetin, associated with the largest number of targets, is the most important component in DXP for EC. Quercetin is a polyhydroxyflavonoid with various biological functions and has antioxidant, anti-inflammatory, and antiviral properties [33]. A New Jersey survey reported that lean women who consume quercetin have a lower risk of developing EC [34]. It has also been found that quercetin inhibits the proliferation of Ishikawa cells (an EC cell line) by downregulating EGF and cyclin D1 [35]. In addition, kaempferol also played a crucial role in the suppression of EC. Kaempferol is a dietary bioflavonoid with anticancer, anti-inflammatory, and antioxidant properties that inhibits human cancer cell proliferation through various mechanisms, including induction of the tumor suppressor p53 and inhibition of Erα [36]. Additionally, Lei’s study has substantiated that kaempferol against EC by inducing cell apoptosis, arresting the cell cycle in G2/M phase, suppressing cell invasion and upregulating the m-TOR/PI3K signaling pathway [37]. These reports combined with the present research provide evidence for the therapeutic potential of quercetin and kaempferol in DXP to cure EC.PPI analysis showed that the hub targets of DXP for EC were TP53, AKT1, INS, IL6, EGFR, CASP3, ESR1, VEGFA, MAPK3, and JUN. To annotate the functions and related pathway of 87 targets, we conducted the GO and KEGG signal pathway analysis. GO analysis results suggested that targets were mainly related to responses to oxidative stress, nutrient levels, steroid hormones, and other biological processes. The KEGG analysis results suggested that the pathways involved in DXP treating EC were mainly related to the IL-17 signaling pathway, TNF signaling pathway, PI3K-AKT signaling pathway, and ERK signaling pathway. Among them, IL-17 is a proinflammatory cytokine and is significantly expressed in a variety of cancers [38–40]. Lai’s study demonstrates that IL-17 may regulate chemokines and cytokines in gynecologic cancers through activation of IκBα and phosphorylation of extracellular signal-regulated kinase 1/2 [41]. PI3K-AKT signaling pathway is a key intracellular signaling pathway that regulates cell proliferation and cell cycle, and is closely related to tumorigenesis [42, 43]. Additionally, numerous studies have shown that the PI3K-AKT signaling pathway is involved in the malignant process of EC, such as picropodophyllin inhibiting EC cell proliferation by disrupting the PI3K-AKT signaling pathway [44]; NLRC5 promoting EC progression through activation of the PI3K/AKT pathway [45]. The ERK signaling pathway is part of the mitogen-activated protein kinase (MAPK) signaling pathway, which is involved in various cellular activities, such as cell proliferation, differentiation, apoptosis, inflammation, and other cellular activities, and is closely related to tumor-related immune escape mechanisms [46, 47]. Studies have reported that the ERK signaling pathway is one of the important pathways involved in endometrial carcinogenesis and is involved in tumor cell metabolism, growth, proliferation, survival, and angiogenesis [48].In vitro experiments were performed to further verify the molecular mechanism of DXP against EC. The cell experiments demonstrated that DXP for the treatment of EC suppressed cell proliferation, arrested the cell cycle in G0/G1 phase, induced cell apoptosis, weakened the ability of cell migration, downregulated RAS expression and RAF, MEK, and ERK phosphorylation levels. These findings verified the predicted result of the network pharmacology approach.


## 5. Conclusion

In summary, we combined network pharmacology andin vitro experiments to verify the potential mechanism of the bioactive components and targets of DXP in treating EC. We demonstrated that DXP could inhibit the proliferation and migration of EC cells and enhance EC cell apoptosis. Moreover, DXP protects against EC malignant progression via downregulating the ERK signaling pathway. This study provides an effective drug for the treatment of EC. Meanwhile, network pharmacology analysis reveals the molecular mechanisms involved in the treatment of EC with DXP, providing a research basis for subsequent studies.



---

*Source: 8330926-2022-06-13.xml*