GENESIS细胞结构建模：从SWC数据到多房室神经元仿真

GENESIS系列写了好几篇，之前一直在聊通道动力学、突触机制这些偏“功能”的东西。这次第七篇，我打算回头把基础打牢，专门讲细胞结构建模。

为什么要专门写这个？因为我发现很多刚开始玩GENESIS的人，喜欢一上来就堆channel、调突触参数，结果仿真出来的波形怎么都不对。问题大概率不是出在通道上，而是细胞结构本身就没建对。在电生理仿真里，细胞结构不是背景道具，它直接决定了电流怎么在空间上传播、树突信号怎么整合、胞体动作电位怎么触发。结构建歪了，后面所有分析都是在沙滩上盖楼。

这篇文章我会从GENESIS里“结构”到底指什么讲起，然后带你从形态学数据开始，一步步算出建模需要的几何参数，再用GENESIS脚本把一个带树突分支的神经元骨架搭出来，最后聊一聊怎么验证结构对不对、怎么把结构图输出成PDF留存。不管你是刚接触计算神经科学的学生，还是已经在跑网络模型但想回头抠单细胞细节的研究者，这篇都能给你一套可以直接复现的流程。

1. 细胞结构建模的核心思路与建模对象

1.1 为什么单房室模型不够用

很多人第一次接触神经元仿真，用的都是单房室模型：把整个神经元当成一个等电位球，所有的离子通道都怼在同一个点上，膜电位处处相同。这个模型对理解离子通道的基本动力学很有用，但一到真实场景就露馅了。

真实神经元的结构非常不均匀：胞体大而圆，树突细长且高度分枝，轴突又细又长。这些不同的结构段，膜面积不一样，轴向电阻不一样，受到突触输入后的局部反应也不一样。如果全部压成一个点，你就完全看不到树突远端输入在向胞体传递过程中的衰减和延迟，也看不到多个树突分支对输入的整合效应。

我在实际做仿真时最深刻的体会是：加入空间结构之后，很多原本需要靠精细通道参数才能调出来的现象，其实单纯靠结构就能自然出现。比如远端树突的EPSP传到胞体时幅值会变小、时程会被拉宽，这不需要任何复杂机制，几何结构本身就决定了。所以，单房室模型适合学原理，但要做真正能和电生理实验对标的研究，多房室建模是绕不开的一步。

1.2 GENESIS眼中的“细胞结构”是什么

GENESIS里对细胞结构的基本抽象特别直接：一切结构上的单元都是一个“compartment”（房室），你可以把它理解成一段圆柱体。胞体是一段比较短粗的圆柱，树突是一段段细长的圆柱，轴突也是一段段圆柱。多个圆柱通过拓扑连接关系串在一起，就组成了一个完整的神经元形态。

每个compartment内部有两个核心属性：一个是它的几何尺寸，也就是length和diameter，这决定了它代表多大的膜面积、内部有多少胞质电阻；另一个是它的生物物理属性，包括膜电容Cm、膜电阻Rm、轴向电阻Ra、静息电位Em等。GENESIS的求解器在做仿真时，会把每一个compartment当成一个独立的电路节点，然后根据compartment之间的连接关系建立方程组，统一求解。

理解了这个抽象方式，你就能明白：在GENESIS里“建结构”，本质上是在做两件事——确定每个compartment放在哪、长什么样（几何），确定compartment之间怎么连（拓扑）。这两件事做扎实了，后面的生物物理参数只有填进去的份。

1.3 结构建模在完整仿真流程中的位置

一个完整的GENESIS神经元模型开发流程，通常分四步：第一步是拿到或重建形态学数据；第二步是把这个形态数据离散成一个个compartment并定义连接；第三步是在不同compartment上添加离子通道、突触等机制；第四步是加刺激、跑仿真、做分析。

结构建模就是第二步，它卡在“数据”和“机制”之间。这一步做得越精细，后续加通道时的位置映射就越准确。比如某些钾通道只在胞体和近端树突表达，某些钙通道只在远端树突表达，你只有在结构建清楚的前提下，才能精确指定这些通道挂在哪个compartment上。

