Exploring Electron/Proton Transfer and Conformational Changes in the Nitrogenase MoFe Protein and FeMo-cofactor Through Cryoreduction/EPR Measurements

Analysis of the Geometric and Electronic Structure of Spin-Coupled Iron-Sulfur Dimers with Broken-Symmetry DFT: Implications for FeMoco

The open-shell electronic structure of iron–sulfur clusters presents considerable challenges to quantum chemistry, with the complex iron–molybdenum cofactor (FeMoco) of nitrogenase representing perhaps the ultimate challenge for either wavefunction or density functional theory. While broken-symmetry density functional theory has seen some success in describing the electronic structure of such cofactors, there is a large exchange–correlation functional dependence in calculations that is not fully understood. In this work, we present a geometric benchmarking test set, FeMoD11, of synthetic spin-coupled Fe–Fe and Mo–Fe dimers, with relevance to the molecular and electronic structure of the Mo-nitrogenase FeMo cofactor. The reference data consists of high-resolution crystal structures of metal dimer compounds in different oxidation states. Multiple density functionals are tested on their ability to reproduce the local geometry, specifically the Fe–Fe/Mo–Fe distance, for both antiferromagnetically coupled and ferromagnetically coupled dimers via the broken-symmetry approach. The metal–metal distance is revealed not only to be highly sensitive to the amount of exact exchange in the functional but also to the specific exchange and correlation functionals. For the antiferromagnetically coupled dimers, the calculated metal–metal distance correlates well with the covalency of the bridging metal–ligand bonds, as revealed via the corresponding orbital analysis, Hirshfeld S/Fe charges, and Fe–S Mayer bond order. Superexchange via bridging ligands is expected to be the dominant interaction in these dimers, and our results suggest that functionals that predict accurate Fe–Fe and Mo–Fe distances describe the overall metal–ligand covalency more accurately and in turn the superexchange of these systems. The best performing density functionals of the 16 tested for the FeMoD11 test set are revealed to be either the nonhybrid functionals r²SCAN and B97-D3 or hybrid functionals with 10–15% exact exchange: TPSSh and B3LYP*. These same four functionals are furthermore found to reproduce the high-resolution X-ray structure of FeMoco well according to quantum mechanics/molecular mechanics (QM/MM) calculations. Almost all nonhybrid functionals systematically underestimate Fe–Fe and Mo–Fe distances (with r²SCAN and B97-D3 being the sole exceptions), while hybrid functionals with >15% exact exchange (including range-separated hybrid functionals) overestimate them. The results overall suggest r²SCAN, B97-D3, TPSSh, and B3LYP* as accurate density functionals for describing the electronic structure of iron–sulfur clusters in general, including the complex FeMoco cluster of nitrogenase.

The electronic structure of FeV-cofactor in vanadium-dependent nitrogenase

The electronic structure of the active-site metal cofactor (FeV-cofactor) of resting-state V-dependent nitrogenase has been an open question, with earlier studies indicating that it exhibits a broad S = 3/2 EPR signal (Kramers state) having g values of ∼4.3 and 3.8, along with suggestions that it contains metal-ions with valencies [1V3+, 3Fe3+, 4Fe2+]. In the present work, genetic, biochemical, and spectroscopic approaches were combined to reveal that the EPR signals previously assigned to FeV-cofactor do not correlate with active VFe-protein, and thus cannot arise from the resting-state of catalytically relevant FeV-cofactor. It, instead, appears resting-state FeV-cofactor is either diamagnetic, S = 0, or non-Kramers, integer-spin (S = 1, 2 etc.). When VFe-protein is freeze-trapped during high-flux turnover with its natural electron-donating partner Fe protein, conditions which populate reduced states of the FeV-cofactor, a new rhombic S = 1/2 EPR signal from such a reduced state is observed, with g = [2.18, 2.12, 2.09] and showing well-defined 51V (I = 7/2) hyperfine splitting, a iso = 110 MHz. These findings indicate a different assignment for the electronic structure of the resting state of FeV-cofactor: S = 0 (or integer-spin non-Kramers state) with metal-ion valencies, [1V3+, 4Fe3+, 3Fe2+]. Our findings suggest that the V3+ does not change valency throughout the catalytic cycle.

