Flux Balance Analysis

Dynamic flux balance analysis (DFBA) provides a platform for detailed design, control and optimization of biochemical process technologies. It is a promising modeling framework that combines genome-scale metabolic network analysis with dynamic simulation of the extracellular environment. Dynamic flux balance analysis assumes that the intracellular species concentrations are in equilibrium with the extracellular environment. The resulting underdetermined stoichiometric model is solved under the assumption of a biochemical objective such as growth rate maximization. The model of the metabolism is coupled with the dynamic mass balance equations of the extracellular environment via expressions for the rates of substrate uptake and product excretion, which imposes additional constraints on the linear program (LP) defined by growth rate maximization of the metabolism. The linear program is embedded into the dynamic model of the bioreactor, and together with the additional constraints this provides an accurate model of the substrate consumption, product secretion, and biomass production during operation. A DFBA model consists of a system of ordinary differential equations for which the evaluation of the right-hand side requires not only function evaluations, but also the solution of one or more linear programs. The numerical tool presented here accurately and efficiently simulates large-scale dynamic flux balance models. The main advantages that this approach has over existing implementation are that the integration scheme has a variable step size, that the linear program only has to be solved when qualitative changes in the optimal flux distribution of the metabolic network occur, and that it can reliably simulate behavior near the boundary of the domain where the model is defined. This is illustrated through large-scale examples taken from the literature. Biotechnol. Bioeng. 2013; 110: 792–802. © 2012 Wiley Periodicals, Inc.

En adoptant, en 1992, la Convention sur la diversité biologique (CDB), la communauté internationale a endossé les objectifs de conservation de la biodiversité, de son utilisation durable et du partage équitable des avantages qui découlent de l’utilisation des ressources génétiques comme un enjeu de gouvernance mondiale.
La diversité biologique joue un rôle essentiel dans le bon fonctionnement des écosystèmes qui est lui-même indispensable au bien-être humain. Malgré son importance fondamentale, reconnue par l’adoption de la CDB au Sommet de la Terre de Rio de Janeiro, la biodiversité continue de s’appauvrir. Reconnaissant le besoin urgent d’agir, l’Assemblée générale des Nations unies a déclaré la période 2011-2020 « Décennie des Nations unies pour la biodiversité ». Les Parties à la Convention ont adopté en 2010 un Plan stratégique pour la biodiversité 2011-2020 comprenant vingt objectifs ambitieux (dits objectifs d’Aichi), qui constitue un « cadre flexible pour la mise en œuvre cohérente d’objectifs nationaux et régionaux ». Consacré au financement du plan stratégique, le dernier objectif s’inscrit dans une Stratégie de mobilisation de(s) ressources dont le chiffrage avait été précisé en octobre 2014, à l’occasion de la prochaine conférence des parties (COP 12) tenue à Pyongchang en Corée du Sud.
A partir du constat de l’insuffisance des seules ressources de l’aide publique au développement (APD) pour la réalisation à l’horizon 2015 des Objectifs du millénaire pour le développement (OMD, 2000) – au rang desquels figure la diminution du taux de perte de biodiversité (objectif 7) – la notion de financements innovants par les budgets nationaux a émergée comme sources ou leviers de mobilisation de ressources complémentaires.
Le Maroc s’est engagé dans la voix de l’actualisation de sa stratégie de la biodiversité et l’élaboration d’un Plan d’action et la mobilisation des ressources pour remplir ses engagements internationaux en tant que Partie à la Convention biodiversité.
Concernant l’intervention de l’économiste, outre la sensibilisation à la valeur de la biodiversité (rapport I) et son intégration dans les systèmes de planification, de comptabilité et de suivi (Rapport II) sont ciblés dans ce rapport III l’évaluation économique de la biodiversité les mécanismes de mobilisation de financement (Les marchés verts, les régimes d’incitations et de subvention, la (sur) compensation).
1. Objet et objectifs de l’étude
Ce dernier et troisième rapport est consacré au plan de mobilisation des ressources pour le financement de la stratégie et des activités du plan d’action national de la biodiversité 2011 - 2020.
L’étude interviendra au niveau des buts D et E de la stratégie nationale pour répondre aux objectifs D4 relatif à la mobilisation des ressources financières et E1 relatif à améliorer l’évaluation des valeurs économiques.
Les données financières pour la période 2006-2010 ont servi à l’établissement d’une valeur de référence (baseline) et les données financières enregistrées après 2010 concernant les besoins et les perspectives de financement à 2020 ont servi à élaborer la note de communication de la CDB formant cadre de suivre des progrès du pays et ce conformément aux indicateurs adoptés dans la décision X/3.

This review focuses on how system biology may assist techniques that are used in pharmacological research, such as high throughput screening, high throughput analytical characterization of biological samples, preclinical and clinical trials, as well as targets and drug validation in order to reach patients at the lowest possible cost in a translational perspective. In signaling networks, targets can be assessed through topological criteria such as their connectivity or/and centrality. In metabolic networks, the relevance of a target for drug development may rather be assessed through some sort of enzymatic specificity resulting from remote homology, analogy or specificity in its strict sense. The concept of specificity is especially valuable in the context of a host-parasite relationship where targeting a protein specific of a parasite compared to its host is expected to minimize the noxious collateral effects of the inhibitor to the host. The relevance of putative molecular target must be proven through bench and animal validations prior to go through clinical trials. Flux balance analysis and other modeling methods of system biology enable to assess whether a molecular target can be considered as pathway's choke or not in a network context, which may facilitate the decision of developing a drugs for it.

