Benchmarking substrate-based kinase activity inference using phosphoproteomic data

A computational tool to infer enzyme activity using post-translational modification profiling data

Enzymes play a pivotal role in orchestrating complex cellular responses to external stimuli and environmental changes through signal transduction pathways. Despite their crucial roles, measuring enzyme activities is typically indirect and performed on a smaller scale, unlike protein abundance measured by high-throughput proteomics. Moreover, it is challenging to derive the activity of enzymes from proteome-wide post-translational modification (PTM) profiling data. To address this challenge, we introduce enzyme activity inference with structural equation modeling under the JUMP umbrella (JUMPsem), a novel computational tool designed to infer enzyme activity using PTM profiling data. We demonstrate that the JUMPsem program enables estimating kinase activities using phosphoproteome data, ubiquitin E3 ligase activities from the ubiquitinome, and histone acetyltransferase (HAT) activities based on the acetylome. In addition, JUMPsem is capable of establishing novel enzyme-substrate relationships through searching motif sequences. JUMPsem outperforms widely used kinase activity tools, such as IKAP and KSEA, in terms of the number of kinases and the computational speed. The JUMPsem program is scalable and publicly available as an open-source R package and user-friendly web-based R/Shiny app. Collectively, JUMPsem provides an improved tool for inferring protein enzyme activities, potentially facilitating targeted drug development.

PhosNetVis: A web-based tool for fast kinase-substrate enrichment analysis and interactive 2D/3D network visualizations of phosphoproteomics data

Protein phosphorylation involves the reversible modification of a protein (substrate) residue by another protein (kinase). Liquid chromatography-mass spectrometry studies are rapidly generating massive protein phosphorylation datasets across multiple conditions. Researchers then must infer kinases responsible for changes in phosphosites of each substrate. However, tools that infer kinase-substrate interactions (KSIs) are not optimized to interactively explore the resulting large and complex networks, significant phosphosites, and states. There is thus an unmet need for a tool that facilitates user-friendly analysis, interactive exploration, visualization, and communication of phosphoproteomics datasets. We present PhosNetVis, a web-based tool for researchers of all computational skill levels to easily infer, generate, and interactively explore KSI networks in 2D or 3D by streamlining phosphoproteomics data analysis steps within a single tool. PhostNetVis lowers barriers for researchers by rapidly generating high-quality visualizations to gain biological insights from their phosphoproteomics datasets. It is available at https://gumuslab.github.io/PhosNetVis/.

PhosNetVis: A web-based tool for fast kinase-substrate enrichment analysis and interactive 2D/3D network visualizations of phosphoproteomics data

Protein phosphorylation involves the reversible modification of a protein (substrate) residue by another protein (kinase). Liquid chromatography-mass spectrometry studies are rapidly generating massive protein phosphorylation datasets across multiple conditions. Researchers then must infer kinases responsible for changes in phosphosites of each substrate. However, tools that infer kinase-substrate interactions (KSIs) are not optimized to interactively explore the resulting large and complex networks, significant phosphosites, and states. There is thus an unmet need for a tool that facilitates user-friendly analysis, interactive exploration, visualization, and communication of phosphoproteomics datasets. We present PhosNetVis, a web-based tool for researchers of all computational skill levels to easily infer, generate and interactively explore KSI networks in 2D or 3D by streamlining phosphoproteomics data analysis steps within a single tool. PhostNetVis lowers barriers for researchers in rapidly generating high-quality visualizations to gain biological insights from their phosphoproteomics datasets. It is available at: https://gumuslab.github.io/PhosNetVis/.

Proteomic Sensors for Quantitative, Multiplexed and Spatial Monitoring of Kinase Signaling

Understanding kinase action requires precise quantitative measurements of their activity in vivo . In addition, the ability to capture spatial information of kinase activity is crucial to deconvolute complex signaling networks, interrogate multifaceted kinase actions, and assess drug effects or genetic perturbations. Here we developed a proteomic kinase activity sensor platform (ProKAS) for the analysis of kinase signaling using mass spectrometry. ProKAS is based on a tandem array of peptide sensors with amino acid barcodes that allow multiplexed analysis for spatial, kinetic, and screening applications. We engineered a ProKAS module to simultaneously monitor the activities of the DNA damage response kinases ATR, ATM, and CHK1 in response to genotoxic drugs, while also uncovering differences between these signaling responses in the nucleus, cytosol, and replication factories. Furthermore, we developed an in silico approach for the rational design of specific substrate peptides expandable to other kinases. Overall, ProKAS is a novel versatile system for systematically and spatially probing kinase action in cells.

ERK synchronizes embryonic cleavages in Drosophila

Extracellular-signal-regulated kinase (ERK) signaling controls development and homeostasis and is genetically deregulated in human diseases, including neurocognitive disorders and cancers. Although the list of ERK functions is vast and steadily growing, the full spectrum of processes controlled by any specific ERK activation event remains unknown. Here, we show how ERK functions can be systematically identified using targeted perturbations and global readouts of ERK activation. Our experimental model is the Drosophila embryo, where ERK signaling at the embryonic poles has thus far only been associated with the transcriptional patterning of the future larva. Through a combination of live imaging and phosphoproteomics, we demonstrated that ERK activation at the poles is also critical for maintaining the speed and synchrony of embryonic cleavages. The presented approach to interrogating phosphorylation networks identifies a hidden function of a well-studied signaling event and sets the stage for similar studies in other organisms.

PhosX: data-driven kinase activity inference from phosphoproteomics experiments

SUMMARY

The inference of kinase activity from phosphoproteomics data can point to causal mechanisms driving signalling processes and potential drug targets. Identifying the kinases whose change in activity explains the observed phosphorylation profiles, however, remains challenging, and constrained by the manually curated knowledge of kinase-substrate associations. Recently, experimentally determined substrate sequence specificities of human kinases have become available, but robust methods to exploit this new data for kinase activity inference are still missing. We present PhosX, a method to estimate differential kinase activity from phosphoproteomics data that combines state-of-the-art statistics in enrichment analysis with kinases' substrate sequence specificity information. Using a large phosphoproteomics dataset with known differentially regulated kinases we show that our method identifies upregulated and downregulated kinases by only relying on the input phosphopeptides' sequences and intensity changes. We find that PhosX outperforms the currently available approach for the same task, and performs better or similarly to state-of-the-art methods that rely on previously known kinase-substrate associations. We therefore recommend its use for data-driven kinase activity inference.

AVAILABILITY AND IMPLEMENTATION

PhosX is implemented in Python, open-source under the Apache-2.0 licence, and distributed on the Python Package Index. The code is available on GitHub (https://github.com/alussana/phosx).

Proteogenomic analysis reveals non-small cell lung cancer subtypes predicting chromosome instability, and tumor microenvironment

Non-small cell lung cancer (NSCLC) is histologically classified into lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LSCC). However, some tumors are histologically ambiguous and other pathophysiological features or microenvironmental factors may be more prominent. Here we report integrative multiomics analyses using data for 229 patients from a Korean NSCLC cohort and 462 patients from previous multiomics studies. Histological examination reveals five molecular subtypes, one of which is a NSCLC subtype with PI3K-Akt pathway upregulation, showing a high proportion of metastasis and poor survival outcomes regardless of any specific NSCLC histology. Proliferative subtypes are present in LUAD and LSCC, which show strong associations with whole genome doubling (WGD) events. Comprehensive characterization of the immune microenvironment reveals various immune cell compositions and neoantigen loads across molecular subtypes, which predicting different prognoses. Immunological subtypes exhibit a hot tumor-enriched state and a higher efficacy of adjuvant therapy.

