登录

需要订阅 JoVE 才能查看此. 登录或开始免费试用。

本文内容

  • 摘要
  • 摘要
  • 引言
  • 研究方案
  • 结果
  • 讨论
  • 披露声明
  • 致谢
  • 材料
  • 参考文献
  • 转载和许可

摘要

该协议演示了 SILICOFCM 平台的工作流程,通过应用心脏的多尺度机电模型,从患者特定的超声图像自动生成左心室的参数模型。该平台支持 计算机临床试验, 旨在减少实际临床试验并最大限度地提高积极的治疗效果。

摘要

SILICOFCM项目主要旨在开发一个用于家族性心肌病(FCM) 的计算机临床试验的 计算平台。该平台的独特之处在于整合了患者特定的生物学、遗传学和临床成像数据。该平台允许测试和优化医疗,以最大限度地提高积极的治疗效果。因此,可以避免不良反应和药物相互作用,可以预防心源性猝死,并且可以缩短药物治疗开始与预期结果之间的时间。本文介绍了通过应用心脏的机电模型从患者特定的超声图像自动生成的左心室的参数模型。通过入口和出口流量的特定边界条件、心电图测量和心肌特性的钙功能来规定药物效应。患者的遗传数据通过心室壁的物质特性被纳入。心尖视图分析涉及使用先前训练的U-net框架分割左心室,并根据舒张期和收缩期中左心室的长度计算边界矩形。M 模式视图分析包括在 M 模式视图中左心室特征区域的边界。提取左心室尺寸后,根据网格选项生成有限元网格,并使用用户提供的入口和出口速度运行有限元分析模拟。用户可以在平台上直接可视化各种模拟结果,例如压力-体积图、压力-应变图和心肌工作时间图,以及位移、压力、速度和剪切应力等不同领域的动画。

引言

近年来信息技术、仿真软件包和医疗设备的快速发展为收集大量临床信息提供了机会。因此,创建全面而详细的计算工具对于从丰富的可用数据中处理特定信息变得至关重要。

从医生的角度来看,区分特定患者的"正常"和"异常"表型对于估计疾病进展、治疗反应和未来风险至关重要。最近的计算模型显著提高了对肥厚型(HCM)和扩张型(DCM)心肌病中心肌行为的综合理解1。使用高分辨率、详细且解剖学上准确的全心电活动模型至关重要,这需要大量的计算时间、专用软件和超级计算机123。最近使用基于Holzapfel实验的线弹性和正交各向异性材料模型开发了真实3D心脏模型的方法,该方法可以准确预测心脏内的电信号传输和位移场4。开发新的综合建模方法可以成为区分多基因疾病患者症状类型和严重程度以及评估正常身体活动受损程度的有效工具。

然而,针对患者的建模存在许多新的挑战。人类心脏的物理和生物特性无法完全确定。非侵入性测量通常包括嘈杂的数据,从中很难估计个体患者的特定参数。大规模计算需要大量时间才能运行,而临床时间框架是有限的。患者个人数据的管理方式应使生成的元数据可以重复使用,而不会影响患者的机密性。尽管存在这些挑战,但多尺度心脏模型可以包含足够的细节水平,以实现紧密遵循观察到的瞬时反应的预测,从而为未来的临床应用提供希望。

然而,尽管多个研究实验室付出了大量的科学努力和大量的资助,目前,只有一个用于多尺度和全心脏模拟的商用软件包,称为 SIMULIA 活体心脏模型5。它包括动态机电仿真、精细的心脏几何形状、血流模型和完整的心脏组织表征,包括被动和主动特性、纤维性质和电通路。该模型旨在用于个性化医疗,但活性材料表征基于Guccione等人引入的现象学模型67。因此,SIMULIA 无法直接准确地翻译在许多心脏病中观察到的收缩蛋白功能特征的变化。这些变化是由分子和亚细胞水平的突变和其他异常引起的6。在临床实践中,SIMULIA 软件在少数应用中的有限使用是当今开发更高级别多尺度人心模型的一个很好的例子。另一方面,它推动了新一代多尺度程序包的开发,这些程序包可以从分子到器官规模追溯突变的影响。

心脏电生理学的主要目的是确定躯干内的信号传播和所有隔室456的特性。SILICOFCM8项目使用患者特定的生物学,遗传学和临床成像数据预测心肌病疾病的发展。它是通过对现实的肉瘤系统、患者的遗传特征、肌纤维方向、流固耦合和电生理耦合进行多尺度建模来实现的。左心室变形、二尖瓣运动和复杂血流动力学的影响给出了特定患者心脏病的详细功能行为。

