Identifying constraints that govern cell behavior: a key to converting conceptual to computational models in biology?

Proteome Regulation Patterns Determine Escherichia coli Wild-Type and Mutant Phenotypes

It is generally recognized that proteins constitute the key cellular component in shaping microbial phenotypes. Due to limited cellular resources and space, optimal allocation of proteins is crucial for microbes to facilitate maximum proliferation rates while allowing a flexible response to environmental changes. To account for the growth condition-dependent proteome in the constraint-based metabolic modeling of Escherichia coli, we consolidated a coarse-grained protein allocation approach with the explicit consideration of enzymatic constraints on reaction fluxes. Besides representing physiologically relevant wild-type phenotypes and flux distributions, the resulting protein allocation model (PAM) advances the predictability of the metabolic responses to genetic perturbations. A main driver of mutant phenotypes was ascribed to inherited regulation patterns in protein distribution among metabolic enzymes. Moreover, the PAM correctly reflected metabolic responses to an augmented protein burden imposed by the heterologous expression of green fluorescent protein. In summary, we were able to model the effects of important and frequently applied metabolic engineering approaches on microbial metabolism. Therefore, we want to promote the integration of protein allocation constraints into classical constraint-based models to foster their predictive capabilities and application for strain analysis and engineering purposes.

IMPORTANCE Predictive metabolic models are important, e.g., for generating biological knowledge and designing microbes with superior performance for target compound production. Yet today’s whole-cell models either show insufficient predictive capabilities or are computationally too expensive to be applied to metabolic engineering purposes. By linking the inherent genotype-phenotype relationship to a complete representation of the proteome, the PAM advances the accuracy of simulated phenotypes and intracellular flux distributions of E. coli. Being equally computationally lightweight as classical stoichiometric models and allowing for the application of established in silico tools, the PAM and related simulation approaches will foster the use of a model-driven metabolic research. Applications range from the investigation of mechanisms of microbial evolution to the determination of optimal strain design strategies in metabolic engineering, thus supporting basic scientists and engineers alike.

A Machine Learning Approach for Efficient Selection of Enzyme Concentrations and Its Application for Flux Optimization

The metabolic engineering of pathways has been used extensively to produce molecules of interest on an industrial scale. Methods like gene regulation or substrate channeling helped to improve the desired product yield. Cell-free systems are used to overcome the weaknesses of engineered strains. One of the challenges in a cell-free system is selecting the optimized enzyme concentration for optimal yield. Here, a machine learning approach is used to select the enzyme concentration for the upper part of glycolysis. The artificial neural network approach (ANN) is known to be inefficient in extrapolating predictions outside the box: high predicted values will bump into a sort of “glass ceiling”. In order to explore this “glass ceiling” space, we developed a new methodology named glass ceiling ANN (GC-ANN). Principal component analysis (PCA) and data classification methods are used to derive a rule for a high flux, and ANN to predict the flux through the pathway using the input data of 121 balances of four enzymes in the upper part of glycolysis. The outcomes of this study are i. in silico selection of optimum enzyme concentrations for a maximum flux through the pathway and ii. experimental in vitro validation of the “out-of-the-box” fluxes predicted using this new approach. Surprisingly, flux improvements of up to 63% were obtained. Gratifyingly, these improvements are coupled with a cost decrease of up to 25% for the assay.

Bayesian Integrative Modeling of Genome-Scale Metabolic and Regulatory Networks

The integration of high-throughput data to build predictive computational models of cellular metabolism is a major challenge of systems biology. These models are needed to predict cellular responses to genetic and environmental perturbations. Typically, this response involves both metabolic regulations related to the kinetic properties of enzymes and a genetic regulation affecting their concentrations. Thus, the integration of the transcriptional regulatory information is required to improve the accuracy and predictive ability of metabolic models. Integrative modeling is of primary importance to guide the search for various applications such as discovering novel potential drug targets to develop efficient therapeutic strategies for various diseases. In this paper, we propose an integrative predictive model based on techniques combining semantic web, probabilistic modeling, and constraint-based modeling methods. We applied our approach to human cancer metabolism to predict in silico the growth response of specific cancer cells under approved drug effects. Our method has proven successful in predicting the biomass rates of human liver cancer cells under drug-induced transcriptional perturbations.

Flux prediction using artificial neural network (ANN) for the upper part of glycolysis

The selection of optimal enzyme concentration in multienzyme cascade reactions for the highest product yield in practice is very expensive and time-consuming process. The modelling of biological pathways is a difficult process because of the complexity of the system. The mathematical modelling of the system using an analytical approach depends on the many parameters of enzymes which rely on tedious and expensive experiments. The artificial neural network (ANN) method has been successively applied in different fields of science to perform complex functions. In this study, ANN models were trained to predict the flux for the upper part of glycolysis as inferred by NADH consumption, using four enzyme concentrations i.e., phosphoglucoisomerase, phosphofructokinase, fructose-bisphosphate-aldolase, triose-phosphate-isomerase. Out of three ANN algorithms, the neuralnet package with two activation functions, “logistic” and “tanh” were implemented. The prediction of the flux was very efficient: RMSE and R² were 0.847, 0.93 and 0.804, 0.94 respectively for logistic and tanh functions using a cross validation procedure. This study showed that a systemic approach such as ANN could be used for accurate prediction of the flux through the metabolic pathway. This could help to save a lot of time and costs, particularly from an industrial perspective. The R-code is available at: https://github.com/DSIMB/ANN-Glycolysis-Flux-Prediction.

