Resolving the structure of the E1 state of Mo nitrogenase through Mo and Fe K-edge EXAFS and QM/MM calculations

Understanding the P-Cluster of Vanadium Nitrogenase: an EPR and XAS Study of the Holo vs. Apo Forms of the Enzyme

The catalytic moiety of nitrogenases contains two complex metalloclusters: The M-cluster (also called cofactor), where the catalytic reduction of substrates takes place, and the [Fe8S7] P-cluster responsible for electron transfer. Due to discrepancies between crystallography and EPR spectroscopy, the exact structure of the P-cluster in the VFe protein remains a topic of debate. Herein, we use an apo-form of VFe (which retains the P-cluster but lacks the FeVco) to study the VFe P-cluster. SDS-PAGE and NativePAGE showed a heterogeneous composition of the VFe and the apo-VFe samples with the presence of α1β2δ2 and α1β2 complexes. The parallel mode EPR measurements of IDS oxidized MoFe, apo-MoFe, and VFe samples reveal a signal at g=12 associated with the two-electron oxidized state of the P-cluster (P2+) for all three samples, albeit with different intensities. In contrast, no P2+ was observed for IDS oxidized apo-VFe. Additionally, comparisons between apo-MoFe, apo-VFe and the model complex (NBu4)2[Fe4S4(SPh)4] via EXAFS measurements showed that apo-VFe does not contain a fully formed [Fe8S7] P-cluster, but rather is comprised of fragmented iron-sulfur clusters. Our results point to a possible variation in the structure of the P-cluster in the different forms of the nitrogenase.

Investigating the Molybdenum Nitrogenase Mechanistic Cycle Using Spectroelectrochemistry

Molybdenum nitrogenase plays a crucial role in the biological nitrogen cycle by catalyzing the reduction of dinitrogen (N2) to ammonia (NH3) under ambient conditions. However, the underlying mechanisms of nitrogenase catalysis, including electron and proton transfer dynamics, remain only partially understood. In this study, we covalently attached molybdenum nitrogenase (MoFe) to gold electrodes and utilized surface-enhanced infrared absorption spectroscopy (SEIRA) coupled with electrochemistry techniques to investigate its catalytic mechanism. Our biohybrid system enabled electron transfer via a mild mediator, likely mimicking the natural electron flow through the P-cluster to FeMoco, the enzyme's active site. For the first time, we experimentally observed both terminal and bridging S-H stretching frequencies, resulting from the protonation of bridging sulfides in FeMoco during turnover conditions providing direct evidence of their role in catalysis. These experimental observations are further supported by QM/MM calculations. Additionally, we investigated CO inhibition, demonstrating both CO binding and unbinding dynamics under electrochemical conditions. These insights not only advance our understanding of the mechanistic cycle of molybdenum nitrogenase but also establish a foundation for studying alternative nitrogenases, including vanadium and iron nitrogenases.

The nitrogenase mechanism: new roles for the dangler?

Dangler sites protruding from a core metallocluster were introduced into the bioinorganic lexicon in 2000 by R.D. Britt and co-workers in an analysis of the tetramanganese oxygen-evolving cluster in photosystem II. In this perspective, we consider whether analogous dangler sites could participate in the mechanism of dinitrogen reduction by nitrogenase. Two possible roles for dynamic danglers in the active site FeMo cofactor are highlighted that might occur transiently during turnover. The first role for a dangler involves the S2B belt sulfur associated with displacement by carbon monoxide and other ligands, while the second dangler role could involve the entire cluster upon displacement of the His-  α  442 side chain to the molybdenum by a free carboxyl group of the homocitrate ligand. To assess whether waters might be able to interact with the cofactor, a survey of small ligands (water and alkali metal ions) contacting [4Fe4S] clusters in synthetic compounds and proteins was conducted. This survey reveals a preference for these sites to pack over the centers of 2Fe2S rhombs. Waters are excluded from the S2B site in the resting state of nitrogenase, suggesting it is unlikely that water molecules coordinate to the FeMo cofactor during catalysis. While alkali metal ions are found to generally influence the properties of catalysts for dinitrogen reduction, no convincing evidence was found that any of the waters near the FeMo cofactor could instead be sodium or potassium ions. Dangler sites, if they exist in the nitrogenase mechanism, are likely formed transiently by localized changes to the resting-state FeMo cofactor structure.

Sulfide release and rebinding in the mechanism for nitrogenase

Nitrogenases are the only enzymes that activate the strong triple bond in N2. The mechanism for the activation has been very difficult to determine in spite of decades of work. In previous modeling studies it has been suggested that the mechanism for nitrogen activation starts out by four pre-activation steps (A0-A4) before catalysis. That suggestion led to excellent agreement with experimental Elecrtron Paramagnetic Resonance (EPR) observations in the step where N2 becomes protonated (E4). An important part of the pre-activation is that a sulfide is released. In the present paper, the details of the pre-activation are modeled, including the release of the sulfide. Several possible transition states for the release have been obtained. An A4(E0) state is reached which is very similar to the E4 state. For completeness, the steps going back from A4(E0) to A0 after catalysis are also modeled, including the insertion of a sulfide.

