Parámetros Agrupados Y Modelado De Elementos Finitos De La Insuficiencia Cardíaca Con Fracción De Eyección Preservada

Summary

Este trabajo introduce dos modelos computacionales de insuficiencia cardíaca con fracción de eyección preservada basados en un enfoque de parámetros agrupados y análisis de elementos finitos. Estos modelos se utilizan para evaluar los cambios en la hemodinámica del ventrículo izquierdo y la vasculatura relacionada inducida por sobrecarga de presión y disminución del cumplimiento ventricular.

Abstract

Los esfuerzos científicos en el campo del modelado computacional de enfermedades cardiovasculares se han centrado en gran medida en la insuficiencia cardíaca con fracción de eyección reducida (HFrEF), pasando por alto ampliamente la insuficiencia cardíaca con fracción de eyección preservada (HFpEF), que más recientemente se ha convertido en una forma dominante de insuficiencia cardíaca en todo el mundo. Motivados por la escasez de representaciones in silico de HFpEF, dos modelos de cómputo distintos se presentan en este papel para simular la hemodinámica de HFpEF resultando de sobrecarga ventricular izquierda de la presión. En primer lugar, se desarrolló un modelo de parámetros agrupados orientado a objetos utilizando un solucionador numérico. Este modelo se basa en una red de dimensión cero (0D) similar a Windkessel, que depende de las propiedades geométricas y mecánicas de los elementos constitutivos y ofrece la ventaja de bajos costos computacionales. En segundo lugar, se utilizó un paquete de software de análisis de elementos finitos (FEA) para la implementación de una simulación multidimensional. El modelo FEA combina modelos multifísicos tridimensionales (3D) de la respuesta cardíaca electromecánica, deformaciones estructurales y hemodinámica basada en cavidades de fluidos y utiliza un modelo simplificado de parámetros agrupados para definir los perfiles de intercambio de flujo entre diferentes cavidades de fluidos. Con cada acercamiento, los cambios hemodinámicos agudos y crónicos en el ventrículo izquierdo y la vasculatura próxima resultando de sobrecarga de la presión fueron simulados con éxito. Específicamente, la sobrecarga de presión se modeló mediante la reducción del área del orificio de la válvula aórtica, mientras que la remodelación crónica se simuló mediante la reducción del cumplimiento de la pared ventricular izquierda. De acuerdo con la literatura científica y clínica de la FEPH, los resultados de ambos modelos muestran (i) una elevación aguda del gradiente de presión transaórtico entre el ventrículo izquierdo y la aorta y una reducción en el volumen del accidente cerebrovascular y (ii) una disminución crónica del volumen ventricular izquierdo diastólico final, indicativo de disfunción diastólica. Finalmente, el modelo FEA demuestra que el estrés en el miocardio HFpEF es notablemente mayor que en el tejido cardíaco sano a lo largo del ciclo cardíaco.

Introduction

La insuficiencia cardíaca es una de las principales causas de muerte en todo el mundo, que ocurre cuando el corazón es incapaz de bombear o llenar adecuadamente para mantenerse al día con las demandas metabólicas del cuerpo. La fracción de eyección, es decir, la cantidad relativa de sangre almacenada en el ventrículo izquierdo que se expulsa con cada contracción se utiliza clínicamente para clasificar la insuficiencia cardíaca en (i) insuficiencia cardíaca con fracción de eyección reducida (HFrEF) y (ii) insuficiencia cardíaca con fracción de eyección preservada (HFpEF), para fracciones de eyección menores o superiores al 45%, respectivamente^1,2,3. Los síntomas de HFpEF se convierten a menudo en respuesta a la sobrecarga ventricular izquierda de la presión, que se puede causar por varias condiciones incluyendo estenosis aórtica, hipertensión, y obstrucción ventricular izquierda de la zona de salida^3,4,5,6,7. La sobrecarga de presión conduce a una cascada de aberraciones moleculares y celulares, lo que conduce al engrosamiento de la pared ventricular izquierda (remodelación concéntrica) y, en última instancia, al endurecimiento de la pared o a la pérdida de cumplimiento^8,9,10. Estos cambios biomecánicos afectan profundamente a la hemodinámica cardiovascular, ya que dan lugar a una elevada relación presión-volumen diastólica final y a una reducción del volumen diastólico final^11.

