ERK3/MAPK6 dictates CDC42/RAC1 activity and ARP2/3-dependent actin polymerization

The actin cytoskeleton is tightly controlled by RhoGTPases, actin binding-proteins and nucleation-promoting factors to perform fundamental cellular functions. We have previously shown that ERK3, an atypical MAPK, controls IL-8 production and chemotaxis (Bogueka et al., 2020). Here, we show in human cells that ERK3 directly acts as a guanine nucleotide exchange factor for CDC42 and phosphorylates the ARP3 subunit of the ARP2/3 complex at S418 to promote filopodia formation and actin polymerization, respectively. Consistently, depletion of ERK3 prevented both basal and EGF-dependent RAC1 and CDC42 activation, maintenance of F-actin content, filopodia formation, and epithelial cell migration. Further, ERK3 protein bound directly to the purified ARP2/3 complex and augmented polymerization of actin in vitro. ERK3 kinase activity was required for the formation of actin-rich protrusions in mammalian cells. These findings unveil a fundamentally unique pathway employed by cells to control actin-dependent cellular functions.


Introduction
Actin is one of the most abundant and highly conserved proteins with over 95% homology among all isoforms (Vandekerckhove and Weber, 1979). In cells, it is present in globular/monomeric form (G-actin), which can polymerize into branched and elongated filamentous forms (F-actin) in a dynamic, spatially, and temporally controlled process (Skruber et al., 2018;Lee and Dominguez, 2010;Carlsson, 2006;Carlier et al., 1997). The cytoplasmic actin network constitutes an important part of the cytoskeleton, which not only mechanically supports the plasma membrane and gives the cell its shape, but fulfills a variety of other functions: It regulates velocity and directionality of cell migration, enables intracellular signaling and transport, supports cell division, and more (Yamaguchi and Condeelis, 2007;dos Remedios et al., 2003). Moreover, by branching and bundling of F-actin filaments the actin cytoskeleton can form unique structures such as lamellipodia and filopodia, which in epithelial cells regulate cell polarization and are crucial for the contact with their environment, while in other cell types these actin rich structures control motility, chemotaxis, and haptotaxis (Vasioukhin et al., 2000;Raich et al., 1999;Millard and Martin, 2008;Khurana and George, 2011;Bahri et al., 2010).
The small RhoGTPases belong to the Ras homologous (Rho) superfamily of GTP-binding proteins (Hodge and Ridley, 2016). They function as binary molecular switches circulating between active (GTP-bound) and inactive (GDP-bound) state. The on-off regulation of the RhoGTPases is under the control of guanine exchange factors (GEFs) and GTPase-activating proteins (GAPs) (Hodge and Ridley, 2016). GTP-bound active CDC42 and RAC1 promote the formation of filopodia and lamellipodia, respectively, by activating the NPFs, WASP, and WAVE (Rohatgi et al., 1999;Carlier et al., 1999). Recent studies revealed a role for kinases in directly activating RhoGTPases mainly through phosphorylation, but the identity of the relevant kinase(s) remained elusive (Tong et al., 2013;Forget et al., 2002;Chang et al., 2011;Tu et al., 2003;Ren et al., 2001). ERK3 (MAPK6) is an atypical member of the mitogen-activated protein kinases (MAPKs). Its physiological roles are tissue-specific and can be both dependent and independent of its catalytic activity: kinase-dependent roles of ERK3 include promoting lung cancer cell migration and invasion, while its function in the regulation of breast cancer cell morphology and migration is shown to be partially kinase-independent (Al-Mahdi et al., 2015;Bogucka et al., 2021;Bogucka et al., 2020;Elkhadragy et al., 2020). Moreover, ERK3 is a labile protein whose rapid turnover has been implicated in cellular differentiation (Bogucka et al., 2021;Coulombe et al., 2003). Although the specific signaling mechanisms involved in the regulation of ERK3 have been studied by us and others and provided tissuespecific insights into this atypical MAPK function (Bogucka et al., 2021;Elkhadragy et al., 2020;El Merahbi et al., 2020;Elkhadragy et al., 2017;Long et al., 2012;Sauma and Friedman, 1996), many aspects of ERK3 biology such as its role in suppressing melanoma cell growth and invasiveness remain elusive (Chen et al., 2019). It is, however, known that ERK3 possesses a single phosphorylation site at serine 189 (S189) within S-E-G activation motif, located in its kinase domain, which was shown to be phosphorylated by group I p-21-activated kinases (PAKs), the direct effectors of RAC1 and CDC42 GTPases (Déléris et al., 2011;De la Mota-Peynado et al., 2011;Coulombe and Meloche, 2007).
ERK3 signaling is required for the motility and migration of different cancer cells, yet its direct role in the regulation of the polarized phenotype of the cells and actin cytoskeleton is lacking. Here, we unveil a multilayered role for ERK3 in regulating actin filament assembly. We demonstrate that ERK3 controls bundling of actin filaments into filopodia via activation of CDC42. Mechanistically, ERK3 directly binds to RAC1 and CDC42 in a nucleotide-independent manner and is required for the activity of both RhoGTPases in vivo. Furthermore, ERK3 could function as a direct GEF for CDC42 in vitro. In addition, ERK3 stimulates ARP2/3-dependent actin polymerization and maintains F-actin abundance in cells via direct binding and phosphorylation of the ARP3 subunit at S418. While the kinase activity is not required for CDC42 activation, RAC1 activity partially depends on ERK3 kinase function. The in vivo relevance of these ERK3 functions was corroborated by analyses of the motility of tumor cells in orthotopic mammary tumors grown in mice. Together, our results establish a hitherto unknown regulatory pathway directly controlled by ERK3 involving the major signaling molecules CDC42, RAC1, and ARP2/3 in the modulation of the actin cytoskeleton and cell migration.

Results
ERK3 is required for motility of mammary epithelial cells ERK3 has been shown to regulate mammary epithelial cancer cell migration and metastasis (Al-Mahdi et al., 2015;Bogucka et al., 2020). To uncover the molecular underpinnings of the role of ERK3 in cell migration, we investigated the intracellular localization of ERK3 by immunocytochemical analyses using a validated antibody (Bogucka et al., 2021;Bogucka et al., 2020). Endogenous ERK3 co-localized to the F-actin-rich protrusions in human mammary epithelial cells (HMECs) ( Figure 1A). Reorganization of actin protrusions at the leading edge of the cells initiates cell migration (Yamaguchi and Condeelis, 2007;Affolter and Weijer, 2005); therefore, we further analyzed the significance of ERK3 in cancer cell motility by loss-of-function studies. Depletion of ERK3 reduced speed and displacement length (distance between the start point and the end position of the cell along each axis) of MDA-MB231 cells . As shown in Figure 1C, depletion of ERK3 strongly reduced random cell motility, which remained low over time, as indicated by the calculated acceleration ( Figure 1D) and concomitantly resulted in a shorter overall displacement length ( Figure 1E) with shorter migration track length and average speed (Figure 1-figure supplement 1A and B). Considering that responses of tumor cells depend on their microenvironment, we tested the in vivo motility of these cells. Intravital imaging of the orthotopic mammary tumors confirmed reduced motility of the ERK3 knockdown MDA-MB231-GFP cells ( Figure 1G and Videos 5 and 6). Together, these data suggest that ERK3 likely controls actin cytoskeleton dynamics, thereby influencing cell shape, motility, and polarized migration.

