陳國泰 鄭艷紅 易 丹 曾巧云

(福建師范大學(xué)數(shù)學(xué)與統(tǒng)計學(xué)院,福州 350117)

(福建師范大學(xué)福建省分析數(shù)學(xué)及應(yīng)用重點實驗室,福州 350117)

(福建師范大學(xué)福建省應(yīng)用數(shù)學(xué)中心,福州 350117)

引言

帕金森病(Parkinson’s disease,PD)是一種神經(jīng)性大腦疾病,其癥狀包括靜止性震顫、運動遲緩、肌強直和姿勢步態(tài)障礙等[1-2].據(jù)統(tǒng)計,全球約有800 萬人患有帕金森病,其中65 歲以上人群的帕金森患病率約為1%～3%[3].我國龐大的人口基數(shù)之下,人口老齡化的問題凸顯,已知的PD 患者和潛在的患病數(shù)給社會帶來了一定的挑戰(zhàn).

有研究表明,PD 的主要病理特征是由大腦黑質(zhì)致密部中多巴胺能神經(jīng)元的大量退化導(dǎo)致的[4].黑質(zhì)致密部的變性進一步影響了整個基底神經(jīng)節(jié)—丘腦—皮質(zhì)(basal ganglia-thalamic-cortex,BGTC)環(huán)路的神經(jīng)元集群活動,特別是基底神經(jīng)節(jié)(BG)中的丘腦底核(subthalamic nucleus,STN)和蒼白球(globus pallidus,GP)神經(jīng)元集群表現(xiàn)出異常同步的β 頻段(13～35 Hz)振蕩[5].由于病理性的β 振蕩動力學(xué)的機制尚不明確,因此需要進行更多的研究來探索帕金森病異常β 振蕩產(chǎn)生的機制.

國內(nèi)外許多學(xué)者利用神經(jīng)元模型探索了BGTC環(huán)路中產(chǎn)生β 振蕩的起源問題[6-11].Holgado 等[12]認為帕金森病的振蕩行為起源于STN 與蒼白球外側(cè)(globus pallidus externa,GPe)兩個神經(jīng)元集群之間的相互作用.他們通過對系統(tǒng)穩(wěn)定性進行分岔分析發(fā)現(xiàn)當(dāng)突觸連接強度增加時,可誘導(dǎo)STN-GPe 回路中異常同步β 振蕩的產(chǎn)生,并發(fā)生Hopf 分岔,使系統(tǒng)在健康和帕金森病這兩個狀態(tài)之間轉(zhuǎn)換.Hu 等[13-14]建立了一個兩STN 神經(jīng)元集群的STN-GPe 網(wǎng)絡(luò)模型,進一步拓展了STN-GPe 神經(jīng)回路,使得模型更加完善.此外,還有學(xué)者認為,除了STN-GPe 環(huán)路,在BG 之外的皮質(zhì)(cortex,CTX)對病理性的β 振蕩也具有重要的調(diào)節(jié)作用[15-19].Pavlides 等[20]基于CTX 和STN-GPe 環(huán)路在PD 中產(chǎn)生β 振蕩的重要作用,提出了CTX-STN-GPe 模型,其中,皮質(zhì)分為興奮性和抑制性兩個神經(jīng)元集群[20],通過此模型成功地重現(xiàn)了Nambu 等[21]的實驗結(jié)果.另外,突觸連接的強度總在神經(jīng)元活動中起著重要的作用[22-24],例如阻斷興奮性突觸傳遞能夠減少神經(jīng)元的放電活動.有研究表明,帕金森病中多巴胺能神經(jīng)元的缺失可能會影響樹突的興奮性[25].劉晨[26]基于CTXSTN-GPe 回路建立了一個雙振蕩神經(jīng)群模型,探討了由多巴胺耗盡引起的突觸連接的變化對回路中各個神經(jīng)元集群的影響.

