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

摘要

该项目的目标是开发一个交互式的, 病人特定的建模管道, 以模拟的影响, 深脑刺激在近实时, 并提供有意义的反馈, 这些设备如何影响神经活动的大脑。

摘要

深部脑刺激 (DBS), 其中包括插入电极, 以提供刺激到局部大脑区域, 是一个既定的治疗运动紊乱, 并正在应用于越来越多的疾病。计算建模已成功地用于预测 DBS 的临床效果;然而, 需要新的建模技术来跟上星展设备日益增长的复杂性。这些模型还需要快速准确地生成预测。该项目的目标是开发一个图像处理管道, 将结构磁共振成像 (MRI) 和扩散加权成像 (DWI) 纳入一个交互式的, 病人的具体模型, 以模拟的影响, DBS。虚拟星展导线可以放置在病人模型内, 以及主动接触和刺激设置, 在那里铅位置或方向的变化产生一个新的有限元网格和解决生物场问题近实时,timespan 大约10秒。该系统还能够在接近接近度的情况下模拟多个引线, 以便在不同的引线上通过不同的阳极和阴极对电流进行转向。本文提出的技术减少了生成和使用计算模型的负担, 同时为研究人员或临床医生提供了对电极位置、电极设计和刺激配置的影响的有意义的反馈。谁可能不会建模专家。

引言

深部脑刺激 (DBS) 是一种已建立的治疗运动障碍, 如必要的震颤¹和帕金森病²。这项治疗也被调查为一个潜在的治疗越来越多的疾病, 包括创伤性脑损伤³, 抽动秽语综合征⁴, 抑郁症⁵。星展系统需要外科植入电极导致在局部脑区进行刺激, 以调节正在进行的神经活动⁶。电极的位置和刺激参数都对调节神经元回路产生治疗效果有一定的影响。小的位置变化可能会影响治疗窗口, 可能会增加不良副作用的可能性, 在治疗效果达到^7,8,9。在实践中, 通常很难预测刺激对神经活动的影响;因此, 这一治疗效果的窗口是在病人的基础上确定的, 因为刺激装置是由临床医生^8,9编程。随着新一代星展设备的推出, 这个过程变得更加复杂。例如, 新的引线设计正在引入更多的联系人^10,11,12, 在某些情况下, 多个潜在顾客被植入近距离¹³。因此, 有必要能够探索和预测 DBS 在一个大的和不断增长的参数空间的影响。

计算建模和分析可用于预测星展集团在患者特定基础上的生理和临床效果。这些模型使用有限元建模 (FEM) 建立三维表示的脑组织和生物物理特性的植入电极。有限元生物场模型已成功地用于预测星展¹⁴的影响, 但到目前为止, 这些都是耗时和计算成本高昂的产生。需要新的建模技术来跟上星展设备日益增长的复杂性。这些病人特定的模型应提供近实时视觉反馈的影响, 星展集团的潜在位置或刺激参数的变化。用户将在几秒钟内获得有关潜在顾客位置和刺激设置的反馈, 从而在几分钟内继续改进铅的放置过程。病人特异性是通过结合病人的解剖, 他们的大脑的形状和大小, 当建立有限元和应用他们的大脑的生物物理特性, 如各向异性组织电导率。各向异性电导率描述了电流将如何通过不同的大脑区域传播, 并且可以非侵入性测量整个大脑类似于典型的磁共振图像 (MRI)。

星展的建模方法不使用病人特定的信息可以提供快速, 但不准确的预测刺激效应, 由于广义几何和电导率的大脑组织。在这种方法中, 对所有患者使用单一的有限元, 预测的神经活动可以提前计算。由于每个个体都建立了一个新的有限元模型, 因此不能对特定的病人进行泛化和预计算。这些模型需要更多的努力来构建, 但可以更准确。有几个因素限制了这些模型的生成和使用速度: 1) 在模型生成管线开始时修改参数, 如电极位置, 需要手动更新所有后续步骤;和 2) 建模管道中的步骤不易相互集成, 需要在多个软件包之间传递数据。通常, 我们要评估许多不同的情况, 如电极位置, 刺激参数, 或电极设计。为了提供有意义的反馈, 这些变化对患者将接受的治疗效果的影响, 这些结果应该是准确的, 并迅速产生。

