Yanfei Zhang, Theresa Kouril, Jacky L. Snoep, Bettina Siebers, Matteo Barberis, Hans V. Westerhoff

Mathematical models are key to systems biology where they typically describe the topology and dynamics of biological networks, listing biochemical entities and their relationships with one another. Some (hyper)thermophilic Archaea contain an enzyme, called non-phosphorylating glyceraldehyde-3-phosphate dehydrogenase (GAPN), which catalyzes the direct oxidation of glyceraldehyde-3-phosphate to 3-phosphoglycerate omitting adenosine 5'-triphosphate (ATP) formation by substrate-level-phosphorylation via phosphoglycerate kinase. In this study we formulate three hypotheses that could explain functionally why GAPN exists in these Archaea, and then construct and use mathematical models to test these three hypotheses. We used kinetic parameters of enzymes of Sulfolobus solfataricus (S. solfataricus) which is a thermo-acidophilic archaeon that grows optimally between 60 and 90 °C and between pH 2 and 4. For comparison, we used a model of Saccharomyces cerevisiae (S. cerevisiae), an organism that can live at moderate temperatures. We find that both the first hypothesis, i.e., that the glyceraldehyde-3-phosphate dehydrogenase (GAPDH) plus phosphoglycerate kinase (PGK) route (the alternative to GAPN) is thermodynamically too much uphill and the third hypothesis, i.e., that GAPDH plus PGK are required to carry the flux in the gluconeogenic direction, are correct. The second hypothesis, i.e., that the GAPDH plus PGK route delivers less than the 1 ATP per pyruvate that is delivered by the GAPN route, is only correct when GAPDH reaction has a high rate and 1,3-bis-phosphoglycerate (BPG) spontaneously degrades to 3PG at a high rate.