|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Departments of 1Biomedical Engineering, 2Surgery, and 3Cellular and Integrative Physiology, Indiana University Purdue University Indianapolis, Indianapolis, Indianapolis
Submitted 17 February 2007 ; accepted in final form 15 May 2007
| ABSTRACT |
|---|
|
|
|---|
endothelium; stress; deformation; mechanotransduction
An arterial wall is typically composed of three functional layers from inner to outer wall surfaces: intima, media, and adventitia. The intima comprises a single layer of endothelial cells and a basal lamina in normal vessels. The media consists mainly of smooth muscle cells, elastin, and collagen fibrils. The adventitia contains mainly collagen fibers, ground substance, fibroblasts, and fibrocytes. Since the intima is very thin and much more compliant than the media and the adventitia, its material property is generally ignored. Consequently, blood vessels are typically viewed either as a single layer of homogeneous material (3, 5, 19, 32) or as a two-layer composite of intima-media and adventitia under normal conditions (16, 26, 27, 34, 39). It has been well accepted that endothelial cells (the major constituent of intima) are very sensitive to changes in flow-induced shear stress, and many biological activities are related to this shear loading (1, 4, 7, 30, 31, 35, 38, 40). For instance, the diameter and thickness of the vessel wall are modulated by flow-induced shear stress, even though the shear stress is several orders of magnitude smaller than the circumferential wall stress. This raises the question of why blood vessels are so sensitive to such small loading. Our hypothesis is that shear strain can be substantial in endothelial cells, despite the small shear stress. Our rationale is that the deformation of endothelium can be underestimated if the mechanical property of intima is not considered individually.
Most previous investigations on the mechanics of vessel wall have focused on the pressure-induced normal stress and strain. The shear strain, however, due to the relatively small magnitude of flow-induced shear stress, has not received significant attention. In this study, we will examine the shear deformation in arterial walls under physiological pressure and flow conditions. We will demonstrate that the flow-induced shear strain in the intima layer can be of the same magnitude as normal (circumferential and axial) strains caused by pressure and axial prestretch. This result may shed light on mechanotransductory mechanisms in blood vessels, as well as arterial growth and remodeling modulated by the mechanical environment of endothelial cells.
Glossary




=
/(
–
)


P

