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     

Discovery and analysis of novel metabolic pathways for the biosynthesis of industrial chemicals: 3‐hydroxypropanoate

Sustainable microbial production of high‐value organic compounds such as 3‐hydroxypropanoate (3HP) is becoming an increasingly attractive alternative to organic syntheses that utilize petrochemical feedstocks. We applied the Biochemical Network Integrated Computational Explorer (BNICE) framework to the automated design and evaluation of novel biosynthetic routes for the production of 3HP from pyruvate. Among the pathways generated by the BNICE framework were all of the known pathways for the production of 3HP as well as numerous novel pathways. The pathways generated by BNICE were ranked based on four criteria: pathway length, thermodynamic feasibility, maximum achievable yield to 3HP from glucose, and maximum achievable activity at which 3HP can be produced. Four pathways emerged from this ranking as the most promising for the biosynthesis of 3HP, and three of these pathways, including the shortest pathways discovered, were novel. We also discovered novel routes for the biosynthesis of 28 commercially available compounds that are currently produced exclusively through organic synthesis. Examination of the optimal pathways for the biosynthesis of these 28 compounds in E. coli revealed pyruvate and succinate to be ideal intermediates for achieving high product yields from glucose. Biotechnol. Bioeng. 2010; 106: 462–473. © 2010 Wiley Periodicals, Inc.     

PathMiner: predicting metabolic pathways by heuristic search

MOTIVATION: Automated methods for biochemical pathway inference are becoming increasingly important for understanding biological processes in living and synthetic systems. With the availability of data on complete genomes and increasing information about enzyme-catalyzed biochemistry it is becoming feasible to approach this problem computationally. In this paper we present PathMiner, a system for automatic metabolic pathway inference. PathMiner predicts metabolic routes by reasoning over transformations using chemical and biological information. RESULTS: We build a biochemical state-space using data from known enzyme-catalyzed transformations in Ligand, including, 2917 unique transformations between 3890 different compounds. To predict metabolic pathways we explore this state-space by developing an informed search algorithm. For this purpose we develop a chemically motivated heuristic to guide the search. Since the algorithm does not depend on predefined pathways, it can efficiently identify plausible routes using known biochemical transformations.     

In silico metabolic engineering of Clostridium ljungdahlii for synthesis gas fermentation

Synthesis gas fermentation is one of the most promising routes to convert synthesis gas (syngas; mainly comprised of H₂ and CO) to renewable liquid fuels and chemicals by specialized bacteria. The most commonly studied syngas fermenting bacterium is Clostridium ljungdahlii, which produces acetate and ethanol as its primary metabolic byproducts. Engineering of C. ljungdahlii metabolism to overproduce ethanol, enhance the synthesize of the native byproducts lactate and 2,3-butanediol, and introduce the synthesis of non-native products such as butanol and butyrate has substantial commercial value. We performed in silico metabolic engineering studies using a genome-scale reconstruction of C. ljungdahlii metabolism and the OptKnock computational framework to identify gene knockouts that were predicted to enhance the synthesis of these native products and non-native products, introduced through insertion of the necessary heterologous pathways. The OptKnock derived strategies were often difficult to assess because increase product synthesis was invariably accompanied by decreased growth. Therefore, the OptKnock strategies were further evaluated using a spatiotemporal metabolic model of a syngas bubble column reactor, a popular technology for large-scale gas fermentation. Unlike flux balance analysis, the bubble column model accounted for the complex tradeoffs between increased product synthesis and reduced growth rates of engineered mutants within the spatially varying column environment. The two-stage methodology for deriving and evaluating metabolic engineering strategies was shown to yield new C. ljungdahlii gene targets that offer the potential for increased product synthesis under realistic syngas fermentation conditions.     

The metabolic potential of plastics as biotechnological carbon sources – Review and targets for the future

The plastic crisis requires drastic measures, especially for the plastics’ end-of-life. Mixed plastic fractions are currently difficult to recycle, but microbial metabolism might open new pathways. With new technologies for degradation of plastics to oligo- and monomers, these carbon sources can be used in biotechnology for the upcycling of plastic waste to valuable products, such as bioplastics and biosurfactants. We briefly summarize well-known monomer degradation pathways and computed their theoretical yields for industrially interesting products. With this information in hand, we calculated replacement scenarios of existing fossil-based synthesis routes for the same products. Thereby, we highlight fossil-based products for which plastic monomers might be attractive alternative carbon sources. Notably, not the highest yield of product on substrate of the biochemical route, but rather the (in-)efficiency of the petrochemical routes (i.e., carbon, energy use) determines the potential of biochemical plastic upcycling. Our results might serve as a guide for future metabolic engineering efforts towards a sustainable plastic economy.     

