Elasticity of 3D networks with rigid filaments and compliant crosslinks

FilamentSensor 2.0: An open-source modular toolbox for 2D/3D cytoskeletal filament tracking

Cytoskeletal pattern formation and structural dynamics are key to a variety of biological functions and a detailed and quantitative analysis yields insight into finely tuned and well-balanced homeostasis and potential pathological alterations. High content life cell imaging of fluorescently labeled cytoskeletal elements under physiological conditions is nowadays state-of-the-art and can record time lapse data for detailed experimental studies. However, systematic quantification of structures and in particular the dynamics (i.e. frame-to-frame tracking) are essential. Here, an unbiased, quantitative, and robust analysis workflow that can be highly automatized is needed. For this purpose we upgraded and expanded our fiber detection algorithm FilamentSensor (FS) to the FilamentSensor 2.0 (FS2.0) toolbox, allowing for automatic detection and segmentation of fibrous structures and the extraction of relevant data (center of mass, length, width, orientation, curvature) in real-time as well as tracking of these objects over time and cell event monitoring.

A computational framework for biomaterials containing three-dimensional random fiber networks based on the affine kinematics

Understanding the structure-function relationship of biomaterials can provide insights into different diseases and advance numerous biomedical applications. This paper presents a finite element-based computational framework to model biomaterials containing a three-dimensional fiber network at the microscopic scale. The fiber network is synthetically generated by a random walk algorithm, which uses several random variables to control the fiber network topology such as fiber orientations and tortuosity. The geometric information of the generated fiber network is stored in an array-like data structure and incorporated into the nonlinear finite element formulation. The proposed computational framework adopts the affine fiber kinematics, based on which the fiber deformation can be expressed by the nodal displacement and the finite element interpolation functions using the isoparametric relationship. A variational approach is developed to linearize the total strain energy function and derive the nodal force residual and the stiffness matrix required by the finite element procedure. Four numerical examples are provided to demonstrate the capabilities of the proposed computational framework, including a numerical investigation about the relationship between the proposed method and a class of anisotropic material models, a set of synthetic examples to explore the influence of fiber locations on material local and global responses, a thorough mesh-sensitivity analysis about the impact of mesh size on various numerical results, and a detailed case study about the influence of material structures on the performance of eggshell-membrane-hydrogel composites. The proposed computational framework provides an efficient approach to investigate the structure-function relationship for biomaterials that follow the affine fiber kinematics.

The elastic properties and deformation mechanisms of actin filament networks crosslinked by filamins

As a substructure of cell cytoskeleton, the crosslinked actin filament networks (CAFNs) play a major role in different cell functions, however, the elastic properties and the deformation mechanisms of CAFNs still remain to be understood. In this paper, a novel three-dimensional (3D) finite element (FE) model has been developed to mimic the mechanical properties of actin filament (F-actin) networks crosslinked by filamin A (FLNA). The simulation results indicate that although the Young's modulus of CAFNs varies in different directions for each random model, the statistical mean value is in-plane isotropic. The crosslinking density and the actin filament volume fraction are found to strongly affect the in-plane shear modulus of CAFNs. The simulation results agree well with the relevant experimental results. In addition, an L-shaped cantilever beam model has been developed for dimensional analysis on the shear stiffness of CAFNs and for quantifying the deformation mechanisms. It has been demonstrated that the in-plane shear modulus of CAFNs is mainly dominated by FLNA (i.e., cross-linkers), and that the bending and torsion deformations of FLNA have almost the same contribution to the stiffness of CAFNs. It has also been found that the stiffness of CAFNs is almost insensitive to the variation of the Poisson's ratios of FLNA and actin filament in the range from 0.29 to 0.499.

Mechanisms of the Complex Thermo-Mechanical Behavior of Polymer Glass Across a Wide Range of Temperature Variations

