C3 photosynthesis in silico

Assessing key parameters in simultaneous simulation of rapid kinetics of chlorophyll a fluorescence and trans-thylakoid electric potential difference

Our study attempts to address the following questions: among numerous photosynthetic modules, which parameters notably influence the rapid chlorophyll fluorescence (ChlF) rise, the so-called O-J-I-P transient, in conjunction with the P515 signal, as these two records are easily obtained and widely used in photosynthesis research, and how are these parameters ranked in terms of their importance? These questions might be difficult to answer solely through experimental assays. Therefore, we employed an established photosynthesis model. Firstly, we utilized the model to simulate the measured rapid ChlF rise and P515 kinetics simultaneously. Secondly, we employed the sensitivity analysis (SA) tool by randomly altering model parameters to observe their effects on model output variables. Thirdly, we systematically identified significant parameters for both or one of the kinetics across various scenarios. A novel aspect of our study is the application of the Morris method, a global SA tool, to simultaneously assess the significance of model parameters in shaping both or one of the kinetics. The Morris SA technique enables the quantification of how much a specific parameter affects O-J-I-P transient during particular time intervals (e.g., J, I, and P steps). This allowed us to theoretically analyze which step is more significantly influenced by the parameter. In summary, our study contributes to the field by providing a comprehensive analysis of photosynthesis kinetics and emphasizing the importance of parameter selection in modelling this process. These findings can inform future research efforts aimed at improving photosynthesis models and advancing our understanding of photosynthetic processes.

From leaf to multiscale models of photosynthesis: applications and challenges for crop improvement

To keep up with the growth of human population and to circumvent deleterious effects of global climate change, it is essential to enhance crop yield to achieve higher production. Here we review mathematical models of oxygenic photosynthesis that are extensively used, and discuss in depth a subset that accounts for diverse approaches providing solutions to our objective. These include models (1) to study different ways to enhance photosynthesis, such as fine-tuning antenna size, photoprotection and electron transport; (2) to bioengineer carbon metabolism; and (3) to evaluate the interactions between the process of photosynthesis and the seasonal crop dynamics, or those that have included statistical whole-genome prediction methods to quantify the impact of photosynthesis traits on the improvement of crop yield. We conclude by emphasizing that the results obtained in these studies clearly demonstrate that mathematical modelling is a key tool to examine different approaches to improve photosynthesis for better productivity, while effective multiscale crop models, especially those that also include remote sensing data, are indispensable to verify different strategies to obtain maximized crop yields.

Screening rate constants in the simulation of rapid kinetics of chlorophyll a fluorescence using the Morris method

Chlorophyll a fluorescence, a sensitive and cost-effective probe, is widely used in photosynthetic research. Its rapid phase, occurring within 1 second under intense illumination, displays complex O-J-I-P transients, providing valuable insights into various aspects of photosynthesis. In addition to employing experimental approaches to measure the rapid Fluorescence Induction (FI) kinetics, mathematical modeling serves as a crucial tool for understanding the underlying mechanisms that drive FI dynamics. However, the significant uncertainty and arbitrary nature of selecting model parameters amplify concerns about the effectiveness of modeling tools in aiding photosynthesis research. Therefore, there is a need to gain a deeper understanding of how these models operate and how arbitrary parameter choices may influence their outcomes. In this study, we employed the Morris method, a global Sensitivity Analysis (SA) tool, to assess the significance of rate constants employed in an existing fluorescence model, particularly those linked to the entire electron transport chain, in shaping the rapid FI dynamics. In summary, utilizing the insights gained from the Morris SA allows for targeted refinement of the photosynthesis model, thereby improving our understanding of the complex processes inherent in photosynthetic systems.

Dynamic modulation of transthylakoid electric potential by chloroplast ATP synthases

The light-induced transthylakoid membrane potential (ΔΨm) can function as a driving force to help catalyzing the formation of ATP molecules, proving a tight connection between ΔΨm and the ATP synthase. Naturally, a question can be raised on the effects of altered functioning of ATP synthases on regulating ΔΨm, which is attractive in the area of photosynthetic research. Lots of findings, when making efforts of solving this difficulty, can offer an in-depth understanding into the mechanism behind. However, the functional network on modulating ΔΨm is highly interdependent. It is difficult to comprehend the consequences of altered activity of ATP synthases on adjusting ΔΨm because parameters that have influences on ΔΨm would themselves be affected by ΔΨm. In this work, a computer model was applied to check the kinetic changes in polarization/depolarization across the thylakoid membrane (TM) regulated by the modified action of ATP synthases. The computing data revealed that under the extreme condition by numerically "switching off" the action of the ATP synthase, the complete inactivation of ATP synthase would markedly impede proton translocation at the cytb6f complex. Concurrently, the KEA3 (CLCe) porter, actively pumping protons into the stroma, further contributes to achieving a sustained low level of ΔΨm. Besides, the quantitative consequences on every particular component of ΔΨm adjusted by the modified functioning of ATP synthases were also explored. By employing the model, we bring evidence from the theoretical perspective that the ATP synthase is a key factor in forming a transmembrane proton loop thereby maintaining a propriate steady-state ΔΨm to meet variable environmental conditions.

Comparative modeling of fluorescence and P700 induction kinetics for alga Scenedesmus sp. obliques and cyanobacterium Synechocystis sp. PCC 6803. Role of state 2-state 1 transitions and redox state of plastoquinone pool

The model of thylakoid membrane system (T-M model) (Belyaeva et al. Photosynth Res 2019, 140:1-19) has been improved in order to analyze the induction data for dark-adapted samples of algal (Scenedesmus obliques) and cyanobacterial (Synechocystis sp. PCC 6803) cells. The fluorescence induction (FI) curves of Scenedesmus were measured at light exposures of 5 min, while FI and P700 redox transformations of Synechocystis were recorded in parallel for 100 s intervals. Kinetic data comprising the OJIP-SMT fluorescence induction and OABCDEF P700+ absorbance changes were used to study the processes underlying state transitions qT2→1 and qT1→2 associated with the increase/decrease in Chl fluorescence emission. A formula with the Hill kinetics (Ebenhöh et al. Philos Trans R Soc B 2014, 369:20130223) was introduced into the T-M model, with a new variable to imitate the flexible size of antenna AntM(t) associated with PSII. Simulations revealed that the light-harvesting capacity of PSII increases with a corresponding decrease for that of PSI upon the qT2→1 transition induced by plastoquinone (PQ) pool oxidation. The complete T-M model fittings were attained on Scenedesmus or Synechocystis fast waves OJIPS of FI, while SMT wave of FI was reproduced at intervals shorter than 5 min. Also the fast P700 redox transitions (OABC) for Synechocystis were fitted exactly. Reasonable sets of algal and cyanobacterial electron/proton transfer (ET/PT) parameters were found. In the case of Scenedesmus, ET/PT traits remained the same irrespective of modeling with or without qT2→1 transitions. Simulations indicated a high extent (20%) of the PQ pool reduction under dark conditions in Synechocystis compared to 2% in Scenedesmus.

Modelling plants across scales of biological organisation for guiding crop improvement

