內皮調節對小動脈管腔運動影響的模型分析1)

唐元梁 賀 纓

(大連理工大學能源與動力學院，大連116024)

內皮調節對小動脈管腔運動影響的模型分析1)

唐元梁 賀 纓2)

(大連理工大學能源與動力學院，大連116024)

微循環是血液和組織之間發生物質交換的主要場所，它可以通過改變管徑實現對血壓、血流量的局部調節.血管內皮層對小動脈運動有重要的調節作用.本文基于連續介質假設，建立了內皮調節過程中主要活性物質在管壁中的擴散--反應動力學模型，并分析了非線性黏彈性血管的徑向運動特性.利用該模型首先得到了內皮舒張因子一氧化氮(NO)、平滑肌細胞內鈣離子(Ca2+)以及磷酸化肌球蛋白(actin-myosin complexes,AMC)在管壁內的的徑向濃度分布；為分析內皮調節的動態過程，進一步對小動脈的被動舒張、血流量發生擾動時的管徑響應分別進行了模擬.研究結果顯示：當沒有活性物質參與調節小動脈被動舒張時，管徑無振蕩發生；當血流量變化引起內皮調節時，內皮舒張因子NO濃度和管徑均出現衰減振蕩，振蕩周期約60s.分析認為內皮調節對壁面剪切力的反饋控制，可能是NO濃度和管徑產生周期性振蕩的原因.內皮調節過程呈現的頻譜特征可以為血管內皮功能障礙的診斷提供幫助.

剪切力,內皮調節,一氧化氮,糖尿病,微循環

引言

人體組織或器官的正常代謝功能需要通過適量的營養和氧的供給，以及代謝廢物的及時排出來維持，生理狀況下人體通過全身性的神經/體液調節以及微循環的局部血流調節，能夠滿足組織或器官的代謝需求[1].血管內皮層結構及功能的完整性對于微循環的局部血流調節非常重要.動脈粥樣硬化、糖尿病等常伴有血管內皮功能障礙，造成人體局部血流調節機制失效.臨床研究表明神經、大血管及微血管病變是兩種類型糖尿病患者常見的慢性并發癥[2].考察血管活性物質濃度變化對血管運動的影響有助于加深對血管管腔內皮調節機理的認識，進一步為相關代謝、血管疾病的早期診斷及治療提供幫助.

1980年Furchgott等[3]發現內皮細胞對血管平滑肌細胞的舒張起重要作用.流動剪切力不僅影響和調控了內皮、平滑肌細胞的形狀和排列取向，而且能刺激內皮細胞分泌多種信號分子影響平滑肌細胞的生理功能[4].依賴流動剪切力的內皮調節反應主要發生在人體的小動脈及比較大的微動脈[5]，該反應可緩解甚至完全消除因肌源性反應造成的下游血管收縮、血流持續下降的問題[6]；此外，由于代謝反應主要引起末梢微動脈的有限舒張，不能完全滿足周圍組織在代謝量增大時的血流需求，因此上游小動脈的內皮舒張反應能極大彌補末梢微動脈代謝調節反應的局限性.

目前一氧化氮(NO)被認為是剪切力誘導內皮細胞釋放的最主要舒血管因子，NO的產生、傳輸及作用是影響血管內皮調節反應機制的關鍵.細胞內重要的第二信使Ca2+在內皮細胞NO合成過程中起到重要作用，剪切力或其他生化因子引起內皮細胞膜上Ca2+通道打開，內皮細胞合成NO增多.Li等[7-9]對剪切流與ATP(adenosine triphosphate，三磷酸腺苷)協同定量調控血管內皮細胞內Ca2+響應進行了理論和實驗研究，通過控制剪切力或ATP信號較好地實現了對內皮細胞內Ca2+信號的定量控制.管壁內皮NO的產生研究方面，Liao等[10]通過實驗獲得了在無流動及有血流情況下血管不同的肌源性收縮反應，并基于實驗數據建立了相關的內皮調節反應經驗公式.Cornelissen等[11]建立了內皮分泌NO速率與壁面剪切力的S形函數關系，并進一步發展了NO濃度與冠狀動脈樹血管緊張度關系的數學模型.關于NO在管壁、管腔及周圍組織的傳輸特性的研究工作方面，Chen等[12]構建了包含小動脈、毛細血管網及組織的三維NO傳輸模型，對NO在血液、血管壁及周圍組織中的分布進行了模擬分析.結果表明壁面剪切力、紅細胞壓積、血管分叉等因素對NO的生成及分布均有影響.Liu等[13]基于流動模型和對流--擴散--反應方程對NO在狹窄血管中的傳輸進行了模擬分析.研究結果表明血管狹窄處的流動失穩對NO的傳輸有顯著的影響，可能是造成動脈粥樣硬化的重要因素.Tsoukias等[14]則對組織總肌紅蛋白、血漿中血紅蛋白對微循環周圍NO的擴散影響進行了模型分析.