Role of DNA damage response in cyclophosphamide-induced premature ovarian failure in mice

AIM

To investigate the DNA damage response (DDR) in a cyclophosphamide (CTX)-induced mouse model of premature ovarian failure (POF).

METHODS

The POF model was established by injecting mice with CTX. The body, ovarian weights, the estrus cycle, and pathological changes of the ovaries were recorded. The serum levels of 17 β-estradiol (E2) and follicle-stimulating hormone (FSH) were measured. The expression of Ki67, β-galactosidase (β-gal), p21, p53, γH2AX, and pATM in ovarian tissues was detected by immunohistochemistry. The expression of β-gal, γH2AX, and pATM was analyzed by immunofluorescence staining of primary cultured granulosa cells (GCs).

RESULTS

The body and ovarian weights decreased, the estrus cycles were erratic, and the FSH level increased, whereas the E2 level decreased in POF mice compared to controls. The pathological consequences of POF revealed an increase in atretic follicles, corpus luteum, and primordial follicles and a decrease in the number of primary, secondary, and tertiary follicles. Ki67 expression was reduced, β-gal, p21, p53, γH2AX, and pATM expression were elevated in the ovaries of POF mice. The expression of β-gal, γH2AX, and pATM increased in GCs with the concentration in a time-dependent manner.

CONCLUSION

In total, CTX induced POF in mice, which was mediated by the DDR pathway of ATM-P53-P21.

Calorie restriction and rapamycin distinctly mitigate aging-associated protein phosphorylation changes in mouse muscles

Calorie restriction (CR) and treatment with rapamycin (RM), an inhibitor of the mTORC1 growth-promoting signaling pathway, are known to slow aging and promote health from worms to humans. At the transcriptome and proteome levels, long-term CR and RM treatments have partially overlapping effects, while their impact on protein phosphorylation within cellular signaling pathways have not been compared. Here we measured the phosphoproteomes of soleus, tibialis anterior, triceps brachii and gastrocnemius muscles from adult (10 months) and 30-month-old (aged) mice receiving either a control, a calorie restricted or an RM containing diet from 15 months of age. We reproducibly detected and extensively analyzed a total of 6960 phosphosites, 1415 of which are not represented in standard repositories. We reveal the effect of these interventions on known mTORC1 pathway substrates, with CR displaying greater between-muscle variation than RM. Overall, CR and RM have largely consistent, but quantitatively distinct long-term effects on the phosphoproteome, mitigating age-related changes to different degrees. Our data expands the catalog of protein phosphorylation sites in the mouse, providing important information regarding their tissue-specificity, and revealing the impact of long-term nutrient-sensing pathway inhibition on mouse skeletal muscle.

Waves of regulated protein expression and phosphorylation rewire the proteome to drive gametogenesis in budding yeast

Sexually reproducing eukaryotes employ a developmentally regulated cell division program-meiosis-to generate haploid gametes from diploid germ cells. To understand how gametes arise, we generated a proteomic census encompassing the entire meiotic program of budding yeast. We found that concerted waves of protein expression and phosphorylation modify nearly all cellular pathways to support meiotic entry, meiotic progression, and gamete morphogenesis. Leveraging this comprehensive resource, we pinpointed dynamic changes in mitochondrial components and showed that phosphorylation of the FoF1-ATP synthase complex is required for efficient gametogenesis. Furthermore, using cryoET as an orthogonal approach to visualize mitochondria, we uncovered highly ordered filament arrays of Ald4ALDH2, a conserved aldehyde dehydrogenase that is highly expressed and phosphorylated during meiosis. Notably, phosphorylation-resistant mutants failed to accumulate filaments, suggesting that phosphorylation regulates context-specific Ald4ALDH2 polymerization. Overall, this proteomic census constitutes a broad resource to guide the exploration of the unique sequence of events underpinning gametogenesis.

Inferring kinase activity from phosphoproteomic data: Tool comparison and recent applications

Aberrant cellular signaling pathways are a hallmark of cancer and other diseases. One of the most important signaling mechanisms involves protein phosphorylation/dephosphorylation. Protein phosphorylation is catalyzed by protein kinases, and over 530 protein kinases have been identified in the human genome. Aberrant kinase activity is one of the drivers of tumorigenesis and cancer progression and results in altered phosphorylation abundance of downstream substrates. Upstream kinase activity can be inferred from the global collection of phosphorylated substrates. Mass spectrometry-based phosphoproteomic experiments nowadays routinely allow identification and quantitation of >10k phosphosites per biological sample. This substrate phosphorylation footprint can be used to infer upstream kinase activities using tools like Kinase Substrate Enrichment Analysis (KSEA), Posttranslational Modification Substrate Enrichment Analysis (PTM-SEA), and Integrative Inferred Kinase Activity Analysis (INKA). Since the topic of kinase activity inference is very active with many new approaches reported in the past 3 years, we would like to give an overview of the field. In this review, an inventory of kinase activity inference tools, their underlying algorithms, statistical frameworks, kinase-substrate databases, and user-friendliness is presented. The most widely-used tools are compared in-depth. Subsequently, recent applications of the tools are described focusing on clinical tissues and hematological samples. Two main application areas for kinase activity inference tools can be discerned. (1) Maximal biological insights can be obtained from large data sets with group comparisons using multiple complementary tools (e.g., PTM-SEA and KSEA or INKA). (2) In the oncology context where personalized treatment requires analysis of single samples, INKA for example, has emerged as tool that can prioritize actionable kinases for targeted inhibition.

Comprehensive multi-omics analysis reveals WEE1 as a synergistic lethal target with hyperthermia through CDK1 super-activation

Hyperthermic intraperitoneal chemotherapy's role in ovarian cancer remains controversial, hindered by limited understanding of hyperthermia-induced tumor cellular changes. This limits developing potent combinatory strategies anchored in hyperthermic intraperitoneal therapy (HIPET). Here, we perform a comprehensive multi-omics study on ovarian cancer cells under hyperthermia, unveiling a distinct molecular panorama, primarily characterized by rapid protein phosphorylation changes. Based on the phospho-signature, we pinpoint CDK1 kinase is hyperactivated during hyperthermia, influencing the global signaling landscape. We observe dynamic, reversible CDK1 activity, causing replication arrest and early mitotic entry post-hyperthermia. Subsequent drug screening shows WEE1 inhibition synergistically destroys cancer cells with hyperthermia. An in-house developed miniaturized device confirms hyperthermia and WEE1 inhibitor combination significantly reduces tumors in vivo. These findings offer additional insights into HIPET, detailing molecular mechanisms of hyperthermia and identifying precise drug combinations for targeted treatment. This research propels the concept of precise hyperthermic intraperitoneal therapy, highlighting its potential against ovarian cancer.

