A Prospective Study on the Fermentation Landscape of Gaseous Substrates to Biorenewables Using Methanosarcina acetivorans Metabolic Model

The abundance of methane in shale gas and of other gases such as carbon monoxide, hydrogen, and carbon dioxide as chemical process byproducts has motivated the use of gas fermentation for bioproduction. Recent advances in metabolic engineering and synthetic biology allow for engineering of microbes metabolizing a variety of chemicals including gaseous feeds into a number of biorenewables and transportation liquid fuels. New computational tools enable the systematic exploration of all feasible conversion alternatives. Here we computationally assessed all thermodynamically feasible ways of co-utilizing CH₄, CO, and CO₂ using ferric as terminal electron acceptor for the production of all key precursor metabolites. We identified the thermodynamically feasible co-utilization ratio ranges of CH₄, CO, and CO₂ toward production of the target metabolite(s) as a function of ferric uptake. A revised version of the iMAC868 genome-scale metabolic model of Methanosarcina acetivorans was chosen to assess co-utilization of CH₄, CO, and CO₂ and their conversion into selected target products using the optStoic pathway design tool. This revised version contains the latest information on electron flow mechanisms by the methanogen while supplied with methane as the sole carbon source. The interplay between different gas co-utilization ratios and the energetics of reverse methanogenesis were also analyzed using the same metabolic model.

Introduction

The global increase in oil production, fossil fuel combustion, biomass burning, and hydraulic fracturing of shale gas and climate change concerns has motivated the reduction of emissions from anthropogenic sources. Mitigation of gaseous emissions (such as methane, carbon dioxide, and carbon monoxide) from the environment and their microbial conversion into useful products provides a sustainable and transformative solution that avoids the “food vs. fuel” dilemma. Methane, the major constituent of natural gas, has the highest oxidation potential amongst carbon dioxide, carbon monoxide, and glucose to be converted into a wide range of products including liquid fuels such as ethanol and butanol. Carbon monoxide, often as synthesis gas with varying levels of carbon dioxide and hydrogen (Aasberg-Petersen et al., 2001), along with methane could yield a variable mixture of gases that can be tapped for microbial conversion.

Existing chemical gas-to-liquid (GTL) technologies (i.e., GTL process using the Fischer-Tropsch method) require high operating temperatures and pressures, involve high CapEx costs, yield generally low carbon conversion efficiency, and cannot directly convert methane into the desired bioproducts (Dry, 2002; Steynberg, 2004; Haynes and Gonzalez, 2014). The biological routes of methane utilization, have received renewed interest because of process simplicity (Lopez et al., 2013), selectivity toward targeted pathways (Haynes and Gonzalez, 2014; Mueller et al., 2015), and recent advancements in the characterization and genetic tools of methanotrophic microbes enabling direct transformation of methane into valuable chemicals and fuel molecules (Coleman et al., 2014; Fei et al., 2014; Strong et al., 2015; Henard et al., 2016). Much of the current industrial applications of methane utilization have been devoted to the use of aerobic methanotrophic bacteria (Fei et al., 2014). In contrast, the global methane cycle is primarily controlled by the syntrophy of microorganisms living in anoxic environments. Although biological methane conversion can occur in oxic habitats (Conrad, 2009; Knittel and Boetius, 2009), more than 80% of methane produced in the world's oceans is estimated to be converted anaerobically (Orphan et al., 2001). In addition, anaerobic routes for methane metabolism offer better carbon and energy efficiency compared with aerobic pathways (Mueller et al., 2015; Nazem-Bokaee et al., 2016). Difficulties in culturing anaerobic methanotrophs in the lab, arising from syntrophy requirements, have hampered their rapid characterization and application. Nonetheless, recent observation of methane utilization by anaerobic methanotrophic archaea (ANME) decoupled from their sulfate-reducing bacteria (SRB) partners in the presence of artificial electron acceptors (Scheller et al., 2016) revealed new avenues for direct anaerobic conversion of methane by ANMEs into useful chemicals. So far there is no microbe capable of AOM utilizing other gaseous substrates at industrial scale. Acetogens has been the workhouse for gas fermentation in industry for over two decades. Anaerobic conversion of carbon monoxide into valuable products such as ethanol, acetate, and 2,3-butanediol at industrial scale has been pursued using different strains of Clostridium (Simpson et al., 2010; Köpke et al., 2011a,b; Tran and Simpson, 2015; Daniell et al., 2016; Martin et al., 2016). A recent study on the co-utilization of carbon dioxide and carbon monoxide or hydrogen to produce acetate using Moorella thermoacetica (Hu et al., 2016) further demonstrates the need for systematic study of co-utilization of various C₁ gases in other potential microbial hosts.

