Diazene (HN=NH) is a substrate for nitrogenase: insights into the pathway of N2 reduction

Investigating the Molybdenum Nitrogenase Mechanistic Cycle Using Spectroelectrochemistry

Molybdenum nitrogenase plays a crucial role in the biological nitrogen cycle by catalyzing the reduction of dinitrogen (N2) to ammonia (NH3) under ambient conditions. However, the underlying mechanisms of nitrogenase catalysis, including electron and proton transfer dynamics, remain only partially understood. In this study, we covalently attached molybdenum nitrogenase (MoFe) to gold electrodes and utilized surface-enhanced infrared absorption spectroscopy (SEIRA) coupled with electrochemistry techniques to investigate its catalytic mechanism. Our biohybrid system enabled electron transfer via a mild mediator, likely mimicking the natural electron flow through the P-cluster to FeMoco, the enzyme's active site. For the first time, we experimentally observed both terminal and bridging S-H stretching frequencies, resulting from the protonation of bridging sulfides in FeMoco during turnover conditions providing direct evidence of their role in catalysis. These experimental observations are further supported by QM/MM calculations. Additionally, we investigated CO inhibition, demonstrating both CO binding and unbinding dynamics under electrochemical conditions. These insights not only advance our understanding of the mechanistic cycle of molybdenum nitrogenase but also establish a foundation for studying alternative nitrogenases, including vanadium and iron nitrogenases.

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.

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.

Reaction Mechanism for CO Reduction by Mo-Nitrogenase Studied by QM/MM

We have studied the conversion of two molecules of carbon monoxide to ethylene catalyzed by nitrogenase. We start from a recent crystal structure showing the binding of two carbon monoxide molecules to nitrogenase and employ the combined quantum mechanics and molecular mechanics approach. Our results indicate that the reaction is possible only if S2B dissociates as H2S (i.e., the charge of the FeMo cluster remains the same as in the E0 state, indicating that the Fe ions are formally reduced two steps when CO binds). Eight electrons and protons are needed for the reaction, and our mechanism suggests that the first four bind alternatively to the two carbon atoms. The C-C bond formation takes place already after the first protonation (in the E3 state). The next two protons bind to the same O atom, which then dissociates as water. In the same state (E8), the second C-O bond is cleaved, forming the ethylene product. The last two electrons and protons are used to form a water molecule that can be exchanged by S2B or by two CO molecules to start a new reaction cycle.

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.

Hole-scavenging in photo-driven N2 reduction catalyzed by a CdS-nitrogenase MoFe protein biohybrid system

The light-driven reduction of dinitrogen (N2) to ammonia (NH3) catalyzed by a cadmium sulfide (CdS) nanocrystal‑nitrogenase MoFe protein biohybrid is dependent on a range of different factors, including an appropriate hole-scavenging sacrificial electron donor (SED). Here, the impact of different SEDs on the overall rate of N2 reduction catalyzed by a CdS quantum dot (QD)-MoFe protein system was determined. The selection of SED was guided by several goals: (i) molecules with standard reduction potentials sufficient to reduce the oxidized CdS QD, (ii) molecules that do not absorb the excitation wavelength of the CdS QD, and (iii) molecules that could be readily reduced by sustainable processes. Earlier studies utilized buffer molecules or ascorbic acid as the SED. The effectiveness of ascorbic acid as SED was compared to dithionite (DT), triethanolamine (TEOA), and hydroquinone (HQ) across a range of concentrations in supporting N2 reduction to NH3 in a CdS QD-MoFe protein photocatalytic system. It was found that TEOA supported N2 reduction rates comparable to those observed for dithionite and ascorbic acid. HQ was found to support significantly higher rates of N2 reduction compared to the other SEDs at a concentration of 50 mM. A comparison of the rates of N2 reduction by the biohybrid complex to the standard reduction potential (Eo) of the SEDs reveals that Eo is not the only factor impacting the efficiency of hole-scavenging. These findings reveal the importance of the SED properties for improving the efficiency of hole-scavenging in the light-driven N2 reduction reaction catalyzed by a CdS QD-MoFe protein hybrid.

