# 基于电子显微镜重建星形胶质细胞分支几何的钙扩散与钙-钙调蛋白反应-扩散过程多平台模拟系统分析

《PLOS Computational Biology》：Modeling reaction–diffusion processes in simplified to complex astroglial branchlet geometries using different simulation tools

　　

神经科学中空间生化过程的计算建模正日益重要，因为实验与成像技术的进步不断提供关于细胞组织结构的精细信息。反应-扩散模型（reaction–diffusion models）常用于研究神经元细胞内离子与分子的扩散及相互作用，并已有多个开放获取的模拟平台被开发以支持此类分析。然而，目前仍缺乏对常用模拟工具（包括MCell、STEPS、NeuroRD和NEURON）的全面分析，尤其缺乏针对真实胶质细胞几何形态的分析，因为这些工具最初是为神经元开发的。在此，研究人员通过在一系列空间域中实现等效的钙扩散与钙-钙调蛋白（calmodulin，CaM）反应-扩散模型，对上述四种模拟环境进行了系统分析。所涉空间域包括简化几何结构（圆柱体和长方体）以及源自电子显微镜（electron microscopy，EM）数据的生物学真实星形胶质细胞形态，涵盖小鼠Bergmann胶质细胞以及小鼠和人皮层星形胶质细胞。研究人员采用匹配的动力学参数和扩散常数，在可比条件下考察各平台的模拟输出与形态学差异。结果发现，在所考虑的系统中，所有工具均产生总体一致的结果，表明每一种工具都能可靠地刻画简化与复杂几何结构中的扩散及反应-扩散动力学。同时，研究人员发现，人星形胶质细胞和小鼠Bergmann胶质细胞分支的结构异质性会产生不同的钙信号，表明精细尺度形态对星形胶质细胞动力学具有关键塑造作用，且必须被模拟框架准确刻画。此外，各工具在几何表示、实现策略、计算层面和易用性方面存在显著差异。通过揭示细胞内生化过程建模中几何保真度、灵活性与计算成本之间的关键权衡，这些发现为选择合适的模拟工具提供了实用指导。更广泛地讲，该结果支持在形态真实的星形胶质细胞环境中应用反应-扩散建模，并有助于实现突触及突触周位点局部信号事件的机制研究。

**一、研究背景与问题**

亚细胞水平时空生化信号解析的计算资源与实验手段不断进步，使脑系统的反应-扩散建模日益受到关注，用以补充和拓展传统基于反应的模型。胶质细胞（如星形胶质细胞）形态复杂，其时空组织的信号需要比当前领域最常用的单房室模型更精细的模型。

现有开放获取模拟工具可分为三类：基于布朗动力学（Brownian dynamics）的粒子型微观模型，如MCell与Smoldyn；基于反应-扩散主方程（reaction–diffusion master equation）的体素型介观模型，如STEPS、NeuroRD与MesoRD；基于偏微分方程（partial differential equations）的房室建模工具，如NEURON、ASTRO与GENESIS/Chemesis；此外还有Virtual Cell、MOOSE等混合方法工具。既往研究或对部分工具进行数值基准测试，或作定性分析，但缺乏在一致的模型与真实细胞形态下对MCell、STEPS、NeuroRD和NEURON rxd模块的联合系统评估，尤其在生物学真实几何中尚未有研究在可比条件下共同应用这四种工具。

体积EM技术可在纳米分辨率下重建细胞形态，MICrONS小鼠皮层数据集、人皮层重建与小鼠小脑数据集等大规模数据通过Neuroglancer、BossDB等平台日益可及；Blender、MeshLab、GAMer 2、NeuroMorphoVis、Ultraliser、Alpha_Mesh_Swc、CellRemorph等工具支持形态的可视化、变换与预处理。然而，将EM衍生几何整合进反应-扩散模拟仍很有限，多数模型依赖理想化几何或聚焦单一星形胶质细胞微域，跨脑区与跨物种的形态多样性很少在空间反应-扩散模拟中探讨，这一缺失在皮层星形胶质细胞与小脑Bergmann胶质细胞之间、以及啮齿类与人星形胶质细胞之间存在显著结构差异的背景下尤为突出。因此，研究人员需要系统分析不同模拟工具如何处理此类形态，以充分利用日益增长的高分辨率EM数据。

