Exploring the diversity of complex metabolic networks

MOTIVATION: Metabolism, the network of chemical reactions that make life possible, is one of the most complex processes in nature. We describe here the development of a computational approach for the identification of every possible biochemical reaction from a given set of enzyme reaction rules that allows the de novo synthesis of metabolic pathways composed of these reactions, and the evaluation of these novel pathways with respect to their thermodynamic properties. RESULTS: We applied this framework to the analysis of the aromatic amino acid pathways and discovered almost 75,000 novel biochemical routes from chorismate to phenylalanine, more than 350,000 from chorismate to tyrosine, but only 13 from chorismate to tryptophan. Thermodynamic analysis of these pathways suggests that the native pathways are thermodynamically more favorable than the alternative possible pathways. The pathways generated involve compounds that exist in biological databases, as well as compounds that exist in chemical databases and novel compounds, suggesting novel biochemical routes for these compounds and the existence of biochemical compounds that remain to be discovered or synthesized through enzyme and pathway engineering. AVAILABILITY: Framework will be available via web interface at http://systemsbiology.northwestern.edu/BNICE (site under construction). CONTACT: vassily@northwestern.edu or broadbelt@northwestern.edu SUPPLEMENTARY INFORMATION: http://systemsbiology.northwestern.edu/BNICE/publications.     

Computational framework for predictive biodegradation

As increasing amounts of anthropogenic chemicals are released into the environment, it is vital to human health and the preservation of ecosystems to evaluate the fate of these chemicals in the environment. It is useful to predict whether a particular compound is biodegradable and if alternate routes can be engineered for compounds already known to be biodegradable. In this work, we describe a computational framework (called BNICE) that can be used for the prediction of novel biodegradation pathways of xenobiotics. The framework was applied to 4‐chlorobiphenyl, phenanthrene, γ‐hexachlorocyclohexane, and 1,2,4‐trichlorobenzene, compounds representing various classes of xenobiotics with known biodegradation routes. BNICE reproduced the proposed biodegradation routes found experimentally, and in addition, it expanded the biodegradation reaction networks through the generation of novel compounds and reactions. The novel reactions involved in the biodegradation of 1,2,4‐trichlorobenzene were studied in depth, where pathway and thermodynamic analyses were performed. This work demonstrates that BNICE can be applied to generate novel pathways to degrade xenobiotic compounds that are thermodynamically feasible alternatives to known biodegradation routes and attractive targets for metabolic engineering. Biotechnol. Bioeng. 2009; 104: 1086–1097. © 2009 Wiley Periodicals, Inc.     

Exploring De Novo metabolic pathways from pyruvate to propionic acid

Industrial biotechnology provides an efficient, sustainable solution for chemical production. However, designing biochemical pathways based solely on known reactions does not exploit its full potential. Enzymes are known to accept non‐native substrates, which may allow novel, advantageous reactions. We have previously developed a computational program named Biological Network Integrated Computational Explorer (BNICE) to predict promiscuous enzyme activities and design synthetic pathways, using generalized reaction rules curated from biochemical reaction databases. Here, we use BNICE to design pathways synthesizing propionic acid from pyruvate. The currently known natural pathways produce undesirable by‐products lactic acid and succinic acid, reducing their economic viability. BNICE predicted seven pathways containing four reaction steps or less, five of which avoid these by‐products. Among the 16 biochemical reactions comprising these pathways, 44% were validated by literature references. More than 28% of these known reactions were not in the BNICE training dataset, showing that BNICE was able to predict novel enzyme substrates. Most of the pathways included the intermediate acrylic acid. As acrylic acid bioproduction has been well advanced, we focused on the critical step of reducing acrylic acid to propionic acid. We experimentally validated that Oye2p from Saccharomyces cerevisiae can catalyze this reaction at a slow turnover rate (10^?3 s^?1), which was unknown to occur with this enzyme, and is an important finding for further propionic acid metabolic engineering. These results validate BNICE as a pathway‐searching tool that can predict previously unknown promiscuous enzyme activities and show that computational methods can elucidate novel biochemical pathways for industrial applications. © 2016 American Institute of Chemical Engineers Biotechnol. Prog., 32:303–311, 2016     

