Mechanism of Erianin anti-triple negative breast cancer based on transcriptomics methods and network pharmacology

Triple negative breast cancer (TNBC) is a highly aggressive illness that lacks effective targeted treatments. Although Erianin has shown potential antitumor properties, its precise mechanism of action and target in TNBC remain unclear, hampering the development of drugs. The present study investigated the underlying mechanism of action of Erianin in treating TNBC by using transcriptomics and network pharmacology approaches. We evaluated Erianin’s bioactivity in TNBC cell lines and xenograft tumor models. The results showed that Erianin significantly inhibited TNBC cell proliferation and impeded tumor growth. A subsequent analysis of transcriptomic and network pharmacological data identified 51 mutual targets. Analysis of protein-protein interactions identified eight hub targets. Furthermore, molecular docking indicated that the PPARA binding energy was the lowest for Erianin among the hub targets, followed by ROCK2, PDGFRB, CCND1, MUC1, and CDK1. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes functional enrichment analysis showed that the common targets were associated with multiple cancer-related signaling pathways, including focal adhesion, PI3K-Akt signaling pathway, Rap1 signaling pathway, microRNAs in cancer, and human papillomavirus infection. The results of the Western blot and immunohistochemistry experiment further showed that Erianin could suppress PI3K/Akt signaling pathway activation. After co-incubation with SC79, the cell inhibition rate of Erianin was decreased, which further confirmed that Erianin inhibits TNBC progression via the PI3K-AKT signaling pathway. In conclusion, our results indicated that Erianin has the potential to inhibit the proliferation of TNBC by downregulating the PI3K/AKT signaling pathway by transcriptomics and network pharmacology. Therefore, Erianin appears to be a promising compound for the effective treatment of TNBC.


INTRODUCTION
Breast cancer (BC) is the most commonly diagnosed cancer in women, posing a significant threat to women's health.The statistics on cancer in 2020 reveal that the incidence of BC is remarkably high, with 2.26 million new cases, exceeding the 2.2 million cases of lung cancer, thus making it the most prevalent cancer across the globe [1].In China, the yearly incidence and mortality rates of BC constitute 12.2% and 9.6% of the global sum, respectively [2].TNBC is a type of BC, comprising 10% to 20% of all invasive BCs [3].Because it does not express estrogen receptor (ER), progesterone receptor (PR) or human epidermal growth factor receptor 2 (HER2) and lacks clear therapeutic targets, endocrine therapy and anti-HER2 therapy are often not effective in TNBC patients.At the same time, TNBC exhibits the traits of heightened metastasis and recurrence [4][5][6].At present, the combination of taxane and anthracycline is the common choice for the clinical treatment of TNBC.However, anthracycline causes irreversible toxic damage to the heart, which is difficult for many patients to tolerate.Moreover, once patients develop chemotherapy drug resistance, the tumor will relapse and metastasize rapidly [7,8].Therefore, research and development of anticancer drugs with high efficiency and low toxicity has become a hot and difficult point in the treatment of TNBC.
Increasing attention is being paid to research on the study and application of substrates derived from natural plants.Erianin (2-methoxy-5-(2-(3,4,5-trimethoxy phenyl) ethyl)-phenol), a bibenzyl compound, is one of the active compounds extracted from Dendrobium.It has numerous biological functions, including programmed cell death induction, angiogenesis inhibition, and antioxidant and antitumor properties.Previous studies have shown that Erianin can inhibit the growth of tumor cells by inducing apoptosis [9,10], autophagy [11,12], ferroptosis [13] and other pathways.Sheng et al. [14] found that Erianin could exert its anti-liver cancer effect by inhibiting the activity of pyruvate carboxylase.Chen et al. [13] found that Erianin could inhibit the proliferation and migration of lung cancer cells via calcium/calmodulin-dependent ferroptosis.At present, little is known about the antitumor effects of Erianin on human BC cells in the BC setting, and it has been reported that Erianin can inhibit the proliferation and migration of T47D cells and can also induce cell apoptosis [15].Erianin and its derivatives (Ecust004) can suppress BC cell growth, invasion, and migration via EMT regulation [16].Erianin induces apoptosis in TNBC cells by inhibiting the PI3K/Akt pathway [17].However, the molecular mechanism of action and the drug targets of Erianin in TNBC are unclear, which limits the further development of this natural anticancer product with significant potential.
In recent years, network pharmacology has been used to study the pharmacological mechanisms of traditional Chinese medicine.It can mechanically associate drugs and diseases, quantitatively represent the key nodes of the network, including key molecules, key pathways or key modules, and predict the relationship between drugs and disease targets.In this study, we utilized transcriptomics and network pharmacology to explore the antiproliferative mechanism of Erianin in TNBC (Figure 1).First, we evaluated the tumor inhibitory Next, common targets, hub genes, biological functions and KEGG signaling pathways of Erianin against TNBC were analyzed by network pharmacology and transcriptomics.Furthermore, molecular docking analysis was used to predict the target genes of Erianin.Finally, the results were validated through biological experiments.Then, the PI3K-AKT signaling pathway was validated.

