Mutant studies of phosphofructo-2-kinases do not reveal an essential role of fructose-2,6-bisphosphate in the regulation of carbon fluxes in yeast cells

Induction of point and structural mutations in engineered yeast Saccharomyces cerevisiae improve carotenoid production

β-Carotene is an attractive compound and that its biotechnological production can be achieved by using engineered Saccharomyces cerevisiae. In a previous study, we developed a technique for the efficient establishment of diverse mutants through the introduction of point and structural mutations into the yeast genome. In this study, we aimed to improve β-carotene production by applying this mutagenesis technique to S. cerevisiae strain that had been genetically engineered for β-carotene production. Point and structural mutations were introduced into β-carotene-producing engineered yeast. The resulting mutants showed higher β-carotene production capacity than the parental strain. The top-performing mutant, HP100_74, produced 37.6 mg/L of β-carotene, a value 1.9 times higher than that of the parental strain (20.1 mg/L). Gene expression analysis confirmed an increased expression of multiple genes in the glycolysis, mevalonate, and β-carotene synthesis pathways. In contrast, expression of ERG9, which functions in the ergosterol pathway competing with β-carotene production, was decreased in the mutant strain. The introduction of point and structural mutations represents a simple yet effective method for achieving mutagenesis in yeasts. This technique is expected to be widely applied in the future to produce chemicals via metabolic engineering of S. cerevisiae.

Gene expression regulates metabolite homeostasis during the Crabtree effect: Implications for the adaptation and evolution of Metabolism

A key issue in both molecular and evolutionary biology has been to define the roles of genes and phenotypes in the adaptation of organisms to environmental changes. The dominant view has been that an organism's metabolic adaptations are driven by gene expression and that gene mutations, independent of the starting phenotype, are responsible for the evolution of new metabolic phenotypes. We propose an alternate hypothesis, in which the phenotype and genotype together determine metabolic adaptation both in the lifetime of the organism and in the evolutionary selection of adaptive metabolic traits. We tested this hypothesis by flux-balance and metabolic-control analysis of the relative roles of the starting phenotype and gene expression in regulating the metabolic adaptations during the Crabtree effect in yeast, when they are switched from a low- to high-glucose environment. Critical for successful short-term adaptation was the ability of the glycogen/trehalose shunt to balance the glycolytic pathway. The role of later gene expression of new isoforms of glycolytic enzymes, rather than flux control, was to provide additional homeostatic mechanisms allowing an increase in the amount and efficiency of adenosine triphosphate and product formation while maintaining glycolytic balance. We further showed that homeostatic mechanisms, by allowing increased phenotypic plasticity, could have played an important role in guiding the evolution of the Crabtree effect. Although our findings are specific to Crabtree yeast, they are likely to be broadly found because of the well-recognized similarities in glucose metabolism across kingdoms and phyla from yeast to humans.

Identifying the Growth Modulon of Corynebacterium glutamicum

The growth rate (μ) of industrially relevant microbes, such as Corynebacterium glutamicum, is a fundamental property that indicates its production capacity. Therefore, understanding the mechanism underlying the growth rate is imperative for improving productivity and performance through metabolic engineering. Despite recent progress in the understanding of global regulatory interactions, knowledge of mechanisms directing cell growth remains fragmented and incomplete. The current study investigated RNA-Seq data of three growth rate transitions, induced by different pre-culture conditions, in order to identify transcriptomic changes corresponding to increasing growth rates. These transitions took place in minimal medium and ranged from 0.02 to 0.4 h^-1 μ. This study enabled the identification of 447 genes as components of the growth modulon. Enrichment of genes within the growth modulon revealed 10 regulons exhibiting a significant effect over growth rate transition. In summary, central metabolism was observed to be regulated by a combination of metabolic and transcriptional activities orchestrating control over glycolysis, pentose phosphate pathway, and the tricarboxylic acid cycle. Additionally, major responses to changes in the growth rate were linked to iron uptake and carbon metabolism. In particular, genes encoding glycolytic enzymes and the glucose uptake system showed a positive correlation with the growth rate.

Mutations in the phosphatase domain of the 6-phosphofructo-2-kinase/fructose 2,6-bisphosphatase result in the transcriptional activation of the alternative oxidase and gluconeogenic pathways in Podospora anserina

By screening suppressors of a respiratory mutant lacking a functional cytochrome pathway in the filamentous fungus Podospora anserina, we isolated a mutation located in the phosphatase domain of the bi-functional enzyme 6-phosphofructo-2-kinase/fructose 2,6-bisphosphatase (PFK-2/FBPase-2). We show that the inactivation of the phosphatase but not of the kinase domain is responsible for the suppressor effect that results from the activation of the RSEs transcription factors that control expression of AOX, an alternative oxidase able to bypass the mitochondria cytochrome pathway of respiration. Remarkably, activation of the RSEs also stimulates the expression of the gluconeogenic enzymes, fructose-1,6 bi-phosphatase (FBPase-1) and phosphoenolpyruvate carboxykinase (PCK-1). We thus reveal in P. anserina an apparently paradoxical situation where the inactivation of the phosphatase domain of PFK-2/FBPase-2, supposed to stimulate glycolysis, is correlated with the transcriptional induction of the gluconeogenic enzymes. Phylogenic analysis revealed the presence of multiple presumed PFK-2/FBPase-2 isoforms in all the species of tested Ascomycetes.

Combination of Three Methods to Reduce Glucose Metabolic Rate For Improving N-Acetylglucosamine Production in Saccharomyces cerevisiae

