On the Natural History of Flavin-Based Electron Bifurcation

Electron bifurcation is here described as a special case of the continuum of electron transfer reactions accessible to two-electron redox compounds with redox cooperativity. We argue that electron bifurcation is foremost an electrochemical phenomenon based on (a) strongly inverted redox potentials of the individual redox transitions, (b) a high endergonicity of the first redox transition, and (c) an escapement-type mechanism rendering completion of the first electron transfer contingent on occurrence of the second one. This mechanism is proposed to govern both the traditional quinone-based and the newly discovered flavin-based versions of electron bifurcation. Conserved and variable aspects of the spatial arrangement of electron transfer partners in flavoenzymes are assayed by comparing the presently available 3D structures. A wide sample of flavoenzymes is analyzed with respect to conserved structural modules and three major structural groups are identified which serve as basic frames for the evolutionary construction of a plethora of flavin-containing redox enzymes. We argue that flavin-based and other types of electron bifurcation are of primordial importance to free energy conversion, the quintessential foundation of life, and discuss a plausible evolutionary ancestry of the mechanism.

Introduction

Life fundamentally depends on the free energy (ΔG) provided by the electrochemical disequilibrium between reduced, electron-donating and oxidized, electron accepting environmental substrates (see Schoepp-Cothenet et al., 2013). One might think that metabolic redox reactions must thus occur within the redox span imposed by these substrates. This clearly is not the case. The acidophilic iron oxidizers from the Acidithiobacillaceae for instance provide a rather extreme counterexample to such redox range limitations. These organisms use Fe²⁺ (with E_m >>0 mV) as environmental reductant but nonetheless manage to maintain their NAD and ferredoxin pools (featuring E_ms below -300 mV) partially reduced. The necessary up-conversion of redox energy in these species (Ingledew, 1982; Ferguson and Ingledew, 2008; reviewed in Nitschke and Bonnefoy, 2016) as well as in a multitude of other organisms facing the same problem is carried out with the help of the chemiosmotic membrane potential which serves to accumulate incremental ΔG from several individual redox reactions. While this scheme is often referred to as “driving electrons uphill,” a more appropriate view is that it uplifts the reducing power of these electrons so that they can flow downhill toward NAD and ferredoxins.

This kind of chemiosmotically driven “reverse electron transfer” was considered the general and unique mechanism to augment the reducing power of environmentally provided reducing equivalents to the levels demanded by metabolic needs. That is, until flavin-based electron bifurcation entered the scene. In 2008, the enzyme ETF (electron-transfer-flavoprotein) was reported to singlehandedly, i.e., in one redox reaction and without the implication of a chemiosmotic membrane system, be able to reduce electron acceptors having electrochemical potentials substantially more reducing than that of the electron donating substrate (Herrmann et al., 2008; Li et al., 2008). Reduction of the low redox potential acceptor, however, was observed to be dependent on the presence and concomitant reduction of a second substrate with much higher redox potential. Likewise, reduction of the high potential substrate was found to depend on the presence and reduction of the low potential acceptor (Chowdhury et al., 2014, 2015, 2016; Demmer et al., 2017). Enzyme-bound flavin, a two-electron redox compound, was proposed as the crucial redox cofactor permitting the coupling between an endergonic one-electron transfer reaction to the low potential acceptor and a (more strongly) exergonic one reducing the high potential acceptor (Chowdhury et al., 2016). The two electrons on the fully reduced flavin are thus thought to bifurcate toward energetically dissimilar pathways heaving one electron onto a more strongly reducing state while still obeying the thermodynamic imperative of a negative ΔG for the entire two-electron redox reaction (Buckel and Thauer, 2013). Since the first proposal of flavin-based electron bifurcation (Herrmann et al., 2008), a multitude of electron bifurcating enzymes containing flavins have been reported to operate in diverse types of redox metabolisms (for reviews see Peters et al., 2016; Buckel and Thauer, 2018a,b; as well as the most recent member of the family described in Duan et al., 2018).

The conceptual similarity of flavin-based electron bifurcation to the so-called Q_o-site reaction in the membrane-integral Rieske/cytb complexes (alias Complex III, bc₁, b₆f, etc.) was pointed out early on (Herrmann et al., 2008). The core functional reaction in Peter Mitchell’s Q-cycle mechanism (Mitchell, 1975, 1976) consists in the coupled redox transitions between a quinone, i.e., a two-electron redox compound, bound at the Q_o-site of the enzyme and a pair of one-electron redox centers one of which features a higher and the other one a lower electrochemical potential than the quinone (see also Wikström and Berden, 1972). The fundamental electrochemical and structural properties allowing the Q_o-site reaction to occur have been worked out over several decades (for comprehensive reviews, see Osyczka et al., 2004; Crofts et al., 2013) and part of this contribution will be dedicated to further explore common features of these two types of electron bifurcation to progress toward a comprehensive molecular and thermodynamic understanding of the flavin-based case (see also Nitschke and Russell, 2012; Peters et al., 2016; Buckel and Thauer, 2018a). At the end of this contribution we will discuss the potential of the electron bifurcation reaction to contribute to the emergence of life on earth.

The Peculiarities of Two-Electron Redox Compounds and the Concept of Redox Cooperativity

As worked out almost a century ago (Michaelis, 1932), certain heteronuclear aromatic two-electron redox compounds (and most prominently quinones and flavins) can exhibit astonishingly versatile redox properties. Under certain conditions they undergo two consecutive redox transitions between the oxidized (ox) and the one-electron- (semi-) reduced and between the 1-electron- and the 2-electron “fully” reduced (red) states. In this case, the redox midpoint potentials of the two individual transitions follow the intuitive rule that putting the second electron on the molecule will be “harder” (i.e., occur at a lower electrochemical potential) than adding the first one. However, under appropriate conditions in protic solvents, the first transition toward the semi-reduced state can render, due to redox-coupled protonation/hydrogen-bonding events (Lim et al., 2013; Tan et al., 2015), the second transition more favorable than the first one so that both electrons will come or go virtually together and the observed redox transition will follow an n = 2-type behavior¹. A redox transition where both electrons are transferred together is commonly called “cooperative.” More generally, the terms “positive” and “negative cooperativity” are used to distinguish the “both together” from the “one-by-one” redox regimes. Negative cooperativity means that the first reduction reaction leads to a decrease in the midpoint potential for the second reduction (i.e., corresponding to the “intuitive” situation mentioned above) and vice versa. As shown by Michaelis (1932), negatively and positively cooperative two-electron redox transitions form a continuum characterized by a smoothly varying spacing ΔE between the redox midpoint potentials of the individual transitions ΔE = E_ox/sr – E_sr/red.

While of course strong negative cooperativity implies that in equilibrium redox titrations the semi-reduced form is stable over a wide range of ambient potentials, increasingly positive cooperativity will be accompanied by a decreasing stability constant K_S² of this semi-reduced state translating into lower percentages of the semi-reduced state observable at equilibrium. The relationship between K_S and ΔE is most clearly illustrated by the so-called “redox seesaw” representation (Mitchell, 1976)³. Figure 1 illustrates the relationship between the redox midpoint potentials of the individual redox transitions and K_S. For K_S < 1 (logK_S < 0), E₂ is more positive than E₁. The logK_S < 0 domain is therefore frequently referred to as the region of “inverted” or “crossed-over” potentials.

FIGURE 1

FIGURE 1. Redox seesaw-type representation of the correlation between the electron transfer function of selected flavoenzymes and the redox cooperativity of their flavin cofactors. The redox midpoint potentials of the first (oxidized to one-electron reduced, red dashes) and the second redox transition (half-reduced to two-electron reduced, blue dashes) is shown as a function of the semiquinone stability constant K_S = [semi-reduced]²/[ox]x[red]. The maximal occupancy of the flavosemiquinone as would be observed in equilibrium redox titrations on the respective flavins is indicated as a violet dashed curve. Note that in this and all following figures more negative potentials are toward the upper part of the graphs.

The Flavin, a Redox Cofactor for All Seasons

As indicated in Figure 1, the specific binding sites in enzymes confer very diverse degrees of cooperativity to the flavin cofactors. The resulting differences in redox behavior allow these two-electron redox compounds to be used by life for a plethora of different purposes. The strictly one-by-one-type flavodoxins thus play the role of one-electron transport proteins with, in general, only one of the two redox transitions being physiologically relevant. A well-studied example are the flavodoxins found in certain cyanobacteria under iron-limiting growth conditions. In these systems, the semi-reduced/reduced transition serves to accept electrons from Photosystem I and to distribute them to the physiological electron acceptors just as [2Fe-2S]-ferredoxins would do under iron-replete conditions (Sétif, 2001). By contrast, flavin cofactors mediating two-electron redox reactions between substrates with strong positive redox-cooperativity, such as NADH and quinone, function optimally if they feature strongly destabilized semi-reduced states, i.e., if they are far into the region of inverted potentials shown in Figure 1. The flavoenzyme NDH-2 for example, a single-polypeptide NADH:quinone oxidoreductase associated to the n- (i.e., negatively charged) side of bioenergetic membranes mediates the two-electron transfer involved in the oxidation of NADH and the reduction of quinone. It does not involve observable semi-reduced states of the flavin (de Vries and Grivell, 1988). Finally, two-electron redox compounds with redox properties in the vicinity of logK_S = 0 in Figure 1 strike a compromise between strictly one- and strictly two-electron transitions and hence function as “gates” mediating between 2-electron and 1-electron segments of electron transfer chains. A prominent example is provided by FNR (ferredoxin NAD(P)⁺ reductase) which deals with both one-electron (ferredoxin) and two-electron (NAD) redox partners (Beirns and Wang, 1972; Batie and Kamin, 1981; Sánchez-Azqueta et al., 2014). In Figure 1 we distinguish the mentioned types of electron transfer reactions by color-coding “single-electron” reactions in blue, “gating”-type ones in violet and “both-together” ones in orange/red. So where does electron bifurcation fit into this scheme?

A Common Scheme for Quinone- and Flavin-Based Electron Bifurcation

