Dihydroxyacetone kinases in Saccharomyces cerevisiae are involved in detoxification of dihydroxyacetone

Identification of propranolol and derivatives that are chemical inhibitors of phosphatidate phosphatase as potential broad-spectrum fungicides

Plant diseases cause enormous economic losses in agriculture and threaten global food security, and application of agrochemicals is an important method of crop disease control. Exploration of disease-resistance mechanisms and synthesis of highly bioactive agrochemicals are thus important research objectives. Here, we show that propranolol, a phosphatidate phosphatase (Pah) inhibitor, effectively suppresses fungal growth, sporulation, sexual reproduction, and infection of diverse plants. The MoPah1 enzyme activity of the rice blast fungus Magnaporthe oryzae is inhibited by propranolol. Alterations in lipid metabolism are associated with inhibited hyphal growth and appressorium formation caused by propranolol in M. oryzae. Propranolol inhibits a broad spectrum of 12 plant pathogens, effectively inhibiting infection of barley, wheat, maize, tomato, and pear. To improve antifungal capacity, we synthesized a series of propranolol derivatives, one of which shows a 16-fold increase in antifungal ability and binds directly to MoPah1. Propranolol and its derivatives can also reduce the severity of rice blast and Fusarium head blight of wheat in the field. Taken together, our results demonstrate that propranolol suppresses fungal development and infection through mechanisms involved in lipid metabolism. Propranolol and its derivatives may therefore be promising candidates for fungicide development.

Yeast osmoregulation - glycerol still in pole position

In response to osmotic dehydration cells sense, signal, alter gene expression, and metabolically counterbalance osmotic differences. The main compatible solute/osmolyte that accumulates in yeast cells is glycerol, which is produced from the glycolytic intermediate dihydroxyacetone phosphate. This review covers recent advancements in understanding mechanisms involved in sensing, signaling, cell-cycle delays, transcriptional responses as well as post-translational modifications on key proteins in osmoregulation. The protein kinase Hog1 is a key-player in many of these events, however, there is also a growing body of evidence for important Hog1-independent mechanisms playing vital roles. Several missing links in our understanding of osmoregulation will be discussed and future avenues for research proposed. The review highlights that this rather simple experimental system—salt/sorbitol and yeast—has developed into an enormously potent model system unravelling important fundamental aspects in biology.

Formate Dehydrogenase Improves the Resistance to Formic Acid and Acetic Acid Simultaneously in Saccharomyces cerevisiae

Bioethanol from lignocellulosic biomass is a promising and sustainable strategy to meet the energy demand and to be carbon neutral. Nevertheless, the damage of lignocellulose-derived inhibitors to microorganisms is still the main bottleneck. Developing robust strains is critical for lignocellulosic ethanol production. An evolved strain with a stronger tolerance to formate and acetate was obtained after adaptive laboratory evolution (ALE) in the formate. Transcriptional analysis was conducted to reveal the possible resistance mechanisms to weak acids, and fdh coding for formate dehydrogenase was selected as the target to verify whether it was related to resistance enhancement in Saccharomyces cerevisiae F3. Engineered S. cerevisiae FA with fdh overexpression exhibited boosted tolerance to both formate and acetate, but the resistance mechanism to formate and acetate was different. When formate exists, it breaks down by formate dehydrogenase into carbon dioxide (CO₂) to relieve its inhibition. When there was acetate without formate, FDH1 converted CO₂ from glucose fermentation to formate and ATP and enhanced cell viability. Together, fdh overexpression alone can improve the tolerance to both formate and acetate with a higher cell viability and ATP, which provides a novel strategy for robustness strain construction to produce lignocellulosic ethanol.

Physiological and Molecular Characterization of an Oxidative Stress-Resistant Saccharomyces cerevisiae Strain Obtained by Evolutionary Engineering

Oxidative stress is a major stress type observed in yeast bioprocesses, resulting in a decrease in yeast growth, viability, and productivity. Thus, robust yeast strains with increased resistance to oxidative stress are in highly demand by the industry. In addition, oxidative stress is also associated with aging and age-related complex conditions such as cancer and neurodegenerative diseases. Saccharomyces cerevisiae, as a model eukaryote, has been used to study these complex eukaryotic processes. However, the molecular mechanisms underlying oxidative stress responses and resistance are unclear. In this study, we have employed evolutionary engineering (also known as adaptive laboratory evolution – ALE) strategies to obtain an oxidative stress-resistant and genetically stable S. cerevisiae strain. Comparative physiological, transcriptomic, and genomic analyses of the evolved strain were then performed with respect to the reference strain. The results show that the oxidative stress-resistant evolved strain was also cross-resistant against other types of stressors, including heat, freeze-thaw, ethanol, cobalt, iron, and salt. It was also found to have higher levels of trehalose and glycogen production. Further, comparative transcriptomic analysis showed an upregulation of many genes associated with the stress response, transport, carbohydrate, lipid and cofactor metabolic processes, protein phosphorylation, cell wall organization, and biogenesis. Genes that were downregulated included those related to ribosome and RNA processing, nuclear transport, tRNA, and cell cycle. Whole genome re-sequencing analysis of the evolved strain identified mutations in genes related to the stress response, cell wall organization, carbohydrate metabolism/transport, which are in line with the physiological and transcriptomic results, and may give insight toward the complex molecular mechanisms of oxidative stress resistance.

The high osmolarity glycerol (HOG) pathway in fungi†

While fungi are widely occupying nature, many species are responsible for devastating mycosis in humans. Such niche diversity explains how quick fungal adaptation is necessary to endow the capacity of withstanding fluctuating environments and to cope with host-imposed conditions. Among all the molecular mechanisms evolved by fungi, the most studied one is the activation of the phosphorelay signalling pathways, of which the high osmolarity glycerol (HOG) pathway constitutes one of the key molecular apparatus underpinning fungal adaptation and virulence. In this review, we summarize the seminal knowledge of the HOG pathway with its more recent developments. We specifically described the HOG-mediated stress adaptation, with a particular focus on osmotic and oxidative stress, and point out some lags in our understanding of its involvement in the virulence of pathogenic species including, the medically important fungi Candida albicans, Cryptococcus neoformans, and Aspergillus fumigatus, compared to the model yeast Saccharomyces cerevisiae. Finally, we also highlighted some possible applications of the HOG pathway modifications to improve the fungal-based production of natural products in the industry.

Identification and biochemical characterization of an ATP-dependent dihydroxyacetone kinase from Trypanosoma cruzi