Previously, the production of N-acetylglucosamine (GlcNAc) in Saccharomyces cerevisiae was improved by deletion of the genes encoding phosphofructokinase 2 (PFK-2) isoforms, which reduced the glycolytic flux by eliminating the pathway to produce fructose-2,6-bisphosphate, an allosteric activator of phosphofructokinase 1 (PFK-1). We further examined the effects of an additional reduction in glucose metabolic rate on N-acetylglucosamine production. Glucose uptake rate was lowered by expressing a gene encoding truncated glucose-sensing regulator ( MTH1-Δ T). In addition, catalytically dead Cas9 (dCas9) was introduced in order to down-regulate the expression levels of PFK-1 and pyruvate kinase-1 (Pyk1). Finally, the three strategies were introduced into S. cerevisiae strains in a combinatorial way; the strain containing all three modules resulted in the highest N-acetylglucosamine production yield. The results showed that the three modules cooperatively reduced the glucose metabolism and improved N-acetylglucosamine production up to 3.0 g/L in shake flask cultivation.

Flux control through protein phosphorylation in yeast

Protein phosphorylation is one of the most important mechanisms regulating metabolism as it can directly modify metabolic enzymes by the addition of phosphate groups. Attributed to such a rapid and reversible mechanism, cells can adjust metabolism rapidly in response to temporal changes. The yeast Saccharomyces cerevisiae, a widely used cell factory and model organism, is reported to show frequent phosphorylation events in metabolism. Studying protein phosphorylation in S. cerevisiae allows for gaining new insight into the function of regulatory networks, which may enable improved metabolic engineering as well as identify mechanisms underlying human metabolic diseases. Here we collect functional phosphorylation events of 41 enzymes involved in yeast metabolism and demonstrate functional mechanisms and the application of this information in metabolic engineering. From a systems biology perspective, we describe the development of phosphoproteomics in yeast as well as approaches to analysing the phosphoproteomics data. Finally, we focus on integrated analyses with other omics data sets and genome-scale metabolic models. Despite the advances, future studies improving both experimental technologies and computational approaches are imperative to expand the current knowledge of protein phosphorylation in S. cerevisiae.

Effect of the cancer specific shorter form of human 6-phosphofructo-1-kinase on the metabolism of the yeast Saccharomyces cerevisiae

Background

At first glance, there appears to be a high degree of similarity between the metabolism of yeast (the Crabtree effect) and human cancer cells (the Warburg effect). At the root of both effects is accelerated metabolic flow through glycolysis which leads to overflows of ethanol and lactic acid, respectively. It has been proposed that enhanced glycolytic flow in cancer cells is triggered by the altered kinetic characteristics of the key glycolytic regulatory enzyme 6-phosphofructo-1-kinase (Pfk). Through a posttranslational modification, highly active shorter Pfk-M fragments, which are resistant to feedback inhibition, are formed after the proteolytic cleavage of the C-terminus of the native human Pfk-M. Alternatively, enhanced glycolysis is triggered by optimal growth conditions in the yeast Saccharomyces cerevisiae.

Results

To assess the deregulation of glycolysis in yeast cells, the sfPFKM gene encoding highly active human shorter Pfk-M fragments was introduced into pfk-null S. cerevisiae. No growth of the transformants with the sfPFKM gene was observed on glucose and fructose. Glucose even induced rapid deactivation of Pfk1 activities in such transformants. However, Pfk1 activities of the sfPFKM transformants were detected in maltose medium, but the growth in maltose was possible only after the addition of 10 mM of ethanol to the medium. Ethanol seemed to alleviate the severely unbalanced NADH/NADPH ratio in the sfPFKM cells. However, the transformants carrying modified Pfk-M enzymes grew faster than the transformants with the human native human Pfk-M enzyme in a narrow ecological niche with a low maltose concentration medium that was further improved by additional modifications. Interestingly, periodic extracellular accumulation of phenylacetaldehyde was detected during the growth of the strain with modified Pfk-M but not with the strain encoding the human native enzyme.

Conclusions

Highly active cancer-specific shorter Pfk-M fragments appear to trigger several controlling mechanisms in the primary metabolism of yeast S. cerevisiae cells. These results suggest more complex metabolic regulation is present in S. cerevisiae as free living unicellular eukaryotic organisms in comparison to metazoan human cells. However, increased productivity under broader growth conditions may be achieved if more gene engineering is performed to reduce or omit several controlling mechanisms.

Electronic supplementary material

The online version of this article (doi:10.1186/s12896-017-0362-5) contains supplementary material, which is available to authorized users.

Global Metabolic Engineering of Glycolytic Pathway via Multicopy Integration in Saccharomyces cerevisiae

The use of renewable feedstocks for producing biofuels and biobased chemicals by engineering metabolic pathways of yeast Saccharomyces cerevisiae has recently become an attractive option. Many researchers attempted to increase glucose consumption rate by overexpressing some glycolytic enzymes because most target biobased chemicals are derived through glycolysis. However, these attempts have met with little success. In this study, to create a S. cerevisiae strain with high glucose consumption rate, we used multicopy integration to develop a global metabolic engineering strategy. Among approximately 350 metabolically engineered strains, YPH499/dPdA3-34 exhibited the highest glucose consumption rate. This strain showed 1.3-fold higher cell growth rate and glucose consumption rate than the control strain. Real-time PCR analysis revealed that transcription levels of glycolysis-related genes such as HXK2, PFK1, PFK2, PYK2, PGI1, and PGK1 in YPH499/dPdA3-34 were increased. Our strategy is thus a promising approach to optimize global metabolic pathways in S. cerevisiae.