As stated above, flavin-based electron bifurcation likely relies on the same thermodynamic concept as its quinone-based counterpart observed at the Q_o-site of bc₁-complexes and related enzymes. Research on the Q_o-site reaction is ahead of the flavin case by several decades and many years of unsuccessful attempts to observe the EPR signature of the semi-reduced state of the quinone in the Q_o-site of the Rieske/cytb complexes under conditions of equilibrium redox titrations suggested a very low K_S value. This led to the conclusion that the Q_o-site quinone must undergo a two-electron transition with strong positive cooperativity during enzyme turnover in line with Peter Mitchell’s original hypothesis (Mitchell, 1976). More recent quantitative determinations yielded indeed astonishingly low values in the region of logK_S = -14 (Zhang et al., 2007; Crofts et al., 2013), that is, even lower than the semiquinone stability in bulk aqueous solutions (logK_S = -9 to -10.5 according to the chemical type of quinone, Hauska et al., 1983; Song and Buettner, 2010). The midpoint potentials for the individual one-electron redox transitions corresponding to such low K_S values yield a redox landscape for electron transfer through the Rieske/cytb complexes as depicted in Figure 2A. Remarkably, the general layout of this redox landscape appears to be conserved in all types of Rieske/cytb complexes analyzed so far (Bergdoll et al., 2016). The most striking feature of this redox scheme is the strongly uphill electron transfer from the quinol/semiquinone couple to the Rieske [2Fe-2S] cluster. This endergonic reaction upends the traditional view which referred to reduction of the Rieske cluster as the exergonic branch of the bifurcating Q_o-site reaction and the reduction of heme b_L as the endergonic one. However, the strong endergonicity (ΔE_m ∼ 200 mV) of electron transfer from the quinol to the Rieske cluster provides a straightforward rationale for the strong coupling observed for the two redox steps in oxidation of the quinol. While the electron tunneling back and forth between the quinol and the Rieske [Fe-S]-cluster resides more than 99.9% of the time on the quinone, the full two-electron oxidation reaction becomes (slightly) exergonic if oxidation of the semiquinone by the second electron acceptor (heme b_L) can also occur. The driving force for oxidation of the semiquinone therefore pulls the thermodynamically impeded “first” electron transfer over the activation barrier. This semiquinone, however is of course only present and thus available for oxidation while this “first” electron resides on the [2Fe-2S]-cluster. The probability of this configuration thus determines the kinetics of the coupled two-electron reaction. The combined reaction pattern has been aptly described as an “escapement” mechanism (Crofts et al., 2013) wherein the entire reaction can only occur if the two individual steps are both occurring. The coupling between the individual reactions is imposed by the specific electrochemical properties (i.e., the strongly inverted individual redox midpoint potentials) conferred to the quinone by the Q_o-binding site. The entire reaction is of course reversible and flavins (or quinones) with crossed-over potentials can be reduced by two electrons from distinct donors with the first reduction step at a lower redox potential than the second one. This process is referred to as electron confurcation and systems proposed to rely on electron confurcation encompass enzymes involved in sulfur metabolism such as QmoA and QmoB (Ramos et al., 2012; Grein et al., 2013) as well as certain hydrogenases (Schut and Adams, 2009; Wang et al., 2013b; Zheng et al., 2014).

FIGURE 2

FIGURE 2. Electrochemical landscapes for the electron bifurcation reactions in (A) the Rieske/cytb complexes (the shown potentials are for the ubiquinone-oxidizing enzymes in mitochondria and alpha-proteobacteria; for enzymes operating with other types of quinones, see Bergdoll et al., 2016), (B) ETF, and (C) Nfn-1.

Being a relatively recently discovered phenomenon, it is not surprising that knowledge of the electrochemical details regarding flavin-based electron bifurcation is still limited. We have depicted the two systems for which the relevant redox parameters have been obtained, i.e., the enzymes NADH-dependent ferredoxin:NADP⁺ oxidoreductase (Nfn-1) and the bifurcating electron-transfer-flavoprotein (ETF) in Figures 2B,C next to that detailed above for the Rieske/cytb complex. It is noteworthy that the exact value for the K_S of the bifurcating flavin in ETF has not been determined so far but has been assumed (Chowdhury et al., 2016; Demmer et al., 2017) in analogy to the Nfn-1 results (Lubner et al., 2017) or to the value obtained for the quinone in the Q_o-site reaction (Osyczka et al., 2004; Crofts et al., 2013). Determination of K_S for the bifurcating flavin in Nfn-1 is itself based on an indirect method involving the redox midpoint potential of the [4Fe-4S]-cluster oxidizing the photochemically generated flavo-semiquinone to the fully oxidized state. Redox midpoint potentials of the two [4Fe-4S]-clusters in Nfn-1 have been obtained via equilibrium redox titrations (Lubner et al., 2017) and the E_m of the cluster accepting the electron from the flavin under conditions of the photochemical experiment are likely to differ from equilibrium titration values as a result of electrostatic interactions in the latter approach. These experimental uncertainties notwithstanding, the bifurcating flavin in ETF can be safely assumed to feature strongly inverted individual redox midpoint potentials since strict two-electron behavior has been observed in equilibrium redox titrations on the enzyme from Megasphaera elsdenii (Sato et al., 2013) and on the FixAB-homolog in Rhodopseudomonas palustris (Duan et al., 2018).

The comparison of the three redox landscapes depicted in Figure 2 reveals an intriguing unicity for the overall thermodynamics of the different types of electron bifurcation. In all cases, extracting the first electron from the fully reduced two-electron redox compound turns out to be a strongly endergonic step which needs the second oxidation step to occur in order to become stabilized. Similarly, in all cases the full oxidation of the two-electron redox compound is slightly exergonic but remains very close to thermodynamic equilibrium.

The picture emerging from the comparison of the admittedly small sample of sufficiently well-characterized flavin-based systems (Figures 2B,C) to the quinone-based Q_o-site reaction (Figure 2A) is thus that of a common mechanistic principle based on strongly crossed-over electrochemical potentials for the individual redox transitions of the two-electron redox compounds and a strong endergonicity of the first oxidation step yielding all the ingredients required for the escapement-type mechanism discussed above.

Intriguingly, in all quinone- and flavin-based electron bifurcating enzymes for which 3D-information is available, a large-scale conformational movement of a protein domain, in most cases that carrying the electron acceptor for the first electron, has been observed or inferred (Zhang et al., 1998; Demmer et al., 2017; Lubner et al., 2017; Wagner et al., 2017). The above described mechanism per se does not present an obvious need for such conformational changes. These conformational movements may serve to optimize electron transfer rates and minimize efficiency-decreasing back- or deleterious side-reactions. A more comprehensive characterization of the flavin-based systems may in the future allow a deeper understanding of the roles of conformational changes in the phenomenon of electron bifurcation.

It is worth noting that electron bifurcation obviously occurs in roughly the same region of inverted redox midpoint potentials as do “both together”-type 2-electron transfer reactions such as for example in NDH-2 (Figure 1). This is not a coincidence since two-electron redox compounds react extremely sluggishly with a one-electron donor or acceptor because of their strongly inverted redox midpoint potential. The K_S for the NADH/NAD⁺ couple has been determined at 10^-20 (Anderson, 1980) yielding redox midpoint potentials for the first oxidation of +340 mV and for the second one of -940 mV. Electron acceptors for NADH typically have redox midpoint potentials well below 0 mV and the first oxidation step is therefore steeply uphill. This endergonicity is even dramatically enhanced if the acceptor is a two-electron redox compound which has itself strongly inverted redox midpoint potentials since its 1^st reduction potential E₁ is substantially lower than its average two-electron potential E_m (see Figure 3). The transfer of the second electron (from -940 mV), however, will occur toward the 2^nd redox transition of the acceptor molecule which is, in the case of inverted redox midpoint potentials, higher than the average two-electron potential E_m. This step will therefore be strongly exergonic and render the combined reaction (gray arrow in Figure 3) favorable whereas in the presence of only single-electron acceptors the electron will remain on the NADH molecule. “Both-together” reactions can therefore be regarded as processes which bifurcate the individual electrons thermodynamically, that is, toward strongly dissimilar reduction potentials, but not spatially. The two redox transitions of the two-electron acceptor with dissimilar (crossed-over) redox midpoint potentials play the same role as the two electrochemically AND spatially distinct electron acceptors in electron bifurcation for rendering the total reaction thermodynamically downhill.

FIGURE 3

FIGURE 3. Electrochemical landscape for the two-electron transfer reaction between NADH and FAD in the framework of the escapement-type mechanism discussed in the text.

The Importance of Electron Bifurcation Based on Strongly Inverted Redox Potentials to Life

The thermodynamic outcome of electron bifurcation as described above fundamentally differs from that of the vast majority of enzyme-mediated catalyzes. While the latter perform reactions which in general uniformly decrease the free energy of the system (a different way of saying that they are exergonic), electron bifurcation certainly is exergonic overall but it generates a product (one out of two) which is even farther from thermodynamic equilibrium than the educt. In ETF for instance, NADH can reduce flavins exergonically in a two-electron reduction reaction, but the ferredoxin reduced as a result of ensuing flavin-based electron bifurcation is substantially more reducing than the initial reductant, NADH. Extending this reasoning to the thermodynamics of life in its redox context, i.e., its source of environmental free energy (Schoepp-Cothenet et al., 2013), an organism growing for example on molecular hydrogen as electron donor augments, via electron bifurcation, the reducing power available to drive metabolic reactions far above the reducing power directly furnished by its habitat. This increase in electrochemical disequilibrium corresponds to a decrease in entropy in the sub-group of redox reactions involving the so-produced reduced ferredoxin molecules. As a pay-back, the other electron emerging from the electron bifurcation reaction features a higher redox midpoint potential resulting in a lower redox disequilibrium between this electron and its acceptor compounds. Electron bifurcations can therefore alternatively be regarded as free energy converters which drive an increase of order in the living system over that of the environment via a concomitant generation of “heat,” i.e., the exergonic, entropy-increasing redox reaction.

The essence of electron bifurcation as encountered in ETF, Nfn-1, or the Q_o-site reaction thus lies in its ability to increase redox energy above that provided by the environment. As detailed above, this necessarily involves strongly inverted redox potentials. Recent articles addressed the question whether electron bifurcation was conceivable without involving strong positive redox cooperativity (Hoben et al., 2017; Zhang et al., 2017). While we do not dispute the general possibility of the proposed scenarios based on uncrossed potentials (Zhang et al., 2017) we would hold that, for the enzymes considered in these articles, such scenarios miss the point since the respective types of electron bifurcation would be useless to biology. In these scenarios the bifurcating cofactor needs to be fully reduced by upstream electron donors with redox potentials commensurate or lower than that of the low potential electron acceptor of the bifurcating reaction. The entropic benefit for the organism would therefore be nil in scenarios based on uncrossed potentials. However, in certain cases a kinetic rather than thermodynamic benefit may be reaped from only spatially bifurcating electrons as we will argue happens in the case of the NuoE/F system of Complex I discussed below. In such circumstances strong inversion of redox midpoint potentials may be dispensable.

What About Geometry?