另外提醒一句，GENESIS的结构建模不是一个一次性的工作。你在跑仿真过程中发现某个分支电阻太大导致信号传不过去，或者某个房室太小导致数值不稳定，都要回来改结构。所以，结构建模不是“建完就算完”，它是整个仿真周期里需要反复迭代的底盘部分。

2. 建模前必须搞定的几何数据与参数计算

2.1 形态学数据从哪来

做细胞结构建模，第一步当然是拿到细胞形态。这里有三条比较常用的路子。

第一条路，去公开数据库下现成的形态学数据。最常用的是Neuromorpho.org，这个库收集了海量已发表的神经元重建数据，支持按物种、脑区、细胞类型筛选，下载格式通常包含SWC。SWC是目前最通用的神经元形态格式，每一行代表一个节点，记录了节点编号、类型、三维坐标、半径，以及与父节点的连接关系。拿到SWC文件，基本上就拿到了建模的“图纸”。

第二条路，从文献里抠参数。如果你研究的细胞类型没有人发过完整的重建形态，那就退而求其次，找文献中报道的胞体直径、树突长度、分支数量等典型值，自己手工搭一个简化结构。这个方法精度低一些，但对于很多机制研究来说完全够用。

第三条路，自己从实验图像里重建。这需要你有共聚焦或双光子图像数据，用半自动追踪工具导出形态。这条路工作量最大，但最贴合你自己的实验数据。个人建议，非必要不走这条路，前期用公开数据跑通流程更重要。

2.2 坐标、长度单位的统一处理

拿到SWC文件后，千万别直接往GENESIS里塞，先做单位统一。SWC文件里的坐标通常以微米为单位，而GENESIS内部计算时用的是标准单位，如果你不换算直接填进去，可能出现长度差好几个数量级的离谱模型。

我自己的习惯是这样：先把所有几何量统一换算成米，再填进GENESIS字段。因为GENESIS里compartment的length字段单位就是米，diameter也是米。如果你从SWC里读到某段树突长度是150.5，数字旁边有注释说明单位是微米，那就得先乘1e-6再填。

这里有个坑特别容易踩：SWC文件的节点代表的是采样点，节点之间的连线才是一段真实的树突片段。所以你在建模时，不能把每个节点当成一个compartment直接建，而是要把相邻两个节点之间那段当成一个compartment。也就是说，compartment的length是相邻节点的欧氏距离，diameter取了两个节点半径的平均值或较小值。我自己一般取平均值，这样更平滑。

2.3 轴向电阻和膜参数的推导计算

结构建模不只是“画形状”，更重要的是给每一段结构赋予物理意义。这里面最有技术含量的是轴向电阻的计算。

轴向电阻描述的是电流在胞质内流动时遇到的阻力，它和两个因素有关：一是胞质的固有电阻率Ra，二是这一段圆柱的几何形状。计算公式是：

R_axial = Ra * L / (π * (D/2)^2)

其中Ra的单位是欧姆·厘米，L是房室长度，D是直径。注意L和D都要换算成厘米才能和Ra对上。很多新手在这里栽跟头，就是因为单位没换算一致，算出来的电阻差了10的8次方，仿真自然没法看。

举个例子，一段树突，长度80微米，直径1.5微米，Ra取100欧姆·厘米。换算之后，长度是0.008厘米，半径是0.000075厘米。带入公式得到：

R = 100 * 0.008 / (3.1416 * (0.000075)^2) ≈ 4.53e8 欧姆

看到这个数值不要慌，膜电阻动辄几百兆欧，轴向电阻在细长结构里本来就是个很大的量级。

膜参数里，Cm通常取1微法每平方厘米，然后用膜面积去换算成具体compartment的电容值。膜面积是2πrL，如果直接用GENESIS的默认单位体系，你只需要保证Cm、Rm填的是单位面积归一化后的值，GENESIS的求解器会结合几何尺寸自动折算。这也是为什么很多人建议用GENESIS自带字段而不自己手算Cm，避免二次换算出错。

