A native phosphoglycolate salvage pathway of the synthetic autotrophic yeast Komagataella phaffii

Synthetic autotrophs can serve as chassis strains for bioproduction from CO2 as a feedstock to take measures against the climate crisis. Integration of the Calvin-Benson-Bassham (CBB) cycle into the methylotrophic yeast Komagataella phaffii (Pichia pastoris) enabled it to use CO2 as the sole carbon source. The key enzyme in this cycle is ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO) catalyzing the carboxylation step. However, this enzyme is error prone to perform an oxygenation reaction leading to the production of toxic 2-phosphoglycolate. Native autotrophs have evolved different recycling pathways for 2-phosphoglycolate. However, for synthetic autotrophs, no information is available for the existence of such pathways. Deletion of CYB2 in the autotrophic K. phaffii strain led to the accumulation of glycolate, an intermediate in phosphoglycolate salvage pathways, suggesting that such a pathway is enabled by native K. phaffii enzymes. 13C tracer analysis with labeled glycolate indicated that the yeast pathway recycling phosphoglycolate is similar to the plant salvage pathway. This orthogonal yeast pathway may serve as a sensor for RuBisCO oxygenation, and as an engineering target to boost autotrophic growth rates in K. phaffii.

The industrial yeast Pichia pastoris is converted from a heterotroph into an autotroph capable of growth on CO2

The methylotrophic yeast Pichia pastoris is widely used in the manufacture of industrial enzymes and pharmaceuticals. Like most biotechnological production hosts, P. pastoris is heterotrophic and grows on organic feedstocks that have competing uses in the production of food and animal feed. In a step toward more sustainable industrial processes, we describe the conversion of P. pastoris into an autotroph that grows on CO2. By addition of eight heterologous genes and deletion of three native genes, we engineer the peroxisomal methanol-assimilation pathway of P. pastoris into a CO2-fixation pathway resembling the Calvin-Benson-Bassham cycle, the predominant natural CO2-fixation pathway. The resulting strain can grow continuously with CO2 as a sole carbon source at a µmax of 0.008 h-1. The specific growth rate was further improved to 0.018 h-1 by adaptive laboratory evolution. This engineered P. pastoris strain may promote sustainability by sequestering the greenhouse gas CO2, and by avoiding consumption of an organic feedstock with alternative uses in food production.

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.

The independent prokaryotic origins of eukaryotic fructose-1, 6-bisphosphatase and sedoheptulose-1, 7-bisphosphatase and the implications of their origins for the evolution of eukaryotic Calvin cycle

BACKGROUND

In the Calvin cycle of eubacteria, the dephosphorylations of both fructose-1, 6-bisphosphate (FBP) and sedoheptulose-1, 7-bisphosphate (SBP) are catalyzed by the same bifunctional enzyme: fructose-1, 6-bisphosphatase/sedoheptulose-1, 7-bisphosphatase (F/SBPase), while in that of eukaryotic chloroplasts by two distinct enzymes: chloroplastic fructose-1, 6-bisphosphatase (FBPase) and sedoheptulose-1, 7-bisphosphatase (SBPase), respectively. It was proposed that these two eukaryotic enzymes arose from the divergence of a common ancestral eubacterial bifunctional F/SBPase of mitochondrial origin. However, no specific affinity between SBPase and eubacterial FBPase or F/SBPase can be observed in the previous phylogenetic analyses, and it is hard to explain why SBPase and/or F/SBPase are/is absent from most extant nonphotosynthetic eukaryotes according to this scenario.

RESULTS

Domain analysis indicated that eubacterial F/SBPase of two different resources contain distinct domains: proteobacterial F/SBPases contain typical FBPase domain, while cyanobacterial F/SBPases possess FBPase_glpX domain. Therefore, like prokaryotic FBPase, eubacterial F/SBPase can also be divided into two evolutionarily distant classes (Class I and II). Phylogenetic analysis based on a much larger taxonomic sampling than previous work revealed that all eukaryotic SBPase cluster together and form a close sister group to the clade of epsilon-proteobacterial Class I FBPase which are gluconeogenesis-specific enzymes, while all eukaryotic chloroplast FBPase group together with eukaryotic cytosolic FBPase and form another distinct clade which then groups with the Class I FBPase of diverse eubacteria. Motif analysis of these enzymes also supports these phylogenetic correlations.

CONCLUSIONS

