血管彎曲程度對內(nèi)源性一氧化氮分布的影響

趙紅俊 韓金濤 劉 聰 Dov Jaron 喬惠婷#*

1(北京航空航天大學(xué)生物與醫(yī)學(xué)工程學(xué)院，北京 100191)2(北京大學(xué)第三醫(yī)院介入血管外科，北京 100191)3(德雷塞爾大學(xué)生物醫(yī)學(xué)工程、理學(xué)與健康系統(tǒng)學(xué)院，美國 費(fèi)城 PA19104)

血管彎曲程度對內(nèi)源性一氧化氮分布的影響

趙紅俊1韓金濤2劉 聰1Dov Jaron3喬惠婷1#*

1(北京航空航天大學(xué)生物與醫(yī)學(xué)工程學(xué)院，北京 100191)2(北京大學(xué)第三醫(yī)院介入血管外科，北京 100191)3(德雷塞爾大學(xué)生物醫(yī)學(xué)工程、理學(xué)與健康系統(tǒng)學(xué)院，美國 費(fèi)城 PA19104)

一氧化氮(NO)在血流調(diào)控中起著重要作用，目前關(guān)于血管彎曲程度對NO分布的影響尚不明確。建立11個(gè)不同彎曲程度的幾何模型，利用多物理場耦合建模仿真的方法，對內(nèi)源性NO在血液、血管壁及周圍組織中的濃度分布進(jìn)行數(shù)值模擬。數(shù)值結(jié)果表明，在血管彎曲處，NO徑向分布不均，內(nèi)外側(cè)濃度存在差異，且隨著彎曲程度的增加，濃度差值呈先上升后緩慢降低的趨勢。在血管投影長度為定值20 mm 的情況下，當(dāng)彎曲高度取2 mm、對應(yīng)曲率為0.04時(shí)，NO在內(nèi)外側(cè)平滑肌區(qū)域的濃度差值達(dá)到最大，外側(cè)邊界比內(nèi)側(cè)高23.74%。研究揭示彎曲程度對內(nèi)源性NO分布的影響，所產(chǎn)生的NO濃度分布差異可能是導(dǎo)致血管形態(tài)變化以及血管病變的重要因素，為進(jìn)一步探究發(fā)病機(jī)理以及患病風(fēng)險(xiǎn)評估提供一定的理論依據(jù)。

一氧化氮； 彎曲曲率； 多物理場仿真； 平滑肌舒張

引言

內(nèi)源性一氧化氮(Nitric Oxide，NO)在血流調(diào)控過程中起著重要的作用[1- 2]，它可以松弛平滑肌調(diào)節(jié)血管的緊張程度，阻止內(nèi)膜增生，抑制血小板聚集防止動(dòng)脈粥樣硬化的發(fā)生[3- 5]。NO主要生成于內(nèi)皮細(xì)胞，由內(nèi)皮型一氧化氮合成酶(eNOS)產(chǎn)生。所產(chǎn)生的NO以一定濃度擴(kuò)散到平滑肌細(xì)胞附近，引起平滑肌的舒張[6]。研究表明，NO的產(chǎn)生速率受壁面剪切力影響，在擴(kuò)散過程中NO濃度受血流流動(dòng)以及反應(yīng)消耗影響，因此NO代謝是一個(gè)受多因素影響的復(fù)雜過程[7- 8]。

目前，關(guān)于NO分布在心腦血管系統(tǒng)中的相關(guān)研究主要利用的是數(shù)值模擬的方法。Deonikar等建立了圓柱形血管的NO分布模型[9]，Chen等建立了血管分叉處NO的濃度分布的相關(guān)模型[10]，Liu等則研究了血管狹窄處的NO濃度分布的情況[11]。然而動(dòng)脈粥樣硬化除了好發(fā)于動(dòng)脈血管的分叉處、血管狹窄處，還易發(fā)生在血管彎曲處[12]。至今針對彎曲血管NO濃度分布的研究較少，盡管已有學(xué)者已經(jīng)開始針對具有生理彎曲的主動(dòng)脈的NO濃度分布進(jìn)行了相應(yīng)的研究[13]，但是他們的模型中沒有考慮血管彎曲程度對NO濃度分布的影響，大多只考慮了單一血管形態(tài)下NO的濃度分布。Bouissou等通過尸檢發(fā)現(xiàn)血管彎曲程度存在差異[14]，Han等基于TOF- MRA技術(shù)對血管形態(tài)進(jìn)行測量也發(fā)現(xiàn)了血管構(gòu)型是存在差異的[15]，并且Rudic等動(dòng)物實(shí)驗(yàn)數(shù)據(jù)表明，NO與血管形態(tài)差異以及動(dòng)脈粥樣硬化等血管疾病相關(guān)[16]。因此，對于血管不同彎曲程度的NO濃度分布的研究非常重要，這將有助于更好地揭示慢性心腦血管疾病發(fā)病機(jī)理，并為不同彎曲程度的血管的患病風(fēng)險(xiǎn)評估提供一定的理論依據(jù)。