In this study a step-wise optimization procedure was developed to predict solvent production using continuous ABE fermentation with immobilized cells. The modeling approach presented here utilizes previously published constraint-based metabolic model for Clostridium acetobutylicum without direct flux constraints. A recently developed flux ratio constraint method was adopted for the model. An experimental data set consisting of 25 experiments using different sugar mixtures as substrates and differing dilution rates was simulated successfully with the modeling approach. Converted to end product concentrations the mean absolute error for acetone was 0.31g/l, for butanol 0.49g/l, and for ethanol 0.17g/l. The modeling approach was validated with another data set from similar experimental setup. The model errors for the validation data set was 0.24g/l, 0.60g/l, and 0.17g/l for acetone, butanol, and ethanol, respectively.

After hundreds of generations of adaptive evolution at exponential growth, Escherichia coli grows as predicted using flux balance analysis (FBA) on genome-scale metabolic models (GEMs). However, it is not known whether the predicted pathway usage in FBA solutions is consistent with gene and protein expression in the wild-type and evolved strains. Here, we report that 498% of active reactions from FBA optimal growth solutions are supported by transcriptomic and proteomic data. Moreover, when E. coli adapts to growth rate selective pressure, the evolved strains upregulate genes within the optimal growth predictions, and downregulate genes outside of the optimal growth solutions. In addition, bottlenecks from dosage limitations of computationally predicted essential genes are overcome in the evolved strains. We also identify regulatory processes that may contribute to the development of the optimal growth phenotype in the evolved strains, such as the downregulation of known regulons and stringent response suppression. Thus, differential gene and protein expression from wild-type and adaptively evolved strains supports observed growth phenotype changes, and is consistent with GEM-computed optimal growth states.

Over the past decade, a growing community of researchers has emerged around the use of constraint-based reconstruction and analysis (COBRA) methods to simulate, analyze and predict a variety of metabolic phenotypes using genome-scale models. The COBRA Toolbox, a MATLAB package for implementing COBRA methods, was presented earlier. Here we present a substantial update of this in silico toolbox. Version 2.0

Background Infections with Salmonella cause significant morbidity and mortality worldwide. Replication of Salmonella typhimurium inside its host cell is a model system for studying the pathogenesis of intracellular bacterial infections. Genome-scale modeling of bacterial metabolic networks provides a powerful tool to identify and analyze pathways required for successful intracellular replication during host-pathogen interaction. Results We have developed and validated a genome-scale metabolic network of Salmonella typhimurium LT2 (iRR1083). This model accounts for 1,083 genes that encode proteins catalyzing 1,087 unique metabolic and transport reactions in the bacterium. We employed flux balance analysis and in silico gene essentiality analysis to investigate growth under a wide range of conditions that mimic in vitro and host cell environments. Gene expression profiling of S. typhimurium isolated from macrophage cell lines was used to constrain the model to predict metabolic pathways that are likely to be operational during infection. Conclusion Our analysis suggests that there is a robust minimal set of metabolic pathways that is required for successful replication of Salmonella inside the host cell. This model also serves as platform for the integration of high-throughput data. Its computational power allows identification of networked metabolic pathways and generation of hypotheses about metabolism during infection, which might be used for the rational design of novel antibiotics or vaccine strains.

Background: Over the last few years a number of methods have been proposed for the phenotype simulation of microorganisms under different environmental and genetic conditions. These have been used as the basis to support the discovery of successful genetic modifications of the microbial metabolism to address industrial goals. However, the use of these methods has been restricted to bioinformaticians or other expert researchers. The main aim of this work is, therefore, to provide a user-friendly computational tool for Metabolic Engineering applications.

