Leading edge maintenance in migrating cells is an emergent property of branched actin network growth
Abstract
Animal cell migration is predominantly driven by the coordinated, yet stochastic, polymerization of thousands of nanometerscale actin filaments across micronscale cell leading edges. It remains unclear how such inherently noisy processes generate robust cellular behavior. We employed highspeed imaging of migrating neutrophillike HL60 cells to explore the finescale shape fluctuations that emerge and relax throughout the process of leading edge maintenance. We then developed a minimal stochastic model of the leading edge that reproduces this stable relaxation behavior. Remarkably, we find lamellipodial stability naturally emerges from the interplay between branched actin network growth and leading edge shape – with no additional feedback required – based on a synergy between membraneproximal branching and lateral spreading of filaments. These results thus demonstrate a novel biological noisesuppression mechanism based entirely on system geometry. Furthermore, our model suggests that the Arp2/3mediated ~70–80° branching angle optimally smooths lamellipodial shape, addressing its longmysterious conservation from protists to mammals.
Editor's evaluation
This paper describes analysis and modeling of leading edge fluctuations in migrating cells driven by a branched Arp2/3 lamellipodial network. A stochastic model shows how branching contributes to shape stability, and reproduces the measured spectrum and dynamics of leading edge fluctuations. Analysis of the model as a function of branching angle suggests that the Arp2/3 branching angle might be selected to smooth lamellipodial shape. This work provides new ideas to a big field of research, including Fourier analysis of leading edge fluctuations.
https://doi.org/10.7554/eLife.74389.sa0eLife digest
In every human cell, there are tens of millions of proteins which work together to control everything from the cell’s shape to its behavior. One of the most abundant proteins is actin, which organizes itself into filaments that mechanically support the cell and help it to move.
These filaments are very dynamic, with individual actin molecules constantly being added or removed. This allows the cell to build large structures with distinct shapes and properties. Many motile cells, for example, have a structure called a lamellipodium which protrudes at their ‘leading edge’ and pushes them forward. The lamellipodium has a very robust shape that does not vary much between different cell types, or change significantly as cells migrate. But how the tens of thousands of actin molecules inside the lamellipodium organize themselves into this large, stable structure is not fully understood.
To investigate, Garner and Theriot used highspeed video microscopy to track the shape of human cells cultured in the laboratory. As the cells crawled along a glass surface, their leading edge undulated like strings being plucked on a guitar. A computer simulation showed that these ripples can be caused by filaments randomly adding and removing actin molecules.
While these random movements could destabilize the structure of the leading edge, the simulation suggests that another aspect of actin filament growth smooths out any fluctuations in the lamellipodium’s shape. Actin networks in the lamellipodium have a branched configuration, with new strands emerging off each other at an angle like branches in a tree. Garner and Theriot found that the specific angle in which new filaments are added smooths out the lamellipodium’s shape, which may explain why this geometry has persisted throughout evolution.
These findings suggest that the way in which actin filaments join together helps to maintain the shape of large cellular structures. In the future, scientists could use this design principle to build molecular machines that can selforganize into microstructures. These engineered constructs could be used to modulate the activity of living cells that have been damaged by disease.
Introduction
Cell migration driven by actin polymerization plays an essential role in countless organisms spanning the eukaryotic tree of life (Pollard and Cooper, 2009; FritzLaylin et al., 2017a; Welch et al., 1997). Across this broad phylogeny, cells have been observed to form a dizzying array of protrusive actin structures, each exhibiting unique physical and biological properties (Svitkina, 2018). In all cases, the fundamental molecular unit of these micronscale structures is the single actin filament, which polymerizes stochastically by addition of single monomers to push the leading edge membrane forward (Mogilner and Oster, 1996; Peskin et al., 1993; Theriot et al., 1992; Prass et al., 2006). Higher order actin structures, and the biological functions they robustly enable, are therefore mediated by the collective action of thousands of stochastically growing filaments (Svitkina, 2013). It remains an open question how cells control for – or leverage – this inherent stochasticity to maintain stable leading edge protrusions over length and time scales more than three orders of magnitude larger than the scales of actin monomer addition (Rafelski and Theriot, 2004; Vavylonis et al., 2005).
Perhaps the archetype of dynamically stable actin structures is the lamellipodium, a flat ‘leaflike’ protrusion that is ~200 nm tall, up to 100 µm wide, and filled with a dense network of dendritically branched actin filaments (Svitkina et al., 1997; Abraham et al., 1999; Laurent et al., 2005; FritzLaylin et al., 2017b). Cell types that undergo lamellipodial migration (most notably fish epidermal keratocytes and vertebrate neutrophils) can maintain a single, stable lamellipodium for minutes to hours, allowing the cells to carry out their biological functions (Tsai et al., 2019; Lacayo et al., 2007). For example, in their in vivo role as first responders of the innate immune system, neutrophils must undergo persistent migration over millimeterscale distances to reach sites of inflammation and infection (de Oliveira et al., 2016; Kolaczkowska and Kubes, 2013). Regardless of the cell type, the origins of this striking stability in the face of stochastic actin filament polymerization remain elusive. It has been widely been assumed for decades that some sort of regulatory or mechanical feedback mechanism must be required for lamellipodial shape stability, with extensive experimental efforts identifying membrane tension (DizMuñoz et al., 2016; Houk et al., 2012; Mueller et al., 2017; Gauthier et al., 2012; Tsujita et al., 2015; Batchelder et al., 2011; Sens and Plastino, 2015), plasma membrane curvaturesensing proteins (Tsujita et al., 2015; Pipathsouk et al., 2019), a competition for membraneassociated free monomers (Mullins et al., 2018), and forcefeedback via directional filament branching (Risca et al., 2012) as potential contributors. The stability of lamellipodia has also been theoretically proposed to depend on the dendritically branched structure of their actin networks, wherein filaments are oriented at an angle relative to the cell’s direction of migration, allowing growing filament tips to spread out laterally along the leading edge as they polymerize (Lacayo et al., 2007; Grimm et al., 2003). Although any acute angle would permit spreading, we note that filament orientation in cells has been experimentally observed to be highly stereotyped, averaging ±35° relative to the membrane normal (Maly and Borisy, 2001; Verkhovsky et al., 2003) – approximately one half of the highly evolutionarilyconserved ~70° branch angle mediated by the Arp2/3 complex (Mullins et al., 1998; Volkmann et al., 2001; Rouiller et al., 2008.) In contrast to the proposed stabilizing role of spreading, several other features of lamellipodial actin are known to impart nonlinearities on network growth, which might amplify stochastic fluctuations. For instance, dendritic branching is an autocatalytic process which can lead to explosive growth (Mullins et al., 2018; Carlsson, 2001). In addition, the growth rates of the actin network are dependent on the velocity of the flexible membrane surface it is pushing, in a manner which imparts hysteresis to the system (Mueller et al., 2017; Parekh et al., 2005). How spreading might interact with these complexities – and the ultimate consequences for maintenance of a stable leading edge – remains unknown.
Seeking to dissect the origins of lamellipodial stability, we pursued complementary experimental and computational methodologies. First, we performed highspeed, highresolution microscopy on migrating human neutrophillike HL60 cells to monitor their leading edge shape dynamics. In contrast to the remarkable overall lamellipodial stability observed over minutes, highspeed imaging revealed that the leading edge shape is extremely dynamic at shorter time and length scales, constantly undergoing finescale fluctuations around the average cell shape. We determined that these shape fluctuations continually dissipate (thereby enabling long time scale lamellipodial maintenance) in a manner quantitatively consistent with viscous relaxation back to the timeaveraged leading edge shape. We next developed a minimal stochastic model of branched actin network growth against a flexible membrane, broadly applicable to a wide variety of cell types, that was able to recapitulate the global leading edge stability and finescale fluctuation relaxation behavior observed in cells. Our model suggests that the suppression of stochastic fluctuations is an intrinsic, emergent property of collective actin dynamics at the leading edge, as branched network geometry alone is necessary and sufficient to generate lamellipodial stability. Moreover, we find that the evolutionarilyconserved geometry, the ~70° branching angle of the Arp2/3 complex, optimally quells shape fluctuations.
Results
Finescale leading edge shape fluctuations revealed at high spatiotemporal resolution
Neutrophils form lamellipodia that are intrinsically lamellar, maintaining a thin, locally flat sheet of actin even in the absence of support structures like the substrate surface (FritzLaylin et al., 2017b). Here, we study the migration of neutrophillike HL60 cells (Spellberg et al., 2005) within quasitwodimensional confinement between a glass coverslip and an agarose pad overlay (Millius and Weiner, 2009). In addition to serving as an excellent in vitro model for neutrophil surveillance of tissues, this assay allows for easy visualization and quantification of lamellipodial dynamics by restraining the lamellipodium to a single imaging plane. Cells in this type of confinement can migrate persistently, maintaining nearlyconstant cell shape, for time scales on the order of minutes to hours (Tsai et al., 2019; Garner et al., 2020). In order to capture leading edge dynamics on time scales more relevant to the stochastic growth of individual filaments, we performed highspeed (20 Hz) imaging of migrating HL60 cells. These experiments revealed dynamic, finescale fluctuations around the average leading edge shape (Figure 1a–c, Video 1, Figure 1—figure supplements 1–2, Materials and methods), where local instabilities in the leading edge emerge, grow, and then relax. Notably, these previouslyunobserved lamellipodial dynamics are phenotypically distinct from – and almost 100fold faster than – the oscillatory protrusionretraction cycles seen in other, slowermoving cell types (e.g. fibroblasts) (Giannone et al., 2004; Ryan et al., 2012; Ma et al., 2018).
We estimate we were able to reliably measure fluctuations with wavelengths as small as ~650 nm, and amplitudes down to ~65 nm, by fitting the phase contrast halo around the leading edge (Figure 1—figure supplements 1–2, Materials and methods). These values should approximately correspond to 25 actin filaments at physiological spacing (Svitkina et al., 1997) and 25 actin monomers assembled into a filament lattice along the direction of motion. While our measurements of shape dynamics cannot resolve polymerization events of individual filaments, our results are consistent with the hypothesis that stochasticity in actin growth at the level of monomer addition – occurring throughout the leading edge actin network – ultimately manifests as the observed micronscale leading edge fluctuations. In particular, kymograph analysis of curvature and velocity (Figure 1b–c) showed that relatively longlived shape fluctuations are formed by the continual timeintegration of seemingly uncorrelated and very shortlived (subsecond, submicron) velocity fluctuations. Because the average cell shape remains constant over time, there must be some form of feedback acting on leading edge curvature to sustain stable lamellipodial growth. These rich, measurable finescale dynamics therefore provide a unique opportunity to directly observe the timeevolution of leading edge maintenance. Taking advantage of our highprecision measurements, we aimed to quantitatively investigate the properties of the observed fluctuations, with the goal of determining the mechanisms by which molecular machinery at the leading edge coordinates the stochastic polymerization of individual actin filaments.
Lamellipodial stability mediated by viscous relaxation of shape fluctuations
The relaxation of finescale shape fluctuations back to the steadystate leading edge shape is essential for the long time scale stability of lamellipodia. As for any physical system, the nature of this relaxation reflects the system’s underlying physical properties; in this case, the characteristics of – and interactions between – actin filaments and the membrane. To provide a framework for exploration of the physical mechanisms underlying stable lamellipodial protrusion, we quantified the relaxation dynamics by performing timeautocorrelation analysis on the leading edge shape (Materials and methods). Applied in this context, this analytical technique calculates the extent to which the lamellipodium contour loses similarity with the shape at previous time points as fluctuations emerge and relax. As most material systems (activelydriven or otherwise) exhibit relaxation behavior with a characteristic wavelengthdependence (e.g. Figure 1d), we performed Fourier decomposition on the leading edge shape to separate out fluctuations at different length scales, and then performed autocorrelation analysis separately on each Fourier mode. We validated our analytical methods using simulations of membrane dynamics, for which there exists a wellestablished analytical theory (Materials and methods, Figure 1—figure supplement 3), and show that our results are not sensitive to an extension of our analysis to longer length and time scales (Materials and methods, Figure 1—figure supplements 4–5). Further, the membrane simulation control nicely demonstrates how visual features of curvature kymographs (e.g. Figure 1b) can be misleading (Materials and methods, Figure 1—figure supplement 3), and motivates the necessity of our more comprehensive technique.
Autocorrelation analysis revealed a monotonic relaxation of shape fluctuations at each wavelength (Figure 1e–f); the decay at every spatial scale is wellfit by an exponential form (Figure 1f), consistent with overdamped viscous relaxation. Importantly, we do not detect any increase in the autocorrelation over time, which would have appeared if there were any sustained, correlated growth of the fluctuations before they decay. This again suggests that the fluctuations arise from uncorrelated stochastic processes, such as fluctuations in actin density. A clear wavelengthdependence is observed, with shorter wavelengths decaying faster and having smaller amplitudes (Figure 1g–h). This general trend is shared by many physical systems with linear elastic constraints, such as idealized membranes (Brown, 2008) and polymers (De Gennes, 2002) freely fluctuating under Brownian motion, but can be contrasted with systems that have a dominant wavelength, as in the case of buckling or wrinkling of materials under compression (Cerda and Mahadevan, 2003). Importantly, these qualitative and quantitative properties of the leading edge fluctuations are not specific to cell type or experimental conditions (e.g. agarose overlay, ECM), as we also observe this phenomenon in fish epidermal keratocytes (Figure 1—figure supplement 6, Video 2).
Leading edge stability as an emergent property of branched actin growth
The rich behavior and quantitative nature of our leading edge shape fluctuation data made them ideal for comparison with physical models. In order to understand how molecularscale actin assembly and biomechanics might give rise to the observed micronscale shape dynamics, we aimed to reproduce this behavior in a stochastic model of branched actin network growth against a membrane (Figure 2a–c, Video 3, Materials and methods). Previous stochastic models of protrusive actinbased forces largely focused on actin polymerization against rigid obstacles (e.g. the bacterial cell wall for the Listeria comet tail Carlsson, 2001; Carlsson, 2003 or a single, flat membrane segment in models of lamellipodia Mueller et al., 2017). Expanding on this general framework, and in an approach conceptually similar to previous work simulating small (< 2 µm) patches of a lamellipodium (Schaus et al., 2007; Schaus and Borisy, 2008), we incorporated a twodimensional leading edge with filaments polymerizing against a flexible membrane, which we modeled as a system of flat membrane segments coupled elastically to each other. The size of the membrane segments was comparable to the spatial resolution of our experimental measurements, allowing us to assay fluctuations over a similar dynamic range of wavelengths. Simulated filaments apply force to the membrane following the classic untethered Brownian ratchet formalism (Mogilner and Oster, 1996), consistent with recent experiments showing that cellular protrusions are formed by largely untethered actin networks (Bisaria et al., 2020). Designed to be as comparable as possible to our experimental data, the model incorporated experimentally measured values from the literature for the membrane tension, membrane bending modulus, and biochemical rate constants (Tables 1–2, Mogilner and Oster, 1996; Lieber et al., 2013). As we were specifically interested in identifying biophysical mechanisms regulating leading edge stability, we minimized the model’s biological complexity by including only the core biochemical elements of actin network growth dynamics: polymerization, depolymerization, branching, and capping. All filament nucleation in the model occurs through dendritic branching observed in cells to be mediated by the Arp2/3 complex (Welch et al., 1997; Svitkina et al., 1997), which catalyzes the nucleation of new ‘daughter’ actin filaments as branches from the sides of preexisting ‘mother’ filaments at a characteristic angle of ~70° (Mullins et al., 1998; Volkmann et al., 2001; Rouiller et al., 2008). By simulating individual filament kinetics, the model captures the evolutionary dynamics of the filament network, allowing us to directly test hypothesized mechanisms for the interplay between actin network properties (e.g. filament orientation) and protrusion dynamics (Mogilner and Oster, 1996; Lacayo et al., 2007; Mueller et al., 2017; Grimm et al., 2003; Maly and Borisy, 2001; Schaus et al., 2007; Figure 2b).
To our great surprise, this very simple model was able to recapitulate stable leading edge fluctuations. Nascent leading edges reach steady state values for filament density, filament length, filament angle, membrane velocity, and (most importantly) membrane fluctuation amplitude within seconds, a biologically realistic time scale (Figure 2d–h). Furthermore, the steady state values obtained are in quantitative agreement with both our own experimental data and previously published measurements, with the model yielding mean values of: 0.3 filaments/nm for filament density (~30 nm filament spacing for a lamellipodium that is 10 filaments tall) (Abraham et al., 1999), ~150 nm for filament length, ~40° for filament angle (with respect to the direction of migration), ~0.35 nm/ms for membrane velocity, and ~50 nm for membrane fluctuation amplitude (Svitkina et al., 1997; Maly and Borisy, 2001; Verkhovsky et al., 2003).
As observed in the experimental measurements, simulated leading edge shape stability is mediated by an exponential decay of shape fluctuations (Figure 2i–j). Furthermore, the minimal model correctly predicts the monotonic trends of fluctuation amplitude and decay time scale with wavelength (Figure 2k–l) in a way that was not sensitive to our choices of simulation time step, membrane segment length, and overall length of the leading edge (Figure 2—figure supplements 1–3). It should be noted that the generation of the simulated data in Figure 2k–l did not involve any curvefitting (and therefore no free parameters that could be fit) to the experimentally measured autocorrelation dynamics in order to parameterize the model. Rather, the simulated fluctuation relaxation behavior, qualitatively reproducing our experimental measurements, emerges directly from the molecularscale actin growth model, in which all biochemical parameters were estimated from measurements in the existing literature (Tables 1–2) – leaving no free simulation parameters.
Predicting effects of drug treatment with Latrunculin B
We were interested in further assaying the predictive power of this minimal stochastic model by determining whether the output of the simulations was congruent with experimental observations under conditions that had not been tested prior to model development. As an example, we elected to test whether the model could correctly predict the response of HL60 cells to treatment with the drug Latrunculin B, which binds to and sequesters actin monomers. Qualitatively, cells treated with Latrunculin B (Video 4) present with enhanced bleb formation and more variable leading edge shapes, in comparison with cells treated with a DMSO vehicle control (Video 5). In our model, addition of this drug can be simulated by reducing the free monomer concentration, which consequently reduces both the polymerization rate and the branching rate. At low effective doses, subtle but measurable changes to leading edge fluctuations were predicted: specifically, an increase in the amplitudes at large wavelengths, and a decrease in the decay rates across all wavelengths (Figure 3a–b). Our experimental results were consistent with these quantitative predictions; Latrunculin Btreated cells exhibited increased fluctuation amplitudes and decreased fluctuation rates over the predicted ranges (Figure 3c–d).
Geometry as the core determinant of simulated leading edge stability
Given the success of the model in reproducing experimental results, we next wanted to determine which features of the simulation were responsible for leading edge stability and relaxation of fluctuations. The simplicity of the model allowed us to determine the stability mechanism by process of elimination, selectively removing elements of the model (Materials and methods) and determining whether stability was retained. To assay the importance of membrane tension and bending rigidity, which has been suggested to be a key factor regulating lamellipodial organization (Batchelder et al., 2011; Sens and Plastino, 2015), we simply removed the forces between the membrane segments (Figure 2a, F_{spring}) from the simulation (Materials and methods). Surprisingly, the coupling between the membrane segments (i.e. the effects of tension and bending at length scales larger than the size of an individual membrane segment) was completely dispensable for leading edge stability (Figure 4a–e).
Following a similar process of elimination, we determined that in fact only two elements were required for stability. First, as reported previously for dendritic actin network polymerization against a single stiff obstacle, it was necessary to constrain branching to occur only within a fixed distance from the leading edge membrane (Figure 2a–b) in order to maintain a steady state actin density (Carlsson, 2001). The molecular motivation for this spatially limited ‘branching window’ is rooted in that fact that activators of the Arp2/3 complex, which render Arp2/3 competent for actin filament nucleation, are typically membraneassociated proteins (Suetsugu, 2013). Second, we found that stability is inherently tied to the ability of filaments to spread laterally to neighboring membrane segments (Figure 2b IIIII, Figure 4f–i). Recall that, because the branched actin network geometry causes filaments to grow, on average, at an angle relative to the membrane normal (Maly and Borisy, 2001; Verkhovsky et al., 2003), polymerizing tips spread laterally along the leading edge (Lacayo et al., 2007; Grimm et al., 2003). Removing filament spreading from the model by fixing filaments to remain associated with their nearest membrane segment at birth (Materials and methods) led to actin density divergence: network regions with low filament density eventually underwent complete depolymerization, while highdensity regions continued to accumulate actin (Figure 4f–g).
These findings lead us to a simple molecular feedback mechanism for leading edge stability, based on a synergy between filament spreading and membraneproximal branching (Figure 4j): To begin with, regions with initially high filament density come to protrude beyond the average position of the rest of the membrane, representing the emergence of a leading edge shape fluctuation (Figure 4j, I, II). This induces asymmetric filament spreading, where filaments from highdensity regions can spread productively into neighboring regions (Figure 2b, II), but filaments spreading from adjacent lowdensity regions cannot keep up with the fast moving membrane segments (Figure 2b, III), and thus are unproductive (Figure 4j, II, III). In this way, the branched geometry inherent to dendritic actin polymerization, as well as its interaction with the shape of the membrane, naturally encodes leading edge stability (Figure 4j, IV). Thus our results directly demonstrate a ‘stabilitythroughspreading’ mechanism that has previously only been assumed in meanfield analytical theories (Lacayo et al., 2007; Grimm et al., 2003). Remarkably, this means that leading edge maintenance is an intrinsic, emergent property of branched actin network growth against a membrane, without requiring any further regulatory governance. Geometrical constraints imposed simply by the nature of membraneproximal actin branching ensure that any small variations in either local actin filament density or growth rate are inherently selfcorrecting to regress toward the mean.
Of note, it has previously been shown that Arp2/3mediated branching is required for lamellipodial formation in a wide variety of cell types; cells with inhibited or depleted Arp2/3 complex exhibit complete disruption of the lamellipodium shape and often switch to a different mode of migration altogether, such as filopodial motility (FritzLaylin et al., 2017b; Henson et al., 2015; Wu et al., 2012; Davidson et al., 2018). Indeed, HL60s treated with the Arp2/3 inhibitor CK666 have extremely variable leading edge shapes, characterized by long, thin filopodialike protrusions (Video 6). Our theoretical results provide a mechanistic interpretation for this striking phenomenon, suggesting that the vital lamellipodial maintenance role of Arp2/3mediated branching stems from its ability to mediate efficient filament spreading and equilibration of actin density fluctuations, purely because the daughter filament always grows at an angle distinct from its mother.
Optimal suppression of fluctuations by the highly conserved ~70° branching angle
Given the essential contribution of branched network geometry to the stability of the simulated leading edges, we reasoned that variations in the branching geometry alone might have a significant effect on leading edge fluctuations. We therefore performed simulations to determine the effects of changing the average branching angle and branching angle variability on filament orientation, filament density, and leading edge fluctuation fit parameters (Figure 5). In this context, we highlight the distinction between the branching angle, θ_{br} (i.e. the angle of a daughter filament relative to its mother), and the filament angle or orientation, θ_{f} (i.e. the angle of a filament relative to the direction of migration) (Figure 5a, inset). Due to the sterotypical branching angle, θ_{br}, there is a direct correspondence between the orientation, θ_{f}^{mother}, of a mother filament and the orientation, θ_{f}^{daughter}, she passes on to her daughter branches. Our simulations are thus, in effect, selection assays, as mother filaments compete to stay within the fixed branching window, spawn daughter branches, and thus pass down their angle to their progeny (Maly and Borisy, 2001; Schaus et al., 2007). For example, when filaments are initialized with a random orientation, and the branching angle is fixed (i.e. there is no variability in the branching angle), only a handful of the initial filament angles (θ_{f}) survive until the end of the simulation (Figure 5a). The surviving, successful filament angles are narrowly and symmetrically distributed around one half of the branching angle (Figure 5a–c). This optimal filament angle allows mother and daughter filaments to branch back and forth symmetrically about the membrane normal, such that mother filaments do not outcompete their progeny (as has been described previously) (Maly and Borisy, 2001; Schaus et al., 2007; Schaus and Borisy, 2008; Figure 2b).
In living cells, branching is mediated by the Arp2/3 complex, which has been experimentally measured to form highly regular and sterotyped branches at ~70° (Mullins et al., 1998; Volkmann et al., 2001; Rouiller et al., 2008). Intriguingly, this protein complex is highly conserved (Welch et al., 1997), with measurements of the branching angle in a wide variety of species, including protists (Mullins et al., 1998; Volkmann et al., 2001; Rouiller et al., 2008; Blanchoin et al., 2000), yeast (Rouiller et al., 2008), mammals (Rouiller et al., 2008; Blanchoin et al., 2000; Cai et al., 2008), and amphibians (Svitkina and Borisy, 1999), using various experimental techniques (platinum replica electron microscopy, cryoelectron microscopy, and total internal reflection microscopy), all falling within the range of 67–78° (± 2–13°). The high degree of conservation hints that this specific angle might carry some functional optimality, but the question has not been addressable experimentally; due to the lack of naturally occuring Arp2/3 variants with a substantially different branching angle, an alternative branching structure would hypothetically have to be designed do novo, presumably by altering the protein interaction interface by which the Arp2/3 complex binds to the side of a mother filament (Volkmann et al., 2001). We thus sought to explore the possible functional significance of this conserved angle using our minimal stochastic model. Excitingly, we found that in simulations with no branch angle variability, a 70–80° branching angle was optimal for minimizing both actin density fluctuations (Figure 5d and i) and leading edge fluctuation amplitudes for wavelengths smaller than ~2 µm (Figure 5f and h). These smoothing effects are therefore predicted to be relevant within the experimentallymeasurable range of wavelengths (between ~0.7 µm and ~2 µm), but are most beneficial for the smallest wavelengths resolved by our simulations (down to ~0.3 µm) – closest to the length scales of individual filament polymerization. Overall, these results provide tantalizing mechanistic insight into the longstanding question of why the characteristic branching angle is so ubiquitous.
Heritability of filament orientation (i.e., the extent to which mother filament orientation determines the orientation of the daughters) is set by the degree of variability in the branching angle (which, in turn, reflects the influence of thermal fluctuations). Perhaps unsurprisingly, decreasing this orientational heritability significantly reduces the dependence of fluctuation amplitude and filament density variability on the branching angle (Figure 5j–q), and thereby counteracts the beneficial effect of the optimal angle on leading edge fluctuations. Introducing a branch angle variability of ±2° (on the lower end of the experimentallymeasured values) broadens the range of nearoptimal branch angles but maintains the optimum at ~70–80° (Figure 5j and l), while introducing a variability of ±10° (on the higher end of the measured range) completely removes the optimum (Figure 5n and p). In both cases, increasing the branch angle variability increases the minimum possible fluctuation amplitude (Figure 5h, l and p – insets) and filament density variability (Figure 5i, m and q), representing a decrease in the noisesuppression capabilities of the system. Overall, these results provide strong support for the idea that actin network geometry is not only essential for leading edge stability, but also plays a major role in determining the fundamental limits of smoothness in lamellipodial shape.
Discussion
The emergence of robust collective behaviors from stochastic elements is an enduring biological mystery which we are only beginning to unravel (Huang et al., 2016; Battich et al., 2015; Chang and Marshall, 2017; Mohapatra et al., 2016; Raj and van Oudenaarden, 2008; Gray et al., 2019; Oates, 2011). The apparent dichotomy in actinbased motility between the random elongation of individual filaments and the stable formation of smooth and persistent higher order actin structures such as lamellipodia exemplifies this enigma, and provides an avenue toward understanding general strategies for noise suppression in biological systems.
In recent years it has become clear that perturbationfree experiments which examine fluctuations around the mean at steady state (in contrast to probing the change in the mean due to a perturbation) can be a powerful tool for understanding noisy systems (Welf and Danuser, 2014). In this work, application of that principle in combination with highprecision measurements, quantitative analytical techniques, and physical modeling led to the surprising revelation that the suppression of stochastic fluctuations naturally emerges from the interactions between a growing actin network and the leading edge membrane, with no additional feedback required. Our insights into the molecular mechanisms mediating lamellipodial stability were largely enabled by experimentally measuring micronscale leading edge shape dynamics and comparing them to a molecularscale actin network growth model that correctly predicts this emergent behavior (as well as many other experimentallymeasured features of lamellipodial actin networks). Ultimately, we hope our results inspire future experimental work to directly measure the nanometerscale interactions predicted by our simulations, which may be accomplished using superresolution imaging of actin dynamics in vivo or in an in vitro reconstitution of lamellipodial protrusion (i.e. branched polymerization of actin driving the motion of a flexible barrier).
The model developed in this work provides an understanding of the basic biophysical mechanisms underlying lamellipodial migration. Of course, living cells are home to array of additional complexities which are likely to further modulate leading edge fluctuations and stability, and our model may provide a framework for future exploration of such effects across diverse cell types and experimental conditions. Cells migrating in vivo inevitably experience a much more challenging and dynamic environment, in which additional feedback mechanisms will almost certainly be required for the maintenance of polarized migration. The simplicity and biophysical realism of our modeling framework should make it particularly wellsuited for future studies focused on predicting and understanding the effects of additional potential feedback mechanisms, including tethering (Mogilner and Oster, 2003; Soo and Theriot, 2005; Alberts and Odell, 2004; Kuo and McGrath, 2000), a limiting monomer concentration (Mullins et al., 2018), forcedependent branching (Risca et al., 2012; Parekh et al., 2005; Chaudhuri et al., 2007), and regulation by curvaturesensing proteins (Tsujita et al., 2015; Zhao et al., 2011). Further, this computational model might be useful for exploring the effects of various extracellular forces, such as those produced by obstacles or variations in matrix density, or intracellular forces, such as those produced by hydrostatics. We also note that a certain degree of biochemical signaling is implicit in our model in the form of biochemical rate constants that are invariant in time and space; this property relies on signaling networks to maintain uniform gradients of actinassociated molecules (Devreotes et al., 2017). How local biochemical control (or lack thereof) over these rate constants might affect leading edge fluctuations remains an interesting avenue for future investigation, both theoretically and experimentally.
The defining characteristic and major advance of our model was the explicit inclusion of both the evolutionary dynamics of the actin network and its interaction with the twodimensional geometry of the leading edge. By selectively removing elements of the model, we determined that lateral filament spreading, combined with a fixed branching window, is indispensable for leading edge stability. This highlights the crucial role in lamellipodial maintenance of the branched structure of actin networks, wherein each daughter filament inherits angular information from its mother. Our further investigations into the evolutionary properties of actin network growth revealed that a ~70–80° branch angle maximally suppresses finescale actin density and leading edge shape fluctuations, showing for the first time that Arp2/3mediated branching imparts optimal functionality, as was long hypothesized based on strong sequence, structural, and functional conservation throughout the eukaryotic tree of life. It is interesting to note that the evolutionarilyconserved branching angle that we find maximally suppresses leading edge fluctuations appears not to be the same angle that optimizes the polymerization velocity of single filaments – predicted to be a broad angle closer to ~90–100° (and quite loaddependent) in the lowload regime (Mogilner and Oster, 1996). This contrast suggests that evolutionary selection acts at the level of actin network properties, rather than force production by individual filaments.
Returning to the broader question of how noisy biological systems control for stochasticity, we find that stability in the case of lamellipodial dynamics is inherently encoded by the geometry of branched actin network growth. It will be interesting to see whether similar principles hold for other cytoskeletal structures with clear geometric constraints, such as endocytic pits, the cytokinetic ring, and the mitotic spindle.
Materials and methods
HL60 cell culture and differentiation
Request a detailed protocolHL60 cells were cultured as described previously (Millius and Weiner, 2009; Garner et al., 2020). In brief, cells were maintained at a density of 0.1–1 x 10^{5} cells/mL by passaging every 2–3 days into fresh RPMI media supplemented with 10% heatinactivated fetal bovine serum and antibiotics/antimycotics. Supplementation with 1.57% DMSO was used to differentiate the cells into a neutrophillike state. Cells were subsequently extracted for experiments at 6 days postdifferentiation. Our HL60 cell line was obtained from Orion Weiner’s lab at UCSF, who originally received them from Henry Bourne’s lab at UCSF. The identity of this suspension cell line was confirmed based on the behavior of the cells, including differentiating into a neutrophillike state upon exposure to DMSO that exhibits characteristic phenotypes for substrate adhesion, rapid migration, and elongated morphology. HL60s are not listed as a misidentified cell line on the Register of Misidentified Cell Lines. The cell line tested negative to mycoplasma contamination.
Underagarose motility assays with HL60s
Request a detailed protocolDifferentiated HL60 cells were plated on fibronectincoated coverslips and then overlaid with a 1% agarose pad containing 1 nM fMLP (to enhance migratory behavior), as described previously (Garner et al., 2020). Microscopy of the migrating cells was performed at 37 °C, using transmitted light to image phase contrast on an epifluorescence microscope at ×100 magnification (100 × 1.45 NA Plan Apo oil objective, Nikon MRD31905). A more detailed description of our microscopy system can be found in previous publications (Garner et al., 2020). For treatment with Latrunculin B or CK666, the drug was embedded into the agarose pad by adding the drug to the unpolymerized agarose pad solution before gelling (at the same time as adding fMLP), such that drug treatment begins when cells are overlaid with the agarose pad and is maintained throughout imaging. Cells were imaged at 45 min postplating. Drugs were first diluted down to 1000 X in DMSO, then added to the agarose solution at a dilution of 1:1,000 (for a final concentration of 30 nM for Latrunculin B and 100 μM for CK666), giving a final DMSO concentration in the pad of 0.1%. Controls were performed by adding 0.1% DMSO to the agarose pad alone.
Keratocyte isolation and motility assays
Request a detailed protocolKeratocytes were cultured from wildtype zebrafish embryos at 2 days postfertilization as described previously (Lou et al., 2015). Briefly, zebrafish embryos were collected at 2 days postfertilization, dechorionated, and anesthetized using tricaine. To dissociate the keratocytes, dechorionated fish were then washed in PBS, incubated in Cell Dissociation Buffer for 30 min at 4 °C, incubated in a solution of 0.25% trypsin and 1 mM EDTA for ~15 min at 28 °C, and then incubated in fetal bovine serum to quench the trypsin. From this point on, cells were maintained in antibiotic and antimycotic to deter microbial growth. The keratocyterich supernatant was then concentrated by centrifugation at 500 g for 3 min. Keratocytes were then plated on collagencoated coverslips and incubated at room temperature for ~1 hr to allow cells to adhere. Once adherent, the supernatant was exchanged for imaging media (10% fetal bovine serum in L15) and allowed to incubate another 15 min at room temperature before imaging. The keratocytes were imaged at 28 °C under similar conditions as those used for HL60 cells. Data included in Figure 1—figure supplement 6 represent cells from a single coverslip. Experiments were approved by University of Washington Institutional Animal Care and Use Committee (protocol 4427–01).
Image segmentation
Request a detailed protocolMost segmentation algorithms penalize curvature in the contour in order to reduce the noise in the fitting algorithm (Seroussi et al., 2012). However, this runs the risk of introducing artificial correlations and structure into the data. For example, a springlike curvature penalty would artificially make fluctuations appear to be stretchdominated. Therefore, we performed the following custom segmentation algorithm to avoid these potential artifacts (Figure 1—figure supplement 1). Phase contrast timelapse videos were manually aligned to the direction of motion of the cell, such that each cell migrates up the yaxis on an xy coordinate system (Figure 1—figure supplement 1a, b). The videos were then cropped to isolate the cell leading edges and exclude the cell body, for easier segmentation. If the cell migrates up the yaxis, this means every image pixel along the xaxis has an associated leading edge position along the yaxis. A leading edge yposition was assigned to each xaxis pixel independently of knowledge about neighboring xaxis pixels, to avoid injecting the artifacts discussed above. A manual segmentation was performed for the first time point in each movie. A custom, automated segmentation algorithm written in MATLAB then performed a line scan of the phase contrast intensity along the yaxis separately for each pixel along the xaxis. For each line scan, the algorithm performs a local search for the leading edge position, constrained to be within a fixed number of pixels from either the manual segmentation (for the first time point) or the previous time point (for subsequent timepoints). The leading edge position is defined as the midpoint between the brightest phase intensity (phase halo) and the point of steepest intensity gradient (transition from phase halo to phasedense cytoplasm).
Preparation of the curvature and velocity kymographs
Request a detailed protocolThe curvature and velocity kymographs (Figure 1b–c, Figure 1—figure supplement 3) were prepared using custom MATLAB code. Curvatures are calculated as the inverse radius of the bestfit circle corresponding to a 30 pixelwide (~1.5 µm) region about each position. The most prominent fluctuation events seen in the curvature kymographs (Figure 1b, Figure 1—figure supplement 3) somewhat correspond to (but do not exactly match) the length of the fitting window. For example, the simulated and experimental data shown in Figure 1—figure supplement 3 were fit using the same ~1.5 µm fitting window and have similar apparent ‘dominant wavemodes’ of ~3 μm, or twice the fitting window. However, despite being fit with the same fitting window, the simulated data has an observably smaller apparent ‘dominant wavemode’ than the experimental data. Velocities were calculated as the distance traveled over 250ms (five 50ms timepoints) nonoverlapping windows.
Processing of segmented cell shapes for autocorrelation analysis
Request a detailed protocolKymographs of curvature and velocity such as those shown in Figure 1b–c, while helpful to obtain a qualitative sense of the fluctuation data, are visually dominated by the largest sizescale features of the leading edge. They thus offer an incomplete description of the shape fluctuations (Ma et al., 2018) – notably deemphasizing the finescale features that are the subject of this study. Further, curvature kymographs emphasize features that are approximately the same size as the fitting window, and fail to pick up fluctuations at different size scales. To perform a quantitative analysis which faithfully captures fluctuations at all size scales, we choose to perform Fourier decomposition on the leading edge shape, and analyze the dynamics of each wavemode separately. As cells migrate, their global leading edge shape undergoes long timescale changes, such as variations in width, largescale curvature, or slight turning of the cell, which can dominate the Fourier amplitudes and the subsequent autocorrelation signal. As we are most interested in extracting the finescale fluctuations, we performed background subtraction on the segmented leading edge shapes (Figure 1—figure supplement 2). To do this, we defined the ‘background’ leading edge shape as the contour after smoothing (by the lowess method, using a span of 7 μm). This rather large smoothing window was chosen specifically to preserve finescale features. The backgroundsubtracted yposition is thus defined as the difference between the segmented leading edge and its smoothed counterpart. This process removed the largescale features of the leading edge. We next wanted to remove the longtimescale features of the leading edge, so we also subtracted the timeaveraged backgroundsubtracted yposition for each xpixel. Altogether, these preprocessing steps maintained the features of interest in the curvature kymograph (Figure 1—figure supplement 3). After performing background subtraction, we still needed to control for changes in leading edge width over time. The wavelengths represented in the Fourier transform are defined as λ = L/n, where L is the length of the leading edge, and n is an integer from 0 to one half the number of pixels. If the leading edge length were to vary over time, then so would the wavelengths, making it impossible to track the behavior of a single wavelength fluctuation over time. We thus cropped the dataset along the xaxis to include only pixels which contain the cell for all timepoints in the video, thereby extracting a fixedlength leading edge subset for further analyses.
Autocorrelation analysis and fitting
Request a detailed protocolAutocorrelation analysis and fitting were performed separately for each cell and simulation. For experimental data, the entire video was analyzed. For simulated data, analysis was only performed on the time points after the simulation had reached steady state, for which we used a conservative cutoff of 10 s (see Figure 2d–h). To separate out fluctuations at different length scales, we first performed a spatial Fourier transform on the leading edge shape. Referencing the coordinate system defined in Figure 2a, the pixels (experiments) or membrane segments (simulations) are equally spaced in the xdirection, allowing us to perform a onedimensional Fourier transform (MATLAB fft() function, which assumes periodic boundary conditions) on the ypositions of a segmented leading edge for each time point. We then normalize the Fourier transform by a factor of dx/$\surd L$, where dx is the pixel/membrane segment size and L is leading edge length. This normalization preserves the variance and accounts for the pixel size. To measure the fluctuation relaxation, we calculated the timeaveraged autocorrelation (A_{n}(τ) = < Y_{n}(t+τ)∙Y_{n}*(t) > _{t}, using nonoverlapping windows in t) of each Fourier mode amplitude. The autocorrelation function extracted from this analysis contains complex elements of the form A = a + ib. We performed all plotting and fitting on the complex magnitude (sqrt (a^{2} + b^{2})) of the autocorrelation function, which is most representative of the total autocorrelation.
Note that because we are plotting the complex magnitude (which is always positive), the autocorrelation plots shown in Figures 1e–f ,–2i–j are expected to decay to some nonzero background noise window, rather than to zero. Indeed, the membrane simulation control (see Validation of autocorrelation analysis implementation), which is predicted to have a purely exponentiallydecaying autocorrelation function, also shows a decay to a noise window at long times (Figure 1—figure supplement 4). We fit each Fourier mode timeautocorrelation to the exponential decay function described in the main text (Figure 1f), fitting ln(A_{n}(τ)) vs τ to a line using MATLAB’s polyfit function, to extract fit parameters for each cell and simulation. Each curve was fit out to a drop in the amplitude by a factor of e/2, or at least 10 points. To average fit parameters over many cells, we controlled for celltocell variability in leading edge length by binning the parameters by spatial frequency, and then calculating summary statistics separately for each bin.
The spatial background subtraction performed on the leading edge shape (discussed in Materials and methods: Processing of segmented cell shapes for autocorrelation analysis) was necessary to extract the finescale shape fluctuations studied in this work. This background subtraction is expected to remove fluctuations with wavelengths larger than ~7 µm (i.e. reduce their amplitude to zero). For this reason, only wavelengths less than 7 µm are plotted in Figure 1e–h. We note that it is possible that fluctuations with wavelengths less than, but near 7 µm might also have slightly reduced measured amplitudes (i.e. the shape of the curve in Figure 1g may artificially level off at low spatial frequency). However, any such effect would be performed uniformly in time, and therefore is not expected to affect the measured temporal dynamics (Figure 1h). Indeed, when we extend the span of our background subtraction by ~50% (up to 10 μm), we find that only the amplitude of the largest mode is altered (slightly increased) and the measured relaxation timescales are not affected (Figure 1—figure supplement 5).
Validation of autocorrelation analysis implementation
Request a detailed protocolTo validate our autocorrelation method, we analyzed control simulations of a membrane freely fluctuating under Brownian motion in the absence of actin, and showed it recapitulates predictions from analytical theory for this system (Figure 1—figure supplement 3). These simulations were performed exactly as in the leading edge simulations, using the same parameters, but without actin. The equation of motion used for membrane segments in the control simulations, as well as a derivation of the associated autocorrelation function, can be found in Appendix 2.4.4: Choice of timestep. Interestingly, the curvature kymographs of these control simulations exhibit striking visual features reminiscent of instabilities, dominant wavemodes, or oscillations – and yet such effects are absent from this system by definition (which we confirm quantitatively using our autocorrelation method, Figure 1—figure supplement 3). This suggests that similar features in the experimentallymeasured curvature kymograph are also not indicative of instabilities, dominant wavemodes, or oscillations, which we confirmed by autocorrelation analysis.
Modeling
Request a detailed protocolHere, we briefly describe the geometry and major elements of the model. Please see the Appendix for a detailed description of the model and Tables 1 and 2 for a list of the chosen parameters. A simulated patch of leading edge was modeled by a branched network of actin filaments stochastically polymerizing towards a 2D strip of membrane, subject to periodic boundary conditions. The 2D strip was discretized as membrane segments that are fixed in position along one axis and move only along the direction of motion of the simulated cell. Stochastic, fixed time step Brownian dynamics simulations, implemented with custom MATLAB code, were performed to update the membrane position and actin network properties. Actin network growth evolved from constant rate Poisson processes for polymerization, depolymerization, branching, and capping. Once polymerized, filaments were fixed in position at their branch point of origin (in the lab frame of reference), and did not undergo retrograde flow (i.e. translation of the filament position opposite the direction of migration) or translational diffusion. The membrane strip was subject to forces of bending and stretching, drag from fluid viscosity, as well as the force of actin (Mogilner and Oster, 1996; Peskin et al., 1993; Carlsson, 2001). Filaments apply force to the membrane segments according to the untethered Brownian ratchet formalism (Mogilner and Oster, 1996), in which filament pointed end positions are assumed to be rigidly connected to the network (via their branch point of origin) and their barbed end positions are able to freely fluctuate. As previously, we ignore the possibility of filament buckling due to the fact that lamellipodial filaments exist in a sufficiently lowload, high branch density regime (Mogilner and Oster, 1996), and experimental evidence shows no indication of buckling (Svitkina et al., 1997). We expanded this formalism (which previously only considered filament fluctuations perpendicular to the filament’s long axis) to include all fluctuations of the filament along the filament’s short and long axes. Each filament pushes the membrane segment that spans the growing tip’s xposition. We note that the filament angles used to determine the filament forces on the membrane and presented throughout this work are always calculated relative to the global average direction of motion of the leading edge, rather than the local average membrane normal (a simplifying approximation necessitated by the discrete geometry and motivated by the shallow curvatures exhibited by cell leading edges). Control simulations were run to verify that the leading edge fluctuation behavior described in this work was not dependent on the temporal discretization (i.e. simulations were run with sufficiently small timesteps to resolve the fastest dynamics, Figure 2—figure supplement 1), spatial discretization (i.e. simulations were run with sufficiently short membrane segments to resolve the smallest length scales at which there is significant bending, Figure 2—figure supplement 2), or leading edge length (i.e. the periodic boundary conditions were implemented correctly, such that a simulated small patch of leading edge behaves identically to an equivalently sized portion of a larger simulated patch of leading edge, Figure 2—figure supplement 3). In cases where membrane tension and bending rigidity were removed, these forces were simply not calculated in the simulation (Figure 4a–e). To remove filament spreading, we modified how the filament position was updated upon addition of a monomer in order to maintain the growing filament tip’s xposition. Addition of monomers contributed only to changes in the barbedend yposition, leaving the xposition intact, while updating the filament length correctly (effectively sliding the pointed end xposition backwards, rather than advancing the barbed end xposition forwards, Figure 4f–i).
Appendix 1
Detailed description of the model
In this section, we outline the main features of the model and reference the Detailed Derivations for more detailed derivations.
Model geometry
Given the flat structure of lamellipodial protrusions, we considered leading edge dynamics in two dimensions. In cartesian coordinates, the cell migrates in the xy plane along the yaxis and maintains a fixed leading edge height of 200 nm along the zaxis. The leading edge membrane was modeled as a 2D strip that restricts bending and stretching in the xy plane and is perfectly flat along the zaxis. We discretized the membrane such that a 20 μm leading edge membrane was modeled as 200 flat, rectangular segments (each 100 nm in length along the xaxis and 200 nm in height along zaxis) whose surface normal is fixed to lie along the yaxis. The membrane is implented using a Monge parameterization, such that these membrane segments are fixed in position along the x and zaxes and move only along the yaxis (the direction of motion of the simulated cell).
Updating the membrane position
The membrane segments are assumed to move in a viscous medium at low Reynolds number, with drag force ${F}_{drag}=\gamma \frac{\partial y}{\partial t}$ and Stoke’s drag coefficient $\gamma =6\pi \eta r$, where r is the membrane segment length (xaxis) and $\eta $ is the dynamic viscosity of water. The membrane also acts under the forces of membrane bending and stretching (${F}_{S/B}$, see Appendix section 1.3), and Brownian ratchet forces by the actin filaments (${F}_{BR}$, see Appendix section 1.4). This gives us the following equation of motion.
Brownian dynamics simulations, implemented with a 4th order RungeKutta algorithm, were performed to update the membrane segment positions. Thermal fluctuations of the membrane were ignored in this implementation, as thermal fluctuations of the (much stiffer) actin filaments dominate the membrane’s motion as well as monomer incorporation into the actin network (Mogilner and Oster, 1996). See Appendix Detailed Derivations IV for a discussion of the numerical approximations included in the simulations.
Forces of membrane bending and stretching
The simulated leading edge membrane acts under the energetic constrains of stretching and bending, characterized by the experimentally measurable parameters of membrane tension ($\sigma$, $pN\cdot n{m}^{1}$) and bending modulus ($\kappa$, $pN\cdot nm$), and using the following energy functional:
where $h$ is the height of the leading edge in the zdimension,$\frac{\partial y}{\partial x}$ is extension of the membrane (i.e. an increase in contour length), and $\frac{{\partial}^{2}y}{\partial {x}^{2}}$ is curvature in (bending of) the membrane. The force of membrane stretching and bending on a single membrane segment (${F}_{S/B}$) is defined as the free energy gained by movement of the segment ($seg$), and takes the following form… (See Detailed Derivations I for a derivation of the functional derivative and resulting force.)
Actin filament Brownian ratchet forces
Actin filaments constantly undulate due to thermal fluctuations, bending and stretching to sample their conformational space. The presence of the membrane restricts fluctuations of the filament past the membrane, presenting an entropic cost and a reduction in the free energy of the filament. The force of the filament exerted on the membrane segment, ${F}_{actin}$, can therefore be calculated as the gain in free energy, G, by an infinitesmal movement of the membrane position, $\xi $:
The partition function, ${Z}_{conf}$, determined by the energetic cost of bending the filament, ${E}_{bend}$, defines the conformational landscape. A few considerations must be made to determine the partition funciton in our system. Each filament applies force only to the membrane segment under which the filament’s barbed (growing) end equilibrium position sits, making the approximation that each membrane segment acts as an infinite wall past which filament fluctuations are blocked. Given this assumption, the membrane only restricts filament fluctuations along the yaxis. The partition function is therefore integrated over all xpositions, but only the subset of ypositions where the filament is not restricted by the membrane. With these considerations in mind, we arrive at our partition function: (See Appendix Detailed Derivations II for a derivation of the filament bending energy and Detailed Derivations III for the derivation of the force).
where y is measured relative to the membrane surface and $\beta =\frac{1}{{k}_{B}T}$. The final equation for the force of a filament on the membrane becomes…
where y_{0} is the equilibrium filament position measured relative to the membrane surface (the filament pokes through the membrane for ${y}_{0}<0$) and
and ${\overline{\kappa}}_{\u27c2}=\frac{3}{{L}^{2}N}$, ${\overline{\kappa}}_{\parallel}=\frac{16}{{L}^{2}{N}^{2}}$, and $N=\frac{L}{{l}_{p}}$ for a filament with length L and persistence length l_{p} at an angle $\theta $ relative to the membrane segment normal.
Actin network dynamics
Actin network dynamics including polymerization, depolymerization, branching, and capping were assumed to be independent, constant rate Poisson processes. The choice of rates are described in their respective sections below. For a given time step of $\mathrm{\Delta}t$ and rate $r$, the probability of an event happening during any given time step is…
where ${e}^{r\cdot \mathrm{\Delta}t}$ is the probability that the event did NOT take place within a time step of $\mathrm{\Delta}t$. A random number generator was used to determine which processes occurred in each time step. In particular, a random number was chosen between 0 and 1. If the random number lied below the probability $p$, then the event occured. For each simulation the random number generator was seeded with a unique, semirandom number, based on the current time.
Polymerization and depolymerization
The rate of polymerization was assumed to be $r}_{on}={k}_{on}\cdot M{\cdot}_{ga{p}^{\prime}$, where $k}_{on$ is the rate of polymerization per free monomer concentration in $\mathrm{m}\mathrm{o}\mathrm{n}\mathrm{o}\mathrm{m}\mathrm{e}\mathrm{r}\mathrm{s}\phantom{\rule{thinmathspace}{0ex}}m{s}^{}1\phantom{\rule{thinmathspace}{0ex}}\mu {M}^{1}$, M is the free monomer concentration in $\mu M$, and $p}_{gap$ is the probability that enough space opens up in between the filament tip and the membrane to add a monomer (such that the polymerization rate far away from the membrane is ${r}_{(y=\mathrm{\infty})}={k}_{on}\cdot M$). Previously, it was determined that this probability is set by thermal fluctuations of the filament (Mogilner and Oster, 1996). In other words, the probability of adding a monomer is determined by the probability that thermal fluctuations, by chance, overcome the bending and stretching energies of the filament and bend the filament tip away from the membrane enough to open up a space of sufficient size to add a monomer. Given a monomer width $\mathrm{\Delta}$, the probability of adding a monomer is…
where ${\kappa}_{\text{Eff}}$ takes into account the thermal energy (${k}_{B}T$) as well as the flexibility, length, and orientation of the filament. Depolymerization was assumed to be a constant rate process with rate ${r}_{off}$ ($\text{monomers}m{s}^{1}$), independent of membrane proximity.
Branching
The rate of branching was calculated for each filament at each time step, such that ${r}_{\text{branch}}={k}_{\text{branch}}\cdot M\cdot {l}_{\text{branch}}$. Here ${k}_{\text{branch}}$ is the branching rate per free monomer concentration per length of mother filament in units of $\text{branches}m{s}^{1}\mu {M}^{1}n{m}^{1}$, M is the free monomer concentration in $\mu M$, and ${l}_{\text{branch}}$ is the length of the filament that sits inside the branching window. A fixedlength branching window is required for model stability, and is wellsupported by experimental evidence that branching activation is localized to the leading edge membrane. Unlike linear polymerization, branching was not inhibited by membrane proximity. For simplicity, new branches were placed on the tip of the mother filament. If this placement caused the new branch tip position to extend past the membrane, then the branches were placed on the side of the mother filament such that the branch tip is flush with the membrane. The angle of the branch relative to the mother filament was randomly selected from a normal distribution with mean $\mu ={\theta}_{branch}$ and standard deviation $\sigma =\mathrm{\Delta}{\theta}_{branch}$ as specified for each simulation. The side of the mother filament on which the branch was placed was random.
Capping
Capping was assumed to be a constant rate process with rate r_{c} ($m{s}^{1}$), independent of membrane proximity. Filaments were not allowed to uncap. Capped filaments were not allowed to polymerize or depolymerize.
Filament deletion
Our simulations were intended to capture only leading edge actin dynamics. We therefore chose, for simulation efficiency, to only keep track of filaments actively applying force to the membrane. Filaments which were both (1) capped and (2) cumulatively provided less than 0.1% of the force on a given membrane segment were deleted from the simulation.
Detailed Derivations
This Detailed Derivations contains detailed derivations and clarifications for the material discussed in the previous section of the Appendix: “Detailed description of the model”.
Detailed Derivations I: Forces of membrane bending and stretching
In this section, we use the energy functional for membrane elasticity defined in section A.3 to calculate the elastic forces on a discrete membrane segment, filling in the steps of Appendix equation (3).
Functional derivative of the membrane stretch/bend energy functional
We can solve the functional derivative $\frac{\delta E}{\delta y}$ for our particular energy function (Appendix equation (2)) using the fact that for a functional of the following type…
The functional derivative is calculated by….
So for this particular energy functional
Note that $\frac{\delta E}{\delta y(x)}$ is a functional derivative (rather than an ordinary derivative), with units of a force per unit length (rather than a force). Finally, we can calculate the total elastic force on a discrete membrane segment of length $\mathrm{\Delta}x$, by integrating this force density over the length of a segment (within which the spatial derivatives of $y$ do not vary).
where ${\sigma}_{\text{eff}}=\sigma h\mathrm{\Delta}x$ and ${\kappa}_{\text{eff}}=\kappa h\mathrm{\Delta}x$.
Detailed Derivations II: 2D thermal fluctuations of actin filaments
In this section, we characterize the thermal fluctuations of an actin filament with a given length and persistence length (and the associated bending modulus). We first decompose the fluctuations into their (small) endtoend and (larger) sidetoside fluctuations, and then further perform a wavelength decomposition on the sidetoside fluctuations. We next apply the equipartition theorem, giving $\frac{1}{2}{k}_{B}T$ to each independent bending mode, to determine the fluctuation magnitude of each of these modes. We then determine the average total sidetoside and endtoend fluctuation magnitudes, allowing us to calculate effective bending coefficients for the sidetoside (${\kappa}_{\u27c2}$) and endtoend (${\kappa}_{\parallel}$) fluctuations. These effective bending energies will then be used in Detailed Derivations III to determine the force of a filament on a membrane segment.
Determining the wavelengthdependence of thermal fluctuations
Assume we have a filament of length $L$, which, in the absence of thermal fluctuations, points vertically upward in the direction $\widehat{i}$. Due to thermal fluctuations, the actin polymer will fluctuate all along its length, as well as along the directions of both the short and long axis of the filament. At any point $s$ along the polymer, we can define the local position vector $\overrightarrow{r}(s)$, and the local tangent vector of the polymer $\overrightarrow{t}(s)=\frac{\partial \overrightarrow{r}}{\partial s}$. Given a persistence length l_{p}, and thus a bending modulus $B={l}_{p}{k}_{B}T$, the bending energy ${E}_{bend}$ of a filament is defined as…
where $\beta =\frac{1}{{k}_{B}T}$ and $\overrightarrow{t}(s=0)=\widehat{i}$. If we assume the thermal fluctuations are small ($\frac{L}{{l}_{p}}<<1$), then we can write the position vector along the polymer as…
Then…
And…
Because $\overrightarrow{t}$ is the unit tangent vector, we know $\overrightarrow{t}\cdot \overrightarrow{t}=1$. This allows us to solve for $\frac{\partial \delta}{\partial s}$ in terms of $\frac{\partial J}{\partial s}$.
Plugging this back into our equation for ${E}_{bend}$…
By decomposing $J(s)$ into its wavemodes, we can determine the relative amplitudes ${A}_{n}$ of thermal fluctuations at different lengthscales, where $n$ refers to the specific wavemode. This choice of wavemode decomposition oscillates around $\sqrt{2}{A}_{n}$ from 0 to $2\sqrt{2}{A}_{N}$ ($\u27e8{A}_{n}\u27e9=0$), where the filament is pinned at zero at $s=0$ and open at the other end.
Because…
Now we can solve the integral for the parts of the function that contain $L$. For the part where $n=m$…
And for $n\ne m$…
where $N=\frac{L}{{l}_{p}}$. By the equipartition theorem ($\beta {E}_{bend}=1/2$ for 1 degree of freedom), we can determine the Fourier coefficients.
We can now calculate a few useful integrals. From the equipartition theorem and Appendix Equation 26:
Because…
Finally…
Because for a Gaussian distribution, the 4th moment is related to the variance in the following way…
and
Determining the average endtoend retraction
The the new effective length of the filament is $I(s=L)$ can be found by integrating $\widehat{i}$ component of the tangent vector $\overrightarrow{t}$ over the arc length of the filament $s$.
From Appendix Equation 28…
Determining the endtoend retraction fluctuations
From Appendix Equation 28 and Equation 30…
Determining the sidetoside fluctuations
Deriving effective stretching constants
We know by the equipartition theorem that for each degree of freedom, $\frac{1}{2}{k}_{B}T=\frac{1}{2}\kappa \u27e8{x}^{2}\u27e9$. So for the $\widehat{x}$ direction…
And for the $\widehat{i}$ direction…
This gives us an effective bending energy as a function of actin filament tip position….
where ${\overline{\kappa}}_{\u27c2}=\frac{3}{{L}^{2}N}$ and ${\overline{\kappa}}_{\parallel}=\frac{16}{{L}^{2}{N}^{2}}$.
Detailed Derivations III: The force of an actin filament on the membrane
In this section, we take the effective sidetoside and endtoend filament bending energies calculated in Detailed Derivations II to determine the force of a filament on a membrane segment. We start by converting from the coordinate system used in Detailed Derivations II (in the frame of the filament long axis) to the reference frame of the membrane segment. We then use the bending energies to evaluate the partition function ${Z}_{conf}$ for a filament which is constrained by the membrane segment  and then use the partition function to determine the free energy $G$ of the filament. Finally, we evaluate the force of the filament on the membrane segment (${F}_{actin}$) as the increase in filament free energy obtained by an incremental movement of the membrane ${\frac{\partial G}{\partial \xi}}_{\xi =0}$.
Converting to the reference frame of the membrane segment
We derived this bending energy function in the reference frame of the filament. However, our partition function, ${Z}_{conf}$, will need to be integrated across the xy coordinate system used in the rest of the paper. The energy function can be extended to an arbitrary filament orientation in a ${\overrightarrow{u}}_{0}$ and tip position $\overrightarrow{r}$ relative to resting filament tip position ${\overrightarrow{r}}_{0}=(0,0)$ in the following way:
Taking a coordinate system centered around the equilibrium filament tip position and aligned with the membrane normal, a filament lying at angle $\theta $ relative to the membrane normal will have filament orientation vector ${\overrightarrow{u}}_{0}=sin(\theta )\widehat{x}+cos(\theta )\widehat{y}$ and filament tip position $\overrightarrow{r}=x\widehat{x}+y\widehat{y}$ will having the following perpendicular and parallel displacements…
and…
We arrive at the following bending energy as a function of the filament angle, relative to the membrane normal.
The x and ycomponents have been separated for ease of integration in the next section.
Simplification of the bending energy function
Currently, the exponential being integrated is of the form ${\int}_{\mathrm{\infty}}^{\mathrm{\infty}}{\int}_{0}^{\mathrm{\infty}}{e}^{a{x}^{2}+bxyc{y}^{2}}dydx$. We can rewrite the integral in the form ${\int}_{\mathrm{\infty}}^{\mathrm{\infty}}{\int}_{0}^{\mathrm{\infty}}{e}^{{\left(dx+fy\right)}^{2}g{y}^{2}}dydx$, to take advantage of the fact that ${\int}_{\mathrm{\infty}}^{\mathrm{\infty}}{e}^{{(x+c)}^{2}}dx={\int}_{\mathrm{\infty}}^{\mathrm{\infty}}{e}^{{x}^{2}}dx$. To do this, we can complete the square:$f(x)=a{x}^{2}+bx+c=a{\left(x+\frac{b}{2a}\right)}^{2}+c\frac{{b}^{2}}{4a}$
Now looking just at the ${y}^{2}$ portion, we can determine an effective bending coefficient in the $\widehat{y}$ direction.
Because $\overline{\kappa}}_{\perp}<{\overline{\kappa}}_{\parallel$, ${\kappa}_{\text{Eff}}$ has a maximum at $\theta =0$ and a minimum at $\theta =\frac{\pi}{2}$. In the limit where $\overline{\kappa}}_{\perp}<<{\overline{\kappa}}_{\parallel$, we find
where ${\kappa}_{\text{Eff}}=\frac{3{l}_{p}{k}_{B}T}{{L}^{3}si{n}^{2}(\theta )}$, agreeing with previous models (Mogilner and Oster, 1996, BiophysJ), except for the fact that we arrive at a prefactor of 3 rather than 4, as we took into account all of the bending modes of the filament, rather than assuming bending of the filament lies along an arc of constant curvature  a difference also arrived at by Dickinson and colleagues (Dickinson, Caro, and Purich, 2004, BiophysJ).
This gives us a final bending energy…
Calculation of the partition function
Using these energies, we can determine the partition function, ${Z}_{conf}$, summing over all possible positions of the filament. In this case, the filament can fluctuate freely along the xaxis, but cannot fluctuate past the membrane position along the yaxis. Here we define the membrane position $d$ relative to the filament tip.
If we then have a change of variables, where we center the system at d, the equation becomes…
where y_{0} is the equilibrium filament position, measured relative to the membrane surface (the filament pokes through the membrane for ${y}_{0}<0$).
Calculation of the force of actin on a membrane
Upon an infinitesimal membrane position perturbation $\xi $, the equilibrium filament tip position becomes ${y}_{0}\to {y}_{0}+\xi $. The partition function becomes…
and the derivative of the partition function with respect to the perturbation is (by Leibniz’s rule)…
Plugging Appendix Equation 51 and Equation 53 into Appendix Equation 5, we get the following force:
Note that this equation for the force assumes that the pointed (nongrowing) ends of the filaments (as well as branches) are rigidly fixed to a stiff and immobile actin network.
Detailed Derivations IV: Numerical approximations
In this section, we describe the various numerical approximations made in the development of our computational model.
Membrane stretch/bend forces
For membrane stretch/bend force calculations, the 2nd and 4th order spatial derivatives were calculated using a central finite difference approximation with 8th order accuracy.
Force of actin on the membrane
Many equations used in this model required performing numerical calculations of Gaussians over halfspace. For this purpose, we used the error function. The denominator of Appendix Equation 54 can be rewritten in terms of the error function in the following way…
The error function appears in the equation as $1+\text{erf}(z)$, which can give inaccurate calculations when $z<<0$ and $1+\text{erf}(z)\approx 0$. Therefore, in the low $z$ regime ($z<0$), we replaced $1+\text{erf}(z)$ with $\text{erfc}(z)$. (Because $\text{erf}(z)=\text{erf}(z)$, $\text{erf}(z(z<0))=\text{erf}(z(z<0))$). It follows that $1+\text{erf}(z(z<0))=1\text{erf}(z(z<0))=\text{erfc}(z(z<0))$. The following cases are listed below for the relevant equations.
When the error function calculations failed (e.g., produced values of 0 or $\mathrm{\infty}$), variable precision accuracy was used.
Probability of adding a monomer
We used similar error function approximations to calculate the probabilities that thermal fluctuations of the filament open up a gap between the filament tip and the membrane of sufficient size to add a monomer: (See Detailed Derivations IV: Numerical approximations, Force of actin on the membrane for rational on using the erf and erfc functions.)
When the error function calculations failed (e.g., produced values of 0 or $\mathrm{\infty}$), variable precision accuracy was used.
Choice of timestep
The timestep was chosen to capture the fastest dynamics in the system, which could either be the membrane stretch/bend relaxation, or actin network growth dynamics. This was implemented as min($\frac{0.1}{{r}_{max}^{actin}}$,$\frac{1}{{r}_{max}^{memb}}$). This timestep was calculated specifically for each simulation, depending on the parameters chosen. The timescales of the actin dynamics are set by the rates of polymerization, depolymerization, branching, and capping. To estimate the timescales of relaxation for membrane elastic forces, we calculate the relaxation timescales of the membrane fluctuating under Brownian thermal forces. (We do this because we do not a priori have an analytical theory describing the shape profile of the actin dynamics. Importantly, the fact that our simulations are not affected by the timestep (Fig. S3) provide evidence that we are sufficiently resolving all system dynamics using this approximation to set the timestep.) For a membrane at low Reynolds number, we have the following equation of motion.
The random Brownian force is Gaussian distributed with mean $\mu =0$ and variance ${\sigma}^{2}=2\gamma {k}_{B}T$ and is defined by its autocorrelation function
We are ultimately interested in the time evolution of wave mode solutions to this equation, so we first take the Fourier transform in order to solve for the amplitude of each wavemode. Using the following properties of the Fourier transform:
we can rewrite the equation in Fourier space as…
where $\alpha =\frac{2\pi k}{L}$. The random Brownian force in Fourier space is Gaussian distributed with mean $\mu =0$ and variance $2\gamma {k}_{B}T\mathrm{\Delta}x$ as defined by its autocorrelation function. For a discrete system, a Fourier transform with normalization $\frac{1}{\sqrt{N}}$, where N is the number of discrete membrane segments in this case, preserves the variance and standard deviation of a vector of normally distrbuted random values. Here, our normalization is $\frac{1}{\sqrt{L}}$ in addition to integrating over the fixed segment length $\mathrm{\Delta}x$, giving a final normalization of $\sqrt{\frac{\mathrm{\Delta}x}{N}}$ for the Brownian force $\mathrm{\Xi}(k,t)$. The factor of $\frac{1}{\sqrt{N}}$ preserves the variance, leaving the variance changed only by the multiplicative factor of ${(\sqrt{\mathrm{\Delta}x})}^{2}=\mathrm{\Delta}x$.
We can solve this equation using a Laplace transform.
and then in inverse Laplace transform
Using the properties
and
and
and
Now that we have $Y(k,t)$, we can solve for for the timeautocorrelation function $\u27e8Y({k}^{\prime},t)Y(k,0)\u27e9$.
We know that $\u27e8Y({k}^{\prime},0)Y(k,0)\u27e9$ is nonzero only if $k={k}^{\prime}$ .
by a variant of Fubini’s theorem
and a property of the Delta function
So finally…
From this equation, the rate of relaxation due to membrane elasticity is…
Data availability
Analysis and modeling code for this paper is available on the Theriot lab Gitlab: <https://gitlab.com/theriot_lab/leadingedgestabilityinmotilecellsisanemergentpropertyofbranchedactinnetworkgrowth> under the MIT license. Figure data are available in the Source Data files. The large size of the raw video microscopy data (865 GB of image files in the Open Microscopy Environment OMETIFF format) and the associated analyzed data (320 GB) prohibits their upload to a public repository. The complete raw and analyzed data files for one example experimental dataset and one example simulated dataset (corresponding to the data shown in Fig. 1af and Fig. 2cj, respectively) are available on Figshare <https://figshare.com/projects/Leading_edge_stability_in_motile_cells_is_an_ emergent_property_of_branched_actin_network_growth/132878>. Code to analyze this data are publicly available on Gitlab as noted above. Requests for additional raw or analyzed data should be sent to the corresponding author by email. Data will be made available in the form of a hard drive shipped by mail. There are no restrictions on who may access the data.
References

The ActinBased Nanomachine at the Leading Edge of Migrating CellsBiophysical Journal 77:1721–1732.https://doi.org/10.1016/S00063495(99)770189

Membraneproximal Factin restricts local membrane protrusions and directs cell migrationScience (New York, N.Y.) 368:1205–1210.https://doi.org/10.1126/science.aay7794

Elastic Modeling of Biomembranes and Lipid BilayersAnnual Review of Physical Chemistry 59:685–712.https://doi.org/10.1146/annurev.physchem.59.032607.093550

Growth of Branched Actin Networks against ObstaclesBiophysical Journal 81:1907–1923.https://doi.org/10.1016/S00063495(01)758420

Growth Velocities of Branched Actin NetworksBiophysical Journal 84:2907–2918.https://doi.org/10.1016/S00063495(03)700186

Geometry and Physics of WrinklingPhysical Review Letters 90:074302.https://doi.org/10.1103/PhysRevLett.90.074302

Organelles  Understanding noise and heterogeneity in cell biology at an intermediate scaleJournal of Cell Science 130:819–826.https://doi.org/10.1242/jcs.181024

The Role of Actin Polymerization in Cell MotilityAnnual Review of Physiology 53:585–605.https://doi.org/10.1146/annurev.ph.53.030191.003101

WASP family proteins and formins compete in pseudopod and blebbased migrationThe Journal of Cell Biology 217:701–714.https://doi.org/10.1083/jcb.201705160

Dynamics of Entangled Polymer SolutionsMacromolecules 9:587–593.https://doi.org/10.1021/ma60052a011

Neutrophil migration in infection and wound repair: Going forward in reverseNature Reviews. Immunology 16:378–391.https://doi.org/10.1038/nri.2016.49

Excitable Signal Transduction Networks in Directed Cell MigrationAnnual Review of Cell and Developmental Biology 33:103–125.https://doi.org/10.1146/annurevcellbio100616060739

WASP and SCAR are evolutionarily conserved in actinfilled pseudopodbased motilityThe Journal of Cell Biology 216:1673–1688.https://doi.org/10.1083/jcb.201701074

Neutrophil‐like HL‐60 cells expressing only GFP‐tagged β‐actin exhibit nearly normal motilityCytoskeleton (Hoboken, N.J.) 77:181–196.https://doi.org/10.1002/cm.21603

Mechanical feedback between membrane tension and dynamicsTrends in Cell Biology 22:527–535.https://doi.org/10.1016/j.tcb.2012.07.005

Neutrophil recruitment and function in health and inflammationNature Reviews. Immunology 13:159–175.https://doi.org/10.1038/nri3399

Myosin light chain kinase regulates cell polarization independently of membrane tension or Rho kinaseThe Journal of Cell Biology 209:275–288.https://doi.org/10.1083/jcb.201409001

Profiling cellular morphodynamics by spatiotemporal spectrum decompositionPLOS Computational Biology 14:e1006321.https://doi.org/10.1371/journal.pcbi.1006321

Chemotaxis in neutrophillike HL60 cellsMethods in Molecular Biology (Clifton, N.J.) 571:167–177.https://doi.org/10.1007/9781607611981_11

Cell motility driven by actin polymerizationBiophysical Journal 71:3030–3045.https://doi.org/10.1016/S00063495(96)794961

Design Principles of Length Control of Cytoskeletal StructuresAnnual Review of Biophysics 45:85–116.https://doi.org/10.1146/annurevbiophys070915094206

From solution to surface to filament: actin flux into branched networksBiophysical Reviews 10:1537–1551.https://doi.org/10.1007/s1255101804695

What’s all the noise about developmental stochasticityDevelopment (Cambridge, England) 138:601–607.https://doi.org/10.1242/dev.059923

Loading history determines the velocity of actinnetwork growthNature Cell Biology 7:1219–1223.https://doi.org/10.1038/ncb1336

Cellular motions and thermal fluctuations: the Brownian ratchetBiophysical Journal 65:316–324.https://doi.org/10.1016/S00063495(93)81035X

Rate constants for the reactions of ATP and ADPactin with the ends of actin filamentsThe Journal of Cell Biology 103:2747–2754.https://doi.org/10.1083/jcb.103.6.2747

Molecular Mechanisms Controlling Actin Filament Dynamics in Nonmuscle CellsAnnual Review of Biophysics and Biomolecular Structure 29:545–576.https://doi.org/10.1146/annurev.biophys.29.1.545

Actin, a Central Player in Cell Shape and MovementScience (New York, N.Y.) 326:1208–1212.https://doi.org/10.1126/science.1175862

Direct measurement of the lamellipodial protrusive force in a migrating cellThe Journal of Cell Biology 174:767–772.https://doi.org/10.1083/jcb.200601159

Crawling Toward a Unified Model of Cell Motility: Spatial and Temporal Regulation of Actin DynamicsAnnual Review of Biochemistry 73:209–239.https://doi.org/10.1146/annurev.biochem.73.011303.073844

The structural basis of actin filament branching by the Arp2/3 complexThe Journal of Cell Biology 180:887–895.https://doi.org/10.1083/jcb.200709092

A review of models of fluctuating protrusion and retraction patterns at the leading edge of motile cellsCytoskeleton (Hoboken, N.J.) 69:195–206.https://doi.org/10.1002/cm.21017

Dynamics of capping protein and actin assembly in vitro: Uncapping barbed ends by polyphosphoinositidesThe Journal of Cell Biology 135:169–179.https://doi.org/10.1083/jcb.135.1.169

Performance of a population of independent filaments in lamellipodial protrusionBiophysical Journal 95:1393–1411.https://doi.org/10.1529/biophysj.107.125005

Membrane tension and cytoskeleton organization in cell motilityJournal of Physics. Condensed Matter 27:273103.https://doi.org/10.1088/09538984/27/27/273103

A phagocytic cell line markedly improves survival of infected neutropenic miceJournal of Leukocyte Biology 78:338–344.https://doi.org/10.1189/jlb.0205072

Activation of nucleation promoting factors for directional actin filament elongation: Allosteric regulation and multimerization on the membraneSeminars in Cell & Developmental Biology 24:267–271.https://doi.org/10.1016/j.semcdb.2013.01.006

Analysis of the actinmyosin II system in fish epidermal keratocytes: Mechanism of cell body translocationThe Journal of Cell Biology 139:397–415.https://doi.org/10.1083/jcb.139.2.397

Ultrastructure of protrusive actin filament arraysCurrent Opinion in Cell Biology 25:574–581.https://doi.org/10.1016/j.ceb.2013.04.003

The actin cytoskeleton and actinbased motilityCold Spring Harbor Perspectives in Biology 10:1–21.https://doi.org/10.1101/cshperspect.a018267

Electron tomography reveals unbranched networks of actin filaments in lamellipodiaNature Cell Biology 12:429–435.https://doi.org/10.1038/ncb2044

Orientational Order of the Lamellipodial Actin Network as Demonstrated in Living Motile CellsMolecular Biology of the Cell 14:4667–4675.https://doi.org/10.1091/mbc.e02100630

Structure of arp2/3 complex in its activated state and in actin filament branch junctionsScience (New York, N.Y.) 293:2456–2459.https://doi.org/10.1126/science.1063025

ParticleTracking Microrheology of Living Cells: Principles and ApplicationsAnnual Review of Biophysics 38:301–326.https://doi.org/10.1146/annurev.biophys.050708.133724

IBAR domain proteins: linking actin and plasma membrane dynamicsCurrent Opinion in Cell Biology 23:14–21.https://doi.org/10.1016/j.ceb.2010.10.005
Decision letter

Alphee MichelotReviewing Editor; Institut de Biologie du Développement, France

Anna AkhmanovaSenior Editor; Utrecht University, Netherlands

Alphee MichelotReviewer; Institut de Biologie du Développement, France
Our editorial process produces two outputs: i) public reviews designed to be posted alongside the preprint for the benefit of readers; ii) feedback on the manuscript for the authors, including requests for revisions, shown below. We also include an acceptance summary that explains what the editors found interesting or important about the work.
Decision letter after peer review:
Thank you for submitting your article "Leading edge maintenance in migrating cells is an emergent property of branched actin network growth" for consideration by eLife. Your article has been reviewed by two peer reviewers, including Alphee Michelot as the Reviewing Editor and Reviewer #1, and the evaluation has been overseen by Anna Akhmanova as the Senior Editor.
The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.
Essential revisions:
1) The main point raised by Reviewer #2 is closely related to the last point raised by previous Reviewer #3. The reviewer found that the analysis of Figure A1, A2 and A3 in the appeal to the other Journal shows that the redblue stripe pattern of Figure 1b is indeed not intrinsic to the system. Presumably the features of Figure 1b reflect the length scale over which the curvature was calculated. Thus we agree that you should include the analysis of these figures in the manuscript and discuss these results to clarify these concerns.
2) Some explanations on spatial resolution would help (Reviewer #1). The spatial scales accessible with the experiments and those described in the model are very different. It would be useful to explain to which scales the experimental data remain informative for the understanding of these molecular mechanisms.
Reviewer #1 (Recommendations for the authors):
I have only one comment below which should be addressed before publication:
A missing point of this manuscript is the spatial resolution of these experiments based on phasecontrast microscopy. I find hard to believe that measured fluctuations arise from stochasticity in actin growth at the level of monomer addition, as claimed by the authors. This does not necessarily precludes all the conclusions from this study, but these seemingly exaggerated claims do not help readers understand the link between the experimental data and the model.
In a revised manuscript, I would suggest the authors to add a section addressing this problem of spatial resolution. This section should provide estimates of typical shape fluctuations that the experimental system is capable of detecting (or not), together with estimates of minimal network growth that can be detected. Next, the authors should indicate over which range of parameters their model can be compared with experimental data, and over which range of parameters their model becomes simply predictive.
Reviewer #2 (Recommendations for the authors):
To resolve my concern, the authors could plot the autocorrelation functions in Figure 1e at longer times and calculate Figure 1g using a smoothing length scale sufficiently larger than 7 microns.
https://doi.org/10.7554/eLife.74389.sa1Author response
Essential revisions:
1) The main point raised by Reviewer #2 is closely related to the last point raised by previous Reviewer #3. The reviewer found that the analysis of Figure A1, A2 and A3 in the appeal to the other Journal shows that the redblue stripe pattern of Figure 1b is indeed not intrinsic to the system. Presumably the features of Figure 1b reflect the length scale over which the curvature was calculated. Thus we agree that you should include the analysis of these figures in the manuscript and discuss these results to clarify these concerns.
2) Some explanations on spatial resolution would help (Reviewer #1). The spatial scales accessible with the experiments and those described in the model are very different. It would be useful to explain to which scales the experimental data remain informative for the understanding of these molecular mechanisms.
We have now addressed both of these essential revisions in our updated manuscript. For our response to point (1), please see our additions to the Results (lines 121128), Methods (lines 427432, 452454, 480482, 497512), and three new supplemental figures (Figure 1 —figure supplement 35). For our response to point (2), please see our expanded discourse on differences in spatial resolution between the experimental and theoretical aspects of our work in the Results (lines 9299, 144145, 153154, 286291) and Discussion (lines 318325).
Reviewer #1 (Recommendations for the authors):
I have only one comment below which should be addressed before publication:
A missing point of this manuscript is the spatial resolution of these experiments based on phasecontrast microscopy. I find hard to believe that measured fluctuations arise from stochasticity in actin growth at the level of monomer addition, as claimed by the authors. This does not necessarily precludes all the conclusions from this study, but these seemingly exaggerated claims do not help readers understand the link between the experimental data and the model.
In a revised manuscript, I would suggest the authors to add a section addressing this problem of spatial resolution. This section should provide estimates of typical shape fluctuations that the experimental system is capable of detecting (or not), together with estimates of minimal network growth that can be detected. Next, the authors should indicate over which range of parameters their model can be compared with experimental data, and over which range of parameters their model becomes simply predictive.
We have now done this. Please see our expanded discourse on this topic in the Results and Discussion.
Reviewer #2 (Recommendations for the authors):
To resolve my concern, the authors could plot the autocorrelation functions in Figure 1e at longer times and calculate Figure 1g using a smoothing length scale sufficiently larger than 7 microns.
We have now done this. Please see additions to the Results, Materials and methods, and in three new supplemental figures (Figure 1 —figure supplement 35).
[Editors' note: we include below the reviews that the authors received from another journal, along with the authors’ responses.]
Reviewer #1:
In their manuscript Garner and Theriot take a perturbation free observational approach in combination with theoretical modelling to address the question why lamellipodia have a smooth structure. They employ neutrophil like cells migrating under agarose pads and obtain movies with very high time resolution and find that the edge fluctuates locally. Autocorrelations show that the local fluctuations show a characteristic relaxation pattern that can be modelled using minimal components of actin polymerisation dynamics and membrane interaction. A striking result was that the relaxation could be reproduced even without including lateral mechanical interaction with the plasma membrane. The authors propose a model where the relaxation is an emergent property of the actin network geometry. Due to the 70{degree sign} branchpoint connections the filaments slide sideways and thereby the adjacent segments catch up, as the filaments sliding in from the lagging segments cannot catch up with the leading membrane and fail to spawn nucleation sites (the only occur in association with the membrane). Varying branching angles they find that the 70{degree sign} Arp2/3 branch is optimal for providing this behavior.
This is very elegant and deep work and an extremely useful framework to think about the crucial question why the branch junction is conserved around 70{degree sign}. Some questions should be addressed before publication.
1. Neutrophils under agarose are certainly develop lamellipodia, however their hydrostatic components seem quite strong and even in the phase contrast movie of the untreated cell a bleb seems to occur. It would be good to control that fluctuations are really driven by polymerisation and not by hydrostatics. Blebbistatin treated cells should still migrate and it would be good to see that the fluctuations are not affected.
The reviewer raises an excellent point that hydrostatics may play a role in mediating leading edge shape fluctuations, and the strength of this effect remains an exciting area of study. Unfortunately, blebbistatin would almost certainly perturb actin turnover and thus free monomer concentration (which would lead to similar phenotype to that observed under Latrunculin B treatment), in addition to hydrostatic pressure. This experiment would therefore be unlikely to shed light on the relative roles of actin polymerization and hydrostatic pressure on leading edge fluctuation dynamics. Future modeling efforts which explicitly include hydrostatic pressure might allow predictions of fluctuation behavior in hydrostaticsdriven or polymerizationdriven regimes, but we consider this to be outside the scope of the current work. We have added clarification in the discussion that other cellular properties, including hydrostatics, may be acting in cells to additionally modulate leading edge shape dynamics, and that subsequent additions to the model could be used to further dissect the relative strength of these effects in diverse contexts (lines 299302, 310).
2. As neutrophils have little to no retrograde flow of actin the following question arises: if the agarose pad has some structural irregularity the load on the lamellipodial tip would fluctuate as well and forcevelocity relation could also drive fluctuations. It would be satisfying to see that the model also holds in a situation without agarose overlay. The Theriot lab should have plenty of movies of e.g. fish keratocytes to check if the fluctuations are comparable.
We thank the reviewer for their intriguing question regarding the effects of agarose and comparative analyses of leading edge fluctuations in HL60s cells and keratocytes. We were sufficiently intrigued by the reviewer’s suggestion that we collected a new, highspeed imaging dataset on keratocytes for the first time. We are happy to report that keratocytes, in the absence of an agarose pad overlay, indeed exhibit leading edge fluctuations qualitatively and quantitatively similar to those of HL60 cells under confinement. We have added a new supplementary companion figure to Figure 1 (Figure S3) describing leading edge fluctuations in keratocytes, and discuss these results in the main text (lines 123125). We believe this fully addresses the reviewer’s question on the effect of agarose on the fluctuations in HL60s. We note the model itself is not cell typespecific, and so should generally apply regardless of experimental condition. Therefore it was indeed satisfying to see that the results of the model are congruent with experimental data from both HL60 cells and keratocytes.
3. The authors propose that filaments laterally sliding in from the lagging segments cannot keep up with the membrane, fall behind and do not contribute to new branches. Therefore the lagging segments that receive additional filaments sliding in from the leading segment catch up while the initially leading segment thins out and slows. This is a plausible result of the model – but I have difficulties understanding that the filaments laterally arriving from the lagging parts cannot keep up. Following a force velocity curve they should speed up and keep contact with the membrane instead of falling behind and being capped.
The reviewer states that filaments following a forcevelocity curve should speed up to keep contact with the membrane. However, this is not the case for lamellipodial actin networks, according to established theoretical and experimental work. Rather, forcevelocity curves of branched actin networks as measured by AFM in keratocyte lamellipodia (Prass, Jacobson, Mogilner, and Radmacher, 2006) and in vitro (Parekh, Chaudhuri, Theriot, and Fletcher, 2005) exhibit a largely loadindependent velocity in low load regimes (<1 nN/µm^{2}). This is largely due to the fact that dendritically branched actin networks increase in actin density (via branching) in response to increased force to maintain a constant amount of force per filament. Lamellipodial actin networks, in which filaments are much shorter than their persistence length and much stiffer than the plasma membrane, operate in this low load regime (Mogilner and Oster, 1996), consistent with electron microscopy images of the keratocyte lamellipodium, in which filaments are observed to remain straight, rather than bend, under the load of the leading edge membrane (Svitkina and Borisy, 1999). Therefore we do not expect significant changes in the polymerization velocity due to small variations in the load (and in particular, we do not expect lagging filaments to speed up). Rather, we expect actin density to decrease in response to a reduced load to maintain a constant polymerization velocity, which is consistent with our model. Thus there is no discrepancy between our model and the known forcevelocity relationship.
4. Can the sideway sliding be directly seen as lateral propagation in in the curvature plots (Figure 1b)? When the red area splits the resulting Y junction should have a characteristic angle. Is this the case and if yes is it in quantitative agreement with the model?
We thank the reviewer for their intriguing observation of lateral propagation in the curvature kymograph. The filament spreading invoked by the model and schematized in Figure 4j should, as far as we know, act uniformly in both directions rather than slide unidirectionally. Although sliding should not be unidirectional, actin density in cells tends to be highest in the center of the leading edge and decreases towards the edges. Thus we might imagine spreading could manifest as curvature events that slide from the center towards the sides of the cell. However, spreading in our model occurs at the singlefilament level. It remains an open area of research how this molecularscale spreading phenomenon should propagate to the larger lengthscales seen in the curvature kymograph – and further what implications that might have for lateral spreading of micronscale curvature events. We remain excited about this future direction but find it to be outside of the scope of the current work.
Reviewer #2:
In this work, the authors study the protrusion of membranes by polymerizing a branched actin network. In particular, the address the question of how a flat membrane configuration can be stable in presence of fluctuations due to molecular noise in the actin polymerization and branching process. They find that high actin density regions, which would locally deform the leading edge, are dissolved because newly nucleated filaments (daughter filaments) branch off under an angle from existing filaments (mothers). Their simulations indicate that the experimentally observed branching angle of 70 degrees is optimal for stabilizing a flat leading edge.
The problem as such is interesting – analyzing the fluctuations of a membrane that is pushed by a polymerizing actin network – and also the mechanism by which fluctuations relax is nicely worked out but the authors, but I have several issues with this work that, in my opinion, preclude its publication in the Journal. First of all, the main motivation for the study is "Because the average cell shape remains constant over time, there must be some form of feedback acting on leading edge curvature to sustain stable lamellipodial growth." (ll 83) This is not necessarily true – and sure enough the authors find that no feedback is necessary, which they sell as a "surprising revelation" (l 278). There is a large body of literature in physics that studies surface growth and there it is very common to see that the roughness of the surface saturates with time, see for example, the ballistic deposition model. Also in the context of a polymerizing branched network this is at least implicitly known, see Weichsel and Schwarz, PNAS 107, 6305 (2010).
While the reviewer correctly states that many general physical models of surface growth exhibit a surface roughness that saturates with time, this is not necessarily the case for all systems, and certainly has not been established for lamellipodial actin. In particular, several features of lamellipodial actin are known to lead to nonlinearities in network growth, which might amplify stochastic fluctuations. First, dendritic branching is an autocatalytic process which can lead to explosive growth. Second, the “deposition” of actin is dependent on the velocity of the flexible membrane surface it is pushing, in a manner which imparts hysteresis to the system (Parekh et al., 2005). Finally, the geometry of actin network growth exhibits nontrivial interactions with the motion and shape of the leading edge membrane. It is far from obvious whether these generic surface growth models would capture each of these complexities. The reviewer cites Weichsel and Schwarz, PNAS 107, 6305 (2010), but this work only investigates actin network growth against a flat rigid surface – in which “surface roughness” is not defined and therefore does not provide a relevant point of comparison. Further, our spreadingnull simulations (Figure 4fi) explicitly demonstrate a case in which actin network growth exhibits instability – providing a direct counterexample to the reviewer’s expectation from generic surface growth models. For all of these reasons, the existence of a class of generic surface growth models that exhibit stability in no way precludes the significance and novelty of our results.
Furthermore, decades of work in cell biology and cell biophysics have investigated the need for extrinsic feedback mechanisms to maintain leading edge stability. This has been considered a major open question in the field – until now. For instance, several influential reviews and perspectives dedicated to the subject have identified the potential roles of membrane tension (DizMuñoz, Fletcher, and Weiner, 2013; Sens and Plastino, 2015), curvaturesensing proteins (Zegers and Friedl, 2015), and a competition for membraneassociated free monomers (Mullins, Bieling, and Fletcher, 2018) in maintaining a single, stable lamellipodium – in addition to several influential papers outlining mechanisms based on forcedependent directional filament branching (Risca et al., 2012) and WAVE complex selforganization (Pipathsouk et al., 2019). The current work shows that none of these mechanisms are required for the generation and maintenance of a stable lamellipodium, and thus represents a major advance in the field as well as a surprising and novel finding.
We have made substantial changes to the Introduction which more clearly set up why our major finding is surprising and a significant contribution to the field, including rearranging the logical flow of the Introduction, as well as: 1) strengthening our language regarding the decades of prior literature assuming some sort of regulatory or mechanical must be required for leading edge stability (lines 4347); and 2) highlighting sources of known nonlinearities in actin network growth that likely amplify stochastic fluctuations (lines 5359).
My second issue is that the argument the authors give for the 70 degrees branching angle are not convincing to me. First of all, the minimum in relaxation times is not very pronounced. Second, why should evolution single out the network structure that minimizes surface roughness? Would there be any way to test the evolutionary advantage experimentally?
We assume that the reviewer’s reference to the “minimum in relaxation times” is meant to refer to the minimum in leading edge shape fluctuation amplitude, seen in Figure 5h, as the manuscript never references a minimum relaxation time. It is difficult to state in absolute terms whether this minimum is or is not “pronounced”. While it is unknown quantitatively how the effect size on leading edge fluctuation amplitude translates to a functional advantage in cell motility, and further unknown how a functional advantage on the cellular level conveys a fitness advantage to the organism, there is no reason to assume the downstream effects of this minimum are small. In general, a small perturbation to one component of a physical system can have large effects the behavior of the system as a whole – and this certainly holds for many cellular processes (e.g. signaling cascades). Specifically in the case of evolutionary selection, it has certainly been established that even phenotypes with very small relative fitness advantages can come to dominate a population, and thus act as a driver of evolution (Crow and Kimura, 1970; Imhof and Schlotterer, 2001) Therefore it is not possible to conclude that an apparently shallow minimum is not functionally or biologically significant. Further, this minimum observed in our work is the only evidence of any kind proposed in the literature addressing the striking evolutionary conservation of the Arp2/3mediated branching angle, and thus represents a significant finding for the field.
The reviewer’s next question regarding why evolution might minimize surface roughness is an excellent one, and one that could only even be asked now that our work has identified this important potential mechanism. As discussed in the previous paragraph, the relationship between molecular mechanism, cellular function, and organismal fitness is extremely nontrivial and currently poorly understood – and thus beyond the scope of our work, but we certainly hope that our thought provoking results will inspire future work to address this important question.
My third issue is with the way the model is set up, notably the membrane part. The membrane bending energy is incorporated in a strange way. By themselves, the springs connecting the membrane segments do not lead to a bending energy. Only by constraining the motion of the segments along the ydirection, it comes into play in some way. But then this is also membrane tension, so it's quite confusing to me.
Furthermore, if I understood correctly, thermal fluctuations of the membrane are not taken into account.
These also depend on the membrane tension, which would thus affect the roughness of the surface. Continuing with the segments, their size corresponds to which physical parameter? Probably something like the membrane's persistence length. When changing the membrane tension, then the persistence length changes and thus the size of the segments. All these points remained unclear to me and I think that the description of the membrane is not physically sound.
The reviewer is incorrect to state that bending energy comes into play within our model implementation “only by constraining the motion of the segments along the ydirection”. Rather, the bending energy is incorporated as the 2^{nd} derivative of the membrane position (curvature) in the energy functional (Supplementary text equation 2). For clarity we now provide a succinct physical interpretation of our energy functional, including where bending appears, to the Supplementary Text corresponding to Equation 2. Perhaps the reviewer interpreted “connected by springs” more literally than intended, as the springs were meant to represent elastic coupling (i.e., resisting both extension and curvature) between the segments. For clarity, we have revised the text to replace “connected by springs” with “coupled elastically to each other” (line 137).
The reviewer is also incorrect to assume that membrane fluctuations would affect “surface roughness”. It has long been established that actin filament fluctuations, rather than membrane fluctuations, drive membrane motion and monomer incorporation into the network (Mogilner and Oster, 1996). Therefore we do not expect thermal fluctuations of the membrane to significantly affect leading edge shape fluctuations – thus it is physically sound to exclude them. We have added clarification to this extent in the supplementary text associated with Equation 1.
Finally, the reviewer is correct that the membrane segment size must be sufficiently short to resolve the smallest length scales at which there is significant bending. Indeed, we took this into account in developing our model. For a freely floating membrane fluctuating under Brownian motion, the appropriate segment size would be at least as small as the membrane’s persistence length to correctly resolve the membrane dynamics. We chose a membrane segment size well below the persistence length of this membrane “strip”. However, what we have modeled is not a membrane fluctuating under Brownian motion, but a membrane being driven by a very stiff actin network. It is possible, even likely, that the system as a whole would have a different effective persistence length than a membrane alone (e.g., the actin growth dynamics might suppress bending of the leading edge below a certain curvature). Lacking analytical theory describing this effective persistence length of the actinmembrane system, the only way to know that the relevant length scales have been sufficiently resolved is to vary the segment size and confirm that there is no effect on the dynamics. We directly demonstrated that the membrane segment size is sufficiently small by varying the membrane segment size and show quantitatively that the leading edge fluctuations are not affected (Figure S5). This was referenced in the main text (lines 165166), and now is described more explicitly in the Methods section on Modeling (lines 478479). Therefore we are convinced that our treatment of the membrane is physically sound.
In conclusion, in my opinion, the authors need to tone down their claims and provide more physical justification of their model. With such revisions, they could submit the work to a more specialized journal.
Minor comments
ll. 79 "reminiscent of an excitable system in which local instabilities in the leading edge emerge, grow, and then relax" This is not sufficient for describing it as an excitable system. Fluctuations always grow initially and – if the system is stable – relax, but not all systems subject to fluctuations are excitable.
The reference to excitable systems has been removed (line 88).
ll. 80 "Consistent with the idea that these fluctuations arise from stochasticity in actin growth, kymograph analysis of curvature and velocity (Figure 1bc) showed that relatively longlived curvature fluctuations are formed by the continual timeintegration of very shortlived (subsecond, submicron) velocity fluctuations." I fear that I do not understand this sentence.
It is unclear which part of the sentence confused the reviewer, but we thank them for drawing our attention to some potential points of confusion which we have clarified in the text as outlined below.
As our manuscript was written for a broad audience, one function of this sentence is to act as a succinct reminder (for the benefit of readers lacking a strong physics or mathematics background) of the mathematical relationship between velocity, position, and curvature. As velocity is defined as the first timederivative of position, the membrane position results from the continual timeintegration of the velocity fluctuations (i.e., the velocity fluctuations add up over time to determine the membrane positional fluctuations). As curvature is defined as the second spatialderivative of the membrane position, curvature fluctuations also arise via a timeintegration of the velocity fluctuations. We have replaced the phrase “curvature fluctuations” with “shape fluctuations” to avoid inadvertently implying that curvature is directly determined as the timeintegral of the velocity (line 93).
Additionally, this sentence serves to emphasize that the small subsecond, submicron scales of the velocity fluctuations – visually appearing uncorrelated to and acting at different scales than the membrane curvature – are consistent with fluctuations arising from stochastic actin polymerization at the level of monomer addition. Several clarifying points along these lines have been added to this sentence (lines 9194).
Figure 1d: I am not sure what to learn from this sketch.
Figure 1d is intended to provide readers without a strong background in physics or mathematics with an intuition for the relaxation of shape fluctuations, a pictographic description of the wavelength, λ, and relaxation timescale, τ, as used throughout the paper. Based on our experience presenting this work to audiences with a biological background, we have found that it is helpful to include a cartoon representation of this commonly observed relationship between wavelength and relaxation timescale for elastic systems (see figure caption for Figure 1d).
Figure 1e,f: What are the wavelengths associated with the chosen Fourier modes?
The spatial frequencies are equally spaced between 0.220.62 µm^{1} in 0.06 µm^{1} intervals, corresponding to wavelengths between 1.64.5 µm. This information has been added to the legend for Figure 1.
Ll 109 "This general trend is shared by many physical systems with linear elastic constraints, such as idealized membranes and polymers freely fluctuating under Brownian motion, but can be contrasted with systems that have a dominant wavelength, as in the case of buckling or wrinkling of materials under compression." I do not understand the relation between buckling and fluctuations. The buckling instability can be discussed within the frame of linear elasticity (LandauLifshitz, vol 7, {section sign}21).
Buckling is invoked here merely as an example to give readers without a strong physics background some intuition for the kinds of systems that exhibit a dominant wavelength. We are not claiming any specific connection between buckling and fluctuations.
L 138 "To our great surprise" Depends on your experience… I would suggest to take it out.
As we have described in our response to the reviewer’s first major concern, the result that a simple actin network growth model can reproduce stable leading edge fluctuations, in the absence of any extrinsic sources of feedback, is indeed a surprising result in the field. Therefore we believe that, while our choice is a stylistic preference, the use of this phrase correctly indicates to readers less familiar with actin biology that this result was highly unexpected given the state of the field.
Ll 151 "It should be noted that the generation of the simulated data in Figure 2kl did not involve any curvefitting (and therefore no free parameters that could be fit) to the experimentallymeasured autocorrelation dynamics." I do not understand. The relaxation times and fluctuation amplitudes should depend on parameters. I think what you want to state is that the parameters you have taken from the literature give results that agree decently with the experimental data.
We did not mean to imply that the simulation results do not depend on the input parameters. Rather we state that we did not use curvefitting on the experimentally measured fluctuations to parameterize the simulations. In modeling of biological systems, it is common to parameterize models by fitting to the experimental data, so we needed to make clear this is not the case. We have clarified these points in the text (line 168).
Figure 3 "Minimal model correctly predicts response of HL60 cells to drug treatment." The model does not predict anything – and "correct" is not defined. You show that the trends and the range of fluctuation amplitude and decay rate are similar.
We thank the reviewer for prompting us to clarify how we defined predictive power in this context. The title for Figure 3 uses the shorthand “correctly predicts”, which is intended to convey that the outputs of the simulation are congruent with experimental observations under conditions that had not been tested prior to model development. (Specifically, our model originally developed using parameters measured for unperturbed cells is able to recapitulate qualitative and quantitative features of the experimental data for latrunculintreated cells, demonstrating how a reduction in free monomer concentration should affect quantitative parameters such as leading edge fluctuation amplitudes and timescales.) This is obviously too long for a figure title; however, to address the reviewer’s concern we now explicitly define how we assess the model’s predictive power in the main text associated with Figure 3 (lines 175177).
Ll 233 "only a handful of mother filament angles (θf) survive until the end of the simulation" I do not understand – there is a welldefined average of θ_{f} and some variation around this angle; what does "handful of mother filament angles" mean? In Figure 5a there is no distinction between mother and daughter filaments, right?
We thank the reviewer for drawing our attention to this potential point of confusion. The reviewer is correct that there is no distinction between mother and daughter filaments in Figure 5a. In this sentence we used “mother filament[s]” to identify the filaments which seeded the simulation. Because there are a finite set of discrete filament angles represented in this randomized starting population, the “surviving” filament angles are descended from a few “ideal” seed filament angles. We agree this is slightly confusing, as any filament which creates a branch at any time in the simulation is technically a mother filament to its branch. For clarity, we have replaced “mother filament angles” with “the initial filament angles” (line 251).
Ll 234 "successful filament angles are narrowly (…) distributed around one half of the branching angle" What does "narrowly" mean? 'Narrow' is a relative notion, right?
We agree with the reviewer that “narrow” is a relative notion. However, the referenced distribution is narrow by almost all reasonable standards. The distribution is narrow relative to the initial population of seed filament angles, relative to the range of total possible angles (0°180°), relative to the branching angle, and relative to the resolution of experimental measurements of lamellipodial filament angles by scanning electron microscopy (SEM).
Figure 5h,l, q: add labels to the xaxes.
It is unclear what the reviewer is referring to. The xaxes on Figure 5 h, l, and q are labeled. Possibly the reviewer is referring to the insets of Figure 5h,l,p? The figure legend states “(h,l,p) […] Insets have identical xaxes to main panels.” Adding duplicate, identical xaxis labels to the insets would only make the figure more cluttered and difficult to interpret.
Ll 293 "The defining characteristic and major advance of our model was the explicit inclusion of both the evolutionary dynamics of the actin network and its interaction with the twodimensional geometry of the leading edge." But this the membrane is not essential, see also Weichsel and Schwarz op cit, who produced a similar network dynamics but without an explicit membrane.
As the reviewer correctly states, Weichsel and Schwarz 2010 did not include an explicit membrane – and therefore could not possibly have explored the interaction between network dynamics and membrane bending leading to the twodimensional leading edge geometry (as we have done here). We have been careful to cite the landmark papers describing the evolutionary aspects of the actin network (Maly and Borisy, 2001; Schaus and Borisy, 2008; Schaus, Taylor, and Borisy, 2007), and we do not claim that our reproduction of these results is novel. Rather, it is the interaction between these evolutionary dynamics and the twodimensional leading edge shape which is novel in our work.
Ll 368 "for which we used a conservative cutoff of 10 sec." I am not sure what this means. Could you present a graph showing the steady state has been reached after 10s?
Our model’s approach to steady state can be seen in three out of our five main text figures (Fid 2dh, Figure 4ae, and Figure 5a), as well as three out of our five supplementary figures (Figure S4cg, Figure S5cg, Figure S6cg), and is directly referenced in the main text: “Nascent leading edges reach steady state values for filament density, filament length, filament angle, membrane velocity, and (most importantly) membrane fluctuation amplitude within seconds, a biologically realistic time scale (Figure 2dh).” (lines 154156). However, for good measure we have added a reference to Figure 2dh at the end of the line referenced by the reviewer (line 423).
Ll 400 "We expanded this formalism to include all fluctuations of the filament along the filament's short and long axes." In ll 397 you write that filaments were fixed in position, now you write about fluctuations along the filament direction and perpendicular to it. Could you clarify?
As described in the work by Mogilner and Oster (1996) cited in this section, “filaments apply force to the membrane segments according to the untethered Brownian ratchet formalism^{5}”. In this formalism, filaments are assumed to be rigidly attached to the actin network at their branch points. The filament thus fluctuates around this rigid attachment. This rigid attachment is fixed in position (in the lab frame of reference), such that filaments do not undergo retrograde flow (i.e. translation of the filament position opposite the direction of migration) or translational diffusion. Clarification of how filaments are translationally fixed in position at one end (lines 461462), and freely fluctuating at the other end (lines 465466) is now included.
Supplementary Material
I do not understand Equation (3). F_{S/B} should be $\frac{{\text{\delta E}}_{S/B}}{\text{\delta y}}$– I am not sure, why you integrate over x. Furthermore, why is the integral over x replaced by multiplication with Δ x? Doesn't the integral extend over the whole membrane? – Now I have seen Equation (14), which clarifies things. I would suggest that you add the integration boundaries to Equation (3) to make it clearer.
As the reviewer correctly understood after having read the appendix referenced in the text corresponding to Equation 3, $\frac{{\text{\delta E}}_{S/B}}{\text{\delta y}}$ is a functional derivative (as opposed to an ordinary derivative) having units of force per unit length, and must be integrated over the size of a discrete membrane segment in order to determine the force on said segment. We have added the integration boundaries as suggested, and added clarifications in the text associated with Eq 3 and Eq 14.
App. II: I still do not understand how connecting segments by springs lead to a bending energy. I guess that it is somehow related to the fact that the segments can only move along the ydirection.
Please see our response to the reviewer’s third major concern (above), where we address this question.
Reviewer #3:
This paper describes analysis and modeling of leading edge fluctuations in migrating cells driven by a branched Arp2/3 lamellipodial network. A stochastic model shows how branching contributes to shape stability, thus explaining the measured spectrum and dynamics of leading edge fluctuations. Analysis of the model as a function of branching angle suggests that the Arp2/3 branching angle might be selected to smooth lamellipodial shape. The topic of actin driven cell motility will be of general interest and appropriate for the Journal. The authors provide new ideas to a big field of research, including Fourier analysis of leading edge fluctuations, which to my knowledge is a novel approach. The modeling methods and model design seem valid and appropriate to me and the paper is well written.
I have however several questions on issues that are unclear to me and are central to the main results of the paper. The following comments focus on the validity of the analysis of experimental images to derive the observed spectrum of fluctuations, the interpretation of the model results, as well as whether the derived optimal Arp2/3 complex value may be dependent on the chosen parameter values.
Major comments
1) I have some questions and concerns regarding the analysis of leading edge behavior in Figure 1.
Firstly, I am worried that looking at deviations of the leading edge contour away from a smoothed profile over 7 microns (as was done in this work) may be erasing the dynamics that cause the most significant fluctuations in the system. The behavior of the large wavelength fluctuations in Figure 1eh may depend on this smoothing. Figure 1B shows a nearly periodic redblue curvature pattern with a scale of order 4 microns that persists over 1015 sec, a time over which the cell advances by a distance of order the size of the lamellipodium. Such a nearly periodic pattern would lead to a peak at the corresponding wavelength, which is not seen however in Figure 1g.
We appreciate the reviewer’s remarks on the details of our spectral analysis, which prompted us to incorporate several additions and clarifications to the manuscript text regarding the methods – including: 1) a detailed description on how the kymographs were generated, 2) a discussion on why kymographs fall short in characterizing finescale leading edge fluctuations, and 3) a discussion on how or whether the background subtraction is expected to affect the measured fluctuation amplitudes and dynamics. We believe these changes, in addition to our pointbypoint response below, fully address the reviewer’s concerns.
With regards to the reviewer’s concern that the background subtraction removes “the most significant fluctuations in the system”, we emphasize that the most obvious fluctuations seen in the curvature kymograph are not necessarily the most significant. Indeed, the most obvious features in a kymograph will inevitably be the largest features occurring over the longest timescales. In contrast, the explicit goal of this work was to understand the finescale fluctuations in the leading edge shape that we discovered, and then took great pains to measure, using highspeed, high resolution imaging (Figure 1a). In fact, our work nicely demonstrates the limitations of kymograph analysis, where the largerscale fluctuations are arbitrarily emphasized over finescale fluctuations – making obvious the need for a more quantitative treatment that encompasses all spatial scales, as exemplified by our Fourier autocorrelation analysis. We include the kymographs in Figure 1 to serve as a visual aid for readers who are most familiar with this kind of analysis to obtain a qualitative sense of the data, before quickly moving on to our more indepth definitive analysis. To clarify these points for the reader, we have added to the Methods: 1) a new section describing the generation of the kymographs (lines 389392) and 2) a discussion of our motivation for using spectral decomposition, including the limitations of kymograph analysis (lines 394401).
In response to the reviewer’s concern that the background subtraction is removing oscillatory dynamics in the data, we do not believe a lack of oscillations in our Fourier analysis means that there is something wrong with our background subtraction technique. Most importantly, we do not see a mechanism by which spatial smoothing/filtering should remove temporal oscillations in our dataset. Even if for some reason the background subtraction were to partially suppress the 4 µm signal, it would do so uniformly in time, and therefore would maintain oscillations in the Fourier autocorrelation results. To elucidate these points, we have expanded our description of the spectrum analysis in the Methods to include an explicit discussion on the expected effects (or lack thereof) of the background subtraction on the measured fluctuation amplitudes and rates (lines 442450).
Finally, we do not see evidence of strong temporal oscillations, as opposed to longtimescale fluctuations which appear and then decay, in our curvature kymograph. In general, it is nontrivial to prove the existence of oscillations in any noisy signal (e.g., see Liu, Oh, Peshkin, and Kirschner, 2020); certainly visual inspection of a kymograph does not suffice. We acknowledge some readers familiar with the field may be primed to look for oscillations due to their familiarity with protrusionretraction patterns – which are not directly relevant to this system (please see our answer to the reviewer’s next comment for further discussion on this topic).
Therefore, our curvature kymograph is very much consistent with our Fourier autocorrelation results, and we have no reason to doubt our use of background subtraction for fluctuations larger than 7 µm.
Leading edge fluctuations over these spatial and temporal scales (microns, lamellipodium turnover time) have been considered in several prior works, including modeling studies by Zimmerman and Falcke, Ryan and Vavylonis, Gov group and others. The fact that there is no connection made to this body of work makes me wonder if the authors may have missed relevant questions answered or raised by these prior works.
We thank the reviewer for drawing attention to a potential point of confusion regarding the previously wellstudied case of lamellipodial protrusionretraction dynamics and their comparison to the leading edge shape fluctuations presented in this paper. The modeling studies mentioned by the reviewer only considered lamellipodia that undergo protrusionretraction cycles, which is a fundamentally different mode of motility compared to that studied in this work – operating with distinct actin dynamics and biomechanics, exhibiting qualitatively different shape dynamics, and evolving over markedly slower (hundredfold longer) timescales.
In particular, these cells cannot be undergoing protrusionretraction cycles, because the leading edge does not exhibit any evidence of retraction (i.e. the velocity is never negative in Figure 1c). A line has been added to the caption for Figure 1c pointing out the uniformly positive velocity shown in the kymograph. In addition, cells that undergo protrusionretractionbased motility, such as fibroblasts, typically migrate inefficiently, moving with average speeds ten to a hundredfold slower than HL60 cells, and the oscillatory period of their protrusionretraction dynamics is typically on the order of tens of minutes, while the fluctuations studied in this work maximally last tens of seconds. Therefore the reviewer is incorrect to say that the referenced prior publications considered leading edge fluctuations over the same temporal scales as our current work.
For all of these reasons, the leading edge fluctuations identified and explored in this work represent a novel class of leading edge dynamics, which are mechanistically distinct from protrusionretraction cycles. To highlight the novelty of the observed dynamics, we have added a sentence to the manuscript clarifying how our measured leading edge fluctuations are distinct from the previously wellstudied case of protrusionretraction cycles, including citations to some of this literature (lines 8991).
I think that the autocorrelation functions in Figure 1e should be shown at longer times to check that they decay to zero (or perhaps those with long λ may show weak periodicity as Figure 1b would suggest). It's not clear to me why none of the examples in Figure 1e do not decay significantly over time.
The reviewer’s assumption that the autocorrelation functions as plotted in Figure 1e should decay to zero is not correct, for reasons we now make more clear in the Methods (lines 432436) and discuss here. The amplitude being plotted in Figure 1e is the complex magnitude of the autocorrelation values (this is stated in the figure legend as well as in the Methods). The Fourier transform and resulting autocorrelation values are complex, having both real and imaginary parts of the form A = a + ib. The complex magnitude of this value, sqrt(a^{2} + b^{2}), is most representative of the total autocorrelation. We of course expect the autocorrelation to decay down to some noise window, which may be wavelengthdependent. Because we are plotting the complex magnitude (which is always positive), the noise window appears as a levelingoff of the autocorrelation curve to some nonzero background value. Therefore we do not expect the complex magnitude to decay down to zero, but rather to some finite background value, which is indeed what we find. The timescales shown in Figure 1e are adequate to demonstrate this decay to background levels.
The leading velocity kymograph in Figure 1c is very noisy. Does it show anything beyond noise around an average value and was there any effort to calculate the velocity over varying time intervals to smooth it out and possibly reveal a pattern?
As alluded to in our answer to the reviewer’s first comment, the kymographs are not central to any conclusions in the paper, but rather serve as a visual aid for readers to obtain a qualitative sense of the data. We did not investigate the velocity kymograph further, because, as mentioned in the paper (lines 9192), a noisy velocity kymograph is consistent with the very simple expectation that leading edge shape fluctuations are driven by stochastic fluctuations in actin network growth. Further, as kymographs are ultimately a qualitative method, we chose to analyze the fluctuations using a superior quantitative method: the Fourier autocorrelation analysis that forms the basis for the rest of the manuscript.
2. A model with fluctuating membrane pushed by a dendritic network has been developed by Schaus et al. (Ref. 45) with a follow up paper in Biophys J. A comparison/discussion I feel is necessary, even if the similarities end up superficial. These authors also varied the branching angle.
While the two papers referenced by the reviewer included a model similar to the current work (in that they incorporated stochastic filament polymerization near a flexible membrane), these papers were mainly focused on selforganization of the actin network, largely ignoring the membrane shape. For example, Schaus et al. (2007) designed their model to study the selforganization of actin filament orientation, and Schaus and Borisy (2008) sought to understand limitations on the optimal “performance” and loadsharing of filaments from an engineering perspective. Importantly, these papers simulated only small segments of a leading edge that were < 2 µm – and thus barely touched on questions of leading edge stability and shape fluctuations. Therefore there are no results that could be compared to the current work, and any comparative discussion would potentially mislead and/or confuse the reader. However, we now acknowledge the conceptual similarity of our approach to the referenced publications in the introduction of our model (lines 134135). Of course, we also reference this work where it is most relevant to our results: in the context of previous literature discussing selforganization of actin filament orientation (lines 246248, 253255) – as was recognized by the reviewer – and now with added emphasis on line 255.
3. Intrinsic stability of branched actin networks seems to be a main message of the paper and it's discussed in Figure 4. But why would one expect instability and why was it a surprise that the simple model would lead to stable fluctuations (page 9, beginning of last paragraph)? From Figure 4g, this it seems the authors may have in mind a density/curvature instability, perhaps with a characteristic wavelength (which is yet what appears to be seen in Figure 1b as I mention on comment 1).
We thank the reviewer for raising this point, which prompted us to make several changes to the Introduction that clarify why our result is surprising. For our complete response, please see our reply to the first comment by Reviewer #2, who had a similar question.
Also, given the emphasis on spectral analysis in Figure 13, in Figure 4 I was expecting to see how membrane tension and bending influence the amplitude and decay constants of the modes (but this is not shown). Do membrane tension and bending energy play any role in the fluctuations, and at what values? I would anticipate the spectral analysis is most relevant in Figure 4, to show the wavelengthdependence of the stability mechanism in Figure 4j (and possibly its relationship to the spectrum of Figure 2).
We are glad the reviewer takes interest in our spectral analysis. Indeed, we also believe a dissection of the exact shape of the fluctuation amplitude and decay rate curves, as well as their dependence on model parameters such as membrane tension, would be an exciting area for future study. However, the major findings of this work are the emergent stability of lamellipodial actin and the optimal angle of the Arp2/3complex – neither of which depend on the details of the spectral analysis. Rather, we use the spectral analysis to quantitatively describe the experimentally observed fluctuations and to validate our model by comparison to the experimental data. While further interrogation of the spectra would be interesting, it is not required for any of our main conclusions and is thus beyond the scope of the work.
With respect to the stability mechanism suggested in Figure 4j, it is not trivial to intuit the associated fluctuation spectra from first principles, and therefore any comparison to the results shown in Figure 2 would be speculative. One could develop an analytical model reflecting our proposed mechanism, which could then predict the wavelengthdependence of the proposed stability mechanism and be compared to our simulations. However, we believe our identification of this novel mechanism, on its own, is already a major advance in the field. The additional development and validation of an associated analytical model would be a substantial undertaking and is left for future studies.
4. (Related to the previous question) Is there a qualitative or quantitative interpretation of the behavior of the amplitude of fluctuations versus wavelength in Figure 2k? For example, does this graph show a crossover between two (powerlaw?) regimes from high amplitude large wavelengths confined by the cell/simulation size to suppressed fluctuations at short wavelengths as a result of the stability mechanism in Figure 4j? I am partly asking this trying to understand the long wavelength behavior as function of the leading edge length in Figure S5 where I would have expected the longest wavelength fluctuations to increase with system size.
As discussed in our answer to the reviewer’s previous comment, we believe an indepth investigation into the exact shape of the fluctuation spectrum curves is beyond the scope of the current work and is not required to inform the major conclusions of our study. The observed curves could possibly be determined by any number of complex interactions between the relevant variables (polymerization velocity, actin density, actin filament orientation, membrane tension/bending rigidity, etc.) – and so even a qualitative interpretation would be speculative at best. Therefore we did not include any speculation in the text of the manuscript.
With regards to the reviewer’s specific reference to Figure S5 (now Figure S6), the fluctuation amplitudes for each wavelength are not, and should not be, constrained by the length of the leading edge. If we simulate a longer leading edge, we can resolve longer wavelength fluctuations, but the amplitude of any particular mode should not be affected by the size of the system (length of the leading edge). The size of the system is a simulation parameter, not a biological parameter, so it should not affect the result (i.e., a simulated small patch of membrane should behave identically to an equivalently sized portion of a larger simulated patch of membrane). Figure S5 (now Figure S6) is therefore a control demonstrating that our results correctly do not depend on the leading edge length simulated. To clarify these points for the reader, we have added an explanation for each of the control simulations to the Methods (lines 477482).
5. Various modeling works have shown that the filament angle distribution at steady state have different preferred orientation, depending on both the branching angle as well as the ratio of protrusion to polymerization rate (which determine if a polymerizing branch can catch up with the leading edge). I assume that as the branching angle is varied in Figure 5, the system may change preferred steady state orientations. As I understand, the analysis in Figure 5 showing that the geneticallyencoded Arp2/3 angle is optimal was performed for fixed polymerization rate. Wouldn't the optimal value change depending on the polymerization rate?
We are excited that the reviewer is curious about the mechanisms by which the Arp2/3mediated branching angle optimally smooths fluctuations – a topic that we believe will be of considerable interest for followup studies based on this work.
The reviewer brings up an interesting question about the role of the polymerization rate in leading edge fluctuations. However, we do not believe there is a trivial relationship between the polymerization rate and the optimal branching angle. While the relationship between branching angle, the preferred filament orientation, and polymerization velocity is well understood, how or whether any of these variables might affect the leading edge fluctuation amplitude remains entirely unclear. If we cannot a priori predict how a change in the polymerization velocity should affect the fluctuation amplitude for a given branching angle, then we certainly cannot predict how a change in the polymerization velocity would affect the minimum amplitude across all possible branching angles. Given these complexities, any investigation or discussion of the effect of the polymerization rate – however interesting – would be beyond the scope of the current work.
In any case, the polymerization rate used in our simulations was specifically chosen to be within a biologicallyrelevant regime and wellgrounded in experimental data. The rate was taken directly from published experimental measurements (Table S1), and produces a simulated leading edge velocity nearly identical to that of our experimentally measured cell velocities (Figure 1c, Figure 2g). Therefore the 7080° branching angle is indeed optimal for the experimentallymeasured cell velocities and actin polymerization rates – which should be most relevant for the question we were most interested in, regarding the evolutionary conservation of the branching angle in cells. We emphasize that this is the first mechanistic evidence of any kind regarding the optimality of the evolutionarily conserved angle, and thus represents a major advance that stands on its own. Any explorations into the roles of the polymerization rate on the optimal branching angle would be outside the scope of this work.
Minor points
6. In Table S1, the actin filament persistence length seems very small, 1 micron. This is an important parameter for this Brownian ratchet model. Is this a typo?
We thank the reviewer for their keen eye – however, our use of 1µm as the persistence length is not a typo. We use the same value for the persistence length as Mogilner and Oster (1996), and as measured by analysis of thermal undulations of actin filaments in Käs et al. (1996) (Käs et al., 1996; Mogilner and Oster, 1996). While we acknowledge that a persistence length of ~1020 µm is a more commonly cited range, all of the known experimentally measured values have been acquired using in vitro systems, rather than in cells. Given recent evidence that microtubules exhibit a persistence length ~100 times smaller in cells than those observed in vitro (Pallavicini et al., 2014), we believe it is not unreasonable to assume a value on the lower end of in vitro experimental measurements for Factin.
7. An outline or plan of the Appendix material would be helpful to understand the purpose of various calculations involved as one starts reading it. The math can also greatly be reduced, for example I don't think showing the details of a simple integration such as Equation (55) is needed, this takes away from focusing on more important equations such as 56.
We thank the reviewer for their suggestion to add an outline of the Appendix material. An introductory paragraph has been added at the beginning of each Appendix in the Supplementary Text.
We respectfully disagree that the detailed derivations should be removed. As this manuscript was written for the broad readership of the Journal, we believe that not all readers will find Supp. Equation 55 to be simple – particularly those without a strong background in math and physics, many of whom will not be familiar with the error function, for example. We find it is far easier for those having familiarity with the material to skip over seemingly trivial details, than for someone with no familiarity to fill in the gaps themselves.
8. A simulation movie at shorter time intervals than the one provided might be more helpful: the current one shows snapshots that differ significantly from frame to frame.
The purpose of this video was to display the simulations at a similar timescale to that shown for the experimental data – in particular, to demonstrate the long timescale stability of the simulated leading edges. Therefore we show the full simulated time (100s) at 0.5 s resolution. Displaying at faster time resolution would quickly make the file size prohibitively large.
9. It may be helpful to list somewhere the values of the wavelength of the modes in Figure 1e,f and elsewhere. At first reading I understood these curves to be the first few modes rather than a sampling of different modes.
We thank the reviewer for prompting us to clarify the wavemodes shown in Figure 1ef. The spatial frequencies are equally spaced between 0.220.62 µm^{1} in 0.06 µm^{1} intervals, corresponding to wavelengths between 1.64.5 µm. This information has been added to the legend for Figure 1.
10. Figure 5: is there a reason to choose branching angles of 17.5, 35 and 70 deg in 5c while 17.5, 70 and 140 degrees in 5d?
A different set of angles were displayed in Figure 5c, as compared to Figure 5d, as the two graphs are making a different set of comparisons. Figure 5c shows not only how the dominant angle is always one half of the branching angle, θ_{br}, but also show additional resonant peaks at θ_{br}/2 + θ for branches growing off of filaments lying at the dominant angle. Given the axes limits of 090° (filaments pointing towards the membrane), these resonant peaks can only be seen for branch angles less than 60°. In contrast, Figure 5de shows that a 70° branching angle optimally minimizes actin density variability, for which angles both above and below 70° must be displayed for comparison.
11) In case this helps, I noticed the κ_{eff} above Equation (49) has the numerical prefactor in Dickinson Biophys J 87:2838 (2004) that they mention is different to Mogilner and Oster, though I am not sure if the calculation here is the same.
We thank the reviewer again for their keen eye in drawing our attention to the fact that Dickinson et al. also arrived at this prefactor. As discussed in Dickinson et al., Mogilner and Oster arrive at a different prefactor because they approximate filament bending to occur only along an arc of constant curvature, rather than treating the full spectrum of bending fluctuations – as performed in this work and in Dickinson et al. We have added a discussion to this extent to the Supplementary Text associated with Equation (49).
References:
Crow, J. F., and Kimura, M. (1970). An introduction to population genetics theory. An introduction to population genetics theory. New York, Evanston and London: Harper and Row, Publishers.
DizMuñoz, A., Fletcher, D. a., and Weiner, O. D. (2013). Use the force: Membrane tension as an organizer of cell shape and motility. Trends in Cell Biology, 23(2), 47–53. https://doi.org/10.1016/j.tcb.2012.09.006
Imhof, M., and Schlotterer, C. (2001). Fitness effects of advantageous mutations in evolving Escherichia coli populations. Proceedings of the National Academy of Sciences, 98(3), 1113–1117. https://doi.org/10.1073/pnas.98.3.1113
Käs, J., Strey, H., Tang, J. X., Finger, D., Ezzell, R., Sackmann, E., and Janmey, P. A. (1996). Factin, a model polymer for semiflexible chains in dilute, semidilute, and liquid crystalline solutions. Biophysical Journal, 70(2 I), 609–625. https://doi.org/10.1016/S00063495(96)796303
Liu, X., Oh, S., Peshkin, L., and Kirschner, M. W. (2020). Computationally enhanced quantitative phase microscopy reveals autonomous oscillations in mammalian cell growth. Proceedings of the National Academy of Sciences of the United States of America, 117(44), 27388–27399. https://doi.org/10.1073/pnas.2002152117
Maly, I. V., and Borisy, G. G. (2001). Selforganization of a propulsive actin network as an evolutionary process. Proceedings of the National Academy of Sciences, 98(20), 11324–11329. https://doi.org/10.1073/pnas.181338798
Mogilner, A., and Oster, G. (1996). Cell motility driven by actin polymerization. Biophysical Journal, 71(6), 3030–3045. https://doi.org/10.1016/S00063495(96)794961
Mullins, R. D., Bieling, P., and Fletcher, D. A. (2018). From solution to surface to filament: actin flux into branched networks. Biophysical Reviews, 10(6), 1537–1551. https://doi.org/10.1007/s1255101804695
Narayanan, J., Xiong, J. Y., and Liu, X. Y. (2006). Determination of agarose gel pore size: Absorbance measurements vis a vis other techniques. Journal of Physics: Conference Series, 28(1), 83–86. https://doi.org/10.1088/17426596/28/1/017
Pallavicini, C., Levi, V., Wetzler, D. E., Angiolini, J. F., Benseñor, L., Despósito, M. A., and Bruno, L. (2014). Lateral motion and bending of microtubules studied with a new singlefilament tracking routine in living cells. Biophysical Journal, 106(12), 2625–2635. https://doi.org/10.1016/j.bpj.2014.04.046
Parekh, S. H., Chaudhuri, O., Theriot, J. A., and Fletcher, D. A. (2005). Loading history determines the velocity of actinnetwork growth. Nature Cell Biology, 7(12), 1219–1223. https://doi.org/10.1038/ncb1336
Pipathsouk, A., Brunetti, R., Town, J., Breuer, A., Pellett, P., Marchuk, K., … Weiner, O. (2019). WAVE complex selforganization templates lamellipodial formation. BioRxiv. https://doi.org/10.1101/836585
Prass, M., Jacobson, K., Mogilner, A., and Radmacher, M. (2006). Direct measurement of the lamellipodial protrusive force in a migrating cell. Journal of Cell Biology, 174(6), 767–772. https://doi.org/10.1083/jcb.200601159
Risca, V. I., Wang, E. B., Chaudhuri, O., Chia, J. J., Geissler, P. L., and Fletcher, D. a. (2012). Actin filament curvature biases branching direction. Proceedings of the National Academy of Sciences, 109(8), 2913–2918. https://doi.org/10.1073/pnas.1114292109
Schaus, T. E., and Borisy, G. G. (2008). Performance of a population of independent filaments in lamellipodial protrusion. Biophysical Journal, 95(3), 1393–1411. https://doi.org/10.1529/biophysj.107.125005
Schaus, T. E., Taylor, E. W., and Borisy, G. G. (2007). Selforganization of actin filament orientation in the dendritic nucleation/arraytreadmilling model. Proceedings of the National Academy of Sciences of the United States of America, 104(17), 7086–7091. https://doi.org/10.1073/pnas.0701943104
Sens, P., and Plastino, J. (2015). Membrane tension and cytoskeleton organization in cell motility. Journal of Physics: Condensed Matter, 27(27), 273103. https://doi.org/10.1088/09538984/27/27/273103
Svitkina, T. M., and Borisy, G. G. (1999). Arp2/3 complex and actin depolymerizing factor/cofilin in dendritic organization and treadmilling of actin filament array in lamellipodia. Journal of Cell Biology, 145(5), 1009–1026. https://doi.org/10.1083/jcb.145.5.1009
Zegers, M. M., and Friedl, P. (2015). Translating Membrane Tension into Cytoskeletal Action by FBP17. Developmental Cell, 33(6), 628–630. https://doi.org/10.1016/j.devcel.2015.06.006
https://doi.org/10.7554/eLife.74389.sa2Article and author information
Author details
Funding
National Science Foundation
 Rikki M Garner