我们的目标是提供新的技术, 以建立病人特定的模型, 利用通用模型获得的速度和自动化许多管道步骤, 以创建一个交互式建模环境, 提供近实时视觉对 DBS 效应的反馈。交互式模拟允许用户在不关注模型构造细节的情况下快速测试预测并获得结果。这是有益的, 当有一个大的参数空间来探索, 这些参数如何影响模拟是不确定的。我们将描述在处理管道中的步骤, 以产生交互式的, 病人特定的有限元模型的磁共振成像 (MRI) 收购。利用本文概述的工具和技术, 可以减少创建有限元生物场模型的时间成本, 并提供一种方法, 使这些模型可供没有建模专家的研究人员和临床医生使用。

该协议描述了如何从获得的 MRI 容量中构建病人特定的有限元模型, 然后模拟星展电极引起的电场。创建这些模型的主要步骤是: 1) 建立一个代表患者大脑和植入电极的有限元模型 (FEM), 2) 增加大脑的生物物理特性, 并将星展设备的刺激参数限定为有限元, 3) 解决在模型中由电极产生的电压。需要两种成像方式来建立星展模拟的病人特定模型。T1 MRI 用于构建脑表面、心室和特定细胞核的分割。扩散加权成像 (DWI) 是一种水扩散率测量方法, 用于估计整个脑组织¹⁵的扩散张量。将扩散张量转换成电导率张数, 以体素为基础, 将组织的非均匀、各向异性生物物理特性量化为¹⁶。通过求解泊松方程, 通过有限元法简化到线性方程组的线性系统, 计算出电极所引起的整个脑部的电压分布, 其中a是刚度矩阵,表示网格的电导率和几何, x是网格中每个节点的电压解, b根据边界条件和电流源进行修改。

研究方案

1. 图像处理

1. T1 MRI 分割
   1. 下载并安装 FreeSurfer¹⁷。假设 T1 MRI 已获得, 并以 DICOM 或 NIFTI 格式。
   2. 键入以下命令以创建病人目录并在 FreeSurfer 中添加他们的 T1 卷:侦察全 patientName-i/全/路径/到/nii
   3. 键入以下命令以运行 FreeSurfer 的自动分割:侦察全 s patientName-所有
   4. 使用mri_convert将 "aseg.auto.mgz" 转换为 FreeSurfer 的专有 mgz 格式。
2. 扩散加权成像 (DWI)
   1. 下载并安装切片器¹⁸。假定已执行了一个 DWI 捕获, 并可作为 DICOM 文件的一叠。
      注意: 执行此重建的一个好的命令行工具是 FMRIB 软件库 (FSL)¹⁹。在该协议中, FSL 用于在重建前执行回波平面、运动和涡流畸变校正。该分析中使用的 DWI 数据是一系列的 DICOM 文件, 用于41种不同的扩散敏感梯度方向。
   2. 使用切片器中的DWIConverter模块将41个单独的收购重构为单个卷。指定在其中的 DWI 数据所在的输入 Dicom 数据目录。选择DicomToNrrd单选按钮并提供输出 DWI 音量的名称。
   3. 用扩散加权体积掩蔽模块创建张量掩码, 以消除 DTI 估计的背景噪声。将阈值参数保留为默认为0.5。
   4. 使用在上一步中创建的掩码作为扩散张量掩码, 将此单一的 dwi 体积转换为扩散张量 (dti) 体积, 并将dwi 估计为 dti 。将估计参数设置为最小二乘法, 并确保检查了移位负特征值选项。
   5. 以 NRRD 格式保存新创建的 DTI 卷。

2. 有限元模型生成

