pith. sign in

arxiv: 2605.28652 · v1 · pith:DV7IG55Tnew · submitted 2026-05-27 · 🧬 q-bio.QM · nlin.CD· q-bio.PE· q-bio.SC

Widespread quasi-steady state assumption in biological interaction modeling mischaracterizes system transitions

Pan-Jun Kim This is my paper

Pith reviewed 2026-06-29 09:00 UTC · model grok-4.3

classification 🧬 q-bio.QM nlin.CDq-bio.PEq-bio.SC
keywords quasi-steady state approximationbiological modelingsystem transitionsrelaxation dynamicsoscillationsfeedback interactionscellular decision-makingmetabolic oscillations
0
0 comments X

The pith

The quasi-steady state assumption misestimates transition durations and oscillation onsets in biological systems by ignoring relaxation dynamics.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper establishes that the widely used quasi-steady state assumption in modeling biological interactions discards important relaxation effects toward those states. Near points where the system changes to a different state, these overlooked processes lead to incorrect estimates of how long transitions take and where oscillations begin. The framework shows that feedback interactions between components create time-delay effects that shift these critical points. A sympathetic reader would care because this assumption underpins models from cellular decisions to ecological cycles, so correcting it could change predicted behaviors in those systems.

Core claim

The paper claims that despite extreme slowdown near transition points, the QSSA alone misestimates transition durations and erroneously predicts the onset point for oscillations, while explicit consideration of relaxation dynamics, driven by common feedback interactions, facilitates or suppresses oscillation onset through a counterintuitive time-delay effect, as verified in simulations of cellular decision-making, metabolic oscillations, and ecological cycles.

What carries the argument

Relaxation processes toward quasi-steady states in systems with feedback interactions between components.

If this is right

  • QSSA misestimates the duration of the transition from one state to another.
  • The QSSA erroneously predicts the transition point for the onset of oscillations.
  • Relaxation dynamics facilitate or suppress the oscillation onset with a time-delay effect.
  • Common feedback interactions between biological components are pivotal to the relaxation effects.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Models of metabolic oscillations may need adjustment to account for delayed onset due to relaxation.
  • Ecological cycle predictions could shift if relaxation is included in near-transition analysis.
  • This framework might extend to other systems with fast-slow dynamics beyond biology.
  • Experimental tests could involve measuring transition timings in cellular systems under controlled conditions.

Load-bearing premise

The assumption that relaxation processes toward quasi-steady states are non-negligible near transition points and that feedback interactions drive the time-delay effects.

What would settle it

A simulation or experiment where the predicted oscillation onset point differs between the QSSA and the relaxation-inclusive model, or where measured transition duration matches one but not the other.

Figures

Figures reproduced from arXiv: 2605.28652 by Pan-Jun Kim.

