Two functionally distinct NADP+-dependent ferredoxin oxidoreductases maintain the primary redox balance of Pyrococcus furiosus

Structures and Electron Transport Paths in the Four Families of Flavin-Based Electron Bifurcation Enzymes

Oxidoreductases facilitating electron transfer between molecules are pivotal in metabolic pathways. Flavin-based electron bifurcation (FBEB), a recently discovered energy coupling mechanism in oxidoreductases, enables the reversible division of electron pairs into two acceptors, bridging exergonic and otherwise unfeasible endergonic reactions. This chapter explores the four distinct FBEB complex families and highlights a decade of structural insights into FBEB complexes. In this chapter, we discuss the architecture, electron transfer routes, and conformational changes across all FBEB families, revealing the structural foundation that facilitate these remarkable functions.

Identification of a deep-branching thermophilic clade sheds light on early bacterial evolution

It has been proposed that early bacteria, or even the last universal common ancestor of all cells, were thermophilic. However, research on the origin and evolution of thermophily is hampered by the difficulties associated with the isolation of deep-branching thermophilic microorganisms in pure culture. Here, we isolate a deep-branching thermophilic bacterium from a deep-sea hydrothermal vent, using a two-step cultivation strategy ("Subtraction-Suboptimal", StS) designed to isolate rare organisms. The bacterium, which we name Zhurongbacter thermophilus 3DAC, is a sulfur-reducing heterotroph that is phylogenetically related to Coprothermobacterota and other thermophilic bacterial groups, forming a clade that seems to represent a major, early-diverging bacterial lineage. The ancestor of this clade might be a thermophilic, strictly anaerobic, motile, hydrogen-dependent, and mixotrophic bacterium. Thus, our study provides insights into the early evolution of thermophilic bacteria.

Optimizing Strategies for Bio-Based Ethanol Production Using Genome-Scale Metabolic Modeling of the Hyperthermophilic Archaeon, Pyrococcus furiosus

A genome-scale metabolic model, encompassing a total of 623 genes, 727 reactions, and 865 metabolites, was developed for Pyrococcus furiosus, an archaeon that grows optimally at 100°C by carbohydrate and peptide fermentation. The model uses subsystem-based genome annotation, along with extensive manual curation of 237 gene-reaction associations including those involved in central carbon metabolism, amino acid metabolism, and energy metabolism. The redox and energy balance of P. furiosus was investigated through random sampling of flux distributions in the model during growth on disaccharides. The core energy balance of the model was shown to depend on high acetate production and the coupling of a sodium-dependent ATP synthase and membrane-bound hydrogenase, which generates a sodium gradient in a ferredoxin-dependent manner, aligning with existing understanding of P. furiosus metabolism. The model was utilized to inform genetic engineering designs that favor the production of ethanol over acetate by implementing an NADPH and CO-dependent energy economy. The P. furiosus model is a powerful tool for understanding the relationship between generation of end products and redox/energy balance at a systems-level that will aid in the design of optimal engineering strategies for production of bio-based chemicals and fuels. IMPORTANCE The bio-based production of organic chemicals provides a sustainable alternative to fossil-based production in the face of today's climate challenges. In this work, we present a genome-scale metabolic reconstruction of Pyrococcus furiosus, a well-established platform organism that has been engineered to produce a variety of chemicals and fuels. The metabolic model was used to design optimal engineering strategies to produce ethanol. The redox and energy balance of P. furiosus was examined in detail, which provided useful insights that will guide future engineering designs.

Evidence for a Putative Isoprene Reductase in Acetobacterium wieringae

Recent discoveries of isoprene-metabolizing microorganisms suggest they might play an important role in the global isoprene budget. Under anoxic conditions, isoprene can be used as an electron acceptor and is reduced to methylbutene. This study describes the proteogenomic profiling of an isoprene-reducing bacterial culture to identify organisms and genes responsible for the isoprene hydrogenation reaction. A metagenome-assembled genome (MAG) of the most abundant (89% relative abundance) lineage in the enrichment, Acetobacterium wieringae, was obtained. Comparative proteogenomics and reverse transcription-PCR (RT-PCR) identified a putative five-gene operon from the A. wieringae MAG upregulated during isoprene reduction. The operon encodes a putative oxidoreductase, three pleiotropic nickel chaperones (2 × HypA, HypB), and one 4Fe-4S ferredoxin. The oxidoreductase is proposed as the putative isoprene reductase with a binding site for NADH, flavin adenine dinucleotide (FAD), two pairs of canonical [4Fe-4S] clusters, and a putative iron-sulfur cluster site in a Cys₆-bonding environment. Well-studied Acetobacterium strains, such as A. woodii DSM 1030, A. wieringae DSM 1911, or A. malicum DSM 4132, do not encode the isoprene-regulated operon but encode, like many other bacteria, a homolog of the putative isoprene reductase (~47 to 49% amino acid sequence identity). Uncharacterized homologs of the putative isoprene reductase are observed across the Firmicutes, Spirochaetes, Tenericutes, Actinobacteria, Chloroflexi, Bacteroidetes, and Proteobacteria, suggesting the ability of biohydrogenation of unfunctionalized conjugated doubled bonds in other unsaturated hydrocarbons. IMPORTANCE Isoprene was recently shown to act as an electron acceptor for a homoacetogenic bacterium. The focus of this study is the molecular basis for isoprene reduction. By comparing a genome from our isoprene-reducing enrichment culture, dominated by Acetobacterium wieringae, with genomes of other Acetobacterium lineages that do not reduce isoprene, we shortlisted candidate genes for isoprene reduction. Using comparative proteogenomics and reverse transcription-PCR we have identified a putative five-gene operon encoding an oxidoreductase referred to as putative isoprene reductase.

Evolving a New Electron Transfer Pathway for Nitrogen Fixation Uncovers an Electron Bifurcating-Like Enzyme Involved in Anaerobic Aromatic Compound Degradation

