<title>Introduction</title> <p>Actin polymerization at the plasma membrane results in the formation of cellular protrusions known as lamellipodia or filopodia, which mediate cell migration ##REF##9013669##[1]##–##REF##15306794##[3]##. The distinct organization and generation of filaments in each structure uses a different mechanism to produce mechanical force. In the lamellipodia, the actin filaments organize into a flat 2D branched network ##REF##12600310##[2]##, ##REF##10352018##[4]## whereas in the filopodia they are assembled into long, parallel, closely packed bundles ##REF##12566431##[5]##–##REF##1447299##[7]##. Different proteins control the assembly of these structures; in the lamellipodia, the branched nucleation is driven by activation of the Arp2/3 complex ##REF##10801131##[8]## by Wiskott-Aldrich syndrome protein family ##REF##9889097##[9]##, ##REF##11329366##[10]## (WASP), followed by filament elongation and barbed-end capping by capping proteins (CP) ##REF##10679358##[11]##. Formin and Ena/VASP proteins are concentrated at the tips of filopodia ##REF##17475772##[12]##–##REF##15908944##[15]## promoting the continuous elongation of filaments ##REF##15306794##[3]##, ##REF##15294161##[16]## and their successive binding by fascin, which crosslinks actin filaments of the same polarity into the bundles constituting the “body” of filopodia ##REF##11948621##[17]##. In both structures, the barbed (i.e., the fast growing) ends of actin filaments point towards the plasma membrane ##REF##10477762##[18]##.</p> <p>Most cultured animal cells assemble both lamellipodia and filopodia. Some cells, like dendritic cells, are dominated by filopodia ##REF##15294161##[16]## while keratocytes grow exclusively lamellipodia ##REF##8968574##[19]##. It is still not fully understood what determines the <italic>in vivo</italic> preference to lamellipodia vs. filopodia, and why certain cells are dominated by only one type of structure. So far, two mechanisms were suggested to explain the formation of filopodia ##REF##17712139##[20]##. In the “<italic>de novo</italic> filopodia nucleation” model, actin filament nucleation and elongation is mediated by Dia2 proteins ##REF##15908944##[15]##, ##REF##16740473##[21]##. According to this mechanism bundles do not emanate from lamellipodia and Arp2/3 complex is not essential for filopodia formation. In the second mechanism, known as the “convergent elongation” model, filopodia emerge from the Arp2/3 lamellipodial network ##REF##12642617##[22]##–##REF##18256280##[25]##. During formation of filopodia, lamellipodial filaments associate at their barbed ends with mDia2 or VASP, elongate continuously and gradually converge into filopodia bundles, by the cross-link of fascin. In contrast to the first mechanism, here filopodia are anchored to lamellipodia. Both systems can be reconstituted <italic>in vitro</italic>. This work primarily focuses on the “convergent elongation” model.</p> <p>Notwithstanding recent advances in understanding the dynamic organization of lamellipodia and filopodia protrusions, it is still not fully understood how cells control the transition between these structures, what directs the localization of filopodia formation along the cell leading edge, and how their thickness is regulated. It is expected that the emergence of filopodia from lamellipodia would be strongly affected by the properties (e.g., the density and length of filaments) of the branched lamellipodial network from which they emanate. The concentration of Arp2/3, which strongly affects these structural properties, was recently observed to have a dramatic influence on the formation of filopodia <italic>in vivo</italic>
##REF##18256280##[25]##. The abundance of fascin, which is the driving force for bundling is also expected to affect bundles formation. Yet, it is still not clear how the structure of lamellipodia and the fascin concentration control filopodia formation.</p> <p>To resolve these issues, it is essential to understand the mechanisms underlying the formation of lamellipodia and filopodia and evaluate the factors controlling their structure and dynamics. To this end, in the present work we study a system containing the major proteins participating in the assembly of lamellipodia and filopodia: a) the constitutively active VCA ##REF##12015607##[26]## (or WA ##REF##10801131##[8]##) domain of WASp; b) Arp2/3 complex; c) fascin, and d) actin monomers. This model system enables us to control the concentrations of the various proteins, and to observe actin organization in bulk. Recently, we have demonstrated the feasibility of using such a simple reconstituted system for studying the formation of lamellipodia and filopodia-like structures <italic>in vitro</italic>
##REF##16549794##[27]##. In this early work, the roles of Arp2/3 complex and fascin in the transition from the Arp2/3 branched network to filopodia-like bundles was addressed only briefly. Also, no attempt was made to model this transition nor to resolve the primary factors controlling the bundling process.</p> <p>Our goal in the present work is to study in detail the self-assembly characteristics of actin in the presence of variable amounts of Arp2/3 complex and fascin. We will show that in the absence of fascin, actin organizes into dense 3D, aster-like, structures composed of branched networks of actin filaments ##REF##16549794##[27]##. The growth of the network advances with the barbed ends of the peripheral actin filaments pointing on average outward, as also demonstrated in the simulations ##REF##16549794##[27]##. Addition of fascin leads to the formation of star-like structures, where actin/fascin bundles are nucleated and emanate radially from the branched network core. The polarity of the actin filaments during the transition is preserved, as in the transition from lamellipodia to filopodia ##REF##16549794##[27]##. Our experiments show that in the early stages of actin polymerization, Arp2/3 mediates the formation of the dense and highly branched aster-like network, whereas fascin is rather passive. Fascin becomes dominant at a later stage of actin polymerization, when the (sparsely branched) filaments in the periphery of the aster become long enough, so that the energy gained due to actin-fascin-actin links between neighboring filaments overcomes the unfavorable bending energy required to bring them into contact. As time evolves the initially thin bundles elongate, thus lowering their bending energy thereby enabling their association into thicker bundles, and so on until no more actin monomers are left.</p> <p>This picture of stepwise bundle elongation and thickening underlies our structural-energetic model of bundle formation and growth. Using reasonable approximations for the bending energies of actin filaments and bundles, and the cohesive energy due to fasin-actin bonds we can explain the dependence of the thickness and length of the mature bundles upon [Arp2/3]. More explicitly, using Kinetic Monte Carlo (KMC) simulations we can model the nucleation and growth of asters and derive their structure as a function of time thus obtaining two central structural parameters: the length, <italic>L</italic>, of filament tips at the aster periphery, and the spacing, <italic>b</italic>, between their origins (i.e., their pointed ends, anchored at Arp2/3 branch points). These parameters, which determine the initial conditions at the onset of bundling, depend sensitively on [Arp2/3]. As we shall see, our theoretical scheme, which combines the KMC and the stepwise bundle assembly model, accounts adequately for the experimentally observed dependence of bundle structure on the initial concentration of Arp2/3 in the system.</p>
<title>Materials and Methods</title> <title>Protein purification</title> <p>Actin was purified from rabbit skeletal muscle acetone powder ##REF##4254541##[33]##. Recombinant WASp-VCA ##REF##12015607##[26]## and fascin (recombinant fascin was prepared by a modification of the method of ##REF##8999969##[34]##) were expressed in <italic>E. coli</italic> as GST fusion proteins. Arp2/3 complex was purchased from Cytoskeleton Co. Actin was labeled on Cys374 with Alexa fluor 568 or 488 (Invitrogen, Co.).</p> <title>Motility assay</title> <p>The motility medium contained 10 mM HEPES, pH = 7.7, 1.7 mM Mg-ATP, 5.5 mM DTT, 0.12 mM Dabco [an anti-bleaching reagent], 0.1 M KCl, 1 mM MgCl<sub>2</sub>, 1% BSA and various concentrations of G-actin, Arp2/3 complex, VCA, and fascin.</p> <title>Microscopy</title> <p>Actin assembly was monitored for about an hour by fluorescence with an Olympus IX-71 microscope. The labeled actin fraction was 1/40 and the temperature was 22°C. Time-lapse images were acquired using an Andor DV887 EMCCD camera (Andor Co., England). Data acquisition and analysis was performed using METAMORPH (Universal Imaging Co.).</p>
<title>Results</title> <p>The first part of this section describes the three types of structures, i.e., asters, stars and network of bundles (a phase of actin-fascin bundles) that were experimentally observed and the transition between them. We also present a phase diagram of the system to illustrate the regions of existence of the three types of structures as a function of the initial [Arp2/3] and [fascin]. In the second part, we focus specifically on the transition from asters to stars, as it mimics the transition from lamellipodia to filopodia in cells. Special attention is given to evaluating the roles of [Arp2/3] and [fascin] on filopodia formation and dimensions. In the last part of this section we describe our theoretical model for bundle assembly and apply it to explain our experimental results.</p> <title>Fascin mediates structural transition of actin self-assembly</title> <p>First, we aimed at studying the self–assembly of actin in the presence of increasing amounts of fascin, at a given [Arp2/3]. At low Arp2/3 complex concentration of 6.25 nM (the [VCA]/[Arp2/3] ratio was kept equal to two in all experiments) two structural transitions are observed, the first, from ‘asters’ to ‘stars’, and the second, from ‘stars’ to ‘network of bundles’ (##FIG##0##Fig. 1##). At a very low fascin concentration (3 nM) the structure of actin network is dominated by Arp2/3 nucleation and branching activity–the system organizes into diffuse highly branched 3D aster-like structures (##FIG##0##Fig. 1a##), similar to those formed in the absence of fascin ##REF##16549794##[27]##. Surprisingly, a minute concentration of fascin ([fascin] = 4 nM) was sufficient to mediate the nucleation and growth of bundles emanating from branched asters' cores. ##FIG##0##Fig. 1b## shows such a star-like structure after reaching final dimensions. The bundles grow with their barbed ends pointing outward, as previously demonstrated ##REF##16549794##[27]##. The time span of this process ranges between few to ten minutes, depending on [Arp2/3] and [fascin] (data not shown). In this work all the experimental data analyses were performed for systems that have reached their final state.</p> <p>Gradual addition of fascin, from 0.005 to 0.2 µM, increases the number of bundles emanating from the asters core as well as their length (##FIG##0##Figs. 1c–1f##). Above 0.2 µM of fascin, the system is dominated by the bundling activity of fascin. Most of the actin appears to be located in actin/fascin bundles residing in the bulk solution (see image background, ##FIG##0##Fig. 1g##; [fascin = 0.5 µM]). The stars that still exist are composed of fewer and shorter bundles. Eventually, at high enough [fascin] (3 µM) stars no longer form; instead, a network of actin/fascin bundles is generated (##FIG##0##Fig. 1h##).</p> <title>[Arp2/3] controls fascin ability to induce structural changes</title> <p>Structural transformations mediated by fascin seem to be facilitated when the initial [Arp2/3] was reduced. The top and bottom panels in ##FIG##1##Fig. 2## distinguish between fluorescence microscopy images corresponding to mixtures containing low ([Arp2/3] = 12.5 nM) and high ([Arp2/3] = 100 nM) initial Arp2/3 contents, respectively. In both cases [fascin] varies in the same manner. Initially, at a very low [fascin] of 3 nM (##FIG##1##Figs. 2a and 2e##), actin organization is fully dominated by Arp2/3 complex activity; diffuse asters are formed, whose number increases with the initial [Arp2/3]. These asters are smaller and their branched actin cores are denser (##FIG##1##Fig. 2e##) compared to those formed at the lower [Arp2/3] (##FIG##1##Fig. 2a##). Increasing the concentration of fascin to 6 nM (blue arrow) induced a transition from asters to stars when [Arp2/3] = 12.5 nM (##FIG##1##Fig. 2b##; ##SUPPL##1##Figs. S1a and S1b## show that the transition occurs between 5 nM and 6 nM of fascin). On the other hand, 6 nM fascin is not sufficient to induce bundle formation in the [Arp2/3] = 100 nM system (##FIG##1##Fig. 2f##). In this latter case, more fascin (7 nM) is required to nucleate bundles from the branched network core (see also ##SUPPL##1##Figs. S1c and S1d##).</p> <p>Further addition of fascin induces the formation of additional and longer bundles, resulting in the appearance of fully developed stars, as shown in ##FIG##1##Figs. 2c and 2g##. The total number of bundles nucleated per star in the low [Arp2/3] (12.5 nM) case appears to be larger than in the high [Arp2/3] (100 nM) case. Moreover, the bundles in the low [Arp2/3] system look both thicker (i.e., they comprise a larger number of actin filaments) and longer. Clearly then, the initial concentration of Arp2/3 in solution plays a major role in determining the conditions for the fascin-mediated assembly of actin bundles, thus (albeit indirectly) controlling bundles thickness and length.</p> <p>Finally, when the amount of fascin was further elevated to 3 µM (red arrow), a transition from ‘stars’ to ‘a network of bundles’ was observed in the 12.5 nM Arp2/3 system (##FIG##1##Fig. 2d##). In contrast, the system with 100 nM Arp2/3 remained unchanged (##FIG##1##Fig. 2h##) and did not show any structural transition up to 6 µM fascin (data not shown). Actually ‘networks of bundles’ did not appear in any of the systems containing more than 25 nM of Arp2/3 complex, regardless of [fascin].</p> <p>In ##FIG##2##Fig. 3## we present a phase diagram of the system, depicting the regions of existence of the three types of experimentally observed structures, as a function of the initial [Arp2/3] and [fascin]. The data shown here corresponds to experiments conducted at [G-actin] = 7 µM; a similar qualitative behavior was observed for the other concentrations of actin we have tested, ranging between 1 µM and 5 µM (data not shown). Careful inspection shows that below 3 nM of fascin only aster-like structures form (circles), regardless of [Arp2/3]. Increasing [fascin], for a given amount of Arp2/3 complex, induces structural transitions into stars (inverted triangles). The amount of fascin required to mediate this transition increases nonlinearly with [Arp2/3], saturating above 40 nM Arp2/3, (the blue curve marks the approximate boundary between the two regions). While only few nano-molars of fascin are sufficient to induce the transition from asters to stars, the subsequent transition to ‘networks’ (squares) requires much larger amounts of fascin (several µM; the red curve marks the approximate boundary between the two regions). We note that a transition to ‘networks’ is observed only for sufficiently low [Arp2/3] (&lt;25 nM), within this region, the amount of fascin required for the nucleation of bundles increases monotonically with Arp2/3 content.</p> <title>Aster to star transition–Experimental results</title> <title>The densities of asters and stars depend on [Arp2/3] but not on [fascin]</title> <p>To test our assumption that each star originates from a preformed aster, we have measured the number of “aggregates” (i.e., asters or stars) per unit area (surface density), within the 1 µm section of the bulk sample observed. The surface density of aggregates <italic>σ</italic> = <italic>d</italic>
<sup>−2</sup> was evaluated by measuring the average distance, <italic>d</italic>, between adjacent aggregates. Keeping a constant ratio [VCA]/[Arp2/3] = 2, we have measured <italic>d</italic>
<sup>−2</sup> for several values of the initial [Arp2/3] and [fascin] (##FIG##3##Fig. 4##). For clarity we present only two concentrations of fascin: 5 nM (triangles) and 0.2 µM (circles). Data points taken at 5 nM fascin correspond to density of asters, while those taken at 0.2 µM fascin are associated with density of stars. In both cases, <italic>d</italic>
<sup>−2</sup> increases monotonically with [Arp2/3]. In contrast, [fascin] does not appear to affect <italic>d</italic>
<sup>−2</sup> (see inset; conditions: [Arp2/3] = 12.5 nM), indicating that the density of aggregates is dictated primarily by the amount of Arp2/3 complex in the system. These results support our previous findings that Arp2/3 complexes are aster nucleators and that stars develop from preformed asters ##REF##16549794##[27]##. This is consistent with our present conclusion that, at early times, actin polymerization is governed by Arp2/3 nucleation and branching. Fascin comes into play and becomes dominant in star formation only later, when the structural properties of the aster enable the onset of filament bundling.</p> <title>The length of bundles in stars depends on Arp2/3 and fascin concentrations</title> <p>For all the concentrations of Arp2/3 that we have tested, the mean final length of bundles, <italic>L<sub>Fin</sub></italic>, was found to increase with [fascin], reaching asymptotically a maximal ([Arp2/3] dependent) value,  (##FIG##4##Fig. 5a##). The system appears more sensitive to fascin addition when [Arp2/3] is lower (as reflected by the steeper slopes in ##FIG##4##Fig. 5a##). Consequently, the amount of fascin required to reach  decreases with the decrease in [Arp2/3]. In contrast, for a fixed concentration of fascin, we found that <italic>L<sub>Fin</sub></italic> decreases monotonously with [Arp2/3], (as seen for [fascin] = 100 nM in ##FIG##4##Fig. 5b##). The inverse correlation between <italic>L<sub>Fin</sub></italic> and [Arp2/3] results from the reduced amount of actin monomers available for bundle growth when [Arp2/3] is increased, as discussed in more detail in the theory section below.</p> <title>Theoretical model and analysis</title> <title>Bundle formation: A simple model and Kinetic Monte Carlo Simulations</title> <p>Our experiments show that upon mixing actin monomers with Arp2/3 complex, VCA and fascin molecules, the first structures to appear in solution are dense and highly branched aster-like networks of actin filaments ##REF##16549794##[27]##. The asters grow radially outwards, and since Arp2/3 joins both old and newly formed filaments, the density of branching points and hence also of polymerized actin is highest in the center of the aster, decreasing gradually towards its periphery, where the average spacing between adjacent branches is relatively large ##REF##16549794##[27]##. To associate into bundles, nearby filaments are generally required to bend towards each other. Inside the aster core, owing primarily to excluded volume interactions between the rigid, dense and highly branched network, filament bending involves a prohibitively large energetic penalty and bundle nucleation is thus highly unlikely.</p> <p>Bending the filaments emanating from the surface of the aster is easier, yet two conditions must be fulfilled to enable the onset of bundle formation. First, <italic>L</italic>, the average length of the filament tips in the aster periphery, should be larger than the distance, <italic>b</italic>, between their pointed ends (which mark their origins). The second requirement is that the energy gained due to fascin links between neighboring filaments should exceed the bending energy associated with bringing them into contact. Both conditions are met at some time <italic>t<sub>i</sub></italic>, in the late stages of aster growth, owing mainly to the decrease in [Arp2/3] and hence the increase in <italic>L</italic>. Qualitatively, this explains why in the early stages of actin polymerization the (dense aster-like) structure of the actin network is dominated by Arp2/3, with fascin playing a rather passive role. However, once the peripheral filaments are long enough, fascin comes into play–linking the filaments into bundles, whereas branching essentially ceases; especially along and within the bundles, because the (already largely depleted) Arp2/3 complexes cannot penetrate the dense bundle.</p> <p>Our goal in this section is to explain the experimental observation that bundles are thicker and longer when the initial [Arp2/3] is low (see ##FIG##1##Figs. 2## and ##FIG##4##5##). Qualitatively, this behavior is understood based on the fact that since Arp2/3 is an aster nucleator, so that for a given concentration of G-actin, higher [Arp2/3] results in the appearance of more asters in solution. Thus, on average, fewer actin monomers participate in the growth of each aster. Also, since the asters formed are denser, fewer monomers are later available for bundle growth, resulting in relatively short bundles, as observed in our experiments (##FIG##1##Figs. 2## and ##FIG##4##5##). It should be noted that in living cells, the plasma membrane also affects the length of filopodia, as was previously analyzed theoretically ##REF##15879474##[23]##, ##REF##16214866##[28]##. Furthermore, in our system, owing to the large energy penalty associated with the bending of short bundles towards each other, and because the bending energy increases rapidly with bundle thickness (see below), when [Arp2/3] is high, association of bundles is energetically costly, and the mature bundles are both shorter and thinner than in the case of low [Arp2/3].</p> <p>The structural properties of the asters at the onset of bundle formation, primarily the distances <italic>b</italic> and <italic>L</italic> defined above, enter our model as input parameters. Numerical estimates of these quantities were derived using our 3D KMC simulations of aster growth. To derive the aster structural characteristics relevant to the experiments reported here, we have used a limited version of the simulation, which includes only actin polymerization and Arp2/3 mediated branching ##REF##16549794##[27]##, (see also ##SUPPL##0##Appendix S1##, in “Supporting information”). A simulation snapshot of an aster is shown in ##FIG##5##Fig. 6a##. It should be noted that (at least presently) these simulations do not account for filament bending, and thus cannot describe bundle formation; they are used here only to provide the initial <italic>L</italic> and <italic>b</italic>.</p> <p>At the onset of the aster-to-star transition the radius of the aster is generally much larger than both <italic>b</italic> and <italic>L</italic>. Since bundles are formed by the association of neighboring filaments whose origins are just a few <italic>b</italic>'s apart, we ignore the curvature of the aster's surface and treat it as being planar. It may be noted, however, that once bundles become much larger than the aster core, surface curvature correction should (and can) be taken into account.</p> <p>Assuming that all the <italic>M</italic> filaments participating initially in bundle formation are of the same length <italic>L</italic>, and organize into <italic>M/N</italic> identical bundles, each of which consists of <italic>N</italic> filaments, the total energy of the system is given by <italic>E<sub>tot</sub></italic> = (<italic>M</italic>/<italic>N</italic>)<italic>E<sub>bundle</sub></italic>(<italic>N</italic>;<italic>L</italic>,<italic>b</italic>). The bundle energy <italic>E<sub>bundle</sub></italic>(<italic>N</italic>;<italic>L</italic>,<italic>b</italic>) includes the cohesive energy due to the fascin-mediated linking of filaments, the bending energy cost of bringing the filaments into contact, and the (unfavorable) surface energy of a finite size bundle. All these contributions depend on <italic>N</italic>, and parametrically on <italic>L</italic> and <italic>b</italic> (Additional contributions, such as the entropy loss experienced by the filaments upon bundling can be neglected). If bundle elongation were slow compared to filament (and bundle) binding-unbinding events, then the optimal bundle size <italic>N</italic>
<sup>*</sup> could be derived by minimizing <italic>E<sub>tot</sub></italic> with respect to <italic>N</italic>. Using reasonable models for the bending, cohesion and surface energies one can show that the resulting <italic>N</italic>
<sup>*</sup> is indeed larger when [Arp2/3] is lower. It must be noted, however, that in this scheme the optimal size and length of bundles is determined by the principles of equilibrium thermodynamics. However, kinetic estimates of the rates of bundle dissociation (based on calculations of filaments dimerization time, ##REF##16565053##[29]##, ##REF##17930069##[30]## (data not shown) indicate that once bundles are formed, their life time is significantly longer than the time scale of actin polymerization, so that filament and bundle associations are practically irreversible on the time scale of our present experiments.</p> <p>Based on this notion we propose here an alternative, stepwise, mechanism whereby single filaments first associate into thin bundles (##FIG##5##Fig. 6b##). The nascent thin bundles are hard to bend, but after elongating their bending and hence fascin-mediated association into thicker bundles becomes easier. This scenario may continue, leading to stepwise bundle thickening. It should be noted that <italic>in vivo</italic> observations of bundles roots also indicate that bundles initiate in a gradual process of elongation and thickening ##REF##12566431##[5]##. In our <italic>in vitro</italic> experiments, this process terminates when there are essentially no more actin monomers in the system.</p> <p>Suppose for simplicity that all possible bundles have a circular cross section consisting of a central straight filament surrounded by <italic>s</italic> concentric shells of filaments; <italic>s</italic> = <italic>R</italic>/<italic>d</italic> where <italic>R</italic> is the radius of the bundle, and <italic>d</italic> is the average distance between adjacent filaments. ##FIG##5##Fig. 6c## shows Δ<italic>E</italic>
<sub>dim</sub>(<italic>s</italic>), the energy change corresponding to the dimerization of two bundles of the same size <italic>s</italic> containing <italic>N</italic>≈<italic>πs</italic>
<sup>2</sup> filaments, into a thicker bundle comprising 2<italic>N</italic> filaments. This simplified calculation ignores many association processes other than bundle doubling (e.g., monomer-dimer, monomer-tetramer binding, etc.) Yet, it captures, at least qualitatively, the kinetic-energetic picture of bundle growth.</p> <p>Results are shown for two of the initial [Arp2/3] studied experimentally (12.5 and 100 nM), and several representative values of <italic>L</italic> which, according to our aster growth simulations, represent relevant filament lengths for the time and concentration scales of the experiments. More specifically, the range of <italic>L</italic> considered is <italic>L<sub>min</sub></italic>([<italic>Arp2</italic>/<italic>3</italic>])&lt;<italic>L</italic>&lt;<italic>L<sub>max</sub></italic>([<italic>Arp2</italic>/<italic>3</italic>]), where <italic>L<sub>min</sub></italic>([<italic>Arp2</italic>/<italic>3</italic>]) is the length of filaments at the time <italic>t<sub>i</sub></italic>, marking the onset of bundling, and <italic>L<sub>max</sub></italic>([<italic>Arp2</italic>/<italic>3</italic>]) is the length of bundles when all G-actin was consumed. Our simulations show that there are about 10<sup>6</sup> free monomers per aster at <italic>t<sub>i</sub></italic> (see <xref ref-type="sec" rid="s4">methods</xref>) when [Arp2/3] = 100 nM, compared to ∼10<sup>8</sup> monomers for 12.5 nM. In both cases, the number of filaments in the aster periphery is about 5000 at <italic>t<sub>i</sub></italic>. Thus <italic>L</italic> ranges between 0.8–1.2 µm for 100 nM Arp2/3. In the 12.5 nM Arp2/3 <italic>L</italic> is much larger: 1.5–30 µm.</p> <p>Association of two adjacent bundles of thickness <italic>s</italic> to form a thicker one is energetically favorable when Δ<italic>E</italic>
<sub>dim</sub>(<italic>s</italic>) is negative. ##FIG##5##Fig. 6c## shows, for instance, that for [Arp2/3] = 12.5 nM, bundles of length 9 µm that are composed of 5 shells of filaments will not associate with each other. However, once their length increases to, say, <italic>L</italic> = 12 µm, they tend to stick to each other since Δ<italic>E</italic>
<sub>dim</sub> becomes negative. From ##FIG##5##Fig. 6c## we note that for both values of [Arp2/3] considered, as <italic>L</italic> increases thicker bundles become energetically favorable; the optimal bundle thickness (corresponding to the minimum of Δ<italic>E</italic>
<sub>dim</sub>(<italic>s</italic>)) increases as well. For example, in the case [Arp2/3] = 12.5 nM, when <italic>L</italic> = 9 µm dimerization is favorable for bundles comprising no more than <italic>s</italic> = 4 shells, and the highest energetic gain is obtained upon dimerization of bundles comprising of <italic>s</italic> = 2 shells. On the other hand, when <italic>L</italic> = 12 µm, the maximum number of shells in bundles capable of dimerization to is <italic>s</italic> = 5, and the optimal thickness is obtained for doubling <italic>s</italic> = 3 bundles. ##FIG##5##Fig. 6c## also reveals the important role of [Arp2/3] on the length and thickness of the bundles. When [Arp2/3] = 100 nM a larger fraction of actin monomers (as compared to the case [Arp2/3] = 12.5 nM) are consumed during aster's assembly, leaving less monomers for bundle growth. This fact, as well as the initial <italic>b</italic> and <italic>L</italic>, are responsible for the shorter and thinner bundles corresponding to the high [Arp2/3] system, (see magnified box). In this case, the typical bundle consists of a small (s∼1) number of shells.</p> <p>To derive the results shown in ##FIG##5##Fig. 6c## we expressed the dimerization energy as a sum of bending and fascin linking energies: Δ<italic>E<sub>dim</sub></italic> = <italic>E<sub>bend</sub></italic>+Δ<italic>E<sub>fascin</sub></italic>, both depending on <italic>s</italic>. To calculate the bending energy of a bundle upon joining another one we have used a simple geometrical model, whereby its bent portion, of length <italic>l</italic>, is represented as composed of two oppositely (and moderately) bent, but otherwise identical arcs, joining smoothly each other to an s-like shape (as may be seen in ##FIG##5##Fig. 6b##). It can be shown that this implies <italic>E<sub>bend</sub></italic> = <italic>γ a</italic>
<sup>2</sup>/<italic>l</italic>
<sup>3</sup> (details will be published elsewhere). Here <italic>a</italic> is the distance between the origins of the bundles (<italic>a</italic> = <italic>bs</italic> for two neighboring bundles of size <italic>s</italic>), <italic>γ</italic> = 32<italic>ξkT</italic> where <italic>ξ</italic> is the persistence length of the bundle, <italic>k</italic> is Boltzmann's constant and <italic>T</italic> is the absolute temperature. Since fascin is abundant in the system, the actin filaments within bundles are tightly bound, and the bundles thus behave as semiflexible rods of thickness <italic>s</italic>, and persistence length <italic>ξ<sub>s</sub></italic> = <italic>π</italic>
<sup>2</sup>
<italic>s</italic>
<sup>4</sup>
<italic>ξ</italic>
<sub>1</sub>, where <italic>ξ</italic>
<sub>1</sub>≈10 <italic>µm</italic> is the persistence length of a single actin filament.</p> <p>The energy change upon linking two bundles is proportional to the length, <italic>L-l</italic>, of their straight and parallel portions (##FIG##5##Fig. 6b##), i.e., Δ<italic>E<sub>fascin</sub></italic> = −(<italic>L</italic>−<italic>l</italic>)<italic>ε</italic>, where <italic>ε</italic> is the gain, per unit length, in the cohesive energy due to the fascin-F-actin bonds. Assuming that all possible actin-fascin contacts are saturated, and that actin filaments are hexagonally packed within a bundle, we find that for two cylindrical bundles which associate into a thicker cylindrical bundle <italic>ε</italic> = 0.2<italic>πsε</italic>
<sub>1</sub>, where <italic>ε</italic>
<sub>1</sub> = 833<italic>kT</italic>/<italic>µm</italic> is the binding energy of fascin to a unit length of a single F-actin. Minimizing Δ<italic>E<sub>dim</sub></italic> = <italic>γb</italic>
<sup>2</sup>
<italic>s</italic>
<sup>2</sup>/<italic>l</italic>
<sup>3</sup>−(<italic>L</italic>−<italic>l</italic>)<italic>ε</italic> with respect to <italic>l</italic>, we find the optimal length of the bent (‘root-like”, ##FIG##5##Fig. 6b##) portion of the bundle: <italic>l<sub>s</sub></italic> = <italic>αb</italic>
<sup>1/2</sup>
<italic>s</italic>
<sup>1/2</sup>, where <italic>α</italic> = (3<italic>γ</italic>/<italic>ε</italic>)<sup>1/4</sup>. Substituting <italic>l<sub>s</sub></italic> into Δ<italic>E</italic>
<sub>dim</sub>(<italic>s</italic>) we get Δ<italic>E<sub>dim</sub></italic>(<italic>s</italic>) = <italic>C<sub>A</sub>s</italic>
<sup>9/4</sup>−<italic>C<sub>B</sub>Ls</italic> where  and <italic>C<sub>B</sub></italic> = 0.2<italic>πε</italic>
<sub>1</sub>.</p> <p>We conclude this analysis by noting that from some point in time onwards, the bundles observed experimentally appear to continue elongating without changing their thickness. Based on our model we may conclude that this happens when the time scale of bundle bending fluctuations which lead to bundle-bundle association becomes long compared to the experimental time scale, or simply long relative to the time it takes to all actin monomers to join the growing bundles.</p>
<title>Discussion</title> <p>In this work we studied in detail the steady state structures formed in a simple <italic>in vitro</italic> system, in which Arp2/3 complex and fascin compete on their binding to actin filaments. Varying mixtures of actin, Arp2/3, and fascin were analyzed in order to elucidate their delicate interplay in determining which of the structures–asters, stars, or networks of bundles–appear in a given concentration regime.</p> <p>We found that in the absence or at very low concentration of fascin the system is dominated by Arp2/3 complex nucleation and branching activities, resulting in the appearance of dense 3D aster-like networks of actin. Increasing fascin concentration induces phase changes, first to stars and then to network of bundles. A star is composed of a dense aster core with actin/fascin bundles emanating radially from its surface. Our experiments show that in the early stages of actin polymerization fascin is passive, while Arp2/3 mediates dense aster-like structures of actin, whose structure is very similar to the one observed in the absence of fascin. Fascin comes into play when actin filaments in the periphery of the aster get long enough, and can thus bend and associate with each other into bundles of parallel filaments, held together by fascin linkers. This is in accordance with <italic>in vivo</italic> experiments that emphasize the importance of filaments length in the creation of filopodia ##REF##15294161##[16]##. The aster-star transition appeared for all Arp2/3 concentrations that we have tested. The second transition, from stars to network of bundles, has only been observed at sufficiently low [Arp2/3]. In conclusion, we find that the system is more sensitive to phase changes when [Arp2/3] is low, and the concentration of fascin necessary to induce structural transformations is lower.</p> <p>The competition between Arp2/3 and fascin is critical in determining the actin structures formed. This is because these two proteins are nucleators of different actin geometries; Arp2/3 complex promotes the formation of branched actin seeds ##REF##10801131##[8]## while fascin initiates unbranched bundle seeds ##UREF##0##[31]##. The ability of fascin to nucleate actin bundles <italic>in vitro</italic>, was recently demonstrated in experiments showing that fascin enhances actin polymerization by promoting the formation of stable disc-like bundle nuclei. These nuclei are composed of short filaments (a few monomers long), which serve as seeds for subsequent actin polymerization ##UREF##0##[31]##. In the presence of Arp2/3, this process competes with branch nucleation of actin. The rates of these two processes are expected to be regulated by the abundance of Arp2/3 complex and fascin.</p> <p>In formulating our model for bundle formation and growth we assumed that only filaments originating at the surface of the aster take part in this process (similar to <italic>in vivo</italic> systems where filopodia originate in the vicinity of the plasma membrane). Filaments from the inner core of the aster cannot significantly contribute to this process, owing to severe excluded volume interactions on the bending of these highly branched, dense, and entangled filaments. We assumed that bundles begin to form when filaments in the aster surface are long enough so that their bending energy penalty is compensated by the gain of actin-fascin linking energy. Since the unbinding time of bundles is significantly longer than the polymerization time of filaments, their association is essentially irreversible, and thus cannot be adequately described based on thermodynamic equilibrium consideration. We have thus proposed a step-wise mechanism of irreversible assembly of single filaments/bundles into thicker bundles, followed by their elongation. The newly formed bundles elongate until their bending energy decreases, and their net association energy into a thicker bundle becomes energetically favorable, etc. This model can also explain the gradual thickening of newly formed bundles <italic>in vivo</italic>, forming Λ-precursors shaped roots ##REF##12566431##[5]##. In our experimental system, the elongation-thickening sequence terminates when all actin monomers are consumed. Based on this picture we could explain why at low [Arp2/3] the bundles emerging from the aster's core are both longer and thicker as compared to those observed in systems with a higher initial [Arp2/3].</p> <p>The evolution of asters to stars is of direct biological relevance because of its similarity to lamellipodia to filopodia transition in cells. Recall that two structural parameters, <italic>L</italic> and <italic>b</italic>, characterize the periphery of aster's network. In our <italic>in vitro</italic> system, the length, <italic>L</italic>, of the (unbranched sections of the) surface filaments increases in time, while the spacing between their origins decreases. In living cells, analogous structural parameters may determine the conditions favoring the formation of filopodia. In the cellular system, filaments length (<italic>L</italic>) and the density of the lamellipodial network (characterized by <italic>b</italic>) is controlled by specific regulatory proteins, such as Arp2/3 complex, VASP and Dia2. In particular, Dia2 was recently shown to be recruited to lamellipodia, where it induces the formation of many long, unbranched filaments (thus locally increasing <italic>L</italic> and reducing <italic>b</italic>), which then gradually converged into filopodial bundles ##REF##18044991##[24]##, ##REF##12086607##[32]##. This behavior is qualitatively consistent with our model, according to which an increasing <italic>L</italic> or decreasing <italic>b</italic>, favors the onset of bundle formation. By recruiting proteins like Dia2 to specific sites along lamellipodia, such a mechanism may enable the cell to control the localization of filopodia along its leading edge.</p> <p>To conclude, using a minimal model system containing only Arp2/3 complex, actin and fascin we were able to mimic complex events which take place at the leading edge of cells. Compared to a cell this system is simple and easily controlled, making it possible to characterize the roles of a limited set of proteins in the higher order assembly of actin filaments. In the future, it will be interesting to measure the size of the actin/fascin bundle using high resolution electron microscopy, and correlate it with the properties of the aster network structure (e.g., density of the actin network, filament length and orientation at the aster periphery). It is worthwhile to test the effects of mDia2 and capping proteins concentrations on bundle formation and dimensions. Finally, quantitative information from <italic>in vivo</italic> studies will be of great value, especially regarding the correlation between filopodia thicknesses and structural properties of the lamellipodial network.</p>
<p>Conceived and designed the experiments: LH ABS AB. Performed the experiments: YI YBK. Analyzed the data: YI YBK. Wrote the paper: YBK ABS AB.</p> <p>During cellular migration, regulated actin assembly takes place at the cell leading edge, with continuous disassembly deeper in the cell interior. Actin polymerization at the plasma membrane results in the extension of cellular protrusions in the form of lamellipodia and filopodia. To understand how cells regulate the transformation of lamellipodia into filopodia, and to determine the major factors that control their transition, we studied actin self-assembly in the presence of Arp2/3 complex, WASp-VCA and fascin, the major proteins participating in the assembly of lamellipodia and filopodia. We show that in the early stages of actin polymerization fascin is passive while Arp2/3 mediates the formation of dense and highly branched aster-like networks of actin. Once filaments in the periphery of an aster get long enough, fascin becomes active, linking the filaments into bundles which emanate radially from the aster's surface, resulting in the formation of star-like structures. We show that the number of bundles nucleated per star, as well as their thickness and length, is controlled by the initial concentration of Arp2/3 complex ([Arp2/3]). Specifically, we tested several values of [Arp2/3] and found that for given initial concentrations of actin and fascin, the number of bundles per star, as well as their length and thickness are larger when [Arp2/3] is lower. Our experimental findings can be interpreted and explained using a theoretical scheme which combines Kinetic Monte Carlo simulations for aster growth, with a simple mechanistic model for bundles' formation and growth. According to this model, bundles emerge from the aster's (sparsely branched) surface layer. Bundles begin to form when the bending energy associated with bringing two filaments into contact is compensated by the energetic gain resulting from their fascin linking energy. As time evolves the initially thin and short bundles elongate, thus reducing their bending energy and allowing them to further associate and create thicker bundles, until all actin monomers are consumed. This process is essentially irreversible on the time scale of actin polymerization. Two structural parameters, <italic>L</italic>, which is proportional to the length of filament tips at the aster periphery and <italic>b</italic>, the spacing between their origins, dictate the onset of bundling; both depending on [Arp2/3]. Cells may use a similar mechanism to regulate filopodia formation along the cell leading edge. Such a mechanism may allow cells to have control over the localization of filopodia by recruiting specific proteins that regulate filaments length (e.g., Dia2) to specific sites along lamellipodia.</p>
<title>Supporting Information</title>
<p>We would like to thank Nir Gov for useful discussions.</p>
<fig id="pone-0003297-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003297.g001</object-id><label>Figure 1</label><caption><title>Phase transition as a function of fascin concentration.</title><p>Conditions are 7 µM G-actin, 6.25 nM Arp2/3, and 12.5 nM GST-VCA. Fascin concentration is: (a) 3 nM, (b) 4 nM, (c) 5 nM, (d) 7 nM, (e) 50 nM, (f) 200 nM, (g) 500 nM, and (h) 3 µM. At very low concentrations of fascin asters are formed. Transition to star-like structures occurs above a fascin concentration of 4 nM (b). The density and the length of the bundles emanating from the star core increase with fascin concentration (b–f). Above a certain fascin concentration of 0.5 µM (g) the size and the number of stars decreases; the stars coexist with an entangled network of actin/fascin bundles, seen in the image background. At 3 µM stars do not form anymore; the system is composed solely of entangled network of actin bundles (h). Bar is 10 µm.</p></caption></fig> <fig id="pone-0003297-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003297.g002</object-id><label>Figure 2</label><caption><title>Transition from aster to star-like to network structures induced by fascin depends on Arp2/3 concentration.</title><p>Conditions are: 7 µM G-actin, upper line (a–d) 12.5 nM Arp2/3 complex, and 25 nM GST-VCA and bottom line (e–h) 100 nM Arp2/3, and 200 nM GST-VCA. Transition from aster to stars and from stars to network: occurs between 3 to 6 nM and between 0.2 to 3 µM at 12.5 nM [Arp2/3] (blue and red arrows, respectively). At 100 nM [Arp2/3] only the transition from aster to star is visible; transition to a network structure is not visible in these fascin concentration ranges. Bars are 10 µm (a–h).</p></caption></fig> <fig id="pone-0003297-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003297.g003</object-id><label>Figure 3</label><caption><title>Phase diagram for actin-fascin-arp2/3 complex.</title><p>Condition is: 7 µM G-actin, Arp2/3 complex and fascin concentration was changed according to the graph. Condition for images is: 3 nM, 200 nM and 3 µM fascin; 6.25 nM, 12.5 nM and 6.25 nM Arp2/3 complex (a, b and c, respectively). Experimental points demonstrate different type of structures formed; these are represented by squares - entangled network, triangles - stars and circles - asters. Phase changes from asters to stars (blue line) and from stars to entangled bundle network (red line). Lines of phase separation are plotted as averages of two points of different structures. For low concentration of Arp2/3 complex phase transition is more abrupt for the same change in fascin concentration. No transition to entangled network was visible for high concentrations of Arp2/3 complex.</p></caption></fig> <fig id="pone-0003297-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003297.g004</object-id><label>Figure 4</label><caption><title>The surface density of objects (asters or stars) is controlled solely by Arp2/3 and not by fascin.</title><p>Conditions are: 7 µM G-actin; fascin concentration of 5 nM (triangles) and 200 nM (circles); and variable amount of Arp2/3 complex: 6.25, 12.5, 25, 40, and 100 nM. The [VCA]/[Arp2/3 complex] = 2 in all experiments. We observe a monotonic growth in density of objects with Arp2/3 complex concentration. The density of objects doesn't show a dependence on fascin concentration (inset, a half log plot is given in order to see clearly all experimental data points). Error bars represent standard deviation from average values.</p></caption></fig> <fig id="pone-0003297-g005" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003297.g005</object-id><label>Figure 5</label><caption><title>Mean star-like bundles' length, <italic>L<sub>Fin</sub></italic>, dependence on the [Arp2/3] and [fascin].</title><p>(a) Conditions are: [G-actin] = 7 µM; [Arp2/3] = 6.25 nM (green dots), 25 nM (blue dots), and 100 nM (red dots). For all three cases, the mean bundles' length, <italic>L<sub>Fin</sub></italic>, increases with [fascin]. (b) Conditions are: 7 µM G-actin; 100 nM fascin and variable amounts of Arp2/3 complex. The data shows that the average of bundles length, <italic>L<sub>Fin</sub></italic>, decrease with [Arp2/3] monotonically (error bars are ±SD).</p></caption></fig> <fig id="pone-0003297-g006" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003297.g006</object-id><label>Figure 6</label><caption><title>Analytic Model and KMC simulations emphasize the role of Arp2/3 in determining the thickness of bundles.</title><p>(a). Left: snapshot of an aster taken from the simulations. Filaments originating in the periphery shell of the aster might bend towards each other and link into bundles. Right: Values of <italic>L<sub>min</sub></italic>, <italic>L<sub>max</sub></italic> and <italic>b</italic> that were taken from the simulations at <italic>t<sub>i</sub></italic>, for the two [Arp2/3] tested. (b). Bundles formation from the branched network occurs by the elongation–association model. The length <italic>l</italic> represents the bent portion of a filament/bundle. (c) The change in energy caused by bending and linking of two adjacent bundles, as a function of their number of shells s (which also corresponds to the bundles' thickness), for several values of <italic>L</italic>. For [Arp2/3] = 12.5 nM: <italic>L</italic> = 6 µm (red), 9 µm (blue), 12 µm (light green), and 15 µm (pink), and <italic>b</italic> = 0.3 µm. For [Arp2/3] = 100 nM (magnified in the black box): <italic>L</italic> = 0.8 µm (light blue), 1.0 µm (brown), and 1.2 µm (green), and <italic>b</italic> = 0.11 µm. In both cases [G-actin] = 7.5 µM.</p></caption></fig>
<inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula>
<supplementary-material content-type="local-data" id="pone.0003297.s001"><label>Appendix S1</label><caption><p>(0.38 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0003297.s002"><label>Figure S1</label><caption><p>Conditions are: 7 µM G-actin, (a–b) 12.5 nM Arp2/3 complex, and 25 nM GST-VCA and (c–d) 100 nM Arp2/3, and 200 nM GST-VCA. The amount of fascin required inducing transition from aster (a and c) to stars (b and d) ranges between 5 to 6 nM and between 6.5 to 7 nM at 12.5 nM and 100 nM Arp2/3 complex, respectively. The insets in b and d show zoom-in images of the bundles emanating from the stars at each [Arp2/3]. At the transition, the bundles emanating from the stars at lower [Arp2/3] (b) are thicker than those originating from the aster core at 100 nM Arp2/3 (d). Bars are 10 µm.</p><p>(2.78 MB TIF)</p></caption></supplementary-material>
<fn-group><fn fn-type="COI-statement"><p><bold>Competing Interests: </bold>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p><bold>Funding: </bold>A.B.G. wishes to thank the Joseph and May Winston Foundation Career Development Chair in Chemical Engineering, the Israel Cancer Association (grant No. 20070020B) and the Israel Science Foundation (grant No. 551/04) for financial support. A.B.S. thanks the US-Israel Binational Science Foundation (BSF grant 2006-401), and the Israel Science Foundation (ISF grant 659/06) for financial support. The Fritz Haber research center, where the calculations were carried out, is supported by the Minerva foundation, Germany.</p></fn></fn-group>
<graphic xlink:href="pone.0003297.g001"/> <graphic xlink:href="pone.0003297.g002"/> <graphic xlink:href="pone.0003297.g003"/> <graphic xlink:href="pone.0003297.g004"/> <inline-graphic xlink:href="pone.0003297.e001.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0003297.e002.jpg" mimetype="image"/> <graphic xlink:href="pone.0003297.g005"/> <graphic xlink:href="pone.0003297.g006"/> <inline-graphic xlink:href="pone.0003297.e003.jpg" mimetype="image"/>
<media xlink:href="pone.0003297.s001.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0003297.s002.tif"><caption><p>Click here for additional data file.</p></caption></media>
[{"label": ["31"], "element-citation": ["\n"], "surname": ["Haviv", "Gov", "Ideses", "Bernheim-Groswasser"], "given-names": ["L", "N", "Y", "A"], "year": ["2007"], "article-title": ["Thickness distribution of actin bundles in vitro."], "source": ["Eur Biophys J "], "comment": ["10.1007/s00249-007-0236-1"]}]
PLoS One. 2008 Sep 29; 3(9):e3297