tTIS阈上刺激可工程化吗？

传统单细胞动作电位阈值是指能够引起单个神经元产生全或无锋电位的最小刺激强度（常用电流密度或电场强度表示）。该阈值由轴突起始段高密度钠通道的再生性去极化决定，并受膜时间常数（典型值 10–20 ms）、刺激波形上升速率及膜电位状态的影响。在离体脑片实验中，低频（<100 Hz）均匀电场引起锥体细胞动作电位的阈值约为 10–30 V/m；而高频（kHz）电场因膜的低通滤波特性，阈值急剧升高至 200 V/m 以上。这一概念是衡量神经元内在兴奋性的金标准，也是理解外源电刺激能否直接诱发神经放电的物理基础。论文中的计算模型（如图1C所示的单室 AdEx 神经元）正是基于这一原理，通过注入电流模拟电场对兴奋性神经元的极化作用。

图1：TI刺激的多尺度研究框架：电场计算、离体海马γ振荡记录及神经网络模型

A1：TI刺激电极放置示意图。两对电极分别置于左右半球（FT7/P7 和 FT8/P8），施加频率略有差异（1 kHz 与 1.005 kHz）的正弦电流。

A2：人体头模型计算出的电场幅值空间分布（前后方向）。颜色越红表示电场越强，可见皮层区域（靠近电极）电场最高（~80 V/m），深部较低。

A3：电场调幅深度（amplitude modulation depth）的空间分布。深部中线区域调幅深度最高（>50%），皮层区域调幅很浅（~15%）。

B1：离体海马脑片实验装置示意图。均匀电场通过两个平行银丝施加于脑片，记录电极位于 CA3a 区，参考电极置于浴液中。

B2：卡巴胆碱诱导的 γ 振荡原始迹线（局部场电位），显示稳定的振荡。

B3：γ 振荡的时频图（频谱随时间变化），证明振荡频率（20–40 Hz）和功率在记录期间保持稳定。

C1：网络模型结构图。包含 800 个兴奋性神经元（红色）和 250 个抑制性神经元（蓝色），连接概率和突触权重依类型而异。

C2：网络模型模拟的 γ 振荡（局部场电位），与离体实验形态相似。

C3：神经元发放光栅图，展示兴奋性（红）和抑制性（蓝）神经元在 γ 振荡周期中的时序关系。

作用：为后续的剂量‑效应实验和多尺度整合提供了物理、实验和计算的基础框架。

HUIYING

网络振荡调制阈值概述

网络振荡调制阈值是指能够显著改变神经元群体节律性场电位（如γ振荡，20–40 Hz）功率或同步性所需的最小外源刺激强度。该阈值并非固定数值，而是依赖于网络状态、刺激波形和统计判据的操作定义。在海马脑片卡巴胆碱诱导的γ振荡模型中（图1B显示实验记录位置及振荡稳定性），研究者通过比较刺激前后局部场电位功率的比值（经 t 检验或 ANOVA）来确定阈值。与单细胞阈值不同，网络阈值可以利用突触传递、兴奋‑抑制平衡及网络共振等涌现机制，在较低刺激强度下即可产生可观测的效应。例如，5 Hz 正弦电场仅需 1–5 V/m 即可显著调制γ振荡的动态变化，远低于该频率下的单细胞动作电位阈值。

HUIYING

静态调制比概述

静态调制比定义为刺激期间（通常取后段稳态）振荡功率与刺激前基线功率的比值，用于量化电刺激引起的平均功率增强或抑制。计算时，对局部场电位进行时频分析，在目标频带（如γ频段）内平均功率，再与基线相除。若比值 >1 表示兴奋，<1 表示抑制。该指标不关注刺激波形内部的动态变化，因此适用于非调幅或直流样刺激。论文中，非调幅 2 kHz 正弦波在 80 V/m 时静态调制比达到 1.15（图2D1展示了不同波形下静态调制随电场强度的变化），表明γ功率平均增加 15%。静态调制主要由持续去极化或超极化驱动，对 GABA慢抑制不敏感，甚至移除 GABA后静态调制反而增强（图3A3）。

HUIYING

动态调制概述

动态调制用于捕捉刺激波形周期性引起的振荡功率节律性波动，通常计算刺激周期内“峰”区间与“谷”区间的功率比值。对于低频正弦波，比较正电场半周期与负电场半周期的γ功率；对于调幅波，比较包络峰值段与谷值段的γ功率（图2A2示意了调幅波包络与γ功率的相位关系）。该指标反映了网络能否跟随刺激的低频包络（如 5 Hz）产生交替的兴奋/抑制，其实现依赖于 GABA等慢适应机制（时间常数 ~50 ms），该机制快于包络周期（200 ms）但慢于载频周期。论文中，5 Hz 调幅 1 kHz 载频在 60 V/m 时动态调制比达 1.33（图2D2），表明γ功率在包络峰时比谷时高 33%。动态调制是 TI 刺激实现“功能性选择性”的核心生理基础。

