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.     

Genome-scale analysis of Mannheimia succiniciproducens metabolism

Mannheimia succiniciproducens MBEL55E isolated from bovine rumen is a capnophilic gram-negative bacterium that efficiently produces succinic acid, an industrially important four carbon dicarboxylic acid. In order to design a metabolically engineered strain which is capable of producing succinic acid with high yield and productivity, it is essential to optimize the whole metabolism at the systems level. Consequently, in silico modeling and simulation of the genome-scale metabolic network was employed for genome-scale analysis and efficient design of metabolic engineering experiments. The genome-scale metabolic network of M. succiniciproducens consisting of 686 reactions and 519 metabolites was constructed based on reannotation and validation experiments. With the reconstructed model, the network structure and key metabolic characteristics allowing highly efficient production of succinic acid were deciphered; these include strong PEP carboxylation, branched TCA cycle, relative weak pyruvate formation, the lack of glyoxylate shunt, and non-PTS for glucose uptake. Constraints-based flux analyses were then carried out under various environmental and genetic conditions to validate the genome-scale metabolic model and to decipher the altered metabolic characteristics. Predictions based on constraints-based flux analysis were mostly in excellent agreement with the experimental data. In silico knockout studies allowed prediction of new metabolic engineering strategies for the enhanced production of succinic acid. This genome-scale in silico model can serve as a platform for the systematic prediction of physiological responses of M. succiniciproducens to various environmental and genetic perturbations and consequently for designing rational strategies for strain improvement.     

Fermentation characterization and flux analysis of recombinant strains of Clostridium acetobutylicum with an inactivated solR gene

The effect of solR inactivation on the metabolism of Clostridium acetobutylicum was examined using fermentation characterization and metabolic flux analysis. The solR-inactivated strain (SolRH) of this study had a higher rate of glucose utilization and produced higher solvent concentrations (by 25%, 14%, and 81%, respectively, for butanol, acetone, and ethanol) compared to the wild type. Strain SolRH(pTAAD), carrying a plasmid-encoded copy of the bifunctional alcohol/aldehyde dehydrogenase gene (aad) used in butanol production, produced even higher concentrations of solvents (by 21%, 45%, and 62%, respectively, for butanol, acetone, and ethanol) than strain SolRH. Clarithromycin used for strain SolRH maintenance during SolRH(pTAAD) fermentations did not alter product formation; however, tetracycline used for pTAAD maintenance resulted in 90% lower solvent production. Journal of Industrial Microbiology & Biotechnology (2001) 27, 322–328. Received 12 September 2000/ Accepted in revised form 21 July 2001     

Metabolic engineering of Clostridium acetobutylicum for the enhanced production of isopropanol‐butanol‐ethanol fuel mixture

Butanol is considered as a superior biofuel, which is conventionally produced by clostridial acetone‐butanol‐ethanol (ABE) fermentation. Among ABE, only butanol and ethanol can be used as fuel alternatives. Coproduction of acetone thus causes lower yield of fuel alcohols. Thus, this study aimed at developing an improved Clostridium acetobutylicum strain possessing enhanced fuel alcohol production capability. For this, we previously developed a hyper ABE producing BKM19 strain was further engineered to convert acetone into isopropanol. The BKM19 strain was transformed with the plasmid pIPA100 containing the sadh (primary/secondary alcohol dehydrogenase) and hydG (putative electron transfer protein) genes from the Clostridium beijerinckii NRRL B593 cloned under the control of the thiolase promoter. The resulting BKM19 (pIPA100) strain produced 27.9 g/l isopropanol‐butanol‐ethanol (IBE) as a fuel alcohols with negligible amount of acetone (0.4 g/l) from 97.8 g/l glucose in lab‐scale (2 l) batch fermentation. Thus, this metabolically engineered strain was able to produce 99% of total solvent produced as fuel alcohols. The scalability and stability of BKM19 (pIPA100) were evaluated at 200 l pilot‐scale fermentation, which showed that the fuel alcohol yield could be improved to 0.37 g/g as compared to 0.29 g/g obtained at lab‐scale fermentation, while attaining a similar titer. To the best of our knowledge, this is the highest titer of IBE achieved and the first report on the large scale fermentation of C. acetobutylicum for IBE production. © 2013 American Institute of Chemical Engineers Biotechnol. Prog., 29:1083–1088, 2013     

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.     

Genome-scale analysis of Saccharomyces cerevisiae metabolism and ethanol production in fed-batch culture

