OptStrain: a computational framework for redesign of microbial production systems

MRE: a web tool to suggest foreign enzymes for the biosynthesis pathway design with competing endogenous reactions in mind

To rationally design a productive heterologous biosynthesis system, it is essential to consider the suitability of foreign reactions for the specific endogenous metabolic infrastructure of a host. We developed a novel web server, called MRE, which, for a given pair of starting and desired compounds in a given chassis organism, ranks biosynthesis routes from the perspective of the integration of new reactions into the endogenous metabolic system. For each promising heterologous biosynthesis pathway, MRE suggests actual enzymes for foreign metabolic reactions and generates information on competing endogenous reactions for the consumption of metabolites. These unique, chassis-centered features distinguish MRE from existing pathway design tools and allow synthetic biologists to evaluate the design of their biosynthesis systems from a different angle. By using biosynthesis of a range of high-value natural products as a case study, we show that MRE is an effective tool to guide the design and optimization of heterologous biosynthesis pathways. The URL of MRE is http://www.cbrc.kaust.edu.sa/mre/.

Photorespiratory Bypasses Lead to Increased Growth in Arabidopsis thaliana: Are Predictions Consistent with Experimental Evidence?

Arguably, the biggest challenge of modern plant systems biology lies in predicting the performance of plant species, and crops in particular, upon different intracellular and external perturbations. Recently, an increased growth of Arabidopsis thaliana plants was achieved by introducing two different photorespiratory bypasses via metabolic engineering. Here, we investigate the extent to which these findings match the predictions from constraint-based modeling. To determine the effect of the employed metabolic network model on the predictions, we perform a comparative analysis involving three state-of-the-art metabolic reconstructions of A. thaliana. In addition, we investigate three scenarios with respect to experimental findings on the ratios of the carboxylation and oxygenation reactions of Ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO). We demonstrate that the condition-dependent growth phenotypes of one of the engineered bypasses can be qualitatively reproduced by each reconstruction, particularly upon considering the additional constraints with respect to the ratio of fluxes for the RuBisCO reactions. Moreover, our results lend support for the hypothesis of a reduced photorespiration in the engineered plants, and indicate that specific changes in CO₂ exchange as well as in the proxies for co-factor turnover are associated with the predicted growth increase in the engineered plants. We discuss our findings with respect to the structure of the used models, the modeling approaches taken, and the available experimental evidence. Our study sets the ground for investigating other strategies for increase of plant biomass by insertion of synthetic reactions.

Synthetic and systems biology for microbial production of commodity chemicals

The combination of synthetic and systems biology is a powerful framework to study fundamental questions in biology and produce chemicals of immediate practical application such as biofuels, polymers, or therapeutics. However, we cannot yet engineer biological systems as easily and precisely as we engineer physical systems. In this review, we describe the path from the choice of target molecule to scaling production up to commercial volumes. We present and explain some of the current challenges and gaps in our knowledge that must be overcome in order to bring our bioengineering capabilities to the level of other engineering disciplines. Challenges start at molecule selection, where a difficult balance between economic potential and biological feasibility must be struck. Pathway design and construction have recently been revolutionized by next-generation sequencing and exponentially improving DNA synthesis capabilities. Although pathway optimization can be significantly aided by enzyme expression characterization through proteomics, choosing optimal relative protein expression levels for maximum production is still the subject of heuristic, non-systematic approaches. Toxic metabolic intermediates and proteins can significantly affect production, and dynamic pathway regulation emerges as a powerful but yet immature tool to prevent it. Host engineering arises as a much needed complement to pathway engineering for high bioproduct yields; and systems biology approaches such as stoichiometric modeling or growth coupling strategies are required. A final, and often underestimated, challenge is the successful scale up of processes to commercial volumes. Sustained efforts in improving reproducibility and predictability are needed for further development of bioengineering.

MapMaker and PathTracer for tracking carbon in genome-scale metabolic models

Constraint-based reconstruction and analysis (COBRA) modeling results can be difficult to interpret given the large numbers of reactions in genome-scale models. While paths in metabolic networks can be found, existing methods are not easily combined with constraint-based approaches. To address this limitation, two tools (MapMaker and PathTracer) were developed to find paths (including cycles) between metabolites, where each step transfers carbon from reactant to product. MapMaker predicts carbon transfer maps (CTMs) between metabolites using only information on molecular formulae and reaction stoichiometry, effectively determining which reactants and products share carbon atoms. MapMaker correctly assigned CTMs for over 97% of the 2,251 reactions in an Escherichia coli metabolic model (iJO1366). Using CTMs as inputs, PathTracer finds paths between two metabolites. PathTracer was applied to iJO1366 to investigate the importance of using CTMs and COBRA constraints when enumerating paths, to find active and high flux paths in flux balance analysis (FBA) solutions, to identify paths for putrescine utilization, and to elucidate a potential CO2 fixation pathway in E. coli. These results illustrate how MapMaker and PathTracer can be used in combination with constraint-based models to identify feasible, active, and high flux paths between metabolites.

Designing intracellular metabolism for production of target compounds by introducing a heterologous metabolic reaction based on a Synechosystis sp. 6803 genome-scale model

Background

Designing optimal intracellular metabolism is essential for using microorganisms to produce useful compounds. Computerized calculations for flux balance analysis utilizing a genome-scale model have been performed for such designs. Many genome-scale models have been developed for different microorganisms. However, optimal designs of intracellular metabolism aimed at producing a useful compound often utilize metabolic reactions of only the host microbial cells. In the present study, we added reactions other than the metabolic reactions with Synechosystis sp. 6803 as a host to its genome-scale model, and constructed a metabolic model of hybrid cells (SyHyMeP) using computerized analysis. Using this model provided a metabolic design that improves the theoretical yield of succinic acid, which is a useful compound.

