Name: 应用交互式、患者特异性模型对神经纤维束进行深部脑刺激治疗的靶向性 (视频) | JoVE Journal
Uploaded: 2018-08-12T12:30:00.000+00:00
Duration: 14 min 14 s
Description: 了解如何使用交互式患者特异性模型模拟深部脑刺激对神经纤维束的影响。该方法利用SCIRun、FreeSurfer和3DSlicer软件实现实时计算建模，为帕金森病等运动障碍疾病的个性化治疗策略提供重要支持。

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

本项目得到了美国国立卫生研究院资助项目UH3、NS095554的支持。技术援助由科学计算与成像研究所的整合生物医学计算中心提供，并部分得益于根据美国国立卫生研究院P41-GM103545项目开发的软件，该软件属于整合生物医学计算中心。
感谢科学计算与成像研究所的Lexie Floor和Nathan Galli在投稿视频制作与编辑方面提供的帮助，同时感谢Theresa Lins在手稿准备过程中的协助。

1. 图像处理
<ol>
	<li>T1 MRI 分割
	<ol>
		<li>下载并安装 FreeSurfer17。假设已获取 T1 MRI 图像，并且其格式为 DICOM 或 NIFTI。</li>
		<li>输入以下命令以在 FreeSurfer 中创建患者目录并添加其 T1 图像：recon-all -s patientName -i /Full/Path/To/nii</li>
		<li>输入以下命令以运行 FreeSurfer 的自动分割：recon-all -s patientName -all</li>
		<li>使用 mri_convert 将 &#34;aseg.auto.mgz&#34; 从 FreeSurfer 专有的 .mgz 格式中转换出来。</li>
	</ol>
	</li>
	<li>弥散加权成像（DWI）
	<ol>
		<li>下载并安装 Slicer18。假设已进行弥散加权成像，并以一系列 DICOM 文件形式提供。 
		注意：一个优秀的命令行工具用于执行此重建是 FMRIB 软件库（FSL）19。本实验方案使用 FSL 在重建前完成平面回波、运动及涡流畸变校正。本分析中使用的 DWI 数据包含 41 个不同弥散敏感梯度方向的 DICOM 文件序列。</li>
		<li>在 Slicer 中使用 DWIConverter 模块将 41 次独立采集重建为单个体积数据。指定 DWI 数据所在的 Input Dicom Data Directory。选择 DicomToNrrd 单选按钮，并提供输出 DWI 体积的名称。</li>
		<li>使用 Diffusion Weighted Volume Masking 模块创建张量掩膜，以消除 DTI 估计中的背景噪声。将阈值参数保持默认值 0.5。</li>
		<li>使用上一步创建的掩膜作为 Diffusion Tensor Mask，通过 DWI to DTI Estimation 将该单个 DWI 体积转换为弥散张量（DTI）体积。将估计参数设置为最小二乘法，并确保勾选 Shift Negative Eigenvalues 选项。</li>
		<li>将新生成的 DTI 体积以 NRRD 格式保存。</li>
	</ol>
	</li>
</ol>
2. 有限元模型生成
<ol>
	<li>SCIRun
	<ol>
		<li>下载并安装 SCIRun 5 版本（http://www.sci.utah.edu/software/scirun.html）。SCIRun 是一个用于建模、仿真和可视化科学问题的通用问题求解环境。本项目中使用该计算工作台来创建、求解并与所开发的模型进行交互。 
		注意：SCIRun 网络通过连接执行特定任务的独立模块序列构建而成。将数据传递通过由简单模块构成的网络，可实现更复杂的处理与仿真。</li>
	</ol>
	</li>
	<li>构建电极几何结构 
	注意：本项目中建模的电极为 Medtronic 3387 DBS 引线20。它是一种圆柱形引线，包含四个触点，每个触点高度为 1.5 mm，直径为 1.27 mm，触点之间相距 1.5 mm。四个触点之间由非导电材料隔开。
	<ol>
		<li>根据所列尺寸，为 DBS 引线中的每个组件创建封闭的表面网格。通过使用两种基本几何体——圆柱体和球体——来构建完整的引线。 
		注意：引线几何结构可在多种 3D 建模软件中创建。</li>
		<li>创建一个厚度为 0.5 mm 的包膜层表面网格，包裹整个引线。</li>
	</ol>
	</li>
	<li>创建全脑网格
	<ol>
		<li>加载第 1 部分中创建的大脑表面，作为有限元法（FEM）的外边界。</li>
		<li>创建两个围绕 DBS 引线的同心盒状表面，以控制电极周围的网格密度。 
		注意：DBS 引线和盒状表面应沿正 z 轴对齐，电极杆的底端位于坐标原点。这对于后续电极的旋转和平移操作至关重要。</li>
		<li>创建一个点云，在网格的每个不同区域中放置一个点。网格的区域包括：杆段、四个触点、盒状表面以及大脑表面。附加在这些点上的数值代表该区域内每个四面体单元的最大体积约束。通过调整这些约束来控制电极内部及周围的网格密度。</li>
		<li>使用 &#34;InterfaceWithTetGen&#34; 模块创建全脑网格。设置输入端口如下： 
		端口 1：大脑表面（来自 1.1） 
		端口 2：区域属性点云，体积约束（来自 2.3.3） 
		端口 3：无 
		&#8203;端口 4：电极模型表面（来自 2.2）</li>
	</ol>
	</li>
	<li>交互式电极放置 
	注意：此步骤的目标是能够在大脑内部任意位置轻松移动电极，并使系统自动更新生成新的网格和仿真结果。
	<ol>
		<li>使用两个级联的 CreateGeometricTransform 模块，在 X 和 Y 平面内旋转电极。点击 Rotate，将第一个模块的 Axis Y 设置为 1.00，第二个模块的 Axis X 设置为 1.00。旋转角度通过底部滑块 Rotate Theta (degrees) 设置。</li>
		<li>为每一个需要变换的变换操作和对象连接一个 TransformMeshWithTransform 模块。</li>
		<li>创建一个 GenerateSinglePointProbeFromField 模块，并将其几何端口连接到 ViewScene 模块。在 ViewScene 窗口中按住 Shift 键并单击左键来移动球体。释放鼠标后，点探针模块将更新球体的位置。 
		注意：该位置用于确定施加于电极几何结构的平移变换。</li>
		<li>使用一个 TransformMeshWithTransform 模块，根据点探针的位置对电极几何结构的每一部分进行平移。</li>
		<li>将两个旋转变换以及来自点探针控件的平移应用于除大脑表面以外的所有对象。 
		注意：这三个变换应在每个对象送入 InterfaceWithTetGen 模块之前依次完成，以便网格生成算法能够知晓需在大脑内部的新位置和方向上对电极几何结构进行网格划分。每次在大脑内移动电极时，都将重建一个新的网格。</li>
	</ol>
	</li>
