Anti-tumor target screening of sea cucumber saponin Frondoside A: a bioinformatics and molecular docking analysis

Cancer remains the leading cause of death worldwide. In spite of significant advances in targeted and immunotherapeutic approaches, clinical outcomes for cancer remain poor. The aim of the present study was to investigate the potential mechanisms and therapeutic targets of Frondoside A for the treatment of liver, pancreatic, and bladder cancers. The data presented in our study demonstrated that Frondoside A reduced the viability and migration of HepG2, Panc02, and UM-UC-3 cancer cell in vitro. Moreover, we utilized the GEO database to screen and identify for differentially expressed genes (DEGs) in liver, pancreatic, and bladder cancers, which resulted in the identification of 714, 357, and 101 DEGs, respectively. Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway annotation were performed using the Metascape database for DEGs that were significantly associated with cancer development. The protein-protein interaction (PPI) networks of the identified DEGs in liver, pancreatic, and bladder cancers were analyzed using Cytoscape 3.9.0 software, and subsequently identified potential key genes that were associated with these networks. Subsequently, their prognostic values were assessed by gene expression level analysis and Kaplan-Meier survival analysis (GEPIA). Furthermore, we utilized TIMER 2.0 to investigate the correlation between the expression of the identified key gene and cancer immune infiltration. Finally, molecular docking simulations were performed to assess the affinity of Frondoside A and key genes. Our results showed a significant correlation between these DEGs and cancer progression. Combined, these analyses revealed that Frondoside A involves in the regulation of multiple pathways, such as drug metabolism, cell cycle in liver cancer by inhibiting the expression of CDK1, TOP2A, CDC20, and KIF20A, and regulates protein digestion and absorption, receptor interaction in pancreatic cancer by down-regulation of ASPM, TOP2A, DLGAP5, TPX2, KIF23, MELK, LAMA3, and ANLN. While in bladder cancer, Frondoside A regulates muscle contraction, complement and coagulation cascade by increase FLNC expression. In conclusion, the present study offers valuable insights into the molecular mechanism underlying the anticancer effects of Frondoside A, and suggests that Frondoside A can be used as a functional food supplement or further developed as a natural anti-cancer drug.


