Requirements for contractility in disordered cytoskeletal bundles

Actomyosin contractility is essential for biological force generation, and is well understood in highly organized structures such as striated muscle. Additionally, actomyosin bundles devoid of this organization are known to contract both in vivo and in vitro, which cannot be described by standard muscle models. To narrow down the search for possible contraction mechanisms in these systems, we investigate their microscopic symmetries. We show that contractile behavior requires non-identical motors that generate large enough forces to probe the nonlinear elastic behavior of F-actin. This suggests a role for filament buckling in the contraction of these bundles, consistent with recent experimental results on reconstituted actomyosin bundles.


I. INTRODUCTION
The large-scale motion of living organisms often depends on their ability to harness the power of nanometersized molecular motors to generate macroscopic displacements. For instance, our striated muscles rely on clusters-or "thick filaments"-of the molecular motor myosin to generate forces. Myosin thick filaments are able to slide directionally towards the barbed end of polar actin filaments (F-actin), and the characteristic organization of striated muscles into periodic sarcomeres arranged in series allows the transfer of this microscopic motion to larger scales [ Fig. 1(a)] [1].
Despite its familiarity, sarcomere-like organization is far from a universal feature of contractile actomyosin assemblies. In some instances, partially periodic arrangements reminiscent of sarcomeres are observed, as in subcellular contractile bundles known as stress fibers [2]. In many other cases, however, no such organization is known to exist. Examples include smooth muscle fibers [3], transverse arcs [4], graded polarity bundles [5], the cell cortex [6] and lamellar networks [7]. Sarcomere-like contraction is unlikely to apply to these systems, and there is no consensus regarding their actual deformation mechanism.
In vitro experiments using purified proteins are useful for understanding contraction in these systems, and have been used to identify the minimum requirements of actomyosin contractility since the 1940s [8]. Modern attempts using dilute actomyosin gels were not able to induce contractility in the presence of actin and myosin * martinlenz@uchicago.edu alone, although adding the actin cross-linker α-actinin did produce observable contraction [9][10][11][12]. However, a more recent study using denser actomyosin bundles shows that F-actin and myosin can induce contractility on their own [13]. Unlike in sarcomeres, in these bundles F-actin lacks polarity ordering and myosins are not  aligned in register, and we thus refer to them as "disordered". Previous theoretical work on disordered actomyosin systems include continuum models focused on length scales much larger than an individual actin filament [14][15][16][17]. In these elegant descriptions, contractility is introduced phenomenologically, which circumvents the question of its emergence from microscopic interactions. Several other studies do however investigate this connection. In simulations without sarcomeric organization or crosslinkers, thick filaments simply sort F-actin by polarity without inducing any overall contraction [18]. To restore contractility, several models assume that thick filaments tend to dwell at the barbed end of F-actin after sliding over its whole length [19][20][21][22], or more generally that their velocity depends on their position relative to the filament [23]. In these models, F-actin tend to have immobilized motors that transiently act as passive cross-linkers at their barbed ends. This essentially introduces a small amount of sarcomere-like organization and results in contractility [18]. However, no direct experimental evidence of thick filaments dwelling at the barbed end of F-actin is available.
Here we investigate the possibility of bundle contraction in the absence of any sarcomeric organization, including motors dwelling at the filament barbed ends. After introducing a general bundle model in Sec. II, we show in Sec. III that underlying symmetries between contraction and extension [ Fig. 1(b)] imply that disordered bundle contraction requires non-identical motors and a non-linear elastic behavior of the filaments. Intuitively, non-identical motors induce mechanical frustration in a disordered actomyosin bundle, thus generating both contractile and extensile stresses. The filament non-linear behavior then allows the former to deform the bundle while resisting the latter, yielding overall contraction. Finally, in Sec. IV we discuss the limitations of our model and propose a minimal model for the contraction of a non-sarcomeric bundle.

II. BUNDLE MODEL
In order to make general statements about a disordered bundle of potentially complex internal geometry, we develop a detailed description of its mechanics without resorting to the simplifying mean-field approximations widely used in previous studies [19][20][21][22][23]. The main assumptions of our model are as follows. First, we assume that the velocity of a motor only depends on the force applied to it by the F-actin to which it is bound. Second, the average deformation of a thermally fluctuating section of filament in the bundle only depends on the force applied longitudinally at its ends. Third, we consider a stabilized bundle where F-actin polymerization and depolymerization do not occur. Finally, our model does not include the attachment-detachment dynamics of motors and cross-linkers to F-actin. The relevance of these choices is discussed in Sec. IV.
We describe a bundle of arbitrary geometry by subdividing it into interacting "units" of three types characterized by their lengths, velocities and tensions (Sec. II A). Introducing general notations to describe the physical connections between different units, we express the relationships between their forces and velocities as a function of their spatial arrangement (Sec. II B). The central physics of bundle mechanics are then described by introducing the filament force-extension relationships and the motor force-velocity relationships (Sec. II C). Conservation equations are then used to derive a compact description of the bundle amenable to further discussion (Sec. II D).

