🧬 biology

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

The paper introduces ECMSim, a high-performance, interactive web application that utilizes WebAssembly and compiled C++ to solve 1.37 million coupled ordinary differential equations in real time, enabling the spatiotemporal simulation of cardiac extracellular matrix remodeling and fibrosis across a 10,000-cell spatial array.

Original authors: Hasi Hays, William J. Richardson

Published 2026-05-26

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Hasi Hays, William J. Richardson

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ⚕️ This is an AI-generated explanation of a preprint that has not been peer-reviewed. It is not medical advice. Do not make health decisions based on this content. Read full disclaimer

Imagine your heart is a bustling city. When the city gets damaged (like after a heart attack), it tries to repair itself by building a "scar." This scar is made of a tough, sticky material called the Extracellular Matrix (ECM), which is like the city's concrete and steel.

Sometimes, the city builds too much concrete. This is called cardiac fibrosis, and it makes the heart stiff and unable to pump blood properly, leading to heart failure.

The problem is that figuring out exactly how and why this happens is incredibly hard. It involves millions of tiny chemical messages passing between cells, changing over time, and spreading across space. Trying to simulate this on a computer usually requires super-computers and complex software that only experts can use.

Enter ECMSim.

Think of ECMSim as a high-tech, interactive video game that runs right in your web browser. It simulates the heart's repair city in real-time, allowing anyone to see how the "concrete" (ECM) builds up and spreads.

Here is how it works, broken down into simple concepts:

1. The Digital City Grid

Imagine a giant checkerboard with 10,000 squares (100 by 100). Each square represents a single heart cell (a fibroblast).

Inside every single square, there is a tiny, complex factory running 137 different chemical processes at once.
The whole system solves 1.37 million math equations simultaneously, every single second, to figure out what each cell is doing.

2. The "Magic Paintbrush"

This is the coolest part. Instead of typing in code, you use a digital paintbrush.

You can "paint" over a specific group of cells on the checkerboard.
When you paint, you can tell those cells to receive a specific chemical signal (like a stress signal or a growth factor).
You can watch in real-time as the painted cells react, send out signals to their neighbors, and start building more "concrete" (ECM).

3. The Ripple Effect (Diffusion)

In the real world, chemicals don't stay in one place; they drift to neighbors.

In ECMSim, when a painted cell produces a signal, it diffuses (spreads out) to the surrounding cells, just like a drop of ink spreading in water.
This allows the simulation to show how a small, localized injury can cause a wave of scarring to spread across the whole tissue.
The app even shows you this spreading as colorful heat maps, 3D hills, and flowing arrows so you can "see" the invisible chemistry.

4. The Engine Under the Hood

Usually, running a simulation this big would require a massive supercomputer. But the creators used a special technology called WebAssembly.

Think of WebAssembly as a "universal translator" that lets a web browser run heavy-duty math code almost as fast as a native computer program.
This means you don't need to install anything. You just open the website, and it crunches 4.84 million calculations every second right on your laptop or tablet.

5. What Can You Do With It?

The paper describes ECMSim as a digital laboratory.

Test Hypotheses: You can ask, "What happens if I add a lot of this specific chemical to just the top-left corner?" and watch the result instantly.
Mix and Match: You can combine different signals (like stress + inflammation) to see how they work together to cause scarring.
Tweak the Rules: You can slide a bar to change how fast the cells react or how fast the chemicals break down, seeing how the "city" changes its behavior immediately.

Why This Matters (According to the Paper)

The paper claims that ECMSim fills a huge gap. Previous tools were either too simple (ignoring how cells talk to neighbors) or too hard to use (requiring special software). ECMSim combines high-speed math, spatial realism (seeing where things happen), and easy-to-use web tools into one package.

It allows researchers to explore "what-if" scenarios about heart scarring and drug effects without needing a wet lab or a supercomputer, making complex biology accessible to anyone with a web browser.

Technical Summary of ECMSim: A High-Performance Interactive Web Application for Real-Time Spatiotemporal Simulation of Cardiac ECM Signaling and Diffusion

Problem Statement
Extracellular matrix (ECM) remodeling is a computationally intensive and mechanistically complex process central to cardiovascular pathophysiology, particularly in cardiac fibrosis and heart failure. Accurately predicting ECM dynamics requires the simultaneous integration of multi-scale reaction-diffusion dynamics, mechanochemical feedback, and spatially heterogeneous cellular responses. Existing computational models face significant limitations: many lack the spatial resolution necessary to model cell-cell interactions and tissue organization; they often require specialized software installation and expertise, limiting accessibility; and they generally lack real-time interactivity for parameter exploration. Furthermore, traditional experimental approaches struggle to capture the full complexity of multi-pathway signaling dynamics combined with spatial molecular diffusion and feedback regulation.

Methodology
The authors introduce ECMSim, a web-based platform designed to simulate heterogeneous matrix remodeling in cardiac scar tissue. The system is built upon a comprehensive, system-level model of the cardiac fibroblast signaling network, extending a previously validated single-cell logic-based ODE model (Zeigler et al.) into a spatially resolved framework.