本文演示了如何使用 SILICOFCM 平台使用具有机电耦合的流结构心脏模型,从患者特定的超声图像自动生成左心室 (LV) 的参数模型。利用深度学习算法生成LV的顶点视图和M模式视图分析。然后,利用网格发生器自动构建有限元模型,模拟左心室收缩9全周期的不同边界条件。在该平台上,用户可以直接可视化压力-体积、压力-应变和心肌工作时间图等仿真结果,以及位移、压力、速度和剪切应力等不同领域的动画。来自特定患者的输入参数是超声图像的几何形状、左心室输入和输出边界流动条件下的速度分布以及特定的药物治疗(例如,entresto、地高辛、mavacamten等)。

研究方案

本研究中的方案于2019年2月6日获得英国国家卫生服务健康研究局东北泰恩威尔南部研究伦理委员会的批准,参考编号为18 / NE / 0318,并被每个参与中心的机构审查委员会采用。该研究是在良好临床实践原则和赫尔辛基宣言之后进行的。从参与研究的所有受试者那里获得了知情同意。患者信息保持匿名。

1. 超声M模式或心尖视图DICOM图像分析和左心室参数提取的工作流程

注意:要开始此协议,用户必须登录到SILICOFCM平台8 并选择适当的工作流程(即使用M模式或顶点视图进行超声图像分析)。例如,超声M型DICOM图像分析和左心室参数提取的工作流程涉及几个步骤。第一步是模板匹配。要匹配的模板包括所有相关的必要边框。应该强调的是,对于每个新数据集,应为对应于特定超声机器的每个数据集手动提取一次模板。工作流图如图 1 所示。在Felsenszwalb高效的基于图形的图像分割算法找到对应于隔膜边界和LV壁的"强"边界后,将从分析的图像中提取与模板匹配的区域。 基于这些边界以及心脏直径最大(对应于舒张直径)和最小(对应于收缩期直径)的地方, 将计算各种LV尺寸。用户应定义视图是 M 模式还是顶视图。

  1. 使用用户名和密码登录平台。在 虚拟人口模块下,选择超声 M 模式或顶端视图工作流。
  2. 在可用工作流程列表中,选择 超声合并工作流程
  3. 在"文件上载"部分中,上载存储在用户计算机上本地的图像和DICOM 文件图2)。
    注意:来自特定患者的模块M模式视图DICOM图像分析和左心室参数提取的输入是测试图像(由用户以DICOM格式提供,代表M模式或顶端视图)和静态文件(图1)。静态文件包括 a) 输入.txt 文件,该文件是用于覆盖模块计算输出度量的模板文件(从该文件创建新的输出文件)和 b) 模板图像文件,用于模板匹配(即 template_dicom_GEMS.jpg)。
  4. 选择"专用"或"公用"文件夹作为文件 的目标文件夹 。在启动工作流之前,在"注释"部分中键入所需的注释或 注释 。选择先前上传的顶端视图和 M 模式左心室超声图像和 DICOM 文件(步骤 1.3)。(图3)。
    注意:超声心脏分割工具由两个子模块组成:心尖视图和M模式视图。这两个模块都可以通过平台访问。分析包括计算相应横截面中可见的特征参数。在这两种情况下,必要的值都会写入一个文件,该文件从用户端作为参数心脏模型的输入(有关详细信息,请参阅讨论)。
  5. 单击 "执行 "按钮。交互式平台会在正在运行的工作流完成时通知用户。
  6. 直接在平台上可视化LV创建的几何形状(屏幕上会自动出现3D模型,可以使用鼠标旋转)。可用选项包括用于可视化的 着色 线框 模型。
    注意:提取LV尺寸后,根据网格选项,生成有限元网格,并使用用户规定的入口和出口速度运行有限元分析模拟(见下文)。

2. 用于流固(FS)仿真的PAK有限元求解器