Howard Hughes Medical Institute
 Julie A Theriot
Washington Research Foundation
 Julie A Theriot
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We dedicate this work to A B Savinov, who was in preparation alongside this manuscript and was born during revisions. We are exceedingly grateful to E Labuz for generously isolating and preparing fish epidermal keratocytes used for the experiments presented in Figure 1—figure supplement 6. We also kindly thank E F Koslover and A J Spakowitz for their thoughtful advice on model conceptualization and data analysis, and E F Koslover, A Mogilner, N M Belliveau, P Radhakrishnan, and A Savinov for helpful comments on the manuscript. Funding: Howard Hughes Medical Institute to JAT, Washington Research Foundation to JAT, Gerald J Lieberman Fellowship to RMG, NSF Graduate Research Fellowship to RMG.
Ethics
Experiments using zebrafish larvae were approved by the University of Washington Institutional Animal Care and Use Committee (protocol 442701).
Senior Editor
 Anna Akhmanova, Utrecht University, Netherlands
Reviewing Editor
 Alphee Michelot, Institut de Biologie du Développement, France
Reviewer
 Alphee Michelot, Institut de Biologie du Développement, France
Publication history
 Preprint posted: August 22, 2020 (view preprint)
 Received: October 2, 2021
 Accepted: March 9, 2022
 Accepted Manuscript published: March 11, 2022 (version 1)
 Version of Record published: April 22, 2022 (version 2)
 Version of Record updated: May 3, 2022 (version 3)
