E4 State Research Articles

Final E5 to E8 Steps in the Nitrogenase Mechanism for Nitrogen Fixation.

Nitrogenase converts nitrogen in the air to ammonia. It is often regarded as the second most important enzyme in nature after photosystem II. The mechanism for how nitrogenase is able to perform the difficult task of cleaving the strong bond in N2 is debated. It is known that for every electron that is donated to N2, two ATP are hydrolyzed. In the experimentally suggested mechanism, the activation occurs after four reductions of the ground state, but there is no suggestion for how the enzyme uses the hydrolysis energy to perform catalysis. In the theoretical mechanism, it is suggested that hydrolysis is used to reduce the electron donor. In previous papers, the steps leading to the activation of N2 in the so-called E4 state has been investigated, using both the experimental and theoretical mechanism, showing that only the theoretical one leads to agreement with EPR observations for E4. In the present paper, the four steps following E4, leading to the release of two ammonia molecules, are described using the same methodology as used in the previous studies.

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.

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.

Computational Model Study of the Experimentally Suggested Mechanism for Nitrogenase.

The mechanism for N2 activation in the E4 state of nitrogenase was investigated by model calculations. In the experimentally suggested mechanism, the E4 state is obtained after four reductions to the ground state. In a recent theoretical study, results for a different mechanism have been found in excellent agreement with available Electron Paramagnetic Resonance (EPR) experiments for E4. The two hydrides in E4 leave as H2 concertedly with the binding of N2. The mechanism suggested differs from the experimentally suggested one by a requirement for four activation steps prior to catalysis. In the present study, the experimentally suggested mechanism is studied using the same methods as those used in the previous study on the theoretical mechanism. The computed results make it very unlikely that a structure obtained after four reductions from the ground state has two hydrides, and the experimentally suggested mechanism does therefore not agree with the EPR experiments for E4. Another structure with only one hydride is here suggested to be the one that has been observed to bind N2 after only four reductions of the ground state.

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

Cryo-annealing of Photoreduced CdS Quantum Dot-Nitrogenase MoFe Protein Complexes Reveals the Kinetic Stability of the E4(2N2H) Intermediate.

A critical step in the mechanism of N2 reduction to 2NH3 catalyzed by the enzyme nitrogenase is the reaction of the four-electron/four-proton reduced intermediate state of the active-site FeMo-cofactor (E4(4H)). This state is a junction in the catalytic mechanism, either relaxing by the reaction of a metal bound Fe-hydride with a proton forming H2 or going forward with N2 binding coupled to the reductive elimination (re) of two Fe-hydrides as H2 to form the E4(2N2H) state. E4(2N2H) can relax to E4(4H) by the oxidative addition (oa) of H2 and release of N2 or can be further reduced in a series of catalytic steps to release 2NH3. If the H2 re/oa mechanism is correct, it requires that oa of H2 be associative with E4(2N2H). In this report, we have taken advantage of CdS quantum dots in complex with MoFe protein to achieve photodriven electron delivery in the frozen state, with cryo-annealing in the dark, to reveal details of the E-state species and to test the stability of E4(2N2H). Illumination of frozen CdS:MoFe protein complexes led to formation of a population of reduced intermediates. Electron paramagnetic resonance spectroscopy identified E-state signals including E2 and E4(2N2H), as well as signals suggesting the formation of E6 or E8. It is shown that in the frozen state when pN2 is much greater than pH2, the E4(2N2H) state is kinetically stable, with very limited forward or reverse reaction rates. These results establish that the oa of H2 to the E4(2N2H) state follows an associative reaction mechanism.

Understanding the Electronic Structure Basis for N2 Binding to FeMoco: A Systematic Quantum Mechanics/Molecular Mechanics Investigation.

