Deriving metabolic engineering strategies from genome-scale modeling with flux ratio constraints

Optimized production of bio-based fuels and chemicals from microbial cell factories is a central goal of systems metabolic engineering. To achieve this goal, a new computational method of using flux balance analysis with flux ratios (FBrAtio) was further developed in this research and applied to five case studies to evaluate and design metabolic engineering strategies. The approach was implemented using publicly available genome-scale metabolic flux models. Synthetic pathways were added to these models along with flux ratio constraints by FBrAtio to achieve increased (i) cellulose production from Arabidopsis thaliana; (ii) isobutanol production from Saccharomyces cerevisiae; (iii) acetone production from Synechocystis sp. PCC6803; (iv) H₂ production from Escherichia coli MG1655; and (v) isopropanol, butanol, and ethanol (IBE) production from engineered Clostridium acetobutylicum. The FBrAtio approach was applied to each case to simulate a metabolic engineering strategy already implemented experimentally, and flux ratios were continually adjusted to find (i) the end-limit of increased production using the existing strategy, (ii) new potential strategies to increase production, and (iii) the impact of these metabolic engineering strategies on product yield and culture growth. The FBrAtio approach has the potential to design “fine-tuned” metabolic engineering strategies in silico that can be implemented directly with available genomic tools.     

Characterization of recombinant strains of the Clostridium acetobutylicum butyrate kinase inactivation mutant: need for new phenomenological models for solventogenesis and butanol inhibition?

Two metabolic engineering tools, namely gene inactivation and gene overexpression, were employed to examine the effects of two genetic modifications on the fermentation characteristics of Clostridium acetobutylicum. Inactivation of the butyrate kinase gene (buk) was examined using strain PJC4BK, while the combined effect of buk inactivation and overexpression of the aad gene-encoding the alcohol aldehyde dehydrogense (AAD) used in butanol formation-was examined using strain PJC4BK(pTAAD). The two strains were characterized in controlled pH > or = 5.0 fermentations, and by a recently enhanced method of metabolic flux analysis. Strain PJC4BK was previously genetically characterized, and fermentation experiments at pH > or = 5.5 demonstrated good, but not exceptional, solvent-production capabilities. Here, we show that this strain is a solvent superproducer in pH > or = 5.0 fermentations producing 225 mM (16.7 g/L) of butanol, 76 mM of acetone (4.4 g/L), and 57 mM (2.6 g/L) of ethanol. Strain PJC4BK(pTAAD) produced similar amounts of butanol and acetone but 98 mM (4.5 g/L) of ethanol. Both strains overcame the 180 mM (13 g/L) butanol toxicity limit, without any selection for butanol tolerance. Work with strain PJC4BK(pTAAD) is the first reported use of dual antibiotic selection in C. acetobutylicum. One antibiotic was used for selection of strain PJC4BK while the second antibiotic selected for the pTAAD presence. Overexpression of aad from pTAAD resulted in increased ethanol production but did not increase butanol titers, thus indicating that AAD did not limit butanol production under these fermentation conditions. Metabolic flux analysis showed a decrease in butyrate formation fluxes by up to 75% and an increase in acetate formation fluxes of up to 100% during early growth. The mean specific butanol and ethanol formation fluxes increased significantly in these recombinant strains, up to 300% and 400%, respectively. Onset of solvent production occurred during the exponential-growth phase when the culture optical density was very low and when total and undissociated butyric acid levels were <1 mM. Butyrate levels were low throughout all fermentations, never exceeding 20 mM. Thus, threshold butyrate concentrations are not necessary for solvent production in these stains, suggesting the need for a new phenomenological model to explain solvent formation.     

Aldehyde-alcohol dehydrogenase and/or thiolase overexpression coupled with CoA transferase downregulation lead to higher alcohol titers and selectivity in Clostridium acetobutylicum fermentations

