Elucidating central metabolic redox obstacles hindering ethanol production in Clostridium thermocellum

Engineering the cellulolytic bacterium, Clostridium thermocellum, to co-utilize hemicellulose

Consolidated bioprocessing (CBP) of lignocellulosic biomass holds promise to realize economic production of second-generation biofuels/chemicals, and Clostridium thermocellum is a leading candidate for CBP due to it being one of the fastest degraders of crystalline cellulose and lignocellulosic biomass. However, CBP by C. thermocellum is approached with co-cultures, because C. thermocellum does not utilize hemicellulose. When compared with a single-species fermentation, the co-culture system introduces unnecessary process complexity that may compromise process robustness. In this study, we engineered C. thermocellum to co-utilize hemicellulose without the need for co-culture. By evolving our previously engineered xylose-utilizing strain in xylose, an evolved clonal isolate (KJC19-9) was obtained and showed improved specific growth rate on xylose by ∼3-fold and displayed comparable growth to a minimally engineered strain grown on the bacteria's naturally preferred substrate, cellobiose. To enable full xylan deconstruction to xylose, we recombinantly expressed three different β-xylosidase enzymes originating from Thermoanaerobacterium saccharolyticum into KJC19-9 and demonstrated growth on xylan with one of the enzymes. This recombinant strain was capable of co-utilizing cellulose and xylan simultaneously, and we integrated the β-xylosidase gene into the KJC19-9 genome, creating the KJCBXint strain. The strain, KJC19-9, consumed monomeric xylose but accumulated xylobiose when grown on pretreated corn stover, whereas the final KJCBXint strain showed significantly greater deconstruction of xylan and xylobiose. This is the first reported C. thermocellum strain capable of degrading and assimilating hemicellulose polysaccharide while retaining its cellulolytic capabilities, unlocking significant potential for CBP in advancing the bioeconomy.

Kinetic characterization of annotated glycolytic enzymes present in cellulose-fermenting Clostridium thermocellum suggests different metabolic roles

BACKGROUND

The efficient production of sustainable biofuels is important for the reduction of greenhouse gas emissions. Clostridium thermocellum ATCC 27405 is a candidate for ethanol production from lignocellulosic biomass using consolidated bioprocessing. Fermentation of cellulosic biomass goes through an atypical glycolytic pathway in this thermophilic bacterium, with various glycolytic enzymes capable of utilizing different phosphate donors, including GTP and inorganic pyrophosphate (PP_i), in addition to or in place of the usual ATP. C. thermocellum contains three annotated phosphofructokinases (PFK) genes, the expression of which have all been detected through proteomics and transcriptomics. Pfp (Cthe_0347) was previously characterized as pyrophosphate dependent with fructose-6-phosphate (F6P) as its substrate.

RESULTS

We now demonstrate that this enzyme can also phosphorylate sedoheptulose-7-phosphate (an intermediate in the pentose phosphate pathway), with the V_max and K_m of F6P being approximately 15 folds higher and 43 folds lower, respectively, in comparison to sedoheptulose-7-phosphate. Purified PfkA shows preference for GTP as the phosphate donor as opposed to ATP with a 12.5-fold difference in K_m values while phosphorylating F6P. Allosteric regulation is a factor at play in PfkA activity, with F6P exhibiting positive cooperativity, and an apparent requirement for ammonium ions to attain maximal activity. Phosphoenolpyruvate and PP_i were the only inhibitors for PfkA determined from the study, which corroborates what is known about enzymes from this subfamily. The activation or inhibition by these ligands lends support to the argument that glycolysis is regulated by metabolites such as PP_i and NH₄⁺ in the organism. PfkB, showed no activity with F6P, but had significant activity with fructose, while utilizing either ATP or GTP, making it a fructokinase. Rounding out the upper glycolysis pathway, the identity of the fructose-1,6-bisphosphate aldolase in the genome was verified and reported to have substantial activity with fructose-1,6-bisphosphate, in the presence of the divalent ion, Zn²⁺.

CONCLUSION

These findings along with previous proteomic data suggest that Pfp, plays a role in both glycolysis and the pentose phosphate pathway, while PfkA and PfkB may phosphorylate sugars in glycolysis but is responsible for sugar metabolism elsewhere under conditions outside of growth on sufficient cellobiose.

Rewiring metabolism of Clostridium thermocellum for consolidated bioprocessing of lignocellulosic biomass poplar to produce short-chain esters

Consolidated bioprocessing (CBP) of lignocellulosic biomass uses cellulolytic microorganisms to enable enzyme production, saccharification, and fermentation to produce biofuels, biochemicals, and biomaterials in a single step. However, understanding and redirecting metabolisms of these microorganisms compatible with CBP are limited. Here, a cellulolytic thermophile Clostridium thermocellum was engineered and demonstrated to be compatible with CBP integrated with a Co-solvent Enhanced Lignocellulosic Fractionation (CELF) pretreatment for conversion of hardwood poplar into short-chain esters with industrial use as solvents, flavors, fragrances, and biofuels. The recombinant C. thermocellum engineered with deletion of carbohydrate esterases and stable overexpression of alcohol acetyltransferases improved ester production without compromised deacetylation activities. These esterases were discovered to exhibit promiscuous thioesterase activities and their deletion enhanced ester production by rerouting the electron and carbon metabolism. Ester production was further improved up to 80-fold and ester composition could be modulated by deleting lactate biosynthesis and using poplar with different pretreatment severity.