Nitrogenase is the key enzyme involved in nitrogen fixation and uses low potential electrons delivered by ferredoxin (Fd) or flavodoxin (Fld) to reduce dinitrogen gas (N₂) to produce ammonia, generating hydrogen gas (H₂) as an obligate product of this activity. Although the phototrophic alphaproteobacterium Rhodopseudomonas palustris encodes multiple proteins that can reduce Fd, the FixABCX complex is the only one shown to support nitrogen fixation, and R. palustris Fix^- mutants grow poorly under nitrogen-fixing conditions. To investigate how native electron transfer chains (ETCs) can be redirected toward nitrogen fixation, we leveraged the strong selective pressure of nitrogen limitation to isolate a suppressor of an R. palustris ΔfixC strain that grows under nitrogen-fixing conditions. We found two mutations were required to restore growth under nitrogen-fixing conditions in the absence of functional FixABCX. One mutation was in the gene encoding the primary Fd involved in nitrogen fixation, fer1, and the other mutation was in aadN, which encodes a homolog of NAD⁺-dependent Fd:NADPH oxidoreductase (Nfn). We present evidence that AadN plays a role in electron transfer to benzoyl coenzyme A reductase, the key enzyme involved in anaerobic aromatic compound degradation. Our data support a model where the ETC for anaerobic aromatic compound degradation was repurposed to support nitrogen fixation in the ΔfixC suppressor strain. IMPORTANCE There is increasing evidence that protein electron carriers like Fd evolved to form specific partnerships with select electron donors and acceptors to keep native electron transfer pathways insulated from one another. This makes it challenging to integrate a Fd-dependent pathway such as biological nitrogen fixation into non-nitrogen-fixing organisms and provide the high-energy reducing power needed to fix nitrogen. Here, we show that amino acid substitutions in an electron donor for anaerobic aromatic compound degradation and an Fd involved in nitrogen fixation enabled electron transfer to nitrogenase. This study provides a model system to understand electron transfer chain specificity and how new electron transfer pathways can be evolved for biotechnologically valuable pathways like nitrogen fixation.

Proteome-wide 3D structure prediction provides insights into the ancestral metabolism of ancient archaea and bacteria

Ancestral metabolism has remained controversial due to a lack of evidence beyond sequence-based reconstructions. Although prebiotic chemists have provided hints that metabolism might originate from non-enzymatic protometabolic pathways, gaps between ancestral reconstruction and prebiotic processes mean there is much that is still unknown. Here, we apply proteome-wide 3D structure predictions and comparisons to investigate ancestorial metabolism of ancient bacteria and archaea, to provide information beyond sequence as a bridge to the prebiotic processes. We compare representative bacterial and archaeal strains, which reveal surprisingly similar physiological and metabolic characteristics via microbiological and biophysical experiments. Pairwise comparison of protein structures identify the conserved metabolic modules in bacteria and archaea, despite interference from overly variable sequences. The conserved modules (for example, middle of glycolysis, partial TCA, proton/sulfur respiration, building block biosynthesis) constitute the basic functions that possibly existed in the archaeal-bacterial common ancestor, which are remarkably consistent with the experimentally confirmed protometabolic pathways. These structure-based findings provide a new perspective to reconstructing the ancestral metabolism and understanding its origin, which suggests high-throughput protein 3D structure prediction is a promising approach, deserving broader application in future ancestral exploration.

An unprecedented function for a tungsten-containing oxidoreductase

Five tungstopterin-containing oxidoreductases were characterized from the hyperthermophile Pyrococcus furiosus. Each enzyme catalyzes the reversible conversion of one or more aldehydes to the corresponding carboxylic acid, but they have different specificities. The physiological functions of only two of these enzymes are known: one, termed GAPOR, is a glycolytic enzyme that oxidizes glyceraldehyde-3-phosphate, while the other, termed AOR, oxidizes multiple aldehydes generated during peptide fermentation. Two of the enzymes have known structures (AOR and FOR). Herein, we focus on WOR5, the fifth tungstopterin enzyme to be discovered in P. furiosus. Expression of WOR5 was previously shown to be increased during cold shock (growth at 72 ℃), although the physiological substrate is not known. To gain insight into WOR5 function, we sought to determine both its structure and identify its intracellular substrate. Crystallization experiments were performed with a concentrated cytoplasmic extract of P. furiosus grown at 72 ℃ and the structure of WOR5 was deduced from the crystals that were obtained. In contrast to a previous report, WOR5 is heterodimeric containing an additional polyferredoxin-like subunit with four [4Fe-4S] clusters. The active site structure of WOR5 is substantially different from that of AOR and FOR and the significant electron density observed adjacent to the tungsten cofactor of WOR5 was modeled as an aliphatic sulfonate. Biochemical assays and product analysis confirmed that WOR5 is an aliphatic sulfonate ferredoxin oxidoreductase (ASOR). A catalytic mechanism for ASOR is proposed based on the structural information and the potential role of ASOR in the cold-shock response is discussed.

Culexarchaeia, a novel archaeal class of anaerobic generalists inhabiting geothermal environments

Geothermal environments, including terrestrial hot springs and deep-sea hydrothermal sediments, often contain many poorly understood lineages of archaea. Here, we recovered ten metagenome-assembled genomes (MAGs) from geothermal sediments and propose that they constitute a new archaeal class within the TACK superphylum, "Candidatus Culexarchaeia", named after the Culex Basin in Yellowstone National Park. Culexarchaeia harbor distinct sets of proteins involved in key cellular processes that are either phylogenetically divergent or are absent from other closely related TACK lineages, with a particular divergence in cell division and cytoskeletal proteins. Metabolic reconstruction revealed that Culexarchaeia have the capacity to metabolize a wide variety of organic and inorganic substrates. Notably, Culexarchaeia encode a unique modular, membrane associated, and energy conserving [NiFe]-hydrogenase complex that potentially interacts with heterodisulfide reductase (Hdr) subunits. Comparison of this [NiFe]-hydrogenase complex with similar complexes from other archaea suggests that interactions between membrane associated [NiFe]-hydrogenases and Hdr may be more widespread than previously appreciated in both methanogenic and non-methanogenic lifestyles. The analysis of Culexarchaeia further expands our understanding of the phylogenetic and functional diversity of lineages within the TACK superphylum and the ecology, physiology, and evolution of these organisms in extreme environments.

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.

A Novel NADP-Dependent Formate Dehydrogenase From the Hyperthermophilic Archaeon Thermococcus onnurineus NA1

