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

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.

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.

Low light intensity elongates period and defers peak time of photosynthesis: a computational approach to circadian-clock-controlled photosynthesis in tomato

Photosynthesis is involved in the essential process of transforming light energy into chemical energy. Although the interaction between photosynthesis and the circadian clock has been confirmed, the mechanism of how light intensity affects photosynthesis through the circadian clock remains unclear. Here, we propose a first computational model for circadian-clock-controlled photosynthesis, which consists of the light-sensitive protein P, the core oscillator, photosynthetic genes, and parameters involved in the process of photosynthesis. The model parameters were determined by minimizing the cost function ( [Formula: see text]), which is defined by the errors of expression levels, periods, and phases of the clock genes (CCA1, PRR9, TOC1, ELF4, GI, and RVE8). The model recapitulates the expression pattern of the core oscillator under moderate light intensity (100 μmol m ^-2 s^-1). Further simulation validated the dynamic behaviors of the circadian clock and photosynthetic outputs under low (62.5 μmol m^-2 s^-1) and normal (187.5 μmol m^-2 s^-1) intensities. When exposed to low light intensity, the peak times of clock and photosynthetic genes were shifted backward by 1-2 hours, the period was elongated by approximately the same length, and the photosynthetic parameters attained low values and showed delayed peak times, which confirmed our model predictions. Our study reveals a potential mechanism underlying the circadian regulation of photosynthesis by the clock under different light intensities in tomato.

Molecular, Brownian, kinetic and stochastic models of the processes in photosynthetic membrane of green plants and microalgae

The paper presents the results of recent work at the Department of Biophysics of the Biological Faculty, Lomonosov Moscow State University on the kinetic and multiparticle modeling of processes in the photosynthetic membrane. The detailed kinetic models and the rule-based kinetic Monte Carlo models allow to reproduce the fluorescence induction curves and redox transformations of the photoactive pigment P700 in the time range from 100 ns to dozens of seconds and make it possible to reveal the role of individual carriers in their formation for different types of photosynthetic organisms under different illumination regimes, in the presence of inhibitors, under stress conditions. The fitting of the model curves to the experimental data quantifies the reaction rate constants that cannot be directly measured experimentally, including the non-radiative thermal relaxation reactions. We use the direct multiparticle models to explicitly describe the interactions of mobile photosynthetic carrier proteins with multienzyme complexes both in solution and in the biomembrane interior. An analysis of these models reveals the role of diffusion and electrostatic factors in the regulation of electron transport, the influence of ionic strength and pH of the cellular environment on the rate of electron transport reactions between carrier proteins. To describe the conformational intramolecular processes of formation of the final complex, in which the actual electron transfer occurs, we use the methods of molecular dynamics. The results obtained using kinetic and molecular models supplement our knowledge of the mechanisms of organization of the photosynthetic electron transport processes at the cellular and molecular levels.

A System Dynamics Approach to Model Photosynthesis at Leaf Level Under Fluctuating Light

Photosynthesis has been mainly studied under steady-state conditions even though this assumption results inadequate for assessing the biochemical responses to rapid variations occurring in natural environments. The combination of mathematical models with available data may enhance the understanding of the dynamic responses of plants to fluctuating environments and can be used to make predictions on how photosynthesis would respond to non-steady-state conditions. In this study, we present a leaf level System Dynamics photosynthesis model based and validated on an experiment performed on two soybean varieties, namely, the wild type Eiko and the chlorophyll-deficient mutant MinnGold, grown in constant and fluctuating light conditions. This mutant is known to have similar steady-state photosynthesis compared to the green wild type, but it is found to have less biomass at harvest. It has been hypothesized that this might be due to an unoptimized response to non-steady-state conditions; therefore, this mutant seems appropriate to investigate dynamic photosynthesis. The model explained well the photosynthetic responses of these two varieties to fluctuating and constant light conditions and allowed to make relevant conclusions on the different dynamic responses of the two varieties. Deviations between data and model simulations are mostly evident in the non-photochemical quenching (NPQ) dynamics due to the oversimplified combination of PsbS- and zeaxanthin-dependent kinetics, failing in finely capturing the NPQ responses at different timescales. Nevertheless, due to its simplicity, the model can provide the basis of an upscaled dynamic model at a plant level.

