A mechanical model for morphological response of endothelial cells under combined wall shear stress and cyclic stretch loadings

The shape and morphology of endothelial cells (ECs) lining the blood vessels are a good indicator for atheroprone and atheroprotected sites. ECs of blood vessels experience both wall shear stress (WSS) and cyclic stretch (CS). These mechanical stimuli influence the shape and morphology of ECs. A few models have been proposed for predicting the morphology of ECs under WSS or CS. In the present study, a mathematical cell population model is developed to simulate the morphology of ECs under combined WSS and CS conditions. The model considers the cytoskeletal filaments, cell–cell interactions, and cell–extracellular matrix interactions. In addition, the reorientation and polymerization of microfilaments are implemented in the model. The simulations are performed for different conditions: without mechanical stimuli, under pure WSS, under pure CS, and under combined WSS and CS. The results are represented as shape and morphology of ECs, shape index, and angle of orientation. The model is validated qualitatively and quantitatively with several experimental studies, and good agreement with experimental studies is achieved. To the best of our knowledge, it is the first model for predicting the morphology of ECs under combined WSS and CS condition. The model can be used to indicate the atheroprone regions of a patient’s artery.     

Effects of biaxial oscillatory shear stress on endothelial cell proliferation and morphology

Wall shear stress (WSS) on anchored cells affects their responses, including cell proliferation and morphology. In this study, the effects of the directionality of pulsatile WSS on endothelial cell proliferation and morphology were investigated for cells grown in a Petri dish orbiting on a shaker platform. Time and location dependent WSS was determined by computational fluid dynamics (CFD). At low orbital speed (50 rpm), WSS was shown to be uniform (0-1 dyne/cm(2)) across the bottom of the dish, while at higher orbital speed (100 and 150 rpm), WSS remained fairly uniform near the center and fluctuated significantly (0-9 dyne/cm(2)) near the side walls of the dish. Since WSS on the bottom of the dish is two-dimensional, a new directional oscillatory shear index (DOSI) was developed to quantify the directionality of oscillating shear. DOSI approached zero for biaxial oscillatory shear of equal magnitudes near the center and approached one for uniaxial pulsatile shear near the wall, where large tangential WSS dominated a much smaller radial component. Near the center (low DOSI), more, smaller and less elongated cells grew, whereas larger cells with greater elongation were observed in the more uniaxial oscillatory shear (high DOSI) near the periphery of the dish. Further, cells aligned with the direction of the largest component of shear but were randomly oriented in low magnitude biaxial shear. Statistical analyses of the individual and interacting effects of multiple factors (DOSI, shear magnitudes and orbital speeds) showed that DOSI significantly affected all the responses, indicating that directionality is an important determinant of cellular responses.     

Tetrahedral vs. polyhedral mesh size evaluation on flow velocity and wall shear stress for cerebral hemodynamic simulation

Haemodynamic factors, in particular wall shear stresses (WSSs) may have significant impact on growth and rupture of cerebral aneurysms. Without a means to measure WSS reliably in vivo, computational fluid dynamic (CFD) simulations are frequently employed to visualise and quantify blood flow from patient-specific computational models. With increasing interest in integrating these CFD simulations into pretreatment planning, a better understanding of the validity of the calculations in respect to computation parameters such as volume element type, mesh size and mesh composition is needed. In this study, CFD results for the two most common aneurysm types (saccular and terminal) are compared for polyhedral- vs. tetrahedral-based meshes and discussed regarding future clinical applications. For this purpose, a set of models were constructed for each aneurysm with spatially varying surface and volume mesh configurations (mesh size range: 5119-258, 481 volume elements). WSS distribution on the model wall and point-based velocity measurements were compared for each configuration model. Our results indicate a benefit of polyhedral meshes in respect to convergence speed and more homogeneous WSS patterns. Computational variations of WSS values and blood velocities are between 0.84 and 6.3% from the most simple mesh (tetrahedral elements only) and the most advanced mesh design investigated (polyhedral mesh with boundary layer).     

Fluid mechanics of the zebrafish embryonic heart trabeculation