In this work, we aim at developing a computational framework allowing for designing overall thermodynamically feasible conversions of mixes of gaseous molecules into selected metabolites and, then, investigating the metabolic capabilities of a selected microorganism in response to introducing new gas mixture combinations. Using optStoic (Chowdhury and Maranas, 2015) we exhaustively identified all thermodynamically feasible optimal conversion stoichiometries making use of a combination of CH₄, CO, and CO₂. Note that there exist many other computational tools for pathway design (Hadadi and Hatzimanikatis, 2015; Long et al., 2015; Nazem-Bokaee and Senger, 2015; Huang et al., 2017). Ten key branch point (precursor) metabolites (Noor et al., 2010) were selected owing to their essentiality for anabolic processes found in all forms of life as well as their crucial role as building blocks for producing many commodity and specialty chemicals listed as top value-added chemicals by the U.S. Department Of Energy (DOE). Maximum uptake of carbon coming from CH₄, CO, or CO₂ and their co-utilization ratios have been assessed as well as the indispensability of ferric ion as an electron acceptor. To analyze metabolic pathway usage at different co-utilization ratios of CH₄, CO, and CO₂ designed by optStoic algorithm, a revised version of the iMAC868 genome-scale metabolic model of the methanogenic archaeon Methanosarcina acetivorans (Nazem-Bokaee et al., 2016) was used allowing for full tracking of carbon and electron flow within the reversal of methanogenesis pathway. Recent studies identified the existence of an electron bifurcating multi-complex enzyme, cytosolic heterodisulfide reductase HdrABC, shedding light into pathways for utilizing methane by M. acetivorans in the presence of ferric to produce useful chemicals such as acetate (Yan et al., 2017; Nazem-Bokaee et al., 2018). It has been shown before that M. acetivorans is capable of growing with carbon monoxide (Rother and Metcalf, 2004; Lessner et al., 2006) and metabolizing carbon dioxide (in the form of bicarbonate) along with methane (Soo et al., 2016), thus, making the archaeon a suitable platform to study the conversion of varying mixtures of these gases into useful products. The computational framework put forth in this study could inform design of novel metabolic engineering strategies for the industrial production of bio-based chemicals and liquid fuels from mixed gaseous feeds.

Methods

Computational Design of Overall Stoichiometries for Gas Co-utilization

To explore optimal overall stoichiometries for conversion of gaseous molecules (i.e., CH₄, CO, and CO₂) into target products, the optStoic procedure (Chowdhury and Maranas, 2015) was implemented in Python so that it can be freely accessible (Supplementary Data Sheet 2). The goal was to design overall stoichiometries informing thermodynamically feasible co-utilization of the gaseous molecules leading to the production of 10 C-mol of products listed in Table 1 (equation 1).

