Evolution of functional diversity in the holozoan tyrosine kinome

An atlas of bacterial serine-threonine kinases reveals functional diversity and key distinctions from eukaryotic kinases

Bacterial serine-threonine kinases (STKs) regulate diverse cellular processes associated with cell growth, virulence, and pathogenicity and are evolutionarily related to the druggable eukaryotic STKs. A deeper understanding of how bacterial STKs differ from their eukaryotic counterparts and how they have evolved to regulate diverse bacterial signaling functions is crucial for advancing the discovery and development of new antibiotic therapies. Here, we classified more than 300,000 bacterial STK sequences from the NCBI RefSeq nonredundant and UniProt protein databases into 35 canonical and seven pseudokinase families on the basis of the patterns of evolutionary constraints in the conserved catalytic domain and flanking regulatory domains. Through statistical comparisons, we identified features distinguishing bacterial STKs from eukaryotic STKs, including an arginine residue in a regulatory helix (C helix) that dynamically couples the ATP- and substrate-binding lobes of the kinase domain. Biochemical and peptide library screens demonstrated that evolutionarily constrained residues contributed to substrate specificity and kinase activation in the Mycobacterium tuberculosis kinase PknB. Together, these findings open previously unidentified avenues for investigating bacterial STK functions in cellular signaling and for developing selective bacterial STK inhibitors.

Redox regulation and dynamic control of brain-selective kinases BRSK1/2 in the AMPK family through cysteine-based mechanisms

In eukaryotes, protein kinase signaling is regulated by a diverse array of post-translational modifications, including phosphorylation of Ser/Thr residues and oxidation of cysteine (Cys) residues. While regulation by activation segment phosphorylation of Ser/Thr residues is well understood, relatively little is known about how oxidation of cysteine residues modulate catalysis. In this study, we investigate redox regulation of the AMPK-related brain-selective kinases (BRSK) 1 and 2, and detail how broad catalytic activity is directly regulated through reversible oxidation and reduction of evolutionarily conserved Cys residues within the catalytic domain. We show that redox-dependent control of BRSKs is a dynamic and multilayered process involving oxidative modifications of several Cys residues, including the formation of intramolecular disulfide bonds involving a pair of Cys residues near the catalytic HRD motif and a highly conserved T-loop Cys with a BRSK-specific Cys within an unusual CPE motif at the end of the activation segment. Consistently, mutation of the CPE-Cys increases catalytic activity in vitro and drives phosphorylation of the BRSK substrate Tau in cells. Molecular modeling and molecular dynamics simulations indicate that oxidation of the CPE-Cys destabilizes a conserved salt bridge network critical for allosteric activation. The occurrence of spatially proximal Cys amino acids in diverse Ser/Thr protein kinase families suggests that disulfide-mediated control of catalytic activity may be a prevalent mechanism for regulation within the broader AMPK family.

A stress-responsive p38 signaling axis in choanoflagellates

Animal kinases regulate cellular responses to environmental stimuli, including cell differentiation, migration, survival, and response to stress, but the ancestry of these functions is poorly understood. Choanoflagellates, the closest living relatives of animals, encode homologs of diverse animal kinases and have emerged as model organisms for reconstructing animal origins. However, efforts to identify key kinase regulators in choanoflagellates have been constrained by the limitations of currently available genetic tools. Here, we report on a framework that combines small molecule-driven kinase discovery with targeted genetics to reveal kinase function in choanoflagellates. To study the physiological roles of choanoflagellate kinases, we established two high-throughput platforms to screen the model choanoflagellate Salpingoeca rosetta with a curated library of human kinase inhibitors. We identified 95 diverse kinase inhibitors that disrupt S. rosetta cell proliferation. By focusing on one inhibitor, sorafenib, we identified a p38 kinase as a regulator of the heat shock response in S. rosetta. This finding reveals a conserved p38 function between choanoflagellates, animals, and fungi. Moreover, this study demonstrates that existing kinase inhibitors can serve as powerful tools to examine the ancestral roles of kinases that regulate modern animal development.

Untangling the Role of MYC in Sarcomas and Its Potential as a Promising Therapeutic Target

MYC plays a pivotal role in the biology of various sarcoma subtypes, acting as a key regulator of tumor growth, proliferation, and metabolic reprogramming. This oncogene is frequently dysregulated across different sarcomas, where its expression is closely intertwined with the molecular features unique to each subtype. MYC interacts with critical pathways such as cell cycle regulation, apoptosis, and angiogenesis, amplifying tumor aggressiveness and resistance to standard therapies. Furthermore, MYC influences the tumor microenvironment by modulating cell-extracellular matrix interactions and immune evasion mechanisms, further complicating therapeutic management. Despite its well-established centrality in sarcoma pathogenesis, targeting MYC directly remains challenging due to its "undruggable" protein structure. However, emerging therapeutic strategies, including indirect MYC inhibition via epigenetic modulators, transcriptional machinery disruptors, and metabolic pathway inhibitors, offer new hope for sarcoma treatment. This review underscores the importance of understanding the intricate roles of MYC across sarcoma subtypes to guide the development of effective targeted therapies. Given MYC's central role in tumorigenesis and progression, innovative approaches aiming at MYC inhibition could transform the therapeutic landscape for sarcoma patients, providing a much-needed avenue to overcome therapeutic resistance and improve clinical outcomes.