Embryonic heart development is a mechanosensitive process, where specific fluid forces are needed for the correct development, and abnormal mechanical stimuli can lead to malformations. It is thus important to understand the nature of embryonic heart fluid forces. However, the fluid dynamical behaviour close to the embryonic endocardial surface is very sensitive to the geometry and motion dynamics of fine-scale cardiac trabecular surface structures. Here, we conducted image-based computational fluid dynamics (CFD) simulations to quantify the fluid mechanics associated with the zebrafish embryonic heart trabeculae. To capture trabecular geometric and motion details, we used a fish line that expresses fluorescence at the endocardial cell membrane, and high resolution 3D confocal microscopy. Our endocardial wall shear stress (WSS) results were found to exceed those reported in existing literature, which were estimated using myocardial rather than endocardial boundaries. By conducting simulations of single intra-trabecular spaces under varied scenarios, where the translational or deformational motions (caused by contraction) were removed, we found that a squeeze flow effect was responsible for most of the WSS magnitude in the intra-trabecular spaces, rather than the shear interaction with the flow in the main ventricular chamber. We found that trabecular structures were responsible for the high spatial variability of the magnitude and oscillatory nature of WSS, and for reducing the endocardial deformational burden. We further found cells attached to the endocardium within the intra-trabecular spaces, which were likely embryonic hemogenic cells, whose presence increased endocardial WSS. Overall, our results suggested that a complex multi-component consideration of both anatomic features and motion dynamics were needed to quantify the trabeculated embryonic heart fluid mechanics.     

Morphological differences between guinea pig aortic and venous endothelial cells in situ

Endothelial cells (ECs) respond to fluid shear stress. They reveal shear stress related morphological changes in both their cell shape and cytoskeletal organization. Little is known about the cytoskeletal organization of ECs in situ. We studied, together with the living ultrasound high resolution imaging system, the distribution of stress fibers (SFs), certain focal adhesion (FA) and signal transduction associated proteins in guinea pig aortic and venous ECs. Although SFs present in the basal portion of venous ECs ran along the direction of the blood flow, their size was smaller and their number was fewer than those of aortic ECs. Venous ECs were elongated to the direction of flow than in aortic ECs exposed over normal shear stress (SS). Since fluid SS in the vein is low, a sustained and uni-directional low SS over a long period might thus cause these structural features observed in venous ECs.     

Epigenetic response of endothelial cells to different wall shear stress magnitudes: A report of new mechano-miRNAs

Endothelial cells (ECs) respond to flow stress via a variety of mechanisms, leading to various intracellular responses that can modulate the vessel wall and lead to diseases if the flow is disturbed. Mechano-microRNAs (miRNAs) are a subset of miRNAs in the ECs that are flow responsive. Mechano-miRNAs were shown to be related to atherosclerosis pathophysiology, and a number of them were identified as pathologic. Here, we exposed human carotid ECs to different wall shear stresses (WSS), high and low, and evaluated the response of miRNAs by microarray and quantitative polymerase chain reaction analysis. We discovered five new mechano-miRNAs that were not reported in that context previously to the best of our knowledge. Moreover, functional pathway analysis revealed that under low WSS conditions, several pathways regulating apoptosis are affected. In addition, KLF2 and KLF4, known atheroprotective genes, were downregulated under low WSS and upregulated under high WSS. KLF2 and VCAM1, both angiogenic, were upregulated under high WSS. NOS3, which is vascular protective, was also upregulated with higher WSS. On the contrary, ICAM-1 and E-selectin, both atherogenic and proinflammatory, were upregulated with high WSS. Collectively, the epigenetic landscape with the gene expression analysis reveals that low WSS is associated with a proapoptotic state, while high WSS is associated with a proliferative and proinflammatory state.     

Numerical simulation of pulsatile flow in a compliant curved tube model of a coronary artery

