训练神经元：构建海马CA1锥体神经元的精细化、最新、多用途生物物理模型

《Frontiers in Computational Neuroscience》：How to train your neuron: developing a detailed, up-to-date, multipurpose model of hippocampal CA1 pyramidal cells

　　

解剖学和生物物理细节丰富的神经元模型已被广泛用于研究这些细胞中的信息处理。大多数研究侧重于理解特定现象，而旨在同时捕获多种细胞过程的通用模型仍然罕见，尽管此类模型对于预测神经元在更复杂、自然条件下的行为是必需的。在本研究中，研究人员旨在开发一种详细的、数据驱动的、多用途的海马CA1锥体神经元生物物理模型。研究人员利用该细胞类型可获得的广泛的形态学、生物物理和生理学数据，并建立了一个系统化的模型构建和验证工作流程，该流程依赖于研究人员最近开发的软件工具。该模型基于高质量形态学重建，并包含一组经过整理的多样化离子通道模型。在纳入现有的离子通道分布约束后，剩余的自由参数使用Neuroptimus工具进行优化，以匹配从体细胞全细胞记录中提取的多种电生理特征。使用HippoUnit进行的验证确认了该模型复制关键电生理特征的能力，包括体细胞对电流输入的电压响应、突触电位和反向传播动作电位的衰减，以及斜向树突中的非线性突触整合。该模型还包含活跃的树突棘，这些树突棘可以显式建模，也可以通过将其生物物理机制并入母树突的机制中来建模。研究人员发现，神经元行为的许多方面不受棘建模细节水平的影响，但准确建模非线性突触整合需要显式建模棘。这种匹配多种实验约束的数据驱动CA1锥体细胞模型，是研究这些细胞活动和可塑性的通用工具，也可以作为海马网络精细模型的可靠组成部分。这种构建和验证多用途模型的系统化方法也应适用于其他细胞类型。

研究背景与目的：解剖学和生物物理细节丰富的神经元模型已广泛用于研究信息处理，但多数已有模型为解释特定现象而构建，缺乏能够同时捕获多种细胞过程的通用模型，且常存在欠约束、参数手工调整、缺乏系统验证等问题。对已有CA1锥体细胞（CA1 pyramidal cell, CA1 PC）模型的系统比较显示，没有任何已有模型能同时复制HippoUnit所覆盖的所有相关行为。因此，研究人员旨在利用大鼠CA1 PC丰富的形态学、生物物理和生理学数据，以及新开发的自动化参数优化和模型验证工具，构建并验证一个数据驱动的多用途生物物理模型，为研究该细胞的活动、可塑性以及海马网络建模提供可靠基础。

研究方法：研究使用Judit Makara提供的大鼠CA1 PC体细胞全细胞电流钳记录（6个细胞，每细胞3次测量，排除1个双峰放电细胞）作为优化目标，利用eFEL提取电生理特征，以Neuroptimus中的CMAES算法优化漏电导、离子通道密度与部分动力学参数；模型基于Megías等的形态学重建并补充Bianchi等的轴突，包含13种电压门控和/或钙门控离子通道以及自定义NMDA受体模型；用HippoUnit的五项测试进行多范式验证，并通过F-factor或显式建模方式处理树突棘。

研究结果：

（1）模型构建流程概述：研究者在每轮拟合中同时运行20次随机初始化优化，最终得到20个优化模型。多数参数未一致贴近边界；与体细胞反应相关的参数在各模型间高度一致，而主要影响树突行为的参数变异性较大，因为优化目标未包含树突特征。

（2）电压门控离子通道模型：研究者系统筛选并校准了13种通道模型，包括体细胞、树突和轴突具有不同动力学的Na⁺通道、多种K⁺通道及钙通道，并根据文献设置非均匀分布。通过逐个阻断通道的实验发现，各通道在模型中的作用与实验预期一致，例如阻断Na⁺通道完全阻止放电，阻断A型K⁺通道消除树突中反向传播动作电位的距离依赖性衰减，阻断M型和D型K⁺通道改变放电模式和兴奋性，阻断BK、SK和慢AHP钾通道则特征性地改变后超极化和放电模式。

