Elastic Behavior of Cross-Linked and Bundled Actin Networks

See allHide authors and affiliations

Science  28 May 2004:
Vol. 304, Issue 5675, pp. 1301-1305
DOI: 10.1126/science.1095087


Networks of cross-linked and bundled actin filaments are ubiquitous in the cellular cytoskeleton, but their elasticity remains poorly understood. We show that these networks exhibit exceptional elastic behavior that reflects the mechanical properties of individual filaments. There are two distinct regimes of elasticity, one reflectingbendingof single filaments and a second reflectingstretchingof entropic fluctuations of filament length. The mechanical stiffness can vary by several decades with small changes in cross-link concentration, and can increase markedly upon application of external stress. We parameterize the full range of behavior in a state diagram and elucidate its origin with a robust model.

Actin is a ubiquitous protein that plays a critical role in virtually all eukaryotic cells. It is a major constituent of the cytoskeletal network that is an essential mechanical component in a variety of cellular processes, including motility, mechanoprotection, and division (1, 2). The mechanical properties, or elasticity, of these cytoskeletal networks are intimately involved in determining how forces are generated and transmitted in living cells. In vitro, actin can polymerize to form long rigid filaments (F-actin), with a diameter D ≈ 7 nm and contour lengths up to L ≈ 20 μm (3). Despite their rigidity, thermal effects nevertheless modify the dynamic configuration of actin filaments, leading to bending fluctuations. The length scale where these fluctuations completely change the direction of the filament is the persistence length, Embedded Image, determined when the bending energy equals the thermal energy; here, κo is the bending modulus of the actin filament. Entangled solutions of F-actin are a well-studied model for semiflexible polymers; they exhibit exceptionally large elasticity at small volume fractions as compared with more traditional flexible polymers at similar concentrations (46). However, in vivo, the cytoskeletal network is regulated and controlled not only by the concentration of F-actin, but also by accessory proteins that bind to F-actin. Indeed, nature provides a host of actin-binding proteins (ABPs) with a wide range of functions: They can align actin filaments into bundles and they can cross-link filaments or bundles into networks (2, 7). However, unlike flexible polymers, where changes in cross-link density do not markedly affect the elasticity (8), small changes in the concentration of cross-linker ABPs dramatically alter the elasticity of F-actin networks (46, 911). A quantitative understanding of this behavior is essential to control and exploit biomimetic materials based on actin or other semiflexible polymer networks. It is also an essential first step in understanding the elasticity of cytoskeletal networks and the pivotal role of ABPs in the mechanical behavior of the cell. However, very little is known about the microscopic mechanisms underlying the elasticity of semiflexible polymer networks and, in particular, the quantitative behavior of cross-linked and bundled F-actin networks.

We study the elasticity of a model composite network of bundled and cross-linked F-actin. We use an ABP that simultaneously bundles and cross-links actin filaments while being sufficiently rigid so as not to contribute to the network compliance. The elastic behavior of such in vitro networks exhibits marked variation as the F-actin and cross-link densities are varied. For example, by changing the cross-link concentration alone, the network elastic modulus can be reliably and precisely tuned over more than three decades (Fig. 1A); this is in marked contrast to conventional flexible polymer materials, where the elasticity is rather insensitive to chemical cross-links (8). In addition, the network can exhibit considerable stiffening with applied stress, with the effective elasticity increasing by more than a decade with virtually no change in strain; thus, the compliance nearly vanishes. By contrast, at lower concentrations of either F-actin or cross-links, the nonlinear strain stiffening completely disappears, and the network elasticity can remain linear for strains as large as unity. Thus, there are two qualitatively distinct regimes of elasticity, highlighted by the different-colored regions of the state diagram in Fig. 1A. A robust model elucidates the fundamental mechanisms of the network elasticity in each regime and rationalizes both the linear and nonlinear behavior of the elasticity captured by the state diagram (1214).

Fig. 1.

