基于计算模拟将精神分裂症前额叶微环路抑制减弱与脑电生物标志物关联研究

《PLOS Computational Biology》：Linking reduced prefrontal microcircuit inhibition in schizophrenia to EEG biomarkers in silico

　　

精神分裂症中由小清蛋白（parvalbumin, PV）阳性中间神经元介导的皮层抑制降低被认为与前额叶皮层（prefrontal cortex, PFC）处理功能受损及脑电信号改变相关，例如在奇异刺激任务中出现Oddball失配负波（mismatch negativity, MMN）降低。近期研究亦提示生长抑素（somatostatin, SST）阳性中间神经元抑制功能存在损失。然而，在人体实验中建立抑制减弱与MMN降低之间的因果关联目前尚不可行。为克服上述挑战，研究人员构建了健康与精神分裂症状态下的人前额叶微环路精细化模型，其中PV和SST中间神经元抑制的降低幅度受死后患者数据约束，并模拟了基线态及Oddball响应期间的神经元放电活动与脑电信号。研究结果表明，PV中间神经元抑制减弱可解释精神分裂症患者中观察到的MMN振幅降低现象，且存在效应较低的阈值，这与高危患者中的表现一致。与之相对，SST中间神经元抑制减弱对MMN振幅无显著影响。进一步分析显示，两种类型抑制的降低均为解释精神分裂症静息态脑电改变所必需：SST中间神经元抑制减弱导致宽带功率增加，而PV与SST中间神经元抑制减弱均导致频谱由α频段向β频段右移。该研究 thus 将精神分裂症中PV和SST中间神经元抑制降低与特异性脑电生物标志物相关联，为利用非侵入性脑信号改善患者分层与早期检测提供了理论基础。

该研究聚焦于精神分裂症神经环路机制与脑电生物标志物之间的关联这一核心科学问题。精神分裂症作为一种严重的精神疾病，其发病机制至今尚未完全阐明。现有研究表明，皮层功能障碍涉及多个脑区的协同活动异常，同时在细胞和微环路层面也存在显著改变。其中，中间神经元介导的抑制性调控障碍被认为是关键病理机制之一。特别值得关注的是，PV阳性中间神经元通过及时抑制锥体（pyramidal, Pyr）神经元活动，对高频γ振荡（20–80 Hz）产生重要调节作用，进而支持前额叶皮层的认知功能；而SST阳性中间神经元则主要通过侧向抑制调节基线活动，对低频θ和α振荡（4–12 Hz）具有调控作用。尽管死后研究已证实精神分裂症患者前额叶皮层存在PV和SST表达降低，但如何在活体人体中建立这些细胞层面改变与脑电生物标志物之间的直接关联，始终面临技术和伦理层面的双重限制。另一方面，脑电图技术作为一种非侵入性检测手段，在精神分裂症的早期诊断和患者分层中展现出重要应用前景，尤其是MMN作为预测精神病发作风险的敏感指标，以及静息态脑电中θ、β功率增加伴α功率降低等特征性改变，但其背后的微环路机制亟待阐明。基于此，研究人员开展了这项计算模拟研究，旨在通过构建精细化的人前额叶微环路模型，桥接细胞机制与脑电表型之间的认知鸿沟，为精神分裂症的机制理解和临床转化提供新视角。

该研究主要由Henriquez、Hjorth、Kozlov、Barbour、Bhalla、Blaabjerg、S?derman和Müller等人合作完成，发表于《PLOS Computational Biology》杂志。研究的技术路线主要涵盖以下关键方法：基于Yao等人先前构建的人皮层第2/3层微环路模型进行适应性改造，依据人前额叶皮层实验数据调整神经元比例（Pyr神经元72%、SST中间神经元5%、PV中间神经元9%、VIP中间神经元14%）；利用NEURON仿真平台、LFPy工具包及高性能并行计算进行神经元放电活动与局部场电位模拟；采用四维球导体模型（脑组织、脑脊液、颅骨、头皮）计算微环路生成的脑电信号；运用Ornstein-Uhlenbeck随机过程模拟背景兴奋性输入；基于死后患者基因表达数据约束PV中间神经元输出机制（突触性和张力性抑制电导降低22%）和输入机制（Pyr至PV的NMDA受体NR2A亚基介导的突触电导降低20%），以及SST中间神经元抑制降低；通过模拟Oddball听觉任务期间的刺激响应重现事件相关电位，并采用FOOOF库对脑电功率谱进行周期性和非周期性成分分解。

研究结果部分包含以下主要内容：

**健康前额叶微环路模型的基线活动与脑电特征**：研究人员构建了包含1000个神经元的人前额叶微环路模型，四种神经元类型均采用详细重建的人源形态学。模型基线放电率与人源在体记录数据基本一致：Pyr神经元为0.94±0.05 Hz（实验值0.66±0.51 Hz），PV中间神经元为7.44±0.42 Hz（实验值2.63±2.55 Hz）。模拟脑电的功率谱密度在α频段（8–12 Hz）呈现显著峰值，并具有1/f特征，与人前额叶静息态脑电谱特性相符。