</ol>
3. 生物电场计算
<ol>
	<li>设置电导率参数
	<ol>
		<li>加载第1步中创建的DTI体积数据，并使用MapFieldDataFromSourceToDestination模块的线性（加权）选项，将数据映射到第2步中通过InterfaceWithTetGen生成的四面体网格上。 
		注意：此步骤中可实现任何各向异性电导率估计。</li>
		<li>使用CalculateFieldData将电极杆区域的电导率设为1e-6，接触区域的电导率设为1e6。设置完所有电导率值后，将场数据输入至BuildFEMatrix模块。</li>
	</ol>
	</li>
	<li>设置电流/电压源及边界条件
	<ol>
		<li>将第2步网格化过程中使用的脑表面数据输入至SetFieldDataToConstantValue模块，并将其值设为0，从而在整个脑表面创建一个电压汇。然后使用InsertVoltageSource模块将电压值施加到由InterfaceWithTetGen输出的四面体网格上。</li>
		<li>对于单极刺激，通过在其中一个接触面使用ReportFieldInfo模块计算FieldCenter值，在该接触位置中心创建一个单一点。将该点的场数据设为-1伏特，然后使用InsertVoltageSource模块将该点源插入四面体网格中。 
		注意：此时存在两个四面体网格副本，一个在脑表面设置了电压值，另一个在接触区域内包含一个点源。</li>
		<li>合并两个InsertVoltageSource模块的输出结果，并将其与BuildFEMatrix的输出一同输入至ApplyFEMVoltageSource模块。</li>
	</ol>
	</li>
	<li>求解正向问题 
	注意：ApplyFEMVoltageSource模块的输出在端口1提供刚度矩阵，在端口2提供右端项矩阵。这两个组分是求解方程组所必需的。
	<ol>
		<li>将ApplyFEMVoltageSource的两个输出连接至SolveLinearSystem模块，选择共轭梯度法作为求解方法，Jacobi作为预条件子。该模块的输出为解矩阵，即四面体网格中每个节点处的电压值。</li>
		<li>使用SetFieldData模块将电压解矩阵赋值到四面体网格，以实现电压分布的三维可视化。</li>
	</ol>
	</li>
	<li>等值面可视化
	<ol>
		<li>将SetFieldData模块的结果输入至ExtractSimpleIsosurface模块。打开用户界面中的Quantity选项卡，将均匀分布的等值数设为10。</li>
		<li>使用ShowField模块并将其连接至ViewScene模块，以透明方式显示等值面，并采用彩虹色图对电压解的最小值和最大值进行缩放显示。</li>
	</ol>
	</li>
</ol>

