《GLIA》：A Computational Framework for Myelin as a Regulator of Temporal Integration, Neuronal Excitability, and Synaptic Plasticity
编辑推荐：
髓鞘不仅通过实现快速传导来调节神经信号传递，还通过动态调整轴突传导速度（conduction velocity, CV）进而塑造锋电位时序。虽然活动依赖性髓鞘形成被认为可以优化神经通信，但髓鞘变化如何转化为单细胞动力学和可塑性的改变仍不清楚。采用计算建模方法，
髓鞘不仅通过实现快速传导来调节神经信号传递，还通过动态调整轴突传导速度（conduction velocity, CV）进而塑造锋电位时序。虽然活动依赖性髓鞘形成被认为可以优化神经通信，但髓鞘变化如何转化为单细胞动力学和可塑性的改变仍不清楚。采用计算建模方法，研究人员表明轴突CV的破坏会重组突触前输入相关性，导致突触后兴奋性、兴奋-抑制平衡和脉冲时序依赖可塑性（spike-timing-dependent plasticity, STDP）的系统性变化。通过涵盖与健康髓鞘化、脱髓鞘和过度髓鞘化相关的CV区间，模拟结果表明髓鞘化不足和过度均可能诱导适应不良的神经元动力学。将放电率稳态（firing-rate homeostasis, FRH）纳入模型进一步提示急性和慢性髓鞘破坏具有不同的神经元动力学后果。最终，这些结果表明轴突CV可能是连接髓鞘可塑性与神经计算的关键机制纽带，为解释健康和疾病状态下髓鞘失调的功能后果提供了原则性框架。
一、研究背景与问题提出
髓鞘包裹轴突膜并通过跳跃式传导显著提高传导速度（conduction velocity, CV），从而补偿不同轴突在长度和直径上的差异。少突胶质细胞（oligodendrocytes, OLs）的髓鞘产生受神经元活动调控，这一现象称为活动依赖性髓鞘形成（activity-dependent myelination, ADM），被认为可通过调整动作电位沿轴突传播的CV来优化神经通信。多个动作电位到达时间的变异性通常对网络运作有害，而ADM对轴突CV的调节可标准化传导延迟，促进突触前动作电位同时、同步地到达靶神经元。
然而，关于髓鞘结构性变化对单个神经元和群体动力学的影响，尤其是突触后兴奋性和突触可塑性如何受到轴突CV变化的影响，目前知之甚少。量化髓鞘变化对单个神经元动力学影响的一种方法是考察轴突CV如何影响其突触前输入的统计相关性。突触前输入相关性反映了传入动作电位时序的同步程度：若锋电位同时或以一致的时间关系到达突触后细胞，则这些锋电位序列更具相关性。无论其起源如何，髓鞘对轴突CV的调节在塑造此类相关性中发挥重要作用，并对神经通信和时间整合具有重要意义。
为了表征髓鞘变化对单细胞动力学的影响，有必要考察在健康和病理状态下产生的不同轴突CV范围内突触前相关性如何变化。在健康脑中，髓鞘对突触前相关性的调节发挥着重要功能：调控神经信号传递、突触后兴奋性并促进突触可塑性。然而在脱髓鞘疾病中，髓鞘减少会抑制相关性，影响时间整合并导致细胞兴奋性的适应不良性转变。此外，通过限制网络调节其CV的能力，脱髓鞘可能阻碍轴突传导延迟的正常化，使脉冲时序依赖可塑性（spike-timing-dependent plasticity, STDP）受损。相反，异常增加的过度髓鞘化可能减少锋电位时序分散，放大突触前相关性和突触后放电率，带来潜在的病理后果。
由于髓鞘化和神经活动的变化发生在截然不同的时间尺度上——神经活动在毫秒到秒级别运作，而髓鞘化则在数小时到数天内发展——因此难以解释这些变化。在较长时间尺度上，细胞代偿机制如放电率稳态（firing-rate homeostasis, FRH）发挥重要作用。FRH指神经元通过一种或多种活动依赖性控制机制适应其兴奋性以正常化自发活动，当神经元持续暴露于异常强度的刺激时被激活。
本研究扩展了现有框架，在细胞水平建模轴突髓鞘化对神经动力学的影响，探索轴突髓鞘化的扰动（以轴突CV的变化来操作化）如何影响时间整合、神经元兴奋性和突触可塑性等一系列神经元层面动力学。研究结果揭示了髓鞘可塑性在细胞尺度上的多种潜在动力学特征，对健康和病理脑状态下的髓鞘可塑性参与具有启示意义。
二、主要研究方法
研究人员采用了稀疏连接的泊松（Poisson）神经元网络模型，设定在波动驱动（fluctuation-driven）区间运行。模型中每个突触后神经元接受来自一群兴奋性或抑制性细胞的汇聚性突触前输入，突触后膜电位（postsynaptic membrane potential, PSP）由无量纲化的随机微分方程描述。髓鞘的影响通过其对轴突平均CV的作用来建模，CV范围基于小鼠、兔和猫皮层及胼胝体有髓纤维的实验数据确定为0.3-4 m/s，低于此范围反映髓鞘化不足，高于此范围反映过度髓鞘化。轴突长度从伽马分布（gamma distribution）中采样，范围约0-30 mm。放电率稳态通过使阈值参数适应放电率变化来实现，阈值遵循特定的适应定律。信号保真度通过计算突触前和突触后锋电位序列之间的互相关（cross-correlation）来量化。脉冲时序依赖可塑性的建模考虑了轴突CV对突触前锋电位到达时间的影响。
三、研究结果
3.1 建模活动依赖性髓鞘形成
研究人员采用此前用于研究有髓神经元网络动力学的理论框架，将轴突髓鞘化变化操作化为轴突CV的变化。模型量化了轴突CV变化如何调节汇聚到共同突触后神经元的突触前输入的时间整合，以及这种时间整合的调节如何反过来塑造突触后动力学。采用纯前馈网络结构，排除了突触后反馈的混杂影响。
3.2 轴突CV调节突触前相关性
通过数学推导，研究人员发现突触后膜电位的方差随轴突CV增加而增大。当突触前输入存在相关性时，方差与突触前输入数量的平方成比例缩放，并依赖于CV。有效配对相关性概念被引入以综合考虑轴突长度和CV的影响。
3.3 轴突CV影响突触后兴奋性和放电率
利用高斯近似的突触后膜电位概率密度函数，研究人员计算了平均突触后放电率。结果显示，轴突CV的增加通过增大突触后膜电位方差来提高突触后兴奋性和自发放电率，而CV降低则产生相反效应。
3.4 轴突CV影响兴奋-抑制平衡
通过分别考虑兴奋性和抑制性突触前输入及其各自的CV，研究人员推导出基于方差的兴奋-抑制比（E/I ratio）。当兴奋性和抑制性轴突的CV相等时达到平衡；当两者不等时，E/I比率向兴奋或抑制方向倾斜，具体取决于ΔCV的正负。
3.5 轴突CV影响信号保真度
通过计算突触前和突触后锋电位序列的互相关，研究发现更快的CV促进突触前锋电位重合，从而放大突触前-突触后协方差，提高信号传递保真度。在纯兴奋性网络中，当突触前输入相关性非零时，互相关随CV增加而增大。
3.6 轴突CV与脉冲时序依赖可塑性（STDP）的联系
研究人员将STDP建模公式重新表述以包含轴突CV的影响。轴突CV的变化通过调整突触前锋电位相对于突触后锋电位的到达时间来驱动突触增强或抑制：CV增加使突触前锋电位提前（可能促进增强），CV降低则造成延迟（可能促进抑制）。
四、讨论与结论
本研究通过计算建模揭示了髓鞘变化如何通过调节轴突CV来影响神经元层面的动力学过程。研究结果表明，许多在神经元水平观察到的过程可能自然源于轴突CV影响突触前相关性的能力。具体而言，模型表明轴突CV可以在急性和慢性两种时间尺度上引起神经元活动、兴奋性、兴奋-抑制平衡以及上游活动模式忠实传递（即信号保真度）的变化。此外，结果为髓鞘变化与STDP调节之间的联系提供了可能的机制通路。急性阶段反映CV变化后、FRH启动前的突触后动力学；慢性阶段则反映FRH稳定后的突触后动力学。两者呈现出不同的动力学后果。
总体而言，模拟结果突出了髓鞘可塑性在细胞尺度的多种潜在动力学特征，对理解健康和病理脑状态下髓鞘可塑性的参与具有启示意义。轴突CV可能作为连接髓鞘可塑性与神经计算的关键机制纽带，为解释髓鞘失调的功能后果提供了原则性框架。