图2：海马脑片γ振荡及网络模型中的调制结果

A2：展示了5 Hz 调幅波（载频 1 kHz）的波形，以及 γ 功率随包络峰谷波动的动态调制概念。

B1：5 Hz 正弦波刺激下，γ 功率的原始迹线及频谱。显示功率随电场极性交替增减（动态调制），但平均功率不变。

B2：100 Hz 非调幅正弦波刺激，γ 功率在刺激开始后逐渐上升并维持（静态调制），无节律性波动。

B3：2 kHz 非调幅正弦波刺激，需更高电场强度（80–100 V/m）才能产生静态调制。

C1：5 Hz 调幅、载频 100 Hz 的刺激结果，可见明显的动态调制（功率随包络节律波动）。

C2：载频1 kHz 的调幅波，动态调制阈值约为 60 V/m。

C3：载频2 kHz 的调幅波，动态调制阈值升至 80 V/m。

D1：静态调制比随电场强度变化的汇总曲线（不同波形对比）。横坐标为电场强度，纵坐标为静态调制比。可见非调幅 100 Hz 和 2 kHz 需较高场强才能产生显著静态调制。

D2：动态调制比随电场强度变化的汇总曲线。调幅波（各载频）均能产生显著动态调制，且载频越低阈值越低；5 Hz 正弦波也产生动态调制，但机制不同（直接跟随电场极性）。

E1（模型结果）：计算模型复现的静态调制比曲线，与实验（D1）趋势一致。

E2（模型结果）：计算模型复现的动态调制比曲线，与实验（D2）一致。

作用：系统展示了波形参数（频率、调幅与否）对γ振荡调制的定量影响，验证了网络模型的预测能力。

HUIYING

敏感性和选择性的物理与生理边界概述

敏感性和选择性的物理与生理边界是指 TI 刺激能够有效且特异性地调控深部脑区所需满足的临界条件。物理边界包括：深部靶区的调幅电场成分需超过网络振荡调制阈值（如 60 V/m），对应头皮总电流约 167 mA；同时，皮层非调幅电场应尽量低于其静态调制阈值以避免脱靶效应。生理边界要求：神经元膜时间常数 ≤1 ms（主要指向轴突）才能对 kHz 载频有效解调，且存在功能性 GABA抑制以实现动态调制。论文通过整合头模型（图1A显示电场和调幅深度的空间分布）与网络模型（图3B展示多尺度预测流程）发现，当深部调幅成分达到阈值时，皮层非调幅电场（~80 V/m）已超过其静态调制阈值（~60 V/m），因此选择性边界无法独立满足——TI 必然同时激活皮层（静态调制）与深部（动态调制），只是调制模式不同。

HUIYING

阈上刺激概述

阈上刺激在经典神经生理学中指刺激强度达到或超过单细胞动作电位阈值，能可靠引发动作电位。在本文语境中，该概念被扩展至网络水平：当刺激强度超过网络振荡调制阈值时，即可称为功能阈上刺激。对于低频正弦波，1–5 V/m 已能显著调制γ振荡，但该强度远低于单细胞动作电位阈值（10–30 V/m），因此属于功能阈上但单细胞阈下。而对于 TI 所需的调幅 kHz 电场（60 V/m），该强度不仅超过网络动态调制阈值，也超过多数皮层轴突的单细胞动作电位阈值（约 10–30 V/m），因此对深部和皮层均属于经典阈上刺激。论文中的计算模型（图3A2显示单位活动在 80 V/m 时显著增加）证实了这一点。使用时必须明确“阈上”所指的效应尺度（单细胞放电 vs 网络振荡功率）和脑区。

HUIYING

临床研究（离体实验与计算模型）

研究方法

离体海马脑片：从 3–5 周龄雄性 Wistar 大鼠制备 400 μm 厚横向海马切片，在界面记录室中用 20 μM 卡巴胆碱诱导 CA3a 区γ振荡（20–40 Hz）。

电刺激施加：通过两个平行银丝在脑片上产生均匀电场，施加 2 秒的测试波形（低频 5 Hz 正弦、非调幅 100 Hz/2 kHz 正弦、5 Hz 调幅且载频 0.1/1/2 kHz 的波形），电场强度逐级递增（图1B1）。

记录与分析：用玻璃微电极记录局部场电位，计算静态调制比（刺激期平均γ功率/基线）和动态调制比（峰区间/谷区间功率比），通过线性混合效应模型和 ANOVA 确定统计显著性阈值（p<0.05）。