Cell culture
The human BC cell line MDA-MB-231 and mouse BC cell line 4T1 were purchased from the Peking Union Medical College Cell Culture Center (Beijing, China).MDA-MB-231 cells were cultured in Dulbecco's modified Eagle's medium (DMEM; C11995500BT, Thermo Fisher Scientific, USA) supplemented with 10% fetal bovine serum (FBS, 10099141, Thermo Fisher Scientific).The 4T1 cells were maintained in Roswell Park Memorial Institute-1640 medium (C11875500BT, Thermo Fisher Scientific) supplemented with 10% FBS.All cells were cultured in a humidified atmosphere of 5% CO2 at 37°C.

CCK-8 assay
The CCK-8 experiment was used to evaluate the effect of drug treatment on the proliferation and viability of MDA-MB-231 and 4T1 cells.The cells were inoculated into 48-well plates at a density of 2 × 10 4 cells per well.When the cells were 70% to 80% confluent, the cells were treated with the drug, and after a certain period of time, 20 µl of CCK-8 solution was added to each well.
Then, the optical density (OD) at 450 nm was measured by an enzyme meter.Cell survival rate (%) = (OD experimental group-OD blank group)/(OD negative control group-OD blank group) × 100%.

Wound healing assay
Cells were seeded (4 × 10 5 cells/well) in a 6-well plate, when the cell density was above 90%, some straight lines were drawn on the bottom of the 6-well plate with a sterile pipette tip.Discarded the old medium, washed it twice with preheated PBS, and added medium containing 1% serum.Then placed the 6-well plate in the incubator (37°C, 5% CO2) 24 hours.After the experiment, ImageJ software was used for data analysis.

