Downloaded from rsif.royalsocietypublishing.org on June 15, 2014
Computational modelling suggests good, bad and ugly roles of glycosaminoglycans in arterial wall mechanics and mechanobiology S. Roccabianca, C. Bellini and J. D. Humphrey J. R. Soc. Interface 2014 11, 20140397, published 11 June 2014
References
This article cites 61 articles, 18 of which can be accessed free
Email alerting service
Receive free email alerts when new articles cite this article - sign up in the box at the top right-hand corner of the article or click here
http://rsif.royalsocietypublishing.org/content/11/97/20140397.full.html#ref-list-1
To subscribe to J. R. Soc. Interface go to: http://rsif.royalsocietypublishing.org/subscriptions
Downloaded from rsif.royalsocietypublishing.org on June 15, 2014
rsif.royalsocietypublishing.org
Computational modelling suggests good, bad and ugly roles of glycosaminoglycans in arterial wall mechanics and mechanobiology S. Roccabianca1, C. Bellini1 and J. D. Humphrey1,2 1
Research Cite this article: Roccabianca S, Bellini C, Humphrey JD. 2014 Computational modelling suggests good, bad and ugly roles of glycosaminoglycans in arterial wall mechanics and mechanobiology. J. R. Soc. Interface 11: 20140397. http://dx.doi.org/10.1098/rsif.2014.0397
Received: 14 April 2014 Accepted: 21 May 2014
Subject Areas: biomechanics
2
Department of Biomedical Engineering, Yale University, New Haven, CT, USA Vascular Biology and Therapeutics Program, Yale School of Medicine, New Haven, CT, USA The medial layer of large arteries contains aggregates of the glycosaminoglycan hyaluronan and the proteoglycan versican. It is increasingly thought that these aggregates play important mechanical and mechanobiological roles despite constituting only a small fraction of the normal arterial wall. In this paper, we offer a new hypothesis that normal aggregates of hyaluronan and versican pressurize the intralamellar spaces, and thereby put into tension the radial elastic fibres that connect the smooth muscle cells to the elastic laminae, which would facilitate mechanosensing. This hypothesis is supported by novel computational simulations using two complementary models, a mechanistically based finite-element mixture model and a phenomenologically motivated continuum hyperelastic model. That is, the simulations suggest that normal aggregates of glycosaminoglycans/proteoglycans within the arterial media may play equally important roles in supporting (i.e. a structural role) and sensing (i.e. an instructional role) mechanical loads. Additional simulations suggest further, however, that abnormal increases in these aggregates, either distributed or localized, may over-pressurize the intralamellar units. We submit that these situations could lead to compromised mechanosensing, anoikis and/or reduced structural integrity, each of which represent fundamental aspects of arterial pathologies seen, for example, in hypertension, ageing and thoracic aortic aneurysms and dissections.
Keywords: hyaluronan, versican, mechanosensing, anoikis, stress analysis
1. Introduction Author for correspondence: J. D. Humphrey e-mail:
[email protected] Glycosaminoglycans (GAGs) compose only 2 –5%, by dry weight, of the wall of normal elastic arteries [1], which may be one reason why their potential mechanical roles have not been well studied in arterial wall mechanics. Nevertheless, increasing evidence shows that GAGs are important contributors to normal arterial structure and function [2–4], and that GAG content tends to increase in hypertension, ageing, atherosclerosis and thoracic aortic aneurysms and dissection [5 –8]. That is, GAGs now appear to be important contributors to arterial morphogenesis, homeostasis and pathogenesis. The three primary classes of GAGs within the arterial wall are hyaluronan, chondroitin/dermatan sulfate and heparan sulfate; hyaluronan and the chondroitin/dermatan sulfate GAGs are synthesized mainly by smooth muscle cells, whereas the heparan sulfate GAGs are synthesized mainly by endothelial cells [9]. Hyaluronan does not associate with a protein core, but chondroitin, dermatan and heparan sulfate GAGs bind to protein cores to form larger molecules called proteoglycans (PGs). The two primary classes of PGs of importance in arterial mechanics and biology are those that aggregate with hyaluronan and contribute to the compressive stiffness and viscoelasticity of the wall (e.g. the chondroitin sulfate-rich PG called versican) and those that participate in collagen fibrillogenesis and the sequestering of growth factors (e.g. the small leucinerich PGs called biglycan and decorin). Herein, however, we are motivated by the potential mechanical roles of aggregates of hyaluronan and versican,
& 2014 The Author(s) Published by the Royal Society. All rights reserved.
Downloaded from rsif.royalsocietypublishing.org on June 15, 2014
Our hypothesis should hold for elastic arteries in general, but owing to the availability of copious data on the murine
3. Finite-element model 3.1. General framework Ateshian et al. [27] developed a mixture-based finite-element approach to describe nonlinear responses of GAG-endowed continua; this formulation was implemented in the open source software FEBio (www.sci.utah.edu/software/40mrl/39-febio). They demonstrated the utility of this approach, for example, by studying residual stresses in rodent aortic tissue having layer-specific differences in fixed charge density (FCD; i.e. GAG-associated negative charges), as described in [11]. Herein, we use their approach as implemented in FEBio. Figure 1 shows half of our baseline axisymmetric finite-element model of an intact artery to emphasize the need to model possible transmural differences in the FCD, that is, the distribution of GAGs. Also shown is an illustrative computational result for the radial cutting of an intact traction-free configuration. A series of such computations of the resulting opening angle F0 served to validate our use of FEBio via a comparison of results with those in [11]. That is, although our focus herein is on the in vivo, not residual, stress field, calculation of residual stress-related quantities proved useful in comparing our two models. As shown below, we also used FEBio to model a small portion (intramural tangent plane) of the arterial wall that contained a focal, rather than diffuse, accumulation of GAGs. Convergence of solutions in both cases was ensured via standard mesh refinement and comparison of results.
2
J. R. Soc. Interface 11: 20140397
2. Illustrative artery of interest
common carotid artery in our laboratory (cf. [19]), we illustrate our ideas based on simulations for this vessel. The normal mouse carotid artery consists of two structurally important layers, the media and adventitia, which occupy approximately 40% and 60% of the wall, respectively. The media consists of two to three medial lamellar units, which in turn consist of smooth muscle, collagen and GAGs/PGs bounded by paired elastic laminae; the adventitia consists of a single layer of fibrillar collagen, with small amounts of elastic fibres. Of particular importance herein, there also exist elastin-associated microfibrils that either connect the smooth muscle cells to the elastic laminae or connect the elastic laminae to each other; these microfibrils are so thin that they are best seen via electron microscopy [20,21]. Reports differ on the overall mass fractions of these primary constituents, perhaps depending on the amount of perivascular tissue/adventitia that is removed during dissection. We use the report that the murine common carotid artery consists of approximately 26% elastic fibres, 21% smooth muscle and 50% collagen, with the remaining 3% mainly medial GAGs/PGs [22]. When transected, these vessels tend to retract approximately 70%, which is thought to be due primarily to recoil of the elastic fibres [23], which are deposited and cross-linked during the perinatal period [4] and subsequently extended elastically during somatic growth [24]. When a traction-free ring is cut radially, a small opening angle Fo occurs owing to the relief of residual stresses (cf. [25]). We note that medial aggregates of hyaluronan –versican would cause an intralamellar swelling which, along with radial gradients in material properties [12] and the remaining recoil of elastic fibres [26], probably affects the measured opening angle (cf. [11]) and hence the overall distribution of wall stress.
rsif.royalsocietypublishing.org
which are found primarily within the extracellular space of the media of elastic arteries [10]. A recent study suggested, for example, that medial GAGs/PGs may contribute to the existence of residual stress in arteries [11], which of course play important roles in homeostasis as well as adaptations to altered haemodynamics [12]. Hyaluronan is a high-molecular-weight biomolecule that sequesters considerable water, and thereby contributes to the compressive stiffness and viscoelastic behaviour of many soft tissues [13]. Versican contains 15–20 chondroitin sulfate side chains, which are highly negatively charged [1]. Hence, versican tends to attract positive ions (e.g. Naþ) to enforce local electroneutrality, which in turn imbibes water and generates an associated local swelling pressure. This so-called Donnan swelling pressure is well known in physical chemistry [14] and is well studied in cartilage mechanics wherein PGs play a fundamental structural role [15]. Of particular relevance herein, research on cartilage revealed that a local swelling pressure can cause collagen II fibres to be in tension despite the existence of an overall compressive stress [16]. Motivated by this finding, our focus is on potential effects of intralamellar aggregates of hyaluronan and versican on the elastic fibre-associated microfibrillar connections between smooth muscle cells and the delimiting elastic laminae in large arteries. In particular, we hypothesize that, in addition to contributing to the compressive load carrying capability of the extracellular matrix, normal intralamellar aggregates of hyaluronan and versican can place radially oriented microfibrillar connections in tension despite the overall compressive radial stress in the arterial wall. This tension, in turn, could provide a means for the smooth muscle cells to sense directly the mechanical loads in the elastic fibre network that is fundamental to normal large artery structure and function. Pathological accumulations of these GAG/PG aggregates could overload these connections, however, thus resulting in compromised or lost mechanosensing. Indeed, in the extreme cases of localized pooling of these aggregates, a complete loss of fibre –cell connections may result in smooth muscle anoikis [17] or cause high-stress concentrations that could adversely affect the structural integrity of the wall and lead to delaminations [18]. In this paper, we use two complementary biomechanical modelling approaches—a mechanistic mixture-based finite-element model, and a phenomenological continuumbased semi-analytic model—to illustrate these potentially good, bad and ugly roles of medial GAGs/PGs in large arteries. Although contributions to normal load bearing and mechanosensing are obviously good, abnormal increases in GAGs/PGs that compromise or sever elastic fibre– smooth muscle connections would be bad, for this would adversely affect mechanosensing or cell survival in ways that further diminish the structural integrity of the wall. Finally, following up on a recent hypothesis [8], we re-emphasize that a pathological local pooling of GAGs/ PGs could lead to the ugly situation of local stress concentrations that may initiate potentially lethal intramural delaminations/dissections.
Downloaded from rsif.royalsocietypublishing.org on June 15, 2014
constitutive parameter:
fixed- charge density: FCDAdv
media
FCDMed1 FCDMed2
j
W kj ¼
F*0 /2 g2
Figure 1. Finite-element model of a segment of a murine common carotid artery, where g2 represents a possible ‘good’ situation with FCDMed2 ¼ 270, FCDMed1 ¼ 240 and FCDAdv ¼ 220 mEq l21. Mechanical properties (i.e. strain energy parameters) were prescribed separately for the media and adventitia and similarly for the swelling related fixed charge density (FCD), with the exception of radial differences in FCD also within the media.
3.2. Constitutive behaviour Consistent with the vast majority of the literature on arterial wall mechanics [12], we used a locally homogenized stored energy function W(C) to describe the material behaviour, albeit separately for the media and adventitia. Here, C ¼ FT F, where C is the right Cauchy– Green tensor, and F is the deformation gradient tensor. To facilitate comparisons between the finite-element and continuum (described below) models, we used a strain energy function W in the finite-element simulations having both neo-Hookean and exponential (four fibre families) contributions. This socalled four fibre family constitutive relation has been shown to describe well the mean biaxial mechanical data of central arteries in mice [19,22]. Relative to the mass-averaged strain energy function that we used for the continuum model (see below), the finite-element model (FEM) accounted for an amorphous matrix consisting of elastin plus GAGs/PGs and modelled as a single material; it also combined the circumferential behaviours of the collagen fibres and smooth muscle. Hence, the overall relation is matrix W FEM ¼ fmatrix þ M,A W
4 X
k
j fM,A W kj ,
j¼1
where fmatrix and fkj are mass fractions that can be prescribed separately in the media (M) and adventitia (A). As is usual (cf. [28]), we let W matrix ¼
cmatrix (IC 3), 2
j 4c2
j
j
2
{exp[c2 (IVC 1) ] 1},
j
where c1 and c2 are material parameters (the former with dimension of stress and the latter non-dimensional) and j IVC ¼ Mj CMj represent squares of the stretch of the jth fibre family. The definition of fibre orientations was consistent with that for the continuum model: family one in the axial direction, family two in the circumferential direction and families three and four representing symmetric-diagonal orientations. It is important to note that we allowed the three primary families of fibres (i.e. axial, circumferential and symmetric-diagonal) to have different material behaviours, hence necessitating estimation of up to seven fibre-related material parameters (see table 1). This material characterization differed between the two models as noted below.
4. Continuum model Although one should prefer a mixture-based finite-element model that directly includes effects of FCD on the local Donnan swelling pressure, information on spatial and temporal changes in the FCD in arteries remains wanting (cf. [18]). Hence, phenomenological continuum models can serve as useful tools for studying parametrically the consequences of changes in Donnan swelling pressures, that is, imposed local changes in volume (detF = 1). Although swelling causes a change in volume, the swollen artery appears to deform isochorically thereafter in response to transient loading; hence, one can use incompressible finite elasticity provided the swelling is included. Towards this end, we used but extended the method of Demirkoparan & Pence [29], who analysed effects of swelling on fibre-reinforced cylindrical tubes with v ¼ detF a constant function of radius. We first studied the effect of a constant swelling on a bilayered arterial model, but then studied effects of more physiological situations wherein the media swelled more than the adventitia. This transmural difference was quantified using a sigmoidal function 100(@MA @) n ¼ na þ qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi , 1 þ [(200=(ni þ na ))(@MA @)]2 where ni and na represent values of swelling in inner (medial side) and outer (advential side) portions of the wall, respectively, and @MA represents the location of the medial– adventitial interface. In this way, we were able to compare, for the first time, the predictions of a mechanistic finite-element model and a phenomenological continuum model.
4.1. Configurations of interest Figure 2 is a schema of general configurations of interest computationally; others can be obtained easily from the ones shown. We again modelled the normal arterial wall as a bilayered structure, noting that the inner (i.e. medial) layer
J. R. Soc. Interface 11: 20140397
j
c1
3
rsif.royalsocietypublishing.org
adventitia
where cmatrix is a material parameter and IC ¼ trC. Note, therefore, that effects of ‘swelling’ associated with the GAGs/PGs were accounted for in the finite-element model by an additional contribution to the tangent stiffness matrix that depends on detF and the FCD (cf. equation 9 in [11]). Similarly, we described the behaviour of the four fibre families (dominated by collagen) using a Fung exponential strain energy function, namely
Downloaded from rsif.royalsocietypublishing.org on June 15, 2014
rP ¼ rp (r), uP ¼ u, zP ¼ z
(from kh to kP ), (r, u, z) ! (@, q, z): @ ¼ @(r), q ¼ u, z ¼
z liv z
(from kh to ktf ), (@, q, z) ! (R, Q, Z): R ¼ R(@), Q ¼
p F0 z q,Z¼ LZ p
(from ktf to kc ) (@, q, z) ! (@ , q , z ):
and
@ ¼ @ (@) , q ¼ q, z ¼ lz z
(from ktf to ktf ): Similarly, consider mappings of the swollen wall, (@ , q , z ) ! (r , u , z ): (from ktf
r ¼ r (@ ), u ¼ q , z ¼ liv z z
to kh ),
(@ , q , z ) ! (R , Q , Z ): R ¼ R (@ ), Q ¼ and
p F0 z q, Z ¼ p LZ
(from ktf to kc )
(r , u , z ) ! (rP , uP , zP ): rP ¼ rP (r ), uP ¼ u , zP ¼ z
(from kh to kP ), where F0 and F0 are the opening angle in the cut, stress-free configuration and the opening angle in the swollen, cut, stress-free configuration, respectively. The deformation
while Fh assumes the same form as F1 tf when r becomes r , @ iv iv becomes @ and lz becomes lz . Similarly, Fc assumes the same form as Fc when R becomes R , @ becomes @ , F0 becomes F0 and Lz becomes LZ , and FP assumes the same form of FP when r becomes r and rP becomes rP . The associated deformation gradients of interest for stress analysis were obtained via multiplicative operations. For example, relative to kh, the deformation gradient associated with kh was F ¼ Fh F Ftf . Of course, this general schema allowed different degrees of swelling (parametrized by det(F ) ¼ n ) as well as different levels of biaxial loading via pressurization and applied axial loads. In general, the value of the relative volume change (equals to either 1 or n ) allowed the radial stretch to be computed, in each deformation process, from the circumferential and axial stretches and the radial position occupied by each material particle that was mapped between different configurations, namely qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi rP ¼ r2Pi þ (r2 r2i ) (from kh to kP ), qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi rP ¼ (rPi )2 þ ((r )2 (ri )2 ) (from kh to kP ), qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 2 2 @ ¼ @2i þ liv z (r ri ) (from kh to ktf ), sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 1 r ¼ (ri )2 þ iv ((@ )2 (@i )2 ) (from ktf to kh ), lz rffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi p R ¼ R2i þ Lz (@2 @2i ) (from ktf to kc ), p F0 rffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi p R ¼ (Ri )2 þ L ((@ )2 (@i )2 ) (from ktf to kc ), p F0 z vffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi u ð@ u 2 u 2 and @ ¼ t(@i ) þ n (@)@ d@ (from ktf to ktf ), lz
@i
where subscripts i denote inner radii.
4.2. Constitutive relation Consistent with the kinematics, our constitutive approach borrowed ideas from the growth and remodelling literature, which allowed us to endow each constituent with not only a different stored energy function, but also a different strain in the homeostatic configuration. Briefly, the elastindominated behaviour was described by a neo-Hookean relation, whereas both the collagen-dominated and the passive smooth muscle behaviours were described by locally parallel fibre families described by Fung exponentials [26]. In addition, however, we added a new term (modified Blatz –Ko material) to account for the GAG/PG-dominated behaviour and related swelling. Although each constituent was constrained to deform with the wall (i.e. the mixture), its individual natural configuration allowed different levels of deformation within any overall configuration. Recalling
4
J. R. Soc. Interface 11: 20140397
(r, u, z) ! (rP , uP , zP ):
gradients F associated with these mappings, as shown in figure 2, are computed via @ rP rP [FP ] ¼ diag 1 , @r r " # @@ @ 1 [Ftf ] ¼ diag , @ r r liv z @ R R p F0 1 [Fc ] ¼ diag @@ @ p Lz @ @ @ lz , and [F ] ¼ diag n (@) @@ @
rsif.royalsocietypublishing.org
is normally swollen slightly relative to the outer (adventitial) layer owing to the presence of hyaluronan and versican in the former [10]. The homeostatic configuration kh (with coordinates r, u, z) at mean arterial pressure (MAP) 93 mmHg and the in vivo axial stretch (liv z 1:7) defined a computationally convenient and biologically relevant reference configuration; note that this choice was different from the classical use of a traction-free reference that was used in the finite-element model, but it proved useful in the continuum model and is more relevant to the physiology [26]. Towards this end, we defined a natural, stress-free configuration kan for each structurally significant constituent a (¼ e, kj, m, PG for elastic fibres, j ¼ 1, 2, 3, 4 families of collagen fibres, smooth muscle cells and GAGs/PGs, respectively). Mappings from these natural configurations to the in vivo homeostatic reference configuration were defined by Ga, the constituent-specific deposition stretch tensors. From the central reference configuration, one can obtain any nonhomeostatic configuration kP at pressure P and liv z (with coordinates rP, uP, zP), the intact, traction-free unloaded configuration ktf (with coordinates @, q, z), and the radially cut, traction-free configuration kc (with coordinates R, Q, Z). To examine consequences of medial swelling, we also let the intact configuration ktf be mapped into a swollen, intact, traction-free configuration ktf (with coordinates @ , q , z ). From this swollen, traction-free configuration ktf one can axially stretch and pressurize the carotid to MAP to obtain a new swollen and physiologically loaded configuration kh (with coordinates r*, u*, z*) or pressurize the vessel to any perturbed pressure P to obtain kP (with coordinates rP , uP , zP ). Of course, as in the finite-element model, one can also introduce a radial cut numerically to obtain a nearly stress-free swollen configuration kc (with coordinates R , Q , Z ), which as noted above was useful in comparing results across models. The mappings of material points among the different configurations shown in figure 2 are given by
Downloaded from rsif.royalsocietypublishing.org on June 15, 2014
SMCs collagen Gm
Gk
rsif.royalsocietypublishing.org
elastin
5
GAGs
kn GPG
Ge
kh
Ftf Fc k220
ktf
F0 kc
F*
F* 0
k *c
k tf*
k 220 *
F*c F*h k *h F*220
Figure 2. Schematic diagram of key configurations that were used either in the nonlinear regression or to test the predictive ability of the constitutive formulations for the mouse carotid artery. kh (kh ), homeostatic (swollen homeostatic) configuration at a mean arterial pressure of 93 mmHg and in vivo axial stretch iv iv livz ¼ 1:70 (liv z ¼ 1:32 to 1:59). kp (kp ), loaded (swollen loaded) configurations at lz (lz ) and any pressure P, with representative values shown for P ¼ 220 mmHg. ktf (ktf ), intact (swollen intact), traction-free configuration. kc (kc ), radially cut (swollen radially cut), traction-free configuration. kn, natural (i.e. stress-free) configurations for individual constituents. that mappings between each natural configuration kan and the homeostatic reference configuration kh were described by the linear transformations Ga, the constituent-specific deformations needed in the stored energy functions were Fa ¼ FGa , where F represents any deformation of the arterial wall relative to the homeostatic reference configuration, as described previously. Hence, F 5 I in the homeostatic (loaded) configuration whereby Fa ¼ Ga, by definition. Moreover, note that Ga account for both the cell-mediated deposition stretch (e.g. when collagen fibres are deposited at a preferred stress in extant matrix) and developmentally induced stretches (e.g. in the elastic fibres, which are deposited during the perinatal period and have a long half-life). The mass-averaged strain energy function for the continuum model (CM) was thus assumed of the form PG PG m m W CM ¼ fPG (F ) þ feM,A W e (Fe ) þ fm M,A W M,A W (l )
þ
4 X j¼1
k
j fM,A W kj (lkj ),
where superscripts PG, e, m and kj refer to GAGs/PGs, elastin, smooth muscle and each of the four families of collagen fibres ( j ¼ 1, 2, 3, 4), and Wa and faM,A are the strain energy density functions and the mass fractions of the a constituents that compose the mixture, whereas a subscript M or A again refers separately to the medial or adventitial layer, respectively. Similar to that in the finite-element model, we modelled the elastin-dominated behaviour using a neo-Hookean model, W e (Fe ) ¼
ce (trCe 3), 2
where c e is a material parameter with the dimension of stress, but C e is now the constituent-specific right Cauchy – Green tensor for elastin, where Ce ¼ (Fe )T Fe and Fe ¼ FGe . Also, similar to the finite-element model, we let the behaviour of the fibrillar collagen and passive smooth muscle in tension be described by Fung exponential strain
J. R. Soc. Interface 11: 20140397
F220
0.14
a0k (deg) c3,4 2
2.19
2.10
(kPa)
2.58 1024
1.09 1024
0.05
2.27 1022
28.25
27.75 28.01 0.05 17.64 one layer
81.47
80.88 0.05 142.10
(kPa) (kPa) (kPa)
19.51
c3,4 1
36.67 19.43
c22 c21
449.44 1.47
c12 c11
1275.97
bz bu ak0 ( deg ) c2k c1k (kPa)
0.11
two layers
where N is the total number of data points and the overbars denote mean values. Goodness of fit, evaluated in terms of the root mean sum-of-the-squares of the error, was always between 0.02 and 0.15. The resulting best-fit parameters as well as values of the mass fractions and distributions of constituents in each layer determined from histology (cf. [26]) are in tables 1 and 2. Note that the parameters bu and bz, which
c
Biaxial data from three different cyclic pressure–diameter protocols at multiple fixed axial stretches [19] were used to find best-fit material parameters for the bilayered constitutive relations used in both the finite-element and continuum models. Specifically, best-fit values were determined by minimizing the error between theoretically predicted (th) and experimentally measured (exp) pressures and axial forces via a nonlinear least-squares regression, namely by minimizing the objective function 2 !2 3 th N exp 2 th exp X P P f f 4 5, e¼ þ exp f exp P i i¼1 i
FEM model
4.3. Parameter estimation
118.10
where J2PG is a function of the second and third invariant of CPG , namely J2PG ¼ IICPG =IIICPG , where IICPG ¼ 1=2[(trCPG )2 tr(CPG )2 ] and IIICPG ¼ detCPG ¼ (n )2 , cPG is a material parameter with the dimension of stress, CPG is the right Cauchy– Green tensor for the GAGs/PGs, evaluated as CPG ¼ (FPG )T FPG , where FPG ¼ FGPG describes the deformation gradient experienced by this amorphous matrix, and n* describes the relative change in volume with respect to the intact traction-free configuration ktf.
matrix
cPG PG ( J 3(n )2=3 ), 2 2
172.18
W PG (FPG ) ¼
J. R. Soc. Interface 11: 20140397
where ck1 and cm 1 are material parameters with dimension of stress, ck2 and cm 2 are non-dimensional material parameters for the collagen-dominated (superscript k) and smooth muscledominated (superscript m) behaviours, respectively. Neither the smooth muscle nor the collagen fibres were endowed with a compressive stiffness. Rather, the compressive behaviour of the mixture was modelled by the mechanical properties of the amorphous matrix. The deformation of each family of fibres pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi kj was evaluated as lkj ¼ Gpkjffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi Mkj CM ffi , for the jth collagen m m m fibre family, and l ¼ Gm M CM for the circumferentially oriented smooth muscle (i.e. Mm ¼ (0, 1, 0)). The collagen fibre orientations were defined by the unit vectors k k k Mkj ¼ (0, sin a0j , cos a0j ), where angles a0j are defined between directions of fibre reinforcement and the axial direction. Hence, axial and circumferential fibres are fixed at ak01 ¼ 0 ( deg) and ak02 ¼ 90 (deg), respectively, whereas symmetric-diagonal fibres are accounted for by ak03 ¼ ak04 ¼ ak0 , which was determined via parameter estimation using in vitro biaxial data. Finally, we modelled the behaviour of the GAGs/PGs using a modified Blatz –Ko model for a material that can swell, namely
c2m
for lm . 1,
c1m (kPa)
and
c e (kPa)
cm 2 m 2 1 {exp[cm 2 ((l ) 1) ] 1}, 4cm 2
for lkj . 1,
c GAGs (kPa)
W m (lm ) ¼
{exp[ck2 ((lkj )2 1)2 ] 1},
6
rsif.royalsocietypublishing.org
W kj (lkj ) ¼
ck1 4ck2
continuum model
energy functions, namely
Table 1. Best-fit values of model parameters for both the continuum and finite-element models. Note that the finite-element model accounts further for the effect of the GAGs via a swelling induced contribution to the tangent stiffness matrix (cf. [11,27]).
Downloaded from rsif.royalsocietypublishing.org on June 15, 2014
Downloaded from rsif.royalsocietypublishing.org on June 15, 2014
collagen GAGs
elastin
SMCs
circumferential
axial
diagonal
area fraction
0.050
0.199
0.279
0.051
0.065
0.342
media (%) adventitia (%)
80 20
97 3
100 0
0 100
15 85
15 85
fibre families FEM
amorphous matrix
circumferential
axial
diagonal
area fraction media (%)
0.250 80
0.280 5
0.230 15
0.230 15
adventitia (%)
20
95
85
85
Table 3. Values of the inner and outer radii, medial –adventitial border and opening angle for four cases of swelling for both the continuum (CM) and finiteelement (FEM) models. Shown, also, is the percentage difference between the model predictions, which represents the first such inter-model comparison of the effects of intramural swelling. g1 inner radius (mm)
medial – adventitial border (mm)
outer radius (mm)
opening angle (deg)
g2
b1
b2
CM
123
129
145
156
FEM % diff
120 2.43
126 2.60
142 2.07
157 0.93
CM
166
173
191
202
FEM % diff
162 2.77
168 2.99
186 2.72
199 1.34
CM FEM
216 216
220 221
231 237
239 248
% diff
0.15
0.62
2.61
3.76
CM FEM
12.21 12.78
31.13 31.62
78.33 76.79
114.35 114.76
% diff
4.46
1.55
2.01
0.36
describe in the continuum model the fraction of collagen fibres-oriented circumferentially and axially, respectively, were considered free parameters in the nonlinear regression. In addition, based on prior work [26], we prescribed the deposition stretches for each constituent in the continuum model. For the elastic network, the deposition stretch tensor Ge ¼ diag[0:15, 2:80, 2:45] for the radial, circumferential and axial directions, respectively. For the collagen fibres, Gk1 ¼ 1:06 for axial fibres, Gk2 ¼ 1:12 for circumferential fibres (located only in the adventitial layer) and Gk3,4 ¼ 1:1 for diagonal fibres. For the smooth muscle, which was directed circumferentially within the media, we assigned a deposition stretch of Gm ¼ 1.03. Finally, the PGs were assumed to be in traction in the radial direction in the homeostatic configuration to mimic the small swelling observed in vivo and, via incompressibility, compressed in the other two directions, thus yielding a deposition stretch tensor of GPG ¼ diag[1:25, 0:90, 0:90] relative to radial, circumferential and axial directions, respectively. Finally, values of FCD in the finite-element model were studied parametrically based on expected mean basal values of
FCD ¼ 250 mEq l21 in the media and 220 mEq l21 in the adventitia under normal conditions [11,18,30]. Consequences of differences in FCD in the media were considered over the range 250 to 2160 mEq l21, which yielded macroscopic changes in opening angle over a range from 12.8 to 114.8 deg (table 3).
5. Results Figure 3 shows data (mean + s.d.) from a pressure–diameter protocol performed on wild-type murine common carotid arteries at the in vivo axial stretch, consistent with that found in [19]. Shown, also, are best-fit results for the aforementioned constitutive relations for the uni- and bilayered finite-element models as well as the continuum model. Results were similar for data from the other protocols, but are not shown for clarity. Table 1 lists the associated best-fit values of the parameters for all three models and table 2 lists layer-specific (medial versus adventitial) ‘mass’ fractions for the different structurally significant constituents for the
J. R. Soc. Interface 11: 20140397
continuum model
7
rsif.royalsocietypublishing.org
Table 2. Distributions of primary load bearing constituents motivated by microstructural information for the mouse common carotid artery [22,26]. Area fractions were obtained by analysing histological samples using a custom code in Matlab; the percentage of the constituents in each layer was inferred from considerations of the macroscopic behaviour of the vessel in vitro.
Downloaded from rsif.royalsocietypublishing.org on June 15, 2014
150
pressure (mmHg)
90 60 30
380
460 540 620 outer diameter (mm)
700
Figure 3. Representative luminal pressure– outer diameter responses for a simulated (lines) and experimental (closed circles) biaxial pressurization and axial extension protocol, consistent with that in Ferruzzi et al. [19]. Shown, specifically, are the biaxial responses for the CM (solid line) and for both a unilayered (FEM 1 layer; dotted-dashed line) and bilayered (FEM 2 layers; dashed line) finite-element model. bilayered finite-element and continuum models; fractions for the unilayer finite-element model of the wall (see below) necessarily account for the global mass fractions. Figure 4 (left column) shows the relative volume change (detF = 1 during swelling, not transient loading) for both the bilayered finite-element (open squares) and continuum (solid line) models for two cases of swelling that generated two different residual stress-related opening angles F0 . In particular, the distribution of swelling in the continuum model was chosen to mirror phenomenologically the macroscopic consequences of prescribed mechanistic variations in FCD in the finite-element model (250 and 270 mEq l21, respectively). Clearly, the two models yielded similar (i.e. within typical experimental uncertainty) results as desired, F0 in this case. In addition, there was good agreement in other measurable metrics: inner and outer radii as well as location of the medial –adventitial border in the swollen configuration differed by less than 4% between the finite-element and continuum computations (table 3). These findings represent, to the best of our knowledge, the first demonstration of model equivalency for these two very different but complementary approaches, and thereby provide confidence in interpreting additional results from the continuum model. Figure 4b,d also shows predicted transmural distributions of circumferential and axial Cauchy stress obtained with the continuum model. Albeit not shown, the radial stress varied continuously through the wall (as required by traction continuity at the medial–adventitial interface) and was compressive, from negative the luminal pressure (P ¼ 93 mmHg) at the inner wall to zero at the outer wall, which satisfied radial traction boundary conditions. Note, also, the predicted jumps in circumferential and axial stress at the medial–adventitial interface as expected (cf. [12]), and that the mean circumferential stress (obtained via the Laplace relation, Pa/h) agreed well with the computed distributions. Consistent with findings by Bellini et al. [26], which did not include effects of swelling, both the circumferential and axial stresses were nearly uniform within each layer, with mean values of both stresses being higher in the media than in the adventitia at mean arterial pressure and the in vivo axial stretch. This result is consistent
8
J. R. Soc. Interface 11: 20140397
0 300
with the thought that the media bears most of the load under normal conditions, with the adventitia serving primarily as a protective sheath that is progressively engaged as the blood pressure increases acutely above normal levels. Of particular note, however, the increased medial swelling (good case g2 relative to g1) tended to give rise to a larger opening angle consistent with findings in [11] as well as slightly lower values of the stress in the media and slightly higher values in the adventitia. This latter finding suggests a ‘good’ consequence of normal swelling, which could help homogenize the transmural distribution of stress with the adventitia protecting the media to some extent. Figures 5 and 6 show two cases (b1 and b2) that we submit could be potentially bad consequences of intralamellar swelling. In particular, figure 5 shows that increased medial swelling could increase the opening angle further, but also increase the transmural gradient in stress in both the media and adventitia while decreasing the mean stress in the former and increasing it in the latter. This potentially remarkable change in the stress field could alter the mechanobiological responses by both the medial smooth muscle cells and adventitial fibroblasts. Indeed, this effect is even clearer in figure 6, which compares transmural circumferential stresses in the unswollen homeostatic condition with all four swelling cases (g1, g2, b1, b2), for two values of pressure: P ¼ 120 mmHg (figure 6a) and P ¼ 220 mmHg (figure 6b). The higher acute pressure mimics exercise conditions independent of any remodelling that would occur if the increased pressure were sustained. The effect of swelling was exacerbated by the higher pressure, suggesting that in the extreme cases it could increase stress in the adventitial layer by approximately 100 kPa. Figure 7 collects together results associated with those of figures 4 and 5, but shows the transmural distribution of radial stretch (lr ¼ detF/lulz) evaluated for the whole mixture for the different cases of swelling as well as for the case of no swelling (solid grey line) in the homeostatic condition (MAP and liv z ). Recall that, with the unswollen ‘homeostatic’ configuration serving as the reference, Fh ¼ I and the radial stretch is equal to 1 in this case. Hence, increased values of radial stretch owing to swelling can be considered as fold-changes. Given these results, one sees that the intralamellar units would necessarily experience radial expansion (i.e. lr . 1, relative to the unswollen in vivo state), not compression, owing to swelling of the media. Under normal conditions, such intralamellar separation could put radially oriented elastic fibres/ microfibrils in tension, which would be advantageous for mechanosensing and thus good (cf. results in figure 4). Extreme swelling, however, could possibly overextend such fibres/fibrils and result in damage, which would be bad from the perspective of mechanosensing despite the possible existence of favourable macroscopic circumferential and axial stresses within the media (cf. results in figure 5). This is, to the best of our knowledge, the first recognition of these potentially competing effects between the mechanosensing and mechanics within the arterial wall. Indeed, overextending all such radial connections could result in a loss of contact with a medial smooth muscle cell and thereby anoikis (i.e. a form of apoptosis owing to lost connections within the matrix; [17]). Although we do not know the extension at which microfibrils would sever, large elastic fibres appear to fail at extensions on the order of 1.8–2.2 [31], which is well within the range predicted here.
rsif.royalsocietypublishing.org
120
FEM 2 layers FEM 1 layer CM experimental data
Downloaded from rsif.royalsocietypublishing.org on June 15, 2014
(b)
2.10
,FEM = 12.78 F* 0 ,CM = 12.21 F* 0
1.66 1.44
,FEM = 31.62 F* 0 ,CM = 31.13 F* 0
1.88 1.66 1.44
g2
1.22
0
0.2
0.4 0.6 0.8 normalized radius
200 axial stress hoop stress Laplace
160 Cauchy stress (kPa)
relative volume change (v*)
g1
120 80 g2
40
1.0
0
0.2
0.4 0.6 0.8 normalized radius
1.0
Figure 4. Volume change (a,c), relative to ktf (cf. figure 2), for the finite-element (open squares) and continuum (solid line) models for two cases of swelling (g1 and g2, that represent ‘good’ situations associated with FCDMed2 ¼ 250 and 270 mEq l21, respectively), which generated two different residual stress-related opening angles F0 . Note the good agreement between the two different modelling approaches. (b,d) Associated transmural distributions of circumferential (solid line) and axial iv (dashed line) Cauchy wall stress calculated using the continuum model for an axial stretch of liv z ¼ 1:59 (top) and lz ¼ 1:54 (bottom) and an internal pressure of 93.3 mmHg. Also shown is the mean circumferential stress obtained via the classical Laplace relation (dotted-dashed line). Finally, figure 8 shows a finite-element solution for a highly localized pooling of GAGs/PGs that caused a localized swelling within a simple planar model (intramural tangent plane) of the arterial wall (cf. [18]). As it can be seen, a local rather than diffuse accumulation of GAGs/ PGs could result in a significant stress concentration near the leading edge of the pool. Indeed, the predicted local tensile radial stresses (approx. 140 kPa tensile, with associated relative radial stretches . 1) are of sufficient magnitude to cause a local tearing of tissue and thus delamination (cf. [32]). This potentially ugly consequence of the GAGs/PGs re-emphasized that which was presented in [18] based on material properties for the human aorta. The predicted stress concentration was higher for the baseline murine model, however, which reminds us that the wall properties and aspect ratio of the pool can also play significant roles. Moreover, noting that the far-field baseline wall stresses (approx. 100 kPa circumferentially and axially) at mean arterial pressure and the in vivo axial stretch were obtained by applying appropriate displacement boundary conditions at the outer edges of the rectangular domain, the more compliant murine artery required greater prescribed displacements than did the stiffer human artery (i.e. the prescribed displacements were 45% and 30% greater herein in the circumferential and axial directions, respectively). This difference proved to be a challenge for the meshing, with nearly twice as many
elements (150 000 in this model versus 75 000 in [18]) needed to achieve convergence.
6. Discussion It has long been known that vascular smooth muscle cells are highly mechanosensitive; they change their gene expression in response to changes in mechanical loading [33–35]. Indeed, smooth muscle cells express multiple integrins (e.g. a1b1, a5b1 and avb3; [36]), which are known to be important mechanosensors [37,38]. Yet, similar to most studies on other vascular cells [39,40], data on smooth muscle mechanobiology have come primarily from studies of cells plated on two-dimensional substrates, not those seeded within three-dimensional matrices, and, consequently, there has yet to be a careful study of potential mechanosensing via the radial elastic fibres/microfibrils that connect smooth muscle cells to the delimiting elastic laminae in central arteries (cf. [2,21,41]). Indeed, even the precise molecular connection between smooth muscle cells and elastic fibres has yet to be agreed upon. Some have postulated the existence of unique elastin–laminin receptors [42,43], whereas others have focused on both direct and indirect connections via elastinassociated glycoproteins such as the fibulins or fibrillins (e.g. via a5b1 and avb3; [44,45]). Regardless, it is intuitive that mechanical connections between smooth muscle cells and the elastic
J. R. Soc. Interface 11: 20140397
(d)
2.10
1.00
80
0
1.00 (c)
120
40
g1
1.22
axial stress hoop stress Laplace
160 Cauchy stress (kPa)
1.88
9
200
rsif.royalsocietypublishing.org
relative volume change (v*)
(a)
Downloaded from rsif.royalsocietypublishing.org on June 15, 2014
(b)
2.10
F*0 ,FEM = 76.79 ,CM = 78.33 F* 0
1.66 1.44
(d ) ,FEM = 114.76 F* 0 ,CM = 114.35 F* 0
1.88 1.66 1.44 1.22
0.2
0.4 0.6 0.8 normalized radius
axial stress hoop stress Laplace
120 80 40
b2
0
200 160
Cauchy stress (kPa)
relative volume change (v*)
b1
0
2.10
1.00
80 40
b1
1.00 (c)
120
J. R. Soc. Interface 11: 20140397
1.22
axial stress hoop stress Laplace
160 Cauchy stress (kPa)
1.88
10
200
rsif.royalsocietypublishing.org
relative volume change (v*)
(a)
1.0
0
b2
0.2
0.4 0.6 0.8 normalized radius
1.0
Figure 5. Volume change (a,c), relative to ktf (cf. figure 1), for the finite-element (open squares) and continuum (solid line) models for two cases of swelling (b1 and b2, that represent ‘bad’ situations associated with FCDMed2 ¼ 2120 and 2160 mEq l21, respectively), which generated two different residual stress-related opening angles F0 . (b,d) Associated transmural distributions of circumferential (solid line) and axial (dashed line) Cauchy wall stress calculated using the continuum iv model for an axial stretch of liv z ¼ 1:40 (top) and lz ¼ 1:32 (bottom) and an internal pressure of 93.3 mmHg. Also shown is the mean circumferential stress obtained via the classical Laplace relation (dotted-dashed line).
laminae may affect both passive and active mechanics of the wall [46] and it is now clear that the presence of functional elastic fibres exerts a strong direct influence on smooth muscle phenotype [46]. Hence, loss of such connections should be expected to result in phenotypic modulation [47] or, in the extreme cases, apoptosis via anoikis [17]. It is these expectations that motivated the present work. We used both a mixture-based finite-element model [27] to examine gross mechanical consequences of transmural differences in GAG-associated FCD in the arterial wall and a semi-analytical continuum model [29] to examine phenomenologically the biomechanical consequences of the resulting transmural differences in intralamellar swelling. Consistent with general findings on the equivalence of short-term soft tissue responses modelled using a biphasic mixture theory and an incompressible finite elasticity theory [48], we also found that predicted gross mechanical metrics were very similar using these two approaches (table 3 and figure 4). The importance of this cross-validation is that one can exploit a mixture theory to examine physico-chemical mechanisms by which GAGs influence arterial mechanics, but because of the continuing lack of detailed information on regional and longitudinal changes in the FCD in normalcy, ageing and disease, one can also use a mechanically equivalent continuum model to quantify structural and biological
implications of histologically observed differences in regional swelling. This comparison also validated, for the first time, use of the in vivo configuration as a reference given constituent-specific deposition stretches and material properties. We submit that both are important findings. We also compared parametrically two possible classes of transmural distributions of FCD and arterial swelling (uniform versus sigmoidal). Consistent with prior findings, however, we focused on the preferential presence of GAGs within the arterial media [9,11]. As revealed by figure 4, modelling the effects of higher concentrations of GAGs in the media was essential to capture empirically observed increases in the residual stress-related opening angle (cf. [11]). Of course, although carotid arteries dissect [49], it is the thoracic aorta that is most vulnerable [50]. Hence, although the present results should hold qualitatively for any elastic artery, our use of data from the mouse carotid was dictated primarily by the availability of relevant data. There is a pressing need, therefore, to perform similar studies based directly on data for the thoracic aorta, which has, for example, a higher basal FCD than the abdominal aorta [30] and a higher propensity to dissect spontaneously [8]. It has long been known from the cartilage literature that local GAG-induced swelling pressures can cause nearby connective tissue fibres to be in tension despite an overall
Downloaded from rsif.royalsocietypublishing.org on June 15, 2014
2.10
at P = 120 mmHg
400
homogen g1 g2 b1 b2
300 200
Hoop Cauchy stress (kPa)
1.38
0.90
at P = 220 mmHg
0
0.2
0.4 0.6 0.8 normalized radius
1.0
Figure 7. Transmural distribution of radial stretch for the different cases of swelling considered in figures 4 and 5: g1, dashed line; g2, dotted line; b1, dotted-dashed line; b2, solid black line and the case of no swelling (solid grey line). Because the homeostatic configuration is taken as a reference, the radial stretches represented here are ‘real’ not ‘relative’ and the numerical values can be interpreted as fold-changes owing to swelling.
400 300 200 100
0
0.2
0.4 0.6 0.8 normalized radius
1.0
Figure 6. Representative transmural distributions of circumferential Cauchy stress for (a) systolic and (b) acutely elevated blood pressure, each evaluated for five different degrees of swelling: good g1 (dashed line) and g2 (dotted line) distributions (see figure 4), bad b1 (dotted-dashed line) and b2 (solid black line) distributions (figure 5), and a non-swollen homogeneous (grey solid line) distribution, i.e. J ¼ 1 at every point of the wall. The effect of the intramural swelling was similar in some ways to an increase in luminal pressure: it increased the tensile stress in the adventitia while decreasing it in the medial layer. The same trend can be seen for the mean arterial pressure (MAP) by comparing circumferential values of the Cauchy stress in figures 4 and 5. compressive state of stress [16], but we are the first to suggest the potential role of this mixture phenomenon in smooth muscle mechanosensing within the arterial wall (cf. figure 6). As revealed by the exquisite transmission electron micrographs reported by Davis [21], the elastic laminae of large arteries connect to the resident smooth muscle cells via elastic fibre/microfibrillar extensions that consist primarily of elastinassociated fibrillin-1. Although mechanobiological roles of these connections have not been elucidated, we submit that they could serve an important role in mechanosensing, particularly because they tend to colocalize with dense plaques (i.e. focal adhesions) and co-align with associated intracellular contractile proteins [21]. Indeed, Bunton et al. [47] hypothesized that ‘loss of cell attachments signals a nonproductive program to synthesize and remodel an elastic matrix’ in Marfan syndrome, which results from mutations in the gene encoding fibrillin-1. This concept of adverse effects of compromised elastic fibre–smooth muscle connections is also consistent with the demonstrated role of elastin in influencing smooth muscle phenotype. Karnik et al. [46] showed using an elastin haploinsufficient mouse that ‘elastin stabilizes the
arterial structure by inducing a quiescent contractile state in vascular smooth muscles’. That is, interactions between elastic fibres and smooth muscle cells are essential determinants of smooth muscle phenotype, thus it follows that altered interactions can result in phenotypic modulation and thus altered intracellular signalling and resulting gene products. Noting that O’Connell et al. [41] reported that 29% of all medial elastin is found within the intralamellar spaces, either connecting laminae to the cells or connecting laminae to laminae, suggests that we must understand better both the structural and the instructional roles of non-laminar elastic fibres. The present simulations suggest, for the first time, that these future studies should include the related effects of the GAGs/PGs in stretching these fibres. That is, as should be the case, theory/ computation should guide experiments. We are not the first to suggest, however, that GAGs/PGs might affect the mechanobiology of arterial smooth muscle cells. Evanko et al. [51] noted that aggregates of hyaluronan and versican tend to localize near cells but away from focal adhesions, and suggested that ‘by altering the tension of the cell-matrix interaction, this ‘replusive’ but relatively malleable quality of the HA–PG (hyaluronan –proteoglycan) matrix could potentially perturb matrix-induced signalling that occurs via integrin –fibronectin or integrin –collagen interactions’. Our computational results are the first to support such speculations, however. Although we agree with these prior suggestions, we submit further that these aggregates could similarly perturb signalling that is mediated directly via the smooth muscle cell –elastic fibre connections, not just integrins targeting fibronectin or collagen. Although each of these effects could be mediated solely via the mechanical effects of aggregated hyaluronan and versican (i.e. Donnan swelling), it should also be noted that excessive GAGs/PGs may also interfere sterically with normal cellular interactions with matrix via focal adhesions [1]. Hence, there is a need to investigate broadly how GAG/PG aggregates affect diverse aspects of mechanosensing. Indeed, notwithstanding our current focus, other mechano- and chemo-mediated mechanisms certainly contribute further to
J. R. Soc. Interface 11: 20140397
500
1.62
11
1.14
100 0
(b)
homogen g1 g2 b1 b2
1.86
radial stretch
500
rsif.royalsocietypublishing.org
Hoop Cauchy stress (kPa)
(a)
Downloaded from rsif.royalsocietypublishing.org on June 15, 2014
(a) FE model
(c)
(d)
y
radial stress (kPa)
140 117 94 71 48 25 2 –21 –44 –67 –90
circumferential stress (kPa)
250 216 182 148 114 80 46 12 –22 –56 –90
axis stress (kPa)
250 216 182 148 114 80 46 12 –22 –56 –90
Figure 8. (a) Final quarter-symmetric finite-element model (approx. 150 000 hexahedral elements) of the overall 250 250 84 mm rectangular portion of the murine carotid wall, including a single centrally located GAG/PG inclusion (volume of 574.7 mm3). The value of the fixed charge density equals 2160 mEq l21 for the GAG inclusion and 250 mEq l21 for the non-affected aortic wall. The far-field stresses correspond to those developed in the wall owing to the homeostatic blood pressure evaluated in a unilayer cylindrical simulation (radial ¼ 26 kPa, circumferential and axial ¼ 100 kPa) and induced here via appropriate displacement boundary conditions. (b–d) Shown is a greyscale representation of multi-axial stresses in and near the GAG/PG inclusion. Note, in particular, the transition of radial stress from compressive to highly tensile near the leading edge of the inclusion.
J. R. Soc. Interface 11: 20140397
(b)
x
12
rsif.royalsocietypublishing.org
z
the integrity of elastic fibre–smooth muscle interactions. As noted by Sherratt [52], ‘mechanical fatigue has also been postulated as a mechanism of age-related elastic fibre failure, with indirect evidence suggesting that the microfibril/elastin interface may be the main weak point in the composite elastic fibre’. Increased proteolysis via elastases (cf. [53]) would also be expected to affect first the thinner elastic fibres and thus the microfibrillar connections (cf. [47]). Indeed, blocking lysyl oxidase in rats revealed that the radial elastic fibres were lost prior to any notable damage to the thicker elastic fibres that constitute the elastic laminae [20]. There are, therefore, multiple means by which smooth muscle cells might lose their mechanosensing capability via compromised microfibrillar connections. Clearly, there is a need to consider roles of radial stresses/strains, not just the traditionally considered correlations between altered circumferential and axial stresses and arterial growth and remodelling (cf. [54]). Finally, we would be remiss not to mention the potential importance of the present findings to Marfan syndrome and related diseases of the thoracic aorta. This connective tissue disorder results from mutations in the gene encoding fibrillin-1, the primary glycoprotein constituting the elastinassociated microfibrils [55]. Like most extracellular matrix constituents, fibrillin-1 plays multiple roles. It appears to confer long-term stability to elastic fibres, that is, to decrease elastic fibre susceptibility to both mechanical damage (fatigue) and chemical insult ( proteolysis); it participates in sequestering latent TGF-b, which thereby affects many cell behaviours, including control of matrix turnover; and it appears to be the main intralamellar connection between smooth muscle cells and the elastic laminae [21], which implies a role in mechanosensing. Although loss of all three functions of microfibrils probably contributes to aortic dilatation and dissection in Marfan syndrome, we recently suggested that associated increases in TGF-b activity may promote the characteristic increased accumulation of GAGs/PGs within the ascending aorta (cf. [56]), because TGF-b is a known promoter of the synthesis of GAGs/PGs by smooth muscle cells [3,57]. We suggested further that localized pooling of GAGs/PGs could result in stress concentrations that initiate delaminations of the media, which in turn could initiate dissections. This hypothesis was supported by simulations reported in [18]. In this paper, we extended our prior findings by suggesting that diffuse accumulations of GAGs/PGs may also perturb or eliminate smooth muscle mechanosensing via microfibrils, whether they are compromised by genetic mutations, fatigue, proteolysis or other factors, and in the extreme cases lead to anoikis. This latter possibility is consistent with the suggestion of Allaire et al. [58] that ‘enhanced TGF-b signalling in MFS (Marfan syndrome) is thought to increase VSMC (vascular smooth muscle cell) apoptosis and impair aortic healing because of hyaluronan up-regulation’. Regarding the control of GAG/PG production within the arterial wall, cyclic mechanical stress, angiotensin-II and TGF-b all promote the synthesis of GAGs/PGs by arterial smooth muscle cells [1,3,59]. Blocking effects of Ang-II via the angiotensin-II type 1 receptor (e.g. via losartan) or blocking TGF-b via neutralizing antibodies could thus decrease the synthesis of GAGs/PGs and, notably, attenuate the incidence of GAG-related dissection. Caution must be exercised in drug treatments, however, because TGF-b is also involved in the maintenance of adventitial collagen, which serves as a
Downloaded from rsif.royalsocietypublishing.org on June 15, 2014
Funding statement. This work was supported, in part, by grants from the NIH (R01 HL086418, R01 HL105297 and U01 HL116323) and the National Marfan Foundation.
References 1.
2.
3.
4.
5.
6.
7.
Wight NT. 1989 Cell biology of arterial proteoglycans. Arteriosclerosis 9, 1 –20. (doi:10. 1161/01.ATV.9.1.1) Dingemans KP, Teeling P, Lagendijk JH, Becker AE. 2000 Extracellular matrix of the human aortic media: an ultrastructural histochemical and immunohistochemical study of the adult aortic media. Anat. Rec. 258, 1–14. (doi:10.1002/ (SICI)1097-0185(20000101)258:1,1::AID-AR1.3. 0.CO;2-7) Kinsella MG, Bressler SL, Wight TN. 2004 The regulated synthesis of versican, decorin, and biglycan: extracellular matrix proteoglycans that influence cellular phenotype. Crit. Rev. Eukary. Gene Express. 14, 203–234. (doi:10.1615/CritRevEukaryot GeneExpr.v14.i3.40) Wagenseil JE, Mecham RP. 2009 Vascular extracellular matrix and arterial mechanics. Physiol. Rev. 89, 957–989. (doi:10.1152/physrev.00041. 2008) Reynertson RH, Parmley RT, Roden L, Oparil S. 1986 Proteoglycans and hypertension: I. A biochemical and ultrastructural study of aorta glycosaminoglycans in spontaneously hypertensive rats. Coll. Relat. Res. 6, 77– 101. (doi:10.1016/ S0174-173X(86)80033-4) Rasmussen LM, Wolf YG, Ruoslahti E. 1995 Vascular smooth muscle cells from injured rat aortas display elevated matrix production associated with transforming growth factor-beta activity. Am. J. Pathol. 147, 1041 –1048. Yang SNY, Burch ML, Tannock LR, Evanko S, Osman N, Little PJ. 2010 Transforming growth factor-beta regulation of proteoglycan synthesis in vascular smooth muscle: contribution to lipid binding and accelerated atherosclerosis in diabetes. J. Diabetes 2, 233– 242. (doi:10.1111/j.1753-0407. 2010.00089.x)
8.
9.
10.
11.
12. 13.
14.
15.
16.
17.
Humphrey JD. 2013 Possible mechanical roles of glycosaminoglycans in thoracic aortic dissection and associations with dysregulated transforming growth factor-b. J. Vasc. Res. 50, 1–10. (doi:10.1159/ 000342436) Lark MW, Yeo T-K, Mar H, Lara S, Hellstrom I, Hellstrom K-E, Wight TN. 1988 Arterial chondroitin sulfate proteoglycan: Localization with a monoclonal antibody. J. Histochem. Cytochem. 36, 1211 –1221. (doi:10.1177/36.10.3047228) Wight NT. 2008 Arterial remodeling in vascular disease: a key role for hyaluronan and versican. Front. Biosci. 13, 4933–4937. (doi:10.2741/3052) Azeloglu EU, Albro MB, Thimmappa VA, Ateshian GA, Costa KD. 2008 Heterogeneous transmural proteoglycan distribution provides a mechanism for regulating residual stress in the aorta. Am. J. Physiol. 294, H1197– H1205. Humphrey JD. 2002 Cardiovascular solid mechanics: cells, tissues, and organs. New York, NY: Springer. Tammi MI, Day AJ, Turley EA. 2002 Hyaluronan and homeostasis: a balancing act. J. Biol. Chem. 277, 4581 –4584. (doi:10.1074/jbc.R100037200) Tinoco I, Sauer K, Wang JC, Puglisi JD. 2002 Physical chemistry, principles and applications in biological sciences. Hoboken, NJ: Prentice Hall. Ateshian GA, Rajan V, Chahine NO, Canal CE, Hung CT. 2009 Modeling the matrix of articular cartilage using a continuous fiber angular distribution predicts many observed phenomena. J. Biomech. Eng. 131, 061003. (doi:10.1115/1.3118773) Maroudas A. 1976 Balance between swelling pressure and collagen tension in normal and degenerate cartilage. Nature 260, 808–809. (doi:10.1038/260808a0) Michel J-B. 2003 Anoikis in the cardiovascular system. Known and unknown extracellular mediators. Arterioscler. Thromb. Vasc. Biol. 23,
18.
19.
20.
21.
22.
23.
24.
25.
26.
2146– 2154. (doi:10.1161/01.ATV.0000099882. 52647.E4) Roccabianca S, Ateshian G, Humphrey JD. 2013 Biomechanical roles of medial pooling of glycosaminoglycans in thoracic aortic dissection. Biomech. Model. Mechanobiol. 13, 13 –25. (doi:10. 1007/s10237-013-0482-3) Ferruzzi J, Bersi MR, Humphrey JD. 2013 Biomechanical phenotyping of central arteries in health and disease: advantages of and methods for murine models. Ann. Biomed. Eng. 41, 1311–1130. (doi:10.1007/s10439-013-0799-1) Nakashima Y, Sueishi K. 1992 Alteration of elastic architecture in the lathyritic rat aorta implies the pathogenesis of aortic dissecting aneurysm. Am. J. Pathol. 140, 959–969. Davis EC. 1993 Smooth muscle cell to elastic lamina connections in developing mouse aorta: role in medial organization. Lab Invest. 68, 89 –99. Bersi M, Collins MJ, Wilson E, Humphrey JD. 2012 Disparate changes in the mechanical properties of murine carotid arteries and aorta in response to chronic infusion of angiotensin-II. Int. J. Adv. Eng. Sci. Appl. Math. 4, 228– 240. (doi:10.1007/s12572012-0052-4) Humphrey JD, Eberth JF, Dye WW, Gleason RL. 2009 Fundamental role of axial stress in compensatory adaptations by arteries. J. Biomech. 42, 1 –8. (doi:10.1016/j.jbiomech.2008.11.011) Davis EC. 1995 Elastic lamina growth in the developing mouse aorta. J. Histochem. Cytochem. 43, 1115 –1123. (doi:10.1177/43.11.7560894) Chuong CJ, Fung YC. 1986 Residual stress in arteries. J. Biomech. Eng. 108, 189 –192. (doi:10. 1115/1.3138600) Bellini C, Ferruzzi J, Roccabianca S, DiMartino E, Humphrey JD. 2014 A microstructurally-motivated model of arterial wall mechanics with
13
J. R. Soc. Interface 11: 20140397
In conclusion, as with so many things in life, too much of a good thing can be a bad thing. For example, Galis & Khatri [64] noted that, at normal concentrations, active matrix metalloproteinases are essential for tissue homeostasis, yet when these concentrations deviate from normal the consequences can be dire. We similarly submit that normal distributions and concentrations of GAGs/PGs within the arterial wall are essential for proper mechanical and biological functions, including homeostasis, yet, when increased above normal levels, these same GAGs/PGs may compromise or eliminate particular aspects of normal mechanosensing, or, if pooled, they may initiate potentially lethal dissections. There is, therefore, a pressing need to increase our understanding of GAGs/PGs and their potentially good, bad and ugly roles in arterial health and disease, which will require increased experimental information and advances in computational modelling.
rsif.royalsocietypublishing.org
protective sheath for the arterial wall [60]. Nevertheless, whereas normal mechanical loading of radial elastic fibres could promote matrix homeostasis within the media, loss of such loading and thus loss of mechanosensing could result in marked changes in matrix, including the early upregulation of hyaluronan and versican as is typical in wound healing [61]. That is, without proper mechanosensing, the resident cells will not be able to maintain, remodel or repair the extracellular matrix that is fundamental to the overall structural integrity of the wall [62]. The complexity of such interactions among cell membrane bound receptors and their signalling pathways is mirrored by the complexity of extracellular matrix interactions [63]. For example, versican not only interacts with hyaluronan, but it also interacts with the fibulins and fibrillins that associate with elastin [55], and it has been suggested that increased versican could compromise effective elastic fibre assembly [3]. Clearly, the potential consequences of altered GAGs/PGs can go well beyond just the mechanics, although ultimately it is the mechanics that determines the fate of an artery (e.g. dissection or rupture).
Downloaded from rsif.royalsocietypublishing.org on June 15, 2014
28.
30.
31.
32.
33.
34.
35.
36.
37.
38.
39.
41.
42.
43.
44.
45.
46.
47.
48.
49.
50.
51.
52.
53.
54.
55.
56.
57.
58.
59.
60.
61.
62.
63.
64.
factor stimulates the formation of versican –hyaluronan aggregates and pericellular matrix expansion in arterial smooth muscle cells. Arch. Biochem. Biophys. 394, 29 –38. (doi:10.1006/ abbi.2001.2507) Sherratt MJ. 2009 Tissue elasticity and the ageing elastic fibre. Age 31, 305 –325. (doi:10.1007/ s11357-009-9103-6) Hornebeck W, Emonard H. 2011 The cell –elastin – elastase(s) interacting triade directs elastolysis. Front. Biosci. 16, 707–722. (doi:10.2741/3714) Humphrey JD. 2008 Vascular adaptation and mechanical homeostasis at tissue, cellular, and subcellular levels. Cell Biochem. Biophys. 50, 53 –78. (doi:10.1007/s12013-007-9002-3) Kielty CM. 2006 Elastic fibres in health and disease. Expert Rev. Mol. Med. 8, 1–17. (doi:10.1017/ S146239940600007X) Nataatmadja M, West J, West M. 2006 Overexpression of transforming growth factor-beta is associated with increased hyaluronan content and impairment of repair in Marfan syndrome aortic aneurysm. Circulation 114, I371 –I377. (doi:10. 1161/circulationaha.105.000927) Chen JK, Hoshi H, McKeehan WL. 1987 Transforming growth factor type beta specifically stimulates synthesis of proteoglycan in human adult arterial smooth muscle cells. Proc. Natl Acad. Sci. USA 84, 5287– 5291. (doi:10.1073/pnas.84.15.5287) Allaire E et al. 2009 New insight in aetiopathogenesis of aortic diseases. Eur. J. Vasc. Endovasc. Surg. 37, 531– 537. (doi:10.1016/j.ejvs. 2009.02.002) Lee RT et al. 2001 Mechanical strain induces specific changes in the synthesis and organization of proteoglycans by vascular smooth muscle cells. J. Biol. Chem. 276, 13 847– 13 851. Li W et al. 2014 Tgfbr2 disruption in postnatal smooth muscle impairs aortic wall homeostasis. J. Clin. Invest. 124, 755 –767. (doi:10.1172/ JCI69942) Wight TN, Potter-Perigo S. 2011 The extracellular matrix: an active or passive player in fibrosis. Am. J. Physiol. 301, G950–G955. Humphrey JD, Milewicz DM, Tellides G, Schwartz MA. 2014 Dysfunctional mechanosensing in aneurysms. Science 344, 477– 479. (doi:10.1126/ science.1253026) Midwood KS, Schwarzbauer JE. 2002 Elastic fibers: building bridges between cells and their matrix. Curr. Biol. 12, R279 –R281. (doi:10.1016/S09609822(02)00800-X) Galis ZS, Khatri JJ. 2002 Matrix metalloproteinases in vascular remodeling and atherogenesis: the good, the bad, and the ugly. Circ. Res. 90, 251 –262.
14
J. R. Soc. Interface 11: 20140397
29.
40.
clinical perspectives. Physiol. Rev. 91, 327 –387. (doi:10.1152/physrev.00047.2009) Chiquet M, Gelman L, Lutz R, Maier S. 2009 From mechanotransduction to extracellular matrix gene expression in fibroblasts. Biochim. Biophys. Acta 1793, 911 –920. (doi:10.1016/j.bbamcr.2009. 01.012) O’Connell MK et al. 2008 The three-dimensional micro- and nanostructure of the aortic medial lamellar unit measured using 3-D confocal and electron microscopy. Matrix Biol. 27, 171– 181. (doi:10.1016/j.matbio.2007.10.008) Spofford CM, Chilian WM. 2001 The elastin – laminin receptor functions as a mechanotransducer in vascular smooth muscle. Am. J. Physiol. Heart Circ. Physiol. 280, 1354–1360. Arribas SM, Hinek A, Gonza´lez MC. 2006 Elastic fibres and vascular structure in hypertension. Pharmacol. Ther. 111, 771–791. (doi:10.1016/j. pharmthera.2005.12.003) Bax DV, Bernard SE, Lomas A, Morgan A, Humphries J, Shuttleworth CA, Humphries MJ, Kielty CM. 2003 Cell adhesion to fibrillin-1 molecules and microfibrils is mediated by a5b1 and avb3 integrins. J. Biol. Chem. 278, 34605–34616. (doi:10.1074/jbc.M303159200) Bezie Y, Lacolley P, Laurent S, Gabella F. 1998 Connection of smooth muscle cells to elastic lamellae in aorta of spontaneously hypertensive rats. Hypertension 32, 166– 169. (doi:10.1161/01. HYP.32.1.166) Karnik SK, Brooke BS, Bayes-Genis A, Sorensen L, Wythe JD, Schwartz RS, Keating MT, Li DY. 2003 A critical role for elastin signaling in vascular morphogenesis and disease. Development 130, 411 –423. (doi:10.1242/dev.00223) Bunton TE, Biery NJ, Myers L, Gayraud B, Ramirez F, Dietz HC. 2001 Phenotypic alteration of vascular smooth muscle cells precedes elastolysis in a mouse model of Marfan syndrome. Circ. Res. 88, 37 –43. (doi:10.1161/01.RES.88.1.37) Ateshian G, Ellis BJ, Weiss JA. 2007 Equivalence between short-time biphasic and incompressible elastic material responses. J. Biomech. Eng. 129, 405 –412. (doi:10.1115/1.2720918) Stevic I, Chan HHW, Chan AKC. 2011 Carotid artery dissections: thrombosis of the false lumen. Thromb. Res. 128, 317–324. (doi:10.1016/j.thromres.2011.06.024) Grond-Ginsbach C, Pjonteck R, Aksay SS, Hyhlik-Durr A, Bockler D, Gross-Weissmann M-L. 2010 Spontaneous arterial dissection: Phenotype and molecular pathogenesis. Cell Mol. Life Sci. 67, 799 –1815. (doi:10.1007/s00018-010-0276-z) Evanko SP, Johnson PY, Braun KR, Underhill CB, Dudhia J, Wight TN. 2001 Platelet-derived growth
rsif.royalsocietypublishing.org
27.
mechanobiological implications. Ann. Biomed. Eng. 42, 488–502. (doi:10.1007/s10439-013-0928-x) Ateshian G, Albro MB, Maas S, Weiss JA. 2011 Finite element implementation of mechanochemical phenomena in neutral deformable porous media under finite deformation. J. Biomech. Eng. 133, 081005. (doi:10.1115/1.4004810) Holzapfel GA, Gasser TC, Ogden RW. 2000 A new constitutive framework for arterial wall mechanics and a comparative study of material models. J. Elast. 61, 1 –48. (doi:10.1023/A:1010835316564) Demirkoparan H, Pence TJ. 2007 Swelling of an internally pressurized nonlinearly elastic tube with fiber reinforcing. Int. J. Solid Struct. 44, 4009– 4029. (doi:10.1016/j.ijsolstr.2006.11.006) Porterfield SP, Calhoon TB, Weiss HS. 1968 Changes in connective tissue colloidal charge density with atherosclerosis and age. Am. J. Physiol. 215, 324–329. Lillie MA, Gosline JM. 2007 Limits to the durability of arterial elastic tissue. Biomaterials 28, 2021–2031. (doi:10.1016/j.biomaterials.2007.01. 016) MacLean NF, Dudek NL, Roach MR. 1999 The role of radial elastic properties in the development of aortic dissection. J. Vasc. Surg. 29, 703–710. (doi:10. 1016/S0741-5214(99)70317-4) Leung DY, Glagov S, Matthews MB. 1976 Cyclic stretching stimulates synthesis of matrix components by arterial smooth muscle cells in vitro. Science 191, 475–477. (doi:10.1126/science. 128820) Li Q, Muragaki Y, Hatamura I, Ueno H, Ooshima A. 1998 Stretch-induced collagen synthesis in cultured smooth muscle cells from rabbit aortic media and a possible involvement of angiotensin II and transforming growth factor-b. J. Vasc. Res. 35, 93 –103. (doi:10.1159/000025570) Li C, Xu Q. 2007 Mechanical stress-initiated signal transduction in vascular smooth muscle cells in vitro and in vivo. Cell Signal. 19, 881 –891. (doi:10.1016/ j.cellsig.2007.01.004) Moiseeva EP. 2001 Adhesion receptors of vascular smooth muscle cells and their functions. Cardiovasc. Res. 52, 372 –386. (doi:10.1016/S0008-6363(01) 00399-6) DuFort CC, Paszek MJ, Weaver VM. 2011 Balancing forces: architectural control of mechanotransduction. Nat. Rev. Mol. Cell Biol. 12, 308– 319. (doi:10.1038/ nrm3112) Schwartz MA. 2010 Integrins and extracellular matrix in mechanotransduction. Cold Spring Harb. Perspect. Biol. 2, a005066. (doi:10.1101/cshperspect. a005066) Chiu J-J, Chien S. 2011 Effects of disturbed flow on vascular endothelium: pathophysiological basis and