A Mechanobiologically Equilibrated Constrained Mixture Model for Growth and Remodeling of Soft Tissues

Computational approaches for mechanobiology in cardiovascular development and diseases

The cardiovascular development in vertebrates evolves in response to genetic and mechanical cues. The dynamic interplay among mechanics, cell biology, and anatomy continually shapes the hydraulic networks, characterized by complex, non-linear changes in anatomical structure and blood flow dynamics. To better understand this interplay, a diverse set of molecular and computational tools has been used to comprehensively study cardiovascular mechanobiology. With the continual advancement of computational capacity and numerical techniques, cardiovascular simulation is increasingly vital in both basic science research for understanding developmental mechanisms and disease etiologies, as well as in clinical studies aimed at enhancing treatment outcomes. This review provides an overview of computational cardiovascular modeling. Beginning with the fundamental concepts of computational cardiovascular modeling, it navigates through the applications of computational modeling in investigating mechanobiology during cardiac development. Second, the article illustrates the utility of computational hemodynamic modeling in the context of treatment planning for congenital heart diseases. It then delves into the predictive potential of computational models for elucidating tissue growth and remodeling processes. In closing, we outline prevailing challenges and future prospects, underscoring the transformative impact of computational cardiovascular modeling in reshaping cardiovascular science and clinical practice.

A Fluid-Solid-Growth Solver for Cardiovascular Modeling

We implement full, three-dimensional constrained mixture theory for vascular growth and remodeling into a finite element fluid-structure interaction (FSI) solver. The resulting "fluid-solid-growth" (FSG) solver allows long term, patient-specific predictions of changing hemodynamics, vessel wall morphology, tissue composition, and material properties. This extension from short term (FSI) to long term (FSG) simulations increases clinical relevance by enabling mechanobioloigcally-dependent studies of disease progression in complex domains.

A homogenized constrained mixture model of cardiac growth and remodeling: analyzing mechanobiological stability and reversal

Cardiac growth and remodeling (G&R) patterns change ventricular size, shape, and function both globally and locally. Biomechanical, neurohormonal, and genetic stimuli drive these patterns through changes in myocyte dimension and fibrosis. We propose a novel microstructure-motivated model that predicts organ-scale G&R in the heart based on the homogenized constrained mixture theory. Previous models, based on the kinematic growth theory, reproduced consequences of G&R in bulk myocardial tissue by prescribing the direction and extent of growth but neglected underlying cellular mechanisms. In our model, the direction and extent of G&R emerge naturally from intra- and extracellular turnover processes in myocardial tissue constituents and their preferred homeostatic stretch state. We additionally propose a method to obtain a mechanobiologically equilibrated reference configuration. We test our model on an idealized 3D left ventricular geometry and demonstrate that our model aims to maintain tensional homeostasis in hypertension conditions. In a stability map, we identify regions of stable and unstable G&R from an identical parameter set with varying systolic pressures and growth factors. Furthermore, we show the extent of G&R reversal after returning the systolic pressure to baseline following stage 1 and 2 hypertension. A realistic model of organ-scale cardiac G&R has the potential to identify patients at risk of heart failure, enable personalized cardiac therapies, and facilitate the optimal design of medical devices.

A multiscale computational model of arterial growth and remodeling including Notch signaling

Blood vessels grow and remodel in response to mechanical stimuli. Many computational models capture this process phenomenologically, by assuming stress homeostasis, but this approach cannot unravel the underlying cellular mechanisms. Mechano-sensitive Notch signaling is well-known to be key in vascular development and homeostasis. Here, we present a multiscale framework coupling a constrained mixture model, capturing the mechanics and turnover of arterial constituents, to a cell-cell signaling model, describing Notch signaling dynamics among vascular smooth muscle cells (SMCs) as influenced by mechanical stimuli. Tissue turnover was regulated by both Notch activity, informed by in vitro data, and a phenomenological contribution, accounting for mechanisms other than Notch. This novel framework predicted changes in wall thickness and arterial composition in response to hypertension similar to previous in vivo data. The simulations suggested that Notch contributes to arterial growth in hypertension mainly by promoting SMC proliferation, while other mechanisms are needed to fully capture remodeling. The results also indicated that interventions to Notch, such as external Jagged ligands, can alter both the geometry and composition of hypertensive vessels, especially in the short term. Overall, our model enables a deeper analysis of the role of Notch and Notch interventions in arterial growth and remodeling and could be adopted to investigate therapeutic strategies and optimize vascular regeneration protocols.

