電場對中間神經(jīng)元亞線性樹突整合的調(diào)制作用

魏熙樂，樊亞琴，趙美佳，伊國勝，王?江

電場對中間神經(jīng)元亞線性樹突整合的調(diào)制作用

魏熙樂，樊亞琴，趙美佳，伊國勝，王?江

(天津大學電氣自動化與信息工程學院，天津 300072)

中間神經(jīng)元的樹突形態(tài)大多呈現(xiàn)出空間對稱性分布，使得電場在其胞體附近引起的膜電位極化作用微弱，因此大多數(shù)研究認為電場對中間神經(jīng)元難以產(chǎn)生調(diào)制作用. 然而，電場在中間神經(jīng)元末端的樹突棘位置能夠誘發(fā)較強的膜電壓極化效應，遠端的樹突極化是否以及如何影響中間神經(jīng)元的樹突整合特性尚不十分清楚. 為此本文針對分散和集中空間分布的AMPA類型突觸輸入，研究了電場對中間神經(jīng)元全局和局部亞線性樹突整合特性的調(diào)制作用. 為了描述電場作用下樹突對AMPA突觸輸入積分的潛在規(guī)律，首先建立了描述中間神經(jīng)元輸入-輸出關系的簡化電路模型，此模型具有多個被動樹突分支且以不均勻細胞外膜電勢表征電場誘發(fā)的極化效應；其次，基于奇異攝動理論分析了電場調(diào)節(jié)下的亞線性樹突整合的動力學特性，推導出了可描述局部和全局樹突整合的閾下輸入-輸出關系的漸近解析表達式. 理論分析結(jié)果表明，電場通過改變AMPA突觸輸入的驅(qū)動力調(diào)節(jié)中間神經(jīng)元局部以及全局亞線性樹突整合特性，調(diào)制效應依賴于電場對遠端樹突的極化程度. 相比于局部樹突整合，電場對全局樹突整合調(diào)制更敏感. 最后，利用電場作用下具有真實樹突形態(tài)的海馬CA1區(qū)中間神經(jīng)元模型仿真驗證了理論分析結(jié)果.

中間神經(jīng)元；電場；極化；亞線性樹突整合

在中樞神經(jīng)系統(tǒng)中，興奮性和抑制性神經(jīng)元的共同激活被認為是各種大腦活動的基本原理[1]．抑制性神經(jīng)元主要為GABA能中間神經(jīng)元，GABA能中間神經(jīng)元通過釋放GABA神經(jīng)遞質(zhì)使得細胞膜超極化而抑制突觸后神經(jīng)元，進而控制和調(diào)節(jié)興奮性神經(jīng)元的活動[2]，因此促進了信息的平穩(wěn)流動以及維護了神經(jīng)回路的穩(wěn)態(tài)．興奮性與抑制性(excitatory/inhibit-tory，E/I)神經(jīng)元比例失衡導致癲癇、精神分裂癥、阿爾茨海默病(Alzheimer’s disease，AD)和帕金森病等神經(jīng)精神疾病[3-6]．大量研究表明，GABA能中間神經(jīng)元已經(jīng)成為了治療神經(jīng)精神疾病的有效靶點．例如，Guzman等[7]對近年來GABA能信號系統(tǒng)作為治療AD疾病的有效藥物靶點的研究進行了概述；Xu等[8]發(fā)現(xiàn)調(diào)節(jié)GABA能系統(tǒng)是一種用于治療精神分裂癥的有效干預方法，并且討論了將GABA能神經(jīng)元作為認知缺陷的潛在治療靶點和藥物治療的可行性．因此研究中間神經(jīng)元有助于為治療神經(jīng)精神疾病提供理論指導．

經(jīng)顱磁刺激[9-10](transcranial magnetic stimula-tion，TMS)是一種常用的有效的神經(jīng)調(diào)控技術(shù)，利用短暫的電流通過電磁線圈所產(chǎn)生的感應電場調(diào)節(jié)神經(jīng)系統(tǒng)，進而治療腦疾病[11-12]．錐體神經(jīng)元不對稱的樹突形態(tài)使得電場在胞體處產(chǎn)生顯著的極化作用，因此之前研究重點關注電場對錐體神經(jīng)元的調(diào)制作用．Bikson等[13]在大鼠海馬切片實驗中，發(fā)現(xiàn)電場通過在近胞體處的極化調(diào)節(jié)了神經(jīng)元的興奮性，改變了神經(jīng)元放電閾值和放電起始位置；Radman等[14]發(fā)現(xiàn)電場可以改變海馬錐體神經(jīng)元的放電時刻．電場在單神經(jīng)元水平引起的放電時刻、放電閾值的變化會經(jīng)過網(wǎng)絡作用得到放大[15]，從而調(diào)節(jié)海馬區(qū)節(jié)律．對于樹突形態(tài)對稱的中間神經(jīng)元而言，電場在胞體誘發(fā)的極化十分微弱[16]，其胞體極化的重要性基本被忽略[17]，因此大多數(shù)研究認為電場對中間神經(jīng)元基本沒有調(diào)制作用．Murphy等[18]發(fā)現(xiàn)，TMS產(chǎn)生的感應電場通過激活GABA能中間神經(jīng)元而抑制皮層深層錐體神經(jīng)元的樹突活動，說明電場可以調(diào)節(jié)GABA能中間神經(jīng)元的放電活動，進而影響神經(jīng)回路的E/I平衡．但是在突觸背景下，電場如何通過調(diào)制中間神經(jīng)元樹突整合而影響胞體響應尚不清楚．

