Elasticity of biopolymer filaments

Nonlinear Poisson effect in affine semiflexible polymer networks

Stretching an elastic material along one axis typically induces contraction along the transverse axes, a phenomenon known as the Poisson effect. From these strains, one can compute the specific volume, which generally either increases or, in the incompressible limit, remains constant as the material is stretched. However, in networks of semiflexible or stiff polymers, which are typically highly compressible yet stiffen significantly when stretched, one instead sees a significant reduction in specific volume under finite strains. This volume reduction is accompanied by increasing alignment of filaments along the strain axis and a nonlinear elastic response, with stiffening of the apparent Young's modulus. For semiflexible networks, in which entropic bending elasticity governs the linear elastic regime, the nonlinear Poisson effect is caused by the nonlinear force-extension relationship of the constituent filaments, which produces a highly asymmetric response of the constituent polymers to stretching and compression. The details of this relationship depend on the geometric and elastic properties of the underlying filaments, which can vary greatly in experimental systems. Here, we provide a comprehensive characterization of the nonlinear Poisson effect in an affine network model and explore the influence of filament properties on essential features of both microscopic and macroscopic response, including strain-driven alignment and volume reduction.

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.

Numerical analysis of the impact of cytoskeletal actin filament density alterations onto the diffusive vesicle-mediated cell transport

The interior of a eukaryotic cell is a highly complex composite material which consists of water, structural scaffoldings, organelles, and various biomolecular solutes. All these components serve as obstacles that impede the motion of vesicles. Hence, it is hypothesized that any alteration of the cytoskeletal network may directly impact or even disrupt the vesicle transport. A disruption of the vesicle-mediated cell transport is thought to contribute to several severe diseases and disorders, such as diabetes, Parkinson’s and Alzheimer’s disease, emphasizing the clinical relevance. To address the outlined objective, a multiscale finite element model of the diffusive vesicle transport is proposed on the basis of the concept of homogenization, owed to the complexity of the cytoskeletal network. In order to study the microscopic effects of specific nanoscopic actin filament network alterations onto the vesicle transport, a parametrized three-dimensional geometrical model of the actin filament network was generated on the basis of experimentally observed filament densities and network geometries in an adenocarcinomic human alveolar basal epithelial cell. Numerical analyzes of the obtained effective diffusion properties within two-dimensional sampling domains of the whole cell model revealed that the computed homogenized diffusion coefficients can be predicted statistically accurate by a simple two-parameter power law as soon as the inaccessible area fraction, due to the obstacle geometries and the finite size of the vesicles, is known. This relationship, in turn, leads to a massive reduction in computation time and allows to study the impact of a variety of different cytoskeletal alterations onto the vesicle transport. Hence, the numerical simulations predicted a 35% increase in transport time due to a uniformly distributed four-fold increase of the total filament amount. On the other hand, a hypothetically reduced expression of filament cross-linking proteins led to sparser filament networks and, thus, a speed up of the vesicle transport.

Many vital processes in our eukaryotic cells and organs require an astonishingly precise routing of intermediate products to various intra- and extracellular destinations using vesicles as transporters. This can be illustrated by numerous examples, such as the production and destruction of proteins, the export of neurotransmitters or insulin to the extracellular domain, etc. However, the inside of a cell is tightly packed with numerous structural scaffoldings (filaments), which serve as obstacles and impede the vesicle motion. It is thought that any disturbances of the vesicle-mediated cell transport contribute to numerous degenerative diseases and disorders, which highlights the clinical relevance for investigating this intracellular transport mechanism by developing computational models and performing experimental studies. In this study, we numerically quantified how different specific alterations of the filament density inside a human lung cell—due to changed mechanical loadings or genetic disorders of proteins being responsible for filament branching—affect the diffusion of vesicles inside the intracellular fluid. Therefore, based on the concept of homogenization, a computationally efficient numerical method was developed and utilized to simulate the diffusion of vesicles inside the whole cell, considering the detailed structural information of the filament network.

On the mechanical response of the actomyosin cortex during cell indentations

A mechanical model is presented to analyze the mechanics and dynamics of the cell cortex during indentation. We investigate the impact of active contraction on the cross-linked actin network for different probe sizes and indentation rates. The essential molecular mechanisms of filament stretching, cross-linking and motor activity, are represented by an active and viscous mechanical continuum. The filaments behave as worm-like chains linked either by passive rigid linkers or by myosin motors. In the first example, the effects of probe size and loading rate are evaluated using the model for an idealized rounded cell shape in which properties are based on the results of parallel-plate rheometry available in the literature. Extreme cases of probe size and indentation rate are taken into account. Afterward, AFM experiments were done by engaging smooth muscle cells with both sharp and spherical probes. By inverse analysis with finite element software, our simulations mimicking the experimental conditions show the model is capable of fitting the AFM data. The results provide spatiotemporal dependence on the size and rate of the mechanical stimuli. The model captures the general features of the cell response. It characterizes the actomyosin cortex as an active solid at short timescales and as a fluid at longer timescales by showing (1) higher levels of contraction in the zones of high curvature; (2) larger indentation forces as the probe size increases; and (3) increase in the apparent modulus with the indentation depth but no dependence on the rate of the mechanical stimuli. The methodology presented in this work can be used to address and predict microstructural dependence on the force generation of living cells, which can contribute to understanding the broad spectrum of results in cell experiments.

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.

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.

Numerical investigation of the influence of pattern topology on the mechanical behavior of PEGDA hydrogels

