Monte Carlo study of multiply crosslinked semiflexible polymer networks

Strain-Controlled Critical Slowing Down in the Rheology of Disordered Networks

Networks and dense suspensions frequently reside near a boundary between soft (or fluidlike) and rigid (or solidlike) regimes. Transitions between these regimes can be driven by changes in structure, density, or applied stress or strain. In general, near the onset or loss of rigidity in these systems, dissipation-limiting heterogeneous nonaffine rearrangements dominate the macroscopic viscoelastic response, giving rise to diverging relaxation times and power-law rheology. Here, we describe a simple quantitative relationship between nonaffinity and the excess viscosity. We test this nonaffinity-viscosity relationship computationally and demonstrate its rheological consequences in simulations of strained filament networks and dense suspensions. We also predict critical signatures in the rheology of semiflexible and stiff biopolymer networks near the strain stiffening transition.

Nonaffine Deformation of Semiflexible Polymer and Fiber Networks

Networks of semiflexible or stiff polymers such as most biopolymers are known to deform inhomogeneously when sheared. The effects of such nonaffine deformation have been shown to be much stronger than for flexible polymers. To date, our understanding of nonaffinity in such systems is limited to simulations or specific 2D models of athermal fibers. Here, we present an effective medium theory for nonaffine deformation of semiflexible polymer and fiber networks, which is general to both 2D and 3D and in both thermal and athermal limits. The predictions of this model are in good agreement with both prior computational and experimental results for linear elasticity. Moreover, the framework we introduce can be extended to address nonlinear elasticity and network dynamics.

Structural Features and Nonlinear Rheology of Self-Assembled Networks of Cross-Linked Semiflexible Polymers

Disordered networks of semiflexible filaments are common support structures in biology. Familiar examples include fibrous matrices in blood clots, bacterial biofilms, and essential components of cells and tissues of plants, animals, and fungi. Despite the ubiquity of these networks in biomaterials, we have only a limited understanding of the relationship between their structural features and their highly strain-sensitive mechanical properties. In this work, we perform simulations of three-dimensional networks produced by the irreversible formation of cross-links between linker-decorated semiflexible filaments. We characterize the structure of networks formed by a simple diffusion-dependent assembly process and measure their associated steady-state rheological features at finite temperature over a range of applied prestrains that encompass the strain-stiffening transition. We quantify the dependence of network connectivity on cross-linker availability and detail the associated connectivity dependence of both linear elasticity and nonlinear strain-stiffening behavior, drawing comparisons with prior experimental measurements of the cross-linker concentration-dependent elasticity of actin gels.

Biomechanics of cells and subcellular components: A comprehensive review of computational models and applications

Cells are a fundamental structural, functional and biological unit for all living organisms. Up till now, considerable efforts have been made to study the responses of single cells and subcellular components to an external load, and understand the biophysics underlying cell rheology, mechanotransduction and cell functions using experimental and in silico approaches. In the last decade, computational simulation has become increasingly attractive due to its critical role in interpreting experimental data, analysing complex cellular/subcellular structures, facilitating diagnostic designs and therapeutic techniques, and developing biomimetic materials. Despite the significant progress, developing comprehensive and accurate models of living cells remains a grand challenge in the 21st century. To understand current state of the art, this review summarises and classifies the vast array of computational biomechanical models for cells. The article covers the cellular components at multi-spatial levels, that is, protein polymers, subcellular components, whole cells and the systems with scale beyond a cell. In addition to the comprehensive review of the topic, this article also provides new insights into the future prospects of developing integrated, active and high-fidelity cell models that are multiscale, multi-physics and multi-disciplinary in nature. This review will be beneficial for the researchers in modelling the biomechanics of subcellular components, cells and multiple cell systems and understanding the cell functions and biological processes from the perspective of cell mechanics.

Non-affinity in multi-material mechanical metamaterials