Toward implementation of carbon-conservation networks in nonmodel organisms

Decarboxylation - the release of carbon dioxide (CO₂) from a substrate - reduces the carbon yield of bioproduced chemicals. When overlaid onto central carbon metabolism, carbon-conservation networks (CCNs) that reroute flux around CO₂ release can theoretically achieve higher carbon yields for products derived from intermediates that traditionally require CO₂ release, such as acetyl-CoA. Recently, CCNs have started to be implemented in model organisms to produce compounds at higher carbon yields. However, implementation of CCNs in nonmodel hosts may have the greatest impact given their ability to assimilate a larger array of feedstocks, greater environmental tolerance, and unique biosynthetic pathways, ultimately enabling access to a wider range of products. Here, we review recent advances in CCNs with a focus on their application to nonmodel organisms. The differences in central carbon metabolism among different nonmodel hosts reveal opportunities to engineer and apply new CCNs.

A detailed genome-scale metabolic model of Clostridium thermocellum investigates sources of pyrophosphate for driving glycolysis

Lignocellulosic biomass is an abundant and renewable source of carbon for chemical manufacturing, yet it is cumbersome in conventional processes. A promising, and increasingly studied, candidate for lignocellulose bioprocessing is the thermophilic anaerobe Clostridium thermocellum given its potential to produce ethanol, organic acids, and hydrogen gas from lignocellulosic biomass under high substrate loading. Possessing an atypical glycolytic pathway which substitutes GTP or pyrophosphate (PP_i) for ATP in some steps, including in the energy-investment phase, identification, and manipulation of PP_i sources are key to engineering its metabolism. Previous efforts to identify the primary pyrophosphate have been unsuccessful. Here, we explore pyrophosphate metabolism through reconstructing, updating, and analyzing a new genome-scale stoichiometric model for C. thermocellum, iCTH669. Hundreds of changes to the former GEM, iCBI655, including correcting cofactor usages, addressing charge and elemental balance, standardizing biomass composition, and incorporating the latest experimental evidence led to a MEMOTE score improvement to 94%. We found agreement of iCTH669 model predictions across all available fermentation and biomass yield datasets. The feasibility of hundreds of PP_i synthesis routes, newly identified and previously proposed, were assessed through the lens of the iCTH669 model including biomass synthesis, tRNA synthesis, newly identified sources, and previously proposed PP_i-generating cycles. In all cases, the metabolic cost of PP_i synthesis is at best equivalent to investment of one ATP suggesting no direct energetic advantage for the cofactor substitution in C. thermocellum. Even though no unique source of PP_i could be gleaned by the model, by combining with gene expression data two most likely scenarios emerge. First, previously investigated PP_i sources likely account for most PP_i production in wild-type strains. Second, alternate metabolic routes as encoded by iCTH669 can collectively maintain PP_i levels even when previously investigated synthesis cycles are disrupted. Model iCTH669 is available at github.com/maranasgroup/iCTH669.

Assessing the impact of substrate-level enzyme regulations limiting ethanol titer in Clostridium thermocellum using a core kinetic model

Clostridium thermocellum is a promising candidate for consolidated bioprocessing because it can directly ferment cellulose to ethanol. Despite significant efforts, achieved yields and titers fall below industrially relevant targets. This implies that there still exist unknown enzymatic, regulatory, and/or possibly thermodynamic bottlenecks that can throttle back metabolic flow. By (i) elucidating internal metabolic fluxes in wild-type C. thermocellum grown on cellobiose via ¹³C-metabolic flux analysis (¹³C-MFA), (ii) parameterizing a core kinetic model, and (iii) subsequently deploying an ensemble-docking workflow for discovering substrate-level regulations, this paper aims to reveal some of these factors and expand our knowledgebase governing C. thermocellum metabolism. Generated ¹³C labeling data were used with ¹³C-MFA to generate a wild-type flux distribution for the metabolic network. Notably, flux elucidation through MFA alluded to serine generation via the mercaptopyruvate pathway. Using the elucidated flux distributions in conjunction with batch fermentation process yield data for various mutant strains, we constructed a kinetic model of C. thermocellum core metabolism (i.e. k-ctherm138). Subsequently, we used the parameterized kinetic model to explore the effect of removing substrate-level regulations on ethanol yield and titer. Upon exploring all possible simultaneous (up to four) regulation removals we identified combinations that lead to many-fold model predicted improvement in ethanol titer. In addition, by coupling a systematic method for identifying putative competitive inhibitory mechanisms using K-FIT kinetic parameterization with the ensemble-docking workflow, we flagged 67 putative substrate-level inhibition mechanisms across central carbon metabolism supported by both kinetic formalism and docking analysis.