Low-Valent Transition Metalate Anions in Synthesis, Small Molecule Activation, and Catalysis

This review surveys the synthesis and reactivity of low-oxidation state metalate anions of the d-block elements, with an emphasis on contributions reported between 2006 and 2022. Although the field has a long and rich history, the chemistry of transition metalate anions has been greatly enhanced in the last 15 years by the application of advanced concepts in complex synthesis and ligand design. In recent years, the potential of highly reactive metalate complexes in the fields of small molecule activation and homogeneous catalysis has become increasingly evident. Consequently, exciting applications in small molecule activation have been developed, including in catalytic transformations. This article intends to guide the reader through the fascinating world of low-valent transition metalates. The first part of the review describes the synthesis and reactivity of d-block metalates stabilized by an assortment of ligand frameworks, including carbonyls, isocyanides, alkenes and polyarenes, phosphines and phosphorus heterocycles, amides, and redox-active nitrogen-based ligands. Thereby, the reader will be familiarized with the impact of different ligand types on the physical and chemical properties of metalates. In addition, ion-pairing interactions and metal-metal bonding may have a dramatic influence on metalate structures and reactivities. The complex ramifications of these effects are examined in a separate section. The second part of the review is devoted to the reactivity of the metalates toward small inorganic molecules such as H2, N2, CO, CO2, P4 and related species. It is shown that the use of highly electron-rich and reactive metalates in small molecule activation translates into impressive catalytic properties in the hydrogenation of organic molecules and the reduction of N2, CO, and CO2. The results discussed in this review illustrate that the potential of transition metalate anions is increasingly being tapped for challenging catalytic processes with relevance to organic synthesis and energy conversion. Therefore, it is hoped that this review will serve as a useful resource to inspire further developments in this dynamic research field.

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

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.

On the Temperature Sensitivity of Electrochemical Reaction Thermodynamics

Herein, we report a method to estimate the thermodynamic potentials of electrochemical reactions at different temperatures. We use a two-term Taylor series approximation of thermodynamic potential as a function of temperature, and we calculate the temperature sensitivity for a family of twenty seven known half reactions. We further analyze pairs of cathode and anode half-cells to pinpoint optimal voltage matches and discuss implications of changes in temperature on overall cell voltages. Using these observations, we look forward to increased interest in temperature and idealized half-reaction pairing as experimental choices for the optimization of electrochemical processes.

Light Alters the NH3 vs N2 H4 Product Profile in Iron-catalyzed Nitrogen Reduction via Dual Reactivity from an Iron Hydrazido (Fe=NNH2 ) Intermediate

Whereas synthetically catalyzed nitrogen reduction (N2 R) to produce ammonia is widely studied, catalysis to instead produce hydrazine (N2 H4 ) has received less attention despite its considerable mechanistic interest. Herein, we disclose that irradiation of a tris(phosphine)borane (P3 B ) Fe catalyst, P3 B Fe+ , significantly alters its product profile to increase N2 H4 versus NH3 ; P3 B Fe+ is otherwise known to be highly selective for NH3 . We posit a key terminal hydrazido intermediate, P3 B Fe=NNH2 , as selectivity-determining. Whereas its singlet ground state undergoes protonation to liberate NH3 , a low-lying triplet excited state leads to reactivity at Nα and formation of N2 H4 . Associated electrochemical and spectroscopic studies establish that N2 H4 lies along a unique product pathway; NH3 is not produced from N2 H4 . Our findings are distinct from the canonical mechanism for hydrazine formation, which proceeds via a diazene (HN=NH) intermediate and showcase light as a tool to tailor selectivity.

Many-body analysis and spectroscopic characterization of diazene oligomers: A theoretical study