The genome of the hyperthermophilic archaeon Thermococcus onnurineus NA1 contains three copies of the formate dehydrogenase (FDH) gene, fdh1, fdh2, and fdh3. Previously, we reported that fdh2, clustered with genes encoding the multimeric membrane-bound hydrogenase and cation/proton antiporter, was essential for formate-dependent growth with H₂ production. However, the functionality of the other two FDH-coding genes has not yet been elucidated. Herein, we purified and characterized cytoplasmic Fdh3 to understand its functionality. The purified Fdh3 was identified to be composed of a tungsten-containing catalytic subunit (Fdh3A), an NAD(P)-binding protein (Fdh3B), and two Fe-S proteins (Fdh3G1 and Fdh3G2). Fdh3 oxidized formate with specific activities of 241.7 U/mg and 77.4 U/mg using methyl viologen and NADP⁺ as electron acceptors, respectively. While most FDHs exhibited NAD⁺-dependent formate oxidation activity, the Fdh3 of T. onnurineus NA1 showed a strong preference for NADP⁺ over NAD⁺ as a cofactor. The catalytic efficiency (k_cat/K_m) of Fdh3 for NADP⁺ was measured to be 5,281 mM⁻¹ s⁻¹, which is the highest among NADP-dependent FDHs known to date. Structural modeling suggested that Arg²⁰⁴ and Arg²⁰⁵ of Fdh3B may contribute to the stabilization of the 2′-phosphate of NADP(H). Fdh3 could also use ferredoxin as an electron acceptor to oxidize formate with a specific activity of 0.83 U/mg. Furthermore, Fdh3 showed CO₂ reduction activity using reduced ferredoxin or NADPH as an electron donor with a specific activity of 0.73 U/mg and 1.0 U/mg, respectively. These results suggest a functional role of Fdh3 in disposing of reducing equivalents by mediating electron transfer between formate and NAD(P)H or ferredoxin.

The Piezo-Hyperthermophilic Archaeon Thermococcus piezophilus Regulates Its Energy Efficiency System to Cope With Large Hydrostatic Pressure Variations

Deep-sea ecosystems share a common physical parameter, namely high hydrostatic pressure (HHP). Some of the microorganisms isolated at great depths have a high physiological plasticity to face pressure variations. The adaptive strategies by which deep-sea microorganisms cope with HHP variations remain to be elucidated, especially considering the extent of their biotopes on Earth. Herein, we investigated the gene expression patterns of Thermococcus piezophilus, a piezohyperthermophilic archaeon isolated from the deepest hydrothermal vent known to date, under sub-optimal, optimal and supra-optimal pressures (0.1, 50, and 90 MPa, respectively). At stressful pressures [sub-optimal (0.1 MPa) and supra-optimal (90 MPa) conditions], no classical stress response was observed. Instead, we observed an unexpected transcriptional modulation of more than a hundred gene clusters, under the putative control of the master transcriptional regulator SurR, some of which are described as being involved in energy metabolism. This suggests a fine-tuning effect of HHP on the SurR regulon. Pressure could act on gene regulation, in addition to modulating their expression.

Genome-Scale Metabolic Model of Caldicellulosiruptor bescii Reveals Optimal Metabolic Engineering Strategies for Bio-based Chemical Production

Metabolic modeling was used to examine potential bottlenecks that could be encountered for metabolic engineering of the cellulolytic extreme thermophile Caldicellulosiruptor bescii to produce bio-based chemicals from plant biomass. The model utilizes subsystems-based genome annotation, targeted reconstruction of carbohydrate utilization pathways, and biochemical and physiological experimental validations. Specifically, carbohydrate transport and utilization pathways involving 160 genes and their corresponding functions were incorporated, representing the utilization of C5/C6 monosaccharides, disaccharides, and polysaccharides such as cellulose and xylan. To illustrate its utility, the model predicted that optimal production from biomass-based sugars of the model product, ethanol, was driven by ATP production, redox balancing, and proton translocation, mediated through the interplay of an ATP synthase, a membrane-bound hydrogenase, a bifurcating hydrogenase, and a bifurcating NAD- and NADP-dependent oxidoreductase. These mechanistic insights guided the design and optimization of new engineering strategies for product optimization, which were subsequently tested in the C. bescii model, showing a nearly 2-fold increase in ethanol yields. The C. bescii model provides a useful platform for investigating the potential redox controls that mediate the carbon and energy flows in metabolism and sets the stage for future design of engineering strategies aiming at optimizing the production of ethanol and other bio-based chemicals. IMPORTANCE The extremely thermophilic cellulolytic bacterium, Caldicellulosiruptor bescii, degrades plant biomass at high temperatures without any pretreatments and can serve as a strategic platform for industrial applications. The metabolic engineering of C. bescii, however, faces potential bottlenecks in bio-based chemical productions. By simulating the optimal ethanol production, a complex interplay between redox balancing and the carbon and energy flow was revealed using a C. bescii genome-scale metabolic model. New engineering strategies were designed based on an improved mechanistic understanding of the C. bescii metabolism, and the new designs were modeled under different genetic backgrounds to identify optimal strategies. The C. bescii model provided useful insights into the metabolic controls of this organism thereby opening up prospects for optimizing production of a wide range of bio-based chemicals.

Iron-sulfur flavoenzymes: the added value of making the most ancient redox cofactors and the versatile flavins work together

Iron-sulfur (Fe-S) flavoproteins form a broad and growing class of complex, multi-domain and often multi-subunit proteins coupling the most ancient cofactors (the Fe-S clusters) and the most versatile coenzymes (the flavin coenzymes, FMN and FAD). These enzymes catalyse oxidoreduction reactions usually acting as switches between donors of electron pairs and acceptors of single electrons, and vice versa. Through selected examples, the enzymes' structure−function relationships with respect to rate and directionality of the electron transfer steps, the role of the apoprotein and its dynamics in modulating the electron transfer process will be discussed.

The Pyrococcus furiosus ironome is dominated by [Fe4S4]2+ clusters or thioferrate-like iron depending on the availability of elemental sulfur