Grain yield improvement in globally important staple crops is critical in the coming decades if production is to keep pace with growing demand; so there is increasing interest in understanding and manipulating plant growth and developmental traits for better crop productivity. However, this is confounded by complex cross-scale feedback regulations and a limited ability to evaluate the consequences of manipulation on crop production. Plant/crop modelling could hold the key to deepening our understanding of dynamic trait-crop-environment interactions and predictive capabilities for supporting genetic manipulation. Using photosynthesis and crop growth as an example, this review summarises past and present experimental and modelling work, bringing about a model-guided crop improvement thrust, encompassing research into: (1) advancing cross-scale plant/crop modelling that connects across biological scales of organisation using a trait dissection-integration modelling principle; (2) improving the reliability of predicted molecular-trait-crop-environment system dynamics with experimental validation; and (3) innovative model application in synergy with cross-scale experimentation to evaluate G×M×E and predict yield outcomes of genetic intervention (or lack of it) for strategising further molecular and breeding efforts. The possible future roles of cross-scale plant/crop modelling in maximising crop improvement are discussed.

An exploratory steady-state redox model of photosynthetic linear electron transport for use in complete modelling of photosynthesis for broad applications

A photochemical model of photosynthetic electron transport (PET) is needed to integrate photophysics, photochemistry, and biochemistry to determine redox conditions of electron carriers and enzymes for plant stress assessment and mechanistically link sun-induced chlorophyll fluorescence to carbon assimilation for remotely sensing photosynthesis. Towards this goal, we derived photochemical equations governing the states and redox reactions of complexes and electron carriers along the PET chain. These equations allow the redox conditions of the mobile plastoquinone pool and the cytochrome b₆ f complex (Cyt) to be inferred with typical fluorometry. The equations agreed well with fluorometry measurements from diverse C₃ /C₄ species across environments in the relationship between the PET rate and fraction of open photosystem II reaction centres. We found the oxidation of plastoquinol by Cyt is the bottleneck of PET, and genetically improving the oxidation of plastoquinol by Cyt may enhance the efficiency of PET and photosynthesis across species. Redox reactions and photochemical and biochemical interactions are highly redundant in their complex controls of PET. Although individual reaction rate constants cannot be resolved, they appear in parameter groups which can be collectively inferred with fluorometry measurements for broad applications. The new photochemical model developed enables advances in different fronts of photosynthesis research.

Prying into the green black-box

Life-long efforts of the Tartu photosynthesis research group have been summarized. The measurements were facilitated by self-designed instruments, distinct in multifunctionality and fastresponse time. The black-box type kinetical analysis on intact leaves has revealed several physiologically significant features of leaf photosynthesis. Rubisco studies reflected competition for the active site between the substrates and products, linearizing in vivo kinetics compared with the low-K_m in vitro responses. Rubisco Activase usually activates only a small part of the Rubisco, making the rest of it a storage protein. Precisely quantifying absorbed photons and the responding transmittance changes, electron flow rates through cytochrome b₆f, plastocyanin and photosystem I were measured, revealing competition between the proton-uncoupled cyclic electron flow from PSI to Cyt b₆f to P700⁺ and the proton-coupled linear flow from PSII to Cyt b₆f to P700⁺. Analyzing responses of O₂ evolution and Chl fluorescence to ms-length light pulses we concluded that explanation of the sigmoidal fluorescence induction by excitonic connectivity between PSII units is a misconception. Each PSII processes excitation from its own antenna, but the sigmoidicity is caused by rise of the fluorescence yield of the Q_A-reduced PSII units after their Q_B site becomes occupied by reduced plastoquinone (or diuron). Unlike respiration, photosynthetic electrons must prepare their acceptor by coupled synthesis of 3ATP/4e^-. Feedback regulation of this ratio leads to oscillations under saturating light and CO₂, when the rate is P_i-limited. The slow oscillations (period 60s) indicate that the magnitudes of the deflections in the 3ATP/4e^- ratio, corrected by regulating cyclic and alternative electron flow (including the Mehler type O₂ reduction), are only a fraction of a per cent. The P_i limitation causes slip in the ATP synthase, slightly increasing the basic 12H⁺/3ATP requirement.

Potential metabolic mechanisms for inhibited chloroplast nitrogen assimilation under high CO2

Improving photosynthesis is considered a major and feasible option to dramatically increase crop yield potential. Increased atmospheric CO2 concentration often stimulates both photosynthesis and crop yield, but decreases protein content in the main C3 cereal crops. This decreased protein content in crops constrains the benefits of elevated CO2 on crop yield and affects their nutritional value for humans. To support studies of photosynthetic nitrogen assimilation and its complex interaction with photosynthetic carbon metabolism for crop improvement, we developed a dynamic systems model of plant primary metabolism, which includes the Calvin-Benson cycle, the photorespiration pathway, starch synthesis, glycolysis-gluconeogenesis, the tricarboxylic acid cycle, and chloroplastic nitrogen assimilation. This model successfully captures responses of net photosynthetic CO2 uptake rate (A), respiration rate, and nitrogen assimilation rate to different irradiance and CO2 levels. We then used this model to predict inhibition of nitrogen assimilation under elevated CO2. The potential mechanisms underlying inhibited nitrogen assimilation under elevated CO2 were further explored with this model. Simulations suggest that enhancing the supply of α-ketoglutarate is a potential strategy to maintain high rates of nitrogen assimilation under elevated CO2. This model can be used as a heuristic tool to support research on interactions between photosynthesis, respiration, and nitrogen assimilation. It also provides a basic framework to support the design and engineering of C3 plant primary metabolism for enhanced photosynthetic efficiency and nitrogen assimilation in the coming high-CO2 world.

Computational Analysis of Alternative Photosynthetic Electron Flows Linked With Oxidative Stress

During photosynthesis, organisms respond to their energy demand and ensure the supply of energy and redox equivalents that sustain metabolism. Hence, the photosynthetic apparatus can, and in fact should, be treated as an integrated supply-demand system. Any imbalance in the energy produced and consumed can lead to adverse reactions, such as the production of reactive oxygen species (ROS). Reaction centres of both photosystems are known sites of ROS production. Here, we investigate in particular the central role of Photosystem I (PSI) in this tightly regulated system. Using a computational approach we have expanded a previously published mechanistic model of C3 photosynthesis by including ROS producing and scavenging reactions around PSI. These include two water to water reactions mediated by Plastid terminal oxidase (PTOX) and Mehler and the ascorbate-glutathione (ASC-GSH) cycle, as a main non-enzymatic antioxidant. We have used this model to predict flux distributions through alternative electron pathways under various environmental stress conditions by systematically varying light intensity and enzymatic activity of key reactions. In particular, we studied the link between ROS formation and activation of pathways around PSI as potential scavenging mechanisms. This work shines light on the role of alternative electron pathways in photosynthetic acclimation and investigates the effect of environmental perturbations on PSI activity in the context of metabolic productivity.

Development of a minimized model structure and a feedback control framework for regulating photosynthetic activities

In this work, the main activities of the plant photosynthesis process are discussed to yield a minimized mathematical model structure with photosystem II (PSII) chlorophyll a fluorescence (ChlF) as a measurable output. After experimental validation of the model structure, we demonstrate that the states of the photosynthetic process may be observed by using this model and the extended Kalman filter method. We then show a feedback control framework that can be used to alter a given photosynthetic activity. The control framework is demonstrated with an example in which PSII ChlF is used as the feedback signal and light intensity is used as a controllable process input to regulate plastoquinone reduction. Although there are caveats, and further research is needed, the results lay the groundwork for further research on novel methods for optimization and regulation of photosynthetic activities, with a goal for sustainability.

Model quantification of the light-induced thylakoid membrane processes in Synechocystis sp. PCC 6803 in vivo and after exposure to radioactive irradiation

