Dynamic flux balance analysis of diauxic growth in Escherichia coli

Genome-based Modeling and Design of Metabolic Interactions in Microbial Communities

Biotechnology research is traditionally focused on individual microbial strains that are perceived to have the necessary metabolic functions, or the capability to have these functions introduced, to achieve a particular task. For many important applications, the development of such omnipotent microbes is an extremely challenging if not impossible task. By contrast, nature employs a radically different strategy based on synergistic combinations of different microbial species that collectively achieve the desired task. These natural communities have evolved to exploit the native metabolic capabilities of each species and are highly adaptive to changes in their environments. However, microbial communities have proven difficult to study due to a lack of suitable experimental and computational tools. With the advent of genome sequencing, omics technologies, bioinformatics and genome-scale modeling, researchers now have unprecedented capabilities to analyze and engineer the metabolism of microbial communities. The goal of this review is to summarize recent applications of genome-scale metabolic modeling to microbial communities. A brief introduction to lumped community models is used to motivate the need for genome-level descriptions of individual species and their metabolic interactions. The review of genome-scale models begins with static modeling approaches, which are appropriate for communities where the extracellular environment can be assumed to be time invariant or slowly varying. Dynamic extensions of the static modeling approach are described, and then applications of genome-scale models for design of synthetic microbial communities are reviewed. The review concludes with a summary of metagenomic tools for analyzing community metabolism and an outlook for future research.

Metabolic modelling in the development of cell factories by synthetic biology

Cell factories are commonly microbial organisms utilized for bioconversion of renewable resources to bulk or high value chemicals. Introduction of novel production pathways in chassis strains is the core of the development of cell factories by synthetic biology. Synthetic biology aims to create novel biological functions and systems not found in nature by combining biology with engineering. The workflow of the development of novel cell factories with synthetic biology is ideally linear which will be attainable with the quantitative engineering approach, high-quality predictive models, and libraries of well-characterized parts. Different types of metabolic models, mathematical representations of metabolism and its components, enzymes and metabolites, are useful in particular phases of the synthetic biology workflow. In this minireview, the role of metabolic modelling in synthetic biology will be discussed with a review of current status of compatible methods and models for the in silico design and quantitative evaluation of a cell factory.

Software applications toward quantitative metabolic flux analysis and modeling

Metabolites and their pathways are central for adaptation and survival. Metabolic modeling elucidates in silico all the possible flux pathways (flux balance analysis, FBA) and predicts the actual fluxes under a given situation, further refinement of these models is possible by including experimental isotopologue data. In this review, we initially introduce the key theoretical concepts and different analysis steps in the modeling process before comparing flux calculation and metabolite analysis programs such as C13, BioOpt, COBRA toolbox, Metatool, efmtool, FiatFlux, ReMatch, VANTED, iMAT and YANA. Their respective strengths and limitations are discussed and compared to alternative software. While data analysis of metabolites, calculation of metabolic fluxes, pathways and their condition-specific changes are all possible, we highlight the considerations that need to be taken into account before deciding on a specific software. Current challenges in the field include the computation of large-scale networks (in elementary mode analysis), regulatory interactions and detailed kinetics, and these are discussed in the light of powerful new approaches.

Optimization of auto-induction medium for G-CSF production by Escherichia coli using artificial neural networks coupled with genetic algorithm

Granulocyte colony-stimulating factor (G-CSF) is a cytokine widely used in cancer patients receiving high doses of chemotherapeutic drugs to prevent the chemotherapy-induced suppression of white blood cells. The production of recombinant G-CSF should be increased to meet the increasing market demand. This study aims to model and optimize the carbon source of auto-induction medium to enhance G-CSF production using artificial neural networks coupled with genetic algorithm. In this approach, artificial neural networks served as bioprocess modeling tools, and genetic algorithm (GA) was applied to optimize the established artificial neural network models. Two artificial neural network models were constructed: the back-propagation (BP) network and the radial basis function (RBF) network. The root mean square error, coefficient of determination, and standard error of prediction of the BP model were 0.0375, 0.959, and 8.49 %, respectively, whereas those of the RBF model were 0.0257, 0.980, and 5.82 %, respectively. These values indicated that the RBF model possessed higher fitness and prediction accuracy than the BP model. Under the optimized auto-induction medium, the predicted maximum G-CSF yield by the BP-GA approach was 71.66 %, whereas that by the RBF-GA approach was 75.17 %. These predicted values are in agreement with the experimental results, with 72.4 and 76.014 % for the BP-GA and RBF-GA models, respectively. These results suggest that RBF-GA is superior to BP-GA. The developed approach in this study may be helpful in modeling and optimizing other multivariable, non-linear, and time-variant bioprocesses.

Reconstruction and comparison of the metabolic potential of cyanobacteria Cyanothece sp. ATCC 51142 and Synechocystis sp. PCC 6803

Cyanobacteria are an important group of photoautotrophic organisms that can synthesize valuable bio-products by harnessing solar energy. They are endowed with high photosynthetic efficiencies and diverse metabolic capabilities that confer the ability to convert solar energy into a variety of biofuels and their precursors. However, less well studied are the similarities and differences in metabolism of different species of cyanobacteria as they pertain to their suitability as microbial production chassis. Here we assemble, update and compare genome-scale models (iCyt773 and iSyn731) for two phylogenetically related cyanobacterial species, namely Cyanothece sp. ATCC 51142 and Synechocystis sp. PCC 6803. All reactions are elementally and charge balanced and localized into four different intracellular compartments (i.e., periplasm, cytosol, carboxysome and thylakoid lumen) and biomass descriptions are derived based on experimental measurements. Newly added reactions absent in earlier models (266 and 322, respectively) span most metabolic pathways with an emphasis on lipid biosynthesis. All thermodynamically infeasible loops are identified and eliminated from both models. Comparisons of model predictions against gene essentiality data reveal a specificity of 0.94 (94/100) and a sensitivity of 1 (19/19) for the Synechocystis iSyn731 model. The diurnal rhythm of Cyanothece 51142 metabolism is modeled by constructing separate (light/dark) biomass equations and introducing regulatory restrictions over light and dark phases. Specific metabolic pathway differences between the two cyanobacteria alluding to different bio-production potentials are reflected in both models.

