Functional diversity of prokaryotic HdrA(BC) modules: Role in flavin-based electron bifurcation processes and beyond

Fe/S proteins in microbial sulfur oxidation

Iron-sulfur clusters serve as indispensable cofactors within proteins across all three domains of life. Fe/S clusters emerged early during the evolution of life on our planet and the biogeochemical cycle of sulfur is one of the most ancient and important element cycles. It is therefore no surprise that Fe/S proteins have crucial roles in the multiple steps of microbial sulfur metabolism. During dissimilatory sulfur oxidation in prokaryotes, Fe/S proteins not only serve as electron carriers in several steps, but also perform catalytic roles, including unprecedented reactions. Two cytoplasmic enzyme systems that oxidize sulfane sulfur to sulfite are of particular interest in this context: The rDsr pathway employs the reverse acting dissimilatory sulfite reductase rDsrAB as its key enzyme, while the sHdr pathway utilizes polypeptides resembling the HdrA, HdrB and HdrC subunits of heterodisulfide reductase from methanogenic archaea. Both pathways involve components predicted to bind unusual noncubane Fe/S clusters acting as catalysts for the formation of disulfide or sulfite. Mapping of Fe/S cluster machineries on the sulfur-oxidizing prokaryote tree reveals that ISC, SUF, MIS and SMS are all sufficient to meet the Fe/S cluster maturation requirements for operation of the sHdr or rDsr pathways. The sHdr pathway is dependent on lipoate-binding proteins that are assembled by a novel pathway, involving two Radical SAM proteins, namely LipS1 and LipS2. These proteins coordinate sulfur-donating auxiliary Fe/S clusters in atypical patterns by three cysteines and one histidine and act as lipoyl synthases by jointly inserting two sulfur atoms to an octanoyl residue. This article is part of a Special Issue entitled: Biogenesis and Function of Fe/S proteins.

Fusion/fission protein family identification in Archaea

The majority of newly discovered archaeal lineages remain without a cultivated representative, but scarce experimental data from the cultivated organisms show that they harbor distinct functional repertoires. To unveil the ecological as well as evolutionary impact of Archaea from metagenomics, new computational methods need to be developed, followed by in-depth analysis. Among them is the genome-wide protein fusion screening performed here. Natural fusions and fissions of genes not only contribute to microbial evolution but also complicate the correct identification and functional annotation of sequences. The products of these processes can be defined as fusion (or composite) proteins, the ones consisting of two or more domains originally encoded by different genes and split proteins, and the ones originating from the separation of a gene in two (fission). Fusion identifications are required for proper phylogenetic reconstructions and metabolic pathway completeness assessments, while mappings between fused and unfused proteins can fill some of the existing gaps in metabolic models. In the archaeal genome-wide screening, more than 1,900 fusion/fission protein clusters were identified, belonging to both newly sequenced and well-studied lineages. These protein families are mainly associated with different types of metabolism, genetic, and cellular processes. Moreover, 162 of the identified fusion/fission protein families are archaeal specific, having no identified fused homolog within the bacterial domain. Our approach was validated by the identification of experimentally characterized fusion/fission cases. However, around 25% of the identified fusion/fission families lack functional annotations for both composite and split states, showing the need for experimental characterization in Archaea.IMPORTANCEGenome-wide fusion screening has never been performed in Archaea on a broad taxonomic scale. The overlay of multiple computational techniques allows the detection of a fine-grained set of predicted fusion/fission families, instead of rough estimations based on conserved domain annotations only. The exhaustive mapping of fused proteins to bacterial organisms allows us to capture fusion/fission families that are specific to archaeal biology, as well as to identify links between bacterial and archaeal lineages based on cooccurrence of taxonomically restricted proteins and their sequence features. Furthermore, the identification of poorly characterized lineage-specific fusion proteins opens up possibilities for future experimental and computational investigations. This approach enhances our understanding of Archaea in general and provides potential candidates for in-depth studies in the future.

Genomic comparison of deep-sea hydrothermal genera related to Aeropyrum, Thermodiscus and Caldisphaera, and proposed emended description of the family Acidilobaceae