**模拟Oddball响应与精神分裂症微环路中的MMN**：研究人员依据灵长类动物实验记录，通过激活Pyr神经元群及PV、SST中间神经元群模拟Oddball响应，重现了100–160 ms刺激后约6–11 Hz的Pyr神经元放电率变化特征。在健康模型中，基线放电率与标准刺激响应无显著差异。当PV中间神经元抑制以10%–50%梯度降低时，基线和Oddball响应放电率呈线性至超线性递增：每降低10%抑制，基线放电率稳步增加约10%，但30%至40%降幅时增幅跃升至15%，40%至50%时更达25%；Oddball响应放电率每降低10%抑制增加约6%，30%至40%及40%至50%降幅时超线性跳跃达10%。值得注意的是，抑制降低对基线活动的影响大于Oddball响应，导致信噪比（signal-to-noise ratio, SNR）随抑制降低而下降。单独分析表明，PV中间神经元突触输出机制对基线放电率、响应放电率和SNR的贡献均大于NMDA输入机制。在事件相关电位方面，MMN振幅随PV中间神经元抑制降低而线性递减，每降低10%抑制约减少5%振幅，但40%至50%降幅时出现超线性跳跃至-11%（2倍效应）。PV中间神经元突触输出机制对事件相关电位的贡献（-11.7%）约为NMDA机制（-5.9%）的两倍。

**PV与SST中间神经元抑制对Oddball响应和MMN的影响比较**：对比50%降低PV或SST中间神经元抑制的效应发现，两者均使基线Pyr神经元放电率显著增加且幅度相近（PV：+86%；SST：+99%），但Oddball响应增加幅度差异显著（PV：+47%；SST：+14%）。信噪比降低幅度以SST减弱更为显著（-43% vs -21%）。关键发现在于，SST中间神经元抑制降低并未改变MMN振幅（Cohen's d = -0.12），即便与PV抑制降低联合作用时，其MMN效应亦与单独PV降低无显著差异。这一结果的机制在于SST抑制降低对基线事件相关电位（-14%）和Oddball事件相关电位（-9%）产生相近影响，在差值波形计算中相互抵消。

**PV与SST中间神经元抑制对静息态脑电特征的影响**：PV中间神经元抑制降低（40%）导致功率谱峰值频率从11.2 Hz右移至12.7 Hz（+13%），低β频段（12–20 Hz）功率显著增加（+41%），但非周期性宽带功率无显著变化。SST中间神经元抑制降低（40%）同样引起峰值频率右移至12.6 Hz（+13%），但同时还导致多频段功率增加：α频段+11%，β频段+125%；其非周期性特征表现为指数降低（-17%）和宽带功率增加（+36%）。PV与SST联合降低时，功率谱呈现两者效应的叠加，峰值频率大幅右移至14.4 Hz（+28%），非周期性指数降低（-23%），宽带功率增加（+21%），但同时伴有δ频段绝对功率降低（-27%）。

在讨论部分，研究人员系统阐释了研究发现的理论与临床意义。首先，研究确立了PV中间神经元抑制降低作为精神分裂症患者MMN振幅降低的主要机制，验证了既往假设；而SST中间神经元抑制降低由于同时影响基线和Oddball响应事件相关电位，其效应在MMN差值计算中被抵消，故对MMN无独立贡献。值得关注的是，MMN降低存在显著的阈值效应：50%抑制降低导致25%的MMN振幅降低，与患者数据一致；而10%–40%的较轻度抑制降低仅引起较轻微效应，与临床高危人群平均15%的MMN降低相符，提示该人群可能存在尚未进展为典型症状的PV中间神经元抑制功能损害。这一发现为基于神经环路机制的患者分层和早期检测提供了量化依据。在静息态脑电方面，PV或SST抑制降低均能解释α至β的频谱右移现象，而SST抑制降低对于重现宽带功率增加不可或缺，这与SST中间神经元调控低频活动的已知功能一致。研究人员还探讨了模型简化与局限：当前模型仅涵盖第2/3层且神经元数量约为实际的1/7，导致模拟脑电和事件相关电位振幅偏小，建议未来可通过归一化生物标志物或振幅缩放因子应用于临床数据；此外，模型中50%的NMDA NR2A亚基表达降低被简化为整体20%的突触NMDA电导降低，且基因表达百分比变化可能转化为更大的生理学效应。研究人员还指出，未来需扩展至多层面模型以研究不同层间中间神经元比例差异的影响，并探索多脑区交互作用的建模。

研究结论：该研究通过精细化计算模拟，机制性地将PV和SST中间神经元抑制降低与精神分裂症的特异性脑电生物标志物相关联。PV中间神经元抑制降低主要解释MMN振幅降低和α至β频谱右移，SST中间神经元抑制降低则主要贡献于静息态宽带功率增加和频谱右移，两者联合作用产生更为广泛的脑电表型改变。这些发现为精神分裂症的非侵入性生物标志物开发、患者分层和早期检测提供了基于神经环路机制的理论框架和可检验的定量预测。