Abstract:In recent years, live-cell imaging has generated detailed spatiotemporal datasets of biochemical networks within cells. These networks often exhibit characteristics of spatially-distributed excitable systems, with propagating waves of signaling activity that govern processes such as cell migration, division, and other essential physiological functions. Traditionally, these reaction-diffusion systems have been modeled using stochastic partial differential equations incorporating spatial Langevin-type dynamics. Although these knowledge-based models have provided valuable insights, they are typically not directly inferred from experimental data. In this study, we introduce and apply two data-driven methodologies for learning the structure and parameters of spatial Langevin equation (SLE) models: 1) Kramers-Moyal regression, an adaptation of an established approach to fit microscale stochastic dynamics, and 2) wave features optimization, a novel approach to match macroscale spatiotemporal patterns. Unlike black-box neural network models that predict system behavior but provide no mechanistic insight, this approach estimates nonlinear model equations that directly relate to system structure and dynamics. As a proof-of-concept, we focus on simulation datasets derived from two stochastic reaction-diffusion models: one based on the FitzHugh-Nagumo equations (FHN model) and another on a biochemically adapted version of the FHN equations (FR model). Our results demonstrate that optimizing stochastic dynamics and wave patterns enables the accurate estimation of SLE model structure and parameters directly from 1D and 2D spatial data. With sparsity enforcement, we also identify novel sparse reaction-diffusion models for excitable systems. We show that this approach effectively approximates system behavior even when working with datasets similar to experimental conditions, with low temporal resolution and unobserved molecular components. By leveraging machine learning techniques for robust estimation of excitable reaction-diffusion models from spatiotemporal data, this work enhances our ability to model and understand the complex systems that regulate cell behavior.

Using a library of chemical reactions to fit systems of ordinary differential equations to agent-based models: a machine learning approach

Learning differential equation models from stochastic agent-based model simulations

Data2Dynamics: a Modeling Environment Tailored to Parameter Estimation in Dynamical Systems.

Bayesian Parameter Estimation for Biochemical Reaction Networks Using Region-Based Adaptive Parallel Tempering

Data-driven model discovery and model selection for noisy biological systems

A Scalable Moment-Closure Approximation for Large-Scale Biochemical Reaction Networks

Spatial-Stochastic Simulation of Reaction-Diffusion Systems

Optimisation of Simulations of Stochastic Processes by Removal of Opposing Reactions

Machine learning sparse reaction-diffusion models from stochastic dynamics and spatiotemporal patterns

A machine learning based approach to reaction rate estimation

Atmospheric chemistry surrogate modeling with sparse identification of nonlinear dynamics

Rapid Bayesian identification of sparse nonlinear dynamics from scarce and noisy data

Learning Interactions in Reaction Diffusion Equations by Neural Networks

The SBML ODE Solver Library: a native API for symbolic and fast numerical analysis of reaction networks

Algorithmic discovery of dynamic models from infectious disease data

Reduced order dynamical models for complex dynamics in manufacturing and natural systems using machine learning

Direct ab initio MD simulation of silver ion diffusion in chalcogenide glasses

Learning black- and gray-box chemotactic PDEs/closures from agent based Monte Carlo simulation data

A general method for numerically simulating the stochastic time evolution of coupled chemical reactions

Hybrid methods in reaction-diffusion equations