Utilization of Monosaccharides by Hungateiclostridium thermocellum ATCC 27405 through Adaptive Evolution

Hungateiclostridium thermocellum ATCC 27405 is a promising bacterium for consolidated bioprocessing with a robust ability to degrade lignocellulosic biomass through a multienzyme cellulosomal complex. The bacterium uses the released cellodextrins, glucose polymers of different lengths, as its primary carbon source and energy. In contrast, the bacterium exhibits poor growth on monosaccharides such as fructose and glucose. This phenomenon raises many important questions concerning its glycolytic pathways and sugar transport systems. Until now, the detailed mechanisms of H. thermocellum adaptation to growth on hexose sugars have been relatively poorly explored. In this study, adaptive laboratory evolution was applied to train the bacterium in hexose sugars-based media, and genome resequencing was used to detect the genes that got mutated during adaptation period. RNA-seq data of the first culture growing on either fructose or glucose revealed that several glycolytic genes in the Embden–Mayerhof–Parnas pathway were expressed at lower levels in these cells than in cellobiose-grown cells. After seven consecutive transfer events on fructose and glucose (~42 generations for fructose-adapted cells and ~40 generations for glucose-adapted cells), several genes in the EMP glycolysis of the evolved strains increased the levels of mRNA expression, accompanied by a faster growth, a greater biomass yield, a higher ethanol titer than those in their parent strains. Genomic screening also revealed several mutation events in the genomes of the evolved strains, especially in those responsible for sugar transport and central carbon metabolism. Consequently, these genes could be applied as potential targets for further metabolic engineering to improve this bacterium for bio-industrial usage.

Engineering promiscuity of chloramphenicol acetyltransferase for microbial designer ester biosynthesis

Robust and efficient enzymes are essential modules for metabolic engineering and synthetic biology strategies across biological systems to engineer whole-cell biocatalysts. By condensing an acyl-CoA and an alcohol, alcohol acyltransferases (AATs) can serve as interchangeable metabolic modules for microbial biosynthesis of a diverse class of ester molecules with broad applications as flavors, fragrances, solvents, and drop-in biofuels. However, the current lack of robust and efficient AATs significantly limits their compatibility with heterologous precursor pathways and microbial hosts. Through bioprospecting and rational protein engineering, we identified and engineered promiscuity of chloramphenicol acetyltransferases (CATs) from mesophilic prokaryotes to function as robust and efficient AATs compatible with at least 21 alcohol and 8 acyl-CoA substrates for microbial biosynthesis of linear, branched, saturated, unsaturated and/or aromatic esters. By plugging the best engineered CAT (CATec3 Y20F) into the gram-negative mesophilic bacterium Escherichia coli, we demonstrated that the recombinant strain could effectively convert various alcohols into desirable esters, for instance, achieving a titer of 13.9 g/L isoamyl acetate with 95% conversion by fed-batch fermentation. The recombinant E. coli was also capable of simulating the ester profile of roses with high conversion (>97%) and titer (>1 g/L) from fermentable sugars at 37 °C. Likewise, a recombinant gram-positive, cellulolytic, thermophilic bacterium Clostridium thermocellum harboring CATec3 Y20F could produce many of these esters from recalcitrant cellulosic biomass at elevated temperatures (>50 °C) due to the engineered enzyme's remarkable thermostability. Overall, the engineered CATs can serve as a robust and efficient platform for designer ester biosynthesis from renewable and sustainable feedstocks.

Modeling coordinated enzymatic control of saccharification and fermentation by Clostridium thermocellum during consolidated bioprocessing of cellulose

Consolidated bioprocessing (CBP) of cellulose is a cost-effective route to produce valuable biochemicals by integrating saccharification, fermentation and cellulase synthesis in a single step. However, the lack of understanding of governing factors of interdependent saccharification and fermentation in CBP eludes reliable process optimization. Here, we propose a new framework that synergistically couples population balances (to simulate cellulose depolymerization) and cybernetic models (to model enzymatic regulation of fermentation) to enable improved understanding of CBP. The resulting framework, named the unified cybernetic-population balance model (UC-PBM), enables simulation of CBP driven by coordinated control of enzyme synthesis through closed-loop interactions. UC-PBM considers two key aspects in controlling CBP: (1) heterogeneity in cellulose properties and (2) cellular regulation of competing cell growth and cellulase secretion. In a case study on Clostridium thermocellum, UC-PBM not only provides a decent fit with various exometabolomic data, but also reveals that: (i) growth-decoupled cellulase-secreting pathways are only activated during famine conditions to promote the production of growth substrates, and (ii) starting cellulose concentration has a strong influence on the overall flux distribution. Equipped with mechanisms of cellulose degradation and fermentative regulations, UC-PBM is practical to explore phenotypic functions for primary evaluation of microorganisms' potential for metabolic engineering and optimal design of bioprocess.

Modeling a co-culture of Clostridium autoethanogenum and Clostridium kluyveri to increase syngas conversion to medium-chain fatty-acids

•