In the postulated overall stoichiometry s₁, s₂, and s₃ are the optimal coefficients of methane, carbon dioxide, and carbon monoxide, respectively. Because the target products listed in Table 1 contain varying number of carbons, fixing the stoichiometry of target product in Equation 1 enables a direct comparison of gaseous feed ratios on a per carbon mol basis. In the optStoic algorithm, water molecules and protons can be taken up or produced as needed so that Equation 1 remains elementally and charge balanced. Furthermore, phosphate, ammonia, and hydrogen sulfide were added to balance Equation 1 when a target product contains phosphorous, nitrogen, and sulfur, respectively. No carbon-containing compound other than methane, carbon dioxide, and carbon monoxide was allowed as an additional substrate. The choice of the products listed in Table 1 is based on their essentiality in the metabolism of almost all forms of life (Noor et al., 2010) and their significance in being used as building blocks of many commodity and specialty chemicals as mentioned in the DOE list of top value-added chemicals. The performance criteria of the overall conversion shown in Equation 1 were to maximize s₁, s₂, or s₃ separately at a specified ferric uptake. To safeguard the thermodynamic feasibility of all conversions, the minimum overall standard ΔG was set to be less than zero. A previously assembled database of metabolites (Chowdhury and Maranas, 2015) was used to explore the optimal combination of reactants and products for any given overall stoichiometry. COBRApy (Ebrahim et al., 2013) with built-in cGLPK (http://www.gnu.org/software/glpk/) solver was used to solve the optimization problems written in Python 2.7.

TABLE 1

Table 1. The key branch point (precursor) metabolites essential for anabolic processes found in all forms of life considered as target products of gaseous fermentation in this study.

Modifications to the iMAC868 Metabolic Model of M. acetivorans

Since the development and release of the iMAC868 metabolic model (Nazem-Bokaee et al., 2016), there have been new experimental studies aimed at better understanding the electron flow mechanisms and biochemistry of M. acetivorans growing on methane (Yan et al., 2017, 2018; Nazem-Bokaee et al., 2018). This provided the impetus for updating the iMAC868 model of this methanogen to catalog these findings. It was recently shown that M. acetivorans expresses a multi-unit cytosolic heterodisulfide reductase complex, HdrA2B2C2, when grown with methane in the presence of ferric (Yan et al., 2017; Yan and Ferry, 2018) that can partition electrons (i.e., bifurcate electrons) coming from cofactor F₄₂₀ (reduced) between ferredoxin (with lower electrode potential) and heterodisulfide (with higher electrode potential) (Figure 1). Therefore, HdrA2B2C2 complex bypasses thermodynamic uphill for direct electron transfer from cofactor F₄₂₀ to ferredoxin during ferric-dependent methanotrophy by M. acetivorans. This important finding introduces a new metabolic capability of M. acetivorans and, therefore, was cataloged in the updated version of the iMAC868 metabolic model. The resulting coenzyme M and coenzyme B are re-used to regenerate heterodisulfide used for activation of methane. The reduced ferredoxin is used to drive the biosynthesis of acetyl-CoA by CO dehydrogenase, Cdh. Therefore, we replaced the previously used electron flow routes in our model with the new route representing the newly elucidated function of HdrA2B2C2 (see Figure 1). We found that the model accommodated the new electron bifurcation mechanism providing new insights about the key role of ferric in the distribution of electrons between major products of methanotrophy as well as on energy conservation mechanisms (Nazem-Bokaee et al., 2018). The model was assembled in a format compatible for flux balance analysis (Orth et al., 2010). FBA optimization problems were solved by GNU Linear Programming Kit (GLPK) (http://www.gnu.org/software/glpk/) solver in Matlab using COBRA toolbox (Schellenberger et al., 2011). Flux variability analysis (FVA) was performed to obtain range of fluxes under optimal growth conditions as described previously (Mahadevan and Schilling, 2003). Both FBA and FVA problems incorporated overall thermodynamic feasibility constraints (overall ΔG≤0).

FIGURE 1

Figure 1. Electron bifurcation mechanism by HdrA2B2C2 complex of M. acetivorans in the presence of external electron acceptor when grown with methane (see Yan et al., 2017 for more details). F₄₂₀: Cofactor F₄₂₀; F₄₂₀H₂: reduced form of cofactor F₄₂₀; Fdx: ferredoxin; Fdx²⁻: reduced form of ferredoxin; HSCoM: coenzyme M; HSCoB: coenzyme B; CoM-S-S-CoB: heterodisulfide.

Results and Discussion

Thermodynamically Feasible Gas Co-utilization Stoichiometries Designed by Optstoic

The thermodynamically feasible ranges of co-utilization of CH₄, CO, and CO₂ for the production of target chemicals listed in Table 1 were predicted by optStoic to be dependent on the level of available ferric. Figure 2 shows this dependency for three selected products with varying degrees of reduction. As the ferric level goes up (i.e., increasing the electron sink capacity), methane usage increases in proportion. However, only some ratio ranges of the CH₄-CO-CO₂ triplet lead to thermodynamically feasible production of the target molecules (Figure 2). Here, increasing ferric levels provides opportunity for CO₂ utilization levels to go up by accepting electrons coming from methane. This increase, however, is at the expense of reduction in CO utilization levels to satisfy stoichiometric and thermodynamics feasibility of the overall conversion.

FIGURE 2

Figure 2. Ternary diagrams showing the contribution of gaseous carbon sources (i.e., CH₄, CO, or CO₂) in the production of 10 C-mol oxaloacetate (A), glyceraldehyde-3-phosphate (B), or acetyl-CoA (C) as selected target products. Colorful symbols on the bottom right of the figure show the range of ferric (Fe³⁺) uptake (in moles) at which the overall gases-to-product conversion shown in Equation 1 is thermodynamically feasible. Each symbol on the ternary plots represents a single independent thermodynamically feasible stoichiometric conversion of gases-to-product simulated by optStoic algorithm. In each simulation, the stoichiometries of the target product, ferric, and one of the gases are fixed and the objective is to maximize the stoichiometries of the other two gases. The moles of CH₄, CO, or CO₂ in the overall stoichiometry are normalized to be between zero and one in the ternary diagram.

The minimum and maximum moles of ferric required to maintain any thermodynamically feasible gas co-utilization are given in Table 2 for all target products listed in Table 1. For example, a minimum of 3.04 mol ferric was required to obtain any feasible conversion of gaseous substrates toward glyceraldehyde-3-phosphate (GAP) while no feasible overall stoichiometry was found with a methane carbon contribution <38%. The overall conversions given in Table 2 also unmask the possibility of designing gas bi-utilization (where either CH₄ and CO or CH₄ and CO₂ can be co-utilized) at varying levels of ferric, which is further explained in the following sections.

TABLE 2

Table 2. optStoic-predicted overall stoichiometric conversions (middle column) for which the stoichiometry of CH₄, CO₂, or CO were maximized independently.

The maximum amount of carbon that can be incorporated to target products from CH₄, CO, or CO₂, depends on the target molecule C/O ratio and reduction level. Figure 3 displays how the choice of target molecule (those listed in Table 1) affects the maximum carbon contributed by the three gaseous feeds. For example, CO could be the top supplier of carbon for oxaloacetate (OXA), as expected, because OXA is highly oxidized. Note that Figure 3 does not directly represent maximum co-utilization of ratios of gases; however, it demarcates the theoretical limits on utilizing any of the gases for the production of each target product. For example, under the defined criteria for optStoic, it would be thermodynamically infeasible to design an overall stoichiometry for pyruvate production in which carbon coming from methane co-utilized with other gases exceeded 50% (see also Table 2 for all stoichiometric designs).

FIGURE 3

Figure 3. optStoic-predicted maximum carbon (shown as percentage on the Y-axis) contribution from carbon monoxide (blue), methane (green), or carbon dioxide (red) for the production of different target products (For abbreviations see Table 1). These maxima are from different independent overall stoichiometry designs predicted by optStoic (see Table 2 for all stoichiometries and performance criteria).

Nonetheless, methane contributes the most carbon at maximum ferric uptake levels. In addition, imposing a more negative requirement for the overall standard free energy of change results in less carbon contributed from CO₂ (since it has the lowest Gibbs free energy of formation among CO and CH₄) leading to a decline in maximum co-utilization ratios of CO₂-to-CH₄. Such information could be useful in designing and/or modifying bioconversions based on varying compositions of industrial gas waste streams (Subramani and Gangwal, 2008; Lackey et al., 2015).

In the following section we describe how overall stoichiometry designs generated by optStoic could be used to inform metabolic engineering strategies through using the updated iMAC868 metabolic model of M. acetivorans as a platform.

Metabolic Capabilities of M. acetivorans During Gas Co-utilization

The ratio of industrial waste gases is often highly variable from stream to stream leading to difficulties in predicting desirable gas stream-to-target product conversions (Williams et al., 2007; Subramani and Gangwal, 2008). The optStoic designs could serve as a guide to estimate feasible conversions using the metabolic model of M. activorans. We selected oxaloacetate (OXA), glyceraldehyde 3-phosphate (GAP), and acetyl-CoA (ACA) (out of hundreds of unique overall stoichiometry designs) based on their distinct differences as shown in Figures 2, 3 (also in Table 1) as well as their importance as building blocks of numerous valuable end products. We chose the stoichiometric ratios of CH₄, CO, and CO₂ at an arbitrary ferric level of 4 moles at which co-utilization of the three gases for the production of OXA, ACA, and GAP was predicted by optStoic to be thermodynamically feasible (Table 3). To implement these stoichiometric ratios in the context of the metabolic model of M. acetivorans, the lower and upper bounds of the reactions corresponding to the uptake of CH₄, CO, and CO₂ in the iMAC868 metabolic model were fixed to the stoichiometric ratios of CH₄, CO, and CO₂ shown in Table 3. Analysis of the flux distribution through the metabolic network confirmed the usage of the reversal of the methanogenesis pathway indicating the incorporation of the gaseous substrates into biomass and cofactor biosynthesis. The iMAC868 metabolic model also predicts the uptake of ammonia, hydrogen sulfide, and phosphate as essential sources of nitrogen, sulfur, and phosphorus, respectively, consistent with the overall optStoic design. Flux variability analysis results in predicting a maximum yield of 2.499 (mol per 10 C-mol of gases) for OXA. This is in agreement with a 2.5 stoichiometric value predicted by optStoic leading to the same ratio of CH₄, CO, and CO₂ co-utilization implying that metabolism remains unaffected even at maximum OXA production yield. The maximum yields of GAP and ACA predicted by the iMAC868 metabolic model are 2.944 and 0.388 (mol per 10 C-mol of gases), respectively, which is only 11.5, and 11% less than the optimal overall stoichiometries obtained by optStoic. This difference is due to the inclusion of many more cofactors and intermediate metabolites in the metabolic network compared to the consideration of one simple overall stoichiometry as that shown in Equation 1. Further analysis of the flux through the formation of biomass, as another product of the metabolic network, reveals a maximum biomass yield of 0.217 at a ferric level of 4.2 (mol per 10 C-mol of gases) when using the gas ratios optimized for ACA production predicted by optStoic (see Table 3). Therefore, the optStoic design could quickly inform potential gas co-utilization ratios at which a certain level of cellular growth can be achieved. It should be noted that there exist other possible gas co-utilization ratios that could end up obtaining similar biomass yields. For example, using the gas ratios optimized for GAP production (see Table 3) results in achieving a maximum biomass yield of 0.224 at a ferric level of 4.5 (mol per 10 C-mol of gases), which is only 3% higher than what could be achieved at a gas composition optimized for ACA production and is slightly richer in CO (see Table 3).

TABLE 3

Table 3. optStoic-designed stoichiometries (mol) of methane (s_CH4), carbon monoxide (s_CO), and carbon dioxide (s_CO2) resulted in the production of 10 C-mol of three selected target products used to constrain the in silico uptake of these gases by the iMAC868 metabolic model of M. acetivorans.

It has been postulated that M. acetivorans reduces ferric at multi-heme c-type cytochromes sites to which electrons are shuttled by membrane-bound methanophenazine (Yan et al., 2017). Depending on the composition of gaseous substrates being used (given in Table 3), the iMAC868 metabolic model predicts that at least 13% (up to 20%) of heterodisulfide has to be reduced through the membrane-bound heterodisulfide reductase (HdrDE) that reduces methanophenazine. The remaining heterodisulfide can be reduced via either the cytosolic HdrA2B2C2 or HdrDE. Reduced cofactor F₄₂₀, which donates electrons to ferredoxin and heterodisulfide at the HdrA2 site, can be regenerated through any of F₄₂₀ dehydrogenase (Fpo), F₄₂₀-dependent methylene-H₄MPT reductase (Mer), F₄₂₀-dependent methylene-H₄MPT dehydrogenase (Mtd), or F₄₂₀-dependent NADP reductase enzyme complexes according to the metabolic model predictions.

To further explore the metabolic capabilities of M. acetivorans, we decided to analyze the theoretical limits of ethanol and butanol co-production during CH₄ and CO co-utilization by the iMAC868 model. The biological co-production of alcohols has been reported in the literature where acetone/butanol/ethanol (ABE) fermentation process by clostridial strains has been studied the most and implemented industrially (Worden et al., 1991; Lee et al., 2008; Tracy et al., 2012; Carlson and Papoutsakis, 2017; Fernandez-Naveira et al., 2017). However, most traditional ABE process suffers from high feedstock costs (Green, 2011) and, thus, the use of cheap sources such as C₁ gas substrates suggests a promising alternate route (Dürre, 2017). Nonetheless, the current C₁ gas fermentation technology is mainly relied on making use of acetogens (De Tissera et al., 2017). Here, the motivation was to study the co-production of alcohols from co-utilization of C₁ gaseous substrates in non-traditional hosts such as M. acetivorans. For that, first, optStoic was used to design overall conversions such that CH₄ and CO co-utilization (using one mole of ferric as basis) results in production of one mole butanol while maximizing the production of ethanol and several other alcohol molecules as co-products (Figure 4). optStoic was also applied to examine how conversion of CH₄ and CO to butanol would vary for different electron acceptors other than ferric. Almost all electron acceptors examined allowed for the same ratio of CH₄ and CO co-utilization except for thrithionate that enabled about three times higher co-utilization ratio (Supplementary Figure S1 in Data Sheet 1). However, ethanol production as a co-product of butanol production when using trithionate/bisulfate as the electron acceptor pair was only 0.3% of that achievable by using ferric/ferrous as electron acceptor pair. Thus, the overall conversion design using ferric as electron acceptor was employed for analyzing metabolic capabilities of M. acetivorans for co-production of ethanol and butanol. The original version of the iMAC868 metabolic model comprises the biosynthetic pathways for ethanol and butanol production (Nazem-Bokaee et al., 2016).

FIGURE 4

Figure 4. optStoic-predicted co-production of selected alcohols (with their number of carbons given in parenthesis) along with butanol in the presence of ferric as electron acceptor. Y-axis indicates that under the design criteria of optStoic, where the only products of CO and CH₄ co-utilization are butanol and one of the shown alcohols, how much (in percent) of the total product could be each alcohol molecule (gray area of the bars show percent butanol of the total). Italic numbers on top of the bars show CO to CH₄ gas co-utilization ratios.

By constraining the lower and upper bounds of the reaction corresponding to the exchange of butanol in the metabolic model to one, and fixing the bounds of reactions corresponding to the uptake of CH₄ and CO to the respective ratio given in Figure 4 (i.e., 0.47), the model predicts that a maximum of 3.779 moles of ethanol per mole of butanol could be produced (Figure 5). The ethanol-to-butanol molar ratio predicted by optStoic at the same gas co-utilization ratio was 3.682, which is only 2.5% different from that predicted by the iMAC868 metabolic model.

FIGURE 5

Figure 5. Predictive capabilities of iMAC868 metabolic model of M. acetivorans during CO and CH₄ co-utilization in the presence of ferric for butanol and ethanol co-production. Top panel: prediction of ethanol and butanol co-production feasibility over a range of ferric reduction levels. Bottom panel: partitioning of methyl-tetrahydrosarcinapterin (CH₃-H₄SPT) flux (denoted as v) between CO₂ pathway (Mer) and acetyl-CoA biosynthesis pathway (Cdh) during reversal of the methanogenesis pathway by M. acetivorans.

Ethanol co-production with butanol was predicted by iMAC868 metabolic model to be feasible over a range of ferric reduction values from 0.68 up to 2.2 (Figure 5). However, ethanol co-production decreases as ferric reduction levels increases because the reducing power for generating acetyl-CoA, the precursor for both ethanol and butanol production, diminishes. The bottom panel of Figure 5 shows that increasing ferric reduction capacity results in re-routing more methane (through Mer) toward the methyotrophic pathway. Thus, acetyl-CoA production via Cdh remains at stoichiometric limits necessary for satisfying fixed amount of butanol production. Nonetheless, the flux through Cdh could never become zero and at ferric levels of 2.2 mol/mol of gases at least 20% of the CH₃-H₄SPT has to be converted to acetyl-CoA to maintain cellular growth. This analysis demonstrate the usefulness of computational tools such as optStoic in guiding metabolic engineering design/analysis for a given bioconversion.

Summary and Conclusion

In this work, we have demonstrated the utility of deploying computational tools such as optStoic along with “genome-scale” metabolic modeling to inform optimal metabolic engineering designs and strategies satisfying an overall desired bioconversion. The optStoic formulation allowed for the exploration of all overall conversions rooting from the co-utilization of low-value C₁ gaseous feedstocks (i.e., CH₄, CO₂, and CO) ending up in the production of precursors used for making high-value biorenewables. We targeted ten key branch point metabolites that have been used extensively as building blocks for the production of many commodity and specialty chemicals such as acetate, terpenoids, and synthetic sugars among others. We showed that the proper choice of an electron acceptor (i.e., ferric) could bypass the thermodynamic barriers for electron flow in the gas-to-chemicals conversions. Furthermore, we showed that there exist well defined gas co-utilization ranges, which are feasible at varying levels of ferric, dependent on the choice of target product. Maximum ferric usage as well as maximum carbon contribution from each of the CH₄, CO₂, and CO was analyzed that could lead to new or improved gas co-utilization designs. Using optStoic designs as a guide, metabolic capacities of M. acetivorans as the model host was examined owing to its diverse substrate utilization abilities and the progress in its genetic engineering tools. Equipped with latest electron flow mechanisms during growth with methane, the iMAC868 metabolic model of M. acetivorans provided information on the partitioning of electrons within the methanogenesis reversal pathway as well as on distribution of carbons coming from co-utilization of mixtures of gases toward selected products. The combined use of optStoic and metabolic modeling presented in this work puts forth an efficient platform for quickly exploring in silico the feasibility and limits of various gaseous substrate utilization options.

Author Contributions

HN-B wrote computer scripts, performed the simulations and analyses, designed and generated the figures and tables, and wrote the manuscript. CM supervised and contributed to the design of the study, wrote the manuscript, and critically revised the manuscript. Both authors read and approved the final manuscript.

Funding

This work was supported by funding from The Advanced Research Projects Agnecy-Energy (ARPA-E), US Department of Energy to CM, Drs. Thomas K. Wood, and James G. Ferry, Grant DE-AR0000431.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.01855/full#supplementary-material

References

Aasberg-Petersen, K., Hansen, J. H. B., Christensen, T. S., Dybkjaer, I., Christensen, P. S., Nielsen, C. S., et al. (2001). Technologies for large-scale gas conversion. Appl. Catal. Gen. 221, 379–387. doi: 10.1016/S0926-860X(01)00811-0

CrossRef Full Text | Google Scholar

Carlson, E. D., and Papoutsakis, E. T. (2017). Heterologous expression of the clostridium carboxidivorans CO dehydrogenase alone or together with the acetyl coenzyme a synthase enables both reduction of CO2 and oxidation of CO by clostridium acetobutylicum. Appl. Environ. Microbiol. 83:e00829-17. doi: 10.1128/AEM.00829-17

PubMed Abstract | CrossRef Full Text

Chowdhury, A., and Maranas, C. D. (2015). Designing overall stoichiometric conversions and intervening metabolic reactions. Sci. Rep. 5:16009. doi: 10.1038/srep16009

PubMed Abstract | CrossRef Full Text | Google Scholar

Coleman, W. J., Vidanes, G. M., Cottarel, G., Muley, S., Kamimura, R., Javan, A. F., et al. (2014). Biological Conversion of Multi-Carbon Compounds from Methane. US Patent US20140273128 A1.

Conrad, R. (2009). The global methane cycle: recent advances in understanding the microbial processes involved. Environ. Microbiol. Rep. 1, 285–292. doi: 10.1111/j.1758-2229.2009.00038.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Daniell, J., Nagaraju, S., Burton, F., Kopke, M., and Simpson, S. D. (2016). Low-carbon fuel and chemical production by anaerobic gas fermentation. Adv. Biochem. Eng. Biotechnol. 156, 293–321. doi: 10.1007/10_2015_5005

PubMed Abstract | CrossRef Full Text | Google Scholar

De Tissera, S., Kopke, M., Simpson, S. D., Humphreys, C., Minton, N. P., and Durre, P. (2017). Syngas biorefinery and syngas utilization. Adv. Biochem. Eng. Biotechnol. doi: 10.1007/10_2017_5. [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

Dry, M. E. (2002). The Fischer-Tropsch process: 1950-2000. Catal Today 71, 227–241. doi: 10.1016/S0920-5861(01)00453-9

CrossRef Full Text | Google Scholar

Dürre, P. (2017). Gas fermentation - a biotechnological solution for today's challenges. Microb. Biotechnol. 10, 14–16. doi: 10.1111/1751-7915.12431

PubMed Abstract | CrossRef Full Text | Google Scholar

Ebrahim, A., Lerman, J. A., Palsson, B. O., and Hyduke, D. R. (2013). COBRApy: COnstraints-based reconstruction and analysis for python. BMC Syst. Biol. 7:74. doi: 10.1186/1752-0509-7-74

PubMed Abstract | CrossRef Full Text | Google Scholar

Fei, Q., Guarnieri, M. T., Tao, L., Laurens, L. M., Dowe, N., and Pienkos, P. T. (2014). Bioconversion of natural gas to liquid fuel: opportunities and challenges. Biotechnol. Adv. 32, 596–614. doi: 10.1016/j.biotechadv.2014.03.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Fernandez-Naveira, A., Veiga, M. C., and Kennes, C. (2017). H-B-E (hexanol-butanol-ethanol) fermentation for the production of higher alcohols from syngas/waste gas. J. Chem. Technol. Biotechnol. 92, 712–731. doi: 10.1002/jctb.5194

CrossRef Full Text | Google Scholar

Green, E. M. (2011). Fermentative production of butanol–the industrial perspective. Curr. Opin. Biotechnol. 22, 337–343. doi: 10.1016/j.copbio.2011.02.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Hadadi, N., and Hatzimanikatis, V. (2015). Design of computational retrobiosynthesis tools for the design of de novo synthetic pathways. Curr. Opin. Chem. Biol. 28, 99–104. doi: 10.1016/j.cbpa.2015.06.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Haynes, C. A., and Gonzalez, R. (2014). Rethinking biological activation of methane and conversion to liquid fuels. Nat. Chem. Biol. 10, 331–339. doi: 10.1038/nchembio.1509

PubMed Abstract | CrossRef Full Text | Google Scholar

Henard, C. A., Smith, H., Dowe, N., Kalyuzhnaya, M. G., Pienkos, P. T., and Guarnieri, M. T. (2016). Bioconversion of methane to lactate by an obligate methanotrophic bacterium. Sci. Rep. 6:21585. doi: 10.1038/srep21585

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, P., Chakraborty, S., Kumar, A., Woolston, B., Liu, H., Emerson, D., et al. (2016). Integrated bioprocess for conversion of gaseous substrates to liquids. Proc. Natl. Acad. Sci. U.S.A. 113, 3773–3778. doi: 10.1073/pnas.1516867113

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, Y., Zhong, C., Lin, H. X., and Wang, J. (2017). A method for finding metabolic pathways using atomic group tracking. PLoS ONE 12:e0168725. doi: 10.1371/journal.pone.0168725

PubMed Abstract | CrossRef Full Text | Google Scholar

Knittel, K., and Boetius, A. (2009). Anaerobic oxidation of methane: progress with an unknown process. Annu. Rev. Microbiol. 63, 311–334. doi: 10.1146/annurev.micro.61.080706.093130

PubMed Abstract | CrossRef Full Text | Google Scholar

Köpke, M., Mihalcea, C., Bromley, J. C., and Simpson, S. D. (2011a). Fermentative production of ethanol from carbon monoxide. Curr. Opin. Biotechnol. 22, 320–325. doi: 10.1016/j.copbio.2011.01.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Köpke, M., Mihalcea, C., Liew, F., Tizard, J. H., Ali, M. S., Conolly, J. J., et al. (2011b). 2,3-butanediol production by acetogenic bacteria, an alternative route to chemical synthesis, using industrial waste gas. Appl. Environ. Microbiol. 77, 5467–5475. doi: 10.1128/AEM.00355-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Lackey, J. C., Peppley, B., Champagne, P., and Maier, A. (2015). Composition and uses of anaerobic digestion derived biogas from wastewater treatment facilities in North America. Waste Manag. Res. 33, 767–771. doi: 10.1177/0734242X15589781

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, S. Y., Park, J. H., Jang, S. H., Nielsen, L. K., Kim, J., and Jung, K. S. (2008). Fermentative butanol production by Clostridia. Biotechnol. Bioeng. 101, 209–228. doi: 10.1002/bit.22003

PubMed Abstract | CrossRef Full Text | Google Scholar

Lessner, D. J., Li, L., Li, Q., Rejtar, T., Andreev, V. P., Reichlen, M., et al. (2006). An unconventional pathway for reduction of CO2 to methane in CO-grown Methanosarcina acetivorans revealed by proteomics. Proc. Natl. Acad. Sci. U.S.A. 103, 17921–17926. doi: 10.1073/pnas.0608833103

PubMed Abstract | CrossRef Full Text | Google Scholar

Long, M. R., Ong, W. K., and Reed, J. L. (2015). Computational methods in metabolic engineering for strain design. Curr. Opin. Biotechnol. 34, 135–141. doi: 10.1016/j.copbio.2014.12.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Lopez, J. C., Quijano, G., Souza, T. S., Estrada, J. M., Lebrero, R., and Munoz, R. (2013). Biotechnologies for greenhouse gases (CH(4), N(2)O, and CO(2)) abatement: state of the art and challenges. Appl. Microbiol. Biotechnol. 97, 2277–2303. doi: 10.1007/s00253-013-4734-z

CrossRef Full Text | Google Scholar

Mahadevan, R., and Schilling, C. H. (2003). The effects of alternate optimal solutions in constraint-based genome-scale metabolic models. Metab. Eng. 5, 264–276. doi: 10.1016/j.ymben.2003.09.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Martin, M. E., Richter, H., Saha, S., and Angenent, L. T. (2016). Traits of selected Clostridium strains for syngas fermentation to ethanol. Biotechnol. Bioeng. 113, 531–539. doi: 10.1002/bit.25827

PubMed Abstract | CrossRef Full Text | Google Scholar

Mueller, T. J., Grisewood, M. J., Nazem-Bokaee, H., Gopalakrishnan, S., Ferry, J. G., Wood, T. K., et al. (2015). Methane oxidation by anaerobic archaea for conversion to liquid fuels. J. Ind. Microbiol. Biotechnol. 42, 391–401. doi: 10.1007/s10295-014-1548-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Nazem-Bokaee, H., Gopalakrishnan, S., Ferry, J. G., Wood, T. K., and Maranas, C. D. (2016). Assessing methanotrophy and carbon fixation for biofuel production by Methanosarcina acetivorans. Microbial Cell Factories 15, 1–13. doi: 10.1186/s12934-015-0404-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Nazem-Bokaee, H., and Senger, R. S. (2015). ToMI-FBA: A genome-scale metabolic flux based algorithm to select optimum hosts and media formulations for expressing pathways of interest. AIMS Bioeng. 2, 335–374. doi: 10.3934/bioeng.2015.4.335

CrossRef Full Text | Google Scholar

Nazem-Bokaee, H., Yan, Z., Maranas, C. D., and Ferry, J. G. (2018). “The biochemistry and physiology of respiratory-driven reversed methanogenesis,” in Methane Biocatalysis: Paving the Way to Sustainability, eds M. G. Kalyuzhnaya and X. H. Xing, (Cham: Springer International Publishing), 183–197.

Google Scholar

Noor, E., Eden, E., Milo, R., and Alon, U. (2010). Central carbon metabolism as a minimal biochemical walk between precursors for biomass and energy. Mol. Cell 39, 809–820. doi: 10.1016/j.molcel.2010.08.031

PubMed Abstract | CrossRef Full Text | Google Scholar

Orphan, V. J., House, C. H., Hinrichs, K. U., McKeegan, K. D., and DeLong, E. F. (2001). Methane-consuming archaea revealed by directly coupled isotopic and phylogenetic analysis. Science 293, 484–487. doi: 10.1126/science.1061338

PubMed Abstract | CrossRef Full Text | Google Scholar

Orth, J. D., Thiele, I., and Palsson, B. O. (2010). What is flux balance analysis? Nat. Biotechnol. 28, 245–248. doi: 10.1038/nbt.1614

PubMed Abstract | CrossRef Full Text | Google Scholar

Rother, M., and Metcalf, W. W. (2004). Anaerobic growth of Methanosarcina acetivorans C2A on carbon monoxide: an unusual way of life for a methanogenic archaeon. Proc. Natl. Acad. Sci. U.S.A. 101, 16929–16934. doi: 10.1073/pnas.0407486101

PubMed Abstract | CrossRef Full Text | Google Scholar

Schellenberger, J., Que, R., Fleming, R. M., Thiele, I., Orth, J. D., Feist, A. M., et al. (2011). Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox v2.0. Nat. Protoc. 6, 1290–1307. doi: 10.1038/nprot.2011.308

PubMed Abstract | CrossRef Full Text | Google Scholar

Scheller, S., Yu, H., Chadwick, G. L., McGlynn, S. E., and Orphan, V. J. (2016). Artificial electron acceptors decouple archaeal methane oxidation from sulfate reduction. Science 351, 703–707. doi: 10.1126/science.aad7154

PubMed Abstract | CrossRef Full Text | Google Scholar

Simpson, S. D., Collet, C., Forster, R. L. S., Cockrem, M. C. M., Oakley, S. D., and Kopke, M. (2010). Carbon Capture in Fermentation. US Patent US20100323417 A1.

Soo, V. W., McAnulty, M. J., Tripathi, A., Zhu, F., Zhang, L., Hatzakis, E., et al. (2016). Reversing methanogenesis to capture methane for liquid biofuel precursors. Microb. Cell Fact. 15:11. doi: 10.1186/s12934-015-0397-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Steynberg, A. P. (2004). Introduction to Fischer-Tropsch technology. Fischer-Tropsch Technol. 152, 1–63. doi: 10.1016/S0167-2991(04)80458-0

CrossRef Full Text | Google Scholar

Strong, P. J., Xie, S., and Clarke, W. P. (2015). Methane as a resource: can the methanotrophs add value? Environ. Sci. Technol. 49, 4001–4018. doi: 10.1021/es504242n

PubMed Abstract | CrossRef Full Text | Google Scholar

Subramani, V., and Gangwal, S. K. (2008). A review of recent literature to search for an efficient catalytic process for the conversion of syngas to ethanol. Energy Fuels 22, 814–839. doi: 10.1021/ef700411x

CrossRef Full Text | Google Scholar

Tracy, B. P., Jones, S. W., Fast, A. G., Indurthi, D. C., and Papoutsakis, E. T. (2012). Clostridia: the importance of their exceptional substrate and metabolite diversity for biofuel and biorefinery applications. Curr. Opin. Biotechnol. 23, 364–381. doi: 10.1016/j.copbio.2011.10.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Tran, L. P., and Simpson, S. D. (2015). Fermentation Process. US Patent US20130065282 A1.

Williams, T. C., Shaddix, C. R., and Schefer, R. W. (2007). Effect of syngas composition and CO2-diluted oxygen on performance of a premixed swirl-stabilized combustor. Combust. Sci. Technol. 180, 64–88. doi: 10.1080/00102200701487061

CrossRef Full Text | Google Scholar

Worden, R. M., Grethlein, A. J., Jain, M. K., and Datta, R. (1991). Production of butanol and ethanol from synthesis gas via fermentation. Fuel 70, 615–619. doi: 10.1016/0016-2361(91)90175-A

CrossRef Full Text | Google Scholar

Yan, Z., and Ferry, J. G. (2018). Electron Bifurcation and Confurcation in Methanogenesis and Reverse Methanogenesis. Front Microbiol. 9:1322. doi: 10.3389/fmicb.2018.01322

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, Z., Joshi, P., Gorski, C. A., and Ferry, J. G. (2018). A biochemical framework for anaerobic oxidation of methane driven by Fe(III)-dependent respiration. Nat. Commun. 9:1642. doi: 10.1038/s41467-018-04097-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, Z., Wang, M., and Ferry, J. G. (2017). A Ferredoxin- and F420H2-dependent, electron-bifurcating, heterodisulfide reductase with homologs in the domains bacteria and archaea. MBio 8:e02285-16. doi: 10.1128/mBio.02285-16

PubMed Abstract | CrossRef Full Text | Google Scholar