Deep-sea hydrothermal vents host archaeal and bacterial thermophilic communities, including taxonomically and functionally diverse Thermoproteota. Despite their prevalence in high-temperature submarine communities, Thermoproteota are chronically under-represented in genomic databases and issues have emerged regarding their nomenclature, particularly within the Aeropyrum-Thermodiscus-Caldisphaera. To resolve some of these problems, we identified 47 metagenome-assembled genomes (MAGs) within this clade, from 20 previously published deep-sea hydrothermal vent and submarine volcano metagenomes, and 24 MAGs from public databases. Using phylogenomic analysis, Genome Taxonomy Database Toolkit (GTDB-Tk) taxonomic assessment, 16S rRNA gene phylogeny, average amino acid identity (AAI) and functional gene patterns, we re-evaluated of the taxonomy of the Aeropyrum-Thermodiscus-Caldisphaera. At least nine genus-level clades were identified with two or more MAGs. In accordance with SeqCode requirements and recommendations, we propose names for three novel genera, viz. Tiamatella incendiivivens, Hestiella acidicharens and Calypsonella navitae. A fourth genus was also identified related to Thermodiscus maritimus, for which no available sequenced genome exists. We propose the novel species Thermodiscus eudorianus to describe our high-quality Thermodiscus MAG, which represents the type genome for the genus. All three novel genera and T. eudorianus are likely anaerobic heterotrophs, capable of fermenting protein-rich carbon sources, while some Tiamatella, Calypsonella and T. eudorianus may also reduce polysulfides, thiosulfate, sulfur and/or selenite, and the likely acidophile, Hestiella, may reduce nitrate and/or perchlorate. Based on phylogenomic evidence, we also propose the family Acidilobaceae be amended to include Caldisphaera, Aeropyrum, Thermodiscus and Stetteria and the novel genera described here.

Combining metabolic flux analysis with proteomics to shed light on the metabolic flexibility: the case of Desulfovibrio vulgaris Hildenborough

Introduction

Desulfovibrio vulgaris Hildenborough is a gram-negative anaerobic bacterium belonging to the sulfate-reducing bacteria that exhibits highly versatile metabolism. By switching from one energy mode to another depending on nutrients availability in the environments" it plays a central role in shaping ecosystems. Despite intensive efforts to study D. vulgaris energy metabolism at the genomic, biochemical and ecological level, bioenergetics in this microorganism remain far from being fully understood. Alternatively, metabolic modeling is a powerful tool to understand bioenergetics. However, all the current models for D. vulgaris appeared to be not easily adaptable to various environmental conditions.

Methods

To lift off these limitations, here we constructed a novel transparent and robust metabolic model to explain D. vulgaris bioenergetics by combining whole-cell proteomic analysis with modeling approaches (Flux Balance Analysis).

Results

The iDvu71 model showed over 0.95 correlation with experimental data. Further simulations allowed a detailed description of D. vulgaris metabolism in various conditions of growth. Altogether, the simulations run in this study highlighted the sulfate-to-lactate consumption ratio as a pivotal factor in D. vulgaris energy metabolism.

Discussion

In particular, the impact on the hydrogen/formate balance and biomass synthesis is discussed. Overall, this study provides a novel insight into D. vulgaris metabolic flexibility.

Mcr-dependent methanogenesis in Archaeoglobaceae enriched from a terrestrial hot spring

The preeminent source of biological methane on Earth is methyl coenzyme M reductase (Mcr)-dependent archaeal methanogenesis. A growing body of evidence suggests a diversity of archaea possess Mcr, although experimental validation of hypothesized methane metabolisms has been missing. Here, we provide evidence of a functional Mcr-based methanogenesis pathway in a novel member of the family Archaeoglobaceae, designated Methanoglobus nevadensis, which we enriched from a terrestrial hot spring on the polysaccharide xyloglucan. Our incubation assays demonstrate methane production that is highly sensitive to the Mcr inhibitor bromoethanesulfonate, stimulated by xyloglucan and xyloglucan-derived sugars, concomitant with the consumption of molecular hydrogen, and causing a deuterium fractionation in methane characteristic of hydrogenotrophic and methylotrophic methanogens. Combined with the recovery and analysis of a high-quality M. nevadensis metagenome-assembled genome encoding a divergent Mcr and diverse potential electron and carbon transfer pathways, our observations suggest methanogenesis in M. nevadensis occurs via Mcr and is fueled by the consumption of cross-fed byproducts of xyloglucan fermentation mediated by other community members. Phylogenetic analysis shows close affiliation of the M. nevadensis Mcr with those from Korarchaeota, Nezhaarchaeota, Verstraetearchaeota, and other Archaeoglobales that are divergent from well-characterized Mcr. We propose these archaea likely also use functional Mcr complexes to generate methane on the basis of our experimental validation in M. nevadensis. Thus, divergent Mcr-encoding archaea may be underestimated sources of biological methane in terrestrial and marine hydrothermal environments.

