Search

Cyanobacteria are photosynthetic prokaryotes with strong potential to be used for industrial terpenoid production. However, the key enzymes forming the principal terpenoid building blocks, called short-chain prenyltransferases (SPTs), are insufficiently characterized. Here, we examined SPTs in the model cyanobacteria Synechococcus elongatus sp. PCC 7942 and Synechocystis sp. PCC 6803. Each species has a single putative SPT (SeCrtE and SyCrtE, respectively). Sequence analysis identified these as type-II geranylgeranyl pyrophosphate synthases (GGPPSs) with high homology to GGPPSs found in the plastids of green plants and other photosynthetic organisms. In vitro analysis demonstrated that SyCrtE is multifunctional, producing geranylgeranyl pyrophosphate (GGPP; C

Triterpenoids are a subgroup of terpenoids and have wide applications in the food, cosmetics, and pharmaceutical industries. The heterologous production of various triterpenoids in

Triterpenoids, one of the most diverse classes of natural products, have been used for centuries as active ingredients in essential oils and Chinese medicines and are of interest for many industrial applications ranging from low-calorie sweeteners to cosmetic ingredients and vaccine adjuvants. However, not only can the extraction from plant material be cumbersome due to low concentrations of the specific triterpenoid, but concerns are also increasing regarding the sustainability of wild plant harvest while meeting market demands. The alternative is to produce triterpenoids with engineered microbes. Here, we present a generally applicable strategy for triterpenoid production in the yeast Saccharomyces cerevisiae based on a modified oxidosqualene cyclase Erg7. The modification reduces the flux into the sterol pathway while increasing the precursor supply for triterpenoid production. The minimally engineered strain was exploited for the exemplary production of the lupane triterpenoids betulin, betulin aldehyde, and betulinic acid at a total titer above 6 g/L, the highest reported so far. To further highlight the chassis concept, squalene, oleanane- and dammarane-type triterpenoids were synthesized to titers at a similar gram scale. We propose the developed baker’s yeast as a host for the thousands of triterpenoid synthesis pathways from plants, reducing the pressure on the natural resources.

Isoprenoids, also known as terpenes or terpenoids, are compounds made of one or more isoprene (C

The yeast

In metabolic engineering, loss-of-function experiments are used to understand and optimise metabolism. A conditional gene inactivation tool is required when gene deletion is lethal or detrimental to growth. Here, we exploit auxin-inducible protein degradation as a metabolic engineering approach in yeast. We demonstrate its effectiveness using terpenoid production. First, we target an essential prenyl-pyrophosphate metabolism protein, farnesyl pyrophosphate synthase (Erg20p). Degradation successfully redirects metabolic flux toward monoterpene (C10) production. Second, depleting hexokinase-2, a key protein in glucose signalling transduction, lifts glucose repression and boosts production of sesquiterpene (C15) nerolidol to 3.5 g L−1 in flask cultivation. Third, depleting acetyl-CoA carboxylase (Acc1p), another essential protein, delivers growth arrest without diminishing production capacity in nerolidol-producing yeast, providing a strategy to decouple growth and production. These studies demonstrate auxin-mediated protein degradation as an advanced tool for metabolic engineering. It also has potential for broader metabolic perturbation studies to better understand metabolism., Loss-of-function experiments are used in metabolic engineering to understand and optimise metabolism. Here, the authors apply auxin inducible protein degradation to test different metabolic engineering strategies for improved terpenoid production in yeast.

Methyl ketones present a group of highly reduced platform chemicals industrially produced from petroleum-derived hydrocarbons. They find applications in the fragrance, flavor, pharmacological, and agrochemical industries, and are further discussed as biodiesel blends. In recent years, intense research has been carried out to achieve sustainable production of these molecules by re-arranging the fatty acid metabolism of various microbes. One challenge in the development of a highly productive microbe is the high demand for reducing power. Here, we engineeredPseudomonas taiwanensisVLB120 for methyl ketone production as this microbe has been shown to sustain exceptionally high NAD(P)H regeneration rates. The implementation of published strategies resulted in 2.1 g Laq-1methyl ketones in fed-batch fermentation. We further increased the production by eliminating competing reactions suggested by metabolic analyses. These efforts resulted in the production of 9.8 g Laq-1methyl ketones (corresponding to 69.3 g Lorg-1in thein situextraction phase) at 53 % of the maximum theoretical yield. This represents a 4-fold improvement in product titer compared to the initial production strain and the highest titer of recombinantly produced methyl ketones reported to date. Accordingly, this study underlines the high potential ofP. taiwanensisVLB120 to produce methyl ketones and emphasizes model-driven metabolic engineering to rationalize and accelerate strain optimization efforts.

Methyl ketones present a group of highly reduced platform chemicals industrially produced from petroleum-derived hydrocarbons. They find applications in the fragrance, flavor, pharmacological, and agrochemical industries, and are further discussed as biodiesel blends. In recent years, intense research has been carried out to achieve sustainable production of these molecules by re-arranging the fatty acid metabolism of various microbes. One challenge in the development of a highly productive microbe is the high demand for reducing power. Here, we engineered Pseudomonas taiwanensis VLB120 for methyl ketone production as this microbe has been shown to sustain exceptionally high NAD(P)H regeneration rates. The implementation of published strategies resulted in 2.1 g Laq−1 methyl ketones in fed-batch fermentation. We further increased the production by eliminating competing reactions suggested by metabolic analyses. These efforts resulted in the production of 9.8 g Laq−1 methyl ketones (corresponding to 69.3 g Lorg−1 in the in situ extraction phase) at 53% of the maximum theoretical yield. This represents a 4-fold improvement in product titer compared to the initial production strain and the highest titer of recombinantly produced methyl ketones reported to date. Accordingly, this study underlines the high potential of P. taiwanensis VLB120 to produce methyl ketones and emphasizes model-driven metabolic engineering to rationalize and accelerate strain optimization efforts.

Microbial production of plant derived, biologically active compounds has the potential to provide economic and ecologic alternatives to existing low productive, plant-based processes. Current production of the pharmacologically active cyclic triterpenoid betulinic acid is realized by extraction from the bark of plane tree or birch. Here, we reengineered the reported betulinic acid pathway into S. cerevisiae and used this novel strain to develop efficient fermentation and product purification methods. Fed-batch cultivations with ethanol excess, using either an ethanol-pulse feed or controlling a constant ethanol concentration in the fermentation medium, significantly enhanced production of betulinic acid and its triterpenoid precursors. The beneficial effect of excess ethanol was further exploited in nitrogen-limited resting cell fermentations, yielding betulinic acid concentrations of 182 mg/L and total triterpenoid concentrations of 854 mg/L, the highest concentrations reported so far. Purification of lupane-type triterpenoids with high selectivity and yield was achieved by solid-liquid extraction without prior cell disruption using polar aprotic solvents such as acetone or ethyl acetate and subsequent precipitation with strong acids This study highlights the potential of microbial production of plant derived triterpenoids in S. cerevisiae by combining metabolic and process engineering. This article is protected by copyright. All rights reserved.

BMC bioinformatics 20(1), 447 (2019). doi:10.1186/s12859-019-2946-7, Published by Springer, Heidelberg

Exploring the dynamic behavior of cellular metabolism requires a standard laboratory method that guarantees rapid sampling and extraction of the cellular content. We propose a versatile sampling technique applicable to cells with different cell wall and cell membrane properties. The technique is based on irreversible electroporation with simultaneous quenching and extraction by using a microfluidic device. By application of electric pulses in the millisecond range, permanent lethal pores are formed in the cell membrane of Escherichia coli and Saccharomyces cerevisiae, facilitating the release of the cellular contents, here demonstrated by the measurement of glucose-6-phosphate and the activity of the enzyme glucose-6-phosphate dehydrogenase. The successful application of this device was demonstrated by pulsed electric field treatment in a flow-through configuration of the microfluidic chip in combination with sampling, inactivation, and extraction of the intracellular content in a few seconds. Minimum electric field strengths of 10 kV/cm for E. coli and 7.5 kV/cm for yeast S. cerevisiae were required for successful cell lysis. The results are discussed in the context of applications in industrial biotechnology, where metabolomics analyses are important.

Pseudomonas putida efficiently utilizes many different carbon sources without the formation of byproducts even under conditions of stress. This implies a high degree of flexibility to cope with conditions that require a significantly altered distribution of carbon to either biomass or energy in the form of NADH. In the literature, co‐feeding of the reduced C1 compound formate to Escherichia coli heterologously expressing the NAD(+)‐dependent formate dehydrogenase of the yeast Candida boidinii was demonstrated to boost various NADH‐demanding applications. Pseudomonas putida as emerging biotechnological workhorse is inherently equipped with an NAD(+)‐dependent formate dehydrogenase encouraging us to investigate the use of formate and its effect on P. putida’s metabolism. Hence, this study provides a detailed insight into the co‐utilization of formate and glucose by P. putida. Our results show that the addition of formate leads to a high increase in the NADH regeneration rate resulting in a very high biomass yield on glucose. Metabolic flux analysis revealed a significant flux rerouting from catabolism to anabolism. These metabolic insights argue further for P. putida as a host for redox cofactor demanding bioprocesses.

Heterologous synthesis of triterpenoids in Saccharomyces cerevisiae from its native metabolite squalene has been reported to offer an alternative to chemical synthesis and extraction from plant material if productivities can be increased.Here, we physiologically characterized a squalene overproducing S. cerevisiae CEN.PK strain to elucidate the effect of cultivation conditions on the production of this central triterpenoid precursor. The maximum achievable squalene concentration was substantially influenced by nutritional conditions, medium composition and cultivation mode. Batch growth on glucose resulted in minimal squalene accumulation, while squalene only significantly accumulated during ethanol consumption (up to 59 mg/gCDW), probably due to increased acetyl-CoA supply on this carbon source. Likewise, low squalene concentrations were observed in glucose-limited chemostat cultivations and improved up to 8-fold upon increasing the ethanol fraction in the feed. In those experiments, a constant, growth-rate-independent specific squalene accumulation rate (2.2 mg/gCDW/h) was recorded resulting in a maximal squalene loading of 30 mg/gCDW at low dilution rates with longer residence times. Coenzyme A availability was identified as possible bottleneck as increased vitamin concentrations, including the Coenzyme A precursor pantothenate, improved squalene titers in batch and chemostat cultivations. This analysis demonstrates that thorough physiologic characterization of production strains is valuable for the identification of bottlenecks already in early stages of strain development and for guiding further optimization efforts.

Rational strain engineering requires solid testing of phenotypes including productivity and ideally contributes thereby directly to our understanding of the genotype-phenotype relationship. Actually, the test step of the strain engineering cycle becomes the limiting step, as ever advancing tools for generating genetic diversity exist. Here, we briefly define the challenge one faces in quantifying phenotypes and summarize existing analytical techniques that partially overcome this challenge. We argue that the evolution of volatile metabolites can be used as proxy for cellular metabolism. In the simplest case, the product of interest is a volatile (e.g., from bulk alcohols to special fragrances) that is directly quantified over time. But also nonvolatile products (e.g., from bulk long-chain fatty acids to natural products) require major flux rerouting that result potentially in altered volatile production. While alternative techniques for volatile determination exist, rather few can be envisaged for medium to high-throughput analysis required for phenotype testing. Here, we contribute a detailed protocol for an ion mobility spectrometry (IMS) analysis that allows volatile metabolite quantification down to the ppb range. The sensitivity can be exploited for small-scale fermentation monitoring. The insights shared might contribute to a more frequent use of IMS in biotechnology, while the experimental aspects are of general use for researchers interested in volatile monitoring.

Volatile organic compounds (VOCs) are small molecular mass substances, which exhibit low boiling points and high-vapour pressures. They are ubiquitous in nature and produced by almost any organism of all kingdoms of life. VOCs are involved in many inter- and intraspecies interactions ranging from antimicrobial or fungal effects to plant growth promotion and human taste perception of fermentation products. VOC profiles further reflect the metabolic or phenotypic state of the living organism that produces them. Hence, they can be exploited for non-invasive medicinal diagnoses or industrial fermentation control. Here, we introduce the reader to these diverse applications associated with the monitoring and analysis of VOC emissions. We also present our vision of real-time VOC analysis enabled by newly developed analytical techniques, which will further broaden the use of VOCs in even wider applications. Hence, we foresee a bright future for VOC research and its associated fields of applications.

Rational strain engineering requires solid testing of phenotypes including productivity and ideally contributes thereby directly to our understanding of the genotype-phenotype relationship. Actually, the test step of the strain engineering cycle becomes the limiting step, as ever advancing tools for generating genetic diversity exist. Here, we briefly define the challenge one faces in quantifying phenotypes and summarize existing analytical techniques that partially overcome this challenge. We argue that the evolution of volatile metabolites can be used as proxy for cellular metabolism. In the simplest case, the product of interest is a volatile (e.g., from bulk alcohols to special fragrances) that is directly quantified over time. But also nonvolatile products (e.g., from bulk long-chain fatty acids to natural products) require major flux rerouting that result potentially in altered volatile production. While alternative techniques for volatile determination exist, rather few can be envisaged for medium to high-throughput analysis required for phenotype testing. Here, we contribute a detailed protocol for an ion mobility spectrometry (IMS) analysis that allows volatile metabolite quantification down to the ppb range. The sensitivity can be exploited for small-scale fermentation monitoring. The insights shared might contribute to a more frequent use of IMS in biotechnology, while the experimental aspects are of general use for researchers interested in volatile monitoring.

Background: Volatile organic compounds (VOCs) are small molecular mass substances, which exhibit high-vapor pressures, low boiling points, and lipophilic character. VOCs are produced by all organisms including eukaryotic microbes like yeast, whose volatile metabolites are for centuries exploited for examples as flavors in bread, beer, and wine. Notably, while the applications of VOCs are many, the knowledge on their biochemical synthesis is still limited.Objective: We review here the current information of yeast volatile metabolites and techniques to further explore the VOC landscape made possible by improvements of the analytical possibilities, regarding sampling frequency, identification, and quantification and the development to computationally interpret (high-throughput) data. Especially possibilities for online and even real-time analysis should trigger new experimental approaches that elucidate the biochemistry as well as the regulation of VOC synthesis. Baker's yeast is here the organism of choice as the genetic inventory can be linked to VOC formation and with this in hand improved applications can be envisaged. The physical, chemical or biological properties make many VOCs interesting targets for different industrial sectors while their natural function as semiochemicals or in defense mechanisms can be exploited to engineer synthetic microbial communities or to develop new antibiotics.Conclusion: VOCs produced by microbes including yeast are a chemical diverse group of compounds with highly different applications. The new analytical techniques briefly summarized here will enable the use of VOCs in even broader applications including human health monitoring and bioprocess control. We envisage a bright future for VOC research and for the resulting applications.

The intracellular reaction rates (fluxes) are the ultimate outcome of the activities of the complete inventory (from DNA to metabolite) and in their sum determine the cellular phenotype. The genotype-phenotype relationship is fundamental in such different fields as cancer research and biotechnology. Here, we summarize the developments in determining metabolic fluxes, inferring major pathways from the DNA-sequence, estimating optimal flux distributions, and how these flux distributions can be achieved in vivo. The technical advances to intervene with the many levels of the cellular architecture allow the implementation of new strategies in for example Metabolic Engineering.

Whole-cell biocatalysis utilizes native or recombinant enzymes produced by cellular metabolism to perform synthetically interesting reactions. Besides hydrolases, oxidoreductases represent the most applied enzyme class in industry. Oxidoreductases are attributed a high future potential, especially for applications in the chemical and pharmaceutical industries, as they enable highly interesting chemistry (e.g., the selective oxyfunctionalization of unactivated C-H bonds). Redox reactions are characterized by electron transfer steps that often depend on redox cofactors as additional substrates. Their regeneration typically is accomplished via the metabolism of whole-cell catalysts. Traditionally, studies towards productive redox biocatalysis focused on the biocatalytic enzyme, its activity, selectivity, and specificity, and several successful examples of such processes are running commercially. However, redox cofactor regeneration by host metabolism was hardly considered for the optimization of biocatalytic rate, yield, and/or titer. This article reviews molecular mechanisms of oxidoreductases with synthetic potential and the host redox metabolism that fuels biocatalytic reactions with redox equivalents. The tools discussed in this review for investigating redox metabolism provide the basis for studies aiming at a deeper understanding of the interplay between synthetically active enzymes and metabolic networks. The ultimate goal of rational whole-cell biocatalyst engineering and use for fine chemical production is discussed.

A key limitation of whole-cell redox biocatalysis for the production of valuable, specifically functionalized products is substrate/product toxicity, which can potentially be overcome by using solvent-tolerant micro-organisms. To investigate the inter-relationship of solvent tolerance and energy-dependent biocatalysis, we established a model system for biocatalysis in the presence of toxic low logPow solvents: recombinant solvent-tolerant Pseudomonas putida DOT-T1E catalyzing the stereospecific epoxidation of styrene in an aqueous/octanol two-liquid phase reaction medium. Using 13C tracer based metabolic flux analysis, we investigated the central carbon and energy metabolism and quantified the NAD(P)H regeneration rate in the presence of toxic solvents and during redox biocatalysis, which both drastically increased the energy demands of solvent-tolerant P. putida. According to the driven by demand concept, the NAD(P)H regeneration rate was increased up to eightfold by two mechanisms: (a) an increase in glucose uptake rate without secretion of metabolic side products, and (b) reduced biomass formation. However, in the presence of octanol, only ∼ 1% of the maximally observed NAD(P)H regeneration rate could be exploited for styrene epoxidation, of which the rate was more than threefold lower compared with operation with a non-toxic solvent. This points to a high energy and redox cofactor demand for cell maintenance, which limits redox biocatalysis in the presence of octanol. An estimated upper bound for the NAD(P)H regeneration rate available for biocatalysis suggests that cofactor availability does not limit redox biocatalysis under optimized conditions, for example, in the absence of toxic solvent, and illustrates the high metabolic capacity of solvent-tolerant P. putida. This study shows that solvent-tolerant P. putida have the remarkable ability to compensate for high energy demands by boosting their energy metabolism to levels up to an order of magnitude higher than those observed during unlimited growth.

(13)C label-based metabolic flux analysis is a powerful technique for the determination of intracellular reaction rates and is used in such different research fields as quantitative physiology, metabolic engineering, and systems biology. Metabolic fluxes can be determined at high quality using (pseudo)-steady-state cultures and advanced mathematical models for data interpretation. Here, we describe a protocol for parallel metabolic flux analysis that consists of downsized microbial (yeast) cultivation, miniaturized sample preparation, and semiautomated analytics and data evaluation. With this protocol dozens of metabolic flux analyses can be carried out in 1 week, thereby enabling for example the analysis of genetic and environmental perturbations on the operation of metabolic networks.

(13)C label-based metabolic flux analysis is a powerful technique for the determination of intracellular reaction rates and is used in such different research fields as quantitative physiology, metabolic engineering, and systems biology. Metabolic fluxes can be determined at high quality using (pseudo)-steady-state cultures and advanced mathematical models for data interpretation. Here, we describe a protocol for parallel metabolic flux analysis that consists of downsized microbial (yeast) cultivation, miniaturized sample preparation, and semiautomated analytics and data evaluation. With this protocol dozens of metabolic flux analyses can be carried out in 1 week, thereby enabling for example the analysis of genetic and environmental perturbations on the operation of metabolic networks.

Adenosine phosphate and NAD cofactors play a vital role in the operation of cell metabolism, and their levels and ratios are carefully regulated in tight ranges. Perturbations of the consumption of these metabolites might have a great impact on cell metabolism and physiology. Here, we investigated the impact of increased ATP hydrolysis and NADH oxidation rates on the metabolism of Pseudomonas putida KT2440 by titration of 2,4-dinitrophenol (DNP) and overproduction of a water-forming NADH oxidase, respectively. Both perturbations resulted in a reduction of the biomass yield and, as a consequence of the uncoupling of catabolic and anabolic activities, in an amplification of the net NADH regeneration rate. However, a stimulation of the specific carbon uptake rate was observed only when P. putida was challenged with very high 2,4-dinitrophenol concentrations and was comparatively unaffected by recombinant NADH oxidase activity. This behavior contrasts with the comparably sensitive performance described, for example, for Escherichia coli or Saccharomyces cerevisiae . The apparent robustness of P. putida metabolism indicates that it possesses a certain buffering capacity and a high flexibility to adapt to and counteract different stresses without showing a distinct phenotype. These findings are important, e.g., for the development of whole-cell redox biocatalytic processes that impose equivalent burdens on the cell metabolism: stoichiometric consumption of (reduced) redox cofactors and increased energy expenditures, due to the toxicity of the biocatalytic compounds.

A key limitation of whole-cell redox biocatalysis for the production of valuable, specifically functionalized products is substrate/product toxicity, which can potentially be overcome by using solvent-tolerant micro-organisms. To investigate the inter-relationship of solvent tolerance and energy-dependent biocatalysis, we established a model system for biocatalysis in the presence of toxic low logP(ow) solvents: recombinant solvent-tolerant Pseudomonas putida DOT-T1E catalyzing the stereospecific epoxidation of styrene in an aqueous/octanol two-liquid phase reaction medium. Using (13)C tracer based metabolic flux analysis, we investigated the central carbon and energy metabolism and quantified the NAD(P)H regeneration rate in the presence of toxic solvents and during redox biocatalysis, which both drastically increased the energy demands of solvent-tolerant P. putida. According to the driven by demand concept, the NAD(P)H regeneration rate was increased up to eightfold by two mechanisms: (a) an increase in glucose uptake rate without secretion of metabolic side products, and (b) reduced biomass formation. However, in the presence of octanol, only approximately 1% of the maximally observed NAD(P)H regeneration rate could be exploited for styrene epoxidation, of which the rate was more than threefold lower compared with operation with a non-toxic solvent. This points to a high energy and redox cofactor demand for cell maintenance, which limits redox biocatalysis in the presence of octanol. An estimated upper bound for the NAD(P)H regeneration rate available for biocatalysis suggests that cofactor availability does not limit redox biocatalysis under optimized conditions, for example, in the absence of toxic solvent, and illustrates the high metabolic capacity of solvent-tolerant P. putida. This study shows that solvent-tolerant P. putida have the remarkable ability to compensate for high energy demands by boosting their energy metabolism to levels up to an order of magnitude higher than those observed during unlimited growth.

Whole-cell redox biocatalysis relies on redox cofactor regeneration by the microbial host. Here, we applied flux balance analysis based on the Escherichia coli metabolic network to estimate maximal NADH regeneration rates. With this optimization criterion, simulations showed exclusive use of the pentose phosphate pathway at high rates of glucose catabolism, a flux distribution usually not found in wild-type cells. In silico, genetic perturbations indicated a strong dependency of NADH yield and formation rate on the underlying metabolic network structure. The linear dependency of measured epoxidation activities of recombinant central carbon metabolism mutants on glucose uptake rates and the linear correlation between measured activities and simulated NADH regeneration rates imply intracellular NADH shortage. Quantitative comparison of computationally predicted NADH regeneration and experimental epoxidation rates indicated that the achievable biocatalytic activity is determined by metabolic and enzymatic limitations including non-optimal flux distributions, high maintenance energy demands, energy spilling, byproduct formation, and uncoupling. The results are discussed in the context of cellular optimization of biotransformation processes and may guide a priori design of microbial cells as redox biocatalysts.

The production of biofuels must be economically and ecologically viable, in particular due to the competition with existing mature fossil-based fuels. There are many alternative products and pathways to generate biofuels. Hence, focusing on the most promising products and process alternatives requires identification of corresponding pathways from biomass to biofuel candidates using scarce data of high uncertainty. This contribution presents a methodology for i) process evaluation of fermentation and downstream processing (DP) in terms of costs and primary energy demand (PED) as main sustainability criteria and ii) identification of the most promising platform chemicals gained by fermentation for a biofuel production. The methodology considers the energy requirement already at an early design stage bridging the gap between performance screenings solely based on reaction stoichiometry and time-consuming process design, enabling an early process analysis, detection of bottlenecks and ranking of various processes. The focus herein is on lignocellulosic biomass.

