Protein tyrosine kinases activate signaling pathways by catalyzing the phosphorylation of tyrosine residues in their substrates. Mounting evidence suggests that, in addition to recognizing phosphorylated tyrosine (pTyr) residues through specific phosphobinding modules, many protein kinases selectively recognize pTyr directly adjacent to the tyrosine residue they phosphorylate and catalyze the formation of twin pTyr-pTyr sites. Here, we demonstrate the importance of this phosphopriming-driven twin pTyr signaling in promoting cell cycle progression through the cell cycle-inhibitory protein p27Kip1. We identify, structurally resolve, and tune two distinct molecular determinants driving the selective recognition of pTyr directly N- and C-terminal to the target phospho-acceptor tyrosine site. We further show structural and biochemical conservation in this recognition, and identify cancer-associated alterations to these determinants that are unable to recognize phosphoprimed substrates. Finally, using an in vivo mouse model of leukemia we show that Bcr-Abl mutants unable to recognize phosphoprimed substrates paradoxically result in enhanced tumor development and progression. These data indicate that Bcr-Abl, like other proto-oncogenes such as Ras or Myc, engages both pro- and anti-oncogenic programs - but in the case of Bcr-Abl, this is accomplished through a mechanism involving traditional and phosphoprimed substrate recognition.

Protein Serine Kinase H1 (PSKH1) was recently identified as a crucial factor in kidney development and is overexpressed in prostate, lung, and kidney cancers. However, little is known about PSKH1 regulatory mechanisms, leading to its classification as a "dark" kinase. Here, we used biochemistry and mass spectrometry to define PSKH1's consensus substrate motif, protein interactors, and how interactors, including Ca2+ sensor proteins, promote or suppress activity. Intriguingly, despite the absence of a canonical Calmodulin binding motif, Ca2+-Calmodulin activated PSKH1 while, in contrast, the ER-resident Ca2+ sensor of the Cab45, Reticulocalbin, Erc55, Calumenin (CREC) family, Reticulocalbin-3, suppressed PSKH1 catalytic activity. In addition to antagonistic regulation of the PSKH1 kinase domain by Ca2+ sensing proteins, we identified UNC119B as a protein interactor that activates PSKH1 via direct engagement of the kinase domain. Our findings identify complementary allosteric mechanisms by which regulatory proteins tune PSKH1's catalytic activity and raise the possibility that different Ca2+ sensors may act more broadly to tune kinase activities by detecting and decoding extremes of intracellular Ca2+ concentrations.

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.

Kinases regulate cellular processes and are essential for understanding cellular function and disease. To investigate the regulatory state of a kinase, numerous methods have been developed to infer kinase activities from phosphoproteomics data using kinase-substrate libraries. However, few phosphorylation sites can be attributed to an upstream kinase in these libraries, limiting the scope of kinase activity inference. Moreover, inferred activities vary across methods, necessitating evaluation for accurate interpretation. Here, we present benchmarKIN, an R package enabling comprehensive evaluation of kinase activity inference methods. Alongside classical perturbation experiments, benchmarKIN introduces a tumor-based benchmarking approach utilizing multi-omics data to identify highly active or inactive kinases. We used benchmarKIN to evaluate kinase-substrate libraries, inference algorithms and the potential of adding predicted kinase-substrate interactions to overcome the coverage limitations. Our evaluation shows most computational methods perform similarly, but the choice of library impacts the inferred activities with a combination of manually curated libraries demonstrating superior performance in recapitulating kinase activities. Additionally, in the tumor-based evaluation, adding predicted targets from NetworKIN further boosts the performance. We then demonstrate how kinase activity inference aids characterize kinase inhibitor responses in cell lines. Overall, benchmarKIN helps researchers to select reliable methods for identifying deregulated kinases.

Hepatic insulin resistance is central to type 2 diabetes (T2D) and metabolic syndrome, but defining the molecular basis of this defect in humans is challenging because of limited tissue access. Utilizing inducible pluripotent stem cells differentiated into hepatocytes from control individuals and patients with T2D and liquid chromatography with tandem mass spectrometry-based (LC-MS/MS-based) phosphoproteomics analysis, we identified a large network of cell-intrinsic alterations in signaling in T2D. Over 300 phosphosites showed impaired or reduced insulin signaling, including losses in the classical insulin-stimulated PI3K/AKT cascade and their downstream targets. In addition, we identified over 500 phosphosites of emergent, i.e., new or enhanced, signaling. These occurred on proteins involved in the Rho-GTPase pathway, RNA metabolism, vesicle trafficking, and chromatin modification. Kinome analysis indicated that the impaired phosphorylation sites represented reduced actions of AKT2/3, PKCθ, CHK2, PHKG2, and/or STK32C kinases. By contrast, the emergent phosphorylation sites were predicted to be mediated by increased action of the Rho-associated kinases 1 and 2 (ROCK1/2), mammalian STE20-like protein kinase 4 (MST4), and/or branched-chain α-ketoacid dehydrogenase kinase (BCKDK). Inhibiting ROCK1/2 activity in T2D induced pluripotent stem cell-derived hepatocytes restored some of the alterations in insulin action. Thus, insulin resistance in the liver in T2D did not simply involve a loss of canonical insulin signaling but the also appearance of new phosphorylations representing a change in the balance of multiple kinases. Together, these led to altered insulin action in the liver and identified important targets for the therapy of hepatic insulin resistance.

