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

The primary carbon metabolism in cyanobacteria and its regulation

Cyanobacteria are the only prokaryotes capable of performing oxygenic photosynthesis. Many cyanobacterial strains can live in different trophic modes, ranging from photoautotrophic and heterotrophic to mixotrophic growth. However, the regulatory mechanisms allowing a flexible switch between these lifestyles are poorly understood. As anabolic fixation of CO2 in the Calvin-Benson-Bassham (CBB) cycle and catabolic sugar-degradation pathways share intermediates and enzymatic capacity, a tight regulatory network is required to enable simultaneous opposed metabolic fluxes. The Entner-Doudoroff (ED) pathway was recently predicted as one glycolytic route, which cooperates with other pathways in glycogen breakdown. Despite low carbon flux through the ED pathway, metabolite analyses of mutants deficient in the ED pathway revealed a distinct phenotype pointing at a strong regulatory impact of this route. The small Cp12 protein downregulates the CBB cycle in darkness by inhibiting phosphoribulokinase and glyceraldehyde 3-phosphate dehydrogenase. New results of metabolomic and redox level analyses on strains with Cp12 variants extend the known role of Cp12 regulation towards the acclimation to external glucose supply under diurnal conditions as well as to fluctuations in CO2 levels in the light. Moreover, carbon and nitrogen metabolism are closely linked to maintain an essential C/N homeostasis. The small protein PirC was shown to be an important regulator of phosphoglycerate mutase, which identified this enzyme as central branching point for carbon allocation from CBB cycle towards lower glycolysis. Altered metabolite levels in the mutant ΔpirC during nitrogen starvation experiments confirm this regulatory mechanism. The elucidation of novel mechanisms regulating carbon allocation at crucial metabolic branching points could identify ways for targeted redirection of carbon flow towards desired compounds, and thus help to further establish cyanobacteria as green cell factories for biotechnological applications with concurrent utilization of sunlight and CO2.

Double blocking of carbon metabolism causes a large increase of Calvin-Benson cycle compounds in cyanobacteria

Carbon-flow-regulator A (CfrA) adapts carbon flux to nitrogen conditions in nondiazotrophic cyanobacteria. Under nitrogen deficiency, CfrA leads to the storage of excess carbon, which cannot combine with nitrogen, mainly as glycogen. cfrA overexpression from the arsenite-inducible, nitrogen-independent ParsB promoter allows analysis of the metabolic effects of CfrA accumulation. Considering that the main consequence of cfrA overexpression is glycogen accumulation, we examined carbon distribution in response to cfrA expression in Synechocystis sp. PCC 6803 strains impaired in synthesizing this polymer. We carried out a comparative phenotypic analysis to evaluate cfrA overexpression in the wild-type strain and in a mutant of ADP-glucose pyrophosphorylase (ΔglgC), which is unable to synthesize glycogen. The accumulation of CfrA in the wild-type background caused a photosynthetic readjustment although growth was not affected. However, in a ΔglgC strain, growth decreased depending on CfrA accumulation and photosynthesis was severely affected. An elemental analysis of the H, C, and N content of cells revealed that cfrA expression in the wild-type caused an increase in the C/N ratio, due to decreased nitrogen assimilation. Metabolomic study indicated that these cells store sucrose and glycosylglycerol, in addition to the previously described glycogen accumulation. However, cells deficient in glycogen synthesis accumulated large amounts of Calvin-Benson cycle intermediates as cfrA was expressed. These cells also showed increased levels of some amino acids, mainly alanine, serine, valine, isoleucine, and leucine. The findings suggest that by controlling cfrA expression, in different conditions and strains, we could change the distribution of fixed carbon, with potential biotechnological benefits.

Machine learning predicts system-wide metabolic flux control in cyanobacteria

