Genome-scale modeling using flux ratio constraints to enable metabolic engineering of clostridial metabolism in silico
© McAnulty et al.; licensee BioMed Central Ltd. 2012
Received: 8 January 2012
Accepted: 14 May 2012
Published: 14 May 2012
Genome-scale metabolic networks and flux models are an effective platform for linking an organism genotype to its phenotype. However, few modeling approaches offer predictive capabilities to evaluate potential metabolic engineering strategies in silico.
A new method called “f lux b alance a nalysis with flux ratio s (FBrAtio)” was developed in this research and applied to a new genome-scale model of Clostridium acetobutylicum ATCC 824 (i CAC490) that contains 707 metabolites and 794 reactions. FBrAtio was used to model wild-type metabolism and metabolically engineered strains of C. acetobutylicum where only flux ratio constraints and thermodynamic reversibility of reactions were required. The FBrAtio approach allowed solutions to be found through standard linear programming. Five flux ratio constraints were required to achieve a qualitative picture of wild-type metabolism for C. acetobutylicum for the production of: (i) acetate, (ii) lactate, (iii) butyrate, (iv) acetone, (v) butanol, (vi) ethanol, (vii) CO2 and (viii) H2. Results of this simulation study coincide with published experimental results and show the knockdown of the acetoacetyl-CoA transferase increases butanol to acetone selectivity, while the simultaneous over-expression of the aldehyde/alcohol dehydrogenase greatly increases ethanol production.
FBrAtio is a promising new method for constraining genome-scale models using internal flux ratios. The method was effective for modeling wild-type and engineered strains of C. acetobutylicum.
KeywordsGenome-scale model clostridia flux ratio flux balance analysis metabolic engineering systems biology
Modeling clostridial metabolism
Butanol is of considerable research interest as a potential biofuel, and its renewable production through fermentation is sought largely from the clostridia. In particular, Clostridium acetobutylicum ATCC 824 has been one of multiple clostridia researched for butanol production over the past few decades. In fact, the first applications of metabolic flux balancing were performed using a model of C. acetobutylicum primary metabolism to understand what caused this organism to produce butanol and the competing metabolic byproducts: (i) acetate, (ii) butyrate, (iii) lactate, (iv) acetone, (v) ethanol, and several others in small amounts [1, 2]. Flux modeling of the primary metabolism of C. acetobutylicum has led to a better understanding of the role cofactor balancing plays in directing global metabolic changes. It has played a significant role in metabolic engineering by identifying bottlenecks and critical flux distributions at metabolic branch points [3–8]. Multiple “genome-scale” metabolic network reconstructions now exist for C. acetobutylicum[9–12]. Similar networks and their corresponding genome-scale models have been reviewed extensively [10, 13–20]. In general, they are used to (i) complete genome annotation , (ii) predict optimal culturing conditions [22–24], (iii) discover genomic regulation [20, 25, 26], (iv) identify essential genes and drug targets [27–33], (v) study strain evolution [34, 35], and (vi) design productive strains [36–38]. Modeling results on the genome-scale have been applied to both “acidogenic” and “solventogenic” programs of clostridial metabolism . The acidogenic program is characterized by high acids (i.e., acetate and butyrate) production and high growth rates, and the solventogenic program (i.e., acetone and butanol production) largely coincides with the stationary growth phase of the culture. During solventogenesis, acetate and butyrate are re-consumed by the culture and converted to acetone and butanol. The genetic program of this metabolic shift between acids and solvents production has been studied in detail . Several insights into C. acetobutylicum metabolism have been gained from “gap filling” the metabolic network by locating previously unknown enzymes and biochemical reactions [11, 12]. The total rate at which a cell produces/consumes protons through the several membrane transport mechanisms is termed the specific proton flux (SPF), and this parameter has shown to significantly reduce the total number of flux “solutions” available for the under-determined genome-scale model of C. acetobutylicum. Reducing the number of degrees of freedom of these genome-scale models through application of genetic regulation and physicochemical constraints has been recognized as a key strategy for generating metabolic flux predictions that coincide with experimental observations .
Engineering clostridial metabolism
Metabolic engineering in silico
The goal of metabolic engineering in silico is to derive (or at least evaluate) potential metabolic engineering strategies prior to constructing them in the laboratory. For example, will a particular gene over-expression or knockout in C. acetobutylicum increase butanol production? Answering questions of this type is one of the potential uses of genome-scale modeling. However, with the initial genome-scale model for C. acetobutylicum[11, 12], these questions could not be addressed without constraints on acid/solvent production. These constraints artificially specified ranges for secretion rates of acid and solvent products. These were necessary due to the large number of degrees of freedom that exist in the under-determined genome-scale model and the high degree of branching in the primary metabolism of clostridia. Simply, too many flux solutions were available if the user was only to define the substrate uptake rate and a proper objective function. The production of products/byproducts by a metabolic network not only completes elemental balances but it also regenerates and balances cofactors. In clostridial metabolism, ATP is regenerated by the production of acetate or butyrate, and NAD+ is produced by the production of (i) lactate, (ii) ethanol, or (iii) butanol. With several options to balance cofactors available, information about enzyme specificity is necessary to achieve reasonable selectivity. If constraints in a genome-scale model are simply placed around secretion of a product or byproduct, the model does not represent the cellular mechanisms that result in proper selection. Thus, an effective metabolic engineering strategy cannot be formulated in silico given these types of constraints.
With the ultimate goal of re-directing metabolic flux through the butanol production pathway in C. acetobutylicum, few tools, with the notable exception of OptKnock , exist for deriving a metabolic engineering strategy. Even with its many successes, OptKnock is restricted to gene knockouts and cannot suggest over-expression and partial gene knockdown strategies to engineer metabolism. However, the recently published OptForce algorithm  provides the capability to identify both gene over-expressions and knockdowns required of a metabolic network to produce a targeted amount of a specified product. Ultimately, methods that target the regulatory network of the cell and re-direct metabolic flux at network branch points will enable even more effective metabolic engineering in silico. The research presented here is a first step to constraining metabolic branching based on enzyme specificity. This approach also enables simulation of gene over-expressions and partial gene knockdowns in addition to gene knockouts.
Considering metabolic flux ratios
The experimental determination of metabolic flux and pathway usage through the use of isotope tracers has significantly contributed to the overall understanding of regulated metabolism. One approach to characterize metabolism is through the use of metabolic flux ratio analysis (METAFoR) [45–47]. This method is used to determine the degree of converging pathway usage to produce a metabolite pool when multiple synthesis routes exist. For example, METAFoR can reveal the relative contributions of anaplerosis and the TCA cycle to the formation of the oxaloacetate pool. Early results revealed the robustness of central carbon metabolism of Escherichia coli[45, 47] as many calculated flux ratios were found impervious to genetic perturbations. Additional computational method development led to the formulation of constraints for flux balancing from measured flux ratios [48, 49]. The resulting algorithm was effective given small metabolic networks of primary metabolism and the use of nonlinear programming methods. Unfortunately, these aspects have limited the applicability to large genome-scale metabolic networks, which often must rely on linear programming.
Genome-scale modeling with flux ratios
Of the several successful (and unsuccessful) metabolic engineering strategies applied to clostridia (many of which are not mentioned here), it was not immediately apparent which design(s) would be successful upon conception. The mutant strains had to be created in the laboratory and analyzed. From these results, hypotheses were formed that guided more advanced designs. The purpose of metabolic engineering in silico is to analyze and optimize engineering strategies a priori so that only the most promising candidates are constructed in the laboratory. While genome-scale modeling has provided the necessary platform for metabolic engineering in silico, the large number of degrees of freedom of these models has been limiting. Here, a new approach called “f lux b alance a nalysis with flux ratio s (FBrAtio)” is developed and applied. One significant advantage of FBrAtio is that flux ratio constraints are built into the stoichiometric matrix directly. This approach allows for multiple flux ratio constraints to be included simultaneously, and the flux balancing problem can be solved using simple linear programming. In particular, FBrAtio is used to show that the butanol to acetone production ratio of C. acetobutylicum increases in the presence of CoAT knockdown by antisense RNA (asRNA). This metabolic engineering strategy is also simulated in the presence of AAD knockdown and over-expression to show this method can predict these published outcomes [43, 50].
A new genome-scale model for C. acetobutylicum ATCC 824 was constructed by expanding the previously published model by Senger and Papoutsakis [11, 12]. The new model is called i CAC490 and contains 707 metabolites involved in 794 biochemical reactions, including 66 membrane transport reactions. The model includes 490 genes from the C. acetobutylicum genome. The newly updated i CAC490 model differs from the original Senger and Papoutsakis model [11, 12] in that it contains 242 more reactions (a 44% increase) and 285 more metabolites (a 68% increase). The new reactions added to create the i CAC794 model were obtained from the KEGG database  and recent literature. The i CAC490 model also contains an updated TCA cycle that operates in both oxidative and reductive directions to succinate, as shown by recent fluxomics studies [52, 53]. The model allows the export of succinate since its metabolic fate has not yet been resolved conclusively. The i CAC490 model is also fully compartmentalized and allows the presence of chemical reactions in the extracellular environment. Thermodynamic reaction reversibility constraints based on Gibb’s free energy calculations from the group contribution method [54, 55] have also been applied. The biomass equation was also updated for the i CAC490 model, using the initial version by Senger and Papoutsakis [11, 12] as a template. A nonlinear optimization procedure was applied (manuscript in preparation) to optimize the biomass equation given specific environmental conditions. The biomass equation derived for exponential growth was used extensively in simulation studies reported here. It was found that the exponential growth biomass equation could result in qualitatively accurate model predictions. It is acknowledged that an updated and dynamic biomass equation will be required to obtain model predictions that are quantitatively accurate. The reconstructed metabolic network of the i CAC490 model is included as Additional file 1. The SBML formatted model is included as Additional file 2.
Flux balance analysis
The i CAC490 genome-scale model was simulated using flux balance analysis through the COBRA toolbox . The open-source GLPK linear programming software was used to solve the flux balance equation (S · v = 0), where S is a stoichiometric coefficient matrix and v is a vector of flux values. Methods related to construction of the stoichiometric matrix and the required steady-state approximation for intracellular metabolite concentrations have been detailed elsewhere . The objective functions used for all simulations were (i) maximizing the specific growth rate of the cell while (ii) minimizing the total flux of the system.
The specific proton flux
The concept of the specific proton flux (SPF) was first introduced by Senger and Papoutsakis  and describes the total rate of proton influx/efflux through all membrane transport mechanisms. This value is negative when protons are leaving the cell and positive when protons are taken-up by the cell. For the case of C. acetobutylicum, the SPF is highly negative during exponential growth (acidogenesis) and turns slightly positive during the stationary phase (solventogenesis). In this research, the SPF was constrained to specific values by constraining the proton exchange reaction (the total flux of protons in/out of the systems boundary). The SPF range was between −30 (proton efflux) and 5 (proton influx), and the limits were chosen from experimental observations.
Next, a new row is added to the stoichiometric matrix. In this new row, two values are added (all other values in the row are zero). In the column representing the reaction catalyzed by THL (Equation 1), the coefficient 1 is added to the matrix. In the column representing the reaction catalyzed by the PTA, the coefficient −2 is added (in the new row). With these additions, when the flux balance equation (S · v = 0) is solved, the ratio of fluxes for the reactions of Equation 1 and Equation 2 will be exactly 2. If the flux ratio chosen leads to an impossible solution of the metabolic network, no solution will be found to the flux balance equation.
The goal of this research was to develop a method of constraining a metabolic network so that metabolic engineering can be performed in silico. Until now, gene knockouts have been the dominant strategy for designing metabolic engineering strategies in silico. However, flux ratios offer the ability to include over-expression and flux re-direction at key branch points in a metabolic network. These results can offer a snapshot of the metabolic potential of the engineered cell and offer the metabolic engineer an experimental target to achieve these results. Simulations performed in this research focus on using FBrAtio to reproduce metabolic engineering strategies that have been experimentally validated in C. acetobutylicum[6, 43, 50]. In particular, simulations were performed with the i CAC490 model in which the glucose uptake rate and the SPF were the only specified membrane transport fluxes. Next, FBrAtio was applied to achieve the experimentally observed wild-type metabolic activity of C. acetobutylicum. The following metabolic characteristics were sought on a qualitative level: (i) at highly negative values of SPF (acidogenesis), acetate and butyrate are produced in high quantities, (ii) high hydrogen production accompanies acidogenesis, (iii) solvents are produced at SPF values close to zero and slightly positive, (iv) hydrogen production decreases during solventogenesis, (v) the maximum growth rate of the culture occurs during acidogenesis, (vi) the production of butyrate is slightly greater than the production of acetate and much greater than the production of lactate, and (vii) the production of butanol is greater than the production of acetone and is much greater than the production of ethanol. Following obtaining a qualitatively accurate simulation of wild-type metabolism, additional flux ratios were applied through FBrAtio in attempt to predict the following experimental observations [43, 50]. Knockdown of the CoAT (by asRNA) resulted in increased butanol to acetone selectivity, but this strategy resulted in decreased ethanol and butanol production . The asRNA was designed against the mRNA of the cftB gene in particular, which is a part of the tricistronic operon (aad-ctfA-ctfB). It was hypothesized that AAD activity was also compromised by this asRNA construct, so aad was over-expressed under its native promoter. Significantly higher ethanol and butanol yields were observed as a result of this metabolic engineering strategy . Flux ratios were designed to (i) knockdown CoAT activity only, (ii) knockdown activity of both CoAT and AAD, and (iii) knockdown CoAT while over-expressing AAD at and above wild-type levels.
Simulations with a minimal set of constraints
The i CAC490 model was simulated with a glucose uptake rate constrained to 10 and the SPF was varied between −30 and 5 . Only thermodynamic reversibility constraints were used initially. Results showed acidogenic and solventogenic metabolic phases that coincided with SPF values . Results also showed a maximum specific growth rate at an SPF value of −10 , which is consistent with previous findings [2, 12]. However, during acidogenesis, acetate was the primary acid produced, and acetone was the primary solvent produced during solventogenesis. Hydrogen (H2) production was also maximized during solventogenesis. These characteristics are not consistent with experimental observations. Since only acetate was produced in acidogenesis, this demonstrates that the network required the generation of ATP. Since butyrate or ethanol was not produced, this means that NAD+ was regenerated in a futile cycle elsewhere in the network. The ability of the network to artificially balance NAD(P)+/NAD(P)H also explains why hydrogen production remained high during solventogenesis. By maximizing the specific growth rate of the cell and minimizing the total flux of the system, flux in longer pathways, such as butyrate/butanol production were minimized in favor of shorter ATP (acetate) and NAD+ (futile cycle) regenerating pathways.
Approximating wild-type metabolism with FBrAtio
Reactions and updated constraints involving the ferredoxins
PFO: pyruvate ferredoxin oxidoreductase
FNO: ferredoxin NAD+ oxidoreducatase
FNPO: ferredoxin NADP+ oxidoreductase
asRNA knockdown of CoAT only
Knockdown of both CoAT and AAD
Over-express AAD and knockdown CoAT
The large number of degrees of freedom in primary clostridial metabolism makes this system challenging to model. Initial efforts [1, 7, 8] relied on experimentally measured data to fit a basic metabolic model and back-calculate pathways fluxes. With the development of genome-scale models, there was initial enthusiasm that this approach would result in a model capable of predicting the metabolic response of the organism to genetic and environmental manipulations. However, this level of prediction was not achieved by the first genome-scale model for C. acetobutylicum[11, 12]. This original genome-scale model was updated in this research with additional reactions and thermodynamic constraints. Even with a more complete model and updated constraints, the number of degrees of freedom of the primary metabolic network proved too large to generate meaningful predictions, even of wild-type metabolism. This is evident from the results shown in Figure 2. To build a truly predictive model, care must be taken when determining how proper constraints are imposed. It is important that these constraints not only lead to accurate representations of metabolism but can be manipulated to mimic genetic and environmental perturbations. For example, a common method is to artificially constrain the glucose uptake rate (as was done in this research). From there, constraints can be imposed on product (e.g., acetate, butyrate, butanol, etc.) secretion fluxes to mimic the wild-type metabolism. This approach is detrimental to metabolic engineering. For example, if constraints are placed on secretion of the end-products, how do these constraints change when a genetic manipulation is made elsewhere in the metabolic network (e.g., at the thiolase enzyme)? There is no clear mathematical relationship between a secretion flux constraint and the metabolic flux through an enzyme elsewhere in the network. Thus, constraints that are imposed to achieve accurate representations of metabolism must be imposed at the metabolic engineering targets themselves. However, this leads to the questions, what is a metabolic engineering target? And, how can constraints be imposed there? This research has focused on “branch points” (or critical nodes) of the metabolic network as potential sites of metabolic engineering. The use of acetyl-CoA in clostridial metabolism is a good example of a metabolic branch point. Acetyl-CoA can be used in the production of (i) acetate, (ii) butyrate/butanol, (iii) ethanol, and (iv) macromolecules required for cell growth. Each of these routes produces/consumes different cofactors, and the balancing of these cofactors ultimately determines the cellular phenotype.
The use of metabolic flux ratio constraints through FBrAtio enabled qualitatively accurate modeling of acidogenic and solventogenic metabolism of C. acetobutylicum using the new i CAC490 genome-scale model. The use of flux ratios allows for constraints to be placed directly at points where metabolic engineering strategies can be applied. For example, flux ratios can be manipulated to achieve a desired result (e.g., maximized butanol production). Then, genetic manipulations such as (i) over-expression, (ii) knockout, and (iii) asRNA knockdown can be applied to achieve the optimum ratios. In this research, flux ratio constraints were implemented to achieve a qualitative picture of metabolism that mimics experimental observations. As a proof of concept, the wild-type and two engineered strains analyzed were consistent with published experimental results. The case of AAD over-expression went a step further and exposed a possible metabolic engineering limit to re-routing flux into the alcohol production pathways. This suggests that the approach of flux ratio constraints is tunable. The flux values obtained here were not converted into concentrations of metabolites and biomass and compared directly to published values. The results obtained here are qualitative (not quantitative) pictures of metabolism. There are several reasons for this. First, a fixed glucose uptake rate of 10 was used for all values of the SPF examined. Previous results  have shown that the glucose uptake rate varies with the SPF. However, the relationship between the glucose uptake rate and the SPF remains uncharacterized. At best, a causal relationship can be established between these two with the current level of knowledge. Next, a single biomass equation was used for all values of the SPF examined. Previous research has shown that the biomass composition, including the maintenance ATP requirement, of C. acetobutylicum changes with the SPF [10, 12]. To obtain quantitatively accurate predictions, one must first understand the relationships that exist between glucose uptake and biomass composition with the SPF. While research is underway to uncover these relationships, the use of parameters associated with exponential growth seemed to be sufficient with the FBrAtio approach.
FBrAtio is a new method to derive metabolic engineering strategies to achieve optimum phenotypes. The concept of using metabolic flux ratios was initially developed with the METAFoR approach [45, 47]. It enabled researchers to determine how multiple biosynthetic pathways contributed to the production of a metabolite pool. This enabled identification of new metabolic pathways and regulatory mechanisms. Since the implementation of FBrAtio accommodates the use of linear programming, flux ratios found with METAFoR can now easily be applied to appropriate genome-scale models using the techniques described in the Methods section (see Equations 14). The FBrAtio approach is different from METAFoR in that it considers how a metabolite pool is distributed as a substrate among competing enzymes. Of course, this process is governed by thermodynamics. This means that enzyme availability and intermediate accumulation downstream (among other factors) are responsible for flux ratios in physical systems. The FBrAtio approach can lead a metabolic engineer to optimum flux ratios, and enzyme availability can be manipulated through gene (i) over-expression, (ii) knockout, or (iii) partial knockdown. However, the FBrAtio approach cannot predict the potential accumulation of downstream intermediates once flux is redirected. This remains a problem for the experimentalist that may be addressed through additional gene over-expression or enzyme engineering.
The FBrAtio approach is presented in detail here and is applied to model previously published metabolic engineering approaches in C. acetobutylicum. Obviously, the full potential of FBrAtio will be realized when it can be used systematically. To do this, algorithms are needed to identify critical nodes (metabolite pools) in the metabolic network where flux ratios can be optimized to produce a desired phenotype. Research is currently underway to address this challenging task. The end result will provide the metabolic engineer with a list of flux ratios that can be manipulated using existing toolsets. Although additional complications may be encountered in some cases due to unforeseen regulatory interactions, the FBrAtio approach has the potential to provide effective “fine-tuned” metabolic engineering strategies.
The FBrAtio approach for -incorporating metabolic flux ratio constraints into a genome-scale metabolic network and generating solutions using simple linear programming was developed in this research. The approach proved effective in modeling wild-type metabolism of C. acetobutylicum. FBrAtio was then applied to metabolically engineered strains, and a high ethanol producing strain was effectively modeled. A nonlinear relationship exists between the flux ratios at a critical node and the resulting phenotype. FBrAtio is capable of capturing these nonlinearities. How flux ratio constraints can be used to design metabolic engineering strategies is currently a subject of much future research, and the developments presented here represent the first steps toward truly predictive genome-scale models that can accurately reflect the impacts of genetic and environmental manipulations.
This research was funded by the USDA AFRI Biobased Products and Bioenergy Production Program (Award Number 2010-65504-20346). BF was supported by the Doctoral Scholars Program of the Institute of Critical Technologies and Applied Science (ICTAS) at Virginia Tech. JY was supported by a grant from the Biodesign and Bioprocessing Research Center at Virginia Tech.
- Papoutsakis ET: Equations and calculations for fermentations of butyric acid bacteria. Biotechnol Bioeng 1984, 26: 174-187. 10.1002/bit.260260210View ArticleGoogle Scholar
- Monot F, Martin JR, Petitdemange H, Gay R: Acetone and butanol production by Clostridium acetobutylicum in a synthetic medium. Appl Environ Microbiol 1982, 44: 1318-1324.Google Scholar
- Girbal L, Soucaille P: Regulation of Clostridium-acetobutylicum metabolism as revealed by mixed-substrate steady-state continuous cultures - Role of Nadh/Nad ratio and Atp Pool. J Bacteriol 1994, 176: 6433-6438.Google Scholar
- Papoutsakis ET, Desai RP, Nielsen LK: Stoichiometric modeling of Clostridium acetobutylicum fermentations with non-linear constraints. J Biotechnol 1999, 71: 191-205. 10.1016/S0168-1656(99)00022-XView ArticleGoogle Scholar
- Calik P, Ozdamar TH: Bioreaction network flux analysis for industrial microorganisms: A review. Rev Chem Eng 2002, 18: 553-596.View ArticleGoogle Scholar
- Sillers R, Chow A, Tracy B, Papoutsakis ET: Metabolic engineering of the non-sporulating, non-solventogenic Clostridium acetobutylicum strain M5 to produce butanol without acetone demonstrate the robustness of the acid-formation pathways and the importance of the electron balance. Metab Eng 2008, 10: 321-332. 10.1016/j.ymben.2008.07.005View ArticleGoogle Scholar
- Desai RP, Harris LM, Welker NE, Papoutsakis ET: Metabolic flux analysis elucidates the importance of the acid-formation pathways in regulating solvent production by Clostridium acetobutylicum . Metab Eng 1999, 1: 206-213. 10.1006/mben.1999.0118View ArticleGoogle Scholar
- Desai RP, Nielsen LK, Papoutsakis ET: Stoichiometric modeling of Clostridium acetobutylicum fermentations with non-linear constraints. J Biotechnol 1999, 71: 191-205. 10.1016/S0168-1656(99)00022-XView ArticleGoogle Scholar
- Lee J, Yun H, Feist AM, Palsson BO, Lee SY: Genome-scale reconstruction and in silico analysis of the Clostridium acetobutylicum ATCC 824 metabolic network. Appl Microbiol Biotechnol 2008, 80: 849-862. 10.1007/s00253-008-1654-4View ArticleGoogle Scholar
- Senger RS: Biofuel production improvement with genome-scale models: The role of cell composition. Biotechnol J 2010, 5: 671-685. 10.1002/biot.201000007View ArticleGoogle Scholar
- Senger RS, Papoutsakis ET: Genome-scale model for Clostridium acetobutylicum : Part I. Metabolic network resolution and analysis. Biotechnol Bioeng 2008, 101: 1036-1052. 10.1002/bit.22010View ArticleGoogle Scholar
- Senger RS, Papoutsakis ET: Genome-scale model for Clostridium acetobutylicum : Part II. Development of specific proton flux states and numerically determined sub-systems. Biotechnol Bioeng 2008, 101: 1053-1071. 10.1002/bit.22009View ArticleGoogle Scholar
- Milne CB, Kim PJ, Eddy JA, Price ND: Accomplishments in genome-scale in silico modeling for industrial and medical biotechnology. Biotechnol J 2009, 4: 1653-1670. 10.1002/biot.200900234View ArticleGoogle Scholar
- Mahadevan R, Palsson BO, Lovley DR: In situ to in silico and back: elucidating the physiology and ecology of Geobacter spp. using genome-scale modelling. Nat Rev Microbiol 2011, 9: 39-50. 10.1038/nrmicro2456View ArticleGoogle Scholar
- Feist AM, Palsson BO: The growing scope of applications of genome-scale metabolic reconstructions using Escherichia coli . Nat Biotechnol 2008, 26: 659-667. 10.1038/nbt1401View ArticleGoogle Scholar
- Feist AM, Herrgard MJ, Thiele I, Reed JL, Palsson BO: Reconstruction of biochemical networks in microorganisms. Nat Rev Microbiol 2009, 7: 129-143.View ArticleGoogle Scholar
- Papin JA, Price ND, Wiback SJ, Fell DA, Palsson BO: Metabolic pathways in the post-genome era. Trends Biochem Sci 2003, 28: 250-258. 10.1016/S0968-0004(03)00064-1View ArticleGoogle Scholar
- Papin JA, Stelling J, Price ND, Klamt S, Schuster S, Palsson BO: Comparison of network-based pathway analysis methods. Trends Biotechnol 2004, 22: 400-405. 10.1016/j.tibtech.2004.06.010View ArticleGoogle Scholar
- Price ND, Papin JA, Schilling CH, Palsson BO: Genome-scale microbial in silico models: the constraints-based approach. Trends Biotechnol 2003, 21: 162-169. 10.1016/S0167-7799(03)00030-1View ArticleGoogle Scholar
- Price ND, Reed JL, Palsson BO: Genome-scale models of microbial cells: evaluating the consequences of constraints. Nat Rev Microbiol 2004, 2: 886-897. 10.1038/nrmicro1023View ArticleGoogle Scholar
- Reed JL, Famili I, Thiele I, Palsson BO: Towards multidimensional genome annotation. Nat Rev Genet 2006, 7: 130-141. 10.1038/nrg1769View ArticleGoogle Scholar
- Teusink B, Wiersma A, Molenaar D, Francke C, de Vos WM, Siezen RJ, Smid EJ: Analysis of growth of Lactobacillus plantarum WCFS1 on a complex medium using a genome-scale metabolic model. J Biol Chem 2006, 281: 40041-40048. 10.1074/jbc.M606263200View ArticleGoogle Scholar
- Gonzalez O, Gronau S, Falb M, Pfeiffer F, Mendoza E, Zimmer R, Oesterhelt D: Reconstruction, modeling & analysis of Halobacterium salinarum R-1 metabolism. Mol Biosyst 2008, 4: 148-159. 10.1039/b715203eView ArticleGoogle Scholar
- Oliveira AP, Nielsen J, Forster J: Modeling Lactococcus lactis using a genome-scale flux model. BMC Microbiol 2005, 5: 39. 10.1186/1471-2180-5-39View ArticleGoogle Scholar
- Gianchandani EP, Papin JA, Price ND, Joyce AR, Palsson BO: Matrix formalism to describe functional states of transcriptional regulatory systems. PLoS Comput Biol 2006, 2: e101. 10.1371/journal.pcbi.0020101View ArticleGoogle Scholar
- Chandrasekaran S, Price ND: Probabilistic integrative modeling of genome-scale metabolic and regulatory networks in Escherichia coli and Mycobacterium tuberculosis . Proc Natl Acad Sci U S A 2010, 107: 17845-17850. 10.1073/pnas.1005139107View ArticleGoogle Scholar
- Reed JL, Vo TD, Schilling CH, Palsson BO: An expanded genome-scale model of Escherichia coli K-12 (iJR904 GSM/GPR). Genome Biol 2003, 4: R54. 10.1186/gb-2003-4-9-r54View ArticleGoogle Scholar
- Archer CT, Kim JF, Jeong H, Park JH, Vickers CE, Lee SY, Nielsen LK: The genome sequence of E. coli W (ATCC 9637): comparative genome analysis and an improved genome-scale reconstruction of E. coli. BMC Genomics 2011, 12: 9. 10.1186/1471-2164-12-9View ArticleGoogle Scholar
- Perumal D, Samal A, Sakharkar KR, Sakharkar MK: Targeting multiple targets in Pseudomonas aeruginosa PAO1 using flux balance analysis of a reconstructed genome-scale metabolic network. J Drug Target 2011, 19: 1-13. 10.3109/10611861003649753View ArticleGoogle Scholar
- Kim HU, Kim TY, Lee SY: Genome-scale metabolic network analysis and drug targeting of multi-drug resistant pathogen Acinetobacter baumannii AYE. Mol Biosyst 2010, 6: 339-348. 10.1039/b916446dView ArticleGoogle Scholar
- AbuOun M, Suthers PF, Jones GI, Carter BR, Saunders MP, Maranas CD, Woodward MJ, Anjum MF: Genome scale reconstruction of a Salmonella metabolic model: comparison of similarity and differences with a commensal Escherichia coli strain. J Biol Chem 2009, 284: 29480-29488. 10.1074/jbc.M109.005868View ArticleGoogle Scholar
- Rambaut A, Pybus OG, Nelson MI, Viboud C, Taubenberger JK, Holmes EC: The genomic and epidemiological dynamics of human influenza A virus. Nature 2008, 453: 615-619. 10.1038/nature06945View ArticleGoogle Scholar
- Heinemann M, Kummel A, Ruinatscha R, Panke S: In silico genome-scale reconstruction and validation of the Staphylococcus aureus metabolic network. Biotechnol Bioeng 2005, 92: 850-864. 10.1002/bit.20663View ArticleGoogle Scholar
- Fong SS, Burgard AP, Herring CD, Knight EM, Blattner FR, Maranas CD, Palsson BO: In silico design and adaptive evolution of Escherichia coli for production of lactic acid. Biotechnol Bioeng 2005, 91: 643-648. 10.1002/bit.20542View ArticleGoogle Scholar
- Kuepfer L, Sauer U, Blank LM: Metabolic functions of duplicate genes in Saccharomyces cerevisiae . Genome Res 2005, 15: 1421-1430. 10.1101/gr.3992505View ArticleGoogle Scholar
- Burgard AP, Pharkya P, Maranas CD: Optknock: a bilevel programming framework for identifying gene knockout strategies for microbial strain optimization. Biotechnol Bioeng 2003, 84: 647-657. 10.1002/bit.10803View ArticleGoogle Scholar
- Pharkya P, Burgard AP, Maranas CD: OptStrain: a computational framework for redesign of microbial production systems. Genome Res 2004, 14: 2367-2376. 10.1101/gr.2872004View ArticleGoogle Scholar
- Ranganathan S, Suthers PF, Maranas CD: OptForce: an optimization procedure for identifying all genetic manipulations leading to targeted overproductions. PLoS Comput Biol 2010, 6: e1000744. 10.1371/journal.pcbi.1000744View ArticleGoogle Scholar
- Jones SW, Paredes CJ, Tracy B, Cheng N, Sillers R, Senger RS, Papoutsakis ET: The transcriptional program underlying the physiology of clostridial sporulation. Genome Biol 2008, 9: R114. 10.1186/gb-2008-9-7-r114View ArticleGoogle Scholar
- Mermelstein LD, Papoutsakis ET: In vivo methylation in Escherichia coli by the Bacillus subtilis phage phi 3 T I methyltransferase to protect plasmids from restriction upon transformation of Clostridium acetobutylicum ATCC 824. Appl Environ Microbiol 1993, 59: 1077-1081.Google Scholar
- Green EM, Boynton ZL, Harris LM, Rudolph FB, Papoutsakis ET, Bennett GN: Genetic manipulation of acid formation pathways by gene inactivation in Clostridium acetobutylicum ATCC 824. Microbiology 1996, 142: 2079-2086. 10.1099/13500872-142-8-2079View ArticleGoogle Scholar
- Harris LM, Desai RP, Welker NE, Papoutsakis ET: Characterization of recombinant strains of the Clostridium acetobutylicum butyrate kinase inactivation mutant: need for new phenomenological models for solventogenesis and butanol inhibition? Biotechnol Bioeng 2000, 67: 1-11. 10.1002/(SICI)1097-0290(20000105)67:1<1::AID-BIT1>3.0.CO;2-GView ArticleGoogle Scholar
- Tummala SB, Junne SG, Papoutsakis ET: antisense RNA downregulation of coenzyme a transferase combined with alcohol-aldehyde dehydrogenase overexpression leads to predominantly alcohologenic Clostridium acetobutylicum fermentations. J Bacteriol 2003, 185: 3644-3653. 10.1128/JB.185.12.3644-3653.2003View ArticleGoogle Scholar
- Sillers R, Al-Hinai MA, Papoutsakis ET: Aldehyde–alcohol dehydrogenase and/or thiolase overexpression coupled with CoA transferase downregulation lead to higher alcohol titers and selectivity in Clostridium acetobutylicum fermentations. Biotechnol Bioeng 2009, 102: 38-49. 10.1002/bit.22058View ArticleGoogle Scholar
- Fischer E, Sauer U: Metabolic flux profiling of Escherichia coli mutants in central carbon metabolism using GC-MS. European journal of biochemistry /FEBS 2003, 270: 880-891. 10.1046/j.1432-1033.2003.03448.xView ArticleGoogle Scholar
- Fischer E, Zamboni N, Sauer U: High-throughput metabolic flux analysis based on gas chromatography–mass spectrometry derived 13 C constraints. Anal Biochem 2004, 325: 308-316. 10.1016/j.ab.2003.10.036View ArticleGoogle Scholar
- Sauer U, Lasko DR, Fiaux J, Hochuli M, Glaser R, Szyperski T, Wuthrich K, Bailey JE: Metabolic flux ratio analysis of genetic and environmental modulations of Escherichia coli central carbon metabolism. J Bacteriol 1999, 181: 6679-6688.Google Scholar
- Buhler B, Park JB, Blank LM, Schmid A: NADH availability limits asymmetric biocatalytic epoxidation in a growing recombinant Escherichia coli strain. Appl Environ Microbiol 2008, 74: 1436-1446. 10.1128/AEM.02234-07View ArticleGoogle Scholar
- Sauer U, Hatzimanikatis V, Bailey JE, Hochuli M, Szyperski T, Wuthrich K: Metabolic fluxes in riboflavin-producing Bacillus subtilis . Nat Biotechnol 1997, 15: 448-452. 10.1038/nbt0597-448View ArticleGoogle Scholar
- Tummala SB, Welker NE, Papoutsakis ET: Design of antisense RNA constructs for downregulation of the acetone formation pathway of Clostridium acetobutylicum . J Bacteriol 2003, 185: 1923-1934. 10.1128/JB.185.6.1923-1934.2003View ArticleGoogle Scholar
- Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res 2000, 28: 27-30. 10.1093/nar/28.1.27View ArticleGoogle Scholar
- Crown SB, Indurthi DC, Ahn WS, Choi J, Papoutsakis ET, Antoniewicz MR: Resolving the TCA cycle and pentose-phosphate pathway of Clostridium acetobutylicum ATCC 824: Isotopomer analysis, in vitro activities and expression analysis. Biotechnol J 2011, 6: 300-305. 10.1002/biot.201000282View ArticleGoogle Scholar
- Amador-Noguez D, Feng XJ, Fan J, Roquet N, Rabitz H, Rabinowitz JD: Systems-level metabolic flux profiling elucidates a complete, bifurcated tricarboxylic acid cycle in Clostridium acetobutylicum . J Bacteriol 2010, 192: 4452-4461. 10.1128/JB.00490-10View ArticleGoogle Scholar
- Jankowski MD, Henry CS, Broadbelt LJ, Hatzimanikatis V: Group contribution method for thermodynamic analysis of complex metabolic networks. Biophys J 2008, 95: 1487-1499. 10.1529/biophysj.107.124784View ArticleGoogle Scholar
- Henry CS, Broadbelt LJ, Hatzimanikatis V: Thermodynamics-based metabolic flux analysis. Biophys J 2007, 92: 1792-1805. 10.1529/biophysj.106.093138View ArticleGoogle Scholar
- Becker SA, Feist AM, Mo ML, Hannum G, Palsson BO, Herrgard MJ: Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox. Nat Protoc 2007, 2: 727-738. 10.1038/nprot.2007.99View ArticleGoogle Scholar
- Edwards JS, Ibarra RU, Palsson BO: In silico predictions of Escherichia coli metabolic capabilities are consistent with experimental data. Nat Biotechnol 2001, 19: 125-130. 10.1038/84379View ArticleGoogle Scholar
- Harris LM, Blank L, Desai RP, Welker NE, Papoutsakis ET: Fermentation characterization and flux analysis of recombinant strains of Clostridium acetobutylicum with an inactivated solR gene. J Ind Microbiol Biotechnol 2001, 27: 322-328. 10.1038/sj.jim.7000191View ArticleGoogle Scholar
- Meyer CL, Papoutsakis ET: Increased levels of ATP and NADH are associated with increased solvent production in continuous cultures of Clostridium acetobutylicum . Appl Environ Microbiol 1989, 30: 450-459.Google Scholar