Integrative modelling of pH-dependent enzyme activity and transcriptomic regulation of the acetone-butanol-ethanol fermentation of Clostridium acetobutylicum in continuous culture

Mathematical modeling of ABE fermentation on glucose substrate with Zn supplementation for enhanced butanol production

Supplementation or limitation of some micronutrients during acetone-butanol-ethanol (ABE) fermentation has led to improvement in butanol yield and productivity. A mechanistic model of ABE fermentation offers insights in understanding these complex interactions and improving productivity through optimal culture conditions. This study proposes a mechanistic kinetic model of ABE fermentation by two Clostridium Acetobutylicum strains, L7 and ATCC 824 using glucose as sole carbon source without zinc and with various zinc doses. The model incorporates enzyme regulation by zinc on several glycolytic, acidogenesis and solventogenesis enzymes. The model was fitted and validated to experimental data collected from the published literature. The simulated results were in compliance with the experimental data, most importantly indicating higher glucose consumption and butanol productivity when supplemented with zinc compared to the control culture. The average squared correlation factor (R2) between the experimental and the simulated results, without and with zinc, were 0.99 and 0.96 for glucose, and 0.89 and 0.95 for butanol, respectively. A sensitivity analysis performed on the fitted and validated model indicated that the glucose consumption and growth parameters most influenced the model outputs. The developed model can be used as a template for modeling ABE fermentation under different combinations of micronutrients that may offer improved butanol yield and productivity.

The pH-dependent lactose metabolism of Lactobacillus delbrueckii subsp. bulgaricus: An integrative view through a mechanistic computational model

The fermentation process of milk to yoghurt using Lactobacillus delbrueckii subsp. bulgaricus in co-culture with Streptococcus thermophilus is hallmarked by the breakdown of lactose to organic acids such as lactate. This leads to a substantial decrease in pH - both in the medium, as well as cytosolic. The latter impairs metabolic activities due to the pH-dependence of enzymes, which compromises microbial growth. To quantitatively elucidate the impact of the acidification on metabolism of L. bulgaricus in an integrated way, we have developed a proton-dependent computational model of lactose metabolism and casein degradation based on experimental data. The model accounts for the influence of pH on enzyme activities as well as cellular growth and proliferation of the bacterial population. We used a machine learning approach to quantify the cell volume throughout fermentation. Simulation results show a decrease in metabolic flux with acidification of the cytosol. Additionally, the validated model predicts a similar metabolic behaviour within a wide range of non-limiting substrate concentrations. This computational model provides a deeper understanding of the intricate relationships between metabolic activity and acidification and paves the way for further optimization of yoghurt production under industrial settings.

Exploring the Influence of pH on the Dynamics of Acetone-Butanol-Ethanol Fermentation

Clostridium acetobutylicum is an anaerobic bacterium that is extensively studied for its ability to produce butanol. Over the past two decades, various genetic and metabolic engineering approaches have been used to investigate the physiology and regulation system of the biphasic metabolic pathway in this organism. However, there has been a relatively limited amount of research focused on the fermentation dynamics of C. acetobutylicum. In this study, we developed a pH-based phenomenological model to predict the fermentative production of butanol from glucose using C. acetobutylicum in a batch system. The model describes the relationship between the dynamics of growth and the production of desired metabolites and the extracellular pH of the media. Our model was found to be successful in predicting the fermentation dynamics of C. acetobutylicum, and the simulations were validated using experimental fermentation data. Furthermore, the proposed model has the potential to be extended to represent the dynamics of butanol production in other fermentation systems, such as fed-batch or continuous fermentation using single and multi-sugars.

pH dependencies of glycolytic enzymes of yeast under in vivo-like assay conditions