Pyrococcus furiosus is a hyperthermophilic anaerobic archaeon whose metabolism depends on whether elemental sulfur is (+S⁰) or is not (-S⁰) included in growth medium. Under +S⁰ conditions, expression of respiratory hydrogenase declines while respiratory membrane-bound sulfane reductase and the putative iron-storage protein IssA increase. Our objective was to investigate the iron content of WT and ΔIssA cells under these growth conditions using Mössbauer spectroscopy. WT-S⁰ cells contained ∼1 mM Fe, with ∼85% present as two spectroscopically distinct forms of S = 0 [Fe₄S₄]²⁺ clusters; the remainder was mainly high-spin Fe^II. WT+S⁰ cells contained 5 to 9 mM Fe, with 75 to 90% present as magnetically ordered thioferrate-like (TFL) iron nanoparticles. TFL iron was similar to chemically defined thioferrates; both consisted of Fe^III ions coordinated by an S₄ environment, and both exhibited strong coupling between particles causing high applied fields to have little spectral effect. At high temperatures with magnetic hyperfine interactions abolished, TFL iron exhibited two doublets overlapping those of [Fe₄S₄]²⁺ clusters in -S⁰ cells. This coincidence arose because of similar coordination environments of TFL iron and cluster iron. The TFL structure was more heterogeneous in the presence of IssA. Presented data suggest that IssA may coordinate insoluble iron sulfides as TFL iron, formed as a byproduct of anaerobic sulfur respiration under high iron conditions, which thereby reduces its toxicity to the cell. This was the first Mössbauer characterization of the ironome of an archaeon, and it illustrates differences relative to the iron content of better-studied bacteria such as Escherichia coli.

Flavins in the electron bifurcation process

The discovery of a new energy-coupling mechanism termed flavin-based electron bifurcation (FBEB) in 2008 revealed a novel field of application for flavins in biology. The key component is the bifurcating flavin endowed with strongly inverted one-electron reduction potentials (FAD/FAD•^- ≪ FAD•^-/FADH^-) that cooperatively transfers in its reduced state one low and one high-energy electron into different directions and thereby drives an endergonic with an exergonic reduction reaction. As energy splitting at the bifurcating flavin apparently implicates one-electron chemistry, the FBEB machinery has to incorporate prior to and behind the central bifurcating flavin 2e-to-1e and 1e-to-2e switches, frequently also flavins, for oxidizing variable medium-potential two-electron donating substrates and for reducing high-potential two-electron accepting substrates. The one-electron carriers ferredoxin or flavodoxin serve as low-potential (high-energy) electron acceptors, which power endergonic processes almost exclusively in obligate anaerobic microorganisms to increase the efficiency of their energy metabolism. In this review, we outline the global organization of FBEB enzymes, the functions of the flavins therein and the surrounding of the isoalloxazine rings by which their reduction potentials are specifically adjusted in a finely tuned energy landscape.

The role of thermodynamic features on the functional activity of electron bifurcating enzymes

Electron bifurcation is a biological mechanism to drive a thermodynamically unfavorable redox reaction through direct coupling with an exergonic reaction. This process allows microorganisms to generate high energy reducing equivalents in order to sustain life and is often found in anaerobic metabolism, where the energy economy of the cell is poor. Recent work has revealed details of the redox energy landscapes for a variety of electron bifurcating enzymes, greatly expanding the understanding of how energy is transformed by this unique mechanism. Here we highlight the plasticity of these emerging landscapes, what is known regarding their mechanistic underpinnings, and provide a context for interpreting their biochemical activity within the physiological framework. We conclude with an outlook for propelling the field toward an integrative understanding of the impact of electron bifurcation.

Overcoming Energetic Barriers in Acetogenic C1 Conversion

Currently one of the biggest challenges for society is to combat global warming. A solution to this global threat is the implementation of a CO₂-based bioeconomy and a H₂-based bioenergy economy. Anaerobic lithotrophic bacteria such as the acetogenic bacteria are key players in the global carbon and H₂ cycle and thus prime candidates as driving forces in a H₂- and CO₂-bioeconomy. Naturally, they convert two molecules of CO₂via the Wood-Ljungdahl pathway (WLP) to one molecule of acetyl-CoA which can be converted to different C2-products (acetate or ethanol) or elongated to C4 (butyrate) or C5-products (caproate). Since there is no net ATP generation from acetate formation, an electron-transport phosphorylation (ETP) module is hooked up to the WLP. ETP provides the cell with additional ATP, but the ATP gain is very low, only a fraction of an ATP per mol of acetate. Since acetogens live at the thermodynamic edge of life, metabolic engineering to obtain high-value products is currently limited by the low energy status of the cells that allows for the production of only a few compounds with rather low specificity. To set the stage for acetogens as production platforms for a wide range of bioproducts from CO₂, the energetic barriers have to be overcome. This review summarizes the pathway, the energetics of the pathway and describes ways to overcome energetic barriers in acetogenic C1 conversion.

Cross-Stress Adaptation in a Piezophilic and Hyperthermophilic Archaeon From Deep Sea Hydrothermal Vent

Hyperthermophiles, living in environments above 80°C and usually coupling with multi-extreme environmental stresses, have drawn great attention due to their application potential in biotechnology and being the primitive extant forms of life. Studies on their survival and adaptation mechanisms have extended our understanding on how lives thrive under extreme conditions. During these studies, the "cross-stress" behavior in various organisms has been observed between the extreme high temperature and other environmental stresses. Despite the broad observation, the global view of the cross-stress behavior remains unclear in hyperthermophiles, leaving a knowledge gap in our understanding of extreme adaptation. In this study, we performed a global quantitative proteomic analysis under extreme temperatures, pH, hydrostatic pressure (HP), and salinity on an archaeal strain, Thermococcus eurythermalis A501, which has outstanding growth capability on a wide range of temperatures (50-100°C), pH (4-9), and HPs (0.1-70 MPa), but a narrow range of NaCl (1.0-5.0 %, w/v). The proteomic analysis (79.8% genome coverage) demonstrated that approximately 61.5% of the significant differentially expressed proteins (DEPs) responded to multiple stresses. The responses to most of the tested stresses were closely correlated, except the responses to high salinity and low temperature. The top three enriched universal responding processes include the biosynthesis and protection of macromolecules, biosynthesis and metabolism of amino acids, ion transport, and binding activities. In addition, this study also revealed that the specific dual-stress responding processes, such as the membrane lipids for both cold and HP stresses and the signal transduction for both hyperosmotic and heat stresses, as well as the sodium-dependent energetic processes might be the limiting factor of the growth range in salinity. The present study is the first to examine the global cross-stress responses in a piezophilic hyperthermophile at the proteomic level. Our findings provide direct evidences of the cross-stress adaptation strategy (33.5% of coding-genes) to multiple stresses and highlight the specific and unique responding processes (0.22-0.63% of coding genes for each) to extreme temperature, pH, salinity, and pressure, which are highly relevant to the fields of evolutionary biology as well as next generation industrial biotechnology (NGIB).