Kinetics of bifurcation reactions, like all electron transfer kinetics, depend on several parameters implemented in Marcus’ theory of outer shell electron transfer (Moser et al., 1992). Beside ΔG and reorganization energy, distances between co-factors strongly influence the kinetic details of such electron transfer reactions. For the majority of flavin-based systems crucial information necessary to apply this rigorous formalism is still missing today or is incomplete. For example, in the 3D structures of both ETF and heterodisulfide reductase (Hdr), a number of electron transfer distances to or from the flavin are too long to permit biologically relevant rate constants and large-scale domain movements must therefore be invoked (Demmer et al., 2017; Wagner et al., 2017). This problem is less striking for Nfn-1, although the distance between the bifurcating flavin and the [2Fe-2S]-cluster on the high potential branch at about 14 Å may be close to the upper limit for physiologically sensible electron transfer (Lubner et al., 2017). Still, the attribution of redox midpoint potentials to individual [4Fe-4S]-centers on the low potential branch in Nfn-1 so far is arbitrary and not substantiated by experimental evidence. As mentioned before, these redox midpoint potentials have furthermore been determined through equilibrium titrations and may be lower than “operating” potentials due to electrostatic repulsion effects.

As a first approximation, we therefore compared the available 3D-structures in order to assess whether general properties of electron transfer pathways between the bifurcating flavins and their electron donors or acceptors can be found. An exhaustive representation of this comparison is presented as Supplementary Figure S1 and a few representative examples are shown in Figure 4. In this figure, the structural arrangements of electron donors and acceptors with respect to the flavin are depicted from two different view angles for the enzymes Nfn-1, Hdr, NDH-2, and SQR (Sulfide quinone reductase). Since in Hdr the presumable electron acceptor on the low potential branch is one of two roughly equidistant [4Fe-4S]-centers (Wagner et al., 2017), we have shown both of these clusters in Figure 4. Furthermore, substantial movement of the domain binding these two centers is likely to occur (Wagner et al., 2017) and the indicated direction of electron transfer is therefore somewhat uncertain. The quinone acceptor was not resolved in the published 3D-structure of NDH-2 and is therefore only indicated for the SQR enzyme. Vice versa, the electron donating NADH was only identified in NDH-2.

FIGURE 4

FIGURE 4. Spatial arrangement with respect to the flavin of electron donating and accepting cofactors in four different flavoenzymes, two of which perform electron bifurcation (Nfn-1 and Hdr), while the other two transfer both electrons together from/to their redox partners. The sideview in (A) renders the cofactors as seen from within a plane perpendicular to the flavin macrocycle (indicated in gray) while in (B) the same cofactors are seen from above this plane. LP and HP stand for low and high potential electron acceptors, respectively.

The shown structural comparison yields the following conclusions. (i) The donor/acceptor substrates NADH/NAD⁺ are accommodated co-planar to the isoalloxazine ring system on the side opposite to the bending direction of the tail leading to the adenine nucleotide, i.e., uncorrelated to the Si- or the Re-side (see Supplementary Figure S1). (ii) Other one- and two-electron redox partners (iron–sulfur centers, quinone) are clustered in or near a plane perpendicular to the isoalloxazine ring system and containing its long axis (Figure 4A and Supplementary Figure S1). (iii) Most single-electron transfer reactions from or to the flavin seem to proceed from the extremities of the conjugated isoalloxazine moiety. Notable exceptions to this rule are provided by Hdr (Figure 4, blue cofactors) or Complex I (see below). In all available structures, however, the two electrons involved in bifurcating/confurcating electron transfer reactions take diametrically opposite directions to/from their acceptor/donor molecules (Figure 4B).

The comparison shown in Figure 4 thus indicates that overall electron transfer geometry does not seem to play a significant role in the type of electron transfer the flavin is involved in (gating, both-together, bifurcating). This comforts our view that electron bifurcation is foremost an electrochemical phenomenon.

The Flavoprotein Construction Kit; Toward a Basic Set of Structural Motifs

As illustrated in Figure 1, the redox properties of diverse flavoproteins cover the whole range allowed by the electrochemistry of two-electron redox compounds and they correspondingly perform as well “one-by-one” as “one/two gating,” “both together,” and “bi- and confurcating” electron transfer reactions. While the vast majority of flavin-binding domains in proteins are made up from Rossmann-fold-type (Rao and Rossman, 1973) basic units, substantially differing variations to this general theme have been observed. We have tried to classify a number of flavin-binding domains relevant to bioenergetic redox reactions into structural motifs based on their 3D structures and have identified three distinguishable architectures which we denote as (a) the NDH-2-type, (b) the FNR-type, and (c) the flavodoxin-type folds (see Figure 5). During our searches of 3D databanks we came across cases which either seemed to be more distantly related to one of these families (Supplementary Figure S2a, right hand part and S2a′), to result from convergent evolution or to represent entirely different folds (Supplementary Figures S2d–f). We aren’t familiar with these enzymes which are at the border of the field of bioenergetics but have nevertheless included these additional cases in the Supplementary Figure S2 and listed their electrochemical properties in Supplementary Table S1.

FIGURE 5

FIGURE 5. Structural comparison of selected flavoenzymes or their flavin-containing domains highlighting conserved folds as discussed in the text. The color-coding indicates the mode of electron transfer reactions as defined in Figure 1. A more extended version of this figure, containing substantially more cases, can be found as Supplementary Figure S2. Detailed electrochemical parameters and relevant references are listed in Supplementary Table S1.

The NDH-2- (Figure 5A) and the FNR-type (Figure 5B) motifs fall into distinct structural domains with the flavin sandwiched in between these domains. In Figures 5A,B, the N-terminal domains are at the bottom of the depicted structures followed by domains colored in a lighter shade and represented in the upper part of the individual structures. The NDH-2 folds feature a third domain (marked in bright orange in Figure 5) that reaches back down toward the N-terminal domain. The C-terminal domain of the FNR fold (Figure 5B, light pink) strongly resembles the flavodoxin fold (Figure 5C). The flavodoxin fold can also be identified in pyruvate oxidases (Supplementary Figure S2g), in NuoF, a subunit of Complex I (Supplementary Figure S2h), and in flavorubredoxin reductase (Supplementary Figure S2i). Pyruvate oxidases and flavorubredoxin reductase are composed of two flavodoxin-like domains. NuoF features the flavodoxin domain at its N-terminus followed by a β-sheet and an α-helical domain, whereas in FNR the flavodoxin domain constitutes the C-terminal part of the protein and is preceded by an entirely β-sheet-containing N-terminal domain. The flavodoxin fold therefore appears to have been fused several times independently to various other domains in order to build the above discussed enzymes in a nice illustration of the construction kit principle used by nature to build novel enzymes (Baymann et al., 2003). Sandwiching the flavin moiety within the contact region of two/three domains results in a relatively uniform accommodation of the isoalloxazine ring system among the different systems. The peripheral location of the flavin in flavodoxin-type folds, by contrast, allows for substantial variability in arrangement of the cofactor in this structure (see Figure 5C).

Several cases, such as for example HdrA, feature extensive sequence insertions within the flavin-binding units resulting in additional, sometimes cofactor-binding domains while nevertheless maintaining the common overall structure of the flavin-domain. While these somewhat freakish insertions are relatively rare, a sequential lining up of separate domains containing extra cofactors or that are even devoid of such cofactors in a single polypeptide chain frequently occurs. Thus, for the sake of clarity, the structures shown in Figure 4 feature only the relevant flavin-domains omitting the (sometimes substantial) remainder of the entire protein.

The NuoF case is part of the folds where an extra cofactor, a [4Fe-4S]-cluster is harbored by an additional domain and is positioned at 8 Å distance from the flavin. A [2Fe-2S]-center is furthermore located at a distance of 13 Å on the neighboring NuoE subunit. This structural arrangement thus looks perfectly suited for electron bifurcation as in the cases mentioned above. However, the electrochemical properties determined for the involved cofactors (Gnandt et al., 2016) shows that the NuoE/F system doesn’t follow the scheme depicted in Figure 2. In the following we will discuss the fold and the electron transfer properties of this intriguing case in more detail.

The Puzzling Case of the NuoE/F-Family

NuoE and NuoF (also denoted Nqo1 and Nqo2) form a compact structural module (Sazanov, 2016) containing a flavin, a [4Fe-4S]- (“N3”) and a [2Fe-2S]-cluster (“N1a”). 3D structures with NAD in close proximity to the flavin have been obtained (Figure 6A) and the module is thought to negotiate between 2-electron oxidation of NADH and 1-electron reduction of the chain of iron–sulfur centers which leads toward the site of quinone reduction. The flavin in this case therefore acts as a gate between 2- and 1-electron redox compounds. Cluster N3 is the entry point to the electron wire to the quinone-binding site while cluster N1a appears as a dead end roughly opposite N3 with respect to the flavin moiety (Figure 6A) and has no obvious connection to other electron transfer chains. In Complex I from Escherichia (E.) coli and Aquifex (A.) aeolicus, cluster N1a features redox potentials of about -250 mV and is observed to go reduced upon addition of NADH to the enzyme. In the mitochondrial and the Thermus thermophilus enzymes it titrates close to -400 mV and does not seem to be reduced by NADH (Gnandt et al., 2016). Cluster N3 titrates in the vicinity of -250 mV in all systems (Gnandt et al., 2016). The flavosemiquinone was studied by EPR and a small but measurable logK_S of 1.34 was determined (Sled et al., 1994). This translates into crossed-over potentials with ΔE ∼-80 mV and individual potentials for the first and second oxidation reactions of the flavin at -336 and -414 mV, respectively. The electrochemical landscape resulting from these redox parameters is depicted in Figure 6B.

FIGURE 6

FIGURE 6. Spatial arrangement (A) and electrochemical landscape (B) of the NuoE/F-module in Complex I. The spatial bifurcation reaction is indicated in (A) by the two continuous arrows while the subsequent electron transfer events down the chain of iron–sulfur centers and toward the quinone (not shown) are marked by dashed and dotted arrows. In (B) two different redox potentials are depicted for cluster N1a, corresponding to the values reported for the mitochondrial and Aquifex aeolicus enzyme (–250 mV) and for the Thermus thermophilus enzyme (–400 mV).

Obviously, this system lacks the endergonicity for the first electron leaving the flavin characteristic for canonical electron bifurcation, since it can flow downhill to cluster N3 (or also cluster N1a for the systems with E_m = -250 mV) with a drop in potential of almost 100 mV. The second electron can then exergonically reduce the remaining cluster or, for the case where E_m of N1a is around -400 mV, equilibrate between the flavosemiquinone and cluster N1a (Figure 6B). The two individual flavin oxidation events therefore aren’t coupled via the described escapement mechanism and no up-conversion of redox energy will take place. In agreement with these electrochemical properties, NADH-oxidation activity by quinone is only marginally affected when the redox potential of cluster N1a is mutationally shifted from -250 to -400 mV (Gnandt et al., 2016) and thus eliminated as an efficient electron acceptor. It therefore seems that the comparatively weak potential inversion of the flavin in the NuoE/F system allows for successive electron transfer down the wire through cluster N3 as well as for simultaneous reduction of the two iron–sulfur clusters N1a and N3. In the latter case the electrons bifurcate spatially but not thermodynamically. To put this in mechanistic terms, the gear-transmission-type coupling of the bifurcation reaction appears to have been replaced by a slip-clutch mechanism.

