Control of the rate of photosynthetic carbon dioxide fixation

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.

Modeling the Calvin-Benson cycle

BACKGROUND

Modeling the Calvin-Benson cycle has a history in the field of theoretical biology. Anyone who intends to model this system will look at existing models to adapt, refine and improve them. With the goal to study the regulation of carbon metabolism, we investigated a broad range of relevant models for their suitability to provide the basis for further modeling efforts. Beyond a critical analysis of existing models, we furthermore investigated the question how adjacent metabolic pathways, for instance photorespiration, can be integrated in such models.

RESULTS

Our analysis reveals serious problems with a range of models that are publicly available and widely used. The problems include the irreproducibility of the published results or significant differences between the equations in the published description of the model and model itself in the supplementary material. In addition to and based on the discussion of existing models, we furthermore analyzed approaches in PGA sink implementation and confirmed a weak relationship between the level of its regulation and efficiency of PGA export, in contrast to significant changes in the content of metabolic pool within the Calvin-Benson cycle.

CONCLUSIONS

In our study we show that the existing models that have been investigated are not suitable for reuse without substantial modifications. We furthermore show that the minor adjacent pathways of the carbon metabolism, neglected in all kinetic models of Calvin-Benson cycle, cannot be substituted without consequences in the mass production dynamics. We further show that photorespiration or at least its first step (O2 fixation) has to be implemented in the model if this model is aimed for analyses out of the steady state.

Molecular cloning of the Arabidopsis thaliana sedoheptulose-1,7-biphosphatase gene and expression studies in wheat and Arabidopsis thaliana

We report here the isolation and nucleotide sequence of genomic clones encoding the chloroplast enzyme sedoheptulose-1,7-bisphosphatase (SBPase) from Arabidopsis thaliana. The coding region of this gene contains eight exons (72-76 bp) and seven introns (75-91 bp) and encodes a polypeptide of 393 amino acids. Unusually, the 5' non-coding region contains two additional AUG codons upstream of the translation initiation codon. A comparison of the deduced Arabidopsis and wheat SBPase polypeptide sequences reveals 78.6%, identity. Expression studies showed that the level of SBPase mRNA in Arabidopsis and wheat is regulated in a light-dependent manner and is also influenced by the developmental stage of the leaf. Although the Arabidopsis SBPase gene is present in a single copy, two hybridizing transcripts were detected in some tissues, suggesting the presence of alternate transcription start sites in the upstream region.

A light- and developmentally-regulated DNA-binding interaction is common to the upstream sequences of the wheat Calvin cycle bisphosphatase genes

We have characterised a DNA-binding interaction common to the upstream sequences of the wheat fructose-1,6-bisphosphatase (FBPase) and sedoheptulose-1,7-bisphosphatase (SBPase) genes. The recognition site for this sequence-specific binding activity, designated wheat FBPase factor (WF-1), is located within 125 bp of the transcription start site of each gene. Within these regions there are no sequence motifs similar to those shown to be important for light-regulated expression in other species. The binding activity was not detected in wheat root nuclear extracts, or in pea leaf extracts. There was a higher level of binding activity in light-grown than in dark-grown wheat leaves. The level was also found to decline when light-grown plants were given an extended dark treatment, but could be reinduced by light. Utilising the gradient of developmental maturity which exists within the wheat leaf it was found that WF-1 activity increases during leaf development.

cDNA and gene sequences of wheat chloroplast sedoheptulose-1,7-bisphosphatase reveal homology with fructose-1,6-bisphosphatases

The nucleotide sequence encoding the chloroplast enzyme, sedoheptulose-1,7-bisphosphatase [Sed(1,7)P2ase], was obtained from wheat cDNA and genomic clones. The transcribed region of the Sed(1,7)P2ase gene has eight exons (72-507 bp) and seven introns (85-626 bp) and encodes a precursor polypeptide of 393 amino acids. Comparison of the deduced amino acid sequence of Sed(1,7)P2ase with those of fructose-1,6-bisphosphatase [Fru(1,6)P2ase] enzymes from a variety of sources reveals 19% identity, rising to 42% if conservative changes are considered. Most importantly, the amino acid residues which form the active site of Fru(1,6)P2ase are highly conserved in the Sed(1,7)P2ase molecule, indicating a common catalytic mechanism. Interestingly, although the activities of both Sed(1,7)P2ase and chloroplast Fru(1,6)P2ase are modulated by light via the thioredoxin system, the amino acid sequence motif identified as having a role in this regulation in chloroplast Fru(1,6)P2ase is not found in the Sed(1,7)P2ase enzyme.

