Entner-Doudoroff pathway in Synechocystis PCC 6803: Proposed regulatory roles and enzyme multifunctionalities.

The Entner-Doudoroff pathway (ED-P) was established in 2016 as the fourth glycolytic pathway in Synechocystis sp. PCC 6803. ED-P consists of two reactions, the first catalyzed by 6-phosphogluconate dehydratase (EDD), the second by keto3-deoxygluconate-6-phosphate aldolase/4-hydroxy-2-oxoglutarate aldolase (EDA). ED-P was previously concluded to be a widespread (∼92%) pathway among cyanobacteria, but current bioinformatic analysis estimated the occurrence of ED-P to be either scarce (∼1%) or uncommon (∼47%), depending if dihydroxy-acid dehydratase (ilvD) also functions as EDD (currently assumed). Thus, the biochemical characterization of ilvD is a task pending to resolve this uncertainty. Next, we have provided new insights into several single and double glycolytic mutants based on kinetic model of central carbon metabolism of Synechocystis. The model predicted that silencing 6-phosphogluconate dehydrogenase (gnd) could be coupled with ∼90% down-regulation of G6P-dehydrogenase, also limiting the metabolic flux via ED-P. Furthermore, our metabolic flux estimation implied that growth impairment linked to silenced EDA under mixotrophic conditions is not caused by diminished carbon flux via ED-P but rather by a missing mechanism related to the role of EDA in metabolism. We proposed two possible, mutually non-exclusive explanations: (i) Δeda leads to disrupted carbon catabolite repression, regulated by 2-keto3-deoxygluconate-6-phosphate (ED-P intermediate), and (ii) EDA catalyzes the interconversion between glyoxylate and 4-hydroxy-2-oxoglutarate + pyruvate in the proximity of TCA cycle, possibly effecting the levels of 2-oxoglutarate under Δeda. We have also proposed a new pathway from EDA toward proline, which could explain the proline accumulation under Δeda. In addition, the presented in silico method provides an alternative to 13C metabolic flux analysis for marginal metabolic pathways around/below the threshold of ultrasensitive LC-MS. Finally, our in silico analysis provided alternative explanations for the role of ED-P in Synechocystis while identifying some severe uncertainties.

A new insight into role of phosphoketolase pathway in Synechocystis sp. PCC 6803.

Phosphoketolase (PKET) pathway is predominant in cyanobacteria (around 98%) but current opinion is that it is virtually inactive under autotrophic ambient CO2 condition (AC-auto). This creates an evolutionary paradox due to the existence of PKET pathway in obligatory photoautotrophs. We aim to answer the paradox with the aid of bioinformatic analysis along with metabolic, transcriptomic, fluxomic and mutant data integrated into a multi-level kinetic model. We discussed the problems linked to neglected isozyme, pket2 (sll0529) and inconsistencies towards the explanation of residual flux via PKET pathway in the case of silenced pket1 (slr0453) in Synechocystis sp. PCC 6803. Our in silico analysis showed: (1) 17% flux reduction via RuBisCO for Δpket1 under AC-auto, (2) 11.2-14.3% growth decrease for Δpket2 in turbulent AC-auto, and (3) flux via PKET pathway reaching up to 252% of the flux via phosphoglycerate mutase under AC-auto. All results imply that PKET pathway plays a crucial role under AC-auto by mitigating the decarboxylation occurring in OPP pathway and conversion of pyruvate to acetyl CoA linked to EMP glycolysis under the carbon scarce environment. Finally, our model predicted that PKETs have low affinity to S7P as a substrate.

Advanced integration of fluid dynamics and photosynthetic reaction kinetics for microalgae culture systems.