Intriguingly, the so-called “electron bifurcating hydrogenases” contain two subunits with substantial sequence similarity to NuoE and NuoF. As is in their name, these hydrogenases have been reported to bifurcate electrons derived from H₂ oxidation toward reduction of NAD⁺ and ferredoxin. Enzymes of this type have been reported from several species, such as Acetobacterium woodii (Schuchmann and Müller, 2012), Thermotoga maritima (Schut and Adams, 2009), Moorella thermoacetica (Wang et al., 2013b), Clostridium autoethanogenum (Wang et al., 2013a), and Ruminococcus albus (Zheng et al., 2014). They all contain, in addition to the NuoE/F orthologs, an [Fe-Fe]-type hydrogenase subunit related to the monomeric [Fe-Fe]-hydrogenase of Clostridium pasteurianum. A subpopulation of these electron bifurcating hydrogenases furthermore contains a thioredoxin-like [Fe-S]-containing subunit which, however, does not seem essential for electron bifurcation (Poudel et al., 2016).

Like for many of the enzymes dealt with in this article, electron transfer direction is readily reversed in these hydrogenases and some of them thus also perform redox confurcation [e.g., the enzymes from Moorella thermoacetica, Clostridium autoethanogenum, and Ruminococcus albus (Wang et al., 2013a,b; Zheng et al., 2014)]. Unfortunately, no 3D structure has so far been solved for any of these enzymes. Based on the detection (Verhagen et al., 1999; Schut and Adams, 2009; Schuchmann and Müller, 2012; Wang et al., 2013b) of a single flavin-binding sequence motif (Hanukoglu and Gutfinger, 1989) in electron bifurcating hydrogenases, initially the presence of only one flavin molecule was proposed. The respective motif in the hydrogenases is at a sequence position equivalent to that of the motif in NuoE in our multiple sequence alignments. Putative ligands for two iron–sulfur centers corresponding to N1a and N3 in Complex I are also found in these alignments. More recently, Buckel and Thauer (2013) have postulated that “two FMN should be bound, one being required for NADH dehydrogenation and a second one for flavin-based electron bifurcation.” This reasoning makes perfect sense since two electrons need to be accumulated by a gating flavin to reduce NAD⁺ in a two-electron reaction, a task impossible to be achieved by the bifurcating flavin which sends its two electrons in opposite directions. The substantial sequence conservation between the two relevant subunits carrying the flavin and the N1a and N3 clusters and their homologs in the hydrogenases, together with the absence of any flagrant second flavin binding motif in both subunits, does not argue for the presence of a second flavin. This conundrum recently prompted the proposal that the bifurcating site in bifurcating hydrogenases may actually be at the active-site H-cluster rather than at the flavin (Peters et al., 2018).

Since in our view electron bifurcation arises from cooperative two-electron chemistry, the phenomenon likely isn’t restricted to redox-active organic molecules such as quinones and flavins but may be expected to also occur on transition metals or clusters thereof provided they feature three or more redox states and sufficiently strong positive redox cooperativity between the respective redox transitions. The case of molybdenum, for which all these properties have been observed empirically, will be discussed below. While the redox spacing between the two transitions (ferrous/ferric and ferric/ferryl) of mononuclear iron is generally found to be discouragingly wide (almost 1 V) and redox potential inversion hence unlikely, the electrochemical properties of the two-iron containing H-cluster may be more conducive to crossing of the individual potentials. However, if the 3D-structure of the [Fe-Fe]-hydrogenase from Clostridium pasteurianum (with almost 40% sequence similarity to the HydA subunit in the bifurcating enzyme from Thermotoga maritima) provides any guidance, we have difficulties seeing a potential second entry pathway (complementing the linear chain of [4Fe-4S]-clusters for electron confurcation) toward the H-cluster which is deeply buried within its parent protein subunit. Nevertheless, the proposed shift of the bifurcating site from the flavin to the H-cluster represents a novel approach to rationalize the enigma of bifurcating hydrogenases and calls for experimental scrutiny.

Evolutionary Use of the Construction Kit

The observation that the overwhelming majority of the flavoenzymes involved in bioenergetic electron transfer seems to be put together from a very limited set of flavin-binding domains together with a host of other elementary, often redox-cofactor containing, building components also found in many other enzymes once more emphasizes the construction-kit principle employed by evolution in the generation of novel bioenergetic enzymes (Friedrich and Weiss, 1997; Baymann et al., 2003; Grein et al., 2013; Grimaldi et al., 2013; Magalon et al., 2016). Rather than evolving novel protein architectures, nature appears to mainly draw from a pre-existing reservoir of structural motifs. Interestingly, these motifs often, but not always, are correlated to specific functional features. For example, enzymes employing the so-called Molybdo/tungsto-bisPGD building component generally catalyze redox conversions of two-electron redox compounds present in the environment (Grimaldi et al., 2013), the membrane-integral diheme cytochromes interconnect lipo-soluble and water-soluble electron transfer chains (Schoepp-Cothenet et al., 2013) and tetracubane iron–sulfur proteins serve as electron wires (for further examples, see Baymann et al., 2003). This begs the question whether the differing structural motifs are correlated to the different positively or negatively cooperative two-electron transfer properties discussed above and indicated in Figure 1. In Supplementary Table S1, we have compiled electrochemical properties of flavin-cofactors in respective enzymes as reported in the literature.

The color-coding employed in Figure 5 correlates the depicted structural modules to electron transfer function (that is, bifurcating and both-together in red/orange, gating in violet and one-by-one-type in blue). A more complete representation can be found in Supplementary Figure S2 and Supplementary Table S1. A clear trend, albeit not a strict rule, emerges from this correlation which admittedly may be biased by limitations in the currently available set of 3D structures. So far, members of the NDH-2 family appear to predominantly perform either both-together or bifurcating electron transfer and accordingly feature strongly inverted potentials (in most cases evidenced by the absence of an observable flavo-semiquinone). The depicted cases structurally resembling FNR all mediate electron transfer between 2-electron and 1-electron redox centers. The group of flavodoxin folds encompasses flavoproteins featuring all types of electron transfer properties (both together, bifurcating, gating and single-electron).

The fact that some correlation between protein structure and electron transfer properties of the flavin cofactor can be observed is astonishing in the light of results obtained on FNR from the cyanobacterium Anabaena. Very subtle modifications of the hydrogen-bonding network involving the N₅-nitrogen of the flavin moiety were found to induce drastic alterations in the relative redox midpoint potentials and to shift the individual potentials by about 100 mV further into the inverted region (Sánchez-Azqueta et al., 2014) which corresponds to a decrease of K_S by roughly two orders of magnitude. Given this subtle dependence of the strength of redox cooperativity on details of the hydrogen-bonding network surrounding the flavin, evolutionary changes likely would be able to easily redox-tune any of the discussed structural motifs to meet all kinds of electrochemical requirements. They indeed do so as for example in the enzyme dihydropyrimidine dehydrogenase (DPD) which features the canonical fold of the NDH-2 family but a mode of electron transfer falling into the realm of gating mechanisms. Supplementary Figure S1 shows several further cases which feature similar folds but strongly dissimilar electron transfer properties. In many cases, however, evolution appears to prefer recruiting for a specific function the protein fold which already performed the respective task earlier. This observed tendency to use ready-made functional modules rather than de novo engineered protein units nicely illustrates François Jacob’s notion of evolution by tinkering (Jacob, 1977).

Two or three distinct flavin-containing folds can now be added to the construction kit for electron transfer enzymes, depending on whether we count the FNR fold as a family of its own right or as a complexified version of the structurally simpler flavodoxin motif. Conserved 3D structures are frequently interpreted to suggest common evolutionary origins and it was for this reason that we refer to groups of enzymes exhibiting related folds as “families” in this contribution. Consequently, the protein or gene sequences of the individual proteins in the families depicted in Figure 3 may contain the imprints of their evolutionary history (Zuckerkandl and Pauling, 1965). Reconstruction of the phylogenetic tree of such a family may, under certain circumstances, allow the inference of evolutionary ancestry of the fold and consequently of its related function. This approach has been applied to Rieske/cytb complexes which harbor the quinone-based Q_o-site electron bifurcation reaction. Evolutionary roots reaching back to the last universal common ancestor (LUCA) of living beings were indeed derived from these analyses (Schütz et al., 2000; Lebrun et al., 2006; Ducluzeau et al., 2009; Ten Brink et al., 2013; Kao and Hunte, 2014). Several of the flavin enzyme families discussed in this article, such as ETF (Garcia Costas et al., 2017), Hdr or NDH-2 (Ojha et al., 2007; Marreiros et al., 2016), are widely distributed over both Bacteria and Archaea and their presence as early as in LUCA is therefore imaginable. A recently published comprehensive analysis of the phylogenetic distribution of ETF reconstructed an unrooted tree based on concatenated sequences (Garcia Costas et al., 2017). Although the lack of a reliable root precludes conclusions on pre- or post-LUCA origins of the observed clades of paralogs, the reported phylogeny does not exclude a pre-LUCA emergence of the entire enzyme superfamily. The Hdr enzyme from methanogens, for which a 3D structure has recently been obtained (Wagner et al., 2017), actually is part of a large superfamily with representatives in many organisms from both the archaeal and the bacterial domains (Grein et al., 2013). Hdr-related enzymes in the form of the Flx-HdrABC complex have notably been postulated to participate in sulfate-reducing bioenergetic pathways and to reduce ferredoxin and disulfides via electron bifurcation in Desulfovibrio species (Ramos et al., 2015). Hdr-related genes have furthermore been identified in several sulfur-oxidizing prokaryotes and a corresponding enzyme has been purified from the chemolithotrophic sulfur-oxidizer Aquifex aeolicus (Boughanemi et al., 2016). The phylogenetic relationships between all these Hdr-like entities remain so far only poorly understood. In general, phylogenetic approaches to the young field of flavin-based electron bifurcation are in their infancies and carry the potential for exciting future discoveries. In the absence of reliable top–down results, that is, conclusions derived from molecular phylogenies of enzymes involved in flavin-based electron bifurcation, we will in the following try to address the problem of the mechanism’s ancestry via plausibility arguments based on thermodynamics and palaeogeochemistry.

On the Evolutionary Emergence of Electron Bifurcation

As argued already more than half a century ago (Schrödinger, 1944), the crucial property defining life is its ability to convert environmental disequilibria into intracellular disequilibria. This argument has since been elaborated to great detail (Cottrell, 1979; Branscomb and Russell, 2013, 2018; Branscomb et al., 2017). Frequently cited examples for extant biological processes performing this fundamental task comprise ATP synthases and the redox-driven ion pumps of bioenergetic systems. The general principle of action in these systems is that endergonic and exergonic partial reactions are strictly coupled so that the overall reaction, while being exergonic, is contingent upon the actual proceeding of the endergonic partial reaction (Branscomb et al., 2017), just as described above for electron bifurcation. Given the outstanding importance of electron bifurcation as an up-converter of redox disequilibria and therefore of a local increase in order, we consider that life could have taken advantage of this process from its very beginnings. Furthermore, since electron bifurcation relies on the redox properties of specific molecules rather than on the specific action of dedicated enzymes, the general mechanism may well have operated in nascent life prior to the emergence of catalysis based on large enzymatic machineries.