樹突整合被定義為突觸后神經(jīng)元將突觸前信息轉(zhuǎn)換為興奮性突觸后電勢(excitatory postsynaptic potential，EPSP)，并傳輸至胞體的過程，體現(xiàn)了對突觸輸入的整合，樹突濾波以及胞體響應特性[19]. NMDA和AMPA作為兩類重要的興奮性突觸，廣泛分布在神經(jīng)元的樹突上[20]．NMDA受體激活使得樹突進行超線性整合[19]，而AMPA受體激活使得樹突進行亞線性整合[21]且有助于NMDA受體激活．非線性整合特性極大豐富了神經(jīng)元計算功能，例如特征識別[22]．突觸電流動態(tài)特性在塑造樹突亞線性整合中起著關鍵作用，對于相同的突觸電導其驅(qū)動力為膜電勢與突觸反電勢的差值，驅(qū)動力飽和導致樹突亞線性計算[23]．盡管加電場對中間神經(jīng)元胞體的極化作用微弱，其在遠端樹突產(chǎn)生的較強極化作用可以改變突觸電流的驅(qū)動力，因此本文提出假設：電場能夠通過改變中間神經(jīng)元遠端樹突膜電位，進而調(diào)節(jié)AMPA亞線性樹突整合特性．

本文建立了一個能夠描述多樹突分支、多突觸輸入的中間神經(jīng)元多間室電路模型，各個間室具有不同的胞外膜電勢從而刻畫了電場的空間極化效應．基于奇異攝動理論建立中間神經(jīng)元動力學與AMPA類型樹突整合之間的聯(lián)系，推導出其輸入-輸出關系的漸近解析表達式．理論分析預測出電場的強度與極性在遠端樹突產(chǎn)生的不同極化對亞線性樹突整合特性的調(diào)制作用，并對比了全局與局部突觸分布對電場調(diào)制敏感性的影響．最后，利用電場作用下樹突形態(tài)呈放射狀對稱的海馬CA1區(qū)中間神經(jīng)元模型進行仿真，驗證理論分析結(jié)果．

1?模型與方法

1.1?中間神經(jīng)元簡化電路模型

1.1.1?單樹突分支模型

不同于錐體神經(jīng)元，大多數(shù)中間神經(jīng)元的樹突形態(tài)呈現(xiàn)空間對稱性[17]，因此不區(qū)分頂樹突和基樹突．本文以樹突形態(tài)呈星狀放射狀對稱的海馬CA1區(qū)中間神經(jīng)元模型[24]作為研究對象，對此進行建模，其三維樹突形態(tài)如圖1所示．由于中間神經(jīng)元樹突分支空間分布對稱且大都直接連接至胞體，因此構(gòu)建的模型樹突形態(tài)對稱，且各個樹突分支的電纜模型結(jié)構(gòu)相同，如圖1(a)所示．該模型由胞體、遠端樹突和近端樹突3部分構(gòu)成，刻畫了真實神經(jīng)元遠端樹突接收突觸信號并經(jīng)近端樹突傳遞至胞體的過程．

圖1?單樹突分支電纜及三間室電路模型

表1?模型被動參數(shù)取值

Tab.1?Passive parameter values of cable model

由單樹突分支電纜模型得到其三間室模型，如圖1(b)所示，該模型由遠端樹突刺激間室、近端樹突間室和胞體間室構(gòu)成，每個間室包含膜電容與漏電導，膜電容取值為1μF/cm2．兩個相鄰的間室通過軸向電導連接，根據(jù)電導計算公式，得到每個間室的漏電導以及連接相鄰兩個間室的軸向電導表達式為

神經(jīng)元的每個間室都可以用一個具有獨立膜電壓的等效電路來表示，設置每個間室胞外膜不同的電勢來體現(xiàn)電場作用下神經(jīng)元不均勻的空間極化效應．圖1(b)中，胞內(nèi)與胞外膜電勢的差值為跨膜電壓，表達式為

式中md1、md2、ms分別為3個間室的跨膜電壓，由圖1(b)列出每個間室跨膜電壓的微分方程為

定義新的狀態(tài)變量為

