In silico modeling of the FR was conducted using previously published methods 47 which potentially allows for the identification of potential autocatalytic pathways based on the open-source Python MØD platform , a graph theory-based chemical reaction modeling software package that can generate CRNRs using user-defined reaction rules (see SI File 1). The FR network explored here was computed using MØD run on a Hewlett-Packard Z820 computer with 128 GB RAM and 16 cores in two Intel® Xeon® E5-2670 2.60 GHz processors, running Ubuntu 20.04 in a virtual machine hosted on a Windows 10 system. MØD uses state-of-the-art graph canonicalization algorithms which are efficient at dealing with large collections of molecular graphs .