Response to Comment on "Structural evidence for a dynamic metallocofactor during N2 reduction by Mo-nitrogenase"

Highly Activated Terminal Carbon Monoxide Ligand in an Iron-Sulfur Cluster Model of FeMco with Intermediate Local Spin State at Fe

Nitrogenases, the enzymes that convert N2 to NH3, also catalyze the reductive coupling of CO to yield hydrocarbons. CO-coordinated species of nitrogenase clusters have been isolated and used to infer mechanistic information. However, synthetic FeS clusters displaying CO ligands remain rare, which limits benchmarking. Starting from a synthetic cluster that models a cubane portion of the FeMo cofactor (FeMoco), including a bridging carbyne ligand, we report a heterometallic tungsten-iron-sulfur cluster with a single terminal CO coordination in two oxidation states with a high level of CO activation (νCO = 1851 and 1751 cm-1). The local Fe coordination environment (2S, 1C, 1CO) is identical to that in the protein making this system a suitable benchmark. Computational studies find an unusual intermediate spin electronic configuration at the Fe sites promoted by the presence the carbyne ligand. This electronic feature is partly responsible for the high degree of CO activation in the reduced cluster.

An orbitally adapted push-pull template for N2 activation and reduction to diazene-diide

A Lewis superacidic bis(borane) C6F4{B(C6F5)2}2 was reacted with tungsten N2-complexes [W(N2)2(R2PCH2CH2PR2)2] (R = Ph or Et), affording zwitterionic boryldiazenido W(ii) complexes trans-[W(L)(R2PCH2CH2PR2)2(N2{B(C6F5)2(C6F4B(C6F5)3})] (L = ø, N2 or THF). These compounds feature only one N-B linkage of the covalent type, as a result of intramolecular boron-to-boron C6F5 transfer. Complex trans-[W(THF)(Et2PCH2CH2PEt2)2(N2{B(C6F5)2C6F4B(C6F5)3})] (5) was shown to split H2, leading to a seven-coordinate complex [W(H)2(Et2PCH2CH2PEt2)2(N2{B(C6F5)2}2C6F4)] (7). Interestingly, hydride storage at the metal triggers backward C6F5 transfer. This reverts the bis(boron) moiety to its bis(borane) state, now doubly binding the distal N, with structural parameters and DFT computations pointing to dative N→B bonding. By comparison with an N2 complex [W(H)2(Et2PCH2CH2PEt2)2(N2{B(C6F5)3}] (10) differing only in the Lewis acid (LA), namely B(C6F5)3, coordinated to the distal N, we demonstrate that two-fold LA coordination imparts strong N2 activation up to the diazene-diide (N22-) state. To the best of our knowledge, this is the first example of a neutral LA coordination that induces reduction of N2.

Nitrogenase beyond the Resting State: A Structural Perspective

Nitrogenases have the remarkable ability to catalyze the reduction of dinitrogen to ammonia under physiological conditions. How does this happen? The current view of the nitrogenase mechanism focuses on the role of hydrides, the binding of dinitrogen in a reductive elimination process coupled to loss of dihydrogen, and the binding of substrates to a binuclear site on the active site cofactor. This review focuses on recent experimental characterizations of turnover relevant forms of the enzyme determined by cryo-electron microscopy and other approaches, and comparison of these forms to the resting state enzyme and the broader family of iron sulfur clusters. Emerging themes include the following: (i) The obligatory coupling of protein and electron transfers does not occur in synthetic and small-molecule iron-sulfur clusters. The coupling of these processes in nitrogenase suggests that they may involve unique features of the cofactor, such as hydride formation on the trigonal prismatic arrangement of irons, protonation of belt sulfurs, and/or protonation of the interstitial carbon. (ii) Both the active site cofactor and protein are dynamic under turnover conditions; the changes are such that more highly reduced forms may differ in key ways from the resting-state structure. Homocitrate appears to play a key role in coupling cofactor and protein dynamics. (iii) Structural asymmetries are observed in nitrogenase under turnover-relevant conditions by cryo-electron microscopy, although the mechanistic relevance of these states (such as half-of-sites reactivity) remains to be established.

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.

Sulfur-Ligated [2Fe-2C] Clusters as Synthetic Model Systems for Nitrogenase

Metal clusters featuring carbon and sulfur donors have coordination environments comparable to the active site of nitrogenase enzymes. Here, we report a series of di-iron clusters supported by the dianionic yldiide ligands, in which the Fe sites are bridged by two μ2-C atoms and four pendant S donors.The [L2Fe2] (L = {[Ph2P(S)]2C}2-) cluster is isolable in two oxidation levels, all-ferrous Fe2II and mixed-valence FeIIFeIII. The mixed-valence cluster displays two peaks in the Mössbauer spectra, indicating slow electron transfer between the two sites. The addition of the Lewis base 4-dimethylaminopyridine to the Fe2II cluster results in coordination with only one of the two Fe sites, even in the presence of an excess base. Conversely, the cluster reacts with 8 equiv of isocyanide tBuNC to give a monometallic complex featuring a new C-C bond between the ligand backbone and the isocyanide. The electronic structure descriptions of these complexes are further supported by X-ray absorption and resonant X-ray emission spectroscopies.

The HD Reaction of Nitrogenase: a Detailed Mechanism

Nitrogenase is the enzyme that converts N₂ to NH₃ under ambient conditions. The chemical mechanism of this catalysis at the active site FeMo-co [Fe₇ S₉ CMo(homocitrate)] is unknown. An obligatory co-product is H₂ , while exogenous H₂ is a competitive inhibitor. Isotopic substitution using exogenous D₂ revealed the N₂ -dependent reaction D₂ +2H⁺ +2e^- →2HD (the 'HD reaction'), together with a collection of additional experimental characteristics and requirements. This paper describes a detailed mechanism for the HD reaction, developed and elaborated using density functional simulations with a 486-atom model of the active site and surrounding protein. First D₂ binds at one Fe atom (endo-Fe6 coordination position), where it is flanked by H-Fe6 (exo position) and H-Fe2 (endo position). Then there is synchronous transfer of these two H atoms to bound D₂ , forming one HD bound to Fe2 and a second HD bound to Fe6. These two HD dissociate sequentially. The final phase is recovery of the two flanking H atoms. These H atoms are generated, sequentially, by translocation of a proton from the protein surface to S3B of FeMo-co and combination with introduced electrons. The first H atom migrates from S3B to exo-Fe6 and the second from S3B to endo-Fe2. Reaction energies and kinetic barriers are reported for all steps. This mechanism accounts for the experimental data: (a) stoichiometry; (b) the N₂ -dependence results from promotional N₂ bound at exo-Fe2; (c) different N₂ binding K_m for the HD reaction and the NH₃ formation reaction results from involvement of two different sites; (d) inhibition by CO; (e) the non-occurrence of 2HD→H₂ +D₂ results from the synchronicity of the two transfers of H to D₂ ; (f) inhibition of HD production at high pN₂ is by competitive binding of N₂ at endo-Fe6; (g) the non-leakage of D to solvent follows from the hydrophobic environment and irreversibility of proton introduction.

Understanding the tethered unhooking and rehooking of S2B in the reaction domain of FeMo-co, the active site of nitrogenase

The active site of the nitrogen fixing enzyme nitrogenase is an Fe₇MoS₉C cluster, and investigations of the enigmatic chemical mechanism of the enzyme have focussed on a pair of Fe atoms, Fe2 and Fe6, and the S2B atom that bridges them. There are three proposals for the status of the Fe2-S2B-Fe6 bridge during the catalytic cycle: one that it remains intact, another that it is completely labile and absent during catalysis, and a third that S2B is hemilabile, unhooking one of its bonds to Fe2 or Fe6. This report examines the tethered unhooking of S2B and factors that affect it, using DFT calculations of 50 geometric/electronic possibilities with a 485 atom model including all relevant parts of surrounding protein. The outcomes are: (a) unhooking the S2B-Fe2 bond is feasible and favourable, but alternative unhooking of the S2B-Fe6 bond is unlikely for steric reasons, (b) energy differences between hooked and unhooked isomers are generally <10 kcal mol^-1, usually with unhooked more stable, (c) ligation at the exo-Fe6 position inhibits unhooking, (d) unhooking of hydrogenated S2B is more favourable than that of bare S2B, (e) hydrogen bonding from the NεH function of His195 to S2B occurs in hooked and unhooked forms, and possibly stabilises unhooking, (f) unhooking is reversible with kinetic barriers ranging 10-13 kcal mol^-1. The conclusion is that energetically accessible reversible unhooking of S2B or S2BH, as an intrinsic property of FeMo-co, needs to be considered in the formulation of mechanisms for the reactions of nitrogenase.

Second and Outer Coordination Sphere Effects in Nitrogenase, Hydrogenase, Formate Dehydrogenase, and CO Dehydrogenase

Gases like H₂, N₂, CO₂, and CO are increasingly recognized as critical feedstock in "green" energy conversion and as sources of nitrogen and carbon for the agricultural and chemical sectors. However, the industrial transformation of N₂, CO₂, and CO and the production of H₂ require significant energy input, which renders processes like steam reforming and the Haber-Bosch reaction economically and environmentally unviable. Nature, on the other hand, performs similar tasks efficiently at ambient temperature and pressure, exploiting gas-processing metalloenzymes (GPMs) that bind low-valent metal cofactors based on iron, nickel, molybdenum, tungsten, and sulfur. Such systems are studied to understand the biocatalytic principles of gas conversion including N₂ fixation by nitrogenase and H₂ production by hydrogenase as well as CO₂ and CO conversion by formate dehydrogenase, carbon monoxide dehydrogenase, and nitrogenase. In this review, we emphasize the importance of the cofactor/protein interface, discussing how second and outer coordination sphere effects determine, modulate, and optimize the catalytic activity of GPMs. These may comprise ionic interactions in the second coordination sphere that shape the electron density distribution across the cofactor, hydrogen bonding changes, and allosteric effects. In the outer coordination sphere, proton transfer and electron transfer are discussed, alongside the role of hydrophobic substrate channels and protein structural changes. Combining the information gained from structural biology, enzyme kinetics, and various spectroscopic techniques, we aim toward a comprehensive understanding of catalysis beyond the first coordination sphere.

Radical SAM-dependent formation of a nitrogenase cofactor core on NifB

Nitrogenase is a versatile metalloenzyme that reduces N₂, CO and CO₂ at its cofactor site. Designated the M-cluster, this complex cofactor has a composition of [(R-homocitrate)MoFe₇S₉C], and it is assembled through the generation of a unique [Fe₈S₉C] core prior to the insertion of Mo and homocitrate. NifB is a radical S-adenosyl-L-methionine (SAM) enzyme that is essential for nitrogenase cofactor assembly. This review focuses on the recent work that sheds light on the role of NifB in the formation of the [Fe₈S₉C] core of the nitrogenase cofactor, highlighting the structure, function and mechanism of this unique radical SAM methyltransferase.

A Mechanism for Nitrogenase Including Loss of a Sulfide

Nitrogenase is the only enzyme in nature that can fix N₂ from the air. The active cofactor of the leading form of this enzyme contains seven irons and one molybdenum connected by sulfide bridges. In several recent experimental studies, it has been suggested that the cofactor is very flexible, and might lose one of its sulfides during catalysis. In this study, the possible loss of a sulfide has been investigated by model calculations. In previous studies, we have shown that there should be four activation steps before catalysis starts, and this study is based on that finding. It was found here that, after the four reductions in the activation steps, a sulfide will become very loosely bound and can be released in a quite exergonic step with a low barrier. The binding of N₂ has no part in that release. In our previous studies, we suggested that the central carbide should be protonated three times after the four activation steps. With the new finding, there will instead be a loss of a sulfide, as the barrier for the loss is much lower than the ones for protonating the carbide. Still, it is suggested here that the carbide will be protonated anyway, but only with one proton, in the E₃ to E₄ step. A very complicated transition state for H₂ formation involving a large structural change was obtained. The combined step, with a loss of H₂ and binding of N₂ , is calculated to be endergonic by +2.3 kcal mol^-1 ; this is in excellent agreement with experiments in which an easily reversible step has been found.

Comparisons of bond valences and distances for CO- and N2-bound clusters of FeMo-cofactors

By comparisons of N2 and isoelectronic substrate CO bound FeMo-cofactors (FeMo-cos) in nitrogenases, we have used a classical bond valence method to calculate the oxidation states of the iron and molybdenum atoms in FeMo-cos.

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

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

The E2 state of FeMoco: Hydride Formation versus Fe Reduction and a Mechanism for H2 Evolution

The iron‐molybdenum cofactor (FeMoco) is responsible for dinitrogen reduction in Mo nitrogenase. Unlike the resting state, E₀, reduced states of FeMoco are much less well characterized. The E₂ state has been proposed to contain a hydride but direct spectroscopic evidence is still lacking. The E₂ state can, however, relax back the E₀ state via a H₂ side‐reaction, implying a hydride intermediate prior to H₂ formation. This E₂→E₀ pathway is one of the primary mechanisms for H₂ formation under low‐electron flux conditions. In this study we present an exploration of the energy surface of the E₂ state. Utilizing both cluster‐continuum and QM/MM calculations, we explore various classes of E₂ models: including terminal hydrides, bridging hydrides with a closed or open sulfide‐bridge, as well as models without. Importantly, we find the hemilability of a protonated belt‐sulfide to strongly influence the stability of hydrides. Surprisingly, non‐hydride models are found to be almost equally favorable as hydride models. While the cluster‐continuum calculations suggest multiple possibilities, QM/MM suggests only two models as contenders for the E₂ state. These models feature either i) a bridging hydride between Fe₂ and Fe₆ and an open sulfide‐bridge with terminal SH on Fe₆ (E₂‐hyd) or ii) a double belt‐sulfide protonated, reduced cofactor without a hydride (E₂‐nonhyd). We suggest both models as contenders for the E₂ redox state and further calculate a mechanism for H₂ evolution. The changes in electronic structure of FeMoco during the proposed redox‐state cycle, E₀→E₁→E₂→E₀, are discussed.

