Genetic wiring maps of single-cell protein states reveal an off-switch for GPCR signalling.

As key executers of biological functions, the activity and abundance of proteins are subjected to extensive regulation. Deciphering the genetic architecture underlying this regulation is critical for understanding cellular signalling events and responses to environmental cues. Using random mutagenesis in haploid human cells, we apply a sensitive approach to directly couple genomic mutations to protein measurements in individual cells. Here we use this to examine a suite of cellular processes, such as transcriptional induction, regulation of protein abundance and splicing, signalling cascades (mitogen-activated protein kinase (MAPK), G-protein-coupled receptor (GPCR), protein kinase B (AKT), interferon, and Wingless and Int-related protein (WNT) pathways) and epigenetic modifications (histone crotonylation and methylation). This scalable, sequencing-based procedure elucidates the genetic landscapes that control protein states, identifying genes that cause very narrow phenotypic effects and genes that lead to broad phenotypic consequences. The resulting genetic wiring map identifies the E3-ligase substrate adaptor KCTD5 (ref. 1) as a negative regulator of the AKT pathway, a key signalling cascade frequently deregulated in cancer. KCTD5-deficient cells show elevated levels of phospho-AKT at S473 that could not be attributed to effects on canonical pathway components. To reveal the genetic requirements for this phenotype, we iteratively analysed the regulatory network linked to AKT activity in the knockout background. This genetic modifier screen exposes suppressors of the KCTD5 phenotype and mechanistically demonstrates that KCTD5 acts as an off-switch for GPCR signalling by triggering proteolysis of Gßγ heterodimers dissociated from the Gα subunit. Although biological networks have previously been constructed on the basis of gene expression, protein-protein associations, or genetic interaction profiles, we foresee that the approach described here will enable the generation of a comprehensive genetic wiring map for human cells on the basis of quantitative protein states.

The α-(1,3)-glucan synthase gene agsE impacts the secretome of Aspergillus niger.

Aspergillus niger is widely used as a cell factory for the industrial production of enzymes. Previously, it was shown that deletion of α-1-3 glucan synthase genes results in smaller micro-colonies in liquid cultures of Aspergillus nidulans. Also, it has been shown that small wild-type Aspergillus niger micro-colonies secrete more protein than large mirco-colonies. We here assessed whether deletion of the agsC or agsE α-1-3 glucan synthase genes results in smaller A. niger micro-colonies and whether this is accompanied by a change in protein secretion. Biomass formation was not affected in the deletion strains but pH of the culture medium had changed from 5.2 in the case of the wild-type to 4.6 and 6.4 for ΔagsC and ΔagsE, respectively. The diameter of the ΔagsC micro-colonies was not affected in liquid cultures. In contrast, diameter of the ΔagsE micro-colonies was reduced from 3304 ± 338 µm to 1229 ± 113 µm. Moreover, the ΔagsE secretome was affected with 54 and 36 unique proteins with a predicted signal peptide in the culture medium of MA234.1 and the ΔagsE, respectively. Results show that these strains have complementary cellulase activity and thus may have complementary activity on plant biomass degradation. Together, α-1-3 glucan synthesis (in)directly impacts protein secretion in A. niger.

Quantifying Positional Isomers (QPI) by Top-Down Mass Spectrometry.

Proteomics has exposed a plethora of posttranslational modifications, but demonstrating functional relevance requires new approaches. Top-down proteomics of intact proteins has the potential to fully characterize protein modifications in terms of amount, site(s), and the order in which they are deposited on the protein; information that so far has been elusive to extract by shotgun proteomics. Data acquisition and analysis of intact multimodified proteins have however been a major challenge, in particular for positional isomers that carry the same number of modifications at different sites. Solutions were previously proposed to extract this information from fragmentation spectra, but these have so far mainly been limited to peptides and have entailed a large degree of manual interpretation. Here, we apply high-resolution Orbitrap fusion top-down analyses in combination with bioinformatics approaches to attempt to characterize multiple modified proteins and quantify positional isomers. Automated covalent fragment ion type definition, detection of mass precision and accuracy, and extensive use of replicate spectra increase sequence coverage and drive down false fragment assignments from 10% to 1.5%. Such improved performance in fragment assignment is key to localize and quantify modifications from fragment spectra. The method is tested by investigating positional isomers of Ubiquitin mixed in known concentrations, which results in quantification of high ratios at very low standard errors of the mean (<5%), as well as with synthetic phosphorylated peptides. Application to multiphosphorylated Bora provides an estimation of the so far unknown stoichiometry of the known set of phosphosites and uncovers new sites from hyperphosphorylated Bora.

