Pre-Bötzinger复合体是哺乳动物产生呼吸节律的脑干区域。磁通控制忆阻器用于连接膜电位和磁通量。在Butera动力学模型的基础上引入忆阻器,研究电流和磁流对神经元放电模式的影响及其产生的动力学机理。结果表明,电流和磁流作用都可以使pre-Bötzinger复合体放电模式由方波簇转迁到混合簇,再转迁为调幅峰。 The pre-Bötzinger complex in the brainstem region can generate respiratory rhythm in mammals. The magnetic flux-controlled memristor is used to bridge the membrane potential and the magnetic flux. By adding the magnetic flow to Butera’s neuron model, we study how the electric current and magnetic flux influence neuronal activities and its dynamic mechanism. The results show that both electric current and magnetic flux can transform the firing patterns of the pre-Bötzinger complex from square wave bursting to mixed bursting, and then amplitude modulation spiking.
Pre-Bötzinger复合体是哺乳动物产生呼吸节律的脑干区域。磁通控制忆阻器用于连接膜电位和磁通量。在Butera动力学模型的基础上引入忆阻器,研究电流和磁流对神经元放电模式的影响及其产生的动力学机理。结果表明,电流和磁流作用都可以使pre-Bötzinger复合体放电模式由方波簇转迁到混合簇,再转迁为调幅峰。
Pre-Bötzinger复合体,忆阻器,磁流,分岔
Wenchao Ji*, Moutian Liu, Huiru Qi
School of Science, North China University of Technology, Beijing
Received: Apr. 1st, 2021; accepted: Apr. 20th, 2021; published: Apr. 27th, 2021
The pre-Bötzinger complex in the brainstem region can generate respiratory rhythm in mammals. The magnetic flux-controlled memristor is used to bridge the membrane potential and the magnetic flux. By adding the magnetic flow to Butera’s neuron model, we study how the electric current and magnetic flux influence neuronal activities and its dynamic mechanism. The results show that both electric current and magnetic flux can transform the firing patterns of the pre-Bötzinger complex from square wave bursting to mixed bursting, and then amplitude modulation spiking.
Keywords:Pre-Bötzinger Complex, Memristor, Magnetic Flux, Bifurcation
Copyright © 2021 by author(s) and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
神经科学是一门综合性学科,可以对动物和人类的神经系统进行生物学研究。随着神经实验和计算技术的发展与创新,神经科学研究取得了长足的进步。为了全面而深入地了解大脑活动的规律性,人们已经深刻认识到建立神经科学理论体系的重要性。Pre-Bötzinger复合体(pre-BötC)是兴奋性耦合的神经元网络。实验发现pre-BötC对呼吸节律的产生至关重要 [
非线性动力学是研究非线性动力系统中各种运动状态的定量和定性规律,特别是运动模式演化行为的科学。一般来说,非线性科学的主体包括混沌、分岔、孤立子和斑图等。其中,分岔和混沌是两个最基本的非线性现象,由于它们的主要特点是在动力学的基础上去揭示事物非线性本质特性,所以它们成为了非线性科学研究领域中的活跃主题。Rinzel [
电活动是普遍存在的神经元生物电现象,它有许多不同的模式来编码生物信息的表达,构成神经元之间信号传播的整个过程。近年来,电流和磁流对神经元放电活动的影响引起科学家的关注 [
为了研究电流和磁流对pre-BötC神经元活动的影响,我们在Butera [
d V d t = ( − I L − I K − I Na − I NaP − I tonice − I CAN + I − k 1 V ρ ( φ ) ) / C (1)
d n d t = ( n ∞ ( V ) − n ) / τ n ( V ) (2)
d h d t = ( h ∞ ( V ) − h ) / τ h ( V ) (3)
d φ d t = V − k 2 φ (4)
钙动力学为:
d [ Ca ] d t = f m ( J ER IN − J ER OUT ) (5)
d l d t = A K d ( 1 − l ) − A [ Ca ] l (6)
在方程(1)中,V是膜电位,C是膜电容,t是时间。IL, IK, INa, INaP, Itonice和ICAN分别代表漏电流,钾电流,快速钠电流,持续钠电流,兴奋性刺激产生的电流和钙激活的非特异性阳离子电流。I是外部刺激的直流电。 k 1 V ρ ( φ ) 可用于描述电磁感应引起的感应电流:参数k1对膜电位上的电流起到调节作用, ρ ( φ ) 是指磁通控制忆阻器的记忆电导,用来描述磁通量 φ 与V之间的关系,其表达式为:
ρ ( φ ) = d q ( φ ) d t = α + 3 β φ 2 ,其中 α , β 是固定的参数值。定义 I a p p = I − k 1 V ρ ( φ ) 来表示由电流和磁流引
起的外激励总和。方程(1)右边各离子电流描述如下:
I L = g L ( V − V L ) , I K = g K n 4 ( V − V K ) ,
I Na = g Na m ∞ 3 ( V ) ( 1 − n ) ( V − V Na ) , I NaP = g NaP m p ∞ ( V ) h ( V − V Na ) ,
I tonice = g tonice ( V − V syn_e ) , I CAN = g CAN f ( [ Ca ] ) ( V − V Na ) ,
f ( [ Ca ] ) = 1 / ( 1 + ( K CAN / [ Ca ] ) n CAN ) 。
方程(2)和(3)中的n和h是电压门控钾和钠通道的门控变量,且这两个方程以及离子电流表达式中出现的平衡函数 n ∞ ( V ) , h ∞ ( V ) , m ∞ ( V ) , m p ∞ ( V ) 可表示为: x ∞ ( V ) = 1 / ( 1 + exp ( V − θ x / σ x ) ) , x ∈ { n , h , m , m p } 。方程(2)和(3)中所涉及的时间尺度函数表达式为: τ x ( V ) = τ ¯ x / cosh [ ( V − θ x ) / 2 σ x ] , x ∈ { n , h } 。
方程(4)中的 φ 表示跨膜的磁通量,可以描述电磁感应和磁场的影响。k2用来抑制磁通量的增加。根据法拉第电磁感应定律和忆阻器的描述,电荷对磁通量的依赖性由记忆电导定义如下:
i ′ = d q ( φ ) d t = d q ( φ ) d φ d φ d t = k 1 V ρ ( φ ) 。方程(1)-(4)构成胞体子系统。
方程(5)和(6)组成细胞内的钙动力系统,可以认为是细胞内的树突子系统。方程(6)中的l表示ER膜中尚未灭活的IP3通道比例,这取决于细胞内钙浓度([Ca])。方程(5)指定[Ca]是由从胞浆( J ER IN )进入内质网的通量和从内质网进入胞浆( J ER OUT )的通量决定的。这些通量受细胞内IP3的浓度[IP3]以及IP3通道门控变量l的调节,并表示为
J ER IN = ( L IP 3 + P IP 3 [ [ IP 3 ] [ Ca ] l ( [ IP 3 ] + K I ) ( [ Ca ] + K a ) ] 3 ) ⋅ ( [ Ca ] ER − [ Ca ] ) ,
J ER OUT = V SERCA [ Ca ] 2 K SERCA 2 + [ Ca ] 2 。
假设细胞内的钙总量([Ca]Tot)是恒定的。在此假设下,[Ca]ER可按如下方式计算:
[ Ca ] ER = [ Ca ] Tot − [ Ca ] σ 。
方程(1)~(6)组成全系统。模型中的参数取值见表1。
参数 | 参数值 | 参数 | 参数值 | 参数 | 参数值 | 参数 | 参数值 |
---|---|---|---|---|---|---|---|
g L | 2.3 nS | V L | −65 mV | σ n | −4 mV | L IP 3 | 0.37 pL∙s−1 |
g K | 11.2 nS | V K | −85 mV | σ h | 5 mV | P IP 3 | 31000 pL∙s−1 |
g Na | 28 nS | V Na | 50 mV | σ m | −5 mV | [ Ca ] Tot | 1.25 μM |
g NaP | 2 nS | V syn_e | 0 mV | σ m p | −6 mV | f m | 0.000025 pL−1 |
g tonice | 0.3 nS | θ n | −29 mV | K I | 1.0 μM | V SERCA | 400 aMol∙s−1 |
g CAN | 0.7 nS | θ h | −48 mV | K a | 0.4 μM | K SERCA | 0.2 μM |
τ ¯ n | 10 ms | θ m | −34 mV | K CAN | 0.74 μM | A | 0.005 μM−1∙s−1 |
τ ¯ h | 10000 ms | θ m p | −40 mV | n CAN | 0.97 | K d | 0.4 μM |
C | 21 μF | σ | 0.185 | [ IP 3 ] | 0.96 μM |
表1. 模型中的参数值
实验证实了持续钠电流(INaP)和钙激活的非特异性阳离子电流(ICAN)在pre-BötC神经元中普遍存在,并在神经元节律脉冲中起着重要的作用 [
除了pre-BötC神经元的呼吸网络外,树突状钙振荡与体细胞之间的相互作用可能发生在大脑区域的各种神经元中,包括丘脑皮质神经元、新皮质神经元和海马体中间神经元以及小脑颗粒细胞。因此,本文的研究可应用于其他类神经元细胞模型中。
当外激励 I app = 0 时,系统的簇放电模式如图1(a)所示。黑色实线是膜电位V随时间t的变化过程,红色虚线是细胞内钙浓度[Ca]随时间t的变化过程。此时系统呈现方波簇。
在系统中,V,n和 φ 是快变量,h和[Ca]是慢变量,l是超慢变量。选择慢变量[Ca]作为快子系统的分岔参数,并对h取平均值(h = 0.4192)进行分岔分析。为方便起见,用 g CANTot = g CAN f ( [ Ca ] ) 代替[Ca]作为分岔参数,其中 f ( [ Ca ] ) 是单调递增函数。快子系统对慢变参数 g CANTot 的分岔如图1(b)所示。平衡点形成S形曲线(黑色实线),曲线的下分支和中支分别由稳定结点和不稳定鞍点组成。曲线的上分支由稳定和不稳定的焦点组成,不稳定焦点经由Hopf分岔(subH)变为稳定焦点,并在Hopf分岔处产生极限环(红色实线)。点F1和F2表示平衡点的鞍结分岔。全系统在( g CANTot ,V)平面上的轨线(绿色曲线)也叠加在分岔图上。当 g CANTot 增大时,S型曲线下支的静息态经由平衡点的鞍结分岔F1跃迁至上支的稳定极限环,并由于极限环的吸引反复振荡,振荡态最终跃迁至下支的静息态,从而完成了一个周期振荡。根据Izhikevich簇放电分类的标准 [
图1. I app = 0 时的簇放电模式和分岔分析。(a) 膜电位V的时间序列;(b) 分岔分析
本节研究电流对混合簇的影响。当 k 1 = 0.1 , k 2 = 3 s − 1 , α = 1 M Ω − 1 , β = 0.00006 M Ω − 1 ⋅ V − 2 ⋅ s − 2 时,不同电流值所对应的簇放电模式如图2所示。黑色实线是膜电位V随时间t的变化过程,红色虚线是细胞内钙浓度[Ca]随时间t的变化过程。
当 I = − 2 μA 时,外激励电流的值比较小,与无外界刺激时的放电状态类似,系统呈现方波簇放电,如图2(a)所示。当电流I的值分别为5 μA, 10 μA, 15 μA时,系统产生胞体-树突簇(somato-dendritic bursting),是一种特殊的混合簇,分别如图2(b)~(d)所示。此种胞体-树突簇与实验中观察到的一种叹息簇 [
图2. 不同电流值下的簇放电模式。(a) I = − 2 μ A ;(b) I = 5 μ A ;(c) I = 10 μ A ;(d) I = 15 μ A ;(e) I = 20 μ A
峰峰间期,即ISI序列能够展现从方波簇到混合簇再到调幅峰的复杂转迁规律。取对数后的ISI序列随电流I的变化情况如图3所示。随着电流的增大,两簇之间由原来静息的状态逐步变为尖峰放电状态,导致前一个簇结束的时间与后一个簇开始的时间之间的间距逐渐变小,直至消失,全部变为峰放电状态,对应最大ISI逐渐降低。
图3. ISI序列随电流I的变化情况
将不同电流值下的簇放电模式在( g CANTot ,V)平面上的轨迹叠加到快子系统的分岔图上,如图4(a)~(e)所示。平衡点形成S形曲线,其他字母及曲线表示的含义与图1(b)相同,且簇放电模式均为“fold/homoclinic”型簇放电。
图4. 不同电流值下的快子系统的分岔,插图是局部放大。(a) I = − 2 μ A ( h = 0.3692 ) ;(b) I = 5 μ A ( h = 0.2669 ) ;(c) I = 10 μ A ( h = 0.2004 ) ;(d) I = 15 μ A ( h = 0.1404 ) ;(e) I = 20 μ A ( h = 0.0919 )
下面我们固定 I = 0 ,设置参数 k 2 = 3 s − 1 , α = 1 M Ω − 1 , β = 0.00006 M Ω − 1 ⋅ V − 2 ⋅ s − 2 ,研究磁流反馈系数 k 1 对簇放电模式的影响。不同 k 1 值所对应的簇放电模式如图5所示。黑色实线是膜电位V随时间t的变化过程,红色虚线是细胞内钙浓度[Ca]随时间t的变化过程。其中,图5(a)是方波簇,图5(b)和图5(c)为胞体–树突簇,图5(d)是调幅峰。与电流作用效果相类似,增大磁流也可以实现由方波簇到混合簇,再到调幅峰的转换。
图5. 不同 k 1 值下的簇放电模式。(a) k 1 = 0.01 ;(b) k 1 = 0.2 ;(c) k 1 = 0.4 ;(d) k 1 = 0.6
取对数后的ISI序列随参数 k 1 的变化情况如图6所示。与ISI序列随电流I的变化情况类似,随着参数 k 1 的增大,两簇之间由原来静息的状态逐步变为尖峰放电状态,导致前一个簇结束的时间与后一个簇开始的时间之间的间距逐渐变小,直至消失,最后全部变为峰放电状态,对应最大ISI逐渐降低。大约在 k 1 = 0.5 时簇放电转变为峰放电。
图6. ISI序列随参数k1的变化情况
神经动力学是研究神经系统的动力学与控制问题的重要理论基础,涉及脑神经电生理、信息和认知活动的动力学建模和行为分析等,是国际前沿研究领域。本文主要应用ISI序列以及快慢变量分岔分析的方法,研究电磁流对神经元放电模式的影响。结果表明,电流和磁流都可以实现由方波簇到调幅峰的转迁。在转迁过程中,随着外激励强度的不断增加,两簇之间的静息状态逐步变为峰放电状态,导致一种特殊的混合簇:胞体–树突簇的出现。外激励强度继续增加会导致原来两簇之间静息态完全消失,全部转变为峰放电状态。由于[Ca]的波动,峰放电之间振幅有所不同,我们称此时的放电状态为调幅峰放电。这一系列过程即是由于外激励强度增加,导致系统放电状态由方波簇到混合簇再到调幅峰的转迁的过程。利用ISI序列展现由方波簇到调幅峰转迁的复杂过程,随着外激励增加,最大ISI序列的逐渐降低对应两簇之间静息态持续时间的逐渐减少。分岔分析图中全系统轨线与时间序列图中的簇放电状态能够很好的对应,进一步解释了这一转迁背后的动力学机理。
调幅峰现象较为新颖,之前的研究极少出现,本文实现了由方波簇到调幅峰的转换,并揭示其背后的动力学机理,为神经动力学行为添加了新的示例,进一步丰富了非线性动力学行为的内容。以往的研究多关注离子电导的影响及两个细胞的耦合、同步等现象,而本研究则关注电磁作用的影响。研究结果有助于了解电磁流作用的动力学机理,为神经系统疾病的预防和治疗提供一定的理论支持。
本文得到国家自然科学基金(项目编号:11872003)及北方工业大学“特色学科建设”项目(项目编号:110052972027/014)的资助。
冀文超,刘谋天,齐会如. 电磁作用实现由方波簇到调幅峰的转迁Realizing the Transition of Electromagnetic Action from Square Wave Bursting to Amplitude-Modulation Spiking[J]. 动力系统与控制, 2021, 10(02): 129-138. https://doi.org/10.12677/DSC.2021.102014