Control and regulation of pathways via negative feedback

The biochemical networks found in living organisms include a huge variety of control mechanisms at multiple levels of organization. While the mechanistic and molecular details of many of these control mechanisms are understood, their exact role in driving cellular behaviour is not. For example, yeast glycolysis has been studied for almost 80 years but it is only recently that we have come to understand the systemic role of the multitude of feedback and feed-forward controls that exist in this pathway. In this article, control theory is discussed as an approach to dissect the control logic of complex pathways. One of the key issues is distinguishing between the terms control and regulation and how these concepts are applied to regulated enzymes such as phosphofructokinase. In doing so, one of the paradoxes in metabolic regulation can be resolved where enzymes such as phosphofructokinase have little control but, nevertheless, possess significant regulatory influence.

Multi-tasking of biosynthetic and energetic functions of glycolysis explained by supply and demand logic

After more than a century of research on glycolysis, we have detailed descriptions of its molecular organization, but despite this wealth of knowledge, linking the enzyme properties to metabolic pathway behavior remains challenging. These challenges arise from multi-layered regulation and the context and time dependence of component functions. However, when viewed as a system that functions according to the principles of supply and demand, a simplifying theoretical framework can be applied to study its regulation logic and to assess the coherence of experimental interpretations. These principles are universally applicable, as they emphasize the common metabolic tasks of glycolysis: the provision of free-energy carriers, and precursors for biosynthesis and stress-related compounds. Here we will review the regulation of multi-tasking by glycolysis and consider how an understanding of this central metabolic pathway can be pursued using general principles, rather than focusing on the biochemical details of constituent components.

Yeast phosphofructokinase-1 subunit Pfk2p is necessary for pH homeostasis and glucose-dependent vacuolar ATPase reassembly

V-ATPases are conserved ATP-driven proton pumps that acidify organelles. Yeast V-ATPase assembly and activity are glucose-dependent. Glucose depletion causes V-ATPase disassembly and its inactivation. Glucose readdition triggers reassembly and resumes proton transport and organelle acidification. We investigated the roles of the yeast phosphofructokinase-1 subunits Pfk1p and Pfk2p for V-ATPase function. The pfk1Δ and pfk2Δ mutants grew on glucose and assembled wild-type levels of V-ATPase pumps at the membrane. Both phosphofructokinase-1 subunits co-immunoprecipitated with V-ATPase in wild-type cells; upon deletion of one subunit, the other subunit retained binding to V-ATPase. The pfk2Δ cells exhibited a partial vma growth phenotype. In vitro ATP hydrolysis and proton transport were reduced by 35% in pfk2Δ membrane fractions; they were normal in pfk1Δ. In vivo, the pfk1Δ and pfk2Δ vacuoles were alkalinized and the cytosol acidified, suggestive of impaired V-ATPase proton transport. Overall the pH alterations were more dramatic in pfk2Δ than pfk1Δ at steady state and after readdition of glucose to glucose-deprived cells. Glucose-dependent reassembly was 50% reduced in pfk2Δ, and the vacuolar lumen was not acidified after reassembly. RAVE-assisted glucose-dependent reassembly and/or glucose signals were disturbed in pfk2Δ. Binding of disassembled V-ATPase (V1 domain) to its assembly factor RAVE (subunit Rav1p) was 5-fold enhanced, indicating that Pfk2p is necessary for V-ATPase regulation by glucose. Because Pfk1p and Pfk2p are necessary for V-ATPase proton transport at the vacuole in vivo, a role for glycolysis at regulating V-ATPase proton transport is discussed.

Control of glycolytic flux in directed biosynthesis of uridine-phosphoryl compounds through the manipulation of ATP availability

Adenosine triphosphate (ATP), the most important energy source for metabolic reactions and pathways, plays a vital role in control of metabolic flux. Considering the importance of ATP in regulation of the glycolytic pathway, the use of ATP-oriented manipulation is a rational and efficient route to regulate metabolic flux. In this paper, a series of efficient ATP-oriented regulation methods, such as changing ambient temperature and altering reduced nicotinamide adenine dinucleotide (NADH), was developed. To satisfy the different demand for ATP at different phases in directed biosynthesis of uridine-phosphoryl compounds, a multiphase ATP supply regulation strategy was also used to enhance to yield of target metabolites.

An efficient xylose-fermenting recombinant Saccharomyces cerevisiae strain obtained through adaptive evolution and its global transcription profile

Factors related to ethanol production from xylose in engineered Saccharomyces cerevisiae that contain an exogenous initial metabolic pathway are still to be elucidated. In the present study, a strain that expresses the xylose isomerase gene of Piromyces sp. Pi-xylA and overexpresses XKS1, RPE1, RKI1, TAL1, and TKL1, with deleted GRE3 and COX4 genes was constructed. The xylose utilization capacity of the respiratory deficiency strain was poor but improved via adaptive evolution in xylose. The μ (max) of the evolved strain in 20 g l(-1) xylose is 0.11 ± 0.00 h(-1), and the evolved strain consumed 17.83 g l(-1) xylose within 72 h, with an ethanol yield of 0.43 g g(-1) total consumed sugars during glucose-xylose cofermentation. Global transcriptional changes and effect of several specific genes were studied. The result revealed that the increased xylose isomerase acivity, the upregulation of enzymes involved in glycolysis and glutamate synthesis, and the downregulation of trehalose and glycogen synthesis, may have contributed to the improved xylose utilization of the strain. Furthermore, the deletion of PHO13 decreased the xylose growth in the respiration deficiency strain although deleting PHO13 can improve the xylose metabolism in other strains.