1. SCIRun
   1. 下载并安装 SCIRun 版本 5 (http://www.sci.utah.edu/software/scirun.html)。SCIRun 是一个解决问题的环境, 用于建模、模拟和可视化科学问题。这是用于创建、解决和与此项目中开发的模型进行交互的计算工作台。
      注意: SCIRun 网络是由执行特定任务的各个模块的连接序列构建的。通过由简单模块组成的网络传递数据, 可以进行更复杂的处理和模拟。
2. 生成电极几何
   注: 在此项目中建模的电极是美敦力3387星展²⁰铅。它是一个圆柱形引线, 四个接触1.5 毫米的高度, 直径1.27 毫米, 间隔1.5 毫米分开。非导电材料将四接触点隔开。
   1. 为星展公司潜在顾客中的每个组件创建一个闭合曲面网格, 以匹配列出的维度。使用两个原始几何图形、圆柱和球体创建整个引线。
      注意: 可以在各种3D 建模程序中创建引线几何图形。
   2. 创建一个表面网格的封装层0.5 毫米厚周围的整个铅。
3. 创建整个脑网
   1. 加载1部分中创建的脑表面, 作为有限元的外边界。
   2. 创建两个同心盒表面围绕星展集团的领导, 以控制网格密度围绕电极。
      注: 星展导线和盒面应沿正 z 轴与位于原点的电极轴的底部对齐。这对于以后电极的旋转和平移是很重要的。
   3. 创建一个点云, 其位于网格的每个不同区域中的单个点。网格的区域有: 轴段、四触点、盒子和脑表面。附加到这些点的值表示该区域网格中每个四面体元素的最大卷约束。这些约束调整, 以控制网格密度在和周围的电极。
   4. 使用 "InterfaceWithTetGen" 模块创建整个脑网。按如下方式设置输入端口:
      端口 1: 脑表面 (从 1.1)
      端口 2: 区域属性点云, 卷约束 (从 2.3.3)
      端口 3: 无
      端口 4: 电极模型表面 (从 2.2)
4. 交互式电极放置
   注意: 这一步的目标是使电极在大脑内的任何地方都能轻松移动, 并让系统自动更新新的网格和模拟。
   1. 使用两个级联CreateGeometricTransform模块在 X 和 Y 平面上旋转电极。单击 "旋转", 并将第一个模块设置为1.00 轴Y , 第二个到 1.00轴 X。旋转角度设置为底部滑块,旋转θ (度数)。
   2. 为需要转换的每个转换和对象连接一个TransformMeshWithTransform模块。
   3. 创建一个GenerateSinglePointProbeFromField模块, 并将其几何端口连接到ViewScene模块。通过按住 shift 键和左键单击, 在ViewScene窗口中移动球体。点探头模块在释放时更新球体的位置。
      注意: 此位置用于确定应用于电极几何的平移变换。
   4. 使用TransformMeshWithTransform模块将电极几何的每个部分转换为点探头位置。
   5. 将两个旋转变换和从点探针构件到除脑表面以外的每个物体的平移。
      注: 所有三的转换都应在每个对象上执行, 然后再管道进入InterfaceWithTetGen模块, 这样啮合算法就知道在新的位置和方向上将大脑内的电极几何网格化。每次在大脑中移动电极时, 都会生成一个新的网格。

3. 生物场计算

1. 设置电导率参数
   1. 加载在步骤1中创建的 DTI 卷, 并将数据从步骤2中创建的InterfaceWithTetGen的四面体网格输出设置为使用线性 ("加权") 选项的MapFieldDataFromSourceToDestination 。
      注意: 任何各向异性电导率估计都可以在这个步骤中实现。
   2. 使用CalculateFieldData将轴区域电导率设置为 1e-6, 并将触点区域设为1e6。一旦设置了所有的电导率值, 将该字段管入BuildFEMatrix。
2. 设置电流/电压源和边界条件
   1. 管在啮合过程中使用的大脑表面从步骤2到SetFieldDataToConstantValue , 并将该值设置为0。这将在整个大脑表面产生一个电压接收器。然后使用InsertVoltageSource将电压值放在InterfaceWithTetGen的四面体网格输出上。
   2. 对于单刺激, 通过计算一个接触面上的ReportFieldInfo模块中的 FieldCenter 值, 在一个联系位置的中心创建一个点。将此点上的字段数据设置为-1 伏。然后使用InsertVoltageSource将点源放置在四面体网格中。
      注: 现在有两个四面体网格的副本。一个有电压值设置在大脑表面和一个点源内的联系人。
   3. 加入两个InsertVoltageSource模块的输出, 并将输出ApplyFEMVoltageSource输出BuildFEMatrix。
3. 解决正向问题
   注: ApplyFEMVoltageSource的输出给出端口1上的刚度矩阵和端口2上的右端矩阵。这是解决方程组所需的两个组成部分。
   1. 将ApplyFEMVoltageSource的两个输出管成SolveLinearSystem。选取共轭梯度解法和雅比法作为预设条件。本模块的输出是求解矩阵, 在四面体网格中每个节点的电压值。
   2. 使用SetFieldData将电压解矩阵放到四面体网格上, 以3D 表示电压分布。
4. 等值面可视化
   1. 把SetFieldData的结果管成ExtractSimpleIsosurface。将 UI 打开到 "数量" 选项卡, 并将平均间隔的 isovalues 数设置为10。
   2. 使用ShowField模块并将其连接到ViewScene模块, 以显示等值面为透明, 并将彩虹色图缩放到电压解决方案的最小和最大值。

结果

在本议定书的末尾提供了创建病人特定模型所需的所有输入: 脑表面、电极几何和电导率张量。还应该创建一个 SCIRun 网络, 将所有输入集成起来, 建立一个有限元模型, 并模拟诱导的生物场。该模拟模型允许在大脑中的电极运动和参数的修改, 如主动接触或刺激幅度。

图 1演示了在多个电极位置模拟和比较光纤激活预测的能力。在该协议中没有显示出从扩散成像中产生的光纤束, 但这种分析可以用任何确定性纤维束成像方法进行。位置2预测, 与位置1相比, 在所有触点上, 目标纤维束在较低的刺激振幅下更健壮地活化。该分析用于星展手术的规划阶段, 确定电极轨迹, 有效地刺激靶纤维通路。

该系统的一个新特点是能够快速切换模型中的电极几何, 并能够在接近彼此的情况下模拟多个引线。这两种方法都用于对电极周围电场的形状和转向提供更大的控制。图 2显示了轴对称和定向星展集团铅之间的纤维活化预测的比较, 同时研究了铅的位置, 以刺激中心丘脑纤维通路和避免附近地区, 如感官核.对于相同的电极轨迹, 我们可以看到定向引线能够引导纤维活化到目标区域, 同时避免刺激不受欢迎的纤维。这种分析可以用来确定什么时候需要改变的设备, 病人将接受星展治疗, 以避免刺激附近地区, 会诱发不必要的副作用。这种方法有助于靶向丘脑核治疗帕金森病的症状, 同时避免附近的内囊。

图 1: 电极位置的比较和目标纤维活化的预测.预测刺激振幅的总纤维束活化 (-0.5 v 到-5.0 v) 和所有四触点为两个不同的电极位置。请单击此处查看此图的较大版本.

图 2: 中心丘脑纤维束活化预测与轴对称和定向星展集团线索的比较.(A)可视化美敦力3387圆柱星展集团铅和智人定向星展集团铅与激活功能计算的目标中央丘脑纤维束为一个单一的单刺激幅度。(B)(i)通过星展集团的铅和纤维束进行横断面切片。(ii)跨截面积的二维投影, 以识别星展集团的铅、靶纤维和纤维以避免。(C)在两个目标纤维和避免纤维之间的活化扩散, 以刺激振幅增加为轴对称和定向星展集团铅。给定刺激振幅的活化纤维以红色显示, 而未激活的纤维以蓝色显示。(D)汇编 C 组中显示的所有模拟, 显示了在光纤束的所有区域内预测的刺激阈值。请单击此处查看此图的较大版本.

讨论

该协议展示了降低星展集团患者特定模型的时间负担和计算成本的技术, 使其能够在一定程度上实现模拟结果的近实时反馈。快速反馈使探索一个大的参数空间, 以更好地了解参数变化如何影响模型结果。这些参数包括引线位置、主动触点的选择、振幅、脉冲宽度和激发波形在这些接触处的频率。所建议的工具的主要特点是: 1) 一个简单的用户界面, 以近实时可视化的方式来调整模型参数, 这些参数如何影响仿真, 以及 2) 从一小部分输入中创建模型的自动化: 大脑表面,脑组织电导率张量和电极几何的表面表示法。这种自动化加速了为许多人创建病人特定的模型, 他们有不同的大脑几何和组织电导率, 以及评估插入不同的电极设计到现有的模型的效果。本协议中描述的图像预处理步骤不是完全自动化的, 可能需要一天的处理时间。但是, 一旦完成了从这些步骤生成的数据被认为是静态的, 这意味着在模拟过程中不会修改这些数据。模型创建的自动化是由系统在没有人工努力的情况下将这些数据应用于特定于病人的有限元的能力。SCIRun 网络执行所有的模型生成、仿真、可视化步骤只需要构建一次。因此, 只有图像预处理步骤需要再次执行, 才能为传入的病人生成特定于病人的模型。

