Emergent structure-dependent relaxation spectra in viscoelastic fiber networks in extension

Does the Extracellular Matrix Support Cell-Cell Communication by Elastic Wave Packets?

The extracellular matrix (ECM) is a fibrous network supporting biological cells and provides them a medium for interaction. Cells modify the ECM by applying traction forces, and these forces can propagate to long ranges and establish a mechanism of mechanical communication between neighboring cells. Previous studies have mainly focused on analysis of static force transmission across the ECM. In this study, we explore the plausibility of dynamic mechanical interaction, expressed as vibrations or abrupt fluctuations, giving rise to elastic waves propagating along ECM fibers. We use a numerical mass-spring model to simulate the longitudinal and transversal waves propagating along a single ECM fiber and across a 2D random fiber network. The elastic waves are induced by an active contracting cell (signaler) and received by a passive neighboring cell (receiver). We show that dynamic wave propagation may amplify the signal at the receiver end and support up to an order of magnitude stronger mechanical cues and longer-ranged communication relative to static transmission. Also, we report an optimal impulse duration corresponding to the most effective transmission, as well as extreme fast impulses, in which the waves are encaged around the active cell and do not reach the neighboring cell, possibly due to the Anderson localization effect. Finally, we also demonstrate that extracellular fluid viscosity reduces, but still allows, dynamic propagation along embedded ECM fibers. Our results motivate future biological experiments in mechanobiology to investigate, on the one hand, the mechanosensitivity of cells to dynamic forces traveling and guided by the ECM and, on the other hand, the impact of ECM architecture and remodeling on dynamic force transmission and its spectral filtering, dispersion, and decay.

REVIEW ARTICLE Engineering bio-inks for 3D bioprinting cell mechanical microenvironment

144Three-dimensional (3D) bioprinting has become a promising approach to constructing functional biomimetic tissues for tissue engineering and regenerative medicine. In 3D bioprinting, bio-inks are essential for the construction of cell microenvironment, thus affecting the biomimetic design and regenerative efficiency. Mechanical properties are one of the essential aspects of microenvironment, which can be characterized by matrix stiffness, viscoelasticity, topography, and dynamic mechanical stimulation. With the recent advances in functional biomaterials, various engineered bio-inks have realized the possibility of engineering cell mechanical microenvironment in vivo. In this review, we summarize the critical mechanical cues of cell microenvironments, review the engineered bio-inks while focusing on the selection principles for constructing cell mechanical microenvironments, and discuss the challenges facing this field and the possible solutions for them.

Elucidating the signal for contact guidance contained in aligned fibrils with a microstructural-mechanical model

Despite its importance in physiological processes and tissue engineering, the mechanism underlying cell contact guidance in an aligned fibrillar network has defied elucidation due to multiple interdependent signals that such a network presents to cells, namely, anisotropy of adhesion, porosity and mechanical behaviour. A microstructural-mechanical model of fibril networks was used to assess the relative magnitudes of these competing signals in networks of varied alignment strength based on idealized cylindrical pseudopods projected into the aligned and orthogonal directions and computing the anisotropy of metrics chosen for adhesion, porosity and mechanical behaviour: cylinder-fibre contact area for adhesion, persistence length of pores for porosity and total force to displace fibres from the cylindrical volume as well as network stiffness experienced upon cylinder retraction for mechanical behaviour. The signals related to mechanical anisotropy are substantially higher than adhesion and porosity anisotropy, especially at stronger network alignments, although their signal to noise (S/N) values are substantially lower. The former finding is consistent with a recent report that fibroblasts can sense fibril alignment via anisotropy of network mechanical resistance, and the model reveals this can be due to either mechanical resistance to pseudopod protrusion or retraction given their signal and S/N values are similar.

Stress relaxation in network materials: the contribution of the network

