arxiv: 2604.15716 · v1 · submitted 2026-04-17 · 🧮 math.DS · q-bio.MN

Recognition: unknown

Mathematical modeling of biochemical signal propagation in many-stage enzymatic pathways

Chathranee Jayathilaka , Mark B. Flegg

Authors on Pith no claims yet

Pith reviewed 2026-05-10 08:09 UTC · model grok-4.3

classification 🧮 math.DS q-bio.MN

keywords travelling wavesenzymatic pathwaysMichaelis-Menten kineticsbiochemical signalingspatial rescalingwavefront propagationheterogeneous networksmodel reduction

0 comments

The pith

A reciprocal-velocity spatial rescaling absorbs kinetic variations to restore uniform wave propagation in heterogeneous enzymatic pathways.

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

The paper develops a framework for travelling waves in feed-forward enzymatic signaling cascades governed by nonlinear Michaelis-Menten kinetics. In uniform pathways it shows that activation bias between nodes acts as the bifurcation parameter that determines whether stable waves exist. Heterogeneous kinetics cause wavefront distortion and fluctuating propagation speeds. A reciprocal-velocity rescaling of the spatial coordinate is introduced to compensate for these local differences, smoothing velocities and preserving profile shape without parameter tuning or continuum approximations. Extreme variations are shown to produce functional fragmentation that justifies model reduction.

Core claim

In canonical feed-forward enzymatic pathways, travelling waves exist when activation bias exceeds a threshold value, and a reciprocal-velocity spatial rescaling can be applied to absorb local kinetic heterogeneity so that wave velocity becomes constant and wavefront profiles are preserved in the transformed coordinate.

What carries the argument

The reciprocal-velocity spatial rescaling, a coordinate transformation in which space is measured inversely to the local wave velocity, which normalizes propagation across varying kinetic parameters.

If this is right

Activation bias between connected nodes functions as the bifurcation parameter that dictates the existence of travelling waves in uniform pathways.
Parameter gradients and random kinetic variations distort wavefronts and produce heavy fluctuations in propagation speed.
The rescaling technique recovers predictable signal transmission by smoothing velocities and preserving profiles without bespoke tuning.
Severe kinetic bottlenecks produce functional pathway fragmentation that supplies a mathematically justified basis for model reduction.

Where Pith is reading between the lines

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

The rescaling could be used to simplify analysis of noisy enzyme kinetics in actual cellular signaling networks.
Fragmentation thresholds might indicate points where signaling fails in disease states or where interventions could be targeted.
The same transformation may apply to wave problems in other heterogeneous media such as reaction-diffusion systems with spatially varying rates.

Load-bearing premise

The pathways consist of strictly feed-forward chains with nonlinear Michaelis-Menten kinetics in which activation bias between nodes is the dominant control on wave formation.

What would settle it

Numerical integration of the pathway equations with randomly drawn kinetic parameters at each node, verifying whether wave speed becomes constant and wavefront shape is preserved after the reciprocal-velocity coordinate change.

read the original abstract

Biochemical signalling cascades transduce extracellular stimuli into cellular responses through sequences of discrete, node-to-node activations. While signal fidelity depends critically on local interaction kinetics, the mechanisms governing information propagation in realistic, highly variable kinetic contexts remain poorly understood. In this paper, we develop a mathematical framework for travelling waves in canonical feed-forward pathways governed by nonlinear Michaelis-Menten-type kinetics. For uniform pathways, we characterise the complete steady-state landscape and demonstrate that activation bias (the contribution of the binary states of each node to downstream activation) between connected nodes acts as a key bifurcation parameter dictating wave existence. Extending this framework to heterogeneous networks, we show how parameter gradients and random kinetic variations distort wavefronts and induce heavy fluctuations in propagation speed. To recover predictable signal transmission, we introduce a novel reciprocal-velocity spatial rescaling technique. We demonstrate that this coordinate transformation inherently absorbs local kinetic variations, effectively smoothing wave velocities and preserving wavefront profiles without requiring bespoke parameter tuning or continuous limits. Finally, by testing the framework's limits against extreme parameter variability, we reveal how severe kinetic bottlenecks lead to functional pathway fragmentation, offering a mathematically justified basis for rational model reduction in complex biochemical networks.

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.

Desk Editor's Note private letter to a colleague

The reciprocal-velocity rescaling is a practical coordinate trick for smoothing waves in heterogeneous enzymatic cascades, but its exactness on finite discrete chains is the part that needs the closest look.

read the letter