Metabolic fluxes and their control mechanisms are fundamental in cellular metabolism, offering insights for the study of biological systems and biotechnological applications. However, quantitative and predictive understanding of controlling biochemical reactions in microbial cell factories, especially at the system level, is limited. In this work, we present ARCTICA, a computational framework that integrates constraint-based modelling with machine learning tools to address this challenge. Using the model cyanobacterium Synechocystis sp. PCC 6803 as chassis, we demonstrate that ARCTICA effectively simulates global-scale metabolic flux control. Key findings are that (i) the photosynthetic bioproduction is mainly governed by enzymes within the Calvin-Benson-Bassham (CBB) cycle, rather than by those involve in the biosynthesis of the end-product, (ii) the catalytic capacity of the CBB cycle limits the photosynthetic activity and downstream pathways and (iii) ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO) is a major, but not the most, limiting step within the CBB cycle. Predicted metabolic reactions qualitatively align with prior experimental observations, validating our modelling approach. ARCTICA serves as a valuable pipeline for understanding cellular physiology and predicting rate-limiting steps in genome-scale metabolic networks, and thus provides guidance for bioengineering of cyanobacteria.

Photorespiration - Rubisco's repair crew

The photorespiratory repair pathway (photorespiration in short) was set up from ancient metabolic modules about three billion years ago in cyanobacteria, the later ancestors of chloroplasts. These prokaryotes developed the capacity for oxygenic photosynthesis, i.e. the use of water as a source of electrons and protons (with O₂ as a by-product) for the sunlight-driven synthesis of ATP and NADPH for CO₂ fixation in the Calvin cycle. However, the CO₂-binding enzyme, ribulose 1,5-bisphosphate carboxylase (known under the acronym Rubisco), is not absolutely selective for CO₂ and can also use O₂ in a side reaction. It then produces 2-phosphoglycolate (2PG), the accumulation of which would inhibit and potentially stop the Calvin cycle and subsequently photosynthetic electron transport. Photorespiration removes the 2-PG and in this way prevents oxygenic photosynthesis from poisoning itself. In plants, the core of photorespiration consists of ten enzymes distributed over three different types of organelles, requiring interorganellar transport and interaction with several auxiliary enzymes. It goes together with the release and to some extent loss of freshly fixed CO₂. This disadvantageous feature can be suppressed by CO₂-concentrating mechanisms, such as those that evolved in C₄ plants thirty million years ago, which enhance CO₂ fixation and reduce 2PG synthesis. Photorespiration itself provided a pioneer variant of such mechanisms in the predecessors of C₄ plants, C₃-C₄ intermediate plants. This article is a review and update particularly on the enzyme components of plant photorespiration and their catalytic mechanisms, on the interaction of photorespiration with other metabolism and on its impact on the evolution of photosynthesis. This focus was chosen because a better knowledge of the enzymes involved and how they are embedded in overall plant metabolism can facilitate the targeted use of the now highly advanced methods of metabolic network modelling and flux analysis. Understanding photorespiration more than before as a process that enables, rather than reduces, plant photosynthesis, will help develop rational strategies for crop improvement.

The novel PII-interactor PirC identifies phosphoglycerate mutase as key control point of carbon storage metabolism in cyanobacteria

Nitrogen limitation imposes a major transition in the lifestyle of nondiazotrophic cyanobacteria that is controlled by a complex interplay of regulatory factors involving the pervasive signal processor PII Immediately upon nitrogen limitation, newly fixed carbon is redirected toward glycogen synthesis. How the metabolic switch for diverting fixed carbon toward the synthesis of glycogen or of cellular building blocks is operated was so far poorly understood. Here, using the nondiazotrophic cyanobacterium Synechocystis sp. PCC 6803 as model system, we identified a novel PII interactor, the product of the sll0944 gene, which we named PirC. We show that PirC binds to and inhibits the activity of 2,3-phosphoglycerate-independent phosphoglycerate mutase (PGAM), the enzyme that deviates newly fixed CO2 toward lower glycolysis. The binding of PirC to either PII or PGAM is tuned by the metabolite 2-oxoglutarate (2-OG), which accumulates upon nitrogen starvation. In these conditions, the high levels of 2-OG dissociate the PirC-PII complex to promote PirC binding to and inhibition of PGAM. Accordingly, a PirC-deficient mutant showed strongly reduced glycogen levels upon nitrogen deprivation, whereas polyhydroxybutyrate granules were overaccumulated compared to wild-type. Metabolome analysis revealed an imbalance in 3-phosphoglycerate to pyruvate levels in the pirC mutant, confirming that PirC controls the carbon flux in cyanobacteria via mutually exclusive interaction with either PII or PGAM.