Under carbon source transitions, the intracellular pH of Saccharomyces cerevisiae is subject to change. Dynamics in pH modulate the activity of the glycolytic enzymes, resulting in a change in glycolytic flux and ultimately cell growth. To understand how pH affects the global behavior of glycolysis and ethanol fermentation, we measured the activity of the glycolytic and fermentative enzymes in S. cerevisiae under in vivo-like conditions at different pH. We demonstrate that glycolytic enzymes exhibit differential pH dependencies, and optima, in the pH range observed during carbon source transitions. The forward reaction of GAPDH shows the highest decrease in activity, 83%, during a simulated feast/famine regime upon glucose removal (cytosolic pH drop from 7.1 to 6.4). We complement our biochemical characterization of the glycolytic enzymes by fitting the V_max to the progression curves of product formation or decay over time. The fitting analysis shows that the observed changes in enzyme activities require changes in V_max , but changes in K_m cannot be excluded. Our study highlights the relevance of pH as a key player in metabolic regulation and provides a large set of quantitative data that can be explored to improve our understanding of metabolism in dynamic environments.

Kinetic Modeling of Saccharomyces cerevisiae Central Carbon Metabolism: Achievements, Limitations, and Opportunities

Central carbon metabolism comprises the metabolic pathways in the cell that process nutrients into energy, building blocks and byproducts. To unravel the regulation of this network upon glucose perturbation, several metabolic models have been developed for the microorganism Saccharomyces cerevisiae. These dynamic representations have focused on glycolysis and answered multiple research questions, but no commonly applicable model has been presented. This review systematically evaluates the literature to describe the current advances, limitations, and opportunities. Different kinetic models have unraveled key kinetic glycolytic mechanisms. Nevertheless, some uncertainties regarding model topology and parameter values still limit the application to specific cases. Progressive improvements in experimental measurement technologies as well as advances in computational tools create new opportunities to further extend the model scale. Notably, models need to be made more complex to consider the multiple layers of glycolytic regulation and external physiological variables regulating the bioprocess, opening new possibilities for extrapolation and validation. Finally, the onset of new data representative of individual cells will cause these models to evolve from depicting an average cell in an industrial fermenter, to characterizing the heterogeneity of the population, opening new and unseen possibilities for industrial fermentation improvement.

Transcriptomic studies of solventogenic clostridia, Clostridium acetobutylicum and Clostridium beijerinckii

Solventogenic clostridia are not a strictly defined group within the genus Clostridium but its representatives share some common features, i.e. they are anaerobic, non-pathogenic, non-toxinogenic and endospore forming bacteria. Their main metabolite is typically 1-butanol but depending on species and culture conditions, they can form other metabolites such as acetone, isopropanol, ethanol, butyric, lactic and acetic acids, and hydrogen. Although these organisms were previously used for the industrial production of solvents, they later fell into disuse, being replaced by more efficient chemical production. A return to a more biological production of solvents therefore requires a thorough understanding of clostridial metabolism. Transcriptome analysis, which reflects the involvement of individual genes in all cellular processes within a population, at any given (sampling) moment, is a valuable tool for gaining a deeper insight into clostridial life. In this review, we describe techniques to study transcription, summarize the evolution of these techniques and compare methods for data processing and visualization of solventogenic clostridia, particularly the species Clostridium acetobutylicum and Clostridium beijerinckii. Individual approaches for evaluating transcriptomic data are compared and their contributions to advancements in the field are assessed. Moreover, utilization of transcriptomic data for reconstruction of computational clostridial metabolic models is considered and particular models are described. Transcriptional changes in glucose transport, central carbon metabolism, the sporulation cycle, butanol and butyrate stress responses, the influence of lignocellulose-derived inhibitors on growth and solvent production, and other respective topics, are addressed and common trends are highlighted.

Crystal Structure and Molecular Mechanism of Phosphotransbutyrylase from Clostridium acetobutylicum

