Physicochemical and metabolic constraints for thermodynamics-based stoichiometric modelling under mesophilic growth conditions.

Metabolic engineering in the post-genomic era is characterised by the development of new methods for metabolomics and fluxomics, supported by the integration of genetic engineering tools and mathematical modelling. Particularly, constraint-based stoichiometric models have been widely studied: (i) flux balance analysis (FBA) (in silico), and (ii) metabolic flux analysis (MFA) (in vivo). Recent studies have enabled the incorporation of thermodynamics and metabolomics data to improve the predictive capabilities of these approaches. However, an in-depth comparison and evaluation of these methods is lacking. This study presents a thorough analysis of two different in silico methods tested against experimental data (metabolomics and 13C-MFA) for the mesophile Escherichia coli. In particular, a modified version of the recently published matTFA toolbox was created, providing a broader range of physicochemical parameters. Validating against experimental data allowed the determination of the best physicochemical parameters to perform the TFA (Thermodynamics-based Flux Analysis). An analysis of flux pattern changes in the central carbon metabolism between 13C-MFA and TFA highlighted the limited capabilities of both approaches for elucidating the anaplerotic fluxes. In addition, a method based on centrality measures was suggested to identify important metabolites that (if quantified) would allow to further constrain the TFA. Finally, this study emphasised the need for standardisation in the fluxomics community: novel approaches are frequently released but a thorough comparison with currently accepted methods is not always performed.

Improved CRISPR/Cas9 Tools for the Rapid Metabolic Engineering of Clostridium acetobutylicum.

Clustered regularly interspaced short palindromic repeats (CRISPR)/Cas (CRISPR-associated proteins)9 tools have revolutionized biology-several highly efficient tools have been constructed that have resulted in the ability to quickly engineer model bacteria, for example, Escherichia coli. However, the use of CRISPR/Cas9 tools has lagged behind in non-model bacteria, hampering engineering efforts. Here, we developed improved CRISPR/Cas9 tools to enable efficient rapid metabolic engineering of the industrially relevant bacterium Clostridium acetobutylicum. Previous efforts to implement a CRISPR/Cas9 system in C. acetobutylicum have been hampered by the lack of tightly controlled inducible systems along with large plasmids resulting in low transformation efficiencies. We successfully integrated the cas9 gene from Streptococcuspyogenes into the genome under control of the xylose inducible system from Clostridium difficile, which we then showed resulted in a tightly controlled system. We then optimized the length of the editing cassette, resulting in a small editing plasmid, which also contained the upp gene in order to rapidly lose the plasmid using the upp/5-fluorouracil counter-selection system. We used this system to perform individual and sequential deletions of ldhA and the ptb-buk operon.

Synthetic Biology on Acetogenic Bacteria for Highly Efficient Conversion of C1 Gases to Biochemicals.

Synthesis gas, which is mainly produced from fossil fuels or biomass gasification, consists of C1 gases such as carbon monoxide, carbon dioxide, and methane as well as hydrogen. Acetogenic bacteria (acetogens) have emerged as an alternative solution to recycle C1 gases by converting them into value-added biochemicals using the Wood-Ljungdahl pathway. Despite the advantage of utilizing acetogens as biocatalysts, it is difficult to develop industrial-scale bioprocesses because of their slow growth rates and low productivities. To solve these problems, conventional approaches to metabolic engineering have been applied; however, there are several limitations owing to the lack of required genetic bioparts for regulating their metabolic pathways. Recently, synthetic biology based on genetic parts, modules, and circuit design has been actively exploited to overcome the limitations in acetogen engineering. This review covers synthetic biology applications to design and build industrial platform acetogens.

Roles of the F-domain in [FeFe] hydrogenase.

The role of accessory Fe-S clusters of the F-domain in the catalytic activity of M3-type [FeFe] hydrogenase and the contribution of each of the two Fe-S surface clusters in the intermolecular electron transfer from ferredoxin are both poorly understood. We designed, constructed, produced and spectroscopically, electrochemically and biochemically characterized three mutants of Clostridium acetobutylicum CaHydA hydrogenase with modified Fe-S clusters: two site-directed mutants, HydA_C100A and HydA_C48A missing the FS4C and the FS2 surface Fe-S clusters, respectively, and a HydA_ΔDA mutant that completely lacks the F-domain. Analysis of the mutant enzyme activities clearly demonstrated the importance of accessory clusters in retaining full enzyme activity at potentials around and higher than the equilibrium 2H+/H2 potential but not at the lowest potentials, where all enzymes have a similar turnover rate. Moreover, our results, combined with molecular modelling approaches, indicated that the FS2 cluster is the main gate for electron transfer from reduced ferredoxin.