<ol><li>Benabid, A. L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Chronic+electrical+stimulation+of+the+ventralis+intermedius+nucleus+of+the+thalamus+as+a+treatment+of+movement+disorders%5BTitle%5D)+AND+%22Journal+of+Neurosurgery%22%5BJournal%5D)">Chronic electrical stimulation of the ventralis intermedius nucleus of the thalamus as a treatment of movement disorders.</a> Journal of Neurosurgery. 84 (2), 203-214 (1996).</li>
<li>Limousin, P., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Effect+of+parkinsonian+signs+and+symptoms+of+bilateral+subthalamic+nucleus+stimulation%5BTitle%5D)+AND+%22The+Lancet%22%5BJournal%5D)">Effect of parkinsonian signs and symptoms of bilateral subthalamic nucleus stimulation.</a> The Lancet. 345 (8942), 91-95 (1995).</li>
<li>Schiff, N. D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Behavioural+improvements+with+thalamic+stimulation+after+severe+traumatic+brain+injury%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">Behavioural improvements with thalamic stimulation after severe traumatic brain injury.</a> Nature. 448 (7153), 600-603 (2007).</li>
<li>Vandewalle, V., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Stereotactic+treatment+of+Gilles+de+la+Tourette+syndrome+by+high+frequency+stimulation+of+thalamus%5BTitle%5D)+AND+%22Lancet%22%5BJournal%5D)">Stereotactic treatment of Gilles de la Tourette syndrome by high frequency stimulation of thalamus.</a> Lancet. 353 (9154), 724(1999).</li>
<li>Mayberg, H. S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Deep+brain+stimulation+for+treatment-resistant+depression%5BTitle%5D)+AND+%22Neuron%22%5BJournal%5D)">Deep brain stimulation for treatment-resistant depression.</a> Neuron. 45 (5), 651-660 (2005).</li>
<li>Hashimoto, T., Elder, C. M., Okun, M. S., Patrick, S. K., Vitek, J. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Stimulation+of+the+subthalamic+nucleus+changes+the+firing+pattern+of+pallidal+neurons%5BTitle%5D)+AND+%22Journal+of+Neuroscience%22%5BJournal%5D)">Stimulation of the subthalamic nucleus changes the firing pattern of pallidal neurons.</a> Journal of Neuroscience. 23 (5), 1916-1923 (2003).</li>
<li>York, M. K., Wilde, E. A., Simpson, R., Jankovic, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Relationship+between+Neuropsychological+Outcome+and+DBS+Surgical+Trajectory+and+Electrode+Location%5BTitle%5D)+AND+%22J.+Neurol.+Sci%22%5BJournal%5D)">Relationship between Neuropsychological Outcome and DBS Surgical Trajectory and Electrode Location.</a> J. Neurol. Sci. 287 (1-2), 159-171 (2009).</li>
<li>Machado, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Deep+brain+stimulation+for+Parkinson%27s+disease%3A+surgical+technique+and+perioperative+management%5BTitle%5D)+AND+%22Movement+Disorders%22%5BJournal%5D)">Deep brain stimulation for Parkinson's disease: surgical technique and perioperative management.</a> Movement Disorders. 21, S247-S258 (2006).</li>
<li>Volkmann, J., Moro, E., Pahwa, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Basic+algorithms+for+the+programming+of+deep+brain+stimulation+in+Parkinson%27s+disease%5BTitle%5D)+AND+%22Movement+Disorders%22%5BJournal%5D)">Basic algorithms for the programming of deep brain stimulation in Parkinson's disease.</a> Movement Disorders. 21 (S14), S284-S289 (2006).</li>
<li>Contarino, M. F., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Directional+steering%3A+A+novel+approach+to+deep+brain+stimulation%5BTitle%5D)+AND+%22Neurology%22%5BJournal%5D)">Directional steering: A novel approach to deep brain stimulation.</a> Neurology. 83 (13), 1163-1169 (2014).</li>
<li>Pollo, C., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Directional+deep+brain+stimulation%3A+An+intraoperative+double-blind+pilot+study%5BTitle%5D)+AND+%22Brain%22%5BJournal%5D)">Directional deep brain stimulation: An intraoperative double-blind pilot study.</a> Brain. 137 (7), 2015-2026 (2014).</li>
<li>Willsie, A. C., Dorval, A. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Fabrication+and+initial+testing+of+the+%C2%B5DBS%3A+a+novel+deep+brain+stimulation+electrode+with+thousands+of+individually+controllable+contacts%5BTitle%5D)+AND+%22Biomedical+Microdevices%22%5BJournal%5D)">Fabrication and initial testing of the µDBS: a novel deep brain stimulation electrode with thousands of individually controllable contacts.</a> Biomedical Microdevices. 17 (56), 9961(2015).</li>
<li>Oliveria, S. F., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Safety+and+efficacy+of+dual-lead+thalamic+deep+brain+stimulation+for+patients+with+treatment-refractory+multiple+sclerosis+tremor%3A+a+single-centre%2C+randomised%2C+single-blind%2C+pilot+trial%5BTitle%5D)+AND+%22The+Lancet+Neurology%22%5BJournal%5D)">Safety and efficacy of dual-lead thalamic deep brain stimulation for patients with treatment-refractory multiple sclerosis tremor: a single-centre, randomised, single-blind, pilot trial.</a> The Lancet Neurology. 16 (9), 691-700 (2017).</li>
<li>Butson, C. R., Cooper, S. E., Henderson, J. M., Wolgamuth, B., Mcintyre, C. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Probabilistic+Analysis+of+Activation+Volumes+Generated+During+Deep+Brain+Stimulation%5BTitle%5D)+AND+%22Neuroimage%22%5BJournal%5D)">Probabilistic Analysis of Activation Volumes Generated During Deep Brain Stimulation.</a> Neuroimage. 54 (3), 2096-2104 (2011).</li>
<li>Basser, P. J., Mattiello, J., LeBihan, D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((MR+diffusion+tensor+spectroscopy+and+imaging%5BTitle%5D)+AND+%22Biophys.+J%22%5BJournal%5D)">MR diffusion tensor spectroscopy and imaging.</a> Biophys. J. 66 (1), 259-267 (1994).</li>
<li>Gullmar, D., Haueisen, J., Reichenbach, J. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Influence+of+anisotropic+electrical+conductivity+in+white+matter+tissue+on+the+EEG%2FMEG+forward+and+inverse+solution.+A+high-resolution+whole+head+simulation+study%5BTitle%5D)+AND+%22NeuroImage%22%5BJournal%5D)">Influence of anisotropic electrical conductivity in white matter tissue on the EEG/MEG forward and inverse solution. A high-resolution whole head simulation study.</a> NeuroImage. 51 (1), 145-163 (2010).</li>
<li>Fischl, B. <a target="_blank" href="https://www.ncbi.nlm.nih.gov/pubmed/22248573">FreeSurfer.</a> Neuroimage. 62 (2), 774-781 (2012).</li>
<li>Fedorov, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((3D+Slicer+as+an+Image+Computing+Platform+for+the+Quantitative+Imaging+Network%5BTitle%5D)+AND+%22Magn.+Reson.+Imaging%22%5BJournal%5D)">3D Slicer as an Image Computing Platform for the Quantitative Imaging Network.</a> Magn. Reson. Imaging. 30 (9), 1323-1341 (2012).</li>
<li>Jenkinson, M., Beckmann, C. F., Behrens, T. E., Woolrich, M. W., Smith, S. M. <a target="_blank" href="https://www.ncbi.nlm.nih.gov/pubmed/21979382">FSL.</a> Neuroimage. 62 (2), 782-790 (2012).</li>
<li><a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=%22Medtronic+DBS+3387%2F3389+Lead+Kit+Manual%22">Medtronic DBS 3387/3389 Lead Kit Manual</a>. , (2018).</li>
<li>Baker, J. L., et al. <a target="_blank" href="https://www.ncbi.nlm.nih.gov/pubmed/27582298">Robust modulation of arousal regulation, performance and frontostriatal activity through central thalamic deep brain stimulation in healthy non-human primates.</a> Journal of Neurophysiology. 116 (5), 2383-2404 (2016).</li>
<li>Fregni, F., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Treatment+of+major+depression+with+transcranial+direct+current+stimulation%5BTitle%5D)+AND+%22Bipolar+Disorders%22%5BJournal%5D)">Treatment of major depression with transcranial direct current stimulation.</a> Bipolar Disorders. 8 (2), 203-204 (2006).</li>
<li>Hodaie, M., Wennberg, R. A., Dostrovsky, J. O., Lozano, A. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Chronic+anterior+thalamus+stimulation+for+intractable+epilepsy%5BTitle%5D)+AND+%22Epilepsia%22%5BJournal%5D)">Chronic anterior thalamus stimulation for intractable epilepsy.</a> Epilepsia. 43 (6), 603-608 (2002).</li>
<li>Rosenow, F., Lüders, H. <a target="_blank" href="https://academic.oup.com/brain/article/124/9/1683/303186">Presurgical evaluation of epilepsy.</a> Brain. 124 (9), 1683-1700 (2001).</li>
</ol>