Metabolic engineering (ME) of Clostridium acetobutylicum has led to increased solvent (butanol, acetone, and ethanol) production and solvent tolerance, thus demonstrating that further efforts have the potential to create strains of industrial importance. With recently developed ME tools, it is now possible to combine genetic modifications and thus implement more advanced ME strategies. We have previously shown that antisense RNA (asRNA)-based downregulation of CoA transferase (CoAT, the first enzyme in the acetone-formation pathway) results in increased butanol to acetone selectivity, but overall reduced butanol yields and titers. In this study the alcohol/aldehyde dehydrogenase (aad) gene (encoding the bifunctional protein AAD responsible for butanol and ethanol production from butyryl-CoA and acetyl-CoA, respectively) was expressed from the phosphotransbutyrylase (ptb) promoter to enhance butanol formation and selectivity, while CoAT downregulation was used to minimize acetone production. This led to early production of high alcohol (butanol plus ethanol) titers, overall solvent titers of 30 g/L, and a higher alcohol/acetone ratio. Metabolic flux analysis revealed the likely depletion of butyryl-CoA. In order to increase then the flux towards butyryl-CoA, we examined the impact of thiolase (THL, thl) overexpression. THL converts acetyl-CoA to acetoacetyl-CoA, the first step of the pathway from acetyl-CoA to butyryl-CoA, and thus, combining thl overexpression with aad overexpression decreased, as expected, acetate and ethanol production while increasing acetone and butyrate formation. thl overexpression in strains with asRNA CoAT downregulation did not significantly alter product formation thus suggesting that a more complex metabolic engineering strategy is necessary to enhance the intracellular butyryl-CoA pool and reduce the acetyl-CoA pool in order to achieve improved butanol titers and selectivity.     

Metabolic engineering of Clostridium acetobutylicum ATCC 824 for isopropanol-butanol-ethanol fermentation

Clostridium acetobutylicum naturally produces acetone as well as butanol and ethanol. Since acetone cannot be used as a biofuel, its production needs to be minimized or suppressed by cell or bioreactor engineering. Thus, there have been attempts to disrupt or inactivate the acetone formation pathway. Here we present another approach, namely, converting acetone to isopropanol by metabolic engineering. Since isopropanol can be used as a fuel additive, the mixture of isopropanol, butanol, and ethanol (IBE) produced by engineered C. acetobutylicum can be directly used as a biofuel. IBE production is achieved by the expression of a primary/secondary alcohol dehydrogenase gene from Clostridium beijerinckii NRRL B-593 (i.e., adh(B-593)) in C. acetobutylicum ATCC 824. To increase the total alcohol titer, a synthetic acetone operon (act operon; adc-ctfA-ctfB) was constructed and expressed to increase the flux toward isopropanol formation. When this engineering strategy was applied to the PJC4BK strain lacking in the buk gene (encoding butyrate kinase), a significantly higher titer and yield of IBE could be achieved. The resulting PJC4BK(pIPA3-Cm2) strain produced 20.4 g/liter of total alcohol. Fermentation could be prolonged by in situ removal of solvents by gas stripping, and 35.6 g/liter of the IBE mixture could be produced in 45 h.     

Switching Clostridium acetobutylicum to an ethanol producer by disruption of the butyrate/butanol fermentative pathway

Solventogenic clostridia are well-known since almost a century due to their unique capability to biosynthesize the solvents acetone and butanol. Based on recently developed genetic engineering tools, a targeted 3-hydroxybutyryl-CoA dehydrogenase (Hbd)-negative mutant of Clostridium acetobutylicum was generated. Interestingly, the entire butyrate/butanol (C₄) metabolic pathway of C. acetobutylicum could be inactivated without a severe growth limitation and indicated the general feasibility to manipulate the central fermentative metabolism for product pattern alteration. Cell extracts of the mutant C. acetobutylicum hbd::int(69) revealed clearly reduced thiolase, Hbd and crotonase but increased NADH-dependent alcohol dehydrogenase enzyme activities as compared to the wildtype strain. Neither butyrate nor butanol were detected in cultures of C. acetobutylicum hbd::int(69), and the formation of molecular hydrogen was significantly reduced. Instead up to 16 and 20 g/l ethanol were produced in glucose and xylose batch cultures, respectively. Further sugar addition in glucose fed-batch fermentations increased the ethanol production to a final titer of 33 g/l, resulting in an ethanol to glucose yield of 0.38 g/g.     

Extracellular glycosyl hydrolase activity of the clostridia producing acetone, butanol, and ethanol