目前階段研究者多采用唯象模型來研究小動脈血流調節反應，即不考慮具體的調節反應過程，只是概括地去總結和提煉實驗結果.如Carlson等[15-17[18]建立了血管平滑肌運動的電化學--化學機械細胞模型，在細胞力學模型中加入了與主動收縮張力元件串聯的黏彈性元件，隨后Yang等[19]加入了平滑肌細胞內NO/cGMP(cyclic guanosine monophosphate,環磷酸鳥苷)反應模型，考察了NO的舒血管作用.Kudryashov等[20]基于連續介質模型，預測了NO、平滑肌細胞內Ca2+以及磷酸化肌球蛋白在管壁的穩態濃度分布，并通過改變血流量考查了模型的自動調節功能以及血管口徑在向新平衡態發展過程中出現的振蕩現象.

在人體的許多組織中，小動脈管徑及血流的自發振蕩現象被觀測到，但是其生理意義及振蕩機理仍不甚清楚.研究表明血流/皮膚溫度振蕩信號與相應調節反應間的頻率具有對應關系，可以從皮膚溫度振蕩信號提取出內皮調節信號(0.0095～0.02Hz)，(交感)神經調節信號(0.02～0.06Hz)及肌源性調節信號(0.06～0.15Hz)等[21-22].Smirnova等[23]對比了健康組、二型糖尿病組以及糖耐量受損患者的手部皮膚溫度信號，發現無血管病變實驗對象在手部受到冷水刺激時內皮、神經調節頻段的溫度波動幅值均減小，冷水刺激結束后波動幅值逐漸恢復；而二型糖尿病和糖耐量受損患者在結束冷水刺激后，相應頻段的溫度波動幅值并沒有及時恢復.其認為這可能是由于微血管內皮功能障礙及神經病變引起的.由于多種調節反應可能共同參與血管的調節過程且互相影響，單純的實驗研究難以準確把握某一調節反應的內在機理，因此結合合適的數學模型對于管徑調節機理的研究會有很大的幫助.

由于唯象模型不能夠對內皮調節反應機制作出解釋，本文將基于Kudryashov等[20]的工作對小動脈內皮調節功能進行模型分析.基于連續介質假設,對內皮調節過程主要介質在管壁中的兩層擴散--反應動力學變化進行了建模，并考慮管壁的非線性黏彈性性質計算血管的徑向運動.利用該模型本文將對給定血流條件下內皮舒張因子NO以及平滑肌細胞內Ca2+和AMC的徑向濃度分布進行計算，進一步模擬血流變化引起的內皮調節反應過程，并對管徑振蕩的力學機制等進行討論.

1 研究方法

1.1 小動脈壁力學性質及依賴剪切力的內皮調節反應生理機制

小動脈管壁由幾種不同的組織構成.小動脈壁內膜表面是單層內皮細胞，中膜含有大量的平滑肌細胞.小動脈管壁的力學性質表現為非線彈性和黏彈性，同時平滑肌在生理控制下會提供收縮張力.

平滑肌細胞伸縮的初始刺激源是細胞內Ca2+濃度的變化.Ca2+引起平滑肌細胞收縮的過程可以描述如下[24]：平滑肌細胞內Ca2+與鈣調蛋白結合；鈣調蛋白與Ca2+結合物進一步與肌球蛋白激酶(一種磷酸化酶)結合且將其激活；在肌球蛋白激酶作用下，肌球蛋白頭調節輕鏈磷酸化；輕鏈發生磷酸化后，肌球蛋白頭對肌動蛋白絲進行間歇性的牽拉，從而引起平滑肌緊張并產生收縮力.

流動剪切力調節血管內皮細胞的代謝，使之合成血管活性物質NO等調節平滑肌細胞內Ca2+濃度及血管的緊張度.依賴剪切力的血管內皮舒張過程如圖1所示：當血流量增加、管壁剪切力增大時，內皮細胞膜上Ca2+通道打開，內皮細胞合成NO增加；NO合成后擴散進入管壁平滑肌層；NO穿越平滑肌細胞膜進入細胞內部，刺激平滑肌細胞產生cGMP，后者能夠引起細胞內Ca2+的外流；如前所述，Ca2+濃度對平滑肌細胞的收縮性和張力起到決定作用，其流失導致細胞內收縮性的肌動--肌球蛋白復合物濃度下降，最終導致平滑肌的松弛和張力的下降.

圖1 血管內皮調節圖示[24]Fig.1 Schematic of endothelial regulation[24]

1.2 數學模型

1.2.1 內皮層NO的產生

內壁面切應力的大小決定了內皮細胞產生NO的速率，軸向血流速度分布采用Hagen--Poiseuille模型，給定血流量Q，則內壁面剪切應力τsh(t)可以表示為

式中，r為徑向坐標，R(t)為血管內徑，μ為血液黏度，vx(r)為軸向流速，Q為血液體積流量.內壁面切應力的大小決定了內皮細胞產生NO的速率，此外血液--特別是血液中的血紅蛋白，會消耗一部分內皮細胞產生的NO，因此內皮層NO濃度的變化可用常微分方程(2)描述

式中ne內皮層NO的濃度，等式右邊第1項為血流對NO的消耗率，ke為NO在血液中的消耗系數；右邊第2項表示剪切力依賴的內皮NO產生率，其與壁面剪切力τsh(t)成正比，k1為產生系數.也有研究者[11]認為NO的產生率與壁面剪切力成S形曲線關系.血液中消耗NO的主要成分是血紅蛋白，其消耗速率與血液的流變學參數相關，在本文中ke與k3均設為常數.將壁面剪切力式(1)代入方程(2)中，并引入管徑相對擾動量η=(R-R0)/R0替代管徑R，其中R0為參考半徑，得到

1.2.2 管壁NO的傳輸

NO在內皮細胞合成后，經過擴散和生物化學反應在管壁中傳輸.內膜和中膜NO濃度變化的動力學過程可以表示為如下方程

其中，j=1,2時分別表示管壁的內膜、中膜，nj為管壁NO濃度.方程右邊第1項為擴散項，Dj為NO在管壁中的擴散系數；第2項為NO在管壁的反應項.與血液對NO的消耗類似，模型認為NO在管壁的反應速率與NO濃度成正比，δj為NO在管壁的反應系數.利用式(3)可以給出內皮層的NO濃度并將其作為管壁內邊界(r=Ri)的邊界條件；內膜與中膜交界面(r=Rm)上NO濃度及擴散量均保持連續；在中膜外邊界(r=Ra)上，假設NO不外擴散.上述邊界條件描述為

式中D1和D2分別為血管內膜、中膜的NO擴散系數.方程(3)、方程(4)、邊界條件(5)及合適的初值條件構成了完整的動脈管壁NO擴散--動力學反應雙層模型.

1.2.3 平滑肌Ca2+與磷酸化肌球蛋白濃度

平滑肌細胞內Ca2+的濃度變化由平滑肌細胞質Ca2+的凈流量決定.平滑肌細胞質中Ca2+源包括細胞間質、細胞內的細胞器--肌質網.與骨骼肌細胞不同，血管平滑肌細胞的肌質網并不發達，平滑肌細胞質Ca2+的主要來源是細胞間質.由于細胞間質Ca2+濃度遠高于平滑肌細胞內(約104倍)，Ca2+可通過平滑肌細胞膜上的Ca2+通道進入細胞內.該動力學過程可以用方程(6)描述

式中,C為平滑肌細胞質內Ca2+濃度，方程右邊第一項為Ca2+向細胞外的主動運輸，主動運輸系數α(n)受平滑肌層NO濃度影響，α(n)=α1+k2n2(r,t)；第二項為Ca2+的被動運輸過程.β為傳輸系數，當跨壁壓差不變時，β取為常數.Ce為細胞間質Ca2+濃度.因Ce?C，可將方程(6)簡化為

其中φ0=βCe.

平滑肌細胞內磷酸化肌球蛋白的濃度直接決定了平滑肌細胞的收縮性及主動張力的大小.其磷酸化程度(即濃度大小)由細胞內Ca2+濃度、閾值濃度及脫磷酸化速率共同決定.該動力學過程可以描述為

式中f為磷酸化肌球蛋白濃度，方程右邊第1項表示肌球蛋白脫磷酸化過程，α2為脫磷酸化系數；第2項表示Ca2+依賴的肌球蛋白磷酸化，Cth為Ca2+閾值濃度，γ為磷酸化系數，θ為Heaviside單位函數

1.2.4 管壁徑向運動

圖2所示為管壁受力狀態，從圖2可以看出，小動脈的管徑由跨壁壓差P0=P-Pe(P和Pe分別為管內、外壓力)、管壁周向張力T決定.

圖2 管壁受力狀態Fig.2 Stress condition on vessel wall

對于質量為Δm、半徑為R的血管單元，根據動量守恒，有

式中Δm=ρw2πRΔxh，ρw為血管密度，h為壁厚，Δx為血管單元的長度，fr和fp分別為管壁周向張力T、跨壁壓差P0產生的徑向力.fr和fp可分別表示為

外力平衡時有fr+fp=0，代入方程(11)可得T=PR0，即周向張力與跨壁壓差滿足拉普拉斯定律.當血管為薄壁管時，管壁周向應力可認為均勻分布，即T=σθθh，其中σθθ為管壁平均周向應力.考慮動脈壁的非線彈性、黏性及平滑肌產生的主動收縮張力等情況，σθθ可以表示為

式中ξ為泊松比，κ1為非線性彈性系數，λ為黏性系數.方程右邊最后一項為平滑肌主動收縮張力，其與平滑肌磷酸化肌球蛋白的體積平均濃度F成正比，k3為比例系數.管壁楊氏模量E(F)的大小也與F有關(E(F)=E0(1+εF/F0),F0為跨壁壓差P0=15.332kPa),Ψ0=0.1μmol/s時磷酸化肌球蛋白的體積平均濃度，E0為不考慮平滑肌主動張力時的楊氏模量.

將式(11)和式(12)代入式(10)，并引入血管半徑的相對擾動量η=(R-R0)/R0,可以得到關于半徑相對擾動量η的方程

1.3 參數選取與數值方法

式 (13)中的參數 λ = 1.5E0t0/R0單位為kg·m-2s-1,其中t0=10s；k3=P0R0/(F0h0)[20].其他參數的取值及來源如表1所示.

表1 模型參數Table 1 Model parameters

引入如下無量綱參數對模型進行無量綱化處理

其中特征參數n0=0.1μmol/L，Cth=0.1μmol/L.R0=2mm為跨壁壓差P0=15.332kPa,Ψ0=0.1μmol/s時的血管內徑.最后得到相應的無量綱模型為

采用有限差分法對上述模型進行數值求解.對方程(16)擴散項進行空間中心差分離散，并進行隱式求解，采用共軛梯度法進行迭代；方程 (15)、方程(17)、方程(18)均采用歐拉隱式格式求解；對方程(19)慣性項、黏性項進行中心差分離散.無量綱時間、空間步長分別為Δt′=3.0×10-5和Δr′=1.2×10-3.

2 結果分析

2.1 內皮調節主要介質濃度的穩態分布

基于上述模型，首先得到了Ψ0=0.1μm/s時內皮調節主要介質濃度在管壁的穩態分布.如圖3所示，由于內皮細胞是管壁NO的唯一產生源，在內皮層的NO濃度最高；沿管徑向外，隨著r′的增大,NO濃度因在管壁的反應消耗而逐漸降低(中膜外層的NO濃度為0.9030).平滑肌細胞內Ca2+濃度沿徑向方向略有上升(中膜內層濃度為1.1296，中膜外層濃度為1.1314)，這反映了舒血管物質NO的濃度變化與Ca2+變化相反，NO濃度的增大會造成平滑肌細胞內Ca2+濃度的下降.由于Ca2+濃度高于敏感閾值(Cth=0.1μmol/L)，Ca2+濃度依賴的肌球蛋白磷酸化作用處于激活態，因此沿徑向磷酸化肌球蛋白濃度也隨Ca2+濃度的上升而略有增大(中膜內層濃度為0.9969，中膜外層濃度為1.0104).

圖3 內皮調節主要介質在管壁內的徑向濃度分布Fig.3 Stationary concentration distributions in the arterial wall of the key agents

2.2 血管的被動舒張

當管壁NO或其他舒血管物質濃度上升時，平滑肌細胞內Ca2+濃度會逐漸下降并可能低于閾值濃度Cth.此時平滑肌完全放松，主動收縮張力消失，動脈壁模型退化為完全被動的黏彈性管，此時管徑相對擾動變化方程(13)變為

在η(0)=0的初始狀態下，管壁周向應力減小使得跨壁壓差與周向應力不平衡，管徑會發生被動擴張.文獻[20]給出了血管被動擴張模型的解析解，其形式為

圖4為該解析解與我們得到的數值解的比較.結果表明管徑從初始狀態η(0)=0經過被動擴張過程后達到新平衡態η=1.2511，數值解與精確解吻合較好.

圖4 管壁的被動擴張變化Fig.4 Passive dilatation of the artery

2.3 血流量變化時的內皮調節過程

由方程(2)可知，當管內血流量變化時，壁面剪切力依賴的內皮NO合成會發生改變，并影響到血管平滑肌的收縮張力、管徑大小.將Ψ0=0.1μmol/s時血管管徑、內皮調節主要介質的穩態濃度作為初始條件，對血流量發生變化時管徑的內皮調節過程進行了模擬分析.

圖5為血流量增加25%時模型的響應情況.管壁NO濃度均從初始態開始經過衰減振蕩，逐漸穩定于新的平衡態濃度.血管管徑也呈現衰減振蕩，最終在新平衡態(η1.25=0.0536)穩定下來，此時跨壁壓差與管壁周向應力達到平衡.結果顯示NO濃度、管徑的振蕩周期均為60s左右，這與依賴內皮的血流調節頻率相一致(0.0095～0.02Hz[21]).

圖5 NO濃度和管徑隨時間的變化Fig.5 Time dependence of the changes of NO concentration and artery radius

表 2列出了不同血流條件下小動脈達到新平衡態時的管徑擾動量和壁面剪切力.從該表可以看到，管徑擾動和壁面剪切應力隨血流量增加呈非線性增加.當血流量Q由0.75Q0增大到1.50Q0時，壁面剪切力τsh只增大了0.20τ0.由式(1)可知，當管徑無變化(即η=0)時，壁面剪切力τsh會隨血流量Q的增大而線性增大；內皮調節引起的管徑增大在很大程度上限制了τsh的快速上升，這對于維持血管內皮受力的相對穩定性非常重要.反之，若壁面剪切力τsh維持不變，由式(1)可得則管徑擾動隨血流量變化曲線斜率隨血流量Q增大而逐漸減小.如圖6實線所示，利用內皮調節模型得到的數據也呈現了類似的特征，說明內皮調節的作用是使血管管壁切應力盡量保持穩定.

表2 不同血流量條件下管徑和壁面剪切力的響應Table 2 Responses of vessel radius and shear stress under dif f erent fl w conditions

圖6 管徑擾動隨血流量的變化Fig.6 Changes of artery radius under dif f erent fl w conditions

3 討論

3.1 管徑振蕩變化的產生機理

通過模型分析可看出，血流量變化時，內皮舒張因子NO和管徑會產生大周期振蕩.為探究該振蕩發生的原因，本文進一步對活性物質動力學特性和血管運動的相互作用進行了分析.

在血管管徑運動方程(13)中我們考慮了血管的慣性、黏性及彈性.這些力學性質的組合可能會產生管壁阻尼振蕩運動.阻尼運動方程的一般形式為

另一方面，若假設血管管徑無變化(η=0)，則根據方程(15)～方程(18)得到的平滑肌細胞內NO，Ca2+及磷酸化肌球蛋白濃度變化如圖7所示.它們均從初始態平緩過渡至新的平衡態，無振蕩出現.

圖7 內皮調節介質濃度隨時間的變化(η=0,Q=1.25Q0)Fig.7 Variations of key agents’concentrations(η=0,Q=1.25Q0)

上述分析表明，血管管徑、內皮調節活性物質濃度均不會獨立地產生振蕩變化，這種振蕩可能是由于對壁面剪切力的負反饋控制作用而產生.圖8描述了這一負反饋控制過程示意圖，本文所采用模型很好地體現了這一特點.

圖8 內皮調節過程的剪切力反饋控制Fig.8 Feedback control of shear stress in the process of endothelial regulation

血管的自發振蕩現象需要管壁平滑肌緊張度的協同變化，而能夠在血管尺度迅速引發這一響應的方式極有可能是電信號在細胞之間的傳輸[28].有研究認為膜電勢及細胞內Ca2+濃度的振蕩變化是引起細胞緊張度變化的原因，這些電信號位于血管分叉處少量的特殊“起搏器”細胞內[29].Secomb[30]總結前人工作指出血管自發振動機理的兩種可能性：血管壁細胞本身固有振動特性，即“起搏器”性質；管壁振蕩運動由細胞反應、管壁力學性質以及血液流動之間的相互作用引起.Kvernmo等[31]認為內皮活動頻率或反映了內皮釋放血管活性物質的速率，或反映了平滑肌細胞對這些物質的響應頻率.本文模型基于宏觀連續介質假設，不能從細胞／細胞群尺度分析血管振蕩的內在機理，但該模型體現了對壁面剪切力的反饋控制特點，使得內皮舒張因子NO、平滑肌內Ca2+、磷酸化肌球蛋白濃度以及血管管徑逐步得到校正.因此我們推測內皮調節過程的反饋控制特點可能是內皮舒張因子NO、管徑周期性振蕩的原因；反饋環節的缺失使得振蕩消失.

3.2 模型的特點和局限性

本文發展的內皮調節模型將活性物質濃度變化與血管周向張力、血管管徑變化關聯起來.相比于一般的唯象模型，本模型能夠體現內皮調節反應過程，而且不涉及過于復雜的單個細胞內的具體反應機理，有助于我們加深對血管管腔內皮調節機理的認識.同時本模型有如下幾個方面需要說明和進一步改進.

本文所采用的活性物質反應動力學模型與Kudryashov等[20]的模型相同，但在分析小動脈徑向運動時，考慮到方程(13)中慣性力項的作用很小，又發現內皮調節引起的管徑大周期波動并不是來自于管壁的機械振蕩，因此在模型分析中略去了該慣性項.相較于Kudryashov等[20]的工作，本模型對血管運動振蕩特性的分析更深入.

在流體介質方面，本模型假設血液為牛頓流體.采用牛頓流體假設會對壁面切應力的大小有一定影響，但不會影響活性物質在管壁的分布規律和管壁運動特性.

目前，本模型假定內皮NO的釋放率與壁面剪切力成正比(k1τsh)，且k1假定為定值.當內皮發生功能障礙時，NO的釋放減少，模型中對應的比例系數k1也會相應減小，最終會影響血管正常的內皮舒張功能和壁面剪切力的維持.類似地，本模型假定血液對NO的消耗率與內皮NO濃度呈正比其中ke為定值.NO在血液中的傳輸和反應實際上與血液流速、紅細胞壓積、紅細胞膜的滲透性等因素有關，這些因素的變化會改變ke的取值.在另外的研究中[32]，建立了內皮NO在血液和組織中的反應—對流擴散模型，并對ke的影響進行了考察，發現ke的增大會使壁面NO濃度降低；同時，也在考察紅細胞運動，變形和壓積對NO生成的影響，這些結果會為確定不同狀態下血液對NO的消耗率提供很好的數據支撐.

4 結論

本文結合擴散、反應動力學模型和黏彈管模型，對剪切力依賴的小動脈內皮調節反應進行了模擬，并對調節引起的振蕩現象進行了分析.基于連續介質假設模型對管壁內皮層NO的產生和NO在管壁中的兩層擴散--反應動力學過程、平滑肌細胞內Ca2+及磷酸化肌球蛋白的反應動力學過程進行了數學描述；考慮管壁的黏彈性并結合內皮調節產生的主動張力變化，建立了血管的徑向運動模型.研究結果表明小動脈被動舒張過程管徑無振蕩發生，而血流量變化引起的內皮調節使內皮舒張因子NO濃度和管徑均出現衰減振蕩.本文認為內皮調節對壁面剪切力的反饋控制，可能是NO濃度—管徑產生周期性振蕩的原因.

內皮調節過程呈現的頻譜特征可以為內皮功能障礙的診斷提供幫助.內皮功能受損常見于糖尿病發病早期.糖尿病內皮功能障礙的確切機制還不明確，高血糖造成的多種因素可能參與了內皮的受損過程，包括NO的產生、傳輸及效應環節，即：NO合成下降，NO擴散至平滑肌細胞層的功能受損，平滑肌細胞對NO反應性下降或內皮衍生的收縮因子抑制了血管舒張.合適的數學模型可以實現對上述環節的合理描述，通過模型分析可有助于對內皮功能受損原因的排查.

1 Davis MJ,Hill MA,Kuo L.Local Regulation of Microvascular Perfusion.Comprehensive Physiology,2010

2 Bohlen,Glenn H.Microvascular Consequences of Obesity and Diabetes.Comprehensive Physiology,2010

3 Furchgott RF,Zawadzki JV.The obligatory role of endothelial cells in the relaxation of arterial smooth muscle by acetylcholine.Nature,1980,288(5789)：373-376

4 任長輝，劉肖，康紅艷等.剪切力條件下血管內皮細胞與平滑肌細胞的相互作用.醫用生物力學,2015,30(2)：185-191(Ren Changhui,Liu Xiao,Kang Hongyan,et al.Interactions between vascular endothelial cells and smooth muscle cells under shear stress.Journal of Medical Biomechanics,2015,30(2)：185-191(in Chinese))

5 Pohl U,De WC,Gloe T.Large arterioles in the control of blood fl w：role of endothelium-dependent dilation.Acta Physiologica Scandinavica,2000,168(168)：505-510

6 Koller A,Kaley G.Endothelial regulation of wall shear stress and blood fl w in skeletal muscle microcirculation.American Journal of Physiology,1991,260(260)：862-868

7 LiLF,ChengX,QinKR.ModelingofTRPV4-C1-mediatedcalcium signaling in vascular endothelial cells induced by flui shear stress and ATP.Biomechanics&Modeling in Mechanobiology,2015, 14(5)：1-15

8 陳宗正,覃開蓉,高爭鳴等.用于分析剪切力和生化因子對細胞聯合作用的微流控裝置.水動力學研究與進展,2015,30(6)：643-649 (Chen Zongzheng,Qin Kairong,Gao Zhengming,et al.A microflu idic device for analyzing cellular behavior in response to shear stress and biochemical factors.Chinese Journal of Hydrodynamics,2015, 30(6)：643-649(in Chinese))

9 高爭鳴,覃開蓉,陳宗正等.剪切力和ATP定量調控血管內皮細胞鈣信號的數值仿真研究.水動力學研究與進展,2015,30(6)：679-684(Gao Zhengming,Qin Kairong,Chen Zongzheng,et al. Numerical simulation research on quantitative regulation of intracellular Ca2+signals in vascular endothelial cells via shear stress and ATP modulation.Chinese Journal of Hydrodynamics,2015,30(6)：679-684(in Chinese))

10 Liao JC,Kuo L.Interaction between adenosine and fl w-induced dilation in coronary microvascular network.AJP Heart&Circulatory Physiology,1997,272(4Pt2)：H1571-H1581

11 Cornelissen AJ,Dankelman J,Vanbavel E,et al.Balance between myogenic,fl w-dependent,and metabolic fl w control in coronary arterial tree：a model study.AJP Heart&Circulatory Physiology, 2002,282(6)：H2224

12 Chen X,Buerk DG,Barbee KA,et al.3D network model of NO transport in tissue.Medical&Biological Engineering&Computing,2011,49(6)：633-647

13 Liu X,Fan Y,Xu XY,et al.Nitric oxide transport in an axisymmetric stenosis.Journal of the Royal Society Interface,2012,9(75)：2468-2478

14 TsoukiasNM,PopelAS.Amodelofnitricoxidecapillaryexchange.Microcirculation,2003,10(6)：479-495

15 Carlson BE,Secomb TW.A theoretical model for the myogenic response based on the length-tension characteristics of vascular smooth muscle.Microcirculation,2005,12(4)：327-338

16 Arciero JC,Carlson BE,Secomb TW.Theoretical model of metabolic blood fl w regulation：roles of ATP release by red blood cells and conducted responses.AJP Heart&Circulatory Physiology,2008,295(4)：H1562-H1571

17 Carlson BE,Arciero JC,Secomb TW.Theoretical model of blood fl w autoregulation：roles of myogenic,shear-dependent,and metabolic responses.AJP Heart&Circulatory Physiology,2008, 295(4)：H1572-H1579

18 Yang J,Jr J C.The myogenic response in isolated rat cerebrovascular arteries：smooth muscle cell model.Medical Engineering&Physics,2013,25(8)：691-709

19 Yang J,Clark JW,Bryan RM,et al.Mathematical modeling of the nitric oxide/cGMP pathway in the vascular smooth muscle cell.AJP Heart&Circulatory Physiology,2005,289(2)：H886-H897

20 Kudryashov NA,Chernyavskii IL.Numerical simulation of the process of autoregulation of the arterial blood fl w.Fluid Dynamics, 2008,43(1)：32-48

21 Kvernmo HD,Stefanovska A,Bracic M,et al.Spectral analysis of the laser doppler perfusion signal in human skin before and after exercise.Microvascular Research,1998,56(3)：173-182

22 Sagaidachnyi AA,Usanov DA,Skripal AV.Correlation of skin temperature and blood fl w oscillations//Proceedings of SPIEThe International Society for Optical Engineering,2011,8337(4)：83370A-83370A-8

23 Elena S,Sergey P,Irina M,et al.Assessment of endothelial dysfunction in patients with impaired glucose tolerance during a cold pressor test.Diabetes&Vascular Disease Research Of ficial Journal of the International Society of Diabetes&Vascular Disease,2013, 10(6)：489-497

24 Hall JE.Textbook of Medical Physiology.13th Edn.Oxford：Elsevier Ltd,2016.203-214

25 Formaggia L,Lamponi D,Quarteroni A.One-dimensional models for blood fl w in arteries.Journal of Engineering Mathematics, 2003,47(3-4)：251-276

26 Regirer SA,Nkh S.Mathematical models of nitric oxide transport in a blood vessel.Biofizik,2005,50(3)：515-536

27 Snow HM,Mcaulif f e SJ,Moors JA,et al.The relationship between blood fl w and diameter in the iliac artery of the anaesthetized dog：the role of endothelium-derived relaxing factor and shear stress.Experimental Physiology,1993,16(3)：278-282

28 Nilsson H,Aalkjaer C.Vasomotion：mechanisms and physiological importance.Molecular Interventions,2003,3(2)：79-89,51

29 Goligorsky MS,Colfles D,Gordienko D,et al.Branching points of renal resistance arteries are enriched in L-type calcium channels and initiate vasoconstriction.American Journal of Physiology,1995, 268(2 Pt 2)：F251

30 Secomb TW.Theoretical models for regulation of blood fl w.Microcirculation,2008,15(8)：765-775

31 Kvernmo HD,Stefanovska A,Kirkeboen KA,et al. Oscillations in the human cutaneous blood perfusion signal modifieby endothelium-dependent and endothelium-independent vasodilators.Microvascular Research,1999,57：298-309

32 Wei YJ,He Y,Tang YL.Finite element analysis on NO transport in system of permeable capillary and tissue//11th Asian Computational Fluid Dynamics Conference.Accepted

MODEL ANALYSIS OF ENDOTHELIUM-DEPENDENT VASOMOTION OF SMALL ARTERY1)