Remodeling of the uterine artery during and early after pregnancy in the mouse

Pregnancy associates with dramatic changes in maternal cardiovascular physiology that ensure that the utero-placental circulation can support the developing fetus. Particularly striking is the marked flow-induced remodeling of uterine arteries during pregnancy and their recovery following birth. Whereas details are available in the literature on alterations in hemodynamics within and changes in the dimensions of uterine arteries during and following pregnancy in mice, we report here the first biaxial biomechanical phenotyping of these arteries during this dynamic period of growth and remodeling (G&R). To gain additional insight into the measured G&R, we also use a computational constrained mixture model to describe and predict findings, including simulations related to complications that may arise during pregnancy. It is found that dramatic pregnancy-induced remodeling of the uterine artery is largely, but not completely, reversed in the postpartum period, which appears to be driven by increases in collagen turnover among other intramural changes. By contrast, data on the remodeling of the ascending aorta, an elastic artery, reveal modest changes that are fully recovered postpartum. There is strong motivation to continue biomechanical studies on this critical aspect of women's health, which has heretofore not received appropriate consideration from the biomechanics community.

Aortic Remodeling Kinetics in Response to Coarctation-Induced Mechanical Perturbations

Background: Coarctation of the aorta (CoA; constriction of the proximal descending thoracic aorta) is among the most common congenital cardiovascular defects. Coarctation-induced mechanical perturbations trigger a cycle of mechano-transduction events leading to irreversible precursors of hypertension including arterial thickening, stiffening, and vasoactive dysfunction in proximal conduit arteries. This study sought to identify kinetics of the stress-mediated compensatory response leading to these alterations using a preclinical rabbit model of CoA. Methods: A prior growth and remodeling (G&R) framework was reformulated and fit to empirical measurements from CoA rabbits classified into one control and nine CoA groups of various severities and durations (n = 63, 5-11/group). Empirical measurements included Doppler ultrasound imaging, uniaxial extension testing, catheter-based blood pressure, and wire myography, yielding the time evolution of arterial thickening, stiffening, and vasoactive dysfunction required to fit G&R constitutive parameters. Results: Excellent agreement was observed between model predictions and observed patterns of arterial thickening, stiffening, and dysfunction among all CoA groups. For example, predicted vascular impairment was not significantly different from empirical observations via wire myography (p-value > 0.13). Specifically, 48% and 45% impairment was observed in smooth muscle contraction and endothelial-dependent relaxation, respectively, which were accurately predicted using the G&R model. Conclusions: The resulting G&R model, for the first time, allows for prediction of hypertension precursors at neonatal ages that is currently challenging to examine in preclinical models. These findings provide a validated computational tool for prediction of persistent arterial dysfunction and identification of revised severity-duration thresholds that may ultimately avoid hypertension from CoA.

Beyond CFD: Emerging methodologies for predictive simulation in cardiovascular health and disease

Physics-based computational models of the cardiovascular system are increasingly used to simulate hemodynamics, tissue mechanics, and physiology in evolving healthy and diseased states. While predictive models using computational fluid dynamics (CFD) originated primarily for use in surgical planning, their application now extends well beyond this purpose. In this review, we describe an increasingly wide range of modeling applications aimed at uncovering fundamental mechanisms of disease progression and development, performing model-guided design, and generating testable hypotheses to drive targeted experiments. Increasingly, models are incorporating multiple physical processes spanning a wide range of time and length scales in the heart and vasculature. With these expanded capabilities, clinical adoption of patient-specific modeling in congenital and acquired cardiovascular disease is also increasing, impacting clinical care and treatment decisions in complex congenital heart disease, coronary artery disease, vascular surgery, pulmonary artery disease, and medical device design. In support of these efforts, we discuss recent advances in modeling methodology, which are most impactful when driven by clinical needs. We describe pivotal recent developments in image processing, fluid-structure interaction, modeling under uncertainty, and reduced order modeling to enable simulations in clinically relevant timeframes. In all these areas, we argue that traditional CFD alone is insufficient to tackle increasingly complex clinical and biological problems across scales and systems. Rather, CFD should be coupled with appropriate multiscale biological, physical, and physiological models needed to produce comprehensive, impactful models of mechanobiological systems and complex clinical scenarios. With this perspective, we finally outline open problems and future challenges in the field.

