[论文解读] A space-fractional cable equation for the propagation of action potentials in myelinated neurons

Drapaca, Corina S., Ozdemir, Sahin+1 more|arXiv (Cornell University)|May 17, 2020

Fractional Differential Equations Solutions被引用 3

一句话总结

本文提出一种基于Caputo分数阶导数的空间分数阶电缆方程，用于模拟髓鞘化轴突上动作电位传播中的非局部空间相互作用。通过将分数阶空间导数引入改进的Hodgkin-Huxley模型，研究表明：由细胞外空间和髓鞘中异常离子扩散所驱动的非局部性，导致膜电位影响范围更广，且动作电位幅度发生变化，特别是当 α = 0.65 时，其对实验观测的拟合效果优于经典局部模型（α = 1）。

ABSTRACT

Myelinated neurons are characterized by the presence of myelin, a multilaminated wrapping around the axons formed by specialized neuroglial cells. Myelin acts as an electrical insulator and therefore, in myelinated neurons, the action potentials do not propagate within the axons but happen only at the nodes of Ranvier which are gaps in the axonal myelination. Recent advancements in brain science have shown that the shapes, timings, and propagation speeds of these so-called saltatory action potentials are controlled by various biochemical interactions among neurons, glial cells, and the extracellular space. Given the complexity of brain's structure and processes, the work hypothesis made in this paper is that non-local effects are involved in the optimal propagation of action potentials. A space-fractional cable equation for the action potentials propagation in myelinated neurons is proposed that involves spatial derivatives of fractional order. The effects of non-locality on the distribution of the membrane potential are investigated using numerical simulations.

研究动机与目标

• 研究非局部空间效应在髓鞘化神经元动作电位传播中的作用。
• 模拟细胞外空间和髓鞘中异常离子扩散对膜电位动力学的影响。
• 建立并数值求解一种扩展经典电缆方程的时空分数阶电缆方程，引入分数阶空间导数。
• 比较分数阶模型（α = 0.65）与经典模型（α = 1）下的时空膜电位分布。
• 探索双向传播特性，以及变化的节点耦合参数（p, q）对非局部动力学的影响。

提出的方法

• 采用阶数 α ∈ (0,1) 的Caputo型分数阶导数，构建空间分数阶电缆方程，以表征非局部空间相互作用。
• 将Hodgkin-Huxley模型改进以包含非局部效应，通过分数阶电缆方程在漏电节间区修改离子通道动力学。
• 应用数值方法求解含分数阶空间导数的偏微分方程组。
• 模拟不同轴突直径和节间长度下的轴突单向及双向动作电位传播。
• 利用参数 α 控制非局部性的程度，其中 α = 0.65 表示强远距离相互作用。
• 在固定时间和位置分析膜电位分布，以比较分数阶与经典模型。

实验结果

研究问题

• RQ1 通过分数阶导数建模的空间非局部性，如何影响髓鞘化轴突中膜电位的分布？
• RQ2 轴突直径如何影响分数阶模型中非局部膜电位传播的范围与形态？
• RQ3 分数阶参数 α 如何影响动作电位诱导的电位变化的幅度与衰减速率？
• RQ4 在双向传播中，不同的耦合参数（p, q）如何影响非局部条件下的时空动力学？
• RQ5 与经典局部模型相比，分数阶模型在多大程度上能更准确地再现实验中观察到的跳跃式传导现象？

主要发现

• 当 α = 0.65 时，膜电位在节间区的更广区域保持非零，表明存在长距离非局部效应；而当 α = 1 时，电位在远离节点处迅速衰减至零。
• 随着轴突直径增大，分数阶模型（α = 0.65）中的非局部效应（远距离电位影响）更加显著，电位在节间区的大部分区域保持非零。
• 在固定位置 x = 700 μm 处，无论 α = 0.65 还是 α = 1，膜电位振荡幅度均随轴突直径增大而增加，且 α = 0.65 时幅度略高。
• 当节间长度较短时，α = 0.65 与 α = 1 的电位分布差异减小，表明非局部性在较长轴突节段中最为显著。
• 在双向传播中，降低耦合参数 p（从 1 降至 0.25）会减少固定位置处膜电位振荡的幅度，且该效应在分数阶模型中更为明显。
• 在 t = 0.1 ms 时，电压分布形状在不同 p 和 q 值下保持相似，但当 p 降低时，振幅系统性减小，表明在非局部框架中有效耦合减弱。