建模管线产生结果的性能增益是由网格生成、生物场计算和解的可视化集成到单个软件环境中的。采用自适应网格生成等现有建模技术, 在电极周围产生较高的网格密度, 使电极的密度更低, 从而缩短了有限元的生成和求解时间。该软件 SCIRun 还能实现网格生成和生物字段计算的自动化。使用交互式小部件的用户定义的电极移动会触发新的网格的构造, 并更新电极位置。这包括修改新电极位置的边界条件和电导率值。

电极几何被视为一个自由运动物体内的大脑体积之前, 其位置被纳入有限元。这种网格构造方法的一个关键含义是可以很容易地将多个电极插入模型中。例如, 电极几何的第二个拷贝可以放置几毫米远, 并且两者将包括在有限元中。在最近的研究中, 两个电极被植入近距离治疗多发性硬化震颤¹³ , 并已用于非人类灵长类实验, 以探索有效的刺激目标²¹。使用多电极的好处是更好地控制在较大区域内组织产生的电场。刺激与活跃接触在两个电极可以引导电流向目标区域和远离区域, 将导致负面的副作用。更精细地控制刺激在一个更大的区域也是有用的探索不同的刺激地点, 当确切的地点的预期目标是未知的, 像这样的情况下, 许多新兴疗法应用的 DBS。然而, 由于已有较大的参数空间的增加, 确定获得治疗性刺激的参数比单个电极更具挑战性。