The role of Cytochrome b6f in the control of steady-state photosynthesis: a conceptual and quantitative model

Here, we present a conceptual and quantitative model to describe the role of the Cytochrome [Formula: see text] complex in controlling steady-state electron transport in [Formula: see text] leaves. The model is based on new experimental methods to diagnose the maximum activity of Cyt [Formula: see text] in vivo, and to identify conditions under which photosynthetic control of Cyt [Formula: see text] is active or relaxed. With these approaches, we demonstrate that Cyt [Formula: see text] controls the trade-off between the speed and efficiency of electron transport under limiting light, and functions as a metabolic switch that transfers control to carbon metabolism under saturating light. We also present evidence that the onset of photosynthetic control of Cyt [Formula: see text] occurs within milliseconds of exposure to saturating light, much more quickly than the induction of non-photochemical quenching. We propose that photosynthetic control is the primary means of photoprotection and functions to manage excitation pressure, whereas non-photochemical quenching functions to manage excitation balance. We use these findings to extend the Farquhar et al. (Planta 149:78-90, 1980) model of [Formula: see text] photosynthesis to include a mechanistic description of the electron transport system. This framework relates the light captured by PS I and PS II to the energy and mass fluxes linking the photoacts with Cyt [Formula: see text], the ATP synthase, and Rubisco. It enables quantitative interpretation of pulse-amplitude modulated fluorometry and gas-exchange measurements, providing a new basis for analyzing how the electron transport system coordinates the supply of Fd, NADPH, and ATP with the dynamic demands of carbon metabolism, how efficient use of light is achieved under limiting light, and how photoprotection is achieved under saturating light. The model is designed to support forward as well as inverse applications. It can either be used in a stand-alone mode at the leaf-level or coupled to other models that resolve finer-scale or coarser-scale phenomena.

Constructing and analysing dynamic models with modelbase v1.2.3: a software update

BACKGROUND

Computational mathematical models of biological and biomedical systems have been successfully applied to advance our understanding of various regulatory processes, metabolic fluxes, effects of drug therapies, and disease evolution and transmission. Unfortunately, despite community efforts leading to the development of SBML and the BioModels database, many published models have not been fully exploited, largely due to a lack of proper documentation or the dependence on proprietary software. To facilitate the reuse and further development of systems biology and systems medicine models, an open-source toolbox that makes the overall process of model construction more consistent, understandable, transparent, and reproducible is desired.

RESULTS AND DISCUSSION

We provide an update on the development of modelbase, a free, expandable Python package for constructing and analysing ordinary differential equation-based mathematical models of dynamic systems. It provides intuitive and unified methods to construct and solve these systems. Significantly expanded visualisation methods allow for convenient analysis of the structural and dynamic properties of models. After specifying reaction stoichiometries and rate equations modelbase can automatically assemble the associated system of differential equations. A newly provided library of common kinetic rate laws reduces the repetitiveness of the computer programming code. modelbase is also fully compatible with SBML. Previous versions provided functions for the automatic construction of networks for isotope labelling studies. Now, using user-provided label maps, modelbase v1.2.3 streamlines the expansion of classic models to their isotope-specific versions. Finally, the library of previously published models implemented in modelbase is growing continuously. Ranging from photosynthesis to tumour cell growth to viral infection evolution, all these models are now available in a transparent, reusable and unified format through modelbase.

CONCLUSION

With this new Python software package, which is written in currently one of the most popular programming languages, the user can develop new models and actively profit from the work of others. modelbase enables reproducing and replicating models in a consistent, tractable and expandable manner. Moreover, the expansion of models to their isotopic label-specific versions enables simulating label propagation, thus providing quantitative information regarding network topology and metabolic fluxes.

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.

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.

Modelling and simulation of chlorophyll fluorescence from PSII of a plant leaf as affected by both illumination light intensities and temperatures

