Geometric characterization and simulation of planar layered elastomeric fibrous biomaterials

Structural descriptor and surrogate modeling for design of biodegradable scaffolds

Biodegradable scaffolds are important to regenerative medicine in that they provide an amicable environment for tissue regrowth. However, establishing structure-property (SP) relationships for scaffold design is challenging due to the complexity of the three-dimensional porous scaffold geometry. The complexity requires high-dimensional geometric descriptors. The training of such a SP surrogate model will need a large amount of experimental or simulation data. In this work, a schema of constructing SP relationship surrogates is developed to predict the degraded mechanical properties from the initial scaffold geometry. A new structure descriptor, the extended surfacelet transform (EST), is proposed to capture important details of pores associated with the degradation of scaffolds. The efficiency is further enhanced with principal component analysis to reduce the high-dimensional EST data into a low-dimensional representation. The schema also includes a kinetic Monte Carlo biodegradation model to simulate the biodegradation of polymer scaffolds and to generate the training data for the formation of SP relationships. The schema is demonstrated with the design of polycaprolactone biodegradable scaffolds by connecting the initial scaffold geometry to the degraded compressive modulus.

Bioinspired 3D microprinted cell scaffolds: Integration of graph theory to recapitulate complex network wiring in lymph nodes

Physical networks are ubiquitous in nature, but many of them possess a complex organizational structure that is difficult to recapitulate in artificial systems. This is especially the case in biomedical and tissue engineering, where the microstructural details of 3D cell scaffolds are important. Studies of biological networks-such as fibroblastic reticular cell (FRC) networks-have revealed the crucial role of network topology in a range of biological functions. However, cell scaffolds are rarely analyzed, or designed, using graph theory. To understand how networks affect adhered cells, 3D culture platforms capturing the complex topological properties of biologically relevant networks would be needed. In this work, we took inspiration from the small-world organization (high clustering and low path length) of FRC networks to design cell scaffolds. An algorithmic toolset was created to generate the networks and process them to improve their 3D printability. We employed tools from graph theory to show that the networks were small-world (omega factor, ω = -0.10 ± 0.02; small-world propensity, SWP = 0.74 ± 0.01). 3D microprinting was employed to physicalize networks as scaffolds, which supported the survival of FRCs. This work, therefore, represents a bioinspired, graph theory-driven approach to control the networks of microscale cell niches.

Dynamical modeling reveals RNA decay mediates the effect of matrix stiffness on aged muscle stem cell fate

Loss of muscle stem cell (MuSC) self-renewal with aging reflects a combination of influences from the intracellular (e.g., post-transcriptional modifications) and extracellular (e.g., matrix stiffness) environment. Whereas conventional single cell analyses have revealed valuable insights into factors contributing to impaired self-renewal with age, most are limited by static measurements that fail to capture nonlinear dynamics. Using bioengineered matrices mimicking the stiffness of young and old muscle, we showed that while young MuSCs were unaffected by aged matrices, old MuSCs were phenotypically rejuvenated by young matrices. Dynamical modeling of RNA velocity vector fields in silico revealed that soft matrices promoted a self-renewing state in old MuSCs by attenuating RNA decay. Vector field perturbations demonstrated that the effects of matrix stiffness on MuSC self-renewal could be circumvented by fine-tuning the expression of the RNA decay machinery. These results demonstrate that post-transcriptional dynamics dictate the negative effect of aged matrices on MuSC self-renewal.

Simultaneous Wide-Field Planar Strain-Fiber Orientation Distribution Measurement Using Polarized Spatial Domain Imaging

