BACKGROUND
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  <dig> 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:  acetate,  butyrate,  lactate,  acetone,  ethanol, and several others in small amounts  <cit> .
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  <cit> .
multiple “genome-scale” metabolic network reconstructions now exist for c. acetobutylicum <cit> .
similar networks and their corresponding genome-scale models have been reviewed extensively  <cit> .
in general, they are used to  complete genome annotation  <cit> ,  predict optimal culturing conditions  <cit> ,  discover genomic regulation  <cit> ,  identify essential genes and drug targets  <cit> ,  study strain evolution  <cit> , and  design productive strains  <cit> .
modeling results on the genome-scale have been applied to both “acidogenic” and “solventogenic” programs of clostridial metabolism  <cit> .
the acidogenic program is characterized by high acids  production and high growth rates, and the solventogenic program  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  <cit> .
several insights into c. acetobutylicum metabolism have been gained from “gap filling” the metabolic network by locating previously unknown enzymes and biochemical reactions  <cit> .
the total rate at which a cell produces/consumes protons through the several membrane transport mechanisms is termed the specific proton flux , 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 <cit> .
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  <cit> .
engineering clostridial metabolism
knowledge of the metabolic pathways of butanol fermentation has allowed for targeted engineering approaches.
as acids  and alcohols  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  <dig> .
an example of this strategy is the knockdown of the butyrate kinase   and phosphotransacetylase   genes in clostridial metabolism  <cit> .
the proximity of the pta gene in the genome to the acetate kinase   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  <cit> .
the aldehyde/alcohol dehydrogenase gene   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  <cit> .
in a separate metabolic engineering strategy, the aad was over-expressed while knocking-down the gene for subunit b of acetoacetyl-coa transferase  .
this led to a strain with similar butanol productivity but the ability to produce extraordinary ethanol concentrations of 200 mm   <cit> .
however, when the aad gene was put under control of the ptb gene promoter  an increase in butanol concentrations to 300 mm  was observed along with faster accumulation of butanol in the culture  <cit> .
metabolic engineering in silico
the goal of metabolic engineering in silico is to derive  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 <cit> , 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  lactate,  ethanol, or  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  <cit> , 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  <cit>  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   <cit> .
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 <cit>  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  <cit> .
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  metabolic engineering strategies applied to clostridia , it was not immediately apparent which design 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 “flux balance analysis with flux ratios ” 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 .
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  <cit> .
methods
genome-scale model
a new genome-scale model for c. acetobutylicum atcc  <dig> was constructed by expanding the previously published model by senger and papoutsakis  <cit> .
the new model is called icac <dig> and contains  <dig> metabolites involved in  <dig> biochemical reactions, including  <dig> membrane transport reactions.
the model includes  <dig> genes from the c. acetobutylicum genome.
the newly updated icac <dig> model differs from the original senger and papoutsakis model  <cit>  in that it contains  <dig> more reactions  and  <dig> more metabolites .
the new reactions added to create the icac <dig> model were obtained from the kegg database  <cit>  and recent literature.
the icac <dig> model also contains an updated tca cycle that operates in both oxidative and reductive directions to succinate, as shown by recent fluxomics studies  <cit> .
the model allows the export of succinate since its metabolic fate has not yet been resolved conclusively.
the icac <dig> 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  <cit>  have also been applied.
the biomass equation was also updated for the icac <dig> model, using the initial version by senger and papoutsakis  <cit>  as a template.
a nonlinear optimization procedure was applied  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 icac <dig> model is included as additional file  <dig>  the sbml formatted model is included as additional file  <dig>
flux balance analysis
the icac <dig> genome-scale model was simulated using flux balance analysis through the cobra toolbox  <cit> .
the open-source glpk linear programming software was used to solve the flux balance equation , 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  <cit> .
the objective functions used for all simulations were  maximizing the specific growth rate of the cell while  minimizing the total flux of the system.
the specific proton flux
the concept of the specific proton flux  was first introduced by senger and papoutsakis  <cit>  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  and turns slightly positive during the stationary phase .
in this research, the spf was constrained to specific values by constraining the proton exchange reaction .
the spf range was between − <dig> mmolh+hr·gdcw  and  <dig> mmolh+hr·gdcw , and the limits were chosen from experimental observations.
fbratio algorithm
metabolic engineering in silico was enabled through the application of flux ratio constraints.
in the fbratio method developed in this research, flux ratio constraints were incorporated into the stoichiometric matrix directly, enabling the flux balancing problem to be solved by simple linear programming.
several flux ratios were investigated, and a critical metabolic branch point was identified around the use of the acetyl-coa metabolite pool.
in clostridial metabolism, acetyl-coa can be used to produce multiple acids, solvents, or macromolecules to produce biomass.
the routes taken by acetyl-coa either regenerate or consume different amounts of atp and nad+.
ultimately, the balancing of these cofactors determines the production of acids and solvents.
the large number of degrees of freedom associated with the genome-scale model of clostridial metabolism allows a large number of acid/solvent production combinations that satisfy cofactor balancing and the overall mass balance, while satisfying the objective functions of the optimization.
thus, flux ratios ultimately reduce the number of degrees of freedom of the system and can be used to define selectivity.
the following example demonstrates the application of a flux ratio related to the consumption of acetyl-coa and the incorporation of this flux ratio into the stoichiometric matrix.
as shown in figure  <dig>  two possibilities for acetyl-coa are  usage by the thiolase enzyme  for conversion to acetoacetyl-coa and  usage by the phosphotransacetylase enzyme  for conversion to acetyl phosphate.
the reactions catalyzed by these enzymes are given equations.
<dig> and  <dig>
2acetyl−coa→thlacetoacetyl−coa+coa
acetyl−coa+orthophosphate→ptaacetyl-phophate+coa
next, a ratio of fluxes for these reactions is assumed.
for this hypothetical example, it is assumed that twice as much flux proceeds through the thl reaction  than the pta reaction .
ultimately, biochemical origins of differing fluxes through competing reactions, where known, can be used to calculate flux ratios.
the flux through the reaction catalyzed by thl is f and the flux through the pta catalyzed reaction is f. this flux ratio is called f:f and is represented as equation  <dig>
frthl:frpta=frthlfrpta= <dig>
to build this flux ratio constraint into the stoichiometric matrix , first equation  <dig> is rearranged to the following.
frthl−2frpta= <dig>
next, a new row is added to the stoichiometric matrix.
in this new row, two values are added .
in the column representing the reaction catalyzed by thl , the coefficient  <dig> is added to the matrix.
in the column representing the reaction catalyzed by the pta, the coefficient − <dig> is added .
with these additions, when the flux balance equation  is solved, the ratio of fluxes for the reactions of equation  <dig> and equation  <dig> will be exactly  <dig>  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 <cit> .
in particular, simulations were performed with the icac <dig> 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:  at highly negative values of spf , acetate and butyrate are produced in high quantities,  high hydrogen production accompanies acidogenesis,  solvents are produced at spf values close to zero and slightly positive,  hydrogen production decreases during solventogenesis,  the maximum growth rate of the culture occurs during acidogenesis,  the production of butyrate is slightly greater than the production of acetate and much greater than the production of lactate, and  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  <cit> .
knockdown of the coat  resulted in increased butanol to acetone selectivity, but this strategy resulted in decreased ethanol and butanol production  <cit> .
the asrna was designed against the mrna of the cftb gene in particular, which is a part of the tricistronic operon .
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  <cit> .
flux ratios were designed to  knockdown coat activity only,  knockdown activity of both coat and aad, and  knockdown coat while over-expressing aad at and above wild-type levels.
RESULTS
simulations with a minimal set of constraints
the icac <dig> model was simulated with a glucose uptake rate constrained to  <dig> mmolhr·gdcw and the spf was varied between − <dig> and  <dig> mmolh+hr·gdcw.
only thermodynamic reversibility constraints were used initially.
results showed acidogenic and solventogenic metabolic phases that coincided with spf values  <cit> .
results also showed a maximum specific growth rate at an spf value of − <dig> mmolh+hr·gdcw, which is consistent with previous findings  <cit> .
however, during acidogenesis, acetate was the primary acid produced, and acetone was the primary solvent produced during solventogenesis.
hydrogen  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+/nadh 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  and nad+  regenerating pathways.
approximating wild-type metabolism with fbratio
obvious futile cycles allowing artificial nad+/nadh balancing were located and corrected, but the problems described above remained.
given the ability of the metabolic network to artificially balance nad+/nadh 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  <dig>  simulation of metabolism, given these and glucose uptake and spf constraints are shown in figure  <dig>  in this simulation, the following results do not coincide with experimental observations:  succinate was produced in high levels during acidogenesis,  the production of acetate and lactate far exceed butyrate production,  only butyrate was re-consumed ,  the amount of butyrate re-consumed exceeded the amount of butyrate produced,  butanol was produced well before acetone, and  hydrogen  production fell to zero and rose during solventogenesis.
these results show that during early production of butanol , butyrate was consumed through the production pathway.
this is consistent with previous findings  <cit> , and the thermodynamic reversibility calculations  <cit>  insist the reactions catalyzed by the butyrate kinase  and the phosphotransbutyrylase   remain reversible.
*the reduced ferredoxin  and the oxidized ferredoxin .
a total of five flux ratio constraints were found necessary to generate simulations of metabolism that were qualitatively consistent with experimental observations.
first, the butyrate to acetate uptake ratio was constrained according to previously published findings  <cit> .
the assumption was made that the extracellular butyrate concentration was twice that of the extracellular acetate concentration.
from the original relationship  <cit> , this meant that the flux of butyrate uptake through the coat, f, to acetate uptake, f, was equal to  <dig> , as shown in equation  <dig>  this flux ratio is referred to as f:f.
frcoat,butyrate:frcoat,acetate=fcoat,butyratefcoat,acetate= <dig>
implementing this flux ratio alone led to high initial production of succinate, and hydrogen production increased during solventogenesis.
over-production of ethanol and lactate was also observed.
so, a flux ratio was installed to direct the conversion of pyruvate to either lactate, through the lactate dehydrogenase , or to acetyl-coa through the pyruvate ferredoxin oxidoreductase  .
this flux ratio is shown as equation  <dig>  this ratio was set equal to  <dig> to coincide with published experimental observations that lactate production is much less than acetate and butyrate production  <cit> .
frpfofrldh= <dig>
simulations with these two flux ratios were characterized by  high acetate and ethanol production,  low butyrate and butanol production,  high initial secretion of succinate, and  low initial production of hydrogen.
the secretion of succinate was investigated next.
this is indicative of high fluxes through the tca cycle.
no significant succinate export has been reported for c. acetobutylicum, so an optimized metabolic model must produce only minimal succinate .
in the icac <dig> model, production of oxaloacetate from pyruvate requires hco- <dig>  when co <dig> is produced, it is either  transported out of the cell or  converted to hco- <dig>  the fate of co <dig> is determined by physicochemical properties of the intracellular environment and has a significant impact on intracellular metabolism.
thus, a ratio constraint for co <dig> export against conversion was derived.
this ratio is shown in equation  <dig> and was set equal to  <dig> to approximate intracellular conditions.
this value was chosen because it led to effective simulations.
the physicochemical nature of this flux ratio constraint is currently under investigation.
fco2exportfco2conversion= <dig>
this flux ratio constraint corrected hydrogen production  and minimized succinate secretion, which is consistent with experimental observations.
this constraint also resulted in increased co <dig> production during solventogenesis, relative to production during acidogenesis.
however, metabolic activity was still characterized by  high acetate and ethanol production and  low butyrate and butanol production .
to address this, flux ratios were constructed around the use of the acetyl-coa metabolite pool in clostridial metabolism .
acetyl-coa can be utilized to  produce atp through acetate production,  regenerate nad+ through ethanol production, or  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  with greater flux than through the phosphotransacetylase  towards acetate or through the bifunctional aldehyde/alcohol dehydrogenase .
to ensure this, additional flux ratio constraints were derived.
the ratio of metabolic flux through the thl relative to the pta, f:f, was set equal to  <dig>  and the flux through the thl relative to the aad, f:f, was set equal to  <dig>  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  <dig> and  <dig>  simulation results are shown in figure  <dig>  due to the constant glucose uptake rate, solutions to the flux balance equation were only possible for spf values greater than − <dig> mmolh+hr·gdcw.
in a previous study  <cit> , much higher glucose uptake rates enabled flux solutions at spf values approaching − <dig> mmolh+hr·gdcw.
in the simulations shown in figure  <dig>  a roughly 1: <dig> production ratio of h2/co <dig> was observed during exponential growth, and the production of h <dig> decreased in solventogenesis, while the production of co <dig> increased.
this is consistent with experimental observations  <cit> .
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  <dig> were applied through fbratio, the icac <dig> 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.
frthl:frpta=frthlfrpta= <dig>
frthl:fraad_1=frthlfraad_1= <dig>
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  <cit> , 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  <dig>  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 .
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: f, where the flux through the coat includes fluxes for both acetate and butyrate re-uptake.
in the wild-type simulations shown in figure  <dig>  given an spf of  <dig> mmolh+hr·gdcw, the flux through the coat was  <dig>  mmolhr·gdcw for acetate re-uptake and  <dig>  mmolhr·gdcw for butyrate re-uptake .
the wild-type flux through bhbd was  <dig>  mmolhr·gdcw to yield a wild-type f: f flux ratio of  <dig> .
the flux ratio is given in equations  <dig>  however, this flux ratio in the wild-type was a function of the spf, and this relationship is shown in figure  <dig>  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  75%,  50%, and  25% of the wild-type value, as shown in figure  <dig>  the second approach was to fix the flux ratio to a specified value for all values of the spf.
the flux ratios chosen were:   <dig>    <dig> ,   <dig> , and   <dig> .
frcoatfrbhbd= <dig>
fbratio results with the f:f flux ratios of figure  <dig> are shown in figure  <dig>  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 − <dig> to  <dig> mmolh+hr·gdcw.
consistent with published experimental results  <cit> , simulated knockdown of the coat resulted in an increased butanol to acetone ratio.
significantly lower acetone production  was predicted with increased butanol production  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:f flux ratio must be represented as a function of the spf .
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  <dig> for flux ratios of:   <dig>    <dig> ,   <dig> , and   <dig> .
similar trends were obtained.
as the f: f flux ratio decreased, metabolic flux was forced through the bhbd enzyme , resulting in  decreased acetate/butyrate re-uptake,  decreased acetone production, and  increased flux through the butanol production pathway.
the exaggerated flux ratio simulations of figure  <dig>  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  <cit> , it was suspected that the aad gene was knocked-down when the asrna for cftb  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  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  <dig>  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  acetyl-coa and  butyryl-coa .
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  <cit> .
to incorporate the knockdown into flux ratio constraints, the f:f flux ratio constraint was increased from a value of  <dig> to  <dig>  to implement the flux ratio constraint at butyryl-coa, a new f:f flux ratio constraint was created.
in wild-type simulations , this flux ratio was unconstrained and had an average value of  <dig>  at an spf value of  <dig> mmolh+hr·gdcw and  <dig> at an spf value of  <dig> mmolh+hr·gdcw .
at highly negative values of the spf, this ratio becomes infinite .
the f:f was constrained to  <dig>   over all values of the spf.
the f:f flux ratio constraint was held constant at  <dig> .
the coat was knocked-down 80% by adjusting the f:f flux ratio to  <dig>  over all values of the spf.
results for solventogenesis  are shown in figure  <dig>  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_ <dig> and aad_ <dig> and is shown in equation  <dig>  this new flux ratio constraint was set to  <dig>  since it is known that acetaldehyde is not produced exclusively through aad_ <dig> .
fbratio results given the f:f 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  as the percent knockdown was increased.
ff= <dig>
over-express aad and knockdown coat
the aad was over-expressed using the f:f and f:f flux ratio constraints.
in this case, the f:f flux ratio constraint was adjusted from  <dig>  to values of   <dig>    <dig> , and   <dig>  to simulate over-expression  by  100%,  200%, and  500% respectively.
the f:f flux ratio constraint was also adjusted accordingly to simulate these levels of over-expression.
fbratio results of solventogenesis  are shown in figure  <dig>  the over-expression of aad under control of its native promoter led to large increases in ethanol production in experimental observations  <cit> .
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% , ethanol production was increased an additional 60% .
this result is possibly explained by the presence of multiple copies of the plasmid present in the cell during the experimental trials.
discussion
the large number of degrees of freedom in primary clostridial metabolism makes this system challenging to model.
initial efforts  <cit>  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 <cit> .
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  <dig>  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 .
from there, constraints can be imposed on product  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 ?
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”  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  acetate,  butyrate/butanol,  ethanol, and  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 icac <dig> 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 .
then, genetic manipulations such as  over-expression,  knockout, and  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  pictures of metabolism.
there are several reasons for this.
first, a fixed glucose uptake rate of  <dig> mmolhr·gdcw was used for all values of the spf examined.
previous results  <cit>  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  <cit> .
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  <cit> .
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 .
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  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  over-expression,  knockout, or  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  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.
CONCLUSIONS
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.
competing interests
the authors declare that they have no competing interests.
authors’ contributions
mm is responsible for constructing the new genome-scale model of c. acetobutylicum, initially developing the fbratio calculations, and writing initial sections of the manuscript.
jy also assisted with developing the fbratio approach and assisted in compiling the genome-scale model.
bf made significant contributions to the background section of the manuscript and assisted in editing.
rs conceived the idea of the fbratio approach and contributed significant writing and editing to the manuscript.
all authors read and approved the final manuscript.
supplementary material
additional file 1
icac490_spreadsheet.
click here for file
additional file 2
icac490_sbml.
click here for file
acknowledgments
this research was funded by the usda afri biobased products and bioenergy production program .
bf was supported by the doctoral scholars program of the institute of critical technologies and applied science  at virginia tech.
jy was supported by a grant from the biodesign and bioprocessing research center at virginia tech.

 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 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: acetate, butyrate, lactate, acetone, ethanol, and several others in small amounts. 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. multiple “genome-scale” metabolic network reconstructions now exist for c. acetobutylicum. similar networks and their corresponding genome-scale models have been reviewed extensively. in general, they are used to complete genome annotation, predict optimal culturing conditions, discover genomic regulation, identify essential genes and drug targets, study strain evolution, and design productive strains. 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 production and high growth rates, and the solventogenic program 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. the total rate at which a cell produces/consumes protons through the several membrane transport mechanisms is termed the specific proton flux, 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 knowledge of the metabolic pathways of butanol fermentation has allowed for targeted engineering approaches. as acids and alcohols 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. an example of this strategy is the knockdown of the butyrate kinase and phosphotransacetylase genes in clostridial metabolism. the proximity of the pta gene in the genome to the acetate kinase 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. the aldehyde/alcohol dehydrogenase gene 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. in a separate metabolic engineering strategy, the aad was over-expressed while knocking-down the gene for subunit b of acetoacetyl-coa transferase. this led to a strain with similar butanol productivity but the ability to produce extraordinary ethanol concentrations of 200 mm. however, when the aad gene was put under control of the ptb gene promoter an increase in butanol concentrations to 300 mm was observed along with faster accumulation of butanol in the culture. metabolic engineering in silico the goal of metabolic engineering in silico is to derive 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, 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 lactate, ethanol, or 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. 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 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. 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 metabolic engineering strategies applied to clostridia, it was not immediately apparent which design 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 “flux balance analysis with flux ratios ” 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. 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. methods genome-scale model a new genome-scale model for c. acetobutylicum atcc was constructed by expanding the previously published model by senger and papoutsakis. the new model is called icac and contains metabolites involved in biochemical reactions, including membrane transport reactions. the model includes genes from the c. acetobutylicum genome. the newly updated icac model differs from the original senger and papoutsakis model in that it contains more reactions and more metabolites. the new reactions added to create the icac model were obtained from the kegg database and recent literature. the icac model also contains an updated tca cycle that operates in both oxidative and reductive directions to succinate, as shown by recent fluxomics studies. the model allows the export of succinate since its metabolic fate has not yet been resolved conclusively. the icac 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 have also been applied. the biomass equation was also updated for the icac model, using the initial version by senger and papoutsakis as a template. a nonlinear optimization procedure was applied 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 icac model is included as additional file the sbml formatted model is included as additional file flux balance analysis the icac 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, 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 maximizing the specific growth rate of the cell while minimizing the total flux of the system. the specific proton flux the concept of the specific proton flux 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 and turns slightly positive during the stationary phase. in this research, the spf was constrained to specific values by constraining the proton exchange reaction. the spf range was between − mmolh+hr·gdcw and mmolh+hr·gdcw, and the limits were chosen from experimental observations. fbratio algorithm metabolic engineering in silico was enabled through the application of flux ratio constraints. in the fbratio method developed in this research, flux ratio constraints were incorporated into the stoichiometric matrix directly, enabling the flux balancing problem to be solved by simple linear programming. several flux ratios were investigated, and a critical metabolic branch point was identified around the use of the acetyl-coa metabolite pool. in clostridial metabolism, acetyl-coa can be used to produce multiple acids, solvents, or macromolecules to produce biomass. the routes taken by acetyl-coa either regenerate or consume different amounts of atp and nad+. ultimately, the balancing of these cofactors determines the production of acids and solvents. the large number of degrees of freedom associated with the genome-scale model of clostridial metabolism allows a large number of acid/solvent production combinations that satisfy cofactor balancing and the overall mass balance, while satisfying the objective functions of the optimization. thus, flux ratios ultimately reduce the number of degrees of freedom of the system and can be used to define selectivity. the following example demonstrates the application of a flux ratio related to the consumption of acetyl-coa and the incorporation of this flux ratio into the stoichiometric matrix. as shown in figure two possibilities for acetyl-coa are usage by the thiolase enzyme for conversion to acetoacetyl-coa and usage by the phosphotransacetylase enzyme for conversion to acetyl phosphate. the reactions catalyzed by these enzymes are given equations.  and 2acetyl−coa→thlacetoacetyl−coa+coa acetyl−coa+orthophosphate→ptaacetyl-phophate+coa next, a ratio of fluxes for these reactions is assumed. for this hypothetical example, it is assumed that twice as much flux proceeds through the thl reaction than the pta reaction. ultimately, biochemical origins of differing fluxes through competing reactions, where known, can be used to calculate flux ratios. the flux through the reaction catalyzed by thl is f and the flux through the pta catalyzed reaction is f. this flux ratio is called f:f and is represented as equation  frthl:frpta=frthlfrpta= to build this flux ratio constraint into the stoichiometric matrix, first equation  is rearranged to the following. frthl−2frpta= next, a new row is added to the stoichiometric matrix. in this new row, two values are added. in the column representing the reaction catalyzed by thl, the coefficient is added to the matrix. in the column representing the reaction catalyzed by the pta, the coefficient − is added. with these additions, when the flux balance equation is solved, the ratio of fluxes for the reactions of equation  and equation  will be exactly 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. in particular, simulations were performed with the icac 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: at highly negative values of spf, acetate and butyrate are produced in high quantities, high hydrogen production accompanies acidogenesis, solvents are produced at spf values close to zero and slightly positive, hydrogen production decreases during solventogenesis, the maximum growth rate of the culture occurs during acidogenesis, the production of butyrate is slightly greater than the production of acetate and much greater than the production of lactate, and 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. knockdown of the coat 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. 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 knockdown coat activity only, knockdown activity of both coat and aad, and knockdown coat while over-expressing aad at and above wild-type levels.  simulations with a minimal set of constraints the icac model was simulated with a glucose uptake rate constrained to mmolhr·gdcw and the spf was varied between − and mmolh+hr·gdcw. 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 − mmolh+hr·gdcw, which is consistent with previous findings. however, during acidogenesis, acetate was the primary acid produced, and acetone was the primary solvent produced during solventogenesis. hydrogen 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+/nadh 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 and nad+ regenerating pathways. approximating wild-type metabolism with fbratio obvious futile cycles allowing artificial nad+/nadh balancing were located and corrected, but the problems described above remained. given the ability of the metabolic network to artificially balance nad+/nadh 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 simulation of metabolism, given these and glucose uptake and spf constraints are shown in figure in this simulation, the following results do not coincide with experimental observations: succinate was produced in high levels during acidogenesis, the production of acetate and lactate far exceed butyrate production, only butyrate was re-consumed, the amount of butyrate re-consumed exceeded the amount of butyrate produced, butanol was produced well before acetone, and hydrogen production fell to zero and rose during solventogenesis. these results show that during early production of butanol, butyrate was consumed through the production pathway. this is consistent with previous findings, and the thermodynamic reversibility calculations insist the reactions catalyzed by the butyrate kinase and the phosphotransbutyrylase remain reversible. *the reduced ferredoxin and the oxidized ferredoxin. a total of five flux ratio constraints were found necessary to generate simulations of metabolism that were qualitatively consistent with experimental observations. first, the butyrate to acetate uptake ratio was constrained according to previously published findings. the assumption was made that the extracellular butyrate concentration was twice that of the extracellular acetate concentration. from the original relationship, this meant that the flux of butyrate uptake through the coat, f, to acetate uptake, f, was equal to, as shown in equation  this flux ratio is referred to as f:f. frcoat,butyrate:frcoat,acetate=fcoat,butyratefcoat,acetate= implementing this flux ratio alone led to high initial production of succinate, and hydrogen production increased during solventogenesis. over-production of ethanol and lactate was also observed. so, a flux ratio was installed to direct the conversion of pyruvate to either lactate, through the lactate dehydrogenase, or to acetyl-coa through the pyruvate ferredoxin oxidoreductase. this flux ratio is shown as equation  this ratio was set equal to to coincide with published experimental observations that lactate production is much less than acetate and butyrate production. frpfofrldh= simulations with these two flux ratios were characterized by high acetate and ethanol production, low butyrate and butanol production, high initial secretion of succinate, and low initial production of hydrogen. the secretion of succinate was investigated next. this is indicative of high fluxes through the tca cycle. no significant succinate export has been reported for c. acetobutylicum, so an optimized metabolic model must produce only minimal succinate. in the icac model, production of oxaloacetate from pyruvate requires hco- when co is produced, it is either transported out of the cell or converted to hco- the fate of co is determined by physicochemical properties of the intracellular environment and has a significant impact on intracellular metabolism. thus, a ratio constraint for co export against conversion was derived. this ratio is shown in equation  and was set equal to to approximate intracellular conditions. this value was chosen because it led to effective simulations. the physicochemical nature of this flux ratio constraint is currently under investigation. fco2exportfco2conversion= this flux ratio constraint corrected hydrogen production 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 high acetate and ethanol production and low butyrate and butanol production. to address this, flux ratios were constructed around the use of the acetyl-coa metabolite pool in clostridial metabolism. acetyl-coa can be utilized to produce atp through acetate production, regenerate nad+ through ethanol production, or 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 with greater flux than through the phosphotransacetylase towards acetate or through the bifunctional aldehyde/alcohol dehydrogenase. to ensure this, additional flux ratio constraints were derived. the ratio of metabolic flux through the thl relative to the pta, f:f, was set equal to and the flux through the thl relative to the aad, f:f, was set equal to 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  and simulation results are shown in figure due to the constant glucose uptake rate, solutions to the flux balance equation were only possible for spf values greater than − mmolh+hr·gdcw. in a previous study, much higher glucose uptake rates enabled flux solutions at spf values approaching − mmolh+hr·gdcw. in the simulations shown in figure a roughly 1: production ratio of h2/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. 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  were applied through fbratio, the icac 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. frthl:frpta=frthlfrpta= frthl:fraad_1=frthlfraad_1= 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, 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 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. 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: f, where the flux through the coat includes fluxes for both acetate and butyrate re-uptake. in the wild-type simulations shown in figure given an spf of mmolh+hr·gdcw, the flux through the coat was mmolhr·gdcw for acetate re-uptake and mmolhr·gdcw for butyrate re-uptake. the wild-type flux through bhbd was mmolhr·gdcw to yield a wild-type f: f flux ratio of. the flux ratio is given in equations  however, this flux ratio in the wild-type was a function of the spf, and this relationship is shown in figure 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 75%, 50%, and 25% of the wild-type value, as shown in figure the second approach was to fix the flux ratio to a specified value for all values of the spf. the flux ratios chosen were:,, and. frcoatfrbhbd= fbratio results with the f:f flux ratios of figure are shown in figure 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 − to mmolh+hr·gdcw. consistent with published experimental results, simulated knockdown of the coat resulted in an increased butanol to acetone ratio. significantly lower acetone production was predicted with increased butanol production 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:f flux ratio must be represented as a function of the spf. 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 for flux ratios of:,, and. similar trends were obtained. as the f: f flux ratio decreased, metabolic flux was forced through the bhbd enzyme, resulting in decreased acetate/butyrate re-uptake, decreased acetone production, and increased flux through the butanol production pathway. the exaggerated flux ratio simulations of figure 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, it was suspected that the aad gene was knocked-down when the asrna for cftb 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 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 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 acetyl-coa and butyryl-coa. 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. to incorporate the knockdown into flux ratio constraints, the f:f flux ratio constraint was increased from a value of to to implement the flux ratio constraint at butyryl-coa, a new f:f flux ratio constraint was created. in wild-type simulations, this flux ratio was unconstrained and had an average value of at an spf value of mmolh+hr·gdcw and at an spf value of mmolh+hr·gdcw. at highly negative values of the spf, this ratio becomes infinite. the f:f was constrained to over all values of the spf. the f:f flux ratio constraint was held constant at. the coat was knocked-down 80% by adjusting the f:f flux ratio to over all values of the spf. results for solventogenesis are shown in figure 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_ and aad_ and is shown in equation  this new flux ratio constraint was set to since it is known that acetaldehyde is not produced exclusively through aad_. fbratio results given the f:f 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 as the percent knockdown was increased. ff= over-express aad and knockdown coat the aad was over-expressed using the f:f and f:f flux ratio constraints. in this case, the f:f flux ratio constraint was adjusted from to values of, and to simulate over-expression by 100%, 200%, and 500% respectively. the f:f flux ratio constraint was also adjusted accordingly to simulate these levels of over-expression. fbratio results of solventogenesis are shown in figure the over-expression of aad under control of its native promoter led to large increases in ethanol production in experimental observations. 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%, ethanol production was increased an additional 60%. this result is possibly explained by the presence of multiple copies of the plasmid present in the cell during the experimental trials. discussion the large number of degrees of freedom in primary clostridial metabolism makes this system challenging to model. initial efforts 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. 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 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. from there, constraints can be imposed on product 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? 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” 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 acetate, butyrate/butanol, ethanol, and 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 icac 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. then, genetic manipulations such as over-expression, knockout, and 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 pictures of metabolism. there are several reasons for this. first, a fixed glucose uptake rate of mmolhr·gdcw 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. 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. 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. 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 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 over-expression, knockout, or 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 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. competing interests the authors declare that they have no competing interests. authors’ contributions mm is responsible for constructing the new genome-scale model of c. acetobutylicum, initially developing the fbratio calculations, and writing initial sections of the manuscript. jy also assisted with developing the fbratio approach and assisted in compiling the genome-scale model. bf made significant contributions to the background section of the manuscript and assisted in editing. rs conceived the idea of the fbratio approach and contributed significant writing and editing to the manuscript. all authors read and approved the final manuscript. supplementary material additional file 1 icac490_spreadsheet. click here for file additional file 2 icac490_sbml. click here for file acknowledgments this research was funded by the usda afri biobased products and bioenergy production program. bf was supported by the doctoral scholars program of the institute of critical technologies and applied science at virginia tech. jy was supported by a grant from the biodesign and bioprocessing research center at virginia tech.
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. five flux ratio constraints were required to achieve a qualitative picture of wild-type metabolism for c. acetobutylicum for the production of: acetate, lactate, butyrate, acetone, butanol, ethanol, co and h 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. a new method called “flux balance analysis with flux ratios ” was developed in this research and applied to a new genome-scale model of clostridium acetobutylicum atcc that contains metabolites and reactions. the method was effective for modeling wild-type and engineered strains of c. acetobutylicum. fbratio is a promising new method for constraining genome-scale models using internal flux ratios. the fbratio approach allowed solutions to be found through standard linear programming. however, few modeling approaches offer predictive capabilities to evaluate potential metabolic engineering strategies in silico.
 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 “flux balance analysis with flux ratios ” was developed in this research and applied to a new genome-scale model of clostridium acetobutylicum atcc that contains metabolites and 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: acetate, lactate, butyrate, acetone, butanol, ethanol, co and h 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. genome-scale modelclostridiaflux ratioflux balance analysismetabolic engineeringsystems biology 