Centipeda minima active components and mechanisms in lung cancer
BMC Complementary Medicine and Therapies volume 23, Article number: 89 (2023)
Traditional Chinese medicine (TCM) has been extensively used for neoplasm treatment and has provided many promising therapeutic candidates. We previously found that Centipeda minima (C. minima), a Chinese medicinal herb, showed anti-cancer effects in lung cancer. However, the active components and underlying mechanisms remain unclear. In this study, we used network pharmacology to evaluate C. minima active compounds and molecular mechanisms in lung cancer.
We screened the TCMSP database for bioactive compounds and their corresponding potential targets. Lung cancer-associated targets were collected from Genecards, OMIM, and Drugbank databases. We then established a drug-ingredients-gene symbols-disease (D-I-G-D) network and a protein–protein interaction (PPI) network using Cytoscape software, and we performed Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses using R software. To verify the network pharmacology results, we then performed survival analysis, molecular docking analysis, as well as in vitro and in vivo experiments.
We identified a total of 21 C. minima bioactive compounds and 179 corresponding targets. We screened 804 targets related to lung cancer, 60 of which overlapped with C. minima. The top three candidate ingredients identified by D-I-G-D network analysis were quercetin, nobiletin, and beta-sitosterol. PPI network and core target analyses suggested that TP53, AKT1, and MYC are potential therapeutic targets. Moreover, molecular docking analysis confirmed that quercetin, nobiletin, and beta-sitosterol, combined well with TP53, AKT1, and MYC respectively. In vitro experiments verified that quercetin induced non-small cell lung cancer (NSCLC) cell death in a dose-dependent manner. GO and KEGG analyses found 1771 enriched GO terms and 144 enriched KEGG pathways, including a variety of cancer related pathways, the IL-17 signaling pathway, the platinum drug resistance pathway, and apoptosis pathways. Our in vivo experimental results confirmed that a C. minima ethanol extract (ECM) enhanced cisplatin (CDDP) induced cell apoptosis in NSCLC xenografts.
This study revealed the key C. minima active ingredients and molecular mechanisms in the treatment of lung cancer, providing a molecular basis for further C. minima therapeutic investigation.
Lung cancer is one of the most common malignancies and the leading cause of all cancer-related deaths worldwide [1, 2]. Non-small cell lung cancer (NSCLC) represents approximately 85% of lung cancer cases. Current treatments for NSCLC mainly include surgery, radiotherapy, chemotherapy, targeted therapy, or some combination of these [3, 4]; however, the therapeutic efficacy of available treatments is not yet satisfactory. The majority of NSCLC patients have advanced disease at diagnosis and are not candidates for surgery [5, 6] and chemotherapy and radiotherapy effectiveness is frequently limited due to the onset of drug resistance and adverse effects. Thus, alternative therapeutic strategies are urgently needed for NSCLC treatment.
Many natural products and derivatives, such as Vinca alkaloids, taxanes, and camptothecins, have become part of the standard neoplasm chemotherapy repertoire . Traditional Chinese medicine (TCM) has been used in cancer treatment for thousands of years, providing a rich basis for anti-carcinogen discovery and development. TCM products have been shown to exert anti-cancer effects through multiple signaling pathways and molecular targets . Centipeda minima (L.)A.Br.et Aschers (C. minima), also called Ebushicao in TCM, belongs to the daisy family Asteraceae , is native to the tropical regions, and is famous for its extraordinary efficacy in treating nasal allergy, whooping cough, asthma, and malaria in TCM . Phytochemical studies revealed that C. minima chemical constituents mainly include terpenoids, flavonoids, sterols, saponins, monophenols, fatty acids, and amides . Crude extracts and purified compounds from C. minima have been shown to exert promising anti-cancer effects in a variety of human malignancies such as lung cancer, multiple myelomanasopharyngeal carcinoma, breast cancer, and colon cancer [12,13,14,15]. Nevertheless, the C. minima ingredients with anti-cancer activities are poorly characterized. Our previous studies found that a C. minima ethanol extract (ECM) combined with cisplatin (CDDP) had synergistic anti-cancer effects both in vitro and in vivo ; however, the underlying molecular mechanisms remain ill-defined.
Here we applied network pharmacology and molecular docking to explore C. minima active ingredients and mechanisms against NSCLC. Network pharmacology, based on the general ideas of systems biology, is considered to be an efficient method for identifying drug mechanisms of action, and molecular docking is a powerful tool used in structural molecular biology and in computer-aided drug design . Furthermore, we performed in vitro and in vivo experiments to verify the obtained results (Fig. 1).
Materials and methods
Screening for C. minima active components and targets
The Traditional Chinese Medicine Systems Pharmacy Database and Analysis Platform (TCMSP, http://lsp.nwu.edu.cn/tcmsp.php) is a unique systems pharmacology platform for Chinese herbal medicines that captures relationships among drugs, targets, and diseases. All ingredients were annotated using the Mol ID which is a unique numbering system for the TCMSP database. To identify C. minima compounds that may influence NSCLC, we preliminarily screened compounds collected from the TCMSP database using oral bioavailability (OB) and drug-likeness (DL) as parameters. OB is an important pharmacokinetic parameter used to measure a drug’s ADME (absorption, distribution, metabolism, excretion) in vivo. DL refers to the structural similarity between herbal ingredients and known drugs, and it plays a major role in new drug development. In this study we selected compounds with OB ≥ 30% and DL ≥ 0.18 as candidate active ingredients for further analysis. We searched for active compound standard names in the PubChem database (https://pubchem.ncbi.nlm.nih.gov/). Chemical constituents’ protein targets were obtained from the TCMSP database with OB ≥ 30% and DL ≥ 0.18 as screening conditions. Chemical ingredient and target databases were constructed using the data extracted from the TCMSP database.
Lung cancer related targets were obtained from the Online Mendelian Inheritance in Man database (OMIM, http://omim.org/), Genecards database (https://www.genecards.org/), and Drugbank database (http://www.drugbank.ca/) using “lung cancer” as the keyword. The chemical constituents’ targets collected above were subsequently mapped to the lung cancer related targets, and the candidate targets related to both C. minima and lung cancer were obtained using VENNY2.1 (https://bioinfogp.cnb.csic.es/ tools/ venny/). All targets were finally identified and normalized using the UniProt database (https://www.uniprot.org/).
Drug-ingredients-gene symbols-disease (D-I-G-D) network construction
Cytoscape  is an open source software platform that visualizes complex networks by integrating any type of attribute information. We imported the identified candidate targets, bioactive compounds, diseases, and herbs into Cytoscape 3.7.2 to construct a D-I-G-D network. In this network, nodes indicate drugs, active ingredients, targets, and diseases, and the edges indicate the interactions between them. The network analyzer Cytoscape plug-in was applied to obtain network topology parameters, such as nodes, degrees, and edge betweenness.
Protein–protein interaction (PPI) network construction
To obtain protein–protein interaction relationships, candidate target gene data were imported into the STRING database (https://string-db.org/)  with “Homo sapiens” as the selected species and the minimum required interaction score limited to a medium confidence score of 0.400. The results were downloaded in TSV format and imported into Cytoscape 3.7.2 to visualize the PPI network interactions. The CytoHubba Cytoscape plugin was used to identify hub genes in the PPI network.
We input candidate targets into R4.0.5 software to perform GO and KEGG enrichment analyses. P-value < 0.05 and Q-value < 0.05 were considered statistically significant. The GO enrichment analysis included cellular components (CC), biological processes (BP), and molecular functions (MF). The top twenty significant terms for each category are presented, and bubble diagrams were drawn using R4.0.5 software. KEGG enrichment analysis clarified C. minima potential mechanisms against lung cancer, and we constructed bubble charts depicting the top 20 significant pathways .
NSCLC gene expression dataset analysis
We used the Kaplan–Meier plotter (http://www.kmplot.com) to determine the correlation between hub gene expression and NSCLC patient prognosis.
Component-target molecular docking
Molecular docking is a method used to study intermolecular interactions and predict their binding modes and affinity . Here, we used it to confirm interactions between C. minima core compounds and hub targets in treating lung cancer and to verify the network pharmacology predictions’ accuracy. The specific process was as follows:
2.7.1 Small ligand molecule file preparation: Three-dimensional (3D) files for C. minima were downloaded in SDF format from the PubChem database and converted to PDB format using OpenBabel. They were saved as PDBQT ligand files using AutoDock Tools 1.5.6 software.
2.7.2 Macromolecular receptor file preparation: Receptor 3D structures were obtained from the Protein Data Bank (PDB, http://www.rcsb.org/, in PDB file format). Pymol 2.5.2 software was then used to extract the target proteins’ original ligand conformation by removing water molecules, co-crystallized ligand, and ions. Subsequently, AutoDock tools 1.5.6 was used to add polar hydrogen and Gasteiger charge to the processed receptors’ structures. The structures were then saved as PDBQT protein receptor files.
2.7.3 Defining docking parameters: The receptor and ligand PDBQT structures were imported into Autodock Tools 1.5.6 and the molecular docking range was defined. The docking center was set on the macromolecule, and the center coordinate (center x/y/z) and box size (size x/y/z) parameters were set to keep the protein completely covered by the docking box.
2.7.4 Docking and visualization: Autodock Vina 1.1.2 was used to conduct molecular docking and calculate docking affinity. The molecular docking results analysis refers to the binding energy (ΔGbind). A binding energy less than -5 kJ∙mol−1 indicates that the target has certain binding activity with the compound . The lower the binding energy, the better the docking effect. Finally, the compound and protein docking results were analyzed and visualized using Pymol 2.5.2.
Chemicals and reagents
CDDP was purchased from Sigma Aldrich (St. Louis, MO, USA). Primary antibodies against Bcl-2 and GAPDH and all secondary antibodies were purchased from Cell Signaling Technology (Beverly, MA, USA). High purity (98.02%) quercetin was obtained from MedChemExpress (Monmouth Junction, NJ, USA).
A549 and H1299 cells, two commonly used NSCLC cell lines, were purchased from the American Type Culture Collection (ATCC, Manassas, VA, USA). A549 cells were cultured in Dulbecco’s Modified Eagle Medium (DMEM) supplemented with 10% fetal bovine serum (FBS). H1299 cells were maintained in RPMI-1640 medium containing 10% FBS. All reagents listed above were purchased from Gibco (CA, USA). Cells were cultured in a humidified incubator with 5% CO2 at 37 °C.
Cell viability assay
A549 and H1299 cells were seeded in 96-well plates at a density of 4000 cells per well, incubated for 24 h, and then treated with quercetin at the indicated concentrations for 48 h. The cells were then incubated with Cell Counting Kit-8 (CCK-8) solution for an additional 1 h at 37 °C and the absorbance was measured at 450 nm using a microplate reader (ThermoFisher, Waltham, MA).
ECM powder was prepared using 95% ethanol as previously described . ECM was dissolved in ethanol to make a stock solution of 200 mg/ml and then diluted in corn oil to final concentrations of 20 and 40 mg/ml.
In vivo animal experiments
BALB/c nude mice were obtained from Guangdong Medical Laboratory Animal Center. All mice were maintained in the specific pathogen free grade laboratory, and were subcutaneously inoculated with A549 cells (1 × 107 cells mixed with 50 μl PBS and 50 μl Matrigel) into the right flank. 14 days after subcutaneous injection, the mice were randomly divided into 5 groups (seven mice per group): (1) vehicle control (solvent: 10% ethanol in normal saline); (2) ECM group (200 mg/kg ECM gavaged daily); (3) CDDP group (2 mg/kg CDDP injected i.p. once every three days); (4) CDDP + ECM-L (100 mg/kg); (5) CDDP + ECM-H (200 mg/kg). Mice were sacrificed at the end of the experiment by decapitation, and the tumors were collected for western blotting and hematoxylin and eosin (HE) staining.
Samples were lysed in RIPA lysis buffer. The bicinchoninic acid assay (BCA assay) was performed to quantify protein concentrations. Equal amounts of protein were separated by SDS-PAGE using 10% polyacrylamide. Next, proteins were transferred to PVDF membranes. Membranes were blocked with 5% fat-free milk, cut into strips, and incubated with primary antibodies and secondary antibodies. Detection was performed using an enhanced chemiluminescent kit (Millipore, MA, USA). Images were acquired using a ChemiDoc XRS system with Quantity One software (Bio-Rad, Richmond, CA, USA).
The tumor tissue was collected and fixed in 4% paraformaldehyde. Then the samples were dehydrated with an alcohol concentration gradient and xylenes and then embedded in paraffin. For HE experiments, Sects. (4 μm) were stained with hematoxylin and eosin, and the slides were observed and photographed under a light microscope.
Statistical analyses were performed using SPSS software 17.0 and GraphPad Prism 7.0. Data shown is from at least three independent experiments (mean ± standard error of the mean (SEM)). To determine differences, unpaired t-tests and one-way ANOVA were performed. P-values < 0.05 were considered statistically significant.
Screening for C. minima active components and related targets
We obtained a total of 108 C. minima compounds from the TCMSP database. After screening for bioactive compounds using specific OB and DL conditions, we selected 21 active ingredients for further analysis, including sesquiterpenoids (plenolin, arnicolide D, arnicolide C, florilenalin angelic acid, florilenalin isobutyrate, senecioylplenolin, helenalin, 2beta-hydroxyl-2,3-dihydrogen-6-O-angeloplenolin, and 2,3,11,13-tetrahydrohelenalin), triterpenoids (taraxasteryl acetate, taraxerol, and (18alpha,19alpha)-5alpha-Urs-20(30)-ene-3beta-ol palmitate), steroids (sitosterol, stigmasterol, and beta-sitosterol), flavonoids (quercetin der., nobiletin, and quercetin), amino and amide derivatives (aurantiamide acetate), and phenols (8-hydroxy-9,10-diisobutyryloxythymol and (Z)-3,3',5,5'-Tetramethoxystilbene). Detailed information about these compounds is listed in Table 1.
Potential target identification
After collecting potential targets from the TCMSP database, converting them into the UniProt database, and deleting redundant items, we obtained 179 known target symbols associated with C. minima. We also screened a total of 804 genes related to lung cancer from the OMIM, Genecards, and Drugbank databases. We obtained a total of 60 overlapping genes (Fig. 2, see Venn diagram) and used these as C. minima candidate targets against lung cancer. Some of the target information is shown in Table 3.
D-I-G-D network analyses
The D–I–G–D network contained 73 nodes (60 genes, 11 chemicals, 1 drug, and 1 disease) and 163 edges. In this network, the blue and red nodes represent C. minima and lung cancer, respectively. The 60 purple nodes represent the overlapping genes between C. minima and lung cancer, and the 11 pink nodes indicate active ingredients in C. minima (10 active ingredients with no common intersection targets were deleted). In addition, the degree of a node was evaluated by the proportion of edges connected to the node, which represents its importance in the network (Fig. 3). Quercetin, nobiletin and beta-sitosterol had the top three degree value, thus we considered them be the major potential C. minima compounds for lung cancer treatment (Table 2).
PPI network analyses
The constructed PPI network contained 60 nodes and 961 edges (Fig. 4A). Nodes and edges represent proteins and the connections between two proteins, respectively. We visualized the String result predictions using Cytoscape software (Fig. 4B). The node colors (light to dark) and sizes (small to large) are proportional to the node degrees. Nodes with a higher degree are considered to be more important in the network. The top 30 identified hub genes with the highest degree of connectivity were: TP53, AKT1, MYC, CASP3, EGFR, CCND1, VEGFA, JUN, ESR1, MAPK8, MAPK1, ERBB2, EGF, PTGS2, CXCL8, TNF, MMP9, BCL2L1, AR, MMP2, CASP8, IL1B, CDKN1A, PPARG, CCL2, RELA, HIF1A, PGR, CASP9, and SPP1. Among these genes, TP53 (degree = 59), AKT1 (degree = 53) and MYC (degree = 52) showed the three highest node degrees (Fig. 4C and Table 3). Survival analysis indicated that expression of top ten genes were closely linked to the lung cancer prognosis (Fig. 5), which was consistent with the results above.
GO function enrichment analysis
We performed three categories of GO functional annotation analysis on the 60 potential targets, and the top 20 significant terms from each category are shown in Fig. 5. Overall, we identified 1771 GO functions, including 1630 BP functions, 30 CC functions, and 111 MF functions. In the CC category, organelle outer membrane, outer membrane, membrane raft, membrane microdomain, and membrane region were enriched (Fig. 6A). In the BP category, cellular response to chemical stress, epithelial cell proliferation, regulation of apoptotic signaling pathway, response to oxidative stress, and response to UV were enriched (Fig. 6B). In the MF category, ubiquitin-like protein ligase binding, DNA-binding transcription factor binding, ubiquitin protein ligase binding, RNA polymerase II-specific DNA-binding transcription factor binding, and nuclear receptor activity were enriched (Fig. 6C).
KEGG pathway enrichment analysis
We identified a total of 144 pathways significantly associated with target genes. The top 20 significant pathways included cancer pathways, platinum drug resistance, apoptosis, and IL-17 signaling pathway (Fig. 7, Table 4). These pathways may be involved in the molecular mechanisms underlying C. minima’s impact on lung cancer.
Compound-target interaction validation by molecular docking
To further confirm the binding capability between active drugs and important targets, we performed molecular docking. The three top targets with the highest degrees from the PPI network (TP53 (PDB ID:7BWN), AKT1 (PDB ID:4EJN), and MYC (PDB ID:7T1Z)) were considered key targets. Quercetin (PubChem CID:5,280,343), nobiletin (PubChem CID:72,344) and beta-sitosterol (PubChem CID:222,284), which had the highest degree in the active components target network, were considered major active compounds. At the same time, Erlotinib was used to be docked with the three proteins as a positive drug . The details of docking results were shown in Fig. 8, and Table 5. All of the binding energy values were smaller than -5 kJ∙mol − 1, suggesting that they have good binding activity. Erlotinib was chosen as positive control, and the binding energies of quercetin, nobiletin and beta-sitosterol were lower than that of Erlotinib, indicating that the binding stability of quercetin, nobiletin and beta-sitosterol was better than that of Erlotinib. Besides, we found that quercetin shows the lowest binding affinity. Molecular docking analysis showed that quercetin mainly forms six hydrogen bonds with residues THR-97, GLU-95, LYS-107, THR-246 and GLN-248 on TP53 protein (Fig. 8A), five hydrogen bonds with residues THR-211, ASP-292, GLN-97, and ASN-54 on AKT1 protein (Fig. 8D), and eleven hydrogen bonds with residues SER-427, TRP-425, ARG-465, ARG-505, MET-587, and GLN-548 on MYC protein (Fig. 8G).
Quercetin induced NSCLC cell death
To investigate how quercetin from C. minima impacts lung cancer cells, we performed a cell viability assay. We found that quercetin treatment reduced NSCLC cell viability in a dose-dependent manner (Fig. 9A). After quercetin exposure, the number of adherent cells decreased, and cells shrank and became round (Fig. 9B). These results demonstrated that quercetin had significant cytotoxic effects on NSCLC cells, which was consistent with the network pharmacology results.
ECM enhances CDDP-induced cell apoptosis in NSCLC xenografts
KEGG pathway enrichment analysis suggested that the platinum drug resistance pathway and apoptosis pathways may be involved in C. minima’s effects on lung cancer cells. To validate these results, we performed in vivo animal studies. CDDP as a most common platinum drug, has been widely used in the first-line treatment of lung cancer. We found that CDDP or ECM administration alone slightly induced apoptosis. However, ECM and CDDP co-treatment robustly increased the apoptosis level (Fig. 10A). Bcl2 belongs to the anti-apoptotic Bcl2 family and is required for cancer cells to survive . Consistently, combined ECM and CDDP treatment reduced Bcl2 expression compared with CDDP treatment alone (Fig. 10B, C). These results demonstrated that ECM enhanced CDDP’s anti-cancer effects, at least partly, due to apoptosis pathway activation, consistent with the network pharmacology results.
TCM is a complementary or alternative medical system used to treat various diseases and is also a potential source of candidate resources for cancer treatment. C. minima, a Chinese herb commonly used in clinical practice, has been shown to have anti-cancer [10, 27], anti-proliferative [15, 28], antioxidant , anti-inflammatory [30,31,32] and antibacterial activities [33, 34]. In our previous study, we confirmed the herb’s anti-cancer effects in NSCLC both in vitro and in vivo , but the underlying mechanism was unclear. The basic network pharmacology theory is consistent with TCM’s holistic principals and is well suited for analyzing multi-targeted Chinese medicine [35,36,37]. Therefore, we used a network pharmacology approach that combined screening for active components, drug targeting, network analysis, and pathway analysis to explore the therapeutic mechanisms underlying lung cancer treatment with C. minima.
We screened 21 potential compounds from C. minima and, eventually, 11 compounds were included in our D-I-G-D network, including sesquiterpenoids (arnicolide D and 2beta-hydroxyl-2,3-dihydrogen-6-O-angeloplenolin), steroids (sitosterol, stigmasterol, and beta-sitosterol), flavonoids (quercetin der., nobiletin, and quercetin), amino and amide derivatives (aurantiamide acetate), and phenols (8-hydroxy-9, 10-diisobutyryloxy-thymol, and 1-[(Z)-2-(3,5-dimethoxyphenyl)vinyl]-3,5-dimethoxybenzene). Our results indicated that these 11 compounds may be related to C. minima’s anti-cancer effects. According to the 11 compounds’ D-I-G-D network degrees and molecular docking analysis, quercetin emerged as a key C. minima active compound against lung cancer. Quercetin is a member of the flavonoid subclass of flavonols and has been found to exhibit significant anti-tumor properties . Quercetin has been shown to inhibit proliferation in several cancer cell lines [39,40,41], and one study found that quercetin had anti-tumor activity in NSCLCs harboring wild-type EGFR by inhibiting AXL and inducing apoptosis . Another study found that quercetin induced pro-apoptotic autophagy in A549 and H1299 cells . Our in vitro experiment confirmed quercetin’s anti-tumor effects in NSCLC.
Further, our results suggest that TP53, AKT1, and MYC may be key targets for C. minima active compound anti-lung cancer activity. TP53 encodes p53 tumor suppressor protein and acts as a critical failsafe mechanism for cellular anti-cancer defenses [44,45,46]. TP53 has been implicated in a number of fundamental biological processes including DNA repair, cell-cycle arrest, apoptosis, autophagy, senescence, metabolism, and aging [47,48,49,50]. A TP53 tumor suppressor gene mutation occurs in approximately 50% of NSCLC cases . Alpha subunit protein kinase B (AKT1), a member of the serine/threonine protein kinase family, has been shown to have a role in a variety of tumors, such as prostate cancer, osteosarcoma, ovarian cancer, and endometrial cancer, by inducing apoptosis, inhibiting tumor cell proliferation, and decreasing tumor cell invasion and metastatic capacity [52,53,54,55]. AKT1 inhibition has been reported to promote migration and invasion of NSCLC cells with KRAS or EGFR mutations in vitro . Wu et al. found that miR-377-5p can inhibit lung cancer cell proliferation, invasion, and cell cycle progression by targeting AKT1 signaling . Thus, AKT1 has emerged as an attractive target for new anti-cancer therapy development. MYC consists of 3 paralogs, C-MYC, N-MYC, and L-MYC, that are frequently deregulated driver genes in human cancers . A previous study found that MYC-induced long non-coding RNAs inhibited NSCLC cell cycle arrest and apoptosis by activating C-MYC and its downstream effectors . In our study, we used AutoDock Vina and PyMOL software to investigate quercetin’s possible binding sites with related proteins and found that quercetin had strong binding activity with TP53, AKT1, and MYC.
CDDP has been widely used in first-line lung cancer treatment . However, the development of chemoresistance remains a major obstacle for treatment success . We used GO enrichment analysis to annotate protein target functions in three main categories: BP, CC, and MF. Our KEGG pathway analysis revealed that these targets were significantly enriched in a variety of cancer pathways, platinum drug resistance, and apoptosis pathways, interesting areas for further research. Notably, recent research found that quercetin could overcome CDDP resistance in cervical cancer cells . Additionally, our previous research found that ECM can sensitize NSCLC cells to CDDP and hypersensitize CDDP-induced apoptosis in NSCLC, but the underlying mechanism was unclear . In the present study, western blotting and HE staining showed that ECM could significantly enhance CDDP-induced apoptosis. All of these experimental results validated our network pharmacology prediction results.
Based on our network pharmacology results and experimental evidence from previous work, C. minima clearly has the potential for use in developing new NSCLC treatments. Additionally, our study demonstrated that network pharmacology predictions can be used to characterize TCM pharmacological mechanisms in detail. Evidence suggests that C. minima functions as an anti-cancer drug through multiple targets, amongst which inducing apoptosis and influencing platinum drug resistance were pathways identified by both in vivo and in vitro experiments. However, the current work has several limitations. First, more comprehensive TCM target gene databases are needed to obtain more reliable results. Furthermore, it is necessary to explore how the three hub genes change in response to single-drug therapy or two-drug combination therapy. Finally, although we determined that quercetin was the most important bioactive compound in C. minima, it may not be completely interchangeable for C. minima. Thus, further studies investigating C. minima’s anti-lung cancer molecular mechanisms are needed.
Here we used network pharmacology along with molecular docking and experimental validation to explore C. minima’s anti-lung cancer material basis and mechanism of action. We found that quercetin was the key active ingredient and that TP53, AKT1, and MYC may be three potential C. minima targets in lung cancer treatment. Additionally, we found that C. minima may exert anti-cancer effects by inducing apoptosis and influencing platinum drug resistance. In conclusion, our research provides a reference for further study of C. minima’s potential clinical application for lung cancer treatment.
Availability of data and materials
The datasets used and/or analyzed in the current study are available from the corresponding author upon reasonable request.
- C. minima :
Centipeda minima (L.)A.Br.et Aschers
Non-small cell lung cancer
Ethanol extract of C. minima
Traditional Chinese medicine
Traditional Chinese Medicine Systems Pharmacy database and analysis platform
Online Mendelian Inheritance in Man
Kyoto Encyclopedia of Genes and Genomes
Cellular tumor antigen p53
RAC-alpha serine/threonine-protein kinase
V-myc myelocytomatosis viral oncogene homolog
Sung H, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. Ca-a Cancer J Clin. 2021;71(3):209–49.
Bade BC, Dela Cruz CS. Lung cancer 2020 epidemiology, etiology, and prevention. Clin Chest Med. 2020;41(1):1.
Jones GS, Baldwin DR. Recent advances in the management of lung cancer. Clin Med. 2018;18(2):S41–6.
Herbst RS, Morgensztern D, Boshoff C. The biology and management of non-small cell lung cancer. Nature. 2018;553(7689):446–54.
Alexander M, Kim SY, Cheng HY. Update 2020: management of non-small cell lung cancer. Lung. 2020;198(6):897–907.
Zang FR, et al. Shikonin suppresses NEAT1 and Akt signaling in treating paclitaxel-resistant non-small cell of lung cancer. Mole Med. 2020;26(1):7.
Efferth T, et al. From traditional Chinese medicine to rational cancer therapy. Trends Mole Med. 2007;13(8):353–61.
Yan ZK, Lai ZJ, Lin JM. Anticancer properties of traditional Chinese medicine. Combinatorial Chem High Throughput Screen. 2017;20(5):423–9.
Guan YM, Zhang NS, Zhang YY. First report of fusarium commune causing root and basal stem rot on centipeda minima in China. Plant Dis. 2016;100(12):2525–6.
Linh NT, et al. Medicinal plant centipeda minima: a resource of bioactive compounds. Mini-Rev Med Chem. 2021;21(3):273–87.
Tan JC, et al. Centipeda minima: an update on its phytochemistry, pharmacology and safety. J Ethnopharmacol. 2022;292:115027.
Wang Y, et al. Inhibition of Nrf2 enhances the anticancer effect of 6-O-angeloylenolin in lung adenocarcinoma. Biochem Pharmacol. 2017;129:43–53.
Liu R, et al. Brevilin A induces cell cycle arrest and apoptosis in nasopharyngeal carcinoma. Front Pharmacol. 2019;10:594.
Lee MML, et al. Anti-cancer activity of centipeda minima extract in triple negative breast cancer via Inhibition of AKT, NF-kappa B, and STAT3 signaling pathways. Front Oncol. 2020;10:491.
You PT, et al. Brevilin A induces apoptosis and autophagy of colon adenocarcinoma cell CT26 via mitochondrial pathway and PI3K/AKT/mTOR inactivation. Biomed Pharmacother. 2018;98:619–25.
Fan XZ, et al. Centipeda minima extract sensitizes lung cancer cells to DNA-crosslinking agents via targeting Fanconi anemia pathway. Phytomedicine. 2021;91:153689.
Saikia S, Bordoloi M. Molecular docking: challenges, advances and its use in drug discovery perspective. Curr Drug Targets. 2019;20(5):501–21.
Shannon P, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
Hsia CW, et al. Analysis of dermal papilla cell interactome using STRING database to profile the ex vivo hair growth inhibition effect of a Vinca alkaloid drug, colchicine. Int J Mol Sci. 2015;16(2):3579–98.
Kanehisa M, et al. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 2023;51:D587.
Chen GL, Seukep AJ, Guo MQ. Recent advances in molecular docking for the research and discovery of potential marine drugs. Mar Drugs. 2020;18(11):545.
Gaillard T. Evaluation of autodock and autoDock Vina on the CASF-2013 benchmark. J Chem Inform Model. 2018;58(8):1697–706.
Wang YJ, et al. Ethanol extract of centipeda minima exerts antioxidant and neuroprotective effects via activation of the Nrf2 signaling pathway. Oxid Med Cell Longev. 2019;2019:9421037.
Midha A, Dearden S, McCormack R. EGFR mutation incidence in non-small-cell lung cancer of adenocarcinoma histology: a systematic review and global map by ethnicity (mutMapII). Am J Cancer Res. 2015;5(9):2892.
Vogler M. Targeting BCL2-proteins for the treatment of solid tumours. Adv Med. 2014;2014:943648.
Zhou W, et al. Systems pharmacology exploration of botanic drug pairs reveals the mechanism for treating different diseases. Sci Rep. 2016;6:36985.
Liu R, et al. Arnicolide D, from the herb centipeda minima, is a therapeutic candidate against nasopharyngeal carcinoma. Molecules. 2019;24(10):1908.
Su MX, et al. Antiproliferative effects of volatile oils from centipeda minima on human nasopharyngeal cancer CNE cells. Nat Prod Commun. 2010;5(1):151–6.
Huang SS, et al. Antioxidant and anti-inflammatory activities of aqueous extract of Centipeda minima. J Ethnopharmacol. 2013;147(2):395–405.
Qin Q, et al. Brevilin A inhibits NLRP3 inflammasome activation in vivo and in vitro by acting on the upstream of NLRP3-induced ASC oligomerization. Mole Immunol. 2021;135:116–26.
Xue PH, et al. Cytotoxic and anti-inflammatory sesquiterpenes from the whole plants of centipeda minima. J Nat Prod. 2021;84(2):247–58.
Chan BD, et al. Centipeda minima extract attenuates dextran sodium sulfate-induced acute colitis in mice by inhibiting macrophage activation and monocyte chemotaxis. Front Pharmacol. 2021;12:738139.
Liang HX, et al. Antibacterial thymol derivatives isolated from Centipeda minima. Molecules. 2007;12(8):1606–13.
Taylor RS, et al. Antibacterial constituents of the Nepalese medicinal herb, Centipeda minima. Phytochemistry. 1998;47(4):631–4.
Chen Y, et al. Anti-endometriosis mechanism of Jiawei Foshou San based on network pharmacology. Front Pharmacol. 2018;9:811.
He SB, et al. A computational toxicology approach to screen the hepatotoxic ingredients in traditional Chinese medicines: polygonum multiflorum thunb as a case study. Biomolecules. 2019;9(10):577.
He SS, et al. Uncovering the molecular mechanism of the Qiang-Xin 1 formula on sepsis-induced cardiac dysfunction based on systems pharmacology. Oxid Med Cell Longev. 2020;2020:3815185.
Kashyap D, et al. Fisetin and Quercetin: promising flavonoids with chemopreventive potential. Biomolecules. 2019;9(5):174.
Jeong JH, et al. Effects of low dose quercetin: cancer cell-specific inhibition of cell cycle progression. J Cell Biochem. 2009;106(1):73–82.
Ren MX, et al. Effect of quercetin on the proliferation of the human ovarian cancer cell line SKOV-3 in vitro. Exp Ther Med. 2015;10(2):579–83.
Granato M, et al. Quercetin induces apoptosis and autophagy in primary effusion lymphoma cells by inhibiting PI3K/AKT/mTOR and STAT3 signaling pathways. J Nutr Biochem. 2017;41:124–36.
Huang KY, et al. Growth suppression in lung cancer cells harboring EGFR-C797S mutation by quercetin. Biomolecules. 2021;11(9):1271.
Guo HJ, et al. Quercetin induces pro-apoptotic autophagy via SIRT1/AMPK signaling pathway in human lung cancer cell lines A549 and H1299 in vitro. Thoracic Cancer. 2021;12(9):1415–22.
Donehower LA, et al. Integrated analysis of TP53 gene and pathway alterations in the cancer genome atlas. Cell Rep. 2019;28(11):3010–3010.
Garancher A, et al. Tumor necrosis factor overcomes immune evasion in P53-mutant medulloblastoma. Neuro Oncol. 2019;21:93–93.
Kastenhuber ER, Lowe SW. Putting p53 in context. Cell. 2017;170(6):1062–78.
Cavic M, et al. TP53 and DNA-repair gene polymorphisms genotyping as a low-cost lung adenocarcinoma screening tool. J Clin Pathol. 2019;72(1):75–80.
Hu W, et al. TP53, TP53 Target Genes (DRAM, TIGAR), and Autophagy. Adv Exp Med Biol. 2019;1206:127–49.
Lacroix M, et al. Metabolic functions of the tumor suppressor p53: Implications in normal physiology, metabolic disorders, and cancer. Mole Metab. 2020;33:2–22.
Mijit M, et al. Role of p53 in the regulation of cellular senescence. Biomolecules. 2020;10(3):420.
Hu B, et al. High expression of CARM1 inhibits lung cancer progression by targeting TP53 by regulating CTNNB1. Lung. 2020;198(2):415–22.
Alwhaibi A, et al. Nodal pathway activation due to Akt1 suppression is a molecular switch for prostate cancer cell epithelial-to-mesenchymal transition and metastasis. Biochem Pharmacol. 2019;168:1–13.
Huo X, et al. Clinical and expression significance of AKT1 by co-expression network analysis in endometrial cancer. Front Oncol. 2019;9:1147.
Wang JJ, et al. MAT1 facilitates the lung metastasis of osteosarcoma through upregulation of AKT1 expression. Life Sci. 2019;234:116771.
Tian XY, et al. Hexokinase 2 promoted cell motility and proliferation by activating Akt1/p-Akt1 in human ovarian cancer cells. J Ovarian Res. 2022;15(1):92.
Rao GH, et al. Inhibition of AKT1 signaling promotes invasion and metastasis of non-small cell lung cancer cells with K-RAS or EGFR mutations. Sci Rep. 2017;7:7066.
Wu H, et al. miR-377-5p inhibits lung cancer cell proliferation, invasion, and cell cycle progression by targeting AKT1 signaling. J Cell Biochem. 2019;120(5):8120–8.
Duffy MJ, et al. MYC as a target for cancer treatment. Cancer Treat Rev. 2021;94:102154.
Chen SJ, et al. Roles of MYC-targeting long non-coding RNA MINCR in cell cycle regulation and apoptosis in non-small cell lung cancer. Respir Res. 2019;20(1):202.
Dho SH, et al. Development of a radionuclide-labeled monoclonal anti-CD55 antibody with theranostic potential in pleural metastatic lung cancer. Sci Rep. 2018;8:8960.
Wakelee H, et al. 50 Years of progress in the systemic therapy of non-small cell lung cancer. Am Soc Clin Oncol Educ Book. 2014;34:177–89.
Ji HH, et al. Prediction of the mechanisms by which quercetin enhances cisplatin action in cervical cancer: a network pharmacology study and experimental validation. Front Oncol. 2022;11:780387.
This work was supported by the Natural Science Foundation of Shandong Province (No. ZR202111290609) and the Science and Technology Project of Binzhou Medical University (No. BY2021KYQD33).
Ethics approval and consent to participate
All methods were carried out in accordance with relevant guidelines and regulations. The protocol established in this research was approved by the Institutional Animal Care and Use Committee of Binzhou Medical University Hospital, Approval Number: 20211008-20. The study was carried out in compliance with the ARRIVE guidelines.
Consent for publication
All authors declared that they have no conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Gao, C., Pan, H., Ma, F. et al. Centipeda minima active components and mechanisms in lung cancer. BMC Complement Med Ther 23, 89 (2023). https://doi.org/10.1186/s12906-023-03915-y
- Centipeda minima
- Lung cancer
- Network pharmacology
- Molecular docking
- Experimental validation
- Molecular mechanism