Intracellular signaling control of mechanical homeostasis in the aorta

Mature arteries exhibit a preferred biomechanical state in health evidenced by a narrow range of intramural and wall shear stresses. When stresses are perturbed by changes in blood pressure or flow, homeostatic mechanisms tend to restore target values via altered contractility and/or cell and matrix turnover. In contrast, vascular disease associates with compromised homeostasis, hence we must understand mechanisms underlying mechanical homeostasis and its robustness. Here, we use a multiscale computational model wherein mechanosensitive intracellular signaling pathways drive arterial growth and remodeling. First, we identify an ensemble of cell-level parameterizations where tissue-level responses are well-regulated and adaptive to hemodynamic perturbations. The responsible mechanism is persistent multiscale negative feedback whereby mechanosensitive signaling drives mass turnover until homeostatic target stresses are reached. This demonstrates how robustness emerges despite inevitable cell and individual heterogeneity. Second, we investigate tissue-level effects of signaling node knockdowns (ATIR, ROCK, TGF[Formula: see text]RII, PDGFR, ERK1/2) and find general agreement with experimental reports of fault tolerance. Robustness against structural changes manifests via low engagement of the node under baseline stresses or compensatory multiscale feedback via upregulation of additional pathways. Third, we show how knockdowns affect collagen and smooth muscle turnover at baseline and with perturbed stresses. In several cases, basal production is not remarkably affected, but sensitivities to stress deviations, which influence feedback strength, are reduced. Such reductions can impair adaptive responses, consistent with previously reported aortic vulnerability despite grossly normal appearances. Reduced stress sensitivities thus form a candidate mechanism for how robustness is lost, enabling transitions from health towards disease.

Neural operator learning of heterogeneous mechanobiological insults contributing to aortic aneurysms

Thoracic aortic aneurysm (TAA) is a localized dilatation of the aorta that can lead to life-threatening dissection or rupture. In vivo assessments of TAA progression are largely limited to measurements of aneurysm size and growth rate. There is promise, however, that computational modelling of the evolving biomechanics of the aorta could predict future geometry and properties from initiating mechanobiological insults. We present an integrated framework to train a deep operator network (DeepONet)-based surrogate model to identify TAA contributing factors using synthetic finite-element-based datasets. For training, we employ a constrained mixture model of aortic growth and remodelling to generate maps of local aortic dilatation and distensibility for multiple TAA risk factors. We evaluate the performance of the surrogate model for insult distributions varying from fusiform (analytically defined) to complex (randomly generated). We propose two frameworks, one trained on sparse information and one on full-field greyscale images, to gain insight into a preferred neural operator-based approach. We show that this continuous learning approach can predict the patient-specific insult profile associated with any given dilatation and distensibility map with high accuracy, particularly when based on full-field images. Our findings demonstrate the feasibility of applying DeepONet to support transfer learning of patient-specific inputs to predict TAA progression.

In vivo development of tissue engineered vascular grafts: a fluid-solid-growth model

Methods of tissue engineering continue to advance, and multiple clinical trials are underway evaluating tissue engineered vascular grafts (TEVGs). Whereas initial concerns focused on suture retention and burst pressure, there is now a pressing need to design grafts to have optimal performance, including an ability to grow and remodel in response to changing hemodynamic loads. Toward this end, there is similarly a need for computational methods that can describe and predict the evolution of TEVG geometry, composition, and material properties while accounting for changes in hemodynamics. Although the ultimate goal is a fluid-solid-growth (FSG) model incorporating fully 3D growth and remodeling and 3D hemodynamics, lower fidelity models having high computational efficiency promise to play important roles, especially in the design of candidate grafts. We introduce here an efficient FSG model of in vivo development of a TEVG based on two simplifying concepts: mechanobiologically equilibrated growth and remodeling of the graft and an embedded control volume analysis of the hemodynamics. Illustrative simulations for a model Fontan conduit reveal the utility of this approach, which promises to be particularly useful in initial design considerations involving formal methods of optimization which otherwise add considerably to the computational expense.