The paper's main new piece is the reciprocal-velocity rescaling: they turn the node index into a continuous coordinate weighted by 1 over the local wave speed from each node's Michaelis-Menten parameters. This absorbs kinetic differences so the wavefront shape stays roughly the same even when rates vary along the chain. For the uniform case they give a full steady-state picture and show that activation bias between nodes is the parameter that decides whether a traveling wave exists at all. That part is clean and useful for anyone who already works with these cascades. They also map out how random or gradient variations distort speed and shape, then test what happens at extreme bottlenecks where the pathway can fragment. Those sections give a concrete basis for deciding when a reduced model is still safe. The soft spot is the claim that the rescaling works exactly without any continuous-limit assumption. On a finite discrete chain the new coordinate is really an interpolation, so the transformed ODEs only decouple into a clean traveling-wave equation if the variations are slow or the number of stages is large. The abstract asserts it holds for arbitrary variations, but the stress-test note is right to flag that this needs an explicit justification or counter-example in the full text; if the derivations only hold asymptotically, the practical payoff shrinks. The math looks formally grounded where they stay inside the uniform or slowly varying regime, and the citations to prior wave work in signaling are standard rather than cherry-picked. This is for readers who build or simplify multi-stage biochemical models and want a quick way to handle heterogeneity without fitting every parameter. It is not going to change the field, but the framework is concrete enough that a serious referee could check the derivations and the numerical tests in a few hours. I would send it to peer review.

Referee Report

1 major / 2 minor

Summary. The paper develops a dynamical-systems framework for traveling waves in feed-forward enzymatic cascades governed by Michaelis-Menten kinetics. For uniform pathways it characterizes the full steady-state landscape and identifies activation bias as the bifurcation parameter controlling wave existence. For heterogeneous pathways it introduces a reciprocal-velocity spatial rescaling that is claimed to absorb arbitrary local kinetic variations, smooth propagation speeds, and preserve wavefront profiles exactly in the discrete finite-N chain without parameter tuning or continuous limits. The work concludes by showing how extreme bottlenecks produce functional fragmentation, thereby justifying rational model reduction.

Significance. If the rescaling is shown to be exact for finite discrete chains, the framework supplies a parameter-free coordinate transformation that converts heterogeneous kinetic networks into effectively uniform traveling-wave problems. This would be a useful addition to the literature on signal propagation in biochemical networks and could guide model reduction when kinetic data are incomplete.

major comments (1)

[heterogeneous-network analysis / rescaling section] The central claim of the reciprocal-velocity rescaling (abstract and the section introducing the heterogeneous-network analysis) is that the transformation “inherently absorbs local kinetic variations … without requiring … continuous limits.” For a finite discrete chain the rescaled node index is an interpolation whose local metric is 1/v_i (v_i obtained from each node’s Michaelis-Menten parameters). The transformed system of ODEs therefore decouples into a uniform traveling-wave equation only if the variations satisfy an unstated slow-variation or large-N condition. Please supply the explicit transformed equations and demonstrate that exact preservation holds for arbitrary finite-N random draws without additional assumptions.

minor comments (2)

[abstract] The abstract states that activation bias “acts as a key bifurcation parameter” but does not specify the precise mathematical definition (e.g., the functional form relating binary node states to downstream activation probability). A short clarifying sentence would help readers before the uniform-pathway analysis.
[throughout] Notation for the local velocity v_local and the rescaled coordinate should be introduced once and used consistently; at present the same symbol appears to be overloaded between the uniform and heterogeneous cases.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading and constructive feedback on our manuscript. We address the single major comment below with a detailed explanation of the rescaling construction and commit to revisions that supply the requested explicit equations and demonstrations.

read point-by-point responses

Referee: [heterogeneous-network analysis / rescaling section] The central claim of the reciprocal-velocity rescaling (abstract and the section introducing the heterogeneous-network analysis) is that the transformation “inherently absorbs local kinetic variations … without requiring … continuous limits.” For a finite discrete chain the rescaled node index is an interpolation whose local metric is 1/v_i (v_i obtained from each node’s Michaelis-Menten parameters). The transformed system of ODEs therefore decouples into a uniform traveling-wave equation only if the variations satisfy an unstated slow-variation or large-N condition. Please supply the explicit transformed equations and demonstrate that exact preservation holds for arbitrary finite-N random draws without additional assumptions.

