丁學(xué)利 賈冰 李玉葉

1) (阜陽職業(yè)技術(shù)學(xué)院基礎(chǔ)教學(xué)部,阜陽 236031)

2) (復(fù)旦大學(xué)生命科學(xué)學(xué)院生理與生物物理系,上海 200438)

3) (赤峰學(xué)院數(shù)學(xué)與統(tǒng)計學(xué)院,赤峰 024000)

1 引 言

生物系統(tǒng)是復(fù)雜的非線性系統(tǒng),一直受到非線性科學(xué)的重點關(guān)注.非線性動力學(xué)在識別神經(jīng)系統(tǒng)中各種復(fù)雜的放電或振蕩模式中起著重要的作用[1?4].神經(jīng)元的電活動包括靜息和放電等.靜息對應(yīng)于動力學(xué)中的穩(wěn)定平衡點.而膜電位處于大幅值振蕩狀態(tài)稱為放電,放電包括峰放電和簇放電等,這與動力學(xué)中穩(wěn)定的極限環(huán)相對應(yīng).神經(jīng)元從靜息到放電以及從放電到靜息的變化,可以用動力學(xué)中的平衡點分岔進行描述[1?4].已知的平衡點分岔包括不變圓上的鞍結(jié)分岔和Hopf分岔等.不變圓上的鞍結(jié)分岔和Hopf分岔會導(dǎo)致神經(jīng)元產(chǎn)生兩類不同的興奮性類型[4?6],即I型興奮性和II型興奮性.隨著去極化電流的增加,I型興奮性的放電頻率從零逐步增加到某個非零值,對應(yīng)于不變圓上的鞍結(jié)分岔; 而Ⅱ型興奮性的放電頻率則從某個非零值開始,對應(yīng)于Hopf分岔.I型和II型興奮性神經(jīng)元對應(yīng)的外界興奮脈沖擾動的相位響應(yīng)曲線(phase response curve,PRC)[7?10]不同.脈沖擾動下的放電相位提前,也就是PRC的值全正,對應(yīng)I型興奮性; 擾動后放電的相位依賴于刺激脈沖相位,可能提前,也可能延遲,也就是PRC的值部分為正、部分為負,對應(yīng)于II型興奮性.這分別對應(yīng)I,II型PRC,在經(jīng)典的二維Morris-Lecar (ML)模型[7,10]中得到了廣泛的研究.

除了脈沖式的調(diào)控,神經(jīng)元電活動還受到噪聲和突觸介導(dǎo)的耦合電流等諸多因素的影響.比如,在很多類的神經(jīng)元,其放電模式的變化會受到自突觸的調(diào)控[11?15].自突觸,即神經(jīng)元通過突觸與自己建立聯(lián)系,這在大腦皮層、新皮層、視覺皮層、紋狀體、海馬和小腦等神經(jīng)系統(tǒng)中都有解剖學(xué)觀察[11?15].通過生物學(xué)實驗發(fā)現(xiàn),抑制性自突觸是存在的[15?19],實驗結(jié)合理論模型的工作確認了慢抑制性突觸的存在性[15,20],可引起精確放電[15].精確放電可用多次相同刺激試驗中依次產(chǎn)生的動作電位的產(chǎn)生時間的標(biāo)準(zhǔn)差(standard deviation,STD)來度量,或通過放電的峰峰間隔(interspike intervals,ISIs)的變異系數(shù)(coefficient of variation,CV)來刻畫[15?20].CV值越小,動作電位產(chǎn)生的時刻越準(zhǔn)時,放電序列越有規(guī)律,意味著放電時刻越精確.抑制性自突觸誘發(fā)的精確放電發(fā)生在大腦皮層的新皮層,并被認為通過影響新皮層的同步行為參與信息處理、行為整合和認知功能.例如,Bacci等[15,16]研究發(fā)現(xiàn)在新皮層上通過抑制性中間神經(jīng)元的自突觸傳輸可增加放電的精確性.當(dāng)Theta和Gamma振蕩出現(xiàn)在海馬和新皮層時,海馬中間神經(jīng)元表現(xiàn)出精確的放電[17,18].Zhao等[19]通過仿真發(fā)現(xiàn)超極化激活非特異性陽離子電流誘導(dǎo)神經(jīng)元模型產(chǎn)生精確放電.數(shù)值仿真還發(fā)現(xiàn),興奮性自突觸和抑制性自突觸對神經(jīng)元放電頻率產(chǎn)生不同的影響[21?27].因為實驗發(fā)現(xiàn)慢抑制性自突觸的存在,即抑制性耦合電流再作用回神經(jīng)元是需要時間的,可以用時滯來模擬這一慢特征.例如,當(dāng)引入時間延遲時,自突觸可以幫助增強神經(jīng)元網(wǎng)絡(luò)的集群活動和誘發(fā)靶波、螺旋波及同步行為[28,29]; 帶有時間延遲的抑制性自突觸可以抑制神經(jīng)元的混沌行為[25]、產(chǎn)生相干共振[30]等.此外,大量的電生理實驗中采用動態(tài)鉗(dynamic clamp)技術(shù)來調(diào)控突觸或反饋電流的幅值或快慢來調(diào)控神經(jīng)元或網(wǎng)絡(luò)的電活動或時空行為[31?33].其中,自反饋調(diào)控也是一種潛在的調(diào)控神經(jīng)元電活動的手段.與現(xiàn)實神經(jīng)系統(tǒng)的突觸的快慢相對不能方便調(diào)控不同,動態(tài)鉗能夠方便地調(diào)控突觸的快慢.例如,當(dāng)抑制性突觸的時間常數(shù)由5 ms調(diào)整到600 ms,交互抑制耦合的兩個神經(jīng)元可以從反相同步達到同相同步[31].因此,從理論上研究較長時滯的自反饋或自突觸的調(diào)控作用對于認識現(xiàn)實的自突觸的特性或者提出能用于動態(tài)鉗實驗的調(diào)控措施有指導(dǎo)意義.

