The first step in the procedure encodes the enzymatic oscillator as a species-species network (Supplementary Information, S2). In here, each edge represents a unique chemical conversion from one species to another represented by nodes. Our procedure transforms the elementary reactions into a network. The species are labeled as reactants, intermediates and products resulting in a total of 10 nodes; three reactants (Tg, Ap, Pro-I), six intermediates (Tr, Int-I, I, Tr•Tg, Tr•Pro-I, Ap•Int-I), and one product (P) that summarizes the species that do not contribute to the overall dynamics (the inhibited trypsin complex and hydrolysed inhibitors). Subsequently, we assigned a directed edge for each pair of species in an elementary reaction, with the direction of the edge indicating whether the species was consumed or produced. We note that multiple chemical conversions could take place between the same pair of species. To account for this multiplicity, we defined weights based on the number of edges assigned to each pair of species. Supplementary Information, Table show the outputs of the algorithm applied on the enzymatic oscillator, where the final column (a so-called edge list) defines the species-species network: A graph with 10 nodes and 19 weighted edges (G1).