Genome-scale metabolic maps can be reconstructed from annotated genome sequence data, biochemical literature, bioinformatic analysis, and strain-specific information. Flux-balance analysis has been useful for qualitative and quantitative analysis of metabolic reconstructions. In the past, FBA has typically been performed in one growth condition at a time, thus giving a limited view of the metabolic capabilities of a metabolic network. We have broadened the use of FBA to map the optimal metabolic flux distribution onto a single plane, which is defined by the availability of two key substrates. A finite number of qualitatively distinct patterns of metabolic pathway utilization were identified in this plane, dividing it into discrete phases. The characteristics of these distinct phases are interpreted using ratios of shadow prices in the form of isoclines. The isoclines can be used to classify the state of the metabolic network. This methodology gives rise to a “phase plane” analysis of the metabolic genotype–phenotype relation relevant for a range of growth conditions. Phenotype phase planes (PhPPs) were generated for Escherichia coli growth on two carbon sources (acetate and glucose) at all levels of oxygenation, and the resulting optimal metabolic phenotypes were studied. Supplementary information can be downloaded from our website (http://epicurus.che.udel.edu). © 2002 John Wiley & Sons, Inc. Biotechnol Bioeng 77: 27–36, 2002.

A genome-scale metabolic model was used to identify new target genes for enhanced biosynthesis of sesquiterpenes in the yeast Saccharomyces cerevisiae. The effect of gene deletions on the flux distributions in the metabolic model of S. cerevisiae was assessed using OptGene as the modeling framework and minimization of metabolic adjustments (MOMA) as objective function.

infections, it remains a major hospital-and community-acquired opportunistic pathogen. With the emergence of multidrug-resistant S. aureus strains, there is an urgent need for the discovery of new antimicrobial drug targets in the organism. To this end, we reconstructed the metabolic networks of multidrug-resistant S. aureus strains using genome annotation, functional-pathway analysis, and comparative genomic approaches, followed by flux balance analysis-based in silico single and double gene deletion experiments. We identified 70 single enzymes and 54 pairs of enzymes whose corresponding metabolic reactions are predicted to be unconditionally essential for growth. Of these, 44 single enzymes and 10 enzyme pairs proved to be common to all 13 S. aureus strains, including many that had not been previously identified as being essential for growth by gene deletion experiments in S. aureus. We thus conclude that metabolic reconstruction and in silico analyses of multiple strains of the same bacterial species provide a novel approach for potential antibiotic target identification.

Isolated liver perfusion systems have been used to characterize intrinsic metabolic changes in liver as a result of various perturbations, including systemic injury, hepatotoxin exposure, and warm ischemia. Most of these studies were done using hyperoxic conditions (95% O 2 ) but without the use of oxygen carriers in the perfusate. Prior literature data do not clearly establish the impact of oxygenation, and in particular that of adding oxygen carriers to the perfusate, on the metabolic functions of the liver. Therefore, herein the effects of oxygen delivery in the perfusion system on liver metabolism were investigated by comparing three modes of oxygenation. Rat livers were perfused via the portal and hepatic veins at a constant flow rate of 3 mL/min/g liver in a recirculating perfusion system. In the first group, the perfusate was equilibrated in a membrane oxygenator with room air (21% O 2 ) before entering the liver. In the second group, the perfusate was equilibrated with a 95% O 2 /5% CO 2 gas mixture. In the third group, the perfusate was supplemented with washed bovine red blood cells (RBCs) at 10% hematocrit and also equilibrated with the 95% O 2 /5% CO 2 gas mixture. Oxygen and CO 2 gradients across the liver were measured periodically with a blood gas analyzer. The rate of change in the concentration of major metabolites in the perfusate was measured over time. Net extracellular fluxes were calculated from these measurements and applied to a stoichiometric-based optimization problem to determine the intracellular fluxes and active pathways in the perfused livers. Livers perfused with RBCs consumed oxygen at twice the rate observed using hyperoxic (95% O 2 ) perfusate without RBCs, and also produced more urea and ketone bodies. At the flow rate used, the oxygen supply in perfusate without RBCs was just sufficient to meet the average oxygen demand of the liver but would be insufficient if it increased above baseline, as is often the case in response to environmental perturbations. Metabolic pathway analysis suggests that significant anaerobic glycolysis occurred in the absence of RBCs even using hyperoxic perfusate. Conversely, when RBCs were used, glucose production from lactate and glutamate, as well as pathways related to energy metabolism were upregulated. RBCs also reversed an increase in PPP fluxes induced by the use of hyperoxic perfusate alone. In conclusion, the use of oxygen carriers is required to investigate the effect of various perturbations on liver metabolism.

Background: Tuberculosis still remains one of the largest killer infectious diseases, warranting the identification of newer targets and drugs. Identification and validation of appropriate targets for designing drugs are critical steps in drug discovery, which are at present major bottle-necks. A majority of drugs in current clinical use for many diseases have been designed without the knowledge of the targets, perhaps because standard methodologies to identify such targets in a high-throughput fashion do not really exist. With different kinds of 'omics' data that are now available, computational approaches can be powerful means of obtaining short-lists of possible targets for further experimental validation.

Synthetic biology, which involves the design of genetic circuits, has been valuable for understanding and engineering biological systems. Such genetic circuits could be valuable in metabolic engineering applications, where the goal is to manipulate the machinery of the organism in order to force or improve the bioengineering objective (which is typically in conflict with the organism's objective governed by evolution, i.e. to optimally allocate the resources for growth).

Genetically modified Saccharomyces cerevisiae strain (YPB-G) which secretes a bifunctional fusion protein that contains both Bacillus subtilis ␣-amylase and Aspergillus awamori glucoamylase activities was used for the direct conversion of starch into ethanol. Starch was either supplied initially to different nutrient media or added instantaneously to the reactor at various discrete time instants (pulse feeding). Stoichiometric modeling was used to investigate the effects of initial substrate concentration and growth rate of the recombinant yeast culture on ethanol production. Reaction stoichiometries describing both the anabolism and catabolism of the microorganism were used as an input to flux balance analysis (FBA), the preferred metabolic modeling approach since the constructed stoichiometric network was underdetermined. Experiments for batch and fed-batch systems at different substrate concentrations were analyzed theoretically in terms of flux distributions using ethanol production rate as the maximization criteria. Calculated ethanol rates were in agreement with experimental measurements, suggesting that this recombinant microorganism is sufficiently evolved to optimize its ethanol production. The function of the main pathways of yeast metabolism (PPP, EMP, TCA) are discussed together with the node analyses of glucose-6-P and pyruvate branch points. Theoretical node analysis revealed that if the split ratio in G6P branch point is changed by genetic manipulations, the ethanol yield would be affected considerably.

over the past decade, a growing community of researchers has emerged around the use of constraint-based reconstruction and analysis (coBra) methods to simulate, analyze and predict a variety of metabolic phenotypes using genome-scale models. the coBra toolbox, a MatlaB package for implementing coBra methods, was presented earlier. Here we present a substantial update of this in silico toolbox. Version 2.0 of the coBra toolbox expands the scope of computations by including in silico analysis methods developed since its original release. new functions include (i) network gap filling, (ii) 13 c analysis, (iii) metabolic engineering, (iv) omics-guided analysis and (v) visualization. as with the first version, the coBra toolbox reads and writes systems biology markup language-formatted models. In version 2.0, we improved performance, usability and the level of documentation. a suite of test scripts can now be used to learn the core functionality of the toolbox and validate results. this toolbox lowers the barrier of entry to use powerful coBra methods.

The elucidation of organism-scale metabolic networks necessitates the development of integrative methods to analyze and interpret the systemic properties of cellular metabolism. A shift in emphasis from single metabolic reactions to systemically defined pathways is one consequence of such an integrative analysis of metabolic systems. The constraints of systemic stoichiometry, and limited thermodynamics have led to the definition of the flux space within the context of convex analysis. The flux space of the metabolic system, containing all allowable flux distributions, is constrained to a convex polyhedral cone in a high-dimensional space. From metabolic pathway analysis, the edges of the highdimensional flux cone are vectors that correspond to systemically defined "extreme pathways" spanning the capabilities of the system. The addition of maximum flux capacities of individual metabolic reactions serves to further constrain the flux space and has led to the development of flux balance analysis using linear optimization to calculate optimal flux distributions. Here we provide the precise theoretical connections between pathway analysis and flux balance analysis allowing for their combined application to study integrated metabolic function. Shifts in metabolic behavior are calculated using linear optimization and are then interpreted using the extreme pathways to demonstrate the concept of pathway utilization. Changes to the reaction network, such as the removal of a reaction, can lead to the generation of suboptimal phenotypes that can be directly attributed to the loss of pathway function and capabilities. Optimal growth phenotypes are calculated as a function of environmental variables, such as the availability of substrate and oxygen, leading to the definition of phenotypic phase planes. It is illustrated how optimality properties of the computed flux distributions can be interpreted in terms of the extreme pathways. Together these developments are applied to an example network and to core metabolism of Escherichia coli demonstrating the connections between the extreme pathways, optimal flux distributions, and phenotypic phase planes. The consequences of changing environmental and internal conditions of the network are examined for growth on glucose and succinate in the face of a variety of gene deletions. The convergence of the calculation of optimal phenotypes through linear programming and the definition of extreme pathways establishes a different perspective for the understanding of how a defined metabolic network is best used under different environmental and internal conditions or, in other words, a pathway basis for the interpretation of the metabolic reaction norm.

Background: Because metabolism is fundamental in sustaining microbial life, drugs that target pathogen-specific metabolic enzymes and pathways can be very effective. In particular, the metabolic challenges faced by intracellular pathogens, such as Mycobacterium tuberculosis, residing in the infected host provide novel opportunities for therapeutic intervention.

A persistent problem hampering our understanding of the dynamics of large-scale metabolic networks is the lack of experimentally determined kinetic parameters that are necessary to build computational models of biochemical processes. To overcome some of the limitations imposed by absent or incomplete kinetic data, structural kinetic modeling (SKM) was proposed recently as an intermediate approach between stoichiometric analysis and a full kinetic description. SKM extends stationary flux-balance analysis (FBA) by a local stability analysis utilizing an appropriate parametrization of the Jacobian matrix. To characterize the Jacobian, we utilize results from robust control theory to determine subintervals of the Jacobian' entries that correspond to asymptotically stable metabolic states. Furthermore, we propose an efficient sampling scheme in combination with methods from computational geometry to sketch the stability region. A glycolytic pathway model comprising 12 uncertain parameters is used to assess the feasibility of the method.

Background: Over the last few years a number of methods have been proposed for the phenotype simulation of microorganisms under different environmental and genetic conditions. These have been used as the basis to support the discovery of successful genetic modifications of the microbial metabolism to address industrial goals. However, the use of these methods has been restricted to bioinformaticians or other expert researchers. The main aim of this work is, therefore, to provide a user-friendly computational tool for Metabolic Engineering applications. Results: OptFlux is an open-source and modular software aimed at being the reference computational application in the field. It is the first tool to incorporate strain optimization tasks, i.e., the identification of Metabolic Engineering targets, using Evolutionary Algorithms/Simulated Annealing metaheuristics or the previously proposed OptKnock algorithm. It also allows the use of stoichiometric metabolic models for (i) phenotype simulation of both wild-type and mutant organisms, using the methods of Flux Balance Analysis, Minimization of Metabolic Adjustment or Regulatory on/off Minimization of Metabolic flux changes, (ii) Metabolic Flux Analysis, computing the admissible flux space given a set of measured fluxes, and (iii) pathway analysis through the calculation of Elementary Flux Modes. OptFlux also contemplates several methods for model simplification and other pre-processing operations aimed at reducing the search space for optimization algorithms. The software supports importing/exporting to several flat file formats and it is compatible with the SBML standard. OptFlux has a visualization module that allows the analysis of the model structure that is compatible with the layout information of Cell Designer, allowing the superimposition of simulation results with the model graph. Conclusions: The OptFlux software is freely available, together with documentation and other resources, thus bridging the gap from research in strain optimization algorithms and the final users. It is a valuable platform for researchers in the field that have available a number of useful tools. Its open-source nature invites contributions by all those interested in making their methods available for the community. Given its plug-in based architecture it can be extended with new functionalities. Currently, several plug-ins are being developed, including network topology analysis tools and the integration with Boolean network based regulatory models.

Metabolic Control Analysis (MCA) is a conceptual and mathematical formalism that models the relative contributions of individual effectors in a pathway to both the flux through the pathway and the concentrations of individual intermediates within it. To exploit MCA in an initial Systems Biology analysis of the eukaryotic cell, two categories of experiments are required. In category 1 experiments, flux is changed and the impact on the levels of the direct and indirect products of gene action is measured. We have measured the impact of changing the flux on the transcriptome, proteome and metabolome of Saccharomyces cerevisiae . In this whole-cell analysis, flux equates to growth rate. In category 2 experiments, the levels of individual gene products are altered, and the impact on the flux is measured. We have used competition analyses between the complete set of heterozygous yeast deletion mutants to reveal genes encoding proteins with high flux control coefficients. These genes may b...

The scope and breadth of genome-scale metabolic reconstructions have continued to expand over the last decade. Herein, we introduce a genome-scale model for a plant with direct applications to food and bioenergy production (i.e., maize). Maize annotation is still underway, which introduces significant challenges in the association of metabolic functions to genes. The developed model is designed to meet rigorous standards on gene-protein-reaction (GPR) associations, elementally and charged balanced reactions and a biomass reaction abstracting the relative contribution of all biomass constituents. The metabolic network contains 1,563 genes and 1,825 metabolites involved in 1,985 reactions from primary and secondary maize metabolism. For approximately 42% of the reactions direct literature evidence for the participation of the reaction in maize was found. As many as 445 reactions and 369 metabolites are unique to the maize model compared to the AraGEM model for A. thaliana. 674 metabolites and 893 reactions are present in Zea mays iRS1563 that are not accounted for in maize C4GEM. All reactions are elementally and charged balanced and localized into six different compartments (i.e., cytoplasm, mitochondrion, plastid, peroxisome, vacuole and extracellular). GPR associations are also established based on the functional annotation information and homology prediction accounting for monofunctional, multifunctional and multimeric proteins, isozymes and protein complexes. We describe results from performing flux balance analysis under different physiological conditions, (i.e., photosynthesis, photorespiration and respiration) of a C4 plant and also explore model predictions against experimental observations for two naturally occurring mutants (i.e., bm1 and bm3). The developed model corresponds to the largest and more complete to-date effort at cataloguing metabolism for a plant species.

Improvement of biological strains through targeted modification of metabolism is essential for successful development of bioprocesses. The computational complexity of optimization procedures routinely used for identifying genetic targets limits their application to genome-scale metabolic networks. In this study, we combined graph theoretic approaches with mixed-integer liner programming (MILP) to reduce the search space and thus reducing computational time. Specifically, we used cut-sets (minimal set of reactions that cuts metabolic networks) as additional constraints to reduce the search space. The efficacy of proposed approach is illustrated by identifying minimal reaction set for Saccharomyces Cerevisiae.

Metabolism is central to cell physiology, and metabolic disturbances play a role in numerous disease states. Despite its importance, the ability to study metabolism at a global scale using genomic technologies is limited. In principle, complete genome sequences describe the range of metabolic reactions that are possible for an organism, but cannot quantitatively describe the behaviour of these reactions. We present a novel method for modeling metabolic states using whole cell measurements of gene expression. Our method, which we call E-Flux (as a combination of flux and expression), extends the technique of Flux Balance Analysis by modeling maximum flux constraints as a function of measured gene expression. In contrast to previous methods for metabolically interpreting gene expression data, E-Flux utilizes a model of the underlying metabolic network to directly predict changes in metabolic flux capacity. We applied E-Flux to Mycobacterium tuberculosis, the bacterium that causes tuberculosis (TB). Key components of mycobacterial cell walls are mycolic acids which are targets for several first-line TB drugs. We used E-Flux to predict the impact of 75 different drugs, drug combinations, and nutrient conditions on mycolic acid biosynthesis capacity in M. tuberculosis, using a public compendium of over 400 expression arrays. We tested our method using a model of mycolic acid biosynthesis as well as on a genome-scale model of M. tuberculosis metabolism. Our method correctly predicts seven of the eight known fatty acid inhibitors in this compendium and makes accurate predictions regarding the specificity of these compounds for fatty acid biosynthesis. Our method also predicts a number of additional potential modulators of TB mycolic acid biosynthesis. E-Flux thus provides a promising new approach for algorithmically predicting metabolic state from gene expression data.

Cell metabolism is the biochemical machinery that provides energy and building blocks to sustain life. Understanding its fine regulation is of pivotal relevance in several fields, from metabolic engineering applications to the treatment of metabolic disorders and cancer. Sophisticated computational approaches are needed to unravel the complexity of metabolism. To this aim, a plethora of methods have been developed, yet it is generally hard to identify which computational strategy is most suited for the investigation of a specific aspect of metabolism. This review provides an up-to-date description of the computational methods available for the analysis of metabolic pathways, discussing their main advantages and drawbacks. In particular, attention is devoted to the identification of the appropriate scale and level of accuracy in the reconstruction of metabolic networks, and to the inference of model structure and parameters, especially when dealing with a shortage of experimental measurements. The choice of the proper computational methods to derive in silico data is then addressed, including topological analyses, constraint-based modeling and simulation of the system dynamics. A description of some computational approaches to gain new biological knowledge or to formulate hypotheses is finally provided.

Poly-3-hydroxybutyrate (PHB) synthesis in Escherichia coli elicits regulatory responses that affect product yield and productivity. We used controlled, steady-state cultures (chemostats) of a genetically stable strain to determine growth-independent metabolic flux regulation. We measured flux and steady-state intracellular metabolite concentrations across different dilution rates (0.05, 0.15, 0.3h(-1)), limitations (glucose, gluconate and nitrogen), and operon copy counts of the PHB pathway (0, 6, 17, and 29). As PHB flux increases, specific substrate consumption and lactate secretion increase while formate and acetate secretion decreases in N-limited, glucose-fed conditions. To understand the regulatory mechanisms that resulted in these macroscopic changes, we used a flux balance analysis model to analyze intracellular redox conditions. Our model shows that under N-limited conditions, synthesis of PHB creates excess reducing equivalents. Cells, under these conditions, secrete more reduced metabolites in order to recycle reducing equivalents. By switching to a more oxidized substrate (gluconate) that decreased excess reducing equivalents, PHB flux yield increased 1.6 fold compared to glucose-fed fermentations. High flux of PHB (~1.2 mmol/g DCWh) was maintained under these steady-state, oxidized conditions. These results imply redox imbalance is a driving force in industrial production of PHB, and substrates that are more oxidized than glucose can increase productivity.

The current state of the art for linear optimization in Flux Balance Analysis has been limited to single objective functions. Since mammalian systems perform various functions, a multiobjective approach is needed when seeking optimal flux distributions in these systems. In most of the available multiobjective optimization methods, there is a lack of understanding of when to use a particular objective, and how to combine and/or prioritize mutually competing objectives to achieve a truly optimal solution. To address these limitations we developed a soft constraints based linear physical programming-based flux balance analysis (LPPFBA) framework to obtain a multiobjective optimal solutions. The developed framework was first applied to compute a set of multiobjective optimal solutions for various pairs of objectives relevant to hepatocyte function (urea secretion, albumin, NADPH, and glutathione syntheses) in bioartificial liver systems. Next, simultaneous analysis of the optimal solutions for three objectives was carried out. Further, this framework was utilized to obtain true optimal conditions to improve the hepatic functions in a simulated bioartificial liver system. The combined quantitative and visualization framework of LPPFBA is applicable to any large-scale metabolic network system, including those derived by genomic analyses.

The reconstruction of genome-scale metabolic models and their applications represent a great advantage of systems biology. Through their use as metabolic flux simulation models, production of industrially-interesting metabolites can be predicted. Due to the growing number of studies of metabolic models driven by the increasing genomic sequencing projects, it is important to conceptualize steps of reconstruction and analysis. We have focused our work in the cyanobacterium Synechococcus elongatus PCC7942, for which several analyses and insights are unveiled. A comprehensive approach has been used, which can be of interest to lead the process of manual curation and genome-scale OPEN ACCESS Metabolites 2014, 4 681 metabolic analysis. The final model, iSyf715 includes 851 reactions and 838 metabolites. A biomass equation, which encompasses elementary building blocks to allow cell growth, is also included. The applicability of the model is finally demonstrated by simulating autotrophic growth conditions of Synechococcus elongatus PCC7942.

Understanding the evolution and dynamics of metabolism in microbial ecosystems is an ongoing challenge in microbiology. A promising approach towards this goal is the extension of genome-scale flux balance models of metabolism to multiple interacting species. However, since the detailed distribution of metabolic functions among ecosystem members is often unknown, it is important to investigate how compartmentalization of metabolites and reactions affects flux balance predictions. Here, as a first step in this direction, we address the importance of compartmentalization in the well characterized metabolic model of the yeast Saccharomyces cerevisiae, which we treat as an "ecosystem of organelles". In addition to addressing the impact that the removal of compartmentalization has on model predictions, we show that by systematically constraining some individual fluxes in a de-compartmentalized version of the model we can significantly reduce the flux prediction errors induced by the removal of compartments. We expect that our analysis will help predict and understand metabolic functions in complex microbial communities. In addition, further study of yeast as an ecosystem of organelles might provide novel insight on the evolution of endosymbiosis and multicellularity.

The reconstruction of genome-scale metabolic models and their applications represent a great advantage of systems biology. Through their use as metabolic flux simulation models, production of industrially-interesting metabolites can be predicted. Due to the growing number of studies of metabolic models driven by the increasing genomic sequencing projects, it is important to conceptualize steps of reconstruction and analysis. We have focused our work in the cyanobacterium Synechococcus elongatus PCC7942, for which several analyses and insights are unveiled. A comprehensive approach has been used, which can be of interest to lead the process of manual curation and genome-scale OPEN ACCESS Metabolites 2014, 4 681 metabolic analysis. The final model, iSyf715 includes 851 reactions and 838 metabolites. A biomass equation, which encompasses elementary building blocks to allow cell growth, is also included. The applicability of the model is finally demonstrated by simulating autotrophic growth conditions of Synechococcus elongatus PCC7942.

Genome-scale metabolic models have proven useful for answering fundamental questions about metabolic capabilities of a variety of microorganisms, as well as informing their metabolic engineering. However, only a few models are available for oxygenic photosynthetic microorganisms, particularly in cyanobacteria in which photosynthetic and respiratory electron transport chains (ETC) share components. We addressed the complexity of cyanobacterial ETC by developing a genomescale model for the diazotrophic cyanobacterium, Cyanothece sp. ATCC 51142. The resulting metabolic reconstruction, iCce806, consists of 806 genes associated with 667 metabolic reactions and includes a detailed representation of the ETC and a biomass equation based on experimental measurements. Both computational and experimental approaches were used to investigate light-driven metabolism in Cyanothece sp. ATCC 51142, with a particular focus on reductant production and partitioning within the ETC. The simulation results suggest that growth and metabolic flux distributions are substantially impacted by the relative amounts of light going into the individual photosystems. When growth is limited by the flux through photosystem I, terminal respiratory oxidases are predicted to be an important mechanism for removing excess reductant. Similarly, under photosystem II flux limitation, excess electron carriers must be removed via cyclic electron transport. Furthermore, in silico calculations were in good quantitative agreement with the measured growth rates whereas predictions of reaction usage were qualitatively consistent with protein and mRNA expression data, which we used to further improve the resolution of intracellular flux values.

Biologically inspired computation has been recently used with mathematical models towards the design of new synthetic organisms. In this work, we use Pareto optimality to optimize these organisms in a multi-objective fashion. We infer the best knockout strategies to perform specific tasks in bacteria, which involve concurrent maximization/minimization of multiple functions (codomain) and optimization of several decision variables (domain). Furthermore, we propose and exploit a mapping between the metabolism and a register machine. We show that optimized bacteria have computational capability and act as molecular Turing machines programmed using a Pareto optimal solution. Finally, we investigate communication between bacteria as a means to evaluate their computational capability. We report that the density and gradient of the Pareto curve are useful tools to compare models and understand their structure, while modelling organisms as computers proves useful to carry out computation using biological machines with specific input-output conditions, as well as to estimate the bacterial computational effort for specific tasks.

Background It is a daunting task to identify all the metabolic pathways of brain energy metabolism and develop a dynamic simulation environment that will cover a time scale ranging from seconds to hours. To simplify this task and make it more practicable, we undertook stoichiometric modeling of brain energy metabolism with the major aim of including the main interacting pathways in and between astrocytes and neurons. Model The constructed model includes central metabolism (glycolysis, pentose phosphate pathway, TCA cycle), lipid metabolism, reactive oxygen species (ROS) detoxification, amino acid metabolism (synthesis and catabolism), the well-known glutamate-glutamine cycle, other coupling reactions between astrocytes and neurons, and neurotransmitter metabolism. This is, to our knowledge, the most comprehensive attempt at stoichiometric modeling of brain metabolism to date in terms of its coverage of a wide range of metabolic pathways. We then attempted to model the basal physiolo...

This review focuses on how system biology may assist techniques that are used in pharmacological research, such as high-throughput screening, high-throughput analytical characterization of biological samples, preclinical and clinical trials, as well as targets and drug validation in order to reach patients at the lowest possible cost in a translational perspective. In signaling networks, targets can be assessed through topological criteria such as their connectivity and/or centrality. In metabolic networks, the relevance of a target for drug development may rather be assessed through some sort of enzymatic specificity resulting from remote homology, analogy, or specificity in its strict sense. The concept of specificity is especially valuable in the context of a host-parasite relationship where targeting a protein specific of a parasite compared to its host is expected to minimize the noxious collateral effects of the inhibitor to the host. The relevance of putative molecular target must be proven through bench and animal validations prior to going through clinical trials. Flux balance analysis and other modeling methods of system biology enable to assess whether a molecular target can be considered as pathway’s choke or not in a network context, which may facilitate the decision of developing drugs for it.

Cell metabolism is the biochemical machinery that provides energy and building blocks to sustain life. Understanding its fine regulation is of pivotal relevance in several fields, from metabolic engineering applications to the treatment of metabolic disorders and cancer. Sophisticated computational approaches are needed to unravel the complexity of metabolism. To this aim, a plethora of methods have been developed, yet it is generally hard to identify which computational strategy is most suited for the investigation of a specific aspect of metabolism. This review provides an up-to-date description of the computational methods available for the analysis of metabolic pathways, discussing their main advantages and drawbacks. In particular, attention is devoted to the identification of the appropriate scale and level of accuracy in the reconstruction of metabolic networks, and to the inference of model structure and parameters, especially when dealing with a shortage of experimental measurements. The choice of the proper computational methods to derive in silico data is then addressed, including topological analyses, constraint-based modeling and simulation of the system dynamics. A description of some computational approaches to gain new biological knowledge or to formulate hypotheses is finally provided.

The identification of a set of genetic manipulations that result in a microbial strain with improved production capabilities of a metabolite with industrial interest is a big challenge in Metabolic Engineering. Evolutionary Algorithms and Simulated Annealing have been used in this task to identify sets of reaction deletions, towards the maximization of a desired objective function. To simulate the cell phenotype for each mutant strain, the Flux Balance Analysis approach is used, assuming organisms have maximized their growth along evolution. In this work, transcriptional information is added to the models using gene-reaction rules. The aim is to find the (near-)optimal set of gene knockouts necessary to reach a given productivity goal. The results obtained are compared with the ones reached using the deletion of reactions, showing that we obtain solutions with similar quality levels and number of knockouts, but biologically more feasible. Indeed, we show that several of the previous solutions are not viable using the provided rules.

Cell metabolism is the biochemical machinery that provides energy and building blocks to sustain life. Understanding its fine regulation is of pivotal relevance in several fields, from metabolic engineering applications to the treatment of metabolic disorders and cancer. Sophisticated computational approaches are needed to unravel the complexity of metabolism. To this aim, a plethora of methods have been developed, yet it is generally hard to identify which computational strategy is most suited for the investigation of a specific aspect of metabolism. This review provides an up-to-date description of the computational methods available for the analysis of metabolic pathways, discussing their main advantages and drawbacks. In particular, attention is devoted to the identification of the appropriate scale and level of accuracy in the reconstruction of metabolic networks, and to the inference of model structure and parameters, especially when dealing with a shortage of experimental measurements. The choice of the proper computational methods to derive in silico data is then addressed, including topological analyses, constraint-based modeling and simulation of the system dynamics. A description of some computational approaches to gain new biological knowledge or to formulate hypotheses is finally provided.

The current state of the art for linear optimization in Flux Balance Analysis has been limited to single objective functions. Since mammalian systems perform various functions, a multiobjective approach is needed when seeking optimal flux distributions in these systems. In most of the available multiobjective optimization methods, there is a lack of understanding of when to use a particular objective, and how to combine and/or prioritize mutually competing objectives to achieve a truly optimal solution. To address these limitations we developed a soft constraints based linear physical programming-based flux balance analysis (LPPFBA) framework to obtain a multiobjective optimal solutions. The developed framework was first applied to compute a set of multiobjective optimal solutions for various pairs of objectives relevant to hepatocyte function (urea secretion, albumin, NADPH, and glutathione syntheses) in bioartificial liver systems. Next, simultaneous analysis of the optimal solutions for three objectives was carried out. Further, this framework was utilized to obtain true optimal conditions to improve the hepatic functions in a simulated bioartificial liver system. The combined quantitative and visualization framework of LPPFBA is applicable to any large-scale metabolic network system, including those derived by genomic analyses.

The human red blood cell (RBC) metabolism is investigated by calculating steady state fluxes using constraint-based stoichiometric modeling approaches. For the normal RBC metabolism, flux balance analysis (FBA) is performed via optimization of various alternative objective functions, and the maximization of production of ATP and NADPH is found to be the primary objective of the RBC metabolism. FBA and two novel approaches, minimization of metabolic adjustment (MOMA) and regulatory on-off minimization (ROOM), which can describe the behavior of the metabolic networks in case of enzymopathies, are applied to observe the relative changes in the flux distribution of the deficient network. The deficiencies in several enzymes in RBC metabolism are investigated and the flux distributions are compared with the non-deficient FBA distribution to elucidate the metabolic changes in response to enzymopathies. It is found that the metabolism is mostly affected by the glucose-6-phosphate dehydrogenase (G6PDH) and phosphoglycerate kinase (PGK) enzymopathies, whereas the effects of the deficiency in DPGM on the metabolism are negligible. These stoichiometric modeling results are found to be in accordance with the experimental findings in the literature related to metabolic behavior of the human red blood cells, showing that human RBC metabolism can be modeled stoichiometrically.