Modeling biological growth and remodeling: Contrasting methods, contrasting needs

Bayesian calibration of a computational model of tissue expansion based on a porcine animal model

Tissue expansion is a technique used clinically to grow skin in situ to correct large defects. Despite its enormous potential, lack of fundamental knowledge of skin adaptation to mechanical cues, and lack of predictive computational models limit the broader adoption and efficacy of tissue expansion. In our previous work, we introduced a finite element model of tissue expansion that predicted key patterns of strain and growth which were then confirmed by our porcine animal model. Here we use the data from a new set of experiments to calibrate the computational model within a Bayesian framework. Four 10×10cm2 patches were tattooed in the dorsal skin of four 12 weeks-old minipigs and a total of six patches underwent successful tissue expander placement and inflation to 60cc for expansion times ranging from 1 h to 7 days. Six patches that did not have expanders implanted served as controls for the analysis. We find that growth can be explained based on the elastic deformation. The predicted area growth rate is k∈[0.02,0.08] [h-1]. Growth is anisotropic and reflects the anisotropic mechanical behavior of porcine dorsal skin. The rostral-caudal axis shows greater deformation than the transverse axis, and the time scale of growth in the rostral-caudal direction is given by rate parameters k1∈[0.04,0.1] [h-1] compared to k2∈[0.01,0.05] [h-1] in the transverse direction. Moreover, the calibration results underscore the high variability in biological systems, and the need to create probabilistic computational models to predict tissue adaptation in realistic settings. STATEMENT OF SIGNIFICANCE: Tissue expansion is a widely used technique in reconstructive surgery because it triggers growth of skin for the correction of large skin lesions and for breast reconstruction after mastectomy. Despite of its potential, complications and undesired outcomes persist due to our incomplete understanding of skin mechanobiology. Here we quantify the deformation and growth fields induced by an expander over 7 days in a porcine animal model and use these data to calibrate a computational model of skin growth using finite element simulations and a Bayesian framework. The calibrated model is a leap forward in our understanding skin growth, we now have quantitative understanding of this process: area growth is anisotropic and it is proportional to stretch with a characteristic rate constant of k∈[0.02,0.08] [h-1].

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.

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.