3. GENESIS脚本实操：搭建一个多房室神经元骨架

3.1 最小区块：compartment的创建与参数写入

在GENESIS里创建一个房室，用的命令是create。最基础的写法是：

create compartment soma

这条命令会创建一个名为soma的compartment对象，放在当前默认路径下。之后你用setfield给它写入几何和生物物理参数，用addmsg把它和其他房室连起来。

setfield soma Cm 0.01 Rm 0.3 Ra 1.0 Em -0.07 setfield soma length 10e-6 diameter 20e-6

这里我故意没有直接用上面算出来的具体值，因为不同细胞类型的参数差异太大。但有一点必须弄清楚：GENESIS里面Cm、Rm、Ra这些字段的单位和普遍文献用的单位习惯不完全一样。Cm的单位是法拉，Rm的单位是欧姆·米平方，Ra的单位是欧姆·米。所以你在填数值之前，先看清自己手里的文献参数用的什么单位，做完换算再填。

我会在模型文件里加注释，把原始文献值和换算后的填值并列写清楚，这样过一个月回来看脚本，还能知道当初每个数是怎么来的。

3.2 分支拓扑的表达与连接规则

有了单独的compartment之后，关键问题来了：怎么把多个compartment串成分支结构？答案是addmsg命令，它用来建立compartment之间的信息传递关系。

GENESIS中最常用的连接消息有几种：AXIAL表示两个房室之间存在轴向电流通路，这是所有相邻房室之间必须有的；VMEQ表示把某个房室的膜电位作为输入传给本房室的机制对象，一般每个房室自己给自己发一条；如果是连接通道相关，还会用到CHANNEL等。

拿一个最简单的二分叉树突来举例：soma分出两个分支dend1和dend2，dend1末端又分出dend1a和dend1b。那么拓扑连接关系是这样的：

addmsg soma dend1 AXIAL addmsg dend1 soma AXIAL addmsg soma dend2 AXIAL addmsg dend2 soma AXIAL addmsg dend1 dend1a AXIAL addmsg dend1a dend1 AXIAL addmsg dend1 dend1b AXIAL addmsg dend1b dend1 AXIAL

注意这里是双向的，每个相邻pair要连两个方向的AXIAL，因为电流既可以从soma流向树突，也可以从树突流回soma。如果只连单边，仿真时电流就会像单向阀门一样，结果完全错误。

3.3 完整脚本示例：一个带两级分支的神经元

下面我给一个可以直接保存运行的完整脚本。这个脚本在GENESIS里构建了一个简化神经元：一个胞体、两根初级树突、其中一根树突继续分出两个次级分支。

// 细胞结构建模示例：soma + 2 primary dendrites + secondary branches create neutral /cell pushe /cell create compartment soma setfield soma Cm 0.01 Rm 0.3 Ra 1.0 Em -0.07 setfield soma length 12e-6 diameter 20e-6 create compartment dend1 setfield dend1 Cm 0.01 Rm 0.3 Ra 1.0 Em -0.07 setfield dend1 length 80e-6 diameter 2e-6 create compartment dend2 setfield dend2 Cm 0.01 Rm 0.3 Ra 1.0 Em -0.07 setfield dend2 length 100e-6 diameter 2e-6 create compartment dend1a setfield dend1a Cm 0.01 Rm 0.3 Ra 1.0 Em -0.07 setfield dend1a length 60e-6 diameter 1.5e-6 create compartment dend1b setfield dend1b Cm 0.01 Rm 0.3 Ra 1.0 Em -0.07 setfield dend1b length 70e-6 diameter 1.5e-6 addmsg soma dend1 AXIAL addmsg dend1 soma AXIAL addmsg soma dend2 AXIAL addmsg dend2 soma AXIAL addmsg dend1 dend1a AXIAL addmsg dend1a dend1 AXIAL addmsg dend1 dend1b AXIAL addmsg dend1b dend1 AXIAL addmsg soma soma VMEQ addmsg dend1 dend1 VMEQ addmsg dend2 dend2 VMEQ addmsg dend1a dend1a VMEQ addmsg dend1b dend1b VMEQ pope

