计算机模型预测缺血和扩散性抑制下的差异性树突易损性

《PLOS Computational Biology》：Computer models predict differential dendritic vulnerability with ischemia and spreading depression

　　

缺血，无论是突发性的还是慢性的，都会限制ATP（三磷酸腺苷）的生产并破坏依赖ATP的稳态机制，导致细胞内和细胞外离子浓度的改变。神经元ATP不足会触发K⁺（钾离子）释放，而细胞外K⁺浓度升高会使神经元去极化，导致更多K⁺释放；这种正反馈现象被称为扩散性去极化（SD，spreading depolarization）。当去极化效应足够强烈时，细胞进入去极化阻滞状态，即扩散性抑制（spreading depression）。过量的细胞外K⁺增加了Na⁺-K⁺泵的能量需求，形成了能量需求增加与能量供应减少的病理性汇聚。由此产生的变化在脑功能的亚细胞、细胞和网络尺度上均有深远影响。研究人员假设缺血或SD稳态失败的后果在亚细胞尺度上存在差异，且海马CA1锥体神经元的不同树突区域之间会有所不同。为了评估形态学与离子浓度之间的相互作用，研究人员采用了一种包含神经元形态学、泵、交换体、电压门控和Ca²⁺敏感离子通道的机制性模拟。在两种情况下，钙积累在基底树突中最为显著，提示这些树突将表现出最明显的兴奋性毒性效应。相比之下，在缺血（而非SD）条件下，远端顶树突暴露于更高的细胞内氯离子浓度，这可能导致树突珠串样改变。

论文解读

缺血性脑损伤是导致死亡和残疾的主要原因之一，其病理生理机制涉及能量代谢障碍、离子稳态失衡以及随后的细胞损伤级联反应。在缺血条件下，氧供应不足限制了ATP的生成，而ATP是维持神经元离子稳态的关键能量来源。Na⁺-K⁺泵作为消耗ATP最多的离子转运体，其功能受损会导致细胞内Na⁺积累和细胞外K⁺升高。细胞外K⁺浓度的升高可触发一种称为扩散性去极化（SD）的正反馈现象，即去极化导致更多K⁺释放，进一步加剧去极化。当去极化足够强烈时，神经元进入去极化阻滞状态，即扩散性抑制。缺血与SD之间存在病理性协同作用：SD会加重局部缺血，而缺血又阻碍SD的恢复。在临床上，SD与创伤性脑损伤、癫痫、肿瘤和缺血性卒中等多种神经系统疾病相关，其频率和持续时间可作为组织代谢状态的生物标志物。

然而，既往研究对缺血和SD在亚细胞水平上如何影响不同树突区域仍缺乏深入理解。神经元具有复杂的形态学结构，不同树突区域（如基底树突和顶树突）在离子通道分布、突触密度以及局部的氧和能量供应方面存在显著差异，这些差异可能导致缺血和SD条件下不同树突区域表现出不同的脆弱性。为此，研究人员构建了一个包含详细形态学信息的海马CA1锥体神经元计算模型，并结合一个简化形态的星形胶质细胞模型和周围细胞外空间，系统模拟了缺血和SD条件下的离子稳态崩溃过程，旨在阐明亚细胞水平上的差异性损伤机制。

该研究采用的主要关键技术方法包括：基于NEURON模拟器的反应-扩散模块构建多尺度三维模型，其中包含507个神经元区室，整合了离子通道、泵、交换体和钙信号通路。能量代谢模块显式模拟了O₂和ATP的浓度动态，Na⁺-K⁺泵、SERCA（肌浆网Ca²⁺-ATPase）和PMCA（质膜Ca²⁺-ATPase）均采用米氏动力学描述其对ATP的依赖性。谷氨酸稳态模型为每个区室配备了一个“平均突触间隙”，以模拟SD期间过量的谷氨酸能驱动。此外，模型还包含了细胞外空间内多种离子和分子的扩散模拟，采用Delaunay三角剖分限定边界，细胞外有效自由体积缩减5倍以近似多神经元共享的生理环境。研究还使用了来自耶鲁大学病理学系的三例62-75岁人脑皮层亚急性梗死组织标本进行神经病理学评估。模型参数通过BluePyOpt工具进行优化，进行了超过15万次模拟，以匹配12项电生理特征并维持离子稳态。