Drug-target identification in COVID-19 disease mechanisms using computational systems biology approaches

Introduction

The COVID-19 Disease Map project is a large-scale community effort uniting 277 scientists from 130 Institutions around the globe. We use high-quality, mechanistic content describing SARS-CoV-2-host interactions and develop interoperable bioinformatic pipelines for novel target identification and drug repurposing.

Methods

Extensive community work allowed an impressive step forward in building interfaces between Systems Biology tools and platforms. Our framework can link biomolecules from omics data analysis and computational modelling to dysregulated pathways in a cell-, tissue- or patient-specific manner. Drug repurposing using text mining and AI-assisted analysis identified potential drugs, chemicals and microRNAs that could target the identified key factors.

Results

Results revealed drugs already tested for anti-COVID-19 efficacy, providing a mechanistic context for their mode of action, and drugs already in clinical trials for treating other diseases, never tested against COVID-19.

Discussion

The key advance is that the proposed framework is versatile and expandable, offering a significant upgrade in the arsenal for virus-host interactions and other complex pathologies.

The genetic and dietary landscape of the muscle insulin signalling network

Metabolic disease is caused by a combination of genetic and environmental factors, yet few studies have examined how these factors influence signal transduction, a key mediator of metabolism. Using mass spectrometry-based phosphoproteomics, we quantified 23,126 phosphosites in skeletal muscle of five genetically distinct mouse strains in two dietary environments, with and without acute in vivo insulin stimulation. Almost half of the insulin-regulated phosphoproteome was modified by genetic background on an ordinary diet, and high-fat high-sugar feeding affected insulin signalling in a strain-dependent manner. Our data revealed coregulated subnetworks within the insulin signalling pathway, expanding our understanding of the pathway's organisation. Furthermore, associating diverse signalling responses with insulin-stimulated glucose uptake uncovered regulators of muscle insulin responsiveness, including the regulatory phosphosite S469 on Pfkfb2, a key activator of glycolysis. Finally, we confirmed the role of glycolysis in modulating insulin action in insulin resistance. Our results underscore the significance of genetics in shaping global signalling responses and their adaptability to environmental changes, emphasising the utility of studying biological diversity with phosphoproteomics to discover key regulatory mechanisms of complex traits.

Mapping protein states and interactions across the tree of life with co-fractionation mass spectrometry

We present CFdb, a harmonized resource of interaction proteomics data from 411 co-fractionation mass spectrometry (CF-MS) datasets spanning 21,703 fractions. Meta-analysis of this resource charts protein abundance, phosphorylation, and interactions throughout the tree of life, including a reference map of the human interactome. We show how large-scale CF-MS data can enhance analyses of individual CF-MS datasets, and exemplify this strategy by mapping the honey bee interactome.

Pharmacological approaches to understanding protein kinase signaling networks

Protein kinases play vital roles in controlling cell behavior, and an array of kinase inhibitors are used successfully for treatment of disease. Typical drug development pipelines involve biological studies to validate a protein kinase target, followed by the identification of small molecules that effectively inhibit this target in cells, animal models, and patients. However, it is clear that protein kinases operate within complex signaling networks. These networks increase the resilience of signaling pathways, which can render cells relatively insensitive to inhibition of a single kinase, and provide the potential for pathway rewiring, which can result in resistance to therapy. It is therefore vital to understand the properties of kinase signaling networks in health and disease so that we can design effective multi-targeted drugs or combinations of drugs. Here, we outline how pharmacological and chemo-genetic approaches can contribute to such knowledge, despite the known low selectivity of many kinase inhibitors. We discuss how detailed profiling of target engagement by kinase inhibitors can underpin these studies; how chemical probes can be used to uncover kinase-substrate relationships, and how these tools can be used to gain insight into the configuration and function of kinase signaling networks.

Proteomic and genetic analyses of influenza A viruses identify pan-viral host targets

Influenza A Virus (IAV) is a recurring respiratory virus with limited availability of antiviral therapies. Understanding host proteins essential for IAV infection can identify targets for alternative host-directed therapies (HDTs). Using affinity purification-mass spectrometry and global phosphoproteomic and protein abundance analyses using three IAV strains (pH1N1, H3N2, H5N1) in three human cell types (A549, NHBE, THP-1), we map 332 IAV-human protein-protein interactions and identify 13 IAV-modulated kinases. Whole exome sequencing of patients who experienced severe influenza reveals several genes, including scaffold protein AHNAK, with predicted loss-of-function variants that are also identified in our proteomic analyses. Of our identified host factors, 54 significantly alter IAV infection upon siRNA knockdown, and two factors, AHNAK and coatomer subunit COPB1, are also essential for productive infection by SARS-CoV-2. Finally, 16 compounds targeting our identified host factors suppress IAV replication, with two targeting CDK2 and FLT3 showing pan-antiviral activity across influenza and coronavirus families. This study provides a comprehensive network model of IAV infection in human cells, identifying functional host targets for pan-viral HDT.

Defining blood-induced microglia functions in neurodegeneration through multiomic profiling

Blood protein extravasation through a disrupted blood-brain barrier and innate immune activation are hallmarks of neurological diseases and emerging therapeutic targets. However, how blood proteins polarize innate immune cells remains largely unknown. Here, we established an unbiased blood-innate immunity multiomic and genetic loss-of-function pipeline to define the transcriptome and global phosphoproteome of blood-induced innate immune polarization and its role in microglia neurotoxicity. Blood induced widespread microglial transcriptional changes, including changes involving oxidative stress and neurodegenerative genes. Comparative functional multiomics showed that blood proteins induce distinct receptor-mediated transcriptional programs in microglia and macrophages, such as redox, type I interferon and lymphocyte recruitment. Deletion of the blood coagulation factor fibrinogen largely reversed blood-induced microglia neurodegenerative signatures. Genetic elimination of the fibrinogen-binding motif to CD11b in Alzheimer's disease mice reduced microglial lipid metabolism and neurodegenerative signatures that were shared with autoimmune-driven neuroinflammation in multiple sclerosis mice. Our data provide an interactive resource for investigation of the immunology of blood proteins that could support therapeutic targeting of microglia activation by immune and vascular signals.

A Chemical Proteomics Approach to Discover Regulators of Innate Immune Signaling

Innate immune pathways are tightly regulated to balance an appropriate response to infectious agents and tolerable levels of inflammation. Dysregulation of innate immune pathways can lead to severe autoinflammatory disorders or susceptibility to infections. Here, we aimed to identify kinases in common cellular pathways that regulate innate immune pathways by combining small-scale kinase inhibitor screening with quantitative proteomics. We found that inhibitors of kinases ATM, ATR, AMPK, and PLK1 reduced the induction of interferon-stimulated gene expression in response to innate immune pathway activation by poly(I:C) transfection. However, siRNA depletion of these kinases did not validate findings with kinase inhibitors, suggesting that off-target effects may explain their activities. We mapped the effects of kinase inhibitors to various stages in innate immune pathways. Determining the mechanisms by which kinase inhibitors antagonize these pathways may illuminate novel mechanisms of innate immune pathway control.

