Skip to main content

Assessing methanotrophy and carbon fixation for biofuel production by Methanosarcina acetivorans



Methanosarcina acetivorans is a model archaeon with renewed interest due to its unique reversible methane production pathways. However, the mechanism and relevant pathways implicated in (co)utilizing novel carbon substrates in this organism are still not fully understood. This paper provides a comprehensive inventory of thermodynamically feasible routes for anaerobic methane oxidation, co-reactant utilization, and maximum carbon yields of major biofuel candidates by M. acetivorans.


Here, an updated genome-scale metabolic model of M. acetivorans is introduced (iMAC868 containing 868 genes, 845 reactions, and 718 metabolites) by integrating information from two previously reconstructed metabolic models (i.e., iVS941 and iMB745), modifying 17 reactions, adding 24 new reactions, and revising 64 gene-protein-reaction associations based on newly available information. The new model establishes improved predictions of growth yields on native substrates and is capable of correctly predicting the knockout outcomes for 27 out of 28 gene deletion mutants. By tracing a bifurcated electron flow mechanism, the iMAC868 model predicts thermodynamically feasible (co)utilization pathway of methane and bicarbonate using various terminal electron acceptors through the reversal of the aceticlastic pathway.


This effort paves the way in informing the search for thermodynamically feasible ways of (co)utilizing novel carbon substrates in the domain Archaea.


Methane, the second most important greenhouse gas, is regulated primarily by microbial processes [1]. A renewed interest in methane as a gas substrate for the production of biofuels is spearheaded by its abundance in shale gas [25]. At the same time, concerns related with methane’s role as a potent greenhouse gas drives the need to mitigate its adverse environmental impact [6]. Advances in the characterization of microbial consortia in anoxic sediments have revealed the potential of transforming methane into various products through biological routes [710].

The global methane cycle is predominantly controlled by anaerobic methanotrophic archaea (ANME) in anoxic environments [11, 12] and aerobic methanotrophic bacteria at the anoxic–oxic interface of habitats [13, 14]. Aerobic methanotrophy [15], proceeds via the oxidation of methane to methanol by a methane monooxygenase and then to formaldehyde by methanol dehydrogenase, which is subsequently integrated into central carbon metabolism through the ribulose monophosphate or the serine pathway [16]. This scheme, however, requires an initial activation cost in the form of NAD(P)H, which is replenished at the expense of carbon efficiency. Shaped by the paucity of available energy, the anaerobic methanotrophy has been shown to exhibit better carbon and energy efficiency [17, 18]. However, in contrast to the aerobic route, anaerobic methanotrophy is relatively poorly characterized as a result of the difficulties in culturing ANMEs in the lab [11] arising from syntrophy requirements. In such environments, the anaerobic methanotroph oxidizes methane and the microbial partner reduces an electron acceptor, often an inorganic ion such as NO3 [19] or SO4 2− [20]. Despite these difficulties, recent metagenomics analysis of ANMEs has partially revealed the methanotrophic pathways, observed in most methanogenic archaea, demonstrating the phylogenetic relationship between ANMEs and methanogens [21, 22]. Of particular interest is the methanogenic archaeon Methanosarcina acetivorans for which trace methane oxidation has been observed [23, 24] implying that it possesses the necessary pathways and electron flow systems to accomplish methanotrophy. However, a complete reversal of methanogenesis pathway is thermodynamically infeasible unless coupled with an electron-accepting pathway [25, 26].

M. acetivorans, a strictly anaerobic marine methanogen possessing one of the largest known archaeal genomes [27], has emerged as a model archaeon owing to the availability of genetic tools [28] and versatility in substrate utilization [29, 30]. While the pathways describing the metabolism on native substrates have been extensively studied [3134], relevant pathways and electron flows for methane oxidation by M. acetivorans remain largely uncharacterized. Two genome-scale metabolic (GSM) models, iVS941 [35] and iMB745 [36], for this organism have been proposed. Both models, however, are not up to date with the current literature on the stoichiometry of ion transport across the membrane and ATP synthesis [3742]. Recent findings on the electron flow mechanisms of M. acetivorans cell extracts grown with methane (unpublished observations, Zhen Yan and James G. Ferry) motivates an update in the existing genome-scale models to incorporate recent findings and to allow for the analysis of methane utilization in silico.

In this paper we make use of a revamped GSM for M. acetivorans to postulate pathways for reversing methanogenesis while maintaining overall thermodynamic feasibility. We first generated an up-to-date GSM model for M. acetivorans by combining information from two earlier models (i.e., iVS941 and iMB745) along with the most recent data from literature and databases. M. acetivorans has transcriptome and proteome profiles that differ depending on growth substrate [31, 34]. We augmented the updated gene-protein-reaction (GPR) associations with regulatory (i.e., −R) switches to incorporate proteomics data to the updated metabolic reconstruction by switching on/off reactions for different substrates. Using the model as a starting point a thermodynamically feasible pathway is proposed for the co-utilization of methane and bicarbonate in the presence of Fe3+, NO3 , SO4 2−, and MnO2 as external electron acceptors. Overall ΔG ≤ 0 is imposed as a constraint to ensure thermodynamic feasibility of methanogenesis reversal in the presence of an external electron acceptor. The interplay between externally supplied electron acceptors and various by-products is analyzed. The feasibility of methanotrophy by resting cells is assessed when all carbons coming from methane and bicarbonate are converted into acetate, formate, CO2, and methyl sulfide, the known byproducts of M. acetivorans’ metabolism [30, 43] some of which were also observed recently by Wood et al. [44] as end products of methanotrophy by the archaeon.

Results and discussion

Updated genome-scale metabolic model reconstruction of M. acetivorans, iMAC868

iMAC868 contains 868 genes, 845 reactions, and 718 metabolites (Additional file 1) and provides better agreement with the observed growth yields on methanol and acetate compared to earlier reconstructions (see Table 1). Improved prediction is due to the correction of charge and mass imbalances of the reactions inherited from the previous metabolic models, incorporation of accurate ion exchange stoichiometries for membrane-bound reactions, and optimization of Na+/H+ ratio for sodium/proton antiporter (Mrp) and ATP synthase. Among the charge re-balanced reactions, those involving cofactor F420 in the methylotrophic pathway and ATP synthase also required proton rebalancing to accurately account for proton exchange across the cell membrane. The number of Na+ pumped out by ferredoxin-dependent methanophenazine reductase (Rnf) was updated from three Na+ in iMB745 [36] to four Na+ per methanophenazine reduced in accordance with experimental findings [38]. In addition to this, the ATP synthase reaction was modified to co-utilize Na+ and H+ gradients [37]. Using the procedure described in Methods section, the optimal Na+/H+ ratio recapitulating the growth yields on native substrates was estimated to be 2:1 for Mrp and 3:1 and 1:2 for ATP synthase, respectively. The two identified solutions for ATP synthase are rendered equivalent by the reversible 2:1 antiport of Na+/H+ across the cell membrane by Mrp, which makes one intracellular H+ equivalent to two extracellular Na+. We also added an F420-dependent NADP reductase to the iMAC868 model (personal communications with James G. Ferry), which functions as a source of NADPH for cell biosynthesis.

Table 1 Growth yield predictions of iMAC868 model of M. acetivorans compared with predictions of previous models iVS941 [35] and iMB745 [36]

