Skip to main content

Tyrosine phosphorylation of myosin heavy chain during skeletal muscle differentiation: an integrated bioinformatics approach

Abstract

Background

Previously it has been shown that insulin-mediated tyrosine phosphorylation of myosin heavy chain is concomitant with enhanced association of C-terminal SRC kinase during skeletal muscle differentiation. We sought to identify putative site(s) for this phosphorylation event.

Results

A combined bioinformatics approach of motif prediction and evolutionary and structural analyses identified tyrosines163 and 1856 of the skeletal muscle heavy chain as the leading candidate for the sites of insulin-mediated tyrosine phosphorylation.

Conclusion

Our work is suggestive that tyrosine phosphorylation of myosin heavy chain, whether in skeletal muscle or in platelets, is a significant event that may initiate cytoskeletal reorganization of muscle cells and platelets. Our studies provide a good starting point for further functional analysis of MHC phosphor-signalling events within different cells.

Introduction

Myosins, actin-based motor proteins, are expressed as multiple isoforms in all eukaryotic cells. They are oligomers consisting of one or two heavy chains to which one or more light chains are non-covalently attached. Myosins have been classified into 18 families based on the amino acid sequence differences in the N-terminal head domains, which contain highly conserved regions including actin- and nucleotide-binding sites [1, 2]. The tail of myosin is the most variable domain and seems to be responsible for the specific role myosin plays in the cell.

Functional activities of most myosins such as actin-dependent ATPase activity or ability to move actin filaments in vitro are regulated in several ways, mainly by phosphorylation of the regulatory light chain, Ca2+-binding, or phosphorylation of the heavy chain [1, 3] It has been previously claimed that the myosin heavy chain (MHC)undergoes tyrosine phosphorylation during insulin-mediated skeletal muscle differentiation, thus linking signal transduction to highly ordered myosin assembly [4]. Insulin modulates an association of myosin with C-terminal SRC kinase (Csk), a tyrosine kinase signalling molecule, and these interactions are fundamental in skeletal muscle differentiation. Although the claims of tyrosine phosphorylation of MHC in vivo remain somewhat controversial, tyrosine phosphorylation of non-muscle MHC IIa has also been implicated as an early event in human platelet activation [5]. To settle this controversy -and establish the role, if any, of MHC tyrosine phosphorylation it is important to identify sites at which such phosphorylation events may occur.

We have mapped potential phosphorylation sites on the skeletal muscle myosin heavy chain utilizing an integrated bioinformatics approach, supporting web-based motif predictions with evolutionary and structural data. Of all the sites analyzed in the bioinformatics approach, the data suggest Y163 and Y1856 as the leading candidates for insulin-mediated tyrosine phosphorylation.

Methods

Tyrosine Phosphorylation Predictions

Tyrosine phosphorylation site predictions were made with two different online resources using the sequences described below. NetPhos 2.0 produces neural network predictions based on sequence and structure [6]. Scansite predicts target motifs for different kinases using a positional selectivity matrix based on peptide library screening data [7, 8] In addition, Scansite predictions were made for known phosphotyrosine recognition motifs for evidence of downstream signalling events. All Scansite predictions were made on the 'Low Stringency' setting to identify as many putative sites as possible. These sites were then supported or rejected on the basis of further analysis as described below.

Evolutionary Analysis

Protein sequences for adult skeletal muscle myosin heavy chains (MYHSA) 1 and 2 were extracted from the SwissProt database [9] MYHSA1 [SwissProt : MYH1_HUMAN, P12882]; MYHSA2 [SwissProt ID: MYH2_HUMAN, Q9UKX2] and used as query sequences to extract closely related homologous proteins. First, BLAST [10] was used to search SwissProt-TrEMBL [9] and the known, novel and Genscan-predicted peptides of five EnsEMBL genome databases (Human, Mouse, Rat, Fugu, Zebrafish) [11] Redundant sequences were removed and ALIGN [12, 13] was used to make pairwise alignments of each homologue with MYH1_HUMAN and to calculate the percentage identity across the entire length of the protein. Vertebrate homologues with at least 60% global identity were processed using an in-house homologue processing tool, HAQESAC [14]. Homologues were aligned using CLUSTALW [15] and badly-aligned sequences eliminated from the dataset. A neighbour-joining tree with 1000 bootstrap replicates was constructed using CLUSTALW and the sequences were grouped into subfamilies of orthologous proteins. The clade corresponding to skeletal muscle myosin heavy chains in Amniota (mammals, reptiles and birds) were then used as sequences for tyrosine phosphorylation motif prediction as described above.