La modelización computacional del sistema cardiovascular ha avanzado en la comprensión de las presiones y flujos sanguíneos tanto en la fisiología como en la enfermedad y ha fomentado el desarrollo de estrategias diagnósticas y terapéuticas^12. Los modelos in silico se clasifican en modelos de baja o alta dimensión, utilizándose los primeros métodos analíticos para evaluar las propiedades hemodinámicas globales con baja demanda computacional y los segundos proporcionando una descripción multiescala y multifísica más extensa de la mecánica cardiovascular y la hemodinámica en el dominio 2D o 3D^13. La representación de Windkessel de parámetro agrupado es la más común entre las descripciones de baja dimensión. Basado en la analogía del circuito eléctrico (ley de Ohm), esto imita el comportamiento hemodinámico general del sistema cardiovascular a través de una combinación de elementos resistivos, capacitivos e inductivos^14. Un estudio reciente de este grupo ha propuesto un modelo de Windkessel alternativo en el dominio hidráulico que permite el modelado de cambios en la geometría y mecánica de grandes recipientes -cámaras de corazón y válvulas- de una manera más intuitiva que los modelos analógicos eléctricos tradicionales. Esta simulación se desarrolla en un solucionador numérico orientado a objetos (ver la Tabla de Materiales)y puede capturar la hemodinámica normal, los efectos fisiológicos del acoplamiento cardiorrespiratorio, el flujo sanguíneo impulsado por las vías respiratorias en la fisiología de un solo corazón y los cambios hemodinámicos debidos a la constricción aórtica. Esta descripción amplía las capacidades de los modelos de parámetros agrupados al ofrecer un enfoque físicamente intuitivo para modelar un espectro de afecciones patológicas, incluida la insuficiencia cardíaca^15.

Los modelos de alta dimensión se basan en FEA para calcular la hemodinámica espaciotemporal y las interacciones fluido-estructura. Estas representaciones pueden proporcionar descripciones detalladas y precisas del campo de flujo sanguíneo local; sin embargo, debido a su baja eficiencia computacional, no son adecuados para estudios de todo el árbol cardiovascular^16,17. Se empleó un paquete de software (ver la Tabla de Materiales)como una plataforma FEA anatómicamente precisa del corazón humano adulto de 4 cámaras, que integra la respuesta electromecánica, las deformaciones estructurales y la hemodinámica basada en cavidades fluidas. El modelo de corazón humano adaptado también comprende un modelo simple de parámetros agrupados que define el intercambio de flujo entre las diferentes cavidades de fluidos, así como una caracterización mecánica completa del tejido cardíaco^18,19.

Se han formulado varios modelos de insuficiencia cardíaca de parámetros agrupados y FEA para capturar anomalías hemodinámicas y evaluar estrategias terapéuticas, particularmente en el contexto de los dispositivos mecánicos de asistencia circulatoria para HFrEF^{20,21,22,23,24.} Por lo tanto, una amplia gama de modelos de parámetros agrupados en 0D de diversas complejidades ha capturado con éxito la hemodinámica del corazón humano en condiciones fisiológicas y HFrEF a través de la optimización de sistemas windkessel analógicos eléctricos de dos o tres^{elementos 20,21,23,24.} La mayoría de estas representaciones son modelos uni o biventriculares basados en la formulación de la elastance variable en el tiempo para reproducir la acción contráctil del corazón y utilizar una relación no lineal de presión-volumen diastólica final para describir el llenado ventricular^25,26,27. Los modelos comprensivos, que capturan la red cardiovascular compleja e imitan la acción de bombeo atrial y ventricular, se han utilizado como plataformas para la prueba del dispositivo. Sin embargo, aunque existe un cuerpo significativo de literatura en torno al campo de la HFrEF, muy pocos modelos in silico de HFpEF se han propuesto^{20,22,28,29,30,31.}

Un modelo de baja dimensión de la hemodinámica de la HFpEF, desarrollado recientemente por Burkhoff et al.³² y Granegger et al.^28,puede capturar los bucles presión-volumen (PV) del corazón de 4 cámaras, recapitulando completamente la hemodinámica de varios fenotipos de HFpEF. Además, utilizan su plataforma in silico para evaluar la viabilidad de un dispositivo circulatorio mecánico para HFpEF, siendo pioneros en la investigación computacional de HFpEF para estudios de fisiología, así como el desarrollo de dispositivos. Sin embargo, estos modelos siguen siendo incapaces de capturar los cambios dinámicos en los flujos sanguíneos y las presiones observadas durante la progresión de la enfermedad. Un estudio reciente de Kadry et al.³⁰ captura los diversos fenotipos de disfunción diastólica ajustando la relajación activa del miocardio y la rigidez pasiva del ventrículo izquierdo en un modelo de baja dimensión. Su trabajo proporciona un análisis hemodinámico completo de la disfunción diastólica basado en las propiedades activas y pasivas del miocardio. Del mismo modo, la literatura de modelos de alta dimensión se ha centrado principalmente en HFrEF^{19,33,34,35,36,37.} Bakir et al.³³ propusieron un modelo FEA de fluido cardíaco-electromecánica totalmente acoplado para predecir el perfil hemodinámico de la HFrEF y la eficacia de un dispositivo de asistencia ventricular izquierda (DAVI). Este modelo biventricular (o de dos cámaras) utilizó un circuito de Windkessel acoplado para simular la hemodinámica del corazón sano, HFrEF y HFrEF con soporte de LVAD^33,37.