Results

Constructing the SyHyMeP model enabled new metabolic designs for producing useful compounds. In the present study, we developed a metabolic design that allowed for improved theoretical yield in the production of succinic acid during glycogen metabolism by Synechosystis sp. 6803. The theoretical yield of succinic acid production using a genome-scale model of these cells was 1.00 mol/mol-glucose, but use of the SyHyMeP model enabled a metabolic design with which a 33 % increase in theoretical yield is expected due to the introduction of isocitrate lyase, adding activations of endogenous tree reactions via D-glycerate in Synechosystis sp. 6803.

Conclusions

The SyHyMeP model developed in this study has provided a new metabolic design that is not restricted only to the metabolic reactions of individual microbial cells. The concept of construction of this model requires only replacement of the genome-scale model of the host microbial cells and can thus be applied to various useful microorganisms for metabolic design to produce compounds.

Electronic supplementary material

The online version of this article (doi:10.1186/s12934-016-0416-8) contains supplementary material, which is available to authorized users.

Applications of Genome-Scale Metabolic Models in Biotechnology and Systems Medicine

Genome-scale metabolic models (GEMs) have become a popular tool for systems biology, and they have been used in many fields such as industrial biotechnology and systems medicine. Since more and more studies are being conducted using GEMs, they have recently received considerable attention. In this review, we introduce the basic concept of GEMs and provide an overview of their applications in biotechnology, systems medicine, and some other fields. In addition, we describe the general principle of the applications and analyses built on GEMs. The purpose of this review is to introduce the application of GEMs in biological analysis and to promote its wider use by biologists.

A systems biology approach to reconcile metabolic network models with application to Synechocystis sp. PCC 6803 for biofuel production

Metabolic network models can be optimized for the production of desired materials like biofuels.

Flux balance analysis of primary metabolism in the diatom Phaeodactylum tricornutum

Diatoms (Bacillarophyceae) are photosynthetic unicellular microalgae that have risen to ecological prominence in oceans over the past 30 million years. They are of interest as potential feedstocks for sustainable biofuels. Maximizing production of these feedstocks will require genetic modifications and an understanding of algal metabolism. These processes may benefit from genome-scale models, which predict intracellular fluxes and theoretical yields, as well as the viability of knockout and knock-in transformants. Here we present a genome-scale metabolic model of a fully sequenced and transformable diatom: Phaeodactylum tricornutum. The metabolic network was constructed using the P. tricornutum genome, biochemical literature, and online bioinformatic databases. Intracellular fluxes in P. tricornutum were calculated for autotrophic, mixotrophic and heterotrophic growth conditions, as well as knockout conditions that explore the in silico role of glycolytic enzymes in the mitochondrion. The flux distribution for lower glycolysis in the mitochondrion depended on which transporters for TCA cycle metabolites were included in the model. The growth rate predictions were validated against experimental data obtained using chemostats. Two published studies on this organism were used to validate model predictions for cyclic electron flow under autotrophic conditions, and fluxes through the phosphoketolase, glycine and serine synthesis pathways under mixotrophic conditions. Several gaps in annotation were also identified. The model also explored unusual features of diatom metabolism, such as the presence of lower glycolysis pathways in the mitochondrion, as well as differences between P. tricornutum and other photosynthetic organisms.

Metabolic engineering of a novel muconic acid biosynthesis pathway via 4-hydroxybenzoic acid in Escherichia coli

cis,cis-Muconic acid (MA) is a commercially important raw material used in pharmaceuticals, functional resins, and agrochemicals. MA is also a potential platform chemical for the production of adipic acid (AA), terephthalic acid, caprolactam, and 1,6-hexanediol. A strain of Escherichia coli K-12, BW25113, was genetically modified, and a novel nonnative metabolic pathway was introduced for the synthesis of MA from glucose. The proposed pathway converted chorismate from the aromatic amino acid pathway to MA via 4-hydroxybenzoic acid (PHB). Three nonnative genes, pobA, aroY, and catA, coding for 4-hydroxybenzoate hydrolyase, protocatechuate decarboxylase, and catechol 1,2-dioxygenase, respectively, were functionally expressed in E. coli to establish the MA biosynthetic pathway. E. coli native genes ubiC, aroF(FBR), aroE, and aroL were overexpressed and the genes ptsH, ptsI, crr, and pykF were deleted from the E. coli genome in order to increase the precursors of the proposed MA pathway. The final engineered E. coli strain produced nearly 170 mg/liter of MA from simple carbon sources in shake flask experiments. The proposed pathway was proved to be functionally active, and the strategy can be used for future metabolic engineering efforts for production of MA from renewable sugars.

In Silico Constraint-Based Strain Optimization Methods: the Quest for Optimal Cell Factories

Shifting from chemical to biotechnological processes is one of the cornerstones of 21st century industry. The production of a great range of chemicals via biotechnological means is a key challenge on the way toward a bio-based economy. However, this shift is occurring at a pace slower than initially expected. The development of efficient cell factories that allow for competitive production yields is of paramount importance for this leap to happen. Constraint-based models of metabolism, together with in silico strain design algorithms, promise to reveal insights into the best genetic design strategies, a step further toward achieving that goal. In this work, a thorough analysis of the main in silico constraint-based strain design strategies and algorithms is presented, their application in real-world case studies is analyzed, and a path for the future is discussed.

Designing overall stoichiometric conversions and intervening metabolic reactions