一般情況下,興奮性作用增加放電,抑制性作用降低放電[34],表現(xiàn)在放電頻率或峰峰間隔ISIs的變化.但是,最近的研究發(fā)現(xiàn)抑制性作用可以增強放電,如抑制性自突觸可引起簇放電的增加[27,35,36]; 而興奮性作用也會降低放電個數(shù)[37].實際上,目前的峰放電的II型PRC曲線都是針對興奮性脈沖的,其特征是興奮性脈沖在一些作用相位下可以使得放電提前,在另外的相位可以使得放電延后.抑制性脈沖激勵下的II型相位響應(yīng)會表現(xiàn)出什么特征,目前并不清晰.抑制性的脈沖激勵是否能使得放電提前也就是增強放電? 假若能,那么抑制性自反饋或自突觸會不會引起峰放電的增加呢? 這都是亟待回答的非線性動力學(xué)問題.此外,抑制性自反饋或自突觸誘發(fā)的精確放電是否能夠在理論模型中獲得解釋,也是對神經(jīng)科學(xué)問題有重要意義的.

因此,本文通過理論模型研究神經(jīng)峰放電在抑制性作用下的PRC,發(fā)現(xiàn)某些相位下放電提前; 進一步研究了峰放電頻率在抑制性自反饋作用下會提高,可以用PRC進行解釋; 最后,較長時滯的自反饋能夠誘發(fā)放電精確性提高,與實驗發(fā)現(xiàn)的慢抑制性自突觸誘發(fā)放電精確性增加相一致.研究結(jié)果揭示了抑制性作用的新功能并給出了非線性動力學(xué)解釋,給出了調(diào)控放電頻率增加的新措施,抑制性反饋而不是興奮性反饋; 這些結(jié)果有助于認識現(xiàn)實神經(jīng)系統(tǒng)的抑制性自突觸的潛在功能.

2 神經(jīng)元模型

2.1 ML模型和自反饋模型

ML模型和自反饋模型[10,20,23]如下:

其中C是膜電容,V(t)和w(t)分別是膜電壓和延遲整流K+電流的激活變量.(1a)式等號右邊的前三項分別表示Ca2+電流、電壓門控延遲整流K+電流和漏電流.參數(shù)gCa,gK和gL分別是鈣電流、鉀電流和漏電流的最大電導(dǎo)率.VCa,VK和VL分別是鈣電流、鉀電流和漏電流的反轉(zhuǎn)電位.I表示外激勵的(去極化)電流;Iaut為自反饋電流,可以對應(yīng)現(xiàn)實神經(jīng)系統(tǒng)的自突觸電流.Iaut的表達式如下:

其中Γ(t?τ)=1/{1+exp[?(V(t?τ)?θsyn)/λ]}.這里qsyn=–20 mV,l=1,t表示時間延遲,gaut是自反饋強度,可以對應(yīng)自突觸的電導(dǎo).當(dāng)Vsyn=–60 mV時,Iaut是抑制性自反饋電流,可以對應(yīng)抑制性自突觸電流; 當(dāng)Vsyn=50 mV時,Iaut是興奮性自反饋電流,可以對應(yīng)興奮性自突觸電流.相關(guān)的函數(shù)的表達式如下:

ML模型產(chǎn)生II型興奮的參數(shù)設(shè)置如下:C=5 mF/cm2,VK=–80 mV,VCa=120 mV,VL=–60 mV,gK=8 mS/cm2,gCa=4 mS/cm2,gL=2 mS/cm2,V1=–1.2 mV,V2=18 mV,V3=4 mV,V4=17.4 mV,f=0.066667.

2.2 PRC的定義

PRC,又稱相位重置曲線[7?10],用于描述神經(jīng)元振子做周期性放電時,外界激勵對峰放電發(fā)放時間或相位偏移的影響.假設(shè)神經(jīng)元模型存在一個周期為T0的穩(wěn)定極限環(huán).外界激勵作用于以峰值時刻開始計時的時刻ts,刺激后的放電的峰值出現(xiàn)在時刻T1,如圖1(a)所示,則刺激相位j=ts/T0,放電的相位為T1/T0,因此,相位重置D(j)定義為

D(j)與j的關(guān)系曲線即相位重置曲線(或PRC).此外,為了方便,也可以定義

本文使用外激勵電流為幅值為A,寬度為d的方波脈沖.例如,當(dāng)I=45.5 μA/cm2時,ML模型的放電周期T0≈ 56.37 ms,如圖1(a)中的點線所示; 在ts=40 ms處,受方波脈沖電流(幅值A(chǔ)=1.65 μA/cm2,寬度d=4.4 ms)激勵,如圖1(a)中的實線所示,方波脈沖擾動后的放電發(fā)生在T1≈52.3 ms.由(3)式可計算在ts=40 ms處的相位重置PR≈ 0.072,對應(yīng)j≈ 0.71的D(j) ≈ 0.072.若計算周期內(nèi)每一點的相位重置,則可得PRC,如圖1(b)所示.在圖1(b)中,大約在ts< 27 ms時,PRC取值為負,表明方波脈沖激勵使峰放電延后,而在ts> 27 ms時,PR取值為正,說明方波脈沖激勵使峰放電提前.這種相位有正有負的PRC即為II型PRC,對應(yīng)Hopf分岔附近的放電[7?10].

2.3 精確放電的度量

由于現(xiàn)實的神經(jīng)元放電受噪聲的影響,因此噪聲會影響神經(jīng)元放電的精確性.為了研究噪聲對神經(jīng)元放電的影響,在模型的方程(1a)的右端加入高斯白噪聲Inoise如下:

其中,x(t)表示高斯白噪聲,D表示噪聲強度.為了估計放電的規(guī)律性,假設(shè)放電序列中的尖峰時刻發(fā)生在ti(0 ≤i≤N),那么放電的峰峰間隔ISIs定義為ISIsi=ti+1–ti.

圖1 當(dāng)I=45.5 μA·cm–2,方波脈沖幅值A(chǔ)=1.65 μA·cm–2,寬度d=4.4 ms時,ML模型的放電序列和PRC (a) 方波脈沖刺激電流(短劃線)、沒有方波脈沖的放電(點線)和有方波脈沖的放電(實線); (b) PRCFig.1.Spiketrains and PRC of the ML model when I=45.5 μA·cm–2,A=1.65 μA·cm–2,and d=4.4 ms: (a) Square pulse disturbance current (dashed line),and spike trains without (dotted line) and with square pulse disturbance (solid line); (b) PRC.

本研究采用ISIs的CV來刻畫放電的精確性.通過如下公式計算CV:

3 仿真結(jié)果

3.1 ML模型的平衡點的亞臨界Hopf分岔