The activating capture of N2 at the active site of Mo-nitrogenase

Dinitrogen is inherently inert. This report describes detailed density functional calculations (with a 485+ atom model) of mechanistic steps by which the enzyme nitrogenase activates unreactive N2 at the intact active site FeMo-co, to form a key intermediate with bound HNNH. This mechanism does not bind N2 first and then add H atoms, but rather captures N2 ('N2-ready') that diffuses in through the substrate channel and enters a strategic gallery of H atom donors in the reaction zone, between Fe2 and Fe6. This occurs at the E4 stage of the complete mechanism. Exploration of possible reactions of N2 in this space leads to the conclusion that the first reaction step is transfer of H on Fe7 to one end of N2-ready, soon followed by Fe-N bond formation, and then a second H transfer from bridging S2BH to the other N. Two H-N bonds and one or two N-Fe bonds are formed, in some cases with a single transition state. The variable positions and orientations of N2-ready lead to various reaction trajectories and products. The favourable products resulting from this capture, judged by the criteria of reaction energies, reaction barriers, and mechanistic competence for further hydrogenation reactions in the nitrogenase cycle, have Fe2-NH-NH bonding. The trajectory of one N2 capture reaction is described in detail, and calculations that separate the H atom component and the 'heavy atom' components of the classical activation energy are described, in the context of possible H atom tunneling in the activation of N2-ready. I present arguments for the activation of N2 by the pathway of concerted hydrogenation and binding of N2-ready, alternative to the commonly assumed pathway of binding N2 first, with subsequent hydrogenation. The active site of nitrogenase is well primed for the thermodynamic and kinetic advantages of N2 capture.

The energetics of N2 reduction by vanadium containing nitrogenase

The main class of nitrogenases has a molybdenum in its cofactor. A mechanism for Mo-nitrogenase has recently been described. In the present study, another class of nitrogenases has been studied, the one with a vanadium instead of a molybdenum in its cofactor. It is generally believed that these classes use the same general mechanism to activate nitrogen. The same methodology has been used here as the one used for Mo-nitrogenase. N2 activation is known to occur after four reductions in the catalytic cycle, in the E4 state. The main features of the mechanism for Mo-nitrogenase found in the previous study are an activation process in four steps prior to catalysis, the release of a sulfide during the activation steps and the formation of H2 from two hydrides in E4, just before N2 is activated. The same features have been found here for V-nitrogenase. A difference is that five steps are needed in the activation process, which explains why the ground state of V-nitrogenase is a triplet (even number) and the one for Mo-nitrogenase is a quartet (odd number). The reason an additional step is needed for V-nitrogenase is that V3+ can be reduced to V2+, in contrast to the case for Mo3+ in Mo-nitrogenase. The fact that V3+ is Jahn-Teller active has important consequences. N2H2 is formed in E4 with reasonably small barriers.

H2 formation from the E2-E4 states of nitrogenase

Nitrogenase is the only enzyme that can cleave the strong triple bond in N2, making nitrogen available for biological lifeforms. The active site is a MoFe7S9C cluster (the FeMo cluster) that binds eight electrons and protons during one catalytic cycle, giving rise to eight intermediate states E0-E7. It is experimentally known that N2 binds to the E4 state and that H2 is a compulsory byproduct of the reaction. However, formation of H2 is also an unproductive side reaction that should be avoided, especially in the early steps of the reaction mechanism (E2 and E3). Here, we study the formation of H2 for various structural interpretations of the E2-E4 states using combined quantum mechanical and molecular mechanical (QM/MM) calculations and four different density-functional theory methods. We find large differences in the predictions of the different methods. B3LYP strongly favours protonation of the central carbide ion and H2 cannot form from such structures. On the other hand, with TPSS, r2SCAN and TPSSh, H2 formation is strongly exothermic for all structures and En and therefore need strict kinetic control to be avoided. For the E2 state, the kinetic barriers for the low-energy structures are high enough to avoid H2 formation. However, for both the E3 and E4 states, all three methods predict that the best structure has two hydride ions bridging the same pair of Fe ions (Fe2 and Fe6) and these two ions can combine to form H2 with an activation barrier of only 29-57 kJ mol-1, corresponding to rates of 7 × 102 to 5 × 107 s-1, i.e. much faster than the turnover rate of the enzyme (1-5 s-1). We have also studied H-atom movements within the FeMo cluster, showing that the various protonation states can quite freely be interconverted (activation barriers of 12-69 kJ mol-1).

Nitrogenase beyond the Resting State: A Structural Perspective