Existing computational tools for de novo metabolic pathway assembly, either based on mixed integer linear programming techniques or graph-search applications, generally only find linear pathways connecting the source to the target metabolite. The overall stoichiometry of conversion along with alternate co-reactant (or co-product) combinations is not part of the pathway design. Therefore, global carbon and energy efficiency is in essence fixed with no opportunities to identify more efficient routes for recycling carbon flux closer to the thermodynamic limit. Here, we introduce a two-stage computational procedure that both identifies the optimum overall stoichiometry (i.e., optStoic) and selects for (non-)native reactions (i.e., minRxn/minFlux) that maximize carbon, energy or price efficiency while satisfying thermodynamic feasibility requirements. Implementation for recent pathway design studies identified non-intuitive designs with improved efficiencies. Specifically, multiple alternatives for non-oxidative glycolysis are generated and non-intuitive ways of co-utilizing carbon dioxide with methanol are revealed for the production of C₂₊ metabolites with higher carbon efficiency.

Flux-sum analysis identifies metabolite targets for strain improvement

Background

Rational design of microbial strains for enhanced cellular physiology through in silico analysis has been reported in many metabolic engineering studies. Such in silico techniques typically involve the analysis of a metabolic model describing the metabolic and physiological states under various perturbed conditions, thereby identifying genetic targets to be manipulated for strain improvement. More often than not, the activation/inhibition of multiple reactions is necessary to produce a predicted change for improvement of cellular properties or states. However, as it is more computationally cumbersome to simulate all possible combinations of reaction perturbations, it is desirable to consider alternative techniques for identifying such metabolic engineering targets.

Results

In this study, we present the modified version of previously developed metabolite-centric approach, also known as flux-sum analysis (FSA), for identifying metabolic engineering targets. Utility of FSA was demonstrated by applying it to Escherichia coli, as case studies, for enhancing ethanol and succinate production, and reducing acetate formation. Interestingly, most of the identified metabolites correspond to gene targets that have been experimentally validated in previous works on E. coli strain improvement. A notable example is that pyruvate, the metabolite target for enhancing succinate production, was found to be associated with multiple reaction targets that were only identifiable through more computationally expensive means. In addition, detailed analysis of the flux-sum perturbed conditions also provided valuable insights into how previous metabolic engineering strategies have been successful in enhancing cellular physiology.

Conclusions

The application of FSA under the flux balance framework can identify novel metabolic engineering targets from the metabolite-centric perspective. Therefore, the current approach opens up a new research avenue for rational design and engineering of industrial microbes in the field of systems metabolic engineering.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-015-0198-3) contains supplementary material, which is available to authorized users.

Rational design of efficient modular cells

The modular cell design principle is formulated to devise modular (chassis) cells. These cells can be assembled with exchangeable production modules in a plug-and-play fashion to build microbial cell factories for efficient combinatorial biosynthesis of novel molecules, requiring minimal iterative strain optimization steps. A modular cell is designed to be auxotrophic, containing core metabolic pathways that are necessary but insufficient to support cell growth and maintenance. To be functional, it must tightly couple with an exchangeable production module containing auxiliary metabolic pathways that not only complement cell growth but also enhance production of targeted molecules. We developed a MODCELL (modular cell) framework based on metabolic pathway analysis to implement the modular cell design principle. MODCELL identifies genetic modifications and requirements to construct modular cell candidates and their associated exchangeable production modules. By defining the degree of similarity and coupling metrics, MODCELL can evaluate which exchangeable production module(s) can be tightly coupled with a modular cell candidate. We first demonstrated how MODCELL works in a step-by-step manner for example metabolic networks, and then applied it to design modular Escherichia coli cells for efficient combinatorial biosynthesis of five alcohols (ethanol, propanol, isopropanol, butanol and isobutanol) and five butyrate esters (ethyl butyrate, propyl butyrate, isopropyl butyrate, butyl butyrate and isobutyl butyrate) from pentose sugars (arabinose and xylose) and hexose sugars (glucose, mannose, and galactose) under anaerobic conditions. We identified three modular cells, MODCELL1, MODCELL2 and MODCELL3, that can couple well with Group 1 of modules (ethanol, isobutanol, butanol, ethyl butyrate, isobutyl butyrate, butyl butyrate), Group 2 (isopropanol, isopropyl butyrate), and Group 3 (propanol, isopropanol), respectively. We validated the design of MODCELL1 for anaerobic production of ethanol, butanol, and ethyl butyrate using experimental data available in literature.

FogLight: an efficient matrix-based approach to construct metabolic pathways by search space reduction

MOTIVATION

A fundamental computational problem in the area of metabolic engineering is finding metabolic pathways between a pair of source and target metabolites efficiently. We present an approach, namely FogLight, for searching metabolic networks utilizing Boolean (AND-OR) operations represented in matrix notation to efficiently reduce the search space. This enables the enumeration of all pathways between metabolites that are too distant for the application of brute-force methods.

RESULTS

Benchmarking tests run with FogLight show that it can reduce the search space by up to 98%, after which the accelerated search for high accurate results is guaranteed. Using FogLight, several pathways between eight given pairs of metabolites are found of which the pathways from CO2 to ethanol are specifically discussed. Additionally, in comparison with three path-finding tools, namely PHT, FMM and RouteSearch, FogLight can find shorter and more pathways for attempted source-target metabolite pairs.

CONTACT

szamani@aut.ac.ir, gholamreza.bidkhori@vtt.fi

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Cell scale host-pathogen modeling: another branch in the evolution of constraint-based methods