Transition Metal Carbide Complexes

Carbide complexes remain a rare class of molecules. Their paucity does not reflect exceptional instability but is rather due to the generally narrow scope of synthetic procedures for constructing carbide complexes. The preparation of carbide complexes typically revolves around generating L_nM-CE_x fragments, followed by cleavage of the C-E bonds of the coordinated carbon-based ligands (the alternative being direct C atom transfer). Prime examples involve deoxygenation of carbonyl ligands and deprotonation of methyl ligands, but several other p-block fragments can be cleaved off to afford carbide ligands. This Review outlines synthetic strategies toward terminal carbide complexes, bridging carbide complexes, as well as carbide-carbonyl cluster complexes. It then surveys the reactivity of carbide complexes, covering stoichiometric reactions where the carbide ligands act as C₁ reagents, engage in cross-coupling reactions, and enact Fischer-Tropsch-like chemistry; in addition, we discuss carbide complexes in the context of catalysis. Finally, we examine spectroscopic features of carbide complexes, which helps to establish the presence of the carbide functionality and address its electronic structure.

The Conversion of Carbon Monoxide and Carbon Dioxide by Nitrogenases

Nitrogenases are the only known family of enzymes that catalyze the reduction of molecular nitrogen (N₂) to ammonia (NH₃). The N₂ reduction drives biological nitrogen fixation and the global nitrogen cycle. Besides the conversion of N₂, nitrogenases catalyze a whole range of other reductions, including the reduction of the small gaseous substrates carbon monoxide (CO) and carbon dioxide (CO₂) to hydrocarbons. However, it remains an open question whether these ‘side reactivities’ play a role under environmental conditions. Nonetheless, these reactivities and particularly the formation of hydrocarbons have spurred the interest in nitrogenases for biotechnological applications. There are three different isozymes of nitrogenase: the molybdenum and the alternative vanadium and iron‐only nitrogenase. The isozymes differ in their metal content, structure, and substrate‐dependent activity, despite their homology. This minireview focuses on the conversion of CO and CO₂ to methane and higher hydrocarbons and aims to specify the differences in activity between the three nitrogenase isozymes.