Nitrogenases have the remarkable ability to catalyze the reduction of dinitrogen to ammonia under physiological conditions. How does this happen? The current view of the nitrogenase mechanism focuses on the role of hydrides, the binding of dinitrogen in a reductive elimination process coupled to loss of dihydrogen, and the binding of substrates to a binuclear site on the active site cofactor. This review focuses on recent experimental characterizations of turnover relevant forms of the enzyme determined by cryo-electron microscopy and other approaches, and comparison of these forms to the resting state enzyme and the broader family of iron sulfur clusters. Emerging themes include the following: (i) The obligatory coupling of protein and electron transfers does not occur in synthetic and small-molecule iron-sulfur clusters. The coupling of these processes in nitrogenase suggests that they may involve unique features of the cofactor, such as hydride formation on the trigonal prismatic arrangement of irons, protonation of belt sulfurs, and/or protonation of the interstitial carbon. (ii) Both the active site cofactor and protein are dynamic under turnover conditions; the changes are such that more highly reduced forms may differ in key ways from the resting-state structure. Homocitrate appears to play a key role in coupling cofactor and protein dynamics. (iii) Structural asymmetries are observed in nitrogenase under turnover-relevant conditions by cryo-electron microscopy, although the mechanistic relevance of these states (such as half-of-sites reactivity) remains to be established.

On the Shoulders of Giants-Reaching for Nitrogenase

Only a single enzyme system-nitrogenase-carries out the conversion of atmospheric N2 into bioavailable ammonium, an essential prerequisite for all organismic life. The reduction of this inert substrate at ambient conditions poses unique catalytic challenges that strain our mechanistic understanding even after decades of intense research. Structural biology has added its part to this greater tapestry, and in this review, I provide a personal (and highly biased) summary of the parts of the story to which I had the privilege to contribute. It focuses on the crystallographic analysis of the three isoforms of nitrogenases at high resolution and the binding of ligands and inhibitors to the active-site cofactors of the enzyme. In conjunction with the wealth of available biochemical, biophysical, and spectroscopic data on the protein, this has led us to a mechanistic hypothesis based on an elementary mechanism of repetitive hydride formation and insertion.

Protonation of Homocitrate and the E1 State of Fe-Nitrogenase Studied by QM/MM Calculations

Nitrogenase is the only enzyme that can cleave the strong triple bond in N2, making nitrogen available for biological life. There are three isozymes of nitrogenase, differing in the composition of the active site, viz., Mo, V, and Fe-nitrogenase. Recently, the first crystal structure of Fe-nitrogenase was presented. We have performed the first combined quantum mechanical and molecular mechanical (QM/MM) study of Fe-nitrogenase. We show with QM/MM and quantum-refinement calculations that the homocitrate ligand is most likely protonated on the alcohol oxygen in the resting E0 state. The most stable broken-symmetry (BS) states are the same as for Mo-nitrogenase, i.e., the three Noodleman BS7-type states (with a surplus of β spin on the eighth Fe ion), which maximize the number of nearby antiferromagnetically coupled Fe-Fe pairs. For the E1 state, we find that protonation of the S2B μ2 belt sulfide ion is most favorable, 14-117 kJ/mol more stable than structures with a Fe-bound hydride ion (the best has a hydride ion on the Fe2 ion) calculated with four different density-functional theory methods. This is similar to what was found for Mo-nitrogenase, but it does not explain the recent EPR observation that the E1 state of Fe-nitrogenase should contain a photolyzable hydride ion. For the E1 state, many BS states are close in energy, and the preferred BS state differs depending on the position of the extra proton and which density functional is used.

Catalysis and structure of nitrogenases

In providing bioavailable nitrogen as building blocks for all classes of biomacromolecules, biological nitrogen fixation is an essential process for all organismic life. Only a single enzyme, nitrogenase, performs this task at ambient conditions and with ATP as an energy source. The assembly of the complex iron-sulfur enzyme nitrogenase and its catalytic mechanism remains a matter of intense study. Recent progress in the structural analysis of the three known isoforms of nitrogenase-differentiated primarily by the heterometal in their active site cofactor-has revealed a degree of structural plasticity of these clusters that suggest two distinct binding sites for substrates and reaction intermediates. A mechanistic proposal based on this finding integrates most of the available experimental data. Furthermore, the first applications of high-resolution cryo-electron microscopy have highlighted further dynamic conformational changes. Structures obtained under turnover conditions support the proposed alternating half-site reactivity in the C2-symmetric nitrogenase complex.

Multiscale QM/MM modelling of catalytic systems with ChemShell

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods are a powerful computational tool for the investigation of all forms of catalysis, as they allow for an accurate description of reactions occurring at catalytic sites in the context of a complicated electrostatic environment. The scriptable computational chemistry environment ChemShell is a leading software package for QM/MM calculations, providing a flexible, high performance framework for modelling both biomolecular and materials catalysis. We present an overview of recent applications of ChemShell to problems in catalysis and review new functionality introduced into the redeveloped Python-based version of ChemShell to support catalytic modelling. These include a fully guided workflow for biomolecular QM/MM modelling, starting from an experimental structure, a periodic QM/MM embedding scheme to support modelling of metallic materials, and a comprehensive set of tutorials for biomolecular and materials modelling.

The E3 state of FeMoco: one hydride, two hydrides or dihydrogen?