The present study reports the many-body analysis and spectroscopic characterization of linear and cyclic diazene oligomers in gas and water solvent states. The oligomers of diazene from monomer to pentamer have been considered for the study. The spectroscopic studies such as geometrical parameters, infrared spectra, electronic absorption spectra, and natural transition orbitals (NTOs) were reported. Many-body analysis techniques have been implemented to study the interactions among the diazene oligomers. These calculations have been performed using exchange and correlation functional (B3LYP) and 6-311++G (d,p) basis set. The geometrical parameters and infrared modes of monomer diazene in the gas state are well-matched with the available experimental determinations at this level of theory. A significant change in vibrational modes of linear and cyclic diazene oligomers has been observed in the gas phase-to-water solvent state. The time-dependent density functional theory (TD-DFT) has been used to calculate the electronic absorption spectra of diazene oligomers. The Wavelength of electronic transitions, oscillator strength, and HOMO to LUMO gap has been reported. Many-body analysis shows that two-, three-, four-, and five-body energies have a remarkable contribution to the binding energy in addition to relaxation energies. All these calculations have been performed using Gaussian 16 program package.

Near ambient N2 fixation on solid electrodes versus enzymes and homogeneous catalysts

The Mo/Fe nitrogenase enzyme is unique in its ability to efficiently reduce dinitrogen to ammonia at atmospheric pressures and room temperature. Should an artificial electrolytic device achieve the same feat, it would revolutionize fertilizer production and even provide an energy-dense, truly carbon-free fuel. This Review provides a coherent comparison of recent progress made in dinitrogen fixation on solid electrodes, homogeneous catalysts and nitrogenases. Specific emphasis is placed on systems for which there is unequivocal evidence that dinitrogen reduction has taken place. By establishing the cross-cutting themes and synergies between these systems, we identify viable avenues for future research.

Quantum Mechanical Calculations of Redox Potentials of the Metal Clusters in Nitrogenase

We have calculated redox potentials of the two metal clusters in Mo-nitrogenase with quantum mechanical (QM) calculations. We employ an approach calibrated for iron-sulfur clusters with 1-4 Fe ions, involving QM-cluster calculations in continuum solvent and large QM systems (400-500 atoms), based on structures from combined QM and molecular mechanics (QM/MM) geometry optimisations. Calculations on the P-cluster show that we can reproduce the experimental redox potentials within 0.33 V. This is similar to the accuracy obtained for the smaller clusters, although two of the redox reactions involve also proton transfer. The calculated P1+/PN redox potential is nearly the same independently of whether P1+ is protonated or deprotonated, explaining why redox titrations do not show any pH dependence. For the FeMo cluster, the calculations clearly show that the formal oxidation state of the cluster in the resting E0 state is MoIIIFe3IIFe4III , in agreement with previous experimental studies and QM calculations. Moreover, the redox potentials of the first five E0-E4 states are nearly constant, as is expected if the electrons are delivered by the same site (the P-cluster). However, the redox potentials are insensitive to the formal oxidation states of the Fe ion (i.e., whether the added protons bind to sulfide or Fe ions). Finally, we show that the later (E4-E8) states of the reaction mechanism have redox potential that are more positive (i.e., more exothermic) than that of the E0/E1 couple.

QM/MM Study of Partial Dissociation of S2B for the E2 Intermediate of Nitrogenase

