Constraint-based Modeling Approach Research Articles

Microarray data can predict diurnal changes of starch content in the picoalga Ostreococcus

BackgroundThe storage of photosynthetic carbohydrate products such as starch is subject to complex regulation, effected at both transcriptional and post-translational levels. The relevant genes in plants show pronounced daily regulation. Their temporal RNA expression profiles, however, do not predict the dynamics of metabolite levels, due to the divergence of enzyme activity from the RNA profiles.Unicellular phytoplankton retains the complexity of plant carbohydrate metabolism, and recent transcriptomic profiling suggests a major input of transcriptional regulation.ResultsWe used a quasi-steady-state, constraint-based modelling approach to infer the dynamics of starch content during the 12 h light/12 h dark cycle in the model alga Ostreococcus tauri. Measured RNA expression datasets from microarray analysis were integrated with a detailed stoichiometric reconstruction of starch metabolism in O. tauri in order to predict the optimal flux distribution and the dynamics of the starch content in the light/dark cycle. The predicted starch profile was validated by experimental data over the 24 h cycle. The main genetic regulatory targets within the pathway were predicted by in silico analysis.ConclusionsA single-reaction description of starch production is not able to account for the observed variability of diurnal activity profiles of starch-related enzymes. We developed a detailed reaction model of starch metabolism, which, to our knowledge, is the first attempt to describe this polysaccharide polymerization while preserving the mass balance relationships. Our model and method demonstrate the utility of a quasi-steady-state approach for inferring dynamic metabolic information in O. tauri directly from time-series gene expression data.

A constraint-based topic modeling approach for name disambiguation

Name ambiguity refers to a problem that different people might be referenced with an identical name. This problem has become critical in many applications, particularly in online bibliography systems, such as DBLP and CiterSeer. Although much work has been conducted to address this problem, there still exist many challenges. In this paper, a general framework of constraint-based topic modeling is proposed, which can make use of user-defined constraints to enhance the performance of name disambiguation. A Gibbs sampling algorithm that integrates the constraints has been proposed to do the inference of the topic model. Experimental results on a real-world dataset show that significant improvements can be obtained by taking the proposed approach.

Understanding Clostridium acetobutylicum ATCC 824 Metabolism Using Genome-Scale Thermodynamics and Metabolomics-based Modeling

Biobutanol has a potential application as a biofuel and can replace gasoline as an alternative fuel. However, the productivity of the biobutanol process has to be increased significantly, before it can be commercialized. Clostridium acetobutylicum has been the primary microbial host used for biobutanol production. Consequently, systems biology based genome-scale metabolic model of C. acetobutylicum metabolism is essential to optimize the biobutanol production rate via genetic and process optimization and metabolic engineering. An updated genome-scale metabolic model of C. acetobutylicum ATCC824 consists of 700 genes, 709 reactions, 29 exchange reactions, and 679 metabolites, which covers 16.3% of ORFs and is more comprehensive compared to two other existing models. This metabolic network was used as a platform for simulating phenotypes using the constraint-based modeling approach. Flux variability analysis shows the presence of alternate carbon and electron sinks that allows for different carbon assimilation patterns and ranges of product formation rates. Incorporation of metabolomics and thermodynamics data into this model results in thermodynamically feasible flux distributions. We apply an NMR based approach to quantify 31 intracellular metabolites in both metabolic phases. The total molar concentration of metabolites was approximately 245mM, where asparagine was the most dominant metabolite.

Constraint-based modelling of Chlamydomonas reinhardtii photoautotrophic metabolism

A constraint-based modelling approach was developed to investigate the metabolic response of the eukaryotic microalgae Chlamydomonas reinhardtii under photoautotrophic conditions. The model explicitly includes thermodynamic and energetic constraints on the functioning metabolism. A mixed integer linear programming method was used to determine the optimal flux distributions with regard to this set of constraints. It enabled us in particular to highlight the existence of a light-driven respiration depending on the incident photon flux density in photobioreactors.

Evaluating the integration of proteomic data for the prediction of intracellular fluxes after knockout experiments

So far, few large scale kinetic models of metabolic networks have been successfully constructed. The main reasons for this are not only the associated mathematical complexity, but also the large number of unknown kinetic parameters required in the rate equations to define the system. In contrast to kinetic models, the constraint-based modelling approach bypasses these difficulties by using basically only stoichiometric information with certain physicochemical constraints to delimit the solution space without large fitted parameter sets. Although these constraint-based models are highly relevant to predict feasible steady-state fluxes under a diverse range of genetic and environmental conditions, the steady-state assumption may oversimplify cellular behaviour and cannot predict time-course profiles. To overcome these problems, combining these two approaches appears as a reasonable alternative to modelling large-scale metabolic networks. On the other hand, several of the experimental data required for model construction are often rare and in this way it is usually assumed that the enzyme concentrations are constant.In this work, we used a central carbon metabolic network of E. coli to investigate whether including high throughput enzyme concentration data into a kinetic model allows improved predictions of metabolic flux distributions in response to single knockouts perturbations. For this purpose, an E. coli model, based on results obtained from flux balance analysis (FBA) and approximate lin-log kinetics was constructed. The intracellular fluxes distributions, obtained using this model, were compared with published in vivo measurements.

