# Actomyosin pulsation and flows in an active elastomer with turnover and network remodeling

Deb Sankar Banerjee, Akankshi Munjal, Thomas Lecuit, Madan Rao

► To cite this version:

Deb Sankar Banerjee, Akankshi Munjal, Thomas Lecuit, Madan Rao. Actomyosin pulsation and flows in an active elastomer with turnover and network remodeling. Nature Communications, 2017, 8 (1), 10.1038/s41467-017-01130-1. hal-01691521

HAL Id: hal-01691521

https://hal.science/hal-01691521

Submitted on 12 Apr 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
ARTICLE

DOI: 10.1038/s41467-017-01130-1

OPEN

# Actomyosin pulsation and flows in an active elastomer with turnover and network remodeling

Deb Sankar Banerjee<fnref n="1" />, Akankshi Munjal², Thomas Lecuit² & Madan Rao³

Tissue remodeling requires cell shape changes associated with pulsation and flow of the actomyosin cytoskeleton. Here we describe the hydrodynamics of actomyosin as a confined active elastomer with turnover of its components. Our treatment is adapted to describe the diversity of contractile dynamical regimes observed in vivo. When myosin-induced contractile stresses are low, the deformations of the active elastomer are *affine* and exhibit spontaneous oscillations, propagating waves, contractile collapse and spatiotemporal chaos. We study the nucleation, growth and coalescence of actomyosin-dense regions that, beyond a threshold, spontaneously move as a spatially localized traveling front. Large myosin-induced contractile stresses lead to nonaffine deformations due to enhanced actin and crosslinker turnover. This results in a transient actin network that is constantly remodeling and naturally accommodates intranetwork flows of the actomyosin-dense regions. We verify many predictions of our study in *Drosophila* embryonic epithelial cells undergoing neighbor exchange during germband extension.

<footnote n="1">¹Raman Research Institute, Bangalore 560080, India. ²IBDM, UMR7288, CNRS-Aix Marseille Universite, Campus de Luminy case 907, 13009 Marseille, France and College de France, Paris, France. ³Simons Centre for the Study of Living Machines, National Centre for Biological Sciences (TIFR), Bangalore 560065, India. Correspondence and requests for materials should be addressed to M.R. (email: madan@ncbs.res.in)</footnote>
Tissue remodeling in diverse developmental contexts such as apical constriction in *Drosophila*¹⁻⁴, or *Caenorhabditis elegans*⁵, cell intercalation during *Drosophila* germband extension⁶⁻⁹, and xenopus extension and convergence¹⁰ have been observed to be associated with pulsation and flows of the medial actomyosin cytoskeleton. For instance, tissue extension in the *Drosophila* embryo proceeds by the intercalation of cells, a so-called T₁-process, initiated by the active shrinkage of a subset of cell junctions called 'vertical junctions' (aligned with the dorsal-ventral axis of the embryo, see Supplementary Fig. 1). The junctional shrinkage events are associated with medial-apical actomyosin pulsation and subsequent flow towards this vertical junction⁶, ⁸, ¹¹. In this paper, we develop a general theory for the dynamics of such spontaneous actomyosin pulsation and symmetry breaking flows, which should be applicable not only in the context of germband elongation, but also during other morphogenetic events. In addition, using germband cells in the *Drosophila* embryo as a model system, we provide experimental justification for the assumptions underlying our theoretical framework and verify many of its key predictions.

The apically located cortical actomyosin cytoskeleton comprises Myosin-II minifilaments, which bind onto a crosslinked actin filament network. The actin mesh is connected to E-cadherin adhesion molecules at the cell junctions via molecular linkers such as α-catenin an actin-binding protein⁵, ¹²⁻¹⁴ and β-catenin, which binds α-catenin and E-cadherin. Here we model the medial actomyosin mesh as an active elastomer embedded in a solvent, subject to active contractile stresses arising from the binding of myosin minifilaments (Fig. 1a)⁵, ⁷, ¹⁵, ¹⁶, and turnover of all components. Viewing the spatially resolved time-lapse images and movies (Supplementary Fig. 1 and Supplementary Movie 1) showing the diversity of dynamical regimes, including nucleation, growth, coalescence and flow of clusters of labeled myosin toward the cell junction in a *Drosophila* germband cell, should immediately convince one of the need for adopting a hydrodynamic approach to describe the spatiotemporal evolution of actomyosin densities.

The hydrodynamic equations for this active elastomer are derived from very general arguments based on symmetry considerations and conservation laws, and include minimal

phenomenological inputs. When the local contractile stress generated by bound myosin is not too large, the deformations are affine. In this regime, we describe the hydrodynamic modes without inertia and obtain phase diagrams by numerically solving the hydrodynamic equations. Similar equations have been written down in¹⁷, ¹⁸, where the primary focus has been on linear analysis, together with an analysis of the leading nonlinear effects. We find that the active affine elastomer exhibits spontaneous oscillations, contractile instabilities and coarsening of clusters enriched in actomyosin. The coarsening often leads to stable actomyosin-dense clusters that, beyond a threshold, acquire a polarity. This results in spontaneous movement as a spatially localized traveling front. Such localized traveling front solutions appear in other excitable systems such as the FitzHugh-Nagumo model¹⁹, to which our affine model bears a close resemblance. On the other hand, large myosin-induced contractile stresses can lead to nonaffine deformations due to actin turnover or network rupture. This results in a transient actin network that exhibits large intermittent strain fluctuations and intranetwork flow of the actomyosin-dense regions as a consequence of filament unbinding and rebinding. Interestingly, both the affine and nonaffine theories predict that the driving force for spontaneous movement comes from the actomyosin-dense region itself and not the cell boundary—we provide robust experimental verification of this. Our hydrodynamic analysis of the affine elastomer could be viewed as being closely related to the spring model of²⁰, though it goes far beyond this in the nonlinear analysis of the traveling front. Our analysis of the nonaffine elastomer, could be viewed as an active generalization of physical gels²¹ driven by local contractile force dipoles. Our general perspective reveals several significantly new aspects and provides a fresh conceptual understanding of this ubiquitous phenomenon.

## Results

**Active elastomer with turnover and network remodeling.** As just stated, we model the medial actomyosin crosslinked mesh as an active elastomer embedded in a solvent, subject to active contractile stresses arising from the binding of myosin minifilaments (Fig. 1a)⁵, ⁷, ¹⁵, ¹⁶, and turnover of all components. We will

<figure><img src="image_2.png" /><figcaption>**Fig. 1** Schematic of apical cortex and linear stability results: a Schematic showing the medial actomyosin cytoskeletal meshwork within the apical region of a cell belonging to the tissue. The actin filaments are attached to the cell junctions via E-cadherin (red dots). Myosin minifilaments bind (unbind) with rates $k_b(k_μ)$ and when bound, apply contractile stresses on the actin filament meshwork—the red circle demarcates a region of higher mesh compression. Both actin filaments and myosin minifilaments undergo turnover. b (I-II) Linear stability phase diagrams in (I) effective elastic stress density vs. contractile stress density at $k=1$ and (II) Effective contractile stress density vs. inverse lifetime of bound myosin at $B=4$. The stresses are normalized by the frictional stress density, $\Gamma k_b l^2$. The phases are described in the legend. Rest of the dimensionless parameters are $\alpha = 0.1$, $c = 0.1$, $D = 0.1$ (see “Methods” section)</figcaption></figure>
see that this allows us to span the diverse dynamical regimes exhibited by this system. Below, we provide a justification for treating this as an active elastomer with turnover, rather than directly as an active fluid, as was done in the context of C. elegans embryo<sup>22, 23</sup>.