Christopher R. Butson 博士曾担任 NeuroPace、Advanced Bionics、Boston Scientific、Intelect Medical、St. Jude Medical 和 Functional Neuromodulation 的顾问。

本方案展示了一系列技术，可显著降低构建用于脑深部电刺激（DBS）的患者特异性模型所需的时间和计算成本，从而实现仿真结果的近实时反馈。快速反馈使得能够探索较大的参数空间，以更深入地理解参数变化如何影响模型结果。这些参数包括电极位置、激活触点的选择，以及在这些触点上的刺激波形的幅度、脉宽和频率。所提出工具的关键特性包括：1）提供一个简洁的用户界面，可调节模型参数，并对这些参数如何影响仿真结果进行近实时可视化；2）通过少量输入自动完成模型构建：包括脑表面、脑组织电导率张量，以及电极几何结构的表面表示。该自动化过程加快了针对不同脑部几何结构和组织电导率的多个个体构建患者特异性模型的速度，同时也便于评估将不同电极设计植入现有模型所产生的影响。本方案中描述的图像预处理步骤尚未完全自动化，可能需要长达一天的处理时间。然而，一旦完成，这些步骤生成的数据被视为静态数据，即在仿真过程中不会被修改。模型构建的自动化源于系统能够将这些数据自动应用于患者特异性的有限元模型（FEM），而无需人工干预。用于执行全部模型生成、仿真和可视化步骤的 SCIRun 网络只需构建一次，因此，为新患者生成患者特异性模型时，仅需重新执行图像预处理步骤即可。
建模流程在生成结果方面的性能提升，得益于将网格生成、生物电场计算以及解的可视化整合到单一软件环境中。采用现有的建模技术（如自适应网格生成），可在电极周围生成较高密度的网格，而在远离电极的区域生成较低密度的网格，从而缩短有限元法（FEM）模型的构建和求解时间。软件 SCIRun 还支持网格生成和生物电场计算的自动化。用户通过交互式控件自定义电极移动时，会触发根据更新后的电极位置重建新网格的过程，包括为新的电极位置调整边界条件和电导率值。
在将电极位置整合到有限元模型（FEM）之前，电极几何结构被视为脑容积内可自由移动的物体。这种建模方法在网格构建中的一个关键优势是可轻松地将多个电极插入模型中。例如，可将同一电极几何结构的第二个副本放置在距离原位置数毫米处，两个电极都将被包含在有限元模型中。在近期研究中，已有将两个电极近距离植入以治疗多发性硬化症震颤的案例13，并在非人灵长类动物实验中用于探索有效的刺激靶点21。使用多个电极的优势在于能够在更大区域内更精确地控制组织中产生的电场。通过在两个电极的激活触点上同时进行刺激，可将电流引导至目标区域，并避开可能引起不良副作用的区域。在目标靶点确切位置未知的情况下，对更大区域实现更精细的刺激控制，对于探索不同刺激位置尤为有用，这正是许多新兴脑深部电刺激（DBS）治疗应用所面临的情况。然而，由于参数空间本身已较大，增加电极数量进一步扩大了参数空间，因此确定实现治疗性刺激的参数比单电极情况更具挑战性。
我们设想，这种交互式建模工具可在脑深部电刺激（DBS）植入术的术前规划过程中提供帮助。关于电刺激在脑组织中作用范围的反馈信息，可使外科医生在手术计划中调整电极位置，从而对目标区域提供有效的治疗性刺激。尽管DBS治疗是开发该工具的主要动因，但本文介绍的技术也可应用于其他有限元生物电场模型，适用于不同的刺激或记录范式。例如，用于治疗抑郁症的经颅直流电刺激22，或用于癫痫治疗的深部电极应用23，在确定最佳刺激位置以实现治疗效果方面，面临着与DBS相同的挑战。皮层脑电图（electrocorticography）是一种在脑表面布设电极阵列以识别癫痫发作起始区域的记录技术，其挑战在于如何确定电极的放置位置，以有效记录脑内目标区域的信号24。所有这些应用均依赖于电极的位置，同时又需应对电流在脑组织中传导路径的不确定性。本文所介绍的技术降低了构建和使用计算模型的负担，同时为非建模专业背景的研究人员和临床医生在使用这些设备时提供了有意义的反馈信息。