Tang Yuanliang He Ying2)
(School of Energy and Power Engineering，Dalian University of Technology，Dalian116024，China)

Metabolic substance exchange between blood and tissue occurs mainly in microcirculation,which can locally regulate blood pressure and blood fl w by changing their diameters.Vascular endothelium plays an important role in the vasomotor regulation of small artery.In this paper,a model was adopted to study the endothelial regulation mechanism. Based on the continuum assumption,two layers of dif f usion&kinetic processes of key agents in endothelial regulation were modeled,and the nonlinear viscoelastic properties of the wall material were considered in the computation of radial motion of small artery.The stationary distributions of nitric oxide(NO),calcium ion(Ca2+)and the contracting actinmyosin complexes(AMC)concentrations in the wall were firstl obtained;then the process of arterial passive dilation and the autoregulation process to the disturbance of blood fl w were simulated.Numerical results showed that there was no oscillation of arterial diameter occurred during the passive dilation process.However,when there was a change in blood fl w,the whole system transferred from the initial state to a new equilibrium state with slow damped oscillations. The oscillating period was about 60s.It is supposed that the endothelial oscillation of artery diameter and NO concentration occurring during the dynamic regulating process is caused by the feedback control of shear stress.This oscillationcharacteristic can provide useful information for the diagnosis of endothelial dysfunction.

shear stress,endothelial regulation,nitric oxide,diabetes,microcirculation

R318.01

A doi：10.6052/0459-1879-16-171

2016-06-13收稿，2016-09-29錄用,2016-10-18網絡版發表.

1)國家自然科學基金(51576033)和大連市科技計劃(2015F11GH092)資助項目.

2)賀纓，教授,主要研究方向：血液動力學，生物傳熱傳質.E-mail：heying@dlut.edu.cn

唐元梁,賀纓.內皮調節對小動脈管腔運動影響的模型分析.力學學報,2017,49(1)：182-190

Tang Yuanliang,He Ying.Model analysis of endothelium-dependent vasomotion of small artery.Chinese Journal of Theoretical and Applied Mechanics,2017,49(1)：182-190