The Sporomusa type Nfn is a novel type of electron-bifurcating transhydrogenase that links the redox pools in acetogenic bacteria

Flavin-based electron bifurcation is a long hidden mechanism of energetic coupling present mainly in anaerobic bacteria and archaea that suffer from energy limitations in their environment. Electron bifurcation saves precious cellular ATP and enables lithotrophic life of acetate-forming (acetogenic) bacteria that grow on H₂ + CO₂ by the only pathway that combines CO₂ fixation with ATP synthesis, the Wood–Ljungdahl pathway. The energy barrier for the endergonic reduction of NADP⁺, an electron carrier in the Wood–Ljungdahl pathway, with NADH as reductant is overcome by an electron-bifurcating, ferredoxin-dependent transhydrogenase (Nfn) but many acetogens lack nfn genes. We have purified a ferredoxin-dependent NADH:NADP⁺ oxidoreductase from Sporomusa ovata, characterized the enzyme biochemically and identified the encoding genes. These studies led to the identification of a novel, Sporomusa type Nfn (Stn), built from existing modules of enzymes such as the soluble [Fe–Fe] hydrogenase, that is widespread in acetogens and other anaerobic bacteria.

Engineering the cellulolytic extreme thermophile Caldicellulosiruptor bescii to reduce carboxylic acids to alcohols using plant biomass as the energy source

Caldicellulosiruptor bescii is the most thermophilic cellulolytic organism yet identified (T_opt 78 °C). It grows on untreated plant biomass and has an established genetic system thereby making it a promising microbial platform for lignocellulose conversion to bio-products. Here, we investigated the ability of engineered C. bescii to generate alcohols from carboxylic acids. Expression of aldehyde ferredoxin oxidoreductase (aor from Pyrococcus furiosus) and alcohol dehydrogenase (adhA from Thermoanaerobacter sp. X514) enabled C. bescii to generate ethanol from crystalline cellulose and from biomass by reducing the acetate produced by fermentation. Deletion of lactate dehydrogenase in a strain expressing the AOR-Adh pathway increased ethanol production. Engineered strains also converted exogenously supplied organic acids (isobutyrate and n-caproate) to the corresponding alcohol (isobutanol and hexanol) using both crystalline cellulose and switchgrass as sources of reductant for alcohol production. This is the first instance of an acid to alcohol conversion pathway in a cellulolytic microbe.

Modification of the glycolytic pathway in Pyrococcus furiosus and the implications for metabolic engineering

The key difference in the modified Embden-Meyerhof glycolytic pathway in hyperthermophilic Archaea, such as Pyrococcus furiosus, occurs at the conversion from glyceraldehyde-3-phosphate (GAP) to 3-phosphoglycerate (3-PG) where the typical intermediate 1,3-bisphosphoglycerate (1,3-BPG) is not present. The absence of the ATP-yielding step catalyzed by phosphoglycerate kinase (PGK) alters energy yield, redox energetics, and kinetics of carbohydrate metabolism. Either of the two enzymes, ferredoxin-dependent glyceraldehyde-3-phosphate ferredoxin oxidoreductase (GAPOR) or NADP⁺-dependent non-phosphorylating glyceraldehyde-3-phosphate dehydrogenase (GAPN), responsible for this "bypass" reaction, could be deleted individually without impacting viability, albeit with differences in native fermentation product profiles. Furthermore, P. furiosus was viable in the gluconeogenic direction (growth on pyruvate or peptides plus elemental sulfur) in a ΔgapnΔgapor strain. Ethanol was utilized as a proxy for potential heterologous products (e.g., isopropanol, butanol, fatty acids) that require reducing equivalents (e.g., NAD(P)H, reduced ferredoxin) generated from glycolysis. Insertion of a single gene encoding the thermostable NADPH-dependent primary alcohol dehydrogenase (adhA) (Tte_0696) from Caldanaerobacter subterraneus, resulted in a strain producing ethanol via the previously established aldehyde oxidoreductase (AOR) pathway. This strain demonstrated a high ratio of ethanol over acetate (> 8:1) at 80 °C and enabled ethanol production up to 85 °C, the highest temperature for bio-ethanol production reported to date.

Small protein folds at the root of an ancient metabolic network

Life on Earth is driven by electron transfer reactions catalyzed by a suite of enzymes that comprise the superfamily of oxidoreductases (Enzyme Classification EC1). Most modern oxidoreductases are complex in their structure and chemistry and must have evolved from a small set of ancient folds. Ancient oxidoreductases from the Archean Eon between ca. 3.5 and 2.5 billion years ago have been long extinct, making it challenging to retrace evolution by sequence-based phylogeny or ancestral sequence reconstruction. However, three-dimensional topologies of proteins change more slowly than sequences. Using comparative structure and sequence profile-profile alignments, we quantify the similarity between proximal cofactor-binding folds and show that they are derived from a common ancestor. We discovered that two recurring folds were central to the origin of metabolism: ferredoxin and Rossmann-like folds. In turn, these two folds likely shared a common ancestor that, through duplication, recruitment, and diversification, evolved to facilitate electron transfer and catalysis at a very early stage in the origin of metabolism.

DsrL mediates electron transfer between NADH and rDsrAB in Allochromatium vinosum

Dissimilatory sulphite reductase DsrAB occurs in sulphate/sulphite-reducing prokaryotes, in sulphur disproportionators and also in sulphur oxidizers, where it functions in reverse. Predictions of physiological traits in metagenomic studies relying on the presence of dsrAB, other dsr genes or combinations thereof suffer from the lack of information on crucial Dsr proteins. The iron-sulphur flavoprotein DsrL is an example of this group. It has a documented essential function during sulphur oxidation and was recently also found in some metagenomes of probable sulphate and sulphite reducers. Here, we show that DsrL and reverse acting rDsrAB can form a complex and are copurified from the phototrophic sulphur oxidizer Allochromatium vinosum. Recombinant DsrL exhibits NAD(P)H:acceptor oxidoreductase activity with a strong preference for NADH over NADPH. In vitro, the rDsrABL complex effectively catalyses NADH-dependent sulphite reduction, which is strongly enhanced by the sulphur-binding protein DsrC. Our work reveals NAD⁺ as suitable in vivo electron acceptor for sulphur oxidation in organisms operating the rDsr pathway and points to reduced nicotinamide adenine dinucleotides as electron donors for sulphite reduction in sulphate/sulphite-reducing prokaryotes that contain DsrL. In addition, dsrL cannot be used as a marker distinguishing sulphate/sulphite reducers and sulphur oxidizers in metagenomic studies without further analysis.