Acetone-butanol-ethanol (ABE) fermentation by the anaerobic bacterium Clostridium acetobutylicum has been considered a promising process of industrial biofuel production. Phosphotransbutyrylase (phosphate butyryltransferase, PTB) plays a crucial role in butyrate metabolism by catalyzing the reversible conversion of butyryl-CoA into butyryl phosphate. Here, we report the crystal structure of PTB from the Clostridial host for ABE fermentation, C. acetobutylicum, (CaPTB) at a 2.9 Å resolution. The overall structure of the CaPTB monomer is quite similar to those of other acyltransferases, with some regional structural differences. The monomeric structure of CaPTB consists of two distinct domains, the N- and C-terminal domains. The active site cleft was formed at the interface between the two domains. Interestingly, the crystal structure of CaPTB contained eight molecules per asymmetric unit, forming an octamer, and the size-exclusion chromatography experiment also suggested that the enzyme exists as an octamer in solution. The structural analysis of CaPTB identifies the substrate binding mode of the enzyme and comparisons with other acyltransferase structures lead us to speculate that the enzyme undergoes a conformational change upon binding of its substrate.

An integrated modular framework for modeling the effect of ammonium on the sialylation process of monoclonal antibodies produced by CHO cells

BACKGROUND

Monoclonal antibodies (mABs) have emerged as one of the most important therapeutic recombinant proteins in the pharmaceutical industry. Their immunogenicity and therapeutic efficacy are influenced by post-translational modifications, specifically the glycosylation process. Bioprocess conditions can influence the intracellular process of glycosylation. Among all the process conditions that have been recognized to affect the mAB glycoforms, the detailed mechanism underlying how ammonium could perturb glycosylation remains to be fully understood. It was shown that ammonium induces heterogeneity in protein glycosylation by altering the sialic acid content of glycoproteins. Hence, understanding this mechanism would aid pharmaceutical manufacturers to ensure consistent protein glycosylation.

METHODS

Three different mechanisms have been proposed to explain how ammonium influences the sialylation process. In the first, the inhibition of CMP-sialic acid transporter, which transports CMP-sialic acid (sialylation substrate) into the Golgi, by an increase in UDP-GlcNAc content that is brought about by the augmented incorporation of ammonium into glucosamine formation. In the second, ammonia diffuses into the Golgi and raises its pH, thereby decreasing the sialyltransferase enzyme activity. In the third, the reduction of sialyltransferase enzyme expression level in the presence of ammonium. We employed these mechanisms in a novel integrated modular platform to link dynamic alteration in mAB sialylation process with extracellular ammonium concentration to elucidate how ammonium alters the sialic acid content of glycoproteins.

RESULTS

Our results show that the sialylation reaction rate is insensitive to the first mechanism. At low ammonium concentration, the second mechanism is the controlling mechanism in mAB sialylation and by increasing the ammonium level (< 8 mM) the third mechanism becomes the controlling mechanism. At higher ammonium concentrations (> 8 mM) the second mechanism becomes predominant again.

CONCLUSION

The presented model in this study provides a connection between extracellular ammonium and the monoclonal antibody sialylation process. This computational tool could help scientists to develop and formulate cell culture media. The model illustrated here can assist the researchers to select culture media that ensure consistent mAB sialylation.

Protein Acetylation and Butyrylation Regulate the Phenotype and Metabolic Shifts of the Endospore-forming Clostridium acetobutylicum

Clostridium acetobutylicum is a strict anaerobic, endospore-forming bacterium, which is used for the production of the high energy biofuel butanol in metabolic engineering. The life cycle of C. acetobutylicum can be divided into two phases, with acetic and butyric acids being produced in the exponential phase (acidogenesis) and butanol formed in the stationary phase (solventogenesis). During the transitional phase from acidogenesis to solventogenesis and latter stationary phase, concentration peaks of the metabolic intermediates butyryl phosphate and acetyl phosphate are observed. As an acyl group donor, acyl-phosphate chemically acylates protein substrates. However, the regulatory mechanism of lysine acetylation and butyrylation involved in the phenotype and solventogenesis of C. acetobutylicum remains unknown. In our study, we conducted quantitative analysis of protein acetylome and butyrylome to explore the dynamic change of lysine acetylation and butyrylation in the exponential phase, transitional phase, and stationary phase of C. acetobutylicum Total 458 lysine acetylation sites and 1078 lysine butyrylation sites were identified in 254 and 373 substrates, respectively. Bioinformatics analysis uncovered the similarities and differences between the two acylation modifications in C. acetobutylicum Mutation analysis of butyrate kinase and the central transcriptional factor Spo0A was performed to characterize the unique role of lysine butyrylation in the metabolic pathway and sporulation process of C. acetobutylicum Moreover, quantitative proteomic assays were performed to reveal the relationship between protein features (e.g. gene expression level and lysine acylation level) and metabolites in the three growth stages. This study expanded our knowledge of lysine acetylation and butyrylation in Clostridia and constituted a resource for functional studies on lysine acylation in bacteria.

