Pith. sign in

REVIEW 5 minor 36 references

ssys turns supported ODE models into exact S-system or GMA form by symbolic lifting, with certificates that dynamics match on the constraint manifold.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · grok-4.5

2026-07-11 05:02 UTC pith:BSQ3UQNO

load-bearing objection Solid software paper that finally turns classical S-system/GMA recasting into a usable, validated open-source tool with real BioModels evidence.

arxiv 2607.05607 v1 pith:BSQ3UQNO submitted 2026-07-06 q-bio.QM

ssys: Exact algebraic recasting of ODE models into S-system or GMA form

William S. Hlavacek This is my paper
classification q-bio.QM
keywords S-systemGeneralized Mass Actionexact recastingsymbolic liftingODE modelsSBMLAntimonysystems biology
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

Classical Biochemical Systems Theory showed that many nonlinear ODEs can be rewritten exactly as S-systems or Generalized Mass Action systems by introducing auxiliary variables, but no general open-source tool did that rewriting for models already written in standard systems-biology formats. ssys fills that gap: it reads Antimony or SBML, lifts composite nonlinearities into power-law products, and emits a mathematically equivalent system whose right-hand sides are pure monomials or differences of two monomials. Exactness is local to trajectories that stay inside the admissible domain and whose auxiliary initial conditions obey the defining relations. A three-test suite (symbolic Jacobian check, pointwise numerical sampling, and trajectory comparison) certifies each recast. Benchmarks on curated textbook models and hundreds of BioModels entries show that the pipeline runs at repository scale and produces validated canonical forms for a large fraction of ODE models that meet its feature limits.

Core claim

ssys produces a mathematically equivalent S-system or GMA representation of supported ODE models: the original and transformed systems have identical dynamics on the invariant constraint manifold defined by the auxiliary-variable definitions, given consistent initial conditions, and this equivalence is certified by independent symbolic, numerical, and trajectory validators.

What carries the argument

Symbolic lifting of composite functions and rational denominators into auxiliary variables whose ODEs are derived by the chain rule, followed by pool-auxiliary construction that factors sums of monomials into products, yielding exact S-system or GMA form on the constraint manifold.

Load-bearing premise

Exactness holds only while trajectories stay inside the admissible real domain (no zero denominators, positive log and power-law bases, valid roots) and auxiliaries start on their defining relations; models that can hit zero or negative values must be preprocessed into the positive orthant by the user.

What would settle it

Take any supported model that ssys claims to recast successfully; if either the symbolic Jacobian identity fails, a random admissible point shows a nonzero residual between original and recast vector fields, or integrated trajectories diverge beyond numerical tolerance under matching initial conditions, the central exactness claim is false.

Watch this falsifier — get emailed when new claim-graph text bears on it.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

0 major / 5 minor

Summary. ssys is a Python package that performs exact algebraic recasting of supported ODE models (Antimony/SBML) into S-system or GMA form by introducing auxiliary variables via symbolic lifting. The recast is claimed to be mathematically equivalent to the original on the invariant constraint manifold defined by the auxiliary definitions (with consistent initial conditions), subject to stated domain restrictions. Correctness is certified by a three-test suite (symbolic Jacobian chain-rule verification, pointwise numerical sampling, and libRoadRunner trajectory comparison). The package includes a CLI, notebook generation, 712 pytest cases covering 117 curated classical models, and a BioModels-scale benchmark (978 candidates after filtering; 848 completed; 739 numerically validated).

Significance. If the claims hold, ssys fills a genuine tooling gap: classical BST recasting (Savageau & Voit 1987) has long been known in theory but lacked a general-purpose, open-source implementation that starts from standard SBML/Antimony models, produces validated S-system/GMA forms, and ships independent certificates of equivalence. Strengths that should be credited explicitly include the three independent validators, the large pytest suite with literature-derived targets, transparent BioModels filtering and counts, and public code/docs. These make exact power-law recasting practical for algebraic steady-state analysis, ground-truth generation for structure-learning methods, and reproducible systems-biology workflows.