(A) Schematic representation of the state diagram of cross-linked and bundled F-actin networks, showing the variations in the elasticity of the networks with changes in actin concentration, cA, or density of cross-linking R = cS/cA, where cS is the scruin concentration. By varying either cA or R, the elastic modulus, Go, can be varied by more than three orders of magnitude, from less than 0.1 Pa to 300 Pa. For large values of cA or R (indicated by the red plane), the network elasticity becomes nonlinear at very small deformations, and the elasticity becomes strain dependent and increases markedly upon increased strain. By contrast, at very low values of cA or R (indicated by the blue plane), the network mechanical response remains linear up to large strains, and no strain stiffening is observed. Transitions between these two regimes of elasticity can be made by changing either R or cA, as indicated by the dashed lines. (B to D) Confocal microscope images of actin networks, cA = 11.9 μM, labeled with Texas Red–phalloidin as a function of R. (B) R = 0, (C) R = 0.07, (D) R = 0.5. As R is increased for a fixed cA, actin filaments are cross-linked into tight bundles. The average bundle thickness is weakly proportional to R between R = 0.03 and R = 1. Because the total filament concentration remains constant, the average mesh size also increases. Bar (B), 3 μm. The bright circles in (B) to (D) are 1-μm-diameter colloidal particles.

The actin cytoskeleton is a highly dynamic and complex network, in which ABPs play a myriad of roles (2, 7). Most ABPs implicated in cross-linking or bundling the network are themselves highly dynamic and compliant, making it extremely difficult to isolate the origins of elasticity of the F-actin network itself (5, 9, 10, 15). To overcome these difficulties, and define a benchmark for the elasticity of cross-linked actin networks, we use scruin, which is found in the sperm cell of the horseshoe crab (16). Scruin-calmodulin dimers decorate individual F-actin filaments, and nonspecific scruin-scruin interactions cross-link and bundle neighboring filaments (17). In vivo, scruin mediates the formation of a single, crystalline bundle of 80 actin filaments; this functions as a mechanochemical spring in the acrosomal process (16, 18). In vitro, actin filaments polymerized in the presence of scruin are cross-linked and bundled by scruin-scruin contacts and form an isotropic, disordered three-dimensional network that is macroscopically homogeneous over a large range of cross-linking concentrations (19). Scruin bonds are irreversible, and scruin itself is not compliant; thus, the elasticity of our networks appears to be due entirely to the F-actin filaments.

In the absence of scruin, an F-actin solution forms an isotropic and homogeneous mesh; at a concentration of 11.9 μM, the mean separation between filaments is ∼300 nm (Fig. 1B). (20). We tune the degree of filament cross-linking and bundling by varying the scruin concentration, cS, for a fixed actin concentration, cA, such that the density of cross-links, R, is cS/cA and varies from 0 to 1; a pelleting assay confirms that the degree of polymerized actin remains constant for a given cA as R is varied (19). The morphology of actin-scruin networks remains similar to that of an entangled F-actin solution with values of R up to 0.03. For R > 0.03, we observe the formation of compact bundles (Fig. 1C); the bundle thickness increases even further as R is increased, as shown for R = 0.5 in Fig. 1D. Using electron microscopy, we measure the average bundle diameter, DB, as a function of R; for R > 0.03, DB varies approximately as DBRx, where x ≃ 0.3 (19). For R = 1, the average bundle thickness is DB ≈ 50 nm, which leads to a significant increase in the rigidity of the actin bundles. The bending rigidity of any rod increases rapidly with its width, varying as D4 (21), and for scruin-mediated actin bundles, κB ≃ κo(DB/D)4, consistent with theory (17). Thus, the average bundle stiffness increases to ∼600κo, greatly increasing the persistence length. We can also directly control the average distance between filaments, ξ, by varying both R and cA (22), Embedded Image(1) Thus, we can finely tune the network morphology.

To measure the elasticity of these networks, we apply a force and measure the resultant deformation; the elastic modulus, G′, is defined as the ratio of the stress, σ, or force per unit area, to the strain, γ; thus, G′ = σ/γ. The networks can also have a viscous or dissipative response, and the resultant stress depends on the strain rate, Embedded Image, defining the loss modulus, Embedded Image. More generally, G′(ω) and G″(ω) depend on the frequency and are measured by applying a small oscillatory stress at a frequency, ω. In these cross-linked networks, the elastic modulus dominates the mechanical response, reaching a frequency-independent plateau, Go, at frequencies less than 1 Hz (fig. S1). The magnitude of Go is highly dependent on the morphology of the network; for a small amount of cross-linking, the composite network is a soft elastic gel, and Go decreases weakly as R decreases below 0.03 (Fig. 2A). However, for R > 0.03, Go increases dramatically with R. Thus, the elasticity is strongly dependent on the scruin-mediated interactions that cross-link and bundle actin filaments, even at the same cA.