注意: 此工具可用于固-流体耦合问题的有限元分析。它支持固体和流体之间的强耦合和松散耦合。元素可以是六面体或四面体,元素中心有或没有附加节点。该求解器具有内置的材料模型,例如Holzapfel模型,猎人肌肉模型等。PAK-FS信息流程图如图 4所示。它从输入文件和 PAK 预处理器开始。PAK 预处理器工具输出一个 DAT 文件,该文件将用作有限元求解器的输入文件。求解器的最终输出是 VTK 文件,其中包含有限元仿真的结果:左心室的速度、压力、应变和应力。

  1. 使用底部的按钮下载入口和出口速度边界条件的模板文件(图3);如果患者特定的流边界条件可用,请下载并使用这些文件。通过单击底部的相应按钮下载网格选项。将这些文件保存到"专用"或"公用"文件夹中。
  2. 以与上传图像类似的方式上传这些文件(图 3)。规定的入口和出口速度模拟药物状况,而网格选项控制有限元网格的密度。要模拟患者特定的条件,请修改压力、流量、材料属性和钙函数的默认值。
  3. 单击 "执行 "按钮。列表中将显示一个新的正在运行的工作流。如果工作流的任何部分不清楚,请单击 "帮助文件 "按钮(图 3)以查看有关如何使用此工作流的详细说明,以及解释结果。
    注意:如果所有内容执行没有错误,则工作流的状态将从"正在运行"更改为"完成正常"。
  4. 使用多个选项查看结果。或者,下载结果;结果文件夹包含 VTK 文件、CSV 文件和动画。
    1. 单击眼睛按钮可查看射血分数和全局工作效率值,以及压力与体积、压力与应变和心肌功与时间的关系图。
    2. 单击相机按钮可预览和播放位移、压力、剪切应力和速度场的动画。
    3. 单击3D可视化按钮,在 ParaView Glance中在线可视化输出。
      1. 加载之前下载的多个 VTK 文件作为结果。查看几个感兴趣的参数,并将字段更改为可视化速度。
      2. 旋转模型或更改配色方案。选择" 带边的曲面 "或 "线框 "作为曲面的表示。对每个加载的 VTK 文件应用相同的方法。

3. 通过心电图测量确定心室激活序列

注意:实施了改进的心脏细胞FitzHugh-Nagumo模型。心前区导联线使用标准的六个电极进行建模。心脏的潜力通过逆心电图得到优化。从窦房结(这是时间的函数)的激活开始,通过心脏和躯干的异质动作电位,用户可以在全躯干模型上获得电活动。PAK-TORSO 信息流程图如图 5 所示。用户提供躯干模型在所有方向(x,y,z)上的缩放和ECG信号功能。然后创建缩放模型,并使用 PAK-FS 求解器模拟其行为。用户在文本文件中提供这些输入值。仿真的输出是一个VTK文件,其中包含躯干嵌入环境中心脏的电活动。

  1. 在主页上,转到执行工作流,然后在可用 工作流 列表中选择 torse-cwl 。在执行工作流之前,在"注释"部分添加注释或 注释
  2. 单击"输入 模板文件" 按钮,将网页上显示的内容保存为 输入.txt 文件,该文件将用于躯干模型。在 "输入文件 "字段中,选择下载的 输入.txt 文件。导入文件后,单击 "执行 "按钮开始计算。
    注意:新工作流将显示在左下角,其状态显示为"正在运行";计算持续约0.5小时,然后状态将更改为"完成正常"。
  3. 单击左下角的眼睛或相机按钮,直接在平台上可视化可用的模拟报告或动画。
    注意:结果包括每个时间步的电场、速度、压力、剪切应力和变形的分布。可以查看每个分布的动画,以显示整个心动周期。
  4. 或者,单击3D可视化按钮,在 ParaView Glance中在线可视化输出。
    1. 选择 "打开文件 "按钮,转到" GIRDER "选项卡,输入用户凭据(如果出现提示),然后打开"专用"文件夹。
    2. 在下一页上,选择 工作流输出 文件夹,然后打开 torso-cwl 文件夹。打开列表中的第一个文件夹。
    3. 请参阅表示仿真结果的 VTK 文件列表。选择一个或多个文件,然后单击 "选择" 按钮以在 ParaView Glance中加载文件。
    4. 使用鼠标操作模型几何图形(即移动、旋转、放大或缩小等)。
    5. 为模型表示选择不同的选项,如下所示。
      1. 选择 "线框 "选项以查看躯干内部,躯干内有一个心。选择 "点 "选项以显示具有全心形网格的躯干模型的虚线表示。
      2. 调整 "磅 值"值以更改显示结果。调整" 不透明度 "值以查看躯干内部并在心形网格内显示结果。
      3. 单击" 颜色依据 "下拉菜单,然后选择所需的选项,例如电势。将默认色阶更改为列出的任何选项。

结果