Non-affine deformations enable mechanical metamaterials to achieve their unusual properties while imposing implications for their structural integrity. The presence of multiple phases with different mechanical properties results in additional non-affinity of the deformations, a phenomenon that has never been studied before in the area of extremal mechanical metamaterials. Here, we studied the degree of non-affinity, [Formula: see text], resulting from the random substitution of a fraction of the struts,[Formula: see text], that make up a lattice structure and are printed using a soft material (elastic modulus = [Formula: see text]) by those printed using a hard material ([Formula: see text]). Depending on the unit cell angle (i.e., [Formula: see text] = 60°, 90°, or 120°), the lattice structures exhibited negative, near-zero, or positive values of the Poisson's ratio, respectively. We found that the auxetic structures exhibit the highest levels of non-affinity, followed by the structures with positive and near-zero values of the Poisson's ratio. We also observed an increase in [Formula: see text] with [Formula: see text] and [Formula: see text] until [Formula: see text] =104 and [Formula: see text]= 75%-90% after which [Formula: see text] saturated. The dependency of [Formula: see text] upon [Formula: see text] was therefore found to be highly asymmetric. The positive and negative values of the Poisson's ratio were strongly correlated with [Formula: see text]. Interestingly, achieving extremely high or extremely low values of the Poisson's ratio required highly affine deformations. In conclusion, our results clearly show the importance of considering non-affinity when trying to achieve a specific set of mechanical properties and underscore the structural integrity implications in multi-material mechanical metamaterials.

Chain length-dependent luminescence in acceptor-doped conjugated polymers

Semiconducting polymers doped with a minority fraction of energy transfer acceptors feature a sensitive coupling between chain conformation and fluorescence emission, that can be harnessed for advanced solution-based molecular sensing and diagnostics. While it is known that chain length strongly affects chain conformation, and its response to external cues, the effects of chain length on the emission patterns in chromophore-doped conjugated polymers remains incompletely understood. In this paper, we explore chain-length dependent emission in two different acceptor-doped polyfluorenes. We show how the binomial distribution of acceptor incorporation, during the probabilistic polycondensation reaction, creates a strong chain-length dependency in the optical properties of this class of luminescent polymers. In addition, we also find that the intrachain exciton migration rate is chain-length dependent, giving rise to additional complexity. Both effects combined, make for the need to develop sensoric conjugated polymers of improved monodispersity and chemical homogeneity, to improve the accuracy of conjugated polymer based diagnostic approaches.

Stress-stabilized subisostatic fiber networks in a ropelike limit

The mechanics of disordered fibrous networks such as those that make up the extracellular matrix are strongly dependent on the local connectivity or coordination number. For biopolymer networks this coordination number is typically between 3 and 4. Such networks are sub-isostatic and linearly unstable to deformation with only central force interactions, but exhibit a mechanical phase transition between floppy and rigid states under strain. The introduction of weak bending interactions stabilizes these networks and suppresses the critical signatures of this transition. We show that applying external stress can also stabilize subisostatic networks with only tensile central force interactions, i.e., a ropelike potential. Moreover, we find that the linear shear modulus shows a power-law scaling with the external normal stress, with a non-mean-field exponent. For networks with finite bending rigidity, we find that the critical stain shifts to lower values under prestress.

Cohesive and adhesive properties of crosslinked semiflexible biopolymer networks

Biomolecular semiflexible polymer networks with persistence lengths well above those of single polymeric chains serve important structural and adhesive roles in biology, biomaterials, food science and many other fields. While relationships between the structure and viscoelasticity of semiflexible polymer networks have been previously investigated, it remains challenging to systematically relate fibril and network properties to cohesive and adhesive properties that govern the function of these materials. To address this issue, here we utilize coarse-grained molecular dynamics simulations to thoroughly elucidate how the work of adhesion of a semiflexible polymer network to a surface depends on crosslink density and fibril persistence length. Two emergent characteristics of the network are its elasticity and its interfacial energy with the surface. Stiff networks that are either highly crosslinked or have high persistence length fibrils tend to have lower interfacial energy, and consequently, lower work of adhesion. For lightly crosslinked networks with flexible fibrils, considerable strain energy must be stored within the adhesive during detachment, which creates an additional penalty to detachment. Increasing persistence length while keeping crosslink density constant leads to porous, low density networks, leading to an optimal fibril persistence length at which maximum work of adhesion per mass density is attained for a given crosslink density. For any given fibril persistence length, increasing crosslink density has a slightly negative effect on network mass density and interfacial energy. A critical crosslink density is found, below which the networks have no significant load-bearing capacity. Lightly crosslinked networks above this threshold absorb more strain energy during desorption and consequently possess greater work of adhesion. The conflict between mass density and stiffness results in a non-monotonic trend between the ratio of work of adhesion to interfacial energy and persistence length. These findings provide physical insight into the adhesive mechanisms of biomaterials based on crosslinked semiflexible polymer networks, and reveal important design guidelines for bio-adhesives.

