Stochastic Analysis of Chemical Reaction Networks Using Linear Noise Approximation

Abstract

Stochastic evolution of Chemical Reactions Networks (CRNs) over time is usually analysed through solving the Chemical Master Equation (CME) or performing extensive simulations. Analysing stochasticity is often needed, particularly when some molecules occur in low numbers. Unfortunately, both approaches become infeasible if the system is complex and/or it cannot be ensured that initial populations are small. We develop a probabilistic logic for CRNs that enables stochastic analysis of the evolution of populations of molecular species. We present an approximate model checking algorithm based on the Linear Noise Approximation (LNA) of the CME, whose computational complexity is independent of the population size of each species and polynomial in the number of different species. The algorithm requires the solution of first order polynomial differential equations. We prove that our approach is valid for any CRN close enough to the thermodynamical limit. However, we show on four case studies that it can still provide good approximation even for low molecule counts. Our approach enables rigorous analysis of CRNs that are not analyzable by solving the CME, but are far from the deterministic limit. Moreover, it can be used for a fast approximate stochastic characterization of a CRN.

Using Spatial Partitioning to Reduce the Bit Error Rate of Diffusion-Based Molecular Communications

“…However, for systems with nonlinear reaction rates, which is the case for our system, it is not possible to solve the RDME analytically because of the moment closure problem [55] which refers to the issue that the computation of the lower order moments requires the values of the higher order moments. An alternative method is the Linear Noise Approximation (LNA) [30], [56] which uses Gaussian distribution to approximate the distribution of the molecular counts when the system is in the steady state. The LNA is known to be a good approximation when the number of the molecules in the system is large.…”

“…We can use similar method to write down the SDEs for other species in the system. With this set of SDEs, we can compute the covariance matrix for vectored Gaussian distribution of the species concentrations by solving a Lyapunov differential equation [30], [56].…”

Bioinformatics Challenges and Potentialities in Studying Extreme Environments

“…by using a network fusion approach, where layers can be weighted and the multi-layer network can be fused to a single-layer network. Finally, although not covered here, another important approach to study metabolic networks is stochastic simulation based on an approach pioneered by Gillespie [33], and relying on molecular counts to simulate the evolution of populations of chemical species [34].…”

Formal Quantitative Analysis of Reaction Networks Using Chemical Organisation Theory

“…We build on their notion of discrete organisation but focus on quantitative analysis of the transitive dynamics among the organisations, which was not considered in [14]. Other approaches for approximate analysis of discrete models of reaction networks include the use of Linear Noise Approximation [4], the Central Limit Approximation [3] and "sliding window" abstractions [17]. the system and a specification of one or more required properties of that system, normally in temporal logic (such as CTL or LTL).…”