神经元兴奋性与Hodgkin-Huxley模型中的参数变异性

《PLOS Computational Biology》：Neuronal excitability and parameter variability in the Hodgkin-Huxley model

　　

生物物理细节丰富的神经元模型通常被构建为一种单向流程：电压钳数据被简化为单一的最优通道参数集，随后这些参数被组合成一个确定性的放电模型。这种做法丢弃了实验观察到的离散性和拟合不确定性，掩盖了兴奋性系统中鲁棒性（robustness）和简并性（degeneracy）产生的机制。本研究中，研究人员将拟合参数的不确定性重新引入Hodgkin-Huxley模型，并将不确定性分析与全局敏感性分析嵌入模型构建过程。研究人员对Hodgkin和Huxley原始文献中的钠和钾速率常数数据进行了数字化处理，并使用bootstrap重采样（bootstrap resampling）来估计拟合电压依赖性动力学参数的不确定性。随后，通过大规模蒙特卡洛模拟（Monte Carlo simulations），将这些不确定性估计传播到一个空间扩展的乌贼轴突电缆模型中，其中每个样本定义了一套完整的动力学、电导、被动和结构参数集。在通道水平，一阶Sobol敏感性指数（first-order Sobol sensitivity indices）显示所有动力学参数均以强烈的时间依赖性方式贡献于输出方差，不同参数分别控制钾和钠电导的瞬态和稳态行为。在神经元兴奋性水平，模拟产生了一个异质的放电行为群体，包括不放电、相位性（phasic）放电、规则放电和自发活动。跨刺激强度下，主导放电模式为刺激起始时的单脉冲放电，这与乌贼巨轴突在快速逃逸行为中的生理作用一致。经典的1952年Hodgkin-Huxley参数集落在规则放电的少数亚群体中，而非代表一个独特或主导的解。在相位性亚群体中，动作电位传播和传导速度变化很大但仍处于实验范围内。最后，放电期间的敏感性分析显示一阶Sobol指数均匀较小而总阶指数（total-order indices）较大，表明兴奋性主要由所有参数之间的强相互作用而非任何子集所控制。这些结果共同支持将Hodgkin-Huxley模型重新框架化为一个实验约束的行为集合体，而非单一的特权参数集，其中具有生理相关性的放电模式从参数空间的结构化区域中涌现。

本研究发表于《PLOS Computational Biology》，旨在探讨生物物理神经元模型构建中的核心方法论问题。当前，神经元建模领域存在一个长期存在的困境：随着模型细节不断增加，从实验继承的参数通常被视为固定值，尽管这些测量值本身受到变异性、不确定性和模型依赖性假设的影响。具体而言，Hodgkin-Huxley模型自1952年建立以来，其构建流程本质上是一种还原主义的单向管道——从电压钳记录到速率常数，再到电压依赖性方程，最终到放电神经元模型，每一步都将可变数据拟合为少量参数，而原始测量中的不确定性被压缩为点估计。这种做法的代价在通道动力学层面尤为突出：一旦电压钳数据被拟合为速率常数参数，这些参数就被视为模型的固定属性，原始测量的不确定性未能继续传递。由此产生的模型可能再现典型的放电模式，但无法自行揭示该模式是稳健的、典型的，还是更广泛可能行为家族中的一个成员。

值得注意的是，这一问题与已建立的离子通道电导简并性（degeneracy）研究有所区别。先前研究已表明，相似的神经元或网络活动可由多种不同的最大电导和其他细胞参数组合产生，参数简并性是传导性模型的特征而非伪迹。然而，在整合的Hodgkin-Huxley型放电模型中，通道动力学通常作为预先规定的模型结构被继承，而推断仅在电导和少量额外膜参数或适应参数上进行。因此，动力学参数的可变性在建模管道的构建阶段和分析阶段均未被充分利用——在构建阶段，通道动力学常被一次性拟合后冻结，作为固定支架；在分析阶段，变异性和敏感性分析通常在模型构建完成后才应用，而非作为模型构建本身的组成部分。

本研究的核心目的在于填补这一空白。研究人员首先利用bootstrap重采样方法，从原始Hodgkin-Huxley文献的图表中估计钠和钾通道动力学参数的不确定性，获得各拟合参数的经验分布。随后，通过大规模蒙特卡洛模拟，将这些不确定性传播到一个空间扩展的Hodgkin-Huxley电缆模型中，其中每个样本定义了包含动力学、电导、被动和结构参数的完整模型。最终，应用基于Sobol分解的全局方差敏感性分析，评估单个参数及其相互作用对通道层面行为和神经元兴奋性的贡献。

研究采用的关键技术方法包括：基于数字化原始文献数据的bootstrap重采样（对六个速率常数各进行4000次重采样）以获得参数的经验分布；构建空间扩展的一维圆柱形电缆模型（长10 cm，半径0.025 cm，80个等电位 compartments），采用密封端边界条件；使用Python（NumPy和JAX库）实现向量化并行的大规模蒙特卡洛模拟；以及采用基于Saltelli扩展的Sobol方法计算全局敏感性指数，包括一阶敏感性指数和总阶敏感性指数，以区分参数的独立效应与交互效应。

研究结果部分呈现如下几个核心发现。

**Bootstrap参数估计与通道水平敏感性分析**