Del mismo modo, Sack et al.³⁵ desarrollaron un modelo biventricular para investigar la disfunción ventricular derecha. Su geometría biventricular se obtuvo a partir de los datos de resonancia magnética (RM) de un paciente, y la malla de elementos finitos del modelo se construyó utilizando la segmentación de imágenes para analizar la hemodinámica de un ventrículo derecho defectuoso apoyado por^{DAV 35.} Se han desarrollado enfoques cardíacos fea de cuatro cámaras para mejorar la precisión de los modelos del comportamiento electromecánico del^{corazón 19,34}. En contraste con las descripciones biventriculares, los modelos de cuatro cámaras del corazón humano derivados de la RMN proporcionan una mejor representación de la anatomía cardiovascular¹⁸. El modelo de corazón empleado en este trabajo es un ejemplo establecido de un modelo FEA de cuatro cámaras. A diferencia de los modelos de parámetros agrupados y fea biventricular, esta representación captura los cambios hemodinámicos a medida que ocurren durante la progresión de la enfermedad^34,37. Genet et al.^34,por ejemplo, utilizaron la misma plataforma para implementar un modelo numérico de crecimiento de la remodelación observada en HFrEF y HFpEF. Sin embargo, estos modelos evalúan los efectos de la hipertrofia cardiaca sobre mecánico estructural solamente y no proporcionan una descripción comprensiva de la hemodinámica asociada.

Para abordar la falta de modelos in silico de HFpEF en este trabajo, el modelo de parámetros agrupados previamente desarrollado por este grupo¹⁵ y el modelo FEA fueron readaptados para simular el perfil hemodinámico de HFpEF. Con este fin, primero se demostrará la capacidad de cada modelo para simular la hemodinámica cardiovascular al inicio del estudio. Los efectos de la sobrecarga ventricular izquierda estenosis-inducida de la presión y de la conformidad ventricular izquierda disminuida debido al remodelado cardiaco-un sello típico de HFpEF-entonces serán evaluados.

Protocol

1. Modelo de parámetros agrupados 0D Configuración de simulaciónNOTA: En el entorno del solucionador numérico (consulte la Tabla de materiales),construya el dominio como se muestra en la Figura 1. Esto se compone del corazón de 4 cámaras, la parte superior del cuerpo, abdominal, parte inferior del cuerpo y compartimentos torácicos, así como la vasculatura proximal, incluyendo la aorta, la arteria pulmonar y las venas cavas superior e inferior. Los eleme…

Representative Results

Los resultados de las simulaciones de línea base se ilustran en la Figura 3. Esto representa las formas de onda de presión y volumen del ventrículo izquierdo y la aorta (Figura 3A) así como el bucle fotovoltaico ventricular izquierdo (Figura 3B). Los dos modelos in silico muestran hemodinámica aórtica y ventricular izquierda similares, que se encuentran dentro del rango fisiológico. Se pueden notar diferencias menores en la r…

Discussion

Las plataformas lumped-parameter y FEA propuestas en este trabajo recapitularon la hemodinámica cardiovascular en condiciones fisiológicas, tanto en la fase aguda de sobrecarga de presión inducida por estenosis como en la FEPH crónica. Al capturar el papel que juega la sobrecarga de presión en las fases aguda y crónica del desarrollo de la FEPH, los resultados de estos modelos están de acuerdo con la literatura clínica de la FEPH, incluyendo la aparición de un gradiente de presión transaórtico debido a la este…

Materials

Abaqus Software	Dassault Systèmes Simulia Corp.		Version used: 2018; FEA simulation software
HETVAL	Dassault Systèmes Simulia Corp.		Version used: 2018
Hydraulic (Isothermal) library	MathWorks		Version used: 2020a
Living Heart Human Model	Dassault Systèmes Simulia Corp.		Version used: V2_1, anatomically accurate FEA platform of 4-chamber adult human heart
MATLAB	MathWorks		Version used: 2020a, object-oriented numerical solver
SIMSCAPE FLUIDS	MathWorks
UAMP	Dassault Systèmes Simulia Corp.		Version used: 2018
VUANISOHYPER	Dassault Systèmes Simulia Corp.		Version used: 2018