Measurements of OJIP-SMT patterns of fluorescence induction (FI) in Synechocystis sp. PCC 6803 (Synechocystis) cells on a time scale up to several minutes were mathematically treated within the framework of thylakoid membrane (T-M) model (Belyaeva et al., Photosynth Res 140:1-19, 2019) that was renewed to account for the state transitions effects. Principles of describing electron transfer in reaction centers of photosystems II and I (PSII and PSI) and cytochrome b₆f complex remained unchanged, whereas parameters for dissipative reactions of non-radiative charge recombination were altered depending on the oxidation state of Q_B-site (neutral, reduced by one electron, empty, reduced by two electrons). According to our calculations, the initial content of plastoquinol (PQH₂) in the total quinone pool of Synechocystis cells adapted to darkness for 10 min ranged between 20 and 40%. The results imply that the PQ pool mediates photosynthetic and respiratory charge flows. The redistribution of PBS antenna units responsible for the increase of Chl fluorescence in cyanobacteria (qT_{2 → 1}) upon state 2 → 1 transition or the fluorescence lowering (qT_{1 → 2}) due to state 1 → 2 transition were described in the model by exponential functions. Parameters of dynamically changed effective cross section were found by means of simulations of OJIP-SMT patterns observed on Synechocystis cells upon strong (3000 μmol photons m^-2s^-1) and moderate (1000 μmol photons m^-2s^-1) actinic light intensities. The corresponding light constant values kL_ΣAnt = 1.2 ms^-1 and 0.4 ms^-1 define the excitation of total antenna pool dynamically redistributed between PSII and PSI reaction centers. Although the OCP-induced quenching of antenna excitation is not involved in the model, the main features of the induction signals have been satisfactorily explained. In the case of strong illumination, the effective cross section decreases by approximately 33% for irradiated Synechocystis cells as compared to untreated cells. Under moderate light, the irradiated Synechocystis cells showed in simulations the same cross section as the untreated cells. The thylakoid model renewed with state transitions description allowed simulation of fluorescence induction OJIP-SMT curves detected on time scale from microseconds to minutes.

A cyanobacterial photorespiratory bypass model to enhance photosynthesis by rerouting photorespiratory pathway in C3 plants

Plants employ photosynthesis to produce sugars for supporting their growth. During photosynthesis, an enzyme Ribulose 1,5 bisphosphate carboxylase/oxygenase (Rubisco) combines its substrate Ribulose 1,5 bisphosphate (RuBP) with CO2 to produce phosphoglycerate (PGA). Alongside, Rubisco also takes up O2 and produce 2-phosphoglycolate (2-PG), a toxic compound broken down into PGA through photorespiration. Photorespiration is not only a resource-demanding process but also results in CO2 loss which affects photosynthetic efficiency in C3 plants. Here, we propose to circumvent photorespiration by adopting the cyanobacterial glycolate decarboxylation pathway into C3 plants. For that, we have integrated the cyanobacterial glycolate decarboxylation pathway into a kinetic model of C3 photosynthetic pathway to evaluate its impact on photosynthesis and photorespiration. Our results show that the cyanobacterial glycolate decarboxylation bypass model exhibits a 10% increase in net photosynthetic rate (A) in comparison with C3 model. Moreover, an increased supply of intercellular CO2 (Ci) from the bypass resulted in a 54.8% increase in PGA while reducing photorespiratory intermediates including glycolate (- 49%) and serine (- 32%). The bypass model, at default conditions, also elucidated a decline in phosphate-based metabolites including RuBP (- 61.3%). The C3 model at elevated level of inorganic phosphate (Pi), exhibited a significant change in RuBP (+ 355%) and PGA (- 98%) which is attributable to the low availability of Ci. Whereas, at elevated Pi, the bypass model exhibited an increase of 73.1% and 33.9% in PGA and RuBP, respectively. Therefore, we deduce a synergistic effect of elevation in CO2 and Pi pool on photosynthesis. We also evaluated the integrative action of CO2, Pi, and Rubisco carboxylation activity (Vcmax) on A and observed that their simultaneous increase raised A by 26%, in the bypass model. Taken together, the study potentiates engineering of cyanobacterial decarboxylation pathway in C3 plants to bypass photorespiration thereby increasing the overall efficiency of photosynthesis.

Photosynthesis: basics, history and modelling

BACKGROUND

With limited agricultural land and increasing human population, it is essential to enhance overall photosynthesis and thus productivity. Oxygenic photosynthesis begins with light absorption, followed by excitation energy transfer to the reaction centres, primary photochemistry, electron and proton transport, NADPH and ATP synthesis, and then CO2 fixation (Calvin-Benson cycle, as well as Hatch-Slack cycle). Here we cover some of the discoveries related to this process, such as the existence of two light reactions and two photosystems connected by an electron transport 'chain' (the Z-scheme), chemiosmotic hypothesis for ATP synthesis, water oxidation clock for oxygen evolution, steps for carbon fixation, and finally the diverse mechanisms of regulatory processes, such as 'state transitions' and 'non-photochemical quenching' of the excited state of chlorophyll a.

SCOPE

In this review, we emphasize that mathematical modelling is a highly valuable tool in understanding and making predictions regarding photosynthesis. Different mathematical models have been used to examine current theories on diverse photosynthetic processes; these have been validated through simulation(s) of available experimental data, such as chlorophyll a fluorescence induction, measured with fluorometers using continuous (or modulated) exciting light, and absorbance changes at 820 nm (ΔA820) related to redox changes in P700, the reaction centre of photosystem I.

CONCLUSIONS

We highlight here the important role of modelling in deciphering and untangling complex photosynthesis processes taking place simultaneously, as well as in predicting possible ways to obtain higher biomass and productivity in plants, algae and cyanobacteria.

Mechanistic modeling of light-induced chemotactic infiltration of bacteria into leaf stomata

Light is one of the factors that can play a role in bacterial infiltration into leafy greens by keeping stomata open and providing photosynthetic products for microorganisms. We model chemotactic transport of bacteria within a leaf tissue in response to photosynthesis occurring within plant mesophyll. The model includes transport of carbon dioxide, oxygen, bicarbonate, sucrose/glucose, bacteria, and autoinducer-2 within the leaf tissue. Biological processes of carbon fixation in chloroplasts, and respiration in mitochondria of the plant cells, as well as motility, chemotaxis, nutrient consumption and communication in the bacterial community are considered. We show that presence of light is enough to boost bacterial chemotaxis through the stomatal opening and toward photosynthetic products within the leaf tissue. Bacterial chemotactic ability is a major player in infiltration, and plant stomatal defense in closing the stomata as a perception of microbe-associated molecular patterns is an effective way to inhibit the infiltration.

Exposure to light can trigger photosynthesis in plant leaves, such as leafy-greens, and increase concentrations of photosynthetic products, such as glucose, within the leaf tissue. Bacteria existing at the leaf surfaces may respond to the available photosynthetic products and migrate into the leaf tissue by chemotaxis toward nutrient concentration gradients. Once the bacteria are inside the leaf tissue, they cannot be washed away, presenting a risk to the consumer. Here, a physics-based model for this light-driven infiltration is presented. This mechanistic model couples transport of bacteria and nutrients, and photosynthesis within a leaf tissue around one stomatal opening. The model shows that the ability of bacteria to transport via chemotaxis is a major factor in infiltration. A moderate intensity light is sufficient to promote chemotactic infiltration of bacteria on a leaf surface into its interior. Infiltration is enhanced in the presence of blue, white and red lights, and for a larger stomatal aperture.

Modeling Plant Metabolism: From Network Reconstruction to Mechanistic Models