Dihydroxyacetone (DHA) can be used as an energy source by many cell types; however, it is toxic at high concentrations. The enzyme dihydroxyacetone kinase (DAK) has shown to be involved in DHA detoxification and osmoregulation. Among protozoa of the genus Trypanosoma, T. brucei, which causes sleeping sickness, is highly sensitive to DHA and does not have orthologous genes to DAK. Conversely, T. cruzi, the etiological agent of Chagas Disease, has two putative ATP-dependent DAK (TcDAKs) sequences in its genome. Here we show that T. cruzi epimastigote lysates present a DAK specific activity of 27.1 nmol/min/mg of protein and that this form of the parasite is able to grow in the presence of 2 mM DHA. TcDAK gene was cloned and the recombinant enzyme (recTcDAK) was expressed in Escherichia coli. An anti-recTcDAK serum reacted with a protein of the expected molecular mass of 61 kDa in epimastigotes. recTcDAK presented maximal activity using Mg⁺², showing a Km of 6.5 μM for DHA and a K_0.5 of 124.7 μM for ATP. As it was reported for other DAKs, recTcDAK activity was inhibited by FAD with an IC₅₀ value of 0.33 mM. In conclusion, TcDAK is the first DAK described in trypanosomatids confirming another divergent metabolism between T. brucei and T. cruzi.

Efficient Conversion of Glycerol to Ethanol by an Engineered Saccharomyces cerevisiae Strain

Glycerol is an eco-friendly solvent that enhances plant biomass decomposition via glycerolysis in many pretreatment methods. Nonetheless, inefficient conversion of glycerol to ethanol by natural Saccharomyces cerevisiae limits its use in these processes. In this study, we have developed an efficient glycerol-converting yeast strain by genetically modifying the oxidation of cytosolic NAD (NADH) by an O₂-dependent dynamic shuttle and abolishing both glycerol phosphorylation and biosynthesis in S. cerevisiae strain D452-2, as well as by vigorous expression of whole genes in the dihydroxyacetone (DHA) pathway (Candida utilis glycerol facilitator, Ogataea polymorpha glycerol dehydrogenase, endogenous dihydroxyacetone kinase, and triosephosphate isomerase). The engineered strain showed conversion efficiencies (CE) up to 0.49 g ethanol/g glycerol (98% of theoretical CE), with a production rate of >1 g liter^-1 h^-1 when glycerol was supplemented in a single fed-batch fermentation in a rich medium. Furthermore, the engineered strain converted a mixture of glycerol and glucose into bioethanol (>86 g/liter) with 92.8% CE. To the best of our knowledge, this is the highest reported titer of bioethanol produced from glycerol and glucose. Notably, we developed a glycerol-utilizing transformant from a parent strain which cannot utilize glycerol as a sole carbon source. The developed strain converted glycerol to ethanol with a productivity of 0.44 g liter^-1 h^-1 on minimal medium under semiaerobic conditions. Our findings will promote the utilization of glycerol in eco-friendly biorefineries and integrate bioethanol and plant oil industries. IMPORTANCE With the development of efficient lignocellulosic biorefineries, glycerol has attracted attention as an eco-friendly biomass-derived solvent that can enhance the dissociation of lignin and cell wall polysaccharides during the pretreatment process. Coconversion of glycerol with the sugars released from biomass after glycerolysis increases the resources for ethanol production and lowers the burden of component separation. However, low conversion efficiency from glycerol and sugars limits the industrial application of this process. Therefore, the generation of an efficient glycerol-fermenting yeast will promote the applicability of integrated biorefineries. Hence, metabolic flux control in yeast grown on glycerol will lead to the generation of cell factories that produce chemicals, which will boost biodiesel and bioethanol industries. Additionally, the use of glycerol-fermenting yeast will reduce global warming and generation of agricultural waste, leading to the establishment of a sustainable society.

Cell-free chemoenzymatic starch synthesis from carbon dioxide

[Figure: see text].

Strategies and challenges with the microbial conversion of methanol to high-value chemicals

As alternatives to traditional fermentation substrates, methanol (CH₃ OH), carbon dioxide (CO₂ ) and methane (CH₄ ) represent promising one-carbon (C1) sources that are readily available at low-cost and share similar metabolic pathway. Of these C1 compounds, methanol is used as a carbon and energy source by native methylotrophs, and can be obtained from CO₂ and CH₄ by chemical catalysis. Therefore, constructing and rewiring methanol utilization pathways may enable the use of one-carbon sources for microbial fermentations. Recent bioengineering efforts have shown that both native and nonnative methylotrophic organisms can be engineered to convert methanol, together with other carbon sources, into biofuels and other commodity chemicals. However, many challenges remain and must be overcome before industrial-scale bioprocessing can be established using these engineered cell refineries. Here, we provide a comprehensive summary and comparison of methanol metabolic pathways from different methylotrophs, followed by a review of recent progress in engineering methanol metabolic pathways in vitro and in vivo to produce chemicals. We discuss the major challenges associated with establishing efficient methanol metabolic pathways in microbial cells, and propose improved designs for future engineering.

Autosomal recessive hypotrichosis with loose anagen hairs associated with TKFC mutations

BACKGROUND

Loose anagen hair is a rare form of impaired hair anchorage in which anagen hairs that lack inner and outer root sheaths can be gently and painlessly plucked from the scalp. This condition usually occurs in children and is often self-limiting. A genetic basis for the disorder has been suggested but not proven. A better understanding the aetiology of loose anagen hair may improve prevention and treatment strategies.

OBJECTIVES

To identify a possible genetic basis of loose anagen hair using next-generation DNA sequencing and functional analysis of variants identified.

METHODS

In this case study, whole-exome sequencing analysis of a pedigree with one affected individual with features of loose anagen hair was performed.

RESULTS

The patient was found to be compound heterozygous for two single-nucleotide substitutions in TKFC resulting in the following missense mutations: c.574G> C (p.Gly192Arg) and c.682C> T (p.Arg228Trp). Structural analysis of human TKFC showed that both mutations are located near the active site cavity. Kinetic assays of recombinant proteins bearing either of these amino acid substitutions showed almost no dihydroxyacetone kinase or D-glyceraldehyde kinase activity, and FMN cyclase activity reduced to just 10% of wildtype catalytic activity.

CONCLUSIONS

TKFC missense mutations may predispose to the development of loose anagen hairs. Identification of this new biochemical pathobiology expands the metabolic and genetic basis of hypotrichosis.

Hog1-induced transcription of RTC3 and HSP12 is robust and occurs in cells lacking Msn2, Msn4, Hot1 and Sko1