Integrating cellular metabolism into a multiscale whole-body model

Cellular metabolism continuously processes an enormous range of external compounds into endogenous metabolites and is as such a key element in human physiology. The multifaceted physiological role of the metabolic network fulfilling the catalytic conversions can only be fully understood from a whole-body perspective where the causal interplay of the metabolic states of individual cells, the surrounding tissue and the whole organism are simultaneously considered. We here present an approach relying on dynamic flux balance analysis that allows the integration of metabolic networks at the cellular scale into standardized physiologically-based pharmacokinetic models at the whole-body level. To evaluate our approach we integrated a genome-scale network reconstruction of a human hepatocyte into the liver tissue of a physiologically-based pharmacokinetic model of a human adult. The resulting multiscale model was used to investigate hyperuricemia therapy, ammonia detoxification and paracetamol-induced toxication at a systems level. The specific models simultaneously integrate multiple layers of biological organization and offer mechanistic insights into pathology and medication. The approach presented may in future support a mechanistic understanding in diagnostics and drug development.

Cellular metabolism is a key element in human physiology. Ideally the metabolic network needs to be considered within the context of the surrounding tissue and organism since the various levels of biological organization are mutually influencing each other. To mechanistically describe the interplay between intracellular space and extracellular environment, we here integrate the genome-scale metabolic network model HepatoNet1 at the cellular scale into physiologically-based pharmacokinetic models at the whole-body level. The resulting multiscale model allows the quantitative description of metabolic behavior in the context of time-resolved metabolite concentration profiles in the body and the surrounding liver tissue. The model has been applied to three case studies covering fundamental aspects of medicine and pharmacology: drug administration, biomarker identification and drug-induced toxication. Most notably, our multiscale approach fosters an improved quantitative understanding of drug action and the impact of metabolic disorders at an organism level, based on a genome-scale representation of cellular metabolism. Computational models such as the one presented include various aspects of human physiology and may therefore significantly support rational approaches in medical diagnostics and pharmaceutical drug development in the future.

Engineering ecosystems and synthetic ecologies

Microbial ecosystems play an important role in nature. Engineering these systems for industrial, medical, or biotechnological purposes are important pursuits for synthetic biologists and biological engineers moving forward. Here we provide a review of recent progress in engineering natural and synthetic microbial ecosystems. We highlight important forward engineering design principles, theoretical and quantitative models, new experimental and manipulation tools, and possible applications of microbial ecosystem engineering. We argue that simply engineering individual microbes will lead to fragile homogenous populations that are difficult to sustain, especially in highly heterogeneous and unpredictable environments. Instead, engineered microbial ecosystems are likely to be more robust and able to achieve complex tasks at the spatial and temporal resolution needed for truly programmable biology.

Mathematical optimization applications in metabolic networks

Genome-scale metabolic models are increasingly becoming available for a variety of microorganisms. This has spurred the development of a wide array of computational tools, and in particular, mathematical optimization approaches, to assist in fundamental metabolic network analyses and redesign efforts. This review highlights a number of optimization-based frameworks developed towards addressing challenges in the analysis and engineering of metabolic networks. In particular, three major types of studies are covered here including exploring model predictions, correction and improvement of models of metabolism, and redesign of metabolic networks for the targeted overproduction of a desired compound. Overall, the methods reviewed in this paper highlight the diversity of queries, breadth of questions and complexity of redesigns that are amenable to mathematical optimization strategies.

MicrobesFlux: a web platform for drafting metabolic models from the KEGG database

Background

Concurrent with the efforts currently underway in mapping microbial genomes using high-throughput sequencing methods, systems biologists are building metabolic models to characterize and predict cell metabolisms. One of the key steps in building a metabolic model is using multiple databases to collect and assemble essential information about genome-annotations and the architecture of the metabolic network for a specific organism. To speed up metabolic model development for a large number of microorganisms, we need a user-friendly platform to construct metabolic networks and to perform constraint-based flux balance analysis based on genome databases and experimental results.

Results

We have developed a semi-automatic, web-based platform (MicrobesFlux) for generating and reconstructing metabolic models for annotated microorganisms. MicrobesFlux is able to automatically download the metabolic network (including enzymatic reactions and metabolites) of ~1,200 species from the KEGG database (Kyoto Encyclopedia of Genes and Genomes) and then convert it to a metabolic model draft. The platform also provides diverse customized tools, such as gene knockouts and the introduction of heterologous pathways, for users to reconstruct the model network. The reconstructed metabolic network can be formulated to a constraint-based flux model to predict and analyze the carbon fluxes in microbial metabolisms. The simulation results can be exported in the SBML format (The Systems Biology Markup Language). Furthermore, we also demonstrated the platform functionalities by developing an FBA model (including 229 reactions) for a recent annotated bioethanol producer, Thermoanaerobacter sp. strain X514, to predict its biomass growth and ethanol production.

Conclusion