脑深部电刺激（DBS）是一种已确立的治疗运动障碍的方法，适用于原发性震颤1和帕金森病2。该疗法目前也正在被研究用于治疗越来越多的其他疾病，包括创伤性脑损伤3、抽动秽语综合征4和抑郁症5。DBS 系统需要通过手术植入电极导线，以在特定脑区施加电刺激，从而调节正在进行的神经活动6。电极的位置和刺激参数都会影响提供治疗效果的神经回路调节。电极位置的微小差异可能会影响治疗窗口，甚至在达到治疗效果之前增加不良副作用发生的风险7,8,9。在实际临床中，通常难以预测刺激对神经活动的影响；因此，这一治疗窗口需由临床医生根据每位患者的具体情况，在编程调试刺激设备的过程中逐一确定8,9。随着新一代 DBS 设备的不断出现，这一过程正变得愈加复杂。例如，新型电极导线设计正在引入更多电触点10,11,12，在某些情况下，多个电极导线被植入彼此邻近的区域13。因此，有必要能够探索并预测 DBS 在日益庞大和复杂的参数空间中的作用效果。
计算建模与分析可用于在个体化基础上预测深部脑刺激（DBS）对生理和临床效果的影响。这些模型采用有限元建模（FEM）方法，构建脑组织以及植入电极生物物理特性的三维表征。FEM生物电场模型已成功用于预测DBS的效果14，但迄今为止，这些模型的生成过程耗时较长且计算成本高昂。随着DBS设备复杂性的不断增加，亟需开发新型建模技术。这些个体化患者模型应能够在改变电极位置或刺激参数时，提供接近实时的视觉反馈。用户可在数秒内获得关于电极位置和刺激设置的反馈信息，从而在几分钟内持续优化电极的放置。个体化特性的实现依赖于在构建FEM时整合患者的解剖结构、大脑的形态与尺寸，并应用其脑组织的生物物理特性，例如各向异性组织电导率。各向异性电导率描述了电流在不同脑区中的扩散方式，可通过非侵入性方法对全脑进行测量，其过程类似于常规磁共振成像（MRI）。
不使用患者特异性信息的脑深部电刺激（DBS）建模方法可以快速预测刺激效应，但由于采用了通用的脑组织几何结构和电导率值，其预测准确性较低。在此类方法中，所有患者均使用同一个有限元模型（FEM），且可预先计算预测的神经活动。而患者特异性模型由于需为每位个体单独构建新的有限元模型，因此无法通用化或预先计算。这类模型构建更为耗时，但准确性更高。目前有多个因素限制了这些模型的构建和使用速度：1）在建模流程初期修改参数（如电极位置）需要人工手动更新后续所有步骤；2）建模流程中的各个步骤之间难以有效集成，需在多个软件平台之间传递数据。通常情况下，我们希望评估多种不同情形，例如电极位置、刺激参数或电极设计。为了对这些变化对患者治疗效果的影响提供有意义的反馈，相关结果必须既准确又能够快速生成。
我们的目标是介绍一种构建患者特异性模型的新技术，该技术充分利用通用模型所获得的速度优势，并自动化处理流程中的多个步骤，从而创建一个交互式建模环境，能够提供关于深部脑刺激（DBS）效果的近实时可视化反馈。交互式仿真允许用户快速测试预测结果并获得输出，而无需专注于模型构建的细节。当需要探索的参数空间较大且这些参数对仿真结果的影响尚不明确时，这种方法尤为有益。我们将描述从磁共振成像（MRI）数据生成交互式、患者特异性有限元法（FEM）模型的处理流程中的各个步骤。利用本文所述的工具与技术，可显著降低构建FEM生物电场模型的时间成本，并为非建模专业背景的研究人员和临床医生提供一种便捷途径来使用这些模型。
本方案描述了如何从获取的磁共振成像（MRI）数据构建患者特异性的有限元模型，并模拟深部脑刺激（DBS）电极诱导的电场。构建此类模型的主要步骤包括：1）建立代表患者脑组织及植入电极的有限元模型（FEM）；2）向该有限元模型中添加脑组织的生物物理特性以及DBS设备的刺激参数；3）求解电极在模型中产生的电压分布。为构建用于DBS仿真的患者特异性模型，需要两种成像模态：T1加权MRI用于构建脑表面、脑室及特定核团的分割图像；弥散加权成像（DWI）——一种测量水分子弥散性的技术——用于估算全脑组织内的弥散张量15。将弥散张量转换为电导率张量，以逐体素方式量化组织的非均质性和各向异性生物物理特性16。通过求解泊松方程计算电极在全脑诱导的电压分布，该方程经有限元法（FEM）处理后简化为线性方程组 Ax = b，其中 A 为表示网格电导率和几何结构的刚度矩阵，x 为网格各节点处的电压解，b 则根据边界条件和电流源进行调整。

<table><tbody><tr><td>FreeSurfer</td><td>哈佛大学 Athinoula A. Martinos 生物医学成像中心</td><td></td><td>&lt;a href=&quot;https://surfer.nmr.mgh.harvard.edu/">https://surfer.nmr.mgh.harvard.edu/&lt;/a&gt;</td></tr><tr><td>3D Slicer</td><td>布莱根妇女医院，哈佛大学</td><td></td><td>&lt;a href=&quot;https://www.slicer.org/">https://www.slicer.org/&lt;/a&gt;</td></tr><tr><td>SCIRun&amp;nbsp;</td><td>犹他大学整合生物医学计算中心</td><td></td><td>&lt;a href=&quot;http://www.sci.utah.edu/cibc-software/scirun.html">http://www.sci.utah.edu/cibc-software/scirun.html&lt;/a&gt;</td></tr></tbody></table>

利用交互式患者特异性模型为深部脑刺激疗法定位神经纤维束

Targeting Neuronal Fiber Tracts for Deep Brain Stimulation Therapy Using Interactive, Patient-Specific Models