因此本研究利用多物理場耦合建模仿真的方法對不同彎曲程度的彎曲血管腔內(nèi)、血管壁，以及周圍組織中的NO濃度分布進(jìn)行數(shù)值模擬與分析，模擬得出不同彎曲程度下的血流動(dòng)力學(xué)特征和NO濃度分布特征，并初步揭示彎曲程度對內(nèi)源性NO分布的影響規(guī)律。

1 方法

1.1 模型構(gòu)建

選取顱內(nèi)動(dòng)脈粥樣硬化高發(fā)的大腦中動(dòng)脈M1段的血管尺寸。利用SolidWorks?構(gòu)建彎曲程度變化的血管幾何模型,并將幾何導(dǎo)入到有限元分析軟件COMSOL Multiphysics?中，通過劃分棱柱網(wǎng)格和四面體網(wǎng)格對模型進(jìn)行計(jì)算域離散化，然后分別設(shè)置了血管腔內(nèi)的流場以及整個(gè)域的物質(zhì)傳輸場，通過添加邊界條件，最終得到11個(gè)具有不同彎曲程度的模型。

基于MRA的實(shí)際測量數(shù)據(jù)，血管投影長度設(shè)為固定值20 mm，血液半徑設(shè)為1.4 mm。血管壁厚度設(shè)為0.34 mm[13]，組織層厚度為3.26 mm。本研究選用曲率來描述血管的彎曲程度，實(shí)際大腦中動(dòng)脈M1段并非一個(gè)均勻的圓弧，即動(dòng)脈中各點(diǎn)的曲率并不相等，在頂點(diǎn)處曲率值最大，在本研究投影長度為定值時(shí)，彎曲高度與頂點(diǎn)的曲率存在線性關(guān)系。因此，本研究動(dòng)脈在中心頂點(diǎn)曲率值(k)的選取結(jié)合了M1段實(shí)際測量的彎曲高度(Curved Height，CH)，分別建立了CH為 0(模型1)、0.5(模型2)、1.0(模型3)、1.5(模型4)、2.0(模型5)、2.5(模型6)、3.0(模型7)、3.5(模型8)、4.0(模型9)、4.5(模型10)和5.0 mm(模型11)的不同彎曲程度的幾何模型，與上述11個(gè)模型對應(yīng)的最大曲率分別為0、0.01、0.02、0.03、0.04、0.05、0.06、0.07、0.08、0.09和0.10。

如圖1所示，幾何模型分為3層，從血液中心向外依次分為血液層(blood，B)、血管壁層(vascular wall，W)和組織層(tissue，T)。血液層主要體現(xiàn)為血液的流動(dòng)和NO的擴(kuò)散、對流以及相應(yīng)的反應(yīng)，血液中的紅細(xì)胞中的血紅蛋白為主要的NO清除物質(zhì)。血管壁層主要包括了靠近血管腔的內(nèi)皮細(xì)胞(endothelium，EN)以及靠近血管壁外側(cè)的平滑肌細(xì)胞(smooth muscle，SMC)，其中平滑肌舒張程度是NO血管舒張功能的重要體現(xiàn)。此外，由于EN區(qū)域的厚度相對于血管壁層可忽略，本研究將EN區(qū)域簡化為血管壁層與血液層的邊界，簡稱W- B邊界。同時(shí)血管壁層與組織層的邊界(W- T邊界)代表了最外側(cè)SMC區(qū)域。

圖1 幾何模型說明Fig.1 Illustration of the geometrical models