Shallow Distance Dependence for Proton-Coupled Tyrosine Oxidation in Oligoproline Peptides

We have explored the kinetic effect of increasing electron transfer (ET) distance in a biomimetic, proton-coupled electron-transfer (PCET) system. Biological ET often occurs simultaneously with proton transfer (PT) in order to avoid the high-energy, charged intermediates resulting from the stepwise transfer of protons and electrons. These concerted proton-electron-transfer (CPET) reactions are implicated in numerous biological ET pathways. In many cases, PT is coupled to long-range ET. While many studies have shown that the rate of ET is sensitive to the distance between the electron donor and acceptor, extensions to biological CPET reactions are sparse. The possibility of a unique ET distance dependence for CPET reactions deserves further exploration, as this could have implications for how we understand biological ET. We therefore explored the ET distance dependence for the CPET oxidation of tyrosine in a model system. We prepared a series of metallopeptides with a tyrosine separated from a Ru(bpy)32+ complex by an oligoproline bridge of increasing length. Rate constants for intramolecular tyrosine oxidation were measured using the flash-quench transient absorption technique in aqueous solutions. The rate constants for tyrosine oxidation decreased by 125-fold with three added proline residues between tyrosine and the oxidant. By comparison, related intramolecular ET rate constants in very similar constructs were reported to decrease by 4-5 orders of magnitude over the same number of prolines. The observed shallow distance dependence for tyrosine oxidation is proposed to originate in part from the requirement for stronger oxidants, leading to a smaller hole-transfer effective tunneling barrier height. The shallow distance dependence observed here and extensions to distance-dependent CPET reactions have potential implications for long-range charge transfers.

The Spectroscopy of Nitrogenases

Nitrogenases are responsible for biological nitrogen fixation, a crucial step in the biogeochemical nitrogen cycle. These enzymes utilize a two-component protein system and a series of iron-sulfur clusters to perform this reaction, culminating at the FeMco active site (M = Mo, V, Fe), which is capable of binding and reducing N2 to 2NH3. In this review, we summarize how different spectroscopic approaches have shed light on various aspects of these enzymes, including their structure, mechanism, alternative reactivity, and maturation. Synthetic model chemistry and theory have also played significant roles in developing our present understanding of these systems and are discussed in the context of their contributions to interpreting the nature of nitrogenases. Despite years of significant progress, there is still much to be learned from these enzymes through spectroscopic means, and we highlight where further spectroscopic investigations are needed.

Spectroscopic Description of the E1 State of Mo Nitrogenase Based on Mo and Fe X-ray Absorption and Mössbauer Studies

Mo nitrogenase (N2ase) utilizes a two-component protein system, the catalytic MoFe and its electron-transfer partner FeP, to reduce atmospheric dinitrogen (N₂) to ammonia (NH₃). The FeMo cofactor contained in the MoFe protein serves as the catalytic center for this reaction and has long inspired model chemistry oriented toward activating N₂. This field of chemistry has relied heavily on the detailed characterization of how Mo N2ase accomplishes this feat. Understanding the reaction mechanism of Mo N2ase itself has presented one of the most challenging problems in bioinorganic chemistry because of the ephemeral nature of its catalytic intermediates, which are difficult, if not impossible, to singly isolate. This is further exacerbated by the near necessity of FeP to reduce native MoFe, rendering most traditional means of selective reduction inept. We have now investigated the first fundamental intermediate of the MoFe catalytic cycle, E₁, as prepared both by low-flux turnover and radiolytic cryoreduction, using a combination of Mo Kα high-energy-resolution fluorescence detection and Fe K-edge partial-fluorescence-yield X-ray absorption spectroscopy techniques. The results demonstrate that the formation of this state is the result of an Fe-centered reduction and that Mo remains redox-innocent. Furthermore, using Fe X-ray absorption and ⁵⁷Fe Mössbauer spectroscopies, we correlate a previously reported unique species formed under cryoreducing conditions to the natively formed E₁ state through annealing, demonstrating the viability of cryoreduction in studying the catalytic intermediates of MoFe.

Short-lived neutral FMN and FAD semiquinones are transient intermediates in cryo-reduced yeast NADPH-cytochrome P450 reductase