Free energy converting reactions performed by transition metals and clusters thereof have been discussed previously as fundamental processes in nascent life. Mineral-harbored metal clusters lend themselves as first catalysts since they are abundant in the geosphere and have versatile electron transfer and ligand binding properties (Nitschke et al., 2013). Electron bifurcation, however, has been discussed in this contribution as a mechanism based on organic molecules (quinones and flavins; however, see the final paragraph of the NuoEF-section) which may at first sight seem to exclude their contribution in nascent life.

We will in the following advocate that firstly, transition metals themselves can perform bifurcating electron transfer reactions just as quinones and flavins do and secondly, that flavins may have been among life’s earliest organic molecules. Tungsten and molybdenum are transition metals which readily perform two-electron transitions. Molybdenum is found in cofactors of various enzymes, notably in the form of molybdo/tungsto-pterins. Enzymes containing the Mo/W-pterin cofactor perform the redox conversion of a plethora of environmental substrates with two-electron chemistry (Grimaldi et al., 2013) and several of them likely are at least as old as the LUCA (Schoepp-Cothenet et al., 2012). The Mo-pterin cofactor in specific enzymes from the superfamily can adopt strongly positive cooperative redox behavior (Hoke et al., 2004; Duval et al., 2016), that is, feature substantially inverted potentials for its individual one-electron transitions (Duval et al., 2016). The Mo/W-pterin cofactor therefore possesses all the properties necessary for electron bifurcating reactions, although a Mo-pterin-based bifurcating enzyme has not been described so far.

These observations from extant biology lead us to propose a specific succession of events occurring at the putative habitats for life’s emergence, that is, submarine alkaline hydrothermal vents, and paving the way from minerals to the biological redox cofactor flavin. Intriguingly, occasional enrichments of molybdenum are found in the 3.8 billion years old Eoarchaean Banded Iron Formation (BIF) in Isua Greenland (Appel, 1980; Dymek and Klein, 1988). This is not so surprising given its solubility in reduced alkaline hydrothermal fluids, variously as [MoO₄]^2-, [MoS₄]^2-, [FeO(OH)MoS₄]^3-, and [(Fe₂S₂)(MoS₄)₂]^4- (Helz et al., 2013; Russell and Nitschke, 2017). These ions would precipitate at the vent on meeting the acidulous and reduced early ocean, as would the double layer oxyhydroxide green rust (GR, ∼[Fe^II _6xFe^III_6(1-x) O₁₂H_2(7-3x)]²⁺ [CO₃.3.3H₂O]^2-) (Antony et al., 2008), the precursor to the iron oxide minerals presently constituting the BIF (Arrhenius, 2003; Russell et al., 2013; Tosca et al., 2016; Halevy et al., 2017; Russell, 2017). In these conditions the carbonates could have been replaced by the molybdate counter-ions proton-caged in the interlayer galleries where they may have enabled the various two-electron transfer processes now carried out by the molybdo-enzymes today (Itaya et al., 1987; Strajescu et al., 1997; Rives and Ulibarri, 1999; Evans and Slade, 2006; Wander et al., 2008; Duan et al., 2011; Duval et al., 2016).

The early history of life was marked by the subjugation, sequestering and the partial dissolution of these mineral nanoengines. Upon the emergence of the first simple organic molecules, pteridines (synthesized today in Archaea and Bacteria from guanosine triphosphate, GTP) or their derivatives, the pterins, may have started to serve as ligands for molybdenum, thereby solubilizing the transition metals out of the green rust mineral and eventually yielding the molybdopterin cofactor of extant enzymes. When the first cells ventured, or were born by currents, beyond the alkaline waters into the wider Hadean or Archaean acidulous and reduced ocean, their survival would be challenged by the dearth of soluble molybdenum anions (Anbar, 2008; Boyd et al., 2011; Zhang et al., 2014). During this period of the organic takeover, pterins may have evolved into the chemically closely related flavins (Basu and Burgmayer, 2011) which mimic the electrochemical properties of the molybdopterins based only on the aromatic ring system. Figure 7 illustrates the structural similarities and differences between pterins and flavins. The similarities mainly trace back to the pteridine moiety from which both pterins and flavins are derived while the third rings of the two systems differ in bond-saturation. This ring is fully conjugated and thus coplanar with the pteridine-part in flavins whereas its higher saturation due to the oxygen heteronucleus leads to substantial puckering highlighted by the 3D-representation at the bottom of the figure. Intriguingly, both pterins and flavins undergo two-electron oxidoreduction reactions (Figure 7). Noteworthily, even in extant molybdoenzymes, redox properties are strongly influenced by the H-bonding network surrounding the pterin moiety (Wu et al., 2015; Duval et al., 2016) just as is the case for isoalloxazine molecule in flavoenzymes. Thus, while molydopterin itself was retained for redox ranges inaccessible to the flavins, the flavoproteins themselves may have evolved to stand in for the absent molybdenum anion in many of its other former roles, especially perhaps in the redox bifurcation of electrons.

FIGURE 7

FIGURE 7. Comparison of the chemical layout (top part) and 3D-conformation (bottom part) of pterins and flavins together with a representation of their precursor-moiety, the pteridine (top). Protons indicated in italics on the flavin molecule feature pK-values in the physiological range of pH.

Author Contributions

All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.

Funding

MR’s research was supported by the National Aeronautics and Space Administration and the Centre National de la Recherche Scientifique, through the NASA Astrobiology Institute under cooperative agreement issued through the Science Mission Directorate (No. NNH13ZDA017C) (Icy Worlds) at the Jet Propulsion Laboratory.

Conflict of Interest Statement

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

Acknowledgments

The authors would like to thank Wolfgang Buckel, Rolf Thauer, and Tristan Wagner (all three MPI for Terrestrial Microbiology, Marburg, Germany) as well as Elbert Branscomb (Department of Physics, University of Illinois, Urbana, IL, United States) for stimulating discussions.

Supplementary Material

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

FIGURE S1 | Structures of flavoenzymes containing at least one substrate molecule, substrate analog, or electron transfer partner of the flavin are superimposed on their flavin ring system (in gray in the figure). Substrates are found to be mainly located on the Re-side of the flavin (above the flavin in the figure) with their ring systems oriented parallel to the isoalloxazine moiety. Cofactors serving as electron transfer partners are mainly found on either side of the isoalloxazine ring and rarely on its Re- or Si-side (i.e., above or below the plane of the isoalloxazine moiety). These spatial preferences are independent of the electron transfer properties of the flavin as evidenced by the color-coding (purple for substrates and/or electron transfer partners of gating type enzymes and orange for both-together or bifurcating types). Pdb-entry codes for the structures used in this comparison: subunit F of Complex I (pdb: 3IAM, Bos taurus), dihydroorotate dehydrogenase B, FAD domain (pdb: 1EP3, Lactococcus lactis); dihydroorotate dehydrogenase B, FMN domain (pdb: 1EP3, Lactococcus lactis), sulfite-quinone reductase (pdb: 3HYW, Aquifex aeolicus), dihydropyrimidine dehydrogenase, FMN domain (pdb: 1GTE, Sus scorfa), dihydropyrimidine dehydrogenase, FAD domain (pdb: 1GTE, Sus scorfa), adenylsulfate reductase (pdb: 1JNR, Archaeoglobus fulgidus), NADH-dependent Ferredoxin:NADP oxidoreductase, small subunit (5JCA_S, Thermotoga maritima), heterodisulfide reductase (pdb: 5ODC, Methanothermococcus), NADH-dependent glutamate synthase alpha subunit, FMN domain (pdb: 2VDC, Azospirillum brasilense), NADH-dependent glutamate synthase alpha subunit, FAD domain (pdb: 2VDC), fumarate quinol oxidoreductase (1QLB, Wolinella succinogenes), NADH-dehydrogenase 2 (pdb 4NWZ, Saccharomyces cerevisiae), sarcosine oxidase (pdb: 1X31, Corynebacterium sp. U-96), Ferredoxin NADH oxidoreductase (pdb: 1QGA, Pisum sativum), Ferredoxin: ferredoxin NADH oxidoreductase complex (1EWY, Anabaena PCC7119), Na⁺-translocating NADH oxidoreductase (pdb: 4P6V, Vibrio cholerae).

FIGURE S2 | Structures of diverse FAD- or FMN-containing enzymes. Only the flavin-binding domains are depicted. The flavin is in red. Shades of orange were used for enzymes that perform both-together electron transfer or electron bifurcation, shades of violet for enzymes that perform the gating-type electron transfer and blue for enzymes that transfer only a single electron. The enzymes are classified according to their fold. See text for a more detailed description. NDH-2, NADH-dehydrogenase 2 (pdb 4NWZ); SQR, sulfite-quinone reductase (pdb: 3HYW); Nfn-1 (L), NADH-dependent Ferredoxin:NADP oxidoreductase, large subunit (5JCA,L); Hdr, heterodisulfide reductase (pdb: 5ODC); TrxR, thioredoxin reductase (pdb: 1f6m); DPD, dihydropyrimidine dehydrogenase (pdb: 1GTE); GR-FNR, glutathione reductase-type NADH-ferredoxin oxidoreductase (pdb: 2LZX); GRT, glutathione reductase (pdb: 2GRT); NROR, NADH rubredoxin reductase (pdb: 3KLJ); AprA, adenylsulfate reductase (pdb: 1JNR); Etf:Q, electron-transfer flavoprotein-ubiquinone oxidoreductase (pdb: 2GMJ); CII (Complex II), succinate quinone oxidoreductase (1QLB); SO, sarcosine oxidase (pdb: 1X31); FNR, ferredoxin NADH oxidoreductase (pdb: 1QUE); Nqr, Na⁺-translocating NADH oxidoreductase (pdb: 4P6V); sMMO, soluble methane monooxygenase (pdb: 1TVC); Nfn-1 (S), NADH-dependent Ferredoxin:NADP oxidoreductase, small subunit (5JCA,S); b5r, cytochrome b5 reductase (pdb: 1IBO); PDR, phthalate dioxygenase reductase (pdb: 2PIA); Fld, flavodoxin (pdb: 2WC1); FIRd, flavo-diiron nitric oxide reductase (pdb: 5LMC); ETF-α, electron-transfer-flavoprotein, alpha subunit (pdb: 2L2I); ETF-β, Electron Transfer Flavoprotein, beta subunit (pdb: 2L2I); DHOdhA_FMNa, dihydroorotate dehydrogenase A, N-ter FMN domain (pdb: 1DOR); DHOdhA_FMNb, dihydroorotate dehydrogenase A, C-ter FMN domain (pdb: 1DOR); DHOdhB_FAD, dihydroorotate dehydrogenase B, FAD domain (pdb: 1EP3); DHOdhB_FMN, dihydroorotate dehydrogenase B, FMN domain (pdb: 1EP3); Fd-Glts, glutamate synthase (pdb: 1OFE); Glts α NADH-dependent glutamate synthase alpha subunit (pdb: 2VDC); Flcytb2, flavocytochrome b2 (pdb: 1FCB); IYD, iodo-tyrosine deiodinase (pdb: 4TTC); Nfr, nitro-flavin-reductase (pdb: 3N2S); Bcd, butyryl CoA dehydrogenase (pdb: 1BUC); SCAD, short chain acyl-CoA dehydrogenase (pdb: SCAD); PyrOx, pyruvate oxidase (pdb: 4FEG); Q-PyrOx, quinone pyruvate oxidase (pdb: 3EY9); NuoF, subunit F of Complex I (pdb: 3IAM).