Stress relaxation in network materials with permanent crosslinks is due to the transport of fluid within the network (poroelasticity), the viscoelasticity of the matrix and the viscoelasticity of the network. While relaxation associated with the matrix was studied extensively, the contribution of the network remains unexplored. In this work we consider two and three-dimensional stochastic fiber networks with viscoelastic fibers and explore the dependence of stress relaxation on network structure. We observe that relaxation has two regimes - an initial exponential regime, followed by a stretched exponential regime - similar to the situation in other disordered materials. The stretch exponent is a function of density, fiber diameter and the network structure, and has a minimum at the transition between the affine and non-affine regimes of network behavior. The relaxation time constant of the first, exponential regime is similar to the relaxation time constant of individual fibers and is independent of network density and fiber diameter. The relaxation time constant of the second, stretched exponential regime is a weak function of network parameters. The stretched exponential emerges from the heterogeneity of relaxation dynamics on scales comparable with the mesh size, with higher heterogeneity leading to smaller stretch exponents. In composite networks of fibers whose relaxation time constant is selected from a distribution with set mean, the stretch exponent decreases with increasing the coefficient of variation of the fiber time constant distribution. As opposed to thermal glass formers and colloids, in these athermal systems the dynamic heterogeneity is introduced by the network structure and does not evolve during relaxation. While in thermal systems the control parameter is the temperature, in this athermal case the control parameter is a non-dimensional structural parameter which describes the degree of non-affinity of the network.

Sensitivity of the shear wave speed-stress relationship to soft tissue material properties and fiber alignment

The use of shear wave propagation to noninvasively measure material properties and loading in tendons and ligaments is a growing area of interest in biomechanics. Prior models and experiments suggest that shear wave speed primarily depends on the apparent shear modulus (i.e., shear modulus accounting for contributions from all constituents) at low loads, and then increases with axial stress when axially loaded. However, differences in the magnitudes of shear wave speeds between ligaments and tendons, which have different substructures, suggest that the tissue's composition and fiber alignment may also affect shear wave propagation. Accordingly, the objectives of this study were to (1) characterize changes in the apparent shear modulus induced by variations in constitutive properties and fiber alignment, and (2) determine the sensitivity of the shear wave speed-stress relationship to variations in constitutive properties and fiber alignment. To enable systematic variations of both constitutive properties and fiber alignment, we developed a finite element model that represented an isotropic ground matrix with an embedded fiber distribution. Using this model, we performed dynamic simulations of shear wave propagation at axial strains from 0% to 10%. We characterized the shear wave speed-stress relationship using a simple linear regression between shear wave speed squared and axial stress, which is based on an analytical relationship derived from a tensioned beam model. We found that predicted shear wave speeds were both in-range with shear wave speeds in previous in vivo and ex vivo studies, and strongly correlated with the axial stress (R² = 0.99). The slope of the squared shear wave speed-axial stress relationship was highly sensitive to changes in tissue density. Both the intercept of this relationship and the apparent shear modulus were sensitive to both the shear modulus of the ground matrix and the stiffness of the fibers' toe-region when the fibers were less well-aligned to the loading direction. We also determined that the tensioned beam model overpredicted the axial tissue stress with increasing load when the model had less well-aligned fibers. This indicates that the shear wave speed increases likely in response to a load-dependent increase in the apparent shear modulus. Our findings suggest that researchers may need to consider both the material and structural properties (i.e., fiber alignment) of tendon and ligament when measuring shear wave speeds in pathological tissues or tissues with less well-aligned fibers.

A Hybrid Microstructural-Continuum Multiscale Approach for Modeling Hyperelastic Fibrous Soft Tissue

The heterogeneous, nonlinear, anisotropic material behavior of biological tissues makes precise definition of an accurate constitutive model difficult. One possible solution to this issue would be to define microstructural elements and perform fully coupled multiscale simulation. However, for complex geometries and loading scenarios, the computational costs of such simulations can be prohibitive. Ideally then, we should seek a method that contains microstructural detail, but leverages the speed of classical continuum-based finite-element (FE) modeling. In this work, we demonstrate the use of the Holzapfel-Gasser-Ogden (HGO) model [1, 2] to fit the behavior of microstructural network models. We show that Delaunay microstructural networks can be fit to the HGO strain energy function by calculating fiber network strain energy and average fiber stretch ratio. We then use the HGO constitutive model in a FE framework to improve the speed of our hybrid model, and demonstrate that this method, combined with a material property update scheme, can match a full multiscale simulation. This method gives us flexibility in defining complex FE simulations that would be impossible, or at least prohibitively time consuming, in multiscale simulation, while still accounting for microstructural heterogeneity.