Mathematical modeling of plant metabolism enables the plant science community to understand the organization of plant metabolism, obtain quantitative insights into metabolic functions, and derive engineering strategies for manipulation of metabolism. Among the various modeling approaches, metabolic pathway analysis can dissect the basic functional modes of subsections of core metabolism, such as photorespiration, and reveal how classical definitions of metabolic pathways have overlapping functionality. In the many studies using constraint-based modeling in plants, numerous computational tools are currently available to analyze large-scale and genome-scale metabolic networks. For ¹³C-metabolic flux analysis, principles of isotopic steady state have been used to study heterotrophic plant tissues, while nonstationary isotope labeling approaches are amenable to the study of photoautotrophic and secondary metabolism. Enzyme kinetic models explore pathways in mechanistic detail, and we discuss different approaches to determine or estimate kinetic parameters. In this review, we describe recent advances and challenges in modeling plant metabolism.

Relationship between irradiance and levels of Calvin-Benson cycle and other intermediates in the model eudicot Arabidopsis and the model monocot rice

Metabolite profiles provide a top-down overview of the balance between the reactions in a pathway. We compared Calvin-Benson cycle (CBC) intermediate profiles in different conditions in Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa) to learn which features of CBC regulation differ and which are shared between these model eudicot and monocot C3 species. Principal component analysis revealed that CBC intermediate profiles follow different trajectories in Arabidopsis and rice as irradiance increases. The balance between subprocesses or reactions differed, with 3-phosphoglycerate reduction being favoured in Arabidopsis and ribulose 1,5-bisphosphate regeneration in rice, and sedoheptulose-1,7-bisphosphatase being favoured in Arabidopsis compared with fructose-1,6-bisphosphatase in rice. Photosynthesis rates rose in parallel with ribulose 1,5-bisphosphate levels in Arabidopsis, but not in rice. Nevertheless, some responses were shared between Arabidopsis and rice. Fructose 1,6-bisphosphate and sedoheptulose-1,7-bisphosphate were high or peaked at very low irradiance in both species. Incomplete activation of fructose-1,6-bisphosphatase and sedoheptulose-1,7-bisphosphatase may prevent wasteful futile cycles in low irradiance. End-product synthesis is inhibited and high levels of CBC intermediates are maintained in low light or in low CO2 in both species. This may improve photosynthetic efficiency in fluctuating irradiance, and facilitate rapid CBC flux to support photorespiration and energy dissipation in low CO2.

Balancing energy supply during photosynthesis - a theoretical perspective

The photosynthetic electron transport chain (PETC) provides energy and redox equivalents for carbon fixation by the Calvin-Benson-Bassham (CBB) cycle. Both of these processes have been thoroughly investigated and the underlying molecular mechanisms are well known. However, it is far from understood by which mechanisms it is ensured that energy and redox supply by photosynthesis matches the demand of the downstream processes. Here, we deliver a theoretical analysis to quantitatively study the supply-demand regulation in photosynthesis. For this, we connect two previously developed models, one describing the PETC, originally developed to study non-photochemical quenching, and one providing a dynamic description of the photosynthetic carbon fixation in C3 plants, the CBB Cycle. The merged model explains how a tight regulation of supply and demand reactions leads to efficient carbon fixation. The model further illustrates that a stand-by mode is necessary in the dark to ensure that the carbon fixation cycle can be restarted after dark-light transitions, and it supports hypotheses, which reactions are responsible to generate such mode in vivo.

Analyzing both the fast and the slow phases of chlorophyll a fluorescence and P700 absorbance changes in dark-adapted and preilluminated pea leaves using a Thylakoid Membrane model

The dark-to-light transitions enable energization of the thylakoid membrane (TM), which is reflected in fast and slow (OJIPSMT or OABCDE) stages of fluorescence induction (FI) and P700 oxidoreduction changes (ΔA810). A Thylakoid Membrane model (T-M model), in which special emphasis has been placed on ferredoxin-NADP+-oxidoreductase (FNR) activation and energy-dependent qE quenching, was applied for quantifying the kinetics of FI and ΔA810. Pea leaves were kept in darkness for 15 min and then the FI and ΔA810 signals were measured upon actinic illumination, applied either directly or after a 10-s light pulse coupled with a subsequent 10-s dark interval. On the time scale from 40 µs to 30 s, the parallel T-M model fittings to both FI and ΔA810 signals were obtained. The parameters of FNR activation and the buildup of qE quenching were found to differ for dark-adapted and preilluminated leaves. At the onset of actinic light, photosystem II (PSII) acceptors were oxidized (neutral) after dark adaptation, while the redox states with closed and/or semiquinone QA(-)QB(-) forms were supposedly generated after preillumination, and did not relax within the 10 s dark interval. In qE simulations, a pH-dependent Hill relationship was used. The rate constant of heat losses in PSII antenna kD(t) was found to increase from the basic value kDconst, at the onset of illumination, to its maximal level kDvar due to lumenal acidification. In dark-adapted leaves, a low value of kDconst of ∼ 2 × 106 s-1 was found. Simulations on the microsecond to 30 s time scale revealed that the slow P-S-M-T phases of the fluorescence induction were sensitive to light-induced FNR activation and high-energy qE quenching. Thus, the corresponding time-dependent rate constants kD(t) and kFNR(t) change substantially upon the release of electron transport on the acceptor side of PSI and during the NPQ development. The transitions between the cyclic and linear electron transport modes have also been quantified in this paper.

Toward Multiscale Models of Cyanobacterial Growth: A Modular Approach

Oxygenic photosynthesis dominates global primary productivity ever since its evolution more than three billion years ago. While many aspects of phototrophic growth are well understood, it remains a considerable challenge to elucidate the manifold dependencies and interconnections between the diverse cellular processes that together facilitate the synthesis of new cells. Phototrophic growth involves the coordinated action of several layers of cellular functioning, ranging from the photosynthetic light reactions and the electron transport chain, to carbon-concentrating mechanisms and the assimilation of inorganic carbon. It requires the synthesis of new building blocks by cellular metabolism, protection against excessive light, as well as diurnal regulation by a circadian clock and the orchestration of gene expression and cell division. Computational modeling allows us to quantitatively describe these cellular functions and processes relevant for phototrophic growth. As yet, however, computational models are mostly confined to the inner workings of individual cellular processes, rather than describing the manifold interactions between them in the context of a living cell. Using cyanobacteria as model organisms, this contribution seeks to summarize existing computational models that are relevant to describe phototrophic growth and seeks to outline their interactions and dependencies. Our ultimate aim is to understand cellular functioning and growth as the outcome of a coordinated operation of diverse yet interconnected cellular processes.

Thylakoid membrane model of the Chl a fluorescence transient and P700 induction kinetics in plant leaves

A new Thylakoid model is presented, which describes in detail the electron/proton transfer reactions between membrane protein complexes including photosystems II and I (PSII, PSI), cytochrome (Cyt) b 6 f, mobile plastoquinone PQ pool in the thylakoid membrane, plastocyanin in lumen and ferredoxin in stroma, reduction of NADP via FNR and cyclic electron transfer. The Thylakoid model parameters were fitted both to Chl fluorescence induction data (FI) and oxido-reductions of P700 (ΔA 810) measured from 20 μs up to 20 s in pea leaves. The two-wave kinetics of FI and ΔA 810 (O(JI)PSM and OABCDE) were described quantitatively, provided that the values of membrane electrochemical potential components ΔΨ(t), pHL(t)/pHS(t) are in physiologically relevant ranges. The time courses on the time scale from nanoseconds to tens of seconds of oxido-reduction changes of ET components as well as concentrations of proton/ions (K+, Cl-) were calculated. We assume a low constant FNR activity over this period. Charge movements across the thylakoid membrane by passive leakage and active ATPase transport and proton buffer reactions are simulated. The dynamics of charge fluxes during photosynthetic induction under low light (PFD 200 μmol photons m-2 s-1) were analyzed. The initial wave of P700 oxidation within 20 ms during independent operation of PSI and PSII was followed after 50 ms by PSI donor-side reduction from reduced PQ pool via Cyt b 6 f site. The Cyt b 6 f reactions contribute to the stabilization of fluxes in the time range 1 s < t < 10 s. The detailed analysis of Chl a fluorescence at the PSM stage (t > 10 s) would need the investigation of FNR activation effect in order to explain the transitions between cyclic and linear electron transport.