A System-wide Approach to Monitor Responses to Synergistic BRAF and EGFR Inhibition in Colorectal Cancer Cells.

Intrinsic and/or acquired resistance represents one of the great challenges in targeted cancer therapy. A deeper understanding of the molecular biology of cancer has resulted in more efficient strategies, where one or multiple drugs are adopted in novel therapies to tackle resistance. This beneficial effect of using combination treatments has also been observed in colorectal cancer patients harboring the BRAF(V600E) mutation, whereby dual inhibition of BRAF(V600E) and EGFR increases antitumor activity. Notwithstanding this success, it is not clear whether this combination treatment is the only or most effective treatment to block intrinsic resistance to BRAF inhibitors. Here, we investigate molecular responses upon single and multi-target treatments, over time, using BRAF(V600E) mutant colorectal cancer cells as a model system. Through integration of transcriptomic, proteomic and phosphoproteomics data we obtain a comprehensive overview, revealing both known and novel responses. We primarily observe widespread up-regulation of receptor tyrosine kinases and metabolic pathways upon BRAF inhibition. These findings point to mechanisms by which the drug-treated cells switch energy sources and enter a quiescent-like state as a defensive response, while additionally compensating for the MAPK pathway inhibition.

Next-generation proteomics: towards an integrative view of proteome dynamics.

Next-generation sequencing allows the analysis of genomes, including those representing disease states. However, the causes of most disorders are multifactorial, and systems-level approaches, including the analysis of proteomes, are required for a more comprehensive understanding. The proteome is extremely multifaceted owing to splicing and protein modifications, and this is further amplified by the interconnectivity of proteins into complexes and signalling networks that are highly divergent in time and space. Proteome analysis heavily relies on mass spectrometry (MS). MS-based proteomics is starting to mature and to deliver through a combination of developments in instrumentation, sample preparation and computational analysis. Here we describe this emerging next generation of proteomics and highlight recent applications.

A central role for TFIID in the pluripotent transcription circuitry.

Embryonic stem (ES) cells are pluripotent and characterized by open chromatin and high transcription levels, achieved through auto-regulatory and feed-forward transcription factor loops. ES-cell identity is maintained by a core of factors including Oct4 (also known as Pou5f1), Sox2, Klf4, c-Myc (OSKM) and Nanog, and forced expression of the OSKM factors can reprogram somatic cells into induced pluripotent stem cells (iPSCs) resembling ES cells. These gene-specific factors for RNA-polymerase-II-mediated transcription recruit transcriptional cofactors and chromatin regulators that control access to and activity of the basal transcription machinery on gene promoters. How the basal transcription machinery is involved in setting and maintaining the pluripotent state is unclear. Here we show that knockdown of the transcription factor IID (TFIID) complex affects the pluripotent circuitry in mouse ES cells and inhibits reprogramming of fibroblasts. TFIID subunits and the OSKM factors form a feed-forward loop to induce and maintain a stable transcription state. Notably, transient expression of TFIID subunits greatly enhanced reprogramming. These results show that TFIID is critical for transcription-factor-mediated reprogramming. We anticipate that, by creating plasticity in gene expression programs, transcription complexes such as TFIID assist reprogramming into different cellular states.

Robust, Sensitive, and Automated Phosphopeptide Enrichment Optimized for Low Sample Amounts Applied to Primary Hippocampal Neurons.

Because of the low stoichiometry of protein phosphorylation, targeted enrichment prior to LC-MS/MS analysis is still essential. The trend in phosphoproteome analysis is shifting toward an increasing number of biological replicates per experiment, ideally starting from very low sample amounts, placing new demands on enrichment protocols to make them less labor-intensive, more sensitive, and less prone to variability. Here we assessed an automated enrichment protocol using Fe(III)-IMAC cartridges on an AssayMAP Bravo platform to meet these demands. The automated Fe(III)-IMAC-based enrichment workflow proved to be more effective when compared to a TiO2-based enrichment using the same platform and a manual Ti(IV)-IMAC-based enrichment workflow. As initial samples, a dilution series of both human HeLa cell and primary rat hippocampal neuron lysates was used, going down to 0.1 µg of peptide starting material. The optimized workflow proved to be efficient, sensitive, and reproducible, identifying, localizing, and quantifying thousands of phosphosites from just micrograms of starting material. To further test the automated workflow in genuine biological applications, we monitored EGF-induced signaling in hippocampal neurons, starting with only 200 000 primary cells, resulting in ∼50 µg of protein material. This revealed a comprehensive phosphoproteome, showing regulation of multiple members of the MAPK pathway and reduced phosphorylation status of two glutamate receptors involved in synaptic plasticity.