一碳气体利用微生物及其基因工程改造

一碳气体主要包括CO、CO_(2)和CH_(4)等,这些气体来源于陆地生物活动、工业废气以及气化合成气等,其中CO_(2)与CH_(4)是温室气体,对全球气候变化有着重要的影响。利用微生物进行一碳气体生物转化既可以解决废气排放的问题,又能生产燃料及多种化学品。近年来,运用CRISPR/Cas9等基因编辑技术对一碳气体利用微生物进行改造,是提高它们的产物得率、增加产物类型的重要途径。本文主要围绕甲烷营养菌、自养乙酸菌、一氧化碳营养菌等一碳气体利用微生物,综述了其生物学特性、好氧和厌氧代谢途径、代谢产物,以及常用的基因编辑技术(利用同源重组的基因中断技术、二类内含子ClosTron法、CRISPR/Cas基因编辑及以噬菌体重组酶介导的DNA大片段引入等)在它们中的应用,为后续相关研究提供参考。     

Production of 1-octanol in Escherichia coli by a high flux thioesterase route

1-octanol is a valuable molecule in the chemical industry, where it is used as a plasticizer, as a precursor in the production of linear low-density polyethylene (LLDPE), and as a growth inhibitor of tobacco plant suckers. Due to the low availability of eight-carbon acyl chains in natural lipid feedstocks and the selectivity challenges in petrochemical routes to medium-chain fatty alcohols,1-octanol sells for the highest price among the fatty alcohol products. As an alternative, metabolic engineers have pursued sustainable 1-octanol production via engineered microbes. Here, we report demonstration of gram per liter titers in the model bacterium Escherichia coli via the development of a pathway composed of a thioesterase, an acyl-CoA synthetase, and an acyl-CoA reductase. In addition, the impact of deleting fermentative pathways was explored E. coli K12 MG1655 strain for production of octanoic acid, a key octanol precursor. In order to overcome metabolic flux barriers, bioprospecting experiments were performed to identify acyl-CoA synthetases with high activity towards octanoic acid and acyl-CoA reductases with high activity to produce 1-octanol from octanoyl-CoA. Titration of expression of key pathway enzymes was performed and a strain with the full pathway integrated on the chromosome was created. The final strain produced 1-octanol at 1.3 g/L titer and a >90% C₈ specificity from glycerol. In addition to the metabolic engineering efforts, this work addressed some of the technical challenges that arise when quantifying 1-octanol produced from cultures grown under fully aerobic conditions where evaporation and stripping are prevalent.     

Production of bulk chemicals via novel metabolic pathways in microorganisms

Metabolic engineering has been playing important roles in developing high performance microorganisms capable of producing various chemicals and materials from renewable biomass in a sustainable manner. Synthetic and systems biology are also contributing significantly to the creation of novel pathways and the whole cell-wide optimization of metabolic performance, respectively. In order to expand the spectrum of chemicals that can be produced biotechnologically, it is necessary to broaden the metabolic capacities of microorganisms. Expanding the metabolic pathways for biosynthesizing the target chemicals requires not only the enumeration of a series of known enzymes, but also the identification of biochemical gaps whose corresponding enzymes might not actually exist in nature; this issue is the focus of this paper. First, pathway prediction tools, effectively combining reactions that lead to the production of a target chemical, are analyzed in terms of logics representing chemical information, and designing and ranking the proposed metabolic pathways. Then, several approaches for potentially filling in the gaps of the novel metabolic pathway are suggested along with relevant examples, including the use of promiscuous enzymes that flexibly utilize different substrates, design of novel enzymes for non-natural reactions, and exploration of hypothetical proteins. Finally, strain optimization by systems metabolic engineering in the context of novel metabolic pathways constructed is briefly described. It is hoped that this review paper will provide logical ways of efficiently utilizing ‘big’ biological data to design and develop novel metabolic pathways for the production of various bulk chemicals that are currently produced from fossil resources.     

Energy dependant plant stress acclimation