Constraint-based models have become popular methods for systems biology as they enable the integration of complex, disparate datasets in a biologically cohesive framework that also supports the description of biological processes in terms of basic physicochemical constraints and relationships. The scope, scale, and application of genome scale models have grown from single cell bacteria to multi-cellular interaction modeling; host-pathogen modeling represents one of these examples at the current horizon of constraint-based methods. There are now a small number of examples of host-pathogen constraint-based models in the literature, however there has not yet been a definitive description of the methodology required for the functional integration of genome scale models in order to generate simulation capable host-pathogen models. Herein we outline a systematic procedure to produce functional host-pathogen models, highlighting steps which require debugging and iterative revisions in order to successfully build a functional model. The construction of such models will enable the exploration of host-pathogen interactions by leveraging the growing wealth of omic data in order to better understand mechanism of infection and identify novel therapeutic strategies.

Exhaustive Analysis of a Genotype Space Comprising 10(15 )Central Carbon Metabolisms Reveals an Organization Conducive to Metabolic Innovation

All biological evolution takes place in a space of possible genotypes and their phenotypes. The structure of this space defines the evolutionary potential and limitations of an evolving system. Metabolism is one of the most ancient and fundamental evolving systems, sustaining life by extracting energy from extracellular nutrients. Here we study metabolism’s potential for innovation by analyzing an exhaustive genotype-phenotype map for a space of 10¹⁵ metabolisms that encodes all possible subsets of 51 reactions in central carbon metabolism. Using flux balance analysis, we predict the viability of these metabolisms on 10 different carbon sources which give rise to 1024 potential metabolic phenotypes. Although viable metabolisms with any one phenotype comprise a tiny fraction of genotype space, their absolute numbers exceed 10⁹ for some phenotypes. Metabolisms with any one phenotype typically form a single network of genotypes that extends far or all the way through metabolic genotype space, where any two genotypes can be reached from each other through a series of single reaction changes. The minimal distance of genotype networks associated with different phenotypes is small, such that one can reach metabolisms with novel phenotypes – viable on new carbon sources – through one or few genotypic changes. Exceptions to these principles exist for those metabolisms whose complexity (number of reactions) is close to the minimum needed for viability. Increasing metabolic complexity enhances the potential for both evolutionary conservation and evolutionary innovation.

Genotype-phenotype mapping is one of the ultimate goals of computational systems biology, and can provide new insights into the function and evolution of biological systems. We present a comprehensive genotype-phenotype map for a space of metabolic genotypes that comprises more than 10¹⁵ central carbon metabolisms. Only one in a million of these metabolisms can sustain life on any one of 10 carbon sources we consider, but these viable metabolisms form connected genotype networks that extend far through genotype space. In addition, they render multiple novel metabolic phenotypes in their immediate neighborhood accessible through small evolutionary changes that require only the alteration of single metabolic reactions. The map we construct reveals an organization of core metabolism that simultaneously facilitates evolutionary conservation of existing metabolic phenotypes, and the origination of novel metabolic traits that allow viability on novel carbon sources. Such metabolic innovation is essential, particularly for organisms that experience unexpected environmental changes, and that explore or invade new habitats.

Computational Methods for Modification of Metabolic Networks

In metabolic engineering, modification of metabolic networks is an important biotechnology and a challenging computational task. In the metabolic network modification, we should modify metabolic networks by newly adding enzymes or/and knocking-out genes to maximize the biomass production with minimum side-effect. In this mini-review, we briefly review constraint-based formalizations for Minimum Reaction Cut (MRC) problem where the minimum set of reactions is deleted so that the target compound becomes non-producible from the view point of the flux balance analysis (FBA), elementary mode (EM), and Boolean models. Minimum Reaction Insertion (MRI) problem where the minimum set of reactions is added so that the target compound newly becomes producible is also explained with a similar formalization approach. The relation between the accuracy of the models and the risk of overfitting is also discussed.

Co-evolution of strain design methods based on flux balance and elementary mode analysis

More than a decade ago, the first genome-scale metabolic models for two of the most relevant microbes for biotechnology applications, Escherichia coli and Saccaromyces cerevisiae, were published. Shortly after followed the publication of OptKnock, the first strain design method using bilevel optimization to couple cellular growth with the production of a target product. This initiated the development of a family of strain design methods based on the concept of flux balance analysis. Another family of strain design methods, based on the concept of elementary mode analysis, has also been growing. Although the computation of elementary modes is hindered by computational complexity, recent breakthroughs have allowed applying elementary mode analysis at the genome scale. Here we review and compare strain design methods and look back at the last 10 years of in silico strain design with constraint-based models. We highlight some features of the different approaches and discuss the utilization of these methods in successful in vivo metabolic engineering applications.

•

Computational strain design methods are divided into two main families.

•

We trace the evolutionary history of these two families.

•

Surveyed successful cases of model-guided strain design for industrial applications.

•

Most proposed methods have not yet been tested in real applications.

•

Agreement between in silico and in vivo results shows potential of tested methods.

Data-driven integration of genome-scale regulatory and metabolic network models

Microbes are diverse and extremely versatile organisms that play vital roles in all ecological niches. Understanding and harnessing microbial systems will be key to the sustainability of our planet. One approach to improving our knowledge of microbial processes is through data-driven and mechanism-informed computational modeling. Individual models of biological networks (such as metabolism, transcription, and signaling) have played pivotal roles in driving microbial research through the years. These networks, however, are highly interconnected and function in concert-a fact that has led to the development of a variety of approaches aimed at simulating the integrated functions of two or more network types. Though the task of integrating these different models is fraught with new challenges, the large amounts of high-throughput data sets being generated, and algorithms being developed, means that the time is at hand for concerted efforts to build integrated regulatory-metabolic networks in a data-driven fashion. In this perspective, we review current approaches for constructing integrated regulatory-metabolic models and outline new strategies for future development of these network models for any microbial system.

Genome-scale strain designs based on regulatory minimal cut sets

MOTIVATION