Catch bonding in the forced dissociation of a polymer endpoint

Applying a force to certain supramolecular bonds may initially stabilize them, manifested by a lower dissociation rate. We show that this behavior, known as catch bonding and by now broadly reported in numerous biophysics bonds, is generically expected when either or both the trapping potential and the force applied to the bond possess some degree of nonlinearity. We enumerate possible scenarios and for each identify the possibility and, if applicable, the criterion for catch bonding to occur. The effect is robustly predicted by Kramers theory and Mean First Passage Time theory and confirmed in direct molecular dynamics simulation. Among the catch scenarios, one plays out essentially any time the force on the bond originates in a polymeric object, implying that some degree of catch bond behavior is to be expected in many settings relevant to polymer network mechanics or optical tweezer experiments.

Linear rheology of reversibly cross-linked biopolymer networks

We suggest a simple model for reversible cross-links, binding, and unbinding to/from a network of semiflexible polymers. The resulting frequency dependent response of the network to an applied shear is calculated via Brownian dynamics simulations. It is shown to be rather complex with the time scale of the linkers competing with the excitations of the network. If the lifetime of the linkers is the longest time scale, as is indeed the case in most biological networks, then a distinct low frequency peak of the loss modulus develops. The storage modulus shows a corresponding decay from its plateau value, which for irreversible cross-linkers extends all the way to the static limit. This additional relaxation mechanism can be controlled by the relative weight of reversible and irreversible linkers.

Bimodal probability density characterizes the elastic behavior of a semiflexible polymer in 2D under compression

We explore the elastic behavior of a wormlike chain under compression in terms of exact solutions for the associated probability densities. Strikingly, the probability density for the end-to-end distance projected along the applied force exhibits a bimodal shape in the vicinity of the critical Euler buckling force of an elastic rod, reminiscent of the smeared discontinuous phase transition of a finite system. These two modes reflect the almost stretched and the S-shaped configuration of a clamped polymer induced by the compression. Moreover, we find a bimodal shape of the probability density for the transverse fluctuations of the free end of a cantilevered polymer as fingerprint of its semiflexibility. In contrast to clamped polymers, free polymers display a circularly symmetric probability density and their distributions are identical for compression and stretching forces.

Hindrances to precise recovery of cellular forces in fibrous biopolymer networks

How cells move through the three-dimensional extracellular matrix (ECM) is of increasing interest in attempts to understand important biological processes such as cancer metastasis. Just as in motion on flat surfaces, it is expected that experimental measurements of cell-generated forces will provide valuable information for uncovering the mechanisms of cell migration. However, the recovery of forces in fibrous biopolymer networks may suffer from large errors. Here, within the framework of lattice-based models, we explore possible issues in force recovery by solving the inverse problem: how can one determine the forces cells exert to their surroundings from the deformation of the ECM? Our results indicate that irregular cell traction patterns, the uncertainty of local fiber stiffness, the non-affine nature of ECM deformations and inadequate knowledge of network topology will all prevent the precise force determination. At the end, we discuss possible ways of overcoming these difficulties.

Self-Consistent Field Lattice Model for Polymer Networks