ERK3 is activated and required for EGF-mediated directional migration in epithelial cells
To further elucidate the physiological role of ERK3 in breast epithelial cell motility, we assessed EGFinduced chemotactic responses to evaluate directional migratory properties of control and ERK3depleted primary HMECs and metastatic MDA-MB231 cells. Interestingly, knockdown of ERK3 significantly attenuated the EGF-induced chemotaxis of both cell types (Figure 1-figure supplement 1C-F, respectively).
Considering that ERK3 participated in EGF-mediated chemotactic responses of both primary and oncogenic mammary epithelial cells, we tested ERK3 protein kinetics in response to the growth factor. EGF induced serine 189 (S189) phosphorylation of ERK3 at early time points in primary epithelial cells (Figure 1-figure supplement 2A and B) and MDA-MB231 cells (Figure 1-figure supple ment 2C and D). A slight increase in ERK3 protein levels was detected in HMECs upon EGF treatment ( Figure 1-figure supplement 2A and B). Further cycloheximide (CHX) chase experiments coupled with RT-PCR analyses suggested that EGF treatment did not increase ERK3 protein half-life ( Figure 1-figure supplement 2E-2H). Taken together, these data confirmed that ERK3 was activated in response to EGF and required for EGF-mediated directional migration of epithelial cells. We then investigated the significance of S189 phosphorylation of ERK3 on F-actin and formation Figure 1. ERK3 localizes to the F-actin-rich protrusions and regulates mammary epithelial cells motility. (A) Confocal analyses of control (shCo) and ERK3 knockdown (shERK3) human mammary epithelial cells (HMECs) co-stained with anti-F-actin (green) and anti-ERK3 antibodies (red). Hoechst staining (blue) was used to visualize cell nuclei. Scale bars 28 µm. Higher magnification images of the boxed areas are shown on the right. Red and white arrows exemplify lamellipodia and filopodia, respectively. (B-F) Single-cell tracking analyses of control (shCo) and ERK3 knockdown (shERK3) MDA-MB231- Figure 1 continued on next page GFP cells. (B) Representative images of the cell migration assay were taken at the start (top panel) and end point (20 hr) (lower panel) of the full track length. Random fields were selected, and boxed areas were magnified on the right. Arrows indicate exemplified cell tracking at the beginning and at the end of the tracking. (C-E) Violin plots present cell distribution according to the analyzed motility parameters, calculated as described in the 'Materials and methods' section. Plotted are median (solid line) and 25%/75% quartiles (dashed lines). (C) Speed (µm/s), (D) acceleration (µm/s 2 ), and (E) displacement length (µm) of 14,659 single cells (n = 14659) are depicted. Significance was determined using nonparametric Mann-Whitney test. *p<0.0332, **p<0.0021, ***p<0.0002, ****p<0.0001. Analyses of the migration tracks length and average speed are depicted in Figure 1-figure supplement 1A and B. (F) Cell migration patterns of four randomly selected control and ERK3 knockdown cells were plotted as x-y trajectories over the 20 hr tracking. (G) Quantification of the intravital imaging of the control (shCo) and ERK3 knockdown (shERK3) orthotopic mammary tumors (MDA-MB231-GFP). The number of GFP-positive motile cells was quantified as described in the 'Materials and methods' section and is presented as mean ± SEM from seven animals (n = 7) per condition. Please see representative Videos 5 and 6. Chemotactic responses of the ERK3-depleted cells were assessed in the presence of EGF and are presented in Figure 1-figure supplement 1C-F. ERK3 protein kinetics and stability in response to EGF are depicted in Figure 1-figure supplement 2. (H-J) Depletion of ERK3 alters leading edge protrusions in primary mammary epithelial cells. (H) Confocal F-actin (green)/Hoechst (blue) and SEM images analyses of control and ERK3 knockdown HMECs are presented after 4 hr starvation. For the confocal analyses, exemplified filopodia fitting was included in the middle panel. (I) Western blot analyses presenting ERK3 knockdown efficiency and total levels of β-actin. Ponceau S was used as a loading control. (J) Filopodia number was quantified as described in the 'Materials and methods' section and is presented as mean ± SEM of 17 single cells (n = 17); *p<0.0332, **p<0.0021, ***p<0.0002, ****p<0.0001, unpaired t-test. Effect of S189 phosphorylation of ERK3 on F-actin assembly is presented in Figure 1-figure supplement 3.
The online version of this article includes the following source data and figure supplement(s) for figure 1: Source data 1. Prism and Excel file for Figure 1C-E.
Source data 2. Excel file for Figure 1F.
Source data 3. Prism and Excel file for Figure 1G.
Source data 4. Full membrane scans for western blot images for Figure 1I.
Source data 5. Prism and Excel file for Figure 1J.             https://elifesciences.org/articles/85167/figures#video2 of actin-rich protrusions. ERK3-depleted HMECs cells (shERK3 -3′ UTR) presented in Figure 1A were transfected with V5-tagged ERK3 wild type (ERK3 WT) or S189 mutants (ERK3 S189A/ERK3 S189D) and expression of exogenous ERK3 was visualized by V5-tag staining. Co-staining of the cells with phalloidin green did not show any significant differences in F-actin organization upon the overexpression of the S189 phosphorylation mutants compared to the wild type ERK3 transfected cells (Figure 1-figure supplement 3). Moreover, we observed that exogenous ERK3 localized predominantly in the cytosol compared to the endogenous ERK3, which we detected at the cell edges and F-actin rich protrusions ( Figure 1A). These results prompted us to prefer ERK3 knockdown cells for the study of ERK3-mediated phenotypes.

ERK3 is required for actin-rich protrusions
In order to migrate, cells acquire morphological asymmetry to drive their locomotion. One of the first steps into breaking spatial symmetry of the cell involves its polarization and formation of membrane extensions (Theriot and Mitchison, 1991). In response to extracellular guidance, cells polarize and extend F-actin-rich protrusions at the leading edge to direct cell migration (Affolter and Weijer, 2005;Ridley et al., 2003;Lauffenburger and Horwitz, 1996). As shown in Figure 1A, ERK3 co-localized with F-actin at the edge of the cell and at the protruding filopodial spikes. Moreover, knockdown of ERK3 clearly reduced F-actin content and actin-rich protrusions in human mammary primary Video 3. Movies representing motility of control (Videos 1 and 2) and ERK3-knockdown (Videos 3 and 4) MDA-MB231 cells. https://elifesciences.org/articles/85167/figures#video3 Video 4. Movies representing motility of control (Videos 1 and 2) and ERK3-knockdown (Videos 3 and 4) MDA-MB231 cells. https://elifesciences.org/articles/85167/figures#video4 Video 5. Representative videos following the motility of the control (shCo) and ERK3 knockdown (shERK3) GFP-expressing MDA-MB231 cells. https://elifesciences.org/articles/85167/figures#video5 Video 6. Representative videos following the motility of the control (shCo) and ERK3 knockdown (shERK3) GFP-expressing MDA-MB231 cells. epithelial cells ( Figure 1A). Further analyses of control and ERK3 knockdown HMECs at the singlecell level confirmed that ERK3-depleted cells had limited, diffuse F-actin, concentrated in the cortical areas of the cells, and a significantly decreased number of filopodia ( Figure 1G-J).