在沒有自突觸的作用下(Iaut=0),以I作為分岔參數(shù),ML模型的平衡點分岔見圖2(a).隨著分岔參數(shù)I的增大,ML模型在處發(fā)生亞臨界的Hopf分岔.通過該分岔點,穩(wěn)定焦點(粗實線)轉(zhuǎn)化為不穩(wěn)定的焦點(虛線),同時產(chǎn)生不穩(wěn)定的極限環(huán)(點線).不穩(wěn)定極限環(huán)和穩(wěn)定極限環(huán)(細實線)在處發(fā)生極限環(huán)的鞍結(jié)分岔.當(dāng)44.65 μA/cm245.23 μA/cm2時,系統(tǒng)只出現(xiàn)穩(wěn)定的極限環(huán),即單穩(wěn)狀態(tài).如當(dāng)I=45.5 μA/cm2時,可得到ML模型的單穩(wěn)放電序列,放電周期如圖2(b)所示.而當(dāng)I< 44.65 μA/cm2時,ML模型處于靜息狀態(tài),如圖2(c)(當(dāng)I=44 μA/cm2時的靜息狀態(tài))所示.在極限環(huán)穩(wěn)定的范圍內(nèi),隨著I的增大,放電周期ISIs(圖2(d)實線)單調(diào)減小,頻率(圖2(d)虛線)單調(diào)增加,但二者在單穩(wěn)放電區(qū)域內(nèi)變化不大,如圖2(d)所示.本文研究的是靠近分岔點的放電/單穩(wěn)極限環(huán).

3.2 抑制性脈沖作用下的II型PRC及推測

3.2.1 抑制性脈沖作用下的II型PRC

以前的研究表明,正向方波脈沖電流可使靠近Hopf分岔點的放電提前或延遲,即II型PRC,如圖1所示.本文研究負向方波脈沖電流作用于放電的不同相位的PRC,發(fā)現(xiàn)該曲線與正向方波脈沖引起的PRC的形狀關(guān)于橫坐標(biāo)(刺激時刻)軸大約是對稱的,其相位也是有正有負的.例如,當(dāng)I=45.5 μA/cm2時,ML模型峰放電靠近Hopf分岔,接受負向脈沖電流的刺激,刺激相位不同時,可使得放電提前或延遲,如圖3所示.圖3中紅色的點線表示沒有方波脈沖刺激,也就是原始的放電序列,虛線表示負向方波脈沖電流,實線表示負向方波脈沖擾動后的放電序列.取負向方波脈沖的幅值A(chǔ)=–0.6 μA/cm2,寬度d=4.9 ms,作用在ts=22 ms處,可得到負向方波脈沖激勵后的放電序列(實線在紅色的點線前)提前了,如圖3(a)所示.同樣地,增大負向方波脈沖幅值的絕對值,如A=–1.65 μA/cm2,d=4.8 ms,負向方波脈沖激勵后的放電序列(實線在紅色的點線前)也提前了,如圖3(b)所示.而當(dāng)負向方波脈沖作用在ts=40 ms處,負向脈沖作用后的峰放電(實線在紅色的點線后)延遲了,如圖3(c)和圖3(d)所示.刺激強度越大,放電的延遲越大,如圖3(d)比圖3(c)的刺激強度大,其延遲就越大.

圖2 (a) ML模型隨I的平衡點分岔; (b) 當(dāng)I=45.5 μA·cm–2時,ML模型的放電序列圖; (c) 當(dāng)I=44 μA·cm–2時,ML模型處于靜息狀態(tài); (d) ISIs和頻率隨I的變化(實線表示ISIs,虛線表示頻率)Fig.2.(a) Bifurcation of ML model with respect to I; (b) spike trains of ML model when I=45.5 μA·cm–2; (c) resting state of ML model when I=44 μA·cm–2; (d) the changes of ISIs (solid line) and frequency (dashed line) with respect to I.

利用圖3中負向方波脈沖的幅值和寬度可分別做出其對應(yīng)的PRC曲線,如圖4所示.圖4(a)中,當(dāng)負向方波脈沖的幅值A(chǔ)=–0.6 μA/cm2,寬度d=4.9 ms時,大約在ts< 27.2 ms時,PR取值為正,表明負向方波脈沖激勵使峰放電提前,而在ts> 27.2 ms時,PR取值為負,說明負向方波脈沖激勵使峰放電延遲,這與圖3(a)和圖3(c)相一致.類似地,當(dāng)負向方波脈沖的幅值A(chǔ)=–1.65 μA/cm2,寬度d=4.8 ms時,其PRC如圖4(b)所示.大約在ts< 27.4 ms時,PR取值為正,而在ts> 27.4 ms時,PR取值為負,這也較好地與圖3(b)和圖3(d)相吻合.圖4所示的負向脈沖PRC與圖1(b)所示的正向脈沖誘發(fā)的PRC的形狀關(guān)于橫坐標(biāo)(刺激時刻)軸大約是對稱的.刺激強度越大,放電的延遲越大,如圖4(b)比圖4(a)的刺激強度大,其延遲就越大.