Upon correcting 64 GPRs based on updated gene annotations (Additional file 2) and implementing proteomics-dependent growth condition-specific R-GPR switches, iMAC868 correctly predicts gene knockout outcomes for 27 out of 28 mutants of M. acetivorans (see Table 2; Additional file 1). The only false prediction by the iMAC868 model is missing the in vivo essentiality of the mutant lacking methanol-specific methyltransferases (ΔmtaA1 ΔmtaCB1 ΔmtaCB2 ΔmtaCB3) growing with acetate due to the unknown role of the enzymes in acetate-grown cells [45]. The model correctly captures the essentiality of mch [46] by identifying the role of the methylotrophic pathway as a source of reduced F420 for NADPH production in acetate-grown cells thereby rendering mch (methenyl-H4SPT cyclohydrolase) essential. Comparing with the iVS941 model, the iMAC868 model correctly predicts the essentiality of rnf, mtr, and the membrane-bound hdr due to the updated GPRs and ion transport stoichiometries included in this model.

Table 2 Gene deletion lethality predictions by iMAC868 model of M. acetivorans compared with predictions of previous models

Model customization to capture methanotrophy by M. acetivorans

In order to allow for methanotrophy, the iMAC868 model was customized to enable three new processes: (1) reversal of methyl-coenzyme M reductase (Mcr) reaction, (2) inclusion of a cytosolic methyltransferase (CmtA), and (3) inclusion of a mechanism enabling electron bifurcation and its subsequent discharge to an external electron acceptor. The methyl-coenzyme M reductase of an anaerobic methanotroph (ANME-MCR), capable of oxidizing methane [47], was appended to the iMAC868 model upon deactivating the native Mcr to prevent methanogenesis based on evidence regarding the reversal of methanogenesis in M. acetivorans [23, 24, 48], reversibility of native Mcr [49], and the confirmed heterologous expression of ANME-MCR in M. acetivorans demonstrating ferric-dependent methanotrophy [44]. CmtA [50] serves as a soluble alternative to membrane-bound Mtr, allowing the conversion of methyl-coenzyme M to methyl-tetrahydrosarcinapterin without drawing on sodium ion gradients across the membrane. Cdh, a key enzyme in the modeled pathway, is dependent on ferredoxin to reduce CO2 that generates the carbonyl group in acetyl-CoA [51]. Two flavin-based mechanism are postulated in which an electron pair derived from oxidation of coenzyme B and coenzyme M (E m = −143 mV) bifurcates yielding a high-potential electron reducing Fe3+ to Fe2+ (E m = +770 mV) and a low-potential electron reducing ferredoxin (E m = −420 mV). Flavin-based electron bifurcation is common among anaerobic microbes including methanogens [52, 53]. Although, both postulated mechanisms depend on delivering electrons to Fe3+ on the outer aspect of the cytoplasmic membrane, the bifurcation event occurs either at the cytoplasm or the membrane. Bifurcation in the membrane depends on the Rnf complex, abundant in acetate-grown M. acetivorans, which interacts with ferredoxin and contains two FMN-bound subunits that are possible sites for electron bifurcation [54]. Oxidation of coenzyme B and coenzyme M is catalyzed by the membrane-bound CoMS-SCoB heterodisulfide reductase (HdrDE) [54]. Methanophenazine (MP) is a quinone-like electron carrier that shuttles electrons between HdrDE and the Rnf complex. Importantly, the reduction of ferredoxin is not dependent on a sodium gradient. Bifurcation in the cytoplasm is postulated to be dependent on the fused HdrA2:MvhD protein shown previously to be present in acetate-grown M. acetivorans [34]. Oxidation of HS-CoB and HS-CoM is catalyzed by the soluble heterodisulfide reductase (HdrB2) that donates electrons to the flavin-containing HdrA2 component where bifurcation takes place reducing ferredoxin and transferring an electron to the membrane where reduction of Fe3+ takes place. Finally, an electron transfer reaction is included in the model to transfer the electrons from reduced methanophenazine to an externally supplied electron acceptor based on its reported in vivo essentiality [19, 20, 55]. The essentiality of this reaction was confirmed by the absence of any in silico external electron acceptor-independent thermodynamically feasible metabolic state despite allowing the production of all reported reduced products such as hydrogen gas [56] and organic acids such as acetate and formate [30]. These additions complete the pathways for the oxidation of methane to various end products such as acetate, formate and CO2. The addition of prospective biofuel molecule production pathways for ethanol [57], butanol [58], and isobutanol [59] to the model allows the exploration of their thermodynamically feasible maximum theoretical yields for different electron acceptors.

Products of electron-acceptor-dependent AOM

The model supports acetate, CO2 and biomass as the main products of methanotrophy using all tested electron acceptors. Methane is oxidized to methyl-H4SPT by ANME-MCR and Mtr (or CmtA), a part of which is oxidized via the methylotrophic pathway to produce intracellular CO2. The remaining methyl-H4SPT is used to produce acetyl-CoA, the primary building block for all biomass precursors. Additional carbon fixation occurs via reductive carboxylation by Cdh and Por. Acetyl-CoA is also converted to acetate, generating ATP via substrate-level phosphorylation (Fig. 1). The electrons released upon activation of methane by ANME-MCR are transferred to ferredoxin via soluble and membrane-bound electron transport chains involving flavin-based electron bifurcation mechanisms. In addition, further oxidation of methyl-H4SPT through methylotrophic pathway generates reducing equivalents in the form of ferredoxin and F420. Fpo and Rnf complexes facilitate the transfer of electrons from reduced F420 and ferredoxin, respectively, to the external electron acceptor via methanophenazine, thereby generating H+ and Na+ gradients across the membrane for chemiosmotic ATP synthesis (see Fig. 1). The primary carbon fixation mechanism via reductive carboxylation prompted a quantitative analysis of the impact of utilizing CO2 as a co-substrate in the form of bicarbonate on acetate and biomass yields.

Fig. 1
figure 1

Proposed methanogenesis reversal pathway supported by the iMAC868 model of M. acetivorans for co-metabolization of methane and bicarbonate in the presence of Fe3+ as external electron acceptor. Soluble and membrane-bound electron bifurcation routes are shown as blue and orange, respectively, and enzymes within ovals. In both routes, electrons originate from coenzyme B and coenzyme M. For enzymes with multiple subunits, only the subunits of soluble Hdr and Rnf involved in electron bifurcation are shown in detail. Numbers in italics next to enzyme ovals denote reaction fluxes (in mmol/gDCW-h) calculated under maximization of acetate production at bicarbonate to methane ratio of 0.44. This ratio corresponds to the maximum thermodynamically feasible value ensuring biomass production at 30 % of its theoretical maximum for Fe3+ as the electron acceptor. The flux towards growth was calculated by assuming that 1 g of biomass contains 36 mmol of carbon. Intracellular proton and water stoichiometries are omitted for the sake of simplicity. Soluble methyltransferase (CmtA) is not present in the network since the minimum possible flux through this reaction is zero. Mcr * putative ANME-like Mcr homolog to methyl-coenzyme M reductase, HdrBC:HdrA:MvhD soluble ferredoxin-dependent heterodisulfide reductase, Mtr methyl-THSPT:coenzyme M methyltransferase, Mer methenyl-THSPT reductase, Mtd methenyl-THSPT dehydrogenase, Mch methenyl-THSPT cyclohydrolase, Ftr formylmethanofuran:THSPT formyltransferase, Fmd formylmethylfuran dehydrogenase, Cdh CO dehydrogenase, Pta phosphotransacetylase, Ack acetate kinase, Por pyruvate synthase, Atps ATP synthase, Mrp sodium/proton antiporter, Rnf methanophenazine reductase, Cyt cytochrome c subunit of Rnf complex, Fpo F420 dehydrogenase, Cam carbonic anhydrase, F4nr F420-dependent NADP reductase, THSPT tetrahydrosarcinapterin, MF methanofuran, MP methanophenazine, MPH 2 reduced methanophenazine, Fd o oxidized ferredoxin, Fd r reduced ferredoxin, F 420 coenzyme F420, F 420 H 2 reduced coenzyme F420

