Time-resolved mass spectrometry of tyrosine phosphorylation sites in the epidermal growth factor receptor signaling network reveals dynamic modules

Ligand binding to cell surface receptors initiates a cascade of signaling events regulated by dynamic phosphorylation events on a multitude of pathway proteins. Quantitative features, including intensity, timing, and duration of phosphorylation of particular residues, may play a role in determining cellular response, but experimental data required for analysis of these features have not previously been available. To understand the dynamic operation of signaling cascades, we have developed a method enabling the simultaneous quantification of tyrosine phosphorylation of specific residues on dozens of key proteins in a time-resolved manner, downstream of epidermal growth factor receptor (EGFR) activation. Tryptic peptides from four different EGFR stimulation time points were labeled with four isoforms of the iTRAQ reagent to enable downstream quantification. After mixing of the labeled samples, tyrosine-phosphorylated peptides were immunoprecipitated with an anti-phosphotyrosine antibody and further enriched by IMAC before LC/MS/MS analysis. Database searching and manual confirmation of peptide phosphorylation site assignments led to the identification of 78 tyrosine phosphorylation sites on 58 proteins from a single analysis. Replicate analyses of a separate biological sample provided both validation of this first data set and identification of 26 additional tyrosine phosphorylation sites and 18 additional proteins. iTRAQ fragment ion ratios provided time course phosphorylation profiles for each site. The data set of quantitative temporal phosphorylation profiles was further characterized by self-organizing maps, which resulted in identification of several cohorts of tyrosine residues exhibiting self-similar temporal phosphorylation profiles, operationally defining dynamic modules in the EGFR signaling network consistent with particular cellular processes. The presence of novel proteins and associated tyrosine phosphorylation sites within these modules indicates additional components of this network and potentially localizes the topological action of these proteins. Additional analysis and modeling of the data generated in this study are likely to yield more sophisticated models of receptor tyrosine kinase-initiated signal transduction, trafficking, and regulation.     

Construction of human activity‐based phosphorylation networks

