A model for arterial adaptation combining microstructural collagen remodeling and 3D tissue growth

Mechanobiological tissue instability induced by stress-modulated growth

The growth of biological tissues, which is regulated by a variety of factors, can induce stresses that may, in turn, destabilize the tissues into diverse patterns. In most previous studies, however, tissue growth was usually assumed as a prescribed parameter independent of stresses, limiting our understanding of the mechanobiological morphogenesis of real tissues. In this paper, we propose a theoretical model to investigate the mechanobiological response of soft tissues undergoing stress-modulated growth. Linear stability analysis is first performed to elucidate the surface instability mechanism induced by stress-modulated volumetric growth. We further conduct finite element simulations to validate the theoretical prediction and, particularly, to capture the post-buckling pattern evolution. Our results show that the non-uniform stresses, which evolve with the tissue growth and morphogenesis, exert mechanical feedback on the growth itself, producing up-down asymmetric surface morphologies as observed in, for example, the gyrification of human brains and brain organoids. It is also revealed that large residual stresses are unnecessary to cause mechanobiological instability and subsequent asymmetric patterning, which has long been believed to be driven by sufficiently high stresses. The present work could help us to understand the morphological changes of biological tissues under physiological and pathological conditions.

Microstructure-based numerical simulation of the mechanical behaviour of ocular tissue

This paper aims to present a novel full-eye biomechanical material model that incorporates the characteristics of ocular tissues at microstructural level, and use the model to analyse the age-related stiffening in tissue behaviour. The collagen content in ocular tissues, as obtained using X-ray scattering measurements, was represented by sets of Zernike polynomials that covered both the cornea and sclera, then used to reconstruct maps of collagen fibril magnitude and orientation on the three-dimensional geometry of the eye globe. Fine-mesh finite-element (FE) models with eye-specific geometry were built and supported by a user-defined material model (UMAT), which considered the regional variation of fibril density and orientation. The models were then used in an iterative inverse modelling study to derive the material parameters that represent the experimental behaviour of ocular tissues from donors aged between 50 and 90 years obtained in earlier ex vivo studies. Sensitivity analysis showed that reducing the number of directions that represented the anisotropy of collagen fibril orientation at each X-ray scattering measurement point from 180 to 16 would have limited and insignificant effect on the FE solution (0.08%). Inverse analysis resulted in material parameters that provided a close match with experimental intraocular pressure-deformation behaviour with a root mean square of error between 3.6% and 4.3%. The results also demonstrated a steady increase in mechanical stiffness in all ocular regions with age. A constitutive material model based on distributions of collagen fibril density and orientation has been developed to enable the accurate representation of the biomechanical behaviour of ocular tissues. The model offers a high level of control of stiffness and anisotropy across ocular globe, and therefore has the potential for use in planning surgical and medical procedures.

Computational modeling of cardiac growth and remodeling in pressure overloaded hearts-Linking microstructure to organ phenotype

Cardiac growth and remodeling (G&R) refers to structural changes in myocardial tissue in response to chronic alterations in loading conditions. One such condition is pressure overload where elevated wall stresses stimulate the growth in cardiomyocyte thickness, associated with a phenotype of concentric hypertrophy at the organ scale, and promote fibrosis. The initial hypertrophic response can be considered adaptive and beneficial by favoring myocyte survival, but over time if pressure overload conditions persist, maladaptive mechanisms favoring cell death and fibrosis start to dominate, ultimately mediating the transition towards an overt heart failure phenotype. The underlying mechanisms linking biological factors at the myocyte level to biomechanical factors at the systemic and organ level remain poorly understood. Computational models of G&R show high promise as a unique framework for providing a quantitative link between myocardial stresses and strains at the organ scale to biological regulatory processes at the cellular level which govern the hypertrophic response. However, microstructurally motivated, rigorously validated computational models of G&R are still in their infancy. This article provides an overview of the current state-of-the-art of computational models to study cardiac G&R. The microstructure and mechanosensing/mechanotransduction within cells of the myocardium is discussed and quantitative data from previous experimental and clinical studies is summarized. We conclude with a discussion of major challenges and possible directions of future research that can advance the current state of cardiac G&R computational modeling. STATEMENT OF SIGNIFICANCE: The mechanistic links between organ-scale biomechanics and biological factors at the cellular size scale remain poorly understood as these are largely elusive to investigations using experimental methodology alone. Computational G&R models show high promise to establish quantitative links which allow more mechanistic insight into adaptation mechanisms and may be used as a tool for stratifying the state and predict the progression of disease in the clinic. This review provides a comprehensive overview of research in this domain including a summary of experimental data. Thus, this study may serve as a basis for the further development of more advanced G&R models which are suitable for making clinical predictions on disease progression or for testing hypotheses on pathogenic mechanisms using in-silico models.