利用坐標變換得到單樹突分支的等效3間室模型(見圖1(c))．將方程(3)與(4)做差值，得到3個間室膜電位變化量的動力學方程為

利用坐標變換，得到電場作用下可刻畫突觸輸入變化所引起的細胞膜電壓變化量的狀態(tài)方程，可用圖1(c)的單樹突分支等效3間室模型描述．

1.1.2?中間神經(jīng)元的多樹突分支模型

模型中多個樹突分支可以看作并聯(lián)關系，再通過軸向電導in_s與胞體間室連接，因此從個刺激分支流出的電流一部分流入胞體間室，一部分流入到非刺激樹突分支，可表示為

由此建立了一個描述多樹突分支、多突觸輸入非對稱分布的中間神經(jīng)元電路模型，該電路模型描述了突觸輸入電流沿遠端樹突-近端樹突-胞體路徑傳輸?shù)桨w的過程．

1.2?具有真實樹突形態(tài)的中間神經(jīng)元電場模型

NEURON仿真平臺[25]是一個可以替代真實神經(jīng)元進行虛擬實驗的工具，其網(wǎng)址為https：//neuron.yale.edu/neuron/．本文在NEURON軟件中的ModelDB數(shù)據(jù)庫中選取具有真實樹突形態(tài)的海馬CA1區(qū)中間神經(jīng)元模型進行仿真實驗(見圖3)，該模型由Saraga等[26]對一只小鼠的海馬CA1區(qū)細胞進行重建得到，由胞體和呈放射狀對稱的三維樹突形態(tài)構(gòu)成，此模型下載網(wǎng)址為：https：//senselab.med. yale.edu/ModelDB/ShowModel？model＝28316．該模型的空間形態(tài)可以通過NEURON中的“點過程”機制實現(xiàn)，其樹突直徑變化范圍從近端樹突1.38μm到遠端樹突0.61μm；樹突長度范圍為70～330μm；軸向電阻為120Ω·cm；膜電容為1μF/cm2；胞體靜息膜電位－65.21mV．由于閾下電場刺激下樹突的主動離子通道活動較弱，本文中樹突僅考慮了被動特性，重點研究電場對遠端突觸輸入整合的調(diào)節(jié)作用，因此樹突中不含主動離子通道．胞體中含有鈉離子通道和鉀離子通道，鈉離子和鉀離子電導分別為0.01S/cm2、0.0319S/cm2，AMPA電導為0.3nS．

基于NEURON軟件中的“胞外機制”建立勻強電場，該建模方式常用在NEURON中模擬在電場刺激下神經(jīng)元的空間極化效應以及其電活動[27]．電場對神經(jīng)元興奮性的調(diào)制作用是通過其引起的空間極化效應實現(xiàn)的，如圖3所示，靠近電場陽極的樹突其膜電位降低而發(fā)生超極化，靠近電場陰極的樹突其膜電位升高而發(fā)生去極化．本文對電場正負極性的定義如下：電場強度＞0V/m時為正極性電場(見圖3(a))，此時胞體去極化；電場強度＜0V/m時為負極性電場，此時胞體超極化．圖3(b)中，當電場強度從－20V/m以5V/m的步長增加至20V/m時，胞體最大極化大約為0.5mV，而錐體神經(jīng)元胞體極化為5～10mV[16]，表明在電場作用下中間神經(jīng)元的胞體極化效應微弱，遠小于錐體神經(jīng)元的胞體極化作用．然而電場在中間神經(jīng)元遠端樹突的膜電位極化作用較大，并且遠端樹突膜電位隨電場強度變化而線性增加．

圖3?海馬CA1區(qū)中間神經(jīng)元模型以及電場空間極化效應

1.3?亞線性樹突計算分析方法

樹突整合特性常用閾下輸入-輸出關系描述，即可以通過比較期望EPSP與觀測EPSP之間的關系描述，其中多個突觸輸入同時激活所產(chǎn)生的胞體EPSP稱為實際觀測EPSP，而相同數(shù)目下依次激活的單個突觸輸入得到的胞體EPSP的代數(shù)和稱為期望EPSP．當實際觀測EPSP大于期望EPSP時，樹突整合表現(xiàn)為超線性計算；當實際觀測EPSP等于期望EPSP時，樹突整合為線性計算；當實際觀測EPSP小于期望EPSP時，樹突整合表現(xiàn)為亞線性計算．

在中間神經(jīng)元的一個樹突分支施加AMPA突觸刺激，突觸位置如圖4(a)所示，每個突觸刺激位置用數(shù)字表示．實測EPSP為突觸同時激活時胞體實際膜電壓與靜息膜電壓的差值，期望EPSP為依次激活的單個突觸激活產(chǎn)生的實測EPSP的代數(shù)和．突觸刺激個數(shù)從1增大到6時得到實測和期望的胞體膜電壓波形，分別如圖4(a)和圖4(b)所示．隨著突觸個數(shù)的增加，胞體膜電壓增大的幅度逐漸減小，這是由于在樹突為被動時對AMPA類型輸入執(zhí)行亞線性樹突計算，如圖4(c)所示．

