Bioenergetics of the formyl-methanofuran dehydrogenase and heterodisulfide reductase reactions in Methanothermobacter thermautotrophicus

An integrated systems biology approach reveals differences in formate metabolism in the genus Methanothermobacter

Methanogenesis allows methanogenic archaea to generate cellular energy for their growth while producing methane. Thermophilic hydrogenotrophic species of the genus Methanothermobacter have been recognized as robust biocatalysts for a circular carbon economy and are already applied in power-to-gas technology with biomethanation, which is a platform to store renewable energy and utilize captured carbon dioxide. Here, we generated curated genome-scale metabolic reconstructions for three Methanothermobacter strains and investigated differences in the growth performance of these same strains in chemostat bioreactor experiments with hydrogen and carbon dioxide or formate as substrates. Using an integrated systems biology approach, we identified differences in formate anabolism between the strains and revealed that formate anabolism influences the diversion of carbon between biomass and methane. This finding, together with the omics datasets and the metabolic models we generated, can be implemented for biotechnological applications of Methanothermobacter in power-to-gas technology, and as a perspective, for value-added chemical production.

Limited Mechanistic Link Between the Monod Equation and Methanogen Growth: a Perspective from Metabolic Modeling

The Monod equation has been widely applied as the general rate law of microbial growth, but its applications are not always successful. By drawing on the frameworks of kinetic and stoichiometric metabolic models and metabolic control analysis, the modeling reported here simulated the growth kinetics of a methanogenic microorganism and illustrated that different enzymes and metabolites control growth rate to various extents and that their controls peak at either very low, intermediate, or very high substrate concentrations. In comparison, with a single term and two parameters, the Monod equation only approximately accounts for the controls of rate-determining enzymes and metabolites at very high and very low substrate concentrations, but neglects the enzymes and metabolites whose controls are most notable at intermediate concentrations. These findings support a limited link between the Monod equation and methanogen growth, and unify the competing views regarding enzyme roles in shaping growth kinetics. The results also preclude a mechanistic derivation of the Monod equation from methanogen metabolic networks and highlight a fundamental challenge in microbiology: single-term expressions may not be sufficient for accurate prediction of microbial growth. IMPORTANCE The Monod equation has been widely applied to predict the rate of microbial growth, but its application is not always successful. Using a novel metabolic modeling approach, we simulated the growth of a methanogen and uncovered a limited mechanistic link between the Monod equation and the methanogen's metabolic network. Specifically, the equation provides an approximation to the controls by rate-determining metabolites and enzymes at very low and very high substrate concentrations, but it is missing the remaining enzymes and metabolites whose controls are most notable at intermediate concentrations. These results support the Monod equation as a useful approximation of growth rates and highlight a fundamental challenge in microbial kinetics: single-term rate expressions may not be sufficient for accurate prediction of microbial growth.

Controls on the isotopic composition of microbial methane

Microbial methane production (methanogenesis) is responsible for more than half of the annual emissions of this major greenhouse gas to the atmosphere. Although the stable isotopic composition of methane is often used to characterize its sources and sinks, strictly empirical descriptions of the isotopic signature of methanogenesis currently limit these attempts. We developed a metabolic-isotopic model of methanogenesis by carbon dioxide reduction, which predicts carbon and hydrogen isotopic fractionations, and clumped isotopologue distributions, as functions of the cell's environment. We mechanistically explain multiple isotopic patterns in laboratory and natural settings and show that these patterns constrain the in situ energetics of methanogenesis. Combining our model with data from environments in which methanogenic activity is energy-limited, we provide predictions for the biomass-specific methanogenesis rates and the associated isotopic effects.

Sulfide and transition metals - A partnership for life

Sulfide and transition metals often came together in Biology. The variety of possible structural combinations enabled living organisms to evolve an array of highly versatile metal-sulfide centers to fulfill different physiological roles. The ubiquitous iron‑sulfur centers, with their structural, redox, and functional diversity, are certainly the best-known partners, but other metal-sulfide centers, involving copper, nickel, molybdenum or tungsten, are equally crucial for Life. This review provides a concise overview of the exclusive sulfide properties as a metal ligand, with emphasis on the structural aspects and biosynthesis. Sulfide as catalyst and as a substrate is discussed. Different enzymes are considered, including xanthine oxidase, formate dehydrogenases, nitrogenases and carbon monoxide dehydrogenases. The sulfide effect on the activity and function of iron‑sulfur, heme and zinc proteins is also addressed.