Temporal constrained objects for modelling neuronal dynamics

BACKGROUND

Several new programming languages and technologies have emerged in the past few decades in order to ease the task of modelling complex systems. Modelling the dynamics of complex systems requires various levels of abstractions and reductive measures in representing the underlying behaviour. This also often requires making a trade-off between how realistic a model should be in order to address the scientific questions of interest and the computational tractability of the model.

METHODS

In this paper, we propose a novel programming paradigm, called temporal constrained objects, which facilitates a principled approach to modelling complex dynamical systems. Temporal constrained objects are an extension of constrained objects with a focus on the analysis and prediction of the dynamic behaviour of a system. The structural aspects of a neuronal system are represented using objects, as in object-oriented languages, while the dynamic behaviour of neurons and synapses are modelled using declarative temporal constraints. Computation in this paradigm is a process of constraint satisfaction within a time-based simulation.

RESULTS

We identified the feasibility and practicality in automatically mapping different kinds of neuron and synapse models to the constraints of temporal constrained objects. Simple neuronal networks were modelled by composing circuit components, implicitly satisfying the internal constraints of each component and interface constraints of the composition. Simulations show that temporal constrained objects provide significant conciseness in the formulation of these models. The underlying computational engine employed here automatically finds the solutions to the problems stated, reducing the code for modelling and simulation control. All examples reported in this paper have been programmed and successfully tested using the prototype language called TCOB. The code along with the programming environment are available at http://github.com/compneuro/TCOB_Neuron.

DISCUSSION

Temporal constrained objects provide powerful capabilities for modelling the structural and dynamic aspects of neural systems. Capabilities of the constraint programming paradigm, such as declarative specification, the ability to express partial information and non-directionality, and capabilities of the object-oriented paradigm especially aggregation and inheritance, make this paradigm the right candidate for complex systems and computational modelling studies. With the advent of multi-core parallel computer architectures and techniques or parallel constraint-solving, the paradigm of temporal constrained objects lends itself to highly efficient execution which is necessary for modelling and simulation of large brain circuits.

Combining inferred regulatory and reconstructed metabolic networks enhances phenotype prediction in yeast

Gene regulatory and metabolic network models have been used successfully in many organisms, but inherent differences between them make networks difficult to integrate. Probabilistic Regulation Of Metabolism (PROM) provides a partial solution, but it does not incorporate network inference and underperforms in eukaryotes. We present an Integrated Deduced And Metabolism (IDREAM) method that combines statistically inferred Environment and Gene Regulatory Influence Network (EGRIN) models with the PROM framework to create enhanced metabolic-regulatory network models. We used IDREAM to predict phenotypes and genetic interactions between transcription factors and genes encoding metabolic activities in the eukaryote, Saccharomyces cerevisiae. IDREAM models contain many fewer interactions than PROM and yet produce significantly more accurate growth predictions. IDREAM consistently outperformed PROM using any of three popular yeast metabolic models and across three experimental growth conditions. Importantly, IDREAM's enhanced accuracy makes it possible to identify subtle synthetic growth defects. With experimental validation, these novel genetic interactions involving the pyruvate dehydrogenase complex suggested a new role for fatty acid-responsive factor Oaf1 in regulating acetyl-CoA production in glucose grown cells.

Application of theoretical methods to increase succinate production in engineered strains

Computational methods have enabled the discovery of non-intuitive strategies to enhance the production of a variety of target molecules. In the case of succinate production, reviews covering the topic have not yet analyzed the impact and future potential that such methods may have. In this work, we review the application of computational methods to the production of succinic acid. We found that while a total of 26 theoretical studies were published between 2002 and 2016, only 10 studies reported the successful experimental implementation of any kind of theoretical knowledge. None of the experimental studies reported an exact application of the computational predictions. However, the combination of computational analysis with complementary strategies, such as directed evolution and comparative genome analysis, serves as a proof of concept and demonstrates that successful metabolic engineering can be guided by rational computational methods.

Predicting network functions with nested patterns

Identifying suitable patterns in complex biological interaction networks helps understanding network functions and allows for predictions at the pattern level: by recognizing a known pattern, one can assign its previously established function. However, current approaches fail for previously unseen patterns, when patterns overlap and when they are embedded into a new network context. Here we show how to conceptually extend pattern-based approaches. We define metabolite patterns in metabolic networks that formalize co-occurrences of metabolites. Our probabilistic framework decodes the implicit information in the networks' metabolite patterns to predict metabolic functions. We demonstrate the predictive power by identifying 'indicator patterns', for instance, for enzyme classification, by predicting directions of novel reactions and of known reactions in new network contexts, and by ranking candidate network extensions for gap filling. Beyond their use in improving genome annotations and metabolic network models, we expect that the concepts transfer to other network types.

Optimal control analysis of the dynamic growth behavior of microorganisms

Understanding the growth behavior of microorganisms using modeling and optimization techniques is an active area of research in the fields of biochemical engineering and systems biology. In this paper, we propose a general modeling framework, based on Monod model, to model the growth of microorganisms. Utilizing the general framework, we formulate an optimal control problem with the objective of maximizing a long-term cellular goal and solve it analytically under various constraints for the growth of microorganisms in a two substrate batch environment. We investigate the relation between long term and short term cellular goals and show that the objective of maximizing cellular concentration at a fixed final time is equivalent to maximization of instantaneous growth rate. We then establish the mathematical connection between the generalized framework and optimal and cybernetic modeling frameworks and derive generalized governing dynamic equations for optimal and cybernetic models. We finally illustrate the influence of various constraints in the cybernetic modeling framework on the optimal growth behavior of microorganisms by solving several dynamic optimization problems using genetic algorithms.

