Photoinduced Electron Transfer Informs on Pathway Coupling in Flavin-Based Electron Bifurcation
Seth A. Wiley, Carolyn E. Lubner

TL;DR
This paper investigates how electron transfer is controlled in a key enzyme from a heat-loving microbe, revealing insights into how high-energy electrons are managed in a process called electron bifurcation.
Contribution
The study introduces a novel method combining low-temperature photoexcitation and EPR to probe unstable intermediates in flavin-based electron bifurcation.
Findings
A radical intermediate and its nearby [4Fe-4S] cluster were observed to grow together during cryogenic illumination.
Photogenerated species remained stable at low temperatures and recombined at higher temperatures.
The study suggests a potential gating mechanism involving residue movement for reversible electron flow.
Abstract
Flavin-based electron bifurcation (FBEB) is an enzymatic mechanism that generates extremely high-energy electrons to drive unfavorable chemical reactions. It is utilized by the NADH-dependent ferredoxin:NADP+-oxidoreductase (Nfn) enzyme in hyperthermophile Pyrococcus furiosus to bifurcate electrons from NADPH into the coupled low-potential (endergonic) and high-potential (exergonic) pathways. This process enables P. furiosus to live in harsh and uninviting environments. Despite its biological importance, the mechanisms used by Nfn to facilitate exceptional directional control over short-lived, high-energy electrons and to prevent undesired transfer, particularly along the low-potential pathway, are still not well understood. To elucidate how the protein environment contributes to electronic control in the low-potential pathway, new techniques must be utilized to probe these unstable…
Genes, proteins, chemicals, diseases, species, mutations and cell lines named across the full text — each resolved to its canonical identifier and authoritative record.
Click any figure to enlarge with its caption.
1
2
3
4
5
6
7- —Basic Energy Sciences10.13039/100006151
- —Chemical Sciences, Geosciences, and Biosciences Division10.13039/100013145
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsLight effects on plants · Photosynthetic Processes and Mechanisms · Photochemistry and Electron Transfer Studies
Introduction
The ability to control the direction of energy flow is a crucial part of how organisms maintain homeostasis and adapt to dynamic or extreme environments, such as hydrothermal vents. ?,? Critically, to survive in such conditions, organisms must control the energetic flux to gain every advantage they can while preventing unnecessary energy loss. Over the past 50+ years ?,? there has been much interest in determining how these energy-efficient biochemical processes occur, and more recently, a lot of attention has been focused on the phenomenon of electron bifurcation. ?,?
Electron bifurcation is a mechanism of energy conservation many organisms use to minimize energy loss during interconversion of cellular currencies for metabolic processes.? Bifurcation promotes the reversible conversion of a middle-energy starting material/substrate to a coupled pair of high- and low-energy reducing equivalents that the cell can use for more energy-demanding processes.? In this mechanism, the overall thermodynamic landscape dictates how the energetic sacrifice in one electron is coupled to the increase of energy in the other, resulting in an energetic redistribution between the electrons. ?,? At the thermodynamic limits of life, many microbes exploit electron bifurcation to redistribute energetic equivalents and increase the overall efficiency of their metabolic processes.? Flavin-Based Electron Bifurcation (FBEB) uses a flavin cofactor as the site of electron bifurcation? which are apt for facilitating both 1- and 2-electron transfers by way of their 3 predominant redox states: oxidized quinone (OX), 1-electron reduced semiquinone (SQ), and 2-electron reduced hydroquinone (HQ). ?,?
In the hyperthermophillic organism Pyrococcus furiosus, FBEB is performed by the NADH-dependent ferredoxin:NADP^+^-oxidoreductase complex I, also known as NfnSL (Figure).? This dimeric enzyme contains two [4Fe-4S] clusters in the large subunit, one [2Fe-2S] cluster in the small subunit, and two flavin adenine dinucleotides (FAD), one in each subunit (FigureA), assigned respectively as the L-FAD in NfnL and the S-FAD in NfnS. Overall, the NfnSL complex carries out NADPH oxidation coupled to the simultaneous reduction of NAD^+^ and ferredoxin (Fd) (eq; FigureC).
NfnSL from Pyrococcus furiosus. Panel A shows the NfnSL complex containing NfnS (green) and NfnL (blue), with NfnS containing the [2Fe-2S] cluster and S-FAD, and NfnL containing the 2x [4Fe-4S] clusters and L-FAD, respectively. The low-potential pathway is located in the NfnL subunit, with the pathway and key residues, R201 and R333, and key distances highlighted in panel B. The 5.4 Å distance was measured between the O4' to C101 sulfur (7.4 Å O4' to Fe). Panel C depicts the overall reaction catalyzed by the NfnSL complex, with the L-FAD as the site of electron bifurcation (center), the high-potential pathway on the bottom in NfnS (blue arrows), and the low-potential pathway on top in NfnL (rust arrows). Once the L-FAD is reduced, one electron transfers “downhill” to the [2Fe-2S] cluster in NfnS, and then the other electron transfers from the L-FAD to the nearby proximal [4Fe-4S] cluster, then to the distal [4Fe-4S] cluster, and ultimately to a partner ferredoxin. NfnL is outlined in panel A as it is the only subunit used in the experiments in this report. PDB:5JFC.
The subunits of the NfnSL complex are effectively separated into the high-potential pathway in NfnS and the low-potential pathway in NfnL, with the bifurcating flavin, L-FAD, located in NfnL at the fork between the two pathways. The bifurcation reaction is triggered by a spontaneous downhill electron transfer (ET) through the high-potential pathway, leaving a transient, unstable anionic semiquinone (ASQ) radical on the L-FAD, which then rapidly undergoes ET via the low-potential pathway to reduce externally located Fd. ?,? Recently, much has been learned about the thermodynamic landscape in NfnSL, particularly regarding how the cofactors are poised for electron transfer.? The midpoint redox potentials of the high- and low-potential pathways in NfnSL have been investigated by square-wave voltammetry, showing that within the low-potential pathway the iron–sulfur cluster closest to the L-FAD, the differentially ligated proximal [4Fe-4S] cluster, has a midpoint potential of ∼ −700 mV, and that the cluster further away, the distal [4Fe-4S] cluster, has a midpoint potential of ∼ −520 mV.? The L-FAD has been found to have crossed-potentials, where under equilibrium conditions only the oxidized (“L-FAD_OX_”) and hydroquinone states (“L-FAD_HQ_”), but not the semiquinone state (“L-FAD_SQ_”), are observed, with a 2 e^–^ midpoint redox potential of −435 mV.? Despite the L-FAD showing only 2 electron chemistry under equilibrium conditions, ultrafast transient absorption spectroscopy (TAS) has shown the existence of a highly unstable ASQ with a lifetime of only 10 ps, that was coincident with an L-FAD_OX_ bleach. ?,?,? Based on the 7.5 Å edge-to-edge distance between the L-FAD and the proximal [4Fe-4S] cluster, the proximal [4Fe-4S] cluster’s midpoint potential, and the lifetime (10 ps) of the ASQ, the L-FAD OX to ASQ midpoint potential of −911 mV was calculated using the Moser-Dutton equation. ?,? These recent findings have revealed the low-potential pathway in NfnL is poised for rapid, downstream electron transfer once the electron pair is separated and the ASQ is generated.?
This rapid electron transfer in the low potential pathway is further controlled by key residues located within the cofactor microenvironments, such as the site-differentiation of the proximal [4Fe-4S] cluster? and the residues surrounding the L-FAD.? The role of site-differentiation in pathway coupling has been investigated, with key studies showing how the high-potential pathway? and the site-differentiated ligand on the proximal [4Fe-4S] cluster ?,? dictate the control of electron transfer in NfnSL. Without the high-potential electron acceptor, NAD^+^, NfnSL is capable of only 3 bifurcation events before stalling due to the full occupancy of the high-potential pathway cofactors (1x [2Fe-2S] and 2x S-FAD).? When the proximal [4Fe-4S] cluster’s unique glutamic acid ligation is substituted to the canonical cysteine, the low- and high-potential pathways appear to decouple, leading to bifurcation activity without the high-potential electron acceptor.? How NfnSL tunes the environments surrounding the redox cofactors to control electron transfer, in particular, the initial high-energy steps of electron bifurcation, is critical for exploiting and expanding this unique chemical transformation mechanism.
Two notable arginine residues, R201 and R333 (FigureB), either directly interface with L-FAD N5 and O4’ (R201)? or are highly conserved only in electron bifurcating enzymes (R333)? and appear to be critical for facilitating, tuning, and stabilizing the transient ASQ intermediate on the L-FAD. How nearby arginines and the [4Fe-4S] site-differentiated glutamate residue influence electronic control and directional bias are of significant interest, but to gain insight into each residue’s contribution to the overall mechanism, methods to interrogate the high-energy electron transfer steps must be realized.
Ultrafast tools, like femtosecond laser spectroscopy, have proved indispensable for probing transient intermediates and enzyme dynamics of light-activated enzymes with high temporal resolution.? In addition, these methods can elucidate the sequence of events following photon absorption, providing the critical temporal context to the steps underlying enzymatic mechanisms postactivation.? While ultrafast TAS experimental evidence allowed for determination of the mechanistically relevant L-FAD intermediate, the short lifetime has complicated the assessment and investigation of the impact of the surrounding redox environment on enabling this ASQ for subsequent electron transfer down the low-potential pathway. How these nearby residues and cofactors within the low-potential pathway interact with the ASQ intermediate and prevent short-circuit electron transfer to the high-potential pathway is currently unknown. To better understand how energy is coupled in electron bifurcating systems, it is necessary to study the high energy electron transfer step from the ASQ to the proximal [4Fe-4S] cluster. As photoexcitation has been previously used to study many photosensitive biomolecules and proteins, ?−? ? we aimed to adapt cryo-photoexcitation methods ?−? ? to an FBEB enzyme and utilize the intrinsic photophysical properties of the L-FAD and cryogenic temperatures to trap and study the mechanistically relevant intermediates by electron paramagnetic resonance (EPR) spectroscopy. EPR spectroscopy allows for monitoring of characteristic spectral features (known as g-values) for a variety of paramagnetic species, such as the [4Fe-4S] clusters and flavin radical species. Tracking these paramagnetic species sheds new light on the mechanism of bifurcation, the low-potential pathway’s electron transfer network, and allows for determining fundamental thermodynamic, physical, and electronic characteristics of key bifurcation intermediates.
In this report, we have observed reversible photoinduced electron transfer, showing the unambiguous generation of both the ASQ and proximal [4Fe-4S] cluster signals in the isolated NfnL subunit in the presence of NADPH. The signals for these photogenerated species increased in intensity with illumination time, and when stored at liquid nitrogen temperatures (∼77 K), the ASQ and proximal cluster species remained stable without losing intensity. Annealing led to pronounced recombination within the low-potential pathway, signifying a physical gating step in the directionality of electron transfer. The outcomes gained from this work expand the model of pathway coupling and electron transfer control in FBEB.
Results
To investigate the low-potential pathway and unambiguously probe the bifurcating L-FAD without interference from the additional S-FAD in NfnS, we performed studies utilizing only the NfnL subunit. We have previously shown that this methodology results in full incorporation of the NfnL cofactors, and the binding with NADP(H), L-FAD intermediate identity (ASQ) and lifetime, and cofactor reduction potentials are unaltered compared to NfnSL. ?,?
L-FAD ASQ and Proximal [4Fe-4S] Cluster Signals Form upon Photoexcitation
Previous experiments detected reduction of only the distal [4Fe-4S] cluster in the presence of NADPH.? To probe the transient ASQ observed via ultrafast TAS experiments, ?,? we photoexcited the oxidized L-FAD present in NfnL in the presence of its native substrate, NADPH. These experiments were conducted at 233 K to separate electron transfer processes from structural changes to have the best chance of trapping high energy intermediates.? As depicted in Figure, upon photoexcitation of NfnL in the presence of 8.8 mM NADPH at 233 K with 2 W 405 nm illumination for 1 h, two new species arose on top of the existing rhombic distal [4Fe-4S] cluster signal with g-values of 2.034, 1.934, and 1.896 (Figures S1 and S2, Table S2): one isotropic, organic radical-like spectrum with a g iso = 2.0033, and another resembling the previously observed proximal [4Fe-4S] cluster with newly reassigned semiaxial g-values corresponding to 1.972, 1.954, and 1.924 (Figures S3–S6, Table S3). These two emergent signals were distinct, and each of the two species had unique relaxation characteristics consistent with an organic radical, assigned to the L-FAD, and the NfnL proximal [4Fe-4S] cluster. The proximal [4Fe-4S] cluster signal was consistent with a fast-relaxing iron–sulfur cluster, losing signal above 30 K, in contrast with the underlying distal cluster signal, which relaxes more slowly and is observable up to 60 K. A reference sample consisting of 100 μM NfnL with 1 mM dithionite (“DTH”) was measured and further confirmed the proximal [4Fe-4S] cluster signal identity (10x DTH; FigureB). At a temperature of 100 K, the radical signal was stable and remained visible even up to 230 K, where the higher temperature allowed the sample to anneal halfway back to starting conditions over the course of 45 min, nearly returning to the preillumination spectrum after ∼ 5 min at 250 K (t 1/2 ≥ 30 min; Figure S7). Not only was this radical longer lived than anticipated, but it also had well-behaved signal decay kinetics over the course of 5 h at 230 K, indicative of a complex mechanism of electron recombination within the low potential pathway, particularly between the radical and proximal [4Fe-4S] cluster (Figure S8, Table S4).
233 K illumination of NfnL in the presence of 8.8 mM NADPH with 2 W 405 nm for 1 h reveals two new species not seen in the preillumination spectrum. Power and modulation amplitude used for these data were set to 10 mW and 10 G, respectively. Panel A shows the characteristic distal [4Fe-4S] cluster before illumination (Dark; lavender) with g-values of 2.035, 1.934, and 1.896. Illumination with 405 nm light reveals two new species with distinct temperature dependencies (Light; navy, blue, light-blue) corresponding to the proximal [4Fe-4S] cluster (g-values of 1.972, 1.954, and 1.924) only observable at 20 K and below, and an ASQ radical (g iso = 2.0033) observable up to 230 K (Figures S4–S7). Upon annealing the sample’s temperature to 250 K for ≤ 5 min (Figure S7), the sample was able to nearly return to the initial Dark condition (Anneal; red). Panel B shows the 20 K photoexcitation data from panel A overlaid with dithionite (DTH) references to confirm the identities of paramagnetic species. The normalized 1x DTH (gray) data overlaps nearly identically to the Dark and Anneal spectra when corrected for intensity differences. With the exception of an organic radical at g = 2.0033, the 10x DTH (cyan) spectrum is consistent with the proximal [4Fe-4S] cluster g-values observed in the illuminated samples, confirming the assignment of the 20 K Light (blue) spectra.
Illuminating oxidized NfnL alone without any NADPH generated only 1% of the total spectral intensity observed in Figure. Further, NfnL, FAD, or NADPH on their own was unable to reproduce the photogenerated species (Figures S9–S11). These control experiments suggested that the photogenerated species originate from the interaction between NADPH and the L-FAD in NfnL and require an oxidized L-FAD to absorb an ∼ 405 nm photon (Figure S12). Additionally, after 1 h of illumination at 233 K with 405 nm, no apparent photodegradation was observed for NfnL alone, implying the photogenerated species are part of a nondestructive and catalytically relevant pathway in NfnL (Figures S13–S15).
We assigned the observed photogenerated radical identity as an ASQ, based on its line width, ?,?−? ? line shape, ?−? ?,? spin–spin interaction with the proximal [4Fe-4S] cluster (Figures, S6), and the consistency with previous TAS data. ?,? We compared the photogenerated radical’s peak-to-peak line width to the established neutral semiquinone (NSQ) signal observed on the S-FAD located in the full NfnSL complex (FigureA,C).? As shown in Figure, the photoexcitated sample generated a flavin radical signal with a peak-to-peak line width of 13.4 G at its narrowest (modulation amplitude at 1 G), in contrast to the narrowest S-FAD NSQ line width of 25.0 G, slightly broader than typical NSQ line widths of 18 to 20 G ?,? due to spin–spin coupling between the [2Fe-2S] and the NSQ. ?,?,? The large difference in peak-to-peak line width between the two radicals confirms the photogenerated radical is distinct from the S-FAD NSQ, and the 13.4 G line width is consistent with previous ASQ assignments. ?,?,?,? The ASQ radical also exhibited minor, but consistent anisotropy, aka “wings”, flanking the main isotropic signal seen in Figure. The “wings” were consistent in all tested spectral parameters (temperature, power, modulation amplitude), suggesting the radical “wings” are a feature of a single radical species and not from two different overlapping radical species. ?−? ?
FigureB depicts the normalized, overlaid spectra of the photogenerated radical taken at the listed temperatures, showing only one species is apparent up to a temperature of 200 K.
*Photogenerated radical line width is consistent with an ASQ and inconsistent with an NSQ. Direct comparison of peak-to-peak line widths for the photogenerated flavin radical (blue traces) to the NSQ in the S-FAD of NfnSL (charcoal traces). Panel A shows the full radical spectra, with reference lines showing the peak measurements. Reducing the modulation amplitude from 10 to 1 G did not narrow the peak-to-peak line width for the NSQ, but did however narrow the photogenerated radical line width by ∼ 2 G. Panel B depicts normalized spectra of the observed radical species with a line width of 13.4 G at a modulation amplitude of 1 G over multiple temperatures (80 to 200 K) and shows the radical “wings” are part of the main radical feature and not from an underlying unique radical species. −
Panel B swept 100 G on either side of the 3350 G center field, from 3250 to 3450 G, resulting in high resolution of the radical species. The microwave power was 10 mW for panels A and B, the temperature was 100 K for panel A, and modulation amplitude is explicit for each trace, for which the intensities have been corrected for. Illumination at 405 nm was 90 min for the photoexcited sample spectra.*
To assess the order in which the paramagnetic signals appeared in response to illumination, a time course was conducted, probing both g ∼ 2 and g ∼ 4 regions at temperatures of 3.6, 20, and 100 K. As seen in Figure, with increasing illumination time, the low potential pathway becomes populated with singly reduced paramagnetic species. Within the first 10 min of illumination, the radical is reduced immediately, preceding the reduction of both the proximal [4Fe-4S] cluster and distal [4Fe-4S] cluster (Figure S16). After 10 min, the constant illumination conditions led to complex rates of formation for the paramagnetic species, although these species all increased in intensity as the illumination time increased (Figure S16). After 90 min of illumination, the total signal intensity corresponding to the distal [4Fe-4S] cluster appeared to be roughly double what it was at the beginning, implying that a noteworthy portion of the preilluminated NfnL sample had no reduced distal [4Fe-4S] cluster (FigureA). This increase in distal [4Fe-4S] cluster intensity underlies the emergence of the radical and reduced proximal [4Fe-4S] cluster and increased with longer illumination times.
Illumination of NfnL with NADPH shows growth of both radical and proximal [4Fe-4S] cluster signals over the course of 90 min. Panel A depicts the overlaid spectra corresponding to the increase in the radical signal at 100 K (above) and the proximal [4Fe-4S] cluster signal at 20 K (below). In addition to the ASQ and proximal [4Fe-4S] clusters, there is an increase in the overall distal [4Fe-4S] cluster signal intensity as well. Panel B shows the low-field region at 3.6 K, where a signal at g ∼ 3.89 arises with increased illumination time (30+ min), due to the interaction between the reduced distal and proximal [4Fe-4S] clusters. The aberrant signal at g ∼ 3.7 for the 90 min low-field spectrum in B is a spectral artifact of unknown origin and likely unrelated to the experimental results (Figure S17). Spectra parameters: A = 10 mW power, 10 G modulation amplitude, 20 and 100 K; B = 10 mW power, 10 G modulation amplitude, 3.6 K.
In addition to the intensity increases in the g ∼ 2 region, we observed coincident increases in the intensity of a signal at g ∼ 3.89 indicative of spin–spin coupling between the proximal and distal [4Fe-4S] clusters once the proximal [4Fe-4S] cluster was sufficiently reduced (FigureB; vide infra). These spectral observations were repeatable between samples as well as with the same protein sample, underscoring that this process is consistent between different preparations and can regenerate the starting material. The feature at g = 3.7 in FigureB is likely a spectral artifact (Table S1, Figures and S17). An extended storage period of 24 days for the 90 min phototreated sample in Figure validated that the photoinduced intermediates were stable in liquid nitrogen (∼77 K) without any significant signal decay (Figure S17).
In order to see the maximal signals possible, the NADPH-treated NfnL sample was subjected to an extended illumination of 4.5 h, with the resulting spectra seen in Figure. After illumination, the photogenerated species were replicated with the greatest intensities observed so far, accounting for roughly 53% of the NfnL having an ASQ and showing nearly double-to-triple the previous amounts seen for the radical and proximal [4Fe-4S] cluster species (Figures and S18). The overall intensity of the proximal [4Fe-4S] cluster in the 4.5 h illuminated sample was proportional to the 10x DTH reference sample (FigureB versus FigureB), showing the illumination process can achieve comparably low potentials to that of excess DTH at pH 8.8. Additionally, the observance of extra isotropic intensity flanking the main ASQ feature at 20 K (FigureA) indicated a spin–spin (vide infra; S = 1/2 + 1/2 = 1) interaction? between the proximal [4Fe-4S] cluster and the ASQ radical, further providing evidence for low-potential pathway reduction with a flavin radical rather than an exogenously produced radical (Figure S6).
4.5 h illumination reveals large radical and proximal [4Fe-4S] cluster signals. Panel A shows the Storage (charcoal), Anneal (red), and Light (blue) conditions at 20 K to see relevant paramagnetic species. Storage spectrum reflects the 90 min illuminated conditions for the 2nd phototitration experiment (Figure S17) after 13 days of storage in liquid nitrogen. The Anneal spectrum was taken after melting the sample at room temperature before refreezing in liquid nitrogen, showing regeneration of the preillumination condition. Light spectrum was taken after 4.5 h of illumination at 405 nm and shows much stronger signals for the radical and proximal [4Fe-4S] cluster species, replicating the initial photochemical experiment in Figure with greater intensity. Panel B depicts the Light spectrum in A at 20 K in comparison to that of the 10x DTH reference sample (cyan). The intensity of the proximal [4Fe-4S] cluster signal appears to have doubled or tripled in the 4.5 h experiment relative to the maximal signal observed in previous illumination spectra.
Low-Field Cluster-to-Cluster Coupling Signal Observed When Both
[4Fe-4S] Clusters are Reduced
While acquiring the illumination data, we consistently observed a low-field, semianisotropic feature located at g ∼ 3.89 in multiple samples, but only once both the distal [4Fe-4S] cluster and proximal [4Fe-4S] cluster were reduced (FiguresB and ?). This g ∼ 3.89 feature was notably absent in the preillumination, annealed, and 1x DTH reference NfnL samples, but was present in variable intensity for any sample with significant proximal [4Fe-4S] cluster signal, especially for the 10x DTH reference sample and the 4.5 h illumination sample (FiguresA and S19). This g ∼ 3.89 signal dependence on the proximal [4Fe-4S] cluster strongly implies this feature is a Δm _ S _ = 2 transition stemming from triplet state coupling between the two S = 1/2 spins (S tot = 1/2 + 1/2 = 1) located on the iron–sulfur clusters, especially since they are only 7.4 Å edge-to-edge away (FigureB). ?−? ? Further, the g ∼ 3.89 signal-to-noise ratio was best resolved at 7 K but quickly became weak due to the rapidly relaxing character of the signal (FigureB). Altogether, the g ∼ 3.89 coupling signal reflects a minor but important electromagnetic communication between the two reduced [4Fe-4S] clusters and may be indicative of a directional preference in electron transfer within the low-potential pathway.
Low-field signals arising from reduced [4Fe-4S] cluster–cluster interactions. Panel A depicts the low field data corresponding to Figure A, showing the g ∼ 3.89 feature arises only in the samples with significant amounts of proximal [4Fe-4S] cluster, whereas the feature is absent in samples without proximal [4Fe-4S] cluster reduction (Figure S19). The 1x DTH (gray) and 10x DTH (cyan) values were normalized to match the “junk iron” signal at g ∼ 4.3 in the Dark (lavender) and Light (blue) spectra. Panel B shows the Curie corrected (SignalTemperature) temperature dependence for the low field signal corresponding to the 4.5 h illumination sample in Figure , with the best signal-to-noise at 7 K. Spectral parameters were 10 mW microwave power, modulation amplitude 10 G (panel A) or 6 G (panel B), and temperatures as specified in each graph.*
Discussion
Interpretating the Significance of the High Energy Intermediates
The photoexcitation method utilized here allows for use of a nondegenerative, in situ process that is sufficient to reduce low-potential intermediates that are unstable and difficult to characterize otherwise. In comparison with the 1x and 10x DTH references, the photogenerated proximal [4Fe-4S] cluster can be clearly identified on top of the preexisting distal [4Fe-4S] cluster signal, and based on line width, temperature, and annealing behavior data, the radical is consistent with an ASQ.
Due to the significant spectral overlap of the three iron–sulfur cofactors found in NfnSL and the difficulty in fully populating the low-potential proximal [4Fe-4S] cluster during the former g-value assignments, the initially reported values described a more conventional rhombic shape with g-values of 2.035, 1.955, 1.916. ?,? In this report, the proximal cluster consistently and repeatedly best fit a semiaxial shape with g-values of 1.972, 1.954, and 1.924, and was unexpected considering typical [4Fe-4S] cluster lineshapes are more commonly rhombic in nature. ?,? Both reduction by DTH and the photoexcitation method produced the proximal cluster line shape, showing the signal is an accurate representation of the true paramagnetic species and not an artifact of DTH. These two sets of g-value parameters were compared and fit to different accuracies, where the semiaxial parameters fit significantly better than the previously reported values, indicating this new semiaxial set is more representative of the proximal [4Fe-4S] cluster’s line shape (Figures S3–S6). The signal line shape, and ultimately the electronic structure, of the proximal [4Fe-4S] cluster in NfnL is unique relative to other site-differentiated iron sulfur clusters observed previously, ?−? ? particularly in reference to those with Fe–O ligations. ?−? ? ? ? These results show the proximal [4Fe-4S] cluster has an unconventionally compact S = 1/2 line shape in X-band (∼9.4 GHz) EPR.? The signal shape of the proximal [4Fe-4S] cluster likely stems from the glutamic acid interaction with the Fe, leading to less electron delocalization over the mixed valence Fe^2.5+^ face resulting in a more electronically isolated spin system. ?−? ? Alternatively, the unique line shape for the proximal [4Fe-4S] cluster could be due to magnetic interactions between the two reduced [4Fe-4S] clusters, and that if there were no distal cluster signal at all then the proximal cluster would be much more rhombic in shape. ?,? Although significant magnetic coupling between the two [4Fe-4S] clusters is plausible, this situation is unsupported by the low temperature data which shows that neither the proximal [4Fe-4S] cluster or the distal [4Fe-4S] cluster have any major line shape splitting at lower temperatures,? nor does the low field indicate any single spin system higher than S = 1/2 (Figure).
We have assigned the photogenerated radical species as an ASQ due to the numerous similarities with previously characterized ASQ species, ?,?,?,? the identification of an ASQ at this site in other photoexcitation experiments (TAS), and its inconsistency with the broad NSQ signal observed in the NfnS S-FAD (Figure). ?,? This radical signal completely disappears when the sample is warmed above ∼240 K without any byproduct generation or apparent damage to the protein (FiguresA, S7, S13, S17), and similarly appears in very minor quantities in NfnL when saturated with NADPH (Figure). ?,? Considering these attributes, as well as the reversibility and coaccumulation of the radical and the proximal [4Fe-4S] cluster, we have assigned this ASQ as a catalytically relevant species in the mechanism of electron bifurcation. The illumination control experiments showing the very minor radical and distal [4Fe-4S] cluster signals (∼1% relative to NADPH samples) formed in the oxidized NfnL control further support the relevance of the NADPH (Figures S9–S11, S14, S15). Over the course of 1 h of illumination, the free FAD control photoreduced to the HQ state without evidence of a radical in situ (Figures S14, S15). This photoactivity of free FAD further solidified the mechanistic relevance of NADPH binding near the L-FAD, as without NfnL or NADPH, there is no significant accumulation of any radical species, displaying the ASQ radical is facilitated specifically by the nearby environment surrounding the flavin in NfnL. For the free FAD control to be reduced at pH 8.8 without exogenous chemical reductants, the electrons are either intramolecularly sourced ?−? ? or pulled from the buffer solution, possibly from HEPES itself, as it has a minor amount of redox activity. ?−? ? The in situ photoreduction of the FAD control suggests that the buffer solution is the source of the minor reduction (1%) of the photoexcited L-FAD observed in the oxidized NfnL control (Figures S9–S11).
The observation of an anisotropic low-field (∼1600 G; g ∼ 3.89) signal at low-temperatures (T_opt_ ∼ 7 K) in reduced NfnL was assigned as iron–sulfur cluster spin–spin coupling and is consistent with the edge-to-edge distance between the proximal and distal [4Fe-4S] clusters of about 7.4 Å (FigureB). This low-field signal is emergent from the spin system having S = 1/2 + 1/2 = 1, and as this S = 1 signal depends on both clusters, it is primarily dependent on the reduction of the proximal [4Fe-4S] cluster. ?−? ? ? This is validated by the results in FigureA, showing that the most intense g ∼ 3.89 features occur with the most intense proximal cluster signals and are completely absent when no proximal [4Fe-4S] cluster is observed. This coupling signal is an insight into the finely tuned intercommunication between the downstream distal [4Fe-4S] cluster and the intermediate proximal [4Fe-4S] cluster. This through space coupling between the clusters may be specifically tuned for electron transfer reversibility as well, especially since annealing the sample at higher temperatures reestablishes the original equilibrium (Figure). ?−? ? ?
The Mechanism of Photoexcited Electron Transfer
Flavin photoexcitation is a well-studied process where the oxidized state generates a transient superoxidizing flavin by exciting an electron within a singlet-state electron pair to an excited triplet-state from the HOMO to the LUMO. ?,?,?,?,? This excited state oxidizes a nearby electron donor to fill the newly unoccupied lower energy state, producing a singly reduced semiquinone. In our experiment, the electron source is the nearby NADPH.? This photoinduced mechanism of ASQ generation in free flavin is well understood and essentially only operates via single electron chemistry.?
The mechanism of photoexcited electron transfer in NfnL is hypothesized to generate an ASQ that is capable of transferring its electron to the proximal [4Fe-4S] cluster and to the distal [4Fe-4S] cluster (if it is oxidized; Figure). The injection of 2 electrons into the low-potential pathway is possible since NADPH is a two-electron donor and could feasibly provide both electrons to the low-potential pathway. If the ASQ is photogenerated from a nearby NADPH, then a highly unstable NADP^•^ would be coincidently formed. As we do not see any obvious evidence of an additional pyridine nucleotide radical (e.g., “NADP^•^”)? in the photogenerated spectra (Figures and S6), we can conclude that there must be a rapid 2 electron mechanism that produces the ASQ which transfers to the proximal [4Fe-4S] cluster at 233 K. Literature on the single electron chemistry of NAD(H) calculates the first midpoint reduction potential of NAD to NAD^•^ as −940 mV vs SHE, with the second midpoint reduction potential of NAD^•^ to NADH as +300 mV vs SHE.? These large single electron redox potentials imply that once one electron has been stripped from the nearby NADPH substrate by the photoexcited L-FAD*, the reducing power of the NADP^•^ facilitates rapid electron transfer within the low-potential pathway to repopulate the ASQ on the L-FAD while simultaneously reducing the proximal [4Fe-4S] cluster.
Hypothesized photogeneration mechanism for the observed species based on the data presented in this report. Starting in the oxidized state (bottom left), the L-FAD is excited by 405 nm light (left, blue), producing a highly transient superoxidizing flavin state (top left). This highly oxidizing flavin strips an electron from the nearby NADPH (Figure S20) to transiently reduce the L-FAD to the ASQ state (top right), producing an unobserved NADP• radical (asterisks, blue box). This transient ASQ rapidly transfers its electron down the low-potential pathway to the proximal [4Fe-4S] cluster while being immediately rereduced by the NADP• radical to produce the main paramagnetic intermediates observed in this report (yellow box). Once the sample is brought above 230 K, the sample is able to recombine the electron pair to form the lower-energy diamagnetic L-FAD HQ state observed in UV–visible experiments and ultimately regenerate the resting equilibrium present at the beginning of the experiment. If the distal [4Fe-4S] cluster is already reduced preceding this photogeneration process, then all 3 redox cofactors are paramagnetic and the coupling between the proximal and distal [4Fe-4S] clusters is observable.
Comparison of Photoexcitation Methods
Significant accumulation of an unstable ASQ is challenging to accomplish in TAS experiments due to spectral overlap with NADPH from 340 to 370 nm and at room temperature (293 K). ?,? Not only is the accumulation of the ASQ difficult but also the proportionally small absorbance changes for [4Fe-4S] clusters make tracking ET between the proximal [4Fe-4S] cluster and the ASQ difficult to assess by TAS. The main distinction between the TAS experiment and this study is the constant 405 nm illumination used to excite the L-FAD, as opposed to the pulsed photoexcitation in TAS, which results in continued accumulation of low-potential pathway intermediates (Figure S16). Regardless, the method used in this study is proposed to utilize the same initial photoexcitation mechanism as the TAS experiment that targets and requires the oxidized L-FAD to produce (and constantly accumulate) the ASQ intermediate (Figure S12), contributing to the refined model depicted in Figure. Consistent with the TAS experimental mechanism, the oxidized L-FAD in NfnL absorbs a 405 nm photon and becomes photoexcited to form the highly transient superoxidizing flavin species (“*[L-FAD OX]**”). ?,?,? The superoxidizing flavin can then accept an electron from the nearby NADPH donor to produce an ASQ (Figure S20). After the initial transfer, the highly reducing NADP^•^ provides sufficient pressure to push further downstream reduction of the proximal [4Fe-4S] cluster as the NADP^•^ simultaneously re-reduces the L-FAD to the ASQ to produce the observed results in this report (Figures S6 and S16). If the distal [4Fe-4S] cluster is already reduced, as some initially are in the presence of NADPH (Figure), then the NfnL has unpaired electrons located on all 3 redox cofactors and exhibits the coupling signal between the proximal and distal [4Fe-4S] clusters, located at g ∼ 3.89 (Figure). However, if the distal [4Fe-4S] cluster is in the oxidized form, there is an additional intramolecular electron transfer from the proximal [4Fe-4S] cluster to the distal [4Fe-4S] cluster (Figure S16), leaving just a radical on the flavin as an ASQ along with the reduced distal [4Fe-4S] cluster (Figure S10A).
Advancing a More Robust Model of Electron Bifurcation
By utilizing photoexcitation to initiate electron transfer at cryogenic temperatures, we were able to populate electrons on all of the cofactors in NfnL, including the accumulation of the unstable ASQ species at the L-FAD active site. The resulting EPR data provide insights into the interactions between the low potential pathway cofactors, which have been difficult to ascertain before now due to the inability to visualize the ASQ and proximal [4Fe-4S] cluster simultaneously and with high yields. Because these intermediates are visible and enrichable, we can observe the downhill electron transfer order from the ASQ to the proximal [4Fe-4S] cluster to the distal [4Fe-4S] cluster. Once these reduced states are produced, electron recombination at higher temperatures reveals critical insights into the electron transfer mechanism of confurcation (i.e., the reverse of bifurcation) in NfnL, which are otherwise extremely difficult to visualize. The observation of coupling, notably the low field g = 3.89 feature between the proximal and distal [4Fe-4S] clusters as well as extra g ∼ 2 intensity flanking the radical at g = 2.0033 due to ASQ and proximal [4Fe-4S] cluster interaction, suggests the low-potential pathway has a network of through-space interactions that appear to prime nearby cofactors to sense pathway redox states.
Comparing the results here with photoexcitation performed at room temperature (i.e., TAS) highlights a key difference in the stability and accumulation of ASQ. At cryogenic temperatures, both ASQ and proximal cluster are generated and steadily increase in concentration over 4.5 h (Figure S16), while in the TAS experiment ASQ is rapidly formed but due to back electron transfer with the proximal cluster it has a lifetime of only ∼10 ps. Furthermore, we observed 3 main radical decay phases when the sample is warmed to 230 K (Figure S8) that suggests a more complex mechanism of electron recombination through the low-potential pathway and may include various switching mechanisms to prevent electron transfer in the confurcating direction. One possibility for this temperature dependent behavior may be movement of the Arg333 residue, which is a newly identified determinant of electron bifurcation ability.? When the β-factors between the NADPH-bound and -unbound NfnSL crystal structures are compared, there are notable changes in the region corresponding to the ∼300 to ∼400 amino acid loop (Figure S21). While the R333 residue is not interacting with the L-FAD in either of the crystal structures (and likely represents the low-temperature state), there is ample room for it to move over the O2’, N3, and O4’ moieties of the L-FAD. As modeled in the crystal structure, the R333 interacts with two glutamic acid residues, E335 and E336 (Figure S22). At room temperature, there may be enough conformational sampling that a change in the interaction of R333 from E335/336 to L-FAD may be possible. At cryogenic temperatures, we hypothesize that it is locked out of position (similar to its position in the crystal structure; FigureB and Figure S21) and when the sample is annealed it can move, affecting the reversibility of the electron transfer events across the low potential pathway and/or the stability of the ASQ (Figure S22). These results suggest a previously unobserved gate of electron transfer along the low potential pathway in NfnL that is likely critical for efficient partitioning of electrons during the electron bifurcation mechanism.
Conclusions
The reversible photogeneration of unstable and high-energy paramagnetic intermediates described here allowed for the unambiguous observation of key electron transfer events within the low-potential pathway of Nfn. The recombination of electrons from the low-potential pathway back to the L-FAD suggested a complex and redox-sensitive switch, representing a previously unknown aspect of this elegant mechanism. These studies have described the fundamental characteristics of key bifurcation intermediates and furthered our understanding of how cofactor environments are tuned to initiate and control high-energy electron transfer. Exploiting photoexcitation further opens the door to studying more subtle mechanistic steps within bifurcation, such as proton-coupled electron transfers and the effects of additional nearby amino acid residues on electron transfer (FiguresB, S21, S22). In addition to expanding our mechanistic understanding of Nfn, we demonstrate the utility of applying this photoexcitation method to other FBEB enzymes and flavin-involved electron transfer systems.
Understanding the physical principles that govern the mechanism of electron transfer in biology is key to building better and more robust synthetic systems capable of highly controlled electron transfer. Additionally, the photogeneration of the low-potential species observed in this report may have broader applications, such as exogenous in vitro reductant upconversion for industrial or biosynthetic applications, providing low-potential electrons sourced from NADPH in a controlled manner and pave a way for greater photochemical applications of this unique and thermostable enzyme.
Materials and Methods
All materials and methods were prepared and performed anoxically unless written otherwise.
Production of Recombinant P.f. NfnL and P.f. NfnSL in E. coli BL21 (DE3)
NfnL and NfnSL were grown, overexpressed, purified, and reconstituted from Escherichia coli (E. coli) as previously described. ?,? Previously constructed streptomycin^R^ pCDFDuet-1 plasmids containing either the Pyrococcus furiosus (P.f.; DSM 3638) nfnL gene with an N-terminal Strep-II polypeptide affinity tag (protein sequence GWSHPQFEK) or the nfnS and nfnL genes with an N-terminal Strep-II polypeptide affinity tag were used in this study (GenScript, USA). Since As-purified NfnL and NfnSL lack complete cofactors (8 Fe, 1 FAD for NfnL; 10 Fe, 2 FAD for NfnSL), FAD and Fe reconstitution was carried out as previously described ?,? to ensure enzymes were cofactor-replete.
Protein Concentration Assessed by Rose Bengal Colorimetric Assay
Protein concentrations were quantified aerobically using an adapted version of the rose bengal dye-based colorimetric assay. ?,?,? For the assays, the following stock solutions were prepared: 1 mg/mL rose bengal dye solution, 50% acetic acid quenching solution, and 10 mM potassium phosphate (KP_i_) at pH 7.0. The setup for the rose bengal colorimetric assay was as follows: Add 1800 μL of 10 mM KP_i_ minus the volume of intended protein sample (2–20 μL), to which the protein sample volume (2–20 μL) was added. To begin the assay, 100 μL of the 1 mg/mL rose bengal solution was added to the 1800 μL protein solution. The protein-dye mixture was allowed to react for 3 min. After the incubation time, the protein-dye mixture was quenched with 100 μL of the 50% acetic acid solution (final mixture volume is 2000 μL), and a UV–visible spectrum focused on 560 nm was taken within 5 min of acid-quenching. Finally, the sample-of-interest’s absorbance was compared to a BSA standard curve (11 points; 0 – 20 μg of BSA, increments of 2 μg) to assess protein concentration in sample-of-interest in μg/mL and corrected for any dilutions.
Cofactor Quantitation of NfnL and NfnSL
Cofactor reconstitution was assessed aerobically by iron-counting as previously described, ?,?,?,? as well as flavin-counting as previously described. ?,?,? For NfnL, iron and FAD loading was determined to be 9.0 ± 0.2 iron atoms per NfnL monomer (112 ± 2%) and 1.21 ± 0.03 FAD per NfnL monomer (121 ± 3%). For NfnSL, iron and FAD loading was determined to be 10.2 ± 0.1 iron atoms per NfnSL unit (102 ± 1%) and 1.81 ± 0.04 FAD per NfnSL unit (91 ± 2%).
Preparation of Electron Paramagnetic Resonance Samples
Electron paramagnetic resonance (EPR) samples were prepared as outlined in Table S1 to final volumes of 220 μL with 150 mM HEPES pH 8.80, 200 mM NaCl, and 5% glycerol (“buffer”). Of the total 220 μL, 200 μL was transferred to a quartz EPR tube for spectral acquisition and photoexcitation experiments, and the remaining 20 μL was used for concurrent preillumination UV–visible measurements and protein concentration confirmation by rose bengal colorimetric assay on a Beckman Coulter DU 800 UV–visible spectrophotometer with deuterium (UV) and tungsten (visible) lamps. NADPH was sourced from Research Products International (RPI; Mt. Prospect, IL) and stock concentrations were quantified at 340 nm via UV–visible measurements using the extinction coefficient (ε_340_) of 6.22 mM^–1^ cm^–1^. ?−? ? Sodium dithionite and flavin adenine dinucleotide (FAD) were sourced from MilliporeSigma (Burlington, MA), where they were weighed prior to the preparation of stock solutions. For photoexcitation experiments (NfnL_PE_), NfnL concentrations were prepared to a final concentration of 140 μM with 8.8 mM NADPH final concentration. The relatively high concentration of 8.8 mM NADPH was selected to maintain consistency with previous studies? and to stave off any depletion of NADPH in the case that it was subject to photodegradation by sustained 405 nm illumination. 140 μM NfnL alone, 200 μM FAD and 8.8 mM NADPH were prepared from stock solutions as photoexcitation control experiments. Reduced NfnL (NfnL_Red_) was prepared with 10x dithionite (1 mM DTH) as a reference for the proximal cluster EPR spectrum for the fully reduced state of NfnL. Partially reduced NfnL (NfnL_Dist_) was prepared with 1x dithionite (0.1 mM DTH) as a reference for the distal cluster EPR spectrum in NfnL. As-isolated NfnSL was prepared to a final concentration of 150 μM as a neutral semiquinone (NSQ) reference for flavin radical line width comparisons for organic radicals observed in this study.
Illumination Conditions
Samples were illuminated using a custom, variable temperature cryogenic illumination system comprising a THORLABS (Newton, NJ) 4P4 – 100 mm Integrating Sphere with 4 modular faces outfitted with 2 opposing THORLABS M405L4 1000 mW minimum LED diodes for 2 W illumination at 405 nm excitation wavelength. All illumination was carried out at 233 K (−40 °C) to facilitate electron transfer ?,? and integrating sphere temperature was confirmed using an external K-type thermocouple (Omega HH306A) placed within an EPR tube. NfnL_PE‑I_, FAD, NADPH, and NfnL_Ox_ were all illuminated at 233 K for 1 h before spectral acquisition. NfnL_PE‑T_ was illuminated at 233 K for various times, at which spectral acquisition was taken at the conclusion of each illumination interval.
EPR Data Acquisition, Parameters, and Data Processing
All EPR spectra were collected on an ELEXSYS E500 continuous-wave (CW) X-band (∼9.4 GHz) spectrometer (Bruker; Billerica, MA) outfitted with a super high-Q perpendicular-mode resonator and an in-cavity, cryogen-free variable temperature helium system (ColdEdge Technologies; Allentown, PA) paired to a Mercury iTC temperature controller (Oxford Instruments; Abingdon, UK). Variable parameters such as temperature, modulation amplitude (MA), and microwave power are specified for relevant traces in each figure and range from 3.5 to 230 K, 1 to 10 G, and 1 to 10 mW, respectively. Parameters applying to all collected spectra were: Modulation frequency of 100 kHz; gain of 60 dB; conversion time of 81.92 ms; 1024 points per spectrum; 5 scan averaging; field sweep 2000 G; and depending on the region of interest, center field of 3350 G (g ∼ 2 region) or 1400 G (g ∼ 4). Microwave frequency was around 9.37 GHz for all spectra, and each spectrum’s g-value plot was corrected by its corresponding microwave frequency recorded during spectral acquisition. EPR data was exported and initially processed using an export program that calculates g-values in MATLAB R2020b. Data were processed in OriginPro 2023, and S = 1/2 (g ∼ 2 region) spectra were simulated using EasySpin’s pepper function in MATLAB R2020b. ?,? ASQ radical spins were quantified in comparison to those of a 100 μM Cu(II) TEA standard at a volume of 200 μL recorded at the same spectral acquisition parameters and temperatures noted for the corresponding radical spectra (Figure S18).
Protein model figures were produced using PyMOL software (The PyMOL Molecular Graphics System, Version 1.2r3pre, Schrödinger, LLC.) which were further processed by layering images in Microsoft PowerPoint. Figures with EPR data were generated by using OriginPro 2023.
Supplementary Material
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Buckel W.Thauer R. K.Energy Conservation via Electron Bifurcating Ferredoxin Reduction and Proton/Na+ Translocating Ferredoxin Oxidation Evol. Asp. Bioenerg. Syst.2013182729411310.1016/j.bbabio.2012.07.00222800682 · doi ↗ · pubmed ↗
- 2Peters J. W.Miller A.-F.Jones A. K.King P. W.Adams M. W.Electron Bifurcation Biocatal. Biotransformation Bioinorg. Chem.20163114615210.1016/j.cbpa.2016.03.00727016613 · doi ↗ · pubmed ↗
- 3Mitchell P.The Protonmotive Q Cycle: A General Formulation FEBS Lett.197559213713910.1016/0014-5793(75)80359-01227927 · doi ↗ · pubmed ↗
- 4Mitchell P.Possible Molecular Mechanisms of the Protonmotive Function of Cytochrome Systems J. Theor. Biol.197662232736710.1016/0022-5193(76)90124-7186667 · doi ↗ · pubmed ↗
- 5Herrmann G.Jayamani E.Mai G.Buckel W.Energy Conservation via Electron-Transferring Flavoprotein in Anaerobic Bacteria J. Bacteriol.2008190378479110.1128/JB.01422-0718039764 PMC 2223574 · doi ↗ · pubmed ↗
- 6Yuly J. L.Zhang P.Beratan D. N.Energy Transduction by Reversible Electron Bifurcation Curr. Opin. Electrochem.20212910076710.1016/j.coelec.2021.100767 · doi ↗
- 7Wise C. E.Ledinina A. E.Yuly J. L.Artz J. H.Lubner C. E.The Role of Thermodynamic Features on the Functional Activity of Electron Bifurcating Enzymes Biochim. Biophys. Acta BBA - Bioenerg.20211862414837710.1016/j.bbabio.2021.14837733453185 · doi ↗ · pubmed ↗
- 8Buckel W.Thauer R. K.Flavin-Based Electron Bifurcation, A New Mechanism of Biological Energy Coupling Chem. Rev.201811873862388610.1021/acs.chemrev.7b 0070729561602 · doi ↗ · pubmed ↗