Hydrides are present in the reduced states of the iron-molybdenum cofactor (FeMoco) of Mo nitrogenase and are believed to play a key mechanistic role in the dinitrogen reduction reaction catalyzed by the enzyme. Two hydrides are present in the E4 state according to 1H ENDOR and there is likely a single hydride in the E2 redox state. The 2-hydride E4 state has been experimentally observed to bind N2 and it has been speculated that E3 may bind N2 as well. However, the E3 state has not been directly observed and very little is known about its molecular and electronic structure or reactivity. In recent computational studies, we have explored the energy surfaces of the E2 and E4 by QM/MM modelling, and found that the most stable hydride isomers contain bridging or partially bridging hydrides with an open protonated belt sulfide-bridge. In this work we systematically explore the energy surface of the E3 redox state, comparing single hydride and two-hydride isomers with varying coordination and bridging vs. terminal sulfhydryl groups. We also include a model featuring a triply protonated carbide. The results are only mildly dependent on the QM-region size. The three most stable E3 isomers at the r2SCAN level of theory have in common: an open belt sulfide-bridge (terminal sulfhydryl group on Fe6) and either 2 bridging hydrides (between Fe2 and Fe6), 1 bridging-1-terminal hydride (around Fe2 and Fe6) or a dihydrogen ligand bound at the Fe2 site. Analyzing the functional dependency of the results, we find that functionals previously found to predict accurate structures of spin-coupled Fe/Mo dimers and FeMoco (TPSSh, B97-D3, r2SCAN, and B3LYP*) are in generally good agreement about the stability of these 3 E3 isomers. However, B3LYP*, similar to its parent B3LYP method, predicts a triply protonated carbide isomer as the most stable isomer, an unlikely scenario in view of the lack of experimental evidence for carbide protonation occurring in reduced FeMoco states. Distinguishing further between the 3 hydride isomers is difficult and this flexible coordination nature of hydrides suggests that multiple hydride isomers could be present during experimental conditions. N2 binding was explored and resulted in geometries with 2 bridging hydrides and N2 bound to either Fe2 or Fe6 with a local low-spin state on the Fe. N2 binding is predicted to be mildly endothermic, similar to the E2 state, and it seems unlikely that the E3 state is capable of binding N2.

Connecting the geometric and electronic structures of the nitrogenase iron-molybdenum cofactor through site-selective 57Fe labelling

Understanding the chemical bonding in the catalytic cofactor of the Mo nitrogenase (FeMo-co) is foundational for building a mechanistic picture of biological nitrogen fixation. A persistent obstacle towards this goal has been that the 57Fe-based spectroscopic data-although rich with information-combines responses from all seven Fe sites, and it has therefore not been possible to map individual spectroscopic responses to specific sites in the three-dimensional structure. Here we have addressed this challenge by incorporating 57Fe into a single site of FeMo-co. Spectroscopic analysis of the resting state informed on the local electronic structure of the terminal Fe1 site, including its oxidation state and spin orientation, and, in turn, on the spin-coupling scheme for the entire cluster. The oxidized resting state and the first intermediate in nitrogen fixation were also characterized, and comparisons with the resting state provided molecular-level insights into the redox chemistry of FeMo-co.

Evolution of Stabilized 1T-MoS2 by Atomic-Interface Engineering of 2H-MoS2 /Fe-Nx towards Enhanced Sodium Ion Storage

Metallic conductive 1T phase molybdenum sulfide (MoS₂ ) has been identified as promising anode for sodium ion (Na⁺ ) batteries, but its metastable feature makes it difficult to obtain and its restacking during the charge/discharge processing result in part capacity reversibility. Herein, a synergetic effect of atomic-interface engineering is employed for constructing 2H-MoS₂ layers assembled on single atomically dispersed Fe-N-C (SA Fe-N-C) anode material that boosts its reversible capacity. The work-function-driven-electron transfer occurs from SA Fe-N-C to 2H-MoS₂ via the Fe-S bonds, which enhances the adsorption of Na⁺ by 2H-MoS₂ , and lays the foundation for the sodiation process. A phase transfer from 2H to 1T/2H MoS₂ with the ferromagnetic spin-polarization of SA Fe-N-C occurs during the sodiation/desodiation process, which significantly enhances the Na⁺ storage kinetics, and thus the 1T/2H MoS₂ /SA Fe-N-C display a high electronic conductivity and a fast Na⁺ diffusion rate.

Can the E1 state in nitrogenase tell if there is an activation process prior to catalysis?

Model calculations have been performed for the singly reduced ground state of Mo-nitrogenase, usually termed E₁. Contradictory conclusions have been reached in two recent experimental studies. In a study based on EPR, it was concluded that there is a bridging hydride in E₁, while in an X-ray study it was concluded that there is no hydride in E₁. Therefore, the EPR study implies that there is an oxidation of the cofactor going from E₀ to E₁, the X-ray study implies a reduction. DFT methods have here been used, which have previously been benchmarked on a set of redox enzymes that led to the conclusion that the accuracy is about 3 kcal mol^-1 in all cases, even for redox transitions. The methodology should therefore be adequate for resolving the question of the hydride presence in E₁. As a comparison, calculations are performed on both Mo- and V-nitrogenase with the same conclusion. The conclusion from the calculations has far reaching consequences for the mechanism of nitrogenase.