TABLE S1 | The table represents a compilation of flavoproteins for which electrochemical properties are available in the literature and a comparison to the relevant redox potentials of free FAD and FMN as well as of the soluble 2-electron redox substrate NAD(P). The color-coding in the column labeled “reaction” is detailed at the bottom of the table and corresponds to the nomenclature for differing types of 2-electron redox reactions as used throughout the main text of this article.

Footnotes

1. ^According to the Nernst equation E_h = E_m + 2.3RT/nF^∗log([ox]/[red] with E_h, solution potential; E_m, redox midpoint potential; n, number of transferred electrons; T, temperature in Kelvin; R, universal gas constant; F, Faraday constant; note that under specific conditions the more appropriate chemical activities “a” rather than concentrations “[c]” need to be considered).
2. ^K_S = [semi-reduced]²/[ox]x[red].
3. ^^∗: for a compilation of pertinent mathematical derivations and correlations as well as graphic representations thereof, see our dedicated website at http://bip.cnrs-mrs.fr/bip09/2electron.html.

References

Anbar, A. D. (2008). Elements and evolution. Science 322, 1481–1483. doi: 10.1126/science.1163100

PubMed Abstract | CrossRef Full Text | Google Scholar

Anderson, R. F. (1980). Energetics of the one-electron steps in the NAD⁺/NADH redox couple. Biochim. Biophys. Acta 590, 277–281. doi: 10.1016/0005-2728(80)90032-8

CrossRef Full Text | Google Scholar

Antony, H., Legrand, L., and Chaussé, A. (2008). Carbonate and sulphate green rusts—mechanisms of oxidation and reduction. Electrochim. Acta 53, 7146–7156. doi: 10.1016/j.electacta.2008.05.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Appel, P. W. U. (1980). On the early Archaean Isua iron-formation, west Greenland. Precambrian Res. 11, 73–87. doi: 10.1016/0301-9268(80)90081-9

CrossRef Full Text | Google Scholar

Arrhenius, G. O. (2003). Crystals and life. Helv. Chim. Acta 86, 1569–1586. doi: 10.1002/hlca.200390135

CrossRef Full Text | Google Scholar

Basu, P., and Burgmayer, S. J. N. (2011). Pterin chemistry and its relationship to the molybdenum cofactor. Coord. Chem. Rev. 255, 1016–1038. doi: 10.1016/j.ccr.2011.02.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Batie, C. J., and Kamin, H. (1981). The relation of pH and oxidation-reduction potential to the association state of the ferredoxin ferredoxin:NADP⁺ reductase complex. J. Biol. Chem. 256, 7756–7763.

Google Scholar

Baymann, F., Lebrun, E., Brugna, M., Schoepp-Cothenet, B., Giudici-Orticoni, M.-T., and Nitschke, W. (2003). The redox protein construction kit; pre-LUCA evolution of energy conserving enzymes. Philos. Trans. R. Soc. Lond. B Biol. Sci. 358, 267–274. doi: 10.1098/rstb.2002.1184

PubMed Abstract | CrossRef Full Text | Google Scholar

Beirns, J. J., and Wang, J. H. (1972). Studies on nicotinamide adenine dinucleotide phosphate reductase of spinach chloroplasts. J. Biol. Chem. 247, 7374–7382.

Google Scholar

Bergdoll, L., Ten Brink, F., Nitschke, W., Picot, D., and Baymann, F. (2016). From low- to high-potential bioenergetic chains: thermodynamic constraints of Q-cycle function. Biochim. Biophys. Acta 1857, 1569–1579. doi: 10.1016/j.bbabio.2016.06.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Boughanemi, S., Lyonnet, J., Infossi, P., Bauzan, M., Kosta, A., Lignon, S., et al. (2016). Microbial oxidative sulfur metabolism: biochemical evidence of the membrane-bound heterodisulfide reductase-like complex of the bacterium Aquifex aeolicus. FEMS Microbiol. Lett. 363:fnw156. doi: 10.1093/femsle/fnw156

PubMed Abstract | CrossRef Full Text | Google Scholar

Boyd, E. S., Anbar, A. D., Miller, S., Hamilton, T. L., Lavin, M., and Peters, J. W. (2011). A late methanogen origin for molybdenum-dependent nitrogenase. Geobiology 9, 221–232. doi: 10.1111/j.1472-4669.2011.00278.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Branscomb, E., Biancalani, T., Goldenfeld, N., and Russell, M. J. (2017). Escapement mechanisms and the conversion of disequilibria; the engines of creation. Phys. Rep. 677, 1–60. doi: 10.1016/j.physrep.2017.02.001

CrossRef Full Text | Google Scholar

Branscomb, E., and Russell, M. J. (2013). Turnstiles and bifurcators: the disequilibrium converting engines that put metabolism on the road. Biochim. Biophys. Acta 1827, 62–78. doi: 10.1016/j.bbabio.2012.10.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Branscomb, E., and Russell, M. J. (2018). Frankenstein vs. “the bosom of the deep”: Who is responsible for abiogenesis? BioEssays 40:1700182.

Buckel, W., and Thauer, R. K. (2013). Energy conservation via electron bifurcating ferredoxin reduction and proton/Na⁺ translocating ferredoxin oxidation. Biochim. Biophys. Acta 1827, 94–113. doi: 10.1016/j.bbabio.2012.07.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Buckel, W., and Thauer, R. K. (2018a). Flavin-based electron bifurcation, a new mechanism of biological energy coupling. Chem. Rev. 118, 3862–3886. doi: 10.1021/acs.chemrev.7b00707

PubMed Abstract | CrossRef Full Text | Google Scholar

Buckel, W., and Thauer, R. K. (2018b). Flavin-based electron bifurcation, ferredoxin, flavodoxin and anaerobic respiration with protons (Ech) or NAD⁺ (Rnf) as electron acceptors: a historical review. Front. Microbiol. 9:401. doi: 10.3389/fmicb.2018.00401

PubMed Abstract | CrossRef Full Text | Google Scholar

Chowdhury, N. P., Kahnt, J., and Buckel, W. (2015). Reduction of ferredoxin or oxygen by flavin-based electron bifurcation in Megasphaera elsdenii. FEBS J. 282, 3149–3160. doi: 10.1111/febs.13308

PubMed Abstract | CrossRef Full Text | Google Scholar

Chowdhury, N. P., Klomann, K., Seubert, A., and Buckel, W. (2016). Reduction of flavodoxin by electron bifurcation and sodium ion-dependent reoxidation by NAD⁺ catalyzed by ferredoxin-NAD⁺ reductase (Rnf). J. Biol. Chem. 291, 11993–12002. doi: 10.1074/jbc.M116.726299

PubMed Abstract | CrossRef Full Text | Google Scholar

Chowdhury, N. P., Mowafy, A. M., Demmer, J. K., Upadhyay, V., Koelzer, S., Jayamani, E., et al. (2014). Studies on the mechanism of electron bifurcation catalyzed by electron transferring flavoprotein (Etf) and butyryl-CoA dehydrogenase (Bcd) of Acidaminococcus fermentans. J. Biol. Chem. 289, 5145–5157. doi: 10.1074/jbc.M113.521013

PubMed Abstract | CrossRef Full Text | Google Scholar

Cottrell, A. (1979). The natural philosophy of engines. Contemp. Phys. 20, 1–10. doi: 10.1080/00107517908227799

CrossRef Full Text | Google Scholar

Crofts, A. R., Hong, S., Wilson, C., Burton, R., Victoria, D., Harrison, C., et al. (2013). The mechanism of ubihydroquinone oxidation at the Q_o-site of the cytochrome bc₁ complex. Biochim. Biophys. Acta 1827, 1362–1377. doi: 10.1016/j.bbabio.2013.01.009

PubMed Abstract | CrossRef Full Text | Google Scholar

de Vries, S., and Grivell, L. A. (1988). Purification and characterization of a rotenone-insensitive NADH:Q6 oxidoreductase from Saccharomyces cerevisiae. Eur. J. Biochem. 176, 377–384. doi: 10.1111/j.1432-1033.1988.tb14292.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Demmer, J. K., Chowdhury, N. P., Selmer, T., Ermler, U., and Buckel, W. (2017). The semiquinone swing in the bifurcating electron transferring flavoprotein/butyryl-CoA dehydrogenase complex from Clostridium difficile. Nat. Commun. 8:1577. doi: 10.1038/s41467-017-01746-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Duan, H. D., Lubner, C. E., Tokmina-Lukaszewska, M., Gauss, G. H., Bothner, B., King, P. W., et al. (2018). Distinct properties underlie flavin-based electron bifurcation in a novel electron transfer flavoprotein FixAB from Rhodopseudomonas palustris. J. Biol. Chem. 293, 4688–4701. doi: 10.1074/jbc.RA117.000707

PubMed Abstract | CrossRef Full Text | Google Scholar

Duan, X., Lu, J., and Evans, D. G. (2011). “Assembly chemistry of anion-intercalated layered materials,” in Modern Inorganic Synthetic Chemistry, eds R. Xu and Y. Xu (ıAmsterdam: Elsevier), 375–404.

Google Scholar

Ducluzeau, A.-L., van Lis, R., Duval, S., Schoepp-Cothenet, B., Russell, M. J., and Nitschke, W. (2009). Was nitric oxide the first deep electron sink? Trends Biochem. Sci. 34, 9–15. doi: 10.1016/j.tibs.2008.10.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Duval, S., Santini, J. M., Lemaire, D., Chaspoul, F., Russell, M. J., Grimaldi, S., et al. (2016). The H-bond network surrounding the pyranopterins modulates redox cooperativity in the molybdenum-bisPGD cofactor in arsenite oxidase. Biochim. Biophys. Acta 1857, 1353–1362. doi: 10.1016/j.bbabio.2016.05.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Dymek, R. F., and Klein, C. (1988). Chemistry, petrology and origin of banded iron-formation lithologies from the 3800 Ma Isua supracrustal belt, West Greenland. Precambrian Res. 39, 247–302. doi: 10.1016/0301-9268(88)90022-8