MediaDB: a database of microbial growth conditions in defined media

Isolating pure microbial cultures and cultivating them in the laboratory on defined media is used to more fully characterize the metabolism and physiology of organisms. However, identifying an appropriate growth medium for a novel isolate remains a challenging task. Even organisms with sequenced and annotated genomes can be difficult to grow, despite our ability to build genome-scale metabolic networks that connect genomic data with metabolic function. The scientific literature is scattered with information about defined growth media used successfully for cultivating a wide variety of organisms, but to date there exists no centralized repository to inform efforts to cultivate less characterized organisms by bridging the gap between genomic data and compound composition for growth media. Here we present MediaDB, a manually curated database of defined media that have been used for cultivating organisms with sequenced genomes, with an emphasis on organisms with metabolic network models. The database is accessible online, can be queried by keyword searches or downloaded in its entirety, and can generate exportable individual media formulation files. The data assembled in MediaDB facilitate comparative studies of organism growth media, serve as a starting point for formulating novel growth media, and contribute to formulating media for in silico investigation of metabolic networks. MediaDB is freely available for public use at https://mediadb.systemsbiology.net.

High temperature stimulates acetic acid accumulation and enhances the growth inhibition and ethanol production by Saccharomyces cerevisiae under fermenting conditions

Cellular responses of Saccharomyces cerevisiae to high temperatures of up to 42 °C during ethanol fermentation at a high glucose concentration (i.e., 100 g/L) were investigated. Increased temperature correlated with stimulated glucose uptake to produce not only the thermal protectant glycerol but also ethanol and acetic acid. Carbon flux into the tricarboxylic acid (TCA) cycle correlated positively with cultivation temperature. These results indicate that the increased demand for energy (in the form of ATP), most likely caused by multiple stressors, including heat, acetic acid, and ethanol, was matched by both the fermentation and respiration pathways. Notably, acetic acid production was substantially stimulated compared to that of other metabolites during growth at increased temperature. The acetic acid produced in addition to ethanol seemed to subsequently result in adverse effects, leading to increased production of reactive oxygen species. This, in turn, appeared to cause the specific growth rate, and glucose uptake rate reduced leading to a decrease of the specific ethanol production rate far before glucose depletion. These results suggest that adverse effects from heat, acetic acid, ethanol, and oxidative stressors are synergistic, resulting in a decrease of the specific growth rate and ethanol production rate and, hence, are major determinants of cell stability and ethanol fermentation performance of S. cerevisiae at high temperatures. The results are discussed in the context of possible applications.

A guide to integrating transcriptional regulatory and metabolic networks using PROM (probabilistic regulation of metabolism)

The integration of transcriptional regulatory and metabolic networks is a crucial step in the process of predicting metabolic behaviors that emerge from either genetic or environmental changes. Here, we present a guide to PROM (probabilistic regulation of metabolism), an automated method for the construction and simulation of integrated metabolic and transcriptional regulatory networks that enables large-scale phenotypic predictions for a wide range of model organisms.

Construction of a genome-scale kinetic model of mycobacterium tuberculosis using generic rate equations

The study of biological systems at the genome scale helps us understand fundamental biological processes that govern the activity of living organisms and regulate their interactions with the environment. Genome-scale metabolic models are usually analysed using constraint-based methods, since detailed rate equations and kinetic parameters are often missing. However, constraint-based analysis is limited in capturing the dynamics of cellular processes. In this paper, we present an approach to build a genome-scale kinetic model of Mycobacterium tuberculosis metabolism using generic rate equations. M. tuberculosis causes tuberculosis which remains one of the largest killer infectious diseases. Using a genetic algorithm, we estimated kinetic parameters for a genome-scale metabolic model of M. tuberculosis based on flux distributions derived from Flux Balance Analysis. Our results show that an excellent agreement with flux values is obtained under several growth conditions, although kinetic parameters may vary in different conditions. Parameter variability analysis indicates that a high degree of redundancy remains present in model parameters, which suggests that the integration of other types of high-throughput datasets will enable the development of better constrained models accounting for a variety of in vivo phenotypes.

Ethanol reduces mitochondrial membrane integrity and thereby impacts carbon metabolism of Saccharomyces cerevisiae

Saccharomyces cerevisiae is an excellent ethanol producer, but is rather sensitive to high concentration of ethanol. Here, influences of ethanol on cellular membrane integrity and carbon metabolism of S. cerevisiae were investigated to rationalize mechanism involved in ethanol toxicity. Addition of 5% (v/v) ethanol did neither significantly change the permeability of the cytoplasmic membrane of the reference strain S. cerevisiae BY4741 nor of the ethanol-tolerant strain iETS3. However, the addition of ethanol resulted in a marked decrease in the mitochondrial membrane potential and in increased concentrations of intracellular reactive oxygen species (ROS). The carbon flux was redistributed under these conditions from mainly ethanol production to the TCA cycle. This redistribution was possibly a result of increased energy demand for cell maintenance that increased from about zero to 20-40 mmol ATP (g(CDW)  h)(-1) . This increase in maintenance energy might be explained by the ethanol-induced reduction of the proton motive force and the required removal of ROS. Thus, the stability of the mitochondrial membrane and subsequently the capacity to keep ROS levels low could be important factors to improve tolerance of S. cerevisiae against ethanol.