Recent advances in kinase signaling network profiling by mass spectrometry

Mass spectrometry-based phosphoproteomics is currently the leading methodology for the study of global kinase signaling. The scientific community is continuously releasing technological improvements for sensitive and fast identification of phosphopeptides, and their accurate quantification. To interpret large-scale phosphoproteomics data, numerous bioinformatic resources are available that help understanding kinase network functional role in biological systems upon perturbation. Some of these resources are databases of phosphorylation sites, protein kinases and phosphatases; others are bioinformatic algorithms to infer kinase activity, predict phosphosite functional relevance and visualize kinase signaling networks. In this review, we present the latest experimental and bioinformatic tools to profile protein kinase signaling networks and provide examples of their application in biomedicine.

Principles of phosphoproteomics and applications in cancer research

Phosphorylation constitutes the most common and best-studied regulatory post-translational modification in biological systems and archetypal signalling pathways driven by protein and lipid kinases are disrupted in essentially all cancer types. Thus, the study of the phosphoproteome stands to provide unique biological information on signalling pathway activity and on kinase network circuitry that is not captured by genetic or transcriptomic technologies. Here, we discuss the methods and tools used in phosphoproteomics and highlight how this technique has been used, and can be used in the future, for cancer research. Challenges still exist in mass spectrometry phosphoproteomics and in the software required to provide biological information from these datasets. Nevertheless, improvements in mass spectrometers with enhanced scan rates, separation capabilities and sensitivity, in biochemical methods for sample preparation and in computational pipelines are enabling an increasingly deep analysis of the phosphoproteome, where previous bottlenecks in data acquisition, processing and interpretation are being relieved. These powerful hardware and algorithmic innovations are not only providing exciting new mechanistic insights into tumour biology, from where new drug targets may be derived, but are also leading to the discovery of phosphoproteins as mediators of drug sensitivity and resistance and as classifiers of disease subtypes. These studies are, therefore, uncovering phosphoproteins as a new generation of disruptive biomarkers to improve personalised anti-cancer therapies.

Pan-Cancer landscape of protein activities identifies drivers of signalling dysregulation and patient survival

Genetic alterations in cancer cells trigger oncogenic transformation, a process largely mediated by the dysregulation of kinase and transcription factor (TF) activities. While the mutational profiles of thousands of tumours have been extensively characterised, the measurements of protein activities have been technically limited until recently. We compiled public data of matched genomics and (phospho)proteomics measurements for 1,110 tumours and 77 cell lines that we used to estimate activity changes in 218 kinases and 292 TFs. Co-regulation of kinase and TF activities reflects previously known regulatory relationships and allows us to dissect genetic drivers of signalling changes in cancer. We find that loss-of-function mutations are not often associated with the dysregulation of downstream targets, suggesting frequent compensatory mechanisms. Finally, we identified the activities most differentially regulated in cancer subtypes and showed how these can be linked to differences in patient survival. Our results provide broad insights into the dysregulation of protein activities in cancer and their contribution to disease severity.

Phosphoproteomics reveals rewiring of the insulin signaling network and multi-nodal defects in insulin resistance

The failure of metabolic tissues to appropriately respond to insulin ("insulin resistance") is an early marker in the pathogenesis of type 2 diabetes. Protein phosphorylation is central to the adipocyte insulin response, but how adipocyte signaling networks are dysregulated upon insulin resistance is unknown. Here we employ phosphoproteomics to delineate insulin signal transduction in adipocyte cells and adipose tissue. Across a range of insults causing insulin resistance, we observe a marked rewiring of the insulin signaling network. This includes both attenuated insulin-responsive phosphorylation, and the emergence of phosphorylation uniquely insulin-regulated in insulin resistance. Identifying dysregulated phosphosites common to multiple insults reveals subnetworks containing non-canonical regulators of insulin action, such as MARK2/3, and causal drivers of insulin resistance. The presence of several bona fide GSK3 substrates among these phosphosites led us to establish a pipeline for identifying context-specific kinase substrates, revealing widespread dysregulation of GSK3 signaling. Pharmacological inhibition of GSK3 partially reverses insulin resistance in cells and tissue explants. These data highlight that insulin resistance is a multi-nodal signaling defect that includes dysregulated MARK2/3 and GSK3 activity.

Computational systems approach towards phosphoproteomics and their downstream regulation

Protein phosphorylation plays an essential role in modulating cell signalling and its downstream transcriptional and translational regulations. Until recently, protein phosphorylation has been studied mostly using low-throughput biochemical assays. The advancement of mass spectrometry (MS)-based phosphoproteomics transformed the field by enabling measurement of proteome-wide phosphorylation events, where tens of thousands of phosphosites are routinely identified and quantified in an experiment. This has brought a significant challenge in analysing large-scale phosphoproteomic data, making computational methods and systems approaches integral parts of phosphoproteomics. Previous works have primarily focused on reviewing the experimental techniques in MS-based phosphoproteomics, yet a systematic survey of the computational landscape in this field is still missing. Here, we review computational methods and tools, and systems approaches that have been developed for phosphoproteomics data analysis. We categorise them into four aspects including data processing, functional analysis, phosphoproteome annotation and their integration with other omics, and in each aspect, we discuss the key methods and example studies. Lastly, we highlight some of the potential research directions on which future work would make a significant contribution to this fast-growing field. We hope this review provides a useful snapshot of the field of computational systems phosphoproteomics and stimulates new research that drives future development.

SARS-CoV-2 replication in airway epithelia requires motile cilia and microvillar reprogramming

How SARS-CoV-2 penetrates the airway barrier of mucus and periciliary mucins to infect nasal epithelium remains unclear. Using primary nasal epithelial organoid cultures, we found that the virus attaches to motile cilia via the ACE2 receptor. SARS-CoV-2 traverses the mucus layer, using motile cilia as tracks to access the cell body. Depleting cilia blocks infection for SARS-CoV-2 and other respiratory viruses. SARS-CoV-2 progeny attach to airway microvilli 24 h post-infection and trigger formation of apically extended and highly branched microvilli that organize viral egress from the microvilli back into the mucus layer, supporting a model of virus dispersion throughout airway tissue via mucociliary transport. Phosphoproteomics and kinase inhibition reveal that microvillar remodeling is regulated by p21-activated kinases (PAK). Importantly, Omicron variants bind with higher affinity to motile cilia and show accelerated viral entry. Our work suggests that motile cilia, microvilli, and mucociliary-dependent mucus flow are critical for efficient virus replication in nasal epithelia.

Resistance to ATR Inhibitors Is Mediated by Loss of the Nonsense-Mediated Decay Factor UPF2