Modeling the optimal central carbon metabolic pathways under feedback inhibition using flux balance analysis

Metabolism is a complex process for energy production for cellular activity. It consists of a cascade of reactions that form a highly branched network in which the product of one reaction is the reactant of the next reaction. Metabolic pathways efficiently produce maximal amount of biomass while maintaining a steady-state behavior. The steady-state activity of such biochemical pathways necessarily incorporates feedback inhibition of the enzymes. This observation motivates us to incorporate feedback inhibition for modeling the optimal activity of metabolic pathways using flux balance analysis (FBA). We demonstrate the effectiveness of the methodology on a synthetic pathway with and without feedback inhibition. Similarly, for the first time, the Central Carbon Metabolic (CCM) pathways of Saccharomyces cerevisiae and Homo sapiens have been modeled and compared based on the above understanding. The optimal pathway, which maximizes the amount of the target product(s), is selected from all those obtained by the proposed method. For this, we have observed the concentration of the product inhibited enzymes of CCM pathway and its influence on its corresponding metabolite/substrate. We have also studied the concentration of the enzymes which are responsible for the synthesis of target products. We further hypothesize that an optimal pathway would opt for higher flux rate reactions. In light of these observations, we can say that an optimal pathway should have lower enzyme concentration and higher flux rates. Finally, we demonstrate the superiority of the proposed method by comparing it with the extreme pathway analysis.

Overexpression of genes encoding glycolytic enzymes in Corynebacterium glutamicum enhances glucose metabolism and alanine production under oxygen deprivation conditions

We previously reported that Corynebacterium glutamicum strain ΔldhAΔppc+alaD+gapA, overexpressing glyceraldehyde-3-phosphate dehydrogenase-encoding gapA, shows significantly improved glucose consumption and alanine formation under oxygen deprivation conditions (T. Jojima, M. Fujii, E. Mori, M. Inui, and H. Yukawa, Appl. Microbiol. Biotechnol. 87:159-165, 2010). In this study, we employ stepwise overexpression and chromosomal integration of a total of four genes encoding glycolytic enzymes (herein referred to as glycolytic genes) to demonstrate further successive improvements in C. glutamicum glucose metabolism under oxygen deprivation. In addition to gapA, overexpressing pyruvate kinase-encoding pyk and phosphofructokinase-encoding pfk enabled strain GLY2/pCRD500 to realize respective 13% and 20% improved rates of glucose consumption and alanine formation compared to GLY1/pCRD500. Subsequent overexpression of glucose-6-phosphate isomerase-encoding gpi in strain GLY3/pCRD500 further improved its glucose metabolism. Notably, both alanine productivity and yield increased after each overexpression step. After 48 h of incubation, GLY3/pCRD500 produced 2,430 mM alanine at a yield of 91.8%. This was 6.4-fold higher productivity than that of the wild-type strain. Intracellular metabolite analysis showed that gapA overexpression led to a decreased concentration of metabolites upstream of glyceraldehyde-3-phosphate dehydrogenase, suggesting that the overexpression resolved a bottleneck in glycolysis. Changing ratios of the extracellular metabolites by overexpression of glycolytic genes resulted in reduction of the intracellular NADH/NAD(+) ratio, which also plays an important role on the improvement of glucose consumption. Enhanced alanine dehydrogenase activity using a high-copy-number plasmid further accelerated the overall alanine productivity. Increase in glycolytic enzyme activities is a promising approach to make drastic progress in growth-arrested bioprocesses.

Control of ATP homeostasis during the respiro-fermentative transition in yeast

Respiring Saccharomyces cerevisiae cells respond to a sudden increase in glucose concentration by a pronounced drop of their adenine nucleotide content ([ATP]+[ADP]+[AMP]=[AXP]). The unknown fate of 'lost' AXP nucleotides represented a long-standing problem for the understanding of the yeast's physiological response to changing growth conditions. Transient accumulation of the purine salvage pathway intermediate, inosine, accounted for the apparent loss of adenine nucleotides. Conversion of AXPs into inosine was facilitated by AMP deaminase, Amd1, and IMP-specific 5'-nucleotidase, Isn1. Inosine recycling into the AXP pool was facilitated by purine nucleoside phosphorylase, Pnp1, and joint action of the phosphoribosyltransferases, Hpt1 and Xpt1. Analysis of changes in 24 intracellular metabolite pools during the respiro-fermentative growth transition in wild-type, amd1, isn1, and pnp1 strains revealed that only the amd1 mutant exhibited significant deviations from the wild-type behavior. Moreover, mutants that were blocked in inosine production exhibited delayed growth acceleration after glucose addition. It is proposed that interconversion of adenine nucleotides and inosine facilitates rapid and energy-cost efficient adaptation of the AXP pool size to changing environmental conditions.

Regulation of gluconeogenesis in Saccharomyces cerevisiae is mediated by activator and repressor functions of Rds2