Computational modeling of redox enzymes

A computational methodology is briefly described, which appears to be able to accurately describe the mechanisms of redox active enzymes. The method is built on hybrid density functional theory where the inclusion of a fraction of exact exchange is critical. Two examples of where the methodology has been applied are described. The first example is the mechanism for water oxidation in photosystem II, and the second one is the mechanism for N₂ activation by nitrogenase. The mechanism for PSII has obtained very strong support from subsequent experiments. For nitrogenase, the calculations suggest that there should be an activation process prior to catalysis, which is still strongly debated.

Quantum Mechanical Calculations of Redox Potentials of the Metal Clusters in Nitrogenase

We have calculated redox potentials of the two metal clusters in Mo-nitrogenase with quantum mechanical (QM) calculations. We employ an approach calibrated for iron-sulfur clusters with 1-4 Fe ions, involving QM-cluster calculations in continuum solvent and large QM systems (400-500 atoms), based on structures from combined QM and molecular mechanics (QM/MM) geometry optimisations. Calculations on the P-cluster show that we can reproduce the experimental redox potentials within 0.33 V. This is similar to the accuracy obtained for the smaller clusters, although two of the redox reactions involve also proton transfer. The calculated P1+/PN redox potential is nearly the same independently of whether P1+ is protonated or deprotonated, explaining why redox titrations do not show any pH dependence. For the FeMo cluster, the calculations clearly show that the formal oxidation state of the cluster in the resting E0 state is MoIIIFe3IIFe4III , in agreement with previous experimental studies and QM calculations. Moreover, the redox potentials of the first five E0-E4 states are nearly constant, as is expected if the electrons are delivered by the same site (the P-cluster). However, the redox potentials are insensitive to the formal oxidation states of the Fe ion (i.e., whether the added protons bind to sulfide or Fe ions). Finally, we show that the later (E4-E8) states of the reaction mechanism have redox potential that are more positive (i.e., more exothermic) than that of the E0/E1 couple.

QM/MM Study of Partial Dissociation of S2B for the E2 Intermediate of Nitrogenase

Nitrogenase is the only enzyme that can cleave the triple bond in N2, making nitrogen available for all lifeforms. Previous computational studies have given widely diverging results regarding the reaction mechanism of the enzyme. For example, some recent studies have suggested that one of the μ2-bridging sulfide ligands (S2B) may dissociate from one of the Fe ions when protonated in the doubly reduced and protonated E2 state, whereas other studies indicated that such half-dissociated states are unfavorable. We have examined how the relative energies of 26 structures of the E2 state depend on details of combined quantum mechanical and molecular mechanical (QM/MM) calculations. We show that the selection of the broken-symmetry state, the basis set, relativistic effects, the size of the QM system, relaxation of the surroundings, and the conformations of the bound protons may affect the relative energies of the various structures by up to 12, 22, 9, 20, 37, and 33 kJ/mol, respectively. However, they do not change the preferred type of structures. On the other hand, the choice of the DFT functional strongly affects the preferences. The hybrid B3LYP functional strongly prefers doubly protonation of the central carbide ion, but such a structure is not consistent with experimental EPR data. Other functionals suggest structures with a hydride ion, in agreement with the experiments, and show that the ion bridges between Fe2 and Fe6. Moreover, there are two structures of the same type that are degenerate within 1-5 kJ/mol, in agreement with the observation of two EPR signals. However, the pure generalized gradient approximation (GGA) functional TPSS favors structures with a protonated S2B also bridging Fe2 and Fe6, whereas r2SCAN (meta-GGA) and TPSSh (hybrid) prefer structures with S2B dissociated from Fe2 (but remaining bound to Fe6). The energy difference between the two types of structure is so small (7-18 kJ/mol) that both types need to be considered in future investigations of the mechanism of nitrogenase.

Review on the QM/MM Methodologies and Their Application to Metalloproteins

The multiscaling quantum mechanics/molecular mechanics (QM/MM) approach was introduced in 1976, while the extensive acceptance of this methodology started in the 1990s. The combination of QM/MM approach with molecular dynamics (MD) simulation, otherwise known as the QM/MM/MD approach, is a powerful and promising tool for the investigation of chemical reactions’ mechanism of complex molecular systems, drug delivery, properties of molecular devices, organic electronics, etc. In the present review, the main methodologies in the multiscaling approaches, i.e., density functional theory (DFT), semiempirical methodologies (SE), MD simulations, MM, and their new advances are discussed in short. Then, a review on calculations and reactions on metalloproteins is presented, where particular attention is given to nitrogenase that catalyzes the conversion of atmospheric nitrogen molecules N₂ into NH₃ through the process known as nitrogen fixation and the FeMo-cofactor.

