Mechanistic and evolutionary insights into isoform-specific ‘supercharging’ in DCLK family kinases
Abstract
Catalytic signaling outputs of protein kinases are dynamically regulated by an array of structural mechanisms, including allosteric interactions mediated by intrinsically disordered segments flanking the conserved catalytic domain. The doublecortin-like kinases (DCLKs) are a family of microtubule-associated proteins characterized by a flexible C-terminal autoregulatory ‘tail’ segment that varies in length across the various human DCLK isoforms. However, the mechanism whereby these isoform-specific variations contribute to unique modes of autoregulation is not well understood. Here, we employ a combination of statistical sequence analysis, molecular dynamics simulations, and in vitro mutational analysis to define hallmarks of DCLK family evolutionary divergence, including analysis of splice variants within the DCLK1 sub-family, which arise through alternative codon usage and serve to ‘supercharge’ the inhibitory potential of the DCLK1 C-tail. We identify co-conserved motifs that readily distinguish DCLKs from all other calcium calmodulin kinases (CAMKs), and a ‘Swiss Army’ assembly of distinct motifs that tether the C-terminal tail to conserved ATP and substrate-binding regions of the catalytic domain to generate a scaffold for autoregulation through C-tail dynamics. Consistently, deletions and mutations that alter C-terminal tail length or interfere with co-conserved interactions within the catalytic domain alter intrinsic protein stability, nucleotide/inhibitor binding, and catalytic activity, suggesting isoform-specific regulation of activity through alternative splicing. Our studies provide a detailed framework for investigating kinome-wide regulation of catalytic output through cis-regulatory events mediated by intrinsically disordered segments, opening new avenues for the design of mechanistically divergent DCLK1 modulators, stabilizers, or degraders.
eLife assessment
This important study expands on current knowledge of allosteric diversity in the human kinome by C-terminal splicing variants using as a paradigm DCLK1. The authors provide convincing evolutionary and some mechanistic evidence how C-terminal isoform specific variants generated by alternative splicing can regulate catalytic activity by means of coupling specific phosphorylation sites to dynamical and conformational changes controlling active site and substrate pocket occupancy, as well as protein-protein interactions. The data will be of interest to researchers in the kinase and signal transduction field.
https://doi.org/10.7554/eLife.87958.3.sa0Introduction
Protein kinases are one of the largest druggable protein families comprising 1.7% of the human genome and play essential roles in regulating diverse eukaryotic cell signaling pathways (Manning et al., 2002). The doublecortin-like kinases (DCLKs) are understudied members of the calcium/calmodulin-dependent kinase (CAMK) clade of serine-threonine kinases (Agulto et al., 2021; Bayer and Schulman, 2019; Gógl et al., 2019). There are three distinct paralogs of DCLK (1, 2, and 3), the last of which is annotated as a ‘dark’ kinase due to the lack of information pertaining to its function (Berginski et al., 2021). Human DCLK1 (also known as DCAMKL1) was initially identified in 1999 (Sossey-Alaoui and Srivastava, 1999), followed by the cloning of human DCLK2 and -3 paralogs (Ohmae et al., 2006). Full-length DCLK proteins contain N-terminal doublecortin-like (DCX) domains, microtubule-binding elements that play a role in microtubule dynamics, neurogenesis, and neuronal migration (Couillard-Despres et al., 2005; Horesh et al., 1999). DCLKs have garnered much interest as disease biomarkers, since they are upregulated in a variety of cancer pathologies (Cheng et al., 2022a; Gao et al., 2016; Westphalen et al., 2017), as well as neurodegenerative disorders such as Huntington’s disease (Galvan et al., 2018). However, the mechanisms by which DCLK activity is autoregulated, and how and why they have diverged from other protein kinases is not well understood.
Like all protein kinases, the catalytic domain of DCLKs adopts a bi-lobal fold (Gógl et al., 2019), with an N-terminal ATP-binding lobe and C-terminal substrate-binding region. Canonical elements within the two lobes include the DFG motif, a Lys-Glu salt bridge that is associated with the active conformation, Gly-rich loop, and ATP-binding pocket, which are all critical elements for catalysis. Many protein kinases, including CAMKs, tyrosine kinases, and AGCs, elaborate on these core elements with unique N-terminal and C-terminal extensions that flank these catalytic lobes (Kannan et al., 2007; Yeon et al., 2016; Nguyen et al., 2015a; Yeung et al., 2021; Kwon et al., 2018), allowing them to function as allosteric regulators of catalytic activity (Gógl et al., 2019). Indeed, CAMKs are archetypal examples of kinases that can exist in an active-like structural conformation, yet still remain catalytically inactive (Gógl et al., 2019). This is in large part due to the presence of unique C-terminal tails that are capable of blocking ATP or substrate binding in well-studied kinases such as CAMK1 and CAMKII. In canonical CAMKs, autoinhibition may be released upon Ca2+/calmodulin (CaM) interaction with the CAMK C-tail, which makes the substrate-binding pocket and enzyme active site accessible (Rellos et al., 2010). In CAMKII, the N- and C-terminal segments flanking the kinase domain are variable in length across different isoforms and the level of kinase autoinhibition or autoactivation has been reported to be dependent on the linker length (Bhattacharyya et al., 2020). The CAMKII C-tail can be organized into an autoregulatory domain and an intrinsically disordered association domain. The autoregulatory domain also serves as a pseudosubstrate, which physically blocks the substrate-binding pocket until it is competed away by CaM (Hudmon and Schulman, 2002). Notably, this autoregulatory pseudosubstrate can be phosphorylated (Rellos et al., 2010), and phosphorylation of the C-tail makes CAMKII insensitive to CaM binding. Across the CAMK group, several other kinases share autoinhibitory activity via interactions between Ca2+/CaM binding domains and the C-terminal tail (Huse and Kuriyan, 2002; Wayman et al., 2008), and a major feature of these kinases is variation in the tail length across the distinct genetic isoforms.
The human genome encodes four distinct DCLK1 isoforms, termed DCLK1.1–1.4 in UniProt (Table 1, Figure 1; Omori et al., 1998), which display differential activity- and tissue-specific expression profiles. Human DCLK1.1 (also known as DCLK1 alpha) is expressed in a variety of tissues, but is enriched in cells derived from the fetal and adult brain, whereas DCLK1.2 (also known as DCLK1 beta) is expressed exclusively in the embryonic brain (Matsumoto et al., 1999). DCLK1.3 and DCLK1.4, which lack tandem microtubule-binding DCX domains (Figure 1) but are otherwise identical to DCLK1.1 and DCLK1.2, respectively, are also highly expressed in the brain. To aid with clarity, the names of the human DCLK1 genes and their isoforms used in this paper are summarized in Table 1. Recent structural and cellular analyses have begun to clarify the mechanisms by which the DCLK1.2 isoform is regulated by the C-tail (Agulto et al., 2021; Patel et al., 2021; Cheng et al., 2022b). Mechanistically, autophosphorylation of Thr 688, which is present only in the C-tail of DCLK1.2 (and DCLK1.4), blunts kinase activity and subsequently inhibits phosphorylation of the N-terminal DCX domain and thus drives DCLK microtubule association in cells (Agulto et al., 2021). Consistently, deletion of the C-tail or mutation of Thr 688 restores DCLK1.2 kinase activity, subsequently leading to DCX domain phosphorylation and the abolition of microtubule binding. The length and sequence of the C-tail varies across the DCLK1 isoforms; however, how these variations contribute to isoform-specific functions and how they emerged during the course of evolution is not known.
Names of doublecortin-like kinase (DCLK) isoforms discussed in this paper, along with their respective isoform number, UniProt identification, and alternate names that have been used.
Name (this Study) | Isoform number | UniProt ID | Alternate names |
---|---|---|---|
DCLK1.1 | 1 | O15075-2 | DCAMKL1 alpha |
DCLK1.2 | 2 | O15075-1 | DCAMKL1 beta |
ΔDCLK1.1 | 3 | O15075-3 | |
ΔDCLK1.2 | 4 | O15075-4 |

Cartoon schematic of Human DCLK orthologs and respective isoforms.
(A) Schematic representation of domain organization for the known isoforms of the three human doublecortin-like kinase (DCLK) paralogs. Domain boundaries are annotated according to the representative amino acid sequences derived from UniProt. (B) DCLK1 isoforms visualized as cartoons, showing key structural differences between the four human DCLK1 isoforms and a DCLK1 catalytic domain with artificially short linker regions (DCLK1cat).
In this paper, we employ an evolutionary systems approach that combines statistical sequence analysis with experimental studies to generate new models of DCLK evolutionary divergence and functional specialization. We identify the C-terminal tail as the hallmark of DCLK functional specialization across the kingdoms of life and propose a refined model in which this regulated tail functions as a highly adaptable ‘Swiss Army Knife’ that can ‘supercharge’ multiple aspects of DCLK signaling output. Notably, a conserved segment of the C-tail functions as an isoform-specific autoinhibitory motif, which mimics ATP functions through direct tail docking to the nucleotide-binding pocket, where it forms an ordered set of interactions that aligns the catalytic (C) spine of the kinase in the absence of ATP binding. Furthermore, molecular modeling demonstrates that a phosphorylated threonine in the C-tail of DCLK1.2, which is absent in DCLK1.1, is positionally poised to competitively mimic the gamma phosphate of ATP, perhaps in a regulated manner. Other segments of the tail function as a pseudosubstrate by occluding the substrate-binding pocket and tethering to key functional regions of the catalytic domain. Thermostability analysis of purified DCLK1 proteins, combined with molecular dynamics (MD) simulations, confirms major differences in thermal and dynamic profiles of the DCLK1 isoforms, while catalytic activity assays reveal how specific variations in the G-loop and C-tail can rescue DCLK1.2 from the autoinhibited conformation. Together, these studies demonstrate that isoform-specific variations in the C-terminal tail co-evolved with residues in the DCLK kinase domain, contributing to regulatory diversification and functional specialization.
Results
Origin and evolutionary divergence of DCLK family members
The human DCLKs’ repertoire is composed of three genes, termed DCLK1, -2, and -3 (Figure 1A, Table 1). The experimental model employed in this study, DCLK1, is composed of multiple spliced variants in human cells. Those full-length proteins that contain N-terminal DCX domains are usually referred to as DCLK1.1 or DCLK1.2 and the variants that lack the DCX domains are termed here (for simplicity) ΔDCLK 1.1 and ΔDCLK1.2 (also referred to as DCLK1.3 and DCLK1.4). The core catalytic domain with minimal flanking regions (DCLK1cat, Figure 1B) is identical in all DCLK1 proteins, whereas the length of the tail, or the presence of the DCX domains, generates considerable diversity from the single human DCLK1 gene (Figure 1B, Figure 1—figure supplement 1). To infer evolutionary relationships of DCLK paralogs, and especially the evolution of the C-terminal tail regions that lie adjacent to the kinase domain (Figure 1), we performed phylogenetic analysis of 36 DCLK sequences with an outgroup of closely related CAMK sequences (Figure 2A, Figure 2—source data 1). These DCLK sequences are from a representative group of holozoans, which consist of multicellular eukaryotes (metazoans) and closely related unicellular eukaryotes (pre-metazoans). The analysis generated four distinct clades: pre-metazoan DCLK, metazoan DCLK3, vertebrate DCLK2, and vertebrate DCLK1. Interestingly, DCLK genes demonstrated significant expansion and diversification within metazoan taxa. The pre-metazoan DCLK sequences were the most ancestral and showed no DCLK diversity, suggesting the DCLK expansion and diversification correlated with the evolution of multicellular organisms. Within the metazoan expansion of DCLK, DCLK3 is the most ancestral and can be broken down into two sub-clades: protostome DCLK3 and deuterostome DCLK3. Within invertebrates, only two DCLK paralogs were present, one that was identified as a DCLK3 ortholog and another that was not clearly defined as either DCLK1 or DCLK2. This suggests that the diversification into DCLK1 and DCLK2 paralogs from an ancestral DCLK1/2-like paralog occurred after the divergence of invertebrates and vertebrates, which is further supported by the monophyletic DCLK1 and DCLK2 clades in vertebrates (bootstrap value: 99).

