Genome-scale modeling using flux ratio constraints to enable metabolic engineering of clostridial metabolism in silico

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 [21], (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 [9]. 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 [39]. 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[12]. 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 [20].

Engineering clostridial metabolism

Knowledge of the metabolic pathways of butanol fermentation has allowed for targeted engineering approaches. As acids (acetate and butyrate) and alcohols (ethanol and butanol) are competing products of fermentation, metabolic engineering strategies have been designed to silence acid producing pathways in attempt to re-direct carbon flow into the alcohol producing pathways. The primary metabolic network of C. acetobutylicum is shown in Figure 1 (adapted from [6, 8]). An example of this strategy is the knockdown of the butyrate kinase (buk) (BK in Figure 1) and phosphotransacetylase (pta) (PTA in Figure 1) genes in clostridial metabolism [40]. The proximity of the pta gene in the genome to the acetate kinase (ak) (AK in Figure 1) gene resulted in the silencing of both genes simultaneously, further decreasing acetate production. The results of this engineering strategy showed that eliminating genes of the acetate pathway had little effect compared to the wild-type, while knockdown of butyrate pathway genes resulted in 10% more butanol and 50% less acetone than the wild-type [41]. The aldehyde/alcohol dehydrogenase gene (aad or adhE1) (AAD in Figure 1) was over-expressed in the presence of a buk knockout, and this strain yielded a 300% increase in butanol production and 400% increase in ethanol production over the wild-type strain [42]. In a separate metabolic engineering strategy, the aad was over-expressed while knocking-down the gene for subunit B of acetoacetyl-CoA transferase (ctfB) (CoAT in Figure 1). This led to a strain with similar butanol productivity but the ability to produce extraordinary ethanol concentrations of 200 mM (23-fold higher than the wild-type) [43]. However, when the aad gene was put under control of the ptb gene promoter (to increase expression during the early acidogenic phase of the culture) an increase in butanol concentrations to 300 mM (a record high) was observed along with faster accumulation of butanol in the culture [44].

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 [36], 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 [38] 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].

Genome-scale model

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 [51] 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 [56]. 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 [57]. 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 [12] 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 $\frac{mmol{H}^{+}}{hr·gDCW}$ (proton efflux) and 5 $\frac{mmol{H}^{+}}{hr·gDCW}$ (proton influx), and the limits were chosen from experimental observations.

FBrAtio algorithm

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.

Simulations performed

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 [50]. 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 [43]. 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 $\frac{mmol}{hr·gDCW}$ and the SPF was varied between −30 and 5 $\frac{mmol{H}^{+}}{hr·gDCW}$. Only thermodynamic reversibility constraints were used initially. Results showed acidogenic and solventogenic metabolic phases that coincided with SPF values [12]. Results also showed a maximum specific growth rate at an SPF value of −10 $\frac{mmol{H}^{+}}{hr·gDCW}$, 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 (H₂) 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

Obvious futile cycles allowing artificial NAD(P)⁺/NAD(P)H balancing were located and corrected, but the problems described above remained. Given the ability of the metabolic network to artificially balance NAD(P)⁺/NAD(P)H without using the acid/solvent production pathways, additional constraints and flux ratio constraints were implemented. First, the reactions involving the ferredoxins were further constrained for irreversibility. These reactions and their updated constraints are given in Table 1. Simulation of metabolism, given these and glucose uptake and SPF constraints are shown in Figure 2. In this simulation, the following results do not coincide with experimental observations: (i) succinate was produced in high levels during acidogenesis, (ii) the production of acetate and lactate far exceed butyrate production, (iii) only butyrate was re-consumed (negative flux values), (iv) the amount of butyrate re-consumed exceeded the amount of butyrate produced, (v) butanol was produced well before acetone, and (vi) hydrogen (H₂) production fell to zero and rose during solventogenesis. These results show that during early production of butanol (before production of acetone), butyrate was consumed through the production pathway. This is consistent with previous findings [7, 58], and the thermodynamic reversibility calculations [54] insist the reactions catalyzed by the butyrate kinase (BK) and the phosphotransbutyrylase (PTB) (see Figure 1) remain reversible.