Biomechanical and structural responses of the aorta to intermittent hypobaric hypoxia in a rat model

High altitude hypoxia is a condition experienced by diverse populations worldwide. In addition, several jobs require working shifts where workers are exposed to repetitive cycles of hypobaric hypoxia and normobaric normoxia. Currently, few is known about the biomechanical cardiovascular responses of this condition. In the present study, we investigate the cycle-dependent biomechanical effects of intermittent hypobaric hypoxia (IHH) on the thoracic aorta artery, in terms of both structure and function. To determine the vascular effects of IHH, functional, mechanical and histological approaches were carried out in the thoracic aorta artery, using uniaxial, pre-stretch, ring opening, myography, and histological tests. Three groups of rats were established: control (normobaric normoxia, NN), 4-cycles of intermittent hypoxia (short-term intermittent hypobaric hypoxia, STH), and 10-cycles of intermittent hypoxia (long-term intermittent hypobaric hypoxia, LTH). The pre-stretch and ring opening tests, aimed at quantifying residual strains of the tissues in longitudinal and circumferential directions, showed that the hypoxia condition leads to an increase in the longitudinal stretch and a marked decrease of the circumferential residual strain. The uniaxial mechanical tests were used to determine the elastic properties of the tissues, showing that a general stiffening process occurs during the early stages of the IH (STH group), specially leading to a significative increase in the high strain elastic modulus ([Formula: see text]) and an increasing trend of low strain elastic modulus ([Formula: see text]). In contrast, the LTH group showed a more control-like mechanical behavior. Myography test, used to assess the vasoactive function, revealed that IH induces a high sensitivity to vasoconstrictor agents as a function of hypoxic cycles. In addition, the aorta showed an increased muscle-dependent vasorelaxation on the LTH group. Histological tests, used to quantify the elastic fiber, nuclei, and geometrical properties, showed that the STH group presents a state of vascular fibrosis, with a significant increase in elastin content, and a tendency towards an increase in collagen fibers. In addition, advanced stages of IH (LTH), showed a vascular remodeling effect with a significant increase of internal and external diameters. Considering all the multidimensional vascular effects, we propose the existence of a long-term passive adaptation mechanism and vascular dysfunction as cycle-dependent effects of intermittent exposures to hypobaric hypoxia.

Mechano-regulated cell-cell signaling in the context of cardiovascular tissue engineering

Cardiovascular tissue engineering (CVTE) aims to create living tissues, with the ability to grow and remodel, as replacements for diseased blood vessels and heart valves. Despite promising results, the (long-term) functionality of these engineered tissues still needs improvement to reach broad clinical application. The functionality of native tissues is ensured by their specific mechanical properties directly arising from tissue organization. We therefore hypothesize that establishing a native-like tissue organization is vital to overcome the limitations of current CVTE approaches. To achieve this aim, a better understanding of the growth and remodeling (G&R) mechanisms of cardiovascular tissues is necessary. Cells are the main mediators of tissue G&R, and their behavior is strongly influenced by both mechanical stimuli and cell-cell signaling. An increasing number of signaling pathways has also been identified as mechanosensitive. As such, they may have a key underlying role in regulating the G&R of tissues in response to mechanical stimuli. A more detailed understanding of mechano-regulated cell-cell signaling may thus be crucial to advance CVTE, as it could inspire new methods to control tissue G&R and improve the organization and functionality of engineered tissues, thereby accelerating clinical translation. In this review, we discuss the organization and biomechanics of native cardiovascular tissues; recent CVTE studies emphasizing the obtained engineered tissue organization; and the interplay between mechanical stimuli, cell behavior, and cell-cell signaling. In addition, we review past contributions of computational models in understanding and predicting mechano-regulated tissue G&R and cell-cell signaling to highlight their potential role in future CVTE strategies.

