Effects of Collagen Heterogeneity on Myocardial Infarct Mechanics in a Multiscale Fiber Network Model

Mechanical Models of Collagen Networks for Understanding Changes in the Failure Properties of Aging Skin

Skin undergoes mechanical alterations due to changes in the composition and structure of the collagenous dermis with aging. Previous studies have conflicting findings, with both increased and decreased stiffness reported for aging skin. The underlying structure-function relationships that drive age-related changes are complex and difficult to study individually. One potential contributor to these variations is the accumulation of nonenzymatic crosslinks within collagen fibers, which affect dermal collagen remodeling and mechanical properties. Specifically, these crosslinks make individual fibers stiffer in their plastic loading region and lead to increased fragmentation of the collagenous network. To better understand the influence of these changes, we investigated the impact of nonenzymatic crosslink changes on the dermal microstructure using discrete fiber networks representative of the dermal microstructure. Our findings suggest that stiffening the plastic region of collagen's mechanical response has minimal effects on network-level stiffness and failure stresses. Conversely, simulating fragmentation through a loss of connectivity substantially reduces network stiffness and failure stress, while increasing stretch ratios at failure.

An image-driven micromechanical approach to characterize multiscale remodeling in infarcted myocardium

Myocardial infarction (MI) is accompanied by the formation of a fibrotic scar in the left ventricle (LV) and initiates significant alterations in the architecture and constituents of the LV free wall (LVFW). Previous studies have shown that LV adaptation is highly individual, indicating that the identification of remodeling mechanisms post-MI demands a fully subject-specific approach that can integrate a host of structural alterations at the fiber-level to changes in bulk biomechanical adaptation at the tissue-level. We present an image-driven micromechanical approach to characterize remodeling, assimilating new biaxial mechanical data, histological studies, and digital image correlation data within an in-silico framework to elucidate the fiber-level remodeling mechanisms that drive tissue-level adaptation for each subject. We found that a progressively diffused collagen fiber structure combined with similarly disorganized myofiber architecture in the healthy region leads to the loss of LVFW anisotropy post-MI, offering an important tissue-level hallmark for LV maladaptation. In contrast, our results suggest that reductions in collagen undulation are an adaptive mechanism competing against LVFW thinning. Additionally, we show that the inclusion of subject-specific geometry when modeling myocardial tissue is essential for accurate prediction of tissue kinematics. Our approach serves as an essential step toward identifying fiber-level remodeling indices that govern the transition of MI to systolic heart failure. These indices complement the traditional, organ-level measures of LV anatomy and function that often fall short of early prognostication of heart failure in MI. In addition, our approach offers an integrated methodology to advance the design of personalized interventions, such as hydrogel injection, to reinforce and suppress native adaptive and maladaptive mechanisms, respectively, to prevent the transition of MI to heart failure. STATEMENT OF SIGNIFICANCE: Biomechanical and architectural adaptation of the LVFW remains a central, yet overlooked, remodeling process post-MI. Our study indicates the biomechanical adaptation of the LVFW post-MI is highly individual and driven by altered fiber network architecture and collective changes in collagen fiber content, undulation, and stiffness. Our findings demonstrate the possibility of using cardiac strains to infer such fiber-level remodeling events through in-silico modeling, paving the way for in-vivo characterization of multiscale biomechanical indices in humans. Such indices will complement the traditional, organ-level measures of LV anatomy and function that often fall short of early prognostication of heart failure in MI.

The non-affine fiber network solver: A multiscale fiber network material model for finite-element analysis