脑深部电刺激（DBS）通过植入电极对特定脑区进行刺激，是治疗运动障碍的成熟疗法，并正逐步应用于越来越多的其他疾病。计算建模已被成功用于预测DBS的临床效应；然而，随着DBS设备复杂性的不断增加，亟需开发新的建模技术。这些模型还需能够快速而准确地生成预测结果。本项目的目标是建立一个图像处理流程，将结构磁共振成像（MRI）和弥散加权成像（DWI）整合到一个可交互的、针对个体患者的模型中，以模拟DBS的作用效果。可在患者模型内放置虚拟DBS电极导线，并设定激活触点和刺激参数；当导线位置或方向发生变化时，系统可在近实时条件下（约10秒内）生成新的有限元网格并求解生物电场问题。该系统还支持模拟多个邻近导线，通过调节不同导线上阳极和阴极的配置实现电流导向。本文所介绍的技术降低了构建和使用计算模型的难度，同时为非建模专业背景的研究人员或临床医生提供了关于电极位置、电极设计及刺激配置影响的有意义反馈。

在本方案结束时，将提供构建患者特异性模型所需的全部输入信息：脑表面、电极几何结构以及电导率张量。同时还应已创建一个 SCIRun 网络，用于整合所有输入数据，构建有限元模型并模拟诱导产生的生物电场。该仿真模型支持电极在脑内的位置调整，并可修改诸如激活触点或刺激幅度等参数。
图1 展示了在多个电极位置下模拟和比较纤维激活预测的能力。本方案中未展示基于扩散成像生成纤维束的过程，但该分析可通过任何确定性示踪方法完成。与位置1相比，位置2在所有触点上均能在较低刺激强度下实现对目标纤维束更有效的激活。该分析在脑深部电刺激（DBS）手术规划阶段用于确定能够有效刺激目标纤维通路的电极植入路径。
该系统的一个新特性是能够快速切换模型中的电极几何结构，并能够模拟多个彼此靠近的电极导线。这两种方法均可用于更精确地控制电极周围电场的形态与导向。图2展示了在研究电极导线位置以刺激中央丘脑纤维通路并避开感觉核等邻近区域时，轴对称型与定向型脑深部刺激（DBS）导线在纤维激活预测方面的比较。对于相同的电极路径，可以看出定向型导线能够将纤维激活引导至目标区域，同时避免刺激非目标纤维。该分析可用于判断在何种情况下需要更换患者所接受的DBS治疗装置，以避免刺激邻近区域引发不必要的副作用。该方法在靶向刺激治疗帕金森病症状的丘脑底核时尤为有用，同时可避免刺激邻近的内囊。
<img alt="神经纤维激活与刺激强度关系图；电极位置；神经科学研究。" class="xfigimg" src="/files/ftp_upload/57292/57292fig1.jpg" data-altupdatedat="1788250015000" data-alt="图 1"> 
图 1：电极位置与目标纤维激活预测的比较。 针对两种不同电极位置，预测在不同刺激强度（-0.5 V 至 -5.0 V）下四个接触点对总纤维束的激活情况。<a href="https://www.jove.com/files/ftp_upload/57292/57292fig1large.jpg" target="_blank">请点击此处查看该图的放大版本。</a>
<img alt="不同电压下深部脑刺激电极激活与纤维激活的对比示意图。" class="xfigimg" src="/files/ftp_upload/57292/57292fig2.jpg" data-altupdatedat="1788250015000" data-alt="Figure 2"> 
图2：轴对称与定向深部脑刺激（DBS）电极对丘脑中央部纤维束激活预测的比较。 (A) Medtronic 3387 圆柱形 DBS 电极与 Sapiens 定向 DBS 电极的可视化图示，显示在单一单极刺激幅度下，目标丘脑中央部纤维束的激活函数计算结果。(B)(i) 穿过 DBS 电极与纤维束的横截面切片。(ii) 横截面区域的二维投影，标示出 DBS 电极、目标纤维以及需避免激活的纤维。(C) 随着刺激幅度增加，轴对称与定向 DBS 电极在目标纤维及需规避纤维上的激活扩展情况。特定刺激幅度下被激活的纤维以红色表示，未被激活的纤维以蓝色表示。(D) 汇总面板 C 中所有模拟结果，显示纤维束各区域的预测刺激阈值幅度。<a href="https://www.jove.com/files/ftp_upload/57292/57292fig2large.jpg">请点击此处查看该图的放大版本。</a>

本项目的目标是开发一个交互式、针对患者个体的建模流程，以近实时地模拟深部脑刺激的效果，并提供有意义的反馈，说明这些设备如何影响大脑中的神经活动。

脑深部电刺激（DBS）是一种通过植入电极对特定脑区进行刺激的成熟疗法，用于……