Improving metabolic flux predictions using absolute gene expression data

Background

Constraint-based analysis of genome-scale metabolic models typically relies upon maximisation of a cellular objective function such as the rate or efficiency of biomass production. Whilst this assumption may be valid in the case of microorganisms growing under certain conditions, it is likely invalid in general, and especially for multicellular organisms, where cellular objectives differ greatly both between and within cell types. Moreover, for the purposes of biotechnological applications, it is normally the flux to a specific metabolite or product that is of interest rather than the rate of production of biomass per se.

Results

An alternative objective function is presented, that is based upon maximising the correlation between experimentally measured absolute gene expression data and predicted internal reaction fluxes. Using quantitative transcriptomics data acquired from Saccharomyces cerevisiae cultures under two growth conditions, the method outperforms traditional approaches for predicting experimentally measured exometabolic flux that are reliant upon maximisation of the rate of biomass production.

Conclusion

Due to its improved prediction of experimentally measured metabolic fluxes, and of its lack of a requirement for knowledge of the biomass composition of the organism under the conditions of interest, the approach is likely to be of rather general utility. The method has been shown to predict fluxes reliably in single cellular systems. Subsequent work will investigate the method’s ability to generate condition- and tissue-specific flux predictions in multicellular organisms.

F2C2: a fast tool for the computation of flux coupling in genome-scale metabolic networks

Background

Flux coupling analysis (FCA) has become a useful tool in the constraint-based analysis of genome-scale metabolic networks. FCA allows detecting dependencies between reaction fluxes of metabolic networks at steady-state. On the one hand, this can help in the curation of reconstructed metabolic networks by verifying whether the coupling between reactions is in agreement with the experimental findings. On the other hand, FCA can aid in defining intervention strategies to knock out target reactions.

Results

We present a new method F2C2 for FCA, which is orders of magnitude faster than previous approaches. As a consequence, FCA of genome-scale metabolic networks can now be performed in a routine manner.

Conclusions

We propose F2C2 as a fast tool for the computation of flux coupling in genome-scale metabolic networks. F2C2 is freely available for non-commercial use at https://sourceforge.net/projects/f2c2/files/.

Modeling metabolic networks for mammalian cell systems: general considerations, modeling strategies, and available tools

Over the past decades, the availability of large amounts of information regarding cellular processes and reaction rates, along with increasing knowledge about the complex mechanisms involved in these processes, has changed the way we approach the understanding of cellular processes. We can no longer rely only on our intuition for interpreting experimental data and evaluating new hypotheses, as the information to analyze is becoming increasingly complex. The paradigm for the analysis of cellular systems has shifted from a focus on individual processes to comprehensive global mathematical descriptions that consider the interactions of metabolic, genomic, and signaling networks. Analysis and simulations are used to test our knowledge by refuting or validating new hypotheses regarding a complex system, which can result in predictive capabilities that lead to better experimental design. Different types of models can be used for this purpose, depending on the type and amount of information available for the specific system. Stoichiometric models are based on the metabolic structure of the system and allow explorations of steady state distributions in the network. Detailed kinetic models provide a description of the dynamics of the system, they involve a large number of reactions with varied kinetic characteristics and require a large number of parameters. Models based on statistical information provide a description of the system without information regarding structure and interactions of the networks involved. The development of detailed models for mammalian cell metabolism has only recently started to grow more strongly, due to the intrinsic complexities of mammalian systems, and the limited availability of experimental information and adequate modeling tools. In this work we review the strategies, tools, current advances, and recent models of mammalian cells, focusing mainly on metabolism, but discussing the methodology applied to other types of networks as well.

Robustness of metabolic networks: a review of existing definitions

Describing the determinants of robustness of biological systems has become one of the central questions in systems biology. Despite the increasing research efforts, it has proven difficult to arrive at a unifying definition for this important concept. We argue that this is due to the multifaceted nature of the concept of robustness and the possibility to formally capture it at different levels of systemic formalisms (e.g., topology and dynamic behavior). Here we provide a comprehensive review of the existing definitions of robustness pertaining to metabolic networks. As kinetic approaches have been excellently reviewed elsewhere, we focus on definitions of robustness proposed within graph-theoretic and constraint-based formalisms.

Systems biology of the qa gene cluster in Neurospora crassa

An ensemble of genetic networks that describe how the model fungal system, Neurospora crassa, utilizes quinic acid (QA) as a sole carbon source has been identified previously. A genetic network for QA metabolism involves the genes, qa-1F and qa-1S, that encode a transcriptional activator and repressor, respectively and structural genes, qa-2, qa-3, qa-4, qa-x, and qa-y. By a series of 4 separate and independent, model-guided, microarray experiments a total of 50 genes are identified as QA-responsive and hypothesized to be under QA-1F control and/or the control of a second QA-responsive transcription factor (NCU03643) both in the fungal binuclear Zn(II)2Cys6 cluster family. QA-1F regulation is not sufficient to explain the quantitative variation in expression profiles of the 50 QA-responsive genes. QA-responsive genes include genes with products in 8 mutually connected metabolic pathways with 7 of them one step removed from the tricarboxylic (TCA) Cycle and with 7 of them one step removed from glycolysis: (1) starch and sucrose metabolism; (2) glycolysis/glucanogenesis; (3) TCA Cycle; (4) butanoate metabolism; (5) pyruvate metabolism; (6) aromatic amino acid and QA metabolism; (7) valine, leucine, and isoleucine degradation; and (8) transport of sugars and amino acids. Gene products both in aromatic amino acid and QA metabolism and transport show an immediate response to shift to QA, while genes with products in the remaining 7 metabolic modules generally show a delayed response to shift to QA. The additional QA-responsive cutinase transcription factor-1β (NCU03643) is found to have a delayed response to shift to QA. The series of microarray experiments are used to expand the previously identified genetic network describing the qa gene cluster to include all 50 QA-responsive genes including the second transcription factor (NCU03643). These studies illustrate new methodologies from systems biology to guide model-driven discoveries about a core metabolic network involving carbon and amino acid metabolism in N. crassa.