This paper aims to explore the mechanisms of the complex thermo-mechanical behavior of polymer glass across a wide range of temperature variations. To this end, the free vibration frequency spectrum of simply supported poly(methyl methacrylate) (PMMA) beams was thoroughly investigated with the aid of the impulse excitation technique. It was found that the amplitude ratio of the multiple peaks in the frequency spectrum is a strongly dependent on temperature, and that the peaks correspond to the multiple vibrational modes of the molecular network of PMMA. At a low temperature, the vibration is dominated by the overall microstructure of PMMA. With increasing the temperature, however, the contribution of the sub-microstructures is retarded by β relaxation. Above 80 °C, the vibration is fully dominated by the microstructure after relaxation. The relaxation time at the transition temperature is of the same order of the vibration period, confirming the contribution of β relaxation. These findings provide a precise method for establishing reliable physical-based constitutive models of polymer glass.

Co-Entangled Actin-Microtubule Composites Exhibit Tunable Stiffness and Power-Law Stress Relaxation

We use optical tweezers microrheology and fluorescence microscopy to characterize the nonlinear mesoscale mechanics and mobility of in vitro co-entangled actin-microtubule composites. We create a suite of randomly oriented, well-mixed networks of actin and microtubules by co-polymerizing varying ratios of actin and tubulin in situ. To perturb each composite far from equilibrium, we use optical tweezers to displace an embedded microsphere a distance greater than the lengths of the filaments at a speed much faster than their intrinsic relaxation rates. We simultaneously measure the force the filaments exert on the bead and the subsequent force relaxation. We find that the presence of a large fraction of microtubules (>0.7) is needed to substantially increase the measured force, which is accompanied by large heterogeneities in force response. Actin minimizes these heterogeneities by reducing the mesh size of the composites and supporting microtubules against buckling. Composites also undergo a sharp transition from strain softening to stiffening when the fraction of microtubules (ϕT) exceeds 0.5, which we show arises from faster poroelastic relaxation and suppressed actin bending fluctuations. The force after bead displacement relaxes via power-law decay after an initial period of minimal relaxation. The short-time relaxation profiles (t < 0.06 s) arise from poroelastic and bending contributions, whereas the long-time power-law relaxation is indicative of filaments reptating out of deformed entanglement constraints. The scaling exponents for the long-time relaxation exhibit a nonmonotonic dependence on ϕT, reaching a maximum for equimolar composites (ϕT = 0.5), suggesting that reptation is fastest in ϕT = 0.5 composites. Corresponding mobility measurements of steady-state actin and microtubules show that both filaments are indeed the most mobile in ϕT = 0.5 composites. This nonmonotonic dependence of mobility on ϕT demonstrates the important interplay between mesh size and filament rigidity in polymer networks and highlights the surprising emergent properties that can arise in composites.

Effect of crosslink torsional stiffness on elastic behavior of semiflexible polymer networks

Networks of semiflexible filaments are building blocks of different biological and structural materials such as cytoskeleton and extracellular matrix. The mechanical response of these systems when subjected to an applied strain at zero temperature is often investigated numerically using networks composed of filaments, which are either rigidly welded or pinned together at their crosslinks. In the latter, filaments during deformation are free to rotate about their crosslinks while the relative angles between filaments remain constant in the former. The behavior of crosslinks in actual semiflexible networks is different than these idealized models and there exists only partial constraint on torques at crosslinks. The present work develops a numerical model in which two intersecting filaments are connected to each other by torsional springs with arbitrary stiffness. We show that fiber networks composed of rigid and freely rotating crosslinks are the limiting case of the present model. Furthermore, we characterize the effects of stiffness of crosslinks on effective Young's modulus of semiflexible networks as a function of filament flexibility and crosslink density. The effective Young's modulus is determined as a function of the mechanical properties of crosslinks and is found to vanish for networks composed of very weak torsional springs. Independent of the stiffness of crosslinks, it is found that the effective Young's modulus is a function of fiber flexibility and crosslink density. In low density networks, filaments primarily bend and the effective Young's modulus is much lower than the affine estimate. With increasing filament bending stiffness and/or crosslink density, the mechanical behavior of the networks becomes more affine and the stretching of filaments depicts itself as the dominant mode of deformation. The torsional stiffness of the crosslinks significantly affects the effective Young's modulus of the semiflexible random fiber networks.

Crosslinking of fibrous hydrogels