Stoichiometric and constraint-based methods of computational strain design have become an important tool for rational metabolic engineering. One of those relies on the concept of constrained minimal cut sets (cMCSs). However, as most other techniques, cMCSs may consider only reaction (or gene) knockouts to achieve a desired phenotype.

RESULTS

We generalize the cMCSs approach to constrained regulatory MCSs (cRegMCSs), where up/downregulation of reaction rates can be combined along with reaction deletions. We show that flux up/downregulations can virtually be treated as cuts allowing their direct integration into the algorithmic framework of cMCSs. Because of vastly enlarged search spaces in genome-scale networks, we developed strategies to (optionally) preselect suitable candidates for flux regulation and novel algorithmic techniques to further enhance efficiency and speed of cMCSs calculation. We illustrate the cRegMCSs approach by a simple example network and apply it then by identifying strain designs for ethanol production in a genome-scale metabolic model of Escherichia coli. The results clearly show that cRegMCSs combining reaction deletions and flux regulations provide a much larger number of suitable strain designs, many of which are significantly smaller relative to cMCSs involving only knockouts. Furthermore, with cRegMCSs, one may also enable the fine tuning of desired behaviours in a narrower range. The new cRegMCSs approach may thus accelerate the implementation of model-based strain designs for the bio-based production of fuels and chemicals.

AVAILABILITY AND IMPLEMENTATION

MATLAB code and the examples can be downloaded at http://www.mpi-magdeburg.mpg.de/projects/cna/etcdownloads.html.

CONTACT

krishna.mahadevan@utoronto.ca or klamt@mpi-magdeburg.mpg.de

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Applications of genome-scale metabolic network model in metabolic engineering

Genome-scale metabolic network model (GEM) is a fundamental framework in systems metabolic engineering. GEM is built upon extensive experimental data and literature information on gene annotation and function, metabolites and enzymes so that it contains all known metabolic reactions within an organism. Constraint-based analysis of GEM enables the identification of phenotypic properties of an organism and hypothesis-driven engineering of cellular functions to achieve objectives. Along with the advances in omics, high-throughput technology and computational algorithms, the scope and applications of GEM have substantially expanded. In particular, various computational algorithms have been developed to predict beneficial gene deletion and amplification targets and used to guide the strain development process for the efficient production of industrially important chemicals. Furthermore, an Escherichia coli GEM was integrated with a pathway prediction algorithm and used to evaluate all possible routes for the production of a list of commodity chemicals in E. coli. Combined with the wealth of experimental data produced by high-throughput techniques, much effort has been exerted to add more biological contexts into GEM through the integration of omics data and regulatory network information for the mechanistic understanding and improved prediction capabilities. In this paper, we review the recent developments and applications of GEM focusing on the GEM-based computational algorithms available for microbial metabolic engineering.

Development and application of efficient pathway enumeration algorithms for metabolic engineering applications

Metabolic Engineering (ME) aims to design microbial cell factories towards the production of valuable compounds. In this endeavor, one important task relates to the search for the most suitable heterologous pathway(s) to add to the selected host. Different algorithms have been developed in the past towards this goal, following distinct approaches spanning constraint-based modeling, graph-based methods and knowledge-based systems based on chemical rules. While some of these methods search for pathways optimizing specific objective functions, here the focus will be on methods that address the enumeration of pathways that are able to convert a set of source compounds into desired targets and their posterior evaluation according to different criteria. Two pathway enumeration algorithms based on (hyper)graph-based representations are selected as the most promising ones and are analyzed in more detail: the Solution Structure Generation and the Find Path algorithms. Their capabilities and limitations are evaluated when designing novel heterologous pathways, by applying these methods on three case studies of synthetic ME related to the production of non-native compounds in E. coli and S. cerevisiae: 1-butanol, curcumin and vanillin. Some targeted improvements are implemented, extending both methods to address limitations identified that impair their scalability, improving their ability to extract potential pathways over large-scale databases. In all case-studies, the algorithms were able to find already described pathways for the production of the target compounds, but also alternative pathways that can represent novel ME solutions after further evaluation.

Bridging the gap: a roadmap to breaking the biological design barrier

This paper presents an analysis of an emerging bottleneck in organism engineering, and paths by which it may be overcome. Recent years have seen the development of a profusion of synthetic biology tools, largely falling into two categories: high-level “design” tools aimed at mapping from organism specifications to nucleic acid sequences implementing those specifications, and low-level “build and test” tools aimed at faster, cheaper, and more reliable fabrication of those sequences and assays of their behavior in engineered biological organisms. Between the two families, however, there is a major gap: we still largely lack the predictive models and component characterization data required to effectively determine which of the many possible candidate sequences considered in the design phase are the most likely to produce useful results when built and tested. As low-level tools continue to mature, the bottleneck in biological systems engineering is shifting to be dominated by design, making this gap a critical barrier to progress. Considering how to address this gap, we find that widespread adoption of readily available analytic and assay methods is likely to lead to rapid improvement in available predictive models and component characterization models, as evidenced by a number of recent results. Such an enabling development is, in turn, likely to allow high-level tools to break the design barrier and support rapid development of transformative biological applications.

Contemplating 3-Hydroxypropionic Acid Biosynthesis in Klebsiella pneumoniae