The yeast MAP kinase Hog1 pathway activates transcription of several hundreds genes. Large-scale gene expression and DNA binding assays suggest that most Hog1-induced genes are regulated by the transcriptional activators Msn2/4, Hot1 and Sko1. These studies also revealed the target genes of each activator and the putative binding sites on their promoters. In a previous study we identified a group of genes, which we considered the bona fide targets of Hog1, because they were induced in response to expression of intrinsically active mutant of Hog1, in the absence of any stress. We previously analyzed the promoter of the most highly induced gene, STL1, and noticed that some promoter properties were different from those proposed by large-scale data. We therefore continue to study promoters individually and present here analyses of promoters of more Hog1's targets, RTC3, HSP12, DAK1 and ALD3. We report that RTC3 and HSP12 promoters are robust and are induced, to different degrees, even in cells lacking all four activators. DAK1 and ALD3 promoters are not robust and fully depend on a single activator, DAK1 on Sko1 and ALD3 on Msn2/4. Most of these observations could not be inferred from the large-scale data. Msn2/4 are involved in regulating all four promoters. It was assumed, therefore, that the promoters are spontaneously active in ras2Δ cells, in which Msn2/4 are known to be de-repressed. Intriguingly, the promoters were not active in BY4741ras2Δ cells, but were de-repressed, as expected, in ras2Δ cells of other genetic backgrounds. This study describes two phenomena. One, some Hog1's target promoters are most robust, backupped by many activators. Second, in contrast to most laboratory strains, the widely used BY4741 strain does not induce Msn2/4 activity when the Ras/cAMP cascade is downregulated.

Bi-allelic Variants in TKFC Encoding Triokinase/FMN Cyclase Are Associated with Cataracts and Multisystem Disease

We report an inborn error of metabolism caused by TKFC deficiency in two unrelated families. Rapid trio genome sequencing in family 1 and exome sequencing in family 2 excluded known genetic etiologies, and further variant analysis identified rare homozygous variants in TKFC. TKFC encodes a bifunctional enzyme involved in fructose metabolism through its glyceraldehyde kinase activity and in the generation of riboflavin cyclic 4',5'-phosphate (cyclic FMN) through an FMN lyase domain. The TKFC homozygous variants reported here are located within the FMN lyase domain. Functional assays in yeast support the deleterious effect of these variants on protein function. Shared phenotypes between affected individuals with TKFC deficiency include cataracts and developmental delay, associated with cerebellar hypoplasia in one case. Further complications observed in two affected individuals included liver dysfunction and microcytic anemia, while one had fatal cardiomyopathy with lactic acidosis following a febrile illness. We postulate that deficiency of TKFC causes disruption of endogenous fructose metabolism leading to generation of by-products that can cause cataract. In line with this, an affected individual had mildly elevated urinary galactitol, which has been linked to cataract development in the galactosemias. Further, in light of a previously reported role of TKFC in regulating innate antiviral immunity through suppression of MDA5, we speculate that deficiency of TKFC leads to impaired innate immunity in response to viral illness, which may explain the fatal illness observed in the most severely affected individual.

Elevated aldolase 1A, retrogene 1 expression induces cardiac apoptosis in rat experimental autoimmune myocarditis model

Acute myocarditis is an unpredictable heart disease that is caused by inflammation-associated cell death. Although viral infection and drug exposure are known to induce acute myocarditis, the molecular basis for its development remains undefined. Using proteomics and molecular analyses in myosin-induced rat experimental autoimmune myocarditis (EAM), we identified that elevated expression of aldolase 1A, retrogene 1 (Aldoart1) is critical to induce mitochondrial dysfunction and acute myocarditis development. Here, we demonstrate that cardiac cell death is associated with increased expressions of proapoptotic genes in addition to high levels of glucose, lactate, and triglyceride in metabolite profiling. The functional protein association network analysis also suggests that Aldoart1 upregulation correlates with high levels of dihydroxyacetone kinase and triglyceride. In H9c2 cardiac cells, lipopolysaccharides (LPS) or high glucose exposure significantly increases the cytochrome c release and the conversion of pro-caspase 3 into the cleaved form of caspase 3. We also found that LPS- or glucose-induced toxicities are almost completely reversed by siRNA-mediated knockdown of Aldoartl, which consequently increases cell viability. Together, our study strongly suggests that Aldoart1 may be involved in inducing mitochondrial apoptotic processes and can be a novel therapeutic target to prevent the onset of acute myocarditis or cardiac apoptosis.

A comprehensive genome-scale model for Rhodosporidium toruloides IFO0880 accounting for functional genomics and phenotypic data

Rhodosporidium toruloides is a red, basidiomycetes yeast that can accumulate a large amount of lipids and produce carotenoids. To better assess this non-model yeast's metabolic capabilities, we reconstructed a genome-scale model of R. toruloides IFO0880's metabolic network (iRhto1108) accounting for 2204 reactions, 1985 metabolites and 1108 genes. In this work, we integrated and supplemented the current knowledge with in-house generated biomass composition and experimental measurements pertaining to the organism's metabolic capabilities. Predictions of genotype-phenotype relations were improved through manual curation of gene-protein-reaction rules for 543 reactions leading to correct recapitulations of 84.5% of gene essentiality data (sensitivity of 94.3% and specificity of 53.8%). Organism-specific macromolecular composition and ATP maintenance requirements were experimentally measured for two separate growth conditions: (i) carbon and (ii) nitrogen limitations. Overall, iRhto1108 reproduced R. toruloides's utilization capabilities for 18 alternate substrates, matched measured wild-type growth yield, and recapitulated the viability of 772 out of 819 deletion mutants. As a demonstration to the model's fidelity in guiding engineering interventions, the OptForce procedure was applied on iRhto1108 for triacylglycerol overproduction. Suggested interventions recapitulated many of the previous successful implementations of genetic modifications and put forth a few new ones.

Beneficial mutations for carotenoid production identified from laboratory-evolved Saccharomyces cerevisiae

Adaptive laboratory evolution (ALE) is a powerful tool used to increase strain fitness in the presence of environmental stressors. If production and strain fitness can be coupled, ALE can be used to increase product formation. In earlier work, carotenoids hyperproducing mutants were obtained using an ALE strategy. Here, de novo mutations were identified in hyperproducers, and reconstructed mutants were explored to determine the exact impact of each mutation on production and tolerance. A single mutation in YMRCTy1-3 conferred increased carotenoid production, and when combined with other beneficial mutations led to further increased β-carotene production. Findings also suggest that the ALE strategy selected for mutations that confer increased carotenoid production as primary phenotype. Raman spectroscopy analysis and total lipid quantification revealed positive correlation between increased lipid content and increased β-carotene production. Finally, we demonstrated that the best combinations of mutations identified for β-carotene production were also beneficial for production of lycopene.

Chemical and Metabolic Controls on Dihydroxyacetone Metabolism Lead to Suboptimal Growth of Escherichia coli

DHA is an attractive triose molecule with a wide range of applications, notably in cosmetics and the food and pharmaceutical industries. DHA is found in many species, from microorganisms to humans, and can be used by Escherichia coli as a growth substrate. However, knowledge about the mechanisms and regulation of this process is currently lacking, motivating our investigation of DHA metabolism in E. coli. We show that under aerobic conditions, E. coli growth on DHA is far from optimal and is hindered by chemical, hierarchical, and possibly allosteric constraints. We show that optimal growth on DHA can be restored by releasing the hierarchical constraint. These results improve our understanding of DHA metabolism and are likely to help unlock biotechnological applications involving DHA as an intermediate, such as the bioconversion of glycerol or C₁ substrates into value-added chemicals.