A conformational role for NifW in the maturation of molybdenum nitrogenase P-cluster

Reduction of dinitrogen by molybdenum nitrogenase relies on complex metalloclusters: the [8Fe:7S] P-cluster and the [7Fe:9S:Mo:C:homocitrate] FeMo-cofactor. Although both clusters bear topological similarities and require the reductive fusion of [4Fe:4S] sub-clusters to achieve their respective assemblies, P-clusters are assembled directly on the NifD2K2 polypeptide prior to the insertion of FeMo-co, which is fully assembled separately from NifD2K2. P-cluster maturation involves the iron protein NifH2 as well as several accessory proteins, whose role has not been elucidated. In the present work, two NifD2K2 species bearing immature P-clusters were isolated from an Azotobacter vinelandii strain in which the genes encoding NifH and the accessory protein NifZ were deleted, and characterized by X-ray absorption spectroscopy and EPR. These analyses showed that both NifD2K2 complexes harbor clusters that are electronically and structurally similar, with each NifDK unit containing two [4Fe:4S]2+/+ clusters. Binding of the accessory protein NifW parallels a decrease in the distance between these clusters, as well as a subtle change in their coordination. These results support a conformational role for NifW in P-cluster biosynthesis, bringing the two [4Fe:4S] precursors closer prior to their fusion, which may be crucial in challenging cellular contexts.

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 r2SCAN 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 r2SCAN 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 r2SCAN, 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.

Enzymatic X-ray absorption spectroelectrochemistry

The ability to observe the changes that occur at an enzyme active site during electrocatalysis can provide very valuable information for understanding the mechanism and ultimately aid in catalyst design. Herein, we discuss the development of X-ray absorption spectroscopy (XAS) in combination with electrochemistry for operando studies of enzymatic systems. XAS has had a long history of enabling geometric and electronic structural insights into the catalytic active sites of enzymes, however, XAS combined with electrochemistry (XA-SEC) has been exceedingly rare in bioinorganic applications. Herein, we discuss the challenges and opportunities of applying operando XAS to enzymatic electrocatalysts. The challenges due to the low concentration of the photoabsorber and the instability of the protein in the X-ray beam are discussed. Methods for immobilizing enzymes on the electrodes, while maintaining full redox control are highlighted. A case study of combined XAS and electrochemistry applied to a [NiFe] hydrogenase is presented. By entrapping the [NiFe] hydrogenase in a redox polymer, relatively high protein concentrations can be achieved on the electrode surface, while maintaining redox control. Overall, it is demonstrated that the experiments are feasible, but require precise redox control over the majority of the absorber atoms and careful controls to discriminate between electrochemically-driven changes and beam damage. Opportunities for future applications are discussed.

Carbon Monoxide Binding to the Iron-Molybdenum Cofactor of Nitrogenase: a Detailed Quantum Mechanics/Molecular Mechanics Investigation

Carbon monoxide (CO) is a well-known inhibitor of nitrogenase activity. Under turnover conditions, CO binds to FeMoco, the active site of Mo nitrogenase. Time-resolved IR measurements suggest an initial terminal CO at 1904 cm^-1 that converts to a bridging CO at 1715 cm^-1, and an X-ray structure shows that CO can displace one of the bridging belt sulfides of FeMoco. However, the CO-binding redox state(s) of FeMoco (E_n) and the role of the protein environment in stabilizing specific CO-bound intermediates remain elusive. In this work, we carry out an in-depth analysis of the CO-FeMoco interaction based on quantum chemical calculations addressing different aspects of the electronic structure. (1) The local electronic structure of the Fe-CO bond is studied through diamagnetically substituted FeMoco. (2) A cluster model of FeMoco within a polarizable continuum illustrates how CO binding may affect the spin-coupling between the metal centers. (3) A QM/MM model incorporates the explicit influence of the amino acid residues surrounding FeMoco in the MoFe protein. The QM/MM model predicts both a terminal and a bridging CO in the E₁ redox state. The scaled calculated CO frequencies (1922 and 1716 cm^-1, respectively) are in good agreement with the experimentally observed IR bands supporting CO binding to the E₁ state. Alternatively, an E₂ state QM/MM model, which has the same atomic structure as the CO-bound X-ray structure, features a semi-bridging CO with a scaled calculated frequency (1718 cm^-1) similar to the bridging CO in the E₁ model.

The E2 state of FeMoco: Hydride Formation versus Fe Reduction and a Mechanism for H2 Evolution