A lattice model based on polymer self-consistent field theory is developed to predict the equilibrium statistics of arbitrary polymer networks. For a given network topology, our approach uses moment propagators on a lattice to self-consistently construct the ensemble of polymer conformations and cross-link spatial probability distributions. Remarkably, the calculation can be performed “in the dark”, without any prior knowledge on preferred chain conformations or cross-link positions. Numerical results from the model for a test network exhibit close agreement with molecular dynamics simulations, including when the network is strongly sheared. Our model captures nonaffine deformation, mean-field monomer interactions, cross-link fluctuations, and finite extensibility of chains, yielding predictions that differ markedly from classical rubber elasticity theory for polymer networks. By examining polymer networks with different degrees of interconnectivity, we gain insight into cross-link entropy, an important quantity in the macroscopic behavior of gels and self-healing materials as they are deformed.

J-shaped stress-strain diagram of collagen fibers: Frame tension of triangulated surfaces with fixed boundaries

We present Monte Carlo data of the stress-strain diagrams obtained using two different triangulated surface models. The first is the canonical surface model of Helfrich and Polyakov (HP), and the second is a Finsler geometry (FG) model. The shape of the experimentally observed stress-strain diagram is called J-shaped. Indeed, the diagram has a plateau for the small strain region and becomes linear in the relatively large strain region. Because of this highly nonlinear behavior, the J-shaped diagram is far beyond the scope of the ordinary theory of elasticity. Therefore, the mechanism behind the J-shaped diagram still remains to be clarified, although it is commonly believed that the collagen degrees of freedom play an essential role. We find that the FG modeling technique provides a coarse-grained picture for the interaction between the collagen and the bulk material. The role of the directional degrees of freedom of collagen molecules or fibers can be understood in the context of FG modeling. We also discuss the reason for why the J-shaped diagram cannot (can) be explained by the HP (FG) model.

Exact solution for the force-extension relation of a semiflexible polymer under compression

Exact solutions for the elastic and thermodynamic properties for the wormlike chain model are elaborated in terms of Mathieu functions. The smearing of the classical Euler buckling instability for clamped polymers is analyzed for the force-extension relation. Interestingly, at strong compression forces the thermal fluctuations lead to larger elongations than for the elastic rod. The susceptibility defined as the derivative of the force-extension relation displays a prominent maximum at a force that approaches the critical Euler buckling force as the persistence length is increased. We also evaluate the excess entropy and heat capacity induced by the compression and find that they vary nonmonotonically with the load. These findings are corroborated by pseudo-Brownian simulations.

Critical behaviour in the nonlinear elastic response of hydrogels

In this paper we study the elastic response of synthetic hydrogels to an applied shear stress. The hydrogels studied here have previously been shown to mimic the behaviour of biopolymer networks when they are sufficiently far above the gel point. We show that near the gel point they exhibit an elastic response that is consistent with the predicted critical behaviour of networks near or below the isostatic point of marginal stability. This point separates rigid and floppy states, distinguished by the presence or absence of finite linear elastic moduli. Recent theoretical work has also focused on the response of such networks to finite or large deformations, both near and below the isostatic point. Despite this interest, experimental evidence for the existence of criticality in such networks has been lacking. Using computer simulations, we identify critical signatures in the mechanical response of sub-isostatic networks as a function of applied shear stress. We also present experimental evidence consistent with these predictions. Furthermore, our results show the existence of two distinct critical regimes, one of which arises from the nonlinear stretch response of semi-flexible polymers.

Two fundamental mechanisms govern the stiffening of cross-linked networks

Biopolymer networks, such as those constituting the cytoskeleton of a cell or biological tissue, exhibit a nonlinear strain-stiffening behavior when subjected to large deformations. Interestingly, rheological experiments on various in vitro biopolymer networks have shown similar strain-stiffening trends regardless of the differences in their microstructure or constituents, suggesting a universal stiffening mechanism. In this article, we use computer simulations of a random network comprised of cross-linked biopolymer-like fibers to substantiate the notion that this universality lies in the existence of two fundamental stiffening mechanisms. After showing that the large strain response is accompanied by the development of a stress path, i.e., a percolating path of axially stressed fibers and cross-links, we demonstrate that the strain stiffening can be caused by two distinctly different mechanisms: 1) the pulling out of stress-path undulations; and 2) reorientation of the stress path. The former mechanism is bending-dominated and can be recognized by a power-law dependence with exponent 3/2 of the shear modulus on stress, whereas the latter mechanism is stretching-dominated and characterized by a power-law exponent 1/2. We demonstrate how material properties of the constituents, as well as the network microstructure, can affect the transition between the two stiffening mechanisms and, as such, control the dominant power-law scaling behavior.