Multiscale mechanical models in biomaterials research have largely relied on simplifying the microstructure in order to make large-scale simulations tractable. The microscale simplifications often rely on approximations of the constituent distributions and assumptions on the deformation of the constituents. Of particular interest in biomechanics are fiber embedded materials, where simplified fiber distributions and assumed affinity in the fiber deformation greatly influence the mechanical behavior. The consequences of these assumptions are problematic when dealing with microscale mechanical phenomena such as cellular mechanotransduction in growth and remodeling, and fiber-level failure events during tissue failure. In this work, we propose a technique for coupling non-affine network models to finite element solvers, allowing for simulation of discrete microstructural phenomena within macroscopically complex geometries. The developed plugin is readily available as an open-source library for use with the bio-focused finite element software FEBio, and the description of the implementation allows for the adaptation to other finite element solvers.

Fourier analysis of collagen bundle orientation in myocardial infarction scars

Collagen bundle orientation (CBO) in myocardial infarct scars plays a major role in scar mechanics and complications after infarction. We aim to compare four histopathological methods for CBO measurement in myocardial scarring. Myocardial infarction was induced in 21 pigs by balloon coronary occlusion. Scar samples were obtained at 4 weeks, stained with Masson’s trichrome, Picrosirius red, and Hematoxylin–Eosin (H&E), and photographed using light, polarized light microscopy, and confocal microscopy, respectively. Masson’s trichrome images were also optimized to remove non-collagenous structures. Two observers measured CBO by means of a semi-automated, Fourier analysis protocol. Interrater reliability and comparability between techniques were studied by the intraclass correlation coefficient (ICC) and Bland–Altman (B&A) plots and limits of agreement. Fourier analysis showed an almost perfect interrater reliability for each technique (ICC ≥ 0.95, p < 0.001 in all cases). CBO showed more randomly oriented values in Masson’s trichrome and worse comparability with other techniques (ICC vs. Picrosirius red: 0.79 [0.47–0.91], p = 0.001; vs. H&E-confocal: 0.70 [0.26–0.88], p = 0.005). However, optimized Masson’s trichrome showed almost perfect agreement with Picrosirius red (ICC 0.84 [0.6–0.94], p < 0.001) and H&E-confocal (ICC 0.81 [0.54–0.92], p < 0.001), as well as these latter techniques between each other (ICC 0.84 [0.60–0.93], p < 0.001). In summary, a semi-automated, Fourier-based method can provide highly reproducible CBO measurements in four different histopathological techniques. Masson’s trichrome tends to provide more randomly oriented CBO index values, probably due to non-specific visualization of non-collagenous structures. However, optimization of Masson’s trichrome microphotographs to remove non-collagenous components provides an almost perfect comparability between this technique, Picrosirius red and H&E-confocal.

Fabrication and characterization methods for investigating cell-matrix interactions in environments possessing spatial orientation heterogeneity

Fibrillar collagen is a ubiquitous structural protein that plays a significant role in determining the mechanical properties of various tissues. The constituent collagen architecture can give direct insight into the respective functional role of the tissue due to the strong structure-function relationship that is exhibited. In such tissues, matrix structure can vary across local subregions contributing to mechanical heterogeneity which can be implicated in tissue function or failure. The post-myocardial infarction scar environment is an example of note where mechanically insufficient collagen can result in impaired cardiac function and possibly tissue rupture due to post-MI cellular response and matrix interactions. In order to further develop the understanding of cell-matrix interactions within heterogeneous environments, we developed a method of heterogeneous collagen gel fabrication which produces a region of randomly oriented fibers directly adjacent to an interconnected region of anisotropic alignment. To fully capture and evaluate the degree of alignment and spatial orientation heterogeneity, several image processing and automated analysis methods were employed. Our analysis revealed the successful fabrication of an interconnected spatially heterogeneous collagen gel possessing distinct regions of random or preferential alignment. Additionally, embedded cell populations were observed to recognize and reorient with their underlying and surrounding architectures through our cell-centric analysis techniques. STATEMENT OF SIGNIFICANCE: Fibrillar collagen is a structural protein that contributes to the architecture-function relationship exhibited by various tissues where mechanically insufficient collagen architecture can lead to tissue failure. One environment where this can occur is the post-myocardial infarction scar environment where too much or too little collagen accumulation coupled with spatial fiber orientation heterogeneity can lead to environments incapable of normal mechanical functionality. While there are methodologies capable of generating aligned constructs, they do so with varying degrees of control and complexity with many producing uniform construct alignment. The presented platform is simple and produces continuous constructs possessing inherent spatial orientation heterogeneity. Coupling this with image processing and automated analysis methods enables the probing of fundamental cell-matrix interactions within heterogeneous environments.