The iron‐molybdenum cofactor (FeMoco) is responsible for dinitrogen reduction in Mo nitrogenase. Unlike the resting state, E₀, reduced states of FeMoco are much less well characterized. The E₂ state has been proposed to contain a hydride but direct spectroscopic evidence is still lacking. The E₂ state can, however, relax back the E₀ state via a H₂ side‐reaction, implying a hydride intermediate prior to H₂ formation. This E₂→E₀ pathway is one of the primary mechanisms for H₂ formation under low‐electron flux conditions. In this study we present an exploration of the energy surface of the E₂ state. Utilizing both cluster‐continuum and QM/MM calculations, we explore various classes of E₂ models: including terminal hydrides, bridging hydrides with a closed or open sulfide‐bridge, as well as models without. Importantly, we find the hemilability of a protonated belt‐sulfide to strongly influence the stability of hydrides. Surprisingly, non‐hydride models are found to be almost equally favorable as hydride models. While the cluster‐continuum calculations suggest multiple possibilities, QM/MM suggests only two models as contenders for the E₂ state. These models feature either i) a bridging hydride between Fe₂ and Fe₆ and an open sulfide‐bridge with terminal SH on Fe₆ (E₂‐hyd) or ii) a double belt‐sulfide protonated, reduced cofactor without a hydride (E₂‐nonhyd). We suggest both models as contenders for the E₂ redox state and further calculate a mechanism for H₂ evolution. The changes in electronic structure of FeMoco during the proposed redox‐state cycle, E₀→E₁→E₂→E₀, are discussed.

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.

Preparation and spectroscopic characterization of lyophilized Mo nitrogenase

Mo nitrogenase is the primary source of biologically fixed nitrogen, making this system highly interesting for developing new, energy efficient ways of ammonia production. Although heavily investigated, studies of the active site of this enzyme have generally been limited to spectroscopic methods that are compatible with the presence of water and relatively low protein concentrations. One method of overcoming this limitation is through lyophilization, which allows for measurements to be performed on solvent free, high concentration samples. This method also has the potential for allowing efficient protein storage and solvent exchange. To investigate the viability of this preparatory method with Mo nitrogenase, we employ a combination of electron paramagnetic resonance, Mo and Fe K-edge X-ray absorption spectroscopy, and acetylene reduction assays. Our results show that while some small distortions in the metallocofactors occur, oxidation and spin states are maintained through the lyophilization process and that reconstitution of either lyophilized protein component into buffer restores acetylene reducing activity.

Construction of Synthetic Models for Nitrogenase-Relevant NifB Biogenesis Intermediates and Iron-Carbide-Sulfide Clusters

The family of nitrogenase enzymes catalyzes the reduction of atmospheric dinitrogen (N2) to ammonia under remarkably benign conditions of temperature, pressure, and pH. Therefore, the development of synthetic complexes or materials that can similarly perform this reaction is of critical interest. The primary obstacle for obtaining realistic synthetic models of the active site iron-sulfur-carbide cluster (e.g., FeMoco) is the incorporation of a truly inorganic carbide. This review summarizes the present state of knowledge regarding biological and chemical (synthetic) incorporation of carbide into iron-sulfur clusters. This includes the Nif cluster of proteins and associated biochemistry involved in the endogenous biogenesis of FeMoco. We focus on the chemical (synthetic) incorporation portion of our own efforts to incorporate and modify C1 units in iron/sulfur clusters. We also highlight recent contributions from other research groups in the area toward C1 and/or inorganic carbide insertion.

Mixed-Valent Diiron μ-Carbyne, μ-Hydride Complexes: Implications for Nitrogenase

Binding of N₂ by the FeMo-cofactor of nitrogenase is believed to occur after transfer of 4 e^- and 4 H⁺ equivalents to the active site. Although pulse EPR studies indicate the presence of two Fe-(μ-H)-Fe moieties, the structural and electronic features of this mixed valent intermediate remain poorly understood. Toward an improved understanding of this bioorganometallic cluster, we report herein that diiron μ-carbyne complex (P₆ArC)Fe₂(μ-H) can be oxidized and reduced, allowing for the first time spectral characterization of two EPR-active Fe(μ-C)(μ-H)Fe model complexes linked by a 2 e^- transfer which bear some resemblance to a pair of E_n and E_n+2 states of nitrogenase. Both species populate S = 1/2 states at low temperatures, and the influence of valence (de)localization on the spectroscopic signature of the μ-hydride ligand was evaluated by pulse EPR studies. Compared to analogous data for the {Fe₂(μ-H)}₂ state of FeMoco (E₄(4H)), the data and analysis presented herein suggest that the hydride ligands in E₄(4H) bridge isovalent (most probably Fe^III) metal centers. Although electron transfer involves metal-localized orbitals, investigations of [(P₆ArC)Fe₂(μ-H)]⁺¹ and [(P₆ArC)Fe₂(μ-H)]^-1 by pulse EPR revealed that redox chemistry induces significant changes in Fe-C covalency (-50% upon 2 e^- reduction), a conclusion further supported by X-ray absorption spectroscopy, ⁵⁷Fe Mössbauer studies, and DFT calculations. Combined, our studies demonstrate that changes in covalency buffer against the accumulation of excess charge density on the metals by partially redistributing it to the bridging carbon, thereby facilitating multielectron transformations.

Quantum Mechanics/Molecular Mechanics Study of Resting-State Vanadium Nitrogenase: Molecular and Electronic Structure of the Iron-Vanadium Cofactor