The limited supply of fossil fuels and the establishment of new environmental policies shifted research in industry and academia towards sustainable production of the 2ndgeneration of biofuels, with Methyl Ethyl Ketone (MEK) being one promising fuel candidate. MEK is a commercially valuable petrochemical with an extensive application as a solvent. However, as of today, a sustainable and economically viable production of MEK has not yet been achieved despite several attempts of introducing biosynthetic pathways in industrial microorganisms. We used BNICE.ch as a retrobiosynthesis tool to discover all novel pathways around MEK. Out of 1’325 identified compounds connecting to MEK with one reaction step, we selected 3-oxopentanoate, but-3-en-2-one, but-1-en-2-olate, butylamine, and 2-hydroxy-2-methyl-butanenitrile for further study. We reconstructed 3’679’610 novel biosynthetic pathways towards these 5 compounds. We then embedded these pathways into the genome-scale model ofE. coli, and a set of 18’622 were found to be most biologically feasible ones based on thermodynamics and their yields. For each novel reaction in the viable pathways, we proposed the most similar KEGG reactions, with their gene and protein sequences, as candidates for either a direct experimental implementation or as a basis for enzyme engineering. Through pathway similarity analysis we classified the pathways and identified the enzymes and precursors that were indispensable for the production of the target molecules. These retrobiosynthesis studies demonstrate the potential of BNICE.ch for discovery, systematic evaluation, and analysis of novel pathways in synthetic biology and metabolic engineering studies.Graphical abstract