Metabolic flexibility of a butyrate pathway mutant of Clostridium acetobutylicum.

Clostridium acetobutylicum possesses two homologous buk genes, buk (or buk1) and buk2, which encode butyrate kinases involved in the last step of butyrate formation. To investigate the contribution of buk in detail, an in-frame deletion mutant was constructed. However, in all the Δbuk mutants obtained, partial deletions of the upstream ptb gene were observed, and low phosphotransbutyrylase and butyrate kinase activities were measured. This demonstrates that i) buk (CA_C3075) is the key butyrate kinase-encoding gene and that buk2 (CA_C1660) that is poorly transcribed only plays a minor role; and ii) strongly suggests that a Δbuk mutant is not viable if the ptb gene is not also inactivated, probably due to the accumulation of butyryl-phosphate, which might be toxic for the cell. One of the ΔbukΔptb mutants was subjected to quantitative transcriptomic (mRNA molecules/cell) and fluxomic analyses in acidogenic, solventogenic and alcohologenic chemostat cultures. In addition to the low butyrate production, drastic changes in metabolic fluxes were also observed for the mutant: i) under acidogenic conditions, the primary metabolite was butanol and a new metabolite, 2-hydroxy-valerate, was produced ii) under solventogenesis, 58% increased butanol production was obtained compared to the control strain under the same conditions, and a very high yield of butanol formation (0.3gg-1) was reached; and iii) under alcohologenesis, the major product was lactate. Furthermore, at the transcriptional level, adhE2, which encodes an aldehyde/alcohol dehydrogenase and is known to be a gene specifically expressed in alcohologenesis, was surprisingly highly expressed in all metabolic states in the mutant. The results presented here not only support the key roles of buk and ptb in butyrate formation but also highlight the metabolic flexibility of C. acetobutylicum in response to genetic alteration of its primary metabolism.

Reactivity of the Excited States of the H-Cluster of FeFe Hydrogenases.

FeFe hydrogenases catalyze H2 oxidation and formation at an inorganic active site (the "H-cluster"), which consists of a [Fe2(CO)3(CN)2(dithiomethylamine)] subcluster covalently attached to a Fe4S4 subcluster. This active site is photosensitive: visible light has been shown to induce the release of exogenous CO (a reversible inhibitor of the enzyme), shuffle the intrinsic CO ligands, and even destroy the H-cluster. These reactions must be understood because they may negatively impact the use of hydrogenase for the photoproduction of H2. Here, we explore in great detail the reactivity of the excited states of the H-cluster under catalytic conditions by examining, both experimentally and using TDDFT calculations, the simplest photochemical reaction: the binding and release of exogenous CO. A simple dyad model can be used to predict which excitations are active. This strategy could be used for probing other aspects of the photoreactivity of the H-cluster.

Electrochemical Measurements of the Kinetics of Inhibition of Two FeFe Hydrogenases by O2 Demonstrate That the Reaction Is Partly Reversible.

The mechanism of reaction of FeFe hydrogenases with oxygen has been debated. It is complex, apparently very dependent on the details of the protein structure, and difficult to study using conventional kinetic techniques. Here we build on our recent work on the anaerobic inactivation of the enzyme [Fourmond et al. Nat. Chem. 2014, 4, 336-342] to propose and apply a new method for studying this reaction. Using electrochemical measurements of the turnover rate of hydrogenase, we could resolve the first steps of the inhibition reaction and accurately determine their rates. We show that the two most studied FeFe hydrogenases, from Chlamydomonas reinhardtii and Clostridium acetobutylicum, react with O2 according to the same mechanism, despite the fact that the former is much more O2 sensitive than the latter. Unlike often assumed, both enzymes are reversibly inhibited by a short exposure to O2. This will have to be considered to elucidate the mechanism of inhibition, before any prediction can be made regarding which mutations will improve oxygen resistance. We hope that the approach described herein will prove useful in this respect.

Stress-induced evolution of Escherichia coli points to original concepts in respiratory cofactor selectivity.