We first note that, when the time scales of turnover of the actin mesh or crosslinkers are much smaller than the time scales of macroscopic processes of interest, then the stress relaxation is set by the turnover time scale, and the system should be treated as a fluid. Alternatively, when the actin mesh or crosslinker turnover time scales are much larger than the time scales of macroscopic processes, then the system should be treated as a dissipative elastic medium. The interpretation in these two extreme regimes is unambiguous. However, Fluorescence recovery after photo-bleaching(FRAP) measurements in germband cells, reveal an actin turnover time of around 10–20 s<sup>6</sup>, which is the same order as the time scale of nucleation and growth, around 10–20 s, and the period of medial actomyosin pulsation, around 50–100 s<sup>6, 8, 11</sup>. Indeed a more detailed look at the movies of myosin (Supplementary Movie 2), show a distribution of time scales and a range of dynamical regimes, starting from the appearance and disappearance of small myosin-rich speckles over a time scale of 5–10 s. This suggests that the appropriate description should span the short-time elastomeric and the longer time fluid-like regimes. As mentioned in the Introduction, we do this by starting out with an elastic description of a mesh, where the local deformations induced by myosin binding and release are affine. Increased myosin binding can lead to rapid turnover of actin and crosslinkers, resulting in loss of network integrity and its fluidization via intranetwork flows; we study this crossover to a nonaffine regime. Indeed, this strategy of going from the short-time elastic to long-time fluidization, as a function of increasing turnover, opens up novel rheological possibilities, such as correlated strain fluctuations and power-law response, which we take up later.

There are several other empirical reasons for starting with an elastomeric description and allowing for turnover of components:

1. Current high-resolution images of medial actin filaments, both in the germband cells and in the amnioserosa, show two distinct populations of actin filaments—a cell-spanning actin filamentous network which appears connected to the cell boundary and a possibly more rapidly turning over pool of shorter actin filaments (T. Lecuit and B. Dehapiot, unpublished observations). The relative levels of these two architectures of actin filaments is likely to be context dependent.

2. The medial actin mesh is connected to the cell junctions; consistent with this, modulating the strength of the coupling of medial actin to the junction via β-catenin affects the pulsation<sup>8</sup>.

3. Finally, the pulsation of medial actomyosin in these systems appears to be correlated with the oscillations in the area of the apical surface<sup>1, 8, 9, 20</sup>.

With this in mind, we shall start our discussion by describing the hydrodynamics of an affine elastomer subject to myosin turnover. We will next explore the consequences of actin turnover and network remodeling, the hydrodynamics of a nonaffine elastomer, giving rise to intranetwork flows and consequent fluid-like behavior.

**Affine elastomer hydrodynamics.** The hydrodynamic variables of the confined active elastomer (Fig. 1a) are the actin mesh displacement field $u$, the actin mesh density $\rho$, the density of the bound myosin minofilaments $\rho_b$ (we assume that the unbound

myosin constitutes a bath), the junctional E-cadherin density $\rho_c$ and the hydrodynamic velocity $v$. In this section we restrict ourselves to deformations that are affine, i.e., homogeneous deformations over a spatially coarse-grained scale. In the absence of flow and activity, the system relaxes to equilibrium governed by a free-energy functional $F[u, \rho, \rho_b, \rho_c] = \int d\mathbf{r}f_B$, where $f_B$ is the bulk elastic free-energy density of the elastomer written in terms of the linearized strain tensor $e = \frac{1}{2}(\nabla u + (\nabla u)^T)$ (see “Methods” section). In principle, we need to include a boundary contribution reflecting the soft anchoring arising from the actin-cadherin linkage; for the present, we simply include $\rho_c$'s contribution in the elastic moduli and hard boundary conditions on the dynamical equations.

In the presence of active processes, the so-called Rouse dynamics<sup>24</sup> without inertia is described by (Supplementary Note 1),

$$ \Gamma \dot{u} = \nabla \cdot (\sigma^e + \sigma^a + \sigma^d) \quad (1) $$

$$ \dot{\rho}_b + \nabla \cdot (\rho_b \dot{u}) = D \nabla^2 \rho_b + S_m \quad (2) $$

$$ \dot{\rho} + \nabla \cdot (\rho \dot{u}) = M \nabla^2 \frac{\delta F}{\delta \rho} + S_a \quad (3) $$

Equation (1) is a balance between frictional force experienced by the mesh (with coefficient $\Gamma$) and the net force acting on the mesh, written in terms of the total stress $\sigma = \sigma^e + \sigma^a + \sigma^d$, a sum of the elastic stress, $\sigma^e = \frac{\delta F}{\delta e}$, dissipative stress due to the viscosity ($\eta$) of the elastomer network, $\sigma^d = \eta \nabla \dot{u}$, and the active stress $\sigma^a$. The active stress should depend on the density of myosin and actin - it should increase with increasing myosin density and then saturate - it is thus reasonable to propose that

$$ \sigma^a = -\zeta(\rho, \rho_b)\Delta\mu I = -\frac{\zeta_1 \rho_b}{1 + \zeta_2 \rho_b} \chi(\rho) \Delta\mu I, \quad (4) $$

where $\Delta\mu$ is the difference in chemical potential during ATP hydrolysis, $\chi(\rho)$ is some smooth, positive valued function of the mesh density, $I$ the identity matrix and the overall negative sign, with parameters $\zeta_1 < 0$<sup>25</sup> and $\zeta_2 > 0$, ensure that the active stress is contractile (this results in a local “negative pressure”, which draws in the surrounding material). Equation (2) describes the dynamics of bound myosin filament density from advection by the filament velocity $\dot{u}$ and turnover $S_m = -k_u(\epsilon)\rho_b + k_b\rho$ by (un) binding. We allow for a possible strain-induced unbinding with a Hill-form, $k_u(\epsilon) = k_{u0}\epsilon^{\alpha}\epsilon^{26, 27}$. The sign of $\alpha$ can be taken to be either positive or negative (Supplementary Fig. 3a): $\alpha > 0$ implies a local extension (compression) of the mesh will increase (decrease) the myosin unbinding, while $\alpha < 0$ implies a local compression (extension) of the mesh will increase (decrease) the myosin unbinding. The choice $\alpha = 0$ implies that the myosin unbinding rate is a constant, independent of mesh deformation. We thus cover all possibilities (for details, see Supplementary Note 3). Similarly, Eq. (3) describes the dynamics of the actin mesh density including advection, permeation (with mobility $M$) and actin turnover, $S_a$.

In our analysis, we have ignored the long-range hydrodynamic interactions of the fluid; this is suggested by experiments that show the movement of the actin mesh relative to the fluid does not advect particles suspended in the fluid unless they are bound to the actomyosin mesh<sup>6, 8</sup>. For completeness, however, we display the full hydrodynamic equations in Supplementary Note 1.

From Eqs. (1–3), we note that there is a separation of time scales, which correspond to the advection time scale, myosin
<figure><img src="image_3.png" /><figcaption>**Fig. 2** Phase diagram obtained from numerical solutions of Eqs. 7 and 8. **a** Effective elastic stress density vs. contractile stress density, with $k=0.1$. **b** Effective contractile stress density vs. inverse-lifetime, with $B=5$. The phases are (i) Stable (yellow), (ii) spontaneous Oscillatory (blue), (iii) spontaneous Moving (grey) and (iv) contractile Collapse (light-green). The corresponding kymographs of the bound myosin density (indicated by the color and symbol on the upper right corner) is shown in **c**, **d**, **f** and **h**, respectively. The regions marked violet and dark-green are the coexistence phases—the oscillatory-moving coexistence (open circle) and the collapse-moving coexistence (open square), with the corresponding kymographs shown in **e** and **g**, respectively. Apart from the new phases, the topology of the phase diagrams are roughly similar to the linear stability diagram (Fig. 1b), except for the upturn of the phase boundaries towards larger active stress in **b**, which arises from the nonlinear strain-dependent unbinding. Symbols are points at which numerical solutions have been obtained. Rest of the dimensionless parameters are, $\alpha=3$, $c=0.1$, $\chi(\rho_0)\zeta_1=-0.5$, and $\zeta_2=0.1$ (see “Methods” section)</figcaption></figure>

turnover time and actin turnover time. The affine description is valid when the actin turnover time is longer than the other times scales. We will consider the case when the actin mesh density is “fast” compared to the time scale associated with the Peclet number $P_c \equiv L_v V/M$ ($L_v$ is the characteristic scale on which the mesh velocity, of typical magnitude $V$, varies); in this limit, the mesh density is slaved to the local compressive strain, thus with $\rho = \rho_0 + \delta\rho$, we arrive at $\delta\rho \propto -\epsilon_{ii}$, i.e., a local compression of the mesh leads to an increase in actin density (see “Methods” section).