Elastic regimes of subisostatic athermal fiber networks

Athermal models of disordered fibrous networks are highly useful for studying the mechanics of elastic networks composed of stiff biopolymers. The underlying network architecture is a key aspect that can affect the elastic properties of these systems, which include rich linear and nonlinear elasticity. Existing computational approaches have focused on both lattice-based and off-lattice networks obtained from the random placement of rods. It is not obvious, a priori, whether the two architectures have fundamentally similar or different mechanics. If they are different, it is not clear which of these represents a better model for biological networks. Here, we show that both approaches are essentially equivalent for the same network connectivity, provided the networks are subisostatic with respect to central force interactions. Moreover, for a given subisostatic connectivity, we even find that lattice-based networks in both two and three dimensions exhibit nearly identical nonlinear elastic response. We provide a description of the linear mechanics for both architectures in terms of a scaling function. We also show that the nonlinear regime is dominated by fiber bending and that stiffening originates from the stabilization of subisostatic networks by stress. We propose a generalized relation for this regime in terms of the self-generated normal stresses that develop under deformation. Different network architectures have different susceptibilities to the normal stress but essentially exhibit the same nonlinear mechanics. Such a stiffening mechanism has been shown to successfully capture the nonlinear mechanics of collagen networks.

Stability and anomalous entropic elasticity of subisostatic random-bond networks

We study the elasticity of thermalized spring networks under an applied bulk strain. The networks considered are subisostatic random-bond networks that, in the athermal limit, are known to have vanishing bulk and linear shear moduli at zero bulk strain. Above a bulk strain threshold, however, these networks become rigid, although surprisingly the shear modulus remains zero until a second, higher, strain threshold. We find that thermal fluctuations stabilize all networks below the rigidity transition, resulting in systems with both finite bulk and shear moduli. Our results show a T(0.66) temperature dependence of the moduli in the region below the bulk strain threshold, resulting in networks with anomalously high rigidity as compared to ordinary entropic elasticity. Furthermore, we find a second regime of anomalous temperature scaling for the shear modulus at its zero-temperature rigidity point, where it scales as T(0.5), behavior that is absent for the bulk modulus since its athermal rigidity transition is discontinuous.

Nonlinearities of biopolymer gels increase the range of force transmission

We present a model of biopolymer gels that includes two types of elastic nonlinearities, stiffening under extension and softening (due to buckling) under compression, to predict the elastic anisotropy induced by both external as well as internal (e.g., due to cell contractility) stresses in biopolymer gels. We show how the stretch-induced anisotropy and the strain-stiffening nonlinearity increase both the amplitude and power-law range of transmission of internal, contractile, cellular forces, and relate this to recent experiments.

Advances in the mechanical modeling of filamentous actin and its cross-linked networks on multiple scales

The protein actin is a part of the cytoskeleton and, therefore, responsible for the mechanical properties of the cells. Starting with the single molecule up to the final structure, actin creates a hierarchical structure of several levels exhibiting a remarkable behavior. The hierarchy spans several length scales and limitations in computational power; therefore, there is a call for different mechanical modeling approaches for the different scales. On the molecular level, we may consider each atom in molecular dynamics simulations. Actin forms filaments by combining the molecules into a double helix. In a model, we replace molecular subdomains using coarse-graining methods, allowing the investigation of larger systems of several atoms. These models on the nanoscale inform continuum mechanical models of large filaments, which are based on worm-like chain models for polymers. Assemblies of actin filaments are connected with cross-linker proteins. Models with discrete filaments, so-called Mikado models, allow us to investigate the dependence of the properties of networks on the parameters of the constituents. Microstructurally motivated continuum models of the networks provide insights into larger systems containing cross-linked actin networks. Modeling of such systems helps to gain insight into the processes on such small scales. On the other hand, they call for verification and hence trigger the improvement of established experiments and the development of new methods.