Evolution of the doublecortin-like kinase (DCLK) family.
(A) Phylogenetic tree showing the divergence and grouping of DCLK sub-families in different taxonomic groups. Bootstrap values are provided for each clade. (B) shows domain annotations for sequences included in the phylogenetic tree. The length of C-terminal tail segment for these sequences is shown as a histogram (green). The original tree generated using IQTREE is provided in Figure 2—source data 1.
-
Figure 2—source data 1
Full doublecortin-like kinase (DCLK) phylogenetic tree with bootstrap values.
- https://cdn.elifesciences.org/articles/87958/elife-87958-fig2-data1-v1.pdf
-
Figure 2—source data 2
Species and Uniprot IDs for each sequence represented in the domain plot of Figure 2.
- https://cdn.elifesciences.org/articles/87958/elife-87958-fig2-data2-v1.xlsx
Interestingly, the expansion of DCLK in metazoans and the diversification of DCLK1 and DCLK2 within vertebrates correlates well with the length and sequence similarity of the C-terminal tail, which also varies between the different DCLK1 splice variants (Figure 1A, Figure 1—figure supplement 1 and Figure 2—figure supplement 1). Within both protostome and deuterostome DCLK3, the length of the C-terminal tail is ~50 residues or less. This is in marked contrast to the tail lengths of vertebrate DCLK1 and DCLK2, which are ~100 residues long. In addition to the C-terminal tail, an analysis of the domain organization of these DCLKs reveals that DCLK3 predominantly contains only a single N-terminal doublecortin domain (DCX), whereas invertebrate DCLKs, and vertebrate DCLK1 and DCLK2 predominantly contain two DCX domains at the N-terminus of the long isoforms (Figure 2B). In addition, we identified a putative active site-binding motif, VSVI, and a phosphorylatable threonine conserved within vertebrate DCLK1 and DCLK2, which is absent in all other DCLK sequences, including invertebrate DCLK1/2. This raises the possibility that the DCLK1/2 tail extensions are employed for vertebrate-specific regulatory functions.
Next, we compared the type of DCLK1 protein sequence encoded by a range of chordate mammalian genomes. The domain organization of each DCLK1 isoform was compared based on annotated sequences from UniProt, demonstrating the presence of at least one DCLK1 protein that lacks the DCX domains in every species examined, with a mixture of ΔDCLK1.1 and ΔDCLK1.2 splice variants. Interestingly, it was only in the human DCLK1 gene that definitive evidence for ΔDCLK1.1 and ΔDCLK1.2 variants were found (Figure 3A). To establish a model for DCLK1 biophysical analysis, we constructed a recombinant hybrid human DCLK1 catalytic domain with a short C-tail sequence that is equivalent to DCLK1.1 amino acids 351–689, containing the catalytic domain with a short C-tail region. As shown in Figure 3B, incubation of size exclusion chromatography (SEC)-purified GST-tagged DCLK1 with 3C protease generated the mature untagged DCLK1 protein for biophysical analysis. Analytical SEC revealed that purified DCLK1.1 and DCLK1.2 isoforms are monomeric in solution (Figure 3—figure supplements 1–3). We evaluated catalytic activity for DCLK1.1351-689 using a validated peptide phosphorylation assay (Figure 3C), which revealed efficient phosphorylation of a DCLK1 substrate peptide. The KM[ATP] for peptide phosphorylation was close to 20 µM in the presence of Mg2+ ions (Figure 3C, left panel), similar to values measured for other Ser/Thr kinases that are autophosphorylated and active after expression from bacteria (Byrne et al., 2020b). DCLK-dependent peptide phosphorylation was completely blocked (Figure 3C, right panel) by prior incubation of the reaction mixture (containing 1 mM ATP) with the chemical inhibitor DCLK1-IN-1, as expected (Ferguson et al., 2020). In addition to enzyme activity, we monitored thermal denaturation of purified, folded, DCLK1351-689 protein in the presence of ATP, either alone or as a Mg:ATP complex, which is required for catalysis. As shown in Figure 3D, DCLK1 was stabilized by 2.1°C upon incubation with an excess of Mg:ATP, and this protective effect was completely blocked by mutation of Asp 533 (of the conserved DFG motif) to Ala, consistent with canonical ATP interaction in the nucleotide-binding site. Finally, we assessed the thermal effects of a panel of DCLK1 inhibitors on the model DCLK1.1351-689 protein. Prior incubation with DCLK1-IN-1, LRRK2-IN-1, the benzopyrimidodiazipinones XMD8-92 and XMD8-85, which have been reported to potently (though not specifically) inhibit DCLK1 activity (Patel et al., 2021), led to marked protection from thermal unfolding (Figure 3—figure supplement 4). Consistently, the negative control compound DCLK1-Neg (Ferguson et al., 2020) was ineffective in stabilizing DCLK1.

Cartoon cladogram of mammalian DCLK1 isoforms and experimental.
(A) Cartoon cladogram of mammalian species showing the domain organization of each DCLK1 isoform from representative annotated sequences from UniProt. UniProt IDs for each sequence are provided in Figure 3—source data 1. (B) SDS-PAGE of 6His-GST-3C-DCLK1.1 (351–689, top) or a D533A mutant in which the DFG Asp is mutated to Ala (middle). Proteins were separated by size exclusion chromatography, and high-purity fractions were pooled. The affinity tag was removed prior to analysis by incubation with 3C protease, leading to a demonstrable shift in mobility (bottom). (C) Evaluation of catalytic activity toward DCLK1 peptide. DCLK1.1 351–689 possesses a KM[ATP]~20 µM in vitro (left) and real-time substrate phosphorylation was inhibited by prior incubation with the small molecule DCLK1-IN-1 (right). (D) Thermal shift assay demonstrating a 2.1°C increase in the stability of DCLK1 351–689 in the presence of Mg:ATP (left), which was absent in the D533A protein (right). Raw data are provided in Figure 3—source data 2.
-
Figure 3—source data 1
Species and Uniprot IDs for each sequence represented in the mammalian DCLK1 isoform species tree.
- https://cdn.elifesciences.org/articles/87958/elife-87958-fig3-data1-v1.xlsx
-
Figure 3—source data 2
Experimental values for thermal shift and activity assays from Figure 3.
- https://cdn.elifesciences.org/articles/87958/elife-87958-fig3-data2-v1.xlsx
-
Figure 3—source data 3
DCLK1 plasmid sequences generated for Figure 3.
- https://cdn.elifesciences.org/articles/87958/elife-87958-fig3-data3-v1.txt
-
Figure 3—source data 4
Figure 3B (top) uncropped gel.
- https://cdn.elifesciences.org/articles/87958/elife-87958-fig3-data4-v1.zip
-
Figure 3—source data 5
Figure 3B (bottom) uncropped gel.
- https://cdn.elifesciences.org/articles/87958/elife-87958-fig3-data5-v1.zip
-
Figure 3—source data 6
Uncropped gel of incubation of size exclusion chromatography (SEC)-purified DCLK1 with 3C protease.
- https://cdn.elifesciences.org/articles/87958/elife-87958-fig3-data6-v1.zip
-
Figure 3—source data 7
PowerPoint file containing labels for Source Data 4–6.
- https://cdn.elifesciences.org/articles/87958/elife-87958-fig3-data7-v1.pptx
Key differences between isoforms in the C-tail of DCLK1 arise from alternative-splicing and different open-reading frames
Higher-order vertebrates have multiple isoforms of DCLK1 and DCLK2, where sequence variations occur in either or both the N- and C-terminal regions attached to the kinase domain. Human DCLK1, for example, has four unique isoforms. Isoforms 1 and 2 differ in C-terminal tail length due to variations in exon splicing (Figure 4A). Further examination of the intron and exon boundaries indicates that human DCLK1.1 contains an additional exon (exon 16) that is not spliced in DCLK1.2. Exon 16 is spliced with exon 17 with a phase 2 intron, which introduces a shift in the reading frame and an earlier translated stop codon (UGA) in exon 17 (Figure 4B). In DCLK1.2, exon 15 is spliced with exon 17, with a non-disruptive phase 0 intron, resulting in the full translation of exon 17. These changes introduce multiple indels (insertions and deletions) and result in the insertion of a phosphorylatable threonine (T688) (Agulto et al., 2021) in DCLK1.2 that is absent in the DCLK1.1 variant (Figure 4B–C), suggesting a possible exon duplication for adaptive regulation of DCLK1 function by phosphorylation. DCLK1.2 is the best-characterized isoform in terms of structure and function (Cheng et al., 2022b), and to compare it with DCLK1.1, we generated a series of C-terminal tail deletion mutants to evaluate how variations in the C-terminal tail contribute to isoform-specific DCLK1 functions.

Transcriptional Variations in the C-terminal Tail of Human DCLK1 Isoforms: Gene Organization, Exon Alternation, and Protein Sequence Alignment.
(A) Gene and intron-exon organization of DCLK1 human isoforms in the C-terminal tail. The DCLK1 gene is present on locus 13q13.3, and isoforms 1 and 3 contain an additional exon (exon 16) in the C-terminal tail that is absent in DCLK1.2. (B) A phase 2 intron results in the alternative transcript of exon 17 in isoform 1, translating a different open-reading frame and early stop codon, resulting in the shorter sequence. (C) Cartoon organization of the C-tail exons (exon 15, 16, and 17) of the DCLK1 isoforms, comparing the translated protein sequence alignment.
-
Figure 4—source data 1
Scipio output for intron-exon analysis of DCLK1.1.
- https://cdn.elifesciences.org/articles/87958/elife-87958-fig4-data1-v1.zip
-
Figure 4—source data 2
Scipio output for intron-exon analysis of DCLK1.2.
- https://cdn.elifesciences.org/articles/87958/elife-87958-fig4-data2-v1.zip
-
Figure 4—source data 3
Scipio output for intron-exon analysis of DCLK1.3.
- https://cdn.elifesciences.org/articles/87958/elife-87958-fig4-data3-v1.zip
-
Figure 4—source data 4
Scipio output for intron-exon analysis of DCLK1.4.
- https://cdn.elifesciences.org/articles/87958/elife-87958-fig4-data4-v1.zip
Isoform-specific variations encode changes in MD, thermostability, and catalytic activity in DCLK1
To study isoform-specific differences in the C-tail, we employed experimental techniques to compare protein stability and catalytic activity between purified DCLK1 proteins alongside MD simulations for DCLK1.1 and DCLK1.2 with different tail lengths. Isoforms 1.1 and 1.2 share identical sequences across the kinase domain and within the first 38 residues of the C-tail, and we used this information to design a new recombinant protein, termed DCLK1cat (residues 351–686). C-terminal to this totally conserved region, both isoforms possess extended tail segments, which includes the putative inhibitory-binding segment (IBS; residues 682–688) and an additional intrinsically disordered segment (IDS; residues 703-end). To study the role of the C-tail in modulating kinase stability and activity, we purified the DCLK1cat, and C-tail containing (long and short) variants of each isoform, each of which lack the N-terminal DCX domains (Figure 5A). SDS-PAGE demonstrated that protein preparations were essentially homogenous after affinity and gel filtration chromatography (Figure 5—figure supplement 1). The short forms of the recombinant proteins (DCLK1.1351-703 and DCLK1.2351-703) possess a partially truncated C-tail and were designed to match the amino acid sequence previously used to solve the structure of DCLK1.2 protein (Cheng et al., 2022b). Notably, these proteins exclude the IDS. The long forms of the DCLK1 proteins include the full-length C-tail for each isoform (DCLK1.1351-729 and DCLK1.2351-740) and incorporate IDS domains. We first performed comparative thermal shift analyses to quantify variance in thermal stability between the different purified proteins. When contrasting DCLK1.1short and DCLK1.2short which do not differ in size or tail length but encode unique sets of amino acids in their partially truncated C-tail as a result of alternative splicing (Figure 4), we observed that DCLK1.2 was some 14°C more stable than DCLK1.1 (Figure 5B). When compared with DCLK1cat, both DCLK1.1 short and long exhibited only subtle changes in thermal stability (Figure 5C and E), whereas both DCLK1.2 proteins (DCLK1.2short and DCLK1.2long) were significantly stabilized (relative ΔTm >16°C, Figure 5D and E).