The FeMo cofactor (FeMoco) of Mo nitrogenase is responsible for reducing dinitrogen to ammonia, but it requires the addition of 3-4 e-/H+ pairs before N2 even binds. A binding site at the Fe2/Fe3/Fe6/Fe7 face of the cofactor has long been suggested based on mutation studies, with Fe2 or Fe6 nowadays being primarily discussed as possibilities. However, the nature of N2 binding to the cofactor is enigmatic as the metal ions are coordinatively saturated in the resting state with no obvious binding site. Furthermore, the cofactor consists of high-spin Fe(II)/Fe(III) ions (antiferromagnetically coupled but also mixed-valence delocalized), which are not known to bind N2. This suggests that an Fe binding site with a different molecular and electronic structure than the resting state must be responsible for the experimentally known N2 binding in the E4 state of FeMoco. We have systematically studied N2 binding to Fe2 and Fe6 sites of FeMoco at the broken-symmetry QM/MM level as a function of the redox-, spin-, and protonation state of the cofactor. The local and global electronic structure changes to the cofactor taking place during redox events, protonation, Fe-S cleavage, hydride formation, and N2 coordination are systematically analyzed. Localized orbital and quasi-restricted orbital analysis via diamagnetic substitution is used to get insights into the local single Fe ion electronic structure in various states of FeMoco. A few factors emerge as essential to N2 binding in the calculations: spin-pairing of dxz and dyz orbitals of the N2-binding Fe ion, a coordination change at the N2-binding Fe ion aided by a hemilabile protonated sulfur, and finally hydride ligation. The results show that N2 binding to E0, E1, and E2 models is generally unfavorable, likely due to the high-energy cost of stabilizing the necessary spin-paired electronic structure of the N2-binding Fe ion in a ligand environment that clearly favors high-spin states. The results for models of E4, however, suggest a feasible model for why N2 binding occurs experimentally in the E4 state. E4 models with two bridging hydrides between Fe2 and Fe6 show much more favorable N2 binding than other models. When two hydrides coordinate to the same Fe ion, an increased ligand-field splitting due to octahedral coordination at either Fe2 or Fe6 is found. This altered ligand field makes it easier for the Fe ion to acquire a spin-paired electronic structure with doubly occupied dxz and dyz orbitals that backbond to a terminal N2 ligand. Within this model for N2 binding, both Fe2 and Fe6 emerge as possible binding site scenarios. Due to steric effects involving the His195 residue, affecting both the N2 ligand and the terminal SH- group, distinguishing between Fe2 and Fe6 is difficult; furthermore, the binding depends on the protonation state of His195.

The mechanism for N2 activation in the E4 - state of nitrogenase.

Nitrogenases take nitrogen from the air and reduce it to ammonia. It has long been known that N2 becomes activated after four reductions in the catalytic cycle, in the E4 state. Several mechanisms for the activation have been suggested. In the present study a previous mechanism has been revised based on recent experimental findings. In the present mechanism N2H2 is formed in E4. As in the previously suggested mechanism, there are four initial reductions before catalysis (the A-states), after which a sulfide is released and the first state in catalysis (E0) is formed. In E4, N2 becomes bound and protonated in the Fe1, Fe2, Fe4 region, in which the hydrides have left two electrons. The rate-limiting step is the formation of N2H by a hydrogen atom transfer from Cys275 to N2 bound to Fe4, concerted with an additional electron transfer from the cofactor. The mechanism fulfills all requirements set by experiments. The activation of N2 is preceded by a formation of H2 from two hydrides, the carbide is kinetically hindered from being protonated, the E4 state is reversible. An important aspect is the presence of a water molecule in the Fe2, Fe6 region. The non-allowed formations of H2 from a hydride and a proton have been investigated and found to have higher barriers than the allowed formation of H2 from two hydrides.

N2 binding to the E0-E4 states of nitrogenase.

