Filipovic Nenad, Ghimire Kedar, Saveljic Igor, Milosevic Zarko, Ruegg Curzio
a Department for Applied Mechanics and Automatic Control, Faculty of Engineering, Center for Bioengineering , University of Kragujevac , Sestre Janjic 6, Kragujevac 34000 , Serbia.
b Department of Medicine, Faculty of Science , University of Fribourg , Fribourg , Switzerland.
Comput Methods Biomech Biomed Engin. 2016;19(6):581-90. doi: 10.1080/10255842.2015.1051973. Epub 2015 Jun 22.
Vascular endothelial cells are continuously exposed to hemodynamic shear stress. Intensity and type of shear stress are highly relevant to vascular physiology and pathology. Here, we modeled shear stress distribution in a tissue culture well (R = 17.5 mm, fill volume 2 ml) under orbital translation using computational fluid dynamics with the finite element method. Free surface distribution, wall shear stress, inclination angle, drag force, and oscillatory index on the bottom surface were modeled. Obtained results predict nonuniform shear stress distribution during cycle, with higher oscillatory shear index, higher drag force values, higher circular component, and larger inclination angle of the shear stress at the periphery of the well compared with the center of the well. The oscillatory index, inclination angle, and drag force are new quantitative parameters modeled in this system, which provide a better understanding of the hydrodynamic conditions experienced and reflect the pulsatile character of blood flow in vivo. Validation experiments revealed that endothelial cells at the well periphery aligned under flow and increased Kruppel-like Factor 4 (KLF-4), cyclooxygenase-2 (COX-2) expression and endothelial nitric oxide synthase (eNOS) phosphorylation. In contrast, endothelial cells at the center of the well did not show clear directional alignment, did not induce the expression of KLF-4 and COX-2 nor increased eNOS phosphorylation. In conclusion, this improved computational modeling predicts that the orbital shaker model generates different hydrodynamic conditions at the periphery versus the center of the well eliciting divergent endothelial cell responses. The possibility of generating different hydrodynamic conditions in the same well makes this model highly attractive to study responses of distinct regions of the same endothelial monolayer to different types of shear stresses thereby better reflecting in vivo conditions.
血管内皮细胞持续暴露于血流动力学剪切应力之下。剪切应力的强度和类型与血管生理及病理密切相关。在此,我们采用有限元法通过计算流体动力学对组织培养孔(半径R = 17.5毫米,填充体积2毫升)在轨道平移下的剪切应力分布进行建模。对自由表面分布、壁面剪切应力、倾斜角度、阻力以及底面的振荡指数进行了建模。所得结果预测,在循环过程中剪切应力分布不均匀,与孔中心相比,孔周边的振荡剪切指数更高、阻力值更高、圆周分量更大且剪切应力的倾斜角度更大。振荡指数、倾斜角度和阻力是该系统中建模的新定量参数,它们有助于更好地理解所经历的流体动力学条件,并反映体内血流的脉动特性。验证实验表明,孔周边的内皮细胞在流动下排列成行,并增加了 Kruppel样因子4(KLF - 4)、环氧化酶 - 2(COX - 2)的表达以及内皮型一氧化氮合酶(eNOS)的磷酸化。相比之下,孔中心的内皮细胞未显示出明显的定向排列,未诱导KLF - 4和COX - 2的表达,也未增加eNOS的磷酸化。总之,这种改进的计算模型预测,轨道振荡器模型在孔周边与中心产生不同的流体动力学条件,从而引发不同的内皮细胞反应。在同一孔中产生不同流体动力学条件的可能性使得该模型对于研究同一内皮单层不同区域对不同类型剪切应力的反应极具吸引力,从而能更好地反映体内情况。