圖4?亞線性樹突計算描述方法

2?結(jié)?果

2.1?亞線性樹突整合動力學分析

奇異攝動系統(tǒng)是帶有小參數(shù)的動態(tài)系統(tǒng)，具有多時間尺度的特點，其本質(zhì)上是時間尺度分析法的一種近似解法．基于圖1(c)單樹突分支3間室等效電路模型，建立奇異攝動系統(tǒng)分析亞線性樹突整合的動力學特性．為了探究中間神經(jīng)元輸入-輸出關系之間的動力學機制，從而對由膜電位組成的狀態(tài)向量引入一個擾動參數(shù)，將AMPA()看作奇異攝動系統(tǒng)的慢變量，Δms()看作快變量，由此構(gòu)建奇異攝動系統(tǒng)為

式中：為神經(jīng)元膜電容；為攝動參數(shù)．依據(jù)變量時間尺度的不同，可由奇異攝動系統(tǒng)得到快子系統(tǒng)和慢子系統(tǒng)，通過把快、慢變量分開可把系統(tǒng)轉(zhuǎn)化成低階模型以方便分析．

AMPA突觸電流表達式為

式中：syn為突觸反電勢，取值為0mV；AMPA(0)為AMPA突觸電導初始狀態(tài)，取值為0mS．AMPA突觸電導隨時間變化的表達式為

將式(14)結(jié)合狀態(tài)方程(6)代入到奇異攝動系統(tǒng)方程(9)中，得到快子系統(tǒng)的狀態(tài)方程為

根據(jù)快子系統(tǒng)的狀態(tài)方程刻畫其平衡點曲線，如圖5(a)所示．在樹突被動時，樹突對AMPA突觸輸入進行亞線性積分，系統(tǒng)均為穩(wěn)定的平衡點．設定AMPA突觸電導的最小值為min，最大值為max，即AMPA()∈[min，max]．快子系統(tǒng)中電導為常數(shù)，圖5(a)中豎直的虛線為快子系統(tǒng)的一維動態(tài)過程，隨著時間的變化膜電壓在虛線上振蕩，當快子系統(tǒng)時間尺度非常小時其收斂于穩(wěn)定的平衡點．在慢子系統(tǒng)中突觸電導是時變的并且隨著時間增加而先增大再減小，當突觸電導從min開始增大時，慢子系統(tǒng)動態(tài)跟隨穩(wěn)定的平衡點從點運動到點，此時電導和胞體膜電壓達到最大值，而后隨著時間的增大電導減小，慢流形沿著漸近平衡點曲線從點運動到點．圖5(b)為亞線性樹突整合時胞體膜電壓的動態(tài)過程，其與圖5(a)中的慢子系統(tǒng)動態(tài)相對應．

圖5?亞線性樹突整合動力學分析

樹突整合動力學特性描述了神經(jīng)元膜電位的動態(tài)曲線，并且突觸電導也是時變的，在變化的過程中突觸電導與膜電壓的最大值為動力學特性的體現(xiàn)．而神經(jīng)計算不考慮樹突整合動力學特性，在描述閾下輸入-輸出關系時，刻畫的突觸電導與膜電壓均為靜態(tài)的，并且對應了樹突整合動力學特性中的突觸電導與膜電壓的最大值．因此基于圖1(c)單樹突分支3間室等效電路模型，通過建立樹突整合動力學與樹突亞線性整合的聯(lián)系來證明可以用漸近平衡點曲線描述中間神經(jīng)元的閾下輸入-輸出關系．

將式(12)代入到系統(tǒng)狀態(tài)方程(6)中，得到系統(tǒng)的動態(tài)微分方程為

令單個AMPA突觸電導＝10-6mS，遠端樹突靜息膜電位為－64mV，采用Matlab求解式(16)，得到胞體膜電壓變化量ΔVms(t)的動態(tài)曲線．如圖6所示，令膜電容C＝1μF/cm2，當ε＝0.2時，由于系統(tǒng)時間常數(shù)較大而造成延遲，因此AMPA突觸電導gAMPA(t)與ΔVms(t)的曲線隨時間t變化呈環(huán)狀變化，gAMPA(t)到達最大值時ΔVms(t)并沒有到達峰值．而在快子系統(tǒng)中慢變量gAMPA(t)為常數(shù)，取時變下突觸電導gAMPA(t)的最大值gmax與ΔVms(t)的峰值ΔVmax作為特征值刻畫快子系統(tǒng)的動態(tài)軌跡，描述攝動參數(shù)ε對神經(jīng)元計算特性的影響(見圖6(a)、(b))．當ε逐漸減小，ΔVms(t)能盡快到達峰值，快子系統(tǒng)越快并且越接近時的極限狀態(tài)(見圖6(c))．圖6(d)中ΔV為ε的不同取值得到的胞體膜電位的變化量與平衡點電位的差值，可見ε越小，動態(tài)軌跡越接近漸近平衡點曲線．因此可以用快子系?統(tǒng)平衡點描述中間神經(jīng)元的閾下亞線性輸入-輸出?關系．

