Constraints-based genome-scale metabolic simulation for systems metabolic engineering

Random mutagenesis and selection approaches used traditionally for the development of industrial strains have largely been complemented by metabolic engineering, which allows purposeful modification of metabolic and cellular characteristics by using recombinant DNA and other molecular biological techniques. As systems biology advances as a new paradigm of research thanks to the development of genome-scale computational tools and high-throughput experimental technologies including omics, systems metabolic engineering allowing modification of metabolic, regulatory and signaling networks of the cell at the systems-level is becoming possible. In silico genome-scale metabolic model and its simulation play increasingly important role in providing systematic strategies for metabolic engineering. The in silico genome-scale metabolic model is developed using genomic annotation, metabolic reactions, literature information, and experimental data. The advent of in silico genome-scale metabolic model brought about the development of various algorithms to simulate the metabolic status of the cell as a whole. In this paper, we review the algorithms developed for the system-wide simulation and perturbation of cellular metabolism, discuss the characteristics of these algorithms, and suggest future research direction.     

Systems Biology for Biotech Applications

Cover illustration: Systems Biology for Biotech Applications is the topic of this BTJ special issue edited by Nathan D. Price and Sang Yup Lee. You will find the latest progress on genome-scale in silico models of metabolism, building the blueprint of life and applications for biofuels and bioprocessing. Cover image: degradation, PHA biosynthesis, central metabolic networks (back) and genome-scale metabolic networks (front) in Pseudomonas putida. Courtesy of Tae Yong Kim, KAIST, Korea.     

Development and experimental verification of a genome-scale metabolic model for Corynebacterium glutamicum

Background  

In silico genome-scale metabolic models enable the analysis of the characteristics of metabolic systems of organisms. In this study, we reconstructed a genome-scale metabolic model of Corynebacterium glutamicum on the basis of genome sequence annotation and physiological data. The metabolic characteristics were analyzed using flux balance analysis (FBA), and the results of FBA were validated using data from culture experiments performed at different oxygen uptake rates.     

Hybrid metabolic flux analysis: combining stoichiometric and statistical constraints to model the formation of complex recombinant products

Background  

Stoichiometric models constitute the basic framework for fluxome quantification in the realm of metabolic engineering. A recurrent bottleneck, however, is the establishment of consistent stoichiometric models for the synthesis of recombinant proteins or viruses. Although optimization algorithms for in silico metabolic redesign have been developed in the context of genome-scale stoichiometric models for small molecule production, still rudimentary knowledge of how different cellular levels are regulated and phenotypically expressed prevents their full applicability for complex product optimization.     

Genome-scale reconstruction and <Emphasis Type="Italic">in silico</Emphasis> analysis of the <Emphasis Type="Italic">Clostridium acetobutylicum</Emphasis> ATCC 824 metabolic network

To understand the metabolic characteristics of Clostridium acetobutylicum and to examine the potential for enhanced butanol production, we reconstructed the genome-scale metabolic network from its annotated genomic sequence and analyzed strategies to improve its butanol production. The generated reconstructed network consists of 502 reactions and 479 metabolites and was used as the basis for an in silico model that could compute metabolic and growth performance for comparison with fermentation data. The in silico model successfully predicted metabolic fluxes during the acidogenic phase using classical flux balance analysis. Nonlinear programming was used to predict metabolic fluxes during the solventogenic phase. In addition, essential genes were predicted via single gene deletion studies. This genome-scale in silico metabolic model of C. acetobutylicum should be useful for genome-wide metabolic analysis as well as strain development for improving production of biochemicals, including butanol. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. J. L. and H. Y. equally contributed to this work.     

Mathematical models of cell factories: moving towards the core of industrial biotechnology

Industrial biotechnology involves the utilization of cell factories for the production of fuels and chemicals. Traditionally, the development of highly productive microbial strains has relied on random mutagenesis and screening. The development of predictive mathematical models provides a new paradigm for the rational design of cell factories. Instead of selecting among a set of strains resulting from random mutagenesis, mathematical models allow the researchers to predict in silico the outcomes of different genetic manipulations and engineer new strains by performing gene deletions or additions leading to a higher productivity of the desired chemicals. In this review we aim to summarize the main modelling approaches of biological processes and illustrate the particular applications that they have found in the field of industrial microbiology.     

