Abstract:The Chemical Master Equation (CME) is a cornerstone of stochastic analysis and simulation of models of biochemical reaction networks. Yet direct solutions of the CME have remained elusive. Although several approaches overcome the infinite dimensional nature of the CME through projections or other means, a common feature of proposed approaches is their susceptibility to the curse of dimensionality, i.e. the exponential growth in memory and computational requirements in the number of problem dimensions. We present a novel approach that has the potential to "lift" this curse of dimensionality. The approach is based on the use of the recently proposed Quantized Tensor Train (QTT) formatted numerical linear algebra for the low parametric, numerical representation of tensors. The QTT decomposition admits both, algorithms for basic tensor arithmetics with complexity scaling linearly in the dimension (number of species) and sub-linearly in the mode size (maximum copy number), and a numerical tensor rounding procedure which is stable and quasi-optimal. We show how the CME can be represented in QTT format, then use the exponentially-converging hp-discontinuous Galerkin discretization in time to reduce the CME evolution problem to a set of QTT-structured linear equations to be solved at each time step using an algorithm based on Density Matrix Renormalization Group (DMRG) methods from quantum chemistry. Our method automatically adapts the "basis" of the solution at every time step guaranteeing that it is large enough to capture the dynamics of interest but no larger than necessary, as this would increase the computational complexity. Our approach is demonstrated by applying it to three different examples from systems biology: independent birth-death process, an example of enzymatic futile cycle, and a stochastic switch model. The numerical results on these examples demonstrate that the proposed QTT method achieves dramatic speedups and several orders of magnitude storage savings over direct approaches.

Rate Constant Matrix Contraction Method for Stiff Master Equations with Detailed Balance

Fast and Numerically Stable Implementation of Rate Constant Matrix Contraction Method

Method of Conditional Moments (MCM) for the Chemical Master Equation

A Kinetics-Based Method for Constraint Selection in Rate-Controlled Constrained Equilibrium

Explicit stabilized multirate method for stiff differential equations

Explicit stabilized multirate methods for the monodomain model in cardiac electrophysiology

On a positive-preserving, energy-stable numerical scheme to mass-action kinetics with detailed balance

The Equilibrium State Method for Hyperbolic Conservation Laws with Stiff Reaction Terms

A low-rank complexity reduction algorithm for the high-dimensional kinetic chemical master equation

Explicit Stabilized Multirate Method for Stiff Stochastic Differential Equations

Diffusion Limited Kinetics In Reactive Systems

Master Equation Approach to Molecular Motors

Order Reduction of the Chemical Master Equation via Balanced Realisation

Direct Solution of the Chemical Master Equation Using Quantized Tensor Trains

Krylov-based Adaptive-Rank Implicit Time Integrators for Stiff Problems with Application to Nonlinear Fokker-Planck Kinetic Models

Contraction Rate Estimates of Stochastic Gradient Kinetic Langevin Integrators

A Multiple Step Active Stiffness Integration Scheme to Couple a Stochastic Cross-Bridge Model and Continuum Mechanics for Uses in Both Basic Research and Clinical Applications of Heart Simulation

The Stabilized Explicit Variable-Load Solver with Machine Learning Acceleration for the Rapid Solution of Stiff Chemical Kinetics

Optimisation of Simulations of Stochastic Processes by Removal of Opposing Reactions