Energy Conservation and Hydrogenase Function in Methanogenic Archaea, in Particular the Genus Methanosarcina

The biological production of methane is vital to the global carbon cycle and accounts for ca. 74% of total methane emissions. The organisms that facilitate this process, methanogenic archaea, belong to a large and phylogenetically diverse group that thrives in a wide range of anaerobic environments. Two main subgroups exist within methanogenic archaea: those with and those without cytochromes. Although a variety of metabolisms exist within this group, the reduction of growth substrates to methane using electrons from molecular hydrogen is, in a phylogenetic sense, the most widespread methanogenic pathway. Methanogens without cytochromes typically generate methane by the reduction of CO₂ with electrons derived from H₂, formate, or secondary alcohols, generating a transmembrane ion gradient for ATP production via an Na⁺-translocating methyltransferase (Mtr). These organisms also conserve energy with a novel flavin-based electron bifurcation mechanism, wherein the endergonic reduction of ferredoxin is facilitated by the exergonic reduction of a disulfide terminal electron acceptor coupled to either H₂ or formate oxidation. Methanogens that utilize cytochromes have a broader substrate range, and can convert acetate and methylated compounds to methane, in addition to the ability to reduce CO₂ Cytochrome-containing methanogens are able to supplement the ion motive force generated by Mtr with an H⁺-translocating electron transport system. In both groups, enzymes known as hydrogenases, which reversibly interconvert protons and electrons to molecular hydrogen, play a central role in the methanogenic process. This review discusses recent insight into methanogen metabolism and energy conservation mechanisms with a particular focus on the genus Methanosarcina.

Methanothermobacter thermautotrophicus modulates its membrane lipids in response to hydrogen and nutrient availability

Methanothermobacter thermautotrophicus strain ΔH is a model hydrogenotrophic methanogen, for which extensive biochemical information, including the complete genome sequence, is available. Nevertheless, at the cell membrane lipid level, little is known about the responses of this archaeon to environmental stimuli. In this study, the lipid composition of M. thermautotrophicus was characterized to verify how this archaeon modulates its cell membrane components during growth phases and in response to hydrogen depletion and nutrient limitation (potassium and phosphate). As opposed to the higher abundance of phospholipids in the stationary phase of control experiments, cell membranes under nutrient, and energy stress were dominated by glycolipids that likely provided a more effective barrier against ion leakage. We also identified particular lipid regulatory mechanisms in M. thermautotrophicus, which included the accumulation of polyprenols under hydrogen-limited conditions and an increased content of sodiated adducts of lipids in nutrient-limited cells. These findings suggest that M. thermautotrophicus intensely modulates its cell membrane lipid composition to cope with energy and nutrient availability in dynamic environments.

Methanogenic archaea: ecologically relevant differences in energy conservation

Most methanogenic archaea can reduce CO(2) with H(2) to methane, and it is generally assumed that the reactions and mechanisms of energy conservation that are involved are largely the same in all methanogens. However, this does not take into account the fact that methanogens with cytochromes have considerably higher growth yields and threshold concentrations for H(2) than methanogens without cytochromes. These and other differences can be explained by the proposal outlined in this Review that in methanogens with cytochromes, the first and last steps in methanogenesis from CO(2) are coupled chemiosmotically, whereas in methanogens without cytochromes, these steps are energetically coupled by a cytoplasmic enzyme complex that mediates flavin-based electron bifurcation.

The past and present of sodium energetics: may the sodium-motive force be with you

All living cells routinely expel Na(+) ions, maintaining lower concentration of Na(+) in the cytoplasm than in the surrounding milieu. In the vast majority of bacteria, as well as in mitochondria and chloroplasts, export of Na(+) occurs at the expense of the proton-motive force. Some bacteria, however, possess primary generators of the transmembrane electrochemical gradient of Na(+) (sodium-motive force). These primary Na(+) pumps have been traditionally seen as adaptations to high external pH or to high temperature. Subsequent studies revealed, however, the mechanisms for primary sodium pumping in a variety of non-extremophiles, such as marine bacteria and certain bacterial pathogens. Further, many alkaliphiles and hyperthermophiles were shown to rely on H(+), not Na(+), as the coupling ion. We review here the recent progress in understanding the role of sodium-motive force, including (i) the conclusion on evolutionary primacy of the sodium-motive force as energy intermediate, (ii) the mechanisms, evolutionary advantages and limitations of switching from Na(+) to H(+) as the coupling ion, and (iii) the possible reasons why certain pathogenic bacteria still rely on the sodium-motive force.