In this work, we shed light on the metabolism of dihydroxyacetone (DHA), a versatile, ubiquitous, and important intermediate for various chemicals in industry, by analyzing its metabolism at the system level in Escherichia coli. Using constraint-based modeling, we show that the growth of E. coli on DHA is suboptimal and identify the potential causes. Nuclear magnetic resonance analysis shows that DHA is degraded nonenzymatically into substrates known to be unfavorable to high growth rates. Transcriptomic analysis reveals that DHA promotes genes involved in biofilm formation, which may reduce the bacterial growth rate. Functional analysis of the genes involved in DHA metabolism proves that under the aerobic conditions used in this study, DHA is mainly assimilated via the dihydroxyacetone kinase pathway. In addition, these results show that the alternative routes of DHA assimilation (i.e., the glycerol and fructose-6-phosphate aldolase pathways) are not fully activated under our conditions because of anaerobically mediated hierarchical control. These pathways are therefore certainly unable to sustain fluxes as high as the ones predicted in silico for optimal aerobic growth on DHA. Overexpressing some of the genes in these pathways releases these constraints and restores the predicted optimal growth on DHA.

IMPORTANCE DHA is an attractive triose molecule with a wide range of applications, notably in cosmetics and the food and pharmaceutical industries. DHA is found in many species, from microorganisms to humans, and can be used by Escherichia coli as a growth substrate. However, knowledge about the mechanisms and regulation of this process is currently lacking, motivating our investigation of DHA metabolism in E. coli. We show that under aerobic conditions, E. coli growth on DHA is far from optimal and is hindered by chemical, hierarchical, and possibly allosteric constraints. We show that optimal growth on DHA can be restored by releasing the hierarchical constraint. These results improve our understanding of DHA metabolism and are likely to help unlock biotechnological applications involving DHA as an intermediate, such as the bioconversion of glycerol or C₁ substrates into value-added chemicals.

PHO15 genes of Candida albicans and Candida parapsilosis encode HAD-type phosphatases dephosphorylating 2-phosphoglycolate

Most of the phosphatases of human fungal pathogens Candida albicans and C. parapsilosis have never been experimentally characterised, although dephosphorylation reactions are central to many biological processes. PHO15 genes of these yeasts have been annotated as the sequences encoding 4-nitrophenyl phosphatase, on the basis of homology to PHO13 gene from the bakers' yeast Saccharomyces cerevisiae. To examine the real function of these potential phosphatases from Candida spp., CaPho15p and CpPho15p were prepared using expression in Escherichia coli and characterised. They share the hallmark motifs of the haloacid dehalogenase superfamily, readily hydrolyse 4-nitrophenyl phosphate at pH 8-8.3 and require divalent cations (Mg2+, Mn2+ or Co2+) as cofactors. CaPho15p and CpPho15p did not dephosphorylate phosphopeptides, but rather hydrolysed molecules related to carbohydrate metabolism. The preferred substrate for the both phosphatases was 2-phosphoglycolate. Among the other molecules tested, CaPho15 showed preference for glyceraldehyde phosphate and ß-glycerol phosphate, while CpPho15 dephosphorylated mainly 1,3-dihydroxyacetone phosphate. This type of substrate specificity indicates that CaPho15 and CpPho15 may be a part of metabolic repair system of C. albicans and C. parapsilosis.

Closure of the Human TKFC Active Site: Comparison of the Apoenzyme and the Complexes Formed with Either Triokinase or FMN Cyclase Substrates

Human triokinase/flavin mononucleotide (FMN) cyclase (hTKFC) catalyzes the adenosine triphosphate (ATP)-dependent phosphorylation of D-glyceraldehyde and dihydroxyacetone (DHA), and the cyclizing splitting of flavin adenine dinucleotide (FAD). hTKFC structural models are dimers of identical subunits, each with two domains, K and L, with an L2-K1-K2-L1 arrangement. Two active sites lie between L2-K1 and K2-L1, where triose binds K and ATP binds L, although the resulting ATP-to-triose distance is too large (≈14 Å) for phosphoryl transfer. A 75-ns trajectory of molecular dynamics shows considerable, but transient, ATP-to-DHA approximations in the L2-K1 site (4.83 Å or 4.16 Å). To confirm the trend towards site closure, and its relationship to kinase activity, apo-hTKFC, hTKFC:2DHA:2ATP and hTKFC:2FAD models were submitted to normal mode analysis. The trajectory of hTKFC:2DHA:2ATP was extended up to 160 ns, and 120-ns trajectories of apo-hTKFC and hTKFC:2FAD were simulated. The three systems were comparatively analyzed for equal lengths (120 ns) following the principles of essential dynamics, and by estimating site closure by distance measurements. The full trajectory of hTKFC:2DHA:2ATP was searched for in-line orientations and short distances of DHA hydroxymethyl oxygens to ATP γ-phosphorus. Full site closure was reached only in hTKFC:2DHA:2ATP, where conformations compatible with an associative phosphoryl transfer occurred in L2-K1 for significant trajectory time fractions.

Minor Isozymes Tailor Yeast Metabolism to Carbon Availability

Gene duplication is one of the main evolutionary paths to new protein function. Typically, duplicated genes either accumulate mutations and degrade into pseudogenes or are retained and diverge in function. Some duplicated genes, however, show long-term persistence without apparently acquiring new function. An important class of isozymes consists of those that catalyze the same reaction in the same compartment, where knockout of one isozyme causes no known functional defect. Here we present an approach to assigning specific functional roles to seemingly redundant isozymes. First, gene expression data are analyzed computationally to identify conditions under which isozyme expression diverges. Then, knockouts are compared under those conditions. This approach revealed that the expression of many yeast isozymes diverges in response to carbon availability and that carbon source manipulations can induce fitness phenotypes for seemingly redundant isozymes. A driver of these fitness phenotypes is differential allosteric enzyme regulation, indicating isozyme divergence to achieve more-optimal control of metabolism.

Isozymes are enzymes that differ in sequence but catalyze the same chemical reactions. Despite their apparent redundancy, isozymes are often retained over evolutionary time, suggesting that they contribute to fitness. We developed an unsupervised computational method for identifying environmental conditions under which isozymes are likely to make fitness contributions. This method analyzes published gene expression data to find specific experimental perturbations that induce differential isozyme expression. In yeast, we found that isozymes are strongly enriched in the pathways of central carbon metabolism and that many isozyme pairs show anticorrelated expression during the respirofermentative shift. Building on these observations, we assigned function to two minor central carbon isozymes, aconitase 2 (ACO2) and pyruvate kinase 2 (PYK2). ACO2 is expressed during fermentation and proves advantageous when glucose is limiting. PYK2 is expressed during respiration and proves advantageous for growth on three-carbon substrates. PYK2’s deletion can be rescued by expressing the major pyruvate kinase only if that enzyme carries mutations mirroring PYK2’s allosteric regulation. Thus, central carbon isozymes help to optimize allosteric metabolic regulation under a broad range of potential nutrient conditions while requiring only a small number of transcriptional states.