There are two evolutionarily distant classes of eubacterial bifunctional F/SBPase. Eukaryotic FBPase and SBPase do not diverge from either of them but have two independent origins: SBPase share a common ancestor with the gluconeogenesis-specific Class I FBPase of epsilon-proteobacteria (or probably originated from that of the ancestor of epsilon-proteobacteria), while FBPase arise from Class I FBPase of an unknown kind of eubacteria. During the evolution of SBPase from eubacterial Class I FBPase, the SBP-dephosphorylation activity was acquired through the transition "from specialist to generalist". The evolutionary substitution of the endosymbiotic-origin cyanobacterial bifunctional F/SBPase by the two light-regulated substrate-specific enzymes made the regulation of the Calvin cycle more delicate, which contributed to the evolution of eukaryotic photosynthesis and even the entire photosynthetic eukaryotes.

Comparative analysis of carboxysome shell proteins

Carboxysomes are metabolic modules for CO(2) fixation that are found in all cyanobacteria and some chemoautotrophic bacteria. They comprise a semi-permeable proteinaceous shell that encapsulates ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO) and carbonic anhydrase. Structural studies are revealing the integral role of the shell protein paralogs to carboxysome form and function. The shell proteins are composed of two domain classes: those with the bacterial microcompartment (BMC; Pfam00936) domain, which oligomerize to form (pseudo)hexamers, and those with the CcmL/EutN (Pfam03319) domain which form pentamers in carboxysomes. These two shell protein types are proposed to be the basis for the carboxysome's icosahedral geometry. The shell proteins are also thought to allow the flux of metabolites across the shell through the presence of the small pore formed by their hexameric/pentameric symmetry axes. In this review, we describe bioinformatic and structural analyses that highlight the important primary, tertiary, and quaternary structural features of these conserved shell subunits. In the future, further understanding of these molecular building blocks may provide the basis for enhancing CO(2) fixation in other organisms or creating novel biological nanostructures.

Inter-species variation in the oligomeric states of the higher plant Calvin cycle enzymes glyceraldehyde-3-phosphate dehydrogenase and phosphoribulokinase

In darkened leaves the Calvin cycle enzymes glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and phosphoribulokinase (PRK) form a regulatory multi-enzyme complex with the small chloroplast protein CP12. GAPDH also forms a high molecular weight regulatory mono-enzyme complex. Given that there are different reports as to the number and subunit composition of these complexes and that enzyme regulatory mechanisms are known to vary between species, it was reasoned that protein-protein interactions may also vary between species. Here, this variation is investigated. This study shows that two different tetramers of GAPDH (an A2B2 heterotetramer and an A4 homotetramer) have the capacity to form part of the PRK/GAPDH/CP12 complex. The role of the PRK/GAPDH/CP12 complex is not simply to regulate the 'non-regulatory' A4 GAPDH tetramer. This study also demonstrates that the abundance and nature of PRK/GAPDH/CP12 interactions are not equal in all species and that whilst NAD enhances complex formation in some species, this is not sufficient for complex formation in others. Furthermore, it is shown that the GAPDH mono-enzyme complex is more abundant as a 2(A2B2) complex, rather than the larger 4(A2B2) complex. This smaller complex is sensitive to cellular metabolites indicating that it is an important regulatory isoform of GAPDH. This comparative study has highlighted considerable heterogeneity in PRK and GAPDH protein interactions between closely related species and the possible underlying physiological basis for this is discussed.

Physiological and molecular changes in Oryza meridionalis Ng., a heat-tolerant species of wild rice

Oryza meridionalis Ng. is a wild relative of Oryza sativa L. found throughout northern Australia where temperatures regularly exceed 35 degrees C in the monsoon growing season. Heat tolerance in O. meridionalis was established by comparing leaf elongation and photosynthetic rates at 45 degrees C with plants maintained at 27 degrees C. By comparison with O. sativa ssp. japonica cv. Amaroo, O. meridionalis was heat tolerant. Elongation rates of the third leaf of O. meridionalis declined by 47% over 24 h at 45 degrees C compared with a 91% decrease for O. sativa. Net photosynthesis was significantly higher in O. sativa at 27 degrees C whereas the two species had the same assimilation rates at 45 degrees C. The leaf proteome and expression levels of individual heat-responsive genes provided insight into the heat response of O. meridionalis. After 24 h of heat exposure, many enzymes involved in the Calvin Cycle were more abundant, while mRNA of their genes generally decreased. Ferredoxin-NADP(H) oxidoreductase, a key enzyme in photosynthetic electron transport had both reduced abundance and gene expression, suggesting light reactions were highly susceptible to heat stress. Rubisco activase was strongly up-regulated after 24 h of heat, with the large isoform having the largest relative increase in protein abundance and a significant increase in gene expression. The protective proteins Cpn60, Hsp90, and Hsp70 all increased in both protein abundance and gene expression. A thiamine biosynthesis protein (THI1), previously shown to act protectively against stress, increased in abundance during heat, even as thiamine levels fell in O. meridionalis.