MicrobesFlux is an installation-free and open-source platform that enables biologists without prior programming knowledge to develop metabolic models for annotated microorganisms in the KEGG database. Our system facilitates users to reconstruct metabolic networks of organisms based on experimental information. Through human-computer interaction, MicrobesFlux provides users with reasonable predictions of microbial metabolism via flux balance analysis. This prototype platform can be a springboard for advanced and broad-scope modeling of complex biological systems by integrating other “omics” data or ¹³ C- metabolic flux analysis results. MicrobesFlux is available at http://tanglab.engineering.wustl.edu/static/MicrobesFlux.html and will be continuously improved based on feedback from users.

Nutritional requirements of the BY series of Saccharomyces cerevisiae strains for optimum growth

Among the vast variety of Saccharomyces cerevisiae strains, the BY family is particularly important because the widely used deletion collections are based on this background. Here we demonstrate that some standard growth media recipes require substantial modifications to provide optimum growth conditions for auxotrophic BY strains and to avoid growth arrest before glucose is depleted. In addition to the essential supplements that are required to satisfy auxotrophic requirements, we found the four amino acids phenylalanine, glutamic acid, serine, and threonine to be indispensable for optimum growth, despite the fact that BY is 'prototrophic' for these amino acids. Interestingly, other widely used S. cerevisiae strains, such as strains of the CEN.PK family, are less sensitive to lack of the described supplements. Furthermore, we found that the concentration of inositol in yeast nitrogen base is too low to support fast proliferation of yeast cultures until glucose is exhausted. Depletion of inositol during exponential growth induces characteristic changes, namely a decrease in glucose uptake and maximum specific growth rate, increased cell size, reduced viability, and accumulation of lipid storage pools. Thus, several of the existing growth media recipes need to be revised to achieve optimum growth conditions for BY-derived strains.

Estimation of dynamic flux profiles from metabolic time series data

Background

Advances in modern high-throughput techniques of molecular biology have enabled top-down approaches for the estimation of parameter values in metabolic systems, based on time series data. Special among them is the recent method of dynamic flux estimation (DFE), which uses such data not only for parameter estimation but also for the identification of functional forms of the processes governing a metabolic system. DFE furthermore provides diagnostic tools for the evaluation of model validity and of the quality of a model fit beyond residual errors. Unfortunately, DFE works only when the data are more or less complete and the system contains as many independent fluxes as metabolites. These drawbacks may be ameliorated with other types of estimation and information. However, such supplementations incur their own limitations. In particular, assumptions must be made regarding the functional forms of some processes and detailed kinetic information must be available, in addition to the time series data.

Results

The authors propose here a systematic approach that supplements DFE and overcomes some of its shortcomings. Like DFE, the approach is model-free and requires only minimal assumptions. If sufficient time series data are available, the approach allows the determination of a subset of fluxes that enables the subsequent applicability of DFE to the rest of the flux system. The authors demonstrate the procedure with three artificial pathway systems exhibiting distinct characteristics and with actual data of the trehalose pathway in Saccharomyces cerevisiae.

Conclusions

The results demonstrate that the proposed method successfully complements DFE under various situations and without a priori assumptions regarding the model representation. The proposed method also permits an examination of whether at all, to what degree, or within what range the available time series data can be validly represented in a particular functional format of a flux within a pathway system. Based on these results, further experiments may be designed to generate data points that genuinely add new information to the structure identification and parameter estimation tasks at hand.

A model of flux regulation in the cholesterol biosynthesis pathway: Immune mediated graduated flux reduction versus statin-like led stepped flux reduction

The cholesterol biosynthesis pathway has recently been shown to play an important role in the innate immune response to viral infection with host protection occurring through a coordinate down regulation of the enzymes catalysing each metabolic step. In contrast, statin based drugs, which form the principle pharmaceutical agents for decreasing the activity of this pathway, target a single enzyme. Here, we build an ordinary differential equation model of the cholesterol biosynthesis pathway in order to investigate how the two regulatory strategies impact upon the behaviour of the pathway. We employ a modest set of assumptions: that the pathway operates away from saturation, that each metabolite is involved in multiple cellular interactions and that mRNA levels reflect enzyme concentrations. Using data taken from primary bone marrow derived macrophage cells infected with murine cytomegalovirus or treated with IFNγ, we show that, under these assumptions, coordinate down-regulation of enzyme activity imparts a graduated reduction in flux along the pathway. In contrast, modelling a statin-like treatment that achieves the same degree of down-regulation in cholesterol production, we show that this delivers a step change in flux along the pathway. The graduated reduction mediated by physiological coordinate regulation of multiple enzymes supports a mechanism that allows a greater level of specificity, altering cholesterol levels with less impact upon interactions branching from the pathway, than pharmacological step reductions. We argue that coordinate regulation is likely to show a long-term evolutionary advantage over single enzyme regulation. Finally, the results from our models have implications for future pharmaceutical therapies intended to target cholesterol production with greater specificity and fewer off target effects, suggesting that this can be achieved by mimicking the coordinated down-regulation observed in immunological responses.

Detection of transcriptional triggers in the dynamics of microbial growth: application to the respiratorily versatile bacterium Shewanella oneidensis

The capacity of microorganisms to respond to variable external conditions requires a coordination of environment-sensing mechanisms and decision-making regulatory circuits. Here, we seek to understand the interplay between these two processes by combining high-throughput measurement of time-dependent mRNA profiles with a novel computational approach that searches for key genetic triggers of transcriptional changes. Our approach helped us understand the regulatory strategies of a respiratorily versatile bacterium with promising bioenergy and bioremediation applications, Shewanella oneidensis, in minimal and rich media. By comparing expression profiles across these two conditions, we unveiled components of the transcriptional program that depend mainly on the growth phase. Conversely, by integrating our time-dependent data with a previously available large compendium of static perturbation responses, we identified transcriptional changes that cannot be explained solely by internal network dynamics, but are rather triggered by specific genes acting as key mediators of an environment-dependent response. These transcriptional triggers include known and novel regulators that respond to carbon, nitrogen and oxygen limitation. Our analysis suggests a sequence of physiological responses, including a coupling between nitrogen depletion and glycogen storage, partially recapitulated through dynamic flux balance analysis, and experimentally confirmed by metabolite measurements. Our approach is broadly applicable to other systems.