Mathematical Framework: The core of ECMSim consists of a system of 1.37 million coupled ordinary differential equations (ODEs). This system models a 100 × 100 spatial grid (10,000 cells), where each cell contains 137 molecular species (115 intracellular signaling species and 17 extracellular ECM species). The network includes 259 regulatory interactions (edges) covering ten input molecules (e.g., TGF-β, AngII, IL6, tension), membrane receptors, intracellular cascades (MAPK, PI3K-Akt-mTOR, Rho GTPase), transcription factors, and ECM outputs (collagens, MMPs, TIMPs).
Spatial Dynamics: Intercellular communication is mediated through the diffusion of five key feedback signaling molecules (TGF-βfb, AngIIfb, IL6fb, ET1fb, and tensionfb) and limited diffusion of ECM molecules. The spatial evolution is modeled using discrete diffusion equations with finite difference approximations on a cellular grid. Periodic boundary conditions are implemented to minimize finite-size effects. ECM molecules diffuse at 20% the rate of feedback molecules to reflect their larger molecular size and tendency for local deposition.
Numerical Integration: The system is solved using the forward Euler method with adaptive time-stepping. To ensure numerical stability and physiological relevance, molecular concentrations are clamped within a normalized range [0, 1].
Computational Implementation: The simulation engine is written in optimized C++ and compiled to WebAssembly using the Emscripten toolchain. This allows the application to run directly in standard web browsers with performance approaching native C++ code (within ~20% of native speed). The architecture utilizes object-oriented design, pre-allocated memory strategies, and cache-locality optimizations to handle the massive computational load.
User Interface: The platform features an interactive interface with a "brush-style" tool for selecting subsets of cells to apply adjustable input molecule concentrations. It includes parameter sliders for real-time adjustment of nine primary rate constants (e.g., input processing, feedback sensitivity, degradation, diffusion rates) and multiple coupled visualizations (heatmaps, line plots, vector fields, and 3D surfaces).

Key Contributions

Scalable Spatial Modeling: ECMSim expands cardiac fibroblast modeling from single-cell representations to a spatially resolved 10,000-cell grid with reaction-diffusion coupling, enabling the study of tissue-level organization and emergent spatial patterns.
Real-Time Interactivity: By leveraging WebAssembly, the platform achieves real-time performance for simulations involving 1.37 million variables, allowing users to interactively explore parameter spaces and observe dynamic responses without the latency typical of high-performance computing clusters.
Accessibility and Democratization: The web-based deployment eliminates the need for specialized software installation, making sophisticated systems biology modeling accessible to a broader range of researchers and educators across different platforms.
Comprehensive Network Integration: The model integrates major signaling pathways implicated in cardiac fibrosis (TGF-β/Smad, MAPK, PI3K-Akt-mTOR, mechanotransduction) into a unified framework that accounts for autocrine and paracrine feedback loops.

Results

Performance Benchmarks: On a MacBook Pro (Apple M1 Max), ECMSim executes approximately 4.84 million operations per time step, including ODE rate evaluations, Euler integration updates, and diffusion neighbor computations. The system maintains real-time interactivity with a total memory footprint of approximately 523 MB, with over 99% of CPU time dedicated to the ODE solver.
Spatiotemporal Dynamics: Simulations demonstrate the platform's ability to model the propagation of fibrotic responses. Localized stimulation (e.g., X-shaped brush selection with TGF-β) triggers intracellular cascades that lead to ECM production, which then diffuses to neighboring cells, creating spatially heterogeneous concentration gradients. The system captures the lag between signal activation and ECM accumulation, as well as the synergistic effects of combinatorial cytokine stimulation (e.g., TGF-β + AngII + IL6).
Parameter Sensitivity: Sensitivity analysis reveals that the ECM production rate ( $k_{production}$ ) is the dominant regulator of ECM remodeling outcomes, followed by activation and degradation rates. The model is robust to moderate variations in other parameters, consistent with biological understanding that matrix accumulation is primarily driven by production rates in activated fibroblasts.

Significance and Claims
The authors claim that ECMSim represents a significant advancement in accessible computational modeling for complex biological systems. By combining mathematical rigor with user-friendly web technologies, the platform addresses critical gaps in current cardiac fibrosis research, specifically the lack of spatial resolution and real-time interactivity.

The paper posits that ECMSim facilitates the investigation of pathological conditions, hypothetical scenarios, and the effects of experimental drugs on target receptors without the temporal and resource constraints of traditional experiments. It serves as a "digital laboratory" for testing hypotheses and validating models. While primarily focused on cardiac fibrosis, the authors note that the underlying reaction-diffusion framework provides a foundation for modeling other applications, including cancer micro-environments, wound healing, and regenerative biology. The platform aims to bridge the gap between molecular-scale experimental observations and tissue-level pathophysiology, fostering interdisciplinary collaboration and accelerating the exploration of novel therapeutic strategies.