Thioredoxin-mediated reversible dissociation of a stromal multiprotein complex in response to changes in light availability

A Calvin cycle multiprotein complex including phosphoribulokinase (PRK), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), and a small protein, CP12, has previously been identified. In this article, we have studied this complex in leaves and have shown that dissociation and reassociation of the PRK/GAPDH/CP12 complex occurs in a time frame of minutes, allowing for rapid regulation of enzyme activity. Furthermore, we have shown that the extent of formation and dissociation of the PRK/GAPDH/CP12 complex correlates with the quantity of light. These data provide evidence linking the status of this complex with the rapid and subtle regulation of GAPDH and PRK activities in response to fluctuations in light availability. We have also demonstrated that dissociation of this complex depends on electron transport chain activity and that the major factor involved in the dissociation of the pea complex was thioredoxin f. We show here that both PRK and GAPDH are present in the reduced form in leaves in the dark, but are inactive, demonstrating the role of the PRK/GAPDH/CP12 complex in deactivating these enzymes in response to reductions in light intensity. Based on our data, we propose a model for thioredoxin f-mediated activation of PRK and GAPDH by two mechanisms: directly through reduction of disulfide bonds within these enzymes and indirectly by mediating the breakdown of the complex in response to changes in light intensity.

The structure of isolated Synechococcus strain WH8102 carboxysomes as revealed by electron cryotomography

Carboxysomes are organelle-like polyhedral bodies found in cyanobacteria and many chemoautotrophic bacteria that are thought to facilitate carbon fixation. Carboxysomes are bounded by a proteinaceous outer shell and filled with ribulose 1,5-bisphosphate carboxylase/oxygenase (RuBisCO), the first enzyme in the CO(2) fixation pathway, but exactly how they enhance carbon fixation is unclear. Here we report the three-dimensional structure of purified carboxysomes from Synechococcus species strain WH8102 as revealed by electron cryotomography. We found that while the sizes of individual carboxysomes in this organism varied from 114 nm to 137 nm, surprisingly, all were approximately icosahedral. There were on average approximately 250 RuBisCOs per carboxysome, organized into three to four concentric layers. Some models of carboxysome function depend on specific contacts between individual RuBisCOs and the shell, but no evidence of such contacts was found: no systematic patterns of connecting densities or RuBisCO positions against the shell's presumed hexagonal lattice could be discerned, and simulations showed that packing forces alone could account for the layered organization of RuBisCOs.

The metabolic significance of octulose phosphates in the photosynthetic carbon reduction cycle in spinach