IMPORTANCE Gene duplication is one of the main evolutionary paths to new protein function. Typically, duplicated genes either accumulate mutations and degrade into pseudogenes or are retained and diverge in function. Some duplicated genes, however, show long-term persistence without apparently acquiring new function. An important class of isozymes consists of those that catalyze the same reaction in the same compartment, where knockout of one isozyme causes no known functional defect. Here we present an approach to assigning specific functional roles to seemingly redundant isozymes. First, gene expression data are analyzed computationally to identify conditions under which isozyme expression diverges. Then, knockouts are compared under those conditions. This approach revealed that the expression of many yeast isozymes diverges in response to carbon availability and that carbon source manipulations can induce fitness phenotypes for seemingly redundant isozymes. A driver of these fitness phenotypes is differential allosteric enzyme regulation, indicating isozyme divergence to achieve more-optimal control of metabolism.

Toxicity of dihydroxyacetone is exerted through the formation of methylglyoxal in Saccharomyces cerevisiae: effects on actin polarity and nuclear division

Dihydroxyacetone (DHA) is the smallest ketotriose, and it is utilized by many organisms as an energy source. However, at higher concentrations, DHA becomes toxic towards several organisms including the budding yeast Saccharomyces cerevisiae. In the present study, we show that DHA toxicity is due to its spontaneous conversion to methylglyoxal (MG) within yeast cells. A mutant defective in MG-metabolizing enzymes (glo1Δgre2Δgre3Δ) exhibited higher susceptibility to DHA. Intracellular MG levels increased following the treatment of glo1Δgre2Δgre3Δ cells with DHA. We previously reported that MG depolarized the actin cytoskeleton and changed vacuolar morphology. We herein demonstrated the depolarization of actin and morphological changes in vacuoles following a treatment with DHA. Furthermore, we found that both MG and DHA caused the morphological change in nucleus, and inhibited the nuclear division. Our results suggest that the conversion of DHA to MG is a dominant contributor to its cytotoxicity.

Functional genomics for the oleaginous yeast Yarrowia lipolytica

Oleaginous yeasts are valuable systems for biosustainable production of hydrocarbon-based chemicals. Yarrowia lipolytica is one of the best characterized of these yeast with respect to genome annotation and flux analysis of metabolic processes. Nonetheless, progress is hampered by a dearth of genome-wide tools enabling functional genomics. In order to remedy this deficiency, we developed a library of Y. lipolytica insertion mutants via transposon mutagenesis. The Hermes DNA transposon was expressed to achieve saturation mutagenesis of the genome. Over 534,000 independent insertions were identified by next-generation sequencing. Poisson analysis of insertion density classified ~ 22% of genes as essential. As expected, most essential genes have homologs in Saccharomyces cerevisiae and Schizosaccharomyces pombe, and the majority of those are also essential. As an obligate aerobe, Y. lipolytica has significantly more respiration - related genes that are classified as essential than do S. cerevisiae and S. pombe. Contributions of non-essential genes to growth in glucose and glycerol carbon sources were assessed and used to evaluate two recent genome-scale models of Y. lipolytica metabolism. Fluorescence-activated cell sorting identified mutants in which lipid accumulation is increased. Our findings provide insights into biosynthetic pathways, compartmentalization of enzymes, and distinct functions of paralogs. This functional genomic analysis of the oleaginous yeast Y. lipolytica provides an important resource for modeling, bioengineering, and design of synthetic minimalized strains of respiratory yeasts.

Identification and characterization of EYK1, a key gene for erythritol catabolism in Yarrowia lipolytica

Erythritol is a four-carbon sugar alcohol synthesized by osmophilic yeasts, such as Yarrowia lipolytica, in response to osmotic stress. This metabolite has application as food additive due to its sweetening properties. Although Y. lipolytica can produce erythritol at a high level from glycerol, it is also able to consume it as carbon source. This ability negatively affects erythritol productivity and represents a serious drawback for the development of an efficient erythritol production process. In this study, we have isolated by insertion mutagenesis a Y. lipolytica mutant unable to grow on erythritol. Genomic characterization of the latter highlighted that the mutant phenotype is directly related to the disruption of the YALI0F01606g gene. Several experimental evidences suggested that the identified gene, renamed EYK1, encodes an erythrulose kinase. The mutant strain showed an enhanced capacity to produce erythritol as compared to the wild-type strain. Moreover, in specific experimental conditions, it is also able to convert erythritol to erythrulose, another compound of biotechnological interest.

The yeast osmostress response is carbon source dependent

Adaptation to altered osmotic conditions is a fundamental property of living cells and has been studied in detail in the yeast Saccharomyces cerevisiae. Yeast cells accumulate glycerol as compatible solute, controlled at different levels by the High Osmolarity Glycerol (HOG) response pathway. Up to now, essentially all osmostress studies in yeast have been performed with glucose as carbon and energy source, which is metabolised by glycolysis with glycerol as a by-product. Here we investigated the response of yeast to osmotic stress when yeast is respiring ethanol as carbon and energy source. Remarkably, yeast cells do not accumulate glycerol under these conditions and it appears that trehalose may partly take over the role as compatible solute. The HOG pathway is activated in very much the same way as during growth on glucose and is also required for osmotic adaptation. Slower volume recovery was observed in ethanol-grown cells as compared to glucose-grown cells. Dependence on key regulators as well as the global gene expression profile were similar in many ways to those previously observed in glucose-grown cells. However, there are indications that cells re-arrange redox-metabolism when respiration is hampered under osmostress, a feature that could not be observed in glucose-grown cells.

Kinase Screening in Pichia pastoris Identified Promising Targets Involved in Cell Growth and Alcohol Oxidase 1 Promoter (PAOX1) Regulation

As one of the most commonly used eukaryotic recombinant protein expression systems, P. pastoris relies heavily on the AOX1 promoter (PAOX1), which is strongly induced by methanol but strictly repressed by glycerol and glucose. However, the complicated signaling pathways involved in PAOX1 regulation when supplemented with different carbon sources are poorly understood. Here we constructed a kinase deletion library in P. pastoris and identified 27 mutants which showed peculiar phenotypes in cell growth or PAOX1 regulation. We analyzed both annotations and possible functions of these 27 targets, and then focused on the MAP kinase Hog1. In order to locate its potential downstream components, we performed the phosphoproteome analysis on glycerol cultured WT and Δhog1 strains and identified 157 differentially phosphorylated proteins. Our results identified important kinases involved in P. pastoris cell growth and PAOX1 regulation, which could serve as valuable targets for further mechanistic studies.

Genome-Wide Transcription Study of Cryptococcus neoformans H99 Clinical Strain versus Environmental Strains

