ECMSim: A high-performance interactive web application for real-time spatiotemporal simulation of cardiac ECM signaling and diffusion

Resumo Técnico do ECMSim: Uma Aplicação Web Interativa de Alto Desempenho para Simulação Espaciotemporal em Tempo Real de Sinalização e Difusão da Matriz Extracelular Cardíaca

Declaração do Problema
O remodelamento da matriz extracelular (MEC) é um processo computacionalmente intensivo e mecanisticamente complexo, central à fisiopatologia cardiovascular, particularmente na fibrose cardíaca e na insuficiência cardíaca. A previsão precisa da dinâmica da MEC requer a integração simultânea de dinâmicas de reação-difusão multiescala, retroalimentação mecanoquímica e respostas celulares espacialmente heterogêneas. Os modelos computacionais existentes enfrentam limitações significativas: muitos carecem da resolução espacial necessária para modelar interações célula-célula e organização tecidual; frequentemente exigem instalação de software especializado e expertise, limitando a acessibilidade; e geralmente carecem de interatividade em tempo real para exploração de parâmetros. Além disso, as abordagens experimentais tradicionais lutam para capturar a complexidade total das dinâmicas de sinalização de múltiplas vias combinadas com difusão molecular espacial e regulação por retroalimentação.

Metodologia
Os autores apresentam o ECMSim, uma plataforma baseada na web projetada para simular o remodelamento heterogêneo da matriz em tecido cicatricial cardíaco. O sistema é construído sobre um modelo abrangente, de nível sistêmico, da rede de sinalização de fibroblastos cardíacos, estendendo um modelo ODE baseado em lógica de célula única previamente validado (Zeigler et al.) para um framework espacialmente resolvido.

• Framework Matemático: O núcleo do ECMSim consiste em um sistema de 1,37 milhão de equações diferenciais ordinárias (EDOs) acopladas. Este sistema modela uma grade espacial de 100 × 100 (10.000 células), onde cada célula contém 137 espécies moleculares (115 espécies de sinalização intracelular e 17 espécies de MEC extracelular). A rede inclui 259 interações regulatórias (arestas) cobrindo dez moléculas de entrada (por exemplo, TGF-β, AngII, IL6, tensão), receptores de membrana, cascatas intracelulares (MAPK, PI3K-Akt-mTOR, GTPase Rho), fatores de transcrição e saídas de MEC (colágenos, MMPs, TIMPs).
• Dinâmica Espacial: A comunicação intercelular é mediada pela difusão de cinco moléculas de sinalização de retroalimentação chave (TGF-βfb, AngIIfb, IL6fb, ET1fb e tensionfb) e difusão limitada de moléculas de MEC. A evolução espacial é modelada usando equações de difusão discretas com aproximações de diferenças finitas em uma grade celular. Condições de contorno periódicas são implementadas para minimizar efeitos de tamanho finito. As moléculas de MEC difundem a 20% da taxa das moléculas de retroalimentação para refletir seu maior tamanho molecular e tendência para deposição local.
• Integração Numérica: O sistema é resolvido usando o método de Euler forward com passo de tempo adaptativo. Para garantir estabilidade numérica e relevância fisiológica, as concentrações moleculares são limitadas dentro de uma faixa normalizada [0, 1].
• Implementação Computacional: O motor de simulação é escrito em C++ otimizado e compilado para WebAssembly usando a ferramenta Emscripten. Isso permite que o aplicativo seja executado diretamente em navegadores web padrão com desempenho aproximando-se do código C++ nativo (dentro de ~20% da velocidade nativa). A arquitetura utiliza design orientado a objetos, estratégias de memória pré-alocada e otimizações de localidade de cache para lidar com a carga computacional massiva.
• Interface do Usuário: A plataforma apresenta uma interface interativa com uma ferramenta no estilo "pincel" para selecionar subconjuntos de células para aplicar concentrações ajustáveis de moléculas de entrada. Inclui controles deslizantes de parâmetro para ajuste em tempo real de nove constantes de taxa primárias (por exemplo, processamento de entrada, sensibilidade de retroalimentação, degradação, taxas de difusão) e múltiplas visualizações acopladas (mapas de calor, gráficos de linha, campos vetoriais e superfícies 3D).