Normal stresses in elastic networks

When loaded in simple shear deformation, polymeric materials may develop so-called normal stresses: stresses perpendicular to the direction of the applied shear. These normal stresses are intrinsically nonlinear: basic symmetry considerations dictate they may only enter at O(γ(2)), with γ the dimensionless shear strain. There is no fundamental restriction on their sign, and normal stresses may be positive (pushing outward) or negative (pulling inward). Most materials tend to dilate in the normal direction, but a wide variety of biopolymer networks including fibrin and actin gels have been reported to present anomalously large, negative normal stresses-a feature which has been ascribed to the intrinsic elastic nonlinearity of semiflexible fibers. In this work, we present analytical results on a model nonlinear network, which we expand to the required nonlinear order to show that due to geometric, rather than elastic, nonlinearities (negative) normal stresses generically arise in filamentous networks-even in networks composed of linear, Hookean springs. We investigate analytically and numerically how the subsequent addition of elastic nonlinearities, nonaffine deformations, and filament persistence through cross-linkers augment this basic behavior.

Fluctuation-stabilized marginal networks and anomalous entropic elasticity

We study the elastic properties of thermal networks of Hookean springs. In the purely mechanical limit, such systems are known to have a vanishing rigidity when their connectivity falls below a critical, isostatic value. In this work, we show that thermal networks exhibit a nonzero shear modulus G well below the isostatic point and that this modulus exhibits an anomalous, sublinear dependence on temperature T. At the isostatic point, G increases as the square root of T, while we find G∝Tα below the isostatic point, where α≃0.8. We show that this anomalous T dependence is entropic in origin.

Elasticity of a filamentous kagome lattice

The diluted kagome lattice, in which bonds are randomly removed with probability 1-p, consists of straight lines that intersect at points with a maximum coordination number of 4. If lines are treated as semiflexible polymers and crossing points are treated as cross-links, this lattice provides a simple model for two-dimensional filamentous networks. Lattice-based effective-medium theories and numerical simulations for filaments modeled as elastic rods, with stretching modulus μ and bending modulus κ, are used to study the elasticity of this lattice as functions of p and κ. At p=1, elastic response is purely affine, and the macroscopic elastic modulus G is independent of κ. When κ=0, the lattice undergoes a first-order rigidity-percolation transition at p=1. When κ>0, G decreases continuously as p decreases below one, reaching zero at a continuous rigidity-percolation transition at p=p(b)≈0.605 that is the same for all nonzero values of κ. The effective-medium theories predict scaling forms for G, which exhibit crossover from bending-dominated response at small κ/μ to stretching-dominated response at large κ/μ near both p=1 and p(b), that match simulations with no adjustable parameters near p=1. The affine response as p→1 is identified with the approach to a state with sample-crossing straight filaments treated as elastic rods.

Rheology of heterotypic collagen networks

Collagen fibrils are the main structural element of connective tissues. In many tissues, these fibrils contain two fibrillar collagens (types I and V) in a ratio that changes during tissue development, regeneration, and various diseases. Here we investigate the influence of collagen composition on the structure and rheology of networks of purified collagen I and V, combining fluorescence and atomic force microscopy, turbidimetry, and rheometry. We demonstrate that the network stiffness strongly decreases with increasing collagen V content, even though the network structure does not substantially change. We compare the rheological data with theoretical models for rigid polymers and find that the elasticity is dominated by nonaffine deformations. There is no analytical theory describing this regime, hampering a quantitative interpretation of the influence of collagen V. Our findings are relevant for understanding molecular origins of tissue biomechanics and for guiding rational design of collagenous biomaterials for biomedical applications.

The influence of disorder on deformations in semiflexible networks

