|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
,||
,¶
,||
,||,**
From the
Ottawa Health Research Institute, Molecular Medicine Program, Ottawa Hospital, Ottawa, Ontario K1H 8L6, Canada,
Institute for Biological Sciences, National Research Council, Ottawa, Ontario K1A 0R6, Canada, and || Department of Cellular and Molecular Medicine and Centre for Neuromuscular Disease, Faculty of Medicine, University of Ottawa, Ottawa, Ontario K1H 8M5, Canada
| ABSTRACT |
|---|
|
|
|---|
The best characterized kinase signaling factors in myogenesis are the p38 mitogen-activated protein kinases. Inhibitors of p38 signaling block muscle-specific gene expression and myotube formation (3, 4). There is evidence that p38 acts through multiple mechanisms including phosphorylation of the MRF cofactor myocyte enhancer factor 2 (46), phosphorylation of the MyoD cofactor E47, and enhanced recruitment of the SWI/SNF chromatin-remodeling complex to MRF-targeted promoters (7). Another key signaling kinase in muscle is AKT1. Stimulation of AKT1 by IGF-1 or insulin induces muscle hypertrophy, whereas inhibition of AKT1 by glucocorticoids promotes muscle atrophy. Hypertrophy is associated with AKT-mediated mTOR (mammalian target of rapamycin) activation and GSK3B repression (8). Conversely AKT1 inhibition leads to dephosphorylation and activation of FOXO transcription factors resulting in atrophy (9). Although some kinase activities are required for normal muscle development, others must be suppressed. JNK1 is normally inactive during myogenesis, and its activation leads to muscle pathology (10). In addition to these kinases numerous other kinase proteins have been implicated as regulators of the differentiation process. For example CDK5 expression and activity increases during early myogenesis, and expression of dominant-negative forms of CDK5 inhibits muscle formation (11, 12). However, an important unresolved question is how these different pathways are integrated into an overall network of kinase-substrate interactions to control myogenic differentiation.
Previously we profiled changes to protein phosphorylation during myogenesis on a proteome-wide scale by using phosphoprotein enrichment and comparative two-dimensional gel electrophoresis (13). Here we extend those findings by mapping sequence-specific sites of protein phosphorylation in growing and differentiating C2C12 cells. To place these results in a systems context, we took advantage of recent advances in the availability of protein-protein interaction databases. Databases of protein-protein interactions have been constructed from yeast, fruit fly, worm, and human cells using yeast two-hybrid, tag/pulldown, and literature search approaches (14). Subsets of the total protein-protein interaction set such as kinase-substrate interaction maps have also been constructed (15). To date, however, relatively few studies have leveraged this type of information to assist in the interpretation of results. Here we used kinase-substrate interaction databases to bioinformatically reconstruct a kinase signaling network based upon our experimentally identified phosphorylation events. This approach yielded a model kinase-substrate network that predicted several known features of myogenic signaling and revealed the potential relevance of newly identified phosphorylation events in relation to known muscle-related biochemical pathways.
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
Phosphoprotein Enrichment
Phosphoprotein enrichment was performed by affinity column purification of total myoblast proteins using PhosphoProtein Purification kits (Qiagen) as described previously (13).
Analysis of Putative Phosphoproteins by One-dimensional Gel Electrophoresis
Phosphoprotein-enriched samples (typically 100150 µg) were washed in water and concentrated using Amicon 10,000 molecular weight cutoff ultrafiltration columns (Millipore). Aliquots of retentate were separated by SDS-PAGE (12.5% acrylamide "mini" format). Gels were stained with Bio-Safe colloidal G-250 Coomassie Blue (Bio-Rad). 23 1-mm gel sections were excised from each lane, and in-gel digest and extraction of proteins was performed.
In-gel Digest
Gel pieces cut into 1-mm cubes were destained in 100 mM ammonium bicarbonate (NH4HCO3) in 30% ACN. Destained gel pieces were washed in dH2O and shrunk in ACN. Gel pieces were incubated for 1 h at 56 °C in 10 mM DTT, 50 mM NH4HCO3; washed in 50 mM NH4HCO3; and incubated for 1 h at room temperature in 55 mM iodoacetamide, 50 mM NH4HCO3. Gel pieces were shrunk with ACN and reswollen in 50 mM NH4HCO3 with Promega modified trypsin (10 ng/µl). Sufficient 50 mM NH4HCO3 was added to cover the gel pieces, and sealed tubes were incubated overnight at 37 °C.
Preparation of Samples for IMAC
Approximately 15 µg of phosphoprotein-enriched protein extract was suspended in 50 mM NH4HCO3 containing Promega modified trypsin (10 ng/µl). After overnight incubation at 37 °C digest solutions were concentrated to
40 µl and desalted through 20-µl GELoader tips (Eppendorf) packed with 20 µl of OligoR3 resin (Applied Biosystems) as described previously (16). Peptides were washed once by 20 µl of 1% acetic acid and eluted by 20 µl of 50% ACN, 1% acetic acid. Elution solvent was removed by SpeedVac and replaced by 40 µl of the appropriate IMAC binding buffer (see below).
Preparation of Metal-chelating Resin for IMAC
A 50% slurry of POROS 20MC resin (Applied Biosystems) or Propac-IMAC resin (Dionex) was prepared in dH2O. Resin was washed twice in water, stripped by 50 mM EDTA, washed in dH2O, and then washed in 1% acetic acid. Resin was incubated in 100 mM iron chloride, 1% acetic acid for 15 min with agitation. Activated resin was rinsed in 30% ACN, 1% acetic acid and then in 1% acetic acid.
IMAC Enrichment Using GELoader Tips
30 µl of prepared IMAC resin was added to the peptide sample in binding buffer (100 mM NaCl in 1% acetic acid) and agitated for 15 min before loading into a GELoader tip. Resin was rinsed once by 20 µl of a solution of 100 mM NaCl in 30% acetonitrile and 1% acetic acid and eluted in two steps by 250 mM ammonium phosphate, pH 9. In the first step, 5 µl was pushed into the resin and left for 5 min to neutralize the remaining acetic acid. In the second step, 15 µl of the elution buffer was pushed through the resin. The 20 µl of phosphopeptide-carrying eluate was evaporated by SpeedVac, and the phosphopeptides were suspended in 20 µl of 5% acetonitrile in 0.1% acetic acid.
IMAC Enrichment Using Pierce Phosphopeptide Isolation Kit
Gallium-chelated columns were used according to the manufacturers recommendation with the following modifications. Using 40 µl of a solution of 5% acetic acid as binding buffer, the tryptic peptides were transferred to a Pierce column, left for 15 min at room temperature, and then centrifuged for 1 min at 3000 x g. The Gallium-chelated disc was rinsed twice with 50 µl of 100 mM NaCl in 1% acetic acid, twice with 30% acetonitrile in 1% acetic acid, and once with dH2O and then eluted by three additions of 20 µl of 250 mM ammonium phosphate, pH 9. Eluates were pooled and evaporated by SpeedVac, and phosphopeptides were suspended in 20 µl of 5% ACN, 0.1% acetic acid before MS analysis.
Q-TOF Analyses
A fraction of each tryptic digest (typically 510 µl) was analyzed by nano-LC-MS/MS using a Q-TOF Ultima hybrid quadrupole time-of-flight mass spectrometer coupled to a CapLC capillary HPLC system (Waters, Milford, MA). Digests were separated on a 75-µm-inner diameter x 150-mm Inertsil ODS 3, 5-µm nano-HPLC column (Dionex/LC Packings, Sunnyvale, CA) using the following gradient conditions: 560% ACN, 0.2% formic acid in 30 min, and 6090% ACN in 5 min. The mass spectrometer was set to automatically acquire MS/MS spectra on doubly, triply, and quadruply charged ions. The MS survey scan range was 4001600 m/z; the MS/MS range was 502000. After the first LC-MS/MS run, an exclusion file was created listing all previously analyzed ions that appeared to represent non-phosphorylated peptides, and a second fraction of each tryptic digest was injected with the mass spectrometer set to ignore the ions selected during the first run.
Data Processing
Database searching was carried out using Mascot DaemonTM (Matrix Science) against the National Center for Biotechnology Information (NCBI) genome sequence database. In addition to carbamidomethylcysteines and oxidized methionine, Mascot was asked to search for the potential presence of phosphate groups on serine, threonine, or tyrosine residues by looking in the MS/MS fragmentation spectra for the mass of a HPO3 group on these residues or for ions resulting from a neutral loss of HPO3 or H3PO4 (79.9663 and 97.9769 Da, respectively). All putative phosphopeptide MS/MS spectra identified by Mascot were manually verified for both the presence of a phosphate group and the peptide sequence. In many cases, the MS/MS spectrum allowed precise identification of the modified residue.
Identification of Phosphopeptides in Existing LC-MS/MS Data
Previously phosphoprotein-enriched samples were separated by 2D-GE, and 190 proteins were identified by LC-MS/MS (13). Many of these proteins showed evidence of differential phosphorylation between the undifferentiated and differentiating states (13). These data was extensively reanalyzed against the latest available NCBInr (non-redundant) database using Mascot server software version 2.0.04 (Matrix Science) allowing up to two missed cleavages and the variable modifications protein N-terminal acetylation, carbamidomethylation, deamidation (NQ), oxidation (Met), pyro-Glu (N-terminal Gln), and phosphorylation (YST). All putative phosphopeptide MS/MS spectra were manually reviewed.
Reconstructing a Pathway Map from Phosphopeptides
For each experimentally determined phosphorylation site, any kinase(s) known to phosphorylate the site in vivo or in vitro was identified using the bioinformatics tools KinaSource (A. Knebel, Dundee, Scotland, UK) and PhosphoSite (Cell Signaling Technology, Danvers, MA). In a small number of cases where indicated, the relevant kinase was predicted using ScanSite (17). Then KinaSource was queried for other known substrates of each kinase. The process of mapping substrates to kinases was then repeated for all proteins in the model, including the kinases themselves. Database-derived results were verified and extended by examination of the primary research literature. The identified substrate-kinase interactions were then diagrammed using the bioinformatics tool HubView (18).
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
|
|
|
DYRK
DYRK interactions were prominent in the model network. Induction of DYRK family kinases is required for muscle differentiation (19) at least in part via DYRK-mediated phosphorylation and regulation of histone deacetylases (20). Here we found that the putative DYRK substrates CRHSP24 and CRMP4 are also phosphorylated in C2C12 (Table I and Fig. 3). Bioinformatics modeling revealed five more DYRK-substrate interactions (p27kip, p21waf, FOXO1A, glycogen synthase, and STAT3), suggesting additional mechanisms by which DYRK may contribute to myogenesis. Two of the phosphorylation sites ascribed to DYRK2 activity are thought to be priming sites for GSK3. In CRMP4 (Dpysl3, Drp3) Ser-522 appears to represent the priming site for a series of GSK3-mediated phosphorylations at serines 518, 514, and 509 (21). In C2C12 cells both site-specific (Table I) and general phosphorylation (13) of CRMP4 was detected only in differentiating cells, whereas Ser-522 is constitutively phosphorylated in brain. Although CRMP4 has been studied almost exclusively in the context of neurons, these data and the observation that the CRMP4 promoter contains a MyoD/myogenin binding site (22) indicate that CRMP proteins likely play an important role in muscle differentiation. eIF2B (eukaryotic translation initiation factor 2B) is also a substrate of both DYRK and GSK3. DYRK-mediated phosphorylation at Ser-539 is thought to permit subsequent phosphorylation of Ser-535 by GSK3 leading to repression of translation (23). Because muscle hypertrophy is associated with GSK3 inhibition and eIF2 derepression, this suggests that promyogenic functions of DYRK could be blocked by GSK3 activity. These observations are consistent with the model that DYRK activity is increased and GSK3 is repressed during myogenesis. In agreement with this principle, several putative GSK3 substrate proteins that were not DYRK substrates, pyruvate dehydrogenase and CRMP2, showed evidence of phosphorylation only in non-differentiated myoblasts, although specific phosphorylation sites were not determined (13).
|
PEA-15
Our previous observation of PEA-15 phosphorylation in C2C12 (13) was confirmed and extended with evidence of a second phosphorylation site (Table I). Ser-116 phosphorylated PEA-15 was repeatedly observed in growing cells but never in differentiating cells. This observation is consistent with the recent report that the Ser-104/Ser-116 dephosphorylated form of PEA-15 promotes cell cycle arrest via sequestration of ERK (27), and indeed, reduced ERK activity is a prerequisite for myoblast differentiation (28).
PKC-
A potential role for PKC-
was inferred in the model (Fig. 2C), and the potential PKC-
substrates 14-3-3, HMGA1 (high mobility group AT-hook 1), and HNRNPK were detected in phosphoprotein enrichment experiments (13). Although site-specific evidence of PKC-
-mediated phosphorylation is lacking, these observations are interesting in light of a recent report that constitutively active PKC-
promotes whereas kinase-inactive PKC-
inhibits IGF-1-mediated proliferation of myoblasts (29). PKC-
, as well as DYRK1a and ERK1/ERK2, can also phosphorylate STAT3, which is reported to promote myoblast proliferation and inhibit MyoD function and muscle differentiation (3032). Because the aforementioned substrates contain multiple phosphorylation sites, additional studies will be required to determine the precise role of PKC.
Limitations to Application and Interpretation
Because the experimental data cover only a very small proportion of the phosphoproteome, it is important to consider the validity of extrapolating a model kinase network under these circumstances. The strategy used relies on three assumptions: that substrates have a limited number of kinases, that kinases have a large number of substrates, and that the overall topology of the kinase-substrate network is scale-free. First, if a substrate has many potential kinases then inference of specific kinase activities is not possible. Second, if a kinase has few substrates then the odds that one of those substrates will be found experimentally is extremely low. The last assumption determines whether the model is likely to encompass major features of the true in vivo network or if it reveals only local details. In a non-scale-free network a small data sample will show only local features, much as mapping a few roads in a city would reveal nothing about the overall layout of the traffic flow. However in a scale-free network, any two points on the network are connected by a very small number of interactions, and the overall picture emerges quickly, just as a sampling of commercial airline flights would quickly reveal the existence of traffic hubs at the major airports. Work on kinase-substrate interactions in yeast supports these three assumptions. An exhaustive evaluation of 1325 phosphoproteins and 87 kinases by Ptacek et al. (15) found that 73% of substrates were associated with only one or two kinases but that most kinases had multiple substrates. Analysis of the kinase-substrate subset of the extensive protein-protein interaction data available for yeast confirmed that this network is scale-free (18).
Another consideration was protein abundance. If myogenesis resulted in altered protein phosphorylation in only a subset of very low abundance proteins against a large backdrop of abundant invariant phosphoproteins then no relevant results would have been obtained. However, we found previously that the subset of the phosphoproteome that is sufficiently abundant for silver stain detection on gels does in fact undergo extensive alteration during early myogenesis (13). The use of data obtained from two-dimensional gels was especially valuable in retaining a biologically relevant focus. Comparative analysis of gels prepared from growing versus differentiating myocytes allowed us to focus the MS analysis on proteins that exhibited evidence of altered post-translational modification during myogenesis. All of these factors contributed to the ability to produce a model network that exhibited relevant features of interest. However, a number of caveats remain. First, because proteome coverage is incomplete, lack of detection of phosphorylation cannot be taken as strong evidence of dephosphorylation. Second, the results are biased toward the detection of relatively abundant peptides. Third, kinase-substrate databases are far from complete, although the amount of protein-protein interaction data available is increasing rapidly. We anticipate that increased coverage should not invalidate or contradict the results presented here. Because kinase networks appear to be scale-free (18), increased data coverage would likely increase the number of substrates more rapidly than the number of kinases, limiting any changes in the overall topology of the network. However, this remains to be experimentally demonstrated. Another challenge remaining to be addressed is the role of regulated phosphatases in the network. At this time inclusion of a comprehensive phosphatase dataset is not possible given the limited experimental observations in the literature combined with the absence of robust screening tools for phosphatase activity or dephosphorylation events.
Conclusion
Experimental phosphopeptide data and bioinformatics databases of kinase-substrate interactions were used to produce a model signaling network for C2C12 myogenesis. From a relatively small number of phosphopeptides, a complex network was readily constructed that inferred the existence of several known features of myogenic signaling. In addition, places for newly described phosphorylations in signaling cascades were predicted. New protein-protein interaction datasets including kinase-substrate datasets are being produced at an increasing rate, and the utility of this approach is certain to increase as more data become available.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
Published, MCP Papers in Press, September 13, 2006, DOI 10.1074/mcp.M600134-MCP200
1 The abbreviations used are: MRF, myogenic regulatory factor; 2D-GE, two-dimensional gel electrophoresis; CDK, cyclin-dependent kinase; CRHSP, calcium-regulated heat stable protein; CRMP, collapsin response mediator protein; dH2O, deionized water; DYRK, dual specificity tyrosine (Y) phosphorylation-regulated kinase; ERK, extracellular signal-regulated kinase; FOXO, forkhead box (transcription factor); GSK, glycogen synthase kinase; IGF, insulin-like growth factor; JNK, c-Jun NH2-terminal kinase; MAPKAPK, mitogen-activated protein kinase-activated protein kinase; PEA, phosphoprotein enriched in astrocytes; PKC, protein kinase C; STAT, signal transducer and activator of transcription; HNRNPK, heterogeneous nuclear ribonucleoprotein K; CK, casein kinase. ![]()
* This work was supported in part by grants from the Heart and Stroke Foundation of Canada and the Muscular Dystrophy Association (United States) (to L. A. M.). The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact. ![]()
S The on-line version of this article (available at http://www.mcponline.org) contains supplemental material. ![]()
¶ Supported by funding from the Genomic and Health Initiative. ![]()
** Holds the Mach-Gaennelsen Chair in Cardiac Research. To whom correspondence should be addressed: Ontario Genomics Innovation Centre, Ottawa Health Research Institute, 501 Smyth Rd., Ottawa, Ontario K1H 8L6, Canada. Tel.: 613-737-8899 (ext. 78618); Fax: 613-737-8803; E-mail: lmegeney{at}ohri.ca
| REFERENCES |
|---|
|
|
|---|
.
Cell. Signal.
18, 1461
1472[CrossRef][Medline]This article has been cited by other articles:
![]() |
W. Li, G. Wu, and Y. Wan The dual effects of Cdh1/APC in myogenesis FASEB J, November 1, 2007; 21(13): 3606 - 3617. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| All ASBMB Journals | Journal of Biological Chemistry |
| Journal of Lipid Research | ASBMB Today |