Opposite Electron-Transfer Dissociation and Higher-Energy Collisional Dissociation Fragmentation Characteristics of Proteolytic K/R(X)n and (X)nK/R Peptides Provide Benefits for Peptide Sequencing in Proteomics and Phosphoproteomics.

A key step in shotgun proteomics is the digestion of proteins into peptides amenable for mass spectrometry. Tryptic peptides can be readily sequenced and identified by collision-induced dissociation (CID) or higher-energy collisional dissociation (HCD) because the fragmentation rules are well-understood. Here, we investigate LysargiNase, a perfect trypsin mirror protease, because it cleaves equally specific at arginine and lysine residues, albeit at the N-terminal end. LysargiNase peptides are therefore practically tryptic-like in length and sequence except that following ESI, the two protons are now both positioned at the N-terminus. Here, we compare side-by-side the chromatographic separation properties, gas-phase fragmentation characteristics, and (phospho)proteome sequence coverage of tryptic (i.e., (X)nK/R) and LysargiNase (i.e., K/R(X)n) peptides using primarily electron-transfer dissociation (ETD) and, for comparison, HCD. We find that tryptic and LysargiNase peptides fragment nearly as mirror images. For LysargiNase predominantly N-terminal peptide ions (c-ions (ETD) and b-ions (HCD)) are formed, whereas for trypsin, C-terminal fragment ions dominate (z-ions (ETD) and y-ions (HCD)) in a homologous mixture of complementary ions. Especially during ETD, LysargiNase peptides fragment into low-complexity but information-rich sequence ladders. Trypsin and LysargiNase chart distinct parts of the proteome, and therefore, the combined use of these enzymes will benefit a more in-depth and reliable analysis of (phospho)proteomes.

Myc coordinates transcription and translation to enhance transformation and suppress invasiveness.

c-Myc is one of the major human proto-oncogenes and is often associated with tumor aggression and poor clinical outcome. Paradoxically, Myc was also reported as a suppressor of cell motility, invasiveness, and metastasis. Among the direct targets of Myc are many components of the protein synthesis machinery whose induction results in an overall increase in protein synthesis that empowers tumor cell growth. At present, it is largely unknown whether beyond the global enhancement of protein synthesis, Myc activation results in translation modulation of specific genes. Here, we measured Myc-induced global changes in gene expression at the transcription, translation, and protein levels and uncovered extensive transcript-specific regulation of protein translation. Particularly, we detected a broad coordination between regulation of transcription and translation upon modulation of Myc activity and showed the connection of these responses to mTOR signaling to enhance oncogenic transformation and to the TGFß pathway to modulate cell migration and invasiveness. Our results elucidate novel facets of Myc-induced cellular responses and provide a more comprehensive view of the consequences of its activation in cancer cells.

Identification of putative substrates for the periplasmic chaperone YfgM in Escherichia coli using quantitative proteomics.

