Stretching Semiflexible Filaments and Their Networks

Microstructure-based nuclear lamina constitutive model

The nuclear lamina is widely recognized as the most crucial component in providing mechanical stability to the nucleus. However, it is still a significant challenge to model the mechanics of this multilayered protein network. We developed a constitutive model of the nuclear lamina network based on its microstructure, which accounts for the deformation phases at the dimer level, as well as the orientational arrangement and density of lamin filaments. Instead of relying on homology modeling in the previous studies, we conducted molecular simulations to predict the force-extension response of a highly accurate lamin dimer structure obtained through X-ray diffraction crystallography experimentation. Furthermore, we devised a semiflexible worm-like chain extension-force model of lamin dimer as a substitute, incorporating phases of initial stretching, uncoiling of the dimer coiled-coil, and transition of secondary structures. Subsequently, we developed a 2D network continuum model for the nuclear lamina by using our extension-force lamin dimer model and derived stress resultants. By comparing with experimentally measured lamina modulus, we found that the lamina network has sharp initial strain-hardening behavior. This also enabled us to carry out finite element simulations of the entire nucleus with an accurate microstructure-based nuclear lamina model. Finally, we conducted simulations of transendothelial transmigration of a nucleus and investigated the impact of varying network density and uncoiling constants on the critical force required for successful transmigration. The model allows us to incorporate the microstructure characteristics of the nuclear lamina into the nucleus model, thereby gaining insights into how laminopathies and mutations affect nuclear mechanics.

Sex Differences in Axonal Dynamic Responses Under Realistic Tension Using Finite Element Models

Existing axonal finite element models do not consider sex morphological differences or the fidelity in dynamic input. To facilitate a systematic investigation into the micromechanics of diffuse axonal injury, we develop a parameterized modeling approach for automatic and efficient generation of sex-specific axonal models according to specified geometrical parameters. Baseline female and male axonal models in the corpus callosum with random microtubule (MT) gap configurations are generated for model calibration and evaluation. They are then used to simulate a realistic tensile loading consisting of both a loading and a recovery phase (to return to an initial undeformed state) generated from dynamic corpus callosum fiber strain in a real-world head impact simulation. We find that MT gaps and the dynamic recovery phase are both critical to successfully reproduce MT undulation as observed experimentally, which has not been reported before. This strengthens confidence in model dynamic responses. A statistical approach is further employed to aggregate axonal responses from a large sample of random MT gap configurations for both female and male axonal models (n = 10,000 each). We find that peak strains in MTs and the Ranvier node and associated neurofilament failures in female axons are substantially higher than those in male axons because there are fewer MTs in the former and also because of the random nature of MT gap locations. Despite limitations in various model assumptions as a result of limited experimental data currently available, these findings highlight the need to systematically characterize MT gap configurations and to ensure a realistic model input for axonal dynamic simulations. Finally, this study may offer fresh and improved insight into the biomechanical basis of sex differences in brain injury, and sets the stage for more systematic investigations at the microscale in the future, both numerically and experimentally.

Mechanics of 3D Cell-Hydrogel Interactions: Experiments, Models, and Mechanisms

Hydrogels are highly water-swollen molecular networks that are ideal platforms to create tissue mimetics owing to their vast and tunable properties. As such, hydrogels are promising cell-delivery vehicles for applications in tissue engineering and have also emerged as an important base for ex vivo models to study healthy and pathophysiological events in a carefully controlled three-dimensional environment. Cells are readily encapsulated in hydrogels resulting in a plethora of biochemical and mechanical communication mechanisms, which recapitulates the natural cell and extracellular matrix interaction in tissues. These interactions are complex, with multiple events that are invariably coupled and spanning multiple length and time scales. To study and identify the underlying mechanisms involved, an integrated experimental and computational approach is ideally needed. This review discusses the state of our knowledge on cell-hydrogel interactions, with a focus on mechanics and transport, and in this context, highlights recent advancements in experiments, mathematical and computational modeling. The review begins with a background on the thermodynamics and physics fundamentals that govern hydrogel mechanics and transport. The review focuses on two main classes of hydrogels, described as semiflexible polymer networks that represent physically cross-linked fibrous hydrogels and flexible polymer networks representing the chemically cross-linked synthetic and natural hydrogels. In this review, we highlight five main cell-hydrogel interactions that involve key cellular functions related to communication, mechanosensing, migration, growth, and tissue deposition and elaboration. For each of these cellular functions, recent experiments and the most up to date modeling strategies are discussed and then followed by a summary of how to tune hydrogel properties to achieve a desired functional cellular outcome. We conclude with a summary linking these advancements and make the case for the need to integrate experiments and modeling to advance our fundamental understanding of cell-matrix interactions that will ultimately help identify new therapeutic approaches and enable successful tissue engineering.