In silico analysis of bioethanol production from glucose/xylose mixtures during fed-batch fermentation of co-culture and mono-culture systems

This study presents a detailed in silico analysis of bioethanol production from glucose/xylose mixtures of various compositions by fed-batch co-culture and mono-culture fermentation of specialized microbes. The mono-culture consists of recombinant Saccharomyces cerevisise that can metabolize both hexose and pentose sugars while the co-culture system consists of substrate-selective microbes. Dynamic flux balance models based on available genome-scale reconstructions of the microorganisms have been used to analyze bioethanol production in fed-batch culture with constant feed rates and the maximization of ethanol productivity is addressed by computing optimal aerobic-anaerobic switching times. The simulation results clearly point to the superior performance of fed-batch fermentation of microbial co-culture against fed-batch fermentation of mono-culture for bioethanol production from glucose/xylose mixtures. A set of potential genetic engineering strategies for enhancement of S. cerevisiae and Escherichia coli strains performance have been identified. Such in silico predictions using genome-scale models provide valuable guidance for conducting in vivo metabolic engineering experiments.     

Genome-scale metabolic models as platforms for strain design and biological discovery

Genome-scale metabolic models (GEMs) have been developed and used in guiding systems’ metabolic engineering strategies for strain design and development. This strategy has been used in fermentative production of bio-based industrial chemicals and fuels from alternative carbon sources. However, computer-aided hypotheses building using established algorithms and software platforms for biological discovery can be integrated into the pipeline for strain design strategy to create superior strains of microorganisms for targeted biosynthetic goals. Here, I described an integrated workflow strategy using GEMs for strain design and biological discovery. Specific case studies of strain design and biological discovery using Escherichia coli genome-scale model are presented and discussed. The integrated workflow presented herein, when applied carefully would help guide future design strategies for high-performance microbial strains that have existing and forthcoming genome-scale metabolic models.     

The genome-scale metabolic network analysis <Emphasis Type="Italic">of Zymomonas mobilis</Emphasis> ZM4 explains physiological features and suggests ethanol and succinic acid production strategies

Background  

Zymomonas mobilis ZM4 is a Gram-negative bacterium that can efficiently produce ethanol from various carbon substrates, including glucose, fructose, and sucrose, via the Entner-Doudoroff pathway. However, systems metabolic engineering is required to further enhance its metabolic performance for industrial application. As an important step towards this goal, the genome-scale metabolic model of Z. mobilis is required to systematically analyze in silico the metabolic characteristics of this bacterium under a wide range of genotypic and environmental conditions.     

Systems-level analysis of genome-scalein silico metabolic models using MetaFluxNet

The systems-level analysis of microbes with myriad of heterologous data generated by omics technologies has been applied to improve our understanding of cellular function and physiology and consequently to enhance production of various bioproducts. At the heart of this revolution residesin silico genome-scale metabolic model. In order to fully exploit the power of genome-scale model, a systematic approach employing user-friendly software is required. Metabolic flux analysis of genome-scale metabolic network is becoming widely employed to quantify the flux distribution and validate model-driven hypotheses. Here we describe the development of an upgraded MetaFluxNet which allows (1) construction of metabolic models connected to metabolic databases, (2) calculation of fluxes by metabolic flux analysis, (3) comparative flux analysis with flux-profile visualization, (4) the use of metabolic flux analysis markup language to enable models to be exchanged efficiently, and (5) the exporting of data from constraints-based flux analysis into various formats. MetaFluxNet also allows cellular physiology to be predicted and strategies for strain improvement to be developed from genome-based information on flux distributions. This integrated software environment promises to enhance our understanding on metabolic network at a whole organism level and to establish novel strategies for improving the properties of organisms for various biotechnological applications.     

In silico genome-scale metabolic analysis of Pseudomonas putida KT2440 for polyhydroxyalkanoate synthesis,degradation of aromatics and anaerobic survival