We model a co-culture of C. autoethanogenum and C. kluyveri for syngas fermentation.

•

Biomass species ratio affects ethanol and acetate profiles in the co-culture.

•

The model predicts that addition of succinate increases caproate production.

•

Genetic interventions  in C. autoethanogenum could increase caproate production.

Microbial fermentation of synthesis gas (syngas) is becoming more attractive for sustainable production of commodity chemicals. To date, syngas fermentation focuses mainly on the use of Clostridium species for the production of small organic molecules such as ethanol and acetate. The co-cultivation of syngas-fermenting microorganisms with chain-elongating bacteria can expand the range of possible products, allowing, for instance, the production of medium-chain fatty acids (MCFA) and alcohols from syngas. To explore these possibilities, we report herein a genome-scale, constraint-based metabolic model to describe growth of a co-culture of Clostridium autoethanogenum and Clostridium kluyveri on syngas for the production of valuable compounds. Community flux balance analysis was used to gain insight into the metabolism of the two strains and their interactions, and to reveal potential strategies enabling production of butyrate and hexanoate.

The model suggests that one strategy to optimize the production of medium-chain fatty-acids from syngas would be the addition of succinate. According to the prediction, addition of succinate would increase the pool of crotonyl-CoA and the ethanol/acetate uptake ratio in C. kluyveri, resulting in a flux of up to 60% of electrons into hexanoate. Another potential way to further optimize butyrate and hexanoate production would be an increase of C. autoethanogenum ethanol production. Blocking either acetaldehyde dehydrogenase or formate dehydrogenase (ferredoxin) activity or formate transport, in the C. autoethanogenum metabolic model could potentially lead to an up to 150% increase in ethanol production.

Development of a Genome-Scale Metabolic Model of Clostridium thermocellum and Its Applications for Integration of Multi-Omics Datasets and Computational Strain Design

Solving environmental and social challenges such as climate change requires a shift from our current non-renewable manufacturing model to a sustainable bioeconomy. To lower carbon emissions in the production of fuels and chemicals, plant biomass feedstocks can replace petroleum using microorganisms as biocatalysts. The anaerobic thermophile Clostridium thermocellum is a promising bacterium for bioconversion due to its capability to efficiently degrade lignocellulosic biomass. However, the complex metabolism of C. thermocellum is not fully understood, hindering metabolic engineering to achieve high titers, rates, and yields of targeted molecules. In this study, we developed an updated genome-scale metabolic model of C. thermocellum that accounts for recent metabolic findings, has improved prediction accuracy, and is standard-conformant to ensure easy reproducibility. We illustrated two applications of the developed model. We first formulated a multi-omics integration protocol and used it to understand redox metabolism and potential bottlenecks in biofuel (e.g., ethanol) production in C. thermocellum. Second, we used the metabolic model to design modular cells for efficient production of alcohols and esters with broad applications as flavors, fragrances, solvents, and fuels. The proposed designs not only feature intuitive push-and-pull metabolic engineering strategies, but also present novel manipulations around important central metabolic branch-points. We anticipate the developed genome-scale metabolic model will provide a useful tool for system analysis of C. thermocellum metabolism to fundamentally understand its physiology and guide metabolic engineering strategies to rapidly generate modular production strains for effective biosynthesis of biofuels and biochemicals from lignocellulosic biomass.

Endogenous carbohydrate esterases of Clostridium thermocellum are identified and disrupted for enhanced isobutyl acetate production from cellulose

Medium-chain esters are versatile chemicals with broad applications as flavors, fragrances, solvents, and potential drop-in biofuels. Currently, these esters are largely produced by the conventional chemical process that uses harsh operating conditions and requires high energy input. Alternatively, the microbial conversion route has recently emerged as a promising platform for sustainable and renewable ester production. The ester biosynthesis pathways can utilize either lipases or alcohol acyltransferase (AAT), but the AAT-dependent pathway is more thermodynamically favorable in an aqueous fermentation environment. Even though a cellulolytic thermophile Clostridium thermocellum harboring an AAT-dependent pathway has recently been engineered for direct conversion of lignocellulosic biomass into esters, the production is not efficient. One potential bottleneck is the ester degradation caused by the endogenous carbohydrate esterases (CEs) whose functional roles are poorly understood. The challenge is to identify and disrupt CEs that can alleviate ester degradation while not negatively affecting the efficient and robust capability of C. thermocellum for lignocellulosic biomass deconstruction. In this study, by using bioinformatics, comparative genomics, and enzymatic analysis to screen a library of CEs, we identified and disrupted the two most critical CEs, Clo1313_0613 and Clo1313_0693, that significantly contribute to isobutyl acetate degradation in C. thermocellum. We demonstrated that an engineered esterase-deficient C. thermocellum strain not only reduced ester hydrolysis but also improved isobutyl acetate production while maintaining effective cellulose assimilation.

Thermodynamic analysis of the pathway for ethanol production from cellobiose in Clostridium thermocellum