The catalytic mechanism of electron-bifurcating electron transfer flavoproteins (ETFs) involves an intermediary complex with NAD<sup/>

Electron bifurcation plays a key role in anaerobic energy metabolism, but it is a relatively new discovery, and only limited mechanistic information is available on the diverse enzymes that employ it. Herein, we focused on the bifurcating electron transfer flavoprotein (ETF) from the hyperthermophilic archaeon Pyrobaculum aerophilum The EtfABCX enzyme complex couples NADH oxidation to the endergonic reduction of ferredoxin and exergonic reduction of menaquinone. We developed a model for the enzyme structure by using nondenaturing MS, cross-linking, and homology modeling in which EtfA, -B, and -C each contained FAD, whereas EtfX contained two [4Fe-4S] clusters. On the basis of analyses using transient absorption, EPR, and optical titrations with NADH or inorganic reductants with and without NAD+, we propose a catalytic cycle involving formation of an intermediary NAD+-bound complex. A charge transfer signal revealed an intriguing interplay of flavin semiquinones and a protein conformational change that gated electron transfer between the low- and high-potential pathways. We found that despite a common bifurcating flavin site, the proposed EtfABCX catalytic cycle is distinct from that of the genetically unrelated bifurcating NADH-dependent ferredoxin NADP+ oxidoreductase (NfnI). The two enzymes particularly differed in the role of NAD+, the resting and bifurcating-ready states of the enzymes, how electron flow is gated, and the two two-electron cycles constituting the overall four-electron reaction. We conclude that P. aerophilum EtfABCX provides a model catalytic mechanism that builds on and extends previous studies of related bifurcating ETFs and can be applied to the large bifurcating ETF family.

Distribution, Evolution, Catalytic Mechanism, and Physiological Functions of the Flavin-Based Electron-Bifurcating NADH-Dependent Reduced Ferredoxin: NADP+ Oxidoreductase

NADH-dependent reduced ferredoxin:NADP⁺ oxidoreductase (Nfn) is an electron-bifurcating enzyme first discovered in the strict anaerobes Clostridium kluyveri and Moorella thermoacetica. In vivo, Nfn catalyzes the endergonic reduction of NADP⁺ with NADH coupled to the exergonic reduction of NADP⁺ with reduced ferredoxin. Most Nfn homologs consist of two subunits, although in certain species Nfn homologs are fused. In contrast to other electron-bifurcating enzymes, Nfn possess a simpler structure. Therefore, Nfn becomes a perfect model to determine the mechanism of flavin-based electron bifurcation, which is a novel energy coupling mode distributed among anaerobic bacteria and archaea. The crystal structures of Nfn from Thermotoga maritima and Pyrococcus furiosus are known, and studies have shown that the FAD molecule of the NfnB (b-FAD) is the site of electron bifurcation, and other cofactors, including a [2Fe2S] cluster, two [4Fe4S] clusters, and the FAD molecule on the NfnA subunit, contribute to electron transfer. Further, the short-lived anionic flavin semiquinone (ASQ) state of b-FAD is essential for electron bifurcation. Nfn homologs are widely distributed among microbes, including bacteria, archaea, and probably eukaryotes, most of which are anaerobes despite that certain species are facultative microbes and even aerobes. Moreover, potential evidence shows that lateral gene transfer may occur in the evolution of this enzyme. Nfn homologs present four different structural patterns, including the well-characterized NfnAB and three different kinds of fused Nfn homologs whose detailed properties have not been characterized. These findings indicate that gene fusion/fission and gene rearrangement may contribute to the evolution of this enzyme. Under physiological conditions, Nfn catalyzes the reduction of NADP⁺ with NADH and reduced ferredoxin, which is then used in certain NADPH-dependent reactions. Deletion of nfn in several microbes causes low growth and redox unbalance and may influence the distribution of fermentation products. It’s also noteworthy that different Nfn homologs perform different functions according to its circumstance. Physiological functions of Nfn indicate that it can be a potential tool in the metabolic engineering of industrial microorganisms, which can regulate the redox potential in vivo.

Sulfur-dependent microbial lifestyles: deceptively flexible roles for biochemically versatile enzymes

A wide group of microbes are able to "make a living" on Earth by basing their energetic metabolism on inorganic sulfur compounds. Because of their range of stable redox states, sulfur and inorganic sulfur compounds can be utilized as either oxidants or reductants in a diverse array of energy-conserving reactions. In this review the major enzymes and basic chemistry of sulfur-based respiration and chemolithotrophy are outlined. The reversibility and versatility of these enzymes, however, means that they can often be used in multiple ways, and several cases are discussed in which enzymes which are considered to be hallmarks of a particular respiratory or lithotrophic process have been found to be used in other, often opposing, metabolic processes. These results emphasize the importance of taking into account the geochemistry, biochemistry and microbiology of an organism and/or environment when trying to interpret the function of a particular sulfur-dependent redox enzyme.

Probing Cascade complex composition and stability using native mass spectrometry techniques

Adaptive prokaryotic immune systems rely on clustered regularly interspaced short palindromic repeats and their associated genes to provide the components necessary to clear infection by foreign genetic elements. These immune systems are based on highly specific nucleases that bind DNA or RNA and, upon sequence recognition, degrade the bound nucleic acid. Because of their specificity, CRISPR-Cas systems are being co-opted to edit genes in eukaryotic cells. While the general function of these systems is well understood, an understanding of mechanistic details to facilitate engineering and application to this new arena remains a topic of intense study. Here, we present two methods that have been successfully used to study the structure and mechanism of the Type IE CRISPR system, Cascade, from Escherichia coli. We provide the protocol for a typical native mass spectrometry experiment which, because it allows for analysis of a protein complex without disruption of the noncovalent interactions within the complex, can be used to determine complex composition, architecture, and relative affinity between subunits. We, also, provide the protocol for intact protein hydrogen-deuterium exchange mass spectrometry, which provides insight into the overall conformational stability of the complex and changes in complex stability based on conditions such as substrate binding. Investigating the solution-phase structure, stability, and dynamics of these complexes improves the overall understanding of the mechanism facilitating engineered adjustments to function or utility.