BACKGROUND: Photosynthetic microalgae have been in the spotlight of biotechnological production (biofuels, lipids, etc), however, current barriers in mass cultivation of microalgae are limiting its successful industrialization. Therefore, a mathematical model integrating both the biological and hydrodynamical parts of the cultivation process may improve our understanding of relevant phenomena, leading to further optimization of the microalgae cultivation. RESULTS: We introduce a unified multidisciplinary simulation tool for microalgae culture systems, particularly the photobioreactors. Our approach describes changes of cell growth determined by dynamics of heterogeneous environmental conditions such as irradiation and mixing of the culture. Presented framework consists of (i) a simplified model of microalgae growth in a culture system (the advection-diffusion-reaction system within a phenomenological model of photosynthesis and photoinhibition), (ii) the fluid dynamics (Navier-Stokes equations), and (iii) the irradiance field description (Beer-Lambert law). To validate the method, a simple case study leading to hydrodynamically induced fluctuating light conditions was chosen. The integration of computational fluid dynamics (ANSYS Fluent) revealed the inner property of the system, the flashing light enhancement phenomenon, known from experiments. CONCLUSION: Our physically accurate model of microalgae culture naturally exhibits features of real system, can be applied to any geometry of microalgae mass cultivation and thus is suitable for biotechnological applications.

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 CO2 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 CO2 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 CO2 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.

Identification of the light-independent phosphoserine pathway as an additional source of serine in the cyanobacterium Synechocystis sp. PCC 6803.

L-serine is one of the proteinogenic amino acids and participates in several essential processes in all organisms. In plants, the light-dependent photorespiratory and the light-independent phosphoserine pathways contribute to serine biosynthesis. In cyanobacteria, the light-dependent photorespiratory pathway for serine synthesis is well characterized, but the phosphoserine pathway has not been identified. Here, we investigated three candidate genes for enzymes of the phosphoserine pathway in Synechocystis sp. PCC 6803. Only the gene for the D-3-phosphoglycerate dehydrogenase is correctly annotated in the genome database, whereas the 3-phosphoserine transaminase and 3-phosphoserine phosphatase (PSP) proteins are incorrectly annotated and were identified here. All enzymes were obtained as recombinant proteins and showed the activities necessary to catalyse the three-step phosphoserine pathway. The genes coding for the phosphoserine pathway were found in most cyanobacterial genomes listed in CyanoBase. The pathway seems to be essential for cyanobacteria, because it was impossible to mutate the gene coding for PSP in Synechocystis sp. PCC 6803 or in Synechococcus elongatus PCC 7942. A model approach indicates a 30-60% contribution of the phosphoserine pathway to the overall serine pool. Hence, this study verified that cyanobacteria, similar to plants, use the phosphoserine pathway in addition to photorespiration for serine biosynthesis.

Multi-level kinetic model explaining diverse roles of isozymes in prokaryotes.

Current standard methods for kinetic and genomic modeling cannot provide deep insight into metabolic regulation. Here, we developed and evaluated a multi-scale kinetic modeling approach applicable to any prokaryote. Specifically, we highlight the primary metabolism of the cyanobacterium Synechococcus elongatus PCC 7942. The model bridges metabolic data sets from cells grown at different CO2 conditions by integrating transcriptomic data and isozymes. Identification of the regulatory roles of isozymes allowed the calculation and explanation of the absolute metabolic concentration of 3-phosphoglycerate. To demonstrate that this method can characterize any isozyme, we determined the function of two glycolytic glyceraldehyde-3-phosphate dehydrogenases: one co-regulates high concentrations of the 3-phosphoglycerate, the other shifts the bifurcation point in hexose regulation, and both improve biomass production. Moreover, the regulatory roles of multiple phosphoglycolate phosphatases were defined for varying (non-steady) CO2 conditions, suggesting their protective role against toxic photorespiratory intermediates.

Phosphoglycerate mutases function as reverse regulated isoenzymes in Synechococcus elongatus PCC 7942.

Phosphoglycerate-mutase (PGM) is an ubiquitous glycolytic enzyme, which in eukaryotic cells can be found in different compartments. In prokaryotic cells, several PGMs are annotated/localized in one compartment. The identification and functional characterization of PGMs in prokaryotes is therefore important for better understanding of metabolic regulation. Here we introduce a method, based on a multi-level kinetic model of the primary carbon metabolism in cyanobacterium Synechococcus elongatus PCC 7942, that allows the identification of a specific function for a particular PGM. The strategy employs multiple parameter estimation runs in high CO2, combined with simulations testing a broad range of kinetic parameters against the changes in transcript levels of annotated PGMs. Simulations are evaluated for a match in metabolic level in low CO2, to reveal trends that can be linked to the function of a particular PGM. A one-isoenzyme scenario shows that PGM2 is a major regulator of glycolysis, while PGM1 and PGM4 make the system robust against environmental changes. Strikingly, combining two PGMs with reverse transcriptional regulation allows both features. A conclusion arising from our analysis is that a two-enzyme PGM system is required to regulate the flux between glycolysis and the Calvin-Benson cycle, while an additional PGM increases the robustness of the system.

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.

