神经元分叉景观中的简并性：跨不同参数集的等效分叉序列

《PLOS Computational Biology》：Degeneracy in neuronal bifurcation landscapes: Equivalent bifurcation sequences across distinct parameter sets

　　

简并性（Degeneracy），即结构不同的系统实现相似功能的能力，是生物体的基本属性，赋予其适应性与韧性。神经元虽在离子通道表达上存在显著变异，却能维持稳定的活动模式，凸显了不同内在构型可产生等效电生理行为。既往研究表明，神经元活动源于内在膜特性与离子浓度动态之间的相互作用。然而，简并性在神经元动态景观中如何组织，以及不同简并构型如何响应外部扰动，仍不清楚。本研究中，研究人员探究了不同的生物物理参数集如何产生等效分叉序列，以及它们的韧性如何在离子调制下发生变化。研究人员采用具有三个可变电导的极简神经元模型，系统探索细胞外钾浓度（代表生理与病理条件）的变化如何重塑神经元动力学。结果表明，不同的电导构型占据参数空间中结构化的区域，支持相同的分叉路径，从而产生等效活动模式，如爆裂放电（bursting）、癫痫样活动（seizure-like activity）和去极化阻滞（depolarization block, DB）。然而，这些简并构型对外部扰动表现出不同的敏感性，揭示出固定条件下的功能等价并不意味着等效的稳健性。为量化这一效应，研究人员引入韧性（resilience）作为度量指标，用以衡量在外部离子变化下，给定电导构型转变为特定活动模式的阈值。结果表明，简并性超越了多种解的共存：简并构型形成了结构化的景观，在环境调制下呈现不同的轨迹。该框架为研究内在变异如何影响神经元在不同动态状态下的脆弱性与韧性提供了一种定量方法。

一、研究背景

简并性（Degeneracy）指结构不同的系统产生相似功能结果的能力，是生物系统的基本属性，使神经元在离子通道表达显著变异时仍能维持稳定放电活动。Eve Marder及其同事的开创性研究表明，单个神经元可补偿内部参数（如离子电导）扰动以维持稳定放电。然而，既往研究集中于稳态或重复放电活动，简并性如何体现在活动模式间的转变——即分叉（Bifurcation）——中仍不清楚。癫痫中，神经元兴奋性常从静息、爆裂放电（bursting）、癫痫样放电（seizure-like activity）突变为去极化阻滞（depolarization block, DB），可数学描述为分叉。尽管癫痫样电生理特征在不同物种和病因中高度守恒，但核心问题未解决：不同离子电导组合能否在环境变化下产生相同分叉序列？细胞外钾浓度（[K]_bath）升高时简并性如何表现？固定条件下的功能等价是否意味着对扰动具有等效稳健性，亦不清楚。

二、研究开展与意义

研究人员采用单房室（single-compartment）神经元模型，含钠、钾、氯三个可调电导，系统改变三个最大电导与细胞外钾浓度，绘制活动模式与分叉序列图谱。研究得出：不同电导组合占据参数空间中结构化的紧凑区域，支持相同分叉路径，产生等效活动模式；但固定条件下功能等价的简并构型对钾浓度变化表现出不同敏感性，表明功能等价并不意味着等效稳健性。研究人员引入韧性（resilience）指标，量化外部离子变化下电导构型转变为特定活动模式的阈值。该发现将简并性分析扩展到动态景观层面，为理解癫痫治疗差异与药物抵抗提供基础，并提示仅依据分子标记分类神经元并不充分。

三、关键技术方法

研究人员开发仅基于膜电位（V_m）的分叉特征确定性分类算法，通过基线跳变、峰电位振幅比等特征区分SNIC、SN-SH、SN-SupH等分叉类型；采用Runge-Kutta法积分四个耦合常微分方程构成的快慢系统，在三个最大电导与细胞外钾浓度构成的参数空间中扫描，生成10,000组模拟实现作为样本队列；计算各参数点诱发目标模式的最小钾浓度与最近邻的差值构建韧性向量，以其归一化投影评估各电导对韧性的贡献。

四、研究结果

2.1 不同参数化下的相似动力学：模型可重现多种癫痫相关爆裂模式，分叉类型可基于膜电位特征区分；通过比较两个参数集的输出并扩展至系统扫描，生成10,000个不同模型实现，研究简并性的涌现与变化。

2.2 简并参数与共享模式：产生相同模式的参数集紧密聚集，形成紧凑简并区域（degeneracy region）；固定钾浓度下观察到多达八种活动模式。随钾浓度升高，简并区域保持紧凑但分布随模式类别而异，新模式出现而另一些消失，表明简并性组织结构对外部条件敏感。

2.3 外部条件塑造神经元响应：随钾浓度升高，各实现表现出不同数量的活动模式，表达相同模式数量的参数集不一定形成紧凑区域。不同电导配置即使在基线条件下产生相同模式，对钾浓度变化的响应也可能不同；即使模式轨迹相同，转变发生的钾浓度水平也可能不同，提示每个电导空间中的点对外部扰动具有独特敏感性。

2.4 韧性反映离子通道构型：氯电导（g_Cl）对韧性贡献最大，钾电导（g_K）次之，均为正贡献；钠电导（g_Na）贡献较小且为负。这归因于：钾浓度动态显式建模使钾反转电位（E_K）随钾升高正向偏移，削弱了钾驱动力近分叉点的作用；氯浓度恒定使氯反转电位（E_Cl）不变，氯电导变化对电压零线产生更大位移。钠变化因电中性假设与钾耦合、方向相反，故贡献为负。

五、讨论部分

讨论部分指出，分类方法与简并区域识别均仅依赖膜电位，对实验数据具有重要意义，因为实验记录通常只能获取膜电位。研究人员将简并性分析扩展到动力学转变层面，发现产生等效分叉序列的电导构型占据参数空间中拓扑紧凑的区域，其中某些电导被紧密约束，其他电导变异较大。升高钾浓度模拟癫痫病理驱动，虽然简并性保持，但支撑各模式的内部电导谱系发生重排。研究表明简并非静态属性，而是与结构化的变异性共存。

六、研究结论

研究人员提出在神经元活动模式中识别和表征简并性的框架，特别是与癫痫动力学相关的模式。该方法仅凭膜电位即可分类爆裂模式，并引入韧性指标量化内部构型对细胞外扰动的影响。通过改变简化模型中的电导，研究人员证实了固定外部条件下多种分叉类型间存在简并性，并证明简并构型在基线条件下支持相似模式，但在钾浓度变化时向不同模式的转变存在差异。这种对模式稳健性与敏感性的双重分析为研究神经元内在变异提供了新途径。这些发现启示：仅凭分子标记分类神经元并不充分，其扰动响应提供更丰富的分类依据；在癫痫中，单神经元韧性的变异可能有助于解释治疗反应异质性和药物抵抗。本方法为这类研究提供原则性基础，并为跨尺度整合神经元变异性的建模迈出一步。