Remarkable reactivity of alkoxide/acetato-bridged binuclear copper(II) complex as artificial carboxylesterase

Bromo-containing binuclear Schiff base copper(II) complex, Cu₂L(OAc), with an alkoxo/acetato-bridged moiety was employed as a model of carboxylesterases to promote the hydrolytic cleavage of p-nitrophenyl picolinate (PNPP). Furthermore, the reactivity of a mononuclear complex (CuHL) was evaluated for comparing it with that of binuclear one. The results reveal that the as-prepared binuclear Cu₂L(OAc) efficiently accelerated the hydrolysis of PNPP, giving rise to excess four orders of magnitude rate enhancement in contrast to the un-catalyzed reaction. Cu₂L(OAc) represented an enzyme-like bell-shaped pH-responsive kinetic behavior. Moreover, the binuclear one is more reactive than its mononuclear analogue (CuHL) by two orders of magnitude. The total efficiency of Cu₂L(OAc) is about 61-fold than that of its mononuclear analogue, CuHL. In addition, a contrast experiment reveals that binuclear Cu₂L(OAc) displayed good activity in the hydrolysis of PNPP as well another active ester, i.e., S-2-benzothiazolyl 2-amino-alpha-(methoxyimino)-4-thiazolethiolacetate (AE-active ester). Noteworthyly, it was found that mononuclear one inspired more obvious rate enhancement in the hydrolysis of AE-active ester relative to PNPP hydrolysis. The estimated pK _a1 of bound water on the binuclear Cu₂L(OAc) using second derivative method (SDM) is relatively smaller than that for CuHL by a gap of about 0.8 pK unit, which facilitates the hydrolysis of PNPP. Four orders of magnitude rate enhancement was observed for the catalytic hydrolysis of p-nitrophenyl picolinate (PNPP) by one μ-alkoxide/acetato-bridged binuclear copper(II) complex under physiological conditions. Substrate specificity of the resulting binuclear complexes was observed for the hydrolysis of PNPP and AE-active ester.

Mathematical modelling of clostridial acetone-butanol-ethanol fermentation

Clostridial acetone-butanol-ethanol (ABE) fermentation features a remarkable shift in the cellular metabolic activity from acid formation, acidogenesis, to the production of industrial-relevant solvents, solventogensis. In recent decades, mathematical models have been employed to elucidate the complex interlinked regulation and conditions that determine these two distinct metabolic states and govern the transition between them. In this review, we discuss these models with a focus on the mechanisms controlling intra- and extracellular changes between acidogenesis and solventogenesis. In particular, we critically evaluate underlying model assumptions and predictions in the light of current experimental knowledge. Towards this end, we briefly introduce key ideas and assumptions applied in the discussed modelling approaches, but waive a comprehensive mathematical presentation. We distinguish between structural and dynamical models, which will be discussed in their chronological order to illustrate how new biological information facilitates the ‘evolution’ of mathematical models. Mathematical models and their analysis have significantly contributed to our knowledge of ABE fermentation and the underlying regulatory network which spans all levels of biological organization. However, the ties between the different levels of cellular regulation are not well understood. Furthermore, contradictory experimental and theoretical results challenge our current notion of ABE metabolic network structure. Thus, clostridial ABE fermentation still poses theoretical as well as experimental challenges which are best approached in close collaboration between modellers and experimentalists.

Kinetic Study of Acetone-Butanol-Ethanol Fermentation in Continuous Culture

