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

Putative reaction mechanism of nitrogenase with a half-dissociated S2B ligand

We have studied whether dissociation of the S2B sulfide ligand from one of its two coordinating Fe ions may affect the later parts of the reaction mechanism of nitrogenase. Such dissociation has been shown to be favourable for the E2-E4 states in the reaction mechanism, but previous studies have assumed that S2B either remains bridging or has fully dissociated from the active-site FeMo cluster. We employ combined quantum mechanical and molecular mechanical (QM/MM) calculations with two density-functional theory methods, r2SCAN and TPSSh. To make dissociation of S2B possible, we have added a proton to this group throughout the reaction. We study the reaction starting from the E4 state with N2H2 bound to the cluster. Our results indicate that half-dissociation of S2B is unfavourable in most steps of the reaction mechanism. We observe favourable half-dissociation of S2B only when NH or NH2 is bound to the cluster, bridging Fe2 and Fe6. However, the former state is most likely not involved in the reaction mechanism and the latter state is only an intermittent intermediate of the E7 state. Therefore, half-dissociation of S2B seems to play only a minor role in the later parts of the reaction mechanism of nitrogenase. Our suggested mechanism with a protonated S2B is alternating (the two N atoms of the substrate is protonated in an alternating manner) and the substrate prefers to bind to Fe2, in contrast to the preferred binding to Fe6 observed when S2B is unprotonated and bridging Fe2 and Fe6.

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.

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.

Mechanism of Nitrogen Reduction to Ammonia in a Diiron Model of Nitrogenase

