CDK5 Functions as a Tumor Promoter in Human Lung Cancer.

Cyclin-dependent kinase 5 (CDK5), an atypical member of the cyclin-dependent kinase family, plays an important role in the nervous system. Recent studies have shown that CDK5 is also associated with tumors. However, few studies have been done to investigate the mechanism underlying the connection between CDK5 and cancers. To explore the role of CDK5 in cancers by using an extensive bioinformatics data mining process. We mined the transcriptional, survival, functions and structure of CDK5 gene through databases and in vitro experiments. We found that higher CDK5 expression levels in most cancer cell lines while lower expression in liver and brain cancer cell lines. High expression of CDK5 was associated with shorter overall survival (OS) in lung cancer. In addition, high expression level of CDK5 promoted lung cancer cells proliferation and metastasis. Inhibited CDK5 decreases CAP1 phosphorylation. CDK5 may prove to be a valid target of anticancer therapies.

Introduction

Cancer is a critical cause of death globally 1. Although standard chemotherapy and better supportive care have improved the length of survival and quality of life, the prognosis of advanced cancer patients remains poor. Thus, it is important to investigate underlying mechanisms of cancer tumorigenesis and tumor progression, and identify potential prognostic biomarkers that could even be used as drug targets. Deregulation of cell cycle is a fundamental process that underlies cancer proliferation 2. Progression through cell cycle is regulated by coordinated actions of CDKs, which including CDK5 as the major one 3.

CDK5 is considered as a neuron-specific kinase in the past decade due to the abundant existence of its activator p35 in post-mitotic neurons 4. Increasing evidence has demonstrated that the karyoplasm localization of CDK5 is important for its multiple pathological and physiological functions, including induction of cell motility, apoptosis, cell cycle progression, neuronal migration, neuronal cell survival, lymphatic system, vascularization, and insulin secretion 5-13. Recent studies have shown that CDK5 also participates in a series of biological and pathological processes in non-neuronal cells, and is generally dysregulated in various cancer cells 14.

In this study, we first explored the role of CDK5 in cancers by using an extensive bioinformatics data mining process to find out the expression of CDK5 miRNA in different cancers, and then predict the connection between the expression level and OS. The protein network of predicted associations for CDK5 and alterations in cancer genomics were analyzed. Some experiments were performed to validate the reliance of database and the connection between CDK5 and lung cancer.

Materials and Methods

Oncomine database analysis

