A homogenized constrained mixture (and mechanical analog) model for growth and remodeling of soft tissue

Most mathematical models of the growth and remodeling of load-bearing soft tissues are based on one of two major approaches: a kinematic theory that specifies an evolution equation for the stress-free configuration of the tissue as a whole or a constrained mixture theory that specifies rates of mass production and removal of individual constituents within stressed configurations. The former is popular because of its conceptual simplicity, but relies largely on heuristic definitions of growth; the latter is based on biologically motivated micromechanical models, but suffers from higher computational costs due to the need to track all past configurations. In this paper, we present a temporally homogenized constrained mixture model that combines advantages of both classical approaches, namely a biologically motivated micromechanical foundation, a simple computational implementation, and low computational cost. As illustrative examples, we show that this approach describes well both cell-mediated remodeling of tissue equivalents in vitro and the growth and remodeling of aneurysms in vivo. We also show that this homogenized constrained mixture model suggests an intimate relationship between models of growth and remodeling and viscoelasticity. That is, important aspects of tissue adaptation can be understood in terms of a simple mechanical analog model, a Maxwell fluid (i.e., spring and dashpot in series) in parallel with a “motor element” that represents cell-mediated mechanoregulation of extracellular matrix. This analogy allows a simple implementation of homogenized constrained mixture models within commercially available simulation codes by exploiting available models of viscoelasticity.     

A 3-D constrained mixture model for mechanically mediated vascular growth and remodeling

In contrast to the widely applied approach to model soft tissue remodeling employing the concept of volumetric growth, microstructurally motivated models are capable of capturing many of the underlying mechanisms of growth and remodeling; i.e., the production, removal, and remodeling of individual constituents at different rates and to different extents. A 3-dimensional constrained mixture computational framework has been developed for vascular growth and remodeling, considering new, microstructurally motivated kinematics and constitutive equations and new stress and muscle activation mediated evolution equations. Our computational results for alterations in flow and pressure, using reasonable physiological values for rates of constituent growth and turnover, concur with findings in the literature. For example, for flow-induced remodeling, our simulations predict that, although the wall shear stress is restored completely, the circumferential stress is not restored employing realistic physiological rate parameters. Also, our simulations predict different levels of thickening on inner versus outer wall locations, as shown in numerous reports of pressure-induced remodeling. Whereas the simulations are meant to be illustrative, they serve to highlight the experimental data currently lacking to fully quantify mechanically mediated adaptations in the vasculature.     

Theoretical study on the effects of pressure-induced remodeling on geometry and mechanical non-homogeneity of conduit arteries

A structure-based mathematical model for the remodeling of arteries in response to sustained hypertension is proposed. The model is based on the concepts of volumetric growth and constitutive modeling of the arterial tissue within the framework of the constrained mixture theory. The major novel result of this study is that remodeling is associated with a local change in the mass fractions of the wall constituents that ultimately leads to mechanical non-homogeneity of the arterial wall. In the new homeostatic state that develops after a sustained increase in arterial pressure, the mass fraction of elastin decreases from the intimal side to the adventitial side of arteries, while the collagen fraction manifests an opposite trend. The results obtained are supported by some experimental observations reported in the literature.     

A constituent-based model of age-related changes in conduit arteries

In the present report, a constituent-based theoretical model of age-related changes in geometry and mechanical properties of conduit arteries is proposed. The model was based on the premise that given the time course of the load on an artery and the accumulation of advanced glycation end-products in the arterial tissue, the initial geometric dimensions and properties of the arterial tissue can be predicted by a solution of a boundary value problem for the governing equations that follow from finite elasticity, structure-based constitutive modeling within the constrained mixture theory, continuum damage theory, and global growth approach for stress-induced structure-based remodeling. An illustrative example of the age-related changes in geometry, structure, composition, and mechanical properties of a human thoracic aorta is considered. Model predictions were in good qualitative agreement with available experimental data in the literature. Limitations and perspectives for refining the model are discussed.     

A growth mixture theory for cartilage with application to growth-related experiments on cartilage explants

In this paper, we present a growth mixture model for cartilage. The main features of this model are illustrated in a simple equilibrium boundary-value problem that is chosen to illustrate how a mechanical theory of cartilage growth may be applied to growth-related experiments on cartilage explants. The cartilage growth mixture model describes the independent growth of the proteoglycan and collagen constituents due to volumetric mass deposition, which leads to the remodeling of the composition and the mechanical properties of the solid matrix. The model developed here also describes how the material constants of the collagen constituent depend on a scalar parameter that may change over time (e.g., crosslink density); this leads to a remodeling of the structural and mechanical properties of the collagen constituent. The equilibrium boundary-value problem that describes the changes observed in cartilage explants harvested at different stages of a growth or a degenerative process is formulated. This boundary-value problem is solved using existing experimental data for developing bovine cartilage explants harvested at three developmental stages. The solution of the boundary-value problem in conjunction with existing experimental data suggest the types of experimental studies that need to be conducted in the future to determine model parameters and to further refine the model.     