The emission of chlorophyll fluorescence (ChlF) from photosystem II (PSII) of plant leaves the couple with photoelectron transduction cascades in photosynthetic reactions and can be used to probe photosynthetic efficiency and plant physiology. Because of population increase, food shortages, and global warming, it is becoming more and more urgent to enhance plant photosynthesis efficiency by controlling plant growth rate. An effective model structure is essential for plant control strategy development. However, there is a lack of reporting on modelling and simulation of PSII activities under the interaction of both illumination light intensities and temperatures, which are the two important controllable factors affecting, plant growth, especially for a greenhouse. In this work, the authors extended their work on modelling photosynthetic activities as affected by light and temperature to cover both the interaction effects of illumination light intensities and temperature on ChlF emission. Experiments on ChlF were performed under different light intensities and temperatures and used to validate the developed model structure. The average relative error between experimental data and model fitting is <0.3%, which shows the effectiveness of the developed model structure. Simulations were performed to show the interaction effect of light and temperature effects on photosynthetic activities.

Modelling and simulation of photosystem II chlorophyll fluorescence transition from dark-adapted state to light-adapted state

Green houses play a vital role in modern agriculture. Artificial light illumination is very important in a green house. While light is necessary for plant growth, excessive light in a green house may not bring more profit and even damages plants. Developing a plant-physiology-based light control strategy in a green house is important, which implies that a state-space model on photosynthetic activities is very useful because modern control theories and techniques are usually developed according to model structures in the state space. In this work, a simplified model structure on photosystem II activities was developed with seven state variables and chlorophyll fluorescence (ChlF) as the observable variable. Experiments on ChlF were performed. The Levenberg-Marquardt algorithm was used to estimate model parameters from experimental data. The model structure can fit experimental data with a small relative error (<2%). ChlF under different light intensities were simulated to show the effect of light intensity on ChlF emission. A simplified model structure with fewer state variables and model parameters will be more robust to perturbations and model parameter estimation. The model structure is thus expected useful in future green-house light control strategy development.

Computational modelling predicts substantial carbon assimilation gains for C3 plants with a single-celled C4 biochemical pump

Achieving global food security for the estimated 9 billion people by 2050 is a major scientific challenge. Crop productivity is fundamentally restricted by the rate of fixation of atmospheric carbon. The dedicated enzyme, RubisCO, has a low turnover and poor specificity for CO2. This limitation of C3 photosynthesis (the basic carbon-assimilation pathway present in all plants) is alleviated in some lineages by use of carbon-concentrating-mechanisms, such as the C4 cycle-a biochemical pump that concentrates CO2 near RubisCO increasing assimilation efficacy. Most crops use only C3 photosynthesis, so one promising research strategy to boost their productivity focuses on introducing a C4 cycle. The simplest proposal is to use the cycle to concentrate CO2 inside individual chloroplasts. The photosynthetic efficiency would then depend on the leakage of CO2 out of a chloroplast. We examine this proposal with a 3D spatial model of carbon and oxygen diffusion and C4 photosynthetic biochemistry inside a typical C3-plant mesophyll cell geometry. We find that the cost-efficiency of C4 photosynthesis depends on the gas permeability of the chloroplast envelope, the C4 pathway having higher quantum efficiency than C3 for permeabilities below 300 μm/s. However, at higher permeabilities the C4 pathway still provides a substantial boost to carbon assimilation with only a moderate decrease in efficiency. The gains would be capped by the ability of chloroplasts to harvest light, but even under realistic light regimes a 100% boost to carbon assimilation is possible. This could be achieved in conjunction with lower investment in chloroplasts if their cell surface coverage is also reduced. Incorporation of this C4 cycle into C3 crops could thus promote higher growth rates and better drought resistance in dry, high-sunlight climates.

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.

Simulation of chlorophyll fluorescence rise and decay kinetics, and P700-related absorbance changes by using a rule-based kinetic Monte-Carlo method