Nitrogenase is the only enzyme that can convert N2 into NH3. The reaction requires the addition of eight electrons and protons to the enzyme and the mechanism is normally described by nine states, E0-E8, differing in the number of added electrons. Experimentally, it is known that three or four electrons need to be added before the enzyme can bind N2. We have used combined quantum mechanical and molecular mechanics methods to study the binding of N2 to the E0-E4 states of nitrogenase, using four different density functional theory (DFT) methods. We test many different structures for the E2-E4 states and study binding both to the Fe2 and Fe6 ions of the active-site FeMo cluster. Unfortunately, the results depend quite strongly on the DFT methods. The TPSS method gives the strongest bonding and prefers N2 binding to Fe6. It is the only method that reproduces the experimental observation of unfavourable binding to the E0-E2 states and favourable binding to E3 and E4. The other three methods give weaker binding, preferably to Fe2. B3LYP strongly favours structures with the central carbide ion triply protonated. The other three methods suggest that states with the S2B ligand dissociated from either Fe2 or Fe6 are competitive for the E2-E4 states. Moreover, such structures with two hydride ions both bridging Fe2 and Fe6 are the best models for E4 and also for the N2-bound E3 and E4 states. However, for E4, other structures are often close in energy, e.g. structures with one of the hydride ions bridging instead Fe3 and Fe7. Finally, we find no support for the suggestion that reductive elimination of H2 from the two bridging hydride ions in the E4 state would enhance the binding of N2.

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.

Mechanistic Insights into Nitrogenase FeMo-Cofactor Catalysis through a Steady-State Kinetic Model.

