A mechanism with severing near barbed ends and annealing explains structure and dynamics of dendritic actin networks
Abstract
Single molecule imaging has shown that part of actin disassembles within a few seconds after incorporation into the dendritic filament network in lamellipodia, suggestive of frequent destabilization near barbed ends. To investigate the mechanisms behind network remodeling, we created a stochastic model with polymerization, depolymerization, branching, capping, uncapping, severing, oligomer diffusion, annealing, and debranching. We find that filament severing, enhanced near barbed ends, can explain the single molecule actin lifetime distribution, if oligomer fragments reanneal to free ends with rate constants comparable to in vitro measurements. The same mechanism leads to actin networks consistent with measured filament, end, and branch concentrations. These networks undergo structural remodeling, leading to longer filaments away from the leading edge, at the +/-35° orientation pattern. Imaging of actin speckle lifetimes at sub-second resolution verifies frequent disassembly of newly-assembled actin. We thus propose a unified mechanism that fits a diverse set of basic lamellipodia phenomenology.
Editor's evaluation
Although studied for decades, the molecular mechanisms involved in the assembly and remodeling of the lamellipodium still pose a number of questions, among which 1/ how are these networks progressively reorganized from short branched filaments to longer ones while maintaining angular order, 2/ by which mechanisms are actin filaments disassembled in these networks (depolymerization, fragmentation and/or "catastrophic" disassembly), and 3/ what is the importance and contribution of filament annealing? To address these questions, the authors develop one of the most detailed stochastic computational models to date. The model takes into account a large number of chemical reactions, including actin polymerization, depolymerization, filament branching by the Arp2/3 complex, capping, uncapping, severing, oligomer diffusion, annealing, and debranching. Close comparison of in silico and cellular actin networks allows them to evaluate the relative contribution of the different reactions. An important finding of this work is that frequent actin filament severing and annealing are phenomena that cannot be neglected to describe lamellipodial dynamics appropriately and although filament annealing in cells is not a new discovery, it is striking that it is not a negligible and inconsequential phenomenon in the cell, but contributes significantly to the reorganization of actin networks.
https://doi.org/10.7554/eLife.69031.sa0Introduction
The force for lamellipodial protrusions is provided by a dendritic network of actin filaments. This dynamic structure is driven by actin filament polymerization, branch generation by the Arp2/3 complex and regulation of filament elongation by capping protein (Pollard and Borisy, 2003; Watanabe, 2010; Blanchoin et al., 2014). Activated by nucleation promoting factors on the cell membrane, the Arp2/3 complex nucleates filament branches at an angle of approximately 70° from filaments that reach the leading edge. These elongating barbed ends add actin monomers from the cytoplasm to push against the cell membrane and generate force for membrane extension or for the retrograde flow of the whole dendritic actin network toward the cell center. This dendritic lamellipodia network structure, evident in electron micrographs of keratocytes (Svitkina et al., 1997) has been quantified by more recent electron tomograms near the leading edge, revealing the number of barbed ends, branches and filaments (Vinzenz et al., 2012; Mueller et al., 2017). Its characteristic pattern with filaments orientated primarily at ±35° with respect to the protrusion axis (Maly and Borisy, 2002; Schaub et al., 2007; Mueller et al., 2017; Koseki et al., 2019) has been interpreted by two-dimensional dendritic network models (Schaus et al., 2007; Maly and Borisy, 2002; Weichsel and Schwarz, 2010; Atilgan et al., 2005; Holz and Vavylonis, 2018).
In parallel to polymerization and branching at the leading edge, lamellipodia maintain their steady state through continuous disassembly and recycling of actin (Pollard and Borisy, 2003; Watanabe, 2010; Blanchoin et al., 2014). Extensive biochemical and biophysical studies have identified critical aspects of the kinetics and thermodynamics of this turnover process, with cofilin and hydrolysis of ATP bound to actin after polymerization, followed by Pi release, playing a central role. However the precise molecular mechanisms of actin turnover in cells have not been fully resolved (Danuser and Waterman-Storer, 2006; Carlsson, 2010; Carlier and Shekhar, 2017).
Single-Molecule Speckle (SiMS) microscopy of fluorescently labeled actin revealed that actin assembly into the dendritic network is transient and not limited to the leading edge (Yamashiro et al., 2018; Watanabe and Mitchison, 2002). In these SiMS experiments, actin subunits incorporated into the actin network appear as single molecule speckles while diffuse actin contributes to background fluorescence. In the lamellipodium, speckle disappearances occur within a few seconds after speckle appearances, a time which is relatively short compared to the time required for actin treadmilling through the entire lamellipodium. Since filament treadmilling cannot explain these dynamics, Miyoshi and Watanabe proposed the hypothesis of frequent filament severing near barbed ends, following by annealing of the oligomeric fragment (Miyoshi and Watanabe, 2013).
Consistent with the frequent severing near ends and annealing hypothesis, in vitro experiments show dissociation of filament fragments from ends of actin filaments in vitro, in the presence of cofilin and co-factors (Wioland et al., 2017; Kueh et al., 2008; Shekhar and Carlier, 2017; Andrianantoandro and Pollard, 2006; McCullough et al., 2008). End-to-end annealing of actin filaments is also well-established in vitro (Sept et al., 1999; Andrianantoandro et al., 2001; Popp et al., 2007) as well as in budding yeast (Okreglak and Drubin, 2010). Cellular factors such as cofilin and Aip1 may indeed allow filament annealing to the barbed end after severing while also restricting it from resuming elongation (Okada et al., 2002; Wioland et al., 2017). Additionally, an independent-particle Monte Carlo model based on actin SiMS data used to model FRAP of actin in lamellipodia (Smith et al., 2013), as well as actin monomer photoactivation experiments (Vitriol et al., 2015), provided better fits with local recycling of slowly diffusing actin back into the network throughout the lamellipodium. A large fraction of slowly diffusing oligomers were indeed observed in fragments of keratocyte lamellipodia (Raz-Ben Aroush et al., 2017); however, the model developed by these authors did not consider or require the presence of a local recycling mechanism.
Distributed turnover through severing and annealing may relate to another puzzle of lamellipodia structure revealed in electron micrographs (Svitkina et al., 1997): while a dense branched brushwork of filaments is observed near the leading edge (approximately within 1 µm), further away from the leading edge (approximately 3–4 µm away), filaments are longer and appear more linear. The mechanism required for this remodeling has yet to be determined.
To test the hypothesis of frequent severing and annealing in distributed turnover and structural remodeling of the actin network, we created a three-dimensional kinetic model of a steady-state lamellipodium based on the dendritic nucleation model. To develop a network with the observed ±35° filament orientation pattern, we systematically examined the self-organizing filament orientation pattern as a function of the relative network growth speed. Using parameter sets matching lamellipodia of the widely studied keratocyte or XTC cell types, we perform a search over parameters describing uniform severing along the actin filament and enhanced severing near barbed ends. The model considers the diffusion and annealing of severed oligomers. We find that a model without annealing cannot reproduce both the filamentous lamellipodial structure and actin SiMS data. With the addition of oligomer diffusion, annealing and enhanced severing near barbed ends, the structure, SiMS data, and overall increase in filament length with distance from the leading edge can be reproduced for optimized parameters. We support this mechanism by performing SiMS of Dylight-labeled actin on XTC cells indicating frequent disassembly of recently polymerized F-actin close to the leading edge. Our study thus supports that frequent severing and annealing is an important mechanism in cellular actin dynamics, motivating further experimental investigations.
Results
Stochastic simulation of dendritic network
The model shown in Figure 1 includes barbed end polymerization, pointed end depolymerization, capping, uncapping, branching near the leading edge, debranching, severing and annealing, without explicitly considering ATP hydrolysis or phosphate release (see Materials and methods and Table 1). We impose a constant network velocity with respect to the leading edge. We selected parameters corresponding to two frequently-studied cell systems, one for fast moving keratocyte cells and one for XTC or fibroblast cells. For keratocytes, corresponds to the rate of cell protrusion, since the actin network is almost stationary with respect to the substrate in experiments (Keren et al., 2008; Mueller et al., 2017; Yamashiro et al., 2018; Schaub et al., 2007). Lamellipodia of XTC cells are frequently studied in cells that do not crawl on the substrate, so provides the magnitude of the retrograde flow speed (Watanabe and Mitchison, 2002; Ryan et al., 2012). We considered uniform severing with rate constant per filament length and enhanced end severing with a rate per filament near the barbed end. If an oligomer is created from one of these severing events, we assume it can diffuse and anneal to a nearby filament end with rate constant .
Planar branching along lamellipodial plane sharpens the filament orientation pattern
Prior models of dendritic networks demonstrated how the ±35° orientation with respect to the direction of protrusion depends on the relationship between filament elongation velocity and relative extension rate . As shown in the results of the 2D model by Weichsel and Schwarz, 2010 in Figure 2A, for low , the dominant orientation pattern has filaments branching at . Filaments oriented at angles larger than , for which , lose contact with the membrane since they are not polymerizing quick enough to catch up. When this critical angle becomes smaller than 70°, the favored pattern is filaments with orientations centered at ±35°: the filament population around 35° can generate daughter branches at -35° and vice versa; thus the population sustains itself even as individual filaments get capped. The was found to resume when exceeds the critical angle for 35°; in this situation, only the 0° filaments polymerize quick enough to remain in the branching region. The above behavior can be quantified by the orientation order parameter, where a value 1 (or –1) indicates all filaments are in a –70°/0°/70° (or ±35°) pattern.
By contrast to the above 2D results, dendritic network models formulated in 3D have provided apparently contradictory results (Holz and Vavylonis, 2018). Atilgan et al., 2005 reported that obtaining the ±35° pattern requires restricting branching to occur primarily along the lamellipodium plane, which they attributed to structural constraints of the branching machinery at the leading edge. Schreiber et al., 2010 and Hu and Papoian, 2010 however did observe the ±35° pattern in 3D simulations, but the role of in determining the pattern was not examined.
As we are interested in structural aspects of lamellipodia, we performed systematic simulations in both 2D and 3D to examine the filament orientation pattern as function of . In the simulations of Figure 2 we consider the region close to the leading edge where severing, annealing and debranching does not influence the resulting structure, and varied the polymerization rate while keeping constant. We also kept the capping and branching rates and fixed as the orientation pattern is robust with respect to their values (Weichsel and Schwarz, 2010).
We reproduce the results of the Weichsel and Schwarz, 2010 simulations as a function of relative network growth speed, by imposing a tight planar branch restriction along the lamellipodial plane (Figure 2A). Note that the transition among different orientation patterns is not abrupt at the critical angles, as a result of allowing fluctuations in branching angle and a finite size of the branching region. This is the reason for the largest difference occurring at high : because only filaments at small angles can keep up with the leading edge under these conditions, our simulation evolves to a narrow comet-like branching structure where all branching is concentrated; this allows branching among ±35° and other orientations (e.g. ) to be maintained through double branching before filaments exit the branching region.
There is little difference in the orientation order parameter between simulations that allow backward branching (angles greater than 80°, chosen to include angles that would be oriented away from the leading edge considering the variation of branching angle for both orientation patterns) and branching limited toward the leading edge (Figure 2A). Filaments that branch backward exit the branching region quicker than filaments that branched forward, decreasing the likelihood of branches off of backward filaments. Even though backward-facing filaments contribute to angles larger than 80°, they do not influence the order parameter that does not measure them. Since backward branching does not affect the filament orientation pattern, and since backward-facing filaments are not seen in electron microscopy images (Vinzenz et al., 2012; Mueller et al., 2017), backward branching was not allowed in the rest of our simulations.
Next, we allowed branching to occur uniformly in 3D (i.e. with equal probability along any azimuthal angle with respect to the axis of the parent filament) and studied the orientation order parameter as function of (Figure 2B). The order parameter was measured using the angles of filaments projected along the lamellipodial plane. Uniform branching led to low and weakly varying order parameter, unlike the sharp orientation pattern with distinct transitions observed in 2D. By comparison, a quasi-2D simulation with branching allowed to occur within 10° off the lamellipodial plane (close to the maximum off-plane angles observed in electron microscopy [Vinzenz et al., 2012] [A. Narita, personal communication, March 2018]) restored the behavior observed in 2D.
In the simulations described so far, filament elongation was assumed to stop when the polymerizing barbed end reached the top or bottom plane. To investigate the influence of filament bending along the membrane, we also performed simulations with filament ‘kinking’, in which filaments were allowed to continue their elongation parallel to the plane representing the top or bottom cell membrane. Allowing kinking lead to a high concentration of filaments along the top and bottom plane of the lamellipodium but did not however significantly influence the orientation order parameter (Figure 2B).
To further visualize the network structure in simulations, Figure 2C and D (and Figure 2—figure supplement 1 for the case with kinking) show snapshots and filament orientation patterns at low and intermediate values of relative network growth speeds. The simulations with quasi-2D branching show clear –70°/0°/70° and ±35° respective orientation patterns. Interestingly, even though no prominent features are observed in simulations with uniform 3D branching at low relative network speeds (Figure 2C), intermediate relative network speeds do show features at ±35° (Figure 2D). The latter histogram is not very different from experimentally-measured distributions (Maly and Borisy, 2002; Vinzenz et al., 2012; Koseki et al., 2019).
We thus conclude that the ±35° pattern does occur within a broad range of relative network growth speeds with uniform 3D branching, however the peaks at ±35° are not very pronounced. Indeed, the parameters used by Atilgan et al., 2005 corresponded to , a parameter set that mostly lies outside the ±35° region; this is likely the reason why the ±35° was not observed in this study. Our results also suggest why Schreiber et al., 2010 who used and Hu and Papoian, 2010 who had , did observe a ±35° with uniform 3D branching.
Considering the simulation results as well as experimental evidence in electron tomograms for filaments oriented primarily along the lamellipodial plane (Vinzenz et al., 2012), for the rest of the simulations we proceed with the quasi-2D case where filament branching occurs within 10° of the lamellipodial plane and relative network growth speeds result in a ±35° orientation pattern.
Estimated parameters
For both cases of keratocytes and XTC cells, which correspond to different values for , we estimated the rates of polymerizing barbed end elongation , branching , and capping , that are needed for a dendritic network with the anticipated concentration, branch length, and ±35° filament orientation pattern (Table 1). We also assume that barbed ends can uncap with rates comparable to those in SiMS lifetime measurements of capping protein (Miyoshi et al., 2006). The availability of such uncapped barbed ends for annealing is an important assumption of this work. We further use debranching rates by considering the lifetime of Arp2/3 complex components in SiMS (Miyoshi et al., 2006) and measurements of Arp2/3 complex profiles in lamellipodia (Iwasa and Mullins, 2007; Lai et al., 2008; Ryan et al., 2012).
We also assumed that uncapped barbed ends away from the leading edge do not elongate or shrink and that free pointed ends depolymerize with a rate /s. The results we present below are robust with respect to small changes of these parameters, as long as the overall filament disassembly rate away from the leading edge is not reaching values comparable to . Maintaining a wide lamellipodium in the latter case would require a global treadmilling mechanism, which would contradict the evidence for distributed turnover. The assumption of slow barbed end dynamics away from the leading edge is consistent with the slow intensity increase in the back of lamellipodium after FRAP of actin (Smith et al., 2013) or after photoactivation of actin at the cell middle (Lai et al., 2008; Vitriol et al., 2015), as well as evidence that cofilin and twinfilin promote both barbed and pointed end depolymerization (Wioland et al., 2017; Johnston et al., 2015; Hakala et al., 2021; Shekhar et al., 2021).
A steady state with a finite lamellipodium depth is reached in the simulation whenever the net rate of depolymerization balances the net rate of polymerization at the leading edge. For example, in the case without enhanced end severing, annealing, or oligomer dissociation, the depth of the lamellipodium is determined by the rates of uniform severing, , and (Figure 1—figure supplement 1): in this case, the fast growth of barbed ends at the leading edge is balanced by the slower depolymerization of a larger number of pointed ends created by severing.
Given the parameters in Table 1, this leaves three main unknown parameters related to oligomer dissociation: , the rate of uniform severing along each filament; , the longest length of a diffusing oligomer; , the enhanced severing rate near the barbed end (within of the end). We treated these three as fitting parameters and considered separately the cases in the presence or absence of annealing. In the absence of annealing, dissociating oligomers do not reincorporate into the network and are thus discarded from the simulation (corresponding to eventual disassembly into monomers, a process that we did not simulate).
Model without annealing cannot reproduce both actin SiMS data and lamellipodial structure
We conducted a parameter search over the maximum oligomer size as well as end and uniform severing rates , for keratocyte- and XTC-like parameters without annealing (Figure 3—figure supplement 1). We classified each parameter set in terms of how well it described the F-actin structure and concentration (including barbed end and branch concentration profiles), actin SiMS data (speckle lifetime distribution, appearance and disappearance profile), and if there was an increase in length between the filaments located in 0–1 µm and 3–4 µm region (see Materials and methods). Within a range of uniform and end severing rates, and short maximum oligomer lengths, the network structure and concentration of our simulation was close to that expected for keratocyte and XTC lamellipodium (Figure 3—figure supplement 1). However, as expected, we did not find any parameter sets that resulted in a length increase since a mechanism for an increase in length and remodeling is not included. We also did not find any parameter set with a sharp peak at short actin speckle lifetimes, though in some cases the actin speckle lifetime distribution has a peak at relatively short times (Figure 3—figure supplement 1).
To further demonstrate that the model without annealing cannot fit the experimental data, Figure 3 and Figure 3—figure supplement 2 contain examples of results for keratocyte and XTC parameters. In this and subsequent figures we color parameter sets black without end severing and red that include enhanced end severing; these curves correspond to the parameters of the scan that are highlighted with a thick frame of same color in Figure 3—figure supplement 1. For the case with moderate uniform severing and no end severing, both the XTC and keratocyte uniform severing cases have a concentration profile that is comparable to lamellipodium of their respective parameter set (Figure 3) as well as barbed end and branch distributions (Figure 3—figure supplement 2). However, the actin speckle lifetimes do not peak at short lifetimes. With the addition of enhanced end severing, as well as increase of uniform severing, a peak at short lifetimes is observed that is closer to the experimental SiMS curves; however, the lamellipodium becomes too narrow, there is a shortage of long speckle lifetimes compared to experiment, and the location of speckle appearances is restricted close the leading edge (Figure 3 and Figure 3—figure supplement 2 A,F). In conclusion, the model without annealing cannot reproduce the distributed turnover, structure of actin network and increase in filament length.
Model with oligomer annealing can reproduce lamellipodium structure, actin speckle dynamics and increase in filament length away from leading edge
Next, we examined if inclusion of oligomer annealing might be able to provide an adequate fit to the structure, speckle, and length increase criteria. We performed another parameter search over the maximum oligomer length, end and uniform severing rates (, , ) for keratocyte and XTC parameter sets. For these simulations, we used an annealing rate constant measured in crowded surfaces in vitro (Popp et al., 2007). Similar to the parameter scan without annealing, a match to F-actin structure and concentration was obtained when uniform and end severing rates were within a certain range (Figure 4—figure supplement 1 and Figure 5—figure supplement 1). The addition of annealing improved the agreement between the speckle lifetimes, appearance and disappearance location profiles compared to SiMS data, resulting in parameter sets that agree with both structure and speckle dynamics (Figure 4—figure supplement 1 and Figure 5—figure supplement 1). We also see that with high enough severing and maximum oligomer lengths, the filament length increases in a region away from the leading edge (Figure 4—figure supplement 1 and Figure 5—figure supplement 1). Parameter sets where all three fitting criteria are satisfied exist for both keratocyte and XTC cases, and all these triple matches have a finite enhanced end severing rate.
Detailed results from our parameter scan for keratocytes are shown in Figure 4. This figure shows a case without enhanced end severing (black curves, double match in structure and speckles for keratocyte parameters) and with enhanced end severing (red, triple match in structure, speckles and filament length increase). Both parameter sets can reproduce the actin speckle dynamics as seen in the speckle lifetime, appearance and disappearance location distributions as well as the the structure as seen in the F-actin, barbed end, and branch concentration profiles (Figure 4A–F). However, only the case with enhanced end severing results in a simultaneous increase in length away from the leading edge (Figure 4G). The increase brings the filament length close to the estimated average of 800 nm in keratocytes (Schaub et al., 2007). The concentration profile of Arp2/3 complex branches (Figure 4F), snapshots (Figure 4H), and Video 1 of the optimized simulations with enhanced end severing clearly show more branches and shorter filaments in the region near the leading edge than away from the leading edge, similar to the electron micrographs of Svitkina et al., 1997. The ±35° filament orientation pattern is preserved throughout the lamellipodium in the presence of severing and annealing. Additional quantification in Figure 4—figure supplement 2 shows the steady state spatial distribution of capped and uncapped barbed ends, the approximate uniform distribution of oligomer sizes between 0 and , and the profile of F-actin according the mechanism of assembly (polymerization as monomer at the leading edge versus annealing). Most of the F-actin at the back of the lamellipodium has undergone severing and annealing (Figure 4—figure supplement 2D.F), similar to an earlier particle model for distributed turnover (Smith et al., 2013).
For XTC cells, enhanced end severing is needed to match the actin speckle lifetime distribution, which further contributes to filament length increase. Results from our parameter scan for XTC cells in Figure 5 show a case without enhanced end severing (black curves, match in structure only) and a case with enhanced end severing (red, triple match in structure, speckles and filament length increase). The case with enhanced end severing provides a good overall fit to SiMS and structure (Figure 5A–G), demonstrating increase in length away from the leading edge (Figure 5G and H and Video 2). Figure 5B and F show that the model reproduces the narrower distribution of Arp2/3 complex as compared to F-actin in XTC cells (Ryan et al., 2012). Additional quantification in Figure 5—figure supplement 2 shows profiles of capped and uncapped end, oligomer size distribution, and origin of F-actin.
When performing the F-actin concentration and structure match for the XTC parameter sets, we used numbers in between those of electron tomograms in fibroblasts by Vinzenz et al., 2012 and the estimated F-actin concentration of 1000 µM for XTC cells (Watanabe, 2010). The branching rate in the simulations of Figure 5 corresponds to 0.05 µM s-1, which is about half of the Arp2/3 complex nucleation rate of 0.11 µM s-1 estimated by SiMS (Watanabe, 2010; Miyoshi et al., 2006). We checked that simulations with doubled the branching rate still provide a good fit to actin SiMS data as well as a length increase away from the leading edge, with F-actin concentration at the leading edge that was around 1,100 µM (Figure 5—figure supplement 3). We also tested that excluding end severing of polymerizing ends in Figures 4 and 5 did not modify our results for the optimized parameter sets (Figure 5—figure supplement 4).
The optimized parameter sets with speckle dynamics similar to SiMS experiments of Figures 4 and 5 also matched another observation from SiMS (Watanabe and Mitchison, 2002): the distribution of actin speckle lifetimes was weakly dependent on location of appearance with respect to the leading edge (Figure 4—figure supplement 3, Figure 5—figure supplement 5).
The simulations of Figures 4 and 5 implement a mechanism of local oligomer rebinding, which is needed to match actin FRAP and photoactivation data (Smith and Liu, 2013; Vitriol et al., 2015): as shown in Figure 4—figure supplement 3B and Figure 5—figure supplement 5B the distance travelled by oligomers before annealing is in the sub-µm range. As a further check of consistency of our simulations with SiMS results, simulated actin SiMS for parameters with enhanced end severing of Figure 4 (keratocytes) and Figure 5 (XTC) do resemble experimental images from Yamashiro et al., 2014 and Watanabe and Mitchison, 2002 (Videos 3 and 4). In Videos 3 and 4, when an appearing and disappearing speckle are near one another, this is typically a reannealing event occurring quickly over short distances. We confirmed that events that might be limited by spatial and temporal resolution or interpreted as blinking in SiMS experiments correspond to a very small fraction of speckle appearances in the simulation.
Finally, we note that even though the results of this section were obtained for a specific value of annealing rate constant and oligomer diffusion coefficient, they remain valid as long as the annealing of oligomers occurs over a sufficiently short distance. We find that this is the case even for annealing rate constants that are lower by nearly two orders of magnitude compared to the values of Table 1 (Figure 4—figure supplement 3C, Figure 5—figure supplement 5C).
Alternative mechanism with frequent barbed end destabilization
An alternative mechanism to explain the short actin SiMS lifetimes (different to enhanced severing and annealing) is barbed end catastrophic disassembly by factors such as twinfillin (Wioland et al., 2017; Johnston et al., 2015; Hakala et al., 2021; Shekhar et al., 2021) or cooperative strand separation in the presence of cofilin, coronin and Aip1, depending on the cofactor concentrations (Kueh et al., 2008; Jansen et al., 2015; Tang et al., 2020).
To test such a mechanism, we implemented a model with stochastic transitions to rapid barbed end depolymerization (see Materials and methods, Video 5). To maintain the F-actin loss by such a disassembly process, a process of rapid regrowth must also be included. Using the insight gained from our parameter searches with severing and annealing, we can show that such a model can be tuned to come close to matching our three main experimental test criteria: agreement with SiMS data, an increase in filament length with distance from the leading edge, and broad F-actin concentration profile (Figure 4—figure supplement 4).
While we cannot fully exclude such a dynamic-instability-like mechanism, we note that: (1) it would require additional controls or homeostatic mechanisms to balance disassembly and reassembly away from the leading edge (while maintenance of F-actin mass is ensured by a severing and annealing mechanism), (2) ATP hydrolysis associated with ATP-actin monomer polymerization would be energetically more costly, and (3) reassembly of rapidly diffusing monomers away from the leading edge to recover bursting would also be less consistent with FRAP or photoactivation experiments (Smith et al., 2013; Vitriol et al., 2015).
Short actin speckle lifetimes provide evidence for rapid disassembly near barbed ends
The actin speckle lifetime distribution of the models with severing and annealing of Figures 4–5 characteristically peaks at . This is a general feature of a barbed-end disassembly mechanism where newly polymerized monomers are the ones that disassemble at higher rates, being closer to the barbed end. The experimental SiMS lifetimes in Figures 3—5 used a temporal resolution of 1 s (Watanabe and Mitchison, 2002; Yamashiro et al., 2018). To probe the kinetics at even shorter lifetimes, SiMS of Dylight actin was repeated on XTC cells using the methods of Yamashiro et al., 2014, with a temporal resolution of 0.1 s. We measured actin speckles appearing near the leading edge, which should represent polymerization of actin monomers rather than oligomer annealing. Approximately 16% and 30% of speckles disappeared within within 0.5 and 1 s, respectively (Figure 6A). This large amount of short speckle lifetimes strongly supports frequent disassembly of newly polymerized F-actin near the barbed end.
To compare to simulations, Figure 6B contains speckle lifetimes within 1 µm of the leading edge for the XTC parameter set with enhanced end severing and annealing of Figure 5 (the uniform severing case of Figure 5 did not provide as good a match to SiMS measurements and is not shown). The speckle distribution extended to short lifetimes, similar to the experiments of Figure 6A, with the percentage of lifetimes within 0.5 s and 1 s being 3% and 7 %, respectively. Doubling the enhanced end severing rate while also choosing parameters that reproduced speckle dynamics, lamellipodial structure and an overall increase in length ( /sub/s; sub, XTC parameters, Table 1), doubled the percentage of speckle short lifetimes within 0.5 s and 1 s to 7% and 14%, respectively, closer to the experimental percentages of Figure 6A. We conclude the short lifetimes of Figure 6A are within the range of what is expected by an enhanced end severing mechanism, which may even occur at a rate 2–4 times faster than the estimate of Figure 5.
Kueh et al., 2008 observed that actin disassembly in lamellipodia was inhibited in the presence of cytochalasin D (CD), a barbed end capper that also inhibits binding of cofilin to F-actin at high concentrations (Shoji et al., 2012). We tested the effect of CD treatment on XTC cells using SiMS (Figure 6C). We confirmed CD’s inhibitory effect on filament disassembly by observing a larger fraction of actin speckles that survive after CD treatment. These results provide further support for rapid disassembly of newly polymerized actin at barbed ends (either through ‘bursting’ as proposed in Kueh et al., 2008 or enhanced severing).
To check if tropomyosin (TPM) might be involved in the reorganization of lamellipodial networks by debranching and stabilizing progressively the longest filaments, we observed SiMS of TPM-EGFP, as described previously in Higashida et al., 2008, in lamellipodia of XTC cells ~30 min after cell spreading. In this early phase, the density of TPM SiMS was substantially lower in lamellipodia than in lamella (Figure 6—figure supplement 1A, Video 6). Within lamellipodia, the density of TPM SiMS becomes gradually higher toward the base of lamellipodia. TPM SiMS are scarce near the leading edge, suggesting that TPM may not have a significant effect on actin turnover near the leading edge.
We also compared the dissociation rate of tropomyosin between lamellipodia and lamella. In lamellipodia of XTC cells 30 min after cell spreading, TPM SiMS dissociated much faster in the lamellipodia than in the lamella (Figure 6—figure supplement 1B). TPM-EGFP was localized to actin stress fibers after cells were grown for 24 hr in the presence of serum (Figure 6—figure supplement 1C). The localization and kinetics of TPM SiMS are in agreement with slow assembly kinetics of TPM along the actin filament in vitro (Schmidt et al., 2015). We thus conclude that at least under our observation conditions, it is unlikely that TPM stabilizes actin filaments near the tip of XTC lamellipodia.
Discussion
We showed that the hypothesis of frequent severing and annealing (Miyoshi and Watanabe, 2013) provides a mechanism for distributed turnover and structural remodeling of the actin network. Using simulations based on the dendritic nucleation model, under conditions that allow self-organization into a ±35° filament orientation pattern, we determined values for uniform and enhanced end severing rates that can simultaneously account for a diverse set of experimental data: (1) actin SiMS measurements (Watanabe and Mitchison, 2002; Smith et al., 2013; Yamashiro et al., 2018), (2) actin photobleaching and photoactivation experiments (Theriot and Mitchison, 1991; Lai et al., 2008; Smith et al., 2013; Vitriol et al., 2015) (since our model incorporates distributed turnover Smith et al., 2013), (3) the presence of uncapped barbed ends through the lamellipodium (Miyoshi et al., 2006; Raz-Ben Aroush et al., 2017), and (4) the change of network structure of the lamellipodium as a function of distance to leading edge (Svitkina and Borisy, 1999). While there is no direct measurement of actin oligomer lengths in cells, we note that our prediction of average length of 40–80 subunits (0.11–0.21 µm, Figure 4—figure supplement 2 and Figure 5—figure supplement 2) is smaller to 0.25 µm, the length of diffusing actin filaments estimated by FCS in the cortex of Hela cells (Gowrishankar et al., 2012) but larger than 13 subunits estimated by FRAP in keratocyte fragments (Raz-Ben Aroush et al., 2017).
Additionally, our simulations are consistent with SiMS of Arp2/3 complex (Millius et al., 2012) and capping protein (Miyoshi et al., 2006; Smith et al., 2011) in XTC cells. Agreement with SiMS of bound Arp2/3 complex lifetimes occurs because the speckle lifetimes correspond to a narrow distribution of branches as compared to the lamellipodium width (McMillen and Vavylonis, 2016), as in Figure 5F. The model is also in agreement with SiMS measurements of capping protein lifetimes in XTC cells (Miyoshi et al., 2006), which were used as input to the uncapping rate constant. This uncapping is important in the model, to allow uncapped barbed ends for annealing. The presence of slowly diffusing oligomers assumed in the model could also be consistent with the presence of slowly diffusing capping proteins in XTC cell lamellipodia (Smith et al., 2011). For the parameters of Figures 4–5, release of capped oligomers through end severing would not contribute significantly to capping protein SiMS lifetimes, since the corresponding rate is 4–8 times slower than uncapping.
It is also interesting to compare our model to SiMS measurement of Aip1 in XTC cells (Tsuji et al., 2009). Assuming that appearance of Aip1 speckles corresponds to filament disruption, the frequency of Aip1-associated filament disruption was estimated to be 1.8 µM s-1 (Tsuji et al., 2009). We calculated the simulated overall effective severing rate by counting the total number of uniform and enhanced end severing within 20 s in steady state within 5 µm from the leading edge (average of 5 simulations). For the enhanced end severing XTC case in Figure 5 and Figure 5—figure supplement 3, the effective severing rate was 0.17 µM s-1 and 0.32 µM s-1, respectively. These values are a few times smaller, yet not too far from the experimental estimate. This difference could indicate an even higher enhanced end severing rate in cells as compared to Figure 5, as also suggested by the comparison of experimental and simulated short actin lifetimes in Figure 6A and B.
Prior measurements of the intensity pattern of phalloidin-stained actin filaments in keratocytes treated with low doses of CD indicates shorter filaments compared to the control case (Schaub et al., 2007). This filament shortening could be due to capping by CD, as suggested (Schaub et al., 2007), however it may also be related to reduced structural remodeling of the lamellipodium through severing/destabilization near barbed ends.
Actin filament annealing, a basic assumption of our model, has been established in vitro (Sept et al., 1999; Andrianantoandro et al., 2001), including on a crowded surface which is similar to lamellipodial conditions (Popp et al., 2007). Our results are robust with respect to the annealing rate constant and oligomer diffusion coefficient, which is why we did not elaborate on the precise length dependence of oligomer diffusion and annealing rates. Specifically, the results are valid for a range of diffusion coefficients and annealing rate constants, as long as annealing is not dominated by reannealing to the same filament, diffusion is fast enough such that severing contributes to speckle disappearance in SiMS, and annealing does not occur further than approximately 1µm from the severing location. For example, the simulations results of Figures 4–5 are nearly identical when reducing by 10 times compared to the value used from Popp et al., 2007. This is because both annealing rates are in a range that allows the oligomers to diffuse away from the same severed filament and large enough for annealing to remodel the network. Future work is needed however to further investigate how diffusion and annealing of filament segments occurs within the dense lamellipodial actin mesh.
Our simulations showed how severing and depolymerization regulate the length of the lamellipodial dendritic network (Figure 2—figure supplement 1), similar to earlier models that have been formulated at various levels of description (Edelstein-Keshet and Ermentrout, 2001; Mogilner and Edelstein-Keshet, 2002; Carlsson, 2007; Michalski and Carlsson, 2010; Ditlev et al., 2009; Berro et al., 2010; Michalski and Carlsson, 2011; Lewalle et al., 2014; Manhart et al., 2019), including use of explicit dendritic network (Schreiber et al., 2010). We also included annealing, a process previously studied using lattice models by Carlsson and Michalski (Carlsson, 2007; Michalski and Carlsson, 2010; Michalski and Carlsson, 2011). In their models, annealing was implemented as reappearance of lattice links and lead to a wider lamellipodium, similar to our findings. Annealing was also included in the partial differential equations model of Ditlev et al., 2009 but the implications of this mechanism were not explicitly explored.
Closer to our work is the model by Huber et al., 2008 who created a detailed 2D kinetic Monte Carlo model of the keratocyte lamellipodium to calculate the filament lengths and concentration profiles of actin and associated proteins. The model included diffusion of free actin monomers, filament nucleation along predefined ±35° orientations to represent branching, permanent barbed end capping, binding of ADF/cofilin and tropomyosin to filaments, as well as ATP hydrolysis and Pi release. Annealing among filaments was also included, without explicit modeling of severed filament diffusion. With this model, two distinct network regions formed (termed the lamellipodium and lamella in that paper): a region with short filaments close to the leading edge, followed by a region with longer filaments starting at about 2 µm further away.
The length increase with distance from the leading edge in the model of Huber et al., 2008 occurred via two different mechanisms: (i) filament annealing, or (ii) polymerization of barbed ends away from the leading edge, created by severing; these ends were assumed to polymerize faster than barbed ends at the leading edge, as a result of the higher G-actin concentration away from the leading edge. The annealing mechanism (i) of Huber et al., 2008 is different to what we described in our work: in their study it involved the joining together of any pair of filaments at the same distance from the leading edge, without accounting of the dendritic network topology. By contrast, we assumed annealing involves diffusing oligomers. Their work also involved annealing rate constants that were one to two orders of magnitude smaller than the values in this work, and, as we understand, it was assumed that annealing occurred even with capped barbed ends. An annealing mechanism as in Huber et al., 2008 would not contribute to actin speckle appearances and disappearances (that were not quantified in their paper) and is likely inconsistent with SiMS data. However we note that the Huber et al., 2008 model did lead to a peak in actin disassembly at 1–2 µm away from the leading edge, similar to Figure 4D, as well as experimental results by qFSM in other cell types (Ponti et al., 2004). It’s also unclear if fast polymerization away from the leading edge (mechanism (ii) in Huber et al., 2008) is consistent with data using FRAP (Lai et al., 2008; Smith et al., 2013) and photoactivation (Vitriol et al., 2015) of actin in lamellipodia of other cell types. These experiments argue against incorporation of fast-diffusing actin at the back of the lamellipodium (Smith et al., 2013). We also note that when we included fast polymerization of barbed ends away from the leading edge in our simulations, we typically obtained a high F-actin concentration peak away from the leading edge, unlike in the simulations of Huber et al., 2008 or in prior experiments.
In a study of actin dynamics that combined experiment and modeling, Raz-Ben Aroush et al., 2017 found evidence for the presence of a large pool of short actin oligomers in keratocyte lamellipodial fragments. Using FRAP on small regions, fluorescence correlation spectroscopy, and phalloidin labeling, the oligomer diffusion coefficient was estimated to be about 5 µm2s-1 for oligomers with an average length of 13 subunits. Raz-Ben Aroush et al., 2017 also report that two-thirds of actin within these fragments are diffuse, with oligomers composed a sizable fraction of this pool. A partial differential equation 1D model (that included polymerizable and non-polymerizable monomers, oligomers and F-actin) provided agreement with the data, assuming actin disassembly into oligomers throughout the lamellipodia and a broad distribution of polymerizing barbed ends. The finding of oligomers, as well as the proposed distributed F-actin turnover, is consistent with main assumptions of our work. However we note that the absence of local reassembly in the mechanism proposed by Raz-Ben Aroush et al., 2017 may not easily explain experimental observations of actin FRAP or photoactivation over large regions of other cell types (Smith et al., 2013; Vitriol et al., 2015; Yamashiro et al., 2018), or the filament length increase across the lamellipodium. We also note that other studies (Smith et al., 2013; Kiuchi et al., 2011) have suggested much smaller concentrations of diffuse actin in lamellipodia compared to Raz-Ben Aroush et al., 2017, although the situation could be different in the faster keratocytes (Yamashiro et al., 2018) and their fragments.
The mechanism of severing and annealing modeled in this work could represent a general feature of actin dendritic networks, including yeast cells where short actin speckle lifetimes have been observed in actin patches of fission yeast (Lacy et al., 2019). It might provide an energetically efficient mechanism for network remodeling matching different mechanical requirements: close to the leading edge, short branched networks would provide rigidity to compressive stresses (resulting from actin polymerization against the membrane) while longer filaments at the back might be better suited for extensional stresses through myosin motors. Future work is however needed to clarify the biochemical basis of the proposed kinetics, taking into account the energetic requirements associated with ATP hydrolysis and Pi release along actin filaments, as well as mechanics and kinetics of actin filament side-binding proteins such as cofilin, GMF, and tropomyosin.
Materials and methods
Simulation
Request a detailed protocolWe developed a three-dimensional stochastic simulation of the actin network within the lamellipodium (Figure 1). Actin filaments are represented as straight lines without excluded volume, within the simulation box. We work in the coordinate system where the leading edge is at rest. For the purposes of this work, we do not consider explicitly the effects of ATP hydrolysis, diffusion of the actin monomer pool, or excluded volume interactions among filaments. The leading edge at y = 0, as well as the lamellipodium top and bottom at and are hard boundaries. There is no boundary at , allowing the mechanisms and rate constants of the system to determine the length of the lamellipodium. Periodic boundary conditions with length at least 1 µm are applied in the x-direction. A constant relative velocity, , is imposed between the network and the leading edge. The network is initialized by filament seeds of 5 subunits in length at random orientations near the leading edge. The probability of a reaction event (polymerization, depolymerization, capping, uncapping, branching, severing, and annealing) is calculated using the corresponding rates or rate constants as described below and in Table 1. The time step was s.
Mechanisms
Polymerization, depolymerization, capping, and uncapping
Request a detailed protocolFree barbed ends, which are created by branching at the leading edge polymerize with rate . Polymerization (elongation) is simulated as stochastic increase of filament length by 2.7 nm, corresponding to addition of one monomer, when this is allowed by the distance to the hard boundaries at , , or . We examined two different scenarios when an elongating filament reaches or : (i) the filaments stop polymerizing, or (ii) undergo ‘kinking’, namely they continue elongation parallel to or , at the same angle along the plane. The latter is implemented to mimic bending of filaments when they come in contact with the membrane. Capping of free barbed ends occurs at a rate , which stops polymerization and does not allow annealing of oligomers. Capped barbed ends become free with uncapping rate, . Free barbed ends away from the leading edge are expected to undergo different polymerization kinetics compared to barbed end at the leading edge, where membrane bound proteins such as Ena/VASP catalyze fast elongation. Recent evidence suggests that ADF/cofilin and twinfillin may assist in the depolymerization of barbed ends away from the leading edge at a slow rate (Wioland et al., 2017; Johnston et al., 2015; Hakala et al., 2021; Shekhar et al., 2021). For simplicity, and accounting for these recent observations, barbed ends formed by uncapping neither polymerize or depolymerize, still allowing however annealing of oligomers. Free pointed ends (created by severing or debranching) depolymerize with rate , implemented by stochastic length decrease by 2.7 nm corresponding to one monomer. Depolymerization stops when the filament completely depolymerizes or meets the qualifications for an oligomer. These mechanisms are depicted in Figure 1B.
Branching and debranching
Request a detailed protocolFilament branches are nucleated at a total rate and placed randomly along parent filaments, in proportion to their segment length within the branching region, a 27 nm region near the leading edge (Figure 1A), approximately the size of an Arp2/3 complex associated to proteins on the cell membrane (Volkmann et al., 2001). New branches form from an existing parent filament at an angle chosen from a Gaussian distribution centered at 70° with a standard deviation of 5° (Weichsel and Schwarz, 2010; Gong et al., 2017). The azimuthal angle of the branch around the axis of the parent filament is either picked from a uniform distribution (uniform 3D branching) or else uniformly but with the additional condition that the angle between the branch and the lamellipodium plane is smaller than a threshold value (typically 10°, case of branching along lamellipodium plane). Unless otherwise indicated, branching at angles larger than 80° with respect to the axis of protrusion (termed “backward” branching) is not allowed; in cases where an invalid orientation is selected, a new branch location and orientation is tried. Debranching occurs at rate , which we assume results in release of the Arp2/3 complex, leading to a free pointed end for the debranched filament. In simulations with , debranching was assumed to occur for branches that become five monomers or smaller and do not contain branches of their own. Debranching also occurs when the pointed end of a parent filament depolymerizes past a branch.
Severing
Request a detailed protocolWe considered uniform severing with constant rate per filament length, and enhanced severing near the barbed end with rate per filament length (proposed as possible explanation for the short actin lifetimes observed by SiMS Miyoshi et al., 2006). Enhanced end severing occurred between the barbed end and , the longest length for an oligomer. In a severing event, a location on the filament is chosen to split the filament into two, creating a new depolymerizing pointed end and a free barbed end. An oligomer, assumed to be diffusing as discussed below, is created for any filament segment that is shorter than and is not a branch or does not contain a branch. Oligomers formed at the barbed end retain the capped or polymerization state of the original filament. To allow a network to form even at high severing rates, we assumed that no severing occurs within 0.1 µm of the leading edge and did not apply enhanced end severing to barbed ends polymerizing at the leading edge (unless otherwise indicated). The latter assumption does not have a significant effect on our final results (Figure 5—figure supplement 4).
Oligomers and annealing
Request a detailed protocolSince oligomers are short and typically anneal after a short time interval, we assumed that they do not branch, sever, or depolymerize, although we implemented capping and uncapping with the same rates as filaments. We assume that the probability of finding the end of a non-annealed oligomer displaced by distance and along the and directions (in the reference frame where the leading edge is at rest), after time , with respect to the location and time of its creation by severing, is given by 2D diffusion: , where is the oligomer diffusion coefficient. This expression neglects the small effect of advection by cytoplasmic fluid flow as well as the boundary condition at the leading edge.
Oligomer annealing to filament ends is calculated as a bimolecular reaction with rate constant . This is implemented by scanning through all pairs of oligomers and available pointed and barbed ends, converting to a local oligomer concentration by assuming a uniform probability along the thin direction, and using , as the distances between the end of the oligomer and the end on the filament that could anneal to one another. The smallest distance in the x-direction is used according to the periodic boundary conditions. If the annealing event is accepted, the length of the filament increases by the size of the oligomer. If the oligomer anneals to a barbed end, the barbed end state of the oligomer is transferred to the filament. If an oligomer did not anneal within 20 s of its creation, it was removed since in this time it likely disassembled. Removal was unlikely to occur at the chosen annealing rate constant since most oligomers annealed within a shorter time.
We checked the time step used was sufficiently small: µM-1s-1 resulted in less than 40% of oligomers annealing per time step, when using the reference parameter values in Table 1. We also checked that the median of the time to anneal, , and the median distance between the severing and annealing events, , were related as expected from a theoretical approximation of these quantities assuming all barbed ends are free: with , where is the average concentration of barbed ends in the body of the simulated lamellipodium (Figure 4—figure supplement 3, Figure 5—figure supplement 5).
Alternative model with frequent transitions to rapid barbed end depolymerization
Request a detailed protocolAs an alternate mechanism to explain the short actin SiMS lifetimes, we implemented a model with frequent transitions of barbed ends to a state of rapid depolymerization with rate . A large value of can mimic barbed end catastrophic disassembly (by factors such as twinfillin) or massive fragmentation (by factors such as Aip1). In the absence of annealing, this disassembly process must be matched by rapid actin repolymerization at rate throughout the lamellipodium, to maintain the F-actin concentration in the lamellipodium.
We assume the same mechanism of polymerization, capping, uncapping, branching, and uniform severing as described in Mechanisms above, without enhanced severing near the barbed end or annealing. A new feature is the assumption that uncapping leads to either fast barbed end depolymerization or fast repolymerization, with probability and , respectively. We thus use the uncapping rate as a parameter to control the rate of catastrophic disassembly and the capping rate as a rescue rate. To decouple barbed end catastrophic disassembly from severing, we now assumed that new barbed ends that form after severing are capped.
The alternate model introduces several new parameters, which are not known experimentally. However, its viability can be examined using a parameter set tuned to give similar effects to the optimized model with enhanced severing near barbed ends and annealing (Figures 4 and 5). Focusing on keratocytes, we keep the same as in Table 1 and uniform through the lamellipodium. We set such that the typical filament length lost in a catastrophic event, , is close to the typical length of the optimized enhanced severing model, in Figure 4. For simplicity, we set and such that repolymerization balances catastrophic disassembly.
Simulated SiMS
Request a detailed protocolTo compare our simulations to the actin speckle appearance and disappearance in SiMS, we tagged and tracked 1% of actin monomers that add to the network through polymerization at the leading edge. Each such polymerization corresponds to a speckle appearance. Disappearance events occur when the tagged monomers depolymerize off a pointed end or when they become part of an oligomer. Annealing of an oligomer carrying a tagged monomer is an appearance. The speckle lifetime is the time between appearance and disappearance events.
Orientation order parameter
Request a detailed protocolWe counted the number of filaments having a segment within the first micrometer of the leading edge and defined the order parameter similar to Weichsel and Schwarz, 2010 and Mueller et al., 2017: where and are the number of filaments oriented between -20° to 20° and 60° to 80° in either direction, or between 25° - 45° in either direction, respectively. The angle is measured with respect to the axis of network growth. For 3D simulations, we used the angles of the filaments projected along the lamellipodial plane.
SiMS imaging experiments
Request a detailed protocolSiMS experiments using Dylight 550 labeled actin introduced to XTC cells by electroporation was performed as in Yamashiro et al., 2014. Cells adhered on a poly-lysine-coated glass coverslip were observed by epi-fluorescence microscopy. For measurements of speckle lifetime distribution, the exposure time was 0.1 s/frame. The lifetime data were normalized for photobleaching as in Watanabe and Mitchison, 2002. In Figure 6A experiments, the leading edge was manually marked by a line which connects the centers of furthermost speckles in several consecutive images. Only actin SiMS that appeared within ~300 nm from the line were analyzed for lifetimes. Lifetimes of 0.1 s were omitted as being beyond the temporal resolution limit. Disassembly of actin speckles in lamellipodia of XTC cells were also observed after treatment with 5 µM cytochalasin D, added after several seconds from the start of observation, at 2 s intervals.
Parameter scan for severing rates and oligomer size
Request a detailed protocolTo estimate parameters for uniform and enhanced end severing rates (, ), we performed a parameter search over these values as well as the maximum oligomer length . To summarize the results, we classified each set of these parameters in terms of their ability to match the actin network structure and concentration, actin speckle dynamics, and if there was an increase in length between the front and back of the lamellipodium (Figure 3—figure supplement 1, Figure 4—figure supplement 1, and Figure 5—figure supplement 1):
Matching actin network structure and concentration. We marked a parameter set as satisfying this condition when the following occurred: (i) The F-actin concentration near the leading edge was 700–1500 µM for keratocytes (Schaub et al., 2007; Mueller et al., 2017) and 400–1300 µM for XTC parameters (Watanabe, 2010; Vinzenz et al., 2012); (ii) The midpoint of the concentration profile fell between 5–15 µm (keratocytes) and 3–10 µm (XTC cells); (iii) The barbed end concentration profiles were within 50% of measurements of keratocyte lamellipodia (Mueller et al., 2017). For XTC cells, the barbed end and branch concentrations were at least as large as those measured in fibroblast lamellipodia (Vinzenz et al., 2012).
Matching actin speckle dynamics. We indicated agreement with actin speckle SiMS (Watanabe and Mitchison, 2002; Yamashiro et al., 2018; Smith et al., 2013) when the following conditions were met: (i) Actin speckle lifetimes peak at short times, satisfied when the probability of speckles with lifetimes between 2 and 4 s was 0.025–0.055 for keratocytes (from the full distribution ranging from 2 to 60 s) and the probability of lifetimes between 4 and 8 s for XTC cells was 0.035–0.08 (from the full distribution between 4–144 s). We did not include lifetimes shorter than 2 or 4 s, respectively, that could be at the limits of experimental resolution; (ii) The lifetime distribution extended to lifetimes longer than 10 s (at least 20%); (iii) The normalized speckle appearance and disappearance profiles as a function of distance from the leading edge were consistent with SiMS. In simulations without annealing, the speckle lifetime distribution did not peak monotonically at the shortest lifetimes, so we separately marked parameter sets with a peak within the first 15 s (keratocytes) or 35 s (XTC).
Increase in length away from leading edge. We considered length increase to occur if the filament length at 50% of the cumulative length distribution of filaments with segments between 3–4 µm away from the leading edge was at least 0.2 µm (keratocytes) or 0.1 µm (XTC) larger compared to 0–1 µm away from the leading edge.
Estimation of branching and capping rates by comparison to prior electron tomograms
Capping and branching rates were determined by comparing to the the barbed end, branch and filament number of electron micographs of fibroblast cells in Vinzenz et al., 2012 (used for the XTC parameter set) and barbed end, pointed end and filament number from Mueller et al., 2017 for the keratocyte parameter set. In all instances of data comparison for the barbed end and pointed end number of Mueller et al., 2017, the average value of two 0.106 µm bins was compared to simulation results of 0.212 µm bins.
We use a simplified 2D dendritic nucleation model, without severing and annealing, to estimate the capping and branching rates for our simulations of keratocyte and XTC cases. We assume linear actin filaments, all oriented ±35° relative to the leading edge with the direction toward the center of the cell. These filaments polymerize, branch, cap and move with respect to the position of the leading edge due to retrograde flow or cell protrusion. Filaments with free barbed ends (uncapped) remain at the leading edge and continue to polymerize. As they cap, they stop polymerization and move toward the center of the cell with retrograde flow. For the purposes of this section, we do not consider uncapping of already capped filaments. Branches form at the leading edge and also move away with retrograde flow. Denoted by the length of a filament with barbed end located at , the system of equations for the number of uncapped , capped barbed ends and branches are:
The first term on the right hand side of Equation 1 represents the branching source term which creates a new filament at 0 length. The last term of Equation 1 is polymerization of the filaments with . The second term of Equation 1 and first term of Equation 2 represents the loss of uncapped filaments and addition of capped filaments. The second term of Equation 2 represents retrograde flow of capped filaments. Equation 3 accounts for the generation of branches and their motion with retrograde flow toward the center of the cell.
In steady-state, the solutions of the uncapped and capped ends to be exponential and the number of branches is:
where . To compare to electron tomogram quantifications we calculate the number of branches, barbed ends and the number of filaments in bins of increasing distance to the leading edge.
Integrating Equation 6 over from the minimum to maximum of the bin we find the number of branches within the bin, .
where is the bin width centered at such that is the bin that includes the leading edge.
The number of barbed ends within a bin centered at depends on the total number of capped and uncapped filaments. We integrate over all filament lengths and the distance from the leading edge :
We find the solution for the number of barbed ends in the the bin to be:
In the electron tomograms, the filament number was measured by counting the number of filaments crossing the middle plane of the measured bins. We calculate the equivalent filament number by integrating the sum of the capped and uncapped filaments that are long enough to cross the distance in at the center of such bins:
Integrating, we find,
which can be evaluated at bin using .
Parameter estimation for keratocyte parameter set
Request a detailed protocolMueller et al., 2017 quantified the number of barbed, pointed ends and filaments in keratocyte lamellipodia near the leading edge. Assuming all pointed ends near the leading edge are at a branch, we can solve Equation 7 for the branching rate using the relative network velocity , and an estimate for the number of branches, . Unlike Equation 7 in the simple model considered in this supplementary text, the number of pointed ends in Mueller et al., 2017 increased over the first (two bins). Using the average number of pointed ends over two bins from the leading edge as an estimate of in Equation 7 leads to . This number is consistent with the value calculated using the number of pointed ends at longer distances, beyond the first two bins.
Using this estimate for , we can estimate in two different ways:
(i) Using Equation 9 with and Table 2 we find the capping rate to be .
(ii) Numerically solving Equation 11, using the number of filaments from Table 2 with and from above, we find (we do not consider since the number of filaments is increasing within the region 0–0.106 µm, unlike in the current model).
In the main text and figure supplements we used branching rate and capping rates similar to the values calculated in (i) and (ii): was required to produce consistent results when depolymerization and severing was included in the simulation.
Parameter estimation for XTC parameter set
Request a detailed protocolVinzenz et al., 2012 studied lamellipodia of fibroblast cells and measured branch, barbed end and filament number close to the leading edge, as well as the filament length distribution Table 2. In this system .
The branching rate can be estimated using Equation 7 and Table 2 to find .
We can estimate the capping rate in three different ways as follows.
(i) Using the calculation for the number of barbed ends of Equation 9. From Equation 9, the value of barbed ends in bins 2 and higher is of the same order as the number of branches in Equation 7, consistent with the measurements in Vinzenz et al., 2012 in Table 2. The ratio of barbed ends in the first bin to the second bin is however a value that depends on but is independent of :
Solving for the capping rate using leads to .
(ii) Using the solution for the number of filaments of Equation 11 and the values in Table 2. We find for , and , .
(iii) Comparison to the average filament length . The median filament branch length from Vinzenz et al., 2012 is approximately 162 nm ≈ 60 sub (a value similar to Bailly et al., 1999 who measured filament lengths to be between 100 and 200 nm near the leading edge of MTLn3 cells). Using this value for , we find .
In the main text and figure supplements we used values of , and /s to produce the same F-actin density at the leading edge while using (comparable to retrograde flow in XTC cells) and also producing consistent results when depolymerization and severing is added.
Data availability
All data reported in this project are present within the published figures and Supplemental Information. The code for simulations is available at https://github.com/vavylonis/LamellipodiumSeverAnneal, (copy archived at swh:1:rev:063372789ca6d92073d7db753e39b8a5051fd3d4) and will allow for all simulation plots to be reproduced. The experimental SiMS data of Figure 6 and Figure 6-supplement 1 have been provided as excel files containing the speckle tracks using the SpeckleTrackerJ ImageJ plugin.
-
GitHubID LamellipodiumSeverAnneal. LamellipodiumSeverAnneal.
References
-
Kinetic mechanism of end-to-end annealing of actin filamentsJournal of Molecular Biology 312:721–730.https://doi.org/10.1006/jmbi.2001.5005
-
Actin dynamics, architecture, and mechanics in cell motilityPhysiological Reviews 94:235–263.https://doi.org/10.1152/physrev.00018.2013
-
Global treadmilling coordinates actin turnover and controls the size of actin networksNature Reviews. Molecular Cell Biology 18:389–401.https://doi.org/10.1038/nrm.2016.172
-
Disassembly of actin networks by filament severingNew Journal of Physics 9:418.https://doi.org/10.1088/1367-2630/9/11/418
-
Actin dynamics: from nanoscale to microscaleAnnual Review of Biophysics 39:91–110.https://doi.org/10.1146/annurev.biophys.093008.131207
-
Quantitative fluorescent speckle microscopy of cytoskeleton dynamicsAnnual Review of Biophysics and Biomolecular Structure 35:361–387.https://doi.org/10.1146/annurev.biophys.35.040405.102114
-
An open model of actin dendritic nucleationBiophysical Journal 96:3529–3542.https://doi.org/10.1016/j.bpj.2009.01.037
-
A model for actin-filament length distribution in A lamellipodJournal of Mathematical Biology 43:325–355.https://doi.org/10.1007/s002850100102
-
Growing actin networks regulated by obstacle size and shapeActa Mechanica Sinica 33:222–233.https://doi.org/10.1007/s10409-016-0628-5
-
G-actin regulates rapid induction of actin nucleation by mDia1 to restore cellular actin polymersJournal of Cell Science 121:3403–3412.https://doi.org/10.1242/jcs.030940
-
Mechano-chemical feedbacks regulate actin mesh growth in lamellipodial protrusionsBiophysical Journal 98:1375–1384.https://doi.org/10.1016/j.bpj.2009.11.054
-
Growing actin networks form lamellipodium and lamellum by self-assemblyBiophysical Journal 95:5508–5523.https://doi.org/10.1529/biophysj.108.134817
-
High-speed depolymerization at actin filament ends jointly catalysed by Twinfilin and Srv2/CAPNature Cell Biology 17:1504–1511.https://doi.org/10.1038/ncb3252
-
Lamellipodium tip actin barbed ends serve as a force sensorGenes to Cells 24:705–718.https://doi.org/10.1111/gtc.12720
-
Actin disassembly by cofilin, coronin, and Aip1 occurs in bursts and is inhibited by barbed-end cappersThe Journal of Cell Biology 182:341–353.https://doi.org/10.1083/jcb.200801027
-
Arp2/3 complex interactions and actin network turnover in lamellipodiaThe EMBO Journal 27:982–992.https://doi.org/10.1038/emboj.2008.34
-
Self-organization of treadmilling microtubules into a polar arrayTrends in Cell Biology 12:462–465.https://doi.org/10.1016/s0962-8924(02)02369-3
-
Cofilin increases the bending flexibility of actin filaments: implications for severing and cell mechanicsJournal of Molecular Biology 381:550–558.https://doi.org/10.1016/j.jmb.2008.05.055
-
A model actin comet tail disassembling by severingPhysical Biology 8:046003.https://doi.org/10.1088/1478-3975/8/4/046003
-
Diffusion, capture and recycling of SCAR/WAVE and Arp2/3 complexes observed in cells by single-molecule imagingJournal of Cell Science 125:1165–1176.https://doi.org/10.1242/jcs.091157
-
Can filament treadmilling alone account for the F-actin turnover in lamellipodia?Cytoskeleton (Hoboken, N.J.) 70:179–190.https://doi.org/10.1002/cm.21098
-
Regulation of actin dynamics in rapidly moving cells: a quantitative analysisBiophysical Journal 83:1237–1258.https://doi.org/10.1016/S0006-3495(02)73897-6
-
Xenopus actin-interacting protein 1 (XAip1) enhances cofilin fragmentation of filaments by capping filament endsThe Journal of Biological Chemistry 277:43011–43016.https://doi.org/10.1074/jbc.M203111200
-
Loss of Aip1 reveals a role in maintaining the actin monomer pool and an in vivo oligomer assembly pathwayThe Journal of Cell Biology 188:769–777.https://doi.org/10.1083/jcb.200909176
-
Two distinct actin networks drive the protrusion of migrating cellsScience (New York, N.Y.) 305:1782–1786.https://doi.org/10.1126/science.1100533
-
Crowded surfaces change annealing dynamics of actin filamentsJournal of Molecular Biology 368:365–374.https://doi.org/10.1016/j.jmb.2007.01.087
-
Actin Turnover in Lamellipodial FragmentsCurrent Biology 27:2963–2973.https://doi.org/10.1016/j.cub.2017.08.066
-
Excitable Actin Dynamics in Lamellipodial Protrusion and RetractionBiophysical Journal 102:1493–1502.https://doi.org/10.1016/j.bpj.2012.03.005
-
Analysis of actin filament network organization in lamellipodia by comparing experimental and simulated imagesJournal of Cell Science 120:1491–1500.https://doi.org/10.1242/jcs.03379
-
Direct observation of tropomyosin binding to actin filamentsCytoskeleton (Hoboken, N.J.) 72:292–303.https://doi.org/10.1002/cm.21225
-
Twinfilin bypasses assembly conditions and actin filament aging to drive barbed end depolymerizationThe Journal of Cell Biology 220:e202006022.https://doi.org/10.1083/jcb.202006022
-
Cytochalasin D acts as an inhibitor of the actin-cofilin interactionBiochemical and Biophysical Research Communications 424:52–57.https://doi.org/10.1016/j.bbrc.2012.06.063
-
Distributed actin turnover in the lamellipodium and FRAP kineticsBiophysical Journal 104:247–257.https://doi.org/10.1016/j.bpj.2012.11.3819
-
Analysis of the actin-myosin 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
-
Catastrophic actin filament bursting by cofilin, Aip1, and coroninJournal of Biological Chemistry 295:13299–13313.https://doi.org/10.1074/jbc.RA120.015018
-
Actin branching in the initiation and maintenance of lamellipodiaJournal of Cell Science 125:2775–2785.https://doi.org/10.1242/jcs.107623
-
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
-
Single-molecule speckle analysis of actin filament turnover in lamellipodiaScience (New York, N.Y.) 295:1083–1086.https://doi.org/10.1126/science.1067470
-
Inside view of cell locomotion through single-molecule: fast F-/G-actin cycle and G-actin regulation of polymer restorationProceedings of the Japan Academy. Series B, Physical and Biological Sciences 86:62–83.https://doi.org/10.2183/pjab.86.62
-
New single-molecule speckle microscopy reveals modification of the retrograde actin flow by focal adhesions at nanometer scalesMolecular Biology of the Cell 25:1010–1024.https://doi.org/10.1091/mbc.E13-03-0162
-
Myosin-dependent actin stabilization as revealed by single-molecule imaging of actin turnoverMolecular Biology of the Cell 29:1941–1947.https://doi.org/10.1091/mbc.E18-01-0061
Article and author information
Author details
Funding
National Institutes of Health (R35GM136372)
- Danielle Holz
- Aaron R Hall
- Eiji Usukura
- Sawako Yamashiro
- Naoki Watanabe
- Dimitrios Vavylonis
National Institutes of Health (R01GM114201)
- Danielle Holz
- Aaron R Hall
- Dimitrios Vavylonis
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank David Rutkowski for helping with network visualization. This work was supported by National Institutes of Health Grant R01GM114201 and R35GM136372. Use of the high-performance computing capabilities of the Extreme Science and Engineering Discovery Environment (XSEDE), which is supported by the National Science Foundation, project no. TG-MCB180021 is also gratefully acknowledged.
Copyright
© 2022, Holz et al.
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
-
- 1,038
- views
-
- 222
- downloads
-
- 4
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
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
- Physics of Living Systems
The endoplasmic reticulum (ER), the largest cellular compartment, harbours the machinery for the biogenesis of secretory proteins and lipids, calcium storage/mobilisation, and detoxification. It is shaped as layered membranous sheets interconnected with a network of tubules extending throughout the cell. Understanding the influence of the ER morphology dynamics on molecular transport may offer clues to rationalising neuro-pathologies caused by ER morphogen mutations. It remains unclear, however, how the ER facilitates its intra-luminal mobility and homogenises its content. It has been recently proposed that intra-luminal transport may be enabled by active contractions of ER tubules. To surmount the barriers to empirical studies of the minuscule spatial and temporal scales relevant to ER nanofluidics, here we exploit the principles of viscous fluid dynamics to generate a theoretical physical model emulating in silico the content motion in actively contracting nanoscopic tubular networks. The computational model reveals the luminal particle speeds, and their impact in facilitating active transport, of the active contractile behaviour of the different ER components along various time–space parameters. The results of the model indicate that reproducing transport with velocities similar to those reported experimentally in single-particle tracking would require unrealistically high values of tubule contraction site length and rate. Considering further nanofluidic scenarios, we show that width contractions of the ER’s flat domains (perinuclear sheets) generate local flows with only a short-range effect on luminal transport. Only contractions of peripheral sheets can reproduce experimental measurements, provided they are able to contract fast enough.
-
- Cell Biology
- Developmental Biology
How cells regulate the size of their organelles remains a fundamental question in cell biology. Cilia, with their simple structure and surface localization, provide an ideal model for investigating organelle size control. However, most studies on cilia length regulation are primarily performed on several single-celled organisms. In contrast, the mechanism of length regulation in cilia across diverse cell types within multicellular organisms remains a mystery. Similar to humans, zebrafish contain diverse types of cilia with variable lengths. Taking advantage of the transparency of zebrafish embryos, we conducted a comprehensive investigation into intraflagellar transport (IFT), an essential process for ciliogenesis. By generating a transgenic line carrying Ift88-GFP transgene, we observed IFT in multiple types of cilia with varying lengths. Remarkably, cilia exhibited variable IFT speeds in different cell types, with longer cilia exhibiting faster IFT speeds. This increased IFT speed in longer cilia is likely not due to changes in common factors that regulate IFT, such as motor selection, BBSome proteins, or tubulin modification. Interestingly, longer cilia in the ear cristae tend to form larger IFT compared to shorter spinal cord cilia. Reducing the size of IFT particles by knocking down Ift88 slowed IFT speed and resulted in the formation of shorter cilia. Our study proposes an intriguing model of cilia length regulation via controlling IFT speed through the modulation of the size of the IFT complex. This discovery may provide further insights into our understanding of how organelle size is regulated in higher vertebrates.