= 1,...,7)
m
, z)
w
w = dw/dR
, z)
jk
jk
| MATERIALS AND METHODS |
|---|
|
|
|---|
characterizes the circumferential residual strain (5) at the no-load state (i.e., a tube with zero pressure and shear). The deformed cross section is schematically shown in Fig. 1C, where the axial (z) direction is pointing out of the r-
plane. We employ a two-layer model in which the intima is considered different from the rest of the wall; i.e., the media-adventitia layer (Fig. 1D). This treatment differs from the classic intima-media and adventitia two-layer model (16, 39), since the significance of intima layer is of interest in this work.
|
plane is zero throughout the arterial wall. Thus the nontrivial equations of equilibrium on the deformed configuration are simplified as:
![]() | (1) |
![]() | (2) |
jk (j, k=r,
, z) denotes Cauchy stress. Using the approaches developed by Humphrey and Na (19) and Rachev et al. (32) while considering the boundary condition
rr (ri) = –P (where P is blood pressure), the following expressions are obtained (see APPENDIX for details):
![]() | (3) |
![]() | (4) |
![]() | (5) |
, z) is the Green strain (see Eq. A4 in APPENDIX),
m (m = r,
, z) are stretch ratios, and
w = dw/dR with w being the displacement that induces shear strain. It is noted that W is different for the arterial layers, and hence two sets of parameters are needed for the intima and media-adventitia, respectively.
If we consider the boundary condition
rz (ri) = –
(where
is the blood shear stress at the inner surface, see Fig. 1A), the solution of Eq. 2 can be obtained as
rz (r) = –
ri/r, which implies that the perivascular tethering tissue must provide a shear force –
ri/ro on the outer surface to balance the flow-induced shear load. The shear stress at given radial coordinate r can be used to compute
w from Eq. A11:
![]() | (6) |
![]() | (7) |
Simulation parameters.
The dimensions of the artery at zero-stress state were taken as Ri = 12.78 mm, Ro = 13.12 mm,
= 163° (Kassab et al., unpublished data corresponding to Ref. 39). The intima (2
5 µm thick) was assumed to occupy 2% of the wall thickness (100
250 µm). The axial stretch ratio was chosen as
z = 1.4 (27). The Fung-model parameters C and b1
b7 in the intact porcine left anterior descending (LAD) artery are listed in Table 1, as previously reported by our group (39). The porcine coronary artery elastic modulus was estimated to be 170 kPa at physiological loading based on experimental data (26). According to Sato et al. (36), the Young's modulus of porcine aortic endothelial cells is on the order of 0.1 kPa. Hence the modulus in endothelium is
1,000 times smaller than the intact artery. Since complete model parameters are not available for intima, the parameter C was taken to be 1% and b
(
= 1 ,... , 7) to be 10% of those of the intact porcine LAD artery reported by Wang et al. (39). To maintain the overall elastic stiffness of the artery at P = 100 mmHg, parameter C in the media-adventitia was increased by 1.65% from the value in intact wall (Table 1). Note that parameter C in the Fung model is a scale factor, and parameter b
reflects the material nonlinearity (curvature of stress-strain relation). Both C and b
have been rescaled to make the elastic modulus (proportional to Cb
) of intima 1,000 times smaller. Sato et al. (36) considered a linear viscoelastic model for endothelial cells. Here we have chosen small b
(Table 1) to ensure nearly linear stress-strain relationship of intima. More details about parameter selection are provided in the APPENDIX.
|
rr(ro) in Eq. 3] and the prescribed pressure boundary condition is <10–8 kPa. A simple test was done to verify that the unloaded state [P = 0,
rr(ro) = 0, and
z = 1.0] can be achieved with the above parameters (the resulting no-load ri = 1.06 mm and ro = 1.39 mm), for which the residual stresses and strains are distributed across the arterial wall in a pattern similar to those shown in Ref. 5; e.g., circumferential components are compressive at the inner surface and tensile at the outer surface.
|
rr(ri) = –13.33 kPa (the average value of systole and diastole pressure, P = 100 mmHg), while the constraint at the outer surface was taken to be
rr(ro) = –8.67 kPa (65 mmHg), which stems from the tethering myocardium (14). The choice of the external pressure was based on the consideration of myocardial constraint and the resulting outer radius of the vessel (1.81 mm at loaded state), consistent with the reported pressure-diameter relation of in situ porcine LAD artery (14). | RESULTS |
|---|
|
|
|---|
Under simultaneous internal pressure P = 100 mmHg, external pressure 65 mmHg, and blood shear stress
= 1.5 Pa (15 dyn/cm2), it is found that the shear strain is significant in the intima (Erz
–0.08), but close to zero in media-adventitia (Fig. 3A). The shear strain increases sharply at the interface of intima and media-adventitia (at normalized wall thickness of 2%, Fig. 3C), so does the radial strain (slightly) due to the contribution of
w in Err (i.e., discontinuous
w results in discontinuous Err, Eq. A4). The radial stress is compressive in the whole vessel wall, but the circumferential and axial stresses change abruptly from compressive to tensile across the interface of intima and media-adventitia layer (Fig. 3B). All normal stresses are nearly equal to the negative blood pressure (13.33 kPa) in the intima layer (Fig. 3D), indicating that the intima is under a hydrostatic state approximately. The shear stress is always small, however, as expected (see right-hand side of Eq. 6).
|
For the same material properties in Table 1, the circumferential and axial strains do not change noticeably from Fig. 3A when the applied shear stress
is increased to 9 Pa. However, the magnitude of shear strain in intima increases almost linearly with
, and it approaches the magnitude of normal strains, e.g., | Erz |
0.4 at
= 9 Pa (Fig. 4A). The radial strain in media-adventitia is also similar to Fig. 3A, whereas its magnitude decreases nonlinearly with
in intima (Fig. 4A), because the increase of the term
w2 in Eq. A4 leads to increasing Err (its magnitude decreases since Err is negative). This result implies that shear and radial strains in endothelial cells strongly depend on the flow conditions (normal flow vs. flow overload).
|
1.5 Pa). Figure 4B shows that the magnitude of shear strain on the inner surface increases (but radial strain decreases) nonlinearly with decreasing b7 for intima (other parameters were kept unchanged). At b7 = 0.05, the magnitude of shear strain (| Erz | = 0.44) approaches that of the circumferential and axial strains. The shear stress, on the other hand, remains several orders of magnitude smaller than the normal stresses, even in flow overload. | DISCUSSION |
|---|
|
|
|---|
Strain or stress as mechanical stimulus. Experiments by Lu and Kassab (25) revealed that the forward flow promotes nitric oxide production in arterial wall under normal flow, but significantly reduces nitric oxide concentration under reversed flow. This phenomenon has been attributed to increased production of superoxide anion in reversed flow. The biological response of endothelial cells to the same magnitude of wall shear stress (but opposite in direction) is very different. It may be that the shear strain on the endothelial cells is very different, depending on the direction, as speculated below.
The residual strain in blood vessels is closely related to the remodeling process (10). There is enormous evidence that circumferential residual strain (characterized by the opening angle) in arteries changes with remodeling, e.g., Liu and Fung (24). Davies (7) pointed out that prolonged steady flow results in reorganization of endothelial surface and significant reduction of peak shear stress compared with no-flow conditions, which implies that a residual shear strain may exist in the arterial wall due to remodeling under chronic flow-induced shear. As illustrated in Fig. 5, if residual shear strain exists (left), the shear deformation might be quite different for the same shear stress under forward and reversed flows (middle and right). The flow-dependent reorientation of endothelial cells (7, 30) could be related to the residual shear strain. Nonetheless, three-dimensional (3D) reconstructions of endothelial cells under various shear loadings are needed to confirm this hypothesis.
|
Rachev (33), Fridez et al. (8), and Gleason and Humphrey (11) have studied the remodeling of carotid artery, including passive and active behaviors, by considering blood flow and pressure loading. The potential residual shear strain may be included in such analyses, as it could play an important role in interpreting flow-related remodeling and mechanotransduction.
Critique of method.
Theoretically, if material property does not change along the axial direction, the arterial diameter will decrease along the flow direction. This is because flow is a result of pressure gradient. Using Poiseuille's law (10), the internal shear stress can be approximated by
![]() | (8) |
is the rate of flow, and L equals the length along axis z when pressure drop is
P. Therefore, all of the variables depend on coordinate z as well as on r if flow exists. Nevertheless, this dependence on z can be viewed as a secondary effect in a typical blood vessel segment of interest (L
1 cm), because the pressure drop is small and the change of lumen radius is insignificant. Thus our assumption that all variables are independent of z is reasonable.
Since the 3D mechanical data for the endothelium do not exist, we assumed a Fung-type exponential energy function, which is known to characterize many biological tissues (10), and estimated the respective material parameters. The one-dimensional (1D) modulus for the endothelium (36) was compared with that of vessel wall (39) and found to be 1,000 times softer. The issue was to transform the 1D experimental data reported for the endothelium to a 3D constitutive relation. The material parameter selection for the estimated 3D endothelial constitutive relation preserves both the linearity (determined by b
) and the magnitude of the modulus (determined by Cb
) obtained from the shear loading endothelial experiments. Although we selected the parameter C as 1% and b
(
= 1 ,... , 7) as 10% of those of the intact artery, C of 10% and b
of 1% would lead to the same conclusion. Similarly, other combinations of C and b
lead to a similar conclusion regarding the shear deformation, as long as we maintain the observed linearity for the endothelium (decrease b
) and modulus (Cb
of 0.1%).
Here the media and adventitia have been regarded as one layer because their mechanical properties are on the same order. The conclusion that shear strain in intima may be significant will not change if media and adventitia layers are treated differently. We assumed that the intima is equal to 2% of the wall thickness, which is reasonable in young and healthy blood vessels. In diseased or aged arteries, the intima becomes thicker and stiffer (17), and its role becomes more important in tensile support (28, 29). The thickened and hardened intima will change the strain distribution, which may affect mechanotransduction, if strain is the stimulus for homeostasis.
The arterial wall has been modeled as an elastic material, and hence the stress and strain distributions are independent of the loading path. If the viscoelastic behavior of blood vessels is considered, the loading sequence must be taken into account. A 3D finite-element model for the viscoelastic deformation of endothelial cells has been considered in the literature (21) and is beyond the scope of the present study. Other issues, such as mechanical property of perivascular tethering, non-Newtonian flow, complex geometry, and heterogeneity of the vessel wall, will make analytical modeling impractical. Some of these effects can be implemented in finite-element method considering fluid-solid interaction, constraint of surrounding tissues, etc. (20, 41).
It should be noted that our modeling has been conducted on the macroscopic level based on continuum mechanics. Therefore, stress and strain in the arterial wall must be interpreted as the average values of various microstructural constituents. This is different from the works of Huang et al. (18) and Karcher et al. (21), as well as other investigations, where details of stresses and strains in the individual cells with realistic shapes and geometries have been considered. Despite the lack of microstructural details, the present model captures the essence of the question raised.
Chaudhry et al. (3) reported that residual stress and strain reduce the circumferential stress gradient in oscillating arterial wall. Cinthio et al. (6) studied the longitudinal movement of arterial wall in cardiac cycle and concluded that the shear strain in artery may be substantial. Based on a general formulation, including dynamics and smooth muscle active terms developed by Humphrey and Na (19), the current model can be extended to investigate arterial behavior involving pulsatile blood flow, smooth muscle cell activity, and other physiological features.
Summary and significance of study. A computational model that takes into account coupled pressure and shear loading was used to analyze stress and strain distributions in the porcine coronary artery. It is found that the flow-induced shear strain in intima can be as large as normal strains caused by internal blood pressure. This finding may have significant implications for interpreting vascular behaviors, such as vasoactivity due to change of blood pressure and flow, arterial wall remodeling, and pathophysiology (atherosclerosis). It may shed light on the mechanical stimulation mechanisms in the endothelial cell.
It is suggested that the intima layer should be considered separately to investigate the role of flow-induced shear, since the deformation in the intima may be overshadowed by the media and adventitia layers, which have much larger elastic moduli (flow-induced shear strain was found to be Erz
–3 x 10–5 in a one-layer model with intact properties listed in Table 1, or a classic intima-media and adventitia two-layer model as mentioned earlier). From the numerical results of the intima and media-adventitia two-layer model, we hypothesize that shear strain may be a stimulus for the endothelium. In this regard, strain may be a good measure to correlate with biochemical response of blood vessels. This hypothesis needs to be explored in future experimental studies.
| APPENDIX |
|---|
|
|
|---|
, z) are related to those in the zero-stress state (R,
, Z) by
![]() | (A1) |
=
/(
–
),
z is the stretch ratio in the axial direction, and w is the relative displacement in the z direction. The deformation gradient matrix with respect to the zero-stress state is
![]() | (A2) |
, z, and
![]() | (A3) |
jk), where l = r,
, z denotes the dummy index for Einstein summation, is
![]() | (A4) |
![]() | (A5) |
![]() | (A6) |
![]() | (A7) |
![]() | (A8) |
![]() | (A9) |
![]() | (A10) |
![]() | (A11) |
Material parameters for Fung model.
The parameters in the normal directions for an intact vessel were taken as the average values reported by Wang et al. (39). For the parameter b7, we assumed that the shear modulus in the rz plane is similar to that in the
z plane. The shear modulus 160 kPa in the
z plane (27) for
z = 1.4, and circumferential stress 90 kPa (
1.6) were used to estimate b7 = 2.89 for Fung model.
The mechanical behavior of endothelium, which is significantly different when the cells are viewed as fluidlike or solidlike (9), is essential to predict the correct stress and strain in intima. Caille et al. (2) adopted an isotropic material model (Mooney-Rivlin law) to fit the experimental force-deformation curves of bovine aortic endothelial cells and determined that their elastic moduli are on the order of 0.5 kPa for spread and round cells and 5.0 kPa for cells with nuclei. Sato et al. (36) reported the linear viscoelastic properties of cultured endothelial cells under shear stress. They estimated that Young's modulus of porcine aortic endothelial cells is on the order of 0.1 kPa. Considering that the porcine coronary artery elastic modulus is
170 kPa (26), it is rationalized that the elastic stiffness of porcine endothelium is
1,000 times smaller than that of the overall vessel. Due to the lack of data, we assumed that the Fung model applies to the intima layer. The properties for intima (basal lamina and subendothelial layer are assumed to have the same property as endothelium), as given in Table 1, were selected to ensure this modulus criterion (modulus is proportional to the product of C and b
).
It should be noted that shear modulus in the rz plane may be different from that in the
z plane. In quantitative modeling, accurate experimental results are desired. Since the endothelial cells are soft and thin, it is difficult to obtain their 3D constitutive equation and material parameters. Hence, direct measurement of stress-strain relationship of intima remains a technical challenge.
| GRANTS |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
| REFERENCES |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |