A Mechanism for Nitrogenase Including Loss of a Sulfide

The nitrogenase mechanism: new roles for the dangler?

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

Understanding non-reducible N2 in the mechanism of Mo-nitrogenase

In my proposed mechanism of Mo-nitrogenase there are two roles for separate N2 molecules. One N2 diffuses into the reaction zone between Fe2 and Fe6 where a strategic gallery of H atoms can capture N2 to form the Fe-bound HNNH intermediate which is then progressively hydrogenated through intermediates containing HNNH2, NH and NH2 entities and then two NH3 in sequence. The second N2 can be parked in an N2-pocket about 3.2 Å from Fe2 or bind end-on at the exo coordination site of Fe2. This second N2 is outside the reaction zone, not exposed to H atom donors, and so is 'non-reducible'. Here density functional calculations using a 485+ atom model describe the thermodynamics for non-reducible N2 moving between the N2-pocket and the exo-Fe2 position, for the resting state and 19 intermediates in the mechanism. The entropy component is estimated and included. The result is that for all intermediates with ligation by H or NHx at the endo-Fe2 position the free energy for association of non-reducible N2 at exo-Fe2 is negative. There remains some uncertainty about the status of exo-Fe2-N2 during the step in which H2 exchanges with the incoming reducible N2, where at least two unbound molecules are present. At Fe2 it is evident that attainment of octahedral coordination stereochemistry dominates the binding thermodynamics for non-reducible N2. Possibilities for experimental support of these computational conclusions are discussed.

The mechanism of Mo-nitrogenase: from N2 capture to first release of NH3

Mo-nitrogenase hydrogenates N2 to NH3. This report continues from the previous paper [I. Dance, Dalton Trans., 2024, 53, 14193-14211] that described how the active site FeMo-co of the enzyme is uniquely able to capture and activate N2, forming a key intermediate with Fe-bound HNNH. Density functional simulations with a 485+ atom model of the active site and its surroundings are used to describe here the further reactions of this HNNH intermediate. The first step is hydrogenation to form HNNH2 bridging Fe2 and Fe6. Then a single-step reaction breaks the N-N bond, generating an Fe2-NH-Fe6 bridge and forming NH3 bound to Fe6. Then NH3 dissociates from Fe6. Reaction potential energies and kinetic barriers for all steps are reported for the most favourable electronic states of the system. The steps that follow the Fe2-NH-Fe6 intermediate, forming and dissociating the second NH3, and regenerating the resting state of the enzyme, are outlined. These results provide an interpretation of the recent steady-state kinetics data and analysis by Harris et al., [Biochemistry, 2022, 61, 2131-2137] who found a slow step after the formation of the HNNH intermediate. The calculated potential energy barriers for the HNNH2 → NH + NH3 reaction (30-36 kcal mol-1) are larger than the potential energy barriers for the N2 → HNNH reaction (19-29 kcal mol-1). I propose that the post-HNNH slow step identified kinetically is the key HNNH2 → NH + NH3 reaction described here. This step and the N2-capture step are the most difficult in the conversion of N2 to 2NH3. The steps in the complete mechanism still to be computationally detailed are relatively straightforward.

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.

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.

CO Oxidation Mechanism of Silver-Substituted Mo/Cu CO-Dehydrogenase - Analogies and Differences to the Native Enzyme

The aerobic oxidation of carbon monoxide to carbon dioxide is catalysed by the Mo/Cu-containing CO-dehydrogenase enzyme in the soil bacterium Oligotropha carboxidovorans, enabling the organism to grow on the small gas molecule as carbon and energy source. It was shown experimentally that silver can be substituted for copper in the active site of Mo/Cu CODH to yield a functional enzyme. In this study, we employed QM/MM calculations to investigate whether the reaction mechanism of the silver-substituted enzyme is similar to that of the native enzyme. Our results suggest that the Ag-substituted enzyme can oxidize CO and release CO2 following the same reaction steps as the native enzyme, with a computed rate-limiting step of 10.4 kcal/mol, consistent with experimental findings. Surprisingly, lower activation energies for C-O bond formation have been found in the presence of silver. Furthermore, comparison of rate constants for reduction of copper- and silver-containing enzymes suggests a discrepancy in the transition state stabilization upon silver substitution. We also evaluated the effects that differences in the water-active site interaction may exert on the overall energy profile of catalysis. Finally, the formation of a thiocarbonate intermediate along the catalytic pathway was found to be energetically unfavorable for the Ag-substituted enzyme. This finding aligns with the hypothesis proposed for the wild-type form, suggesting that the creation of such species may not be necessary for the enzymatic catalysis of CO oxidation.

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

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.

Enzymatic Fischer-Tropsch-Type Reactions

The Fischer-Tropsch (FT) process converts a mixture of CO and H2 into liquid hydrocarbons as a major component of the gas-to-liquid technology for the production of synthetic fuels. Contrary to the energy-demanding chemical FT process, the enzymatic FT-type reactions catalyzed by nitrogenase enzymes, their metalloclusters, and synthetic mimics utilize H+ and e- as the reducing equivalents to reduce CO, CO2, and CN- into hydrocarbons under ambient conditions. The C1 chemistry exemplified by these FT-type reactions is underscored by the structural and electronic properties of the nitrogenase-associated metallocenters, and recent studies have pointed to the potential relevance of this reactivity to nitrogenase mechanism, prebiotic chemistry, and biotechnological applications. This review will provide an overview of the features of nitrogenase enzymes and associated metalloclusters, followed by a detailed discussion of the activities of various nitrogenase-derived FT systems and plausible mechanisms of the enzymatic FT reactions, highlighting the versatility of this unique reactivity while providing perspectives onto its mechanistic, evolutionary, and biotechnological implications.

How Protons Move in Enzymes─The Case of Nitrogenase

When moving protons in enzymes, water molecules are often used as intermediates. The water molecules used are not necessarily seen in the crystal structures if they move around at high rates. In a different situation, for metal containing cofactors in enzymes, it is sometimes necessary to move protons on the cofactor from the position they enter the cofactor to another position where the energy is lower. That is, for example, the situation in nitrogenase. In recent studies on that enzyme, prohibitively high barriers were sometimes found for transferring protons, and that was used as a strong argument against mechanisms where a sulfide is lost in the mechanism. A high barrier could be due to nonoptimal distances and angles at the transition state. In the present study, possibilities are investigated to use water molecules to reduce these barriers. The study is very general and could have been done for many other enzymes. The effect of water was found to be very large in the case of nitrogenase with a lowering of one barrier from 15.6 kcal/mol down to essentially zero. It is concluded that the effect of water molecules must be taken into account for meaningful results.

The HD Reaction of Nitrogenase: a Detailed Mechanism

Nitrogenase is the enzyme that converts N2 to NH3 under ambient conditions. The chemical mechanism of this catalysis at the active site FeMo-co [Fe7 S9 CMo(homocitrate)] is unknown. An obligatory co-product is H2 , while exogenous H2 is a competitive inhibitor. Isotopic substitution using exogenous D2 revealed the N2 -dependent reaction D2 +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 D2 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 D2 , 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 N2 -dependence results from promotional N2 bound at exo-Fe2; (c) different N2 binding Km for the HD reaction and the NH3 formation reaction results from involvement of two different sites; (d) inhibition by CO; (e) the non-occurrence of 2HD→H2 +D2 results from the synchronicity of the two transfers of H to D2 ; (f) inhibition of HD production at high pN2 is by competitive binding of N2 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.

Computational modeling of redox enzymes

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