Uncovering Redox Non-innocent Hydrogen-Bonding in Cu(I)-Diazene Complexes

The life-sustaining reduction of N₂ to NH₃ is thermoneutral yet kinetically challenged by high-energy intermediates such as N₂H₂. Exploring intramolecular H-bonding as a potential strategy to stabilize diazene intermediates, we employ a series of [^xHetTpCu]₂(μ-N₂H₂) complexes that exhibit H-bonding between pendant aromatic N-heterocycles (^xHet) such as pyridine and a bridging trans-N₂H₂ ligand at copper(I) centers. X-ray crystallography and IR spectroscopy clearly reveal H-bonding in [^pyMeTpCu]₂(μ-N₂H₂) while low-temperature ¹H NMR studies coupled with DFT analysis reveals a dynamic equilibrium between two closely related, symmetric H-bonded structural motifs. Importantly, the ^xHet pendant negligibly influences the electronic structure of ^xHetTpCu^I centers in ^xHetTpCu(CNAr^2,6-Me₂) complexes that lack H-bonding as judged by nearly indistinguishable ν(CN) frequencies (2113-2117 cm^-1). Nonetheless, H-bonding in the corresponding [^xHetTpCu]₂(μ-N₂H₂) complexes results in marked changes in ν(NN) (1398-1419 cm^-1) revealed through resonance Raman studies. Due to the closely matched N-H BDEs of N₂H₂ and the pyH⁰ cation radical, the aromatic N-heterocyclic pendants may encourage partial H-atom transfer (HAT) from N₂H₂ to ^xHet through redox-non-innocent H-bonding in [^xHetTpCu]₂(μ-N₂H₂). DFT studies reveal modest thermodynamic barriers for concerted transfer of both H-atoms of coordinated N₂H₂ to the ^xHet pendants to generate tautomeric [^xHetHTpCu]₂(μ-N₂) complexes, identifying metal-assisted concerted dual HAT as a thermodynamically favorable pathway for N₂/N₂H₂ interconversion.

