{'fraction_non_alphanumeric': 0.03606864818893517, 'fraction_numerical': 0.02950476481235144, 'mean_word_length': 4.961778471138846, 'pattern_counts': {'":': 0, '<': 1, '<?xml version=': 0, '>': 6, 'https://': 1, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
TLK1‐mediated MK5‐S354 phosphorylation drives prostate cancer cell motility and may signify distinct pathologies

Metastases account for the majority of prostate cancer (PCa) deaths, and targeting them is a major goal of systemic therapy. We identified a novel interaction between two kinases: tousled‐like kinase 1 (TLK1) and MAP kinase‐activated protein kinase 5 (MK5) that promotes PCa spread. In PCa progression, TLK1–MK5 signalling appears to increase following antiandrogen treatment and in metastatic castration‐resistant prostate cancer (mCRPC) patients. Determinations of motility rates (2D and 3D) of different TLK1‐ and MK5‐perturbed cells, including knockout (KO) and knockdown (KD), as well as the use of specific inhibitors, showed the importance of these two proteins for in vitro dissemination. We established that TLK1 phosphorylates MK5 on three residues (S160, S354 and S386), resulting in MK5 activation, and additionally, mobility shifts of MK5 also supported its phosphorylation by TLK1 in transfected HEK 293 cells. Expression of MK5‐S354A or kinase‐dead MK5 in MK5‐depleted mouse embryonic fibroblast (MEF) cells failed to restore their motility compared with that of wild‐type (WT) MK5‐rescued MK5−/− MEF cells. A pMK5‐S354 antiserum was used to establish this site as an authentic TLK1 target in androgen‐sensitive human prostate adenocarcinoma (LNCaP) cells, and was used in immunohistochemistry (IHC) studies of age‐related PCa sections from TRAMP (transgenic adenocarcinoma of the mouse prostate) mice and to probe a human tissue microarray (TMA), which revealed pMK5‐S354 level is correlated with disease progression (Gleason score and nodal metastases). In addition, The Cancer Genome Atlas (TCGA) analyses of PCa expression and genome‐wide association study (GWAS) relations identify TLK1 and MK5 as potential drivers of advanced PCa and as markers of mCRPC. Our work suggests that TLK1–MK5 signalling is functionally involved in driving PCa cell motility and clinical features of aggressiveness; hence, disruption of this axis may inhibit the metastatic spread of PCa.

Metastases account for the majority of prostate cancer (PCa) deaths, and targeting them is a major goal of systemic therapy. We identified a novel interaction between two kinases: tousled-like kinase 1 (TLK1) and MAP kinase-activated protein kinase 5 (MK5) that promotes PCa spread. In PCa progression, TLK1-MK5 signalling appears to increase following antiandrogen treatment and in metastatic castration-resistant prostate cancer (mCRPC) patients. Determinations of motility rates (2D and 3D) of different TLK1-and MK5-perturbed cells, including knockout (KO) and knockdown (KD), as well as the use of specific inhibitors, showed the importance of these two proteins for in vitro dissemination. We established that TLK1 phosphorylates MK5 on three residues (S160, S354 and S386), resulting in MK5 activation, and additionally, mobility shifts of MK5 also supported its phosphorylation by TLK1 in transfected HEK 293 cells. Expression of MK5-S354A or kinase-dead MK5 in MK5-depleted mouse embryonic fibroblast (MEF) cells failed to restore their motility compared with that of wild-type (WT) MK5-rescued MK5 À/À MEF cells. A pMK5-S354 antiserum was used to establish this site as an authentic TLK1 target in androgen-sensitive human prostate adenocarcinoma (LNCaP) cells, and was used in immunohistochemistry (IHC) studies of age-related PCa sections from TRAMP (transgenic adenocarcinoma of the mouse prostate) mice and to probe a human tissue microarray (TMA), which revealed pMK5-S354 level is correlated with disease progression (Gleason score and nodal metastases). In addition, The Cancer Genome Atlas (TCGA) analyses of PCa expression and genome-wide association study (GWAS) relations identify TLK1 and MK5 as potential drivers of advanced PCa and as markers of mCRPC. Our work suggests that TLK1-MK5 signalling is functionally involved in driving PCa cell motility and clinical features of aggressiveness; hence, disruption of this axis may inhibit the metastatic spread of PCa.

Introduction
Prostate cancer (PCa) is the 2nd leading death causing malignancy in men in the United States [1]. Like other malignancies, metastasis has been the primary cause of PCa-related deaths. PCa patients with distant organ metastasis have a 5-year survival rate of only around 30% [2], and currently, no curative therapy is available for the treatment of metastatic PCa [3]. Androgen deprivation therapy (ADT) has been the gold standard method to treat PCa. However, after 1-2 years of treatment, the tumour relapses with more aggressive phenotype, known as metastatic castration-resistant PCa (mCRPC). Cellular motility and invasion are the prerequisites for the cancer cells to migrate and localize in the adjacent and distant organs. Cell migration and invasion are multifactor processes involving numerous signalling pathways. Thus, identification of these factors and underlying signalling cascades is required to devise a therapeutic intervention, which may reduce the spread of the cancer cells.
Tousled kinase (TSL) was first discovered in Arabidopsis thaliana and found crucial for leaf and flower development of the plant [4,5]. Later, homologues of TSL were identified in most metazoans, which are named as tousled-like kinase (TLK) [6]. TLKs are serine/threonine kinases, with maximal activity in the Sphase of the cell cycle, play roles in replication, transcription, chromosomal segregation, DNA damage response and repair. In mammals, two genes are present in the tousled-like kinase family-Tlk1 and Tlk2, which share 89% homology in their entire amino acid sequence and 94% similarity in their C-terminal kinase domain (rev. in [7]). Each gene encodes several isoforms of which TLK1B, a spliced variant with long untranslated region and two regulatory uORFs, is well characterized by our group [8]. Histone H3, histone chaperone ASF1, RAD9 and NEK1 are some of the well-validated substrates of TLK1, which can be directly phosphorylated by TLK1/1B. The activating phosphorylation of NEK1 on T141 residue by TLK1 results in the initiation of DNA damage response signalling through TLK1 > NEK1 > ATR > CHK1 axis that contributes to castration-resistant progression of PCa [9,10]. In addition, recent findings suggest that TLK1-mediated activation of NEK1 also leads to apoptotic prevention and YAP stabilization, which facilitates cancer cell survival and conversion to androgen independence [11,12]. Recently, two research groups independently reported that TLK2 enhances the migration rate and invasiveness of breast cancer and glioblastoma cell lines and hence, can be considered a factor for aggressiveness [13,14]. Another study revealed that the gene encoding tousled-like kinase (Tlk) in Drosophila (Drosophila has only one Tlk gene) is required for the collective migration of border cells by activating JAK/STAT signalling pathway in the ovary, a phenomenon that is reminiscent of cancer migration/metastasis [15]. Since TLK1 and TLK2 have higher sequence similarity and function similarly, we hypothesized that TLK1 may also play role in the dissemination of the cancer cells.
Here, we describe the potential implications of TLK1-MK5 signalling in prostate cancer cell motility and whether disruption of this axis can inhibit PCa cell migration. Our data suggest that MK5 is a direct substrate of TLK1, which can be phosphorylated on several serine residues. We also establish MK5 as a genuine promotility factor in PCa cell lines and suggest that TLK1 phosphorylation of MK5 may enhance its promotility activity. Our bioinformatic study also suggested a correlation of TLK1 and MK5 expression with the aggressiveness of prostate tumours, which hints towards a clinical relevance of TLK1-MK5 signalling in PCa metastasis. In addition, IHC analysis of a PCa TMA revealed that high pMK5-S354 stain correlated with the presence of local nodal involvement in biopsies.

Cell culture
Human embryonic kidney 293 (HEK 293) and mouse embryonic fibroblast (MEF) cells were cultured in DMEM supplemented with 10% FBS and 1% penicillin/streptomycin. PCa cell lines-LNCaP, C4-2B, PC3, DU145 and 22Rv1-were cultured in RPMI 1640 supplemented with 10% FBS and 1% penicillin/ streptomycin. VCaP cells were cultured in 1% fibronectin-coated flask in DMEM supplemented with 10% FBS and 1% penicillin/streptomycin. MK5 À/À MEF cells were a kind gift from M. Gaestel, University of Hannover, Germany [43]. All the other cell lines were purchased from American Type Culture Collection (ATCC, Manassas, VA, USA), and all were authenticated within the past three years. All the cells were maintained in humidified incubator at 37°C with 5% CO 2 .

Cell transfection
Wild-type MEF (MEF WT) and LNCaP cells were transfected with eGFP-MK5 by Lipofectamine 3000 (Thermo Scientific, Waltham, MA, USA, cat# L3000-015) following the manufacturer's protocol, and GFPpositive cells were sorted by flow cytometry. MK5 À/À MEF cells were transfected with eGFP-MK5 or TLK1 mammalian expression plasmid by Lipofectamine 3000 and sort-gated by flow cytometry or selected by 1-2 µgÁmL À1 puromycin treatment for seven days. LNCaP cells were also transfected with kinase-dead TLK1 plasmid and selected by 400-500 µgÁmL À1 G418 for seven days. HEK293 cells were first transfected with eGFP-MK5 and selected by flow cytometry. MK5-positive cells were further transfected with either TLK1 or kinase-dead TLK1 plasmid and selected by antibiotic treatment. To knock down TLK1, MEF WT cells were plated as 40 000 cellsÁwell À1 in a 6-well plate and grown until 60-70% confluency. Cells were transfected twice with 100 nM of mouse-specific TLK1 siRNA (Thermo Scientific, Waltham, MA, USA, cat# AM16708) by Lipofectamine 3000. 72 h after the 2nd dose of siRNA treatment, cells were split and harvested for scratch wound repair assay and WB analysis.

Scratch wound repair assay
Different MEF derivatives and PCa (C4-2B, PC3, DU145 and 22Rv1) cells were seeded as 40 000 cellsÁwell À1 in a 96-well ImageLock plate (Essen Biosciences, Inc., Ann Arbor, MI, USA, cat# 4379) in 10% FBS-containing media. LNCaP cells were seeded as 80 000-90 000 cellsÁwell À1 in a 96-well ImageLock plate. After 48 h when 90-100% confluency was achieved, cells were scratched by a wound maker (Essen Biosciences, Inc., cat# 4493), washed twice with serum free media and supplemented with low serum (0.5-1% FBS) media with or without treatment. Media were supplemented with either DMSO as vehicle control (VC) or 10 µM J54 [44] or 5-20 µM of GLPG0259 (Medkoo Biosciences, Inc., Morrisville, NC, USA, cat# 561481). The plates were then transferred to a IncuCyte Zoom incubator with 5% CO 2 at 37°C, and the migration of the cells was monitored by IncuCyte Zoom live cell imaging system (Essen Biosciences, Inc., Ann Arbor, MI, USA). Machinecontrolled images were taken every 4 h to determine the relative wound density of each well. Mean rate (MR) was determined by calculating the slopes of the curve.

3D chemotactic trans-well migration assay
Both sides of the porous membrane of the top insert of IncuCyte ClearView 96-well chemotaxis plate (Essen Biosciences, Inc., cat# 4582) were coated with either fibronectin (1 mgÁmL À1 ), or Matrigel (50 µgÁmL À1 ), or collagen I (50 µgÁmL À1 ) 3 h before seeding the cells. Cells were seeded in the top chamber as 1000 cellsÁwell À1 in low serum (0.5% FBS) media. Chemotactic gradient was created by putting 10% FBScontaining media in the bottom reservoir. Cells were transferred to a IncuCyte Zoom incubator with 5% CO 2 at 37°C, and the migration of the cells towards chemotactic gradient was monitored by IncuCyte Zoom live cell analysis system. Images were taken at every 4 h of both top insert and bottom reservoir, and total cell occupancy in the bottom reservoir normalized to initial top value was determined. Mean rate (MR) was determined by calculating the slopes of the curve.

Proliferation assay
Cell proliferation was determined by IncuCyte Zoom live cell analysis system from Essen Bioscience. 2000 cells /well were seeded in a 96-well plate with 10% FBS-containing media. Cells were incubated in the IncuCyte Zoom incubator at 37°C with 5% CO 2 . Images were taken at every 4 h, and cell occupancy percentage over time was determined. MTS proliferation assay was conducted by seeding 3000 cells/well in a 96-well plate following the manufacturer's protocol (Promega, Madison, WI, USA, cat# G3580).

Co-immunoprecipitation (co-IP)
200 µg of cell lysate was incubated with TLK1-or GFP-specific antibodies for 4 h on ice. 50 µL of preequilibrated protein A/G agarose beads (50% slurry) (SCBT, Dallas, TX, USA, cat# sc-2003) was added and incubated overnight at 4°C with rotation. Beads were washed three times with cell lysis buffer, and bound proteins were eluted with 1X Laemmli buffer. The eluted sample was further subjected to WB analysis.

Protein purification
Recombinant full-length his-tagged TLK1B and histagged MK5 protein were purified according to the published literature [12,40].

In vitro kinase assay
Radioactive in vitro kinase (IVK) assay was conducted by incubating HSP27 with either TLK1B or MK5 or with both TLK1B and MK5, as previously described [12]. For mass spectrometric (MS) analysis, another IVK assay was conducted by incubating either MK5 or TLK1 with nonradiolabelled ATP. MK5 gel bands were excised and sent to the Kentucky MS facility.

Identification of MK5 phosphorylation by mass spectrometry
Identification of the phosphoresidues of MK5 by TLK1 was conducted at the University of Kentucky Proteomics Core Facility as previously described [12]. . GFP signal distribution was monitored using a Nikon Ti2 fluorescent microscope.

Immunohistochemistry
For TRAMP mouse IHC staining, prostates were harvested, formalin-fixed, processed and paraffinembedded. Prostate TMA was prepared as previously described [10], and its layout is outlined in Table TS3. Tumour tissues were cut in serial 5-µm-thin sections, deparaffinized in xylene, rehydrated with decreased concentrations of ethanol and boiled in unmasking solution (10 mM sodium citrate + 1 mM EDTA) for 15 min for antigen retrieval. Cellular peroxidase was quenched with 3% hydrogen peroxide. Afterwards, tissue sections were blocked in 3% BSA for 1 h followed by the incubation of primary antibody against pMK5 Ser354 (1:200 dilutions) overnight at 4°C. Sections were washed three times and incubated in secondary antibody (Vector Laboratories, Burlingame, CA, USA, cat# PK-6200) for 1 h at room temperature. After washing, the sections were then incubated in Vectastain ABC reagent (Vector Laboratories, Burlingame, CA, USA, cat# PK-6200) for 30 min. The sections were washed again and incubated in DAB substrates for 2 min (Vector Laboratories, Burlingame, CA, USA, cat# SK-4105). Slides were then washed, counterstained with haematoxylin, dehydrated and coverslipped. Imaging and quantification were done by our pathologist, J. King.

Western blotting
Cell were harvested and lysed in 1X RIPA buffer by sonication. Protein concentration was measured using Pierce BCA Protein Assay Kit (Thermo Scientific, Waltham, MA, USA, cat# 23225). Samples were run in SDS/PAGE and transferred to a PVDF membrane. After blocking in 5% nonfat dry milk, the membranes were incubated in primary antibodies overnight at 4°C followed by the incubation in HRP-conjugated secondary antibodies for 1 h at room temperature. The blots were incubated in ECL substrates (Thermo Scientific, Waltham, MA, USA, cat# 32106), and imaging was conducted using Bio-Rad ChemiDoc Imaging System (Bio-Rad, Hercules, CA, USA, cat# 12003154).

Bioinformatics analysis
Genomic amplification and mRNA expression of both TLK1 and MK5 based on regional lymph node metastases and Gleason scores were generated interrogating The Cancer Genome Atlas Prostate Adenocarcinoma (TCGA-PRAD) database using cBioPortal (accessed on 7 February 2020) and UALCAN online platform (accessed on 7 February 2020) [45]. TCGA-PRAD data set can also be downloaded from TCGA data portal (https://tcga-data.nci.nih.gov/tcga/tcgaHome2. jsp). The Kaplan-Meier survival plot of PCa patients based on high and low TLK1 expression along with different Gleason scores was generated using UAL-CAN online platform (accessed on 24 November 2021).

Statistical analysis
Statistical calculations were performed using GRAPH-PAD PRISM 9 and MICROSOFT EXCEL software. Data quantifications are expressed as meanAE standard error of the mean (SEM). Statistical significance was calculated by 2-tailed Student's t-test when comparing the mean between two groups, or by one-way ANOVA followed by Tukey's post hoc analysis when comparing more than two groups. P-values < 0.05 were considered significant.

Ethics approval and consent to participate
Construction of the PCa TMA was made at LSUHSC-Shreveport following the Institutional Review Board (IRB) approval. The experiments were undertaken with the understanding and written consent of each subject. Furthermore, the identity of each donor was deleted from the construction of the block. The study methodologies conformed to the standards set by the Declaration of Helsinki. Animal studies were approved by our institutional animal care (IACUC) for our previous work published in Ref. [10], and only archival tumour tissue from TRAMP mice [46] was used. No animals were used for the purpose of this work.

Genetic depletion and pharmacologic inhibition of both MK5 and TLK1 result in reduced cell migration rate
We have utilized a novel proteomic screening assay from Life Technologies to identify the interactome of TLK1B [47]. TLK1B is a spliced variant of TLK1 that lacks the first 238 amino acids in the N-terminal domain of TLK1. TLK1B is translationally regulated, shares the identical C-terminal kinase domain and is believed to have similar substrate specificity with TLK1 [8,41]. Briefly, biotinylated TLK1B was used to hybridize 9000 human full-length proteins spotted in duplicates on glass slides, after which fluorophore-conjugated streptavidin was used to generate the signals from the interactions between TLK1 and arrayed proteins. This assay identified 165 human proteins that interact with TLK1 with high confidence that are involved in various cellular functions such as cell cycle regulation, DNA damage repair, DNA replication and cell motility. MK5 was identified as one of the top interactor proteins (supplementary section of Ref. [47]). Since MK5 promotes cellular motility in both nonmalignant (e.g. PC12 [17], HUVEC [16]) and neoplastic (e.g. HeLa [20]) cells, we tested whether MK5 overexpression or depletion alters the motility rate of mouse embryonic fibroblast (MEF). MK5 was overexpressed in wild-type (WT) MEF cells using eGFP-MK5 mammalian construct, and exogenous MK5 level was determined by western blotting (Fig. S1A). We compared the scratch healing rate among MK5 À/À MEF- [43], WT MEF-and MK5overexpressing MEF cells by scratch repair assay. Scratch repair appeared to be the fastest in the MK5overexpressing MEF and slowest in the MK5 À/À MEF cells ( Fig. 1A and Fig. S1B). It is noteworthy that MK5 À/À MEF cells were generated in a mixed genetic background (129/Ola X C57BL/ 6), whereas wild-type MEF cells were isolated from C57BL/6 mice [43]. To rule out the possibility that the discrepancy in the genetic background of MK5 À/À MEF and WT MEF cells can affect their migration rate in scratch wound repair assay, we reconstituted MK5 in the MK5 À/À MEF cells (Fig. S2A). In addition, we tested whether TLK1 alone in the absence of MK5 can enhance their migration rate, by stably overexpressing TLK1 in the MK5 À/À MEF Cells (Fig. S2B). When we compared the wound healing rates among these three cell lines, namely, MK5-rescued MEF, MK5 À/À MEF and TLK1-overexpressing MK5 À/À MEF cells in a scratch repair assay, we observed that while MK5-  rescued cells migrated significantly faster, there was no significant difference in the migration rate between TLK1-overexpressing MK5 À/À MEF and MK5 À/À MEF cells (Fig. 1B and Fig. S2C). The latter results indicate that without MK5, TLK1 cannot exert its function in motility promotion. To confirm our observations from the scratch assay, we utilized the 3D chemotactic trans-well migration assay in the IncuCyte machine using three different extracellular matrices (ECMs). We observed that in fibronectin, MK5-rescued cells migrated significantly faster and no difference in the migration rate of TLK1-overexpressing MK5 À/À vs. MK5 À/À MEF cells (Fig. 1C and Fig. S3). Similar trends in migration were also observed using Matrigel and collagen I (Fig. S4 and Fig. S5). To examine whether reconstitution of MK5 increases the proliferation rate of the cells, which might contribute to enhanced migration rate that we observed in the scratch and chemotactic migration assays, we conducted a proliferation assay between MEF MK5 À/À and MEF MK5 À/À GFP-MK5 cells. MK5 reconstitution did not significantly increase the proliferation rate compared with MK5 À/À MEF cells ( Fig. 1D and Fig. S6). We conducted another scratch experiment with LNCaP cells, an androgen-dependent PCa cell line. LNCaP cells were transfected with either eGFP-MK5 or TLK1 kinase-dead (KD) dominant negative mammalian construct (Fig. S7A) [41]. We observed significantly faster wound healing in the MK5-overexpressing LNCaP cells compared with the TLK1 KD and WT LNCaP cells ( Fig. 1E and Fig. S7B). Overall, these findings suggest MK5 as a promotility factor; in addition, TLK1 alone in the absence of MK5 cannot increase the migration rate of the cells. Emerging evidence suggests that TLK2 can promote cancer invasiveness by enhancing the migration and invasion of breast cancer and glioblastoma cell. Since TLK1 and TLK2 share high homology in the amino acid sequence and are believed to have partly redundant functions, we hypothesize that TLK1 can also regulate cellular motility. To examine whether TLK1 acts a potential factor to regulate cell motility, we took both genetic and pharmacologic approaches. Both siRNA-mediated knockdown and pharmacologic inhibition of TLK1 using small molecule inhibitor (J54) resulted in the reduced scratch healing rate in WT MEF cells compared with the control cells (Fig. S8A,  Fig. 1F and Fig. S8B). J54 is a potent secondgeneration phenothiazine derivative specific for TLK1 [44]. Scratch healing rate of TLK1-depleted and TLK1-inhibited WT MEF cells was similar, and there is no statistical difference between these two experimental groups. Therefore, we opted to use J54 for our next experiments. J54 treatment also significantly reduces the wound healing rate of MK5-rescued MK5 À/À MEF cells compared with the vehicle (DMSO)-treated MK5-rescued cells without affecting proliferation or viability (Fig. 1G, Fig. S8C and Fig. S9). These findings establish both TLK1 and MK5 as promotility factors, and we hypothesize that MK5 may function as a downstream effector of TLK1 to promote cellular motility. Unexpectedly, TLK1 inhibition by J54 also resulted in reduced wound closure in MK5 À/À MEF cells, which suggest that TLK1 may promote cellular migration through additional downstream effectors even in the absence of MK5 (Fig. S10A and Fig. S10B). We also expressed two mutant variants of eGFP-MK5: S354A (a target phosphorylation by TLK1-see below) and K51E kinasedead variant [42] in MK5 À/À MEF cells. Expression of neither protein could restore motility to MK5 À/À MEF cells (Fig. 1H, Fig. S11A,B).

TLK1 interacts and phosphorylates MK5 both in vitro and in cultured cells and increases its catalytic activity
Results obtained from our previously conducted protoarray assay demonstrated that TLK1B interacts with MK5 in vitro [47]. To test whether TLK1 interacts with MK5 in cultured cells, we conducted coimmunoprecipitation (co-IP) experiments targeting both endogenously and exogenously expressed proteins. Immunoprecipitations of endogenous TLK1 using TLK1-specific antibodies brought down endogenous MK5 in LNCaP cells ( Fig. 2A). We overexpressed MK5 in HEK 293 cells using eGFP-MK5 mammalian construct and immunoprecipitated MK5 using GFPspecific antibody, which also precipitated TLK1 (Fig. 2B). These data suggest that TLK1 and MK5 form a complex in cells. We confirmed this observation performing His and GST pulldown experiments independently. Incubation of purified recombinant histagged TLK1B with MK5-overexpressing HEK cell lysates precipitated both endogenous MK5 and exogenously expressed GFP-MK5, as determined by WB (Fig. 2C). Similarly, incubation of purified recombinant GST-tagged MK5 with TLK1 and MK5 coexpressing cell lysate precipitated both TLK1 (Fig. 2D, left panel) and GFP-MK5 (Fig. 2D, right panel). Together, these results indicate that TLK1 interacts directly with MK5. The observation that endogenous MK5 was also pulled down in both co-IP (Fig. 2B) and GST pulldown assay (Fig. 2D, right panel) suggests that endogenous MK5 may dimerize with exogenously expressed GFP-MK5 and/or GST-MK5. We examined whether TLK1 possesses the ability to phosphorylate MK5 in vitro and in cultured cells. We first purified recombinant his-tagged TLK1B and recombinant his-tagged MK5 following the previously published protocol (Fig. S12A,B) [12,40]. We confirmed the catalytic activity of MK5 using MK5specific protein and peptide substrate, HSP27 (Fig. S12C) and PRAKtide (Fig. S12D), respectively, by ADP hunter assay. Incubation of purified recombinant his-tagged TLK1B with both purified recombinant his-tagged MK5 and HSP27 in an IVK assay synergistically increases the phosphorylation level of MK5 and intrinsic catalytic activity of MK5 towards its known substrate, HSP27 (Fig. 2E,F). To determine whether TLK1 can phosphorylate MK5 in cultured cells, we transfected HEK 293 cells either with GFP-MK5 alone or cotransfected the cells either with GFP-MK5 +TLK1 or with GFP-MK5+TLK1 kinase-dead (KD) expression construct. WB analysis revealed a hyperphosphorylated form of GFP-MK5 only in MK5+ TLK1-coexpressed cell lysate that migrated slower than the main form of GFP-MK5 (Fig. 2G). Lambda protein phosphatase (LPP) treatment reduces this slower migrating band, which suggests that TLK1 can phosphorylate MK5 in cultured cells. We also conducted another IVK assay by incubating both recombinant His-TLK1B and GST-MK5 in the presence of nonradioactive cold ATP. The reactions were run in an SDS/PAGE gel, and the corresponding MK5 bands were excised for mass spectrometric analysis.
3.3. TLK1 phosphorylates MK5 in vitro on three unique residues, which are mapped to the functional domains of MK5 To determine the phosphopeptides of MK5 by TLK1, MK5 bands were digested by trypsin and subjected to shot-gun proteomic analysis using an LTQ-Orbitrap mass spectrometer. MS data sets are searched with MASCOT against a custom database containing human MAPKAPK5. The summary of the analysis is as follows: (a) MAP kinase-activated protein kinase 5 (MAPKAPK5 or MK5) was detected with high protein score (4686-8482) and good peptides coverage (70.40-80.55%) in both samples. (b) Phosphorylated sites were detected in both samples. After comparing MK5 mock and MK5+TLK1B samples, a few unique phosphorylation sites were detected in MK5+TLK1B sample: S160, S348, S354 and S386 ( Fig. S13 and Table S1). The constitutive phosphorylation sites present in both MK5 mock and MK5+TLK1B samples can be explained as the autophosphorylation sites of MK5 (Fig. S14). (c) S348 and S354 are within the same peptide with single phosphorylation event. After comparing the spectra, S354 was determined to be the more likely phosphorylation site. (d) Thus, we conclude that there are three unique phosphorylation sites in MK5+TLK1B sample: S160, S354 and S386. We mapped these serine residues to the domains of MK5 and identified that they are located in the kinase domain (in the activation loop near the residue Thr182, which needs to be phosphorylated for the catalytic activation of MK5 [48,49]), in the nuclear localization signal (NLS) domain and in the ERK3 binding domain (rev. in [37]) (Fig. S15).

Antiandrogen treatment increases pMK5 Ser354 level and can be inhibited by small molecule inhibitor specific to TLK1
Antiserum raised against one of the phosphoresidue of MK5 (Ser354), one of the sites that we identified as being target of TLK1, has recently become commercially available, and we tested it to determine whether this specific MK5 phosphoprotein is present in several androgen-dependent and androgen-independent PCa cell lines (Fig. 3A and D). MK5 À/À MEF cells were used as a negative control. We also detected TLK1 in those PCa cell lines (Fig. 3A). Western blotting analysis revealed the presence of pMK5-S354 level relative to the total MK5 in all PCa cell lines tested, while no pMK5-S354 was detected in MK5 À/À MEF cells. We also previously reported that treatment of LNCaP and other androgen-sensitive PCa cells with bicalutamide results in increased expression of TLK1B [9,10]. We now show that antiandrogen treatment also results in a corresponding dose-dependent increase in pMK5 (S354) (Fig. 3B and E), suggesting that TLK1/1B may be the kinase responsible for this phosphorylation. Moreover, we treated LNCaP cells with two specific small molecule inhibitors of TLK1 (THD or J54), which resulted in the reduction of pMK5-S354 level compared with the DMSO-treated control LNCaP cells, which further supports TLK1 role in MK5 Ser354 phosphorylation (Fig. 3C and F).

pMK5
Ser354 is elevated in PCa progression of TRAMP mice and in patients with high-grade metastatic prostate cancer We wanted to study whether the presence of pMK5 presented with a pattern of progressive intensity during progression to invasive PCa in the TRAMP model, given that we have already reported the expression of TLK1B increases greatly from 12 to 30 weeks of age, particularly after castration [10]. At 12 weeks,  TRAMP mice present feature of early carcinogenic progression of the prostate, with a mix of hyperproliferative lesions such as prostatic intraepithelial neoplasia (PIN), to well-differentiated adenocarcinomas, and in some case some areas of invasiveness with breakdown of the basement membrane. These features are remarkably well represented in Fig. 4A (top panels), which clearly shows that the areas of PIN, confined PCa and locally invasive PCa are stained progressively darker upon pMK5-IHC, compared with normal acini (Fig. 4A and B). At 26 weeks, the PCa has progressed to fully invasive, and the cancer cells are spread thought the glands (Fig. 4A, bottom panels). We interrogated our institutional PCa TMA (with African American prevalence) to establish possible clinicopathologic correlations based on the Gleason score/s (GS), regional lymph node metastasis and the levels of pMK5 Ser354 (Table S2). pMK5 Ser354 staining intensity was progressively higher in tumour samples with increasing Gleason scores and lymph nodes metastases. For instance, our representative IHC images revealed higher positive staining of pMK5-S354 of several GS 8-9 tumours compared with GS 6 (Fig. 4C). Similar increase was also observed in the tumours with 1-3 lymph node metastasis (N1). A very notable feature was that regardless of GS (the key token of pathology determination to date), the nuclear-to-cytoplasm ratio of pMK5 was decreased in several samples, including some N1 of not the highest GS grade (Table S3). MK5 is generally localized to the nuclei, but it shuttles to the cytoplasm upon MAPK activation [42,50], where it is believed to regulate actin cytoskeletal reorganization, although the role of stress-responsive p38-MAPK in activating MK5 and its redistribution is debated [42,43]. We have thus investigated the possible role of TLK1-overexpressing cells in the cellular sublocalization of eGFP-MK5 and found that at steady state and in the absence of cellular stress, the distribution of eGFP-MK5 was strongly shifted to the cytoplasm in contrast to the nucleus of cells with normal level TLK1 (Fig. S16). While our TMA does not report on patients' outcome, it would be intriguing to study this in the future to establish whether the pMK5 nuclear/cytoplasm ratio can be used as a prognostic indicator. Generally, tumours with GS≤ 6 are considered as low-risk tumours with well-differentiated PCa cells, GS=7 is considered as moderate risk, and GS≥ 8 is considered as high-risk tumours with poorly differentiated PCa cells that have higher migratory and invasive potentials. Higher phosphorylation of MK5 may suggest TLK1-MK5 signalling as a key driver of metastasis in PCa subjects.

Pharmacologic inhibition of MK5 strongly reduces wound healing of both androgendependent and androgen-independent PCa cells
To establish that MK5 is indeed important for cell motility, we tested a panel of PCa cell lines, including LNCaP, C4-2B, PC3, DU145 and 22RV1 with the specific MK5 inhibitor GLPG0259 (Fig. 5 and Fig. S17) using the scratch repair assay via the IncuCyte. In all cell lines tested, GLPG resulted in a strong immobilization effect with no apparent toxicity in terms of loss of viable cells. This effect was noticeable within few hours after the addition of the drug, with already maximal effect at 5 lM, which is close to the established IC 50 for GLPG0259 [51,52]. We should emphasize that GLPG0259, at the concentrations used in these assays (which are similar to those reported in human studies [51]), did not affect viability and growth. While we have not yet tested GLPG0259 on a cell line such as PC3 in a xenograft model of metastasis, such work was reported recently with very exciting results of metastasis suppression in xenografts of B16 and MDA-MB-231 cells [53].

Bioinformatic analysis revealed genomic amplification and upregulation of both TLK1 and MK5 in metastatic tumours
Previous studies reported TLK1 as one of the major drivers of CRPC progression after ADT [9,10,12]. Androgen ablation also results in the activation of some MAPK pathways, which will further increase tumour aggressiveness. We interrogated publicly available databases to determine TLK1 and MK5 status in actual PCa patients with advanced tumours. Analysis of TCGA, SU2C, Broad/Cornell and other PCa data sets revealed genomic copy-number increase of MK5 in CRPC     patients compared with localized tumours (Fig. 6A). We analysed the mRNA expression of both TLK1 and MK5 of PCa tumours with different nodal metastatic status using the UALCAN online bioinformatic tool. Tumours with 1-3 regional nodal metastatic lesions (N1) showed slightly higher expression of TLK1, but significant upregulation of MK5 compared to tumours with no nodal metastatic lesion (N0) (Fig. 6B and C) was observed. Similarly, consistent upregulation of both TLK1 and MK5 was observed in tumours with increasing Gleason scores ( Fig. 6D and E), although our driving hypothesis is that it is the activity of MK5 (via upregulation by TLK1) that is really more significant. Importantly, the TLK1 gene was identified by coexpression analysis using WGCNA as a key driver of PCa, highly enriched among candidate genes collected from expression quantitative trait loci (eQTL), somatic copynumber alterations (SCNA) and prognostic analyses [54]. Further analysis of TLK1 stratified as high or low expression indicated that it is a predictor of poor outcome, with earlier death time, when the Gleason score is lower and thus deemed as more favourable (Fig. S18).

Discussion
Prostate cancer (PCa) is a slow-progressing disease where the mortality rate peaks with the metastasis of cancer cells to other vital organs. While organconfined PCa is more clinically manageable and has a nearly 100% five-year survival rate, PCa metastasis drastically decreases the survival of the patients. The median survival rate of the prostate cancer patients with liver, lung, bone and lymph node metastasis is approximately 14, 19, 21 and 32 months, respectively [55]. Failure of ADT leads to more aggressive malignancy and drives metastatic progression of the PCa. However, the molecular mechanism of this metastatic dissemination is not fully understood yet. Motility and invasion are two essential cellular processes for cancer cells to get detached from the primary tumour site and disseminated throughout the body. Motility and invasion are regulated by the concerted action of several factors, coordination of multiple signal transduction pathways, gene expression, cytoskeletal changes and remodelling of surrounding extracellular matrices. These allow the cells to invade and migrate through new tissues. Our discovery of TLK1B-MK5 interaction through protoarray assay prompted us to investigate the role of TLK1 and MK5 in motility in both nonmalignant and neoplastic cell lines as TLK1 paralog, TLK2 and another kinase, MK5, were reported to increase cellular motility and invasion [13,14,16,18,20]. Both genetic depletion and pharmacologic inhibition of MK5 and TLK1 independently resulted in significantly reduced wound healing, whereas overexpression/rescue of MK5 enhances the wound closure rate (Figs. 1 and 5). However, overexpression of TLK1 alone in MK5-depleted cells did not promote the migration, which suggests TLK1 must exert its promotility function primarily through MK5 (Fig. 1B,C). We further confirmed TLK1-MK5 interaction by co-IP, His and GST pulldown assay and elucidated that TLK1 can phosphorylate MK5 in three serine residues (S160, S354 and S386), which are located on some functional domains of MK5 (Figs. 2 and S13-S15). Thus, MK5 is a direct substrate of TLK1, which increases its intrinsic kinase activity in vitro (Fig. 2E and 2G). We reported the presence of TLK1 and MK5 phosphoprotein (pMK5-S354) in all major androgen-dependent and androgenindependent PCa cell lines and demonstrated that antiandrogen (bicalutamide) treatment increases the pMK5-S354 level in a dose-dependent fashion (Fig. 3A, B,D and E). Pharmacologic inhibition of TLK1 reduces this pMK5 level, which suggests TLK1 as an authentic pMK5-S354 kinase ( Fig. 3C and 3F). The importance of S354 phosphorylation for the activation of MK5 was then demonstrated by the lack of motility rescue of MK5 À/À MEF cells complemented with a MK5-S354A mutant, which was no better than the K51E kinasedead variant (Fig. 1H and Fig. S11). In this respect, one could argue that TLK1 appears to play a critical role in the activation of MK5 at least for some of its functions. In fact, overexpressing TLK1 resulted in the relocalization of GFP-MK5 to cytoplasm in the absence of a cellular stress (e.g. osmotic) that might involve the activation of p38-MAPK [42], or by invoking the action of PKA (see below).
In PCa, TLK1/1B-MK5 signalling axis might be crucial as the ADT-related stress increases TLK1/1B level and activates compensatory MAPK signalling, which may result in the activation and stabilization of MK5 [9,56]. Through our bioinformatics study using TCGA and other publicly available data sets, we also observed that both TLK1 and MK5 are upregulated in metastatic and advanced PCa (Fig. 6). Cellular abundance of these kinases may increase their interaction frequencies and start a novel signalling transduction pathway to promote cellular migration. When the tumour grows sufficiently large, cells may shed from the primary tumour site and find a secondary location to avoid hypoxic stress. We propose that TLK1-MK5 signalling aids this process by increasing the motility and invasion of the cells. Interrogation of a PCa TMA and IHC of TRAMP mice of different stages of PCa progression leading to metastatic spread is supportive of the notion that pMK5-S354 is potential marker of disease progression and in particular of invasive PCa (Fig. 4). Notably to date, pathologists mainly use H&E staining in the attempt to determine PCa features that predict worse outcome. However, a relatively mild GS is often no predictor of better outcome, particularly if accompanied by the known presence of nodal involvement (N1). A reliable IHC marker, with an established cellular function, that better predicts the intrinsic invasive potential of the PCa cells would be a significant addition to available tools for pathology determinations.
Another potential mechanism could be focal adhesion complex modifications by TLK1-MK5 signalling that promotes cell motility. Since both TLKs and MK5 are known to form complexes with FAK [16,22], Src [13,14,22] and paxillin [22], it is likely that these complexes depend on TLK1-MK5 interaction. TLK2 is reported to complex with Src and activate EGFR/ Src/FAK signalling pathway to promote the migration and invasion in breast adenocarcinoma and glioblastoma cells [13,14]. Yoshizuka et al. demonstrated that MK5 increases motility in endothelial cells by phosphorylating FAK on Tyr397 residue and localizes it to focal adhesions in advanced skin carcinogenesis [16]. Additionally, another group reported that recombinant MK5 can phosphorylate FAK, Src and paxillin in IVK assays and form complexes with FAK and Src and that endogenous MK5 localizes to focal adhesions [22]. We suggest that TLK1 interacts with MK5 and recruits it to focal adhesions where both TLK1 and MK5 may form complex with either FAK or Src or paxillin, phosphorylate and activate them and transduce downstream signals. FAK activation by TLK1-MK5 axis can regulate the affinity and avidity of integrins towards ECM and initiate a feedforward loop of MK5 activation through Rac1 > PAK > ERK3/ 4 > MK5 signalling [65][66][67][68].
MK5 can be activated by both typical and atypical MAPKs, p38 [48] and ERK3/4 [69][70][71], respectively. MK5 full activation needs the phosphorylation of Thr182 residue, which is located in its activation loop [48,70]. One of the phosphorylation of MK5 by TLK1 is mapped to Ser160, and this site is located on the catalytic domain near to the activation loop of MK5. This phosphorylation may cause a conformational change of MK5 and increase its affinity towards its substrates. On the contrary, TLK1 phosphorylation of MK5 may also enhance nucleocytoplasmic localization of MK5, while MK5 is predominantly nuclear in the resting cells. This is, in agreement with the findings of others that phosphorylation of MK5 by its interaction partners causes the redistribution of MK5 from nucleus to cytoplasm [42,50,[72][73][74]. TLK1-mediated phosphorylation on Ser354 residue, which lies on the NLS domain of MK5, may mask the NLS and promote its nuclear export to exert its functions on motility enhancement.

Conclusion
Implementation of ADT and novel antiandrogens has made an impact in the early survival of PCa patients, but the presence of metastasis remains a significant poor prognostic indicator. Inhibition of metastatic dissemination seems to be the first and foremost challenge in the treatment of mCRPC patients. We identified a novel signalling axis (TLK1-MK5), which is active in the majority of the PCa cell lines and demonstrated that disruption of this signalling significantly reduces PCa cell migration. Further research is needed to determine whether the inhibition of TLK1-MK5 signalling can reduce PCa metastasis in xenograft mouse model using TLK1 and/or MK5 inhibitor. Successful completion of such work may result in the application of these inhibitors in combination with other treatment strategies to delay PCa metastatic burden and lethal progression. fellowship to MIK. We like to thank Dr. Erik First (LSUHSC-Shreveport) for his assistance with MK5 protein purification. We like to thank the INLET facility of LSUHSC-Shreveport, especially Dr. Ana maria Dragoi, Dr. Jennifer Carroll-Duhon and Brian Latimer for their assistance in working with the IncuCyte machines. We also like to thank Dr. Jing Chen and Dr. Haining Zhu from the University of Kentucky Proteomics Core Facility for conducting the mass spectrometric analysis.

Supporting information
Additional supporting information may be found online in the Supporting Information section at the end of the article. Fig. S1. MK5 overexpression and image representation of the scratch wound repair assay among different MEF cells. Fig. S2. MK5 rescue and TLK1 overexpression and image representation of the scratch wound repair assay among different MEF cells. Fig. S3. Image representation of 3D chemotactic trans-well migration assay among different MEF cells. Fig. S4. 3D chemotactic trans-well migration assay among different MEF cells using Matrigel coating. Fig. S5. 3D chemotactic trans-well migration assay among different MEF cells using collagen coating. Fig. S6. Image representation of the proliferation assay between MK5 rescued and MK5 À/À MEF cells. Fig. S7. MK5 or kinase dead TLK1 overexpression and image representation of the scratch wound repair assay among different LNCaP cells. Fig. S8. TLK1 knockdown and image representation of scratch wound repair assays among different MEF cells. Fig. S9. Cell proliferation determination among MEF MK5 À/À GFP-MK5 cells by MTS assay. Fig. S10. Wound healing rate between treated or untreated MEF MK5 À/À cells determined by scratch wound repair assay. Fig. S11. Wild-type or mutant variant of MK5 reconstitution and image representation of the scratch wound repair assay among different MEF MK5 À/À cells. Fig. S12. Recombinant protein purification and MK5 kinase activity determination by ADP hunter assay. Fig. S13. Spectral analysis of MK5 phosphoresidue determination by MS. Fig. S14. MK5 phosphopeptide determination by MS. Fig. S15. Schematic diagram of MK5 domains with TLK1 mediated phosphoresidues.    Table S1. Peptides containing unique phosphorylated sites detected in the TLK1-MK5 samples. Table S2. Layout of PCa TMA. Table S3. Clinicopathological features of the TMA tumour samples.