Simulating biochemical reactions: The Linear Noise Approximation can capture non-linear dynamics
Abstract
There is a plethora of highly stochastic non-linear dynamical systems in fields such as molecular biology, chemistry, epidemiology, and ecology. Yet, none of the currently available stochastic models are both accurate and computationally efficient for long-term predictions of large systems. The Linear Noise Approximation (LNA) model for biochemical reaction networks is analytically tractable, which makes it computationally efficient for simulation, analysis, and inference. However, it is only accurate for linear systems and short-time transitions. Other methods can achieve greater accuracy across a wider range of systems, including non-linear ones, but lack analytical tractability. This paper seeks to challenge the prevailing view by demonstrating that the Linear Noise Approximation can indeed capture non-linear dynamics after certain modifications. We introduce a new framework that utilises centre manifold theory allowing us to identify simple interventions to the LNA that do not significantly compromise its computational efficiency. We develop specific algorithms for systems that exhibit oscillations or bi-stability and demonstrate their accuracy and computational efficiency across multiple examples.