A dynamic flux balance model based on a genome-scale metabolic network reconstruction is developed for in silico analysis of Saccharomyces cerevisiae metabolism and ethanol production in fed-batch culture. Metabolic engineering strategies previously identified for their enhanced steady-state biomass and/or ethanol yields are evaluated for fed-batch performance in glucose and glucose/xylose media. Dynamic analysis is shown to provide a single quantitative measure of fed-batch ethanol productivity that explicitly handles the possible tradeoff between the biomass and ethanol yields. Productivity optimization conducted to rank achievable fed-batch performance demonstrates that the genetic manipulation strategy and the fed-batch operating policy should be considered simultaneously. A library of candidate gene insertions is assembled and directly screened for their achievable ethanol productivity in fed-batch culture. A number of novel gene insertions with ethanol productivities identical to the best metabolic engineering strategies reported in previous studies are identified, thereby providing additional targets for experimental evaluation. The top performing gene insertions were substrate dependent, with the highest ranked insertions for glucose media yielding suboptimal performance in glucose/xylose media. The analysis results suggest that enhancements in biomass yield are most beneficial for the enhancement of fed-batch ethanol productivity by recombinant xylose utilizing yeast strains. We conclude that steady-state flux balance analysis is not sufficient to predict fed-batch performance and that the media, genetic manipulations, and fed-batch operating policy should be considered simultaneously to achieve optimal metabolite productivity.     

Genome-scale modeling and in silico analysis of ethanologenic bacteria Zymomonas mobilis

Bioethanol has been recognized as a potential alternative energy source. Among various ethanol-producing microbes, Zymomonas mobilis has acquired special attention due to its higher ethanol yield and tolerance. However, cellular metabolism in Z. mobilis remains unclear, hindering its practical application for bioethanol production. To elucidate such physiological characteristics, we reconstructed and validated a genome-scale metabolic network (iZM363) of Z. mobilis ATCC31821 (ZM4) based on its annotated genome and biochemical information. The phenotypic behaviors and metabolic states predicted by our genome-scale model were highly consistent with the experimental observations of Z. mobilis ZM4 strain growing on glucose as well as NMR-measured intracellular fluxes of an engineered strain utilizing glucose, fructose, and xylose. Subsequent comparative analysis with Escherichia coli and Saccharomyces cerevisiae as well as gene essentiality and flux coupling analyses have also confirmed the functional role of pdc and adh genes in the ethanologenic activity of Z. mobilis, thus leading to better understanding of this natural ethanol producer. In future, the current model could be employed to identify potential cell engineering targets, thereby enhancing the productivity of ethanol in Z. mobilis.     

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.     

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

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

Metabolic engineering of Clostridium acetobutylicum ATCC 824 for increased solvent production by enhancement of acetone formation enzyme activities using a synthetic acetone operon

The ability to genetically alter the product-formation capabilities of Clostridium acetobutylicum is necessary for continued progress toward industrial production of the solvents butanol and acetone by fermentation. Batch fermentations at pH 4.5, 5.5, or 6.5 were conducted using C. acetobutylicum ATCC 824 (pFNK6). Plasmid pFNK6 contains a synthetic operon (the "ace operon") in which the three homologous acetone-formation genas (adc, ctfA, and ctfB) are transcribed from the adc promoter. The corresponding enzymes (acetoacetate decarboxylase and CoA-transferase) were best expressed in pH 4.5 fermentations. However, the highest levels of solvents were attained at pH 5.5. Relative to the plasmid-free control strain at pH 5.5, ATCC 824 (pFNK6) produced 95%, 37%, and 90% higher final concentrations of acetone, butanol, and ethanol, respectively; a 50% higher yield (g/g) of solvents on glucose; and a 22-fold lower mass of residual carboxylic acids. At all pH values, the acetone-formation enzymes were expressed earlier with ATCC 824 (pFNK6) than in control fermentations, leading to earlier induction of acetone formation. Furthermore, strain ATCC 824 (pFNK6) produced butanol significantly earlier in the fermentation and produced significant levels of solvents at pH 6.5. Only trace levels of solvents were produced by strain ATCC 824 at pH 6.5. Compared with ATCC 824, a plasmid-control strain containing a vector without the ace operon also produced higher levels of solvents [although lower than those of strain ATCC 824 (pFNK6)] and lower levels of acids. Strains containing plasmid-borne derivatives of the ace operon, in which either the acetoacetate decarboxylase or CoA-transferase alone were expressed at elevated levels, produced acids and solvents at levels similar to those of the plasmid-control strain. (c) 1993 John Wiley & Sons, Inc.     

Efficient butanol production under aerobic conditions by coculture of Clostridium acetobutylicum and Nesterenkonia sp. strain F