Why does yeast ferment? A flux balance analysis study

Advances in biological techniques have led to the availability of genome-scale metabolic reconstructions for yeast. The size and complexity of such networks impose limits on what types of analyses one can perform. Constraint-based modelling overcomes some of these restrictions by using physicochemical constraints to describe the potential behaviour of an organism. FBA (flux balance analysis) highlights flux patterns through a network that serves to achieve a particular objective and requires a minimal amount of data to make quantitative inferences about network behaviour. Even though FBA is a powerful tool for system predictions, its general formulation sometimes results in unrealistic flux patterns. A typical example is fermentation in yeast: ethanol is produced during aerobic growth in excess glucose, but this pattern is not present in a typical FBA solution. In the present paper, we examine the issue of yeast fermentation against respiration during growth. We have studied a number of hypotheses from the modelling perspective, and novel formulations of the FBA approach have been tested. By making the observation that more respiration requires the synthesis of more mitochondria, an energy cost related to mitochondrial synthesis is added to the FBA formulation. Results, although still approximate, are closer to experimental observations than earlier FBA analyses, at least on the issue of fermentation.

Towards a genome-scale kinetic model of cellular metabolism

Background

Advances in bioinformatic techniques and analyses have led to the availability of genome-scale metabolic reconstructions. The size and complexity of such networks often means that their potential behaviour can only be analysed with constraint-based methods. Whilst requiring minimal experimental data, such methods are unable to give insight into cellular substrate concentrations. Instead, the long-term goal of systems biology is to use kinetic modelling to characterize fully the mechanics of each enzymatic reaction, and to combine such knowledge to predict system behaviour.

Results

We describe a method for building a parameterized genome-scale kinetic model of a metabolic network. Simplified linlog kinetics are used and the parameters are extracted from a kinetic model repository. We demonstrate our methodology by applying it to yeast metabolism. The resultant model has 956 metabolic reactions involving 820 metabolites, and, whilst approximative, has considerably broader remit than any existing models of its type. Control analysis is used to identify key steps within the system.

Conclusions

Our modelling framework may be considered a stepping-stone toward the long-term goal of a fully-parameterized model of yeast metabolism. The model is available in SBML format from the BioModels database (BioModels ID: MODEL1001200000) and at http://www.mcisb.org/resources/genomescale/.

Metabolic flux analysis in plants: coping with complexity

Theory and experience in metabolic engineering both show that metabolism operates at the network level. In plants, this complexity is compounded by a high degree of compartmentation and the synthesis of a very wide array of secondary metabolic products. A further challenge to understanding and predicting plant metabolic function is posed by our ignorance about the structure of metabolic networks even in well-studied systems. Metabolic flux analysis (MFA) provides tools to measure and model the functioning of metabolism, and is making significant contributions to coping with their complexity. This review gives an overview of different MFA approaches, the measurements required to implement them and the information they yield. The application of MFA methods to plant systems is then illustrated by several examples from the recent literature. Next, the challenges that plant metabolism poses for MFA are discussed together with ways that these can be addressed. Lastly, new developments in MFA are described that can be expected to improve the range and reliability of plant MFA in the coming years.

Quantifying epistatic interactions among the components constituting the protein translation system

In principle, the accumulation of knowledge regarding the molecular basis of biological systems should allow the development of large-scale kinetic models of their functions. However, the development of such models requires vast numbers of parameters, which are difficult to obtain in practice. Here, we used an in vitro translation system, consisting of 69 defined components, to quantify the epistatic interactions among changes in component concentrations through Bahadur expansion, thereby obtaining a coarse-grained model of protein synthesis activity. Analyses of the data measured using various combinations of component concentrations indicated that the contributions of larger than 2-body inter-component epistatic interactions are negligible, despite the presence of larger than 2-body physical interactions. These findings allowed the prediction of protein synthesis activity at various combinations of component concentrations from a small number of samples, the principle of which is applicable to analysis and optimization of other biological systems. Moreover, the average ratio of 2- to 1-body terms was estimated to be as small as 0.1, implying high adaptability and evolvability of the protein translation system.

Flux balance analysis: a geometric perspective

Advances in the field of bioinformatics have led to reconstruction of genome-scale networks for a number of key organisms. The application of physicochemical constraints to these stoichiometric networks allows researchers, through methods such as flux balance analysis, to highlight key sets of reactions necessary to achieve particular objectives. The key benefits of constraint-based analysis lie in the minimal knowledge required to infer systemic properties. However, network degeneracy leads to a large number of flux distributions that satisfy any objective; moreover, these distributions may be dominated by biologically irrelevant internal cycles. By examining the geometry underlying the problem, we define two methods for finding a unique solution within the space of all possible flux distributions; such a solution contains no internal cycles, and is representative of the space as a whole. The first method draws on typical geometric knowledge, but cannot be applied to large networks because of the high computational complexity of the problem. Thus a second method, an iteration of linear programs which scales easily to the genome scale, is defined. The algorithm is run on four recent genome-scale models, and unique flux solutions are found. The algorithm set out here will allow researchers in flux balance analysis to exchange typical solutions to their models in a reproducible format. Moreover, having found a single solution, statistical analyses such as correlations may be performed.