Mass spectrometry-based phosphoproteomics offers a comprehensive view of protein phosphorylation, yet our limited knowledge about the regulation and function of most phosphosites hampers the extraction of meaningful biological insights. To address this challenge, we integrate machine learning with phosphoproteomic data from 1195 tumor specimens spanning 11 cancer types to construct CoPheeMap, a network that maps the co-regulation of 26,280 phosphosites. By incorporating network features from CoPheeMap into a second machine learning model, namely CoPheeKSA, we achieve superior performance in predicting kinase-substrate associations. CoPheeKSA uncovers 24,015 associations between 9399 phosphosites and 104 serine/threonine kinases, shedding light on many unannotated phosphosites and understudied kinases. We validate the accuracy of these predictions using experimentally determined kinase-substrate specificities. Through the application of CoPheeMap and CoPheeKSA to phosphosites with high computationally predicted functional significance and those associated with cancer, we demonstrate their effectiveness in systematically elucidating phosphosites of interest. These analyses unveil dysregulated signaling processes in human cancer and identify understudied kinases as potential therapeutic targets.

Chromosome instability (CIN) is a prevalent hallmark of aggressive cancers that is largely absent in normal cells. Targeting tumors with CIN is a therapeutic strategy that is rapidly advancing. Volastra Therapeutics has undertaken multiple orthogonal approaches to target CIN as a cancer vulnerability. One such approach that relies on the Cancer Dependency Map (DepMap) identified BUB1B (BUB1 Mitotic Checkpoint Serine/Threonine Kinase B) as a top scoring gene required for survival of cell lines with increased genomic instability. BUB1B, a core component of the spindle assembly checkpoint that functions in cell cycle progression and chromosome segregation, was also highlighted in published screens as a major cancer vulnerability in cells with aneuploidy and whole-genome duplication. To corroborate the dependency of chromosomally unstable cancer cells on BUB1B, we used three independent cell models: 1) an isogenic CIN system based on modulating kinetochore-microtubule dynamics, 2) a Cas9-engineered system to induce dicentric chromosomes through telomere deprotection and 3) a non-isogenic cell line model based on differential levels of aneuploidy and micronuclei. Data in these three orthogonal systems confirmed that chromosomally unstable cells rely heavily on BUB1B for survival. We next explored the development of BUB1B small molecule inhibitors. The C-terminus of BUB1B is responsible for chromosome congression and contains the putative kinase domain, suggesting that catalytic activity may be responsible for CIN suppression. Whether BUB1B possesses enzymatic activity has been controversial, with some groups reporting low but inhibitable BUB1B kinase activity, while others indicate that BUB1B lacks key kinase residues and is catalytically inactive. To probe the function of BUB1B, we generated full-length and kinase domain recombinant proteins. Both constructs showed very weak ATPase activity comparable to published and commercially available proteins. Efforts to boost the kinase activity through improved substrate specificity or enhanced phosphorylation were not successful, and additional purification resulted in further loss of kinase activity. Furthermore, modeling of the BUB1B kinase domain suggests that ATP binding to the ATP cleft is physically hindered by bulky residues. Finally, substrate sequence mapping did not match sequences predicted for BUB1B. Together the data suggest that BUB1B is not enzymatically active and that the measured low ATPase activity may be attributable to a co-purifying kinase. Our work confirms chromosomally unstable cells have greater dependency on BUB1B. However, the protein lacks catalytic activity and thus cannot be targeted with a kinase inhibitor strategy. BUB1B’s function as a protein interaction hub with specific partners could permit alternative therapeutic strategies that may result in effective treatments for patients with chromosomally unstable cancers. Nazario Bosco, Rafael S. Depetris, Jasmina Prabhakara, Akanksha Verma, Jared L. Johnson, Lewis C. Cantley, Scott Drutman, Shinsan M. Su, Rumin Zhang, Derek A. Cogan, Christina H. Eng. BUB1B is an inactive kinase required for survival of chromosomally unstable cells [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2025; Part 1 (Regular Abstracts); 2025 Apr 25-30; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2025;85(8_Suppl_1):Abstract nr 2983.

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.

Despite the successes of immunotherapy in cancer treatment over recent decades, less than <10%-20% cancer cases have demonstrated durable responses from immune checkpoint blockade. To enhance the efficacy of immunotherapies, combination therapies suppressing multiple immune evasion mechanisms are increasingly contemplated. To better understand immune cell surveillance and diverse immune evasion responses in tumor tissues, we comprehensively characterized the immune landscape of more than 1,000 tumors across ten different cancers using CPTAC pan -cancer proteogenomic data. We identified seven distinct immune subtypes based on integrative learning of cell type compositions and pathway activities. We then thoroughly categorized unique genomic, epigenetic, transcriptomic, and proteomic changes associated with each subtype. Further leveraging the deep phosphoproteomic data, we studied kinase activities in different immune subtypes, which revealed potential subtype -specific therapeutic targets. Insights from this work will facilitate the development of future immunotherapy strategies and enhance precision targeting with existing agents.

The peptide hormone glucagon is a fundamental metabolic regulator that is also being considered as a pharmacotherapeutic option for obesity and type 2 diabetes. Despite this, we know very little regarding how glucagon exerts its pleiotropic metabolic actions. Given that the liver is a chief site of action, we performed in situ time-resolved liver phosphoproteomics to reveal glucagon signaling nodes. Through pathway analysis of the thousands of phosphopeptides identified, we reveal "membrane trafficking" as a dominant signature with the vesicle trafficking protein SEC22 Homolog B (SEC22B) S137 phosphorylation being a top hit. Hepatocyte-specific loss- and gain-of-function experiments reveal that SEC22B was a key regulator of glycogen, lipid and amino acid metabolism, with SEC22B-S137 phosphorylation playing a major role in glucagon action. Mechanistically, we identify several protein binding partners of SEC22B affected by glucagon, some of which were differentially enriched with SEC22B-S137 phosphorylation. In summary, we demonstrate that phosphorylation of SEC22B is a hepatocellular signaling node mediating the metabolic actions of glucagon and provide a rich resource for future investigations on the biology of glucagon action.

CK1 enzymes phosphorylate substrates that participate in a variety of signaling pathways, yet the mechanisms by which these kinases are regulated remain elusive. Human and yeast CK1s autophosphorylate their C-terminal non-catalytic tails, and C-terminal truncation or dephosphorylation of these enzymes has been found to increase substrate phosphorylation in vitro. This has led to a model in which the autophosphorylated C-terminus acts as a pseudosubstrate to inhibit catalytic activity. Here, we have directly evaluated this proposed mechanism in vitro and in vivo using the CK1 homologues Schizosaccharomyces pombe Hhp1 and human CK1 delta and CK1 epsilon. First, we identified the C-terminal autophosphorylation sites, confirming 6 on Hhp1 and 15 on CK1 epsilon. In vivo, Hhp1 autophosphorylated its C-terminus on the sites we identified, and the occupancy of these sites varied under different cellular conditions, indicating that autophosphorylation is dynamic in cells. In vitro, peptides corresponding to the Hhp1 and CK1 epsilon C-termini directly interacted with the kinase domains, but only when these sites were phosphorylated. Mutating these sites to non-phosphorylatable residues increased the activity of both enzymes towards their substrates. Interestingly, we found that substrates have a higher affinity for the kinase domain than the autophosphorylated C-termini, and substrates competitively inhibited tail binding even when the tail was fully autophosphorylated. Mutating the basic residues on the kinase domain that comprise the substrate binding groove prevented both tail binding and substrate phosphorylation, supporting the model that tail and substrate compete for binding to the substrate binding groove. Furthermore, we found that the C-termini contribute to substrate specificity, in combination with the previously described autophosphorylation site at T220 that resides in the catalytic domain. Truncation of the CK1 delta C-terminus relaxed the kinase's preference for an upstream phosphoprimed residue in the linear substrate motif, allowing additional specificity determinants to be observed. Autophosphorylation at T220 affected the likelihood of tail autophosphorylation and binding, demonstrating that these two types of autophosphorylation sites work in concert. Both C-terminal and T220 autophosphorylation affected the specificity constants measured in vitro for different protein substrates of yeast and human CK1s. Based on these data, we propose a displacement specificity model to explain how CK1 autophosphorylation may regulate phosphorylation of specific substrates in different pathways and under different conditions. This work was funded by the National Institutes of Health, R35-GM131799 to KLG and T32-CA119925 to SNC.

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.

PTEN plays a crucial role in preventing the development of glioblastoma (GBM), a severe and untreatable brain cancer. In GBM, most PTEN deficiencies are missense mutations that have not been thoroughly examined. Here, we leveraged genetically modified mice and isogenic astrocyte cell cultures to investigate the role of clinically relevant mutations (G36E, L42R, C105F, and R173H) in the development of EGFR-driven GBM. We report that the loss of tumor suppression from these mutants is unrelated to their lipid phosphatase activity and rather relate to elevated localization at the cell membrane. Moreover, expression of these PTEN mutations heightened EGFR activity by sequestering EGFR within endomembranes longer and affected its signaling behavior. Through comprehensive studies on global protein phosphorylation and kinase library analyses in cells with the G36E and L42R PTEN mutations, we identified distinct cancer-promoting pathways activated by EGFR, offering targets for treating GBM with these PTEN alterations.

Dysregulation of intracellular signaling networks underpins cancer. However, a systems-level elucidation of how signaling networks within distinct cell subpopulations drive cancer progression in vivo has been unattainable due to technical limitations. We developed INSIGHT (INvestigating SIGnaling network of specific cell subpopulation in Heterogeneous Tissue), a new platform technology combining fluorescence-activated cell sorting with ultra-sensitive mass spectrometry to enable phosphoproteomic characterization of rare and discrete cell subpopulations from fixed tissues. We demonstrated the broad utility of INSIGHT by analyzing the oligodendroglial cell-specific signaling network in the mouse brain. We then applied INSIGHT to investigate the rare, disseminated tumor cell subpopulation in glioblastoma patient-derived xenograft models. INSIGHT uncovered a global rewiring of signaling networks with tumor cell dissemination, marked by a transition from proliferation-associated signaling in the primary tumor cells to signaling associated with postsynapse, neuronal migration, and ion homeostasis in disseminated tumor cells. We reveal interconnections between signaling circuitries within the networks, with numerous proteins, including GluA2, exhibiting altered phosphorylation without protein expression changes, emphasizing the role of post-translational modifications in glioblastoma dissemination. We validated key phosphorylation changes and inferred differentially active kinases with tumor spread to offer new systems-level insights into glioblastoma dissemination mechanisms in vivo. INSIGHT is generally applicable to a wide range of biological systems without genetic engineering and provides quantitative phosphorylation and protein expression data for selected cell subpopulations from heterogeneous tissues. ### Competing Interest Statement The authors have declared no competing interest.

While ultraviolet (UV) radiation damages DNA, eliciting the DNA damage response (DDR), it also damages RNA, triggering transcriptome-wide ribosomal collisions and eliciting a ribotoxic stress response (RSR). However, the relative contributions, timing, and regulation of these pathways in determining cell fate is unclear. Here we use time-resolved phosphoproteomic, chemical-genetic, single-cell imaging, and biochemical approaches to create a chronological atlas of signaling events activated in cells responding to UV damage. We discover that UV-induced apoptosis is mediated by the RSR kinase ZAK and not through the DDR. We identify two negative-feedback modules that regulate ZAK-mediated apoptosis: (1) GCN2 activation limits ribosomal collisions and attenuates ZAK-mediated RSR and (2) ZAK activity leads to phosphodegron autophosphorylation and its subsequent degradation. These events tune ZAK's activity to collision levels to establish regimes of homeostasis, tolerance, and death, revealing its key role as the cellular sentinel for nucleic acid damage.

T cell activation is a complex biological process of naive cells maturing into effector cells. Proteomic and phospho-proteomic approaches have provided critical insights into this process, yet it is not always clear how changes in individual proteins or phosphorylation sites have functional significance. Here, we developed the Phosphorylation Integrated Thermal Shift Assay (PITSA) that combines the measurement of protein or phosphorylation site abundance and thermal stability into a single tandem mass tags experiment and apply this method to study T cell activation. We quantified the abundance and thermal stability of over 7500 proteins and 5000 phosphorylation sites and identified significant differences in chromatin-related, TCR signaling, DNA repair, and proliferative phosphoproteins. PITSA may be applied to a wide range of biological contexts to generate hypotheses as to which proteins or phosphorylation sites are functionally regulated in a given system as well as the mechanisms by which this regulation may occur.

Although genomic anomalies in glioblastoma (GBM) have been well studied for over a decade, its 5-year survival rate remains lower than 5%. We seek to expand the molecular landscape of high-grade glioma, composed of IDH-wildtype GBM and IDH-mutant grade 4 astrocytoma, by integrating proteomic, metabolomic, lipidomic, and post-translational modifications (PTMs) with genomic and transcriptomic measurements to uncover multi-scale regulatory interactions governing tumor development and evolution. Applying 14 proteogenomic and metabolomic platforms to 228 tumors (212 GBM and 16 grade 4 IDH-mutant astrocytoma), including 28 at recurrence, plus 18 normal brain samples and 14 brain metastases as comparators, reveals heterogeneous upstream alterations converging on common downstream events at the proteomic and metabolomic levels and changes in protein-protein interactions and glycosylation site occupancy at recurrence. Recurrent genetic alterations and phosphorylation events on PTPN11 map to important regulatory domains in three dimensions, suggesting a central role for PTPN11 signaling across high-grade gliomas.