The electron configuration of flavin cofactors, FMN and FAD, is a critical factor governing the reactivity of NADPH-cytochrome P450 reductase (CPR). The current view of electron transfer by the mammalian CPR, based on equilibrium redox potentials of the flavin cofactors, is that the two electron-reduced FMN hydroquinone (FMNH2), rather than one electron-reduced FMN semiquinone, serves as electron donor to the terminal protein acceptors. However, kinetic and thermodynamic studies on the CPR species originated from different organisms have shown that redox potentials measured at distinct electron transfer steps differ from redox potentials determined by equilibrium titration. Collectively, previous observations suggest that the short-lived transient semiquinone species may carry electrons in diflavin reductases. In this work, we have investigated spectroscopic properties of the CPR-bound FAD and FMN reduced at 77 K by radiolytically-generated thermalized electrons. Using UV-vis spectroscopy, we demonstrated that upon cryo-reduction of oxidized yeast CPR (yCPR) containing an equimolar ratio of both FAD and FMN, or FAD alone, neutral semiquinones were trapped at 77 K. During annealing at the elevated temperatures, unstable short-lived neutral semiquinones relaxed to spectroscopically distinct air-stable neutral semiquinones. This transition was independent of pH within the 6.0-10.7 range. Our data on yeast CPR are in line with the previous observations of others that the flavin short-lived transient semiquinone intermediates may have a role in the electron transfer by CPR at physiological conditions.

X-ray Magnetic Circular Dichroism Spectroscopy Applied to Nitrogenase and Related Models: Experimental Evidence for a Spin-Coupled Molybdenum(III) Center

Nitrogenase enzymes catalyze the reduction of atmospheric dinitrogen to ammonia utilizing a Mo‐7Fe‐9S‐C active site, the so‐called FeMoco cluster. FeMoco and an analogous small‐molecule (Et₄N)[(Tp)MoFe₃S₄Cl₃] cubane have both been proposed to contain unusual spin‐coupled Mo^III sites with an S(Mo)=1/2 non‐Hund configuration at the Mo atom. Herein, we present Fe and Mo L₃‐edge X‐ray magnetic circular dichroism (XMCD) spectroscopy of the (Et₄N)[(Tp)MoFe₃S₄Cl₃] cubane and Fe L_2,3‐edge XMCD spectroscopy of the MoFe protein (containing both FeMoco and the 8Fe‐7S P‐cluster active sites). As the P‐clusters of MoFe protein have an S=0 total spin, these are effectively XMCD‐silent at low temperature and high magnetic field, allowing for FeMoco to be selectively probed by Fe L_2,3‐edge XMCD within the intact MoFe protein. Further, Mo L₃‐edge XMCD spectroscopy of the cubane model has provided experimental support for a local S(Mo)=1/2 configuration, demonstrating the power and selectivity of XMCD.

Ultrafast Spectroscopic Dynamics of Quinacrine-Riboflavin Binding Protein Interactions

Redox active cofactors play a dynamic role inside protein binding active sites because the amino acids responsible for binding participate in electron transfer (ET) reactions. Here, we use femtosecond transient absorption (FsTA) spectroscopy to examine the ultrafast ET between quinacrine (Qc), an antimalarial drug with potential anticancer activity, and riboflavin binding protein (RfBP) with a known K_d = 264 nM. Steady-state absorption reveals a ∼ 10 nm red-shift in the ground state when QcH₃²⁺ is titrated with RfBP, and a Stern-Volmer analysis shows ∼84% quenching and a blue-shift of the QcH₃²⁺ photoluminescence to form a 1:1 binding ratio of the QcH₃²⁺-RfBP complex. Upon selective photoexcitation of QcH₃²⁺ in the QcH₃²⁺-RfBP complex, we observe charge separation in 7 ps to form ¹[QcH_{3_red}^•+-RfBP^•+], which persists for 138 ps. The FsTA spectra show the spectroscopic identification of QcH_{3_red}^•+, determined from spectroelectrochemical measurements in DMSO. We correlate our results to literature and report lifetimes that are 10-20× slower than the natural riboflavin, Rf-RfBP, complex and are oxygen independent. Driving force (ΔG) calculations, corrected for estimated dielectric constants for protein hydrophobic pockets, and Marcus theory depict a favorable one-electron ET process between QcH₃²⁺ and nearby redox active tyrosine (Tyr) or tryptophan (Trp) residues.