How proteins are trafficked, folded, and assembled into functional units in the cell envelope of Gram-negative bacteria is of significant interest. A number of chaperones have been identified, however, the molecular roles of these chaperones are often enigmatic because it has been challenging to assign substrates. Recently we discovered a novel periplasmic chaperone, called YfgM, which associates with PpiD and the SecYEG translocon and operates in a network that contains Skp and SurA. The aim of the study presented here was to identify putative substrates of YfgM. We reasoned that substrates would be incorrectly folded or trafficked when YfgM was absent from the cell, and thus more prone to proteolysis (the loss-of-function rationale). We therefore used a comparative proteomic approach to identify cell envelope proteins that were lower in abundance in a strain lacking yfgM, and strains lacking yfgM together with either skp or surA. Sixteen putative substrates were identified. The list contained nine inner membrane proteins (CusS, EvgS, MalF, OsmC, TdcB, TdcC, WrbA, YfhB, and YtfH) and seven periplasmic proteins (HdeA, HdeB, AnsB, Ggt, MalE, YcgK, and YnjE), but it did not include any lipoproteins or outer membrane proteins. Significantly, AnsB (an asparaginase) and HdeB (a protein involved in the acid stress response), were lower in abundance in all three strains lacking yfgM. For both genes, we ruled out the possibility that they were transcriptionally down-regulated, so it is highly likely that the corresponding proteins are misfolded/mistargeted and turned-over in the absence of YfgM. For HdeB we validated this conclusion in a pulse-chase experiment. The identification of HdeB and other cell envelope proteins as potential substrates will be a valuable resource for follow-up experiments that aim to delineate molecular the function of YfgM.

Proteomic analyses uncover a new function and mode of action for mouse homolog of Diaphanous 2 (mDia2).

mDia2 is an auto-inhibited Formin influencing actin dynamics upon conversion to the active conformation. mDia2 regulates actin-based protrusions and cell invasion, cell differentiation, vesicle trafficking, and cytokinesis. However, whether mDia2 has additional functions and how its action is functionally specified remain unknown. Here we draw the interactome of auto-inhibited and constitutively active mDia2 to address these issues. We embed mDia2 in protein networks accounting for its attributed functions and unexpectedly link it to the Ubiquitin Proteasome System. Taking FBXO3 as a test case, we show that mDia2 binds FBXO3 and p53, and regulates p53 transcriptional activity in an actin-nucleation-independent and conformation-insensitive manner. Increased mDia2 and FBXO3 levels elevate p53 activity and expression thereby sensitizing cells to p53-dependent apoptosis, whereas their decrease produces opposite effects. Thus, we discover a new role of mDia2 in p53 regulation suggesting that the closed conformation is biologically active and an FBXO3-based mechanism to functionally specify mDia2's activity.

Assessment of SRM, MRM(3) , and DIA for the targeted analysis of phosphorylation dynamics in non-small cell lung cancer.

Hypothesis-driven MS-based targeted proteomics has gained great popularity in a relatively short timespan. Next to the widely established selected reaction monitoring (SRM) workflow, data-independent acquisition (DIA), also referred to as sequential window acquisition of all theoretical spectra (SWATH) was introduced as a high-throughput targeted proteomics method. DIA facilitates increased proteome coverage, however, does not yet reach the sensitivity obtained with SRM. Therefore, a well-informed method selection is crucial for designing a successful targeted proteomics experiment. This is especially the case when targeting less conventional peptides such as those that contain PTMs, as these peptides do not always adhere to the optimal fragmentation considerations for targeted assays. Here, we provide insight into the performance of DIA, SRM, and MRM cubed (MRM(3) ) in the analysis of phosphorylation dynamics throughout the phosphoinositide 3-kinase mechanistic target of rapamycin (PI3K-mTOR) and mitogen-activated protein kinase (MAPK) signaling network. We observe indeed that DIA is less sensitive when compared to SRM, however demonstrates increased flexibility, by postanalysis selection of alternative phosphopeptide precursors. Additionally, we demonstrate the added benefit of MRM(3) , allowing the quantification of two poorly accessible phosphosites. In total, targeted proteomics enabled the quantification of 42 PI3K-mTOR and MAPK phosphosites, gaining a so far unachieved in-depth view mTOR signaling events linked to tyrosine kinase inhibitor resistance in non-small cell lung cancer.

Quantitative Proteomics Illuminates a Functional Interaction between mDia2 and the Proteasome.

Formin mDia2 is a cytoskeleton-regulatory protein that switches reversibly between a closed, autoinhibited and an open, active conformation. Although the open conformation of mDia2 induces actin assembly thereby controlling many cellular processes, mDia2 possesses also actin-independent and conformation-insensitive scaffolding roles related to microtubules and p53, respectively. Thus, we hypothesize that mDia2 may have other unappreciated functions and regulatory modes. Here we identify and validate proteasome and Ubiquitin as mDia2-interacting partners using stable isotope labeling with amino acids in cell culture-based quantitative proteomics and biochemistry, respectively. Although mDia2 is ubiquitinated, binds ubiquitinated proteins and free Ubiquitin, it is not a proteasome substrate. Surprisingly, knockdown of mDia2 increases the activity of the proteasome in vitro, whereas mDia2 overexpression has opposite effects only when it adopts the open conformation and cannot induce actin assembly. Consistently, a combination of candidate and unbiased proteome-wide analyses indicates that mDia2 regulates the cellular levels of proteasome substrate ß-catenin and a number of ubiquitinated actin-regulatory proteins. Hence, these findings add more complexity to the mDia2 activity cycle by showing that the open conformation may control actin dynamics also through actin-independent regulation of the proteasome.