Kinetics of the conformational transition of the spinach chloroplast fructose-1,6-bisphosphatase induced by fructose 2,6-bisphosphate

The activation of oxidized chloroplast fructose-1,6-bisphosphatase by fructose 2,6-bisphosphate and magnesium previously described at pH 7.5 [Soulié et al. (1988) Eur. J. Biochem. 176, 111-117] has now been studied at pH 8, the pH which prevails under light conditions in the chloroplast stroma. The process obeys a hysteretic mechanism but the rate of activation is considerably increased with half-times down to 50 s and the apparent dissociation constant of fructose 2,6-bisphosphate from the enzyme is lowered from 1 mM at pH 7.5 to 3.3 microM at pH 8. The process is strictly metal-dependent with a half-saturation concentration of 2.54 mM for magnesium. The conformational transition postulated in our hysteretic model has been investigated through both the spectrophometric and chemical modification approaches. The activation of the enzyme by fructose 2,6-bisphosphate in the presence of magnesium results in a slow modification of the ultraviolet absorption spectrum of the enzyme with an overall increase of 3% at 290 nm. The same treatment leads to the protection of two free sulfhydryls and an increased reactivity of one sulfhydryl group/enzyme monomer to modification by 5,5'-dithiobis(2-nitrobenzoic acid). The titration of the exposed cysteinyl residue prevents the relaxation of enzyme species induced by fructose 2,6-bisphosphate to the native form. The activation of chloroplast fructose-1,6-bisphosphatase by fructose 2,6-bisphosphate is discussed both with respect to the understanding of the overall regulation properties of the enzyme and to a possible physiological significance of this process.

Control analysis of photosynthate partitioning : Impact of reduced activity of ADP-glucose pyrophosphorylase or plastid phosphoglucomutase on the fluxes to starch and sucrose inArabidopsis thaliana (L.) Heynh

Experiments were carried out to investigate the contribution of ADP-glucose pyrophosphorylase and the plastid phosphoglucosemutase to the control of starch synthesis. Mutants ofArabidopsis thaliana (L.) Heyhn. were constructed with 50% and 7% of the wild-type adenosine 5'-diphosphoglucose pyrophosphorylase (ADPGlc-PPase), or 50% and null plastid phosphoglucomutase (PGM). The changes in the steady-state rates of sucrose synthesis, starch synthesis and CO2 fixation were measured in saturating CO2 in low (75 μmol·m(-2)·s(-1)) and high (600 μmol·m(-2)·s(-1)) irradiance. In low irradiance, a 50% decrease of PGM had no significant effect on fluxes, while a 50% and 93% decrease of ADPGlc-PPase led to a 23% and 74% inhibition of starch synthesis. Decreased ADPGlc-PPase led to an increase of hexose phosphates, triose phosphates and fructose-1,6-bisphosphate. Fixation of CO2 was not inhibited because the inhibition of starch synthesis was matched by a stimulation of sucrose synthesis. In high irradiance, a 50% decrease of PGM led to a 20% inhibition of starch synthesis. A 50% and 93% decrease of ADPGlc-PPase led to a 39% and 90% inhibition of starch synthesis. Sucrose synthesis was also inhibited, and the rate of photosynthesis was decreased. Decreased ADPGlc-PPase led to an increase of hexose phosphates but triose phosphates and fructose-1,6-bisphosphate did not increase. These results are used to estimate flux-control coefficients for these enzymes for starch synthesis. Firstly, the flux to starch is only controlled by ADPGlc-PPase in low irradiance, but control is redistributed to other enzymes in the pathway when a rapid flux is imposed, e.g. in high irradiance and CO2. Secondly, reducing the rate of starch synthesis by decreasing the activity of enzymes in this pathway does not always lead to a compensating increase in the rate of sucrose synthesis. Thirdly, decreasing the activity of an enzyme by a factor of two compared to the remainder of the pathway often leads to it exerting very considerable control. Fourthly, each enzyme starts to exert considerable control when only a fraction of its Vmax activity is being utilised in vivo, for example the maximum flux at ADPGlc-PPase never exceeded 20% of the Vmax activity. The summation theory is also applied to check whether additional major control sites are required. In low irradiance, the efficiency of light harvesting will exert considerable control over the rate of starch synthesis.

Control analysis of photosynthetic CO2 fixation