Genome-scale constraint-based modeling of Geobacter metallireducens

BackgroundGeobacter metallireducens was the first organism that can be grown in pure culture to completely oxidize organic compounds with Fe(III) oxide serving as electron acceptor. Geobacter species, including G. sulfurreducens and G. metallireducens, are used for bioremediation and electricity generation from waste organic matter and renewable biomass. The constraint-based modeling approach enables the development of genome-scale in silico models that can predict the behavior of complex biological systems and their responses to the environments. Such a modeling approach was applied to provide physiological and ecological insights on the metabolism of G. metallireducens.ResultsThe genome-scale metabolic model of G. metallireducens was constructed to include 747 genes and 697 reactions. Compared to the G. sulfurreducens model, the G. metallireducens metabolic model contains 118 unique reactions that reflect many of G. metallireducens' specific metabolic capabilities. Detailed examination of the G. metallireducens model suggests that its central metabolism contains several energy-inefficient reactions that are not present in the G. sulfurreducens model. Experimental biomass yield of G. metallireducens growing on pyruvate was lower than the predicted optimal biomass yield. Microarray data of G. metallireducens growing with benzoate and acetate indicated that genes encoding these energy-inefficient reactions were up-regulated by benzoate. These results suggested that the energy-inefficient reactions were likely turned off during G. metallireducens growth with acetate for optimal biomass yield, but were up-regulated during growth with complex electron donors such as benzoate for rapid energy generation. Furthermore, several computational modeling approaches were applied to accelerate G. metallireducens research. For example, growth of G. metallireducens with different electron donors and electron acceptors were studied using the genome-scale metabolic model, which provided a fast and cost-effective way to understand the metabolism of G. metallireducens.ConclusionWe have developed a genome-scale metabolic model for G. metallireducens that features both metabolic similarities and differences to the published model for its close relative, G. sulfurreducens. Together these metabolic models provide an important resource for improving strategies on bioremediation and bioenergy generation.

Analysis of Growth of Lactobacillus plantarum WCFS1 on a Complex Medium Using a Genome-scale Metabolic Model

A genome-scale metabolic model of the lactic acid bacterium Lactobacillus plantarum WCFS1 was constructed based on genomic content and experimental data. The complete model includes 721 genes, 643 reactions, and 531 metabolites. Different stoichiometric modeling techniques were used for interpretation of complex fermentation data, as L. plantarum is adapted to nutrient-rich environments and only grows in media supplemented with vitamins and amino acids. (i) Based on experimental input and output fluxes, maximal ATP production was estimated and related to growth rate. (ii) Optimization of ATP production further identified amino acid catabolic pathways that were not previously associated with free-energy metabolism. (iii) Genome-scale elementary flux mode analysis identified 28 potential futile cycles. (iv) Flux variability analysis supplemented the elementary mode analysis in identifying parallel pathways, e.g. pathways with identical end products but different co-factor usage. Strongly increased flexibility in the metabolic network was observed when strict coupling between catabolic ATP production and anabolic consumption was relaxed. These results illustrate how a genome-scale metabolic model and associated constraint-based modeling techniques can be used to analyze the physiology of growth on a complex medium rather than a minimal salts medium. However, optimization of biomass formation using the Flux Balance Analysis approach, reported to successfully predict growth rate and by product formation in Escherichia coli and Saccharomyces cerevisiae, predicted too high biomass yields that were incompatible with the observed lactate production. The reason is that this approach assumes optimal efficiency of substrate to biomass conversion, and can therefore not predict the metabolically inefficient lactate formation.

Analysis of enzymopathies in the human red blood cells by constraint-based stoichiometric modeling approaches

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.

Characterization of Metabolism in the Fe(III)-Reducing Organism Geobacter sulfurreducens by Constraint-Based Modeling

Geobacter sulfurreducens is a well-studied representative of the Geobacteraceae, which play a critical role in organic matter oxidation coupled to Fe(III) reduction, bioremediation of groundwater contaminated with organics or metals, and electricity production from waste organic matter. In order to investigate G. sulfurreducens central metabolism and electron transport, a metabolic model which integrated genome-based predictions with available genetic and physiological data was developed via the constraint-based modeling approach. Evaluation of the rates of proton production and consumption in the extracellular and cytoplasmic compartments revealed that energy conservation with extracellular electron acceptors, such as Fe(III), was limited relative to that associated with intracellular acceptors. This limitation was attributed to lack of cytoplasmic proton consumption during reduction of extracellular electron acceptors. Model-based analysis of the metabolic cost of producing an extracellular electron shuttle to promote electron transfer to insoluble Fe(III) oxides demonstrated why Geobacter species, which do not produce shuttles, have an energetic advantage over shuttle-producing Fe(III) reducers in subsurface environments. In silico analysis also revealed that the metabolic network of G. sulfurreducens could synthesize amino acids more efficiently than that of Escherichia coli due to the presence of a pyruvate-ferredoxin oxidoreductase, which catalyzes synthesis of pyruvate from acetate and carbon dioxide in a single step. In silico phenotypic analysis of deletion mutants demonstrated the capability of the model to explore the flexibility of G. sulfurreducens central metabolism and correctly predict mutant phenotypes. These results demonstrate that iterative modeling coupled with experimentation can accelerate the understanding of the physiology of poorly studied but environmentally relevant organisms and may help optimize their practical applications.