Wide range of metabolic adaptations to the acquisition of the Calvin cycle revealed by comparison of microbial genomes

Knowledge of the genetic basis for autotrophic metabolism is valuable since it relates to both the emergence of life and to the metabolic engineering challenge of incorporating CO2 as a potential substrate for biorefining. The most common CO2 fixation pathway is the Calvin cycle, which utilizes Rubisco and phosphoribulokinase enzymes. We searched thousands of microbial genomes and found that 6.0% contained the Calvin cycle. We then contrasted the genomes of Calvin cycle-positive, non-cyanobacterial microbes and their closest relatives by enrichment analysis, ancestral character estimation, and random forest machine learning, to explore genetic adaptations associated with acquisition of the Calvin cycle. The Calvin cycle overlaps with the pentose phosphate pathway and glycolysis, and we could confirm positive associations with fructose-1,6-bisphosphatase, aldolase, and transketolase, constituting a conserved operon, as well as ribulose-phosphate 3-epimerase, ribose-5-phosphate isomerase, and phosphoglycerate kinase. Additionally, carbohydrate storage enzymes, carboxysome proteins (that raise CO2 concentration around Rubisco), and Rubisco activases CbbQ and CbbX accompanied the Calvin cycle. Photorespiration did not appear to be adapted specifically for the Calvin cycle in the non-cyanobacterial microbes under study. Our results suggest that chemoautotrophy in Calvin cycle-positive organisms was commonly enabled by hydrogenase, and less commonly ammonia monooxygenase (nitrification). The enrichment of specific DNA-binding domains indicated Calvin-cycle associated genetic regulation. Metabolic regulatory adaptations were illustrated by negative correlation to AraC and the enzyme arabinose-5-phosphate isomerase, which suggests a downregulation of the metabolite arabinose-5-phosphate, which may interfere with the Calvin cycle through enzyme inhibition and substrate competition. Certain domains of unknown function that were found to be important in the analysis may indicate yet unknown regulatory mechanisms in Calvin cycle-utilizing microbes. Our gene ranking provides targets for experiments seeking to improve CO2 fixation, or engineer novel CO2-fixing organisms.

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 CO₂ 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.

Metabolic model guided strain design of cyanobacteria

Cyanobacteria are oxygenic photoautotrophs that serve as potential platforms for the production of biochemicals from cheap and renewable raw materials - sunlight, water, and carbon dioxide. Systems level analysis of the metabolic network of these organisms could enable the successful engineering of these organisms for the enhanced production of target chemicals. Metabolic modeling techniques including both stoichiometric and kinetic modeling with a genome-wide coverage enable a global assessment of metabolic capabilities. Recent studies guided by such modeling techniques have engineered strains for the enhanced production of valuable chemicals such as ethanol, n-butanol, 1,3-propanediol, glycerol, limonene, and isoprene from CO₂.

Kinetic modeling of the Calvin cycle identifies flux control and stable metabolomes in Synechocystis carbon fixation

Biological fixation of atmospheric CO2 via the Calvin-Benson-Bassham cycle has massive ecological impact and offers potential for industrial exploitation, either by improving carbon fixation in plants and autotrophic bacteria, or by installation into new hosts. A kinetic model of the Calvin-Benson-Bassham cycle embedded in the central carbon metabolism of the cyanobacterium Synechocystis sp. PCC 6803 was developed to investigate its stability and underlying control mechanisms. To reduce the uncertainty associated with a single parameter set, random sampling of the steady-state metabolite concentrations and the enzyme kinetic parameters was employed, resulting in millions of parameterized models which were analyzed for flux control and stability against perturbation. Our results show that the Calvin cycle had an overall high intrinsic stability, but a high concentration of ribulose 1,5-bisphosphate was associated with unstable states. Low substrate saturation and high product saturation of enzymes involved in highly interconnected reactions correlated with increased network stability. Flux control, that is the effect that a change in one reaction rate has on the other reactions in the network, was distributed and mostly exerted by energy supply (ATP), but also by cofactor supply (NADPH). Sedoheptulose 1,7-bisphosphatase/fructose 1,6-bisphosphatase, fructose-bisphosphate aldolase, and transketolase had a weak but positive effect on overall network flux, in agreement with published observations. The identified flux control and relationships between metabolite concentrations and system stability can guide metabolic engineering. The kinetic model structure and parameterizing framework can be expanded for analysis of metabolic systems beyond the Calvin cycle.

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.