Fig. 2.

The elastic modulus, Go, and maximum strain, γmax, as a function of R and cA. (A) Go as a function of R for cA = 11.9 μM. The solid line indicates a scaling of GoR2. (B) γmax versus R for cA = 11.9 μM. The solid line indicates γmaxR–0.6. (Inset) Typical strain stiffening response for a composite network with cA = 11.9 μM and R = 0.03, indicating the onset of nonlinear elastic response, γcrit, and the maximum strain, γmax, before the network breaks. (C) Go as a function of cA for R = 0.03 (◼), 0.13 (△), and 0.3 (⚫). The solid line shows a scaling of Embedded Image and the error bars indicate our sample-to-sample reproducibility. (D) γmax as a function of c for the same values of R. The solid line shows a scaling of Embedded Image.

The magnitude of Go is also highly dependent on actin concentration for a fixed value of R. For R = 0.12, we measure Embedded Image, where α ≈ 2.5 as shown by the open triangles in Fig. 2C, which is consistent with all proposed theoretical network models to within our experimental uncertainty (12, 23, 24). We see similar scaling for all values of R, as shown in Fig. 2C for weakly cross-linked networks, R = 0.03 (squares), and thickly bundled networks, R = 0.3 (circles). Thus, the elastic modulus is also strongly dependent on the filament density.

The elasticity of the networks can also exhibit a marked dependence on applied strain, γ. In the regime where Go increases rapidly with R (Fig. 1A), the elastic modulus is linear below a critical strain, γcrit ∼ 10%. However, above this strain, Go increases dramatically over a very small change in γ up to a maximum strain, γmax, whereupon the network breaks irreversibly (Fig. 2B, inset). This strain stiffening is completely reversible and the network can be cycled through strains up to γmax without any history dependence of the elastic response. These composite networks constitute an exceptional material, whose modulus increases significantly with virtually no change in γ; thus, they have a nearly vanishing compliance for γ>γcrit.

The origin of the elasticity of entangled F-actin solutions is purely entropic, even though the filaments are semiflexible, with a rather large bending rigidity (25, 26). Thus, the elasticity arises because an applied strain reduces the accessible fluctuations; this is an entropic elasticity, similar to that of flexible polymers such as rubber. However, when filaments are chemically cross-linked into networks such that the distance between cross-links, Embedded Imagec, is less than Embedded Imagep, or when the persistence length is increased as a result of bundling, it is unclear whether the network elasticity originates from the actual bending of filaments, which is an enthalpic elasticity, or from stretching out thermal fluctuations of the filaments, which is an entropic elasticity (12, 23, 24).

To delineate the role of entropic effects in network elasticity, we model the network as a collection of thermally fluctuating, semiflexible filament segments of length, Embedded Image, the average distance between cross-links, such that Embedded Image. Thermally driven transverse fluctuations reduce the end-to-end filament extension and lead to an entropic spring; for small extensions, the force, F, required to extend this spring by length Embedded Image is Embedded Image, where kBT is the thermal energy (12) (fig. S2); this is the wormlike chain model that also describes the behavior of DNA (27). Accounting for the concentration of chain segments through the mesh size, we can calculate the resultant stress (8), σ = F2, and thereby determine the linear elastic modulus of the network by applying a strain, Embedded Image. We find (12) (fig. S3) Embedded Image(2) We assume that, for a fixed R, Embedded Imagec is directly proportional to the distance between filament entanglements (12, 28) (fig. S3). Thus, Goc11/5A, which depends strongly on actin concentration, consistent with that observed for both weakly cross-linked networks, R = 0.03, and for thickly bundled networks, R = 0.3. We conclude that the elastic stiffness of both weakly cross-linked and thickly bundled actin networks is entropic in origin, originating from the stretching out of thermal fluctuations of individual actin filaments (Fig. 3).

Fig. 3.