Hermitian splines for modeling biological soft tissue systems that exhibit nonlinear force-elongation curves

Numerical simulation of soft tissue mechanical properties is a critical step in developing valuable biomechanical models of live organisms. A cubic Hermitian spline optimization routine is proposed in this paper to model nonlinear experimental force-elongation curves of soft tissues, in particular when modeled as lumped elements. Boundary conditions are introduced to account for the positive definiteness and the particular curvature of the experimental curve to be fitted. The constrained least-square routine minimizes user intervention and optimizes fitting of the experimental data across the whole fitting range. The routine provides coefficients of a Hermitian spline or corresponding knots that are compatible with a number of constraints that are suitable for modeling soft tissue tensile curves. These coefficients or knots may become inputs to user-defined component properties of various modeling software. Splines are particularly advantageous over the well-known exponential model to account for the traction curve flatness at low elongations and to allow for more flexibility in the fitting process. This is desirable as soft tissue models begin to include more complex physical phenomena.     

Computational modeling of volumetric soft tissue growth: application to the cardiac left ventricle

As an initial step to investigate stimulus–response relations in growth and remodeling (G&R) of cardiac tissue, this study aims to develop a method to simulate 3D-inhomogeneous volumetric growth. Growth is regarded as a deformation that is decomposed into a plastic component which describes unconstrained growth and an elastic component to satisfy continuity of the tissue after growth. In current growth models, a single reference configuration is used that remains fixed throughout the entire growth process. However, considering continuous turnover to occur together with growth, such a fixed reference is unlikely to exist in reality. Therefore, we investigated the effect of tissue turnover on growth by incrementally updating the reference configuration. With both a fixed reference and an updated reference, strain-induced cardiac growth in magnitude of 30% could be simulated. However, with an updated reference, the amplitude of the stimulus for growth decreased over time, whereas with a fixed reference this amplitude increased. We conclude that, when modeling volumetric growth, the choice of the reference configuration is of great importance for the computed growth.     

Locally constrained mixture representation of dynamic imaging data from PET and MR studies

Dynamic positron emission tomography (PET) studies provide measurements of the kinetics of radiotracers in living tissue. This is a powerful technology which can play a major role in the study of biological processes, potentially leading to better understanding and treatment of disease. Dynamic PET data relate to complex spatiotemporal processes and its analysis poses significant challenges. In previous work, mixture models that expressed voxel-level PET time course data as a convex linear combination of a finite number of dominant time course characteristics (called sub-TACs) were introduced. This paper extends that mixture model formulation to allow for a weighted combination of scaled sub-TACs and also considers the imposition of local constraints in the number of sub-TACs that can be active at any one voxel. An adaptive 3D scaled segmentation algorithm is developed for model initialization. Increases in the weighted residual sums of squares is used to guide the choice of the number of segments and the number of sub-TACs in the final mixture model. The methodology is applied to five data sets from representative PET imaging studies. The methods are also applicable to other contexts in which dynamic image data are acquired. To illustrate this, data from an echo-planar magnetic resonance (MR) study of cerebral hemodynamics are considered. Our analysis shows little indication of departure from a locally constrained mixture model representation with at most two active components at any voxel. Thus, the primary sources of spatiotemporal variation in representative dynamic PET and MR imaging studies would appear to be accessible to a substantially simplified representation in terms of the generalized locally constrained mixture model introduced.     

A goal function approach to remodeling of arteries uncovers mechanisms for growth instability

A novel, goal function-based formulation for the growth dynamics of arteries is introduced and used for investigating the development of growth instability in blood vessels. Such instabilities would lead to abnormal growth of the vessel, reminiscent of an aneurysm. The blood vessel is modeled as a thin-walled cylindrical tube, and the constituents that form the vessel wall are assumed to deform together as a constrained mixture. The growth dynamics of the composite material of the vessel wall are described by an evolution equation, where the effective area of each constituent changes in the direction of steepest descent of a goal function. This goal function is formulated in such way that the constituents grow toward a target potential energy and a target composition. The convergence of the simulated response of the evolution equation toward a target homeostatic state is investigated for a range of isotropic and orthotropic material models. These simulations suggest that elastin-deficient vessels are more prone to growth instability. Increased stiffness of the vessel wall, on the other hand, gives a more stable growth process. Another important finding is that an increased rate of degradation of materials impairs growth stability.     