Constrained Mixture Models of Soft Tissue Growth and Remodeling - Twenty Years After

Soft biological tissues compromise diverse cell types and extracellular matrix constituents, each of which can possess individual natural configurations, material properties, and rates of turnover. For this reason, mixture-based models of growth (changes in mass) and remodeling (change in microstructure) are well-suited for studying tissue adaptations, disease progression, and responses to injury or clinical intervention. Such approaches also can be used to design improved tissue engineered constructs to repair, replace, or regenerate tissues. Focusing on blood vessels as archetypes of soft tissues, this paper reviews a constrained mixture theory introduced twenty years ago and explores its usage since by contrasting simulations of diverse vascular conditions. The discussion is framed within the concept of mechanical homeostasis, with consideration of solid-fluid interactions, inflammation, and cell signaling highlighting both past accomplishments and future opportunities as we seek to understand better the evolving composition, geometry, and material behaviors of soft tissues under complex conditions.

A Chemomechanobiological Model of the Long-Term Healing Response of Arterial Tissue to a Clamping Injury

Vascular clamping often causes injury to arterial tissue, leading to a cascade of cellular and extracellular events. A reliable in silico prediction of these processes following vascular injury could help us to increase our understanding thereof, and eventually optimize surgical techniques or drug delivery to minimize the amount of long-term damage. However, the complexity and interdependency of these events make translation into constitutive laws and their numerical implementation particularly challenging. We introduce a finite element simulation of arterial clamping taking into account acute endothelial denudation, damage to extracellular matrix, and smooth muscle cell loss. The model captures how this causes tissue inflammation and deviation from mechanical homeostasis, both triggering vascular remodeling. A number of cellular processes are modeled, aiming at restoring this homeostasis, i.e., smooth muscle cell phenotype switching, proliferation, migration, and the production of extracellular matrix. We calibrated these damage and remodeling laws by comparing our numerical results to in vivo experimental data of clamping and healing experiments. In these same experiments, the functional integrity of the tissue was assessed through myograph tests, which were also reproduced in the present study through a novel model for vasodilator and -constrictor dependent smooth muscle contraction. The simulation results show a good agreement with the in vivo experiments. The computational model was then also used to simulate healing beyond the duration of the experiments in order to exploit the benefits of computational model predictions. These results showed a significant sensitivity to model parameters related to smooth muscle cell phenotypes, highlighting the pressing need to further elucidate the biological processes of smooth muscle cell phenotypic switching in the future.

From Transcript to Tissue: Multiscale Modeling from Cell Signaling to Matrix Remodeling

Tissue-level biomechanical properties and function derive from underlying cell signaling, which regulates mass deposition, organization, and removal. Here, we couple two existing modeling frameworks to capture associated multiscale interactions-one for vessel-level growth and remodeling and one for cell-level signaling-and illustrate utility by simulating aortic remodeling. At the vessel level, we employ a constrained mixture model describing turnover of individual wall constituents (elastin, intramural cells, and collagen), which has proven useful in predicting diverse adaptations as well as disease progression using phenomenological constitutive relations. Nevertheless, we now seek an improved mechanistic understanding of these processes; we replace phenomenological relations in the mixture model with a logic-based signaling model, which yields a system of ordinary differential equations predicting changes in collagen synthesis, matrix metalloproteinases, and cell proliferation in response to altered intramural stress, wall shear stress, and exogenous angiotensin II. This coupled approach promises improved understanding of the role of cell signaling in achieving tissue homeostasis and allows us to model feedback between vessel mechanics and cell signaling. We verify our model predictions against data from the hypertensive murine infrarenal abdominal aorta as well as results from validated phenomenological models, and consider effects of noisy signaling and heterogeneous cell populations.

[Review of studies on the biomechanical modelling of the coupling effect between stent degradation and blood vessel remodeling]

The dynamic coupling of stent degradation and vessel remodeling can influence not only the structural morphology and material property of stent and vessel, but also the development of in-stent restenosis. The research achievements of biomechanical modelling and analysis of stent degradation and vessel remodeling were reviewed; several noteworthy research perspectives were addressed, a stent-vessel coupling model was developed based on stent damage function and vessel growth function, and then concepts of matching ratio and risk factor were established so as to evaluate the treatment effect of stent intervention, which may lay the scientific foundation for the structure design, mechanical analysis and clinical application of biodegradable stent.