本项目的主要目标是开发一个交互式、患者特异性的建模流程，以近实时地模拟深部脑刺激对神经纤维通路的影响。深部脑刺激是一种已获美国食品药品监督管理局（FDA）批准的治疗干预手段，十余年来一直用于帕金森病等运动障碍类疾病的治疗，并已在多种其他神经及精神疾病中得到研究。深部脑刺激的计算建模作为一种研究工具，已被证明具有重要价值，其应用日益广泛，有助于更深入地理解并预测个体受试者中刺激所产生的效应。计算建模面临的一个长期问题是构建和运行计算模型实验所需的时间过长。本实验方案的目的是介绍几种实现近实时且更具交互性的计算建模方法，这些方法通过若干创新得以实现。其中最主要的一项创新是能够对个体受试者的完整头部进行交互式网格划分。本视频中的方案将介绍如何从获取的 MRI 图像数据构建患者特异性的有限元模型，并模拟由植入式深部脑刺激电极诱导的电场。创建这些模型的主要步骤包括：第一，构建代表患者大脑及植入电极的有限元模型；第二，在有限元网格中添加脑组织的生物物理特性以及深部脑刺激设备的刺激参数；第三，求解电极在模型中产生的电压分布。本项目中建模的电极为美敦力 3387 型 DBS 电极。然而，任何电极几何形状均可与此方法配合使用。3387 是一种圆柱形导联，带有四个由非导电材料隔开的触点。在网格划分过程中，将为这些组件中的每一个创建表面表示。为四个触点、五个轴段、封装层以及两个用于控制电极周围网格密度的外围框创建了封闭曲面。创建触点的高密度表面网格非常重要，以保持圆柱几何形状并准确求解生物电场问题。通过求解泊松方程来计算电极在脑组织中诱导的电压分布，借助有限元建模方法，该方程可简化为线性方程组 Ax = B，其中 A 为表示网格电导率和几何特性的刚度矩阵，x 为网格各节点上的电压解，B 则根据边界条件和电流源进行调整。根据屏幕上表格中显示的材料特性，必须为封装层、轴段以及四个触点分配相应的电导率值。采用固定边界条件来计算生物电场。在四个DBS触点中每一个的中心位置设置点表面，可通过任意组合施加单极或双极刺激。将大脑外表面设为零伏，以模拟远处的回路电极。这些条件代表了有限元模型中已知的电压。我们已知四个触点以及回路电极（即植入式脉冲发生器）处的电压。采用共轭梯度法求解线性方程组，以计算模型在一组边界条件下的电压解。SCIRun 是一个先进的科学问题求解环境，用于对多类科学问题进行建模、仿真和可视化。其可视化编程语法允许用户将软件模型相互连接，以处理和可视化科学数据。本方案使用 SCIRun 创建大脑及深部脑刺激电极的有限元网格，并在交互式环境中求解生物电场问题。SCIRun 的一个关键特性是能够在处理流程的任意阶段对数据进行交互式可视化，这一点我们将在下文进行演示。构建有限元模型的第一步是确定大脑的结构解剖。为了执行此步骤，我们选择使用 FreeSurfer，这是一款用于处理和分析人类大脑 MRI 图像数据的开源软件套件。首先，将患者的 T1 加权 MRI 以 DICOM 数据序列的形式输入 FreeSurfer，并运行以下命令。接着，使用以下命令启动自动化的处理与分割流程。FreeSurfer 会创建一个去除颅骨后的新脑组织图像体数据。对该新体数据进行强度阈值分割和等值面提取，以获得脑组织的表面表示，该表面表示随后用于定义有限元模型的外边界。FreeSurfer 流程还可对白质、灰质、皮层下核团以及脑室进行分割。每个脑区被分配一个唯一的整数标签。下一步预处理是从扩散加权成像数据中估算扩散张量，并将该体积数据与T1加权磁共振成像进行共配准。扩散加权成像（DWI）是一种用于测量水分子在脑组织中多个方向上扩散性的技术。在这一步中，我们选择使用 3DSlicer，这是一个用于医学图像处理和成像数据三维可视化的免费开源软件平台。首先，在 Slicer 中使用 DWI 转换模块将 DICOM 序列转换为单个 DWI 体积数据。为了在张量重建过程中消除背景和噪声，需使用扩散加权体积掩膜模块创建掩膜。这将使得张量仅在脑组织内进行计算。下一步是使用DWI转DTI估计模块来生成扩散张量。将张量估计算法设置为最小二乘法，并勾选“校正负特征值”。这将去除与噪声或采集误差相关的具有负特征值的张量。电极模型的三维表面几何结构可以使用任何三维建模软件构建，本视频中不会演示该过程。此处我们为每个电极触点、轴段、封装层以及用于控制网格密度的两个周围框体创建了封闭的表面网格。将第一部分创建的脑表面与电极几何结构结合，使用网格生成工具 TetGen 构建包含植入电极的脑部有限元模型。将脑表面和电极几何结构输入至 InterfaceWithTetGen 模块中。根据屏幕上表格所示，对模型中每个区域的四面体单元生成施加体积约束。该技术在电极-组织界面附近（此处电场梯度较高）生成具有自适应尺寸且网格密度更高的网格。结果得到一个嵌入脑组织中的电极模型的四面体有限元网格。SCIRun 可通过几何变换模块和查看器窗口中的交互式控件实现对象的移动。要改变电极几何结构的角度，请结合使用 CreateGeometricTransform 模块和 TransformMesh 模块。将 X 旋转轴设为 1，并设置 theta 的旋转角度以控制角度。将 GenerateSinglePointProbe 模块直接连接到 ViewScene，以创建一个球形控件，同时连接到另一个 TransformMesh 模块以控制电极尖端的位置。在视图场景中，按住 Shift 键并单击拖动控件，以移动电极尖端的位置。这些变换会在电极几何结构传递给 TetGen 之前应用，因此更新后的电极位置和方向将被整合到有限元模型中。下一步是为有限元模型分配电导率值。加载在 3DSlicer 中处理过的 DTI 数据，并将电极模型和脑组织的电导率值赋给上一节中创建的四面体网格。将电导率值应用到网格后，将场数据输入 BuildFEMatrix。这是建立生物电场计算方程组的第一步。对于电压控制刺激，我们现在将在有限元网格上设置已知的电压边界条件。将用作网格外边界的脑表面设置为零伏特，以模拟远处的回路电极。然后使用 InsertVoltageSource 在网格上施加此边界条件。对于单极刺激，在底部电极中心位置的一个点上分配一个电压值。再次使用 InsertVoltageSource 将此边界条件应用于网格。将 BuildFEMatrix 和 InsertVoltageSource 的输出结果合并后输入至 SolveLinearSystem。使用共轭梯度求解方法并执行模块以计算电压解。利用 SetFieldData 模块将方程组的解重新赋值到有限元网格的节点上。电压解可直接显示在网格节点上，但无法直观反映其在脑组织中的分布情况。更直观的可视化方法是创建类似于地图上等高线的等电位面，这也能提供关于电场方向和强度的信息。将有限元网格的电压解传递给 ExtractSimpleIsosurface，并将数量设置为 10。这将在电极周围生成 10 个等间距的等电位面。利用该技术，我们可以直观地观察将各向异性电导张量引入模型后所产生的效应，以及该参数对脑组织内电压分布的影响。将有限元构建、生物电场计算和结果可视化整合到同一环境中，能够近乎实时地反馈电极位置和刺激强度等参数变化对生物电场与周围组织相互作用的影响。在调整电极位置并重新计算生物电场后，将所得解映射到目标纤维束上。通过计算激活函数来预测在给定刺激强度下纤维束的激活情况。我们收集所有触点及多种刺激强度下的纤维束激活预测结果，快速确定在特定电极位置下激活目标纤维最有效的刺激配置。随后，在脑深部电刺激（DBS）手术规划阶段，我们比较多个电极位置的激活模式，以确定能够有效刺激目标纤维通路的电极植入路径。在此可以看到，与位置一相比，位置二在所有触点上的刺激幅度更低，却能更有效地激活目标纤维束。该方法可用于模拟更为复杂的场景。该系统的一个新特性是能够快速切换模型中的电极几何形状，并能够模拟多个彼此靠近的电极导线。本研究旨在探索用于刺激参与觉醒调节的中央丘脑纤维通路的铅置入位置，以期为创伤性脑损伤患者提供潜在治疗靶点，同时避开感觉核团等邻近区域。针对固定的电极位置，我们模拟了改用定向电极模型后产生的效应，以实现对纤维激活位置更精确的控制，并避免刺激其他纤维束。通过观察电极周围横截面区域，可见定向电极能够有效刺激目标纤维，同时将激活信号扩散至需避免的纤维中的程度降至最低。该方法还可用于丘脑底核靶向治疗，以缓解帕金森病症状，同时避免刺激邻近的内囊。在此，我们模拟植入三根独立的深部脑刺激电极，并预测在非人灵长类动物实验的手术规划中，类似的中央丘脑纤维通路的激活情况。使用多个电极的优势在于能够对更大范围的区域进行刺激控制。方向性电极和多通道导线均增加了确定电极刺激参数的难度。然而，采用本视频中介绍的方法，能够更快速、更有针对性地探索复杂的参数空间，从而确定具有治疗效果的刺激参数设置。我们今天向您展示的患者特异性建模流程，涵盖了从患者影像数据到交互式模型的各个步骤，将有助于弥合工程师与临床医生之间的差距。具有超过四个触点的新型脑深部刺激（DBS）设备的发展，包括插入多个电极，为一个本已复杂的问题增添了更多难度。能够快速可视化不同电极设计在刺激配置下的作用效果，可为深入理解这些设备在大脑中的工作方式提供有价值的反馈信息。

