- Research article
- Open Access
Predicting biological system objectives de novo from internal state measurements
BMC Bioinformatics volume 9, Article number: 43 (2008)
Optimization theory has been applied to complex biological systems to interrogate network properties and develop and refine metabolic engineering strategies. For example, methods are emerging to engineer cells to optimally produce byproducts of commercial value, such as bioethanol, as well as molecular compounds for disease therapy. Flux balance analysis (FBA) is an optimization framework that aids in this interrogation by generating predictions of optimal flux distributions in cellular networks. Critical features of FBA are the definition of a biologically relevant objective function (e.g., maximizing the rate of synthesis of biomass, a unit of measurement of cellular growth) and the subsequent application of linear programming (LP) to identify fluxes through a reaction network. Despite the success of FBA, a central remaining challenge is the definition of a network objective with biological meaning.
We present a novel method called Biological Objective Solution Search (BOSS) for the inference of an objective function of a biological system from its underlying network stoichiometry as well as experimentally-measured state variables. Specifically, BOSS identifies a system objective by defining a putative stoichiometric "objective reaction," adding this reaction to the existing set of stoichiometric constraints arising from known interactions within a network, and maximizing the putative objective reaction via LP, all the while minimizing the difference between the resultant in silico flux distribution and available experimental (e.g., isotopomer) flux data. This new approach allows for discovery of objectives with previously unknown stoichiometry, thus extending the biological relevance from earlier methods. We verify our approach on the well-characterized central metabolic network of Saccharomyces cerevisiae.
We illustrate how BOSS offers insight into the functional organization of biochemical networks, facilitating the interrogation of cellular design principles and development of cellular engineering applications. Furthermore, we describe how growth is the best-fit objective function for the yeast metabolic network given experimentally-measured fluxes.
Systems-based approaches coupled with experimental data have facilitated greater understanding of large-scale biological systems [1, 2]. For example, optimization procedures have recently been used to characterize systemic properties in biology, including phenotypic properties like growth rates and effects of gene knockouts [3–7]. One quantitative measure of a biological phenotype is the set of fluxes through all reactions within a biochemical network . Specifically, flux balance analysis (FBA) is a constraints-based approach that calculates steady-state flux distributions [9–11].
FBA has traditionally been based on the premise that prokaryotes such as Escherichia coli have maximized their growth performance as a response to selective pressure . Consequently, a common objective function in FBA of metabolic networks is the maximization of the rate of synthesis of biomass, a unit of measurement of cellular growth. However, as other types of networks and higher-order systems are interrogated, other objectives may be more accurate in predicting phenotypes. For example, other objective functions that have been previously considered in FBA include optimization of energy production or consumption  and byproduct synthesis . By inferring objective functions of biological systems, cellular design principles may be studied and systems may be exploited for engineering of metabolic byproducts of commercial or medical value [15–20].
In silico frameworks for determining a most-likely objective function have previously been proposed. One such tool, named ObjFind, attempts to identify weightings, termed coefficients of importance (CoIs), on reaction fluxes within a network while minimizing the difference between the resultant flux distribution and known experimental fluxes . In the ObjFind framework, a high CoI indicates a reaction that is more likely a component of the cellular objective function, given available experimental fluxes. However, ObjFind is unable to a priori define objectives, since in FBA the objective function is defined as a single reaction within the system (and represented within the stoichiometric matrix) and not a weighting on multiple reactions (i.e., a set of CoIs). For example, if the true objective reaction has not been experimentally characterized and is not included within the network reconstruction, ObjFind is unable to assign the highest CoI to it and instead chooses an alternate (and consequently suboptimal) reaction or set of reactions as constituting the objective function. Two recent efforts have further attempted to identify the most probable objective of a metabolic system from a set of possible objectives, in one case via a Bayesian-based probability ranking  and in the other case using an Euclidean metric . However, like ObjFind, each of these methods requires that the stoichiometric network reconstruction include the true objective function as an existing reaction in order to yield meaningful predictions.
We present a novel framework, Biological Objective Solution Search (BOSS), for identifying objective functions of biological systems based on the stoichiometry of the underlying biochemical network(s) and known experimental flux data. In this framework, the biological objective function is a de novo reaction (column) that is added to the matrix S representing the stoichiometry of the underlying system. Subsequently, the flux through this particular objective reaction is optimized (maximized) as the objective function using the standard FBA approach. Notably, the objective reaction is not confined to be one of a subset of existing reactions, but rather is allowed to take on any form (e.g., an existing reaction, a combination of existing reactions, or a previously uncharacterized reaction) as specified by an optimization procedure. It is, however, confined to be one linear stoichiometric reaction whose coefficients are determined by the framework.
To illustrate this novelty of BOSS over existing methods (including the three described above), consider the simple system drawn in Figure 1(a). This system is comprised of three components and five reactions. The objective reaction is denoted as robj. Assume that a network reconstruction based on available literature includes all three components and four of the five reactions but excludes the actual system objective, robj, governing the resultant flux distribution. As shown in Figure 1(b), analyzing the network with an approach like ObjFind generates a rank-ordered list of the reactions that best describe the objective function. However, since robj is unknown, the calculated objective will be some combination of weightings on r1 through r4, and the approach will fundamentally fail to capture the actual objective of the system robj. However, as shown in Figure 1(c), by generating a de novo stoichiometric reaction as its objective function, BOSS would in theory be able to recapitulate the actual system objective.
We applied BOSS to the previously reconstructed S. cerevisiae central metabolic network [22, 23] (see Additional file 1 for details of the network) to evaluate which objective reaction it infers for that system given a set of isotopomer flux data [23, 24]. We compared the BOSS-derived objective reaction to the hypothesized system objective of precursor biomass synthesis. We considered two cases, one in which the hypothesized objective reaction (i.e., precursor biomass synthesis) was excluded from the system and another in which it was included within the set of known stoichiometric reactions. We further assessed how BOSS handles noise in experimental flux distributions. Finally, we compared our results to a negative control in which flux distributions comprised of randomly generated values were inputted into BOSS. Ultimately, we illustrate how BOSS extends existing objective function identification tools by inferring the objective reaction of a biochemical network de novo from internal state measurements, and how this tool therefore facilitates future study of cellular design principles and cellular engineering approaches.
The framework that we present for identifying a putative objective function for a given biological system constitutes a single-level optimization problem. Specifically, BOSS integrates network stoichiometry, physico-chemical constraints such as reaction bounds, and experimental flux data to generate a novel stoichiometric reaction corresponding to the most likely objective function of the system. The framework and implementation strategies are described here.
Flux Balance Analysis (FBA)
A key feature in the application of FBA is the optimization of an objective function subject to fundamental constraints on cellular function [9, 11, 25]. FBA requires a stoichiometric network reconstruction, usually represented as a stoichiometric matrix, S, in which components are delineated as rows, component interactions as columns, and stoichiometric relationships as coefficients (see Figure 1 for examples of biochemical networks and the corresponding S matrices). This reconstruction, coupled with reaction fluxes, comprises the principal constraint in FBA, i.e., at steady-state, for each component, the sum of the stoichiometric coefficients multiplied by the corresponding reaction fluxes must equal zero to ensure that mass is balanced within the system.
Commonly used objective functions (i.e., stoichiometric objectives or "objective reactions") in FBA of metabolic networks include biomass production [26, 27], energy production or consumption , and byproduct production . The objective function used in FBA is assumed to be linear by a first-order approximation, since this form simplifies computation and reduces the number of parameters to be defined in the objective experimentally. Changes in a linear objective under varying growth conditions have been observed .
FBA-predicted flux distributions have displayed agreement with experimentally-measured flux data in some cases . Furthermore, FBA has yielded insights into optimal targets for metabolic engineering  and adaptive evolution of E. coli strains , among other applications . Additionally, FBA has successfully predicted metabolic phenotypes of cellular systems by hypothesizing biomass production as the objective function [7, 9, 10, 12, 19, 27, 30–36]. However, other objectives may be equally successful at predicting phenotypes for particular conditions . Furthermore, a metabolic objective may change at different stages of an organism's life cycle . Some studies suggest a greater role for environmental factors in determining cellular objectives than is assumed when FBA is employed with a biomass production objective . Ultimately, as additional network reconstructions are completed and experimental flux data are obtained, there is increasing interest in determining realistic objective functions to explain systemic behavior [18–20, 39].
Formulation of the Objective Function-Finding Algorithm BOSS
The BOSS framework initially takes the form of a bi-level optimization problem that minimizes the sum-squared error between experimentally-measured (in vivo) fluxes and framework-computed (in silico) fluxes ("outer problem") (see Figure 2(a), line 1), subject to the condition that a cellular objective is simultaneously maximized ("inner problem") (lines 2–5). The inner problem takes the canonical form of a FBA problem, wherein an objective reaction introduced to the network is maximized (line 2) subject to thermodynamic (line 3), mass balance (line 4), and uptake (line 5) constraints. The coefficients of the objective reaction are unknown and part of the solution space.
In order to make this bi-level optimization problem computationally tractable, we reformulated it as a single-level optimization problem via the duality theorem of LP, as previously described [3, 40] (see Figure 2(b)). Specifically, the revised form is comprised of an objective that aims to minimize the sum-squared error between experimentally-measured and framework-computed fluxes (line 1) subject to a set of primal (lines 3 through 5) and dual constraints (lines 7 through 9). The objective in line 1 is equivalent to the outer objective of the bi-level optimization problem (see Figure 2(a), line 1), and the primal constraints in lines 3 through 5 are equivalent to the constraints in the original bi-level optimization problem (see Figure 2(a), lines 3 through 5). Two additional constraints are included in the single-level optimization framework. Specifically, the value of the primal and dual problems are set equivalent to one another (line 2), and the flux distribution corresponding to the new objective reaction is normalized to a specific value (line 6). As described in the caption for Figure 2, the decision variables in this optimization are: (1) the stoichiometric coefficients of the objective reaction, Si, objectivewhere i ∈ N, or the set of metabolites, and "objective" denotes the new objective reaction; (2) the framework-computed fluxes vexperimental; (3) the dual variable g associated with the uptake constraint; and (4) the dual variables u indicating shadow prices on the mass balance constraints for each metabolite in the system.
Due to the large number of decision variables, large solution space, existence of multiple local optima, and inherent non-convexity of the problem, we instituted a multiple restart approach, wherein the optimization is run multiple times, each time with different randomly-selected starting values for the decision variables (within a specified range) (see Figure 3(b), item 3). The output of this approach is a series of putative objective reactions, one for each restart. Objective reactions that yield a poor value in the outer optimization (i.e., for which the sum-squared error between experimental and calculated fluxes is relatively high) are removed from the solution set. The remaining putative objective reactions are then clustered into groups based on similarity, and the most populous cluster is chosen as the consensus objective reaction (see Figure 3(b), item 4, as well as Additional file 2). This last step seeks to eliminate solutions to the outer problem that constitute suboptimal local minima. It is assumed that the global minimum to the outer problem (i.e., the best match between BOSS-derived and experimental fluxes, thereby leading to the likely objective reaction) will draw from a larger portion of the initial state space than any local minimum, and thus will gather the greatest proportion of the restarts when solutions are clustered.
This single-level optimization with multiple restarts, followed by clustering and subsequent averaging of the most populous cluster, ultimately yields a stoichiometrically-weighted reaction for which the metabolic network is optimized. This reaction represents a best hypothesis for the network objective function based on available data about the system. Determination of the objective reaction in this manner offers a means to gain insight about network and sub-network objectives and behavior.
Biological System Evaluated: The S. cerevisiae Central Metabolic Network
The previously reconstructed S. cerevisiae central metabolic network [22, 23] was used to assess the BOSS framework (see Additional file 2). This network is a subset of the genome-scale S. cerevisiae metabolic reconstruction [41, 42] and comprises the major carbohydrate metabolism pathways of glycolysis, pentose phosphate, and the citrate cycle, the principal energy metabolism pathway of oxidative phosphorylation, and a precursor biomass synthesis reaction. The network is comprised of 60 metabolites participating in 62 reactions, including six exchange reactions, 55 intracellular reactions, and one precursor biomass synthesis reaction.
A flux distribution for the central metabolic network was previously characterized experimentally [22, 23]. This distribution was obtained by GC-MS tracing of 13C-labeled isotopomers, and corresponds to yeast cells cultivated in batch culture with glucose as the limiting substrate and at a maximum specific growth rate, μmax, of 0.37 h-1. (See  for complete experimental details.) The BOSS framework was applied to this data set.
In addition, the precursor biomass synthesis reaction, which balances the major metabolites from central metabolism that contribute to biomass, was hypothesized to be the "true" objective function of the system for testing the BOSS framework. This hypothesis was based on the underlying premise described previously that organisms have maximized their growth performance through natural selection . The precursor biomass synthesis reaction was previously characterized experimentally [22, 23]. It is important to note that BOSS does not require an assumption of an objective function, and we demonstrate below the success of BOSS at generating a hypothesized objective function de novo.
The general implementation scheme for BOSS is illustrated in Figure 3. Specifically, we inputted the stoichiometric network reconstruction and isotopomer flux data for the S. cerevisiae central metabolic network into BOSS (see Figure 3(a)) and assessed the objective reaction that it derived (see Figure 3(b)). We evaluated two conditions, one in which the hypothesized objective function of precursor biomass synthesis was included in the network reconstruction inputted into BOSS and another in which it was excluded. In the second case, BOSS yielded a zero-vector for the objective reaction, suggesting that the true objective was already a part of the network reconstruction. Therefore, we altered BOSS slightly to identify which of the reactions was the likely objective: we removed, one by one, each reaction flux from the set of experimental flux data (vexperimental in Figure 2), and identified in each case the sum-squared error between the BOSS-derived objective reaction and the reaction corresponding to the removed flux. The reaction that exhibited the smallest sum-squared error through this method was identified as the objective function. This approach was utilized to ensure that, if the actual system objective is already in the stoichiometric network reconstruction, the outer optimization problem of BOSS does not force flux to go through that reaction and rather allows the flux corresponding to the new objective function column to be maximized as part of BOSS.
We considered one additional in silico experiment to evaluate how well our framework handles the kind of variability that is characteristic of actual flux measurements. Specifically, we assessed how a varying amount of noise introduced to the set of experimental flux data affects the objective function derived by BOSS. (See Additional file 3 as well for a discussion of how BOSS handles a paucity of flux data, as applied to a prototypic system.) To perform this test, additional steps were performed during data preparation (see Figure 3(a), item 2). First, each experimental flux for the S. cerevisiae metabolic network was varied randomly via a normal distribution, with standard deviation equal to a set percent of the initial value of the flux. Different percent variances, ranging from five percent to 25 percent in increments of five percent, were considered. Fluxes with an initial value of zero were varied with a standard deviation equal to that of the smallest nonzero flux in the system. BOSS was fed these varied fluxes as "experimental" fluxes. To avoid biasing results toward one particular variant of the "experimental" fluxes, we repeated multiple times this process of randomly varying the experimental flux distribution. In other words, for any given percent variance, we tested multiple noisy flux data sets (n = 20). Consequently, when considering the results of any BOSS run on noisy flux data, we evaluated the coefficients of the objective reaction for each noisy flux data set independently.
Assessing the Results of BOSS
To assess the accuracy and efficacy of our BOSS-derived results, we considered two statistical measures. First, we computed the sum-squared error between the coefficients of the BOSS-derived objective reaction and the coefficients of the hypothesized objective reaction (in this case, precursor biomass synthesis), normalized to the magnitude of the hypothesized objective reaction column vector. This term is abbreviated SSE S throughout the manuscript. Second, we computed the sum-squared error between the reaction fluxes outputted by BOSS and the experimental fluxes inputted to the framework and specified as constraints. This term, which is equivalent to the function being minimized in the BOSS outer optimization problem, is abbreviated SSE v throughout the manuscript. Ideally, smaller values of SSE v should correspond to smaller values of SSE S , as smaller values of SSE S imply more accurate objective reactions.
Simulations of BOSS were run with GAMS v.22.5 on Microsoft Windows XP Professional and Linux (Centos5) machines. The MINOS5 NLP (non-linear programming) solver was utilized. Results were analyzed using code written in MATLAB v.7.5 (part of the MathWorks R2007b release package), with clustering functionality provided by the MATLAB Statistics Toolbox (as detailed in Additional file 2).
Results and Discussion
The S. cerevisiae central metabolic network was used as an experimental system to evaluate the ability of BOSS to identify objective functions of biological systems. Additionally, the network was used to evaluate how well the framework handles noise in flux distributions, which is a characteristic of actual experimental measurements. These results were contrasted with randomly-generated flux data for the S. cerevisiae central metabolic network as a further validation of the BOSS framework. We summarize the results here.
Validation of Objective Function Identification
In inferring the objective function for the S. cerevisiae metabolic network with BOSS, two conditions were evaluated. In the first one, the hypothesized objective reaction (i.e., precursor biomass synthesis) was excluded from the stoichiometric network reconstruction. BOSS identified coefficients approximately equal to that of the precursor biomass synthesis reaction (see Figure 4(a)). The sum-squared error between the BOSS-computed objective reaction and the experimentally-derived precursor biomass synthesis reaction, normalized to the magnitude of the precursor biomass synthesis reaction (SSE S ), was 8.242 × 10-29. The magnitude of this SSE S suggests that BOSS derives an objective reaction equivalent to the precursor biomass synthesis reaction, well within the limits of numerical tolerance.
We also implemented BOSS on the complete S. cerevisiae central metabolic network, i.e., without removing the precursor biomass synthesis reaction. When the precursor biomass synthesis reaction was thus included in the set of stoichiometric reactions, BOSS yielded a zero-vector for the objective reaction (SSE S = 655.0), suggesting that the actual system objective was already a part of the stoichiometric network reconstruction (see Figure 4(b)). Subsequently, we individually removed each of the reaction fluxes from the pool of experimental fluxes defined in the outer optimization problem of BOSS (i.e., vexperimental in Line 1 of Figure 2(a)) and observed the resultant stoichiometric coefficients that BOSS derived for the objective function. SSE S values between the reaction whose flux was removed from vexperimental and the BOSS-derived objective reaction given removal of that flux were compared for each reaction in the system (see Figure 4(c)). The reaction with the smallest SSE S (between the BOSS-derived objective and the experimentally-characterized reaction whose flux was excluded) was the hypothesized system objective of precursor biomass synthesis (SSE S = 4.210 × 10-4). This result further validated the ability of BOSS to identify objective reactions with previously known as well as unknown stoichiometry (see Figures 4(c) and 4(d)). Interestingly, the reaction with the next smallest SSE S (= 42.20) was the ATP maintenance reaction, a reaction that has also been evaluated as an objective of metabolic networks under some conditions [3, 21, 43]. All other reactions in the network had SSE S values that were at least an order of magnitude higher. This key result is also a direct validation of the hypothesis that biomass production is a reasonable objective function for many metabolic networks; with the experimentally-generated flux data, BOSS derived an objective function highly correlated with biomass production although there was no such assumption a priori.
It is important to note that application of FBA to the S. cerevisiae metabolic network does not yield the experimental flux distribution measured by . Because FBA problems have large convex solution spaces, many different flux distributions yield equally optimal results. Therefore, it is significant that BOSS is able to infer the correct coefficients for the precursor biomass synthesis reaction in spite of this difference. Furthermore, although maximization of biomass synthesis is hypothesized as the objective of a metabolic network, the experimental flux data are calculated based on isotopomer measurements and the underlying system stoichiometry but is not biased toward this objective. Thus, the characterization of biomass as the system objective by BOSS is a novel validation of this theoretical assumption.
Accounting for Limitations in Experimental Measurements: Noisy Flux Data
To account for current limitations in experimental (isotopomer) technologies for measuring reaction fluxes, the S. cerevisiae central metabolic network was evaluated by BOSS when different levels of noise were introduced into the experimental flux data. Different percentages of deviation for the flux data were evaluated to determine how quickly the accuracy of the BOSS-computed objective reaction degrades with noisy data.
The results for the S. cerevisiae network when varying levels of noise are introduced into the individual flux values are illustrated in Figure 5. For the purposes of this evaluation, based on the results described above, the hypothesized objective reaction of precursor biomass synthesis was included as part of the underlying network stoichiometry inputted into BOSS, and the set of experimental fluxes consisted of all fluxes except for the flux corresponding to this reaction. Figure 5(a) illustrates the normalized sum-squared error between the BOSS-derived objective reaction and the hypothesized objective function of precursor biomass synthesis (SSE S ) for different percent variances in experimental flux data, ranging from five percent to 25 percent in increments of five percent. (As described above, for each increment, 20 unique "initial" flux distributions, each with the same level of "noise," were computed, and the results for each of these 20 flux distributions were treated independently of the others. Therefore, the plot in Figure 5(a) illustrates the mean SSE S across these 20 distributions, as well as the associated standard deviation within SSES across these 20 distributions.) Up to about 15 percent variance in the flux data, SSE S < 1000, i.e., BOSS is able to generate reconstructions of the hypothesized objective reaction orders of magnitude more accurately than when random flux data are inputted into BOSS (see "Control Case: Random Flux Data" below). As expected, with increasing noise, the results of BOSS exhibit increasing mean SSE S values and the variability of SSE S across different initial flux sets is higher. These data support the utility of our framework in identifying objective functions for biological systems given the current state of experimental flux measurements.
To further highlight BOSS's performance on noisy flux data, snapshots of the objective reactions that it generated at different flux variances are illustrated. Specifically, Figures 5(b), 5(c), and 5(d) correspond to the coefficients identified by BOSS (blue error bars) overlaid on the precursor biomass synthesis reaction (gray bars) for flux variances of 5, 15, and 25 percent, respectively. (Note that, for each flux variance, 20 unique "initial" flux distributions were computed. Consequently, the values for the coefficients derived by BOSS correspond to a representative result for one of these 20 initial flux sets.) For smaller flux variances, the pattern of the coefficients identified by BOSS follows that of the hypothesized network objective of precursor biomass synthesis that was used to generate the noisy flux data (see Figure 5(b)). As expected, the more noisy the flux data, the less precise the objective reaction identified by the framework (see Figure 5(c)).
It is important to note that this analysis involves artificially adding noise to an experimental flux distribution. Although the original experimental fluxes are based on a steady-state metabolic model of yeast central metabolism, it is likely that there is some degree of noise already within the data set as part of the experimental protocol. Consequently, the studies with noise presented here include even higher degrees of noise than the five to 25 percent spectrum evaluated. Nevertheless, the results suggest that BOSS performs well, to an extent, with noisy experimental data.
Control Case: Random Flux Data
To place these results in context, randomly-generated flux values coupled with the underlying network stoichiometry for the S. cerevisiae central metabolic network were inputted into the BOSS framework as a negative control. These fluxes were generated by specifying a normal distribution statistically similar to the characteristics of the experimental flux distribution. Specifically, the randomly-generated fluxes were centered at a mean of 6.697 mmol/g·h with a standard deviation of 9.521 mmol/g·h, identical to the mean and standard deviation of the experimental fluxes reported in [23, 24]. For this test of random flux data, no reaction was removed from the stoichiometric network reconstruction inputted into BOSS to avoid biasing results in any way toward a particular objective. However, to ensure that flux was driven through the objective reaction, thus minimizing SSE v (see discussion above pertaining to Figures 4(b), 4(c), and 4(d)), the flux for the precursor biomass synthesis objective reaction was not included in the pool of experimental flux data inputted into BOSS, as described earlier. A total of 100 randomly-generated flux distributions, each with mean and standard deviation consistent with the experimental flux data, were evaluated, and the SSE S between the hypothesized objective reaction of precursor biomass synthesis and the BOSS-derived objective reaction was 2.076 × 104. In addition, the SSE v was 962.9. These results are orders of magnitude higher than when actual experimental flux data, or even noisy experimental data, are specified (see "Validation of Objective Function Identification" above), thus further confirming the utility of BOSS in inferring objective functions of biological systems with internal state measurements.
We present a novel framework called Biological Objective Solution Search (BOSS) that integrates network stoichiometry and experimental flux data to determine the most likely objective function for a given biological system. We illustrate the utility of BOSS on a model of S. cerevisiae central metabolism with a variety of input conditions, including experimental (isotopomer) flux data, varied experimental flux data meant to mimic noise in experimental measurements, and randomly-generated fluxes. When the underlying network stoichiometry and all reaction fluxes are specified, BOSS identifies an objective reaction with a normalized sum-squared error between the computed objective and the hypothesized objective (i.e., precursor biomass synthesis) (SSE S ) of 4.210 × 10-4. Additionally, when noise is introduced into the flux data to simulate the types of error observed in experimental flux measurements, BOSS identifies the objective reaction with a SSE S of less than 1000 for up to 15 percent experimental flux noise. Thus, we show how BOSS can integrate network stoichiometry and internal state measurements, including in the case of noisy data, to predict the objective function of a biological system.
Interestingly, researchers have theorized that, through evolution, metabolic systems have optimized for biomass production . Here we illustrate that our framework derives this hypothesized objective reaction, precursor biomass synthesis, based on the underlying network stoichiometry and a set of experimentally-measured flux data. This observation validates the BOSS framework and it further confirms the precursor biomass synthesis reaction as an objective of yeast central metabolism.
Current challenges remain in the implementation of the BOSS algorithm. The transformation of a bi-level to a single-level optimization via the strong duality theorem (see "Methods") introduces a high degree of nonlinearity into the problem and non-convexity into the solution space, rendering the BOSS algorithm difficult for nonlinear solvers to tackle. Further, the structure of the inner and outer problems suggests that the solution space for a BOSS problem is not smooth, but rather likely contains discontinuous rifts in the value of the optimization parameter (SSE v ) as one moves smoothly through the state space. These factors complicate the ability of BOSS to determine a global optimum with a small number of restarts, resulting in several hours to days of simulation time depending on the system and input conditions. Perhaps stronger nonlinear solvers, more powerful dedicated computers, and more run-time would yield better predictions of objective functions. These challenges will need to be addressed as BOSS is scaled up to genome-scale applications.
Another hurdle to the meaningful application of BOSS is the paucity of large-scale experimentally-derived flux sets in literature. Isotopomer studies, as well as 13C-constrained fluxomics studies, tend to use simplified models of central metabolism that lump reactions together and thus fail to define a distinct flux for each reaction in the system [44–47]. These simplifications have been necessary due to experimental and computational hurdles, but they limit the degree to which experimental flux data can be used by BOSS. As more powerful fluxomics techniques are developed, more complete flux distribution maps will become available for analysis by BOSS .
BOSS is distinguished from preexisting objective-searching algorithms in its ability to define an objective reaction absent previous knowledge of the stoichiometric structure of the objective function [3, 21]. This attribute could prove useful in analyzing signaling and transcriptional regulatory networks, as their objectives are likely to be more loosely defined and poorly understood than those of metabolic networks [32, 49]. Likewise, studying objectives of multiple metabolically-interacting species or compartmentalized metabolic processes (e.g., mitochondrial metabolism) offers fruitful challenges that cannot be described as simply biomass production. Such analyses could lead to a greater appreciation of the driving evolutionary forces that govern these cellular processes, and could help explain how the whole-cell objective originates. A framework that predicts objectives of a biological system based on known properties of the system (e.g., the stoichiometry of the system and associated reaction fluxes) thus has the potential to address many critical and timely questions in biology.
Stelling J: Mathematical models in microbial systems biology. Curr Opin Microbiol 2004, 7(5):513–518. 10.1016/j.mib.2004.08.004
Klamt S, Stelling J: Two approaches for metabolic pathway analysis? Trends Biotechnol 2003, 21(2):64–69. 10.1016/S0167-7799(02)00034-3
Burgard AP, Maranas CD: Optimization-based framework for inferring and testing hypothesized metabolic objective functions. Biotechnol Bioeng 2003, 82(6):670–677. 10.1002/bit.10617
Heinrich R, Schuster S, Holzhutter HG: Mathematical analysis of enzymic reaction systems using optimization principles. Eur J Biochem 1991, 201(1):1–21. 10.1111/j.1432-1033.1991.tb16251.x
Maynard Smith J: Optimization Theory in Evolution. Ann Rev Ecol Syst 1978, 9: 31–56. 10.1146/annurev.es.09.110178.000335
Deutscher D, Meilijson I, Kupiec M, Ruppin E: Multiple knockout analysis of genetic robustness in the yeast metabolic network. Nat Genet 2006, 38(9):993–998. 10.1038/ng1856
Mahadevan R, Edwards JS, Doyle FJ 3rd: Dynamic flux balance analysis of diauxic growth in Escherichia coli. Biophys J 2002, 83(3):1331–1340.
Sauer U: Metabolic networks in motion: 13C-based flux analysis. Mol Syst Biol 2006, 2: 62. 10.1038/msb4100109
Lee JM, Gianchandani EP, Papin JA: Flux balance analysis in the era of me-tabolomics. Brief Bioinform 2006, 7(2):140–150. 10.1093/bib/bbl007
Kauffman KJ, Prakash P, Edwards JS: Advances in flux balance analysis. Curr Opin Biotechnol 2003, 14(5):491–496. 10.1016/j.copbio.2003.08.001
Price ND, Reed JL, Palsson BO: Genome-scale models of microbial cells: evaluating the consequences of constraints. Nat Rev Microbiol 2004, 2(11):886–897. 10.1038/nrmicro1023
Segre D, Vitkup D, Church GM: Analysis of optimality in natural and perturbed metabolic networks. Proc Natl Acad Sci USA 2002, 99(23):15112–15117. 10.1073/pnas.232349399
Ramakrishna R, Edwards JS, McCulloch A, Palsson BO: Flux-balance analysis of mitochondrial energy metabolism: consequences of systemic stoichiometric constraints. Am J Physiol Regul Integr Comp Physiol 2001, 280(3):R695–704.
Varma A, Boesch BW, Palsson BO: Stoichiometric interpretation of Escherichia coli glucose catabolism under various oxygenation rates. Appl Environ Microbiol 1993, 59(8):2465–2473.
Bro C, Regenberg B, Forster J, Nielsen J: In silico aided metabolic engineering of Saccharomyces cerevisiae for improved bioethanol production. Metab Eng 2006, 8(2):102–111. 10.1016/j.ymben.2005.09.007
Ro DK, Paradise EM, Ouellet M, Fisher KJ, Newman KL, Ndungu JM, Ho KA, Eachus RA, Ham TS, Kirby J, Chang MC, Withers ST, Shiba Y, Sarpong R, Keasling JD: Production of the antimalarial drug precursor artemisinic acid in engineered yeast. Nature 2006, 440(7086):940–943. 10.1038/nature04640
Raman K, Rajagopalan P, Chandra N: Flux balance analysis of mycolic acid pathway: targets for anti-tubercular drugs. PLoS Comput Biol 2005, 1(5):e46. 10.1371/journal.pcbi.0010046
Bailey JE: Toward a science of metabolic engineering. Science 1991, 252(5013):1668–1675. 10.1126/science.2047876
Patil KR, Akesson M, Nielsen J: Use of genome-scale microbial models for metabolic engineering. Curr Opin Biotechnol 2004, 15(1):64–69. 10.1016/j.copbio.2003.11.003
Schuetz R, Kuepfer L, Sauer U: Systematic evaluation of objective functions for predicting intracellular fluxes in Escherichia coli. Mol Syst Biol 2007, 3: 119. 10.1038/msb4100162
Knorr AL, Jain R, Srivastava R: Bayesian-based selection of metabolic objective functions. Bioinformatics 2006.
Wang L, Hatzimanikatis V: Metabolic engineering under uncertainty – II: analysis of yeast metabolism. Metab Eng 2006, 8(2):142–159. 10.1016/j.ymben.2005.11.002
Gombert AK, Moreira dos Santos M, Christensen B, Nielsen J: Network identification and flux quantification in the central metabolism of Saccharomyces cerevisiae under different conditions of glucose repression. J Bacteriol 2001, 183(4):1441–1451. 10.1128/JB.183.4.1441-1451.2001
dos Santos MM, Gombert AK, Christensen B, Olsson L, Nielsen J: Identification of in vivo enzyme activities in the cometabolism of glucose and acetate by Saccharomyces cerevisiae by using 13C-labeled substrates. Eukaryot Cell 2003, 2(3):599–608. 10.1128/EC.2.3.599-608.2003
Papin JA, Palsson BO: The JAK-STAT signaling network in the human B-cell: an extreme signaling pathway analysis. Biophys J 2004, 87(1):37–46. 10.1529/biophysj.103.029884
Varma A, Palsson BO: Stoichiometric flux balance models quantitatively predict growth and metabolic by-product secretion in wild-type Escherichia coli W3110. Appl Environ Microbiol 1994, 60(10):3724–3731.
Edwards JS, Palsson BO: The Escherichia coli MG1655 in silico metabolic genotype: its definition, characteristics, and capabilities. Proc Natl Acad Sci USA 2000, 97(10):5528–5533. 10.1073/pnas.97.10.5528
Fong SS, Burgard AP, Herring CD, Knight EM, Blattner FR, Maranas CD, Pals-son BO: In silico design and adaptive evolution of Escherichia coli for production of lactic acid. Biotechnol Bioeng 2005, 91(5):643–648. 10.1002/bit.20542
Trawick JD, Schilling CH: Use of constraint-based modeling for the prediction and validation of antimicrobial targets. Biochem Pharmacol 2006, 71(7):1026–1035. 10.1016/j.bcp.2005.10.049
Edwards JS, Ibarra RU, Palsson BO: In silico predictions of Escherichia coli metabolic capabilities are consistent with experimental data. Nat Biotechnol 2001, 19(2):125–130. 10.1038/84379
Edwards JS, Palsson BO: Robustness analysis of the Escherichia coli metabolic network. Biotechnol Prog 2000, 16(6):927–939. 10.1021/bp0000712
Klamt S, Stelling J, Ginkel M, Gilles ED: FluxAnalyzer: exploring structure, pathways, and flux distributions in metabolic networks on interactive flux maps. Bioinformatics 2003, 19(2):261–269. 10.1093/bioinformatics/19.2.261
Henry CS, Jankowski MD, Broadbelt LJ, Hatzimanikatis V: Genome-scale thermodynamic analysis of Escherichia coli metabolism. Biophys J 2006, 90(4):1453–1461. 10.1529/biophysj.105.071720
Almaas E, Kovacs B, Vicsek T, Oltvai ZN, Barabasi AL: Global organization of metabolic fluxes in the bacterium Escherichia coli. Nature 2004, 427(6977):839–843. 10.1038/nature02289
Bonarius H, Schmid G, Tramper J: Flux analysis of underdetermined metabolic networks: The quest for the missing constraints. TRENDS BIOTECHNOL 1997, 15(8):308–314. 10.1016/S0167-7799(97)01067-6
Blank LM, Kuepfer L, Sauer U: Large-scale 13C-flux analysis reveals mechanistic principles of metabolic network robustness to null mutations in yeast. Genome Biol 2005, 6(6):R49. 10.1186/gb-2005-6-6-r49
Gupta S: Parasite immune escape: new views into host-parasite interactions. Curr Opin Microbiol 2005, 8(4):428–433. 10.1016/j.mib.2005.06.011
Pfeiffer T, Schuster S: Game-theoretical approaches to studying the evolution of biochemical systems. Trends Biochem Sci 2005, 30(1):20–25. 10.1016/j.tibs.2004.11.006
Stelling J, Klamt S, Bettenbrock K, Schuster S, Gilles ED: Metabolic network structure determines key aspects of functionality and regulation. Nature 2002, 420(6912):190–193. 10.1038/nature01166
Intriligator MD: Mathematical Optimization and Economic Theory. 1st edition. Philadelphia: Society for Industrial and Applied Mathematics; 2002.
Herrgard MJ, Lee BS, Portnoy V, Palsson BO: Integrated analysis of regulatory and metabolic networks reveals novel regulatory mechanisms in Saccharomyces cerevisiae. Genome Res 2006, 16(5):627–635. 10.1101/gr.4083206
Duarte NC, Herrgard MJ, Palsson BO: Reconstruction and validation of Saccharomyces cerevisiae iND750, a fully compartmentalized genome-scale metabolic model. Genome Res 2004, 14(7):1298–1309. 10.1101/gr.2250904
Ibarra RU, Edwards JS, Palsson BO: Escherichia coli K-12 undergoes adaptive evolution to achieve in silico predicted optimal growth. Nature 2002, 420(6912):186–189. 10.1038/nature01149
Emmerling M, Dauner M, Ponti A, Fiaux J, Hochuli M, Szyperski T, Wuthrich K, Bailey JE, Sauer U: Metabolic flux responses to pyruvate kinase knockout in Escherichia coli. J Bacteriol 2002, 184(1):152–164. 10.1128/JB.184.1.152-164.2002
Li M, Ho PY, Yao S, Shimizu K: Effect of lpdA gene knockout on the metabolism in Escherichia coli based on enzyme activities, intracellular metabolite concentrations and metabolic flux analysis by 13C-labeling experiments. J Biotechnol 2006, 122(2):254–266. 10.1016/j.jbiotec.2005.09.016
Zhao J, Baba T, Mori H, Shimizu K: Effect of zwf gene knockout on the metabolism of Escherichia coli grown on glucose or acetate. Metab Eng 2004, 6(2):164–174. 10.1016/j.ymben.2004.02.004
Iwatani S, Van Dien S, Shimbo K, Kubota K, Kageyama N, Iwahata D, Miyano H, Hirayama K, Usuda Y, Shimizu K, Matsui K: Determination of metabolic flux changes during fed-batch cultivation from measurements of intracellular amino acids by LC-MS/MS. J Biotechnol 2007, 128(1):93–111. 10.1016/j.jbiotec.2006.09.004
Reed JHQ, Palsson B: Sensitivity of 13C Isotopomer Based Flux Calculations to Metabolic Network Scope. In Review 2007.
Klamt S, Saez-Rodriguez J, Lindquist JA, Simeoni L, Gilles ED: A methodology for the structural and functional analysis of signaling and regulatory networks. BMC Bioinformatics 2006, 7: 56. 10.1186/1471-2105-7-56
We thank Dr. Kyongbum Lee and Ryan P. Nolan for kindly sharing code with us, as well as Dr. Jong Min Lee, Dr. Kenneth Crowther, Arvind Chavali, Drake Guenther, and Michael Ellis for insightful comments and feedback. Additionally, we thank the National Institutes of Health (NIH) (GM08715/NIH Biotechnology Training Grant to EPG and MAO), the National Science Foundation (CAREER Award to JP), and the Whitaker Foundation for support.
EG and MO conceived of the study and participated in its design and implementation. JP conceived of the study and participated in its design and coordination. AB and CM helped with key analyses of the data. All authors read and approved of the final manuscript.
Erwin P Gianchandani, Matthew A Oberhardt contributed equally to this work.