Looking for Options to Sustainably Fixate Nitrogen. Are Molecular Metal Oxides Catalysts a Viable Avenue?

Fast and reliable industrial production of ammonia (NH₃) is fundamentally sustaining modern society. Since the early 20^th Century, NH₃ has been synthesized via the Haber-Bosch process, running at conditions of around 350-500°C and 100-200 times atmospheric pressure (15-20 MPa). Industrial ammonia production is currently the most energy-demanding chemical process worldwide and contributes up to 3% to the global carbon dioxide emissions. Therefore, the development of more energy-efficient pathways for ammonia production is an attractive proposition. Over the past 20 years, scientists have imagined the possibility of developing a milder synthesis of ammonia by mimicking the nitrogenase enzyme, which fixes nitrogen from the air at ambient temperatures and pressures to feed leguminous plants. To do this, we propose the use of highly reconfigurable molecular metal oxides or polyoxometalates (POMs). Our proposal is an informed design of the polyoxometalate after exploring the catabolic pathways that cyanobacteria use to fix N₂ in nature, which are a different route than the one followed by the Haber-Bosch process. Meanwhile, the industrial process is a "brute force" system towards breaking the triple bond N-N, needing high pressure and high temperature to increase the rate of reaction, nature first links the protons to the N₂ to later easier breaking of the triple bond at environmental temperature and pressure. Computational chemistry data on the stability of different polyoxometalates will guide us to decide the best design for a catalyst. Testing different functionalized molecular metal oxides as ammonia catalysts laboratory conditions will allow for a sustainable reactor design of small-scale production.