Multi-Scale Modeling of Vision-Guided Remodeling and Age-Dependent Growth of the Tree Shrew Sclera During Eye Development and Lens-Induced Myopia

The sclera uses unknown mechanisms to match the eye's axial length to its optics during development, producing eyes with good focus (emmetropia). A myopic eye is too long for its own optics. We propose a multi-scale computational model to simulate eye development based on the assumption that scleral growth is controlled by genetic factors while scleral remodeling is driven by genetic factors and the eye's refractive error. We define growth as a mechanism that changes the tissue volume and mass while remodeling involves internal micro-deformations that are volume-preserving at the macroscale. The model was fitted against longitudinal refractive measurements in tree shrews of different ages and exposed to three different visual conditions: (i) normal development; (ii) negative lens wear to induce myopia; and (iii) recovery from myopia by removing the negative lens. The model was able to replicate the age- and vision-dependent response of the tree shrew experiments. Scleral growth ceased at younger age than scleral remodeling. The remodeling rate decreased as the eye emmetropized but increased at any age when a negative lens was put on. The predictive power of the model was investigated by calculating the susceptibility to scleral remodeling and the response to form deprivation myopia in tree shrews. Both predictions were in good agreement with experimental data that were not used to fit the model. We propose the first model that distinguishes scleral growth from remodeling. The good agreement of our results with experimental data supports the notion that scleral growth and scleral remodeling are two independently controlled mechanisms during eye development.

Growth Description for Vessel Wall Adaptation: A Thick-Walled Mixture Model of Abdominal Aortic Aneurysm Evolution

(1) Background: Vascular tissue seems to adapt towards stable homeostatic mechanical conditions, however, failure of reaching homeostasis may result in pathologies. Current vascular tissue adaptation models use many ad hoc assumptions, the implications of which are far from being fully understood; (2) Methods: The present study investigates the plausibility of different growth kinematics in modeling Abdominal Aortic Aneurysm (AAA) evolution in time. A structurally motivated constitutive description for the vessel wall is coupled to multi-constituent tissue growth descriptions; Constituent deposition preserved either the constituent’s density or its volume, and Isotropic Volume Growth (IVG), in-Plane Volume Growth (PVG), in-Thickness Volume Growth (TVG) and No Volume Growth (NVG) describe the kinematics of the growing vessel wall. The sensitivity of key modeling parameters is explored, and predictions are assessed for their plausibility; (3) Results: AAA development based on TVG and NVG kinematics provided not only quantitatively, but also qualitatively different results compared to IVG and PVG kinematics. Specifically, for IVG and PVG kinematics, increasing collagen mass production accelerated AAA expansion which seems counterintuitive. In addition, TVG and NVG kinematics showed less sensitivity to the initial constituent volume fractions, than predictions based on IVG and PVG; (4) Conclusions: The choice of tissue growth kinematics is of crucial importance when modeling AAA growth. Much more interdisciplinary experimental work is required to develop and validate vascular tissue adaption models, before such models can be of any practical use.

Emergence of Collagen Orientation Heterogeneity in Healing Infarcts and an Agent-Based Model

Spatial heterogeneity of matrix structure can be an important determinant of tissue function. Although bulk properties of collagen structure in healing myocardial infarcts have been characterized previously, regional heterogeneity in infarct structure has received minimal attention. Herein, we quantified regional variations of collagen and nuclear orientations over the initial weeks of healing after infarction in rats, and employed a computational model of infarct remodeling to test potential explanations for the heterogeneity we observed in vivo. Fiber and cell orientation maps were generated from infarct samples acquired previously at 1, 2, 3, and 6 weeks postinfarction in a rat ligation model. We analyzed heterogeneity by calculating the dot product of each fiber or cell orientation vector with every other fiber or cell orientation vector, and plotting that dot product versus distance between the fibers or cells. This analysis revealed prominent regional heterogeneity, with alignment of both fibers and cell nuclei in local pockets far exceeding the global average. Using an agent-based model of fibroblast-mediated collagen remodeling, we found that similar levels of heterogeneity can spontaneously emerge from initially isotropic matrix via locally reinforcing cell-matrix interactions. Specifically, cells that sensed fiber orientation at a distance or remodeled fibers at a distance by traction-mediated reorientation or aligned deposition gave rise to regionally heterogeneous structures. However, only the simulations in which cells deposited collagen fibers aligned with their own orientation reproduced experimentally measured patterns of heterogeneity across all time points. These predictions warrant experimental follow-up to test the role of such mechanisms in vivo and identify opportunities to control heterogeneity for therapeutic benefit.