Physiological potential and evolutionary trajectories of syntrophic sulfate-reducing bacterial partners of anaerobic methanotrophic archaea

Sulfate-coupled anaerobic oxidation of methane (AOM) is performed by multicellular consortia of anaerobic methanotrophic archaea (ANME) in obligate syntrophic partnership with sulfate-reducing bacteria (SRB). Diverse ANME and SRB clades co-associate but the physiological basis for their adaptation and diversification is not well understood. In this work, we used comparative metagenomics and phylogenetics to investigate the metabolic adaptation among the 4 main syntrophic SRB clades (HotSeep-1, Seep-SRB2, Seep-SRB1a, and Seep-SRB1g) and identified features associated with their syntrophic lifestyle that distinguish them from their non-syntrophic evolutionary neighbors in the phylum Desulfobacterota. We show that the protein complexes involved in direct interspecies electron transfer (DIET) from ANME to the SRB outer membrane are conserved between the syntrophic lineages. In contrast, the proteins involved in electron transfer within the SRB inner membrane differ between clades, indicative of convergent evolution in the adaptation to a syntrophic lifestyle. Our analysis suggests that in most cases, this adaptation likely occurred after the acquisition of the DIET complexes in an ancestral clade and involve horizontal gene transfers within pathways for electron transfer (CbcBA) and biofilm formation (Pel). We also provide evidence for unique adaptations within syntrophic SRB clades, which vary depending on the archaeal partner. Among the most widespread syntrophic SRB, Seep-SRB1a, subclades that specifically partner ANME-2a are missing the cobalamin synthesis pathway, suggestive of nutritional dependency on its partner, while closely related Seep-SRB1a partners of ANME-2c lack nutritional auxotrophies. Our work provides insight into the features associated with DIET-based syntrophy and the adaptation of SRB towards it.

DsrC is involved in fermentative growth and interacts directly with the FlxABCD-HdrABC complex in Desulfovibrio vulgaris Hildenborough

DsrC is a key protein in dissimilatory sulfur metabolism, where it works as co-substrate of the dissimilatory sulfite reductase DsrAB. DsrC has two conserved cysteines in a C-terminal arm that are converted to a trisulfide upon reduction of sulfite. In sulfate-reducing bacteria, DsrC is essential and previous works suggested additional functions beyond sulfite reduction. Here, we studied whether DsrC also plays a role during fermentative growth of Desulfovibrio vulgaris Hildenborough, by studying two strains where the functionality of DsrC is impaired by a lower level of expression (IPFG07) and additionally by the absence of one conserved Cys (IPFG09). Growth studies coupled with metabolite and proteomic analyses reveal that fermentation leads to lower levels of DsrC, but impairment of its function results in reduced growth by fermentation and a shift towards more fermentative metabolism during sulfate respiration. In both respiratory and fermentative conditions, there is increased abundance of the FlxABCD-HdrABC complex and Adh alcohol dehydrogenase in IPFG09 versus the wild type, which is reflected in higher production of ethanol. Pull-down experiments confirmed a direct interaction between DsrC and the FlxABCD-HdrABC complex, through the HdrB subunit. Dissimilatory sulfur metabolism, where sulfur compounds are used for energy generation, is a key process in the ecology of anoxic environments, and is more widespread among bacteria than previously believed. Two central proteins for this type of metabolism are DsrAB dissimilatory sulfite reductase and its co-substrate DsrC. Using physiological, proteomic and biochemical studies of Desulfovibrio vulgaris Hildenborough and mutants affected in DsrC functionality, we show that DsrC is also relevant for fermentative growth of this model organism and that it interacts directly with the soluble FlxABCD-HdrABC complex that links the NAD(H) pool with dissimilatory sulfite reduction.