In vivo xenograft assay
Six-to eight-week-old female Balb/c normal mice and Balb/c nude mice were purchased from Kunming Medical University.The mice were housed under a specific pathogen-free (SPF)-grade experimental system.Each Balb/c mouse was injected with 3 × 10 5 4T1 cells, and each Balb/c nude mouse was injected with 3 × 10 5 MDA-MB-231 cells in the second left mammary fat pad.When the tumor size reached approximately 15 mm 3 , all mice were randomly divided into two groups: the experimental group and the control group (five mice in each group).The mice in the experimental group were injected with 4 mg/kg Erianin (based on previous studies [9,12,[18][19][20][21][22][23] and our experimental validation), while the control group received an injection of the same volume of 1% DMSO.Every two days, the tumor volume was measured.The volume was calculated using the following formula: V = 0.5 × length × width 2 .Nude mice were executed after 21 days and Balb/c mice after 17 days and the tumours were photographed and recorded.All protocols were approved by the Animal Ethics Committee of Kunming University of Science and Technology (PZWH K2019-0005).

Transcriptome sequencing and analysis
Cells in the logarithmic growth phase were passaged.After 80% cell confluence, one dish was randomly selected as the experimental group and treated with 40 nM of Erianin.Similarly, another dish was randomly selected as the control group and treated with an equivalent volume of solvent (1‰ DMSO).The cells were then collected after 24 hours, and three samples were replicated three times.Subsequently, the cells were sent to Biomarker Technologies for transcriptome sequencing.Cell samples were sent to Bemac Biotechnology Ltd., for transcriptome sequencing.According to the method described by the company, RNA was initially extracted from the samples.RNA purity and concentration were assessed using a NanoDrop 2000 spectrophotometer.After accurate detection of RNA integrity by Agilent 2100/LabChip GX, library construction occurred upon passing the test.Library quality control was carried out employing the Qubit 3.0 fluorescence quantification instrument, Qsep400, and Q-PCR method.PCR methods were utilized for quality control of the library before sequencing using the Illumina NovaSeq 6000 sequencing platform.The differentially expressed genes (DEGs) were screened for expression differences using the condition | log2Fold Change | >1.

Drug and disease target acquisition
Exploring the genes associated with targeting a certain disease is an important foundation in modern small molecule drug discovery and development.To explore the genes associated with the anti-TNBC progression of Erianin, we analyzed the target genes of Erianin against TNBC by using web-based pharmacology.First, we obtained the 2D structure and 3D structure of Erianin from PubChem through BATMAN-TCM [24] (http://bionet.ncpsb.org.cn/batman-tcm/index.php), PharmMapper [25] (http://www.lilab-ecust.cn/pharmmapper/), Swiss Target Prediction (http://swisstarget prediction.ch/),TargetNet [26] (http://targetnet.scbdd.com/), PubChem (https://pubchem.ncbi.nlm.nih.gov/) and CTD [27] (https://ctdbase.org/)databases to obtain the targets of action of Erianin.Where the BATMAN-TCM database was obtained for Erianinrelated targets and screened for "Score cutoff >2.608" targets.Through the UniProt database, target names were standardized, and duplicates were removed and merged to create the Erianin target dataset."Triple negative breast cancer" or "TNBC" were used as the search terms, and the data were obtained from the GeneCards (https://www.genecards.org/),OMIM (https://www.omim.org/)and DisGenet [28] (https://www.disgenet.org/)databases.Among them, TNBC-related gene information was retrieved from the GeneCards database, and the genes with GeneCards GIFtS greater than 40 were screened.Subsequently, target names were standardized, and duplicates were removed and merged to establish the TNBC-related gene dataset through the UniProt database.

Construction of the protein-protein interaction and hub gene network
The online network analysis platform, STRING 11.5 [29] (http://string-db.org/) was used to analyze the network topology.The common targets were inputted into the STRING 11.5 database for creation, with the selection of "Homo sapiens" as the species.Network was then visualized and analyzed using Cytoscape 3.8.0.The Analyze Network and CytaHubba were then used to screen the hub genes based on degree in the PPI network.

Molecular docking analysis
Molecular docking can effectively predict the binding abilities between small molecule compounds and target genes.The 2D structure of Erianin was downloaded from the PubChem database.The 2D structure was then converted to a 3D structure in Chem3D software.The 3D crystal structures of 8 hub genes were retrieved from the PDB database.AutoDockTools 1.5.6 software was utilized to remove all ligands from the protein receptors and add hydrogens and charges before molecular docking simulation.Ligands and protein receptors were recorded in PDBQT format.Then, the binding sites of Erianin to the target gene receptor protein were examined using PyMOL 2.4.1.The binding energy was determined from the affinity.The higher the absolute affinity value, the stronger the binding affinity of Erianin with proteins.

Gene expression and survival analysis
The BC gene expression dataset, comprising tumor tissue and paraneoplastic normal tumor samples, was procured from the TCGA database, along with patientrelated clinical data.Patients were screened based on clinical information to obtain cases of TNBC with PRnegative, ER-negative, and Her-2-negative conditions.R software was utilized to analyze the expression matrices of both the normal and tumor groups.R software was utilized to extract expression data of potential target genes of Erianin acting on TNBC in the TCGA database.The edge R and DEGseq packages were employed for analyzing and screening DEGs in TNBC tissue samples and normal tissue samples.Histograms were plotted with loglFCl >1 and p < 0.05 thresholds.Additionally, we performed a survival analysis on possible target genes of Erianin in TNBC using the Kaplan-Meier Plotter web platform.We selected the best nodes while setting negative ER, PR, and Her-2 statuses.We regarded the genes with p < 0.05 in the results as influencing TNBC patient survival.

Gene function and pathway enrichment analysis
Gene Ontology (GO) [30] is internationally used to analyze the functional enrichment of genes and proteins, while Kyoto Encyclopedia of Genes and Genomes (KEGG) [31] pathway enrichment analysis provides a deeper explanation of different gene and protein functions.We used the R 4.3.1 programming language packages "clusterProfiler", "org.Hs.eg.db", "enrichplot" and "ggplot2" to analyze the GO and KEGG enrichment of target genes.The top 10 biological processes (BPs), cellular localizations (CCs) and molecular functions (MFs) in GO analysis and the top 15 pathways in KEGG analysis were screened at p < 0.05.

Western blotting (WB)
BC cells were cultured in 6-well plates at a density of 1 × 10 6 /ml per well.After treatment with Erianin for 24 hours, cells were lysed in precooled RIPA buffer.The protein concentrations were determined by a BCA protein assay kit (Beyotime, P0009).Equal amounts of protein were separated by SDS-PAGE and transferred to 0.45 μm PVDF membranes (Millipore, IPVH00010).After blocking with 2% BSA for 2 hours at room temperature, the membranes were incubated with primary antibodies at 4°C overnight.Then, the membranes were washed with TBST buffer 3 times and incubated with HRP-conjugated secondary antibody.The results were developed and recorded by a chemiluminescence analysis system (Tanon-5200, Tanon Science and Technology, China).

Immunohistochemistry (IHC)
Tumors were fixed for 24 hours in 4% paraformaldehyde and rinsed in running water.Then, paraffin-embedded tissue was cut into 5 μm thick sections.Sodium citrate buffer was used to boil the slides for 10 minutes to retrieve the antigens.After blocking with goat serum solution for 30 min, the sections were incubated overnight at 4°C with primary antibodies.A secondary antibody was added the next day, DAB chromogenic reagent was added, and the sections were counterstained with hematoxylin.Scores were assigned based on the intensity of the immunohistochemical signal and the area of the tissue sections stained positively.Finally, Image-Pro Plus 6.0 was used for quantitative analysis of the picture.

Statistical processing
All statistical data were analyzed using GraphPad Prism 8.0 software and R 4.3.1.The data were first tested for normal distribution, followed differences between the two groups were measured by T-test, and comparisons differences between the two groups were evaluated by one-way analysis of variance (ANOVA).The experimental data are expressed as the mean ± SD.Statistical significance was defined as p < 0.05.* p < 0.05, *** p < 0.001, **** p < 0.0001.

Erianin inhibits the progression of TNBC cells
To investigate the antiproliferation of Erianin in TNBC cells, MDA-MB-231 cells and 4T1 cells were treated with different concentrations of Erianin (0, 20, 40, 80, 120, 160 nM) for 24 hours, and the viability of the cells was measured by CCK-8 assay.The results indicated that Erianin had a significant inhibitory effect on the proliferation of MDA-MB-231 and 4T1 cells in a dosedependent manner (Figure 2A, 2B).Furthermore, we evaluated the effect of Erianin on cell migration.Considering the killing effect of Erianin on TNBC cells, we used 10 nM and 16 nM Erianin to treat MDA-MB-231 and 4T1 cells, respectively.In wound healing assay, Erianin significantly inhibited the migratory ability of MDA-MB-231 and 4T1 cells, compared with the control group (Figure 2C, 2D).In addition, two xenograft transplantation models were used to investigate the antitumor effect of Erianin in vivo.MDA-MB-231 cells and 4T1 cells were injected into Balb/c nude mice and Balb/c mice in situ.When the tumor size reached approximately 15 mm 3 , the mice were separated randomly into two experimental groups: the control group (DMSO) and the Erianin group.We found that 4 mg/kg Erianin treatment significantly inhibited the tumours growth of MDA-MB-231 and 4T1 cells (Figure 2E-2H).

Transcriptomic sequencing analysis of TNBC cells after Erianin treatment
By comparing the gene expression profile data of the Erianin group and Control group, we identified 1350 DEGs out of 13629 transcriptome annotated genes, following the threshold standard of |log2-fold change| > 1 and p < 0.05.The analysis revealed 1124 upregulated genes and 226 downregulated genes.The volcano plot and heatmap of all DEGs were shown in Figure 3A, 3B.To facilitate analysis, we excluded genes that were not expressed in any group from the transcriptome sequencing expression matrix.The top 10 upregulated and downregulated genes were shown in Figure 3C, 3D.

Network pharmacological screening of Erianin anti-TNBC target genes
We analyzed the target genes of Erianin against TNBC by using network pharmacology.First, we obtained the 2D structure and 3D structure of Erianin from PubChem (Figure 4A, 4B).Then, we obtained a total of 1002 Erianin-related targets from six drug databases and 5156 TNBC-related targets from the GeneCards, OMIM and DisGenet databases (Figure 4C), and the intersection of these two databases led to the identification of 622 Erianin targets in the TNBCrelated target genes (Figure 4D).Furthermore, we combined the potential target genes obtained from the network pharmacological analysis and the differential genes analyzed by transcriptome sequencing, and a total of 51 intersections were obtained (Figure 4E).We carried out subsequent analyses of these 51 common targets.

PPI network analysis of common targets
The clarity of drug targets is crucial for the in-depth development of drugs, and although our preliminary studies have shown that Erianin plays an important regulatory role in anti-TNBC progression, its target of action in TNBC has not yet been reported in the literature.Therefore, using the STRING database, we explored the interactions between the 51 common targets obtained from the above transcriptomics and network pharmacology to mine the interactions, selecting a confidence score >0.15.A complete interaction network was constructed using Cytoscape, and a total of 51 nodes and 314 edges were obtained.The darker the color, the higher the corresponding degree value.In the PPI network, we employed the degree method to identify the top 8 hub genes.These include CCND1, APOE, PDGFRB, PLAU, CDK1, PPARA, ROCK2 and MUC1 (Figure 5A).Most of these proteins had significant expression differences in TNBC patients (Figure 5B).For example, APOE, PLAU, CDK1 and MUC1 were significantly higher in TNBC tissues than in cancer-parasite tissues, whereas PDGFRB, RET, ROCK2, CCND1 and PPARA were significantly lower in TNBC tissues than in cancerparasite tissues (Figure 5C-5J).In addition, high expression of ROCK2 and PLAU in TNBC patients predicted poorer survival, whereas low expression of PDGFRB and CCND1 in TNBC patients predicted poorer survival (Figure 5K-5N).
To further clarify the specific targets of Erianin anti-TNBC, we performed molecular docking of these hub genes as target targets with Erianin, and the specific docking binding results were shown in Figure 6A.As shown in the table, Erianin has the lowest binding energy to PPARA, followed by ROCK2, PDGFRB, CCND1, MUC1 and CDK1.The lower docking binding energy indicates that Erianin has a higher binding affinity to these targets, suggesting that it has a greater possibility to act on these targets, affecting the structural changes of the targets and then regulating the corresponding signaling pathways.Docking structures of Erianin with the target proteins PPARA, ROCK2, PDGFRB, CCND1, MUC1 and CDK1 were shown in Figure 6B-6G.

Functional and pathway enrichment analyses of common targets
Furthermore, we analyzed GO with KEGG enrichment for 51 common targets and found that there were 2427 results related to BPs, of which 1007 results were significantly different (p < 0.05), 220 results related to CCs, of which 72 results were significantly different (p < 0.05), and 298 results related to MFs, of which 95 results were significantly different (p < 0.05).The GO results are shown in Figure 7A (TOP10).In terms of BP, these genes were enriched in positive regulation of ion transport, regulation of body fluid levels, response to oxygen levels, positive regulation of kinase activity and positive regulation of MAPK cascade.CCs alterations were mainly involved in neuronal cell bodies, cyclin-dependent protein kinase holoenzyme complexes, synaptic clefts, desmosomes and serine/ threonine protein kinase complexes.The MFs included peptide binding, amide binding, alpha-catenin binding, amyloid-beta binding and sulfur compound binding.The KEGG results showed that these potential target genes were enriched in a total of 206 pathways, of which 60 pathways were significantly different (p < 0.05), including focal adhesion, the PI3K-Akt signaling pathway, the Rap1 signaling pathway, microRNAs in cancer and human papillomavirus infection (Figure 7B).Among them, focal adhesion and the PI3K-AKT signaling pathway were significantly enriched.In addition, we conducted KEGG enrichment analysis on the differentially expressed genes (DEGs) obtained from transcriptome sequencing and found that the PI3K-AKT signaling pathway was significantly enriched among the DEGs (Figure 7C).The PI3K-AKT signaling pathway holds significance in responding to extracellular signals and regulating biological processes, including cell proliferation, apoptosis, and migration.Therefore, we explored the regulatory effects of Erianin on the PI3K-AKT signaling pathway by WB and immunohistochemical experiments.We treated two types of TNBC cells with concentrations of Erianin at 40 nM and 80 nM, respectively, which correspond to the half maximal inhibitory concentration for each cell type.Western blot results showed that Erianin significantly inhibited the expression of PI3K proteins and the phosphorylation of AKT after treating MDA-MB-231 and 4T1 cells for 24 hours (Figure 8A-8C).Furthermore, the expression of key proteins in the PI3K-AKT signaling pathway was also detected by immunohistochemistry in both transplantation tumor models, and the results were consistent with those of Western blot analysis (Figure 8D-8F).To further investigate whether Erianin can regulate the progression of TNBC through the PI3K-AKT signaling pathway, rescue experiments were conducted.We incubated TNBC cells with the Akt agonist SC79, CCK-8 assay revealed that 2 ug/ml SC79 treatment significantly reduced the inhibitory effect of Erianin on the proliferation of MDA-MB-231 and 4T1 cells, suggesting that the activation of the PI3K/ Akt pathway was associated with the antitumor effect of Erianin (Figure 8G, 8H).Western blot experiments also showed that SC79 significantly enhanced AKT phosphorylation and attenuated the inhibitory effect of Erianin on AKT phosphorylation (Figure 8I-8K).The above results indicated that Erianin may inhibit TNBC cell proliferation by affecting the PI3K-AKT pathway, which confirmed the findings of network pharmacology and transcriptome sequencing.

DISCUSSION
TNBC, as a complex heterogeneous disease, does not express ER, PR and HER2 and is characterized by early onset, high invasiveness, poor prognosis, early local recurrence and distant metastasis.TNBC has higher mortality and shorter median survival than non-TNBC.Surgery combined with chemotherapy and radiotherapy  is the main treatment for this kind of cancer, and effective targeted therapy strategies are urgently needed.
Compared with conventional chemotherapy, natural products have obvious advantages.They can interact with multiple targets and have high therapeutic value and low systemic toxicity.As a new anticancer drug, they have attracted widespread attention.Erianin is a small molecular dibenzyl compound extracted from Dendrobium.Previous studies have demonstrated that Erianin exerts a variety of pharmacological effects, including antiviral [32], antiaging, antidiabetic, and anti-inflammatory activities [33] and inhibits angiogenesis [34,35].Especially in the treatment of tumors, it has significant antitumor activity, including liver cancer [10,36], lung cancer [13,18,34], nasopharyngeal cancer [37], cervical cancer [38], bladder cancer [9], osteosarcoma [12], etc., However, the molecular mechanism and cellular targets through which Erianin acts against TNBC are poorly systematically defined, which is important for understanding its potential anticancer mechanisms.In this study, we elucidated the effects of Erianin on the proliferation of TNBC, verified the mechanism predicted by transcriptome analysis and network pharmacology, and explored the potential target of Erianin both in vitro and in vivo.
First, MDA-MB-231 and 4T1 cells were selected to determine the anti-proliferative activity of Erianin by CKK-8.The results showed that Erianin significantly inhibited the proliferation of TNBC cells in a dose-dependent manner.Meanwhile, Erianin could significantly inhibit the migratory ability of MDA-MB-231 and 4T1 cells, respectively.Further in vivo experiments also showed that Erianin treatment notably suppressed the growth of transplanted tumors in mice, indicating that Erianin had strong anti-TNBC proliferation activity.Next, 1385 DEGs were obtained by transcriptome sequencing analysis, and volcano and heatmaps were drawn.Meanwhile, to explore the genes related to Erianin anti-TNBC progression, we analyzed the target genes of Erianin anti-TNBC using network pharmacology and identified 622 genes shared by Erianin and TNBC.Furthermore, we combined the potential target genes obtained from the network pharmacology analysis and the differential genes analyzed by transcriptome sequencing, and a total of 51 intersections were obtained.To investigate the molecular mechanism and potential targets of Erianin in inhibiting the proliferation of TNBC.We constructed a PPI network on these 51 genes.PPI network analysis screened eight central targets, including CCND1, APOE, PDGFRB, PLAU, CDK1, PPARA, ROCK2 and MUC1.In previous literature, most of these genes have been associated with the onset and progression of TNBC.For example, CCND1 played a significant role in TNBC proliferation, migration, and treatment [39][40][41].PDGFRB was found to be significantly upregulated in aggressive TNBC tumor cells and the tumor microenvironment.Inhibition of its expression has been shown to effectively reduce tumor growth and lung metastasis [42,43].The cell cycle protein CDK1 was overexpressed in TNBC and may serve as a novel biomarker [41,44] and therapeutic target for the diagnosis and treatment of TNBC [45].The PPARA pathway exhibited downregulation in TNBC.Targeting this pathway could serve as a therapeutic strategy for treating TNBC [46].MUC1 was expressed in TNBC cell lines as a metabolic regulator in TNBC and promotes metabolic reprogramming of glutamine utilization that affects TNBC tumor growth [47].Targeting tumor MUC1 glycoprotein inhibited the growth of TNBC [48], among others.Meanwhile, these proteins have significant expression differences in TNBC patients; for example, APOE, PLAU, CDK1 and MUC1 were significantly higher in TNBC patients than in paraneoplastic tissues, while PDGFRB, RET, ROCK2, CCND1 and PPARA were significantly lower in TNBC patients than in paraneoplastic tissues.In addition, high expression of ROCK2 and PLAU in TNBC patients predicted poorer survival, while low expression of PDGFRB and CCND1 in TNBC patients predicted poorer survival.
To specify the potential targets of Erianin for TNBC, we performed docking of these hub genes with Erianin.Our results indicated that PPARA had the lowest binding energy to Erianin, followed by ROCK2, PDGFRB, CCND1, MUC1, and CDK1.The lower docking binding energy implies that Erianin exhibits higher binding affinity toward these targets, indicating its greater potential to act upon and affect the target structure, ultimately regulating the corresponding signaling pathways.The PPARA pathwaywas suppressed in breast cancer expression, specifically in subtypes of TNBC [46,49].It played a vital role in the invasion and metastasis of TNBC [50,51].PPAR-α can also serve as a target for TNBC treatment [52,53].A case in point was fenofibrate, a PPAR-α agonist that can induce apoptosis in TNBC cells by activating the NF-κB pathway [53].
Further KEGG enrichment analyses showed that the 51 potential targets may exert antitumor effects by modulating multiple signaling pathways, such as focal adhesion, the PI3K-Akt signaling pathway, the Rap1 signaling pathway, microRNAs in cancer and human papillomavirus infection.Notably, focal adhesion and the PI3K-Akt signaling pathway were significantly enriched.Meanwhile, transcriptome analysis results suggested that the PI3K-Akt signaling pathway may be crucial for cell death induced by Erianin, indicating the significant role of the PI3K-AKT signaling pathway in Erianin anti-TNBC progression.The PI3K/AKT signaling pathway is a crucial pathway that reacts to external signals and facilitates metabolism, proliferation, cell viability, enlargement and angiogenesis.It is considered a classic pathway related to various cancers [54].Previous studies have shown that Erianin can inhibit the proliferation of BC [17], gastric cancer [55], lung cancer [37], liver cancer [10] and other tumor cells by inhibiting the PI3K/Akt signaling pathway.Therefore, we detected the activity of the PI3K/AKT signaling pathway after Erianin treatment in TNBC.Consistent with the above studies, the expression levels of PI3K and p-AKT were significantly decreased following treatment with Erianin both in vivo and in vitro.In the rescue experiment, the inhibitory effect of Erianin on AKT phosphorylation and TNBC proliferation was significantly reduced by SC79, an AKT phosphorylation activator.Thus, the PI3K/AKT signaling pathway may play an important role in the growth inhibition of TNBC cells and tissues induced by Erianin.
In general, our study utilized network pharmacology, molecular docking, transcriptome analysis, and in vitro and in vivo experiments to comprehensively investigate the mechanism by which Erianin can treat TNBC.Our research indicates that Erianin can potentially serve as a viable medication for treating TNBC.This is due to its ability to restrain the growth of TNBC by inhibiting the PI3K/AKT signaling pathway.Meanwhile, it offered a new viewpoint for examining the function of natural small molecular compounds in the therapy of malignant tumors.

Figure 1 .
Figure 1.Workflow diagram of the research.

Figure 2 .
Figure 2. Erianin inhibits the progression of TNBC cells in vitro and in vivo.MDA-MB-231 (A) and 4T1 cells (B) were treated with different concentrations of Erianin for 24 hours, and then the viability was determined using the CCK-8 assay.(C) The migration capacity of MDA-MB-231 and 4T1 cell with Erianin was detected by wound-healing assay.(D) Statistical analysis of wound-healing assay.(E) MDA-MB-231 cells were injected into Balb/c nude mice in situ, and xenograft tumors were imaged after Erianin treatment.(F) Tumor volume in each group (N = 5).(G) 4T1 cells were injected into Balb/c normal mice in situ, and xenograft tumors were imaged after Erianin treatment.(H) Tumor volume in each group (N = 5).

Figure 3 .
Figure 3. RNA-seq data analysis.(A) A volcano plot of 1350 DEGs, including 226 downregulated genes and 1124 upregulated genes in response to Erianin treatment, was drawn.(B) Heatmap of 1350 DEGs was created.(C) Heatmap of the top 10 upregulated or downregulated genes based on the ranking of expression changes.(D) Table of the top 10 upregulated or downregulated genes based on the ranking of expression changes.

Figure 4 .
Figure 4. Screening of target genes for anti-TNBC activity of Erianin.(A) 2D structure of the Erianin molecule.(B) 3D structure of the Erianin molecule; (C) Flowchart of the screening process of target genes of Erianin anti-TNBC obtained by network pharmacological analysis.(D) Venn diagram of target genes of Erianin in anti-TNBC obtained by network pharmacological analysis; (E) Venn diagram of target genes identified in network pharmacology analysis and differentially expressed genes (DEGs).

Figure 5 .
Figure 5. Bioinformatic analysis of PPI-enriched and hub genes.(A) PPI of potential hub target genes.(B) The expression of hub genes in TNBC and the correlation analysis with the overall survival time of patients.(C-J) Expression maps of MUC1, APOE, PLAU, CDK1, PDGFRB, CCND1, ROCK2, and PPARA in TNBC patients and normal tissues.(K-N) Impact of ROCK2, PLAU, PDGFRB, and CCND1 expression on the survival of TNBC patients.

Figure 6 .
Figure 6.Molecular docking results of Hub protein and Erianin.(A) Table of molecular docking results of Erianin and hub genes.(B-G) The molecular docking visualization results of Erianin binding to PPARA, ROCK2, PDGFRB, CCND1, MUC1 and CDK1.

Figure 7 .
Figure 7.The results of GO and KEGG pathway enrichment.(A) Top 10 results of GO for target genes of Erianin against TNBC.(B) Top 15 enriched KEGG pathways of target genes of Erianin against TNBC.(C) KEGG results for DEGs obtained by transcriptome sequencing.

Figure 8 .
Figure 8. Erianin inhibits the PI3K-AKT signaling pathway in TNBC.(A) Western blot detection of PI3K and AKT expression as well as AKT phosphorylation after Erianin treatment of MDA-MB-231 and 4T1 cells.(B) Quantitative graphs of PI3K expression after Erianin treatment of MDA-MB-231 and 4T1 cells.(C) Quantitative graphs of p-AKT expression after Erianin treatment of MDA-MB-231 and 4T1 cells.(D) Immunohistochemistry detection of PI3K and AKT expression as well as AKT phosphorylation in two transplantation tumor models.(E) Quantitative graphs of PI3K expression after Erianin treatment of two transplantation tumor models.(F) Immunohistochemical detection of PI3K and AKT expression as well as AKT phosphorylation in two transplantation tumor models.(G) MDA-MB-21 cells were treated with 40 nM of Erianin and 2 ug/ml of SC79 for 24 hours, and then the viability was determined using the CCK-8 assay.(H) 4T1 cells were treated with 80 nM of Erianin and 2 ug/ml of SC79 for 24 hours, and then the viability was determined using the CCK-8 assay.(I) Western blot detection of AKT expression and AKT phosphorylation after Erianin and SC79 treatment of MDA-MB-231 and 4T1 cells.(J) Quantitative graphs of PI3K expression in MDA-MB-231 and 4T1 cells.(K) Quantitative graphs of AKT phosphorylation expression in MDA-MB-231 and 4T1 cells.