Atlas of the Bacterial Serine-Threonine Kinases expands the functional diversity of the kinome

Bacterial serine-threonine protein kinases (STKs) regulate diverse cellular processes associated with cell growth, virulence, and pathogenicity. They are evolutionarily related to the druggable eukaryotic STKs. However, an incomplete knowledge of how bacterial STKs differ from their eukaryotic counterparts and how they have diverged to regulate diverse bacterial signaling functions presents a bottleneck in targeting them for drug discovery efforts. Here, we classified over 300,000 bacterial STK sequences from the NCBI RefSeq non-redundant and UniProt protein databases into 35 canonical and seven non-canonical (pseudokinase) families based on the patterns of evolutionary constraints in the conserved catalytic domain and flanking regulatory domains. Through statistical comparisons, we identified distinguishing features of bacterial STKs, including a distinctive arginine residue in a regulatory helix (C-Helix) that dynamically couples ATP and substrate binding lobes of the kinase domain. Biochemical and peptide-library screens demonstrated that constrained residues contribute to substrate specificity and kinase activation in the Mycobacterium tuberculosis kinase PknB. Collectively, these findings open new avenues for investigating bacterial STK functions in cellular signaling and for the development of selective bacterial STK inhibitors.

Evolutionary sequence and structural basis for the distinct conformational landscapes of Tyr and Ser/Thr kinases

Protein kinases are molecular machines with rich sequence variation that distinguishes the two main evolutionary branches - tyrosine kinases (TKs) from serine/threonine kinases (STKs). Using a sequence co-variation Potts statistical energy model we previously concluded that TK catalytic domains are more likely than STKs to adopt an inactive conformation with the activation loop in an autoinhibitory folded conformation, due to intrinsic sequence effects. Here we investigate the structural basis for this phenomenon by integrating the sequence-based model with structure-based molecular dynamics (MD) to determine the effects of mutations on the free energy difference between active and inactive conformations, using a thermodynamic cycle involving many (n = 108) protein-mutation free energy perturbation (FEP) simulations in the active and inactive conformations. The sequence and structure-based results are consistent and support the hypothesis that the inactive conformation DFG-out Activation Loop Folded, is a functional regulatory state that has been stabilized in TKs relative to STKs over the course of their evolution via the accumulation of residue substitutions in the activation loop and catalytic loop that facilitate distinct substrate binding modes in trans and additional modes of regulation in cis for TKs.

Evolutionary sequence and structural basis for the distinct conformational landscapes of Tyr and Ser/Thr kinases

Protein kinases are molecular machines with rich sequence variation that distinguishes the two main evolutionary branches - tyrosine kinases (TKs) from serine/threonine kinases (STKs). Using a sequence co-variation Potts statistical energy model we previously concluded that TK catalytic domains are more likely than STKs to adopt an inactive conformation with the activation loop in an autoinhibitory "folded" conformation, due to intrinsic sequence effects. Here we investigated the structural basis for this phenomenon by integrating the sequence-based model with structure-based molecular dynamics (MD) to determine the effects of mutations on the free energy difference between active and inactive conformations, using a novel thermodynamic cycle involving many (n=108) protein-mutation free energy perturbation (FEP) simulations in the active and inactive conformations. The sequence and structure-based results are consistent and support the hypothesis that the inactive conformation "DFG-out Activation Loop Folded", is a functional regulatory state that has been stabilized in TKs relative to STKs over the course of their evolution via the accumulation of residue substitutions in the activation loop and catalytic loop that facilitate distinct substrate binding modes in trans and additional modes of regulation in cis for TKs.

Evolutionary sequence and structural basis for the distinct conformational landscapes of Tyr and Ser/Thr kinases