Genome-scale metabolic models have been appearing with increasing frequency and have been employed in a wide range of biotechnological applications as well as in biological studies. With the metabolic model as a platform, engineering strategies have become more systematic and focused, unlike the random shotgun approach used in the past. Here we present the genome-scale metabolic model of the versatile Gram-negative bacterium Pseudomonas putida, which has gained widespread interest for various biotechnological applications. With the construction of the genome-scale metabolic model of P. putida KT2440, PpuMBEL1071, we investigated various characteristics of P. putida, such as its capacity for synthesizing polyhydroxyalkanoates (PHA) and degrading aromatics. Although P. putida has been characterized as a strict aerobic bacterium, the physiological characteristics required to achieve anaerobic survival were investigated. Through analysis of PpuMBEL1071, extended survival of P. putida under anaerobic stress was achieved by introducing the ackA gene from Pseudomonas aeruginosa and Escherichia coli.     

Systems metabolic engineering: Genome-scale models and beyond

The advent of high throughput genome-scale bioinformatics has led to an exponential increase in available cellular system data. Systems metabolic engineering attempts to use data-driven approaches – based on the data collected with high throughput technologies – to identify gene targets and optimize phenotypical properties on a systems level. Current systems metabolic engineering tools are limited for predicting and defining complex phenotypes such as chemical tolerances and other global, multigenic traits. The most pragmatic systems-based tool for metabolic engineering to arise is the in silico genome-scale metabolic reconstruction. This tool has seen wide adoption for modeling cell growth and predicting beneficial gene knockouts, and we examine here how this approach can be expanded for novel organisms. This review will highlight advances of the systems metabolic engineering approach with a focus on de novo development and use of genome-scale metabolic reconstructions for metabolic engineering applications. We will then discuss the challenges and prospects for this emerging field to enable model-based metabolic engineering. Specifically, we argue that current state-of-the-art systems metabolic engineering techniques represent a viable first step for improving product yield that still must be followed by combinatorial techniques or random strain mutagenesis to achieve optimal cellular systems.     

Flux analysis and metabolomics for systematic metabolic engineering of microorganisms

Rational engineering of metabolism is important for bio-production using microorganisms. Metabolic design based on in silico simulations and experimental validation of the metabolic state in the engineered strain helps in accomplishing systematic metabolic engineering. Flux balance analysis (FBA) is a method for the prediction of metabolic phenotype, and many applications have been developed using FBA to design metabolic networks. Elementary mode analysis (EMA) and ensemble modeling techniques are also useful tools for in silico strain design. The metabolome and flux distribution of the metabolic pathways enable us to evaluate the metabolic state and provide useful clues to improve target productivity. Here, we reviewed several computational applications for metabolic engineering by using genome-scale metabolic models of microorganisms. We also discussed the recent progress made in the field of metabolomics and ¹³C-metabolic flux analysis techniques, and reviewed these applications pertaining to bio-production development. Because these in silico or experimental approaches have their respective advantages and disadvantages, the combined usage of these methods is complementary and effective for metabolic engineering.     

Simulation-based model checking approach to cell fate specification during Caenorhabditis elegans vulval development by hybrid functional Petri net with extension

Background  

Model checking approaches were applied to biological pathway validations around 2003. Recently, Fisher et al. have proved the importance of model checking approach by inferring new regulation of signaling crosstalk in C. elegans and confirming the regulation with biological experiments. They took a discrete and state-based approach to explore all possible states of the system underlying vulval precursor cell (VPC) fate specification for desired properties. However, since both discrete and continuous features appear to be an indispensable part of biological processes, it is more appropriate to use quantitative models to capture the dynamics of biological systems. Our key motivation of this paper is to establish a quantitative methodology to model and analyze in silico models incorporating the use of model checking approach.     

Evaluation of a Genome-Scale In Silico Metabolic Model for Geobacter metallireducens by Using Proteomic Data from a Field Biostimulation Experiment