許多學(xué)者在研究丘腦底核和蒼白球環(huán)路時,并未明確區(qū)分蒼白球外側(cè)(GPe)和蒼白球內(nèi)側(cè)(globus pallidus-interna,GPi),而是將兩者視為整體.Plenz等[27]通過觀察基底神經(jīng)節(jié)的振蕩特性,發(fā)現(xiàn)STN,GPe 以及GPi 之間復(fù)雜的交互關(guān)系能夠誘導(dǎo)振蕩的產(chǎn)生.GPi 作為基底神經(jīng)節(jié)的輸出核團,在皮層和丘腦之間的信號傳遞過程中發(fā)揮著關(guān)鍵的作用,Pavlides 等[20]提出的CTX-STN-GPe 生理學(xué)共振模型尚未考慮蒼白球內(nèi)部之間的關(guān)聯(lián),而僅僅關(guān)注了蒼白球與其外部神經(jīng)元集群的聯(lián)系.因此,本文將GPi 考慮進去,建立更加完善的皮質(zhì)-丘腦底核-蒼白球外側(cè)-蒼白球內(nèi)側(cè)模型(Cortex-STN-GPe-GPi),這更具有生理意義.并在第二節(jié)具體研究了與GPi 有關(guān)的參數(shù)對其產(chǎn)生振蕩的影響.該模型有助于對基底神經(jīng)節(jié)重要輸出核團GPi 的研究,有助于理解PD 的致病機理.

最近許多學(xué)者的研究更進一步強調(diào)了突觸傳輸時滯和不同神經(jīng)元集群之間的突觸連接強度會引起異常β 振蕩的產(chǎn)生[13,23,28],陳亞倩等[29]通過皮質(zhì)-基底神經(jīng)節(jié)模型研究了病理性β 振蕩的產(chǎn)生與抑制機理,發(fā)現(xiàn)突觸傳輸時滯的增加能夠?qū)е孪到y(tǒng)發(fā)生超臨界Hopf 分岔.此外,研究還發(fā)現(xiàn),不同神經(jīng)元集群的突觸連接強度會引起β 振蕩的產(chǎn)生.例如,從STN 到GPe 的強興奮連接能夠誘導(dǎo)GPe 產(chǎn)生β 振蕩.以上研究表明,突觸傳輸時滯和不同神經(jīng)元集群之間的突觸連接強度在PD 患者的健康狀態(tài)和帕金森病狀態(tài)之間起到了積極的調(diào)節(jié)作用.

因此,本文主要分為三個部分: 首先,使用突觸傳遞時滯作為分岔參數(shù),通過分岔理論分析了病理性β 振蕩的Hopf 分岔產(chǎn)生機制;其次,通過數(shù)值模擬闡述了突觸傳輸時滯和突觸連接強度能夠誘導(dǎo)異常β 振蕩的產(chǎn)生;最后,研究與GPi 有關(guān)的參數(shù)對其產(chǎn)生振蕩的影響.希望本文的研究結(jié)果能夠?qū)ε两鹕〉闹委煼椒ㄌ峁┮恍﹨⒖?

1 模型與方法

1.1 模型描述

基于Pavlides 等[20]提出的CTX-STN-GPe 模型,本文構(gòu)造了圖1 的CTX-STN-GPe-GPi 生理學(xué)共振模型.模型由兩部分神經(jīng)回路組成,其中,第一部分為CTX 的興奮性神經(jīng)元集群(excitatory pyramidal neuron,E)和抑制性神經(jīng)元集群(inhibitory pyramidal neurons,I),第二部分為STN-GPe-GPi 組成的神經(jīng)回路.此外,E 通過神經(jīng)元軸突將興奮性谷氨酸能傳遞到STN 神經(jīng)元集群,進一步地,STN 將興奮信號傳入GPe 和GPi,同時 STN 還接收到 GPe 的抑制性GABA 能輸入.紋狀體(striatum,STR) 將抑制性GABA 能傳給GPe 和GPi,同時GPi 也接受來自GPe 的抑制信號.

圖1 所示的模型由式(1)描述其動力學(xué)行為

圖1 皮質(zhì)-丘腦底核-蒼白球外側(cè)-蒼白球內(nèi)側(cè)共振模型示意圖Fig.1 The schematic diagram of Cortex-STN-GPe-GPi