In Saccharomyces cerevisiae, RDS2 encodes a zinc cluster transcription factor with unknown function. Here, we unravel a key function of Rds2 in gluconeogenesis using chromatin immunoprecipitation-chip technology. While we observed that Rds2 binds to only a few promoters in glucose-containing medium, it binds many additional genes when the medium is shifted to ethanol, a nonfermentable carbon source. Interestingly, many of these genes are involved in gluconeogenesis, the tricarboxylic acid cycle, and the glyoxylate cycle. Importantly, we show that Rds2 has a dual function: it directly activates the expression of gluconeogenic structural genes while it represses the expression of negative regulators of this pathway. We also show that the purified DNA binding domain of Rds2 binds in vitro to carbon source response elements found in the promoters of target genes. Finally, we show that upon a shift to ethanol, Rds2 activation is correlated with its hyperphosphorylation by the Snf1 kinase. In summary, we have characterized Rds2 as a novel major regulator of gluconeogenesis.

Increase in Fru-2,6-P(2) levels results in altered cell division in Schizosaccharomyces pombe

Mitogenic response to growth factors is concomitant with the modulation they exert on the levels of Fructose 2,6-bisphosphate (Fru-2,6-P2), an essential activator of the glycolytic flux. In mammalian cells, decreased Fru-2,6-P2 concentration causes cell cycle delay, whereas high levels of Fru-2,6-P2 sensitize cells to apoptosis. In order to analyze the cell cycle consequences due to changes in Fru-2,6-P2 levels, the bisphosphatase-dead mutant (H258A) of 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase enzyme was over-expressed in Schizosaccharomyces pombe cells and the variation in cell phenotype was studied. The results obtained demonstrate that the increase in Fru-2,6-P2 levels results in a defective division of S. pombe, as revealed by an altered multisepted phenotype. The H258A-expressing cells showed impairment of cytokinesis, but normal nuclear division. In order to identify cellular mediators responsible for this effect, we transformed different S. pombe strains and observed that the cytokinetic defect was absent in cells defective for Wee1 kinase function. Therefore, in S. pombe, Wee1 integrates the metabolic signal emerging from changes in Fru-2,6-P2 content, thus coupling metabolism with cell proliferation. As the key regulators of the cell cycle checkpoints are conserved throughout evolution, these results may help to understand the experimental evidences obtained by manipulation of Fru-2,6-P2 levels in mammalian cells.

Functional genomic analysis of commercial baker's yeast during initial stages of model dough-fermentation

Gene expression profiles of baker's yeast during initial dough-fermentation were investigated using liquid fermentation (LF) media to obtain insights at the molecular level into rapid adaptation mechanisms of baker's yeast. Results showed that onset of fermentation caused drastic changes in gene expression profiles within 15 min. Genes involved in the tricarboxylic acid (TCA) cycle were down-regulated and genes involved in glycolysis were up-regulated, indicating a metabolic shift from respiration to fermentation. Genes involved in ethanol production (PDC genes and ADH1), in glycerol synthesis (GPD1 and HOR2), and in low-affinity hexose transporters (HXT1 and HXT3) were up-regulated at the beginning of model dough-fermentation. Among genes up-regulated at 15 min, several genes classified as transcription were down-regulated within 30 min. These down-regulated genes are involved in messenger RNA splicing and ribosomal protein biogenesis and in transcriptional regulator (SRB8, MIG1). In contrast, genes involved in amino acid metabolism and in vitamin metabolism, such as arginine biosynthesis, riboflavin biosynthesis, and thiamin biosynthesis, were subsequently up-regulated after 30 min. Interestingly, the genes involved in the unfolded protein response (UPR) pathway were also subsequently up-regulated. Our study presents the first overall description of the transcriptional response of baker's yeast during dough-fermentation, and will thus help clarify genomic responses to various stresses during commercial fermentation processes.

Metabolic functions of duplicate genes in Saccharomyces cerevisiae

The roles of duplicate genes and their contribution to the phenomenon of enzyme dispensability are a central issue in molecular and genome evolution. A comprehensive classification of the mechanisms that may have led to their preservation, however, is currently lacking. In a systems biology approach, we classify here back-up, regulatory, and gene dosage functions for the 105 duplicate gene families of Saccharomyces cerevisiae metabolism. The key tool was the reconciled genome-scale metabolic model iLL672, which was based on the older iFF708. Computational predictions of all metabolic gene knockouts were validated with the experimentally determined phenotypes of the entire singleton yeast library of 4658 mutants under five environmental conditions. iLL672 correctly identified 96%-98% and 73%-80% of the viable and lethal singleton phenotypes, respectively. Functional roles for each duplicate family were identified by integrating the iLL672-predicted in silico duplicate knockout phenotypes, genome-scale carbon-flux distributions, singleton mutant phenotypes, and network topology analysis. The results provide no evidence for a particular dominant function that maintains duplicate genes in the genome. In particular, the back-up function is not favored by evolutionary selection because duplicates do not occur more frequently in essential reactions than singleton genes. Instead of a prevailing role, multigene-encoded enzymes cover different functions. Thus, at least for metabolism, persistence of the paralog fraction in the genome can be better explained with an array of different, often overlapping functional roles.

6-Phosphofructo-2-kinase and fructose-2,6-bisphosphatase in Trypanosomatidae. Molecular characterization, database searches, modelling studies and evolutionary analysis

