1. Introduction
Colorectal cancer (CRC) is one of the deadliest diseases globally, ranking third in cancer morbidity and second in cancer mortality worldwide [
1]. Despite the benefits of early screening, surgery, and other therapeutic interventions, the five-year relative survival rate of CRC patients has not changed significantly in the past decades [
2]. Currently, treatment of colorectal cancer mainly includes surgery, radiotherapy, chemotherapy, and targeted therapies, such as checkpoint inhibitors and anti-angiogenesis therapies [
3]; nonetheless, the clinical outcomes are far from being satisfactory due to the considerable side effects and drug resistance. Therefore, there is an urgent need to search for effective and safe agents for CRC prevention and treatment.
Natural products have been shown to be an excellent and reliable source for the development of new drugs. The plant-derived small-molecule compounds, such as terpenoids, carotenoids, anthocyanidins, and flavonoids, have been extensively investigated for antitumor activities and proven to regulate gene expression and be involved in crucial biological processes such as cell proliferation, differentiation, apoptosis, and autophagy [
4]. 3
β-hydroxy-12-oleanen-27-oic acid (ATA), a cytotoxic oleanane triterpenoid, was isolated from the rhizome of
Astilbe chinensis (Maxim.) Franch. et Savat. [
5]. ATA has the same molecular weight and similar structures to the known antitumor pentacyclic triterpene compounds oleanolic acid (OA) and ursolic acid (UA), with the main difference being the position of the COOH and methyl groups (
Figure 1). ATA has been previously reported to exhibit distinctive cytotoxicity towards tumor cells, including human ovarian carcinoma HO-8910 cells, cervical squamous carcinoma HeLa cells, leukemic HL-60 cells, colorectal carcinoma COLO-205 cells, and human hepatoma HepG2 cells [
6,
7,
8], as well as to inhibit the growth of transplanted B16-F10 melanoma in mice [
8]. ATA was also reported to be considerably more cytotoxic towards HeLa and COLO-205 cells than its congener OA, implying its favorable structure and superior antitumor effect. In our previous experiment, ATA significantly inhibited the proliferation of a serial of human tumor cell lines, such as bladder carcinoma J82, T24, SW1710, and UM-UC-3, breast cancer MCF-7 and
T47D, colorectal cancer DLD-1 and HCT15, hepatoma HepG2, non-small cell lung cancer NCI-H460, ovary cancer Caov-3, prostate cancer LNCap (clone FGC), and renal carcinoma Caki-1, as well as normal human umbilical vein endothelial cells (HUVECs) and lung fibroblast MRC-5 cells via the CellTiter-Glo
® luminescent cell viability assay (CTG) (
Table S1). Notably, the cytotoxic effects of ATA were significantly stronger than its analogue compound OA to all fifteen cell lines. Among them, human colorectal cancer cells DLD-1 and HCT15 were indicated to be more sensitive to ATA than other cells. These results emphasized that ATA might be an excellent antitumor agent, especially for colorectal cancer, and that its unique structure might be responsible for its higher cytotoxic efficiency. Although ATA could inhibit the proliferation of neoplastic cells by inducing apoptosis, more investigations are needed to uncover the in-depth molecular mechanism of ATA to facilitate its development and applications.
High-throughput RNA sequencing (RNA-seq), a promising approach to profile gene expression, has been widely used to unravel molecular mechanisms and explore therapeutic targets of traditional Chinese medicine [
9]. Network pharmacology can establish a “compound-target-disease” network by integrating target prediction, statistical analysis network-based algorithms, and biological information analysis, as well as reveal the regulation principles of small molecules in combination with in vivo and vitro experimental validation [
10,
11]. In the present study, the cytotoxic effects of ATA were first confirmed against a serial of CRC cell lines, which revealed that HCT116 cells were the most sensitive to ATA. Meanwhile, ATA was proved to induce autophagy and ferroptosis in HCT116 cells. Furthermore, the underlying mechanisms of cytotoxic action and the potential targets of ATA were explored via integrative analysis of transcriptomics and network pharmacology in combination with in vitro and in vivo experimental validations.
3. Discussion
ATA has been previously proven to possess in vitro and in vivo antitumor activities and might be a promising antitumor agent with a unique chemical structure, high efficiency, and low toxicity [
5,
6,
7,
8]. In the pre-experiment, the inhibitory effects of ATA on the proliferation of 15 tested cell lines were found to not only be superior to its congener OA, as well as the two other currently widely researched and developed plant-derived compounds 23-hydroxybetulinic acid and
β-elemene, but equivalent to the positive control drugs cisplatin and oxaliplatin (
Table S1). In view of that, human colorectal cancer DLD1 and HCT15 cells seemed to be more sensitive to ATA than other cell lines among the tested tumor cell lines; thus, the cytotoxic effects of ATA against a serial of CRC cell lines were first detected. HCT116 cells were found to be the most sensitive to ATA (
Figure 2A). In this study, the antitumor effect of ATA against colorectal cancer and its molecular mechanism were explored via transcriptomic analysis and network pharmacology using HCT116 cells.
ATA significantly suppressed the proliferation of HCT116 cells in a concentration- and time-dependent manner (
Figure 2B). Meanwhile, ATA-treated HCT116 cells showed typical morphological changes of apoptosis, autophagy, and ferroptosis after specific fluorescence staining (
Figure 2C). FCM analysis also indicated that ATA induces cell cycle arrest at the G0/G1 phase (
Figure 2D,E) and apoptosis in HCT116 cells (
Figure 2F,G). These initial phenotypic validations of ATA against HCT116 cells were consistent with our previous results [
6,
7,
8], and suggested that ATA might inhibit cell proliferation through inducing cell cycle arrest and various types of cell death involving apoptosis, autophagy, and ferroptosis in CRC cells.
The level of transcription initiation is one of the most pivotal patterns in the regulation of gene expression. Transcriptomic analysis confirmed that ATA-induced DEGs in the HCT116 cells are predominately involved in the cell cycle, ferroptosis, and apoptosis, along with the Wnt and p53 signaling pathways, two well-established drivers of colon cancer (
Figure 5B). Since apoptosis resistance is a well-known characteristic of many cancer types, there is a high need for new therapeutic strategies [
16]. In this study, it was, for the first time, concerned that ATA induced cell autophagy and ferroptosis, two emerging types of non-apoptotic cell death currently playing important roles in colorectal cancer suppression [
17,
18], which was further evidenced via the involvement of the typical autophagic mTOR pathway and the positive regulatory PI3K-AKT pathway at the transcriptome level [
19]. Moreover, ATA induced endoplasmic reticulum (ER) stress in HCT116 cells, represented via the activation of UPR (
Figure 5D), which has been overwhelmingly evidenced to be a critical upstream event for apoptosis [
20], autophagy [
21], and ferroptosis [
22]. Among eleven hub genes, HSPA5, as one of the most important ER chaperones, mediates distinct ER stress signaling cascades [
23], while the heat shock proteins HSPA8 and HSP90AA1 are involved in protein refolding and chaperone-mediated autophagy, and HMOX1 takes part in the ferroptosis pathway (
Figure 5E). The above results suggested that ATA induces several promising types of cell death through ER stress as the upstream events.
To further identify the key target of ATA for its cytotoxic effects towards HCT116 cells, the integrating analysis of the network pharmacology computational prediction and disease targets from the transcriptomic data was conducted to reveal 15 potential targets of ATA (
Figure 6A). These 15 overlapping targets were involved in the oxygen level response, DNA damage, and mitotic cell cycle checkpoint (
Figure 6B), suggesting that ATA induces the cell cycle arrest and death of HCT116 cells through regulating these targets. Furthermore, the molecular docking results showed that ATA successfully docked with the three proteins FDFT1, PPARα, and PPARγ (
Figure 6C). However, only the FDFT1 inhibitor significantly blocked the cytotoxicity of ATA towards HCT116 cells (
Figure 6D–F). More notably, FDFT1 was exactly presented as one of the unique targets for ATA distinct from its two analogs UA and OA (
Figure 6G). Although having the same molecular formula and a similar structure, ATA exhibited higher levels of cytotoxic activity on the HCT116 cells than UA and OA (
Figure 6H and
Table S1). FDFT1 encodes farnesyl-diphosphate farnesyltransferase 1, a membrane-associated enzyme acting at a branch point in the mevalonate pathway and playing an important role in cholesterol biosynthesis [
24]. The genes responsible for cholesterol biosynthesis in the ATA-treated HCT116 cells were down-regulated and enriched in the critical module MCODE2 (
Figure 5C,D). In this study, ATA significantly down-regulated the mRNA and protein expression levels of FDFT1 in several CRC cell lines, including the HCT116, HT-29, HCT15, and LOVO cells (
Figure 7A–C and
Figure S5). Moreover, the inhibitory effects of ATA on the mRNA and protein expression levels of FDFT1 in four CRC cell lines were consistent with its cytotoxicity towards these cell lines (
Figure 1A). These findings suggested that FDFT1 acting at the beginning of the steroid biosynthesis pathway might be exactly the ATA-specific target biologically responsible for its distinctive antitumor mechanism. However, the inhibitory effects of ATA were not directly related to the basal expression level of FDFT in these cells (
Figure S5B), suggesting that FDFT1 might not be the only target of ATA against CRC cell lines.
Recently it has been demonstrated that disordered cholesterol homeostasis is a carcinogenic factor, thus leading to cancer progression [
25]. Changes in cholesterol biosynthesis, one of the important manifestations of lipid metabolism reprogramming, are regarded as a hallmark of a variety of cancers and are indispensable for cancer cell survival and proliferation [
26]. The abnormal expression of FDFT1 occurs in various cancers, which might be a novel candidate biomarker and a putative new target for cancer therapy [
27]. At present, there is little research on the role of FDFT1 in CRC, and the results are also controversial. The gene expression profiling interactive analysis (GEPIA) database (
http://gepia.cancer-pku.cn/, accessed on 11 February 2023) revealed that FDFT1 is highly expressed in colon cancer tissues (
Figure S6A) and that FDFT1 expression was associated with a better prognosis (
p = 0.018, log-rank test) (
Figure S6B) based on 275 tumor samples and 349 normal tissue samples from TCGA normal and GTEx data. Although it was also reported that high FDFT1 expression levels predicted better prognoses for patients with CRC in the Fudan University Shanghai Cancer Center (FUSCC) cohort (
p = 0.0238, log-rank test), the expression levels of the FDFT1 gene were significantly lower in CRC tissues than in normal tissues in the GDS2609 and GDS4382 datasets [
28]. Luo et al. reported that FDFT1 expression was higher in CRC tissues than normal tissues (
p < 0.05), and immunohistochemical results also confirmed that FDFT1 was highly expressed in CRC tissues. FDFT1 silencing obviously inhibited the tumor growth in both intraperitoneal implantation and subcutaneous xenograft models of CRC in C57BL/6 mice (
p < 0.05), but the colorectal peritoneal metastasis model of NCG mice presented no difference in tumor growth after silencing FDFT1 (
p > 0.05), suggesting that FDFT1 promoted the growth of CRC in vivo with dependence on the tumor immune microenvironment [
29].
Recently, FDFT1 has been identified as a ferroptosis-related gene and is considered an important gene for prognosis prediction in colorectal cancer patients [
30,
31]. Gathering evidence has indicated that FDFT1 expression increases in cells undergoing proliferation, suggesting that FDFT1 is implicated in proliferative signaling in cancer cells. In fact, FDFT1 regulates the cell cycle, where the inhibition of FDFT1 significantly impedes cells in the S-phase. FDFT1 could activate the NF-kB pathways, leading to an increase in the levels of anti-apoptotic proteins, such as Bcl-xL, Bcl-2, and Bax, and a decrease in the levels of pro-apoptotic proteins, such as caspase-3, thereby blocking apoptosis signaling [
32]. In addition, the high expression levels of FDFT1 increased the intracellular squalene levels, which protected the cell membrane from lipid peroxidation by reactive oxygen species (ROS) and further prevented the cells from entering the ferroptosis pathway [
33]. FDFT1 also participated in the biological process of ferroptosis through the Akt signaling pathway [
34], and the inhibition of FDFT1 could increase endogenous geranylgeranoic acid levels, resulting in incomplete autophagy, along with a reduction in cholesterol levels also inducing autophagy [
35]. The inhibition assay further confirmed that inhibition of FDFT1 could reverse the effect of ATA in the regulation of cellular proliferation and induction of autophagy and apoptosis (
Figure 7D,G–I). Meanwhile, FDFT1 knock-down significantly weakened the anti-proliferative effect of ATA on the HCT116 cells (
Figure 7M). More than that, ATA exhibited potent antitumor activity in the HCT116 xenograft mice (
Figure 8A–C). Meanwhile, the protein expression levels of FDFT1 was notably down-regulated, accompanied by the alteration of critical biomarkers of the cell cycle, autophagy, apoptosis, and ferroptosis in the tumor tissues from the ATA-treated mice (
Figure 8E,F). These results indicated that ATA exerted an antiproliferative effect against the HCT116 cells by targeting FDFT1 as a tumorigenesis gene in the present study (
Figure 9).
In this study, ATA significantly down-regulated the mRNA and protein expression levels of FDFT1 in HCT116 cells. Due to the lack of existing test kits, the impact of ATA on FDFT1 enzyme activity has not been tested. However, YM-53601, an inhibitor of FDFT1 enzyme activity [
36], significantly blocked the cytotoxicity of ATA against HCT116 cells at a non-cytotoxic concentration. The drugs with the same mechanism of action may exhibit competitive antagonism when combined due to varying intrinsic activities and affinities. From this fact, it was thereby speculated that ATA may also inhibit FDFT1 enzyme activity. Therefore, whether the specific mechanism is due to transcriptional blockade, protein degradation, and/or enzyme activity inhibition remains to be further explored.
This study demonstrated that ATA possessed the cytotoxic activities against HCT116 cells by inducing cell apoptosis, autophagy, and ferroptosis in vitro and in vivo. In mechanisms, ATA exerted the antiproliferative effect against HCT116 cells through targeting FDFT1. These data further expanded current knowledge on the mechanism of antitumor action of triterpenoid compounds and provided a new target for the research and design of anticancer agents.
4. Materials and Methods
4.1. Materials and Reagents
For this study, 3-(4,5-dimethyl thiazolyl-2)-2,5-diphenyltetrazolium bromide (MTT), acridine orange (AO), and dimethyl sulfoxide (DMSO) were purchased from Sigma-Aldrich, St. Louis, MO, USA; cell culture medium, trypsin, penicillin, streptomycin, and fetal bovine serum (FBS) were obtained from Gibco, Burlington, MA, USA; TRIzol was purchased from Invitrogen, Carlsbad, CA, USA; reverse transcriptase, oligo(dT)18, and ribonuclease inhibitor were obtained from Shanghai Sangon Biotech Co., Ltd., Shanghai, China; RIPA lysis buffer, BCA protein assay kit, anti-mouse actin mAbs, anti-caspase3 antibody (AF0081), horseradish peroxidase (HRP)-conjugated goat anti-rabbit and anti-mouse IgG (H + L), enhanced chemiluminescence (ECL) kit, and reactive oxygen species (ROS) assay kit were purchased from Beyotime, Shanghai, China; phosphatase inhibitor and protease inhibitor cocktails were obtained from Biotool, Houston, TX, USA; anti-FDFT1 antibody (ab195046) was purchased from Abcam, Cambridge, UK; anti-cyclinD1 antibody (ET1601-31) was obtained from HUABIO, Shanghai, China; anti-GPX4 (67763-1-Ig) antibody was obtained from Proteintech, Rosemont, IL, USA; the antibodies against LC3A/B (4108S) and p62 (39749S) were purchased from CST, Danvers, MA, USA; a monodansylcadaverine (MDC) sensor kit was obtained from KeyGEN BioTECH Corp., Ltd., Nanjing, China; an annexin V-FITC apoptosis detection kit and propidium iodide (PI)/RNase staining buffer were obtained from BD Pharmingen, San Diego, CA, USA; mRFP-GFP-LC3 lentivirus was obtained from Genechem, Shanghai, China.
3
β-Hydroxy-12-oleanen-27-oic acid (ATA, C
30H
48O
3, MW: 456.3594,
Figure 1) was previously isolated from the rhizomes of
A. chinensis [
5]. The purity of ATA was determined to be 98.9% through HPLC. Ursolic acid (UA, ≥ 98%), oleanolic acid (OA, ≥ 98%), and 23-hydroxybetulinic acid (23-HA, ≥ 98%) were purchased from Chengdu Lemeitian Pharmaceutical Technology Co., Ltd., Sichuan, China;
β-elemene (≥ 98%) was obtained from Shanghai Yuanye Bio-Technology Co., Ltd., Shanghai, China; cisplatin, oxaliplatin, PPARγ inhibitor GW9662, and PPARα inhibitor GW6471 were purchased from Selleck Chemicals, Houston, TX, USA; the FDFT1 inhibitor YM-53601 was obtained from MedChemExpress Technology, Monmouth Junction, NJ, USA. The stock solutions in DMSO were prepared and diluted as desired with the cell culture medium. The final concentration of DMSO in all experiments was less than 0.1% and did not show any detectable effect on cell growth or apoptosis.
4.2. Cell Lines and Culture
Murine colon carcinoma CT26, human bladder carcinoma J82, T24, SW1710, and UM-UC-3, breast cancer MCF-7 and T47D, colorectal cancer DLD-1, HCT15, HT29, Lovo, and HCT116, hepatoma HepG2, non-small cell lung cancer NCI-H460, ovary cancer Caov-3, prostate cancer LNCap clone FGC, and renal carcinoma Caki-1, as well as normal HUVEC, intestinal epithelial FHs74Int, and lung fibroblast MRC-5 cell lines were obtained from the National Collection of Authenticated Cell Cultures. Murine colorectal adenocarcinoma MC38 cell lines were purchased from MingzhouBio, Ningbo, Zhejiang, China. The cells were maintained in the logarithmic phase of growth in DMEM, MEM, McCoy’5a, or RPMI medium supplemented with 10% heat-inactivated FBS, 100 IU/mL penicillin, and 100 μg/mL streptomycin in an incubator at 37 °C with a humidified atmosphere of 5% CO2. Cell line identity was validated through short tandem repeat profiling, and routine mycoplasma testing was negative for contamination.
4.3. Experimental Animals
BALB/c-nu nude mice (male, five weeks of age) were purchased from Sino-British SIPPR/BK Laboratory Animal Co., Ltd. (Shanghai, China) and housed in a specific pathogen-free (SPF) environment at the Laboratory Animal Research Center of Zhejiang Chinese Medical University. All animal experiments were performed according to the procedures approved by the Institutional Animal Care and Use Committee of Zhejiang Chinese Medical University.
4.4. Cell Viability Assay
HCT116 cells were seeded at 1.0 × 10
4 cells/well in a 96-well plate and incubated at 37 °C in a humidified atmosphere with 5% CO
2. After 24 h, various concentrations of the tested samples were added into each well, and then cells were incubated for an indicated time. For the inhibition assay, HCT116 cells were pre-treated with YM-53601, GW9662, or GW6471 for 30 min, followed by co-incubation with the tested compounds for an additional 48 h. Each concentration was repeated for four wells. Four hours before the end, the extent of cell proliferation was detected using the MTT assay, as previously described [
37].
4.5. Cell Morphological Observation
The HCT116 cells were seeded at 1 × 10
5 cells/mL into a 24-well plate and incubated at 37 °C in a humidified atmosphere with 5% CO
2. After 24 h, the cells were either treated with ATA (10 μM and 20 μM) or MEM medium for 24 h. For cell apoptotic observation, cells were stained with AO (5 μM) in the culture medium at room temperature for 15 min. For cell autophagy detection, the cells were exposed to MDC (50 μM) at 37 °C for 30 min in the dark for staining using the KEYGEN MDC Sensor kit. The intracellular ROS levels were detected using the ROS assay kit [
37]. The cells were incubated with 10 μM 2′,7′-dichlorofluorescein diacetate (DCFH-DA) at 37 °C for 20 min. The stained cells were visualized via fluorescence microscopy (Axiover200, Zeiss, Oberkochen, Germany).
4.6. siRNA Transfection
FDFT1 siRNA and negative control siRNA were synthesized by GenePharma (Shanghai, China) and were transfected into HCT116 cells for 24 h or 48 h using Lipofectamine 3000 (Invitrogen, L3000015), according to the manufacturer’s protocol. The siRNA duplexes targeting FDFT1 indicated: siFDFT1-1: 5′-GUGCCUGAAUGAACUUAUATT-3′; siFDFT1-2: 5′-CCAUCUACCUGUCGUUUGUTT-3′.
4.7. Real-Time Quantitative Polymerase Chain Reaction (RT-qPCR)
The HCT116 cells were seeded at 1 × 10
5 cells/mL into a 24-well plate and incubated at 37 °C in a humid air with 5% CO
2. After 24 h, the cells were either treated with ATA (20 μM) or MEM medium for 3 h, 6 h, 12 h, 18 h, and 24 h. The total RNA was isolated with TRIzol reagent and reverse transcription was performed. The PCR was conducted on a Bio-Rad CFX 96 system using the SYBR Green qPCR Master Mix. The specific primers for RT-qPCR were synthesized by Sangon Biotech Co., Ltd. (Shanghai, China) and the sequences are listed in
Table S5. Primer amplification efficiency and specificity were verified for each set of primers. The mRNA expression levels of the tested genes relative to GAPDH were determined using the 2
-ΔΔCt method.
4.8. Flow Cytometry (FCM)
The HCT116 cells were seeded at 2 × 105 cells/well in 6-well plates and then cultured at 37 °C in a humidified atmosphere with 5% CO2 for 24 h. After treatment with ATA for 24 h, the cells were harvested and washed twice with PBS. For cell cycle detection, the collected cells were fixed in ice-cold 70% ethanol for overnight at −20 °C and then centrifugated. The pellet was resuspended in 0.5 mL of PI/RNase staining buffer and incubated for 20 min at room temperature in the dark. For the apoptosis assay, the cells were stained with the Annexin V-FITC/PI Apoptosis Kit according to the instructions. The analysis was performed on the flow cytometer (CytoFLEX, Beckman Coulter, CA, USA) using CytExpert 2.4 software.
4.9. RNA Sequencing (RNA-Seq)
The HCT116 cells treated with and without ATA (20 μM) for 3, 6, 12, 18, and 24 h were subjected to RNA-seq. The total RNA was isolated with TRIzol reagent. RNA quality was assessed using the NanoDrop UV–VIS spectrophotometer and Agilent 2100 Bioanalyzer. After digestion with DNase to remove the genomic DNA, mRNA purification and fragmentation were performed using the TruSeq Stranded mRNA LT Sample Prep Kit. The fragmented RNAs were used as templates to synthesize cDNA with reverse transcriptase. The double-stranded cDNA was purified and enriched using the Agencourt AMPure XP kit. Sequencing was conducted using an Illumina HiSeqTM2500 sequencer to afford paired-end data of 125 bp in length after constructing the library. Trimmomatic software was used to control the raw data quality, remove the linkers, and obtain high-quality clean reads [
38]. The obtained clean reads were matched to the human genome (GRCh38) using HISAT2to afford bam files [
39]. The fragments per kilobase of exon per million reads (FPKM) were quantified using cufflinks software (version 0.16) [
40,
41].
4.10. Differentially Expressed Genes (DEGs) Analysis
Significant difference analysis between samples was calculated using the DESeq R package. Differentially expressed genes (DEGs) were screened based on
p < 0.05 and fold change (FC) > 2 or < 0.5 compared with control samples calculated on the three replicates. The resulting significance scores were corrected for multiple testing using the Benjamini–Hochberg correction. Principal component analysis (PCA) was carried out to assess the variability and repeatability of samples using the normalized RNA-seq counts. A volcano plot was generated with the average FC and
p-value using online software (
http://sangerbox.com/, accessed on 10 March 2023). A Venn diagram was produced using online software (
http://bioinformatics.psb.ugent.be/webtools/Venn/, accessed on 13 March 2023). The heat maps were plotted using
https://www.omicstudio.cn/tool (accessed on 12 March 2023). The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were performed using Metascape software (
http://metascape.org/, accessed on 11 March 2023) [
42].
The protein–protein interaction (PPI) network of DEGs was constructed using the Search Tool for the Retrieval of Interacting Genes (STRING) online database (version 11.5) (
https://cn.string-db.org/, accessed on 12 March 2023), and then was visualized using Cytoscape software (v3.9.1) [
43,
44]. The plug-in molecular complex detection technology (MCODE) in Cytoscape software was used to analyze pivotal functional modules by the following criteria: K-core = 2, degree cutoff = 2, max depth = 100, and node score cutoff = 0.2 [
45]. The hub genes were identified using the cytoHubba plug-in of Cytoscape software with eight common algorithms (betweenness, bottleneck, MNC, degree, closeness, stress, EPC, and radiality) [
46]. The intersecting genes of the top twenty genes from all eight approaches of CytoHubba were considered as hub genes. The upstream transcription factors (TFs) of the hub genes were predicted using the Transcriptional Regulatory Relationships Unraveled by Sentence-based Text mining (TRRUST) database (
https://www.grnpedia.org/trrust/, accessed on 13 March 2023), and an adjusted
p-value < 0.05 was considered significant [
47].
4.11. Network Pharmacological Analysis
The structural information of ATA, OA, and UA were downloaded from the PubChem Database (
https://pubchem.ncbi.nlm.nih.gov/, accessed on 13 October 2022). The PubChem, Swiss target prediction (
http://www.swisstargetprediction.ch/, accessed on 13 October 2022), and SuperPred databases (
https://prediction.charite.de/s, accessed on 13 October 2022) were used to predict the potential targets of ATA, OA, and UA. DEGs and TFs regulating the hub genes from transcriptomic analysis were merged as disease targets of ATA. Then, the OmicStudio Venn tool (
https://www.omicstudio.cn/tool/6, accessed on 13 March 2023) was applied to produce Venn diagrams and derive the cross genes between the predicted targets and disease targets, which were considered as the candidate targets of ATA. Subsequently, a co-expression network of these key targets was constructed using GeneMANIA (
http://www.genemania.org/, accessed on 13 March 2023.) [
48]. Further, the unique genes of the predicted targets of the three compounds were respectively obtained using Venn diagrams.
4.12. Molecular Docking
The molecular docking simulation was performed using AutoDockTools (4.2.3) to verify the credibility of the hub targets, as previously described. The mechanical structure of ATA was optimized via Chem3D (21.0). The 3D structures of the proteins were obtained from the PDB database (
http://www.rcsb.org/pdb/home/home.do, accessed on 2 April 2023). The ligand–receptor binding property was analyzed. Binding energies < −7.0 kcal/mol denoted the docking, and the results were visualized using Pymol software (v2.5).
4.13. Western Blot Analysis
The HCT116 cells were seeded at 1 × 106 cells into a 6 cm dish and then incubated at 37 °C for 24 h in a humidified atmosphere with 5% CO2. After being treated with ATA for 24 h, the cells were washed twice with cold PBS. The tumor tissues (ca. 75 mg) or the collected cell pellets were lyzed with RIPA lysis buffer containing PMSF using the KZ-II high speed tissue grinder (Servicebio, Wuhan, China). The supernatants were collected through centrifugation at 12,000 rpm for 10 min at 4 °C. The protein concentrations in the supernatants were detected using the BCA protein assay kit. The denatured proteins were separated on a 10% SDS-polyacrylamide gel via electrophoresis and transferred to a polyvinylidene difluoride (PVDF) membrane. After blocking the membrane with 5% skim milk in Tris-buffered saline containing 0.1% Tween-20 (TBST) for 1 h at 37 °C, the blot was incubated with the primary antibodies overnight at 4 °C. Subsequently, the membranes were washed with TBST and incubated with HRP-conjugated goat anti-rabbit IgG (H + L) for 1 h. After washing the membrane with TBST for three times, the signal was visualized with the ECL kit on the iBright™ CL1500 Imaging System.
4.14. The mRFP-GFP-LC3 Assay
The mRFP and GFP expressed in mRFP-GFP-LC3 tandem fluorescent protein lentivirus were used to track LC3 [
49,
50]. The HCT116 cells were plated in confocal dishes and incubated overnight. After being transfected with mRFP-GFP-LC3 lentivirus for 48 h, the cells were selected with 2 μg/mL puromycin for 48 h. The cells stably expressing mRFP-GFP-LC3 were seeded into 6-well plates at a density of 2 × 10
5 cells in 2 mL of serum-containing medium. After 24 h, the cells were treated with ATA (20 μM) in the presence or absence of YM-53601 (4 μM), or with Torin (1 μM) for 24 h. The cells were visualized using fluorescence microscopy (Axiover200, Zeiss, Germany).
4.15. Xenograft Model in Nude Mice
HCT116 cells (1 × 107 cells in 200 μL of PBS) were injected subcutaneously into the flanks of the mice. When the tumor volume reached 80−100 mm3 on average, the mice were randomly divided into four groups, each consisting of four mice. The mice were orally administered with ATA at the doses of 40 mg/kg and 60 mg/kg once daily or injected i.p. with oxaliplatin (Oxa) at a dose of 6 mg/kg in saline twice a week for three weeks. The dose volume was 0.2 mL/10 g body weight. The model control (MC) groups received the same volume of saline. Body weights and tumor sizes were measured every two days. Tumor size was measured using a digital slide caliper and volumes (mm3) were calculated as follows: Tumor volume (mm3) = L × W2/2, where L is the length of the tumor and W is the width of the tumor. Animals were sacrificed after the final drug administration and tumors were carefully excised and weighed. The tumor tissues were snap frozen in liquid nitrogen and stored at −80 °C for subsequent analysis.
4.16. Statistical Analysis
Data were presented as mean ± SD and examined for their statistically significant differences with the analysis of variance (ANOVA) and Student’s t-test. The p-values of less than 0.05 were statistically significant. The calculations and graphs were performed using GraphPad Prism 9.0 software (GraphPad Software, San Diego, CA, USA).