The endothelial cells (ECs) lining a blood vessel wall are exposed to both the wall shear stress (WSS) of blood flow and the circumferential strain (CS) of pulsing artery wall motion. These two forces and their interaction are believed to play a role in determining remodeling of the vessel wall and development of arterial disease (atherosclerosis). This study focused on the WSS and CS dynamic behavior in a compliant model of a coronary artery taking into account the curvature of the bending artery and physiological radial wall motion. A three-dimensional finite element model with transient flow and moving boundaries was set up to simulate pulsatile flow with physiological pressure and flow wave forms characteristic of the coronary arteries. The characteristic coronary artery curvature and flow conditions applied to the simulation were: aspect ratio (lambda) = 10, diameter variation (DV) = 6 percent, mean Reynolds number (Re) = 150, and unsteadiness parameter (alpha) = 3. The results show that mean WSS is about 50 percent lower on the inside wall than the outside wall while WSS oscillation is stronger on the inside wall. The stress phase angle (SPA) between CS and WSS, which characterizes the dynamics of the mechanical force pattern applied to the endothelial cell layer, shows that CS and WSS are more out of phase in the coronaries than in any other region of the circulation (-220 deg on the outside wall, -250 deg on the inside wall). This suggests that in addition to WSS, SPA may play a role in localization of coronary atherosclerosis.     

MRI and CFD studies of pulsatile flow in healthy and stenosed carotid bifurcation models

Pulsatile flow was studied in physiologically realistic models of a normal and a moderately stenosed (30% diameter reduction) human carotid bifurcation. Time-resolved velocity measurements were made using magnetic resonance imaging, from which wall shear stress (WSS) vectors were calculated. Velocity measurements in the inflow and outflow regions were also used as boundary conditions for a computational fluid dynamics (CFD) model. Experimental flow patterns and derived WSS vectors were compared qualitatively with the corresponding CFD predictions. In the stenosed phantom, flow in the bulb region of the "internal carotid artery" was concentrated along the outer wall, with a region of low and recirculating flow near the inner wall. In the normal phantom, the converse was found, with a low flow region near the outer wall of the bulb. Time-averaged WSS and oscillatory shear index were also markedly different for the two phantoms.     

Approximating hemodynamics of cerebral aneurysms with steady flow simulations

Computational fluid dynamics (CFD) simulations can be employed to gain a better understanding of hemodynamics in cerebral aneurysms and improve diagnosis and treatment. However, introduction of CFD techniques into clinical practice would require faster simulation times. The aim of this study was to evaluate the use of computationally inexpensive steady flow simulations to approximate the aneurysm's wall shear stress (WSS) field. Two experiments were conducted. Experiment 1 compared for two cases the time-averaged (TA), peak systole (PS) and end diastole (ED) WSS field between steady and pulsatile flow simulations. The flow rate waveform imposed at the inlet was varied to account for variations in heart rate, pulsatility index, and TA flow rate. Consistently across all flow rate waveforms, steady flow simulations accurately approximated the TA, but not the PS and ED, WSS field. Following up on experiment 1, experiment 2 tested the result for the TA WSS field in a larger population of 20 cases covering a wide range of aneurysm volumes and shapes. Steady flow simulations approximated the space-averaged WSS with a mean error of 4.3%. WSS fields were locally compared by calculating the absolute error per node of the surface mesh. The coefficient of variation of the root-mean-square error over these nodes was on average 7.1%. In conclusion, steady flow simulations can accurately approximate the TA WSS field of an aneurysm. The fast computation time of 6 min per simulation (on 64 processors) could help facilitate the introduction of CFD into clinical practice.     

A model for studying the effect of shear stress on interactions between vascular endothelial cells and smooth muscle cells

Vascular endothelial cells (ECs) are constantly subjected to blood flow-induced shear stress and the influences of neighboring smooth muscle cells (SMCs). In the present study, a coculture flow system was developed to study the effect of shear stress on EC-SMC interactions. ECs and SMCs were separated by a porous membrane with only the EC side subjected to the flow condition. When ECs were exposed to a shear stress of 12 dynes/cm2 for 24 h, the cocultured SMCs tended to orient perpendicularly to the flow direction. This perpendicular orientation of the cocultured SMCs to flow direction was not observed when ECs were exposed to a shear stress of 2 dynes/cm2. Under the static condition, long and parallel actin bundles were observed in the central regions of the cocultured SMCs, whereas the actin filaments localized mainly at the periphery of the cocultured ECs. After 24 h of flow application, the cocultured ECs displayed very long, well-organized, parallel actin stress fibers aligned with the flow direction in the central regions of the cells. Immunostaining of platelet endothelial cell adhesion molecule-1 confirmed the elongation and alignment of the cocultured ECs with the flow direction. Coculture with SMCs under static condition induced EC gene expressions of growth-related oncogene-alpha and monocyte chemotactic protein-1, and shear stress was found to abolish these SMC-induced gene expressions. Our results suggest that shear stress may serve as a down-regulator for the pathophysiologically relevant gene expression in ECs cocultured with SMCs.     