minor comments (5)
  1. In the Summary, the distinction between relaxed S-system form (zero-valued RHS terms allowed) and canonical form (strictly positive coefficients) is clear, but a short forward pointer to how the simplified vs. canonical output modes implement this would help readers who jump from the equations to the CLI.
  2. BioModels benchmark paragraph: the cascade 1644 → 978 → 848 → 739 → 366 is informative; a one-sentence breakdown of the main failure modes among the non-validated remainder (unsupported features vs. timeouts vs. domain violations) would make the applicability claim even more transparent without changing the result.
  3. Statement of need notes that variables that can be zero or negative must be preprocessed into the positive orthant. A brief cross-reference to RECASTING.md (or a short example of a typical shift/offset) would reduce user friction for models outside the positive orthant.
  4. Software design: the pool-auxiliary construction for factoring sums of monomials is mentioned only briefly; one additional sentence or a pointer to the corresponding section of RECASTING.md would clarify how multi-term GMA equations are rewritten into two-term S-system form.
  5. Minor typography: 'V alidation infrastructure' and 'A vailability' appear to have stray spaces after the initial capital; fix for the camera-ready version.

Circularity Check

0 steps flagged

No significant circularity: classical algebraic recasting implemented as software and checked by independent validators, not by construction or self-fit.

full rationale

The manuscript is a software paper that implements the classical Savageau–Voit exact recasting construction (auxiliary-variable lifting of elementary nonlinearities into power-law form) for Antimony/SBML ODEs. Equivalence is claimed only on the invariant constraint manifold under consistent initial conditions and stated domain restrictions; it is certified by three independent checks (symbolic Jacobian/chain-rule verification, pointwise numerical sampling, and trajectory comparison via external libRoadRunner simulation). No parameters are fitted to data and then re-presented as predictions; BioModels and curated-model counts are descriptive success rates, not self-predictions. Citations to Savageau & Voit (1987) and related BST literature supply the known algebraic method rather than a load-bearing uniqueness theorem authored by the present author. There is therefore no self-definitional loop, no fitted-input-as-prediction, and no self-citation chain that forces the central claim. Score 0 is the honest finding.

Axiom & Free-Parameter Ledger

0 free parameters · 3 axioms · 0 invented entities

The paper rests on the classical Savageau–Voit recasting theorem and on standard symbolic-computation and SBML semantics; it introduces no free parameters that are fitted to force equivalence and no new physical entities. The main load-bearing assumptions are domain restrictions (positivity, non-vanishing denominators) that the user must enforce.

axioms (3)
  • standard math Systems built from sums, products and compositions of elementary functions can be recast into S-system form by introduction of auxiliary variables (Savageau & Voit 1987 theorem).
    Invoked throughout Summary and Statement of need as the theoretical foundation that ssys automates.
  • domain assumption Lifted variables used as power-law bases remain strictly positive, denominators avoid zero, and log/root arguments stay in their real domains along trajectories of interest.
    Stated explicitly as the condition for local exactness; models violating it are rejected or require user preprocessing.
  • domain assumption SBML/Antimony models that pass the package’s feature filters are faithfully represented by the SymPy expressions obtained via libSBML/Antimony.
    Parsing path described in Software design; unsupported semantics are rejected rather than silently accepted.

pith-pipeline@v1.1.0-grok45 · 11791 in / 2268 out tokens · 22886 ms · 2026-07-11T05:02:33.087136+00:00 · methodology

0 comments
read the original abstract

ssys is a Python package for exact algebraic recasting of supported ODE models into S-system or Generalized Mass Action form. It reads Antimony and SBML models, introduces auxiliary variables through symbolic lifting, and validates transformed systems using symbolic, numerical, and trajectory-based checks. The package provides command-line workflows, notebook generation, and benchmark evidence across curated models and BioModels examples, making classical power-law recasting practical for reproducible systems biology modeling.

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

