- Research article
- Open Access
In silico metabolic network analysis of Arabidopsis leaves
BMC Systems Biology volume 10, Article number: 102 (2016)
During the last decades, we face an increasing interest in superior plants to supply growing demands for human and animal nutrition and for the developing bio-based economy. Presently, our limited understanding of their metabolism and its regulation hampers the targeted development of desired plant phenotypes. In this regard, systems biology, in particular the integration of metabolic and regulatory networks, is promising to broaden our knowledge and to further explore the biotechnological potential of plants.
The thale cress Arabidopsis thaliana provides an ideal model to understand plant primary metabolism. To obtain insight into its functional properties, we constructed a large-scale metabolic network of the leaf of A. thaliana. It represented 511 reactions with spatial separation into compartments. Systematic analysis of this network, utilizing elementary flux modes, investigates metabolic capabilities of the plant and predicts relevant properties on the systems level: optimum pathway use for maximum growth and flux re-arrangement in response to environmental perturbation. Our computational model indicates that the A. thaliana leaf operates near its theoretical optimum flux state in the light, however, only in a narrow range of photon usage. The simulations further demonstrate that the natural day-night shift requires substantial re-arrangement of pathway flux between compartments: 89 reactions, involving redox and energy metabolism, substantially change the extent of flux, whereas 19 reactions even invert flux direction. The optimum set of anabolic pathways differs between day and night and is partly shifted between compartments. The integration with experimental transcriptome data pinpoints selected transcriptional changes that mediate the diurnal adaptation of the plant and superimpose the flux response.
The successful application of predictive modelling in Arabidopsis thaliana can bring systems-biological interpretation of plant systems forward. Using the gained knowledge, metabolic engineering strategies to engage plants as biotechnological factories can be developed.
Increasing world population, shortage of arable land and the resulting growing demand for food, feed and raw materials are major drivers to create plant lines with increased performance, e.g. better resistance to disease and drought . In addition, plants play a significant role for the developing bio-economy  and emerge as platforms for sustainable production of therapeutics, renewable chemicals and biofuels, purely from sunlight and carbon dioxide , which adds impetus to the growing interest in enhanced crops [3, 4]. Admittedly, plant metabolic engineering and plant systems metabolic engineering is hampered by our still limited understanding of the underlying metabolism and its regulation [5, 6], involving only little systems-level understanding of the effect of genetic modifications [7, 8]. It is therefore easy to understand that the interest in methods to drive plant metabolic engineering is high [6, 7, 9]. One promising approach to support plant metabolic engineers, obviously, is the use of in silico metabolic modelling. Firstly, in silico metabolic modelling has proven impressively successful in guiding systems metabolic engineering of other biological systems, including bacteria and fungi [10–14]. Secondly, the over 50, fully sequenced plant genomes, mainly crops,  provide a valuable source of information to create plant metabolic models [16–18]. And thirdly, a powerful collection of modelling approaches, currently available to model and simulate stoichiometric metabolic networks, can be adapted to plant metabolic networks in a straightforward manner . Among the available modelling approaches are in silico based analyses, such as elementary flux mode analysis [20, 21] and extreme pathway analysis [22, 23], as well as analyses, which rely on experimental data to deliver necessary constraints, such as 13C-metabolic flux analysis (13C-MFA) [24–26], high-throughput isotope based metabolic screening , flux balance analysis  and metabolic control analysis . Among the approaches, in silico simulation appears particularly interesting due to its high speed, given e.g. the relatively long time of experiments with growing plants.
Here, we conduct detailed in silico analysis of plant central carbon metabolism. We use A. thaliana as a widely applied model plant to study plant biology and biotechnology, reflected also by the fact that it was the first plant sequenced . In short, a metabolic network model of the A. thaliana leaf was constructed, which represented 511 reactions and 1567 metabolic genes. Its reactions and pathways were localized in different subcellular compartments: cytosol, plastid, mitochondrion and peroxisome. Through simulation, relevant physiological scenarios of plant metabolism were studied. This involved the detailed investigation of the diurnal metabolism through comparison of metabolic capabilities in the light and in the dark. In addition, the modelling results were integrated with previous experimental fluxomic and transcriptomic data of A. thaliana leaves in order to explore the systems-wide regulation of the metabolism. This provides an enhanced understanding of plant metabolic pathways and their contribution to growth and product formation.
Metabolic network construction
The metabolic network used in this study comprises the central carbon metabolism of A. thaliana leaves. For re-construction, the initial draft of the network built on recent genome-scale models of the organism [30–32]. The derived pathway bibliography was then carefully checked and updated with recent findings, collected in metabolic pathway databases: Kyoto Encyclopedia of Genes and Genomes  (http://www.genome.jp/kegg), MetaCrop  (http://metacrop.ipk-gatersleben.de/) and AraCyc  (http://pmn.plantcyc.org). This provided state-of-art gross information on the genomic pathway repertoire. Where needed, the network was updated with experimental data and primary literature as described in detail below. Individual additions and specifications considered enzyme localization, cofactor usage, and inter-compartmental transport. Subsequently, the model was condensed, which, however, did not reduce its information content (Additional file 1: Information S1). Furthermore, the applied toolbox for the enumeration of elementary flux modes required in silico external exchange reactions for unbalanced metabolites (biomass, ATP for cellular maintenance, inorganic phosphate, photons and starch), which were added. Both models, the detailed version and the condensed version, used for the simulations, are provided as SBML files in the Additional files 2 and 3.
Computation of elementary flux modes
Elementary flux modes (EFMs) were calculated with efmtool, based on the null space approach and recursive enumeration with bit pattern trees . The EFM matrix, computed by the algorithm, comprises information on all thermodynamically and stoichiometrically possible pathways in the cell, which reduce metabolism into all feasible, unique, non-decomposable biochemical pathways . All simulations were conducted on a quad core personal computer. Normalization of the EFM matrix and subsequent data interpretation was conducted as described previously . In short, relative fluxes were normalized to their respective substrate uptake and the theoretical biomass production of each elementary flux mode was determined. The respective flux modes were expressed in (C-mol) (C-mol)−1. The relative contribution of a particular pathway to anabolic precursor formation was obtained by dividing the underlying pathway flux by the sum of all fluxes, which formed this anabolic precursor. Significant changes between two conditions were identified on the basis of a two-sample t-test (95 % significance level, p-value <0.05) and an absolute log2-value > 0.5.
Transcriptome data processing
Experimental transcriptome data were taken from a recent study on the shift of gene expression between day and night in A. thaliana rosettes . In this previous work, the amplitude limit of gene expression (log2 value) during the diurnal cycle was quantified by ATH1 arrays, and changes in absolute gene expression level were identified with a cut off value of 0.8. From the published data set, we extracted the genes encoding proteins of the central carbon metabolism. The obtained raw data were then inspected to identify the genes, which exhibited a diurnal expression change, i.e. revealed an unambiguous increase or decrease in expression during illumination and the opposite change during the dark period.
Fluxome data processing
The experimentally measured metabolic fluxes of an illuminated A. thaliana rosette  were converted into (C-mol) (C-mol)−1 to enable a straightforward comparison with the respective flux modes, obtained in this work, also given in (C-mol) (C-mol)−1.
Quantum yield data processing
The quantum yield of photosynthesis was derived from previous measurement of the rate of photosynthesis of Arabidopsis and of other C3 plants under ambient atmospheric conditions at different light intensity [39–43]. Due to the fact, that not all harvested quanta are converted into chemical energy, as some are lost through absorption by pigments, unable to contribute their excitation energy to photosynthesis, the experimental quantum yield values were corrected, assuming that 47 % of photons are outside the photosynthetically active range . This provided a direct correlation between assimilated carbon dioxide and properly absorbed photons.
Metabolic network topology
The created metabolic network of A. thaliana’s core carbon pathways reflected 511 metabolic conversions in total (Fig. 1). In total, 511 metabolic reactions are accounted for: 82 cytosolic, 348 plastidic, 27 mitochondrial and 6 peroxisomal reactions (see Additional file 2). Additionally, 41 reactions describe intracellular transport and 5 reactions mediate the uptake of nutrients. The derived network is descriptive for about 1567 metabolic genes, all known to exist in A. thaliana (http://pmn.plantcyc.org). The chosen topology is specific for green leaf tissue. The network contained the Calvin-Benson-Bassham (CBB) cycle for CO2 assimilation and photosynthesis, the gluconeogenesis, the glycolytic Embden-Meyerhof-Parnas (EMP) pathway, the oxidative and the non-oxidative pentose phosphate (PP) pathway, the tricarboxylic acid (TCA) cycle, the reactions of photorespiration, of starch biosynthesis and degradation, of energy metabolism, as well as the anabolic pathways for biomass synthesis (Figs. 1 and 2). The latter considered compartment-specific supply of the individual precursors (Additional file 1: Table S1 and S2) [17, 30]. Internal starch, degraded in green tissue in the dark, and atmospheric CO2, assimilated in the light, were included as natural carbon sources. The network model specifically included metabolic conversions that contribute to leaf growth, whereas growth-unrelated metabolic processes, e.g. chemical defense, were omitted. For the simulations, the model was condensed. Single reactions in linear pathways and stoichiometrically identical parallel pathways were lumped together. The condensation was crucial to enable the computationally demanding simulation of the entire set of elementary flux modes. In this way, the total equation system to describe the plant metabolism could be represented in a compact form by only 129 reactions and 70 metabolites (Fig. 1, Additional file 1: Table S3 and S4, Additional file 3). The condensed model remained fully descriptive of the entire set of metabolic conversions, as demonstrated by additional simulations (Additional file 1: information S1). Among the reactions, 34 reactions belonged to inter-compartmental and to extracellular transport. Based on the information provided in the above-mentioned databases and in the literature, 63 reactions were constrained as irreversible (Additional file 1: Table S5). The condensed model topology allowed the computationally demanding calculation of all feasible metabolic routes in the leaf, through enumeration of elementary flux modes. The model was validated and fine structured with regard to the presence, location and thermodynamic properties of individual reactions, based on experimental evidence, as described below.
The network reflected the four major compartments of plant leaves that contribute to biochemical conversions: cytosol, peroxisome, mitochondrion and plastid (http://metacrop.ipk.gatersleben.de) (http://pmn.plantcyc.org) [31, 32, 45–48]. The cytosol comprised the reactions of the EMP pathway , the oxidative part of the PP pathway  and the reactions of starch degradation from maltose . The plastid contained the photosynthetic CBB cycle, a second copy of the EMP pathway , the oxidative and the non-oxidative PP pathway  and the starch metabolism . The TCA cycle was assigned to the mitochondrion . The photo-respiratory system, known to be a multi-compartment process, was distributed accordingly across plastid, peroxisome and mitochondrion . Additionally, each compartment contained malate dehydrogenase (http://metacrop.ipk.gatersleben.de) . Pyruvate kinase was considered as cytosolic and as plastidic reaction , whereas the pyruvate dehydrogenase complex was assigned to the mitochondrion and to the plastid . The supply of cytosolic acetyl-CoA was attributed to ATP-citrate lyase in the cytoplasm, which uses citrate as a substrate . Malic enzyme, specific for photosynthetic tissue, was integrated into the plastid .
Inter-compartmental and external transport
The separation of the metabolic routes in distinct organelles requires the translocation of specific compounds across cellular membranes. Based on experimental evidence, unidirectional or bidirectional transport between cytosol and mitochondrion was assumed for pyruvate, malate, inorganic phosphate, glycine and serine, respectively, whereas antiporters were considered for malate/oxaloacetate, citrate/oxaloacetate and ATP/ADP [34, 52–55]. Transport reactions between the cytosol and the plastid were implemented for 3-phosphoglycerate, glycerate, glycolate, malate/oxaloacetate, pyruvate, phosphoenolpyruvate, xylulose 5-phosphate, glucose 6-phosphate, dihydroxyacetone phosphate, maltose, ATP/ADP and inorganic phosphate [34, 45, 56–59]. Hereby, the translocation of phosphorylated carbohydrates across the plastid membranes was linked to the simultaneous antiport of inorganic phosphate . In addition, active peroxisomal membrane transfer of malate/oxaloacetate, glycerate, glycolate, glycine and serine was considered [49, 50]. So far, a transporter for acetyl-CoA has not been discovered and was therefore not incorporated . CO2 was assumed to freely diffuse within the cell , photons were capable of penetrating both the extracellular and plastidic membranes  and inorganic phosphate was available from the vacuole .
Redox, energy and phosphate metabolism were compartmentalized across the different organelles. This included the confinement of the photosynthetic light reactions to the plastid and of the oxidative phosphorylation system to the mitochondrion. A vacuolar storage pool for inorganic phosphate was considered . Inorganic phosphate could be transported by specific carriers between the cytosol and both the plastid and the mitochondrion. ATP/ADP antiporters were incorporated into the mitochondrial and the plastidic membrane, respectively. In addition, the exchange of reducing equivalents between cytosol and plastid, mitochondrion and peroxisome was attributed to malate dehydrogenase, which was coupled to malate and oxaloacetate channeling across the organelle barriers. To account for widely abundant isoenzymes, capable of utilizing either NADPH or NADH or both molecules as cofactor, an oxidoreductase for interconversion of NADPH and NADH was included in the cytosol . This compartmented consideration is more realistic. Additional simulations of a simplified network without this strongly compartment-specific energy, redox and phosphate metabolism showed that the maximum biomass formation was not affected by the compartmentalized energy, redox and phosphate acquisition (data not shown). This demonstrates that the network allows for full equilibration of energy, redox and phosphate among the compartments. The following considerations were additionally included to implement energy efficiency. The ratio of ATP: NADPH, produced in a photosynthetic cell, depends on the generally accepted plasticity of the photosynthetic light reactions for energy production (Fig. 2) [65–69]. As it is still unresolved, how this mechanism functions exactly, an average ATP to NADPH ratio of 1.5 was usually chosen, which accounts for cyclic electron flow around photosystem I of two photons and of non-cyclic electron flow caused by eight photons . In a set of additional simulations, the photosynthetic plasticity was investigated by varying the ratio between cyclic and non-cyclic electron flow, i.e. the ATP to NADPH ratio (see below). In order to provide sufficient ATP for maintenance, a surplus of ATP was set as constraint for the modelling . This excluded unrealistic scenarios, which would have produced less ATP than needed for growth, but allowed for scenarios, which produced an apparent excess of ATP to serve for maintenance purposes.
The biochemical composition of A. thaliana leaves was collected from previous work (Additional file 1: Table S1). Nearly half of cellular carbon is stored in the cell wall [71–73], whereas one third is contained in proteins . Based on experimental data, the remaining carbon was distributed among lipids [75–77], carbohydrates [71, 74, 78], porphyrins [79, 80] and other biomass components [74, 81–83]. The anabolic pathways form the most important carbon sink during growth. As most of these peripheral biosynthetic pathways, are linear and a metabolic steady-state is assumed, they can easily be summarized into a single, lumped biomass equation, starting from 12 central metabolic precursor metabolites (Additional file 1: Table S2), without noteworthy degeneration of information content (Additional file 1: Information S1). The intracellular localization of the individual biosynthetic enzymes determined from which organelle a particular precursor originated. For instance, aromatic amino acids are synthesized from phosphoenolpyruvate and erythrose 4-phosphate in the plastid , whereas cellulose originates from hexose 6-phosphate in the cytosol .
Plant metabolism involves highly efficient carbon assimilation and conversion
Plants are subjected to changing environmental conditions, most prominently the light–dark shift. During the day, light is available as copious source of energy, allowing photosynthetic carbon assimilation, whereas during the night, the breakdown of internal starch delivers the necessary energetic power . From a metabolic engineering perspective, knowledge on both physiological states is essential to optimize plants in a way that carbon is channeled optimally towards desired traits throughout the diurnal cycle. Therefore, these two fundamental growth states, light and dark metabolism, were now studied using elementary flux mode analysis.
The full sets of elementary flux modes, which span the entire space of feasible flux distributions, were calculated both for metabolism in the light using CO2 as sole carbon source (light-metabolism), and for metabolism in the dark using internal starch (dark-metabolism). The solution space consisted of 1.2 million elementary flux modes for the light-metabolism and of 5.7 million elementary flux modes for the dark-metabolism (Table 1). The theoretical maximum growth yields were 27.4 and 28.6 (g biomass) (C-mol substrate)−1 for the two scenarios. Assuming equal contribution of day and night metabolism, the resulting mean value of 28.1 (g biomass) (C-mol substrate)−1 closely resembles experimental values for A. thaliana leaves of about 25.2 (g biomass) (C-mol substrate)−1 . The maximum theoretical biomass formation, predicted by the model, revealed a high carbon efficiency of plant metabolism: 86.8 to 90.6 % of the assimilated carbon was incorporated into biomass (Table 2). For prolonged illumination phases, the simulated values further approached the measured ones (Additional file 1: Table S6). This was a first indication that in vivo plant metabolism closely approaches theoretical optimum performance with regard to stoichiometric capacity, which seems a particularly useful characteristic with regard to use in production. When the light influx was omitted from the model, no modes resulted for CO2 assimilation, deducing that degradation of a heterotrophic substrate, such as starch, is the only feasible phenotype in the dark.
A narrow range of absorbed photons supports optimal plant growth
The computed set of elementary modes was used to investigate the impact of the light influx on growth in more detail. First, energetically inefficient modes that included massive cycling of resources were eliminated. Although stoichiometrically possible, these modes involve futile cycling of carbon across the internal membranes to dissipate energy (Additional file 1: Figure S1). They are an emergent property, which arises from the compartmentation of metabolic pathways into intracellular organelles. Because such massive futile cycling is energetically inefficient, they are considered both evolutionary and physiologically implausible [88, 89]. Furthermore, the high reaction rates, observed in such elementary modes, would require the additional synthesis of large amounts of protein with further energetic burden, further justifying the elimination of these modes from the subsequent evaluation. This identified a physiologically plausible subset of 848,629 individual flux modes. When these modes were mapped against the corresponding acquired light influx for each of the modes (Fig. 3f), the production of biomass was found highest in a defined range between 8.6 and 11.5 (mol photons) (mol CO2)−1. This corresponds to a quantum yield between 0.09 - 0.12 (mol CO2) (mol photons)−1, obviously supporting optimal photosynthesis.
Flux rearrangement between optimum day and night metabolism requires reversible translocation of carbohydrates
Next, day-time and night-time metabolism were evaluated on the level of intracellular fluxes. The flux distributions, which resemble optimum growth, revealed large differences between the two conditions (Figs. 4 and 5). Shortly, dark metabolism involved the degradation of starch into maltose, which was subsequently transported across the plastid membrane (Fig. 4). In the cytosol, maltose was hydrolyzed into two glucose molecules, followed by phosphorylation to glucose 6-phosphate. The sugar phosphate was then channeled into the cytosolic and the plastidic EMP pathways. The glucose 6-phosphate translocator mediated the distribution of carbon between the different compartments. In this way, both malate and pyruvate were generated through malate dehydrogenase and through the EMP pathways, respectively, and were subsequently imported into the mitochondrion to fuel the TCA cycle. The light-metabolism assimilated atmospheric CO2 via the CBB cycle in the plastid (Fig. 4). This involved a high flux through the CBB cycle and the non-oxidative PP pathway, as both routes are strongly entwined. Mainly, the triose-phosphate translocators supplied carbon-building blocks into the cytosol, whereby the glucose 6-phosphate translocator recycled carbon back into the plastid. Mostly, malate was built from initial carbon conversions and it was then imported into the mitochondrion. Malic enzyme provided the TCA cycle with pyruvate. In this regard, the translocation of pyruvate between the different compartments was one prominent example of a fully reversed flux, as described above. To conclude, in the light, especially the plastidic metabolism was highly active and the triose phosphate carriers transported carbon out of the plastid. Mainly malate was imported into the mitochondrion, whereas in the dark both malate and pyruvate were channeled across the mitochondrial membrane and the cytosol showed slightly more activity than the plastid. The main carbon export from the plastid occurred in the form of maltose. The most significant differences between light and dark metabolism were found for the plastid metabolism, for the energy metabolism and for specific inter-compartmental transporters (Fig. 5). In contrast, especially the mitochondrial and the peroxisomal reactions remained rather unchanged. The simulation data predict an on/off shift for five previously postulated diurnal switches, i.e. ribulose 1,5-bisphosphate carboxylase, phosphoribulokinase, fructose 1,6-bisphosphate phosphatase, NADP-dependent glyceraldehyde 3-phosphate dehydrogenase and sedoheptulose 1,7-bisphosphate phosphatase . In addition, many more changes were observed. A reversal in flux direction resulted for no less than 19 reactions (Additional file 1: Table S7, Fig. 5), whereas 70 reactions kept the direction of flux, but showed significant changes in flux value. From the latter, 66 fluxes were significantly decreased for the dark metabolism as compared to the light metabolism, whereas 23 reactions showed a significant increase (Fig. 5). Given the fact, that the analysis was on the averaged flux behavior of the top 1 % biomass producing modes, it cannot be excluded that alternative solutions exist within the flux solution space that do not behave in the same way as the here described flux changes (Additional file 1: Figure S3).
Simulation quantifies differences in energy and redox supply between light and dark metabolism
As shown in the previous section, the metabolism in the dark recruited all compartments to supply ATP and reducing power, whereas the plastid was the dominant compartment in the light (Fig. 4). Next, we analyzed in more detail, which metabolic processes are responsible for the supply of energy and reducing power. The analysis again considered the optimal dark and light fluxes, i.e. the average of the top 1 % biomass producing modes. Illuminated leaves largely delivered NADPH and ATP through the photosynthetic light reactions in the plastid (Fig. 6), whereas the NADH stemmed mostly from cytosolic malate dehydrogenase, i.e. was exported as NADH from the plastid . During the dark, ATP was produced predominantly via oxidative phosphorylation in the mitochondrion (60 %), and reduced NADPH was generated by all compartments. Secondly, the formation and assimilation of carbon dioxide was different between the two physiological states (Fig. 6). Generally, the total net flux of CO2 production was higher for the dark-metabolism (10 %), as compared to light-metabolism (5 %), whereby overall, both values were rather low. It was interesting to note that also the origin of the released CO2 differed. Under both conditions, the carbon loss by photorespiration was negligible, whereas the largest contributors to CO2 release were the TCA cycle and pyruvate and malate dehydrogenase. During the dark, additional carbon was released during the generation of NADPH in the oxidative PP pathway, whereas during illumination, previously assimilated CO2 was lost via PEP carboxykinase.
Optimum anabolic pathway use switches between day and night
The striking differences in the localization of energy and reducing power supply among the cellular compartments now suggested to also inspecting anabolic metabolism. Therefore, the metabolic pathways, responsible for the synthesis of each anabolic precursor, were investigated and compared between the optimum growth modes in the light and in the dark (average flux values of the top 1 % biomass producing modes). Interestingly, the biosynthetic origin of many of the twelve anabolic precursors strongly depended on the physiological growth mode, i.e. light and carbon acquisition (Fig. 7). As example, the majority of 3-phosphoglycerate in the light was supplied via the CBB cycle, whereas, in the dark, the precursor exclusively stemmed from the EMP pathway. In addition, ribose 5-phosphate was produced through the non-oxidative PP pathway under illumination, whereas, in the dark, it was mainly generated by the oxidative PP pathway. For other precursors, such as α-ketoglutaric acid, the metabolic network recruited isoenzymes, which differed in cofactor use, depending on the illumination conditions.
Pathway fluxes remain stable upon photosynthetic plasticity
It is generally accepted that cyclic and non-cyclic electron flow both contribute to the in vivo photosynthetic light reaction [65–67]. However, it is unclear at which ratio cyclic and non-cyclic electron flow co-operate. It is hypothesized that the observed flexibility might allow the modulation of the ATP: NADPH ratio to match the demand of metabolism under changing environmental conditions [66, 67], as non-cyclic electron flow allows for simultaneous NADPH and ATP synthesis, whereas cyclic electron flow solely generates ATP. It is prevalent in the literature that under environmental stress, such as cold, light or osmotic stress, certain metabolic changes arise, that modulate the cellular ATP: NADPH demand [91, 92]. In this regard, it has been acknowledged, that this modulation requires fine-tuned regulation of the contribution of cyclic electron flow to photosynthesis . Reversely, sudden light fluctuations, such as shading by a cloud or overhanging leaf during full sunlight at noon, directly cause changes in the contribution of cyclic electron flow to photosynthesis . In particular, an increase in cyclic electron flow appears important in protection of photosystem I against photo inhibition [94, 95]. It was now interesting to investigate the interplay between the metabolic pathways in meeting different ratios for NADPH and ATP, as expected under sudden fluctuations in light. For this purpose, optimal flux distributions (In average flux values of the top 1 % biomass producing modes) of light-metabolism were calculated for 15 scenarios with increasing contribution of cyclic electron flow to the overall electron flow, i.e. an increasing ATP:NADPH ratio. Ratios between non-cyclic and cyclic electron flow of 14:0, 13:1, 12:2, 11:3, 10:4, 9:5, 8:6, 7:7, 6:8, 5:9, 4:10, 3:11, 2:12, 1:13 and 0:14 were considered (Additional file 1: Table S8). For all ratios, except for 100 % cyclic electron flow, the same theoretical maximum for biomass formation was found. Omission of non-cyclic electron flow did not allow for any biomass formation. Three general trends could be observed when plotting the absolute flux values for optimal growth against the ratio between cyclic and non-cyclic electron flow (Fig. 8). Either, (i) optimal flux values were not influenced by the type of phosphorylation (e.g. ribose 5-phosphate isomerase) or an exponential increase at either (ii) extremely low or (iii) extremely high contribution of cyclic electron flow was observed (e.g. mitochondrial phosphate import and glyceraldehyde 3-phosphate dehydrogenase, respectively). Between 20 and 70 % cyclic electron flow contribution all fluxes remained rather unchanged, which can be taken as an indication for a highly flexible plant metabolism, easily capable of adjusting to different ATP:NADPH ratios. The specific metabolic flux changes that allow to cope with an increase in cyclic electron flow, constitute futile cycling of metabolites to actively dissipate an excess in ATP (Additional file 1: Figure S2). Examples include, substrate cycling between phosphoenolpyruvate and pyruvate, substrate cycling between fructose 1,6-bisphosphate and fructose 6-phosphate, cycling across the plastid membrane mediated by 3-phosphoglycerate and glucose 6-phospate transport and the activity of the cytosolic NADPH-dependent glyceraldehyde 3-phosphate dehydrogenase. Furthermore, the metabolic fluxes that are most influenced by a decreased ATP: NADPH ratio, appear to handle an excess of redox power by channeling NADPH through the malate/oxaloacetate shuttle into the mitochondrion, where ATP is produced by oxidative phosphorylation (Additional file 1: Figure S2). This supports the hypothesis that mitochondria can act as a sink for reduced NADPH [96–98]. In addition, plastidic triose phosphate is increasingly exported to the cytosolic EMP pathway, thus effectively reducing the plastidic ATP requirement.
Systematic analysis of the elementary mode solution space unravels the light stress response
The constructed metabolic network enabled the adequate description of fundamental physiological traits, previously found experimentally for A. thaliana, such as growth and photosynthetic efficiency, which can be taken as indication of its high quality and validity (Fig. 3 and Table 2). Obviously, Arabidopsis rosettes operate close to optimum. This becomes evident from the resemblance between the predicted theoretical maximum and the experimental growth yield (Table 2) and the predicted optimum range of photon absorbance, which rather exactly match with the range observed in vivo (Fig. 3) [39–43].
Optimal growth only results for a narrow range of photon absorbance, while higher and lower light intensities reduce growth efficiency, respectively (Fig. 3). Possible molecular reasons for this phenomenon can be directly extracted from the flux mode solution space. The direct product of photo-reduction is NADPH (Fig. 3b). Under high light intensity, its production soon exceeds the biosynthetic need. A physiologically feasible way to cope with such NADPH excess is an increased flux through the photo-respiratory pathway, as photorespiration consumes NADPH. Interestingly, a clear relationship between the flux through the photo respiratory system (glycerate kinase) and phosphoribulokinase, and the quantum requirement under high light intensity suggests that indeed Arabidopsis deals with light stress by increasing photorespiration (Fig. 3d-e). A direct consequence of such an up-regulation is reduced growth, because photorespiration is associated with a net loss of carbon as CO2 [99–101]. Recently, an activated photorespiration has been experimentally recognized in plants as an important light stress response to dissipate excess reducing equivalents and energy . Alternatively, NADPH excess could also be handled by an increased consumption of ATP, e.g. through futile cycles. From a metabolic viewpoint, this would involve a conversion of excess NADPH into NADH by a transhydrogenase-like reaction, and fueling the phosphorylation of ADP into ATP. In this way, growth could be maintained at its optimal rate, however in an energetically less efficient way.
The intracellular pathways of Arabidopsis rosettes operate in vivo near the predicted optimum flux distribution
The integration of in silico and in vivo fluxes allows evaluation of the observed cellular physiology within the overall feasible flux space, as proven valuable for different microorganisms [10, 103]. Fortunately, recent 13C flux analysis data from Arabidopsis rosettes  provide an excellent opportunity to conduct such an integration for the first time in plant cells. For this purpose, the determined flux values were now mapped with the predicted fluxes for optimum growth (Fig. 9). In this context, ‘optimum growth’ refers to the average flux value of the top 1 % biomass producing modes. This provided a striking agreement between the in vivo and optimum in silico fluxome for the entire plant. Considering the similarity of measured and predicted biomass and quantum yield, we conclude that illuminated thale cress is represented adequately by the proposed metabolic model, and its leaves operate very close to their maximal potential with regard to stoichiometric capacity. It should be noted that this may not hold for suboptimal growth conditions, such as high light acclimation , where the stress conditions seem to induce a metabolic burden. The observed high metabolic efficacy of plants is particularly interesting in the light of their biotechnological potential as only little carbon is ‘wasted’ on respiration and fueling of metabolism. A similar comparison for heterotrophic leaf metabolism could not be drawn due to a lack of experimental data. Unfortunately, whole plant in vivo flux studies in the dark are currently still unattainable .
Synthesis of compositional traits could be driven by dynamic metabolic engineering
Clearly, our simulations indicate that a reversible transhydrogenase-like function is crucial to rearrange the metabolism from day to night (Fig. 2). It is interesting to note that dark and light metabolism involves a rather diverse set of anabolic pathways for optimum precursor supply (Fig. 7). This also involves isoenzymes with different prevalence for NADPH and NADH, which vary between the growth regimes. So far, their role is still unclear, however functional diversification  and increased metabolic robustness  have been postulated as possible purposes of isoenzymes. Based on our simulations, the ubiquitous presence of isoenzymes with different cofactor usage in plants seems key for flexible handling of specific anabolic demands of day and night metabolism. The extended knowledge on flux re-arrangement between day and night physiologies is particularly interesting for the development of plants with biotechnologically interesting traits. Plants recruit distinct pathways for the synthesis of anabolic precursors and energy during the day and at night (Figs. 6 and 7). This information seems most helpful to optimize precursor, reducing power and energy supply towards biotechnologically interesting traits. On a first glance, constitutive synthesis would permit continuous accumulation of the desired compound. However, biosynthesis might require much more energy or reducing equivalents during the night, as compared to the day. In such cases, it might be more desirable to connect biotechnological syntheses to genes that are tightly regulated throughout the diurnal cycle. Novel approaches, which allow dynamic programming of metabolism , seem straightforward to exploit this fundamental plant characteristic.
Day-night flux rearrangement is superimposed by selected transcriptional changes
As shown, the optimum metabolism of Arabidopsis recruits a different set of pathways during day and night: diurnal cycling is linked to a substantial redistribution of flux. It was now interesting to see, how regulatory circuits of the plant superimpose the obvious fine-adjustment on the level of the metabolic network, i.e. to which extent transcriptional and post-transcriptional control is involved, and which genes are part of the overarching regulatory network. Our simulation data enabled the integration of fluxes with the expression of encoding genes. This integration of transcriptomic knowledge with the observed flux changes under optimum growth conditions offered a systems level view of the day-night shift (Fig. 10). Reactions, such as the ones catalyzed by mitochondrial and plastidic pyruvate dehydrogenase and citrate synthase were found stable in either fluxome or transcriptome. Others displayed clear diurnal trends, as for instance cytosolic and plastidic pyruvate kinase, phosphoribulokinase, fructose 1,6-bisphosphate phosphatase and many others. Diurnal changes in both transcripts and fluxes were observed for the previously identified on/off switches, i.e. for ribulose 1,5-bisphosphate carboxylase, phosphoribulokinase, fructose 1,6-bisphosphate phosphatase, NADP-dependent glyceraldehyde 3-phosphate dehydrogenase and sedoheptulose 1,7-bisphosphate phosphatase . In addition, the triose phosphate translocator between plastid and cytosol is known to follow a diurnal path for carbon export . This supports the predicted high export flux for 3-phosphoglycerate, phosphoenolpyruvate and dihydroxyacetone phosphate during the light period. Additionally, light/dark modulation of pyruvate-orthophosphate (Pi) dikinase through phosphorylation in C3 leaves  is confirmed by an off switch during the night. Even for those genes that appear to have unaltered expression levels, although fluxomic changes are observed, isogenes exist that are still to be tested, including reactions of the oxidative PP pathway, fumarase or hexose isomerase. Overall, systems level correlation occurs plentiful and seems well distributed over the entire central carbon metabolism, reflecting the necessity of readjusting the core metabolism during day/night-transition.
Dynamics in metabolic pathway fluxes and gene expression
A more detailed view on selected transcriptional changes is obtained by overlaying the time resolved expression of metabolic genes from Arabidopsis grown with day night cycles  with the corresponding flux changes of the encoded reactions, as predicted from the simulations under optimal growth (Fig. 10). In a number of cases, transcription-based and flux-based changes exhibit a close connection, which indicates that the plant strongly recruits regulatory mechanisms to drive core metabolism. Interestingly, the genes that showed analogous behavior in transcript and flux, were located around three main controlling points. Firstly, the fluxes between fructose 6-phosphate and triose-phosphate appeared somewhat regulated. At least, fructose 1,6-bisphosphatase and fructose bisphosphate aldolase showed similar behavior in transcript and flux. Secondly, many enzymes involving malate and pyruvate had a diurnal pattern, indicating malate and/or pyruvate might be a second controlling point in metabolism. Both the fructose 1,6-bisphosphatase and pyruvate kinase controlling points have previously been postulated as diurnal regulators of metabolism . In addition, the complex regulation pattern around the malate node is reflected in the local high network complexity. Finally, genes associated with the TCA cycle seem superimposed by transcriptional regulation.
Besides identifying clusters of potential transcriptional control, this type of analysis can be used to identify, which isoenzymes might be linked to flux changes. For example, for three pyruvate kinase-associated genes (At5g08570, At5g56350 and At3g49160) a diurnal expression pattern was identified, however only the former two are linked to the observed flux changes. This indicates that At3g49160 likely has a functionally different purpose from mediating the pyruvate phosphorylation in central carbon metabolism. For instance, pyruvate kinase could also catalyze several metabolic conversions in ribonucleotide and nucleotide biosynthesis. As isoenzymes might have distinct functions , correlation of flux and transcript could pose a valuable method for identification of iso-enzymatic function and might provide first evidence on their biochemical function. This is particularly useful in cases, where a high number of isoenzymes potentially catalyze a certain reaction. Furthermore, flux-transcript correlation can assist in the confirmation of genes with putative function. Putative genes that do not show significant linkage to flux change possibly do not control flux on a transcriptional level (At3g15020, At2g26080, At1g58150, At3g49160, At5g11670). However, those that do correlate might present candidates involved in fructose bisphosphate aldolase (At2g36460), fructose 1,6-bisphosphatase (At1g43670) and phosphoenolpyruvate carboxylase (At1g53310). One should keep in mind for these interpretations, that not in all cases transcript changes will immediately lead to protein change, because translation in plants can be damped or delayed .
Flux-homeostasis and biotechnological impact of plasticity in photophosphorylation
Plants adapt to changes in light intensity through photo protection and optimization of energy conversion . In this way, also the output ratio of ATP: NADPH can be influenced significantly by the light environment, which has been investigated extensively [65–67, 96, 101]. How cellular metabolism copes with the changed ATP: NADPH supply on the level of intracellular fluxes is investigated here. For instance, it seems that the protection of photosystem I against photo-inhibition through an increase in cyclic electron flow, only poses a small metabolic burden, as only little modulation is necessary to deal with the imposed increase in ATP (Fig. 8). The short time frame, in which metabolism needs to react to such sudden fluctuations in light, is also reflected in the fact that the quick onset of large differences in ATP, only requires small flux changes in response. On the short term, these changes might not require modulation of transcript or protein levels, and could therefore warrant an instantaneous short-term response. Additionally, the observed stability of biomass formation, with only few distinct flux changes, across a wide range of ATP: NADPH ratios (Fig. 8), indicates a certain robustness of metabolism to environmental changes. The metabolic fluxes that are most influenced by small changes around the assumed in vivo ratio of non-cyclic and cyclic electron flow of 12:2, involve specific modulation of the NADPH status and of the dissipation of ATP through futile cycling. Possibly, these specific flux changes permit homeostasis of net flux and adenylate/redox status. A homeostasis of the adenylate status during photosynthesis in a fluctuating environment has previously been indicated, however, here, homeostasis was supposedly attributed to changes in pathway usage . Here, we observe that the net flux distribution remains unchanged under abruptly changing light environment, however, strong differences occur in substrate cycling. Such substrate cycling has previously been identified in vivo in plant tissues using experiments  and in silico through modelling . Likely, the observed metabolic adaptation through futile cycling occurs in addition to changes in pathway use, such as increased photorespiration. This improves the plants’ capacity to cope with a constantly changing light environment, both on the short-term through e.g. ATP dissipation and redox modulation and on the long-term through e.g. changes in pathway usage. Furthermore, it has recently been proposed that the excess in redox power could be directed towards light-driven production of biotechnological compounds through cytochrome P450s-mediated reactions . When an organism is engineered to produce large amounts of a biotechnologically interesting product, its molecular flux patterns change. These metabolic changes engage a different demand for ATP and NADPH that needs to be accustomed by the cell. Due to the observed plasticity in non-cyclic and cyclic electron flow, plants are capable to adapt to such modulated energetic requirements. In addition, our simulations now show that such adaptations do not impede growth, granting plants a high potential for the production of biotechnological products, especially for those compounds requiring much redox power. This emphasizes the potential of photosynthetic light reactions in biotechnology.
The created condensed network model adequately described leaf physiology, and thus, because of its reduced network size as compared to genome-scale models, permits straightforward, comprehensive analysis of the entire elementary flux mode solution space. Taken together, the metabolic simulations provide detailed molecular insights into plant functioning. Arabidopsis can operate close to theoretical pathway optimum and that this is mediated by a fine-adjustment of metabolic flux, strongly under transcriptional control. In this light, the present work is one of the very few examples so far, which link in-vivo with in-silico flux data to a higher-level understanding [103, 114, 115]. It seems straightforward to extend this to other plant systems and to more specific models that address specific plant tissues, which are formed during plant development. The knowledge gained from our systems-biological approach, together with the high potential of plants as biotechnological production platforms, especially for compounds requiring much redox power, will contribute to establish plants as biotechnological factories.
Yan L, Kerr PS. Genetically engineered crops: their potential use for improvement of human nutrition. Nutr Rev. 2002;60(5 Pt 1):135–41.
Dyer JM, Stymne S, Green AG, Carlsson AS. High-value oils from plants. Plant J. 2008;54(4):640–55.
Saha S, Ramachandran S. Genetic improvement of plants for enhanced bio-ethanol production. Recent Pat DNA Gene Seq. 2013;7(1):36–44.
Rajasekaran A, Kalaivani M. Designer foods and their benefits: A review. J Food Sci Technol. 2013;50(1):1–16.
Junker BH. Flux analysis in plant metabolic networks: increasing throughput and coverage. Curr Opin Biotechnol. 2014;26:183–8.
Shachar-Hill Y. Metabolic network flux analysis for engineering plant systems. Curr Opin Biotechnol. 2013;24(2):247–55.
Cusido RM, Onrubia M, Sabater-Jara AB, Moyano E, Bonfill M, Goossens A, Angeles Pedreño M, Palazon J. A rational approach to improving the biotechnological production of taxanes in plant cell cultures of Taxus spp. Biotechnol Adv. 2014;32(6):1157–67.
Sweetlove LJ, Last RL, Fernie AR. Predictive metabolic engineering: a goal for systems biology. Plant Physiol. 2003;132(2):420–5.
Lotz K, Hartmann A, Grafahrend-Belau E, Schreiber F, Junker BH. Elementary flux modes, flux balance analysis, and their application to plant metabolism. Methods Mol Biol. 2014;1083:231–52.
Becker J, Zelder O, Häfner S, Schröder H, Wittmann C. From zero to hero--design-based systems metabolic engineering of Corynebacterium glutamicum for L-lysine production. Metab Eng. 2011;13(2):159–68.
Poblete-Castro I, Binger D, Rodrigues A, Becker J. Martins Dos Santos VA, Wittmann C: In-silico-driven metabolic engineering of Pseudomonas putida for enhanced production of poly-hydroxyalkanoates. Metab Eng. 2013;15:113–23.
Kim B, Park H, Na D, Lee SY. Metabolic engineering of Escherichia coli for the production of phenol from glucose. Biotechnol J. 2014;9(5):621–9.
Trinh CT. Elucidating and reprogramming Escherichia coli metabolisms for obligate anaerobic n-butanol and isobutanol production. Appl Microbiol Biotechnol. 2012;95(4):1083–94.
Hwang KS, Kim HU, Charusanti P, Palsson BØ, Lee SY. Systems biology and biotechnology of Streptomyces species for the production of secondary metabolites. Biotechnol Adv. 2014;32(2):255–68.
Michael TP, Jackson S. The First 50 Plant Genomes. Plant Genome. 2013;6. doi:10.3835/plantgenome2013.03.0001in.
Grafahrend-Belau E, Schreiber F, Koschützki D, Junker BH. Flux balance analysis of barley seeds: a computational approach to study systemic properties of central metabolism. Plant Physiol. 2009;149(1):585–98.
de Oliveira Dal’Molin CG, Quek LE, Palfreyman RW, Brumbley SM, Nielsen LK. AraGEM, a genome-scale reconstruction of the primary metabolic network in Arabidopsis. Plant Physiol. 2010;152(2):579–89.
Saha R, Suthers PF, Maranas CD. Zea mays iRS1563: a comprehensive genome-scale metabolic reconstruction of maize metabolism. PLoS ONE. 2011;6(7):e21784.
Dersch LM, Beckers V, Wittmann C. Green pathways: Metabolic network analysis of plant systems. Metab Eng. 2016;34:1-24.
Schuster S, Dandekar T, Fell DA. Detection of elementary flux modes in biochemical networks: a promising tool for pathway analysis and metabolic engineering. Trends Biotechnol. 1999;17(2):53–60.
Terzer M, Stelling J. Large-scale computation of elementary flux modes with bit pattern trees. Bioinformatics. 2008;24(19):2229–35.
Papin JA, Price ND, Palsson BØ. Extreme pathway lengths and reaction participation in genome-scale metabolic networks. Genome Res. 2002;12(12):1889–900.
Price ND, Papin JA, Palsson BØ. Determination of redundancy and systems properties of the metabolic network of Helicobacter pylori using genome-scale extreme pathway analysis. Genome Res. 2002;12(5):760–9.
Sauer U. Metabolic networks in motion: 13C-based flux analysis. Mol Syst Biol. 2006;2:62.
Wittmann C. Fluxome analysis using GC-MS. Microb Cell Fact. 2007;6:6.
Young JD, Shastri AA, Stephanopoulos G, Morgan JA. Mapping photoautotrophic metabolism with isotopically nonstationary 13C flux analysis. Metab Eng. 2011;13(6):656–65.
Dersch LM, Beckers V, Rasch D, Melzer G, Bolten C, Kiep K, Becker H, Blasing OE, Fuchs R, Ehrhardt T, et al. Novel Approach for High-Throughput Metabolic Screening of Whole Plants by Stable Isotopes. Plant Physiol. 2016;171(1):25–41.
Wang L, Birol I, Hatzimanikatis V. Metabolic control analysis under uncertainty: framework development and case studies. Biophys J. 2004;87(6):3750–63.
Van Norman JM, Benfey PN. Arabidopsis thaliana as a model organism in systems biology. Wiley Interdiscip Rev Syst Biol Med. 2009;1(3):372–9.
Mintz-Oron S, Meir S, Malitsky S, Ruppin E, Aharoni A, Shlomi T. Reconstruction of Arabidopsis metabolic network models accounting for subcellular compartmentalization and tissue-specificity. Proc Natl Acad Sci U S A. 2012;109(1):339–44.
Orzechowski S. Starch metabolism in leaves. Acta Biochim Pol. 2008;55(3):435–45.
Fernie AR, Carrari F, Sweetlove LJ. Respiratory metabolism: glycolysis, the TCA cycle and mitochondrial electron transport. Curr Opin Plant Biol. 2004;7(3):254–61.
Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, Kawashima S, Okuda S, Tokimatsu T, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36(Database issue):D480–4.
Grafahrend-Belau E, Weise S, Koschutzki D, Scholz U, Junker BH, Schreiber F. MetaCrop: a detailed database of crop plant metabolism. Nucleic Acids Res. 2008;36(Database issue):D954–8.
Mueller LA, Zhang P, Rhee SY. AraCyc: a biochemical pathway database for Arabidopsis. Plant Physiol. 2003;132(2):453–60.
Melzer G, Esfandabadi ME, Franco-Lara E, Wittmann C. Flux Design: In silico design of cell factories based on correlation of pathway fluxes to desired properties. BMC Syst Biol. 2009;3:120.
Bläsing OE, Gibon Y, Günther M, Höhne M, Morcuende R, Osuna D, Thimm O, Usadel B, Scheible WR, Stitt M. Sugars and circadian regulation make major contributions to the global regulation of diurnal gene expression in Arabidopsis. Plant Cell. 2005;17(12):3257–81.
Szecowka M, Heise R, Tohge T, Nunes-Nesi A, Vosloh D, Huege J, Feil R, Lunn J, Nikoloski Z, Stitt M, et al. Metabolic fluxes in an illuminated Arabidopsis rosette. Plant Cell. 2013;25(2):694–714.
Ehleringer J, Bjorkman O. Quantum Yields for CO2 Uptake in C3 and C4 Plants: Dependence on Temperature, CO2, and O2 Concentration. Plant Physiol. 1977;59(1):86–90.
Ehleringer J, Pearcy RW. Variation in Quantum Yield for CO2 Uptake among C3 and C4 Plants. Plant Physiol. 1983;73(3):555–9.
Long SP, Postl WF, Bolharnordenkampf HR. Quantum Yields for Uptake of Carbon-Dioxide in C-3 Vascular Plants of Contrasting Habitats and Taxonomic Groupings. Planta. 1993;189(2):226–34.
Osborne BA, Garrett MK. Quantum yield for CO2 uptake in some diploid and tetraploid plant species. Plant Cell Environ. 1983;6:135–44.
Skillman JB. Quantum yield variation across the three pathways of photosynthesis: not yet out of the dark. J Exp Bot. 2008;59(7):1647–61.
Hall DO, Rao K. Photosynthesis. Cambridge/New York: Cambridge University Press; 1999.
Kruger NJ, von Schaewen A. The oxidative pentose phosphate pathway: structure and organisation. Curr Opin Plant Biol. 2003;6(3):236–46.
Plaxton WC. The Organization and Regulation of Plant Glycolysis. Annu Rev Plant Physiol Plant Mol Biol. 1996;47:185–214.
Schnarrenberger C, Martin W. Evolution of the enzymes of the citric acid cycle and the glyoxylate cycle of higher plants. A case study of endosymbiotic gene transfer. Eur J Biochem. 2002;269(3):868–83.
Sweetlove LJ, Fernie AR. The spatial organization of metabolism within the plant cell. Annu Rev Plant Biol. 2013;64:723–46.
Raghavendra AS, Reumann S, Heldt HW. Participation of mitochondrial metabolism in photorespiration. Reconstituted system of peroxisomes and mitochondria from spinach leaves. Plant Physiol. 1998;116(4):1333–7.
Visser WF, van Roermund CW, Ijlst L, Waterham HR, Wanders RJ. Metabolite transport across the peroxisomal membrane. Biochem J. 2007;401(2):365–75.
Fatland BL, Nikolau BJ, Wurtele ES. Reverse genetic characterization of cytosolic acetyl-CoA generation by ATP-citrate lyase in Arabidopsis. Plant Cell. 2005;17(1):182–203.
Laloi M. Plant mitochondrial carriers: an overview. Cell Mol Life Sci. 1999;56(11–12):918–44.
Klingenberg M. The ADP and ATP transport in mitochondria and its carrier. Biochim Biophys Acta. 2008;1778(10):1978–2021.
Haferkamp I, Fernie AR, Neuhaus HE. Adenine nucleotide transport in plants: much more than a mitochondrial issue. Trends Plant Sci. 2011;16(9):507–15.
Picault N, Hodges M, Palmieri L, Palmieri F. The growing family of mitochondrial carriers in Arabidopsis. Trends Plant Sci. 2004;9(3):138–46.
Furumoto T, Yamaguchi T, Ohshima-Ichie Y, Nakamura M, Tsuchida-Iwata Y, Shimamura M, Ohnishi J, Hata S, Gowik U, Westhoff P, et al. A plastidial sodium-dependent pyruvate transporter. Nature. 2011;476(7361):472–5.
Fischer K. The import and export business in plastids: transport processes across the inner envelope membrane. Plant Physiol. 2011;155(4):1511–9.
Eicks M, Maurino V, Knappe S, Flugge UI, Fischer K. The plastidic pentose phosphate translocator represents a link between the cytosolic and the plastidic pentose phosphate pathways in plants. Plant Physiol. 2002;128(2):512–22.
Facchinelli F, Weber AP. The metabolite transporters of the plastid envelope: an update. Front Plant Sci. 2011;2:50.
Schwender J, Shachar-Hill Y, Ohlrogge JB. Mitochondrial metabolism in developing embryos of Brassica napus. J Biol Chem. 2006;281(45):34040–7.
Kaldenhoff R, Kai L, Uehlein N. Aquaporins and membrane diffusion of CO2 in living organisms. Biochim Biophys Acta. 2014;1840(5):1592–5.
Terashima I, Fujita T, Inoue T, Chow WS, Oguchi R. Green light drives leaf photosynthesis more efficiently than red light in strong white light: revisiting the enigmatic question of why leaves are green. Plant Cell Physiol. 2009;50(4):684–97.
Rausch C, Bucher M. Molecular mechanisms of phosphate transport in plants. Planta. 2002;216(1):23–37.
Cheung CY, Williams TC, Poolman MG, Fell DA, Ratcliffe RG, Sweetlove LJ. A method for accounting for maintenance costs in flux balance analysis improves the prediction of plant cell metabolic phenotypes under stress conditions. Plant J. 2013;75(6):1050–61.
Allen JF. Cyclic, pseudocyclic and noncyclic photophosphorylation: new links in the chain. Trends Plant Sci. 2003;8(1):15–9.
Cruz JA, Avenson TJ, Kanazawa A, Takizawa K, Edwards GE, Kramer DM. Plasticity in light reactions of photosynthesis for energy production and photoprotection. J Exp Bot. 2005;56(411):395–406.
Kramer DM, Avenson TJ, Edwards GE. Dynamic flexibility in the light reactions of photosynthesis governed by both electron and proton transfer reactions. Trends Plant Sci. 2004;9(7):349–57.
Foyer CH, Neukermans J, Queval G, Noctor G, Harbinson J. Photosynthetic control of electron transport and the regulation of gene expression. J Exp Bot. 2012;63(4):1637–61.
Munekage YN, Genty B, Peltier G. Effect of PGR5 impairment on photosynthesis and growth in Arabidopsis thaliana. Plant Cell Physiol. 2008;49(11):1688–98.
Krömer JO, Wittmann C, Schröder H, Heinzle E. Metabolic pathway analysis for rational design of L-methionine production by Escherichia coli and Corynebacterium glutamicum. Metab Eng. 2006;8(4):353–69.
Reiter WD, Chapple C, Somerville CR. Mutants of Arabidopsis thaliana with altered cell wall polysaccharide composition. Plant J. 1997;12(2):335–45.
Zablackis E, Huang J, Muller B, Darvill AG, Albersheim P. Characterization of the cell-wall polysaccharides of Arabidopsis thaliana leaves. Plant Physiol. 1995;107(4):1129–38.
Herrero J, Fernández-Pérez F, Yebra T, Novo-Uzal E, Pomar F, Pedreño MA, Cuello J, Guéra A, Esteban-Carrasco A, Zapata JM. Bioinformatic and functional characterization of the basic peroxidase 72 from Arabidopsis thaliana involved in lignin biosynthesis. Planta. 2013;237(6):1599–612.
Tschoep H, Gibon Y, Carillo P, Armengaud P, Szecowka M, Nunes-Nesi A, Fernie AR, Koehl K, Stitt M. Adjustment of growth and central metabolism to a mild but sustained nitrogen-limitation in Arabidopsis. Plant Cell Environ. 2009;32(3):300–18.
Fan J, Yan C, Zhang X, Xu C. Dual role for phospholipid:diacylglycerol acyltransferase: enhancing fatty acid synthesis and diverting fatty acids from membrane lipids to triacylglycerol in Arabidopsis leaves. Plant Cell. 2013;25(9):3506–18.
Shen W, Li JQ, Dauk M, Huang Y, Periappuram C, Wei Y, Zou J. Metabolic and transcriptional responses of glycerolipid pathways to a perturbation of glycerol 3-phosphate metabolism in Arabidopsis. J Biol Chem. 2010;285(30):22957–65.
Stahl U, Carlsson AS, Lenman M, Dahlqvist A, Huang B, Banas W, Banas A, Stymne S. Cloning and functional characterization of a phospholipid:diacylglycerol acyltransferase from Arabidopsis. Plant Physiol. 2004;135(3):1324–35.
Saxena SC, Salvi P, Kaur H, Verma P, Petla BP, Rao V, Kamble N, Majee M. Differentially expressed myo-inositol monophosphatase gene (CaIMP) in chickpea (Cicer arietinum L.) encodes a lithium-sensitive phosphatase enzyme with broad substrate specificity and improves seed germination and seedling growth under abiotic stresses. J Exp Bot. 2013;64(18):5623–39.
Nowicka B, Strzalka W, Strzalka K. New transgenic line of Arabidopsis thaliana with partly disabled zeaxanthin epoxidase activity displays changed carotenoid composition, xanthophyll cycle activity and non-photochemical quenching kinetics. J Plant Physiol. 2009;166(10):1045–56.
Tardy F, Havaux M. Photosynthesis, chlorophyll fluorescence, light-harvesting system and photoinhibition resistance of a zeaxanthin-accumulating mutant of Arabidopsis thaliana. J Photochem Photobiol B. 1996;34(1):87–94.
Murray MG, Thompson WF. Rapid isolation of high molecular weight plant DNA. Nucleic Acids Res. 1980;8(19):4321–5.
Arnqvist L, Persson M, Jonsson L, Dutta PC, Sitbon F. Overexpression of CYP710A1 and CYP710A4 in transgenic Arabidopsis plants increases the level of stigmasterol at the expense of sitosterol. Planta. 2008;227(2):309–17.
Suzuki Y, Kawazu T, Koyama H. RNA isolation from siliques, dry seeds, and other tissues of Arabidopsis thaliana. Biotechniques. 2004;37(4):542. 544.
Tzin V, Galili G. The Biosynthetic Pathways for Shikimate and Aromatic Amino Acids in Arabidopsis thaliana. Arabidopsis Book. 2010;8:e0132.
McFarlane H, Döring A, Persson S. The Cell Biology of Cellulose Synthesis. Annu Rev Plant Biol. 2014;65:69–94.
Smith AM, Stitt M. Coordination of carbon supply and plant growth. Plant Cell Environ. 2007;30(9):1126–49.
Sulpice R, Flis A, Ivakov AA, Apelt F, Krohn N, Encke B, Abel C, Feil R, Lunn JE, Stitt M. Arabidopsis coordinates the diurnal regulation of carbon allocation and growth across a wide range of photoperiods. Mol Plant. 2014;7(1):137–55.
Parsons PA. Energetic efficiency under stress underlies positive genetic correlations between longevity and other fitness traits in natural populations. Biogerontology. 2007;8(1):55–61.
Yun AJ, Lee PY, Doux JD, Conley BR. A general theory of evolution based on energy efficiency: its implications for diseases. Med Hypotheses. 2006;66(3):664–70.
Taiz L, Zeiger E. Plant Physiology. 4th ed. Sunderland: Sinauer Associates, Inc.; 2006.
Paredes M, Quiles MJ. The Effects of Cold Stress on Photosynthesis in Hibiscus Plants. PLoS ONE. 2015;10(9):e0137472.
He Y, Fu J, Yu C, Wang X, Jiang Q, Hong J, Lu K, Xue G, Yan C, James A, et al. Increasing cyclic electron flow is related to Na + sequestration into vacuoles for salt tolerance in soybean. J Exp Bot. 2015;66(21):6877–89.
Strand DD, Livingston AK, Satoh-Cruz M, Froehlich JE, Maurino VG, Kramer DM. Activation of cyclic electron flow by hydrogen peroxide in vivo. Proc Natl Acad Sci U S A. 2015;112(17):5539–44.
Kono M, Noguchi K, Terashima I. Roles of the cyclic electron flow around PSI (CEF-PSI) and O(2)-dependent alternative pathways in regulation of the photosynthetic electron flow in short-term fluctuating light in Arabidopsis thaliana. Plant Cell Physiol. 2014;55(5):990–1004.
Roach T, Krieger-Liszkay A. Regulation of photosynthetic electron transport and photoinhibition. Curr Protein Pept Sci. 2014;15(4):351–62.
Cheung CY, Ratcliffe RG, Sweetlove LJ. A Method of Accounting for Enzyme Costs in Flux Balance Analysis Reveals Alternative Pathways and Metabolite Stores in an Illuminated Arabidopsis Leaf. Plant Physiol. 2015;169(3):1671–82.
Miyake C. Alternative electron flows (water-water cycle and cyclic electron flow around PSI) in photosynthesis: molecular mechanisms and physiological functions. Plant Cell Physiol. 2010;51(12):1951–63.
Sweetlove LJ, Lytovchenko A, Morgan M, Nunes-Nesi A, Taylor NL, Baxter CJ, Eickmeier I, Fernie AR. Mitochondrial uncoupling protein is required for efficient photosynthesis. Proc Natl Acad Sci U S A. 2006;103(51):19587–92.
Poolman MG, Kundu S, Shaw R, Fell DA. Responses to light intensity in a genome-scale model of rice metabolism. Plant Physiol. 2013;162(2):1060–72.
Ma F, Jazmin LJ, Young JD, Allen DK. Isotopically nonstationary 13C flux analysis of changes in Arabidopsis thaliana leaf metabolism due to high light acclimation. Proc Natl Acad Sci U S A. 2014;111(47):16967–72.
Nogales J, Gudmundsson S, Knight EM, Palsson BØ, Thiele I. Detailing the optimality of photosynthesis in cyanobacteria through systems biology analysis. Proc Natl Acad Sci U S A. 2012;109(7):2678–83.
Voss I, Sunil B, Scheibe R, Raghavendra AS. Emerging concept for the role of photorespiration as an important part of abiotic stress response. Plant Biol (Stuttg). 2013;15(4):713–22.
Driouch H, Melzer G, Wittmann C. Integration of in vivo and in silico metabolic fluxes for improvement of recombinant protein production. Metab Eng. 2012;14(1):47–58.
Dersch LM, Beckers V, Wittmann C. Green pathways: Metabolic network analysis of plant systems. Metab Eng. 2016;34:1–24.
Costenoble R, Picotti P, Reiter L, Stallmach R, Heinemann M, Sauer U, Aebersold R. Comprehensive quantitative analysis of central carbon and amino-acid metabolism in Saccharomyces cerevisiae under multiple conditions by targeted proteomics. Mol Syst Biol. 2011;7:464.
Blank LM, Kuepfer L, Sauer U. Large-scale 13C-flux analysis reveals mechanistic principles of metabolic network robustness to null mutations in yeast. Genome Biol. 2005;6(6):R49.
McNerney MP, Watstein DM, Styczynski MP. Precision metabolic engineering: The design of responsive, selective, and controllable metabolic systems. Metab Eng. 2015;31:123–31.
Chastain CJ, Fries JP, Vogel JA, Randklev CL, Vossen AP, Dittmer SK, Watkins EE, Fiedler LJ, Wacker SA, Meinhover KC, et al. Pyruvate, orthophosphate dikinase in leaves and chloroplasts of C3 plants undergoes light-/dark-induced reversible phosphorylation. Plant Physiol. 2002;128(4):1368–78.
Steer BT. Control of diurnal variations in photosynthetic products: I Carbon metabolism. Plant Physiol. 1974;54(5):758–61.
Baerenfaller K, Massonnet C, Walsh S, Baginsky S, Buhlmann P, Hennig L, Hirsch-Hoffmann M, Howell KA, Kahlau S, Radziejwoski A, et al. Systems-based analysis of Arabidopsis leaf growth reveals adaptation to water deficit. Mol Syst Biol. 2012;8:606.
Noctor G, Foyer CH. Homeostasis of adenylate status during photosynthesis in a fluctuating environment. J Exp Bot. 2000;51 Spec No:347–56.
Alonso AP, Vigeolas H, Raymond P, Rolin D, Dieuaide-Noubhani M. A new substrate cycle in plants. Evidence for a high glucose-phosphate-to-glucose turnover from in vivo steady-state and pulse-labeling experiments with [13C]glucose and [14C]glucose. Plant Physiol. 2005;138(4):2220–32.
Lassen LM, Nielsen AZ, Ziersen B, Gnanasekaran T, Moller BL, Jensen PE. Redirecting photosynthetic electron flow into light-driven synthesis of alternative products including high-value bioactive natural compounds. ACS Synth Biol. 2014;3(1):1–12.
van Duuren JB, Puchalka J, Mars AE, Bucker R, Eggink G, Wittmann C, Dos Santos VA. Reconciling in vivo and in silico key biological parameters of Pseudomonas putida KT2440 during growth on glucose under carbon-limited condition. BMC Biotechnol. 2013;13:93.
Moisset P, Vaisman D, Cintolesi A, Urrutia J, Rapaport I, Andrews BA, Asenjo JA. Continuous modeling of metabolic networks with gene regulation in yeast and in vivo determination of rate parameters. Biotechnol Bioeng. 2012;109(9):2325–39.
The authors acknowledge Björn Juncker (University of Halle) for helpful discussion on the network topology.
VB created the metabolic models, designed the simulation experiments, performed the simulation studies, analyzed the results, drafted the figures and drafted the manuscript. LD, KL, GM contributed with database and literature screening to the validation of metabolic model topology. RF, TE, OB contributed by discussion on the manuscript. CW designed and supervised the work, designed the simulation experiments, analyzed the results and drafted the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.