Flux balance analysis: interrogating genome-scale metabolic networks

Flux balance analysis (FBA) is a computational method to analyze reconstructions of biochemical networks. FBA requires the formulation of a biochemical network in a precise mathematical framework called a stoichiometric matrix. An objective function is defined (e.g., growth rate) toward which the system is assumed to be optimized. In this chapter, we present the methodology, theory, and common pitfalls of the application of FBA.

Genome-scale model for Clostridium acetobutylicum: Part II. Development of specific proton flux states and numerically determined sub-systems

A regulated genome-scale model for Clostridium acetobutylicum ATCC 824 was developed based on its metabolic network reconstruction. To aid model convergence and limit the number of flux-vector possible solutions (the size of the phenotypic solution space), modeling strategies were developed to impose a new type of constraint at the endo-exo-metabolome interface. This constraint is termed the specific proton flux state, and its use enabled accurate prediction of the extracellular medium pH during vegetative growth of batch cultures. The specific proton flux refers to the influx or efflux of free protons (per unit biomass) across the cell membrane. A specific proton flux state encompasses a defined range of specific proton fluxes and includes all metabolic flux distributions resulting in a specific proton flux within this range. Effective simulation of time-course batch fermentation required the use of independent flux balance solutions from an optimum set of specific proton flux states. Using a real-coded genetic algorithm to optimize temporal bounds of specific proton flux states, we show that six separate specific proton flux states are required to model vegetative-growth metabolism and accurately predict the extracellular medium pH. Further, we define the apparent proton flux stoichiometry per weak acids efflux and show that this value decreases from approximately 3.5 mol of protons secreted per mole of weak acids at the start of the culture to approximately 0 at the end of vegetative growth. Calculations revealed that when specific weak acids production is maximized in vegetative growth, the net proton exchange between the cell and environment occurs primarily through weak acids efflux (apparent proton flux stoichiometry is 1). However, proton efflux through cation channels during the early stages of acidogenesis was found to be significant. We have also developed the concept of numerically determined sub-systems of genome-scale metabolic networks here as a sub-network with a one-dimensional null space basis set. A numerically determined sub-system was constructed in the genome-scale metabolic network to study the flux magnitudes and directions of acetylornithine transaminase, alanine racemase, and D-alanine transaminase. These results were then used to establish additional constraints for the genome-scale model.

Dynamic analysis of integrated signaling, metabolic, and regulatory networks

Extracellular cues affect signaling, metabolic, and regulatory processes to elicit cellular responses. Although intracellular signaling, metabolic, and regulatory networks are highly integrated, previous analyses have largely focused on independent processes (e.g., metabolism) without considering the interplay that exists among them. However, there is evidence that many diseases arise from multifunctional components with roles throughout signaling, metabolic, and regulatory networks. Therefore, in this study, we propose a flux balance analysis (FBA)–based strategy, referred to as integrated dynamic FBA (idFBA), that dynamically simulates cellular phenotypes arising from integrated networks. The idFBA framework requires an integrated stoichiometric reconstruction of signaling, metabolic, and regulatory processes. It assumes quasi-steady-state conditions for “fast” reactions and incorporates “slow” reactions into the stoichiometric formalism in a time-delayed manner. To assess the efficacy of idFBA, we developed a prototypic integrated system comprising signaling, metabolic, and regulatory processes with network features characteristic of actual systems and incorporating kinetic parameters based on typical time scales observed in literature. idFBA was applied to the prototypic system, which was evaluated for different environments and gene regulatory rules. In addition, we applied the idFBA framework in a similar manner to a representative module of the single-cell eukaryotic organism Saccharomyces cerevisiae. Ultimately, idFBA facilitated quantitative, dynamic analysis of systemic effects of extracellular cues on cellular phenotypes and generated comparable time-course predictions when contrasted with an equivalent kinetic model. Since idFBA solves a linear programming problem and does not require an exhaustive list of detailed kinetic parameters, it may be efficiently scaled to integrated intracellular systems that incorporate signaling, metabolic, and regulatory processes at the genome scale, such as the S. cerevisiae system presented here.

Cellular systems comprise many diverse components and component interactions spanning signal transduction, transcriptional regulation, and metabolism. Although signaling, metabolic, and regulatory activities are often investigated independently of one another, there is growing evidence that considerable interplay occurs among them, and that the malfunctioning of this interplay is associated with disease. The computational analysis of integrated networks has been challenging because of the varying time scales involved as well as the sheer magnitude of such systems (e.g., the numbers of rate constants involved). To this end, we developed a novel computational framework called integrated dynamic flux balance analysis (idFBA) that generates quantitative, dynamic predictions of species concentrations spanning signaling, regulatory, and metabolic processes. idFBA extends an existing approach called flux balance analysis (FBA) in that it couples “fast” and “slow” reactions, thereby facilitating the study of whole-cell phenotypes and not just sub-cellular network properties. We applied this framework to a prototypic integrated system derived from literature as well as a representative integrated yeast module (the high-osmolarity glycerol [HOG] pathway) and generated time-course predictions that matched with available experimental data. By extending this framework to larger-scale systems, phenotypic profiles of whole-cell systems could be attained expeditiously.

Individual-based modelling: an essential tool for microbiology