The potential of control analysis to aid our understanding of regulation and control of photosynthetic carbon metabolism is investigated. Methods of metabolic control analysis are used to determine flux control coefficients of photosynthetic reactions from enzyme elasticities. Equations expressing control coefficients symbolically by enzyme elasticities are derived, and general properties of these expressions are analysed. Suggestions for experimental determination of flux control coefficients from enzyme elasticities are given. A simplified model of the Calvin-Benson cycle is used to illustrate interrelations between patterns of photosynthetic metabolites and that of control coefficients.

Metabolic control therapy and biochemical systems theory: different objectives, different assumptions, different results

The claim by Savageau et al. (1987 a, b, Math. Biosci. 86, 127-145, 147-167) that the theory of metabolic control associated with Kacser & Burns (1973, Symp. Soc. Exp. Biol. 27, 65-104) and with Heinrich & Rapoport (1974, Eur. J. Biochem. 42, 89-102) is no more than a special case of the biochemical systems theory of Savageau and colleagues is examined. It is shown to be based on a misconception of the objectives and assumptions of metabolic control theory. In particular, the control and elasticity coefficients that play a central role in metabolic control theory are not constants and cannot be treated as constants. Consequently they cannot in general be equated with the kinetic orders that appear in biochemical systems theory, though they do correspond at the point where the two theories are tangential to one another.

The matrix method of metabolic control analysis: its validity for complex pathway structures

The sensitivities of the variables of a metabolic system (such as fluxes and concentrations) to variations in enzyme concentration are expressed in metabolic control analysis as control coefficients. The matrix method is a system of writing matrix equations that generate expressions for the control coefficients in terms of the characteristics of the components (principally the enzymes). Previously, the matrix method has been considered in terms of simple pathway structures; here we justify its applicability to complex pathways, such as those with multiple branches. It is shown that this requires modification of the branch point relationship to take account of changes of flux along the limbs of the branch and of stoichiometric factors. The method of deriving the flux control coefficients with respect to different fluxes in the system is extended to cope with these circumstances.

Fractional control of photosynthesis by the QB protein, the cytochrome f/b 6 complex and other components of the photosynthesic apparatus

In order to obtain information on fractional control of photosynthesis by individual catalysts, catalytic activities in photosynthetic electron transport and carbon metabolism were modified by the addition of inhibitors, and the effect on photosynthetic flux was measured using chloroplasts of Spinacia oleracea L. In thylakoids with coupled electron transport, light-limited electron flow to ferricyanide was largely controlled by the QB protein of the electron-transport chain. Fractional control by the cytochrome f/b 6 complex was insignificant under these conditions. Control by the cytochrome f/b 6 complex dominated at high energy fluence rates where the contribution to control of the QB protein was very small. Uncoupling shifted control from the cytochrome f/b 6 complex to the QB protein. Control of electron flow was more complex in assimilating chloroplasts than in thylakoids. The contributions of the cytochrome f/b 6 complex and of the QB protein to control were smaller in intact chloroplasts than in thylakoids. Thus, even though the transit time for an electron through the electron-transport chain may be below 5 ms in leaves, oxidation of plastohydroquinone was only partially responsible for limiting photosynthesis under conditions of light and CO2 saturation. The energy fluence rate influenced control coefficients. Fractional control of photosynthesis by the ATP synthetase, the cytochrome f/b 6 complex and by ribulose-1,5-bisphosphate carboxylase increased with increasing fluence rates, whereas the contributions of the QB protein and of enzymes sensitive to SH-blocking agents decreased. The results show that the burdens of control are borne by several components of the photosynthetic apparatus, and that burdens are shifted as conditions for photosynthesis change.

Metabolic control and its analysis. Extensions to the theory and matrix method

The matrix algebra procedure for determining the flux control coefficients of enzymes in metabolic pathways has been extended to allow determination of the concentration control coefficients. Although it is shown that the procedure is essentially unchanged in most cases, the presence of moiety-conserved cycles in a pathway places additional limitations on the form of the equations that can be used in the matrix formulation for concentration control coefficients. In the case of branched pathways, a new coefficient has been defined, the branch distribution control coefficient, which can be obtained via the matrix procedure. Thus a single matrix equation permits calculation or algebraic evaluation of the control coefficients for flux, concentration and distribution of flux at branches, so that the complete response of a pathway to alteration of enzyme content, or to modulation by an effector, can be determined. The relationships have been determined between flux control coefficients in isolated sections of metabolic pathways and the coefficients for the same enzymes when part of a larger metabolic system. It is shown that the control analysis of the isolated system provides useful information towards determining the control properties of the extended system.