Evolutionary traces decode molecular mechanism behind fast pace of myosin XI
© Syamaladevi and Sowdhamini; licensee BioMed Central Ltd. 2011
Received: 20 March 2011
Accepted: 26 September 2011
Published: 26 September 2011
Cytoplasmic class XI myosins are the fastest processive motors known. This class functions in high-velocity cytoplasmic streaming in various plant cells from algae to angiosperms. The velocities at which they process are ten times faster than its closest class V homologues.
To provide sequence determinants and structural rationale for the molecular mechanism of this fast pace myosin, we have compared the sequences from myosin class V and XI through Evolutionary Trace (ET) analysis. The current study identifies class-specific residues of myosin XI spread over the actin binding site, ATP binding site and light chain binding neck region. Sequences for ET analysis were accumulated from six plant genomes, using literature based text search and sequence searches, followed by triple validation viz. CDD search, string-based searches and phylogenetic clustering. We have identified nine myosin XI genes in sorghum and seven in grape by sequence searches. Both the plants possess one gene product each belonging to myosin type VIII as well. During this process, we have re-defined the gene boundaries for three sorghum myosin XI genes using fgenesh program.
Molecular modelling and subsequent analysis of putative interactions involving these class-specific residues suggest a structural basis for the molecular mechanism behind high velocity of plant myosin XI. We propose a model of a more flexible switch I region that contributes to faster ADP release leading to high velocity movement of the algal myosin XI.
Myosin is an actin based motor protein that generates motion using chemical energy released through ATP hydrolysis. Myosins play many important roles within plant cells such as organelle trafficking [1, 2], remodelling [3, 4], and inheritance . They are also known to be involved in development of various plant parts like root hairs, pollen etc[6, 7]. Though there are around 24 classes of myosins reported in eukaryota , only three classes - class VIII, XI and XIII are seen in plants. The similarity of plant myosin sequences with animal and fungal class V myosins , derived from phylogenetic analysis, suggests a common ancestor from which plants and Opisthokonts  might have evolved.
Myosins, in general, function through an ATP hydrolysis cycle by converting the hydrolysis energy to allosteric conformational changes within the motor head as well as the neck region leading to motion. Swinging cross-bridge hypothesis, proposed by H. E. Huxley, has been the most popular model to explain the molecular mechanism of energy transduction in myosins. Numerous biochemical and biophysical experiments thereafter helped to improve it to the present day swinging lever arm model , according to which, immediately after the release of ADP from the previous cycle, ATP binds to myosin head that is in an actin-bound post-stroke conformation. Upon hydrolysis of ATP, myosin head gets transformed to pre-stroke conformation which is actually actin unbound form. Upon rebinding to the next actin molecule, Pi is released first and the ADP bound myosin head changes conformation from weak to strong actin binding state. This will be followed by ADP release and conformational changes at the head domain back to post-stroke state where myosin strongly interacts with actin.
In plants, from algae to angiosperms, the high velocity cytoplasmic streaming (of the range 40-60 μm/s) is driven by myosins . Myosin XI, found in plants, is the fastest known motor and it moves processively along actin filaments towards the plus-end, performing cellular functions like cytoplasmic streaming and vesicle transport [1, 2]. Myosin XI is architecturally similar to class V myosin in animals, with a motor head followed by six IQ motifs, a coiled-coil and a globular cargo-binding domain called DIL . Myosin XI, just like myosin V, functions as a dimer formed through coiled-coil interaction between the α-helical tail regions of the monomers. Both the classes have comparable step size (an average of 35nm) and same directionality towards the plus end of actin filament.
Due to the difficulties in crystallization of actin-myosin complex, the actual actin binding residues are not known even in well-characterized myosins like myosin II and V. However, docking studies have revealed the actomyosin interface residues (see ref.  for a review). Actomyosin interface is extensive in the rigour state because of the interaction of a single head with regions on two adjacent actin molecules. Rigour state contact between actin and myosin head can be divided into four regions: a large primary binding site on the face of actin, which on three sides, is flanked by three additional sites from surface loops . In this study, using the Evolutionary Trace method, we have identified crucial residues at the actin binding site, at the ATP binding domain and at the beginning of neck region that could contribute to the fast release of ADP and the high velocity. During the process of accumulation of myosin XI sequences from plants, we have recognised nine Myosin XI sequences from sorghum and seven from grape through genome-wide survey and gene prediction.
Algal myosin XI, isolated from Chara corallina, slides F-actin in vitro at a speed equivalent to cytoplasmic streaming speed of 40- 60 μm/s, which is 10 times the speed generated by myosin V [16, 17]. Studies with tobacco myosin XI heavy isoform, by Tominaga and coworkers, revealed that single myosin XI molecules move at velocity 7 μm/s along the actin and generate relatively smaller force, in the order of 0.5 pN  much smaller than the force generated by muscle myosin II  and by myosin V [19–21]. In an attempt to elucidate the mechanism of this fast movement of Chara myosin, Ito and co-workers measured its kinetic properties . The rate constant of ADP dissociation from actin-motor domain complex of Chara myosin XI was estimated to be more than 2800 s-1 and the rate constant of ATP-induced dissociation of motor domain from the actin was 2200 s-1 at a physiological concentration of ATP. The estimated time spent on actin, in strongly bound state, was estimated to be <0.82ms. This value is the shortest among known values for various myosins and it has a duty ratio of <0.3 and a Vmax of actin-activated ATPase activity of 390 s-1. ADP release, which is the rate limiting step in all other myosin types, is dramatically accelerated in this plant myosin. Most of the myosins possess positively charged residues on loop 2 where as Chara myosin has no net positive charge on loop2. Instead, positively charged residues are harboured on loop 3. Ito and coworkers investigated the effect of positively charged residues on the loops at the actin binding region and provided evidence for its partial role on the high velocity movement through mutation studies . Still, the sequence signatures that lead to the specialization of these myosins as the fastest motors and the actual molecular mechanism of such a rapid process are not known completely. In this study, based on sequence analysis and molecular modelling, we propose that the sequence signatures at the switch I region, the ATP-binding site and the neck region as partly responsible for the observed high rate of ADP release, which in turn lead to the specialization of these myosins as the fastest ones.
Sequence search and validation
Evolutionary Trace (ET) analysis
Class-specific residues were identified using ET method for myosin V and XI members. Multiple sequence alignment for 30 sequences from 10 genomes was done using CLUSTALW . ET analysis was initially guided by Trace Suite II . Sequences were divided into three parts based on domain architecture: 1) N-terminal SH3-like domain 2) motor domain and 3) light chain binding neck region. Class-specific residues for the motor domain and light chain binding neck domains were mapped on to myosin V crystal structure (PDB ID: 2DFS) and clustering of conserved residues were analysed.
Molecular modelling of algal myosin XI
Homology model of the algal myosin XI (Chara corallina) was obtained using MODELLER . Myosin V from chicken (PDB code: 1W7J) was used as a template . Sequence identity between algal and chicken myosin head domains was 41%. Further, the predicted secondary structures of algal myosin head were comparable with the secondary structures identified in the template (chicken myosin head), except for fine-tuned boundaries. The best energy model from a set of 20 models, generated by MODELLER, was passed through various validation processes viz. WHATIF  and VERIFY-3D . The model was energy minimized to relieve any inherent strain in the model and to obtain the most favourable conformation including side chains. The energy minimized structures were used to visualize the class-specific residues using PYMOL .
Results and discussion
Sequence search and validation identify putative myosins in soghum and grape
Details of confirmed myosins from sorghum and grape after multi-fold validation
1597682 - 1630481
10420215 - 10531891
8436820 - 8456586.
1674607 - 1696216.
1332352 - 1348606
360176 - 377881
6130176 - 6182090
3770898 - 3797194
186074 - 205207
7057755 - 7066246
14576408 - 14597301
52499134 - 52510404
61226792 - 61232762
64687149 - 64704955
66880218 - 66891041
1 to 566
XP_002440148 & XP_002440149
56017722 - 56023175 & 56024467-56037464
8061068 & 8061069
610 & 1399
8229550 - 8241900.
19620879 - 19630164
Both sorghum and grape genomes contain only one copy of myosin VIII type in chromosome 1 and 19, respectively (Table 1).
Class-specific residues in myosin XI lines actin binding site
It was recently shown that the ATP-induced motor head dissociation in myosin XI occurs at a very fast rate of 2200s-1, resulting in stunningly fast movement of myosin XI . In our analysis, we identified two myosin XI-specific residues in the primary binding site (helix-loop-helix) that corresponds to K502 and K514 of chicken myosin V. In myosin XI, K502 and K514 are mutated in a class-specific manner to proline and methionine, respectively. This may introduce a structural kink and local rigidity to the loop and disrupt the charged interaction of the primary binding site to actin. This can help the myosin XI motor head to get unbound sooner from the actin upon ATP binding after the release of Pi and ADP in the previous ATP-hydrolysis cycle.
Apart from the three key actin binding sites, surrounding residues may also get involved in actin binding. Few of the class-specific residues observed near the actin binding site may also contribute to binding strength (yellow spheres at the actin binding -region in Figure 5). The second actin binding site on myosin XI also harbours a class-specific residue: Lysine, but the mutation R542K seems to be too subtle to introduce a difference in the interaction with actin.
Putative involvement of two myosin XI specific residues in switch I flexibility
The switch I region, which drives an open-close conformational change in the myosin head during ATP hydrolysis , has two class-specific residues (shown in green in Figure 5). One of these is Asn 214 (algal myosin XI numbering), corresponding to Asp215 of myosin V. In myosin V, Asp215 has a favourable charged interaction with Lys668, whereas in myosin XI, due to changes of Asp215 to Asn214 and Lys688 to Leu640, this stabilizing interaction is likely abolished which could lead to a more flexible switch I region. These class-specific residues are mapped on the three-dimensional model of algal myosin head (Figure 6A). Since the movement at the switch I region is intimately associated with γ-phosphate binding, coordination of Mg2+ and there by ADP binding and release, we propose a more flexible switch I region in myosin XI that might facilitate an early or easy release of ADP.
Shutter residues might make the ATP binding cavity closure inefficient in Chara coralina
A closer investigation of the ATP binding site (Figure 6C & 6D) of the myosin V crystal structure, reveal two mutations in plant myosins that would probably decide the opening and closing of the small cleft at the ATP binding site. Asp297 is located in an upper loop at the face of the cleft and has a potential polar interaction partner Gln115 located in a lower loop of myosin V. In Figure 6A, this interaction is highlighted through a dotted line in the near-rigour-open conformation, as seen in the crystal structure of myosin V. This interaction can form a vertical shutter for the ATP binding cavity. These 'shutter' residues are Ala296 and Arg111 in myosin XI from Chara coralina, which are stabilised by weak dispersion forces, and hence the closure of the cleft may become inefficient in Chara myosin.
Weakened interaction between the neck and the converter domain affects lever arm tilting
We have also examined the neck region (also termed as 'lever arm') for class-specific residues and the conservation of amino acid types. We identified a cluster of three class-specific residues, two of which (Ile762 and Lys766, chicken myosin V numbering) have potential interactions with the converter domain in the case of myosin V (coloured yellow in Figure 5, Additional files 1 and 2).
Ile762 (in chicken myosin V, please see Additional file 2) at the beginning of the neck region, is embedded between two hydrophobic residues - a highly conserved Leu714 and a conserved hydrophobic residue at position 735 - of the converter domain forming a "hydrophobic triad". The corresponding residue of Ile762 in plant myosin XI neck region is a charged arginine that possibly weaken the hydrophobic interaction of the N-terminal neck with the converter domain.
Spatially close to this "hydrophobic triad" in myosin V (please see Additional files), the side chain of Lys766, the second class specific residue, interacts with the backbone oxygen of Val713 located at the converter domain by forming a hydrogen bond. In plant myosin XI, Lys766 is replaced by a Valine in a class specific manner. As a result of this change, a weak hydrophobic interaction alone could prevail between the neck and the converter domains of myosin XI (Additional file 1). Such amino acid differences could provide higher conformational flexibility to plant myosins in tilting the lever arm (neck region) with ease, which in turn is associated with a faster ADP release and actin dissociation. The third class-specific residue (not shown in Additional file 1) is likely to interact with SH3 domain.
From this comparative sequence analysis, we could recognise most crucial residue changes that may be responsible for decreased/altered affinity of myosin XI to the actin and rapid ADP release in comparison with myosin V. A molecular mechanism of fast movement of myosin XI could be proposed: As long as ADP is bound to the myosin head, actin binding is stabilized through numerous weak interactions. Due to the enhanced flexibility of the switch I region, ADP gets released very fast after hydrolysis. Release of ADP permits the next ATP to bind and subsequently dissociate the acto-myosin complex through suitable conformational changes at the actin-myosin interface, which would further be enhanced due to the absence of few stabilizing charged interactions. Meanwhile, since the neck region in myosin XI interacts only poorly with the converter domain, the neck flexibility is also high which might assist in faster catalysis and binding events at the head domain.
DPS's stay at NCBS is financially supported by NCBS and Sugarcane Breeding Institute, Coimbatore. This work is supported by Human Frontier Science Program. We also thank NCBS (TIFR) for infrastructural support.
- Boevink P, Oparka KJ, Cruz SS, Martin B, Betteridge A, Hawes C: Stacks on tracks: the plant Golgi apparatus traffics on an actin/ER network. Plant J 1998, 15: 441–447. 10.1046/j.1365-313X.1998.00208.xView ArticlePubMedGoogle Scholar
- Kim H, Park M, Kim SJ, Hwang I: Actin filaments play a critical role in vacuolar trafficking at the Golgi complex in plant cells. Plant Cell 2005, 17: 888–902. 10.1105/tpc.104.028829PubMed CentralView ArticlePubMedGoogle Scholar
- Runions J, Brach T, Kuhner S, Hawes C: Photoactivation of GFP reveals protein dynamics within the endoplasmic reticulum membrane. J Exp Bot 2006, 57: 43–50.View ArticlePubMedGoogle Scholar
- Nebenfuhr A, Gallagher LA, Dunahay TG, Frohlick JA, Mazurkiewicz AM, Meehl JB, Staehelin L: Stop-and-go movements of plant Golgi stacks are mediated by the acto myosin system. Plant Physiol 1999, 121: 1127–1141. 10.1104/pp.121.4.1127PubMed CentralView ArticlePubMedGoogle Scholar
- Sheahan MB, Rose RJ, McCurdy DW: Organelle inheritance in plant cell division: the actin cytoskeleton is required for unbiased inheritance of chloroplasts, mitochondria and endoplasmic reticulum in dividing protoplasts. Plant J 2004, 37: 379–390.View ArticlePubMedGoogle Scholar
- Ojangu EL, Jarve K, Paves H, Truve E: Arabidopsis thaliana myosin XIK is involved in root hair as well as trichome morphogenesis on stems and leaves. Protoplasma 2007, 230: 193–202. 10.1007/s00709-006-0233-8View ArticlePubMedGoogle Scholar
- Jiang SY, Cai M, Ramachandran S: Oriza sativa myosin XI B controls pollen development by photoperiod-sensitive protein localizations. Dev Biol 2007, 304: 579–592. 10.1016/j.ydbio.2007.01.008View ArticlePubMedGoogle Scholar
- Tominaga M, Kojima H, Yokota E, Orii H, Nakamori R, Katayama E, Anson M, Shimmen T, Oiwa K: Higher plant myosin XI moves processively on actin with 35 nm steps at high velocity. The EMBO Journal 2003, 22(6):1263–1272. 10.1093/emboj/cdg130PubMed CentralView ArticlePubMedGoogle Scholar
- Desnos C, Huet S, Darchen F: 'Should I stay or should I go?': myosin V function in organelle trafficking. Biol Cell 2007, 99: 411–423. 10.1042/BC20070021View ArticlePubMedGoogle Scholar
- Foth BJ, Goedecke MC, Soldati D: New insights into myosin evolution and classification. Proc Natl Acad Sci USA 2006, 103: 3681–3686. 10.1073/pnas.0506307103PubMed CentralView ArticlePubMedGoogle Scholar
- Holmes KC: The swinging lever-arm hypothesis of muscle contraction. Current Biology 1997, 7(2):112–118. 10.1016/S0960-9822(06)00055-8View ArticleGoogle Scholar
- Nagai R, Hayama T: Ultra structure of the edoplasmic factor responsible for cytoplasmic streaming in Chara internodal cells. J Cell Sci 1989, 36: 121–136.Google Scholar
- Milligan RA: protein _protein interactions in the rigour actomyosin complex. Proc Natl Acad Sci USA 1996, 93: 21–26. 10.1073/pnas.93.1.21PubMed CentralView ArticlePubMedGoogle Scholar
- Reyment I, Holden HM, Whittaker M, Yohn CB, Lorenz M, Holmes KC, Milligan RA: Structure of the actin-myosin complex and its implications for muscle contraction. Science 1993, 261: 58–65. 10.1126/science.8316858View ArticleGoogle Scholar
- Yamamoto K, Kikuyama M, Sutoh-Yamamoto N, Kamitsubo E: Purification of actin based motor protein from Chara corallina . Proc Japan Acad Series B 1994, 70: 175–180. 10.2183/pjab.70.175View ArticleGoogle Scholar
- Anson M: Temperature dependence and Arrhenius activation energy of F-actin velocity generated invitro by skeletal myosin. J Mol Biol 1992, 224: 1029–1038. 10.1016/0022-2836(92)90467-XView ArticlePubMedGoogle Scholar
- Pierobon P, Achouri S, Courty S, Dunn AR, Spudich JA, Dahan M, Cappello G: Velocity, processivity, and individual steps of single myosin V molecules in live cells. Biophys J 2009, 96(10):4268–4275. 10.1016/j.bpj.2009.02.045PubMed CentralView ArticlePubMedGoogle Scholar
- Guilford WH, Dupuis DE, Kennedy G, Wu J, Patlak JB, Warshaw DM: Smooth muscle and skeletal muscle myosins produce similar unitary forces and displacements in the laser trap. Biophys J 1997, 72: 1006–1021. 10.1016/S0006-3495(97)78753-8PubMed CentralView ArticlePubMedGoogle Scholar
- Rief M, Rock RS, Mehta AD, Mooseker MS, Cheney RE, Spudich JA: Myosin-V stepping kinetics: A molecular model for processivity. Proc Natl Acad Sci USA 2000, 97: 9482–9486.PubMed CentralView ArticlePubMedGoogle Scholar
- Veigel C, Wang F, Bartoo ML, Sellers JR, Molloy JE: The gated gait of theprocessive molecular motor, myosin V. Nature Cell Biol 2002, 4: 59–65. 10.1038/ncb732View ArticlePubMedGoogle Scholar
- Higayashi-fujime S, Ishikawa R, Iwasawa H, Kagami O, Kuri-moto E, Kohama K, Hozumi T: The fastest-actin-based motor protein from the green algae, Chara and its distinct mode of interaction with actin. FEBS lett 1995, 375: 151–154. 10.1016/0014-5793(95)01208-VView ArticleGoogle Scholar
- Ito K, Ikebe M, Kashiyama T, Mogai T, Kon T, Yamamoto K: Kinetic Mechanism of the Fastest Motor Protein, Chara Myosin. Journal of Biological Chemistry 2007, 282(27):19534–19545. 10.1074/jbc.M611802200View ArticlePubMedGoogle Scholar
- Ito K, Yamaguchi Y, Yanase K, Ichikawa Y, Yamamoto K: Unique charge distribution in surface loops confers high velocity on the fast motor protein Chara myosin. Proc Natl Acad Sci USA 2009, 106(51):21585–21590. 10.1073/pnas.0910787106PubMed CentralView ArticlePubMedGoogle Scholar
- Finn RD, Mistry J, Tate J, Coggill P, Heger A, Pollington JE, Gavin OL, Gunesekaran P, Ceric G, Forslund K, Holm L, Sonnhammer EL, Eddy SR, Bateman A: The Pfam protein families database. Nucleic Acids Res 2010, 38: D211–222. 10.1093/nar/gkp985PubMed CentralView ArticlePubMedGoogle Scholar
- Marchler-Bauer A, Bryant SH: CD-Search: protein domain annotations on the fly. Nucleic Acids Res 2004, 32: W327–31. 10.1093/nar/gkh454PubMed CentralView ArticlePubMedGoogle Scholar
- Bryson K, McGuffin LJ, Marsden RL, Ward JJ, Sodhi JS, Jones DT: Protein structure prediction servers at University College London. Nucl Acids Res 2005, 33: W36–38. 10.1093/nar/gki410PubMed CentralView ArticlePubMedGoogle Scholar
- Lupas A, Dyke MV, Stock J: Predicting Coled Coils from Protein Sequences. Science 1991, 252: 1162–1164. 10.1126/science.252.5009.1162View ArticlePubMedGoogle Scholar
- Salamov AA, Solovyev VV: Ab initio Gene Finding in Drosophila Genomic DNA. Genome Res 2000, 10: 516–522. 10.1101/gr.10.4.516PubMed CentralView ArticlePubMedGoogle Scholar
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, Thompson JD, Gibson TJ, Higgins DG: Clustal W and Clustal X version 2.0. Bioinformatics 2007, 23(21):2947–2948. 10.1093/bioinformatics/btm404View ArticlePubMedGoogle Scholar
- Felsenstein J: PHYLIP - Phylogeny Inference Package (Version 3.2). Cladistics 1989, 5: 164–166.Google Scholar
- Innis CA, Shi J, Blundell TL: Evolutionary trace analysis of TGF and related growth factors: implications for site-directed mutagenesis. Protein Eng 2000 13(12):839–847.
- Sali A, Blundell TL: Comparative protein modelling by satisfaction of spatial restraints. J Mol Biol 1993, 234: 779–815. 10.1006/jmbi.1993.1626View ArticlePubMedGoogle Scholar
- Coureux PD, Sweeney HL, Houdusse A: Three myosin V structures delineate essential features of chemo-mechanical transduction. Embo J 2004, 23: 4527. 10.1038/sj.emboj.7600458PubMed CentralView ArticlePubMedGoogle Scholar
- WHATIF server[http://swift.cmbi.ru.nl/servers/html/index.html]
- Eisenberg D, Lüthy R, Bowie JU: VERIFY3D: assessment of protein models with three-dimensional profiles. Methods in enzymology 1997, 277: 396–404.View ArticlePubMedGoogle Scholar
- The PyMOL Molecular Graphics System, Version 1.2r3pre Schrödinger, LLC;
- Geli MI, Lombardi R, Schmelzl B, Riezman H: An intact SH3 domain is required for myosin I-induced actin polymerization. The EMBO Journal 2000, 19: 4281–4291. 10.1093/emboj/19.16.4281PubMed CentralView ArticlePubMedGoogle Scholar
- Reubold TF, Eschenburg S, Becker A, Kull FJ, Manstein DJ: A structural model for actin induced nucleotide release in myosin. Nature Structural Biology 2003, 10(10):826–830. 10.1038/nsb987View ArticlePubMedGoogle Scholar