Nitrogenase is the only enzyme that can cleave the triple bond in N2, making nitrogen available for all lifeforms. Previous computational studies have given widely diverging results regarding the reaction mechanism of the enzyme. For example, some recent studies have suggested that one of the μ2-bridging sulfide ligands (S2B) may dissociate from one of the Fe ions when protonated in the doubly reduced and protonated E2 state, whereas other studies indicated that such half-dissociated states are unfavorable. We have examined how the relative energies of 26 structures of the E2 state depend on details of combined quantum mechanical and molecular mechanical (QM/MM) calculations. We show that the selection of the broken-symmetry state, the basis set, relativistic effects, the size of the QM system, relaxation of the surroundings, and the conformations of the bound protons may affect the relative energies of the various structures by up to 12, 22, 9, 20, 37, and 33 kJ/mol, respectively. However, they do not change the preferred type of structures. On the other hand, the choice of the DFT functional strongly affects the preferences. The hybrid B3LYP functional strongly prefers doubly protonation of the central carbide ion, but such a structure is not consistent with experimental EPR data. Other functionals suggest structures with a hydride ion, in agreement with the experiments, and show that the ion bridges between Fe2 and Fe6. Moreover, there are two structures of the same type that are degenerate within 1-5 kJ/mol, in agreement with the observation of two EPR signals. However, the pure generalized gradient approximation (GGA) functional TPSS favors structures with a protonated S2B also bridging Fe2 and Fe6, whereas r2SCAN (meta-GGA) and TPSSh (hybrid) prefer structures with S2B dissociated from Fe2 (but remaining bound to Fe6). The energy difference between the two types of structure is so small (7-18 kJ/mol) that both types need to be considered in future investigations of the mechanism of nitrogenase.

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 E₄ state, we start from a state where N₂ is bound to the cluster and is protonated to N₂H₂, after dissociation of H₂. We show that the reaction follows an alternating mechanism with HNNH (possibly protonated to HNNH₂) and H₂NNH₂ as intermediates and the two NH₃ products dissociate at the E₇ and E₈ levels. For all intermediates, coordination to Fe6 is preferred, but for the E₄ and E₈ intermediates, binding to Fe2 is competitive. For the E₄, E₅ and E₇ intermediates we find that the substrate may abstract a proton from the hydroxy group of the homocitrate ligand of the FeMo cluster, thereby forming HNNH₂, H₂NNH₂ and NH₃ 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.

Identification of M-NH2 -NH2 Intermediate and Rate Determining Step for Nitrogen Reduction with Bioinspired Sulfur-Bonded FeW Catalyst

The multimetallic sulfur-framework catalytic site of biological nitrogenases allows the efficient conversion of dinitrogen (N2 ) to ammonia (NH3 ) under ambient conditions. Inspired by biological nitrogenases, a bimetallic sulfide material (FeWSx @FeWO4 ) was synthesized as a highly efficient N2 reduction (NRR) catalyst by sulfur substitution of the surface of FeWO4 nanoparticles. Thus prepared FeWSx @FeWO4 catalysts exhibit a relatively high NH3 production rate of 30.2 ug h-1  mg-1cat and a Faraday efficiency of 16.4 % at -0.45 V versus a reversible hydrogen electrode in a flow cell; these results have been confirmed via purified 15 N2 -isotopic labeling experiments. In situ Raman spectra and hydrazine reduction kinetics analysis revealed that the reduction of undissociated hydrazine intermediates (M-NH2 -NH2 ) on the surface of the bimetallic sulfide catalyst is the rate-determing step for the NRR process. Therefore, this work can provide guidance for elucidating the structure-activity relationship of NRR catalysts.

Reactivity, Mechanism, and Assembly of the Alternative Nitrogenases

Biological nitrogen fixation is catalyzed by the enzyme nitrogenase, which facilitates the cleavage of the relatively inert triple bond of N2. Nitrogenase is most commonly associated with the molybdenum-iron cofactor called FeMoco or the M-cluster, and it has been the subject of extensive structural and spectroscopic characterization over the past 60 years. In the late 1980s and early 1990s, two "alternative nitrogenase" systems were discovered, isolated, and found to incorporate V or Fe in place of Mo. These systems are regulated by separate gene clusters; however, there is a high degree of structural and functional similarity between each nitrogenase. Limited studies with the V- and Fe-nitrogenases initially demonstrated that these enzymes were analogously active as the Mo-nitrogenase, but more recent investigations have found capabilities that are unique to the alternative systems. In this review, we will discuss the reactivity, biosynthetic, and mechanistic proposals for the alternative nitrogenases as well as their electronic and structural properties in comparison to the well-characterized Mo-dependent system. Studies over the past 10 years have been particularly fruitful, though key aspects about V- and Fe-nitrogenases remain unexplored.

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.

