膜电位·动作电位模拟器 返回
生物工程

膜电位·能斯特方程式·动作电位模拟器

快速解答
电池的标准电动势为 E°cell = E°cathode − E°anode,E°cell > 0 时反应自发(原电池)。浓度效应由能斯特方程 E = E°cell − (0.0592/n)·log₁₀Q 描述(25℃,n=转移电子数、Q=反应商);放电过程中 Q 增大、E 下降,Q→K 时 E→0(达到平衡)。

使用能斯特方程式和GHK方程式计算离子平衡电位和膜电位。使用霍奇金-赫胥黎简化模型实时模拟动作电位波形和阈值效应。

离子浓度·透过性
细胞类型
K⁺ 浓度 (mM)
[K]_in
mM
[K]_out
mM
Na⁺ 浓度 (mM)
[Na]_in
mM
[Na]_out
mM
Cl⁻ 浓度 (mM)
[Cl]_in
mM
[Cl]_out
mM
透过性比 P_Na/P_K
透过性比 P_Cl/P_K
温度 T
°C
注入电流 I_stim
μA/cm²
计算结果
E_K [mV]
E_Na [mV]
E_Cl [mV]
静息膜电位 V_m [mV]
动作电位峰值 [mV]
动作电位持续时间 [ms]
阈值电流 [μA/cm²]
细胞膜示意图
动作电位波形 V_m(t)
理论·主要公式

$$E_{ion} = \frac{RT}{zF}\ln\frac{[\text{ion}]_{out}}{[\text{ion}]_{in}}$$

能斯特方程式。R:气体常数、T:绝对温度 [K]、z:离子价数、F:法拉第常数(96485 C/mol)

$$V_m = \frac{RT}{F}\ln\frac{P_K[K]_{out}+P_{Na}[Na]_{out}+P_{Cl}[Cl]_{in}}{P_K[K]_{in}+P_{Na}[Na]_{in}+P_{Cl}[Cl]_{out}}$$

Goldman-Hodgkin-Katz(GHK)方程式。P_K·P_Na·P_Cl:各离子的膜透过系数(相对比)

$$C_m\frac{dV}{dt} = -g_{Na}(V-E_{Na}) - g_K(V-E_K) - g_L(V-E_L) + I_{stim}$$

H-H简化模型。C_m:膜电容 [F/m²]、g:导纳 [S/m²]、I_stim:刺激电流 [A/m²]

膜电位·动作电位基础

🙋
「静息膜电位」是什么?教科书上说「细胞内侧比外侧负约-70mV」,但这个电位是怎么形成的呢?
🎓
简单说,是因为细胞膜内外的离子浓度不同。比如,钾离子(K⁺)在细胞内多,钠离子(Na⁺)在细胞外多。试试调动左边的「[K]_in」和「[K]_out」滑块。当浓度差变大时,计算出的「K⁺的平衡电位」的绝对值也会变大。
🙋
但是K⁺和Na⁺都存在,为什么整体的膜电位会接近K⁺的平衡电位(约-90mV)而不是-70mV呢?
🎓
好问题!其实细胞膜对K⁺有高度选择透过性,而对Na⁺几乎不透。这个「通透程度」就是「透过性比 P_Na/P_K」。试试把P_Na/P_K设置为0(只有K⁺通过),膜电位就会等于K⁺的平衡电位。但实际上有一点Na⁺也透过(P_Na/P_K=0.01左右),所以膜电位从-90mV向正方向偏移到-70mV左右。调动中间的透过性滑块试试。
🙋
那「动作电位」那急剧的电位上升,是不是这个透过性瞬间反向了?
🎓
完全正确!试试把下面的「注入电流 I_stim」调大,把膜电位推向阈值(约-55mV)。那时候,电压依赖性Na⁺通道会一下子打开,P_Na/P_K会临时跳到100以上。这样膜电位会快速接近Na⁺的平衡电位(约+60mV),形成那个尖锐的脉冲(除极化)。之后K⁺通道打开,膜电位恢复(复极化)。这整个过程在我们的模拟器上会实时演示出来。

常见问题

能斯特方程式和GHK方程式的结果为什么不同?
能斯特方程式计算的是单个离子的平衡电位,而GHK方程式考虑了多个离子的透过性。实际膜电位受多种离子影响,所以GHK方程式更符合生理学。可以调整透过性比(P_K:P_Na等)来比较。

模拟中动作电位的阈值为什么会变化?
阈值取决于膜的初始状态和离子通道参数(最大导纳、时间常数等)。比如,升高外部K⁺浓度会使静息膜电位除极化,更容易达到阈值。调整各参数观察波形变化。