Nitrogenase is a fascinating enzyme in biology that reduces dinitrogen from air to ammonia through stepwise reduction and protonation. Despite it being studied in detail by experimental and computational groups, there are still many unknown factors in the catalytic cycle of nitrogenase, especially related to the addition of protons and electrons and their order. A recent biomimetic study characterized a potential dinitrogen-bridged diiron cluster as a synthetic model of nitrogenase. Using strong acid and reductants, the dinitrogen was converted into ammonia molecules, but details of the mechanism remains unknown. In particular, it was unclear from the experimental studies whether the proton and electron transfer steps are sequential or alternating. Moreover, the work failed to establish what the function of the diiron core is and whether it split into mononuclear iron fragments during the reaction. To understand the structure and reactivity of the biomimetic dinitrogen-bridged diiron complex [(P2P'PhFeH)2(μ-N2)] with triphenylphosphine ligands, we performed a density functional theory study. Our computational methods were validated against experimental crystal structure coordinates, Mössbauer parameters, and vibrational frequencies and show excellent agreement. Subsequently, we investigated the alternating and consecutive addition of electrons and protons to the system. The calculations identify a number of possible reaction channels, namely, same-site protonation, alternating protonation, and complex dissociation into mononuclear iron centers. The calculations show that the overall mechanism is not a pure sequential set of electron and proton transfers but a mixture of alternating and consecutive steps. In particular, the first reaction steps will start with double proton transfer followed by an electron transfer, while thereafter, there is another proton transfer and a second electron transfer to give a complex whereby ammonia can split off with a low energetic barrier. The second channel starts with alternating protonation of the two nitrogen atoms, whereafter the initial double proton transfer, electrons and protons are added sequentially to form a hydrazine-bound complex. The latter split off ammonia spontaneously after further protonation. The various reaction channels are analyzed with valence bond and orbital diagrams. We anticipate the nitrogenase enzyme to operate with mixed alternating and consecutive protonation and electron transfer steps.

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.

Dynamic effects on ligand field from rapid hydride motion in an iron(ii) dimer with an S = 3 ground state

Hydride complexes are important in catalysis and in iron-sulfur enzymes like nitrogenase, but the impact of hydride mobility on local iron spin states has been underexplored. We describe studies of a dimeric diiron(ii) hydride complex using X-ray and neutron crystallography, Mössbauer spectroscopy, magnetism, DFT, and ab initio calculations, which give insight into the dynamics and the electronic structure brought about by the hydrides. The two iron sites in the dimer have differing square-planar (intermediate-spin) and tetrahedral (high-spin) iron geometries, which are distinguished only by the hydride positions. These are strongly coupled to give an S total = 3 ground state with substantial magnetic anisotropy, and the merits of both localized and delocalized spin models are discussed. The dynamic nature of the sites is dependent on crystal packing, as shown by changes during a phase transformation that occurs near 160 K. The change in dynamics of the hydride motion leads to insight into its influence on the electronic structure. The accumulated data indicate that the two sites can trade geometries by rotating the hydrides, at a rate that is rapid above the phase transition temperature but slow below it. This small movement of the hydrides causes large changes in the ligand field because they are strong-field ligands. This suggests that hydrides could be useful in catalysis not only due to their reactivity, but also due to their ability to rapidly modulate the local electronic structure and spin states at metal sites.

The binding of reducible N2 in the reaction domain of nitrogenase

The binding of N₂ to FeMo-co, the catalytic site of the enzyme nitrogenase, is central to the conversion to NH₃, but also has a separate role in promoting the N₂-dependent HD reaction (D₂ + 2H⁺ + 2e^- → 2HD). The protein surrounding FeMo-co contains a clear channel for ingress of N₂, directly towards the exo-coordination position of Fe2, a position which is outside the catalytic reaction domain. This led to the hypothesis [I. Dance, Dalton Trans., 2022, 51, 12717] of 'promotional' N₂ bound at exo-Fe2, and a second 'reducible' N₂ bound in the reaction domain, specifically the endo-coordination position of Fe2 or Fe6. The range of possibilities for the binding of reducible N₂ in the presence of bound promotional N₂ is described here, using density functional simulations with a 486 atom model of the active site and surrounding protein. The pathway for ingress of the second N₂ through protein, past the first N₂ at exo-Fe2, and tumbling into the binding domain between Fe2 and Fe6, is described. The calculations explore 24 structures involving 6 different forms of hydrogenated FeMo-co, including structures with S2BH unhooked from Fe2 but tethered to Fe6. The calculations use the most probable electronic states. End-on (η¹) binding of N₂ at the endo position of either Fe2 or Fe6 is almost invariably exothermic, with binding potential energies ranging up to -18 kcal mol^-1. Many structures have binding energies in the range -6 to -14 kcal mol^-1. The relevant entropic penalty for N₂ binding from a diffusible position within the protein is estimated to be 4 kcal mol^-1, and so the binding free energies for reducible N₂ are suitably negative. N₂ binding at endo-Fe2 is stronger than at endo-Fe6 in three of the six structure categories. In many cases the reaction domain containing reducible N₂ is expanded. These results inform computational simulation of the subsequent steps in which surrounding H atoms transfer to reducible N₂.

The HD Reaction of Nitrogenase: a Detailed Mechanism

Nitrogenase is the enzyme that converts N₂ to NH₃ under ambient conditions. The chemical mechanism of this catalysis at the active site FeMo-co [Fe₇ S₉ CMo(homocitrate)] is unknown. An obligatory co-product is H₂ , while exogenous H₂ is a competitive inhibitor. Isotopic substitution using exogenous D₂ revealed the N₂ -dependent reaction D₂ +2H⁺ +2e^- →2HD (the 'HD reaction'), together with a collection of additional experimental characteristics and requirements. This paper describes a detailed mechanism for the HD reaction, developed and elaborated using density functional simulations with a 486-atom model of the active site and surrounding protein. First D₂ binds at one Fe atom (endo-Fe6 coordination position), where it is flanked by H-Fe6 (exo position) and H-Fe2 (endo position). Then there is synchronous transfer of these two H atoms to bound D₂ , forming one HD bound to Fe2 and a second HD bound to Fe6. These two HD dissociate sequentially. The final phase is recovery of the two flanking H atoms. These H atoms are generated, sequentially, by translocation of a proton from the protein surface to S3B of FeMo-co and combination with introduced electrons. The first H atom migrates from S3B to exo-Fe6 and the second from S3B to endo-Fe2. Reaction energies and kinetic barriers are reported for all steps. This mechanism accounts for the experimental data: (a) stoichiometry; (b) the N₂ -dependence results from promotional N₂ bound at exo-Fe2; (c) different N₂ binding K_m for the HD reaction and the NH₃ formation reaction results from involvement of two different sites; (d) inhibition by CO; (e) the non-occurrence of 2HD→H₂ +D₂ results from the synchronicity of the two transfers of H to D₂ ; (f) inhibition of HD production at high pN₂ is by competitive binding of N₂ at endo-Fe6; (g) the non-leakage of D to solvent follows from the hydrophobic environment and irreversibility of proton introduction.

13C ENDOR Characterization of the Central Carbon within the Nitrogenase Catalytic Cofactor Indicates That the CFe6 Core Is a Stabilizing "Heart of Steel"

Substrates and inhibitors of Mo-dependent nitrogenase bind and react at Fe ions of the active-site FeMo-cofactor [7Fe-9S-C-Mo-homocitrate] contained within the MoFe protein α-subunit. The cofactor contains a CFe6 core, a carbon centered within a trigonal prism of six Fe, whose role in catalysis is unknown. Targeted 13C labeling of the carbon enables electron-nuclear double resonance (ENDOR) spectroscopy to sensitively monitor the electronic properties of the Fe-C bonds and the spin-coupling scheme adopted by the FeMo-cofactor metal ions. This report compares 13CFe6 ENDOR measurements for (i) the wild-type protein resting state (E0; α-Val70) to those of (ii) α-Ile70, (iii) α-Ala70-substituted proteins; (iv) crystallographically characterized CO-inhibited "hi-CO" state; (v) E4(4H) Janus intermediate, activated for N2 binding/reduction by accumulation of 4[e-/H+]; (vi) E4(2H)* state containing a doubly reduced FeMo-cofactor without Fe-bound substrates; and (vii) propargyl alcohol reduction intermediate having allyl alcohol bound as a ferracycle to FeMo-cofactor Fe6. All states examined, both S = 1/2 and 3/2 exhibited near-zero 13C isotropic hyperfine coupling constants, Ca = [-1.3 ↔ +2.7] MHz. Density functional theory computations and natural bond orbital analysis of the Fe-C bonds show that this occurs because a (3 spin-up/3 spin-down) spin-exchange configuration of CFe6 Fe-ion spins produces cancellation of large spin-transfers to carbon in each Fe-C bond. Previous X-ray diffraction and DFT both indicate that trigonal-prismatic geometry around carbon is maintained with high precision in all these states. The persistent structure and Fe-C bonding of the CFe6 core indicate that it does not provide a functionally dynamic (hemilabile) "beating heart"─instead it acts as "a heart of steel", stabilizing the structure of the FeMo-cofactor-active site during nitrogenase catalysis.

Understanding the tethered unhooking and rehooking of S2B in the reaction domain of FeMo-co, the active site of nitrogenase

The active site of the nitrogen fixing enzyme nitrogenase is an Fe₇MoS₉C cluster, and investigations of the enigmatic chemical mechanism of the enzyme have focussed on a pair of Fe atoms, Fe2 and Fe6, and the S2B atom that bridges them. There are three proposals for the status of the Fe2-S2B-Fe6 bridge during the catalytic cycle: one that it remains intact, another that it is completely labile and absent during catalysis, and a third that S2B is hemilabile, unhooking one of its bonds to Fe2 or Fe6. This report examines the tethered unhooking of S2B and factors that affect it, using DFT calculations of 50 geometric/electronic possibilities with a 485 atom model including all relevant parts of surrounding protein. The outcomes are: (a) unhooking the S2B-Fe2 bond is feasible and favourable, but alternative unhooking of the S2B-Fe6 bond is unlikely for steric reasons, (b) energy differences between hooked and unhooked isomers are generally <10 kcal mol^-1, usually with unhooked more stable, (c) ligation at the exo-Fe6 position inhibits unhooking, (d) unhooking of hydrogenated S2B is more favourable than that of bare S2B, (e) hydrogen bonding from the NεH function of His195 to S2B occurs in hooked and unhooked forms, and possibly stabilises unhooking, (f) unhooking is reversible with kinetic barriers ranging 10-13 kcal mol^-1. The conclusion is that energetically accessible reversible unhooking of S2B or S2BH, as an intrinsic property of FeMo-co, needs to be considered in the formulation of mechanisms for the reactions of nitrogenase.

The One-Electron Reduced Active-Site FeFe-Cofactor of Fe-Nitrogenase Contains a Hydride Bound to a Formally Oxidized Metal-Ion Core

The nitrogenase active-site cofactor must accumulate 4e^-/4H⁺ (E₄(4H) state) before N₂ can bind and be reduced. Earlier studies demonstrated that this E₄(4H) state stores the reducing-equivalents as two hydrides, with the cofactor metal-ion core formally at its resting-state redox level. This led to the understanding that N₂ binding is mechanistically coupled to reductive-elimination of the two hydrides that produce H₂. The state having acquired 2e^-/2H⁺ (E₂(2H)) correspondingly contains one hydride with a resting-state core redox level. How the cofactor accommodates addition of the first e^-/H⁺ (E₁(H) state) is unknown. The Fe-nitrogenase FeFe-cofactor was used to address this question because it is EPR-active in the E₁(H) state, unlike the FeMo-cofactor of Mo-nitrogenase, thus allowing characterization by EPR spectroscopy. The freeze-trapped E₁(H) state of Fe-nitrogenase shows an S = 1/2 EPR spectrum with g = [1.965, 1.928, 1.779]. This state is photoactive, and under 12 K cryogenic intracavity, 450 nm photolysis converts to a new and likewise photoactive S = 1/2 state (denoted E₁(H)*) with g = [2.009, 1.950, 1.860], which results in a photostationary state, with E₁(H)* relaxing to E₁(H) at temperatures above 145 K. An H/D kinetic isotope effect of 2.4 accompanies the 12 K E₁(H)/E₁(H)* photointerconversion. These observations indicate that the addition of the first e^-/H⁺ to the FeFe-cofactor of Fe-nitrogenase produces an Fe-bound hydride, not a sulfur-bound proton. As a result, the cluster metal-ion core is formally one-electron oxidized relative to the resting state. It is proposed that this behavior applies to all three nitrogenase isozymes.

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.