Experimental and modeling analysis of Synechocystis sp. PCC 6803 growth

BACKGROUND/AIMS

The influence of different parameters such as temperature, irradiance, nitrate concentration, pH, and an external carbon source on Synechocystis PCC 6803 growth was evaluated.

METHODS

4.5-ml cuvettes containing 2 ml of culture, a high-throughput system equivalent to batch cultures, were used with gas exchange ensured by the use of a Parafilm™ cover. The effect of the different variables on maximum growth was assessed by a multi-way statistical analysis.

RESULTS

Temperature and pH were identified as the key factors. It was observed that Synechocystis cells have a strong influence on the external pH. The optimal growth temperature was 33°C while light-saturating conditions were reached at 40 µE·m⁻²·s⁻¹.

CONCLUSION

It was demonstrated that Synechocystis exhibits a marked difference in behavior between autotrophic and glucose-based mixotrophic conditions, and that nitrate concentrations did not have a significant influence, probably due to endogenous nitrogen reserves. Furthermore, a dynamic metabolic model of Synechocystis photosynthesis was developed to gain insights on the underlying mechanism enabling this cyanobacterium to control the levels of external pH. The model showed a coupled effect between the increase of the pH and ATP production which in turn allows a higher carbon fixation rate.

Dynamic regulatory on/off minimization for biological systems under internal temporal perturbations

Background

Flux balance analysis (FBA) together with its extension, dynamic FBA, have proven instrumental for analyzing the robustness and dynamics of metabolic networks by employing only the stoichiometry of the included reactions coupled with adequately chosen objective function. In addition, under the assumption of minimization of metabolic adjustment, dynamic FBA has recently been employed to analyze the transition between metabolic states.

Results

Here, we propose a suite of novel methods for analyzing the dynamics of (internally perturbed) metabolic networks and for quantifying their robustness with limited knowledge of kinetic parameters. Following the biochemically meaningful premise that metabolite concentrations exhibit smooth temporal changes, the proposed methods rely on minimizing the significant fluctuations of metabolic profiles to predict the time-resolved metabolic state, characterized by both fluxes and concentrations. By conducting a comparative analysis with a kinetic model of the Calvin-Benson cycle and a model of plant carbohydrate metabolism, we demonstrate that the principle of regulatory on/off minimization coupled with dynamic FBA can accurately predict the changes in metabolic states.

Conclusions

Our methods outperform the existing dynamic FBA-based modeling alternatives, and could help in revealing the mechanisms for maintaining robustness of dynamic processes in metabolic networks over time.

Use of an uncertainty analysis for genome-scale models as a prediction tool for microbial growth processes in subsurface environments

The integration of genome-scale, constraint-based models of microbial cell function into simulations of contaminant transport and fate in complex groundwater systems is a promising approach to help characterize the metabolic activities of microorganisms in natural environments. In constraint-based modeling, the specific uptake flux rates of external metabolites are usually determined by Michaelis-Menten kinetic theory. However, extensive data sets based on experimentally measured values are not always available. In this study, a genome-scale model of Pseudomonas putida was used to study the key issue of uncertainty arising from the parametrization of the influx of two growth-limiting substrates: oxygen and toluene. The results showed that simulated growth rates are highly sensitive to substrate affinity constants and that uncertainties in specific substrate uptake rates have a significant influence on the variability of simulated microbial growth. Michaelis-Menten kinetic theory does not, therefore, seem to be appropriate for descriptions of substrate uptake processes in the genome-scale model of P. putida. Microbial growth rates of P. putida in subsurface environments can only be accurately predicted if the processes of complex substrate transport and microbial uptake regulation are sufficiently understood in natural environments and if data-driven uptake flux constraints can be applied.

Modelling cyanobacteria: from metabolism to integrative models of phototrophic growth

Cyanobacteria are phototrophic microorganisms of global importance and have recently attracted increasing attention due to their capability to convert sunlight and atmospheric CO(2) directly into organic compounds, including carbon-based biofuels. The utilization of cyanobacteria as a biological chassis to generate third-generation biofuels would greatly benefit from an increased understanding of cyanobacterial metabolism and its interplay with other cellular processes. In this respect, metabolic modelling has been proposed as a way to overcome the traditional trial and error methodology that is often employed to introduce novel pathways. In particular, flux balance analysis and related methods have proved to be powerful tools to investigate the organization of large-scale metabolic networks-with the prospect of predicting modifications that are likely to increase the yield of a desired product and thereby to streamline the experimental progress and avoid futile avenues. This contribution seeks to describe the utilization of metabolic modelling as a research tool to understand the metabolism and phototrophic growth of cyanobacteria. The focus of the contribution is on a mathematical description of the metabolic network of Synechocystis sp. PCC 6803 and its analysis using constraint-based methods. A particular challenge is to integrate the description of the metabolic network with other cellular processes, such as the circadian clock, the photosynthetic light reactions, carbon concentration mechanism, and transcriptional regulation-aiming at a predictive model of a cyanobacterium in silico.

Model-based confirmation of alternative substrates of mitochondrial electron transport chain