Bacterial metabolism is characterized by a remarkable capacity to rapidly adapt to environmental changes. We restructured the central metabolic network in Escherichia coli to force a higher production of NADPH, and then grew this strain in conditions favoring adaptive evolution. A six-fold increase in growth capacity was attained that could be attributed in multiple clones, after whole genome mutation mapping, to a specific single mutation. Each clone had an evolved NuoF*(E183A) enzyme in the respiratory complex I that can now oxidize both NADH and NADPH. When a further strain was constructed with an even higher degree of NADPH stress such that growth was impossible on glucose mineral medium, a solid-state screening for mutations restoring growth, led to two different types of NuoF mutations in strains having recovered growth capacity. In addition to the previously seen E183A mutation other clones showed a E183G mutation, both having NADH and NADPH oxidizing ability. These results demonstrate the unique solution used by E. coli to overcome the NADPH stress problem. This solution creates a new function for NADPH that is no longer restricted to anabolic synthesis reactions but can now be also used to directly produce catabolic energy.

Chromosomal integration of the pSOL1 megaplasmid of Clostridium acetobutylicum for continuous and stable advanced biofuels production.

Biofuel production by Clostridium acetobutylicum is compromised by strain degeneration due to loss of its pSOL1 megaplasmid. Here we used engineering biology to stably integrate pSOL1 into the chromosome together with a synthetic isopropanol pathway. In a membrane bioreactor continuously fed with glucose mineral medium, the final strain produced advanced biofuels, n-butanol and isopropanol, at high yield (0.31 g g-1), titre (15.4 g l-1) and productivity (15.5 g l-1 h-1) without degeneration.

Steady-state catalytic wave-shapes for 2-electron reversible electrocatalysts and enzymes.

Using direct electrochemistry to learn about the mechanism of electrocatalysts and redox enzymes requires that kinetic models be developed. Here we thoroughly discuss the interpretation of electrochemical signals obtained with adsorbed enzymes and molecular catalysts that can reversibly convert their substrate and product. We derive analytical relations between electrochemical observables (overpotentials for catalysis in each direction, positions, and magnitudes of the features of the catalytic wave) and the characteristics of the catalytic cycle (redox properties of the catalytic intermediates, kinetics of intramolecular and interfacial electron transfer, etc.). We discuss whether or not the position of the wave is determined by the redox potential of a redox relay when intramolecular electron transfer is slow. We demonstrate that there is no simple relation between the reduction potential of the active site and the catalytic bias of the enzyme, defined as the ratio of the oxidative and reductive limiting currents; this explains the recent experimental observation that the catalytic bias of NiFe hydrogenase depends on steps of the catalytic cycle that occur far from the active site [Abou Hamdan et al., J. Am. Chem. Soc. 2012, 134, 8368]. On the experimental side, we examine which models can best describe original data obtained with various NiFe and FeFe hydrogenases, and we illustrate how the presence of an intramolecular electron transfer chain affects the voltammetry by comparing the data obtained with the FeFe hydrogenases from Chlamydomonas reinhardtii and Clostridium acetobutylicum, only one of which has a chain of redox relays. The considerations herein will help the interpretation of electrochemical data previously obtained with various other bidirectional oxidoreductases, and, possibly, synthetic inorganic catalysts.

Metabolic engineering of Clostridium acetobutylicum ATCC 824 for the high-yield production of a biofuel composed of an isopropanol/butanol/ethanol mixture.

Clostridium acetobutylicum was metabolically engineered to produce a biofuel consisting of an isopropanol/butanol/ethanol mixture. For this purpose, different synthetic isopropanol operons were constructed and introduced on plasmids in a butyrate minus mutant strain (C. acetobutylicum ATCC 824 Δcac15ΔuppΔbuk). The best strain expressing the isopropanol operon from the thl promoter was selected from batch experiments at pH 5. By further optimizing the pH of the culture, a biofuel mixture with almost no by-products was produced at a titer, a yield and productivity never reached before, opening the opportunities to develop an industrial process for alternative biofuels with Clostridial species. Furthermore, by performing in vivo and in vitro flux analysis of the synthetic isopropanol pathway, this flux was identified to be limited by the [acetate](int) and the high Km of CoA-transferase for acetate. Decreasing the Km of this enzyme using a protein engineering approach would be a good target for improving isopropanol production and avoiding acetate accumulation in the culture medium.

Improvement of the Genome Editing Tools Based on 5FC/5FU Counter Selection in Clostridium acetobutylicum.