Viscoelastic characterization of human descending thoracic aortas under cyclic load

Experiments were carried out on 15 human descending thoracic aortas from heart-beating healthy donors who donated organs for transplant. The aortas were kept refrigerated in organ preservation solution and tested were completed within 48 hours from explant. Donors' age was comprised between 25 and 70 years, with an average of 51.7 ± 12.8 years. Quasi-static and dynamic uniaxial tensile test were carried out in thermally controlled physiological saline solution in order to characterize the viscoelastic behavior. Strips were tested under harmonic deformation of different frequency, between 1 and 11 Hz, at three initial pre-stretches. Cyclic deformations of two different amplitudes were used: a physiological one and a small one, the latter one for comparison purposes to understand the accuracy limits of viscoelastic models. Aortic strips in circumferential and longitudinal directions were cut from each aorta. Some strips were dissected to separate the three layers: intima, media and adventitia. They were tested individually in order to obtain layer-specific data. However, strips of the intact wall were also tested. Therefore, 8 strips per donors were tested. Viscoelastic parameters are accurately evaluated from the hysteresis loops. Results show that small-amplitude cyclic strain over-estimate the storage modulus and under-estimate the loss-factor. Therefore, cyclic deformation of physiological amplitude is necessary to obtain correct viscoelastic data of aortic tissue. The value of the applied pre-stretch is significant on the dynamic stiffness ratio (storage modulus divided by the corresponding quasi-static stiffness), while it is less significant for the loss factor. The median of the dynamic stiffness ratios, in physiological conditions, varies between 1.14 and 1.33 for the different layers and the intact wall; the corresponding median of the loss factors varies between 0.050 and 0.066. The lowest dynamic stiffness ratios and loss factors were obtained from donors of the youngest age group. STATEMENT OF SIGNIFICANCE: There is an increasing interest in replacing traditional Dacron grafts used to repair thoracic aortas after acute dissection and aneurysm, with grafts in innovative biomaterials that mimic the mechanical properties and the dynamic behavior of the aorta. The human aorta is a complex laminated structure with hyperelastic and viscoelastic material properties and residual stresses. This study aims to characterize the nonlinear viscoelastic properties of ex-vivo human descending thoracic aortas by measuring hysteresis loops of physiological amplitude under harmonic strain. Results show the necessity to characterize the viscoelastic material properties of the aorta under physiological conditions, as well as the necessity to introduce improved models that take better into account the influence of the initial pre-stretch and amplitude of the cyclic load.

Evaluation of affine fiber kinematics in porcine tricuspid valve leaflets using polarized spatial frequency domain imaging and planar biaxial testing

Collagen fibers are the primary load-bearing microstructural constituent of bodily soft tissues, and, when subjected to external loading, the collagen fibers reorient, uncrimp, and elongate. Specific to the atrioventricular heart valve leaflets, the collagen fiber kinematics form the basis of many constitutive models; however, some researchers claim that modeling the affine fiber kinematics (AFK) are sufficient for accurately predicting the macroscopic tissue deformations, while others state that modeling the non-affine kinematics (i.e., fiber uncrimping together with elastic elongation) is required. Experimental verification of the AFK theory has been previously performed for the mitral valve leaflets in the left-side heart; however, this same evaluation has yet to be performed for the morphologically distinct tricuspid valve (TV) leaflets in the right-side heart. In this work, we, for the first time, evaluated the AFK theory for the TV leaflets using an integrated biaxial testing-polarized spatial frequency domain imaging device to experimentally quantify the load-dependent collagen fiber reorientations for comparison to the AFK theory predictions. We found that the AFK theory generally underpredicted the fiber reorientations by 3.1°, on average, under the applied equibiaxial loading with greater disparity when the tissue was subjected to the applied non-equibiaxial loading. Furthermore, increased AFK errors were observed with increasing collagen fiber reorientations (Pearson coefficient r = -0.36, equibiaxial loading), suggesting the AFK theory is better suited for relatively smaller reorientations. Our findings suggest the AFK theory may require modification for more accurate predictions of the collagen fiber kinematics in the TV leaflets, which will be useful in refining modeling efforts for more accurate TV simulations.