Methyl sulfide production by a novel carbon monoxide metabolism in Methanosarcina acetivorans

We observed dimethyl sulfide and methanthiol production in pure incubations of the methanogen Methanosarcina acetivorans when carbon monoxide (CO) served as the only electron donor. Energy conservation likely uses sodium ion gradients for ATP synthesis. This novel metabolism permits utilization of CO by the methanogen, resulting in quantitative sulfide methylation.

On the origin of biochemistry at an alkaline hydrothermal vent

A model for the origin of biochemistry at an alkaline hydrothermal vent has been developed that focuses on the acetyl-CoA (Wood-Ljungdahl) pathway of CO2 fixation and central intermediary metabolism leading to the synthesis of the constituents of purines and pyrimidines. The idea that acetogenesis and methanogenesis were the ancestral forms of energy metabolism among the first free-living eubacteria and archaebacteria, respectively, stands in the foreground. The synthesis of formyl pterins, which are essential intermediates of the Wood-Ljungdahl pathway and purine biosynthesis, is found to confront early metabolic systems with steep bioenergetic demands that would appear to link some, but not all, steps of CO2 reduction to geochemical processes in or on the Earth's crust. Inorganically catalysed prebiotic analogues of the core biochemical reactions involved in pterin-dependent methyl synthesis of the modern acetyl-CoA pathway are considered. The following compounds appear as probable candidates for central involvement in prebiotic chemistry: metal sulphides, formate, carbon monoxide, methyl sulphide, acetate, formyl phosphate, carboxy phosphate, carbamate, carbamoyl phosphate, acetyl thioesters, acetyl phosphate, possibly carbonyl sulphide and eventually pterins. Carbon might have entered early metabolism via reactions hardly different from those in the modern Wood-Ljungdahl pathway, the pyruvate synthase reaction and the incomplete reverse citric acid cycle. The key energy-rich intermediates were perhaps acetyl thioesters, with acetyl phosphate possibly serving as the universal metabolic energy currency prior to the origin of genes. Nitrogen might have entered metabolism as geochemical NH3 via two routes: the synthesis of carbamoyl phosphate and reductive transaminations of alpha-keto acids. Together with intermediates of methyl synthesis, these two routes of nitrogen assimilation would directly supply all intermediates of modern purine and pyrimidine biosynthesis. Thermodynamic considerations related to formyl pterin synthesis suggest that the ability to harness a naturally pre-existing proton gradient at the vent-ocean interface via an ATPase is older than the ability to generate a proton gradient with chemistry that is specified by genes.

Adaptations to energy stress dictate the ecology and evolution of the Archaea

The three domains of life on Earth include the two prokaryotic groups, Archaea and Bacteria. The Archaea are distinguished from Bacteriabased on phylogenetic and biochemical differences, but currently there is no unifying ecological principle to differentiate these groups. The ecology of the Archaea is reviewed here in terms of cellular bioenergetics. Adaptation to chronic energy stress is hypothesized to be the crucial factor that distinguishes the Archaea from Bacteria. The biochemical mechanisms that enable archaea to cope with chronic energy stress include low-permeability membranes and specific catabolic pathways. Based on the ecological unity and biochemical adaptations among archaea, I propose the hypothesis that chronic energy stress is the primary selective pressure governing the evolution of the Archaea.

Isolation and characterization of an amiloride-resistant mutant of Methanothermobacter thermautotrophicus possessing a defective Na+/H+ antiport

A spontaneous mutant of Methanothermobacter thermautotrophicus resistant to the Na+/H+ antiporter inhibitor amiloride was isolated. The Na+/H+ exchanger activity in the mutant cells was remarkably decreased in comparison with wild-type cells. Methanogenesis rates in the mutant strain were higher than wild-type cells and resistant to the inhibitory effect of 2 mM amiloride. In contrast, methanogenesis in wild-type cells was completely inhibited by the same amiloride concentration. ATP synthesis driven by methanogenic electron transport or by an electrogenic potassium efflux in the presence of sodium ions was significantly enhanced in the mutant cells. ATP synthesis driven by potassium diffusion potential was profoundly inhibited in wild-type cells by the presence of uncoupler 3,3',4',5- tetrachlorosalicylanilide and sodium ions, whereas c. 50% inhibition was observed in the mutant cells under the same conditions.

