Metabolic model of central carbon and energy metabolisms of growing Arabidopsis thaliana in relation to sucrose translocation
© The Author(s). 2016
Received: 22 December 2015
Accepted: 5 August 2016
Published: 28 December 2016
Sucrose translocation between plant tissues is crucial for growth, development and reproduction of plants. Systemic analysis of these metabolic and underlying regulatory processes allow a detailed understanding of carbon distribution within the plant and the formation of associated phenotypic traits. Sucrose translocation from ‘source’ tissues (e.g. mesophyll) to ‘sink’ tissues (e.g. root) is tightly bound to the proton gradient across the membranes. The plant sucrose transporters are grouped into efflux exporters (SWEET family) and proton-symport importers (SUC, STP families). To better understand regulation of sucrose export from source tissues and sucrose import into sink tissues, there is a need for a metabolic model that takes in account the tissue organisation of Arabidopsis thaliana with corresponding metabolic specificities of respective tissues in terms of sucrose and proton production/utilization. An ability of the model to operate under different light modes (‘light’ and ‘dark’) and correspondingly in different energy producing modes is particularly important in understanding regulatory modules.
Here, we describe a multi-compartmental model consisting of a mesophyll cell with plastid and mitochondrion, a phloem cell, as well as a root cell with mitochondrion. In this model, the phloem was considered as a non-growing transport compartment, the mesophyll compartment was considered as both autotrophic (growing on CO2 under light) and heterotrophic (growing on starch in darkness), and the root was always considered as heterotrophic tissue dependent on sucrose supply from the mesophyll compartment. In total, the model includes 413 balanced compounds interconnected by 400 transformers. The structured metabolic model accounts for central carbon metabolism, photosynthesis, photorespiration, carbohydrate metabolism, energy and redox metabolisms, proton metabolism, biomass growth, nutrients uptake, proton gradient generation and sucrose translocation between tissues. Biochemical processes in the model were associated with gene-products (742 ORFs). Flux Balance Analysis (FBA) of the model resulted in balanced carbon, nitrogen, proton, energy and redox states under both light and dark conditions. The main H+-fluxes were reconstructed and their directions matched with proton-dependent sucrose translocation from ‘source’ to ‘sink’ under any light condition.
The model quantified the translocation of sucrose between plant tissues in association with an integral balance of protons, which in turn is defined by operational modes of the energy metabolism.
KeywordsEnergy metabolism Multi-compartment metabolic model Central carbon metabolism Sucrose metabolism Sucrose transport Flux balance analysis Diurnal growth
The aim of this research was to build a multi-compartmental metabolic model of growing Arabidopsis thaliana. The integrated model should describe biomass growth both in light and in dark phases with corresponding formation and consumption of starch and sucrose. Furthermore, the structured metabolic model should take into account major pathways of primary metabolism such as sugar metabolism, central carbon metabolisms, photosynthesis, photorespiration, energy and redox metabolism, proton turnover, sucrose translocation from source to sink tissues and biomass growth. In future, we will use the model to recognize cause-effect relationships and describe regulatory processes in carbon metabolism and transport.
In growing plants, sucrose is the most widespread sugar used to supply both carbon and energy from ‘source’ tissues (e.g. autotrophic mesophyll) to ‘sink’ tissues (e.g. heterotrophic roots, growing shoots or reproductive organs) to build up a biomass . During photosynthesis in plastids of mesophyll cells, triose phosphates (GAP, DHAP) are synthesised and exported into the cytoplasm to support formation of sucrose and biomass. During growth in the light, starch is formed and accumulated in the plastids and becomes a part of the biomass [2, 3]. Starch is a repository of carbon which is later used during the dark phase as the primary carbon source for biomass formation  and fuelling of sucrose biosynthesis and its transport. The diurnal dynamics of starch accumulation is generally well documented in plants, and particularly in Arabidopsis thaliana  this process was even subjected to the analysis of regulatory patterns by means of dynamic mathematical modelling . Perturbation of these tightly regulated metabolic processes results in growth phenotypes of the plants. For example, disruption of the plant’s ability to invest carbon into the day-time storage of starch in the pgm mutant  results in higher cytosolic sucrose levels, higher respiration rates, retarded growth , low seed yield , and slow root growth at night .
In most sink tissues, sucrose is primarily used as a carbon source to support growth and build up biomass. However, sucrose can also serve as metabolite invested into storage compounds in root tissues. It can be converted to starch as in potato tubers or it can directly be stored in the vacuole as in sugar beet or sugar cane. Thus, the whole plant growth is tightly dependent on regulation of sucrose metabolism and transport. The cause-effect relationship of sucrose transport between tissues and phenotypic traits of plants is an important area of current plant research .
The analysis of metabolic networks with a rational approach is an efficient tool for engineering of plant systems [23, 24]. Model-based approaches are increasingly used also in plant biology to gain functional insights and even prediction of metabolic processes [25, 26]. In general, network modelling involves several stages: (i) collection of a priori knowledge and construction of the model by network reconstruction [27, 28]; (ii) a priori check of the proposed metabolic network model, such as topological analysis of stoichiometric matrix to eliminate structural gaps; (iii) addition of constraints and formulation of objective function to the model; (iv) network optimization using constraint-based analyses such as flux balance analysis (FBA); (v) a posteriori consistency check of the identified model by metabolic flux analysis (MFA); (vi) validation of the model with a newly generated experimental dataset [29, 30]. Based on this workflow the successful metabolic models of Arabidopsis thaliana of different scales were already suggested [23, 25, 28, 31–34] including multi-tissue genome scale metabolic models [32, 35, 36]. However, each of these existing metabolic models is dedicated to a specific biological phenomenon to be better understood, such as photosynthesis  or the Calvin-Benson cycle . On a more global scale, the metabolic costs for all amino acids and proteins in a given network were calculated based on the flux balance analysis of genome-scale metabolic network of Arabidopsis thaliana in light and dark conditions . Flux Balance Analysis combined with turnover measurements of 13C-labeled metabolites can provide deeper insights into the underlying metabolic processes. Stability of metabolic fluxes in central metabolism of Arabidopsis thaliana root cells was tested against environmental variation of oxygen using such approach [40, 41].
The stoichiometric models can in future be further developed into a dynamic model based on kinetic expressions of particular biochemical reactions and their integration [4, 37, 42, 43]. Kinetic modelling of plant metabolism was used to unravel local and global system features, such as flux and concentration control coefficients and regulation patterns  in relation to different external or internal states, stimuli and conditions . A kinetic model of sugar metabolism throughout tomato Solanum lycopersicum fruit development revealed importance of different enzymes on different development stages, as well as importance of sugar accumulation in vacuole together with organic acids to enable osmotic-driven vacuole expansion during the cell division . Also, different regulatory scenarios of starch turnover by the circadian clock through dynamic adjustment of starch turnover to changing environmental conditions were suggested based on mathematical modelling .
The multi-compartmental metabolic network of Arabidopsis thaliana was reconstructed and optimized in order to explain growth stoichiometry of the plant both in light and in dark conditions. Balances and turnover of energy (ATP/ADP) and redox (NAD(P)H/NAD(P)) metabolites as well as proton in different compartments were estimated. The model showed that in light conditions, the plastid ATP balance depended on the relationship between fluxes through photorespiration and photosynthesis including both cyclic and non-cyclic electron flow. The ATP balance in plastid depended on the ratio between these processes, and therefore can be either deficient, self-supported or producing a surplus. The excess of redox potential from the photosynthetic light system was translocated to the mitochondrion via the malate/oxaloacetate shuttle. The model showed that the mitochondrion consumed protons under both ‘light’ and ‘dark’ conditions and provided ATP to the cytoplasm. Matching the proton fluxes with proton-dependent translocation of sucrose between tissues from source to sink in light and dark conditions corresponded well to the known molecular mechanism of sucrose transport.
Results and discussion
In our model, we simplified the complex plant tissue organization to four principal compartments: (i) the super-compartment ‘plant’, which includes (ii) a growing autotrophic ‘mesophyll’, (iii) a heterotrophic growing ‘root’ and (iv) a non-growing ‘phloem’ compartment. Mesophyll, phloem and root were interconnected through inner space of the super-compartment ‘plant’, which played the functional role of ‘apoplast’ (Figs. 3 and 4). To model the exchange of solutes with the external environment, the ‘root’ was set to acquire water, protons, nutrients and transport them up to the ‘mesophyll’ via the ‘phloem’, while the ‘mesophyll’ was exchanging CO2, O2, excreting an excess of water, and providing sucrose in counter-flux to the ‘root’ via the ‘phloem’ (Fig. 4). The ‘mesophyll‘compartment included sub-compartments: plastid and mitochondrion, while the ‘root’ compartment contained only mitochondria. Additionally, the model included a set of reactions that imitated the vacuole, as a virtual sink (only accumulating) compartment. We assigned accumulation of nitrogen, orthophosphate, sulphur, and carbon in form of ash to the virtual vacuole. In addition, since sucrose reserves were also defined as a biomass constituent, therefore sucrose contributing to growth was also irreversibly stored in the virtual vacuole.
The network was formulated to focus only on the major relevant metabolic processes from the central carbon metabolism that contribute to the outlined phenomena of sucrose formation, translocation and degradation along with proton formation in metabolism and energy/redox reactions. The elemental and charge balances of each metabolic reaction in the model were verified based on elemental composition of compounds and their known charges.
The model was developed in an iterative way of (i) network reconstruction, formulation, correction, (ii) network setup analysis (to remove inconsistencies), (iii) topological analysis (to remove structural gaps and false-positive effects), and (iv) Flux Balance Analysis to optimize objective function and reach biological consistency. The outcomes of the predicted macroscopic fluxes were compared with generally accepted directions of macroscopic input/output fluxes of the plant metabolites, which are well documented for ‘light’ and ‘dark’ growth phases (Table 1). Additionally, directions of intracellular fluxes predicted by FBA were checked for the biological consistency and in case of inconsistency or biological irrelevance, the network setup was re-designed in order to eliminate these inconsistencies. Then cycle of the analyses was repeated until the model reflected the biologically relevant picture.
The main uptake path of proton from environment to the phloem cell in the model was in symport with nutrients such as orthophosphates, nitrates, or sulphates (Fig. 4). Nevertheless, there is also ATP dependent efflux of protons from cells through the activity of the plasma membrane H+-ATPase in a stoichiometric ratio 1:1 . Thus, the model takes in account simultaneous and independent influx/efflux of protons, but with overall net-flux of protons from the environment to the plant cells as determined experimentally from measured alkalinisation of the growth media for growing Arabidopsis roots. Under standard growth conditions in hydroponic setup , the pH value of the ½MS-medium was found to increase from 5.7 to 6.2 over a period of two weeks.
Special attention was paid to modelling of energy metabolism and photosynthetic light reactions. The photosynthesis light reactions were formulated in accordance with descriptions in AraCyc  and an existing model . There are two processes that can significantly influence the energetic efficiency of photosynthesis, namely (i) photorespiration (oxidative photosynthetic carbon cycle)  and (ii) cyclic electron flow through the photosynthesis light reactions [50, 51]. Photorespiration: Although in some cases photorespiration is considered as a ‘security valve’ to reduce the consequences of over-oxygenation in plastids, it is nevertheless considered as the essential part of the photosynthesis , which potentially reduces photosynthetic output by 20–30 % in C3 plants . Some reactions of photorespiration occur in peroxisome and mitochondrion, however to avoid complication of the model with introduction of an additional compartment (i.e. peroxisome) and introduction of additional pathways to mitochondria, the corresponding peroxisomal reactions were written within the cytoplasm of the mesophyll. In this way, photorespiration was connected to the cytoplasmic folate metabolism. Cyclic electron transport (CET): The cyclic electron transport was reported to occur mainly at high irradiances in combination with very low CO2 concentration , nevertheless the CET generally is reported as an essential element of the photosynthesis system [54–56]. Thus, the current version of the model includes photosynthesis with both non-cyclic and cyclic electron transports through photosynthetic light reactions and photorespiration. The proton motive force for ATP synthase in plastid was formed by the cytochrome b 6 /f complex and the proton/ATP stoichiometry of ATP synthase was fixed to 4.0 [57, 58].
Plastid metabolism was restricted to photosynthetic light reactions, Calvin-Benson cycle (CBC), pentose-phosphate pathway (PPP), sulphate and nitrite reduction pathways, synthesis/degradation of starch. Sulphate and nitrite reduction pathways were interconnected through conserved moiety of ferredoxin(red)/ ferredoxin(ox) with photosynthetic light reactions, as well as to NADPH/NADP+ conserved moieties. Under light conditions, reduced ferredoxin was mainly consumed by ferredoxin-NADP+ oxidoreductase (FNR), which provided NADPH for the Calvin-Benson cycle, whereas the oxidative part of PPP was inactive. However, under dark conditions, when photosynthesis light reactions were inactive, the oxidative part of PPP became active and, in the model, it was considered as the only provider of NADPH in plastids during darkness. Formed NADPH was used by FNR to reduce ferredoxin and therefore to feed sulphate and nitrite reduction during dark phase.
We did not consider other central carbon metabolic pathways in plastid (fatty acids synthesis, amino acids and glutathione biosynthesis), due to increasing complexity of the model and since these pathways did not contribute to a better reflection of the molecular relationships of the phenomenon under study. Thus, the elaborated pathways considered in the plastid were just enough to model photosynthesis, energy and redox balance, proton balance, CO2 fixation, triose phosphate synthesis as well as starch synthesis and degradation.
During the light-growth-phase, starch was synthesized and stored as biomass constituent. During the dark-growth-phase it became the only carbon and energy source. The model considered its consumption as a virtually external metabolite.
Folate (or C1) metabolism
The tetrahydrofolate transformations are important part of C1-metabolism . This pathway connected the transfer of C1-group with amino acid metabolism, nitrogen metabolism , and with photorespiration when it was active.
Nitrogen and sulphate metabolism
In the mesophyll, nitrate was reduced to nitrite in the cytoplasm and then transported to the plastid for further reduction to glutamate/glutamine, which were returned back to the cytoplasm, and further used in transamination reactions for synthesis of others amino acids. Sulphate was transported to the plastid, where it was reduced to H2S, which was returned back to cytoplasm and used further to synthesize cysteine . Nitrate and sulphate reduction are metabolic processes that demand reduced equivalents and overall utilize free protons . Therefore, these reactions were additional redox-loads for the plastidic redox-balance. In roots, the nitrogen and sulphate metabolism was modelled in the cytoplasm and were fed by NADPH produced by the oxidative part of cytoplasmic PPP. This was a compromise to avoid modelling amyloplasts in root [61, 62].
To ensure functional metabolite exchange of metabolites and redistribution of energy and redox load between organelles (plastid, mitochondrion) and cytoplasm so-called metabolite translocators and shuttles were included in the model. It is well documented that the energy-, redox-, sulphur-, nitrogen- and carbon- metabolisms are interdependent via such exchange processes [12, 61, 63]. We integrated the following exchange processes in the model: ATP/ADP translocation among plastid/cytoplasm/mitochondrion , redox exchange via malate/oxaloacetate shuttle (dicarboxylate translocators) between plastid/cytoplasm/mitochondrion [64, 65], carbon (DHAP, GAP) exchange (triose phosphate/phosphate translocator) between plastid/cytoplasm [61, 63, 66], amino acids and ammonia exchange between plastid/cytoplasm [61, 63, 65], organic acid exchange between cytoplasm/mitochondrion [1, 65].
The final version of the model (Additional file 1: Figure S1) was composed of 400 transformers, among which were 229 metabolic reactions from 41 different metabolic pathways, 155 transporters and 16 polymerization steps. The transformers connected 423 compounds, among which 413 were balanced compounds. In the model, external metabolites such as gases, nutrients, biomass were unbalanced compounds, because they were considered as being an infinite source or sink. All biochemical reactions were associated with 742 genes, whose products performed the corresponding biochemical reactions.
Metrics and topological indicators of the stoichiometric model of Arabidopsis thaliana
Degree of freedom (df)
There are 28 conserved moieties across all compartments of the model, among them there are only a biologically determined moieties:
- NAD(P)H/NAD(P) [all compartments]
- CoQH2/CoQ [mitochondria]
- THF/5-Formyl-THF/5,10-Methylene-THF/5-Methyl-THF [tetrahydrofolate : cytoplasms]
- CoA/Acetyl-CoA [CoA : cytoplasms]
- CoA/Acetyl-CoA/Succinyl-CoA [CoA : mitochondria]
- CytC(red)/ CytC(ox) [cytochrome c : mitochondria]
- Fd(red)/Fd(ox) [ferredoxin : plastid]
- PQH2/PQ [plastoquinone : plastid]
- PC(red)/PC(ox) [plastocyanine : plastid]
- ATP/ADP [mitochondria]
- ADP-glc/ATP/ADP/AMP/APS [plastid]
15 solution spaces
Parallel routes & cycles
Analysis of conserved moieties in all compartments revealed only those with biological relevance (Table 2). Interestingly, conserved moiety ATP/ADP/AMP was not in the list of cytoplasmic conserved moieties unlike in the plastid, since cytoplasmic AMP is de novo synthesized and interconnected to the ATP/ADP pool via adenylate kinase, and, at the same time in the model, AMP was a monomer for RNA/DNA polymerization. Thus, ATP/ADP/AMP was an open moiety in the cytoplasm.
Determination of the nullspace of a stoichiometric matrix of a model is an important topological criterion to assess the feasible space containing all possible solutions of the linear equations S × v = 0 without taking into account thermodynamic irreversibility. Thus, the nullspace analysis is used to reveal non-functional regions of the network by finding solution spaces of all sub-networks that are able to operate at steady state. In the model there were 15 overlapping sub-networks that all together cover all the compounds and reactions in the suggested network.
Flux Balance Analysis (FBA)
Table of flux constraints used for Flux Balance Analysis (FBA) of the stoichiometric model of Arabidopsis thaliana
Interpretation of the constraints
T.Biomass.ext < = 5000
SK = 0
0.25*RPC_plastide - RPC2_plastide = 0
T.starch.ext = 0
T.CO2.ext = 0
T.CO2.ext_rev > = 0
FQR = 0
T.ADP.plastid - T.ATP.plastid = 0
(0 … 0.5)*FNR - FQR = 0
RPC2_plastide = 0
GLYK = 0
Maximization of biomass formation under light assumes that:
• A plant consumes CO2 as the carbon source (T.CO2.ext = 0, T.CO2.ext_rev > = 0) there are two constraints because CO2 transport is reversible
• All formed starch becomes a part of the biomass (T.Starch.ext = 0)
• There is no starch degradation under light and therefore starch kinase is inactive (SK = 0)
• Photorespiration flux was fixed at 20 % of flux through RuBisCo (0.25*RPC_plastide - RPC2_plastide = 0) at all tested conditions
Additionally, the impact of different degree of cyclic electron flow through photosynthesis light reactions (particularly through ferredoxin-plastoquinone reductase; FQR) was estimated by means of add-on constraints:
• Either … cyclic electron flow is inactive (FQR = 0)
• Or … flux through ATP/ADP translocator between plastid and cytoplasm is inactive (T.ADP.plastid – T.ATP.plastid = 0), thus plastid’s ATP balance becomes self-sufficient
• Or … flux through FQR varies relative to the non-cyclic electron flow through ferredoxin-NADP+-oxidoreductase (FNR) ((0 … 0.5)*FNR – FQR = 0)
• Or … flux through photorespiration pathway is zero (RPC2_plastide = 0; GLYK = 0), but cyclic electron flow through FQR is subjected to optimization
Resulted growth stoichiometry for ‘light’ conditions (the stoichiometric coefficients [mol i/mol X] are normalized per biomass):
[PR* = 0.25; FQR = 0]:
[PR = 0.25; FQR/FNR = 0.37]:
[PR = 0.25; FQR/FNR = 0.5]:
[PR = 0; FQR/FNR = 0.1]:
3.26e7 × CO 2 + 2.12e7 × H 2 O + 1.00e6 × HPO 4 2 − + 1.99e4 × SO 4 2 − + 1.52e6 × H + + 5.04e6 × NO 3 − + 5.26e8 × hv = > 3.98e7 × O 2 + Biomass_plant
3.26e7 × CO 2 + 2.12e7 × H 2 O + 1.00e6 × HPO 4 2 − + 1.99e4 × SO 4 2 − + 1.52e6 × H + + 5.04e6 × NO 3 − + 5.58e8 × hv = > 3.98e7 × O 2 + Biomass_plant
3.26e7 × CO 2 + 2.12e7 × H 2 O + 1.00e6 × HPO 4 2 − + 1.99e4 × SO 4 2 − + 1.08e6 × H + + 5.04e6 × NO 3 − + 5.92e8 × hv = > 3.98e7 × O 2 + Biomass_plant
3.26e7 × CO 2 + 2.12e7 × H 2 O + 1.00e6 × HPO 4 2 − + 1.99e4 × SO 4 2 − + 5.84e6 × H + + 5.04e6 NO 3 − × + 3.90e8 × hv = > 3.98e7 × O 2 + Biomass_plant
T.Biomass.ext < = 5000
PGM3_plastid = 0
RPC_plastide = 0
RPC2_plastide = 0
GLYK = 0
GDC = 0
T.starch.ext > = 0
T.hv.ext = 0
Maximization of the biomass formation in darkness assumes that:
• There is no light (T.hv.ext = 0)
• Therefore RuBisCo is inactive (RPC_plastide = 0)
• Correspondingly photorespiration is also inactive (RPC2_plastide = 0, GLYK = 0, GDC = 0)
• Consequently there is no synthesis of new starch, therefore phosphoglucomutase is inactive (PGM3_plastid = 0)
• Starch is the carbon source for biomass formation (T.Starch.ext > = 0), which previously has been deposited in course of the light phase
Resulted growth stoichiometry for ‘dark’ conditions (the stoichiometric coefficients [mol i/mol X] are normalized per biomass):
8.48e6 × O 2 + 1.00e6 HPO 4 2 − × + 1.99e4 × SO 4 2 − + 3.84e6 × H ++ 5.04e6 × NO 3 − + 1.34e4 × Starch = > 1.55e7 × CO 2 + 1.92e7 × H 2 O + Biomass_plant
Modelling of the metabolic activity of the mesophyll was the most complex task due to light-dependent activity of energy-, redox- and carbon-metabolism in plastid, in terms of CO2 and O2 consumption/production, starch formation/degradation, triose phosphate export/import, or under conditions of photorespiration on/off. The proton balance in the mesophyll was especially influence by energy and redox metabolism particularly water photolysis and oxygen reduction (Figs. 4 and 5). In contrast, the root metabolism was almost invariant under light and dark growth condition, being the constant sink for sucrose.
Chosen stoichiometry of the anabolic reactions leading to the biomass formation predicted an elemental composition of the Arabidopsis thaliana biomass as CH1.592O0.834N0.144P0.033 (MW x = 29.88 [g dw /Cmol]), which only slightly differed from averaged elemental biomass compositions of microbial species (CH1.596O0.396N0.216P0.017; MW x = 24.59 [g dw /Cmol]) [69–72]. We cannot judge about significance of the observed differences, since there is no confident information published on elemental composition of Arabidopsis thaliana.
ATP and NAD(P)H balances in mesophyll
While interpreting of the FBA results, main attention was placed on the metabolic activity of the mesophyll, since in the model the mesophyll constituted 85 % (w/w) of the mass fraction within the plant biomass.
Photosynthetic light reactions (PLR) in the plastid were modelled in details in order to provide maximal plasticity in ATP and NADPH allocations. It is well documented, that fixation of CO2 to yield triose phosphate (GAP or DHAP) in the Calvin-Benson cycle requires a maximal theoretical ATP/NADPH ratio of 1.5 [1, 63]. However, non-cyclic photosynthetic electron transport (in the way as it is written in AraCyc ) provides an ATP/NADPH ratio of 1.0 and experimentally measured values range up to 1.3 . Thus, the demand for ATP in plastid can exceed the level of ATP synthesis provided by non-cyclic electron transport in photosynthetic light reactions. To overcome this limitation, the C3 plants use several pathways to adjust the ATP/NADPH ratio in the plastid stroma and to finally ensure CO2-fixation in the Calvin-Benson cycle: (i) increased cyclic electron transfer in photosynthetic light reactions to increase ATP yield without concomitant NADPH yield ; (ii) increased ATP consumption through photorespiration; (iii) importing/exporting ATP from/to cytoplasm via ATP/ADP translocator [61, 63]; (iv) NADPH consumption in reductive biosynthetic pathways in the plastid stroma such as nitrite- and sulphate-reduction ; (v) shuttling of reduced equivalents from plastid to mitochondrion via the malate/oxaloacetate shuttle  (Fig. 5); (vi) shuttling of ATP and reduced equivalents to the cytoplasm through DHAP/GAP shuttle  (which was not included into the model, for the reasoning see above).
To explore the contributions of all pathways to plastid ATP/NADPH balance, the FBA was correspondingly constrained in a series of independent runs (Table 3). Growth in the light meant that CO2 was the carbon source, photosynthesis (including both cyclic and non-cyclic electron flow) and also photorespiration were active, NADPH was generated by photosynthesis, ATP was generated by the plastidic ATP-synthase, but dependent on the actual ATP/NADPH balance it was additionally imported/exported from/to cytoplasm, triosephosphates were formed in the Calvin-Benson cycle and used both for export to cytoplasm and starch synthesis. Growth in the darkness meant that carbon source was starch and in plastid the photosystem I and II, RuBisCO, and therefore photorespiration were inactive, ATP was completely imported from the cytosol and NADPH was generated by partially active Calvin-Benson cycle and PPP fed by triosephosphates obtained from starch degradation.
Redox-potentials are required for nitrogen- and sulphur-reduction pathways and are provided in form of reduced ferredoxin, either by PLR under the light or by PPP in the dark. In the model under light conditions, activity of these pathways withdrew electron flux and correspondingly reduced overall NADPH yield of photosynthesis. NADPH formed by ferredoxin-NADP-reductase was mainly utilized by the Calvin-Benson cycle in order to drive CO2 fixation and to form triose phosphates. The excess of NADPH was diverted from plastid to mitochondrion via malate/oxaloacetate shuttle . In mitochondria, the exported reduced equivalents were used for ATP formation via oxidative phosphorylation (Fig. 5). According to FBA, in the dark, the plastid consumed imported GAP from the cytoplasm to provide carbon flux through the CBC and further to oxidative part of PPP in order to generate NADPH required for reductive processes of the CBC itself, as well as for N- and S-reduction. A by-product of the oxidative part of PPP was CO2, which was partially utilized in CBC itself, but was also excreted from the plastid.
The ATP synthesis in the plastid by ATP synthase (with stoichiometry H in + /ATP = 4.0 [57, 58]), was driven by the H in + − motive force between thylakoid lumen and the plastid stroma (Fig. 5). Therefore, addition of cyclic electron transfer by ferredoxin-plastoquinone reductase (FQR) to the photosynthetic light reactions contribute to an increase in ATP yield without corresponding NADPH yield [55, 56]. Additionally, we also investigated potential impact of H in + − leakage from lumen back to stroma, which increases H in + /ATP stoichiometry. As it appears, within the framework of this model, it was impossible to provide reliable constraints to estimate the contribution of this process, therefore this process was eliminated. Overall, the FBA of the current form of the model predicted that ATP can be imported or exported from/to the cytosol to fulfil the ATP balance in photosynthetically active plastids. The overall balance of ATP in the plastid was regulated by levering photorespiration and cyclic electron flow through the photosynthesis light reactions.
Comparison of metabolite turnover predicted by FBA with known values for Arabidopsis thaliana
Ratio of precursor turnover [mol i/mol j]
FQR f) = 0 PR e) = 0.25
FQR / FNR f) = 0.37 PR e) = 0.25
FQR / FNR f) = 0.5 PR e) = 0.25
FQR / FNR f) = 0.1 PR e) = 0
Photon/NO 3 −
CO2/NO 3 −
2.5 – 5
1.3 – 1.5
up to 3
H in + /ATP
Metabolite turnover [mol i/mol X]
ATP (× 1011)
NADPH (× 1011)
H+-turnover (× 1012)
The efficiency of the photosynthetic reactions could be evaluated based on the precursor ratios (Table 4). For example, the achieved hv/CO2 ratio in the model was around 16, which was close to the published hv/CO2 ratio of 9.55 for Arabidopsis thaliana . Also other metabolic indicators were predicted to be within the range of published values (Table 4), suggesting that despite the limitations outlined above, in general the stoichiometry of the plastid metabolism was correctly presented.
The mitochondrion in our model contained only two pathways, the TCA cycle and the oxidative phosphorylation. All intermediates of TCA cycle were exchangeable with the cytoplasm, except succinyl-CoA, which was a part of the mitochondrial conserved moiety CoA/acetyl-CoA/succinyl-CoA. Therefore, in the ‘light’ scenario the reactions of NAD-dependent isocitrate dehydrogenase (EC 220.127.116.11), 2-oxoglutarate dehydrogenase and succinate-CoA ligase (EC 18.104.22.168) were inactive, since 2-oxoglutarate was exchanged into the cytoplasm and through anaplerotic reactions was converted into succinate, which returned into the TCA cycle bypassing the reactions of NAD-dependent isocitrate dehydrogenase, 2-oxoglutarate dehydrogenase, and succinate-CoA ligase.
However, in the ‘dark’ scenario, these reactions were active due to differentially loaded glycolysis and consequently different input with primary metabolites into the anaplerotic reactions. The TCA cycle produced its own NADH, but also accepted cumulative redox-potential translocated into the mitochondrion through the malate/oxaloacetate shuttle from the plastid and cytoplasm (Fig. 5). Thus, in the light, the mitochondrial MDH produced 99 % of the NADH (9 % due to own TCA flux and 91 % due to import of L-malate from the cytoplasm), pyruvate dehydrogenase produced 1 % of mitochondrial NADH, and NAD-dependent isocitrate dehydrogenase was predicted to be inactive. CO2 production by mitochondria was very low in the light conditions due to inactivity of NAD-dependent isocitrate dehydrogenase and 2-oxoglutarate dehydrogenase, although the ATP production was quite high. Consequently, the predicted ATP/CO2 ratio was 65.9. In the dark, the NADH balance was complemented by activities of all reactions of the TCA cycle and correspondingly CO2 production was high leading to an ATP/CO2 ratio of 6.
In the model, all formed NADH in mitochondrion was consumed by complex I of the oxidative phosphorylation. There were no assigned ATP-consuming reactions in the mitochondrion, therefore all formed ATP from complex V of the oxidative phosphorylation was exported to the cytoplasm. The efficiency of oxidative phosphorylation was assessed by the amount of ATP formed per amount of oxygen consumed (P/O ratio). In reality, the P/O ratio is a variable parameter due to substrate availability and proton leakage. The theoretical value of P/O is 5.0, but practically measured values are generally greater than 2.5 . The FBA analysis suggested that the P/O ratio of the model was greater than 4.6 in both mesophyll and root mitochondria and in both ‘light’ and ‘dark’ conditions (Table 4). Thus, we assume that the balances of redox and energy reactions accepted in the sub-compartment ‘mitochondria’ of the model depict biological reality.
Major ATP producing/consuming processes in cytoplasm of the mesophyll quantified by FBA
Major ATP …
‘Light’ FQR / FNR = 0.37 photoresp. = 0.25
83.2 % import from mitochondrion 16.8 % by glycolysis and other processes
84.3 % import from mitochondrion 15.7 % by glycolysis and other processes
74.4 % by cata-/anabolic reactions 25.6 % by H+-ATPase
74.1 % by cata-/anabolic reactions 7.7 % by H+-ATPase 18.2 % export to plastid
Relative ATP turnover [mol ATP/mol X]
The non-growth associated ATP expenses of the plasma membrane H+-ATPase can be considered as maintenance. According to the FBA, the plasma membrane H+-ATPase consumed 25.6 % of the cytoplasmic ATP in light conditions (in case of FQR/FNR = 0.37, photorespiration = 0.25) and only 7.7 % of the cytoplasmic ATP in dark conditions (Table 5).
The NADH balance in the cytoplasm was comprised by major producing reactions (anplerotic reactions, glycolysis, C1-metabolism, pectin biosynthesis, nucleotides and amino acids biosynthesis, photorespiration) and consuming (cytoplasmic nitrate reduction, C1-metabolism and photorespiration in light conditions) reactions. The FBA of the model predicted active export of reduced equivalents from the plastid to the mitochondrion via the malate/oxaloacetate shuttle in form of malate and oxaloacetate during light conditions (Fig. 5). However, in light conditions photorespiration (particularly hydroxypyruvate reductase) became an additional NADH sink. In the dark, there was no export of reduced equivalents from the plastid, but additionally relatively higher glycolytic flux provided a surplus of NADH, which was translocated into the mitochondrion via the malate/oxaloacetate shuttle to be converted into ATP (Fig. 5). Thus, in the light, the mitochondrion was the sink for excess reduced equivalents formed in the plastid during photosynthesis, and photorespiration was additional sink for cytoplasmatic NADH. In the dark, the mitochondrion was the only sink for excess (up to 35 %) of NADH formed in the cytosol by relatively high glycolytic flux. At the same time, the mitochondrion was the major ATP supplier for cytoplasmic reactions both in light and dark growth conditions.
Proton balance & sucrose translocation
Numbers of H+-exchanging processes (reactions and transport) accounted for in the model
Exchange with env.
In the model, protons entered Arabidopsis thaliana from the environment in symport with nutrients into the root (e.g. nitrate, sulphate, phosphate; Fig. 4 and Additional file 1: Figure S1). However, the amount of protons entered in that way exceeds the amount protons going to transport the nutrients, therefore the excess of protons was excreted through H+-ATPase back to environment (Fig. 7) resulting in overall consumption of protons by the plant and accompanied alkalinisation of the growth medium. In the model, the efflux of protons from compartments was carried out by the plasma membrane H+-ATPase (Fig. 7). The turnover of protons predicted by the FBA in each compartment was normalized relative to the turnover of protons in the root cytoplasm (Fig. 7). The metabolic activity and the proton turnover of the root compartment in the model remained almost invariant under light and dark conditions and therefore this compartment was selected as a reference point for data normalization.
In both growing tissues (mesophyll and root), the major consumer of cytoplasmic protons was the mitochondrion (Fig. 7). For example, mesophyll mitochondria consumed 50 % of cytoplasmic protons under light growth conditions (in case of FQR/FNR = 0.37, photorespiration = 0.25). In mitochondria, major proton utilizing processes were oxygen reduction in oxidative phosphorylation. Thus, in order to maintain trans-membrane proton gradients, mitochondria constantly needed additional proton input, which came from the cytoplasm.
Also many metabolic processes resulted in proton production/utilization and therefore additionally contributed to the proton balance in the cytoplasm. Particularly in the light (in case of FQR/FNR = 0.37, photorespiration = 0.25), the plastid provided GAP, which entered glycolysis downstream of sugar phosphorylation reactions, which are producing substantial amounts of protons (Fig. 7). In the plastid, under light conditions, ATP synthase and photosynthesis light reactions produced a surplus of protons due to water photolysis, which was completely utilized by the CBC and nitrogen/sulphur reduction. Thus, turnover of proton in the photosynthetically active plastid was very high relative to the other compartments (Fig. 7), and the plastdal proton balance is self-sufficient in this conditions (FQR/FNR = 0.37, photorespiration = 0.25).
Under dark conditions, starch was decomposed to glucose, which involved phosphorylation by glucokinase (GLK; EC 22.214.171.124) and a release of protons contributing to approximately 17 % of the cytoplasmatic proton pool. This proton excess was mainly utilized by the fully functional TCA cycle in mitochondria. Correspondingly, the mesophyll in dark conditions pumped less protons to the apoplast (Fig. 7). The overall proton consumption by mitochondria was higher in dark than in light growth conditions (Fig. 7, Table 3).
FBA of proton fluxes allowed quantifying the net-fluxes of protons from phloem to mesophyll and root. Further, taking in account pH differences between compartments, the protons enter the root cytoplasm from the phloem along their concentration gradient allowing symport of solutes including sucrose via SUC and STP transporters. The sucrose flux to the root (as the sink tissue) was in accordance with proton symport mechanisms of SUCs and STPs and the proton balance of the tissues involved (Figs. 3, 4 and 7). Furthermore, the directions of proton fluxes between mesophyll and phloem did not disturb sucrose efflux from the mesophyll to the phloem due to independence of the SWEET transporters from the proton motive force. Under both scenarios (‘light’ and ‘dark’) the proton fluxes estimated by the FBA (Fig. 7) matched specific requirements by the molecular mechanisms of sucrose transporters (SWEET and SUC, STP) in the different tissues (Figs. 3, 4 and 7).
In this work we present a multi-compartmental metabolic model of growing Arabidopsis thaliana. The flux balance analysis (FBA) of the model quantified sugar metabolism, central carbon metabolism, photosynthesis, energy and redox metabolism, proton turnover, sucrose translocation from mesophyll to root and biomass growth under both dark- and light-growth conditions with corresponding growth either on CO2 (in the light) or on starch (in darkness). The model showed that in light conditions, interplay between photorespiration and photosynthesis including both cyclic and non-cyclic electron flow defined the ATP balance of the plastid. The plastid was found to be either deficient, self-supported or producing a surplus of ATP. The excess of redox potential from the photosynthetic light system was translocated to the mitochondrion via the malate/oxaloacetate shuttle. Thus, photosynthetically active plastids could achieve the ATP/NADPH ratio required for CO2 fixation. Also, the FBA predicted that the mitochondria were the main ATP provider for cytoplasmic processes together with glycolysis under both light- and dark-growth conditions. At the same time, the mitochondria were the main sink for reduced equivalents translocated both from plastid and from cytoplasm as well as for the protons.
The model described all main metabolic processes participating in the proton metabolism. The translocation of sucrose among plant tissues was associated with an integral balance of protons, which in turn was partially defined by operational modes of the energy metabolism (photosynthesis, respiration). The proton fluxes predicted by FBA generally corresponded to the molecular mechanisms and functional peculiarity of the sucrose transporters SWEET and SUC/STP proton-symporters and net-flux of sucrose from source to sink tissue. Thus, our multi-compartmental model adequately described the growth stoichiometry and sucrose translocation processes of Arabidopsis thaliana in light and dark conditions. The FBA was used to optimize the network structure and to prepare it for further analysis of the distribution of the intracellular fluxes based on experimental measurements by Metabolic Flux Analysis (MFA).
The primary information on the metabolic reactions regarding stoichiometry, direction, metabolite elemental formula, charge state, and associated genes was collected from AraCyc  using the Pathway Tool v19.0  and additionally were verified with TAIR [77, 78], KEGG  and ChEBI  databases.
The model reconstruction was performed with Insilico Discovery™ (Insilico Biotechnology AG, Stuttgart, Germany). The structured metabolic network of Arabidopsis thaliana was reconstructed based on the decision to take in account only those metabolic processes, which contributed to better understanding of sugar metabolism, central carbon metabolism, energy metabolism, proton turnover, biomass growth and sucrose translocation among tissues in Arabidopsis thaliana. Thus, metabolic reactions from the following pathways were used in the reconstructed network: 2-oxoglutorate decarboxylation to succinyl-CoA; adenosine nucleotides de novo biosynthesis; amino acids biosynthesis; aspartate degradation; Calvin-Beson cycle; chorismate biosynthesis; fatty acids biosynthesis; folate metabolism; glutamine biosynthesis; glycolysis; homoserine biosynthesis; inosine-5’-phosphate biosynthesis; maintenance; malate/oxaloacetate shuttle; nitrate reduction; ornithine biosynthesis; oxidative phosphorylation; PRPP biosynthesis; pentose-phosphate pathway; phosphorus metabolism; photosynthesis light reactions; photorespiration; purine nucleotides de novo biosynthesis; pyrimidine ribonucleotides interconversions; pyruvate decarboxylation to acetyl-CoA; pyruvate fermentation; starch biosynthesis; starch degradation; sucrose biosynthesis; sucrose degradation; sulphate reduction; TCA cycle; UDP-glucose biosynthesis; UDP-glucoronate biosynthesis; uridine-5’-phospahte biosynthesis. Many biosynthetic reactions concatenated in linear pathways (e.g. nucleotide or amino acid biosynthesis) were lumped in order to reduce model complexity. All reactions in the model were manually transferred from the database, and manually curated.
The biomass was assumed to be composed of the following polymers: proteins, lipids, RNA, carbohydrates (i.e. starch, sucrose, cellulose, pectin) and ash [33, 40]. The amino acids composition of a protein was accepted from literature . The lipid composition was assumed to be similar to yeast one and considered as a polymer of 9.6 % hexadecanoate, 20.6 % cis-octadecenoate, 34.7 % cis-hexadecenoate, 2.4 % octadecanoate, 0.006 % Linoleate and 32.4 % glycerol. RNA composition was derived by reversed sequence analysis of protein amino acid sequences in FASTA format available from TAIR: 32 % A, 26 % U, 9 % C and 33 % G. Based on these data, the biomass macromolecular composition was considered to be (w/w) 20 % cellulose, 20 % pectin, 10 % starch, 12 % sucrose, 25 % proteins, 9 % lipids, 1 % RNA and 3 % ash. The mass ratio between shoot and root was accepted as 85/15 (w/w) based on own measurements of fresh weight of 21 days old plants grown under long-day (16 h light) conditions in hydroponic culture.
In accordance with concept of sucrose emitting tissues and sucrose demanding tissues as well as under consideration of the molecular mechanism of sucrose translocation mechanisms (Fig. 3), under different energy producing modes (‘light’ and ‘dark’) we generalized the following compartments: a super-compartment ‘plant’ which included (i) the autotrophic sub-compartment ‘mesophyll’ with sub-compartments ‘plastid’ and ‘mitochondrion’, (ii) the heterotrophic sub-compartment ‘root’ with only one sub-compartment ‘mitochondrion’ and (iii) non-growing transport compartment ‘phloem’; whereas the inner space of the super-compartment ‘plant’ served as the extracellular compartment ‘apoplast’. Sub-compartment ‘plastid’ was required to separate photosynthetic light reactions, sulphate and nitrite reduction, Calvin-Benson cycle with associated triose phosphates, hexose phosphates, NADPH/NADP moiety, ATP synthesis and starch synthesis (Additional file 1: Figure S1). Sub-compartment ‘mitochondrion’ was required to separate TCA cycle and associated with its activity NADH/NAD moiety and ATP synthesis (Additional file 1: Figure S1). Compartmentalization was required in order to separate pools of similar metabolites between compartments with different metabolic specificities and to orchestrate their exchange through transport steps . Correspondingly, ‘mesophyll’ and ‘root’ were considered growing compartments, whereas ‘phloem’ was assumed a non-growing connecting compartment. The exchange processes between tissues were routed through the apoplast.
In the model, each biochemical reaction (classified with E.C. number) was referred to corresponding KEGG reaction (KeggID). Additionally, reactions were annotated with accession numbers of genes whose products perform this biochemical process (GeneID). In case of lumped reactions, all corresponding GeneIDs were associated with lumped reaction. This annotation was later used to match the model with mass-spectrometry based proteomic data  and to validate the presence of reactions in the model.
The final version of the model was exported in SBML and MATLAB format (Supplemental data) as well as the stoichiometric matrix of the model and high-resolution printout (JPG) of the model with all marked reactions. Corresponding databases of the compounds/metabolites, transformers/reactions and genes used in the model are provided in in ASCII format (Supplemental data).
Model assumptions and biologically based constraints
Sucrose was considered as the only carbohydrate that is exchanged among tissues.
The anatomy of the whole plant was reduced down only to three tissues: autotrophic mesophyll, heterotrophic root, interconnected by transport phloem. The model was restricted to the mesophyll as always being the ‘source’ tissue of sucrose both in light conditions (during photosynthesis) and in dark conditions (non-photosynthetic condition, carbon derived from starch). The root was always defined as the ‘sink’ tissue for sucrose. The root took up sucrose and used it for biomass formation, but did not store it with a purpose to emit it back to phloem. Thus, we have consciously constrained the exchange of sucrose to a unidirectional net-flux from autotrophic mesophyll to heterotrophic root (Figs. 3 and 4).
Sucrose was considered to be also a part of the formed biomass. Nevertheless, we did not consider the possibility that deposited sucrose can be returned to sucrose pool into the phloem for its exchange. Thus, sucrose, which was directed to biomass, was not considered as available for growth or exchange through the free sucrose pool.
Starch was considered as a metabolite that can be formed and deposited during light growth phase and used as the carbon and energy source during dark growth phase. We assumed that starch was accumulated only in autotrophic tissues (mesophyll), but not in heterotrophic tissues (root).
Starch was formed only during the light growth phase, and stored as a part of biomass. It was not degraded even partially during light growth phase.
Starch was considered as the only carbon source for growth of the plant during dark conditions.
In this version of the model, we did not consider exchange of amino acids and organic acids (glutamate, acetic, citric acids) between tissues.
Fatty acid biosynthesis was considered to be cytoplasmic to reduce complexity of redox balance in plastid.
ATP/ADP was made exchangeable between compartments (Fig. 4) through the adenine nucleotide translocator.
NAD(P)H/NAD(P) could not be exchanged between compartments. However, the reduced equivalent could be transferred between compartments through malate/oxaloacetate shuttle.
Phloem played only a role in transport processes and deliver nutrients and water from root to mesophyll, and sucrose from mesophyll to root as unidirectional net-fluxes (Fig. 4).
Nutrient uptake (nitrate, orthophosphate and sulphate) followed a proton-symport mechanism.
The only exchangeable metabolites from central carbon metabolism between plastid and cytoplasm in the mesophyll were D-glyceraldehyde 3-phosphate (GAP) and maltose. Other potentially exchangeable metabolites (i.e. other triose phosphates, hexose phosphates) were not considered in the model to reduce the complexity.
Gas exchange between plant and environment was considered as being passive.
Before proceeding to advanced network analysis, the detailed consistency check of the network stoichiometry through elemental and charge balance analysis was performed with Insilico Discovery™ (Insilico Biotechnology AG, Stuttgart, Germany). In addition to known reaction stoichiometry, this analysis relied on known elemental composition and charge (at specific pH) of metabolites. Consistency check of elemental and charge balances allowed identifying and revising/eliminating: (i) compounds for which elemental composition and charge were unknown and (ii) inconsistent reactions, polymerizations and transport steps with unbalanced elements or charges. A Check for dead-ends and unused transformers allowed identifying and eliminating transformers and network regions which were not capable of operating at steady state. A check for irreversible transformers allowed visualizing and approving those transformers, which were irreversible from thermodynamic point of view. Irreversibility of a transformer was important constraint for the Flux Balance Analysis.
Topological analysis of the network
Topological analysis of the stoichiometric matrix S of the network is very powerful tool in optimization of the metabolic network structure, because it identifies and illustrates how substances flows within a network. This analysis also allowed identification of false or/and unused network structures. The degree of freedom estimated the dimensions of the network, and further requirements to solve it with experimentally derived measurements. The lower degree of freedom, the easier solutions can be found. Analysis of conserved moiety revealed sets of compounds, whose sum always remained the same even under dynamic conditions. The number of conserved moieties equalled the number of linearly dependent balance equations. Analysis of conserved moieties can reveal false-positive conserved moieties which had accidentally appeared due to designed network structure. Accidental false-positive conserved moieties were eliminated and only biologically relevant conserved moieties (e.g. NAD(P)H/NAD(P), acetyl-CoA/CoA/succinyl-CoA) were remain in the network. Analysis of blocked transformers identified regions of the network, which were not capable of operating at steady state due to the irreversibility constraints of the transformer equation. Computation of the base vector of the nullspace revealed “functional” and “non-functional” regions of a metabolic network. Nullspace of S is a set of all solutions of the linear equation system S × v = 0. S was the stoichiometric matrix of the network (Supplemental data) and v was the rate vector of metabolic fluxes. Computation of parallel routes and cycles revealed a set of reactions/transporters/polymerizators which (i) were capable of maintaining a steady state, (ii) could not be decomposed, and (iii) in total did not consume or produce an external substrate/product. From condition (iii) followed that the net reaction of parallel routes and reaction cycles was zero. In the complex networks, there was a high chance for the accidental formation of parallel routes and cycles which increased internal degree of freedom. Parallel routes and cycles substantially increase the number of elementary flux modes without leading to new phenotypic behaviour. Therefore, their number had to be minimized in order to use elementary mode analysis for evaluation of maximal product yield. Computation of elementary modes revealed a unique set of smallest sub-networks that allowed a reconstructed network to function in steady state. Elementary mode analysis took into account stoichiometry and thermodynamics when evaluating whether a particular metabolic route or (sub-)network was feasible and likely for a given set of proteins/enzymes. The topological analysis of the model was performed with Insilico Discovery™ (Insilico Biotechnology AG, Stuttgart, Germany).
Flux Balance Analysis (FBA)
D-glycerate 3-phosphate (ChEBI:58272)
Dihydroxyacetone phosphate (ChEBI:57642)
D-glyceraldehyde 3-phosphate (ChEBI:59776)
- H+ :
- hv :
Light photon (ChEBI:302012)
is a nitrogen source, in the model it is NO 3 −
is a phosphorus source, in the model it is HPO 4 (2−) , i.e. orthophosphate
- Pi or HPO 4 (2−) :
is a sulphur source, in the model it is SO 4 2 −
Fructose-bisphosphate aldolase (EC 126.96.36.199)
Fructose 1,6-bisphosphatase (EC 188.8.131.52)
Ferredoxin-NADP-reductase (EC 184.108.40.206)
Glyceraldehyde-3-phosphate dehydrogenase (EC 220.127.116.11/13)
Glucokinase (EC 18.104.22.168)
Malate dehydrogenase (EC 22.214.171.124)
6-phosphofructokinase (EC 126.96.36.199)
Phosphoglycerate kinase (EC 188.8.131.52)
Phosphoglucomutase phosphoglucomutase (EC 184.108.40.206)
Phosphoribulokinase (EC 220.127.116.11)
Ribulose-bisphosphate carboxylase, RuBisCo (EC 18.104.22.168)
- SUC and STP:
Sucrose efflux transporter (Sugars Will Eventually be Exported Transporters)
Transketolase (EC 22.214.171.124)
Triosephosphate isomerase (EC 126.96.36.199)
Photosynthesis light reactions
Pentose phosphate pathway
- gdw :
Dry weight of biomass [g]
Flux balance analysis
Biological process- sucrose transport – The directed movement of sucrose into, out of or within a cell, or between cells by means of some agent such as a transporter or pore
The authors would like to thank Azam Salimi (University of Hohenheim, Stuttgart, Germany) for the assistance in curation of reaction and metabolite databases relevant for the model.
This work was supported by a research grant of the German Research foundation (DFG) to WXS.
Availability of data and materials
The metabolic multi-compartmental model of Arabidopsis thaliana in SBML and MATLAB formats.
High-resolution image of the modelled network [JPG].
Model documentation: detailed description of identifiers, stoichiometric equations, association with compartments, pathways and genes [1 file in ASCII format].
Exported stoichiometric matrix of the network [1 file in ASCII format].
FBA constraints for ‘light’ and ‘dark’ growth scenarios [2 files in ASCII formats].
FBA solutions for ‘light’ and ‘dark’ growth scenarios [6 files in ASCII formats].
Information on used in the model metabolites (e.g. elemental composition, charge, molecular weight, external database reference, etc.), transformers (e.g. identifier, trivial name, stoichiometric equation, belonging to pathway, EC number of the reaction, external database reference, etc.) and corresponding genes (gene ID in ATG format, catalysed reaction, gene definition, gene locus, etc.) [3 files in ASCII formats].
WS – idea generator, results interpretation and manuscript compilation; MZ – modelling, analysis, results interpretation and manuscript compilation; IM – database and model curation; YO – discussion of the modelling tools and results; IA – data visualization and results discussion; OK- data and model standards-compliant curation, annotation and database management. All authors have read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Taiz L, Zeiger E. Plant physiology, 5th edn. Sunderland: Sinauer Associates, Inc.; 2010.Google Scholar
- Scialdone A, Mugford ST, Feike D, Skeffington A, Borrill P, Graf A, Smith AM, Howard M. Arabidopsis plants perform arithmetic division to prevent starvation at night. eLife. 2013;2:e00669.PubMedPubMed CentralView ArticleGoogle Scholar
- Chapin FS, Schulze ED, Mooney HA. The ecology and economics of storage in plants. Annu Rev Ecol Syst. 1990;21:423–47.View ArticleGoogle Scholar
- Seaton DD, Ebenhöh O, Millar AJ, Pokhilko A. Regulatory principles and experimental approaches to the circadian control of starch turnover. J R Soc Interface. 2013;11(91):20130979.Google Scholar
- Caspar T, Huber SC, Somerville C. Alterations in growth, photosynthesis, and respiration in a starchless mutant of Arabidopsis thaliana (L.) deficient in chloroplast phosphoglucomutase activity. Plant Physiol. 1985;79:1–7.View ArticleGoogle Scholar
- Schulze W, Stitt M, Schulze ED, Neuhaus HE, Fichtner K. A quantification of the significance of assimilatory starch for growth of Arabidopsis thaliana L. Heynh. Plant Physiol. 1991;95:890–5.PubMedPubMed CentralView ArticleGoogle Scholar
- Schulze W, Schulze ED, Stader J, Heilmeier H, Stitt M, Mooney HA. Growth and reproduction of Arabidopsis thaliana in relation to storage of starch and nitrate in the wild-type and in starch-deficient and nitrate-uptake-deficient mutants. Plant Cell Environ. 1994;17:795–809.View ArticleGoogle Scholar
- Yazdanbakhsh N, Sulpice R, Graf A, Stitt M, Fisahn J. Circadian control of root elongation and C partitioning in Arabidopsis thaliana. Plant Cell Environ. 2001;34:877–94.View ArticleGoogle Scholar
- Lemoine R, La Camera S, Atanassova R, Dédaldéchamp F, Allario T, Pourtau N, Bonnemain J-L, Laloi M, Coutos-Thévenot P, Maurousset L et al. Source to sink transport and regulation by environmental factors. Front Plant Sci. 2013;4:272.Google Scholar
- Ludewig F, Flügge U-I. Role of metabolite transporters in source-sink carbon allocation. Front Plant Sci. 2013;4:231.Google Scholar
- Kühn C, Grof CPL. Sucrose transporters of higher plants. Curr Opin Plant Biol. 2010;13(3):287–97.View ArticleGoogle Scholar
- Flügge U-I, Häusler RE, Ludewig F, Gierth M. The role of transporters in supplying energy to plant plastids. J Exp Bot. 2011;62(7):2381–92.PubMedView ArticleGoogle Scholar
- Chen L-Q. SWEET sugar transporters for phloem transport and pathogen nutrition. New Phytol. 2014;201(4):1150–5.PubMedView ArticleGoogle Scholar
- Shiratake K. Genetics of sucrose transporter in plants. In. Genes, Genomes and Genomics. vol. 1. Ikenobe: Global Science Books, Ltd.; 2007:73–80.Google Scholar
- Boorer KJ, Loo DDF, Frommer WB, Wright EM. Transport mechanism of the cloned potato H+/sucrose cotransporter StSUT1. J Biol Chem. 1996;271:25139–44.PubMedView ArticleGoogle Scholar
- Gahrtz M, Stolz J, Sauer N. A phloem-specific sucrose-H+ symporter from Plantago major L. supports the model of apoplastic phloem loading. Plant J. 1994;6:697–706.PubMedView ArticleGoogle Scholar
- Schulze W, Weise A, Frommer WB, Ward JM. Function of the cytosolic N-terminus of sucrose transporter AtSUT2 in substrate affinity. FEBS Lett. 2000;485:189–94.PubMedView ArticleGoogle Scholar
- Gottwald JR, Krysan PJ, Young JC, Evert RF, Sussman MR. Genetic evidence for the in planta role of phloem-specific plasma membrane sucrose transporters. Proc Natl Acad Sci U S A. 2000;97(25):13979–84.PubMedPubMed CentralView ArticleGoogle Scholar
- Stadler R, Truernit E, Gahrtz M, Sauer N. The AtSUC1 sucrose carrier may represent the osmotic driving force for anther dehiscence and pollen tube growth in Arabidopsis. Plant J. 1999;19:269–78.PubMedView ArticleGoogle Scholar
- Chen LQ, Qu XQ, Hou BH, Sosso D, Osorio S, Fernie AR, Frommer WB. Sucrose efflux mediated by SWEET proteins as a key step for phloem transport. Science. 2012;35:207–11.View ArticleGoogle Scholar
- Zimmermann P, Hirsch-Hoffman M, Hennig L, Gruissem W. GENEVESTIGATOR. Arabidopsis microarray database and analysis toolbox. Plant Physiol. 2004;136:2621–32.PubMedPubMed CentralView ArticleGoogle Scholar
- Milne RJ, Byrt CS, Patrick JW, Grof CP. Are sucrose transporter expression profiles linked with patterns of biomass partitioning in Sorghum phenotypes? Front Plant Sci. 2013;4:223.Google Scholar
- Shachar-Hill Y. Metabolic network flux analysis for engineering plant systems. Curr Opin Biotechnol. 2013;24(2):247–55.PubMedView ArticleGoogle Scholar
- Sweetlove LJ, Obata T, Fernie AR. Systems analysis of metabolic phenotypes: what have we learnt? Trends Plant Sci. 2014;19(4):222–30.PubMedView ArticleGoogle Scholar
- Baghalian K, Hajirezaei M-R, Schreiber F. Plant Metabolic Modeling: Achieving New Insight into Metabolism and Metabolic Engineering. Plant Cell Online. 2014;26(10):3847–66.View ArticleGoogle Scholar
- Poolman MG, Assmus HE, Fell DA. Applications of metabolic modelling to plant metabolism. J Exp Bot. 2004;55(400):1177–86.PubMedView ArticleGoogle Scholar
- Arnold A, Nikoloski Z. Bottom-up Metabolic Reconstruction of Arabidopsis and Its Application to Determining the Metabolic Costs of Enzyme Production. Plant Physiol. 2014;165(3):1380–91.PubMedPubMed CentralView ArticleGoogle Scholar
- de Oliveira Dal’Molin CG, Quek L-E, 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Vanrolleghem PA, Heijnen J. Metabolic network modelling: Improving predictions of microbial metabolism by maximal incorporation of knowledge on biochemical reaction stoichiometry. In: Proceedings 9th Forum for Applied Biotechnology: 1995; Med. Fac. Landbouww. Univ. Gent, 60/4a, 1933–1940. 1995.Google Scholar
- Thiele I, Palsson BO. A protocol for generating a high-quality genome-scale metabolic reconstruction. Nat Protocols. 2010;5(1):93–121.PubMedView ArticleGoogle Scholar
- Stitt M, Sulpice R, Keurentjes J. Metabolic Networks: How to Identify Key Components in the Regulation of Metabolism and Growth. Plant Physiol. 2010;152(2):428–44.PubMedPubMed CentralView ArticleGoogle Scholar
- 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. 2012;109(1):339–44.PubMedView ArticleGoogle Scholar
- Poolman MG, Miguet L, Sweetlove LJ, Fell DA. A genome-scale metabolic model of Arabidopsis and some of Its properties. Plant Physiol. 2009;151(3):1570–81.PubMedPubMed CentralView ArticleGoogle Scholar
- Cheung CYM, 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.Google Scholar
- Gomes De Oliveira Dal’molin C, Quek L-E, Saa PA, Nielsen LK. A multi-tissue genome-scale metabolic modelling framework for the analysis of whole plant systems. Front Plant Sci. 2015;6:4.Google Scholar
- Grafahrend-Belau E, Junker A, Eschenröder A, Müller J, Schreiber F, Junker BH. Multiscale Metabolic Modeling: Dynamic Flux Balance Analysis on a Whole-Plant Scale. Plant Physiol. 2013;163(2):637–47.PubMedPubMed CentralView ArticleGoogle Scholar
- Poolman MG, Fell DA, Thomas S. Modelling photosynthesis and its control. J Exp Bot. 2000;51 suppl 1:319–28.PubMedView ArticleGoogle Scholar
- Arnold A, Nikoloski Z. A quantitative comparison of Calvin–Benson cycle models. Trends Plant Sci. 2011;16(12):676–83.PubMedView ArticleGoogle Scholar
- Sajitz-Hermstein M, Nikoloski Z. A novel approach for determining environment-specific protein costs: the case of Arabidopsis thaliana. Bioinformatics. 2010;26(18):i582–8.PubMedPubMed CentralView ArticleGoogle Scholar
- Williams TCR, Miguet L, Masakapalli SK, Kruger NJ, Sweetlove LJ, Ratcliffe RG. Metabolic Network Fluxes in Heterotrophic Arabidopsis Cells: Stability of the Flux Distribution under Different Oxygenation Conditions. Plant Physiol. 2008;148(2):704–18.PubMedPubMed CentralView ArticleGoogle Scholar
- Masakapalli SK, Le Lay P, Huddleston JE, Pollock NL, Kruger NJ, Ratcliffe RG. Subcellular flux analysis of central metabolism in a heterotrophic Arabidopsis cell suspension using steady-state stable isotope labeling. Plant Physiol. 2010;152(2):602–19.PubMedPubMed CentralView ArticleGoogle Scholar
- Nägele T, Henkel S, Hörmiller I, Sauter T, Sawodny O, Ederer M, Heyer AG. Mathematical Modeling of the Central Carbohydrate Metabolism in Arabidopsis Reveals a Substantial Regulatory Influence of Vacuolar Invertase on Whole Plant Carbon Metabolism. Plant Physiol. 2010;153(1):260–72.PubMedPubMed CentralView ArticleGoogle Scholar
- Schallau K, Junker BH. Simulating Plant Metabolic Pathways with Enzyme-Kinetic Models. Plant Physiol. 2010;152(4):1763–71.PubMedPubMed CentralView ArticleGoogle Scholar
- Rohwer JM. Kinetic modelling of plant metabolic pathways. J Exp Bot. 2012;63(6):2275–92.PubMedView ArticleGoogle Scholar
- Beauvoit BP, Colombié S, Monier A, Andrieu M-H, Biais B, Bénard C, Chéniclet C, Dieuaide-Noubhani M, Nazaret C, Mazat J-P, et al. Model-Assisted Analysis of Sugar Metabolism throughout Tomato Fruit Development Reveals Enzyme and Carrier Properties in Relation to Vacuole Expansion. Plant Cell. 2014;26(8):3224–42.PubMedPubMed CentralView ArticleGoogle Scholar
- Briskin DP, Reynolds-Niesman I. Determination of H+/ATP Stoichiometry for the Plasma Membrane H + −ATPase from Red Beet (Beta vulgaris L.) Storage Tissue. Plant Physiol. 1991;95(1):242–50.PubMedPubMed CentralView ArticleGoogle Scholar
- Schlesier B, Breton F, Mock HP. A hydroponic culture system for growing Arabidopsis thaliana plantlets under sterile conditions. Plant Mol Biol Report. 2003;21:449–56.View ArticleGoogle Scholar
- Zhang P, Foerster H, Tissier CP, Mueller L, Paley S, Karp PD, Rhee SY. MetaCyc and AraCyc. Metabolic Pathway Databases for Plant Research. Plant Physiol. 2005;138(1):27–37.PubMedPubMed CentralView ArticleGoogle Scholar
- Bauwe H, Hagemann M, Fernie AR. Photorespiration: players, partners and origin. Trends Plant Sci. 2010;15(6):330–6.PubMedView ArticleGoogle Scholar
- DalCorso G, Pesaresi P, Masiero S, Aseeva E, Schünemann D, Finazzi G, Joliot P, Barbato R, Leister D. A Complex Containing PGRL1 and PGR5 Is Involved in the Switch between Linear and Cyclic Electron Flow in Arabidopsis. Cell. 2008;132(2):273–85.Google Scholar
- Cleland RE, Bendall DS. Photosystem I cyclic electron transport: Measurement of ferredoxin-plastoquinone reductase activity. Photosynth Res. 1992;34(3):409–18.PubMedView ArticleGoogle Scholar
- Heber U, Bligny R, Streb P, Douce R. Photorespiration is Essential for the Protection of the Photosynthetic Apparatus of C3 Plants Against Photoinactivation Under Sunlight. Botanica Acta. 1996;109(4):307–15.View ArticleGoogle Scholar
- Sharkey TD. Estimating the rate of photorespiration in leaves. Physiol Plant. 1988;73(1):147–52.View ArticleGoogle Scholar
- Harbinson J, Foyer CH. Relationships between the Efficiencies of Photosystems I and II and Stromal Redox State in CO2-Free Air : Evidence for Cyclic Electron Flow in Vivo. Plant Physiol. 1991;97(1):41–9.PubMedPubMed CentralView ArticleGoogle Scholar
- Munekage Y, Hashimoto M, Miyake C, Tomizawa K-I, Endo T, Tasaka M, Shikanai T. Cyclic electron flow around photosystem I is essential for photosynthesis. Nature. 2004;429(6991):579–82.PubMedView ArticleGoogle Scholar
- Shikanai T. Cyclic Electron Transport Around Photosystem I: Genetic Approaches. Annu Rev Plant Biol. 2007;58(1):199–217.PubMedView ArticleGoogle Scholar
- Pänke O, Rumberg B. Energy and entropy balance of ATP synthesis. Biochim Biophys Acta Bioenergetics. 1997;1322(2–3):183–94.View ArticleGoogle Scholar
- Pänke O, Rumberg B. Kinetic modeling of rotary CF0F1-ATP synthase: storage of elastic energy during energy transduction. Biochim Biophys Acta Bioenergetics. 1999;1412(2):118–28.View ArticleGoogle Scholar
- Loizeau K, Gambonnet B, Zhang G-F, Curien G, Jabrin S, Van Der Straeten D, Lambert WE, Rébeillé F, Ravanel S. Regulation of One-Carbon Metabolism in Arabidopsis: The N-Terminal Regulatory Domain of Cystathionine γ-Synthase Is Cleaved in Response to Folate Starvation. Plant Physiol. 2007;145(2):491–503.PubMedPubMed CentralView ArticleGoogle Scholar
- Bohrer A-S, Kopriva S, Takahashi H. Plastid-cytosol partitioning and integration of metabolic pathways for APS/PAPS biosynthesis in Arabidopsis thaliana. Front Plant Sci. 2015;5:751.Google Scholar
- Facchinelli F, Weber APM. The metabolite transporters of the plastid envelope: an update. Front Plant Sci. 2011;2:Article 50.View ArticleGoogle Scholar
- Miernyk JA. The Isolation and Characterization on Nongreen Plastids. In. Cell Components. Edited by Linskens H-F, Jackson J, vol. 1. Berlin: Springer Berlin Germany; 1985:259–95.Google Scholar
- Hoefnagel MHN, Atkin OK, Wiskich JT. Interdependence between chloroplasts and mitochondria in the light and the dark. Biochim Biophys Acta Bioenergetics. 1998;1366(3):235–55.View ArticleGoogle Scholar
- Scheibe R. Malate valves to balance cellular energy supply. Physiol Plant. 2004;120(1):21–6.PubMedView ArticleGoogle Scholar
- Atwell BJ, Kriedemann PE, Turnbull CGN. Plants in Action: Adaptation in Nature, Performance in Cultivation. Sydney: Macmillan Education AU; 1999.Google Scholar
- Flügge UI, Heldt HW. The phosphate-triose phosphate-phosphoglycerate translocator of the chloroplast. Trends Biochem Sci. 1984;9(12):530–3.View ArticleGoogle Scholar
- Flügge U-I. Phosphate translocators in plastids. Annu Rev Plant Physiol Plant Mol Biol. 1999;50(1):27–45.PubMedView ArticleGoogle Scholar
- Flügge U-L. Phosphate translocation in the regulation of photosynthesis. J Exp Bot. 1995;46(special issue):1317–23.View ArticleGoogle Scholar
- Roels JA. Energetics and kinetics in biotechnology. Amsterdam: Elsevier Biomedical Press; 1983.Google Scholar
- Stephanopoulos GN, Aristidou AA, Nielsen J. Metabolic engineering: Principles and methodologies. San Diego: Academic; 1998.Google Scholar
- Nielsen J, Villadsen J, Liden G. Bioreaction engineernig principles, vol. 2nd. New York: Kluwer Academic/Plenum Publishers; 2003.View ArticleGoogle Scholar
- Villadsen J, Nielsen J, Liden G. Bioreaction engineering principles, vol. 3rd. Berlin: Springer; 2011.Google Scholar
- Portis AR, McCarty RE. Quantitative relationships between phosphorylation, electron flow, and internal hydrogen ion concentrations in spinach chloroplasts. J Biol Chem. 1976;251(6):1610–7.PubMedGoogle Scholar
- Badger MR. Photosynthetic oxygen exchange. Annu Rev Plant Physiol. 1985;36(1):27–53.View ArticleGoogle Scholar
- Lee CP, Gu Q, Xiong Y, Mitchell RA, Ernster L. P/O ratios reassessed: mitochondrial P/O ratios consistently exceed 1.5 with succinate and 2.5 with NAD-linked substrates. FASEB J. 1996;10(2):345–50.PubMedGoogle Scholar
- Karp PD, Paley SM, Krummenacker M, Latendresse M, Dale JM, Lee TJ, Kaipa P, Gilham F, Spaulding A, Popescu L, et al. Pathway Tools version 13.0: integrated software for pathway/genome informatics and systems biology. Brief Bioinform. 2010;11(1):40–79.PubMedView ArticleGoogle Scholar
- Huala E, Dickerman AW, Garcia-Hernandez M, Weems D, Reiser L, LaFond F, Hanley D, Kiphart D, Zhuang M, Huang W, et al. The Arabidopsis Information Resource (TAIR): a comprehensive database and web-based information retrieval, analysis, and visualization system for a model plant. Nucleic Acids Res. 2001;29(1):102–5.PubMedPubMed CentralView ArticleGoogle Scholar
- Swarbreck D, Wilks C, Lamesch P, Berardini TZ, Garcia-Hernandez M, Foerster H, Li D, Meyer T, Muller R, Ploetz L et al. The Arabidopsis Information Resource (TAIR): gene structure and function annotation. Nucleic Acids Res. 2008;36(D1009–1014).Google Scholar
- Kanehisa M, Goto S, Sato Y, Kawashima M, Furumichi M, Tanabe M. Data, information, knowledge and principle: back to metabolism in KEGG. Nucleic Acids Res. 2014;42(D1):D199–205.PubMedView ArticleGoogle Scholar
- Hastings J, de Matos P, Dekker A, Ennis M, Harsha B, Kale N, Muthukrishnan V, Owen G, Turner S, Williams M, et al. The ChEBI reference database and ontology for biologically relevant chemistry: enhancements for 2013. Nucleic Acids Res. 2013;41(D1):D456–63.PubMedView ArticleGoogle Scholar
- Piques M, Schulze WX, Höhne M, Usadel B, Gibon Y, Rohwer J, Stitt M. Ribosome and transcript copy numbers, polysome occupancy and enzyme dynamics in Arabidopsis. Mol Syst Biol. 2009;5(1):D456–463.Google Scholar
- Wolstencroft K, Owen S, Krebs O, Nguyen Q, Stanford NJ, Golebiewski M, Weidemann A, Bittkowski M, An L, Shockley D, et al. SEEK: a systems biology data and model management platform. BMC Syst Biol. 2015;9(1):1–12.View ArticleGoogle Scholar
- Roberts I, Smith S, De Rybel B, Van Den Broeke J, Smet W, De Cokere S, Mispelaere M, De Smet I, Beeckman T. The CEP family in land plants: evolutionary analyses, expression studies, and role in Arabidopsis shoot development. J Exp Bot. 2013;64(17):5371–81.PubMedView ArticleGoogle Scholar
- Schmid M, Davison TS, Henz SR, Pape UJ, Demar M, Vingron M, Scholkopf B, Weigel D, Lohmann JU. A gene expression map of Arabidopsis thaliana development. Nat Genet. 2005;37(5):501–6.PubMedView ArticleGoogle Scholar
- Kilian J, Whitehead D, Horak J, Wanke D, Weinl S, Batistic O, D’Angelo C, Bornberg-Bauer E, Kudla J, Harter K. The AtGenExpress global stress expression data set: protocols, evaluation and model data analysis of UV-B light, drought and cold stress responses. Plant J. 2007;50(2):347–63.PubMedView ArticleGoogle Scholar