Fluid shear stress modulates endothelial cell invasion into three-dimensional collagen matrices

Endothelial cells are subjected to biochemical and mechanical stimuli, which regulate their angiogenic potential. We determined the synergistic effects of sphingosine-1-phosphate (S1P) and fluid wall shear stress (WSS) on a previously established model of human umbilical vein endothelial cell invasion into three-dimensional collagen matrices. Collagen matrices were incorporated into a parallel-plate flow chamber to apply controlled WSS to the surface of endothelial monolayers over a period of 24 h. Cell invasion required the presence of S1P, with the effects of S1P being enhanced by shear stress to an extent comparable with S1P combined with angiogenic growth factor stimulation. The number of invading cells depended on the magnitude of shear stress, with a maximal induction at a shear stress of approximately 5 dyn/cm2, whereas the invasion distance was proportional to the magnitude of shear stress. The enhancement of invasion by 5.3 dyn/cm2 shear stress coincided with elevated phosphorylation of Akt and matrix metalloproteinase (MMP)-2 activation. Furthermore, invasion induced by the combined application of WSS and S1P was attenuated by inhibitors of MMPs (GM6001) and the phosphatidylinositol 3-kinase/Akt signaling pathway (wortmannin). These results provide evidence that shear stress is a positive modulator of S1P-induced endothelial cell invasion into collagen matrices through enhanced Akt and MMP-2 activation.     

Optimization of cardiovascular stent design using computational fluid dynamics

Coronary stent design affects the spatial distribution of wall shear stress (WSS), which can influence the progression of endothelialization, neointimal hyperplasia, and restenosis. Previous computational fluid dynamics (CFD) studies have only examined a small number of possible geometries to identify stent designs that reduce alterations in near-wall hemodynamics. Based on a previously described framework for optimizing cardiovascular geometries, we developed a methodology that couples CFD and three-dimensional shape-optimization for use in stent design. The optimization procedure was fully-automated, such that solid model construction, anisotropic mesh generation, CFD simulation, and WSS quantification did not require user intervention. We applied the method to determine the optimal number of circumferentially repeating stent cells (N(C)) for slotted-tube stents with various diameters and intrastrut areas. Optimal stent designs were defined as those minimizing the area of low intrastrut time-averaged WSS. Interestingly, we determined that the optimal value of N(C) was dependent on the intrastrut angle with respect to the primary flow direction. Further investigation indicated that stent designs with an intrastrut angle of approximately 40 deg minimized the area of low time-averaged WSS regardless of vessel size or intrastrut area. Future application of this optimization method to commercially available stent designs may lead to stents with superior hemodynamic performance and the potential for improved clinical outcomes.     

Blood flow in stented arteries: a parametric comparison of strut design patterns in three dimensions

BACKGROUND: Restenosis after stent implantation varies with stent design. Alterations in secondary flow patterns and wall shear stress (WSS) can modulate intimal hyperplasia via their effects on platelet and inflammatory cell transport toward the wall, as well as direct effects on the endothelium. METHOD OF APPROACH: Detailed flow characteristics were compared by estimating the WSS in the near-strut region of realistic stent designs using three-dimensional computational fluid dynamics (CFD), under pulsatile high and low flow conditions. The stent geometry employed was characterized by three geometric parameters (axial strut pitch, strut amplitude, and radius of curvature), and by the presence or lack of the longitudinal connector. RESULTS: Stagnation regions were localized around stent struts. The regions of low WSS are larger distal to the strut. Under low flow conditions, the percentage restoration of mean axial WSS between struts was lower than that for the high flow by 10-12%. The largest mean transverse shear stresses were 30-50% of the largest mean axial shear stresses. The percentage restoration in WSS in the models without the longitudinal connector was as much as 11% larger than with the connector The mean axial WSS restoration between the struts was larger for the stent model with larger interstrut spacing. CONCLUSION: The results indicate that stent design is crucial in determining the fluid mechanical environment in an artery. The sensitivity of flow characteristics to strut configuration could be partially responsible for the dependence of restenosis on stent design. From a fluid dynamics point of view, interstrut spacing should be larger in order to restore the disturbed flow; struts should be oriented to the flow direction in order to reduce the area of flow recirculation. Longitudinal connectors should be used only as necessary, and should be parallel to the axis. These results could guide future stent designs toward reducing restenosis.     