A validated software application to measure fiber organization in soft tissue

The mechanical behavior of soft connective tissue is governed by a dense network of fibrillar proteins in the extracellular matrix. Characterization of this fibrous network requires the accurate extraction of descriptive structural parameters from imaging data, including fiber dispersion and mean fiber orientation. Common methods to quantify fiber parameters include fast Fourier transforms (FFT) and structure tensors; however, information is limited on the accuracy of these methods. In this study, we compared these two methods using test images of fiber networks with varying topology. The FFT method with a band-pass filter was the most accurate, with an error of [Formula: see text] in measuring mean fiber orientation and an error of [Formula: see text] in measuring fiber dispersion in the test images. The accuracy of the structure tensor method was approximately five times worse than the FFT band-pass method when measuring fiber dispersion. A free software application, FiberFit, was then developed that utilizes an FFT band-pass filter to fit fiber orientations to a semicircular von Mises distribution. FiberFit was used to measure collagen fibril organization in confocal images of bovine ligament at magnifications of [Formula: see text] and [Formula: see text]. Grayscale conversion prior to FFT analysis gave the most accurate results, with errors of [Formula: see text] for mean fiber orientation and [Formula: see text] for fiber dispersion when measuring confocal images at [Formula: see text]. By developing and validating a software application that facilitates the automated analysis of fiber organization, this study can help advance a mechanistic understanding of collagen networks and help clarify the mechanobiology of soft tissue remodeling and repair.

A 3D MODEL FOR MURAL-CELL-MEDIATED DESTRUCTIVE REMODELING DURING EARLY DEVELOPMENT OF A CEREBRAL ANEURYSM

Development of a diagnostic tool for predicting the behavior of cerebral aneurysms was the inspiration of many research groups in recent years. In the present study a fluid–solid-growth (FSG) model for the early development of a cerebral aneurysm was presented in a 3D model of the internal carotid artery (ICA). This model is the result of two parallel mechanisms: first, defining arterial wall as a living tissue with the ability of degradation, growth and remodeling and second, full coupling of the wall and the blood flow. Taking into account the shear dependent nature of elastin degradation and mural-cell-mediated destructive activities, here, the degradation process has been linked to high effective stress of the vascular wall. The evolving properties of the elastinous and collagenous constituents have been predicted during the early development of the aneurysm and the code is applicable to more complicated aneurismal growth models.

Bio-Chemo-Mechanical Models of Vascular Mechanics

Models of vascular mechanics are necessary to predict the response of an artery under a variety of loads, for complex geometries, and in pathological adaptation. Classic constitutive models for arteries are phenomenological and the fitted parameters are not associated with physical components of the wall. Recently, microstructurally-linked models have been developed that associate structural information about the wall components with tissue-level mechanics. Microstructurally-linked models are useful for correlating changes in specific components with pathological outcomes, so that targeted treatments may be developed to prevent or reverse the physical changes. However, most treatments, and many causes, of vascular disease have chemical components. Chemical signaling within cells, between cells, and between cells and matrix constituents affects the biology and mechanics of the arterial wall in the short- and long-term. Hence, bio-chemo-mechanical models that include chemical signaling are critical for robust models of vascular mechanics. This review summarizes bio-mechanical and bio-chemo-mechanical models with a focus on large elastic arteries. We provide applications of these models and challenges for future work.

Microstructural modelling of cerebral aneurysm evolution through effective stress mediated destructive remodelling

Recently, researchers have shown an increased interest in the biomechanical modelling of cerebral aneurysm development. In the present study a fluid-solid-growth model for the formation of a fusiform aneurysm has been presented in an axi-symmetric geometry of the internal carotid artery. This model is the result of two parallel mechanisms: first, defining arterial wall as a living tissue with the ability of degradation, growth and remodelling and second, full coupling of the wall and the blood flow. Here for the first time the degradation of elastin has been defined as a function of vascular wall effective stress to take into account the shear dependent nature of degradation and the mural-cell-mediated destructive activities. The model has been stabilized in size and mechanical properties and is consistent with other computational or clinical studies. Furthermore, the evolving microstructural properties of the wall during the evolution process have been predicted.

Continuum mixture models of biological growth and remodeling: past successes and future opportunities

Biological growth processes involve mass exchanges that increase, decrease, or replace material that constitutes cells, tissues, and organs. In most cases, such exchanges alter the structural makeup of the material and consequently affect associated mechanobiological responses to applied loads. Given that the type and extent of changes in structural integrity depend on the different constituents involved (e.g., particular cytoskeletal or extracellular matrix proteins), the continuum theory of mixtures is ideally suited to model the mechanics of growth and remodeling. The goal of this review is twofold: first, to highlight a few illustrative examples that show diverse applications of mixture theory to describe biological growth and/or remodeling; second, to identify some open problems in the fields of modeling soft-tissue growth and remodeling.

Perspectives on biomechanical growth and remodeling mechanisms in glaucoma()

Glaucoma is a blinding diseases in which damage to the axons results in loss of retinal ganglion cells. Experimental evidence indicates that chronic intraocular pressure elevation initiates axonal insult at the level of the lamina cribrosa. The lamina cribrosa is a porous collagen structure through which the axons pass on their path from the retina to the brain. Recent experimental studies revealed the extensive structural changes of the lamina cribrosa and its surrounding tissues during the development and progression of glaucoma. In this perspective paper we review the experimental evidence for growth and remodeling mechanisms in glaucoma including adaptation of tissue anisotropy, tissue thickening/thinning, tissue elongation/shortening and tissue migration. We discuss the existing predictive computational approaches that try to elucidate the potential biomechanical basis of theses growth and remodeling mechanisms and highlight open questions, challenges, and avenues for further development.

Constrained Mixture Models as Tools for Testing Competing Hypotheses in Arterial Biomechanics: A Brief Survey

Hypothesis testing via numerical models has emerged as a powerful tool which permits the verification of theoretical frameworks against canonical experimental and clinical observations. Cleverly designed computational experiments also inspire new methodologies by elucidating important biological processes and restricting parametric spaces. Constrained mixture models of arterial growth and remodeling (G&R) can facilitate the design of computational experiments which can bypass technical limitations in the laboratory, by considering illustrative special cases. The resulting data may then inform the design of focused experimental techniques and the development of improved theories. This work is a survey of computational hypothesis-testing studies, which exploit the unique abilities offered by the constrained mixture theory of arterial G&R. Specifically, we explore the core hypotheses integrated in these models, review their basic mathematical conceptualizations, and recapitulate their most salient and illuminating findings. We then assess how a decade's worth of constrained mixture models have contributed to a lucid, emerging picture of G&R mechanisms.

Mechanics, mechanobiology, and modeling of human abdominal aorta and aneurysms

Biomechanical factors play fundamental roles in the natural history of abdominal aortic aneurysms (AAAs) and their responses to treatment. Advances during the past two decades have increased our understanding of the mechanics and biology of the human abdominal aorta and AAAs, yet there remains a pressing need for considerable new data and resulting patient-specific computational models that can better describe the current status of a lesion and better predict the evolution of lesion geometry, composition, and material properties and thereby improve interventional planning. In this paper, we briefly review data on the structure and function of the human abdominal aorta and aneurysmal wall, past models of the mechanics, and recent growth and remodeling models. We conclude by identifying open problems that we hope will motivate studies to improve our computational modeling and thus general understanding of AAAs.

A Multilayered Wall Model of Arterial Growth and Remodeling

Adaptations of large arteries to sustained alterations in hemodynamics that cause changes in both caliber and stiffness are increasingly recognized as important initiators or indicators of cardiovascular risk to high flow, low resistance organs such as the brain, heart, and kidney. There is, therefore, a pressing need to understand better the underlying causes of geometric and material adaptations by large arteries and the associated time courses. Although such information must ultimately come from well designed experiments, mathematical models will continue to play a vital role in the design of these experiments and their interpretation. In this paper, we present a new multilayered model of the time course of basilar artery growth and remodeling in response to sustained alterations in blood pressure and flow. We show, for example, that single- and multi-layered models consistently predict similar changes in caliber and wall thickness, but multilayered models provide additional insight into other important metrics such as the residual stress related opening angle and the axial prestress, both of which are fundamental to arterial homeostasis and responses to injury or insult.