其中S(t),G1(t),G2(t),E(t)和I(t) 分別表示STN,GPe,GPi,皮質(zhì)興奮性神經(jīng)元集群和抑制性神經(jīng)元集群的放電率.τi(i=S,G1,G2,E,I)表示神經(jīng)元集群i的膜時間常數(shù).Tij和 ωij分別表示神經(jīng)元集群i和j之間的突觸傳遞時滯和連接強度.C表示對E 的興奮性輸入常量,STR表示紋狀體的外界輸入.

FX(X=S,G1,G2,E,I)表示激活函數(shù),滿足下式[30]

其中,BX,MX為每個神經(jīng)元集群的基礎(chǔ)放電率和最大放電率.激活函數(shù)以及激活函數(shù)導(dǎo)數(shù)的圖像如圖2所示.本文中的所有參數(shù),除突觸傳遞時滯和突觸連接權(quán)重之外,都列在表1 中,更多細節(jié)可以在參考文獻[20,31-32]中查閱.

圖2 (a) 激活函數(shù)曲線圖和(b) 激活函數(shù)導(dǎo)數(shù)曲線圖Fig.2 (a) Curves of the activation functions and (b) derivatives curves for activation functions

表1 本文參數(shù)表Table 1 The parameter values used in this paper

1.2 Hopf 分岔的存在性

由于該模型的復(fù)雜性,為了簡化,假設(shè)膜時間常數(shù)和突觸傳遞時滯都相同,記為 τ和T,令 τ=10 ms.因此,可以得到如下動力學(xué)公式

所以系統(tǒng)(3)所對應(yīng)的特征方程為

其中

下面根據(jù)T值的不同分兩種情況討論.

情形1:T=0

當(dāng)T=0 時,式(5)可化為

其中

由Routh-Hurwitz 判據(jù)[33]可得

定理1.1系統(tǒng)(3)在平衡點P處是局部漸進穩(wěn)定的,當(dāng)且僅當(dāng)式(8)的所有根具有負實部,即式(8)的系數(shù)滿足下列條件

另一方面,如果式(5) 有一對復(fù)共軛純虛根λ=±iω(ω ＞0),則系統(tǒng)(3) 在平衡點處的穩(wěn)定性會發(fā)生變化,并在某些條件下會從原點處分岔出一系列小振幅的周期解,即在平衡點處發(fā)生了Hopf 分岔現(xiàn)象.因此,需要找出使h1(λ)和h2(λ) 有一對純虛根的條件.

情形2:T＞0

當(dāng)T＞0時,由式(6)可知,h1(λ) 有一個負實根λ=-1/τ,于是轉(zhuǎn)而研究

根的情況.(1)假設(shè) λ=iω+(ω+＞0)是h3(λ) 的根,可以得到

分離式(10)的實部和虛部,有

因此,把上述式(11) 和式(12) 的兩端平方相加,可得

令z=(ω+)2,則

如果條件

成立,那么式(10)的所有根都具有負實部,因此系統(tǒng)(3)的平衡點 (0,0,0,0,0) 是漸進穩(wěn)定的.

于是

(2) 假設(shè) λ=iω-(ω-＞0) 是h2(λ) 的根當(dāng)且僅當(dāng)λ=iω-滿足以下條件

顯然,同(1)推導(dǎo),可以得到

下面,驗證相應(yīng)的橫截條件.先對式(9)的兩邊關(guān)于突觸傳輸時滯T求導(dǎo),有

引理1.1對于超越方程[34-35]

當(dāng) (τ1,τ2,···,τm) 變動時,當(dāng)且僅當(dāng)有零根出現(xiàn)在虛軸或者穿過虛軸時,在右半平面的零點重數(shù)之和才發(fā)生變化.

由以上分析以及Ruan 和Wei 的推理[32],可以得到定理1.2.

定理1.2對于系統(tǒng)(3),假設(shè)條件H(1),H(3) 和H(4)滿足,則如下三條結(jié)論成立:

(1)當(dāng)T＜T0時,系統(tǒng)(3)在平衡點P處是漸進穩(wěn)定的;

(2) 當(dāng)T=T0時,系統(tǒng)(3) 在平衡點P處產(chǎn)生Hopf 分岔;