Discrimination of metabolic models based on high throughput metabolomics data, reflecting various internal and external perturbations, is essential for identifying the components that contribute to the emerging behavior of metabolic processes. Here, we investigate 12 different models of the mitochondrial electron transport chain (ETC) in Arabidopsis thaliana during dark-induced senescence in order to elucidate the alternative substrates to this metabolic pathway. Our findings demonstrate that the coupling of the proposed computational approach, based on dynamic flux balance analysis, with time-resolved metabolomics data results in model-based confirmations of the hypotheses that, during dark-induced senescence in Arabidopsis, (i) under conditions where the main substrate for the ETC are not fully available, isovaleryl-CoA dehydrogenase and 2-hydroxyglutarate dehydrogenase are able to donate electrons to the ETC, (ii) phytanoyl-CoA does not act even as an indirect substrate of the electron transfer flavoprotein/electron-transfer flavoprotein:ubiquinone oxidoreductase complex, and (iii) the mitochondrial γ-aminobutyric acid transporter has functional significance in maintaining mitochondrial metabolism. Our study provides a basic framework for future in silico studies of alternative pathways in mitochondrial metabolism under extended darkness whereby the role of its components can be computationally discriminated based on available molecular profile data.

OptCom: a multi-level optimization framework for the metabolic modeling and analysis of microbial communities

Microorganisms rarely live isolated in their natural environments but rather function in consolidated and socializing communities. Despite the growing availability of high-throughput sequencing and metagenomic data, we still know very little about the metabolic contributions of individual microbial players within an ecological niche and the extent and directionality of interactions among them. This calls for development of efficient modeling frameworks to shed light on less understood aspects of metabolism in microbial communities. Here, we introduce OptCom, a comprehensive flux balance analysis framework for microbial communities, which relies on a multi-level and multi-objective optimization formulation to properly describe trade-offs between individual vs. community level fitness criteria. In contrast to earlier approaches that rely on a single objective function, here, we consider species-level fitness criteria for the inner problems while relying on community-level objective maximization for the outer problem. OptCom is general enough to capture any type of interactions (positive, negative or combinations thereof) and is capable of accommodating any number of microbial species (or guilds) involved. We applied OptCom to quantify the syntrophic association in a well-characterized two-species microbial system, assess the level of sub-optimal growth in phototrophic microbial mats, and elucidate the extent and direction of inter-species metabolite and electron transfer in a model microbial community. We also used OptCom to examine addition of a new member to an existing community. Our study demonstrates the importance of trade-offs between species- and community-level fitness driving forces and lays the foundation for metabolic-driven analysis of various types of interactions in multi-species microbial systems using genome-scale metabolic models.

Integrating flux balance analysis into kinetic models to decipher the dynamic metabolism of Shewanella oneidensis MR-1

Shewanella oneidensis MR-1 sequentially utilizes lactate and its waste products (pyruvate and acetate) during batch culture. To decipher MR-1 metabolism, we integrated genome-scale flux balance analysis (FBA) into a multiple-substrate Monod model to perform the dynamic flux balance analysis (dFBA). The dFBA employed a static optimization approach (SOA) by dividing the batch time into small intervals (i.e., ∼400 mini-FBAs), then the Monod model provided time-dependent inflow/outflow fluxes to constrain the mini-FBAs to profile the pseudo-steady-state fluxes in each time interval. The mini-FBAs used a dual-objective function (a weighted combination of "maximizing growth rate" and "minimizing overall flux") to capture trade-offs between optimal growth and minimal enzyme usage. By fitting the experimental data, a bi-level optimization of dFBA revealed that the optimal weight in the dual-objective function was time-dependent: the objective function was constant in the early growth stage, while the functional weight of minimal enzyme usage increased significantly when lactate became scarce. The dFBA profiled biologically meaningful dynamic MR-1 metabolisms: 1. the oxidative TCA cycle fluxes increased initially and then decreased in the late growth stage; 2. fluxes in the pentose phosphate pathway and gluconeogenesis were stable in the exponential growth period; and 3. the glyoxylate shunt was up-regulated when acetate became the main carbon source for MR-1 growth.

Metabolic flux analysis and visualization

One of the ultimate goals of systems biology research is to obtain a comprehensive understanding of the control mechanisms of complex cellular metabolisms. Metabolic Flux Analysis (MFA) is a important method for the quantitative estimation of intracellular metabolic flows through metabolic pathways and the elucidation of cellular physiology. The primary challenge in the use of MFA is that many biological networks are underdetermined systems; it is therefore difficult to narrow down the solution space from the stoichiometric constraints alone. In this tutorial, we present an overview of Flux Balance Analysis (FBA) and (13)C-Metabolic Flux Analysis ((13)C-MFA), both of which are frequently used to solve such underdetermined systems, and we demonstrate FBA and (13)C-MFA using the genome-scale model and the central carbon metabolism model, respectively. Furthermore, because such comprehensive study of intracellular fluxes is inherently complex, we subsequently introduce various pathway mapping and visualization tools to facilitate understanding of these data in the context of the pathways. Specific visualization of MFA results using the BioCyc Omics Viewer and Pathway Projector are shown as illustrative examples.

Metabolic pathway determination and flux analysis in nonmodel microorganisms through 13C-isotope labeling

C-isotope labeling is a commonly used technique for determining and quantifying pathways in microorganisms under various growth conditions. The experimental protocol consists of feeding the cell with a composition-defined substrate and measuring isotopic labeling patterns in the synthesized metabolites (often the amino acids). Not only can the labeling information be cross-referenced with genomic information to identify the novel pathways, but it can also be used to decipher absolute carbon fluxes through the metabolic network of interest. This technique can be widely used for functional characterization of nonmodel microbial species, and thus we provide a (13)C-pathway and flux analysis protocol. The five key procedures are: (1) growing cells using labeled substrates, (2) measuring extracellular metabolite and biomass component, (3) analyzing isotopic labeling patterns in amino acids and central metabolites using gas chromatography-mass spectrometry, (4) tracing (13)C carbon transitions in metabolites and discovering new pathways, and (5) estimating flux distributions based on isotopomer constraints. This protocol provides complementary information to the recently published protocol for (13)C-based metabolic flux analysis of the model species Escherichia coli (Nat Protoc 4:878-892, 2009).