We present a model that assesses the different elastic responses of a semiflexible network, which either (i) is constrained to deform in an affine way or (ii) is permitted to thermally fluctuate and deviate from affine response. The thermal, non-affine response of the network is achieved using a Metropolis Monte Carlo algorithm with dynamic step size. We find that non-affine deformations soften the network dramatically at low strains and make the eventual nonlinear strain stiffening far more pronounced. We show that the effect of these non-affine deformations are very sensitive to the degree variation in the lengths of filaments connecting cross-links. Where there is high variation, non-affine deformations allow internal stresses to relax, giving rise to a smaller range of tensile forces in filaments and a dramatic reduction of network stiffness. This highlights that non-affine deformations are crucial in small strain response of stiff polymer networks.

Internal stresses, normal modes, and nonaffinity in three-dimensional biopolymer networks

We numerically investigate deformations and modes of networks of semiflexible biopolymers as a function of crosslink coordination number z and strength of bending and stretching energies. In equilibrium filaments are under internal stress, and the networks exhibit shear rigidity below the Maxwell isostatic point. In contrast to two-dimensional networks, ours exhibit nonaffine bending-dominated response in all rigid states, including those near the maximum of z=4 when bending energies are less than stretching ones.

Frequency-dependent stiffening of semiflexible networks: a dynamical nonaffine to affine transition

By combining the force-extension relation of single semiflexible polymers with a Langevin equation to capture the dissipative dynamics of chains moving through a viscous medium we study the dynamical response of cross-linked biopolymer materials. We find that at low frequencies the network deformations are highly nonaffine, and show a low plateau in the modulus. At higher frequencies, this nonaffinity decreases while the elastic modulus increases. With increasing frequency, more and more nonaffine network relaxation modes are suppressed, resulting in a stiffening. This effect is fundamentally different from the high-frequency stiffening due to the single-filament relaxation modes [F. Gittes and F. C. MacKintosh, Phys. Rev. E 58, R1241 (1998)], not only in terms of its mechanism but also in its resultant scaling: G'(ω) ∼ ω(α) with α > 3/4. This may determine nonlinear material properties at low, physiologically relevant frequencies.

Semiflexible filamentous composites

Inspired by the ubiquity of composite filamentous networks in nature, we investigate models of biopolymer networks that consist of interconnected floppy and stiff filaments. Numerical simulations carried out in three dimensions allow us to explore the microscopic partitioning of stresses and strains between the stiff and floppy fractions cs and cf and reveal a nontrivial relationship between the mechanical behavior and the relative fraction of stiff polymer: when there are few stiff polymers, nonpercolated stiff "inclusions" are protected from large deformations by an encompassing floppy matrix, while at higher fractions of stiff material the stiff network is independently percolated and dominates the mechanical response.

Actin filament length tunes elasticity of flexibly cross-linked actin networks

Networks of the cytoskeletal biopolymer actin cross-linked by the compliant protein filamin form soft gels that stiffen dramatically under shear stress. We demonstrate that the elasticity of these networks shows a strong dependence on the mean length of the actin polymers, unlike networks with small, rigid cross-links. This behavior is in agreement with a model of rigid filaments connected by multiple flexible linkers.

Structural hierarchy governs fibrin gel mechanics

Fibrin gels are responsible for the mechanical strength of blood clots, which are among the most resilient protein materials in nature. Here we investigate the physical origin of this mechanical behavior by performing rheology measurements on reconstituted fibrin gels. We find that increasing levels of shear strain induce a succession of distinct elastic responses that reflect stretching processes on different length scales. We present a theoretical model that explains these observations in terms of the unique hierarchical architecture of the fibers. The fibers are bundles of semiflexible protofibrils that are loosely connected by flexible linker chains. This architecture makes the fibers 100-fold more flexible to bending than anticipated based on their large diameter. Moreover, in contrast with other biopolymers, fibrin fibers intrinsically stiffen when stretched. The resulting hierarchy of elastic regimes explains the incredible resilience of fibrin clots against large deformations.