The functional impact of myofiber macroscopic organization and disarray in computational models of the murine heart

Título: El impacto funcional de la organización macroscópica y el desorden de las fibras miocárdicas en modelos computacionales del corazón murino

1. Problema

La reconstrucción precisa de la arquitectura de las fibras miocárdicas es un desafío crítico en los modelos computacionales de electromecánica cardíaca.

• Limitación de medición: Las mediciones in vivo de la orientación de las fibras son inviables con la tecnología actual.
• Uso de sustitutos: Debido a esto, se utilizan métodos basados en reglas (Rule-Based Methods, RBM), como los métodos Laplace-Dirichlet (LDRBM), que generan campos de fibras ideales basados en la geometría del órgano.
• El problema del desorden (Disarray): Estos métodos suelen ignorar el "desorden de las fibras", que son las desviaciones locales de la orientación media. La literatura carece de validaciones exhaustivas que comparen modelos con RBM frente a datos experimentales de alta resolución, y no está claro si el desorden microscópico afecta significativamente la función cardíaca o si las estrategias actuales para modelarlo (como la activación cruzada de fibras) son adecuadas.

2. Metodología

Los autores desarrollaron un modelo electromecánico biventricular de alta fidelidad de un corazón de ratón y aplicaron una estrategia innovadora para separar la arquitectura macroscópica del desorden microscópico.

• Datos de Referencia (Ground Truth): Se utilizó un campo de fibras de alta resolución (96 µm) obtenido mediante imágenes ópticas mesoscópicas (mesoSPIM) en un corazón de ratón ex vivo. Este campo sirve como referencia experimental.
• Desacoplamiento Frecuencial: Se introdujo una estrategia de filtrado espacial (filtro de Helmholtz) para descomponer el campo de fibras medido en dos componentes:
  1. Arquitectura macroscópica: La tendencia de baja frecuencia (organización helicoidal general).
  2. Desorden microscópico: La señal residual de alta frecuencia (desviaciones locales).
  • Se varió el radio de regularización ($\ell$) para suavizar progresivamente el campo, aislando el efecto del desorden.
• Modelado Comparativo: Se compararon tres escenarios:
  1. Campo experimental: Con desorden natural.
  2. Campos regularizados: Con desorden eliminado pero manteniendo la arquitectura media.
  3. LDRBM: Fibras generadas por métodos basados en reglas (sin desorden y con organización helicoidal idealizada).
• Estrategias de Modelado del Desorden: Se evaluaron dos enfoques comunes para simular el desorden en modelos con fibras ideales:
  1. Reducción de la rigidez efectiva de los puentes cruzados ($a_{XB}$) en el modelo de tensión activa.
  2. Introducción de componentes de tensión activa en direcciones ortogonales a la fibra (direcciones de lámina y normal).
• Simulación: Se resolvió un modelo de electromecánica acoplado con circulación sanguínea (0D) utilizando el framework lifex. Se analizaron la propagación eléctrica, la mecánica pasiva y la contracción activa.

3. Contribuciones Clave

• Validación de Alta Resolución: Es uno de los primeros estudios que integra datos de fibras a resolución celular (96 µm) en un modelo electromecánico completo para validar métodos de sustitución (RBM).
• Metodología de Desacoplamiento: Propone un método matemático riguroso para separar cuantitativamente el impacto de la organización global de las fibras frente al desorden local, algo que no se había hecho antes de esta manera sistemática.
• Evaluación Crítica de Estrategias: Proporciona una evaluación directa de la eficacia de las estrategias actuales para modelar el desorden (reducción de contractilidad vs. tensión ortogonal) comparándolas con un "ground truth" experimental.

4. Resultados

• Electrofisiología: La propagación del impulso eléctrico es poco sensible al desorden de las fibras. Aunque el desorden hace la propagación ligeramente más isotrópica, los tiempos de activación globales y los patrones de activación se mantienen similares entre el campo experimental, el regularizado y el LDRBM.
• Mecánica Pasiva: La respuesta pasiva (inflado del ventrículo) es débilmente afectada por el desorden. La arquitectura macroscópica influye en la rigidez, pero la eliminación del desorden microscópico solo produce cambios mínimos en el volumen (aprox. 0.4%).
• Mecánica Activa (Contracción): Aquí se observó la mayor sensibilidad.
  • Arquitectura Macroscópica: Es el factor dominante. La organización helicoidal es crucial para una contracción eficiente.
  • Efecto del Desorden: Una regularización moderada (eliminación de desorden, $\ell \approx 0.25$ mm) mejora la eficiencia de bombeo (fracción de eyección) al permitir que las fibras actúen de manera más coordinada. Sin embargo, una regularización excesiva (que destruye la organización helicoidal) o el uso de LDRBM (demasiado idealizado) conduce a contracciones poco fisiológicas (demasiado fuertes o ineficientes) y fallos numéricos.
• Evaluación de Estrategias de Sustitución:
  • Reducción de Contractilidad: Ajustar la rigidez de los puentes cruzados ($a_{XB}$) para reducir la contractilidad global reprodujo con mayor precisión los patrones de deformación (strain) locales y los bucles de presión-volumen (PV) del caso experimental.
  • Tensión Ortogonal: Introducir componentes de tensión activa en direcciones cruzadas (lámina/normal) logró reproducir los indicadores hemodinámicos globales (bucles PV), pero falló en capturar los patrones locales de deformación, generando distorsiones espurias y no reproduciendo la heterogeneidad observada experimentalmente.

5. Significado e Implicaciones

• Jerarquía de Importancia: La arquitectura macroscópica de las fibras es el determinante principal de la mecánica activa, mientras que el desorden microscópico tiene un impacto secundario pero significativo en la eficiencia de la contracción y la distribución local de la deformación.
• Limitaciones de los Modelos Actuales: Los métodos basados en reglas (RBM) son suficientes para predecir la propagación eléctrica y la mecánica pasiva global, pero son insuficientes para modelar con precisión la mecánica activa y la deformación tisular local si no se incorporan estrategias adecuadas para el desorden.
• Mejora de Gemelos Digitales: Para crear "gemelos digitales" cardíacos precisos, es prioritario capturar la organización macroscópica de las fibras. Si se utilizan fibras ideales (RBM), la estrategia de reducción de la contractilidad efectiva es superior a la adición de tensiones ortogonales para simular el desorden, ya que ofrece una mejor concordancia con la deformación tisular real.
• Futuro: El estudio sugiere que se necesitan procedimientos de escalado rigurosos que vinculen la distribución microscópica de fibras con formulaciones macroscópicas efectivas de tensión activa, en lugar de depender de ajustes fenomenológicos simples.

En resumen, el trabajo demuestra que, aunque el desorden de las fibras no altera drásticamente la señal eléctrica, es crucial para la eficiencia mecánica y la deformación local, y que las estrategias actuales de modelado deben priorizar la fidelidad de la arquitectura macroscópica y el uso de reducciones de contractilidad para compensar la falta de datos de desorden real.