Side-binding proteins modulate actin filament dynamics
Abstract
Actin filament dynamics govern many key physiological processes from cell motility to tissue morphogenesis. A central feature of actin dynamics is the capacity of filaments to polymerize and depolymerize at their ends in response to cellular conditions. It is currently thought that filament kinetics can be described by a single rate constant for each end. In this study, using direct visualization of single actin filament elongation, we show that actin polymerization kinetics at both filament ends are strongly influenced by the binding of proteins to the lateral filament surface. We also show that the pointed-end has a non-elongating state that dominates the observed filament kinetic asymmetry. Estimates of flexibility as well as effects on fragmentation and growth suggest that the observed kinetic diversity arises from structural alteration. Tuning elongation kinetics by exploiting the malleability of the filament structure may be a ubiquitous mechanism to generate a rich variety of cellular actin dynamics.
https://doi.org/10.7554/eLife.04599.001eLife digest
Actin is one of the most abundant proteins in cells. It forms networks of filaments that provide structural support and generate the forces needed for cell movement, division, and many other processes in cells.
Filaments of actin continuously change in length as actin molecules are added or removed at the ends. One end of an actin filament—called the barbed-end—grows much faster than the other, known as the pointed-end. Many other proteins also help the actin filaments to form. Some of these proteins bind to the ends of the filaments, where they directly control the growth of the filaments. Other proteins bind along the length of the filaments, but how these ‘side-binding’ proteins influence the growth of filaments is not clear.
In this study, Crevenna et al. used a technique called ‘total internal reflection fluorescence (TIRF) microscopy’ to study how several side-binding proteins affect the growth of actin filaments in an artificial system. The growth of the barbed-ends was strongly influenced by which side-binding protein was interacting with the filament. For example, the barbed-end grew rapidly when a protein called VASP was present but grew more slowly in the presence of the protein α-actinin. Although the growth at the pointed-end was generally slow and sporadic, the side-binding proteins also had noticeable effects.
Crevenna et al. found that when the side-binding proteins were present at low levels, filament growth was similar for all proteins studied. It was only when the proteins were present at higher levels that the growth of the actin filaments was altered depending on the specific side-binding protein present. One side-binding protein called α-actinin also altered the shape of the actin filament so that when it was present at high levels, the filaments curved in a particular direction. Together, these results suggest that the growth, structure, and flexibility of actin filaments can be strongly influenced by the various proteins that bind along the length of the filaments.
The next challenges are to understand the precise details of how these side-binding proteins are able to alter the growth and shape of actin and investigate how they influence other processes that control the structure of actin networks in cells.
https://doi.org/10.7554/eLife.04599.002Introduction
Central cellular processes such as cell migration, cytokinesis, endocytosis, and mechanosensation depend critically on actin-based force generation and actin filament turnover (Pollard and Borisy, 2003; Lecuit et al., 2011). The molecular basis of actin filament turnover derives from the association and dissociation of monomers from each filament end and depends on the nucleotide (ATP, ADP · Pi, or ADP) bound to the actin monomer (Pollard, 1986). The filament is kinetically asymmetric, where one end (called the barbed-end) is observed to grow an order of magnitude faster than the other end (the pointed-end) (Pollard, 1986). In addition, the critical concentration for polymerization is different for the two ends. The origin of the asymmetry is not fully understood. Measurements of filament elongation as a function of solution viscosity (Drenckhahn and Pollard, 1986) and particle-analysis from cryo-electron microscopy (Narita et al., 2011) suggest the existence of a non-elongating state at the pointed-end. Although growth pauses have been previously observed during filament elongation measured using total internal reflection fluorescence (TIRF) microscopy (Kuhn and Pollard, 2005; Fujiwara et al., 2007), these pauses were attributed to artifacts and were not characterized further. The dynamics of the pointed-end plays an important role in both the origin of the differences in critical concentration observed at the two ends in the presence of ATP (Pollard, 1986; Fujiwara et al., 2007); and in filament treadmilling, where, barbed-end growth and pointed-end shrinking occur simultaneously (Bugyi and Carlier, 2010). Thus, we have focused on performing an accurate and detailed analysis of both barbed-end and pointed-end dynamics using TIRF microscopy.
In cells, a large number of proteins interact with actin filaments, either at the ends or with the lattice. End-binding proteins regulate actin dynamics by limiting elongation (at the barbed-end) or serving as anchor points (for the pointed-end). Side-binding proteins, on the other hand, are much more diverse encompassing myosin motors, cross-linkers or bundlers as well as severing proteins. The interaction of the actin filament with a particular subset of proteins defines the molecular composition, architecture, and overall turnover of sub-cellular arrays such as stress fibers and filopodia. Some of these arrays are tightly packed (Jasnin et al., 2013) and dynamics of the filaments are influenced by the local environment. The mechanisms of how some proteins are recruited to these structures while others are excluded are a subject of intense research (Cai et al., 2008; Hansen et al., 2013). Although the overall filament dynamics have been thought to be sensitive to the concentration of the side-binding protein (Breitsprecher et al., 2009), it is not understood how and to what extent side-binding proteins alter filament kinetics, structure, and flexibility.
In this study, we used TIRF microscopy to study the effect of side-binding proteins on the dynamics of actin filament growth in vitro. We chose three cross-linking proteins and one motor protein to represent the large variety of interacting proteins and used them to tether filaments directly to the surface of a glass slide for visualization. We used the chemically inactivated myosin II motor protein (NEM-myosin) as it is the standard choice for this type of assay (Kuhn and Pollard, 2005). The filamin protein (Kueh et al., 2008) was used, which is an important player in cellular mechanosensing that is evolutionary-conserved (Razinia et al., 2012), as its use as a tether has recently generated some debate (Mullins, 2012; Niedermayer et al., 2012). Additionally, we selected α-actinin, a molecule that, together with myosin II, forms stress fibers (Langanger et al., 1986), and VASP, a protein that localizes to areas of dynamic actin reorganization such as filopodia and the lamelipodium (Rottner et al., 1999). By carrying out these assays with several proteins that bind to the side of actin filaments, we were able to explore the possible range of modulation available to actin filament dynamics and delineate intrinsic filament properties.
Results
Kinetic modulation at the barbed-end
Fluorescently labeled actin was used to visualize the growth of actin filaments (Figure 1A–B) using TIRF microscopy. In this technique, single actin filaments are tethered to a glass surface via a side-binding protein and their growth and/or shrinkage is monitored in real time (Figure 1A–B). From each frame, the filament is extracted and a kymograph is constructed (Figure 1—figure supplement 1). The position of each end of the filament was then determined by fitting an error function (Demchouk et al., 2011) to each line of the kymograph (see ‘Materials and methods’ and Figure 1—figure supplement 1 for details). This end-detection method provides a more accurate determination of the filament length and thereby a more reliable estimate of the instantaneous elongation velocity compared to methodologies used previously (Figure 1—figure supplement 1).
The single-filament elongation experiments showed that the barbed-end grew at a constant velocity with occasional pauses for all constructs measured while the barbed-end elongation velocity varied depending on the particular side-binding protein used (Figure 1C,D). The elongation velocity ‘E’ at the barbed-end was the fastest with VASP and the slowest with α-actinin (Figure 1C). By varying the free actin concentration from 0.3 to 2 μM, we estimated the barbed-end association and dissociation rates, kon and koff, respectively (Figure 1D), using only the periods of elongation (i.e., E > 1.5 sub·s−1, referred to hereafter as the ‘kinetically active’ phases). Compared to the previously reported value of 11.6 sub·μM−1·s−1 for actin only in the absence of tethering proteins (Pollard, 1986), we found a higher value of kon in the presence of VASP, a similar value for actin alone and for NEM-myosin, and lower values when α-actinin or filamin were used (Figure 1D and Table 1). Extrapolating the elongation velocity as a function of actin concentration to zero actin provides an estimate of the dissociation rate, koff, of ATP-actin at the barbed-end (Table 1). In the presence of filamin, koff is indistinguishable from zero, whereas in the presence of VASP, koff increased compared to the value in the presence of NEM-myosin (1.6 ± 0.5 s−1). The estimated koff we measured in the presence of NEM-myosin was in agreement with the previously reported value of 1.4 s−1 (Pollard, 1986), whereas in the presence of α-actinin, koff was lower than 1.4 s−1. The ratio of inferred dissociation rates to the calculated association rate (i.e., koff/kon) is the critical concentration at which polymerization will occur and has been estimated to be ∼150 nM for the barbed-end (Pollard, 1986). We find a similar value (∼0.2 μM) for filaments elongated in the presence of VASP, α-actinin and NEM-myosin, but close to zero for filamin. The use of VASP induced the largest change in the measured kinetics. In contrast to the other three proteins measured, the kinetics were enhanced. VASP is known to act as a polymerase at the barbed-end by delivering subunits to the growing end (Hansen and Mullins, 2010; Breitsprecher et al., 2011). It achieves this function through its two actin-binding domains: a G-actin-binding domain that delivers monomers (via G-actin binding or GAB domain) while remaining attached to the filament via the F-actin-binding or FAB domain (Breitsprecher et al., 2008, 2011; Hansen and Mullins, 2010). To rule out the polymerase activity of VASP as the cause of enhanced kinetics, we also tested a VASP construct that lacks the GAB domain but retains its capacity to interact with the filament lattice (Breitsprecher et al., 2008). We continued to observe fast polymerization using this VASP-ΔGAB protein as a tether (Figure 1) in agreement with previous reports (Breitsprecher et al., 2008). The measured kinetic rates were 70 ± 30 sub·μM−1·s−1 for kon and 14 ± 9 s−1 for koff, about half of those determined using the full-length VASP construct. Therefore, the effect of immobilized VASP on filament kinetics is not only due to recruitment of monomers to the growing filament. These results demonstrate that ATP-actin kinetics at the barbed-end are sensitive to the particular side-binding protein interacting with the filament.
Kinetic modulation at the pointed-end
Pointed-end association and dissociation rates were estimated in the same manner as those for the barbed-end (Figure 2A). Both the estimated association rates and dissociation rates varied according to the associated side-binding protein used as a tether (Figure 2A and Table 1). The presence of filamin increased the by a factor of ∼5 (from 0.8 in the presence of NEM-myosin to 2.8 sub·μM−1·s−1). The for α-actinin was 0.9 sub·μM−1·s−1, while, when using VASP or VASP ΔGAB, the rate was 44 sub·μM−1·s−1 and 16 sub·μM−1·s−1, respectively. On the other hand, the presence of filamin also increased the inferred by almost an order of magnitude from 0.4 (in the presence of NEM-myosin) to 2.6 s−1. The inferred rates were 0.7 s−1, 8 s−1, and 5 s−1 with α-actinin, VASP, and VASP ΔGAB, respectively (Table 1).
Unlike the barbed-end where there were occasional pauses (Figure 1C), the pointed-end displayed mostly a kinetically inactive phase or paused state and only grew sporadically (Figure 2B,C). Such kinetically inactive phases were observed for all free actin concentrations tested (250 nM–2 μM). Above the pointed-end critical concentration (e.g., using a free actin concentration of 1 μM), we observed a discontinuous (i.e., growth-pause) behavior for all side-binding proteins (Figure 2B). In the presence of VASP or filamin, pointed-end elongation was readily observed. Pointed-end elongation was much more difficult to visualize when using NEM-myosin and α-actinin (Figure 2B) where elongation occurred for brief periods of time and with slower rates. The elongation velocity during kinetically active phases was influenced strongly by the different tethering proteins used (Figure 2A). Elongation velocity followed the order of VASP > VASP ΔGAB > filamin > α-actinin > NEM-myosin (Figure 2A,B). On the other hand, at 300 nM free actin monomer concentration, i.e., below the pointed-end critical concentration of ∼600 nM (Pollard, 1986), we observed barbed-end growth (Figure 1D) and pointed-end depolymerization (Figure 2C), i.e., treadmilling, in the presence of filamin as a tethering protein (Figure 2C). Treadmilling was also present using NEM-myosin and α-actinin, albeit with slower rates, since pointed-end depolymerization establishes the overall treadmilling rate. In contrast to our expectations, there was no shrinkage at the pointed-end below the critical concentration but polymerization in the presence of VASP or VASP ΔGAB (Figure 2). These results suggest that side-binding proteins can also determine actin filament pointed-end growth and depolymerization dynamics. Additionally, these results show that observed effects at one end do not necessarily represent effects at both ends. For example, filamin reduces only the dissociation rate (and therefore the critical concentration) at the barbed-end although it alters both the association and dissociation rate at the pointed-end.
The elongation rate varies with occupancy of the side-binding proteins
Next, we studied how sensitive filament dynamics are to the presence of each of the proteins tested. Therefore, we measured the elongation rates and pausing as a function of the side-binding protein surface density (Figure 3). For this, we varied the total protein concentration that was allowed to adsorb to the glass surface, therefore changing the number of tethering proteins that interact with a single filament. We estimated the lattice-binding protein surface density from the protein concentration, the sample volume (∼10 μL) and the surface to which the sample was adsorbed (a flow cell of 5 mm × 20 mm, giving 100 mm2) as done previously (Howard et al., 1989; Crevenna et al., 2008). All protein in solution was assumed to adsorb on the upper and lower glass surfaces. To achieve consecutive lower tether densities, the total protein concentration was serially diluted. At low tethering protein concentrations, individual filaments swiveled around distinctive attachment points indicating that they are bound to single tethering molecules as observed previously (Howard et al., 1989; Crevenna et al., 2008). To estimate the density in an alternative manner, we measured the average number of pivot points per micron of filament at the two lowest protein concentrations and divided that by the average area covered during swiveling. Assuming a linear scaling with protein concentration, this estimate results in a lower density (by a factor of 2) compared to those reported in Figure 3 and throughout the text. Estimated densities ranged from ∼5 up to ∼18,000 molecules·μm−2, which are equivalent to values between 0.1 and 110 tethers per micron of filament (Figure 3).
At low tether densities (5–200 μm−2 or 0.025 to 1.0 molecules per micron of filament), the dynamics were independent of the tethering protein used. As an example, filamin is shown in Figure 3A–D. Elongating actin filaments (at a free actin monomer concentration of 1 μM) showed mostly kinetically active phases (Figure 3A), and the elongation velocity distributions were centered around ∼9 subunits·s−1 (Figure 3 and Figure 3—figure supplement 1). At high tether densities (600–18,000 μm−2 or 3–110 molecules per micron of filament), each side-binding protein tested generated a particular elongation behavior (Figure 3E and Figure 3—figure supplement 1). Using filamin, increasing the surface tether density decreased the mean elongation velocity of kinetically active phases (Figure 3B,D) and increased the fraction of time the filament spent in a paused state, i.e., the pausing probability ‘Pp’ (Figure 3B,D,F). In contrast, increasing the VASP or the VASP-ΔGAB density increased the elongation velocity while VASP also increased the Pp (Figure 3 and Figure 3—figure supplement 1). Higher surface concentrations of α-actinin or NEM-myosin had also an effect on the elongation velocity (Figure 3—figure supplement 1) and, in addition, the density of NEM-myosin had a strong effect on the Pp (Figure 3F and Figure 3—figure supplement 1).
One possible explanation for these results could be geometric and/or mechanical constrains imposed on the filament by the high density of the surface-immobilized side-binding protein used. To investigate this possibility, we carried out experiments where a very low density of NEM-myosin was used to tether filaments to the surface while a second, side-binding protein was present in solution. We tested the effects of VASP, VASP-ΔGAB, and filamin on filament growth. The influence of all three side-binding proteins on the elongation rate of both the barbed- and pointed-end were similar to what we observed when using them to immobilize the actin filaments to the surface (Figure 3—figure supplement 2 and Table 1). We also performed these experiments with human VASP and the elongation rates measured for both Dd VASP and Hs VASP agree with previously reported results (Breitsprecher et al., 2008; Hansen and Mullins, 2010) (Figure 3—figure supplement 3). These results suggest that a variety of elongation kinetics can arise from the specific interaction of actin filaments with the particular associated side-binding protein.
Intrinsic filament dynamics
To further verify that the observed kinetic changes and pauses originate from the particular side-binding protein used as a tether, we investigated the intrinsic properties of filament elongation and controlled for artifacts. Single elongating filaments were measured at the lowest protein surface density possible that still allowed filament visualization. At the lowest α-actinin tether density used (5 molecules μm−2, which corresponds to 1 tether molecule every 5–10 microns along the filament), the ends swiveled around their tethering site due to Brownian motion and were clearly free of the surface (Figure 4A). Under these conditions, barbed-ends showed continuous elongation (Figure 4B) while pointed-ends were typically in the kinetically inactive state (only 2 of 50 filaments showed growth or depolymerization, Figure 4C–D). For the other tethering proteins, only the paused state was observed on freely swiveling pointed-ends (Figure 4—figure supplement 1). Using only the pause-free elongation velocities for each actin concentration tested, we estimated association and dissociation rates (slopes in Figure 4E, Table 1). Our estimated values for the pause-free elongation kinetics agree well with those previously obtained by EM (Table 1), which were measured on the 20–60 s time scale. When we convolute our pointed-end pause-free elongation rate with the pausing probability, our results are comparable to the kinetics estimated by TIRF experiments (Kuhn and Pollard, 2005). One possible explanation for this discrepancy is that continuous pointed-end growth occurs at the beginning of filament assembly, which is suggested from our data (Figure 4C) and is the time scale on which the EM data was acquired (Pollard, 1986). Moreover, the pausing probability, Pp, at either end was insensitive to the actin concentration used (Figure 4F).
The low density used for these experiments and the observed pauses on freely swiveling actin filaments (pointed-end only) rules out surface effects (Kuhn and Pollard, 2005) as the determining cause for the pauses at the ends. Another possible source of pauses is light-induced photo-dimerization. From the work of Niedermayer et al. (2012), it is possible to quantitatively predict the accumulated fraction of filaments where depolymerization has been paused as a consequence of exposure to light (Figure 4—figure supplement 2). In contrast with this prediction, we observed all swiveling filament pointed-ends, under depolymerizing conditions, to be in a kinetically inactive state at the beginning of image acquisition (N = 40, Figure 4—figure supplement 2). Only in the presence of a medium to high density of tethering proteins did we observe depolyermization of pointed-ends (12 of 55, Figure 4—figure supplement 2).
As an additional test to rule out any tether, surface or light-induced effect of the pausing, we used a two-color solution assay to investigate pointed-end growth. Here, a small seed (formed with atto565-labeled actin) was allowed to grow in solution for 15 min in the presence of atto488-labeled monomers, followed by stabilization, dilution, and visualization of the filaments (Figure 4—figure supplement 3). At a free actin concentration of 1 μM, the concentration used in solution to allow filament elongation, all pointed-ends are expected to grow at an average rate of ∼0.5 sub/s (Pollard, 1986). In contrast to this expectation, we observed that only 20% of the seeds grew at the pointed-end (N = 1000, Figure 4—figure supplement 3). This percentage is higher than we observe in the surface-based experiments, which could be due to annealing of filaments in solution (Sept et al., 1999; Andrianantoandro et al., 2001) or due to lack of the tethering protein. What is clear is that the non-elongating or paused state is not due to either surface or light-induced effects. Taken together, these results show that a single rate constant describes filament elongation kinetics from ATP-monomers in the absence of side-binding proteins and that the pointed-end has an intrinsic kinetically inactive state.
Structural effects of side-binding proteins on filaments
During the course of filament elongation analysis as a function of side-binding protein density on the surface (Figure 3), we noticed that filaments appeared more bent as the tether density increased. To quantify this curviness, we estimated an apparent persistence length ‘’ of individual filaments associated with different side-binding proteins (see ‘Materials and methods’ for details). The persistence length Lp (Boal, 2012) reflects the material properties of the filament, which are related to its structure (Chu and Voth, 2005, 2006; Pfaendtner et al., 2010), and has already been shown to be tunable by side-binding proteins (such as myosin or cofilin [McCullough et al., 2008; Murrell and Gardel, 2012; Bengtsson et al., 2013]). At the lowest side-binding protein density (∼10 molecules μm−2 or ∼0.1 molecules per filament micron), actin filaments had an of ∼18 μm and was independent of the associated protein (Figure 5A). At the highest densities (∼16,000 molecules·μm−2 or ∼100 molecules per filament micron), the presence of NEM-myosin decreased the to 4 ± 1 μm while it was reduced to 3 ± 1 μm, 5 ± 1 μm, and 2.2 ± 0.3 μm when using filamin, VASP and α-actinin, respectively (N > 50 for each condition, Figure 5A). Estimates for the persistence length of surface adsorbed filaments are consistent with what has been determined for freely fluctuating filaments (McCullough et al., 2008; Graham et al., 2014). We also tested the mechanical effect of the side binding proteins on the actin filament when the actin filaments were attached to the surface with a low density of NEM-myosin and the side-binding protein was present in solution. Again, we observed a decrease in the persistence length of about 30% for filamin, α-actinin, and NEM-myosin, whereas the effect was about 50% in the presence of VASP or VASP ΔGAB (Figure 5A).
Two other interesting phenomena were observed in the presence of side-binding proteins at high densities. First, the presence of filamin increased the spontaneous fragmentation of filaments (20 out of 197 filaments vs less than 1 fragmentation even per 200 filaments) (Figure 5B). Second, barbed-end elongation when tethered with α-actinin had a preference to grow in a counterclockwise direction (Figure 5C). This counterclockwise elongation observed when α-actinin was present was independent of the length the filament had when it landed on the surface. These observations suggest an influence of the tethering protein on the structural properties of the filament.
Changes in the structure of the filament by binding proteins are known to be able to propagate over several subunits (Orlova et al., 1995). We hypothesize that structural alteration might be the origin of kinetic modulation. To test how well this interpretation could explain our results, we constructed a simple model based on long-range structural alterations to describe elongation as a function of tether density (Figure 3E). In our model, the interaction between a growing filament and a tethering protein gives rise to a modified association rate (i.e., ) at the site of interaction, which then propagates over a certain distance LC (Figure 3—figure supplement 4). From the interaction site, there is a linear decay of the modified association rate, as a linear decay would be expected for the release of torsional stress. Using a Monte Carlo method (see ‘Materials and methods’ for details), we calculated α and LC (which are the only free parameters) by comparing the simulated elongation rates to the experimental values and minimizing the χ2 (Table 2). This simple model satisfactorily described our experimental results (Figure 3E, solid lines) and provides an estimate for the propagation length. We also considered the possibility that VASP acts via a ‘local increase in monomer concentration’ similar to Breitsprecher et al. (2011) (see ‘Materials and methods’ for details). This local concentration model did not account for the tether density dependence of elongation velocity (Figure 3—figure supplement 5), nor it can account for the effect observed with VASP ΔGAB that does not have the capacity to bind actin monomers.
Our estimate of the propagation length should be considered as a lower limit since additional factors that could potentially influence our assay such as tether unbinding (i.e., which would lower the effective number of interacting side-binding proteins) and/or alternative tether density calculations (see above) would result in longer propagation lengths. The induced effects are local-to-short-ranged (∼2-11 monomers) for α-actinin, filamin, and NEM-myosin while they are long-ranged when using VASP and VASP-ΔGAB (160 and 76 monomers, respectively) (Table 2). Collectively, these results suggest that the side-binding proteins tested alter the structure of the filament.
Discussion
We have shown that asymmetry in filament elongation is a consequence of a non-elongating state at the pointed-end and that the general versatility of actin dynamics may be a response to the binding of various proteins. Through accurate measurements of pointed-end association kinetics, we have observed that experiments performed in the absence of tethering proteins and in the presence of VASP or VASP ΔGAB yielded equivalent critical concentrations for both ends (∼0.2 μM, Table 1). This implies that, under these conditions, detailed balance at equilibrium is fulfilled, i.e., (Hill, 1987). We propose that the existing discrepancy in estimated critical concentrations at both ends (Fujiwara et al., 2007) originates from the presence of a previously uncharacterized kinetically inactive or non-elongating state at the pointed-end. This kinetically inactive state is consistent with a non-elongating structural conformation observed by cryo-electron microscopy (Narita et al., 2011). The kinetic asymmetry of the pure actin filament may be low (∼3) and such non-elongating or closed conformation at the pointed-end would reinforce the effective filament asymmetry. The transition at the pointed-end from the open to the closed state may be coupled to ATP hydrolysis or phosphate release at the terminal subunit. The presence of this open-to-closed transition at the pointed-end would explain why the terminal subunit has an estimated different rate of phosphate release compared to the filament lattice (Fujiwara et al., 2007).
Actin filaments in association with any of the five proteins tested displayed a change in elongation velocity, an increase in pausing, and a change in filament flexibility. Therefore, it is possible that these three characteristics have a common origin. For three of these filament-binding proteins (myosin, α-actinin, and filamin), the binding interface to the actin filament is formed by two consecutive monomers along the same strand (Galkin et al., 2008, 2010; Lorenz and Holmes, 2010). These side-binding proteins might directly occlude the binding site for the next monomer either partially (reducing the elongation velocity) or completely (giving rise to an elongation pause). Partial distortion of the filament could turn into a defect that propagates along the lattice decreasing the observed filament stiffness and impacting the association rate. In this respect, side-binding proteins could be thought of as allosteric regulators of actin filament kinetics. Indeed, actin filaments are known to be subject to allosteric regulation by other associated proteins (Egelman and Orlova, 1995; Galkin et al., 2012). In particular, myosin (Prochniewicz and Thomas, 1997), cofilin (Galkin et al., 2001; Prochniewicz et al., 2005), dystrophin (Prochniewicz et al., 2009), and utrophin (Prochniewicz et al., 2009) are known to induce structural changes in the actin filament. Similar to filamin and α-actinin, dystrophin and utrophin bind actin through calponin-homology (CH) domains (Galkin et al., 2010). Moreover, binding to the filament is cooperative for cofilin (De La Cruz, 2005), αE-catenin (Hansen et al., 2013), and myosin (Orlova and Egelman, 1997). The basis for this allosteric regulation could originate from the stabilization of an existing structural state of the filament (Galkin et al., 2001), given that the actin filament is structurally polymorphic (Galkin et al., 2010). Therefore, it is possible that the observed elongation kinetics and pauses arise from direct modulation of the filament structure. In line with this hypothesis, two other proteins, the actin-binding domain of αE-catenin (Hansen et al., 2013) and an N-WASP construct (Khanduja and Kuhn, 2014), have recently been shown to alter filament kinetics and one of them, the actin-binding domain of αE-catenin, also influences filament structure (Hansen et al., 2013). Although atomically accurate simulations and more high resolution experiments are required to understand the molecular basis of monomer association and dissociation from the filament ends, our results provide evidence that lattice structural changes affect actin filament growth kinetics. The influence of different side-binding proteins on the growth kinetics was found to persist over different length scales. Although we do not currently know the mechanism of this difference, it is interesting to note that actin-binding proteins with globular binding domains (α-actinin, filamin, and NEM-myosin) have short-range affects whereas VASP and VASP ΔGAB, which have an unstructured binding motif, have more long-range affects.
Our experimental approach of using tethers immobilized on a solid surface imposes geometric and/or mechanical constrains on filament growth. As actin filaments form part of the cell cortex (Biro et al., 2013) and focal adhesions (Kanchanawong et al., 2010) where they assemble into oligomeric membrane-anchored complexes with many actin-binding proteins tethered to the plasma membrane surface, our studies may not be too far from the biologically relevant situation in living cells. Moreover, the cell interior is very crowded (Luby-Phelps, 2000) and some sub-cellular actin arrays are tightly packed (Jasnin et al., 2013). Both of these conditions may lead to the immobilization of actin-binding proteins and generate similar constrains during filament growth. In addition, the presence of side-binding proteins in solution is sufficient for altering the filament kinetics and mechanics. Depending on the local cross-linker protein abundance in the cell, turnover kinetics on the order of 1 μm of filament within ∼1 min can be achieved, a rate at which treadmilling could become a contributing factor to cellular retrograde flow in the lamellipodium (Watanabe and Mitchison, 2002; Ponti et al., 2004). Additionally, filament structural changes generated by side-binding proteins may also play a more active role in the identity and turnover of actin-based sub-cellular structures than previously thought, by regulating processes such as branching and fragmentation (Hansen et al., 2013) or network mechanics (Jensen et al., 2014). Given the vast number of side-binding proteins, kinetic modulation via structural alteration may be a general regulatory mechanism of actin dynamics.
Materials and methods
Proteins
Actin was obtained from chicken muscle using the method of acetone powder. Actin was extracted by one round of polymerization and pelleting by centrifugation (Spudich and Watt, 1971). The resulting pellet was depolymerized in G-buffer (1 mM Tris–HCl pH 7.8, 2 mM ATP, 2 mM CaCl2, 2 mM DTT) overnight at 4°C followed by gel filtration on a Sephacryl S-300 column. Myosin was purified and chemically inactivated with N-Ethyl-Maleimide according to the published protocol (Breitsprecher et al., 2009). Atto488-actin, α-actinin, and filamin were purchased from Hypermol (Bielefeld, Germany). Alternatively, actin was labeled with succinimidyl ester atto488 (ATTO-TEC GmbH, Germany) on random lysine residues. Actin labeling was performed under polymerization conditions (50 mM KCl and 2 mM MgCl2) followed by depolymerization and gel filtration in G-buffer. The functionality of 1:1 dye:protein lysine-modified actin was found to be unaffected by the labeling as has been previously characterized using pyrene polymerization assays, TIRF elongation, EM and FCS experiments (Crevenna et al., 2013). Unlabeled and labeled actin were mixed to yield a final ratio of 2:1 unlabeled:labeled actin molecules. The actin mixture (20 μl) was snap frozen and stored at −80°C until further use. Before use, an actin aliquot was centrifuged to remove possible aggregates.
A plasmid containing the gene of Dd VASP was kindly provided by J Faix, (Hanover, Germany). VASP was expressed using a pCoofy plasmid in Sf9 cells with a MBP-tag and purified following standard methods as described previously (Scholz et al., 2013). MBP-VASP was used without cleavage, since removal of the tag resulted in protein aggregation and degradation. For VASP ΔGAB purification, the VASP coding sequence without residues 198–220 was amplified using the pCoofy28-full-length VASP as a template, forward primer 5′-GCGCTTTTATCAACACCGCCACCTGCGGCTGG-3′ and the reverse primer 5′-GCAGGTGGCGGTGTTGATAAAAGCGCTGGTGTACCAACAAAAAC-3′. Then the VASP(delta198-220) coding sequence was further cloned into a pEC-GST vector and expressed using E. coli BL21(DE3) as reported previously (Wang et al., 2014). Briefly, the E. coli strain was grown at 37°C in 2 L of ZY auto-induction medium for 5 hr and then the temperature was reduced to 18°C overnight. Cells were harvested and resuspended in 50 mM Tris pH7.5, 500 mM NaCl, 1 mM DTT supplemented with protease inhibitors and the cells were disrupted using sonication. The protein was purified from clarified cell lysate using a 5-ml GSTrap FF column (GE Healthcare, Germany) with elution buffer 50 mM Tris, pH7.5, 150 mM NaCl, 20 mM reduced glutathione, 1 mM DTT and further purified using size exclusion chromatography (Superdex 75, GE Healthcare, Germany) with buffer 50 mM Tris, pH7.5, 500 mM NaCl, 1 mM DTT.
Imaging
Request a detailed protocolFlow cells were made as a sandwich of a cover slip (20 × 20 mm), parafilm with an approximate 5-mm wide channel and a glass slide. The surfaces of the flow cells were passivated to avoid adsorption of actin to the sample holder by incubating them with 10% (wt/vol) of BSA in PBS for 10 min. Flow cells were washed three times with 90 μl of G-buffer. The tethering protein was then applied for 5 min and the flow cell was then washed again three times with 90 μl of G-buffer. Actin (33% atto488-actin) was incubated 5 min on ice with 1/10 volume of 10x ME buffer (400 μM MgCl2 and 2 mM EGTA) to exchange Ca2+ for Mg2+. The actin-containing solution was mixed with imaging buffer (catalase, β-mercaptoethanol, glucose oxidase, 0.8% [vol/vol] D-glucose, 0.25% [wt/vol] methylcellulose, and 1/10 volume of 10x KMEI buffer [500 mM KCl, 20 mM MgCl2, 20 mM EGTA, and 300 mM imidazole], with a final pH of 7.1) and introduced into the flow cell. TIRF microscopy was performed using a TILL photonics inverted microscope (FEI Munich GmbH, Germany). A single actin aliquot was used within 12 hr.
The lattice-binding protein surface density was estimated from the protein concentration, the sample volume (∼10 μl) and the surface to which the sample was adsorbed (a flow cell of 5 mm × 20 mm, giving 100 mm2) as done previously (Howard et al., 1989; Crevenna et al., 2008). All protein in solution was assumed to adsorb on the upper and lower glass surfaces. To achieve consecutive lower tether densities, the total protein concentration was serially diluted. At low tethering protein concentrations, individual filaments swiveled around distinctive attachment points indicating that they are bound to single tethering molecules as observed previously (Howard et al., 1989; Crevenna et al., 2008). To estimate the density in an alternative manner, we measured the average number of pivot points per micron of filament at the two lowest protein concentrations and divided that by the average area covered during swiveling. Assuming a linear scaling with protein concentration, this estimation resulted in a slightly lower density (by a factor of 2) compared to those reported in Figure 3. The concentration-based estimated densities represent an upper limit and are easy to reproduce. Hence, we report both in the text.
Two-color solution assay
Request a detailed protocolFilaments were formed using atto565-labeled actin in G-buffer by addition of 1/10 volume of 10x KMEI buffer. After more than 2 hr of polymerization at room temperature, filaments were fragmented by shearing and subsequently mixed with atto488-labeled monomers and allowed to elongate for 15 min. Filaments were then stabilized with unlabeled phalloidin and diluted for imaging on an Epi-Fluorescent Microscope (Axiovert 200, Zeiss, Germany).
Data analysis
Request a detailed protocolRaw movies were corrected for x- and y- stage drift by first calculating its magnitude via image correlation spectroscopy (Hebert et al., 2005), and secondly, correcting the drift by bicubic interpolation. Drift estimation and correction were implemented in custom programs written in LabView and MATLAB (The MathWorks, MA). Kymographs of single filaments were made using Metamorph or Image J, while further analysis was carried out using MATLAB. Filament analysis tools are available at: http://www.cup.uni-muenchen.de/pc/lamb/actin_filament_dynamics.html. The position of the filament tip, per line in the kymograph, was estimated by fitting an error function as previously described (Demchouk et al., 2011). More than 20 filaments were analyzed per condition. To estimate the first pause distribution, we used the model described by Niedermayer et al. (2012) with ω = 2 × 106. The light intensity for treadmilling experiments ranged from 0.74 to 0.92 mW·mm−2. Growth orientation was assessed manually with the following criteria: Barbed-end filament growth direction was classified as straight/not-defined, clockwise or counterclockwise from experiments at the highest surface tether density.
Monte Carlo simulations
Request a detailed protocolFor the model presented in Figure 3, we used a Monte Carlo method to simulate the polymerization of actin filaments at the barbed-end. For each condition, a 105 monomer long actin filament was polymerized, and the instantaneous elongation rate was calculated for every point and then averaged over the total length of the filament. The average elongation rate was calculated from the length of the polymer over time for each condition. The effect of the tethering protein was simulated by a change in the effective kon at the site and vicinity of the tethering protein (Figure 3—figure supplement 4). Over the 105 monomer sites, Nb side-binding proteins were randomly placed corresponding to the desired side-binding protein density. This gives on the order of 104 side binding proteins for the highest densities computed. First, the tethering protein's positions were randomly chosen according to the protein density. The surface density was calculated assuming that all added tethering proteins adsorbed to the surface and were functionally active. To convert from surface density to fractional occupancy, we used the area occupied by 1 μm of actin filament (0.006 μm2), using a value of 370 subunits per micron of filament. We neglected tether dissociation from the filament, as this would only reduce the effective tether density. During the simulated polymerization, the effective was changed to at the position of the tethering protein and decreased linearly until reaching the free actin value of after a characteristic length (LC) counted in monomers of actin. The values of 11 μM-1s−1 and 2 s−1 for barbed-end and koff, respectively, were taken from literature (Pollard, 1986), and the average elongation rate was taken from measurements at the lowest tethering protein density (Figure 3E). A first round of simulations was performed to roughly estimate the optimal interval for the parameters (α and LC). A second set of simulations over restrained intervals, with a better resolution on α and LC, yielded a 3D space (α, LC, v). The 2D χ2 was calculated, the minimum value gave the best (α, LC) and the confidence intervals were taken by Δχ2 = 1 (68%) and Δχ2 = 4 (95%). The only free parameters were α and LC, which were determined for each curve (i.e., elongation as a function of lattice-binding protein density) by comparing the simulated elongation rates to the experimental values and minimizing the χ2 using:
where v(xi) is the experimental velocity for density i, vsim(xi, α, LC) is the corresponding simulated velocity with parameters α and LC, and std(xi) is the experimental standard deviation for this data point. Hence, the only free parameters are α and LC and these were determined for each curve (i.e., elongation as a function of lattice-binding protein density). All simulations were done in MATLAB.
Persistence length calculation
Request a detailed protocolIndividual filaments were extracted from the measured data using the algorithm of ‘open active contours’ within JFilamin, a plug in for Image J (Li et al., 2009; Smith et al., 2010). The filament persistence length ‘Lp’ was determined by calculating the angular correlation (Isambert et al., 1995):
where the brackets represent the average correlation function of the tangent θ, measured along the contour length s. The point spacing used to reconstruct a single filament was between 6 and 10 points per micron to avoid artifacts in the Lp estimation (Isambert et al., 1995; McCullough et al., 2008; Smith et al., 2010). All data analysis was done in MATLAB.
Local change in actin concentration model
Request a detailed protocolThe strong increase in the elongation rate in the presence of, for example, VASP as a tethering protein was first thought to originate from the multiple actin monomer binding sites on each VASP protein (Breitsprecher et al., 2008, Hansen et al., 2010, Breitsprecher et al., 2011). Theoretically, the polymerization kinetics is expected to be inhomogeneous along the actin filament and only to be locally enhanced through a higher local concentration of free actin monomers due to the presence of VASP. The growth at one end of the polymer can be written as follows:
where E(x) is the elongation rate at the position x along the filament axis, c(x) the local concentration of globular actin at this position, and kon and koff are the association and dissociation rates. The average elongation rate along the filament length is given by,
where 〈c〉=c0(1+4d), c0 is the free actin concentration in solution and d is the density of VASP protein at the surface. Given that one VASP protein can bind 4 actin monomers, the local concentration of actin monomers available for polymerization can be significantly increased at the site of a tethering protein. The relationship between the average elongation rate and the protein surface density was expected to be linear as follows:
Experimental points do not show a linear dependency on the protein surface density as expected from this model (Figure 3—figure supplement 5). An alternative mode of operation for VASP has been recently postulated where the protein not only increases the local concentration but also transfers monomers from its monomer binding domains to the filament tip (Breitsprecher et al., 2011). The surface density dependency of this alternative model would, nonetheless, predict a linear behavior as well, albeit with a different slope. In addition, this increased local concentration model would not explain the effect observed when using filamin, α-actinin, or with the VASP ΔGAB construct, which is unable to bind monomers.
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
-
Actomyosin interactions and different structural states of actin filamentsBiophysical Journal 104:480a–481a.https://doi.org/10.1016/j.bpj.2012.11.2655
-
Analysis of actin assembly by in vitro TIRF microscopyMethods in molecular biology 571:401–415.https://doi.org/10.1007/978-1-60761-198-1_27
-
Molecular mechanism of Ena/VASP-mediated actin-filament elongationThe EMBO Journal 30:456–467.https://doi.org/10.1038/emboj.2010.348
-
Control of actin filament treadmilling in cell motilityAnnual Review of Biophysics 39:449–470.https://doi.org/10.1146/annurev-biophys-051309-103849
-
Allostery of actin filaments: molecular dynamics simulations and coarse-grained analysisProceedings of the National Academy of Sciences of the United States of America 102:13111–13116.https://doi.org/10.1073/pnas.0503732102
-
Coarse-grained modeling of the actin filament derived from atomistic-scale simulationsBiophysical Journal 90:1572–1582.https://doi.org/10.1529/biophysj.105.073924
-
Electrostatics control actin filament nucleation and elongation kineticsThe Journal of Biological Chemistry 288:12102–12113.https://doi.org/10.1074/jbc.M113.456327
-
Cofilin binding to muscle and non-muscle actin filaments: isoform-dependent cooperative interactionsJournal of Molecular Biology 346:557–564.https://doi.org/10.1016/j.jmb.2004.11.065
-
Microtubule tip tracking and tip structures at the nanometer scale using digital fluorescence microscopyCellular and molecular bioengineering 4:192–204.https://doi.org/10.1007/s12195-010-0155-6
-
Elongation of actin filaments is a diffusion-limited reaction at the barbed end and is accelerated by inert macromoleculesThe Journal of Biological Chemistry 261:12754–12758.
-
Allostery, cooperativity, and different structural states in F-actinJournal of Structural Biology 115:159–162.https://doi.org/10.1006/jsbi.1995.1040
-
Polymerization kinetics of ADP- and ADP-Pi-actin determined by fluorescence microscopyProceedings of the National Academy of Sciences of the United States of America 104:8827–8832.https://doi.org/10.1073/pnas.0702510104
-
High-resolution cryo-EM structure of the F-actin-fimbrin/plastin ABD2 complexProceedings of the National Academy of Sciences of the United States of America 105:1494–1498.https://doi.org/10.1073/pnas.0708667105
-
Actin filaments as tension sensorsCurrent Biology 22:R96–R101.https://doi.org/10.1016/j.cub.2011.12.010
-
Actin depolymerizing factor stabilizes an existing state of F-actin and can change the tilt of F-actin subunitsThe Journal of Cell Biology 153:75–86.https://doi.org/10.1083/jcb.153.1.75
-
Opening of tandem calponin homology domains regulates their affinity for F-actinNature Structural & Molecular Biology 17:614–616.https://doi.org/10.1038/nsmb.1789
-
Structural polymorphism in F-actinNature Structural & Molecular Biology 17:1318–1323.https://doi.org/10.1038/nsmb.1930
-
αE-catenin actin-binding domain alters actin filament conformation and regulates binding of nucleation and disassembly factorsMolecular Biology of the Cell 24:3710–3720.https://doi.org/10.1091/mbc.E13-07-0388
-
VASP is a processive actin polymerase that requires monomeric actin for barbed end associationThe Journal of Cell Biology 191:571–584.https://doi.org/10.1083/jcb.201003014
-
Flexibility of actin filaments derived from thermal fluctuations. Effect of bound nucleotide, phalloidin, and muscle regulatory proteinsThe Journal of Biological Chemistry 270:11437–11444.https://doi.org/10.1074/jbc.270.19.11437
-
Three-dimensional architecture of actin filaments in Listeria monocytogenes comet tailsProceedings of the National Academy of Sciences of the United States of America 110:20521–20526.https://doi.org/10.1073/pnas.1320155110
-
Mechanism of calponin stabilization of cross-linked actin networksBiophysical Journal 106:793–800.https://doi.org/10.1016/j.bpj.2013.12.042
-
Processive acceleration of actin barbed-end assembly by N-WASPMolecular Biology of the Cell 25:55–65.https://doi.org/10.1091/mbc.E12-11-0781
-
Dynamic stabilization of actin filamentsProceedings of the National Academy of Sciences of the United States of America 105:16531–16536.https://doi.org/10.1073/pnas.0807394105
-
The molecular organization of myosin in stress fibers of cultured cellsThe Journal of Cell Biology 102:200–209.https://doi.org/10.1083/jcb.102.1.200
-
Force generation, transmission, and integration during cell and tissue morphogenesisAnnual Review of Cell and Developmental Biology 27:157–184.https://doi.org/10.1146/annurev-cellbio-100109-104027
-
Automated actin filament segmentation, tracking and tip elongation measurements based on open active contour modelsProceedings/IEEE International Symposium on Biomedical Imaging 2009:1302–1305.https://doi.org/10.1109/ISBI.2009.5193303
-
The actin-myosin interfaceProceedings of the National Academy of Sciences of the United States of America 107:12529–12534.https://doi.org/10.1073/pnas.1003604107
-
Cytoarchitecture and physical properties of cytoplasm: volume, viscosity, diffusion, intracellular surface areaInternational review of cytology 192:189–221.https://doi.org/10.1016/S0074-7696(08)60527-6
-
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
-
The instability of stabilizationProceedings of the National Academy of Sciences of the United States of America 109:10743–10744.https://doi.org/10.1073/pnas.1207980109
-
F-actin buckling coordinates contractility and severing in a biomimetic actomyosin cortexProceedings of the National Academy of Sciences of the United States of America 109:20820–20825.https://doi.org/10.1073/pnas.1214753109
-
Structural basis for the slow dynamics of the actin filament pointed endThe EMBO Journal 30:1230–1237.https://doi.org/10.1038/emboj.2011.48
-
Intermittent depolymerization of actin filaments is caused by photo-induced dimerization of actin protomersProceedings of the National Academy of Sciences of the United States of America 109:10769–10774.https://doi.org/10.1073/pnas.1121381109
-
Cooperative rigor binding of myosin to actin is a function of F-actin structureJournal of Molecular Biology 265:469–474.https://doi.org/10.1006/jmbi.1996.0761
-
Structural dynamics of F-actin: II. Cooperativity in structural transitionsJournal of Molecular Biology 245:598–607.https://doi.org/10.1006/jmbi.1994.0049
-
Structure and dynamics of the actin filamentJournal of Molecular Biology 396:252–263.https://doi.org/10.1016/j.jmb.2009.11.034
-
Rate constants for the reactions of ATP- and ADP-actin with the ends of actin filamentsThe Journal of Cell Biology 103:2747–2754.https://doi.org/10.1083/jcb.103.6.2747
-
Dystrophin and utrophin have distinct effects on the structural dynamics of actinProceedings of the National Academy of Sciences of the United States of America 106:7822–7827.https://doi.org/10.1073/pnas.0812007106
-
Cofilin increases the torsional flexibility and dynamics of actin filamentsJournal of Molecular Biology 353:990–1000.https://doi.org/10.1016/j.jmb.2005.09.021
-
Filamins in mechanosensing and signalingAnnual Review of Biophysics 41:227–246.https://doi.org/10.1146/annurev-biophys-050511-102252
-
VASP dynamics during lamellipodia protrusionNature Cell Biology 1:321–322.https://doi.org/10.1038/13040
-
The regulation of rabbit skeletal muscle contraction. I. Biochemical studies of the interaction of the tropomyosin-troponin complex with actin and the proteolytic fragments of myosinThe Journal of Biological Chemistry 246:4866–4871.
-
Structural basis for the extended CAP-Gly domains of p150(glued) binding to microtubules and the implication for tubulin dynamicsProceedings of the National Academy of Sciences of the United States of America 111:11347–11352.https://doi.org/10.1073/pnas.1403135111
Article and author information
Author details
Funding
Deutsche Forschungsgemeinschaft (SPP 1464)
- Roland Wedlich-Söldner
- Don C Lamb
Deutsche Forschungsgemeinschaft (SFB 1035)
- Oliver F Lange
Consejo Nacional de Ciencia y Tecnología (CONACYT DAAD- A/09/7253)
- Marcelino Arciniega
Narodowe Centrum Nauki (2011/01/B/NZ1/00031)
- Kaja Kowalska
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank D Kovar for the Image J macros's to perform the filament analysis, J Faix for the DdVASP expression plasmid and J Dominguez for valuable discussions. This work was supported by the Max Planck Society and by grants from the Deutsche Forschungsgemeinschaft through the SFB 1035 (to OFL), the SPP 1464 (to RWS and DCL), the Excellent Clusters Nanosystems Initative Munich (NIM) and Center for Integrated Protein Science Munich (CIPSM), the Cells-in-Motion Cluster of Excellence (EXC1003 – CiM) of the University of Münster, a CONACYT-DAAD scholarship A/09/7253 (to MA) and the Ludwig-Maximilians-Universität München (via the LMUInnovativ BioImaging Network (BIN) and the Center for NanoScience (CeNS)). This research has also been supported (to KK) by the 2011/01/B/NZ1/00031 grant from the National Science Centre, Poland.
Copyright
© 2015, Crevenna 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
-
- 3,841
- views
-
- 772
- downloads
-
- 25
- 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
-
- Computational and Systems Biology
- Physics of Living Systems
Explaining biodiversity is a fundamental issue in ecology. A long-standing puzzle lies in the paradox of the plankton: many species of plankton feeding on a limited variety of resources coexist, apparently flouting the competitive exclusion principle (CEP), which holds that the number of predator (consumer) species cannot exceed that of the resources at a steady state. Here, we present a mechanistic model and demonstrate that intraspecific interference among the consumers enables a plethora of consumer species to coexist at constant population densities with only one or a handful of resource species. This facilitated biodiversity is resistant to stochasticity, either with the stochastic simulation algorithm or individual-based modeling. Our model naturally explains the classical experiments that invalidate the CEP, quantitatively illustrates the universal S-shaped pattern of the rank-abundance curves across a wide range of ecological communities, and can be broadly used to resolve the mystery of biodiversity in many natural ecosystems.
-
- Computational and Systems Biology
- Physics of Living Systems
Planar cell polarity (PCP) – tissue-scale alignment of the direction of asymmetric localization of proteins at the cell-cell interface – is essential for embryonic development and physiological functions. Abnormalities in PCP can result in developmental imperfections, including neural tube closure defects and misaligned hair follicles. Decoding the mechanisms responsible for PCP establishment and maintenance remains a fundamental open question. While the roles of various molecules – broadly classified into “global” and “local” modules – have been well-studied, their necessity and sufficiency in explaining PCP and connecting their perturbations to experimentally observed patterns have not been examined. Here, we develop a minimal model that captures the proposed features of PCP establishment – a global tissue-level gradient and local asymmetric distribution of protein complexes. The proposed model suggests that while polarity can emerge without a gradient, the gradient not only acts as a global cue but also increases the robustness of PCP against stochastic perturbations. We also recapitulated and quantified the experimentally observed features of swirling patterns and domineering non-autonomy, using only three free model parameters - the rate of protein binding to membrane, the concentration of PCP proteins, and the gradient steepness. We explain how self-stabilizing asymmetric protein localizations in the presence of tissue-level gradient can lead to robust PCP patterns and reveal minimal design principles for a polarized system.