A model of primary photosynthetic reactions in the thylakoid membrane was developed and its validity was tested by simulating three types of experimental kinetic curves: (1) the light-induced chlorophyll a fluorescence rise (OJIP transients) reflecting the stepwise transition of the photosynthetic electron transport chain from the oxidized to the fully reduced state; (2) the dark relaxation of the flash-induced fluorescence yield attributed to the Q_A^- oxidation kinetics in PSII; and (3) the light-induced absorbance changes near 820 or 705 nm assigned to the redox transitions of P₇₀₀ in PSI. A model was implemented by using a rule-based kinetic Monte-Carlo method and verified by simulating experimental curves under different treatments including photosynthetic inhibitors, heat stress, anaerobic conditions, and very high light intensity.

Impact assessment of leaf pigments in selected landscape plants exposed to roadside dust

Continuous addition of undesired effluents to the environment affects foliar surface of leaf, changes their morphology, stomata, photosynthetic pigments, and biochemical constituents which result in massive damage due to persistent nature of the pollutant. In persistent hostile environment, plants fail to grow and develop, and the effects are often extensive. In current study, landscape plants were exposed to different levels of road dust to analyze the effect on various photosynthetic pigments. Dry roadside sediments were collected through a vacuum pump and passed through filters to get fine particles less than 100 μm and sprinkled on Euphorbia milii (EM), Gardenia jasminoides (GJ), and Hibiscus rosa-sinensis (HRs) by using a hand pump, twice daily at T₁ (control), T₂, T₃, and T₄ (0, 2, 4, and 6 g/plant, respectively) for a period of 3 months in green house. Road sediment significantly reduces leaf pigments in landscape plants population and the effects were more severe in high level of dust deposition. Individual response of EM, GJ, and HRs to different levels of road dust was variable; however, road sediment significantly reduces leaf pigments at high dose of roadside dust deposition. EM plants exposed to 2 g/plant roadside dust showed higher chlorophyll-a, chlorophyll-b, total chlorophyll, chlorophyllide-b, and polar carotenoid contents as compared to GJ and HRs. Leaf chlorophyll-a, chlorophyll-b, total chlorophyll, carotenoid, and polar carotenoid contents of EM were higher than GJ and HRs in T₃ and T₄ treatments. However HRs showed significantly higher protochlorophyllide, chlorophyllide-a, and pheophytin-b contents of leaf in T₄ group. EM was found as tolerant landscape plant followed by HRs. GJ was most vulnerable to road dust stress. Present study concludes that the entire biosynthesis of leaf pigments is in chain and interlinked together where effect of road dust on one pigment influences other pigments and their derivatives. Salient features of the present study provide useful evidence to estimate roadside dust as a major risk factor for plant pigments, and plants in green belt along roadside suffer retarded growth and fail to establish and develop.

Mechanisms of Photodamage and Protein Turnover in Photoinhibition

Rapid protein degradation and replacement is an important response to photodamage and a means of photoprotection by recovering proteostasis. Protein turnover and translation efficiency studies have discovered fast turnover subunits in cytochrome b₆f and the NAD(P)H dehydrogenase (NDH) complex, in addition to PSII subunit D1. Mutations of these complexes have been linked to enhanced photodamage at least partially via cyclic electron flow. Photodamage and photoprotection involving cytochrome b₆f, NDH complex, cyclic electron flow, PSI, and nonphotochemical quenching proteins have been reported. Here, we propose that the rapid turnover of specific proteins in cytochrome b₆f and the NDH complex need to be characterised and compared with the inhibition of PSII by excess excitation energy and PSI by excess electron flux to expand our understanding of photoinhibition mechanisms.

A four state parametric model for the kinetics of the non-photochemical quenching in Photosystem II

The phenomenon of non-photochemical quenching (NPQ) was studied in spinach chloroplasts using pulse amplitude modulated (PAM) fluorometry. We present a new analysis method which describes the observed fluorescence quantum yield as the sum of the product of four different states of PSII and their corresponding quantum yields. These four distinct states are PSII in the quenched or unquenched state, and with its reaction center either open or closed depending upon the reduction of the Q_A site. With this method we can describe the dynamics of the NPQ induction and recovery as well as quantify the percentage of photoinactivated RC throughout the measurement. We show that after one cycle of quenching followed by a period of recovery, approximately 8-9% of the RC are photoinactivated, after two cycles of illumination this number becomes 15-17%. The recovery from the quenching appeared with rates of (50s)^-1 and (1h)^-1. The new analysis method presented here is flexible, allowing it to be applied to any type of PAM fluorometry protocol. The method allows to quantitatively compare qualitatively different PAM curves on the basis of statistically relevant fitting parameters and to quantify quenching dynamics and photoinactivation. Moreover, the results presented here demonstrate that the analysis of a single PAM fluorometry quenching experiment can already provide information on the relative quantum yield of the four different states of PSII for the intact chloroplasts - something no other form of spectroscopy could provide in a single measurement.