例如,超声M型DICOM图像分析和左心室参数提取的工作流程如图 1所示。M模式和顶视图可以单独测试,也可以一个接一个地测试,具体取决于感兴趣的参数。如果一个接一个地测试,结果将附加到一个通用文件中(分别用于收缩期和舒张期)。如果只测试一个视图,则未知参数的值取自默认文件 输入.txt图 1)。SILICOFCM工具根据视图模式分析图?...

讨论

SILICOFCM项目是一个 计算机 临床试验平台,用于设计虚拟患者群体以进行风险预测,测试药物治疗的效果,并减少动物实验和人体临床试验。使用规定的入口/出口边界流动条件、钙功能和材料壁特性对药物治疗的效果进行测试。该平台将肉瘤水平的多尺度方法与全心性能和功能优化水平相结合,用于预测心肌病疾病进展期间患者特定疾病的风险。

在SILICOFCM

披露声明

作者没有利益冲突。

致谢

这项研究得到了欧盟地平线2020研究和创新计划的支持,该计划根据SILICOFCM 777204赠款协议和塞尔维亚共和国教育,科学和技术发展部通过第451-03-68 / 2022-14 / 200107号合同提供支持。本文仅反映作者的观点。欧盟委员会对可能使用本文所含信息的行为概不负责。

材料

NameCompanyCatalog NumberComments
SILICOFCM projectwww.silicofcm.euopen access for registered users

参考文献

  1. Gibbons Kroeker, C. A., Adeeb, S., Tyberg, J. V., Shrive, N. G. A 2D FE model of the heart demonstrates the role of the pericardium in ventricular deformation. American Journal of Physiology. 291 (5), 2229-2236 (2006).
  2. Pullan, A. J., Buist, M. L., Cheng, L. K. . Mathematically Modelling the Electrical Activity of the Heart - From Cell To Body Surface and Back Again. , (2005).
  3. Trudel, M. -. C., Dub´e, B., Potse, M., Gulrajani, R. M., Leon, L. J. Simulation of QRST integral maps with a membrane based computer heart model employing parallel processing. IEEE Transactions on Biomedical Engineering. 51 (8), 1319-1329 (2004).
  4. Kojic, M., et al. Smeared multiscale finite element models for mass transport and electrophysiology coupled to muscle mechanics. Frontiers in Bioengineering and Biotechnology. 7, 381 (2019).
  5. Baillargeon, B., Rebelo, N., Fox, D. D., Taylor, R. L., Kuhl, E. The Living Heart Project: A robust and integrative simulator for human heart function. European Journal of Mechanics - A/Solids. 48, 38-47 (2014).
  6. Guccione, J. M., McCulloch, A. D. Mechanics of active contraction in cardiac muscle: Part I--Constitutive relations for fiber stress that describe deactivation. TheJournal of Biomechanical Engineering. 115, 72-81 (1993).
  7. Guccione, J. M., Waldman, L. K., McCulloch, A. D. Mechanics of active contraction in cardiac muscle: Part II--Cylindrical models of the systolic left ventricle. The Journal of Biomechanical Engineering. 115, 82-90 (1993).
  8. Kojic, M., et al. Smeared multiscale finite element model for electrophysiology and ionic transport in biological tissue. Computers in Biology and Medicine. 108, 288-304 (2019).
  9. Wang, Y., Rudy, Y. Application of the method of fundamental solutions to potential-based inverse electrocardiography. Annals of Biomedical Engineering. 34 (8), 1272-1288 (2006).
  10. Van Oosterom, A. The use of the spatial covariance in computing pericardial potentials. IEEE Transactions on Biomedical Engineering. 46 (7), 778-787 (1999).
  11. Van Oosterom, A. The spatial covariance used in computing the pericardial potential distribution. Computational Inverse Problems in Electrocardiography. , 1-50 (2001).
  12. Van Oosterom, A. Source models in inverse electrocardiography. International Journal of Bioelectromagnetism. 5, 211-214 (2003).
  13. Van Oosterom, A. The equivalent double layer: source models for repolarization. Comprehensive Electrocardiology. , 227-246 (2010).

转载和许可

请求许可使用此 JoVE 文章的文本或图形

请求许可

探索更多文章

183

This article has been published

Video Coming Soon

JoVE Logo

政策

使用条款

隐私

科研

教育

关于 JoVE

版权所属 © 2025 MyJoVE 公司版权所有,本公司不涉及任何医疗业务和医疗服务。