Next Article in Journal
Computational and Experimental Evidence for Templated Macrocyclization: The Role of a Hydrogen Bond Network in the Quantitative Dimerization of 24-Atom Macrocycles
Next Article in Special Issue
Prediction of Phage Virion Proteins Using Machine Learning Methods
Previous Article in Journal
A Review of the Regulatory Mechanisms of N-Myc on Cell Cycle
Previous Article in Special Issue
Proteomic Analysis of the Effect of Salmonella Challenge on Broiler Chicken
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Targeted Quantification of Protein Phosphorylation and Its Contributions towards Mathematical Modeling of Signaling Pathways

Biological Sciences Division and Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory, Richland, WA 99352, USA
*
Author to whom correspondence should be addressed.
Molecules 2023, 28(3), 1143; https://doi.org/10.3390/molecules28031143
Submission received: 18 November 2022 / Revised: 12 January 2023 / Accepted: 18 January 2023 / Published: 23 January 2023
(This article belongs to the Special Issue Functional Proteomics in Cell Biology and Beyond)

Abstract

:
Post-translational modifications (PTMs) are key regulatory mechanisms that can control protein function. Of these, phosphorylation is the most common and widely studied. Because of its importance in regulating cell signaling, precise and accurate measurements of protein phosphorylation across wide dynamic ranges are crucial to understanding how signaling pathways function. Although immunological assays are commonly used to detect phosphoproteins, their lack of sensitivity, specificity, and selectivity often make them unreliable for quantitative measurements of complex biological samples. Recent advances in Mass Spectrometry (MS)-based targeted proteomics have made it a more useful approach than immunoassays for studying the dynamics of protein phosphorylation. Selected reaction monitoring (SRM)—also known as multiple reaction monitoring (MRM)—and parallel reaction monitoring (PRM) can quantify relative and absolute abundances of protein phosphorylation in multiplexed fashions targeting specific pathways. In addition, the refinement of these tools by enrichment and fractionation strategies has improved measurement of phosphorylation of low-abundance proteins. The quantitative data generated are particularly useful for building and parameterizing mathematical models of complex phospho-signaling pathways. Potentially, these models can provide a framework for linking analytical measurements of clinical samples to better diagnosis and treatment of disease.

1. Introduction

Protein phosphorylation is the most common and widely studied post-translational modification (PTM) in eukaryotes. Phosphorylation is involved in numerous signal transduction processes including proliferation, differentiation, apoptosis, and metabolic processes that maintain cellular homeostasis and direct physiological responses to stimuli [1]. Phosphorylation regulates function by multiple mechanisms, including conformational changes that alter protein structure and creating docking sites for the assembly of multi-protein complexes [2]. Kinases introduce negatively charged phosphate groups which interact with neutral hydroxyl groups on Serine (Ser), Threonine (Thr), and Tyrosine (Tyr) residues, resulting in a polar anionic moiety. Conversely, phosphatases reverse this biochemical process with the removal of the phosphate group. The dynamics of phosphorylation by these two counteracting enzymes are tightly controlled to properly propagate signal transduction events and modulate protein activity in each biochemical cascade [1,2,3,4,5]. An aberration of protein phosphorylation within any network can cause adverse effects that manifest in many disease states including cancer [6,7,8], diabetes [9,10,11], Alzheimer’s disease [12,13], rheumatoid arthritis [14,15], muscular dystrophy [16] and cardiovascular diseases [17,18]. Furthermore, multiple treatment strategies target protein kinases such as: Imatinib, a tyrosine kinase inhibitor for chronic myelogenous leukemia [19]; Duvelisib, a PI3K inhibitor for chronic lymphocytic leukemia and follicular lymphoma [20,21]; Larotrectinib, a NTRK inhibitor for solid tumors [22]; and Vemurafenib, a BRAF serine/threonine kinase inhibitor for melanoma [23], demonstrating the importance of signal transduction in cancer [24]. In addition, measuring protein phosphorylation is becoming centrally important to understanding healthy states, disease pathogenesis, and identify treatment strategies. Hence, identifying and quantifying protein phosphorylation within the context of complex phosphorylation-mediated signaling cascades has garnered substantial interest over the last decade.
Although substantial progress has been made in elucidating signaling events and uncovering therapeutic strategies, many challenges remain in understanding the relationship between signaling dynamics and disease progression. Specifically, the functions of signaling networks are based on the spatiotemporal dynamics that emerge from extensive protein groups and interactions collectively, rather than the simple presence or absence of individual proteins or interactions [25,26,27]. Thus, understanding signaling dynamics requires the quantitative measurements of multiple sites on multiple proteins through time and space. Advancements in the field of phosphoproteomics have allowed identification and quantification of thousands of phosphorylation sites [28,29,30,31]. Targeted phosphoproteomics has enabled the identification of site-specificity, localization of phosphorylation, the extent of phosphorylation, and quantification of absolute abundance with high throughput. This information can help to not only elucidate the biological relevance of proteins but also provide the high-resolution data needed to parameterize predictive mathematical models [32].
Computational modeling is also becoming increasingly important in understanding the collective behaviors of signaling proteins which are difficult to comprehend intuitively. Many different modeling paradigms and inference algorithms have been proposed to facilitate modeling, simulation, and inference of different types of biological systems [33]. In particular, the kinetic modeling of reaction networks and its derivative modeling method, rule-based modeling, are designed for signaling networks [34,35]. Modeling and experimental approaches are complementary and need to be integrated so that phosphorylation data can inform model structures and parameters while models, in turn, can guide experimental designs. Therefore, it is necessary to consider the information required for different types of models when designing phosphoproteomics experiments. In this review, we discuss the challenges facing targeted phosphoproteomics, the analytical advancements made to address those challenges, and the prospects of using targeted phosphoproteomics data to build predictive models of signaling pathways for biomedical applications.

2. The Case for Measuring Protein Phosphorylation by Mass Spectrometry

Current research in protein phosphorylation is focused on identifying phosphorylation sites, quantifying their abundance, and characterizing functional impact. From an analytical standpoint, there are inherent challenges in phosphorylation analyses due to the labile nature of the modification and sub-stoichiometric concentrations [36]. While different biochemical approaches have been developed and utilized to study protein phosphorylation, they all have inherent limitations. A few common examples are phosphoimaging and immunoassays. Historically, phosphoimaging initially depended on the use of radioisotopes, such as the incorporation of the 32P-labeled ATP to monitor kinase activity towards specific substrates [37]. More recently, fluorescent-based dyes have been used to stain for phospho-residues, selectively or universally, by direct imaging on blotting membranes [38,39]. Alternatively, antibody-based immunoassays like immunoblotting and enzyme-linked immunosorbent assay (ELISA) have been used as a substitute to radiolabeling techniques and can allow for signal amplification [40]. Unfortunately, the above-mentioned methods are generally not site-specific and do not provide stoichiometry information of individual phosphosites [40]. Consequently, the emerging mass spectrometry (MS)-based technologies aim to address these limitations of sensitivity, specificity, and rigor in measuring protein phosphorylation in complex biological samples.

MS Analysis of Protein Phosphorylation

The development of Mass Spectrometry (MS)-based proteomics has revolutionized the study of protein phosphorylation. Conceptual and technical breakthroughs in MS-based phosphoproteomics have enabled the detection and quantification of tens of thousands of phosphosites per sample [30,41,42]. Generally, proteins are digested to produce peptides that are separated by reversed-phase liquid chromatography and analyzed by mass spectrometry (LC-MS). Shotgun proteomics allows the discovery of phosphorylation sites through a data-dependent acquisition (DDA) mode. The DDA mode relies on a full-scan MS spectrum determination of peptide species present, then MS/MS acquisition of a list of most intense peptide precursor ions [43]. DDA approaches are suitable for exploring complex biological samples and identifying large numbers of phosphopeptides without prior knowledge [44]. Therefore, DDA is responsible for the identification of hundreds of thousands of phosphorylation sites across multiple specimens and species [45]. Such site-specific phosphorylation information has been collated into publicly available resource databases for use by the research community. PhosphoSitePlus [46] is the most comprehensive database with over 290,000 non-redundant phosphosites, while PHOSIDA [47,48], Phospho.ELM [49,50], and PhosphoPEP [51,52] contain over 70,000, 40,000, and 10,000 phosphosites, respectively, and counting. However, DDA is inherently biased towards more abundant phosphopeptides, which means that a phosphopeptide of interest is often not detected. In contrast, targeted methods of phosphoproteomics provide a more reliable way to monitor a specified number of ‘target’ phosphoproteins within and across biological samples. This facilitates hypothesis-driven research while simultaneously enabling high specificity and accurate quantification.

3. General Workflow for Targeted Phosphoproteomics

The general workflow for targeted phosphoproteomics is outlined in Figure 1. Typically, it involves the proteolytic digestion of proteins and post-digestion manipulations, such as enrichment to enhance recovery of phosphopeptides followed by MS analysis using a reference internal standard.

3.1. Surrogate ‘Target’ Phosphopeptide Selection and Enzymatic Digestion

The careful selection of target peptides is a crucial step in the identification and quantification of target proteins. The primary criteria for selection of surrogate peptides are that they must uniquely represent the protein of interest and be experimentally detectable within a proteome. They are often referred to as ‘proteotypic peptides’ [55,56]. Proteotypic peptides should be quantifiable (or quantotypic) and correlate with levels of the parent protein [55,57]. However, the selection of quantitative surrogate phosphopeptides for targeted analysis can be challenging because it must include the phosphorylated sites of interest. The typical approach of protein digestion using proteases can complicate the situation. Ideally, the enzyme(s) of choice should have efficient access to the cleavage sites. However, when phosphopeptides contain phosphorylated sites adjacent to the cleavage residues, they can form salt bridges which inhibit the activity of the enzyme and results in missed cleavages [58]. Trypsin is the most commonly used enzyme, which cleaves specifically at Lys and Arg C-terminal residues (except with the presence of Pro at the position after the cut site). Hence, the options for surrogate phosphopeptides can be limited. There are a few options to select alternative peptides. First, utilizing other proteolytic enzymes such as LysC (cleaves Lys C-terminal residue), AspN (cleaves N-terminal aspartic and glutamic acid residues), ArgC (cleaves C-terminal of Arg residues, including sites next to Pro), and GluC (cleaves C-terminal aspartic and glutamic acid residues) provides different options for cleavage sites. Multiple enzyme digests can be performed in parallel or sequentially to diversify peptide pools and increase coverage of the phosphoproteome. For example, sequential digestion using Glu-C and Lys-C enzymes prior to trypsin digestion has been shown to increase phosphosite identification by 83% and 40%, respectively [59,60]. Optimization of digestion efficiency, including the enzyme-to-substrate concentration ratio, optimal buffer conditions, and digestion time will be helpful with achieving reproducibility and accurate quantification.

3.2. Enrichment

Ongoing analytical challenges for targeted analyses include high peptide lability, sub-stoichiometric concentrations, heterogenous iso-forms, and low mass spectrometric response to protein phosphorylation [36]. It is estimated that at least 75% of the proteome can be phosphorylated, with kinases and phosphatases accounting for 1–2% of the genome [42,61]. In addition, the approximate relative abundance for pSer/pThr/pTyr species is 1800:200:1 [28]. Thus, phospho-tyrosine residues, which are critical signatures for tyrosine receptor kinase activation, constitute a very small fraction of all generated peptides. Hence, the efficient enrichment of pSer, pThr, and pTyr residues is crucial for robust and reliable phosphoproteome analyses. Towards this end, chemical resin- and antibody-based enrichment strategies have been developed and applied to enrich either total or tyrosine-phosphorylated phosphopeptides prior to MS analysis.

3.2.1. Immobilized Metal Affinity Chromatography (IMAC)

The most common phosphoenrichment technique introduced over three decades ago is IMAC. It is based on the principle that positively charged metal ions such as Fe3+, Ga3+, Ti4+, and Zr4+ can chelate and capture negatively charged phosphopeptides [62,63,64,65]. These metal ions are immobilized on resins that are commonly coated with iminodiacetic acid (IDA) or nitrilotriacetic acid (NTA) [66]. However, IMAC workflows can show reduced sensitivity caused by buffer interference [67] and its specificity can be reduced by acidic peptides that contain multiple aspartic acids and glutamic acids [68]. Though methyl esterification of the acidic amino acid residues has been used to prevent non-specific binding, it introduced more biological complexity as those modified residues can be prone to spontaneous deamidation and subsequent methylation [69,70]. An alternative approach is the complete protonation of acidic amino acid residues with acid at a low pH (below 2) [62]. It is also important to remove all traces of nucleic acids during sample protein extraction. These can competitively bind to the IMAC matrix, thus reducing IMAC efficiency [71].
To increase specificity and improve robustness of phosphoenrichment, several modifications have been applied to IMAC methodologies. Elements like Nickel have been linked to NTA-based resins to improve specificity [72]. Metal (IV) phosphate/phosphonate chemistry has revealed stronger coordination of metals like Zr4+ and Ti4+ with phosphates and carboxylic acid groups, and this has been implemented in approaches like the Ti4+-IMAC technology [63,73]. The preferential affinity of different metals towards specific phosphopeptides has led to sequential enrichment strategies. For example, Ti4+-IMAC has been shown to be able to enrich multiple basic amino acid residues in addition to acidic residues [74]. Another example is the sequential use of Ga3+-Fe3+-IMAC where Ga3+ is initially applied to enrich multi-phosphorylated peptides followed by Fe3+ to capture monophosphorylated peptides [75]. The extensive optimalization of IMAC-based enrichment strategies has significantly improved their usability for targeted quantification.

3.2.2. Metal Oxide Affinity Chromatography (MOAC)