Numerical knockouts-In silico assessment of factors predisposing to thoracic aortic aneurysms

Myriad risk factors–including uncontrolled hypertension, aging, and diverse genetic mutations–contribute to the development and enlargement of thoracic aortic aneurysms. Detailed analyses of clinical data and longitudinal studies of murine models continue to provide insight into the natural history of these potentially lethal conditions. Yet, because of the co-existence of multiple risk factors in most cases, it has been difficult to isolate individual effects of the many different factors or to understand how they act in combination. In this paper, we use a data-informed computational model of the initiation and progression of thoracic aortic aneurysms to contrast key predisposing risk factors both in isolation and in combination; these factors include localized losses of elastic fiber integrity, aberrant collagen remodeling, reduced smooth muscle contractility, and dysfunctional mechanosensing or mechanoregulation of extracellular matrix along with superimposed hypertension and aortic aging. In most cases, mild-to-severe localized losses in cellular function or matrix integrity give rise to varying degrees of local dilatations of the thoracic aorta, with enlargement typically exacerbated in cases wherein predisposing risk factors co-exist. The simulations suggest, for the first time, that effects of compromised smooth muscle contractility are more important in terms of dysfunctional mechanosensing and mechanoregulation of matrix than in vessel-level control of diameter and, furthermore, that dysfunctional mechanobiological control can yield lesions comparable to those in cases of compromised elastic fiber integrity. Particularly concerning, therefore, is that loss of constituents such as fibrillin-1, as in Marfan syndrome, can compromise both elastic fiber integrity and mechanosensing.

Aneurysms are local dilatations of the arterial wall that are responsible for significant disability and death. Detailed analyses of clinical data continue to provide insight into the natural history of these potentially lethal conditions, with myriad risk factors–including uncontrolled hypertension, aging, and diverse genetic mutations–contributing to their development and enlargement. Yet, because of the co-existence of these risk factors in most cases, it has been difficult to isolate individual effects or to understand how they act in combination. In this paper, we use a computational model of the initiation and progression of thoracic aortic aneurysms to contrast key predisposing factors both in isolation and in combination as well as with superimposed hypertension and aging. The present study recovers many findings from mouse models but with new and important observations that promise to guide in vivo and ex vivo studies as we seek to understand and eventually better treat these complex, multi-factorial lesions, with data-informed patient-specific computations eventually the way forward.

Vascular adaptation in the presence of external support - A modeling study

Vascular grafts have long been used to replace damaged or diseased vessels with considerable success, but a new approach is emerging where native vessels are merely supported, not replaced. Although external supports have been evaluated in diverse situations - ranging from aneurysmal disease to vein grafts or the Ross operation - optimal supports and procedures remain wanting. In this paper, we present a novel application of a growth and remodeling model well suited for parametrically exploring multiple designs of external supports while accounting for mechanobiological and immunobiological responses of the supported native vessel. These results suggest that a load bearing external support can reduce vessel thickening in response to pressure elevation. Results also suggest that the final adaptive state of the vessel depends on the structural stiffness of the support via a mechano-driven adaptation, although luminal encroachment may be a complication in the presence of chronic inflammation. Finally, the supported vessel can stiffen (structurally and materially) along circumferential and axial directions, which could have implications on overall hemodynamics and thus subsequent vascular remodeling. The proposed framework can provide valuable insights into vascular adaptation in the presence of external support, accelerate rational design, and aid translation of this emerging approach.

Next-generation tissue-engineered heart valves with repair, remodelling and regeneration capacity