Enzyme Name	Reaction	Lower Bound	Upper Bound
PFO: pyruvate ferredoxin oxidoreductase	$\begin{array}{l}CoA+pyruvate+F{d}_{\mathit{Ox}}\to \\ C{O}_2+acetyl-CoA+H^{+}+F{d}_{\mathit{Red}}\end{array}$	0	1000
HYDA: hydrogenase	$2H^{+}+F{d}_{\mathit{Red}}\to H_2+F{d}_{\mathit{Ox}}$	0	1000
FNO: ferredoxin NAD+ oxidoreducatase	$\begin{array}{l}NA{D}^{+}+H^{+}+F{d}_{\mathit{Red}}\to \\ NADH+F{d}_{\mathit{Ox}}\end{array}$	0	1000
FNPO: ferredoxin NADP+ oxidoreductase	$\begin{array}{l}NAD{P}^{+}+H^{+}+F{d}_{\mathit{Red}}\to \\ NADPH+F{d}_{\mathit{Ox}}\end{array}$	0	1000

*The reduced ferredoxin (Fd_Red) and the oxidized ferredoxin (Fd_Ox).

This flux ratio constraint corrected hydrogen production (i.e., high hydrogen production during acidogenesis and reduced during solventogenesis) and minimized succinate secretion, which is consistent with experimental observations. This constraint also resulted in increased CO₂ production during solventogenesis, relative to production during acidogenesis. However, metabolic activity was still characterized by (i) high acetate and ethanol production and (ii) low butyrate and butanol production (results not shown). To address this, flux ratios were constructed around the use of the acetyl-CoA metabolite pool in clostridial metabolism (see Figure 1). Acetyl-CoA can be utilized to (i) produce ATP through acetate production, (ii) regenerate NAD⁺ through ethanol production, or (iii) balance both ATP and NAD⁺ by producing butyrate. However, the shorter metabolic pathways result in acetate and ethanol production and accommodate minimizing the total flux of the metabolic network. To approximate wild-type clostridial metabolism, flux must proceed from acetyl-CoA through the thiolase (THL) with greater flux than through the phosphotransacetylase (PTA) towards acetate or through the bifunctional aldehyde/alcohol dehydrogenase (AAD). To ensure this, additional flux ratio constraints were derived. The ratio of metabolic flux through the THL relative to the PTA, f(rTHL):f(rPTA), was set equal to 2, and the flux through the THL relative to the AAD, f(rTHL):f(rAAD_1), was set equal to 10. This ensured the majority of acetyl-CoA was sent to the butyrate producing pathway while a greater amount of acetyl-CoA was converted to acetate than was converted to ethanol. These flux ratios are shown in Equations 8 and 9. Simulation results are shown in Figure 3. Due to the constant glucose uptake rate, solutions to the flux balance equation were only possible for SPF values greater than −17 $\frac{mmol{H}^{+}}{hr·gDCW}$. In a previous study [12], much higher glucose uptake rates enabled flux solutions at SPF values approaching −55 $\frac{mmol{H}^{+}}{hr·gDCW}$. In the simulations shown in Figure 3, a roughly 1:1 production ratio of H₂/CO₂ was observed during exponential growth, and the production of H₂ decreased in solventogenesis, while the production of CO₂ increased. This is consistent with experimental observations [59]. Butyrate was produced in greater amounts than acetate, and both acids were taken up during solventogenesis. The production of butanol was greater than that of ethanol and was similar to that of acetone. The production of lactate was minimal, as was the secretion of succinate. The maximum specific growth rate of the culture occurred towards the end of acidogenesis but prior to solventogenesis onset. Thus, when the flux ratio constraints of Equations 59 were applied through FBrAtio, the i CAC490 genome-scale model was able to capture the major properties of wild-type metabolism. This was done for the first time without constraining acids and solvents production rates directly using constraints on transporters or exchange fluxes.

$f\left(rTHL\right):f\left(rPTA\right)=\frac{f\left(rTHL\right)}{f\left(rPTA\right)}=2$

(8)