Over one million cases of gastric cancer are diagnosed each year globally, and the metastatic disease continues to have a poor prognosis. A significant proportion of gastric tumors have defects in the DNA damage response pathway, creating therapeutic opportunities through synthetic lethal approaches. Several small-molecule inhibitors of ATR, a key regulator of the DNA damage response, are now in clinical development as targeted agents for gastric cancer. Here, we performed a large-scale CRISPR interference screen to discover genetic determinants of response and resistance to ATR inhibitors (ATRi) in gastric cancer cells. Among the top hits identified as mediators of ATRi response were UPF2 and other components of the nonsense-mediated decay (NMD) pathway. Loss of UPF2 caused ATRi resistance across multiple gastric cancer cell lines. Global proteomic, phosphoproteomic, and transcriptional profiling experiments revealed that cell-cycle progression and DNA damage responses were altered in UPF2-mutant cells. Further studies demonstrated that UPF2-depleted cells failed to accumulate in G1 following treatment with ATRi. UPF2 loss also reduced transcription-replication collisions, which has previously been associated with ATRi response, thereby suggesting a possible mechanism of resistance. Our results uncover a novel role for NMD factors in modulating response to ATRi in gastric cancer, highlighting a previously unknown mechanism of resistance that may inform the clinical use of these drugs.

SIGNIFICANCE

Loss of NMD proteins promotes resistance to ATR inhibitors in gastric cancer cells, which may provide a combination of therapeutic targets and biomarkers to improve the clinical utility of these drugs.

Systematic discovery of biomolecular condensate-specific protein phosphorylation

Reversible protein phosphorylation is an important mechanism for regulating (dis)assembly of biomolecular condensates. However, condensate-specific phosphosites remain largely unknown, thereby limiting our understanding of the underlying mechanisms. Here, we combine solubility proteome profiling with phosphoproteomics to quantitatively map several hundred phosphosites enriched in either soluble or condensate-bound protein subpopulations, including a subset of phosphosites modulating protein-RNA interactions. We show that multi-phosphorylation of the C-terminal disordered segment of heteronuclear ribonucleoprotein A1 (HNRNPA1), a key RNA-splicing factor, reduces its ability to locate to nuclear clusters. For nucleophosmin 1 (NPM1), an essential nucleolar protein, we show that phosphorylation of S254 and S260 is crucial for lowering its partitioning to the nucleolus and additional phosphorylation of distal sites enhances its retention in the nucleoplasm. These phosphorylation events decrease RNA and protein interactions of NPM1 to regulate its condensation. Our dataset is a rich resource for systematically uncovering the phosphoregulation of biomolecular condensates.

KSTAR: An algorithm to predict patient-specific kinase activities from phosphoproteomic data

Kinase inhibitors as targeted therapies have played an important role in improving cancer outcomes. However, there are still considerable challenges, such as resistance, non-response, patient stratification, polypharmacology, and identifying combination therapy where understanding a tumor kinase activity profile could be transformative. Here, we develop a graph- and statistics-based algorithm, called KSTAR, to convert phosphoproteomic measurements of cells and tissues into a kinase activity score that is generalizable and useful for clinical pipelines, requiring no quantification of the phosphorylation sites. In this work, we demonstrate that KSTAR reliably captures expected kinase activity differences across different tissues and stimulation contexts, allows for the direct comparison of samples from independent experiments, and is robust across a wide range of dataset sizes. Finally, we apply KSTAR to clinical breast cancer phosphoproteomic data and find that there is potential for kinase activity inference from KSTAR to complement the current clinical diagnosis of HER2 status in breast cancer patients.

Functional Characterization of Co-Phosphorylation Networks

MOTIVATION

Protein phosphorylation is a ubiquitous regulatory mechanism that plays a central role in cellular signaling. According to recent estimates, up to 70% of human proteins can be phosphorylated. Therefore, characterization of phosphorylation dynamics is critical for understanding a broad range of biological and biochemical processes. Technologies based on mass spectrometry are rapidly advancing to meet the needs for high-throughput screening of phosphorylation. These technologies enable untargeted quantification of thousands of phosphorylation sites in a given sample. Many labs are already utilizing these technologies to comprehensively characterize signaling landscapes by examining perturbations with drugs and knockdown approaches, or by assessing diverse phenotypes in cancers, neuro-degerenational diseases, infectious diseases, and normal development.

RESULTS