為探究不同彎曲程度對于NO濃度的影響，除了考慮上述不同彎曲程度的幾何模型外，還需要考慮到血液中的流場以及整個(gè)域的物質(zhì)傳輸?shù)鸟詈戏抡婺M。

1.2 流場

在血液層，假設(shè)血液為不可壓縮牛頓流體，血管壁為固壁、無滑移?？刂品匠滩捎枚ǔＡ?、不可壓縮納維- 斯托克斯(Navier- Stokes)方程組，有

(1)

式中：ρ表示密度，kg·m-3;u表示速度，m·s-1;p表示壓力，Pa；μ表示黏度，Pa·s。

入口邊界條件采用層流流入，平均速度0.6 ms-1[17]；選用壓力出口，壓力為15 993 Pa(120 mm·Hg)[18]。

1.3 物質(zhì)傳輸場

在血液層，控制方程使用對流- 擴(kuò)散- 反應(yīng)方程，有

DNOΔcB-u·▽cB-kBcB=0

(2)

式中：DNO為擴(kuò)散系數(shù)，取3.3×10-9m2·s-1[8]；cB為NO在血液層的濃度；kB為NO與紅細(xì)胞中血紅蛋白的反應(yīng)速系數(shù)，取23 s-1[13]。

在血管壁層及組織層，主要體現(xiàn)了NO的擴(kuò)散以及消耗情況，控制方程使用如下形式：

DNOΔci-kici=0

(3)

式中,i=W、T分別為血管壁層和組織層，其中kW=kT=0.01 s-1[13]。

內(nèi)皮細(xì)胞上的NO生成主要通過內(nèi)皮通量不連續(xù)的邊界條件來體現(xiàn)，采用如下方程：

NB·nB-NW·nB=RNOH

(4)

式中,NB·nB為從內(nèi)皮細(xì)胞傳輸?shù)窖鼙谥蠳O的通量，-NW·nB為從內(nèi)皮細(xì)胞傳輸?shù)窖芮恢蠳O的通量，RNO為NO的生成率，H為內(nèi)皮厚度，設(shè)定為2 μm[13]。

此外，邊界條件還設(shè)置了初始濃度為0，邊界濃度連續(xù)。

1.4 流場與物質(zhì)傳輸?shù)鸟詈霞胺抡嬗?jì)算

流場中的壁面剪切力(WSS)直接影響NO的產(chǎn)生率RNO，基于Andrews[19]的實(shí)時(shí)測量結(jié)果，RNO與流場中的WSS呈雙曲線的關(guān)系，有

(5)

式中,RNO,basal為基本生成率，Rmax為最大生產(chǎn)率，單位為nM·s-1,τw為壁面剪切力，a為常數(shù)。

流場的速度對血液層NO的對流過程有著直接的影響，如式(2)中u為流場的空間速度。

模型計(jì)算使用有限元分析軟件COMSOL Multiphysics?。網(wǎng)格主要使用棱柱型網(wǎng)格以及四面體網(wǎng)格，為保證網(wǎng)格的獨(dú)立性，本模型已通過了網(wǎng)格敏感性分析，NO濃度的平均差值已小于1%。本研究中的11個(gè)模型除幾何彎曲程度不同，其余均在相同條件下進(jìn)行的計(jì)算。

本研究方法已得到了驗(yàn)證，利用該方法所使用的COMSOL平臺(tái)，構(gòu)建了與文獻(xiàn)[6]幾何尺寸、控制方程、邊界條件設(shè)置均相同的模型，仿真復(fù)現(xiàn)了相同條件下NO濃度分布情況，復(fù)現(xiàn)的仿真結(jié)果與文獻(xiàn)數(shù)據(jù)一致。

2 結(jié)果

2.1 彎曲程度對流場的影響

流場數(shù)值計(jì)算結(jié)果顯示：在直管中，WSS呈均勻分布，且血管中心徑向截面速度呈對稱分布；在彎曲血管中, WSS分布不均，且血管彎曲中心徑向截面速度分布呈非對稱性，體現(xiàn)為由中心向外側(cè)壁偏離；在彎曲處，隨著血管在中心頂點(diǎn)曲率k的增加，WSS逐漸增加，且高剪切力區(qū)域主要出現(xiàn)在幾何中心的外側(cè)中下游(見圖2(a))；同時(shí)，隨著k的增加，中心截面速度向外側(cè)偏離程度呈現(xiàn)逐步增加趨勢(見圖2(b))。這一結(jié)果表示，不同彎曲程度的血管血流動(dòng)力學(xué)特征不同。

