用于比较帕金森病中130 Hz丘脑底核深部脑刺激与10 Hz基础脉冲周期性视觉输入的皮质-基底节-丘脑（CBT）in silico测试平台

《Frontiers in Computational Neuroscience》：An in silico cortico-basal ganglia-thalamic testbed for comparing visual and deep brain stimulation in Parkinson’s disease

　　

引言：深部脑刺激（DBS）是帕金森病（PD）的既定干预措施，而采用10 Hz基础脉冲的规定性周期视觉输入的网络水平效应仍不确定。
方法：研究人员开发了一个基于电导的皮质-基底节-丘脑（CBT）模型，包含11个Hodgkin-Huxley神经元群体和一个规定性周期视觉输入单元（RET），用于在同一校准PD网络上比较130 Hz丘脑底核（STN）DBS与存档的周期视觉输入条件。主要条件各使用N = 10次独立运行。附加分析包括非配对0–5 mA/cm² DBS振幅序列、匹配种子的内侧苍白球至丘脑皮质中继（GPi→TC）敏感性分析（在0.8、1.0和1.2 × 名义耦合下，每尺度N = 10），以及匹配种子的三参数±20%局部敏感性分析（每参数-条件单元N = 5）。
结果：DBS将模型内部STN突触电流beta指标从0.0184 ± 0.0060降至0.0018 ± 0.0004 a.u.（Holm校正 p = 1.1 × 10^–5），而周期视觉输入增加了丘脑皮质中继（TC）和网状丘脑（RE）放电，并增强皮质锥体（PY）gamma（30–80 Hz）和窄带38–42 Hz功率。皮质PY放电在PD中为4.86 ± 0.84 Hz，DBS下为5.28 ± 0.87 Hz（与PD相比p = 0.285），周期视觉输入下为10.33 ± 0.85 Hz（与PD相比Holm校正 p = 1.18 × 10^–10）。振幅序列显示1 mA/cm²下STN反应非单调；从2至5 mA/cm²，平均STN beta值处于0.00186–0.00212 a.u.的狭窄范围。在GPi→TC分析中，TC自发放电在30/30次运行中保持0 Hz，所有750,000个保存的电压样本均为有限且阈下（–72.06至–54.74 mV；无0-mV交叉）。在三参数分析中，DBS原始STN-beta抑制保持在88.9%–95.8%，周期视觉输入使TC放电相对于名义PD增加3.70–4.29倍。
讨论：这些发现显示在所测试剂量、GPi→TC和三参数局部域中具有持续性，但不能确定诱发中继保真度、全局参数稳健性、生物学机制或临床疗效。

论文解读：

研究背景：帕金森病（Parkinson’s disease, PD）是进行性神经退行性疾病，黑质纹状体多巴胺丢失会重组直接通路和间接通路活动，并促进皮质-基底节-丘脑（cortico-basal ganglia-thalamic, CBT）网络异常同步。过量beta频段活动（13–30 Hz）与僵直和运动迟缓相关，可作为帕金森状态的电路水平标志。左旋多巴仍是主要药物治疗，但长期使用可能并发运动波动和异动症。对于晚期PD，丘脑底核（subthalamic nucleus, STN）高频深部脑刺激（deep brain stimulation, DBS）是既定手术选择，可改变基底节输出并抑制病理beta活动。节律性感觉刺激也被研究用于影响分布式神经活动，但采用10 Hz基础脉冲的规定性周期视觉输入的网络水平效应仍不确定。计算模型可跨互补层次描述PD，因此研究人员整合皮质、丘脑、纹状体和基底节，构建可重构Hodgkin-Huxley（HH）CBT测试平台，在同一PD网络上比较130 Hz STN-DBS与10 Hz基础脉冲规定性周期视觉输入，评估皮质下频谱、丘脑和皮质放电及频谱输出，目的是生成可定量检验的模型假设，而非确立临床治疗原则。

研究开展与总体结论：研究人员开发了多状态、基于电导的CBT测试平台，并在清醒PD状态下比较DBS与周期视觉输入。结果表明，DBS与较低的模型内部STN beta功率相关，并伴随TC自发放电降低至0 Hz；周期视觉输入则增加TC、RE和皮质PY活动及38–42 Hz功率，但未抑制STN beta代理指标。这些结果定义了两类扰动如何重分配网络活动的可检验、模型条件假设，但不确立诱发中继能力、临床疗效、治疗优越性、全局参数稳健性或通用神经调控设计原则。

关键技术方法：研究人员在Python 3.9中实现基于电导的HH CBT测试平台，含11个神经元群体及规定性周期视觉输入单元RET；采用后向欧拉与Newton-Raphson迭代，时间步0.01 ms，单次运行3,000 ms并排除前500 ms；用130 Hz、0.06 ms/相、电荷平衡双相脉冲模拟STN-DBS，主要振幅5 mA/cm²，并做0–5 mA/cm²非配对序列；以N = 10随机运行和Welch t检验、Holm校正、Mann-Whitney U检验比较；进行匹配种子GPi→TC 0.8、1.0、1.2 × 名义耦合及三参数±20%敏感性分析。无样本队列来源，全部为计算模拟。