Authors: We appreciate the referee drawing attention to the need for explicit detail. The reciprocal-velocity rescaling is defined directly on the discrete chain by the cumulative coordinate s_i = sum_{j=1}^i (1/v_j), where v_j > 0 is the local characteristic speed obtained from the Michaelis-Menten parameters at node j (explicitly, v_j is the speed that would be observed if the entire pathway were uniform with those parameters). The traveling-wave ansatz is then imposed in the rescaled coordinate: the state vector at node i is a function of (s_i - t) with unit speed. Substituting this ansatz into the original system of N coupled ODEs and using the chain rule for the discrete differences yields an equation for the profile that is identical to the uniform-case traveling-wave ODE, because the local speed variations have been absorbed into the nonuniform spacing of the s_i points. The derivation requires no continuum limit, no slow-variation assumption, and holds for arbitrary finite N and arbitrary positive v_i (including random draws). We will add the full set of transformed equations, the step-by-step substitution, and the proof of equivalence to the revised manuscript. We will also include numerical verification for several finite-N random-parameter realizations (N = 5, 10, 20) confirming that the wavefront shape is preserved exactly once the coordinate is rescaled. revision: yes

Circularity Check

1 steps flagged

Reciprocal-velocity rescaling absorbs variations by construction

specific steps

self definitional [Abstract]
"we introduce a novel reciprocal-velocity spatial rescaling technique. We demonstrate that this coordinate transformation inherently absorbs local kinetic variations, effectively smoothing wave velocities and preserving wavefront profiles without requiring bespoke parameter tuning or continuous limits."

The rescaling is defined by setting the new spatial metric to 1/v_local where v_local is computed directly from each node's Michaelis-Menten kinetics. By this definition the transformed propagation speed is forced to uniformity, so the 'inherent absorption' and profile preservation are equivalent to the input choice of coordinate rather than a derived property of the finite discrete system.

full rationale

The paper's central novel technique redefines the spatial coordinate using the reciprocal of locally computed wave velocities derived from the model's own Michaelis-Menten parameters. This makes the claimed absorption of kinetic variations and preservation of wavefront profiles hold tautologically in the transformed system, reducing the result to a reparameterization rather than an independent derivation from the discrete ODE chain. The uniform-pathway bifurcation analysis stands separately, but the heterogeneous extension relies on this definitional step. No other circular reductions were identified in the provided derivation chain.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The framework rests on standard Michaelis-Menten kinetics for enzymatic reactions and the assumption of feed-forward cascade structure; no new entities are postulated.

free parameters (1)

activation bias
Treated as the key bifurcation parameter whose value dictates wave existence; its specific values are not derived from first principles in the abstract.

axioms (1)

domain assumption Canonical feed-forward pathways are governed by nonlinear Michaelis-Menten-type kinetics
Explicitly stated as the governing kinetics for the traveling-wave analysis.

pith-pipeline@v0.9.0 · 5507 in / 1256 out tokens · 26498 ms · 2026-05-10T08:09:22.590306+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

20 extracted references · 3 canonical work pages

[1]

Cell100(1), 113–127 (2000)

Hunter, T.: Signaling—2000 and beyond. Cell100(1), 113–127 (2000)

2000
[2]

In: Othmer, H.G., Adler, F.R., Lewis, M.A., Dallon, J.C

Othmer, H.G.: Signal transduction and second messenger systems. In: Othmer, H.G., Adler, F.R., Lewis, M.A., Dallon, J.C. (eds.) Case Studies in Mathematical Modeling: Ecology, Phys- iology, and Cell Biology, pp. 101–128. Prentice Hall, Englewood Cliffs, NJ (1997). Classic mathematical treatment of signal transduction dynamics with multiple qualitative behaviors

1997
[3]

Bio-protocol11(14), 4089 (2021) https://doi.org/10.21769/BioProtoc.4089

Schwen, P., Geiger, T.: Modeling the nonlinear dynamics of intracellular signaling networks. Bio-protocol11(14), 4089 (2021) https://doi.org/10.21769/BioProtoc.4089 . Protocol on non- linear ODE and PDE models of intracellular signaling dynamics, including multistability and oscillations

work page doi:10.21769/bioprotoc.4089 2021
[4]

Biophysical journal52(5), 807–828 (1987)

Martiel, J.-L., Goldbeter, A.: A model based on receptor desensitization for cyclic amp signaling in dictyostelium cells. Biophysical journal52(5), 807–828 (1987)

1987
[5]

Journal of mathematical biology 79(5), 1831–1883 (2019)

Hurtado, P.J., Kirosingh, A.S.: Generalizations of the ‘linear chain trick’: incorporating more flexible dwell time distributions into mean field ode models. Journal of mathematical biology 79(5), 1831–1883 (2019)

