[论文解读] A Multidomain Model for Ionic Electrodiffusion and Osmosis with an Application to Cortical Spreading Depression

Yoichiro Mori|arXiv (Cornell University)|Oct 30, 2014

Electrochemical Analysis and Applications参考文献 36被引用 8

一句话总结

本文提出了一种热力学一致的多域模型，用于描述生物组织中的离子电扩散和渗透性水流动，该模型被表述为一组守恒自由能的偏微分方程。将该模型应用于皮层扩散性抑制（cortical spreading depression）时，成功计算出实验观测到的细胞外电位负直流偏移，并表明间隙连接耦合不太可能是驱动SD传播速度的主要因素。

ABSTRACT

Ionic electrodiffusion and osmotic water flow are central processes in many physiological systems. We formulate a system of partial differential equations that governs ion movement and water flow in biological tissue. A salient feature of this model is that it satisfies a free energy identity, ensuring the thermodynamic consistency of the model. A numerical scheme is developed for the model in one spatial dimension and is applied to a model of cortical spreading depression, a propagating breakdown of ionic and cell volume homeostasis in the brain.

研究动机与目标

• 开发一种用于多相生物组织中离子电扩散和渗透性水流动的热力学一致模型。
• 通过引入离子浓度动力学和水流动，扩展双域框架，确保通过自由能恒等式实现能量守恒。
• 将该模型应用于皮层扩散性抑制（cortical spreading depression），这是一种与偏头痛和中风相关的病理脑现象。
• 将细胞外电位偏移（负直流偏移）作为模型的自然输出，这是SD的关键实验特征。
• 研究间隙连接耦合和细胞外氯离子浓度在调节SD传播速度中的作用。

提出的方法

• 为多个相互穿插的组织区室中的离子传输和水流动制定一组时间依赖的偏微分方程系统。
• 使用带有源项的能斯特-普朗克方程描述跨膜离子通量，使用水质量守恒方程描述水平衡。
• 通过施加自由能恒等式确保热力学一致性，借鉴昂萨格互惠关系及软物质物理领域的前期工作。
• 在一维空间中应用数值格式求解该模型，初始条件来自稳态预处理模拟。
• 将细胞外电位作为系统的自然输出计算，避免对电中性或电容关系的假设。
• 采用两区室（细胞内和细胞外）框架，显式建模离子通量、膜电位和体积分数动力学。

实验结果

研究问题

• RQ1 热力学一致的多域模型能否准确描述脑组织中的离子电扩散和渗透性水流动？
• RQ2 间隙连接耦合在决定皮层扩散性抑制传播速度中的作用是什么？
• RQ3 细胞外氯离子浓度如何影响皮层扩散性抑制的动力学？
• RQ4 该模型能否定量再现SD实验中观测到的细胞外电位负直流偏移？
• RQ5 快速电流传导与缓慢的离子浓度梯度如何在统一且能量守恒的框架下相互作用？

主要发现

• 该模型成功将细胞外电位的负直流偏移作为自然输出计算出来，这是皮层扩散性抑制的关键实验特征，据作者所知，这是首次实现此类计算。
• 预处理初始状态在离子浓度变化上的相对误差小于10−12时收敛至稳态，表明数值稳定性和一致性。
• 间隙连接耦合被发现不太可能是SD传播速度的主要驱动因素，因为在该模型的参数化下，其对速度的影响微乎其微。
• 细胞外氯离子浓度被证明会影响SD传播速度，表明其在病理过程中具有调节作用。
• 该模型的热力学一致性确保了电扩散和渗透效应以物理上合理的方式耦合，从而能够准确模拟离子和体积稳态的破坏。
• 在无间隙连接耦合的情况下，可基于离子通量与电导率之间的平衡，推导出细胞外电位的简化表达式，但该推导仅在严格条件下成立（一维、无通量边界）。