36 extracted references · 24 canonical work pages

  1. [1]

    Anghel, M., Milano, F., & Papachristodoulou, A. (2013). Algorithmic construction of lyapunov functions for power system stability analysis. IEEE Transactions on Circuits and Systems I: Regular Papers, 60(9), 2533--2546. https://doi.org/10.1109/tcsi.2013.2246233

    doi:10.1109/tcsi.2013.2246233 2013

  2. [2]

    Anguelov, R., Borisov, M., Iliev, A., Kyurkchiev, N., & Markov, S. (2018). On the chemical meaning of some growth models possessing gompertzian‐type property. Mathematical Methods in the Applied Sciences, 41(18), 8365--8376. https://doi.org/10.1002/mma.4539

    work page doi:10.1002/mma.4539 2018

  3. [3]

    J., Keating, S

    Bornstein, B. J., Keating, S. M., Jouraku, A., & Hucka, M. (2008). libSBML : An API library for SBML . Bioinformatics, 24(6), 880--881. https://doi.org/10.1093/bioinformatics/btn051

    work page doi:10.1093/bioinformatics/btn051 2008

  4. [4]

    C., & Nemenman, I

    Daniels, B. C., & Nemenman, I. (2015a). Automated adaptive inference of phenomenological dynamical models. Nature Communications, 6, 8133. https://doi.org/10.1038/ncomms9133

    doi:10.1038/ncomms9133

  5. [5]

    C., & Nemenman, I

    Daniels, B. C., & Nemenman, I. (2015b). Efficient inference of parsimonious phenomenological models of cellular dynamics using S -systems and alternating regression. PLOS ONE, 10(3), e0119821. https://doi.org/10.1371/journal.pone.0119821

    work page doi:10.1371/journal.pone.0119821

  6. [6]

    R., Millman, K

    Harris, C. R., Millman, K. J., Walt, S. J. van der, Gommers, R., Virtanen, P., Cournapeau, D., Wieser, E., Taylor, J., Berg, S., Smith, N. J., Kern, R., Picus, M., Hoyer, S., Kerkwijk, M. H. van, Brett, M., Haldane, A., Fernández del Río, J., Wiebe, M., Peterson, P., Oliphant, T. E. (2020). Array programming with NumPy . Nature, 585(7825), 357--362. https...

    doi:10.1038/s41586-020-2649-2 2020

  7. [7]

    Hernández-Bermejo, B., Fairén, V., & Brenig, L. (1998). Algebraic recasting of nonlinear systems of ODEs into universal formats. Journal of Physics A: Mathematical and General, 31(10), 2415--2430. https://doi.org/10.1088/0305-4470/31/10/016

    work page doi:10.1088/0305-4470/31/10/016 1998

  8. [8]

    L., & Sauro, H

    Heydarabadipour, A., Smith, L., Hellerstein, J. L., & Sauro, H. M. (2026). Antimony 3: Extending human-readable model definitions for SBML level 3 core and packages. bioRxiv. https://doi.org/10.64898/2026.04.07.717118

    work page doi:10.64898/2026.04.07.717118 2026

  9. [9]

    Irvine, D. H. (1988). Efficient solution of nonlinear models expressed in S -system canonical form. Mathematical and Computer Modelling, 11, 123--128. https://doi.org/10.1016/0895-7177(88)90466-9

    work page doi:10.1016/0895-7177(88)90466-9 1988

  10. [10]

    M., Waltemath, D., König, M., Zhang, F., Dräger, A., Chaouiya, C., Bergmann, F

    Keating, S. M., Waltemath, D., König, M., Zhang, F., Dräger, A., Chaouiya, C., Bergmann, F. T., Finney, A., Gillespie, C. S., Helikar, T., Hoops, S., Malik-Sheriff, R. S., Moodie, S. L., Moraru, I. I., Myers, C. J., Naldi, A., Olivier, B. G., Sahle, S., Schaff, J. C., SBML Level 3 Community members. (2020). SBML level 3: An extensible format for the excha...

    work page doi:10.15252/msb.20199110 2020

  11. [11]

    Kluyver, T., Ragan-Kelley, B., Pérez, F., Granger, B., Bussonnier, M., Frederic, J., Kelley, K., Hamrick, J., Grout, J., Corlay, S., Ivanov, P., Avila, D., Abdalla, S., Willing, C., & Jupyter Development Team. (2016). Jupyter notebooks--a publishing format for reproducible computational workflows. In Positioning and power in academic publishing: Players, ...

    doi:10.3233/978-1-61499-649-1-87 2016

  12. [12]

    S., Glont, M., Nguyen, T

    Malik-Sheriff, R. S., Glont, M., Nguyen, T. V. N., Tiwari, K., Roberts, M. G., Xavier, A., Vu, M. T., Men, J., Maire, M., Kananathan, S., Fairbanks, E. L., Meyer, J. P., Arankalle, C., Varusai, T. M., Knight-Schrijver, V., Li, L., Dueñas-Roca, C., Dass, G., Keating, S. M., Hermjakob, H. (2020). BioModels ---15 years of sharing computational models in life...

    doi:10.1093/nar/gkz1055 2020

  13. [13]

    O., Gonzalez-Alcon, C., & Torres, N

    Marin-Sanguino, A., Voit, E. O., Gonzalez-Alcon, C., & Torres, N. V. (2007). Optimization of biotechnological systems through geometric programming. Theoretical Biology and Medical Modelling, 4(1), 38. https://doi.org/10.1186/1742-4682-4-38

    work page doi:10.1186/1742-4682-4-38 2007

  14. [14]

    P., Paprocki, M., Čertík, O., Kirpichev, S

    Meurer, A., Smith, C. P., Paprocki, M., Čertík, O., Kirpichev, S. B., Rocklin, M., Kumar, A., Ivanov, S., Moore, J. K., Singh, S., Rathnayake, T., Vig, S., Granger, B. E., Muller, R. P., Bonazzi, F., Gupta, H., Vats, S., Johansson, F., Pedregosa, F., Scopatz, A. (2017). SymPy : Symbolic computing in P ython. PeerJ Computer Science, 3, e103. https://doi.or...

    doi:10.7717/peerj-cs.103 2017

  15. [15]

    Papachristodoulou, A., & Prajna, S. (2005). Analysis of non-polynomial systems using the sum of squares decomposition. In Positive polynomials in control (Vol. 312, pp. 23--43). Springer Berlin Heidelberg. https://doi.org/10.1007/10997703_2

    work page doi:10.1007/10997703_2 2005

  16. [16]

    Pozo, C., Marín-Sanguino, A., Alves, R., Guillén-Gosálbez, G., Jiménez, L., & Sorribas, A. (2011). Steady-state global optimization of metabolic non-linear dynamic models through recasting into power-law canonical models. BMC Systems Biology, 5(1), 137. https://doi.org/10.1186/1752-0509-5-137

    work page doi:10.1186/1752-0509-5-137 2011

  17. [17]

    F., & Voit, E

    Rust, P. F., & Voit, E. O. (1990). Statistical densities, cumulatives, quantiles, and power obtained by S -system differential equations. Journal of the American Statistical Association, 85(410), 572--578. https://doi.org/10.1080/01621459.1990.10476237

    doi:10.1080/01621459.1990.10476237 1990

  18. [18]

    Savageau, M. A. (1971). Parameter sensitivity as a criterion for evaluating and comparing the performance of biochemical systems. Nature, 229(5286), 542--544. https://doi.org/10.1038/229542a0

    work page doi:10.1038/229542a0 1971

  19. [19]

    Savageau, M. A. (1976). Biochemical systems analysis: A study of function and design in molecular biology. Addison-Wesley. ISBN: 0-201-06739-0

    1976

  20. [20]

    Savageau, M. A. (1988). Introduction to S -systems and the underlying power-law formalism. Mathematical and Computer Modelling, 11, 546--551. https://doi.org/10.1016/0895-7177(88)90553-5

    work page doi:10.1016/0895-7177(88)90553-5 1988

  21. [21]

    Savageau, M. A. (1993). Finding multiple roots of nonlinear algebraic equations using S -system methodology. Applied Mathematics and Computation, 55(2--3), 187--199. https://doi.org/10.1016/0096-3003(93)90020-f

    work page doi:10.1016/0096-3003(93)90020-f 1993

  22. [22]

    A., & Voit, E

    Savageau, M. A., & Voit, E. O. (1987). Recasting nonlinear differential equations as S -systems: A canonical nonlinear form. Mathematical Biosciences, 87(1), 83--115. https://doi.org/10.1016/0025-5564(87)90035-6

    work page doi:10.1016/0025-5564(87)90035-6 1987

  23. [23]

    P., Bergmann, F

    Smith, L. P., Bergmann, F. T., Chandran, D., & Sauro, H. M. (2009). Antimony: A modular model definition language. Bioinformatics, 25(18), 2452--2454. https://doi.org/10.1093/bioinformatics/btp401

    work page doi:10.1093/bioinformatics/btp401 2009

  24. [24]

    T., Bouteiller, J.-M., Glazier, J

    Somogyi, E. T., Bouteiller, J.-M., Glazier, J. A., König, M., Medley, J. K., Swat, M. H., & Sauro, H. M. (2015). libRoadRunner : A high performance SBML simulation and analysis library. Bioinformatics, 31(20), 3315--3321. https://doi.org/10.1093/bioinformatics/btv363

    work page doi:10.1093/bioinformatics/btv363 2015

  25. [25]

    Vadhin, S., & Varner, J. D. (2026). BSTModelKit.jl : A Julia package for constructing, solving, and analyzing biochemical systems theory models. arXiv Preprint arXiv:2603.19115. https://arxiv.org/abs/2603.19115

    work page arXiv 2026

  26. [26]

    Vera, J., Sun, C., Oertel, Y., & Wolkenhauer, O. (2007). PLMaddon : A power-law module for the Matlab™ SBToolbox . Bioinformatics, 23(19), 2638--2640. https://doi.org/10.1093/bioinformatics/btm245

    work page doi:10.1093/bioinformatics/btm245 2007

  27. [27]

    Villaverde, A. F. (2019). Observability and structural identifiability of nonlinear biological systems. Complexity, 2019, 8497093. https://doi.org/10.1155/2019/8497093

    doi:10.1155/2019/8497093 2019

  28. [28]

    Voit, E. O. (1988a). New nonlinear methodologies for modelling molecular and cellular systems. Proceedings of the 1st IFAC Symposium on Modelling and Control in Biomedical Systems, Venice, Italy, 217--228. https://doi.org/10.1016/s1474-6670(17)57559-x

    work page doi:10.1016/s1474-6670(17)57559-x

  29. [29]

    Voit, E. O. (1988b). Recasting nonlinear models as S -systems. Mathematical and Computer Modelling, 11, 140--145. https://doi.org/10.1016/0895-7177(88)90469-4

    work page doi:10.1016/0895-7177(88)90469-4

  30. [30]

    Voit, E. O. (1990). S -system modelling of endemic infections. Computers & Mathematics with Applications, 20(4--6), 161--173. https://doi.org/10.1016/0898-1221(90)90325-e

    work page doi:10.1016/0898-1221(90)90325-e 1990

  31. [31]

    Voit, E. O. (1992). Symmetries of S -systems. Mathematical Biosciences, 109(1), 19--37. https://doi.org/10.1016/0025-5564(92)90050-7

    work page doi:10.1016/0025-5564(92)90050-7 1992

  32. [32]

    Voit, E. O. (1993). S -system modelling of complex systems with chaotic input. Environmetrics, 4(2), 153--186. https://doi.org/10.1002/env.3170040203

    work page doi:10.1002/env.3170040203 1993

  33. [33]

    Voit, E. O. (2000). Computational analysis of biochemical systems: A practical guide for biochemists and molecular biologists (1st ed.). Cambridge University Press. ISBN: 978-0521780872

    2000

  34. [34]

    O., & Kemp, M

    Voit, E. O., & Kemp, M. L. (2025). A first course in systems biology (3rd ed.). Garland Science. ISBN: 978-1040165478

    2025

  35. [35]

    Welsh, C., Xu, J., Smith, L., König, M., Choi, K., & Sauro, H. M. (2023). libRoadRunner 2.0: A high performance SBML simulation and analysis library. Bioinformatics, 39(1), btac770. https://doi.org/10.1093/bioinformatics/btac770

    work page doi:10.1093/bioinformatics/btac770 2023

  36. [36]

    Zhang, Z., Schuerhuber, R., Fickert, L., Friedl, K., Chen, G., & Zhang, Y. (2022). Domain of attraction's estimation for grid connected converters with phase-locked loop. IEEE Transactions on Power Systems, 37(2), 1351--1362. https://doi.org/10.1109/tpwrs.2021.3098960 CSLReferences document

    doi:10.1109/tpwrs.2021.3098960 2022