Copyright
© 2022, Garner and Theriot
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics

 3,191
 Page views

 413
 Downloads

 5
 Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading

 Cell Biology
Mitochondrial dysfunction has been reported in obesity and insulin resistance, but primary genetic mitochondrial dysfunction is generally not associated with these, arguing against a straightforward causal relationship. A rare exception, recently identified in humans, is a syndrome of lower body adipose loss, leptindeficient severe upper body adipose overgrowth, and insulin resistance caused by the p.Arg707Trp mutation in MFN2, encoding mitofusin 2. How the resulting selective form of mitochondrial dysfunction leads to tissue and adipose depotspecific growth abnormalities and systemic biochemical perturbation is unknown. To address this, Mfn2^{R707W/R707W} knockin mice were generated and phenotyped on chow and high fat diets. Electron microscopy revealed adiposespecific mitochondrial morphological abnormalities. Oxidative phosphorylation measured in isolated mitochondria was unperturbed, but the cellular integrated stress response was activated in adipose tissue. Fat mass and distribution, body weight, and systemic glucose and lipid metabolism were unchanged, however serum leptin and adiponectin concentrations, and their secretion from adipose explants were reduced. Pharmacological induction of the integrated stress response in wildtype adipocytes also reduced secretion of leptin and adiponectin, suggesting an explanation for the in vivo findings. These data suggest that the p.Arg707Trp MFN2 mutation selectively perturbs mitochondrial morphology and activates the integrated stress response in adipose tissue. In mice, this does not disrupt most adipocyte functions or systemic metabolism, whereas in humans it is associated with pathological adipose remodelling and metabolic disease. In both species, disproportionate effects on leptin secretion may relate to cell autonomous induction of the integrated stress response.

 Cell Biology
 Microbiology and Infectious Disease
