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

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.