In contrast to most synthetic hydrogels, biological gels are made of fibrous networks. This architecture gives rise to unique properties, like low concentration, high porosity gels with a high mechanical responsiveness as a result of strain-stiffening. Here, we used a synthetic polymer model system, based on polyisocyanides, that we crosslinked selectively inside the bundles. This approach allows us to lock in the fibrous network present at the crosslinking conditions. At minimum crosslink densities, we are able to freeze in the architecture, as well as the associated mechanical properties. Rheology and X-ray scattering experiments show that we able to accurately tailor network mechanics, not by changing the gel composition or architecture, but rather by tuning its (thermal) history. Selective crosslinking is a crucial step in making biomimetic networks with a controlled architecture.

Unlike synthetic hydrogels, biological gels are made of fibrous networks which give rise to unique properties, such as high porosity and mechanical responsiveness. Here the authors use polyisocyanide-based gels and selectively crosslink inside the bundles to lock the fibrous network and thus control the architecture and the mechanics.

Topology determines force distributions in one-dimensional random spring networks

Networks of elastic fibers are ubiquitous in biological systems and often provide mechanical stability to cells and tissues. Fiber-reinforced materials are also common in technology. An important characteristic of such materials is their resistance to failure under load. Rupture occurs when fibers break under excessive force and when that failure propagates. Therefore, it is crucial to understand force distributions. Force distributions within such networks are typically highly inhomogeneous and are not well understood. Here we construct a simple one-dimensional model system with periodic boundary conditions by randomly placing linear springs on a circle. We consider ensembles of such networks that consist of N nodes and have an average degree of connectivity z but vary in topology. Using a graph-theoretical approach that accounts for the full topology of each network in the ensemble, we show that, surprisingly, the force distributions can be fully characterized in terms of the parameters (N,z). Despite the universal properties of such (N,z) ensembles, our analysis further reveals that a classical mean-field approach fails to capture force distributions correctly. We demonstrate that network topology is a crucial determinant of force distributions in elastic spring networks.

Topology Counts: Force Distributions in Circular Spring Networks

Filamentous polymer networks govern the mechanical properties of many biological materials. Force distributions within these networks are typically highly inhomogeneous, and, although the importance of force distributions for structural properties is well recognized, they are far from being understood quantitatively. Using a combination of probabilistic and graph-theoretical techniques, we derive force distributions in a model system consisting of ensembles of random linear spring networks on a circle. We show that characteristic quantities, such as the mean and variance of the force supported by individual springs, can be derived explicitly in terms of only two parameters: (i) average connectivity and (ii) number of nodes. Our analysis shows that a classical mean-field approach fails to capture these characteristic quantities correctly. In contrast, we demonstrate that network topology is a crucial determinant of force distributions in an elastic spring network. Our results for 1D linear spring networks readily generalize to arbitrary dimensions.

A new model linking elastic properties and ionic conductivity of mixed network former glasses

A new statistical thermodynamic model has been developed to describe the activated process of cation hopping in mixed network former glasses based on the systematic comparison between the adiabatic elastic moduli measured using Brillouin light scattering and the ionic conductivity measured using dielectric impedance spectroscopy.

Force percolation of contractile active gels

Living systems provide a paradigmatic example of active soft matter. Cells and tissues comprise viscoelastic materials that exert forces and can actively change shape. This strikingly autonomous behavior is powered by the cytoskeleton, an active gel of semiflexible filaments, crosslinks, and molecular motors inside cells. Although individual motors are only a few nm in size and exert minute forces of a few pN, cells spatially integrate the activity of an ensemble of motors to produce larger contractile forces (∼nN and greater) on cellular, tissue, and organismal length scales. Here we review experimental and theoretical studies on contractile active gels composed of actin filaments and myosin motors. Unlike other active soft matter systems, which tend to form ordered patterns, actin-myosin systems exhibit a generic tendency to contract. Experimental studies of reconstituted actin-myosin model systems have long suggested that a mechanical interplay between motor activity and the network's connectivity governs this contractile behavior. Recent theoretical models indicate that this interplay can be understood in terms of percolation models, extended to include effects of motor activity on the network connectivity. Based on concepts from percolation theory, we propose a state diagram that unites a large body of experimental observations. This framework provides valuable insights into the mechanisms that drive cellular shape changes and also provides design principles for synthetic active materials.

Elastic properties and short-range structural order in mixed network former glasses

A new statistical thermodynamic model has been developed to describe the speciation of network former elements in ternary oxide glasses, which uses data from NMR spectroscopy and the adiabatic elastic moduli measured using Brillouin light scattering as input.