Clostridium thermocellum is a candidate for consolidated bioprocessing by carrying out both cellulose solubilization and fermentation. However, despite significant efforts the maximum ethanol titer achieved to date remains below industrially required targets. Several studies have analyzed the impact of increasing ethanol concentration on C. thermocellum's membrane properties, cofactor pool ratios, and altered enzyme regulation. In this study, we explore the extent to which thermodynamic equilibrium limits maximum ethanol titer. We used the max-min driving force (MDF) algorithm (Noor et al., 2014) to identify the range of allowable metabolite concentrations that maintain a negative free energy change for all reaction steps in the pathway from cellobiose to ethanol. To this end, we used a time-series metabolite concentration dataset to flag five reactions (phosphofructokinase (PFK), fructose bisphosphate aldolase (FBA), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), aldehyde dehydrogenase (ALDH) and alcohol dehydrogenase (ADH)) which become thermodynamic bottlenecks under high external ethanol concentrations. Thermodynamic analysis was also deployed in a prospective mode to evaluate genetic interventions which can improve pathway thermodynamics by generating minimal set of reactions or elementary flux modes (EFMs) which possess unique genetic variations while ensuring mass and redox balance with ethanol production. MDF evaluation of all generated (336) EFMs indicated that, i) pyruvate phosphate dikinase (PPDK) has a higher pathway MDF than the malate shunt alternative due to limiting CO₂ concentrations under physiological conditions, and ii) NADPH-dependent glyceraldehyde-3-phosphate dehydrogenase (GAPN) can alleviate thermodynamic bottlenecks at high ethanol concentrations due to cofactor modification and reduction in ATP generation. The combination of ATP linked phosphofructokinase (PFK-ATP) and NADPH linked alcohol dehydrogenase (ADH-NADPH) with NADPH linked aldehyde dehydrogenase (ALDH-NADPH) or ferredoxin: NADP ​+ ​oxidoreductase (NADPH-FNOR) emerges as the best intervention strategy for ethanol production that balances MDF improvements with ATP generation, and appears to functionally reproduce the pathway employed by the ethanologen Thermoanaerobacterium saccharolyticum. Expanding the list of measured intracellular metabolites and improving the quantification accuracy of measurements was found to improve the fidelity of pathway thermodynamics analysis in C. thermocellum. This study demonstrates even before addressing an organism's enzyme kinetics and allosteric regulations, pathway thermodynamics can flag pathway bottlenecks and identify testable strategies for enhancing pathway thermodynamic feasibility and function.

Rex in Caldicellulosiruptor bescii: Novel regulon members and its effect on the production of ethanol and overflow metabolites

Rex is a global redox-sensing transcription factor that senses and responds to the intracellular [NADH]/[NAD+ ] ratio to regulate genes for central metabolism, and a variety of metabolic processes in Gram-positive bacteria. We decipher and validate four new members of the Rex regulon in Caldicellulosiruptor bescii; a gene encoding a class V aminotransferase, the HydG FeFe Hydrogenase maturation protein, an oxidoreductase, and a gene encoding a hypothetical protein. Structural genes for the NiFe and FeFe hydrogenases, pyruvate:ferredoxin oxidoreductase, as well as the rex gene itself are also members of this regulon, as has been predicted previously in different organisms. A C. bescii rex deletion strain constructed in an ethanol-producing strain made 54% more ethanol (0.16 mmol/L) than its genetic parent after 36 hr of fermentation, though only under nitrogen limited conditions. Metabolomic interrogation shows this rex-deficient ethanol-producing strain synthesizes other reduced overflow metabolism products likely in response to more reduced intracellular redox conditions and the accumulation of pyruvate. These results suggest ethanol production is strongly dependent on the native intracellular redox state in C. bescii, and highlight the combined promise of using this gene and manipulation of culture conditions to yield strains capable of producing ethanol at higher yields and final titer.

Single mutation at a highly conserved region of chloramphenicol acetyltransferase enables isobutyl acetate production directly from cellulose by Clostridium thermocellum at elevated temperatures

BACKGROUND

Esters are versatile chemicals and potential drop-in biofuels. To develop a sustainable production platform, microbial ester biosynthesis using alcohol acetyltransferases (AATs) has been studied for decades. Volatility of esters endows high-temperature fermentation with advantageous downstream product separation. However, due to the limited thermostability of AATs known, the ester biosynthesis has largely relied on use of mesophilic microbes. Therefore, developing thermostable AATs is important for ester production directly from lignocellulosic biomass by the thermophilic consolidated bioprocessing (CBP) microbes, e.g., Clostridium thermocellum.

RESULTS

In this study, we engineered a thermostable chloramphenicol acetyltransferase from Staphylococcus aureus (CAT_Sa) for enhanced isobutyl acetate production at elevated temperatures. We first analyzed the broad alcohol substrate range of CAT_Sa. Then, we targeted a highly conserved region in the binding pocket of CAT_Sa for mutagenesis. The mutagenesis revealed that F97W significantly increased conversion of isobutanol to isobutyl acetate. Using CAT_Sa F97W, we demonstrated direct conversion of cellulose into isobutyl acetate by an engineered C. thermocellum at elevated temperatures.