Catalytic N2-to-NH3 (or -N2H4) Conversion by Well-Defined Molecular Coordination Complexes

Nitrogen fixation, the six-electron/six-proton reduction of N₂, to give NH₃, is one of the most challenging and important chemical transformations. Notwithstanding the barriers associated with this reaction, significant progress has been made in developing molecular complexes that reduce N₂ into its bioavailable form, NH₃. This progress is driven by the dual aims of better understanding biological nitrogenases and improving upon industrial nitrogen fixation. In this review, we highlight both mechanistic understanding of nitrogen fixation that has been developed, as well as advances in yields, efficiencies, and rates that make molecular alternatives to nitrogen fixation increasingly appealing. We begin with a historical discussion of N₂ functionalization chemistry that traverses a timeline of events leading up to the discovery of the first bona fide molecular catalyst system and follow with a comprehensive overview of d-block compounds that have been targeted as catalysts up to and including 2019. We end with a summary of lessons learned from this significant research effort and last offer a discussion of key remaining challenges in the field.

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.

Reduction of Substrates by Nitrogenases

Nitrogenase is the enzyme that catalyzes biological N₂ reduction to NH₃. This enzyme achieves an impressive rate enhancement over the uncatalyzed reaction. Given the high demand for N₂ fixation to support food and chemical production and the heavy reliance of the industrial Haber-Bosch nitrogen fixation reaction on fossil fuels, there is a strong need to elucidate how nitrogenase achieves this difficult reaction under benign conditions as a means of informing the design of next generation synthetic catalysts. This Review summarizes recent progress in addressing how nitrogenase catalyzes the reduction of an array of substrates. New insights into the mechanism of N₂ and proton reduction are first considered. This is followed by a summary of recent gains in understanding the reduction of a number of other nitrogenous compounds not considered to be physiological substrates. Progress in understanding the reduction of a wide range of C-based substrates, including CO and CO₂, is also discussed, and remaining challenges in understanding nitrogenase substrate reduction are considered.

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.

Aerobic nitrogen-fixing bacteria for hydrogen and ammonium production: current state and perspectives

Biological nitrogen fixation (BNF) is accomplished through the action of the oxygen-sensitive enzyme nitrogenase. One unique caveat of this reaction is the inclusion of hydrogen gas (H₂) evolution as a requirement of the reaction mechanism. In the absence of nitrogen gas as a substrate, nitrogenase will reduce available protons to become a directional ATP-dependent hydrogenase. Aerobic nitrogen-fixing microbes are of particular interest, because these organisms have evolved to perform these reactions with oxygen-sensitive enzymes in an environment surrounded by oxygen. The ability to maintain a functioning nitrogenase in aerobic conditions facilitates the application of these organisms under conditions where most anaerobic nitrogen fixers are excluded. In recent years, questions related to the potential yields of the nitrogenase-derived products ammonium and H₂ have grown more approachable to experimentation based on efforts to construct increasingly more complicated strains of aerobic nitrogen fixers such as the obligate aerobe Azotobacter vinelandii. This mini-review provides perspectives of recent and historical efforts to understand and quantify the yields of ammonium and H₂ that can be obtained through the model aerobe A. vinelandii, and outstanding questions that remain to be answered to fully realize the potential of nitrogenase in these applications with model aerobic bacteria.

Extremely large differences in DFT energies for nitrogenase models

