- Research article
- Open Access
Modeling the ascorbate-glutathione cycle in chloroplasts under light/dark conditions
BMC Systems Biologyvolume 10, Article number: 11 (2016)
Light/dark cycles are probably the most important environmental signals that regulate plant development. Light is essential for photosynthesis, but an excess, in combination with the unavoidable presence of atmospheric oxygen inside the chloroplast, leads to excessive reactive oxygen species production. Among the defense mechanisms that activate plants to cope with environmental stress situations, it is worth noting the ascorbate-glutathione cycle, a complex metabolic pathway in which a variety of photochemical, chemical and enzymatic steps are involved.
We herein studied the dynamic behavior of this pathway under light/dark conditions and for several consecutive days. For this purpose, a mathematical model was developed including a variable electron source with a rate law proportional to the intensity of solar irradiance during the photoperiod, and which is continuously turned off at night and on again the next day. The model is defined by a nonlinear system of ordinary differential equations with an on/off time-dependent input, including a parameter to simulate the fact that the photoperiod length is not constant throughout the year, and which takes into account the particular experimental kinetics of each enzyme involved in the pathway. Unlike previous models, which have only provided steady-state solutions, the present model is able to simulate diurnal fluctuations in the metabolite concentrations, fluxes and enzymatic rates involved in the network.
The obtained results are broadly consistent with experimental observations and highlight the key role played by ascorbate recycling for plants to adapt to their surrounding environment. This approach provides a new strategy to in vivo studies to analyze plant defense mechanisms against oxidative stress induced by external changes, which can also be extrapolated to other complex metabolic pathways to constitute a useful tool to the scientific community in general.
The Earth takes approximately 24 h to make one complete turn around its own axis, which gives rise to the succession of day and night. In the daytime, the amount of sunlight that reaches the surface of our planet is not constant, but varies according to different factors such as time of day, season, altitude, latitude and atmospheric composition. The length of the light and darkness periods is not the same, and neither is constant throughout the year. Living organisms have evolved to coordinate their activities with light/dark cycles, which greatly influence many aspects of their metabolism, physiology, and even their behavior.
In plants, chloroplasts are the organelles that capture sunlight energy and store it as chemical energy to be used in photosynthesis. Plants have evolved different mechanisms to cope with natural fluctuations in light intensity to not only be able to harvest light optimally, but to also protect themselves from excess light . When the light absorbed/light used for CO2 fixation ratio is above 1, the generation of reactive oxygen species (ROS) is greatly accelerated, which will lead to the inhibition of photosynthetic machinery. Under such conditions, plants can activate different protective mechanisms to dissipate excess photon energy, including photorespiration, down-regulation of photosystem II (PSII) through proton gradient generation across thylakoid membranes, cyclic electron flow around the PSI or the water-water cycle [2, 3]. In addition, plants can also dissipate excess excitation energy in PSII antenna as heat through nonphotochemical processes, which involves the xanthophyll cycle using ascorbate (ASC) as the reducing agent .
ROS are generally toxic and can cause damage to different biomolecules. So their fast removal is crucial for plant survival. For this purpose, plants have developed a set of antioxidants such as glutathione (GSH) and ASC, as well as antioxidant enzymes such as superoxide dismutase (SOD) and ascorbate peroxidase (APX). GSH, ASC and APX work in the so-called ASC-GSH cycle together with other enzymes, which include glutathione reductase (GR), dehydroascorbate reductase (DHAR) and monodehydroascorbate reductase (MDAR) [2, 3]. The metabolic pathway consists in two coupled substrate cycles  (Fig. 1a) to achieve high sensitivity and amplification of the response [6, 7]. The dynamic behavior of this pathway has been previously analyzed by computer simulation [5, 8], which makes it possible to calculate concentrations of the species and fluxes of the steps involved in the proposed mathematical models in both unstressed and stressed chloroplasts. However, the reaction schemes used in these reports included a constant electron source, so they did not address the effect of the variations of solar irradiance intensity and were only valid under fixed light conditions. Nevertheless, it has been shown that the antioxidants concentrations and enzymatic activities involved in the detoxifying pathway are subject to diurnal fluctuations, with significant changes relating to changes in light intensity over the 24 h period [9–11]. Furthermore, it has also been reported that ROS responsive genes are not under circadian control, but their transcription is induced upon stress, and that the activity of many ROS-scavenging enzymes depends on the redox state of the plant cell . Understanding all the simultaneous metabolic responses against the different light/dark conditions is currently one of the major challenges in plant research to improve crop productivity under a changing global climate.
In the present paper, our goal was to develop a mathematical model based on ordinary differential equations (ODEs) that would allow us to simulate the dynamic behavior of the species involved in the ASC-GSH pathway for several consecutive days when an electron flux is introduced into the model by means of a rate law proportional to the intensity of solar irradiance, which is switched off under dark conditions. Light is one of the most important environmental factors to affect plant growth. Hence it should be considered in mathematical models to gain insights into the dynamics of plant physiological processes. The model developed herein represents an open system that allows the analysis, from a quantitative viewpoint, of the effect of such an environmentally relevant factor like sunlight on the temporal dynamics of the metabolites and fluxes involved in a metabolic network, particularly the H2O2-detoxifying pathway in chloroplasts. This tool is straightforwardly extendible to other light-dependent biological processes that occur in plants and other living organisms.
The metabolic pathway under study is shown in Fig. 1. First of all, the model includes an electron source, PSI, which is active during the photoperiod, but is switched off at night. Therefore, the electron flow from PSI, F 1 , is an input into the system which has been included as an on/off switch. It is a time-dependent variable that has been divided into three competitive routes, which obviously do not work under dark conditions. They are the following:
photoreduction of O2 to O 2 · − (one-electron), whose flux is F 11
photoreduction of NADP+ to NADPH (two-electrons), whose flux is F 12
photoreduction of monodehydroascorbate (MDA) to ASC (one-electron), whose flux is F 13 ,
where F 1 = F 11 + F 12 + F 13 . To express this fact mathematically, three distribution coefficients were defined, c 11 , c 12 and c 13 , respectively, whose sum equals unity; i.e. c 11 + c 12 + c 13 = 1. These coefficients are the following:
In this way, the expression for each electron flow is as follows:
Therefore the distribution of the electron flow among the three branches is not always the same, but varies throughout the photoperiod, and could even differ from day to day, depending on the redox conditions of the chloroplast (particularly on the MDA and NADP+ concentrations). This confers the model much flexibility. The O2 and CO2 concentrations were considered constant in the model.
To establish the most appropriate mathematical function that describes variation in F 1 with time during the photoperiod, the following assumptions were made: F 1 is zero at the beginning and the end of the photoperiod, it shows its maximum at noon, and it follows a rate law proportional to the solar irradiance intensity. So experimental data of average daily global solar irradiance for different months and cities were taken from  and were statistically fitted to different kinds of mathematical equations. The next sinusoidal function was chosen because it provided a good enough fit to experimental data (an example of the goodness-of-fit of real average daily global solar irradiance data to Eq. (11) is illustrated in Additional file 1, for an inland city of Spain):
Furthermore, the sinusoidal function offers the additional advantage of a clearly established physical meaning of the parameters (y 0 is a parameter to displace the sin function in the ordinate axis to avoid negative values, a i is the amplitude, b is the inverse of frequency and c is the phase). For greater comprehensibility, we used the condition y 0 = a i (a i > 0), so that the maximum value of F 1 was 2a i and the minimum value was 0. Therefore in our model, 2a i was the maximum value of F 1 at noon, b was the duration of the photoperiod in hours (if Δ = 0 and/or i = 0, see below) and c = 3π/2. A day counter, i, was included in the model to write the equations in a more simplified manner. Its expression is the following:
The a i value can be the same for every day. However, the intensity of solar irradiance is not the same every day; so F 1 varies from day to day. For this reason, a different value of a i was included in the model for each day, which was randomly obtained by a discrete uniform distribution in the COPASI software :
We must also take into account that the length of the photoperiod is not the same every day because it depends on the season of the year under study. This fact was included in the model by the parameter Δ, whose value can be positive or negative (or null) depending on whether the length of the photoperiod became longer or shorter (or constant) from one solstice to the other. It is important to note here that full day duration (light and dark conditions) was considered to be 24 h, where (b + iΔ) is the length of the photoperiod and 24 – (b + iΔ) the length of night. Simulation started and ended each photoperiod under dark conditions, i.e. F 1 (24i) = 0, F 1 (24i + b + i∆) = 0 and 2a i > F 1 (Time) > 0 if 24i < Time < 24i + b + i∆. After bearing all this in mind, F 1 was defined as follows:
One important feature of the pathway under study is the existence of three moiety-conserved cycles (metabolic structures interconverting different forms of a chemical moiety, while the sum of these forms remains constant [15–17]), which are the sums (NADPH + NADP+), (GSH + 2GSSG) and (ASC + DHA + MDA) (Fig. 1a). NADPH is photoproduced from NADP+ by light-dependent reactions (F 12 ); the antioxidant power of NADPH is then used to regenerate GSH and ASC by the GR- and MDAR-catalyzed reactions through the reduction of GSSG and MDA, respectively; an additional step of NADPH consumption by the Calvin-Benson cycle (and other electron-consuming reactions) with flux F N was also added to the model, as described in . The equation that describes F N is the following:
and k N,cte is the apparent rate constant for the flow F N .
GSH and GSSG are interconverted each other by the action of the enzymes GR and DHAR. MDA radicals are produced following ASC oxidation by the APX-catalyzed H2O2 reduction, and spontaneously undergo disproportionation into dehydroascorbate (DHA) and ASC. DHA is recycled into ASC by DHAR. ASC is also recovered following MDA reduction by both F 13 and MDAR. Given that MDA has been shown to be mainly photoreduced via ferredoxin (F 13 ), but not via NAD(P)H with MDAR, at least in the thylakoidal scavenging system , MDAR was only included into the model when indicated in the text, and NADPH was employed as the electron donor.
The four enzymes involved in the metabolic pathway, GR, DHAR, MDAR and APX, were considered active 24 h a day whenever their respective substrates were available. SOD was also included under the same conditions, but given its high catalytic efficiency , it was evident that it would work only under light conditions when superoxide radicals are generated in the model (F 11 ). Each enzyme was introduced into the model after taking into account its particular mechanism of action based on previously reported experimental data (for further details see ), namely for GR  and MDAR , the ping-pong mechanism:
For DHAR, a bi uni uni uni ping-pong mechanism :
For SOD, after taking into account its remarkably high activity as an O 2 · − scavenger, and the fact that the O 2 · − levels inside chloroplasts are well below K m (350 μM) , first-order kinetics were assumed:
In the particular case of APX, a more detailed mechanism was introduced with the next steps (Fig. 1b): 1) reaction of native enzyme with H2O2 to yield Compound I of APX (CoI); 2) oxidation of ASC to MDA by CoI; 3) oxidation of ASC to MDA by Compound II of APX (CoII); 4) H2O2-induced inactivation of the enzyme through CoI at low concentrations of ASC ; 5) continuous input of APX into the system at a non-constant rate, directly proportional to the difference between the initial enzyme concentration and the concentration of the active forms in each instant (see Eqn. (30)); the aim of this step was to obtain a stable level of APX activity under nonstress conditions. The APX rate equation is defined by the two steps of ASC oxidation to MDA, and then:
Other nonenzymatic reactions that participate in the pathway under study and which were included in our model are indicated in Fig. 1c (Eqs. (1)-(6)) [5, 8]. Briefly, spontaneous dismutation of MDA is necessary to supply DHA to DHAR to thus close the enzymatic cycle between APX and DHAR. Superoxide radicals spontaneously dismutate to O2 and H2O2 quite rapidly. DHA and superoxide are the most relevant molecules to contribute to the uncatalyzed GSSG production in vivo . The spontaneous oxidation of ASC by H2O2 and O 2 · − is also relevant. There are some reactions that proceed both spontaneously and enzymatically-driven, and therefore compete for the same substrates (SOD with Eq. (2), DHAR with Eq. (3), APX with Eq. (5)); in these cases, we numerically checked that the enzymatic reaction rate was higher than the rate of the corresponding chemical reaction under nonstress conditions. The rate constants corresponding to these reactions were considered apparent bimolecular rate constants. Most have been taken from the literature and are the same as those previously described (Table 1). However, the values for the macroscopic kinetic constants corresponding to a set of processes (i.e. min and max values for a i , k 11, k 12, k 13, k N,cte , k 5 APX ) were adapted to the present model in order to obtain reasonable results. The k 3 -value published , led to reaction (3) rates much higher than V DHAR , which considerably weakened the role of this enzyme in the pathway. For this reason, its value was adapted in the model. This possibility has been previously suggested  after considering that the rate constant for the reaction (3) in tissues would very likely be significantly lower than that determined in vitro in aqueous media.
Based on all these considerations, the complete system of ODEs corresponding to the model shown in Fig. 1 is the following:
Simulated progress curves were obtained by numerical solutions of this nonlinear set of ODEs with the initial conditions indicated in Table 2. Numerical integration was performed with the help of the COPASI 4.7 software (Build 34)  using a deterministic algorithm (LSODA) that is able to deal with stiff ODEs . COPASI is one of the most popular software applications for the simulation and analysis of biochemical networks. The resulting ODE model (15 days) consists in 13 species and 67 global quantities (kinetic parameters, enzymatic rates and fluxes).
Results and discussion
Model performance under nonstress conditions
Light-dependent reactions in chloroplasts use solar power to reduce NADP+ to NADPH, which is necessary for CO2 fixation. Photoreduction of molecular oxygen to O 2 · − and MDA to ASC also occurs competitively. Then the detoxifying efficiency of the ASC-GSH pathway can be evaluated by the relative distribution of the electron flux F 1 among the three competitive processes (F 11 , F 12 and F 13 ). Figure 2 shows the results obtained herein for F j (j = 1,11,12,13) (Fig. 2a and b) and c j (j = 11,12,13) (Fig. 2c) by computational simulation under nonstress conditions while running the model for 15 consecutive days, by beginning with a photoperiod of 11 h (b) and increasing this latter (Δ) by 0.25 h/day (which would correspond approximately to a time period from February to June in Spain ). The a i parameters were randomly obtained from Eq. (13) within the 3.5 x 106 and 4.5 x 106 μM h−1 range (which is equivalent to a maximum value of F 1 in the range ~1.9-2.5 mM s−1 ), so that each run of the model is different. Under these conditions, the distribution of the electron flow among the three routes throughout the photoperiod was not constant (Fig. 2b and c), and the photochemical reduction of NADP+ (F 12 ) was predominant despite a relatively low value of k 12 was chosen (Table 1). Therefore, NADPH was synthesized in the model upon electron transport activation at sunrise, and the NADPH/NADP+ ratio was sufficiently high under the light conditions (see also Fig. 3). Thus it is observed that c 12 was the most significant distribution coefficient throughout the photoperiod (c 11 = ~10 %, c 12 = ~78 %, c 13 = ~12 % at noon) (Fig. 2c). The O 2 · − photogeneration rate (F 11 ) was only slightly higher than F 13 at dawn and dusk (Fig. 2b), while during the remaining photoperiod, the photoreduction of ASC was higher than F 11 to maintain ascorbate homeostasis. This is in agreement with the experimental observations made in leaves and chloroplasts under normal physiological conditions . Under these parametric conditions, the maxima superoxide production rates were in the range 7.1 x 105-8.8 x 105 μM h−1 (~200-250 μM s−1), which are in the same order of magnitude as data previously reported [2, 8, 22].
Figure 3 shows the simulated progress curves that correspond to the chemical species and enzymatic activities involved in the model shown in Fig. 1 under the same parametric conditions as in Fig. 2. As expected according to the fact that midday electron fluxes were higher (in agreement with higher light intensities and temperatures), the concentrations and enzymatic activities involved in the ASC-GSH cycle peaked at noon, and were higher (or lower for NADP+, GSH, ASC and APX) than those in the early morning or near sunset. Two kinds of dynamic behavior can be observed from these plots. First, there were some chemical species like NADPH and NADP+, ASC, DHA and MDA, and enzymatic activities like V DHAR , whose levels were similar from day to day, regardless of solar irradiance intensity. Second, there was another group of chemical species and enzymatic rates, among which GSH and GSSG, H2O2, O 2 · − , APX, the Calvin-Benson cycle rate (F N ), V SOD , V GR and V APX were included, whose levels were more closely related to F 1 since a parallel variation was obtained. It is worth noting here that V GR and V DHAR should be equal in a stationary situation as a substrate cycle is established between them [6, 7]. However, the presence of the chemical steps in the model (Fig. 1c) led to a different situation, which was checked by simulations when these steps were absent (data not shown).
The numerical results indicated that ASC and GSH concentrations are maxima at night when ROS generation lowered. In contrast, the highest NADPH concentration was achieved during light conditions because its photoproduction from NADP+ predominated over its consumption for H2O2-detoxification (F 12 > F 11, see Fig. 2a and b). The H2O2 concentration was low early in the morning and peaked at noon before lowering again throughout the afternoon and evening. This situation agrees with the experimental data measured in mangrove leaves  and tobacco leaves . Greater APX inactivation was also obtained at noon as a consequence of the rise in the H2O2 levels (Fig. 3). At sunset, APX recovered through the de novo synthesis to reach to its original value under dark conditions. When a new photoperiod started, APX inactivation once again began. The amount of inactive APX (APXi) formed described a ladder-shaped plot where the greater F 1 was, the more pronounced steps were. Under these parametric conditions, APX was the most active enzyme in the cycle.
Model performance under high light conditions
In a continuously changing environment, plants find themselves under many different biotic and abiotic conditions from day to day, such as pathogen challenges, salinity, drought, temperature extremes, heavy metal toxicity, ozone, high-light intensity, etc., which can disrupt cell redox homeostasis and lead to substantial losses in crop yields and quality. Among these factors, light is an essential source of energy, but also a major source of abiotic stress which leads to ROS overgeneration  which, when combined with other stresses, is exacerbated . Therefore, in order to study the dynamic behavior of our model under stress conditions, the a i values were increased and randomly obtained within the range 4 x 106 to 5 x 106 μM h−1 from Eq. (13) (it is equivalent to maxima values of F 1 in the range ~2.2-2.8 mM s−1). To increase the stress, the distribution of the electron flow among the three routes involved in the model was also changed so that a higher production of O 2 · − was attained, along with greater APX photoinactivation, for which the rate constant corresponding to its de novo biosynthesis (k 5 APX ) was also decreased (the values of the kinetic parameters which have been changed compared to those indicated in Table 1 have been included in the Fig. 4 legend). Figure 4 shows the results thus obtained for 15 consecutive days, beginning with a photoperiod of 11 h (b) with an increase of this latter of 0.25 h/day (as indicated in Fig. 2). The data obtained indicated a predominant electron flow for the photochemical reduction of NADP+ (F 12 ), even though the k 12 -value was lowered to increase stress (c 11 = ~24 %, c 12 = ~55 %, c 13 = ~21 % at noon). Under these parametric conditions, the maxima superoxide production rates were in the range 2.0 x 106 – 2.4 x 106 μM h−1 (~550-660 μM s−1), which are in the order of magnitude as data previously reported for stressed chloroplasts .
Distribution coefficients evolved differently to the previous parametric conditions. At t = 0, F 1 was always split between F 11 and F 12 , with c 13 = 0 because it was assumed that [MDA]0 = 0 (Table 2). However under high light conditions, as F 11 worked very early in the day to generate O 2 · − , c 12 sharply increased to compensate c 11 , so c 12 reached a maximum and c 11 a minimum. Afterward, c 13 began to increase because MDA was produced at the expense of ASC oxidation, so c 12 decreased and c 11 remained approximately constant. In the afternoon, the system recovered until the production of free radicals stopped at night. Under these conditions, c 11 was always above c 13 .
Yet despite the system being subject to extremely stressful conditions, the NADPH levels were well buffered and regulated from day to day to peak at noon, although its concentration was lower than under the nonstress conditions (see Fig. 3). The differences between the maxima levels attained each day at noon were quite small. This result is interesting because it has been shown that chloroplastic NADPH concentrations change little after adjustment to light [31, 32], even in chloroplasts stressed with 10 μM paraquat , and also it is true that the NADPH levels are higher in the light than in the dark . The system became more stressed as the consumption of antioxidants GSH and ASC was greater, with higher concentrations of GSSG, DHA, MDA, APXi, H2O2 and O 2 · − attained at noon. The MDA levels achieved at noon remained approximately constant daily regardless of a i , which once again illustrates the system’s regulatory capacity. The enzymatic activities involved in the cycle were greater than under the previous conditions, which represents the defense mechanism of chloroplasts against increased ROS production. In fact, the overexpression of these enzymes has been shown to be linked to increased plant tolerance against different environmental stresses .
Effects on the system of varying antioxidant enzyme activities
Although chloroplastic APX is the primary target of photooxidative damage , the activity of the other enzymes in the pathway could also be decreased, particularly in combination with other stressful environmental factors. Therefore, computer simulations were also run under the above high-light conditions by decreasing the catalytic constant of each enzyme to investigate the effect produced on the metabolites and fluxes of the metabolic network.
Figure 5 shows the results obtained under GR-limiting conditions. The most marked effect was a substantial drop in the GSH/GSSG ratio at noon, which agrees with the physiological role of GR , and led to greater ASC consumption and higher H2O2 concentrations for the maximum solar irradiance hours (see Fig. 4). As a result, on day 14 (t = 320.4 h in the run shown herein), the H2O2 levels were well above a life-compatible level (in the order of mM), which gave rise to total APX inactivation and ASC depletion. The amount of formed APXi was larger than that obtained under the previous conditions, and the O 2 · − levels were higher (c 11 shows an increase at noon under these stressful conditions), as was SOD activity. However less V DHAR was observed, and even showed a minimum at noon as a result of lower GSH availability, which was accompanied by a subsequent increase in the noncatalyzed parallel reaction (3) rate (V (3) ) to counteract this drop in V DHAR .
MDAR is another important enzyme in the H2O2-detoxifying pathway to help regenerate ASC from MDA at the expense of NAD(P)H. Therefore, and given the importance of ASC availability to keep APX active, the enzyme was included in the model at this point under these highly stressful parametric conditions. The results yielded by the computer after numerous runs indicated the vital role of ASC recycling in response to extremely hard environmental stress conditions, since the numerical data obtained (see Additional file 2) came much closer to those obtained under the nonstress conditions; e.g., the MDA maximum levels at noon lowered to 3–4.5 μM and antioxidants consumption was much lower. Furthermore, ASC recovery through MDAR also led to a slower electron flow for the photochemical MDA reduction (F 13 ), along with a subsequent increase in F 12 to regenerate the excess of NADPH consumed by MDAR.
When k SOD was lowered in the system, the most meaningful effect was a considerable drop in the NADPH/NADP+ ratio (Fig. 6), higher MDA levels attained at noon and an increase in the rates of the uncatalyzed reactions (4) and (6) to remove the excess of O 2 · − at the expense of GSH and ASC. The response of the system was to slow down the electron flow for superoxide photogeneration, with a minimum in c 11 at noon (c 11 = ~20 %, c 12 = ~59.5 %, c 13 = 20.5 % at noon). An increase in the rates of the other enzymes, DHAR, GR and APX was also observed, which also helped compensate the poorer catalytic activity of SOD.
A decrease in k cat DHAR led to a greater consumption of antioxidants ASC and NADPH, although the GSH levels remained approximately constant (Fig. 7). As in the previous case, the rate of the spontaneous reaction (3), which is parallel to that catalyzed by DHAR, became higher. Increased ASC consumption led to greater APX inactivation at noon, while H2O2 and O 2 · − levels were higher (see Fig. 4), although the system showed sufficient capacity to remove it.
For these two enzymes, SOD and DHAR, the catalytic constants could be decreased by several orders of magnitude without a significant amount of H2O2 being accumulated due to the counteracting effect of the spontaneous parallel steps.
Model performance under NADP+-limiting conditions
In water or salt stress situations, a reduction in stomatal conductance occurs to limit water loss [36, 37]. Stomatal closure lowers the CO2 concentration available for the Calvin-Benson cycle, and consequently NADP+ availability to accept the electrons from PSI . This means that more electrons from the electron transport chain in the chloroplast are used for O2 photoreduction, which leads to greater ROS generation (redox poising) [38, 39]. This was simulated in the model by decreasing the apparent rate constant for the Calvin-Benson cycle, k N,cte , and k cat GR , both under the standard conditions indicated in Table 1, and under the stressful high-light conditions (the parametric conditions indicated in Fig. 5), in the absence and presence of MDAR.
Under the standard conditions indicated in Table 1, and in the absence of MDAR, a change in the distribution of the electron flow from PSI was observed, which was in agreement to that previously indicated as F 13 > F 11 > F 12 (Fig. 8). NADP+ levels were in fact considerably lower than those obtained under nonstress conditions (see Fig. 3). The most important differences with respect to the data obtained under nonstress conditions were higher GSH consumption, greater H2O2 and O 2 · − concentrations at noon, F N considerably lowered, and V SOD and V APX were higher. However, MDA levels were in the same order of magnitude.
Under high-light conditions and in the absence of MDAR, the distribution of the electron flow was in the order of F 11 > F 13 > F 12 (i.e. the superoxide production rate was higher than the ASC and MDA photoreduction rates, except early in the morning) (Fig. 9), which resulted in extremely high H2O2 concentrations always being attained on the first day in each run. NADP+concentration increased early in the morning (in parallel with c 12 ), but quickly decreased. There is an instant (t = 2.84 h in the run here shown) in which APX was totally inactivated and ASC was totally depleted (although it is possible to observe an increase in the ASC concentration after that, it has no biological sense. It is due to the spontaneous redox steps involved in the model; however, it must be taken into account that all of these recoveries observed after APX inactivation are produced at H2O2 concentrations above a life-compatible level). In this way, simulations comprehensibly revealed the importance of a tight regulation of the NADPH/NADP+ ratio to ensure sufficient NADP+ levels available to accept electrons in the electron transport chain .
However, if MDAR was included under these fatal parametric conditions, the distribution of the electron flow changed to: F 12 > F 11 > F 13 (Fig. 10a). Once again, a drop in F 13 was observed at the expense of the increase in F 12 as a result of ASC recycling through MDAR. This situation led to physiologically more adequate levels of H2O2, O 2 · − and MDA at noon, and to a less significant consumption of GSH and ASC (although the system is still quite stressed, see Additional file 3). However, when this in silico experiment was repeated with a lower MDAR concentration, the distribution of the electron flow differed, and even differed from day to day at the time of maximum solar irradiance (Fig. 10b). In this scenario, the electron flow deviated mainly toward O 2 · − generation (F 11 ) every day, but at noon F 12 > F 11 was obtained for the days with lower solar irradiance (2, 3, 5, 6, 11, 15 in this run), whereas F 11 > F 12 for the days with greater solar irradiance (1, 4, 7, 8, 10, 12, 13, 14). ASC and GSH consumption was more meaningful (see Additional file 4), which led to higher levels of MDA, DHA and GSSG at noon.
All these results indicate that plants displaying increased activity in the Calvin-Benson cycle will better resist excessive light, i.e., the Calvin-Benson cycle, besides providing sugars, is also regenerating electron acceptors (NADP+ in the model) to thus represent chloroplast self-protection. In fact this is what actually happens. C4 and CAM plants are the best adapted to arid, hot, high-light environments because they possess greater photosynthetic efficiency than C3 plants . It has also been proposed that one of the reasons for the successful adaptation of Antarctic vascular plants to high light/low temperature conditions is the robustness of their CO2 assimilation machinery . This is also the role played by photorespiration in plants as a way to use the excess ATP and NADPH generated in thylakoids, and therefore to dissipate excess energy to prevent over-reduction of the photosynthetic electron chain and consequent damage in the photosynthetic apparatus . As a matter of fact, it is well-known that the Calvin-Benson cycle enzymes subject to regulation by reduced thioredoxin are activated by reduction in the light and deactivated by oxidation in the dark. This mechanism has been demonstrated for several chloroplast enzymes, including fructose-1,6-bisphosphatase, sedoheptulose-1,7-bisphosphatase , phosphoribulokinase, NADP+-glyceraldehyde-3-phosphate dehydrogenase, Rubisco activase and ATPsynthase . Evidently, this activation has a limit after which the electron transport chain becomes saturated and electrons are deflected toward dioxygen reduction.
By last, the mathematical model developed herein is a very useful tool to understand the elements that determine diurnal fluctuations in the chemical species involved in ROS generation-detoxification in chloroplasts. This study can be extended in different ways. One important aspect has been mentioned above, which would allow to make predictions of the adaptability of plants to different solar irradiance conditions. However, this would require knowledge of the model’s kinetic parameters per plant type. Another obvious aspect is to add new steps to the model to study the relation with other nearby metabolic pathways. The ASC-GSH pathway is closely related to the electron transport chain in the chloroplast, which is involved in ATP production for photosynthesis. Ongoing research looks at the relation between photosynthesis and the ASC-GSH pathway. Another very interesting aspect would be to consider the diffusion and compartmentalization processes of chemical species to model the asynchronous metabolite supplies to enzymes. Previously, we modeled the adenylate energy system and ATP production at the systemic cellular level with a system of delay-differential equations to take into account different time scales within the cell . Finally, the present model can be extrapolated to other metabolic pathways to quantitatively analyze the effect of sunlight on their metabolism.
In this paper, a mathematical model able to simulate not only the distribution of the electron flow from PSI in the chloroplast, but also the dynamics of the chemical and enzymatic reactions involved in the ASC-GSH pathway has been developed, which takes into account for the first time the succession of days and nights. The metabolic processes involved in the network have been described by a nonlinear system of ordinary differential equations in which the enzymatic rate equations and biochemical kinetic parameters have been retrieved from previously reported experimental data. One important novelty in the numerical results here shown is that metabolite concentrations and enzymatic activities do not evolve toward a steady-state level, but display oscillatory behavior, which is dependent on the time of day. The results here obtained clearly highlight the importance of the distribution of electron fluxes through the system for the detoxifying efficiency of the ASC-GSH pathway. The model can help to understand the elements that determine diurnal fluctuations in the chemical species involved in ROS generation-detoxification in chloroplasts. This approach can provide strategies to analyze plant defense mechanisms against oxidative stress and can be extrapolated to other metabolic pathways to quantitatively analyze the effect of sunlight on the metabolism.
Availability of supporting data
- APXi :
ascorbate peroxidase inactive
ordinary differential equations
reactive oxygen species
Müller-Moulé P, Golan T, Niyogi KK. Ascorbate-deficient mutants of Arabidopsis grow in high light despite chronic photooxidative stress. Plant Physiol. 2004;134:1163–72.
Asada K. The water-water cycle in chloroplasts: scavenging of active oxygens and dissipation of excess photons. Annu Rev Plant Physiol Plant Mol Biol. 1999;50:601–39.
Asada K. Production and scavenging of reactive oxygen species in chloroplasts and their functions. Plant Physiol. 2006;141:391–6.
Maxwell K, Johnson GN. Chlorophyll fluorescence-a practical guide. J Exp Bot. 2000;51:659–68.
Valero E, González-Sánchez MI, Macià H, García-Carmona F. Computer simulation of the dynamic behavior of the glutathione-ascorbate redox cycle in chloroplasts. Plant Physiol. 2009;149:1958–69.
Valero E, Varón R, García-Carmona F. Kinetic study of an enzymatic cycling system coupled to an enzymatic step. Determination of alkaline phosphatase activity. Biochem J. 1995;309:181–5.
Valero E, Varón R, García-Carmona F. Kinetic analysis of a model for double substrate cycling: highly amplified ADP (and/or ATP) quantification. Biophys J. 2004;86:3598–606.
Polle A. Dissecting the superoxide dismutase-ascorbate-glutathione pathway in chloroplasts by metabolic modeling: computer simulations as a step towards flux analysis. Plant Physiol. 2001;126:445–62.
Chang AC, Yang TY, Riskowski GL. Ascorbic acid, nitrate, and nitrite concentration relationship to the 24 h light/dark cycle for spinach grown in different conditions. Food Chem. 2013;138:382–8.
Cheeseman JM. Hydrogen peroxide concentrations in leaves under natural conditions. J Exp Bot. 2006;57:2435–44.
Noctor G, Arisi ACM, Jouanin L, Valadier MH, Roux Y, Foyer CH. Light-dependent modulation of foliar glutathione synthesis and associated amino acid metabolism in poplar overexpressing γ-glutamylcysteine synthetase. Planta. 1997;202:357–69.
Sánchez A, Shin J, Davis SJ. Abiotic stress and the plant circadian clock. Plant Signal Behav. 2011;6:223–31.
http://re.jrc.ec.europa.eu/pvgis/apps4/pvest.php. Accessed 10 Feb 2014.
http://www.copasi.org. Accessed 20 Mar 2015.
Reich JG, Sel’kov EE. Energy metabolism of the cell. New York: Academic; 1981.
Hofmeyr J-H S, Kacser H, Merwe KJ. Metabolic control analysis of moiety-conserved cycles. Eur J Biochem. 1986;155:631–41.
Valero E, Varón R, García-Carmona F. A kinetic study of a ternary cycle between adenine nucleotides. FEBS J. 2006;273:3598–613.
Bowler C, Van Montagu M, Inze D. Superoxide dismutase and stress tolerance. Annu Rev Plant Physiol Plant Mol Biol. 1992;43:83–116.
Vanoni MA, Wong KK, Ballou DP, Blanchard JS. Glutathione reductase: comparison of steady-state and rapid reaction primary kinetic isotope effects exhibited by the yeast, spinach, and Escherichia coli enzymes. Biochemistry. 1990;29:5790–6.
Hossain MA, Asada K. Monodehydroascorbate reductase from cucumber is a flavin adenine dinucleotide enzyme. J Biol Chem. 1985;260:12920–6.
Shimaoka T, Miyake C, Yokota A. Mechanism of the reaction catalyzed by dehydroascorbate reductase from spinach chloroplasts. Eur J Biochem. 2003;270:921–8.
Asada K, Takahashi M. Production and scavenging of active oxygen species in photosynthesis. In: Kyle D, Osmond C, Arntzen C, editors. Photoinhibition. New York: Elsevier Science Publishers; 1987. p. 227–87.
Miyake C, Asada K. Inactivation mechanism of ascorbate peroxidase at low concentrations of ascorbate; hydrogen peroxide decomposes compound I of ascorbate peroxidase. Plant Cell Physiol. 1996;37:423–30.
Rahantaniaina MS, Tuzet A, Mhamdi A, Noctor G. Missing links in understanding redox signaling via thiol/disulfide modulation: how is glutathione oxidized in plants? Front Plant Sci. 2013;4:477.
Hausladen A, Kunert KJ. Effects of artificially enhanced levels of ascorbate and glutathione on the enzymes monodehydroascorbate reductase, dehydroascorbate reductase and glutathione reductase in spinach (Spinacia oleracea). Physiol Plant. 1990;79:384–8.
Hoops S, Sahle S, Gauges R, Lee C, Pahle J, Simus N, et al. COPASI – A COmplex PAthway SImulator. Bioinformatics. 2006;22:3067–74.
Smirnoff N. Ascorbic acid: metabolism and functions of a multifaceted molecule. Curr Opin Plant Biol. 2011;3:229–35.
Chen Z, Gallie DR. The ascorbic acid redox state controls guard cell signaling and stomatal movement. Plant Cell. 2004;16:1143–62.
Oelze ML, Vogel MO, Alsharafa K, Kahmann U, Viehhauser A, Maurino VG, et al. Efficient acclimation of the chloroplast antioxidant defence of Arabidopsis thaliana leaves in response to a 10-or 100-fold light increment and the possible involvement of retrograde signals. J Exp Bot. 2012;63:1297–313.
Suzuki N, Rivero RM, Shulaev V, Blumwald E, Mittler R. Abiotic and biotic stress combinations. New Phytol. 2014;203:32–43.
Bielawski W, Joy KW. Reduced and oxidised glutathione and glutathione reductase activity in tissues of Pisum sativum. Planta. 1986;169:267–72.
Gerst U, Schönknecht G, Heber U. ATP and NADPH as the driving force of carbon reduction in leaves in relation to thylakoid energization by light. Planta. 1994;193:421–9.
Holfgrefe S, Backhausen JE, Kitzmann C, Scheibe R. Regulation of steady-state photosynthesis in isolated intact chloroplasts under constant light: responses of carbon fluxes, metabolite pools and enzyme-activation states to changes of electron pressure. Plant Cell Physiol. 1997;38:1207–16.
Berg JM, Tymoczko JL, Stryer L. The activity of the Calvin cycle depends on environmental conditions. In: Biochemistry. New York: WH Freeman and Company; 2002. p. 597–600.
Noctor G, Foyer CH. Ascorbate and glutathione: keeping active oxygen under control. Annu Rev Plant Physiol Plant Mol Biol. 1998;49:249–79.
Faize M, Burgos L, Faize L, Piqueras A, Nicolás E, Barba-Espín G, et al. Involvement of cytosolic ascorbate peroxidase and Cu/Zn-superoxide dismutase for improved tolerance against drought. J Exp Bot. 2011;62:2599–613.
Barba-Espín G, Clemente-Moreno MJ, Álvarez S, García-Legaz MF, Hernández JA, Díaz-Vivancos P. Salicylic acid negatively affects the response to salt stress in pea plants: effects on PR1b and MAPK expression. Plant Biol. 2011;13:909–17.
Robinson JM. Does O2 photoreduction occur within chloroplasts in vivo? Physiol Plant. 1988;72:666–80.
Mullineaux P, Karpinski S. Signal transduction in response to excess light: getting out of the chloroplast. Curr Opin Plant Biol. 2002;5:43–8.
Ashraf M, Harris PJC. Photosynthesis under stressful environments: An overview. Photosynthetica. 2013;51:163–90.
Pérez-Torres E, Bascuñán L, Sierra A, Bravo LA, Corcuera LJ. Robustness of activity of Calvin cycle enzymes after high light and low temperature conditions in Antarctic vascular plants. Polar Biol. 2006;29:909–16.
Asada K, Endo T, Mano J, Miyake C. Molecular mechanism for relaxation of and protection from light stress. In: Satoh K, Murata N, editors. Stress responses of photosynthetic organisms. Amsterdam: Elsevier; 1998. p. 37–52.
Harrison EP, Willingham NM, Lloyd JC, Raines CA. Reduced sedoheptulose-1,7-bisphosphatase levels in transgenic tobacco lead to decreased photosynthetic capacity and altered carbohydrate accumulation. Planta. 1998;204:27–36.
Buchanan BB, Gruissem W, Jones RL. Biochemistry and molecular biology of plants. Monona Drive, Rockville, MD, USA: American Society of Plant Physiologist; 2000.
De la Fuente I, Cortés JM, Valero E, Desroches M, Rodrigues S, Malaina I, et al. On the dynamics of the adenylate energy system: homeorhesis vs homeostasis. PLoS One. 2014;9(10), e108676. doi:10.1371/journal.pone.0108676.
Li C, Donizelli M, Rodriguez N, Dharuri H, Endler L, Chelliah V, et al. BioModels Database: An enhanced, curated and annotated resource for published quantitative kinetic models. BMC Syst Biol. 2010;4:92.
https://www.ebi.ac.uk/biomodels-main/ Accessed 5 Oct 2015.
Bielski BHJ, Comstock DA, Bowen RA. Ascorbic acid free radicals. I. Pulse radiolysis study of optical absorption and kinetic properties. J Am Chem Soc. 1971;93:5624–9.
Fridovich I. Oxygen toxicity: a radical explanation. J Exp Biol. 1998;201:1203–9.
Polle A, Junkermann W. Inhibition of apoplastic and symplastic peroxidase activity from Norway spruce by the photooxidant hydroxymethyl hydroperoxide. Plant Physiol. 1996;104:617–23.
This work was funded by the Spanish Ministry of Economy and Competitiveness (MINECO, http://www.mineco.gob.es/portal/site/mineco/idi), Project No. BFU2013-44095-P (cofunded with FEDER funds, EU), and by the Junta de Comunidades de Castilla-La Mancha (http://www.educa.jccm.es/idiuniv/es), Project No. PEII-2014-014-A (cofunded with FEDER funds, EU). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
The authors declare that they have no competing interests.
Conceived and designed the experiments: EV, HM, MIGS and FGC. Performed the experiments: EV, HM, MIGS. Analyzed the data: EV, HM, IMDF, JAH and FGC. Wrote the paper: EV, HM, IMDF, JAH and FGC. Supervised the research: EV, IMDF, JAH and FGC. Developed the computational model: EV, HM and FGC. All authors read and approved the final manuscript.
Fitting the average daily global clear-sky solar irradiance data to Eq. ( 11 ). Data were taken from  after considering these parameters: geographic coordinates = 40° 25' 0'' North, 3° 42' 1'' West (Madrid, Spain), month = September, inclination of plane = 35° and orientation (azimuth) of plane = 0°. Dots represent the real solar irradiance data (adapted so that the photoperiod starts at time = 0) and the line corresponds to the nonlinear regression analysis. Data were fitted by the SigmaPlot Scientific Graphing Software for Windows, version 13.0 (2014, Systat Software, Inc.). (TIF 187 kb)
Simulated progress curves under high-light and GR-limiting conditions in the presence of 2 μM MDAR. Other parametric conditions as indicated in Fig. 5. (TIF 1969 kb)
Simulated progress curves under NADP + -limiting and high-light conditions in the presence of 2 μM MDAR. Parametric conditions as indicated in Fig. 10A. (TIF 1806 kb)