Accurately predicting the interactions between microbial metabolism and the physical subsurface environment is necessary to enhance subsurface energy development, soil and groundwater cleanup, and carbon management. This study was an initial attempt to confirm the metabolic functional roles within an in silico model using environmental proteomic data collected during field experiments. Shotgun global proteomics data collected during a subsurface biostimulation experiment were used to validate a genome-scale metabolic model of Geobacter metallireducens—specifically, the ability of the metabolic model to predict metal reduction, biomass yield, and growth rate under dynamic field conditions. The constraint-based in silico model of G. metallireducens relates an annotated genome sequence to the physiological functions with 697 reactions controlled by 747 enzyme-coding genes. Proteomic analysis showed that 180 of the 637 G. metallireducens proteins detected during the 2008 experiment were associated with specific metabolic reactions in the in silico model. When the field-calibrated Fe(III) terminal electron acceptor process reaction in a reactive transport model for the field experiments was replaced with the genome-scale model, the model predicted that the largest metabolic fluxes through the in silico model reactions generally correspond to the highest abundances of proteins that catalyze those reactions. Central metabolism predicted by the model agrees well with protein abundance profiles inferred from proteomic analysis. Model discrepancies with the proteomic data, such as the relatively low abundances of proteins associated with amino acid transport and metabolism, revealed pathways or flux constraints in the in silico model that could be updated to more accurately predict metabolic processes that occur in the subsurface environment.     

What CHO is made of: Variations in the biomass composition of Chinese hamster ovary cell lines

BackgroundCell line-specific, genome-scale metabolic models enable rigorous and systematic in silico investigation of cellular metabolism. Such models have recently become available for Chinese hamster ovary (CHO) cells. However, a key ingredient, namely an experimentally validated biomass function that summarizes the cellular composition, was so far missing. Here, we close this gap by providing extensive experimental data on the biomass composition of 13 parental and producer CHO cell lines under various conditions.ResultsWe report total protein, lipid, DNA, RNA and carbohydrate content, cell dry mass, and detailed protein and lipid composition. Furthermore, we present meticulous data on exchange rates between cells and environment and provide detailed experimental protocols on how to determine all of the above. The biomass composition is converted into cell line- and condition-specific biomass functions for use in cell line-specific, genome-scale metabolic models of CHO. Finally, flux balance analysis (FBA) is used to demonstrate consistency between in silico predictions and experimental analysis.ConclusionsOur study reveals a strong variability of the total protein content and cell dry mass across cell lines. However, the relative amino acid composition is independent of the cell line and condition and thus needs not be explicitly measured for each new cell line. In contrast, the lipid composition is strongly influenced by the growth media and thus will have to be determined in each case. These cell line-specific variations in biomass composition have a small impact on growth rate predictions with FBA, as inaccuracies in the predictions are rather dominated by inaccuracies in the exchange rate spectra. Cell-specific biomass variations only become important if the experimental errors in the exchange rate spectra drop below twenty percent.     

Genome-scale in silico aided metabolic analysis and flux comparisons of Escherichia coli to improve succinate production

In the post-genome era, it is one challenge to understand the cellular metabolism at the systematic levels. Mathematical modeling of microorganisms and subsequent computer simulation are effective tools for systems biology. In this paper, based on the genome-scale Escherichia coli stoichiometric model iJR904, through the GAMS linear programming package, the in silico maximal succinate yield was estimated to be 1.714 mol/mol glucose. When another two constraints were added, the maximal succinate yield dropped to 1.60 mol/mol glucose. Further analysis substantiated the uniqueness of the flux distribution under such constraints. After comparisons with the metabolic flux analysis (MFA) results computed from the wet experimental data of the three kinds of E. coli, three potential improvement target sites, the glucose phosphotransferase transport system, the pyruvate carboxylase, and the glyoxylate shunt, were identified and selected for the genetic modifications. All the three genetic modified strains showed increased succinate yield. The final strain TUQ19/pQZ6 had a high yield of 1.29 mol succinate/mol glucose and high productivity. The success of the above experiments proved that this in silico optimal succinate production pathway is reasonable and practical. This method may also be used as a general strategy to help enhance the yields of other favorable metabolites in E. coli.     

Nuclear staining and relative distance for quantifying epidermal differentiation in biomarker expression profiling

Background  

The epidermal physiology results from a complex regulated homeostasis of keratinocyte proliferation, differentiation and death and is tightly regulated by a specific protein expression during cellular maturation. Cellular in silico models are considered a promising and inevitable tool for the understanding of this complex system. Hence, we need to incorporate the information of the differentiation dependent protein expression in cell based systems biological models of tissue homeostasis. Such methods require measuring tissue differentiation quantitatively while correlating it with biomarker expression intensities.