(3)當(dāng)T＞T0時,系統(tǒng)(3)在平衡點P處是不穩(wěn)定的.

2 結(jié)果與分析

2.1 突觸傳輸時滯誘導(dǎo)異常beta 振蕩

為了闡明突觸傳輸時滯和突觸連接強度對帕金森病異常β 振蕩產(chǎn)生的影響,本文借助MATLAB軟件對模型進行數(shù)值模擬.模型中使用的參數(shù)均與參考文獻[34]一致.即 ωCS=6.6 ,=3.22 ,=2.56,=2.56 ,=0.9 ,ωCC=3.08 .

越來越多的研究表明,突觸傳輸時滯在誘導(dǎo)帕金森病異常β 振蕩的產(chǎn)生中起著關(guān)鍵作用[29,36,37]根據(jù)1.2 節(jié)的分析推導(dǎo),得到分岔點T0=2.2 ms..

本節(jié)對CTX-STN-GPe 生理學(xué)共振模型進行了分岔分析并繪制了共振模型放電速率的時間歷程圖,相圖以及STN 放電速率關(guān)于傳輸時滯的分岔圖,見圖3(a)～圖3(d).

圖3(a)是關(guān)于突觸傳輸時滯的分岔圖,其中,藍色的點線表示穩(wěn)定的狀態(tài),紅色的點線表示不穩(wěn)定的狀態(tài).HB 表示Hopf 分岔點.從圖中可知,當(dāng)0 ＜T＜2.2ms 時共振模型表現(xiàn)出唯一的穩(wěn)定平衡點(對應(yīng)健康狀態(tài)),在任何初值狀態(tài)下STN 的放電速率都收斂于一個恒定的值.當(dāng)T=2.2 ms 時,共振模型經(jīng)過HB 點,產(chǎn)生極限環(huán)振蕩.當(dāng)T＞2.2 ms 時系統(tǒng)出現(xiàn)不穩(wěn)定的平衡點,產(chǎn)生振蕩(對應(yīng)帕金森病狀態(tài)).可見,共振模型在T0=2.2 ms 處發(fā)生了Hopf 分岔,這與1.2 節(jié)的理論結(jié)果一致.兩種動力學(xué)狀態(tài)所對應(yīng)時間序列圖以及相圖,分別見圖3(b)～圖3(d).在相圖中,紅色曲線代表健康狀態(tài),藍色曲線代表帕金森病狀態(tài).由此可見,改變突觸傳輸時滯,可以誘導(dǎo)異常β 振蕩的產(chǎn)生.

圖3 (a) STN 放電速率關(guān)于傳輸時滯的分岔圖.共振模型時間歷程圖: (b) 健康狀態(tài),(c)帕金森病狀態(tài)和(d)共振模型相圖Fig.3 (a) Bifurcation diagram of the firing rate of the STN population STN against the transmission delay.Time series of the resonance model:(b) the healthy state,(c) the Parkinson state and (d) phase diagram of the resonance model

2.2 突觸連接強度誘導(dǎo)異常beta 振蕩

為了研究突觸連接強度的改變對異常β 振蕩的產(chǎn)生和振蕩的振幅的影響,本節(jié)分別討論了突觸連接強度 ωCS單獨變化對STN 神經(jīng)元集群動力學(xué)行為的影響以及突觸連接強度 ωCS和ωG1S同時變化對STN 神經(jīng)元集群動力學(xué)行為的影響.

首先探討興奮性神經(jīng)元集群E 到STN 的興奮輸入連接強度 ωCS對STN 動力學(xué)行為的影響.ωCS關(guān)于STN 放電速率的分岔圖如圖4(a)所示,由分岔圖可知,當(dāng) 0＜ωCS＜80 時,STN 的動力學(xué)行為表現(xiàn)為振蕩狀態(tài),當(dāng) ωCS＞80 時,STN 從之前的振蕩狀態(tài)轉(zhuǎn)變?yōu)閱畏€(wěn)態(tài).因此,只要 ωCS在大于0 的參數(shù)區(qū)域內(nèi)足夠大,可以有效地阻止異常β 振蕩的產(chǎn)生.