3-Hydroxypropionic acid (3-HP) is a commercially valuable platform chemical from which an array of C3 compounds can be generated. Klebsiella pneumoniae has been considered a promising species for biological production of 3-HP. Despite a wealth of reports related to 3-HP biosynthesis in K. pneumoniae, its commercialization is still in infancy. The major hurdle hindering 3-HP overproduction lies in the poor understanding of glycerol dissimilation in K. pneumoniae. To surmount this problem, this review aims to portray a picture of 3-HP biosynthesis, involving 3-HP-synthesizing strains, biochemical attributes, metabolic pathways and key enzymes. Inspired by the state-of-the-art advances in metabolic engineering and synthetic biology, here we advocate protocols for overproducing 3-HP in K. pneumoniae. These protocols range from cofactor regeneration, alleviation of metabolite toxicity, genome editing, remodeling of transport system, to carbon flux partition via logic gate. The feasibility for these protocols was also discussed.

Reconstruction and analysis of a genome-scale metabolic model of the oleaginous fungus Mortierella alpina

Background

Mortierella alpina is an oleaginous fungus used in the industrial scale production of arachidonic acid (ARA). In order to investigate the metabolic characteristics at a systems level and to explore potential strategies for enhanced lipid production, a genome-scale metabolic model of M. alpina was reconstructed.

Results

This model included 1106 genes, 1854 reactions and 1732 metabolites. On minimal growth medium, 86 genes were identified as essential, whereas 49 essential genes were identified on yeast extract medium. A series of sequential desaturase and elongase catalysed steps are involved in the synthesis of polyunsaturated fatty acids (PUFAs) from acetyl-CoA precursors, with concomitant NADPH consumption, and these steps were investigated in this study. Oxygen is known to affect the degree of unsaturation of PUFAs, and robustness analysis determined that an oxygen uptake rate of 2.0 mmol gDW⁻¹ h⁻¹ was optimal for ARA accumulation. The flux of 53 reactions involving NADPH was significantly altered at different ARA levels. Of these, malic enzyme (ME) was confirmed as a key component in ARA production and NADPH generation. When using minimization of metabolic adjustment, a knock-out of ME led to a 38.28% decrease in ARA production.

Conclusions

The simulation results confirmed the model as a useful tool for future research on the metabolism of PUFAs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-014-0137-8) contains supplementary material, which is available to authorized users.

Genome-scale modeling for metabolic engineering

We focus on the application of constraint-based methodologies and, more specifically, flux balance analysis in the field of metabolic engineering, and enumerate recent developments and successes of the field. We also review computational frameworks that have been developed with the express purpose of automatically selecting optimal gene deletions for achieving improved production of a chemical of interest. The application of flux balance analysis methods in rational metabolic engineering requires a metabolic network reconstruction and a corresponding in silico metabolic model for the microorganism in question. For this reason, we additionally present a brief overview of automated reconstruction techniques. Finally, we emphasize the importance of integrating metabolic networks with regulatory information-an area which we expect will become increasingly important for metabolic engineering-and present recent developments in the field of metabolic and regulatory integration.

Computational methods in metabolic engineering for strain design

Metabolic engineering uses genetic approaches to control microbial metabolism to produce desired compounds. Computational tools can identify new biological routes to chemicals and the changes needed in host metabolism to improve chemical production. Recent computational efforts have focused on exploring what compounds can be made biologically using native, heterologous, and/or enzymes with broad specificity. Additionally, computational methods have been developed to suggest different types of genetic modifications (e.g. gene deletion/addition or up/down regulation), as well as suggest strategies meeting different criteria (e.g. high yield, high productivity, or substrate co-utilization). Strategies to improve the runtime performances have also been developed, which allow for more complex metabolic engineering strategies to be identified. Future incorporation of kinetic considerations will further improve strain design algorithms.

Succinate Overproduction: A Case Study of Computational Strain Design Using a Comprehensive Escherichia coli Kinetic Model

Computational strain-design prediction accuracy has been the focus for many recent efforts through the selective integration of kinetic information into metabolic models. In general, kinetic model prediction quality is determined by the range and scope of genetic and/or environmental perturbations used during parameterization. In this effort, we apply the k-OptForce procedure on a kinetic model of E. coli core metabolism constructed using the Ensemble Modeling (EM) method and parameterized using multiple mutant strains data under aerobic respiration with glucose as the carbon source. Minimal interventions are identified that improve succinate yield under both aerobic and anaerobic conditions to test the fidelity of model predictions under both genetic and environmental perturbations. Under aerobic condition, k-OptForce identifies interventions that match existing experimental strategies while pointing at a number of unexplored flux re-directions such as routing glyoxylate flux through the glycerate metabolism to improve succinate yield. Many of the identified interventions rely on the kinetic descriptions that would not be discoverable by a purely stoichiometric description. In contrast, under fermentative (anaerobic) condition, k-OptForce fails to identify key interventions including up-regulation of anaplerotic reactions and elimination of competitive fermentative products. This is due to the fact that the pathways activated under anaerobic condition were not properly parameterized as only aerobic flux data were used in the model construction. This study shed light on the importance of condition-specific model parameterization and provides insight on how to augment kinetic models so as to correctly respond to multiple environmental perturbations.

Computational prediction of essential metabolic genes using constraint-based approaches

In this chapter, we describe the application of constraint-based modeling to predict the impact of gene deletions on a metabolic phenotype. The metabolic reactions taking place inside cells form large networks, which have been reconstructed at a genome-scale for several organisms at increasing levels of detail. By integrating mathematical modeling techniques with biochemical principles, constraint-based approaches enable predictions of metabolite fluxes and growth under specific environmental conditions or for genetically modified microorganisms. Similar to the experimental knockout of a gene, predicting the essentiality of a metabolic gene for a phenotype further allows to generate hypotheses on its biological function and design of genetic engineering strategies for biotechnological applications. Here, we summarize the principles of constraint-based approaches and provide a detailed description of the procedure to predict the essentiality of metabolic genes with respect to a specific metabolic function. We exemplify the approach by predicting the essentiality of reactions in the citric acid cycle for the production of glucose from fatty acids.