$f\left(rTHL\right):f\left(rAAD\_1\right)=\frac{f\left(rTHL\right)}{f\left(rAAD\_1\right)}=10$

(9)

asRNA knockdown of CoAT only

Previous research found that knockdown of the CoA transferase using asRNA technology resulted in a more favorable butanol to acetone selectivity. Although, overall reduced butanol yields were observed [50], possibly due to simultaneous AAD knockdown. To determine if FBrAtio can predict of these findings, a flux ratio constraint was derived to simulate knockdown of the CoAT. It was assumed this knockdown led to a decreased flux of acetate and butyrate re-uptake. As shown in Figure 1, the CoAT converts acetoacetyl-CoA to acetoacetate while transporting acetate or butyrate into the cell and converting them to acetyl-CoA or butyryl-CoA, respectively. Acetoacetyl-CoA, on the other hand, can also be converted to β-hydroxybutyryl-CoA by the β-hydroxybutyryl-CoA dehydrogenase (BHBD). So, a flux ratio constraint was derived to specify the distribution of β-hydroxybutyryl-CoA utilized by the CoAT and the BHBD. This flux ratio is called f(rCoAT): f(rBHBD), where the flux through the CoAT includes fluxes for both acetate and butyrate re-uptake. In the wild-type simulations shown in Figure 3, given an SPF of 5 $\frac{mmol{H}^{+}}{hr·gDCW}$, the flux through the CoAT was 2.03 $\frac{mmol}{hr·gDCW}$ for acetate re-uptake and 1.28 $\frac{mmol}{hr·gDCW}$ for butyrate re-uptake (total CoAT flux of 3.31 $\frac{mmol}{hr·gDCW}$). The wild-type flux through BHBD was 3.04 $\frac{mmol}{hr·gDCW}$ to yield a wild-type f(rCoAT): f(rBHBD) flux ratio of 1.09. The flux ratio is given in Equations 10. However, this flux ratio in the wild-type was a function of the SPF, and this relationship is shown in Figure 4. To simulate knockdown of the CoAT, two separate approaches were attempted. The first approach continued to treat this flux ratio as a function of the SPF. The flux ratio was constrained to values of (i) 75%, (ii) 50%, and (iii) 25% of the wild-type value, as shown in Figure 4. The second approach was to fix the flux ratio to a specified value for all values of the SPF. The flux ratios chosen were: (i) 1, (ii) 0.5, (iii) 0.1, and (iv) 0.01.

$\frac{f\left(rCoAT\right)}{f\left(rBHBD\right)}=1.09$

(10)

FBrAtio results with the f(rCoAT):f(rBHBD) flux ratios of Figure 4 are shown in Figure 5. In particular, results of acetone, butanol, and ethanol predictions are shown to focus on the impact of CoAT asRNA knockdown on solvent production. The SPF range shown is from −10 to 5 $\frac{mmol{H}^{+}}{hr·gDCW}$. Consistent with published experimental results [43, 50], simulated knockdown of the CoAT resulted in an increased butanol to acetone ratio. Significantly lower acetone production (~50% reduction) was predicted with increased butanol production (~25% increase) when CoAT activity was down-regulated by 75%. Significant about these simulation results is that the impact of “fine-tuning” an asRNA construct can be observed through simulations using flux ratio constraints.

It was important to determine whether the f(rCoAT):f(rBHBD) flux ratio must be represented as a function of the SPF (Figure 4). In further simulations, this flux ratio constraint was artificially set and held constant for all values of the SPF. FBrAtio results are shown in Figure 6 for flux ratios of: (i) 1, (ii) 0.5, (iii) 0.1, and (iv) 0.01. Similar trends were obtained. As the f(rCoAT): f(rBHBD) flux ratio decreased, metabolic flux was forced through the BHBD enzyme (rather than through CoAT), resulting in (i) decreased acetate/butyrate re-uptake, (ii) decreased acetone production, and (iii) increased flux through the butanol production pathway. The exaggerated flux ratio simulations of Figure 6 (ratios of 0.1 and 0.01) show the potential of effective asRNA or gene knockout of the CoAT. However, even though CoAT knockdown simulations resulted in predicted phenotypes with improved butanol to acetone selectivity, these phenotypes did not show overall reduced ethanol and butanol production. This suggests AAD knockdown may be required to impact alcohol yields, which is discussed in the next section.