A. Linkers, filaments and junctions
Here we consider a single bundle constituted of F-actin, thick filaments of myosin molecular motors, and optionally passive actin cross-linkers. F-actin are assumed to be aligned in the z direction, and for simplicity we refer to the direction of positive (negative) z as the "right" ("left") in the following. To describe the bundle, we decompose it so as to distinguish three types of "units" (Fig. 2): • linker units, representing a whole myosin thick filament or an passive actin cross-linker (passive cross-linkers are equivalent to immobile motors). The total number of linker units in the system is denoted as n .
• junction units, representing the point of contact between a myosin thick filament or passive actin cross-linker on the one hand, and an F-actin on the other. The total number of junction units in the system is denoted as n .
• filament units, representing a portion of an Factin comprised between two junction units, or between one junction unit and an F-actin free end. We do not consider freely floating F-actin (filament units with two free ends). The total number of filament units in the system is denoted as n.
We label the filament units by i = 1, . . . , n, and denote by f i the tension of filament unit i (f i > 0 for a filament unit under extension). We allow filament units to bend away from the z-axis while maintaining their overall z orientation, and thus introduce filament unit i's contour length and end-to-end length as two independent variables L i and i , respectively. Finally, we introduce the velocity v r i of the rightmost monomer of filament unit i, and the velocity v l i of its leftmost monomer. These two velocities are Eulerian velocities; in that respect, one might think of junction 1 of Fig. 2 as a bridge and of the actin as the water that flows under it. Then the velocities v r 1 and v l 2 associated with filament units 1 and 2 are the velocities of the water just upstream and downstream (Color online) Schematic of a model bundle comprising three F-actin, one myosin thick filament bound in three different sites (linker 1) and one bound only once (linker 2). Colors indicate the division into linker units (red ), junction units (blue) and filament units (green). Our model does not impose any restrictions on the number of junctions associated with a linker, whether they involve one or more F-actin, or the polarity of the F-actin. from the bridge, as opposed to being the velocities of a specific fluid element. The notation is summarized in Fig. 3

(a).
We label junction units by i = 1, . . . , n , and linker units by i = 1, . . . , n . We introduce the velocity v i of linker unit i , i.e., the velocity of the bridge itself in our previous analogy. We choose to work in the reference frame where the center-of-mass of all linker units is motionless, which reads Note that the velocities v r i and v l i are absolute velocities defined in this reference frame, as opposed to being velocities relative to the neighboring linker. We delay labeling the other forces and velocities involved in junction and linker units until the introduction of convenient notations in the next section.