Introduction
Cancer remains the predominant cause of mortality globally.In 2020, according to the GLOBOCAN report, the annual incidence of new cancer cases is estimated to be 19.3 million, with nearly 10.0 million cancer-related deaths (1).The projected incidence of new cancer cases in China and the USA for the year 2022 is around 4,820,000 and 2,370,000, respectively, with estimated cancer-related deaths of approximately 3,210,000 and 640,000 in these two countries (2).In the past decades, many basic and clinical studies have investigated underlying mechanisms of cancer formation and progression leading to advances in the molecular biology and diagnosis of cancers.However, despite these advancements, the incidence and mortality rates of cancer continue to exhibit high levels, and the available treatment options remain limited in their efficacy and scope (3).While there are many treatment options currently available for cancer patientsand many are clinically effectivethese treatment options are known to have significant side effects (4,5).Consequently, the development of novel and efficacious therapeutic agents for the treatment of cancer constitutes an urgent requirement.In general, natural products have served as a principal reservoir of compounds for the management of a diverse spectrum of cancers, presenting encouraging prospects for examining not only novel anticancer agents but also unexplored and potentially significant mechanisms of action (6).Marine organisms have emerged as a recent research focus for the identification of novel compounds that hold potential in the prevention and treatment of cancer (7), and there is growing interest in exploiting the diverse and complex marine organisms for rational drug discovery (8).The advancement of technology and extensive research on marine natural products have culminated in the identification of a novel cohort of anti-cancer drugs, and a number of marine biological extracts are currently being studied in clinical trials after showing anti-cancer activity in various preclinical studies (7,9,10).So far, 20 drugs from marine sources are already in clinical use (10).
To date, the exploration for anticancer agents in diverse marine phyla has unveiled numerous active compounds (11).The echinoderms, which exclusively inhabit the marine milieu and encompass sea stars, sea urchins, sand dollars, sea cucumbers, and sea lilies, have garnered significant attention as a subject of interest, amongst which sea cucumbers are of particularly interesting as they have been used as dietary remedies for cancer, inflammation, and other diseases for centuries (12).For example, triterpenoid glycosides from diverse species of sea cucumbers have been recognized for their potential anti-cancer properties (13).Frondoside A, a specific triterpenoid saponin, has recently gained significant attention due to its demonstrated potent anti-cancer effects against various types of solid malignancies and leukemias (14).However, few studies have attempted either Frondoside A as a potential treatment option for liver, pancreatic, and bladder cancers or investigated their mechanisms of action.In this study, we analyzed the in vitro antitumor activity of Frondoside A and its potential targets, and mechanisms of action related to the treatment of liver, pancreatic and bladder cancers using in vitro cellular assays, bioinformatics and molecular docking techniques.The full analytical work flow of this study is shown in Figure 1.buffer solution (PBS, pH 7.4) were purchased from Thermo Fisher Scientific (Waltham, MA.Finally, penicillin (100 Units/mL), streptomycin (100 µg/mL), 0.25% Trypsin (×1), and Nonessential Amino Acid Solution (NAAS, 100X) were obtained fromGibco (Grand Island, NY) were used in the study Chemotherapeutic agent Epirubicin hydrochloride (EPI), a commonly used anti-tumor drug to treat solid tumors, for instance liver cancer, breast cancer, and bladder cancer, was purchased from Aladdin.

Cell lines and cell cultures
Three cancer cell lines were used in this study, including Human liver cancer cell HepG2 (ATCC ® No. HB-8065), Human pancreatic cancer cell Panc02 (ATCC ® No. CRL-2553) and Human bladder cancer cell UM-UC-3 (ATCC ® No. CRL-1749).Both HepG2 and Panc02 cells were incubated with DMEM medium supplemented with 10% FBS, and UM-UC-3 cells were cultured in MEM medium supplemented 1% NAAS.Moreover, three cell lines were supplemented with 10% FBS, 100 Units/mL of penicillin, and 100 mg/mL, and incubated at 37°C with 5% CO 2 atmosphere and 100% humidity.

Cell viability assay
To evaluate the cytotoxicity of drugs of interest, HepG2, Panc02, and UM-UC-3 cells were seeded in 96-well plates at a density of 5000 cells per well and incubated overnight.Triplicate treatments with increasing concentrations of Frondoside A were administered to the cells for 48 hours, while EPI (10 mM) and PBS (pH 7.4) were used as positive and blank controls, respectively.Subsequently, PBS (pH 7.4) was used to carefully wash the cells, after which 10% CCK-8 reagent was added, and the cultures were incubated for 3 hours.The plates were analyzed at 450 nm utilizing a microplate reader (Bio-Rad, Hercules, CA).Cell viabilities were expressed as proportional viabilities (%) normalized against the blank control and assumed 100% survival.The half maximal inhibitory concentration (IC 50 ) values of Frondoside A for HepG2, Panc02, and UM-UC-3 cell lines were also determined.

Cell migration assay
Confluent cell monolayers of HepG2, Panc02, and UM-UC-3 cells were cultivated in 6-well tissue culture plates.A plastic micropipette tip was used to carefully scratch the confluent cell monolayers, and the cell cultures were subsequently washed twice with PBS (pH 7.4).The cells were cultured in their respective cell culture media that contained 1% FBS and supplemented with Frondoside A at their respective IC 50 concentrations for 24 h.The widths of the wound in the cultures were documented using an Olympus inverted microscope at 0, 12, and 24 h after the wound.Finally, the scratched area was measured using a NIH Image J (Bethesda, MD), and the degree of wound closure was calculated.

Microarray data collection
The Gene Expression Omnibus (GEO) database (https:// www.ncbi.nlm.nih.gov/geo/) is a publicly available functional genomics resource that provides high-throughput gene expression data, chips, and microarrays (15).For the relevant tumors, gene expression profile data was acquired from the GEO database.The inclusion criteria were: (1) The specimens were obtained from patients with human liver cancer, pancreatic cancer, and bladder cancer; (2) The study population was comprised of both tumor patients and normal controls; (3) The study design was "expression profiling by array".Finally, we downloaded liver cancer-related datasets (GSE14520 and GSE60502), pancreatic cancer datasets (GSE16515 and GSE28735), and bladder cancer datasets (GSE13507, GSE23732, and GSE37815) from the GEO database.

DEGs identification and data visualization
The GEO2R (https://www.ncbi.nlm.nih.gov/geo/geo2r/) is an online analytical tool that is utilized for the identification of differently expressed genes (DEGs) between tumor and healthy controls.The data was categorized into tumor and normal groups and analyzed based on |log 2 FC| ≥ 1.0 and adj.p value < 0.05, respectively.Following the screening of the DEGs within each dataset, the Venny2.0online tool (https://bioinfogp.cnb.csic.es/tools/venny/) was employed to identify common targets within the relevant datasets, and the overlapping DEGs were retained for further analysis.

Gene ontology terms and KEGG enrichment analysis of DEGs
In order to annotate the DEGs identified from the aforementioned comparison groups in a functional manner, GO term annotation and KEGG pathway analysis were conducted with the use of Metascape.Metascape (https://metascape.org/gp/index.html#/main/step1)(16) is a comprehensive analytical website that incorporates functional enrichment, gene annotation, interactome analysis, and membership search, leveraging more than 40 independent knowledge bases.KEGG is a database resource used for elucidating high-level features and effects of biological systems, while GO is a widely used bioinformatics program for high-quality functional gene annotation based on biological processes (BP), molecular functions (MF), and cellular components (CC).Metascape was utilized with a minimum overlap of 3 and a minimum richness of 1.5 as screening criteria to determine the characteristics of DEGs.A significance level of p < 0.01 was considered statistically significant, and GraphPad Prism 9.0 software was used to generate bar graphs to visually represent the obtained results.

Protein-protein interaction network construction and key gene identification
The intersection targets that were obtained were subjected to analysis and construction of Protein-Protein Interaction (PPI) networks using the STRING database (https://string-db.org/).The STRING database is one of the most informative databases for studying protein interactions at this stage, which evaluates and integrates information on known and predicted protein interactions to form a comprehensive protein network covering > 1,100 organisms (17).The STRING database utilizes probability calculations of various evidence channels, including the correction of random interaction probability, to construct PPI networks.The resulting composite score is assigned based on the level of confidence, with a score of 0.4 indicating medium confidence, 0.7 indicating high confidence, and 0.9 indicating the highest level of confidence.In this study, a PPI network with a score of 0.9, i.e., the highest confidence network, was used to import potential targets into the STRING database and select "Homo Sapiens" species to obtain target protein interaction relationships.

Expression analysis of key genes
GEPIA2.0 (http://gepia2.cancer-pku.cn/#index) is a newly developed web server for cancer and normal gene expression analysis and the interactive analysis with rich analysis functions, such as tumor/normal differential expression analysis, survival analysis, gene correlation analysis, and downscaling analysis (18).This study will use the GEPIA2.0database to analyze the obtained key genes and initially validate their expression differences between tumor and normal tissues.

Prognostic analysis of survival of key genes
In this study, GEPIA2.0 was utilized to examine the correlation between the expression of key genes in tumors and overall survival (OS) (18).The statistical analysis was performed by calculating risk ratios (HR) and p values with 95% interval confidence (IC) using a Log-rank test.Finally, genes with insignificant differences were excluded for further analysis to screen for key genes with a significant expression.

Analysis of key genes and tumor immune invasion level
The Tumor Immunity Estimation Resource (TIMER) database contains 10,897 samples from The Cancer Genome Atlas (TCGA) for 32 cancer types (19).In this study, the correlation between key genes and the degree of immune cell infiltration, including B cells, CD8 + T cells, CD4 + T cells, neutrophils, macrophages, and dendritic cells, was investigated using the TIMER2 database (http://timer.cistrome.org/).The tumor purity of PCs was adjusted prior to analysis.Key genes were input into the "gene module" to generate scatter plots, which were used to observe the relationship between gene expression and the level of tumor immune invasion.

Molecular docking
To confirm the binding of Frondoside A to the predicted key targets, the structural formula of Frondoside A was retrieved and d o w n l o a d e d f r o m t h e P u b C h e m d a t a b a s e ( h t t p s : / / www.ncbi.nlm.nih.gov/pccompound) in mol format.The structure was optimized using Chem3D 14.0 software with MM2 force field for energy minimization and saved in pdb format.Threedimensional structures of CDK1 (PDB ID: 4YC6) (20), ASPM (PDB ID: 3QBT) (21), BUB1 (PDB ID: 4QPM) ( 22), KIF20A (PDB ID: 5LEF) (23), TOP2A (PDB ID:5NNE) (24), TPX2 (PDB ID:4C3P) (25), KIF23 (PDB ID: 3VHX) (26), MELK (PDB ID:5TWL) (27), and FLNC (PDB ID: 7OUU) were obtained from RCSB PDB database (https://www.rcsb.org/)(28).The proteins were preprocessed with PyMOL software for dehydration, hydrogenation, and removal of irrelevant ligands.The Grid Box parameters were set with the receptor as the center, and molecular docking analysis was performed with AutoDock Vina 1.1.2.The conformation with the lowest free energy was selected and visualized with PyMOL (29).A binding free energy of < -5 kcal/ mol indicates that the target has some binding activity to the compound (30), and the lower the binding energy, the more stable the receptor-ligand binding (31).

Statistical analysis
The experiments were conducted in triplicate, and the data were presented as mean ± standard deviation (SD).Statistical analysis was performed using unpaired Student's t-test and one-way analysis of variance (ANOVA).Results with a p-value of less than 0.05 were considered statistically significant.

In vitro inhibition of tumor cells by Frondoside A
As shown in Figures 2A to C, Frondoside A exhibited a dosedependent reduction in the cell viability in HepG2, Panc02, and UM-UC-3 cell lines.The IC 50 values of Frondoside A for HepG2, Panc02, and UM-UC-3 were 1.5 mM, 1.5 mM, and 1 mM, respectively.
As shown in Figures 2D to E, both Frondoside A and EPI treatments significantly inhibited cell migration and reduced the proliferation of HepG2, Panc02, and UM-UC-3.In addition, Frondoside A exhibited a more conspicuous suppression of cell migration than EPI (Figures 2G-I) The migration rate of HepG2, Panc02, and UM-UC-3, respectively.

Identification of DEGs
GEO is a free database of microarray/gene profiles and nextgeneration sequencing, from which liver cancer and normal or adjacent liver tissue gene expression profiles of GSE14520 and GSE60502 were obtained.Using adj.p < 0.05 and |log 2 FC|> 1 as cut-off criteria, from the microarray data of GSE14520 and GSE60502, 1238 (573 up-regulated and 665 down-regulated) and 1556 (794 up-regulated and 762 down-regulated) DEGs were extracted, respectively (Figures 3A, B).After integrated bioinformatical analysis, a total of 714 genes were identified from the two profile datasets (Figure 3C).
Similarly, mRNA expression profiles in human bladder cancer tissues and adjacent non-cancerous tissues were also analyzed using the GEO databases (GSE13507, GSE23732, and GSE37815).For GSE13507, there were 459 DEGs (i.e., 63 up-regulated and 396 down-regulated) with significant changes in bladder cancer tissues (Figure 3G).For GSE23732, there were 611 DEGs (i.e., 102 upregulated and 509 down-regulated) significant changes in bladder cancer tissues (Figure 3H).For GSE37815, there were 668 DEGs (i.e., 190 up-regulated and 478 down-regulated) significant changes in bladder urothelial carcinoma tissues (Figure 3I).As a result of the analysis above, a total of 101 genes were used for further investigation (Figure 3J).

Enrichment analysis
To investigate the biological functions of DEGs, functional annotation and pathway enrichment analysis were conducted using the Metascape database.The results were visualized by generating bar graphs with GraphPad Prism 9.0 software.GO (Gene Ontology) analyses are shown in Figures 4A-C, and KEGG analyses are shown in Figures 4D-F.The results of GO analysis of DEG in liver cancer are shown in Figure 4A.The most important biological processes (BPs) of DEGs were "monocarboxylic acid metabolic process", "response to xenobiotic stimulus", "small molecule catabolic process", "response to the hormone", and "mitotic cell cycle process".The molecular functions (MFs) are mainly enriched in "oxidoreductase activity", "protein homodimerization activity", "kinase binding", "ATP-dependent activity", and "amide binding".The cellular component (CC) analysis revealed that the differentially expressed genes (DEGs) were predominantly associated with cellular components such as "blood microparticle", "mitochondrial matrix", "secretory granule Volcano and Venn diagrams of differentially expressed genes.lumen", "spindle", and "external side of the plasma membrane".KEGG pathway analysis shows that the DEGs are significantly involved in several key pathways, including the "pathway in cancer", "drug metabolism-cytochrome P450 pathway", "complement and coagulation cascades", "cell cycle", and "carbon metabolism" (Figure 4D).The results of GO analysis of DEG in pancreatic cancer are shown in Figure 4B, wherein the most important BPs include "regulation of cell adhesion", "positive regulation of cell migration", "tissue morphogenesis", "monocarboxylic acid metabolic process", and "extracellular matrix organization".The main enrichment of MFs is observed in "calcium ion binding", "extracellular matrix structural constituent", "lipid binding", "serine-type peptidase activity", and "oxidoreductase activity".The CC consists mainly of the "extracellular matrix", "apical part of the cell", "endoplasmic reticulum lumen", "perinuclear region of cytoplasm", and "lytic vacuole".KEGG pathway analysis reveals that the DEGs plays a crucial in the "protein digestion and absorption", "retinol metabolism", "ECM-receptor interaction", "transcriptional mis-regulation in cancer", and "complement and coagulation cascades" (Figure 4E).
The results of GO analysis of DEG in bladder cancer are shown in Figure 4C, and its most important BPs include "supramolecular fiber organization", "regulation of supramolecular fiber organization", "cellular response to cytokine stimulus", "positive regulation of cell death", and "regulation of MAPK cascade".The MFs are mainly enriched in "structural molecule activity", "actin binding", "lipid binding", "protein homodimerization activity", and "heparin-binding".The CC consists mainly of the "extracellular matrix", "contractile fiber", "polymeric cytoskeletal fiber", "focal adhesion", and "receptor complex".KEGG pathway analysis reveals that the DEGs plays a significant role in the "vascular smooth muscle contraction", "cardiac muscle contraction", "complement and coagulation cascades", and "malaria" (Figure 4F).

PPI network analysis
using the STRING database, a protein-protein interaction (PPI) network consisting of 714 conserved genes in hepatocellular carcinoma was constructed.The network comprised of 656 nodes and 1461 edges (Figure 5A).By considering the degree of connectivity in the PPI network, CytoNCA identified CDK1, TOP2A, ASPM, BUB1, NUSAP1, CDC20, DLGAP5, BUB1B, KIF20A, and CCNB2 as the top ten genes.These genes were tentatively regarded as potential key genes associated with hepatocellular carcinoma, emphasizing the need for further investigation (Figure 5B).
Similarly, a PPI network was established based on 357 conserved genes in pancreatic cancer, revealing a network with 86 nodes and 1133 edges (Figure 5C).Applying the degree of connectivity analysis, CytoNCA identified ASPM, TOP2A, DLGAP5, TPX2, CENPF, KIF23, MELK, LAMC2, LAMA3, and ANLN as the top ten genes.These genes were provisionally considered as potential key genes associated with pancreatic cancer, warranting further investigation (Figure 5D).
Furthermore, a PPI network was constructed based on 101 conserved genes in bladder cancer, consisting of 101 nodes and 131 edges (Figure 5E).Based on the degree of connectivity, CytoNCA identified TPM2, MYH11, MYL9, TAGLN, CNN1, ACTC1, FLNC, LMOD1, CSRP1, and VIM as the top ten genes.These genes were preliminarily regarded as potential key genes associated with bladder cancer, necessitating further investigation (Figure 5F).

Expression levels of key genes
The GEPIA2.0 database was used to verify the mRNA expression levels of key genes in tumor and normal tissues, and the results are shown in Figure 6.The expression levels of relevant potential key genes were markedly elevated in liver cancer and pancreatic cancer tissues compared to normal tissues (Figures 6A,  B); conversely, the expression levels of these genes were significantly decreased in bladder cancer tissues compared to normal tissues (Figure 6C).

Expression of key genes and their correlation with prognosis
To investigate whether the variation in mRNA expression levels of key genes affects the correlation between overall survival (OS) in cancer patients, further analysis was conducted.The GEPIA2.0 online platform was utilized to conduct Kaplan-Meier survival analysis, to explore the potential associations between key genes and overall survival (OS) in tumor patients.The analysis revealed that a majority of the genes exhibited a significant decrease in survival rate (p < 0.05, Figure 7).High levels of expression of CDK1 (p = 0.00022), TOP2A (p = 0.003), ASPM (p= 0.00072), BUB1 (p = 0.0012), NUSAP1 (p = 0.0067), CDC20 (p = 6.8E-06),DLGAP5 (p= 0.00049), BUB1B (p = 0.0031), and KIF20A (p = 0.0037) were found to be significantly associated with the OS in patients with liver cancer (p < 0.05); while CCNB2 (p = 0.053) expression was not relevant to survival, as shown in Figure 7A.

Key gene expressions are correlated with immune infiltration and immune cells
The study investigated the association between the expression of key genes and six types of infiltrating immune cells, namely B cells, CD8 + T cells, CD4 + T cells, macrophages, neutrophils, and dendritic cells, in order to determine their potential as independent predictors of cancer treatment and prognosis.

Molecular docking
Molecular docking can be used to explore the optimal binding mode between compounds and targets.Therefore, to further investigate interactions between Frondoside A with its target molecules, CDK1, ASPM, TOP2A, BUB1, CDC20, KIF20A, TPX2, KIF23, MELK, and FLNC molecules were used as molecular docking targets, and the molecular binding free energies listed in Table 2. Our results show that the interaction between Frondoside A saponins and key targets via hydrogenbonding indicated high stability and the binding free energy of the Frondoside A and the ten key proteins were all < -5 kcal/mol.Of these, the affinity by CDK1 performed the best (-10.7 kcal/mol).The 2D docking diagrams of Frondoside A with the receptor targets were visualized using PyMOL software as shown in Relationship between target gene expression and survival of tumor patients.(A-C) are the total survival time of hepatocellular carcinoma, pancreatic cancer, and bladder cancer, respectively.

Discussion
Conventional chemotherapeutic drugs currently in use have given rise to drug resistance due to diverse mechanisms of action.Therefore, new and effective therapeutic agents are urgently needed to treat cancer (14).The targeting of cancer growth, survival, migration, and metastasis pathways using drugs with minimal or no toxicity to normal cells is of significant importance (32).Frontside A has been shown to have anticancer activities in many cancer models, such as bladder cancer, lung cancer, breast cancer (33)(34)(35)(36).In the present study, Frondoside A has been demonstrated to exhibit concentration-dependent inhibitory effects on the viability of cancer cells.Frondoside A has been observed to significantly reduce cell migration of HepG2, Panc02, and UM-UC-3 cells in a time-dependent manner at IC 50 concentrations.Moreover, the inhibitory effects of Frondoside A on cell migration were found to be superior to those of EPI.It is important to note that the development of anti-tumor therapeutic agents that induce apoptosis and inhibit tumor invasion or metastasis are highly desirable (37).The effects of Frondoside A on cell viability or proliferation have been tested using multiple different methods in many different cancers (14,38,39).Frondoside A has demonstrated potent growth inhibitory effects on human pancreatic cancer cells, exhibiting an EC50 of approximately 1 µM (38).The observed inhibition of proliferation is accompanied by a significant increase in apoptosis (14,38).Frondoside A exhibits significant cytotoxicity against urothelial carcinoma cells, with IC 50 values ranging from 0.55 to 2.33 mM (33, 40).Consistent with our findings, Frondoside A, a marine-derived compound, has demonstrated potential as a therapeutic agent for the treatment of liver, pancreatic, and bladder cancers.While the precise mechanism underlying the anti-cancer effects of Frondoside A remains unclear (14), and the effect of Frondoside A on cell cycle and apoptosis of HepG2, Panc02, and UM-UC-3 require further investigation.
Recently, a study has reported that metformin can induce G2/M arrest and significantly inhibit the proliferation of HCC cells (41).Additionally, metformin has been found to effectively downregulate the expression of CDK1 (41).Several studies have reported significant upregulation of TOP2A mRNA and protein expression in hepatocellular carcinoma (HCC), suggesting that TOP2A is overexpressed in this cancer type and may serve as a potential biomarker for HCC (42).ASPM overexpression has been recognized as a molecular marker that correlates with heightened invasive and metastatic potential in hepatocellular carcinoma (HCC) (43,44).In contrast, upregulation of BUB1 significantly promoted cell proliferation, while downregulation of BUB1 expression inhibited the proliferation of liver cancer cell lines Analysis of key gene expression in various immune cells in liver cancer.The correlation between the expression of (A) CDK1, (B) TOP2A, (C) ASPM, (D) BUB1, (E) NUSAP1, (F) CDC20, (G) DLGAP5, (H) BUB1B, and (I) KIF20A and the degree of immune invasion in liver cancer was investigated using the TIMER database (http://timer.cistrome.org).The Rho value, which indicates Pearson's correlation coefficient, was used to evaluate the relationship between the genes and immune cells.When |Rho| > 0.1 and p < 0.05, it was considered that a correlation existed between the genes and immune cells.In general, the shape of the curve varied depending on the Rho value.When Rho < 0.5, the curve was elliptical; when Rho = 0.5, the curve was parabolic; and when Rho > 0.5, the curve was hyperbolic.These findings provide insights into the potential roles of these key genes in immune invasion in liver cancer.(45)(46)(47).Nucleolar and spindle-associated protein 1 (NUSAP1) is a member of the NUSP1 family of nucleolar-spindle-associated proteins and is involved in spindle microtubule organization (48).Li (49, 50) et al. studied the impact of CDC20 on the progression of hepatocellular carcinoma (HCC) and found that CDC20 expression was elevated in HCC samples.Transfection with CDC20 small interfering RNA in HCC cells led to reduced cellular proliferation and increased cell numbers in the G2/M phase.DLGAP5 (also known as HURP or KI-AA0008) is a cell-cycle-regulated protein that plays a role in tumor development (51).BUB1B, a crucial mitotic spindle checkpoint, is overexpressed in adrenocortical carcinomas (52) and promotes tumor proliferation while inducing radio resistance in glioblastoma (53,54).KIF20A accumulates in the nucleus during the G2 phase of the cell cycle and promotes both normal and pathological hepatocyte proliferation (55).ASPM-iII has been shown to selectively regulate cyclin E expression levels and cell cycle progression in PDAC cells (56).TOP2A expression is upregulated in pancreatic cancer tissues Analysis of key gene expression in various immune cells in pancreatic cancer.The relationship between the expression of (A) ASPM, (B) TOP2A, (C) DLGAP5, (D) TPX2, (E) KIF23, (F) MELK, and (G) ANLN and the degree of immune invasion in pancreatic cancer was analyzed using the TIMER database (http://timer.cistrome.org).The Rho value, which indicates Pearson's correlation coefficient, was used to evaluate the relationship between the genes and immune cells.When |Rho| > 0.1 and p < 0.05, it was considered that a correlation existed between the genes and immune cells.The shape of the curve varied depending on the Rho value, with smaller Rho values resulting in smoother curves, while larger Rho values resulted in fuller curves.Elliptical curves were observed when Rho < 0.5, parabolic curves when Rho = 0.5, and hyperbolic curves when Rho > 0.5.These findings provide insights into the potential roles of these key genes in immune invasion in pancreatic cancer.
compared to non-tumor tissues, and its upregulation is significantly associated with tumor metastasis and shorter survival in pancreatic cancer patients (57-59).Knockdown of TOP2A in pancreatic cancer cell lines inhibits cell proliferation and migration.DLGAP5 expression is significantly elevated in pancreatic cancer tissues and is correlated with patients' survival and progression-free survival (60).Knockdown of DLGAP5 inhibits pancreatic cancer cell proliferation, invasion, and migration (61).TPX2 has been identified as a prognostic biomarker of KARS-mutant PDAC (62).Knockdown of KIF23 can inhibit pancreatic cell proliferation (63).Silencing of MELK significantly reduces pancreatic cancer development (64), as MELK promotes CDK1 involvement in the cell cycle and cell progression in cancers (65).ANLN expression is significantly upregulated in pancreatic cancer tissues and cell lines and is associated with tumor size, differentiation, TNM stage, lymph node metastasis, distant metastasis, and poor prognosis in pancreatic cancer.ANLN knockdown inhibits several cell-cell adhesions-related genes, including the gene encoding LIM and SH3 protein 1 (LASP1).LASP1 upregulation partially reverses the tumor-suppressive effect of ANLN downregulation on pancreatic cancer cell progression.ANLN contributes to pancreatic cancer progression by regulating the EZH2/miR-218-5p/LASP1 signaling axis (66, 67).
Consistent with our study, previous studies have reported these key genes (TAGLN, CNN1, ACTC1, LMOD1) play important roles in bladder cancer.TAGLN is a recognized actin-binding protein that modulates the dynamics of the actin cytoskeleton (68-70).CNN1 was significantly lowly expressed in BC tissues and cells (70).ACTC1 is responsible for encoding cardiac actin, and a c.G301A mutation in this gene has been shown to be associated with Analysis of key gene expression in various immune cells in bladder cancer.The relationship between the expression of (A) TAGLN, (B) CNN1, (C) ACTC1, (D) FLNC, and (E) LMOD1 and the degree of immune invasion in bladder cancer was analyzed using the TIMER database (http://timer.cistrome.org).The Rho value, which indicates Pearson's correlation coefficient, was used to evaluate the relationship between the genes and immune cells.When |Rho| > 0.1 and p < 0.05, it was considered that a correlation existed between the genes and immune cells.The shape of the curve varied depending on the Rho value, with smaller Rho values resulting in smoother curves, while larger Rho values resulted in fuller curves.Elliptical curves were observed when Rho < 0.5, parabolic curves when Rho = 0.5, and hyperbolic curves when Rho > 0.5.These findings provide insights into the potential roles of these key genes in immune invasion in bladder cancer.hypertrophic cardiomyopathy, which, in some instances, leads to sudden cardiac death (71).Previous studies that ACTC1 is related to cancer prognosis (72, 73), ACTC1 expression is significantly upregulated in glioblastoma and inhibits migration of cancer cells (74), and ACTC1 correlates with the prognosis of glioblastoma and can be used as a novel prognostic marker in glioma (72).LMOD1 belongs to the LMOD family of proteins, which exhibit a striking similarity to actin-capping proteins referred to as Tropomodulins (TMODs) (75).LMOD1, also known as Leiomodin 1, displays widespread expression in a majority of tissues, with particularly elevated levels in the thyroid, skeletal muscle, eye muscle, and ovary (76).Dysregulated expression of LMOD1 may potentially be linked to various pathological conditions (77).LMOD1 was found to be a new gastric cancer biomarker and therapeutic target that induces EMT by regulating the FAK-Akt/mTOR pathway (78).
Molecular docking was used in this study to investigate the mechanism and provide valuable guidance for drug screening and design in future experiments.Based on our molecular investigations, the free energies of Frondoside A and the critical proteins were observed to be below -5 kcal/mol, indicating favorable binding of the aforementioned active ingredients to their respective targets.Studies have shown that Frondoside A can down-regulate CDC20 gene expression, which inhibits the proliferation and migration of cancer cells and promotes cycle arrest and cell apoptosis (79,80).In addition, Zhang and colleagues have shown that TOP2A deletion can significantly inhibit the proliferation and migration of cancer cells and induce cell apoptosis (81-83).Furthermore, down-regulated expression of CDK1 not only inhibits proliferation and migration of hepatocellular carcinoma cells, but also induces cycle arrest and specific apoptosis of hepatocellular carcinoma cells (84, 85).Finally, overexpression of CNN1 has been shown to inhibit the proliferation, invasion and metastasis of bladder cancer cells (86).Therefore, it is likely that that Frondoside A inhibits the proliferation and metastasis of cancer cells by down-regulating CDK1, CDC20, TOP2A and upregulating CNN1.It should be noted that the free binding energy  between Frondoside A and CDK1 was the lowest, -10.7 kcal/mol.This indicates that Frondoside A may inhibit the development of liver cancer through its effect on the cell cycle (39).The specific mechanisms need to be further investigated.
In addition, Frondoside A may affect immune infiltration by modulating certain targets.These targets may involve activation, recruitment, or inhibition of immune cells.For example, Frondoside A may inhibit the expression of certain proinflammatory cytokines in order to reduce the inflammatory response and immune cell infiltration.Alternatively, it may act directly on immune cells, such as T cells or macrophages, affecting their activation and function.This effect may help to improve the tumor microenvironment and enhance anti-tumor effects.However, further studies are needed to validate this hypothesis and determine the specific molecular mechanisms.

Conclusions
The integration of bioinformatics and molecular docking analyses has enhanced our comprehension of the underlying molecular mechanisms through which Frondoside A exerts its effects in hepatocellular carcinoma, pancreatic cancer, and bladder cancer.Our findings suggest that Frondoside A holds great promise as a potential therapeutic candidate for the treatment of these cancers due to its multi-targeted and multimechanistic approach.
FIGURE 2 Inhibition of tumor cell proliferation and migration by Frondoside A. (A-C) Logarithmic growth HepG2, Panc02, and UM-UC-3 cells were treated with EPI and Frondoside A for 48 h.(D-F) Figures of cell migration of HepG2, Panc02, and UM-UC-3 respectively were confluent mono-layers cultured in the presence or absence (control) with the IC50 concentrations of Frondoside A for 0, 12, and 24 h, and EPI (10 mM) as a positive group.(G-I) The migration rate of HepG2, Panc02, and UM-UC-3, respectively.Determination of cells as described in Materials and methods.All experiments were repeated at least three times.Data for SD ± means, *p < 0.05, * *p < 0.01, * * * p < 0.001, * * * *p < 0.0001.

4
FIGURE 4    GO and KEGG enrichment analysis of DEGs.(A-C) GO enrichment in liver cancer, pancreatic cancer, and bladder cancer, respectively; and (D-F) KEGG pathway analysis of liver, pancreatic, and bladder cancers, respectively.

5 PPI
FIGURE 5 PPI and key target screening of differentially expressed genes.(A, C, E) are the PPI figures for liver cancer, pancreatic cancer, and bladder cancer, respectively; (B, D, F) are the key targets of liver cancer, pancreatic cancer, and bladder cancer.

Figure 11 . 6
FIGURE 6 Differential expression analysis of key genes in tumors in the GEPIA2.0database.(A) Liver cancer tissue; (B) Pancreatic cancer tissue; and (C) Bladder cancer tissue.

TABLE 1
Analysis of key gene expression in various immune cells in liver, pancreatic, and bladder cancers.

TABLE 2
Binding Free energy of Frondoside A and targets.