Knockdown of both CoAT and AAD

In published research [43, 50], it was suspected that the aad gene was knocked-down when the asRNA for cftB (CoAT) was applied. The resulting strain produced very little ethanol or butanol. Since these two genes reside in the same operon, this hypothesis is valid. The previous results (Figures 5 and 6) showed that the butanol to acetone ratio was increased with CoAT knockdown. Simulations of the knockdown of both CoAT and AAD to defined levels using FBrAtio are shown in Figure 7. The AAD was knocked-down by the use of two flux ratio constraints. In these simulations, knockdown of AAD must be accommodated at two metabolic branch points (i) acetyl-CoA and (ii) butyryl-CoA (see Figure 1). To simulate the AAD knockdown, the ratio of flux diverted through AAD, relative to its other choices, was decreased by 80%. This factor was chosen based on the published effectiveness of the originally designed asRNA [50]. To incorporate the knockdown into flux ratio constraints, the f(rTHL):f(rAAD_1) flux ratio constraint was increased from a value of 10 to 50. To implement the flux ratio constraint at butyryl-CoA, a new f(rPTB):f(rAAD_2) flux ratio constraint was created. In wild-type simulations (see Figure 3), this flux ratio was unconstrained and had an average value of 0.52 at an SPF value of 0 $\frac{mmol{H}^{+}}{hr·gDCW}$ and 0 at an SPF value of 5 $\frac{mmol{H}^{+}}{hr·gDCW}$ (assumed value of 0.25 during solventogenesis). At highly negative values of the SPF, this ratio becomes infinite (since no butanol is produced). The f(rPTB):f(rAAD_2) was constrained to 1.25 (80% knockdown over 0.25) over all values of the SPF. The f(rTHL):f(rAAD_1) flux ratio constraint was held constant at 50 (also 80% knockdown). The CoAT was knocked-down 80% by adjusting the f(rCoAT):f(rBHBD) flux ratio to 0.2 over all values of the SPF. Results for solventogenesis (SPF = 5) are shown in Figure 7. Initial simulations revealed a large ethanol production that was not consistent with experimental findings. A closer inspection of genome-wide metabolic fluxes revealed several reactions were involved in futile cycles to produce excess acetaldehyde. To correct this, a new flux ratio constraint was created between AAD_1 and AAD_3 and is shown in Equation 11. This new flux ratio constraint was set to 1.2 since it is known that acetaldehyde is not produced exclusively through AAD_1 (although the exact ratio has not been measured). FBrAtio results given the f(rAAD_3):f(rAAD_1) flux ratio constraint yielded significantly reduced ethanol production. Results with 80% knockdown of CoAT and AAD still showed significant butanol production; however, this was significantly reduced (towards zero) as the percent knockdown was increased.

$\frac{f(rAAD\_3)}{f(rAAD\_1)}=1.2$

(11)

Over-express AAD and knockdown CoAT

The AAD was over-expressed using the f(rTHL):f(rAAD_1) and f(rPTB):f(rAAD_2) flux ratio constraints. In this case, the f(rTHL):f(rAAD_1) flux ratio constraint was adjusted from 10 (wild-type) to values of (i) 5, (ii) 2.5, and (iii) 0.3125 to simulate over-expression (over wild-type levels) by (i) 100%, (ii) 200%, and (iii) 500% respectively. The f(rPTB):f(rAAD_2) flux ratio constraint was also adjusted accordingly to simulate these levels of over-expression. FBrAtio results of solventogenesis (SPF = 5) are shown in Figure 8. The over-expression of AAD under control of its native promoter led to large increases in ethanol production in experimental observations [43]. However, in this simulation study, an increase in AAD expression of 500% was required to see this dramatic increase. As AAD over-expression reached 1000% (results not shown), ethanol production was increased an additional 60% (relative to the 500% AAD over-expression level). This result is possibly explained by the presence of multiple copies of the plasmid present in the cell during the experimental trials.