CONCLUSIONS

This study highlights that CAT is a potential thermostable AAT that can be harnessed to develop the thermophilic CBP microbial platform for biosynthesis of designer bioesters directly from lignocellulosic biomass.

Redox cofactor engineering in industrial microorganisms: strategies, recent applications and future directions

NAD and NADP, a pivotal class of cofactors, which function as essential electron donors or acceptors in all biological organisms, drive considerable catabolic and anabolic reactions. Furthermore, they play critical roles in maintaining intracellular redox homeostasis. However, many metabolic engineering efforts in industrial microorganisms towards modification or introduction of metabolic pathways, especially those involving consumption, generation or transformation of NAD/NADP, often induce fluctuations in redox state, which dramatically impede cellular metabolism, resulting in decreased growth performance and biosynthetic capacity. Here, we comprehensively review the cofactor engineering strategies for solving the problematic redox imbalance in metabolism modification, as well as their features, suitabilities and recent applications. Some representative examples of in vitro biocatalysis are also described. In addition, we briefly discuss how tools and methods from the field of synthetic biology can be applied for cofactor engineering. Finally, future directions and challenges for development of cofactor redox engineering are presented.

Targeted redox and energy cofactor metabolomics in Clostridium thermocellum and Thermoanaerobacterium saccharolyticum

BACKGROUND

Clostridium thermocellum and Thermoanaerobacterium saccharolyticum are prominent candidate biocatalysts that, together, can enable the direct biotic conversion of lignocellulosic biomass to ethanol. The imbalance and suboptimal turnover rates of redox cofactors are currently hindering engineering efforts to achieve higher bioproductivity in both organisms. Measuring relevant intracellular cofactor concentrations will help understand redox state of these cofactors and help identify a strategy to overcome these limitations; however, metabolomic determinations of these labile metabolites have historically proved challenging.

RESULTS

Through our validations, we verified the handling and storage stability of these metabolites, and verified extraction matrices and extraction solvent were not suppressing mass spectrometry signals. We recovered adenylate energy charge ratios (a main quality indicator) above 0.82 for all extractions. NADH/NAD+ values of 0.26 and 0.04 for an adhE-deficient strain of C. thermocellum and its parent, respectively, reflect the expected shift to a more reduced redox potential when a species lacks the ability to re-oxidize NADH by synthesizing ethanol. This method failed to yield reliable results with C. bescii and poor-growing strains of T. saccharolyticum.

CONCLUSIONS

Our validated protocols demonstrate and validate the extraction and analysis of selected redox and energy-related metabolites from two candidate consolidated bioprocessing biocatalysts, C. thermocellum and T. saccharolyticum. This development and validation highlights the important, but often neglected, need to optimize and validate metabolomic protocols when adapting them to new cell or tissue types.

Overflow metabolism and growth cessation in Clostridium thermocellum DSM1313 during high cellulose loading fermentations

As a model thermophilic bacterium for the production of second-generation biofuels, the metabolism of Clostridium thermocellum has been widely studied. However, most studies have characterized C. thermocellum metabolism for growth at relatively low substrate concentrations. This outlook is not industrially relevant, however, as commercial viability requires substrate loadings of at least 100 g/L cellulosic materials. Recently, a wild-type C. thermocellum DSM1313 was cultured on high cellulose loading batch fermentations and reported to produce a wide range of fermentative products not seen at lower substrate concentrations, opening the door for a more in-depth analysis of how this organism will behave in industrially relevant conditions. In this work, we elucidated the interconnectedness of overflow metabolism and growth cessation in C. thermocellum during high cellulose loading batch fermentations (100 g/L). Metabolic flux and thermodynamic analyses suggested that hydrogen and formate accumulation perturbed the complex redox metabolism and limited conversion of pyruvate to acetyl-CoA conversion, likely leading to overflow metabolism and growth cessation in C. thermocellum. Pyruvate formate lyase (PFL) acts as an important redox valve and its flux is inhibited by formate accumulation. Finally, we demonstrated that manipulation of fermentation conditions to alleviate hydrogen accumulation could dramatically alter the fate of pyruvate, providing valuable insight into process design for enhanced C. thermocellum production of chemicals and biofuels. Biotechnol. Bioeng. 2017;114: 2592-2604. © 2017 Wiley Periodicals, Inc.

Development of a core Clostridium thermocellum kinetic metabolic model consistent with multiple genetic perturbations

BACKGROUND

Clostridium thermocellum is a Gram-positive anaerobe with the ability to hydrolyze and metabolize cellulose into biofuels such as ethanol, making it an attractive candidate for consolidated bioprocessing (CBP). At present, metabolic engineering in C. thermocellum is hindered due to the incomplete description of its metabolic repertoire and regulation within a predictive metabolic model. Genome-scale metabolic (GSM) models augmented with kinetic models of metabolism have been shown to be effective at recapitulating perturbed metabolic phenotypes.

RESULTS