Synthetic biology triggers new era of antibiotics development

As a discipline to design and construct organisms with desired properties, synthetic biology has generated rapid progresses in the last decade. Combined synthetic biology with the traditional process, a new universal workflow for drug development has been becoming more and more attractive. The new methodology exhibits more efficient and inexpensive comparing to traditional methods in every aspect, such as new compounds discovery & screening, process design & drug manufacturing. This article reviews the application of synthetic biology in antibiotics development, including new drug discovery and screening, combinatorial biosynthesis to generate more analogues and heterologous expression of biosynthetic gene clusters with systematic engineering the recombinant microbial systems for large scale production.

Parameter estimation using meta-heuristics in systems biology: a comprehensive review

This paper gives a comprehensive review of the application of meta-heuristics to optimization problems in systems biology, mainly focussing on the parameter estimation problem (also called the inverse problem or model calibration). It is intended for either the system biologist who wishes to learn more about the various optimization techniques available and/or the meta-heuristic optimizer who is interested in applying such techniques to problems in systems biology. First, the parameter estimation problems emerging from different areas of systems biology are described from the point of view of machine learning. Brief descriptions of various meta-heuristics developed for these problems follow, along with outlines of their advantages and disadvantages. Several important issues in applying meta-heuristics to the systems biology modelling problem are addressed, including the reliability and identifiability of model parameters, optimal design of experiments, and so on. Finally, we highlight some possible future research directions in this field.

Genome-scale consequences of cofactor balancing in engineered pentose utilization pathways in Saccharomyces cerevisiae

Biofuels derived from lignocellulosic biomass offer promising alternative renewable energy sources for transportation fuels. Significant effort has been made to engineer Saccharomyces cerevisiae to efficiently ferment pentose sugars such as D-xylose and L-arabinose into biofuels such as ethanol through heterologous expression of the fungal D-xylose and L-arabinose pathways. However, one of the major bottlenecks in these fungal pathways is that the cofactors are not balanced, which contributes to inefficient utilization of pentose sugars. We utilized a genome-scale model of S. cerevisiae to predict the maximal achievable growth rate for cofactor balanced and imbalanced D-xylose and L-arabinose utilization pathways. Dynamic flux balance analysis (DFBA) was used to simulate batch fermentation of glucose, D-xylose, and L-arabinose. The dynamic models and experimental results are in good agreement for the wild type and for the engineered D-xylose utilization pathway. Cofactor balancing the engineered D-xylose and L-arabinose utilization pathways simulated an increase in ethanol batch production of 24.7% while simultaneously reducing the predicted substrate utilization time by 70%. Furthermore, the effects of cofactor balancing the engineered pentose utilization pathways were evaluated throughout the genome-scale metabolic network. This work not only provides new insights to the global network effects of cofactor balancing but also provides useful guidelines for engineering a recombinant yeast strain with cofactor balanced engineered pathways that efficiently co-utilizes pentose and hexose sugars for biofuels production. Experimental switching of cofactor usage in enzymes has been demonstrated, but is a time-consuming effort. Therefore, systems biology models that can predict the likely outcome of such strain engineering efforts are highly useful for motivating which efforts are likely to be worth the significant time investment.

Dynamic flux balance modeling of S. cerevisiae and E. coli co-cultures for efficient consumption of glucose/xylose mixtures

Current researches into the production of biochemicals from lignocellulosic feedstocks are focused on the identification and engineering of individual microbes that utilize complex sugar mixtures. Microbial consortia represent an alternative approach that has the potential to better exploit individual species capabilities for substrate uptake and biochemical production. In this work, we construct and experimentally validate a dynamic flux balance model of a Saccharomyces cerevisiae and Escherichia coli co-culture designed for efficient aerobic consumption of glucose/xylose mixtures. Each microbe is a substrate specialist, with wild-type S. cerevisiae consuming only glucose and engineered E. coli strain ZSC113 consuming only xylose, to avoid diauxic growth commonly observed in individual microbes. Following experimental identification of a common pH and temperature for optimal co-culture batch growth, we demonstrate that pure culture models developed for optimal growth conditions can be adapted to the suboptimal, common growth condition by adjustment of the non-growth associated ATP maintenance of each microbe. By comparing pure culture model predictions to co-culture experimental data, the inhibitory effect of ethanol produced by S. cerevisiae on E. coli growth was found to be the only interaction necessary to include in the co-culture model to generate accurate batch profile predictions. Co-culture model utility was demonstrated by predicting initial cell concentrations that yield simultaneous glucose and xylose exhaustion for different sugar mixtures. Successful experimental validation of the model predictions demonstrated that steady-state metabolic reconstructions developed for individual microbes can be adapted to develop dynamic flux balance models of microbial consortia for the production of renewable chemicals.

Applications of systems biology towards microbial fuel production

Harnessing the immense natural diversity of biological functions for economical production of fuel has enormous potential benefits. Inevitably, however, the native capabilities for any given organism must be modified to increase the productivity or efficiency of a biofuel bioprocess. From a broad perspective, the challenge is to sufficiently understand the details of cellular functionality to be able to prospectively predict and modify the cellular function of a microorganism. Recent advances in experimental and computational systems biology approaches can be used to better understand cellular level function and guide future experiments. With pressure to quickly develop viable, renewable biofuel processes a balance must be maintained between obtaining depth of biological knowledge and applying that knowledge.

Flux-balance modeling of plant metabolism