The infection of Cryptococcus neoformans is acquired through the inhalation of desiccated yeast cells and basidiospores originated from the environment, particularly from bird’s droppings and decaying wood. Three environmental strains of C. neoformans originated from bird droppings (H4, S48B and S68B) and C. neoformans reference clinical strain (H99) were used for intranasal infection in C57BL/6 mice. We showed that the H99 strain demonstrated higher virulence compared to H4, S48B and S68B strains. To examine if gene expression contributed to the different degree of virulence among these strains, a genome-wide microarray study was performed to inspect the transcriptomic profiles of all four strains. Our results revealed that out of 7,419 genes (22,257 probes) examined, 65 genes were significantly up-or down-regulated in H99 versus H4, S48B and S68B strains. The up-regulated genes in H99 strain include Hydroxymethylglutaryl-CoA synthase (MVA1), Mitochondrial matrix factor 1 (MMF1), Bud-site-selection protein 8 (BUD8), High affinity glucose transporter 3 (SNF3) and Rho GTPase-activating protein 2 (RGA2). Pathway annotation using DAVID bioinformatics resource showed that metal ion binding and sugar transmembrane transporter activity pathways were highly expressed in the H99 strain. We suggest that the genes and pathways identified may possibly play crucial roles in the fungal pathogenesis.

Gene expression associated with intersterility in Heterobasidion

Intersterility (IS) is thought to prevent mating compatibility between homokaryons that belong to different species. Although IS in Heterobasidion is regulated by the genes located at the IS loci, it is not yet known how the IS genes influence sexual compatibility and heterokaryon formation. To increase our understanding of the molecular events underlying IS, we studied mRNA abundance changes during IS compatible and incompatible interactions over time. The clustering of the transcripts into expression profiles, followed by the application of Gene Ontology (GO) enrichment pathway analysis of each of the clusters, allowed inference of biological processes participating in IS. These analyses identified events involved in mating and sexual development (i.e., linked with IS compatibility), which included processes associated with cell-cell adhesion and recognition, cell cycle control and signal transduction. We also identified events potentially involved in overriding mating between individuals belonging to different species (i.e., linked with IS incompatibility), which included reactive oxygen species (ROS) production, responses to stress (especially to oxidative stress), signal transduction and metabolic biosynthesis. Our findings thus enabled detection and characterization of gene expression changes associated with IS in Heterobasidion, as well as identification of important processes and pathways associated with this phenomenon. Overall, the results of this study increase current knowledge regarding the molecular mechanisms underpinning IS in Heterobasidion and allowed for the establishment of a vital baseline for further studies.

Bifunctional homodimeric triokinase/FMN cyclase: contribution of protein domains to the activities of the human enzyme and molecular dynamics simulation of domain movements

Mammalian triokinase, which phosphorylates exogenous dihydroxyacetone and fructose-derived glyceraldehyde, is neither molecularly identified nor firmly associated to an encoding gene. Human FMN cyclase, which splits FAD and other ribonucleoside diphosphate-X compounds to ribonucleoside monophosphate and cyclic X-phosphodiester, is identical to a DAK-encoded dihydroxyacetone kinase. This bifunctional protein was identified as triokinase. It was modeled as a homodimer of two-domain (K and L) subunits. Active centers lie between K1 and L2 or K2 and L1: dihydroxyacetone binds K and ATP binds L in different subunits too distant (≈ 14 Å) for phosphoryl transfer. FAD docked to the ATP site with ribityl 4'-OH in a possible near-attack conformation for cyclase activity. Reciprocal inhibition between kinase and cyclase reactants confirmed substrate site locations. The differential roles of protein domains were supported by their individual expression: K was inactive, and L displayed cyclase but not kinase activity. The importance of domain mobility for the kinase activity of dimeric triokinase was highlighted by molecular dynamics simulations: ATP approached dihydroxyacetone at distances below 5 Å in near-attack conformation. Based upon structure, docking, and molecular dynamics simulations, relevant residues were mutated to alanine, and kcat and Km were assayed whenever kinase and/or cyclase activity was conserved. The results supported the roles of Thr(112) (hydrogen bonding of ATP adenine to K in the closed active center), His(221) (covalent anchoring of dihydroxyacetone to K), Asp(401) and Asp(403) (metal coordination to L), and Asp(556) (hydrogen bonding of ATP or FAD ribose to L domain). Interestingly, the His(221) point mutant acted specifically as a cyclase without kinase activity.

Genome sequence and analysis of methylotrophic yeast Hansenula polymorpha DL1

BACKGROUND

Hansenula polymorpha DL1 is a methylotrophic yeast, widely used in fundamental studies of methanol metabolism, peroxisome biogenesis and function, and also as a microbial cell factory for production of recombinant proteins and metabolic engineering towards the goal of high temperature ethanol production.

RESULTS

We have sequenced the 9 Mbp H. polymorpha DL1 genome and performed whole-genome analysis for the H. polymorpha transcriptome obtained from both methanol- and glucose-grown cells. RNA-seq analysis revealed the complex and dynamic character of the H. polymorpha transcriptome under the two studied conditions, identified abundant and highly unregulated expression of 40% of the genome in methanol grown cells, and revealed alternative splicing events. We have identified subtelomerically biased protein families in H. polymorpha, clusters of LTR elements at G + C-poor chromosomal loci in the middle of each of the seven H. polymorpha chromosomes, and established the evolutionary position of H. polymorpha DL1 within a separate yeast clade together with the methylotrophic yeast Pichia pastoris and the non-methylotrophic yeast Dekkera bruxellensis. Intergenome comparisons uncovered extensive gene order reshuffling between the three yeast genomes. Phylogenetic analyses enabled us to reveal patterns of evolution of methylotrophy in yeasts and filamentous fungi.

CONCLUSIONS

Our results open new opportunities for in-depth understanding of many aspects of H. polymorpha life cycle, physiology and metabolism as well as genome evolution in methylotrophic yeasts and may lead to novel improvements toward the application of H. polymorpha DL-1 as a microbial cell factory.

Engineering of the glycerol decomposition pathway and cofactor regulation in an industrial yeast improves ethanol production

Glycerol is a major by-product of industrial ethanol production and its formation consumes up to 4 % of the sugar substrate. This study modified the glycerol decomposition pathway of an industrial strain of Saccharomyces cerevisiae to optimize the consumption of substrate and yield of ethanol. This study is the first to couple glycerol degradation with ethanol formation, to the best of our knowledge. The recombinant strain overexpressing GCY1 and DAK1, encoding glycerol dehydrogenase and dihydroxyacetone kinase, respectively, in glycerol degradation pathway, exhibited a moderate increase in ethanol yield (2.9 %) and decrease in glycerol yield (24.9 %) compared to the wild type with the initial glucose concentration of 15 % under anaerobic conditions. However, when the mhpF gene, encoding acetylating NAD⁺-dependent acetaldehyde dehydrogenase from Escherichia coli, was co-expressed in the aforementioned recombinant strain, a further increase in ethanol yield by 5.5 % and decrease in glycerol yield by 48 % were observed for the resultant recombinant strain GDMS1 when acetic acid was added into the medium prior to inoculation compared to the wild type. The process outlined in this study which enhances glycerol consumption and cofactor regulation in an industrial yeast is a promising metabolic engineering strategy to increase ethanol production by reducing the formation of glycerol.