3.2.2 基于II型PRC推測抑制性自反饋電流的作用

圖4所示的II型PRC有助于理解抑制性自突觸作用下神經(jīng)元電活動的動力學(xué)性質(zhì)和頻率變化.對于具有抑制性自突觸的神經(jīng)元,抑制性自突觸電流的表達式為Iaut=?gaut(V(t)?Vsyn)Γ(t?τ) ,其時間歷程如圖5中的虛線所示,在很大程度上與圖3中的激勵脈沖相似.自突觸的時滯t可以調(diào)節(jié)抑制性自突觸電流脈沖的作用相位,對應(yīng)于刺激相位ts; 而抑制性自突觸脈沖電流的幅值由自突觸強度gaut決定,對應(yīng)于刺激脈沖電流的幅值A(chǔ).因此,自突觸電流作用下的放電與PRC曲線有一定的對應(yīng)關(guān)系.合適的t值下,抑制性自突觸電流會引起動作電位間隔變小,也就是放電峰峰間期變小,放電頻率增加,即放電增強; 另外的t值下,抑制性自突觸電流會引起動作電位間隔變大,放電頻率降低,即放電減弱.不同的自突觸電導(dǎo)gaut,對應(yīng)于刺激強度A值的不同,放電ISIs或頻率的變化程度是不同的.

圖3 當(dāng)I=45.5 μA·cm–2時,負向方波脈沖電流(虛線)作用在不同相位的放電序列(實線)和無方波脈沖作用的放電序列(紅色的點線) (a) ts=22 ms,A=–0.6 μA·cm–2,d=4.9 ms; (b) ts=22 ms,A=–1.65 μA·cm–2,d=4.8 ms; (c) ts=40 ms,A=–0.6 μA·cm–2,d=4.9 ms; (d) ts=40 ms,A=–1.65 μA·cm–2,d=4.8 msFig.3.Spike trains induced by square pulse current applied at different phases when I=45.5 μA·cm–2.The spike trains (solid line)influenced by negative square pulse current (dashed line) and the trains (red dotted line) without negative square pulse current.(a) ts=22 ms,A=–0.6 μA·cm–2,d=4.9 ms; (b) ts=22 ms,A=–1.65 μA·cm–2,d=4.8 ms; (c) ts=40 ms,A=–0.6 μA·cm–2,d=4.9 ms; (d) ts=40 ms,A=–1.65 μA·cm–2,d=4.8 ms.

圖4 當(dāng)I=45.5 μA·cm–2時,ML模型在Hopf分岔點附近的負向脈沖刺激誘發(fā)的PRC (a) A=–0.6 μA·cm–2,d=4.9 ms;(b) A=–1.65 μA·cm–2,d=4.8 msFig.4.PRC induced by negative square pulse current near the Hopf bifurcation point in the ML model when I=45.5 μA·cm–2:(a) A=–0.6 μA·cm–2,d=4.9 ms; (b) A=–1.65 μA·cm–2,d=4.8 ms.

圖5 當(dāng)I=45.5 μA·cm–2,gaut=0.04 mS·cm–2,具有抑制性自突觸ML模型的放電模式(實線)與抑制性自突觸電流(短劃線)(a) t=0 mS; (b) t=10 mS; (c) t=20 mS; (d) t=30 mS; (e) t=40 mS; (f) t=50 mSFig.5.Inhibitory autapse current (dashed line) and spike trains (solid line) of the ML model with inhibitory autapse when I=45.5 μA·cm–2 and gaut=0.04 mS·cm–2: (a) t=0 mS; (b) t=10 mS; (c) t=20 mS; (d) t=30 mS; (e) t=40 mS; (f) t=50 mS.

3.3 抑制性自反饋對神經(jīng)元模型峰放電模式的影響

3.3.1 抑制性自反饋的時滯t對ML模型峰放電的影響