星形胶质细胞发出细突起包绕并部分包裹突触，构成三联突触（tripartite synapse）的星形胶质成分，支持高度局域化的信号事件。超微结构分析显示，这些分支为极细、具有高表面积体积比（surface-to-volume ratio）的结构，细星形胶质突起报道值约为19–33 μm^-1，小脑Bergmann胶质突起也存在类似高比值微域。突触周星形胶质突起可通过代谢型谷氨酸受体（metabotropic glutamate receptors）等神经递质激活受体检测突触活动，突触诱发信号可触发突起内局域钙升高。纳米尺度成像进一步揭示分支内显著的局部异质性，包括由更细的柄连接的球状腔室以及收缩、扩张和核心等结构基序。这些几何特征可促进生化区室化并影响胞内钙信号在相邻腔室间的传播。尽管近期已有多个人皮层组织EM数据集可用，纳米尺度星形胶质细胞分支形态的定量分析仍有限。这些观察提示，分支口径、体积与连接性的局部变化可能强烈塑造胞内扩散与反应动力学，凸显真实分支几何在星形胶质细胞钙信号及下游钙依赖生化过程模型中的重要性。

据此，研究人员在多个空间域（圆柱体、长方体等简化几何与EM重建的真实星形胶质细胞形态）中实现等效的钙扩散与钙-钙调蛋白反应-扩散模型，并在各平台支持范围内开展模拟。一维情形聚焦NEURON与NeuroRD，三维模型主要考察MCell、STEPS与NEURON，并在长方体中补充NeuroRD模拟。

**二、主要关键技术方法**

第一，选用四种开放获取反应-扩散模拟工具：粒子型蒙特卡罗（Monte Carlo）工具MCell（经CellBlender界面）、基于四面体网格的随机与确定性求解器STEPS（Tetexact、TetODE）、基于异步tau-leaping的NeuroRD以及NEURON的rxd模块，并构建统一可比模型。第二，样本队列来源方面，形态数据取自小鼠小脑皮层序列EM数据集（Bergmann胶质细胞）、MICrONS小鼠视觉皮层数据集与人颞叶皮层petavoxel EM重建（人星形胶质细胞），共8种形态；使用CellRemorph、Neuroglancer与CloudVolume提取与筛选分支。第三，统一的网格制备流程：构建水密多边形表面网格，经MultiCompMesher生成四面体体积网格供STEPS使用，再回提表面网格供MCell；对NEURON则将表面网格体素化、骨架化为SWC形态并作直径缩放以匹配原体积。第四，建立测地距离（geodesic distance，基于Dijkstra算法）等体积分箱框架，将不同离散化下的三维钙分布映射为可跨平台比较的一维剖面。第五，采用RMSD、NRMSD与皮尔逊相关系数（Pearson Correlation Coefficient）定量比较各工具结果，随机模拟取20次运行的均值与标准差。

**三、研究结果**

（一）简化与复杂星形胶质细胞几何中的钙动力学：钙扩散模拟。1D简化几何中，NeuroRD与NEURON在25个亚体积下结果高度相似，符合大亚体积且等截面积几何的预期，尽管前者随机、后者确定。3D简化几何中，四种工具在模拟末段表现一致，同一工具的圆柱与长方体结果相近；但工具间存在差异，圆柱结果常比长方体更相似；NeuroRD亚体积最少，可能贡献更快的扩散，确定性STEPS与随机STEPS单次实现及其他工具结果不同。真实细胞形态的3D扩散中，小鼠Bergmann胶质分支形态2里MCell与STEPS高度一致，确定性STEPS落在随机模拟变异内；NEURON呈现更快的扩散，尤以后三个时间点明显，研究人员认为可能源于其SWC树重建与立方体素离散对薄膜结构与环状区域的平滑或额外连通；但NRMSD均低于5%、相关系数显示各工具高度一致。小鼠Bergmann胶质分支形态1结果同样一致，但更快达到稳态。小鼠皮层星形胶质分支形态1中，NEURON结果更多落在随机模拟标准差内，因其环状少、更呈树状，与NEURON几何表示更兼容；确定性STEPS仍更接近随机STEPS而非确定性NEURON。人皮层星形胶质分支形态1中MCell与STEPS一致，而NEURON因体素离散平滑细窄缩颈而显示更快扩散；人分支形态2与皮层形态1行为类似，但扩散发生更快。形态学可视化显示钙扩散受分支几何约束、沿形态而非均匀分布扩散，工具间差异在该时间尺度上较细微。

（二）反应-扩散模拟。1D简化几何中NEURON与NeuroRD高度一致，NeuroRD的随机性可见。3D简化几何中四工具结果高度一致，STEPS与MCell中Ca₄CaM峰值约90–100，NeuroRD略低于90、NEURON略低于80；由于STEPS确定性TetODE在原精细网格上运行时间不可行，研究人员在保持几何特征的粗化网格上验证了其确定性与随机求解结果一致。真实形态中，两种Bergmann胶质分支下工具表现高度一致，NEURON在形态1中钙与钙调蛋白平台略高、Ca₄CaM基线略低；20次运行的变异相对较小，反映参与反应分子数多；NRMSD低于8%。小鼠皮层星形胶质分支形态1中，NEURON钙与钙调蛋白基线略高、Ca₄CaM峰值约50，而MCell与STEPS约80–90，Ca₄CaM的NRMSD超过7%；对简化后的该形态，STEPS确定性与随机结果一致。人皮层星形胶质分支形态1中NEURON的Ca₄CaM峰值约100，低于MCell与STEPS的约140；形态2中NEURON也能跟随MCell与STEPS曲线。形态学可视化显示三工具反应-扩散行为一致，但MCell在Bergmann胶质分支形态2中Ca₄CaM分布略宽；人分支形态2中所有工具形成的Ca₄CaM较少。