Integrating computational methods to retrofit enzymes to synthetic pathways

Microbial production of desired compounds provides an efficient framework for the development of renewable energy resources. To be competitive to traditional chemistry, one requirement is to utilize the full capacity of the microorganism to produce target compounds with high yields and turnover rates. We use integrated computational methods to generate and quantify the performance of novel biosynthetic routes that contain highly optimized catalysts. Engineering a novel reaction pathway entails addressing feasibility on multiple levels, which involves handling the complexity of large-scale biochemical networks while respecting the critical chemical phenomena at the atomistic scale. To pursue this multi-layer challenge, our strategy merges knowledge-based metabolic engineering methods with computational chemistry methods. By bridging multiple disciplines, we provide an integral computational framework that could accelerate the discovery and implementation of novel biosynthetic production routes. Using this approach, we have identified and optimized a novel biosynthetic route for the production of 3HP from pyruvate.     

Exploring biochemical pathways for mono-ethylene glycol (MEG) synthesis from synthesis gas

Mono-ethylene glycol (MEG) is an important petrochemical with widespread use in numerous consumer products. The current industrial MEG-production process relies on non-renewable fossil fuel-based feedstocks, such as petroleum, natural gas, and naphtha; hence, it is useful to explore alternative routes of MEG-synthesis from gases as they might provide a greener and more sustainable alternative to the current production methods. Technologies of synthetic biology and metabolic engineering of microorganisms can be deployed for the expression of new biochemical pathways for MEG-synthesis from gases, provided that such promising alternative routes are first identified. We used the BNICE.ch algorithm to develop novel and previously unknown biological pathways to MEG from synthesis gas by leveraging the Wood-Ljungdahl pathway of carbon fixation of acetogenic bacteria. We developed a set of useful pathway pruning and analysis criteria to systematically assess thousands of pathways generated by BNICE.ch. Published genome-scale models of Moorella thermoacetica and Clostridium ljungdahlii were used to perform the pathway yield calculations and in-depth analyses of seven (7) newly developed biological MEG-producing pathways from gases, including CO₂, CO, and H₂. These analyses helped identify not only better candidate pathways, but also superior chassis organisms that can be used for metabolic engineering of the candidate pathways. The pathway generation, pruning, and detailed analysis procedures described in this study can also be used to develop biochemical pathways for other commodity chemicals from gaseous substrates.     

Recent advances in the metabolic engineering of microorganisms for the production of 3-hydroxypropionic acid as C3 platform chemical

Development of sustainable technologies for the production of 3-hydroxypropionic acid (3HP) as a platform chemical has recently been gaining much attention owing to its versatility in applications for the synthesis of other specialty chemicals. Several proposed biological synthesis routes and strategies for producing 3HP from glucose and glycerol are reviewed presently. Ten proposed routes for 3HP production from glucose are described and one of which was recently constructed successfully in Escherichia coli with malonyl–Coenzyme A as a precursor. This resulted in a yield still far from the required level for industrial application. On the other hand, strategies employing engineered E. coli and Klebsiella pneumoniae capable of producing 3HP from glycerol are also evaluated. The titers produced by these recombinant strains reached around 3 %. At its current state, it is evident that a bulk of engineering works is yet to be done to acquire a biosynthesis route for 3HP that is acceptable for industrial-scale production.     

Combination of Entner-Doudoroff Pathway with MEP Increases Isoprene Production in Engineered Escherichia coli

