Finally, in Step IV of the procedure, the time-dependent dynamics of the network measures is analyzed. This is done, since reactions take place all the time but at specific moments during the reaction trajectory certain ones are more prevalent (depending on the change in the concentration of the species involved in the reaction). Therefore, we simulated the CRN using the model described by ODEs. Fig. shows the characteristic behaviour of the oscillator in time, wherein the concentration of key intermediates Tr, [Tr], and Int-I, [Int-I], oscillate with the same frequency. The simulated time series are used as input to represent the temporal evolution of the species to determine the addition or removal of edges in G2. Briefly, we considered that an outgoing edge can only exist when the associated node is present, and used a threshold to determine the presence of the nodes at each time interval (for details, see Supplementary Information S3.2). To ensure that the sustained oscillations had established, our analysis starts from 50 h onwards.