- 1Laboratory of Environmental Pollution and Integrative Omics, Guilin Medical University, Guilin, China
- 2Guilin Center for Disease Control and Prevention, Guilin, China
- 3Guangxi Key Laboratory of Tumor Immunology and Microenvironmental Regulation, Guilin Medical University, Guilin, China
- 4Department of Psychiatry, Icahn School of Medicine at Mount Sinai, New York, NY, United States
- 5State Key Laboratory of Agrobiotechnology, School of Life Sciences, The Chinese University of Hong Kong, Shatin, Hong Kong SAR, China
- 6Department of Pharmacy, Guilin Medical University, Guilin, China
The coronavirus disease 2019 (COVID-19) pandemic has led to 4,255,892 deaths worldwide. Although COVID-19 vaccines are available, mutant forms of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have reduced the effectiveness of vaccines. Patients with cancer are more vulnerable to COVID-19 than patients without cancer. Identification of new drugs to treat COVID-19 could reduce mortality rate, and traditional Chinese Medicine(TCM) has shown potential in COVID-19 treatment. In this study, we focused on lung adenocarcinoma (LUAD) patients with COVID-19. We aimed to investigate the use of curcumol, a TCM, to treat LUAD patients with COVID-19, using network pharmacology and systematic bioinformatics analysis. The results showed that LUAD and patients with COVID-19 share a cluster of common deregulated targets. The network pharmacology analysis identified seven core targets (namely, AURKA, CDK1, CCNB1, CCNB2, CCNE1, CCNE2, and TTK) of curcumol in patients with COVID-19 and LUAD. Clinicopathological analysis of these targets demonstrated that the expression of these targets is associated with poor patient survival rates. The bioinformatics analysis further highlighted the involvement of this target cluster in DNA damage response, chromosome stability, and pathogenesis of LUAD. More importantly, these targets influence cell-signaling associated with the Warburg effect, which supports SARS-CoV-2 replication and inflammatory response. Comparative transcriptomic analysis on in vitro LUAD cell further validated the effect of curcumol for treating LUAD through the control of cell cycle and DNA damage response. This study supports the earlier findings that curcumol is a potential treatment for patients with LUAD and COVID-19.
Key points
• COVID-19-/LUAD-associated and curcumol targets were identified.
• Prognostic value of curcumol against LUAD and COVID-19 was characterized.
• We identified seven core pharmacological targets of curcumol, namely, AURKA, CDK1, CCNB1, CCNB2, CCNE1, CCNE2, and TTK, in treating LUAD and COVID-19.
• Comparative transcriptomic analysis specified the effects of curcumol for treating LUAD through control of cell cycle, DNA damage response, and cell apoptosis.
Introduction
Coronavirus disease 2019 (COVID-19) is a novel disease characterized by high infectivity and rapid spread. Widespread community transmission of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has led to a global pandemic. As of 4 March 2022, 440,807,756 confirmed cases of COVID-19 and 5,978,096 deaths had been reported worldwide, according to the World Health Organization (WHO) (https://worldhealthorg.shinyapps.io/covid/). A considerable proportion of patients with COVID-19 related critical illness have comorbidities, which are associated with increased mortality.
The prospective cohort studies have demonstrated that patients with COVID-19 with underlying malignancies have a higher mortality rate than those without cancer (1, 2). The recent studies have shown that the levels of angiotensin-I-converting enzyme 2 (ACE2) and transmembrane serine protease 2 (TMPRSS2) in patients with lung adenocarcinoma (LUAD) are significantly increased, and the increased levels of these enzymes are associated with susceptibility of patients with LUAD to SARS-CoV-2 (3), because ACE2 serves as an important binding site for SARS-CoV-2, leading to facilitate viral entry into target host cells (4). Lung adenocarcinoma is the most common subtype of lung cancer. The development of LUAD is stepwise, beginning with atypical adenomatous hyperplasia (AAH), and progressing to adenocarcinoma in situ (AIS) and then to minimally invasive adenocarcinoma (MIA) (5, 6). Most of the patients are diagnosed with advanced disease and have poor prognosis. Gene mutations (such as EGFR, KRAS, and BRAF mutations) and tumor inflammatory microenvironment are strongly associated with LUAD pathogenesis (7, 8).
A meta-analysis of patients with lung cancer and COVID-19, which included 13 studies, showed that the pooled mortality of patients with lung cancer and COVID-19 (up to 25–42%) was significantly higher than the mortality of patients with other cancers (9, 10). This may be due to different pathophysiological factors, such as pulmonary compromise and smoking history, in patients with lung cancer, compared with other cancers (11). Therefore, there is an urgent need to identify drugs to treat patients with lung cancer and COVID-19. In addition to helping these patients, identification of such drugs will also relieve the pressure on respiratory healthcare services.
The current research findings have shown that traditional Chinese medicine (TCM), such as Lianhua Qingwen Keli, Honeysuckle Flower Cold-Relieving Granules, and Xuebijing injection can be effective in preventing COVID-19 and relieving clinical symptoms of COVID-19. These TCM were officially recommended by the National Medical Products Administration, as adjunctive therapy, in the treatment of COVID-19 (12). TCM can potentially be used for the treatment of comorbidities associated with COVID-19, and has received worldwide attention. Curcumol, a sesquiterpenoid isolated from Curcumae rhizoma, has been shown to have various therapeutic effects, including anticancer, antioxidant, antimicrobial, and anti-inflammatory effects (13). An in vitro study by Li et al. showed that curcumol suppresses proliferation of the LUAD cells, A549 and H460, by arresting the cell cycle, altering the expression of apoptosis signaling pathways and inducing tumor cell apoptosis (14). In chronic asthmatic mice, curcumol was found to reduce pulmonary inflammation and airway remodeling by decreasing cytokine levels (15). In addition, curcumol was reported to inhibit LUAD growth and metastasis and overcome tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) resistance in lung cancer (16). The results of the aforementioned studies show that curcumol may be a potential treatment for patients with LUAD and COVID-19.
In the current study, we used network pharmacology, comparative transcriptome, and systematic bioinformatics analysis to investigate the use of curcumol for COVID-19 and LUAD treatment. We aimed to identify possible therapeutic targets and to unfold the molecular mechanisms underlying the therapeutic effects of curcumol in COVID-19 and LUAD, using clinicopathological analysis, gene ontology, KEGG enrichment analysis, ingenuity pathway analysis (IPA) and molecular docking.
Materials and Methods
Identification of Common Deregulated Targets Between COVID-19 and LUAD
For the identification of COVID-19-associated targets, the keywords “coronavirus COVID-19,” “coronavirus Disease 2019,” “severe acute respiratory syndrome coronavirus 2,” and “COVID-19” were subjected to different databases, including the Genecards database (17), Online Mendelian Inheritance in Man (OMIM) database (https://omim.org/), Therapeutic Target Database (TTD) (18), Comparative Toxicogenomics Database (CTD) (19), and National Center for Biotechnology Information (NCBI) (https://www.ncbi.nlm.nih.gov/), genes with the relevance score >1 were obtained from the databases. To identify LUAD-associated targets, the transcriptome data of patients with LUAD were obtained from The Cancer Genome Atlas (TCGA) database (https://portal.gdc.cancer.gov/) on 26 July 2021. Using the limma package in R on Bioconductor software, genes with FDR < 0.05, and |logfold change| > 2 were considered as differentially expressed genes (DEGs) (20).
Identification of Pharmacological Targets of Curcumol in LUAD and COVID-19 Treatment
The pharmacological targets of curcumol were determined using various online tools and databases, including Swiss Target Prediction database and Bioinformatics Analysis Tool for Molecular mechANism of TCM (BATMAN-TCM) (21). The target genes were subjected to UniProt for human database correction. The common deregulated genes between COVID-19 and LUAD, that were previously identified, were intersected and compared with the targets of curcumol. Interactions between common targets were analyzed using the STRING database (version 11.0) (22) and Cytoscape software (version 3.6.1) (23). The Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.8 was used for the gene ontology (GO) function enrichment analysis and KEGG enrichment analysis, to understand the functional roles of targets and signaling pathways controlled by the genes.
Binding of Curcumol to Predicted Targets
Molecular docking analysis was used to investigate the possible binding between curcumol and its predicted targets. The protein structures of the core targets were searched for from the Protein Data Bank (PDB) database (24). Identified protein structures were then docked with curcumol using the AutoDock Vina program and docking analysis was conducted (25).
The Roles of Curcumol Target Genes in LUAD Pathogenesis
To determine the pathological roles of the core targets in LUAD, Cox proportional hazards models were applied in univariate analysis of survival as a function of clinical variables and gene expression.
Cell Culture
The human lung adenocarcinomic cell lines of A549 were incubated with high glucose dulbecco's modified eagle medium (DMEM) medium (ThermoFisher, Cat. No. 11965118), supplemented with 0.5% penicillin–streptomycin (ThermoFisher, Cat. No. 15140122) and 5% fetal bovine serum (ThermoFisher, Cat. No. 10082147) under 5% CO2 at 37°C.
Cell Proliferation Assays
The cells were seeded in a 96-well plate at a cell density of 2 × 104 cells per well, with eight replicate wells. The cells were treated with different concentrations of curcumol (0.1–100 μM) for 48 h. After the incubation, the cell viability was measured by the CCK-8 assay (Data Inventory Biotechnology) as described previously (26). The colorimetric product formed was measured at an absorbance of 450 nm and 600 nm, ΔOD = OD450nm – OD600nm.
The RNA Sequencing
After the treatment of the cell with 100 μM curcumol for 48 h, the total RNA of the cell was extracted using Trizol reagent (Thermofisher) following the manufacturer's instruction. The RNA quality and quantity were assessed by using Bioanalyzer 2,100 and RNA 6,000 Nano LabChip Kit (Agilent), high-quality RNA samples with RNA integrity number (RIN) number higher than 7.0 were used to construct sequencing library. The average insert size for the final complementary DNA (Cdna) library was about 300 bp. Then 2 × 150 bp paired-end sequencing (PE150) was performed on an Illumina Novaseq™. The high-quality clean reads were mapped to the Human genome reference (Homo sapiens Ensembl v96) using HISAT2 software (version: hisat2-2.0.4) (27). StringTie and ballgown were used to determine the gene expression level (28). The genes with a 1.5 < fold change (treatment/control) < 0.75 and –log10 (q-value) > 1.3 were considered as DEGs. The DEGs were subjected to the DAVID v6.8 analysis (29) and IPA (https://www.qiagenbioinformatics.com/products/ingenuity-pathway-analysis) to delineate the molecular mechanism underlying the effect of curcumol for treating LUAD.
Results
Identification of Pharmacological Targets of Curcumol in COVID-19 and LUAD Treatment
Using the relevant databases, we identified a total of 8,339 targets associated with COVID-19 (Figure 1A). Using the TCGA database, we found 5,538 differential expressed genes associated with LUAD (Figure 1A). When we compared the COVID-19- and LUAD-associated targets, we found 882 shared targets (Figure 1A), among which 216 were downregulated and 666 were upregulated, in patients with LUAD (Figure 1B). To understand the pharmacology of curcumol, a network pharmacology analysis was conducted. We identified 151 curcumol-associated targets using the mentioned databases (Figure 1A) and, after comparison of the curcumol-associated targets with the COVID-19/ LUAD-associated targets, we found 28 targets shared by curcumol, COVID-19, and LUAD (Figure 1A). The molecular network analysis using Cytoscape highlighted seven core targets of curcumol, namely, AURKA, CDK1, CCNB1, CCNB2, CCNE1, CCNE2, and TTK in COVID-19 and LUAD (Figure 1C and Table 1).
Figure 1. Identification of curcumol-targeted COVID-19- LUAD-associated targets. (A) Venn diagram showing the intersecting targets of curcumol /COVID-19/LUAD. (B) Volcano plot showing the expression level of differential expressed COVID-19- LUAD-associated genes in patients with LUAD. The genes with |log2 (fold change)| > 1 and –log10(FDR) > 1.3 were considered as differential expressed genes. (C) Protein–protein interaction analysis of curcumol/COVID-19/LUAD-intersecting genes using STRING tool.
Binding of Curcumol to Targets CDK1, TTK, and AURKA
Three protein structures, CDK1, TTK, and AURKA, out of the seven core targets were available on the PDB database (AURKA, ID:2J50; CDK1, ID:5HQ0, and TTK, ID:6N6O). These protein structures were subjected to docking analysis with curcumol using the AutoDock Vina program. The results were displayed using PyMOL (version 2.3), which showed that curcumol formed a hydrogen bond with LYS-162 (3.2 Å) of AURKA (PDB ID: 2J50) (Figure 2A), and the binding affinity of curcumol for AURKA was −6.4 kcal/mol. A similar bindings were observed between curcumol and the amino acid residue LEU-83 (2.2 Å) of CDK1 (PDB ID: 5HQ0) (Figure 2B) and between curcumol and the amino acid residue LYS-529 (3.0 Å) of TTK (PDB ID: ID:6N6O) (Figure 2C). The binding affinities of curcumol for CDK1 and TTK were −3.2 kcal/mol and −4.7 kcal/mol, respectively. Our data suggested the direct binding of curcumol to its targets.
Figure 2. The direct binding of curcumol to CDK1, TTK, and AURKA. The protein structures of CDK1, TTK, and AURKA, obtained from the PDB database. Using the programs Autodock and Visual Molecular Dynamics for visualization, hydrogen bonds can be seen between curcumol and (A) the amino acid residue LYS-162 (3.2 Å) of AURKA (PDB ID: 2J50), (B) the amino acid residue LEU-83 (2.2 Å) of CDK1 (PDB ID: 5HQ0), and (C) the amino acid residue LYS-529 (3.0 Å) of TTK (PDB ID: ID:6N6O).
The Roles of Curcumol Target Genes in LUAD Pathogenesis
The results of the hazards models showed that the expression of curcumol targets in COVID-19 and LUAD was significantly associated with the relative risk of survival [AURKA (p = 0.001, hazard ratio, 1.086–1.399); CDK1 (p < 0.001, hazard ratio, 1.099–1.401); CCNB1 (p < 0.001, hazard ratio, 1.144–1.499); CCNB2 (p < 0.001, hazard ratio, 1.089–1.388); CCNE1 (p = 0.012, hazard ratio, 1.032–1.283); CCNE2 (p = 0.026, hazard ratio, 1.019–1.341); TTK (p = 0.004, hazard ratio, 1.053–1.307)] in patients with LUAD (Figure 3A). The results of survival analysis using the Kaplan–Meier estimator also showed that patients with LUAD with greater expression of the genes, AURKA, CDK1, CCNB1, CCNB2, CCNE2, and TTK, had poorer overall survival rates (Figure 3B). In addition, the correlation analysis highlighted that increased expression of AURKA was associated with advanced stages of LUAD and the increased number of lymph nodes containing tumor (Figure 3C).
Figure 3. Clinicopathologic analysis of curcumol-targeted COVID-19/LUAD-associated targets. (A) The Univariate Cox proportional hazards models showing that the expression of curcumol target genes in COVID-19 and LUAD was significantly associated with relative risk of survival in the patients with LUAD. (B) The survival analysis using Kaplan–Meier estimator showing that the patients with LUAD with higher expression of AURKA, CDK1, CCNB1, CCNB2, CCNE2, and TTK had poorer overall survival rates. (C) The higher expressions of AURKA in patients with LUAD are associated with advanced stages of LUAD. T, Staging of tumor; N, Number of lymph nodes containing tumor.
Curcumol Target Genes Mediated DNA Damage Response and Cell Cycle Control
To further understand the biological roles of curcumol in the treatment of COVID-19 and LUAD, the identified core targets were subjected to GO enrichment analysis and KEGG pathway enrichment analysis. The results of GO analysis showed that curcumol targets (AURKA, CDK1, and CCNB1) could mediate biological processes related to DNA damage response by controlling DNA integrity and the DNA damage checkpoint (Figure 4A and Supplementary Table 1), or by regulating cyclin-dependent protein kinase activity, which could lead to cell cycle arrest at different cell cycle checkpoints, such as cell cycle G1/S phase transition and cell cycle G2/M phase transition (Figure 4B and Supplementary Table 1). These responses are controlled by the seven curcumol targets (Figure 4B and Supplementary Table 1). The other possible outcome of mediation of biological processes related to the DNA damage response by these targets is the alteration of cellular assembly processes, such as organelle fission, mitotic nuclear envelope disassembly, chromosome segregation, and mitotic spindle organization (Figure 4C and Supplementary Table 1). In addition, curcumol targets were found to regulate many important processes involved in LUAD carcinogenesis, such as histone phosphorylation, oxidative phosphorylation, and cellular respiration (Figure 4D and Supplementary Table 1).
Figure 4. Functional characterization of curcumol-targeted COVID-19/ LUAD-associated targets. (A) The GO enrichment analysis highlighted the biological processes related to DNA damage response that are controlled by curcumol-targeted COVID-19/LUAD-associated targets. The circos plot showing the involvement of AUKRA, CDK1, and CCNB1 in the enriched biological processes. (B) The GO enrichment analysis highlighted the biological processes related to cell cycle control regulated by curcumol-targeted COVID-19/LUAD-associated targets. Circos plot showing the involvement of AUKRA, CDK1, CCNB1, CCNB2, CCNE1, CCNE2, and TTK in the enriched biological processes. (C) The GO enrichment analysis highlighted the biological processes related to cellular assembly processes controlled by curcumol targeted COVID-19/LUAD-associated targets. Circos plot showing the involvement of AUKRA, CDK1, CCNB1, CCNB2, CCNE1, CCNE2, and TTK in the enriched biological processes. (D) The GO enrichment analysis highlighted the biological processes related to carcinogenesis of LUAD. Circos plot showing the involvement of AUKRA, CDK1, CCNB1, and TTK in the enriched biological processes. The size of the dot represents the number of targets. The color intensity of the dot represents the significance of the processes.
The GO cellular components and molecular functions analyses showed that the curcumol targets play roles in many enzymatic complexes related to cell cycle control, such as in cyclin-dependent protein kinase holoenzyme complex and in serine/threonine protein kinase complex (Figure 5A and Supplementary Table 2). More importantly, the targets are involved in chromosome organization, specifically, in organization of centromeric region of chromosomes, kinetochores, mitotic spindle pores, and telomeric regions (Figure 5A and Supplementary Table 2). The results of the molecular function analysis further highlighted the effects of curcumol on cyclin-dependent protein serine/threonine kinase regulator activity, histone kinase activity, and cyclin binding (Figure 5A and Supplementary Table 2). KEGG pathway analysis was used to investigate the pharmacological effect of curcumol on the regulation of cell-signaling pathways. The results showed that curcumol could target many cell-signaling pathways related to cancer development (Figure 5B and Supplementary Table 3), including the p53 signaling pathway, FoxO signaling pathway, and PI3K–Akt signaling pathway (Figure 5B and Supplementary Table 3). The curcumol targets were also found to play a role in immune response to viral infections, like T-cell leukemia virus 1 infection, papillomavirus infection, and immunodeficiency virus 1 infection (Figure 5B and Supplementary Table 3).
Figure 5. Molecular functions and signaling pathways controlled by curcumol-targeted COVID-19/LUAD-associated targets. (A) The GO enrichment analysis highlighted the involved enzymatic complexes, chromosome compartments, and molecular functions in curcumol-targeted COVID-19/LUAD-associated targets. The size of the dot represents the number of targets. The color intensity of the dot represents the significance of the terms. The Circos plot shows the involvement of AUKRA, CDK1, CCNB1, CCNB2, CCNE1, CCNE2, and TTK in the enriched terms. (B) The KEGG analysis highlighted the cell-signaling pathways related to cancer development, mediated by curcumol, against COVID-19 and LUAD. The size of dot represents the number of targets. The color intensity of the dot represents the significance of the pathways. The Circos plot shows the involvement of AUKRA, CDK1, CCNB1, CCNB2, CCNE1, CCNE2, and TTK in the enriched signaling.
Curcumol Inhibited the Cell Proliferation of LUAD Through the Control of Cell Cycle and DNA Damage
To investigate the effect of curcumol on lung cancer, an in vitro LUAD model A549 was used. Our result showed that the treatment of curcumol caused a significant dose-dependent inhibition of cell proliferation in LUAD as compared to the control group (Figure 6A). Then comparative transcriptomic analysis was conducted to delineate the corresponding molecular mechanism. By comparing the gene expression profile of control and curcunol treatment group, we found 348 DEGs, including 206 upregulated genes and 142 downregulated genes (Figure 6B and Supplementary Table 4). The DEGs were used for the DAVID and IPA analysis to understand the biological alteration and gene network mediated by curcumol treatment in LUAD. The result of GO enrichment showed that the curcumol treatment controlled biological processes related to cell cycle and gene transcription of LUAD (Figure 6C and Supplementary Table 5). More importantly, curcumol could trigger DNA damage response, leading to cell death and cell apoptosis in LUAD (Figure 6C and Supplementary Table 5). These results further supported the above findings from network pharmacology. Finally, gene networking of IPA highlighted the involvement of transcription and translation factors, enzymes, kinases, phosphatases, and receptors in curcumol-mediated cell cycle and DNA response (Figure 6D and Table 2).
Figure 6. Curcumol-targeted genes involved in cell cycle and DNA damage of LUAD cell model A549. (A) The curcumol treatment (100 μM) inhibited cell proliferation of LUAD cell in a dose-dependent manner; N = 8 and three experiments. (B) The comparative transcriptomic analysis showed the differential gene expression in LUAD cell caused by curcumol treatment (100 μM). Genes with a 1.5 < fold change (treatment/control) < 0.75 and –log10 (q-value) > 1.3 were considered as DEGs. Blue dots represented downregulated genes and red dots represented upregulated genes. (C) The GO enrichment analysis highlighted the importance of curcumol-dysregulated genes in cell cycle, DNA damage, and cell apoptosis of LUAD cell. The size of bubble represented the number of DEGs involved in the processes, the color of bubble represented the significance of the processes. (D) Gene networking of IPA highlighted the involvement of DEGs in cell cycle and DNA damage response. Red color represented the upregulated genes, green color represented downregulated genes. Orange represented the predicted activated molecules and blue color represented the predicted inhibited molecules in the signaling. * represented p < 0.05, as compared to control group.
Discussion
Using the network pharmacology analysis, we identified seven core targets of curcumol against COVID-19 and LUAD. These seven targets, namely, CCNB1, CCNB2, CCNE1, CCNE2, AURKA, CDK1, and TTK, have been reported to play critical roles in the carcinogenesis and development of LUAD. These targets mediate signaling pathways, such as PI3K/AKT and p53, that are associated with the Warburg effect, which supports SARS-CoV-2 replication and inflammatory response (30, 31). Four cyclin family members (CCNB1, CCNB2, CCNE1, and CCNE2) were identified as core targets in our analysis. Also, cyclin B1 (CCNB1) is associated with poor prognosis in LUAD (32). It is a downstream effector of monoacylglycerol lipase (MGLL), a key enzyme in lipid metabolism which plays an oncogenic role in LUAD progression and metastasis (33). The cyclin family member CCNB1 also contributes to lung inflammation and oxidative stress (34, 35); CCNB2 is an independent predictor of the prognosis of patients with LUAD (36). The functional characterization further highlighted the involvement of CCNB2 in inducing cell cycle arrest and apoptosis in LUAD cells (37). Moreover, the high levels of CCNB2 activate inflammation-induced motility in LUAD (38). Additionally, a study by Ma et al. revealed that knockdown of CCNB2 suppressed proliferation of the LUAD cell line (39). In addition to cyclin B members, our results show that cyclin E members, such as CCNE1 and CCNE2, are also targeted by curcumol; CCNE1 plays a role in progression, cell proliferation, and cell cycle arrest of lung cancer cells (40–42).
Using molecular docking, we found that curcumol binds directly to a group of kinases, including CDK1, TTK, and Aurora kinase A (AURKA). AURKA, a cell cycle kinase, is associated with many cancer types (43). An in vitro study of human LUAD cell lines demonstrated that AURKA plays an important role in the proliferation of LUAD cells, through the regulation of multiple downstream effectors, such as RAF-1, CCND2, CCND3, CDK4, PAK4, and EGFR (44). In addition to its function as a cell cycle kinase, AURKA is considered as an inhibitor which prevents the chromatin assembly of functional replisomes, leading to sensitization of cancer cells to combination therapy (45). Further functional characterization studies of AURKA report that AURKA suppression enhances the radiosensitivity of lung cancer and its response to EGFR inhibitors (46, 47). In addition, the downregulating AURKA inhibits docetaxel chemoresistance in LUAD (48), suggesting that AURKA is a promising target for LUAD therapy. Many studies have also demonstrated the autoimmune and inflammatory roles of AURKA via regulation of M1 macrophage polarization (49, 50).
Our results showed that the curcumol targets the other cell cycle gatekeeper kinases, such as cyclin-dependent kinase 1 (CDK1); CDK1 is a potential prognostic biomarker of/and target for lung cancer (51); CDK1 activity is critical for JAK/STAT3 signaling activation, and the inhibition of CDK1 can suppress lung cancer (52). In addition, an in vitro study of LUAD cells showed that reduced CDK1 activity led to cell cycle arrest and promotion of apoptosis in LUAD (53, 54). CDK1 controls many effectors involved in cell cycle regulation, such as FOXM1, TRAP1, and GCN1 (55–57). Furthermore, CDK1 plays a role in the DNA damage response. For example, CDK1 ensures optimal Fun30 phosphorylation and checkpoint activation at DNA double-strand breaks and plays an important role in the DNA damage response by preventing the formation of lagging chromosomes (58, 59). Also, CDK1 ensures accurate chromosomal segregation via the activity of acetyltransferase TIP60 and chromatin remodeller RSF1 (60, 61). Cumulative studies have reported that DNA damage is associated with higher mortality in patients with COVID-19 (62). Kinase threonine tyrosine kinase (TTK) is a critical component of the spindle assembly checkpoint (63). Also, TTK is a biomarker for prognosis of Non-small cell lung cancer (64), and the upregulation of TTK increases the cancer progression in lung cancer (65). In addition, the TTK antagonism has marked antineoplastic effects against LUAD (66). The targeting of CDK1 and TTK by curcomol provides an opportunity to treat patients with LUAD and COVID-19. The results were further validated by using comparative transcriptomic analysis on LUAD cell. The treatment of curcumol could inhibited cell proliferation of LUAD through its control on cell cycle and DNA damage response.
Conclusions
In conclusion, we identified the pharmacological targets and the therapeutic mechanisms of curcumol in the treatment of COVID-19 and LUAD, including immune response, DNA damage response, and cell cycle arrest, and regulation of cell-signaling pathways such as the p53 signaling pathway, FoxO signaling pathway, and PI3K-Akt signaling pathway. The results were further supported by the comparative transcriptomic analysis on in vitro LUAD cell, suggesting that curcumol has potential for treating patients with LUAD and COVID-19. However, further Pre-clinical study is needed to warrant the findings of the present study before the clinical use.
Data Availability Statement
Sequencing data of transcriptome sequencing that support the findings of this study have been deposited in the NCBI BioProject database (https://www.ncbi.nlm.nih.gov/bioproject) with the BioProject accession code PRJNA793079.
Author Contributions
RL, KL, and XC contributed to the conception, design of the manuscript, drafted this manuscript, and revised this manuscript. LY, HX, XLi, YL, HZ, XLin, and TC contributed to the acquisition, analysis, and interpretation of data in this manuscript. All authors agree to be accountable for all aspects of work ensuring integrity and accuracy. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Natural Science Foundation of China (Nos. 82160282, 82160768, and 82002822) and Natural Science Foundation of Guangxi Province (No. GuikeZD20302006).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnut.2022.870370/full#supplementary-material
Abbreviations
LUAD, lung adenocarcinoma; COVID-19, coronavirus disease 2019; TCGA, The Cancer Genome Atlas; GO, Gene Ontology; BP, Biological process; KEGG, Kyoto Encyclopedia of Genes and Genomes; IPA, Ingenuity Pathway Analysis.
References
1. Assaad S, Avrillon V, Fournier ML, Mastroianni B, Russias B, Swalduz A, et al. High mortality rate in cancer patients with symptoms of COVID-19 with or without detectable SARS-COV-2 on RT-PCR. Eur J Cancer. (2020) 135:251–9. doi: 10.1016/j.ejca.2020.05.028
2. Lee LYW, Cazier JB, Starkey T, Briggs SEW, Arnold R, Bisht V, et al. COVID-19 prevalence and mortality in patients with cancer and the effect of primary tumour subtype and patient demographics: a prospective cohort study. Lancet Oncol. (2020) 21:1309–16. doi: 10.1016/S1470-2045(20)30442-3
3. Kong Q, Xiang Z, Wu Y, Gu Y, Guo J, Geng F. Analysis of the susceptibility of lung cancer patients to SARS-CoV-2 infection. Mol Cancer. (2020) 19:80. doi: 10.1186/s12943-020-01209-2
4. Dhaundiyal A, Kumari P, Jawalekar SS, Chauhan G, Kalra S, Navik U. Is highly expressed ACE 2 in pregnant women “a curse” in times of COVID-19 pandemic? Life Sci. (2021) 264:118676. doi: 10.1016/j.lfs.2020.118676
5. Cha YJ, Shim HS. Biology of invasive mucinous adenocarcinoma of the lung. Transl Lung Cancer Res. (2017) 6:508–12. doi: 10.21037/tlcr.2017.06.10
6. Butnor KJ. Controversies and challenges in the histologic subtyping of lung adenocarcinoma. Transl Lung Cancer Res. (2020) 9:839–46. doi: 10.21037/tlcr.2019.12.30
7. Lou Y, Diao L, Cuentas ER, Denning WL, Chen L, Fan YH, et al. Epithelial-mesenchymal transition is associated with a distinct tumor microenvironment including elevation of inflammatory signals and multiple immune checkpoints in lung adenocarcinoma. Clin Cancer Res. (2016) 22:3630–42. doi: 10.1158/1078-0432.CCR-15-1434
8. Inamura K. Clinicopathological characteristics and mutations driving development of early lung adenocarcinoma: tumor initiation and progression. Int J Mol Sci. (2018) 19:1259. doi: 10.3390/ijms19041259
9. Luo J, Rizvi H, Preeshagul IR, Egger JV, Hoyos D, Bandlamudi C, et al. COVID-19 in patients with lung cancer. Ann Oncol. (2020) 31:1386–96. doi: 10.1016/j.annonc.2020.06.007
10. Lei H, Yang Y, Zhou W, Zhang M, Shen Y, Tao D, et al. Higher mortality in lung cancer patients with COVID-19? A systematic review and meta-analysis. Lung Cancer. (2021) 157:60–5. doi: 10.1016/j.lungcan.2021.05.002
11. Passaro A, Bestvina C, Velez Velez M, Garassino MC, Garon E, Peters S. Severity of COVID-19 in patients with lung cancer: evidence and challenges. J Immunother Cancer. (2021) 9:e002266. doi: 10.1136/jitc-2020-002266
12. Ni L, Chen L, Huang X, Han C, Xu J, Zhang H, et al. Combating COVID-19 with integrated traditional Chinese and Western medicine in China. Acta Pharm Sin B. (2020) 10:1149–62. doi: 10.1016/j.apsb.2020.06.009
13. Wei W, Rasul A, Sadiqa A, Sarfraz I, Hussain G, Nageen B, et al. Curcumol: from plant roots to cancer roots. Int J Biol Sci. (2019) 15:1600–9. doi: 10.7150/ijbs.34716
14. Li S, Zhou G, Liu W, Ye J, Yuan F, Zhang Z. Curcumol inhibits lung adenocarcinoma growth and metastasis via inactivation of PI3K/AKT and Wnt/ß-catenin pathway. Oncol Res. (2020) 28:685–700. doi: 10.3727/096504020X15917007265498
15. Jia S, Guo P, Lu J, Huang X, Deng L, Jin Y, et al. Curcumol ameliorates lung inflammation and airway remodeling via inhibiting the abnormal activation of the Wnt/β-catenin pathway in chronic asthmatic mice. Drug Des Devel Ther. (2021) 15:2641–51. doi: 10.2147/DDDT.S292642
16. Zhang J, Zhou Y, Li N, Liu WT, Liang JZ, Sun Y, et al. Curcumol overcomes TRAIL resistance of non-small cell lung cancer by targeting NRH:Quinone Oxidoreductase 2 (NQO2). Adv Sci. (2020) 7:2002306. doi: 10.1002/advs.202002306
17. Stelzer G, Rosen N, Plaschkes I, Zimmerman S, Twik M, Fishilevich S, et al. The GeneCards suite: from gene data mining to disease genome sequence analyses. Curr Protoc Bioinformatics. (2016) 54:1.30.1–1.30.33. doi: 10.1002/cpbi.5
18. Wang Y, Zhang S, Li F, Zhou Y, Zhang Y, Wang Z, et al. Therapeutic target database 2020: enriched resource for facilitating research and early development of targeted therapeutics. Nucleic Acids Res. (2020) 48:1031–41. doi: 10.1093/nar/gkz981
19. Davis AP, Grondin CJ, Johnson RJ, Sciaky D, McMorran R, Wiegers J, et al. The comparative toxicogenomics database: update 2019. Nucleic Acids Res. (2019) 47:948–54. doi: 10.1093/nar/gky868
20. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. (2015) 43:e47. doi: 10.1093/nar/gkv007
21. Liu Z, Guo F, Wang Y, Li C, Zhang X, Li H, et al. BATMAN-TCM: a bioinformatics analysis tool for molecular mechANism of traditional Chinese medicine. Sci Rep. (2016) 6:21146. doi: 10.1038/srep21146
22. Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, et al. STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res. (2015) 43:447–52. doi: 10.1093/nar/gku1003
23. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. (2003) 13:2498–504. doi: 10.1101/gr.1239303
24. Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, et al. The protein data bank. Nucleic Acids Res. (2000) 28:235–42. doi: 10.1093/nar/28.1.235
25. Trott O, Olson AJ. AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J Comput Chem. (2010) 31:455–61. doi: 10.1002/jcc.21334
26. Wu K, Guo C, Yang B, Wu X, Wang W. Antihepatotoxic benefits of Poria cocos polysaccharides on acetaminophen-lesioned livers in vivo and in vitro. J Cell Biochem. (2019) 120:7482–8. doi: 10.1002/jcb.28022
27. Kim D, Paggi JM, Park C, Bennett C, Salzberg SL. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat Biotechnol. (2019) 37:907–15. doi: 10.1038/s41587-019-0201-4
28. Pertea M, Kim D, Pertea GM, Leek JT, Salzberg SL. Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat Protoc. (2016) 11:1650–67. doi: 10.1038/nprot.2016.095
29. Huang da W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. (2009) 4:44–57. doi: 10.1038/nprot.2008.211
30. Icard P, Lincet H, Wu Z, Coquerel A, Forgez P, Alifano M, et al. The key role of Warburg effect in SARS-CoV-2 replication and associated inflammatory response. Biochimie. (2021) 180:169–77. doi: 10.1016/j.biochi.2020.11.010
31. Lodhi N, Singh R, Rajput SP, Saquib Q. SARS-CoV-2: understanding the transcriptional regulation of ACE2 and TMPRSS2 and the role of single nucleotide polymorphism (SNP) at Codon 72 of p53 in the innate immune response against virus infection. Int J Mol Sci. (2021) 22:8660. doi: 10.3390/ijms22168660
32. Lou X, Ning J, Liu W, Li K, Qian B, Xu D, et al. YTHDF1 promotes cyclin B1 translation through m6A modulation and contributes to the poor prognosis of lung adenocarcinoma with KRAS/TP53 co-mutation. Cells. (2021) 10:1669. doi: 10.3390/cells10071669
33. Zhang L, He M, Zhu W, Lv X, Zhao Y, Yan Y, et al. Identification of a panel of mitotic spindle-related genes as a signature predicting survival in lung adenocarcinoma. J Cell Physiol. (2020) 235:4361–75. doi: 10.1002/jcp.29312
34. Sun XQ, Wu C, Qiu YB, Wu YX, Chen JL, Huang JF, et al. Heme oxygenase-1 attenuates seawater drowning-induced acute lung injury through a reduction in inflammation and oxidative stress. Int Immunopharmacol. (2019) 74:105634. doi: 10.1016/j.intimp.2019.05.019
35. Fountain MD, McLellan LA, Smith NL, Loughery BF, Rakowski JT, Tse HY, et al. Isoflavone-mediated radioprotection involves regulation of early endothelial cell death and inflammatory signaling in Radiation-Induced lung injury. Int J Radiat Biol. (2020) 96:245–56. doi: 10.1080/09553002.2020.1683642
36. Zhang H, Guo W, Zhang F, Li R, Zhou Y, Shao F, et al. Monoacylglycerol lipase knockdown inhibits cell proliferation and metastasis in lung adenocarcinoma. Front Oncol. (2020) 10:559568. doi: 10.3389/fonc.2020.559568
37. Cheng Y, Hou K, Wang Y, Chen Y, Zheng X, Qi J, et al. Identification of prognostic signature and gliclazide as candidate drugs in lung adenocarcinoma. Front Oncol. (2021) 11:665276. doi: 10.3389/fonc.2021.665276
38. Zhou H, Wang L, Huang J, Jiang M, Zhang X, Zhang L, et al. High EGFR_1 inside-out activated inflammation-induced motility through SLC2A1-CCNB2-HMMR-KIF11-NUSAP1-PRC1-UBE2C. J Cancer. (2015) 6:519–24. doi: 10.7150/jca.11404
39. Ma X, Yang S, Jiang H, Wang Y, Xiang Z. Transcriptomic analysis of tumor tissues and organoids reveals the crucial genes regulating the proliferation of lung adenocarcinoma. J Transl Med. (2021) 19:368. doi: 10.1186/s12967-021-03043-6
40. She Y, Han Y, Zhou G, Jia F, Yang T, Shen Z. hsa_circ_0062389 promotes the progression of non-small cell lung cancer by sponging miR-103a-3p to mediate CCNE1 expression. Cancer Genet. (2020) 241:12–9. doi: 10.1016/j.cancergen.2019.12.004
41. Zheng C, Li X, Ren Y, Yin Z, Zhou B. Long noncoding RNA RAET1K enhances CCNE1 expression and cell cycle arrest of lung adenocarcinoma cell by sponging miRNA-135a-5p. Front Genet. (2020) 10:1348. doi: 10.3389/fgene.2019.01348
42. Yao Y, Luo J, Sun Q, Xu T, Sun S, Chen M, et al. HOXC13 promotes proliferation of lung adenocarcinoma via modulation of CCND1 and CCNE1. Am J Cancer Res. (2017) 7:1820–34.
43. Melichar B, Adenis A, Lockhart AC, Bennouna J, Dees EC, Kayaleh O, et al. Safety and activity of alisertib, an investigational aurora kinase A inhibitor, in patients with breast cancer, small-cell lung cancer, non-small-cell lung cancer, head and neck squamous-cell carcinoma, and gastro-oesophageal adenocarcinoma: a five-arm phase 2 study. Lancet Oncol. (2015) 16:395–405. doi: 10.1016/S1470-2045(15)70051-3
44. Zhong N, Shi S, Wang H, Wu G, Wang Y, Ma Q, et al. Silencing Aurora-A with siRNA inhibits cell proliferation in human lung adenocarcinoma cells. Int J Oncol. (2016) 49:1028–38. doi: 10.3892/ijo.2016.3605
45. Guarino Almeida E, Renaudin X, Venkitaraman AR. A kinase-independent function for AURORA-A in replisome assembly during DNA replication initiation. Nucleic Acids Res. (2020) 48:7844–55. doi: 10.1093/nar/gkaa570
46. Liu N, Wang YA, Sun Y, Ecsedy J, Sun J, Li X, et al. Inhibition of Aurora A enhances radiosensitivity in selected lung cancer cell lines. Respir Res. (2019) 20:230. doi: 10.1186/s12931-019-1194-8
47. Shah KN, Bhatt R, Rotow J, Rohrberg J, Olivas V, Wang VE, et al. Aurora kinase A drives the evolution of resistance to third-generation EGFR inhibitors in lung cancer. Nat Med. (2019) 25:111–8. doi: 10.1038/s41591-018-0264-7
48. Cao J, Geng J, Chu X, Wang R, Huang G, Chen L. miRNA-885-3p inhibits docetaxel chemoresistance in lung adenocarcinoma by downregulating Aurora A. Oncol Rep. (2019) 41:1218–30. doi: 10.3892/or.2018.6858
49. Ding L, Gu H, Gao X, Xiong S, Zheng B. Aurora kinase a regulates m1 macrophage polarization and plays a role in experimental autoimmune encephalomyelitis. Inflammation. (2015) 38:800–11. doi: 10.1007/s10753-014-9990-2
50. Katsha A, Soutto M, Sehdev V, Peng D, Washington MK, Piazuelo MB, et al. Aurora kinase A promotes inflammation and tumorigenesis in mice and human gastric neoplasia. Gastroenterology. (2013) 145:1312–22. doi: 10.1053/j.gastro.2013.08.050
51. Li M, He F, Zhang Z, Xiang Z, Hu D. CDK1 serves as a potential prognostic biomarker and target for lung cancer. J Int Med Res. (2020) 48:300060519897508. doi: 10.1177/0300060519897508
52. Kuang Y, Guo W, Ling J, Xu D, Liao Y, Zhao H, et al. Iron-dependent CDK1 activity promotes lung carcinogenesis via activation of the GP130/STAT3 signaling pathway. Cell Death Dis. (2019) 10:297. doi: 10.1038/s41419-019-1528-y
53. Chu J, Wu S, Xing D. Survivin mediates self-protection through ROS/cdc25c/CDK1 signaling pathway during tumor cell apoptosis induced by high fluence low-power laser irradiation. Cancer Lett. (2010) 297:207–19. doi: 10.1016/j.canlet.2010.05.013
54. Tang Q, Li W, Zheng X, Ren L, Liu J, Li S, et al. MELK is an oncogenic kinase essential for metastasis, mitotic progression, and programmed death in lung carcinoma. Signal Transduct Target Ther. (2020) 5:279. doi: 10.1038/s41392-020-00288-3
55. Saldivar JC, Hamperl S, Bocek MJ, Chung M, Bass TE, Cisneros-Soberanis F, et al. An intrinsic S/G2 checkpoint enforced by ATR. Science. (2018) 361:806–10. doi: 10.1126/science.aap9346
56. Sisinni L, Maddalena F, Condelli V, Pannone G, Simeon V, Li Bergolis V, et al. TRAP1 controls cell cycle G2-M transition through the regulation of CDK1 and MAD2 expression/ubiquitination. J Pathol. (2017) 243:123–34. doi: 10.1002/path.4936
57. Yamazaki H, Kasai S, Mimura J, Ye P, Inose-Maruyama A, Tanji K, et al. Ribosome binding protein GCN1 regulates the cell cycle and cell proliferation and is essential for the embryonic development of mice. PLoS Genet. (2020) 16:e1008693. doi: 10.1371/journal.pgen.1008693
58. Chen X, Niu H, Yu Y, Wang J, Zhu S, Zhou J, et al. Enrichment of Cdk1-cyclins at DNA double-strand breaks stimulates Fun30 phosphorylation and DNA end resection. Nucleic Acids Res. (2016) 44:2742–53. doi: 10.1093/nar/gkv1544
59. Kabeche L, Nguyen HD, Buisson R, Zou L. A mitosis-specific and R loop-driven ATR pathway promotes faithful chromosome segregation. Science. (2018) 359:108–14. doi: 10.1126/science.aan6490
60. Mo F, Zhuang X, Liu X, Yao PY, Qin B, Su Z, et al. Acetylation of Aurora B by TIP60 ensures accurate chromosomal segregation. Nat Chem Biol. (2016) 12:226–32. doi: 10.1038/nchembio.2017
61. Lee HS, Park YY, Cho MY, Chae S, Yoo YS, Kwon MH, et al. The chromatin remodeller RSF1 is essential for PLK1 deposition and function at mitotic kinetochores. Nat Commun. (2015) 6:7904. doi: 10.1038/ncomms8904
62. Lorente L, Martín MM, González-Rivero AF, Pérez-Cejas A, Cáceres JJ, Perez A, et al. DNA and RNA oxidative damage and mortality of patients with COVID-19. Am J Med Sci. (2021) 361:585–90. doi: 10.1016/j.amjms.2021.02.012
63. Uitdehaag JCM, de Man J, Willemsen-Seegers N, Prinsen MBW, Libouban MAA, Sterrenburg JG, et al. Target residence time-guided optimization on TTK kinase results in inhibitors with potent anti-proliferative activity. J Mol Biol. (2017) 429:2211–30. doi: 10.1016/j.jmb.2017.05.014
64. Chen J, Wu R, Xuan Y, Jiang M, Zeng Y. Bioinformatics analysis and experimental validation of TTK as a biomarker for prognosis in non-small cell lung cancer. Biosci Rep. (2020) 40:BSR20202711. doi: 10.1042/BSR20202711
65. Tsai YM, Wu KL, Chang YY, Hung JY, Chang WA, Chang CY, et al. Upregulation of Thr/Tyr kinase increases the cancer progression by neurotensin and dihydropyrimidinase-like 3 in lung cancer. Int J Mol Sci. (2020) 21:1640. doi: 10.3390/ijms21051640
Keywords: lung adenocarcinoma, COVID-19, curcumol, network pharmacology, transcriptome
Citation: Yang L, Xiong H, Li X, Li Y, Zhou H, Lin X, Chan TF, Li R, Lai KP and Chen X (2022) Network Pharmacology and Comparative Transcriptome Reveals Biotargets and Mechanisms of Curcumol Treating Lung Adenocarcinoma Patients With COVID-19. Front. Nutr. 9:870370. doi: 10.3389/fnut.2022.870370
Received: 06 February 2022; Accepted: 16 March 2022;
Published: 19 April 2022.
Edited by:
Kuo-Cheng Lu, Fu Jen Catholic University Hospital, TaiwanReviewed by:
Akshay Narkar, United States Food and Drug Administration, United StatesUma Shanker Navik, Central University of Punjab, India
Copyright © 2022 Yang, Xiong, Li, Li, Zhou, Lin, Chan, Li, Lai and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Rong Li, bGlyb25nMTI3OCYjeDAwMDQwOzE2My5jb20=; Keng Po Lai, Z2xtdV9rZW5ncGxhaSYjeDAwMDQwO3llYWgubmV0; Xu Chen, Y2hlbnh1JiN4MDAwNDA7Z2xtYy5lZHUuY24=; Z2xtdV9jaGVueHUmI3gwMDA0MDt5ZWFoLm5ldA==
†These authors have contributed equally to this work