作者：秦杨，刘一辉，杨晶，王文挺，邢俊玲，胡三觉
【摘要】 目的：利用非线性动力学 理论 的 分析 手段阐明三叉神经中脑核神经元放电滞后现象的产生机制. 方法 ： 应用 脑片全细胞膜片钳技术记录大鼠三叉神经中脑核（mes v）神经元在三角波电流刺激下的放电特征，运用动力学方法对其进行分析. 结果：ⅱ型mesv神经元放电表现出滞后现象，加入4ap转型后滞后现象消失；非线性动力学分岔理论可以对以上现象做出解释. 结论：ⅱ型神经元的亚临界霍普夫分岔使得其存在滞后现象.
【关键词】 三叉神经中脑核；神经元；电生 理学 ；分岔；转型；i4ap
【abstract】 aim: to elucidate the discharge hysteresis phenomenon in mesencephalic v (mes v) neurons using nonlinear dynamics theory. methods: wholecell patch clamp technique was used for firing recording and nonlinear dynamic theory was used to analyze the producing mechanism of discharge hysteresis phenomenon in mes v neurons. results: hysteresis of discharge was observed in type 2 mes v neurons. after application of 4ap, type 2 firing pattern was transferred to type 1, in which hysteresis disappeared. the nonlinear dynamic bifurcation theory could provide explanation for this phenomenon. conclusion: it is the character of subcritical andronovhopf bifurcation that leads to the discharge hysteresis in type 2 mes v neurons.
【keywords】 mes v; neurons; electrophysiology; bifurcation; transition; i4ap
0 引言
hodgkin等根据单个神经元动作电位发放模式可把膜兴奋性分为两种“类型”［1］：ⅰ型兴奋性神经元放电阈值较低，并且随着刺激强度的增大，放电频率逐渐增大；ⅱ型神经元的放电频率相对集中在一个较小的范围内，并不随着刺激强度的变化而有显著的变化. 三叉神经中脑核（mesencephalic v,mes v）神经元是位于中枢的初级感觉传入神经元，是传递和调控口面部本体感觉信息的重要部位［2］. 我们在对三叉神经中脑核的 研究 中观察到了两种类型放电模式的存在. 同时发现在三角波电流刺激下，ⅱ型兴奋性神经元的放电模式存在一种物理学上所谓的滞后现象，表现为序列放电中,起始动作电位的阈强度大于结束时动作电位的阈强度［3］. 传统的生物学理论对以上现象难以做出满意的解释.
非线性动力学分析具有定量准确的特点，已广泛用于神经元放电，包括对mesv神经元放电的研究［4-7］，我们希望以此对所记录到的现象进行合理的分析. 分岔在非线性动力学理论中是指系统行为的定性变化［8］. 对于神经元来说，处于静息电位时是一种状态，当细胞持续放电，即膜电位大幅周期性振荡是另一种状态，细胞从静息状态向放电状态的转变，也就是膜电位从保持不变向周期振荡的转变，就是一种神经元性质的定性变化，反之亦然. 因此，可以将细胞所经历的这种转变称为细胞所经历的分岔. 本研究运用分岔理论对mesv神经元滞现象进行分析讨论，以突破传统电生理理论的局限.
1 材料和方法
1．1 材料 两周龄新生sd大鼠，体质量22～30 g，雌雄不限（第四军医大学实验动物中心）；1.4 mm内径记录电极玻璃毛坯（南京泉水教学实验器材厂）；戊巴比妥钠，nacl，kcl，cacl2，mgcl2，nah2po4，nahco3，kglucose和glucose（天津市医药公司）；乳酸（lactic acid）,四乙酸(egta)，羟乙基哌嗪乙磺酸(hepes)，k2atp, nagtph和4氨基吡啶（4ap）（美国sigma公司）； 渗透压仪（model 501, 美国precision system公司）；振动切片机(1000 plus, 美国vibratome公司)；正置微分干涉相差显微镜(bx51wi, 日本olympus公司)； 膜片钳放大器（multiclamp 700a）,模拟数字转换器(digidata 1322a，美国axon公司)；微电极拉制仪(pp83型， 日本narishige instruments公司)；恒流泵(hl2型，上海沪西仪器厂).
细胞外液（人工脑脊液，artificial cerebrospinal fluid, acsf）配方为（mmol/l）：124 nacl, 26 nahco3, 3 kcl, 1.25 nah2po4, 2 cacl2, 2 mgcl2, 10 glucose；切片用液（cuting solution）配方为（mmol/l）：126 nacl, 26 nahco3, 3 kcl, 1.25 nah2po4, 1 cacl2, 5 mgcl2, 10 glucose，以上两种液体用1 mol/l hcl将ph值调为7.3左右，渗透压为290～310 mosm/l；记录内液配方（mmol/l）：9 nacl, 115 kglucose，1 mgcl2，10 hepes，0.2 egta，3k2atp，naatp用naoh将ph值调为7.2左右，渗透压为298 mosm/l，离心管分装冻存备用；4氨基吡啶（4ap）用acsf配成0.5 mmol/l分装离心管冻存，用时用acsf稀释10倍.
1．2 方法
1．2．1 大鼠脑片制作 用戊巴比妥钠（40 mg/kg,ip）麻醉后，将大鼠迅速断头，剪开头皮，咬骨钳咬开颅骨，取出全脑置于含950 ml/l o2和50 ml/l co2饱和的冰冷的切片用液. 2 min后取出放在浸有acsf的滤纸上将全脑修块，然后用502胶水粘在切片槽的固定位置. 槽中放有冰冷的acsf（小于4℃），并不断通入含950 ml/l o2和50 ml/l co2混合气体. 用振动切片机切出含mesv在内的300 mm的冠状脑片，切好的脑片置于25℃下含950 ml/l o2和50 ml/l co2饱和气体的acsf中孵育30 min［7］.
1．2．2 脑片固定与灌流 将脑片移至记录用灌流槽内，浸没于液面下l～2 mm，用铂金丝制成u形框架，将尼龙丝粘于其上，制成盖网轻覆于脑片上. 使用恒流泵持续向脑片浴槽灌注含950 ml/l o2和50 ml/l co2饱和气体的acsf，流速为2～3 ml/min，灌流槽温度保持在26～30℃.
1．2．3 全细胞膜片钳记录 玻璃微管洗净烘干后由微电极拉制仪拉制成尖端直径约为1～2 μm，充灌电极内液后的入水阻抗为3～5 mω的玻璃电极. 实验在正置红外微分干涉显微镜下进行，用5×物镜观察定位mes v，40×水镜确认细胞，选取大于35 mm者进行可视巨阻封接. 成功后采用口吸破膜形成全细胞模式.电容和液接电位（约5 mv）线下补偿，记录膜被动参数，选取静息电位大于-55 mv,动作电位幅值大于70 mv进行实验.刺激所用三角波电流用clampex 9.2编写，刺激强度为600 pa，持续时间600 ms.
统计学处理: 实验中所有电刺激均从胞内给予，模拟数字转换器采样频率10 khz，低通滤波频率为3 khz；数据采集、记录采用axon laboratory提供的软件包clampex 9.2进行，分析软件采用clampfit 9.2，绘图和统计软件采用origin 6.1. 所有数据用x±s表示，用spss 11.0进行配对t检验.
2 结果
随机选取34例mes v神经元统计基本电生理参数：静息膜电位vm=（60.3±2.5） mv，膜电阻rm=（82.7±26.5） mω,膜电容cm=（68.0±12.9） pf.
2.1 两时象斜波刺激下mesv神经元放电的滞后现象 选取ⅱ型放电模式神经元，在三角波电流刺激下此类神经元表现出明显的滞后现象，即起始动作电位阈强度i1大于结束动作电位的阈强度i2，结束的动作电位的时刻较起始动作电位的时刻滞后（图1）. 根据统计学 分析 此类神经元i1和i2，阈电位v1和v2的差异有统计学意义(p<0.01, 表1). 4ap能够抑制mes v神经元上低电位激活的钾电流(i4ap)从而改变细胞的放电类型，灌流加入4ap(50 μmol/l)神经元的放电类型从ⅱ型转型为ⅰ型，同样的刺激方式在已转为ⅰ型放电模式的神经元上没有出现滞后现象（图2）. 根据统计学分析此类神经元i1和i2，v1和v2的差异无统计学意义（p>0.05, 表1）.
图1 缓慢三角波电流刺激使ⅱ型兴奋性神经元放电出现滞后现象（略）
表1 ⅰ型和ⅱ型放电神经元起始与结束放电阈强度及阈电位水平的比较（略）
bp<0.01 vs起始放电.
2.2 两种类型放电神经元动力学分岔图 利用已有的三叉神经中脑核神经元放电的模型资料［2］，可得ⅱ型（图3a）和ⅰ型（图3b）放电神经元放电的动力学分岔图.
图2 4ap作用下神经元转变为一类放电模式，滞后现象消失（略）
a: 为ⅱ型神经元的亚临界霍普夫分岔图，方框内为双稳态区; b: 为ⅰ型神经元的鞍结不变环分岔图. hb为分岔点，实线表示稳定的不动点，虚线表示不稳定的不动点，实心圆圈表示稳定周期解的最大值与最小值，虚心圆圈表示不稳定周期解的最大值与最小值.
图3 两类放电神经元动力学分岔图（略）
3 讨论
传统电生理 理论 把神经元兴奋性定义为产生动作电位的能力，并从阈值、幅度、频率、时程、动作电位的具体形态等来对神经元的放电进行描述和分类. 本文中的分类方式即是采用频率对刺激强度的反应作为标准. 然而此种理论和 方法 仅能对丰富的放电现象进行定性的描述. 例如传统理论认为可兴奋细胞存在一阈值，当刺激强度达到阈值时爆发动作电位，低于阈值不会爆发动作电位. 本实验中观察到的两种类型兴奋性神经元对相同类型刺激产生不同的放电模式，其中ⅱ型兴奋性神经元所产生的滞后现象用传统的定性的电生理理论关于阈值的定义则无法做出解释.
guttman等［3］利用hh方程 研究 枪乌贼巨轴突重复放电时，从解得的极限环图中曾预测轴突的放电存在滞后现象（hysteresis）,并在实验中得到证实.本实验在mes v神经元上也观察了此种现象，并且观察到加入4ap使神经元兴奋类型转变后滞后现象消失.
非线性动力学理论认为：产生不同兴奋性类型的原因是系统经历了不同类型的分岔［9］. 当系统经历鞍结不变环分岔时，产生ⅰ型兴奋性；而当系统经历霍普夫分岔或鞍结分岔时，产生ⅱ型兴奋性. 我们对mesv神经元的研究发现mesv神经元的放电属于亚临界霍普夫分岔，并且在施加药物4ap后这种亚临界霍普夫分岔的性质会被掩盖而转型为鞍结不变环分岔［2］. 利用分岔理论对滞后现象的产生和消失进行解释十分直观.
从分岔图上（图3a）可以很清楚的看出，ⅱ型兴奋性神经元在放电的起始区有一个很明显的双稳态区，在刺激电流由小到大和由大到小连续缓慢变化的过程中，在两个稳态之间跳跃的临界点不同，也就是所达到的刺激电流的大小不同. 当刺激电流由小到大时，膜电位只有到达hb点时产生跳跃（分岔），而当刺激电流由大到小时膜电位因仍处于稳态区会越过hb在跳跃（分岔）. 表现出动作电位起始的阈强度大于其结束时的阈强度，于是产生了滞后现象. 也可以说滞后现象就是亚临界霍普夫分岔的重要性质. 在ⅱ型兴奋性神经元上加入钾电流阻断剂4ap后，由于i4ap的减小它的动力学过程由ⅱ型转变为ⅰ型. 其动力学分岔图也发生了改变（图3b）. 从ⅰ型兴奋性神经元的动力学分叉图上可以看出，不存在一个双稳态区，膜电位会随着刺激电流的增大见效在hb点处跳跃（分岔），不能产生滞后现象. 借助数学模型，通过动力学分析，传统理论无法解释的神经元放电现象得以了合理的说明.
【 参考 文献 】
［1］ hodgkin al. the local electric changes associated with repetitive action in a nonmedullated axon［j］. j physiol, 1948,107:165-181.
［2］ liu y, yang j, hu s. transition between two excitabilities in mesencephalic v neurons［j］. j comput neurosci, 2007：/content/g7378651363q3776/fulltext.pdf
［3］ guttman r, lewis s, rinzel j. control of repetitive firing in squid axon membrane as a model for a neuroneoscillator［j］. j physiol, 1980,305:377-395.
［4］ del negro ca, chandler sh. physiological and theoretical analysis of k+ currents controlling discharge in neonatal rat mesencephalic trigeminal neurons［j］. j neurophysiol, 1997,77:537-553.
［5］ tanaka s, wu n, hsaio cf, et al. development of inward rectification and control of membrane excitability in mesencephalic v neurons［j］. j neurophysiol, 2003,89:1288-1298.
［6］ wu n, enomoto a, tanaka s, et al. persistent sodium currents in mesencephalic v neurons participate in burst generation and control of membrane excitability［j］. j neurophysiol, 2005,93: 2710-2722.
［7］ wu n, hsiao cf, chandler sh. membrane resonance and subthreshold membrane oscillations in mesencephalic v neurons: participants in burst generation［j］. j neurosci, 2001,21:3729-3739.
［8］ alligood k, sauer t, yorke j, et al. an introduction to dynamical systems［m］.new york:springerverlag,1996: 17-21.
［9］ koch c, segev i. methods in neuronal modeling: from ions to networks［m］. london: the mit press cambridge,1998：251-291.
