1
|
Panchenko AY, Tchaicheeyan O, Berinskii IE, Lesman A. Does the Extracellular Matrix Support Cell-Cell Communication by Elastic Wave Packets? ACS Biomater Sci Eng 2022; 8:5155-5170. [PMID: 36346743 DOI: 10.1021/acsbiomaterials.2c01049] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
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.
Collapse
Affiliation(s)
- Artem Y Panchenko
- School of Mechanical Engineering, Faculty of Engineering, Tel Aviv University, Tel Aviv69978, Israel
| | - Oren Tchaicheeyan
- School of Mechanical Engineering, Faculty of Engineering, Tel Aviv University, Tel Aviv69978, Israel
| | - Igor E Berinskii
- School of Mechanical Engineering, Faculty of Engineering, Tel Aviv University, Tel Aviv69978, Israel
| | - Ayelet Lesman
- School of Mechanical Engineering, Faculty of Engineering, Tel Aviv University, Tel Aviv69978, Israel.,The Center for the Physics and Chemistry of Living Systems, Tel Aviv University, Tel Aviv69978, Israel
| |
Collapse
|
2
|
Yang Y, Jia Y, Yang Q, Xu F. REVIEW ARTICLE Engineering bio-inks for 3D bioprinting cell mechanical microenvironment. Int J Bioprint 2022; 9:632. [PMID: 36844247 PMCID: PMC9947384 DOI: 10.18063/ijb.v9i1.632] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2022] [Accepted: 08/29/2022] [Indexed: 11/05/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Yanshen Yang
- The Key Laboratory of Biomedical Information Engineering of Ministry of Education, School of Life Science and Technology, Xi’an Jiaotong University, Xi’an, Shaanxi, 710049, P.R. China,Bioinspired Engineering and Biomechanics Center (BEBC), Xi’an Jiaotong University, Xi’an, Shaanxi, 710049, P.R. China
| | - Yuanbo Jia
- The Key Laboratory of Biomedical Information Engineering of Ministry of Education, School of Life Science and Technology, Xi’an Jiaotong University, Xi’an, Shaanxi, 710049, P.R. China,Bioinspired Engineering and Biomechanics Center (BEBC), Xi’an Jiaotong University, Xi’an, Shaanxi, 710049, P.R. China
| | - Qingzhen Yang
- The Key Laboratory of Biomedical Information Engineering of Ministry of Education, School of Life Science and Technology, Xi’an Jiaotong University, Xi’an, Shaanxi, 710049, P.R. China,Bioinspired Engineering and Biomechanics Center (BEBC), Xi’an Jiaotong University, Xi’an, Shaanxi, 710049, P.R. China
| | - Feng Xu
- The Key Laboratory of Biomedical Information Engineering of Ministry of Education, School of Life Science and Technology, Xi’an Jiaotong University, Xi’an, Shaanxi, 710049, P.R. China,Bioinspired Engineering and Biomechanics Center (BEBC), Xi’an Jiaotong University, Xi’an, Shaanxi, 710049, P.R. China,Corresponding author: Feng Xu ()
| |
Collapse
|
3
|
Isogai T, Hayashi M. Critical Effects of Branch Numbers at the Cross-Link Point on the Relaxation Behaviors of Transesterification Vitrimers. Macromolecules 2022. [DOI: 10.1021/acs.macromol.2c00560] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Taketo Isogai
- Department of Life Science and Applied Chemistry, Graduate School of Engineering, Nagoya Institute of Technology, Gokiso-cho, Showa-ku, Nagoya 466-8555, Japan
| | - Mikihiro Hayashi
- Department of Life Science and Applied Chemistry, Graduate School of Engineering, Nagoya Institute of Technology, Gokiso-cho, Showa-ku, Nagoya 466-8555, Japan
| |
Collapse
|
4
|
Bersie-Larson LM, Lai VK, Dhume RY, Provenzano PP, Barocas VH, Tranquillo RT. Elucidating the signal for contact guidance contained in aligned fibrils with a microstructural-mechanical model. J R Soc Interface 2022; 19:20210951. [PMID: 35582810 DOI: 10.1098/rsif.2021.0951] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Lauren M Bersie-Larson
- Department of Biomedical Engineering, University of Minnesota - Twin Cities, Minneapolis, MN 55455, USA
| | - Victor K Lai
- Department of Chemical Engineering, University of Minnesota - Duluth, Duluth, MN, USA
| | - Rohit Y Dhume
- Department of Mechanical Engineering, University of Minnesota - Twin Cities, Minneapolis, MN 55455, USA
| | - Paolo P Provenzano
- Department of Biomedical Engineering, University of Minnesota - Twin Cities, Minneapolis, MN 55455, USA
| | - Victor H Barocas
- Department of Biomedical Engineering, University of Minnesota - Twin Cities, Minneapolis, MN 55455, USA
| | - Robert T Tranquillo
- Department of Biomedical Engineering, University of Minnesota - Twin Cities, Minneapolis, MN 55455, USA.,Department of Chemical Engineering and Materials Science, University of Minnesota - Twin Cities, Minneapolis, MN 55455, USA
| |
Collapse
|
5
|
Amjad SN, Picu RC. Stress relaxation in network materials: the contribution of the network. SOFT MATTER 2022; 18:446-454. [PMID: 34913052 DOI: 10.1039/d1sm01546j] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
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.
Collapse
Affiliation(s)
- S N Amjad
- Department of Mechanical, Aerospace and Nuclear Engineering, Rensselaer Polytechnic Institute, Troy, NY 12180, USA.
| | - R C Picu
- Department of Mechanical, Aerospace and Nuclear Engineering, Rensselaer Polytechnic Institute, Troy, NY 12180, USA.
| |
Collapse
|
6
|
Blank JL, Thelen DG, Allen MS, Roth JD. Sensitivity of the shear wave speed-stress relationship to soft tissue material properties and fiber alignment. J Mech Behav Biomed Mater 2021; 125:104964. [PMID: 34800889 PMCID: PMC8666097 DOI: 10.1016/j.jmbbm.2021.104964] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2021] [Revised: 10/26/2021] [Accepted: 11/06/2021] [Indexed: 01/03/2023]
Abstract
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 (R2 = 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.
Collapse
Affiliation(s)
- Jonathon L. Blank
- Department of Mechanical Engineering, University of Wisconsin-Madison, Madison, WI, USA
| | - Darryl G. Thelen
- Department of Mechanical Engineering, University of Wisconsin-Madison, Madison, WI, USA,Department of Biomedical Engineering, University of Wisconsin-Madison, Madison, WI, USA
| | - Matthew S. Allen
- Department of Mechanical Engineering, Brigham Young University, Provo, UT, USA
| | - Joshua D. Roth
- Department of Mechanical Engineering, University of Wisconsin-Madison, Madison, WI, USA,Department of Orthopedics and Rehabilitation, University of Wisconsin-Madison, Madison, WI, USA
| |
Collapse
|
7
|
Nikpasand M, Mahutga RR, Bersie-Larson LM, Gacek E, Barocas VH. A Hybrid Microstructural-Continuum Multiscale Approach for Modeling Hyperelastic Fibrous Soft Tissue. JOURNAL OF ELASTICITY 2021; 145:295-319. [PMID: 36380845 PMCID: PMC9648697 DOI: 10.1007/s10659-021-09843-7] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/30/2020] [Accepted: 05/19/2021] [Indexed: 06/16/2023]
Abstract
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.
Collapse
Affiliation(s)
- Maryam Nikpasand
- Department of Mechanical Engineering, University of Minnesota – Twin Cities, Minneapolis, MN, USA
| | - Ryan R. Mahutga
- Department of Biomedical Engineering, University of Minnesota – Twin Cities, Minneapolis, MN, USA
| | - Lauren M. Bersie-Larson
- Department of Biomedical Engineering, University of Minnesota – Twin Cities, Minneapolis, MN, USA
| | - Elizabeth Gacek
- Department of Biomedical Engineering, University of Minnesota – Twin Cities, Minneapolis, MN, USA
| | - Victor H. Barocas
- Department of Biomedical Engineering, University of Minnesota – Twin Cities, Minneapolis, MN, USA
| |
Collapse
|
8
|
Franchini G, Breslavsky ID, Holzapfel GA, Amabili M. Viscoelastic characterization of human descending thoracic aortas under cyclic load. Acta Biomater 2021; 130:291-307. [PMID: 34082105 DOI: 10.1016/j.actbio.2021.05.025] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2021] [Revised: 05/21/2021] [Accepted: 05/26/2021] [Indexed: 12/20/2022]
Abstract
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.
Collapse
Affiliation(s)
- Giulio Franchini
- Department of Mechanical Engineering, McGill University, Montreal, Canada
| | - Ivan D Breslavsky
- Department of Mechanical Engineering, McGill University, Montreal, Canada
| | - Gerhard A Holzapfel
- Institute of Biomechanics, Graz University of Technology, Austria; Department of Structural Engineering, Norwegian University of Science and Technology, Trondheim, Norway
| | - Marco Amabili
- Department of Mechanical Engineering, McGill University, Montreal, Canada; Dipartimento di Ingegneria e Architettura, University of Parma, Parma, Italy.
| |
Collapse
|
9
|
Evaluation of affine fiber kinematics in porcine tricuspid valve leaflets using polarized spatial frequency domain imaging and planar biaxial testing. J Biomech 2021; 123:110475. [PMID: 34004393 DOI: 10.1016/j.jbiomech.2021.110475] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2021] [Revised: 04/16/2021] [Accepted: 04/19/2021] [Indexed: 11/22/2022]
Abstract
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.
Collapse
|
10
|
Amabili M, Balasubramanian P, Breslavsky I. Anisotropic fractional viscoelastic constitutive models for human descending thoracic aortas. J Mech Behav Biomed Mater 2019; 99:186-197. [DOI: 10.1016/j.jmbbm.2019.07.010] [Citation(s) in RCA: 39] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2019] [Revised: 06/11/2019] [Accepted: 07/18/2019] [Indexed: 11/26/2022]
|