Oral inflammatory diseases such as apical periodontitis are common bacterial infectious diseases that may affect the periapical alveolar bone tissues. A protective process occurs simultaneously with the inflammatory tissue destruction, in which mesenchymal stem cells (MSCs) play a primary role. However, a systematic and precise description of the cellular and molecular composition of the microenvironment of bone affected by inflammation is lacking. In this study, we created a single cell atlas of cell populations that compose alveolar bone in healthy and inflammatory disease states. We investigated changes in expression frequency and patterns related to apical periodontitis, as well as the interactions between MSCs and immunocytes. Our results highlight an enhanced selfsupporting network and osteogenic potential within MSCs during apical periodontitisassociated inflammation. MSCs not only differentiated towards osteoblast lineage cells, but also expressed higher levels of osteogenic related markers, including Sparc and Col1a1. This was confirmed by lineage tracing in transgenic mouse models and human samples from oral inflammatoryrelated alveolar bone lesions. In summary, the current study provides an indepth description of the microenvironment of MSCs and immunocytes in both healthy and disease states. We also identified key apical periodontitisassociated MSC subclusters and their biomarkers, which could further our understanding of the protective process and the underlying mechanisms of oral inflammatoryrelated bone disease. Taken together, these results enhance our understanding of heterogeneity and cellular interactions of alveolar bone cells under pathogenic and inflammatory conditions. We provide these data as a tool for investigators not only to better appreciate the repertoire of progenitors that are stress responsive but importantly to help design new therapeutic targets to restore bone lesions caused by apical periodontitis and other inflammatoryrelated bone diseases.