We comprehensively investigate the concept of "co-phosphorylation", defined as the correlated phosphorylation of a pair of phosphosites across various biological states. We integrate nine publicly available phosphoproteomics datasets for various diseases (including breast cancer, ovarian cancer and Alzheimer's disease) and utilize functional data related to sequence, evolutionary histories, kinase annotations, and pathway annotations to investigate the functional relevance of co-phosphorylation. Our results across a broad range of studies consistently show that functionally associated sites tend to exhibit significant positive or negative co-phosphorylation. Specifically, we show that co-phosphorylation can be used to predict with high precision the sites that are on the same pathway or that are targeted by the same kinase. Overall, these results establish co-phosphorylation as a useful resource for analyzing phosphoproteins in a network context, which can help extend our knowledge on cellular signaling and its dysregulation.

Accurate, high-coverage assignment of in vivo protein kinases to phosphosites from in vitro phosphoproteomic specificity data

Phosphoproteomic experiments routinely observe thousands of phosphorylation sites. To understand the intracellular signaling processes that generated this data, one or more causal protein kinases must be assigned to each phosphosite. However, limited knowledge of kinase specificity typically restricts assignments to a small subset of a kinome. Starting from a statistical model of a high-throughput, in vitro kinase-substrate assay, I have developed an approach to high-coverage, multi-label kinase-substrate assignment called IV-KAPhE (“In vivo-Kinase Assignment for Phosphorylation Evidence”). Tested on human data, IV-KAPhE outperforms other methods of similar scope. Such computational methods generally predict a densely connected kinase-substrate network, with most sites targeted by multiple kinases, pointing either to unaccounted-for biochemical constraints or significant cross-talk and signaling redundancy. I show that such predictions can potentially identify biased kinase-site misannotations within families of closely related kinase isozymes and they provide a robust basis for kinase activity analysis.

Proteins can pass around information inside cells about changes in the environment. This process, called intracellular signaling, helps to trigger appropriate cellular responses to environmental changes. One of the main ways information is passed to proteins is through chemical “tagging,” called phosphorylation, by enzymes called protein kinases. We can measure the phosphorylation state of practically all proteins in a cell at any moment. Starting from known cases of phosphorylation by a kinase, many computational methods have been developed to predict if the kinase might tag a certain spot on another protein or if an observed tag was attached by the kinase, with different models for each kinase. I have developed a new method that instead uses a single model to assign one or more kinases to each observed tag, built from the latest large-scale experimental data. This change in focus and unbiased training data allows my method to be significantly more accurate than past methods. I also explored useful applications for my method. For example, I used it to show that much of our knowledge about which kinase is responsible for each tag is probably inaccurately biased towards the commonly studied ones.

Global post-translational modification profiling of HIV-1-infected cells reveals mechanisms of host cellular pathway remodeling

Viruses must effectively remodel host cellular pathways to replicate and evade immune defenses, and they must do so with limited genomic coding capacity. Targeting post-translational modification (PTM) pathways provides a mechanism by which viruses can broadly and rapidly transform a hostile host environment into a hospitable one. We use mass spectrometry-based proteomics to quantify changes in protein abundance and two PTM types-phosphorylation and ubiquitination-in response to HIV-1 infection with viruses harboring targeted deletions of a subset of HIV-1 genes. PTM analysis reveals a requirement for Aurora kinase activity in HIV-1 infection and identified putative substrates of a phosphatase that is degraded during infection. Finally, we demonstrate that the HIV-1 Vpr protein inhibits histone H1 ubiquitination, leading to defects in DNA repair.

decoupleR: ensemble of computational methods to infer biological activities from omics data

Summary

Many methods allow us to extract biological activities from omics data using information from prior knowledge resources, reducing the dimensionality for increased statistical power and better interpretability. Here, we present decoupleR, a Bioconductor and Python package containing computational methods to extract these activities within a unified framework. decoupleR allows us to flexibly run any method with a given resource, including methods that leverage mode of regulation and weights of interactions, which are not present in other frameworks. Moreover, it leverages OmniPath, a meta-resource comprising over 100 databases of prior knowledge. Using decoupleR, we evaluated the performance of methods on transcriptomic and phospho-proteomic perturbation experiments. Our findings suggest that simple linear models and the consensus score across top methods perform better than other methods at predicting perturbed regulators.

Availability and implementation

decoupleR's open-source code is available in Bioconductor (https://www.bioconductor.org/packages/release/bioc/html/decoupleR.html) for R and in GitHub (https://github.com/saezlab/decoupler-py) for Python. The code to reproduce the results is in GitHub (https://github.com/saezlab/decoupleR_manuscript) and the data in Zenodo (https://zenodo.org/record/5645208).

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

CFIm-mediated alternative polyadenylation remodels cellular signaling and miRNA biogenesis

The mammalian cleavage factor I (CFIm) has been implicated in alternative polyadenylation (APA) in a broad range of contexts, from cancers to learning deficits and parasite infections. To determine how the CFIm expression levels are translated into these diverse phenotypes, we carried out a multi-omics analysis of cell lines in which the CFIm25 (NUDT21) or CFIm68 (CPSF6) subunits were either repressed by siRNA-mediated knockdown or over-expressed from stably integrated constructs. We established that >800 genes undergo coherent APA in response to changes in CFIm levels, and they cluster in distinct functional classes related to protein metabolism. The activity of the ERK pathway traces the CFIm concentration, and explains some of the fluctuations in cell growth and metabolism that are observed upon CFIm perturbations. Furthermore, multiple transcripts encoding proteins from the miRNA pathway are targets of CFIm-dependent APA. This leads to an increased biogenesis and repressive activity of miRNAs at the same time as some 3′ UTRs become shorter and presumably less sensitive to miRNA-mediated repression. Our study provides a first systematic assessment of a core set of APA targets that respond coherently to changes in CFIm protein subunit levels (CFIm25/CFIm68). We describe the elicited signaling pathways downstream of CFIm, which improve our understanding of the key role of CFIm in integrating RNA processing with other cellular activities.

Evolution of enhanced innate immune evasion by SARS-CoV-2

The emergence of SARS-CoV-2 variants of concern suggests viral adaptation to enhance human-to-human transmission1,2. Although much effort has focused on the characterization of changes in the spike protein in variants of concern, mutations outside of spike are likely to contribute to adaptation. Here, using unbiased abundance proteomics, phosphoproteomics, RNA sequencing and viral replication assays, we show that isolates of the Alpha (B.1.1.7) variant3 suppress innate immune responses in airway epithelial cells more effectively than first-wave isolates. We found that the Alpha variant has markedly increased subgenomic RNA and protein levels of the nucleocapsid protein (N), Orf9b and Orf6-all known innate immune antagonists. Expression of Orf9b alone suppressed the innate immune response through interaction with TOM70, a mitochondrial protein that is required for activation of the RNA-sensing adaptor MAVS. Moreover, the activity of Orf9b and its association with TOM70 was regulated by phosphorylation. We propose that more effective innate immune suppression, through enhanced expression of specific viral antagonist proteins, increases the likelihood of successful transmission of the Alpha variant, and may increase in vivo replication and duration of infection4. The importance of mutations outside the spike coding region in the adaptation of SARS-CoV-2 to humans is underscored by the observation that similar mutations exist in the N and Orf9b regulatory regions of the Delta and Omicron variants.

SARS-CoV-2 infects human pancreatic β cells and elicits β cell impairment

Emerging evidence points toward an intricate relationship between the pandemic of coronavirus disease 2019 (COVID-19) and diabetes. While preexisting diabetes is associated with severe COVID-19, it is unclear whether COVID-19 severity is a cause or consequence of diabetes. To mechanistically link COVID-19 to diabetes, we tested whether insulin-producing pancreatic β cells can be infected by SARS-CoV-2 and cause β cell depletion. We found that the SARS-CoV-2 receptor, ACE2, and related entry factors (TMPRSS2, NRP1, and TRFC) are expressed in β cells, with selectively high expression of NRP1. We discovered that SARS-CoV-2 infects human pancreatic β cells in patients who succumbed to COVID-19 and selectively infects human islet β cells in vitro. We demonstrated that SARS-CoV-2 infection attenuates pancreatic insulin levels and secretion and induces β cell apoptosis, each rescued by NRP1 inhibition. Phosphoproteomic pathway analysis of infected islets indicates apoptotic β cell signaling, similar to that observed in type 1 diabetes (T1D). In summary, our study shows SARS-CoV-2 can directly induce β cell killing.

Evolution of enhanced innate immune evasion by the SARS-CoV-2 B.1.1.7 UK variant

Emergence of SARS-CoV-2 variants, including the globally successful B.1.1.7 lineage, suggests viral adaptations to host selective pressures resulting in more efficient transmission. Although much effort has focused on Spike adaptation for viral entry and adaptive immune escape, B.1.1.7 mutations outside Spike likely contribute to enhance transmission. Here we used unbiased abundance proteomics, phosphoproteomics, mRNA sequencing and viral replication assays to show that B.1.1.7 isolates more effectively suppress host innate immune responses in airway epithelial cells. We found that B.1.1.7 isolates have dramatically increased subgenomic RNA and protein levels of Orf9b and Orf6, both known innate immune antagonists. Expression of Orf9b alone suppressed the innate immune response through interaction with TOM70, a mitochondrial protein required for RNA sensing adaptor MAVS activation, and Orf9b binding and activity was regulated via phosphorylation. We conclude that B.1.1.7 has evolved beyond the Spike coding region to more effectively antagonise host innate immune responses through upregulation of specific subgenomic RNA synthesis and increased protein expression of key innate immune antagonists. We propose that more effective innate immune antagonism increases the likelihood of successful B.1.1.7 transmission, and may increase in vivo replication and duration of infection.

Combining Mass Spectrometry-Based Phosphoproteomics with a Network-Based Approach to Reveal FLT3-Dependent Mechanisms of Chemoresistance

FLT3 mutations are the most frequently identified genetic alterations in acute myeloid leukemia (AML) and are associated with poor clinical outcome, relapse and chemotherapeutic resistance. Elucidating the molecular mechanisms underlying FLT3-dependent pathogenesis and drug resistance is a crucial goal of biomedical research. Given the complexity and intricacy of protein signaling networks, deciphering the molecular basis of FLT3-driven drug resistance requires a systems approach. Here we discuss how the recent advances in mass spectrometry (MS)-based (phospho) proteomics and multiparametric analysis accompanied by emerging computational approaches offer a platform to obtain and systematically analyze cell-specific signaling networks and to identify new potential therapeutic targets.

PHONEMeS: Efficient Modeling of Signaling Networks Derived from Large-Scale Mass Spectrometry Data

Post-translational modifications of proteins play an important role in the regulation of cellular processes. The mass spectrometry analysis of proteome modifications offers huge potential for the study of how protein inhibitors affect the phosphosignaling mechanisms inside the cells. We have recently proposed PHONEMeS, a method that uses high-content shotgun phosphoproteomic data to build logical network models of signal perturbation flow. However, in its original implementation, PHONEMeS was computationally demanding and was only used to model signaling in a perturbation context. We have reformulated PHONEMeS as an Integer Linear Program (ILP) that is orders of magnitude more efficient than the original one. We have also expanded the scenarios that can be analyzed. PHONEMeS can model data upon perturbation on not only a known target but also deregulated pathways upstream and downstream of any set of deregulated kinases. Finally, PHONEMeS can now analyze data sets with multiple time points, which helps us to obtain better insight into the dynamics of the propagation of signals. We illustrate the value of the new approach on various data sets of medical relevance, where we shed light on signaling mechanisms and drug modes of action.

Robust inference of kinase activity using functional networks

Mass spectrometry enables high-throughput screening of phosphoproteins across a broad range of biological contexts. When complemented by computational algorithms, phospho-proteomic data allows the inference of kinase activity, facilitating the identification of dysregulated kinases in various diseases including cancer, Alzheimer’s disease and Parkinson’s disease. To enhance the reliability of kinase activity inference, we present a network-based framework, RoKAI, that integrates various sources of functional information to capture coordinated changes in signaling. Through computational experiments, we show that phosphorylation of sites in the functional neighborhood of a kinase are significantly predictive of its activity. The incorporation of this knowledge in RoKAI consistently enhances the accuracy of kinase activity inference methods while making them more robust to missing annotations and quantifications. This enables the identification of understudied kinases and will likely lead to the development of novel kinase inhibitors for targeted therapy of many diseases. RoKAI is available as web-based tool at http://rokai.io.

Kinases drive fundamental changes in cell state, but predicting kinase activity based on substrate-level changes can be challenging. Here the authors introduce a computational framework that utilizes similarities between substrates to robustly infer kinase activity.

Proteomics-based interrogation of the kinome and its implications for precision oncology

The identification of specific protein kinases as oncogenic drivers in a variety of cancer types, coupled with the clinical success of particular kinase-directed targeted therapies, has cemented the human kinome as an attractive source of "actionable" targets for cancer therapy. However, "mining" of the human kinome for precision oncology applications has yet to yield its full potential. This reflects a variety of issues, including oncogenic kinase dysregulation at levels not detectable by genomic sequencing and the uncharacterized nature of a considerable fraction of the kinome. In addition, selective therapeutic targeting of specific kinases requires efficient mapping of total kinome space impacted by candidate small molecule drugs. Fortunately, recent developments in proteomics techniques, particularly in mass spectrometry-based phosphoproteomics and kinomics, provide the necessary technology platforms to address these impediments. Moreover, initiatives such as the Clinical Proteomic Tumour Analysis Consortium have enabled the generation, deposition and integration of genomic, transcriptomic and (phospho)proteomic data for many cancer types, providing unprecedented insights into oncogenic kinases and cancer cell signalling generally. These multi-omic data are identifying novel therapeutic targets, highlighting opportunities for drug re-purposing, and helping assign optimal therapies to specific tumour subtypes, heralding a new era of "enhanced" precision oncology.

Proteomic Approaches to Study SARS-CoV-2 Biology and COVID-19 Pathology

The novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of coronavirus disease 2019 (COVID-19), was declared a pandemic infection in March 2020. As of December 2020, two COVID-19 vaccines have been authorized for emergency use by the U.S. Food and Drug Administration, but there are no effective drugs to treat COVID-19, and pandemic mitigation efforts like physical distancing have had acute social and economic consequences. In this perspective, we discuss how the proteomic research community can leverage technologies and expertise to address the pandemic by investigating four key areas of study in SARS-CoV-2 biology. Specifically, we discuss how (1) mass spectrometry-based structural techniques can overcome limitations and complement traditional structural approaches to inform the dynamic structure of SARS-CoV-2 proteins, complexes, and virions; (2) virus-host protein-protein interaction mapping can identify the cellular machinery required for SARS-CoV-2 replication; (3) global protein abundance and post-translational modification profiling can characterize signaling pathways that are rewired during infection; and (4) proteomic technologies can aid in biomarker identification, diagnostics, and drug development in order to monitor COVID-19 pathology and investigate treatment strategies. Systems-level high-throughput capabilities of proteomic technologies can yield important insights into SARS-CoV-2 biology that are urgently needed during the pandemic, and more broadly, can inform coronavirus virology and host biology.

Accurate prediction of kinase-substrate networks using knowledge graphs

Phosphorylation of specific substrates by protein kinases is a key control mechanism for vital cell-fate decisions and other cellular processes. However, discovering specific kinase-substrate relationships is time-consuming and often rather serendipitous. Computational predictions alleviate these challenges, but the current approaches suffer from limitations like restricted kinome coverage and inaccuracy. They also typically utilise only local features without reflecting broader interaction context. To address these limitations, we have developed an alternative predictive model. It uses statistical relational learning on top of phosphorylation networks interpreted as knowledge graphs, a simple yet robust model for representing networked knowledge. Compared to a representative selection of six existing systems, our model has the highest kinome coverage and produces biologically valid high-confidence predictions not possible with the other tools. Specifically, we have experimentally validated predictions of previously unknown phosphorylations by the LATS1, AKT1, PKA and MST2 kinases in human. Thus, our tool is useful for focusing phosphoproteomic experiments, and facilitates the discovery of new phosphorylation reactions. Our model can be accessed publicly via an easy-to-use web interface (LinkPhinder).

LinkPhinder is a new approach to prediction of protein signalling networks based on kinase-substrate relationships that outperforms existing approaches. Phosphorylation networks govern virtually all fundamental biochemical processes in cells, and thus have moved into the centre of interest in biology, medicine and drug development. Fundamentally different from current approaches, LinkPhinder is inherently network-based and makes use of the most recent AI developments. We represent existing phosphorylation data as knowledge graphs, a format for large-scale and robust knowledge representation. Training a link prediction model on such a structure leads to novel, biologically valid phosphorylation network predictions that cannot be made with competing tools. Thus our new conceptual approach can lead to establishing a new niche of AI applications in computational biology.

The Global Phosphorylation Landscape of SARS-CoV-2 Infection

The causative agent of the coronavirus disease 2019 (COVID-19) pandemic, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has infected millions and killed hundreds of thousands of people worldwide, highlighting an urgent need to develop antiviral therapies. Here we present a quantitative mass spectrometry-based phosphoproteomics survey of SARS-CoV-2 infection in Vero E6 cells, revealing dramatic rewiring of phosphorylation on host and viral proteins. SARS-CoV-2 infection promoted casein kinase II (CK2) and p38 MAPK activation, production of diverse cytokines, and shutdown of mitotic kinases, resulting in cell cycle arrest. Infection also stimulated a marked induction of CK2-containing filopodial protrusions possessing budding viral particles. Eighty-seven drugs and compounds were identified by mapping global phosphorylation profiles to dysregulated kinases and pathways. We found pharmacologic inhibition of the p38, CK2, CDK, AXL, and PIKFYVE kinases to possess antiviral efficacy, representing potential COVID-19 therapies.

Using phosphoproteomics data to understand cellular signaling: a comprehensive guide to bioinformatics resources

Mass spectrometry-based phosphoproteomics is becoming an essential methodology for the study of global cellular signaling. Numerous bioinformatics resources are available to facilitate the translation of phosphopeptide identification and quantification results into novel biological and clinical insights, a critical step in phosphoproteomics data analysis. These resources include knowledge bases of kinases and phosphatases, phosphorylation sites, kinase inhibitors, and sequence variants affecting kinase function, and bioinformatics tools that can predict phosphorylation sites in addition to the kinase that phosphorylates them, infer kinase activity, and predict the effect of mutations on kinase signaling. However, these resources exist in silos and it is challenging to select among multiple resources with similar functions. Therefore, we put together a comprehensive collection of resources related to phosphoproteomics data interpretation, compared the use of tools with similar functions, and assessed the usability from the standpoint of typical biologists or clinicians. Overall, tools could be improved by standardization of enzyme names, flexibility of data input and output format, consistent maintenance, and detailed manuals.

Substrate-based kinase activity inference identifies MK2 as driver of colitis

Inflammatory bowel disease (IBD) is a chronic and debilitating disorder that has few treatment options due to a lack of comprehensive understanding of its molecular pathogenesis. We used multiplexed mass spectrometry to collect high-content information on protein phosphorylation in two different mouse models of IBD. Because the biological function of the vast majority of phosphorylation sites remains unknown, we developed Substrate-based Kinase Activity Inference (SKAI), a methodology to infer kinase activity from phosphoproteomic data. This approach draws upon prior knowledge of kinase-substrate interactions to construct custom lists of kinases and their respective substrate sites, termed kinase-substrate sets that employ prior knowledge across organisms. This expansion as much as triples the amount of prior knowledge available. We then used these sets within the Gene Set Enrichment Analysis framework to infer kinase activity based on increased or decreased phosphorylation of its substrates in a dataset. When applied to the phosphoproteomic datasets from the two mouse models, SKAI predicted largely non-overlapping kinase activation profiles. These results suggest that chronic inflammation may arise through activation of largely divergent signaling networks. However, the one kinase inferred to be activated in both mouse models was mitogen-activated protein kinase-activated protein kinase 2 (MAPKAPK2 or MK2), a serine/threonine kinase that functions downstream of p38 stress-activated mitogen-activated protein kinase. Treatment of mice with active colitis with ATI450, an orally bioavailable small molecule inhibitor of the MK2 pathway, reduced inflammatory signaling in the colon and alleviated the clinical and histological features of inflammation. These studies establish MK2 as a therapeutic target in IBD and identify ATI450 as a potential therapy for the disease.

Reconstructing kinase network topologies from phosphoproteomics data reveals cancer-associated rewiring

Understanding how oncogenic mutations rewire regulatory-protein networks is important for rationalizing the mechanisms of oncogenesis and for individualizing anticancer treatments. We report a chemical phosphoproteomics method to elucidate the topology of kinase-signaling networks in mammalian cells. We identified >6,000 protein phosphorylation sites that can be used to infer >1,500 kinase-kinase interactions and devised algorithms that can reconstruct kinase network topologies from these phosphoproteomics data. Application of our methods to primary acute myeloid leukemia and breast cancer tumors quantified the relationship between kinase expression and activity, and enabled the identification of hitherto unknown kinase network topologies associated with drug-resistant phenotypes or specific genetic mutations. Using orthogonal methods we validated that PIK3CA wild-type cells adopt MAPK-dependent circuitries in breast cancer cells and that the kinase TTK is important in acute myeloid leukemia. Our phosphoproteomic signatures of network circuitry can identify kinase topologies associated with both phenotypes and genotypes of cancer cells.

Phosphoproteomics Enables Molecular Subtyping and Nomination of Kinase Candidates for Individual Patients of Diffuse-Type Gastric Cancer

The diffuse-type gastric cancer (DGC) constitutes a subgroup of gastric cancer with poor prognosis and no effective molecular therapies. Here, we report a phosphoproteomic landscape of DGC derived from 83 tumors together with their nearby tissues. Based on phosphorylation, DGC could be classified into three molecular subtypes with distinct overall survival (OS) and chemosensitivity. We identified 16 kinases whose activities were associated with poor OS. These activated kinases covered several cancer hallmark pathways, with the MTOR signaling network being the most frequently activated. We proposed a patient-specific strategy based on the hierarchy of clinically actionable kinases for prioritization of kinases for further clinical evaluation. Our global data analysis indicates that in addition to finding activated kinase pathways in DGC, large-scale phosphoproteomics could be used to classify DGCs into subtypes that are associated with distinct clinical outcomes as well as nomination of kinase targets that may be inhibited for cancer treatments.

Footprint-based functional analysis of multiomic data

Omic technologies allow us to generate extensive data, including transcriptomic, proteomic, phosphoproteomic and metabolomic. These data can be used to study signal transduction, gene regulation and metabolism. In this review, we summarise resources and methods to analysis these types of data. We focus on methods developed to recover functional insights using footprints. Footprints are signatures defined by the effect of molecules or processes of interest. They integrate information from multiple measurements whose abundances are under the influence of a common regulator. For example, transcripts controlled by a transcription factor or peptides phosphorylated by a kinase. Footprints can also be generalised across multiple types of omic data. Thus, we also present methods to integrate multiple types of omic data and features (such as the ones derived from footprints) together. We highlight some examples of studies that leverage such approaches to discover new biological mechanisms.

•

Functional information on signalling pathways, metabolism and gene regulation can be found across multiple types of omic data.

•

One way to extract such information is to consider these data as the footprint of the activity of enzymes and pathways.

•

Information on enzyme/pathway activities and omic data can be integrated together to contextualise multi-scale networks.

•

Such an approach can lead to the discovery of regulatory events spanning across multiple biological processes.

Proteomics and phosphoproteomics in precision medicine: applications and challenges

Recent advances in proteomics allow the accurate measurement of abundances for thousands of proteins and phosphoproteins from multiple samples in parallel. Therefore, for the first time, we have the opportunity to measure the proteomic profiles of thousands of patient samples or disease model cell lines in a systematic way, to identify the precise underlying molecular mechanism and discover personalized biomarkers, networks and treatments. Here, we review examples of successful use of proteomics and phosphoproteomics data sets in as well as their integration other omics data sets with the aim of precision medicine. We will discuss the bioinformatics challenges posed by the generation, analysis and integration of such large data sets and present potential reasons why proteomics profiling and biomarkers are not currently widely used in the clinical setting. We will finally discuss ways to contribute to the better use of proteomics data in precision medicine and the clinical setting.