Several genetic tools have been developed for genome engineering in Clostridium acetobutylicum utilizing 5-fluorouracil (5FU) or 5-fluorocytosine (5FC) resistance as a selection method. In our group, a method based on the integration, by single crossing over, of a suicide plasmid (pCat-upp) followed by selection for the second crossing over using a counter-selectable marker (the upp gene and 5FU resistance) was recently developed for genome editing in C. acetobutylicum. This method allows genome modification without leaving any marker or scar in a strain of C. acetobutylicum that is ∆upp. Unfortunately, 5FU has strong mutagenic properties, inducing mutations in the strain's genome. After numerous applications of the pCat-upp/5FU system for genome modification in C. acetobutylicum, the CAB1060 mutant strain became entirely resistant to 5FU in the presence of the upp gene, resulting in failure when selecting on 5FU for the second crossing over. It was found that the potential repressor of the pyrimidine operon, PyrR, was mutated at position A115, leading to the 5FU resistance of the strain. To fix this problem, we created a corrective replicative plasmid expressing the pyrR gene, which was shown to restore the 5FU sensitivity of the strain. Furthermore, in order to avoid the occurrence of the problem observed with the CAB1060 strain, a preventive suicide plasmid, pCat-upp-pyrR*, was also developed, featuring the introduction of a synthetic codon-optimized pyrR gene, which was referred to as pyrR* with low nucleotide sequence homology to pyrR. Finally, to minimize the mutagenic effect of 5FU, we also improved the pCat-upp/5FU system by reducing the concentration of 5FU from 1 mM to 5 µM using a defined synthetic medium. The optimized system/conditions were used to successfully replace the ldh gene by the sadh-hydG operon to convert acetone into isopropanol.

Exploitation of a Type 1 Toxin-Antitoxin System as an Inducible Counter-Selective Marker for Genome Editing in the Acetogen Eubacterium limosum.

Targeted mutations in the anaerobic methylotroph Eubacterium limosum have previously been obtained using CRISPR-based mutagenesis methods. In this study, a RelB-family toxin from Eubacterium callanderi was placed under the control of an anhydrotetracycline-sensitive promoter, forming an inducible counter-selective system. This inducible system was coupled with a non-replicative integrating mutagenesis vector to create precise gene deletions in Eubacterium limosum B2. The genes targeted in this study were those encoding the histidine biosynthesis gene hisI, the methanol methyltransferase and corrinoid protein mtaA and mtaC, and mtcB, encoding an Mttb-family methyltransferase which has previously been shown to demethylate L-carnitine. A targeted deletion within hisI brought about the expected histidine auxotrophy, and deletions of mtaA and mtaC both abolished autotrophic growth on methanol. Deletion of mtcB was shown to abolish the growth of E. limosum on L-carnitine. After an initial selection step to isolate transformant colonies, only a single induction step was required to obtain mutant colonies for the desired targets. The combination of an inducible counter-selective marker and a non-replicating integrative plasmid allows for quick gene editing of E. limosum.

Covalent attachment of FeFe hydrogenases to carbon electrodes for direct electron transfer.

Direct electron transfer between enzymes and electrodes is now commonly achieved, but obtaining protein films that are very stable may be challenging. This is particularly crucial in the case of hydrogenases, the enzymes that catalyze the biological conversion between dihydrogen and protons, because the instability of the hydrogenase films may prevent the use of these enzymes as electrocatalysts of H(2) oxidation and production in biofuel cells and photoelectrochemical cells. Here we show that two different FeFe hydrogenases (from Chamydomonas reinhardtii and Clostridium acetobutylicum) can be covalently attached to functionalized pyrolytic graphite electrodes using peptidic coupling. In both cases, a surface patch of lysine residues makes it possible to favor an orientation that is efficient for fast, direct electron transfer. High hydrogen-oxidation current densities are maintained for up to one week, the only limitation being the intrinsic stability of the enzyme. We also show that covalent attachment has no effect on the catalytic properties of the enzyme, which means that this strategy can also used be for electrochemical studies of the catalytic mechanism.

Relating diffusion along the substrate tunnel and oxygen sensitivity in hydrogenase.

In hydrogenases and many other redox enzymes, the buried active site is connected to the solvent by a molecular channel whose structure may determine the enzyme's selectivity with respect to substrate and inhibitors. The role of these channels has been addressed using crystallography and molecular dynamics, but kinetic data are scarce. Using protein film voltammetry, we determined and then compared the rates of inhibition by CO and O2 in ten NiFe hydrogenase mutants and two FeFe hydrogenases. We found that the rate of inhibition by CO is a good proxy of the rate of diffusion of O2 toward the active site. Modifying amino acids whose side chains point inside the tunnel can slow this rate by orders of magnitude. We quantitatively define the relations between diffusion, the Michaelis constant for H2 and rates of inhibition, and we demonstrate that certain enzymes are slowly inactivated by O2 because access to the active site is slow.