Secondary Structure Prediction

Secondary structure predictions were made for MYH1_HUMAN using the PSIPRED V2.3 website [16]. Because of the length of the protein, it was submitted in two overlapping chunks: residues 1–814 and 800 +.

3D Structure Analysis

3D structures were obtained from the Protein Data Bank (PDB) [17] and viewed with the RasMol viewer [18]. Three myosin heavy chain structures were identified: 2MYS, Chicken adult skeletal muscle myosin heavy chain; 1BR2, chicken gizzard smooth muscle myosin heavy chain; and 1B7T, Aequipecten irradians (Bay scallop) striated muscle myosin heavy chain. The corresponding SwissProt sequences [Swiss -Prot :2MYS: MYSS_CHICK, P13538]; [Swiss-Prot1BR2: MYHB_CHICK, P10587]; [Swiss -Prot1B7T: MYS_AEQIR, P24733] were downloaded and aligned with Human MYH1_HUMAN and MYH2_HUMAN using CLUSTALW. This alignment was used with the skeletal muscle myosin heavy chains (above) to assign putative tyrosine phosphorylation sites to their corresponding residues in the homologous 3D structures. Visualisation with RasMol and DSSP solvent accessibility data [19] was then used to infer whether potential sites of tyrosine phosphorylation were surface-exposed or buried.

Results

In total, twenty-three myosin heavy chain sequences were used for tyrosine phosphorylation motif prediction, which were divided into five groups of orthologous sequences (Figure 1). Important motifs are likely to be conserved during evolution and so we considered only those sites that were predicted to be phosphorylation motifs in all the sequences of at least one orthologous group. Because phosphorylation site predictors have a tendency to over-predict, we increased stringency by accepting only those motifs that received a NetPhos score of 0.8 or higher, or were predicted by both NetPhos and Scansite, in at least one human adult skeletal myosin heavy chain. This yielded fourteen putative sites (Table 1). Of these, six were predicted by both methods, including two motifs that were conserved across all sequences (MYH1_HUMAN Y163 and Y1856).

Figure 1
figure 1

Neighbour-joining phylogeny of MHC homologues, with bootstrap support. PDB structure 2MYS is marked with a black diamond.

Table 1 Summary of predicted tyrosine phosphorylation sites.

To be phosphorylated, tyrosine residues must be accessible on the surface of the protein. Although the three-dimensional conformations of homologous myosin molecules will not be identical, the high degree of sequence conservation between human adult skeletal muscle myosin heavy chains and the three myosin sequences present in PDB allowed the inference of solvent accessibility. This was confirmed by the generally good agreement in surface accessibility measures both between models and between the different myosin chains of 1BR2 (data not shown). From these data, two sites (Y286 and Y435) were buried while a further two (Y313 and Y504) had very low solvent accessibility (Table 1). 3D data were not available for the four tyrosines in the C-terminal of the protein.

If tyrosine phosphorylation of MMHC-II is part of a signalling cascade, it is likely that some other protein will interact with the phosphotyrosine. We used Scansite to look for phosphotyrosine interaction motifs and found three SH2 domain recognition motifs that matched potentially exposed phosphorylation sites. Because Scansite also identifies the interacting protein, we interrogated the Gene Cards [20] entry for each kinase and SH2 domain for expression patterns. Only four kinases and two SH2 domains had evidence from UniGene [21] or SAGE [22] of expression in skeletal muscle, while only one kinase (INSR) and one SH2 domain (PIK3R1) had evidence from both (Table 2). Interestingly, both of the latter pair were predicted to interact with the same, totally conserved, motif (MYH1_HUMAN Y163). Furthermore, both are involved in insulin-mediated pathways (see Discussion). Two kinases expressed in skeletal muscle were predicted to interact with Y1856. These were an SRC kinase and ABL1, which interacts with SORBS1 following insulin stimulation [23].