B. Bundle geometry and topology
To describe the physical connections between junction and filament units, we define the n × n matrices ρ and λ by For instance, the bundle represented in Fig (3) The usefulness of these matrices is illustrated by introducing the notation f = (f 1 , . . . , f n ) for the vector of all filament tensions, as well as similar notations L, , v r , v l , v . The matrix product ρf is a vector of length n whose i th component (ρf ) i is the tension of the filament unit that is the right-hand neighbor of junction unit i . This neighbor thus exerts a force (ρf ) i on junction unit i , while its left-hand neighbor exerts −(λf ) i . The sum of these two forces is equal and opposite to the force applied to junction unit i by its linker unit. This last force thus reads Further use of ρ and λ indicates that the velocities of the rightmost and leftmost actin monomers involved in junction unit i are (ρv l ) i and (λv r ) i , respectively The reasonings used here can be generalized in the following way: if x i is a quantity associated with filament i, then (ρx) i is associated to the right-hand neighbor of junction unit i . For instance, in the two previous paragraphs we considered x i = f i and x i = v l i or v r i , respectively. Just as this statement relates junction units to quantities associated with the neighboring filament units, we can conversely relate filament units to quantities associated with their junction unit neighbors as follows. If x i is associated with junction unit i , then (ρ T x ) i is associated with the junction unit that is the left neighbor of filament unit i if it has one, or is equal to zero if it does not (the superscript T denotes the matrix transpose). A similar statement holds for λ T . Combining this with the fact that each junction unit has exactly one right-hand and one left-hand neighbor, we find that (ρρ T x ) i = x i always. More generally, where 1l denotes the identity matrix. We further note that filament units may have one or two neighbors, which implies This means that matrix ρ T ρ is a projector onto the subspace of filament units that have a left-hand neighbor. A similar statement holds for λ T λ. This discussion implies that n < n.
To describe the polarity of the filament units, we introduce the diagonal n × n matrix Π, where Π ii = 1 if the pointed end of filament unit i points to the right, and Π ii = −1 if it points to the left. As an example, Fig. 2 has A similar diagonal polarity matrix Π is associated with junction units, and since neighboring filaments and junction units have the same polarity we have Turning to the linker units, we define the n ×n matrix γ by For instance, Fig. 2 has Therefore, the velocity of the linker unit connected to Each junction unit is connected to a linker unit, but one linker unit can be connected to several junction units, which implies n n . To obtain the forces associated with the linker units, we reason that if i experiences a force [(λ−ρ)f ] i from its linker unit, then the i exerts an equal and opposite force Fig. 3(c)]. Force balance imposes that the sum of the forces applied on any linker unit vanishes, and thus or, in vector notation: Although this condition comprises n scalar equations, those equations are not all independent. This can be seen by considering the mechanical subsystem formed by all junction and filament units (but not including the linker units). As inertia and friction against the background fluid are negligible, the sum of all forces applied to this system by junction units must vanish: where the first equality follows from the fact that γ has exactly one element equal to 1 per column and zeros everywhere else. As Eq. (13) is always trivially true, Eq. (12) expresses only n −1 linearly independent scalar conditions.

C. Filament elasticity and motor operation
Having defined notations for all lengths, forces and velocities in our system (Fig. 3), as well as enforced velocity continuity and force balance conditions, we turn to characterizing the more substantial physics of the junction and filament units.
Filament units are sections of semiflexible polymers shorter than or with a length comparable to their persistence length. We thus assume that the force required to hold a filament unit of given contour length L i in mechanical equilibrium is uniquely determined by specifying its end-to-end length i , which defines the force-extension relationship F : In the case of a thermally fluctuating polymer, i denotes the end-to-end length averaged over thermal fluctuations. Thermal bending of the filament unit moreover implies that i is smaller than L i even when f i = 0. Since filament units are shorter than the filament persistence length, we expect deformations of this kind ranging from zero to ≈ 20%. In vector notation, we write Motor operation at junction i is described by a functional relationship between the local velocity of the linker relative to the F-actin and the force applied to i : where the functionṼ i a priori depends on the polarity Π i i of the junction unit. This dependence can be explicitly determined by noting that the force-velocity relationship must not depend on our arbitrary choice of the direction of positive z. Reversing this choice is equivalent to reversing the sign of all velocities, forces and polarities. The only way for Eq. (16) to be invariant under this transformation is to write where function V i is the force-velocity relationship of motor i , an a priori nonlinear function independent of Π .
In vector notation,

D. Dynamical equations for bundle deformation
To provide a kinematic description of bundle contraction and extension, we write the conservation of F-actin contour length. The rate of change of the contour length of a filament unit is directly related to the velocity at which the neighboring junction units slide relative to F-actin. If i has a neighboring junction unit on its right-hand side, the velocity of the linker unit there is (λ T γ T v ) i , and its sliding velocity relative to the actin is (λ T γ T v − v r ) i . We rewrite this sliding velocity as (λ T γ T v − λ T λv r ) i , which is equal to the previous expression if i has a right-hand neighbor, and to zero if it does not. Using a similar reasoning for the left-hand side, we obtain The first (second) term on the right-hand side of this equation accounts for actin-linker sliding on the righthand (left-hand) side of the filaments units, and vanishes for the filaments units that do not have a right-hand (lefthand) neighbor junction unit. We describe the evolution of the end-to-end length of a filament unit in different ways depending on whether it has two junction unit neighbors or has a free end. In the former case, where the two terms in the right-hand side are the absolute velocities of the right and left neighbors of i, respectively. To describe a filament unit with one free end, we first note that it has vanishing tension. If the contour length of the filament unit is known, its end-to-end length is given by its force-extension relation Eq. (14) with f i = 0. Differentiating with respect to time, we obtain Defining the diagonal matrices where δ ij is the Kronecker delta, we rewrite Eq. (21) as for a filament unit with a free end [24]. Making use of Eq. (6) and its analog for λ T λ, we write an equation that describes filament units whether they have one or two neighboring junction units: We finally use Eqs. (4), (15) and (18) to eliminate f , ρv l and λv r in Eqs. (19) and (24). This yields where the nonlinear vector function v s ( , L) is the solution of the linear (in v ) system of equations formed by Eq. (1) and the following vector equation, obtained by combining Eqs. (12) and (15), then differentiating with respect to time and inserting Eqs. (25) into the result: Equation (26), just like Eq. (12), has only n − 1 independent scalar equations, and supplementing it with Eq. (1) thus results in a complete set of equations for v s . Finally, supplementing Eqs. (1), (25) and (26) with an initial condition [ (t = 0), L(t = 0)] completely specifies the dynamics of a bundle of arbitrary geometry.