Effects of phenyl acids on different degradation phases during thermophilic anaerobic digestion

Aromatic compounds like phenyl acids (PA) can accumulate during anaerobic digestion (AD) of organic wastes due to an increased entry of lignocellulose, secondary plant metabolites or proteins, and thermodynamic challenges in degrading the benzene ring. The effects of aromatic compounds can be various – from being highly toxic to be stimulating for methanogenesis – depending on many parameters like inoculum or molecular characteristics of the aromatic compound. To contribute to a better understanding of the consequences of PA exposure during AD, the aim was to evaluate the effects of 10 mM PA on microbial communities degrading different, degradation phase–specific substrates in thermophilic batch reactors within 28  days: Microcrystalline cellulose (MCC, promoting hydrolytic to methanogenic microorganisms), butyrate or propionate (promoting syntrophic volatile fatty acid (VFA) oxidisers to methanogens), or acetate (promoting syntrophic acetate oxidisers to methanogens). Methane production, VFA concentrations and pH were evaluated, and microbial communities and extracellular polymeric substances (EPS) were assessed. The toxicity of PA depended on the type of substrate which in turn determined the (i) microbial diversity and composition and (ii) EPS quantity and quality. Compared with the respective controls, methane production in MCC reactors was less impaired by PA than in butyrate, propionate and acetate reactors which showed reductions in methane production of up to 93%. In contrast to the controls, acetate concentrations were high in all PA reactors at the end of incubation thus acetate was a bottle-neck intermediate in those reactors. Considerable differences in EPS quantity and quality could be found among substrates but not among PA variants of each substrate. Methanosarcina spp. was the dominant methanogen in VFA reactors without PA exposure and was inhibited when PA were present. VFA oxidisers and Methanothermobacter spp. were abundant in VFA assays with PA exposure as well as in all MCC reactors. As MCC assays showed higher methane yields, a higher microbial diversity and a higher EPS quantity and quality than VFA reactors when exposed to PA, we conclude that EPS in MCC reactors might have been beneficial for absorbing/neutralising phenyl acids and keeping (more susceptible) microorganisms shielded in granules or biofilms.

Substrate-Dependent Conformational Switch of the Noncubane [4Fe-4S] Cluster in Heterodisulfide Reductase HdrB

The noncubane [4Fe-4S] cluster identified in the active site of heterodisulfide reductase (HdrB) displays a unique geometry among Fe-S cofactors found in metalloproteins. Here we employ resonance Raman (RR) spectroscopy and density functional theory (DFT) calculations to probe structural, electronic, and vibrational properties of the noncubane cluster in HdrB from a non-methanogenic Desulfovibrio vulgaris (Dv) Hildenborough organism. The immediate protein environment of the two neighboring clusters in DvHdrB is predicted using homology modeling. We demonstrate that in the absence of substrate, the oxidized [4Fe-4S]³⁺ cluster adopts a "closed" conformation. Upon substrate coordination at the "special" iron center, the cluster core translates to an "open" structure, facilitated by the "supernumerary" cysteine ligand switch from iron-bridging to iron-terminal mode. The observed RR fingerprint of the noncubane cluster, supported by Fe-S vibrational mode analysis, will advance future studies of enzymes containing this unusual cofactor.

A metabolic puzzle: Consumption of C1 compounds and thiosulfate in Hyphomicrobium denitrificans XT

Many obligately heterotrophic methylotrophs oxidize thiosulfate as an additional electron source during growth on C₁ compounds. Although two different pathways of thiosulfate oxidation are implemented in Hyphomicrobium denitrificans X^T, a pronounced negative effect on growth rate is observed when it is cultured in the simultaneous presence of methanol and thiosulfate. In this model organism, periplasmic thiosulfate dehydrogenase TsdA catalyzes formation of the dead-end product tetrathionate. By reverse genetics we verified the second pathway that also starts in the periplasm where SoxXA catalyzes the oxidative fusion of thiosulfate to SoxYZ, from which sulfate is released by SoxB. Sulfane sulfur is then further oxidized in the cytoplasm by the sulfur-oxidizing heterodisulfide reductase-like system (sHdr) which is produced constitutively in a strain lacking the transcriptional repressor sHdrR. When exposed to thiosulfate, the ΔshdrR strain exhibited a strongly reduced growth rate even without thiosulfate in the pre-cultures. When grown on methanol, cells exhibit significantly increased NAD⁺/NADH ratios in the presence of thiosulfate. In contrast, thiosulfate did not exert any negative effect on growth rate or increase NAD⁺ levels during growth on formate. On both C₁ substrates, excretion of up to 0.5 mM sulfite as an intermediate of thiosulfate (2 mM) oxidation was recorded. Sulfite is known to form adducts with pyrroloquinoline quinone, the cofactor of periplasmic methanol dehydrogenase. We rationalize that this causes specific inhibition of methanol degradation in the presence of thiosulfate while formate metabolism in the cytoplasm remains unaffected.