具有自反饋的ML模型的數(shù)值模擬結(jié)果驗證了上述的推測結(jié)果.下面以I=45.5 μA/cm2,gaut=0.04 mS/cm2為例介紹ML模型在抑制性自反饋電流作用下,隨著時滯t的增大放電模式的變化.虛線表示抑制性自反饋電流,實線表示神經(jīng)元模型放電序列,如圖5所示.圖5中,抑制性自反饋電流的幅值和寬度大約分別與圖3中抑制性脈沖的幅值和寬度相對應(yīng).當(dāng)t=0 ms時,放電周期為56.48 ms,大于單神經(jīng)元的放電周期; 當(dāng)t=10和20 ms時,放電周期分別為56.31和55.95 ms,小于單神經(jīng)元的放電周期,即峰放電提前,這與圖4(b)中大約在ts< 27.4 ms時,PR取值為正相吻合;當(dāng)t=30,40和50 ms時,放電周期分別為57.14,63.95和65.41 ms,大于單神經(jīng)元的放電周期,即峰放電延遲,這與圖4(b)中大約在ts> 27.4 ms時,PR取值為負相一致.

3.3.2 不同耦合強度gaut下的歸一化ISIs隨時滯t的變化

在抑制性自反饋電流作用下,神經(jīng)元模型可能不會立即引起尖峰,但可以改變ISIs,即放電周期.利用時滯t調(diào)節(jié)抑制性自反饋電流作用在神經(jīng)元模型放電周期的不同相位,可以測量ISIs隨時滯t的變化.為了更好地揭示神經(jīng)元模型放電周期如何隨時滯t變化,將ISIs進行歸一化處理,即將ISIs除以沒有自反饋電流的放電周期T≈ 56.37 ms.若ISIs/T> 1,則放電周期延長; 若ISIs/T< 1,則放電周期縮短.例如,當(dāng)gaut=0.01 mS/cm2時,歸一化的ISIs隨時滯t的變化如圖6所示.圖6(a)中,粗實線表示歸一化的ISIs隨時滯t的變化,細實線表示其放電頻率隨時滯t的變化.圖6(a)和圖6(b)中ISIs/T=1對應(yīng)的t值大約分別是27.2 ms和27.4 ms.當(dāng)時滯t< 27.2 ms時,ISIs/T幾乎都小于1,這與圖4(a)中PR值大于0相對應(yīng),說明放電周期縮短,其對應(yīng)的放電頻率增加.抑制性自反饋可以誘發(fā)放電頻率增加是一個不同于抑制性作用降低放電頻率的傳統(tǒng)觀念的新觀點.而在時滯t> 27.2 ms時,ISIs/T幾乎都大于1,這與圖4(a)中PR值小于0相對應(yīng),說明放電周期延遲,其對應(yīng)的放電頻率降低.同樣地,當(dāng)gaut=0.04 mS/cm2時,也可以得到類似的結(jié)論,這與圖4(b)中的PRC相對應(yīng),如圖6(b)所示.圖6(a)中頻率的變化范圍大約從17.06—17.78 Hz,圖6(b)中頻率的變化范圍大約從13.9—17.9 Hz; 刺激強度越大,對PR的影響越大.

圖6 不同耦合強度gaut下歸一化的ISIs (粗實線)和放電頻率(細實線)隨時滯t的變化 (a) gaut=0.01 mS·cm–2; (b) gaut=0.04 mS·cm–2Fig.6.Change of normalized ISIs (boldsolid line) and firing frequency (thin solid line) with respect to time delay t: (a) gaut=0.01 mS·cm–2; (b) gaut=0.04 mS·cm–2.

3.3.3 在噪聲擾動下的ISIs變化和PRC特征

當(dāng)噪聲強度D=0.5 μA/cm2時,平均ISIs對時滯t和耦合強度gaut的依賴關(guān)系如圖7所示.紅色區(qū)域表示平均ISIs小于56.37 ms,即平均ISIs降低; 綠色區(qū)域表示平均ISIs大于56.37 ms,即平均ISIs增加.說明在t約小于31 ms時放電是提前的,t約大于31 ms時放電是延遲的.與圖3所示的確定性模型的結(jié)果相一致,說明PRC曲線的特征在噪聲作用下有一定的魯棒性.

圖7 當(dāng)噪聲強度D=0.5 μA·cm–2時,平均ISIs對時滯t和耦合強度gaut的依賴關(guān)系Fig.7.Dependence of average ISIs on time delay t and coupling strength gaut when D=0.5 μA·cm–2.

3.4 抑制性自反饋引起精確放電