Flux-balance modeling of plant metabolic networks provides an important complement to (13)C-based metabolic flux analysis. Flux-balance modeling is a constraints-based approach in which steady-state fluxes in a metabolic network are predicted by using optimization algorithms within an experimentally bounded solution space. In the last 2 years several flux-balance models of plant metabolism have been published including genome-scale models of Arabidopsis metabolism. In this review we consider what has been learnt from these models. In addition, we consider the limitations of flux-balance modeling and identify the main challenges to generating improved and more detailed models of plant metabolism at tissue- and cell-specific scales. Finally we discuss the types of question that flux-balance modeling is well suited to address and its potential role in metabolic engineering and crop improvement.

Fifteen years of large scale metabolic modeling of yeast: developments and impacts

Since the first large-scale reconstruction of the Saccharomyces cerevisiae metabolic network 15 years ago the development of yeast metabolic models has progressed rapidly, resulting in no less than nine different yeast genome-scale metabolic models. Here we review the historical development of large-scale mathematical modeling of yeast metabolism and the growing scope and impact of applications of these models in four different areas: as guide for metabolic engineering and strain improvement, as a tool for biological interpretation and discovery, applications of novel computational framework and for evolutionary studies.

Systems-biology approaches for predicting genomic evolution

Is evolution predictable at the molecular level? The ambitious goal to answer this question requires an understanding of the mutational effects that govern the complex relationship between genotype and phenotype. In practice, it involves integrating systems-biology modelling, microbial laboratory evolution experiments and large-scale mutational analyses - a feat that is made possible by the recent availability of the necessary computational tools and experimental techniques. This Review investigates recent progresses in mapping evolutionary trajectories and discusses the degree to which these predictions are realistic.

Bacterial growth rate reflects a bottleneck in resource allocation

BACKGROUND

Growth rate management in fast-growing bacteria is currently an active research area. In spite of the huge progress made in our understanding of the molecular mechanisms controlling the growth rate, fundamental questions concerning its intrinsic limitations are still relevant today. In parallel, systems biology claims that mathematical models could shed light on these questions.

METHODS

This review explores some possible reasons for the limitation of the growth rate in fast-growing bacteria, using a systems biology approach based on constraint-based modeling methods.

RESULTS

Recent experimental results and a new constraint-based modelling method named Resource Balance Analysis (RBA) reveal the existence of constraints on resource allocation between biological processes in bacterial cells. In this context, the distribution of a finite amount of resources between the metabolic network and the ribosomes limits the growth rate, which implies the existence of a bottleneck between these two processes. Any mechanism for saving resources increases the growth rate.

GENERAL SIGNIFICANCE

Consequently, the emergence of genetic regulation of metabolic pathways, e.g. catabolite repression, could then arise as a means to minimise the protein cost, i.e. maximising growth performance while minimising the resource allocation. This article is part of a Special Issue entitled Systems Biology of Microorganisms.

Optimization-based metabolic control analysis

In this work, a novel optimization-based metabolic control analysis (OMCA) method is introduced for reducing data requirement for metabolic control analysis (MCA). It is postulated that using the optimal control approach, the fluxes in a metabolic network are correlated to metabolite concentrations and enzyme activities as a state-feedback control system that is optimal with respect to a homeostasis objective. It is then shown that the optimal feedback gains are directly related to the elasticity coefficients (ECs) of MCA. This approach requires determination of the relative "importance" of metabolites and fluxes for the system, which is possible with significantly reduced experimental data, as compared with typical MCA requirements. The OMCA approach is applied to a top-down control model of glycolysis in hepatocytes. It is statistically demonstrated that the OMCA model is capable of predicting the ECs observed experimentally with few exceptions. Further, an OMCA-based model reconciliation study shows that the modification of four assumed stoichiometric coefficients in the model can explain most of the discrepancies, with the exception of elasticities with respect to the NADH/NAD ratio.

Expanding a dynamic flux balance model of yeast fermentation to genome-scale

Background

Yeast is considered to be a workhorse of the biotechnology industry for the production of many value-added chemicals, alcoholic beverages and biofuels. Optimization of the fermentation is a challenging task that greatly benefits from dynamic models able to accurately describe and predict the fermentation profile and resulting products under different genetic and environmental conditions. In this article, we developed and validated a genome-scale dynamic flux balance model, using experimentally determined kinetic constraints.

Results

Appropriate equations for maintenance, biomass composition, anaerobic metabolism and nutrient uptake are key to improve model performance, especially for predicting glycerol and ethanol synthesis. Prediction profiles of synthesis and consumption of the main metabolites involved in alcoholic fermentation closely agreed with experimental data obtained from numerous lab and industrial fermentations under different environmental conditions. Finally, fermentation simulations of genetically engineered yeasts closely reproduced previously reported experimental results regarding final concentrations of the main fermentation products such as ethanol and glycerol.

Conclusion

A useful tool to describe, understand and predict metabolite production in batch yeast cultures was developed. The resulting model, if used wisely, could help to search for new metabolic engineering strategies to manage ethanol content in batch fermentations.

Rapid identification of bacteria with a disposable colorimetric sensing array

Rapid identification of both species and even specific strains of human pathogenic bacteria grown on standard agar has been achieved from the volatiles they produce using a disposable colorimetric sensor array in a Petri dish imaged with an inexpensive scanner. All 10 strains of bacteria tested, including Enterococcus faecalis and Staphylococcus aureus and their antibiotic-resistant forms, were identified with 98.8% accuracy within 10 h, a clinically important time frame. Furthermore, the colorimetric sensor arrays also proved useful as a simple research tool for the study of bacterial metabolism and as an easy method for the optimization of bacterial production of fine chemicals or other fermentation processes.