Coupling of Methanothermobacter thermautotrophicus methane formation and growth in fed-batch and continuous cultures under different H2 gassing regimens

In nature, H2- and CO2-utilizing methanogenic archaea have to couple the processes of methanogenesis and autotrophic growth under highly variable conditions with respect to the supply and concentration of their energy source, hydrogen. To study the hydrogen-dependent coupling between methanogenesis and growth, Methanothermobacter thermautotrophicus was cultured in a fed-batch fermentor and in a chemostat under different 80% H(2)-20% CO2 gassing regimens while we continuously monitored the dissolved hydrogen partial pressures (pH2). In the fed-batch system, in which the conditions continuously changed the uptake rates by the growing biomass, the organism displayed a complex and yet defined growth behavior, comprising the consecutive lag, exponential, and linear growth phases. It was found that the in situ hydrogen concentration affected the coupling between methanogenesis and growth in at least two respects. (i) The microorganism could adopt two distinct theoretical maximal growth yields (YCH4 max), notably approximately 3 and 7 g (dry weight) of methane formed mol-1, for growth under low (pH2 < 12 kPa)- and high-hydrogen conditions, respectively. The distinct values can be understood from a theoretical analysis of the process of methanogenesis presented in the supplemental material associated with this study. (ii) The in situ hydrogen concentration affected the "specific maintenance" requirements or, more likely, the degree of proton leakage and proton slippage processes. At low pH2 values, the "specific maintenance" diminished and the specific growth yields approached YCH4 max, indicating that growth and methanogenesis became fully coupled.

Modeling methanogenesis with a genome-scale metabolic reconstruction of Methanosarcina barkeri

We present a genome-scale metabolic model for the archaeal methanogen Methanosarcina barkeri. We characterize the metabolic network and compare it to reconstructions from the prokaryotic, eukaryotic and archaeal domains. Using the model in conjunction with constraint-based methods, we simulate the metabolic fluxes and resulting phenotypes induced by different environmental and genetic conditions. This represents the first large-scale simulation of either a methanogen or an archaeal species. Model predictions are validated by comparison to experimental growth measurements and phenotypes of M. barkeri on different substrates. The predicted growth phenotypes for wild type and mutants of the methanogenic pathway have a high level of agreement with experimental findings. We further examine the efficiency of the energy-conserving reactions in the methanogenic pathway, specifically the Ech hydrogenase reaction, and determine a stoichiometry for the nitrogenase reaction. This work demonstrates that a reconstructed metabolic network can serve as an analysis platform to predict cellular phenotypes, characterize methanogenic growth, improve the genome annotation and further uncover the metabolic characteristics of methanogenesis.

Protein complexes in the archaeon Methanothermobacter thermautotrophicus analyzed by blue native/SDS-PAGE and mass spectrometry

Methanothermobacter thermautotrophicus is a thermophilic archaeon that produces methane as the end product of its primary metabolism. The biochemistry of methane formation has been extensively studied and is catalyzed by individual enzymes and proteins that are organized in protein complexes. Although much is known of the protein complexes involved in methanogenesis, only limited information is available on the associations of proteins involved in other cell processes of M. thermautotrophicus. To visualize and identify interacting and individual proteins of M. thermautotrophicus on a proteome-wide scale, protein preparations were separated using blue native electrophoresis followed by SDS-PAGE. A total of 361 proteins, corresponding to almost 20% of the predicted proteome, was identified using peptide mass fingerprinting after MALDI-TOF MS. All previously characterized complexes involved in energy generation could be visualized. Furthermore the expression and association of the heterodisulfide reductase and methylviologen-reducing hydrogenase complexes depended on culture conditions. Also homomeric supercomplexes of the ATP synthase stalk subcomplex and the N5-methyl-5,6,7,8-tetrahydromethanopterin:coenzyme M methyltransferase complex were separated. Chemical cross-linking experiments confirmed that the multimerization of both complexes was not experimentally induced. A considerable number of previously uncharacterized protein complexes were reproducibly visualized. These included an exosome-like complex consisting of four exosome core subunits, which associated with a tRNA-intron endonuclease, thereby expanding the constituency of archaeal exosomes. The results presented show the presence of novel complexes and demonstrate the added value of including blue native gel electrophoresis followed by SDS-PAGE in discovering protein complexes that are involved in catabolic, anabolic, and general cell processes.