Poly(ethylene glycol) diacrylate (PEGDA) hydrogels can be potentially used as scaffold material for tissue engineered heart valves (TEHVs) due to their good biocompatibility and biomechanical tunability. The photolithographic patterning technique is an effective approach to pattern PEGDA hydrogels to mimic the mechanical behavior of native biological tissues that are intrinsically anisotropic. The material properties of patterned PEGDA hydrogels largely depend on the pattern topology. In this paper, we adopt a newly proposed computational framework for fibrous biomaterials to numerically investigate the influence of pattern topology, including pattern ratio, orientation and waviness, on the mechanical behavior of patterned PEGDA hydrogels. The material parameters for the base hydrogel and the pattern stripes are directly calibrated from published experimental data. Several experimental observations reported in the literature are captured in the simulation, including the nonlinear relationship between pattern ratio and material linear modulus, and the decrease of material anisotropy when pattern ratio increases. We further numerically demonstrate that a three-region (toe-heel-linear) stress-strain relationship typically exhibited by biological tissues can be obtained by tuning the pattern waviness and the relative stiffness between the base hydrogel and pattern stripes. The numerical strategy and simulation results presented here can provide helpful guidance to optimize pattern design of PEGDA hydrogels toward the targeted material mechanical properties, therefore advance the development of TEHVs.

STATEMENT OF SIGNIFICANCE

Poly(ethylene glycol) diacrylate (PEGDA) hydrogels can be used as scaffold material for tissue engineered heart values (TEHVs) providing a promising alternative to generate suitable heart valve replacement method. The patterning of PEGDA hydrogels using photolithographic techniques creates materials that mimic the mechanical behavior of native heart valve tissues. However, targeted material properties are obtained via a trial-and-error process. Depending on experiments alone to explore the influence of pattern topology is expensive and time-consuming. We combine a newly proposed computational framework with published experimental data to numerically investigate the influence of pattern geometry on the mechanical behavior of patterned PEGDA hydrogels. The numerical strategy and simulation results presented here can provide guidance to optimize the design of PEGDA hydrogels with targeted material properties, therefore advance the development of TEHVs.

Mechanical properties of gray and white matter brain tissue by indentation

The mammalian brain is composed of an outer layer of gray matter, consisting of cell bodies, dendrites, and unmyelinated axons, and an inner core of white matter, consisting primarily of myelinated axons. Recent evidence suggests that microstructural differences between gray and white matter play an important role during neurodevelopment. While brain tissue as a whole is rheologically well characterized, the individual features of gray and white matter remain poorly understood. Here we quantify the mechanical properties of gray and white matter using a robust, reliable, and repeatable method, flat-punch indentation. To systematically characterize gray and white matter moduli for varying indenter diameters, loading rates, holding times, post-mortem times, and locations we performed a series of n=192 indentation tests. We found that indenting thick, intact coronal slices eliminates the common challenges associated with small specimens: it naturally minimizes boundary effects, dehydration, swelling, and structural degradation. When kept intact and hydrated, brain slices maintained their mechanical characteristics with standard deviations as low as 5% throughout the entire testing period of five days post mortem. White matter, with an average modulus of 1.89 5kPa ± 0.592 kPa, was on average 39% stiffer than gray matter, p<0.01, with an average modulus of 1.389 kPa ± 0.289 kPa, and displayed larger regional variations. It was also more viscous than gray matter and responded less rapidly to mechanical loading. Understanding the rheological differences between gray and white matter may have direct implications on diagnosing and understanding the mechanical environment in neurodevelopment and neurological disorders.

From single fiber to macro-level mechanics: A structural finite-element model for elastomeric fibrous biomaterials

In the present work, we demonstrate that the mesoscopic in-plane mechanical behavior of membrane elastomeric scaffolds can be simulated by replication of actual quantified fibrous geometries. Elastomeric electrospun polyurethane (ES-PEUU) scaffolds, with and without particulate inclusions, were utilized. Simulations were developed from experimentally-derived fiber network geometries, based on a range of scaffold isotropic and anisotropic behaviors. These were chosen to evaluate the effects on macro-mechanics based on measurable geometric parameters such as fiber intersections, connectivity, orientation, and diameter. Simulations were conducted with only the fiber material model parameters adjusted to match the macro-level mechanical test data. Fiber model validation was performed at the microscopic level by individual fiber mechanical tests using AFM. Results demonstrated very good agreement to the experimental data, and revealed the formation of extended preferential fiber orientations spanning the entire model space. We speculate that these emergent structures may be responsible for the tissue-like macroscale behaviors observed in electrospun scaffolds. To conclude, the modeling approach has implications for (1) gaining insight on the intricate relationship between fabrication variables, structure, and mechanics to manufacture more functional devices/materials, (2) elucidating the effects of cell or particulate inclusions on global construct mechanics, and (3) fabricating better performing tissue surrogates that could recapitulate native tissue mechanics.

An affine continuum mechanical model for cross-linked F-actin networks with compliant linker proteins

Cross-linked actin networks are important building blocks of the cytoskeleton. In order to gain deeper insight into the interpretation of experimental data on actin networks, adequate models are required. In this paper we introduce an affine constitutive network model for cross-linked F-actin networks based on nonlinear continuum mechanics, and specialize it in order to reproduce the experimental behavior of in vitro reconstituted model networks. The model is based on the elastic properties of single filaments embedded in an isotropic matrix such that the overall properties of the composite are described by a free-energy function. In particular, we are able to obtain the experimentally determined shear and normal stress responses of cross-linked actin networks typically observed in rheometer tests. In the present study an extensive analysis is performed by applying the proposed model network to a simple shear deformation. The single filament model is then extended by incorporating the compliance of cross-linker proteins and further extended by including viscoelasticity. All that is needed for the finite element implementation is the constitutive model for the filaments, the linkers and the matrix, and the associated elasticity tensor in either the Lagrangian or Eulerian formulation. The model facilitates parameter studies of experimental setups such as micropipette aspiration experiments and we present such studies to illustrate the efficacy of this modeling approach.

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.