Characterization of membrane-bound sulfane reductase: A missing link in the evolution of modern day respiratory complexes

Hyperthermophilic archaea contain a hydrogen gas-evolving,respiratory membrane-bound NiFe-hydrogenase (MBH) that is very closely related to the aerobic respiratory complex I. During growth on elemental sulfur (S°), these microorganisms also produce a homologous membrane-bound complex (MBX), which generates H₂S. MBX evolutionarily links MBH to complex I, but its catalytic function is unknown. Herein, we show that MBX reduces the sulfane sulfur of polysulfides by using ferredoxin (Fd) as the electron donor, and we rename it membrane-bound sulfane reductase (MBS). Two forms of affinity-tagged MBS were purified from genetically engineered Pyrococcus furiosus (a hyperthermophilic archaea species): the 13-subunit holoenzyme (S-MBS) and a cytoplasmic 4-subunit catalytic subcomplex (C-MBS). S-MBS and C-MBS reduced dimethyl trisulfide (DMTS) with comparable K_m (∼490 μm) and V _max values (12 μmol/min/mg). The MBS catalytic subunit (MbsL), but not that of complex I (NuoD), retains two of four NiFe-coordinating cysteine residues of MBH. However, these cysteine residues were not involved in MBS catalysis because a mutant P. furiosus strain (MbsLC85A/C385A) grew normally with S°. The products of the DMTS reduction and properties of polysulfides indicated that in the physiological reaction, MBS uses Fd (E _o' = -480 mV) to reduce sulfane sulfur (E _o' -260 mV) and cleave organic (RS _n R, n ≥ 3) and anionic polysulfides (S _n ^2-, n ≥ 4) but that it does not produce H₂S. Based on homology to MBH, MBS also creates an ion gradient for ATP synthesis. This work establishes the electrochemical reaction catalyzed by MBS that is intermediate in the evolution from proton- to quinone-reducing respiratory complexes.

Origin and Evolution of Flavin-Based Electron Bifurcating Enzymes

Twelve evolutionarily unrelated oxidoreductases form enzyme complexes that catalyze the simultaneous coupling of exergonic and endergonic oxidation–reduction reactions to circumvent thermodynamic barriers and minimize free energy loss in a process known as flavin-based electron bifurcation. Common to these 12 bifurcating (Bf) enzymes are protein-bound flavin, the proposed site of bifurcation, and the electron carrier ferredoxin. Despite the documented role of Bf enzymes in balancing the redox state of intracellular electron carriers and in improving the efficiency of cellular metabolism, a comprehensive description of the diversity and evolutionary history of Bf enzymes is lacking. Here, we report the taxonomic distribution, functional diversity, and evolutionary history of Bf enzyme homologs in 4,588 archaeal, bacterial, and eukaryal genomes and 3,136 community metagenomes. Bf homologs were primarily detected in the genomes of anaerobes, including those of sulfate-reducers, acetogens, fermenters, and methanogens. Phylogenetic analyses of Bf enzyme catalytic subunits (oxidoreductases) suggest they were not a property of the Last Universal Common Ancestor of Archaea and Bacteria, which is consistent with the limited and unique taxonomic distributions of enzyme homologs among genomes. Further, phylogenetic analyses of oxidoreductase subunits reveal that non-Bf homologs predate Bf homologs. These observations indicate that multiple independent recruitments of flavoproteins to existing oxidoreductases enabled coupling of numerous new electron Bf reactions. Consistent with the role of these enzymes in the energy metabolism of anaerobes, homologs of Bf enzymes were enriched in metagenomes from subsurface environments relative to those from surface environments. Phylogenetic analyses of homologs from metagenomes reveal that the earliest evolving homologs of most Bf enzymes are from subsurface environments, including fluids from subsurface rock fractures and hydrothermal systems. Collectively, these data suggest strong selective pressures drove the emergence of Bf enzyme complexes via recruitment of flavoproteins that allowed for an increase in the efficiency of cellular metabolism and improvement in energy capture in anaerobes inhabiting a variety of subsurface anoxic habitats where the energy yield of oxidation-reduction reactions is generally low.

The Role of Mass Spectrometry in Structural Studies of Flavin-Based Electron Bifurcating Enzymes

For decades, biologists and biochemists have taken advantage of atomic resolution structural models of proteins from X-ray crystallography, nuclear magnetic resonance spectroscopy, and more recently cryo-electron microscopy. However, not all proteins relent to structural analyses using these approaches, and as the depth of knowledge increases, additional data elucidating a mechanistic understanding of protein function is desired. Flavin-based electron bifurcating enzymes, which are responsible for producing high energy compounds through the simultaneous endergonic and exergonic reduction of two intercellular electron carriers (i.e., NAD⁺ and ferredoxin) are one class of proteins that have challenged structural biologists and in which there is great interest to understand the mechanism behind electron gating. A limited number of X-ray crystallography projects have been successful; however, it is clear that to understand how these enzymes function, techniques that can reveal detailed in solution information about protein structure, dynamics, and interactions involved in the bifurcating reaction are needed. In this review, we cover a general set of mass spectrometry-based techniques that, combined with protein modeling, are capable of providing information on both protein structure and dynamics. Techniques discussed include surface labeling, covalent cross-linking, native mass spectrometry, and hydrogen/deuterium exchange. We cover how biophysical data can be used to validate computationally generated protein models and develop mechanistic explanations for regulation and performance of enzymes and protein complexes. Our focus will be on flavin-based electron bifurcating enzymes, but the broad applicability of the techniques will be showcased.

Flavin-Based Electron Bifurcation, A New Mechanism of Biological Energy Coupling