研究结果部分，首先在人类皮层梗死组织中观察到，神经丝免疫组化染色显示在梗死核心和半暗带均存在神经突珠串样变（beading），而邻近未受累皮层则无此现象，验证了该病理改变的临床相关性。进一步模拟发现，在正常氧条件下，神经元对500 ms的1.0 nA电流注入产生动作电位爆发后能迅速恢复静息膜电位。当细胞外K⁺升高至10 mM时，神经元在毫秒时间尺度内快速去极化并进入持续的SD样状态；而缺氧条件下，Na⁺-K⁺泵电流在2.95秒内下降至初始值的1%，导致延迟数十秒的去极化，且去极化幅度更大，伴随细胞外K⁺升高和细胞内Na⁺的持续替换。

离子浓度分布分析显示，在K⁺升高和缺氧两种条件下，细胞内Ca²⁺浓度均在基底树突中达到峰值，且浓度从基础的60 nM上升至mM范围。基底树突中更高的Ca²⁺积累归因于NMDAR（N-甲基-D-天冬氨酸受体）介导的Ca²⁺内流与Na⁺-Ca²⁺交换清除之间的失衡。紧密排列的基底树突竞争有限的O₂，导致局部ATP耗竭更为严重，残余Na⁺-K⁺泵活性降低，进一步限制了Na⁺-Ca²⁺交换介导的Ca²⁺清除。与此同时，ATP限制的Ca²⁺泵和SERCA活性也下降，共同导致基底树突Ca²⁺的大量积累。在Cl^-分布方面，缺氧条件下细胞内Cl^-浓度随距胞体距离增加而线性升高（平均1.61 mM/mm），在远端顶树突中达到最大值，这主要归因于KCC2（K⁺-Cl^-共转运体2）的反向转运。KCC2在缺氧条件下的反向转运是由于K⁺和Cl^-的跨膜梯度崩溃所致，而顶树突中较低的Na⁺-K⁺泵残余活性使其更易发生Cl^-积累。相反，在K⁺升高条件下，Cl^-浓度在基底树突中更高，且随距胞体距离增加。

在氧依赖性分析中，膜电位表现出全或无的转变特征：当O₂浓度≥0.03 mM时，神经元维持静息膜电位，而当O₂浓度≤0.02 mM时，神经元发生延迟去极化，去极化时间随O₂浓度降低而缩短（从9.5秒到28.3秒）。在维持膜电位的中间氧水平下（0.03 mM），虽然ATP显著下降且部分丧失离子梯度，但未产生大幅度的Ca²⁺和Cl^-病理性积累。此外，缺氧条件下Na⁺-K⁺泵活性在胞体和基底树突中显著下降（10%残余量分别出现在2.39秒和2.07秒），而顶树突在60秒时仍保持62.4%的泵活性，反映了形态学差异对能量竞争的显著影响。

在讨论部分，研究人员指出模型预测了缺氧性稳态丧失的三个阶段：（1）细胞内Na⁺替换K⁺而膜电位变化较小；（2）Na⁺通道开放导致快速去极化；（3）持续性离子浓度崩溃。该模型强调了基底树突的高脆弱性源于多种因素的协同作用：谷氨酸清除减少导致NMDAR活化介导的Ca²⁺内流增加、局部K⁺和Na⁺浓度的升高以及O₂的竞争性匮乏。与此相反，顶树突在缺血条件下由于KCC2失衡而出现更显著的Cl^-积累，这与远端树突先出现珠串样变、随后向近端扩展的临床观察相吻合。研究人员还讨论了模型的局限性，包括未纳入电中性约束和体积动力学、未模拟糖酵解和乳酸穿梭等无氧ATP生成途径，以及未考虑神经血管耦联等因素。尽管如此，该模型是迄今为止在亚细胞水平上对离子失衡和缺血综合效应进行的最为生物物理详尽的模拟，为理解不同病理类型下的树突损伤模式和开发靶向性神经保护策略提供了重要的理论框架。

研究结论表明，缺血和SD虽均导致离子稳态崩溃，但产生了不同的亚细胞离子分布模式：基底树突对Ca²⁺相关的兴奋性毒性最为敏感，而远端顶树突在缺血条件下更易发生Cl^-相关的珠串样改变。这些发现提示，不同树突区域具有不同的病理易损性特征，可为预测组织损伤模式、分层评估预后以及指导针对性治疗提供依据。该成果发表于《PLOS Computational Biology》。