Coral: Clear and Customizable Visualization of Human Kinome Data

Structure-Based Drug Design of 2-Amino-[1,1'-biphenyl]-3-carboxamide Derivatives as Selective PKMYT1 Inhibitors for the Treatment of CCNE1-Amplified Breast Cancer

CCNE1 amplification occurs in breast cancer and currently lacks effective therapies. PKMYT1 as a synthetic lethal target for CCNE1 amplification holds promise for the treatment of CCNE1-amplified breast cancer. Herein, we discover a series of 2-amino-[1,1'-biphenyl]-3-carboxamide derivatives as potent and selective PKMYT1 inhibitors using structure-based drug design. The representative compound 8ma exhibited excellent potency against PKMYT1, while sparing WEE1. It also suppressed proliferation of the CCNE1-amplified HCC1569 breast cancer cell line and showed synergistic cytotoxicity in combination with gemcitabine. PKMYT1 X-ray cocrystallography confirmed that introduction of key binding interactions between the inhibitors and residues Asp251 and Tyr121 of PKMYT1 greatly enhanced the potency and selectivity of the compounds.

PTMoreR-enabled cross-species PTM mapping and comparative phosphoproteomics across mammals

To support PTM proteomic analysis and annotation in different species, we developed PTMoreR, a user-friendly tool that considers the surrounding amino acid sequences of PTM sites during BLAST, enabling a motif-centric analysis across species. By controlling sequence window similarity, PTMoreR can map phosphoproteomic results between any two species, perform site-level functional enrichment analysis, and generate kinase-substrate networks. We demonstrate that the majority of real P-sites in mice can be inferred from experimentally derived human P-sites with PTMoreR mapping. Furthermore, the compositions of 129 mammalian phosphoproteomes can also be predicted using PTMoreR. The method also identifies cross-species phosphorylation events that occur on proteins with an increased tendency to respond to the environmental factors. Moreover, the classic kinase motifs can be extracted across mammalian species, offering an evolutionary angle for refining current motifs. PTMoreR supports PTM proteomics in non-human species and facilitates quantitative phosphoproteomic analysis.

Strategic Fluorination to Achieve a Potent, Selective, Metabolically Stable, and Orally Bioavailable Inhibitor of CSNK2

The host kinase casein kinase 2 (CSNK2) has been proposed to be an antiviral target against β-coronaviral infection. To pharmacologically validate CSNK2 as a drug target in vivo, potent and selective CSNK2 inhibitors with good pharmacokinetic properties are required. Inhibitors based on the pyrazolo[1,5-a]pyrimidine scaffold possess outstanding potency and selectivity for CSNK2, but bioavailability and metabolic stability are often challenging. By strategically installing a fluorine atom on an electron-rich phenyl ring of a previously characterized inhibitor 1, we discovered compound 2 as a promising lead compound with improved in vivo metabolic stability. Compound 2 maintained excellent cellular potency against CSNK2, submicromolar antiviral potency, and favorable solubility, and was remarkably selective for CSNK2 when screened against 192 kinases across the human kinome. We additionally present a co-crystal structure to support its on-target binding mode. In vivo, compound 2 was orally bioavailable, and demonstrated modest and transient inhibition of CSNK2, although antiviral activity was not observed, possibly attributed to its lack of prolonged CSNK2 inhibition.

Coordinated Targeting of S6K1/2 and AXL Disrupts Pyrimidine Biosynthesis in PTEN-Deficient Glioblastoma

Intrinsic resistance to targeted therapeutics in PTEN-deficient glioblastoma (GBM) is mediated by redundant signaling networks that sustain critical metabolic functions. Here, we demonstrate that coordinated inhibition of the ribosomal protein S6 kinase 1 (S6K1) and the receptor tyrosine kinase AXL using LY-2584702 and BMS-777607 can overcome network redundancy to reduce GBM tumor growth. This combination of S6K1 and AXL inhibition suppressed glucose flux to pyrimidine biosynthesis. Genetic inactivation studies to map the signaling network indicated that both S6K1 and S6K2 transmit growth signals in PTEN-deficient GBM. Kinome-wide ATP binding analysis in inhibitor-treated cells revealed that LY-2584702 directly inhibited S6K1, and substrate phosphorylation studies showed that BMS-777607 inactivation of upstream AXL collaborated to reduce S6K2-mediated signal transduction. Thus, combination targeting of S6K1 and AXL provides a kinase-directed therapeutic approach that circumvents signal transduction redundancy to interrupt metabolic function and reduce growth of PTEN-deficient GBM.

SIGNIFICANCE

Therapy for glioblastoma would be advanced by incorporating molecularly targeted kinase-directed agents, similar to standard of care strategies in other tumor types. Here, we identify a kinase targeting approach to inhibit the metabolism and growth of glioblastoma.

A first-in-class selective inhibitor of EGFR and PI3K offers a single-molecule approach to targeting adaptive resistance

Despite tremendous progress in precision oncology, adaptive resistance mechanisms limit the long-term effectiveness of molecularly targeted agents. Here we evaluated the pharmacological profile of MTX-531 that was computationally designed to selectively target two key resistance drivers, epidermal growth factor receptor and phosphatidylinositol 3-OH kinase (PI3K). MTX-531 exhibits low-nanomolar potency against both targets with a high degree of specificity predicted by cocrystal structural analyses. MTX-531 monotherapy uniformly resulted in tumor regressions of squamous head and neck patient-derived xenograft (PDX) models. The combination of MTX-531 with mitogen-activated protein kinase kinase or KRAS-G12C inhibitors led to durable regressions of BRAF-mutant or KRAS-mutant colorectal cancer PDX models, resulting in striking increases in median survival. MTX-531 is exceptionally well tolerated in mice and uniquely does not lead to the hyperglycemia commonly seen with PI3K inhibitors. Here, we show that MTX-531 acts as a weak agonist of peroxisome proliferator-activated receptor-γ, an attribute that likely mitigates hyperglycemia induced by PI3K inhibition. This unique feature of MTX-531 confers a favorable therapeutic index not typically seen with PI3K inhibitors.

More than an Amide Bioisostere: Discovery of 1,2,4-Triazole-containing Pyrazolo[1,5-a]pyrimidine Host CSNK2 Inhibitors for Combatting β-Coronavirus Replication

The pyrazolo[1,5-a]pyrimidine scaffold is a promising scaffold to develop potent and selective CSNK2 inhibitors with antiviral activity against β-coronaviruses. Herein, we describe the discovery of a 1,2,4-triazole group to substitute a key amide group for CSNK2 binding present in many potent pyrazolo[1,5-a]pyrimidine inhibitors. Crystallographic evidence demonstrates that the 1,2,4-triazole replaces the amide in forming key hydrogen bonds with Lys68 and a water molecule buried in the ATP-binding pocket. This isosteric replacement improves potency and metabolic stability at a cost of solubility. Optimization for potency, solubility, and metabolic stability led to the discovery of the potent and selective CSNK2 inhibitor 53. Despite excellent in vitro metabolic stability, rapid decline in plasma concentration of 53 in vivo was observed and may be attributed to lung accumulation, although in vivo pharmacological effect was not observed. Further optimization of this novel chemotype may validate CSNK2 as an antiviral target in vivo.

Memory effects of prior subculture may impact the quality of multiomic perturbation profiles

Mass spectrometry-based omics technologies are increasingly used in perturbation studies to map drug effects to biological pathways by identifying significant molecular events. Significance is influenced by fold change and variation of each molecular parameter, but also by multiple testing corrections. While the fold change is largely determined by the biological system, the variation is determined by experimental workflows. Here, it is shown that memory effects of prior subculture can influence the variation of perturbation profiles using the two colon carcinoma cell lines SW480 and HCT116. These memory effects are largely driven by differences in growth states that persist into the perturbation experiment. In SW480 cells, memory effects combined with moderate treatment effects amplify the variation in multiple omics levels, including eicosadomics, proteomics, and phosphoproteomics. With stronger treatment effects, the memory effect was less pronounced, as demonstrated in HCT116 cells. Subculture homogeneity was controlled by real-time monitoring of cell growth. Controlled homogeneous subculture resulted in a perturbation network of 321 causal conjectures based on combined proteomic and phosphoproteomic data, compared to only 58 causal conjectures without controlling subculture homogeneity in SW480 cells. Some cellular responses and regulatory events were identified that extend the mode of action of arsenic trioxide (ATO) only when accounting for these memory effects. Controlled prior subculture led to the finding of a synergistic combination treatment of ATO with the thioredoxin reductase 1 inhibitor auranofin, which may prove useful in the management of NRF2-mediated resistance mechanisms.

Ikaros sets the threshold for negative B-cell selection by regulation of the signaling strength of the AKT pathway

Inhibitory phosphatases, such as the inositol-5-phosphatase SHIP1 could potentially contribute to B-cell acute lymphoblastic leukemia (B-ALL) by raising the threshold for activation of the autoimmunity checkpoint, allowing malignant cells with strong oncogenic B-cell receptor signaling to escape negative selection. Here, we show that SHIP1 is differentially expressed across B-ALL subtypes and that high versus low SHIP1 expression is associated with specific B-ALL subgroups. In particular, we found high SHIP1 expression in both, Philadelphia chromosome (Ph)-positive and ETV6-RUNX1-rearranged B-ALL cells. As demonstrated by targeted knockdown of SHIP1 by RNA interference, proliferation of B-ALL cells in vitro and their tumorigenic spread in vivo depended in part on SHIP1 expression. We investigated the regulation of SHIP1, as an important antagonist of the AKT signaling pathway, by the B-cell-specific transcription factor Ikaros. Targeted restoration of Ikaros and pharmacological inhibition of the antagonistic casein kinase 2, led to a strong reduction in SHIP1 expression and at the same time to a significant inhibition of AKT activation and cell growth. Importantly, the tumor suppressive function of Ikaros was enhanced by a SHIP1-dependent additive effect. Furthermore, our study shows that all three AKT isoforms contribute to the pro-mitogenic and anti-apoptotic signaling in B-ALL cells. Conversely, hyperactivation of a single AKT isoform is sufficient to induce negative selection by increased oxidative stress. In summary, our study demonstrates the regulatory function of Ikaros on SHIP1 expression in B-ALL and highlights the relevance of sustained SHIP1 expression to prevent cells with hyperactivated PI3K/AKT/mTOR signaling from undergoing negative selection.

Proteomic Changes Induced by the Immunosuppressant Everolimus in Human Podocytes

mTOR inhibitors (mTOR-Is) may induce proteinuria in kidney transplant recipients through podocyte damage. However, the mechanism has only been partially defined. Total cell lysates and supernatants of immortalized human podocytes treated with different doses of everolimus (EVE) (10, 100, 200, and 500 nM) for 24 h were subjected to mass spectrometry-based proteomics. Support vector machine and partial least squares discriminant analysis were used for data analysis. The results were validated in urine samples from 28 kidney transplant recipients receiving EVE as part of their immunosuppressive therapy. We identified more than 7000 differentially expressed proteins involved in several pathways, including kinases, cell cycle regulation, epithelial-mesenchymal transition, and protein synthesis, according to gene ontology. Among these, after statistical analysis, 65 showed an expression level significantly and directly correlated with EVE dosage. Polo-Like Kinase 1 (PLK1) content was increased, whereas osteopontin (SPP1) content was reduced in podocytes and supernatants in a dose-dependent manner and significantly correlated with EVE dose (p < 0.0001, FDR < 5%). Similar results were obtained in the urine of kidney transplant patients. This study analyzed the impact of different doses of mTOR-Is on podocytes, helping to understand not only the biological basis of their therapeutic effects but also the possible mechanisms underlying proteinuria.

Determining the ERK-regulated phosphoproteome driving KRAS-mutant cancer

To delineate the mechanisms by which the ERK1 and ERK2 mitogen-activated protein kinases support mutant KRAS-driven cancer growth, we determined the ERK-dependent phosphoproteome in KRAS-mutant pancreatic cancer. We determined that ERK1 and ERK2 share near-identical signaling and transforming outputs and that the KRAS-regulated phosphoproteome is driven nearly completely by ERK. We identified 4666 ERK-dependent phosphosites on 2123 proteins, of which 79 and 66%, respectively, were not previously associated with ERK, substantially expanding the depth and breadth of ERK-dependent phosphorylation events and revealing a considerably more complex function for ERK in cancer. We established that ERK controls a highly dynamic and complex phosphoproteome that converges on cyclin-dependent kinase regulation and RAS homolog guanosine triphosphatase function (RHO GTPase). Our findings establish the most comprehensive molecular portrait and mechanisms by which ERK drives KRAS-dependent pancreatic cancer growth.

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.

DCLK1-Mediated Regulation of Invadopodia Dynamics and Matrix Metalloproteinase Trafficking Drives Invasive Progression in Head and Neck Squamous Cell Carcinoma

Head and neck squamous cell carcinoma (HNSCC) is a major health concern due to its high mortality from poor treatment responses and locoregional tumor invasion into life sustaining structures in the head and neck. A deeper comprehension of HNSCC invasion mechanisms holds the potential to inform targeted therapies that may enhance patient survival. We previously reported that doublecortin like kinase 1 (DCLK1) regulates invasion of HNSCC cells. Here, we tested the hypothesis that DCLK1 regulates proteins within invadopodia to facilitate HNSCC invasion. Invadopodia are specialized subcellular protrusions secreting matrix metalloproteinases that degrade the extracellular matrix (ECM). Through a comprehensive proteome analysis comparing DCLK1 control and shDCLK1 conditions, our findings reveal that DCLK1 plays a pivotal role in regulating proteins that orchestrate cytoskeletal and ECM remodeling, contributing to cell invasion. Further, we demonstrate in TCGA datasets that DCLK1 levels correlate with increasing histological grade and lymph node metastasis. We identified higher expression of DCLK1 in the leading edge of HNSCC tissue. Knockdown of DCLK1 in HNSCC reduced the number of invadopodia, cell adhesion and colony formation. Using super resolution microscopy, we demonstrate localization of DCLK1 in invadopodia and colocalization with mature invadopodia markers TKS4, TKS5, cortactin and MT1-MMP. We carried out phosphoproteomics and validated using immunofluorescence and proximity ligation assays, the interaction between DCLK1 and motor protein KIF16B. Pharmacological inhibition or knockdown of DCLK1 reduced interaction with KIF16B, secretion of MMPs, and cell invasion. This research unveils a novel function of DCLK1 within invadopodia to regulate the trafficking of matrix degrading cargo. The work highlights the impact of targeting DCLK1 to inhibit locoregional invasion, a life-threatening attribute of HNSCC.

Mevalonate kinase-deficient THP-1 cells show a disease-characteristic pro-inflammatory phenotype

Objective

Bi-allelic pathogenic variants in the MVK gene, which encodes mevalonate kinase (MK), an essential enzyme in isoprenoid biosynthesis, cause the autoinflammatory metabolic disorder mevalonate kinase deficiency (MKD). We generated and characterized MK-deficient monocytic THP-1 cells to identify molecular and cellular mechanisms that contribute to the pro-inflammatory phenotype of MKD.

Methods

Using CRISPR/Cas9 genome editing, we generated THP-1 cells with different MK deficiencies mimicking the severe (MKD-MA) and mild end (MKD-HIDS) of the MKD disease spectrum. Following confirmation of previously established disease-specific biochemical hallmarks, we studied the consequences of the different MK deficiencies on LPS-stimulated cytokine release, glycolysis versus oxidative phosphorylation rates, cellular chemotaxis and protein kinase activity.

Results

Similar to MKD patients' cells, MK deficiency in the THP-1 cells caused a pro-inflammatory phenotype with a severity correlating with the residual MK protein levels. In the MKD-MA THP-1 cells, MK protein levels were barely detectable, which affected protein prenylation and was accompanied by a profound pro-inflammatory phenotype. This included a markedly increased LPS-stimulated release of pro-inflammatory cytokines and a metabolic switch from oxidative phosphorylation towards glycolysis. We also observed increased activity of protein kinases that are involved in cell migration and proliferation, and in innate and adaptive immune responses. The MKD-HIDS THP-1 cells had approximately 20% residual MK activity and showed a milder phenotype, which manifested mainly upon LPS stimulation or exposure to elevated temperatures.

Conclusion

MK-deficient THP-1 cells show the biochemical and pro-inflammatory phenotype of MKD and are a good model to study underlying disease mechanisms and therapeutic options of this autoinflammatory disorder.

Development of Selective Pyrido[2,3-d]pyrimidin-7(8H)-one-Based Mammalian STE20-Like (MST3/4) Kinase Inhibitors

Mammalian STE20-like (MST) kinases 1-4 play key roles in regulating the Hippo and autophagy pathways, and their dysregulation has been implicated in cancer development. In contrast to the well-studied MST1/2, the roles of MST3/4 are less clear, in part due to the lack of potent and selective inhibitors. Here, we re-evaluated literature compounds, and used structure-guided design to optimize the p21-activated kinase (PAK) inhibitor G-5555 (8) to selectively target MST3/4. These efforts resulted in the development of MR24 (24) and MR30 (27) with good kinome-wide selectivity and high cellular potency. The distinct cellular functions of closely related MST kinases can now be elucidated with subfamily-selective chemical tool compounds using a combination of the MST1/2 inhibitor PF-06447475 (2) and the two MST3/4 inhibitors developed. We found that MST3/4-selective inhibition caused a cell-cycle arrest in the G1 phase, whereas MST1/2 inhibition resulted in accumulation of cells in the G2/M phase.

Illumination of understudied ciliary kinases

Cilia are cellular signaling hubs. Given that human kinases are central regulators of signaling, it is not surprising that kinases are key players in cilia biology. In fact, many kinases modulate ciliogenesis, which is the generation of cilia, and distinct ciliary pathways. Several of these kinases are understudied with few publications dedicated to the interrogation of their function. Recent efforts to develop chemical probes for members of the cyclin-dependent kinase like (CDKL), never in mitosis gene A (NIMA) related kinase (NEK), and tau tubulin kinase (TTBK) families either have delivered or are working toward delivery of high-quality chemical tools to characterize the roles that specific kinases play in ciliary processes. A better understanding of ciliary kinases may shed light on whether modulation of these targets will slow or halt disease onset or progression. For example, both understudied human kinases and some that are more well-studied play important ciliary roles in neurons and have been implicated in neurodevelopmental, neurodegenerative, and other neurological diseases. Similarly, subsets of human ciliary kinases are associated with cancer and oncological pathways. Finally, a group of genetic disorders characterized by defects in cilia called ciliopathies have associated gene mutations that impact kinase activity and function. This review highlights both progress related to the understanding of ciliary kinases as well as in chemical inhibitor development for a subset of these kinases. We emphasize known roles of ciliary kinases in diseases of the brain and malignancies and focus on a subset of poorly characterized kinases that regulate ciliary biology.

Illuminating function of the understudied druggable kinome

The human kinome, with more than 500 proteins, is crucial for cell signaling and disease. Yet, about one-third of kinases lack in-depth study. The Data and Resource Generating Center for Understudied Kinases has developed multiple resources to address this challenge including creation of a heavy amino acid peptide library for parallel reaction monitoring and quantitation of protein kinase expression, use of understudied kinases tagged with a miniTurbo-biotin ligase to determine interaction networks by proximity-dependent protein biotinylation, NanoBRET probe development for screening chemical tool target specificity in live cells, characterization of small molecule chemical tools inhibiting understudied kinases, and computational tools for defining kinome architecture. These resources are available through the Dark Kinase Knowledgebase, supporting further research into these understudied protein kinases.

Multiplex single-cell chemical genomics reveals the kinase dependence of the response to targeted therapy

Chemical genetic screens are a powerful tool for exploring how cancer cells' response to drugs is shaped by their mutations, yet they lack a molecular view of the contribution of individual genes to the response to exposure. Here, we present sci-Plex-Gene-by-Environment (sci-Plex-GxE), a platform for combined single-cell genetic and chemical screening at scale. We highlight the advantages of large-scale, unbiased screening by defining the contribution of each of 522 human kinases to the response of glioblastoma to different drugs designed to abrogate signaling from the receptor tyrosine kinase pathway. In total, we probed 14,121 gene-by-environment combinations across 1,052,205 single-cell transcriptomes. We identify an expression signature characteristic of compensatory adaptive signaling regulated in a MEK/MAPK-dependent manner. Further analyses aimed at preventing adaptation revealed promising combination therapies, including dual MEK and CDC7/CDK9 or nuclear factor κB (NF-κB) inhibitors, as potent means of preventing transcriptional adaptation of glioblastoma to targeted therapy.

Glucose restriction enhances oxidative fiber formation: A multi-omic signal network involving AMPK and CaMK2

Skeletal muscle is a highly plastic organ that adapts to different metabolic states or functional demands. This study explored the impact of permanent glucose restriction (GR) on skeletal muscle composition and metabolism. Using Glut4m mice with defective glucose transporter 4, we conducted multi-omics analyses at different ages and after low-intensity treadmill training. The oxidative fibers were significantly increased in Glut4m muscles. Mechanistically, GR activated AMPK pathway, promoting mitochondrial function and beneficial myokine expression, and facilitated slow fiber formation via CaMK2 pathway. Phosphorylation-activated Perm1 may synergize AMPK and CaMK2 signaling. Besides, MAPK and CDK kinases were also implicated in skeletal muscle protein phosphorylation during GR response. This study provides a comprehensive signaling network demonstrating how GR influences muscle fiber types and metabolic patterns. These insights offer valuable data for understanding oxidative fiber formation mechanisms and identifying clinical targets for metabolic diseases.

Discovery of pyrazolopyrimidines that selectively inhibit CSF-1R kinase by iterative design, synthesis and screening against glioblastoma cells

Glioblastoma multiforme (GBM) is the most aggressive type of brain cancer in adults, with an average life expectancy under treatment of approx. 15 months. GBM is characterised by a complex set of genetic alterations that results in significant disruption of receptor tyrosine kinase (RTK) signaling. We report here an exploration of the pyrazolo[3,4-d]pyrimidine scaffold in search for antiproliferative compounds directed to GBM treatment. Small compound libraries were synthesised and screened against GBM cells to build up structure-antiproliferative activity-relationships (SAARs) and inform further rounds of design, synthesis and screening. 76 novel compounds were generated through this iterative process that found low micromolar potencies against selected GBM lines, including patient-derived stem cells. Phenomics analysis demonstrated preferential activity against glioma cells of the mesenchymal subtype, whereas kinome screening identified colony stimulating factor-1 receptor (CSF-1R) as the lead's target, a RTK implicated in the tumourigenesis and progression of different cancers and the immunoregulation of the GBM microenvironment.

Identification of Central Nervous System Oncologic Disease Biomarkers in EVs from Cerebrospinal Fluid (CSF) of Pediatric Patients: A Pilot Neuro-Proteomic Study

Cerebrospinal fluid (CSF) is a biochemical-clinical window into the brain. Unfortunately, its wide dynamic range, low protein concentration, and small sample quantity significantly limit the possibility of using it routinely. Extraventricular drainage (EVD) of CSF allows us to solve quantitative problems and to study the biological role of extracellular vesicles (EVs). In this study, we implemented bioinformatic analysis of our previous data of EVD of CSF and its EVs obtained from congenital hydrocephalus with the aim of identifying a comprehensive list of potential tumor and non-tumor biomarkers of central nervous system diseases. Among all proteins identified, those enriched in EVs are associated with synapses, synaptosomes, and nervous system diseases including gliomas, embryonal tumors, and epilepsy. Among these EV-enriched proteins, given the broad consensus present in the recent scientific literature, we validated syntaxin-binding protein 1 (STXBP1) as a marker of malignancy in EVD of CSF and its EVs from patients with pilocytic astrocytoma and medulloblastoma. Our results show that STXBP1 is negatively enriched in EVs compared to non-tumor diseases and its downregulation correlates with adverse outcomes. Further experiments are needed to validate this and other EV markers in the blood of pediatric patients for translational medicine applications.

High-throughput chemogenetic drug screening reveals PKC-RhoA/PKN as a targetable signaling vulnerability in GNAQ-driven uveal melanoma

Uveal melanoma (UM) is the most prevalent cancer of the eye in adults, driven by activating mutation of GNAQ/GNA11; however, there are limited therapies against UM and metastatic UM (mUM). Here, we perform a high-throughput chemogenetic drug screen in GNAQ-mutant UM contrasted with BRAF-mutant cutaneous melanoma, defining the druggable landscape of these distinct melanoma subtypes. Across all compounds, darovasertib demonstrates the highest preferential activity against UM. Our investigation reveals that darovasertib potently inhibits PKC as well as PKN/PRK, an AGC kinase family that is part of the "dark kinome." We find that downstream of the Gαq-RhoA signaling axis, PKN converges with ROCK to control FAK, a mediator of non-canonical Gαq-driven signaling. Strikingly, darovasertib synergizes with FAK inhibitors to halt UM growth and promote cytotoxic cell death in vitro and in preclinical metastatic mouse models, thus exposing a signaling vulnerability that can be exploited as a multimodal precision therapy against mUM.

A novel antifolate suppresses growth of FPGS-deficient cells and overcomes methotrexate resistance

Cancer cells make extensive use of the folate cycle to sustain increased anabolic metabolism. Multiple chemotherapeutic drugs interfere with the folate cycle, including methotrexate and 5-fluorouracil that are commonly applied for the treatment of leukemia and colorectal cancer (CRC), respectively. Despite high success rates, therapy-induced resistance causes relapse at later disease stages. Depletion of folylpolyglutamate synthetase (FPGS), which normally promotes intracellular accumulation and activity of natural folates and methotrexate, is linked to methotrexate and 5-fluorouracil resistance and its association with relapse illustrates the need for improved intervention strategies. Here, we describe a novel antifolate (C1) that, like methotrexate, potently inhibits dihydrofolate reductase and downstream one-carbon metabolism. Contrary to methotrexate, C1 displays optimal efficacy in FPGS-deficient contexts, due to decreased competition with intracellular folates for interaction with dihydrofolate reductase. We show that FPGS-deficient patient-derived CRC organoids display enhanced sensitivity to C1, whereas FPGS-high CRC organoids are more sensitive to methotrexate. Our results argue that polyglutamylation-independent antifolates can be applied to exert selective pressure on FPGS-deficient cells during chemotherapy, using a vulnerability created by polyglutamylation deficiency.

Proteomic and Phosphoproteomic Reprogramming in Epithelial Ovarian Cancer Metastasis

Epithelial ovarian cancer (EOC) is a high-risk cancer presenting with heterogeneous tumors. The high incidence of EOC metastasis from primary tumors to nearby tissues and organs is a major driver of EOC lethality. We used cellular models of spheroid formation and readherence to investigate cellular signaling dynamics in each step toward EOC metastasis. In our system, adherent cells model primary tumors, spheroid formation represents the initiation of metastatic spread, and readherent spheroid cells represent secondary tumors. Proteomic and phosphoproteomic analyses show that spheroid cells are hypoxic and show markers for cell cycle arrest. Aurora kinase B abundance and downstream substrate phosphorylation are significantly reduced in spheroids and readherent cells, explaining their cell cycle arrest phenotype. The proteome of readherent cells is most similar to spheroids, yet greater changes in the phosphoproteome show that spheroid cells stimulate Rho-associated kinase 1 (ROCK1)-mediated signaling, which controls cytoskeletal organization. In spheroids, we found significant phosphorylation of ROCK1 substrates that were reduced in both adherent and readherent cells. Application of the ROCK1-specific inhibitor Y-27632 to spheroids increased the rate of readherence and altered spheroid density. The data suggest ROCK1 inhibition increases EOC metastatic potential. We identified novel pathways controlled by Aurora kinase B and ROCK1 as major drivers of metastatic behavior in EOC cells. Our data show that phosphoproteomic reprogramming precedes proteomic changes that characterize spheroid readherence in EOC metastasis.

Discovery and Synthesis of a Naturally Derived Protein Kinase Inhibitor that Selectively Inhibits Distinct Classes of Serine/Threonine Kinases

The DNAJB1-PRKACA oncogenic gene fusion results in an active kinase enzyme, J-PKAcα, that has been identified as an attractive antitumor target for fibrolamellar hepatocellular carcinoma (FLHCC). A high-throughput assay was used to identify inhibitors of J-PKAcα catalytic activity by screening the NCI Program for Natural Product Discovery (NPNPD) prefractionated natural product library. Purification of the active agent from a single fraction of an Aplidium sp. marine tunicate led to the discovery of two unprecedented alkaloids, aplithianines A (1) and B (2). Aplithianine A (1) showed potent inhibition against J-PKAcα with an IC50 of ∼1 μM in the primary screening assay. In kinome screening, 1 inhibited wild-type PKA with an IC50 of 84 nM. Further mechanistic studies including cocrystallization and X-ray diffraction experiments revealed that 1 inhibited PKAcα catalytic activity by competitively binding to the ATP pocket. Human kinome profiling of 1 against a panel of 370 kinases revealed potent inhibition of select serine/threonine kinases in the CLK and PKG families with IC50 values in the range ∼11-90 nM. An efficient, four-step total synthesis of 1 has been accomplished, enabling further evaluation of aplithianines as biologically relevant kinase inhibitors.

The complex impact of cancer-related missense mutations on the stability and on the biophysical and biochemical properties of MAPK1 and MAPK3 somatic variants

Mitogen-activated protein kinases 1 and 3 (MAPK1 and MAPK3), also called extracellular regulated kinases (ERK2 and ERK1), are serine/threonine kinase activated downstream by the Ras/Raf/MEK/ERK signal transduction cascade that regulates a variety of cellular processes. A dysregulation of MAPK cascade is frequently associated to missense mutations on its protein components and may be related to many pathologies, including cancer. In this study we selected from COSMIC database a set of MAPK1 and MAPK3 somatic variants found in cancer tissues carrying missense mutations distributed all over the MAPK1 and MAPK3 sequences. The proteins were expressed as pure recombinant proteins, and their biochemical and biophysical properties have been studied in comparison with the wild type. The missense mutations lead to changes in the tertiary arrangements of all the variants. The thermodynamic stability of the wild type and variants has been investigated in the non-phosphorylated and in the phosphorylated form. Significant differences in the thermal stabilities of most of the variants have been observed, as well as changes in the catalytic efficiencies. The energetics of the catalytic reaction is affected for all the variants for both the MAPK proteins. The stability changes and the variation in the enzyme catalysis observed for most of MAPK1/3 variants suggest that a local change in a residue, distant from the catalytic site, may have long-distance effects that reflect globally on enzyme stability and functions.

Circulating neutrophils from patients with early breast cancer have distinct subtype-dependent phenotypes

PURPOSE

An elevated number of circulating neutrophils is a poor prognostic factor for breast cancer, where evidence of bone marrow cancer-dependent priming is found. However, how early this priming is detectable remains unclear.

PATIENTS AND METHODS

Here, we investigate changes in circulating neutrophils from newly diagnosed breast cancer patients before any therapeutic interventions. To do this, we assessed their lifespan and their broader intracellular kinase network activation states by using the Pamgene Kinome assay which measures the activity of neutrophil kinases.

RESULTS

We found sub-type specific L-selectin (CD62L) changes in circulating neutrophils as well as perturbations in their overall global kinase activity. Strikingly, breast cancer patients of different subtypes (HR+, HER2+, triple negative) exhibited distinct neutrophil kinase activity patterns indicating that quantifiable perturbations can be detected in circulating neutrophils from early breast cancer patients, that are sensitive to both hormonal and HER-2 status. We also detected an increase in neutrophils lifespan in cancer patients, independently of tumour subtype.

CONCLUSIONS

Our results suggest that the tumour-specific kinase activation patterns in circulating neutrophils may be used in conjunction with other markers to identify patients with cancer from those harbouring only benign lesions of the breast. Given the important role neutrophil in breast cancer progression, the significance of this sub-type of specific priming warrants further investigation.

Coevolution of the ileum with Brk/Ptk6 family kinases confers robustness to ileal homeostasis

Brk/Ptk6, Srms, and Frk constitute a Src-related but distinct family of tyrosine kinases called Brk family kinases (BFKs) in higher vertebrates. To date, however, their biological roles have remained largely unknown. In this study, we generated BFK triple-knockout (BFK/TKO) mice lacking all BFK members using CRISPR/Cas9-mediated genome editing. BFK/TKO mice exhibited impaired intestinal homeostasis, represented by a reduced stem/progenitor cell population and defective recovery from radiation-induced severe mucosal damage, specifically in the ileum, which is the most distal segment of the small intestine. RNA-seq analysis revealed that BFK/TKO ileal epithelium showed markedly elevated IL-22/STAT3 signaling, resulting in the aberrant activation of mucosal immune response and altered composition of the ileal microbiota. Since single- or double-knockout of BFK genes did not elicit such abnormalities, BFKs may redundantly confer robust homeostasis to the ileum, the most recently added intestinal segment that plays crucial roles in nutrient absorption and mucosal immunity. Given that BFK diversification preceded the appearance of the ileum in vertebrate phylogeny, the present study highlights the coevolution of genes and organs, the former of which shapes up the latter in higher vertebrates.

Investigation of the Host Kinome Response to Coronavirus Infection Reveals PI3K/mTOR Inhibitors as Betacoronavirus Antivirals

Host kinases play essential roles in the host cell cycle, innate immune signaling, the stress response to viral infection, and inflammation. Previous work has demonstrated that coronaviruses specifically target kinase cascades to subvert host cell responses to infection and rely upon host kinase activity to phosphorylate viral proteins to enhance replication. Given the number of kinase inhibitors that are already FDA approved to treat cancers, fibrosis, and other human disease, they represent an attractive class of compounds to repurpose for host-targeted therapies against emerging coronavirus infections. To further understand the host kinome response to betacoronavirus infection, we employed multiplex inhibitory bead mass spectrometry (MIB-MS) following MERS-CoV and SARS-CoV-2 infection of human lung epithelial cell lines. Our MIB-MS analyses revealed activation of mTOR and MAPK signaling following MERS-CoV and SARS-CoV-2 infection, respectively. SARS-CoV-2 host kinome responses were further characterized using paired phosphoproteomics, which identified activation of MAPK, PI3K, and mTOR signaling. Through chemogenomic screening, we found that clinically relevant PI3K/mTOR inhibitors were able to inhibit coronavirus replication at nanomolar concentrations similar to direct-acting antivirals. This study lays the groundwork for identifying broad-acting, host-targeted therapies to reduce betacoronavirus replication that can be rapidly repurposed during future outbreaks and epidemics. The proteomics, phosphoproteomics, and MIB-MS datasets generated in this study are available in the Proteomics Identification Database (PRIDE) repository under project identifiers PXD040897 and PXD040901.

Uncovering the dynamics of cellular responses induced by iron-carbohydrate complexes in human macrophages using quantitative proteomics and phosphoproteomics

Iron-carbohydrate complexes are widely used to treat iron deficiencies. Macrophages play a crucial role in the uptake and fate of these nanomedicines, however, how complexed iron carbohydrates are taken up and metabolized by macrophages is still not fully understood. Using a (phospho-)proteomics approach, we assessed differences in protein expression and phosphorylation in M2 macrophages triggered by iron sucrose (IS). Our results show that IS alters the expression of multiple receptors, indicative of a complex entry mechanism. Besides, IS induced an increase in intracellular ferritin, the loss of M2 polarization, protective mechanisms against ferroptosis, and an autophagic response. These data indicate that macrophages can use IS as a source of iron for its storage and later release, however, the excess of iron can cause oxidative stress, which can be successfully regulated by the cells. When comparing IS with ferric carboxymaltose (FCM) and iron isomaltoside-1000 (IIM), complexes with a higher carbohydrate ligand stability, we observed that FCM and IIM are metabolized at a slower rate, and trigger M2 polarization loss to a lower extent. These results indicate that the surface characteristics of the iron-carbohydrate complexes may influence the cell responses. Our data show that the application of (phospho-)proteomics can lead to a better understanding of metabolic processes, including the uptake, biodegradation and bioavailability of nanomedicines.

A multi-omics integrative approach unravels novel genes and pathways associated with senescence escape after targeted therapy in NRAS mutant melanoma

Therapy Induced Senescence (TIS) leads to sustained growth arrest of cancer cells. The associated cytostasis has been shown to be reversible and cells escaping senescence further enhance the aggressiveness of cancers. Chemicals specifically targeting senescent cells, so-called senolytics, constitute a promising avenue for improved cancer treatment in combination with targeted therapies. Understanding how cancer cells evade senescence is needed to optimise the clinical benefits of this therapeutic approach. Here we characterised the response of three different NRAS mutant melanoma cell lines to a combination of CDK4/6 and MEK inhibitors over 33 days. Transcriptomic data show that all cell lines trigger a senescence programme coupled with strong induction of interferons. Kinome profiling revealed the activation of Receptor Tyrosine Kinases (RTKs) and enriched downstream signaling of neurotrophin, ErbB and insulin pathways. Characterisation of the miRNA interactome associates miR-211-5p with resistant phenotypes. Finally, iCell-based integration of bulk and single-cell RNA-seq data identifies biological processes perturbed during senescence and predicts 90 new genes involved in its escape. Overall, our data associate insulin signaling with persistence of a senescent phenotype and suggest a new role for interferon gamma in senescence escape through the induction of EMT and the activation of ERK5 signaling.

Differentiating Benign from Malignant Thyroid Tumors by Kinase Activity Profiling and Dabrafenib BRAF V600E Targeting

Differentiated non-medullary thyroid cancer (NMTC) can be effectively treated by surgery followed by radioactive iodide therapy. However, a small subset of patients shows recurrence due to a loss of iodide transport, a phenotype frequently associated with BRAF V600E mutations. In theory, this should enable the use of existing targeted therapies specifically designed for BRAF V600E mutations. However, in practice, generic or specific drugs aimed at molecular targets identified by next generation sequencing (NGS) are not always beneficial. Detailed kinase profiling may provide additional information to help improve therapy success rates. In this study, we therefore investigated whether serine/threonine kinase (STK) activity profiling can accurately classify benign thyroid lesions and NMTC. We also determined whether dabrafenib (BRAF V600E-specific inhibitor), as well as sorafenib and regorafenib (RAF inhibitors), can differentiate BRAF V600E from non-BRAF V600E thyroid tumors. Using 21 benign and 34 malignant frozen thyroid tumor samples, we analyzed serine/threonine kinase activity using PamChip®peptide microarrays. An STK kinase activity classifier successfully differentiated malignant (26/34; 76%) from benign tumors (16/21; 76%). Of the kinases analyzed, PKC (theta) and PKD1 in particular, showed differential activity in benign and malignant tumors, while oncocytic neoplasia or Graves' disease contributed to erroneous classifications. Ex vivo BRAF V600E-specific dabrafenib kinase inhibition identified 6/92 analyzed peptides, capable of differentiating BRAF V600E-mutant from non-BRAF V600E papillary thyroid cancers (PTCs), an effect not seen with the generic inhibitors sorafenib and regorafenib. In conclusion, STK activity profiling differentiates benign from malignant thyroid tumors and generates unbiased hypotheses regarding differentially active kinases. This approach can serve as a model to select novel kinase inhibitors based on tissue analysis of recurrent thyroid and other cancers.

Combined CDK4/6 and ERK1/2 Inhibition Enhances Antitumor Activity in NF1-Associated Plexiform Neurofibroma

PURPOSE

Plexiform neurofibromas (PNF) are peripheral nerve sheath tumors that cause significant morbidity in persons with neurofibromatosis type 1 (NF1), yet treatment options remain limited. To identify novel therapeutic targets for PNF, we applied an integrated multi-omic approach to quantitatively profile kinome enrichment in a mouse model that has predicted therapeutic responses in clinical trials for NF1-associated PNF with high fidelity.

EXPERIMENTAL DESIGN

Utilizing RNA sequencing combined with chemical proteomic profiling of the functionally enriched kinome using multiplexed inhibitor beads coupled with mass spectrometry, we identified molecular signatures predictive of response to CDK4/6 and RAS/MAPK pathway inhibition in PNF. Informed by these results, we evaluated the efficacy of the CDK4/6 inhibitor, abemaciclib, and the ERK1/2 inhibitor, LY3214996, alone and in combination in reducing PNF tumor burden in Nf1flox/flox;PostnCre mice.

RESULTS

Converging signatures of CDK4/6 and RAS/MAPK pathway activation were identified within the transcriptome and kinome that were conserved in both murine and human PNF. We observed robust additivity of the CDK4/6 inhibitor, abemaciclib, in combination with the ERK1/2 inhibitor, LY3214996, in murine and human NF1(Nf1) mutant Schwann cells. Consistent with these findings, the combination of abemaciclib (CDK4/6i) and LY3214996 (ERK1/2i) synergized to suppress molecular signatures of MAPK activation and exhibited enhanced antitumor activity in Nf1flox/flox;PostnCre mice in vivo.

CONCLUSIONS

These findings provide rationale for the clinical translation of CDK4/6 inhibitors alone and in combination with therapies targeting the RAS/MAPK pathway for the treatment of PNF and other peripheral nerve sheath tumors in persons with NF1.

Dynamic applicability domain (dAD): compound-target binding affinity estimates with local conformal prediction

MOTIVATION

Increasing efforts are being made in the field of machine learning to advance the learning of robust and accurate models from experimentally measured data and enable more efficient drug discovery processes. The prediction of binding affinity is one of the most frequent tasks of compound bioactivity modelling. Learned models for binding affinity prediction are assessed by their average performance on unseen samples, but point predictions are typically not provided with a rigorous confidence assessment. Approaches, such as the conformal predictor framework equip conventional models with a more rigorous assessment of confidence for individual point predictions. In this article, we extend the inductive conformal prediction framework for interaction data, in particular the compound-target binding affinity prediction task. The new framework is based on dynamically defined calibration sets that are specific for each testing pair and provides prediction assessment in the context of calibration pairs from its compound-target neighbourhood, enabling improved estimates based on the local properties of the prediction model.

RESULTS

The effectiveness of the approach is benchmarked on several publicly available datasets and tested in realistic use-case scenarios with increasing levels of difficulty on a complex compound-target binding affinity space. We demonstrate that in such scenarios, novel approach combining applicability domain paradigm with conformal prediction framework, produces superior confidence assessment with valid and more informative prediction regions compared to other 'state-of-the-art' conformal prediction approaches.

AVAILABILITY AND IMPLEMENTATION

Dataset and the code are available on GitHub (https://github.com/mlkr-rbi/dAD).

Systematic analysis of the impact of phosphorylation and O-GlcNAcylation on protein subcellular localization

The subcellular localization of proteins is critical for their functions in eukaryotic cells and is tightly correlated with protein modifications. Here, we comprehensively investigate the nuclear-cytoplasmic distributions of the phosphorylated, O-GlcNAcylated, and non-modified forms of proteins to dissect the correlation between protein distribution and modifications. Phosphorylated and O-GlcNAcylated proteins have overall higher nuclear distributions than non-modified ones. Different distributions among the phosphorylated, O-GlcNAcylated, and non-modified forms of proteins are associated with protein size, structure, and function, as well as local environment and adjacent residues around modification sites. Moreover, we perform site-mutagenesis experiments using phosphomimetic and phospho-null mutants of two proteins to validate the proteomic results. Additionally, the effects of the OGT/OGA inhibition on glycoprotein distribution are systematically investigated, and the distribution changes of glycoproteins are related to their abundance changes under the inhibitions. Systematic investigation of the relationship between protein modification and localization advances our understanding of protein functions.

SHIP1 Is Present but Strongly Downregulated in T-ALL, and after Restoration Suppresses Leukemia Growth in a T-ALL Xenotransplantation Mouse Model

Acute lymphoblastic leukemia (ALL) is the most common cause of cancer-related death in children. Despite significantly increased chances of cure, especially for high-risk ALL patients, it still represents a poor prognosis for a substantial fraction of patients. Misregulated proteins in central switching points of the cellular signaling pathways represent potentially important therapeutic targets. Recently, the inositol phosphatase SHIP1 (SH2-containing inositol 5-phosphatase) has been considered as a tumor suppressor in leukemia. SHIP1 serves as an important negative regulator of the PI3K/AKT signaling pathway, which is frequently constitutively activated in primary T-ALL. In contrast to other reports, we show for the first time that SHIP1 has not been lost in T-ALL cells, but is strongly downregulated. Reduced expression of SHIP1 leads to an increased activation of the PI3K/AKT signaling pathway. SHIP1-mRNA expression is frequently reduced in primary T-ALL samples, which is recapitulated by the decrease in SHIP1 expression at the protein level in seven out of eight available T-ALL patient samples. In addition, we investigated the change in the activity profile of tyrosine and serine/threonine kinases after the restoration of SHIP1 expression in Jurkat T-ALL cells. The tyrosine kinase receptor subfamilies of NTRK and PDGFR, which are upregulated in T-ALL subgroups with low SHIP1 expression, are significantly disabled after SHIP1 reconstitution. Lentiviral-mediated reconstitution of SHIP1 expression in Jurkat cells points to a decreased cellular proliferation upon transplantation into NSG mice in comparison to the control cohort. Together, our findings will help to elucidate the complex network of cell signaling proteins, further support a functional role for SHIP1 as tumor suppressor in T-ALL and, much more importantly, show that full-length SHIP1 is expressed in T-ALL samples.

Loss of hepatic PPARα in mice causes hypertension and cardiovascular disease

The leading cause of death in patients with nonalcoholic fatty liver disease (NAFLD) is cardiovascular disease (CVD). However, the mechanisms are unknown. Mice deficient in hepatocyte proliferator-activated receptor-α (PPARα) (PparaHepKO) exhibit hepatic steatosis on a regular chow diet, making them prone to manifesting NAFLD. We hypothesized that the PparaHepKO mice might be predisposed to poorer cardiovascular phenotypes due to increased liver fat content. Therefore, we used PparaHepKO and littermate control mice fed a regular chow diet to avoid complications with a high-fat diet, such as insulin resistance and increased adiposity. After 30 wk on a standard diet, male PparaHepKO mice exhibited elevated hepatic fat content compared with littermates as measured by Echo MRI (11.95 ± 1.4 vs. 3.74 ± 1.4%, P < 0.05), hepatic triglycerides (1.4 ± 0.10 vs. 0.3 ± 0.01 mM, P < 0.05), and Oil Red O staining, despite body weight, fasting blood glucose, and insulin levels being the same as controls. The PparaHepKO mice also displayed elevated mean arterial blood pressure (121 ± 4 vs. 108 ± 2 mmHg, P < 0.05), impaired diastolic function, cardiac remodeling, and enhanced vascular stiffness. To determine mechanisms controlling the increase in stiffness in the aorta, we used state-of-the-art PamGene technology to measure kinase activity in this tissue. Our data suggest that the loss of hepatic PPARα induces alterations in the aortas that reduce the kinase activity of tropomyosin receptor kinases and p70S6K kinase, which might contribute to the pathogenesis of NAFLD-induced CVD. These data indicate that hepatic PPARα protects the cardiovascular system through some as-of-yet undefined mechanism.

Multiplex single-cell chemical genomics reveals the kinase dependence of the response to targeted therapy

Chemical genetic screens are a powerful tool for exploring how cancer cells' response to drugs is shaped by their mutations, yet they lack a molecular view of the contribution of individual genes to the response to exposure. Here, we present sci-Plex-Gene-by-Environment (sci-Plex-GxE), a platform for combined single-cell genetic and chemical screening at scale. We highlight the advantages of large-scale, unbiased screening by defining the contribution of each of 522 human kinases to the response of glioblastoma to different drugs designed to abrogate signaling from the receptor tyrosine kinase pathway. In total, we probed 14,121 gene-by-environment combinations across 1,052,205 single-cell transcriptomes. We identify an expression signature characteristic of compensatory adaptive signaling regulated in a MEK/MAPK-dependent manner. Further analyses aimed at preventing adaptation revealed promising combination therapies, including dual MEK and CDC7/CDK9 or NF-kB inhibitors, as potent means of preventing transcriptional adaptation of glioblastoma to targeted therapy.

Large-Scale Modeling of Sparse Protein Kinase Activity Data

Protein kinases are a protein family that plays an important role in several complex diseases such as cancer and cardiovascular and immunological diseases. Protein kinases have conserved ATP binding sites, which when targeted can lead to similar activities of inhibitors against different kinases. This can be exploited to create multitarget drugs. On the other hand, selectivity (lack of similar activities) is desirable in order to avoid toxicity issues. There is a vast amount of protein kinase activity data in the public domain, which can be used in many different ways. Multitask machine learning models are expected to excel for these kinds of data sets because they can learn from implicit correlations between tasks (in this case activities against a variety of kinases). However, multitask modeling of sparse data poses two major challenges: (i) creating a balanced train-test split without data leakage and (ii) handling missing data. In this work, we construct a protein kinase benchmark set composed of two balanced splits without data leakage, using random and dissimilarity-driven cluster-based mechanisms, respectively. This data set can be used for benchmarking and developing protein kinase activity prediction models. Overall, the performance on the dissimilarity-driven cluster-based split is lower than on random split-based sets for all models, indicating poor generalizability of models. Nevertheless, we show that multitask deep learning models, on this very sparse data set, outperform single-task deep learning and tree-based models. Finally, we demonstrate that data imputation does not improve the performance of (multitask) models on this benchmark set.

Automated Enrichment of Phosphotyrosine Peptides for High-Throughput Proteomics

Phosphotyrosine (pY) enrichment is critical for expanding the fundamental and clinical understanding of cellular signaling by mass spectrometry-based proteomics. However, current pY enrichment methods exhibit a high cost per sample and limited reproducibility due to expensive affinity reagents and manual processing. We present rapid-robotic phosphotyrosine proteomics (R2-pY), which uses a magnetic particle processor and pY superbinders or antibodies. R2-pY can handle up to 96 samples in parallel, requires 2 days to go from cell lysate to mass spectrometry injections, and results in global proteomic, phosphoproteomic, and tyrosine-specific phosphoproteomic samples. We benchmark the method on HeLa cells stimulated with pervanadate and serum and report over 4000 unique pY sites from 1 mg of peptide input, strong reproducibility between replicates, and phosphopeptide enrichment efficiencies above 99%. R2-pY extends our previously reported R2-P2 proteomic and global phosphoproteomic sample preparation framework, opening the door to large-scale studies of pY signaling in concert with global proteome and phosphoproteome profiling.

Comparative biochemical kinase activity analysis identifies rivoceranib as a highly selective VEGFR2 inhibitor

Vascular endothelial growth factor receptor 2 (VEGFR2), a key regulator of tumor angiogenesis, is highly expressed across numerous tumor types and has been an attractive target for anti-cancer therapy. However, clinical application of available VEGFR2 inhibitors has been challenged by limited efficacy and a wide range of side effects, potentially due to inadequate selectivity for VEGFR2. Thus, development of potent VEGFR2 inhibitors with improved selectivity is needed. Rivoceranib is an orally administered tyrosine kinase inhibitor that potently and selectively targets VEGFR2. A comparative understanding of the potency and selectivity of rivoceranib and approved inhibitors of VEGFR2 is valuable to inform rationale for therapy selection in the clinic. Here, we performed biochemical analyses of the kinase activity of VEGFR2 and of a panel of 270 kinases to compare rivoceranib to 10 FDA-approved kinase inhibitors ("reference inhibitors") with known activity against VEGFR2. Rivoceranib demonstrated potency within the range of the reference inhibitors, with a VEGFR2 kinase inhibition IC₅₀ value of 16 nM. However, analysis of residual kinase activity of the panel of 270 kinases showed that rivoceranib displayed greater selectivity for VEGFR2 compared with the reference inhibitors. Differences in selectivity among compounds within the observed range of potency of VEGFR2 kinase inhibition are clinically relevant, as toxicities associated with available VEGFR2 inhibitors are thought to be partly due to their effects against kinases other than VEGFR2. Together, this comparative biochemical analysis highlights the potential for rivoceranib to address clinical limitations associated with off-target effects of currently available VEGFR2 inhibitors.

Establishment of extensive artificial intelligence models for kinase inhibitor prediction: Identification of novel PDGFRB inhibitors

Identifying hit compounds is an important step in drug development. Unfortunately, this process continues to be a challenging task. Several machine learning models have been generated to aid in simplifying and improving the prediction of candidate compounds. Models tuned for predicting kinase inhibitors have been established. However, an effective model can be limited by the size of the chosen training dataset. In this study, we tested several machine learning models to predict potential kinase inhibitors. A dataset was curated from a number of publicly available repositories. This resulted in a comprehensive dataset covering more than half of the human kinome. More than 2,000 kinase models were established using different model approaches. The performances of the models were compared, and the Keras-MLP model was determined to be the best performing model. The model was then used to screen a chemical library for potential inhibitors targeting platelet-derived growth factor receptor-β (PDGFRB). Several PDGFRB candidates were selected, and in vitro assays confirmed four compounds with PDGFRB inhibitory activity and IC₅₀ values in the nanomolar range. These results show the effectiveness of machine learning models trained on the reported dataset. This report would aid in the establishment of machine learning models as well as in the discovery of novel kinase inhibitors.

The Stiffness-Sensitive Transcriptome of Human Tendon Stromal Cells

Extracellular matrix stiffness is a major regulator of cellular states. Stiffness-sensing investigations are typically performed using cells that have acquired "mechanical memory" through prolonged conditioning in rigid environments, e.g., tissue culture plastic (TCP). This potentially masks the full extent of the matrix stiffness-driven mechanosensing programs. Here, a biomaterial composed of 2D mechanovariant silicone substrates with simplified and scalable surface biofunctionalization chemistry is developed to facilitate large-scale cell culture expansion processes. Using RNA sequencing, stiffness-mediated mechano-responses of human tendon-derived stromal cells are broadly mapped. Matrix elasticity (E.) approximating tendon microscale stiffness range (E. ≈ 35 kPa) distinctly favors transcriptional programs related to chromatin remodeling and Hippo signaling; whereas compliant stiffnesses (E. ≈ 2 kPa) are enriched in processes related to cell stemness, synapse assembly, and angiogenesis. While tendon stromal cells undergo dramatic phenotypic drift on conventional TCP, mechanovariant substrates abrogate this activation with tenogenic stiffnesses inducing a transcriptional program that strongly correlates with established tendon tissue-specific expression signature. Computational inference predicts that AKT1 and ERK1/2 are major stiffness-sensing signaling hubs. Together, these findings highlight how matrix biophysical cues may dictate the transcriptional identity of tendon cells, and how matrix mechano-reciprocity regulates diverse sets of previously underappreciated mechanosensitive processes in tendon fibroblasts.

Phosphoproteomic analysis of metformin signaling in colorectal cancer cells elucidates mechanism of action and potential therapeutic opportunities

BACKGROUND

The biguanide drug metformin is a safe and widely prescribed drug for type 2 diabetes. Interestingly, hundreds of clinical trials have been set to evaluate the potential role of metformin in the prevention and treatment of cancer including colorectal cancer (CRC). However, the "metformin signaling" remains controversial.

AIMS AND METHODS

To interrogate cell signaling induced by metformin in CRC and explore the druggability of the metformin-rewired phosphorylation network, we performed integrative analysis of phosphoproteomics, bioinformatics, and cell proliferation assays on a panel of 12 molecularly heterogeneous CRC cell lines. Using the high-resolute data-independent analysis mass spectrometry (DIA-MS), we monitored a total of 10,142 proteins and 56,080 phosphosites (P-sites) in CRC cells upon a short- and a long-term metformin treatment.

RESULTS AND CONCLUSIONS

We found that metformin tended to primarily remodel cell signaling in the long-term and only minimally regulated the total proteome expression levels. Strikingly, the phosphorylation signaling response to metformin was highly heterogeneous in the CRC panel, based on a network analysis inferring kinase/phosphatase activities and cell signaling reconstruction. A "MetScore" was determined to assign the metformin relevance of each P-site, revealing new and robust phosphorylation nodes and pathways in metformin signaling. Finally, we leveraged the metformin P-site signature to identify pharmacodynamic interactions and confirmed a number of candidate metformin-interacting drugs, including navitoclax, a BCL-2/BCL-xL inhibitor. Together, we provide a comprehensive phosphoproteomic resource to explore the metformin-induced cell signaling for potential cancer therapeutics. This resource can be accessed at https://yslproteomics.shinyapps.io/Metformin/.

Finnish-specific AKT2 gene variant leads to impaired insulin signalling in myotubes

Finnish-specific gene variant p.P50T/AKT2 (minor allele frequency (MAF) = 1.1%) is associated with insulin resistance and increased predisposition to type 2 diabetes. Here, we have investigated in vitro the impact of the gene variant on glucose metabolism and intracellular signalling in human primary skeletal muscle cells, which were established from 14 male p.P50T/AKT2 variant carriers and 14 controls. Insulin-stimulated glucose uptake and glucose incorporation into glycogen were detected with 2-[1,2-3H]-deoxy-D-glucose and D-[14C]-glucose, respectively, and the rate of glycolysis was measured with a Seahorse XFe96 analyzer. Insulin signalling was investigated with Western blotting. The binding of variant and control AKT2-PH domains to phosphatidylinositol (3,4,5)-trisphosphate (PI(3,4,5)P3) was assayed using PIP StripsTM Membranes. Protein tyrosine kinase and serine-threonine kinase assays were performed using the PamGene® kinome profiling system. Insulin-stimulated glucose uptake and glycogen synthesis in myotubes in vitro were not significantly affected by the genotype. However, the insulin-stimulated glycolytic rate was impaired in variant myotubes. Western blot analysis showed that insulin-stimulated phosphorylation of AKT-Thr308, AS160-Thr642 and GSK3β-Ser9 was reduced in variant myotubes compared to controls. The binding of variant AKT2-PH domain to PI(3,4,5)P3 was reduced as compared to the control protein. PamGene® kinome profiling revealed multiple differentially phosphorylated kinase substrates, e.g. calmodulin, between the genotypes. Further in silico upstream kinase analysis predicted a large-scale impairment in activities of kinases participating, for example, in intracellular signal transduction, protein translation and cell cycle events. In conclusion, myotubes from p.P50T/AKT2 variant carriers show multiple signalling alterations which may contribute to predisposition to insulin resistance and T2D in the carriers of this signalling variant.

Prenatal methadone exposure selectively alters protein expression in primary motor cortex: Implications for synaptic function

As problematic opioid use has reached epidemic levels over the past 2 decades, the annual prevalence of opioid use disorder (OUD) in pregnant women has also increased 333%. Yet, how opioids affect the developing brain of offspring from mothers experiencing OUD remains understudied and not fully understood. Animal models of prenatal opioid exposure have discovered many deficits in the offspring of prenatal opioid exposed mothers, such as delays in the development of sensorimotor function and long-term locomotive hyperactivity. In attempt to further understand these deficits and link them with protein changes driven by prenatal opioid exposure, we used a mouse model of prenatal methadone exposure (PME) and preformed an unbiased multi-omic analysis across many sensoriomotor brain regions known to interact with opioid exposure. The effects of PME exposure on the primary motor cortex (M1), primary somatosensory cortex (S1), the dorsomedial striatum (DMS), and dorsolateral striatum (DLS) were assessed using quantitative proteomics and phosphoproteomics. PME drove many changes in protein and phosphopeptide abundance across all brain regions sampled. Gene and gene ontology enrichments were used to assess how protein and phosphopeptide changes in each brain region were altered. Our findings showed that M1 was uniquely affected by PME in comparison to other brain regions. PME uniquely drove changes in M1 glutamatergic synapses and synaptic function. Immunohistochemical analysis also identified anatomical differences in M1 for upregulating the density of glutamatergic and downregulating the density of GABAergic synapses due to PME. Lastly, comparisons between M1 and non-M1 multi-omics revealed conserved brain wide changes in phosphopeptides associated with synaptic activity and assembly, but only specific protein changes in synapse activity and assembly were represented in M1. Together, our studies show that lasting changes in synaptic function driven by PME are largely represented by protein and anatomical changes in M1, which may serve as a starting point for future experimental and translational interventions that aim to reverse the adverse effects of PME on offspring.

Suppressing Hepatic UGT1A1 Increases Plasma Bilirubin, Lowers Plasma Urobilin, Reorganizes Kinase Signaling Pathways and Lipid Species and Improves Fatty Liver Disease

Several population studies have observed lower serum bilirubin levels in patients with non-alcoholic fatty liver disease (NAFLD). Yet, treatments to target this metabolic phenotype have not been explored. Therefore, we designed an N-Acetylgalactosamine (GalNAc) labeled RNAi to target the enzyme that clears bilirubin from the blood, the UGT1A1 glucuronyl enzyme (GNUR). In this study, male C57BL/6J mice were fed a high-fat diet (HFD, 60%) for 30 weeks to induce NAFLD and were treated subcutaneously with GNUR or sham (CTRL) once weekly for six weeks while continuing the HFD. The results show that GNUR treatments significantly raised plasma bilirubin levels and reduced plasma levels of the bilirubin catabolized product, urobilin. We show that GNUR decreased liver fat content and ceramide production via lipidomics and lowered fasting blood glucose and insulin levels. We performed extensive kinase activity analyses using our PamGene PamStation kinome technology and found a reorganization of the kinase pathways and a significant decrease in inflammatory mediators with GNUR versus CTRL treatments. These results demonstrate that GNUR increases plasma bilirubin and reduces plasma urobilin, reducing NAFLD and inflammation and improving overall liver health. These data indicate that UGT1A1 antagonism might serve as a treatment for NAFLD and may improve obesity-associated comorbidities.

Modulation of RNA splicing enhances response to BCL2 inhibition in leukemia

Therapy resistance is a major challenge in the treatment of cancer. Here, we performed CRISPR-Cas9 screens across a broad range of therapies used in acute myeloid leukemia to identify genomic determinants of drug response. Our screens uncover a selective dependency on RNA splicing factors whose loss preferentially enhances response to the BCL2 inhibitor venetoclax. Loss of the splicing factor RBM10 augments response to venetoclax in leukemia yet is completely dispensable for normal hematopoiesis. Combined RBM10 and BCL2 inhibition leads to mis-splicing and inactivation of the inhibitor of apoptosis XIAP and downregulation of BCL2A1, an anti-apoptotic protein implicated in venetoclax resistance. Inhibition of splicing kinase families CLKs (CDC-like kinases) and DYRKs (dual-specificity tyrosine-regulated kinases) leads to aberrant splicing of key splicing and apoptotic factors that synergize with venetoclax, and overcomes resistance to BCL2 inhibition. Our findings underscore the importance of splicing in modulating response to therapies and provide a strategy to improve venetoclax-based treatments.

Automated Enrichment of Phosphotyrosine Peptides for High-Throughput Proteomics

Phosphotyrosine (pY) enrichment is critical for expanding fundamental and clinical understanding of cellular signaling by mass spectrometry-based proteomics. However, current pY enrichment methods exhibit a high cost per sample and limited reproducibility due to expensive affinity reagents and manual processing. We present rapid-robotic phosphotyrosine proteomics (R2-pY), which uses a magnetic particle processor and pY superbinders or antibodies. R2-pY handles 96 samples in parallel, requires 2 days to go from cell lysate to mass spectrometry injections, and results in global proteomic, phosphoproteomic and tyrosine specific phosphoproteomic samples. We benchmark the method on HeLa cells stimulated with pervanadate and serum and report over 4000 unique pY sites from 1 mg of peptide input, strong reproducibility between replicates, and phosphopeptide enrichment efficiencies above 99%. R2-pY extends our previously reported R2-P2 proteomic and global phosphoproteomic sample preparation framework, opening the door to large-scale studies of pY signaling in concert with global proteome and phosphoproteome profiling.

Evolutionary divergence in the conformational landscapes of tyrosine vs serine/threonine kinases

Inactive conformations of protein kinase catalytic domains where the DFG motif has a "DFG-out" orientation and the activation loop is folded present a druggable binding pocket that is targeted by FDA-approved 'type-II inhibitors' in the treatment of cancers. Tyrosine kinases (TKs) typically show strong binding affinity with a wide spectrum of type-II inhibitors while serine/threonine kinases (STKs) usually bind more weakly which we suggest here is due to differences in the folded to extended conformational equilibrium of the activation loop between TKs vs. STKs. To investigate this, we use sequence covariation analysis with a Potts Hamiltonian statistical energy model to guide absolute binding free-energy molecular dynamics simulations of 74 protein-ligand complexes. Using the calculated binding free energies together with experimental values, we estimated free-energy costs for the large-scale (~17-20 Å) conformational change of the activation loop by an indirect approach, circumventing the very challenging problem of simulating the conformational change directly. We also used the Potts statistical potential to thread large sequence ensembles over active and inactive kinase states. The structure-based and sequence-based analyses are consistent; together they suggest TKs evolved to have free-energy penalties for the classical 'folded activation loop' DFG-out conformation relative to the active conformation, that is, on average, 4-6 kcal/mol smaller than the corresponding values for STKs. Potts statistical energy analysis suggests a molecular basis for this observation, wherein the activation loops of TKs are more weakly 'anchored' against the catalytic loop motif in the active conformation and form more stable substrate-mimicking interactions in the inactive conformation. These results provide insights into the molecular basis for the divergent functional properties of TKs and STKs, and have pharmacological implications for the target selectivity of type-II inhibitors.

Targeting Gatekeeper Mutations for Kinase Drug Discovery

Clinically acquired resistance is a major challenge in cancer therapies with small-molecule kinase inhibitors (SMKIs). Gatekeeper mutations in the ATP-binding pocket of kinases are the most common mutations leading to acquired resistance. To date, seven new-generation kinase inhibitors targeting gatekeeper mutations have been approved by the FDA; however, the clinical need is still unmet. Here, we systematically summarize the types of gatekeeper mutations across the kinase family, the structural basis for acquired resistance, and newly developed SMKIs targeting gatekeeper mutations as well as highlight the opportunities and challenges of kinase drug discovery for targeting gatekeeper mutations.