Production of acetone, butanol, ethanol, acetic acid, and butyric acid by three strains of anaerobic bacteria, which we identified as Clostridium acetobutylicum, was studied. The yield of acetone and alcohols in 6% flour medium amounted to 12.7-15 g/l with butanol constituting 51.0-55.6%. Activities of these strains towards xylan, beta-glucan, carboxymethylcellulose, and crystalline and amorphous celluloses were studied. C. acertobutylicum 6, C. acetoburylicum 7, and C. acertobutylicum VKPM B-4786 produced larger amounts of acetone and alcohols and displayed higher cellulase and hemicellulase activities than the type strain C. acetobutylicum ATCC 824. It was demonstrated that starch in the medium could be partially substituted with plant biomass.     

丙酮丁醇发酵菌的分子遗传改造

丙酮丁醇梭菌及拜氏梭菌是重要的ABE（丙酮、丁醇和乙醇）工业生产菌株,其发酵产物中的丙酮和丁醇均为重要的化工原料,汽车发动机试验证明丁醇还是一种性能优于乙醇的极具潜力的生物燃料和燃料添加剂。随着新生物技术的不断发展及工业生产的需求,遗传工程改造不断应用于丙酮丁醇生产菌株。在前人研究及工业实践的基础上,对丙酮丁醇生产菌株的遗传特性及其分子遗传改造取得的进展进行了详细概述。     

限制葡萄糖、葡萄糖/乙酸双底物条件下自由控制丙丁梭菌ABE发酵丙酮浓度和丙酮/丁醇比

提出一种可以提高和自由控制丙丁梭菌ABE发酵丙酮浓度与丙酮/丁醇比的方法。（1）通过控制糖化酶用量、反应时间和温度调节玉米培养基初始葡萄糖浓度,使发酵进入到产溶剂期后,残留葡萄糖浓度降至接近于0 g/L的水平;（2）在葡萄糖受限的条件下,诱导丙丁梭菌合成分泌糖化酶,分解寡糖,将葡萄糖维持于低浓度,进而限制梭菌胞内糖酵解途径的碳代谢和NADH生成速度。与此同时,外添乙酸形成葡萄糖/乙酸双底物环境。在能量代谢基本不受破坏、丁醇未达到抑制浓度的条件下,适度抑制丁醇生产,有效地利用外添乙酸强化丙酮合成;（3）在外添乙酸的基础上,添加适量酿酒酵母,形成丙丁梭菌/酿酒酵母混合发酵体系,提高梭菌对高丁醇浓度的耐受能力。整个发酵体系可以将丙酮浓度和丙酮/丁醇比自由控制在5~12 g/L和0.5~1.0的水平,最大丙酮浓度和丙酮/丁醇比达到11.74 g/L和1.02,并可维持丁醇浓度于10~14 g/L的正常水平,充分满足工业ABE发酵对于丙酮和丁醇产品的不同需求。     

Genome-scale model for Clostridium acetobutylicum: Part II. Development of specific proton flux states and numerically determined sub-systems

A regulated genome-scale model for Clostridium acetobutylicum ATCC 824 was developed based on its metabolic network reconstruction. To aid model convergence and limit the number of flux-vector possible solutions (the size of the phenotypic solution space), modeling strategies were developed to impose a new type of constraint at the endo-exo-metabolome interface. This constraint is termed the specific proton flux state, and its use enabled accurate prediction of the extracellular medium pH during vegetative growth of batch cultures. The specific proton flux refers to the influx or efflux of free protons (per unit biomass) across the cell membrane. A specific proton flux state encompasses a defined range of specific proton fluxes and includes all metabolic flux distributions resulting in a specific proton flux within this range. Effective simulation of time-course batch fermentation required the use of independent flux balance solutions from an optimum set of specific proton flux states. Using a real-coded genetic algorithm to optimize temporal bounds of specific proton flux states, we show that six separate specific proton flux states are required to model vegetative-growth metabolism and accurately predict the extracellular medium pH. Further, we define the apparent proton flux stoichiometry per weak acids efflux and show that this value decreases from approximately 3.5 mol of protons secreted per mole of weak acids at the start of the culture to approximately 0 at the end of vegetative growth. Calculations revealed that when specific weak acids production is maximized in vegetative growth, the net proton exchange between the cell and environment occurs primarily through weak acids efflux (apparent proton flux stoichiometry is 1). However, proton efflux through cation channels during the early stages of acidogenesis was found to be significant. We have also developed the concept of numerically determined sub-systems of genome-scale metabolic networks here as a sub-network with a one-dimensional null space basis set. A numerically determined sub-system was constructed in the genome-scale metabolic network to study the flux magnitudes and directions of acetylornithine transaminase, alanine racemase, and D-alanine transaminase. These results were then used to establish additional constraints for the genome-scale model.     

