Abstract Purpose Bile acids promote the progression of colon adenocarcinoma (COAD), and ursodeoxycholic acid (UDCA) is a key drug in promoting bile acid excretion, but its role in COAD unclear. Our study aims to investigate the relationship between COAD and bile acid metabolism and to assess the feasibility of UDCA for the treatment of COAD. Methods Firstly, biological targets closely related to COAD were identified: Based on the cancer genome atlas (TCGA) database, the core genes of COAD were obtained by differential expression analysis and weighted gene-coexpression network analysis (WGCNA), and subjected to gene ontology (GO) function annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. Secondly, finding a drug by target, after identifying UDCA as a candidate drug, the feasibility of UDCA in treating COAD was verified in reverse: Using databases to collect potential targets for COAD and UDCA, and the intersecting genes were the potential targets for UDCA to exert anti-tumor effects. Then Autodock was used for molecular docking to analyze the interaction between UDCA and core target proteins. Finally, experimental validation was performed: MTT assay, wound healing, transwell migration, and angiogenesis assays were used to detect the effects of UDCA on cell proliferation, migration, invasion, and neovascularization. Results 2064 differential genes were screened from TCGA. WGCNA obtained the module most relevant to CRC, containing 493 genes. KEGG analysis found that overlapping genes were mainly concentrated in bile acid metabolic pathways. A total of 26 UDCA anti-tumor targets were obtained in database, and 5 core targets were selected by STRING database and Cytoscape software: TNF, CYP27B1, MDM2, MMP2, CASP3. Molecular docking results showed that UDCA had good binding activity with the core targets. In vitro experiment showed UDCA effectively inhibited the proliferation, migration, invasion and neovascularization in colon cancer cells. Conclusion The antitumor activity of ursodeoxycholic acid may be related to cell apoptosis, proliferation, migration and vascular neogenesis. Keywords: ursodeoxycholic acid, colon adenocarcinoma treatment, bile acid metabolism, cell death Background Colon adenocarcinoma is at the forefront of the incidence and mortality of malignant tumors. According to the 2023 Global Cancer Data report, the incidence of colon adenocarcinoma ranks the third among malignant tumors (10%) and the mortality ranks the second among malignant tumors (9.4%),[48]1,[49]2 which seriously endangers public life and health. The pathogenesis o f colon adenocarcinoma is complex, and the only way to treat the disease with drugs is intravenous infusion of fluorouracil drugs, such as the combination of 5-fluorouracil and folinic acid (5 Fu/LV), which has high toxicity, heavy adverse reactions and lack of specificity,[50]3 and with the increase of treatment course, many patients are faced with the challenge of drug resistance. Although radiotherapy can accurately target the tumor, it inevitably causes damage to the surrounding normal tissue during irradiation, resulting in inflammatory reactions such as radiation enteritis. In addition, the high-energy radiation used in radiotherapy will increase the risk of recurrence of tumors, especially for young colorectal adenocarcinoma patients, pelvic radiotherapy may affect the function of the ovary or testis, causing long-term damage to the patient’s reproductive system. Recently, some studies have combined bioinformatics methods to clarify the key roles of metabolism, inflammation, and bile acid secretion in colon adenocarcinoma.[51]4,[52]5 However, none of these studies applied weighted gene co-expression network analysis (WGCNA) to look for Hub genes in colon adenocarcinoma. At present, researchers have developed drugs targeting metabolism and inflammation, such as STF-31, V-9302, and achieved good clinical benefits. However, there are few studies on the application of drugs targeting bile acid metabolism in the treatment of colon adenocarcinoma.[53]6,[54]7 Therefore, in order to deeply understand the molecular mechanism of colon adenocarcinoma, In addition, to verify whether improving bile acid metabolism can alleviate colon adenocarcinoma (CRC) progression, we innovatively combined WGCNA and network pharmacology in this study to explore the feasibility of ursodeoxycholic acid in the treatment of colon adenocarcinoma. Network pharmacology is an emerging research method to analyze the mechanism of drug multi-target treatment of diseases from the system level, which can effectively locate drug targets and establish a “drug-target-disease” biological network to reveal the regulation principle of small molecules in a high-throughput manner.[55]8,[56]9 Network pharmacology explores the mechanism and pathway of drug action, which helps to improve the therapeutic effect of drugs, reduce toxic side effects, and improve the success rate of clinical trials of new drugs. Ursodeoxycholic acid (UDCA) is a natural hydrophilic bile acid that is actively used in the treatment and prevention of cholelithiasis by inhibiting the activity of cholesterol 7A-hydroxylase and cholesterol esterase in bile, reducing cholesterol synthesis and esterification, and thus reducing cholesterol saturation.[57]10 In addition, UDCA can improve the function and expression of Cl-/HCO3- cotransporter AE2 (Anion2), which is a key transporter for bile outflow,[58]11 and the expression of bile acid transporters in the capillary bile ducts of hepatocytes is upregulated,[59]12 such as BSEP, MRP2, thereby promoting bile, bile acid outflow, widely used clinically in primary biliary cholangitis (PBC) and primary sclerosing cholangitis (PSC). In recent years, the potential of UDCA in cancer adjuvant therapy has attracted much attention. Narisawa et al[60]13 found that ursodeoxycholic acid (UDCA) could reduce the incidence of colon adenocarcinoma in animal models of colon adenocarcinoma induced by N-methylnitrosourea, possibly because the concentration and composition of bile acids in feces affected the abundance and distribution of intestinal flora, thus affecting the occurrence and development of colon adenocarcinoma. This study investigated the mechanism of ursodeoxycholic acid in reducing the occurrence of colon adenocarcinoma through network pharmacological research, and provided a new treatment plan for colon adenocarcinoma patients, so as to reduce the incidence and mortality of colon adenocarcinoma. Methods and Materials Data Source and Pre-Processing In this study, the cancer type was selected as Colon adenocarcinoma (TCGA-COAD) and the file data type was selected as RNA-seq in the Cancer Genome Atlas (TCGA) database ([61]https://portal.gdc.cancer.gov/). Directly download the expression data and patient-related clinical data, and use Perl script and python script to integrate the data extraction. RNA sequencing data corresponded to 41 normal tissue samples and 455 tumor samples from 456 colon cancer patients. If the first two principal components determined by principal component analysis could not distinguish tumor tissue from normal tissue, we exclude the sample. Patient samples with incomplete clinical data should also be excluded. The research workflow is shown in [62]Figure 1. Figure 1. [63]Figure 1 [64]Open in a new tab Study flow chart. Ethics Approval This study involving human data, including publicly available data, was reviewed and approved by The Second Affiliated Hospital of Nanjing Medical University Ethics Committee (KY-2020-092). DEGs Screening Software (edgeR package.) was used to analyze the COUNT data of mRNA in tumor tissue and normal tissue to obtain cancer-related differentially expressed genes (DEGs) in COAD. The thresholds for differential expression screening were |logFC|>2 and FDR<0.05. The DEGs volcano map was drawn using “ggplo2” package in the R. Weighted Gene Co-Expression Network Analysis R-Software (WGCNA package) was used to perform gene module clustering for all the genes obtained in 1.1, and the correlation between module eigengene (ME) and clinical characterization was evaluated. The soft threshold β was set to 7 to ensure a scale-free network, and the blockwiseModules function was used to construct a clustering dendrogram, and the branches of the clustering tree were merged into different gene modules by the dynamic tree clipping method, with the minimum number of genes set to 30 and the heterogeneity threshold set to 0.25. We evaluated the association of different modules with the onset of colon adenocarcinoma and selected the most relevant modules as WGCNA-derived hub genes, and visualized the association of clinical characterization data with gene modules using R-package “pheatmap”. Go/KEGG Analysis Gene ontology (GO) and KEGG enrichment analysis were performed using DAVID ([65]https://david.ncifcrf.gov/)[66]14 in order to facilitate our comprehension of the potential mechanism of COAD. Prediction of Ursodeoxycholic Acid and Its Targets Ursodeoxycholic acid related targets were obtained through the search of Drugbank database. In order to increase the reliability of the data, it is necessary to consider the combination of multiple databases and text mining tools to supplement the collection of other potential targets. Firstly, the SMILE and MOL format files of drugs were searched through Drugbank, and then imported into Swiss Target Prediction (STITCH) ([67]http://www.swisstargetprediction.ch/) and Similarity ensemble approach (SEA) to obtain the potential targets of the active ingredients. At the same time, using PharmMapper ([68]http://www.lilab-ecust.Cn/PharmMapper/submitfile.html) to further supplement the potential targets of ursodeoxycholic acid using the reverse pharmacophore matching method. Subsequently, the predicted targets mentioned above were merged and imported into the UniProt database ([69]http://www.Uniprot.org/), limiting the species to Homo Sapiens, and corrected for the official name of the gene. Collection of Targets Related to Colon Adenocarcinoma First, the intersection of DEGs and WGCNA-derived hub genes was obtained. Then enter the keyword “Colon adenocarcinoma” and select “Homo sapiens” as the species at GeneCards ([70]http://www.genecards.Org/),[71]15 Drugbank ([72]https://go.drugbank.com/),[73]16 DisGenet ([74]https//:www.disgenet.org/),[75]17 OMIM[76]18 ([77]https://omim./) and TDD ([78]https://db.idrblab.net/ttd/) data set, searching the disease-related genes. Finally, the above two parts are merged. Drawing Wayne Diagram By introducing compound targets and disease targets into R-software, overlapping genes are potential targets for bioactive compounds that interact with the body in patients. Constructing Protein-Protein Interaction Network (PPI) The STRING ([79]https://stringdb.org) database collects a large number of known or predicted results of protein-protein interactions (PPI).[80]19 The overlapping genes were uploaded to the database, the species selected “Homo sapiens”, set the confidence level (0.40) as the minimum interaction score requirement, and removed the discrete targets. Plot a network of interactions between targets. The evaluation index of network node importance is expressed in terms of betweenness centrality, closeness centrality and degree. Cytoscape 3.8.0 software was used to analyze the PPI network. Molecular Docking Results Open the ursodeoxycholic acid 3D structure file downloaded by PubChem ([81]https://pubchem.ncbi.nlm.nih.gov/). Download the human target protein structure corresponding to ursodeoxycholic acid and key target gene of COAD from RCSB protein database ([82]https://www.rcsb.org/). PyMOL1.0, AutoDockTools1.5.6 and other software were used for docking and visualization of key target molecules. Before docking, root mean square deviation (RMSD) of ligand molecules was calculated to evaluate the rationality of docking parameter setting, and RMSD<2 was selected as the optimal conformation ligand molecule. Binding energy less than −5.0kcal/mol indicates good binding activity. The smaller the binding energy, the more stable the molecular docking. Cell Culture Human colorectal cancer cell lines HCT116 and SW620 were purchased from the Chinese Academy of Sciences (Shanghai) Cell Bank. Human umbilical vascular endothelial cells were purchased from Beijing Huizhiheyuan Biotechnology Co. HCT116, SW620 and HUVEC cells were cultured in DMEM medium (Gibco, USA) with 10% fetal bovine serum (FBS) (Gibco, USA) and 1% penicillin/streptomycin (Thermo, USA). All cell lines were incubated in a 37°C incubator with 5% CO2 according to the manufacturer’s instructions. Cell line morphology was evaluated regularly to ensure cell viability and quality. Drug Ursodeoxycholic acid (UDCA), specification: 1 g/ bottle (Quality standard: >98%, BR.), manufactured by Nanjing Best Biotechnology Co., LTD., Lot No.: A04088.0.1g of UDCA powder was dissolved in 2mL ethanol, which was stock1. 400ul stock1 and 1ul sterile dimethyl sulfoxide (DMSO, Sigma-Aldrich, USA) were added to 600ul PBS to form stock2. stock2 was diluted with DMEM to different concentrations (0, 0.00625, 0.0125, 0.025, 0.05, 0.1mg/mL) so that the final concentration of DMSO in the medium did not exceed 0.1% in all treatments. Cell Viability Assay Cell proliferation and cytotoxicity were detected by MTT colorimetry. SW620 and HCT116 cells of logarithmic growth stage were inoculated in 96-well plates with a density of 1×104 cells/well and cultured in an incubator at 5%CO2 and 37°C for 24h. After the cells grew to 80% under a microscope, the supernatant was discarded. Media containing different concentrations of UDCA (0,0.00625,0.0125,0.025,0.05,0.1mg/mL) were added to each well (each group has 5 multiple holes). After 24h culture, 10μL solution containing 10% MTT was added to each well and cultured for 4h away from light. The supernatant was discarded and I00μL DMSO was added to each well for 15 minutes. The absorbance was read at 490 nm on the enzyme-labeler of Bio-Rad Laboratories (United States). Flatbed Cloning Experiment SW620 cells of logarithmic growth stage were inoculated in 6-well plates with a density of 1000 cells/well, cultured in an incubator at 5%CO2 and 37°C for 24h, then the supernatant was discarded, cleaned with lxPBS, and DMEM and culture medium containing 0.02mg/mL UDCA were added. Two weeks later, the supernatant was discarded, cleaned with 1xP BS, and 1mL of pre-cooled paraformaldehyde was added to fix it for 15min. After cleaning with ultra-pure water, lmL of crystal violet dyeing solution containing 0.1% was added to stain for 15min. Then, it was gently rinsed and dried for photo taking and analysis. Cell Migration and Invasion Assays Both migration and invasion tests were performed using 24-well plates. Serum-free culture-supplemented cells were inoculated in the upper chamber of Transwell. The subchamber of Transwell was filled with 600 μL medium containing different concentrations (0.02mg/mL) of UDCA and 10% fetal bovine serum. After incubation at 37°C for 24h, the cells were fixed with cold methanol for 20 minutes and then stained with 0.1% crystal violet for 20 minutes. The cells were cleaned with lxPBS and fixed with pre-cooled paraformaldehyde for 15min. After washing, 1mL of 0.1% crystal violet staining solution was added to the chamber for 15min. The cells were washed again and dried to capture cell images in the field of view of each filter under the microscope. The only difference between migration and invasion tests was whether the upper surface of the transwell membrane was pre-coated with Matrigel matrix (BD, USA). Scratch Experiment SW620 and HCT116 cells of logarithmic growth stage were inoculated into 6-well plates and incubated until the cell density was about 90–95%. A vertical line was drawn vertically in the center of each well with a 10μL gun tip, and washed twice with IxPBS. Culture medium with different concentrations (0, 0.01, 0.02, 0.04mg/mL) of UDCA was added to each well. They were observed and photographed under microscope at 0 and 24h, respectively. Angiogenesis Experiment HUVEC cells were treated with different concentrations of UDCA (0, 0.01, 0.02, 0.04 mg/mL) for 24 h. After adding 30 μL of Matrigel Matrix Gel to a 24-well plate and incubating in a 37°C incubator for half an hour, each group of HUVEC cells was made into a cell suspension at a cell density of 3×10^5 cells/mL, and then inoculated in 24-well plates and continued to be incubated in the incubator. Photos were taken at 2h, 6h and 12 h, respectively. Results and Analysis DEGs Screening We first obtained the data of colon adenocarcinoma from the TCGA database. Based on principal component analysis, it was found that tumors and normal samples could be well distinguished. Then, differential expression threshold =|logFC|>2 and FDR<0.05 were used as screening criteria to screen differentially expressed genes between tumor samples and near-normal tissue samples. A total of 2064 differential genes (1000 up-regulated and 1064 down-regulated) were identified ([83]Figure 2A and [84]B). Figure 2. [85]Figure 2 [86]Open in a new tab Identification of DEGs between 41 normal and 455 COAD samples. (A) Principal component analysis. (B) Volcano plot. Blue dots represent genes down-regulated in COAD; gray dots, genes not differing significantly between COAD and normal tissues; and red dots, genes up-regulated in COAD. Construction of WGCNA Network and Identification of Colon Adenocarcinoma Related Modules To identify the potential gene modules with the highest association with colon adenocarcinoma, we performed WGCNA analysis of all candidate genes from the TCGA colon adenocarcinoma data set ([87]Figure 3A) and obtained correlation coefficients. Finally, the Brown module was the most associated with tumor among all modules (p=0.91) ([88]Figure 3B). Figure 3. [89]Figure 3 [90]Open in a new tab Identification of CRC-associated gene modules in TCGA dataset using WGCNA. (A) Dendrogram of all genes in the TCGA dataset was clustered on the basis of a topological overlap matrix (1-TOM). Each branch in the clustering tree represents a gene, while co-expression modules were constructed in different colors. (B) Module-trait heatmap of the correlation between the clustering gene module and COAD in the TCGA dataset. Each module contains the corresponding correlation coefficient and p value. Go/KEGG Analysis to Find Critical Paths We ultimately screened 128 overlapping genes as candidate hub genes that may play an important role in the development and progression of colon adenocarcinoma ([91]Figure 4A). GO and KEGG analyses were then performed to further explore the potential role of these 128 overlapping genes ([92]Figure 4B and [93]C). GO enrichment analysis showed that overlapping genes mainly affect biological functions such as metabolism, exosomes, and various hydrolase activities. KEGG enrichment analysis showed that overlapping genes mainly affected metabolism and bile acid secretion pathways. Figure 4. [94]Figure 4 [95]Open in a new tab Candidate hub genes were screened and validated. (A) Venn diagram revealed 128 overlapping candidate hub genes (B) GO enrichment analysis of candidate hub genes. (C) KEGG pathway analysis of candidate hub genes. Prediction of Ursodeoxycholic Acid and Its Targets Through database analysis, 164 ursodeoxycholic acid-related target genes and 1202 colon adenocarcinoma targets were identified. At the intersection of the two types of targets, 26 potential intersection targets of ursodeoxycholic acid anti-tumor were obtained ([96]Figure 5A). Figure 5. [97]Figure 5 [98]Open in a new tab Ursodeoxycholic acid anti-tumor network diagram. (A) Venn diagram revealed 26 intersection targets of ursodeoxycholic acid drug targets and COADdisease targets. (B) Ursodeoxycholic acid anti-tumor protein-protein interaction (PPI) network diagram. (C) Diagram of the core anti-tumor targets of ursodeoxycholic acid. Construction of Ursodeoxycholic Acid-Colon Adenocarcinoma Related PPI Network and Screening of Core Targets The String database was applied to mine ursodeoxycholic acid-colon adenocarcinoma related PPI network ([99]Figure 5B). Subsequently, secondary screening was performed using the CytoNCA plug-in in Cytoscape, and the core network diagram was finally obtained ([100]Figure 5C), and TNF, CYP27B1, CASP3, MDM2, MMP2 may have a vital function in the treatment of this disease. Molecular Docking Molecular docking of five representative core targets: TNF, CYP27B1, MDM2, MMP2 and CASP3 obtained in the previous step were performed with ursodeoxycholic acid respectively. Before molecular docking, the root mean squared deviation (RMSD) was evaluated, molecules with RMSD<2 were selected for docking and the key target docking information table was drawn. After molecular docking, the binding mode diagram of ursodeoxycholic acid and core target protein receptors was drawn ([101]Figure 6). The results showed that the binding energy of ursodeoxycholic acid and 5 key target proteins were all <-5.0kcal/mol, indicating that it had good binding activity with all 5 key target proteins ([102]Table 1). Figure 6. [103]Figure 6 [104]Open in a new tab Molecular docking diagram of ursodeoxycholic acid and key targets of CRC. (A) Ursodeoxycholic acid and TNF; (B) Ursodeoxycholic acid and CYP27B1; (C) Ursodeoxycholic acid and MDM2; (D) Ursodeoxycholic acid and MMP2; (E) Ursodeoxycholic acid and CASP3. Table 1. Results of Ursodeoxycholic Acid Docking Scores with Key Targets in CRC Key Targets in CRC PDB Protein Binding-energy/kcal·mol–1 TNF 5INL −5.33 CYP27B1 CYP27B1 −5.73 MDM2 2P54 −5.31 MMP2 8H78 −5.04 CASP3 3kjf −5.21 [105]Open in a new tab UDCA Can Inhibit the Proliferation of Colorectal Cancer Cells in vitro In order to evaluate the effect of UDCA on the proliferation of SW620 and HCT116 cells, SW620 and HCT116 cells were treated with different concentrations of UDCA (0, 0.00625, 0.0125, 0.025, 0.05, 0.1mg/mL), respectively, and the cell growth was detected by MTT assay. In [106]Figure 7A, UDCA significantly inhibited cell proliferation, indicating that UDCA had potential anti-proliferation effects. IC50 calculations after 24 h of treatment showed that the IC50 values of SW620 and HCT116 cells were 0.02mg/mL and 0.016mg/mL, respectively. The clonogenesis experiment results showed that compared with the control group, the clonogenesis ability of UDCA treated cells was significantly decreased ([107]Figure 7B). Figure 7. [108]Figure 7 [109]Open in a new tab UDCA inhibits the proliferation of colorectal cancer cells. (A) SW620 and HCT116 cells were treated with UDCA at different concentrations (0, 0.00625, 0.0125, 0.025, 0.05, 0.1mg/mL) for 24 h, and cell viability was determined by MTT method and inhibition rate was calculated. (B) The effect of UDCA on the clonogenesis of SW620 and HCT116 cells. Data were expressed as mean ±SEM (n = 6). UDCA Inhibits the Migration and Invasion of Colorectal Cancer in vitro We studied the effect of UDCA on colorectal cancer cell migration in vitro by using scratch assay. As shown in [110]Figure 8A and [111]B, in SW620 and HCT116 cells, the degree of wound healing was gradually reduced after 0, 12, 24 and 36 hours of UDCA treatment compared with the control group. The cell migration and matrix gel invasion tests were also performed. The results in [112]Figure 8C and [113]D showed that the number of UDCA-treated colorectal cancer cells crossing the membrane was significantly reduced compared with the control group. The results showed that UDCA could effectively inhibit cell migration and invasion. Figure 8. [114]Figure 8 [115]Open in a new tab UDCA suppressed the migration and invasion in SW620 and HCT116 cells. (A) SW620 cells and (B) HCT116 cells were scraped with a pipette tip and then cultured with UDCA at the indicated concentrations for different times. The scratch wound cell migration assays were detected by microscope. (C) SW620 cells and (D) HCT116 cells treated with or not with UDCA were performed for migration and invasion assays. The cells that invaded across the membrane of the transwell were stained with crystal violet. UDCA Inhibits the Formation of Blood Vessels in Colorectal Cancer in vitro Angiogenesis experiment was used to study the effect of UDCA on the new blood vessels of colorectal cancer cells. As shown in [116]Figure 9, in HUVEC cells, compared with the control group, new vessels were significantly reduced in the UDCA-treated group after 2h, 6h, 12h, among which the inhibition effect was most obvious in the 0.02mg/mL and 0.04mg/mL treatment groups. Figure 9. [117]Figure 9 [118]Open in a new tab UDCA inhibits angiogenesis in HUVECs. HUVEC cells were cultured with UDCA at specified concentration, and the changes of angiogenesis capacity in each group were observed under microscope at 2h,6h and 12h. Discussion Colon adenocarcinoma is a malignant tumor occurring in the mucosal epithelium of the colon. At present, surgery combined with chemotherapy is a comprehensive treatment method that is widely used, but it is not suitable for advanced and older patients. Targeted therapeutic drugs such as cetuximab and bevacizumab are expensive, so this paper aims to find an adjuvant drug that can effectively prolong the survival of cancer patients. In this work, we used a public database to obtain the differential genes between colon adenocarcinoma and normal tissue samples, and then used WGCNA to identify the modules most associated with colon adenocarcinoma, where the two intersected to get the core genes, which greatly narrowed the range of genes to screen for. Further enrichment analysis showed that differential genes were mainly enriched in metabolic pathways and bile acid metabolic pathways. Firstly, the genes enriched in the metabolic pathway include AHCYL2, ADH1C, LDHD, PDE9A, AKR1B10, PLAAT2, etc. In previous studies, cyclin-dependent kinase 7 (CDK7) phosphorylated nuclear Yes-associated protein 1 (YAP) at S127 and S397, enhancing its transcriptional function and promoting the expression of LDHD protein. CDK7-YAP-LDHD enables cancer cells to escape from iron death induced by D-lactic acid, and targets LDHD to inhibit cancer cell proliferation and renewal.[119]19 PDE9A encodes phosphodiesterase 9A, a family of enzymes that catalyze the hydrolysis of cyclic nucleotides, hydrolyzing cyclic nucleotides, thereby reducing intracellular cAMP (cyclic adenosine monophosphate) and cGMP (cyclic guanosine monophosphate) levels. PDE9A specifically targets cGMP. Phosphodiesterase inhibitors activate caspase, inhibit Wnt/b-catenin pathway and TCF transcription through the increase of cGMP, thereby inhibiting CDK and survivin’s anticancer effects.[120]20,[121]21 In addition, inhibiting the expression of AKR1B10 and PLAAT2 can affect the production of lipid signaling molecules, which contribute to the proliferation, survival and invasion of cancer cells.[122]22,[123]23,[124]24 Inhibitors targeting LDHD, PDE9A, and AKR have been shown to negatively regulate cancer cells and are being explored as potential therapeutics. In addition, many drugs targeting tumor metabolism have been applied in clinical practice and achieved good survival benefits.[125]25 The authors’ KEGG results showed that, in addition to the enrichment of metabolic pathways, differential genes were the next most enriched in bile acid metabolic pathways. However, there are few reports on the relationship between drugs targeting bile acid metabolism and colon adenocarcinoma. We found that many UGT family molecules were enriched in the core genes of bile acid metabolism pathway: UGT2B15, UGT2A3, UGT2B17, etc. By searching the literature, we found that ursodeoxycholic acid can regulate bile acid metabolism through UGT3A1 and UGT2B7,[126]26–28 and bile acid metabolism has a regulatory effect on the composition and abundance of intestinal flora.[127]29,[128]30 The role of intestinal flora and its metabolites in colon cancer and other cancers has been widely confirmed.[129]31,[130]32 So we hypothesized that ursodeoxycholic acid could alleviate colon adenocarcinoma progression. Next, we explored the feasibility of ursodeoxycholic acid as an adjuvant therapy for colon adenocarcinoma based on network pharmacology. According to the database, 164 ursodeoxychacid-related target genes and 1202 colon adenocarcinoma targets were identified, including 26 common targets. Five representative core targets including TNF, CASP3, CYP27B1, MDM2 and MMP2 were further selected by String and Cytoscape software. The TNF gene, which encodes tumor necrosis factor α (TNF-α), is a cytokine involved in inflammation and immune system regulation. The inflammatory microenvironment is closely related to tumor invasion and metastasis, and chronic inflammation may be one of the important reasons leading to tumor invasion and metastasis. TNF-αwas found to activate MAPKKKs through TRAF2 aggregation and JNK and p38 through MKK4 and MKK7 (MAP ERK kinases). Min-Kyung Choo et al found that TNF-α can enhance the lung metastasis ability of colon adenocarcinoma cells by activating the ERK signaling pathway.[131]33 Simone et al found that TNF can promote the growth of colon cancer cells by activating STAT3 and NF-kB.[132]34 Caspase3 encoded by CASP3 plays a role as a core molecule in cell apoptosis. Previous studies have believed that caspase-3 can eliminate cancerous cells and is a widespread tumor suppressor.[133]35 Surprisingly, Liu et al found that Caspase-3 can promote the genetic instability of cells during radiation and chemotherapy. It plays an important role in tumor transformation.[134]36,[135]37 The CYP27B1 encoded hydroxylase CYP27B1 is primarily responsible for converting 25-hydroxyvitamin D (calcifediol) to its active form 1, 25-dihydroxyvitamin D (calcitriol), After calcitriol binds to vitamin D receptor (VDR) or retinoid-related orphan receptors α or γ (RORα and RORγ), it binds to Vitamin D response elements (VDRE) located in the promoter region to regulate vitamin D related genes. A large number of cell and animal experiments have confirmed that vitamin D can inhibit colon cancer through a variety of ways. Kaler et al found that vitamin D can inhibit the growth of colon adenocarcinoma cells by inhibiting the production of IL-I-β by macrophages and thus blocking the Wnt signaling pathway.[136]38 It was subsequently found that IL-IB could also protect colon adenocarcinoma cells from apoptosis mediated by tumor necrosis factor-associated apoptosis-inducing ligand (TRAIL).[137]39 A recent study also found that vitamin D may promote iron death in colon adenocarcinoma stem cells by regulating SLC7A11 expression.[138]40 MMP2 is a widely studied gene closely related to cell migration.[139]41 MDM2 is an E3 ligase, and studies have found that MDM2 interacts with p53 and inhibits the transcriptional activity of p53[140]42 P53 is not only a tumor suppressor gene, but also a key gene in the regulation of cellular radiosensitivity and chemotherapy sensitivity. Overexpression of MDM2 can inhibit p53-mediated ability to detect DNA damage in G1 phase and block the cell cycle.[141]43 Inhibition of MDM2 expression in glioma cells can also play an important role in resistance to temozolomide through regulation of p53/MDM2 signaling pathway.[142]44 Therefore, it is speculated that ursodeoxycholic acid acts on the above target proteins, which may play a role in adjuvant colon adenocarcinoma therapy through the regulation of cell apoptosis, iron death, cell proliferation and migration and other multicellular biological processes. The results of molecular docking showed that the binding energy of the above target proteins to ursodeoxycholic acid was less than −5 kcal/mol-1, indicating stable ligand-receptor binding. This result further verifies the reliability of network pharmacological predictive analysis. Goossens et al found that UDCA can protect epithelial cells from injury and apoptosis, while inducing inhibition of cancer cell proliferation and apoptosis and/or autophagy death.[143]45 Similar results were obtained in our study, and the clonogenesis experiment showed that the clonogenesis ability of UDCA-treated cells was significantly reduced compared with the control group. In addition to the optimal inhibitory concentration (0.02mg/mL), we also set low concentration (0.01mg/mL) and high concentration group (0.04mg/mL) for scratch test, which results showed that the healing ability of UDCA treated cells was significantly weakened compared with the control group. After 36 hours, the inhibitory effect of 0.01mg/mL and 0.02mg/mL groups was still good, but the inhibitory effect of 0.04mg/mlUDCA group was worse, which may be because the cells were treated with higher concentration of drugs for a long time, and the sensitivity of cells to drugs changed to a certain extent. The Transwell chamber experiment showed that the number of UDCA-treated colorectal cancer cells passing through the chamber was significantly reduced compared to the control group, and consistent results were obtained with the addition of matrix glue. The angiogenesis experiment showed that compared with the control group, the angiogenesis ability of HUVEC cells after UDCA treatment was significantly inhibited, especially in the medium concentration and high concentration treatment groups. These results confirm that UDCA has an inhibitory effect on the proliferation, migration and invasion of colorectal cancer. In summary, this study was conducted based on the bile acid pathway enriched by differential genes between colon adenocarcinoma and normal people. It was found that this pathway enriched a large number of UGT family molecules, and a drug targeting this molecule, namely ursodeoxycholic acid, was found. Subsequently, network pharmacology and molecular docking technology were used to verify its feasibility in the treatment of colon adenocarcinoma at the molecular level, in vitro experiments were performed to validate its inhibitory effect on proliferation, migration and invasion of colorectal cancer cells at the cellular level. However, this study has some limitations. First, data on drug and disease targets are obtained from literature and databases, so the reliability and accuracy of the method’s predictions depend on the quality of the data. Secondly, this study is based on data mining and bioinformatics analysis, and further clinical trials and animal experiments need to be verified. Conclusion Based on the sequencing data, this study innovatively searched for COAD enrichment pathways and key targets, identified ursodeoxycholic acid as the target drug, and then used network pharmacology methods to explore the molecular mechanism and feasibility of its treatment of COAD. Through the screening and re-study of known drugs, new uses of the drug were realized, and it was found that it may have effects on the key targets of apoptosis, proliferation and migration, and has good anti-tumor activity, providing new ideas and new hopes for anti-tumor drug research and development. To clarify the anti-tumor potential of UDCA, more animal and human studies are needed in the future. Funding Statement There is no funding to report. Abbreviations CRC, Colorectal cancer; COAD, Colon adenocarcinoma; DEGs, differentially expressed genes; DMSO, dimethyl sulfoxide; GO, gene ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; PPI, Protein-Protein Interaction Network; RMSD, root mean squared deviation; SEA, Similarity ensemble approach; TCGA, the cancer genome atlas; TOM, topological overlap matrix; UDCA, ursodeoxycholic acid; WGCNA, weighted gene co-expression network analysis. Data Sharing Statement The original contributions presented in the study are included in the article. The dataset used in this study is available in an online database and the name and registration number of the dataset can be found in the article. Further inquiries can be directed to the corresponding authors. Author Contributions All authors made a significant contribution to the work reported, whether that is in the conception, study design, execution, acquisition of data, analysis and interpretation, or in all these areas; took part in drafting, revising or critically reviewing the article; gave final approval of the version to be published; have agreed on the journal to which the article has been submitted; and agree to be accountable for all aspects of the work. Disclosure All authors declare that there exists no potential competing interest. No benefits in any form have been received or will be received from a commercial party related directly or indirectly to the subject of this article. References