Nucleotide degradation and ribose salvage in yeast

Nucleotide degradation is a universal metabolic capability. Here we combine metabolomics, genetics and biochemistry to characterize the yeast pathway. Nutrient starvation, via PKA, AMPK/SNF1, and TOR, triggers autophagic breakdown of ribosomes into nucleotides. A protein not previously associated with nucleotide degradation, Phm8, converts nucleotide monophosphates into nucleosides. Downstream steps, which involve the purine nucleoside phosphorylase, Pnp1, and pyrimidine nucleoside hydrolase, Urh1, funnel ribose into the nonoxidative pentose phosphate pathway. During carbon starvation, the ribose-derived carbon accumulates as sedoheptulose-7-phosphate, whose consumption by transaldolase is impaired due to depletion of transaldolase's other substrate, glyceraldehyde-3-phosphate. Oxidative stress increases glyceraldehyde-3-phosphate, resulting in rapid consumption of sedoheptulose-7-phosphate to make NADPH for antioxidant defense. Ablation of Phm8 or double deletion of Pnp1 and Urh1 prevent effective nucleotide salvage, resulting in metabolite depletion and impaired survival of starving yeast. Thus, ribose salvage provides means of surviving nutrient starvation and oxidative stress.

Adaptive response to acetic acid in the highly resistant yeast species Zygosaccharomyces bailii revealed by quantitative proteomics

Zygosaccharomyces bailii is the most tolerant yeast species to acetic acid-induced toxicity, being able to grow in the presence of concentrations of this food preservative close to the legal limits. For this reason, Z. bailii is the most important microbial contaminant of acidic food products but the mechanisms behind this intrinsic resistance to acetic acid are very poorly characterized. To gain insights into the adaptive response and tolerance to acetic acid in Z. bailii, we explored an expression proteomics approach, based on quantitative 2DE, to identify alterations occurring in the protein content in response to sudden exposure or balanced growth in the presence of an inhibitory but nonlethal concentration of this weak acid. A coordinate increase in the content of proteins involved in cellular metabolism, in particular, in carbohydrate metabolism (Mdh1p, Aco1p, Cit1p, Idh2p, and Lpd1p) and energy generation (Atp1p and Atp2p), as well as in general and oxidative stress response (Sod2p, Dak2p, Omp2p) was registered. Results reinforce the concept that glucose and acetic acid are coconsumed in Z. bailii, with acetate being channeled into the tricarboxylic acid cycle. When acetic acid is the sole carbon source, results suggest the activation of gluconeogenic and pentose phosphate pathways, based on the increased content of several proteins of these pathways after glucose exhaustion.

Characterization of GCY1 in Saccharomyces cerevisiae by metabolic profiling

AIMS

The analytical study of intracellular (IC) metabolites has developed with advances in chromatography-linked mass spectrometry and fast sampling procedures. We applied the IC metabolite analysis to characterize the role of GCY1 in the glycerol (GLY) catabolic pathway in Saccharomyces cerevisiae.

METHODS AND RESULTS

Strains with disrupted or overexpressing GLY catabolic genes such as GCY1, DAK1 and DAK2 were constructed. The strains were cultivated under different aeration conditions and quickly quenched using a novel rapid sampling port. IC concentrations of GLY, dihydroxyacetone (DHA), glycerol 3-phosphate and dihydroxyacetone phosphate were analysed in the strains by gas chromatography/mass spectrometry. DHA was not detected in the gcy1 gene-disrupted strain but accumulated 225.91 μmol g DCW(-1) in a DHA kinase gene-deficient strain under micro-aerobic conditions. Additionally, a 16.1% increase in DHA occurred by overexpressing GCY1 in the DHA kinase-deficient strain.

CONCLUSIONS

Metabolic profiling showed that the GCY1 gene product functions as a GLY dehydrogenase in S. cerevisiae, particularly under micro-aerobic conditions.

SIGNIFICANCE AND IMPACT OF THE STUDY

Metabolic profiling of the GLY dissimilation pathway was successfully demonstrated in S. cerevisiae, and the function of GCY1 was explained by the results.

A constraint-based model analysis of the metabolic consequences of increased NADPH oxidation in Saccharomyces cerevisiae

Controlling the amounts of redox cofactors to manipulate metabolic fluxes is emerging as a useful approach to optimizing byproduct yields in yeast biotechnological processes. Redox cofactors are extensively interconnected metabolites, so predicting metabolite patterns is challenging and requires in-depth knowledge of how the metabolic network responds to a redox perturbation. Our aim was to analyze comprehensively the metabolic consequences of increased cytosolic NADPH oxidation during yeast fermentation. Using a genetic device based on the overexpression of a modified 2,3-butanediol dehydrogenase catalyzing the NADPH-dependent reduction of acetoin into 2,3-butanediol, we increased the NADPH demand to between 8 and 40-fold the anabolic demand. We developed (i) a dedicated constraint-based model of yeast fermentation and (ii) a constraint-based modeling method based on the dynamical analysis of mass distribution to quantify the in vivo contribution of pathways producing NADPH to the maintenance of redox homeostasis. We report that yeast responds to NADPH oxidation through a gradual increase in the flux through the PP and acetate pathways, providing 80% and 20% of the NADPH demand, respectively. However, for the highest NADPH demand, the model reveals a saturation of the PP pathway and predicts an exchange between NADH and NADPH in the cytosol that may be mediated by the glycerol-DHA futile cycle. We also reveal the contribution of mitochondrial shuttles, resulting in a net production of NADH in the cytosol, to fine-tune the NADH/NAD(+) balance. This systems level study helps elucidate the physiological adaptation of yeast to NADPH perturbation. Our findings emphasize the robustness of yeast to alterations in NADPH metabolism and highlight the role of the glycerol-DHA cycle as a redox valve, providing additional NADPH from NADH under conditions of very high demand.

Kinetic analysis of a Saccharomyces cerevisiae strain adapted for improved growth on glycerol: Implications for the development of yeast bioprocesses on glycerol

Glycerol is an agro-industrial residue generated in high amounts during the biodiesel production. The growing production of biodiesel is creating a worldwide glycerol surplus. Therefore, replacing sugar-based feedstock in bioprocesses by glycerol could be potentially attractive. Saccharomyces cerevisiae is one of the most commonly used microorganisms in the agri-food industry and therefore currently produced in large quantities from sugar-based feedstock. Unfortunately, growth of S. cerevisiae strains on glycerol is very low with reported μmax around 0.01 h(-1). This study demonstrates that successive growth of the S. cerevisiae CBS 8066, CEN.PK 113-7 D and Ethanol Red on glycerol as sole carbon source considerably improved the μmax from 0.01 up to 0.2 h(-1). The "adapted strain" CBS 8066-FL20 was kinetically characterized during aerobic and oxygen-limited cultivation in bioreactor and the results discussed in terms of their implication for developing glycerol-based S. cerevisiae bioprocesses.