In the present study, we demonstrate that soft tissue fiber architectural maps captured using polarized spatial frequency domain imaging (pSFDI) can be utilized as an effective texture source for DIC-based planar surface strain analyses. Experimental planar biaxial mechanical studies were conducted using pericardium as the exemplar tissue, with simultaneous pSFDI measurements taken. From these measurements, the collagen fiber preferred direction [Formula: see text] was determined at the pixel level over the entire strain range using established methods ( https://doi.org/10.1007/s10439-019-02233-0 ). We then utilized these pixel-level [Formula: see text] maps as a texture source to quantify the deformation gradient tensor [Formula: see text] as a function of spatial position [Formula: see text] within the specimen at time t. Results indicted that that the pSFDI approach produced accurate deformation maps, as validated using both physical markers and conventional particle based method derived from the DIC analysis of the same specimens. We then extended the pSFDI technique to extract the fiber orientation distribution [Formula: see text] as a function of [Formula: see text] from the pSFDI intensity signal. This was accomplished by developing a calibration procedure to account for the optical behavior of the constituent fibers for the soft tissue being studied. We then demonstrated that the extracted [Formula: see text] was accurately computed in both the referential (i.e. unloaded) and deformed states. Moreover, we noted that the measured [Formula: see text] agreed well with affine kinematic deformation predictions. We also demonstrated this calibration approach could also be effectively used on electrospun biomaterials, underscoring the general utility of the approach. In a final step, using the ability to simultaneously quantify [Formula: see text] and [Formula: see text], we examined the effect of deformation and collagen structural measurements on the measurement region size. For pericardial tissues, we determined a critical length of [Formula: see text] 8 mm wherein the regional variations sufficiently dissipated. This result has immediate potential in the identification of optimal length scales for meso-scale strain measurement in soft tissues and fibrous biomaterials.

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.

Anisotropic elastic behavior of a hydrogel-coated electrospun polyurethane: Suitability for heart valve leaflets

Although xenograft biomaterials have been used for decades in replacement heart valves, they continue to face multiple limitations, including limited durability, mineralization, and restricted design space due to their biological origins. These issues necessitate the need for novel replacement heart valve biomaterials that are durable, non-thrombogenic, and compatible with transcatheter aortic valve replacement devices. In this study, we explored the suitability of an electrospun poly(carbonate urethane) (ES-PCU) mesh coated with a poly(ethylene glycol) diacrylate (PEGDA) hydrogel as a synthetic biomaterial for replacement heart valve leaflets. In this material design, the mesh provides the mechanical support, while the hydrogel provides the required surface hemocompatibility. We conducted a comprehensive study to characterize the structural and mechanical properties of the uncoated mesh as well as the hydrogel-coated mesh (composite biomaterial) over the estimated operational range. We found that the composite biomaterial was functionally robust with reproducible stress-strain behavior within and beyond the functional ranges for replacement heart valves, and was able to withstand the rigors of mechanical evaluation without any observable damage. In addition, the composite biomaterial displayed a wide range of mechanical anisotropic responses, which were governed by fiber orientation of the mesh, which in turn, was controlled with the fabrication process. Finally, we developed a novel constitutive modeling approach to predict the mechanical behavior of the composite biomaterial under in-plane extension and shear deformation modes. This model identified the existence of fiber-fiber mechanical interactions in the mesh that have not previously been reported. Interestingly, there was no evidence of fiber-hydrogel mechanical interactions. This important finding suggests that the hydrogel coating can be optimized for hemocompatibility independent of the structural mechanical responses required by the leaflet. This initial study indicated that the composite biomaterial has mechanical properties well-suited for replacement heart valve applications and that the electrospun mesh microarchitecture and hydrogel biological properties can be optimized independently. It also reveals that the structural mechanisms contributing to the mechanical response are more complicated than what was previously established and paves the pathway for more detailed future studies.

Why Are Viscosity and Nonlinearity Bound to Make an Impact in Clinical Elastographic Diagnosis?

The adoption of multiscale approaches by the biomechanical community has caused a major improvement in quality in the mechanical characterization of soft tissues. The recent developments in elastography techniques are enabling in vivo and non-invasive quantification of tissues' mechanical properties. Elastic changes in a tissue are associated with a broad spectrum of pathologies, which stems from the tissue microstructure, histology and biochemistry. This knowledge is combined with research evidence to provide a powerful diagnostic range of highly prevalent pathologies, from birth and labor disorders (prematurity, induction failures, etc.), to solid tumors (e.g., prostate, cervix, breast, melanoma) and liver fibrosis, just to name a few. This review aims to elucidate the potential of viscous and nonlinear elastic parameters as conceivable diagnostic mechanical biomarkers. First, by providing an insight into the classic role of soft tissue microstructure in linear elasticity; secondly, by understanding how viscosity and nonlinearity could enhance the current diagnosis in elastography; and finally, by compounding preliminary investigations of those elastography parameters within different technologies. In conclusion, evidence of the diagnostic capability of elastic parameters beyond linear stiffness is gaining momentum as a result of the technological and imaging developments in the field of biomechanics.

Contribution of nascent cohesive fiber-fiber interactions to the non-linear elasticity of fibrin networks under tensile load

Fibrin is a viscoelastic proteinaceous polymer that determines the deformability and integrity of blood clots and fibrin-based biomaterials in response to biomechanical forces. Here, a previously unnoticed structural mechanism of fibrin clots' mechanical response to external tensile loads is tested using high-resolution confocal microscopy and recently developed three-dimensional computational model. This mechanism, underlying local strain-stiffening of individual fibers as well as global stiffening of the entire network, is based on previously neglected nascent cohesive pairwise interactions between individual fibers (crisscrossing) in fibrin networks formed under tensile load. Existence of fiber-fiber crisscrossings of reoriented fibers was confirmed using 3D imaging of experimentally obtained stretched fibrin clots. The computational model enabled us to study structural details and quantify mechanical effects of the fiber-fiber cohesive crisscrossing during stretching of fibrin gels at various spatial scales. The contribution of the fiber-fiber cohesive contacts to the elasticity of stretched fibrin networks was characterized by changes in individual fiber stiffness, the length, width, and alignment of fibers, as well as connectivity and density of the entire network. The results show that the nascent cohesive crisscrossing of fibers in stretched fibrin networks comprise an underappreciated important structural mechanism underlying the mechanical response of fibrin to (patho)physiological stresses that determine the course and outcomes of thrombotic and hemostatic disorders, such as heart attack and ischemic stroke. STATEMENT OF SIGNIFICANCE: Fibrin is a viscoelastic proteinaceous polymer that determines the deformability and integrity of blood clots and fibrin-based biomaterials in response to biomechanical forces. In this paper, a novel structural mechanism of fibrin clots' mechanical response to external tensile loads is tested using high-resolution confocal microscopy and newly developed computational model. This mechanism, underlying local strain-stiffening of individual fibers as well as global stiffening of the entire network, is based on previously neglected nascent cohesive pairwise interactions between individual fibers (crisscrossing) in fibrin networks formed under tensile load. Cohesive crisscrossing is an important structural mechanism that influences the mechanical response of blood clots and which can determine the outcomes of blood coagulation disorders, such as heart attacks and strokes.

Mechanical behavior of nonwoven non-crosslinked fibrous mats with adhesion and friction

We present a study of the mechanical behavior of planar fibrous mats stabilized by inter-fiber adhesion. Fibers of various degrees of tortuosity and of infinite and finite length are considered in separate models. Fibers are randomly distributed, are not cross-linked, and interact through adhesion and friction. The variation of structural parameters such as the mat thickness and the mean segment length between contacts along given fibers with the strength of adhesion is determined. These systems are largely dissipative in that most of the work performed during deformation is dissipated frictionally and only a small fraction is stored as strain energy. The response of the mats to tensile loading has three regimes: a short elastic regime in which no sliding at contacts is observed, a well-defined sliding regime characterized by strain hardening, and a rapid stiffening regime at larger strains. The third regime is due to the formation of stress paths after the fiber tortuosity is pulled out and is absent in mats of finite length fibers. Networks of finite length fibers lose stability during the second regime of deformation. The scaling of the yield stress, which characterizes the transition between the first and the second regimes, and of the second regime's strain hardening modulus, with system parameters such as the strength of adhesion and friction and the degree of fiber tortuosity are determined. The strength of mats of finite length fibers is also determined as a function of network parameters. These results are expected to become useful in the design of electrospun mats and other planar fibrous non-cross-linked networks.

Meso-scale topological cues influence extracellular matrix production in a large deformation, elastomeric scaffold model

Physical cues are decisive factors in extracellular matrix (ECM) formation and elaboration. Their transduction across scale lengths is an inherently symbiotic phenomenon that while influencing ECM fate is also mediated by the ECM structure itself. This study investigates the possibility of enhancing ECM elaboration by topological cues that, while not modifying the substrate macro scale mechanics, can affect the meso-scale strain range acting on cells incorporated within the scaffold. Vascular smooth muscle cell micro-integrated, electrospun scaffolds were fabricated with comparable macroscopic biaxial mechanical response, but different meso-scale topology. Seeded scaffolds were conditioned on a stretch bioreactor and exposed to large strain deformations. Samples were processed to evaluate ECM quantity and quality via: biochemical assay, qualitative and quantitative histological assessment and multi-photon analysis. Experimental evaluation was coupled to a numerical model that elucidated the relationship between the scaffold micro-architecture and the strain acting on the cells. Results showed an higher amount of ECM formation for the scaffold type characterized by lowest fiber intersection density. The numerical model simulations associated this result with the differences found for the change in cell nuclear aspect ratio and showed that given comparable macro scale mechanics, a difference in material topology created significant differences in cell-scaffold meso-scale deformations. These findings reaffirmed the role of cell shape in ECM formation and introduced a novel notion for the engineering of cardiac tissue where biomaterial structure can be designed to both mimick the organ level mechanics of a specific tissue of interest and elicit a desirable cellular response.

Structural modeling reveals microstructure-strength relationship for human ascending thoracic aorta

High lethality of aortic dissection necessitates accurate predictive metrics for dissection risk assessment. The not infrequent incidence of dissection at aortic diameters <5.5 cm, the current threshold guideline for surgical intervention (Nishimura et al., 2014), indicates an unmet need for improved evidence-based risk stratification metrics. Meeting this need requires a fundamental understanding of the structural mechanisms responsible for dissection evolution within the vessel wall. We present a structural model of the repeating lamellar structure of the aortic media comprised of elastic lamellae and collagen fiber networks, the primary load-bearing components of the vessel wall. This model was used to assess the role of these structural features in determining in-plane tissue strength, which governs dissection initiation from an intimal tear. Ascending aortic tissue specimens from three clinically-relevant patient populations were considered: non-aneurysmal aorta from patients with morphologically normal tricuspid aortic valve (CTRL), aneurysmal aorta from patients with tricuspid aortic valve (TAV), and aneurysmal aorta from patients with bicuspid aortic valve (BAV). Multiphoton imaging derived collagen fiber organization for each patient cohort was explicitly incorporated in our model. Model parameters were calibrated using experimentally-measured uniaxial tensile strength data in the circumferential direction for each cohort, while the model was validated by contrasting simulated tissue strength against experimentally-measured strength in the longitudinal direction. Orientation distribution, controlling the fraction of loaded collagen fibers at a given stretch, was identified as a key feature governing anisotropic tissue strength for all patient cohorts.

A 2.5D approach to the mechanics of electrospun fibre mats

In this paper, a discrete random network modelling approach specific to electrospun networks is presented. Owing to the manufacturing process, fibres in these materials systems have an enormous length, as compared to their diameters, and form sparse networks since fibre contact over thickness is limited to a narrow range. Representative volume elements are generated, in which fibres span the entire domain, and a technique is developed to apply computationally favourable periodic boundary conditions despite the initial non-periodicity of the networks. To capture sparsity, a physically motivated method is proposed to distinguish true fibre cross-links, in which mechanical interaction takes place, from mere fibre intersections. The model is exclusively informed by experimentally accessible parameters, demonstrates excellent agreement with the mechanical response of electrospun fibre mats, captures typical microscopic deformation patterns, and provides information on the kinematics of fibres and pores. This ability to address relevant mechanisms of deformation at both micro- and macroscopic length scales, together with the moderate computational cost, render the proposed modelling approach a highly qualified tool for the computer-based design and optimization of electrospun networks.

Large strain stimulation promotes extracellular matrix production and stiffness in an elastomeric scaffold model

Mechanical conditioning of engineered tissue constructs is widely recognized as one of the most relevant methods to enhance tissue accretion and microstructure, leading to improved mechanical behaviors. The understanding of the underlying mechanisms remains rather limited, restricting the development of in silico models of these phenomena, and the translation of engineered tissues into clinical application. In the present study, we examined the role of large strip-biaxial strains (up to 50%) on ECM synthesis by vascular smooth muscle cells (VSMCs) micro-integrated into electrospun polyester urethane urea (PEUU) constructs over the course of 3 weeks. Experimental results indicated that VSMC biosynthetic behavior was quite sensitive to tissue strain maximum level, and that collagen was the primary ECM component synthesized. Moreover, we found that while a 30% peak strain level achieved maximum ECM synthesis rate, further increases in strain level lead to a reduction in ECM biosynthesis. Subsequent mechanical analysis of the formed collagen fiber network was performed by removing the scaffold mechanical responses using a strain-energy based approach, showing that the denovo collagen also demonstrated mechanical behaviors substantially better than previously obtained with small strain training and comparable to mature collagenous tissues. We conclude that the application of large deformations can play a critical role not only in the quantity of ECM synthesis (i.e. the rate of mass production), but also on the modulation of the stiffness of the newly formed ECM constituents. The improved understanding of the process of growth and development of ECM in these mechano-sensitive cell-scaffold systems will lead to more rational design and manufacturing of engineered tissues operating under highly demanding mechanical environments.

A structural finite element model for lamellar unit of aortic media indicates heterogeneous stress field after collagen recruitment

Incorporation of collagen structural information into the study of biomechanical behavior of ascending thoracic aortic (ATA) wall tissue should provide better insight into the pathophysiology of ATA. Structurally motivated constitutive models that include fiber dispersion and recruitment can successfully capture overall mechanical response of the arterial wall tissue. However, these models cannot examine local microarchitectural features of the collagen network, such as the effect of fiber disruptions and interaction between fibrous and non-fibrous components, which may influence emergent biomechanical properties of the tissue. Motivated by this need, we developed a finite element based three-dimensional structural model of the lamellar units of the ATA media that directly incorporates the collagen fiber microarchitecture. The fiber architecture was computer generated utilizing network features, namely fiber orientation distribution, intersection density and areal concentration, obtained from image analysis of multiphoton microscopy images taken from human aneurysmal ascending thoracic aortic media specimens with bicuspid aortic valve (BAV) phenotype. Our model reproduces the typical J-shaped constitutive response of the aortic wall tissue. We found that the stress state in the non-fibrous matrix was homogeneous until the collagen fibers were recruited, but became highly heterogeneous after that event. The degree of heterogeneity was dependent upon local network architecture with high stresses observed near disrupted fibers. The magnitude of non-fibrous matrix stress at higher stretch levels was negatively correlated with local fiber density. The localized stress concentrations, elucidated by this model, may be a factor in the degenerative changes in aneurysmal ATA tissue.

On the presence of affine fibril and fiber kinematics in the mitral valve anterior leaflet

In this study, we evaluated the hypothesis that the constituent fibers follow an affine deformation kinematic model for planar collagenous tissues. Results from two experimental datasets were utilized, taken at two scales (nanometer and micrometer), using mitral valve anterior leaflet (MVAL) tissues as the representative tissue. We simulated MVAL collagen fiber network as an ensemble of undulated fibers under a generalized two-dimensional deformation state, by representing the collagen fibrils based on a planar sinusoidally shaped geometric model. The proposed approach accounted for collagen fibril amplitude, crimp period, and rotation with applied macroscopic tissue-level deformation. When compared to the small angle x-ray scattering measurements, the model fit the data well, with an r(2) = 0.976. This important finding suggests that, at the homogenized tissue-level scale of ∼1 mm, the collagen fiber network in the MVAL deforms according to an affine kinematics model. Moreover, with respect to understanding its function, affine kinematics suggests that the constituent fibers are largely noninteracting and deform in accordance with the bulk tissue. It also suggests that the collagen fibrils are tightly bounded and deform as a single fiber-level unit. This greatly simplifies the modeling efforts at the tissue and organ levels, because affine kinematics allows a straightforward connection between the macroscopic and local fiber strains. It also suggests that the collagen and elastin fiber networks act independently of each other, with the collagen and elastin forming long fiber networks that allow for free rotations. Such freedom of rotation can greatly facilitate the observed high degree of mechanical anisotropy in the MVAL and other heart valves, which is essential to heart valve function. These apparently novel findings support modeling efforts directed toward improving our fundamental understanding of tissue biomechanics in healthy and diseased conditions.