Valvular heart disease is a major cause of morbidity and mortality worldwide. Surgical valve repair or replacement has been the standard of care for patients with valvular heart disease for many decades, but transcatheter heart valve therapy has revolutionized the field in the past 15 years. However, despite the tremendous technical evolution of transcatheter heart valves, to date, the clinically available heart valve prostheses for surgical and transcatheter replacement have considerable limitations. The design of next-generation tissue-engineered heart valves (TEHVs) with repair, remodelling and regenerative capacity can address these limitations, and TEHVs could become a promising therapeutic alternative for patients with valvular disease. In this Review, we present a comprehensive overview of current clinically adopted heart valve replacement options, with a focus on transcatheter prostheses. We discuss the various concepts of heart valve tissue engineering underlying the design of next-generation TEHVs, focusing on off-the-shelf technologies. We also summarize the latest preclinical and clinical evidence for the use of these TEHVs and describe the current scientific, regulatory and clinical challenges associated with the safe and broad clinical translation of this technology.

Computational modeling for cardiovascular tissue engineering: the importance of including cell behavior in growth and remodeling algorithms

Understanding cardiovascular growth and remodeling (G&R) is fundamental for designing robust cardiovascular tissue engineering strategies, which enable synthetic or biological scaffolds to transform into healthy living tissues after implantation. Computational modeling, particularly when integrated with experimental research, is key for advancing our understanding, predicting the in vivo evolution of engineered tissues, and efficiently optimizing scaffold designs. As cells are ultimately the drivers of G&R and known to change their behavior in response to mechanical cues, increasing efforts are currently undertaken to capture (mechano-mediated) cell behavior in computational models. In this selective review, we highlight some recent examples that are relevant in the context of cardiovascular tissue engineering and discuss the current and future biological and computational challenges for modeling cell-mediated G&R.

Cell signaling model for arterial mechanobiology

Arterial growth and remodeling at the tissue level is driven by mechanobiological processes at cellular and sub-cellular levels. Although it is widely accepted that cells seek to promote tissue homeostasis in response to biochemical and biomechanical cues—such as increased wall stress in hypertension—the ways by which these cues translate into tissue maintenance, adaptation, or maladaptation are far from understood. In this paper, we present a logic-based computational model for cell signaling within the arterial wall, aiming to predict changes in extracellular matrix turnover and cell phenotype in response to pressure-induced wall stress, flow-induced wall shear stress, and exogenous sources of angiotensin II, with particular interest in mouse models of hypertension. We simulate a number of experiments from the literature at both the cell and tissue level, involving single or combined inputs, and achieve high qualitative agreement in most cases. Additionally, we demonstrate the utility of this modeling approach for simulating alterations (in this case knockdowns) of individual nodes within the signaling network. Continued modeling of cellular signaling will enable improved mechanistic understanding of arterial growth and remodeling in health and disease, and will be crucial when considering potential pharmacological interventions.

Biological soft tissues are characterized by continuous production and removal of material, which endows them with a remarkable ability to adapt to changes in their biochemical and biomechanical environments. For arteries, mechanical stimuli result primarily from changes in blood pressure or flow, and biochemical changes are induced by multiple factors, including pharmacological intervention. In order to understand how arterial properties are maintained in health, or how they adapt or fail to adapt in disease, we must understand better how these diverse stimuli affect material turnover. Extracellular matrix is tightly regulated by mechano-sensing and mechano-regulation, and therefore cell signaling, thus we present a computational model of relevant signaling pathways within the vascular wall, with the aim of predicting changes in wall composition and function in response to three main inputs: pressure-induced wall stress, flow-induced wall shear stress, and exogenous angiotensin II. We obtain qualitative agreement with a range of experimental studies from the literature, and provide illustrative examples demonstrating how such models can be used to further our understanding of arterial remodeling.

Fast, Rate-Independent, Finite Element Implementation of a 3D Constrained Mixture Model of Soft Tissue Growth and Remodeling

Constrained mixture models of soft tissue growth and remodeling can simulate many evolving conditions in health as well as in disease and its treatment, but they can be computationally expensive. In this paper, we derive a new fast, robust finite element implementation based on a concept of mechanobiological equilibrium that yields fully resolved solutions and allows computation of quasi-equilibrated evolutions when imposed perturbations are slow relative to the adaptive process. We demonstrate quadratic convergence and verify the model via comparisons with semi-analytical solutions for arterial mechanics. We further examine the enlargement of aortic aneurysms for which we identify new mechanobiological insights into factors that affect the nearby non-aneurysmal segment as it responds to the changing mechanics within the diseased segment. Because this new 3D approach can be implemented within many existing finite element solvers, constrained mixture models of growth and remodeling can now be used more widely.

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.