A constrained mixture model for developing mouse aorta

Mechanical stresses influence the structure and function of adult and developing blood vessels. When these stresses are perturbed, the vessel wall remodels to return the stresses to homeostatic levels. Constrained mixture models have been used to predict remodeling of adult vessels in response to step changes in blood pressure, axial length and blood flow, but have not yet been applied to developing vessels. Models of developing blood vessels are complicated by continuous and simultaneous changes in the mechanical forces. Understanding developmental growth and remodeling is important for treating human diseases and designing tissue-engineered blood vessels. This study presents a constrained mixture model for postnatal development of mouse aorta with multiple step increases in pressure, length and flow. The baseline model assumes that smooth muscle cells (SMCs) in the vessel wall immediately constrict or dilate the inner radius after a perturbation to maintain the shear stress and then remodel the wall thickness to maintain the circumferential stress. The elastin, collagen and SMCs have homeostatic stretch ratios and passive material constants that do not change with developmental age. The baseline model does not predict previously published experimental data. To approximate the experimental data, it must be assumed that the SMCs dilate a constant amount, regardless of the step change in mechanical forces. It must also be assumed that the homeostatic stretch ratios and passive material constants change with age. With these alterations, the model approximates experimental data on the mechanical properties and dimensions of aorta from 3- to 30-day-old mice.     

Finite element implementation of a new model of slight compressibility for transversely isotropic materials

Modelling transversely isotropic materials in finite strain problems is a complex task in biomechanics, and is usually addressed by using finite element (FE) simulations. The standard method developed to account for the quasi-incompressible nature of soft tissues is to decompose the strain energy function (SEF) into volumetric and deviatoric parts. However, this decomposition is only valid for fully incompressible materials, and its use for slightly compressible materials yields an unphysical response during the simulation of hydrostatic tension/compression of a transversely isotropic material. This paper presents the FE implementation as subroutines of a new volumetric model solving this deficiency in two FE codes: Abaqus and FEBio. This model also has the specificity of restoring the compatibility with small strain theory. The stress and elasticity tensors are first derived for a general SEF. This is followed by a successful convergence check using a particular SEF and a suite of single-element tests showing that this new model does not only correct the hydrostatic deficiency but may also affect stresses during shear tests (Poynting effect) and lateral stretches during uniaxial tests (Poisson's effect). These FE subroutines have numerous applications including the modelling of tendons, ligaments, heart tissue, etc. The biomechanics community should be aware of specificities of the standard model, and the new model should be used when accurate FE results are desired in the case of compressible materials.     

Collagen fiber alignment does not explain mechanical anisotropy in fibroblast populated collagen gels

Many load-bearing soft tissues exhibit mechanical anisotropy. In order to understand the behavior of natural tissues and to create tissue engineered replacements, quantitative relationships must be developed between the tissue structures and their mechanical behavior. We used a novel collagen gel system to test the hypothesis that collagen fiber alignment is the primary mechanism for the mechanical anisotropy we have reported in structurally anisotropic gels. Loading constraints applied during culture were used to control the structural organization of the collagen fibers of fibroblast populated collagen gels. Gels constrained uniaxially during culture developed fiber alignment and a high degree of mechanical anisotropy, while gels constrained biaxially remained isotropic with randomly distributed collagen fibers. We hypothesized that the mechanical anisotropy that developed in these gels was due primarily to collagen fiber orientation. We tested this hypothesis using two mathematical models that incorporated measured collagen fiber orientations: a structural continuum model that assumes affine fiber kinematics and a network model that allows for nonaffine fiber kinematics. Collagen fiber mechanical properties were determined by fitting biaxial mechanical test data from isotropic collagen gels. The fiber properties of each isotropic gel were then used to predict the biaxial mechanical behavior of paired anisotropic gels. Both models accurately described the isotropic collagen gel behavior. However, the structural continuum model dramatically underestimated the level of mechanical anisotropy in aligned collagen gels despite incorporation of measured fiber orientations; when estimated remodeling-induced changes in collagen fiber length were included, the continuum model slightly overestimated mechanical anisotropy. The network model provided the closest match to experimental data from aligned collagen gels, but still did not fully explain the observed mechanics. Two different modeling approaches showed that the level of collagen fiber alignment in our uniaxially constrained gels cannot explain the high degree of mechanical anisotropy observed in these gels. Our modeling results suggest that remodeling-induced redistribution of collagen fiber lengths, nonaffine fiber kinematics, or some combination of these effects must also be considered in order to explain the dramatic mechanical anisotropy observed in this collagen gel model system.