The slow phase of chlorophyll a fluorescence induction in silico: Origin of the S-M fluorescence rise

In higher plants, algae, and cyanobacteria, chlorophyll (Chl) a fluorescence induction (ChlFI) has a fast (under a second) increasing OJIP phase and a slow (few minutes) PS(M)T phase, where O is for origin, the minimum fluorescence, J and I for intermediate levels, P for peak, S for a semi-steady state, M for a maximum (which is sometimes missing), and T for the terminal steady-state level. We have used a photosynthesis model of Ebenhöh et al. (Philos Trans R Soc B, 2014, doi: 10.1098/rstb.2013.0223 ) in an attempt to simulate the slow PS(M)T phase and to determine the origin of the S-M rise in Chlamydomonas (C.) reinhardtii cells. Our experiments in silico show that a slow fluorescence S-M rise (as that observed, e.g., by Kodru et al. (Photosynth Res 125:219-231, 2015) can be simulated only if the photosynthetic samples are initially in a so-called "state 2," when the absorption cross section (CS) of Photosystem II (PSII) is lower than that of PSI, and Chl a fluorescence is low (see, e.g., a review by Papageorgiou and Govindjee (J Photochem Photobiol B 104:258-270, 2011). In this case, simulations show that illumination induces a state 2 (s2) to state 1 (s1) transition (qT₂₁), and a slow S-M rise in the simulated ChlFI curve, since the fluorescence yield is known to be higher in s1, when CS of PSII is larger than that of PSI. Additionally, we have analyzed how light intensity and several photosynthetic processes influence the degree of this qT₂₁, and thus the relative amplitude of the simulated S-M phase. A refinement of the photosynthesis model is, however, necessary in order to obtain a better fit of the simulation data with the measured ChlFI curves.

The influence of leaf anatomy on the internal light environment and photosynthetic electron transport rate: exploration with a new leaf ray tracing model

Leaf photosynthesis is determined by biochemical properties and anatomical features. Here we developed a three-dimensional leaf model that can be used to evaluate the internal light environment of a leaf and its implications for whole-leaf electron transport rates (J). This model includes (i) the basic components of a leaf, such as the epidermis, palisade and spongy tissues, as well as the physical dimensions and arrangements of cell walls, vacuoles and chloroplasts; and (ii) an efficient forward ray-tracing algorithm, predicting the internal light environment for light of wavelengths between 400 and 2500nm. We studied the influence of leaf anatomy and ambient light on internal light conditions and J The results show that (i) different chloroplasts can experience drastically different light conditions, even when they are located at the same distance from the leaf surface; (ii) bundle sheath extensions, which are strips of parenchyma, collenchyma or sclerenchyma cells connecting the vascular bundles with the epidermis, can influence photosynthetic light-use efficiency of leaves; and (iii) chloroplast positioning can also influence the light-use efficiency of leaves. Mechanisms underlying leaf internal light heterogeneity and implications of the heterogeneity for photoprotection and for the convexity of the light response curves are discussed.

A mathematical model of non-photochemical quenching to study short-term light memory in plants

Plants are permanently exposed to rapidly changing environments, therefore it is evident that they had to evolve mechanisms enabling them to dynamically adapt to such fluctuations. Here we study how plants can be trained to enhance their photoprotection and elaborate on the concept of the short-term illumination memory in Arabidopsis thaliana. By monitoring fluorescence emission dynamics we systematically observe the extent of non-photochemical quenching (NPQ) after previous light exposure to recognise and quantify the memory effect. We propose a simplified mathematical model of photosynthesis that includes the key components required for NPQ activation, which allows us to quantify the contribution to photoprotection by those components. Due to its reduced complexity, our model can be easily applied to study similar behavioural changes in other species, which we demonstrate by adapting it to the shadow-tolerant plant Epipremnum aureum. Our results indicate that a basic mechanism of short-term light memory is preserved. The slow component, accumulation of zeaxanthin, accounts for the amount of memory remaining after relaxation in darkness, while the fast one, antenna protonation, increases quenching efficiency. With our combined theoretical and experimental approach we provide a unifying framework describing common principles of key photoprotective mechanisms across species in general, mathematical terms.

Different strategies of metabolic regulation in cyanobacteria: from transcriptional to biochemical control

Cyanobacteria Synechococcus sp. PCC 7942 and Synechocystis sp. PCC 6803 show similar changes in the metabolic response to changed CO₂ conditions but exhibit significant differences at the transcriptomic level. This study employs a systems biology approach to investigate the difference in metabolic regulation of Synechococcus sp. PCC 7942 and Synechocystis sp. PCC 6803. Presented multi-level kinetic model for Synechocystis sp. PCC 6803 is a new approach integrating and analysing metabolomic, transcriptomic and fluxomics data obtained under high and ambient CO₂ levels. Modelling analysis revealed that higher number of different isozymes in Synechocystis 6803 improves homeostatic stability of several metabolites, especially 3PGA by 275%, against changes in gene expression, compared to Synechococcus sp. PCC 7942. Furthermore, both cyanobacteria have the same amount of phosphoglycerate mutases but Synechocystis 6803 exhibits only ~20% differences in their mRNA levels after shifts from high to ambient CO₂ level, in comparison to ~500% differences in the case of Synechococcus sp. PCC 7942. These and other data imply that the biochemical control dominates over transcriptional regulation in Synechocystis 6803 to acclimate central carbon metabolism in the environment of variable inorganic carbon availability without extra cost carried by large changes in the proteome.

Systems approach to excitation-energy and electron transfer reaction networks in photosystem II complex: model studies for chlorophyll a fluorescence induction kinetics

Photosystem II (PS II) is a protein complex which evolves oxygen and drives charge separation for photosynthesis employing electron and excitation-energy transfer processes over a wide timescale range from picoseconds to milliseconds. While the fluorescence emitted by the antenna pigments of this complex is known as an important indicator of the activity of photosynthesis, its interpretation was difficult because of the complexity of PS II. In this study, an extensive kinetic model which describes the complex and multi-timescale characteristics of PS II is analyzed through the use of the hierarchical coarse-graining method proposed in the authors׳ earlier work. In this coarse-grained analysis, the reaction center (RC) is described by two states, open and closed RCs, both of which consist of oxidized and neutral special pairs being in quasi-equilibrium states. Besides, the PS II model at millisecond scale with three-state RC, which was studied previously, could be derived by suitably adjusting the kinetic parameters of electron transfer between tyrosine and RC. Our novel coarse-grained model of PS II can appropriately explain the light-intensity dependent change of the characteristic patterns of fluorescence induction kinetics from O-J-I-P, which shows two inflection points, J and I, between initial point O and peak point P, to O-J-D-I-P, which shows a dip D between J and I inflection points.

Modeling of the redox state dynamics in photosystem II of Chlorella pyrenoidosa Chick cells and leaves of spinach and Arabidopsis thaliana from single flash-induced fluorescence quantum yield changes on the 100 ns-10 s time scale

The time courses of the photosystem II (PSII) redox states were analyzed with a model scheme supposing a fraction of 11-25 % semiquinone (with reduced [Formula: see text]) RCs in the dark. Patterns of single flash-induced transient fluorescence yield (SFITFY) measured for leaves (spinach and Arabidopsis (A.) thaliana) and the thermophilic alga Chlorella (C.) pyrenoidosa Chick (Steffen et al. Biochemistry 44:3123-3132, 2005; Belyaeva et al. Photosynth Res 98:105-119, 2008, Plant Physiol Biochem 77:49-59, 2014) were fitted with the PSII model. The simulations show that at high-light conditions the flash generated triplet carotenoid (3)Car(t) population is the main NPQ regulator decaying in the time interval of 6-8 μs. So the SFITFY increase up to the maximum level [Formula: see text]/F 0 (at ~50 μs) depends mainly on the flash energy. Transient electron redistributions on the RC redox cofactors were displayed to explain the SFITFY measured by weak light pulses during the PSII relaxation by electron transfer (ET) steps and coupled proton transfer on both the donor and the acceptor side of the PSII. The contribution of non-radiative charge recombination was taken into account. Analytical expressions for the laser flash, the (3)Car(t) decay and the work of the water-oxidizing complex (WOC) were used to improve the modeled P680(+) reduction by YZ in the state S 1 of the WOC. All parameter values were compared between spinach, A. thaliana leaves and C. pyrenoidosa alga cells and at different laser flash energies. ET from [Formula: see text] slower in alga as compared to leaf samples was elucidated by the dynamics of [Formula: see text] fractions to fit SFITFY data. Low membrane energization after the 10 ns single turnover flash was modeled: the ∆Ψ(t) amplitude (20 mV) is found to be about 5-fold smaller than under the continuous light induction; the time-independent lumen pHL, stroma pHS are fitted close to dark estimates. Depending on the flash energy used at 1.4, 4, 100 % the pHS in stroma is fitted to 7.3, 7.4, and 7.7, respectively. The biggest ∆pH difference between stroma and lumen was found to be 1.2, thus pH- dependent NPQ was not considered.

Modeling chlorophyll a fluorescence transient: relation to photosynthesis

To honor Academician Alexander Abramovitch Krasnovsky, we present here an educational review on the relation of chlorophyll a fluorescence transient to various processes in photosynthesis. The initial event in oxygenic photosynthesis is light absorption by chlorophylls (Chls), carotenoids, and, in some cases, phycobilins; these pigments form the antenna. Most of the energy is transferred to reaction centers where it is used for charge separation. The small part of energy that is not used in photochemistry is dissipated as heat or re-emitted as fluorescence. When a photosynthetic sample is transferred from dark to light, Chl a fluorescence (ChlF) intensity shows characteristic changes in time called fluorescence transient, the OJIPSMT transient, where O (the origin) is for the first measured minimum fluorescence level; J and I for intermediate inflections; P for peak; S for semi-steady state level; M for maximum; and T for terminal steady state level. This transient is a real signature of photosynthesis, since diverse events can be related to it, such as: changes in redox states of components of the linear electron transport flow, involvement of alternative electron routes, the build-up of a transmembrane pH gradient and membrane potential, activation of different nonphotochemical quenching processes, activation of the Calvin-Benson cycle, and other processes. In this review, we present our views on how different segments of the OJIPSMT transient are influenced by various photosynthetic processes, and discuss a number of studies involving mathematical modeling and simulation of the ChlF transient. A special emphasis is given to the slower PSMT phase, for which many studies have been recently published, but they are less known than on the faster OJIP phase.

The benefits of photorespiratory bypasses: how can they work?

Bypassing the photorespiratory pathway is regarded as a way to increase carbon assimilation and, correspondingly, biomass production in C3 crops. Here, the benefits of three published photorespiratory bypass strategies are systemically explored using a systems-modeling approach. Our analysis shows that full decarboxylation of glycolate during photorespiration would decrease photosynthesis, because a large amount of the released CO2 escapes back to the atmosphere. Furthermore, we show that photosynthesis can be enhanced by lowering the energy demands of photorespiration and by relocating photorespiratory CO2 release into the chloroplasts. The conductance of the chloroplast membranes to CO2 is a key feature determining the benefit of the relocation of photorespiratory CO2 release. Although our results indicate that the benefit of photorespiratory bypasses can be improved by increasing sedoheptulose bisphosphatase activity and/or increasing the flux through the bypass, the effectiveness of such approaches depends on the complex regulation between photorespiration and other metabolic pathways.

Hierarchical modularization of biochemical pathways using fuzzy-c means clustering

Biological systems that are representative of regulatory, metabolic, or signaling pathways can be highly complex. Mathematical models that describe such systems inherit this complexity. As a result, these models can often fail to provide a path toward the intuitive comprehension of these systems. More coarse information that allows a perceptive insight of the system is sometimes needed in combination with the model to understand control hierarchies or lower level functional relationships. In this paper, we present a method to identify relationships between components of dynamic models of biochemical pathways that reside in different functional groups. We find primary relationships and secondary relationships. The secondary relationships reveal connections that are present in the system, which current techniques that only identify primary relationships are unable to show. We also identify how relationships between components dynamically change over time. This results in a method that provides the hierarchy of the relationships among components, which can help us to understand the low level functional structure of the system and to elucidate potential hierarchical control. As a proof of concept, we apply the algorithm to the epidermal growth factor signal transduction pathway, and to the C3 photosynthesis pathway. We identify primary relationships among components that are in agreement with previous computational decomposition studies, and identify secondary relationships that uncover connections among components that current computational approaches were unable to reveal.

Computer modeling of electron and proton transport in chloroplasts

Photosynthesis is one of the most important biological processes in biosphere, which provides production of organic substances from atmospheric CO2 and water at expense of solar energy. In this review, we contemplate computer models of oxygenic photosynthesis in the context of feedback regulation of photosynthetic electron transport in chloroplasts, the energy-transducing organelles of the plant cell. We start with a brief overview of electron and proton transport processes in chloroplasts coupled to ATP synthesis and consider basic regulatory mechanisms of oxygenic photosynthesis. General approaches to computer simulation of photosynthetic processes are considered, including the random walk models of plastoquinone diffusion in thylakoid membranes and deterministic approach to modeling electron transport in chloroplasts based on the mass action law. Then we focus on a kinetic model of oxygenic photosynthesis that includes key stages of the linear electron transport, alternative pathways of electron transfer around photosystem I (PSI), transmembrane proton transport and ATP synthesis in chloroplasts. This model includes different regulatory processes: pH-dependent control of the intersystem electron transport, down-regulation of photosystem II (PSII) activity (non-photochemical quenching), the light-induced activation of the Bassham-Benson-Calvin (BBC) cycle. The model correctly describes pH-dependent feedback control of electron transport in chloroplasts and adequately reproduces a variety of experimental data on induction events observed under different experimental conditions in intact chloroplasts (variations of CO2 and O2 concentrations in atmosphere), including a complex kinetics of P700 (primary electron donor in PSI) photooxidation, CO2 consumption in the BBC cycle, and photorespiration. Finally, we describe diffusion-controlled photosynthetic processes in chloroplasts within the framework of the model that takes into account complex architecture of chloroplasts and lateral heterogeneity of lamellar system of thylakoids. The lateral profiles of pH in the thylakoid lumen and in the narrow gap between grana thylakoids have been calculated under different metabolic conditions. Analyzing topological aspects of diffusion-controlled stages of electron and proton transport in chloroplasts, we conclude that along with the NPQ mechanism of attenuation of PSII activity and deceleration of PQH2 oxidation by the cytochrome b6f complex caused by the lumen acidification, the intersystem electron transport may be down-regulated due to the light-induced alkalization of the narrow partition between adjacent thylakoids of grana. The computer models of electron and proton transport described in this article may be integrated as appropriate modules into a comprehensive model of oxygenic photosynthesis.

Short-term acclimation of the photosynthetic electron transfer chain to changing light: a mathematical model

Photosynthetic eukaryotes house two photosystems with distinct light absorption spectra. Natural fluctuations in light quality and quantity can lead to unbalanced or excess excitation, compromising photosynthetic efficiency and causing photodamage. Consequently, these organisms have acquired several distinct adaptive mechanisms, collectively referred to as non-photochemical quenching (NPQ) of chlorophyll fluorescence, which modulates the organization and function of the photosynthetic apparatus. The ability to monitor NPQ processes fluorometrically has led to substantial progress in elucidating the underlying molecular mechanisms. However, the relative contribution of distinct NPQ mechanisms to variable light conditions in different photosynthetic eukaryotes remains unclear. Here, we present a mathematical model of the dynamic regulation of eukaryotic photosynthesis using ordinary differential equations. We demonstrate that, for Chlamydomonas, our model recapitulates the basic fluorescence features of short-term light acclimation known as state transitions and discuss how the model can be iteratively refined by comparison with physiological experiments to further our understanding of light acclimation in different species.

Characterization of wave phenomena in the relaxation of flash-induced chlorophyll fluorescence yield in cyanobacteria

Fluorescence yield relaxation following a light pulse was studied in various cyanobacteria under aerobic and microaerobic conditions. In Synechocystis PCC 6803 fluorescence yield decays in a monotonous fashion under aerobic conditions. However, under microaerobic conditions the decay exhibits a wave feature showing a dip at 30-50 ms after the flash followed by a transient rise, reaching maximum at ~1s, before decaying back to the initial level. The wave phenomenon can also be observed under aerobic conditions in cells preilluminated with continuous light. Illumination preconditions cells for the wave phenomenon transiently: for few seconds in Synechocystis PCC 6803, but up to one hour in Thermosynechocystis elongatus BP-1. The wave is eliminated by inhibition of plastoquinone binding either to the QB site of Photosystem-II or the Qo site of cytochrome b6f complex by 3-(3',4'-dichlorophenyl)-1,1-dimethylurea or 2,5-dibromo-3-methyl-6-isopropyl-p-benzoquinone, respectively. The wave is also absent in mutants, which lack either Photosystem-I or the NAD(P)H-quinone oxidoreductase (NDH-1) complex. Monitoring the redox state of the plastoquinone pool revealed that the dip of the fluorescence wave corresponds to transient oxidation, whereas the following rise to re-reduction of the plastoquinone pool. It is concluded that the unusual wave feature of fluorescence yield relaxation reflects transient oxidation of highly reduced plastoquinone pool by Photosystem-I followed by its re-reduction from stromal components via the NDH-1 complex, which is transmitted back to the fluorescence yield modulator primary quinone electron acceptor via charge equilibria. Potential applications of the wave phenomenon in studying photosynthetic and respiratory electron transport are discussed. This article is part of a special issue entitled: photosynthesis research for sustainability: keys to produce clean energy.

Regulatory principles and experimental approaches to the circadian control of starch turnover

In many plants, starch is synthesized during the day and degraded during the night to avoid carbohydrate starvation in darkness. The circadian clock participates in a dynamic adjustment of starch turnover to changing environmental condition through unknown mechanisms. We used mathematical modelling to explore the possible scenarios for the control of starch turnover by the molecular components of the plant circadian clock. Several classes of plausible models were capable of describing the starch dynamics observed in a range of clock mutant plants and light conditions, including discriminating circadian protocols. Three example models of these classes are studied in detail, differing in several important ways. First, the clock components directly responsible for regulating starch degradation are different in each model. Second, the intermediate species in the pathway may play either an activating or inhibiting role on starch degradation. Third, the system may include a light-dependent interaction between the clock and downstream processes. Finally, the clock may be involved in the regulation of starch synthesis. We discuss the differences among the models’ predictions for diel starch profiles and the properties of the circadian regulators. These suggest additional experiments to elucidate the pathway structure, avoid confounding results and identify the molecular components involved.

e-Photosynthesis: a comprehensive dynamic mechanistic model of C3 photosynthesis: from light capture to sucrose synthesis

Photosynthesis is arguably the most researched of all plant processes. A dynamic model of leaf photosynthesis that includes each discrete process from light capture to carbohydrate synthesis, e-photosynthesis, is described. It was developed by linking and extending our previous models of photosystem II (PSII) energy transfer and photosynthetic C3 carbon metabolism to include electron transfer processes around photosystem I (PSI), ion transfer between the lumen and stroma, ATP synthesis and NADP reduction to provide a complete representation. Different regulatory processes linking the light and dark reactions are also included: Rubisco activation via Rubisco activase, pH and xanthophyll cycle-dependent non-photochemical quenching mechanisms, as well as the regulation of enzyme activities via the ferredoxin-theoredoxin system. Although many further feedback and feedforward controls undoubtedly exist, it is shown that e-photosynthesis effectively mimics the typical kinetics of leaf CO₂ uptake, O₂ evolution, chlorophyll fluorescence emission, lumen and stromal pH, and membrane potential following perturbations in light, [CO₂] and [O₂] observed in intact C3 leaves. The model provides a framework for guiding engineering of improved photosynthetic efficiency, for evaluating multiple non-invasive measures used in emerging phenomics facilities, and for quantitative assessment of strengths and weaknesses within the understanding of photosynthesis as an integrated process.

Approximating subcellular organisation of carbohydrate metabolism during cold acclimation in different natural accessions of Arabidopsis thaliana

Accessions of Arabidopsis thaliana originating from climatically different habitats show different levels of cold acclimation when exposed to low temperatures. The central carbohydrate metabolism plays a crucial role during this acclimation. Subcellular distribution of carbohydrates over the compartments cytosol, vacuole and plastids, and putative interactions of the compartments, are analyzed in three differentially cold-tolerant accessions of Arabidopsis thaliana, originating from the Iberian Peninsula (C24), Russia (Rschew) and Scandinavia (Tenela), respectively. Subcellular carbohydrate concentrations were determined by applying the nonaqueous fractionation technique. Mathematical modeling and steady-state simulation was used to analyse the metabolic homeostasis during cold exposure. In all accessions, the initial response to cold exposure was a significant increase of plastidial and cytosolic sucrose concentrations. Raffinose accumulated in all cellular compartments of cold-tolerant accessions with a delay of 3 d, indicating that raffinose accumulation is a long-term component of cold acclimation. Minimal rates of metabolite transport permitting steady-state simulations of metabolite concentrations correlated with cold tolerance, indicating an important role of subcellular re-distribution of metabolites during cold acclimation. A highly regulated interplay of enzymatic reactions and intracellular transport processes appears to be a prerequisite for maintaining carbohydrate homeostasis during cold exposure and allowing cold acclimation in Arabidopsis thaliana.

Chlorophyll a fluorescence induction: a personal perspective of the thermal phase, the J-I-P rise

The fast (up to 1 s) chlorophyll (Chl) a fluorescence induction (FI) curve, measured under saturating continuous light, has a photochemical phase, the O-J rise, related mainly to the reduction of Q(A), the primary electron acceptor plastoquinone of Photosystem II (PSII); here, the fluorescence rise depends strongly on the number of photons absorbed. This is followed by a thermal phase, the J-I-P rise, which disappears at subfreezing temperatures. According to the mainstream interpretation of the fast FI, the variable fluorescence originates from PSII antenna, and the oxidized Q(A) is the most important quencher influencing the O-J-I-P curve. As the reaction centers of PSII are gradually closed by the photochemical reduction of Q(A), Chl fluorescence, F, rises from the O level (the minimal level) to the P level (the peak); yet, the relationship between F and [Q(A) (-)] is not linear, due to the presence of other quenchers and modifiers. Several alternative theories have been proposed, which give different interpretations of the O-J-I-P transient. The main idea in these alternative theories is that in saturating light, Q(A) is almost completely reduced already at the end of the photochemical phase O-J, but the fluorescence yield is lower than its maximum value due to the presence of either a second quencher besides Q(A), or there is an another process quenching the fluorescence; in the second quencher hypothesis, this quencher is consumed (or the process of quenching the fluorescence is reversed) during the thermal phase J-I-P. In this review, we discuss these theories. Based on our critical examination, that includes pros and cons of each theory, as well mathematical modeling, we conclude that the mainstream interpretation of the O-J-I-P transient is the most credible one, as none of the alternative ideas provide adequate explanation or experimental proof for the almost complete reduction of Q(A) at the end of the O-J phase, and for the origin of the fluorescence rise during the thermal phase. However, we suggest that some of the factors influencing the fluorescence yield that have been proposed in these newer theories, as e.g., the membrane potential ΔΨ, as suggested by Vredenberg and his associates, can potentially contribute to modulate the O-J-I-P transient in parallel with the reduction of Q(A), through changes at the PSII antenna and/or at the reaction center, or, possibly, through the control of the oxidation-reduction of the PQ-pool, including proton transfer into the lumen, as suggested by Rubin and his associates. We present in this review our personal perspective mainly on our understanding of the thermal phase, the J-I-P rise during Chl a FI in plants and algae.

Elements of a dynamic systems model of canopy photosynthesis

Improving photosynthesis throughout the full canopy rather than photosynthesis of only the top leaves of the canopy is central to improving crop yields. Many canopy photosynthesis models have been developed from physiological and ecological perspectives, however most do not consider heterogeneities of microclimatic factors inside a canopy, canopy dynamics and associated energetics, or competition among different plants, and most models lack a direct linkage to molecular processes. Here we described the rationale, elements, and approaches necessary to build a dynamic systems model of canopy photosynthesis. A systems model should integrate metabolic processes including photosynthesis, respiration, nitrogen metabolism, resource re-mobilization and photosynthate partitioning with canopy level light, CO(2), water vapor distributions and heat exchange processes. In so doing a systems-based canopy photosynthesis model will enable studies of molecular ecology and dramatically improve our insight into engineering crops for improved canopy photosynthetic CO(2) uptake, resource use efficiencies and yields.

Systematic analysis of stability patterns in plant primary metabolism

Metabolic networks are characterized by complex interactions and regulatory mechanisms between many individual components. These interactions determine whether a steady state is stable to perturbations. Structural kinetic modeling (SKM) is a framework to analyze the stability of metabolic steady states that allows the study of the system Jacobian without requiring detailed knowledge about individual rate equations. Stability criteria can be derived by generating a large number of structural kinetic models (SK-models) with randomly sampled parameter sets and evaluating the resulting Jacobian matrices. Until now, SKM experiments applied univariate tests to detect the network components with the largest influence on stability. In this work, we present an extended SKM approach relying on supervised machine learning to detect patterns of enzyme-metabolite interactions that act together in an orchestrated manner to ensure stability. We demonstrate its application on a detailed SK-model of the Calvin-Benson cycle and connected pathways. The identified stability patterns are highly complex reflecting that changes in dynamic properties depend on concerted interactions between several network components. In total, we find more patterns that reliably ensure stability than patterns ensuring instability. This shows that the design of this system is strongly targeted towards maintaining stability. We also investigate the effect of allosteric regulators revealing that the tendency to stability is significantly increased by including experimentally determined regulatory mechanisms that have not yet been integrated into existing kinetic models.

Modelling cyanobacteria: from metabolism to integrative models of phototrophic growth

Cyanobacteria are phototrophic microorganisms of global importance and have recently attracted increasing attention due to their capability to convert sunlight and atmospheric CO(2) directly into organic compounds, including carbon-based biofuels. The utilization of cyanobacteria as a biological chassis to generate third-generation biofuels would greatly benefit from an increased understanding of cyanobacterial metabolism and its interplay with other cellular processes. In this respect, metabolic modelling has been proposed as a way to overcome the traditional trial and error methodology that is often employed to introduce novel pathways. In particular, flux balance analysis and related methods have proved to be powerful tools to investigate the organization of large-scale metabolic networks-with the prospect of predicting modifications that are likely to increase the yield of a desired product and thereby to streamline the experimental progress and avoid futile avenues. This contribution seeks to describe the utilization of metabolic modelling as a research tool to understand the metabolism and phototrophic growth of cyanobacteria. The focus of the contribution is on a mathematical description of the metabolic network of Synechocystis sp. PCC 6803 and its analysis using constraint-based methods. A particular challenge is to integrate the description of the metabolic network with other cellular processes, such as the circadian clock, the photosynthetic light reactions, carbon concentration mechanism, and transcriptional regulation-aiming at a predictive model of a cyanobacterium in silico.

Fungicide impacts on photosynthesis in crop plants

Fungicides are widely used to control pests in crop plants. However, it has been reported that these pesticides may have negative effects on crop physiology, especially on photosynthesis. An alteration in photosynthesis might lead to a reduction in photoassimilate production, resulting in a decrease in both growth and yield of crop plants. For example, a contact fungicide such as copper inhibits photosynthesis by destroying chloroplasts, affecting photosystem II activity and chlorophyll biosynthesis. Systemic fungicides such as benzimidazoles, anilides, and pyrimidine are also phytotoxic, whereas azoles stimulate photosynthesis. This article focuses on the available information about toxic effects of fungicides on photosynthesis in crop plants, highlighting the mechanisms of perturbation, interaction, and the target sites of different classes of fungicides.

The size of the lumenal proton pool in leaves during induction and steady-state photosynthesis

This report describes a new method to measure the chloroplastic lumenal proton pool in leaves (tobacco and sunflower). The method is based on measurement of CO(2) outbursts from leaves caused by the shift in the CO(2) + H(2)O ↔ HCO(3)(-) + H(+) equilibrium in the chloroplast stroma as protons return from the lumen after darkening. Protons did not accumulate in the lumen to a significant extent when photosynthesis was light-limited, but a large pool of >100 μmol H(+) m(-2) accumulated in the lumen as photosynthesis became light-saturated. During thylakoid energization in the light, large amounts of protons are moved from binding sites in the stroma to binding sites in the lumen. The transthylakoidal difference in the chemical potential of free protons (ΔpH) is largely based on the difference in the chemical potential of bound protons in the lumenal and stromal compartments (pK). Over the course of the dark-light induction of photosynthesis protons accumulate in the lumen during reduction of 3-phosphoglycerate. The accumulation of electrons in reduced compounds of the stroma and cytosol is the natural cause for accumulation of a stoichiometric pool of lumenal protons during this transient event.

A quantitative comparison of Calvin-Benson cycle models

The Calvin-Benson cycle (CBC) provides the precursors for biomass synthesis necessary for plant growth. The dynamic behavior and yield of the CBC depend on the environmental conditions and regulation of the cellular state. Accurate quantitative models hold the promise of identifying the key determinants of the tightly regulated CBC function and their effects on the responses in future climates. We provide an integrative analysis of the largest compendium of existing models for photosynthetic processes. Based on the proposed ranking, our framework facilitates the discovery of best-performing models with regard to metabolomics data and of candidates for metabolic engineering.