In this effort, we first update a second-generation genome-scale metabolic model (iCth446) for C. thermocellum by correcting cofactor dependencies, restoring elemental and charge balances, and updating GAM and NGAM values to improve phenotype predictions. The iCth446 model is next used as a scaffold to develop a core kinetic model (k-ctherm118) of the C. thermocellum central metabolism using the Ensemble Modeling (EM) paradigm. Model parameterization is carried out by simultaneously imposing fermentation yield data in lactate, malate, acetate, and hydrogen production pathways for 19 measured metabolites spanning a library of 19 distinct single and multiple gene knockout mutants along with 18 intracellular metabolite concentration data for a Δgldh mutant and ten experimentally measured Michaelis-Menten kinetic parameters.

CONCLUSIONS

The k-ctherm118 model captures significant metabolic changes caused by (1) nitrogen limitation leading to increased yields for lactate, pyruvate, and amino acids, and (2) ethanol stress causing an increase in intracellular sugar phosphate concentrations (~1.5-fold) due to upregulation of cofactor pools. Robustness analysis of k-ctherm118 alludes to the presence of a secondary activity of ketol-acid reductoisomerase and possible regulation by valine and/or leucine pool levels. In addition, cross-validation and robustness analysis allude to missing elements in k-ctherm118 and suggest additional experiments to improve kinetic model prediction fidelity. Overall, the study quantitatively assesses the advantages of EM-based kinetic modeling towards improved prediction of C. thermocellum metabolism and develops a predictive kinetic model which can be used to design biofuel-overproducing strains.

Specialized activities and expression differences for Clostridium thermocellum biofilm and planktonic cells

Clostridium (Ruminiclostridium) thermocellum is a model organism for its ability to deconstruct plant biomass and convert the cellulose into ethanol. The bacterium forms biofilms adherent to lignocellulosic feedstocks in a continuous cell-monolayer in order to efficiently break down and uptake cellulose hydrolysates. We developed a novel bioreactor design to generate separate sessile and planktonic cell populations for omics studies. Sessile cells had significantly greater expression of genes involved in catabolism of carbohydrates by glycolysis and pyruvate fermentation, ATP generation by proton gradient, the anabolism of proteins and lipids and cellular functions critical for cell division consistent with substrate replete conditions. Planktonic cells had notably higher gene expression for flagellar motility and chemotaxis, cellulosomal cellulases and anchoring scaffoldins, and a range of stress induced homeostasis mechanisms such as oxidative stress protection by antioxidants and flavoprotein co-factors, methionine repair, Fe-S cluster assembly and repair in redox proteins, cell growth control through tRNA thiolation, recovery of damaged DNA by nucleotide excision repair and removal of terminal proteins by proteases. This study demonstrates that microbial attachment to cellulose substrate produces widespread gene expression changes for critical functions of this organism and provides physiological insights for two cells populations relevant for engineering of industrially-ready phenotypes.

Improved growth rate in Clostridium thermocellum hydrogenase mutant via perturbed sulfur metabolism

BACKGROUND

Metabolic engineering is a commonly used approach to develop organisms for an industrial function, but engineering aimed at improving one phenotype can negatively impact other phenotypes. This lack of robustness can prove problematic. Cellulolytic bacterium Clostridium thermocellum is able to rapidly ferment cellulose to ethanol and other products. Recently, genes involved in H2 production, including the hydrogenase maturase hydG and NiFe hydrogenase ech, were deleted from the chromosome of C. thermocellum. While ethanol yield increased, the growth rate of ΔhydG decreased substantially compared to wild type.

RESULTS

Addition of 5 mM acetate to the growth medium improved the growth rate in C. thermocellum ∆hydG, whereas wild type remained unaffected. Transcriptomic analysis of the wild type showed essentially no response to the addition of acetate. However, in C. thermocellum ΔhydG, 204 and 56 genes were significantly differentially regulated relative to wild type in the absence and presence of acetate, respectively. Genes, Clo1313_0108-0125, which are predicted to encode a sulfate transport system and sulfate assimilatory pathway, were drastically upregulated in C. thermocellum ΔhydG in the presence of added acetate. A similar pattern was seen with proteomics. Further physiological characterization demonstrated an increase in sulfide synthesis and elimination of cysteine consumption in C. thermocellum ΔhydG. Clostridium thermocellum ΔhydGΔech had a higher growth rate than ΔhydG in the absence of added acetate, and a similar but less pronounced transcriptional and physiological effect was seen in this strain upon addition of acetate.

CONCLUSIONS

Sulfur metabolism is perturbed in C. thermocellum ΔhydG strains, likely to increase flux through sulfate reduction to act either as an electron sink to balance redox reactions or to offset an unknown deficiency in sulfur assimilation.

Genome-Scale Modeling of Thermophilic Microorganisms

Thermophilic microorganisms are of increasing interest for many industries as their enzymes and metabolisms are highly efficient at elevated temperatures. However, their metabolic processes are often largely different from their mesophilic counterparts. These differences can lead to metabolic engineering strategies that are doomed to fail. Genome-scale metabolic modeling is an effective and highly utilized way to investigate cellular phenotypes and to test metabolic engineering strategies. In this review we chronicle a number of thermophilic organisms that have recently been studied with genome-scale models. The microorganisms spread across archaea and bacteria domains, and their study gives insights that can be applied in a broader context than just the species they describe. We end with a perspective on the future development and applications of genome-scale models of thermophilic organisms.