Unique attributes of cyanobacterial metabolism revealed by improved genome-scale metabolic modeling and essential gene analysis

The model cyanobacterium, Synechococcus elongatus PCC 7942, is a genetically tractable obligate phototroph that is being developed for the bioproduction of high-value chemicals. Genome-scale models (GEMs) have been successfully used to assess and engineer cellular metabolism; however, GEMs of phototrophic metabolism have been limited by the lack of experimental datasets for model validation and the challenges of incorporating photon uptake. Here, we develop a GEM of metabolism in S. elongatus using random barcode transposon site sequencing (RB-TnSeq) essential gene and physiological data specific to photoautotrophic metabolism. The model explicitly describes photon absorption and accounts for shading, resulting in the characteristic linear growth curve of photoautotrophs. GEM predictions of gene essentiality were compared with data obtained from recent dense-transposon mutagenesis experiments. This dataset allowed major improvements to the accuracy of the model. Furthermore, discrepancies between GEM predictions and the in vivo dataset revealed biological characteristics, such as the importance of a truncated, linear TCA pathway, low flux toward amino acid synthesis from photorespiration, and knowledge gaps within nucleotide metabolism. Coupling of strong experimental support and photoautotrophic modeling methods thus resulted in a highly accurate model of S. elongatus metabolism that highlights previously unknown areas of S. elongatus biology.

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.

Perspectives for a better understanding of the metabolic integration of photorespiration within a complex plant primary metabolism network

Photorespiration is an essential high flux metabolic pathway that is found in all oxygen-producing photosynthetic organisms. It is often viewed as a closed metabolic repair pathway that serves to detoxify 2-phosphoglycolic acid and to recycle carbon to fuel the Calvin-Benson cycle. However, this view is too simplistic since the photorespiratory cycle is known to interact with several primary metabolic pathways, including photosynthesis, nitrate assimilation, amino acid metabolism, C1 metabolism and the Krebs (TCA) cycle. Here we will review recent advances in photorespiration research and discuss future priorities to better understand (i) the metabolic integration of the photorespiratory cycle within the complex network of plant primary metabolism and (ii) the importance of photorespiration in response to abiotic and biotic stresses.

Identification and Profiling of microRNAs Expressed in Elongating Cotton Fibers Using Small RNA Deep Sequencing

Plant microRNAs (miRNAs) have been shown to play essential roles in the regulation of gene expression. In this study, small RNA deep sequencing was applied to explore novel miRNAs expressed in elongating cotton fibers. A total of 46 novel and 96 known miRNAs, primarily derived from the corresponding specific loci in genome of Gossypium arboreum, were identified. 64 miRNAs were shown to be differentially expressed during the fiber elongation process; 16 were predicted to be novel miRNAs while the remaining 48 belong to known miRNA families. Furthermore, RLM-5' RACE (RNA ligase-mediated rapid amplification of 5'-cDNA ends) experiments identified the targets of eight important miRNAs, and the expression levels of these target genes were confirmed to be negatively correlated with the expression patterns of their corresponding miRNAs. We propose a potential functional network mediated through these eight miRNAs to illustrate their important functions in fiber elongation. Our study provides novel insights into the dynamic profiles of these miRNAs and a basis for investigating the regulatory mechanisms involved in the elongation of cotton fibers.

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 CO₂ 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) CO₂ conditions, suggesting their protective role against toxic photorespiratory intermediates.