On the approaches applied in formulation of a kinetic model of photosystem II: Different approaches lead to different simulations of the chlorophyll alpha fluorescence transients.

Chlorophyll a fluorescence rise (O-J-I-P transient) was in literature simulated using models describing reactions occurring solely in photosystem II (PSII) and plastoquinone (PQ) pool as well as using complex models which described, in addition to the above, also subsequent electron transport occurring beyond the PQ pool. However, there is no consistency in general approach how to formulate a kinetic model and how to describe particular reactions occurring even in PSII only. In this work, simple kinetic PSII models are considered always with the same electron carriers and same type of reactions but some reactions are approached in different ways: oxygen evolving complex is considered bound to PSII or "virtually" separated from PSII; exchange of doubly reduced secondary quinone PSII electron acceptor, Q(B), with PQ molecule from the PQ pool is described by one second order reaction or by two subsequent reactions; and all possible reactions or only those which follow in logical order are considered. By combining all these approaches, eight PSII models are formulated which are used for simulations of the chlorophyll a fluorescence transients. It is shown that the different approaches can lead to qualitatively different results. The approaches are compared with other models found elsewhere in the literature and therefore this work can help the readers to better understand the other models and their results.

Impact of dimeric organization of enzyme on its function: the case of photosynthetic water splitting.

MOTIVATION: It is a question of whether the supramolecular organization of the protein complex has an impact on its function, or not. In the case of the photosystem II (PSII), water splitting might be influenced by cooperation of the PSIIs. Since PSII is the source of the atmospheric oxygen and because better understanding of the water splitting may contribute to the effective use of water as an alternative energy source, possible cooperation should be analyzed and discussed. RESULTS: We suggest that the dimeric organization of the PSII induces cooperation in the water splitting. We show that the model of monomeric PSII is unable to produce the oxygen after the second short flash (associated with the double turnover of the PSII), in contrast to the experimental data and model of dimeric PSII with considered cooperation. On the basis of this fact and partially from the support from other studies, we concluded that the double turnover of the PSII induced by short flashes might be caused by the cooperation in the water splitting. We further discuss a possibility that the known pathway of the electron transport through the PSII might be incomplete and besides D1-Y161, other cofactor which is able to oxidize the special chlorophyll pair (P680) must be considered in the monomeric PSII to explain the oxygen production after the second short flash. AVAILABILITY: Commented SBML codes (.XML files) of the monomeric and dimeric PSII models will be available (at the time of publication) in the BioModels database (www.ebi.ac.uk/biomodels).

Evidence for intermediate S-states as initial phase in the process of oxygen-evolving complex oxidation.

We have analyzed flash-induced period-four damped oscillation of oxygen evolution and chlorophyll fluorescence with the aid of a kinetic model of photosystem II. We have shown that, for simulation of the period-four oscillatory behavior of oxygen evolution, it is essential to consider the so-called intermediate S-state as an initial phase of each of the S(n)-S(n+1), (n = 0, 1, 2, 3) transitions. The intermediate S-states are defined as [S(n)Y(Z)(ox)]-states (n = 0, 1, 2, 3) and are formed with rate constant k(iSn) approximately 1.5 x 10(6) s(-1), which was determined from comparison of theoretical predictions with experimental data. The assumed intermediate S-states shift the equilibrium in reaction P680(+)Y(Z)<-->P680Y(Z)(ox) more to the right and we suggest that kinetics of the intermediate S-states reflects a relaxation process associated with changes of the redox equilibrium in the above reaction. The oxygen oscillation is simulated without the miss and double-hit parameters, if the intermediate S-states, which are not the source of the misses or the double-hits, are included in the simulation. Furthermore, we have shown that the intermediate S-states, together with S(2)Q(A)(-) charge recombination, are prerequisites for the simulation of the period-four oscillatory behavior of the chlorophyll fluorescence.