（3）模型验证：HippoUnit测试显示，20个模型在体细胞特征测试中误差较小，能复现动作电位形状、放电频率、电压sag等特征；反向传播动作电位测试中所有模型均表现为强传播型，与实验中的强传播细胞群体一致；突触后电位衰减测试中，需要将轴向电阻设为50 Ohm×cm才能匹配实验数据，传统100 Ohm×cm会导致过度衰减；去极化阻滞测试中多数模型未通过，因为高电流注入下轴突仍产生小幅度spikelet，使膜电位未达到HippoUnit所定义的平坦阻滞标准，但这一现象与部分实验记录一致，且需要更精细的轴突模型来解释。

（4）使用斜向树突整合测试确定NMDA受体模型的正确电压依赖性：研究者系统比较了8种已发表的NMDA受体电压依赖性，发现它们均不能完全匹配Losonczy和Magee实验中的斜向树突整合特征。通过调整sigmoid函数的半激活电压和斜率，研究人员自定义了NMDA受体模型（半激活约?6 mV，斜率参数18 mV），使模型在5个同步输入时出现由树突Na⁺尖峰驱动的跳跃式反应，并正确复现了非线性整合程度和EPSP幅度。

（5）特定离子通道的作用：通过完全阻断通道和±20%密度扰动实验，研究者发现Na⁺通道密度和性质对多数测试特征具有主导影响；其他扰动多具有选择性效应。N型钙通道与BK钾通道、L型钙通道与慢AHP钾通道之间因钙浓度耦联而表现出类似阻断效应。

（6）显式建模树突棘的影响：研究者比较了无棘F-factor模型、仅显式建模接受突触输入的棘、以及全棘显式模型。在体细胞特征、反向传播动作电位、PSP衰减和去极化阻滞测试中，三种模型表现几乎一致；但在斜向树突整合测试中，突触靶向棘头时NMDA成分幅度更高、时程更长，产生实验观察到的超线性整合，而突触直接靶向树突干时非线性不足。仅建模突触棘与全棘模型结果一致，说明采用混合策略可在较低计算成本下准确模拟突触整合。

讨论与结论：研究者强调，该系统化建模过程使模型可适用于多种应用，也可作为海马网络模型的关键组件。经过更新的离子通道模型集合本身是有价值的资源，但需要更多高质量、细胞类型特异的离子通道数据。自动化参数拟合仍需大量人工设定，包括机制选择、参数边界和固定参数；仅优化体细胞特征不足以约束树突相关参数。HippoUnit测试有助于模型约束，但存在局限：去极化阻滞测试未能区分轴突spikelet与真实去极化阻滞；斜向整合测试的某些特征可能因NMDA慢成分而高估误差。树突棘方面，F-factor方法对非突触信号有效，但突触输入整合需要显式建模棘。研究结论可翻译为：研究人员希望建模和实验学界能从这一新的CA1锥体神经元多用途模型中获益。根据验证结果，该模型能够捕捉这些细胞在体细胞电流注入下的行为（去极化阻滞可能除外）、树突上被动和主动信号的传播（包括朝向和离开胞体），并正确整合斜向树突上的突触输入。需要这些特性的计算机实验可以安全地在该模型上进行。另一方面，通过使用更详细的轴突模型并纳入更多特定实验约束，该模型的轴突行为可能会得到改进；该模型在在体条件下以及作为网络模型组件的行为也需要评估。研究人员的多数工作流程应可推广至其他细胞类型甚至其他神经系统；构建忠实数据驱动模型的最关键要求和潜在瓶颈是获得所有必要类型的实验数据，因为这些数据的质量和数量决定模型构建和验证的成功。