我们设想, 这种交互式建模工具可以为星展公司的植入前规划提供好处。对脑组织中刺激程度的反馈可以让外科医生在其手术计划中修改电极位置, 以便为靶区提供治疗性刺激。星展疗法是开发这一工具的主要动力因素, 本文所提出的技术可以应用于任何具有不同刺激或记录范式的有限元生物场模型。刺激疗法, 如经颅直接电流刺激抑郁症²²或使用深度电极治疗癫痫²³与 DBS 在确定最佳刺激位置, 以实现治疗结果。皮层脑电图, 一种记录技术, 其表面的电极阵列, 以确定癫痫发作的地区, 有一个挑战, 确定在哪里放置电极记录从目标地区的大脑²⁴。所有这些应用都依赖于电极位置, 同时处理电流在脑组织中流动的不确定性。本文所提供的技术减少了生成和使用计算模型的负担, 同时为研究人员和临床医生提供了有意义的反馈, 这些设备不是模型专家。

材料

Name	Company	Catalog Number	Comments
FreeSurfer	Athinoula A. Martinos Center for Biomedical Imaging		https://surfer.nmr.mgh.harvard.edu/
3D Slicer	BWH, Harvard University		https://www.slicer.org/
SCIRun	University of Utah Center for Integrative Biomedical Computing		http://www.sci.utah.edu/cibc-software/scirun.html