MOAC enrichment is based on the Lewis acid-base interaction principle where metal oxides and hydroxides interact with oxygen atoms of phosphate groups to form complexes. Compounds like TiO2, ZrO2, Fe3O4, Ga2O3, Nb2O3, SnO2, Ta2O5, HfO2, and Al(OH)3 have been tested as matrixes for enrichment of phosphopeptides [76]. The first application of MOAC trapped phosphopeptides under acidic conditions on a TiO2-based precolumn. These were subsequently eluted in alkaline buffers for LC-MS/MS analysis [77]. The buffer conditions used in MOAC significantly influence its selectivity for phosphopeptides. For example, the addition of 2.5-dihydroxybenzoic acid (DHB) to loading and washing buffers reduces nonspecific binding of non-phosphorylated peptides to TiO2 and thus, increased the binding affinity for phosphopeptides [70,78]. However, residual DHB interferes with LC-MS/MS analysis. Alternatively, aliphatic hydroxy acid modified metal oxide chromatography (HAMMOC) using lactic acid and β-hydroxypropanoic acid in place of DHB improved phosphopeptide selectivity [78]. Furthermore, sequential elution on a pH gradient has been used to increase enrichment efficiency. For example, trapped mono-phosphopeptides in a TiO2 microcolumn can be first eluted with acidic methylphosphonate buffer followed by successive basic elution of multi-phosphopeptides with disodium hydrogen phosphate, ammonium hydroxide, and pyrrolidine [79].
Since IMAC and MOAC possess different affinities for phosphopeptides, it can be beneficial to combine both enrichment strategies. IMAC and MOAC typically favor multi-phosphorylated peptides and mono-phosphorylated peptides, respectively [67,80]. One combined approach is sequential enrichment with immobilized metal affinity chromatography (SIMAC). Here, the sample is first enriched with Fe3+-IMAC and multi-phosphopeptides are eluted in alkaline conditions, followed by TiO2 enrichment of the unbound acidic fraction for mono-phosphopeptides [81,82]. Like SIMAC, sequential enrichment using metal oxide affinity chromatography (SMOAC) offers a reverse approach. The samples are first enriched by TiO2 resin followed by Fe3+-IMAC [83]. Combining IMAC and MOAC strengthens binding selectivity and reduces the amount of starting material while deriving two distinct phosphopeptide fractions for analysis.

3.2.3. Polymer-Based Metal Affinity Capture (PolyMAC)

PolyMAC is another adaptation of IMAC but using soluble nanopolymers instead of solely metal ions. The nanopolymers are dentrimers multifunctionalized with metal ions that bind to phosphopeptides in solution. The PolyMAC-phosphopeptide complexes are subsequently recovered on a solid phase by coupling to hydrazine-agarose gels, followed by washing and phosphopeptide elution. The PolyMAC-Ti formulation has been well studied and appears to display high recovery, reproducibility, and enrichment compared to IMAC and TiO2 methods [84,85,86].

3.2.4. Antibody-Based Enrichment

Many traditional biochemical separation approaches are based on antibodies binding to specific epitopes. This has been applied to phosphoproteomics using antibodies raised against phosphorylated serine, threonine, and tyrosine residues. The metal-based methods discussed above are efficient for pSer and pThr residues, but not pTyr due to their significantly low relative abundance [87]. Fortunately, highly specific phosphotyrosine antibodies have proven to be effective in capturing the tyrosine phosphorylated peptides and proteins [88,89,90,91]. In a first application, selective immunoprecipitation of tyrosine phosphorylated proteins using anti-phosphotyrosine antibodies was able to identify known tyrosine phosphorylated proteins as well as novel sites in the epidermal growth factor (EGF) receptor signaling pathway [92]. However, most anti-phosphoserine and threonine antibodies do not bind efficiently to their respective sites and hence are not widely used in phosphoproteomics. Moreover, the need for site-specific antibodies makes immunoprecipitation unsuitable for wide coverage studies of the phosphoproteome. Nevertheless, the specificity and selectivity of antibodies provides a useful strategy when targeting specific low abundant proteins, as in the case of tyrosine-phosphorylated proteins. Furthermore, antibody- and IMAC-based strategies can be combined in a single workflow for deeper coverage of the phosphoproteome [93].

3.3. Fractionation

Fractionation is a step employed to reduce sample complexity either prior to phosphoenrichment or MS analysis. These strategies separate complex peptide samples into fractions that are more suitable for further enrichment or to enhance coverage depth of the phosphoproteome [94].
Strong cation exchange (SCX) and strong anion exchange (SAX) are chromatography techniques that fractionate peptides according to their charge. In SCX, positively charged cations interact with the negatively charged SCX matrix. The pH of the loading and elution buffers determine the binding and elution of different peptides. Acid buffers (pH 2.7) are used to protonate the C-terminal arginine and lysine residues and N-terminal amino acid groups, separating phosphorylated peptides from non-phosphorylated peptides by net charge states [95]. Bound peptides are eluted by increasing the ionic strength and/or pH to maximize collection of phosphopeptides in early fractions [96]. In contrast, SAX chromatography depends on the interaction of negatively charged peptides with a positively charged matrix. It is carried out using neutral-to-alkaline buffers so that the C-terminal carboxylic group and glutamate and aspartate residues are deprotonated. Most phosphopeptides are collected in later fractions by increasing ionic strength or decreasing the pH of the eluting buffer [97]. The applications of SCX and SAX are complementary as SCX favors retainment of mono-phosphorylated peptides while SAX favors multi-phosphorylated peptides [95,98]. However, selectively is an issue as other peptides with the same net charge states cannot be distinguished from phosphopeptides [94].
High pH reversed-phase liquid chromatography (HpH RPLC) is another fractionation strategy based on hydrophobic interactions [99]. It is carried out in a high pH environment and orthogonal separation is achieved by low to high pH gradient in the mobile phase [100]. HpH RPLC offers a higher resolving power compared to SCX and SAX and can be coupled with low pH LC-MS/MS [101]. In offline mode RPLC also offers the option of pooling fractions from different parts of the gradient, thus reducing the number of instrument runs for analysis [101,102]. To advance this strategy our group developed PRISM (high-pressure, high-resolution separation with intelligent selection and multiplexing) [103]. Here, high pH reversed-phase capillary LC is used to fractionate samples. Then an online SRM system monitors heavy isotope-labeled internal standards to select target fractions for downstream LC-SRM. Target fractions eluted at different times can be concatenated to increase throughput. The PRISM workflow provides a more-sensitive targeted proteomics workflow that does not require affinity enrichment [104]. RPLC fractionation prior to phosphopeptide enrichment by IMAC or MOAC has yielded extensive coverage of the phosphoproteome [105,106,107,108].
Other fractionation strategies like hydrophilic interaction chromatography (HILIC) and electrostatic repulsion-hydrophilic interaction (ERLIC) are based on hydrophilic interaction [109,110]. HILIC relies on hydrogen bond interaction between a neutral, hydrophilic stationary phase and peptides. More hydrophilic peptides are retained longer and eluted by increasing the polarity of the mobile phase, which is opposite to RPLC [109]. Alternatively, ERLIC is a subset of HILIC which utilizes both hydrophilic interactions and electrostatic forces. In the presence of high-percentage organic mobile phase at low pH, hydrophilic interactions dominate over electrostatic forces, making phosphopeptides less repulsed (or more retained) by the column compared to non-phosphorylated peptides. However, the absence of hydrophilic interactions leads to electrostatic repulsion of phosphopeptides at low pH [110,111].

4. Targeted MS Analysis of Phosphorylation

Current targeted methods of phosphoproteomics are based on selected reaction monitoring (SRM, also called multiple reaction monitoring (MRM)) [55] to provide precise quantification, reliable reproducibility, and robust sensitivity of ‘target’ proteins/phosphoproteins in biological samples. As technological advancements in instrumentation continue, there are targeted acquisition methods like parallel reaction monitoring (PRM) and data independent acquisition (DIA)-MS strategies that are gaining relevance. Due to the sub-stoichiometric nature of phosphorylation events, the continuous gain in sensitivity of targeted MS acquisition tools is useful to quantitatively and consistently measure target phosphopeptides. The current methods will be discussed subsequently.

4.1. SRM/MRM

Selected Reaction Monitoring (SRM), also known as Multiple Reaction Monitoring (MRM), is an approach for reliably measuring phosphoproteins in a highly reproducible manner within and across laboratories and instrument platforms [112,113]. Implementation of SRM/MRM to quantify phosphorylation requires that the mass shifts introduced by phosphorylation are known. MS/MS spectra are isolated and detected by linear collision-induced dissociation (CID) on traditional triple quadrupole (QQQ) mass spectrometer platforms. The first quadrupole (Q1) filters and selects predefined precursors, the second quadrupole (Q2) induces fragmentation of the species via CID, and the third quadrupole (Q3) detects the resulting product/fragment ions with high speed and sensitivity [114]. Fragment ion analysis occurs in different scan modes. In the precursor ion scan mode, only phosphate (PO3)-modified fragments are detected in Q3 by the addition of 79 Da, which is only ideal for highly stable phosphopeptides. Alternatively, neutral loss scan mode detects the loss of a neutral group (H3PO4) from a phosphorylated peptide after low-energy CID. The loss of H3PO4 results in a mass shift of 98 Da, 49 Da and 32.7 Da for singly, doubly, and triply charged phosphopeptides, respectively. However, tyrosine-specific phosphopeptides are more stable and show no neutral loss, thus Q3 scans for an immonium ion (216.042 Da) marker [115]. The acquired mass-to-charge (m/z) values of precursor-fragment ion pairs monitored by Q1 and Q3 are called ‘transitions.’
In establishing SRM assays, selecting optimal transitions for each peptide is important. First, it is best practice to select the most intense fragments because the intensities of fragment ions derived from one precursor differ substantially. Low-energy CID induces the generation of b/y series fragment ions. In this fragmentation environment, y-ion signals are more common, robust, and stable compared to b-ion signals [116]. Even though y-ions are the preferred fragments for building transitions, b-ions can also be good candidates depending on the peptide. It is also good practice to exclude low m/z fragment ions due to low specificity and ensure that chosen fragments have m/z values above the precursor for high selectivity and spectral noise reduction [55]. In general, at least three transitions per phosphopeptide are usually selected for SRM analysis, to allow for room to include more targets in one assay. Once transitions are finalized, hundreds of phosphopeptides can be targeted in a single LC-MS/MS analysis. Targeted data analysis is typically performed using the Skyline (Figure 1), the most commonly used software for targeted proteomics data analysis, and other tools such as MRMer [117,118].
The first utilization of SRM for large-scale targeted phosphoproteomics was reported by White’s group in 2007 [119]. They were able to quantify phosphorylation changes in the Epidermal Growth Factor Receptor (EGFR) signaling network upon Epidermal Growth Factor (EGF) stimulation. Using a two-step enrichment procedure with selective antibody immunoprecipitation of phosphotyrosine-containing peptides followed by IMAC, they quantified 222 phosphorylated peptides across seven time points, in a single two-hour analysis.
For absolute quantification workflows, addition of chemically synthesized heavy isotope-labeled peptides containing the phosphorylated sites of interest is the gold standard for analyzing phosphopeptides and their stoichiometries [55,112,120]. These peptide analogs enable the identification of phosphopeptide targets with high confidence with reduced technical variation and false positives that could arise from co-eluting isoforms [44,53]. The retention time and ion fragmentation pattern are identical between heavy isotope-labeled and endogenous phosphopeptides, though distinguishable by a set mass difference. The intensity of the endogenous peptide relative to the added heavy isotope-labeled standard is the basis for quantification [55]. It is important to adjust the concentrations of spiked-in heavy peptides in order that the relative ratios of endogenous-to-heavy signals within a reasonable range for reliable quantification. To quantify differences between native and phosphorylated peptides, two internal standard peptides can be used: a non-phosphorylated and a phosphorylated peptide. The relative intensities of non-phosphorylated and phosphorylated peptides determine the percent change in phosphorylation at a given site which, in turn, can provide a measure of stoichiometry. Gerber et al. demonstrated this concept by quantifying the cell cycle-dependent Ser-1126 phosphorylation of Human Separase [121]. This approach can be integrated as a normalization strategy because changes in up to 25% of phosphopeptides can be attributable to varying protein levels [122]. With the complexities of biological mixtures, normalization strategies to correct for systematic biases are also often critical [123].
In general, standard-based quantification for targeted phosphoproteomics is gaining relevance with increasing experimental demonstrations by different groups (Table 1). The largest scale utilization of this approach so far has been by Carr’s group, quantifying 284 phosphosites in a multiplexed quantitative assay called SigPath [93]. Our group has developed assays to measure absolute abundances of proteins and phosphorylation sites over time in the Epidermal Growth Factor Receptor-Mitogen-Activated Protein Kinase (EGFR-MAPK) pathway to better understand rewiring in drug resistant cancers [32,53]. We have found that current AQUA tools can reliably measure changes in the phosphoproteome, and hold promise for more widespread research and modeling applications.

4.2. PRM

In the ongoing effort to avoid limiting the number of targets, technical advancements in high-resolution mass-measuring instruments like the hybrid quadrupole-orbitrap mass spectrometer has enabled the development of Parallel Reaction Monitoring (PRM) [130]. Like SRM, precursor ions of interest are selected and fragmented in Q1 and Q2, respectively. However, unlike SRM the Orbitrap replaces the Q3 to scan all fragment ions at high resolution, enabling the selection of the best transitions post-acquisition [131,132]. The Quadrupole-Orbitrap instrumentation setup also offers trapping capabilities which are valuable in analyzing low abundance events like phosphorylation. This approach reduces potential background interference usually encountered with complex biological samples, offering both qualitative capabilities and quantitative analysis [130].
PRM data acquisition takes advantage of certain Quadrupole-Orbitrap instrumentation parameters that can be optimized to improve analysis of complex samples. The quadrupole isolation width can be narrowed to directly control the selectivity of precursor ions and exclude interferences, and the trapping fill time can be longer to enhance sensitivity and widen dynamic detection range [130]. PRM has been used in a highly reproducible manner to measure 96 reduced-representation probes in ~200 samples across cell types, treatment conditions, and timepoints [126]. There are also a few small-scale studies that have used PRM to confirm and quantify protein targets from previous experimental observations [133,134,135]. A more recent data acquisition method, called internal standard triggered-parallel reaction monitoring (IS-PRM), uses internal standard peptides to set parameters and drive spectral acquisition. This real-time scheme dynamically monitors internal standard peptides, performs spectral matching, and then triggers quantitative analysis of heavy and endogenous peptides [136]. The versatility of PRM is a powerful approach for combining exploration and analysis in quantitative phosphoproteomics.