研究结果：
3.1 模型表征与证据检查：报告三个互补证据水平，包括选定参考行为复现、丙泊酚下方向性外部比较及随机实现中的数值和运行间行为；这些检查不确立未测试条件性能或全局稳健性。
3.1.1 选定参考行为的复现：健康态纹状体快放电中间神经元（sFS）、直接通路中等棘神经元（dMSN）和间接通路中等棘神经元（iMSN）放电率分别为5.53 ± 1.13、2.27 ± 0.28和1.83 ± 0.32 Hz，符合低放电纹状体状态；多巴胺耗竭使GPi原始beta功率增加3.8倍、iMSN增加4.0倍，dMSN放电从2.3升至5.4 Hz、iMSN从1.8升至9.1 Hz，但STN电流代理未显示经典临床STN局部场电位增加。
3.1.2 方向性跨状态检查：丙泊酚相关低频增强：在PD网络中，模拟丙泊酚使原始STN delta/theta功率（1–8 Hz）从1.25 ± 0.30 × 10^–3 a.u.升至68.2 ± 31.9 × 10^–3 a.u.，约54倍变化，方向与所报道的STN低频活动增加一致。
3.1.3 数值行为与运行间离散：在0.01 ms隐式时间步下，分析运行保持数值有限；RET生成器产生100 ms间隔的1 ms脉冲（10 Hz）；N = 10时，原始STN beta功率变异系数在PD、DBS和周期视觉输入下分别为32.8%、21.1%和19.4%，TC放电变异系数在非零放电条件下≤4%。
3.2 病理beta振荡的详细复现：校准多巴胺耗竭参数后，GPi和间接通路iMSN的beta频段功率增加，dMSN和iMSN放电增加；模型内部STN输出电流代理未显示经典临床STN-LFP增加，说明这些为选定基准检查。
3.3 耦合过程的数值表征：对“shifting and stacking”耦合程序进行数值表征，时间平均合成输入近似参考静态驱动并保留时变波动，支持耦合实现数值兼容性，但不表示重现生物轴突延迟统计。
3.4 DBS与周期视觉输入的比较模型反应：在5 mA/cm²、130 Hz STN-DBS下，原始STN beta功率从0.0184 ± 0.0060降至0.0018 ± 0.0004 a.u.，描述性beta抑制比为90.2%，TC自发放电从8.6 ± 0.3 Hz降至0 Hz；周期视觉输入未抑制STN beta代理，原始STN beta功率为0.0423 ± 0.0082 a.u.，但使TC放电从8.6 ± 0.3 Hz升至35.3 ± 0.6 Hz，并招募RE放电从0至6.56 ± 4.15 Hz。
3.5 皮质输出与靶向敏感性分析：皮质PY平均放电在健康、PD、DBS和周期视觉输入下分别为5.83 ± 0.56、4.86 ± 0.84、5.28 ± 0.87和10.33 ± 0.85 Hz；周期视觉输入增加PY放电、beta、30–80 Hz gamma和38–42 Hz功率，DBS未显著改变这些皮质指标。0–5 mA/cm²振幅序列显示非单调STN反应，1 mA/cm²时STN beta代理升高而TC放电降至0 Hz；2–5 mA/cm²时平均STN beta处于0.00186–0.00212 a.u.窄范围，TC自发放电均0 Hz。GPi→TC局部敏感性分析中，TC放电在30/30次运行中保持0 Hz，750,000个电压样本均有限且阈下，原始STN beta在各尺度为0.001931 ± 0.000333 a.u.。三参数±20%分析中，DBS原始STN-beta抑制为88.9%–95.8%，周期视觉输入使原始STN beta增加76.2%–230.9%，TC放电增加3.70–4.29倍。

讨论部分总结：4.1 同一CBT架构包含觉醒、非快速眼动睡眠和丙泊酚麻醉参数化；丙泊酚条件提供方向性跨状态检查，DBS与视觉比较限于清醒PD。4.2 DBS表现为模型内部STN beta代理降低和状态依赖性TC静息；全分辨率诊断和GPi→TC敏感性分析未支持持续性去极化阻滞，而更符合阈下自发静息。4.2.2 周期视觉输入表现为丘脑和皮质招募，但未使PY beta正常化或抑制STN beta代理。4.3 局限性包括：靶向敏感性未覆盖全局多变量稳健性；皮质高beta同步较临床预期更强；简化NREM未重现睡眠相关病理beta；STN beta代理方向性与人类STN-LFP经典方向不同；周期视觉输入增加STN频谱能量并非人类视觉刺激降低STN-LFP beta的证据；RET绕过显式外侧膝状核动态，且模型电流密度不等同临床DBS设备电流。

研究结论翻译：研究人员开发了一个多状态、基于电导的CBT测试平台，并用其比较130 Hz STN-DBS与采用10 Hz基础脉冲的规定性周期视觉输入在同一清醒PD网络上的效应。在所检查操作点，DBS与较低的模型内部STN beta功率相关，而周期视觉输入增加丘脑和皮质活动但未抑制该指标。TC = 0输出在1–5 mA/cm²以及匹配种子GPi→TC耦合±20%局部变化中持续存在；所有保存的TC电压保持阈下，该模式不支持持续性去极化阻滞，并与受测模拟中的自发静息一致。皮质分析显示视觉输入增加PY放电和38–42 Hz功率，而DBS未显著恢复PY放电相对于PD。匹配种子±20% GABA_A衰减、STN–GPe耦合和STN电流分析在受测局部域内保留了所有四个预设结局。这些结果定义了关于两种扰动如何重新分配网络活动的可检验、模型条件假设；它们不确立诱发中继能力、临床疗效、治疗优越性、全局参数稳健性或通用神经调控设计原则，这些问题需要更广泛的多变量敏感性分析和实验验证。