III. SITUATIONS WITHOUT TELESCOPIC DEFORMATION
We now ask under what conditions contraction occurs. We are interested in bundles much longer than the size of any single one of their constituents (F-actin or motor). Significant contraction of such a bundle requires that it contracts throughout its length, as opposed to, e.g., contracting at its extremities while the bulk of the bundle retains a constant length. We thus focus on "telescopic deformation", whereby the end-to-end velocity of a bundle contracting under vanishing external load is proportional to its length. This constitutes the standard behavior of contractile actomyosin structures in vivo [1,[25][26][27] and in vitro [13]. Telescopic deformation is characteristic of systems formed of a serial arrangement of many independently deforming elements, often referred to as "contractile units" [13,26,27].
Here we demonstrate two requirements for telescopic deformation. In Sec. III A we show that it cannot arise if the motors all have identical force-velocity relationships. Sec. III B then tackles situations where motors with different force-velocity relationships are present. In that case, we prove that bundles lacking polarity organization (e.g., sarcomeres) comprised of linearly elastic (i.e., rigid) filaments do not undergo telescopic deformation. Such rigid filaments units represent situations where the filament persistence length is very large (e.g., bundles of microtubules and kinesin oligomers), or when the filament units themselves are very short (e.g., strongly crosslinked bundles, where the spacing between two junctions is small).
To completely determine the bundle dynamics, we need to specify an initial condition [ (t = 0), L(t = 0)]. We thus choose an arbitrary vector L 0 of length n and impose L(t = 0) = L 0 . To avoid confusion between the effects of motor-generated stresses, which are relevant for contractility, and those of bundle prestress, which are not, we consider bundles that are initially stress-free: This imposes (t = 0) = 0 , where 0 is the solution of the equation F ( 0 , L 0 ) = 0, and thus completely specifies the bundle initial condition.

A. Motors with identical force-velocity relationships
In a bundle where the motors have identical forcevelocity relationships, the junction units have identical spontaneous sliding velocities v * in the absence of applied force. Defining v * as the length n vector with all its components equal to v * , we can thus write We now demonstrate that under this assumption, all linker units are immobile throughout the dynamics and all right-pointing filaments undergo a uniform translation with constant velocity v * , while left-pointing filaments translate with −v * . We define the set of functions [ * (t), L * (t)] by and by defining * (t) as the solution of F [ * (t), L * (t)] = 0. This set manifestly satisfies the initial condition [ 0 , L 0 ] chosen above. Inserting [ * (t), L * (t)] into Eqs. (25), we further verify that these functions satisfies the equations of motion, implying that they describe the dynamics of the bundle. Using Eqs. (18) and (26), we find that v = 0 and ρv l = λv r = Π v * for all times, thus confirming that linker units are immobile and that the velocities associated with right-and left-pointing actin are v * and −v * , respectively.
In this regime, the maximum relative speed between any two actin filament units is 2|v * | irrespective of bundle geometry. Defining the contraction velocity of the bundle as the difference between the velocity of its leftmost and rightmost filaments, this implies that the contraction velocity cannot exceed the constant 2|v * |. It is thus impossible for the contraction velocity to scale linearly with bundle length, and telescopic contractility does not occur. Instead, filaments are segregated according to polarity, as observed experimentally in Ref. [28]. Note that this result does not require the bundle to be disordered.

B. Disordered bundles with linearly elastic filaments
In a bundle with arbitrary force-velocity relationships, the reasoning of the previous section does not apply, and some amount of contraction or extension is generally present. We thus ask whether bundles contract on average, and find that they do only if filament polarities are organized across the bundles, or if the filaments display nonlinear elastic behavior.
To prove this statement, we first give a mathematical description of bundles devoid of both polarity organization and filament nonlinear elastic behavior (Sec. III B 1), and then use Eqs. (25) and (26) to show that such bundles do not contract (Sec. III B 2).

Mathematical formulation
Consider an arbitrary bundle, which we denote by B. Bundle B is fully characterized by specifying matrices λ, ρ, Π, γ, the force-velocity functions V and the initial condition L 0 . Therefore, a population of bundles is fully characterized by specifying the distribution P(B) = P(λ, ρ, Π, γ, V , L 0 ) of the frequencies at which any possible bundle B arises in the population. For a population without any polarization organization, P(B) must be independent of Π, implying in particular that it is invariant under polarity reversal: This relationship clearly does not apply to a population of sarcomeres. Indeed, in a sarcomere static cross-linkers are restricted to the barbed ends of F-actin while active, mobile motors are found at the pointed ends. Thus the polarity-reversed image of a sarcomere is not a sarcomere; actually, inverting the polarities of filaments in Fig. 1(a) results in an extensile, not contractile, structure. Assuming in the following that Eq. (30) holds thus excludes sarcomeric contractility from our discussion. We further assume that filament units exhibit linear elastic behavior, which reads where ∂F /∂ and ∂F /∂L are constants, i.e. d dt This assumption is a good description of very stiff filaments, where i = L i for any filament unit i. Indeed, such filaments can be described by choosing the forceextension relationship which satisfies Eqs. (31) provided that K is a constant, and enforcing the limit K → +∞.

Proof of the property
To determine whether a bundle B = {λ, ρ, Π, γ, V , L 0 } is contractile or extensile, we imagine labeling its leftmost and rightmost points, and ask whether the distance L B between these two labels tends to increase or decrease with time. To calculate L B , we choose a path along the bundle's linker and filament units going from the left label to the right label as pictured in Fig. 4. We define B i as equal to 1 if the path considered crosses filament unit i from left to right, to −1 if it crosses it from right to left, and to 0 otherwise (see the caption of Fig. 4). The contraction velocity of the bundle can then be defined as: We use the following notation to refer to the solution of the equations of motion for bundle B: We now introduce bundleB as the polarity-reversed image of B, i.e.,B = {λ, ρ, −Π, γ, V , L 0 }. Substituting Eqs. (31) into Eqs. (25) and (26), we find that the dynamics ofB satisfies Combining this with Eqs. (33) and (34) while using the path B i = B i to assess the contraction ofB, we find We finally calculate the average contraction velocity over a population of bundles as where the sum runs over all possible bundles. Reorganizing this sum, we find where Eqs. (30) and (36) are used to derive the second and third equalities, respectively. Eq. (38) demonstrates that bundles without polarity organization or nonlinear elastic behavior do not contract or extend on average. This result does not depend on bundle structure or the form of the motor force-velocity relationships, and can easily be generalized to bundles pinned to a rigid substrate, or to include friction of the linker or filament units with the solvent. Mean-field modeling of dilute actomyosin gels with rigid filaments suggest that this symmetry-based reasoning could have a threedimensional counterpart [23]. However, geometrical nonlinearities in two or more dimensions can take on the role played by elastic nonlinearities in one-dimensional bundles, thus enabling contraction in disordered networks of rigid filaments [29].

IV. DISCUSSION
In this paper we consider bundles of sliding motors and filaments with arbitrary geometries and motor force-velocity relationships and show that their contractility requires 1. non-identical motors 2. and (a) either polarity organization (b) or nonlinear elastic response of the filaments.
While our model is framed in term of F-actin and myosin for clarity, our results are much more general and could equally apply to bundles comprised of other types of motors and filaments (e.g., kinesins and microtubules). Our description includes as a special case the well-understood contractility of striated muscle sarcomeres [ Fig. 1(a)]. Their architecture includes both identical myosin thick filaments and passive cross-linkers (which are mathematically equivalent to motors with velocity zero), thus satisfying condition (1). They moreover have a distinctive polarity organization, and thereby fulfill condition (2a). Similarly, bundles with motors whose velocities depend on their position relative to the filaments [19][20][21][22][23] generically break polarity-reversal symmetry, which results in polarity organization. Besides establishing the requirements for contractility, the formalism presented here can describe the dynamics of a wide range of contractile bundles. For instance, straightforward numerical simulations of Eqs. (25) and (26) could be used to describe bundle deformation as a function of the initial arrangement of the filament and motors. Comparing these predictions to experimental observations while varying these initial parameters could yield insight into the architecture of the bundles, which is currently lacking. Although such a study is beyond the scope of our work, a simplified version of our formalism still successfully predicts the onset of their contraction [30].
While the model used here is designed to describe a large class of contractile bundles, some of our assumptions are especially appropriate to describe the reconstituted bundles of Ref. [13]. F-actin is phalloidin-stabilized in this system, implying that no actin polymerizationdepolymerization takes place; myosin thick filaments are much shorter ( 300 nm) and thicker ( 50 nm) than F-actin ( 5 µm and 5 nm, respectively), justifying our assumption that they behave as rigid objects; and myosin thick filaments do not detach from the bundle on the time scales relevant for contraction, suggesting that the filament-motor attachment-detachment dynamics are inessential to contractility. Indeed, we show elsewhere that attachment-detachment can limit contractility under low myosin conditions [30].
Going beyond the specifics of the system studied in Ref. [13], it is interesting to discuss bundles where this attachment-detachment dynamics is not negligible. In the simplest such situation, motors undergo attachment and detachment at a constant rate, which tends to randomize their distribution in a filament polarityindependent manner. This fails to break the polarity- The Fenn effect impedes bundle contractility. (a) Simple contractile bundle comprised of fast (light) and slow (dark ) motors. Double-headed arrows indicate local contraction or expansion due to the movement of the motors relative to the filaments [ Fig. 1(b)]. (b) As tension builds following contraction, the motors detach and reorganize and contractility diminishes.
reversal symmetry discussed in this paper, and the resulting requirements for contractility are unchanged.
A more interesting question is to ask whether contractility could arise from the load-dependent detachment of myosin motors. Specifically, the detachment rate of myosin motors decreases under increasing load, a tendency known as the "Fenn effect" [31]. To assess its influence on contractility, we consider the simple bundle of rigid filaments illustrated in Fig. 5(a). This bundle comprises a contracting and an expanding region similar to those of Fig. 1(b). When faster motors happen to be concentrated in the contracting region, this bundle is contractile. Assuming that its ends are fixed, the bundle comes under tensile force during contraction. As a consequence, its expanding region comes under negative load and the motors there tend to detach. Conversely, the motors in the contracting region experience a positive load and tend to hold on to the filaments. On average, detached motors thus tend to diffuse from the expanding to the contracting region [ Fig. 5(b)]. The reduced number of motors in the expanding region implies that its resistance to the applied tension decreases, thus increasing its expansion rate. As the contracting region tends to recruit the slower motors from the expanding region, its contraction rate decreases. Overall, the Fenn effect thus tends to suppress contractile configurations rather than amplify them. This suggests that the Fenn effect cannot generate contractility in the absence of polarity organization or filament nonlinear elastic response.
The requirements derived in this paper offer insight into the contractility of actomyosin bundles devoid of sarcomere-like organization, for which no established contraction mechanism exists. As they do not satisfy condition (2a), we propose that they contract by fulfilling conditions (1) and (2b) [30]. Consider two antiparallel filaments interacting through several different motors with distinct speeds [ Fig. 6(a)]. As motors start to move relative to the filaments, stresses build in sections of the filament flanked by motors with different speeds. When the flanking motor proximal to the barbed end is faster than that proximal to the pointed end, compression arises. When it is slower, tension arises. F-actin responds nonlinearly to these stresses by buckling under Figure 6. Filament buckling as a mechanism for bundle contractility. (a) The presence of fast (white) and slow (gray) motors generically induce compressive (red ) and extensile (blue) stresses in filaments. (b) Buckling of the compressed filaments leads to an overall shortening of the bundle. Buckling of filament units of length Li = 500 nm occurs for compressive forces of order kBT p/L 2 i 0.16 pN, where p 10 µm is the F-actin persistence length. Myosin thick filaments typically exert forces of several piconewtons [34], which is sufficient to induce buckling, but remains smaller than the tension required to break F-actin ( 100 pN [35]).
compression while resisting extension, which has previ-ously been proposed to play a role in actomyosin contraction [32]. Following buckling of the compressed filament sections, fast motors are free to move quickly while the others move slowly. This results in the growth of the compressed sections and shrinkage of the extended ones, and thus in overall bundle contraction [ Fig. 6(b)]. Experimental observations suggest that this mechanism could be at the origin of contraction in reconstituted actomyosin bundles [30], and that F-actin buckling can occur in cells [33]. These results offer an interesting new perspective on mechanisms underlying actomyosin bundle contraction in vivo.