Strain-driven criticality underlies nonlinear mechanics of fibrous networks

Networks with only central force interactions are floppy when their average connectivity is below an isostatic threshold. Although such networks are mechanically unstable, they can become rigid when strained. It was recently shown that the transition from floppy to rigid states as a function of simple shear strain is continuous, with hallmark signatures of criticality [Sharma et al., Nature Phys. 12, 584 (2016)1745-247310.1038/nphys3628]. The nonlinear mechanical response of collagen networks was shown to be quantitatively described within the framework of such mechanical critical phenomenon. Here, we provide a more quantitative characterization of critical behavior in subisostatic networks. Using finite-size scaling we demonstrate the divergence of strain fluctuations in the network at well-defined critical strain. We show that the characteristic strain corresponding to the onset of strain stiffening is distinct from but related to this critical strain in a way that depends on critical exponents. We confirm this prediction experimentally for collagen networks. Moreover, we find that the apparent critical exponents are largely independent of the spatial dimensionality. With subisostaticity as the only required condition, strain-driven criticality is expected to be a general feature of biologically relevant fibrous networks.

Sulfo-SMCC Prevents Annealing of Taxol-Stabilized Microtubules In Vitro

Microtubule structure and functions have been widely studied in vitro and in cells. Research has shown that cysteines on tubulin play a crucial role in the polymerization of microtubules. Here, we show that blocking sulfhydryl groups of cysteines in taxol-stabilized polymerized microtubules with a commonly used chemical crosslinker prevents temporal end-to-end annealing of microtubules in vitro. This can dramatically affect the length distribution of the microtubules. The crosslinker sulfosuccinimidyl 4-(N-maleimidomethyl)cyclohexane-1-carboxylate, sulfo-SMCC, consists of a maleimide and an N-hydroxysuccinimide ester group to bind to sulfhydryl groups and primary amines, respectively. Interestingly, addition of a maleimide dye alone does not show the same interference with annealing in stabilized microtubules. This study shows that the sulfhydryl groups of cysteines of tubulin that are vital for the polymerization are also important for the subsequent annealing of microtubules.

Elasticity of fibrous networks under uniaxial prestress

We present theoretical and experimental studies of the elastic response of fibrous networks subjected to uniaxial strain. Uniaxial compression or extension is applied to extracellular networks of fibrin and collagen using a shear rheometer with free water in/outflow. Both uniaxial stress and the network shear modulus are measured. Prior work [van Oosten, et al., Sci. Rep., 2015, 6, 19270] has shown softening/stiffening of these networks under compression/extension, together with a nonlinear response to shear, but the origin of such behaviour remains poorly understood. Here, we study how uniaxial strain influences the nonlinear mechanics of fibrous networks. Using a computational network model with bendable and stretchable fibres, we show that the softening/stiffening behaviour can be understood for fixed lateral boundaries in 2D and 3D networks with comparable average connectivities to the experimental extracellular networks. Moreover, we show that the onset of stiffening depends strongly on the imposed uniaxial strain. Our study highlights the importance of both uniaxial strain and boundary conditions in determining the mechanical response of hydrogels.

Nonlinear Mechanics of Athermal Branched Biopolymer Networks

Naturally occurring biopolymers such as collagen and actin form branched fibrous networks. The average connectivity in branched networks is generally below the isostatic threshold at which central force interactions marginally stabilize the network. In the submarginal regime, for connectivity below this threshold, such networks are unstable toward small deformations unless stabilized by additional interactions such as bending. Here we perform a numerical study on the elastic behavior of such networks. We show that the nonlinear mechanics of branched networks is qualitatively similar to that of filamentous networks with freely hinged cross-links. In agreement with a recent theoretical study,1 we find that branched networks also exhibit nonlinear mechanics consistent with athermal critical phenomena controlled by strain. We obtain the critical exponents capturing the nonlinear elastic behavior near the critical point by performing scaling analysis of the stiffening curves. We find that the exponents evolve with the connectivity in the network. We show that the nonlinear mechanics of disordered networks, independent of the detailed microstructure, can be characterized by a strain-driven second-order phase transition, and that the primary quantitative differences among different architectures are in the critical exponents describing the transition.

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.