圖2 不同彎曲程度下的流場結(jié)果。(a)壁面剪切力；(b)速度分布Fig.2 The results of hemodynamics in different curvature.(a) Wall shear stress; (b) Velocity

2.2 彎曲對NO徑向濃度分布的影響

NO徑向濃度分布數(shù)值計(jì)算結(jié)果顯示：在直管中，NO徑向分布關(guān)于血管中心對稱；在彎管中，NO徑向分布呈非對稱性，體現(xiàn)為外側(cè)濃度大于內(nèi)側(cè)。這一結(jié)果表示，彎曲與否對NO徑向濃度分布存在影響。

同時(shí)，不論是直管還是彎管，NO均在W- B邊界濃度達(dá)到峰值，產(chǎn)生于內(nèi)皮細(xì)胞的NO在向內(nèi)擴(kuò)散到血管腔內(nèi)的過程中由于受到紅細(xì)胞的快速清除作用及血液的流動(dòng)作用，NO的濃度急速下降；NO由內(nèi)皮細(xì)胞向外擴(kuò)散到血管壁及周圍組織的過程中，濃度呈緩慢下降的趨勢(見圖3)，NO濃度分布曲線的該特點(diǎn)與文獻(xiàn)[6,9]相符。

圖3 NO徑向分布特點(diǎn)。(a)整體分布；(b)外側(cè)血管壁局部NO濃度分布放大；(c)內(nèi)側(cè)血管壁局部NO濃度分布放大Fig.3 The radical distributions of NO(a) The whole image; (b)The local concentration distribution of NO in outervascular walls; (c)The local concentration distribution of NO in innervascular walls

2.3 彎曲程度對NO濃度分布的影響

物質(zhì)傳輸?shù)臄?shù)值計(jì)算結(jié)果顯示：在軸向方向，NO濃度最高值出現(xiàn)在W- B邊界，即EN區(qū)域，同時(shí)除直管(k=0)，NO的高濃度區(qū)主要出現(xiàn)在外側(cè)中下游及內(nèi)側(cè)上游(見圖4(b))；在徑向方向，NO在平滑肌細(xì)胞所在的血管壁層的濃度均從W- B邊界到W- T邊界逐漸減少，同時(shí)直管NO濃度關(guān)于中心對稱，在曲率為0.03的模型中，NO外側(cè)血管壁層的濃度略高于內(nèi)側(cè)，在曲率為0.06和0.09的模型中，NO在外側(cè)血管壁層的濃度明顯高于內(nèi)側(cè)(見圖4(c))。這一結(jié)果表示彎曲程度對NO濃度在軸向及徑向分布均存在影響。

圖4 不同彎曲程度下的NO濃度分布。(a)內(nèi)外側(cè)；(b)軸向方向；(c)徑向方向Fig.4 The distributions of NO in different curvature. (a) Inner and outer diagram; (b)Axial direction; (c) Radical direction

此外，表1顯示了11個(gè)彎曲程度不同的模型中NO在W- B邊界的分布情況。在模型1中，NO在內(nèi)外側(cè)W- B邊界的濃度相同；隨著彎曲程度的增加，內(nèi)側(cè)NO的濃度呈現(xiàn)先降低后逐漸增加的趨勢，而外側(cè)NO濃度呈先增加后趨于平緩的趨勢；同時(shí)外側(cè)W- B邊界NO的濃度始終大于內(nèi)側(cè)，差值在曲率為0.04時(shí)達(dá)到最大，外側(cè)比內(nèi)側(cè)高28.02%。表2顯示了11個(gè)彎曲程度不同的模型中NO在W- T邊界的分布情況。在模型1中，NO在內(nèi)外側(cè)W- T邊界的濃度相同；隨著彎曲程度的增加，內(nèi)側(cè)NO的濃度呈現(xiàn)先降低后逐漸增加的趨勢，而外側(cè)NO濃度呈先增加后緩慢降低的趨勢；外側(cè)W- T邊界的濃度始終大于內(nèi)側(cè)濃度，差值同樣是在曲率為0.04時(shí)達(dá)到最大，外側(cè)比內(nèi)側(cè)高23.74%。這一結(jié)果表示，彎曲程度對NO在內(nèi)外兩側(cè)濃度及濃度差值均存在影響。