Flexibility in crosstalk between H2B ubiquitination and H3 methylation in vivo.

Histone H2B ubiquitination is a dynamic modification that promotes methylation of histone H3K79 and H3K4. This crosstalk is important for the DNA damage response and has been implicated in cancer. Here, we show that in engineered yeast strains, ubiquitins tethered to every nucleosome promote H3K79 and H3K4 methylation from a proximal as well as a more distal site, but only if in a correct orientation. This plasticity indicates that the exact location of the attachment site, the native ubiquitin-lysine linkage and ubiquitination cycles are not critical for trans-histone crosstalk in vivo. The flexibility in crosstalk also indicates that other ubiquitination events may promote H3 methylation.

Single-step enrichment by Ti4+-IMAC and label-free quantitation enables in-depth monitoring of phosphorylation dynamics with high reproducibility and temporal resolution.

Quantitative phosphoproteomics workflows traditionally involve additional sample labeling and fractionation steps for accurate and in-depth analysis. Here we report a high-throughput, straightforward, and comprehensive label-free phosphoproteomics approach using the highly selective, reproducible, and sensitive Ti(4+)-IMAC phosphopeptide enrichment method. We demonstrate the applicability of this approach by monitoring the phosphoproteome dynamics of Jurkat T cells stimulated by prostaglandin E2 (PGE2) over six different time points, measuring in total 108 snapshots of the phosphoproteome. In total, we quantitatively monitored 12,799 unique phosphosites over all time points with very high quantitative reproducibility (average r > 0.9 over 100 measurements and a median cv < 0.2). PGE2 is known to increase cellular cAMP levels, thereby activating PKA. The in-depth analysis revealed temporal regulation of a wide variety of phosphosites associated not only with PKA, but also with a variety of other classes of kinases. Following PGE2 stimulation, several pathways became only transiently activated, revealing that in-depth dynamic profiling requires techniques with high temporal resolution. Moreover, the large publicly available dataset provides a valuable resource for downstream PGE2 signaling dynamics in T cells, and cAMP-mediated signaling in particular. More generally, our method enables in-depth, quantitative, high-throughput phosphoproteome screening on any system, requiring very little sample, sample preparation, and analysis time.

Phosphoproteome dynamics in onset and maintenance of oncogene-induced senescence.

Expression of the BRAF(V600E) oncoprotein is known to cause benign lesions, such as melanocytic nevi (moles). Despite the oncogenic function of mutant BRAF, these lesions are arrested by a cell-autonomous mechanism called oncogene-induced senescence. Infrequently, nevi can progress to malignant melanoma, through mechanisms that are incompletely understood. To gain more insight into this vital tumor-suppression mechanism, we performed a mass-spectrometry-based screening of the proteome and phosphoproteome in cycling and senescent cells and in cells with abrogated senescence. Proteome analysis of senescent cells revealed the up-regulation of established senescence biomarkers, including specific cytokines, but also several proteins not previously associated with senescence, including extracellular matrix-interacting. Using both general and targeted phosphopeptide enrichment by Ti(4+)-IMAC and phosphotyrosine antibody enrichment, we identified over 15,000 phosphorylation sites. Among the regulated phosphorylation sites we encountered components of the interleukin, BRAF/MAPK, and CDK-retinoblastoma pathways and several other factors. The extensive proteome and phosphoproteome dataset of BRAF(V600E)-expressing senescent cells provides molecular clues as to how oncogene-induced senescence is initiated, maintained, or evaded, serving as a comprehensive proteomic basis for functional validation.

PhosphoPath: Visualization of Phosphosite-centric Dynamics in Temporal Molecular Networks.