Figure 1
Figure 1. Figure 1: Genetic switch and induction kinetics. (a) Protein production mechanism with positive autoregulation and induction. (b) Bifurcation diagram of the simulated protein level as a function of inducer level 𝜂 at 𝜎ൌ 0.015 (basal transcription rate). The same 𝜎 is used for (c)–(e) as well, but the observation here remains valid for other 𝜎 values (e.g., Fig. S1). The steady state is plotted as 𝜂 increases (solid … view at source ↗
Figure 2
Figure 2. Figure 2: Glycolytic enzyme activity and oscillation onset. [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Consumer cooperativity and population oscillation. (a) Time-series of the simulated, QSSA￾based consumer level (solid line) with consumer cooperativity 𝜙ൌ 1.58, 1.63, 2.1, 2.22, or 2.27 (top to bottom) at 𝑎ൌ 4.4 and 𝑏ൌ 4.25. The analytical steady state from Eqs. (18)–(20) at 𝜙ൌ𝜙୕ୡଵ (dashed line) or 𝜙ൌ𝜙୕ୡଶ (dotted line) is presented for visual guidance. (b) Bifurcation diagram of the simulated, QSSA-based c… view at source ↗
Figure 4
Figure 4. Figure 4: Time-delay and amplitude reduction effects of relaxation. [PITH_FULL_IMAGE:figures/full_fig_p016_4.png] view at source ↗
read the original abstract

From molecular, cellular, to ecological systems, the modeling of biological processes often stands on the assumption that fast components immediately reach the equilibrium at each moment (quasi-steady state) and only slow components govern the relevant system dynamics. This quasi-steady state approximation (QSSA) simplifies the modeling but discards the effects of the relaxation towards each quasi-steady state. Unclear is the QSSA's suitability around the transition point, a specific condition where the system changes to a qualitatively different state. In this regard, we here derived a theoretical framework for the near-transition dynamics of biological systems, explicitly considering the relaxation processes overlooked by the QSSA. Numerical simulations verify our predictions for cellular decision-making, metabolic oscillations, and ecological cycles. Despite the extreme slowdown near the transition point, the QSSA alone misestimates the duration of the transition from one state to another. Moreover, the QSSA erroneously predicts the transition point itself for the onset of oscillations, while the relaxation dynamics facilitates or suppresses the oscillation onset with a counterintuitive time-delay effect. Common feedback interactions between biological components are pivotal to those relaxation effects. Our study provides an analytical foundation to understand the rich transient or rhythmic dynamics of interacting biological components near the transitions.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 1 minor

Summary. The paper claims that the quasi-steady state assumption (QSSA) commonly used in biological modeling discards relaxation processes to quasi-steady states and thereby mischaracterizes transitions, including misestimating their duration and incorrectly predicting the onset of oscillations. The authors derive an analytical framework for near-transition dynamics that retains these relaxation effects, show that common feedback interactions produce counterintuitive time-delay phenomena, and verify the predictions via numerical simulations in cellular decision-making, metabolic oscillations, and ecological cycles.

Significance. If the central claims hold under clearly stated conditions, the work would be significant for systems biology and ecology because the QSSA is foundational; identifying its systematic errors near transitions could alter how transients and rhythms are modeled. The explicit inclusion of relaxation dynamics and the multi-domain simulation verification constitute concrete strengths that allow direct testing of the framework.

major comments (1)
  1. [theoretical framework (near-transition dynamics derivation)] The central claim that relaxation toward quasi-steady states remains relevant and produces measurable shifts in transition duration and oscillation threshold even near the QSSA breakdown point requires an explicit inequality relating the fast relaxation eigenvalue, the distance to the transition, and the feedback strength. No such condition is supplied in the theoretical framework, leaving open whether the claimed mischaracterization occurs in the parameter regimes of the cellular, metabolic, and ecological examples.
minor comments (1)
  1. [Abstract] The abstract states that 'numerical simulations verify our predictions' but does not indicate the range of parameters or initial conditions explored; a brief statement of the tested regimes would improve clarity.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the constructive and detailed report. We address the single major comment below and will revise the manuscript to incorporate an explicit condition as requested.

read point-by-point responses

  1. Referee: The central claim that relaxation toward quasi-steady states remains relevant and produces measurable shifts in transition duration and oscillation threshold even near the QSSA breakdown point requires an explicit inequality relating the fast relaxation eigenvalue, the distance to the transition, and the feedback strength. No such condition is supplied in the theoretical framework, leaving open whether the claimed mischaracterization occurs in the parameter regimes of the cellular, metabolic, and ecological examples.

    Authors: We agree that an explicit inequality would strengthen the presentation of the regime of validity. Although the near-transition derivation already encodes the relevant scaling through the fast eigenvalue and the distance-to-transition parameter, we acknowledge that a standalone inequality was not stated. In the revised manuscript we will derive and insert the condition |λ_fast| ≫ |δ|·β (with λ_fast the fast relaxation eigenvalue, δ the distance to the transition, and β the feedback strength) and will verify that the three example parameter sets satisfy it, thereby confirming that the reported mischaracterization of transition duration and oscillation onset is indeed operative in those regimes. revision: yes

Circularity Check

0 steps flagged

No circularity; framework derived independently then verified by simulation

full rationale

The provided abstract states that a theoretical framework is derived first for near-transition dynamics explicitly including relaxation processes, after which numerical simulations are used to verify the predictions for cellular, metabolic, and ecological examples. No equations, parameter fits, or self-citations appear in the text that would reduce any claimed prediction to an input by construction. The central assertions (QSSA misestimates transition duration and oscillation onset, with relaxation producing time-delay effects) are presented as outputs of the derivation rather than tautological restatements of fitted quantities or prior self-referential results. This matches the default expectation of a self-contained derivation chain.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

No free parameters, axioms, or invented entities are specified in the abstract; the framework is described as a theoretical derivation verified by simulations.

pith-pipeline@v0.9.1-grok · 5756 in / 1069 out tokens · 26126 ms · 2026-06-29T09:00:34.611953+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

53 extracted references

  1. [1]

    Held, Simplicity amid complexity, Science 343, 1206–1207 (2014)

    I. Held, Simplicity amid complexity, Science 343, 1206–1207 (2014)

    2014

  2. [2]

    Born and R

    M. Born and R. Oppenheimer, Zur Quantentheorie der Molekeln, Ann. Phys. 389, 457–484 (1927)

    1927

  3. [3]

    G. E. Briggs and J. B. S. Haldane, A note on the kinetics of enzyme action , Biochem. J. 19, 338–339 (1925)

    1925

  4. [4]

    K. J. Laidler, Theory of the transient phase in kine tics, with special reference to enzyme systems, Can. J. Chem. 33, 1614–1624 (1955)

    1955

  5. [5]

    Kuramoto, Chemical Oscillations, Waves, and Turbulence (Springer-Verlag, Berlin, 1984)

    Y. Kuramoto, Chemical Oscillations, Waves, and Turbulence (Springer-Verlag, Berlin, 1984)

    1984

  6. [6]

    J. D. Murrary, Mathematical Biology. I. An Introduction (Springer-Verlag, New York, 2002)

    2002

  7. [7]

    Momeni, L

    B. Momeni, L. Xie, and W. Shou, Lotka–Volterra pairwise modeling fails to capture diverse pairwise microbial interactions, eLife 6, e25051 (2017)

    2017

  8. [8]

    Kim, T.-W

    P.-J. Kim, T.-W. Ko, H. Jeong, K. J. Lee, and S. K. Han, Emergence of chaotic itinerancy in simple ecological systems, Phys. Rev. E 76, 065201(R) (2007)

    2007

  9. [9]

    M. W. Covert, E. M. Knight, J. L. Reed, M. J. Herrgård, and B. Ø. Palsson, Integrating high- throughput and computational data elucidates bacterial networks, Nature 429, 92–96 (2004)

    2004

  10. [10]

    Haken, Slaving principle revisited, Physica D 97, 95–103 (1996)

    H. Haken, Slaving principle revisited, Physica D 97, 95–103 (1996)

    1996

  11. [11]

    C. M. Bender and S. A. Orszag, Advanced Mathematicai Methods for Scientists and Engineers (McGraw-Hill, New York, 1978)

    1978

  12. [12]

    J. K. Kim and J. J. Tyson, Misuse of the Michaelis–Menten rate law for protein interaction networks and its remedy , PLoS Comput. Biol. 16, e1008258 (2020)

    2020

  13. [13]

    L. A. Segel and M. Slemrod, The quasi-steady-state assu mption: A case study in perturbation , SIAM Rev. 31, 446–477 (1989)

    1989

  14. [14]

    F. G. Heineken, H. M. Tsuchiya, and R. Aris, On the mathematical status of the pseudo- steady state hypothesis of biochemical kinetics, Math. Biosci. 1, 95–113 (1967)

    1967

  15. [15]

    A. N. Tikhonov, Системы дифференциальных уравнений, содержащие малые параметры при производных, Mat. Sb. 31, 575–586 (1952). 20

    1952

  16. [16]

    A. B. Vasil’eva, Asymptotic behaviour of solutions to certain problems involving non-linear differential equations containing a small parameter multiplying the highest derivatives, Russ. Math. Surv. 18, 13–84 (1963)

    1963

  17. [17]

    D. R. Smith, Singular Perturbation Theory (Cambridge University Press, Cambridge, 1985)

    1985

  18. [18]

    L.-Y. Chen, N. Goldenfeld, and Y. Oono, Renormalization group and singular perturbations: Multiple scales, boundary layers, and reductive perturbation theory, Phys. Rev. E 54, 376–394 (1996)

    1996

  19. [19]

    K. C. Chen, A. Csikasz-Nagy, B. Gyorffy, J. Val, B. Novak, and J. J. Tyson, Kinetic analysis of a molecular model of the budding yeast cell cycle , Mol. Biol. Cell 11, 369–391 (2000)

    2000

  20. [20]

    J. K. Kim and D. B. Forger, A mechanism for robust circadian timekeeping via stoichiometric balance , Mol. Syst. Biol. 8, 630 (2012)

    2012

  21. [21]

    M. Foo, D. E. Somers, and P.-J. Kim, Kernel architecture of the genetic circuitry of the Arabidopsis circadian system , PLoS Comput. Biol. 12, e1004748 (2016)

    2016

  22. [22]

    Lechón-Alonso, S

    P. Lechón-Alonso, S. Kundu, P. Lemos-Costa, J. A. Capitán, and S. Allesina, Robust coexistence in competitive ecological communities , Nat. Commun. 17, 2637 (2026)

    2026

  23. [23]

    J. E. Purvis and G. Lahav, Encoding and decoding cellular information through signaling dynamics , Cell 152, 945–956 (2013)

    2013

  24. [24]

    Betz and B

    A. Betz and B. Chance, Phase relationship of glycolytic intermediates in yeast cells with oscillatory metabolic control, Arch. Biochem. Biophys. 109, 585–594 (1965)

    1965

  25. [25]

    S. H. Strogatz, Nonlinear Dynamics and Chaos: With Applications to Physics, Biology, Chemistry, and Engineering (Addison-Wesley, Reading, Massachusetts, 1994)

    1994

  26. [26]

    Wiggins, Introduction to Applied Nonlinear Dynamical Systems and Chaos (Springer- Verlag, New York, 1990)

    S. Wiggins, Introduction to Applied Nonlinear Dynamical Systems and Chaos (Springer- Verlag, New York, 1990)

    1990

  27. [27]

    Nazarimehr, S

    F. Nazarimehr, S. Jafari, M. Perc, and J. C. Sprott, Critical slowing down indicators, Europhys. Lett. 132, 18001 (2020)

    2020

  28. [28]

    M. I. Maturana, C. Meisel, K. Dell, P. J. Karoly, W. D’Souza, D. B. Grayden, A. N. Burkitt, P. Jiruska, J. Kudlacek, J. Hlinka et al., Critical slowing down as a biomarker for seizure susceptibility, Nat. Commun. 11, 1–12 (2020)

    2020

  29. [29]

    Scheffer, Foreseeing tipping points, Nature 467, 411–412 (2010)

    M. Scheffer, Foreseeing tipping points, Nature 467, 411–412 (2010)

    2010

  30. [30]

    J. Chae, R. Lim, T. L. P. Martin, C.-M. Ghim, and P.-J. Kim, Enlightening the blind spot of the Michaelis–Menten rate law: The role of relaxation dynamics in molecular complex formation, J. Theor. Biol. 597, 111989 (2025)

    2025

  31. [31]

    J. R. Reitz, F. J. Milford, and R. W. Christy, Foundations of Electromagnetic Theory (Addison-Wesley, Reading, Massachusetts, 1979). 21

    1979

  32. [32]

    J. K. Hale and S. M. Verduyn Lunel, Introduction to Functional Differential Equations (Springer-Verlag, New York, 1993)

    1993

  33. [33]

    T. J. Perkins and P. S. Swain, Strategies for cellular decision-making , Mol. Syst. Biol. 5, 326 (2009)

    2009

  34. [34]

    R. Lim, T. L. P. Martin, J. Chae, W. J. Kim, C.-M. Ghim, and P.-J. Kim, Generalized Michaelis–Menten rate law with time-varying molecular concentrations , PLoS Comput. Biol. 19, e1011711 (2023)

    2023

  35. [35]

    Y. T. Maeda and M. Sano, Regulatory dynamics of synthetic gene networks with positive feedback, J. Mol. Biol. 359, 1107–1124 (2006)

    2006

  36. [36]

    Jeynes-Smith and R

    C. Jeynes-Smith and R. P. Araujo, Ultrasensitivity and bistability in covalent-modification cycles with positive autoregulation, Proc. R. Soc. A 477, 20210069 (2021)

    2021

  37. [37]

    E. K. Pye, Biochemical mechanisms underlying the metabolic oscillations in yeast, Can. J. Botany 47, 271–285 (1969)

    1969

  38. [38]

    Gustavsson, D

    A.-K. Gustavsson, D. D. van Niekerk, C. B. Adiels, B. Kooi, M. Goksör, and J. L. Snoep, Allosteric regulation of phosphofructokinase controls the emergence of glycolytic oscillations in isolated yeast cells, FEBS J. 281, 2784–2793 (2014)

    2014

  39. [39]

    E. E. Sel’kov, Self-oscillations in glycolysis. 1. A simple kinetic model, European J. Biochem. 4, 79–86 (1968)

    1968

  40. [40]

    Cezairliyan and F

    B. Cezairliyan and F. M. Ausubel, Investment in secreted enzymes during nutrient-limited growth is utility dependent, Proc. Natl. Acad. Sci. USA 114, E7796–E7802 (2017)

    2017

  41. [41]

    T. D. Seeley and P. K. Visscher, Assessing the benefits of cooperation in honeybee foraging: Search costs, forage quality, and competitive ability, Behav. Ecol. Sociobiol. 22, 229–237 (1988)

    1988

  42. [42]

    D. W. Macdonald, S. Creel, and M. G. L. Mills, in The Biology and Conservation of Wild Canids, edited by D. W. Macdonald and C. Sillero-Zubiri (Oxford University Press, New York, 2004), pp. 85–106

    2004

  43. [43]

    Ordaz-Rivas, A

    E. Ordaz-Rivas, A. Rodriguez-Liñan, and L. Torres-Treviño, Flock of robots with self- cooperation for prey–predator task, J. Intel. Robot. Syst. 101, 1–16 (2021)

    2021

  44. [44]

    D. V. Ramana Reddy, A. Sen, and G. L. Johnston, Time delay induced death in coupled limit cycle oscillators , Phys. Rev. Lett. 80, 5109–5112 (1998)

    1998

  45. [45]

    F. M. Atay, Distributed delays facilitate amplitude death of coupled oscillators , Phys. Rev. Lett. 91, 094101 (2003)

    2003

  46. [46]

    T. Song, H. C. Park, and K.-H. Ahn, Proposal for high sensitivit y force sensor inspired by auditory hair cells , Appl. Phys. Lett. 95, 013702 (2009). 22

    2009

  47. [47]

    T. Song, W. S. Lee, and K.-H. Ahn, Physical limits to auditory transduction of hair-cell bundles probed by a biomimetic system , Sci. Rep. 5, 11470 (2015)

    2015

  48. [48]

    A. K. Manrai and J. Gunawardena, The geometry of multisite phosphorylation , Biophys. J. 95, 5533–5543 (2008)

    2008

  49. [49]

    Kim and N

    P.-J. Kim and N. D. Price, Macroscopic kinetic effect of cell-to-cell variation in biochemical reactions , Phys. Rev. Lett. 104, 148103 (2010)

    2010

  50. [50]

    A" to "I

    K. S. Brown and J. P. Sethna, Statistical mechanical approaches to models with many poorly known parameters , Phys. Rev. E 68, 021904 (2003). 1 Supplementary Information (SI) for Widespread quasi-steady state assumption in biological interaction modeling mischaracterizes system transitions Pan-Jun Kim Technical details of numerical simulation and statisti...

    2003

  51. [51]

    2(b)–(d), S2(a)–(c) and Section J)

    Our numerical simulations show the presence of the supercritical Hopf bifurcation, whether based on the original QSSA model or the current extended model with 𝜇୭ିଵ൐0 (Figs. 2(b)–(d), S2(a)–(c) and Section J). From Eqs. (12), (13), (E1), and (E2), 𝐺ଵ൫𝑥ሺ𝑡ሻ, 𝑦ଵሺ𝑡ሻ, 𝑦ଶሺ𝑡ሻ൯ൌ1 െ𝑥ሺ𝑡ሻ, (H18) 𝐺ଶ൫𝑥ሺ𝑡ሻ, 𝑦ଵሺ𝑡ሻ, 𝑦ଶሺ𝑡ሻ൯ൌ𝛼൫𝑥ሺ𝑡ሻെ𝑦ଶሺ𝑡ሻ൯. (H19) At the fixed point ൫𝑥ሺ𝑡ሻ, 𝑦ଵ...

  52. [52]

    (H1)–(H3) to Eqs

    involves the partial, linear reduction of Eqs. (H1)–(H3) to Eqs. (H6)–(H8) (Section H). Therefore, the initial conditions for the simulation were chosen near the fixed point, ensuring the physically sensible, non-negative solutions. Population cycle with cooperative foraging For each parameter set, we examined the presence of sustained oscillations (limit...

  53. [53]

    (I3) (Section I)

    involves the linear reduction of 𝑋ሺ𝑡ሻ in Eq. (I3) (Section I). Therefore, the initial conditions for the simulation were chosen near the fixed point, ensuring the physically sensible, non-negative solutions. 27 Eq. (A7) captures the relaxation’s time-delay and amplitude reduction effects on 𝑥ሺ𝑡ሻ. To examine their contributions to the bifurcation point shi...

    1979