Cyanobacteria as photosynthetic biocatalysts: a systems biology perspective

A review of cyanobacterial biocatalysts highlighting their metabolic features that argues for the need for systems-level metabolic engineering.

Expanding Metabolic Engineering Algorithms Using Feasible Space and Shadow Price Constraint Modules

While numerous computational methods have been developed that use genome-scale models to propose mutants for the purpose of metabolic engineering, they generally compare mutants based on a single criteria (e.g., production rate at a mutant׳s maximum growth rate). As such, these approaches remain limited in their ability to include multiple complex engineering constraints. To address this shortcoming, we have developed feasible space and shadow price constraint (FaceCon and ShadowCon) modules that can be added to existing mixed integer linear adaptive evolution metabolic engineering algorithms, such as OptKnock and OptORF. These modules allow strain designs to be identified amongst a set of multiple metabolic engineering algorithm solutions that are capable of high chemical production while also satisfying additional design criteria. We describe the various module implementations and their potential applications to the field of metabolic engineering. We then incorporated these modules into the OptORF metabolic engineering algorithm. Using an Escherichia coli genome-scale model (iJO1366), we generated different strain designs for the anaerobic production of ethanol from glucose, thus demonstrating the tractability and potential utility of these modules in metabolic engineering algorithms.

•

Added modules to impose multiple design criteria for engineering algorithms.

•

Examples are provided to eliminate by-product secretion.

•

Examples are provided to control coupling between product and biomass formation.

•

Modules are tractable for genome-scale design.

CONSTRICTOR: constraint modification provides insight into design of biochemical networks

Advances in computational methods that allow for exploration of the combinatorial mutation space are needed to realize the potential of synthetic biology based strain engineering efforts. Here, we present Constrictor, a computational framework that uses flux balance analysis (FBA) to analyze inhibitory effects of genetic mutations on the performance of biochemical networks. Constrictor identifies engineering interventions by classifying the reactions in the metabolic model depending on the extent to which their flux must be decreased to achieve the overproduction target. The optimal inhibition of various reaction pathways is determined by restricting the flux through targeted reactions below the steady state levels of a baseline strain. Constrictor generates unique in silico strains, each representing an “expression state”, or a combination of gene expression levels required to achieve the overproduction target. The Constrictor framework is demonstrated by studying overproduction of ethylene in Escherichia coli network models iAF1260 and iJO1366 through the addition of the heterologous ethylene-forming enzyme from Pseudomonas syringae. Targeting individual reactions as well as combinations of reactions reveals in silico mutants that are predicted to have as high as 25% greater theoretical ethylene yields than the baseline strain during simulated exponential growth. Altering the degree of restriction reveals a large distribution of ethylene yields, while analysis of the expression states that return lower yields provides insight into system bottlenecks. Finally, we demonstrate the ability of Constrictor to scan networks and provide targets for a range of possible products. Constrictor is an adaptable technique that can be used to generate and analyze disparate populations of in silico mutants, select gene expression levels and provide non-intuitive strategies for metabolic engineering.

Computational approaches for microalgal biofuel optimization: a review

The increased demand and consumption of fossil fuels have raised interest in finding renewable energy sources throughout the globe. Much focus has been placed on optimizing microorganisms and primarily microalgae, to efficiently produce compounds that can substitute for fossil fuels. However, the path to achieving economic feasibility is likely to require strain optimization through using available tools and technologies in the fields of systems and synthetic biology. Such approaches invoke a deep understanding of the metabolic networks of the organisms and their genomic and proteomic profiles. The advent of next generation sequencing and other high throughput methods has led to a major increase in availability of biological data. Integration of such disparate data can help define the emergent metabolic system properties, which is of crucial importance in addressing biofuel production optimization. Herein, we review major computational tools and approaches developed and used in order to potentially identify target genes, pathways, and reactions of particular interest to biofuel production in algae. As the use of these tools and approaches has not been fully implemented in algal biofuel research, the aim of this review is to highlight the potential utility of these resources toward their future implementation in algal research.

DRUM: a new framework for metabolic modeling under non-balanced growth. Application to the carbon metabolism of unicellular microalgae

Metabolic modeling is a powerful tool to understand, predict and optimize bioprocesses, particularly when they imply intracellular molecules of interest. Unfortunately, the use of metabolic models for time varying metabolic fluxes is hampered by the lack of experimental data required to define and calibrate the kinetic reaction rates of the metabolic pathways. For this reason, metabolic models are often used under the balanced growth hypothesis. However, for some processes such as the photoautotrophic metabolism of microalgae, the balanced-growth assumption appears to be unreasonable because of the synchronization of their circadian cycle on the daily light. Yet, understanding microalgae metabolism is necessary to optimize the production yield of bioprocesses based on this microorganism, as for example production of third-generation biofuels. In this paper, we propose DRUM, a new dynamic metabolic modeling framework that handles the non-balanced growth condition and hence accumulation of intracellular metabolites. The first stage of the approach consists in splitting the metabolic network into sub-networks describing reactions which are spatially close, and which are assumed to satisfy balanced growth condition. The left metabolites interconnecting the sub-networks behave dynamically. Then, thanks to Elementary Flux Mode analysis, each sub-network is reduced to macroscopic reactions, for which simple kinetics are assumed. Finally, an Ordinary Differential Equation system is obtained to describe substrate consumption, biomass production, products excretion and accumulation of some internal metabolites. DRUM was applied to the accumulation of lipids and carbohydrates of the microalgae Tisochrysis lutea under day/night cycles. The resulting model describes accurately experimental data obtained in day/night conditions. It efficiently predicts the accumulation and consumption of lipids and carbohydrates.