The landscape of human phosphorylation networks has not been systematically explored, representing vast, unchartered territories within cellular signaling networks. Although a large number of in vivo phosphorylated residues have been identified by mass spectrometry (MS)‐based approaches, assigning the upstream kinases to these residues requires biochemical analysis of kinase‐substrate relationships (KSRs). Here, we developed a new strategy, called CEASAR, based on functional protein microarrays and bioinformatics to experimentally identify substrates for 289 unique kinases, resulting in 3656 high‐quality KSRs. We then generated consensus phosphorylation motifs for each of the kinases and integrated this information, along with information about in vivo phosphorylation sites determined by MS, to construct a high‐resolution map of phosphorylation networks that connects 230 kinases to 2591 in vivo phosphorylation sites in 652 substrates. The value of this data set is demonstrated through the discovery of a new role for PKA downstream of Btk (Bruton's tyrosine kinase) during B‐cell receptor signaling. Overall, these studies provide global insights into kinase‐mediated signaling pathways and promise to advance our understanding of cellular signaling processes in humans.     

RAFTK/Pyk2 mediates LPA-induced PC12 cell migration

The phospholipid lysophosphatidic acid (LPA) is a normal constituent of serum that functions as a lipid growth factor and intracellular signaling molecule. In this report, we have investigated the signaling mechanism and function of the tyrosine kinase RAFTK/Pyk2 in LPA-induced cell migration. Analysis of tyrosine phosphorylation upon LPA stimulation in neuroendocrine PC12 cells revealed 6 major tyrosine-phosphorylated proteins with estimated sizes of 180, 120, 115, 68, 44, and 42 kDa. These proteins were identified as epidermal growth factor receptor (EGFR), focal adhesion kinase, RAFTK/Pyk2, paxillin, Erk 1, and Erk 2, respectively. Using specific pharmacological inhibitors, we found that the tyrosine phosphorylation of RAFTK/Pyk2 was intracellular Ca2+-dependent, but not EGFR-dependent, during LPA stimulation of these cells. Moreover, the cytoskeletal and signal scaffolding protein, paxillin, associated with and was regulated by RAFTK/Pyk2 in a Ca2+-dependent manner. Characterization of LPA receptors showed that LPA1 (Edg2) and LPA2 (Edg4) are major receptors for LPA, while LPA3 receptor (Edg7) expression was limited. Upon using the LPA1/LPA3 receptor-specific antagonist VPC 32179, we observed that inhibition of the LPA1/LPA3 receptors had no effect on the LPA-induced phosphorylation of RAFTK, strongly suggesting that the LPA2 receptor is a key mediator of RAFTK phosphorylation. Furthermore, LPA induced PC12 cell migration, which was subsequently blocked by the dominant-negative form of FAK, FRNK. Expression of a dominant-negative form of the small GTPase Ras also blocked LPA-induced cell migration and RAFTK phosphorylation. Taken together, these results indicate that RAFTK is a key signaling molecule that mediates LPA-induced PC12 cell migration in a Ras-dependent manner.     

Phosphorylation Time‐Course Study of the Response during Adenovirus Type 2 Infection

PTMs such as phosphorylations are usually involved in signal transduction pathways. To investigate the temporal dynamics of phosphoproteome changes upon viral infection, a model system of IMR‐90 cells infected with human adenovirus type 2 (Ad2) is used in a time‐course quantitative analysis combining titanium dioxide (TiO₂) particle enrichment and SILAC‐MS. Quantitative data from 1552 phosphorylated sites clustered the highly altered phosphorylated sites to the signaling by rho family GTPases, the actin cytoskeleton signaling, and the cAMP‐dependent protein kinase A signaling pathways. Their activation is especially pronounced at early time post‐infection. Changes of several phosphorylated sites involved in the glycolysis pathway, related to the activation of the Warburg effect, point at virus‐induced energy production. For Ad2 proteins, 32 novel phosphorylation sites are identified and as many as 52 phosphorylated sites on 17 different Ad2 proteins are quantified, most of them at late time post‐infection. Kinase predictions highlighted activation of PKA, CDK1/2, MAPK, and CKII. Overlaps of kinase motif sequences for viral and human proteins are observed, stressing the importance of phosphorylation during Ad2 infection.     

Characterization of epidermal growth factor receptor function in lysophosphatidic acid signaling in PC12 cells

Lysophosphatidic acid (LPA) is a lipid metabolite that induces the activation of mitogen-activated protein kinase (MAPK) through binding to the G protein-coupled receptor in a number of cell lines and cultures. Recent studies have revealed that LPA is able to rapidly induce the phosphorylation of MAPK through an epidermal growth factor (EGF) receptor-dependent pathway. We investigated the role of the EGF receptor in the signaling pathway initiated by LPA stimulation in nerve growth factor (NGF)-responsive PC12 cells well known to transiently retract their own neurites upon LPA stimulation. LPA-stimulated MAPK signaling was suppressed by the selective EGF receptor inhibitor and in the dominant negative mutant EGF receptor cell line. As in the EGF signaling pathway, the complex of EGF receptor with adapter proteins Shc and Sos was formed in response to LPA stimulation, suggesting there is an intracellular mechanism for transactivation. A neurite retraction assay was also performed to examine the role of the EGF receptor in PC12 cell differentiation, which related to the involvement of LPA-induced neurite retraction. These results suggest that the receptor tyrosine kinase can be activated in a ligand-independent manner through intracellular crosstalk between the signaling pathways.     

Proteomic,functional and motif‐based analysis of C‐terminal Src kinase‐interacting proteins

C‐terminal Src kinase (Csk) that functions as an essential negative regulator of Src family tyrosine kinases (SFKs) interacts with tyrosine‐phosphorylated molecules through its Src homology 2 (SH2) domain, allowing it targeting to the sites of SFKs and concomitantly enhancing its kinase activity. Identification of additional Csk‐interacting proteins is expected to reveal potential signaling targets and previously undescribed functions of Csk. In this study, using a direct proteomic approach, we identified 151 novel potential Csk‐binding partners, which are associated with a wide range of biological functions. Bioinformatics analysis showed that the majority of identified proteins contain one or several Csk‐SH2 domain‐binding motifs, indicating a potentially direct interaction with Csk. The interactions of Csk with four proteins (partitioning defective 3 (Par3), DDR1, SYK and protein kinase C iota) were confirmed using biochemical approaches and phosphotyrosine 1127 of Par3 C‐terminus was proved to directly bind to Csk‐SH2 domain, which was consistent with predictions from in silico analysis. Finally, immunofluorescence experiments revealed co‐localization of Csk with Par3 in tight junction (TJ) in a tyrosine phosphorylation‐dependent manner and overexpression of Csk, but not its SH2‐domain mutant lacking binding to phosphotyrosine, promoted the TJ assembly in Madin‐Darby canine kidney cells, implying the involvement of Csk‐SH2 domain in regulating cellular TJs. In conclusion, the newly identified potential interacting partners of Csk provided new insights into its functional diversity in regulation of numerous cellular events, in addition to controlling the SFK activity.     

Quantitative phosphoproteome analysis of lysophosphatidic acid induced chemotaxis applying dual-step (18)O labeling coupled with immobilized metal-ion affinity chromatography

Reversible protein phosphorylation is a central cellular regulatory mechanism in modulating protein activity and propagating signals within cellular pathways and networks. Development of more effective methods for the simultaneous identification of phosphorylation sites and quantification of temporal changes in protein phosphorylation could provide important insights into molecular signaling mechanisms in various cellular processes. Here we present an integrated quantitative phosphoproteomics approach and its application for comparative analysis of Cos-7 cells in response to lysophosphatidic acid (LPA) gradient stimulation. The approach combines trypsin-catalyzed (16)O/ (18)O labeling plus (16)O/ (18)O-methanol esterification for quantitation, a macro-immobilized metal-ion affinity chromatography trap for phosphopeptide enrichment, and LC-MS/MS analysis. LC separation and MS/MS are followed by neutral loss-dependent MS/MS/MS for phosphopeptide identification using a linear ion trap (LTQ)-FT mass spectrometer. A variety of phosphorylated proteins were identified and quantified including receptors, kinases, proteins associated with small GTPases, and cytoskeleton proteins. A number of hypothetical proteins were also identified as differentially expressed followed by LPA stimulation, and we have shown evidence of pseudopodia subcellular localization of one of these candidate proteins. These results demonstrate the efficiency of this quantitative phosphoproteomics approach and its application for rapid discovery of phosphorylation events associated with LPA gradient sensing and cell chemotaxis.     

Dynamic profiling of the post-translational modifications and interaction partners of epidermal growth factor receptor signaling after stimulation by epidermal growth factor using Extended Range Proteomic Analysis (ERPA)

In a recent report, we introduced Extended Range Proteomic Analysis (ERPA), an intermediate approach between top-down and bottom-up proteomics, for the comprehensive characterization at the trace level (fmol level) of large and complex proteins. In this study, we extended ERPA to determine quantitatively the temporal changes that occur in the tyrosine kinase receptor, epidermal growth factor receptor (EGFR), upon stimulation. Specifically A 431 cells were stimulated with epidermal growth factor after which EGFR was immunoprecipitated at stimulation times of 0, 0.5, 2, and 10 min as well as 4 h. High sequence coverage was obtained (96%), and methods were developed for label-free quantitation of phosphorylation and glycosylation. A total of 13 phosphorylation sites were identified, and the estimated stoichiometry was determined over the stimulation time points, including Thr(P) and Ser(P) sites in addition to Tyr(P) sites. A total of 10 extracellular domain N-glycan sites were also identified, and major glycoforms at each site were quantitated. No change in the extent of glycosylation with stimulation was observed as expected. Finally potential binding partners to EGFR were identified based on changes in the amount of protein pulled down with EGFR as a function of time of stimulation. Many of the 19 proteins identified are known binding partners of EGFR. This work demonstrates that comprehensive characterization provides a powerful tool to aid in the study of important therapeutic targets. The detailed molecular information will prove useful in future studies in tissue.     

Temporal analysis of phosphotyrosine-dependent signaling networks by quantitative proteomics

To study the global dynamics of phosphotyrosine-based signaling events in early growth factor stimulation, we developed a mass spectrometric method that converts temporal changes to differences in peptide isotopic abundance. The proteomes of three cell populations were metabolically encoded with different stable isotopic forms of arginine. Each population was stimulated by epidermal growth factor for a different length of time, and tyrosine-phosphorylated proteins and closely associated binders were affinity purified. Arginine-containing peptides occurred in three forms, which were quantified; we then combined two experiments to generate five-point dynamic profiles. We identified 81 signaling proteins, including virtually all known epidermal growth factor receptor substrates, 31 novel effectors and the time course of their activation upon epidermal growth factor stimulation. Global activation profiles provide an informative perspective on cell signaling and will be crucial to modeling signaling networks in a systems biology approach.     

Convergence of multiple signaling cascades at glycogen synthase kinase 3: Edg receptor-mediated phosphorylation and inactivation by lysophosphatidic acid through a protein kinase C-dependent intracellular pathway

Lysophosphatidic acid (LPA) is a natural phospholipid with multiple biological functions. We show here that LPA induces phosphorylation and inactivation of glycogen synthase kinase 3 (GSK-3), a multifunctional serine/threonine kinase. The effect of LPA can be reconstituted by expression of Edg-4 or Edg-7 in cells lacking LPA responses. Compared to insulin, LPA stimulates only modest phosphatidylinositol 3-kinase (PI3K)-dependent activation of protein kinase B (PKB/Akt) that does not correlate with the magnitude of GSK-3 phosphorylation induced by LPA. PI3K inhibitors block insulin- but not LPA-induced GSK-3 phosphorylation. In contrast, the effect of LPA, but not that of insulin or platelet-derived growth factor (PDGF), is sensitive to protein kinase C (PKC) inhibitors. Downregulation of endogenous PKC activity selectively reduces LPA-mediated GSK-3 phosphorylation. Furthermore, several PKC isotypes phosphorylate GSK-3 in vitro and in vivo. To confirm a specific role for PKC in regulation of GSK-3, we further studied signaling properties of PDGF receptor beta subunit (PDGFRbeta) in HEK293 cells lacking endogenous PDGF receptors. In clones expressing a PDGFRbeta mutant wherein the residues that couple to PI3K and other signaling functions are mutated with the link to phospholipase Cgamma (PLCgamma) left intact, PDGF is fully capable of stimulating GSK-3 phosphorylation. The process is sensitive to PKC inhibitors in contrast to the response through the wild-type PDGFRbeta. Therefore, growth factors, such as PDGF, which control GSK-3 mainly through the PI3K-PKB/Akt module, possess the ability to regulate GSK-3 through an alternative, redundant PLCgamma-PKC pathway. LPA and potentially other natural ligands primarily utilize a PKC-dependent pathway to modulate GSK-3.     

Characterization of early autophagy signaling by quantitative phosphoproteomics

Under conditions of nutrient shortage autophagy is the primary cellular mechanism ensuring availability of substrates for continuous biosynthesis. Subjecting cells to starvation or rapamycin efficiently induces autophagy by inhibiting the MTOR signaling pathway triggering increased autophagic flux. To elucidate the regulation of early signaling events upon autophagy induction, we applied quantitative phosphoproteomics characterizing the temporal phosphorylation dynamics after starvation and rapamycin treatment. We obtained a comprehensive atlas of phosphorylation kinetics within the first 30 min upon induction of autophagy with both treatments affecting widely different cellular processes. The identification of dynamic phosphorylation already after 2 min demonstrates that the earliest events in autophagy signaling occur rapidly after induction. The data was subjected to extensive bioinformatics analysis revealing regulated phosphorylation sites on proteins involved in a wide range of cellular processes and an impact of the treatments on the kinome. To approach the potential function of the identified phosphorylation sites we performed a screen for MAP1LC3-interacting proteins and identified a group of binding partners exhibiting dynamic phosphorylation patterns. The data presented here provide a valuable resource on phosphorylation events underlying early autophagy induction.     

Signal characteristics of G protein-transactivated EGF receptor.

The epidermal growth factor receptor (EGFR) tyrosine kinase recently was identified as providing a link to mitogen-activated protein kinase (MAPK) in response to G protein-coupled receptor (GPCR) agonists in Rat-1 fibroblasts. This cross-talk pathway is also established in other cell types such as HaCaT keratinocytes, primary mouse astrocytes and COS-7 cells. Transient expression of either Gq- or Gi-coupled receptors in COS-7 cells allowed GPCR agonist-induced EGFR transactivation, and lysophosphatidic acid (LPA)-generated signals involved the docking protein Gab1. The increase in SHC tyrosine phosphorylation and MAPK stimulation through both Gq- and Gi-coupled receptors was reduced strongly upon selective inhibition of EGFR function. Inhibition of phosphoinositide 3-kinase did not affect GPCR-induced stimulation of EGFR tyrosine phosphorylation, but inhibited MAPK stimulation, upon treatment with both GPCR agonists and low doses of EGF. Furthermore, the Src tyrosine kinase inhibitor PP1 strongly interfered with LPA- and EGF-induced tyrosine phosphorylation and MAPK activation downstream of EGFR. Our results demonstrate an essential role for EGFR function in signaling through both Gq- and Gi-coupled receptors and provide novel insights into signal transmission downstream of EGFR for efficient activation of the Ras/MAPK pathway.     

Ligand regulates epidermal growth factor receptor kinase specificity: activation increases preference for GAB1 and SHC versus autophosphorylation sites

The epidermal growth factor receptor (EGFR) kinase catalyzes phosphorylation of tyrosines in its C terminus and in other cellular targets upon epidermal growth factor (EGF) stimulation. Here, by using peptides derived from EGFR autophosphorylation sites and cellular substrates, we tested the hypothesis that ligand may function to regulate EGFR kinase specificity by modulating the binding affinity of peptide sequences to the active site. Measurement of the steady-state kinetic parameters, K(m) and k(cat), revealed that EGF did not affect the binding of EGFR peptides but increased the binding affinity for peptides corresponding to the major EGFR-mediated phosphorylation sites of the adaptor proteins Gab1 (Tyr-627) and Shc (Tyr-317), and for peptides containing the previously identified optimal EGFR kinase substrate sequence EEEEYFELV (3-7-fold). Conversely, EGF stimulation increased k(cat) approximately 5-fold for all peptides. Thus, ligand changed the relative preference of the EGFR kinase for substrates as evidenced by EGF increases of approximately 5-fold in the specificity constants (k(cat)/K(m)) for EGFR peptides, whereas approximately 15-40-fold increases were observed for other peptides, such as Gab1 Tyr-627. Furthermore, we demonstrate that EGF (i) increased the binding affinity of EGFR to Gab1 Tyr-627 and Shc Tyr-317 sites in purified GST fusion proteins approximately 4-6-fold, and (ii) EGF significantly enhanced the phosphorylation of these sites, relative to EGFR autophosphorylation, in cell lysates containing the full-length Gab1 and Shc proteins. Analysis of peptides containing amino acid substitutions indicated that residues C-terminal to the target tyrosine were critical for EGF-stimulated increases in substrate binding and regulation of kinase specificity. To our knowledge, this represents the first demonstration that ligand can alter specificity of a receptor kinase toward physiologically relevant targets.     

Platelet-derived growth factor and lysophosphatidic acid inhibit growth hormone binding and signaling via a protein kinase C-dependent pathway

Growth hormone (GH) regulates body growth and metabolism. GH exerts its biological action by stimulating JAK2, a GH receptor (GHR)-associated tyrosine kinase. Activated JAK2 phosphorylates itself and GHR, thus initiating multiple signaling pathways. In this work, we demonstrate that platelet-derived growth factor (PDGF) and lysophosphatidic acid (LPA) down-regulate GH signaling via a protein kinase C (PKC)-dependent pathway. PDGF substantially reduces tyrosyl phosphorylation of JAK2 induced by GH but not interferon-gamma or leukemia inhibitory factor. PDGF, but not epidermal growth factor, decreases tyrosyl phosphorylation of GHR (by approximately 90%) and the amount of both total cellular GHR (by approximately 80%) and GH binding (by approximately 70%). The inhibitory effect of PDGF on GH-induced tyrosyl phosphorylation of JAK2 and GHR is abolished by depletion of 4beta-phorbol 12-myristate 13-acetate (PMA)-sensitive PKCs with chronic PMA treatment and is severely inhibited by GF109203X, an inhibitor of PKCs. In contrast, extracellular signal-regulated kinases 1 and 2 and phosphatidylinositol 3-kinase appear not to be involved in this inhibitory effect of PDGF. LPA, a known activator of PKC, also inhibits GH-induced tyrosyl phosphorylation of JAK2 and GHR and reduces the number of GHR. We propose that ligands that activate PKC, including PDGF, LPA, and PMA, down-regulate GH signaling by decreasing the number of cell surface GHR through promoting GHR internalization and degradation and/or cleavage of membrane GHR and release of the extracellular domain of GHR.     

Identification of targets of c-Src tyrosine kinase by chemical complementation and phosphoproteomics

The cellular proto-oncogene c-Src is a nonreceptor tyrosine kinase involved in cell growth and cytoskeletal regulation. Despite being dysregulated in a variety of human cancers, its precise functions are not fully understood. Identification of the substrates of c-Src remains a major challenge, because there is no simple way to directly stimulate its activity. Here we combine the chemical rescue of mutant c-Src and global quantitative phosphoproteomics to obtain the first high resolution snapshot of the range of tyrosine phosphorylation events that occur in the cell immediately after specific c-Src stimulation. After enrichment by anti-phosphotyrosine antibodies, we identified 29 potential novel c-Src substrate proteins. Tyrosine phosphopeptide mapping allowed the identification of 382 nonredundant tyrosine phosphopeptides on 213 phosphoproteins. Stable isotope labeling of amino acids in cell culture-based quantitation allowed the detection of 97 nonredundant tyrosine phosphopeptides whose level of phosphorylation is increased by c-Src. A large number of previously uncharacterized c-Src putative protein targets and phosphorylation sites are presented here, a majority of which play key roles in signaling and cytoskeletal networks, particularly in cell adhesion. Integrin signaling and focal adhesion kinase signaling pathway are two of the most altered pathways upon c-Src activation through chemical rescue. In this context, our study revealed the temporal connection between c-Src activation and the GTPase Rap1, known to stimulate integrin-dependent adhesion. Chemical rescue of c-Src provided a tool to dissect the spatiotemporal mechanism of activation of the Rap1 guanine exchange factor, C3G, one of the identified potential c-Src substrates that plays a role in focal adhesion signaling. In addition to unveiling the role of c-Src in the cell and, specifically, in the Crk-C3G-Rap1 pathway, these results exemplify a strategy for obtaining a comprehensive understanding of the functions of nonreceptor tyrosine kinases with high specificity and kinetic resolution.     

Temporal proteomic analysis of IGF‐1R signalling in MCF‐7 breast adenocarcinoma cells

Dysregulation of the insulin‐like growth factor 1 receptor signalling network is implicated in tumour growth and resistance to chemotherapy. We explored proteomic changes resulting from insulin‐like growth factor 1 stimulation of MCF‐7 adenocarcinoma cells as a function of time. Quantitative analysis using iTRAQ? reagents and 2‐D LC‐MS/MS analysis of three biological replicates resulted in the identification of 899 proteins (p≤0.05) with an estimated mean false‐positive rate of 2.6%. Quantitative protein expression was obtained from 681 proteins. Further analysis by supervised k‐means clustering identified five temporal clusters, which were submitted to the FuncAssociate server to assign overrepresented gene ontology terms. Proteins associated with vesicle transport were significantly overrepresented. We further analyzed our data set for proteins showing temporal significance using the software, extraction and analysis of differential gene expression, resulting in 20 significantly and temporally changing proteins (p≤0.1). These significant proteins play roles in, among others, altered glucose metabolism (lactate dehydrogenase A and pyruvate kinase M1/M2) and cellular stress (nascent polypeptide‐associated complex subunit α and heat shock (HSC70) proteins). We used multiple reaction monitoring to validate these interesting proteins and have revealed several differences in relative peptide expression corresponding to protein isoforms and variants.