Embden-Meyerhof pathway (EMP) in tandem with 2-C-methyl-D-erythritol 4-phosphate pathway (MEP) is commonly used for isoprenoid biosynthesis in E. coli. However, this combination has limitations as EMP generates an imbalanced distribution of pyruvate and glyceraldehyde-3-phosphate (G3P). Herein, four glycolytic pathways—EMP, Entner-Doudoroff Pathway (EDP), Pentose Phosphate Pathway (PPP) and Dahms pathway were tested as MEP feeding modules for isoprene production. Results revealed the highest isoprene production from EDP containing modules, wherein pyruvate and G3P were generated simultaneously; isoprene titer and yield were more than three and six times higher than those of the EMP module, respectively. Additionally, the PPP module that generates G3P prior to pyruvate was significantly more effective than the Dahms pathway, in which pyruvate production precedes G3P. In terms of precursor generation and energy/reducing-equivalent supply, EDP+PPP was found to be the ideal feeding module for MEP. These findings may launch a new direction for the optimization of MEP-dependent isoprenoid biosynthesis pathways.     

Biosynthesis of isoprenoids in higher plant chloroplasts proceeds via a mevalonate-independent pathway

Isopentenyl diphosphate (IPP) is the biological C₅ precursor of isoprenoids. By labeling experiments using [1-¹³C]glucose, higher plants were shown to possess two distinct biosynthetic routes for IPP biosynthesis: while the cytoplasmic sterols were formed via the acetate/mevalonate pathway, the chloroplast-bound isoprenoids (β-carotene, lutein, prenyl chains of chlorophylls and plastoquinone-9) were synthesized via a novel IPP biosynthesis pathway (glyceraldehyde phosphate/pyruvate pathway) which was first found in eubacteria and a green alga. The dichotomy in isoprenoid biosynthesis in higher plants allows a reasonable interpretation of previous odd and inconclusive results concerning the biosynthesis of chloroplast isoprenoids, which so far had mainly been interpreted in the frame of models using compartmentation of the mevalonate pathway.     

Auxin-mediated induction of GAL promoters by conditional degradation of Mig1p improves sesquiterpene production in Saccharomyces cerevisiae with engineered acetyl-CoA synthesis

The yeast Saccharomyces cerevisiae uses the pyruvate dehydrogenase-bypass for acetyl-CoA biosynthesis. This relatively inefficient pathway limits production potential for acetyl-CoA-derived biochemical due to carbon loss and the cost of two high-energy phosphate bonds per molecule of acetyl-CoA. Here, we attempted to improve acetyl-CoA production efficiency by introducing heterologous acetylating aldehyde dehydrogenase and phosphoketolase pathways for acetyl-CoA synthesis to enhance production of the sesquiterpene trans-nerolidol. In addition, we introduced auxin-mediated degradation of the glucose-dependent repressor Mig1p to allow induced expression of GAL promoters on glucose so that production potential on glucose could be examined. The novel genes that we used to reconstruct the heterologous acetyl-CoA pathways did not sufficiently complement the loss of endogenous acetyl-CoA pathways, indicating that superior heterologous enzymes are necessary to establish fully functional synthetic acetyl-CoA pathways and properly explore their potential for nerolidol synthesis. Notwithstanding this, nerolidol production was improved twofold to a titre of ˜ 900 mg l⁻¹ in flask cultivation using a combination of heterologous acetyl-CoA pathways and Mig1p degradation. Conditional Mig1p depletion is presented as a valuable strategy to improve the productivities in the strains engineered with GAL promoters-controlled pathways when growing on glucose.     

Understanding the central role of citrate in the metabolism of cancer cells

Cancers cells strongly stimulate glycolysis and glutaminolysis for their biosynthesis. Pyruvate derived from glucose is preferentially diverted towards the production of lactic acid (Warburg effect). Citrate censors ATP production and controls strategic enzymes of anabolic and catabolic pathways through feedback reactions. Mitochondrial citrate diffuses in the cytosol to restore oxaloacetate and acetyl-CoA. Whereas acetyl-CoA serves de novo lipid synthesis and histone acetylation, OAA is derived towards lactate production via pyruvate and / or a vicious cycle reforming mitochondrial citrate. This cycle allows cancer cells to burn their host's lipid and protein reserves in order to sustain their own biosynthesis pathways. In vitro, citrate has demonstrated anti-cancer properties when administered in excess, sensitizing cancer cells to chemotherapy. Understanding its central role is of particular relevance for the development of new strategies for counteracting cancer cell proliferation and overcoming chemoresistance.     