14C-Labelled octulose phosphates were formed during photosynthetic 14CO2 fixation and were measured in spinach leaves and chloroplasts. Because mono- and bisphosphates of D: -glycero- D: -ido-octulose are the active 8-carbon ketosugar intermediates of the L-type pentose pathway, it was proposed that they may also be reactants in a modified Calvin-Benson-Bassham pathway reaction scheme. This investigation therefore initially focussed only on the ido-epimer of the octulose phosphates even though 14C-labelled D: -glycero- D: -altro-octulose mono- and bisphosphates were also identified in chloroplasts and leaves. 14CO2 predominantly labelled positions 5 and 6 of D: -glycero- D: -ido-octulose 1,8-P2 consistent with labelling predictions of the modified scheme. The kinetics of 14CO2 incorporation into ido-octulose was similar to its incorporation into some traditional intermediates of the path of carbon, while subsequent exposure to 12CO2 rapidly displaced the 14C isotope label from octulose with the same kinetics of label loss as some of the confirmed Calvin pathway intermediates. This is consistent with octulose phosphates having the role of cyclic intermediates rather than synthesized storage products. (Storage products don't rapidly exchange isotopically labelled carbons with unlabelled CO2.)A spinach chloroplast extract, designated stromal enzyme preparation (SEP), catalysed and was used to measure rates of CO(2) assimilation with Calvin cycle intermediates and octulose and arabinose phosphates. Only pentose (but not arabinose) phosphates and sedoheptulose 7-phosphate supported CO2 fixation at rates in excess of 120 micromol h(-1) mg(-1) Chl. Rates for octulose, sedoheptulose and fructose bisphosphates, octulose, hexose and triose monophosphates were all notably less than the above rate and arabinose 5-phosphate was inactive. Altro-octulose phosphates were more active than phosphate esters of the ido-epimer. The modified scheme proposed a specific phosphotransferase and SEP unequivocally catalysed reversible phosphate transfer between sedoheptulose bisphosphate and D: -glycero- D: -ido-octulose 8-phosphate. It was also initially hypothesized that arabinose 5-phosphate, an L-Type pentose pathway reactant, may have a role in a modified Calvin pathway. Arabinose 5-phosphate is present in spinach chloroplasts and leaves. Radiochromatography showed that 14C-arabinose 5-phosphate with SEP, but only in the presence of an excess of unlabelled ribose 5-phosphate, lightly labelled ribulose 5-phosphate and more heavily labelled hexose and sedoheptulose mono- and bisphosphates. However, failure to demonstrate any CO2 fixation by arabinose 5-phosphate as sole substrate suggested that the above labelling may have no metabolic significance. Despite this arabinose and ribose 5-phosphates are shown to exhibit active roles as enzyme co-factors in transaldolase and aldolase exchange reactions that catalyse the epimeric interconversions of the phosphate esters of ido- and altro-octulose. Arabinose 5-phosphate is presented as playing this role in a New Reaction Scheme for the path of carbon, where it is concluded that slow reacting ido-octulose 1,8 bisphosphate has no role. The more reactive altro-octulose phosphates, which are independent of the need for phosphotransferase processing, are presented as intermediates in the new scheme. Moreover, using the estimates of phosphotransferase activity with altro-octulose monophosphate as substrate allowed calculation of the contributions of the new scheme, that ranged from 11% based on the intact chloroplast carboxylation rate to 80% using the carboxylation rate required for the support of octulose phosphate synthesis and its role in the phosphotransferase reaction.

Glutathionylation of chloroplast thioredoxin f is a redox signaling mechanism in plants

Thioredoxin f (TRXf) is a key factor in the redox regulation of chloroplastic carbon fixation enzymes, whereas glutathione is an important thiol buffer whose status is modulated by stress conditions. Here, we report specific glutathionylation of TRXf. A conserved cysteine is present in the TRXf primary sequence, in addition to its two active-site cysteines. The additional cysteine becomes glutathionylated when TRXf is exposed to oxidized glutathione or to reduced glutathione plus oxidants. No other chloroplastic TRX, from either Arabidopsis or Chlamydomonas, is glutathionylated under these conditions. Glutathionylation decreases the ability of TRXf to be reduced by ferredoxin-thioredoxin reductase and results in impaired light activation of target enzymes in a reconstituted thylakoid system. Although several mammalian proteins undergoing glutathionylation have already been identified, TRXf is among the first plant proteins found to undergo this posttranslational modification. This report suggests that a crosstalk between the TRX and glutathione systems mediates a previously uncharacterized form of redox signaling in plants in stress conditions.

Structure of Escherichia coli ribose-5-phosphate isomerase: a ubiquitous enzyme of the pentose phosphate pathway and the Calvin cycle

Ribose-5-phosphate isomerase A (RpiA; EC 5.3.1.6) interconverts ribose-5-phosphate and ribulose-5-phosphate. This enzyme plays essential roles in carbohydrate anabolism and catabolism; it is ubiquitous and highly conserved. The structure of RpiA from Escherichia coli was solved by multiwavelength anomalous diffraction (MAD) phasing, and refined to 1.5 A resolution (R factor 22.4%, R(free) 23.7%). RpiA exhibits an alpha/beta/(alpha/beta)/beta/alpha fold, some portions of which are similar to proteins of the alcohol dehydrogenase family. The two subunits of the dimer in the asymmetric unit have different conformations, representing the opening/closing of a cleft. Active site residues were identified in the cleft using sequence conservation, as well as the structure of a complex with the inhibitor arabinose-5-phosphate at 1.25 A resolution. A mechanism for acid-base catalysis is proposed.