Plants may live and grow under suboptimal environmental conditions having certain biochemical and metabolic adaptations that facilitate their survival. Plant “metabolic flexibility” consists of the accomplishment of the same step in a metabolic pathway in a variety of different ways. Pyrophosphate which works as an energy donor when cellular ATP pools become diminished during stresses, alternative glycolytic reactions that bypass ATP-requiring steps, additional pathways for electron transport in plant mithocondria and the salvage pathways are some of the aspects related to “energetic flexibility”. This key feature that plays an important role in plant acclimation to stress can be an important target for engineering enhanced stress tolerance in crop plants.     

Metabolic engineering of Escherichia coli for high production of 1,5-pentanediol via a cadaverine-derived pathway

1,5-Pentanediol (1,5-PDO) is a high value-added chemical which is widely used as a monomer in the polymer industry. There are no natural organisms that could directly produce 1,5-PDO from renewable carbon sources. In this study, we report metabolic engineering of Escherichia coli for high-level production of 1,5-PDO from glucose via a cadaverine-derived pathway. In the newly proposed pathway, cadaverine can be converted to 1,5-PDO via 5-hydroxyvalerate (5-HV) by introducing only one heterologous enzyme in E. coli. Different endogenous genes of E. coli were screened and heterologous carboxylic acid reductase genes were tested to build a functional pathway. Compared to the previously reported pathways, the engineered cadaverine-based pathway has a higher theoretical yield (0.70 mol/mol glucose) and higher catalytic efficiency. By further combining strategies of pathway engineering and process engineering, we constructed an engineered E. coli strain that could produce 2.62 g/L 1,5-PDO in shake-flask and 9.25 g/L 1,5-PDO with a yield of 0.28 mol/mol glucose in fed-batch fermentation. The proposed new pathway and engineering strategies reported here should be useful for developing biological routes to produce 1,5-PDO for real application.     

Alternative pathways of galactose assimilation: could inverse metabolic engineering provide an alternative to galactosemic patients?

The galactose assimilation pathway has been extensively studied as an example of a genetic regulatory switch. Besides the importance of this pathway as a tool in basic biological research, unraveling its structure and regulation is also of major medical importance. Impairment of galactose assimilation is the cause of the genetic metabolic disease known as "galactosemia", while the in vivo activity of the pathway affects the production of glycans. The latter have been connected to tumor metastasis, anti-cancer drug resistance and various cardiovascular diseases. Despite the vast amount of studies, however, galactose assimilation and its interaction with other parts of the metabolic network have not been fully elucidated yet. In yeast and higher eukaryotes, it is still being studied as comprising only the linear Leloir pathway. Recent observations, however, indicate that alternative pathways of galactose assimilation identified in prokaryotes and fungi might also be present in yeast. Such a result is valuable per se, because it could lead to the discovery of these pathways in humans. Even more importantly, these pathways provide alternative phenotypes with known genetic fingerprints that can be used in the context of classical and inverse metabolic engineering to examine and treat the mechanisms of defects of galactose assimilation.     

Functional replacement of isoprenoid pathways in Rhodobacter sphaeroides

Advances in synthetic biology and metabolic engineering have proven the potential of introducing metabolic by-passes within cell factories. These pathways can provide a more efficient alternative to endogenous counterparts due to their insensitivity to host's regulatory mechanisms. In this work, we replaced the endogenous essential 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway for isoprenoid biosynthesis in the industrially relevant bacterium Rhodobacter sphaeroides by an orthogonal metabolic route. The native 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway was successfully replaced by a heterologous mevalonate (MVA) pathway from a related bacterium. The functional replacement was confirmed by analysis of the reporter molecule amorpha-4,11-diene after cultivation with [4-¹³C]glucose. The engineered R. sphaeroides strain relying exclusively on the MVA pathway was completely functional in conditions for sesquiterpene production and, upon increased expression of the MVA enzymes, it reached even higher sesquiterpene yields than the control strain coexpressing both MEP and MVA modules. This work represents an example where substitution of an essential biochemical pathway by an alternative, heterologous pathway leads to enhanced biosynthetic performance.     

The population genetics of alleles affecting enzyme activity.

