|Home | About | Journals | Submit | Contact Us | Français|
The role of the actin cytoskeleton in relation to mitochondria function and dynamics is only recently beginning to be recognized. Myo19 is an actin-based motor that is bound to the outer mitochondrial membrane and promotes the localization of mitochondria to filopodia in response to glucose starvation. However, how glucose starvation induces mitochondria localization to filopodia, what are the dynamics of this process and which enzymatic adaptation allows the translocation of mitochondria to filopodia are not known. Here we show that reactive oxygen species (ROS) mimic and mediate the glucose starvation induced phenotype. In addition, time-lapse fluorescent microscopy reveals that ROS-induced Myo19 motility is a highly dynamic process which is coupled to filopodia elongation and retraction. Interestingly, Myo19 motility is inhibited by back-to-consensus-mutation of a unique residue of class XIX myosins in the motor domain. Kinetic analysis of the purified mutant Myo19 motor domain reveals that the duty ratio (time spent strongly bound to actin) is highly compromised in comparison to that of the WT motor domain, indicating that Myo19 unique motor properties are necessary to propel mitochondria to filopodia tips. In summary, our study demonstrates the contribution of actin-based motility to the mitochondrial localization to filopodia by specific cellular cues.
Mitochondria participate in multiple cellular functions including ATP synthesis, calcium buffering, lipid synthesis, signaling and apoptosis1. The distribution of mitochondria throughout the cell is regulated and affected by external stimuli such as stress or growth factors2, 3. Consistent with the critical roles of the mitochondria, aberrant distribution of these organelles is implicated in neurodegenerative diseases4–6. Mitochondria move primarily on microtubule (MT) tracks via kinesin and dynein motors, although some motility is mediated via the actin cytoskeleton7–9.
Myosins are actin-based motors that participate in diverse cellular functions including muscle contractions, cytoskeleton dynamics, membrane tension and cargo transport10, 11. Their involvement in human diseases such as cancer is coming to be appreciated at the molecular level12, 13. The myosin architecture contains a highly conserved motor domain that binds actin and ATP, followed by a ‘neck domain’ that possesses the light–chain-binding motifs, and a highly variable C-terminal domain that interacts with diverse cargo, proteins and membranes14. Despite the high degree of conservation of the motor structure and the conserved ATPase cycle across all myosin classes and isoforms, the rate and equilibrium constants of the biochemical transitions vary greatly, yielding a unique enzymatic adaptation tailored for each individual myosin’s cellular function15–17.
Myosin 19 (Myo19) is a mitochondrial actin-based molecular motor involved in mitochondrial motility, segregation during cell division, and localization to filopodia tips18–20. Myo19 is embedded in the outer mitochondrial membrane via a highly specific and stable interaction mediated by a unique ~45 amino acids motif in its C-terminal domain20, 21.
Filopodia, thin actin-based protrusions in the cell membrane that allow the cell to probe the environment, function in a variety of cellular activities such as cell motility, wound healing, and phagocytosis, and are also precursors of dendritic spines22, 23. Filopodia form by linear polymerization of G-actin-ATP at their barbed–ends mediated by formins at the tip, a process regulated by small Rho family GTPases and promoted by oxidative stress and reactive oxygen species (ROS)23–25. ROS exert deleterious effects on the cell by damaging lipids, proteins and DNA26. However, a physiological level of oxidative stress is important for normal cellular functions, and ROS also serve as secondary messengers in signaling cascades that affect actin cytoskeleton rearrangement, proliferation, differentiation, cell motility and apoptosis. ROS affect signaling pathways by oxidizing reactive cysteins found in the active site of phosphatases and kinases, thereby inhibiting phosphatases and modulating kinases activities26, 27.
We extended our studies of the effect of glucose starvation on Myo19, and found that the localization of Myo19 and mitochondria to filopodia tips is mediated by ROS. Tracking Myo19 localization during ROS-induced filopodia formation, revealed its inducible motor activity that strongly supports active and directional motility of Myo19-bound mitochondria towards filopodia tips. A highly conserved amino acid in the motor domain at the end of P-loop and beginning of Loop 1 (position 140) is unique to class XIX28. Mutation of this position to the consensus amino acid of Human myosins (W140V) affected its motility to the filopodia tips, although its localization to the mitochondria remained unchanged. To study the protein’s motor properties, we co-expressed and purified the mutant Myo19W140V-3IQ with calmodulin as a light chain in a suspension adapted, serum-free, human HEK293SF-3F6 cell line. Steady-state and transient kinetics, as well as kinetic modeling, allowed us to identify key features that differed from the WT and were affected by the mutation. Strikingly, gentle modulation of Myo19 ATPase cycle parameters strongly affected its cellular function as revealed by the Myo19W140V mutant phenotype.
Glucose starvation increases the cellular levels of ROS by increasing ROS generation by both mitochondria and NADPH oxidase (NOX) enzymes29. NOX enzymes reside on the plasma membrane and oxidize extracellular oxygen to generate ROS that enter the cell30, 31. The resultant increase in oxidative stress inhibits phosphatases and induces a specific phospho-tyrosine phosphorylation signature associated with focal adhesions29.
We sought to investigate the connection between the effects of glucose starvation (and elevated ROS) and Myo19 activation, as well as determine whether the starvation-induced localization of Myo19 is indeed mediated by ROS. To visualize filopodia under starvation conditions with or without the addition of antioxidants, we transfected U2OS cells with the fluorescently labeled actin binding protein Emerald-LifeAct (emLifeAct) (Fig. 1a). Consistent with results obtained previously using fascin as a filopodia marker, starvation of the emLifeAct-expressing cells resulted in filopodia formation (17.6±3.9 filopodia per cell, 4.77±2.1µm in length, N=100 from three independent experiments, Fig. 1a,b). These cells formed fewer filopodia than cells overexpressing fascin, likely due to the bundling activity resulting from fascin overexpression20, 32. Interestingly, addition of the extracellular H2O2 scavenger catalase to the starvation media strongly inhibited filopodia formation, indicating that the signal for filopodia formation is mediated by extracellular H2O2 (Fig. 1a, left panel: 2, Fig. 1b). The addition of the common antioxidant propyl gallate (PG) to starvation media had a synergistic effect on filopodia formation, increasing their number but not their length, but had no effect when added to non-starved cells (Fig. 1a, right panel: 3, Fig. 1b). This surprising and counterintuitive finding can be explained by the fact that PG is a dismutase mimic and can increase the rate at which extracellular non-permeable ROS such as superoxide is converted to H2O2, enters the cell, and integrates into signaling pathways33. Consistent with this, the synergistic effect was completely inhibited by supplementation with catalase (Fig. 1a, left panel: 5, Fig. 1b). Taken together, these results support a mechanism in which glucose starvation induces an increase in cellular oxidative stress, which in turn induces filopodia formation. This prompted us to check whether direct addition of H2O2 to the media was sufficient to induce filopodia. Remarkably, H2O2 was able to mimic the glucose starvation–induced formation of filopodia (Fig. 1a, right panel: 2, Fig. 1b), which contained mitochondria (detected via antibody against ATP-synthase subunit alpha – ATP5A) and endogenous Myo19 (magenta and green, respectively, marked with arrow) at their tips (Fig. 1c).
Myo19 is a plus end–directed myosin that slides on actin filaments towards their barbed ends, which are located at filopodia tips23, 34. This property suggests that, without regulation of its motility, Myo19 would reach the tips of all actin protrusions. To determine whether Myo19 localization to filopodia is regulated, we investigated whether Myo19 localizes to tips of filopodia induced by expression of either a constitutively active form of the formin mDia2, mDia2ΔDAD, or Myo1035, 36. We co-transfected U2OS cells with Ruby-tagged Myo19 and GFP-tagged forms of mDia2ΔDAD or Myo10, and followed the localization of Myo19 in these cells. Both mDia2ΔDAD- and Myo10-expressing cells spontaneously formed filopodia in complete media; however, Myo19 localized to filopodia tips only in cells expressing mDia2ΔDAD, but not in those expressing Myo10 (Fig. 2). The selective localization of Myo19 to mDia2ΔDAD-induced filopodia is likely due to differences in the modes of filopodia induction. Dimerization of Myo10 motor domains is sufficient to promote filopodia formation, relying on the bundling activity of the two motor heads37. By contrast, mDia2 interacts with numerous proteins, such as the tyrosine kinase Src, that may initiate a signaling cascade resulting in activation of Myo1938. These results suggest that Myo19 localization to filopodia tips is not the default preference, and that Myo19 localization is limited to certain actin-based protrusions.
Next, we considered two competing possibilities: Myo19 localization to the filopodia tip could be coupled to filopodia elongation, or instead occur after the filopodia has formed. To distinguish between these two mechanisms and reveal more about the dynamics of Myo19 localization, we performed time-lapse fluorescence microscopy. Myo19 and actin were visualized by transfecting U2OS cells with full-length Emerald Green–tagged Myo19 and Ruby-tagged LifeAct. We stimulated the cells with 0.2mM H2O2 and followed the dynamics of Myo19 and actin over time (Movie S1). Initially, we observed ruffling of the cell membrane, followed by formation of filopodia. The filopodia formation was highly dynamic, with multiple formation and retraction events during the first hour, whereas formation of stable filopodia began after 1hour. We observed formation of filopodia and subsequent directional movement of Myo19 towards their tips, indicating that Myo19 actively migrates on actin filaments and tracks the tips (Fig. 3a). Several filopodia with Myo19 at their tips were present at the beginning of imaging. Some of the filopodia (2/4) retracted back to the cell, and this was preceded by Myo19 transport to the cell body (Fig. 3b). Overall, these results reveal that Myo19 can move both to filopodia tips and back to the cell body. The forward movement is motor-dependent, as the ATPase-deficient Myo19G135R mutant cannot localize to filopodia tips20, 28, and Myo19 is a barbed end–directed myosin34. The retrograde motion can be explained by cessation of actin polymerization at the tip and retraction of filopodia back to the cell body with Myo1939, 40. Towards the end of the time-lapse movie (~1.5h), some filopodia seem to be without Myo19 at their tips, but this was most likely due to bleaching, as ~100% of filopodia contained Myo19 at their tips (Fig. 4c). Unfortunately, use of antioxidants to reduce photo-bleaching resulted in inhibition of ROS-induced filopodia formation.
The results of the time-lapse microscopy suggested that the motility of Myo19 towards the barbed end of the actin protrusions is a directional process dependent on the activity of the Myo19 motor. Moreover, if Myo19 only anchored mitochondria to filopodia actin bundles, it would remain in the shaft of filopodia, as the protrusions are elongated by actin polymerization at the tip.
The results described above suggested that the enzymatic activity of the motor propels Myo19 and mitochondria towards filopodia tips. Therefore, we searched for amino acids specific to Myo19 that might provide more insights into the enzymology of the motor domain. To identify candidate amino acids, we performed sequence alignment of human Myo19 against other human myosins, as well as against Myo19 homologues from 40 species41. We and others28 found that W140 is a unique amino acid in Myo19 (Fig. 4a; other myosin classes contain E [18/38] or V [15/38] at this position). This residue is positioned at the edge of the highly conserved P-loop (ATP-binding motif) and at the beginning of the helix preceding Loop 1, a key determinant of several enzymatic properties such as ATPase and motility rate, and nucleotide-binding kinetics42. To determine the role of the unique Trp on Myo19, we performed structural homology modeling using the Phyre2 webserver and compared the result with solved myosin structures43. The modeling suggested that W140 is in close proximity to the bound nucleotide, and thus might modulate some of the parameters related to nucleotide-binding kinetics (Fig. 4b).
To test the importance of W140, we mutated this residue to a consensus Val and investigated whether Myo19W140V would localize to filopodia tips in response to H2O2. As with WT Myo19, Myo19W140V localized to mitochondria under normal growth conditions, but failed to localize to H2O2-induced filopodia tips (Fig. 4c,d). To further assess the importance of W140, we introduced several other point mutations, examining the effect of the residue properties and volume on H2O2-induced filopodia localization of Myo19. When W140 was mutated to Leu, Glu, Ala, or Gly, the mutant protein failed to show any significant localization to filopodia tips. Myo19W140F was the only mutant indistinguishable from the WT, stressing the critical role of a bulky aromatic residue in position 140 for Myo19 function (Fig. S1). Myo19W140V-3IQ, which lacks the C-terminal outer mitochondrial membrane–binding domain, localized to the shaft region of H2O2-induced filopodia (Fig. S2). This observation strongly indicates that Myo19-3IQW140V is still capable of actin motility when not bound to the mitochondria, i.e., under no load. To further assess the contribution of motor activity to transport of mitochondria to filopodia tips, we performed rescue experiments by co-overexpressing Myo19W140V with WT Myo19. Expression of the WT protein with the W140V mutant at a 1:1 ratio resulted in localization of both myosins to filopodia tips, indicating that the WT protein rescued the mutant phenotype, and strongly supporting that Myo19 functions as an ensemble of motors (Fig. 4e).
This result also indicated that the Myo19W140V mutant does not act in a dominant-negative manner when in stoichiometric proportions with the WT. The underlying mechanism for the inability of Myo19W140V to transport the mitochondria as a cargo could be understood by studying its ATPase cycle to reveal key enzymatic properties that were lost or modified. To this end, we co-expressed and purified the motor domain with the lever arm consisting of three IQ calmodulin-binding motifs, as performed previously for the WT (Myo19-3IQ and Myo19W140V-3IQ, Fig. S3).
To determine whether the W140V point mutation had any effect on Myo19W140V enzymology, we measured the actin-dependent ATPase steady state activity of Myo19W140V-3IQ (Fig. 5). The steady-state parameters of Myo19W140V-3IQ revealed that, although the mutant’s k cat (ATP turnover per myosin head) was slightly faster than the WT’s, its K ATPase, (the actin concentration at the half maximum the ATPase activity), was significantly (5-fold) weaker (Table 1). This suggests that the intrinsic rate constants of Myo19W140V-3IQ ATPase cycle were significantly affected by the point mutation.
In this study, we present the kinetic analysis of Myo19W140V-3IQ within the framework of the WT Myo19-3IQ ATPase cycle that is fully characterized in the companion paper (Ušaj & Henn, parallel submission). The key rate and equilibrium constants that are compared between the WT and mutant protein are presented in Table 1. In the following sections, we present our analysis of ATP and ADP binding according to the nomenclature shown in the kinetic scheme, Fig. 6. The kinetic scheme in Fig. 6 shows the determined rate and equilibrium constants of both protein constructs. Myo19W140V-3IQ construct did not exhibit any FRET signal from the tryptophan(s) (λex, 280 nm) to the fluorescent mant moiety (λem,>400 nm) upon binding of mant nucleotides, whereas the WT exhibits such signal (Fig. S4). The loss of signal upon mant nucleotide binding limited the measurements of the mutant to direct measurements of nucleotide binding. However, this is quite interesting since it clearly shows that position 140 is very sensitive to nucleotide binding. Thus, it may have an effect on the interaction of Myo19 with nucleotide binding despite being situated between, and not within, the two-well-known nucleotide interacting motifs, P-loop and Loop I.
We measured ATP binding to Acto·Myo19W140V-3IQ by monitoring actin dissociation upon ATP binding. Similarly to WT Myo19-3IQ, but in contrast to many other studied myosins, Myo19W140V-3IQ binding to F-actin labeled with common actin-fluorescent probes (i.e., pyrene, coumarin, or IAEDANS) did not yield any signal. Therefore, we monitored ATP binding to Acto·Myo19W140V-3IQ by measuring changes in light scattering. ATP binding to Acto·Myo19W140V-3IQ induced dissociation from actin (Fig. 7a), exhibiting double-exponential kinetics similar to those observed for the WT (Figs 7b, S5). The fast phase is due to rapid ATP binding to readily accessible actomyosin conformation complex, while the slow phase results from the isomerization of the nucleotides binding pocket from the ‘close’ (AMC) state to the ‘open’ (AM°) state preceding ATP binding. K′α is the isomerization equilibrium constant described by the slow phase. The fast and slow observed rate constants (k obs,fast, k obs,slow) of Acto·Myo19W140V-3IQ exhibited hyperbolic dependence on [ATP] (Fig. 7c,d). Thus, the mechanism of ATP binding to Acto·Myo19W140V-3IQ is the same for both the mutant and the WT, with isomerization of Acto·Myo19W140V-3IQ preceding ATP. The ATP-binding mechanism is presented in more details in our co-submitted paper (Ušaj & Henn, parallel submission). Despite having a similar binding scheme, the rate and equilibrium constants for ATP binding to Acto·Myo19W140V-3IQ were strongly altered in comparison with those of WT-Myo19-3IQ (Table 1). Fitting the data of the mutant’s K′obs,fast yielded an equilibrium constant (1/K′1T) for ATP binding of 758±217µM, with a K′+2T (maximum rate of isomerization) (K′+2T≈K′diss) to the dissociated state of 381±72s−1. The fit of the curve for the mutant Myo19 is distinct from the origin (K′−2T=1.4±0.3), suggesting that ATP binding is slightly reversible compared with that observed for the WT (K′+2T» K′−2T ~0) (Fig. 7c). The k obs,fast hyperbolic dependence on [ATP] yields a second-order association rate constant for ATP binding of Acto·Myo19W140V-3IQ (K′1T K′+2T=0.5±0.2µM−1·s−1) that was approximately 4-fold smaller than that of the WT (Table 1). Overall, the rate constants describing ATP binding are significantly different between the WT and the mutant Myo19-3IQ constructs (Table 1), indicating that the W140V mutation strongly affects ATP-binding kinetics.
The slow phase observed upon ATP binding exhibited hyperbolic dependence on [ATP] (Fig. 7d). As with the WT, the slow phase arises from an isomerization between the “open” (AM°) and “closed” (AMC) nucleotide-binding states. Notably, the hyperbolic dependence of the slow phase was more pronounced for the mutant than for the WT. In this regard, the mutant resembles Myo1c, in which the slow phase has larger amplitudes and exhibits clear hyperbolic dependence on [ATP]. Overall, the maximal k obs, slow=K′+α for Myo19W140V-3IQ (29.3±4.3s−1) was 2-fold faster than that for the WT (Fig. 7d). The ratio of the amplitudes of the fast and slow phases (Afast/Aslow) defines the equilibrium constant K α for the AM° to AMC states42. To estimate the mutant’s K α, we averaged the Afast/Aslow ratios at the three highest [ATP] values tested, based on a previously reported approach44. The K α of the WT was 2-fold larger than that of the mutant (6.8±1.6 versus 2.9±1.3, respectively); hence the AMC state was more populated in the mutant than in the WT. As [ATP] increased, the observed transients of ATP binding to Acto·Myo19W140V-3IQ were dominated by the fast phase (Fig. 7e). This is because the product of k +1T×[AM°]×[T] is much greater than K′+α, resulting in the faster depletion of the AM° state at high [ATP] than its replenishment by isomerization of the AMC to AM° state, with a maximal rate of K′+α=29.3±4.3s−1 (Table 1). Finally, knowledge of K′α and K′+α permitted us to calculate K′−α, which was ~4-fold faster for the mutant (10.1±6.0s−1) than for the WT (Table 1).
The rate of ADP release and ADP’s affinity towards actomyosin can be determined by measuring the rate of ATP-induced dissociation of myosin from actin as a function of increasing prebound ADP. Such transients exhibit double exponentials, with the fast phase reporting ATP binding to the fraction of ADP-free actomyosin and the slow phase reporting the release of ADP. The amplitudes of the fast phase decrease as [ADP] increases until the fast phase is no longer observed, and actomyosin dissociation exhibits only a single (slow) phase. The observed rate constants for Acto·Myo19-3IQ-W140V ADP release were similar to those for WT Acto·Myo19-3IQ·ADP (Usaj&Henn), despite having a much lower signal. Figure 8a shows time courses of Acto·Myo19-3IQ-W140V·ADP upon mixing with excess ATP (125µM). The collection of data for the mutant was very challenging due to the weak signal; consequently, Fig. 8a shows only limited time courses. Nevertheless, we can extrapolate from the observed rate constants to obtain the slow-phase rate of k obs,slow≈2.9s−1, which differs only slightly from that of the WT (Fig. 8b). The affinity of ADP for Acto·Myo19-3IQ-W140V is determined from the hyperbolic dependence of the amplitudes of the slow phase (Aslow) versus [ADP]45, and can be fitted according to Equation 1:
The W140V mutation had a 10-fold weaker ADP affinity than the WT (1/K′D=1.6±0.3µM; Fig. 8c, Table 1), but nonetheless maintained a relatively strong affinity towards ADP. Myosin classes with high affinities towards ADP require much higher [ATP] for examination of the true ADP dissociation rate constant46. Therefore, we performed ATP-induced dissociation of Acto·Myo19W140V-3IQ in the presence of saturating [ADP]=25µM with increasing [ATP]. The k obs of the single exponential phase exhibited hyperbolic dependence on [ATP], saturating at K′−1D=12.7±1.6s−1, which is a bit faster than that of the WT (Fig. 8d, Table 1), but in accordance with a faster k cat observed for the mutant.
WT Acto·Myo19-3IQ·ADP dissociation at saturated [ADP] exhibited a well-resolved lag upon mixing with ATP (Ušaj & Henn, parallel submission). By contrast, Acto·Myo19W140V-3IQ·ADP exhibited a much shorter lag phase, even at lower temperatures. This suggests that the mutation has a mild effect on ADP binding and dissociation kinetics (k +1D, k −1D, Fig. 6), but a more pronounced effect on isomerization kinetics (k +2D, k −2D, Fig. 6).
Similarly to WT Myo19, the mutant’s ADP release was faster than the observed k cat but in the range to significantly contribute to the catalytic turnover of the ATPase cycle of Myo19W140V; hence the overall mechanism is conserved (Ušaj & Henn, parallel submission). Comparison of the WT and mutant nucleotide binding and isomerization steps revealed that at least two transitions limit the ATPase cycle under no external load (Fig. 6). Using the rates for ADP release (Fig. 6), the steady-state ATPase rate (k cat) for the mutant would be ~8.6s−1 head−1, twice as fast as our measured k cat (4.1s−1 head−1). Because a similar ratio was observed for WT (Ušaj & Henn, parallel submission), we argue in favor of the existence of an additional isomerization step (AMDC to AMD°) before ADP is completely released. Thus, the mutation inhibited Myo19 biological function by altering the intrinsic kinetic rate constants. The faster observed ADP transitions for Myo19W140V suggest lower efficiency of cargo translocation. This may be a consequence of the shorter dwelling time in the load-bearing ADP states; however, this proposal needs to be tested by force measurements and/or single-molecule experiments in the future.
The simulated ATPase cycle for WT permitted us to integrate the modified equilibrium and rate constants that we measured for Myo19W140V-3IQ, and to test whether they are sufficient to account for the observed weaker K ATPase and a slightly faster k cat for the mutant. We found that weakening of k +APi (defined as the second order rate constant for MADP·Pi binding to actin) by 6-fold compared with the WT resulted in K ATPase and k cat of 65.0±0.7µM (measured 56.7±10.9µM) and 6.5±0.02s−1 head−1 (measured 4.1±0.4s−1 head−1; calculated from the two slowest phases of ADP isomerization and its release ≈8.6s−1 head−1), respectively. The K DR (the [actin] at 0.5 duty ratio) for Myo19W140V-3IQ was 92.7±0.6µM, suggesting that the changes in the measured and estimated rates could account for the mutant motor’s ~7-fold weaker K DR (Fig. 9). Similarly, the duty ratio can be calculated according to De La Cruz et al.47 by the following parameters: k +APi, K H, and K −D (~0.15µM−1s−1, 1.07 and 8.6s−1, respectively), the K DR is 111µM. Myo19W140V-3IQ will require an actin concentration of ~400µM to reach a duty ratio of ~0.8 under saturating [ATP], concluding that Myo19W140V-3IQ is a very low duty ratio motor.
Our results show that the starvation-induced localization of Myo19 to filopodia is mediated by ROS, and that Myo19 can be recruited to filopodia tips by direct addition of H2O2 to the media. We also showed that mDia2ΔDAD, a factor that induces filopodia growth, could promote a similar localization of Myo19 to filopodia tips. By contrast, Myo10, which also promotes filopodia formation, failed to trigger Myo19 localization to filopodia. These observations suggest that Myo19 activation can be mediated by both intracellular and extracellular signals. Time-lapse microscopy revealed that Myo19 pauses during the formation and growth of filopodia prior to tracking their tips, and can also undergo retrograde movement coupled to filopodia retraction. Finally, we connected the enzymology of Myo19 to its filopodial activity by studying in detail the enzymology of a class-specific mutation in a highly-conserved residue, W140. We propose that modulation of the duty ratio of the mutant is partially responsible for its inability to perform the filopodia transport activity.
Glucose starvation increases cellular ROS levels and promotes a specific phosphorylation pattern involving focal adhesion proteins29. We hypothesized that the filopodial localization of Myo19 is mediated in a similar way. The inhibition of glucose starvation–induced filopodia localization by extracellular catalase and the synergistic effect of PG indicated that cells excrete ROS into the environment, which then enter the cell and exert their effects29. The mimicking of this phenotype by direct addition of H2O2 supports this mechanism. ROS participate in and affect signaling pathways by altering the activity of kinases or inhibiting phosphatases, and may directly modify Myo19 or alter the activity of an interacting protein48, 49. For example, EGF signaling involves cascades of kinase phosphorylation, and it is well established that EGF stimulation leads to filopodia formation50, 51. Thus, such signaling pathways may be also linked to Myo19 filopodial activity.
The majority of the motility we observed is directional and persistent towards the filopodia tips. This suggests a possible regulatory mechanism mediated either by affecting the motor’s ability to propel mitochondrial movement or by additional motility licensing factors that are regulated during filopodia elongation. The movement of Myo19 along the filament tracks appears to be deterministic, and is most likely carried out by an ensemble of motors. This is supported by the fact that we can observe motility of foci that must consist of more than a single fluorophore arising from tagged Myo19, which can be easily observed when they wrap around an organelle such as the mitochondria. Additional evidence supporting the idea that Myo19 acts through an ensemble of motors is provided by the observation that co-expression of WT Myo19 and the mutant Myo19W140V rescued the mutant phenotype, resulting in the localization of both WT and mutant motors at filopodia tips. This may arise from the presence of both WT and mutant motors on the same membrane, or from some as-yet-unknown interaction between several Myo19 motors.
Several other myosins localize to filopodia and affect their formation, likely by supplying the necessary components for filopodia formation. These myosins bring various cargos to the filopodia, such as actin bundling proteins, and proteins involved in cell adhesion such as β-integrins and VE-Cadherin52–56. Myo10 possesses unique forward and rearward movements within the filopodia of HeLa cells and promotes filopodia formation when ectopically expressed35, 57–59. Dictyostelium Myo7 fulfills a function similar to that of Myo10 within the filopodia56, 59, 60. Myo3 brings Espin1, an actin bundler, to filopodia and promotes their elongation54, 61. A barbed end–directed mutant of Myo6 can also promote filopodia formation, likely via the bundling activity of the dimer62. Myo15 localizes to filopodia tips, and Myo15-null mice have short stereocilia in comparison with WT animals63.
Myo19 and Myo19W140V pose a great challenge for in vitro characterization, mainly because of their lack of quenching of fluorescent actin, as well as their silent fluorescence upon ATP and ADP binding. However, the clear cellular phenotype allows us to test how the motor enzymology supports the protein’s activity. Most importantly, we were able to identify key parameters in the ATPase cycle of Myo19W140V that may account for the observed phenotype. First, the W140V mutant had a slightly higher k cat but a ~5-fold weaker K ATPase in comparison with the WT. Second, both the k α transition (AMC↔AM°), and the k αD transition (AMDC↔AMD°) were faster in the mutant. Third, the lag preceding ADP release was much more prominent in the WT than in the mutant. Overall, these changes indicate that, although the mechanism remained similar, the rate and hence the equilibrium constants of the biochemical transitions became faster in the mutant. Because we established a strong link between the K ATPase and the k +APi of Myo19, we can model the entire ATPase cycle of Myo19W140V-3IQ. Despite not having all the cycle parameters at hand, we predict that the K DR of the mutant was reduced by ~7-fold compared with the WT. It is reasonable to speculate that a multi-motor ensemble of Myo19 on the mitochondria would compensate for the reduction in the duty ratio of the motors, especially in a dense actin mesh. However, it is also reasonable to assume that unique load-dependence transitions are also important for this motor function; hence the weakened ADP affinity may also affect the force-sustained states. Previous authors have proposed that myosins that undergo the K α transition are highly force-sensitive64. Studies performed on Myo1c on the effect of R156W mutation which has been shown to increases K α 65 may suggest some parallel link to explain how W140V is hindering Myo19 transport of mitochondria. The frictional loading assays indicated that load-induced inhibition of actin gliding is less pronounced in the R156W Myo1c mutant65. It has been proposed that an increase in K α (i.e., a stabilization of the nucleotide-accessible AM° state) might be an indication of reduced tension sensitivity65. If this hypothesis is true, then in our case the W140V mutation (decreases K α) may produce a Myo19 motor with higher sensitivity to load. Such motors would translocate cargo less efficiently. Future work should seek to test this hypothesis.
The physiological role of mitochondria in filopodia formation and elongation remains unknown. Recent studies show that mitochondria migrate into damaged axons, and that a high density of mitochondria in a damaged axon is correlated with improved regeneration. When an axon regenerates, it forms a growth cone, which is an actin-based structure whose size is correlated with the density of mitochondria within it. Collectively, these studies reveal that the major role of the mitochondria in the damaged axon is to produce ATP that is essential for repair, possibly by replenishing the local pool of G-Actin–ATP, which is essential for filopodia elongation by formins66–68. In addition, axon branching occurs preferentially near stalled mitochondria, and the latter act as hot spots of protein synthesis69. These remarkable findings suggest that mitochondrial motility and localization is strongly coupled to ATP production and actin dynamics. We propose that Myo19 plays a similar role, increasing mitochondrial density to supply ATP at filopodia tips.
All chemicals and reagents were of the highest purity commercially available. ATP was purchased from Roche Applied Science (Penzberg, Germany) and ADP was purchased from Bio Basic (Markham, ON, Canada). Nucleotide concentrations were determined by absorbance at 259nm using ε259 of 15,400M−1cm−1. In all experiments one molar equivalent of MgCl2 was added to nucleotide solutions immediately before use.
3-(N-Morpholino)propanesulfonic acid (MOPS), Ethylene glycol-bis(2-aminoethylether)-N,N,N’,N’-tetraacetic acid (EGTA), Apyrase (potato grade VII), phalloidin were purchased from Sigma-Aldrich (St. Louis, Missouri, USA). Alternatively, the same form of phalloidin was purchased also from Setareh Biotech (Eugene, OR, USA). MgCl26H2O came from Bio Basic (Markham, ON, Canada) and KCl from Merk (Darmstadt, Germany).
Emerald tagged Myo19 (emMyo19) and Ruby tagged Myo19 (Myo19-Ruby) were generated by swapping the eGFP from Myo19-eGFP20 with Emerald or Ruby from mEmerald-Lifeact-7 or mRuby-LifeAct, respectively, using the restriction enzymes AgeI and NotI. emMyo19, full length Myo19-Halo, or Myo19-3IQ-Halo20 (Usaj and Henn, parallel submission) were used as templets for generation of Myo19 mutants and Myo19-3IQ mutant constructs (Table S1). For expression and purification, Myo19W140V-3IQ was cloned into the pF14A Flexi Vector using the restriction enzymes NheI and PspXI (Promega). Human calmodulin (Myo19 light chain according to18) was cloned to pF4A CMV Flexi Vector (Promega) according to the manufacturer’s protocol. mEmerald-Lifeact-7 was a gift from Michael Davidson (Addgene, Cambridge, MA, USA, plasmid # 54148). Plasmid expressing Myo10 was a kind gift of Prof. Alexander Bershadsky (MBI, NUS). eGFP-tagged full-length mDia2/Drf3 was a kind gift from Drs. Klemens Rottner and Jan Faix (TU Braunschweig and HMS, Germany, respectively). mRuby-LifeAct was a kind gift from Prof. Arie Admon (Technion - IIT, IL). Primers pairs harboring a mutation (Table S1) were phosphorylated using T4 PNK (NEB, Ipswich, MA, USA) according to the manufacturer’s instructions. The primers were used to amplify the template using Phusion High fidelity DNA polymerase (Thermo Fisher Scientific, Waltham, MA, USA). Ligations using T4 ligase (Thermo Fisher Scientific, Waltham, MA, USA) were carried out according to the manufacturer’s instructions.
All materials were purchased from Sigma unless otherwise indicated. Fetal calf serum, L-Glutamine, HEPES-KOH pH 7.4, penicillin, streptomycin, amphotericin B were purchased from Biological Industries, Beit Haemek, Israel. Paraformaldehyde (PFA) was purchased from EMS (Hatfield, PA, USA).
U2OS cells were grown at 37°C and 5% CO2 in DMEM supplemented with 10% Fetal calf serum (FCS, Biological industries, Beit Haemk, Israel), 4.5g/L glucose, 2 mM L-Glutamine, 20mM HEPES-KOH pH 7.4, 100 U/ml penicillin, 100µg/ml streptomycin and 0.25µg/ml amphotericin B. Starvation conditions: Cells were rinsed twice with PBS and incubated in starvation medium (Glucose free DMEM supplemented with 20mM HEPES-KOH pH 7.4, 5mg/ml BSA, 100 U/ml penicillin, 100µg/ml streptomycin and 0.25µg/ml amphotericin B) for two hours. Antioxidants were prepared according to the manufacturer’s directions and diluted in media prior to being added to the cells.
Cells grown on coverslips were fixed with 4% PFA in CBS (10mM MES pH 6.1, 138mM KCl, 3mM MgCl2, 2mM EGTA, 7.8% sucrose) for 10min. The PFA was removed by washing with 0.3% TBST (0.3% Tween-20 in Tris buffered saline), and the cells were blocked for one hr with 10% BSA in TBST. Primary antibody was added at a dilution of 1:250 in TBST for one hr (Myo19 -Abcam ab174286, ATP5a – Abcam ab14748). The coverslips were washed three times for five mins in TBST and a secondary antibody was added at a dilution of 1:250 in TBST for one hour (Alexa-fluor 488 conjugated anti rabbit (Life Technologies) and Alexa fluor 633 conjugated anti mouse (Life technologies). The coverslips were washed three times for five mins in TBST (once in the presence of Hoechst) and mounted on slides using Fluoromount-G. Actin in fixed samples was visualized by incubating the cells with either Phalloidin 633 or Phalloiding 546 (Life technologies).
Transfections were performed using Polyethylenimine (PEI) purchased from Weis Scientific – PolySciences (Warrington, PA, USA). Adherent U2OS cells were plated a day before transfection on glass bottom dishes purchased from In-vitro Scientific (Mountain View, CA, USA) and allowed to adhere overnight. Plasmid DNA and PEI were diluted separately in 150mM NaCl, combined and complex formation was allowed for 25min at RT before addition to the cells and incubation overnight. Hoechst 33342 (0.75µg/ml) was added 15min prior to imaging. To visualize Halo tagged Myo19, HaloTag TMR ligand (Promega) was added to cells transfected with Halo tagged Myo19 at a final concentration of 50nM and incubated overnight. Images were taken using LSM710 at ×40 magnification or INCELL2000 at ×60 for live cells in an environmental chamber or using LSM710 at ×63 for fixed samples.
The images that were acquired using the INCELL2000 were 2048×2048 pixels in size at 100 DPI, 16bit in bit depth and were acquired at 37°C in an environmental chamber in DMEM (See starvation conditions section). Manufacturer’s settings were used for the excitation and emission. The images that were acquired using the LSM710 had a scaling of ~0.076µm (i.e 1 pixel=0.076µm), 8bit in bit depth and were acquired at 37°C in an environmental chamber in DMEM (See starvation conditions section). The excitation wavelengths used: 405, 488, 543 and 639nm. We used GaAS detectors for detection of the green and red signal using 525/50 and 594/46 filter sets. Blue emission was detected between 410—500nm and far-red emission was detected between 644–797nm.
All of the images that were acquired by the dedicated instrument’s software were processed in FIJI. Images were adjusted for brightness and smoothed. Filopodia length was measured manually in FIJI using the “freehand line”.
Expression of constructs were performed in suspension adapted HEK293SF-3F670 (ATCC) cells grown in serum free media EX-CELL (Sigma-Aldrich) or in house made proprietary cell culture media. Myo19W140V-3IQ mutant was transiently co-transfected with calmodulin (ratio used 2:1). The transfected cells (initial concentration of 1–2×106cells/ml) were allowed to grow for 2–3 days at 37°C while shaking (reaching 2–4×106cells/ml) and were then pelleted at 350×g for 10min. All protein purification steps were carried out on 4°C. On the day of purification cell pellets were re-suspended in HaloTag system compatible lysis buffer (20mM MOPS, 130mM KCl, 5mM MgCl2, 2mM EGTA, 2mM ATP, 0.5mM DTT, 5% Glycerol, 0.05% NP-40, 2µM calmodulin, 1 X Promega protease inhibitor cocktail, pH 7.5) to achieve cell density≈40×106cells/ml and lysed by homogenization in a Dounce homogenizer with a glass pestle by ≈60 strokes. Lysates were supplemented with additional 1mM MgATP and cleared by ultra-centrifugation at 100,000×g for 1h (4°C). From this point on, the purification followed HaloTag purification system manuals and protocols provided by Promega. Final elutions were concentrated using Millipore (Merck, Darmstadt, Germany) 3kDa amicon ultra centrifugal filter unit and dialyzed into storage buffer (20mM MOPS, 75mM KCl, 5mM MgCl2, 2mM EGTA, 0.5mM DTT, 50% Glycerol, pH 7.5). The purified constructs were stored at −20°C. Myo19-3IQ and Myo19-3IQW140V concentration were determined using predicted extinction coefficient at 280nm (λex.coff, WT=105,770M−1cm−1, λex.coff, W140V=100207M−1 cm−1, ExPASy ProtParam). The fraction of [active heads] was obtained as for the WT assuming similar efficiency of the expression and purification system to produce WT and mutant constructs (Usaj and Henn, parallel submission).
Actin was purified from rabbit or chicken skeletal muscle labeled and gel-filtered over Sephacryl S-300 HR71. Ca2+-actin monomers were converted to Mg2+- actin monomers with 0.2mM EGTA and 40µM MgCl2 (excess over [actin]) immediately prior to polymerization by dialysis against KMg50 buffer (20mM MOPS, 50mM KCl, 2mM MgCl2, 0.2mM EGTA, 1mM DTT, pH 7.3 at 25°C). Phalloidin (1:1 molar ratio) was used to stabilize actin filaments. 7-Diethylamino-3-((((2-maleimidyl) ethyl)amino)carbonyl)coumarin-labeled phosphate-binding protein (PiBP) was expressed, purified, and labeled as described72. Calmodulin was expressed and purified in bacteria as described73.
Steady-state ATPase Activity - The actin-activated steady-state ATPase activity of Myo19-3IQW140V was measured at 25±0.1°C in KMg50 buffer supplemented with 2mM MgATP using the NADH coupled assay46 by monitoring changes in absorption at 340nm. The Myo19-3IQW140V concentration was 10–50nM but constant during each individual experiment.
All experiments were performed in KMg50 buffer with the Hi-Tech Scientific SF-61DX2 stopped-flow apparatus (TgK Scientific Limited, Bradford-on-Avon, UK) at 25±0.1°C. The concentrations stated through all the text are always final concentrations after mixing (i.e. in the observation cell) unless noted otherwise. Light scattering was measured at 90 ° with excitation at 313nm. Most time courses shown are of individual, 2000-point transients collected with the instrument in oversampling mode, where the intrinsic time constant for data acquisition is ≈64 µs. Typically, multiple (3 to 5) time courses were averaged before analysis. Time courses displayed fast and slow phases were collected on a logarithmic or split time scale. Using Kinetic Studio software provided with the instrument or with Origin (OriginLab Corporation, Northampton, Massachusetts, USA), time courses of signal (fluorescence, light scattering) change were fitted to a sum of exponentials (Equation 2),
where F(t) is the signal at time t, F∞ is the final signal value, A i is the amplitude, k i is the observed rate constant characterizing the i-th relaxation process, and n is the total number of observed relaxations. The value of n was either one (single exponential) or two (double exponential). For mant-nucleotides binding photobleaching was not affecting time courses inside fitting windows (<5s). The dead time of the instrument determined from the reduction of 2,6-dichlorophenolindophenol with ascorbic acid in absorbance mode was 1 ms. Fitting was limited to data beyond 1 ms to account for the instrument dead time and to exclude data acquired during the continuous flow phase of mixing, as recommended by the manufacturer.
Uncertainties are reported as standard errors in the fits unless stated otherwise and were propagated using the general formula (Equation 3),
where the experimental measurements x1, x2… xn have uncertainties dx1, dx2… dxn and a is a function of x1, x2… xn.
Time courses of nucleotide binding were acquired under pseudo first-order conditions with [nucleotide]>>[myosin or actomyosin]. Actomyosin samples were prepared by mixing equal molar amount of Myo19-3IQ with actin filaments or, where specified, with [actin]» [myosin]. Myosin and actomyosin samples were treated with apyrase (0.01 unit/ml) and equilibrated on ice for 10min before measurements. The final apyrase concentration after mixing was 0.005 unit/ml used to deplete ATP and ADP from Myo19-3IQ and ActoMyo19-3IQ46.
Kinetic simulations and modeling of reaction time courses were performed using KinTek Explorer74, 75. Simulations of the steady-state reaction time courses were performed according rate constants in Table 1. For the mutant, all the rate constants parameters that were measured were incorporated into the model. K′−APi, was the only free parameter that was allowed to be determined by the simulated model until the curve converge to the observed [Actin]-dependence steady-state ATPase. To plot the duty ratio as a function of actin, the fraction of bound and unbound myosin states to actin were directly summed based on the whole reaction simulation mechanism.
Statistical significance between mock and treated groups was performed using two tailed Student’s T-test. We assumed equal variance since the samples originated from the same source. Statistical significance was scored as P≤0.05 *P≤0.01 **P≤0.001***. All presented error bars indicate standard deviation.
No datasets were generated or analysed during the current study.
We are grateful for Dr. Nitzan Dahan’s assistance with the fluorescence microscopy imaging from Life Sciences and Engineering Infrastructure Unit at the Technion. AH acknowledges Marie Curie Career Integration Grant [grant number 1403705/11] and the Israeli Scientific Foundation [grant No. 296/13]. MU acknowledges support from the Lady Davis Post-Doctoral Fellowship and Israeli Government MFA Scholarship program. Special thanks to Mrs. Hadar Levi for her technical assistance. We would like to thank the group members for their helpful comments on the manuscript.
B.S., M.U., and A.H. designed and performed experiments, analyzed data, and wrote the paper; N.W.M. and R.R. performed experiments.
The authors declare that they have no competing interests.
Electronic supplementary material
Supplementary information accompanies this paper at doi:10.1038/s41598-017-11002-9
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.