计算头模型：基于 1 mm 分辨率人体头部模型，有限元法模拟 TI 电极对（FT7/P7 与 FT8/P8）施加 167 mA、1 kHz 和 1.005 kHz 电流时全脑电场分布，提取每个体素的峰值电场强度和调幅深度（图1A）。

网络模型：构建 1050 个单室 AdEx 神经元（800 兴奋性 + 250 抑制性，含 GABA和 GABA 突触），用头模型输出的电场作为刺激输入，模拟γ振荡的静态和动态调制（图1C）。改变膜时间常数（0.4/1/3 ms）和 GABA存在与否，观察对敏感性和选择性的影响。

研究结果

波形依赖性阈值：

5 Hz 正弦：动态调制阈值 >1 V/m，静态调制不显著。

非调幅 100 Hz：静态调制阈值 7–10 V/m；非调幅 2 kHz：静态调制阈值 80–100 V/m。

5 Hz 调幅波：载频 100 Hz 时动态调制阈值 5 V/m；载频 1 kHz 时 60 V/m；载频 2 kHz 时 80 V/m。载频越高，敏感性越低（图2D1–D2）。

关键生理参数：

膜时间常数 τm = 1 ms 可拟合实验数据；τm = 0.4 ms 时敏感性大幅提高（所需电流降至 83 mA）；τm = 3 ms 时几乎不响应（图3A1）。

GABA是动态调制和选择性的必要条件：移除后调幅波的动态调制减弱，而非调幅波的静态调制增强（图3A3）。

图3：泛化网络模型及TI刺激的多尺度预测

A1：在 5 Hz 调幅、1 kHz 载频刺激下，改变膜时间常数（τm = 0.4, 1, 3 ms）和是否包含 GABAb抑制对动态调制的影响。蓝色实线（τm=1 ms, 有 GABAb）拟合实验数据；红色（τm=0.4 ms）敏感性更高；绿色（τm=3 ms）几乎不响应。虚线表示移除 GABAb后动态调制大幅减弱。

A2：不同电场强度下（10, 40, 80 V/m）兴奋性神经元的平均发放率变化。显示在80 V/m 时发放率显著增加，证明该强度已是单细胞阈上。

A3：在非调幅 1 kHz 正弦波刺激下，改变 τm和 GABAb对静态调制的影响。移除 GABAb（虚线）反而增强静态调制，说明 GABAb 抑制了持续的去极化效应。

B：多尺度模型工作流程图。从左到右：头模型计算每个体素的电场峰值与调幅深度 → 输入网络模型 → 输出静态调制与动态调制图 → 与实验阈值比较。

C1：τm = 0.4 ms 时的全脑预测图。上行：动态调制空间分布（深部热点明显）；下行：静态调制空间分布（皮层全域高值）。

C2：τm= 1 ms（实验拟合值）的预测图。深部动态调制热点仍存在，但强度略降；皮层静态调制依然显著。

C3：τm = 3 ms 的预测图。几乎无动态调制热点，静态调制也极弱。 每列左右分别代表有/无 GABAb。结果显示：无论参数如何，深部动态调制高的区域，皮层静态调制也高，无法实现空间选择性。

作用：揭示了 TI 刺激敏感性和选择性的生理决定因素（τm 和 GABAb），并定量预测了人体中无法避免的皮层‑深部共存效应。

人体 TI 预测：

使深部调幅成分达到 60 V/m 需要头皮总电流 ~167 mA，此时皮层非调幅电场 ~80 V/m，超过皮层静态调制阈值（~60 V/m）。

无论参数如何变化（τm或 GABA 调整），深部动态调制高的区域总是伴随皮层静态调制显著（图3C1–C3）。

TI 无法实现真正的空间选择性（深部激活而皮层无效应），只能产生调制模式的空间分化（深部动态、皮层静态）。

HUIYING

总结

本研究通过离体海马γ振荡实验、人体头模型电流场模拟及生物物理网络模型，系统阐明了时间干涉（TI）刺激的敏感性和选择性边界。敏感性（所需头皮电流）由神经元膜时间常数决定，轴突（τm ≤1 ms）是 TI 解调 kHz 载频的关键靶点，人体实现有效深部调制需上百毫安级电流，远高于传统 tACS。选择性依赖于 GABA 介导的慢适应机制，该机制使网络能跟随 5 Hz 包络产生动态调制，而皮层由于近场单电极对主导，形成非调幅 kHz 电场，仅产生静态调制。然而，电场叠加的物理规律决定了当深部调幅达到阈值时，皮层非调幅强度必然超过其静态调制阈值，因此 TI 无法实现深部单独激活——其本质是“调制模式选择性”而非“空间选择性”。这些发现为优化电极蒙太奇、利用脑区异质性参数以及设定临床合理预期提供了理论依据。未来若能将 TI 与更复杂的波形设计或多阵列电极结合，或可进一步改善深部与皮层的效应分离。