Acetone-butanol-ethanol (ABE) fermentation by clostridia has shown promise for industrial-scale production of biobutanol. However, the continuous ABE fermentation suffers from low product yield, titer, and productivity. Systems analysis of the continuous ABE fermentation will offer insights into its metabolic pathway as well as into optimal fermentation design and operation. For the ABE fermentation in continuous Clostridium acetobutylicum culture, this paper presents a kinetic model that includes the effects of key metabolic intermediates and enzymes as well as culture pH, product inhibition, and glucose inhibition. The kinetic model is used for elucidating the behavior of the ABE fermentation under the conditions that are most relevant to continuous cultures. To this end, dynamic sensitivity analysis is performed to systematically investigate the effects of culture conditions, reaction kinetics, and enzymes on the dynamics of the ABE production pathway. The analysis provides guidance for future metabolic engineering and fermentation optimization studies.

System-level modeling of acetone-butanol-ethanol fermentation

Acetone-butanol-ethanol (ABE) fermentation is a metabolic process of clostridia that produces bio-based solvents including butanol. It is enabled by an underlying metabolic reaction network and modulated by cellular gene regulation and environmental cues. Mathematical modeling has served as a valuable strategy to facilitate the understanding, characterization and optimization of this process. In this review, we highlight recent advances in system-level, quantitative modeling of ABE fermentation. We begin with an overview of integrative processes underlying the fermentation. Next we survey modeling efforts including early simple models, models with a systematic metabolic description, and those incorporating metabolism through simple gene regulation. Particular focus is given to a recent system-level model that integrates the metabolic reactions, gene regulation and environmental cues. We conclude by discussing the remaining challenges and future directions towards predictive understanding of ABE fermentation.

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.

Redox-switch regulatory mechanism of thiolase from Clostridium acetobutylicum

Thiolase is the first enzyme catalysing the condensation of two acetyl-coenzyme A (CoA) molecules to form acetoacetyl-CoA in a dedicated pathway towards the biosynthesis of n-butanol, an important solvent and biofuel. Here we elucidate the crystal structure of Clostridium acetobutylicum thiolase (CaTHL) in its reduced/oxidized states. CaTHL, unlike those from other aerobic bacteria such as Escherichia coli and Zoogloea ramegera, is regulated by the redox-switch modulation through reversible disulfide bond formation between two catalytic cysteine residues, Cys88 and Cys378. When CaTHL is overexpressed in wild-type C. acetobutylicum, butanol production is reduced due to the disturbance of acidogenic to solventogenic shift. The CaTHL(V77Q/N153Y/A286K) mutant, which is not able to form disulfide bonds, exhibits higher activity than wild-type CaTHL, and enhances butanol production upon overexpression. On the basis of these results, we suggest that CaTHL functions as a key enzyme in the regulation of the main metabolism of C. acetobutylicum through a redox-switch regulatory mechanism.

Redox-switch regulatory mechanism of thiolase from Clostridium acetobutylicum

Thiolase is the first enzyme catalysing the condensation of two acetyl-coenzyme A (CoA) molecules to form acetoacetyl-CoA in a dedicated pathway towards the biosynthesis of n-butanol, an important solvent and biofuel. Here we elucidate the crystal structure of Clostridium acetobutylicum thiolase (CaTHL) in its reduced/oxidized states. CaTHL, unlike those from other aerobic bacteria such as Escherichia coli and Zoogloea ramegera, is regulated by the redox-switch modulation through reversible disulfide bond formation between two catalytic cysteine residues, Cys88 and Cys378. When CaTHL is overexpressed in wild-type C. acetobutylicum, butanol production is reduced due to the disturbance of acidogenic to solventogenic shift. The CaTHL(V77Q/N153Y/A286K) mutant, which is not able to form disulfide bonds, exhibits higher activity than wild-type CaTHL, and enhances butanol production upon overexpression. On the basis of these results, we suggest that CaTHL functions as a key enzyme in the regulation of the main metabolism of C. acetobutylicum through a redox-switch regulatory mechanism.

Coenzyme A-transferase-independent butyrate re-assimilation in Clostridium acetobutylicum-evidence from a mathematical model

