Coordinate-independent model reductions of chemical reaction networks based on geometric singular perturbation theory
Abstract
The quasi-steady-state approximation (QSSA) is a standard technique for reducing the complexity of chemical reaction networks (CRNs). The validity of any QSSA-based model is restricted to specific parameter regimes. Selecting the appropriate reduction is not always straightforward. At times, QSSAs are misused outside of their validity regions and, even when a particular QSSA is considered valid in a given parameter regime, other QSSAs may be simultaneously valid, creating ambiguity. Here, we employ a more powerful alternative: a constructive model reduction framework based on coordinate-independent geometric singular perturbation theory (ci-GSPT) and the parametrization method. A key advantage of this approach is its ability to derive reduced models independent of a clear timescale separation in the variables for a specific parameter configuration. We demonstrate our approach on two benchmark systems. For the Michaelis-Menten (MM) reaction, we show that the framework provides a systematic approach by exploring parameter configurations across three orders of magnitude: asymptotically large, small, and `order one'. A consequence of this systematic analysis is a geometric classification that categorizes the resulting model reductions and provides a point of comparison between our approach and common QSSA variants in the literature. For the more complex Kim-Forger model, we show that this approach successfully produces a reduction without the need for a coordinate transformation, showcasing its applicability to larger systems.
Turn this paper into a full lesson
ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.