表1 模型1～11 W- B邊界NO濃度分布Tab.1 The distributions of NO in W- B boundary

表2 模型1～11 W- T邊界NO濃度分布Tab.2 The distributions of NO in W- T boundary

3 討論

血管彎曲處為動(dòng)脈粥樣硬化的好發(fā)部位，除血流動(dòng)力學(xué)因素外，彎曲處的局部物質(zhì)傳輸在血管疾病的發(fā)生過程中亦產(chǎn)生重要影響，如內(nèi)源性NO，其濃度分布是血流調(diào)控作用的關(guān)鍵因素。本研究創(chuàng)新性在于不但對彎曲血管中血流動(dòng)力學(xué)因素進(jìn)行了研究，更重要的是探討了血管彎曲程度對局部NO濃度分布的影響。

從流場的結(jié)果可以看出，不同彎曲程度的血管，其血流動(dòng)力學(xué)特征不同。隨著彎曲程度的增加，筆者發(fā)現(xiàn)除模型1外，其他模型速度均偏向外側(cè)并存在二次流現(xiàn)象，這與國外James等[20]以及國內(nèi)研究者喬愛科[21]等關(guān)于彎曲血管血流動(dòng)力學(xué)的研究結(jié)果一致。這就意味著僅考慮血流速度影響，NO在血液層的濃度分布也是十分復(fù)雜的；同時(shí)，隨著彎曲程度的變化，WSS分布不均，這種WSS的不均勻分布會(huì)引起NO的產(chǎn)生量不同。因此，單一的物理場已無法詮釋復(fù)雜的NO濃度分布問題，故采用多物理場耦合是必要且更接近真實(shí)的。

已有研究表明，動(dòng)脈粥樣硬化斑塊形成的標(biāo)志為依賴于內(nèi)皮細(xì)胞的NO的相關(guān)血管舒張功能受損，而內(nèi)皮細(xì)胞功能紊亂發(fā)生之前常表現(xiàn)為血管壁結(jié)構(gòu)的改變[5]。Cheng等的研究也指出，NO在內(nèi)外側(cè)血管壁中起著重要的維持血管形態(tài)的作用[22]。因此本研究關(guān)于不同彎曲程度對NO濃度分布影響的研究，尤其是針對血管壁層的研究就相當(dāng)重要。本研究觀察到在彎曲處，血管壁內(nèi)外兩側(cè)NO濃度存在差異，且外側(cè)的濃度總是高于內(nèi)側(cè)，這與文獻(xiàn)中主動(dòng)脈彎曲處外側(cè)壁濃度高于內(nèi)側(cè)壁濃度的結(jié)果[13]一致。同時(shí)，本研究發(fā)現(xiàn)，隨著曲率的增加，在曲率小于0.04前，內(nèi)外側(cè)NO濃度差呈上升趨勢，即兩側(cè)血管緊張程度的差異是逐漸上升的。當(dāng)曲率在0.04～0.1。范圍內(nèi)，兩側(cè)血管緊張程度的差異是隨著曲率的增加而緩慢降低的，這可能與血管自身調(diào)節(jié)有關(guān)，這種NO濃度分布的差異可能是導(dǎo)致血管形態(tài)改變以及血管病變的重要因素。

此外，本研究觀察到在彎曲處，內(nèi)外兩側(cè)W- B邊界(即EN區(qū)域)，NO濃度存在差異，外側(cè)濃度高于內(nèi)側(cè)，這種濃度差異隨曲率的增加呈先增加后緩慢下降的趨勢，其趨勢轉(zhuǎn)折出現(xiàn)在曲率為0.04時(shí)。Cozzi 等通過熒光探針實(shí)時(shí)監(jiān)測發(fā)現(xiàn)，NO的產(chǎn)生依賴于WSS，且隨著WSS的增加而增加[23]。Andrews等也通過實(shí)時(shí)測量，得出了WSS與NO產(chǎn)生的雙曲線關(guān)系[19]。然而，本研究中發(fā)現(xiàn)，在EN區(qū)域NO并沒有隨著WSS的升高而持續(xù)升高，因此筆者猜測在EN區(qū)域，NO的濃度分布除了受生成影響外，對流及消耗的影響也是不容忽視的，但其影響程度還需進(jìn)一步研究。