Protein phosphorylation is an essential post-translational modification (PTM) regulating many biological processes at the cellular and multicellular level. Continuous improvements in phosphoproteomics technology allow the analysis of this PTM in an expanding biological content, yet up until now proteome data visualization tools are still very gene centric, hampering the ability to comprehensively map and study PTM dynamics. Here we present PhosphoPath, a Cytoscape app designed for the visualization and analysis of quantitative proteome and phosphoproteome data sets. PhosphoPath brings knowledge into the biological network by importing publically available data and enables PTM site-specific visualization of information from quantitative time series. To showcase PhosphoPath performance we use a quantitative proteomics data set comparing patient-derived melanoma cell lines grown in either conventional cell culture or xenografts.

Signal Transduction Reaction Monitoring Deciphers Site-Specific PI3K-mTOR/MAPK Pathway Dynamics in Oncogene-Induced Senescence.

We report a straightforward strategy to comprehensively monitor signal transduction pathway dynamics in mammalian systems. Combining targeted quantitative proteomics with highly selective phosphopeptide enrichment, we monitor, with great sensitivity, phosphorylation dynamics of the PI3K-mTOR and MAPK signaling networks. Our approach consists of a single enrichment step followed by a single targeted proteomics experiment, circumventing the need for labeling and immune purification while enabling analysis of selected phosphorylation nodes throughout signaling pathways. The need for such a comprehensive pathway analysis is illustrated by highlighting previously uncharacterized phosphorylation changes in oncogene-induced senescence, associated with diverse biological phenotypes and pharmacological intervention of the PI3K-mTOR pathway.

Diet-induced neuropeptide expression: feasibility of quantifying extended and highly charged endogenous peptide sequences by selected reaction monitoring.

Understanding regulation and action of endogenous peptides, especially neuropeptides, which serve as inter- and intracellular signal transmitters, is key in understanding a variety of functional processes, such as energy balance, memory, circadian rhythm, drug addiction, etc. Therefore, accurate and reproducible quantification of these bioactive endogenous compounds is highly relevant. The biosynthesis of endogenous peptides, involving multiple possible trimming and modification events, hinders the de novo prediction of the active peptide sequences, making MS-based measurements very valuable in determining the actual active compounds. Here, we report an extended selected reaction monitoring (SRM)-based strategy to reproducibly and quantitatively monitor the abundances of a set of 15 endogenously occurring peptides from Rattus norvegicus hypothalamus. We demonstrate that SRM can be extended toward reproducible detection and quantification of peptides, bearing characteristics very different from tryptic peptides. We show that long peptide sequences, producing precursors with up to five and MS2 fragment ions with up to three charges, can be targeted by SRM on a triple quadrupole instrument. Using this approach to quantify endogenous peptide levels in hypothalami of animals subjected to different diets revealed several significant changes, most notably the significant upregulation of VGF-derived signaling peptide AQEE-30 upon high caloric feeding.

Benchmarking multiple fragmentation methods on an orbitrap fusion for top-down phospho-proteoform characterization.

Top-down analysis of intact proteins by mass spectrometry provides an ideal platform for comprehensive proteoform characterization, in particular, for the identification and localization of post-translational modifications (PTM) co-occurring on a protein. One of the main bottlenecks in top-down proteomics is insufficient protein sequence coverage caused by incomplete protein fragmentation. Based on previous work on peptides, increasing sequence coverage and PTM localization by combining sequential ETD and HCD fragmentation in a single fragmentation event, we hypothesized that protein sequence coverage and phospho-proteoform characterization could be equally improved by this new dual fragmentation method termed EThcD, recently been made available on the Orbitrap Fusion. Here, we systematically benchmark the performance of several (hybrid) fragmentation methods for intact protein analysis on an Orbitrap Fusion, using as a model system a 17.5 kDa N-terminal fragment of the mitotic regulator Bora. During cell division Bora becomes multiply phosphorylated by a variety of cell cycle kinases, including Aurora A and Plk1, albeit at distinctive sites. Here, we monitor the phosphorylation of Bora by Aurora A and Plk1, analyzing the generated distinctive phospho-proteoforms by top-down fragmentation. We show that EThcD and ETciD on a Fusion are feasible and capable of providing richer fragmentation spectra compared to HCD or ETD alone, increasing protein sequence coverage, and thereby facilitating phosphosite localization and the determination of kinase specific phosphorylation sites in these phospho-proteoforms. Data are available via ProteomeXchange with identifier PXD001845.