Contribuições Principais

1. Modelagem Espacial Escalável: O ECMSim expande a modelagem de fibroblastos cardíacos de representações de célula única para uma grade espacialmente resolvida de 10.000 células com acoplamento de reação-difusão, permitindo o estudo da organização em nível tecidual e padrões espaciais emergentes.
2. Interatividade em Tempo Real: Ao aproveitar o WebAssembly, a plataforma alcança desempenho em tempo real para simulações envolvendo 1,37 milhão de variáveis, permitindo que os usuários explorem interativamente espaços de parâmetros e observem respostas dinâmicas sem a latência típica de clusters de computação de alto desempenho.
3. Acessibilidade e Democratização: A implantação baseada na web elimina a necessidade de instalação de software especializado, tornando a modelagem sofisticada de biologia de sistemas acessível a uma gama mais ampla de pesquisadores e educadores em diferentes plataformas.
4. Integração Abrangente de Redes: O modelo integra as principais vias de sinalização implicadas na fibrose cardíaca (TGF-β/Smad, MAPK, PI3K-Akt-mTOR, mecanotransdução) em um framework unificado que considera loops de retroalimentação autócrina e parácrina.

Resultados

• Benchmarks de Desempenho: Em um MacBook Pro (Apple M1 Max), o ECMSim executa aproximadamente 4,84 milhões de operações por passo de tempo, incluindo avaliações de taxa de EDO, atualizações de integração de Euler e computações de vizinhança de difusão. O sistema mantém interatividade em tempo real com uma pegada de memória total de aproximadamente 523 MB, com mais de 99% do tempo de CPU dedicado ao solucionador de EDO.
• Dinâmicas Espaciotemporais: As simulações demonstram a capacidade da plataforma de modelar a propagação de respostas fibróticas. Estímulo localizado (por exemplo, seleção em forma de X com TGF-β) desencadeia cascatas intracelulares que levam à produção de MEC, que então difunde para células vizinhas, criando gradientes de concentração espacialmente heterogêneos. O sistema captura o atraso entre a ativação do sinal e o acúmulo de MEC, bem como os efeitos sinérgicos da estimulação combinatória de citocinas (por exemplo, TGF-β + AngII + IL6).
• Sensibilidade de Parâmetros: A análise de sensibilidade revela que a taxa de produção de MEC ($k_{production}$) é o regulador dominante dos resultados do remodelamento da MEC, seguida pelas taxas de ativação e degradação. O modelo é robusto a variações moderadas em outros parâmetros, consistente com o entendimento biológico de que o acúmulo de matriz é impulsionado principalmente pelas taxas de produção em fibroblastos ativados.

Significado e Alegações
Os autores afirmam que o ECMSim representa um avanço significativo na modelagem computacional acessível para sistemas biológicos complexos. Ao combinar rigor matemático com tecnologias web amigáveis ao usuário, a plataforma aborda lacunas críticas na pesquisa atual de fibrose cardíaca, especificamente a falta de resolução espacial e interatividade em tempo real.

O artigo postula que o ECMSim facilita a investigação de condições patológicas, cenários hipotéticos e os efeitos de drogas experimentais em receptores-alvo sem as restrições temporais e de recursos de experimentos tradicionais. Serve como um "laboratório digital" para testar hipóteses e validar modelos. Embora focado principalmente na fibrose cardíaca, os autores observam que o framework subjacente de reação-difusão fornece uma base para modelar outras aplicações, incluindo microambientes de câncer, cicatrização de feridas e biologia regenerativa. A plataforma visa preencher a lacuna entre observações experimentais em escala molecular e fisiopatologia em nível tecidual, fomentando colaboração interdisciplinar e acelerando a exploração de novas estratégias terapêuticas.