圖4 (a)興奮性神經(jīng)元集群對STN 放電率的影響. ωCS和 ωGeS對STN (b) 振幅和 (c)放電率的影響Fig.4 (a) The effect of excitatory pyramidal neuron on STN firing rate.The effect of ωCSand ωG1Son (b) amplitude and (c) firing rate of STN

為了更全面地探討突觸連接強度對STN 產(chǎn)生β 振蕩的影響,下面討論雙參數(shù) ωCS和 ωG1S變化下STN 的振幅變化,進而分析 ωCS和 ωG1S的整體改變對STN 動力學(xué)行為的影響.圖4(b)給出了雙參數(shù)ωCS和 ωG1S變化對STN 振幅的影響.從圖4(b)可以看出,當(dāng)參數(shù) ωCS和 ωG1S取自藍色區(qū)域時,STN 振蕩消失,系統(tǒng)達到健康狀態(tài),當(dāng)參數(shù) ωCS和 ωG1S取在黃色區(qū)域內(nèi),STN 產(chǎn)生振蕩,對應(yīng)為帕金森病狀態(tài),振蕩的振幅在黃色區(qū)域內(nèi)較大.為了進一步刻畫這種現(xiàn)象,下面分別選取一些典型的參數(shù)來模擬STN的動力學(xué)行為.當(dāng) ωCS與 ωG1S兩者的作用逐漸增強時,STN 容易發(fā)生振蕩,且振幅也逐漸增大.當(dāng)ωG1S=3.22,ωCS=3 時,系統(tǒng)產(chǎn)生較小的振蕩,當(dāng)ωG1S,ωCS逐漸增加,比如當(dāng)ωG1S=8,ωCS=35 時,系統(tǒng)產(chǎn)生高振幅的振蕩,見圖4(c).

2.3 GPi 的相關(guān)參數(shù)對其振蕩的影響

由于GPi 受GPe 突觸連接強度和突觸傳輸時滯的共同作用,故本節(jié)通過考慮振蕩的振幅,來討論GPe 突觸連接強度和突觸傳輸時滯對GPi 振蕩的影響.從圖5(a)中可以看出,參數(shù)區(qū)域分為三部分,當(dāng)參數(shù)在A 區(qū)域變化時,GPi 振蕩,當(dāng)參數(shù)在B 區(qū)域時,振蕩的振幅減小,而在區(qū)域C 時,GPi 振蕩的振幅較小甚至消失.因此,當(dāng)較小的GPe 突觸連接強度和較大的突觸傳輸時滯共同作用時,更容易使得GPi 發(fā)生振蕩,且振幅越來越大.圖5(b)為圖5(a)對應(yīng)的時間序列圖,以便更簡潔明了的通過振幅觀察參數(shù)對GPi 振蕩的影響.從圖5(b)可知較小的和較大的會促使GPi 產(chǎn)生振蕩.

圖5 和 對GPi (a)振幅和(b)放電率的影響Fig.5 The effect of and on (a) amplitude and (b) firing rate of GPi

圖5 和 對GPi (a)振幅和(b)放電率的影響 (續(xù))Fig.5 The effect of and on (a) amplitude and (b) firing rate of GPi (continued)

3 結(jié)論

本文結(jié)合PD 相關(guān)生理理論,分岔理論,從建模分析的角度出發(fā),對皮質(zhì)-丘腦底核-蒼白球外側(cè)共振模型進行拓展,通過分岔分析和數(shù)值模擬揭示了突觸傳輸時滯和突觸連接強度對PD 的動力學(xué)機制的影響.結(jié)果表明,突觸傳輸時滯的增加可以誘導(dǎo)異常β 振蕩的產(chǎn)生.其次,發(fā)現(xiàn)丘腦底核同時受到興奮性神經(jīng)元集群和蒼白球外側(cè)較強的促進作用時更容易發(fā)生振蕩.最后,通過分析β 振蕩的振幅對GPi 相關(guān)參數(shù)的依賴性,發(fā)現(xiàn)較小的GPe 突觸連接強度和較大的突觸傳輸時滯使得GPi 產(chǎn)生振蕩.希望本文的研究結(jié)果能夠為研究帕金森病的機理提供一些參考.