2019
[6]

Nature Communications12(1), 1788 (2021) https://doi

Glass, D.S., Jang, S.H., Fraser, H.B.: Nonlinear delay differential equations and their application to modeling biological network motifs. Nature Communications12(1), 1788 (2021) https://doi. org/10.1038/s41467-021-21700-8

work page doi:10.1038/s41467-021-21700-8 2021
[7]

I: An Introduction, 3rd edn

Murray, J.D.: Mathematical Biology. I: An Introduction, 3rd edn. Interdisciplinary Applied Mathematics, vol. 17. Springer, New York (2002)

2002
[8]

Proceedings of the National Academy of Sciences of the United States of America92(11), 5042–5046 (1995) https://doi.org/10.1073/pnas.92.11.5042

McKeithan, T.W.: Kinetic proofreading in t-cell receptor signal transduction. Proceedings of the National Academy of Sciences of the United States of America92(11), 5042–5046 (1995) https://doi.org/10.1073/pnas.92.11.5042

work page doi:10.1073/pnas.92.11.5042 1995
[9]

Endocrine reviews22(2), 153–183 (2001)

Pearson, G., Robinson, F., Beers Gibson, T., Xu, B.-e., Karandikar, M., Berman, K., Cobb, M.H.: Mitogen-activated protein (map) kinase pathways: regulation and physiological functions. Endocrine reviews22(2), 153–183 (2001)

2001
[10]

Cell149(6), 1192–1205 (2012)

Clevers, H., Nusse, R.: Wnt/β-catenin signaling and disease. Cell149(6), 1192–1205 (2012)

2012
[11]

Proceedings of the National Academy of Sciences93(19), 10078–10083 (1996)

Huang, C.-Y., Ferrell Jr, J.E.: Ultrasensitivity in the mitogen-activated protein kinase cascade. Proceedings of the National Academy of Sciences93(19), 10078–10083 (1996)

1996
[12]

European journal of biochemistry267(6), 1583–1588 (2000) 18

Kholodenko, B.N.: Negative feedback and ultrasensitivity can bring about oscillations in the mitogen-activated protein kinase cascades. European journal of biochemistry267(6), 1583–1588 (2000) 18

2000
[13]

SIAM Journal on Numerical Analysis62(4), 1736–1758 (2024)

Scarabel, F., Vermiglio, R.: Equations with infinite delay: pseudospectral discretization for numerical stability and bifurcation in an abstract framework. SIAM Journal on Numerical Analysis62(4), 1736–1758 (2024)

2024
[14]

Rios, K.I.: Systematic characterisation of network structures underlying adaptive resistance and susceptibility to sequential combination treatment in cancer (2021) https://doi.org/10.26180/ 17209520.v1

2021
[15]

Journal of Mathematical Biology89(2), 26 (2024)

Jayathilaka, C., Araujo, R., Nguyen, L., Flegg, M.: Two wrongs do not make a right: the assump- tion that an inhibitor acts as an inverse activator. Journal of Mathematical Biology89(2), 26 (2024)

2024
[16]

Annals of eugenics7(4), 355–369 (1937)

Fisher, R.A.: The wave of advance of advantageous genes. Annals of eugenics7(4), 355–369 (1937)

1937
[17]

Applicable Mathematics to Non-physical Phenomena (1982)

Kolmogorov, A.: Study of the diffusion equation with growth of the quantity of matter and its applications to biological problems. Applicable Mathematics to Non-physical Phenomena (1982)

1982
[18]

Cell144(1), 119–131 (2011)

O’Shaughnessy, E.C., Palani, S., Collins, J.J., Sarkar, C.A.: Tunable signal processing in synthetic map kinase cascades. Cell144(1), 119–131 (2011)

2011
[19]

In: Partial Differential Equations and Related Topics: Ford Foundation Sponsored Program at Tulane University, January to May, 1974, pp

Aronson, D.G., Weinberger, H.F.: Nonlinear diffusion in population genetics, combustion, and nerve pulse propagation. In: Partial Differential Equations and Related Topics: Ford Foundation Sponsored Program at Tulane University, January to May, 1974, pp. 5–49. Springer, ??? (2006)

1974
[20]

Molecular cell9(5), 957–970 (2002) 19

Heinrich, R., Neel, B.G., Rapoport, T.A.: Mathematical models of protein kinase signal transduction. Molecular cell9(5), 957–970 (2002) 19

2002