Synthesis and accumulation of aromatic aldehydes in an engineered strain of Escherichia coli

Aromatic aldehydes are useful in numerous applications, especially as flavors, fragrances, and pharmaceutical precursors. However, microbial synthesis of aldehydes is hindered by rapid, endogenous, and redundant conversion of aldehydes to their corresponding alcohols. We report the construction of an Escherichia coli K-12 MG1655 strain with reduced aromatic aldehyde reduction (RARE) that serves as a platform for aromatic aldehyde biosynthesis. Six genes with reported activity on the model substrate benzaldehyde were rationally targeted for deletion: three genes that encode aldo-keto reductases and three genes that encode alcohol dehydrogenases. Upon expression of a recombinant carboxylic acid reductase in the RARE strain and addition of benzoate during growth, benzaldehyde remained in the culture after 24 h, with less than 12% conversion of benzaldehyde to benzyl alcohol. Although individual overexpression results demonstrated that all six genes could contribute to benzaldehyde reduction in vivo, additional experiments featuring subset deletion strains revealed that two of the gene deletions were dispensable under the conditions tested. The engineered strain was next investigated for the production of vanillin from vanillate and succeeded in preventing formation of the byproduct vanillyl alcohol. A pathway for the biosynthesis of vanillin directly from glucose was introduced and resulted in a 55-fold improvement in vanillin titer when using the RARE strain versus the wild-type strain. Finally, synthesis of the chiral pharmaceutical intermediate L-phenylacetylcarbinol (L-PAC) was demonstrated from benzaldehyde and glucose upon expression of a recombinant mutant pyruvate decarboxylase in the RARE strain. Beyond allowing accumulation of aromatic aldehydes as end products in E. coli, the RARE strain expands the classes of chemicals that can be produced microbially via aldehyde intermediates.

Computational approaches to metabolic engineering utilizing systems biology and synthetic biology

Metabolic engineering modifies cellular function to address various biochemical applications. Underlying metabolic engineering efforts are a host of tools and knowledge that are integrated to enable successful outcomes. Concurrent development of computational and experimental tools has enabled different approaches to metabolic engineering. One approach is to leverage knowledge and computational tools to prospectively predict designs to achieve the desired outcome. An alternative approach is to utilize combinatorial experimental tools to empirically explore the range of cellular function and to screen for desired traits. This mini-review focuses on computational systems biology and synthetic biology tools that can be used in combination for prospective in silico strain design.

A computational study of the Warburg effect identifies metabolic targets inhibiting cancer migration

Over the last decade, the field of cancer metabolism has mainly focused on studying the role of tumorigenic metabolic rewiring in supporting cancer proliferation. Here, we perform the first genome-scale computational study of the metabolic underpinnings of cancer migration. We build genome-scale metabolic models of the NCI-60 cell lines that capture the Warburg effect (aerobic glycolysis) typically occurring in cancer cells. The extent of the Warburg effect in each of these cell line models is quantified by the ratio of glycolytic to oxidative ATP flux (AFR), which is found to be highly positively associated with cancer cell migration. We hence predicted that targeting genes that mitigate the Warburg effect by reducing the AFR may specifically inhibit cancer migration. By testing the anti-migratory effects of silencing such 17 top predicted genes in four breast and lung cancer cell lines, we find that up to 13 of these novel predictions significantly attenuate cell migration either in all or one cell line only, while having almost no effect on cell proliferation. Furthermore, in accordance with the predictions, a significant reduction is observed in the ratio between experimentally measured ECAR and OCR levels following these perturbations. Inhibiting anti-migratory targets is a promising future avenue in treating cancer since it may decrease cytotoxic-related side effects that plague current anti-proliferative treatments. Furthermore, it may reduce cytotoxic-related clonal selection of more aggressive cancer cells and the likelihood of emerging resistance.

Integration of bioinformatics to biodegradation

Bioinformatics and biodegradation are two primary scientific fields in applied microbiology and biotechnology. The present review describes development of various bioinformatics tools that may be applied in the field of biodegradation. Several databases, including the University of Minnesota Biocatalysis/Biodegradation database (UM-BBD), a database of biodegradative oxygenases (OxDBase), Biodegradation Network-Molecular Biology Database (Bionemo) MetaCyc, and BioCyc have been developed to enable access to information related to biochemistry and genetics of microbial degradation. In addition, several bioinformatics tools for predicting toxicity and biodegradation of chemicals have been developed. Furthermore, the whole genomes of several potential degrading bacteria have been sequenced and annotated using bioinformatics tools.

Elementary Flux Mode Analysis of Acetyl-CoA Pathway in Carboxydothermus hydrogenoformans Z-2901

Carboxydothermus hydrogenoformans is a carboxydotrophic hydrogenogenic bacterium species that produces hydrogen molecule by utilizing carbon monoxide (CO) or pyruvate as a carbon source. To investigate the underlying biochemical mechanism of hydrogen production, an elementary mode analysis of acetyl-CoA pathway was performed to determine the intermediate fluxes by combining linear programming (LP) method available in CellNetAnalyzer software. We hypothesized that addition of enzymes necessary for carbon monoxide fixation and pyruvate dissimilation would enhance the theoretical yield of hydrogen. An in silico gene knockout of pyk, pykC, and mdh genes of modeled acetyl-CoA pathway allows the maximum theoretical hydrogen yield of 47.62 mmol/gCDW/h for 1 mole of carbon monoxide (CO) uptake. The obtained hydrogen yield is comparatively two times greater than the previous experimental data. Therefore, it could be concluded that this elementary flux mode analysis is a crucial way to achieve efficient hydrogen production through acetyl-CoA pathway and act as a model for strain improvement.