Site-Differentiated Iron–Sulfur Cluster Ligation Affects Flavin-Based Electron Bifurcation Activity

Electron bifurcation is an elegant mechanism of biological energy conversion that effectively couples three different physiologically relevant substrates. As such, enzymes that perform this function often play critical roles in modulating cellular redox metabolism. One such enzyme is NADH-dependent reduced-ferredoxin: NADP⁺ oxidoreductase (NfnSL), which couples the thermodynamically favorable reduction of NAD⁺ to drive the unfavorable reduction of ferredoxin from NADPH. The interaction of NfnSL with its substrates is constrained to strict stoichiometric conditions, which ensures minimal energy losses from non-productive intramolecular electron transfer reactions. However, the determinants for this are not well understood. One curious feature of NfnSL is that both initial acceptors of bifurcated electrons are unique iron–sulfur (FeS) clusters containing one non-cysteinyl ligand each. The biochemical impact and mechanistic roles of site-differentiated FeS ligands are enigmatic, despite their incidence in many redox active enzymes. Herein, we describe the biochemical study of wild-type NfnSL and a variant in which one of the site-differentiated ligands has been replaced with a cysteine. Results of dye-based steady-state kinetics experiments, substrate-binding measurements, biochemical activity assays, and assessments of electron distribution across the enzyme indicate that this site-differentiated ligand in NfnSL plays a role in maintaining fidelity of the coordinated reactions performed by the two electron transfer pathways. Given the commonality of these cofactors, our findings have broad implications beyond electron bifurcation and mechanistic biochemistry and may inform on means of modulating the redox balance of the cell for targeted metabolic engineering approaches.

An Abundant and Diverse New Family of Electron Bifurcating Enzymes With a Non-canonical Catalytic Mechanism

Microorganisms utilize electron bifurcating enzymes in metabolic pathways to carry out thermodynamically unfavorable reactions. Bifurcating FeFe-hydrogenases (HydABC) reversibly oxidize NADH (E′∼−280 mV, under physiological conditions) and reduce protons to H2 gas (E°′−414 mV) by coupling this endergonic reaction to the exergonic reduction of protons by reduced ferredoxin (Fd) (E′∼−500 mV). We show here that HydABC homologs are surprisingly ubiquitous in the microbial world and are represented by 57 phylogenetically distinct clades but only about half are FeFe-hydrogenases. The others have replaced the hydrogenase domain with another oxidoreductase domain or they contain additional subunits, both of which enable various third reactions to be reversibly coupled to NAD+ and Fd reduction. We hypothesize that all of these enzymes carry out electron bifurcation and that their third substrates can include hydrogen peroxide, pyruvate, carbon monoxide, aldehydes, aryl-CoA thioesters, NADP+, cofactor F420, formate, and quinones, as well as many yet to be discovered. Some of the enzymes are proposed to be integral membrane-bound proton-translocating complexes. These different functionalities are associated with phylogenetically distinct clades and in many cases with specific microbial phyla. We propose that this new and abundant class of electron bifurcating enzyme be referred to as the Bfu family whose defining feature is a conserved bifurcating BfuBC core. This core contains FMN and six iron sulfur clusters and it interacts directly with ferredoxin (Fd) and NAD(H). Electrons to or from the third substrate are fed into the BfuBC core via BfuA. The other three known families of electron bifurcating enzyme (abbreviated as Nfn, EtfAB, and HdrA) contain a special FAD that bifurcates electrons to high and low potential pathways. The Bfu family are proposed to use a different electron bifurcation mechanism that involves a combination of FMN and three adjacent iron sulfur clusters, including a novel [2Fe-2S] cluster with pentacoordinate and partial non-Cys coordination. The absolute conservation of the redox cofactors of BfuBC in all members of the Bfu enzyme family indicate they have the same non-canonical mechanism to bifurcate electrons. A hypothetical catalytic mechanism is proposed as a basis for future spectroscopic analyses of Bfu family members.