Nitrogenase is the only enzyme that can cleave the triple bond in N2, making nitrogen avaiable for other organisms. It contains a complicated MoFe7S9C(homocitrate) cluster in its active site. Many computational studies with density-functional theory (DFT) of the nitrogenase enzyme have been presented, but they do not show any consensus - they do not even agree where the first four protons should be added, forming the central intermediate E4. We show that the prime reason for this is that different DFT methods give relative energies that differ by almost 600 kJ mol-1 for different protonation states. This is 4-30 times more than what is observed for other systems. The reason for this is that in some structures, the hydrogens bind to sulfide or carbide ions as protons, whereas in other structures they bind to the metals as hydride ions, changing the oxidation state of the metals, as well as the Fe-C, Fe-S and Fe-Fe distances. The energies correlate with the amount of Hartree-Fock exchange in the method, indicating a variation in the amount of static correlation in the structures. It is currently unclear which DFT method gives the best results for nitrogenase. We show that non-hybrid DFT functionals and TPSSh give the most accurate structures of the resting active site, whereas B3LYP and PBE0 give the best H2 dissociation energies. However, no DFT method indicates that a structure of E4 with two bridging hydride ions is lowest in energy, as spectroscopic experiments indicate.

The Mechanism of N-N Double Bond Cleavage by an Iron(II) Hydride Complex

The use of hydride species for substrate reductions avoids strong reductants, and may enable nitrogenase to reduce multiple bonds without unreasonably low redox potentials. In this work, we explore the N═N bond cleaving ability of a high-spin iron(II) hydride dimer with concomitant release of H2. Specifically, this diiron(II) complex reacts with azobenzene (PhN═NPh) to perform a four-electron reduction, where two electrons come from H2 reductive elimination and the other two come from iron oxidation. The rate law of the H2 releasing reaction indicates that diazene binding occurs prior to H2 elimination, and the negative entropy of activation and inverse kinetic isotope effect indicate that H-H bond formation is the rate-limiting step. Thus, substrate binding causes reductive elimination of H2 that formally reduces the metals, and the metals use the additional two electrons to cleave the N-N multiple bond.

Insight into the Iron-Molybdenum Cofactor of Nitrogenase from Synthetic Iron Complexes with Sulfur, Carbon, and Hydride Ligands

Nitrogenase enzymes are used by microorganisms for converting atmospheric N2 to ammonia, which provides an essential source of N atoms for higher organisms. The active site of the molybdenum-dependent nitrogenase is the unique carbide-containing iron-sulfur cluster called the iron-molybdenum cofactor (FeMoco). On the FeMoco, N2 binding is suggested to occur at one or more iron atoms, but the structures of the catalytic intermediates are not clear. In order to establish the feasibility of different potential mechanistic steps during biological N2 reduction, chemists have prepared iron complexes that mimic various structural aspects of the iron sites in the FeMoco. This reductionist approach gives mechanistic insight, and also uncovers fundamental principles that could be used more broadly for small-molecule activation. Here, we discuss recent results and highlight directions for future research. In one direction, synthetic iron complexes have now been shown to bind N2, break the N-N triple bond, and produce ammonia catalytically. Carbon- and sulfur-based donors have been incorporated into the ligand spheres of Fe-N2 complexes to show how these atoms may influence the structure and reactivity of the FeMoco. Hydrides have been incorporated into synthetic systems, which can bind N2, reduce some nitrogenase substrates, and/or reductively eliminate H2 to generate reduced iron centers. Though some carbide-containing iron clusters are known, none yet have sulfide bridges or high-spin iron atoms like the FeMoco.

Pentamethylcyclopentadienyl Half-Sandwich Diazoalkane Complexes of Ruthenium: Preparation and Reactivity