Modeling lamellar disruption within the aortic wall using a particle-based approach

Aortic dissections associate with medial degeneration, thus suggesting a need to understand better the biophysical interactions between the cells and matrix that constitute the middle layer of the aortic wall. Here, we use a recently extended "Smoothed Particle Hydrodynamics" formulation to examine potential mechanisms of aortic delamination arising from smooth muscle cell (SMC) dysfunction or apoptosis, degradation of or damage to elastic fibers, and pooling of glycosaminoglycans (GAGs), with associated losses of medial collagen in the region of the GAGs. First, we develop a baseline multi-layered model for the healthy aorta that delineates medial elastic lamellae and intra-lamellar constituents. Next, we examine stress fields resulting from the disruption of individual elastic lamellae, lost SMC contractility, and GAG production within an intra-lamellar space, focusing on the radial transferal of loading rather than on stresses at the tip of the delaminated tissue. Results suggest that local disruptions of elastic lamellae transfer excessive loads to nearby intra-lamellar constituents, which increases cellular vulnerability to dysfunction or death. Similarly, lost SMC function and accumulations of GAGs increase mechanical stress on nearby elastic lamellae, thereby increasing the chance of disruption. Overall these results suggest a positive feedback loop between lamellar disruption and cellular dropout with GAG production and lost medial collagen that is more pronounced at higher distending pressures. Independent of the initiating event, this feedback loop can catastrophically propagate intramural delamination.

Growth and remodelling of living tissues: perspectives, challenges and opportunities

One of the most remarkable differences between classical engineering materials and living matter is the ability of the latter to grow and remodel in response to diverse stimuli. The mechanical behaviour of living matter is governed not only by an elastic or viscoelastic response to loading on short time scales up to several minutes, but also by often crucial growth and remodelling responses on time scales from hours to months. Phenomena of growth and remodelling play important roles, for example during morphogenesis in early life as well as in homeostasis and pathogenesis in adult tissues, which often adapt to changes in their chemo-mechanical environment as a result of ageing, diseases, injury or surgical intervention. Mechano-regulated growth and remodelling are observed in various soft tissues, ranging from tendons and arteries to the eye and brain, but also in bone, lower organisms and plants. Understanding and predicting growth and remodelling of living systems is one of the most important challenges in biomechanics and mechanobiology. This article reviews the current state of growth and remodelling as it applies primarily to soft tissues, and provides a perspective on critical challenges and future directions.

Mechanobiological Stability of Biological Soft Tissues

Like all other materials, biological soft tissues are subject to general laws of physics, including those governing mechanical equilibrium and stability. In addition, however, these tissues are able to respond actively to changes in their mechanical and chemical environment. There is, therefore, a pressing need to understand such processes theoretically. In this paper, we present a new rate-based constrained mixture formulation suitable for studying mechanobiological equilibrium and stability of soft tissues exposed to transient or sustained changes in material composition or applied loading. These concepts are illustrated for canonical problems in arterial mechanics, which distinguish possible stable versus unstable mechanobiological responses. Such analyses promise to yield insight into biological processes that govern both health and disease progression.

Critical roles of time-scales in soft tissue growth and remodeling

Most soft biological tissues exhibit a remarkable ability to adapt to sustained changes in mechanical loads. These macroscale adaptations, resulting from mechanobiological cellular responses, are important determinants of physiological behaviors and thus clinical outcomes. Given the complexity of such adaptations, computational models can significantly increase our understanding of how contributions of different cell types or matrix constituents, and their rates of turnover and evolving properties, ultimately change the geometry and biomechanical behavior at the tissue level. In this paper, we examine relative roles of the rates of tissue responses and external loading and present a new rate-independent approach for modeling the evolution of soft tissue growth and remodeling. For illustrative purposes, we also present numerical results for arterial adaptations. In particular, we show that, for problems defined by particular characteristic times, this approximate theory captures well the predictions of a fully general constrained mixture theory at a fraction of the computational cost.