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.

Physics-based machine learning for modeling stochastic IP3-dependent calcium dynamics

Physics-informed machine learning for automatic model reduction in chemical reaction networks

An efficient reduced-lattice model of IP3R for probing Ca2+ dynamics

Discrete stochastic modeling of calcium channel dynamics

Learning Moment Closure in Reaction-Diffusion Systems with Spatial Dynamic Boltzmann Distributions

Physics-aware, probabilistic model order reduction with guaranteed stability

Physics-guided machine learning from simulated data with different physical parameters

Particle-Based Multiscale Modeling of Calcium Puff Dynamics

Encoding physics to learn reaction–diffusion processes

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

Mesoscopic behavior from microscopic Markov dynamics and its application to calcium release channels

Physics-informed machine learning

Evolutionary probability density reconstruction of stochastic dynamic responses based on physics-aided deep learning

A nonintrusive hybrid neural-physics modeling of incomplete dynamical systems: Lorenz equations

Physics-constrained neural differential equations for learning multi-ionic transport

Chemistry-informed Machine Learning Explains Calcium-binding Proteins Fuzzy Shape for Communicating Changes in the Atomic States of Calcium Ions

Mathematical model for IP$_{3}$ dependent calcium oscillations and mitochondrial associate membranes in non-excitable cells

Law of Mass Action, Detailed Balance, and the Modeling of Calcium Puffs

Machine learning in biological physics: From biomolecular prediction to design

Physics-guided Deep Markov Models for Learning Nonlinear Dynamical Systems with Uncertainty

Towards Model Discovery Using Domain Decomposition and PINNs