Computer-aided synthesis of biochemical pathways

The synthesis of biochemical pathways satisfying stoichiometric constraints is discussed. Stoichiometric constraints arise primarily from designating compounds as required or allowed reactants, and required or allowed products of the pathways; they also arise from similar restrictions on intermediate metabolites and bioreactions participating in the pathways. An algorithm for the complete and correct solution of the problem is presented; the algorithm satisfies each constraint by recursively transforming a base-set of pathways. The algorithm is applied to the problem of lysine synthesis from glucose and ammonia. In addition to the established synthesis routes, the algorithm constructs several alternative pathways that bypass key enzymes, such as malate dehydrogenase and pyruvate dehydrogenase. Apart from the construction of pathways with desired characteristics, the systematic synthesis of pathways can also uncover fundamental constraints in a particular problem, by demonstrating that no pathways exist to meet certain sets of specifications. In the case of lysine, the algorithm shows that oxaloacetate is a necessary intermediate in all pathways leading to lysine from glucose, and that the yield of lysine over glucose cannot exceed 67% in the absence of enzymatic recovery of carbon dioxide.     

In silico assessment of the metabolic capabilities of an engineered functional reversal of the β-oxidation cycle for the synthesis of longer-chain (C≥4) products

The modularity and versatility of an engineered functional reversal of the β-oxidation cycle make it a promising platform for the synthesis of longer-chain (C≥4) products. While the pathway has recently been exploited for the production of n-alcohols and carboxylic acids, fully capitalizing on its potential for the synthesis of a diverse set of product families requires a system-level assessment of its biosynthetic capabilities. To this end, we utilized a genome scale model of Escherichia coli, in combination with Flux Balance Analysis and Flux Variability Analysis, to determine the key characteristics and constraints of this pathway for the production of a variety of product families under fermentative conditions. This analysis revealed that the production of n-alcohols, alkanes, and fatty acids of lengths C3–C18 could be coupled to cell growth in a strain lacking native fermentative pathways, a characteristic enabling product synthesis at maximum rates, titers, and yields. While energetic and redox constraints limit the production of target compounds from alternative platforms such as the fatty acid biosynthesis and α-ketoacid pathways, the metabolic efficiency of a β-oxidation reversal allows the production of a wide range of products of varying length and functionality. The versatility of this platform was investigated through the simulation of various termination pathways for product synthesis along with the use of different priming molecules, demonstrating its potential for the efficient synthesis of a wide variety of functionalized compounds. Overall, specific metabolic manipulations suggested by this systems-level analysis include deletion of native fermentation pathways, the choice of priming molecules and specific routes for their synthesis, proper choice of termination enzymes, control of flux partitioning at the pyruvate node and the pentose phosphate pathway, and the use of an NADH-dependent trans-enoyl-CoA reductase instead of a ferredoxin-dependent enzyme.     

A computational framework for integration of lipidomics data into metabolic pathways

Lipids are important compounds for human physiology and as renewable resources for fuels and chemicals. In lipid research, there is a big gap between the currently available pathway-level representations of lipids and lipid structure databases in which the number of compounds is expanding rapidly with high-throughput mass spectrometry methods.In this work, we introduce a computational approach to bridge this gap by making associations between metabolic pathways and the lipid structures discovered increasingly thorough lipidomics studies. Our approach, called NICELips (Network Integrated Computational Explorer for Lipidomics), is based on the formulation of generalized enzymatic reaction rules for lipid metabolism, and it employs the generalized rules to postulate novel pathways of lipid metabolism. It further integrates all discovered lipids in biological networks of enzymatic reactions that consist their biosynthesis and biodegradation pathways.We illustrate the utility of our approach through a case study of bis(monoacylglycero)phosphate (BMP), a biologically important glycerophospholipid with immature synthesis and catabolic route(s). Using NICELips, we were able to propose various synthesis and degradation pathways for this compound and several other lipids with unknown metabolism like BMP, and in addition several alternative novel biosynthesis and biodegradation pathways for lipids with known metabolism. NICELips has potential applications in designing therapeutic interventions for lipid-associated disorders and in the metabolic engineering of model organisms for improving the biobased production of lipid-derived fuels and chemicals.     