The nitrogenase enzymes are responsible for all biological nitrogen reduction. How this is accomplished at the atomic level, however, has still not been established. The molybdenum-dependent nitrogenase has been extensively studied and is the most active catalyst for dinitrogen reduction of the nitrogenase enzymes. The vanadium-dependent form, on the other hand, displays different reactivity, being capable of CO and CO2 reduction to hydrocarbons. Only recently did a crystal structure of the VFe protein of vanadium nitrogenase become available, paving the way for detailed theoretical studies of the iron-vanadium cofactor (FeVco) within the protein matrix. The crystal structure revealed a bridging 4-atom ligand between two Fe atoms, proposed to be either a CO32- or NO3- ligand. Using a quantum mechanics/molecular mechanics model of the VFe protein, starting from the 1.35 Å crystal structure, we have systematically explored multiple computational models for FeVco, considering either a CO32- or NO3- ligand, three different redox states, and multiple broken-symmetry states. We find that only a [VFe7S8C(CO3)]2- model for FeVco reproduces the crystal structure of FeVco well, as seen in a comparison of the Fe-Fe and V-Fe distances in the computed models. Furthermore, a broken-symmetry solution with Fe2, Fe3, and Fe5 spin-down (BS7-235) is energetically preferred. The electronic structure of the [VFe7S8C(CO3)]2- BS7-235 model is compared to our [MoFe7S9C]- BS7-235 model of FeMoco via localized orbital analysis and is discussed in terms of local oxidation states and different degrees of delocalization. As previously found from Fe X-ray absorption spectroscopy studies, the Fe part of FeVco is reduced compared to FeMoco, and the calculations reveal Fe5 as locally ferrous. This suggests resting-state FeVco to be analogous to an unprotonated E1 state of FeMoco. Furthermore, V-Fe interactions in FeVco are not as strong compared to Mo-Fe interactions in FeMoco. These clear differences in the electronic structures of otherwise similar cofactors suggest an explanation for distinct differences in reactivity.

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.

Reduction of Substrates by Nitrogenases

Nitrogenase is the enzyme that catalyzes biological N₂ reduction to NH₃. This enzyme achieves an impressive rate enhancement over the uncatalyzed reaction. Given the high demand for N₂ fixation to support food and chemical production and the heavy reliance of the industrial Haber-Bosch nitrogen fixation reaction on fossil fuels, there is a strong need to elucidate how nitrogenase achieves this difficult reaction under benign conditions as a means of informing the design of next generation synthetic catalysts. This Review summarizes recent progress in addressing how nitrogenase catalyzes the reduction of an array of substrates. New insights into the mechanism of N₂ and proton reduction are first considered. This is followed by a summary of recent gains in understanding the reduction of a number of other nitrogenous compounds not considered to be physiological substrates. Progress in understanding the reduction of a wide range of C-based substrates, including CO and CO₂, is also discussed, and remaining challenges in understanding nitrogenase substrate reduction are considered.

A model for dinitrogen binding in the E4 state of nitrogenase

Molybdenum nitrogenase is one of the most intriguing metalloenzymes in nature, featuring an exotic iron-molybdenum-sulfur cofactor, FeMoco, whose mode of action remains elusive. In particular, the molecular and electronic structure of the N2-binding E4 state is not known. In this study we present theoretical QM/MM calculations of new structural models of the E4 state of molybdenum-dependent nitrogenase and compare to previously suggested models for this enigmatic redox state. We propose two models as possible candidates for the E4 state. Both models feature two hydrides on the FeMo cofactor, bridging atoms Fe2 and Fe6 with a terminal sulfhydryl group on either Fe2 or Fe6 (derived from the S2B bridge) and the change in coordination results in local lower-spin electronic structure at Fe2 and Fe6. These structures appear consistent with the bridging hydride proposal put forward from ENDOR studies and are calculated to be lower in energy than other proposed models for E4 at the TPSSh-QM/MM level of theory. We critically analyze the DFT method dependency in calculations of FeMoco that has resulted in strikingly different proposals for this state. Importantly, dinitrogen binds exothermically to either Fe2 or Fe6 in our models, contrary to others, an effect rationalized via the unique ligand field (from the hydrides) at the Fe with an empty coordination site. A low-spin Fe site is proposed as being important to N2 binding. Furthermore, the geometries of these states suggest a feasible reductive elimination step that could follow, as experiments indicate. Via this step, two electrons are released, reducing the cofactor to yield a distorted 4-coordinate Fe2 or Fe6 that partially activates N2. We speculate that stabilization of an N2-bound Fe(i) at Fe6 (not found for Fe2 model) via reductive elimination is a crucial part of N2 activation in nitrogenases, possibly aided by the apical heterometal ion (Mo or V). By using protons from the sulfhydryl group (to regenerate the sulfide bridge between Fe2 and Fe6) and the nearby homocitrate hydroxy group, we calculate a plausible route to yield a diazene intermediate. This is found to be more favorable with the Fe6-bound model than the Fe2-bound model; however, this protonation is uphill in energy, suggesting protonation of N2 might occur later in the catalytic cycle or via another mechanism.