A summary of the interpretation of results presented in this paper. (A) Networks are polymerized between two plates of area, A, and separation, h. By varying the filament or cross-link concentration, we vary the network microstructure; low filament and cross-link density are shown in the left cartoon, whereas high filament and cross-link density are shown in the right cartoon. (B) We measure the mechanical properties of the networks by applying a force, F, per unit area, A, or stress and measuring the deformation or strain γ ≡ x/h. The microscopic distribution of the strain differs in the networks. Our data show two distinct mechanical regimes, distinguished by their nonlinear response. One of these is consistent with macroscopic mechanical properties due to enthalpic filament bending, as expected for low cross-link or filament density. This leads to a highly inhomogeneous distribution of strain, as indicated by the red arrows. By contrast, the mechanical properties of dense networks are consistent with stretching of thermally induced filament fluctuations and an entropic elasticity. The resultant strain is uniform throughout the sample, as indicated by the red arrows. (C) The table summarizes the essential differences in the elasticity of entropic and enthalpic networks.

A critical test of this picture for the origin of the elasticity is the strain dependence; when an entropic, semiflexible network is extended, thermally induced transverse fluctuations are pulled out. At sufficiently large extension, the response ceases to be linear, consistent with our observed strain stiffening above γcrit. Entropic, semiflexible networks break at an applied strain that is expected to decrease weakly with increasing filament concentration (12) (fig. S3); this is also consistent with our observations, as shown in Fig. 2D. We can quantitatively test our model for the elasticity by determining the divergence of the stress response in the nonlinear elastic regime. Here, the force required to extend a single semiflexible filament, Embedded Image, diverges dramatically as full extension is approached, where Embedded Image (27, 29). However, the extreme nonlinearity of the modulus makes precise oscillatory strain measurements virtually impossible, because the measured stress waveforms become highly nonsinusoidal. To overcome this limitation, we superpose a small oscillatory stress of magnitude δσ, for a constant applied stress, σo, with δσ/σo < 0.1 at 0.1 Hz, and determine the differential elastic modulus, K′(σo) = [δσ/δγ]|σo as a function of σo. This measurement applies a constant prestress while measuring the response to a small oscillatory stress, allowing us to accurately probe the divergence of the response as a function of applied prestress.

When the elastic modulus is strain-independent, the differential modulus is the same as the elastic modulus, K′(σo) = Go. However, above the critical stress, σcrit = Goγcrit, K′ increases markedly until the network breaks, as shown for a large range of actin concentrations at fixed R in Fig. 4. In this stress stiffening regime, K′(σo) ∼ σ 3/2o, consistent with that predicted for a single semiflexible polymer, dF/dEmbedded ImageF3/2 (27, 29). The theoretical form of K′(σo) for these networks is predicted by using the full force-extension relation of a single actin filament and taking into account variations in Embedded Imagec and cA, as well as spatial averaging of the orientation of chain segments (19). By scaling the theory to fit the data for one concentration, cA = 29.4 μM and R = 0.03, we determine Embedded Imagec and calibrate the absolute stress; the theory is then in excellent accord with the remaining data sets, with no additional adjustable parameters, as shown by the curves in Fig. 4. Moreover, the theory suggests that the functional form of the data should be the same for all values of cA and R. Consistent with this prediction, all the data can by scaled onto a single master curve (Fig. 4, inset). This provides conclusive support for our hypothesis that the origin of both the linear and the nonlinear elasticity of the network is due to stretching of thermal fluctuations of single filaments. Thus, the elastic stiffening of the network results from the nonlinear force-extension behavior of individual filaments; this further supports our hypothesis that the elasticity of cross-linked actin networks is entropic in origin.

Fig. 4.

The differential elastic modulus, K′, as a function of applied steady shear stress, σo, for R = 0.03 and several values of cA; the lines through the data indicate theoretical predictions for each concentration determined from the single-filament response. The values of cA are 29.4 μM(▲—), 21.4 μM (⚫ - ·), 11.9 μM (◼ - -), and 8.33 μM (♢ - ··). The stress stiffening response of a bundled network, R = 0.5, at cA = 7 μM is also shown (◯). The dotted line at the right indicates σo3/2o. (Inset) The data sets rescaled by σcrit and Go showing the universal form of the stress stiffening response; the rescaled theory is indicated by the solid line.