Structural Analysis, Thermal Stability, and Molecular Dynamics of DCLK1 Variants: Insights into the IBS and IDS Domains.
(A) Cartoons of DCLK1 construct used in our assays, portraying the locations of the inhibitory-binding segment (IBS) and the intrinsically disordered segment (IDS). (B–E) Differential scanning fluorimetry (DSF) thermal denaturation profiling of the purified DCLK1 core catalytic domain, or tail-matched DCLK1.1 and DCLK1.2 proteins. Unfolding curves and changes in Tm values (ΔTm) for each protein relative to WT DCLK1cat are indicated. (F–H) B-factor structural representations of DCLK1short proteins shown in (A). The width of the region indicates the extent of flexibility based on averaged root mean square fluctuation (RMSF) data from three 1 µs molecular dynamics (MD) replicates. (I) DSSP analysis of three replicates of 1 µs MD simulations showing the residues surrounding the IBS in the C-tail of DCLK1.1short and DCLK1.2short. Blue indicates the presence of a beta-sheet or beta-bridge secondary structures and red indicates the presence of alpha-helical structures.
-
Figure 5—source data 1
Root mean square fluctuation (RMSF) values for each mutant from all replicates of molecular dynamics (MD) simulations.
- https://cdn.elifesciences.org/articles/87958/elife-87958-fig5-data1-v1.xlsx
We next performed MD simulations to study the dynamics within the distinct DCLK1 C-tails that might explain the observed difference in protein stability. The crystal structure of DCLK1.2 (PDB: 6KYQ) was employed for the DCLK1.2short model and AlphaFold2 was used to model the other proteins (DCLK1cat and DCLK1.1short). Comparison of the root mean square fluctuations (RMSF) of the two isoforms in three different replicates of MD simulations indicates strikingly different thermal fluctuations in the C-terminal tails and catalytic domains (Figure 5—source data 1). In particular, the IBS (between 682 and 688) is stably docked in the ATP-binding pocket in DCLK1.2 and an alpha-helical conformation is maintained during the microsecond time scale across different replicates (Figure 5I, bottom). In contrast, the IBS is more unstable in DCLK1.1, as indicated by high thermal fluctuations and a lack of secondary structure propensity (Figure 5I, top). A caveat to bear in mind is that DCLK1.1 is an AlphaFold2 model, which will also account for increased RMSF. Analysis of sequence variations and structural interactions provides additional insights into the differential dynamics of the two isoforms. The helical conformation of the IDS in DCK1.2 is maintained during the simulation due, in part, to a capping interaction with Thr 687, which is absent in DCLK1.1 due to the alternative splicing event detailed above. Likewise, another key residue in DCLK1.2, Lys 692, anchors the tail to the catalytic domain through directional salt bridges with the conserved aspartates (Asp 511 and Asp 533) in the HRD and DFG motifs (Figure 5, Figure 5—figure supplement 2A). These interactions are not observed in DCLK1.1 simulations because Lys 692 is substituted to a histidine (His 689), which is unable to form a corresponding interaction with the catalytic domain (Figure 5, Figure 5—figure supplement 2B). We also evaluated the effects of T688A (non-phosphorylated) or T688E (phosphomimetic) mutations through DCLK1 MD simulations and found that the two mutations slightly destabilize the tail relative to WT. Three replicates of the two mutants show increased RMSF of the tail region relative to WT DCLK1.1 (Figure 5, Figure 5—figure supplement 3). Either mutation was not sufficiently destabilizing on its own to unlatch the C-tail, and we hypothesize that other residues in addition to T688 are also likely to be important for contributing to conformational regulation of the kinase domain by the C-terminal tail. The variable docking of the C-tail within the kinase domains of the two DCLK1 isoforms, and the extent to which this contributes to more transient or stable autoinhibited states are explored in more detail in the next section.
Residues contributing to the co-evolution and unique tethering of the C-terminal tail to the DCLK catalytic domain
To identify specific residues that contribute to the unique modes of DCLK regulation by the C-terminal tail, we performed statistical analysis of the evolutionary constraints acting on DCLK and related CAMK family sequences. We aligned the catalytic domain of DCLK and related CAMK sequences from diverse organisms and employed the Bayesian partitioning with pattern selection (BPPS) method (Neuwald, 2014) to identify residues that most distinguish DCLK sequences (foreground alignment in Figure 6B) from CAMK sequences (background alignment). Beyond the catalytic domain, DCLKs share sequence and structural similarities in the first helix of the tail (αR1 in CAMK1) (Rosenberg et al., 2005; Goldberg et al., 1996), with other CAMKs but share no detectable sequence similarity beyond this helical segment. DCLKs also share a CAMK-specific insert segment located between F and G helices in the catalytic domain, although the nature of residues conserved within the insert is unique to individual CAMK families (Figure 6—figure supplement 1). BPPS analysis revealed DCLK-specific constraints in different regions of the kinase domain, most notably, the ATP-binding G-loop, N terminus of the C-helix, the activation loop, and C-terminus of the F-helix (Figure 6—figure supplement 2).

Identification of doublecortin-like kinase (DCLK)-specific constraints.
(A) Cartoon of DCLK1.2 and the intrinsically disordered segment (IDS) with evolutionary constraints mapped to the kinase domain and C-tail. (B) Sequence constraints that distinguish DCLK1/2/3 sequences from closely related calcium calmodulin kinase (CAMK) sequences are shown in a contrast hierarchical alignment (CHA). The CHA shows DCLK1/2/3 sequences from diverse organisms as the display alignment. The foreground consists of DCLK sequences while the background alignment contains related CAMK sequences. The foreground and background alignments are shown as residue frequencies below the display alignment in integer tenths (Manning et al., 2002; Agulto et al., 2021; Bayer and Schulman, 2019; Gógl et al., 2019; Berginski et al., 2021; Sossey-Alaoui and Srivastava, 1999; Ohmae et al., 2006; Couillard-Despres et al., 2005; Horesh et al., 1999). The histogram (red) indicates the extent to which distinguishing residues in the foreground diverge from the corresponding position in the background alignment. Black dots indicate the alignment positions used by the Bayesian partitioning with pattern selection (BPPS) (Neuwald, 2014) procedure when classifying DCLK sequences from related CAMK sequences. Alignment number is based on the human DCLK1.2 sequence (UniProt ID: O15075-2). (C) Sequence alignment of human DCLK1 isoforms.
Some of the most significant DCLK-specific residue constraints map to the ATP-binding G-loop (GDGNFA motif) (Figure 6A–B). In particular, Asp 398, Asn 400, and Ala 402 are unique to DCLKs as the corresponding residues are strikingly different in other CAMKs. Asp 398 is typically a charged residue (K/R) in other CAMKs while Asn 400 and Ala 402 are typically hydrophobic and polar residues, respectively (see residue frequencies in background alignment; Figure 6B). Notably, both Asn 400 and Asp 398 make direct interactions with residues in the C-tail either in the crystal structure or MD simulations (see below). Likewise, DCLK conserved residues in the C-helix and activation loop tether the C-terminal tail to functional regions of the kinase core, suggesting co-option of the DCLK catalytic domain to uniquely interact with the flanking cis-regulatory tail.
An autoinhibitory ATP-mimic completes the C-spine and mimics the gamma phosphate of ATP
The most stable segment of the C-tail based on the B-factor and RMSF fluctuations in MD simulations is a unique region (682–688 in DCLK1.2) that docks into the ATP-binding pocket through both hydrophobic and hydrogen-bonding interactions. Remarkably, this peptide segment mimics the physiological ATP ligand, and stabilizes the catalytic domain through ‘completion’ of the hydrophobic catalytic spine (Figure 7A; Zheng et al., 1993). The residues that mimic adenosine and complete the C-spine of DCLK1 are Val 682, Val 684, and Ile 685 (PDB: 6KYQ), which are part of the alpha-helix that docks to the ATP-binding pocket (Figure 7B). Interestingly, based on our BPPS analyses, these C-tail residues are uniquely vertebrate DCLK1-specific pattern constraints. At the tail end of this alpha-helix are two Thr residues, Thr 687 and Thr 688. As previously noted, these Thr residues mark the beginning of exon 17, and are one of the key variations between human DCLK1 isoforms, found only in DCLK1.2 variants.