Geometrical nonlinear elasticity of axon under tension: A coarse-grained computational study

Axon bundles cross-linked by microtubule (MT) associate proteins and bounded by a shell skeleton are critical for normal function of neurons. Understanding effects of the complexly geometrical parameters on their mechanical properties can help gain a biomechanical perspective on the neurological functions of axons and thus brain disorders caused by the structural failure of axons. Here, the tensile mechanical properties of MT bundles cross-linked by tau proteins are investigated by systematically tuning MT length, axonal cross-section radius, and tau protein spacing in a bead-spring coarse-grained model. Our results indicate that the stress-strain curves of axons can be divided into two regimes, a nonlinear elastic regime dominated by rigid-body like inter-MT sliding, and a linear elastic regime dominated by affine deformation of both tau proteins and MTs. From the energetic analyses, first, the tau proteins dominate the mechanical performance of axons under tension. In the nonlinear regime, tau proteins undergo a rigid-body like rotating motion rather than elongating, whereas in the nonlinear elastic regime, tau proteins undergo a flexible elongating deformation along the MT axis. Second, as the average spacing between adjacent tau proteins along the MT axial direction increases from 25 to 125 nm, the Young's modulus of axon experiences a linear decrease whereas with the average space varying from 125 to 175 nm, and later reaches a plateau value with a stable fluctuation. Third, the increment of the cross-section radius of the MT bundle leads to a decrease in Young's modulus of axon, which is possibly attributed to the decrease in MT numbers per cross section. Overall, our research findings offer a new perspective into understanding the effects of geometrical parameters on the mechanics of MT bundles as well as serving as a theoretical basis for the development of artificial MT complexes potentially toward medical applications.

Single-molecule stretching experiments of flexible (wormlike) chain molecules in different ensembles: Theory and a potential application of finite chain length effects to nick-counting in DNA

We propose a formalism for deriving force-elongation and elongation-force relations for flexible chain molecules from analytical expressions for their radial distribution function, which provides insight into the factors controlling the asymptotic behavior and finite chain length corrections. In particular, we apply this formalism to our previously developed interpolation formula for the wormlike chain end-to-end distance distribution. The resulting expression for the asymptotic limit of infinite chain length is of similar quality to the numerical evaluation of Marko and Siggia's variational theory and considerably more precise than their interpolation formula. A comparison to numerical data suggests that our analytical finite chain length corrections achieve a comparable accuracy. As an application of our results, we discuss the possibility of inferring the time-dependent number of nicks in single-molecule stretching experiments on double-stranded DNA from the accompanying changes in the effective chain length.

A NanoFE simulation-based surrogate machine learning model to predict mechanical functionality of protein networks from live confocal imaging

Sub-cellular mechanics plays a crucial role in a variety of biological functions and dysfunctions. Due to the strong structure-function relationship in cytoskeletal protein networks, light can be shed on their mechanical functionality by investigating their structures. Here, we present a data-driven approach employing a combination of confocal live imaging of fluorescent tagged protein networks, in silico mechanical experiments and machine learning to investigate this relationship. Our designed image processing and nanoFE mechanical simulation framework resolves the structure and mechanical behaviour of cytoskeletal networks and the developed gradient boosting surrogate models linking network structure to its functionality. In this study, for the first time, we perform mechanical simulations of Filamentous Temperature Sensitive Z (FtsZ) complex protein networks with realistic network geometry depicting its skeletal functionality inside organelles (here, chloroplasts) of the moss Physcomitrella patens. Training on synthetically produced simulation data enables predicting the mechanical characteristics of FtsZ network purely based on its structural features (R2⩾0.93), therefore allowing to extract structural principles enabling specific mechanical traits of FtsZ, such as load bearing and resistance to buckling failure in case of large network deformation.

A mesoscopic theory to describe the flexibility regulation in F-actin networks: An approach of phase transitions with nonlinear elasticity