Does the Heart Want What It Wants? A Case for Self-Adapting, Mechano-Sensitive Therapies After Infarction

There is a critical need for interventions to control the development and remodeling of scar tissue after myocardial infarction. A significant hurdle to fibrosis-related therapy is presented by the complex spatial needs of the infarcted ventricle, namely that collagenous buildup is beneficial in the ischemic zone but detrimental in the border and remote zones. As a new, alternative approach, we present a case to develop self-adapting, mechano-sensitive drug targets in order to leverage local, microenvironmental mechanics to modulate a therapy's pharmacologic effect. Such approaches could provide self-tuning control to either promote fibrosis or reduce fibrosis only when and where it is beneficial to do so.

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.

Biomechanics of infarcted left ventricle: a review of modelling

Mathematical modelling in biomechanics of infarcted left ventricle (LV) serves as an indispensable tool for remodelling mechanism exploration, LV biomechanical property estimation and therapy assessment after myocardial infarction (MI). However, a review of mathematical modelling after MI has not been seen in the literature so far. In the paper, a systematic review of mathematical models in biomechanics of infarcted LV was established. The models include comprehensive cardiovascular system model, essential LV pressure-volume and stress-stretch models, constitutive laws for passive myocardium and scars, tension models for active myocardium, collagen fibre orientation optimization models, fibroblast and collagen fibre growth/degradation models and integrated growth-electro-mechanical model after MI. The primary idea, unique characteristics and key equations of each model were identified and extracted. Discussions on the models were provided and followed research issues on them were addressed. Considerable improvements in the cardiovascular system model, LV aneurysm model, coupled agent-based models and integrated electro-mechanical-growth LV model are encouraged. Substantial attention should be paid to new constitutive laws with respect to stress-stretch curve and strain energy function for infarcted passive myocardium, collagen fibre orientation optimization in scar, cardiac rupture and tissue damage and viscoelastic effect post-MI in the future.

Ex Vivo Mechanical Tests and Multiscale Computational Modeling Highlight the Importance of Intramural Shear Stress in Ascending Thoracic Aortic Aneurysms

Ascending thoracic aortic aneurysms (ATAAs) are anatomically complex in terms of architecture and geometry, and both complexities contribute to unpredictability of ATAA dissection and rupture in vivo. The goal of this work was to examine the mechanism of ATAA failure using a combination of detailed mechanical tests on human tissue and a multiscale computational model. We used (1) multiple, geometrically diverse, mechanical tests to characterize tissue properties; (2) a multiscale computational model to translate those results into a broadly usable form; and (3) a model-based computer simulation of the response of an ATAA to the stresses generated by the blood pressure. Mechanical tests were performed in uniaxial extension, biaxial extension, shear lap, and peel geometries. ATAA tissue was strongest in circumferential extension and weakest in shear, presumably because of the collagen and elastin in the arterial lamellae. A multiscale, fiber-based model using different fiber properties for collagen, elastin, and interlamellar connections was specified to match all of the experimental data with one parameter set. Finally, this model was used to simulate ATAA inflation using a realistic geometry. The predicted tissue failure occurred in regions of high stress, as expected; initial failure events involved almost entirely interlamellar connections, consistent with arterial dissection-the elastic lamellae remain intact, but the connections between them fail. The failure of the interlamellar connections, paired with the weakness of the tissue under shear loading, is suggestive that shear stress within the tissue may contribute to ATAA dissection.