2.2?亞線性輸入-輸出關系的漸近表達式

在快子系統(tǒng)中突觸電導為常數(shù)，因此針對圖2中的中間神經(jīng)元多樹突分支模型，式(7)中AMPA突觸電流可以表示為

由于所有樹突分支看作并聯(lián)關系，因此每個分支的近端樹突間室膜電壓的變化量相同，并且由圖1(c)可以得到近端樹突與胞體膜電位變化量的關系，表示為

由式(8)、式(18)得到近端與遠端樹突膜電位變化量的關系為

此處考慮勻質(zhì)情況，因此每個樹突分支上的AMPA突觸電導和遠端樹突靜息膜電位相同，即

因此得到每個分支的遠端樹突膜電位變化量與分布突觸刺激分支的電流保持不變，則有

由式(17)～式(21)得到胞體膜電位變化量與遠端樹突膜電位的變化量以及突觸電流輸入的表達式分別為

式中md、V11、V12、V2為常數(shù)，且其表達式分別為

由式(14)得到快子系統(tǒng)中AMPA突觸電流表達式為

將式(24)代入到式(22)中，得到AMPA類型的中間神經(jīng)元亞線性輸入-輸出關系漸近解析表達式為

設單個AMPA突觸的最大電導為AMPA_s，單樹突分支的突觸電導可表示為

則多樹突分支上的多個突觸輸入激活時，中間神經(jīng)元胞體實際峰值膜電位的變化量為

式中nmsyn為激活突觸的全部數(shù)量，并且樹突的總分支數(shù)，突觸輸入激活的樹突分支數(shù)和nmsyn滿足的條件和關系表達式為

根據(jù)式(24)得到單樹突分支上的單個突觸輸入被激活時，中間神經(jīng)元胞體EPSP為

單樹突分支上的個AMPA突觸輸入激活時，中間神經(jīng)元胞體實際觀測的EPSP表達式為

在單個樹突分支上激活多個突觸輸入時，中間神經(jīng)元胞體期望的EPSP為

在多個樹突分支上的多個突觸刺激下，中間神經(jīng)元胞體期望的EPSP為

需要指出的是，分散的突觸輸入的樹突末端在電場刺激下的極化作用與角度相關，因此本文得出的表達式可以描述分散分布的突觸輸入不對稱情況下的樹突整合．

2.3?樹突靜息電位對亞線性樹突整合特性的影響

電場的空間極化效應使得中間神經(jīng)元遠端樹突的靜息膜電位發(fā)生改變，因此在簡化電路模型中通過刻畫中間神經(jīng)元遠端樹突靜息膜電位的變化來體現(xiàn)電場極化作用．根據(jù)上文得到的中間神經(jīng)元的閾下亞線性輸入-輸出關系的漸近表達式，采用Matlab仿真研究電場極化對中間神經(jīng)元樹突亞線性整合特性的影響．模型中考慮全局和局部兩種突觸分布形式，如圖7(a)與(b)所示，設置模型總分支數(shù)＝10，遠端樹突靜息膜電位為－64mV，AMPA的突觸電導為10-6mS，兩種分布下突觸個數(shù)nmsyn均從0變化到80個．局部式突觸分布為突觸輸入激活的遠端樹突的數(shù)量＝1；全局式突觸分布為＝7．