The DCLK1 C-tail 'completes’ the regulatory C-spine (green).
(A) PKA crystal structure (PDB: 1ATP) with bound ATP in red and Mg2+ in purple. The C-spine is completed by the adenine ring of ATP. The gamma phosphate of ATP hydrogen bonds with the second glycine of the G-loop. (B) DCLK1.2 crystal structure (PDB: 6KYQ) showing how the C-tail (red) docks underneath the pocket and mimics the ATP structure. The C-spine is completed by V682 and V684 in the C-tail and helical segments defined using DSSP are shown. T687 is also depicted making multiple hydrogen bonds with the backbone of V684 and I685 (dashed lines). (C) DCLK1.1 AlphaFold2 model showing an unstructured loop in the C-tail docking into the ATP-binding pocket, where V684 and I685 are predicted to complete the C-spine. The average per-residue confidence of the C-tail is 49%. (D–F) Zoomed-out versions of A–C, demonstrating how the DCLK1 C-tail docks into the ATP-binding cleft, akin to ATP in PKA.
Structural analysis and MD simulations reveal that Thr 687 in DCLK1.2 caps the stable alpha-helix that extends the C-spine (Figure 7B, Figure 7—figure supplement 1B). In comparison, the same region in DCLK1.1, which lacks Thr 687, is predicted to be unstructured. Upon phosphorylation, Thr 688 in DCLK1.2 can mimic the gamma phosphate of ATP by maintaining a stable hydrogen bonding distance with the backbone of the second glycine of the G-loop (G399) (Figure 7B, Figure 7—figure supplement 1C). We additionally observe that the sidechain of Asn 400, a DCLK-specific G-loop constraint, further stabilizes the phosphate group in pThr 688 through hydrogen bonding. As previously described, Thr 688 is unique to DCLK1.2. The lack of this functional site in DCK1.1 is correlated with increased RMSF and instability of the ATP-mimic segment in isoform 1 MDs (Figure 5C–D, Figure 7—figure supplement 1B). Comparatively, MD analysis of DCLK1.2 and a phosphothreonine-containing DCLK1.2 demonstrates reduced C-tail fluctuations, suggesting the potential regulatory involvement of Thr 688 phosphorylation for further modulation of the autoinhibited conformation (Figure 7—figure supplement 1C), consistent with previous findings (Agulto et al., 2021).
Mutational analysis support isoform-specific allosteric control of catalytic activity by the C-terminal tail
To evaluate how sequence differences between DCLK1.1 and DCLK1.2 affected both thermal stability and catalytic potential, we generated targeted mutations at contact residues within the Gly-rich loop and C-tail of DCK1.1 and DCLK1.2 (at the indicted residues depicted in Figure 8A) which we predicted would disrupt or destabilize C-tail docking within the domain. All proteins were purified to near homogeneity by IMAC and SEC (Figure 5—figure supplement 1, Figure 5—figure supplement 1—source data 1), and the thermal stability of a panel of DCLK1.2 mutant and WT proteins were compared side-by-side with the DCLK1cat (Figure 8B). The ΔTm values obtained (Figure 8C) demonstrate that mutation of Asp 398 or Asn 400 in the Gly-rich loop are by themselves insufficient to destabilize DCLK1.2. In marked contrast, dual mutation of the hydrophobic pair of Val 682 and Val 684 residues to Thr, or mutation of the acidic tail residue Asp 691, resulted in a pronounced reduction in DCLK1.2 thermal stability. Moreover, the recorded Tm values for these latter two mutations quite closely resembled the Tm of DCLK1cat (which lacks the C-tail entirely), which is consistent with the uncoupling of the C-tail and a commensurate decrease in thermal stability associated with loss of this interaction. We next determined the catalytic activity of our recombinant DCLK1.1 and DCLK1.2 proteins side-by-side (Figure 8D, Figure 8—figure supplement 1). Although partially diminished in relation to DCLK1cat, both DCLK1.1short (351–703) and DCLK1.1long (351–729) variants possess robust catalytic activity. This suggested ineffective ATP-competitive autoinhibition mediated by the C-tail segment of DCLK1.1 and is consistent with their closely matched Tm values to DCLK1cat (Figure 5C). Interestingly, both C-tail containing variants of DCLK1.1 (and particularly DCLK1.1351-729) exhibited lower affinity for ATP (inferred from KM[ATP] for peptide phosphorylation), which is consistent with partial occlusion of the ATP-binding pocket (Figure 8—figure supplement 1). In marked contrast, the detectable kinase activity for short (351–703) or long (351–740) DCLK1.2 proteins was significantly blunted compared to DCLK1cat, exhibiting just ~5% of the activity of the catalytic domain alone, and consistently, the calculated KM[ATP] was ~4-fold higher compared to the catalytic domain lacking the C-tail. We also utilized autophosphorylation as a proxy for overall kinase activity. Quantitative tandem mass spectrometry (MS/MS) analysis of site-specific autophosphorylation within DCLK1.1short and DCLK1.2short demonstrate a marked reduction in the site-specific abundance of phosphate in DCLK1.2 when compared to DCLK1.1 at two separate sites that could be directly and accurately quantified by MS (S438 and S660, DCLK1.1 relative abundance set to 1, Figure 8—figure supplement 2). LC-MS/MS also indicated that several autophosphorylation sites identified in isoform 1 were absent in DCLK1.2 (Ser 683 and Thr 692, the latter of which is an amino acid that is unique to the C-tail of DCLK1.1, Figure 8—figure supplement 2). Interestingly, amino acid substitutions in the G-loop or the C-tail of DCLK1.2 designed to subvert C-tail and ATP site interactions also had major effects on DCLK1.2 phosphorylation and catalytic activity. DCLK1.2 D398A was activated some fivefold when compared to the WT form, whereas DCLK1.2 N400A was almost as active as the DCLK1cat. Consistently, DCLK1.2 V682T/V684T and D691A proteins were also much more active than the WT form of DCLK1.2. Kinetic analysis also confirmed higher Vmax (but broadly similar KM[ATP]) values for DCLK1.2 D398A and V682T/V684T relative to the WT protein (Figure 8—figure supplement 1). Moreover, comprehensive LC-MS/MS phosphosite mapping revealed a marked increase in the total number of phosphorylated amino acids in all of the mutant DCLK1.2 proteins, consistent with the enhanced catalytic activity of these proteins when compared to WT DCLK1.2 (Figure 8—figure supplement 2). Together, these observations confirm that targeted mutations are sufficient to relieve ATP-competitive C-tail autoinhibition by physical tail uncoupling; this model is also strongly supported by marked changes in the biophysical stability of the mutant proteins, particularly for V682/684T and D691A (Figure 8B).

Structural and Functional Analysis of DCLK1.2 Variants: Insights into Amino Acid Modifications, Thermal Stability, and Kinase Activity.
(A) Structural depiction of DCLK1.2 (PDB: 6KYQ) showing the location of modified DCLK1 amino acids on the G-loop (purple) or C-tail (red). (B–C) Differential scanning fluorimetry assays depicting thermal denaturation profiles of each protein along with the calculated Tm value. (D) Kinase assays. DCLK1-dependent phosphate incorporation (pmol/min) into the DCLK1 peptide substrate was calculated for DCLK1cat, long and short DCLK1.1, and the indicated DCLK1.2 variants. (E) Thermal stability analysis in the presence of ATP or DCLK1-IN-1 for DCLK1 proteins. For DCLK1.2, all proteins were generated in the DCLK1.2 short background.
To expand on this finding, we investigated the interaction of Mg:ATP or DCLK1-IN-1 to our panel of DCLK1.1 and -1.2 proteins (Figure 8E), using changes in thermal stability as a reporter of ligand binding. DCLK1cat (351–686), DCLK1.1short (351–703) and DCLK1.1long (351–729) proteins all behaved similarly in the presence of either Mg:ATP or DCLK1-IN-1, inducing marked stabilization. In contrast, DCLK1.2short (351–703) or DCLK1.2long (351–740) proteins registered negligible thermal shifts in the presence of the same concentration of either ligand, which is in line with the C-tail tightly occupying the ATP-binding site and obstructing their binding. Remarkably, D398A and N400A DCLK1.2, whose high basal Tm values (compared to DCLK1cat) are consistent with stabilization by docking of the C-tail within the kinase domain, were markedly destabilized in the presence of either ligand. This suggests that incorporation of these G-loop mutations in isolation is insufficient to dislodge the C-tail, but rather that the stability of the interaction is compromised to the extent that either ATP or DCLK-IN-1 can competitively dislodge the bound tail from the ATP active site, resulting in a net destabilization caused by lack of tail engagement. This observation is corroborated by the results of our kinase assays, where both D398A and N400A mutants were more active than WT DCLK1.2, confirming appropriate ATP binding (a prerequisite for catalysis). Finally, DCLK1.2 V682T/V684T and D691A, which exhibit lower basal Tm values than the WT protein (indicating a loss of tail interaction), were both stabilized in a ligand-dependent manner to a similar degree to that observed for DCLK1cat and DCLK1.1 proteins (Figure 8E). Collectively, these observations clearly demonstrate that the C-tail section of DCLK1.2 can both stabilize the canonical DCLK kinase domain and inhibit kinase activity (by impeding the binding and structural coordination of Mg:ATP) much more effectively than that of DCLK1.1. Our targeted mutational analysis of key contact residues in DCLK1.2 also clearly shows that this is a consequence of specific amino acid interactions that are absent in DCLK1.1 due to alternative-splicing and subsequent sequence variation.
Classification of DCLK regulatory segments
We synthesized our experimental findings by classifying the DCLK C-tail into six functional segments, based on interactions with different regions of the catalytic domain and their conservation between the C-tail splice variants in our analysis (Figure 9): First is an ATP-mimetic peptide segment (residues 682–688 in DCLK1.2) that readily mimics physiological ATP binding by completing the C-spine in the nucleotide-binding site. The inhibitory peptide also contains a phosphorylatable Thr residue, which sits adjacent to the highly characteristic Gly-rich loop (GDGNFA, residues 396–402). Second, we define a pseudosubstrate mimic (residues 692–701), which interacts with the acidic HRD and DFG Asp sidechains and docks in the substrate pocket occluding substrate access. Third, at the C-terminus of the tail lies an IDS (residues 702–749, Figure 9—figure supplement 1), which packs dynamically against DCLK conserved residues in the kinase activation loop. Fourth, at the beginning of the C-tail lies a CAMK-tether (residues 654–664), a set of residues that pack against a CAMK-specific insert in the C-lobe. In many CAMK crystal structures, this insert makes multiple contacts with the F-helix and C-tail (Figure 6—figure supplement 1). Fifth, this is followed by a highly dynamic pseudosubstrate region (residues 672–678) that occludes the substrate pocket and will thus interfere with substrate phosphorylation. Sixth, a transient beta-strand is formed in DCLK1.2 through amino acid-specific sequences that help modulate and potentially strengthen binding of the C-tail in this isoform (Figure 7—figure supplement 1A–B). Collectively, these segments and their associated interaction sites demonstrate that co-evolution of the unique C-tail with the catalytic domain is the central hallmark of DCLK functional divergence, and that changes in these segments possess the ability to ‘supercharge’ catalytic output of the kinase. In particular, the variable C-terminal segments of the tail might contribute to isoform-specific functional specialization. The combinatorial diversity of events that modulate C-tail function may allow DCLKs to nimbly coordinate various tasks including ATP-binding, substrate-based phosphorylation, regulation of DCX domain phosphorylation and structural disposition, kinase autoinhibition, and allosteric regulation. Isoform-specific variability provides additional nuance to regulatory and catalytic signaling events and may even contribute to differences in cellular localization (e.g. cytoplasm or nucleoplasm) and tissue-specific activity, enabling contextual DCLK regulation through these modular sequence segments.