Fructose 2,6-bisphosphate is a potent allosteric activator of trypanosomatid pyruvate kinase and thus represents an important regulator of energy metabolism in these protozoan parasites. A 6-phosphofructo-2-kinase, responsible for the synthesis of this regulator, was highly purified from the bloodstream form of Trypanosoma brucei and kinetically characterized. By searching trypanosomatid genome databases, four genes encoding proteins homologous to the mammalian bifunctional enzyme 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase (PFK-2/FBPase-2) were found for both T. brucei and the related parasite Leishmania major and four pairs in Trypanosoma cruzi. These genes were predicted to each encode a protein in which, at most, only a single domain would be active. Two of the T. brucei proteins showed most conservation in the PFK-2 domain, although one of them was predicted to be inactive due to substitution of residues responsible for ligating the catalytically essential divalent metal cation; the two other proteins were most conserved in the FBPase-2 domain. The two PFK-2-like proteins were expressed in Escherichia coli. Indeed, the first displayed PFK-2 activity with similar kinetic properties to that of the enzyme purified from T. brucei, whereas no activity was found for the second. Interestingly, several of the predicted trypanosomatid PFK-2/FBPase-2 proteins have long N-terminal extensions. The N-terminal domains of the two polypeptides with most similarity to mammalian PFK-2s contain a series of tandem repeat ankyrin motifs. In other proteins such motifs are known to mediate protein-protein interactions. Phylogenetic analysis suggests that the four different PFK-2/FBPase-2 isoenzymes found in Trypanosoma and Leishmania evolved from a single ancestral bifunctional enzyme within the trypanosomatid lineage. A possible explanation for the evolution of multiple monofunctional enzymes and for the presence of the ankyrin-motif repeats in the PFK-2 isoenzymes is presented.

Role of transcriptional regulation in controlling fluxes in central carbon metabolism of Saccharomyces cerevisiae. A chemostat culture study

In contrast to batch cultivation, chemostat cultivation allows the identification of carbon source responses without interference by carbon-catabolite repression, accumulation of toxic products, and differences in specific growth rate. This study focuses on the yeast Saccharomyces cerevisiae, grown in aerobic, carbon-limited chemostat cultures. Genome-wide transcript levels and in vivo fluxes were compared for growth on two sugars, glucose and maltose, and for two C2-compounds, ethanol and acetate. In contrast to previous reports on batch cultures, few genes (180 genes) responded to changes of the carbon source by a changed transcript level. Very few transcript levels were changed when glucose as the growth-limiting nutrient was compared with maltose (33 transcripts), or when acetate was compared with ethanol (16 transcripts). Although metabolic flux analysis using a stoichiometric model revealed major changes in the central carbon metabolism, only 117 genes exhibited a significantly different transcript level when sugars and C2-compounds were provided as the growth-limiting nutrient. Despite the extensive knowledge on carbon source regulation in yeast, many of the carbon source-responsive genes encoded proteins with unknown or incompletely characterized biological functions. In silico promoter analysis of carbon source-responsive genes confirmed the involvement of several known transcriptional regulators and suggested the involvement of additional regulators. Transcripts involved in the glyoxylate cycle and gluconeogenesis showed a good correlation with in vivo fluxes. This correlation was, however, not observed for other important pathways, including the pentose-phosphate pathway, tricarboxylic acid cycle, and, in particular, glycolysis. These results indicate that in vivo fluxes in the central carbon metabolism of S. cerevisiae grown in steadystate, carbon-limited chemostat cultures are controlled to a large extent via post-transcriptional mechanisms.

The extent to which ATP demand controls the glycolytic flux depends strongly on the organism and conditions for growth

Using molecular genetics we have introduced uncoupled ATPase activity in two different bacterial species, Escherichia coli and Lactococcus lactis, and determined the elasticities of the growth rate and glycolytic flux towards the intracellular [ATP]/[ADP] ratio. During balanced growth in batch cultures of E. coli the ATP demand was found to have almost full control on the glycolytic flux (FCC=0.96) and the flux could be stimulated by 70%. In contrast to this, in L. lactis the control by ATP demand on the glycolytic flux was close to zero. However, when we used non-growing cells of L. lactis (which have a low glycolytic flux) the ATP demand had a high flux control and the flux could be stimulated more than two fold. We suggest that the extent to which ATP demand controls the glycolytic flux depends on how much excess capacity of glycolysis is present in the cells.

Genome-wide analysis of the response to cell wall mutations in the yeast Saccharomyces cerevisiae

Perturbations of the yeast cell wall trigger a repair mechanism that reconfigures its molecular structure to preserve cell integrity. To investigate this mechanism, we compared the global gene expression in five mutant strains, each bearing a mutation (i.e. fks1, kre6, mnn9, gas1, and knr4 mutants) that affects in a different manner the cell wall construction. Altogether, 300 responsive genes were kept based on high stringency criteria during data processing. Functional classification of these differentially expressed genes showed a substantial subset of induced genes involved in cell wall construction and an enrichment of metabolic, energy generation, and cell defense categories, whereas families of genes belonging to transcription, protein synthesis, and cellular growth were underrepresented. Clustering methods isolated a single group of approximately 80 up-regulated genes that could be considered as the stereotypical transcriptional response of the cell wall compensatory mechanism. The in silico analysis of the DNA upstream region of these co-regulated genes revealed pairwise combinations of DNA-binding sites for transcriptional factors implicated in stress and heat shock responses (Msn2/4p and Hsf1p) with Rlm1p and Swi4p, two PKC1-regulated transcription factors involved in the activation genes related to cell wall biogenesis and G1/S transition. Moreover, this computational analysis also uncovered the 6-bp 5'-AGCCTC-3' CDRE (calcineurin-dependent response element) motif in 40% of the co-regulated genes. This motif was recently shown to be the DNA binding site for Crz1p, the major effector of calcineurin-regulated gene expression in yeast. Taken altogether, the data presented here lead to the conclusion that the cell wall compensatory mechanism, as triggered by cell wall mutations, integrates three major regulatory systems: namely the PKC1-SLT2 mitogen-activated protein kinase-signaling module, the "global stress" response mediated by Msn2/4p, and the Ca2+/calcineurin-dependent pathway. The relative importance of these regulatory systems in the cell wall compensatory mechanism is discussed.