Dynamic flux balance modeling of microbial co-cultures for efficient batch fermentation of glucose and xylose mixtures

Sequential uptake of pentose and hexose sugars that compose lignocellulosic biomass limits the ability of pure microbial cultures to efficiently produce value-added bioproducts. In this work, we used dynamic flux balance modeling to examine the capability of mixed cultures of substrate-selective microbes to improve the utilization of glucose/xylose mixtures and to convert these mixed substrates into products. Co-culture simulations of Escherichia coli strains ALS1008 and ZSC113, engineered for glucose and xylose only uptake respectively, indicated that improvements in batch substrate consumption observed in previous experimental studies resulted primarily from an increase in ZSC113 xylose uptake relative to wild-type E. coli. The E. coli strain ZSC113 engineered for the elimination of glucose uptake was computationally co-cultured with wild-type Saccharomyces cerevisiae, which can only metabolize glucose, to determine if the co-culture was capable of enhanced ethanol production compared to pure cultures of wild-type E. coli and the S. cerevisiae strain RWB218 engineered for combined glucose and xylose uptake. Under the simplifying assumption that both microbes grow optimally under common environmental conditions, optimization of the strain inoculum and the aerobic to anaerobic switching time produced an almost twofold increase in ethanol productivity over the pure cultures. To examine the effect of reduced strain growth rates at non-optimal pH and temperature values, a break even analysis was performed to determine possible reductions in individual strain substrate uptake rates that resulted in the same predicted ethanol productivity as the best pure culture.

Prospects for systems biology and modeling of the gut microbiome

Abundant microorganisms that inhabit the human intestine are implicated in health and disease. The gut microbiome has been studied with metagenomic tools, and over 3 million genes have been discovered, constituting a 'parts list' of this ecosystem; further understanding requires studies of the interacting parts. Mouse models have provided a glimpse into the microbiota and host interactions at metabolic and immunologic levels; however, to provide more insight, there is a need to generate mathematical models that can reveal genotype-phenotype relationships and provide scaffolds for integrated analyses. To this end, we propose the use of genome-scale metabolic models that have successfully been used in studying interactions between human hosts and microbes, as well as microbes in isolation and in communities.

The application of flux balance analysis in systems biology

An increasing number of genome-scale reconstructions of intracellular biochemical networks are being generated. Coupled with these stoichiometric models, several systems-based approaches for probing these reconstructions in silico have been developed. One such approach, called flux balance analysis (FBA), has been effective at predicting systemic phenotypes in the form of fluxes through a reaction network. FBA employs a linear programming (LP) strategy to generate a flux distribution that is optimized toward a particular 'objective,' subject to a set of underlying physicochemical and thermodynamic constraints. Although classical FBA assumes steady-state conditions, several extensions have been proposed in recent years to constrain the allowable flux distributions and enable characterization of dynamic profiles even with minimal kinetic information. Furthermore, FBA coupled with techniques for measuring fluxes in vivo has facilitated integration of computational and experimental approaches, and is allowing pursuit of rational hypothesis-driven research. Ultimately, as we will describe in this review, studying intracellular reaction fluxes allows us to understand network structure and function and has broad applications ranging from metabolic engineering to drug discovery.

Genome-scale metabolic networks

During the last decade, models have been developed to characterize cellular metabolism at the level of an entire metabolic network. The main concept that underlies whole-network metabolic modeling is the identification and mathematical definition of constraints. Here, we review large-scale metabolic network modeling, in particular, stoichiometric- and constraint-based approaches. Although many such models have been reconstructed, few networks have been extensively validated and tested experimentally, and we focus on these. We describe how metabolic networks can be represented using stoichiometric matrices and well-defined constraints on metabolic fluxes. We then discuss relatively successful approaches, including flux balance analysis (FBA), pathway analysis, and common extensions or modifications to these approaches. Finally, we describe techniques for integrating these approaches with models of other biological processes.

Bridging the gap between fluxomics and industrial biotechnology

Metabolic flux analysis is a vital tool used to determine the ultimate output of cellular metabolism and thus detect biotechnologically relevant bottlenecks in productivity. ¹³C-based metabolic flux analysis (¹³C-MFA) and flux balance analysis (FBA) have many potential applications in biotechnology. However, noteworthy hurdles in fluxomics study are still present. First, several technical difficulties in both ¹³C-MFA and FBA severely limit the scope of fluxomics findings and the applicability of obtained metabolic information. Second, the complexity of metabolic regulation poses a great challenge for precise prediction and analysis of metabolic networks, as there are gaps between fluxomics results and other omics studies. Third, despite identified metabolic bottlenecks or sources of host stress from product synthesis, it remains difficult to overcome inherent metabolic robustness or to efficiently import and express nonnative pathways. Fourth, product yields often decrease as the number of enzymatic steps increases. Such decrease in yield may not be caused by rate-limiting enzymes, but rather is accumulated through each enzymatic reaction. Fifth, a high-throughput fluxomics tool hasnot been developed for characterizing nonmodel microorganisms and maximizing their application in industrial biotechnology. Refining fluxomics tools and understanding these obstacles will improve our ability to engineer highlyefficient metabolic pathways in microbial hosts.

Whole-genome metabolic network reconstruction and constraint-based modeling

With the advent of modern high-throughput genomics, there is a significant need for genome-scale analysis techniques that can assist in complex systems analysis. Metabolic genome-scale network reconstructions (GENREs) paired with constraint-based modeling are an efficient method to integrate genomics, transcriptomics, and proteomics to conduct organism-specific analysis. This text explains key steps in the GENRE construction process and several methods of constraint-based modeling that can help elucidate basic life processes and development of disease treatment, bioenergy solutions, and industrial bioproduction applications.