K-Cone Analysis: Determining All Candidate Values for Kinetic Parameters on a Network Scale

The absence of comprehensive measured kinetic values and the observed inconsistency in the available in vitro kinetic data has hindered the formulation of network-scale kinetic models of biochemical reaction networks. To meet this challenge we present an approach to construct a convex space, termed the k-cone, which contains all the allowable numerical values of the kinetic constants in large-scale biochemical networks. The definition of the k-cone relies on the incorporation of in vivo concentration data and a simplified approach to represent enzyme kinetics within an established constraint-based modeling approach. The k-cone approach was implemented to define the allowable combination of numerical values for a full kinetic model of human red blood cell metabolism and to study its correlated kinetic parameters. The k-cone approach can be used to determine consistency between in vitro measured kinetic values and in vivo concentration and flux measurements when used in a network-scale kinetic model. k-Cone analysis was successful in determining whether in vitro measured kinetic values used in the reconstruction of a kinetic-based model of Saccharomyces cerevisiae central metabolism could reproduce in vivo measurements. Further, the k-cone can be used to determine which numerical values of in vitro measured parameters are required to be changed in a kinetic model if in vivo measured values are not reproduced. k-Cone analysis could identify what minimum number of in vitro determined kinetic parameters needed to be adjusted in the S. cerevisiae model to be consistent with the in vivo data. Applying the k-cone analysis a priori to kinetic model development may reduce the time and effort involved in model building and parameter adjustment. With the recent developments in high-throughput profiling of metabolite concentrations at a whole-cell scale and advances in metabolomics technologies, the k-cone approach presented here may hold the promise for kinetic characterization of metabolic networks as well as other biological functions at a whole-cell level.

The effects of alternate optimal solutions in constraint-based genome-scale metabolic models

Genome-scale constraint-based models of several organisms have now been constructed and are being used for model driven research. A key issue that may arise in the use of such models is the existence of alternate optimal solutions wherein the same maximal objective (e.g., growth rate) can be achieved through different flux distributions. Herein, we investigate the effects that alternate optimal solutions may have on the predicted range of flux values calculated using currently practiced linear (LP) and quadratic programming (QP) methods. An efficient LP-based strategy is described to calculate the range of flux variability that can be present in order to achieve optimal as well as suboptimal objective states. Sample results are provided for growth predictions of E. coli using glucose, acetate, and lactate as carbon substrates. These results demonstrate the extent of flux variability to be highly dependent on environmental conditions and network composition. In addition we examined the impact of alternate optima for growth under gene knockout conditions as calculated using QP-based methods. It was observed that calculations using QP-based methods can show significant variation in growth rate if the flux variability among alternate optima is high. The underlying biological significance and general source of such flux variability is further investigated through the identification of redundancies in the network (equivalent reaction sets) that lead to alternate solutions. Collectively, these results illustrate the variability inherent in metabolic flux distributions and the possible implications of this heterogeneity for constraint-based modeling approaches. These methods also provide an efficient and robust method to calculate the range of flux distributions that can be derived from quantitative fermentation data.

Mechanical systems and assemblies modeling using knowledge-intensive Petri nets formalisms

This paper presents a novel knowledge-based Petri net approach to mechanical systems and assemblies modeling within a design with objects environment. A new unified class of object-oriented knowledge Petri nets, which can incorporate a knowledge-based system with ordinary Petri nets, is defined and used for the unified representations of assembly design and modeling. The object knowledge Petri nets, as a graphical language and a new knowledge-based description scheme, can be used to express the qualitative and quantitative aspects of the assembly design and modeling process in an interactive and integrated way. The four-level hierarchy model is proposed and constructed in terms of function-behaviors, structures, geometries, and features. The function-behavior-structure description is built on more abstract concepts so that it can match well top-down design. The static and dynamic characteristics in the design of assembly can also be captured. With the help of fuzzy logic, the incomplete, imprecise knowledge and uncertainty in the design process can also be dealt with. Therefore, the hybrid design object model can incorporate product data model, top-down design process, and assembly process model using an object-oriented, knowledge-based, feature-based, parametric, and constraint-based modeling approach, and can provide a more accurate and more flexible representation. To verify and demonstrate the effective use of the proposed hybrid design object model, a prototype system has been developed. This research provides a knowledge-intensive framework for intelligent assembly design and modeling.