Glycolysis without pyruvate kinase in Clostridium thermocellum

The metabolism of Clostridium thermocellum is notable in that it assimilates sugar via the EMP pathway but does not possess a pyruvate kinase enzyme. In the wild type organism, there are three proposed pathways for conversion of phosphoenolpyruvate (PEP) to pyruvate, which differ in their cofactor usage. One path uses pyruvate phosphate dikinase (PPDK), another pathway uses the combined activities of PEP carboxykinase (PEPCK) and oxaloacetate decarboxylase (ODC). Yet another pathway, the malate shunt, uses the combined activities of PEPCK, malate dehydrogenase and malic enzyme. First we showed that there is no flux through the ODC pathway by enzyme assay. Flux through the remaining two pathways (PPDK and malate shunt) was determined by dynamic ¹³C labeling. In the wild-type strain, the malate shunt accounts for about 33±2% of the flux to pyruvate, with the remainder via the PPDK pathway. Deletion of the ppdk gene resulted in a redirection of all pyruvate flux through the malate shunt. This provides the first direct evidence of the in-vivo function of the malate shunt.

Engineering electron metabolism to increase ethanol production in Clostridium thermocellum

The NfnAB (NADH-dependent reduced ferredoxin: NADP⁺ oxidoreductase) and Rnf (ion-translocating reduced ferredoxin: NAD⁺ oxidoreductase) complexes are thought to catalyze electron transfer between reduced ferredoxin and NAD(P)⁺. Efficient electron flux is critical for engineering fuel production pathways, but little is known about the relative importance of these enzymes in vivo. In this study we investigate the importance of the NfnAB and Rnf complexes in Clostridium thermocellum for growth on cellobiose and Avicel using gene deletion, enzyme assays, and fermentation product analysis. The NfnAB complex does not seem to play a major role in metabolism, since deletion of nfnAB genes had little effect on the distribution of fermentation products. By contrast, the Rnf complex appears to play an important role in ethanol formation. Deletion of rnf genes resulted in a decrease in ethanol formation. Overexpression of rnf genes resulted in an increase in ethanol production of about 30%, but only in strains where the hydG hydrogenase maturation gene was also deleted.

Exploring complex cellular phenotypes and model-guided strain design with a novel genome-scale metabolic model of Clostridium thermocellum DSM 1313 implementing an adjustable cellulosome

BACKGROUND

Clostridium thermocellum is a gram-positive thermophile that can directly convert lignocellulosic material into biofuels. The metabolism of C. thermocellum contains many branches and redundancies which limit biofuel production, and typical genetic techniques are time-consuming. Further, the genome sequence of a genetically tractable strain C. thermocellum DSM 1313 has been recently sequenced and annotated. Therefore, developing a comprehensive, predictive, genome-scale metabolic model of DSM 1313 is desired for elucidating its complex phenotypes and facilitating model-guided metabolic engineering.

RESULTS

We constructed a genome-scale metabolic model iAT601 for DSM 1313 using the KEGG database as a scaffold and an extensive literature review and bioinformatic analysis for model refinement. Next, we used several sets of experimental data to train the model, e.g., estimation of the ATP requirement for growth-associated maintenance (13.5 mmol ATP/g DCW/h) and cellulosome synthesis (57 mmol ATP/g cellulosome/h). Using our tuned model, we investigated the effect of cellodextrin lengths on cell yields, and could predict in silico experimentally observed differences in cell yield based on which cellodextrin species is assimilated. We further employed our tuned model to analyze the experimentally observed differences in fermentation profiles (i.e., the ethanol to acetate ratio) between cellobiose- and cellulose-grown cultures and infer regulatory mechanisms to explain the phenotypic differences. Finally, we used the model to design over 250 genetic modification strategies with the potential to optimize ethanol production, 6155 for hydrogen production, and 28 for isobutanol production.

CONCLUSIONS

Our developed genome-scale model iAT601 is capable of accurately predicting complex cellular phenotypes under a variety of conditions and serves as a high-quality platform for model-guided strain design and metabolic engineering to produce industrial biofuels and chemicals of interest.

Metabolic modeling of clostridia: current developments and applications

Anaerobic Clostridium spp. is an important bioproduction microbial genus that can produce solvents and utilize a broad spectrum of substrates including cellulose and syngas. Genome-scale metabolic (GSM) models are increasingly being put forth for various clostridial strains to explore their respective metabolic capabilities and suitability for various bioconversions. In this study, we have selected representative GSM models for six different clostridia (Clostridium acetobutylicum, C. beijerinckii, C. butyricum, C. cellulolyticum, C. ljungdahlii and C. thermocellum) and performed a detailed model comparison contrasting their metabolic repertoire. We also discuss various applications of these GSM models to guide metabolic engineering interventions as well as assessing cellular physiology.