Genome Sequence of Eubacterium limosum B2 and Evolution for Growth on a Mineral Medium with Methanol and CO2 as Sole Carbon Sources.

Eubacterium limosum is an acetogen that can produce butyrate along with acetate as the main fermentation end-product from methanol, a promising C1 feedstock. Although physiological characterization of E. limosum B2 during methylotrophy was previously performed, the strain was cultured in a semi-defined medium, limiting the scope for further metabolic insights. Here, we sequenced the complete genome of the native strain and performed adaptive laboratory evolution to sustain growth on methanol mineral medium. The evolved population significantly improved its maximal growth rate by 3.45-fold. Furthermore, three clones from the evolved population were isolated on methanol mineral medium without cysteine by the addition of sodium thiosulfate. To identify mutations related to growth improvement, the whole genomes of wild-type E. limosum B2, the 10th, 25th, 50th, and 75th generations, and the three clones were sequenced. We explored the total proteomes of the native and the best evolved clone (n°2) and noticed significant differences in proteins involved in gluconeogenesis, anaplerotic reactions, and sulphate metabolism. Furthermore, a homologous recombination was found in subunit S of the type I restriction-modification system between both strains, changing the structure of the subunit, its sequence recognition and the methylome of the evolved clone. Taken together, the genomic, proteomic and methylomic data suggest a possible epigenetic mechanism of metabolic regulation.

Molecular characterization of the missing electron pathways for butanol synthesis in Clostridium acetobutylicum.

Clostridium acetobutylicum is a promising biocatalyst for the renewable production of n-butanol. Several metabolic strategies have already been developed to increase butanol yields, most often based on carbon pathway redirection. However, it has previously demonstrated that the activities of both ferredoxin-NADP+ reductase and ferredoxin-NAD+ reductase, whose encoding genes remain unknown, are necessary to produce the NADPH and the extra NADH needed for butanol synthesis under solventogenic conditions. Here, we purify, identify and partially characterize the proteins responsible for both activities and demonstrate the involvement of the identified enzymes in butanol synthesis through a reverse genetic approach. We further demonstrate the yield of butanol formation is limited by the level of expression of CA_C0764, the ferredoxin-NADP+ reductase encoding gene and the bcd operon, encoding a ferredoxin-NAD+ reductase. The integration of these enzymes into metabolic engineering strategies introduces opportunities for developing a homobutanologenic C. acetobutylicum strain.

Molecular characterization of the glycerol-oxidative pathway of Clostridium butyricum VPI 1718.

The glycerol oxidative pathway of Clostridium butyricum VPI 1718 plays an important role in glycerol dissimilation. We isolated, sequenced, and characterized the region coding for the glycerol oxidation pathway. Five open reading frames (ORFs) were identified: dhaR, encoding a putative transcriptional regulator; dhaD (1,142 bp), encoding a glycerol dehydrogenase; and dhaK (995 bp), dhaL (629 bp), and dhaM (386 bp), encoding a phosphoenolpyruvate (PEP)-dependent dihydroxyacetone (DHA) kinase enzyme complex. Northern blot analysis demonstrated that the last four genes are transcribed as a 3.2-kb polycistronic operon only in glycerol-metabolizing cultures, indicating that the expression of this operon is regulated at the transcriptional level. The transcriptional start site of the operon was determined by primer extension, and the promoter region was deduced. The glycerol dehydrogenase activity of DhaD and the PEP-dependent DHA kinase activity of DhaKLM were demonstrated by heterologous expression in different Escherichia coli mutants. Based on our complementation experiments, we proposed that the HPr phosphoryl carrier protein and His9 residue of the DhaM subunit are involved in the phosphoryl transfer to dihydroxyacetone-phosphate. DhaR, a potential regulator of this operon, was found to contain conserved transmitter and receiver domains that are characteristic of two-component systems present in the AraC family. To the best of our knowledge, this is the first molecular characterization of a glycerol oxidation pathway in a Gram-positive bacterium.

CO disrupts the reduced H-cluster of FeFe hydrogenase. A combined DFT and protein film voltammetry study.

Carbon monoxide is often described as a competitive inhibitor of FeFe hydrogenases, and it is used for probing H(2) binding to synthetic or in silico models of the active site H-cluster. Yet it does not always behave as a simple inhibitor. Using an original approach which combines accurate electrochemical measurements and theoretical calculations, we elucidate the mechanism by which, under certain conditions, CO binding can cause permanent damage to the H-cluster. Like in the case of oxygen inhibition, the reaction with CO engages the entire H-cluster, rather than only the Fe(2) subsite.