（三）反应-扩散框架的适用性。MCell的CellBlender界面直观、可处理任意水密三角网格，但文档陈旧（描述MCell3与Blender 2.79）；STEPS经Python界面使用顺畅、可直接在四面体网格上计算，但需网格预处理，且确定性求解计算代价高；NeuroRD不支持一般1D模拟需借2D强制实现，复杂形态实现困难，故未实现真实形态；NEURON经Python与rxd使用友好，但需形态转换与外部可视化，SWC重建后体积偏差多数约1–6%、小鼠皮层形态1约17%；1D电缆形式不能表示环状拓扑，而3D rxd可借共享体素恢复部分环状特征。简化与真实形态的表面积体积比分别约9–10 μm^-1与约9–20 μm^-1。

**四、讨论与结论（译文）**

本研究分析了四种主要用于神经元反应-扩散系统建模的先进开放获取工具，并考察星形胶质细胞分支形态对反应-扩散建模的意义。研究人员发现，人星形胶质细胞与小鼠Bergmann胶质分支的结构异质性会产生定性不同的钙信号状态，表明精细尺度形态关键地塑造星形胶质细胞动力学；钙、钙调蛋白与Ca₄CaM等关键模型变量依局部几何呈现显著不同的信号行为，凸显星形胶质细胞钙动力学对窄通道、复杂分支模式与环状结构等微观结构特征的强敏感性，也对模拟框架准确解析形态依赖的扩散与反应提出高要求。

选择这四种工具是因为据研究人员所知它们尚未在单一研究中被联合考察，尤其缺乏真实星形胶质细胞形态背景。目标是评估其在未来复杂星形胶质分支（包括小鼠与人皮层星形胶质分支及Bergmann胶质分支）钙动力学研究中的适用性，而非产生完全优化或工具特定的基准实现。结果显示，钙扩散与钙-钙调蛋白反应-扩散在不同工具间高度一致，尽管其确定性/随机算法不同。既往研究表明，低拷贝数下随机方法可刻画离散分子相互作用引起的涨落，但其后果依系统而定；本模型缺乏放大的可兴奋或正反馈机制，故预期随机性主要表现为运行间变异而非定性的不同信号状态，与观察一致。由于1D模型未含表面积依赖过程，长方体与圆柱表面积体积比差异未导致动力学差异。某些情形出现差异，尤其是NeuroRD长方体3D扩散与NEURON在更复杂星形胶质形态中的结果；反应-扩散中早期可见Ca₄CaM瞬态空间异质性，最初集中在注射点附近并在部分形态中局限于被窄分支分隔的特定子区域，反映几何收缩下的扩散限制，类似树突棘细颈的暂时限制作用；随模拟推进，Ca₄CaM在各工具中趋于均匀分布。未来可纳入多钙源与内质网等胞内区室。四种工具均适于此模型，但在易用性、形态格式（MCell多边形表面、STEPS体积网格、NeuroRD立方网格、NEURON骨架）、是否支持EM真实形态、划分亚体积方式、输入形式与运行时间（MCell约10分钟，STEPS数小时；NEURON由0.17秒至数小时，NeuroRD由0.014秒至约2分钟）等方面有明显差异；STEPS的效率差异应结合模型与求解器理解，其优势预期在稀疏反应、强区室化或通道介导钙流情形。研究还提及Smoldyn、MesoRD、Virtual Cell、PyRID、NEST与ARACHNE、Brian2、Uncertainpy及物理信息神经网络（physics-informed neural networks，PINNs）等相关工具与方向。

结论：本研究分析了MCell、STEPS、NeuroRD与NEURON rxd四种广泛使用模拟工具在扩散与反应-扩散方案中的易用性、几何处理、实现投入、灵活性与跨模拟器一致性，并考察其如何表示星形胶质分支的复杂亚细胞几何。分析凸显各平台的优点与局限，为其在未来网络级模拟协同仿真框架中整合反应-扩散过程提供洞见。结果表明，无论简化几何还是真实星形胶质分支形态，工具间均高度一致；由EM数据重建的部分星形胶质几何使系统评估形态复杂性如何影响跨平台模拟结果成为可能。研究强调形态对时空依赖的钙扩散与钙-钙调蛋白反应-扩散系统的影响；进一步研究可借助更复杂生化模型更深入揭示形态的作用。总体而言，这些结果阐明了当前工具在复杂细胞几何中反应-扩散建模的能力与局限，并为其未来使用提供指导。