title: Layered decomposition for the model order reduction of timescale separated biochemical reaction networks.
abstract: Biochemical reaction networks tend to exhibit behaviour on more than one timescale and they are inevitably modelled by stiff systems of ordinary differential equations. Singular perturbation is a well-established method for approximating stiff systems at a given timescale. Standard applications of singular perturbation partition the state variable into fast and slow modules and assume a quasi-steady state behaviour in the fast module. In biochemical reaction networks, many reactants may take part in both fast and slow reactions; it is not necessarily the case that the reactants themselves are fast or slow. Transformations of the state space are often required in order to create fast and slow modules, which thus no longer model the original species concentrations. This paper introduces a layered decomposition, which is a natural choice when reaction speeds are separated in scale. The new framework ensures that model reduction can be carried out without seeking state space transformations, and that the effect of the fast dynamics on the slow timescale can be described directly in terms of the original species.
authors.affiliation: Department of Engineering Science, University of Oxford, Parks Road, Oxford OX1 3PJ, United Kingdom; Life Sciences Interface Doctoral Training Centre, University of Oxford, Parks Road, Oxford OX1 3QU, United Kingdom.
authors.forename: Thomas P
authors.initials: TP
authors.lastname: Prescott
authors.affiliation: Department of Engineering Science, University of Oxford, Parks Road, Oxford OX1 3PJ, United Kingdom. Electronic address: antonis@eng.ox.ac.uk.
authors.forename: Antonis
authors.initials: A
authors.lastname: Papachristodoulou
mesh_terms: D008954:Models, Biological
keywords: Nonlinear ODEs; Reaction kinetics; Singular perturbation; Wnt pathway