CrossRef Full Text | Google Scholar

Evans, D. G., and Slade, R. C. (2006). “Structural aspects of layered double hydroxides,” in Layered Double Hydroxides, eds X. Duan and D. G. Evans (Berlin: Springer), 1–87.

Google Scholar

Ferguson, S. J., and Ingledew, W. J. (2008). Energetic problems faced by micro-organisms growing or surviving on parsimonious energy sources and at acidic pH: I. Acidithiobacillus ferrooxidans as a paradigm. Biochim. Biophys. Acta 1777, 1471–1479. doi: 10.1016/j.bbabio.2008.08.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Friedrich, T., and Weiss, H. (1997). Modular evolution of the respiratory NADH:ubiquinone oxidoreductase and the origin of its modules. J. Theor. Biol. 187, 529–540. doi: 10.1006/jtbi.1996.0387

PubMed Abstract | CrossRef Full Text | Google Scholar

Garcia Costas, A. M., Poudel, S., Miller, A. F., Schut, G. J., Ledbetter, R. N., Fixen, K. R., et al. (2017). Defining electron bifurcation in the electron transferring flavoprotein family. J. Bacteriol. 199:e00440-17. doi: 10.1128/JB.00440-17

PubMed Abstract | CrossRef Full Text | Google Scholar

Gnandt, E., Dörner, K., Strampraad, M. F. J., de Vries, S., and Friedrich, T. (2016). The multitude of iron–sulfur clusters in respiratory complex I. Biochim. Biophys. Acta 1857, 1068–1072. doi: 10.1016/j.bbabio.2016.02.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Grein, F., Ramos, A. R., Venceslau, S. S., and Pereira, I. A. C. (2013). Unifying concepts in anaerobic respiration: insights from dissimilatory sulfur metabolism. Biochim. Biophys. Acta 1827, 145–160. doi: 10.1016/j.bbabio.2012.09.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Grimaldi, S., Cécaldi, P., Schoepp-Cothenet, B., Guigliarelli, B., and Magalon, A. (2013). The prokaryotic Mo/W-bisPGD enzymes family: a catalytic workhorse in bioenergetics. Biochim. Biophys. Acta 1827, 1048–1085. doi: 10.1016/j.bbabio.2013.01.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Halevy, I., Alesker, M., Schuster, E. M., Popovitz-Biro, R., and Feldman, Y. (2017). A key role for green rust in the Precambrian oceans and the genesis of iron formations. Nat. Geosci. 10, 135–139. doi: 10.1038/ngeo2878

CrossRef Full Text | Google Scholar

Hanukoglu, I., and Gutfinger, T. (1989). cDNA sequence of adrenodoxin reductase. Identification of NADP-binding sites in oxidoreductases. Eur. J. Biochem. 180, 479–484. doi: 10.1111/j.1432-1033.1989.tb14671.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Hauska, G., Hurt, E., Gabellini, N., and Lockau, W. (1983). Comparative aspects of quinol-cytochrome c/plastocyanin oxidoreductases. Biochim. Biophys. Acta 726, 97–133. doi: 10.1016/0304-4173(83)90002-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Helz, G. R., Erickson, B. E., and Vorlicek, T. P. (2013). Stabilities of thiomolybdate complexes of iron: implications for retention of essential trace elements (Fe, Cu, Mo) in sulfidic waters. Metallomics 6, 1131–1140. doi: 10.1039/C3MT00217A

PubMed Abstract | CrossRef Full Text | Google Scholar

Herrmann, G., Jayamani, E., Mai, G., and Buckel, W. (2008). Energy conservation via electron-transferring flavoprotein in anaerobic bacteria. J. Bacteriol. 190, 784–791. doi: 10.1128/JB.01422-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoben, J. P., Lubner, C. E., Ratzloff, M. W., Schut, G. J., Nguyen, D. M. N., Hempel, K. W., et al. (2017). Equilibrium and ultrafast kinetic studies manipulating electron transfer: a short-lived flavin semiquinone is not sufficient for electron bifurcation. J. Biol. Chem. 292, 14039–14049. doi: 10.1074/jbc.M117.794214

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoke, K. R., Cobb, N., Armstrong, F. A., and Hille, R. (2004). Electrochemical studies of arsenite oxidase: an unusual example of a highly cooperative two-electron molybdenum center. Biochemistry 43, 1667–1674. doi: 10.1021/bi0357154

PubMed Abstract | CrossRef Full Text | Google Scholar

Ingledew, W. J. (1982). Thiobacillus ferrooxidans. The bioenergetics of an acidophilic chemolithotroph. Biochim. Biophys. Acta 683, 89–117. doi: 10.1016/0304-4173(82)90007-6

CrossRef Full Text | Google Scholar

Itaya, K., Chang, H. C., and Uchida, I. (1987). Anion-exchanged clay (hydrotalcite-like compounds) modified electrodes. Inorg. Chem. 26, 624–626. doi: 10.1021/ic00251a028

CrossRef Full Text | Google Scholar

Jacob, F. (1977). Evolution and tinkering. Science 196, 1161–1166. doi: 10.1126/science.860134

CrossRef Full Text | Google Scholar

Kao, W. C., and Hunte, C. (2014). The molecular evolution of the Qo motif. Genome Biol. Evol. 6, 1894–1910. doi: 10.1093/gbe/evu147

CrossRef Full Text | Google Scholar

Lebrun, E., Santini, J. M., Brugna, M., Ducluzeau, A.-L., Ouchane, S., Schoepp-Cothenet, B., et al. (2006). The Rieske protein; a case study on the pitfalls of multiple sequence alignments and phylogenetic reconstruction. Mol. Biol. Evol. 23, 1180–1191. doi: 10.1093/molbev/msk010

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, F., Hinderberger, J., Seedorf, H., Zhang, J., Buckel, W., and Thauer, R. K. (2008). Coupled ferredoxin and crotonyl coenzyme A (CoA) reduction with NADH catalyzed by the butyryl-CoA dehydrogenase/Etf complex from Clostridium kluyveri. J. Bacteriol. 190, 843–850. doi: 10.1128/JB.01417-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Lim, Z. H., Lay Khim Chng, E., Hui, Y., and Webster, R. D. (2013). The hydrogen-bonded dianion of Vitamin K1 produced in aqueous-organic solutions exists in equilibrium with its hydrogen-bonded semiquinone anion radical. J. Phys. Chem. B 117, 2396–2402. doi: 10.1021/jp4003966

PubMed Abstract | CrossRef Full Text | Google Scholar

Lubner, C. E., Jennings, D. P., Mulder, D. W., Schut, G. J., Zadvornyy, O. A., Hoben, J. P., et al. (2017). Mechanistic insights into energy conservation by flavin-based electron bifurcation. Nat. Chem. Biol. 13, 655–659. doi: 10.1038/nchembio.2348

PubMed Abstract | CrossRef Full Text | Google Scholar

Magalon, A., Ceccaldi, P., and Schoepp-Cothenet, B. (2016). “The prokaryotic Mo/W-bisPGD enzymes family,” in Molybdenum and Tungsten Enzymes, ed. R. Hille (Cambridge: RSC Publishing).

Google Scholar

Marreiros, B. C., Sena, F. V., Sousa, F. M., Batista, A. P., and Pereira, M. M. (2016). Type II NADH:quinone oxidoreductase family: phylogenetic distribution, structural diversity and evolutionary divergences. Environ. Microbiol. 18, 4697–4709. doi: 10.1111/1462-2920.13352

PubMed Abstract | CrossRef Full Text | Google Scholar

Michaelis, L. (1932). Theory of the reversible two-step oxidation. J. Biol. Chem. 96, 703–715. doi: 10.1021/ic201601q

PubMed Abstract | CrossRef Full Text | Google Scholar

Mitchell, P. (1975). The protonmotive Q cycle: a general formulation. FEBS Lett. 59, 137–139. doi: 10.1016/0014-5793(75)80359-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Mitchell, P. (1976). Possible molecular mechanisms of protonmotive function of cytochrome systems. J. Theor. Biol. 62, 327–367. doi: 10.1016/0022-5193(76)90124-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Moser, C. C., Keske, J. M., Warncke, K., Farid, R. S., and Dutton, P. L. (1992). Nature of biological electron transfer. Nature 355, 796–802. doi: 10.1038/355796a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Nitschke, W., and Bonnefoy, V. (2016). “Chapter 2 Energy acquisition in low-pH environments,” in Acidophiles, eds R. Quatrini and D. B. Johnson (Poole: Caister Academic Press).

Nitschke, W., McGlynn, S. E., Milner-White, E. J., and Russell, M. J. (2013). On the antiquity of metalloenzymes and their substrates in bioenergetics. Biochim. Biophys. Acta 1827, 871–881. doi: 10.1016/j.bbabio.2013.02.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Nitschke, W., and Russell, M. J. (2012). Redox bifurcations: mechanisms and importance to life now, and at its origin: a widespread means of energy conversion in biology unfolds. Bioessays 34, 106–109. doi: 10.1002/bies.201100134

PubMed Abstract | CrossRef Full Text | Google Scholar

Ojha, S., Meng, E. C., and Babbitt, P. C. (2007). Evolution of function in the “two dinucleotide binding domains” flavoproteins. PLoS Comput. Biol. 3:e121. doi: 10.1371/journal.pcbi.0030121

PubMed Abstract | CrossRef Full Text | Google Scholar

Osyczka, A., Moser, C. C., Daldal, F., and Dutton, P. L. (2004). Reversible redox energy coupling in electron transfer chains. Nature 427, 607–612. doi: 10.1038/nature02242

PubMed Abstract | CrossRef Full Text | Google Scholar

Peters, J. W., Beratan, D. N., Schut, G. J., and Adams, M. W. (2018). On the nature of organic and inorganic centers that bifurcate electrons, coupling exergonic and endergonic oxidation-reduction reactions. Chem. Commun. 54, 4091–4099. doi: 10.1039/c8cc01530a

PubMed Abstract | CrossRef Full Text | Google Scholar