温度滑块改变膜电位的原因是什么?
能斯特方程式和GHK方程式都含有绝对温度T。温度升高会改变RT/F项,导致同样浓度梯度下平衡电位的绝对值增大。霍奇金-赫胥黎模型也会因温度改变通道开闭速率。

为什么Cl⁻离子的价数设为-1?
Cl⁻是阴离子,电荷为-1,所以在能斯特方程式中要用z=-1。这会导致对数项符号反转,使得外部浓度高时平衡电位为负值。选择离子时务必确认价数正确。

实际应用

神经工程·深部脑刺激(DBS):帕金森病等疾病的治疗中,在脑深部放置电极进行电刺激。利用像本模拟器这样的神经细胞电学模型,设计最优刺激波形(频率、幅值、脉冲宽度),既能抑制症状又能最小化副作用。

人工内耳·视网膜芯片:将声音或光信息转换成电脉冲直接刺激听神经或视网膜神经。通过本工具的基础模型,可验证什么样的刺激模式最能可靠且自然地诱发神经活动电位。

心脏电生理·心律不齐模型:心肌细胞也会产生动作电位。通过建立包含Na⁺、K⁺、Ca²⁺等通道动态的模型,可用计算机重现多个细胞级联兴奋产生搏动的过程,以及节律失常的机制。

药物开发·药理试验:很多药物(尤其是心律不齐治疗药和局部麻醉药)作用于离子通道。通过膜片钳实验测得的电流数据与本模拟器模型对照,可推断药物对哪个通道的透过性或动态特性的影响。

常见误区和注意事项

开始使用这个模拟器时,容易踩的坑有几个。首先是「静息膜电位总是比K⁺平衡电位更正」这个思维定式。在神经元确实如此,但试试把P_Cl/P_K调大试试。Cl⁻透过性变高时,膜电位会被E_Cl(通常-70mV左右)拉动。有时候计算出的静息膜电位甚至可能比E_K还负。用工具玩一玩,打破「应该是这样」的先入为主。

其次是滑块变化对膜电位的「灵敏度」不均匀。比如,在静息状态(P_Na/P_K=0.01)时,把[K]_out从5mM改到10mM,膜电位会从-80mV大幅跳到-65mV。但同样改变[Na]_out却几乎没动静。这是因为膜对K⁺选择性高。调参数时要时刻意识到哪个离子起支配作用。

最后,「动作电位波形不随注入电流强度改变」这个错误认知。实际试试就能看出,I_stim恰好在阈值时,脉冲发火延迟,波形平缓。强电流一下子打进去,发火迅速,脉冲尖锐。虽然模拟器用的是简化模型,但这个现象和真实神经中一样。通过工具亲身体验,「阈值」并不是绝对的值,而是取决于刺激方式的「事件」。

使用指南

  1. 输入钾(K)和钠(Na)的细胞内外浓度。在kInNum/kOutNum、数值输入框/naOutNum中设定细胞内外的mM浓度
  2. 使用能斯特方程式 E = (RT/zF)×ln([离子]外/[离子]内) 自动计算各离子的平衡电位。体温37°C、F=96485 C/mol条件下进行计算
  3. 用GHK方程式求出静息膜电位,通过霍奇金-赫胥黎模型参数(Na/K-ATP酶活动、Na·K·Cl透过性)开始模拟。实时显示动作电位波形、峰值、持续时间

具体计算示例

神经元的标准条件:细胞内K 140mM、细胞外K 5mM、细胞内Na 10mM、细胞外Na 145mM。计算出能斯特平衡电位 E_K = -90mV、E_Na = +60mV。将相对透过性比设为Pk:PNa:PCl = 1:0.03:0.1,用GHK方程式得到静息膜电位 V_m = -70mV。在霍奇金-赫胥黎模型中注入0.5μA/cm² 刺激电流,观察到动作电位峰值上升到+30mV,持续时间约1.2ms

实务中的注意事项

  1. 低K外液(2mM以下)条件模拟时,E_K超过-110mV会导致数值计算不稳定,需将时间步设为0.01ms以下
  2. 温度与37°C差异大时,RT/F系数会改变(25°C时59mV/decade,37°C时61.5mV/decade),需准确设定实验条件
  3. 膜片钳或全膜片记录实验测得膜电位与计算值有偏差时,应检查Shilling补正和Leak补正影响。使用KCl填充电极时还需加上液接合电位补正(通常+15mV)