{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,6,2]],"date-time":"2026-06-02T22:47:49Z","timestamp":1780440469854,"version":"3.54.1"},"update-to":[{"DOI":"10.1371\/journal.pcbi.1011039","type":"new_version","label":"New version","source":"publisher","updated":{"date-parts":[[2023,4,25]],"date-time":"2023-04-25T00:00:00Z","timestamp":1682380800000}}],"reference-count":44,"publisher":"Public Library of Science (PLoS)","issue":"4","license":[{"start":{"date-parts":[[2023,4,13]],"date-time":"2023-04-13T00:00:00Z","timestamp":1681344000000},"content-version":"vor","delay-in-days":0,"URL":"http:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"DOI":"10.13039\/501100010446","name":"Institute for Basic Science","doi-asserted-by":"crossref","award":["IBS-R029-C3"],"award-info":[{"award-number":["IBS-R029-C3"]}],"id":[{"id":"10.13039\/501100010446","id-type":"DOI","asserted-by":"crossref"}]},{"name":"NSF","award":["DMS-2213390"],"award-info":[{"award-number":["DMS-2213390"]}]}],"content-domain":{"domain":["www.ploscompbiol.org"],"crossmark-restriction":false},"short-container-title":["PLoS Comput Biol"],"abstract":"<jats:p>The long-term behaviors of biochemical systems are often described by their steady states. Deriving these states directly for complex networks arising from real-world applications, however, is often challenging. Recent work has consequently focused on network-based approaches. Specifically, biochemical reaction networks are transformed into weakly reversible and deficiency zero generalized networks, which allows the derivation of their analytic steady states. Identifying this transformation, however, can be challenging for large and complex networks. In this paper, we address this difficulty by breaking the complex network into smaller independent subnetworks and then transforming the subnetworks to derive the analytic steady states of each subnetwork. We show that stitching these solutions together leads to the analytic steady states of the original network. To facilitate this process, we develop a user-friendly and publicly available package, COMPILES (COMPutIng anaLytic stEady States). With COMPILES, we can easily test the presence of bistability of a CRISPRi toggle switch model, which was previously investigated via tremendous number of numerical simulations and within a limited range of parameters. Furthermore, COMPILES can be used to identify absolute concentration robustness (ACR), the property of a system that maintains the concentration of particular species at a steady state regardless of any initial concentrations. Specifically, our approach completely identifies all the species with and without ACR in a complex insulin model. Our method provides an effective approach to analyzing and understanding complex biochemical systems.<\/jats:p>","DOI":"10.1371\/journal.pcbi.1011039","type":"journal-article","created":{"date-parts":[[2023,4,13]],"date-time":"2023-04-13T18:31:04Z","timestamp":1681410664000},"page":"e1011039","update-policy":"https:\/\/doi.org\/10.1371\/journal.pcbi.corrections_policy","source":"Crossref","is-referenced-by-count":9,"title":["A framework for deriving analytic steady states of biochemical reaction networks"],"prefix":"10.1371","volume":"19","author":[{"ORCID":"https:\/\/orcid.org\/0000-0002-5632-216X","authenticated-orcid":true,"given":"Bryan S.","family":"Hernandez","sequence":"first","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-7551-2911","authenticated-orcid":true,"given":"Patrick Vincent N.","family":"Lubenia","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Matthew D.","family":"Johnston","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-7842-2172","authenticated-orcid":true,"given":"Jae Kyoung","family":"Kim","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"340","published-online":{"date-parts":[[2023,4,13]]},"reference":[{"key":"pcbi.1011039.ref001","volume-title":"An introduction to systems biology: design principles of biological circuits","author":"U Alon","year":"2007"},{"key":"pcbi.1011039.ref002","volume-title":"Mathematical Modeling in Systems Biology: An Introduction","author":"BP Ingalls","year":"2013"},{"key":"pcbi.1011039.ref003","doi-asserted-by":"crossref","DOI":"10.1017\/CBO9781139854610","volume-title":"Systems Biology: Constraint-based Reconstruction and Analysis","author":"B Palsson","year":"2015"},{"key":"pcbi.1011039.ref004","doi-asserted-by":"crossref","first-page":"365","DOI":"10.1016\/j.copbio.2007.07.009","article-title":"Biochemical and statistical network models for systems biology","volume":"18","author":"N Price","year":"2007","journal-title":"Curr Opin Biotechnol"},{"key":"pcbi.1011039.ref005","volume-title":"Industrial Chemistry","author":"BK Sharma","year":"2000"},{"key":"pcbi.1011039.ref006","first-page":"35","article-title":"Studies Concerning Affinity","author":"CM Guldberg","year":"1864","journal-title":"C M Forhandlinger: Videnskabs-Selskabet i Chistiana"},{"key":"pcbi.1011039.ref007","unstructured":"Feinberg M. Lectures on chemical reaction networks; 1979. Available from: https:\/\/crnt.osu.edu\/LecturesOnReactionNetworks."},{"key":"pcbi.1011039.ref008","doi-asserted-by":"crossref","first-page":"81","DOI":"10.1007\/BF00251225","article-title":"General mass action kinetics","volume":"47","author":"F Horn","year":"1972","journal-title":"Arch Ration Mech Anal"},{"key":"pcbi.1011039.ref009","doi-asserted-by":"crossref","first-page":"172","DOI":"10.1007\/BF00255664","article-title":"Necessary and sufficient conditions for complex balancing in chemical kinetics","volume":"49","author":"F Horn","year":"1972","journal-title":"Arch Ration Mech Anal"},{"key":"pcbi.1011039.ref010","doi-asserted-by":"crossref","first-page":"187","DOI":"10.1007\/BF00255665","article-title":"Complex balancing in general kinetic systems","volume":"49","author":"M Feinberg","year":"1972","journal-title":"Arch Ration Mech Anal"},{"issue":"11","key":"pcbi.1011039.ref011","doi-asserted-by":"crossref","first-page":"1551","DOI":"10.1016\/j.jsc.2008.08.006","article-title":"Toric Dynamical Systems","volume":"44","author":"G Craciun","year":"2009","journal-title":"J Symbolic Comput"},{"key":"pcbi.1011039.ref012","doi-asserted-by":"crossref","first-page":"1926","DOI":"10.1137\/110847056","article-title":"Generalized mass action systems: Complex balancing equilibria and sign vectors of the stoichiometric and kinetic-order subspaces","volume":"72","author":"S M\u00fcller","year":"2012","journal-title":"SIAM J Appl Math"},{"key":"pcbi.1011039.ref013","doi-asserted-by":"crossref","first-page":"1081","DOI":"10.1007\/s11538-014-9947-5","article-title":"Translated chemical reaction networks","volume":"76","author":"MD Johnston","year":"2014","journal-title":"Bull Math Biol"},{"key":"pcbi.1011039.ref014","doi-asserted-by":"crossref","first-page":"1389","DOI":"10.1126\/science.1183372","article-title":"Structural sources of robustness in biochemical reaction networks","volume":"327","author":"G Shinar","year":"2010","journal-title":"Science"},{"issue":"1","key":"pcbi.1011039.ref015","doi-asserted-by":"crossref","first-page":"39","DOI":"10.1016\/j.mbs.2011.02.012","article-title":"Design principles for robust biochemical reaction networks: What works, what cannot work, and what might almost work","volume":"231","author":"G Shinar","year":"2011","journal-title":"Math Biosci"},{"key":"pcbi.1011039.ref016","doi-asserted-by":"crossref","first-page":"2306","DOI":"10.1007\/s11538-018-0458-7","article-title":"Network translation and steady-state properties of chemical reaction systems","volume":"80","author":"E Tonello","year":"2018","journal-title":"Bull Math Biol"},{"key":"pcbi.1011039.ref017","doi-asserted-by":"crossref","first-page":"1027","DOI":"10.1007\/s11538-011-9685-x","article-title":"Chemical Reaction Systems with Toric Steady States","volume":"74","author":"M P\u00e9rez Mill\u00e1n","year":"2012","journal-title":"Bull Math Biol"},{"issue":"10","key":"pcbi.1011039.ref018","doi-asserted-by":"crossref","first-page":"e1005751","DOI":"10.1371\/journal.pcbi.1005751","article-title":"Identifying parameter regions for multistationarity","volume":"13","author":"C Conradi","year":"2017","journal-title":"PLOS Comp Bio"},{"key":"pcbi.1011039.ref019","doi-asserted-by":"crossref","first-page":"1065","DOI":"10.1007\/s11538-015-0077-5","article-title":"A computational approach to steady state correspondence of regular and generalized mass action systems","volume":"77","author":"MD Johnston","year":"2015","journal-title":"Bull Math Biol"},{"key":"pcbi.1011039.ref020","doi-asserted-by":"crossref","first-page":"1143","DOI":"10.1007\/s11538-018-00562-0","article-title":"A deficiency-based approach to parametrizing positive equilibria of biochemical reaction systems","volume":"81","author":"MD Johnston","year":"2019","journal-title":"Bull Math Biol"},{"key":"pcbi.1011039.ref021","doi-asserted-by":"crossref","first-page":"1613","DOI":"10.1007\/s11538-019-00579-z","article-title":"Computing weakly reversible deficiency zero network translations using elementary flux modes","volume":"81","author":"MD Johnston","year":"2019","journal-title":"Bull Math Biol"},{"key":"pcbi.1011039.ref022","first-page":"1","article-title":"Prevalence of deficiency zero reaction networks in an Erd\u00f6s\u2013R\u00e9nyi framework","author":"DF Anderson","year":"2022","journal-title":"J Appl Probab"},{"key":"pcbi.1011039.ref023","doi-asserted-by":"crossref","DOI":"10.1137\/22M150469X","article-title":"Computational translation framework identifies biochemical reaction networks with special topologies and their long-term dynamics","author":"H Hong","year":"2023","journal-title":"SIAM J Appl Math"},{"key":"pcbi.1011039.ref024","doi-asserted-by":"crossref","first-page":"2229","DOI":"10.1016\/0009-2509(87)80099-4","article-title":"Chemical reaction network structure and the stability of complex isothermal reactors I: The deficiency zero and deficiency one theorems","volume":"42","author":"M Feinberg","year":"1987","journal-title":"Chem Eng Sci"},{"key":"pcbi.1011039.ref025","doi-asserted-by":"crossref","DOI":"10.1007\/978-3-030-03858-8","volume-title":"Foundations of Chemical Reaction Network Theory","author":"M Feinberg","year":"2019"},{"key":"pcbi.1011039.ref026","doi-asserted-by":"crossref","first-page":"2746","DOI":"10.1038\/s41467-020-16574-1","article-title":"Multistable and dynamic CRISPRi-based synthetic circuits","volume":"11","author":"J Santos-Moreno","year":"2020","journal-title":"Nat Commun"},{"key":"pcbi.1011039.ref027","doi-asserted-by":"crossref","first-page":"E1084","DOI":"10.1152\/ajpendo.00571.2001","article-title":"A mathematical model of metabolic insulin signaling pathways","volume":"283","author":"AR Sedaghat","year":"2002","journal-title":"Am J Physiol Endocrinol Metab"},{"key":"pcbi.1011039.ref028","doi-asserted-by":"crossref","DOI":"10.1017\/CBO9780511609589","volume-title":"Enumerative Combinatorics","author":"R Stanley","year":"1999"},{"key":"pcbi.1011039.ref029","doi-asserted-by":"crossref","first-page":"129","DOI":"10.1007\/s11538-022-01087-3","article-title":"Reaction Network Analysis of Metabolic Insulin Signaling","volume":"84","author":"PVN Lubenia","year":"2022","journal-title":"Bull Math Biol"},{"key":"pcbi.1011039.ref030","doi-asserted-by":"crossref","first-page":"623","DOI":"10.1007\/s10013-021-00524-5","article-title":"Absolute Concentration Robustness in Networks with Low-Dimensional Stoichiometric Subspace","volume":"50","author":"N Meshkat","year":"2022","journal-title":"Vietnam J Math"},{"key":"pcbi.1011039.ref031","doi-asserted-by":"crossref","first-page":"76","DOI":"10.1109\/TAC.2015.2427691","article-title":"New approach to the stability of chemical reaction networks: piecewise linear in rates Lyapunov functions","volume":"61","author":"M Ali Al-Radhawi","year":"2016","journal-title":"IEEE Trans Automat Contr"},{"key":"pcbi.1011039.ref032","doi-asserted-by":"crossref","first-page":"707","DOI":"10.3934\/dcds.2019029","article-title":"Planar S-systems: Global stability and the center problem","volume":"39","author":"B Boros","year":"2019","journal-title":"Discrete Contin Dyn Syst Ser A"},{"key":"pcbi.1011039.ref033","doi-asserted-by":"crossref","first-page":"899","DOI":"10.1137\/120898486","article-title":"Local and global stability of equilibria for a class of chemical reaction networks","volume":"12","author":"P Donnell","year":"2013","journal-title":"SIAM J Appl Dyn Syst"},{"key":"pcbi.1011039.ref034","doi-asserted-by":"crossref","first-page":"1253","DOI":"10.1007\/s00285-019-01394-9","article-title":"Discrepancies between extinction events and boundary equilibria in reaction networks","volume":"79","author":"D Anderson","year":"2019","journal-title":"J Math Biol"},{"key":"pcbi.1011039.ref035","doi-asserted-by":"crossref","first-page":"140","DOI":"10.1016\/j.mbs.2008.09.001","article-title":"Homotopy methods for counting reaction network equilibria","volume":"216","author":"G Craciun","year":"2008","journal-title":"Math Biosci"},{"key":"pcbi.1011039.ref036","doi-asserted-by":"crossref","first-page":"191","DOI":"10.1016\/j.amc.2017.12.012","article-title":"Inheritance of oscillation in chemical reaction networks","volume":"325","author":"M Banaji","year":"2018","journal-title":"Appl Math Comput"},{"key":"pcbi.1011039.ref037","doi-asserted-by":"crossref","first-page":"442","DOI":"10.3934\/mbe.2020024","article-title":"Complex-balanced equilibria of generalized mass-action systems: necessary conditions for linear stability","volume":"17","author":"B Boros","year":"2020","journal-title":"Math Biosci Eng"},{"key":"pcbi.1011039.ref038","unstructured":"M\u00fcller S, Regensburger G. Sufficient conditions for linear stability of complex-balanced equilibria in generalized mass-action systems. arXiv 2022."},{"key":"pcbi.1011039.ref039","doi-asserted-by":"crossref","first-page":"1947","DOI":"10.1007\/s11538-010-9517-4","article-title":"Product-form stationary distributions for deficiency zero chemical reaction networks","volume":"72","author":"DF Anderson","year":"2010","journal-title":"Bull Math Biol"},{"issue":"6","key":"pcbi.1011039.ref040","doi-asserted-by":"crossref","first-page":"e1005571","DOI":"10.1371\/journal.pcbi.1005571","article-title":"Reduction of multiscale stochastic biochemical reaction networks using exact moment derivation","volume":"13","author":"JK Kim","journal-title":"PLoS Comput Biol"},{"key":"pcbi.1011039.ref041","doi-asserted-by":"crossref","first-page":"620","DOI":"10.1038\/s42003-021-02117-x","article-title":"Derivation of stationary distributions of biochemical reaction networks via structure transformation","volume":"4","author":"H Hong","year":"2021","journal-title":"Commun Biol"},{"issue":"7","key":"pcbi.1011039.ref042","doi-asserted-by":"crossref","first-page":"76","DOI":"10.1007\/s11538-021-00906-3","article-title":"Independent decompositions of chemical reaction networks","volume":"83","author":"BS Hernandez","year":"2021","journal-title":"Bull Math Biol"},{"issue":"2","key":"pcbi.1011039.ref043","doi-asserted-by":"crossref","first-page":"367","DOI":"10.46793\/match.87-2.367H","article-title":"Independent, incidence independent and weakly reversible decompositions of chemical reaction networks","volume":"87","author":"BS Hernandez","year":"2022","journal-title":"MATCH Commun Math Comput Chem"},{"key":"pcbi.1011039.ref044","first-page":"302","volume-title":"Computer Algebra in Scientific Computing. CASC 2014. Lecture Notes in Computer Science","author":"S M\u00fcller","year":"2014"}],"updated-by":[{"DOI":"10.1371\/journal.pcbi.1011039","type":"new_version","label":"New version","source":"publisher","updated":{"date-parts":[[2023,4,25]],"date-time":"2023-04-25T00:00:00Z","timestamp":1682380800000}}],"container-title":["PLOS Computational Biology"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/dx.plos.org\/10.1371\/journal.pcbi.1011039","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2023,12,10]],"date-time":"2023-12-10T17:33:28Z","timestamp":1702229608000},"score":1,"resource":{"primary":{"URL":"https:\/\/dx.plos.org\/10.1371\/journal.pcbi.1011039"}},"subtitle":[],"editor":[{"given":"Mark","family":"Alber","sequence":"first","affiliation":[],"role":[{"vocabulary":"crossref","role":"editor"}]}],"short-title":[],"issued":{"date-parts":[[2023,4,13]]},"references-count":44,"journal-issue":{"issue":"4","published-online":{"date-parts":[[2023,4,13]]}},"URL":"https:\/\/doi.org\/10.1371\/journal.pcbi.1011039","relation":{"new_version":[{"id-type":"doi","id":"10.1371\/journal.pcbi.1011039","asserted-by":"object"}]},"ISSN":["1553-7358"],"issn-type":[{"value":"1553-7358","type":"electronic"}],"subject":[],"published":{"date-parts":[[2023,4,13]]}}}