钠的正反馈与两种刹车
把三个门控方程和膜方程联立,得到四个变量的常微分方程组。电压稍有去极化,快速的钠激活就跟上来,钠电流往里流,电压上升更快。这一步没有刹车,电压一路冲向钠的反转电位。
刹车来自两个慢过程。钠失活变量 在一毫秒左右把钠通道关上,钾激活变量 同时把钾通道打开。电压于是回落,并且冲过静息电位,在静息以下约 11 mV 处停一会儿,再慢慢回到静息。
Hodgkin 和 Huxley 按 Hartree 的方法一步步数值积分这组方程(论文第 523 页),上升段的步长是 0.02 ms。我们用四阶 Runge–Kutta 积分,步长减半时峰值变化小于 。
阈值在 6 和 7 nC/cm² 之间
Hodgkin 和 Huxley 在 6 °C 下算了一组初始去极化:6 mV 不放电,7 mV 放电(论文图 12、图 21)。我们的模型在 6.3 °C 下用二分法找到阈值 ;另一套按论文原始符号独立写的积分器给出 6.504。
超过阈值以后,电击再大,峰值也差不多。论文表 4 列出:7 mV 起点给出 102.1 mV 的峰,15 mV 给出 105.4 mV,90 mV 给出 108.5 mV。起点差了十几倍,峰值只差 6%。这就是“全或无”。
阈值不是一个固定的电压。刚放过一次电,钠通道大多处于失活状态,钾通道还开着,再来一次同样的电击也不放电,这是不应期(论文图 20)。反过来,先把膜超极化 30 mV 再突然放开,钠失活解除了一部分,钾通道又关着,膜会自己放一次电,这叫阳极断开(论文图 22)。
完整例题:一次 15 nC/cm² 电击
题目。6.3 °C 下给静息膜一次 15 nC/cm² 的瞬时电击。求峰值、最大上升速率,以及这次冲动额外带进的钠和带出的钾,并与 1952 年表 4、表 5 比较。把电容充到峰值需要多少钠,与实际进入的钠相比如何?
建模。从 出发积分式 (26)。流量按论文的做法积分到峰后第三次回到静息电平为止,再减去同一时段内静息时本来就有的流量;减去静息流量这一步是我们根据数值推断的,论文没有明说。
计算。峰值 ,表 4 印 105.4;最大上升速率 ,印 311;电导峰值 ,印 37.0。额外进入的钠 ,表 5 印 14.46;额外流出的钾 14.33,印 14.32。把 的电容充到 105 mV 只需 。
解释。实际进入的钠是充电所需的 13 倍。上升沿后半段钾通道已经打开,钠流进来的电荷大部分被同时流出的钾抵消,只有一小部分留在电容上。峰电位在电学上很省,在离子上却很费,后面泵要把这些钠钾搬回去。
检验。把同样的计算放到 18.5 °C,模型给出 96.9 mV 和 563 V/s,表 4 印 96.8 和 564;额外的钠、钾为 4.01 和 4.10 pmol/cm²,表 5 印 3.99 和 4.07。阳极断开那一行,钾的差别最大,模型 16.59、论文 16.77,差 1.1%。