Survey of the Geometric and Electronic Structures of the Key Hydrogenated Forms of FeMo-co, the Active Site of the Enzyme Nitrogenase: Principles of the Mechanistically Significant Coordination Chemistry

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.

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.

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.

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.

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.

Biological nitrogen fixation in theory, practice, and reality: a perspective on the molybdenum nitrogenase system

Nitrogenase is the sole enzyme responsible for the ATP-dependent conversion of atmospheric dinitrogen into the bioavailable form of ammonia (NH₃ ), making this protein essential for the maintenance of the nitrogen cycle and thus life itself. Despite the widespread use of the Haber-Bosch process to industrially produce NH₃ , biological nitrogen fixation still accounts for half of the bioavailable nitrogen on Earth. An important feature of nitrogenase is that it operates under physiological conditions, where the equilibrium strongly favours ammonia production. This biological, multielectron reduction is a complex catalytic reaction that has perplexed scientists for decades. In this review, we explore the current understanding of the molybdenum nitrogenase system based on experimental and computational research, as well as the limitations of the crystallographic, spectroscopic, and computational techniques employed. Finally, essential outstanding questions regarding the nitrogenase system will be highlighted alongside suggestions for future experimental and computational work to elucidate this essential yet elusive process.

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.

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

The chemical reactions occurring at the Fe₇MoS₉C(homocitrate) cluster, FeMo-co, the active site of the enzyme nitrogenase (N₂ → NH₃), are enigmatic. Experimental information collected over a long period reveals aspects of the roles of N₂ and H₂, each with more than one type of reactivity. This paper reports investigations of the binding of H₂ and N₂ at intact FeMo-co, using density functional simulations of a large 486 atom relevant portion of the protein, resulting in 27 new structures containing H₂ and/or N₂ bound at the exo and endo coordination sites of the participating Fe atoms, Fe2 and Fe6. Binding energies and transition states for association/dissociation are determined, and trajectories for the approach, binding and separation of H₂/N₂ are described, including diffusion of these small molecules through proximal protein. Influences of surrounding amino acids are identified. FeMo-co deforms geometrically when binding H₂ or N₂, and a procedure for calculating the energy cost involved, the adaptation energy, is introduced here. Adaptation energies, which range from 7 to 36 kcal mol^-1 for the reported structures, are influenced by the protonation state of the His195 side chain. Seven N₂ structures and three H₂ structures have negative binding free energies, which include the estimated entropy penalties for binding of N₂, H₂ from proximal protein. These favoured structures have N₂ bound end-on at exo-Fe2, exo-Fe6 and endo-Fe2 positions of FeMo-co, and H₂ bound at the endo-Fe2 position. Various postulated structures with N₂ bridging Fe2 and Fe6 revert to end-on-N₂ at endo positions. The structures are also assessed via the calculated potential energy barriers for association and dissociation. Barriers to the binding of H₂ range from 1 to 20 kcal mol^-1 and barriers to dissociation of H₂ range from 3 to 18 kcal mol^-1. Barriers to the binding of N₂, in either side-on or end-on mode, range from 2 to 18 kcal mol^-1, while dissociation of bound N₂ encounters barriers of 3 to 8 kcal mol^-1 for side-on bonding and 7 to 18 kcal mol^-1 for end-on bonding. These results allow formulation of mechanisms for the H₂/N₂ exchange reaction, and three feasible mechanisms for associative exchange and three for dissociative exchange are identified. Consistent electronic structures and potential energy surfaces are maintained throughout. Changes in the spin populations of Fe2 and Fe6 connected with cluster deformation and with metal-ligand bond formation are identified.

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.

On the Use of Normalized Metrics for Density Sensitivity Analysis in DFT

In the past years, there has been a discussion about how the errors in density functional theory might be related to errors in the self-consistent densities obtained from different density functional approximations. This, in turn, brings up the discussion about the different ways in which we can measure such errors and develop metrics that assess the sensitivity of calculated energies to changes in the density. It is important to realize that there cannot be a unique metric in order to look at this density sensitivity, simultaneously needing size-extensive and size-intensive metrics. In this study, we report two metrics that are widely applicable to any density functional approximation. We also show how they can be used to classify different chemical systems of interest with respect to their sensitivity to small variations in the density.

Quantum-refinement studies of the bidentate ligand of V‑nitrogenase and the protonation state of CO-inhibited Mo‑nitrogenase