Thermodynamic feasibility of methanotrophy is ensured only when the free energy of reduction (ΔGred) of the supplied electron acceptor is less than 50.5 kJ/electron-pair (Fig. 2), corresponding to the maximum free energy equivalents generated by CO2 production (see Table 3). Using methane as the sole carbon source, maximum biomass yield is constrained by thermodynamic feasibility when ΔGred of the electron acceptor is greater than −20 kJ/electron-pair. Sulfate-dependent methanotrophy falls within this regime, in which thermodynamic coupling with an exergonic pathway such as acetate or CO2 production (Table 3) drives only partial conversion of methane to biomass. In contrast, biomass production is limited only by stoichiometry during ferric-dependent methanotrophy due to the far greater free energy equivalents produced by the reduction of Fe3+ to Fe2+ (ΔG = −140.44 kJ/electron-pair) compared to SO4 2− reduction (ΔG = 44.53 kJ/electron-pair). This thermodynamic advantage for Fe3+ allows for the co-utilization of bicarbonate up to a maximum HCO3 /CH4 ratio of 0.44 with complete incorporation of all substrate carbons into biomass. However, the endergonic nature of bicarbonate uptake disallows HCO3 /CH4 co-utilization for increasing biomass yield under sulfate-dependent methanotrophy. Methanotrophy using NO3 and MnO2 also allows for co-utilization of bicarbonate as both electron acceptors have ΔGred greater than −20 kJ/electron-pair (see Table 3).

Fig. 2
figure 2

Biomass yield (per 10 mmol methane) as a function of the ΔG of external electron acceptor reduction (kJ/electron-pair) predicted by the iMAC868 model of M. acetivorans. Solid line methane as the sole carbon source; dashed line bicarbonate and methane (at a ratio of HCO3 /CH4 = 0.44) as carbon sources. Vertical dotted lines show the ΔG (kJ/electron-pair) of reduction for Fe3+ (−140.44), MnO2 (−77.65), NO3 (−68.15), and SO4 2− (+44.53). A magnified insert shows the maximum biomass yields for sulfate-dependent methanotrophy. All ΔG values were evaluated at pH of 7, 25 °C, and an ionic strength of 0.25 M as described by Alberty [64]

Table 3 Oxidation half reactions of methane to various products and reduction half reactions of various electron acceptors

The model predicts a maximum acetate production (0.5 mol/mol-methane), constrained only by stoichiometry for both Fe3+ and SO4 2− during growth on only methane. This yield is further increased to 0.94 mol/mol-methane at an optimal HCO3 /CH4 ratio of 0.88 for ferric-dependent methanotrophy, and 0.68 mol/mol-methane at an optimal HCO3 /CH4 ratio of 0.36 during sulfate-dependent methanotrophy (Fig. 3a). The improvement in acetate yield arises from the reduction in the fraction of methane oxidized via the methylotrophic pathway from 50 to 6 % and 32 % during ferric- and sulfate-dependent methanotrophy, respectively. However, a complete reversal of the aceticlastic pathway with a co-utilization ratio of one could not be achieved using either electron acceptor due to thermodynamic restrictions during sulfate-dependent methanotrophy and reduced ferredoxin availability during ferric-dependent methanotrophy. Under sulfate-dependent methanotrophy, the minimum essential flux through the methylotrophic pathway enables thermodynamic coupling with sulfate reduction for the generation of free energy equivalents. Mandatory channeling of electrons towards ferric ions by the electron bifurcation mechanism decreases available reduced ferredoxin for acetate synthesis during ferric-dependent methanotrophy. Despite the exergonic nature and the ATP generating capability of the acetate production pathway, it is never essential (minimum acetate production is always zero) at any HCO3 /CH4 ratio due to the fact that there exist other competing products and an electron acceptor-driven chemiosmotic ATP synthesis.

Fig. 3
figure 3

Acetate and carbon dioxide production yields as a function of bicarbonate to methane ratio (a, b) and biomass yield (c, d) using Fe3+ (solid lines) or SO4 2− (dashed lines) as external electron acceptors. Methane was the sole carbon source for generating the plots shown in panels (c) and (d). All yields are per 10 mmol of methane

CO2 production remains non-essential during ferric-dependent methanotrophy as revealed by the model (Fig. 3b) due to the fact that reductive carboxylation of acetyl-CoA allows the production of many different thermodynamically feasible products. In contrast, CO2 production for sulfate-dependent methanotrophy beyond a HCO3 /CH4 ratio of 0.36 becomes mandatory. Beyond this ratio, increased CO2 production via the methylotrophic pathway serves to offset the free energy increase associated with uptake of bicarbonate. Carbon channeling towards the methylotrophic pathway leads to increased CO2 production thus decreasing methane flow towards other major products, thereby adversely affecting acetate and biomass yields at HCO3 /CH4 ratios beyond 0.36. The trade-off plot between the products of AOM and biomass did not reveal any thermodynamic restrictions in the solution space during ferric-dependent methanotrophy (Fig. 3c, d). However, the model predicts that acetate becomes thermodynamically constrained beyond a biomass yield of 0.018 for sulfate-dependent methanotrophy. Up to this yield value, the minimum required CO2 production remains zero due to the fact that either acetate or CO2 production pathways can generate the necessary free energy equivalents, ATP and reducing equivalents for biomass production. At biomass yields above 0.018, CO2 production becomes mandatory.

The production of ethanol, butanol, isobutanol, and methanol is thermodynamically feasible through both ferric-dependent and sulfate-dependent methanotrophy (Fig. 4). However, complete carbon conversion of methane to candidate biofuel molecules is thermodynamically feasible only for ferric-dependent methanotrophy with methane as the sole carbon source (Fig. 4a) due to the favorable thermodynamics of coupling the biofuel production pathways by ferric reduction (see Table 3). Upon co-utilization of methane and bicarbonate, electron bifurcation limits the availability of reduced ferredoxin for fixing CO2 by Cdh to produce acetyl-CoA (biofuel precursor), thereby restricting maximum achievable biofuel yield (Fig. 4a). Moreover, biofuel production pathways require additional energy in the form of NAD(P)H necessitating elevated amounts of reduced F420 at increasing bicarbonate to methane ratios which is also controlled by electron bifurcation. Incorporation of bicarbonate into methanol occurred via the CO2 reduction pathway (reversal of the methylotrophic pathway) as opposed to CO2 reduction by acetyl-CoA synthesis, causing all electrons to be generated by the ANME-MCR. During sulfate-dependent methanotrophy, none of the products could be produced with the complete carbon conversion efficiency due to the fact that coupling biofuel production with SO4 2− reduction remains thermodynamically infeasible (see Table 3) requiring the co-production of by-products such as acetate or CO2. As a consequence of this, co-utilization of bicarbonate and methane is not supported (see Fig. 4b). Although both oxidation of methane to methanol and reduction of sulfate to sulfide are thermodynamically infeasible, methanol can be still produced with SO4 2− due to coupling with the concomitant production of CO2. The lower bound for the production of all biofuel molecules is zero indicating that their production is not growth-coupled when methane is either the sole carbon source or co-utilized with bicarbonate.