Flow rates in perfusion bioreactors to maximise mineralisation in bone tissue engineering in vitro

In bone tissue engineering experiments, fluid-induced shear stress is able to stimulate cells to produce mineralised extracellular matrix (ECM). The application of shear stress on seeded cells can for example be achieved through bioreactors that perfuse medium through porous scaffolds. The generated mechanical environment (i.e. wall shear stress: WSS) within the scaffolds is complex due to the complexity of scaffold geometry. This complexity has so far prevented setting an optimal loading (i.e. flow rate) of the bioreactor to achieve an optimal distribution of WSS for stimulating cells to produce mineralised ECM. In this study, we demonstrate an approach combining computational fluid dynamics (CFD) and mechano-regulation theory to optimise flow rates of a perfusion bioreactor and various scaffold geometries (i.e. pore shape, porosity and pore diameter) in order to maximise shear stress induced mineralisation. The optimal flow rates, under which the highest fraction of scaffold surface area is subjected to a wall shear stress that induces mineralisation, are mainly dependent on the scaffold geometries. Nevertheless, the variation range of such optimal flow rates are within 0.5–5 mL/min (or in terms of fluid velocity: 0.166–1.66 mm/s), among different scaffolds. This approach can facilitate the determination of scaffold-dependent flow rates for bone tissue engineering experiments in vitro, avoiding performing a series of trial and error experiments.     

Influence of surface model extraction parameter on computational fluid dynamics modeling of cerebral aneurysms

Threshold image intensity for reconstructing patient-specific vascular models is generally determined subjectively. We assessed the effects of threshold image intensity differences on computational fluid dynamics (CFD) using a simple method of threshold determination. This study included 11 consecutive patients with internal carotid artery aneurysms collected retrospectively between April 2009 and March 2010. In 3-dimensional rotational angiography image data, we set a line probe across the coronal cross-section of the parent internal carotid artery, and calculated a profile curve of the image intensity along this line. We employed the threshold coefficient (C(thre)) value in this profile curve, in order to determine the threshold image intensity objectively. We assessed the effects of C(thre) value differences on vascular model configuration and the wall shear stress (WSS) distribution of the aneurysm. The threshold image intensity increased as the C(thre) value increased. The frequency of manual editing increased as the C(thre) value decreased, while disconnection of the posterior communicating artery occurred more frequently as the C(thre) value increased. The volume of the vascular model decreased and WSS increased according to the C(thre) value increase. The pattern of WSS distribution changed remarkably in one case. Threshold image intensity differences can produce profound effects on CFD. Our results suggest the uniform setting of C(thre) value is important for objective CFD.     

Flow patterns and wall shear stress distribution in human internal carotid arteries: the geometric effect on the risk for stenoses

It has been widely observed that atherosclerotic stenosis occurs at sites with complex hemodynamics, such as arteries with high curvature or bifurcations. These regions usually have very low or highly oscillatory wall shear stress (WSS). In the present study, 3D sinusoidally pulsatile blood flow through the models of internal carotid artery (ICA) with different geometries was investigated with computational simulation. Three preferred sites of stenoses were found along the carotid siphon with low and highly oscillatory WSS. The risk for stenoses at these sites was scaled with the values of time-averaged WSS and oscillating shear index (OSI). The local risk for stenoses at every preferred site of stenoses was found different between 3 types of ICA, indicating that the geometry of the blood vessel plays significant roles in the atherogenesis. Specifically, the large curvature and planarity of the vessel were found to increase the risk for stenoses, because they tend to lower WSS and elevate OSI. Therefore, the geometric study makes it possible to estimate the stenosis location in the ICA siphon as long as the shape of ICA was measured.     

