Este trabalho introduz dois modelos computacionais de insuficiência cardíaca com fração de ejeção preservada baseada em uma abordagem de parâmetros e análise de elementos finitos. Esses modelos são utilizados para avaliar as alterações na hemodinâmica do ventrículo esquerdo e vasculatura relacionada induzida pela sobrecarga de pressão e diminuição da conformidade ventricular.
Os esforços científicos no campo da modelagem computacional de doenças cardiovasculares têm se concentrado em grande parte na insuficiência cardíaca com redução da fração de ejeção (HFrEF), amplamente negligenciando a insuficiência cardíaca com a fração de ejeção preservada (HFpEF), que mais recentemente se tornou uma forma dominante de insuficiência cardíaca em todo o mundo. Motivados pela escassez de HFpEF em representações de silico, dois modelos computacionais distintos são apresentados neste artigo para simular a hemodinâmica do HFpEF resultante da sobrecarga de pressão ventricular esquerda. Primeiro, um modelo de parâmetros amontoados orientados a objeto foi desenvolvido usando um solucionador numérico. Este modelo é baseado em uma rede derivada de Windkessel em dimensões zero (0D), que depende das propriedades geométricas e mecânicas dos elementos constitutivos e oferece a vantagem de baixos custos computacionais. Em segundo lugar, foi utilizado um pacote de software de análise de elementos finitos (FEA) para a implementação de uma simulação multidimensional. O modelo FEA combina modelos multifísicos tridimensionais (3D) da resposta cardíaca eletromecânica, deformações estruturais e hemodinâmica à base de cavidade fluida e utiliza um modelo simplificado de parâmetros fixos para definir os perfis de troca de fluxo entre diferentes cavidades fluidas. Através de cada abordagem, tanto as alterações hemodinâmicas agudas e crônicas no ventrículo esquerdo quanto na vasculatura proximal resultantes da sobrecarga de pressão foram simuladas com sucesso. Especificamente, a sobrecarga de pressão foi modelada pela redução da área orifício da válvula aórtica, enquanto a remodelagem crônica foi simulada pela redução da conformidade da parede ventricular esquerda. Consistente com a literatura científica e clínica da HFpEF, os resultados de ambos os modelos mostram (i) uma elevação aguda do gradiente de pressão transaréltica entre o ventrículo esquerdo e a aorta e uma redução no volume do AVC e (ii) uma diminuição crônica no volume ventricular esquerdo terminante, indicativo de disfunção diastólica. Finalmente, o modelo FEA demonstra que o estresse no miocárdio HFpEF é notavelmente maior do que no tecido cardíaco saudável durante todo o ciclo cardíaco.
A insuficiência cardíaca é uma das principais causas de morte em todo o mundo, que ocorre quando o coração é incapaz de bombear ou preencher adequadamente para acompanhar as demandas metabólicas do corpo. A fração de ejeção, ou seja, a quantidade relativa de sangue armazenada no ventrículo esquerdo que é ejetada a cada contração é usada clinicamente para classificar a insuficiência cardíaca em (i) insuficiência cardíaca com fração de ejeção reduzida (HFrEF) e (ii) insuficiência cardíaca com fração de ejeção preservada (HFpEF), para frações de ejeção inferiores ou superiores a 45%,respectivamente 1,2,3. Os sintomas do HFpEF frequentemente se desenvolvem em resposta à sobrecarga de pressão ventricular esquerda, que pode ser causada por várias condições, incluindo estenose aórtica, hipertensão e obstrução do trato de saída ventricular esquerda3,4,5,6,7. A sobrecarga de pressão impulsiona uma cascata de aberrações moleculares e celulares, levando ao espessamento da parede ventricular esquerda (remodelagem concêntrica) e, finalmente, ao enrijecimento da parede ou perda de conformidade8,9,10. Essas alterações biomecânicas afetam profundamente a hemodinâmica cardiovascular, pois resultam em uma relação elevada de aumento do volume de pressão diastólica e na redução do volume diastólico final11.
A modelagem computacional do sistema cardiovascular avançou na compreensão das pressões e fluxos sanguíneos tanto na fisiologia quanto na doença e tem fomentado o desenvolvimento de estratégias diagnósticas e terapêuticas12. Em modelos silico são classificados em modelos de baixa ou alta dimensão, com o primeiro utilizando métodos analíticos para avaliar propriedades hemodinâmicas globais com baixa demanda computacional e este último fornecendo uma descrição multiescala e multifísica mais extensa da mecânica cardiovascular e hemodinâmica no domínio 2D ou 3D13. A representação de Windkessel com parâmetros é a mais comum entre as descrições de baixa dimensão. Com base na analogia do circuito elétrico (lei de Ohm), isso imita o comportamento hemodinâmico geral do sistema cardiovascular através de uma combinação de elementos resistivos, capacitivos e indutivos14. Um estudo recente deste grupo propôs um modelo alternativo windkessel no domínio hidráulico que permite a modelagem de mudanças na geometria e mecânica de grandes vasos-câmaras cardíacas e válvulas – de forma mais intuitiva do que os modelos analógicos elétricos tradicionais. Esta simulação é desenvolvida em um solucionador numérico orientado a objetos (ver a Tabela de Materiais) e pode capturar a hemodinâmica normal, efeitos fisiológicos do acoplamento cardiorrespiratório, fluxo sanguíneo respiratório na fisiologia uni-cardíaca e alterações hemodinâmicas devido à constrição aórtica. Esta descrição expande-se sobre as capacidades de modelos de parâmetros agrupados, oferecendo uma abordagem fisicamente intuitiva para modelar um espectro de condições patológicas, incluindo insuficiência cardíaca15.
Modelos de alta dimensão são baseados na FEA para calcular hemodinâmica esposteral e interações de estrutura fluida. Essas representações podem fornecer descrições detalhadas e precisas do campo de fluxo sanguíneo local; no entanto, devido à sua baixa eficiência computacional, eles não são adequados para estudos de toda a árvore cardiovascular16,17. Um pacote de software (ver a Tabela de Materiais) foi empregado como uma plataforma FEA anatomicamente precisa do coração humano adulto de 4 câmaras, que integra a resposta eletromecânica, deformações estruturais e hemodinâmica baseada em cavidade fluida. O modelo de coração humano adaptado também compreende um modelo simples de parâmetros amontoados que define a troca de fluxo entre as diferentes cavidades fluidas, bem como uma caracterização mecânica completa do tecido cardíaco18,19.
Vários modelos de metametro e FEA de insuficiência cardíaca foram formulados para capturar anormalidades hemodinâmicas e avaliar estratégias terapêuticas, particularmente no contexto de dispositivos de assistência circulatória mecânica para HFrEF20,21,22,23,24. Uma ampla gama de modelos de parâmetros 0D de várias complexidades capturou com sucesso a hemodinâmica do coração humano em condições fisiológicas e HFrEF através da otimização de dois ou três elementos sistemas analógicos elétricos Windkessel20,21,23,24. A maioria dessas representações são modelos uni ou biventriculares baseados na formulação de tempo variando-elastância para reproduzir a ação contratil do coração e usar uma relação de pressão-volume não linear para descrever o preenchimento ventricular25,26,27. Modelos abrangentes, que capturam a complexa rede cardiovascular e imitam a ação de bombeamento atrial e ventricular, têm sido usados como plataformas para testes de dispositivos. No entanto, embora exista um corpo significativo de literatura no campo da HFrEF, muito poucos em modelos silico de HFpEF foram propostos20,22,28,29,30,31.
Um modelo de baixa dimensão da hemodinâmica HFpEF, recentemente desenvolvido por Burkhoff et al.32 e Granegger et al.28,pode capturar os laços de volume de pressão (PV) do coração de 4 câmaras, recapitulando totalmente a hemodinâmica de vários fenótipos de HFpEF. Além disso, utilizam sua plataforma em silico para avaliar a viabilidade de um dispositivo circulatório mecânico para hfpef, pesquisa computacional pioneira do HFpEF para estudos de fisiologia, bem como desenvolvimento de dispositivos. No entanto, esses modelos permanecem incapazes de captar as mudanças dinâmicas nos fluxos sanguíneos e pressões observadas durante a progressão da doença. Um estudo recente de Kadry et al.30 captura os vários fenótipos de disfunção diastólica ajustando o relaxamento ativo do miocárdio e a rigidez passiva do ventrículo esquerdo em um modelo de baixa dimensão. Seu trabalho fornece uma análise hemodinâmica abrangente da disfunção diastólica com base nas propriedades ativas e passivas do miocárdio. Da mesma forma, a literatura de modelos de alta dimensão tem focado principalmente em HFrEF19,33,34,35,36,37. Bakir et al.33 propuseram um modelo FEA de fluido cardíaco totalmente acoplado para prever o perfil hemodinâmico HFrEF e a eficácia de um dispositivo de assistência à esquerda-ventricular (LVAD). Este modelo biventricular (ou duas câmaras) utilizou um circuito windkessel acoplado para simular a hemodinâmica do coração saudável, HFrEF e HFrEF com suporte LVAD33,37.
Da mesma forma, Sack et al.35 desenvolveram um modelo biventricular para investigar disfunção ventricular direita. Sua geometria biventricular foi obtida a partir dos dados de ressonância magnética (RM) do paciente, e a malha de elemento finito do modelo foi construída usando segmentação de imagem para analisar a hemodinâmica de um ventrículo direito com falha suportada pelo VAD35. Abordagens cardíacas FEA de quatro câmaras foram desenvolvidas para melhorar a precisão dos modelos do comportamento eletromecânico do coração19,34. Em contraste com as descrições biventriculares, modelos de quatro câmaras derivadas da ressonância magnética do coração humano fornecem uma melhor representação da anatomia cardiovascular18. O modelo cardíaco empregado neste trabalho é um exemplo estabelecido de um modelo FEA de quatro câmaras. Ao contrário dos modelos FEA de parâmetros e biventriculares, essa representação captura alterações hemodinâmicas à medida que ocorrem durante a progressão da doença34,37. Genet et al.34, por exemplo, utilizaram a mesma plataforma para implementar um modelo de crescimento numérico da remodelagem observada no HFrEF e HFpEF. No entanto, esses modelos avaliam os efeitos da hipertrofia cardíaca apenas na mecânica estrutural e não fornecem uma descrição abrangente da hemodinâmica associada.
Para suprida a falta de HFpEF em modelos de silico neste trabalho, o modelo de parâmetros agrupados anteriormente desenvolvido por este grupo15 e o modelo FEA foram readaptados para simular o perfil hemodinâmico da HFpEF. Para isso, a capacidade de cada modelo de simular hemodinâmica cardiovascular na linha de base será demonstrada pela primeira vez. Os efeitos da sobrecarga de pressão ventricular esquerda induzida por estenose e da diminuição da conformidade ventricular esquerda devido à remodelagem cardíaca – uma marca típica do HFpEF – serão então avaliados.
As plataformas de parâmetros e FEA propostas neste trabalho recapitularam a hemodinâmica cardiovascular em condições fisiológicas, tanto na fase aguda da sobrecarga de pressão induzida por estenose quanto na HFpEF crônica. Ao captar o papel que a sobrecarga de pressão desempenha nas fases agudas e crônicas do desenvolvimento da HFpEF, os resultados desses modelos estão de acordo com a literatura clínica da HFpEF, incluindo o surgimento de um gradiente de pressão transatróptica devido à estenose aórtica, au…
The authors have nothing to disclose.
Reconhecemos o financiamento do programa de Ciências e Tecnologia do Instituto de Tecnologia de Harvard-Massachusetts e do Prêmio Fundação SITA do Instituto de Engenharia Médica e Ciência.
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 |