Metabolic engineering of Clostridium acetobutylicum M5 for highly selective butanol production

To improve butanol selectivity, Clostridium acetobutylicum M5(pIMP1E1AB) was constructed by adhE1-ctfAB complementation of C. acetobutylicum M5, a derivative strain of C. acetobutylicum ATCC 824, which does not produce solvents due to the lack of megaplasmid pSOL1. The gene products of adhE1-ctfAB catalyze the formation of acetoacetate and ethanol/butanol with acid re-assimilation in solventogenesis. Effects of the adhE1-ctfAB complementation of M5 were studied by batch fermentations under various pH and glucose concentrations, and by flux balance analysis using a genome-scale metabolic model for this organism. The metabolically engineered M5(pIMP1E1AB) strain was able to produce 154 mM butanol with 9.9 mM acetone at pH 5.5, resulting in a butanol selectivity (a molar ratio of butanol to total solvents) of 0.84, which is much higher than that (0.57 at pH 5.0 or 0.61 at pH 5.5) of the wild-type strain ATCC 824. Unlike for C. acetobutylicum ATCC 824, a higher level of acetate accumulation was observed during fermentation of the M5 strain complemented with adhE1 and/or ctfAB. A plausible reason for this phenomenon is that the cellular metabolism was shifted towards acetate production to compensate reduced ATP production during the largely growth-associated butanol formation by the M5(pIMP1E1AB) strain.     

SpoIIE regulates sporulation but does not directly affect solventogenesis in Clostridium acetobutylicum ATCC 824

Using gene expression reporter vectors, we examined the activity of the spoIIE promoter in wild-type and spo0A-deleted strains of Clostridium acetobutylicum ATCC 824. In wild-type cells, the spoIIE promoter is active in a transient manner during late solventogenesis, but in strain SKO1, where the sporulation initiator spo0A is disrupted, no spoIIE promoter activity is detectable at any stage of growth. Strains 824(pMSpo) and 824(pASspo) were created to overexpress spoIIE and to decrease spoIIE expression via antisense RNA targeted against spoIIE, respectively. Some cultures of strains 824(pMSpo) degenerated during fermentations by losing the pSOL1 megaplasmid and hence did not produce the solvents ethanol, acetone, and butanol. The frequent degeneration event was shown to require an intact copy of spoIIE. Nondegenerate cultures of 824(pMSpo) exhibited normal growth and solvent production. Strain 824(pASspo) exhibited prolonged solventogenesis characterized by increased production of ethanol (225%), acetone (43%), and butanol (110%). Sporulation in strains harboring pASspo was significantly delayed, with sporulating cells exhibiting altered morphology. These results suggest that SpoIIE has no direct effect on the control of solventogenesis and that the changes in solvent production in spoIIE-downregulated cells are mediated by effects on the cell during sporulation.     

Metabolome remodeling during the acidogenic-solventogenic transition in Clostridium acetobutylicum

The fermentation carried out by the biofuel producer Clostridium acetobutylicum is characterized by two distinct phases. Acidogenesis occurs during exponential growth and involves the rapid production of acids (acetate and butyrate). Solventogenesis initiates as cell growth slows down and involves the production of solvents (butanol, acetone, and ethanol). Using metabolomics, isotope tracers, and quantitative flux modeling, we have mapped the metabolic changes associated with the acidogenic-solventogenic transition. We observed a remarkably ordered series of metabolite concentration changes, involving almost all of the 114 measured metabolites, as the fermentation progresses from acidogenesis to solventogenesis. The intracellular levels of highly abundant amino acids and upper glycolytic intermediates decrease sharply during this transition. NAD(P)H and nucleotide triphosphates levels also decrease during solventogenesis, while low-energy nucleotides accumulate. These changes in metabolite concentrations are accompanied by large changes in intracellular metabolic fluxes. During solventogenesis, carbon flux into amino acids, as well as flux from pyruvate (the last metabolite in glycolysis) into oxaloacetate, decreases by more than 10-fold. This redirects carbon into acetyl coenzyme A, which cascades into solventogenesis. In addition, the electron-consuming reductive tricarboxylic acid (TCA) cycle is shutdown, while the electron-producing oxidative (clockwise) right side of the TCA cycle remains active. Thus, the solventogenic transition involves global remodeling of metabolism to redirect resources (carbon and reducing power) from biomass production into solvent production.     