Hydrogen concentrations in methane-forming cells probed by the ratios of reduced and oxidized coenzyme F420

Coenzyme F420 is the central low-redox-potential electron carrier in methanogenic metabolism. The coenzyme is reduced under hydrogen by the action of F420-dependent hydrogenase. The standard free-energy change at pH 7 of F420 reduction was determined to be -15 kJ mol(-1), irrespective of the temperature (25-65 degrees C). Experiments performed with methane-forming cell suspensions of Methanothermobacter thermautotrophicus incubated under various conditions demonstrated that the ratios of reduced and oxidized F420 were in thermodynamic equilibrium with the gas-phase hydrogen partial pressures. During growth in a fed-batch fermenter, ratios changed in connection with the decrease in dissolved hydrogen. For most of the time, the changes were as expected for thermodynamic equilibrium between the oxidation state of F420 inside the cells and extracellular hydrogen. Also, methanol-metabolizing, but not acetate-converting, cells of Methanosarcina barkeri maintained the ratios of reduced and oxidized coenzyme F420 in thermodynamic equilibrium with external hydrogen. The results of the study demonstrate that F420 is a useful probe to assess in situ hydrogen concentrations in H2-metabolizing methanogens.

The membrane-bound electron transport system of Methanosarcina species

Members of the genus Methanosarcina are strictly anaerobic archaea that derive their metabolic energy from the conversion of a restricted number of substrates to methane. H2 + CO2 are converted to CH4 via the CO2-reducing pathway, while methanol and methylamines are metabolized by the methylotrophic pathway. Two novel electron transport systems are involved in the process of methanogenesis. Both systems are able to use a heterodisulfide as electron acceptor and either H2 or F420H2 as electron acceptors and generate a proton-motive force by redox potential-driven H(+)-translocation. The H2:heterodisulfide oxidoreductase is composed of an F420-nonreducing hydrogenase and the heterodisulfide reductase. The latter protein is also part of the F420H2:heterodisulfide oxidoreductase system. The second component of this system is referred to as F420H2 dehydrogenase. The archaeal protein is a homologue of complex I of the respiratory chain from bacteria and mitochondria. This review focuses on the biochemical and genetic characteristics of the three energy-transducing enzymes and on the mechanisms of ion translocation.

Bioenergetics of the formyl-methanofuran dehydrogenase and heterodisulfide reductase reactions in Methanothermobacter thermautotrophicus

The synthesis of formyl-methanofuran and the reduction of the heterodisulfide (CoM-S-S-CoB) of coenzyme M (HS-CoM) and coenzyme B (HS-CoB) are two crucial, H2-dependent reactions in the energy metabolism of methanogenic archaea. The bioenergetics of the reactions in vivo were studied in chemostat cultures and in cell suspensions of Methanothermobacter thermautotrophicus metabolizing at defined dissolved hydrogen partial pressures ( pH2). Formyl-methanofuran synthesis is an endergonic reaction (DeltaG degrees ' = +16 kJ.mol-1). By analyzing the concentration ratios between formyl-methanofuran and methanofuran in the cells, free energy changes under experimental conditions (DeltaG') were found to range between +10 and +35 kJ.mol-1 depending on the pH2 applied. The comparison with the sodium motive force indicated that the reaction should be driven by the import of a variable number of two to four sodium ions. Heterodisulfide reduction (DeltaG degrees ' = -40 kJ.mol-1) was associated with free energy changes as high as -55 to -80 kJ.mol-1. The values were determined by analyzing the concentrations of CoM-S-S-CoB, HS-CoM and HS-CoB in methane-forming cells operating under a variety of hydrogen partial pressures. Free energy changes were in equilibrium with the proton motive force to the extent that three to four protons could be translocated out of the cells per reaction. Remarkably, an apparent proton translocation stoichiometry of three held for cells that had been grown at pH2<0.12 bar, whilst the number was four for cells grown above that concentration. The shift occurred within a narrow pH2 span around 0.12 bar. The findings suggest that the methanogens regulate the bioenergetic machinery involved in CoM-S-S-CoB reduction and proton pumping in response to the environmental hydrogen concentrations.