Schemes of flux control in a model of Saccharomyces cerevisiae glycolysis

We used parameter scanning to emulate changes to the limiting rate for steps in a fitted model of glucose-derepressed yeast glycolysis. Three flux-control regimes were observed, two of which were under the dominant control of hexose transport, in accordance with various experimental studies and other model predictions. A third control regime in which phosphofructokinase exerted dominant glycolytic flux control was also found, but it appeared to be physiologically unreachable by this model, and all realistically obtainable flux control regimes featured hexose transport as a step involving high flux control.

Genetic manipulation of 6-phosphofructo-1-kinase and fructose 2,6-bisphosphate levels affects the extent to which benzoic acid inhibits the growth of Saccharomyces cerevisiae

The mechanisms by which the weak acid preservative benzoic acid inhibits the growth of Saccharomyces cerevisiae have been investigated. A reduction in the pyruvate kinase level, which decreases glycolytic flux, did not increase the sensitivity of yeast to benzoic acid. However, a decrease in 6-phosphofructo-1-kinase (PF1K), which does not affect glycolytic flux, did increase sensitivity to benzoic acid. Also, resistance was increased by elevating PF1K levels. Hence, resistance to benzoic acid was not dependent upon optimum glycolytic flux, but upon an adequate PF1K activity. Benzoic acid was shown to depress fructose 2,6-bisphosphate levels in YKC14, a mutant with low PF1K levels. This effect was partially suppressed by overexpressing constitutively active 6-phosphofructo-2-kinase (Pfk26(Asp644)) or by inactivating fructose-2,6-bisphosphatase (in a Deltafbp26 mutant). The inactivation of PF2K (in a Deltapfk26 Deltapfk27 mutant) increased benzoic acid sensitivity. Therefore, the antimicrobial effects of benzoic acid can be relieved, at least in part, by the genetic manipulation of PF1K or fructose 2,6-bisphosphate levels.

Simultaneous overexpression of enzymes of the lower part of glycolysis can enhance the fermentative capacity of Saccharomyces cerevisiae

Recombinant S. cerevisiae strains, with elevated levels of the enzymes of lower glycolysis (glyceraldehyde-3-phosphate dehydrogenase, phosphoglycerate mutase, phosphoglycerate kinase, enolase, pyruvate kinase, pyruvate decarboxylase and alcohol dehydrogenase) were physiologically characterized. During growth on glucose the enzyme levels in the recombinant strains (YHM4 and YHM7) were 1.1-3.4-fold higher than in the host strain (CEN.PK.K45). The recombinant strains were grown in aerobic or anaerobic batch cultures on glucose or a mixture of glucose and galactose. The specific ethanol production rates in the recombinant strains were the same as for the host strain and the physiological behaviour of the recombinant strains and the host strain was similar. When the cellular demand for ATP was increased by means of glucose pulses (final concentrations of 3.9 g/l or 2.0 g/l, respectively) to aerobic chemostat cultures maintained at a dilution rate of 0.08/h, the specific carbon dioxide production rate (qCO(2)) of CEN.PK.K45 accelerated at 6x10(-3) mmol/g/min(2) during the first 15 min, whereas during the same time period the qCO(2) of YHM7 accelerated twice as fast at 12x10(-3) mmol/g/min(2), indicating a higher fermentative capacity in the recombinant strain.

Reconstitution of ethanolic fermentation in permeabilized spheroplasts of wild-type and trehalose-6-phosphate synthase mutants of the yeast Saccharomyces cerevisiae

In the yeast Saccharomyces cerevisiae, TPS1-encoded trehalose-6-phosphate synthase (TPS) exerts an essential control on the influx of glucose into glycolysis, presumably by restricting hexokinase activity. Deletion of TPS1 results in severe hyperaccumulation of sugar phosphates and near absence of ethanol formation. To investigate whether trehalose 6-phosphate (Tre6P) is the sole mediator of hexokinase inhibition, we have reconstituted ethanolic fermentation from glucose in permeabilized spheroplasts of the wild-type, tps1Delta and tps2Delta (Tre6P phosphatase) strains. For the tps1Delta strain, ethanol production was significantly lower and was associated with hyperaccumulation of Glu6P and Fru6P. A tps2Delta strain shows reduced accumulation of Glu6P and Fru6P both in intact cells and in permeabilized spheroplasts. These results are not consistent with Tre6P being the sole mediator of hexokinase inhibition. Reconstitution of ethanolic fermentation in permeabilized spheroplasts with glycolytic intermediates indicates additional target site(s) for the Tps1 control. Addition of Tre6P partially shifts the ethanol production rate and the metabolite pattern in permeabilized tps1Delta spheroplasts to those of the wild-type strain, but only with glucose as substrate. This is observed at a very high ratio of glucose to Tre6P. Inhibition of hexokinase activity by Tre6P is less efficiently counteracted by glucose in permeabilized spheroplasts compared to cell extracts, and this effect is largely abolished by deletion of TPS2 but not TPS1. In permeabilized spheroplasts, hexokinase activity is significantly lower in a tps2Delta strain compared to a wild-type strain and this difference is strongly reduced by additional deletion of TPS1. These results indicate that Tps1-mediated protein-protein interactions are important for control of glucose influx into yeast glycolysis, that Tre6P inhibition of hexokinase might not be competitive with respect to glucose in vivo and that also Tps2 appears to play a role in the control of hexokinase activity.