targeting-neuronal-fiber-tracts-for-deep-brain-stimulation-therapy

Research

科研

JoVE Journal

Bioengineering

生物工程

应用交互式、患者特异性模型对神经纤维束进行深部脑刺激治疗的靶向性

深部脑刺激（DBS）手术治疗作用背后的机制仍需进一步研究。本文所述方法介绍了一种实验性方法，通过分析可能促进DBS术后神经发生的候选基因的基因表达谱，来探讨DBS所触发的细胞事件。...

丘脑前核深部脑刺激后大鼠海马区基因表达变化的分析

本文概述了以磁共振成像引导的重复经颅磁刺激作用于背内侧前额叶皮层，作为重度抑郁症实验性治疗方法的操作步骤。...

MRI引导的背内侧前额叶皮层重复经颅磁刺激治疗难治性重度抑郁症

本方案展示了一种家庭多模态数据采集平台的原型，该平台支持针对神经系统运动障碍患者优化自适应深部脑刺激（aDBS）的研究。我们还展示了将该平台部署于一名帕金森病患者家中超过一年所获得的关键研究结果。...

将诊所带回家：一种支持自适应深部脑刺激的家庭多模态数据采集生态系统

本方案描述了在功能磁共振成像（fMRI）期间进行局灶性经颅直流电刺激（tDCS）时，采用神经导航电极定位的方法。...

在功能磁共振成像期间施用的神经导航聚焦经颅直流电刺激

一种结合多种磁共振成像（MRI）数据模态指导下的个性化经颅磁刺激与脑电图（TMS-EEG）皮层网络定位的新方案。...

神经影像引导的经颅磁刺激–脑电图用于实时皮层网络定位

本文介绍了在缺乏导航设备的情况下，如何定位重复经颅磁刺激干预或治疗的功能特异性靶点。...

在缺乏导航设备的情况下定位经颅磁刺激的功能特异性靶点

本文介绍了一种使用重复经颅磁刺激（nTMS）结合基于扩散张量成像（DTI）的皮质脊髓束（CST）重建进行运动区定位的标准化操作流程。该方案具有可重复性，临床可行性高，可轻松整合到常规临床工作流程中，为运动通路评估、神经可塑性研究以及康复方案制定提供了可靠且有价值的技术框架。...

使用导航经颅磁刺激进行功能性运动定位的标准化方案

本方案描述了机器人辅助的fMRI引导个性化经颅磁刺激治疗难治性抑郁症的工作流程。采用本方案进行临床实施是可行的，开放标签研究结果支持该方法在真实世界中的有效性。...

机器人递送的fMRI引导个性化经颅磁刺激治疗难治性抑郁症

在此，我们描述了一种实验方案 离体 从接受治疗的肌张力障碍仓鼠模型的急性小脑脑片中进行高密度微电极阵列记录 体内...

利用高密度微电极阵列研究苍白球深部脑刺激的全网络作用机制

在许多生物学和临床情境中，研究细胞过程在其天然微环境中的动态演变具有显著优势。本文介绍了低成本光纤显微镜的组装及其应用，该显微镜可实现实时成像，适用于细胞培养、动物实验以及临床患者研究。...

高分辨率光纤微内窥镜技术用于 原位 细胞成像

了解如何使用交互式患者特异性模型模拟深部脑刺激对神经纤维束的影响。该方法利用SCIRun、FreeSurfer和3DSlicer软件实现实时计算建模，为帕金森病等运动障碍疾病的个性化治疗策略提供重要支持。

利用交互式患者特异性模型为深部脑刺激疗法定位神经纤维束

探索更多视频