Mo-nitrogenase catalyzes the challenging N2-to-NH3 reduction. This complex reaction proceeds through a series of intermediate states (En) of its active site FeMo-cofactor. An understanding of the kinetics of the conversion between En states is central to defining the mechanism of nitrogenase. Here, rate constants of key steps have been determined through a steady-state kinetic model with fits to experimental data. The model reveals that the rate for H2 formation from the early electron populated state E2(2H) is much slower than that from the more reduced E4(4H) state. Further, it is found that the competing reactions of H2 formation and N2 binding at the E4(4H) state occur with equal rate constants. The H2-dependent reverse reaction of the N2 binding step is found to have a rate constant of 5.5 ± 0.2 (atm H2)-1 s-1 (7.2 ± 0.3 (mM H2)-1 s-1). Importantly, the reduction of N2 bound to FeMo-cofactor proceeds with a rate constant of 1 ± 0.1 s-1, revealing a previously unrecognized slow step in the Mo-nitrogenase catalytic cycle associated with the chemical transformation of N2 to 2 NH3. Finally, the populations of En states under different reaction conditions are predicted, providing a powerful tool to guide the spectroscopic and mechanistic studies of Mo-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+ (E4(4H) state) before N2 can bind and be reduced. Earlier studies demonstrated that this E4(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 N2 binding is mechanistically coupled to reductive-elimination of the two hydrides that produce H2. The state having acquired 2e-/2H+ (E2(2H)) correspondingly contains one hydride with a resting-state core redox level. How the cofactor accommodates addition of the first e-/H+ (E1(H) state) is unknown. The Fe-nitrogenase FeFe-cofactor was used to address this question because it is EPR-active in the E1(H) state, unlike the FeMo-cofactor of Mo-nitrogenase, thus allowing characterization by EPR spectroscopy. The freeze-trapped E1(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 E1(H)*) with g = [2.009, 1.950, 1.860], which results in a photostationary state, with E1(H)* relaxing to E1(H) at temperatures above 145 K. An H/D kinetic isotope effect of 2.4 accompanies the 12 K E1(H)/E1(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.

Thermodynamically Favourable States in the Reaction of Nitrogenase without Dissociation of any Sulfide Ligand.

We have used combined quantum mechanical and molecular mechanical (QM/MM) calculations to study the reaction mechanism of nitrogenase, assuming that none of the sulfide ligands dissociates. To avoid the problem that there is no consensus regarding the structure and protonation of the E4 state, we start from a state where N2 is bound to the cluster and is protonated to N2H2, after dissociation of H2. We show that the reaction follows an alternating mechanism with HNNH (possibly protonated to HNNH2) and H2NNH2 as intermediates and the two NH3 products dissociate at the E7 and E8 levels. For all intermediates, coordination to Fe6 is preferred, but for the E4 and E8 intermediates, binding to Fe2 is competitive. For the E4, E5 and E7 intermediates we find that the substrate may abstract a proton from the hydroxy group of the homocitrate ligand of the FeMo cluster, thereby forming HNNH2, H2NNH2 and NH3 intermediates. This may explain why homocitrate is a mandatory component of nitrogenase. All steps in the suggested reaction mechanism are thermodynamically favourable compared to protonation of the nearby His‐195 group and in all cases, protonation of the NE2 atom of the latter group is preferred.

The active E4 structure of nitrogenase studied with different DFT functionals.

The present study concerns the technical aspects of obtaining the energetics for the E4 state of nitrogenase, the enzyme that fixes N2 in nature. EPR experiments have shown that the critical E4 structure that activates N2 should contain two bridging hydrides in the FeMo‐cofactor. It is furthermore in equilibrium with a structure where the two hydrides have been released and N2 binds. These observations led to the suggestion that E4 should have two bridging hydrides and two protonated sulfides. It is important to note that the structure for E4 has not been determined, but only suggested. For a long time, no DFT study led to the suggested structure, independent of which functional was used. However, in two recent DFT studies a good agreement with the experimental suggestion was claimed to have been obtained. In one of them the TPSS functional was used. That was the first out of 11 functionals tried that led to the experimentally suggested structure. In the second of the recent DFT studies, a similar conclusion was reached using the TPSSh functional. The conclusions in the recent studies have here been studied in detail, by calculating a critical energetic value strongly implied by the same EPR experiments. Both the TPSS and TPSSh functionals have been used. The present calculations suggest that those DFT functionals would not lead to agreement with the experimental EPR results either.

What Is the Structure of the E4 Intermediate in Nitrogenase?

Nitrogenase is the only enzyme that can cleave the strong triple bond in N2. The active site contains a complicated MoFe7S9C cluster. It is believed that it needs to accept four protons and electrons, forming the E4 state, before it can bind N2. However, there is no consensus on the atomic structure of the E4 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 E4 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.

Reactivity toward Unsaturated Small Molecules of Thiolate-Bridged Diiron Hydride Complexes.

In the presence of 1 equiv of tBuNC, the homolytic cleavage of the FeIII-H bond in the diiron terminal hydride complex [Cp*Fe( t-H)(μ-η2:η4-bdt)FeCp*][BF4] (1[BF4]) smoothly took place to release 1/2 H2, followed by binding of a tBuNC group to the unsaturated FeII center. Interestingly, upon exposure of 1[BF4] to 1 atm of acetylene, the isomerization process of the hydride ligand from the terminal to bridging coordination site was unaffected. Upon treatment of the diiron hydride bridged complex 2[BF4] with acetylene at 30 °C, two FeIII-H bonds were broken, and then an acetylene molecule was coordinated to the diiron centers in a novel μ-η2:η2 side-on fashion. In the above reaction system, the hydride ligands whether terminal or bridging all play a role as the electron donor for the reduction of the diiron centers from FeIIIFeIII to FeIIIFeII. These reaction patterns are reminiscent of the vital E4 state responsible for N2 binding and H2 liberation in the catalytic cycle of nitrogenase, which contains two {Fe-H-Fe} motifs as electron reservoirs for the reduction of the iron centers. Differently, when treating 1[BF4] with TMSN3, the terminal hydride ligand was inserted into the azide group to give a diiron amide complex 4[BF4] in moderate yield.

Protonation and Reduction of the FeMo Cluster in Nitrogenase Studied by Quantum Mechanics/Molecular Mechanics (QM/MM) Calculations.

We have performed a systematic computational study of the relative energies of possible protonation states of the FeMo cluster in nitrogenase in the E0-E4 states, i.e., the resting state and states with 1-4 electrons and protons added but before N2 binds. We use the combined quantum mechanics and molecular mechanics (QM/MM) approach, including the complete solvated heterotetrameric enzyme in the calculations. The QM system consisted of 112 atoms, i.e., the full FeMo cluster, as well all groups forming hydrogen bonds to it within 3.5 Å. It was treated with either the TPSS-D3 or B3LYP-D3 methods with the def2-SV(P) or def2-TZVPD basis sets. For each redox state, we calculated relative energies of at least 50 different possible positions for the proton, added to the most stable protonation state of the level with one electron less. We show quite conclusively that the resting E0 state is not protonated using quantum refinement and by comparing geometries to the crystal structure. The E1 state is protonated on S2B, in agreement with most previous computational studies. However, for the E2-E4 states, the two QM methods give diverging results, with relative energies that differ by over 300 kJ/mol for the most stable E4 states. TPSS favors hydride ions binding to the Fe ions. The first bridges Fe2 and Fe6, whereas the next two bind terminally to either Fe4, Fe5, or Fe6 with nearly equal energies. On the other hand, B3LYP disfavors hydride ions and instead suggests that 1-3 protons bind to the central carbide ion.

Kinetic Understanding of N2 Reduction versus H2 Evolution at the E4(4H) Janus State in the Three Nitrogenases.

The enzyme nitrogenase catalyzes the reduction of N2 to ammonia but also that of protons to H2. These reactions compete at the mechanistically central 'Janus' intermediate, denoted E4(4H), which has accumulated 4e-/4H+ as two bridging Fe-H-Fe hydrides on the active-site cofactor. This state can lose e-/H+ by hydride protonolysis (HP) or become activated by reductive elimination ( re) of the two hydrides and bind N2 with H2 loss, yielding an E4(2N2H) state that goes on to generate two NH3 molecules. Thus, E4(4H) represents the key branch point for these competing reactions. Here, we present a steady-state kinetic analysis that precisely describes this competition. The analysis demonstrates that steady-state, high-electron flux turnover overwhelmingly populates the E4 states at the expense of less reduced states, quenching HP at those states. The ratio of rate constants for E4(4H) hydride protonolysis ( kHP) versus reductive elimination ( kre) provides a sensitive measure of competition between these two processes and thus is a central parameter of nitrogenase catalysis. Analysis of measurements with the three nitrogenase variants (Mo-nitrogenase, V-nitrogenase, and Fe-nitrogenase) reveals that at a fixed N2 pressure their tendency to productively react with N2 to produce two NH3 molecules and an accompanying H2, rather than diverting electrons to the side reaction, HP production of H2, decreases with their ratio of rate constants, k re/ kHP: Mo-nitrogenase, 5.1 atm-1; V-nitrogenase, 2 atm-1; and Fe-nitrogenase, 0.77 atm-1 (namely, in a 1:0.39:0.15 ratio). Moreover, the lower catalytic effectiveness of the alternative nitrogenases, with more H2 production side reaction, is not caused by a higher kHP but by a significantly lower k re.

Is there computational support for an unprotonated carbon in the E4 state of nitrogenase?

In the key enzyme for nitrogen fixation in nature, nitrogenase, the active site has a metal cluster with seven irons and one molybdenum bound by bridging sulfurs. Surprisingly, there is also a carbon in the center of the cluster, with a role that is not known. A mechanism has been suggested experimentally, where two hydrides leave as a hydrogen molecule in the critical E4 state. A structure with two hydrides, two protonated sulfurs and an unprotonated carbon has been suggested for this state. Rather recently, DFT calculations supported the experimental mechanism but found an active state where the central carbon is protonated all the way to CH3 . Even more recently, another DFT study was made that instead supported the experimentally suggested structure. To sort out the origin of these quite different computational results, additional calculations have here been performed using different DFT functionals. The conclusion from these calculations is very clear and shows no computational support for an unprotonated carbon in E4 . © 2017 Wiley Periodicals, Inc.