Expression of plasmid-encoded aad in Clostridium acetobutylicum M5 restores vigorous butanol production.

Mutant M5 of Clostridium acetobutylicum ATCC 824, which produces neither butanol nor acetone and is deficient in butyraldehyde dehydrogenase (BYDH), acetoacetate decarboxylase, and acetoacetyl-coenzyme A:acetate/butyrate:coenzyme A-transferase activities, was transformed with plasmid pCAAD, which carries the gene aad (R. V. Nair, G. N. Bennett, and E. T. Papoutsakis, J. Bacteriol, 176:871-885, 1994). In batch fermentation studies, aad expression restored butanol formation (84 mM) in mutant M5 without any acetone formation or any significant increase in ethanol production. The corresponding protein (AAD) appeared as a ca. 96-kDa band in a denaturing protein gel. Expression of AAD in M5 resulted in restoration of BYDH activity and small increases in the activities of acetaldehyde dehydrogenase, butanol dehydrogenase, and ethanol dehydrogenase. These findings suggest that BYDH activity in C. acetobutylicum ATCC 824 resides largely in AAD, and that AAD's primary role is in the formation of butanol rather than of ethanol.     

Use of immobilized microbial membrane fragments to remove oxygen and favor the acetone-butanol fermentation.

Oxygen-reducing membrane fragments obtained from Escherichia coli were used with Clostridium acetobutylicum (C. acetobutylicum) to provide an oxygen-free microenvironment for the conversion of glucose to acetone, butanol, and ethanol (ABE). The batch fermentation of suspended C. acetobutylicum NRRL-B-643 and its ability to produce solvents in the presence of membranes as the oxygen-elimination agent are described and compared with the conventional sparging technique used to maintain anaerobiosis. The use of membrane fragments to remove oxygen for fermentation by C. acetobutylicum was successful and gave slightly improved results over the use of sparing with regard to lag, biomass, and solvent production (e.g., final butanol concentration of 3.25 and 2.7 g/L, respectively). Solvent production is also reported for a continuous columnar reactor with coimmobilized cells and membranes in kappa-carrageenan gel beads and air-saturated liquid feed.     

New insights into the butyric acid metabolism of Clostridium acetobutylicum

Biosynthesis of acetone and n-butanol is naturally restricted to the group of solventogenic clostridia with Clostridium acetobutylicum being the model organism for acetone-butanol-ethanol (ABE) fermentation. According to limited genetic tools, only a few rational metabolic engineering approaches were conducted in the past to improve the production of butanol, an advanced biofuel. In this study, a phosphotransbutyrylase-(Ptb) negative mutant, C. acetobutylicum ptb::int(87), was generated using the ClosTron methodology for targeted gene knock-out and resulted in a distinct butyrate-negative phenotype. The major end products of fermentation experiments without pH control were acetate (3.2?g/l), lactate (4.0?g/l), and butanol (3.4?g/l). The product pattern of the ptb mutant was altered to high ethanol (12.1?g/l) and butanol (8.0?g/l) titers in pH?≥?5.0-regulated fermentations. Glucose fed-batch cultivation elevated the ethanol concentration to 32.4?g/l, yielding a more than fourfold increased alcohol to acetone ratio as compared to the wildtype. Although butyrate was never detected in cultures of C. acetobutylicum ptb::int(87), the mutant was still capable to take up butyrate when externally added during the late exponential growth phase. These findings suggest that alternative pathways of butyrate re-assimilation exist in C. acetobutylicum, supposably mediated by acetoacetyl-CoA:acyl-CoA transferase and acetoacetate decarboxylase, as well as reverse reactions of butyrate kinase and Ptb with respect to previous studies.     