Errors in power-law estimations of inflow rates for intracranial aneurysm CFD

Patient-specific inflow rates are rarely available for computational fluid dynamics (CFD) studies of intracranial aneurysms. Instead, inflow rates are often estimated from parent artery diameters via power laws, i.e. Q ∝ Dⁿ, reflecting adaptation of conduit arteries to demanded flow. The present study aimed to validate the accuracy of these power laws. Internal carotid artery (ICA) flow rates were measured from 25 ICA aneurysm patients via 2D phase contrast MRI. ICA diameters, derived from 3D segmentation of rotational angiograms, were used to estimate inflow rates via power laws from the aneurysm CFD literature assuming the same inlet wall shear stress (WSS) (n = 3), velocity (n = 2) or flow rate (n = 0) for all cases. To illustrate the potential impact of errors in flow rate estimates, pulsatile CFD was carried out for four cases having large errors for at least one power law. Flow rates estimated by n = 3 and n = 0 power laws had significant (p < 0.01) mean biases of −22% to +32%, respectively, but with individual errors ranging from −78% to +120%. The n = 2 power law had no significant bias, but had non-negligible individual errors of −58% to +71%. CFD showed similarly large errors for time-averaged sac WSS; however, these were reduced after normalizing by parent artery WSS. High frequency WSS fluctuations, evident in 2/4 aneurysms, were also sensitive to inflow rate errors. Care should therefore be exercised in the interpretation of aneurysm CFD studies that rely on power law estimates of inflow rates, especially if absolute (vs. normalized) WSS, or WSS instabilities, are of interest.     

Flow pattern analysis in a highly stenotic patient-specific carotid bifurcation model using a turbulence model

The aim of this study is to investigate the blood flow pattern in carotid bifurcation with a high degree of luminal stenosis, combining in vivo magnetic resonance imaging (MRI) and computational fluid dynamics (CFD). A newly developed two-equation transitional model was employed to evaluate wall shear stress (WSS) distribution and pressure drop across the stenosis, which are closely related to plaque vulnerability. A patient with an 80% left carotid stenosis was imaged using high resolution MRI, from which a patient-specific geometry was reconstructed and flow boundary conditions were acquired for CFD simulation. A transitional model was implemented to investigate the flow velocity and WSS distribution in the patient-specific model. The peak time-averaged WSS value of approximately 73 Pa was predicted by the transitional flow model, and the regions of high WSS occurred at the throat of the stenosis. High oscillatory shear index values up to 0.50 were present in a helical flow pattern from the outer wall of the internal carotid artery immediately after the throat. This study shows the potential suitability of a transitional turbulent flow model in capturing the flow phenomena in severely stenosed carotid arteries using patient-specific MRI data and provides the basis for further investigation of the links between haemodynamic variables and plaque vulnerability. It may be useful in the future for risk assessment of patients with carotid disease.     

Computational fluid dynamics in abdominal aorta bifurcation: non-Newtonian versus Newtonian blood flow in a real case study

Hemodynamic in abdominal aorta bifurcation was investigated in a real case using computational fluid dynamics. A Newtonian and non-Newtonian (Walburn-Schneck) viscosity models were compared. The geometrical model was obtained by 3D reconstruction from CT-scan and hemodynamic parameters obtained by laser-Doppler. Blood was assumed incompressible fluid, laminar flow in transient regime and rigid vessel wall. Finite volume-based was used to study the velocity, pressure, wall shear stress (WSS) and viscosity throughout cardiac cycle. Results obtained with Walburn-Schneck’s model, during systole, present lower viscosity due to shear thinning behavior. Furthermore, there is a significant difference between the results obtained by the two models for a specific patient. During the systole, differences are more pronounced and are preferably located in the tortuous regions of the artery. Throughout the cardiac cycle, the WSS amplitude between the systole and diastole is greater for the Walburn-Schneck’s model than for the Newtonian model. However, the average viscosity along the artery is always greater for the non-Newtonian model, except in the systolic peak. The hemodynamic model is crucial to validate results obtained with CFD and to explore clinical potential.