Micro-organisms play a central role in every ecosystem and in the global biomass cycle. They are strongly involved in many fields of human interest, from medicine to the food industry and waste control. Nevertheless, most micro-organisms remain almost unknown, and nearly 99% of them have not yet been successfully cultured in vitro. Therefore, new approaches and new tools must be developed in order to understand the collective behaviour of microbial communities in any natural or artificial setting. In particular, theoretical and practical methodologies to deal with such systems at a mesoscopic level of description (covering the range from 100 to 10(8) cells) are required. Individual-based modelling (IBM) has become a widely used tool for describing complex systems made up of autonomous entities, such as ecosystems and social networks. Individual-based models (IBMs) provide some advantages over the traditional whole-population models: (a) they are bottom-up approaches, so they describe the behaviour of a system as a whole by establishing procedural rules for the individuals and for their interactions, and thus allow more realistic assumptions for the model of the individuals than population models do; (b) they permit the introduction of randomness and individual variability, so they can reproduce the diversity found in real systems; and (c) they can account for individual adaptive behaviour to their environmental conditions, so the evolution of the whole system arises from the dynamics that govern individuals in their pursuit of optimal fitness. However, they also present some drawbacks: they lack the clarity of continuous models and may easily become rambling, which makes them difficult to analyse and communicate. All in all, IBMs supply a holistic description of microbial systems and their emerging properties. They are specifically appropriate to deal with microbial communities in non-steady states, and spatially explicit IBMs are particularly appropriate to study laboratory and natural microbiological systems with spatial heterogeneity. In this paper, we review IBM methodology applied to microbiology. We also present some results obtained from the application of Individual Discrete Simulations, an IBM of ours, to the study of bacterial communities, yeast cultures and Plasmodium falciparum-infected erythrocytes in vitro cultures of Plasmodium falciparum-infected erythrocytes.

Impact of the solvent capacity constraint on E. coli metabolism

Background

Obtaining quantitative predictions for cellular metabolic activities requires the identification and modeling of the physicochemical constraints that are relevant at physiological growth conditions. Molecular crowding in a cell's cytoplasm is one such potential constraint, as it limits the solvent capacity available to metabolic enzymes.

Results

Using a recently introduced flux balance modeling framework (FBAwMC) here we demonstrate that this constraint determines a metabolic switch in E. coli cells when they are shifted from low to high growth rates. The switch is characterized by a change in effective optimization strategy, the excretion of acetate at high growth rates, and a global reorganization of E. coli metabolic fluxes, the latter being partially confirmed by flux measurements of central metabolic reactions.

Conclusion

These results implicate the solvent capacity as an important physiological constraint acting on E. coli cells operating at high metabolic rates and for the activation of a metabolic switch when they are shifted from low to high growth rates. The relevance of this constraint in the context of both the aerobic ethanol excretion seen in fast growing yeast cells (Crabtree effect) and the aerobic glycolysis observed in rapidly dividing cancer cells (Warburg effect) should be addressed in the future.

Predicting gene essentiality using genome-scale in silico models

Genome-scale metabolic models of organisms can be reconstructed using annotated genome sequence information, well-curated databases, and primary research literature. The metabolic reaction stoichiometry and other physicochemical factors are incorporated into the model, thus imposing constraints that represent restrictions on phenotypic behavior. Based on this premise, the theoretical capabilities of the metabolic network can be assessed by using a mathematical technique known as flux balance analysis (FBA). This modeling framework, also known as the constraint-based reconstruction and analysis approach, differs from other modeling strategies because it does not attempt to predict exact network behavior. Instead, this approach uses known constraints to separate the states that a system can achieve from those that it cannot. In recent years, this strategy has been employed to probe the metabolic capabilities of a number of organisms, to generate and test experimental hypotheses, and to predict accurately metabolic phenotypes and evolutionary outcomes. This chapter introduces the constraint-based modeling approach and focuses on its application to computationally predicting gene essentiality.

Metabolic reconstruction and analysis for parasite genomes

With the completion of sequencing projects for several parasite genomes, efforts are ongoing to make sense of this mass of information in terms of the gene products encoded and their interactions in the growth, development and survival of parasites. The emerging science of systems biology aims to explain the complex relationship between genotype and phenotype by using network models. One area in which this approach has been particularly successful is in the modeling of metabolism. With an accurate picture of the set of metabolic reactions encoded in a genome, it is now possible to identify enzymes or transporters that might be viable targets for new drugs. Because these predictions greatly depend on the quality and completeness of the genome annotation, there are substantial efforts in the scientific community to increase the numbers of metabolic enzymes identified. In this review, we discuss the opportunities for using metabolic reconstruction and analysis tools in parasitology research, and their applications to protozoan parasites.

Something from nothing: bridging the gap between constraint-based and kinetic modelling

Two divergent modelling methodologies have been adopted to increase our understanding of metabolism and its regulation. Constraint-based modelling highlights the optimal path through a stoichiometric network within certain physicochemical constraints. Such an approach requires minimal biological data to make quantitative inferences about network behaviour; however, constraint-based modelling is unable to give an insight into cellular substrate concentrations. In contrast, kinetic modelling aims to characterize fully the mechanics of each enzymatic reaction. This approach suffers because parameterizing mechanistic models is both costly and time-consuming. In this paper, we outline a method for developing a kinetic model for a metabolic network, based solely on the knowledge of reaction stoichiometries. Fluxes through the system, estimated by flux balance analysis, are allowed to vary dynamically according to linlog kinetics. Elasticities are estimated from stoichiometric considerations. When compared to a popular branched model of yeast glycolysis, we observe an excellent agreement between the real and approximate models, despite the absence of (and indeed the requirement for) experimental data for kinetic constants. Moreover, using this particular methodology affords us analytical forms for steady state determination, stability analyses and studies of dynamical behaviour.

Toward whole cell modeling and simulation: comprehensive functional genomics through the constraint-based approach

The increasing availability of various system-level, or so-called 'omics', datasets, in concert with existing data from the primary research literature, is facilitating the development of genome-scale metabolic models for many organisms. By incorporating the metabolic reaction stoichiometry as well as other physicochemical properties into systemic network reconstructions, these models account for the constraints that restrict an organism's phenotypic behavior. Accordingly, unlike many contemporary modeling strategies, this constraint-based modeling approach does not attempt to predict network behavior exactly; rather, it seeks to clearly distinguish those network states that a system can achieve from those that it cannot. A variety of analytical tools have been designed and developed to probe these models, thus enabling studies that investigate the metabolic capabilities of a number of organisms, that generate and test experimental hypotheses, and that predict accurately metabolic phenotypes and evolutionary outcomes. This chapter introduces the concepts that underlie the constraint-based modeling approach, and describes several of its applications with an emphasis on those potentially relevant to the drug development field. In addition, while this chapter focuses on the primary application of the constraint-based approach to date, namely in modeling metabolic networks, the latter sections of the chapter discuss its relatively recent application to modeling other cellular systems. Finally, the chapter concludes with an assessment of future directions focusing on the efforts that will be required to utilize the constraint-based approach in generating a holistic model of a viable organism.

Genome-scale thermodynamic analysis of Escherichia coli metabolism

Genome-scale metabolic models are an invaluable tool for analyzing metabolic systems as they provide a more complete picture of the processes of metabolism. We have constructed a genome-scale metabolic model of Escherichia coli based on the iJR904 model developed by the Palsson Laboratory at the University of California at San Diego. Group contribution methods were utilized to estimate the standard Gibbs free energy change of every reaction in the constructed model. Reactions in the model were classified based on the activity of the reactions during optimal growth on glucose in aerobic media. The most thermodynamically unfavorable reactions involved in the production of biomass in E. coli were identified as ATP phosphoribosyltransferase, ATP synthase, methylene-tetra-hydrofolate dehydrogenase, and tryptophanase. The effect of a knockout of these reactions on the production of biomass and the production of individual biomass precursors was analyzed. Changes in the distribution of fluxes in the cell after knockout of these unfavorable reactions were also studied. The methodologies and results discussed can be used to facilitate the refinement of the feasible ranges for cellular parameters such as species concentrations and reaction rate constants.

Uniform sampling of steady-state flux spaces: means to design experiments and to interpret enzymopathies

Reconstruction of genome-scale metabolic networks is now possible using multiple different data types. Constraint-based modeling is an approach to interrogate capabilities of reconstructed networks by constraining possible cellular behavior through the imposition of physicochemical laws. As a result, a steady-state flux space is defined that contains all possible functional states of the network. Uniform random sampling of the steady-state flux space allows for the unbiased appraisal of its contents. Monte Carlo sampling of the steady-state flux space of the reconstructed human red blood cell metabolic network under simulated physiologic conditions yielded the following key results: 1), probability distributions for the values of individual metabolic fluxes showed a wide variety of shapes that could not have been inferred without computation; 2), pairwise correlation coefficients were calculated between all fluxes, determining the level of independence between the measurement of any two fluxes, and identifying highly correlated reaction sets; and 3), the network-wide effects of the change in one (or a few) variables (i.e., a simulated enzymopathy or fixing a flux range based on measurements) were computed. Mathematical models provide the most compact and informative representation of a hypothesis of how a cell works. Thus, understanding model predictions clearly is vital to driving forward the iterative model-building procedure that is at the heart of systems biology. Taken together, the Monte Carlo sampling procedure provides a broadening of the constraint-based approach by allowing for the unbiased and detailed assessment of the impact of the applied physicochemical constraints on a reconstructed network.

Genomic and mechanistic insights into the biodegradation of organic pollutants

Several new methodologies have enabled recent studies on the microbial biodegradation mechanisms of organic pollutants. Culture-independent techniques for analysis of the genetic and metabolic potential of natural and model microbial communities that degrade organic pollutants have identified new metabolic pathways and enzymes for aerobic and anaerobic degradation. Furthermore, structural studies of the enzymes involved have revealed the specificities and activities of key catabolic enzymes, such as dioxygenases. Genome sequencing of several biodegradation-relevant microorganisms have provided the first whole-genome insights into the genetic background of the metabolic capability and biodegradation versatility of these organisms. Systems biology approaches are still in their infancy, but are becoming increasingly helpful to unravel, predict and quantify metabolic abilities within particular organisms or microbial consortia.

Genome-scale models of microbial cells: evaluating the consequences of constraints

Microbial cells operate under governing constraints that limit their range of possible functions. With the availability of annotated genome sequences, it has become possible to reconstruct genome-scale biochemical reaction networks for microorganisms. The imposition of governing constraints on a reconstructed biochemical network leads to the definition of achievable cellular functions. In recent years, a substantial and growing toolbox of computational analysis methods has been developed to study the characteristics and capabilities of microorganisms using a constraint-based reconstruction and analysis (COBRA) approach. This approach provides a biochemically and genetically consistent framework for the generation of hypotheses and the testing of functions of microbial cells.