ERK3-dependent regulation of RAC1 and CDC42 activity
RhoGTPases are the major membrane signaling transmitters that drive polarized cell asymmetry by locally regulating the formation of F-actin-rich protrusions (Hall, 1998;Nobes and Hall, 1995). EGF signaling triggers activation of RAC1 and CDC42, which link the signal to the actin cytoskeleton leading to the formation of lamellipodia-and filopodia-like protrusions, respectively (Figure 2A; Kurokawa et al., 2004).
Considering that loss of ERK3 significantly affected F-actin distribution and filopodia formation in primary mammary epithelial cells ( Figure 1A and G-J), we tested the activity of the two major regulators of actin assembly, CDC42 and RAC1 ( Figure 2B-G). ERK3 knockdown significantly decreased the levels of both basal and EGF-induced GTP-bound CDC42 and RAC1 in primary (HMEC) ( Figure 2B, D, and F) and triple-negative breast cancer (MDA-MB231) mammary epithelial cells ( Figure 2C, E, and G), respectively. HMECs were cultured in growth factor/hormone-enriched medium that also contained EGF. Therefore, we withdrew the supplements to test the effect of EGF alone on the activity of RAC1 and CDC42 in these cells and compared the responses with those of the cancer cells. Interestingly, we detected high RAC1 and CDC42 levels at steady state in HMEC cells, which underwent a rapid activation-inactivation cycle after starvation and restimulation with EGF alone. The inhibitory effect of ERK3 knockdown on RAC1 was not as striking as with CDC42 ( Figure 2D-G). Intriguingly, we could observe co-precipitation of ERK3 protein with the active RhoGTPases at endogenous levels, suggesting that these proteins exist in a multimeric complex ( Figure 2B and C). We expanded these analyses to different cell types and consistently found that depletion of ERK3 reduced the activity of RAC1 and CDC42 and that the active RhoGTPases specifically co-precipitated endogenous ERK3 (Figure 2-figure supplements 1 and 2). Interestingly, the effect of ERK3 on RAC1 and CDC42 activity is either EGF-dependent (MCF7 cells [ Figure

ERK3 directly binds to RAC1 and CDC42
In light of these observations, we tested whether any direct interaction existed between these RhoGTPases and ERK3. By employing purified full-length recombinant proteins, we found that ERK3 directly bound to RAC1 and CDC42 in a nucleotide-independent manner ( Figure 2H and I). The nucleotide-loading status of the purified RAC1 and CDC42 proteins was simultaneously assessed by GST-PAK1-PBD pull-down assay followed by immunoblots ( Figure 2-figure supplement 3). To further assess the affinity of the binding, we performed concentration-dependent protein binding assays with ELISA as described in the 'Materials and methods' section. Binding of full-length ERK3 to RAC1 and CDC42 could be detected at a low nanomolar range (5 nM) ( Figure 2J and K). Next, we tested whether the kinase domain of ERK3 (amino acids [aa] 9-327) would be sufficient for this interaction. In vitro GST-pull-down assays confirmed that the kinase domain of ERK3 was sufficient for its binding to RAC1 and CDC42 ( Figure 2L and

ERK3 functions as a GEF for CDC42
The intrinsic GDP-GTP exchange of RhoGTPases is a slow process, which is stimulated in vivo by GEFs ( Figure 3A; Hodge and Ridley, 2016). To examine the ability of the ERK3 kinase domain to stimulate GDP-GTP exchange, we incubated GDP-loaded RhoGTPases with non-hydrolysable GTPγS in the presence or absence of the ERK3 kinase domain and quantified the final amount of GTP-bound CDC42 ( Figure 3B and C) and RAC1 ( Figure 3D and E) using GST-PAK1-PBD fusion beads. The ERK3 kinase domain failed to interact directly with the GST-PAK1-PBD protein but bound to active CDC42 in the same assay confirming that conformationally stable protein was employed in these assays (Figure 2figure supplement 4). These experiments showed that although ERK3 regulated the activity of both   Figure 3B and C). We further corroborated these results by measuring the GEF activity of the ERK3 kinase domain in vitro by employing a second fluorophore-based assay. The results confirmed those of the previous assay as addition of ERK3 stimulated GTP-loading on CDC42 ( Figure 3F) but not on RAC1 ( Figure 3G). Notably, ERK3 exerted a more potent and long-lasting effect than Dbl's Big Sister (Dbs), a well-established GEF for CDC42 and RhoA proteins ( Figure 3F; Baumeister et al., 2006).
To assess whether full-length ERK3 could potentiate the effect in terms of activation of CDC42, we compared GEF activity of the ERK3 kinase domain with that of the full-length protein. Interestingly, full-length ERK3 exhibited weaker GEF activity in vitro than the kinase domain of ERK3 (Figure 3figure supplement 1A). The observed in vitro discrepancy between the kinase domain and the fulllength protein could be attributed to differences in posttranslational modifications due to different production methods of the two recombinant proteins, which might affect conformation and/or activity. The full-length protein employed in these assays was purified from Sf9 cells, while the kinase domain was purified from bacterial cells and therefore lacked -among other things -Ser189 phosphorylation and thus kinase activity in vitro ( Figure 3-figure supplement 1B). We next examined whether ERK3 is required for the localization and activity of CDC42 at the plasma membrane (PM). Cellular fractionation assays revealed that knockdown of ERK3 does not disrupt the localization of CDC42 and RAC1 to the PM ( Figure 4A-C). On the contrary, the PM fraction from ERK3-depleted HMECs had more ERK3 with active RAC1 and CDC42. Relative levels of (D, E) active CDC42 and (F, G) RAC1 were calculated with respect to the total protein levels and are presented as mean ± SEM from minimum three (n = 3) independent experiments; *p<0.0332, **p<0.0021, ***p<0.0002, ****p<0.0001, one-way ANOVA, Tukey's post-test. ERK3-dependent regulation of the RAC1 and CDC42 activity was assessed in multiple cell types and data are presented in Figure 2-figure supplements 1 and 2. (H, I) In vitro interaction between ERK3 and GDP/GTP bound (H) RAC1 and (I) CDC42 was assessed, GST was used as a negative control. Pull-down efficiency was assessed with GST antibody and levels of bound ERK3 were verified. The nucleotide-loading status of the purified RAC1 and CDC42 proteins is presented in Figure 2-figure supplement 3. (J, K) Concentration-dependent binding affinity of ERK3-GST protein to (J) RAC1 and (K) CDC42 was determined. Interacting ERK3/GST proteins were used at 5, 10, 20, and 40 nM concentrations. Data from three independent experiments run in triplicates are presented as mean ± SEM. (L) ERK3 kinase domain (aa 9-327) binds to RAC1 and CDC42. Representative analysis of the in vitro GST pull-down of RAC1 and CDC42 and the interaction with the ERK3 kinase domain recombinant protein. Binding affinity of ERK3 (aa 9-327) with PAK1-PBD was verified and is presented in Figure 2-figure supplement 4.
The online version of this article includes the following source data and figure supplement(s) for figure 2: Source data 1. Full membrane scans for western blot images for Figure 2B, C, H, I, and L.
Source data 2. Prism and Excel file for Figure 2D.
Source data 3. Prism and Excel file for Figure 2E.
Source data 4. Prism and Excel file for Figure 2F.
Source data 5. Prism and Excel file for Figure 2G.
Source data 6. Prism and Excel file for Figure 2J.
Source data 7. Prism and Excel file for Figure 2K.              total CDC42 and RAC1 than control cells. Interestingly, subsequent pull-down of active CDC42/RAC1 from the PM fraction revealed that although both RhoGTPases are present at the plasma membrane in the absence of ERK3, their activity is significantly reduced ( Figure 4B and C). These results were corroborated by the colocalization of ERK3 with CDC42 to the protrusions at the leading edge of the cell ( Figure 4D-G) with mean Pearson's correlation coefficient (PCC (r)) of 0.6638 ± 0.02946 and Spearman's rank correlation coefficient (SRCC (ρ)): 0.7122 ± 0.02586 ( Figure 4F and G, respectively).
As ERK3 protein co-precipitated with active RAC1 and CDC42 ( Figure 2B and C), we further investigated whether also components of the ARP2/3 complex precipitate with RAC1/CDC42. Interestingly, we readily detected the ARP2/3 complex subunits ARP3, ARP2, and ARPC1A, as well as ERK3 by immunoblots in active RAC1/CDC42 pull-downs ( Figure 5B). Consistent with the reduction in active RAC1 and CDC42 levels, we observed that knockdown of ERK3 by CRISPR/Cas9 or with si/ shRNAs led to a reduction in F-actin staining in mammary epithelial cells as shown in Figure 5C and Figure 5-figure supplement 1A. To further quantify these effects and corroborate the role of ERK3 in polymerization of actin filaments, we evaluated the ratio of cytoskeleton-incorporated filamentous (F-) and globular/monomeric actin (G-actin) in control and ERK3-depleted mammary epithelial cells by ultracentrifugation. Knockdown of ERK3 significantly shifted total F-actin to G-actin in both primary HMECs ( Figure 5D and E) and MDA-MB231 cancer cells ( Figure 5-figure supplement 1B and C). Moreover, in the absence of CDC42, we could still detect colocalization between endogenous ERK3 and the ARP3 subunit of the ARP2/3 complex ( Figure 5-figure supplement 2A-D). These data prompted us to test whether ERK3 directly controls ARP2/3-dependent actin polymerization in the absence of RAC1 and CDC42.
In vitro, the ARP2/3 complex expresses low intrinsic actin nucleating activity. In vivo, its activity is tightly controlled by NPFs, which stimulate conformational changes (Mullins et al., 1998;Welch et al., 1998). To determine whether ERK3 had any direct effect on the assembly of actin filaments, we used a pyrene-actin polymerization assay to monitor fluorescence kinetics over time as the actin filaments assembled in the presence of the purified ARP2/3 complex. We used the VCA domain blot analyses of active RAC1 and CDC42 pull-down, respectively, using PAK1-PBD fusion beads. Levels of active RhoGTPases were detected using RAC1-and CDC42-specific antibodies. Levels of PAK1-PBD protein were detected by Ponceau S staining and used for the quantification presented in (C) and (E). Fold change in GTP-loading was calculated by normalization of the signal obtained for the samples with GTPγS in the presence of ERK3 with samples incubated with GTPγS alone and is presented as mean ± SEM from (D) five (n = 5) independent experiments; *p<0.0332, **p<0.0021, ***p<0.0002, ****p<0.0001, unpaired t-test. (F, G) In vitro RhoGEF activity assay was performed to assess guanine nucleotide exchange activity on (F) CDC42 and (G) RAC1 in the presence and absence of recombinant ERK3 protein (9-327 aa). After six initial readings, Dbs-GEF or ERK3 protein were added at 0.5 µM final concentration. GEF activity was expressed as mean relative fluorescence units (RFUs) from at least three independent experiments. RhoGEF Dbs protein was used as a positive control. In vitro GEF activity of full-length ERK3 vs. ERK3 kinase domain towards CDC42 was compared as is presented in Figure 3-figure supplement 1A.
The online version of this article includes the following source data and figure supplement(s) for figure 3: Source data 1. Full membrane scans for western blot images for Figure 3B and D.
Source data 2. Prism and Excel file for Figure 3C.
Source data 3. Prism and Excel file for Figure 3E.
Source data 4. Prism and Excel file for Figure 3F.
Source data 5. Prism and Excel file for Figure 3G.     of WASP protein as a positive control. Indeed, recombinant ERK3 stimulated ARP2/3-dependent actin polymerization at a nanomolar concentration (5 nM) in the presence of 10 nM of the ARP2/3 complex ( Figure 5F). Moreover, the nucleation-promoting activity of ERK3 was comparable to WASP (VCA), a well-known NPF. ERK3 protein alone did not exert any stimulatory effect on actin nucleation ( Figure 5F). These data suggest that ERK3 could function as nucleation promoting factor that stimulates ARP2/3-dependent actin polymerization. Additionally, we measured the effect of having both ERK3 and WASP (VCA) present on ARP2/3-dependent actin polymerization. Using the same concentration as for the initial screening (ERK3 4.8 nM and WASP [VCA] at 400 nM), we did not detect any additive effect when both proteins were combined (

ERK3 directly binds to ARP3
The ARP2/3 complex is composed of seven subunits: the actin-related proteins ARP2 and ARP3 and five associated proteins (ARPC1-5) that sequester ARP2 and ARP3 in the complex's inactive conformation (Robinson et al., 2001;Nolen et al., 2004; Figure 6A). In vitro, NPFs such as WASP stimulate actin polymerization by directly binding to and activating ARP2/3. Considering that ERK3 exerted a similar effect on the ARP2/3-dependent actin polymerization as NPFs, we further investigated the mode of interaction between ERK3 and the ARP2/3 protein complex by employing purified components. Indeed, using quantitative ELISAs we found that full-length ERK3 bound to the ARP2/3 complex with high-affinity in vitro ( Figure 6B).
ARP2 and ARP3 are essential for the nucleating activity of the ARP2/3 complex and occupy the central position in actin nucleation. Structural modeling has suggested that although ARP2 and ARP3 interact with each other in the inactive conformation, they are not in the right configuration to act as a seed for actin nucleation. After activation by NPF binding, the heterodimer mimics two actin subunits. Binding of an actin monomer then initiates the formation of the actin nucleus and the polymerization process (Mullins et al., 1998;Gournier et al., 2001;Kelleher et al., 1995). Extensive characterization of the individual subunits also revealed that ARP3 is crucial for the stimulation and activity of the ARP2/3 complex and nucleation in vitro (Gournier et al., 2001). Considering that ARP3 is a structural component regulating the nucleating properties of the ARP2/3 complex, we further determined whether ERK3 can directly bind to ARP3. Through GST pull-down experiments with purified recombinant proteins, we could indeed detect a concentration-dependent direct binding between full-length ERK3 and ARP3 ( Figure 6C). These results were further confirmed with ELISA analyses suggesting a high-affinity binding between these two proteins ( Figure 6D). The ERK3 kinase domain binds also directly to the Arp3 protein, albeit with low affinity in comparison to the full-length protein ( Figure 6-figure supplement 1A). Finally, we also detected ERK3 co-precipitating with ARP3 and ARP2 in HMECs at endogenous levels ( Figure 6E). instructions (Cat# 16118/19, Thermo Fisher) and as described in the 'Materials and methods' section. (B) Cellular fractionation was performed. Expression levels of ERK3, RAC1, and CDC42 were assessed in nuclear (N), cytosolic (C), and PM fractions. Na/K ATPases and GAPDH were used as controls for the analyzed fractions. (C) Active RAC1/CDC42 pull-down was performed from the isolated PM fraction. (D-G) Colocalization of ERK3 and CDC42 in polarized cells. Human mammary epithelial cells (HMECs) were seeded and cultured on cover slips. When cells became around 70% confluent, scratch wounds were introduced to the cover slip using a 200 µl tip, medium was exchanged to supplement-free, and cells were cultured for additional 6 hr. Afterward, cells were fixed and subjected to IF staining as described in the 'Materials and methods' section with anti-CDC42 (red) (secondary antibody: anti-mouse Cy3; Cat# A10521, Thermo Fisher Scientific) and anti-ERK3 (green) (secondary antibody: Alexa Fluor 488; Cat# A11008, Thermo Fisher Scientific) antibodies. Merged images show the colocalization of both proteins in yellow. Magnification of the boxed regions is shown on the right for better visualization. In (D), a group of cells at the scratch site is presented and ERK3-CDC42 colocalization is marked with red arrows at the cell protrusions and with white arrows at the cell body. (E) Images representing ERK3-CDC42 colocalization at a single-cell level at the scratch site. Scale bars 28 µm. (F, G) Colocalization of ERK3 and CDC42 was analyzed as described in the 'Materials and methods' section and values for (F) Pearson's correlation coefficient as well as the (G) Spearman's rank correlation coefficient are presented for eight randomly selected cells (n = 8). Scores above 0 indicate a tendency towards colocalization with a perfect colocalization with a score of 1.
The online version of this article includes the following source data for figure 4: Source data 1. Full membrane scans for western blot images for Figure 4A-C.
Source data 2. Prism and Excel file for Figure 4F and G.

ERK3 phosphorylates ARP3 at S418
Conformational changes activating ARP2/3 are induced by binding to NPFs, actin, and ATP. Additionally, phosphorylation of various ARP2/3 residues has been reported previously to play a crucial role in the conformational rearrangements within the ARP2/3 complex (Vadlamudi et al., 2004;Singh et al., 2003;Narayanan et al., 2011;LeClaire et al., 2008;Choi et al., 2013). Most of the studies focused on phosphorylation of the ARP2 subunit, which induced the conformational repositioning of this subunit toward ARP3, further allowing binding of NPFs and full activation of the ARP2/3 complex leading to the formation of the ARP2-ARP3 heterodimer (Narayanan et al., 2011;LeClaire et al., 2008;Choi et al., 2013). From the seven components of the ARP2/3 complex ARP3, ARP2 and ARPC1A purified from Acanthamoeba castellani were shown to be phosphorylated (LeClaire et al., 2008). Since ERK3 is a kinase, we assessed phosphorylation of the ARP2/3 complex by ERK3 by employing in vitro kinase assays and subsequent analysis of phosphopeptides by mass spectrometry. Initial analyses detected phosphorylation of the ARP3 subunit at serine 418 (S418) (peptide sequence: HNPVFGVMS, please see Supplementary file 1). We further validated the significance of the detected phosphorylation in vivo by overexpressing the non-phosphorylatable (S418A) and phospho-mimicking (S418D) mutants of ARP3 in primary mammary epithelial cells and analyzing the cell morphology and F-actin abundance. Cells transduced with wild type (WT) ARP3 or empty vector (EV) were used as reference controls ( Figure 6F). We were able to validate the phosphorylation of ARP3 at Ser418 in ARP3overexpressing cells ( Figure 6G) and detect a decrease in the S418 phosphorylation in ERK3-depleted the V-verpolin-like motif binds actin monomer (G-actin), C-central and A-acidic domains bind and activate the ARP2/3 complex. Conformational changes induced by the binding of the ARP2/3 complex promote its binding to the actin filament, which is strengthened by the additional interaction of the ARP2/3 complex with WASP (VCA)-G-actin. Further conformational changes will secure the ARP2/3 complex on the filament and allow its binding to the actin monomer and the polymerization of the newly nucleated filament. Actin polymerizes at the fast-growing/barbed end, elongating toward the plasma membrane and the ARP2/3 complex would cross-link newly polymerizing filament to the existing filament. (B) ERK3 co-precipitates with active RAC1 and CDC42 in complex with ARP2/3. Active RAC1/CDC42 pull-down was performed using control and ERK3 knockdown human mammary epithelial cells (HMECs). Levels of the active RAC1 and CDC42 were assessed as well as the co-immunoprecipitation levels of ERK3, ARP2, ARP3, and ARPC1A. Levels of the total protein expression were evaluated in the total cell lysates (TCL) and Ponceau S staining was used as a loading control. (C-F) ERK3 regulates F-actin levels in vitro and in vivo. (C) Western blot analyses of control (CRISPR Co) and ERK3-depleted (CRISPR ERK3) HMECs are presented alongside with representative confocal images of F-actin staining. (D, E) In vivo analysis of F-and G-actin levels in HMECs upon ERK3 knockdown. (D) Representative western blot analyses of the enriched F-and G-actin fractions as well as the ERK3 knockdown validation and total actin levels in the TCL are presented. (E) F-and G-actin levels were quantified, and ratios were calculated from five (n = 5) independent experiments and are presented as mean ± SEM; *p<0.0332, **p<0.0021, ***p<0.0002, ****p<0.0001, unpaired t-test. Analyses of ERK3-dependent regulation of F-actin levels in cancerous MDA-MB231 cells is presented in Figure 5-figure supplement 1. Cellular colocalization between endogenous ERK3 and the ARP2/3 was assessed in the absence of CDC42 and is presented in Figure 5-figure supplement 2. (F) Effect of full-length ERK3 on ARP2/3-dependent pyrene actin polymerization was assessed using a pyrene actin polymerization assay. Polymerization induced by the VCA domain of WASP that served as a positive control (green) as well as the ARP2/3 (orange) and ERK3 protein alone (blue) are shown for reference. Actin alone (black) was used to establish a baseline of polymerization. Fluorescence at 360/415 was measured over time and is presented as mean fold change from at least three independent experiments after normalization to the first time point within the respective group. ARP2/3-dependent actin polymerization was measured in the presence of both ERK3 and WASP (VCA) domain, and the results are depicted in Figure 5-figure supplement 3.
The online version of this article includes the following source data and figure supplement(s) for figure 5: Source data 1. Full membrane scans for western blot images for Figure 5B and C and Figure 4D.
Source data 2. Prism and Excel file for Figure 5E.
Source data 3. Prism and Excel file for Figure 5F.           Figure 6H) and in the kinase dead (KD) ERK3-overexpressing cells (Figure 6-figure supplement 1B). Of note, we observed that expression of exogenous ARP3 reduced the protein levels of the endogenous ARP3 ( Figure 6G), possibly due to disruption of the protein complex, thereby affecting the stability of the endogenous proteins. Strikingly, expression of the phosphorylation-mimicking ARP3 mutant (S418D) under these settings predominantly induced actin filament formation, as indicated by the intense phalloidin staining ( Figure 6F) and enhanced the F/G-actin ratio in primary mammary epithelial cells ( Figure 6I and J). Moreover, S418D-overexpressing cells exhibited F-actin-rich protrusions ( Figure 6F, lower panel). In contrast, most of the S418A-expressing cells were smaller in size and had a round morphology ( Figure 6F, upper panel). Quantification of the F-actin levels revealed that despite the morphological distortion, cells expressing non-phosphorylatable ARP3 (ARP3 S418A) had no significant effect on overall F-actin levels ( Figure 6I and J). To further determine the relevance of the S418 phosphorylation in ERK3-regulated ARP2/3-dependent actin polymerization in cells, we depleted ERK3 in the ARP3 S418D-overexpressing HMECs using shRNA ( Figure 6K-M). Knockdown by ELISA as described in the 'Materials and methods' section and mean absorbance (Abs) ± SEM from three independent experiments is presented. (E) Co-immunoprecipitation (IP) of ARP2/3 protein complex and ERK3 was performed in HMECs using ARP3 antibody. Levels of precipitated ARP3 as well as co-IP of ARP2 and ERK3 were assessed. IgG control was included to determine the specificity of the interaction. Total cell lysate (TCL) was included to present expression levels of the verified interacting partners. Ponceau S staining was used as a loading control. (F, G) Actin phenotype of the human mammary epithelial cells (HMECs) was validated upon stable overexpression of the ARP3 non-phosphorylatable (S418A) and the phosphomimicking (S418D) mutant, respectively. Wild type (WT) ARP3 was used as a control for the mutants and empty vector (EV) served negative control for the overexpression itself. (F) F-actin expression and organization in the negative (S418A) and phospho-mimicking (S418D) ARP3 mutant was visualized by green phalloidin and merged with the Hoechst staining of the nuclei. Four representative confocal images are presented. Images of EV-transfected and ARP3 WT-overexpressing HMECs are presented as controls. (G) Western blot validation of the overexpression efficiency and phosphorylation of ARP3 at S418. Anti-V5-tag antibody was used to detect levels of exogenous ARP3 WT, S418A, and S418D. Expression levels of the endogenous ARP3 were assessed as well as the phosphorylation at S418, total actin was validated. Ponceau S staining was used as a loading control. (H) Detection of the S418 phosphorylation of ARP3 in CRISPR ERK3 HMECs presented in Figure 5C and D. (I, J) Effect of the ARP3 mutant overexpression on F-actin levels was quantified using F/G actin in vivo assay. (I) Representative western blot analyses of F-and G-actin levels detected in fractions obtained from EV, ARP3 WT, S418A, S418D HMECs. (J) Quantification of the F/G actin ratios was performed for three (n = 3) independent experiments and is presented as mean ± SEM; *p<0.0332, **p<0.0021, ***p<0.0002, ****p<0.0001, one-way ANOVA, Tukey's post-test. (K-M) Effect of ERK3 depletion on dense F-actin phenotype of the ARP3 S418D-overexpressing HMECs. HMECs stably overexpressing ARP3 S418D were transduced with lentiviral particles targeting ERK3 (shERK3) and stable knockdown was established as described in the 'Materials and methods' section. Cells were further subjected to analyses of the F-actin levels. (K) IF staining with Oregon Green Phalloidin 488 to visualize F-actin levels and organization. Scale bars 28 µm. (L, M) Effect of the ERK3 knockdown on F-actin levels was quantified in the ARP3 S418D mutant overexpressing HMECs using F/G actin in vivo assay. (L) Representative western blot analyses of F/G actin levels. ARP3 S418D-(V5-tagged) overexpression and ERK3 knockdown efficiency were validated in TCL. Actin and Ponceau S staining were used as loading controls. (M) Calculated ratios of F/G actin are presented as mean ± SEM from three (n = 3) independent experiments; *p<0.0332, **p<0.0021, ***p<0.0002, ****p<0.0001, paired t-test. Colocalization of endogenous ERK3 with endogenous and exogenous ARP3 mutant (S418D) was verified, and further effect of the ERK3 depletion on the RAC1 and CDC42 activity was assessed in ARP3 S418Doverexpressing HMECs and presented in Figure 6-figure supplement 2.
The online version of this article includes the following source data and figure supplement(s) for figure 6: Source data 1. Full membrane scans for western blot images for Figure 6A, C, E, G, I, and L.
Source data 2. Prism and Excel file for Figure 6B.
Source data 3. Prism and Excel file for Figure 6D.
Source data 4. Prism and Excel file for Figure 6J.
Source data 5. Prism and Excel file for Figure 6M.         of ERK3 led to a significant reduction of the F-actin levels in the ARP3 S418D-overexpressing cells ( Figure 6K-M). These data suggested that although phosphorylation at S418 of ARP3 promoted actin polymerization and thus the F-actin content it is not absolutely essential like ERK3. A concomitant active RAC1/CDC42 pull-down revealed a significant decrease in RAC1 activity, with almost no effect on CDC42 (Figure 6-figure supplement 2A and B). We could still detect ARP3 and ARP2 in complex with CDC42 under these settings ( Figure 6-figure supplement 2A). Using immunofluorescence analysis, we observed significant colocalization between endogenous ERK3 and ARP3 in HMECs ( Figure 6-figure supplement 2C) and between endogenous ERK3 and overexpressed ARP3 S418D (Figure 6-figure supplement 2D-F), respectively. Further, depletion of ERK3, as expected, reduced the migration of ARP3 S418D-expressing cells (Figure 6-figure supplement 2G and H).
Kinase activity of ERK3 is necessary for membrane protrusions in primary mammary epithelial cells ERK3 controls both F-actin levels by regulation of the new filament assembly via ARP3-binding and its branching and bundling into the actin-rich protrusions by binding to CDC42/RAC1. We further investigated whether kinase activity of ERK3 was required for the formation of actin-rich protrusions and F-actin levels in mammary epithelial cells. Control (shCo) and ERK3-depleted (shERK3, 3′ UTR) HMECs were reconstituted with WT or kinase dead (KD) (K49A/K50A) ERK3, and EV was introduced as a control ( Figure 7A). Knockdown of ERK3 led to a decrease in actin-rich protrusions and overall F-actin staining ( Figure 7B) with concomitant decrease in F-actin abundance ( Figure 7C and D). Complementation with ERK3 WT rescued overall F-actin levels and the protrusive phenotype of the HMECs ( Figure 7B and D). Interestingly, KD ERK3 recovered cytoskeletal F-actin levels to the same extent as ERK3 WT ( Figure 7C and D). However, although the abundance of F-actin is rescued upon KD ERK3 overexpression, cells formed a dense meshwork of actin membrane ruffles which did not protrude into filopodia but rather curved around the edges of the cell forming a tangled web ( Figure 7B). These results indicate that kinase activity of ERK3 is not required for the polymerization of actin, but rather for the bundling and/or branching of the rapidly polymerizing actin filaments. To further corroborate these observations, the WT and KD ERK3 proteins were expressed in the rabbit reticulocyte lysate (RRL) system and then employed to stimulate ARP2/3-dependent actin polymerization in vitro. Kinase activity of ERK3 did not seem to affect the ARP2/3-dependent actin polymerization in vitro (Figure 7-figure supplement 1). We performed active RAC1/CDC42 pull-down from HMECs shERK3 (3′ UTR) cells reconstituted with either WT or the KD ERK3 and found that the kinase activity of ERK3 was not required for interactions or activation of CDC42. However, RAC1 activation was partially dependent on the kinase activity of ERK3 ( Figure 7E-I).

Discussion
The regulation of the actin cytoskeleton is an intensely investigated area because of its fundamental role in many basic cellular functions. The signaling machinery that controls the actin skeleton dynamics is deregulated in cancer contributing to metastasis. Kinases are major drivers of signaling events and form the largest part of the druggable genome. Several drugs have successfully been developed to target deregulated kinases in cancer. Of the more than 500 kinases encoded by the human genome, many studies have been focused on the conventional MAPKs, while the pathophysiological significance of atypical MAPKs remains underexplored. We have recently demonstrated that ERK3 directly contributed to AP-1/c-Jun activation, IL-8 production, and chemotaxis (Bogucka et al., 2020). Emerging studies have further shown that ERK3 functions in a context-and tissue-specific manner in controlling tumorigenesis and metastasis (Al-Mahdi et al., 2015;Bogucka et al., 2021;Elkhadragy et al., 2020;Elkhadragy et al., 2017;Choi et al., 2013;Elkhadragy et al., 2018;Alshammari et al., 2021). Mice expressing a catalytically inactive ERK3 mutant survive to adulthood and are fertile, but the kinase activity of ERK3 is necessary for optimal postnatal growth. In this study, we aimed to elucidate how ERK3 influences cell shape and actin cytoskeleton and unexpectedly revealed a pivotal role for this atypical MAPK in the regulation of RhoGTPases as GEF, of ARP2/3-dependent polymerization of actin probably as NPF, and of cell migration in mammalian cells (Figure 8). Further, our studies unveiled an evolutionarily conserved role for ERK3 in the formation of actin-rich protrusions.

ERK3 is a GEF for CDC42
The activation of both RAC1 and CDC42 was compromised in the absence of ERK3 in primary mammary epithelial cells and oncogenic MDA-MB-231 cells. While the activation of RAC1 was partially rescued in HMECs upon EGF stimulation in ERK3-depleted cells, this was not the case in MDA-MB231 cells ( Figure 2B, C, F, and G). Furthermore, while ERK3 regulated the activity of both RAC1 and CDC42 in vivo ( Figure 2B-G), the kinase domain of ERK3, which lacks kinase activity, directly facilitated GDP-GTP nucleotide exchange specifically for CDC42 ( Figure 3B, C, and G), which is implicated as a main regulator of filopodia assembly. In some cell types, we often detect a reduction in the total protein levels of RAC1 and CDC42 after depletion of ERK3. Whether inactivation of these RhoGTPases prompts them for proteasomal degradation under these settings deserves further analysis. Using several experimental procedures, we demonstrated that the interaction between ERK3 and RAC1 or CDC42 is direct and nucleotide-independent, which suggests that ERK3 is probably targeting the C-terminus of these RhoGTPases. Interestingly, it has been reported that RAC1 possesses a C-terminal docking (D site) for the canonical ERK ( 183 KKRKRKCLLL 192 ) (Tong et al., 2013). Whether the same site is being exploited by ERK3 for its binding to RAC1 or CDC42 deserves further studies. It is interesting to note that while the kinase domain of ERK3 promoted GTP binding to CDC42, it slightly attenuated the same process with RAC1 ( Figure 3B-G). While the kinase activity of ERK3 is not required for CDC42 activation, RAC1 activity partially depends on it. Further structural, biochemical, and biophysical studies are clearly warranted to unveil the molecular basis of this interaction and the subsequent functional implications.
In cells, regulation of RhoGTPases and that of the actin assembly is a much more complex and integrative process (Kurokawa et al., 2004). For example, a cell type-dependent, hierarchical crosstalk between RhoGTPases was observed in the regulation of actin-rich protrusions and CDC42 activity was shown to initiate RAC1-dependent lamellipodia protrusions (Nobes and Hall, 1995;Zamudio-Meza et al., 2009;Nobes and Hall, 1999). While ERK3 might activate RAC1 through its GEF-like function toward CDC42, it could also function as an intermediate link between RAC1 and its GEF by scaffolding both proteins and/or activating the GEF, thus indirectly contributing to RAC1 activity.
It is highly interesting that PAK family kinases, the effector kinases of RAC1 and CDC42, directly phosphorylate the SEG motif in the activation loop of ERK3 (Déléris et al., 2011;De la Mota-Peynado et al., 2011). The observations presented here strongly suggest a positive feedback loop, where ERK3 is required for the GTP loading of CDC42 and RAC1 to induce their interaction with PAKs, which in turn phosphorylate and activate ERK3. We show that active ERK3 kinase is required for the eventual formation of actin-regulated membrane protrusions. The functional and physical uncoupling of this loop could possibly be controlled by a phosphatase or by inducing the rapid turnover of ERK3 protein.
The constituents of this dynamic complex need further evaluation and characterization.
representing cell phenotype and magnified regions of the cell edges on the right show the actin distribution. (C, D) Levels of F-actin in HMECs analyzed by phalloidin staining were quantified using a F/G actin in vivo assay. (C) Calculated ratios of F/G actin are depicted as mean ± SEM from three (n = 3) independent experiments; *p<0.0332, **p<0.0021, ***p<0.0002, ****p<0.0001, one-way ANOVA, Tukey's post-test. (D) Representative western blot analyses of F-and G-actin levels. (E-I) Active RAC1 and CDC42 pull-down assays were performed from control (shCo) and ERK3-depleted (shERK3, 3′ UTR) HMEC reconstituted with either WT or KD ERK3 as described in the 'Materials and methods' section. PAK1-PBD was used to capture the active form of CDC42 and RAC1 in the respective cell lysates. Levels of active (GTP-bound) CDC42 and RAC1, as well as the total protein expression were assessed, ERK3 knockdown efficiency and overexpression was verified in the total cell lysate (TCL) using ERK3 antibody and V5-tag antibody to detect the exogenous WT and KD version of ERK3. (E) Western blot analyses are presented. (F, G) Relative levels of active CDC42 and RAC1 were calculated with respect to the total protein levels and are presented as ratio of GTP/Total RhoGTPases. (H, I) Additionally, the ratio of GTP/total RAC1 and CDC42 was normalized to the control cells (shERK3+EV) and is presented as fold change in activation.
The online version of this article includes the following source data and figure supplement(s) for figure 7: Source data 1. Full membrane scans for western blot images for Figure 7A, D, and E.
Source data 2. Prism and Excel file for Figure 7C.     . Schematic summary of the ERK3-dependent mechanisms regulating actin cytoskeleton and cell motility. ERK3 directly binds and activates the ARP2/3 protein complex as well as the CDC42 and RAC1 Rho GTPases. Activation of the ARP2/3 complex and RAC1/CDC42 is required for nucleation of the new actin filaments, elongation, and branching into the lamellipodia and filopodia. ERK3 regulates actin-rich protrusions, which play a direct role in cell motility.

ERK3 binds ARP3 to regulate its function in actin polymerization
Cells rely on several mechanisms to regulate the assembly of actin filaments. To complicate matters further, regulatory pathways involved in the filopodia formation seem to be cell type-specific, and CDC42-, WASP-, and ARP2/3-independent mechanisms of membrane spikes formation have been proposed (Czuchra et al., 2005;Snapper et al., 2001;Steffen et al., 2006). In this study, we demonstrate that ERK3 directly binds the ARP2/3 complex and contributes to actin polymerization (Figures 5  and 6). Conformational changes in ARP2/3 complex assembly are crucial for the initial heterodimer formation between ARP2 and ARP3 subunits (Machesky et al., 1999;Gournier et al., 2001;Welch et al., 1997). In the unstimulated complex, ARP2 and ARP3 subunits exist in a splayed confirmation and upon stimulation by NPFs the two proteins align into a side-by-side position creating a surface mimicking the first two subunits of the new actin filament (Robinson et al., 2001;Mullins and Pollard, 1999). Binding of these two crucial subunits by NPFs could reassure an active conformation. Interestingly, phosphorylation of certain subunits of ARP2/3 complex has been shown to be indispensable for the destabilization of its inactive conformation prior the full activation by the NPFs (Vadlamudi et al., 2004;Narayanan et al., 2011;LeClaire et al., 2008;Choi et al., 2013). LeClaire et al. further reported that Nck-interacting kinase (NIK) regulates ARP2 phosphorylation and grooms the ARP2/3 complex for the activation by NPFs. Phosphorylation events within the ARP2/3 complex could regulate nucleation activity by affecting several properties, such as complex formation itself, conformational changes in the complex assembly, or the affinity of the ARP2/3 to NPFs. Further structural simulation studies are needed to test whether the interaction of ERK3 with the ARP2/3 complex and ARP3 phosphorylation at S418 inflicts key conformational changes to support actin polymerization. While we uncover a role for ERK3 in the formation of filopodia and ARP2/3-dependent actin polymerization, the possible role for ERK3 in the regulation of other actin nucleation factors including formins cannot be ruled out.
While the kinase activity of ERK3 is clearly required for the formation of actin-rich protrusions, expression of KD constructs in ERK3-depleted cells failed to inhibit F-actin enrichment. These data suggest that perhaps ERK3 binding to ARP2/3 and possibly other hitherto unknown factors per se contribute to actin polymerization, while its kinase activity is probably still required for ARP2/3dependent and CDC42-induced filopodia formation, polarization, and migration of mammalian epithelial cells.
Overall, our studies shed new light on the multilayered regulation of the actin cytoskeleton by an understudied MAPK (Figure 8). ERK3 is druggable, and kinase inhibitors have already been developed for clinical use, which could also serve as tools to decipher kinase-dependent and -independent functions. Our observations will not only enhance the current understanding of actin cytoskeleton regulation, but also instigate new lines of investigations on the molecular, structural, and functional characterization of RhoGTPases and the ARP2/3 complex.

Materials and methods
Cell culture Medical Center of the Johannes Gutenberg University Mainz) and were cultured on 0.2% gelatin coating in Endothelial Cell Growth Medium (Cat# C-22010, PromoCell). HT-29 cells were purchased from ATCC (ATCC HTB-38) and cultured in McCoy's medium supplemented with 10% heat-inactivated FBS. Calu-1 cells were obtained from Sigma and cultured in DMEM supplemented with 10% heatinactivated FBS.
All cell lines used in this study were authenticated cell lines obtained from ATCC or DSMZ. All used cells were periodically tested for Mycoplasma contamination with negative results.

Stimulation of cells
HMECs and MDA-MB231 cells were seeded in 12-well plates at an initial density of 2 × 10 5 cells/ well. After cells reached 70% confluence, medium was exchanged to MEGM with no supplements for HMECs and DMEM minus FBS for MDA-MB231 cells, 4 hr prior treatment with recombinant human epidermal growth factor (EGF) (Cat# RP-10927, Invitrogen) at 5, 10, 15, and 30 min. Afterward, cells were subjected to western blot analyses.

Cycloheximide (CHX) chase experiments
To determine ERK3 half-life in mammary epithelial cells and its alteration upon EGF treatment, HMECs were seeded in 12-well plates at an initial density of 2 × 10 5 cells/well and cultured until 70% confluent. Medium was exchanged to MEGM with no supplements, 4 hr before cells were treated with human recombinant EGF at 100 µg/ml for 15 min, followed by treatment with protein synthesis inhibitor CHX (Cat# C-7698, Sigma, stock 100 mg/ml in DMSO) at 0.5, 1, 2, 4, and 6 hr. Cells were subjected to western blot analyses and quantification. Fold change in ERK3 protein levels was calculated with respect to the untreated cells (-EGF, 0 hr) or to the respective control in each group (0 hr) for unstimulated (-EGF) and EGF-stimulated (+EGF) cells using ImageJ software.  Figure 6G). ARP3 (Cat# ab151729, Abcam), ARP2 (Cat# ab128934, Abcam), and beta-actin HRP-conjugated antibody (Cat# ab49900, Abcam). HRP-conjugated secondary antibodies for rabbit and mouse IgG were obtained from Invitrogen (Cat# A16096 and A16066, respectively).

siRNA transfection
Cells were seeded 1 d before transfection at an initial density of 2 × 10 5 cells/well in 12-well plates or at 3 × 10 5 cells/well (6-well plates). Cells were transfected using SAINT-sRNA transfection reagent (SR-2003, Synvolux) according to the manufacturer's instructions. Cells were analyzed 48 hr posttransfection, and knockdown efficiency was verified by western blot using target-specific antibodies.
If EGF stimulation was included in the experimental settings, 48 hr post-transfection medium was exchanged for serum/supplements-free medium for 4 hr prior the EGF (100 ng/ml) treatment for 15 min.
The images of both GFP fluorescence and transmission of the live cells were acquired with a Leica SP8 confocal microscope (Leica, Mannheim, Germany) performed with a 10 × 0.3 NA objective, with 488 nm excitation (at approximately 150 µW), with an emission window of 500-590 nm for GFP detection and with scanning Differential Interference Contrast transmission imaging in a 1500 µm × 1500 µm frame format with 400 lines per second, 0.71 µm/pixel (2048 × 2048 pixel per frame) and with two times averaging per line with a frame acquisition of every 30 min per selected position within the chamber.
Phosphorylation of the ARP2/3 protein complex subunits by ERK3 was detected by mass spectrometry analyses after in vitro kinase assay. Proteins were in-gel-digested and subsequently enriched for phosphopeptides using TiO 2 beads.

RNA isolation, cDNA synthesis, and RT-PCR analysis
For gene expression analyses, cells were washed with cold PBS and total RNA was extracted using Trizol (Cat# 15596018, Ambion) according to the manufacturer's instructions. Isolated RNA was then used as a template for cDNA synthesis with the RevertAid First Strand cDNA synthesis kit (Cat# K1621, Thermo Fisher Scientific) and random hexamer primers.
Proteins size and purity were assessed throughout the purification procedures.
Synthesis of WT and KD ERK3 proteins using RRL ERK3 WT or ERK3 KD (K49A K50A) in pcDNA3/V5-Dest40 vector were expressed in T7 RRLs using TNT T7 Quick Coupled transcription/Translation System (Cat# TM045, Promega) following the manufacturer's protocol. Expression levels were assessed by western blot, using V5-tag specific antibody.
Lysates were further used in actin polymerization assay as recombinant proteins.
Endogenous pull-down of active RAC1/CDC42 Levels of GTP-bound RAC1 and CDC42 were determined using either active RAC1/CDC42 pull-down and detection kit (Cat# 16118/19, Thermo Fisher Scientific) according to he manufacturer's instructions and buffer composition or purified GST-PAK1-PBD fusion beads. HMEC and MDA-MB231 cells were seeded in 6-well plates at an initial density of 3 × 10 5 cells/well. After cells reached 70% confluence, medium was exchanged to MEGM (no supplements) for HMEC or DMEM without FBS for MDA-MB231 cells, 4 hr prior 15 min stimulation with 100 ng/ml of recombinant human EGF. Afterward, cells were subjected to active RAC1/CDC42 pull-down. The precipitated samples and total cell lysates were subjected to western blot analyses and ImageJ quantification. Relative levels of active RAC1 and CDC42 were determined by calculating the ratio of active (GTP loaded) RhoGTPases with respective total protein levels in TCL.

Endogenous pull-down of ARP3
HMECs were seeded in 6-well plate at an initial density of 2 × 10 5 cells per well and cultured until 80% confluent. For the immunoprecipitation (IP) of the endogenous ARP3, cells were washed with ice-cold PBS and lysed with ice-cold IP buffer (10 mM HEPES pH 7.4; 150 mM NaCl, 1% Triton X-100, plus protease inhibitor cocktail Set I-Calbiochem 1:100 [Cat# 539131, Merck Millipore], 1 mM Na 3 VO 4 , and 1mM NaF). Cell lysates were incubated with Protein A/G-Agarose beads (Cat# 11 134 515 001/11 243 233 001, Roche), and either ARP3 or normal rabbit IgG antibody for 2 hr at 4°C with rotating. After the incubation, beads were washed with IP buffer and analyzed by immunoblot. Levels of the immunoprecipitated ARP3 as well as the co-immunoprecipitation of ARP2 and ERK3 were assessed. Total cell lysates were used as a control.

Binding interaction-GST pull-down assays
Purified recombinant GST-fusion RAC1 and CDC42 proteins immobilized on the beads were used for in vitro GST pull-down experiments.
To verify the relevance of the nucleotide binding in the interaction of the GTPases with ERK3 protein, RAC1 and CDC42 proteins were loaded with non-hydrolyzable GTPγS or GDP (components of the active RAC1/CDC42 pull-down and detection kit; Cat# 16118/19, Thermo Fisher Scientific).
Recombinant GST-fusion GTPases protein beads were incubated with gentle shaking for 15 min at 30°C in 100 µl of 25 mM Tris-HCl, pH 7.2, 150 mM NaCl, 5 mM MgCl 2 , 1% NP-40, and 5% glycerol binding/wash buffer containing 0.1 mM of GTPγS or 1 mM GDP in the presence of 10 mM EDTA pH 8.0 to facilitate nucleotide exchange. Reaction was terminated by placing samples on ice and addition of 60 mM of MgCl 2 . GTPγS-/GDP-loaded RAC1 and CDC42 beads were centrifuged and supernatants were removed. Beads were subjected to GST pull-down ERK3 binding assay.
GTPγS/GDP-loaded GST-RAC1 and GST-CDC42 protein beads were incubated with recombinant human ERK3 protein (Cat# OPCA01714, Aviva Systems Biology) in 100 µl of binding/wash buffer supplemented with protease inhibitor with protease inhibitor cocktail Set I-Calbiochem 1:100 (Cat# 539131, Merck Millipore) and phosphatase inhibitors: 1 mM sodium orthovanadate (Na 3 VO 4 ), 1 mM sodium fluoride (NaF) for 1 hr at 4°C with rotation. Recombinant GST protein bound to glutathione beads was used as a negative control. After the incubation, beads were washed three times with 400 µl binding/wash buffer, centrifuged each time at 2000 rpm. Samples were eluted with 4× Laemmli buffer supplemented with 100 mM DTT and 5 min boiling at 95°. Samples were subjected to western blot analyses.
For the in vitro interaction between GST-RAC1/ GST-CDC42 and ERK3 kinase domain, RAC1/ CDC42 WT beads and ERK3 (amino acids (aa) 9-327) (Crelux) Escherichia coli purification was employed. Additional information on the Crelux protein can be viewed in the Supplementary file 1.
For ERK3-ARP3 binding studies, human recombinant GST-fusion ERK3 protein (SignalChem) and E. coli purified ARP3. Glutathione beads were incubated with 37 nM of GST-ERK3 or GST alone and indicated concentrations of ARP3 for 2 hr at 4°C in the binding buffer (5 mM Tris-HCl, pH 8.0, 0.2 mM CaCl 2 ) supplemented with protease and phosphatases inhibitors. Afterward, beads were washed three times with 400 µl binding buffer supplemented with 1% NP-40. Samples were eluted with 4× Laemmli buffer supplemented with 50 mM DTT and 5 min boiling at 95°. Immobilized protein complexes were detected by GST and ARP3-specific antibodies.

Determination of GTP loading on Rho GTPases
In vitro active RAC1/CDC42 pull-down assay The efficiency of the nucleotide loading under chosen conditions of the GTP-loading status of the purified RAC1 and CDC42 proteins was determined using GST-PAK1-PBD fusion beads.
Recombinant RAC1 and CDC42 proteins (no tag) were loaded with GTPγS or GDP by incubation for 15 min at 30°C in 100 µl of 25 mM Tris-HCl, pH 7.2, 150 mM NaCl, 5 mM MgCl 2 , 1% NP-40, and 5% glycerol binding/wash buffer containing 0.1 mM of GTPγS or 1 mM GDP in the presence of 10 mM EDTA pH 8.0. Reaction was terminated by placing samples on ice and addition of 60 mM of MgCl 2 . GTPγS-/GDP-loaded RAC1 and CDC42 as well as the native (WT) proteins were subjected to active RAC1/CDC42 pull-down assay using GST-PAK1-PBD fusion beads. Levels of GTP-loaded GTPases were detected using RAC1 or CDC42-specific antibodies.
GDP/GTP nucleotide exchange assays GDP/GTP exchange assay was performed using recombinant RAC1 and CDC42 proteins in binding buffer containing: 20 mM Tris-HCl (pH 7.5), 100 mM NaCl, and protease/phosphatases inhibitors. Firstly, RhoGTPases were stripped of nucleotides by incubation with binding buffer containing 10 mM EDTA for 10 min at RT. Afterward, reactions were supplemented with 50 mM MgCl 2 and 500 µM GDP (Active RAC1/CDC42 pull-down kit, Thermo Fisher, Cat# 16118/19) and incubated for 15 min at 37°C with gentle agitation. Samples of GDP-loaded RhoGTPases were kept on ice for further analyses as controls. The remaining GDP-loaded RAC1 and CDC42 were then incubated with 500 µM GTPγS (Thermo Fisher kit Cat# 16118/19) in the presence or absence of 2 nM of recombinant ERK3 protein (aa 9-327) (Crelux) for 30 min at 37°C with gentle agitation. All reactions were terminated by addition of 60 mM MgCl 2 on ice.
To isolate the active (GTP-bound) RAC1 and CDC42, the specific GST-PAK1-PBD-fusion beads were used. Nucleotide exchange reaction were incubated with the beads for 1 hr at 4°C, rotating. Afterward, beads were washed three times with binding/wash buffer: 25 mM Tris-HCl, pH 7.2, 150 mM NaCl, 5 mM MgCl 2 , 1% NP-40, and 5% glycerol. Samples were eluted with 4× sample buffer supplemented with 100 mM DTT. Levels of active RAC1 and CDC42 were detected using RAC1 and CDC42specific antibodies from the active RAC1/CDC42 pull-down kit (Cat# 16118/19, Thermo Fisher). Levels of GST-PAK1-PBD protein were detected using Ponceau S staining of the membrane. Efficiency of the GDP-GTP exchange was calculated as the ratio of GTPγS+ERK3 samples and GTPγS alone and presented as fold change where GTPγS samples were used as a reference.
Guanine nucleotide exchange (GEF) activity of ERK3 (9-327 aa) was further evaluated using RhoGEF Exchange Assay (Cat# BK100, Cytoskeleton). Assay was performed according to the manufacturer's instructions and included Dbs GEF as a positive control.

In vivo F-actin/G-actin assay
HMECs and MDA-MB231 cells were seeded in 12-well or 6-well plates for F/G actin assay and parallel western blot analyses. When cells reached about 70% confluency, plates were subjected to either western blot analyses or G-actin/F-actin in vivo assay using cytoskeleton kit (Cat# BK037) according to the manufacturer's instructions. Briefly, cells were washed with 1× PBS and lysed in appropriate volume of lysis and F-actin stabilization buffer supplemented with 1 mM ATP and protease inhibitor mixture for F-actin stabilization (Cytoskeleton). Lysates were centrifuged at 100,000 × g at 37°C for 1 hr to separate the G-actin fraction in the supernatants and pelleted F-actin. F-actin pellets were further depolymerized for 1 hr on ice using F-actin depolymerizing buffer. Both sets of samples were mixed with 4× SDS-PAGE sample buffer. Equal volumes of G-actin and F-actin lysates were analyzed by SDS-PAGE and immunoblotting with anti-actin rabbit polyclonal antibody (Cat# AANO1, Cytoskeleton). Densitometric analyses of the G-actin and F-actin levels were performed using ImageJ software. Ratios of F-actin in respect to G-actin levels were calculated.

In vitro actin polymerization assay
Actin polymerization assay was performed using actin polymerization biochem kit (Cat# BK003, Cytoskeleton) according to the provided protocol. Pyrene-labeled rabbit skeletal muscle actin (Cat# AP05-A, Cytoskeleton) (2.3 µM per reaction) was diluted with general actin buffer (5mM Tris-HCl pH 8.0, 0.2 mM CaCl 2 ) supplemented with 0.2 mM ATP and 1 mM DTT and incubated for 1 hr on ice to depolymerize actin oligomers. Actin stock was centrifuged for 30 min at 14,000 rpm at 4°C. Pyrenelabeled actin was incubated with recombinant ARP2/3 protein complex alone (Cat# RP01P, Cytoskeleton) or along with the human recombinant proteins: WASP-VCA domain protein (Cat# VCG03, Cytoskeleton) or ERK3 protein (M31-34G, SignalChem). ARP2/3 protein complex, WASP-VCA domain protein, and full-length ERK3 recombinant protein were used at 10 nM, 400 nM, and 4.8 nM final concentration, respectively. Actin alone was used to establish a baseline of polymerization rate. Actin polymerization was induced by addition of 1.5× actin polymerization buffer (Cat# BSA02) (10× buffer: 20 mM MgCl 2 , 500 mM KCl, 10 mM ATP, 50 mM guanidine carbonate in 100 mM Tris-HCl, pH 7.5). Actin polymerization was measured by fluorescence emission at 415 nm over 1 hr (30-60 s interval time) at RT with 360 nm excitation wavelength using multimode microplate reader Tecan SPARK. Buffer background signal at each interval was subtracted and relative fluorescence units (RFUs) are depicted.

Immunofluorescence (IF) and confocal analyzes
Cells cultured on coverslips were fixed in 3.7% formaldehyde (Cat# CP10.1, Roth) for 15 min, followed by washing with PBS pH 7.5 and 3 min permeabilization using 0.1% Triton X-100 (AppliChem). After washing twice with PBS, cells were blocked with 1% BSA in PBS for 15 min and washed once with PBS. Filamentous actin was labeled with Oregon Green 488 Phalloidin (Cat# O7466, Invitrogen) or alternatively Rhodamine Phalloidin (Cat# R415, Invitrogen) in blocking buffer for 1 hr at RT in the dark. Nuclei were stained with 10 µg/ml of DNA dye (Hoechst 33342; Cat# H3570, Invitrogen). For co-staining of endogenous ERK3 and F-actin, cells were incubated in 1:250 dilution of anti-ERK3 antibody (mouse Cat# MAB3196, R&D or rabbit Cat# ab53277, Abcam) in blocking solution for 1 hr at RT. For the ARP3 and CDC42 staining, the anti-ARP3 antibody (Cat# ab151729, Abcam) and anti-CDC42 antibody (Cat# 610929, BD Transduction) were used in blocking buffer for 1 hr at RT in the dark. Nuclei were stained with 10 µg/ml of DNA dye (Hoechst 33342; Cat# H3570, Invitrogen). Afterward, cells were washed with PBS and incubated with secondary anti-mouse IgG-Cyanine3 (Cat# A10521, Thermo Fisher Scientific) at 5 µg/ml, DNA dye (Hoechst 33342), and Green Phalloidin in blocking solution for 1 hr at RT in the dark. Samples were washed twice with PBS and cells were mounted onto glass slides using Moviol (+DABCO) (Sigma). For exogenous ERK3 detection upon overexpression of the ERK3 WT or serine 189 phosphorylation mutants (S189A/S189D) in shERK3 HMECs, V5-tagspecific antibody was used at 1:100 dilution (Cat # R960-25, Invitrogen). Cells were imaged using a Leica DMi8 confocal microscope (×63, oil immersion objective). ERK3-depleted cells were used as a control for the ERK3 staining.
For the colocalization analyses, a mask was created to select excluding nuclei where no colocalization was observed. The ARP3 and ERK3 images were scaled and thresholded in the same way to maximize the dynamic scale and then had the colocalization of the Pearson's correlation and Spearman's rank correlation coefficients (PCC and SCC) on a pixel-wise basis with the 'Coloc 2' algorithm in ImageJ. For both the Pearson and Spearman correlation coefficients, a value of +1 is perfectly correlated or colocalized, a value of 0 means that the two signals are randomly localizing, and a value of -1 indicates that the two signals are perfectly separated.

Scanning electron microscopy (SEM)
Control and ERK3-depleted HMECs were seeded in 12-well plate. Cells were fixed overnight with 2.5% glutaraldehyde, rinsed with PBS, and then again fixed and stained with 2% osmium tetroxide. Afterward, samples were rinsed with distilled water, frozen, and then freeze-dried with Crist Alpha LSC Plus Freeze Drier. The SEM scans were performed with a Philipps ESEM XL30 scanning electron microscope.

Filopodia quantification and analyses
To detect, quantify, and analyze the number and lengths of the filopodia of the fixed and actinphalloidin-stained control and ERK3 knockdown HMECs, the FiloQuant open access software and routines (Jacquemet et al., 2017;Jacquemet et al., 2019) were applied in ImageJ and in some cases were applied in Imaris version 9.3.1 (Imaris, RRID:SCR_007370) for additional verification. Briefly, the raw cell images were linearly adjusted for brightness and contrast for optimal filopodial observation and also had a mask applied to just analyze the filopodial region of interest in each cell and applied to each cell in the same way. Additionally, the following parameters were initially applied to each cell in the same way: Cell Edge Threshold = 20, Number of Iterations = 10, Number of Erode Cycles = 0, Fill Holes on Edges = checked, Filopodia Threshold = 25, Filopodia Minimum Size = 10 pixel. The results of the detected filopodial filaments were the overlayed in white over the actin-phalloidin green fluorescence of the filapodia for visual verification in ImageJ and in some cases in Imaris version 9.3.1.
The images of both GFP fluorescence and transmission of the live cells were acquired with a Leica SP8 confocal microscope using a 10 × 0.3 NA objective, with 488 nm excitation (at approximately 150 µW) and emission window of 500-590 nm for GFP detection and with scanning differential interference contrast transmission imaging in a 1500 µm × 1500 µm frame format with 400 lines per second, 0.71 µm/pixel (2048 × 2048 pixel per frame), and with two times averaging per line with a frame acquisition of every 30 min per selected position within the chamber. The images were first acquired in multiple regions of each cell type for 20 h.
The image sequences were imported into Imaris version 9.3.1 and detected automatically by fluorescence with both the whole-cell spot and whole-cell surface analysis. As the surface and spot analysis center positions differed in control analysis by less than 1%, the whole cell spot automated analysis was applied to all images in the same way with a 16 µm per cell diameter estimate. The automated tracking also occurred in the same way for all image sequences within Imaris with the autoregressive motion algorithm, with a maximum average distance of 60 µm per step and with zero step gap applied. The speed (or cell speed) was determined by dividing each x, y step by 1800s (or 30 min). The acceleration (or cell acceleration) was determined by subtracting a step speed from the previous step speed and dividing by 1800s (or 30 min). The displacement length was determined by subtracting the initial x, y position from the final x, y position and then to determine the difference vector length for each track over 20 hr of acquisition. The track length added each absolute x, y vector step for an entire track over 20 hr of acquisition. The track mean speed averaged the speed of the steps for individual tracks.
The graphs were created with GraphPad Prism (RRID:SCR_002798), and cell distribution according to the analyzed parameter was visualized by violin plots. Significance was determined using nonparametric Mann-Whitney test.

Transwell cell migration assay
Migratory properties of control and ERK3-depleted HMECs and MDA-MB231 cells were assessed using two-chamber Transwell system (Cat# 3422, Corning). Cells were seeded in 12-well plates. HMECs were deprived of supplements from the medium 24 hr prior the migration. HMECs and MDA-MB231 cells were then trypsinized and resuspended in supplements-free or serum-free medium, respectively. 500 µl of supplements-free/serum-free medium was mixed with 100 ng/ ml of EGF and added into the bottom chamber and 1 × 10 5 cells in 120 µl of supplements-free/ serum-free medium were added into each insert. Plates were incubated at 37°C for 24 hr. To quantify the migrated cells, cells were removed from the upper surface of the insert using cotton swabs and inserts were washed with PBS. Afterward, cells that migrated to the lower side were fixed within 3.7% formaldehyde for 10 min and stained with Hoechst solution in PBS for 20 min at 37°C. Cells were then visualized using Leica DMi8 microscope (×5 dry objective), images of 3-4 regions per each experimental condition were taken. Number of cells was quantified by Fiji/ImageJ software (Fiji, RRID:SCR_002285) using particle analyses for each field of view and averaged for each membrane. Directional migration was presented and the percentage of the respective control (expressed as 100%).

Phosphorylation site identification on ARP3
Phosphorylation of the ARP2/3 protein complex subunits by full-length ERK3 was detected by in vitro kinase assay and in-gel digestion of the respective ARP2/3 complex proteins followed by subsequent spectrometry analyses.

Phosphopeptide enrichment
Dried tryptic peptide samples were dissolved in loading buffer (1 M glycolic acid, 6% trifluoroacetic acid, 5% glycerol, and 80% acetonitrile) under continuous shaking. TiO 2 beads (Titansphere, TiO 2 , GL Sciences Inc) were washed in loading buffer three times before transferring them to the dissolved tryptic peptide samples. After 1 hr of continuous shaking, the supernatant was collected and transferred to a new tube containing freshly washed TiO 2 beads for a second incubation. The TiO 2 beads were collected separately and gently washed with 200 μl of loading buffer, 200 μl 80% acetonitrile/2% trifluoroacetic acid, 200 mM ammonium glutamate, and 200 μl of 50% acetonitrile/1% trifluoroacetic acid, respectively. The TiO 2 beads were dried and bound peptides were eluted sequentially in 10 min at first with 50 μl of 10% ammonium hydroxide, pH 11.7, then with 50 μl of 15% ammonium hydroxide/60% acetonitrile, and finally with 50 μl of 1% pyrrolidine. Eluted peptides were acidified by adding 75 μl 50% formic acid and cleaned up using OMIX C18, 10 µl SPE tips (Agilent).

LC-MS analysis
The samples were solved in 10 µl 0.1% formic acid and 5 µl was analyzed by LC-MS using a timsTOF Pro (Bruker Daltonik, Bremen, Germany), which was coupled online to a nanoElute nanoflow liquid chromatography system (Bruker Daltonik) via a CaptiveSpray nanoelectrospray ion source. The peptides were separated on a reversed-phase C18 column (25 cm × 75 µm, 1.6 µm, IonOpticks; Fitzroy, VIC, Australia). Mobile phase A contained water with 0.1% (vol/vol) formic acid, and acetonitrile with 0.1% (vol/vol) formic acid was used as mobile phase B. The peptides were separated by a gradient from 0 to 35% mobile phase B over 25 min at a flow rate of 300 nl/min at a column temperature of 50°C. MS acquisition was performed in DDA-PASEF mode.

ERK3-dependent tumor cell motility in vivo
All in vivo experiments were performed in accordance with the Swiss animal welfare ordinance and approved by the cantonal veterinary office Basel-Stadt. Female NSG mice were maintained in the Department of Biomedicine animal facilities in accordance with Swiss guidelines on animal experimentation (license number 2464). NSG mice are from in-house colonies. Mice were maintained in a sterile-controlled environment (a gradual light-dark cycle with light from 7:00 to 17:00, 21-25°C, 45-65% humidity).

Intravital imaging
Mice were anesthetized with Attane Isofluran (Provet AG) and anesthesia was maintained throughout the experiment with a nose cone. Tumors were exposed by skin flap surgery on a Nikon Ti2 A1plus multiphoton microscope and imaged at 880 nm with an Apochromat ×25/1.1 NA water immersion objective at a resolution of 1.058 µm per pixel. Cell motility was monitored by time-lapse imaging over 30 min in 2 min cycles, where a 100 mm Z-stack at 5-mm increments was recorded for each field of view starting at the tumor capsule. Three-dimensional time-lapse videos were analyzed using ImageJ. Images were registered using demon algorithm in MATLAB to correct for breathing movement (Kroon, 2023). Tumor cell motility was quantified manually. A tumor cell motility event was defined as a protrusion of half a cell length or more over the course of a 30 min video.

Statistical analyses
All experiments were repeated at least three times and exact replicate number (n) is specified for each figure. GraphPad Prism 9 was used to analyze data. Analyses performed for each data set, including statistical test and post-test, were specified for each figure. Where applicable, data are presented as mean ± SEM of at least three independent experiments. Significance levels are displayed in GP (GraphPad) style: *p<0.0332, **p<0.0021, ***p<0.0002, ****p<0.0001.