The origin of the network elasticity appears to remains entropic, even as κ increases as a result of bundling, or as Embedded Imagec decreases as a result of increased cross-linking. We demonstrate this by tuning these two parameters simultaneously, by varying R while keeping cA fixed. The resultant networks exhibit variations in both Go and γmax with R that are consistent with our model (fig. S3). They also exhibit similar strain stiffening, and the qualitative form of the divergence of the incremental modulus remains unchanged. For a given cA, the modulus increases as GoR2, and the critical strain simultaneously decreases as γmaxR–0.6 (Fig. 2). Interestingly, two networks can have the same Go for different actin concentrations, provided that R varies; however, the onset of nonlinear behavior occurs at a different σcrit, as shown by the open circles in Fig. 4. Thus, changes in the network morphology caused by variation in the filament cross-linking and bundle thickness dramatically affect both the linear and nonlinear elastic behavior.

At very low values of cA or R, the qualitative features of the elastic response change dramatically. This is most pronounced at very low values of R, where Go exhibits virtually no dependence on R (Fig. 2A). Furthermore, these networks do not exhibit any strain stiffening but instead deform easily upon increased applied stress, and the elastic response remains linear for strains as large as γ∼ 1. This suggests that there is a second regime of elasticity with a qualitatively different origin. To capture the full variety of elastic behavior, we summarize all the data in the cA-R state diagram (Fig. 5). The elastic modulus varies continuously over four orders of magnitude as both cA and R are varied. We delineate the two distinct regimes with different symbols. In the first regime, for sufficiently large cA or R, Go is highly sensitive to both cA and R and exhibits pronounced strain stiffening (denoted by the plus signs in Fig. 5). The origin of the mechanical response is quantitatively explained by our model of entropic elasticity, which reflects the stretching and compression of individual filaments. This model implicitly assumes that the deformations of the network are affine, or distributed uniformly throughout the sample, and thus the strain is homogeneous at all length scales (Fig. 3). There is a sharp transition, denoted by the dashed line in Fig. 5, to the second regime of elasticity, where Go is much less sensitive to both cA and R and exhibits no strain stiffening (denoted by the open circles in Fig. 5). This behavior is consistent with a model that suggests that the network deformation at lower filament or cross-link concentrations may be nonaffine; the elastic connections of the network are sparse, and its response is dominated by bending of individual filaments in isolated locations, leading to a highly inhomogeneous strain of the network (13, 14) (Fig. 3). In this model, the linear modulus is predicted to be Goc2A (23, 24). Furthermore, because network elasticity is determined by filament bending, which is rather small on the molecular scale, there is no mechanism for strain stiffening. Thus, as summarized in Fig. 5, a network of rigidly cross-linked semiflexible polymers can exhibit a rich variety of elastic behaviors simply by varying the filament concentration, cross-link density, or bundle thickness.

Fig. 5.

The R-cA state diagram of actin:scruin networks detailing the tunability of Go; colors indicate a range from 0.03 Pa (purple) to 300 Pa (red), as shown by the legend. Symbols differentiate networks that stress stiffen (+) from those that do not (◯).

We speculate that there may also be a third regime of elasticity, in the limit of very high values of cA and R; here, the response will again be affine but will be totally enthalpic in nature, and the compliance will be dominated by the mechanical stretching and compression of filaments or bundles, with no universal mechanism for strain stiffening. Our results also have important implications for how a cell may regulate its mechanical response, with small variations in either filament or cross-link concentration to make rapid and precise transitions to qualitatively alter the mechanical strength of the cytoskeletal network. Furthermore, the mechanical properties of the cytoskeletal actin network may also be varied as a function of external prestress, providing a potential mechanism for mechanosensing and mechanoprotection in the cell by allowing small variations in force to significantly increase its mechanical strength. Finally, our data highlight the importance of understanding single-molecule mechanics to interpret the mechanical properties of networks. Because we can so precisely parameterize the elastic behavior of cross-linked actin networks, it is now possible to be quantitatively predictive for in vivo physiological conditions. Other physiological cross-linkers, unlike scruin, are neither rigid nor inextensible; as a result, the single-molecule elasticity and dynamics of the cross-linker must also be included when modeling these more complex networks. However, our data using the rigid cross-linker scruin provide an essential benchmark for the elasticity of the cytoskeletal actin network.

Supporting Online Material

Materials and Methods

Figs. S1 to S3


References and Notes

View Abstract

Stay Connected to Science

Navigate This Article