A systems-wide understanding of photosynthetic acclimation in algae and higher plants

The ability of phototrophs to colonise different environments relies on robust protection against oxidative stress, a critical requirement for the successful evolutionary transition from water to land. Photosynthetic organisms have developed numerous strategies to adapt their photosynthetic apparatus to changing light conditions in order to optimise their photosynthetic yield, which is crucial for life on Earth to exist. Photosynthetic acclimation is an excellent example of the complexity of biological systems, where highly diverse processes, ranging from electron excitation over protein protonation to enzymatic processes coupling ion gradients with biosynthetic activity, interact on drastically different timescales from picoseconds to hours. Efficient functioning of the photosynthetic apparatus and its protection is paramount for efficient downstream processes, including metabolism and growth. Modern experimental techniques can be successfully integrated with theoretical and mathematical models to promote our understanding of underlying mechanisms and principles. This review aims to provide a retrospective analysis of multidisciplinary photosynthetic acclimation research carried out by members of the Marie Curie Initial Training Project, AccliPhot, placing the results in a wider context. The review also highlights the applicability of photosynthetic organisms for industry, particularly with regards to the cultivation of microalgae. It intends to demonstrate how theoretical concepts can successfully complement experimental studies broadening our knowledge of common principles in acclimation processes in photosynthetic organisms, as well as in the field of applied microalgal biotechnology.

Frequently asked questions about chlorophyll fluorescence, the sequel

Using chlorophyll (Chl) a fluorescence many aspects of the photosynthetic apparatus can be studied, both in vitro and, noninvasively, in vivo. Complementary techniques can help to interpret changes in the Chl a fluorescence kinetics. Kalaji et al. (Photosynth Res 122:121-158, 2014a) addressed several questions about instruments, methods and applications based on Chl a fluorescence. Here, additional Chl a fluorescence-related topics are discussed again in a question and answer format. Examples are the effect of connectivity on photochemical quenching, the correction of F V /F M values for PSI fluorescence, the energy partitioning concept, the interpretation of the complementary area, probing the donor side of PSII, the assignment of bands of 77 K fluorescence emission spectra to fluorescence emitters, the relationship between prompt and delayed fluorescence, potential problems when sampling tree canopies, the use of fluorescence parameters in QTL studies, the use of Chl a fluorescence in biosensor applications and the application of neural network approaches for the analysis of fluorescence measurements. The answers draw on knowledge from different Chl a fluorescence analysis domains, yielding in several cases new insights.

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.

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.

A reductionist approach to model photosynthetic self-regulation in eukaryotes in response to light

Along with the development of several large-scale methods such as mass spectrometry or micro arrays, genome wide models became not only a possibility but an obvious tool for theoretical biologists to integrate and analyse complex biological data. Nevertheless, incorporating the dynamics of photosynthesis remains one of the major challenges while reconstructing metabolic networks of plants and other photosynthetic organisms. In this review, we aim to provide arguments that small-scale models are still a suitable choice when it comes to discovering organisational principles governing the design of biological systems. We give a brief overview of recent modelling efforts in understanding the interplay between rapid, photoprotective mechanisms and the redox balance within the thylakoid membrane, discussing the applicability of a reductionist approach in modelling self-regulation in plants and outline possible directions for further research.

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.

Sharing light between two photosystems: mechanism of state transitions