Genome-enabled determination of amino acid biosynthesis in Xanthomonas campestris pv. campestris and identification of biosynthetic pathways for alanine, glycine, and isoleucine by 13C-isotopologue profiling

To elucidate the biosynthetic pathways for all proteinogenic amino acids in Xanthomonas campestris pv. campestris, this study combines results obtained by in silico genome analysis and by (13)C-NMR-based isotopologue profiling to provide a panoramic view on a substantial section of bacterial metabolism. Initially, biosynthesis pathways were reconstructed from an improved annotation of the complete genome of X. campestris pv. campestris B100. This metabolic reconstruction resulted in the unequivocal identification of biosynthesis routes for 17 amino acids in total: arginine, asparagine, aspartate, cysteine, glutamate, glutamine, histidine, leucine, lysine, methionine, phenylalanine, proline, serine, threonine, tryptophan, tyrosine, and valine. Ambiguous pathways were reconstructed from the genome data for alanine, glycine, and isoleucine biosynthesis. (13)C-NMR analyses supported the identification of the metabolically active pathways. The biosynthetic routes for these amino acids were derived from the precursor molecules pyruvate, serine, and pyruvate, respectively. By combining genome analysis and isotopologue profiling, a comprehensive set of biosynthetic pathways covering all proteinogenic amino acids was unraveled for this plant pathogenic bacterium, which plays an important role in biotechnology as a producer of the exopolysaccharide xanthan. The data obtained lay ground for subsequent functional analyses in post-genomics and biotechnology, while the innovative combination of in silico and wet lab technology described here is promising as a general approach to elucidate metabolic pathways.     

Carbohydrate and amino acid metabolism in Tuber borchii mycelium during glucose utilization: a (13)C NMR study

The metabolism of [1-13C]glucose in the vegetative mycelium of the ectomycorrhizal ascomycete Tuber borchii was studied in order to characterize the biochemical pathways for the assimilation of glucose and amino acid biosynthesis. The pathways were characterized using nuclear magnetic resonance spectroscopy in conjunction with [1-13C]glucose labeling. The enzymes of mannitol cycle and ammonium assimilation were also evaluated. The majority of the 13C label was incorporated into mannitol and this polyol was formed via a direct route from absorbed glucose. Amino acid biosynthesis was also an important sink of assimilated carbon and 13C was mainly incorporated into alanine and glutamate. From this intramolecular 13C enrichment, it is concluded that pyruvate, arising from [1-13C]glucose catabolism, was used by alanine aminotransferase, pyruvate dehydrogenase and pyruvate carboxylase before entering the Krebs cycle. The transfer of 13C-labeled mycelium on [12C]glucose showed that mannitol, alanine, and glutamate carbon were used to synthesize glutamine and arginine that likely play a storage role.     

Revisiting the Central Metabolism of the Bloodstream Forms of Trypanosoma brucei: Production of Acetate in the Mitochondrion Is Essential for Parasite Viability

Background

The bloodstream forms of Trypanosoma brucei, the causative agent of sleeping sickness, rely solely on glycolysis for ATP production. It is generally accepted that pyruvate is the major end-product excreted from glucose metabolism by the proliferative long-slender bloodstream forms of the parasite, with virtually no production of succinate and acetate, the main end-products excreted from glycolysis by all the other trypanosomatid adaptative forms, including the procyclic insect form of T. brucei.

Methodology/Principal Findings

A comparative NMR analysis showed that the bloodstream long-slender and procyclic trypanosomes excreted equivalent amounts of acetate and succinate from glucose metabolism. Key enzymes of acetate production from glucose-derived pyruvate and threonine are expressed in the mitochondrion of the long-slender forms, which produces 1.4-times more acetate from glucose than from threonine in the presence of an equal amount of both carbon sources. By using a combination of reverse genetics and NMR analyses, we showed that mitochondrial production of acetate is essential for the long-slender forms, since blocking of acetate biosynthesis from both carbon sources induces cell death. This was confirmed in the absence of threonine by the lethal phenotype of RNAi-mediated depletion of the pyruvate dehydrogenase, which is involved in glucose-derived acetate production. In addition, we showed that de novo fatty acid biosynthesis from acetate is essential for this parasite, as demonstrated by a lethal phenotype and metabolic analyses of RNAi-mediated depletion of acetyl-CoA synthetase, catalyzing the first cytosolic step of this pathway.