4.3. DIA/SWATH

Sequential window acquisition of all theoretical mass spectra (SWATH-MS) operating in a data-independent acquisition (DIA) mode is an approach designed to combine the advantages of both shotgun and targeted proteomics to achieve reproducibility and consistency in a high-throughput manner [137,138]. This method, like PRM, typically utilizes a Quadrupole-Orbitrap/TOF hybrid as a mass analyzer. In principle, the DIA step acquires the full m/z range of a selected precursor, followed by fragmentation in sequential small m/z windows. These acquisitions are continuously repeated to collect information on all detectable precursor and fragment ions within the specified windows [139]. Assay development with DIA/SWATH eliminates the requirement for identifying predetermined sets of ions as in SRM/PRM. Instead, defining mass ranges, precursor isolation window width, and number of MS2 scans per cycle are sufficient for data acquisition [137]. A recent study developed a hybrid DIA workflow with spectral libraries from DDA to quantify 36,350 phosphosites in lung cancer cell lines within 2 h [140]. Although this approach can quantify thousands of phosphopeptides, software development for data analysis is still undergoing improvement to extract phosphopeptide-specific information from available empirical knowledge. A few examples of currently available tools include OpenSWATH and DIA-NN [137,141,142].

5. Recent Targeted MS Applications to Elucidate Phosphorylation-Driven Mechanisms

Altered phosphorylation drives processes underlying numerous diseases. Disease-associated phosphorylation cascades can be restricted to specific pathways or, in many cases, involves crosstalk with other pathways. Hence, altered phosphorylation-mediated signaling can include tens to hundreds of proteins. These can potentially be followed and analyzed by SRM-based phosphoproteomics (Table 1). One prime example is the elucidation of phosphorylation-driven mechanisms in the EGFR-MAPK-PI3K-mTOR pathways that can be mapped to individual functional phosphosites (Figure 1B,C) [53,54]. For example, in the MAPK signal transduction cascade, MAPK1/ERK2 and MAPK3/ERK1 are two closely related and often indistinguishable protein isoforms, with corresponding phosphorylation sites at T185/Y187 and T202/Y204, respectively. The rapid activation of MAPK1/3 by dual phosphorylation shows similar trends in phosphorylation across these isoforms and indicates a level of redundancy between MAPK1 and MAPK 3. Ribosomal protein S6 (RPS6) is downstream of the MAPK pathway, integrating proliferating signals from both the MAPK and PI3K-mTOR pathways via S235 and S236 phosphorylation from ribosomal S6 and p70S6 kinases, respectively [54,143]. Consequently, phosphorylated RPS6 is a potential biomarker for the growth state of cancer cells and their drug resistance [144,145]. The ability to resolve RPS6 phosphorylation by two different kinases demonstrates the power of SRM assays.

6. Quantitative Phosphoproteomics in Modeling of Signaling Pathways

MS-based instrumentations have revealed extensive phosphorylation modifications in the proteome [29,30,69,95,146,147,148,149]. Though with limited data so far, the quantitative capabilities of SRM/MRM-based phosphoproteomics have demonstrated the ability to reliably measure phosphorylation abundances in signaling pathways (Table 1). These changes in phosphorylation in response to defined perturbations can provide useful data to different research areas from cellular biology to modelling. With the need to further understand how PTMs such as phosphorylation shape cellular function, mathematical models are a promising way to associate specific mechanisms with cellular perturbations [150].
The complexity of phosphorylation-mediated signaling and the subsequent heterogeneity make systems-level models most useful. To highlight some relevant sources of cellular complexities: (1) the presence of potential phosphorylation sites on receptors does not mean all sites can be phosphorylated at the same time, or ever at all [151]. Sometimes, spatially restricted subsets of receptors are not in the same compartment as effector proteins [152,153]; (2) phosphorylated sites can bind to adaptor proteins that protect them from dephosphorylation [154]; (3) there are feedback loops that can modulate pathway phosphorylation activity [155]; and (4) many protein–protein interactions are regulated by phosphorylation states, and this binding interaction is dependent on the concentration of binding partners and interaction affinity [156]. Multiple kinases are often coregulated, rather than acting independently, with patterns that modulate molecular switches which contribute toward the cellular response of a perturbation [29,150,157]. Considering the multifaceted nature of phosphorylation in building mathematical models, we need analytical methods that provide information on the state of individual molecules rather than the average state of receptors. Recent advancements in MS-based targeted phosphoproteomics provide the platform to obtain such relevant experimental information. We can also quantify phosphorylation on receptors and its impact on signaling molecules with high throughput. In addition, quantitative approaches can compartmentalize these changes in organellar fractions such as the mitochondria to establish localization indices for phosphorylation modifications [150].
Predictive modeling in cellular signaling aims to develop mathematical tools that are capable of computing expected outputs from given inputs under defined conditions [158]. With a model, one can either explore the possible outcomes with different structures and parameter regimes or optimize the model structures and parameters towards defined functional phenotypes [25]. However, to understand the dynamics of real biological systems at normal or disease states, one needs to determine and/or estimate the model structure and parameter space under different conditions. The parameter spaces under which real signaling networks operate are only a small subset of the essentially infinite parameter space in which a model could theoretically operate. To constrain model parameters to a relevant set, we can use both biophysical properties (e.g., diffusion and spatial scales) and experimental measurements. For example, for signaling networks the range of stimulation should be constrained to physiologically measured limits and relevant timescales [159]. Phosphoproteomics data can also be useful in determining which processes are operational in a cell [160] and how the activation of different proteins in a pathway are corelated either temporally or through a stimulation dose-range. These correlations can provide a useful constraint for the parameter space of a model [161].
A few studies have undertaken modeling approaches using different approaches of generating phosphorylation data (Table 2). Our group was one of the pioneers of using targeted phosphoproteomics datasets for mathematical modeling applications [104]. In a recent paper, quantitative targeted phosphoproteomics was used to determine the absolute abundance and phosphorylation states of proteins in the MAPK pathway, uncovering the negative feedback regulatory mechanisms of the cascade [53]. In turn, these data were suitable for parametrizing a model to predict the resistance of receptor driven ERK signaling to RAF and MEK inhibition [32].

Modeling Approaches to Capture Phosphorylation Signaling Complexities