Effects of introducing heterologous pathways on microbial metabolism with respect to metabolic optimality

Although optimality of microbial metabolism under genetic and environmental perturbations is well studied, the effects of introducing heterologous reactions on the overall metabolism are not well understood. This point is important in the field of metabolic engineering because heterologous reactions are more frequently introduced into various microbial hosts. The genome-scale metabolic simulations of Escherichia coli strains engineered to produce 1,4-butanediol, 1,3-propanediol, and amorphadiene suggest that microbial metabolism shows much different responses to the introduced heterologous reactions in a strain-specific manner than typical gene knockouts in terms of the energetic status (e.g., ATP and biomass generation) and chemical production capacity. The 1,4-butanediol and 1,3-propanediol producers showed greater metabolic optimality than the wild-type strains and gene knockout mutants for the energetic status, while the amorphadiene producer was metabolically less optimal. For the optimal chemical production capacity, additional gene knockouts were most effective for the strain producing 1,3-propanediol, but not for the one producing 1,4-butanediol. These observations suggest that strains having heterologous metabolic reactions have metabolic characteristics significantly different from those of the wild-type strain and single gene knockout mutants. Finally, comparison of the theoretically predicted and ¹³C-based flux values pinpoints pathways with non-optimal flux values, which can be considered as engineering targets in systems metabolic engineering strategies. To our knowledge, this study is the first attempt to quantitatively characterize microbial metabolisms with different heterologous reactions. The suggested potential reasons behind each strain’s different metabolic responses to the introduced heterologous reactions should be carefully considered in strain designs.     

Genomic analysis of carbon monoxide utilization and butanol production by Clostridium carboxidivorans strain P7

Increasing demand for the production of renewable fuels has recently generated a particular interest in microbial production of butanol. Anaerobic bacteria, such as Clostridium spp., can naturally convert carbohydrates into a variety of primary products, including alcohols like butanol. The genetics of microorganisms like Clostridium acetobutylicum have been well studied and their solvent-producing metabolic pathways characterized. In contrast, less is known about the genetics of Clostridium spp. capable of converting syngas or its individual components into solvents. In this study, the type of strain of a new solventogenic Clostridium species, C. carboxidivorans, was genetically characterized by genome sequencing. C. carboxidivorans strain P7(T) possessed a complete Wood-Ljungdahl pathway gene cluster, involving CO and CO(2) fixation and conversion to acetyl-CoA. Moreover, with the exception of an acetone production pathway, all the genetic determinants of canonical ABE metabolic pathways for acetate, butyrate, ethanol and butanol production were present in the P7(T) chromosome. The functionality of these pathways was also confirmed by growth of P7(T) on CO and production of CO(2) as well as volatile fatty acids (acetate and butyrate) and solvents (ethanol and butanol). P7(T) was also found to harbour a 19 Kbp plasmid, which did not include essential or butanol production related genes. This study has generated in depth knowledge of the P7(T) genome, which will be helpful in developing metabolic engineering strategies to improve C. carboxidivorans's natural capacity to produce potential biofuels from syngas.     

Semi-automated Curation of Metabolic Models via Flux Balance Analysis: A Case Study with Mycoplasma gallisepticum

Primarily used for metabolic engineering and synthetic biology, genome-scale metabolic modeling shows tremendous potential as a tool for fundamental research and curation of metabolism. Through a novel integration of flux balance analysis and genetic algorithms, a strategy to curate metabolic networks and facilitate identification of metabolic pathways that may not be directly inferable solely from genome annotation was developed. Specifically, metabolites involved in unknown reactions can be determined, and potentially erroneous pathways can be identified. The procedure developed allows for new fundamental insight into metabolism, as well as acting as a semi-automated curation methodology for genome-scale metabolic modeling. To validate the methodology, a genome-scale metabolic model for the bacterium Mycoplasma gallisepticum was created. Several reactions not predicted by the genome annotation were postulated and validated via the literature. The model predicted an average growth rate of 0.358±0.12, closely matching the experimentally determined growth rate of M. gallisepticum of 0.244±0.03. This work presents a powerful algorithm for facilitating the identification and curation of previously known and new metabolic pathways, as well as presenting the first genome-scale reconstruction of M. gallisepticum.