The transcription level of CDK5 gene in different cancers was identified by Oncomine database (https://www.oncomine.org/resource/main.html) 15. The mRNA expression level in cancer tissue compared to the normal was obtained as the parameters of p-value<1E-4, fold change>2 and top gene rank 10%, the analyses were summarized in Table 1.

Table 1

CDK5 expression in lung cancer

Cancer subtype	P-value	Fold change	Rank(10%)	Sample	Reference
Lung Carcinoid Tumor	2.72E-7	7.946	10	37	16
Large Cell Lung Carcinoma	1.00E-8	2.122	10	84	17
Lung Adenocarcinoma	2.30E-7	2.168	10	39	18
Lung Adenocarcinoma	9.26E-6	2.203	10	57	19

PrognoScan database analysis

The relationship between CDK5 expression and survival in various types of cancers was analyzed by PrognoScan database (http://www.prognoscan.org/) 20. The threshold was set as cox p-value<0.05 and the analyses were summarized in Fig. 2, Table 2 and Table S2.

Figure 2

CDK5 gene in lung, skin and blood cancer (Prognoscan Database). (A, B) Lung cancer, (C) skin cancer, (D) blood cancer was plotted from PrognoScan database as the threshold of cox p-value < 0.05. The survival curve comparing the patient with high (red) and low (blue) expression was plotted from PrognoScan database.

Table 2

The relationship between transcription level of CDK5 and OS in cancer patients

Cancer type	N	Cox p-value	HR	Endpoint	Dataset	Probe ID
Lung	104	0.000846	2.09	OS	HG_U133A	204247_S_at
	204	0.004384	1.15	OS	HG_U133_Plus_2	204247_S_at
Skin	38	0.001214	1.30	OS	HG_U133_Plus_2	204247_S_at
Blood	53	0.015104	-0.97	OS	HG_U133A	204247_S_at

Kaplan-Meier plotter database analysis

The Kaplan Meier plotter is capable to assess the effect of 54,675 genes on survival using 10,461 cancer samples. These include 5,143 breast, 1,816 ovarian, 2,437 lung and 1,065 gastric cancer patients. The correlation of CDK5 expression and survival in breast, gastric, ovarian and lung cancer was analyzed by Kaplan-Meier plotter database (http://kmplot.com/analysis/) 21 and the result was showed in Fig. 3. The hazard ratio with 95% confidence intervals and log-rank p-value was also computed.

Figure 3

CDK5 gene in ovarian, breast, lung, gastric cancer (Kaplan-Meier Plotter). (A) Ovarian, (B) Breast, (C) Lung, (D) Gastric cancer. The survival curve comparing the patient with high (red) and low (black) expression was plotted from Kaplan-Meier Plotter database. The p-value was 0.73, 0.43, 0.0013 and 0.32 respectively.

STRING database analysis

We analyzed the functional protein association network by STRING database (https://string-db.org/) 22. Protein name: CDK5 and organism; homo sapiens were analyzed. The results showed in Fig. 4.

Figure 4

Identification of known and predicted structural proteins essential for CDK5 function (STRING). Interacting nodes are displayed in colored circles using String. Predicted functional partners of CDK5 are shown based upon peer reviewed published data and curated database entries.

cBioPortal database analysis

The mutations and copy number changes were analyzed by cBioPortal database (http://www.cbioportal.org/) 23. We selected 13 small cell lung cancer and non-small cell lung cancer datasets. Data Type Priority: Mutation and CAN, Gene Set: CDK5, CDK5R1, CDK5R2, MAPT were analyzed. The results showed in Table 3 and Fig. 5-6.

Table 3

The alteration frequency of summary for CDK5, CDK5R1, CDK5R2 and MAPT genes in lung cancer subtypes

Cancer	Data source	N	Frequency (%)	Amplification (%)	Deletion (%)	Mutation (%)	Multiple alteration (%)
Lung squ	TCGA, Provisional	178	7.9	2.2	2.2	2.8	0.6
Lung squ	TCGE,Nature 2012	178	6.2	0.6	1.7	3.4	0.6
NSCLC	TCGA,Nat Genet 2016	1144	5.9	2.2	0.4	3.1	0.2
Lung adeno	TCGA,Nat Genet 2014	230	5.7	3.9		1.3	0.4
Lung adeno	TCGA, Provisional	230	5.7	3.9		1.3	0.4

Figure 5

Copy number alteration of CDK5 gene and lung cancer (cBioPortal). (A) The alteration frequency of CDK5, CDK5R1, CDK5R2 and MAPT genes in lung cancer were analyzed by cBioPortal database. (B) The percentage of alteration in CDK5, CDK5R1, CDK5R2 and MAPT in lung squamous cell carcinoma. (C) CDK5 mutations in lung squamous cell carcinoma. Minimal percentage of alteration was set as 1%. The alteration frequency included amplification (red), deep deletions (blue), mutation (green) or multiple alterations (grey).

Figure 6

Genetic analysis of CDK5, p35, p39 and MAPT in lung cancer subtypes (cBioPortal). Network view of the CDK5 neighborhood in lung cancer. The network contains 54 nodes, including 4 query genes and the 50 most frequently altered neighbor genes (out of a total of 360). The depth of color represents the degree of alteration.

Cell lines and cell culture reagents

Human lung carcinoma cell lines A549 and PC9 were obtained from FDCC (Shanghai, China) and maintained in RPMI 1640 medium or DMEM (High Glucose) (Thermo, Suzhou, China) supplemented with 10% FBS (Gemini, USA) and 1% antibiotics (penicillin-streptomycin, Gibco, USA). Human brain cancer cell line U87 were obtained from ATCC (USA) and maintained in DMEM (High Glucose) (Thermo, Suzhou, China) supplemented with 10% FBS (Gemini, USA) and 1% antibiotics (penicillin-streptomycin, Gibco, USA). Human liver cancer cell line HepG2, breast cancer cell line MCF-7, bladder cancer cell line EJ were obtained from SIBS (Shanghai, China) and maintained in RPMI 1640 or DMEM (High Glucose) (Thermo, Suzhou, China) medium supplemented with 10% FBS (Gemini, USA) and 1% antibiotics (penicillin-streptomycin, Gibco, USA). Human retinal pigment epithelium (ARPE-19) cells were bought from iCell Bioscience (Shanghai, China) and maintained in DMEM/F12 cell medium (Hyclone, USA) supplemented with 10% FBS (Gemini, USA) and 1% antibiotics (penicillin-streptomycin, Gibco, USA). MCF-10A were obtained from SIBS (Shanghai, China) and maintained in DMEM (Invitrogen, USA) supplemented with 20% Horse Serum (Invitrogen, USA), 0.02% EGF (Peprotech, USA) , 0.05% Hydrocortisone (Sigma, USA), 0.1%Insulin (Sigma, USA) and 1% antibiotics (penicillin-streptomycin, Gibco, USA). Roscovitine was from Selleck (Shanghai, China).

Transfection

CDK5 siRNAs were purchased from Sangon Biotech (Shanghai, China). The sequences of the S1, S2 and NC sense strand were 5'- CAGAACCUUCUGAAGUGUAAC-3', 5'- GGGCUGGGAUUCUGUCAUATT-3', and 5'-UUCUCCGAACGUGUCACGUTT-3' respectively. Cells were seeded into six-well plates (Corning, USA) at 1x105/ well and cultured overnight to 40-50% confluence. The cells were transfected with CDK5 siRNA by Lipofectamine 2000 (Invitrogen, USA) according to the manufacturer's instructions. Silencing of CDK5 was examined by Western blots after transfected for 48h.

Antibodies and Western blotting

The A549 cells were homogenized in RIPA lysis buffer containing protease inhibitor mixture (meilunbio, Shanghai, China) on ice for 30 min, then removed from six-well plates (Corning, USA) using a cell scraper and centrifuged for 10 min at 12,000 rpm in a 4℃ Eppendorf microfuge (Germany). Loading buffer was added to the protein lysate and boiled for 10 minutes. Equal amounts of protein were loaded, separated by 10% SDS-PAGE, and transferred to NC membranes (Merck Millipore, Germany). The membranes were incubated first in a 5% powdered milk solution (BBI, Shanghai, China) for 1 hour, then incubated with primary antibodies overnight 4℃, followed by appropriate secondary antibodies. The following primary antibodies were used by CAP1 (1:2000; Santa Cruz, USA), pCAP1 (1:500, donated by Professor Field 24), CDK5 (1:1000, Cell Signaling Technology, USA), β-Actin (1:1000, Sigma, USA) and GAPDH (1:2000, Bioworld, USA).

MTT cell proliferation assay

Approximately 1000 cells were seeded into each well of a 96-well plate (Corning, USA). At each time point, the cells were incubated with 10μl Dye solution (Solarbio, Beijing, China) for 4 hours at 37 °C followed by cytolysis. The released formazan product was detected with an Epoch Microplate Spectrometer (BioTek, USA) by absorbance reading at 490 nM.

Colony formation assay

For colony formation assay, 1000 transfected cells were seeded into 6-well plates and cultured for 10 days at RPMI 1640 medium or DMEM (High Glucose) (Thermo, Suzhou, China) supplemented with 10% FBS (Gemini, USA) and 1% antibiotics (penicillin-streptomycin, Gibco, USA). The colonies were imaged and counted after fixed with 75% ethanol and staining with 0.1% crystal violet solution.

Wound healing and transwell assays

Cell migration was detected by wound healing and transwell assays. Cells were seeded into 6-well plates and cultured until the confluence reached 95%. A sterile 10 μl pipette tip was used to generate a scratch through each well. The wound closure was observed after 0h and 24h and photographed under a microscope (Olympus, Tokyo, Japan).

For migration assays, 3x104 transfected cells in serum-free medium were seeded in the upper well of the chambers equipped with 8 μm membranes were used (Corning Incorporated, Corning, NY, USA), and 500 μl of RPMI 1640 medium or DMEM (High Glucose) (Thermo, Suzhou, China) supplemented with 10% FBS (Gemini, USA) and 1% antibiotics (penicillin-streptomycin, Gibco, USA) were added into the lower chambers. After incubated at 37°C for 24 hours, cells in the upper chamber were removed with a cotton tip, the chambers were washed by 1xPBS for 3 times, the cells migrating to the lower surface of the chamber were fixed with 95% ethanol for 10 min, stained with 0.1% crystal violet solution for 10 min, washed for three times, air dried, photographed and counted in five randomly selected fields for each well by a light microscope (Olympus, Tokyo, Japan), cells were counted by ImageJ software (National Institutes of Health, USA).

Statistical analysis

All experiments were performed at least three times. Data are presented as mean ± standard deviation (SD). Statistical analysis was performed with GraphPad Prism (version 6.01, GraphPad Software, USA). Statistical significance was assumed if p≤0.05.

Results

The transcription level of CDK5 in different cancers

Firstly, the mRNA expression of CDK5 in cancer was analyzed. The main functions of Oncomine database are gene expression differential analysis, gene expression and clinical relevance analysis and multi-gene co-expression analysis. Oncomine database was used to detect the expression of CDK5 mRNA in cancer and normal clinical specimens, and the expression of CDK5 mRNA in lung cancer are shown in Table 1, other cancers results are shown in Table S1. The threshold was set as: p-value<1E-4, fold change>2 and top gene rank 10%. Eventually, 162 datasets, including 18,295 samples, were selected. Compared with normal tissues, the mRNA expression of CDK5 was significantly up-regulated in some cancers, and down-regulated in others (Fig. 1). CDK5 was hyper-expressed in bladder, breast, colorectal, head and neck, lung, lymphoma, melanoma, myeloma, ovarian, uterine corpus leiomyoma, yolk sac tumor and seminoma cancers, and hypo-expressed in brain and central nervous system (CNS), leukemia and liver cancers, suggesting that the transcription level of CDK5 mainly depended on the cancer type.

Figure 1

The transcription level of CDK5 in various cancer types (Oncomine Database). (A) The left column in red showed the number of datasets with CDK5 overexpression and the right column in blue show the under expression datasets number after compared clinical specimen of cancer vs. normal tissue. The box plot comparing specific CDK5 expression in normal (left plot) and cancer tissue (right plot) was derived from Oncomine database. The analysis was shown in (B) lung carcinoid tumor, (C) large cell lung carcinoma, (D, E) lung adenocarcinoma relative to normal lung. The threshold set as follow: p-value: 1E-4, fold change: 2, gene rank: 10%.

The relationship between transcription level of CDK5 and prognosis of cancer

The results of the previous database showed that the expression level of CDK5 mRNA in most cancers including bladder, breast, colorectal, head and neck, lung, lymphoma, melanoma, myeloma, ovarian, uterine corpus leiomyoma, yolk sac tumor and seminoma cancers were higher than that in the normal group and lower in brain and CNS, leukemia and liver cancers. Thus, it is necessary to explore the relationship between the transcription level of CDK5 and prognosis of cancer, and clarify whether CDK5 is a cancer promoter or a tumor suppressor. PrognoScan searches the relation between gene expression and patient prognosis across a large collection of publicly available cancer microarray datasets, provides a platform for evaluating potential tumor markers and therapeutic targets.

The PrognoScan database included 13 types of cancers including bladder, blood, brain, breast, colorectal, esophageal, eye, head and neck, lung, ovarian, prostate, skin and soft tissue cancers. The result of analysis revealed that cancers with statistical significance regarding OS were lung, skin and blood cancers, the p-value was 0.000846, 0.004384, 0.001214 and 0.015104 respectively (Table 2). The high transcription level of CDK5 showed poor prognosis in lung and skin cancer patients, the HR were 8.05, 3.15 and 3.68 respectively (Fig. 2 A-C), whereas good prognosis in blood cancer the HR was 0.38 (Fig. 2D). The Kaplan Meier plotter is capable to assess the effect of 54,675 genes on survival using 10,461 cancer samples, which included data from breast, ovarian, lung and gastric cancers. In order to verify the search results of Prognoscan, we searched the Kaplan-Meier database to predict the relationship between the transcription level of CDK5 and prognosis of cancer. Only in lung cancer, the relationship between CDK5 and prognosis was statistically significant (p=0.0013; HR: 1.23), overexpression of CDK5 related to poor OS (Fig. 3C). There was no statistically significant difference in ovarian, breast and gastric cancer, the p-value was 0.73, 0.43 and 0.32 respectively (Fig. 3A.3B.3D). Both of the two databases indicated that CDK5 was associated with poor prognosis in lung cancer.

CDK5 predicted protein-protein interaction analysis

STRING is a database known to be used to predict protein-protein interactions. The co-expression analysis revealed that CDK5 was co-expressed with cyclin-dependent kinase 5 regulatory subunit 1(CDK5R1/p35), cyclin-dependent kinase 5 regulatory subunit 2(CDK5R2/p39), microtubule-associated protein tau (MAPT), cyclin B1(CCNB1), protein phosphatase 1(PPP1R1B), CDK5 and Abl enzyme substrate 1(CABLES1), tumor protein p53 (TP53), cyclin B2(CCNB2), nudE nuclear distribution E homolog(A. nidulans)-like 1(NDEL1) and dihydropyrimidinase-like 2(DPYSL2), whose correlation score was 0.999, 0.994, 0.991, 0.988, 0.986, 0.986, 0.983, 0.982, 0.978 and 0.976 respectively (Fig. 4). p35 and p39 are widely recognized as activators of CDK5, and the p35/CDK5 complex is essential for the development of the nervous system. MAPT is thought to be associated with the formation and maintenance of the nervous system polarity. We found that the above three proteins (CDK5R1, CDK5R2 and MAPT) had the highest correlation with CDK5, and therefore require more in-depth analysis.

CDK5, p35, p39 and MAPT genomic mutation and CNA in lung cancer subtypes

The cBioPortal for Cancer Genomics provides visualization and analysis of large-scale cancer genomics datasets. Genetic alterations in cancer mainly include mutation, amplification, deletion and multiple alterations. Both PrognoScan and Kaplan-Meier plotter database showed that there was a connection between CDK5 and lung cancer. In order to provide a theoretical basis for further research experiments, CDK5, p35, p39 and MAPT genomic changes in lung cancer subtypes (4104 samples) were studied using 13 datasets of cBioPortal database. Minimal percentage altered samples were set as 1%, and the results are shown in Fig. 5A. Five datasets including lung squamous cell carcinoma (lung squ), non-small cell lung cancer (NSCLC) and lung adenocarcinoma (lung adeno) samples showed gene alteration rates were varied from 5.7% to 7.9%. The sample sizes for these datasets were 178, 178, 1144, 230 and 230 respectively. The detailed information is shown in Table 3.

The OncoPrint database was applied to explore the specific alteration in each gene (Fig. 5A). For example, alteration's percentage in CDK5, CDK5R1, CDK5R2 and MAPT gene among lung squamous cell carcinoma varied from 1.1% to 3.0% in individual genes (CDK5, 2.2%; CDK5R1, 1.1%; CDK5R2, 2.2%; MAPT, 3.0%). The alterations included amplification, deep deletion, truncating mutation and missense mutation (Fig. 5B). Mutations page showed that the 48th amino acid in CDK5 was susceptible to mutations, and the mutation type belonged to missense mutations (Fig. 5C).

Genetic analysis of CDK5, p35, p39 and MAPT in lung cancer subtype

The gene network, which could further analyze genes that interacted with CDK5, p35, p39 and MAPT in lung cancer. The network contains 54 nodes, including 4 query genes and the 50 most frequently altered neighbor genes (out of a total of 360). The depth of color represents the degree of alteration. Darker red nodes including TP53, protein phosphatase 1 regulatory inhibitor subunit 2 (PPP1R2), protein kinase C iota(PRKCI), p21 (RAC1) activated kinase 2(PAK2), titin(TTN) and TRAF2 and NCK interacting kinase(TNIK) gene showed an increased frequency of alteration in lung squamous cell carcinoma (Fig. 6).

Expression of CDK5 in human cancer cell lines

The expression of CDK5 was detected in 8 cell lines, which include 6 human cancer cell lines and 2 normal cell lines (Fig.7A). Analysis showed that expression of CDK5 was up-regulated in most cancer cell lines MCF-7, EJ, A549, PC9, while obviously down-regulated in HepG2 and miner down-regulated in U87 cell line, and nearly no expression of CDK5 was detected in normal cell lines. This was in line with the research outcome of Oncomine database. According to Oncomine CDK5 was hyper-expressed in bladder, breast, colorectal, head and neck, lung, lymphoma, melanoma, myeloma, ovarian, uterine corpus leiomyoma, yolk sac tumor and seminoma cancers, and hypo-expressed in brain and central nervous system (CNS), leukemia and liver cancers. This set of data points to the reliability of database and CDK5 was high expression in human lung cancer cell lines.

Figure 7

Expression of CDK5 in human cancer cell lines (A) and CDK5 promoted the proliferation ability of lung cancer cells The expression of CDK5 was detected in six cancer cell lines and two normal cell lines by western blot. β-actin was used as internal control. (B) Proliferation assay of MTT experiments was performed in A549 cells transfected with NC siRNA or CDK5 siRNA, NC: Normal control. Error bars represent the mean ± s.d. of three independent experiments. (C) The growth inhibitory rate of A549 after treated with roscovitine at the concentration for 10um, 20um, 40um. Error bars represent the mean ± s.d. of three independent experiments. (D) Colony-formation assays was used for detecting the proliferation ability in A549 cell line after treated with roscovitine. (# p<0.05; * p<0.01; ** p<0.001).

CDK5 promoted the proliferation ability of lung cancer cells in vitro

To investigate the biological function of CDK5 in lung cancer cells, we transfected A549 and PC9 cells with CDK5 siRNA and pharmacological inhibited CDK5 by small molecule roscovitine. MTT and colony formation assays were performed to validate the impact of CDK5 in proliferation. Experiments showed that S2 interference effect was better than S1, without specific notification CDK5 siRNA means S2.MTT assay showed that after down-regulated the expression of CDK5 by CDK5 siRNA significantly inhibited the proliferation (Fig.7B, Fig.S1A). Roscovitine inhibited the proliferation of A549 and PC9, and showed the IC50 was between 20um and 40um (Fig.7C, Fig.S1B). Colony formation assays showed the similar results, roscovitine significantly inhibited the colony formation. (Fig.7D, Fig.S1C)

CDK5 promoted lung cancer cell migration in vitro

To determine the effect of CDK5 on metastasis of lung cancer cells. A549 and PC9 cells were transfected with CDK5 siRNA or treated with CDK5 specific inhibitor roscovitine then performed transwell invasion and wound healing assay respectively in indicated time. The results showed that transfected lung cancer cells by CDK5 siRNA will inhibit cells migration (Fig.8A, Fig.S2A). A concentration-dependent decrease of migration by Cdk5 inhibitor roscovitine was observed (Fig.8B, Fig.S2B).

Figure 8

CDK5 promoted lung cancer cell migration . (A) Transwell assay was performed in A549 cells transfected with NC siRNA and CDK5 siRNA. NC: Normal control, S1: si-RNA1, S2: si-RNA2. (B) Wound healing assay was performed in A549 cells after treated with roscovitine at the concentration of 20um or 40um. Error bars represent the mean ± s.d. of three independent experiments. (# p<0.05; * p<0.01; ** p<0.001).

Inhibition of CDK5 decreases CAP1 phosphorylation

Recent findings suggested CDK5 as a kinase that phosphorylates the S307/S309 regulatory site on cyclase-associated protein 1 (CAP1; H. Zhang & G.L. Zhou, unpublished results). We tested the impact of CDK5 inhibition on the phosphorylation of CAP1 in lung cancer cell line A549 by performing different assays using the CDK5 inhibitor roscovitine and CDK5 siRNA. We found that phosphorylated CAP1 (pCAP1) was decreased after treatment with CDK5 inhibitor roscovitine and knockdown CDK5 by siRNA in lung cancer cell (Fig. 9). These results support that CDK5 regulates CAP1 phosphorylation in lung cancer cells as well.

Figure 9

CDK5 regulates CAP1 phosphorylation. Western blots analysis showed the protein of pCAP1 was decreased after treatment with CDK5 inhibitor roscovitine and knockdown CDK5 by siRNA in lung cancer cell A549. NC: Normal control, S1: si-RNA1, S2: si-RNA2.

Discussion

The present study aimed at exploring the role of CDK5 in cancers by using an extensive bioinformatics data mining process to find out the expression of CDK5 mRNA in different cancers, and then predict the connection between the expression level and OS. The protein network of predicted associations for CDK5 and alterations in cancer genomics were analyzed. Western blots, cell proliferation and migration experiments were performed to validate the connection between CDK5 and lung cancer.

CDKs, a group of proline-directed serine/ threonine protein kinases, are conserved throughout evolution and can be found in species from saccharomyces cerevisia to humans. In humans, there are 13 different CDKs (CDK1-CDK13) highly expressed in mitotic cells 25. But CDK5 is unusual because it is not typically activated upon binding with a cyclin and does not require T-loop phosphorylation for activation; instead, it is activated through specific binding with p35 or p39 protein, or their respective cleaved counterparts p25 and p29 26. CDK5 is well characterized for its role in the CNS. While increased levels of CDK5 target proteins are being considered as possible biomarkers of specific cancers. Zhou et al. found that CDK5 activated FAK/AKT signaling pathway to generated vasculogenic mimicry in NSCLC, and identification of CDK5 as a key factor regulating migration and metastasis by regulated actin 27. Zhuang et al. revealed that CDK5 promotes human colorectal cancer via ERK5-AP-1 axis 28. Several researchers established that CDK5 enhanced cell-cycle progression and facilitated medullary thyroid carcinoma (MTC) proliferation by phosphorylated retinoblastoma protein (Rb) at Ser807/811 29. R. Dixon Dorand found that after injecting tumor cells that block the activity of CDK5 into mice, more than half of the mice survived, while almost all the mice in the control group died, which verified the CDK5 associated with tumor development 30. The regulation of CDK5 activity is now emerging as a candidate therapeutic target. The inhibition or knockdown of CDK5 has been proven to play an anti-cancer role through various mechanisms, and can synergize the killing effect of chemotherapeutics. The results are similar to our study. Our analysis predicted CDK5 were higher expressed in most cancer patients compared to normal by Oncomine database. In addition, to further explore the relationship between CDK5 and OS in various cancer types by PrognoScan and Kaplan-Meier Plotter databases. However, the results between OS and other cancers were not statistically significant except lung cancer. Those results suggest that a high transcription level of CDK5 was correlated with poor prognosis in lung cancer patients.

The cause of human cancer is mainly due to irreparable structural mutations in cells. These mutations at very specific genomic locations can alter the function of the gene and DNA copy number 31, 32. Pollack et al. found that copy number alterations (CNAs) uncover all gene expression, which could be a critical element in the tumor development 33. Here in identifying CNAs is a new method for linking CNA with the disease phenotype. The present study aimed to determine if the CNAs of the CDK5 correlate with aggressive cancer sub-types, based on the cBioPortal 34, 35. The mutations alteration (2.8-3.4%) was the most common gene alteration in lung squamous cell carcinoma, while in lung adenocarcinoma patients amplifications (3.9%) were the most common gene alteration. Meanwhile, the 48th amino acid in CDK5 was susceptible to missense mutations (Fig. 5C).

Furthermore, the cBioPortal can be used for interactive analysis and visualization of altered networks. According to STRING analysis, CDK5 is mainly activated by p35/p25, p39/p29 (Fig. 4). In addition, CDK5 is also co-expressed with MAPT, CCNB1, PPP1R1B, CABLES1, TP53, CCNB2, NDEL1 and DPYSL2 in lung cancer (Fig. 4). MAPT and NDEL1 related to microtubule, DPYSL2 related to metastasis which subsequent remodeling of the cytoskeleton 36-38. Fig. 6 showed that the network view of the CDK5 in lung squamous cell carcinoma, those results will help researchers to better understand the molecular mechanisms of CDK5 in lung cancer. Serine/threonine protein kinases encoded by PRKCI gene is a member of the protein kinase C family. It plays a role in microtubule dynamics in the early secretory pathway, which promotes metastasis of esophageal cancer via PKCι-SKP2-AKT pathway 39. The p21 activated kinases (PAK) are critical effectors that link Rho GTPases to cytoskeleton reorganization and nuclear signaling 40. TTN gene encodes a large number of muscle related proteins 41. Protein encoded by TNIK gene is a serine/threonine kinase that functions as an activator of the Wnt signaling pathway which plays important roles in carcinogenesis 42. These data analysis results show that CDK5 is related to cell migration and skeleton, providing a direction for further research.

In order to verify the reliance of databases we performed some experiments. Western blots showed that the expression level of CDK5 were up-regulated in most cancer cell lines MCF-7, EJ, A549, PC9, while obviously down-regulated in HepG2 and miner down-regulated in U87 cell line, and nearly no expression of CDK5 was detected in normal cell lines. This was in line with the research outcome of Oncomine database and some previous reports 43-45. In addition, we also found that knockdown of CDK5 inhibited lung cancer cells proliferation and migration (Fig. 7.8). These results indicated that CDK5 may function as a tumor activator in lung cancer. CDK5 has been reported up-regulated and to act as a tumor activator in many human cancers, including colorectal cancer 43, breast cancer 44, pancreatic cancer45. These results showed that CDK5 can be a promising anticancer target.

Caldesmon phosphorylation and actin polymerization can be observed in CDK5 activity 46. CAP1, a protein encoded by the CAP gene, plays an important role in cell movement and morphological changes by acting synergistically with cofilin to regulate cytoskeleton movement 47. As is known to all tumor cell growth, differentiation, apoptosis and motility are accompanied by abnormal phosphorylation of the protein in cell signal transduction 48, 49. Our group has done a long study of CAP1 and found that CAP1 is abnormal in lung cancer cells and lung cancer tissues, and may be related to poor prognosis and the signal transduction of cytoskeleton and movement 50. Compare with normal tissue, CAP1 is overexpressed in lung cancer tissue, particularly in the metastasis state. Clinical and experimental results are consistent with the results of the database analysis 51, 52.

Zhou et al. found that there were 9 CAP1 phosphorylation sites; one of the phosphorylation sites S309 of CAP1 can recognize the sequence of CDKs (S/T-P-X-K/R/H) 53. CDKs are a large family of serine/threonine kinases. Research showed that this sequence is a specific sequence of phosphorylated substrate of CDK5 54. Therefore CDK5 may be potential kinase of CAP1. After treated with the inhibitor of CDK5, the level of phosphorylated CAP1 was decreased (Fig. 9). Therefore, the phosphorylation of CAP1 is regulated by CDK5.

Database is a powerful tool can help us to screen some target molecules that are worth studying and provide more basis for research, which can predict the direction of the experiment and save research funding. Our study aimed to better understand the molecular mechanisms by extensive oncogenic databases. This study provides a theoretical basis for clinicians and researchers, and finds out a new signaling pathway or biomarker in cancer, which may help develop novel therapeutic approaches for early intervention in cancer prevention. However, this study still has some limitations, the role of CDK5 needed to be further clarified. We intend to further verify the role and function of CDK5 in lung cancer by regulating the phosphorylation site S309 of CAP1.

Conclusion

In summary, this study explored the role of CDK5 in cancers by using an extensive bioinformatics data mining process and seeks new strategies for targeting CDK5 and its downstream mechanisms. CDK5 may prove to be a valid target for anticancer therapies.

Supplementary figures and tables.

Click here for additional data file.