这段脚本跑起来，你就在GENESIS里拥有了一个可仿真的空间非均匀神经元模型。注意这里我还加了每个compartment给自己的VMEQ消息，这一步是让机制对象能读取本房室的电压，很多刚入门的人容易漏掉。

3.4 结构建模里的几个隐蔽坑

第一个坑是重复创建compartment。如果你在一个脚本里多次执行create compartment soma，而且当前路径没有切走，GENESIS不会报错，而是会生成一个带后缀的新对象或者直接覆盖。这个问题在大型模型里特别隐蔽，因为报错信息不明显，结果就是仿真结果漂移。我的习惯是，在每个脚本开头先用delete /cell清理旧对象，或者切换到新的neutral路径下创建。

第二个坑是pushe和pope的配对问题。pushe是把路径切到某个object里面，pope是弹回来。如果你push了忘记pop，后面的setfield可能全部写到错误的对象里，排查起来非常痛苦。我在脚本里习惯在每个区块末尾检查一下路径层级。

第三个坑是直径的单位。很多形态学重建软件给出的直径单位并不统一，有的给微米，有的给毫米，甚至有直接给像素的。在填入diameter字段之前，一定要确认这个数字代表的是直径而不是半径。我在一次建模里就因为把半径当直径填了，导致所有树突的输入阻抗差了整整一倍，后来查了好久才发现是这里。

4. 结构建模的检查与可视化：从仿真验证到PDF导出

4.1 仿真第一课：结构和参数对不对，先看静息电位

结构搭好了，别急着加刺激，先跑一个没有任何输入的静息状态仿真。这一步的目的是验证全局的数值稳定性，以及结构是否处于合理的电学状态。

正常的做法是设置一个足够长的仿真时间，比如100毫秒，然后记录soma的电压。如果结构搭建正确，各compartment的静息电位最终都应该稳定在Em附近，不会有趋势性的漂移，更不会发散到离谱的值。

如果发现soma电压稳在了一个明显不是Em的值，或者曲线呈锯齿状震荡，基本可以确定结构或者参数有问题。常见的可能性有：AXIAL消息漏连、某个compartment的Rm填错量级、Cm单位不一致。这些都值得逐一排查，不要急着加电流刺激掩盖问题。

4.2 电流注入测试：结构对信号传导的影响

静息稳定之后，我会在胞体注入一个短促的方波电流，观察动作电位或者被动衰减的传播情况。这一步能够直观看出结构对信号传导的影响。

比如在soma注入一个400皮安、2毫秒的方波，你可能会关注以下几个问题：soma电压变化幅度多少？dend1远端的电压变化是不是明显变小、变慢？dend2因为更长，衰减是不是比dend1更严重？这些现象都是结构本身带来的，不需要任何通道机制参与就能看到。

我拿到一个新模型时，经常用这个方法来快速判断结构是否“说得通”。如果soma的输入阻抗异常大，或者远端信号几乎没有衰减，那我就要怀疑是不是某个compartment的几何尺寸填错了。

4.3 结构可视化：把模型“画”出来检查

脚本建出来的结构毕竟是抽象的，肉眼很难发现问题。所以一定要可视化检查。GENESIS本身带有简单的绘图工具，比如xgraph，可以实时画出某个compartment的电压曲线。但这只是数值可视化，不是结构形态可视化。

如果你想真正“看到”建出来的细胞长什么样，我常用的做法是把结构参数导出到文本文件，再用外部工具画。你可以用print命令把每个compartment的编号、type、坐标、直径、父节点信息输出成一个SWC文件，然后扔给支持SWC可视化的工具去看。

这个步骤特别值得做一次。因为你会发现，很多在脚本里感觉合理的结构，画出来之后可能长得很奇怪。比如有的分支角度不合理，有的compartment长度和直径比例严重失调，这些都会影响仿真精度。

4.4 打印PDF命令的使用方法详解