Protein kinases are molecular machines with rich sequence variation that distinguishes the two main evolutionary branches - tyrosine kinases (TKs) from serine/threonine kinases (STKs). Using a sequence co-variation Potts statistical energy model we previously concluded that TK catalytic domains are more likely than STKs to adopt an inactive conformation with the activation loop in an autoinhibitory "folded" conformation, due to intrinsic sequence effects. Here we investigated the structural basis for this phenomenon by integrating the sequence-based model with structure-based molecular dynamics (MD) to determine the effects of mutations on the free energy difference between active and inactive conformations, using a novel thermodynamic cycle involving many (n=108) protein-mutation free energy perturbation (FEP) simulations in the active and inactive conformations. The sequence and structure-based results are consistent and support the hypothesis that the inactive conformation "DFG-out Activation Loop Folded", is a functional regulatory state that has been stabilized in TKs relative to STKs over the course of their evolution via the accumulation of residue substitutions in the activation loop and catalytic loop that facilitate distinct substrate binding modes in trans and additional modes of regulation in cis for TKs.

The intrinsic substrate specificity of the human tyrosine kinome

Phosphorylation of proteins on tyrosine (Tyr) residues evolved in metazoan organisms as a mechanism of coordinating tissue growth1. Multicellular eukaryotes typically have more than 50 distinct protein Tyr kinases that catalyse the phosphorylation of thousands of Tyr residues throughout the proteome1-3. How a given Tyr kinase can phosphorylate a specific subset of proteins at unique Tyr sites is only partially understood4-7. Here we used combinatorial peptide arrays to profile the substrate sequence specificity of all human Tyr kinases. Globally, the Tyr kinases demonstrate considerable diversity in optimal patterns of residues surrounding the site of phosphorylation, revealing the functional organization of the human Tyr kinome by substrate motif preference. Using this information, Tyr kinases that are most compatible with phosphorylating any Tyr site can be identified. Analysis of mass spectrometry phosphoproteomic datasets using this compendium of kinase specificities accurately identifies specific Tyr kinases that are dysregulated in cells after stimulation with growth factors, treatment with anti-cancer drugs or expression of oncogenic variants. Furthermore, the topology of known Tyr signalling networks naturally emerged from a comparison of the sequence specificities of the Tyr kinases and the SH2 phosphotyrosine (pTyr)-binding domains. Finally we show that the intrinsic substrate specificity of Tyr kinases has remained fundamentally unchanged from worms to humans, suggesting that the fidelity between Tyr kinases and their protein substrate sequences has been maintained across hundreds of millions of years of evolution.

Redox Regulation of Brain Selective Kinases BRSK1/2: Implications for Dynamic Control of the Eukaryotic AMPK family through Cys-based mechanisms

In eukaryotes, protein kinase signaling is regulated by a diverse array of post-translational modifications (PTMs), including phosphorylation of Ser/Thr residues and oxidation of cysteine (Cys) residues. While regulation by activation segment phosphorylation of Ser/Thr residues is well understood, relatively little is known about how oxidation of cysteine residues modulate catalysis. In this study, we investigate redox regulation of the AMPK-related Brain-selective kinases (BRSK) 1 and 2, and detail how broad catalytic activity is directly regulated through reversible oxidation and reduction of evolutionarily conserved Cys residues within the catalytic domain. We show that redox-dependent control of BRSKs is a dynamic and multilayered process involving oxidative modifications of several Cys residues, including the formation of intramolecular disulfide bonds involving a pair of Cys residues near the catalytic HRD motif and a highly conserved T-Loop Cys with a BRSK-specific Cys within an unusual CPE motif at the end of the activation segment. Consistently, mutation of the CPE-Cys increases catalytic activity in vitro and drives phosphorylation of the BRSK substrate Tau in cells. Molecular modeling and molecular dynamics simulations indicate that oxidation of the CPE-Cys destabilizes a conserved salt bridge network critical for allosteric activation. The occurrence of spatially proximal Cys amino acids in diverse Ser/Thr protein kinase families suggests that disulfide mediated control of catalytic activity may be a prevalent mechanism for regulation within the broader AMPK family.

Mechanistic and evolutionary insights into isoform-specific 'supercharging' in DCLK family kinases

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.

Molecular Dynamics and Machine Learning Give Insights on the Flexibility-Activity Relationships in Tyrosine Kinome

Tyrosine kinases are a subfamily of kinases with critical roles in cellular machinery. Dysregulation of their active or inactive forms is associated with diseases like cancer. This study aimed to holistically understand their flexibility-activity relationships, focusing on pockets and fluctuations. We studied 43 different tyrosine kinases by collecting 120 μs of molecular dynamics simulations, pocket and residue fluctuation analysis, and a complementary machine learning approach. We found that the inactive forms often have increased flexibility, particularly at the DFG motif level. Noteworthy, thanks to these long simulations combined with a decision tree, we identified a semiquantitative fluctuation threshold of the DGF+3 residue over which the kinase has a higher probability to be in the inactive form.

Cross-kingdom analyses of transmembrane protein kinases show their functional diversity and distinct origins in protists