Table 2 Interacting enzymes predicted by Scansite.

Discussion

Bioinformatics alone cannot identify a functional motif; supporting experiments will always be needed for conclusive evidence. Nevertheless, while other sites cannot be categorically excluded, the combined data presented here identify Y163 and Y1856 as the most likely sites for tyrosine phosphorylation events in skeletal muscle. Both NetPhos and Scansite predicted these motifs for all mammalian adult skeletal myosin heavy chain sequences analysed, indicating strong evolutionary conservation (Figure 1). A kinase predicted to be responsible for phosphorylation of each site is expressed in skeletal muscle, as was an SH2 domain protein that was predicted to interact with a phosphotyrosine at Y163. Analysis of predicted secondary structures and homologous 3D structures indicates that these sites may be accessible on the protein surface. The position of Y163 as part of an alpha helix within the globular myosin head domain does raise concerns that it is potentially difficult to phosphorylate, even though it is on the surface of the domain. Nevertheless, depending on the relative conformations of the solved in vitro chicken myosin structures compared to in vivo human myosin, Y163 might still be available for phosphorylation. Y1856 is in region of low predicted secondary structure (Table 1), indicative of a flexible loop region more usually associated with phosphorylation sites.

Myosin heavy chain (MHC) undergoes tyrosine phosphorylation during insulin-mediated differentiation in skeletal muscles and the degree of phosphorylation increases in line with differentiation [4]. Interestingly, for the strongest candidate tyrosine phosphorylation site, Y163, both the kinase and interacting SH2 domain predicted by Scansite are involved in insulin-mediated pathways. The kinase INSR is a transmembrane receptor that binds insulin [24] while the SH2 domain protein PIK3R1 is necessary for the insulin-stimulated increase in glucose uptake and glycogen synthesis in insulin-sensitive tissues [25]. We can therefore conclude that Y163 remains a strong candidate site for insulin-mediated tyrosine phosphorylation of myosin heavy chain, despite concerns over accessibility.

As phosphorylation sites are often in the tails of proteins, the tyrosines outside the main globular domains, namely Y1379, Y1492 and Y1856, are also potential candidates for phosphorylation. The strongest of these is Y1856, which is both C-terminal and predicted to be phosphorylated by the kinases SRC and ABL1, which are found in skeletal muscle (Table 2). Perhaps of most interest is ABL1, a protein known to be associated with "Sorbin and SH3 domain containing 1" (SORBS1) during insulin signalling in other cell lines [23]. SORBS1 is highly expressed in skeletal muscle (data not shown) and is involved in formation of actin stress fibres and focal adhesions; its orthologue, CAP, has been identified as an important adaptor during insulin signalling in mice [26–28]. Furthermore, the SORBS1 gene has been implicated in the pathogenesis of human disorders with insulin resistance [29].

Csk has been shown to be associated with the hormone 1, 25(OH)2-vitamin D3 resulting in the stimulation of the growth-related mitogen-activated protein kinase (MAPK). The phosphorylated form of MAPK is then translocated to the nucleus where it induces the expression of c-myc oncoprotein associated with skeletal muscle proliferation [30]. In addition, Csk has also been implicated in the regulation of integrins and the control of cell attachment and shape [31] Goel et al. showed that insulin can phosphorylate myosin, leading to an association with Csk and thus to a decrease in c-Src activity. This has also been shown in fibroblast cell lines following stimulation of the insulin-like growth factor-I receptor [32]. This demonstrates the potential for skeletal muscle differentiation after phosphorylation of Y163 and/or Y1856 of the MHC.