In general, mathematical models are generally useful for: (1) aiding experimental designs to test both simple and complex hypotheses; (2) reconciling discrepancies in experimental data; and (3) consolidating knowledge about a system. Simple models like the partial least squares regression (PLSR) and multiple linear regression (MLR) are based on linear regressions [166]. The data are separated into sets of input and output, and a linear solution that correlates the input to the output is determined. However, these models are unable to capture nonlinear phenomena and coupled effects [167], which can be more pronounced in the disease state.
In complex biological systems, reducing the biochemistry (reactions) of the system to construct a network of ordinary differential equations (ODE) or partial differential equations (PDE) is a common approach [168]. Beginning with an equation for each species reaction written in terms of rate of change, empirical data then provide estimates of values of unknown parameters such as kinetic constants [169]. To simplify these models, certain assumptions are frequently made, for instance multiple phosphorylation sites on a peptide/protein having a similar impact [170]. Including all the phosphorylation states of a protein in a model can lead to a combinatorial explosion [171,172] and an exponential increase in reactions and subsequent model size. For example, a study that modeled signaling in the ERK and Akt network that included 7 receptor dimers and 28 proteins in response to two ligands generated a total of 499 differential equations, 201 reaction rates, and 28 non-zero initial conditions [2. Biological dynamics complicate parameter optimization of ODE/PDE models. One way around the complexity of ODE/PDE models is to simplify the system by “coarse-graining” which reduces the number of equations to just those describing molecular patterns that operate interdependently and are relevant to the system [158,173,174].
Rule-based models offer an approach which factors in the combinatorial complexities of cellular signaling mechanisms [174]. These models are based on the modularity of interactions, representing biomolecular interactions in terms of rules instead of equations [158,175]. These rules then are used to generate reactions and species. The rules are comprised of recognition patterns, reactant-product maps, and rate laws [174]. Consequently, the rules can specify the features of proteins (condition), such as the modification state of a phosphorylation site, which are required for or affected by a particular interaction (action) which can then be executed as simulations. As such, multiple species may qualify as reactants as defined by a rule, and multiple reactions may be generated as defined by the same characteristic law rate [176]. The common languages to write and simulate rule-based models are BioNetGen Language (BNGL) [34] and Kappa [35]. These algorithms run their simulations as: (1) indirect methods—where the state of the system is defined by traditional population variables such as concentrations; (2) direct methods—where the state of the system is defined by the collective states of individual sites or particles; and (3) hybrid particle/population method—where a subset of species can be treated as population variables rather than particles [34,177,178,179]. An example of an application of rule-based models was reported by Chylek et al. A rule-based modeling approach was developed which used temporal global phosphoproteomics data to characterize reshaping of the T-cell phosphoproteome in response to TCR/CD28 co-stimulation. The implicated signaling mechanisms were dependent on specific site phosphorylations and phosphoforms of receptors. The model detected over 100 pTyr sites with greater than two-fold changes in abundance and identified novel negative regulatory sites that were subsequently validated experimentally [165]. Rule-based modeling has also been combined within a Python-based programming language, PySB, which promotes model extensibility and reuse [180].
Rule-based and ODE/PDE models rely heavily, and in some cases entirely, on generated experimental knowledge of well-studied systems [158], which is an opportunity to leverage MS-based phosphoproteomics. MS-based technologies can generate datasets containing site and abundance information that can provide a basis, or “rules” to build biologically realistic models. The ongoing generation of high-quality proteomics and phosphoproteomics data across different experimental sources (such as different cell lines and tissue types) will be essential for the generation of cell-type specific models, and to understand the basis by which such specificity arises [181].
Previous studies have shown the utility of experimental data for constraining models and estimating model parameters. In considering the power of different experimental methods we argue that absolute quantification methods, such as targeted phosphoproteomics, provide much greater advantage in parameter estimation than relative quantification of phosphorylation levels. Absolute quantification can greatly reduce parameter space by directly providing values for protein abundance (concentration) and phosphorylation levels. Most kinetic models use “lumped” rate constants that combine both protein abundance and activity. Since protein abundance varies across cell types [181], knowing this parameter provides a way to infer real differences in protein activity. In addition, functional protein abundance is frequently dependent on protein phosphorylation state. For example, phosphorylation of SOS1 allows binding of 14-3-3 proteins that removes SOS1 from the pool that can participate in signaling [182]. Knowing the functional abundance of proteins through a combinate of targeted proteomics and phosphoproteomics reduces the uncertainty in model parameters and significantly improves model identifiability (Figure 2).

7. Future Perspectives

Targeted MS-based quantitative proteomics is becoming increasing popular across scientific and clinical research areas due to significant advancements in assay sensitivity, specificity, and reproducibility. These exceptional characteristics make them suitable for multiple applications such as biomarker discovery in basic and translational research. Classic SRM/MRM-based proteomics has been applied to verify protein biomarkers for different types of diseases such as melanoma, prostate cancer, lung cancer, Parkinson’s disease, cardiovascular diseases, and type 1 diabetes using different specimens like plasma, serum, urine, cerebrospinal fluid, and tissue [183]. Clinical applications could benefit immensely from the implementation of targeted phosphoproteomics methods to identify comprehensive site-specific phosphorylation that are of diagnostic and prognostic assessment value, such as mechanisms of cancer resistance mediated by the PI3K/AKT/mTOR signaling network. The constant improvement and application of SRM-based phosphoproteomics will advance the detection and quantification of biologically relevant phosphorylation sites, which could serve as biomarkers.
However, the broad utility of phosphoproteomics is not without existing challenges. Despite great improvements in proteomics workflows, sample throughput in MS-based proteomics is still relatively low compared to antibody-based techniques, such as Reverse Phase Protein Arrays [184]. Phosphorylation sites that can serve as biomarkers are typically of very low abundance and can be obscured by irrelevant peptides. Fortunately, fast, high-resolution, and highly sensitive measuring instruments like the orbitrap are being integrated into acquisition workflows as in the case of PRM applications. An added benefit with the orbitrap is the potential for simultaneous discovery and targeted phosphoproteomics as demonstrated by DIA/SWATH. Such approaches can significantly reduce instrumentation time, increase throughput, and enable measurements across many samples. DIA-based phosphoproteomics platforms could also expedite deep mining of phosphorylation events in a targeted fashion. Furthermore, the gravitation towards artificial intelligence capabilities like machine learning proposes the utility of phosphoproteomics data to achieve more precise and scalable models in big data analytics. The applications of artificial intelligence to phosphoproteomics data fall in two folds. On one hand, the unsupervised machine learning algorithms where the program is allowed to independently find patterns and trends. On the other hand, the supervised machine learning algorithms can be used to better identify the spectrum from raw MS data as well as interpolate and extrapolate the unmeasured time points and doses with a high fidelity [185]. In either choice, it is important to consider the quality of data that the algorithm learns, because noisy data will not generate meaningful models. Worth noting is that the recently developed physics-informed machine learning methods [186], such as neural differential equations [187], sparse regression [188], and symbolic regression [189], can be further developed and applied in proteomic-data-driven systems identifications and predictions. These applications require more accurate information from measurements, where biophysical rules to be discovered can only operate with absolute quantification of variables in the machine learning models. Besides phosphoproteomics, the algorithms can be designed to integrate other data, such as genomics, epigenomics, transcriptomics, metabolomics, etc. [190].
Another challenge is that due to the low abundance levels of phosphorylation sites, a large amount of input material is usually required for enrichment and subsequent analysis. This could pose a limitation for applications that generate only a small amount of sample/specimen. With the introduction of single-cell proteomics, there is the possibility for pushing the sensitivity of targeted phosphoproteomics to measure target analytes in small amounts of samples, such as those generated by flow cytometry. In terms of specificity, mass spectrometers cannot always differentiate all target phosphopeptides, especially if multiple phosphorylations result in similar mass or structural changes to the peptides. Consequently, utilizing recent tools like structure for lossless ion manipulation (SLIM) in resolution-type analysis, which can distinguish targets by size and shape in addition to mass, may help to resolve this issue [191]. Furthermore, bioinformatic tools are being developed that can identify functionally important phosphorylation sites within a context of a multi-phosphorylated protein. These are based on deep-learning programs like AlphaFold2 [192] that can potentially predict protein structural changes arising from phosphorylation events. Combining advanced MS and bioinformatic tools could benefit SRM applications for site specificity and help overcome ambiguous site-specific identifications using traditional MS approaches.

Concluding Remarks

As biomedical research gravitates towards a systems approach to answer important and increasingly complex questions, the integration of experimental and computational tools will continue to increase. To support this work, targeted phosphoproteomics has provided a platform to obtain high quality data with scalable throughput. It also provides important information for building predictive models that can assist in experimental designs and clinical decision making. Experimental and computational integration hold promise to better decipher phosphorylation mechanisms underlying disease pathogenesis and identify potential biomarker/therapeutic candidates for use in traditional and precision medicine.

Author Contributions

P.P.D., S.F., H.S.W., W.-J.Q. drafted the manuscript and prepared figures; T.S., J.M.J., and W.-J.Q. provided conceptual and technical expertise. P.P.D., S.F., T.S., J.M.J., H.S.W., and W.-J.Q. edited and revised manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

Portions of this work were supported by National Institutes of Health grants U01 CA227544, U54 DA049116, U01 DK124020, and P41 GM103493.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created or analyzed in this review. Data sharing is not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hunter, T. Signaling—2000 and beyond. Cell 2000, 100, 113–127. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Pawson, T.; Gish, G.D.; Nash, P. SH2 domains, interaction modules and cellular wiring. Trends Cell Biol. 2001, 11, 504–511. [Google Scholar] [CrossRef] [PubMed]
  3. Pawson, T.; Scott, J.D. Protein phosphorylation in signaling—50 years and counting. Trends Biochem. Sci. 2005, 30, 286–290. [Google Scholar] [CrossRef]
  4. Tarrant, M.K.; Cole, P.A. The chemical biology of protein phosphorylation. Annu. Rev. Biochem. 2009, 78, 797–825. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Cohen, P. Signal integration at the level of protein kinases, protein phosphatases and their substrates. Trends Biochem. Sci. 1992, 17, 408–413. [Google Scholar] [CrossRef]
  6. Ruprecht, B.; Lemeer, S. Proteomic analysis of phosphorylation in cancer. Expert Rev. Proteom. 2014, 11, 259–267. [Google Scholar] [CrossRef]
  7. Blume-Jensen, P.; Hunter, T. Oncogenic kinase signalling. Nature 2001, 411, 355–365. [Google Scholar] [CrossRef]
  8. Zanivan, S.; Meves, A.; Behrendt, K.; Schoof, E.M.; Neilson, L.J.; Cox, J.; Tang, H.R.; Kalna, G.; van Ree, J.H.; van Deursen, J.M.; et al. In vivo SILAC-based proteomics reveals phosphoproteome changes during mouse skin carcinogenesis. Cell Rep. 2013, 3, 552–566. [Google Scholar] [CrossRef] [Green Version]
  9. Meyerovitch, J.; Backer, J.M.; Kahn, C.R. Hepatic phosphotyrosine phosphatase activity and its alterations in diabetic rats. J. Clin. Investig. 1989, 84, 976–983. [Google Scholar] [CrossRef] [Green Version]
  10. Ortsater, H.; Grankvist, N.; Honkanen, R.E.; Sjoholm, A. Protein phosphatases in pancreatic islets. J. Endocrinol. 2014, 221, R121–R144. [Google Scholar] [CrossRef]
  11. Stanley, W.J.; Litwak, S.A.; Quah, H.S.; Tan, S.M.; Kay, T.W.; Tiganis, T.; de Haan, J.B.; Thomas, H.E.; Gurzov, E.N. Inactivation of Protein Tyrosine Phosphatases Enhances Interferon Signaling in Pancreatic Islets. Diabetes 2015, 64, 2489–2496. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Neddens, J.; Temmel, M.; Flunkert, S.; Kerschbaumer, B.; Hoeller, C.; Loeffler, T.; Niederkofler, V.; Daum, G.; Attems, J.; Hutter-Paier, B. Phosphorylation of different tau sites during progression of Alzheimer’s disease. Acta Neuropathol. Commun. 2018, 6, 52. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Drummond, E.; Pires, G.; MacMurray, C.; Askenazi, M.; Nayak, S.; Bourdon, M.; Safar, J.; Ueberheide, B.; Wisniewski, T. Phosphorylated tau interactome in the human Alzheimer’s disease brain. Brain 2020, 143, 2803–2817. [Google Scholar] [CrossRef]
  14. Nie, H.; Zheng, Y.; Li, R.; Guo, T.B.; He, D.; Fang, L.; Liu, X.; Xiao, L.; Chen, X.; Wan, B.; et al. Phosphorylation of FOXP3 controls regulatory T cell function and is inhibited by TNF-alpha in rheumatoid arthritis. Nat. Med. 2013, 19, 322–328. [Google Scholar] [CrossRef]
  15. Kuuliala, K.; Kuuliala, A.; Hamalainen, M.; Koivuniemi, R.; Kautiainen, H.; Moilanen, E.; Repo, H.; Leirisalo-Repo, M. Impaired Akt Phosphorylation in Monocytes of Patients with Rheumatoid Arthritis. Scand. J. Immunol. 2017, 85, 155–161. [Google Scholar] [CrossRef] [Green Version]
  16. Kaliman, P.; Llagostera, E. Myotonic dystrophy protein kinase (DMPK) and its role in the pathogenesis of myotonic dystrophy 1. Cell. Signal. 2008, 20, 1935–1941. [Google Scholar] [CrossRef]
  17. El-Armouche, A.; Pohlmann, L.; Schlossarek, S.; Starbatty, J.; Yeh, Y.H.; Nattel, S.; Dobrev, D.; Eschenhagen, T.; Carrier, L. Decreased phosphorylation levels of cardiac myosin-binding protein-C in human and experimental heart failure. J. Mol. Cell. Cardiol. 2007, 43, 223–229. [Google Scholar] [CrossRef]
  18. Walker, L.A.; Fullerton, D.A.; Buttrick, P.M. Contractile protein phosphorylation predicts human heart disease phenotypes. Am. J. Physiol. Heart Circ. Physiol. 2013, 304, H1644–H1650. [Google Scholar] [CrossRef] [Green Version]
  19. Sawyers, C.L.; Hochhaus, A.; Feldman, E.; Goldman, J.M.; Miller, C.B.; Ottmann, O.G.; Schiffer, C.A.; Talpaz, M.; Guilhot, F.; Deininger, M.W.; et al. Imatinib induces hematologic and cytogenetic responses in patients with chronic myelogenous leukemia in myeloid blast crisis: Results of a phase II study. Blood 2002, 99, 3530–3539. [Google Scholar] [CrossRef] [PubMed]
  20. O’Brien, S.; Patel, M.; Kahl, B.S.; Horwitz, S.M.; Foss, F.M.; Porcu, P.; Jones, J.; Burger, J.; Jain, N.; Allen, K.; et al. Duvelisib, an oral dual PI3K-delta, gamma inhibitor, shows clinical and pharmacodynamic activity in chronic lymphocytic leukemia and small lymphocytic lymphoma in a phase 1 study. Am. J. Hematol. 2018, 93, 1318–1326. [Google Scholar] [CrossRef]
  21. Flinn, I.W.; Miller, C.B.; Ardeshna, K.M.; Tetreault, S.; Assouline, S.E.; Mayer, J.; Merli, M.; Lunin, S.D.; Pettitt, A.R.; Nagy, Z.; et al. DYNAMO: A Phase II Study of Duvelisib (IPI-145) in Patients with Refractory Indolent Non-Hodgkin Lymphoma. J. Clin. Oncol. 2019, 37, 912–922. [Google Scholar] [CrossRef] [PubMed]
  22. Drilon, A.; Laetsch, T.W.; Kummar, S.; DuBois, S.G.; Lassen, U.N.; Demetri, G.D.; Nathenson, M.; Doebele, R.C.; Farago, A.F.; Pappo, A.S.; et al. Efficacy of Larotrectinib in TRK Fusion-Positive Cancers in Adults and Children. N. Engl. J. Med. 2018, 378, 731–739. [Google Scholar] [CrossRef]
  23. Sosman, J.A.; Kim, K.B.; Schuchter, L.; Gonzalez, R.; Pavlick, A.C.; Weber, J.S.; McArthur, G.A.; Hutson, T.E.; Moschos, S.J.; Flaherty, K.T.; et al. Survival in BRAF V600-mutant advanced melanoma treated with vemurafenib. N. Engl. J. Med. 2012, 366, 707–714. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. National Institute of Diabetes and Digestive and Kidney Diseases. Protein Kinase Inhibitors. In LiverTox: Clinical and Research Information on Drug-Induced Liver Injury; National Institute of Diabetes and Digestive and Kidney Diseases: Bethesda, MD, USA, 2012. [Google Scholar]
  25. Feng, S.; Ollivier, J.F.; Swain, P.S.; Soyer, O.S. BioJazz: In silico evolution of cellular networks with unbounded complexity using rule-based modeling. Nucleic Acids Res. 2015, 43, e123. [Google Scholar] [CrossRef] [Green Version]
  26. Bowsher, C.G.; Swain, P.S. Environmental sensing, information transfer, and cellular decision-making. Curr. Opin. Biotechnol. 2014, 28, 149–155. [Google Scholar] [CrossRef]
  27. Perkins, T.J.; Swain, P.S. Strategies for cellular decision-making. Mol. Syst. Biol. 2009, 5, 326. [Google Scholar] [CrossRef]
  28. Chan, C.Y.; Gritsenko, M.A.; Smith, R.D.; Qian, W.J. The current state of the art of quantitative phosphoproteomics and its applications to diabetes research. Expert Rev. Proteom. 2016, 13, 421–433. [Google Scholar] [CrossRef] [Green Version]
  29. Olsen, J.V.; Blagoev, B.; Gnad, F.; Macek, B.; Kumar, C.; Mortensen, P.; Mann, M. Global, in vivo, and site-specific phosphorylation dynamics in signaling networks. Cell 2006, 127, 635–648. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Olsen, J.V.; Vermeulen, M.; Santamaria, A.; Kumar, C.; Miller, M.L.; Jensen, L.J.; Gnad, F.; Cox, J.; Jensen, T.S.; Nigg, E.A.; et al. Quantitative phosphoproteomics reveals widespread full phosphorylation site occupancy during mitosis. Sci. Signal. 2010, 3, ra3. [Google Scholar] [CrossRef]
  31. Ishihama, Y.; Sato, T.; Tabata, T.; Miyamoto, N.; Sagane, K.; Nagasu, T.; Oda, Y. Quantitative mouse brain proteomics using culture-derived isotope tags as internal standards. Nat. Biotechnol. 2005, 23, 617–621. [Google Scholar] [CrossRef]
  32. Gerosa, L.; Chidley, C.; Frohlich, F.; Sanchez, G.; Lim, S.K.; Muhlich, J.; Chen, J.Y.; Vallabhaneni, S.; Baker, G.J.; Schapiro, D.; et al. Receptor-Driven ERK Pulses Reconfigure MAPK Signaling and Enable Persistence of Drug-Adapted BRAF-Mutant Melanoma Cells. Cell Syst. 2020, 11, 478–494.e9. [Google Scholar] [CrossRef] [PubMed]
  33. Liu, F.; Heiner, M.; Gilbert, D. Hybrid modelling of biological systems: Current progress and future prospects. Brief. Bioinform. 2022, 23, bbac081. [Google Scholar] [CrossRef] [PubMed]
  34. Faeder, J.R.; Blinov, M.L.; Hlavacek, W.S. Rule-based modeling of biochemical systems with BioNetGen. Methods Mol. Biol. 2009, 500, 113–167. [Google Scholar]
  35. Feret, J.; Danos, V.; Krivine, J.; Harmer, R.; Fontana, W. Internal coarse-graining of molecular systems. Proc. Natl. Acad. Sci. USA 2009, 106, 6453–6458. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Tsai, C.-F.; Wang, Y.-T.; Lin, P.-Y.; Chen, Y.-J. Phosphoproteomics by Highly Selective IMAC Protocol. In Neuroproteomics; Li, K.W., Ed.; Human Press: Totowa, NJ, USA, 2011; Volume 57, pp. 181–196. [Google Scholar]
  37. Karra, A.S.; Stippec, S.; Cobb, M.H. Assaying Protein Kinase Activity with Radiolabeled ATP. J. Vis. Exp. 2017, 123, 55504. [Google Scholar] [CrossRef] [PubMed]
  38. Akita, S.; Umezawa, N.; Kato, N.; Higuchi, T. Array-based fluorescence assay for serine/threonine kinases using specific chemical reaction. Bioorg. Med. Chem. 2008, 16, 7788–7794. [Google Scholar] [CrossRef]
  39. Martin, K.; Steinberg, T.H.; Cooley, L.A.; Gee, K.R.; Beechem, J.M.; Patton, W.F. Quantitative analysis of protein phosphorylation status and protein kinase activity on microarrays using a novel fluorescent phosphorylation sensor dye. Proteomics 2003, 3, 1244–1255. [Google Scholar] [CrossRef]
  40. Blackburn, K.; Goshe, M.B. Challenges and strategies for targeted phosphorylation site identification and quantification using mass spectrometry analysis. Brief. Funct. Genom. Proteom. 2009, 8, 90–103. [Google Scholar] [CrossRef] [Green Version]
  41. Huttlin, E.L.; Jedrychowski, M.P.; Elias, J.E.; Goswami, T.; Rad, R.; Beausoleil, S.A.; Villen, J.; Haas, W.; Sowa, M.E.; Gygi, S.P. A tissue-specific atlas of mouse protein phosphorylation and expression. Cell 2010, 143, 1174–1189. [Google Scholar] [CrossRef] [Green Version]
  42. Sharma, K.; D’Souza, R.C.; Tyanova, S.; Schaab, C.; Wisniewski, J.R.; Cox, J.; Mann, M. Ultradeep human phosphoproteome reveals a distinct regulatory nature of Tyr and Ser/Thr-based signaling. Cell Rep. 2014, 8, 1583–1594. [Google Scholar] [CrossRef] [Green Version]
  43. Aebersold, R.; Goodlett, D.R. Mass spectrometry in proteomics. Chem. Rev. 2001, 101, 269–295. [Google Scholar] [CrossRef] [PubMed]
  44. Osinalde, N.; Aloria, K.; Omaetxebarria, M.J.; Kratchmarova, I. Targeted mass spectrometry: An emerging powerful approach to unblock the bottleneck in phosphoproteomics. J. Chromatogr. B Analyt. Technol. Biomed. Life Sci. 2017, 1055–1056, 29–38. [Google Scholar] [CrossRef] [PubMed]
  45. Von Stechow, L.; Francavilla, C.; Olsen, J.V. Recent findings and technological advances in phosphoproteomics for cells and tissues. Expert Rev. Proteom. 2015, 12, 469–487. [Google Scholar] [CrossRef] [Green Version]
  46. Hornbeck, P.V.; Chabra, I.; Kornhauser, J.M.; Skrzypek, E.; Zhang, B. PhosphoSite: A bioinformatics resource dedicated to physiological protein phosphorylation. Proteomics 2004, 4, 1551–1561. [Google Scholar] [CrossRef]
  47. Gnad, F.; Gunawardena, J.; Mann, M. PHOSIDA 2011: The posttranslational modification database. Nucleic Acids Res. 2011, 39, D253–D260. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Gnad, F.; Ren, S.; Cox, J.; Olsen, J.V.; Macek, B.; Oroshi, M.; Mann, M. PHOSIDA (phosphorylation site database): Management, structural and evolutionary investigation, and prediction of phosphosites. Genome Biol. 2007, 8, R250. [Google Scholar] [CrossRef] [Green Version]
  49. Dinkel, H.; Chica, C.; Via, A.; Gould, C.M.; Jensen, L.J.; Gibson, T.J.; Diella, F. Phospho.ELM: A database of phosphorylation sites--update 2011. Nucleic Acids Res. 2011, 39, D261–D267. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Diella, F.; Cameron, S.; Gemund, C.; Linding, R.; Via, A.; Kuster, B.; Sicheritz-Ponten, T.; Blom, N.; Gibson, T.J. Phospho.ELM: A database of experimentally verified phosphorylation sites in eukaryotic proteins. BMC Bioinform. 2004, 5, 79. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Bodenmiller, B.; Campbell, D.; Gerrits, B.; Lam, H.; Jovanovic, M.; Picotti, P.; Schlapbach, R.; Aebersold, R. PhosphoPep—A database of protein phosphorylation sites in model organisms. Nat. Biotechnol. 2008, 26, 1339–1340. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Bodenmiller, B.; Malmstrom, J.; Gerrits, B.; Campbell, D.; Lam, H.; Schmidt, A.; Rinner, O.; Mueller, L.N.; Shannon, P.T.; Pedrioli, P.G.; et al. PhosphoPep—A phosphoproteome resource for systems biology research in Drosophila Kc167 cells. Mol. Syst. Biol. 2007, 3, 139. [Google Scholar] [CrossRef]
  53. Yi, L.; Shi, T.; Gritsenko, M.A.; X’Avia Chan, C.Y.; Fillmore, T.L.; Hess, B.M.; Swensen, A.C.; Liu, T.; Smith, R.D.; Wiley, H.S.; et al. Targeted Quantification of Phosphorylation Dynamics in the Context of EGFR-MAPK Pathway. Anal. Chem. 2018, 90, 5256–5263. [Google Scholar] [CrossRef] [PubMed]
  54. De Graaf, E.L.; Kaplon, J.; Mohammed, S.; Vereijken, L.A.; Duarte, D.P.; Redondo Gallego, L.; Heck, A.J.; Peeper, D.S.; Altelaar, A.F. Signal Transduction Reaction Monitoring Deciphers Site-Specific PI3K-mTOR/MAPK Pathway Dynamics in Oncogene-Induced Senescence. J. Proteome Res. 2015, 14, 2906–2914. [Google Scholar] [CrossRef] [PubMed]
  55. Lange, V.; Picotti, P.; Domon, B.; Aebersold, R. Selected reaction monitoring for quantitative proteomics: A tutorial. Mol. Syst. Biol. 2008, 4, 222. [Google Scholar] [CrossRef] [PubMed]
  56. Hoofnagle, A.N.; Whiteaker, J.R.; Carr, S.A.; Kuhn, E.; Liu, T.; Massoni, S.A.; Thomas, S.N.; Townsend, R.R.; Zimmerman, L.J.; Boja, E.; et al. Recommendations for the Generation, Quantification, Storage, and Handling of Peptides Used for Mass Spectrometry-Based Assays. Clin. Chem. 2016, 62, 48–69. [Google Scholar] [CrossRef] [Green Version]
  57. Brownridge, P.; Holman, S.W.; Gaskell, S.J.; Grant, C.M.; Harman, V.M.; Hubbard, S.J.; Lanthaler, K.; Lawless, C.; O’Cualain, R.; Sims, P.; et al. Global absolute quantification of a proteome: Challenges in the deployment of a QconCAT strategy. Proteomics 2011, 11, 2957–2970. [Google Scholar] [CrossRef]
  58. Dickhut, C.; Feldmann, I.; Lambert, J.; Zahedi, R.P. Impact of digestion conditions on phosphoproteomics. J. Proteome Res. 2014, 13, 2761–2770. [Google Scholar] [CrossRef]
  59. Bian, Y.; Ye, M.; Song, C.; Cheng, K.; Wang, C.; Wei, X.; Zhu, J.; Chen, R.; Wang, F.; Zou, H. Improve the coverage for the analysis of phosphoproteome of HeLa cells by a tandem digestion approach. J. Proteome Res. 2012, 11, 2828–2837. [Google Scholar] [CrossRef]
  60. Wisniewski, J.R.; Mann, M. Consecutive proteolytic digestion in an enzyme reactor increases depth of proteomic and phosphoproteomic analysis. Anal. Chem. 2012, 84, 2631–2637. [Google Scholar] [CrossRef]
  61. Delom, F.; Chevet, E. Phosphoprotein analysis: From proteins to proteomes. Proteome Sci. 2006, 4, 15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Posewitz, M.C.; Tempst, P. Immobilized gallium(III) affinity chromatography of phosphopeptides. Anal. Chem. 1999, 71, 2883–2892. [Google Scholar] [CrossRef]
  63. Zhou, H.; Ye, M.; Dong, J.; Han, G.; Jiang, X.; Wu, R.; Zou, H. Specific phosphopeptide enrichment with immobilized titanium ion affinity chromatography adsorbent for phosphoproteome analysis. J. Proteome Res. 2008, 7, 3957–3967. [Google Scholar] [CrossRef] [PubMed]
  64. Feng, S.; Ye, M.; Zhou, H.; Jiang, X.; Jiang, X.; Zou, H.; Gong, B. Immobilized zirconium ion affinity chromatography for specific enrichment of phosphopeptides in phosphoproteome analysis. Mol. Cell. Proteom. 2007, 6, 1656–1665. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Andersson, L.; Porath, J. Isolation of phosphoproteins by immobilized metal (Fe3+) affinity chromatography. Anal. Biochem. 1986, 154, 250–254. [Google Scholar] [CrossRef]
  66. Neville, D.C.; Rozanas, C.R.; Price, E.M.; Gruis, D.B.; Verkman, A.S.; Townsend, R.R. Evidence for phosphorylation of serine 753 in CFTR using a novel metal-ion affinity resin and matrix-assisted laser desorption mass spectrometry. Protein Sci. 1997, 6, 2436–2445. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Jensen, S.S.; Larsen, M.R. Evaluation of the impact of some experimental procedures on different phosphopeptide enrichment techniques. Rapid Commun. Mass. Spectrom. 2007, 21, 3635–3645. [Google Scholar] [CrossRef] [PubMed]
  68. Seeley, E.H.; Riggs, L.D.; Regnier, F.E. Reduction of non-specific binding in Ga(III) immobilized metal affinity chromatography for phosphopeptides by using endoproteinase glu-C as the digestive enzyme. J. Chromatogr. B Analyt. Technol. Biomed. Life Sci. 2005, 817, 81–88. [Google Scholar] [CrossRef] [PubMed]
  69. Ficarro, S.B.; McCleland, M.L.; Stukenberg, P.T.; Burke, D.J.; Ross, M.M.; Shabanowitz, J.; Hunt, D.F.; White, F.M. Phosphoproteome analysis by mass spectrometry and its application to Saccharomyces cerevisiae. Nat. Biotechnol. 2002, 20, 301–305. [Google Scholar] [CrossRef]
  70. Larsen, M.R.; Thingholm, T.E.; Jensen, O.N.; Roepstorff, P.; Jorgensen, T.J. Highly selective enrichment of phosphorylated peptides from peptide mixtures using titanium dioxide microcolumns. Mol. Cell. Proteom. 2005, 4, 873–886. [Google Scholar] [CrossRef] [Green Version]
  71. Li, Y.; Luo, Y.; Wu, S.; Gao, Y.; Liu, Y.; Zheng, D. Nucleic acids in protein samples interfere with phosphopeptide identification by immobilized-metal-ion affinity chromatography and mass spectrometry. Mol. Biotechnol. 2009, 43, 59–66. [Google Scholar] [CrossRef]
  72. Ficarro, S.B.; Adelmant, G.; Tomar, M.N.; Zhang, Y.; Cheng, V.J.; Marto, J.A. Magnetic bead processor for rapid evaluation and optimization of parameters for phosphopeptide enrichment. Anal. Chem. 2009, 81, 4566–4575. [Google Scholar] [CrossRef] [Green Version]
  73. Zhou, H.; Xu, S.; Ye, M.; Feng, S.; Pan, C.; Jiang, X.; Li, X.; Han, G.; Fu, Y.; Zou, H. Zirconium phosphonate-modified porous silicon for highly specific capture of phosphopeptides and MALDI-TOF MS analysis. J. Proteome Res. 2006, 5, 2431–2437. [Google Scholar] [CrossRef] [PubMed]
  74. Lai, A.C.; Tsai, C.F.; Hsu, C.C.; Sun, Y.N.; Chen, Y.J. Complementary Fe(3+)- and Ti(4+)-immobilized metal ion affinity chromatography for purification of acidic and basic phosphopeptides. Rapid Commun. Mass. Spectrom. 2012, 26, 2186–2194. [Google Scholar] [CrossRef] [PubMed]
  75. Tsai, C.F.; Hsu, C.C.; Hung, J.N.; Wang, Y.T.; Choong, W.K.; Zeng, M.Y.; Lin, P.Y.; Hong, R.W.; Sung, T.Y.; Chen, Y.J. Sequential phosphoproteomic enrichment through complementary metal-directed immobilized metal ion affinity chromatography. Anal. Chem. 2014, 86, 685–693. [Google Scholar] [CrossRef]
  76. Leitner, A. Phosphopeptide enrichment using metal oxide affinity chromatography. TrAC Trends Anal. Chem. 2010, 29, 177–185. [Google Scholar] [CrossRef]
  77. Pinkse, M.W.; Uitto, P.M.; Hilhorst, M.J.; Ooms, B.; Heck, A.J. Selective isolation at the femtomole level of phosphopeptides from proteolytic digests using 2D-NanoLC-ESI-MS/MS and titanium oxide precolumns. Anal. Chem. 2004, 76, 3935–3943. [Google Scholar] [CrossRef]
  78. Sugiyama, N.; Masuda, T.; Shinoda, K.; Nakamura, A.; Tomita, M.; Ishihama, Y. Phosphopeptide enrichment by aliphatic hydroxy acid-modified metal oxide chromatography for nano-LC-MS/MS in proteomics applications. Mol. Cell. Proteom. 2007, 6, 1103–1109. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Wakabayashi, M.; Kyono, Y.; Sugiyama, N.; Ishihama, Y. Extended Coverage of Singly and Multiply Phosphorylated Peptides from a Single Titanium Dioxide Microcolumn. Anal. Chem. 2015, 87, 10213–10221. [Google Scholar] [CrossRef] [PubMed]
  80. Yue, X.; Schunter, A.; Hummon, A.B. Comparing multistep immobilized metal affinity chromatography and multistep TiO2 methods for phosphopeptide enrichment. Anal. Chem. 2015, 87, 8837–8844. [Google Scholar] [CrossRef] [Green Version]
  81. Thingholm, T.E.; Jensen, O.N.; Robinson, P.J.; Larsen, M.R. SIMAC (sequential elution from IMAC), a phosphoproteomics strategy for the rapid separation of monophosphorylated from multiply phosphorylated peptides. Mol. Cell. Proteom. 2008, 7, 661–671. [Google Scholar] [CrossRef] [Green Version]
  82. Carrascal, M.; Ovelleiro, D.; Casas, V.; Gay, M.; Abian, J. Phosphorylation analysis of primary human T lymphocytes using sequential IMAC and titanium oxide enrichment. J. Proteome Res. 2008, 7, 5167–5176. [Google Scholar] [CrossRef]
  83. Choi, J.; Snovida, S.; Bomgarden, R.D.; Rogers, J.C. Sequential enrichment from Metal Oxide Affinity Chromatography (SMOAC), a phosphoproteomics strategy for the separation of multiply phosphorylated from monophosphorylated peptides. Am. Soc. Mass Spectrom. Conf. 2017, WP601. [Google Scholar]
  84. Iliuk, A.; Jayasundera, K.; Wang, W.H.; Schluttenhofer, R.; Geahlen, R.L.; Tao, W.A. In-Depth Analyses of B Cell Signaling Through Tandem Mass Spectrometry of Phosphopeptides Enriched by PolyMAC. Int. J. Mass. Spectrom. 2015, 377, 744–753. [Google Scholar] [CrossRef]
  85. Iliuk, A.B.; Martin, V.A.; Alicie, B.M.; Geahlen, R.L.; Tao, W.A. In-depth analyses of kinase-dependent tyrosine phosphoproteomes based on metal ion-functionalized soluble nanopolymers. Mol. Cell. Proteom. 2010, 9, 2162–2172. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Iliuk, A.; Jayasundera, K.; Schluttenhofer, R.; Tao, W.A. Functionalized soluble nanopolymers for phosphoproteome analysis. Methods Mol. Biol. 2011, 790, 277–285. [Google Scholar]
  87. Ding, S.J.; Qian, W.J.; Smith, R.D. Quantitative proteomic approaches for studying phosphotyrosine signaling. Expert Rev. Proteom. 2007, 4, 13–23. [Google Scholar] [CrossRef] [PubMed]
  88. Rush, J.; Moritz, A.; Lee, K.A.; Guo, A.; Goss, V.L.; Spek, E.J.; Zhang, H.; Zha, X.M.; Polakiewicz, R.D.; Comb, M.J. Immunoaffinity profiling of tyrosine phosphorylation in cancer cells. Nat. Biotechnol. 2005, 23, 94–101. [Google Scholar] [CrossRef] [PubMed]
  89. Rikova, K.; Guo, A.; Zeng, Q.; Possemato, A.; Yu, J.; Haack, H.; Nardone, J.; Lee, K.; Reeves, C.; Li, Y.; et al. Global survey of phosphotyrosine signaling identifies oncogenic kinases in lung cancer. Cell 2007, 131, 1190–1203. [Google Scholar] [CrossRef] [Green Version]
  90. Blagoev, B.; Ong, S.E.; Kratchmarova, I.; Mann, M. Temporal analysis of phosphotyrosine-dependent signaling networks by quantitative proteomics. Nat. Biotechnol. 2004, 22, 1139–1145. [Google Scholar] [CrossRef]
  91. Heibeck, T.H.; Ding, S.J.; Opresko, L.K.; Zhao, R.; Schepmoes, A.A.; Yang, F.; Tolmachev, A.V.; Monroe, M.E.; Camp, D.G., 2nd; Smith, R.D.; et al. An extensive survey of tyrosine phosphorylation revealing new sites in human mammary epithelial cells. J. Proteome Res. 2009, 8, 3852–3861. [Google Scholar] [CrossRef] [Green Version]
  92. Steen, H.; Kuster, B.; Fernandez, M.; Pandey, A.; Mann, M. Tyrosine phosphorylation mapping of the epidermal growth factor receptor signaling pathway. J. Biol. Chem. 2002, 277, 1031–1039. [Google Scholar] [CrossRef] [Green Version]
  93. Keshishian, H.; McDonald, E.R., 3rd; Mundt, F.; Melanson, R.; Krug, K.; Porter, D.A.; Wallace, L.; Forestier, D.; Rabasha, B.; Marlow, S.E.; et al. A highly multiplexed quantitative phosphosite assay for biology and preclinical studies. Mol. Syst. Biol. 2021, 17, e10156. [Google Scholar] [CrossRef] [PubMed]
  94. Fila, J.; Honys, D. Enrichment techniques employed in phosphoproteomics. Amino Acids 2012, 43, 1025–1047. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  95. Beausoleil, S.A.; Jedrychowski, M.; Schwartz, D.; Elias, J.E.; Villen, J.; Li, J.; Cohn, M.A.; Cantley, L.C.; Gygi, S.P. Large-scale characterization of HeLa cell nuclear phosphoproteins. Proc. Natl. Acad. Sci. USA 2004, 101, 12130–12135. [Google Scholar] [CrossRef]
  96. Gruhler, A.; Olsen, J.V.; Mohammed, S.; Mortensen, P.; Faergeman, N.J.; Mann, M.; Jensen, O.N. Quantitative phosphoproteomics applied to the yeast pheromone signaling pathway. Mol. Cell. Proteom. 2005, 4, 310–327. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  97. Nie, S.; Dai, J.; Ning, Z.B.; Cao, X.J.; Sheng, Q.H.; Zeng, R. Comprehensive profiling of phosphopeptides based on anion exchange followed by flow-through enrichment with titanium dioxide (AFET). J. Proteome Res. 2010, 9, 4585–4594. [Google Scholar] [CrossRef]
  98. Han, G.; Ye, M.; Zhou, H.; Jiang, X.; Feng, S.; Jiang, X.; Tian, R.; Wan, D.; Zou, H.; Gu, J. Large-scale phosphoproteome analysis of human liver tissue by enrichment and fractionation of phosphopeptides with strong anion exchange chromatography. Proteomics 2008, 8, 1346–1361. [Google Scholar] [CrossRef] [PubMed]
  99. Rosenbaum, D. Quantitative analysis of a multicomponent analgesic product containing butalbital, using high speed reverse-phase liquid chromatography. Anal. Chem. 1974, 46, 2226–2228. [Google Scholar] [CrossRef] [PubMed]
  100. Gilar, M.; Olivova, P.; Daly, A.E.; Gebler, J.C. Orthogonality of separation in two-dimensional liquid chromatography. Anal. Chem. 2005, 77, 6426–6434. [Google Scholar] [CrossRef]
  101. Wang, Y.; Yang, F.; Gritsenko, M.A.; Wang, Y.; Clauss, T.; Liu, T.; Shen, Y.; Monroe, M.E.; Lopez-Ferrer, D.; Reno, T.; et al. Reversed-phase chromatography with multiple fraction concatenation strategy for proteome profiling of human MCF10A cells. Proteomics 2011, 11, 2019–2026. [Google Scholar] [CrossRef] [Green Version]
  102. Song, C.; Ye, M.; Han, G.; Jiang, X.; Wang, F.; Yu, Z.; Chen, R.; Zou, H. Reversed-phase-reversed-phase liquid chromatography approach with high orthogonality for multidimensional separation of phosphopeptides. Anal. Chem. 2010, 82, 53–56. [Google Scholar] [CrossRef]
  103. Shi, T.; Fillmore, T.L.; Sun, X.; Zhao, R.; Schepmoes, A.A.; Hossain, M.; Xie, F.; Wu, S.; Kim, J.S.; Jones, N.; et al. Antibody-free, targeted mass-spectrometric approach for quantification of proteins at low picogram per milliliter levels in human plasma/serum. Proc. Natl. Acad. Sci. USA 2012, 109, 15395–15400. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  104. Shi, T.; Gao, Y.; Gaffrey, M.J.; Nicora, C.D.; Fillmore, T.L.; Chrisler, W.B.; Gritsenko, M.A.; Wu, C.; He, J.; Bloodsworth, K.J.; et al. Sensitive targeted quantification of ERK phosphorylation dynamics and stoichiometry in human cells without affinity enrichment. Anal. Chem. 2015, 87, 1103–1110. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  105. Mertins, P.; Qiao, J.W.; Patel, J.; Udeshi, N.D.; Clauser, K.R.; Mani, D.R.; Burgess, M.W.; Gillette, M.A.; Jaffe, J.D.; Carr, S.A. Integrated proteomic analysis of post-translational modifications by serial enrichment. Nat. Methods 2013, 10, 634–637. [Google Scholar] [CrossRef] [PubMed]
  106. Mertins, P.; Yang, F.; Liu, T.; Mani, D.R.; Petyuk, V.A.; Gillette, M.A.; Clauser, K.R.; Qiao, J.W.; Gritsenko, M.A.; Moore, R.J.; et al. Ischemia in tumors induces early and sustained phosphorylation changes in stress kinase pathways but does not affect global protein levels. Mol. Cell. Proteom. 2014, 13, 1690–1704. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  107. Yue, X.S.; Hummon, A.B. Combination of multistep IMAC enrichment with high-pH reverse phase separation for in-depth phosphoproteomic profiling. J. Proteome Res. 2013, 12, 4176–4186. [Google Scholar] [CrossRef] [PubMed]
  108. Villen, J.; Gygi, S.P. The SCX/IMAC enrichment approach for global phosphorylation analysis by mass spectrometry. Nat. Protoc. 2008, 3, 1630–1638. [Google Scholar] [CrossRef]
  109. McNulty, D.E.; Annan, R.S. Hydrophilic interaction chromatography reduces the complexity of the phosphoproteome and improves global phosphopeptide isolation and detection. Mol. Cell. Proteom. 2008, 7, 971–980. [Google Scholar] [CrossRef] [Green Version]
  110. Alpert, A.J. Electrostatic repulsion hydrophilic interaction chromatography for isocratic separation of charged solutes and selective isolation of phosphopeptides. Anal. Chem. 2008, 80, 62–76. [Google Scholar] [CrossRef] [Green Version]
  111. Chien, K.Y.; Liu, H.C.; Goshe, M.B. Development and application of a phosphoproteomic method using electrostatic repulsion-hydrophilic interaction chromatography (ERLIC), IMAC, and LC-MS/MS analysis to study Marek’s Disease Virus infection. J. Proteome Res. 2011, 10, 4041–4053. [Google Scholar] [CrossRef]
  112. Method of the Year 2012. Nat. Methods 2013, 10, 1. [CrossRef]
  113. Addona, T.A.; Abbatiello, S.E.; Schilling, B.; Skates, S.J.; Mani, D.R.; Bunk, D.M.; Spiegelman, C.H.; Zimmerman, L.J.; Ham, A.J.; Keshishian, H.; et al. Multi-site assessment of the precision and reproducibility of multiple reaction monitoring-based measurements of proteins in plasma. Nat. Biotechnol. 2009, 27, 633–641. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  114. Picotti, P.; Aebersold, R. Selected reaction monitoring-based proteomics: Workflows, potential, pitfalls and future directions. Nat. Methods 2012, 9, 555–566. [Google Scholar] [CrossRef] [PubMed]
  115. Friedrich Lottspeich, J.E. Analysis of Post-translational Modifications: Phosphorylation and Acetylation of Proteins. In Bioanalytics: Analytical Methods and Concepts in Biochemistry and Molecular Biology; Friedrich Lottspeich, J.E., Ed.; Wiley-VCH: Weinheim, Germany, 2018; pp. 645–662. [Google Scholar]
  116. Steen, H.; Mann, M. The ABC’s (and XYZ’s) of peptide sequencing. Nat. Rev. Mol. Cell Biol. 2004, 5, 699–711. [Google Scholar] [CrossRef] [PubMed]
  117. MacLean, B.; Tomazela, D.M.; Shulman, N.; Chambers, M.; Finney, G.L.; Frewen, B.; Kern, R.; Tabb, D.L.; Liebler, D.C.; MacCoss, M.J. Skyline: An open source document editor for creating and analyzing targeted proteomics experiments. Bioinformatics 2010, 26, 966–968. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  118. Martin, D.B.; Holzman, T.; May, D.; Peterson, A.; Eastham, A.; Eng, J.; McIntosh, M. MRMer, an interactive open source and cross-platform system for data extraction and visualization of multiple reaction monitoring experiments. Mol. Cell. Proteom. 2008, 7, 2270–2278. [Google Scholar] [CrossRef] [Green Version]
  119. Wolf-Yadlin, A.; Hautaniemi, S.; Lauffenburger, D.A.; White, F.M. Multiple reaction monitoring for robust quantitative proteomic analysis of cellular signaling networks. Proc. Natl. Acad. Sci. USA 2007, 104, 5860–5865. [Google Scholar] [CrossRef] [Green Version]
  120. Kirkpatrick, D.S.; Gerber, S.A.; Gygi, S.P. The absolute quantification strategy: A general procedure for the quantification of proteins and post-translational modifications. Methods 2005, 35, 265–273. [Google Scholar] [CrossRef]
  121. Gerber, S.A.; Rush, J.; Stemman, O.; Kirschner, M.W.; Gygi, S.P. Absolute quantification of proteins and phosphoproteins from cell lysates by tandem MS. Proc. Natl. Acad. Sci. USA 2003, 100, 6940–6945. [Google Scholar] [CrossRef] [Green Version]
  122. Wu, R.; Dephoure, N.; Haas, W.; Huttlin, E.L.; Zhai, B.; Sowa, M.E.; Gygi, S.P. Correct interpretation of comprehensive phosphorylation dynamics requires normalization by protein expression changes. Mol. Cell. Proteom. 2011, 10, M111009654. [Google Scholar] [CrossRef] [Green Version]
  123. O’Rourke, M.B.; Town, S.E.L.; Dalla, P.V.; Bicknell, F.; Koh Belic, N.; Violi, J.P.; Steele, J.R.; Padula, M.P. What is Normalization? The Strategies Employed in Top-Down and Bottom-Up Proteome Analysis Workflows. Proteomes 2019, 7, 29. [Google Scholar] [CrossRef] [Green Version]
  124. Narumi, R.; Murakami, T.; Kuga, T.; Adachi, J.; Shiromizu, T.; Muraoka, S.; Kume, H.; Kodera, Y.; Matsumoto, M.; Nakayama, K.; et al. A strategy for large-scale phosphoproteomics and SRM-based validation of human breast cancer tissue samples. J. Proteome Res. 2012, 11, 5311–5322. [Google Scholar] [CrossRef]
  125. Schmidlin, T.; Garrigues, L.; Lane, C.S.; Mulder, T.C.; van Doorn, S.; Post, H.; de Graaf, E.L.; Lemeer, S.; Heck, A.J.; Altelaar, A.F. Assessment of SRM, MRM(3), and DIA for the targeted analysis of phosphorylation dynamics in non-small cell lung cancer. Proteomics 2016, 16, 2193–2205. [Google Scholar] [CrossRef] [PubMed]
  126. Abelin, J.G.; Patel, J.; Lu, X.; Feeney, C.M.; Fagbami, L.; Creech, A.L.; Hu, R.; Lam, D.; Davison, D.; Pino, L.; et al. Reduced-representation Phosphosignatures Measured by Quantitative Targeted MS Capture Cellular States and Enable Large-scale Comparison of Drug-induced Phenotypes. Mol. Cell. Proteom. 2016, 15, 1622–1641. [Google Scholar] [CrossRef] [PubMed]
  127. Kennedy, J.J.; Yan, P.; Zhao, L.; Ivey, R.G.; Voytovich, U.J.; Moore, H.D.; Lin, C.; Pogosova-Agadjanyan, E.L.; Stirewalt, D.L.; Reding, K.W.; et al. Immobilized Metal Affinity Chromatography Coupled to Multiple Reaction Monitoring Enables Reproducible Quantification of Phospho-signaling. Mol. Cell. Proteom. 2016, 15, 726–739. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  128. Whiteaker, J.R.; Zhao, L.; Yan, P.; Ivey, R.G.; Voytovich, U.J.; Moore, H.D.; Lin, C.; Paulovich, A.G. Peptide Immunoaffinity Enrichment and Targeted Mass Spectrometry Enables Multiplex, Quantitative Pharmacodynamic Studies of Phospho-Signaling. Mol. Cell. Proteom. 2015, 14, 2261–2273. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  129. Whiteaker, J.R.; Zhao, L.; Saul, R.; Kaczmarczyk, J.A.; Schoenherr, R.M.; Moore, H.D.; Jones-Weinert, C.; Ivey, R.G.; Lin, C.; Hiltke, T.; et al. A Multiplexed Mass Spectrometry-Based Assay for Robust Quantification of Phosphosignaling in Response to DNA Damage. Radiat. Res. 2018, 189, 505–518. [Google Scholar] [CrossRef] [PubMed]
  130. Bourmaud, A.; Gallien, S.; Domon, B. Parallel reaction monitoring using quadrupole-Orbitrap mass spectrometer: Principle and applications. Proteomics 2016, 16, 2146–2159. [Google Scholar] [CrossRef] [PubMed]
  131. Gallien, S.; Duriez, E.; Crone, C.; Kellmann, M.; Moehring, T.; Domon, B. Targeted proteomic quantification on quadrupole-orbitrap mass spectrometer. Mol. Cell. Proteom. 2012, 11, 1709–1723. [Google Scholar] [CrossRef] [Green Version]
  132. Peterson, A.C.; Russell, J.D.; Bailey, D.J.; Westphall, M.S.; Coon, J.J. Parallel reaction monitoring for high resolution and high mass accuracy quantitative, targeted proteomics. Mol. Cell. Proteom. 2012, 11, 1475–1488. [Google Scholar] [CrossRef] [Green Version]
  133. Kang, T.; Jensen, P.; Huang, H.; Lund Christensen, G.; Billestrup, N.; Larsen, M.R. Characterization of the Molecular Mechanisms Underlying Glucose Stimulated Insulin Secretion from Isolated Pancreatic beta-cells Using Post-translational Modification Specific Proteomics (PTMomics). Mol. Cell. Proteom. 2018, 17, 95–110. [Google Scholar] [CrossRef] [Green Version]
  134. Morales Betanzos, C.; Federspiel, J.D.; Palubinsky, A.M.; McLaughlin, B.; Liebler, D.C. Dynamic Phosphorylation of Apoptosis Signal Regulating Kinase 1 (ASK1) in Response to Oxidative and Electrophilic Stress. Chem. Res. Toxicol. 2016, 29, 2175–2183. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  135. Osinalde, N.; Mitxelena, J.; Sanchez-Quiles, V.; Akimov, V.; Aloria, K.; Arizmendi, J.M.; Zubiaga, A.M.; Blagoev, B.; Kratchmarova, I. Nuclear Phosphoproteomic Screen Uncovers ACLY as Mediator of IL-2-induced Proliferation of CD4+ T lymphocytes. Mol. Cell. Proteom. 2016, 15, 2076–2092. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  136. Gallien, S.; Kim, S.Y.; Domon, B. Large-Scale Targeted Proteomics Using Internal Standard Triggered-Parallel Reaction Monitoring (IS-PRM). Mol. Cell. Proteom. 2015, 14, 1630–1644. [Google Scholar] [CrossRef] [PubMed]
  137. Ludwig, C.; Gillet, L.; Rosenberger, G.; Amon, S.; Collins, B.C.; Aebersold, R. Data-independent acquisition-based SWATH-MS for quantitative proteomics: A tutorial. Mol. Syst. Biol. 2018, 14, e8126. [Google Scholar] [CrossRef] [PubMed]
  138. Law, K.P.; Lim, Y.P. Recent advances in mass spectrometry: Data independent analysis and hyper reaction monitoring. Expert Rev. Proteom. 2013, 10, 551–566. [Google Scholar] [CrossRef] [PubMed]
  139. Gillet, L.C.; Navarro, P.; Tate, S.; Rost, H.; Selevsek, N.; Reiter, L.; Bonner, R.; Aebersold, R. Targeted data extraction of the MS/MS spectra generated by data-independent acquisition: A new concept for consistent and accurate proteome analysis. Mol. Cell. Proteom. 2012, 11, O111016717. [Google Scholar] [CrossRef] [Green Version]
  140. Kitata, R.B.; Choong, W.K.; Tsai, C.F.; Lin, P.Y.; Chen, B.S.; Chang, Y.C.; Nesvizhskii, A.I.; Sung, T.Y.; Chen, Y.J. A data-independent acquisition-based global phosphoproteomics system enables deep profiling. Nat. Commun. 2021, 12, 2539. [Google Scholar] [CrossRef]
  141. Rost, H.L.; Rosenberger, G.; Navarro, P.; Gillet, L.; Miladinovic, S.M.; Schubert, O.T.; Wolski, W.; Collins, B.C.; Malmstrom, J.; Malmstrom, L.; et al. OpenSWATH enables automated, targeted analysis of data-independent acquisition MS data. Nat. Biotechnol. 2014, 32, 219–223. [Google Scholar] [CrossRef] [Green Version]
  142. Demichev, V.; Messner, C.B.; Vernardis, S.I.; Lilley, K.S.; Ralser, M. DIA-NN: Neural networks and interference correction enable deep proteome coverage in high throughput. Nat. Methods 2020, 17, 41–44. [Google Scholar] [CrossRef]
  143. Roux, P.P.; Shahbazian, D.; Vu, H.; Holz, M.K.; Cohen, M.S.; Taunton, J.; Sonenberg, N.; Blenis, J. RAS/ERK signaling promotes site-specific ribosomal protein S6 phosphorylation via RSK and stimulates cap-dependent translation. J. Biol. Chem. 2007, 282, 14056–14064. [Google Scholar] [CrossRef] [Green Version]
  144. Corcoran, R.B.; Rothenberg, S.M.; Hata, A.N.; Faber, A.C.; Piris, A.; Nazarian, R.M.; Brown, R.D.; Godfrey, J.T.; Winokur, D.; Walsh, J.; et al. TORC1 suppression predicts responsiveness to RAF and MEK inhibition in BRAF-mutant melanoma. Sci. Transl. Med. 2013, 5, 196ra98. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  145. Qian, Z.R.; Ter-Minassian, M.; Chan, J.A.; Imamura, Y.; Hooshmand, S.M.; Kuchiba, A.; Morikawa, T.; Brais, L.K.; Daskalova, A.; Heafield, R.; et al. Prognostic significance of MTOR pathway component expression in neuroendocrine tumors. J. Clin. Oncol. 2013, 31, 3418–3425. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  146. Bodenmiller, B.; Mueller, L.N.; Mueller, M.; Domon, B.; Aebersold, R. Reproducible isolation of distinct, overlapping segments of the phosphoproteome. Nat. Methods 2007, 4, 231–237. [Google Scholar] [CrossRef] [PubMed]
  147. Dephoure, N.; Zhou, C.; Villen, J.; Beausoleil, S.A.; Bakalarski, C.E.; Elledge, S.J.; Gygi, S.P. A quantitative atlas of mitotic phosphorylation. Proc. Natl. Acad. Sci. USA 2008, 105, 10762–10767. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  148. Guo, A.; Villen, J.; Kornhauser, J.; Lee, K.A.; Stokes, M.P.; Rikova, K.; Possemato, A.; Nardone, J.; Innocenti, G.; Wetzel, R.; et al. Signaling networks assembled by oncogenic EGFR and c-Met. Proc. Natl. Acad. Sci. USA 2008, 105, 692–697. [Google Scholar] [CrossRef] [Green Version]
  149. Choudhary, C.; Olsen, J.V.; Brandts, C.; Cox, J.; Reddy, P.N.; Bohmer, F.D.; Gerke, V.; Schmidt-Arras, D.E.; Berdel, W.E.; Muller-Tidow, C.; et al. Mislocalized activation of oncogenic RTKs switches downstream signaling outcomes. Mol. Cell 2009, 36, 326–339. [Google Scholar] [CrossRef]
  150. Cox, J.; Mann, M. Quantitative, high-resolution proteomics for data-driven systems biology. Annu. Rev. Biochem. 2011, 80, 273–299. [Google Scholar] [CrossRef]
  151. Perez Verdaguer, M.; Zhang, T.; Paulo, J.A.; Gygi, S.; Watkins, S.C.; Sakurai, H.; Sorkin, A. Mechanism of p38 MAPK-induced EGFR endocytosis and its crosstalk with ligand-induced pathways. J. Cell Biol. 2021, 220, e202102005. [Google Scholar] [CrossRef]
  152. Pinilla-Macua, I.; Watkins, S.C.; Sorkin, A. Endocytosis separates EGF receptors from endogenous fluorescently labeled HRas and diminishes receptor signaling to MAP kinases in endosomes. Proc. Natl. Acad. Sci. USA 2016, 113, 2122–2127. [Google Scholar] [CrossRef] [Green Version]
  153. Kuwada, S.K.; Lund, K.A.; Li, X.F.; Cliften, P.; Amsler, K.; Opresko, L.K.; Wiley, H.S. Differential signaling and regulation of apical vs. basolateral EGFR in polarized epithelial cells. Am. J. Physiol. 1998, 275, C1419–C1428. [Google Scholar] [CrossRef]
  154. Salazar-Cavazos, E.; Nitta, C.F.; Mitra, E.D.; Wilson, B.S.; Lidke, K.A.; Hlavacek, W.S.; Lidke, D.S. Multisite EGFR phosphorylation is regulated by adaptor protein abundances and dimer lifetimes. Mol. Biol. Cell 2020, 31, 695–708. [Google Scholar] [CrossRef] [PubMed]
  155. Cirit, M.; Wang, C.C.; Haugh, J.M. Systematic quantification of negative feedback mechanisms in the extracellular signal-regulated kinase (ERK) signaling network. J. Biol. Chem. 2010, 285, 36736–36744. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  156. Kiel, C.; Verschueren, E.; Yang, J.S.; Serrano, L. Integration of protein abundance and structure data reveals competition in the ErbB signaling network. Sci. Signal. 2013, 6, ra109. [Google Scholar] [CrossRef] [PubMed]
  157. Hlavacek, W.S.; Faeder, J.R.; Blinov, M.L.; Perelson, A.S.; Goldstein, B. The complexity of complexes in signal transduction. Biotechnol. Bioeng. 2003, 84, 783–794. [Google Scholar] [CrossRef]
  158. Saez-Rodriguez, C.T.J. Modeling Signaling Networks Using High-throughput Phospho-proteomics. In Advances in Systems Biology; Goryachev, I.I.G.A.B., Ed.; Springer: Berlin/Heidelberg, Germany, 2011; Volume 736, pp. 19–57. [Google Scholar]
  159. Pinilla-Macua, I.; Grassart, A.; Duvvuri, U.; Watkins, S.C.; Sorkin, A. EGF receptor signaling, phosphorylation, ubiquitylation and endocytosis in tumors in vivo. eLife 2017, 6, e31993. [Google Scholar] [CrossRef]
  160. Babur, O.; Luna, A.; Korkut, A.; Durupinar, F.; Siper, M.C.; Dogrusoz, U.; Vaca Jacome, A.S.; Peckner, R.; Christianson, K.E.; Jaffe, J.D.; et al. Causal interactions from proteomic profiles: Molecular data meet pathway knowledge. Patterns 2021, 2, 100257. [Google Scholar] [CrossRef]
  161. Koch, Y.; Wolf, T.; Sorger, P.K.; Eils, R.; Brors, B. Decision-tree based model analysis for efficient identification of parameter relations leading to different signaling states. PLoS One 2013, 8, e82593. [Google Scholar] [CrossRef] [Green Version]
  162. Chen, W.W.; Schoeberl, B.; Jasper, P.J.; Niepel, M.; Nielsen, U.B.; Lauffenburger, D.A.; Sorger, P.K. Input-output behavior of ErbB signaling pathways as revealed by a mass action model trained against dynamic data. Mol. Syst. Biol. 2009, 5, 239. [Google Scholar] [CrossRef]
  163. Stepanov, A.; Karelina, T.; Markevich, N.; Demin, O.; Nicholas, T. A mathematical model of multisite phosphorylation of tau protein. PLoS One 2018, 13, e0192519. [Google Scholar]
  164. Santos, S.D.; Verveer, P.J.; Bastiaens, P.I. Growth factor-induced MAPK network topology shapes Erk response determining PC-12 cell fate. Nat. Cell Biol. 2007, 9, 324–330. [Google Scholar] [CrossRef]
  165. Chylek, L.A.; Akimov, V.; Dengjel, J.; Rigbolt, K.T.; Hu, B.; Hlavacek, W.S.; Blagoev, B. Phosphorylation site dynamics of early T-cell receptor signaling. PLoS One 2014, 9, e104240. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  166. Leonidas, G.; Alexopoulos, J.S.-R.; Espelin, C.W. High-Throughput Protein-Based Technologies and Computational Models for Drug Development, Efficacy, and Toxicity. In Drug Efficacy, Safety, and Biologics Discovery: Emerging Technologies and Tools; Sean Ekins, J.J.X., Ed.; John Wiley and Sons, Inc.: Hoboken, NJ, USA, 2008; pp. 29–52. [Google Scholar]
  167. Nelander, S.; Wang, W.; Nilsson, B.; She, Q.B.; Pratilas, C.; Rosen, N.; Gennemark, P.; Sander, C. Models from experiments: Combinatorial drug perturbations of cancer cells. Mol. Syst. Biol. 2008, 4, 216. [Google Scholar] [CrossRef] [PubMed]
  168. Aldridge, B.B.; Burke, J.M.; Lauffenburger, D.A.; Sorger, P.K. Physicochemical modelling of cell signalling pathways. Nat. Cell Biol. 2006, 8, 1195–1203. [Google Scholar] [CrossRef]
  169. Birtwistle, M.R.; Kholodenko, B.N. Endocytosis and signalling: A meeting with mathematics. Mol. Oncol. 2009, 3, 308–320. [Google Scholar] [CrossRef] [PubMed]
  170. Kamioka, Y.; Yasuda, S.; Fujita, Y.; Aoki, K.; Matsuda, M. Multiple decisive phosphorylation sites for the negative feedback regulation of SOS1 via ERK. J. Biol. Chem. 2010, 285, 33540–33548. [Google Scholar] [CrossRef] [Green Version]
  171. Mayer, B.J.; Blinov, M.L.; Loew, L.M. Molecular machines or pleiomorphic ensembles: Signaling complexes revisited. J. Biol. 2009, 8, 81. [Google Scholar] [CrossRef]
  172. Borisov, N.M.; Markevich, N.I.; Hoek, J.B.; Kholodenko, B.N. Signaling through receptors and scaffolds: Independent interactions reduce combinatorial complexity. Biophys. J. 2005, 89, 951–966. [Google Scholar] [CrossRef] [Green Version]
  173. Sneddon, M.W.; Faeder, J.R.; Emonet, T. Efficient modeling, simulation and coarse-graining of biological complexity with NFsim. Nat. Methods 2011, 8, 177–183. [Google Scholar] [CrossRef]
  174. Chylek, L.A.; Harris, L.A.; Faeder, J.R.; Hlavacek, W.S. Modeling for (physical) biologists: An introduction to the rule-based approach. Phys. Biol. 2015, 12, 045007. [Google Scholar] [CrossRef] [Green Version]
  175. Chylek, L.A.; Stites, E.C.; Posner, R.G.; Hlavacek, W.S. Innovations of the Rule-Based Modeling Approach. In Systems Biology: Integrative Biology and Simulation Tools; Csukás, A.P.B., Ed.; Springer: Berlin/Heidelberg, Germany, 2013; Volume 1, pp. 273–300. [Google Scholar]
  176. Chylek, L.A.; Harris, L.A.; Tung, C.S.; Faeder, J.R.; Lopez, C.F.; Hlavacek, W.S. Rule-based modeling: A computational approach for studying biomolecular site dynamics in cell signaling systems. Wiley Interdiscip. Rev. Syst. Biol. Med. 2014, 6, 13–36. [Google Scholar] [CrossRef] [Green Version]
  177. Hogg, J.S.; Harris, L.A.; Stover, L.J.; Nair, N.S.; Faeder, J.R. Exact hybrid particle/population simulation of rule-based models of biochemical systems. PLoS Comput. Biol. 2014, 10, e1003544. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  178. Colvin, J.; Monine, M.I.; Faeder, J.R.; Hlavacek, W.S.; Von Hoff, D.D.; Posner, R.G. Simulation of large-scale rule-based models. Bioinformatics 2009, 25, 910–917. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  179. Colvin, J.; Monine, M.I.; Gutenkunst, R.N.; Hlavacek, W.S.; Von Hoff, D.D.; Posner, R.G. RuleMonkey: Software for stochastic simulation of rule-based models. BMC Bioinform. 2010, 11, 404. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  180. Lopez, C.F.; Muhlich, J.L.; Bachman, J.A.; Sorger, P.K. Programming biological models in Python using PySB. Mol. Syst. Biol. 2013, 9, 646. [Google Scholar] [CrossRef]
  181. Shi, T.; Niepel, M.; McDermott, J.E.; Gao, Y.; Nicora, C.D.; Chrisler, W.B.; Markillie, L.M.; Petyuk, V.A.; Smith, R.D.; Rodland, K.D.; et al. Conservation of protein abundance patterns reveals the regulatory architecture of the EGFR-MAPK pathway. Sci. Signal. 2016, 9, rs6. [Google Scholar] [CrossRef] [Green Version]
  182. Saha, M.; Carriere, A.; Cheerathodi, M.; Zhang, X.; Lavoie, G.; Rush, J.; Roux, P.P.; Ballif, B.A. RSK phosphorylates SOS1 creating 14-3-3-docking sites and negatively regulating MAPK activation. Biochem. J. 2012, 447, 159–166. [Google Scholar] [CrossRef] [Green Version]
  183. Shi, T.; Song, E.; Nie, S.; Rodland, K.D.; Liu, T.; Qian, W.J.; Smith, R.D. Advances in targeted proteomics and applications to biomedical research. Proteomics 2016, 16, 2160–2182. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  184. Boellner, S.; Becker, K.F. Reverse Phase Protein Arrays-Quantitative Assessment of Multiple Biomarkers in Biopsies for Clinical Use. Microarrays 2015, 4, 98. [Google Scholar] [CrossRef] [Green Version]
  185. Rahmani, A.M.; Azhir, E.; Ali, S.; Mohammadi, M.; Ahmed, O.H.; Yassin Ghafour, M.; Hasan Ahmed, S.; Hosseinzadeh, M. Artificial intelligence approaches and mechanisms for big data analytics: A systematic study. PeerJ Comput. Sci. 2021, 7, e488. [Google Scholar] [CrossRef]
  186. Karniadakis, G.E.; Kevrekidis, I.G.; Lu, L.; Perdikaris, P.; Wang, S.; Yang, L. Physics-informed machine learning. Nat. Rev. Phys. 2021, 3, 422–440. [Google Scholar] [CrossRef]
  187. Raissi, M.; Perdikaris, P.; Karniadakis, G.E. Physics-informed neural networks: A deep learning framework for solving forward and inverse problems involving nonlinear partial differential equations. J. Comput. Phys. 2019, 378, 686–707. [Google Scholar] [CrossRef]
  188. Brunton, S.L.; Proctor, J.L.; Kutz, J.N. Discovering governing equations from data by sparse identification of nonlinear dynamical systems. Proc. Natl. Acad. Sci. USA 2016, 113, 3932–3937. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  189. Udrescu, S.M.; Tegmark, M. AI Feynman: A physics-inspired method for symbolic regression. Sci. Adv. 2020, 6, eaay2631. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  190. Biswas, N.; Chakrabarti, S. Artificial Intelligence (AI)-Based Systems Biology Approaches in Multi-Omics Data Analysis of Cancer. Front. Oncol. 2020, 10, 588221. [Google Scholar] [CrossRef] [PubMed]
  191. Ibrahim, Y.M.; Hamid, A.M.; Deng, L.; Garimella, S.V.; Webb, I.K.; Baker, E.S.; Smith, R.D. New frontiers for mass spectrometry based upon structures for lossless ion manipulations. Analyst 2017, 142, 1010–1021. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  192. Jumper, J.; Evans, R.; Pritzel, A.; Green, T.; Figurnov, M.; Ronneberger, O.; Tunyasuvunakool, K.; Bates, R.; Zidek, A.; Potapenko, A.; et al. Highly accurate protein structure prediction with AlphaFold. Nature 2021, 596, 583–589. [Google Scholar] [CrossRef] [PubMed]
Figure 1. (A) Standard workflow for targeted analysis of protein phosphorylation. Selected studies have demonstrated the use of SRM to measure phosphorylation dynamics of signaling molecules complexes such as the (B) EGFR-MAPK pathway [53] and (C) PI3K-mTOR and MAPK pathway [54].
Figure 1. (A) Standard workflow for targeted analysis of protein phosphorylation. Selected studies have demonstrated the use of SRM to measure phosphorylation dynamics of signaling molecules complexes such as the (B) EGFR-MAPK pathway [53] and (C) PI3K-mTOR and MAPK pathway [54].
Molecules 28 01143 g001
Figure 2. Absolute quantifications increase the identifiability of model parameters.
Figure 2. Absolute quantifications increase the identifiability of model parameters.
Molecules 28 01143 g002
Table 1. Select studies that have demonstrated the use of targeted phosphoproteomics to quantify phosphorylation dynamics of signaling pathways.
Table 1. Select studies that have demonstrated the use of targeted phosphoproteomics to quantify phosphorylation dynamics of signaling pathways.
PathwayGoalSample TypeStarting AmountEnrichment StrategyQuantified PhosphositesReference
Breast cancerClassify high- and low-risk patient groupsBreast tumor500 µgFe3+-IMAC19[124]
PI3K-mTOR and MAPKOncogene-induced senescence in the presence of BEZ235, a dual PI3K/mTOR inhibitorHuman diploid fibrolast (HDF) cell line200 µgTi4+-IMAC89[54]
PI3K-mTOR and MAPKPhosphorylation dynamics of rapamycin mechanistic targetsPC9 and H1975 NSCLC cell line300 µgTi4+-IMAC + HpH42[125]
PI3K-mTOR and Cell CycleSensitive disruptions MAPK, PIK3/mTOR, and cell cycle signaling pathways using 26 inhibitorsMCF7, PC3, and HL60 cell lines500 µgFe3+-IMAC92[126]
DNA damage response (DDR)Understanding DDR functionality, in response to ionizing radiation and methyl methanesulfonateMCF10A cell line and human PBMCs200 µgFe3+-IMAC107[127]
DNA damage response (DDR)Understanding DDR functionality, in response ionizing radiationMCF10A cell line and human breast tumors500 µgAntibody29[128]
DNA damage response (DDR)Understanding DDR functionality, in response ionizing radiation and 4-nitroquinoline 1-oxideLCL cell line and human PBMCs500 µgAntibody25[129]
EGFR-MAPKPhosphorylation dynamics in response to EGF perturbationMCF7 and Hs578T cell lines25–100 µgFe3+-IMAC and TiO234[53]
MAPKRewiring of MAPK signaling in drug-resistant BRAFV600E melanomasA375 cell line10 and 100 µgFe3+-IMAC22[32]
Multiple cancer pathwaysCharacterize dynamic signaling across diverse cancer pathwaysPatient-derived xenografts of triple negative breast cancer and human medulloblastoma tumors500 µgAntibody and Fe3+-IMAC284[93]
Table 2. Selected applications of phosphorylation data in mathematical modeling.
Table 2. Selected applications of phosphorylation data in mathematical modeling.
PathwayData SourceExperimental ConditionsModelReadoutsLiterature
ErbB signaling pathwayWestern blotting, literature10 timepoints, 2 ligands, 2 concentrations, 4 cell linesODE3 phosphorylated proteins[162]
Multisite Tau phosphorylationWestern blotting5 timepoints, 2 ligandsODE10 phosphosites from 2 phosphorylated proteins[163]
MAPK signalingWestern blotting2 timepoints, 2 ligands, 3 gene knockdownsModular Response Analysis3 phosphorylated proteins[164]
T Cell receptor signalingGlobal Phosphoproteomics5 timepoints, 2 stimulatorsNFSim (rule-based model)22 phosphosites from 17 phosphorylated proteins[165]
MAPK signalingELISA and targeted phosphoproteomics1 ligand, 2 inhibitors at 10 concentrationsMARM1 (rule-based model)2 phosphorylated proteins[32]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Dakup, P.P.; Feng, S.; Shi, T.; Jacobs, J.M.; Wiley, H.S.; Qian, W.-J. Targeted Quantification of Protein Phosphorylation and Its Contributions towards Mathematical Modeling of Signaling Pathways. Molecules 2023, 28, 1143. https://doi.org/10.3390/molecules28031143

AMA Style

Dakup PP, Feng S, Shi T, Jacobs JM, Wiley HS, Qian W-J. Targeted Quantification of Protein Phosphorylation and Its Contributions towards Mathematical Modeling of Signaling Pathways. Molecules. 2023; 28(3):1143. https://doi.org/10.3390/molecules28031143

Chicago/Turabian Style

Dakup, Panshak P., Song Feng, Tujin Shi, Jon M. Jacobs, H. Steven Wiley, and Wei-Jun Qian. 2023. "Targeted Quantification of Protein Phosphorylation and Its Contributions towards Mathematical Modeling of Signaling Pathways" Molecules 28, no. 3: 1143. https://doi.org/10.3390/molecules28031143

APA Style

Dakup, P. P., Feng, S., Shi, T., Jacobs, J. M., Wiley, H. S., & Qian, W. -J. (2023). Targeted Quantification of Protein Phosphorylation and Its Contributions towards Mathematical Modeling of Signaling Pathways. Molecules, 28(3), 1143. https://doi.org/10.3390/molecules28031143

Article Metrics

Back to TopTop