生物學(xué)實驗中將神經(jīng)元發(fā)放的精確性定義為動作電位峰值時刻的抖動(jitter),量化為多次相同刺激實驗中順序產(chǎn)生的動作電位的尖峰時間的標(biāo)準(zhǔn)差STD[15,16],也可量化為動作電位平均ISIs的CV[15,16].放電精確性衡量神經(jīng)系統(tǒng)的電活動在存在外界擾動比如噪聲的情況下的響應(yīng)特性.通過神經(jīng)元模型的平均ISIs的標(biāo)準(zhǔn)差和CV來反映抑制性自反饋引起的精確放電現(xiàn)象.

當(dāng)I=45.5 μA/cm2,噪聲強度D=0.5 μA/cm2時,放電的平均ISIs的標(biāo)準(zhǔn)差STD在 (t,gaut)平面的分布,也就是對時滯和耦合強度的依賴關(guān)系,如圖8(a)所示.平均ISIs的標(biāo)準(zhǔn)差有以下幾個特點: 首先,STD的低值主要出現(xiàn)在(t,gaut)平面中間的區(qū)域,即放電精確性高的產(chǎn)生區(qū)域; 如在gaut=0.25—1 mS/cm2和t=25—31 ms,STD的變化范圍大約是3—6.其次,當(dāng)gaut< 0.25 mS/cm2,或t< 25 ms,或t> 31 ms時,STD的值都較大(大于6).最后,在延遲較小時(t=1—2 ms),隨著耦合強度的增大(gaut> 0.35 mS/cm2),STD出現(xiàn)較大的波動.

當(dāng)I=45.5 μA/cm2,噪聲強度D=0.5 μA/cm2時,放電的平均ISIs的CV在(t,gaut)平面的分布,也就是對時滯和耦合強度的依賴關(guān)系,如圖8(b)所示.CV的值越小,表明神經(jīng)元模型放電具有較高的放電精確性.從圖8(b)可以看出,時滯t大約在25—35 ms的范圍內(nèi)、及在耦合強度gaut>0.2 mS/cm2且時滯t大約在35—48 ms的范圍內(nèi)神經(jīng)元模型放電具有較高的放電精確性.上述結(jié)果說明精確放電發(fā)生在時滯相對較大的情況下,這與實驗發(fā)現(xiàn)抑制性自突觸是慢突觸有一定的對應(yīng)性[15].

圖8 當(dāng)I=45.5 μA·cm–2,D=0.5 μA·cm–2時,(a) ISIs的STD對時滯t和耦合強度gaut的依賴關(guān)系,(b) ISIs的CV對時滯t和耦合強度gaut的依賴關(guān)系Fig.8.(a) Dependence of standard deviation of ISIs (STD) on time delay t and coupling strength gaut; (b) the dependence of coefficient of variation of ISIs (CV) on delay t and coupling strength gaut.The parameter values are I=45.5 μA·cm–2 and D=0.5 μA·cm–2.

在圖8(b)中,若固定時滯t,可得到ISIs的CV隨著耦合強度gaut的變化,如圖9所示.當(dāng)D=0.5 μA/cm2時,在時滯t=27 ms,隨著耦合強度gaut的增大,CV的值逐漸變小,說明放電越來越精確,如圖9(a)所示.在時滯t=30,40,46 ms時,耦合強度越大,放電精確性越高,如圖9(b)—(d).因此當(dāng)固定時滯t時,隨著耦合強度gaut的增大,峰放電精確性越高,這與實驗現(xiàn)象相一致[15].

與圖9類似,若固定耦合強度gaut,可得到ISIs的CV隨著時滯t的變化,如圖10所示.在圖10(a)中,隨著時滯t的增大,CV的值先變小后逐漸變大,大約在t=28 ms時,CV達到最小值.同樣地,在gaut=0.61 mS/cm2,t=30 ms時(圖10(b)),CV達到最小值.

從神經(jīng)元模型的放電序列圖也可以看出時滯的變化對精確放電的影響.例如當(dāng)I=45.5 μA/cm2,D=0.5 μA/cm2,gaut=0.61 mS/cm2時,在圖11中展示了不同時滯下精確放電的差異.當(dāng)時滯t=1,10,30,50 ms時,其變異系數(shù)CV的值約分別為0.72,0.16,0.06和0.22,即在t=30 ms時,CV的值最小.說明在t=30 ms時抑制性自反饋可較好地使神經(jīng)元模型放電精確,如圖11(c)所示.