Conclusions/Significance

Acetate produced in the mitochondrion from glucose and threonine is synthetically essential for the long-slender mammalian forms of T. brucei to feed the essential fatty acid biosynthesis through the “acetate shuttle” that was recently described in the procyclic insect form of the parasite. Consequently, key enzymatic steps of this pathway, particularly acetyl-CoA synthetase, constitute new attractive drug targets against trypanosomiasis.     

Understanding the central role of citrate in the metabolism of cancer cells

Cancers cells strongly stimulate glycolysis and glutaminolysis for their biosynthesis. Pyruvate derived from glucose is preferentially diverted towards the production of lactic acid (Warburg effect). Citrate censors ATP production and controls strategic enzymes of anabolic and catabolic pathways through feedback reactions. Mitochondrial citrate diffuses in the cytosol to restore oxaloacetate and acetyl-CoA. Whereas acetyl-CoA serves de novo lipid synthesis and histone acetylation, OAA is derived towards lactate production via pyruvate and / or a vicious cycle reforming mitochondrial citrate. This cycle allows cancer cells to burn their host's lipid and protein reserves in order to sustain their own biosynthesis pathways. In vitro, citrate has demonstrated anti-cancer properties when administered in excess, sensitizing cancer cells to chemotherapy. Understanding its central role is of particular relevance for the development of new strategies for counteracting cancer cell proliferation and overcoming chemoresistance.     

The production of organic acids.

The production of organic acids covers two aspects: first, the metabolic pathways involved in the biosynthesis, and, second, the industrial process strategy adopted. The review seeks to show the underlying biochemical similarities in the biosynthesis of organic acids and the resulting similarities in the commercial processes. Two groups of acids are defined, those with a "long" biosynthetic path from glucose, involving much of the glycolytic pathway and the tricarboxylic acid cycle, and those acids with a "short pathway", essentially a biotransformation of glucose. The regulation of the pathways and the future developments in metabolic control theory and genetic manipulations relating to them are considered. The organisms used industrially are also limited, Aspergillus sp. and Candida yeasts; again the underlying metabolic similarities lead to similar strategies for all the acids discussed.     

Generation of an atlas for commodity chemical production in Escherichia coli and a novel pathway prediction algorithm,GEM-Path

The production of 75% of the current drug molecules and 35% of all chemicals could be achieved through bioprocessing (Arundel and Sawaya, 2009). To accelerate the transition from a petroleum-based chemical industry to a sustainable bio-based industry, systems metabolic engineering has emerged to computationally design metabolic pathways for chemical production. Although algorithms able to provide specific metabolic interventions and heterologous production pathways are available, a systematic analysis for all possible production routes to commodity chemicals in Escherichia coli is lacking. Furthermore, a pathway prediction algorithm that combines direct integration of genome-scale models at each step of the search to reduce the search space does not exist. Previous work (Feist et al., 2010) performed a model-driven evaluation of the growth-coupled production potential for E. coli to produce multiple native compounds from different feedstocks. In this study, we extended this analysis for non-native compounds by using an integrated approach through heterologous pathway integration and growth-coupled metabolite production design. In addition to integration with genome-scale model integration, the GEM-Path algorithm developed in this work also contains a novel approach to address reaction promiscuity. In total, 245 unique synthetic pathways for 20 large volume compounds were predicted. Host metabolism with these synthetic pathways was then analyzed for feasible growth-coupled production and designs could be identified for 1271 of the 6615 conditions evaluated. This study characterizes the potential for E. coli to produce commodity chemicals, and outlines a generic strain design workflow to design production strains.