In the thylakoid membrane, the two photosystems act in series to promote linear electron flow, with the concomitant production of ATP and reducing equivalents such as NADPH. Photosystem I, which is preferentially activated in far-red light, also energizes cyclic electron flow which generates only ATP. Thus, changes in light quality and cellular metabolic demand require a rapid regulation of the activity of the two photosystems. At low light intensities, this is mediated by state transitions. They allow the dynamic allocation of light harvesting antennae to the two photosystems, regulated through protein phosphorylation by a kinase and phosphatase pair that respond to the redox state of the electron transfer chain. Phosphorylation of the antennae leads to remodeling of the photosynthetic complexes.

Quantitative proteomic analysis of thylakoid from two microalgae (Haematococcus pluvialis and Dunaliella salina) reveals two different high light-responsive strategies

Under high light (HL) stress, astaxanthin-accumulating Haematococcus pluvialis and β-carotene-accumulating Dunaliella salina showed different responsive patterns. To elucidate cellular-regulating strategies photosynthetically and metabolically, thylakoid membrane proteins in H. pluvialis and D. salina were extracted and relatively quantified after 0 h, 24 h and 48 h of HL stress. Proteomic analysis showed that three subunits of the cytochrome b6/f complex were greatly reduced under HL stress in H. pluvialis, while they were increased in D. salina. Additionally, the major subunits of both photosystem (PS) II and PSI reaction center proteins were first reduced and subsequently recovered in H. pluvialis, while they were gradually reduced in D. salina. D. salina also showed a greater ability to function using the xanthophyll-cycle and the cyclic photosynthetic electron transfer pathway compared to H. pluvialis. We propose a reoriented and effective HL-responsive strategy in H. pluvialis, enabling it to acclimate under HL. The promising metabolic pathway described here contains a reorganized pentose phosphate pathway, Calvin cycle and glycolysis pathway participating in carbon sink formation under HL in H. pluvialis. Additionally, the efficient carbon reorientation strategy in H. pluvialis was verified by elevated extracellular carbon assimilation and rapid conversion into astaxanthin.

Chloroplast remodeling during state transitions in Chlamydomonas reinhardtii as revealed by noninvasive techniques in vivo

Plants respond to changes in light quality by regulating the absorption capacity of their photosystems. These short-term adaptations use redox-controlled, reversible phosphorylation of the light-harvesting complexes (LHCIIs) to regulate the relative absorption cross-section of the two photosystems (PSs), commonly referred to as state transitions. It is acknowledged that state transitions induce substantial reorganizations of the PSs. However, their consequences on the chloroplast structure are more controversial. Here, we investigate how state transitions affect the chloroplast structure and function using complementary approaches for the living cells of Chlamydomonas reinhardtii. Using small-angle neutron scattering, we found a strong periodicity of the thylakoids in state 1, with characteristic repeat distances of ∼ 200 Å, which was almost completely lost in state 2. As revealed by circular dichroism, changes in the thylakoid periodicity were paralleled by modifications in the long-range order arrangement of the photosynthetic complexes, which was reduced by ∼ 20% in state 2 compared with state 1, but was not abolished. Furthermore, absorption spectroscopy reveals that the enhancement of PSI antenna size during state 1 to state 2 transition (∼ 20%) is not commensurate to the decrease in PSII antenna size (∼ 70%), leading to the possibility that a large part of the phosphorylated LHCIIs do not bind to PSI, but instead form energetically quenched complexes, which were shown to be either associated with PSII supercomplexes or in a free form. Altogether these noninvasive in vivo approaches allow us to present a more likely scenario for state transitions that explains their molecular mechanism and physiological consequences.

Changing the light environment: chloroplast signalling and response mechanisms

Light is an essential environmental factor required for photosynthesis, but it also mediates signals to control plant development and growth and induces stress tolerance. The photosynthetic organelle (chloroplast) is a key component in the signalling and response network in plants. This theme issue of Philosophical Transactions of the Royal Society of London B: Biology provides updates, highlights and summaries of the most recent findings on chloroplast-initiated signalling cascades and responses to environmental changes, including light and biotic stress. Besides plant molecular cell biology and physiology, the theme issue includes aspects from the cross-disciplinary fields of environmental adaptation, ecology and agronomy.