There are two types of electron bifurcation (EB), either quinone- or flavin-based (QBEB/FBEB), that involve reduction of a quinone or flavin by a two-electron transfer and two reoxidations by a high- and low-potential one-electron acceptor with a reactive semiquinone intermediate. In QBEB, the reduced low-potential acceptor (cytochrome b) is exclusively used to generate ΔμH⁺. In FBEB, the "energy-rich" low-potential reduced ferredoxin or flavodoxin has dual function. It can give rise to ΔμH⁺/Na⁺ via a ferredoxin:NAD reductase (Rnf) or ferredoxin:proton reductase (Ech) or conducts difficult reductions such as CO₂ to CO. The QBEB membrane complexes are similar in structure and function and occur in all domains of life. In contrast, FBEB complexes are soluble and occur only in strictly anaerobic bacteria and archaea (FixABCX being an exception). The FBEB complexes constitute a group consisting of four unrelated families that contain (1) electron-transferring flavoproteins (EtfAB), (2) NAD(P)H dehydrogenase (NuoF homologues), (3) heterodisulfide reductase (HdrABC) or HdrABC homologues, and (4) NADH-dependent ferredoxin:NADP reductase (NfnAB). The crystal structures and electron transport of EtfAB-butyryl-CoA dehydrogenase and NfnAB are compared with those of complex III of the respiratory chain (cytochrome bc₁), whereby unexpected common features have become apparent.

Flavin-Based Electron Bifurcation, Ferredoxin, Flavodoxin, and Anaerobic Respiration With Protons (Ech) or NAD+ (Rnf) as Electron Acceptors: A Historical Review

Flavin-based electron bifurcation is a newly discovered mechanism, by which a hydride electron pair from NAD(P)H, coenzyme F₄₂₀H₂, H₂, or formate is split by flavoproteins into one-electron with a more negative reduction potential and one with a more positive reduction potential than that of the electron pair. Via this mechanism microorganisms generate low- potential electrons for the reduction of ferredoxins (Fd) and flavodoxins (Fld). The first example was described in 2008 when it was found that the butyryl-CoA dehydrogenase-electron-transferring flavoprotein complex (Bcd-EtfAB) of Clostridium kluyveri couples the endergonic reduction of ferredoxin (E₀′ = −420 mV) with NADH (−320 mV) to the exergonic reduction of crotonyl-CoA to butyryl-CoA (−10 mV) with NADH. The discovery was followed by the finding of an electron-bifurcating Fd- and NAD-dependent [FeFe]-hydrogenase (HydABC) in Thermotoga maritima (2009), Fd-dependent transhydrogenase (NfnAB) in various bacteria and archaea (2010), Fd- and H₂-dependent heterodisulfide reductase (MvhADG-HdrABC) in methanogenic archaea (2011), Fd- and NADH-dependent caffeyl-CoA reductase (CarCDE) in Acetobacterium woodii (2013), Fd- and NAD-dependent formate dehydrogenase (HylABC-FdhF2) in Clostridium acidi-urici (2013), Fd- and NADP-dependent [FeFe]-hydrogenase (HytA-E) in Clostridium autoethanogrenum (2013), Fd(?)- and NADH-dependent methylene-tetrahydrofolate reductase (MetFV-HdrABC-MvhD) in Moorella thermoacetica (2014), Fd- and NAD-dependent lactate dehydrogenase (LctBCD) in A. woodii (2015), Fd- and F₄₂₀H₂-dependent heterodisulfide reductase (HdrA2B2C2) in Methanosarcina acetivorans (2017), and Fd- and NADH-dependent ubiquinol reductase (FixABCX) in Azotobacter vinelandii (2017). The electron-bifurcating flavoprotein complexes known to date fall into four groups that have evolved independently, namely those containing EtfAB (CarED, LctCB, FixBA) with bound FAD, a NuoF homolog (HydB, HytB, or HylB) harboring FMN, NfnB with bound FAD, or HdrA harboring FAD. All these flavoproteins are cytoplasmic except for the membrane-associated protein FixABCX. The organisms—in which they have been found—are strictly anaerobic microorganisms except for the aerobe A. vinelandii. The electron-bifurcating complexes are involved in a variety of processes such as butyric acid fermentation, methanogenesis, acetogenesis, anaerobic lactate oxidation, dissimilatory sulfate reduction, anaerobic- dearomatization, nitrogen fixation, and CO₂ fixation. They contribute to energy conservation via the energy-converting ferredoxin: NAD⁺ reductase complex Rnf or the energy-converting ferredoxin-dependent hydrogenase complex Ech. This Review describes how this mechanism was discovered.

On the nature of organic and inorganic centers that bifurcate electrons, coupling exergonic and endergonic oxidation–reduction reactions

Bifurcating electrons to couple endergonic and exergonic electron-transfer reactions has been shown to have a key role in energy conserving redox enzymes.

H/D exchange mass spectrometry and statistical coupling analysis reveal a role for allostery in a ferredoxin-dependent bifurcating transhydrogenase catalytic cycle

Recent investigations into ferredoxin-dependent transhydrogenases, a class of enzymes responsible for electron transport, have highlighted the biological importance of flavin-based electron bifurcation (FBEB). FBEB generates biomolecules with very low reduction potential by coupling the oxidation of an electron donor with intermediate potential to the reduction of high and low potential molecules. Bifurcating systems can generate biomolecules with very low reduction potentials, such as reduced ferredoxin (Fd), from species such as NADPH. Metabolic systems that use bifurcation are more efficient and confer a competitive advantage for the organisms that harbor them. Structural models are now available for two NADH-dependent ferredoxin-NADP⁺ oxidoreductase (Nfn) complexes. These models, together with spectroscopic studies, have provided considerable insight into the catalytic process of FBEB. However, much about the mechanism and regulation of these multi-subunit proteins remains unclear. Using hydrogen/deuterium exchange mass spectrometry (HDX-MS) and statistical coupling analysis (SCA), we identified specific pathways of communication within the model FBEB system, Nfn from Pyrococus furiosus, under conditions at each step of the catalytic cycle. HDX-MS revealed evidence for allosteric coupling across protein subunits upon nucleotide and ferredoxin binding. SCA uncovered a network of co-evolving residues that can provide connectivity across the complex. Together, the HDX-MS and SCA data show that protein allostery occurs across the ensemble of iron‑sulfur cofactors and ligand binding sites using specific pathways that connect domains allowing them to function as dynamically coordinated units.