A Third Way of Energy Conservation in Acetogenic Bacteria

Acetogenic bacteria are a group of strictly anaerobic bacteria that make a living from acetate formation from two molecules of CO₂ via the Wood-Ljungdahl pathway (WLP). The free energy change of this reaction is very small and allows the synthesis of only a fraction of an ATP. How this pathway is coupled to energy conservation has been an enigma since its discovery ~90 years ago. Here, we describe an electron transport chain in the cytochrome- and quinone-containing acetogen Sporomusa ovata that leads from molecular hydrogen as an electron donor to an intermediate of the WLP, methylenetetrahydrofolate (methylene-tetrahydrofolate [THF]), as an electron acceptor. The catalytic site of the hydrogenase is periplasmic and likely linked cytochrome b to the membrane. We provide evidence that the MetVF-type methylenetetrahydrofolate reductase is linked proteins MvhD and HdrCBA to the cytoplasmic membrane. Membrane preparations catalyzed the H₂-dependent reduction of methylene-THF to methyl-THF. In our model, a transmembrane electrochemical H⁺ gradient is established by both scalar and vectorial protons that leads to the synthesis of 0.5 mol ATP/mol methylene-THF by a H⁺-F₁F_o ATP synthase. This H₂- and methylene-THF-dependent electron transport chain may be present in other cytochrome-containing acetogens as well and represents a third way of chemiosmotic energy conservation in acetogens, but only in addition to the well-established respiratory enzymes Rnf and Ech.

IMPORTANCE Acetogenic bacteria grow by making acetate from CO₂ and are considered the first life forms on Earth since they couple CO₂ reduction to the conservation of energy. How this is achieved has been an enigma ever since. Recently, two respiratory enzymes, a ferredoxin:NAD⁺ oxidoreductase (Rnf) and a ferredoxin:H⁺ oxidoreductase (Ech), have been found in cytochrome-free acetogenic model bacteria. However, some acetogens contain cytochromes in addition, and there has been a long-standing assumption of a cytochrome-containing electron transport chain in those acetogens. Here, we provide evidence for a respiratory chain in Sporomusa ovata that has a cytochrome-containing hydrogenase as the electron donor and a methylenetetrahydrofolate reductase as the terminal electron acceptor. This is the third way of chemiosmotic energy conservation found in acetogens.

HMS-S-S: a tool for the identification of sulfur metabolism-related genes and analysis of operon structures in genome and metagenome assemblies

Sulfur compounds are used in a variety of biological processes including respiration and photosynthesis. Sulfide and sulfur compounds of intermediary oxidation state can serve as electron donors for lithotrophic growth while sulfate, thiosulfate and sulfur are used as electron acceptors in anaerobic respiration. The biochemistry underlying the manifold transformations of inorganic sulfur compounds occurring in sulfur metabolizing prokaryotes is astonishingly complex and knowledge about it has immensely increased over the last years. The advent of next-generation sequencing approaches as well as the significant increase of data availability in public databases has driven focus of environmental microbiology to probing the metabolic capacity of microbial communities by analysis of this sequence data. To facilitate these analyses, we created HMS-S-S, a comprehensive equivalogous hidden Markov model (HMM)-supported tool. Protein sequences related to sulfur compound oxidation, reduction, transport and intracellular transfer are efficiently detected and related enzymes involved in dissimilatory sulfur oxidation as opposed to sulfur compound reduction can be confidently distinguished. HMM search results are coupled to corresponding genes, which allows analysis of co-occurrence, synteny and genomic neighborhood. The HMMs were validated on an annotated test dataset and by cross-validation. We also proved its performance by exploring meta-assembled genomes isolated from samples from environments with active sulfur cycling, including members of the cable bacteria, novel Acidobacteria and assemblies from a sulfur-rich glacier, and were able to replicate and extend previous reports.