Fig. 4
figure 4

Biofuel yields using methane (shaded in black) as the sole carbon source and HCO3 /CH4 (shaded in gray) with a ratio of 0.44 for ferric-dependent (a) and sulfate-dependent (b) methanotrophy. The dashed bars denote the excess carbons that could not be incorporated into the product due to thermodynamic restrictions. Ferric-dependent methanotrophy predicted increased yields at higher HCO3 /CH4 ratios, but the carbon conversion was less than 100 %. Sulfate-dependent methanotrophy cannot achieve 100 % carbon conversion at any HCO3 /CH4 ratio due to mandatory co-production of thermodynamically feasible by-products. Product yields are in mol per mol methane

Interplay between electron acceptors and by-products of AOM at no growth

The interplay between the external electron acceptor choice and various products of AOM is pictorially illustrated (see Fig. 5) using feasible production envelopes for growth-arrested cells. We constrained the model for zero growth, ATP production for only maintenance requirements, and bicarbonate to methane ratio of 0.44. Analysis of the product profiles predicted by the iMAC868 model, based on the imposed constraints, identifies acetate as the main product of co-utilization of methane with bicarbonate along with the possible production of formate, CO2, and methyl sulfide (Fig. 5). A minimum Fe3+ uptake (i.e., 0.5 mol/mol-methane) is necessary to maintain thermodynamic feasibility of the observed solution spaces (Fig. 5a–d) at which methyl sulfide is found to be essential (Fig. 5a) due to the fact that it is the least oxidized by-product of methanotrophy by M. acetivorans. The maximum methyl sulfide yield at this Fe3+ uptake exceeds methane uptake indicating that bicarbonate is reduced via the methylotrophic pathway. The reversal of the methylotrophic pathway, however, is limited by the availability of reducing equivalents provided by Mcr, thereby resulting in an incomplete conversion of substrate carbons (methane and bicarbonate) to methyl sulfide. Increasing Fe3+ uptake allows more flux through the methylotrophic pathway, thereby generating additional intracellular CO2 for an increased acetate production of up to a maximum of 0.71 mol/mol-methane at an Fe3+ uptake of 2.2 mol/mol-methane (Fig. 5b). At this uptake rate, acetate can be produced as the sole product of methanotrophy resulting in non-essentiality of methyl sulfide production. Beyond this Fe3+ uptake rate, acetate production decreases due to the paucity of methyl-coenzyme M arising from increased flux through the methylotrophic pathway and channeling of electrons towards Fe3+ reduction via the membrane-bound electron transport chain. This increase in methylotrophic pathway flux also increases the yield of formate, an intermediate of this pathway. Maximum formate yield is found to be 1.44 mol/mol-methane at a Fe3+ uptake of 5.1 mol/mol-methane (Fig. 5c) where all taken up carbons are converted to formate. Beyond this Fe3+ uptake rate, CO2 production becomes essential so as to generate sufficient electrons for reduction of Fe3+ (Fig. 5d). A consequence of the essentiality of CO2 is the reduction of maximum formate yield. At a maximum Fe3+ uptake of 8 mol/mol-methane, only CO2 is produced due to the fact that it is the most oxidized form of carbon that can be produced by M. acetivorans.

Fig. 5
figure 5

Thermodynamically feasible production envelope (highlighted in grey) of methyl sulfide (a), acetate (b), formate (c), and carbon dioxide (d) as a function of Fe3+ uptake predicted by the iMAC868 model under no growth. All product yields and Fe3+ uptake are defined as mol per mol methane


An updated genome-scale metabolic reconstruction (iMAC868) for the archaeon Methanosarcina acetivorans that integrates the latest literature findings and provides complete pathways and electron flow systems for reversing methanogenesis is introduced. Upon improving growth and gene-knockout outcome prediction for M. acetivorans grown on its native substrates, the model was used in a prospective mode for assessing thermodynamically feasible methanotrophic pathways leading to the production of biofuel candidate molecules such as methanol, ethanol, butanol, and isobutanol. We found that anaerobic methanotrophy favored the production of acetate and CO2 as they provide free energy equivalents to support growth. Co-utilization of CO2 (in the form of bicarbonate) and methane was feasible for certain ratios leading to improved carbon yields for acetate and biofuel molecules. Re-routing of a fraction of carbon towards CO2 was found to be a recurring mechanism for driving growth and production within thermodynamically constrained metabolic states. Finally, the availability in excess of inorganic electron acceptors resulted in a switch between substrate-level and chemiosmotic ATP synthesis. Thermodynamic constraints were often the limiting factor in product yields. Unsurprisingly, the inability to completely reverse the aceticlastic pathway in the absence of an external electron acceptor was confirmed.

The theoretical limits of external electron acceptor utilization possibilities to drive the reversal of methanogenesis were thoroughly explored. Key challenges that confound the fidelity of model predictions are still unknown sodium gradient requirements, substrate-dependent regulation and the detailed mechanism of electron transport from internal electron carriers to their external counterparts. Shedding light to these questions will require systematic experimental investigations to confirm or refute electron flow paths guided by the rapidly expanding modeling infrastructure.


Model assembly and growth/product formation predictions

