《PLOS Computational Biology》：Coupled beta and high-frequency oscillations emerge from synchronized bursting in a minimal model of the parkinsonian subthalamic nucleus
编辑推荐：
帕金森病（Parkinson's disease，PD）丘脑底核（subthalamic nucleus，STN）记录到的局部场电位（local field potentials，LFPs）表现出独特的跨尺度频谱特征：β频段振荡（13–30 Hz）增强，并与高
帕金森病（Parkinson's disease，PD）丘脑底核（subthalamic nucleus，STN）记录到的局部场电位（local field potentials，LFPs）表现出独特的跨尺度频谱特征：β频段振荡（13–30 Hz）增强，并与高频振荡（high-frequency oscillations，HFOs，200–400 Hz）相耦合，HFO振幅锁定于β周期相位。这种相位-振幅耦合（phase-amplitude coupling，PAC）已被认定为帕金森病状态的一种有前景的生物标志物，然而尚无生物物理模型解释其产生机制、HFO频率的決定因素，或解释为何在用药状态下的STN中HFO可在无β调制的情况下存在。本研究表明，具有循环耦合的异质性兴奋性Izhikevich神经元群体可产生三种动力学状态：（i）异步紧张性放电；（ii）异步爆发放电，其中神经元各自爆发产生宽带HFO功率，但缺乏群体水平的相干PAC；（iii）同步爆发放电，由此产生β-HFO PAC。这些状态由两个具有生物物理可解释性的参数控制，它们分别反映多巴胺耗竭的互补效应：一个反映内在神经元兴奋性的变化，另一个反映突触耦合强度的变化。在该模型中，从异步爆发到同步爆发的转变模拟了帕金森病状态下STN病理性神经元活动的出现。HFO峰值频率在二维参数空间中连续变化，提示了临床上观察到的用药状态间从慢HFO（200–300 Hz）向快HFO（300–400 Hz）转变的可能机制。同步化转变的性质取决于基线兴奋性，从低兴奋性下爆发与同步的急剧共同出现，到中等兴奋性下的解耦转变——此时爆发比例趋于饱和而群体同步化随耦合继续增强。STN群体中的异步-同步爆发转变及β-HFO PAC同样出现在模型的扩展版本中，该版本包含与外侧苍白球（globus pallidus externa，GPe）样抑制性神经元群体的交互耦合，β频率随两群体间突触耦合延迟的变化而变化。该模型为未来临床和实验研究提供了可检验的预测，提供了中尺度LFP特征如何映射到微观神经元动力学的数值解析，并为未来环路水平模型提供了计算构建模块，这些模型可能为针对患者特异性STN动力学状态的脑刺激策略提供信息。
论文解读
一、研究背景与问题
帕金森病（Parkinson's disease，PD）患者丘脑底核（subthalamic nucleus，STN）的局部场电位（local field potentials，LFPs）具有独特的跨尺度频谱特征：β频段（13–30 Hz）振荡显著增强，并与200–400 Hz的高频振荡（high-frequency oscillations，HFOs）形成相位-振幅耦合（phase-amplitude coupling，PAC），即HFO振幅随β相位系统性变化。临床研究已证实，STN中β振荡功率与运动症状严重程度相关，而PAC强度与运动迟缓和强直严重程度相关。值得注意的是，在用药（ON）状态下，HFO虽然存在但不再受β相位调制，表现为"从β调制中解放出来"的功能性促运动特征。临床文献将STN HFO分为慢成分（sHFO，200–300 Hz）和快成分（fHFO，300–400 Hz），多巴胺能药物引起功率从sHFO向fHFO的特征性转移，且sHFO/fHFO功率比与临床运动评分相关。
尽管这些临床证据充分，但三个关键问题悬而未决：第一，什么生物物理机制能在STN样神经元群体中产生β调制的HFO？第二，什么决定了HFO频率，为何其在用药状态间发生转移？第三，在用药状态下HFO为何能在无β调制的情况下存在，何种神经元放电动力学的转变构成了病理性β-HFO PAC的基础？既往基底节计算模型主要关注通过多区域环路产生β振荡，但均未涉及HFO生成或跨频率耦合问题。基于二次整合-放电（quadratic integrate-and-fire，QIF）神经元的精确平均场模型虽能产生β嵌套快振荡，但复位电位在宏观方程中消失，无法区分内在爆发放电群体与紧张性放电群体——这正是帕金森病STN的核心区分点。因此，有必要在放电神经元群体水平上研究这一问题。
二、研究内容与核心发现
研究人员构建了一个由异质性兴奋性Izhikevich（自适应QIF）神经元组成的单群体最小模型，包含慢适应和循环耦合，以两个具有多巴胺底物的参数刻画多巴胺耗竭的效应：中位复位电位（c?，控制内在兴奋性）和突触耦合强度（J）。模型产生三种动力学状态：异步紧张性放电、异步爆发放电（无群体水平PAC的宽带HFO功率）和同步爆发放电（产生β-HFO PAC）。从异步爆发到同步爆发的转变模拟了从用药状态到帕金森病状态的转变。HFO峰值频率在二维参数空间中连续变化，为临床上观察到的sHFO/fHFO转移提供了机制解释。模型还复现了Meidahl等人报告的HFO包络振幅尖峰触发平均（spike-triggered average，STA）调制模式，与单单元记录一致。扩展的STN-GPe环路模型证实，即使缺乏STN内循环兴奋，β-HFO PAC机制仍可涌现，β频率随环路延迟变化而HFO频率保持不变。
三、主要技术方法
研究采用Izhikevich（自适应QIF）神经元模型构建STN单群体网络（5000个神经元），以洛伦兹分布设定复位电位异质性，采用稀疏随机AMPA型兴奋性突触连接。通过数值模拟系统扫描（c?，J）二维参数空间，计算平均放电率变异系数（CV）作为同步性指标、相对β功率、Tort调制指数（modulation index，MI）量化PAC强度，并采用替代数据检验显著性。模型扩展版本引入STN-GPe双向兴奋-抑制环路，包含传导延迟参数。所有分析基于合成LFP信号和尖峰序列的频谱分析。
四、研究结果
单神经元现象学：紧张性放电向爆发放电的转变。 复位电位c控制锋电位后超极化的深度，决定神经元以紧张性单峰放电还是爆发模式放电。低c值（?65 mV）产生规则单峰放电；c增大超过约?57 mV后，神经元产生爆发。每个神经元内在具有快慢双时间尺度：爆发内放电频率处于HFO范围（200–400 Hz），爆发间间隔处于β范围（13–30 Hz）。两个时间尺度均随c连续变化，建立了群体水平连续HFO频率景观的基础。
网络动力学：从紧张性放电到异步和同步爆发放电。 同步化转变的性质取决于群体相对于单神经元紧张性-爆发边界的分布。低中位兴奋性（c?=?64 mV）时，群体大部分低于爆发阈值，耦合作为有效去极化驱动将神经元推向爆发边界，爆发与同步在临界耦合处共同出现，无异步爆发中间状态。中等兴奋性（c?=?58 mV）时，约37%神经元为内在爆发放电者，β功率随耦合增加而上升，但爆发比例和PAC保持基线水平；更高耦合下更多神经元被招募进入爆发并锁定在一起，产生强β振荡和β-HFO PAC，爆发比例饱和而同步化继续增强。高兴奋性（c?=?53 mV）时，几乎所有神经元均为内在爆发放电者，无耦合时各自产生HFO频率爆发内放电但爆发起始时间不相关，LFP显示HFO功率峰但无β峰或PAC；耦合增加逐步相关化爆发时间，同步化后产生尖锐β峰和强PAC。
兴奋性-耦合相图。 二维参数空间映射显示，β-HFO PAC在同步化转变处涌现，其出现跟随CV=0.5边界；HFO中心频率在整个参数空间连续变化，主要由内在兴奋性决定，更负的c?产生200–300 Hz HFO，较不负的c?将中心移向300–400 Hz。临床上的sHFO/fHFO分类可能对应连续频谱表面的不同区域而非两种独立振荡机制。
跨同步化转变的频谱度量。 慢振荡峰值频率随耦合增加而降低，下降速率和起始频率取决于c?。HFO中心频率主要由c?决定，耦合仅起次要调制作用。PAC调制指数在转变阈值以下接近零，超过阈值后急剧上升。
从无PAC的快HFO到有PAC的慢HFO：双参数转变。 异步爆发状态（c?=?52 mV，J=0.050）对应用药状态，具有fHFO功率但无PAC；同步爆发状态（c?=?56 mV，J=0.166）对应帕金森病状态，β功率和β-HFO PAC显著，HFO频率移至sHFO范围。两状态平均放电率约44 Hz和40 Hz，处于STN神经元在体记录范围。同步状态中HFO包络振幅的STA显示清晰的β频率调制，与Meidahl等人的实验发现一致。爆发内ISI峰对应HFO频谱峰，爆发间ISI峰对应β频谱峰，证实HFO是独立生成的节律而非β谐波。
滞后与双稳态。 在c?=?62和?54 mV时，耦合上扫和下扫遵循不同分支，转变显示滞后现象，对应异步或部分同步稳定态与强同步稳定态共存的双稳态区间；c?=?58 mV时两方向重合，无滞后窗口。
STN-GPe环路中的候选机制。 在无STN内循环兴奋的STN-GPe环路模型中，增加耦合强度驱动STN群体从异步爆发转向同步爆发，β功率和β-HFO PAC在同步化边界共同上升，HFO中心频率仍由内在兴奋性决定。β振荡频率随传导延迟从3 ms增至10 ms而从约27 Hz降至约18 Hz，而HFO频率保持在约200 Hz不变。GPe群体通过速率调制提供同步化反馈，自身不爆发，HFO及其与β的耦合仍为爆发性STN群体的产物。
五、讨论与结论
该模型为帕金森病STN的β-HFO频谱特征提供了最小尖峰网络机制解释。核心结论是：具有慢适应和循环兴奋耦合的异质性QIF尖峰神经元群体产生三种动力学状态，由两个生物物理可解释参数控制；从异步爆发到同步爆发的转变模拟了从用药状态到帕金森病状态的转变。HFO频率由单神经元爆发内放电率决定并随内在兴奋性连续变化，解释了临床上sHFO/fHFO的区分可能是连续频谱表面的采样而非两种独立机制。β-HFO PAC仅在群体爆发同步化时出现，HFO无β调制标记异步爆发状态而非爆发的缺失。模型预测：HFO中心频率在患者间应形成连续体而非两个离散群组；PAC需要同步化爆发，异步爆发状态（有HFO无PAC）是模型特有的；沿帕金森病方向移动预测β峰值频率下移及频谱变化的有序序列；若β由网络反馈产生，其频率应随传导延迟变化而HFO频率不变；双稳态窗口附近的去同步化刺激可产生超刺激持续效应且依赖基线兴奋性；β同步性是爆发的同步而非单峰放电的同步，这一区别以嵌套HFO及其相位耦合为标志。该框架为基于患者特异性动力学状态的脑深部刺激（deep brain stimulation，DBS）策略提供了计算基础，提示β爆发反映的是嵌套的双时间尺度同步爆发放电模式而非单纯的群体放电率增加，刺激效果可能取决于其相对于同步爆发的时间位置。