The diazoalkane complexes [Ru(η(5)-C5Me5)(N2CAr1Ar2){P(OR)3}L]BPh4 (1-4) [R = Me, L = P(OMe)3 (1); R = Et, L = P(OEt)3 (2); R = Me, L = PPh3 (3); R = Et, L = PPh3 (4); Ar1 = Ar2 = Ph (a); Ar1 = Ph, Ar2 = p-tolyl (b); Ar1Ar2 = C12H8 (c); Ar1 = Ph, Ar2 = PhC(O) (d)] and [Ru(η(5)-C5Me5){N2C(C12H8)}{PPh(OEt)2}(PPh3)]BPh4 (5c) were prepared by allowing chloro-compounds RuCl(η(5)-C5Me5)[P(OR)3]L to react with the diazoalkane Ar1Ar2CN2 in the presence of NaBPh4. Treatment of complexes 1-4 with H2O afforded 1,2-diazene derivatives [Ru(η(5)-C5Me5)(η(2)-NH═NH){P(OR)3}L]BPh4 (6-9) and ketone Ar1Ar2CO. A reaction path involving nucleophilic attack by H2O on the coordinated diazoalkane is proposed and supported by density functional theory calculations. The complexes were characterized spectroscopically (IR and (1)H, (31)P, (13)C, (15)N NMR) and by X-ray crystal structure determination of [Ru(η(5)-C5Me5)(N2CC12H8){P(OEt)3}2]BPh4 (2c) and [Ru(η(5)-C5Me5)(η(2)-NH═NH){P(OEt)3}2]BPh4 (7).

Facet-dependent solar ammonia synthesis of BiOCl nanosheets via a proton-assisted electron transfer pathway

Under the pressure of a fossil fuels shortage and global climate change, solar ammonia synthesis and the need to develop N2 fixation under mild conditions is becoming more urgent need; however, their intrinsic mechanisms still remain unclear. Herein, we demonstrate that the kinetic inertia of N2 can be overcome using oxygen vacancies (OVs) of BiOCl as the catalytic centers to create lower energy molecular steps, which are amendable for the solar light driven N-N triple bond cleavage via a proton-assisted electron transfer pathway. Moreover, the distinct structures of OVs on different BiOCl facets strongly determine the N2 fixation pathways by influencing both the adsorption structure and the activation level of N2. The fixation of terminal end-on bound N2 on the OVs of BiOCl {001} facets follows an asymmetric distal mode by selectively generating NH3, while the reduction of side-on bridging N2 on the OVs of BiOCl {010} facets is more energetically favorable in a symmetric alternating mode to produce N2H4 as the main intermediate.

Mechanistic insights into nitrogen fixation by nitrogenase enzymes

Biological nitrogen fixation by nitrogenase enzymes is a process that activates dinitrogen (N2) one of the most inert molecules in nature, within the confines of a living organism and at ambient conditions. Despite decades of study, there are still no complete explanations as to how this is possible. Here we describe a model of N2 reduction using the Mo-containing nitrogenase (FeMoco) that can explain the reactivity of the active site via a series of electrochemical steps that reversibly unseal a highly reactive Fe edge site. Our model can explain the 8 proton-electron transfers involved in biological ammonia synthesis within the kinetic scheme of Lowe and Thorneley, the obligatory formation of one H2 per N2 reduced, and the behavior of known inhibitors.

Hydrolysis of coordinated diazoalkanes to yield side-on 1,2-diazene derivatives

Diazoalkane complexes [Ru(η(5)-C5Me5)(N2CAr1Ar2)(PPh3){P(OR)3}]BPh4 [R = Me (1), Et (2); Ar1 = Ar2 = Ph (a); Ar1 = Ph, Ar2 = p-tolyl (b); Ar1Ar2 = C12H8 (c)] were prepared by allowing chloro complexes RuCl(η(5)-C5Me5)(PPh3)[P(OR)3] to react with diazoalkane Ar1Ar2CN2 in ethanol. The treatment of compounds 1 and 2 with H2O afforded 1,2-diazene derivatives [Ru(η(5)-C5Me5)(η(2)-NH═NH)(PPh3){P(OR)3}]BPh4 (3 and 4) and ketone Ar1Ar2CO. A reaction path involving nucleophilic attack by H2O on the coordinated diazoalkane is proposed. The complexes were characterized spectroscopically (IR and NMR) and by X-ray crystal structure determination of [Ru(η(5)-C5Me5)(η(2)-NH═NH)(PPh3){P(OMe)3}]BPh4 (3).

Misconception of reductive elimination of H2, in the context of the mechanism of nitrogenase

Calculated atom partial charges reveal misconceptions of reductive elimination of H₂.