Figure 9: Simulation of flux distribution under CO2 elevation. Simulating flux distribution at any key branch-points within the network infers a formal metabolic route response in the different starch formation under CO2 elevation (grey arrows): fluxes through (a) aldolase (V5) and triose-phosphate transporter (V31-GAP and V31a-DHAP); (b) starch phosphorylase (V19), beta-amylase (V20), disproportional enzyme (V21 and V23), and alpha-amylase (V22); (c) phosphoglucomutase (V33), phosphoglucoisomerase (V38) and sucrose phosphate synthase (V40); (d) sucrose synthase (V42) and invertase (V45); (e) tuber phosphoglucomutase (V44) and cytosolic ADPG pyrophosphorylase (V52); (f) G6P transporter (V55a-ADPG) and plastidic ADPG pyrophosphorylase (V57); (g) granule-bound starch synthase (V58) and starch synthase (V59). CO2 concentrations (x-axis) and fluxes (y-axis) are simulated in mM and mM/s, respectively. Clear symbols are results plotted against the left y-axis, while filled symbols are those plotted against the right y-axis. Grey arrows superimposed on each cascade represent the preferential metabolic route of plants, suggested by flux distribution (summarised in Figure 10).