Peters, J. W., Miller, A. F., Jones, A. K., King, P. W., and Adams, M. W. (2016). Electron bifurcation. Curr. Opin. Chem. Biol. 31, 146–152. doi: 10.1016/j.cbpa.2016.03.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Poudel, S., Tokmina-Lukaszewska, M., Colman, D. R., Refai, M., Schut, G. J., King, P. W., et al. (2016). Unification of [FeFe]-hydrogenases into three structural and functional groups. Biochim. Biophys. Acta 1860, 1910–1921. doi: 10.1016/j.bbagen.2016.05.034

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramos, A. R., Grein, F., Oliveira, G. P., Venceslau, S. S., Keller, K. L., Wall, J. D., et al. (2015). The FlxABCD-HdrABC proteins correspond to a novel NADH dehydrogenase/heterodisulfide reductase widespread in anaerobic bacteria and involved in ethanol metabolism in Desulfovibrio vulgaris Hildenborough. Environ. Microbiol. 17, 2288–2305. doi: 10.1111/1462-2920.12689

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramos, A. R., Keller, K. L., Wall, J. D., and Pereira, I. A. (2012). The membrane QmoABC complex tnteracts directly with the dissimilatory adenosine 5′-phosphosulfate reductase in sulfate reducing bacteria. Front. Microbiol. 3:137. doi: 10.3389/fmicb.2012.00137

CrossRef Full Text | Google Scholar

Rao, S. T., and Rossman, M. G. (1973). Comparisons of super-secondary structures in proteins. J. Mol. Biol. 76, 241–256. doi: 10.1016/0022-2836(73)90388-4

CrossRef Full Text | Google Scholar

Rives, V., and Ulibarri, M. A. (1999). Layered double hydroxides (LDH) intercalated with metal coordination compounds and oxometalates. Coord. Chem. Rev. 181, 61–120. doi: 10.1016/S0010-8545(98)00216-1

CrossRef Full Text | Google Scholar

Russell, M. J. (2017). Life is a verb, not a noun. Geology 45, 1143–1144. doi: 10.1130/focus112017.1

CrossRef Full Text | Google Scholar

Russell, M. J., and Nitschke, W. (2017). Methane: fuel or exhaust at the emergence of life. Astrobiology 17, 1053–1066. doi: 10.1089/ast.2016.1599

PubMed Abstract | CrossRef Full Text | Google Scholar

Russell, M. J., Nitschke, W., and Branscomb, E. (2013). The inevitable journey to being. Philos. Trans. R. Soc. Lond. B Biol. Sci. 368:20120254. doi: 10.1098/rstb.2012.0254

PubMed Abstract | CrossRef Full Text | Google Scholar

Sánchez-Azqueta, A., Herguedas, B., Hurtado-Guerrero, R., Hervás, M., Navarro, J. A., Martínez-Júlvez, M., et al. (2014). A hydrogen bond network in the active site of Anabaena ferredoxin-NADP⁺ reductase modulates its catalytic efficiency. Biochim. Biophys. Acta 1837, 251–263. doi: 10.1016/j.bbabio.2013.10.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Sato, K., Nishina, Y., and Shiga, K. (2013). Interaction between NADH and electron-transferring flavoprotein from Megasphaera elsdenii. J. Biochem. 153, 565–572. doi: 10.1093/jb/mvt026

PubMed Abstract | CrossRef Full Text | Google Scholar

Sazanov, L. A. (2016). Structure of bacterial respiratory complex I. Biochim. Biophys. Acta 1857, 892–901. doi: 10.1016/j.bbabio.2016.01.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Schoepp-Cothenet, B., van Lis, R., Atteia, A., Baymann, F., Capowiez, L., Ducluzeau, A.-L., et al. (2013). On the universal core of bioenergetics. Biochim. Biophys. Acta 1827, 79–93. doi: 10.1016/j.bbabio.2012.09.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Schoepp-Cothenet, B., van Lis, R., Philippot, P., Magalon, A., Russell, M. J., and Nitschke, W. (2012). The ineluctable requirement for the trans-iron elements molybdenum and/or tungsten in the origin of life. Sci. Rep. 2:263. doi: 10.1038/srep00263

PubMed Abstract | CrossRef Full Text | Google Scholar

Schrödinger, E. (1944). What Is Life? The Physical Aspect of the Living Cell and Mind. Cambridge: Cambridge University Press.

Google Scholar

Schuchmann, K., and Müller, V. (2012). A bacterial electron-bifurcating hydrogenase. J. Biol. Chem. 287, 31165–31171. doi: 10.1074/jbc.M112.395038

PubMed Abstract | CrossRef Full Text | Google Scholar

Schut, G. J., and Adams, M. W. W. (2009). The Iron-hydrogenase of Thermotoga maritima utilizes ferredoxin and NADH synergistically: a new perspective on anaerobic hydrogen production. J. Bacteriol. 191, 4451–4457. doi: 10.1128/JB.01582-08

PubMed Abstract | CrossRef Full Text | Google Scholar

Schütz, M., Brugna, M., Lebrun, E., Baymann, F., Huber, R., Stetter, K.-O., et al. (2000). Early evolution of cytochrome bc complexes. J. Mol. Biol. 300, 663–676. doi: 10.1006/jmbi.2000.3915

PubMed Abstract | CrossRef Full Text | Google Scholar

Sétif, P. (2001). Ferredoxin and flavodoxin reduction by photosystem I. Biochim. Biophys. Acta 1507, 161–179. doi: 10.1016/S0005-2728(01)00205-5

CrossRef Full Text | Google Scholar

Sled, V. D., Rudnitzky, N. I., Hatefi, Y., and Ohnishi, T. (1994). Thermodynamic analysis of flavin in mitochondrial NADH:ubiquinone oxidoreductase (Complex I). Biochemistry 33, 10069–10075. doi: 10.1021/bi00199a034

CrossRef Full Text | Google Scholar

Song, Y., and Buettner, G. R. (2010). Thermodynamic and kinetic considerations for the reaction of semiquinone radicals to form superoxide and hydrogen peroxide. Free Radic. Biol. Med. 49, 919–962. doi: 10.1016/j.freeradbiomed.2010.05.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Strajescu, M., Gabrus, R., Bratean, C., and Jeleriu, S. (1997). Dinitrogen reduction in hydroxide systems of Mo3+, Fe2+, Ti3+, and Mg2+. Rev. Roum. Chim. 42, 435–438.

Google Scholar

Tan, S. L. J., Novianti, M. L., and Webster, R. D. (2015). Effects of low to intermediate water concentrations on proton-coupled electron transfer (PCET) reactions of flavins in aprotic solvents and a comparison with the PCET reactions of quinones. J. Phys. Chem. B 119, 14053–14064. doi: 10.1021/acs.jpcb.5b07534

PubMed Abstract | CrossRef Full Text | Google Scholar

Ten Brink, F., Schoepp-Cothenet, B., van Lis, R., Nitschke, W., and Baymann, F. (2013). Multiple Rieske/cytb complexes in a single organism. Biochim. Biophys. Acta 1827, 1392–1406. doi: 10.1016/j.bbabio.2013.03.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Tosca, N. J., Guggenheim, S., and Pufahl, P. K. (2016). An authigenic origin for Precambrian greenalite: implications for iron formation and the chemistry of ancient seawater. Geol. Soc. Am. Bull. 128, 511–530. doi: 10.1130/B31339.1

CrossRef Full Text | Google Scholar

Verhagen, M. F., O’Rourke, T., and Adams, M. W. (1999). The hyperthermophilic bacterium, Thermotoga maritima, contains an unusually complex iron-hydrogenase: amino acid sequence analyses versus biochemical characterization. Biochim. Biophys. Acta 1412, 212–229. doi: 10.1016/S0005-2728(99)00062-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Wagner, T., Koch, J., Ermler, U., and Shima, S. (2017). Methanogenic heterodisulfide reductase (HdrABC-MvhAGD) uses two noncubane [4Fe-4S] clusters for reduction. Science 357, 699–703. doi: 10.1126/science.aan0425

PubMed Abstract | CrossRef Full Text | Google Scholar

Wander, M. C., Kubicki, J. D., and Schoonen, M. A. (2008). Reduction of N2 by Fe²⁺ via homogeneous and heterogeneous reactions part 2: the role of metal binding in activating N₂ for reduction; a requirement for both pre-biotic and biological mechanisms. Orig. Life Evol. Biospheres 38, 195–209. doi: 10.1007/s11084-008-9133-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, S., Huang, H., Kahnt, J., Mueller, A. P., Köpke, M., and Thauer, R. K. (2013a). NADP-specific electron-bifurcating [FeFe]-hydrogenase in a functional complex with formate dehydrogenase in Clostridium autoethanogenum grown on CO. J. Bacteriol. 195, 4373–4386. doi: 10.1128/JB.00678-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, S., Huang, H., Kahnt, J., and Thauer, R. K. (2013b). A reversible electron-bifurcating ferredoxin- and NAD-dependent [FeFe]-hydrogenase (HydABC) in Moorella thermoacetica. J. Bacteriol. 195, 1267–1275. doi: 10.1128/JB.02158-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Wikström, M. K., and Berden, J. A. (1972). Oxidoreduction of cytochrome b in the presence of antimycin. Biochim. Biophys. Acta 283, 403–420. doi: 10.1016/0005-2728(72)90258-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, S. Y., Rothery, R. A., and Weiner, J. H. (2015). Pyranopterin coordination controls molybdenum electrochemistry in Escherichia coli nitrate reductase. J. Biol. Chem. 290, 25164–25173. doi: 10.1074/jbc.M115.665422

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, H., Osyczka, A., Dutton, P. L., and Moser, C. C. (2007). Exposing the Complex III Qo semiquinone radical. Biochim. Biophys. Acta 1767, 883–887. doi: 10.1016/j.bbabio.2007.04.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, P., Yuly, J. L., Lubner, C. E., Mulder, D. W., King, P. W., Peters, J. W., et al. (2017). Electron bifurcation: thermodynamics and kinetics of two-electron brokering in biological redox chemistry. Acc. Chem. Res. 50, 2410–2417. doi: 10.1021/acs.accounts.7b00327

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, X., Sigman, D. M., Morel, F. M., and Kraepiel, A. M. (2014). Nitrogen isotope fractionation by alternative nitrogenases and past ocean anoxia. Proc. Natl. Acad. Sci. U.S.A. 111, 4782–4787. doi: 10.1073/pnas.1402976111

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., Huang, L., Shulmeister, V. M., Chi, Y., Kim, K. K., Hung, L. W., et al. (1998). Electron transfer by domain movement in cytochrome bc1. Nature 392, 677–684. doi: 10.1038/33612

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, Y., Kahnt, J., Kwon, I. H., Mackie, R. I., and Thauer, R. K. (2014). Hydrogen formation and its regulation in Ruminococcus albus: involvement of an electron-bifurcating [FeFe]-hydrogenase, of a non-electron-bifurcating [FeFe]-hydrogenase, and of a putative hydrogen-sensing [FeFe]-hydrogenase. J. Bacteriol. 196, 3840–3852. doi: 10.1128/JB.02070-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Zuckerkandl, E., and Pauling, L. (1965). Molecules as documents of evolutionary history. J. Theor. Biol. 8, 357–366. doi: 10.1016/0022-5193(65)90083-4

PubMed Abstract | CrossRef Full Text | Google Scholar