Clostridium acetobutylicum is widely used for the microbial production of butanol in a process known as acetone–butanol–ethanol (ABE) fermentation. However, this process suffers from several disadvantages including high oxygen sensitivity of the bacterium which makes the process complicated and necessitate oxygen elimination in the culture medium. Nesterenkonia sp. strain F has attracted interests as the only known non-Clostridia microorganism with inherent capability of butanol production even in the presence of oxygen. This bacterium is not delimited by oxygen sensitivity, a challenge in butanol biosynthesis, but the butanol titer was far below Clostridia. In this study, Nesterenkonia sp. strain F was cocultivated with C. acetobutylicum to form a powerful “coculture” for butanol production thereby eliminating the need for oxygen removal before fermentation. The response surface method was used for obtaining optimal inoculation amount/time and media formulation. The highest yield, 0.31 g/g ABE (13.6 g/L butanol), was obtained by a coculture initiated with 1.5 mg/L Nesterenkonia sp. strain F and inoculated with 15 mg/L C. acetobutylicum after 1.5 hr in a medium containing 67 g/L glucose, 2.2 g/L yeast extract, 4 g/L peptone, and 1.4% (vol/vol) P2 solution. After butanol toxicity assessment, where Nesterenkonia sp. strain F showed no butanol toxicity, the coculture was implemented in a 2 L fermenter with continual aeration leading to 20 g/L ABE.     

Complete activity profile of Clostridium acetobutylicum [FeFe]-hydrogenase and kinetic parameters for endogenous redox partners

In Clostridium acetobutylicum, [FeFe]-hydrogenase is involved in hydrogen production in vivo by transferring electrons from physiological electron donors, ferredoxin and flavodoxin, to protons. In this report, by modifications of the purification procedure, the specific activity of the enzyme has been improved and its complete catalytic profile in hydrogen evolution, hydrogen uptake, proton/deuterium exchange and para-H2/ortho-H2 conversion has been determined. The major ferredoxin expressed in the solvent-producing C. acetobutylicum cells was purified and identified as encoded by ORF CAC0303. Clostridium acetobutylicum recombinant holoflavodoxin CAC0587 was also purified. The kinetic parameters of C. acetobutylicum [FeFe]-hydrogenase for both physiological partners, ferredoxin CAC0303 and flavodoxin CAC0587, are reported for hydrogen uptake and hydrogen evolution activities.     

Parallel labeling experiments validate Clostridium acetobutylicum metabolic network model for 13C metabolic flux analysis

In this work, we provide new insights into the metabolism of Clostridium acetobutylicum ATCC 824 obtained using a systematic approach for quantifying fluxes based on parallel labeling experiments and ¹³C-metabolic flux analysis (¹³C-MFA). Here, cells were grown in parallel cultures with [1-¹³C]glucose and [U-¹³C]glucose as tracers and ¹³C-MFA was used to quantify intracellular metabolic fluxes. Several metabolic network models were compared: an initial model based on current knowledge, and extended network models that included additional reactions that improved the fits of experimental data. While the initial network model did not produce a statistically acceptable fit of ¹³C-labeling data, an extended network model with five additional reactions was able to fit all data with 292 redundant measurements. The model was subsequently trimmed to produce a minimal network model of C. acetobutylicum for ¹³C-MFA, which could still reproduce all of the experimental data. The flux results provided valuable new insights into the metabolism of C. acetobutylicum. First, we found that TCA cycle was effectively incomplete, as there was no measurable flux between α-ketoglutarate and succinyl-CoA, succinate and fumarate, and malate and oxaloacetate. Second, an active pathway was identified from pyruvate to fumarate via aspartate. Third, we found that isoleucine was produced exclusively through the citramalate synthase pathway in C. acetobutylicum and that CAC3174 was likely responsible for citramalate synthase activity. These model predictions were confirmed in several follow-up tracer experiments. The validated metabolic network model established in this study can be used in future investigations for unbiased ¹³C-flux measurements in C. acetobutylicum.     

基因组规模代谢网络模型构建及其应用

微生物制造产业的发展迫切需要进一步提高认识、设计和改造微生物细胞代谢的能力,以推动工业生物技术快速发展。随着微生物全基因组序列等高通量数据的不断积聚和生物信息学策略的持续涌现,使全局性、系统化地解析、设计、调控微生物生理代谢功能成为可能。而基于基因组序列注释和详细生化信息整合的基因组规模代谢网络模型(GSMM)构建为全局理解和理性调控微生物生理代谢功能提供了最佳平台。以下在详述GSMM的应用基础上,描述了如何构建一个高精确度的GSMM,并展望了未来的发展方向。