The hetero-dimeric CoA-transferase CtfA/B is believed to be crucial for the metabolic transition from acidogenesis to solventogenesis in Clostridium acetobutylicum as part of the industrial-relevant acetone-butanol-ethanol (ABE) fermentation. Here, the enzyme is assumed to mediate re-assimilation of acetate and butyrate during a pH-induced metabolic shift and to faciliate the first step of acetone formation from acetoacetyl-CoA. However, recent investigations using phosphate-limited continuous cultures have questioned this common dogma. To address the emerging experimental discrepancies, we investigated the mutant strain Cac-ctfA398s::CT using chemostat cultures. As a consequence of this mutation, the cells are unable to express functional ctfA and are thus lacking CoA-transferase activity. A mathematical model of the pH-induced metabolic shift, which was recently developed for the wild type, is used to analyse the observed behaviour of the mutant strain with a focus on re-assimilation activities for the two produced acids. Our theoretical analysis reveals that the ctfA mutant still re-assimilates butyrate, but not acetate. Based upon this finding, we conclude that C. acetobutylicum possesses a CoA-tranferase-independent butyrate uptake mechanism that is activated by decreasing pH levels. Furthermore, we observe that butanol formation is not inhibited under our experimental conditions, as suggested by previous batch culture experiments. In concordance with recent batch experiments, acetone formation is abolished in chemostat cultures using the ctfa mutant.

Gene duplication and the environmental regulation of physiology and development

When different life stages have different environmental tolerances, development needs to be regulated so that each life stage experiences environmental conditions that are suitable for it, if fitness is to be maintained. Restricting the timing of developmental transitions to occur under specific combinations of environmental conditions is therefore adaptively important. However, impeding development can itself incur demographic and fitness costs. How do organisms regulate development and physiological processes so that they occur under the broadest range of permissive conditions? Gene duplication offers one solution: Multiple genes contribute to the same downstream process, but do so under distinct combinations of environmental conditions. We present a simple model to examine how environmental sensitivities of genes and how gene duplication influence the distribution of environmental conditions under which an end process will proceed. The model shows that the duplication of genes that retain their downstream function but diverge in environmental sensitivities can allow an end process to proceed under more than one distinct combination of environmental conditions. The outcomes depend on how upstream genes regulate downstream components, which genes in the pathway have diversified in their sensitivities, and the structure of the pathway.

Proteomic analyses of the phase transition from acidogenesis to solventogenesis using solventogenic and non-solventogenic Clostridium acetobutylicum strains

The fermentation carried out by the solvent-producing bacterium, Clostridium acetobutylicum, is characterized by two distinct phases: acidogenic and solventogenic phases. Understanding the cellular physiological changes occurring during the phase transition in clostridial fermentation is important for the enhanced production of solvents. To identify protein changes upon entry to stationary phase where solvents are typically produced, we herein analyzed the proteomic profiles of the parental wild type C. acetobutylicum strains, ATCC 824, the non-solventogenic strain, M5 that has lost the solventogenic megaplasmid pSOL1, and the synthetic simplified alcohol forming strain, M5 (pIMP1E1AB) expressing plasmid-based CoA-transferase (CtfAB) and aldehyde/alcohol dehydrogenase (AdhE1). A total of 68 protein spots, corresponding to 56 unique proteins, were unambiguously identified as being differentially present after the phase transitions in the three C. acetobutylicum strains. In addition to changes in proteins known to be involved in solventogenesis (AdhE1 and CtfB), we identified significant alterations in enzymes involved in sugar transport and metabolism, fermentative pathway, heat shock proteins, translation, and amino acid biosynthesis upon entry into the stationary phase. Of these, four increased proteins (AdhE1, CAC0233, CtfB and phosphocarrier protein HPr) and six decreased proteins (butyrate kinase, ferredoxin:pyruvate oxidoreductase, phenylalanyl-tRNA synthetase, adenylosuccinate synthase, pyruvate kinase and valyl-tRNA synthetase) showed similar patterns in the two strains capable of butanol formation. Interestingly, significant changes of several proteins by post-translational modifications were observed in the solventogenic phase. The proteomic data from this study will improve our understanding on how cell physiology is affected through protein levels patterns in clostridia.