Structures and reaction dynamics of N2 and H2 binding at FeMo-co, the active site of nitrogenase

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).

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.

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.

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.

Calculating the chemical mechanism of nitrogenase: new working hypotheses

The enzyme nitrogenase converts N₂ to NH₃ with stoichiometry N₂ + 8H⁺ + 8e^- → 2NH₃ + H₂. The mechanism is chemically complex with multiple steps that must be consistent with much accumulated experimental information, including exchange of H₂ and N₂ and the N₂-dependent hydrogenation of D₂ to HD. Previous investigations have developed a collection of working hypotheses that guide ongoing density functional investigations of mechanistic steps and sequences. These include (i) hypotheses about the serial provision of protons and their conversion to H atoms bonded to S and Fe atoms of the FeMo-co catalytic site, (ii) the migration of H atoms over the surface of FeMo-co, (iii) the roles of His195, (iv) identification of three protein channels, one for the ingress of N₂, a separate pathway for the passage of exogenous H₂ (D₂) and product H₂ (HD), and a hydrophilic pathway for egress of product NH₃. Two additional working hypotheses are described in this paper. N₂ passing along the N₂ channel approaches and binds end-on to the exo coordination position of Fe2, with favourable energetics when FeMo-co is pre-hydrogenated. This exo-Fe2-N₂ is apparently not reduced but has a promotional role by expanding the reaction zone. A second N₂ can enter via the N₂ ingress channel and bind at the endo-Fe6 position, where it is surrounded by H atom donors suitable for the N₂ → NH₃ conversion. It is proposed that this endo-Fe6 position is also the binding site for H₂ (generated or exogenous), accounting for the competitive inhibition of N₂ reduction by H₂. The HD reaction occurs at the endo-Fe6 site, promoted by N₂ at the exo-Fe2 site. The second hypothesis concerns the most stable electronic states of FeMo-co with ligands bound at Fe2 and Fe6, and provides a protocol for management of electronic states in mechanism calculations.