Chemical reaction networks (CRNs) are essential for modeling and analyzing complex systems across fields, from biochemistry to economics. Autocatalytic reaction network -- networks where certain species catalyze their own production -- are particularly significant for understanding self-replication dynamics in biological systems and serve as foundational elements in formalizing the concept of a circular economy. In a previous study, we developed a mixed-integer linear optimization-based procedure to enumerate all minimal autocatalytic subnetworks within a network. In this work, we define the maximum growth factor (MGF) of an autocatalytic subnetwork, develop mathematical optimization approaches to compute this metric, and explore its implications in the field of economics and dynamical systems. We develop exact approaches to determine the MGF of any subnetwork based on an iterative procedure with guaranteed convergence, which allows for identifying autocatalytic subnetworks with the highest MGF. We report the results of computational experiments on synthetic CRNs and two well-known datasets, namely the Formose and E. coli reaction networks, identifying their autocatalytic subnetworks and exploring their scientific ramifications. Using advanced optimization techniques and interdisciplinary applications, our framework adds an essential resource to analyze complex systems modeled as reaction networks.
翻译:暂无翻译