本研究通過多物理場耦合的仿真方法，對血管彎曲程度對內(nèi)源性NO濃度分布的影響有了更深入的了解。初步得到了隨彎曲程度變化的NO濃度分布規(guī)律。本研究主要考慮彎曲程度的影響，模型采用單一血管半徑、固定流速、血液黏度一定，即雷諾數(shù)Re為508.8固定值，因此在本研究討論的模型中Dean數(shù)對流動(dòng)的影響被簡化為曲率半徑對流動(dòng)的影響，然而對于其他血管尺寸以及非定常流因素的影響還需要進(jìn)一步研究。

4 結(jié)論

通過多物理場耦合仿真方法，發(fā)現(xiàn)血管彎曲程度對NO濃度分布有著重要的影響。彎曲血管的外側(cè)NO濃度明顯高于內(nèi)側(cè)，并且在一定范圍內(nèi)，濃度差值是逐漸上升的，當(dāng)彎曲程度超過臨界值后，濃度差值開始緩慢下降。該研究揭示了彎曲程度對內(nèi)源性NO分布的影響規(guī)律，這種NO濃度分布的差異可能是導(dǎo)致血管形態(tài)變化以及血管病變的重要因素，為進(jìn)一步探究慢性心血管疾病的發(fā)病機(jī)理，以及對不同彎曲程度的血管患病風(fēng)險(xiǎn)評估提供一定的理論依據(jù)。

[1] F?rstermann U, Sessa WC. Nitric oxide synthases: regulation and function [J]. European Heart Journal, 2012, 33(7): 829- 837.

[2] Iadecola C. Bright and dark sides of nitric oxide in ischemic brain injury [J]. Trends in Neurosciences, 1997, 20(3): 132- 139.

[3] 薛永亮,唐寧,華曉東,等. 一氧化氮與動(dòng)脈粥樣硬化 [J]. 中國動(dòng)脈硬化雜志, 2009, 17(8): 698- 701.

[4] Napoli C, Nigris F, Williams- Ignarro S, et al. Nitric oxide and atherosclerosis: an update[J]. Nitric Oxide, 2006, 15(4): 265- 279.

[5] Plank MJ, Wall DJN, David T. The role of endothelial calcium and nitric oxide in the localisation of atherosclerosis [J]. Mathematical Biosciences, 2007, 207(1): 26- 39.

[6] Chen Xuewen, Jaron D, Barbee KA, et al. The influence of radial RBC distribution, blood velocity profiles, and glycocalyx on coupled NO/O2 transport [J]. Journal of Applied Physiology, 2006, 100(2): 482- 492.

[7] Butler AR, Megson IL, Wright PG. Diffusion of nitric oxide and scavenging by blood in the vasculature [J]. Biochim Biophys Acta, 1998, 1425(1):168- 176.

[8] Lancaster JR. A tutorial on the diffusibility and reactivity of free nitric oxide [J]. Nitric Oxide, 1997, 1(1): 18- 30.

[9] Deonikar P, Kavdia M. A computational model for nitric oxide, nitrite and nitrate biotransport in the microcirculation: effect of reduced nitric oxide consumption by red blood cells and blood velocity [J]. Microvascular Research, 2010, 80(3): 464- 476.

[10] Chen Xuewen, Buerk DG, Barbee KA, et al. 3D network model of NO transport in tissue [J]. Medical & Biological Engineering & Computing, 2011, 49(6): 633- 647.

[11] Liu Xiao, Fan Yubo, Deng Xiaoyan, et al. Nitric oxide transport in an axisymmetric stenosis [J]. Journal of the Royal Society Interface, 2012, 9(75):2468- 78.

[12] 王中群, 袁偉, 楊永宗,等. 動(dòng)脈粥樣硬化病變位點(diǎn)特異性研究新進(jìn)展[J]. 中華老年心腦血管病雜志, 2013, 15(04):436- 437.

[13] Liu Xiao, Wang Zhenze, Zhao Ping, et al. Nitric oxide transport in normal human thoracic aorta: effects of hemodynamics and nitric oxide scavengers [J]. PLoS ONE, 2014, 9(11):e112395.