Nitrogenase is the only enzyme that can cleave the triple bond in N₂, making nitrogen available to plants (although the enzyme itself is strictly microbial). It has been studied extensively with both experimental and computational methods, but many details of the reaction mechanism are still unclear. X-ray crystallography is the main source of structural information for biomacromolecules, but it has problems to discern hydrogen atoms or to distinguish between elements with the same number of electrons. These problems can sometimes be alleviated by introducing quantum chemical calculations in the refinement, providing information about the ideal structure (in the same way as the empirical restraints used in standard crystallographic refinement) and comparing different interpretations of the structure with normal crystallographic and quantum mechanical quality measures. We have performed such quantum-refinement calculations to address two important issues for nitrogenase. First, we show that the bidentate ligand of the active-site FeV cluster in V‑nitrogenase is carbonate, rather than bicarbonate or nitrate. Second, we study the CO-inhibited structure of Mo‑nitrogenase. CO binds to a reduced and protonated state of the enzyme by replacing one of the sulfide ions (S2B) in the active-site FeMo cluster. We examined if it is possible to deduce from the crystal structure the location of the protons. Our results indicates that the crystal structure is best modelled as fully deprotonated.

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.

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.

N2H2 binding to the nitrogenase FeMo cluster studied by QM/MM methods

We have made a systematic combined quantum mechanical and molecular mechanical (QM/MM) investigation of possible structures of the N₂ bound state of nitrogenase. We assume that N₂ is immediately protonated to a N₂H₂ state, thereby avoiding the problem of determining the position of the protons in the cluster. We have systematically studied both end-on and side-on structures, as well as both HNNH and NNH₂ states. Our results indicate that the binding of N₂H₂ is determined more by interactions and steric clashes with the surrounding protein than by the intrinsic preferences of the ligand and the cluster. The best binding mode with both the TPSS and B3LYP density-functional theory methods has trans-HNNH terminally bound to Fe2. It is stabilised by stacking of the substrate with His-195 and Ser-278. However, several other structures come rather close in energy (within 3-35 kJ/mol) at least in some calculations: The corresponding cis-HNNH structure terminally bound to Fe2 is second best with B3LYP. A structure with HNNH₂ terminally bound to Fe6 is second most stable with TPSS (where the third proton is transferred to the substrate from the homocitrate ligand). Structures with trans-HNNH, bound to Fe4 or Fe6, or cis-HNNH bound to Fe6 are also rather stable. Finally, with the TPSS functional, a structure with cis-HNNH side-on binding to the Fe3-Fe4-Fe5-Fe7 face of the cluster is also rather low in energy, but all side-on structures are strongly disfavoured by the B3LYP method.

What Is the Structure of the E4 Intermediate in Nitrogenase?

Nitrogenase is the only enzyme that can cleave the strong triple bond in N₂. The active site contains a complicated MoFe₇S₉C cluster. It is believed that it needs to accept four protons and electrons, forming the E₄ state, before it can bind N₂. However, there is no consensus on the atomic structure of the E₄ state. Experimental studies indicate that it should contain two hydride ions bridging two pairs of Fe ions, and it has been suggested that both hydride ions as well as the two protons bind on the same face of the cluster. On the other hand, density functional theory (DFT) studies have indicated that it is energetically more favorable with either three hydride ions or with a triply protonated carbide ion, depending on the DFT functional. We have performed a systematic combined quantum mechanical and molecular mechanical (QM/MM) study of possible E₄ states with two bridging hydride ions. Our calculations suggest that the most favorable structure has hydride ions bridging the Fe2/6 and Fe3/7 ion pairs. In fact, such structures are 14 kJ/mol more stable than structures with three hydride ions, showing that pure DFT functionals give energetically most favorable structures in agreement with experiments. An important reason for this finding is that we have identified a new type of broken-symmetry state that involves only two Fe ions with minority spin, in contrast to the previously studied states with three Fe ions with minority spin. The energetically best structures have the two hydride ions on different faces of the FeMo cluster, whereas better agreement with ENDOR data is obtained if they are on the same face; such structures are only 6–22 kJ/mol less stable.

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

Biological nitrogen fixation is predominately accomplished through Mo nitrogenase, which utilizes a complex MoFe7S9C catalytic cluster to reduce N2 to NH3. This cluster requires the accumulation of three to four reducing equivalents prior to binding N2; however, despite decades of research, the intermediate states formed prior to N2 binding are still poorly understood. Herein, we use Mo and Fe K-edge X-ray absorption spectroscopy and QM/MM calculations to investigate the nature of the E1 state, which is formed following the addition of the first reducing equivalent to Mo nitrogenase. By analyzing the extended X-ray absorption fine structure (EXAFS) region, we provide structural insight into the changes that occur in the metal clusters of the protein when forming the E1 state, and use these metrics to assess a variety of possible models of the E1 state. The combination of our experimental and theoretical results supports that formation of E1 involves an Fe-centered reduction combined with the protonation of a belt-sulfide of the cluster. Hence, these results provide critical experiment and computational insight into the mechanism of this important enzyme.

The mechanism for nitrogenase including all steps

Nitrogen in the air is turned into biologically useful ammonia by the nitrogenase enzyme. The leading member of this group has a cofactor with one molybdenum and seven irons linked together by sulfurs. The structure that binds N₂ has a triply protonated carbide and a rotated homocitrate. Both these structural changes are necessary for the activation.