利用式(30)與式(31)得到遠端樹突靜息膜電位對局部亞線性樹突整合的影響；根據(jù)式(27)和式(32)得到遠端樹突靜息膜電位對非對稱全局亞線性樹突整合的影響，仿真結(jié)果如圖7所示．依據(jù)胞體期望EPSP與實際觀測EPSP的曲線關系，得出隨著遠端樹突靜息膜電位的增大，胞體飽和膜電位sat的變化規(guī)律．仿真結(jié)果表明AMPA突觸輸入下樹突的全局與局部整合均隨著突觸刺激的增強，胞體膜電位最終達到飽和．全局與局部樹突整合中發(fā)現(xiàn)相似的規(guī)律：中間神經(jīng)元遠端樹突在正向電場超極化時靜息膜電位降低，增強突觸電流的驅(qū)動力進而弱化樹突亞線性整合能力，使胞體輸出增大；而在負向電場中遠端樹突去極化，靜息膜電位升高，減小突觸驅(qū)動力進而增強樹突亞線性整合能力，使胞體輸出減?。虼苏{(diào)制作用取決于中間神經(jīng)元遠端樹突靜息膜電位的變化，也就是依賴于電場的極性與強度改變．

通過對比全局與局部突觸分布，研究電場極化對中間神經(jīng)元全局與局部亞線性樹突整合的敏感性，利用參數(shù)d衡量其調(diào)制的敏感性，d的計算方法為：對圖7(c)與(d)中的直線求解其斜率的絕對值，表示單位遠端樹突膜電位變化下EPSP的變化量．電場極化作用主要體現(xiàn)在遠端樹突膜電位，并且電場強度與遠端樹突膜電位的變化關系近似線性，因此d越大說明電場極化對其調(diào)制效果越明顯．仿真結(jié)果如圖7(e)所示，表明全局樹突整合得到的d遠遠大于局部樹突整合得到的d，這是由于電場對中間神經(jīng)元遠端樹突膜電位的極化呈同樣趨勢變化，即同時去極化同時超極化，使得電場對分散性突觸輸入的調(diào)制具有全局效應，對全局樹突整合的調(diào)節(jié)更敏感．

圖7?遠端樹突靜息電位對亞線性樹突整合特性的影響

2.4?真實樹突形態(tài)的中間神經(jīng)元NEURON仿真

為了驗證模型中電場極化對亞線性樹突整合的理論預測，利用在海馬區(qū)中間神經(jīng)元電場模型中加入AMPA突觸刺激進行分析．突觸分布形式如圖8(a)與(b)所示，局部式突觸分布為突觸隨機分布在神經(jīng)元的一個樹突分支上，而全局式分布是指，突觸隨機分布在一部分樹突上．在兩種分布下突觸個數(shù)從0增加到80個且在距離胞體大于60μm處同時激活，sat表示80個突觸輸入對應的實測EPSP(見圖8(a)).為了衡量電場對中間神經(jīng)元全局與局部亞線性樹突整合調(diào)制的敏感性．將圖8(c)中的點擬合成直線，Δ為＝20V/m與＝－20V/m對應的胞體EPSP峰值之差，Δ為電場變化量40V/m，樹突增益dend的計算公式為

式中dend表示單位電場變化下胞體EPSP的變化量，樹突增益越大說明電場對其調(diào)制效果越明顯．以此方法計算單位電場變化時所引起的胞體極化膜電壓的變化量，將圖3(b)中的曲線擬合成直線并計算該直線的斜率，記為soma，結(jié)果見圖8(d)中的藍色柱體．

圖8(a)和(b)中，當＝20V/m時，電場使刺激端樹突膜電位超極化而降低，增強了突觸電流驅(qū)動力進而使AMPA類型的亞線性樹突整合能力減弱，樹突輸入輸出關系更加線性化，故sat增大；當＝??-20V/m時，電場作用下刺激端樹突膜電位去極化而升高，降低突觸電流驅(qū)動力使AMPA亞線性整合能力增強，樹突輸入輸出關系亞線性化顯著，故sat減?。虼穗妶鰪姸葟模?0V/m以10V/m的步長增加至20V/m，sat隨著電場強度的變化而增加(見圖8(a))．觀察樹突增益發(fā)現(xiàn)，突觸為全局分布時得到的樹突增益遠大于局部分布的樹突增益，因此電場對中間神經(jīng)元全局亞線性樹突整合的調(diào)制更敏感，與理論簡化模型結(jié)論一致．將突觸分布不對稱的情況(見圖8(b))與對稱情況(見圖8(e))相比較，結(jié)果表明：當突觸空間分布呈對稱時，樹突增益有較明顯的下降(見圖8(f))，說明對稱性突觸輸入的整合對電場調(diào)節(jié)不敏感這一特性，與之前的研究表明電場對中間神經(jīng)元調(diào)節(jié)效果不顯著相一致．

將樹突增益與單位電場引起的胞體極化進行比較(見圖8(d))，發(fā)現(xiàn)單位電場下的樹突增益大約為胞體膜電位極化程度的2倍，因此電場通過極化遠端樹突膜電位引起的胞體膜電位變化量大于電場對胞體膜電位的極化，表明電場可以顯著調(diào)節(jié)中間神經(jīng)元的樹突整合從而影響其胞體響應，并且分散的突觸激活位置越不對稱效果越顯著．

目前研究表明同步激活的突觸信號更有利于神經(jīng)元信號整合[28]，事實上突觸并不是同時激活的．圖8(g)和(h)分析了突觸激活時間方差對電場調(diào)節(jié)亞線性樹突整合的影響，其中asy表示激活時間方差，其值分別為0ms、20ms和100ms．由此發(fā)現(xiàn)突觸激活時間方差相較突觸分布對電場的調(diào)節(jié)作用影響微弱，但激活方差越大使突觸輸入之間的交互減弱，樹突亞線性增強，因此神經(jīng)元更難放電．與文獻[28]結(jié)論?一致.

3?結(jié)?語

本文針對中間神經(jīng)元單分支的多間室模型建立了奇異攝動系統(tǒng)，對AMPA類型的局部或全局樹突亞線性整合進行了動力學分析，得到可描述中間神經(jīng)元多樹突分支的閾下輸入-輸出關系的漸近解析表達式．結(jié)果表明電場通過在中間神經(jīng)元遠端樹突膜電位的極化作用，影響了AMPA突觸電流驅(qū)動力，進而調(diào)節(jié)了中間神經(jīng)元的全局與局部亞線性樹突整合特性，調(diào)節(jié)作用依賴于電場的強度與極性，并且相比于局部樹突整合，電場對全局樹突整合的調(diào)制更敏感．利用真實樹突形態(tài)的中間神經(jīng)元模型進行NEURON仿真，從而驗證模型結(jié)果．本文研究表明，即使中間神經(jīng)元胞體極化作用微弱，但在突觸激活位置不對稱的情況下，電場調(diào)節(jié)下的中間神經(jīng)元的樹突積分會顯著影響其胞體輸出，從而影響神經(jīng)元計算，此結(jié)果有助于理解電場調(diào)制中間神經(jīng)元輸入-輸出功能的潛在機制，為闡釋磁刺激技術(shù)調(diào)制腦功能的機理提供理論?基礎．

［1］ Isaacson J S，Scanziani M. How inhibition shapes cortical activity[J]. Neuron，2011，72(2)：231-243.

［2］ Shetty A K，Bates A. Potential of GABA-ergic cell therapy for schizophrenia，neuropathic pain，and Alzheimer’s and Parkinson’s diseases[J]. Brain Research，2016，1638：74-87.

［3］ Southwell D G，Nicholas C R，Basbaum A I，et al. Interneurons from embryonic development to cell-based therapy[J]. Science，2014，344(6180)：167-187.

［4］ Amatniek J C，Hauser W A，DelCastillo-Castaneda C，et al. Incidence and predictors of seizures in patients with Alzheimer’s disease[J]. Epilepsia，2006，47(5)：867-872.

［5］ Shetty A K，Turner D A. Fetal hippocampal grafts containing CA3 cells restore host hippocampal glutamate decarboxylase-positive interneuron numbers in a rat model of temporal lobe epilepsy[J]. Journal Neuroscience，2020，20(23)：8788-8801.

［6］ Tyson J A，Anderson S A. GABAergic interneuron transplants to study development and treat disease[J]. Trends Neurosciences，2014，37(3)：169-177.

［7］ Guzman B C F，Vinnakota C，Govindpani K，et al. The GABAergic system as a therapeutic target for Alzheimer’s disease[J]. Journal of Neurochemistry，2018，146(6)：649-669.

［8］ Xu M Y，Wong A H C. GABAergic inhibitory neurons as therapeutic targets for cognitive impairment in schizophrenia[J]. Acta Pharmacologica Sinica，2018，39(5)：733-753.

［9］ Hallett M，Di lorio R，Rossini P M，et al. Contribution of transcranial magnetic stimulation to assessment of brain connectivity and networks[J]. Clinical Neurophysiology，2017，128(11)：2125-2139.

［10］ Ladenbauer J，Ladenbauer J，Kuelzow N，et al. Promoting sleep oscillations and their functional coupling by transcranial stimulation enhances memory consolidation in mild cognitive impairment[J]. The Journal of Neuroscience，2017，37(30)：7111-7124.

［11］ Huang Y，Liu A A，Lafon B，et al. Measurements and models of electric fields in thehuman brain during transcranial electric stimulation[J]. ELife，2017，6：1-26.

［12］ 于?凱，王?江，鄧?斌，等. 磁刺激作用下感生電場對網(wǎng)絡同步的影響[J]. 天津大學學報(自然科學與工程技術(shù)版)，2013，46(8)：726-736.

Yu Kai，Wang Jiang，Deng Bin，et al. Effects of induced electric field on network synchrony under magnetic stimulation[J]. Journal of Tianjin University(Sci-ence and Technology)，2013，46(8)：726-736(in Chinese).

［13］ Bikson M，Inoue M，Akiyama H，et al. Effect of uniform extracellular DC electric fields on excitability in rat hippocampal slices[J]. Journal of Physiology-London，2004，557(1)：175-190.

［14］ Radman T，Su Y Z，An J H，et al. Spike timing amplifies the effect of electric fields on neurons：Implications for endogenous field effects[J]. Journal of Neuroscience，2007，27(11)：3030-3036.

［15］ Reato D，Rahman A，Bikson M，et al. Low-intensity electrical affects network dynamics by modulating population rate and spike timing[J]. Journal of Neuroscience，2010，30(45)：15067-15079.

［16］ Radman T，Ramos R L，Brumberg J C，et al. Role of cortical cell type and morphology in subthreshold and suprathreshold uniform electric field stimulation[J]. Brain Stimulation，2009，2(4)：215-228.

［17］ Jackson M P，Rahman A，Lafon B，et al. Animal models of transcranial direct current stimulation：Methods and mechanisms[J]. Clinical Neurophysiol-ogy，2016，127(11)：3425-3454.

［18］ Murphy S C，Palmer L M，Nyffeler T，et al. Transcranial magnetic stimulation(TMS) inhibits cortical dendrites[J]. Elife，2016，5：1-12.

［19］ Stuart G J，Spruston N. Dendritic integration：60 years of progress[J]. Nature Neuroscience，2015，18(12)：1713-1721.

［20］ Stafford J，Brownlow M L，Qualley A，et al. AMPA receptor translocation and phosphorylation are induced by transcranial direct current stimulation in rats[J]. Neu-robiology of Learning and Memory，2018，150：36-41.

［21］ Abrahamsson T，Cathala L，Matsui K，et al. Thin dendrites of cerebellar interneurons confer sublinear synaptic integration and a gradient of short-term plasticity[J]. Neuron，2012，73(6)：1159-1172.

［22］ Vervaeke K，Lorincz A，Nusser Z，et al. Gap junctions compensate for sublinear dendritic integration in an inhibitory network[J]. Science，2012，335(6076)：1624-1628.

［23］ Tran-Van-Minh A，Caze R D，Abrahamsson T，et al. Contribution of sublinear and supralinear dendritic integration to neuronal computations[J]. Frontiers in Cellular Neuroscience，2015，9：1-15.

［24］ Poleg-Polsky A. Effects of neural morphology and input distribution on synaptic processing by global and focal NMDA-spikes[J]. PLoS One，2015，10(10)：1-21.

［25］ Hines M L，Carnevale N T. The NEURON simulation environment[J]. Neural Computation，1997，9(6)：1179-1209.

［26］ Saraga F，Wu C P，Zhang L，et al. Active dendrites and spike propagation in multi-compartment models of oriens-lacunosum/moleculare hippocampal interneurons[J]. Journal of Physiology-London，2003，552(3)：673-689.

［27］ 伊國勝，王?江，魏熙樂，等. 無創(chuàng)式腦調(diào)制的神經(jīng)效應研究進展[J]. 科學通報，2016，61(8)：819-834.

Yi Guosheng，Wang Jiang，Wei Xile，et al. Developments of neural effects induced by noninvasive brain modulation[J]. Chinese Science Bulletin，2016，61(8)：819-834(in Chinese).

［28］ Li Xiumin. Signal integration on the dendrites of a pyramidal neuron model[J]. Cognitive Neurodynamics，2014，8(1)：81-85.

Modulation Effect of Electric Fields on Sublinear Dendritic Integration of the Interneuron

Wei Xile，F(xiàn)an Yaqin，Zhao Meijia，Yi Guosheng，Wang Jiang

(School of Electrical and Information Engineering，Tianjin University，Tianjin 300072，China)

Most studies have suggested that the modulatory effect of the electrical field (EF) on the interneuron can be ignored due to the EF-induced weak somatic polarization that results from the symmetric distribution of dendritic morphology. However，EF induces strong polarization at distal dendritic spines of the interneuron. Whether and how these polarizations at the distal dendrites partake in the dendritic integration of alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) synapses remain unclear. In this regard，a simplified circuit model with multiple passive dendritic branches was established to capture the dendritic input-output relationship of the interneuron receiving asymmetrically distributed AMPA synapses，in which the inhomogeneous extracellular potentials characterize EF-induced polarizations. The dynamics of sublinear dendritic integration of interneuron were analyzed based on the singular perturbation theory. Our theoretical results showed that the EFs regulated the local and global sublinear dendritic integration properties of the interneuron by altering the driving force of AMPA synaptic inputs，depending on the degree of polarization at the distal dendrites. Additionally，the EF-regulated effect on global dendritic integration has higher efficiency compared to local dendritic integration. The simulation of an interneuron model with realistic dendritic morphology at the hippocampal CA1 in the case of EF stimulation verified the results of theoretical analysis.

interneuron；electric field；polarization；sublinear dendritic integration

10.11784/tdxbz202012009

Q64

A

0493-2137(2022)02-0199-12

2020-12-04；

2021-04-06.

魏熙樂（1975—??），男，博士，教授，xilewei@tju.edu.cn．Email：m_bigm@tju.edu.cn

伊國勝，guoshengyi@tju.edu.cn.

國家自然科學基金資助項目(61771330，62071324)；天津市自然科學基金資助項目(19JCQNJC01200).

the National Natural Science Foundation of China(No.61771330，No.62071324)，the Natural Science Foundation of Tianjin，China(No. 19JCQNJC01200).

(責任編輯：孫立華)