研究人员对Potassium电导的五个动力学参数和Sodium电导的十个动力学参数进行了bootstrap估计。几乎所有15个参数均呈钟形分布，但控制Sodium电导激活反向速率的两个参数因原始数据点少且离散大而呈现异常。基于这些估计，对钾电导进行的电压钳模拟显示，五个动力学参数各自以独特的时间依赖性方式影响输出方差。一阶Sobol敏感性指数表明，参数A_α,n、V_1/2,α,n、z_α,n、A_β,n和z_β,n均对钾通道开放概率的变异性有贡献，且这种贡献随时间动态变化。类似地，Sodium电导的激活和失活过程的十个参数也表现出时间依赖性敏感性，激活相关参数在电压阶跃时呈现瞬态变化，而失活参数因过程慢十倍而表现出更慢的敏感性动态。

**放电行为群体的结构**

将不确定性传播至空间扩展的电缆模型后，300,000组参数的蒙特卡洛模拟产生了一个异质的放电行为群体。基于电压响应，这些轨迹被分为四个亚群体：不放电（仅亚阈值去极化）、相位性放电（刺激起始时单一动作电位）、规则放电（刺激期间多个动作电位）和自发放电（无注入电流时发生动作电位）。结果显示，7.2%的模拟表现出自发放电；随着电流增加，相位性和规则放电亚群体比例增长，亚阈值响应在8 μA时近乎消失。值得注意的是，在8 μA刺激下，规则放电亚群体仅占约14%，而大多数模拟（77.7%）产生相位性响应。这一结果出人意料，因为经典Hodgkin-Huxley模型通常产生规则放电模式，但更符合乌贼巨轴突在快速逃逸行为中的生理角色。

**规则放电亚群体的放电率变异性**

对规则放电亚群体的进一步分析显示，不同电流水平下该亚群体产生广泛的规则放电频率分布，其均值和标准差随注入电流变化。关键发现是：使用原始Hodgkin-Huxley参数集的模拟结果（红色符号）与规则放电亚群体的平均放电率高度吻合。这表明原始参数集并非不确定性约束群体中的主导行为，而是落在规则放电少数群体中，其经典放电模式仅代表参数空间的一个结构化区域，而非典型行为。

**相位性动作电位的传播与变异性**

鉴于相位性响应的主导地位及其生理相关性，研究人员进一步分析了该亚群体。500个代表性模拟显示，部分动作电位在电流注入位点成功起始但未能传播至远端记录位点（8 cm处），导致远端幅度分布广泛。成功传播定义为远端动作电位幅度达70 mV。传导速度（conduction velocity）计算为注入位点与远端记录位点间距离除以动作电位峰值时间延迟，结果范围为约5–15 m/s，与实验数据一致。

**兴奋性的全局敏感性结构**

对相位性放电亚群体的敏感性分析（基于182,386个成功发放动作电位的模拟）揭示了关键发现。一阶Sobol敏感性指数在动作电位期间（t=12 ms）和随后的持续膜电位期间（t=70 ms）均较小，表明单个参数在耦合微分方程积分后的独立效应有限。相对的，总阶Sobol敏感性指数在几乎所有参数上均较大，表明参数间存在强烈的相互依赖性。两个钾电导参数V_1/2,α,n和A_β,n对持续电位有可观的独立影响，这与它们在初始动作电位后去极化阻滞中的直接作用一致。这一模式在动作电位特征和持续膜电位两个时间点上均成立，表明交互作用主导的敏感性是整合兴奋性的稳健特征。

研究讨论部分深入分析了上述结果的科学与方法论意义。研究人员指出，参数敏感性强烈依赖于模型被 interrogated 的层次：在电压钳条件下，单个动力学参数对通道行为有显著的一阶效应；但当这些电导嵌入兴奋性膜并驱动发放时，情况发生质变——低阶敏感性并不意味弱影响，而是表明参数的影响依赖于其他参数的值。在Hodgkin-Huxley模型中，这种耦合通过膜电位介导：一个动力学、电导或被动参数的改变影响电压，进而改变所有电压依赖性速率，再反馈至电压。因此，从孤立通道动力学到整合兴奋性的转变，将参数影响从相对直接的控制转变为交互作用主导的控制。

关于放电行为群体的生理相关性，研究强调相位性放电的主导地位与乌贼巨轴突已知的生理功能高度一致——该结构触发快速、全或无的 mantle 收缩以产生逃逸行为，在自然条件下通常对感觉刺激产生单一动作电位。这与经典Hodgkin-Huxley参数集通常产生重复放电形成对比，表明原始参数集仅代表更广泛、生理相关的参数景观中的一个区域。

该研究与参数据有学（database approach）和模拟推断（simulation-based inference）等互补性工作形成对话。Prinz等建立的参数据方法表明相似的神经活动可从不同的最大电导组合产生，确立了神经元建模中多解性的概念；而模拟推断使用似然自由贝叶斯推断和神经密度估计来恢复生物物理参数的后验分布。本研究将同样的群体逻辑扩展至动力学参数本身——这一区分至关重要，因为最大电导主要作为比例因子，而动力学参数决定塑造发放起始、不应期和传播的时序和非线性状态转换。

研究结论部分明确指出：Hodgkin-Huxley模型仍然是理解神经元兴奋性的独特强有力基础，但当实验约束的拟合速率参数的不确定性被保留并传播时，该模型产生了一个以相位性响应为主导的结构化放电行为群体，这与乌贼巨轴突的生理角色一致，而经典的1952年参数集作为规则放电的少数群体置于其中。该群体中的兴奋性由动力学、电导和被动参数之间的相互作用而非任何单一参数所控制。因此，Hodgkin-Huxley模型最好被理解为一个实验约束的行为集合体，而非单一的特权参数集。这一发现支持将生物物理建模工作流重新框架化：首先执行不确定性分析，然后生成作为结构化群体分析的大型模型集合，其中生理行为对应于参数空间的区域而非孤立点。选择单一"最优"模型不可避免地丢弃了这一结构，类似于过拟合。