A survey of yeast genomic assays for drug and target discovery

Over the past decade, the development and application of chemical genomic assays using the model organism Saccharomyces cerevisiae has provided powerful methods to identify the mechanism of action of known drugs and novel small molecules in vivo. These assays identify drug target candidates, genes involved in buffering drug target pathways and also help to define the general cellular response to small molecules. In this review, we examine current yeast chemical genomic assays and summarize the potential applications of each approach.

The gld1+ gene encoding glycerol dehydrogenase is required for glycerol metabolism in Schizosaccharomyces pombe

The budding yeast Saccharomyces cerevisiae is able to utilize glycerol as the sole carbon source via two pathways (glycerol 3-phosphate pathway and dihydroxyacetone [DHA] pathway). In contrast, the fission yeast Schizosaccharomyces pombe does not grow on media containing glycerol as the sole carbon source. However, in the presence of other carbon sources such as galactose and ethanol, S. pombe could assimilate glycerol and glycerol was preferentially utilized over ethanol and galactose. No equivalent of S. cerevisiae Gcy1/glycerol dehydrogenase has been identified in S. pombe. However, we identified a gene in S. pombe, SPAC13F5.03c (gld1 (+)), that is homologous to bacterial glycerol dehydrogenase. Deletion of gld1 caused a reduction in glycerol dehydrogenase activity and prevented glycerol assimilation. The gld1 Delta cells grew on 50 mM DHA as the sole carbon source, indicating that the glycerol dehydrogenase encoded by gld1 (+) is essential for glycerol assimilation in S. pombe. Strains of S. pombe deleted for dak1 (+) and dak2 (+) encoding DHA kinases could not grow on glycerol and showed sensitivity to a higher concentration of DHA. The dak1 Delta strain showed a more severe reduction of growth on glycerol and DHA than the dak2 Delta strain because the expression of dak1 (+) mRNA was higher than that of dak2 (+). In wild-type S. pombe, expression of the gld1 (+), dak1 (+), and dak2 (+) genes was repressed at a high concentration of glucose and was derepressed during glucose starvation. We found that gld1 (+) was regulated by glucose repression and that it was derepressed in scr1 Delta and tup12 Delta strains.

Computing the shortest elementary flux modes in genome-scale metabolic networks

MOTIVATION

Elementary flux modes (EFMs) represent a key concept to analyze metabolic networks from a pathway-oriented perspective. In spite of considerable work in this field, the computation of the full set of elementary flux modes in large-scale metabolic networks still constitutes a challenging issue due to its underlying combinatorial complexity.

RESULTS

In this article, we illustrate that the full set of EFMs can be enumerated in increasing order of number of reactions via integer linear programming. In this light, we present a novel procedure to efficiently determine the K-shortest EFMs in large-scale metabolic networks. Our method was applied to find the K-shortest EFMs that produce lysine in the genome-scale metabolic networks of Escherichia coli and Corynebacterium glutamicum. A detailed analysis of the biological significance of the K-shortest EFMs was conducted, finding that glucose catabolism, ammonium assimilation, lysine anabolism and cofactor balancing were correctly predicted. The work presented here represents an important step forward in the analysis and computation of EFMs for large-scale metabolic networks, where traditional methods fail for networks of even moderate size.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Reduced formation of advanced glycation endproducts via interactions between glutathione peroxidase 3 and dihydroxyacetone kinase 1

Dihydroxyacetone (DHA) induces the formation of advanced glycation endproducts (AGEs), which are involved in several diseases. Earlier, we identified dihydroxyacetone kinase 1 (Dak1) as a candidate glutathione peroxidase 3 (Gpx3)-interacting protein in Saccharomyces cerevisiae. This finding is noteworthy, as no clear evidence on the involvement of oxidative stress systems in DHA-induced AGE formation has been found to date. Here, we demonstrate that Gpx3 interacts with Dak1, alleviates DHA-mediated stress by upregulating Dak activity, and consequently suppresses AGE formation. Based on these results, we propose that defense systems against oxidative stress and DHA-induced AGE formation are related via interactions between Gpx3 and Dak1.

Proteome profiling of aging in mouse models: differential expression of proteins involved in metabolism, transport, and stress response in kidney

Aging is a time-dependent complex biological phenomenon observed in various organs and organelles of all living organisms. To understand the molecular mechanism of age-associated functional loss in aging kidneys, we have analyzed the expression of proteins in the kidneys of young (19-22 wk) and old (24 months) C57/BL6 male mice using 2-DE followed by LC-MS/MS. We found that expression levels of 49 proteins were upregulated (p < or = 0.05), while that of only ten proteins were downregulated (p < or = 0.05) due to aging. The proteins identified belong to three broad functional categories: (i) metabolism (e.g., aldehyde dehydrogenase family, ATP synthase beta-subunit, malate dehydrogenase, NADH dehydrogenase (ubiquinone), hydroxy acid oxidase 2), (ii) transport (e.g., transferrin), and (iii) chaperone/stress response (e.g., Ig-binding protein, low density lipoprotein receptor-related protein associated protein 1, selenium-binding proteins (SBPs)). Some proteins with unknown functions were also identified as being differentially expressed. ATP synthase beta subunit, transferrin, fumarate hydratase, SBPs, and albumin are present in multiple forms, possibly arising due to proteolysis or PTMs. The above functional categories suggest specific mechanisms and pathways for age-related kidney degeneration.

Dihydroxyacetone-induced death is accompanied by advanced glycation endproduct formation in selected proteins ofSaccharomyces cerevisiaeandCaenorhabditis elegans

Advanced glycation endproduct (AGE) formation is an important mechanism for protein deterioration during diabetic complications and ageing. The effects on AGE formation following dihydroxyacetone (DHA) stress were studied in two model organisms, the yeast Saccharomyces cerevisiae and the nematode Caenorhabditis elegans. Total protein AGEs, detected using an anti-N(epsilon)-carboxyalkyllysine-specific monoclonal antibody, displayed a strong correlation to DHA-induced yeast cell mortality in the wild-type and hypersensitive as well as resistant mutant strains. During DHA-induced cell death we also detected AGEs as the formation of acidic protein modifications by 2-D PAGE. Furthermore, we confirmed AGE targets immunologically on 2-D gel-separated protein extracted from DHA-treated cells. AGE modification of several metabolic enzymes (Eno2p, Adh1p, Met6 and Pgk1p) and actin (Act1p) displayed a strong correlation to DHA-induced cell death. DHA was toxic to C. elegans even at low concentration and also in this organism AGE formation accompanied death. We propose the use of DHA as a model AGE-generating substance for its apparent lack of a clear oxidative stress connection, and yeast and worm as model organisms to identify genetic determinants of protein AGE formation.