The synthetic actin network arouses great interest as bio-material due to its soft and wet nature that mimics many biological scaffolding structures. Inside the cell, the actin network is regulated by tens of actin-binding proteins (ABP's), which make for a highly complex system with several emergent behaviors. In particular, calponin is an ABP that was identified as an actin stabiliser, but whose mechanism is still poorly understood. Recent experiments using an in vitro model system of cross-linked actin with calponin and large deformation bulk rheology, found that networks with exhibited a delayed onset and were able to withstand a higher maximal strain before softening. In this work, we show that at network scale the actin network with calponin furthermore the reduction of the persistence length allows: (i) The reduction in the network pre-strain. (ii) The increment of the crosslinks adhesion energy. We verify these effects theoretically using nonlinear continuum mechanics for the semiflexible and crosslinked network. In addition, the alterations over the microstructure are described by the definition of an interaction parameter Γ according the formalism of Landau for phase transitions. According to this model we demonstrates that the interaction parameter can describe the experimental observations following a scaling exponent as Γ~|c-c_cr|^1/2, where c is the ratio between concentration of calponin and actin. This result provides interesting feedback to improve our understanding of several mechano-biological pathways.

Effects of cosolvent partitioning on conformational transitions and chain flexibility of thermoresponsive microgels

The conformational collapse of polymers in mixtures of two individually good solvents is an intriguing yet puzzling phenomenon termed cononsolvency. In this paper, the concept of the preferential adsorption of the cosolvent is combined with mean-field approaches to elaborate the cononsolvency effect of dimethylformamide (DMF) on the thermoresponsive poly(N-isopropylacrylamide) (PNIPAM) microgels in aqueous solutions. We give a quantitative description concerning the effects of DMF preferential adsorption and partitioning on the reentrant transition of PNIPAM microgels below the lower critical solution temperature (LCST) of PNIPAM. While the DMF cononsolvency incurs the conformational collapse, the affinity of DMF molecules to PNIPAM chains becomes increasingly stronger, which reveals that the conformational collapse is decoupled from the solvent quality of DMF-water mixtures. Considering the chain elasticity, spatial constraints, and surface charge of microgels, we explore the cononsolvency effect on the persistence length quantifying the PNIPAM flexibility. Our analysis elucidates that, depending on chain length and temperature, the DMF cononsolvency-induced collapse of PNIPAM microgels leads to a remarkable increase in the persistent length below LCST, which is comparable to the experimental data regarding suspension mechanical properties of PNIPAM microgels in water above LCST.

Multiscale FEM simulations of cross-linked actin network embedded in cytosol with the focus on the filament orientation

The present contribution focuses on the application of the multiscale finite element method to the modeling of actin networks that are embedded in the cytosol. These cell components are of particular importance with regard to the cell response to external stimuli. The homogenization strategy chosen uses the Hill-Mandel macrohomogeneity condition for bridging 2 scales: the macroscopic scale that is related to the cell level and the microscopic scale related to the representative volume element. For the modeling of filaments, the Holzapfel-Ogden β-model is applied. It provides a relationship between the tensile force and the caused stretches, serves as the basis for the derivation of the stress and elasticity tensors, and enables a novel finite element implementation. The elements with the neo-Hookean constitutive law are applied for the simulation of the cytosol. The results presented corroborate the main advantage of the concept, namely, its flexibility with regard to the choice of the representative volume element as well as of macroscopic tests. The focus is particularly placed on the study of the filament orientation and of its influence on the effective behavior.

Concentration-dependent swelling and structure of ionic microgels: simulation and theory of a coarse-grained model

We study swelling and structural properties of ionic microgel suspensions within a comprehensive coarse-grained model that combines the polymeric and colloidal natures of microgels as permeable, compressible, charged spheres governed by effective interparticle interactions. The model synthesizes the Flory-Rehner theory of cross-linked polymer gels, the Hertz continuum theory of effective elastic interactions, and a theory of density-dependent effective electrostatic interactions. Implementing the model using Monte Carlo simulation and thermodynamic perturbation theory, we compute equilibrium particle size distributions, swelling ratios, volume fractions, net valences, radial distribution functions, and static structure factors as functions of concentration. Trial Monte Carlo moves comprising particle displacements and size variations are accepted or rejected based on the total change in elastic and electrostatic energies. The theory combines first-order thermodynamic perturbation and variational free energy approximations. For illustrative system parameters, theory and simulation agree closely at concentrations ranging from dilute to beyond particle overlap. With increasing concentration, as microgels deswell, we predict a decrease in the net valence and an unusual saturation of pair correlations. Comparison with experimental data for deionized, aqueous suspensions of PNIPAM particles demonstrates the capacity of the coarse-grained model to predict and interpret measured swelling behavior.

Gels with sense: supramolecular materials that respond to heat, light and sound

Advances in the field of supramolecular chemistry have made it possible, in many situations, to reliably engineer soft materials to address a specific technological problem. Particularly exciting are "smart" gels that undergo reversible physical changes on exposure to remote, non-invasive environmental stimuli. This review explores the development of gels which are transformed by heat, light and ultrasound, as well as other mechanical inputs, applied voltages and magnetic fields. Focusing on small-molecule gelators, but with reference to organic polymers and metal-organic systems, we examine how the structures of gelator assemblies influence the physical and chemical mechanisms leading to thermo-, photo- and mechano-switchable behaviour. In addition, we evaluate how the unique and versatile properties of smart materials may be exploited in a wide range of applications, including catalysis, crystal growth, ion sensing, drug delivery, data storage and biomaterial replacement.

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.

Effects of disorder and chain stiffening on the elasticity of flexible polymer networks

We examine how the distribution of contour lengths and the high-stretch stiffening of individual chain segments affect the macroscopic shear modulus of flexible polymer gels, using a two-dimensional numerical model in which polymer segments form a triangular network and disorder is introduced by varying their contour lengths. We show that, in the relevant parameter range: (i) the nonaffine contribution to the shear modulus is negligible, i.e., the Born approximation is satisfactory, and (ii) the shear modulus is dominated by the contribution originating from equilibrium chain tensions. Moreover, mechanical equilibration at the nodes induces specific correlations between the end-to-end distances and contour lengths of chain segments, which must be properly accounted for to construct reasonable estimates of chain pressure and shear modulus.

Theory of Semiflexible Filaments and Networks

We briefly review the recent developments in the theory of individual semiflexible filaments, and of a crosslinked network of such filaments, both permanent and transient. Starting from the free energy of an individual semiflexible chain, models on its force-extension relation and other mechanical properties such as Euler buckling are discussed. For a permanently crosslinked network of filaments, theories on how the network responds to deformation are provided, with a focus on continuum approaches. Characteristic features of filament networks, such as nonlinear stress-strain relation, negative normal stress, tensegrity, and marginal stability are discussed. In the new area of transient filament network, where the crosslinks can be dynamically broken and re-formed, we show some recent attempts for understanding the dynamics of the crosslinks, and the related rheological properties, such as stress relaxation, yield stress and plasticity.

Conformational Properties of Active Semiflexible Polymers

The conformational properties of flexible and semiflexible polymers exposed to active noise are studied theoretically. The noise may originate from the interaction of the polymer with surrounding active (Brownian) particles or from the inherent motion of the polymer itself, which may be composed of active Brownian particles. In the latter case, the respective monomers are independently propelled in directions changing diffusively. For the description of the polymer, we adopt the continuous Gaussian semiflexible polymer model. Specifically, the finite polymer extensibility is taken into account, which turns out to be essential for the polymer conformations. Our analytical calculations predict a strong dependence of the relaxation times on the activity. In particular, semiflexible polymers exhibit a crossover from a bending elasticity-dominated dynamics to the flexible polymer dynamics with increasing activity. This leads to a significant activity-induced polymer shrinkage over a large range of self-propulsion velocities. For large activities, the polymers swell and their extension becomes comparable to the contour length. The scaling properties of the mean square end-to-end distance with respect to the polymer length and monomer activity are discussed.

Nonlinear elasticity of semiflexible filament networks

We develop a continuum theory for equilibrium elasticity of a network of crosslinked semiflexible filaments, spanning the full range between flexible entropy-driven chains to stiff athermal rods. We choose the 3-chain constitutive model of network elasticity over several plausible candidates, and derive analytical expressions for the elastic energy at arbitrary strain, with the corresponding stress-strain relationship. The theory fits well to a wide range of experimental data on simple shear in different filament networks, quantitatively matching the differential shear modulus variation with stress, with only two adjustable parameters (which represent the filament stiffness and the pre-tension in the network, respectively). The general theory accurately describes the crossover between the positive and negative Poynting effect (normal stress on imposed shear) on increasing the stiffness of filaments forming the network. We discuss the network stability (the point of marginal rigidity) and the phenomenon of tensegrity, showing that filament pre-tension on crosslinking into the network determines the magnitude of linear modulus G0.

F1 rotary motor of ATP synthase is driven by the torsionally-asymmetric drive shaft

F₁F₀ ATP synthase (ATPase) either facilitates the synthesis of ATP in a process driven by the proton moving force (pmf), or uses the energy from ATP hydrolysis to pump protons against the concentration gradient across the membrane. ATPase is composed of two rotary motors, F₀ and F₁, which compete for control of their shared γ -shaft. We present a self-consistent physical model of F₁ motor as a simplified two-state Brownian ratchet using the asymmetry of torsional elastic energy of the coiled-coil γ -shaft. This stochastic model unifies the physical concepts of linear and rotary motors, and explains the stepped unidirectional rotary motion. Substituting the model parameters, all independently known from recent experiments, our model quantitatively reproduces the ATPase operation, e.g. the ‘no-load’ angular velocity is ca. 400 rad/s anticlockwise at 4 mM ATP. Increasing the pmf torque exerted by F₀ can slow, stop and overcome the torque generated by F₁, switching from ATP hydrolysis to synthesis at a very low value of ‘stall torque’. We discuss the motor efficiency, which is very low if calculated from the useful mechanical work it produces - but is quite high when the ‘useful outcome’ is measured in the number of H⁺ pushed against the chemical gradient.

Torsional behavior of axonal microtubule bundles

Axonal microtubule (MT) bundles crosslinked by microtubule-associated protein (MAP) tau are responsible for vital biological functions such as maintaining mechanical integrity and shape of the axon as well as facilitating axonal transport. Breaking and twisting of MTs have been previously observed in damaged undulated axons. Such breaking and twisting of MTs is suggested to cause axonal swellings that lead to axonal degeneration, which is known as "diffuse axonal injury". In particular, overstretching and torsion of axons can potentially damage the axonal cytoskeleton. Following our previous studies on mechanical response of axonal MT bundles under uniaxial tension and compression, this work seeks to characterize the mechanical behavior of MT bundles under pure torsion as well as a combination of torsional and tensile loads using a coarse-grained computational model. In the case of pure torsion, a competition between MAP tau tensile and MT bending energies is observed. After three turns, a transition occurs in the mechanical behavior of the bundle that is characterized by its diameter shrinkage. Furthermore, crosslink spacing is shown to considerably influence the mechanical response, with larger MAP tau spacing resulting in a higher rate of turns. Therefore, MAP tau crosslinking of MT filaments protects the bundle from excessive deformation. Simultaneous application of torsion and tension on MT bundles is shown to accelerate bundle failure, compared to pure tension experiments. MAP tau proteins fail in clusters of 10-100 elements located at the discontinuities or the ends of MT filaments. This failure occurs in a stepwise fashion, implying gradual accumulation of elastic tensile energy in crosslinks followed by rupture. Failure of large groups of interconnecting MAP tau proteins leads to detachment of MT filaments from the bundle near discontinuities. This study highlights the importance of torsional loading in axonal damage after traumatic brain injury.

Cavity approach for modeling and fitting polymer stretching

The mechanical properties of molecules are today captured by single molecule manipulation experiments, so that polymer features are tested at a nanometric scale. Yet devising mathematical models to get further insight beyond the commonly studied force-elongation relation is typically hard. Here we draw from techniques developed in the context of disordered systems to solve models for single and double-stranded DNA stretching in the limit of a long polymeric chain. Since we directly derive the marginals for the molecule local orientation, our approach allows us to readily calculate the experimental elongation as well as other observables at wish. As an example, we evaluate the correlation length as a function of the stretching force. Furthermore, we are able to fit successfully our solution to real experimental data. Although the model is admittedly phenomenological, our findings are very sound. For single-stranded DNA our solution yields the correct (monomer) scale and yet, more importantly, the right persistence length of the molecule. In the double-stranded case, our model reproduces the well-known overstretching transition and correctly captures the ratio between native DNA and overstretched DNA. Also in this case the model yields a persistence length in good agreement with consensus, and it gives interesting insights into the bending stiffness of the native and overstretched molecule, respectively.

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.

Mechanics of biological networks: from the cell cytoskeleton to connective tissue

From the cell cytoskeleton to connective tissues, fibrous networks are ubiquitous in metazoan life as the key promoters of mechanical strength, support and integrity. In recent decades, the application of physics to biological systems has made substantial strides in elucidating the striking mechanical phenomena observed in such networks, explaining strain stiffening, power law rheology and cytoskeletal fluidisation - all key to the biological function of individual cells and tissues. In this review we focus on the current progress in the field, with a primer into the basic physics of individual filaments and the networks they form. This is followed by a discussion of biological networks in the context of a broad spread of recent in vitro and in vivo experiments.

Wrinkling of random and regular semiflexible polymer networks

We investigate wrinkling of two-dimensional random and triangular semiflexible polymer networks under shear. Both types of semiflexible networks exhibit wrinkling above a small critical shear angle, which scales with an exponent of the bending modulus between 1.9 and 2.0. Random networks exhibit hysteresis at the wrinkling threshold. Wrinkling lowers the total elastic energy by up to 20% and strongly affects the elastic properties of all semiflexible networks such as the crossover between bending and stretching dominated behavior. In random networks, we also find evidence for metastable wrinkled configurations. While the disordered microstructure of random networks affects the scaling behavior of wrinkle amplitudes, it has little effect on wrinkle wavelength. Therefore, wrinkles represent a robust, microstructure-independent assay of shear strain or elastic properties.

Viscoelasticity of cross-linked actin networks: experimental tests, mechanical modeling and finite-element analysis

Filamentous actin is one of the main constituents of the eukaryotic cytoskeleton. The actin cortex, a densely cross-linked network, resides underneath the lipid bilayer. In the present work we propose a continuum mechanical formulation for describing the viscoelastic properties of in vitro actin networks, which serve as model systems for the cortex, by including the microstructure, i.e. the behavior of a single filament and its spatial arrangement. The modeling of the viscoelastic response in terms of physically interpretable parameters is conducted using a multiscale approach consisting of two steps: modeling of the single filament response of F-actin by a worm-like chain model including the extensibility of the filament, and assembling the three-dimensional biopolymer network by using the microsphere model which accounts for filaments equally distributed in space. The viscoelastic effects of the network are taken into account using a generalized Maxwell model. The Cauchy stress and elasticity tensors are obtained within a continuum mechanics framework and implemented into a finite-element program. The model is validated on the network level using large strain experiments on reconstituted actin gels. Comparisons of the proposed model with rheological experiments recover reasonable values for the material parameters. Finite-element simulations of the indentation of a sphere on a network slab and the aspiration of a droplet in a micropipette allow for further insights of the viscoelastic behavior of actin networks.

Elasticity of biopolymer filaments

Within the general one-dimensional theory of nonlinear elasticity we analyze the elasticity of biopolymer filaments. The approach adopted is purely mechanical but is reconciled with statistical physics approaches and allows for a proper formulation of boundary-value problems. By specializing the general framework we obtain force-extension relations for inextensible filaments and show how previous work on the biophysics of filaments fits within the framework. On the other hand, within the same framework, the theory of extensible filaments, which is appropriate for semi-flexible filaments such as F-actin, enables us to fit representative F-actin data. The specific formulas derived are relatively simple and the parameters involved have direct mechanical interpretations and are immediately related to the filament properties, including the initial end-to-end length, contour length and persistence length.

Microrheology of highly crosslinked microtubule networks is dominated by force-induced crosslinker unbinding

We determine the time- and force-dependent viscoelastic responses of reconstituted networks of microtubules that have been strongly crosslinked by biotin-streptavidin bonds. To measure the microscale viscoelasticity of such networks, we use a magnetic tweezers device to apply localized forces. At short time scales, the networks respond nonlinearly to applied force, with stiffening at small forces, followed by a reduction in the stiffening response at high forces, which we attribute to the force-induced unbinding of crosslinks. At long time scales, force-induced bond unbinding leads to local network rearrangement and significant bead creep. Interestingly, the network retains its elastic modulus even under conditions of significant plastic flow, suggesting that crosslinker breakage is balanced by the formation of new bonds. To better understand this effect, we developed a finite element model of such a stiff filament network with labile crosslinkers obeying force-dependent Bell model unbinding dynamics. The coexistence of dissipation, due to bond breakage, and the elastic recovery of the network is possible because each filament has many crosslinkers. Recovery can occur as long as a sufficient number of the original crosslinkers are preserved under the loading period. When these remaining original crosslinkers are broken, plastic flow results.

A new approach to model cross-linked actin networks: multi-scale continuum formulation and computational analysis

The mechanical properties of a cell are defined mainly by the cytoskeleton. One contributor within this three-dimensional structure is the actin cortex which is located underneath the lipid bilayer. It forms a nearly isotropic and densely cross-linked protein network. We present a continuum mechanical formulation for describing the mechanical properties of in vitro model systems based on their micro-structure, i.e. the behavior of a single filament and its spatial arrangement. The network is considered elastic, viscous effects being neglected. Filamentous actin is a biopolymer with a highly nonlinear force-stretch relationship. This can be well described by a worm-like chain model that includes extensibility of the filament, which we call the β-model. A comparison with experimental data shows good agreement with values for the physically interpretable parameters. To make these properties applicable to three dimensions we used a non-affine micro-sphere network, which accounts for filaments, equally distributed in space. The assembled model results in a strain-energy density which is a function of the deformation gradient, and it is validated with experimental data from rheological experiments of in vitro reconstituted actin networks. The Cauchy stress and elasticity tensors are obtained within the continuum mechanics framework and implemented into a finite element program to solve boundary-value problems.

Computational modeling of axonal microtubule bundles under tension

Microtubule bundles cross-linked by tau protein serve a variety of neurological functions including maintaining mechanical integrity of the axon, promoting axonal growth, and facilitating cargo transport. It has been observed that axonal damage in traumatic brain injury leads to bundle disorientation, loss of axonal viability, and cognitive impairment. This study investigates the initial mechanical response of axonal microtubule bundles under uniaxial tension using a discrete bead-spring representation. Mechanisms of failure due to traumatic stretch loading and their impact on the mechanical response and stability are also characterized. This study indicates that cross-linked axonal microtubule bundles in tension display stiffening behavior similar to a power-law relationship from nonaffine network deformations. Stretching of cross-links and microtubule bending were the primary deformation modes at low stresses. Microtubule stretch was negligible up to tensile stresses of ∼1 MPa. Bundle failure occurred by failure of cross-links leading to pull-out of microtubules and loss of bundle integrity. This may explain the elongation, undulation, and delayed elasticity of axons following traumatic stretch loading. More extensively cross-linked bundles withstood higher tensile stresses before failing. The bundle mechanical behavior uncovered by these computational techniques should guide future experiments on stretch-injured axons.

Polymer physics of the cytoskeleton

The cytoskeleton is generally visualized by light or electron microscopy as a meshwork of protein filaments that spans the space between the nuclear envelope and the plasma membrane. In most cell types, this meshwork is formed by a three dimensional composite network of actin filaments, microtubules (MT), and intermediate filaments (IF) together with the host of proteins that bind to the sides or ends of these linear polymers. Cytoskeletal binding proteins regulate filament length, crosslink filaments to each other, and apply forces to the filaments. One approach to modeling the mechanical properties of the cytoskeleton and of cell in general is to consider the elements of the cytoskeleton as polymers, using experimental methods and theoretical models developed for traditional polymers but modified for the much larger, stiffer, and fragile biopolymers comprising the cytoskeleton. The presence of motor proteins that move actin filaments and microtubules also creates a new class of active materials that are out of thermodynamic equilibrium, and unconstrained by limitations of the fluctuation-dissipation theorem. These active materials create rich opportunities for experimental design and theoretical developments. The degree to which the mechanics of live cells can usefully be modeled as highly complex polymer networks is by no means certain, and this article will discuss recent progress in quantitatively measuring cytoskeletal polymer systems and relating them to the properties of the cell.

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.

Stretching weakly bending filaments with spontaneous curvature in two dimensions

Some important biomolecules (for instance, bacterial FtsZ and eukaryotic DNA) are known to posses spontaneous (intrinsic) curvature. Using a simple extension of the wormlike chain model, we study the response of a weakly bending filament in two dimensions to a pulling force applied at its ends (a configuration common in classical in-vitro experiments and relevant to several in-vivo cell cases). The spontaneous curvature of such a chain or filament can in general be arc-length dependent and we study a case of sinusoidal variation, from which an arbitrary case can be reconstructed via Fourier transformation. We obtain analytic results for the force-extension relationship and the width of transverse fluctuations. We show that spontaneous-curvature undulations can affect the force-extension behavior even in relatively flexible filaments with a persistence length smaller than the contour length.