Lamina Cribrosa Thickening in Early Glaucoma Predicted by a Microstructure Motivated Growth and Remodeling Approach

Glaucoma is among the leading causes of blindness worldwide. The ocular disease is characterized by irreversible damage of the retinal ganglion cell axons at the level of the lamina cribrosa (LC). The LC is a porous, connective tissue structure whose function is believed to provide mechanical support to the axons as they exit the eye on their path from the retina to the brain. Early experimental glaucoma studies have shown that the LC remodels into a thicker, more posterior structure which incorporates more connective tissue after intraocular pressure (IOP) elevation. The process by which this occurs is unknown. Here we present a microstructure motivated growth and remodeling (G&R) formulation to explore a potential mechanism of these structural changes. We hypothesize that the mechanical strain experienced by the collagen fibrils in the LC stimulates the G&R response at the micro-scale. The proposed G&R algorithm controls collagen fibril synthesis/degradation and adapts the residual strains between collagen fibrils and the surrounding tissue to achieve biomechanical homeostasis. The G&R algorithm was applied to a generic finite element model of the human eye subjected to normal and elevated IOP. The G&R simulation underscores the biomechanical need for a LC at normal IOP. The numerical results suggest that IOP elevation leads to LC thickening due to an increase in collagen fibril mass, which is in good agreement with experimental observations in early glaucoma monkey eyes. This is the first study to demonstrate that a biomechanically-driven G&R mechanism can lead to the LC thickening observed in early experimental glaucoma.

Modeling of cardiac growth and remodeling of myofiber orientation

The heart has the ability to respond to long-term changes in its environment through changes in mass (growth), shape (morphogenesis) and tissue properties (remodeling). For improved quantitative understanding of cardiac growth and remodeling (G&R) experimental studies need to be complemented by mathematical models. This paper reviews models for cardiac growth and remodeling of myofiber orientation, as induced by mechanical stimuli. A distinction is made between optimization models, that focus on the end stage of G&R, and adaptation models, that aim to more closely describe the mechanistic relation between stimulus and effect. While many models demonstrate qualitatively promising results, a lot of questions remain, e.g. with respect to the choice of the stimulus for G&R or the long-term stability of the outcome of the model. A continued effort combining information on mechanotransduction at the cellular level, experimental observations on G&R at organ level, and testing of hypotheses on stimulus-effect relations in mathematical models is needed to answer these questions on cardiac G&R. Ultimately, models of cardiac G&R seem indispensable for patient-specific modeling, both to reconstruct the actual state of the heart and to assess the long-term effect of potential interventions.

A structural constitutive model considering angular dispersion and waviness of collagen fibres of rabbit facial veins

Background

Structural constitutive models of vascular wall integrate information on composition and structural arrangements of tissue. In blood vessels, collagen fibres are arranged in coiled and wavy bundles and the individual collagen fibres have a deviation from their mean orientation. A complete structural constitutive model for vascular wall should incorporate both waviness and orientational distribution of fibres. We have previously developed a model, for passive properties of vascular wall, which considers the waviness of collagen fibres. However, to our knowledge there is no structural model of vascular wall which integrates both these features.

Methods

In this study, we have suggested a structural strain energy function that incorporates not only the waviness but also the angular dispersion of fibres. We studied the effect of parameters related to the orientational distribution on macro-mechanical behaviour of tissue during inflation-extension tests. The model was further applied on experimental data from rabbit facial veins.

Results

Our parametric study showed that the model is less sensitive to the orientational dispersion when fibres are mainly oriented circumferentially. The macro-mechanical response is less sensitive to changes in the mean orientation when fibres are more dispersed. The model accurately fitted the experimental data of veins, while not improving the quality of the fit compared to the model without dispersion. Our results showed that the orientational dispersion of collagen fibres could be compensated by a less abrupt and shifted to higher strain collagen engagement pattern. This should be considered when the model is fitted to experimental data and model parameters are used to study structural modifications of collagen fibre network in physiology and disease.

Conclusions

The presented model incorporates structural features related to waviness and orientational distribution of collagen fibres and thus offers possibilities to better understand the relation between structure and function in the vascular wall. Also, the model can be used to further study mechanically induced collagen remodelling in vascular tissue in health and disease.