圖9 固定時滯t在不同水平下,ISIs的CV隨著耦合強度gaut的變化 (a) t=27 ms; (b) t=30 ms; (c) t=40 ms; (d) t=46 msFig.9.Changes of coefficient of variation (CV) of ISIs with respect to coupling strength gaut when time delay t is fixed at different values: (a) t=27 ms; (b) t=30 ms; (c) t=40 ms; (d) t=46 ms.

圖10 固定耦合強度gaut在不同水平下,ISIs的變異系數(shù)CV隨著時滯t的變化 (a) gaut=0.31 mS·cm–2; (b) gaut=0.61 mS·cm–2Fig.10.Changes of coefficient of variation (CV) of ISIs with respect to time delay t when coupling strength gaut is fixed at different levels: (a) gaut=0.31 mS·cm–2; (b) gaut=0.61 mS·cm–2.

圖11 當(dāng)I=45.5 μA·cm–2,D=0.5 μA·cm–2,耦合強度gaut=0.61 mS·cm–2時,時滯t對神經(jīng)元模型的精確放電的影響 (a) t=1 ms; (b) t=10 ms; (c) t=30 ms; (d) t=50 msFig.11.Effect of time delay t on spike-timing precision of neuron model when I=45.5 μA·cm–2,D=0.5 μA·cm–2,and gaut=0.61 mS·cm–2: (a) t=1 ms; (b) t=10 ms; (c) t=30 ms; (d) t=50 ms.

4 結(jié) 論

本文研究了抑制性自反饋增強放電活動這一新現(xiàn)象,并可以用PRC進行解釋,提示了自突觸的新的潛在功能,主要結(jié)論如下.

1)抑制性自反饋的負反饋電流能引起峰放電頻率增加,在噪聲作用下這一特性仍舊能夠維持,這是不同于傳統(tǒng)結(jié)果—抑制性作用引起頻率降低的新發(fā)現(xiàn),也是反常的神經(jīng)動力學(xué)和神經(jīng)電活動現(xiàn)象.以前抑制性刺激引起的反常神經(jīng)電活動研究包括抑制性刺激誘發(fā)靜息變?yōu)榉烹奫26]、抑制性耦合引起簇放電頻率增強[27]和抑制性耦合增強同相同步[35,36],而本文則揭示了抑制性作用增強峰放電的放電頻率,是一個不同于前兩個結(jié)果的抑制性刺激增強電活動的新例證,豐富了神經(jīng)動力學(xué)關(guān)于反常電活動的內(nèi)涵.

2)抑制性作用增強電活動的反常現(xiàn)象可以用PRC進行解釋,揭示了抑制性電流脈沖激勵引起峰放電增強的動力學(xué)機制.以前的研究主要關(guān)注了正向方波脈沖電流可使靠近Hopf分岔點的放電提前或延遲,即II型PRC,并用PRC揭示了多類神經(jīng)電活動的動力學(xué)特征如興奮性和同步等[7?9,38?40],這些結(jié)果有利于認識興奮性作用的特性和功能.而本文研究了負向方波脈沖電流作用于放電的不同相位的PRC,發(fā)現(xiàn)該曲線與正向方波脈沖引起的PRC的形狀關(guān)于刺激時刻的坐標(biāo)軸大約是對稱的,其相位也是有正有負的,即放電提前或延遲,從而揭示了抑制性電流脈沖激勵引起峰放電增強的動力學(xué)機制.因此,該II型PRC將有助于進一步理解-抑制性作用下神經(jīng)元電活動的動力學(xué)性質(zhì)和功能.

3)本文還揭示了抑制性自反饋可以增強放電精確性,這是在較長的時滯下產(chǎn)生的,與實驗發(fā)現(xiàn)的慢抑制性自突觸誘發(fā)精確放電相一致[15]; 進一步,抑制性耦合越強,精確放電越高,也與實驗結(jié)果相一致[15].研究結(jié)果理論仿真了抑制性自突觸誘發(fā)和增強放電精確性的實驗結(jié)果,這就通過理論模型揭示了抑制性自突觸新的潛在功能.而放電精確性是神經(jīng)系統(tǒng)實現(xiàn)生理功能的重要保證.

因此,本研究揭示了新的神經(jīng)電活動的反?，F(xiàn)象和相應(yīng)的動力學(xué)機制,提供了調(diào)控神經(jīng)電活動和生理功能的新手段,對于進一步認識抑制性作用的動力學(xué)機制和生理功能有重要的意義.