The updated genome-scale metabolic model for M. acetivorans, iMAC868, was constructed by appending genes and reactions from iVS941 to iMB745. The stoichiometric coefficients of the ions associated with the membrane-bound electron transport chain were updated based on recent findings. This includes Rnf, ATP synthase and the Na+/H+ antiporter Mrp. All reactions in the model were checked for mass and charge balances and corrected if necessary. Five reactions from amino acid, two from tRNA, one from cofactor biosynthesis pathways, three from methanogenesis, and three metabolite transport reactions required elemental and charge rebalancing. The biomass equation formulation of iMB745 model was adopted in the new iMAC868 model as reported previously [36]. It was ensured that all biomass precursors could be independently produced in the model to avoid feasibiltiy tolerance related errors [60]. The iMAC868 model is available in Excel format in Additional file 1. All reaction fluxes are in mmol/gDCW-h except for the reaction representing cell biomass formation that is expressed in h−1. The medium composition was assumed to be a defined high-salt medium [29]. The model was assembled in a format compatible for flux balance analysis [61]. FBA optimization problems were solved by GNU Linear Programming Kit (GLPK) ( and Gurobi ( solvers in Matlab using COBRA toolbox [62]. Flux variability analysis (FVA) was performed to obtain range of fluxes under optimal growth conditions as described previously [63]. Both FBA and FVA problems incorporated overall thermodynamic feasibility constraints (overall ΔG ≤ 0). FBA was performed by solving the following Linear Programming (LP) problem:

$$\begin{array}{*{20}c} {Maximize } & {v_{biomass} } & {} & {} \\ {Subject \, to} & {\mathop \sum \limits_{j} S_{ij} v_{j} = 0,} & {\forall i \in I, j \in J} & {(1)} \\ {} & {\mathop \sum \limits_{j} \varDelta G_{j} v_{j} \le 0,} & {\forall j \in J^{ex} \mathop \cup \nolimits \, \{ biomass\} } & {(2)} \\ {} & {LB_{j} \le v_{j} \le UB_{j} ,} & {\forall j \in J\backslash \{ Ex - methane, \,ATPM\} } & {(3)} \\ {} & {v_{Ex - methane} = - 10} & {} & {(4)} \\ {} & {v_{ATPM} \ge 2.5} & {} & {(5)} \\ {} & {\varDelta G_{j} = \varDelta G_{j}^{met} , } & {\forall j \in J^{ex} \mathop \cup \nolimits \, \{ biomass\} } & {(6)} \\ {} & {v_{j} \in {\mathbb{R}}} & {\forall j \in J} & {} \\ \end{array}$$

where sets, variables, and parameters are defined as follows:


  • I = {i|i = 1,2,…,M} = Set of metabolites in the stoichiometric model

  • J = {j|j = 1,2,…,N} = Set of reactions in the stoichiometric model

  • J ex = {j|j = 1,2,…,N ex} = Set of exchange reactions in the stoichiometric model


  • v j  = Flux of reaction \(j \in J\)

  • v biomass  = Flux of the biomass formation reaction


  • S ij  = Stoichiometric coefficient of metabolite \(i \in I\) in reaction \(j \in J\)

  • UB j  = Upper bound for the flux of reaction \(j \in J\)

  • LB j  = Lower bound for the flux of reaction \(j \in J\)

  • \(v_{Ex - methane}\) = Flux of the methane exchange reaction

  • \(v_{ATPM}\) = Flux of the non-growth associated maintenance ATP reaction

  • \(\varDelta G_{j}^{met}\) = The \(\varDelta G\) of formation of the metabolite associated with exchange reaction \(j \in J^{ex} \mathop \cup \nolimits \{ biomass\}\)

All standard transformed ΔG values were calculated at pH of 7, temperature of 25 °C and ionic concentration of 0.25 M [64] listed in Additional file 1. The upper bound of the free energy of biomass formation (ΔGbiomass) is estimated to be 3750 J/gDCW. This value is identified so as the overall stoichiometries for growth on acetate and methanol (see below) remain thermodynamically feasible:

$$10 \,CH_{3} COOH \to 9.3\, CH_{4} + 9.3 \,CO_{2} + 0.026\,\, biomass$$
$$10 \,CH_{3} OH \to 6.4 \,CH_{4} + 1.5\, CO_{2} + 6 \,H_{2} O + 0.052 \,\,biomass$$

The iMAC868 model arrived at these overall conversion stoichiometries by solving the FBA problem subject to constraints (1), (3), and (5) only.

In the above LP problem, the flux of the biomass is maximized subject to the constraints of stoichiometry (1), thermodynamics (2), metabolic network fluxes (3), fixed uptake of methane (4), minimum requirements of maintenance ATP of 2.5 mmol/gDCW-h (5), and fixed ΔG values of input/output metabolites to/from the system (6). Flux ranges for target products were obtained by iteratively solving the above LP problem to minimize and maximize all v j separately subject to stoichiometric and thermodynamic constraints.

Formulation of R-GPR to integrate ‘-omics’ data into the metabolic model

Gene-protein-reactions (GPRs) associations in the iMAC868 model were thoroughly assessed and 64 GPRs were corrected using a list of 781 newly revised gene annotations (Additional file 2) along with database entries from KEGG [65], MetaCyc [66], BRENDA [67], and TransportDB [68]. Following this, R-GPR switches were implemented using a dataset of quantitative protein levels for over 250 genes of M. acetivorans grown with acetate and methanol [34]. This dataset is given in Additional file 1. The R-GPR approach allows for the incorporation of ‘omics’ data for conditional switching on/off of reactions allowing improved gene-knockout predictions by providing an insight into the likelihood of a reaction to be active or inactive under specific growth conditions. The following systematic procedure elaborates this conditional activation/inactivation of reactions by the R-GPR switches:

Step 1:

Calculate the ratio of protein abundance under different growth substrates. For each gene k in the total gene set K, the value c k was calculated as the ratio of protein abundance for cells grown on substrate 1 to cells grown on substrate 2 where substrate 1 and 2 can be any of acetate or methanol.

Step 2:

Compare the ratio c k with a pre-defined cutoff value (i.e., equal to 25 %). If the ratio c k is below or equal to the cutoff value then gene k is added to a candidate list G for which the feasibility of removing the corresponding reactions is evaluated.

Step 3:

Re-evaluate GPRs within list G. The GPR for each reaction is re-evaluated assuming that all genes in list G are eliminated. If re-evaluation of the GPR reveals no associated gene then the reaction is added to the set J exp that contains the candidate reactions for removal.

Step 4:

The following mathematical formulation is used to identify the maximum number of reactions in the set J exp that can be removed from the model without dropping the biomass yield below the experimental value:

$$\begin{array}{*{20}c} {Minimize} & {\mathop \sum \limits_{{j \in J^{exp} }} y_{j}} & {} & {} \\ {Subject \, to} & {\mathop \sum \limits_{j} S_{ij} v_{j} = 0,} & {\forall i \in I, j \in J} & {(7)} \\ {} & {y_{j} LB_{j} \le v_{j} \le y_{j} UB_{j} ,} & {\forall j \in J^{exp} \backslash \left( {J^{on} \mathop \cup \nolimits J^{off} } \right)} & {(8)} \\ \begin{aligned} \hfill \\ \hfill \\ \hfill \\ \hfill \\ \end{aligned} & \begin{aligned} LB_{j} \le v_{j} \le UB_{j} , \hfill \\ v_{j} = 0 \hfill \\ v_{biomass} \ge v_{biomass,exp} \hfill \\ v_{j} \in {\mathbb{R}},\text{ }y_{j} \in \{ 0,1\} \hfill \\ \end{aligned} & \begin{aligned} \forall j \in \left( {J\backslash J^{exp} } \right)\mathop \cup \nolimits J^{on} \hfill \\ \forall j \in J^{off} \hfill \\ \hfill \\ \forall j \in J \hfill \\ \end{aligned} & \begin{aligned} (9) \hfill \\ (10) \hfill \\ (11) \hfill \\ \hfill \\ \hfill \\ \end{aligned} \\ \end{array}$$

where sets, variables, and parameters in this MILP problem have the same definition as those defined earlier for the LP problem except for the followings:


  • J exp = {j|j = 1, 2, …, N exp} = Set of reactions whose GPRs are evaluated due to availability of experimental data (i.e., proteomic data)

  • J on = {j|j = 1, 2, …, N on} = Set of reactions for which the evaluation of their GPR indicates active reactions

  • J off = {j|j = 1, 2, …, N off} = Set of reactions for which the evaluation of their GPR indicates inactive reactions


  • v biomass  = Flux of the reaction producing biomass

  • y j  = Binary variable associated with flux of v j


  • v biomass,exp  = The experimentally measured biomass yield for the given substrate

Here, the sum of binary variables y j is minimized subject to the limitations on stoichiometry (7), the flux of the reactions for which proteomic data is available but their GPRs are not evaluated because c k is greater than the cutoff value (8), the flux of all other metabolic network reactions which do not belong to set J exp along with those belong to set J on (9), and the flux reactions which belong to set J off (10), and the flux of the reaction producing biomass being equal or greater than the experimentally measured biomass yield for the given substrate (11).

Representation of external electron acceptor in the metabolic network of iMAC868

Electron acceptor reactions are modeled using “electron acceptor equivalents (EAE)”, which serves to drain electron pairs from the model. Each electron pair is drained from reduced methanophenazine. The general form of the electron acceptor reaction used in this model is:

$$Reduced\,\,methanophenazine \,+\, aH^{ + } [c] \to Oxidized \,\,methanophenazine \,+\, EAE \,+\, bH^{ + } [e]$$

here, a and b correspond to the excess protons drained from the cytosol and secreted into the extracellular medium, respectively. These values are electron acceptor specific and can be obtained from the corresponding electron acceptor reduction reactions described in Table 3. In order to make EAE a balanced metabolite, an exchange reaction describing the draining of EAE is added to the model, with ΔGEAE corresponding to the ΔGred (kJ/electron-pair) of the specific electron acceptor described in Table 3.

Optimization of Na+/H+ ratios of ATP synthase and Mrp antiporter

To find the optimal Na+/H+ stoichiometric ratios of ATP synthase and Mrp matching best the observed growth yields, Flux Balance Analysis (FBA) was used as described earlier in this section. FBA optimization problems were solved at varying Na+/H+ ratios of 3:1, 1:1, 1:3, 2:1, and 1:2 for ATP synthase and at varying Na+/H+ ratios of 3:1, 2:1, 1:1, 1:2, and 1:3 for Mrp. These ratios were selected based on recent findings involving archaeal ATP synthase and their dependence on Mrp. Welte and Deppenmeier [42] estimated that 3–4 translocated ions (H+ or Na+) are required to synthesize one molecule of ATP based on the measured electrochemical ion gradients available for a limited number of methanogens. Jasso-Chavez et al. [40] proposed that Mrp is used for optimization of the thermodynamic efficiency of the ATP synthase in M. acetivorans. The objective function was the maximization of cellular growth. This procedure was repeated for both acetate and methanol as sole carbon sources. The sum squared error (SSE) between the predicted and observed growth yields were calculated and results are shown in Table 4.

Table 4 SSE between predicted (by iMAC868) and observed growth yields for acetate- and methanol-grown cells at different ratios of Na+ and H+ exchange by Mrp antiporter and varying ratios of Na+ and H+ uptake by ATP synthase



anaerobic methanotrophic archaea


anaerobic oxidation of methane


genome-scale metabolic model


in silico Methanosarcina acetivorans metabolic model containing 868 genes


gene-protein-reaction associations


regulatory gene-protein-reaction associations


electron acceptor equivalent


  1. Conrad R. The global methane cycle: recent advances in understanding the microbial processes involved. Environ Microbiol Rep. 2009;1:285–92.

    Article  CAS  Google Scholar 

  2. Wang Q, Chen X, Jha AN, Rogers H. Natural gas from shale formation—the evolution, evidences and challenges of shale gas revolution in United States. Renew Sustain Energy Rev. 2014;30:1–28.

    Article  Google Scholar 

  3. McGlade C, Speirs J, Sorrell S. Unconventional gas—a review of regional and global resource estimates. Energy. 2013;55:571–84.

    Article  Google Scholar 

  4. Finkel ML, Law A. The rush to drill for natural gas: a public health cautionary tale. Am J Public Health. 2011;101:784–5.

    Article  Google Scholar 

  5. Hu D, Xu S. Opportunity, challenges and policy choices for China on the development of shale gas. Energy Policy. 2013;60:21–6.

    Article  Google Scholar 

  6. Pierrehumbert RT. Short-lived climate pollution. Annu Rev Earth Planet Sci. 2014;42:341–79.

    Article  CAS  Google Scholar 

  7. Wegener G, Niemann H, Elvert M, Hinrichs KU, Boetius A. Assimilation of methane and inorganic carbon by microbial communities mediating the anaerobic oxidation of methane. Environ Microbiol. 2008;10:2287–98.

    Article  CAS  Google Scholar 

  8. Zehnder AJ, Brock TD. Methane formation and methane oxidation by methanogenic bacteria. J Bacteriol. 1979;137:420–32.

    CAS  Google Scholar 

  9. Hinrichs KU, Hayes JM, Sylva SP, Brewer PG, DeLong EF. Methane-consuming archaebacteria in marine sediments. Nature. 1999;398:802–5.

    Article  CAS  Google Scholar 

  10. Hoehler TM, Alperin MJ, Albert DB, Martens CS. Field and laboratory studies of methane oxidation in an anoxic marine sediment: evidence for a methanogen-sulfate reducer consortium. Global Biogeochem Cycles. 1994;8:451–63.

    Article  CAS  Google Scholar 

  11. Knittel K, Boetius A. Anaerobic oxidation of methane: progress with an unknown process. Annu Rev Microbiol. 2009;63:311–34.

    Article  CAS  Google Scholar 

  12. Thauer RK, Kaster AK, Seedorf H, Buckel W, Hedderich R. Methanogenic archaea: ecologically relevant differences in energy conservation. Nat Rev Microbiol. 2008;6:579–91.

    Article  CAS  Google Scholar 

  13. Nazaries L, Murrell JC, Millard P, Baggs L, Singh BK. Methane, microbes and models: fundamental understanding of the soil methane cycle for future predictions. Environ Microbiol. 2013;15:2395–417.

    Article  CAS  Google Scholar 

  14. Conrad R. Soil microorganisms as controllers of atmospheric trace gases (H2, CO, CH4, OCS, N2O, and NO). Microbiol Rev. 1996;60:609–40.

    CAS  Google Scholar 

  15. Hakemian AS, Rosenzweig AC. The biochemistry of methane oxidation. Annu Rev Biochem. 2007;76:223–41.

    Article  CAS  Google Scholar 

  16. Lieberman RL, Rosenzweig AC. Biological methane oxidation: regulation, biochemistry, and active site structure of particulate methane monooxygenase. Crit Rev Biochem Mol Biol. 2004;39:147–64.

    Article  CAS  Google Scholar 

  17. Haynes CA, Gonzalez R. Rethinking biological activation of methane and conversion to liquid fuels. Nat Chem Biol. 2014;10:331–9.

    Article  CAS  Google Scholar 

  18. Mueller TJ, Grisewood MJ, Nazem-Bokaee H, Gopalakrishnan S, Ferry JG, Wood TK, Maranas CD. Methane oxidation by anaerobic archaea for conversion to liquid fuels. J Ind Microbiol Biotechnol. 2015;42:391–401.

    Article  CAS  Google Scholar 

  19. Haroon MF, Hu S, Shi Y, Imelfort M, Keller J, Hugenholtz P, Yuan Z, Tyson GW. Anaerobic oxidation of methane coupled to nitrate reduction in a novel archaeal lineage. Nature. 2013;500:567–70.

    Article  CAS  Google Scholar 

  20. Valentine DL, Reeburgh WS. New perspectives on anaerobic methane oxidation. Environ Microbiol. 2000;2:477–84.

    Article  CAS  Google Scholar 

  21. Wang FP, Zhang Y, Chen Y, He Y, Qi J, Hinrichs KU, Zhang XX, Xiao X, Boon N. Methanotrophic archaea possessing diverging methane-oxidizing and electron-transporting pathways. ISME J. 2014;8:1069–78.

    Article  CAS  Google Scholar 

  22. Meyerdierks A, Kube M, Kostadinov I, Teeling H, Glockner FO, Reinhardt R, Amann R. Metagenome and mRNA expression analyses of anaerobic methanotrophic archaea of the ANME-1 group. Environ Microbiol. 2010;12:422–39.

    Article  CAS  Google Scholar 

  23. Moran JJ, House CH, Freeman KH, Ferry JG. Trace methane oxidation studied in several Euryarchaeota under diverse conditions. Archaea. 2005;1:303–9.

    Article  CAS  Google Scholar 

  24. Moran JJ, House CH, Thomas B, Freeman KH. Products of trace methane oxidation during nonmethyltrophic growth by Methanosarcina. J Geophys Res Biogeosci 2007;112:G02011.

    Article  Google Scholar 

  25. Thauer RK. Anaerobic oxidation of methane with sulfate: on the reversibility of the reactions that are catalyzed by enzymes also involved in methanogenesis from CO2. Curr Opin Microbiol. 2011;14:292–9.

    Article  CAS  Google Scholar 

  26. Thauer RK, Shima S. Methane as fuel for anaerobic microorganisms. Ann N Y Acad Sci. 2008;1125:158–70.

    Article  CAS  Google Scholar 

  27. Galagan JE, Nusbaum C, Roy A, Endrizzi MG, Macdonald P, FitzHugh W, Calvo S, Engels R, Smirnov S, Atnoor D, et al. The genome of M. acetivorans reveals extensive metabolic and physiological diversity. Genome Res. 2002;12:532–42.

    Article  CAS  Google Scholar 

  28. Kohler PR, Metcalf WW. Genetic manipulation of Methanosarcina spp. Front Microbiol. 2012;3:259.

    Article  Google Scholar 

  29. Sowers KR, Baron SF, Ferry JG. Methanosarcina acetivorans sp. nov., an acetotrophic methane-producing bacterium isolated from marine sediments. Appl Environ Microbiol. 1984;47:971–8.

    CAS  Google Scholar 

  30. Rother M, Metcalf WW. Anaerobic growth of Methanosarcina acetivorans C2A on carbon monoxide: an unusual way of life for a methanogenic archaeon. Proc Natl Acad Sci USA. 2004;101:16929–34.

    Article  CAS  Google Scholar 

  31. Li Q, Li L, Rejtar T, Karger BL, Ferry JG. Proteome of Methanosarcina acetivorans Part II: comparison of protein levels in acetate- and methanol-grown cells. J Proteome Res. 2005;4:129–35.

    Article  CAS  Google Scholar 

  32. Li Q, Li L, Rejtar T, Karger BL, Ferry JG. Proteome of Methanosarcina acetivorans Part I: an expanded view of the biology of the cell. J Proteome Res. 2005;4:112–28.

    Article  CAS  Google Scholar 

  33. Lessner DJ, Li L, Li Q, Rejtar T, Andreev VP, Reichlen M, Hill K, Moran JJ, Karger BL, Ferry JG. An unconventional pathway for reduction of CO2 to methane in CO-grown Methanosarcina acetivorans revealed by proteomics. Proc Natl Acad Sci USA. 2006;103:17921–6.

    Article  CAS  Google Scholar 

  34. Li L, Li Q, Rohlin L, Kim U, Salmon K, Rejtar T, Gunsalus RP, Karger BL, Ferry JG. Quantitative proteomic and microarray analysis of the archaeon Methanosarcina acetivorans grown with acetate versus methanol. J Proteome Res. 2007;6:759–71.

    Article  CAS  Google Scholar 

  35. Satish Kumar V, Ferry JG, Maranas CD. Metabolic reconstruction of the archaeon methanogen Methanosarcina Acetivorans. BMC Syst Biol. 2011;5:28.

    Article  Google Scholar 

  36. Benedict MN, Gonnerman MC, Metcalf WW, Price ND. Genome-scale metabolic reconstruction and hypothesis testing in the methanogenic archaeon Methanosarcina acetivorans C2A. J Bacteriol. 2012;194:855–65.

    Article  CAS  Google Scholar 

  37. Schlegel K, Leone V, Faraldo-Gomez JD, Muller V. Promiscuous archaeal ATP synthase concurrently coupled to Na+ and H+ translocation. Proc Natl Acad Sci USA. 2012;109:947–52.

    Article  CAS  Google Scholar 

  38. Schlegel K, Welte C, Deppenmeier U, Muller V. Electron transport during aceticlastic methanogenesis by Methanosarcina acetivorans involves a sodium-translocating Rnf complex. FEBS J. 2012;279:4444–52.

    Article  CAS  Google Scholar 

  39. Schlegel K, Muller V. Evolution of Na(+) and H(+) bioenergetics in methanogenic archaea. Biochem Soc Trans. 2013;41:421–6.

    Article  CAS  Google Scholar 

  40. Jasso-Chavez R, Apolinario EE, Sowers KR, Ferry JG. MrpA functions in energy conversion during acetate-dependent growth of Methanosarcina acetivorans. J Bacteriol. 2013;195:3987–94.

    Article  CAS  Google Scholar 

  41. Suharti S, Wang M, de Vries S, Ferry JG. Characterization of the RnfB and RnfG subunits of the Rnf complex from the archaeon Methanosarcina acetivorans. PLoS One. 2014;9:e97966.

    Article  Google Scholar 

  42. Welte C, Deppenmeier U. Bioenergetics and anaerobic respiratory chains of aceticlastic methanogens. Biochim Biophys Acta. 2014;1837:1130–47.

    Article  CAS  Google Scholar 

  43. Moran JJ, House CH, Vrentas JM, Freeman KH. Methyl sulfide production by a novel carbon monoxide metabolism in Methanosarcina acetivorans. Appl Environ Microbiol. 2008;74:540–2.

    Article  CAS  Google Scholar 

  44. Wood Tk, Soo V, McAnulty M, Tripathy A, Zhu F, Zhang L, Hatzakis E, Smith P, Agrawal S, Nazem-Bokaee H, et al. Reversing methanogenesis to capture methane for liquid biofuels precursors. Microb Cell Fact. 2015. doi:10.1186/s12934-015-0397-z.

  45. Bose A, Pritchett MA, Metcalf WW. Genetic analysis of the methanol- and methylamine-specific methyltransferase 2 genes of Methanosarcina acetivorans C2A. J Bacteriol. 2008;190:4017–26.

    Article  CAS  Google Scholar 

  46. Guss AM, Mukhopadhyay B, Zhang JK, Metcalf WW. Genetic analysis of mch mutants in two Methanosarcina species demonstrates multiple roles for the methanopterin-dependent C-1 oxidation/reduction pathway and differences in H(2) metabolism between closely related species. Mol Microbiol. 2005;55:1671–80.

    Article  CAS  Google Scholar 

  47. Shima S, Krueger M, Weinert T, Demmer U, Kahnt J, Thauer RK, Ermler U. Structure of a methyl-coenzyme M reductase from Black Sea mats that oxidize methane anaerobically. Nature. 2012;481:98–101.

    Article  CAS  Google Scholar 

  48. Moran JJ, Beal EJ, Vrentas JM, Orphan VJ, Freeman KH, House CH. Methyl sulfides as intermediates in the anaerobic oxidation of methane. Environ Microbiol. 2008;10:162–73.

    CAS  Google Scholar 

  49. Scheller S, Goenrich M, Boecher R, Thauer RK, Jaun B. The key nickel enzyme of methanogenesis catalyses the anaerobic oxidation of methane. Nature. 2010;465:606–8.

    Article  CAS  Google Scholar 

  50. Vepachedu VR, Ferry JG. Role of the fused corrinoid/methyl transfer protein CmtA during CO-dependent growth of Methanosarcina acetivorans. J Bacteriol. 2012;194:4161–8.

    Article  CAS  Google Scholar 

  51. Terlesky KC, Ferry JG. Ferredoxin requirement for electron transport from the carbon monoxide dehydrogenase complex to a membrane-bound hydrogenase in acetate-grown Methanosarcina thermophila. J Biol Chem. 1988;263:4075–9.

    CAS  Google Scholar 

  52. Buckel W, Thauer RK. Energy conservation via electron bifurcating ferredoxin reduction and proton/Na(+) translocating ferredoxin oxidation. Biochim Biophys Acta. 2013;1827:94–113.

    Article  CAS  Google Scholar 

  53. Kaster AK, Moll J, Parey K, Thauer RK. Coupling of ferredoxin and heterodisulfide reduction via electron bifurcation in hydrogenotrophic methanogenic archaea. Proc Natl Acad Sci USA. 2011;108:2981–6.

    Article  CAS  Google Scholar 

  54. Li Q, Li L, Rejtar T, Lessner DJ, Karger BL, Ferry JG. Electron transport in the pathway of acetate conversion to methane in the marine archaeon Methanosarcina acetivorans. J Bacteriol. 2006;188:702–10.

    Article  CAS  Google Scholar 

  55. Beal EJ, House CH, Orphan VJ. Manganese- and iron-dependent marine methane oxidation. Science. 2009;325:184–7.

    Article  CAS  Google Scholar 

  56. Lovley DR, Ferry JG. Production and Consumption of H(2) during Growth of Methanosarcina spp. on Acetate. Appl Environ Microbiol. 1985;49:247–9.

    CAS  Google Scholar 

  57. Lee SK, Chou H, Ham TS, Lee TS, Keasling JD. Metabolic engineering of microorganisms for biofuels production: from bugs to synthetic biology to fuels. Curr Opin Biotechnol. 2008;19:556–63.

    Article  CAS  Google Scholar 

  58. Atsumi S, Cann AF, Connor MR, Shen CR, Smith KM, Brynildsen MP, Chou KJ, Hanai T, Liao JC. Metabolic engineering of Escherichia coli for 1-butanol production. Metab Eng. 2008;10:305–11.

    Article  CAS  Google Scholar 

  59. Atsumi S, Li Z, Liao JC. Acetolactate synthase from Bacillus subtilis serves as a 2-ketoisovalerate decarboxylase for isobutanol biosynthesis in Escherichia coli. Appl Environ Microbiol. 2009;75:6306–11.

    Article  CAS  Google Scholar 

  60. Chindelevitch L, Trigg J, Regev A, Berger B. An exact arithmetic toolbox for a consistent and reproducible structural analysis of metabolic network models. Nat Commun. 2014;5:4893.

    Article  CAS  Google Scholar 

  61. Orth JD, Thiele I, Palsson BO. What is flux balance analysis? Nat Biotechnol. 2010;28:245–8.

    Article  CAS  Google Scholar 

  62. Schellenberger J, Que R, Fleming RM, Thiele I, Orth JD, Feist AM, Zielinski DC, Bordbar A, Lewis NE, Rahmanian S, et al. Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox v2.0. Nat Protoc. 2011;6:1290–307.

    Article  CAS  Google Scholar 

  63. Mahadevan R, Schilling CH. The effects of alternate optimal solutions in constraint-based genome-scale metabolic models. Metab Eng. 2003;5:264–76.

    Article  CAS  Google Scholar 

  64. Alberty RA. Calculation of standard transformed Gibbs energies and standard transformed enthalpies of biochemical reactants. Arch Biochem Biophys. 1998;353:116–30.

    Article  CAS  Google Scholar 

  65. Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M. KEGG: Kyoto Encyclopedia of genes and genomes. Nucleic Acids Res. 1999;27:29–34.

    Article  CAS  Google Scholar 

  66. Caspi R, Altman T, Dreher K, Fulcher CA, Subhraveti P, Keseler IM, Kothari A, Krummenacker M, Latendresse M, Mueller LA, et al. The MetaCyc database of metabolic pathways and enzymes and the BioCyc collection of pathway/genome databases. Nucleic Acids Res. 2012;40:D742–53.

    Article  CAS  Google Scholar 

  67. Schomburg I, Chang A, Placzek S, Sohngen C, Rother M, Lang M, Munaretto C, Ulas S, Stelzer M, Grote A, et al. BRENDA in 2013: integrated reactions, kinetic data, enzyme function data, improved disease classification: new options and contents in BRENDA. Nucleic Acids Res. 2013;41:D764–72.

    Article  CAS  Google Scholar 

  68. Ren Q, Chen K, Paulsen IT. TransportDB: a comprehensive database resource for cytoplasmic membrane transport systems and outer membrane channels. Nucleic Acids Res. 2007;35:D274–9.

    Article  CAS  Google Scholar 

  69. Summer H. Improved approach for transferring and cultivating Methanosarcina acetivorans C2A (DSM 2834). Lett Appl Microbiol. 2009;48:786–9.

    CAS  Google Scholar 

  70. Buan NR, Metcalf WW. Methanogenesis by Methanosarcina acetivorans involves two structurally and functionally distinct classes of heterodisulfide reductase. Mol Microbiol. 2010;75:843–53.

    Article  CAS  Google Scholar 

  71. Oelgeschlager E, Rother M. In vivo role of three fused corrinoid/methyl transfer proteins in Methanosarcina acetivorans. Mol Microbiol. 2009;72:1260–72.

    Article  Google Scholar 

  72. Mahapatra A, Srinivasan G, Richter KB, Meyer A, Lienard T, Zhang JK, Zhao G, Kang PT, Chan M, Gottschalk G, et al. Class I and class II lysyl-tRNA synthetase mutants and the genetic encoding of pyrrolysine in Methanosarcina spp. Mol Microbiol. 2007;64:1306–18.

    Article  CAS  Google Scholar 

  73. Rother M, Boccazzi P, Bose A, Pritchett MA, Metcalf WW. Methanol-dependent gene expression demonstrates that methyl-coenzyme M reductase is essential in Methanosarcina acetivorans C2A and allows isolation of mutants with defects in regulation of the methanol utilization pathway. J Bacteriol. 2005;187:5552–9.

    Article  CAS  Google Scholar 

Download references

Author's contributions

HN-B reconstructed the metabolic network, formulated the LP problems, performed the simulations and analyses, designed and generated the figures and tables, and wrote the manuscript. SG contributed to the reconstruction and analysis of the metabolic network and interpretation of simulation results, participated in designing the figures and tables, and revised the manuscript. JGF provided the data for Additional File 2, contributed to the envisioning and generation of Figure 1, proposed the electron bifurcation mechanism and interpretation of the corresponding simulations results, and contributed in writing of the manuscript. TKW provided insights to the interpretation of the simulation results. CDM supervised and contributed to the design of the study and critically revised the manuscript. All authors read and approved the final manuscript.


This work was supported by funding from The Advanced Research Projects Agnecy-Energy (ARPA-E), U S Department of Energy to Drs. Costas D. Maranas, Thomas K. Wood, and James G. Ferry, Grant DE-AR0000431. The authors would like to thank Anupam Chowdhury for his useful inputs during the preparation of the manuscript.

Competing interests

The authors declare that they have no competing interests.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Costas D. Maranas.

Additional files

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Nazem-Bokaee, H., Gopalakrishnan, S., Ferry, J.G. et al. Assessing methanotrophy and carbon fixation for biofuel production by Methanosarcina acetivorans . Microb Cell Fact 15, 10 (2016).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: