A systems biology approach to investigate the effect of pH-induced gene regulation on solvent production by Clostridium acetobutylicum in continuous culture
- Sylvia Haus†1,
- Sara Jabbari†2Email author,
- Thomas Millat1,
- Holger Janssen3,
- Ralf-Jörg Fischer3,
- Hubert Bahl3,
- John R King2 and
- Olaf Wolkenhauer1, 4
© Haus et al; licensee BioMed Central Ltd. 2011
Received: 20 July 2010
Accepted: 19 January 2011
Published: 19 January 2011
Clostridium acetobutylicum is an anaerobic bacterium which is known for its solvent-producing capabilities, namely regarding the bulk chemicals acetone and butanol, the latter being a highly efficient biofuel. For butanol production by C. acetobutylicum to be optimized and exploited on an industrial scale, the effect of pH-induced gene regulation on solvent production by C. acetobutylicum in continuous culture must be understood as fully as possible.
We present an ordinary differential equation model combining the metabolic network governing solvent production with regulation at the genetic level of the enzymes required for this process. Parameterizing the model with experimental data from continuous culture, we demonstrate the influence of pH upon fermentation products: at high pH (pH 5.7) acids are the dominant product while at low pH (pH 4.5) this switches to solvents. Through steady-state analyses of the model we focus our investigations on how alteration in gene expression of C. acetobutylicum could be exploited to increase butanol yield in a continuous culture fermentation.
Incorporating gene regulation into the model of solvent production by C. acetobutylicum enables an accurate representation of the pH-induced switch to solvent production to be obtained and theoretical investigations of possible synthetic-biology approaches to be pursued. Steady-state analyses suggest that, to increase butanol yield, alterations in the expression of single solvent-associated genes are insufficient; a more complex approach targeting two or more genes is required.
A renewed interest in the development of biofuels is emerging as a result of a variety of factors including dwindling crude oil reserves, concerns over the environmental impact of fossil fuels and threats to national security potentially limiting access to resources . In recent years, biofuels have been predominantly sourced from crops, resulting in competition for limited food resources and land [2–5]; bacterial fermentation has been considered a possible answer to this problem . One of the best-studied bacteria for biofuel production is Clostridium acetobutylicum. Clostridial bacteria are strictly anaerobic, Gram-positive and form highly-resistant spores. Many of the clostridial species, such as Clostridium difficile and Clostridium botulinum, are highly pathogenic and cause devastating diseases. Some, however, like C. acetobutylicum which was first isolated from corn in 1912 by Chaim Weizmann , are harmless to humans, animals and plants and make a wide range of useful chemicals .
The metabolism of C. acetobutylicum is characterized by the so called acetone-butanol fermentation (AB fermentation) which is also referred to as acetone-butanol-ethanol (ABE) fermentation. Since butanol is a more efficient biofuel than many other solvents such as ethanol, much research is currently focused on this bacterium [1, 9]. The metabolic pathway of AB fermentation comprises two characteristic phases: acidogenesis and solventogenesis. During the transition phase the generation of the solvents acetone and butanol is induced and these become the dominant fermentation products. This switch is called the metabolic (or solventogenic) shift. Though the metabolic pathways leading to solvent and acid production are clearly defined , the mechanisms governing the shift are not well understood. It has been shown, however, that under continuous culture conditions in a chemostat, the external pH is a crucial prerequisite for the induction of this metabolic change [11, 12].
Experiments showed that in continuous culture C. acetobutylicum shifts its metabolism in response to the external pH: cells produce predominantly acetate and butyrate above a pH of 5.0 and acetone and butanol when the pH is lower than 5.0 [11, 12, 16, 17]. It is postulated that the switch to solventogenesis represents one aspect of a survival strategy in response to external pH [1, 10]. Concerning this aspect it is noteworthy that in contrast to many other organisms, clostridial bacteria seem to be unable to maintain the internal pH value at a more or less constant level. Instead they generate an approximately constant pH gradient across the cell membrane which results in an internal pH value of about one unit higher than the external one . Consequently, the internal pH follows the external one without significant delay .
While many research groups work with C. acetobutylicum in batch culture, we focus on continuous culture experiments in order to investigate the effect of pH on the AB metabolism . The 'forward' dynamic-shift experiment runs in the following way: after approximately 65 hours (requiring five volume changes, each after 13.3 hours) the culture is in steady state in the acid producing phase. This status can be kept for up to ten days. After this period the mega plasmid pSol1 which is essential for solvent production will be lost and, with it, the ability to produce solvents . To initiate the metabolic shift, control of the pH of the culture via KOH is turned off. After approximately 65 hours the bacterial population is in solventogenic steady state and can, in principle, be held permanently in this state.
In this study we introduce a model at the metabolic, proteomic, and genomic levels to explain the pH-induced changes in AB metabolism and to elucidate the key genes governing this shift. We perform steady-state analyses of solventogenesis to study the effect on metabolic outputs of artificially altering gene regulation with the aim of predicting targets for optimal genetically engineered production strains with respect to butanol. The model is parameterized using data from our dynamic-shift experiments. Our findings indicate that mutations to, or otherwise altered expression of, single genes in the solvent-producing pathway are insufficient to increase significantly the butanol fermentation; a more complicated approach involving alterations in the expression of multiple genes is required.
A further important aspect of the AB fermentation is the flow of energetic and electron carriers and their availability which was investigated in the 1980 s and 1990 s by several experimental groups [19, 20]. The results suggest that ATP and NADH levels might be important for the regulation of AB fermentation. Nevertheless, for this first dynamic model presented in this manuscript we do not consider this aspect in order to simplify the investigations of the pH-induced metabolic switch and its optimization with respect to butanol production.
Experimental setup in continuous culture
In contrast to many other groups handling batch cultures of C. acetobutylicum, we used continuous cultures in order to study the effect of pH on solvent production, in particular on the production of butanol. A fermentation process that can be operated continuously has several advantages over a batch process in industrial and biotechnological manufacturing: only one series of precultures is needed for a long production period, the "dead season" necessary for the filling, sterilization, cooling and clearing of the equipment is largely diminished and the volume of the fermenter vessel can be reduced without a loss of production capacity . Stable solvent production can be maintained for much longer in a synthetic medium under phosphate limitation in a chemostat culture than in the traditional batch process .
Our model is based around continuous culture dynamic shift experiments. Experiments were performed according to a standardized experimental setup [16, 18] and using standard operating procedures for extracting and handling different types of samples. The strain C. acetobutylicum ATCC824 was grown under anaerobic conditions at 37°C and the precultures were prepared as previously described . The phosphate-limited chemostat experiments were performed in a BiostatB 1.5-l fermenter system (BBI, Melsungen, Germany) with 0.5 mM KH2PO4 and 4% (wt/vol) glucose in the supplying medium  and a dilution rate of D = 0.075 h-1. The external pH in the culture medium was adjusted to and kept constant at pH 5.7 and pH 4.5 by automatic addition of 2 M KOH. The analysis of the fermentation products was accomplished as described previously . Three individual experiments were performed shifting the culture from pH 5.7 to pH 4.5 (which we call the 'forward shift' experiments) and one from pH 4.5 to pH 5.7 (the 'reverse shift' experiment). Data were taken over the full length of the observation time. The biological data is given in Additional file 1: Tables S1 - S4.
pH-dependent modeling of the AB fermentation
As mentioned before, the metabolism of AB fermentation in C. acetobutylicum displays two characteristic metabolic phases, acidogenesis and solventogenesis. A detailed representation of the corresponding metabolic pathways has been published by Jones and Woods . Up to the present, only a few metabolic models describing this fermentative process have been published. Papoutsakis  developed a stoichiometric model that could be used to estimate the rates of reactions occurring within the AB fermentation pathways of several (AB-producing) clostridial bacteria in batch culture. Importantly, however, these results are not transferable to continuous cultivations like the chemostat cultures used in this study where the cells are growing exponentially throughout. Votruba et al.  formulated a model of the fermentation process for batch cultures without including the intermediate metabolites, restricting the variables to biomass, glucose and the end products and this model captures well the two phases. Metabolic flux analysis has been applied since to the pathway [25–27]. Existing kinetic models describing the dynamics of ABE fermentation do not seek to capture the effect of pH upon the metabolic network: in Shinto et al.  (which actually considers Clostridium saccharoperbutylacetonicum) the switch is assumed to be glucose-dependent and enzyme levels are taken to be constant, while in Li et al.  enzyme activity is incorporated, but regulation is fed in directly to the model from experimental data (rather than being allowed to vary freely within the model or be influenced by other model components). To our knowledge, therefore, the model presented here is the first to consider the effect of pH upon the metabolic network.
The reduced metabolic network
Reactions that form the metabolic pathways of acidogenesis and solventogenesis and the enzymes which are required to catalyze the solvent-associated reactions, i.e. 3-7 and 9.
glucose → acetyl-CoA
acetyl-CoA → acetate
acetate + acetoacetyl-CoA → acetoacetate + acetyl-CoA
acetyl-CoA → acetoacetyl-CoA
acetyl-CoA → ethanol
butyrate + acetoacetyl-CoA → acetoacetate + butyryl-CoA
acetoacetate → acetone
butyryl-CoA → butyrate
butyryl-CoA → butanol
acetoacetyl-CoA → butyryl-CoA
When gene regulation is not included explicitly in the representation of a reaction, we employ Michaelis-Menten expressions. In the following section we explain our approach when gene regulation is included explicitly.
Incorporating gene regulation
The reactions required for solventogenesis are tightly regulated at the genetic level: production of the enzymes required to catalyze these reactions can be switched on or off (or, more generally, increased or decreased) at the transcriptional level by regulatory proteins binding at the appropriate DNA sites. The levels of the specific regulatory proteins are adjusted in response to both internal and external signals (for example pH) that are transmitted through the cell. It is therefore expected that the switch from acidogenesis to solventogenesis (or, indeed, vice versa), can be explained, at least partially, via pH-regulation of the enzyme-associated genes.
The principal enzymes involved in solventogenesis are encoded by the genes adc, adhE, bdhA/B, ctfA/B and thlA - see Table 1. We note that there are two adhE genes in C. acetobutylicum; the one we incorporate into our model is sometimes referred to as adhE1. For reasons outlined below, we include the influence of only three of these genes in our dynamic model, namely adc, adhE and ctfA/B.
Expression of the ctfA/B gene results in the acetoacetyl-CoA: acetate/butyrate:CoA-transferase  (or simply CoA-transferase), which is the enzyme responsible for converting acetoacetyl-CoA and the previously secreted acetate or butyrate into acetoacetate and acetyl-CoA or butyryl-CoA, respectively. The adhE gene is part of the same operon as ctfA/B, namely the sol operon . The dehydrogenase AdhE catalyzes the conversion of butyryl-CoA and acetyl-CoA into butanol and ethanol, respectively [31, 32].
The adc gene encodes acetoacetate decarboxylase (Adc), which is the enzyme responsible for the decarboxylation of acetoacetate into acetone and CO2[14, 15]. There are two known bdh genes: bdhA and bdhB; the products of both genes are butanol dehydrogenases  which (in addition to AdhE) convert butyryl-CoA into butanol. Since their roles are similar to that of AdhE their effects can be absorbed into the parameter choice for AdhE. Thus, we neglected them from our dynamic model. The gene product of thlA is a thiolase which catalyzes the conversion of acetyl-CoA into acetoacetyl-CoA and is therefore needed for both acidogenesis and solventogenesis . It is speculated that expression of the thlA gene is constitutive  and recent experiments indicate that the difference in transcription levels of thlA between acidogenesis and solventogenesis is much less marked than that of ctfA/B, adc, and adhE. Therefore, ThlA was also neglected from our dynamic model as we assume ThlA levels remain relatively constant throughout the experiments as found in . The adc, adhE and ctfA/B genes are regulated in accordance with the metabolic phase, with their transcription increasing in the solventogenic phase . Thus, they are likely to be induced via a pH-related signal and we therefore incorporate each of these genes into the model. We assume that they are each transcribed at two distinct rates: basal (r E for enzyme E, see the dashed downward arrow in (1) below) when the pH is sufficiently high (indeed, low levels of these transcripts exist before the onset of solventogenesis, see for example [37, 38]) and at a higher rate ( , see dashed upward arrow in (1)) when the pH of the environment is low.
where p* represents the threshold pH level around which the switch occurs and n dictates the steepness of the smooth switch function.
to the pH data, where c1,c2 and c3 are constants, to gain a distinct function representing pH for each simulation; this function is then fed into the model via the switch function F(p).
Combining the pH-dependent metabolic network with gene regulation
where we have added to each equation an out-flow term which is the product of the dilution rate with the concentration of the corresponding metabolite because we have a constant out-flow of both extracellular and intracellular (as a result of cell out-flow) products through the chemostat.
We fit the model to experimental data as described in the following section.
Data normalization, parameter estimation and model simulation
Estimated parameter values.
Steady-state analysis of the metabolic network
Wild-type-associated production rates (to three significant figures) for each enzyme in the steady-state investigations.
Results and Discussion
Modeling and Simulation of the dynamic shift experiments
Here we apply our dynamic process-oriented model of AB fermentation by C. acetobutylicum that combines the metabolic network with the regulation of pH-induced proteome concentration (see Methods: pH-dependent Modeling of AB Fermentation), the pH value being the key factor for the metabolic changes in C. acetobutylicum. We fitted the model to experimental data of the dynamic shift experiments through the use of a switch function F(p) (see Equation (3)) to demonstrate the influence of the pH upon the fermentation products: acetate and butyrate are produced at a high pH value and acetone, butanol, and ethanol at a low pH value. The resulting numerical simulations provide a good fit to the data, see Figure 2. Thus, our model is able to explain the metabolic shift of C. acetobutylicum in continuous culture.
For the benefit of the reader, the SBML implementations of the model for each of the four experiments have been included as additional files. The first 'forward shift' experiment is encoded in Additional file 2, the second 'forward shift' experiment in Additional file 3, and the third 'forward shift' experiment in Additional File 4. Finally, the 'reverse shift' experiment is provided in Additional file 5.
The model is able to describe four independent experiments: three 'forward shift' experiments (Figure 2(a) - Figure 2(c)) and one 'reverse shift' experiment (Figure 2(d)). In the 'forward' experiments the culture was shifted from pH 5.7 to 4.5. In the first of these ('Forward' 1, Figure 2(a)), acidogenesis at pH 5.7 was maintained for 137 hours, after which the pH control was stopped, allowing the natural metabolic shift to the production of the solvents to begin. In this case, the shift took 22 hours, meaning that 159 hours after the start of the experiment, the cultures reached a pH of 4.5, where solventogenesis occured. The final measurement was taken at 215 hours. At this point, however, the continuous culture had not yet reached the steady state of this phase (in general C. acetobutylicum needs five complete volume changes to be in a stable state which would be after approximately 224 hours). Consequently, two repetitions of this experiment were performed, ensuring that both of these attained the steady state in each phase. In the first repetition of the 'forward' experiment ('Forward' 2, Figure 2(b)), the pH control was switched off after 137.5 hours, the system needed around 33.5 hours for the shift before reaching the solvent producing phase after 171 hours. Figure 2(b) illustrates that steady state was reached after roughly 236 hours. In the second repetition ('Forward' 3, Figure 2(c)), following 121 hours in acidogenesis, the shift lasted approximately 29 hours. From Figure 2(d), we see that the solventogenic steady state was reached after around 215 hours. In the 'reverse' experiment (Figure 2(d)) the culture was shifted from pH 4.5 to 5.7. In this case the pH control was switched off after the culture was kept for 129 hours in the solvent producing phase. The metabolic shift took around 17 hours before reaching the acid producing phase.
In each case the two distinct pH-dependent phases are evident: at pH 5.7 the main fermentation products are always the acids acetate and butyrate, while at pH 4.5 these are acetone and butanol. Due to its position in the pathway, ethanol is produced at low levels in both phases.
For the variables which are measured experimentally, initial conditions for each simulation are the corresponding first data points, and zero for all other variables. Importantly, only the 'forward' data were used to estimate parameters: no further fitting was undertaken to replicate the 'reverse' experiment in Figure 2(d). In all instances, the model describes well the experimental behavior. There are two principal discrepancies between the simulations and the data. Firstly, the model does not capture the rise in butyrate concentration which occurs during acidogenesis. It is possible that this is caused by our choice of initial conditions for the variables which are not measured experimentally (i.e. all the intermediate metabolites) since these will undoubtedly affect the dynamics of the system before reaching steady state (we will see in the subsequent section that the system is monostable meaning that the steady states in either phase will be achieved regardless of initial conditions). Secondly, rather than maintaining a relatively constant level of ethanol throughout the time course, the simulation demonstrates a switch in the amount of this metabolite being produced, the model attaining higher ethanol levels during the solventogenic phase (though the relative errors between the simulations and the data for ethanol concentration are comparable to those of the other end-products and, indeed, smaller than those of the other solvents). This second discrepancy may be more complex since the disagreement between the simulation and the data concerns the steady-state level in acid phase, rather than the dynamics. We focused on AdhE1, the major enzyme for butanol production, and did not consider AdhE2 in the model. However, recent experiments have shown that the transcript of the adhE2 gene is significantly upregulated during acidogenesis and antagonistically regulated with adhE1. Thus, it might be that this enzyme is responsible for ethanol production at higher pH values.
Our model suggests that a simultaneous change in regulation of the genes adc, ctfA/B, adhE is sufficient to induce the change in phenotype involving the switch between acidogenesis and solventogenesis in either direction since these are the only genes included explicitly in the model given by (5) which matches well the experimental data in Figure 2. Thus the model has helped in the characterization of the metabolic shift by identifying key genes involved in the process.
A systems biology study of genetic engineering of the metabolic network: steady-state analysis of butanol yield
In order to investigate means of enhancing butanol yield, we use the previous model to study the effects of altered expression levels of the genes associated with solventogenesis, i.e. adc, ctfA/B, adhE, bdhA/B and thlA, on the steady-state butanol concentration ( ). A number of recent experimental studies investigate the effects of altering gene regulation on product yields [43–46], but all of these consider C. acetobutylicum in batch culture, rather than continuous culture, meaning that parameters and product concentrations are not directly transferable between the two experimental setups. For instance, reaction rates can be dependent upon, amongst other things, pH, and, unlike in our experiments, the pH of a batch culture is continuously evolving (in our work pH is constant in each phase). Thus, though we will make comparisons between the qualitative results of the experimental studies and our theoretical work, it must be remembered that there are fundamental differences between the two.
In order to investigate the consequence of altering gene expression, we examine the dependence of the steady states of the system on the enzyme production rates. We display the results for butanol yield in Figure 3. The results for other metabolites can be found in Additional file 1, Figures S1 - S6.
Firstly we consider , the parameter controlling production of the enzyme Adc. Adc is responsible for the conversion of acetoacetate into acetone and the analysis shows that variation in Adc levels should have a negligible effect upon butanol production, see Figure 3(a). This is consistent with results of  which suggested that, although required for acetone production, downregulation of adc has no effect upon acetone yield. These results, however, have been recently contradicted in  where it was demonstrated that disruption of the adc gene in a hyperbutanol-producing strain of C. acetobutylicum (EA2018) could, in fact, cause a decrease in acetone production (again in batch culture). Nevertheless, the absolute amount of butanol was not significantly changed, but consequently the ratio between butanol and acetone yield was altered. The strain variation could at least partially account for the differences in the qualitative results, in addition to a different gene manipulation technique being used (TargeTron , rather than antisense RNA). However, for our model to reproduce such behavior we would need to include the reverse of the CoA-transferase-mediated reaction so that the build-up of acetoacetate caused by downregulation of adc could be consumed in the production of end-products other than acetone. To our knowledge, the activity of the reverse reaction has not been shown to be significant (and this is why it is excluded from our model) but the results of  are a hint that this reverse reaction may indeed arise in certain circumstances. There are also many more factors which should be taken into account such as the evolving pH in batch culture (and the dependence of the parameters in the model upon this pH), the effect of the genetic alterations on growth rate (the disruption in the adc gene caused a significant inhibition on growth of the EA2018 cells) and the effect on redox balance.
The CoA-transferase (CtfA/B) synthesizes the conversion of acetate and acetoacetyl-CoA into acetoacetate, acetyl-CoA and butyryl-CoA, i.e. CtfA/B regulates the conversion of acids into solvent precursors (acetoacetate for acetone, acetyl-CoA for ethanol and butyryl-CoA for butanol). Increasing causes an increase in the production of acetone and ethanol, while changes in butanol production are negligible, see Figure 3(b). In , the authors investigated downregulating ctfA/B (again in batch culture) with the expectation that this should lower acetone production and possibly therefore give rise to increased levels of butanol. However, it was found that downregulation of ctfA/B lowered levels of both acetone and butanol and this is consistent with our results, see Figure 3(c). It is speculated in  that the reduced butanol levels could be as a result of the antisense RNA causing unintended downregulation of adhE. While this is, indeed, a valid explanation, our model suggests that, at least in continuous culture, downregulation of ctfA/B alone could lower butanol levels because less butyryl-CoA is available from which butanol can be formed (see Additional file 1: Figure S2).
AdhE is required for two solvent-producing reactions, the conversion of acetyl-CoA into ethanol and of butyryl-CoA into butanol. We would therefore expect that increasing (the rate of AdhE production) would increase both ethanol and butanol output. However, Figure 3(d) illustrates that this is not necessarily the case. While increasing correlates with increased levels of ethanol, butanol levels are initially raised (coinciding with a decrease in acids) before plateauing over a range of around the value seen experimentally in the wild-type strain (the wild-type associated value of lies just to the left of the plateau). In this region of we observe a sharp rise in acetate production and a smaller one in butyrate. If is increased further, butanol levels decrease along with both acids (though the acetate levels remain high). To our knowledge, no experimental studies have been published considering the effect on product yield of adhE in isolation but the batch experiments reported in [44, 48] examine simultaneous suppression of ctfA/B and overexpression of adhE. In this case, overexpression of adhE causes an increase in all solvents (though butanol levels are not increased above wild-type levels: the adhE overexpression compensates for ctfA/B downregulation). This is consistent with the lower range of discussed above. The dual role of AdhE in solvent production obstructs the possibility of exploiting overexpression of adhE to increase butanol yield.
Thus we turn to BdhA and BdhB, which play a similar role to AdhE with respect to butanol production alone, and not of ethanol. Figure 3(e) demonstrates the effect of varying r B : butanol levels do increase with larger values of r B but the rate of increase decelerates with increasing r B , the levels saturating at a yield slightly higher than that seen in the model under default parameter values, i.e. those representing the wild-type cells.
The last enzyme that we consider is ThlA, which is required for the flow through the network downstream of acetyl-CoA. Figure 3(f) illustrates the consequences of varying this parameter. It is evident that, although a thlA mutant should be deficient in butanol, overexpressing thlA cannot push butanol yield significantly above that observed in the wild-type strain; this is consistent with experimental results in .
We are primarily concerned with genetic manipulation to maximize butanol production, this solvent being the most sought after for the manufacture of biofuels. It would be logical to predict that overexpressing the genes responsible for producing the enzymes involved in the solvent pathways of Figure 1 would increase butanol production. However, this is only true to a limited extent. The most effective way of achieving enhanced butanol yield is, according to our model, via overexpression of the bdhA/B genes, but the maximum levels of butanol attained via this method are not markedly higher than those observed in our wild-type experiments. Though a small increase in adhE transcription will induce a higher butanol yield (but, again, not much higher than that observed in the wild-type), any further increase has a detrimental effect with the increased solvent production being associated only with ethanol (shown in Additional file 1: Figure S4).
The model presented here, however, provides an efficient and simple means by which more complex approaches involving alteration to the expression of multiple genes simultaneously can be investigated. Such a theoretical approach should generate experimentally-testable hypotheses; future publications will address this.
The dependence of the metabolic phase of C. acetobutylicum on the pH of the environment is confirmed by our experimental results: at low pH the dominant fermentation products are solvents, with acids at high pH. Incorporating gene regulation into the mathematical model of the fermentation process furnishes a mechanistic representation of this pH-induced switch between the two phases. Though the model captures well the observed behavior, including additional regulation may improve the fit. For instance, we have assumed, for simplicity, that each solvent-associated gene is activated at the same threshold pH value; assigning distinct threshold values to each gene may improve the accuracy of the model.
Systems biology approaches involving alterations to gene expression levels are frequently adopted to investigate and exploit bacteria, for example in increasing the yield of valuable fermentation products specific to the bacterium of interest, which in the case of C. acetobutylicum is butanol. By incorporating gene regulation into our model we are able to perform steady-state analyses to make predictions upon the likely efficacy of potential systems biology approaches. By varying gene expression level in silico, we infer that, to increase butanol yield, alterations in the expression of single solvent-associated genes are insufficient; a more complex approach targeting two or more genes simultaneously is required.
Kinetic models such as that presented here are a useful tool for examining gene regulation. However, the number of parameters to be measured or estimated can limit the optimal size of such a model before the level of unknown information becomes too great. Accordingly, components which could affect the behavior of the model in certain circumstances can be neglected in favor of reducing the number of parameters. Genome-scale metabolic models such as  and  can take into account the overall flux through the whole cell, but they lack key regulatory information. To gain full benefit from a systems biology approach, both types of model should be used to complement the other, e.g. the feasibility of useful mutants predicted by kinetic modeling can be tested by flux balance analysis in the genome-scale models and these genome-scale models can be used to assess the most important supplementary components to be included in the kinetic model or to identify additional sinks for existing variables. Thus, drawn together, the two modeling approaches can be mutually beneficial. The model presented in this work is therefore a stepping stone from which a theoretical exploration of gene manipulation in C. acetobutylicum can be conducted to generate experimentally-testable hypotheses.
List of abbreviations
Acetoacetate decarboxylase (Adc)
Aldehyde dehydrogenase (AdhE)
CoA transferase (CtfA/B)
Limiting rate Michaelis-Menten constant .
The authors acknowledge support by the German Federal Ministry for Education and Research (BMBF) (FKZ 0313981 D, 0315782 D, 0313978F, and 0315784E) and the Biotechnology and Biological Sciences Research Council (BBSRC) (BBF 0033821 and BBI 0045131) as part of the European Transnational Network - Systems Biology of Microorganism (SysMO) - within the COSMIC and BaCell consortia. JRK also acknowledges the funding of the Royal Society and Wolfson Foundation.
- Dürre P: Biobutanol: An attractive biofuel. Biotechnol J. 2007, 2: 1525-1534.View ArticlePubMedGoogle Scholar
- Fargione J, Hill J, Tilman D, Polasky S, Hawthorne P: Land clearing and the biofuel carbon debt. Science. 2008, 319: 1235-1238. 10.1126/science.1152747View ArticlePubMedGoogle Scholar
- Searchinger T, Heimlich R, Houghton R, Dong F, Elobeid A, Fabiosa J, Tokgoz S, Hayes D, Yu T: Use of US croplands for biofuels increases greenhouse gases through emissions from land-use change. Science. 2008, 319: 1238-1240. 10.1126/science.1151861View ArticlePubMedGoogle Scholar
- Kline K, Dale V: Biofuels: effects on land and fire. Science. 2008, 321: 199- 10.1126/science.321.5886.199View ArticlePubMedGoogle Scholar
- Koh L, Tan H, Sodhi N: Biofuels: waste not want not. Science. 2008, 320: 1418-10.1126/science.320.5882.1419a.View ArticleGoogle Scholar
- Sheridan C: Making green. Nat Biotechnol. 2009, 27: 1074-1076. 10.1038/nbt1209-1074View ArticlePubMedGoogle Scholar
- Dürre P, : Handbook on clostridia. 2005, Boca Raton: CRC PressGoogle Scholar
- Woods D: The clostridia and biotechnology. 1993, Butterworth-Heinemann, BostonGoogle Scholar
- Lee S, Park J, Jang S, Nielsen L, Kim J, Jung K: Fermentative butanol production by clostridia. Biotechnol Bioeng. 2008, 101: 209-228. 10.1002/bit.22003View ArticlePubMedGoogle Scholar
- Jones D, Woods D: Acetone-Butanol Fermentation Revisited. Microbiol Rev. 1986, 50: 484-524.PubMed CentralPubMedGoogle Scholar
- Bahl H, Andersch W, Braun K, Gottschalk G: Effect of pH and butyrate concentration on the production of acetone and butanol by Clostridium acetobutylicum grown in continuous culture. Appl Microbiol Biotechnol. 1982, 14: 17-20. 10.1007/BF00507998.View ArticleGoogle Scholar
- Gottwald M, Gottschalk G: The internal pH of Clostridium acetobutylicum and its effect on the shift from acid to solvent formation. Arch Microbiol. 1985, 143: 42-46. 10.1007/BF00414766.View ArticleGoogle Scholar
- Segel I: Enzyme Kinetics: Behavior and Analysis of Rapid Equilibrium and Steady-State Enzyme Systems. 1993, New York: Wiley-InterscienceGoogle Scholar
- Gerischer U, Dürre P: Cloning, sequencing, and molecular analysis of the acetoacetate decarboxylase gene region from Clostridium acetobutylicum. J Bacteriol. 1990, 172: 6907-6918.PubMed CentralPubMedGoogle Scholar
- Petersen D, Bennett G: Purification of acetoacetate decarboxylase from Clostridium acetobutylicum ATCC 824 and cloning of the acetoacetate decarboxylase gene in Escherichia coli. Appl Environ Microb. 1990, 56: 3491-3498.Google Scholar
- Janssen H, Döring C, Ehrenreich A, Voigt B, Hecker M, Bahl H, Fischer RJ: A proteomic and transcriptional view of acidogenic and solventogenic steady state cells of Clostridium acetobutylicum in a chemostat culture. Appl Microbiol Biotechnol. 2010, 87: 2209-2226. 10.1007/s00253-010-2741-xPubMed CentralView ArticlePubMedGoogle Scholar
- Bahl H, Andersch W, Gottschalk G: Continuous production of acetone and butanol by Clostridium acetobutylicum in a two-stage phosphate limited chemostat. Appl Microbiol Biotechnol. 1982, 15: 201-205. 10.1007/BF00499955.View ArticleGoogle Scholar
- Fischer RJ, Oehmcke S, Meyer U, Mix M, Schwarz K, Fiedler T, Bahl H: Transcription of the pst operon of Clostridium acetobutylicum is dependent on phosphate concentration and pH. J Bacteriol. 2006, 188: 5469-5478. 10.1128/JB.00491-06PubMed CentralView ArticlePubMedGoogle Scholar
- Grupe H, Gottschalk G: Physiological events in Clostridium acetobutylicum during the shift from acidogenesis to solventogenesis in continuous culture and presentation of a model for shift induction. Appl Environ Mircobiol. 1992, 58: 3896-3902.Google Scholar
- Meyer C, Papoutsakis E: Increased levels of ATP and NADH are associated with increased solvent production in continuous culture of Clostridium acetobutylicum. Appl Microbiol Biotechnol. 1989, 30: 450-459. 10.1007/BF00263849.View ArticleGoogle Scholar
- Jones D, van der Westhuizen A, Long S, Allcock E, Reid S, Woods D: Solvent production and morphological changes in Clostridium acetobutylicum. Appl Environ Microb. 1982, 43: 1434-1439.Google Scholar
- Fiedler T, Mix M, Meyer U, Mikkat S, Glocker M, Bahl H, Fischer RJ: The two-component system PhoPR of Clostridium acetobutylicum is involved in phosphate-dependent gene regulation. J Bacteriol. 2008, 190: 6559-6567. 10.1128/JB.00574-08PubMed CentralView ArticlePubMedGoogle Scholar
- Papoutsakis E: Equations and calculations for fermentations of butyric acid bacteria. Biotechnol Bioeng. 1984, 26: 174-187. 10.1002/bit.260260210View ArticlePubMedGoogle Scholar
- Votruba J, Volesky B, Yerushalmi L: Mathematical model of a batch acetone-butanol fermentation. Biotechnol Bioeng. 1986, 28: 247-255. 10.1002/bit.260280215View ArticlePubMedGoogle Scholar
- Desai R, Harris L, Welker N, Papoutsakis E: 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 ArticlePubMedGoogle Scholar
- Senger R, Papoutsakis E: Genome-scale model for Clostridium acetobutylicum: Part I. Metabolic network resolution and analysis. Biotechnol Bioeng. 2008, 101: 1036-1052. 10.1002/bit.22010PubMed CentralView ArticlePubMedGoogle Scholar
- Senger R, Papoutsakis E: 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.22009PubMed CentralView ArticlePubMedGoogle Scholar
- Shinto H, Tashiro Y, Yamashita M, Kobayashi G, Sekiguchi T, Hanai T, Kuriya Y, Okamoto M, Sonomoto K: Kinetic modeling and sensitivity analysis of acetone-butanol-ethanol production. J Biotechnol. 2007, 131: 45-56. 10.1016/j.jbiotec.2007.05.005View ArticlePubMedGoogle Scholar
- Li RD, Li YY, Lu LY, Ren C, Li YX, Liu L: An improved kinetic model for acetone-butanol-ethanol pathway of Clostridium acetobutylicum and model-based perturbation analysis. The Fourth Internation Conference on Computational Systems Biology. 2010, 62-71.Google Scholar
- Fischer RJ, Helms J, Dürre P: Cloning, sequencing, and molecular analysis of the sol operon of Clostridium acetobutylicum, a chromosomal locus involved in solventogenesis. J Bacteriol. 1993, 175: 6959-6969.PubMed CentralPubMedGoogle Scholar
- Nair R, Bennett G, Papoutsakis E: Molecular characterization of an aldehyde/alcohol dehydrogenase gene from Clostridium acetobutylicum ATCC 824. J Bacteriol. 1994, 176: 871-885.PubMed CentralPubMedGoogle Scholar
- Thormann K, Feustel L, Lorenz K, Nakotte S, Dürre P: Control of butanol formation in Clostridium acetobutylicum by transcriptional activation. J Bacteriol. 2002, 184: 1966-1973. 10.1128/JB.184.7.1966-1973.2002PubMed CentralView ArticlePubMedGoogle Scholar
- Walter K, Bennett G, Papoutsakis E: Molecular characterization of two Clostridium acetobutylicum ATCC 824 butanol dehydrogenase isozyme genes. J Bacteriol. 1992, 174: 7149-7158.PubMed CentralPubMedGoogle Scholar
- Winzer K, Lorenz K, Zickner B, Dürre P: Differential regulation of two thiolase genes from Clostridium acetobutylicum DSM 792. J Mol Microbiol Biotechnol. 2000, 2: 531-541.PubMedGoogle Scholar
- Dürre P, Fischer RJ, Kuhn A, Lorenz K, Schreiber W, Stürzenhofecker B, Ullmann S, Winzer K, Sauer U: Solventogenic enzymes of Clostridium acetobutylicum: catalytic properties, genetic organization, and transcriptional regulation. FEMS Microbiol Rev. 1995, 17: 251-262.PubMedGoogle Scholar
- Grimmler C, Janssen H, Krauβe D, Fischer RJ, Bahl H, Dürre P, Liebl W, Ehrenreich A: Genome-wide gene expression analysis of the switch between acidogenesis and solventogenesis in continuous cultures of Clostridium acetobutylicum. J Mol Microbiol Biotechnol. 2011, 20: 1-15. 10.1159/000320973View ArticlePubMedGoogle Scholar
- Feustel L, Nakotte S, Dürre P: Characterization and development of two reporter gene systems for Clostridium acetobutylicum. Appl Environ Microbiol. 2004, 70: 798-803. 10.1128/AEM.70.2.798-803.2004PubMed CentralView ArticlePubMedGoogle Scholar
- Harris L, Welker N, Papoutsakis E: Northern, morphological, and fermentation analysis of spo0A inactivation and overexpression in Clostridium acetobutylicum ATCC 824. J Bacteriol. 2002, 184: 3586-3597. 10.1128/JB.184.13.3586-3597.2002PubMed CentralView ArticlePubMedGoogle Scholar
- Cornish-Bowden A: Fundamentals of Enzyme Kinetics. 2004, Colchester: Portland Press, Third revisedGoogle Scholar
- Millat T, Bullinger E, Rohwer J, Wolkenhauer O: Approximations and their consequences for dynamic modelling of signal transduction pathways. Math Biosci. 2007, 207: 40-57. 10.1016/j.mbs.2006.08.012View ArticlePubMedGoogle Scholar
- MathWorks (Nattick, MA, USA): MATLAB - The Language Of Technical Computing. 2010, http://www.mathworks.comGoogle Scholar
- Schmidt H, Jirstrand M: Systems Biology Toolbox for MATLAB: a computational platform for research in systems biology. Bioinformatics. 2006, 22: 514-515. 10.1093/bioinformatics/bti799View ArticlePubMedGoogle Scholar
- Sillers R, Al-Hinai M, Papoutsakis E: 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 ArticlePubMedGoogle Scholar
- Tummala S, Junne S, Papoutsakis E: 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.2003PubMed CentralView ArticlePubMedGoogle Scholar
- Tummala S, Welker N, Papoutsakis E: 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.2003PubMed CentralView ArticlePubMedGoogle Scholar
- Jiang Y, Xu C, Dong F, Yang Y, Jiang W, Yang S: Disruption of the acetoacetate decarboxylase gene in solvent-producing Clostridium acetobutylicum increases the butanol ratio. Met Eng. 2009, 11: 284-291. 10.1016/j.ymben.2009.06.002.View ArticleGoogle Scholar
- Shao L: Targeted gene disruption by use of a group II intron (targetron) vector in Clostridium acetobutylicum. Cell Res. 2007, 17: 963-965. 10.1038/cr.2007.91View ArticlePubMedGoogle Scholar
- Tummala S, Junne S, Paredes C, Papoutsakis E: Transcriptional analysis of product-concentration driven changes in cellular programs of recombinant Clostridium acetobutylicum strains. Biotechnol Bioeng. 2003, 84: 842-854. 10.1002/bit.10851View ArticlePubMedGoogle Scholar