很多同行问我，GENESIS跑的图怎么导出成PDF放到论文里。这里我说一下我常用的操作流程。

GENESIS本身不直接生成PDF，它的图形输出更多是靠仿真完成后保存数据，再用外部绘图工具出图。我的标准做法是：先用GENESIS把仿真数据保存下来，再用gnuplot或Python重新绘图，最后导出为PDF。

在GENESIS脚本里，保存数据的命令是print配合重定向。示例：

str file = "sim_result.dat" openfile {file} w print {file} {soma.Vm} closefile {file}

或者更简单一些，用print -i soma.Vm -t 0 -v 0 -o sim_result.dat把指定时间范围内的电压波形直接输出到文件。

拿到数据之后，在gnuplot里画图并输出PDF的命令很直接：

set terminal pdf set output "soma_voltage.pdf" plot "sim_result.dat" using 1:2 with lines

这里有一个小细节：gnuplot的PDF终端生成的是矢量图，放进论文里缩放不糊。如果要在Linux命令行直接一条命令出图，可以写成：

gnuplot -e "set terminal pdf; set output 'soma_voltage.pdf'; plot 'sim_result.dat' using 1:2 with lines"

把它写成一个shell脚本，以后每次跑完仿真，数据一出来图就出来了，非常省事。

需要说明的是，有些老版本的GENESIS系统里打包了基于X11的绘图组件，可以直接在屏幕上截图存成位图，但画质远不如重新出图。我个人强烈建议走“数据落盘+外部绘图”这条路，既灵活又清晰。

5. 高频问题排查与实战经验速查

5.1 仿真发散：多半是结构几何出问题

仿真一跑就发散的经典场景是：数值在几步之内就直接涨到NaN。这个问题出现在结构建模阶段，十有八九是某个compartment太细太长，导致轴向电阻过大，求解器数值不稳定。

解决办法有几个方向：第一，把过长的compartment细分成多段，保证每一段的长度不要超过直径的几倍，通常建议长细比在5以内；第二，检查Ra的值是否合理，有些默认值在文献里是欧姆·厘米，在GENESIS里要换算成欧姆·米，换算错了会差100倍；第三，减小仿真步长，但这不是根治办法，只能临时绕过。

5.2 信号断流：检查拓扑连接

如果仿真不发散，但电流注入后电压毫无反应，或者某段树突电压恒等于静息值，那大概率是拓扑连接出了问题。优先检查所有相邻compartment之间的AXIAL消息是否双向都加上了，其次检查是不是有compartment只在命名空间上被创建，但从来没有和任何邻居连接。

我可以给你一个排查命令：在GENESIS里用showmsg查看某个compartment所有消息，快速确认连接是否齐全。

showmsg dend1

如果输出里只有VMEQ而没有AXIAL，那这个compartment就是“电隔离”的，信号自然传不过去。

5.3 数值异常：单位与浮点精度

还有一种情况是仿真能跑，但波形出现小幅高频振荡或者幅度异常。这种问题最烦人，因为不容易定位。我见过最多的原因是单位换算不一致。判断技巧是：把模型里所有长度和半径值统一列出来，看它们的量级是否符合形态学常识。比如一个树突直径填了2e-6米是2微米，合理；如果填了2e-9，那就是纳米量级，显然错了。

另外，GENESIS里的浮点精度默认是单精度，如果你某个compartment的几何量极其小，单精度可能会被截断成0。这种情况建议把几何量量级控制在1e-3到1e-12之间，避免极端值。

5.4 高频问题速查表

我把上面这些内容整理成一个速查表，运行时遇到问题可以先对照一遍。

最后再分享一个我个人的习惯，也是踩过不少坑之后总结出来的：不管项目多急，建完结构之后，一定先花十分钟做一个最小化测试——不加任何通道，不调任何参数，只给胞体打一针短促电流，看整个细胞的被动电学行为是否符合基本常识。这个测试能过滤掉至少八成的结构问题。等你以后跑大网络模型的时候就会发现，一个干净可靠的单细胞结构，比一堆精心调过的参数值值钱得多。