[14] Bouissou H, Emery MC, Sorbara R. Age related changes of the middle cerebral artery and a comparison with the radial and coronary artery [J]. Angiology, 1975, 26(3): 257- 268.

[15] Han Jintao, Qiao Huiting, Li Xuan, et al. The three- dimensional shape analysis of the M1 segment of the middle cerebral artery using MRA at 3T [J]. Neuroradiology, 2014, 56(11): 995- 1005.

[16] Rudic RD, Shesely EG, Maeda N, et al. Direct evidence for the importance of endothelium- derived nitric oxide in vascular remodeling [J]. Journal of Clinical Investigation, 1998, 101(4): 731.

[17] Shojima M, Oshima M, Takagi K, et al. Magnitude and role of wall shear stress on cerebral aneurysm computational fluid dynamic study of 20 middle cerebral artery aneurysms [J]. Stroke, 2004, 35(11): 2500- 2505.

[18] Leng X, Scalzo F, Ip HL, et al. Computational fluid dynamics modeling of symptomatic intracranial atherosclerosis may predict risk of stroke recurrence [J]. PLoS ONE, 2014, 9(5): e97531.

[19] Andrews AM, Jaron D, Buerk DG, et al. Direct, real- time measurement of shear stress- induced nitric oxide produced from endothelial cells in vitro [J]. Nitric Oxide, 2010, 23(4): 335- 342.

[20] Moore JE, Weydahl ES, Santamarina A. Frequency dependence of dynamic curvature effects on flow through coronary arteries [J]. Journal of Biomechanical Engineering, 2001, 123(2): 129- 133.

[21] 喬愛科, 劉有軍, 伍時(shí)桂. 彎曲動(dòng)脈的血流動(dòng)力學(xué)數(shù)值分析[J]. 計(jì)算力學(xué)學(xué)報(bào), 2003, 20(2): 155- 163.

[22] Cheng C, Tempel D, Oostlander A, et al. Rapamycin modulates the eNOS vs. shear stress relationship [J]. Cardiovascular Research, 2008, 78(1): 123- 129.

[23] Cozzi MR, Guglielmini G, Battiston M, et al. Visualization of nitric oxide production by individual platelets during adhesion in flowing blood [J]. Blood, 2015, 125(4): 697- 705.

The Effect of Different Curvature on the Distribution of Nitric Oxide Concentration in the Curved Arterial Segment

Zhao Hongjun1Han Jintao2Liu Cong1Dov Jaron3Qiao Huiting1#*

1(School of Biomedical Engineering, Beihang University, Beijing 100191, China)2(Department of Interventional Radiology and Vascular Surgery, Peking University Third Hospital, Beijing 100191，China)3(School of Biomedical Engineering, Science and Health Systems, Drexel University, Philadelphia, PA 19104, USA)

Nitric oxide (NO) produced by the endothelium plays an important role in regulating vascular tone. However, the effect of different curvature on the distribution of NO concentration in the curved arterial segment has not been clear. We developed eleven mass transport models with different curvature to simulate the distributions for NO in blood, vascular walls and surrounding tissue. Results showed that the radical distribution of NO was uneven in the curved segment. The outer concentration was higher than the inner one. And the concentration difference increased and then decreased with the increased curvature. In the smooth muscle region, the concentration difference reached the maximum when the curved height was increased to 2mm under the condition of the vessel projection length was a constant (20 mm),where the curvature was increased to 0.04,the outer boundary concentration was 23.74% higher than the inner one. This study revealed the effect of different curvature on NO concentration, and suggested that the uneven distribution of NO might cause vascular morphology changes and lead to vascular disease. Those models may provide a theoretical support for the further research and vascular disease risk assessment.

nitric oxide; curvature; multiphysics simulation; smooth muscle relaxation

10.3969/j.issn.0258- 8021. 2017. 03.008

2016-06-07，， 錄用日期:2016-08-27

國家重點(diǎn)研發(fā)計(jì)劃(2016YFC1304300)

R318

A

0258- 8021(2017) 03- 0316- 06

# 中國生物醫(yī)學(xué)工程學(xué)會(huì)高級會(huì)員(Senior member, Chinese Society of Biomedical Engineering)

*通信作者(Corresponding author)，E- mail: qht@buaa.edu.cn