Harney et al. have shown that non-muscle myosin heavy chain type IIA in platelets undergoes tyrosine phosphorylation and subsequent dephosphorylation in a time-dependent manner [5]. In common with other cells, the cytoskeleton of platelets comprises actin filaments, microtubules and myosin molecules. Myosins form rings within the platelet that maintain a spherical shape and several lines of evidence suggest that these rings reorient following platelet activation to permit spreading [33, 34]. While myosin function therefore appears critical to platelet spreading, studies using cytoskeleton inhibitors have shown that at least the early events of platelet activation are not dependent on the cytoskeletal changes [35]. Our work is suggestive that tyrosine phosphorylation of myosin heavy chain, whether in skeletal muscle or in platelets, is a significant event that may initiate cytoskeletal reorganization of muscle cells and platelets. Our studies provide a good starting point for further functional analysis of MHC phosphor-signalling events within different cells.

Supplementary Information

Full prediction results, sequence alignments and links to in-house software used can be found at: http://www.bioinformatics.rcsi.ie/~redwards/phos/

References

  1. Korn ED: Coevolution of head, neck, and tail domains of myosin heavy chains. Proc Natl Acad Sci U S A. 2000, 97: 12559-12564. 10.1073/pnas.230441597.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  2. Yamashita RA, Sellers JR, Anderson JB: Identification and analysis of the myosin superfamily in Drosophila: a database approach. J Muscle Res Cell Motil. 2000, 21: 491-505. 10.1023/A:1026589626422.

    Article  CAS  PubMed  Google Scholar 

  3. Sellers J: Myosins (Protein Profile S). 1999, Oxford University Press, Oxford

    Google Scholar 

  4. Goel HL, Dey CS: Insulin-mediated tyrosine phosphorylation of myosin heavy chain and concomitant enhanced association of C-terminal SRC kinase during skeletal muscle differentiation. Cell Biol Int. 2002, 26: 557-561. 10.1006/cbir.2002.0878.

    Article  CAS  PubMed  Google Scholar 

  5. Harney D, Cagney G, Thiede B, Treumann A, Fitzgerald DJ, Maguire PB: Proteomic analysis of platelets detects the May-Hegglin gene, non-muscle myosin heavy chain A in platelets and demonstrates its tyrosine phosphorylation following thrombin stimulation. [abstract]. Circulation. 2003, 108: s200-

    Google Scholar 

  6. Blom N, Gammeltoft S, Brunak S: Sequence and structure-based prediction of eukaryotic protein phosphorylation sites. J Mol Biol. 1999, 294: 1351-1362. 10.1006/jmbi.1999.3310.

    Article  CAS  PubMed  Google Scholar 

  7. Songyang Z, Blechner S, Hoagland N, Hoekstra MF, Piwnica-Worms H, Cantley LC: Use of an oriented peptide library to determine the optimal substrates of protein kinases. Curr Biol. 1994, 4: 973-982. 10.1016/S0960-9822(00)00221-9.

    Article  CAS  PubMed  Google Scholar 

  8. Obenauer JC, Cantley LC, Yaffe MB: Scansite 2.0: Proteome-wide prediction of cell signaling interactions using short sequence motifs. Nucleic Acids Res. 2003, 31: 3635-3641. 10.1093/nar/gkg584.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Boeckmann B, Bairoch A, Apweiler R, Blatter MC, Estreicher A, Gasteiger E, Martin MJ, Michoud K, O'Donovan C, Phan I, Pilbout S, Schneider M: The SWISS-PROT protein knowledgebase and its supplement TrEMBL in 2003. Nucleic Acids Res. 2003, 31: 365-370. 10.1093/nar/gkg095.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1996, 215: 403-410. 10.1006/jmbi.1990.9999.

    Article  Google Scholar 

  11. Clamp M, Andrews D, Barker D, Bevan P, Cameron G, Chen Y, Clark L, Cox T, Cuff J, Curwen V, Down T, Durbin R, Eyras E, Gilbert J, Hammond M, Hubbard T, Kasprzyk A, Keefe D, Lehvaslaiho H, Iyer V, Melsopp C, Mongin E, Pettett R, Potter S, Rust A, Schmidt E, Searle S, Slater G, Smith J, Spooner W, Stabenau A, Stalker J, Stupka E, Ureta-Vidal A, Vastrik I, Birney E: Ensembl 2002: accommodating comparative genomics. Nucleic Acids Res. 2003, 31: 38-42. 10.1093/nar/gkg083.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  12. Pearson WR: Using the FASTA program to search protein and DNA sequence databases. Methods Mol Biol. 1994, 24: 307-331.

    CAS  PubMed  Google Scholar 

  13. Pearson WR: Flexible sequence similarity searching with the FASTA3 program package. Methods Mol Biol. 2000, 132: 185-219.

    CAS  PubMed  Google Scholar 

  14. Edwards RJ: HAQESAC: Homologue Alignment Quality, Establishment of Subfamilies and Ancestor Construction. 2004,http://www.bioinformatics.rcsi/~redwards/haqesac/

    Google Scholar 

  15. Higgins DG, Sharp PM: CLUSTAL: a package for performing multiple sequence alignment on a microcomputer. Gene. 1988, 73: 237-244. 10.1016/0378-1119(88)90330-7.

    Article  CAS  PubMed  Google Scholar 

  16. McGuffin LJ, Bryson K, Jones DT: The PSIPRED protein structure prediction server. Bioinformatics. 2000, 16: 404-405. 10.1093/bioinformatics/16.4.404.

    Article  CAS  PubMed  Google Scholar 

  17. Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE: The Protein Data Bank. Nucleic Acids Res. 2000, 28: 235-242. 10.1093/nar/28.1.235.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  18. Sayle RA, Milner-White EJ: RASMOL: biomolecular graphics for all. Trends Biochem Sci. 1995, 20: 374-10.1016/S0968-0004(00)89080-5.

    Article  CAS  PubMed  Google Scholar 

  19. Kabsch W, Sander C: Dictionary of protein secondary structure: pattern recognition of hydrogen-bonded and geometrical features. Biopolymers. 1983, 22: 2577-2637. 10.1002/bip.360221211.

    Article  CAS  PubMed  Google Scholar 

  20. Rebhan M, Chalifa-Caspi V, Prilusky J, Lancet D: GeneCards: encyclopedia for genes, proteins and diseases. 1997,http://bioinformatics.weizmann.ac.il/cards

    Google Scholar 

  21. Schuler GD, Boguski MS, Stewart EA, Stein LD, Gyapay G, Rice K, White RE, Rodriguez-Tome P, Aggarwal A, Bajorek E, Bentolila S, Birren BB, Butler A, Castle AB, Chiannilkulchai N, Chu A, Clee C, Cowles S, Day PJ, Dibling T, Drouot N, Dunham I, Duprat S, East C, Hudson TJ: A gene map of the human genome. Science. 1996, 274: 540-546. 10.1126/science.274.5287.540.

    Article  CAS  PubMed  Google Scholar 

  22. Velculescu VE, Zhang L, Vogelstein B, Kinzler KW: Serial analysis of gene expression. Science. 1995, 270: 484-487.

    Article  CAS  PubMed  Google Scholar 

  23. Lin WH, Huang CJ, Liu MW, Chang HM, Chen YJ, Tai TY, Chuang LM: Cloning, mapping, and characterization of the human sorbin and SH3 domain containing 1 (SORBS1) gene: a protein associated with c-Abl during insulin signaling in the hepatoma cell line Hep3B. Genomics. 2001, 74: 12-20. 10.1006/geno.2001.6541.

    Article  CAS  PubMed  Google Scholar 

  24. White MF: The IRS-signaling system: a network of docking proteins that mediate insulin and cytokine action. Recent Prog Horm Res. 1998, 53: 119-138.

    CAS  PubMed  Google Scholar 

  25. Hers I, Bell CJ, Poole AW, Jiang D, Denton RM, Schaefer E, Tavare JM: Reciprocal feedback regulation of insulin receptor and insulin receptor substrate tyrosine phosphorylation by phosphoinositide 3-kinase in primary adipocytes. Biochem J. 2002, 368: 875-884. 10.1042/BJ20020903.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. Ribon V, Herrera R, Kay BK, Saltiel AR: A role for CAP, a novel, multifunctional Src homology 3 domain-containing protein in formation of actin stress fibers and focal adhesions. J Biol Chem. 1998, 273: 4073-4080. 10.1074/jbc.273.7.4073.

    Article  CAS  PubMed  Google Scholar 

  27. Ribon V, Printen JA, Hoffman NG, Kay BK, Saltiel AR: A novel, multifuntional c-Cbl binding protein in insulin receptor signaling in 3T3-L1 adipocytes. Mol Cell Biol. 1998, 18: 872-879.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. Baumann CA, Ribon V, Kanzaki M, Thurmond DC, Mora S, Shigematsu S, Bickel PE, Pessin JE, Saltiel AR: CAP defines a second signalling pathway required for insulin-stimulated glucose transport. Nature. 2002, 407: 202-207.

    Google Scholar 

  29. Lin WH, Chiu KC, Chang HM, Lee KC, Tai TY, Chuang LM: Molecular scanning of the human sorbin and SH3-domain-containing-1 (SORBS1) gene: positive association of the T228A polymorphism with obesity and type 2 diabetes. Hum Mol Genet. 2001, 10: 1753-1760. 10.1093/hmg/10.17.1753.

    Article  CAS  PubMed  Google Scholar 

  30. Buitrago C, Boland R, de Boland AR: The tyrosine kinase c-Src is required for 1,25(OH)2-vitamin D3 signalling to the nucleus in muscle cells. Biochim Biophys Acta. 2001, 1541: 179-87. 10.1016/S0167-4889(01)00142-2.

    Article  CAS  PubMed  Google Scholar 

  31. Bergman M, Joukov V, Virtanen I, Alitalo K: Overexpressed Csk tyrosine kinase is localized in focal adhesions, causes reorganization of alpha v beta 5 integrin, and interferes with HeLa cell spreading. Mol Cell Biol. 1995, 15: 711-722.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Arbet-Engels C, Tartare-Deckert S, Eckhart W: C-terminal Src kinase associates with ligand-stimulated insulin-like growth factor-I receptor. J Biol Chem. 1999, 274: 5422-5428. 10.1074/jbc.274.9.5422.

    Article  CAS  PubMed  Google Scholar 

  33. Nachmias VT: Platelet and megakaryocyte shape change: triggered alterations in the cytoskeleton. Semin Hematol. 1983, 20: 261-281.

    CAS  PubMed  Google Scholar 

  34. White JG, Radha E, Krumwiede M: Isolation of circumferential microtubules from platelets without simultaneous fixation. Blood. 1986, 67: 873-877.

    CAS  PubMed  Google Scholar 

  35. White JG, Rao GH: Influence of a microtubule stabilizing agent on platelet structural physiology. Am J Pathol. 1983, 112: 207-217.

    PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors would like to thank Dr G. Cagney (GC) and Dr Patricia Maguire (PBM) for many helpful comments during the analysis and preparation of the manuscript. In addition this work was supported in part by a fellowship from Enterprise Ireland (PBM), the Health Research Board of Ireland (PBM) and the Higher Education Authority of Ireland (GC) and by a Science Foundation Ireland award (grant no. 02/IN.1/B117).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to DF Harney.

Additional information

Competing interests

The author(s) declare that they have no competing interests.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Harney, D., Butler, R. & Edwards, R. Tyrosine phosphorylation of myosin heavy chain during skeletal muscle differentiation: an integrated bioinformatics approach. Theor Biol Med Model 2, 12 (2005). https://doi.org/10.1186/1742-4682-2-12

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1742-4682-2-12

Keywords