It is possible to predict the population genetics of allozymes by assuming that fitness is proportional to flux through a biochemical pathway. The model presented here extends previous work by incorporating two additional features of biological realism. Firstly, that more than one biochemical route may exist between any two metabolites. The major routes have been identified as the classical biochemical pathways but in the event of a mutation blocking a major route, minor routes become significant. These minor routes are named "bypass fluxes" and have profound effects on the population genetics of allozymes. Secondly, recent work has suggested that a metabolic cost is associated with enzyme synthesis; this will constitute an additional selective pressure on alleles which affect the amount of enzyme synthesized. The model generates a fitness curve which predicts the fitness associated with any level of enzyme activity. It can utilize data on null or near-null, structural or regulatory, mutations in the presence or absence of bypass fluxes. When data from natural populations of Drosophila are investigated, it is concluded that selection pressures acting on enzyme variants may be much higher than previously thought.     

PathAligner

MOTIVATION: Analysis of metabolic pathways is a central topic in understanding the relationship between genotype and phenotype. The rapid accumulation of biological data provides the possibility of studying metabolic pathways at both the genomic and the metabolic levels. Retrieving metabolic pathways from current biological data sources, reconstructing metabolic pathways from rudimentary pathway components, and aligning metabolic pathways with each other are major tasks. Our motivation was to develop a conceptual framework and computational system that allows the retrieval of metabolic pathway information and the processing of alignments to reveal the similarities between metabolic pathways. RESULTS: PathAligner extracts metabolic information from biological databases via the Internet and builds metabolic pathways with data sources of genes, sequences, enzymes, metabolites etc. It provides an easy-to-use interface to retrieve, display and manipulate metabolic information. PathAligner also provides an alignment method to compare the similarity between metabolic pathways. AVAILABILITY: PathAligner is available at http://bibiserv.techfak.uni-bielefeld.de/pathaligner.     

Rethinking glycolysis: on the biochemical logic of metabolic pathways

Metabolic pathways may seem arbitrary and unnecessarily complex. In many cases, a chemist might devise a simpler route for the biochemical transformation, so why has nature chosen such complex solutions? In this review, we distill lessons from a century of metabolic research and introduce new observations suggesting that the intricate structure of metabolic pathways can be explained by a small set of biochemical principles. Using glycolysis as an example, we demonstrate how three key biochemical constraints--thermodynamic favorability, availability of enzymatic mechanisms and the physicochemical properties of pathway intermediates--eliminate otherwise plausible metabolic strategies. Considering these constraints, glycolysis contains no unnecessary steps and represents one of the very few pathway structures that meet cellular demands. The analysis presented here can be applied to metabolic engineering efforts for the rational design of pathways that produce a desired product while satisfying biochemical constraints.     

Metabolic engineering of Escherichia coli for direct production of 1,4-butanediol

1,4-Butanediol (BDO) is an important commodity chemical used to manufacture over 2.5 million tons annually of valuable polymers, and it is currently produced exclusively through feedstocks derived from oil and natural gas. Herein we report what are to our knowledge the first direct biocatalytic routes to BDO from renewable carbohydrate feedstocks, leading to a strain of Escherichia coli capable of producing 18 g l(-1) of this highly reduced, non-natural chemical. A pathway-identification algorithm elucidated multiple pathways for the biosynthesis of BDO from common metabolic intermediates. Guided by a genome-scale metabolic model, we engineered the E. coli host to enhance anaerobic operation of the oxidative tricarboxylic acid cycle, thereby generating reducing power to drive the BDO pathway. The organism produced BDO from glucose, xylose, sucrose and biomass-derived mixed sugar streams. This work demonstrates a systems-based metabolic engineering approach to strain design and development that can enable new bioprocesses for commodity chemicals that are not naturally produced by living cells.     

Engineering primary metabolic pathways of industrial micro-organisms

Metabolic engineering is a powerful tool for the optimisation and the introduction of new cellular processes. This is mostly done by genetic engineering. Since the introduction of this multidisciplinary approach, the success stories keep accumulating. The primary metabolism of industrial micro-organisms has been studied for long time and most biochemical pathways and reaction networks have been elucidated. This large pool of biochemical information, together with data from proteomics, metabolomics and genomics underpins the strategies for design of experiments and choice of targets for manipulation by metabolic engineers. These targets are often located in the primary metabolic pathways, such as glycolysis, pentose phosphate pathway, the TCA cycle and amino acid biosynthesis and mostly at major branch points within these pathways. This paper describes approaches taken for metabolic engineering of these pathways in bacteria, yeast and filamentous fungi.