Transmembrane kinases (TMKs) are important mediators of cellular signaling cascades. The kinase domains of most metazoan and plant TMKs belong to the serine/threonine/tyrosine kinase (S/T/Y-kinase) superfamily. They share a common origin with prokaryotic kinases and have diversified into distinct subfamilies. Diverse members of the eukaryotic crown radiation such as amoebae, ciliates, and red and brown algae (grouped here under the umbrella term "protists") have long diverged from higher eukaryotes since their ancient common ancestry, making them ideal organisms for studying TMK evolution. Here, we developed an accurate and high-throughput pipeline to predict TMKomes in cellular organisms. Cross-kingdom analyses revealed distinct features of TMKomes in each grouping. Two-transmembrane histidine kinases constitute the main TMKomes of bacteria, while metazoans, plants, and most protists have a large proportion of single-pass TM S/T/Y-kinases. Phylogenetic analyses classified most protist S/T/Y-kinases into three clades, with clades II and III specifically expanded in amoebae and oomycetes, respectively. In contrast, clade I kinases were widespread in all protists examined here, and likely shared a common origin with other eukaryotic S/T/Y-kinases. Functional annotation further showed that most non-kinase domains were grouping-specific, suggesting that their recombination with the more conserved kinase domains led to the divergence of S/T/Y-kinases. However, we also found that protist leucine-rich repeat (LRR)- and G-protein-coupled receptor (GPCR)-type TMKs shared similar sensory domain architectures with respective plant and animal TMKs, despite that they belong to distinct kinase subfamilies. Collectively, our study revealed the functional diversity of TMKomes and the distinct origins of S/T/Y-kinases in protists.

Mechanistic and evolutionary insights into isoform-specific 'supercharging' in DCLK family kinases

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 auto-regulation 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.

Crystal structure of the kinase domain of a receptor tyrosine kinase from a choanoflagellate, Monosiga brevicollis

Genomic analysis of the unicellular choanoflagellate, Monosiga brevicollis (MB), revealed the remarkable presence of cell signaling and adhesion protein domains that are characteristically associated with metazoans. Strikingly, receptor tyrosine kinases, one of the most critical elements of signal transduction and communication in metazoans, are present in choanoflagellates. We determined the crystal structure at 1.95 Å resolution of the kinase domain of the M. brevicollis receptor tyrosine kinase C8 (RTKC8, a member of the choanoflagellate receptor tyrosine kinase C family) bound to the kinase inhibitor staurospaurine. The chonanoflagellate kinase domain is closely related in sequence to mammalian tyrosine kinases (~ 40% sequence identity to the human Ephrin kinase domain EphA3) and, as expected, has the canonical protein kinase fold. The kinase is structurally most similar to human Ephrin (EphA5), even though the extracellular sensor domain is completely different from that of Ephrin. The RTKC8 kinase domain is in an active conformation, with two staurosporine molecules bound to the kinase, one at the active site and another at the peptide-substrate binding site. To our knowledge this is the first example of staurospaurine binding in the Aurora A activation segment (AAS). We also show that the RTKC8 kinase domain can phosphorylate tyrosine residues in peptides from its C-terminal tail segment which is presumably the mechanism by which it transmits the extracellular stimuli to alter cellular function.

Higher-order interactions of Bcr-Abl can broaden chronic myeloid leukemia (CML) drug repertoire

Bcr-Abl, a nonreceptor tyrosine kinase, is associated with leukemias, especially chronic myeloid leukemia (CML). Deletion of Abl's N-terminal region, to which myristoyl is linked, renders the Bcr-Abl fusion oncoprotein constitutively active. The substitution of Abl's N-terminal region by Bcr enables Bcr-Abl oligomerization. Oligomerization is critical: it promotes clustering on the membrane, which is essential for potent MAPK signaling and cell proliferation. Here we decipher the Bcr-Abl specific, step-by-step oligomerization process, identify a specific packing surface, determine exactly how the process is structured and identify its key elements. Bcr's coiled coil (CC) domain at the N-terminal controls Bcr-Abl oligomerization. Crystallography validated oligomerization via Bcr-Abl dimerization between two Bcr CC domains, with tetramerization via tight packing between two binary assemblies. However, the structural principles guiding Bcr CC domain oligomerization are unknown, hindering mechanistic understanding and drugs exploiting it. Using molecular dynamics (MD) simulations, we determine that the binary complex of the Bcr CC domain serves as a basic unit in the quaternary complex providing a specific surface for dimer-dimer packing and higher-order oligomerization. We discover that the small α1-helix is the key. In the binary assembly, the helix forms interchain aromatic dimeric packing, and in the quaternary assembly, it contributes to the specific dimer-dimer packing. Our mechanism is supported by the experimental literature. It offers the key elements controlling this process which can expand the drug discovery strategy, including by Bcr CC-derived peptides, and candidate residues for small covalent drugs, toward quenching oligomerization, supplementing competitive and allosteric tyrosine kinase inhibitors.