The importance of ATP as a regulator of glycolytic flux in Saccharomyces cerevisiae

The control of glycolytic flux in the yeast Saccharomyces cerevisiae was studied by using permeabilized cells. Cells were harvested from chemostat cultures and, after removal of the cell wall, nystatin was used to permeabilize the spheroplasts. By this method it is possible to study the performance and regulation of a complete and functional metabolic pathway and not only a single enzymatic step. The results showed that ATP has a strong negative effect on glycolytic activity affecting several of the glycolytic enzymes. However, the main targets for ATP inhibition was phosphofructokinase and pyruvate kinase. Phospofructokinase was inhibited by ATP concentrations starting at about 1-2 mM, while pyruvate kinase required ATP levels above 2.5 mM before any inhibition was visible. These ATP concentrations were in the same range as measured for nitrogen- and glucose-limited cells cultivated in chemostat cultures. Other potential candidates as enzymes susceptible to ATP inhibition included hexokinase and enolase. The ATP:ADP ratio, as well as trehalose-6-phosphate levels, did not seem to influence the glycolytic activity.

Simultaneous genomic overexpression of seven glycolytic enzymes in the yeast Saccharomyces cerevisiae

Fusions of the glycolytic genes TPI1, PGK1, ENO1, PYK1, PDC1, and ADH1 with the lacZ reporter gene of Escherichia coli and a lacZ fusion construct of a 390-bp fragment from the promoter of the HXT7 gene were assayed for beta-galactosidase activity. The glycolytic promoters were induced after addition of glucose to ethanol-grown cells, whereas the HXT7 promoter fragment showed a constitutive beta-galactosidase expression on both carbon sources. The genes coding for the seven enzymes of lower glycolysis Tdh, Pgk, Gpm, Eno, Pyk, Pdc, and Adh were simultaneously put under the control of the same strong promoter, a truncated HXT7 promoter that is constitutively active on ethanol as well as on glucose medium. Genomic expression of the glycolytic genes under the control of this promoter, resulted in an at least 2-fold overexpression. The gene MSG5 was isolated, coding for a protein phosphatase normally involved in cell cycle regulation, as a factor that possibly influences the expression of the HXT7 gene. However, overexpression of MSG5 had no effect on the expression of the HXT7/lacZ fusion, whereas a deletion of this gene resulted in a decreased expression of beta-galactosidase.

A potential role of the cytoskeleton of Saccharomyces cerevisiae in a functional organization of glycolytic enzymes

Numerous individual enzymes participate in a given synthetic or degradative pathway in which the product of one reaction becomes the substrate for the subsequent enzyme. This raises the question of whether the product of one 'soluble' enzyme diffuses freely through the available cell volume, where it accidentally collides with the subsequent 'soluble' enzyme. Alternatively, enzymes acting in a given pathway may be organized in ordered structures, metabolons. Certain glycolytic enzymes have been shown to co-localize with the cytoskeleton in mammalian cells. We deleted genes coding for proteins associated with the cytoskeleton of Saccharomyces cerevisiae: TPM1 coding for tropomyosin, SAC6 for fimbrin and CIN1 for a microtubule-associated protein. Single deletions or deletions of two such genes had no effect on the specific activities of glycolytic enzymes, or on the rates of glucose consumption and ethanol production. However, the concentrations of glycolytic metabolites during a switch from a gluconeogenic mode of metabolism, growth on an ethanol medium, to glycolysis after glucose addition showed transient deviations from the normal change in metabolite concentrations, as observed in wild type cells. However, all metabolites in mutant strains reached wild-type levels within 2-4 h after the shift. Only ATP levels remained low in all but the tmp1-Delta-sac6-Delta double mutant strains. These observations can be interpreted to mean that metabolic reorganization from a gluconeogenic to a glycolytic metabolism is facilitated by an intact cytoskeleton in yeast.

Fructose-2,6-bisphosphate and control of carbohydrate metabolism in eukaryotes

Fructose-2,6-bisphosphate is an important intracellular biofactor in the control of carbohydrate metabolic fluxes in eukaryotes. It is generated from ATP and fructose-6-phosphate by 6-phosphofructo-2-kinase and degraded to fructose-6-phosphate and phosphate ion by fructose-2,6-bisphosphatase. In most organisms these enzymatic activities are contained in a single polypeptide. The reciprocal modulation of the kinase and bisphosphatase activities by post-translational modifications places the level of the biofactor under the control of extra-cellular signals. In general, these signals are generated in response to changing nutritional states, therefore, fructose-2,6-bisphosphate plays a role in the adaptation of organisms, and the tissues within them, to changes in environmental and metabolic states. Although the specific mechanism of fructose-2,6-bisphosphate action varies between species and between tissues, most involve the allosteric activation of 6-phosphofructo-1-kinase and inhibition of fructose-1,6-bisphosphatase. These highly conserved enzymes regulate the fructose-6-phosphate/fructose-1,6-bisphosphate cycle, and thereby, determine the carbon flux. It is by reciprocal modulation of these activities that fructose-2,6-bisphosphate plays a fundamental role in eukaryotic carbohydrate metabolism.

Starting up yeast glycolysis

Successfully igniting the yeast glycolytic flux during the transition from gluconeogenic to fermentative growth seems to be a matter of balance and coordination between a multitude of events. The contours of the sugar sensing and signalling pathways that regulate this transition are only beginning to emerge.