刺激·参数设置
$\quad\quad\quad\quad - g_K n^4(V-E_K) - g_L(V-E_L) + I$
$m,h,n$: 通道门变量
基于霍奇金-郝克斯利模型的动作电位数值模拟。通过改变刺激电流和通道导电性,实时计算离子通道开关和膜电位变化。
本模拟器基于霍奇金-郝克斯利模型,用数值积分描述神经细胞的兴奋性,计算膜电位随时间的演变。膜容量为\(C_m\),各离子通道的导电性为\(g_{Na}, g_K, g_L\),相应的平衡电位为\(E_{Na}, E_K, E_L\),膜电位\(V\)的变化由以下方程给出: $$ C_m \frac{dV}{dt} = I_{ext} - g_{Na}m^3h(V-E_{Na}) - g_K n^4(V-E_K) - g_L(V-E_L) $$ 其中\(m, h, n\)是门变量,它们遵循由电压依赖速率常数\(\alpha, \beta\)确定的以下微分方程: $$ \frac{dm}{dt} = \alpha_m(1-m) - \beta_m m $$ 用户可以用滑块调整刺激电流\(I_{ext}\)和最大导电性值,通道开关概率和膜电位波形会实时更新。这使得用户能够直观地理解动作电位的发生阈值和不应期的机制。
产业中的实际应用案例
制药行业使用神经元·动作电位模拟器在开发局部麻醉药(利多卡因等)和抗癫痫药(苯妥英等)时,提前评估钠通道阻断效果。医疗设备制造商(如美敦力公司)利用它来优化神经刺激装置(脊髓刺激疗法)的参数,通过模拟刺激电流和通道反应的关系,实现了动物实验的减少和开发周期的缩短。
研究和教育应用
在大学生理学和神经科学课程中,学生通过改变通道导电性和刺激强度,实时观察膜电位变化,直观理解动作电位的发生机制(去极化、复极化、不应期)。在研究中,它被用于再现帕金森病和癫痫的异常放电模式,有助于探索新的治疗靶点。
与CAE分析的联系和在实务中的定位
本模拟器作为生体电气现象的CAE(计算机辅助工程)分析工具而得到定位。在实务中,它被纳入神经界面设计(人工耳蜗和脑深部刺激电极)的预验证流程,通过数值分析电极形状、位置与神经响应的关系。此外,它还与心脏电生理模拟器相连联,用于药物诱发QT延长综合征的风险评估等,作为多尺度生体模拟的基础而被应用。
人们通常认为"增大刺激电流,动作电位的振幅也会成比例增加",但实际上在霍奇金-郝克斯利模型中,动作电位的振幅由钠通道的平衡电位决定,一旦刺激电流超过阈值,振幅基本保持不变。过大的电流只会增加放电频率,而振幅会饱和。需要注意这一点。
人们通常认为"改变通道导电性只是简单地加快膜电位的时间变化",但实际上钠和钾导电性的比例会复杂地影响动作电位的波形和不应期。特别是过度增加钾导电性会导致复极化过快,动作电位变小;反之降低会导致不应期异常延长。参数之间的平衡很重要。
人们通常认为"模拟结果准确再现了生体本身",但实际上霍奇金-郝克斯利模型是基于鱿鱼巨轴突的近似模型,人类神经细胞和心肌细胞的通道类型和动态是不同的。此外,温度和细胞外离子浓度的影响也是固定的,因此定量解释时需要谨慎。
对以直径1mm的牛蛙巨轴突为对象,在iext=7.0μA/cm²、ton=10ms、tstim=1.0ms的刺激条件下,V_m初值-65mV在约13ms时因Na⁺内流达到-30mV阈值,动作电位发生。此时m值(Na⁺激活)从0.05急升至0.8,h值(Na⁺失活)从0.6下降至0.1。其后因K⁺外流,V_m复极化至-90mV,n值增加产生过极化。无刺激时,V_m收敛至-65mV的静息膜电位。