A DCLK1 C-tail can act as a multi-functional Swiss Army Knife, using six distinct segments for a variety of regulatory functions including mimicking ATP binding/association, stabilizing the G-loop, occluding the substrate-binding pocket, and packing against the kinase activation loop.
Discussion
The kinase domain is a conserved switch for phosphorylation-based catalytic regulation. Yet the complexity of cell signaling pathways demands other nuanced forms of regulation beyond binary ‘on’ or ‘off’ switch-based mechanisms. For many Ser/Thr kinases, including AGC and CAMK families, these distinct regulatory functions come from segments which flank the kinase domain, N- and C-terminal regions, which serve to modulate activity through allosteric activation, inhibition, or rheostatic behaviors that change based on environmental conditions (Gógl et al., 2019). In this study, we expand on our knowledge of allosteric diversity in the human kinome by revealing how alternative splicing of the DCLK1 C-tail contributes to isoform-specific behaviors, coupling regulation of catalytic output, phosphorylation, protein dynamics and stability, substrate binding, and protein-protein interactions. Our ‘Swiss Army Knife’ model for DCLK1 expands our view of allosteric regulation as not just a dynamic process facilitated by proteins, but one where adaptive genetic mechanisms, like differential splicing, dexterously tune isoform-specific functions for specific cellular signaling roles; in the case of DCLK1.1, this allows ‘supercharging’ of catalysis between splice variants due to key amino acid differences in the C-tail that are lacking in the DCLK1.2 isoform.
Multiple members of the human kinome have independently evolved C-tail regions that dock to the N- or C-lobe of the kinase domain in cis. In the case of the AGC kinases, the C-terminal tail is a very well-studied in cis modulatory element that serves to explain a variety of regulatory properties in this kinase sub-family (Kannan et al., 2007; Romano et al., 2009; Baffi and Newton, 2022; Taylor and Kornev, 2011). Classical deletion studies with members of the CAMK family have also revealed a cis-acting inhibitory element lying C-terminal to the catalytic domain of both CAMK1 (Yokokura et al., 1995) and CAMK2 (Yang and Schulman, 1999). More recent examples of C-tail functional diversity in CAMK family members are presented by the human pseudokinases TRIB1 and TRIB2, which employ C-tail sequences to either latch (and structurally restrict) the atypical kinase domain or to bind competitively to the ubiquitin E3 ligase COP1 (Eyers, 2015; Eyers et al., 2017). Functional disengagement of the TRIB1 or TRIB2 tail through deletion, mutagenesis, or small molecule binding has marked effects on pseudokinase conformation, intrinsic stability, and cellular transformation (Foulkes et al., 2018; Harris et al., 2022; Keeshan et al., 2010; Murphy et al., 2015).
A novel pseudosubstrate region encoded by DCLK1
In addition to the marked differences between DCLK1 splice variants relevant to nucleotide binding, small molecule interactions, and catalysis, our work also reveals two unique pseudosubstrate segments present before and after the IBS. Before the IBS, we observe the formation of an anti-parallel transient beta-sheet with the beta1 strand in the catalytic domain (Figure 5I, Figure 7—figure supplement 1A–B). During the formation of this transient structure, the C-tail dynamically occludes part of the substrate-binding pocket. A beta-sheet is observed in all three MD replicates of DCLK1.1, but only in a single replicate of DCLK1.2 dynamic analysis. At the other end of the IBS is another pseudosubstrate segment whose structure and dynamics change as a result of alternative exon splicing. In DCLK1.2, the pseudosubstrate segment is stable, with an average RMSF of 1.8 Å, facilitated by key interactions from Lys 692, which coordinates acidic residues in the HRD and DFG motifs (Figure 5—figure supplement 2A). In DCLK1.1, Lys 692 is replaced by a His, which weakly coordinates with the HRD and DFG motifs, resulting in increased dynamics of the segment and an RMSF of 3.1 Å (Figure 5, Figure 5—figure supplement 2B). Together, residue variation between isoforms contributes to differences in stability and alteration of dynamics of the tail. HPCAL1 was recently proposed as a possible substrate that activates DCLK1 in a calcium-dependent manner, but it is unclear how it may bind DCLK1 (Patel et al., 2021). Because only exon 15 of the C-tail is conserved between the isoforms, it is possible that the location of binding occurs in this dynamic pseudosubstrate segment prior to the IBS, where increased flexibility and occlusion of the substrate pocket is reflective of the absence of HPCAL1, or a similar CaM-like substrate.
Discovery of the DCLK1 ATP-mimic region, a splice-variant specific regulatory module
Our structural analysis of DCLK1 reveals a remarkable structural mimic of ATP located in the C-tail, which differs markedly between DCLK1.1 and DLCK1.2 splice variants. We note for the first time that a set of three residues in the ATP-mimic, Val 682, Val 684, and Ile 685, are conserved across all isoforms of DCLK1 and DCLK2 (Figure 1C) and serve to extend the kinase C-spine. Mutation of these residues in DCLK1.2 uncouples tail binding and activates the kinase. Proximal to these residues are two Thr residues (Thr 687 and Thr 688), which are present in DCLK1.2, but absent in DCLK1.1. Based on published phosphoproteomics data, both Thr residues can be phosphorylated (Agulto et al., 2021) and are thus likely to be regulatory in DCLK1.2. Although we could not detect phosphorylation of either of these predicted regulatory sites in the WT form of DCLK1.2, we consistently observed pThr 688 in activated mutant DCLK1.2 variants (Figure 8—figure supplement 2). By analyzing DCLK1 dynamics using MD simulations, we observed multiple key interactions between the G-loop and the C-tail in DCLK1, such as dipole interactions with the second glycine in the G-loop by the phosphothreonine. In addition, Thr 687 contributes to increased stabilization of DCLK1.2 tail by forming a C-cap with the helical ATP-mimic segment. We aligned the intensively-studied protein kinase (PKA, PDB: 1ATP), a DCLK1.2 structure (PDB: 6KYQ), and frames of our MD trajectory, which demonstrate remarkable overlap of the ATP gamma phosphate and C-tail phosphothreonine, which seemingly acts as a mimic for the ATP co-factor. As phosphorylation is reported to lead to DCLK1 inhibition, this suggests a complex mechanism of regulation, in which the DCLK-specific constraints in the G-loop, the intrinsic flexibility of the C-tail, and threonine phosphorylation, by cis- or trans-mediated modification, systematically prevent hyperphosphorylation of the doublecortin domains and cellular effects. Somewhat paradoxically, we could only identify phosphorylated Thr 688 in activated DCLK1.2 mutants, but not in the autoinhibited (WT) version form. This suggests that the selected mutations exhibit a regulatory hierarchical dominance over inhibitory Thr 688 phosphorylation and are sufficient to liberate DCLK1.2 from its autoinhibited, C-tail-bound state. This also implies that phosphorylation of Thr 688 may only be minimally required for autoinhibition, especially given its association with the hyperactive variants obtained by mutagenesis (Figure 8—figure supplements 1–2).
Finally, we have evaluated the terminal residues in the DCLK1 C-tail, which are predicted to be intrinsically disordered. Side-by-side analysis of DCLK1.1 and DCLK1.2 in which this region is added from a common core terminating at residue 703 shows that increasing the length of the tail in both DCLK1.1 and DCLK1.2 has little additional effect on the inhibitory or stabilizing effects driven by the highly ordered tail regions that precede it. Kinase domains are regulated by IDRs in a multitude of ways, but the CAMK family is specifically enriched for adapted C-terminal extensions that, as we show here, can block the ATP and substrate binding, and enzymatically inactivate the kinase domain by occlusion of the activation loop through a flexible helical IDS on their C-tail (Gógl et al., 2019). We note that the DCLK1.2 kinase domain crystallizes in an ‘active’ closed conformation, despite binding of the C-tail in an autoinhibitory manner (Cheng et al., 2022b). Repeated packing motions of the IDS against the activation loop in all replicates of DCLK1.1 MD simulations suggest that tail may occlude the activation loop, similar to other CAMKs, possibly pointing to a mode of cis autoregulation. Conversely, AlphaFold2 predicts the placement of the DCX domains as adjacent to the IDS in both DCLK1.1 and DCLK1.2 (Figure 1—figure supplement 1). It is possible that, like other CAMKs, the IDS facilitates protein binding, whether to the DCX domains, calcium-modulated proteins, or other kinases.
Evolutionary divergence and functional specialization of DCLKs
For the DCLK family as a whole, we discovered phylogenetic divergence between DCLK1 and -2 as a relatively recent event (Figure 1A) in which metazoan DCLK3 is the more ancestral DCLK gene from which DCLK1 and -2 emerged after duplication. Because of shared evolutionary constraints and the recent divergence between DCLK1 and -2, we surmise the functional specialization of the DCLK1 tail is shared between these paralogs. Moreover, we quantify key differences between human DCLK1.1 and -1.2 activity that are impacted by amino acid changes that contribute to the function of the C-tail. The differences between DCLK1 isoforms 1 and 2 are generated by variations in exon splicing, which change both the C-tail protein sequence, and introduce or exclude potential phosphorylation sites. Expression of the highly autoinhibitable DCLK1.2 isoform is believed to be predominant in the brain during embryogenesis, although DCLK1.1 is also thought to be present in the adult brain (Burgess and Reiner, 2002). It is therefore possible that an altered ratio in DCLK1.1/1.2 expression, accompanied by changes in the requirement for DCLK1 autoregulation, is relevant for early neurogenesis. The overall sequence similarity at the protein level, despite the loss of a pair of putative phosphorylation sites, suggests a possible exon duplication in the C-tail, whereby polymorphisms have allowed for adaptive regulation during development and proliferation. Indeed, we also speculate that the induced expression of DCLK1.2 splice variants in multiple cancer subtypes (Qu et al., 2019) is likely to be indicative of a survival and drug resistance role that could be targetable with new types of small molecule. Although nanomolar DCLK1 ATP site inhibitors such as DCLK1-IN-1 have been developed that can bind tightly to the DCLK1 ATP site (Figure 3 and Figure 3—figure supplement 1), the ‘problematic’ existence of human DCLK1.1 and DCLK1.2 splice variants with distinct autoinhibitory properties may present a challenge to compound engagement in the cell, where relief of autoinhibition through C-tail undocking in DCLK1.2 is likely to require a high concentration of compound in order to compete and disengage interactions at the ATP site. Indeed, although potent chemical DCLK1 inhibitors such as DCLK1-IN-1 are known to influence DCLK1 autophosphorylation and cell motility, they have relatively modest effects in cells in terms of cytotoxicity (Ferguson et al., 2020; Ding et al., 2021). Therefore, we propose that the dual inhibitory effects of the C-tail and the transmission of this information to adjacent DCX domains, which control adaptive cellular phenotypes such as EMT in cancer cells (Major et al., 1985), may make allosteric classes of DCLK1 inhibitor a preferred therapeutic option, especially if they can be tailored specifically toward DCLK1.1 or DCLK1.2, whose autoregulation is different in terms of the varied molecular details we have uncovered here.
Materials and methods
Ortholog identification
Request a detailed protocolTo identify orthologs, we used the software KinOrtho (Huang et al., 2021) to query one-to-one orthologous relationships for DCLK1/2/3 across the proteome. After collation of the various orthologs, we parsed the sequence data for taxonomic information and classified each sequence by family. We further separated human DCLK1 into each unique isoform and aligned them.
Phylogenetic analysis
Request a detailed protocolWe identified diverse DCLK orthologs from the UniProt database (UniProt Consortium, 2021) using an profile-based approach (Neuwald, 2009). From this dataset, we manually curated a taxonomically diverse set of DCLK orthologs composed of 36 sequences spanning 16 model organisms. These sequences were used to generate a maximum-likelihood phylogenetic tree using IQTREE version 1.6.12 (Nguyen et al., 2015b). Branch support values were generated using ultrafast bootstrap with 1000 resamples (Hoang et al., 2018). The consensus tree was selected as the final tree. The optimal substitution model for our final topology was determined to be LG (Le and Gascuel, 2008) with invariable sites and discrete gamma model (Gu et al., 1995) based on the Bayesian information criterion as determined by ModelFinder (Kalyaanamoorthy et al., 2017). We rooted our final tree against an outgroup of 17 closely related human CAMK kinases using ETE Toolkit version 3.1.2 (Huerta-Cepas et al., 2016).
Sequence and structure analysis
Request a detailed protocolMAFFT (Katoh et al., 2002) generated multiple sequence alignments were fed into the BPPS tool to determine evolutionarily conserved and functionally significant residues (Neuwald, 2014). Constraints mapped onto ΑlphaFold-predicted structures were visualized in PyMOL to analyze biochemical interactions.
Rosetta loop modeling
Request a detailed protocolLoop modeling was performed on the crystal structure (6KYQ) using the kinematic loop modeling protocol (Mandell et al., 2009) to model missing residues. Following this, the structure underwent five cycles of rotamer repacking and minimization using the Rosetta Fast-relax protocol (Tyka et al., 2011).
DSSP analysis
Request a detailed protocolTo analyze changes in secondary structure over our MDs, we employed the DSSP (Kabsch and Sander, 1983) command in GROMACS. This produces an output that contains an array of secondary structure values against each residue. The MDAnalysis python module (Michaud-Agrawal et al., 2011) was used to plot these values.
Molecular dynamics
Request a detailed protocolPDB constructs were generated by retrieving structural models RCSB and the ΑlphaFold2 database. Post-translational modifications were performed in PyMOL using the PyTMs plugin. All structures were solvated using the TIP3P water model (Jorgensen et al., 1983). Energy minimization was run for a maximum of 10,000 steps, performed using the steepest descent algorithm, followed by the conjugate gradient algorithm. The system was heated from 0 K to a temperature of 300 K. After two equilibration steps that each lasted 20 ps, 1-µs-long simulations were run at a 2 fs timestep. Long-range electrostatics were calculated via particle mesh Ewald algorithms using the GROMACS MD engine (Pronk et al., 2013). We utilized the CHARMM36 force field (Huang and MacKerell, 2013). The resulting output was visualized using VMD 1.9.3 (Humphrey et al., 1996). All MD analysis was conducted using scripts coded in Python using the MDAnalysis module (Michaud-Agrawal et al., 2011).
Computational mutational analysis
Request a detailed protocolCartesian ddG in Rosetta (Park et al., 2016) was utilized to predict potential stabilizing and destabilizing mutations in the enzyme structure. We performed three replicates per mutation and averaged the Rosetta energies. All mutant energies were then subtracted by the WT Rosetta energy to generate a panel of ddG values relative to WT. Combined with our sequence analyses, we mutated kinase and DCLK-specific constraints to identify destabilizing interactions in the C-tail.
Exon-intron boundary mapping
Request a detailed protocolThe precise gene structure of DCLK1 isoforms were mapped onto the human genome with each isoform used as a query protein sequence in order to generate exon-intron borders. This was achieved using Scipio (version 1.4.1) (Keller et al., 2008) with default settings. Exons were numbered based on Ensembl annotations (Cunningham et al., 2022). The translation of each annotated gene sequence to protein sequence was provided with the output file (Figure 4—source data 1–4).
DCLK1 cloning and recombinant protein expression
Request a detailed protocol6His-DCLK1 catalytic domain (351–686), DCLK1.1 351–703 (short C-tail) or 351–729 (full C-tail), DCLK1.2 351–703 (short C-tail) or 351–740 (full C-tail), and DCLK1.2 351–703 containing D398A, N400A, V682T/V684T, or D691A substitutions were synthesized by Twist Biosciences in pET28a. 6His-GST-(3C) DCLK1.1 351–689 was amplified by PCR and cloned into pOPINJ. Kinase dead, D533A 6His-GST-(3C) DCLK1.1 351–689 was generated by PCR-based site-directed mutagenesis (Figure 3—source data 3). All plasmids were sequenced prior to their use in protein expression studies. All proteins, including 6His-GST-(3C) DCLK1 351–689, with a 3C protease cleavable affinity tag, were expressed in BL21(DE3)pLysS Escherichia coli (Novagen) and purified by affinity and SEC. The short N-terminal 6His affinity tag present on all other DCLK1 proteins described in this paper was left in situ on recombinant proteins, since it does not appear to interfere with differential scanning fluorimetry (DSF), biochemical interactions, or catalysis. For analytical SEC (Figure 3—figure supplement 1—source data 1, Figure 3—figure supplement 2—source data 1), 1 mg of each DCLK1 protein was assayed on a Superdex 200 Increase 10/300 GL (Cytiva), and the eluted fractions were also analyzed by SDS-PAGE and Coomassie blue staining to confirm composition. The molecular weight standards were loaded in a mixture of 200 µg of bovine serum albumin, carbonic anhydrase, and alcohol dehydrogenase each.
Mass spectrometry
Request a detailed protocolPurified DCLK1 proteins (5 μg) were diluted (~40-fold) in 100 mM ammonium bicarbonate pH 8.0 and reduced (DTT) and alkylated iodoacetamide, as previously described (Ferries et al., 2017), and digested with a 25:1 (wt/wt) trypsin gold (Promega) at 37°C for 18 hr with gentle agitation. Digests were then subjected to strong cation exchange chromatography using in-house packed stage tip clean-up (Daly et al., 2021). Dried tryptic peptides were solubilized in 20 μl of 3% (vol/vol) acetonitrile and 0.1% (vol/vol) TFA in water, sonicated for 10 min, and centrifuged at 13,000 × g for 10 min at 4°C and supernatant collected. LC-MS/MS separation was performed using an Ultimate 3000 nano system (Dionex), over a 60 min gradient (Ferries et al., 2017). Briefly, samples were loaded at a rate of 12 μl/min onto a trapping column (PepMap100, C18, 300 μm×5 mm) in loading buffer (3% [vol/vol] acetonitrile, 0.1% [vol/vol] TFA). Samples were then resolved on an analytical column (Easy-Spray C18 75 μm×500 mm, 2 μm bead diameter column) using a gradient of 97% A (0.1% [vol/vol] formic acid): 3% B (80% [vol/vol] acetonitrile, 0.1% [vol/vol] formic acid) to 80% B over 30 min at a flow rate of 300 nl/min. All data acquisition was performed using a Fusion Lumos Tribrid mass spectrometer (Thermo Scientific). Samples were injected twice with either higher-energy C-trap dissociation (HCD) fragmentation (set at 32% normalized collision energy [NCE]) or electron transfer dissociation with supplemental 30% NCE HCD (EThcD) for 2+ to 4+ charge states using a top 3 s top speed mode. MS1 spectra were acquired at a 120 K resolution (at 200 m/z), over a range of 300–2000 m/z, normalized AGC target = 50%, maximum injection time = 50 ms. MS2 spectra were acquired at a 30 K resolution (at 200 m/z), AGC target = standard, maximum injection time = dynamic. A dynamic exclusion window of 20 s was applied at a 10 ppm mass tolerance. Data was analyzed by Proteome Discoverer 2.4 in conjunction with the MASCOT search engine using a custom database of the UniProt E. coli reviewed database (updated January 2023) with the DCLK1 mutant variant amino acid sequences manually added, and using the search parameters: fixed modification = carbamidomethylation (C), variable modifications = oxidation (M) and phospho (S/T/Y), MS1 mass tolerance = 10 ppm, MS2 mass tolerance = 0.01 Da, and the ptmRS node on; set to a score >99.0. For HCD data, instrument type = electrospray ionization-Fourier transform ion cyclotron resonance, for EThcD data, instrument type = EThcD. For label-free relative quantification of phosphopeptide abundances of the different DCLK1 variants, the minora feature detector was active and set to calculate the area under the curve for peptide m/z ions. Abundance of phosphopeptide ions were normalized against the total protein abundance (determined by the HI3 method [Silva et al., 2006], as in the minora feature detector node) to account for potential protein load variability during analysis.
DCLK1 DSF
Request a detailed protocolThermal shift assays were performed using DSF in a StepOnePlus Real-Time PCR machine (Life Technologies) in combination with Sypro-Orange dye (Invitrogen) and a thermal ramping protocol (0.3°C per minute between 25°C and 94°C). Recombinant DCLK1 proteins were assayed at a final concentration of 5 μM in 50 mM Tris-HCl (pH 7.4) and 100 mM NaCl in the presence or absence of the indicated concentrations of ligand (ATP or Mg:ATP) or DCLK1 inhibitor compounds, with final DMSO concentrations never higher than 4% (vol/vol). Thermal melting data were processed using the Boltzmann equation to generate sigmoidal denaturation curves, and average Tm/ΔTm values were calculated as described using GraphPad Prism software, as previously described, from three technical repeats (Byrne et al., 2020a).
DCLK1 kinase assays
Request a detailed protocolDCLK1 peptide-based enzyme assays (Byrne et al., 2016; Omar et al., 2023) were carried out using the LabChip EZ Reader platform, which monitors and quantifies real-time phosphorylation-induced changes in the mobility of the fluorescently labeled DCLK1 peptide substrate 5-FAM-KKALRRQETVDAL-CONH2. To assess DCLK1 catalytic domains, or DCLK1.1 or DCLK1.2 variants, 100 ng of purified protein were incubated with a high (1 mM) concentration of ATP (to mimic cellular levels of nucleotide) and 2 μM of the fluorescent substrate in 25 mM HEPES (pH 7.4), 5 mM MgCl2, and 0.001% (vol/vol) Brij 35. DCLK1-IN-1 and DCLK1-NEG (kind gifts from Dr Fleur Ferguson, UCSF) enzyme inhibition was quantified under identical assay conditions in the presence of 10 μM of each compound. Assays are either reported as rates (pmoles/min phosphate incorporation) during linear phosphate incorporation (e.g. total substrate phosphorylation limited to <20–30% to prevent ATP depletion and to ensure assay linearity) or presented as time-dependent percentage substrate phosphorylation (kinetic mode). Rates of substrate phosphorylation (pmol phosphate incorporation per minute) were determined using a fixed amount of kinase and linear regression analysis with GraphPad Prism software; Vmax and KM[ATP] values were calculated at 2 µM substrate peptide concentration, as previously described (McSkimming et al., 2016). Rates are normalized to enzyme concentration and all enzyme rate and kinetic data are presented as mean and SD of four technical replicates.
Data availability
All data generated in this study are included within the manuscript. Source data are provided for each figure. MD simulations and associated data have been deposited to Dryad (https://doi.org/10.5061/dryad.8931zcrxb).
-
Dryad Digital RepositoryMolecular Dynamics Simulations and associated data for Mechanistic and evolutionary insights into isoform-specific 'supercharging' in DCLK family kinases.https://doi.org/10.5061/dryad.8931zcrxb
References
-
mTOR Regulation of AGC Kinases: New Twist to an Old TailMolecular Pharmacology 101:213–218.https://doi.org/10.1124/molpharm.121.000310
-
Alternative splice variants of doublecortin-like kinase are differentially expressed and have different kinase activitiesThe Journal of Biological Chemistry 277:17696–17705.https://doi.org/10.1074/jbc.M111981200
-
Research Progress of DCLK1 Inhibitors as Cancer TherapeuticsCurrent Medicinal Chemistry 29:2261–2273.https://doi.org/10.2174/0929867328666210709110721
-
Doublecortin expression levels in adult brain reflect neurogenesisThe European Journal of Neuroscience 21:1–14.https://doi.org/10.1111/j.1460-9568.2004.03813.x
-
TRIBBLES: a twist in the pseudokinase tailStructure 23:1974–1976.https://doi.org/10.1016/j.str.2015.10.003
-
Discovery of a selective inhibitor of doublecortin like kinase 1Nature Chemical Biology 16:635–643.https://doi.org/10.1038/s41589-020-0506-0
-
DCLK1 is up-regulated and associated with metastasis and prognosis in colorectal cancerJournal of Cancer Research and Clinical Oncology 142:2131–2140.https://doi.org/10.1007/s00432-016-2218-0
-
Disordered protein kinase regions in regulation of kinase domain coresTrends in Biochemical Sciences 44:300–311.https://doi.org/10.1016/j.tibs.2018.12.002
-
Maximum likelihood estimation of the heterogeneity of substitution rate among nucleotide sitesMolecular Biology and Evolution 12:546–557.https://doi.org/10.1093/oxfordjournals.molbev.a040235
-
Analysis of human Tribbles 2 (TRIB2) pseudokinaseMethods in Enzymology 667:79–99.https://doi.org/10.1016/bs.mie.2022.03.025
-
UFBoot2: improving the ultrafast bootstrap approximationMolecular Biology and Evolution 35:518–522.https://doi.org/10.1093/molbev/msx281
-
Doublecortin, a stabilizer of microtubulesHuman Molecular Genetics 8:1599–1610.https://doi.org/10.1093/hmg/8.9.1599
-
CHARMM36 all-atom additive protein force field: validation based on comparison to NMR dataJournal of Computational Chemistry 34:2135–2145.https://doi.org/10.1002/jcc.23354
-
Structure-function of the multifunctional Ca2+/calmodulin-dependent protein kinase IIThe Biochemical Journal 364:593–611.https://doi.org/10.1042/BJ20020228
-
ETE 3: reconstruction, analysis, and visualization of phylogenomic dataMolecular Biology and Evolution 33:1635–1638.https://doi.org/10.1093/molbev/msw046
-
VMD: visual molecular dynamicsJournal of Molecular Graphics 14:33–38.https://doi.org/10.1016/0263-7855(96)00018-5
-
Comparison of simple potential functions for simulating liquid waterThe Journal of Chemical Physics 79:926–935.https://doi.org/10.1063/1.445869
-
MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transformNucleic Acids Research 30:3059–3066.https://doi.org/10.1093/nar/gkf436
-
Coupled regulation by the juxtamembrane and sterile α motif (SAM) linker is a hallmark of ephrin tyrosine kinase evolutionThe Journal of Biological Chemistry 293:5102–5116.https://doi.org/10.1074/jbc.RA117.001296
-
An improved general amino acid replacement matrixMolecular Biology and Evolution 25:1307–1320.https://doi.org/10.1093/molbev/msn067
-
KinView: a visual comparative sequence analysis tool for integrated kinome researchMolecular bioSystems 12:3651–3665.https://doi.org/10.1039/c6mb00466k
-
MDAnalysis: a toolkit for the analysis of molecular dynamics simulationsJournal of Computational Chemistry 32:2319–2327.https://doi.org/10.1002/jcc.21787
-
A Bayesian sampler for optimization of protein domain hierarchiesJournal of Computational Biology 21:269–286.https://doi.org/10.1089/cmb.2013.0099
-
IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogeniesMolecular Biology and Evolution 32:268–274.https://doi.org/10.1093/molbev/msu300
-
Molecular identification and characterization of a family of kinases with homology to Ca2+/calmodulin-dependent protein kinases I/IVThe Journal of Biological Chemistry 281:20427–20439.https://doi.org/10.1074/jbc.M513212200
-
Classification of Cushing’s syndrome PKAc mutants based upon their ability to bind PKIThe Biochemical Journal 480:875–890.https://doi.org/10.1042/BCJ20230183
-
Expression and chromosomal localization of KIAA0369, a putative kinase structurally related to DoublecortinJournal of Human Genetics 43:169–177.https://doi.org/10.1007/s100380050063
-
Simultaneous optimization of biomolecular energy functions on features from small molecules and macromoleculesJournal of Chemical Theory and Computation 12:6201–6212.https://doi.org/10.1021/acs.jctc.6b00819
-
A chimeric mechanism for polyvalent trans-phosphorylation of PKA by PDK1Protein Science 18:1486–1497.https://doi.org/10.1002/pro.146
-
Absolute quantification of proteins by LCMSE: a virtue of parallel MS acquisitionMolecular & Cellular Proteomics 5:144–156.https://doi.org/10.1074/mcp.M500230-MCP200
-
Protein kinases: evolution of dynamic regulatory proteinsTrends in Biochemical Sciences 36:65–77.https://doi.org/10.1016/j.tibs.2010.09.006
-
Alternate states of proteins revealed by detailed energy landscape mappingJournal of Molecular Biology 405:607–618.https://doi.org/10.1016/j.jmb.2010.11.008
-
UniProt: the universal protein knowledgebase in 2021Nucleic Acids Research 49:D480–D489.https://doi.org/10.1093/nar/gkaa1100
-
Structural examination of autoregulation of multifunctional calcium/calmodulin-dependent protein kinase IIThe Journal of Biological Chemistry 274:26199–26208.https://doi.org/10.1074/jbc.274.37.26199
-
Evolution of functional diversity in the holozoan tyrosine kinomeMolecular Biology and Evolution 38:5625–5639.https://doi.org/10.1093/molbev/msab272
-
The regulatory region of calcium/calmodulin-dependent protein kinase I contains closely associated autoinhibitory and calmodulin-binding domainsThe Journal of Biological Chemistry 270:23851–23859.https://doi.org/10.1074/jbc.270.40.23851
-
2.2 A refined crystal structure of the catalytic subunit of cAMP-dependent protein kinase complexed with MnATP and A peptide inhibitorActa Crystallographica. Section D, Biological Crystallography 49:362–365.https://doi.org/10.1107/S0907444993000423
Peer review
Reviewer #1 (Public Review):
In the study by Venkat et al. the authors expand the current knowledge of allosteric diversity in the human kinome by c-terminal splicing variants using as a paradigm DCLK1. In this work, the authors provide evolutionary and some mechanistic evidence about how c-terminal isoform specific variants generated by alternative splicing can regulate catalytic activity by means of coupling specific phosphorylation sites to dynamical and conformational changes controlling active site and substrate pocket occupancy, as well as interfering with protein-protein interacting interfaces that altogether provides evidence of c-terminal isoform specific regulation of the catalytic activity in protein kinases.
The paper is overall well written, the rationale and the fundamental questions are clear and well explained, the evolutionary and MD analyses are very detailed and well explained. Overall I think this is a study that brings some new aspects and concepts that are important for the protein kinase field, in particular the allosteric regulation of the catalytic core by c-terminal segments, and how evolutionary cues generate more sophisticated mechanisms of allosteric control in protein kinases.
Current submission: I have read and gone through the revised manuscript and the rebuttal letter and I confirm that the authors did an excellent job answering all the comments satisfactorily.
https://doi.org/10.7554/eLife.87958.3.sa1Reviewer #2 (Public Review):
In this study, the authors explore the structure/function of the DCLK kinases, most specifically DCLK1 as it is the most studied to date. Recently, the C-terminal domain has garnered attention as it was found to regulate the kinase domain, however, the different isoforms retain additional amino acid sequences with as-yet-undefined functions. The authors provide an evolutionary and biochemical characterization of these regions and provide evidence for some functionality for these additional C-terminal sequences. While these experiments are informative they do require that the protein is soluble and not membrane-bound as has been suggested to be important for functionality in other studies. Still, this is a major contribution to understanding the structure/function of these proteins that will be important in future experimental designs.
https://doi.org/10.7554/eLife.87958.3.sa2Author response
The following is the authors’ response to the original reviews.
eLife assessment
This important study expands on current knowledge of allosteric diversity in the human kinome by C-terminal splicing variants using as a paradigm DCLK1. The authors provide solid evolutionary and some mechanistic evidence how C-terminal isoform specific variants generated by alternative splicing can regulate catalytic activity by means of coupling specific phosphorylation sites to dynamical and conformational changes controlling active site and substrate pocket occupancy, as well as protein-protein interactions. The data will be of interest to researchers in the kinase and signal transduction field.
We thank the editor for coordinating the review of our manuscript and the reviewers for their valuable feedback. We have significantly revised the manuscript in response to the reviewer’s comments. Our point-by-point response to each comment is present below. We have uploaded both a clean draft of our revised manuscript as well as a version with the revisions highlighted in yellow. We hope the revised manuscript is now acceptable for publication in eLife. We have additionally updated the preprint on bioRxiv and have included the link: We thank the editor for coordinating the review of our manuscript and the reviewers for their valuable feedback. We have significantly revised the manuscript in response to the reviewer’s comments. Our point-by-point response to each comment is present below. We have uploaded both a clean draft of our revised manuscript as well as a version with the revisions highlighted in yellow. We hope the revised manuscript is now acceptable for publication in eLife. We have additionally updated the preprint on biorxiv and have included the link here: https://www.biorxiv.org/content/10.1101/2023.03.29.534689v2.
Reviewer #1
Summary
In the study by Venkat et al. the authors expand the current knowledge of allosteric diversity in the human kinome by c-terminal splicing variants using as a paradigm DCLK1. In this work, the authors provide evolutionary and some mechanistic evidence about how c-terminal isoform specific variants generated by alternative splicing can regulate catalytic activity by means of coupling specific phosphorylation sites to dynamical and conformational changes controlling active site and substrate pocket occupancy, as well as interfering with protein-protein interacting interfaces that altogether provides evidence of c-terminal isoform specific regulation of the catalytic activity in protein kinases.
The paper is overall well written, the rationale and the fundamental questions are clear and well explained, the evolutionary and MD analyses are very detailed and well explained. The methodology applied in terms of the biochemical and biophysical tools falls a bit short in some places and some comments and suggestions are given in this respect. If the authors could monitor somehow protein auto-phosphorylation as a functional readout would be very useful by means of using phospho-specific antibodies to monitor activity. Overall I think this is a study that brings some new aspects and concepts that are important for the protein kinase field, in particular the allosteric regulation of the catalytic core by c-terminal segments, and how evolutionary cues generate more sophisticated mechanisms of allosteric control in protein kinases. However a revision would be recommended.
Major Comments
The authors explain in the introduction the role of T688 autophosphorylation site in the function of DCLK1.2. This site when phosphorylated have a detrimental impact on catalytic activity and inhibits phosphorylation of the DCX domain. allowing the interaction with microtubules. In the paper they show how this site is generated by alternative splicing and intron skipping in DCLK1.2. However there is no further functional evidence along the functional experiments presented in this study.
1. What is the effect of a non-phosphorylable T688 mutant in terms of stability and enzymatic activity? What would be the impact of this mutant in the overall auto-phosphorylation reaction?
The role of T688 phosphorylation on DCLK1 functions has been explored in previous studies (Agulto et al, 2020: PMID: 34310279), although only relevant to DCLK1.2 splice variants, since this site is lacking in DCLK1.1. These studies showed that mutation of T688 to an alanine increases total kinase autophosphorylation (ie autoactivity) and the subsequent phosphorylation of DCX domains, which in turn decreases microtubule binding. Given this information, our goal was to use an evolutionary perspective to investigate this, alongside less-well characterized aspects of DCLK autoregulation, including co-conserved residues in the catalytic domain and C-terminal tail. However, to address the reviewers question of a non-phosphorylatable T688 mutant, we performed MD simulations of T688A and T688E (a phosphomimic) mutant and include a new supplementary figure (Figure 5-supplement 3) which show the two mutants slightly destabilize the C-tail relative to wt (1 and 2 angstrom increase in RMSF for T688E and T688A respectively), but by themselves cannot dislodge the C-tail from the ATP binding pocket. Thus, other co-conserved interactions as revealed by our analysis, are likely to contribute to the autoregulation of the kinase domain by the C-terminal tail. We have incorporated these observations into the revised results section.
Furthermore, to address the reviewer’s question in terms of site-specific autophosphorylation as a marker of DCLK1.2 activity, we have now performed a much-more detailed phosphoproteomic analysis of a panel of purified DCLK1.2 proteins after purification from E. coli (Figure 8-figure supplement 2). This showed that we are only able to detect Thr 688 phosphorylated in our ‘activated’ DCLK1.2 mutants, and not in the autoinhibited WT DCLK1.2 version of the protein. This apparent contradiction does not necessary discount Thr 688 as an important regulatory hotspot, but, together with the MD simulations, may imply a decreased contribution of pThr 688 in facilitating/maintaining DCLK1.2 auto-inhibition than previously anticipated, especially in the context of the numerous other stabilizing amino acid contacts that we describe between the C-tail and the ATP-binding pocket. We do, however, propose a mechanism for pThr688 as a potential ATP mimic based on MD analysis. However, we only found MS-based evidence for phosphorylation at this (and other sites in the same peptide) in highly active DCLK1.2 mutants, in which the C-tail remains uncoupled from the ATP-binding site, even in the presence of this regulatory PTM. We acknowledge that better understanding of DCLK biology will require a detailed appraisal of how the DCLK auto-inhibited states are subsequently physiologically regulated (PTMs, protein-protein interaction etc.), but this is beyond the scope of our current evolutionary investigation, and the absence of phosphospecific antibodies makes this challenging currently. We intend to expand upon our current work by assessing the relative contribution of multiple DCLK phosphorylation sites (including, but not limited to, Thr 688) with regard to cellular DCLK auto-regulation in future studies, in part by generating such site-specific phospho-antibodies.
1. Have the authors made an equivalent T687/688 tanden in DCLK1.1 instead of the two prolines?
This is a good point. We have not considered introducing a T687/688 tandem mutation into DCLK1.1 (at the equivalent position to that of DCLK1.2), primarily because the amino acid composition of their respective C-tail domains are so highly divergent across the tail (due to alternative splicing, as discussed in our paper). As discussed in our present study, there are numerous contacts made between specific amino acids in the regulatory C-tail and the kinase domain of DCLK1.2, which functionally occlude ATP binding, and thus change catalytic output. It is these contacts, which are determined by the specific amino acid sequence identity, and not the extended length of the DCLK1.2 C-tail per se, that drives autoinhibition. The alternate amino acid sequence identity of the C-tail of DCLK1.1 does not enable such contacts to form, which we believe explains the different activities of the two isoforms.
Furthermore, our mutational analysis reveals clearly that Thr688 and several other sites are more highly autophosphorylated in the artificially activated DCLK1.2 constructs than WT DCLK1.2, and as such it remains our hypothesis that introduction of the tandem phosphorylation sites into DCLK1.1 is unlikely to be sufficient to impose an auto-inhibitory conformation of the enzyme.
1. Could T688 autophosphorylation be used as a functional readout to evaluate DCLK1.2 activity?
We agree with the reviewer’s suggestion about using autophosphorylation (including potentially Thr688 for DCLK1.2) as a functional read out for DCLK1 activity. In our present study, we identify phosphorylated peptides containing pThr688 only in the mutationally activated DCLK1.2 variants. We have now taken this analytical approach further and performed a detailed comparative phosphoproteomic characterisation of all of our DCLK1 constructs, where we observe marked differences in the overall phosphorylation profiles of the mutant DCLK1.2 (and DCLK1.1) proteins relative to the less phosphorylated WT DCLK1.2 kinase. This manifests as a depletion in the total number of confidently assigned phosphorylation sites within the kinase domain and C-tail of WT DCLK1.2, and also as a depletion in the abundance of phosphorylated peptides for a given site. To help visualise this, individual phosphorylation sites have been schematically mapped onto DCLK1, which has been included as a new extended supplementary figure (Figure 8-figure supplement 2). For comparative analysis of phosphosite abundance, we could only select peptides that could be directly compared between all mutants (identical amino acid sequences) and those found to be phosphorylated in all proteins (these are Ser660 and Thr438); these are now shown in figure supplement 2 as a table. These site occupancies follow what we see with respect to the increased catalytic activity between DCLK1.1 and DCLK1.2 mutants versus DCLK1.2. We also detect increased phosphorylation of DCLK1.1 and activated DCLK1.2 mutants in comparison to (autoinhibited) DCLK1.2, supporting the hypothesis that these mutants are relieving the autoinhibited conformation.
1. What are the evidences of the here described c-terminal specific interactions to be intra-molecular rather than inter-molecular? Have the authors looked at the monodispersion and molecular mass in solution of the different protein evaluated in this study? Basically, are the proteins in solutions monomers or dimers/oligomers?
Analysis of symmetry mates in the crystal structure of DCLK1.2 (PDB ID: 6KYQ) provide no evidence for inter-molecular interactions. Furthermore, to evaluate oligomerization status in solution, we conducted an analytical size exclusion chromatography (SEC) and our analysis reveals that both DCLK1.1 and DCLK1.2 predominantly exist as monomers in solution (Figure 3-Supplements 1-3). These results suggest that the C-terminal tail interactions are primarily intra-molecular.
1. (Figure 3) Did the authors look at the mono-dispersion of the protein preparation? The sec profile did result in one single peak or multiple peaks? Could the authors show the chromatogram? how many species do you have in solution? Was the tag removed from the recombinant proteins or not?
Yes, as mentioned above, the SEC profile resulted in a single peak for both DCLK1.1 and DCLK1.2, which was confirmed as DCLK1 by subsequent SDS-PAGE. We have included the chromatogram and gels in supporting data (Figure 3-supplements 1-3) in the revised manuscript and updated the Methods section. ‘The short N-terminal 6-His affinity tag present on all other DCLK1 proteins described in this paper was left in situ on recombinant proteins, since it does not appear to interfere with DSF, biochemical interactions or catalysis.’
1. Authors should do Michaelis-Menten saturation kinetics as shown in Figure 3C with the WT when comparing all the functional variant analysed in the study. So we can compared the catalytic rates and enzymatic constants (depicted in a table also) kcat, Km and catalytic efficiency constants (kcat/Km)
Thank you for your suggestion. We have performed the requested comparative kinetics analyses for selected functional DCLK1 variants at the same concentration as suggested, using our real-time assay to determine Vmax for peptide phosphorylation as a function of ATP, but at a fixed substrate concentration (we are unable to assess Vmax above 5 µM peptide for technical reasons). The results of these analyses have been included in the revised version of Figure 8-Supplement 1, where they support differences in both Vmax and Km[ATP]; the ratio of these values very clearly points to differences in activities falling into ‘low’ or ‘high’. This kinetic analysis fully supports our initial activity assays, where mutations predicted to uncouple the auto-inhibitory C-tail rescue DCLK1.2 activity to levels similar to DCLK1.1 towards a common substrate.
Minor Comments
It is very interesting how the IBS together with the pT688 mimics ATP in the case of DCLK1.2 to reach full occupancy of the active site. On Figure 8 you evaluate residues of the GRL and IBS interface to probe such interactions.
1. Did the authors look at the T688 non-phosphorylable mutant?
See our response to Major Comment 1 above. In addition, due to the absence of T688 in DCLK1.1, we did not look at the T688A mutant of DCLK1.2 biochemically, partially because it has been characterized in previous studies, but partially because this site is preceeded by another Thr residue. The lack of a selective antibody towards this site makes it difficult to evaluate the role of T688 phosphorylation specifically with respect to DCLK cellular functions and interactions. Therefore, we focused our in vitro efforts to understand how mutations in the IBS impact the catalytic activity of DCLK1.2 by comparing different variants to DCLK1.1.
1. Classification of DCLK C-terminal regulatory elements.
It would be useful to connect the different regulatory elements described in this study to a specific functional and biological setting where these different switches play a role e.g. microtubule interactions and dynamics, cell cycle, cancer, etc..
While the primary focus of our paper is on the mechanism of allosteric regulation of DCLK1, we have indeed touched upon the potential implications of the various regulatory elements of the tail on functions such as microtubule binding and phenotypic effects like cancer progression. However, we acknowledge that a comprehensive understanding of these effects would necessitate a more detailed investigation. This could potentially involve the integration of RNA-seq data with extensive cell assays to evaluate phenotypic effects. We believe that such a future study would be a valuable extension of our current work and could provide further insights into the functional roles of DCLK1.
1. (Figure 3) Could the authors explain the differences in yield between the WT and the D531A mutant. Apparently, it [the yield] does not appear to be caused by a lower stability as indicated by the Tm. Could the authors comment on this?It is important to compare different samples in parallel, in the same experiment and side by side. This applies to the thermal shift data comparing WT and a D531A mutant on panel D and also on panel C a comparison between WT and D531A as negative control should be shown.
WT and D533A (kinase-dead) were indeed analysed in parallel, but have been split in two panels to make the data easier to interpret. The modest differences in yield is likely explained by experimental prep-to-prep variations. Our experience shows that many protein kinase yields vary between kinase and kinase-dead variants, likely due to bacterial toxicity related to enzyme activity. In regards to thermal stability, we would like to emphasize that Differential Scanning Fluorimetry (DSF) is to our mind a more informative and quantitative measure of protein stability than yield from bacteria, because both assess purified proteins at the same concentration. We believe that the DSF data provide a more accurate representation of the real stability differences between the WT and D533A mutant.
https://doi.org/10.7554/eLife.87958.3.sa3Article and author information
Author details
Funding
National Institute of General Medical Sciences (GM139656)
- Natarajan Kannan
North West Cancer Research (CR1208)
- Leonard A Daly
- Claire E Eyers
- Patrick A Eyers
- Dominic P Byrne
Biotechnology and Biological Sciences Research Council (BB/S018514/1)
- Dominic P Byrne
- Leonard A Daly
- Patrick A Eyers
- Claire E Eyers
Biotechnology and Biological Sciences Research Council (BB/N021703/1)
- Dominic P Byrne
- Leonard A Daly
- Patrick A Eyers
- Claire E Eyers
Biotechnology and Biological Sciences Research Council (BB/X002780/1)
- Dominic P Byrne
- Leonard A Daly
- Patrick A Eyers
- Claire E Eyers
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
Funding from NK (grant no. R35 GM139656) is acknowledged. PAE acknowledges funding from a University of Liverpool BBSRC IAA award. AV acknowledges funding from ARCS Foundation. DPB, LAD, PAE, and CEE also acknowledge BBSRC grants BB/S018514/1, BB/N021703/1, and BB/X002780/1 and North West Cancer Research (NWCR) grant CR1208. EEF thanks the MRC for a DiMeN DTP studentship (No. 1961582). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
Senior and Reviewing Editor
- Volker Dötsch, Goethe University, Germany
Version history
- Preprint posted: March 29, 2023 (view preprint)
- Sent for peer review: March 29, 2023
- Preprint posted: June 1, 2023 (view preprint)
- Preprint posted: September 5, 2023 (view preprint)
- Version of Record published: October 26, 2023 (version 1)
Cite all versions
You can cite all versions using the DOI https://doi.org/10.7554/eLife.87958. This DOI represents all versions, and will always resolve to the latest one.
Copyright
© 2023, Venkat, Watterson, Byrne et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 181
- Page views
-
- 19
- Downloads
-
- 0
- Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Computational and Systems Biology
- Physics of Living Systems
Diabetes is caused by the inability of electrically coupled, functionally heterogeneous -cells within the pancreatic islet to provide adequate insulin secretion. Functional networks have been used to represent synchronized oscillatory [Ca2+] dynamics and to study -cell subpopulations, which play an important role in driving islet function. The mechanism by which highly synchronized -cell subpopulations drive islet function is unclear. We used experimental and computational techniques to investigate the relationship between functional networks, structural (gap-junction) networks, and intrinsic -cell dynamics in slow and fast oscillating islets. Highly synchronized subpopulations in the functional network were differentiated by intrinsic dynamics, including metabolic activity and KATP channel conductance, more than structural coupling. Consistent with this, intrinsic dynamics were more predictive of high synchronization in the islet functional network as compared to high levels of structural coupling. Finally, dysfunction of gap junctions, which can occur in diabetes, caused decreases in the efficiency and clustering of the functional network. These results indicate that intrinsic dynamics rather than structure drive connections in the functional network and highly synchronized subpopulations, but gap junctions are still essential for overall network efficiency. These findings deepen our interpretation of functional networks and the formation of functional sub-populations in dynamic tissues such as the islet.
-
- Computational and Systems Biology
- Immunology and Inflammation
T cells are required to clear infection, and T cell motion plays a role in how quickly a T cell finds its target, from initial naive T cell activation by a dendritic cell to interaction with target cells in infected tissue. To better understand how different tissue environments affect T cell motility, we compared multiple features of T cell motion including speed, persistence, turning angle, directionality, and confinement of T cells moving in multiple murine tissues using microscopy. We quantitatively analyzed naive T cell motility within the lymph node and compared motility parameters with activated CD8 T cells moving within the villi of small intestine and lung under different activation conditions. Our motility analysis found that while the speeds and the overall displacement of T cells vary within all tissues analyzed, T cells in all tissues tended to persist at the same speed. Interestingly, we found that T cells in the lung show a marked population of T cells turning at close to 180o, while T cells in lymph nodes and villi do not exhibit this “reversing” movement. T cells in the lung also showed significantly decreased meandering ratios and increased confinement compared to T cells in lymph nodes and villi. These differences in motility patterns led to a decrease in the total volume scanned by T cells in lung compared to T cells in lymph node and villi. These results suggest that the tissue environment in which T cells move can impact the type of motility and ultimately, the efficiency of T cell search for target cells within specialized tissues such as the lung.