The conceptual features of the affine dynamics are captured by a simple scalar version with one-elastic constant, $\sigma^e = B\epsilon$. It helps to make the equations dimensionless by choosing length and time in units of the screening length $l = \sqrt{\eta/\Gamma}$ and inverse binding rate $k_b^{-1}$, respectively; in these units, $u/l \to u$, $\rho_b/\rho_{b0} \to \rho_b$, $B/\Gamma k_b l^2 \to B$, $\zeta_1 \Delta\mu \rho_{b0}/k_b \Gamma l^2 \to \zeta_1 \Delta\mu$, $D/k_b l^2 \to D$ and $k_{u0}/k_b \to k$ are dimensionless. See “Methods” for parameter values.

**Linear analysis.** For low levels of bound myosin, the mesh deformation is small. It is thus appropriate to analyze the linear stability about the unstrained, homogeneous steady state ($u=0$, $\rho_b = \rho_{b0}$). To this order, the active stress reduces to $\sigma^a = -\zeta_1 \Delta\mu (1 + \zeta' \rho) \rho_b$, with $\zeta_1 < 0$ for contractility and $\zeta' > 0$.

Keeping in mind that $-\zeta_1 \Delta\mu$, $B$, and $D$ are the dimensionless active stress, elastic modulus and myosin diffusion, respectively, our linear stability analysis demonstrates that: (i) when the active stress is smaller than the elastic stiffness, $-\zeta_1 \Delta\mu < (B + D(1 + \sqrt{k/D})^2)/2$, the elastic mesh is stable with dispersion $\omega \sim q^2$, (ii) as the contractile strength exceeds $(B + D(1 + \sqrt{k/D})^2)/2$, the elastomer undergoes unstable oscillations with an amplitude that increases exponentially with time, (iii) at a threshold boundary, the elastomer supports a traveling wave solution with a speed $v_c^* = \sqrt[4]{(kDB^2/4)}$, (iv) beyond an active stress $-\zeta_1 \Delta\mu \sim B$, the elastomer contracts indefinitely. The phase diagrams shown in Fig. 1b reflect these four phases. Since this kind of linear analysis appears in¹⁷, ¹⁸, we simply state our results here and refer to Supplementary Note 2 and Supplementary Fig. 2 for details of calculations and dispersion relations.

In addition, we find that the qualitative features of these transitions remain unaltered when the strain-dependent unbinding parameter $\alpha$ is varied, as long as $\alpha \le \alpha_{max}(B, \zeta_1 \Delta\mu)$. This is discussed in detail in Supplementary Note 3 and Supplementary Fig. 3.

In spite of the simplifying nature of the linear analysis, it captures some gross features and offers some useful hints: (a) although unstable, the oscillatory behavior or pulsation (which will be stabilized by nonlinearities) requires advection and myosin
turnover<sup>8, 28</sup>, (b) the contractile instability is promoted by lowering elastic stiffness (by reducing levels of β-catenin/cadherin)<sup>5, 7, 15</sup> or reducing myosin unbinding rates<sup>8, 28, 29</sup>. The linear analysis, however, fails to capture finite-amplitude sustained oscillations and the phenomenology of moving configurations of actomyosin.

**Leading order nonlinearities.** It is possible that the oscillatory instability obtained in the linear analysis is tempered by the nonlinear terms in the dynamics, arising from strain-dependent unbinding $k_u(\epsilon)$, advection and active stress. Indeed, we find that by taking into account the nonlinear contributions to leading order, we can recast the equations as an excitatory-inhibitory dynamical system that exhibit sustained spontaneous oscillations. To see this, we work in 1-dim, and after fourier transforming Eqs. (1),(2) in a finite domain [0, L] with Neumann boundary conditions, retain only the smallest wavenumber, a 1-mode Galerkin truncation, which accommodates the strain-dependent unbinding nonlinearity. The resulting coupled ODEs upon re-scaling describe a generalized van der Pol oscillator<sup>30</sup> with linear damping and cubic nonlinearities. This admits a limit cycle through a supercritical-Hopf bifurcation for $-\zeta_1\Delta\mu \ge \frac{B}{2} + \frac{1}{\pi^2}$—a signature of the appearance of sustained spontaneous oscillations. At the onset of bifurcation, we use a fluctuation analysis to obtain the time period of oscillation $T = 2\pi(\eta/k_{u0}(B + \zeta_1\Delta\mu\rho_{b0}))^{1/2}$. To see the effects of the advective nonlinearity, we need to extend the above mode-truncation analysis to 2-modes. The resulting 3-dimensional dynamical system exhibits, in addition to limit cycles, temporal chaos as seen by the algebraic decay of the power-spectrum, positive Lyapunov exponent and denseness of the Poincare section (D.B., manuscript in preparation). However, our attempt to include nonlinear effects using low order mode-truncation, fails to capture the phenomenology of moving configurations of actomyosin.

Note that wave-like dispersion relations obtained from this truncated model are not a response to an external perturbation, but are self-generated. This cell-autonomous behavior is consistent with observations of pulsatile dynamics in medial actomyosin<sup>6, 11</sup>.

**Full nonlinear analysis shows traveling front solutions.** We now study the full nonlinear theory, which includes the effect of the nonlinearity arising from the active stress. This gives rise to a new set of solutions, namely the spatially localized traveling front solutions. To see this, we perform a numerical analysis of the full nonlinear equations in 1-dim. We first Taylor expand $\chi(\rho)$ in Eq. (4) about $\rho_0$, the mesh density in the unstrained configuration, and recast the active stress as (see “Methods” section)

$$ \sigma^a = \frac{-\zeta_1 \Delta \mu \rho_b}{1 + \zeta_2 \rho_b} (\chi(\rho_0) - c \chi'(\rho_0) \epsilon + c^2 \chi''(\rho_0) \epsilon^2 + \dots), \quad (5) $$

where c is a positive constant (see “Methods” section). Separating out the terms dependent on $\epsilon$ and only on $\rho_b$, and combining the former with the elastic stress $\sigma^a = B\epsilon$ in Eq. 1, leads to an effective “elastic free-energy”,

$$ \Phi(\epsilon) = \frac{1}{2} K_2(\rho_b, \rho_0) \epsilon^2 + \frac{1}{3} K_3(\rho_b, \rho_0) \epsilon^3 + \frac{1}{4} K_4(\rho_b, \rho_0) \epsilon^4, \quad (6) $$

where $K_i$ ($i=1, 2, 3$) are density dependent coefficients (see “Methods” section), and the quartic term with $K_4 > 0$ ensures that the local compressive strain does not grow without bound, as a consequence of steric hinderance, filament rigidity or crosslinking myosin. The $\Phi(\epsilon)$ that emerges as a consequence of activity, has 3 new features: (i) for weak active contractile stress, the minima at $\epsilon=0$ gets shallower, indicating that the elastic stiffness B

decreases, (ii) as we increase the active stress, there appears another minimum at $\epsilon = \epsilon_0$ (iii) for large active stresses, the $\epsilon = 0$ state can be unstable, with the effective $B < 0$ (Supplementary Fig. 4). The final 1-dim equations of motion are given by,

$$ \Gamma \dot{u} = \partial_x \Phi'(\epsilon) + \partial_x \sigma^a(\rho_b) \quad (7) $$

$$ \dot{\rho}_b = -\partial_x(\rho_b \dot{u}) + D \partial_x^2 \rho_b + S_m(\epsilon, \rho_b), \quad (8) $$

where $\Phi' \equiv \frac{\delta\Phi}{\delta\epsilon}$, $\sigma^a(\rho_b) = \frac{-\zeta_1 \Delta\mu \rho_b}{1 + \zeta_2 \rho_b} \chi(\rho_0)$ and myosin turnover $S_m(\epsilon, \rho_b) = -k_{u0} e^{\alpha\epsilon} \rho_b + k_b(1 - c\epsilon)$.

These equations are numerically solved with either periodic or Neumann boundary conditions using a finite difference scheme (see “Methods” section). Initial conditions are small amplitude random fluctuations about the homogeneous unstrained state. The numerical phase diagram, displayed in Fig. 2 shows several new features compared to the linear phase diagram, which we discuss below.

**Steady-state phase diagram.** The two features that are expected to arise from nonlinear effects, namely, the tempering of the linear instabilities to obtain both finite-amplitude oscillatory and finite-amplitude contractile collapse phases at intermediate and high contractile stresses, respectively, show up in the steady-state phase diagram, Fig. 2a, b. The corresponding kymographs in the bound myosin density (Fig. 2d, h) show the appearance of these steady state at late times. The time development of configurations in these phases can be summarized as follows—starting from a generic state with small random fluctuations about the homogeneous unstrained state, the configuration quickly results in a spatially heterogeneous (un)binding of myosin filaments onto the actin mesh, transiently generating localized compression. This will increase the local concentration of actin, which in turn will facilitate more myosin recruitment and hence more compression. This local compression will be resisted by an elastic restoring force, and the resulting strain can lead to an enhanced myosin unbinding. If it does, this will lead to a relaxation of the compressed region, to be followed by another round of binding-compression-unbinding, leading to the observed oscillations. In this spontaneous oscillating phase, the frequency gets smaller with increasing the active stress or decreasing unbinding rate<sup>8</sup>. On the other hand, if myosin unbinding does not occur fast enough, the elastomer will undergo a contractile instability, to be eventually stabilized by nonlinear effects such as steric hinderance and filament rigidity.

Additionally, there is a wholly unexpected feature that emerges from a numerical solution of the full nonlinear equations. In the parameter regime between the oscillatory and the contractile collapse phases, there appears a moving phase (Fig. 2), where spatially localized actomyosin-dense regions (which we later identify as traveling fronts) spontaneously move to either the left or right boundary. In the regimes between the pure moving phase and the oscillatory and collapse phases lie the coexistence phases where the moving phase coexists with oscillations and collapse, respectively. The corresponding kymographs in the bound myosin density (Fig. 2e-g) show the appearance of these steady state at late times. We may understand the occurrence of these phase transitions using a simple argument based on the relative time scales of these dynamical events, as displayed in Supplementary Fig. 5.

Several qualitative assertions follow immediately from the affine theory, such as: (i) the existence of bounded (finite-amplitude) oscillations requires both strain-dependent unbinding and turnover of myosin, (ii) the coexisting oscillation-moving and collapse-moving phases cannot be obtained in the absence of
<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /><figcaption>**Fig. 3** Kymographs of bound myosin density from theory and experiment: **a, b** Kymograph of the spatial profile of bound myosin density from theory, shows nucleation and growth (0 < *t* < 0.15), followed by coalescence (0.15 < *t* < 0.7) and movement (*t* > 0.7). Here *B* = 6, −ζ₁Δμ = 5.5, *k* = 0.5, *α* = 1 and *D* = 0.15. Rest of the parameters as in Fig. 2. **c** Kymograph of the spatial profile of labeled myosin from experiment, shows nucleation and growth (0 < *t* < 5 s), followed by coalescence (5 < *t* < 10 s) and **d** eventual movement of the formed actomyosin-rich region (*t* > 10 s)</figcaption></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /><figcaption>**Fig. 4** Time evolution of ρ_b, the density profile of bound myosin in a cluster prior to movement: **a** Prior to movement the myosin density profile is seen to be symmetric, following which we compute the instantaneous left-right fluxes, *J*<sub>L,R</sub> of myosin drawn into it. **b** Time evolution of the algebraic sum *J*<sub>L</sub> + *J*<sub>R</sub>, shows that after a while, there develops a net flux in one direction as a precursor to the asymmetric traveling front. **c** Emergence of the asymmetric traveling front which moves towards the right in a shape-invariant manner. **d–f** Myosin intensity profiles from experiments, which shows how an initial stationary symmetric myosin profile at *t* = 0, finally evolves to an asymmetric profile at *t* = 4.5 s, which then travels to the right. The degree of shape asymmetry of the profiles is described by the skewness *S*, the standard error of mean(s.e.m.) reported is due to projection of the images to one dimension</figcaption></figure>
<figure><img src="image_14.png" /></figure>

**Fig. 5** Anatomy of the traveling front in the co-moving frame: **a** Spatial profile of excess bound-myosin density (black) and strain $\epsilon$ (red) profile, **b** Spatial profile of Myosin density and derivative of the effective elastic free energy $\Phi'(\epsilon)$ (red), **c** Spatial profile of myosin density and active force (red). Horizontal axis is distance from centre of mass position, $x_{cm}$. Here $B = 6.0$, $-\zeta_1\Delta\mu = 4.8$, $k = 0.2$, $\alpha = 1.0$ and $D = 0.25$. **d** Theory predicts that the traveling front velocity is proportional to the net active force integrated over the front profile across the moving front. We demonstrate this fact from a numerical solution of the dynamical equations by varying the parameters of the active stress (circle) and the elastic stress (triangle). for different values of $B$ (with $-\zeta_1\Delta\mu = 6.0$, $k = 0.2$ and $D = 0.15$ fixed) and $-\zeta_1\Delta\mu$ (with $B = 8.0$, $k = 0.2$ and $D = 0.15$ fixed). Rest of the parameters as in Fig. 2. The color bar shows the magnitude of these stresses in dimensionless units. (inset) Skewness of bound myosin profile in the traveling front vs. velocity, obtained by varying the contractile stress $-\zeta_1\Delta\mu$ from 3-4 shows a linear increase followed by saturation. The other dimensionless parameters are: $B = 4$, $k = 0.2$, $D = 0.1$, $\alpha = 1$ and $c = 0.1$. The error bars are calculated as s.e.m

strain-dependent unbinding, (iii) advection is a necessary condition for front movement. We now look more closely at the dynamics of actomyosin-dense regions and compare the results of the affine theory with early time dynamics of myosin dense regions in vivo.

**Nucleation, growth and coalescence.** In the moving regime, the effective ‘elastic free-energy’ functional $\Phi(\epsilon)$ develops a second minima at $\epsilon = \epsilon_0$ corresponding to a local compression due to contractility (Supplementary Fig. 4). Initiation of movement starts with the nucleation of actomyosin-dense regions, which grow and coalesce to form larger actomyosin-dense regions. This is best seen using a space-time analysis of Eqs. 7–8 with initial conditions stated above. Kymographs of the spatial profile of bound myosin density, calculated from the theory show nucleation and growth ($0 < t < 0.15$), followed by coalescence ($0.15 < t < 0.7$) and eventual movement ($t > 0.7$) (Fig. 3a, b). This space-time behavior accurately recapitulates the early time dynamics of medial myosin in vivo as seen from the experimental kymographs, Fig. 3c, d.

**Asymmetric actomyosin profile of traveling front.** Here we investigate the origins of the spontaneous movement of the actomyosin-dense region. We study the configuration of the

localized actomyosin-dense region just prior to movement, and find that it assumes a symmetric localized profile (Fig. 4a) within which the strain $\epsilon = \epsilon_0$ (the second minimum) where $\Phi'(\epsilon) = 0$. The active stress within this region is higher than outside, the resulting gradient in stress should induce inflowing myosin currents from either side of it. We verify this by monitoring the fluxes $J_L$ and $J_R$, coming from the left and right of this symmetric profile. Over time, owing to stochasticity either in the initial conditions or the dynamics, there is a net flux, $J_L + J_R$, from either the left or the right (Fig. 4b), leading an asymmetric profile (Fig. 4c), and hence a gradient of $\rho_b$ across the profile. This marks the onset of the traveling front. This feature also appears to be present in the early time dynamics of the moving myosin profiles observed in vivo, as seen in Fig. 4d–f.

The affine theory predicts that the myosin density profile is asymmetric and moves as a traveling front, with a constant velocity while maintaining its shape (as long as there are no further coalescence events). We confirm this using a variety of initial conditions of the bound myosin density $\rho_b$, including starting with a single symmetric gaussian profile. We analyze the asymmetric profile of the traveling front by transforming to the co-moving frame $x \pm vt$ (Fig. 5).

We find that within the traveling front, the strain takes a value slightly more compressed relative to $\epsilon_0$, the value of the strain at the second minima, where $\Phi' = 0$. The traveling front is stably
<figure><img src="image_15.png" /><figcaption>**Fig. 6** Possible mechanisms for the movement of actomyosin-dense structures in an active affine elastomer: **a** moving deformation of the actomyosin mesh without turnover, implying a traveling front, and **b** moving deformation of the actomyosin mesh with differential myosin binding unbind rates at the leading and trailing edges of the front</figcaption></figure>

compressed in a force-free state (Fig. 5a, b). The asymmetric myosin profile gives rise to a gradient in the active stress (Fig. 5c), which provides the propulsion force for the traveling front to move to the right in Fig. 5c.

**Dynamics of traveling front.** The dynamics of the traveling front that emerges from the affine theory is local, its propulsion is therefore independent of the boundary or the distance from the boundary. We calculate the velocity of the traveling front by integrating Eqs. (7 and 8) across the scale $\Omega$ of the traveling front in the co-moving frame. This leads to the formula, $v = \Gamma^{-1} \int_{\Omega} \partial_{\tau} \sigma^a \equiv \Gamma^{-1} \int_{\Omega} f_{act}$, which states that the velocity depends only on the shape asymmetry of the front; if the shape is maintained over time, then the velocity is a constant. In Fig. 5d, we plot the $f_{act}$ vs. the velocity and show that they are proportional to each other over a large range of active and elastic stresses. We will later make connection with experiments, where it is convenient to use the shape asymmetry or skewness of the myosin profile, $S = \frac{\int_{\Omega} (x-x_{CM})^3 \rho_b(x)}{\int_{\Omega} (x-x_{CM})^2 \rho_b(x)^{3/2}}$, as a proxy for the driving force. We find that the traveling front velocity increases with increasing skewness before saturating at larger velocities (inset of Fig. 5d).

**Affine theory predicts moving deformation.** It is important to note that the movement of the actomyosin-dense region arising from affine deformations of the active elastomer is a moving deformation of the actomyosin mesh, and once established, is not contingent on myosin turnover, as shown in Fig. 6a. One could also sustain a traveling front or moving deformation of the actomyosin mesh by ensuring a differential myosin binding and unbinding rates at the leading and trailing edges of the front, in a kind of treadmilling movement, Fig. 6b. None of these however is associated with mass flow of actin and myosin. Indeed, the dynamical equations describing the active affine elastomer, Eqs. 7–8, bear a close resemblance to the generalized FitzHugh-Nagumo model, an excitable system which is known to exhibit traveling front solutions¹⁹.

**Nonaffine elastomer with network remodeling and intranetwork flow.** As time progresses and the myosin-dense clusters grow, the local contractile stresses can become large. Such large contractile stresses may dramatically enhance actin turnover and crosslinker unbinding. In Fig. 7a, we give a schematic plot of the

actin turnover time as a function of contractile stress, and compare it with the time scales of oscillation, front propagation and contractile collapse obtained from the affine theory. Indeed there is experimental evidence that actin turnover times at first increase with increasing contractility³¹⁻³³ - this is likely due to the fact that bound myosin might occupy the binding sites of actin remodeling proteins such as cofilin. There is also evidence that for large levels of contractility, the actin turnover rates are large, possibly because of destabilization of actin or unbinding of crosslinkers at large contractile stress. We make the plausible assumption that the qualitative change in turnover rates is sudden at a stress scale $\sigma^*$ (Fig. 7a), beyond which the deformation of the mesh can no longer be considered affine.

We use this insight to arrive at a nonaffine description of the active elastomer. Consider a disordered mesh comprising actin filaments linked to each other by crosslinkers such as $\alpha$-actinin and myosin (we will assume that this is an unentangled network). The bound myosin locally compresses the mesh here and there, recruiting more myosin in the process. When the local bound-myosin concentration goes beyond a threshold (so that the configuration now samples the second minimum of the effective free energy, $\Phi(e)$), the local compression is high and the mesh surrounding this myosin-dense region gets significantly stretched. This could lead to a tearing or ripping of the mesh, either via the unbinding of crosslinkers, destabilization of actin or by the sliding and slipping of filaments past each other. This mesh breakage subsequently healing by the rebinding of crosslinkers or actin itself.

With this picture in mind, we refer to several seminal studies on thermally activated reversibly crosslinked networks in the context of the dynamic properties of physical gels²¹, ³⁴, ³⁵. The most dramatic feature of such reversible networks is its internal fluidity, where each chain can diffuse through the entire network due to the finiteness of the crosslinker life time, in spite of being partially connected to the macroscopic network structure in the course of movement. These systems thus flow under an external stress on time scales longer than the crosslinker dissociation time.

In the context of actomyosin networks, the disrupting influence of filament turnover and the ultimate fluidization has been the subject of some study³⁶. In a very recent submission, the healing effect of turnover has also been investigated, and long-range network flows have been demonstrated in simulations of a model actomyosin cortex³⁷.

In the present context, when the local bound myosin concentration rises beyond a threshold and attains an asymmetric
<figure><img src="image_16.png" /><figcaption>**Fig. 7** Enhanced actin turnover results in nonaffine deformation of the elastomer: **a** Qualitative behavior of the actin turnover time as a function of contractile stress, which is consistent with data from refs. 31-33. This should be compared with the time scales of oscillation, front propagation and contractile collapse obtained from the affine theory (see, Supplementary Fig. 5). Based on the discussion in the text, we have placed the crossover stress $\sigma^*$ in the moving regime, thus implying that the crossover to the nonaffine description occurs in this regime. **b** Schematic showing the intranetwork flow of an actomyosin-dense region (enclosed within the yellow circle) resulting from active stress induced unbinding and rapid turnover of the actin in a transient actomyosin network</figcaption></figure>

density profile, it induces mesh breakage in its surrounding regions. This actomyosin-dense region can move through the entire network due to the finiteness of the crosslinker life time, in spite of being partially connected to the macroscopic network structure in the course of movement. These systems should thus exhibit flow under an internal active contractile stresses on time scales longer than the crosslinker dissociation time. This is depicted in Fig. 7b.

To describe this mathematically, it is convenient to define physical quantities coarse-grained over the scale of the actomyosin-dense region $\Omega$, such as $\bar{\rho}_b = \Omega^{-1} \int_{\Omega} \rho_b$ and $\mathbf{p} = \Omega^{-1} \int_{\Omega} \nabla \cdot \sigma^a$, the net force-dipole associated with the anisotropy of the myosin profile. We may decompose myosin density configuration $\bar{\rho}_b$ into a sum of actomyosin-dense clumps (contributing to nonaffine deformations) with volume fraction $\phi$ and a background (contributing to affine deformations). The equation for the $\bar{\rho}_b$ may now be written as,

$$ \bar{\rho}_b = -\nabla \cdot (\bar{\rho}_b \mathbf{v}) + D(1 - \phi)\nabla^2 \bar{\rho}_b + S_m(\bar{\rho}_b), \quad (9) $$

where, $\mathbf{v} = \phi\beta(\mathbf{p} + \gamma\mathbf{p} \cdot \mathbf{e}) + (1 - \phi)\dot{u}$, and $\beta$ is the mesh breakage probability and $\gamma$ is a strain-alignment parameter.

Flow described here is a consequence of internal active deformations in a transient actomyosin network. We believe this is a new physical phenomenon; unlike the flow observed in

physical gels under external load, this active flow is generated by internally generated stresses even in the absence of an external load.

In general, the interplay between actively generated contractile stresses and stress-dependent turnover of components that shows up in many cellular contexts, promises a rich phenomenology with novel rheological consequences. The constant remodeling and turnover of the actin mesh could drive the system from an elastic regime to a fluid-like behavior via a critical elastic state, which is characterized by correlated strain fluctuations, which might be intermittent. A more complete hydrodynamic theory of nonaffine deformations of a actomyosin network with turnover is a task for the future.

**Comparison with experiments.** In making comparisons of the theory presented here with experiments in germband cells in vivo, it is important to demarcate the affine and nonaffine regimes of the elastomer. The difficulty in doing this, is that, we do not know the detailed physical mechanism that would allow us to compute the form of the actin turnover time as a function of contractile stress (Fig. 7a). Nevertheless, it is clear that starting from an unstrained elastomer, the early time dynamics should be described by the linear and leading nonlinear analysis of the affine theory, as described earlier. The dynamical behaviors described by this analysis include oscillations and contractile collapse, from
<figure><img src="image_17.png" /><figcaption>**Fig. 8** Spatial profile and movement of myosin dense regions: Excess bound myosin density in the moving frame. **a–c** Three separate examples of the spatial profile of myosin intensity in the co-moving frame of the flowing actomyosin-dense cluster displayed at different times. Note that the actomyosin-dense cluster hardly changes its shape as it flows towards a junction</figcaption></figure>

<figure><img src="image_18.png" /><figcaption>**Fig. 9** Velocity and profile of myosin dense regions: **a** Histogram of the number of flows that move to the right junction *R* starting from either the right *R* (*P*(*R*IR*) or the left *L* (*P*(*R*IL*) region of the cell (inset shows schematic). The data collected from 24 actomyosin-dense clusters over 18 cells. The fact that the histograms are similar is consistent with the theoretical prediction that the flow is spontaneous and not driven by the cell boundaries. **b** Velocity of an isolated flowing actomyosin-dense cluster monitored over time shows that it is a constant, as predicted by theory. **c** Here we present average skewness values of 28 different pulses plotted against average velocities of the respective pulses. The linear behavior is quite clear from the data and a linear fit produced a slope value of 2.6 ± 0.5. This intensity data set was prepared from 20 myosin-GFP-tagged germband cells (N = 20). Error bars indicate s.e.m. of fluctuation of skewness and velocity values of a pulse in time</figcaption></figure>

which we make the following qualitative assertions: (i) the existence of bounded (finite-amplitude) oscillations requires advection, strain-dependent unbinding and turnover of myosin, (ii) the coexisting oscillation-moving and collapse-moving phases cannot be obtained in the absence of strain-dependent unbinding. As the kymographs in Fig. 2 show, these dynamical behaviors start emerging at $t \approx 2$ – which translates to a real time of $\approx 6$ s (see “Methods” section), which is smaller than the actin turnover time of around $10–20$ s⁶.

The stable, oscillatory and the collapse phases of the affine theory compare favorably with in vivo experiments in germband cell⁸. First, our phase diagram showing the stable, oscillatory and contractile collapse phase is grossly consistent with the experimental phase diagram that appears in ref. ⁸. Pulsatory solutions are obtained over a wide range of parameters, which include the pure oscillatory and the coexistence phases (Fig. 2d–g). In addition, we compare the finer aspects of the oscillatory phase with the pulsation seen in experiments (Supplementary Fig. 10). Consistent with ref. ⁸, we see that advection is crucial to obtain oscillations of bound myosin, both locally (Supplementary Fig. 10a, b) and cell-averaged (Supplementary Fig. 10(c)). We find a strong correlation between convergent (divergent) advection velocities and increased (decreased) myosin density (Supplementary Fig. 10a–c). Moreover, the amplitude of the oscillation decreases when we reduce actin density $\rho_0$ (Supplementary Fig. 10d), which is consistent with the actin perturbation experiments in ref. ⁸.

Now the nucleation, growth, and coalescence events of the myosin-dense clusters, precursors to the eventual moving phase of the affine theory, also take place over these time scales and compare well with the in vivo kymographs (Fig. 3). The establishment of the asymmetric profile of the myosin-dense cluster at the onset of movement also appears to be present in the early time dynamics of the moving myosin profiles observed in vivo, Fig. 4d–f.

The moving (traveling front) solution is the one significant prediction of the full nonlinear affine theory, and appears at later times. By this stage, the local contractile stresses are high, opening up the possibility of a crossover to the nonaffine regime. Both the affine and nonaffine theories predict movement—the affine theory predicts a moving deformation or a traveling front, while the movement in the nonaffine theory is associated with mass flow as a consequence of the steady turnover of actin filaments.

To test which of these pictures is true in vivo, we appeal to FRAP experiments performed on a small region within the actomyosin-dense cluster⁶. These preliminary experiments show loss of recovery upon FRAP, suggesting the possibility of actual mass flow, consistent with the predictions of the nonaffine theory (see Supplementary Information in ref. ⁶).

The fact that the actomyosin-dense region commences to move with a constant velocity when it has attained an asymmetric profile, is common to both the affine and nonaffine theory. It moves with a constant velocity in the direction where the leading edge has the smoother slope and maintains its asymmetric shape
as it moves. This appears to be consistent with the situation in vivo, Fig. 8a-c.

The driving force for this movement is established within the medial actomyosin-dense cluster and not the cell boundary. Its propulsion is therefore independent of the boundary or the distance from the boundary. We ask whether this is true of the moving actomyosin-dense regions in vivo. Figure 9a shows that irrespective of its initial position at the commencement of the flow, the moving actomyosin-dense region travels to the left or right cell boundary with equal probability. Further, Fig. 9b shows that the moving actomyosin-dense region travels with a constant velocity as it moves towards a given cell boundary, its speed does not depend on the distance from the cell boundary.

From the myosin intensity, we compute the shape asymmetry via the skewness, $S = \frac{\int_{\Omega} (x-x_{CM})^3 \rho_b(x)}{\int_{\Omega} (x-x_{CM})^2 \rho_b(x)^{3/2}}$, of the myosin profile, and find that it is proportional to the speed of the moving myosin-dense cluster (Fig. 9c).

To our mind, this establishes unambiguously that the flow toward the junctions is spontaneous with the driving force coming from the gradient in myosin established within the front. The boundary does not affect the flow speed, at best, weak asymmetries that may arise at the boundary (for instance due to an asymmetry in functional cadherin) may bias the direction of the flow$^7$.

One consequence of this is that a cluster moving to the right, might reverse its direction following a coalescence with another moving cluster; such reversals are observed in in vivo experiments (Supplementary Fig. 7 and Supplementary Fig. 9).

## Discussion

Although simplified, in that it has completely ignored the coupling of actomyosin dynamics to local chemical signaling such as Rho$^8$, we believe this active elastomer model with strain-dependent turnover of components, admitting both affine and nonaffine deformations, captures the essential physics of actomyosin pulsation and flows observed in a wide variety of tissue remodeling contexts such as *Drosophila* germband extension and dorsal closure in the amnioserosa. The minimal ingredients for actomyosin pulsation and flow are mesh-elasticity, actomyosin contractility, advection and turnover of both myosin and actin.

In this study, we have modeled the medial actin mesh during apical constriction and germband extension in the *Drosophila* embryo as an active elastomer embedded in a solvent$^5, 7, 15, 16$, which undergoes turnover of all its components. Our description goes from the hydrodynamics of an affine elastomer to a nonaffine elastomer, the latter incorporating network rupture and remodeling, resulting in intranetwork flows. It thus goes from an elastomer to a fluid-like description. Together, the affine and nonaffine active elastomer model captures the range of dynamical regimes exhibited in this system.

To make a detailed comparison of the spatiotemporal actomyosin patterns with experiments generated using quantitative imaging, we will need to extend this numerical study to 2-dimensions, using appropriate (anisotropic) boundary conditions and allowing for shear. The nucleation and growth of the actomyosin-rich domains are similar to that seen in 1-dim, with the difference being that domains can move around each other in 2-dim and can exhibit anisotropic movement.

In future, we would like to extend this framework to understand the dynamical coupling of the medial actomyosin with degrees of freedom (concentration of E-cadherin) attached to a deformable cell junction. Though the emergence of actomyosin flows does not depend on specific boundary conditions, cell boundaries may directionally bias the intrinsic ability of actomyosin networks to generate flow, as proposed before$^7$.

## Methods

**Free-energy functional.** The bulk elastic free-energy density of the elastomer in terms of the linearized strain tensor $\epsilon_{ij}$ is given by $f_B = \frac{\mu}{2}\epsilon_{ij}\epsilon_{ij} + \frac{1}{2}\epsilon_{ii}\epsilon_{jj} + \dots + C\delta\rho\epsilon_{ii} + \frac{1}{2}\delta\rho^2 + \dots$, where the ... represent higher-order terms to stabilize a possible activity induced contracted state and to prevent a runway density increase ($\delta\rho$ is the deviation of the local mesh density from its equilibrium value, $\delta\rho = -(C/A)\epsilon_{ii}$). It is this $C/A \equiv c > 0$ that makes its appearance in all the figure captions and in Eq. 5.

The coefficients that appear in the effective elastic free-energy functional $\Phi(\epsilon)$ (Eq. 6) conceived from renormalization of elastic stress by the active stress $-\zeta(\rho, \rho_b)$ $\Delta\mu I$ are given by, $K_2 = B + \frac{\zeta_1\Delta\mu\rho_b}{1+\zeta_2\rho_b}c\chi'(\rho_0)$, $K_3 = -\frac{\zeta_1\Delta\mu\rho_b}{1+\zeta_2\rho_b}c^2\chi''(\rho_0)$ and $K_4 = \frac{\zeta_1\Delta\mu\rho_b}{1+\zeta_2\rho_b}c^3\chi'''(\rho_0)$. Here $\zeta_1 < 0$ and $\chi'(\rho_0) > 0$, so $K_2$ goes from being positive to negative as contractility increases. The signs of the other constants are as follows: $\zeta_2 > 0, \chi''(\rho_0) > 0$ and $\chi'''(\rho_0) < 0$, so that $K_3$ and $K_4$ are always positive.

**Parameter values.** Here we relate the values of the dimensionless parameters to real units extracted from a variety of experimental measurements. For the unit of length, $l = \sqrt{\eta/\Gamma}$ or the actin mesh size, we take 0.5 µm (consistent with the rough estimates in ref. 6). The unit of time, $k_c^{-1}$, can be estimated from the myosin FRAP data$^8$, we find $k_u = 0.2 \pm 0.08$ s$^{-1}$, and taking the ratio $k = k_u/k_b$ to be 1.0, we obtain a binding rate, $k_b = 0.2$ s$^{-1}$. The viscosity of the mesh is taken to be 50 Pa s$^{36}$.

We can now convert all the dimensionless values into real values, and check for consistency with other experimental estimates. Thus, a dimensionless value of the bulk modulus $B = 5$ translates to $B = 42$ Pa (consistent with what can be estimated from ref. 38). Similarly, a dimensionless value of the magnitude of the active stress, $|\zeta_1\Delta\mu| = 5$ translates to $|\zeta_1\Delta\mu\rho_b| = 42$ Pa (roughly the order of magnitude estimated from ref. 38). Finally, the dimensionless diffusion coefficient $D = 0.25$ implies a real value of $D = 0.01$ µm$^2$/s.

This implies that a dimensionless front velocity $v = 1$ (Fig. 9c) translates to a real velocity of $v = 0.08$ µm/s (consistent with the flow velocity of actomyosin reported in ref. 6, also see Supplementary Fig.8). Likewise, the time period of oscillation obtained in the section on “Leading order nonlinearities”, $T = 2\pi(\eta/k_{u0}(B + \zeta_1\Delta\mu\rho_b))^{1/2}$ translates to be ~20 s and the active propulsion force $f_{act}$ is estimated at 30–60 pN (consistent with ref. 38).

**Numerical methods.** Since the dynamical equations have nonlinear advection and diffusion, care must be taken in evaluating the flux due to dissipative and dispersive errors arising from spatial discretization. We use the finite volume method for spatial discretization$^{39}$, which has been found to be useful for nonlinear advection equations$^{40}$.

We calculate the numerical flux using the Van-Leer’s flux limiter, which uses a different formula to calculate the spatial derivative depending on how sharply the $\rho_b$ profile changes in space. When the profile changes very fast, the scheme implements the upwind method$^{41}$, which reduces the dispersion error through numerical diffusion. When the profile changes smoothly the scheme implements a second order accurate method called Lax-Wendroff method$^{42}$. In our numerical scheme, the density ($\rho_b$ or $\rho$) flux on the interface between $i^{th}$ and $(i-1)^{th}$ node is computed as,

$$f_{i-\frac{1}{2}}^{n+\frac{1}{2}} = \frac{1}{2} v_{i-\frac{1}{2}} \left[ (1 + \theta_{i-\frac{1}{2}}) \rho_{i-1}^n + \left(1 - \theta_{i-\frac{1}{2}}\right) \rho_i^n \right] + \frac{1}{2} |v_{i-1}| \left(1 - \left|\frac{v_{i-1} + \Delta x}{\Delta x}\right|\right) \phi_{i-\frac{1}{2}}^n \left(\rho_i^n - \rho_{i-1}^n\right) \quad (10)$$

Here, $v_{i-\frac{1}{2}} = \frac{v_{i-1} - v_{i-2}}{2}$ and

$$\theta_{i-\frac{1}{2}} = \begin{cases} +1, & \text{if } v_{i-\frac{1}{2}} > 0 \\ -1, & \text{if } v_{i-\frac{1}{2}} \le 0 \end{cases} \quad (11)$$

$$r_{i-\frac{1}{2}}^{n} = \begin{cases} \frac{\rho_{i-1}^{n} - \rho_{i-2}^{n}}{\rho_{i}^{n} - \rho_{i-1}^{n}}, & \text{if } v_{i-\frac{1}{2}} > 0 \\ \frac{\rho_{i-1}^{n} - \rho_{i}^{n}}{\rho_{i}^{n} - \rho_{i-1}^{n}}, & \text{if } v_{i-\frac{1}{2}} \le 0 \end{cases} \quad (12)$$

The function $\phi(r)$ is the Van-Leer flux limiter

$$\phi(r) = \frac{r + |r|}{1 + |r|}. \quad (13)$$

The time integration is done with a total variation diminishing (TVD) 3rd order Runge-Kutta method$^{43}$. All other derivative terms were discretized using a simple finite difference. The initial conditions were chosen from a uniform random
distribution of fixed width about a uniform, unstrained configuration. We used periodic boundary conditions throughout and a time-space discretization, $\Delta t = 10^{-4}$ and $\Delta x = 5 \times 10^{-2}$.

**Experimental methods and data analysis.** For fluorescence time-lapse imaging, embryos at stage 7 were dechorionated with 100% bleach and mounted on No. 1 coverslip with halocarbon oil<sup>44</sup>. A 100×, 1.4 N.A oil immersion objective with Nikon spinning disc Eclipse Ti inverted microscope was used<sup>44</sup>. The system acquires images using the MetaMorph software. Starting from the most apical plane, 4–7 z-sections 0.5 µm apart were acquired every 1.5–4 s (depending on the experiment) using a single camera. Sum-intensity z-projection of slices was used for all quantifications, followed by a background subtraction using the available plugin in Fiji.

All the data analysis are done using customized code written in MATLAB. For realization of pulse shape, velocity and skewness pulses were projected on the line of their movement (for details see Supplementary Fig. 6).

**Data availability.** All the relevant data are available from the authors upon request.

Received: 5 April 2016 Accepted: 22 August 2017

Published online: 24 October 2017

## References

1. Martin, A. C., Kaschube, M. & Wieschaus, E. F. Pulsed contractions of an actin-myosin network drive apical constriction. *Nature* **457**, 495–499 (2009).

2. Solon, J., Kaya-Copur, A., Colombelli, J. & Brunner, D. Pulsed forces timed by a ratchet-like mechanism drive directed tissue movement during dorsal closure. *Dev. Cell* **13**, 1331–1342 (2009).

3. Meghana, C. et al. Integrin adhesion drives the emergent polarization of active cytoskeletal stresses to pattern cell delamination. *Proc. Natl Acad. Sci.* **108**, 9107–9112 (2011).

4. Blanchard, G. B., Murugesu, S., Adams, R. J., Martinez-Arias, A. & Gorfinkiel, N. Cytoskeletal dynamics and supracellular organisation of cell shape fluctuations during dorsal closure. *Development* **137**, 2743–2752 (2010).

5. Roh-Johnson, M. et al. Triggering a cell shape change by exploiting preexisting actomyosin contractions. *Science* **335**, 1232–1235 (2012).

6. Rauzi, M., Lenne, P. F. & Lecuit, T. Planar polarized actomyosin contractile flows control epithelial junction remodelling. *Nature* **468**, 1110–1114 (2010).

7. Levayer, R. & Lecuit, T. Oscillation and polarity of E-cadherin asymmetries control actomyosin flow patterns during morphogenesis. *Dev. Cell* **26**, 162–175 (2013).

8. Munjal, A., Philippe, J.-M., Munro, E. & Lecuit, T. A self-organized biomechanical network drives shape changes during tissue morphogenesis. *Nature* **524**, 351–355 (2015).

9. Fernandez-Gonzalez, R. & Zallen, J. A. Oscillatory behaviors and hierarchical assembly of contractile structures in intercalating cells. *Phys. Biol.* **8**, 045005 (2011).

10. Keller, R. et al. Mechanisms of convergence and extension by cell intercalation. *Phil. Trans. R. Soc. Lond* **355**, 897–922 (2000).

11. Rauzi, M., Verant, P., Lecuit, T. & Lenne, P. F. Nature and anisotropy of cortical forces orienting drosophila tissue morphogenesis. *Nat. Cell Biol.* **10**, 1401–1410 (2008).

12. Cavey, M., Rauzi, M., Lenne, P. F. & Lecuit, T. A two-tiered mechanism for stabilization and immobilization of E-cadherin. *Nature* **453**, 751–756 (2008).

13. Yonemura, S., Wada, Y., Watanabe, T., Nagafuchi, A. & Shibata, M. alpha-Catenin as a tension transducer that induces adherens junction development. *Nat. Cell Biol.* **12**, 533–542 (2010).

14. Buckley, C. D. et al. The minimal cadherin-catenin complex binds to actin filaments under force. *Science* **346**, 1254211 (2014).

15. Martin, A. C., Gelbart, M., Fernandez-Gonzalez, R., Kaschube, M. & Wieschaus, E. F. Integration of contractile forces during tissue invagination. *J. Cell Biol.* **188**, 735–749 (2010).

16. Wu, S. K. et al. Cortical F-actin stabilization generates apical-lateral patterns of junctional contractility that integrate cells into epithelia. *Nat. Cell Biol.* **16**, 167–178 (2014).

17. Banerjee, S., Liverpool, T. B. & Marchetti, M. C. Generic phases of cross-linked active gels: relaxation, oscillation and contractility. *Europhys. Lett.* **96**, 58004 (2011).

18. Banerjee, S. & Marchetti, M. C. Instabilities and oscillations in isotropic active gels. *Soft Matter* **7**, 463 (2011).

19. J. Keener and J. Sneyd, In *Mathematical Physiology, I: Cellular Physiology*, (2nd Ed., Springer, 2009).

20. Dierkes, K., Sumi, A., Solon, J. & Salbreux, G. Spontaneous oscillations of elastic contractile materials with turnover. *Phys. Rev. Lett.* **113**, 148102 (2014).

21. Tanaka, F. & Edwards, S. F. Viscoelastic properties of physically crosslinked networks: Transient network theory. *Macromolecules* **25**, 1516–1523 (1992).

22. Mayer, M., Depken, M., Bois, J. S., Julicher, F. & Grill, S. W. Anisotropies in cortical tension reveal the physical basis of polarizing cortical flows. *Nature* **467**, 617–621 (2010).

23. Goehring, N. W. et al. Polarizing PAR proteins by advective triggering of a pattern-forming system. *Science* **334**, 1137–1141 (2011).

24. Doi, M. & Edwards, S. F. In *The Theory of Polymer Dynamics* (Oxford University Press, 1986).

25. Marchetti, M. C. et al. Hydrodynamics of soft active matter. *Rev. Mod. Phys.* **85**, 1143 (2013).

26. Fernandez-Gonzalez, R. et al. Myosin II dynamics are regulated by tension in intercalating cells. *Dev. Cell* **17**, 736–743 (2009).

27. Kovacs, M., Thirumurugan, K., Knight, P. J. & Sellers, J. R. Load-dependent mechanism of nonmuscle myosin 2. *Proc. Natl Acad. Sci.* **104**, 9994–9999 (2007).

28. Vasquez, C. G., Tworoger, M. & Martin, A. C. Dynamic myosin phosphorylation regulates contractile pulses and tissue integrity during epithelial morphogenesis. *J. Cell Biol.* **206**, 435–450 (2014).

29. Kasza, K. E., Farrell, D. L. & Zallen, J. A. Spatiotemporal control of epithelial remodeling by regulated myosin phosphorylation. *Proc. Natl Acad. Sci.* **111**, 11732–11737 (2014).

30. Strogatz, S. In *Nonlinear Dynamics and Chaos* (2nd Ed., Perseus, 2014).

31. Murthy, K. & Wadsworth, P. Myosin-II-dependent localization and dynamics of f-actin during cytokinesis. *Curr. Biol.* **15**, 724–731 (2005).

32. Engl, W., Arasi, B., Yap, L. L., Thiery, J. P. & Viasnoff, V. Actin dynamics modulate mechanosensitive immobilization of E-cadherin at adherens junctions. *Nat. Cell Biol.* **16**, 584–591 (2014).

33. Guha, M., Zhou, M. & Wang, Y. Cortical actin turnover during cytokinesis requires myosin II. *Curr. Biol.* **15**, 732–736 (2005).

34. Murthy, A. K. & Muthukumar, M. Dynamic mechanical properties of poly (γ-benzyl L-α-glutamate) gels in benzyl alcohol. *Macromolecules* **20**, 564–569 (1987).

35. Leibler, L., Rubinstein, M. & Colby, R. H. Dynamics of reversible networks. *Macromolecules* **24**, 4701–4707 (1991).

36. Kim, T., Gardel, M. L. & Munro, E. Determinants of fluidlike behavior and effective viscosity in cross-linked actin networks. *Biophys. J.* **106**, 526–534 (2014).

37. McFadden, W. M., McCall, P. M. & Munro, E. M. Filament turnover is essential for continuous long range contractile flow in a model actomyosin cortex. Preprint arXiv:1612.07430 [q-bio.SC] (2016).

38. Bambardekar, K., Clement, R., Blanc, O., Chardes, C. & Lenne, P. F. Direct laser manipulation reveals the mechanics of cell contacts in vivo. *Proc. Natl Acad. Sci.* **112**, 1416–1421 (2015).

39. Sweby, P. K. High resolution schemes using flux limiters for hyperbolic conservation laws. *SIAM J. Num. Anal.* **21**, 995–1011 (1984).

40. R.J. LeVeque In *Numerical Methods for Conservation Laws*, Lectures in Mathematics (2nd edition, Birkhauser 1992).

41. Courant, R., Isaacson, E. & Rees, M. On the solution of nonlinear hyperbolic differential equations by finite differences. *Comm. Pure App. Math.* **5**, 243–255 (1952).

42. Lax, P. & Wendroff, B. Systems of conservation laws. *Comm. Pure App. Math.* **13**, 217–237 (1960).

43. Gottlieb, S. & Shu, C.-W. Total variation diminishing Runge-Kutta scheme. *Math. Comp* **67**, 73–85 (1988).

44. Cavey, M. & Lecuit, T. Imaging cellular and molecular dynamics in live embryos using fluorescent proteins, Meth. *Mol. Biol.* **420**, 219–238 (2008).

## Acknowledgements

We thank the CNRS (Laboratoire International Associé between IBDM and NCBS), AMIDEX and the Labex INFORM (ANR-11-LABX-0054), for support. A.M. and T.L. were supported by the ERC (Biomecamorph 323027). D.S.B. and M.R. thank K. Vijaykumar, S. Ramaswamy and R. Morris for extremely useful discussions. D.S.B. thanks R. Mandal, K. Husain and A. Das for numerous helpful discussions on numerics and data analysis.

## Author contributions

Theory done by D.S.B. and M.R. inspired by experiments done by A.M. and T.L. Numerics done by D.S.B. All authors contributed in writing the paper.

## Additional information

*Supplementary Information* accompanies this paper at doi:10.1038/s41467-017-01130-1.

**Competing interests:** The authors declare no competing financial interests.
**Reprints and permission** information is available online at http://npg.nature.com/
reprintsandpermissions/

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_19.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing,

adaptation, distribution and reproduction in any medium or format, as long as you give

appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2017