Comparative assessment of the composition and charge state of nitrogenase FeMo-cofactor

Controllable Construction of a Mo2C/MoO2 Interface with an Ideal Mo2C/MoO2 Ratio for Efficient Electrocatalytic Nitrogen Reduction to Ammonia

Electrocatalytic nitrogen reduction reaction (NRR) is considered to be a viable contender for the production of NH3. However, due to the sluggish adsorption and activation of the electrocatalyst toward inert N2 molecules, there is an urgent need for developing effective catalysts to facilitate the reaction. Inspired by natural nitrogenase, in which Mo atoms are the active centers, Mo-based electrocatalysts have received considerable attention, but further exploration is still necessary. Interface-engineered electrocatalysts can effectively optimize the absorption and activation of the catalytic active center for N2 and thus improve the electrocatalytic activity of NRR. However, the lack of studies for controllably constructing an optimal ratio of two phases at the interface hinders the development of NRR electrocatalysts. Herein, a series of Mo2C/MoO2 interface-engineered electrocatalysts with various Mo2C/MoO2 ratios were constructed by controlling the Y dosages. The controlled experimental results verified that the catalytic activity of NRR, the dosage of Y, and the ratio of Mo2C/MoO2 were strongly correlated. Density functional theory calculations show that the C-Mo-O coordination at the Mo2C/MoO2 interface can optimize the reaction path and reduce the energy barrier of the reaction intermediates, thereby enhancing the reaction kinetics of NRR.

Deciphering Iron-Catalyzed C-H Amination with Organic Azides: N2 Cleavage from a Stable Organoazide Complex

Catalytic C-N bond formation by direct activation of C-H bonds offers wide synthetic potential. En route to C-H amination, complexes with organic azides are critical precursors towards the reactive nitrene intermediate. Despite their relevance, α-N coordinated organoazide complexes are scarce in general, and elusive with iron, although iron complexes are by far the most active catalysts for C-H amination with organoazides. Herein, we report the synthesis of a stable iron α-N coordinated organoazide complex from [Fe(N(SiMe3 )2 )2 ] and AdN3 (Ad=1-adamantyl) and its crystallographic, IR, NMR and zero-field 57 Fe Mössbauer spectroscopic characterization. These analyses revealed that the organoazide is in fast equilibrium between the free and coordinated state (Keq =62). Photo-crystallography experiments showed gradual dissociation of N2 , which imparted an Fe-N bond shortening and correspond to structural snapshots of the formation of an iron imido/nitrene complex. Reactivity of the organoazide complex in solution showed complete loss of N2 , and subsequent formation of a C-H aminated product via nitrene insertion into a C-H bond of the N(SiMe3 )2 ligand. Monitoring this reaction by 1 H NMR spectroscopy indicates the transient formation of the imido/nitrene intermediate, which was supported by Mössbauer spectroscopy in frozen solution.

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.

Bioinspired Framework Catalysts: From Enzyme Immobilization to Biomimetic Catalysis

Enzymatic catalysis has fueled considerable interest from chemists due to its high efficiency and selectivity. However, the structural complexity and vulnerability hamper the application potentials of enzymes. Driven by the practical demand for chemical conversion, there is a long-sought quest for bioinspired catalysts reproducing and even surpassing the functions of natural enzymes. As nanoporous materials with high surface areas and crystallinity, metal-organic frameworks (MOFs) represent an exquisite case of how natural enzymes and their active sites are integrated into porous solids, affording bioinspired heterogeneous catalysts with superior stability and customizable structures. In this review, we comprehensively summarize the advances of bioinspired MOFs for catalysis, discuss the design principle of various MOF-based catalysts, such as MOF-enzyme composites and MOFs embedded with active sites, and explore the utility of these catalysts in different reactions. The advantages of MOFs as enzyme mimetics are also highlighted, including confinement, templating effects, and functionality, in comparison with homogeneous supramolecular catalysts. A perspective is provided to discuss potential solutions addressing current challenges in MOF catalysis.

Connecting the geometric and electronic structures of the nitrogenase iron-molybdenum cofactor through site-selective 57Fe labelling

Understanding the chemical bonding in the catalytic cofactor of the Mo nitrogenase (FeMo-co) is foundational for building a mechanistic picture of biological nitrogen fixation. A persistent obstacle towards this goal has been that the 57Fe-based spectroscopic data-although rich with information-combines responses from all seven Fe sites, and it has therefore not been possible to map individual spectroscopic responses to specific sites in the three-dimensional structure. Here we have addressed this challenge by incorporating 57Fe into a single site of FeMo-co. Spectroscopic analysis of the resting state informed on the local electronic structure of the terminal Fe1 site, including its oxidation state and spin orientation, and, in turn, on the spin-coupling scheme for the entire cluster. The oxidized resting state and the first intermediate in nitrogen fixation were also characterized, and comparisons with the resting state provided molecular-level insights into the redox chemistry of FeMo-co.

Two Shared Icosahedral Metallacarboranes through Iron: A Joint Experimental and Theoretical Refinement of Mössbauer Spectrum in [Fe(1,2-C2B9H11)2]Cs

Mössbauer and X-ray photoelectron spectroscopies (XPS) are complemented with high-level quantum-chemical computations in the study of the geometric and electronic structure of the paramagnetic salt of the metallacarborane sandwich complex [Fe(1,2-C₂B₉H₁₁)₂]Cs = FeSanCs. Experimental ⁵⁷Fe isomer shifts and quadrupole splitting parameters are compared with the theoretical prediction, with good agreement. The appearance of two sets of Cs(3d) doublets in the XPS spectrum, separated by 2 eV, indicates that Cs has two different chemical environments due to ease of the Cs⁽⁺⁾ cation moving around the sandwich complex with low-energy barriers, as confirmed by quantum-chemical computations. Several minimum-energy geometries of the FeSanCs structure with the corresponding energies and Mössbauer parameters are discussed, in particular the atomic charges and spin population and the surroundings of the Fe atom in the complex. The Mössbauer spectra were taken at different temperatures showing the presence of a low-spin Fe atom with S = 1/2 and thus confirming a paramagnetic Fe^III species.

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

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

Carbon Monoxide Binding to the Iron-Molybdenum Cofactor of Nitrogenase: a Detailed Quantum Mechanics/Molecular Mechanics Investigation

Carbon monoxide (CO) is a well-known inhibitor of nitrogenase activity. Under turnover conditions, CO binds to FeMoco, the active site of Mo nitrogenase. Time-resolved IR measurements suggest an initial terminal CO at 1904 cm^-1 that converts to a bridging CO at 1715 cm^-1, and an X-ray structure shows that CO can displace one of the bridging belt sulfides of FeMoco. However, the CO-binding redox state(s) of FeMoco (E_n) and the role of the protein environment in stabilizing specific CO-bound intermediates remain elusive. In this work, we carry out an in-depth analysis of the CO-FeMoco interaction based on quantum chemical calculations addressing different aspects of the electronic structure. (1) The local electronic structure of the Fe-CO bond is studied through diamagnetically substituted FeMoco. (2) A cluster model of FeMoco within a polarizable continuum illustrates how CO binding may affect the spin-coupling between the metal centers. (3) A QM/MM model incorporates the explicit influence of the amino acid residues surrounding FeMoco in the MoFe protein. The QM/MM model predicts both a terminal and a bridging CO in the E₁ redox state. The scaled calculated CO frequencies (1922 and 1716 cm^-1, respectively) are in good agreement with the experimentally observed IR bands supporting CO binding to the E₁ state. Alternatively, an E₂ state QM/MM model, which has the same atomic structure as the CO-bound X-ray structure, features a semi-bridging CO with a scaled calculated frequency (1718 cm^-1) similar to the bridging CO in the E₁ model.

Structure, reactivity, and spectroscopy of nitrogenase-related synthetic and biological clusters

The reduction of dinitrogen (N2) is essential for its incorporation into nucleic acids and amino acids, which are vital to life on earth. Nitrogenases convert atmospheric dinitrogen to two ammonia molecules (NH3) under ambient conditions. The catalytic active sites of these enzymes (known as FeM-cofactor clusters, where M = Mo, V, Fe) are the sites of N2 binding and activation and have been a source of great interest for chemists for decades. In this review, recent studies on nitrogenase-related synthetic molecular complexes and biological clusters are discussed, with a focus on their reactivity and spectroscopic characterization. The molecular models that are discussed span from simple mononuclear iron complexes to multinuclear iron complexes and heterometallic iron complexes. In addition, recent work on the extracted biological cofactors is discussed. An emphasis is placed on how these studies have contributed towards our understanding of the electronic structure and mechanism of nitrogenases.

Mössbauer Spectroscopic and Computational Investigation of An Iron Cyclopentadienone Complex

(Cyclopentadienone)iron carbonyl complexes have recently received particular attention for their use as catalysts in hydrogenation or transfer hydrogenation reactions including the N-alkylation of amines with alcohols. This is due to their easy synthesis from simple and cheap materials, air and water stabilities, and the crucial metal-ligand cooperation giving rise to unique catalytic properties. Here, we report a Mössbauer spectroscopic and computational investigation of such a complex and its corresponding activated species for dehydrogenation and hydrogenation reactions. This study affords a deeper understanding of the species formed by the reaction with Me₃NO and their distribution upon the added amount of an oxidant.

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

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

Postbiosynthetic modification of a precursor to the nitrogenase iron-molybdenum cofactor

Nitrogenases utilize Fe-S clusters to reduce N₂ to NH₃ The large number of Fe sites in their catalytic cofactors has hampered spectroscopic investigations into their electronic structures, mechanisms, and biosyntheses. To facilitate their spectroscopic analysis, we are developing methods for incorporating ⁵⁷Fe into specific sites of nitrogenase cofactors, and we report herein site-selective ⁵⁷Fe labeling of the L-cluster-a carbide-containing, [Fe₈S₉C] precursor to the Mo nitrogenase catalytic cofactor. Treatment of the isolated L-cluster with the chelator ethylenediaminetetraacetate followed by reconstitution with ⁵⁷Fe²⁺ results in ⁵⁷Fe labeling of the terminal Fe sites in high yield and with high selectivity. This protocol enables the generation of L-cluster samples in which either the two terminal or the six belt Fe sites are selectively labeled with ⁵⁷Fe. Mössbauer spectroscopic analysis of these samples bound to the nitrogenase maturase Azotobacter vinelandii NifX reveals differences in the primary coordination sphere of the terminal Fe sites and that one of the terminal sites of the L-cluster binds to H35 of Av NifX. This work provides molecular-level insights into the electronic structure and biosynthesis of the L-cluster and introduces postbiosynthetic modification as a promising strategy for studies of nitrogenase cofactors.

Assignment of protonated R-homocitrate in extracted FeMo-cofactor of nitrogenase via vibrational circular dichroism spectroscopies

Protonation of FeMo-cofactor is important for the process of substrate hydrogenation. Its structure has been clarified as Δ-Mo*Fe₇S₉C(R-homocit*)(cys)(Hhis) for the efforts of nearly 30 years, while it remains controversial whether FeMo-cofactor is protonated or deprotonated with chelated ≡C-O(H) homocitrate. We have used protonated molybdenum(V) lactates 1 and its enantiomer as model compounds for R-homocitrate in FeMo-cofactor of nitrogenase. Vibrational circular dichroism (VCD) spectrum of 1 at 1051 cm^-1 is attributed to ≡C-O_H vibration, and molybdenum(VI) R-lactate at 1086 cm^-1 is assigned as ≡C-O _α-alkoxy vibration. These vibrations set up labels for the protonation state of coordinated α-hydroxycarboxylates. The characteristic VCD band of NMF-extracted FeMo-cofactor is assigned to ν(C-O_H), which is based on the comparison of molybdenum(VI) R-homocitrate. Density Functional Theory calculations are in consistent with these assignments. To the best of our knowledge, this is the first time that protonated R-homocitrate in FeMo-cofactor is confirmed by VCD spectra.

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.

Structural evidence for a dynamic metallocofactor during N2 reduction by Mo-nitrogenase

The enzyme nitrogenase uses a suite of complex metallocofactors to reduce dinitrogen (N2) to ammonia. Mechanistic details of this reaction remain sparse. We report a 1.83-angstrom crystal structure of the nitrogenase molybdenum-iron (MoFe) protein captured under physiological N2 turnover conditions. This structure reveals asymmetric displacements of the cofactor belt sulfurs (S2B or S3A and S5A) with distinct dinitrogen species in the two αβ dimers of the protein. The sulfur-displaced sites are distinct in the ability of protein ligands to donate protons to the bound dinitrogen species, as well as the elongation of either the Mo-O5 (carboxyl) or Mo-O7 (hydroxyl) distance that switches the Mo-homocitrate ligation from bidentate to monodentate. These results highlight the dynamic nature of the cofactor during catalysis and provide evidence for participation of all belt-sulfur sites in this process.

Bond-valence analyses of the crystal structures of FeMo/V cofactors in FeMo/V proteins

The bond-valence method has been used for valence calculations of FeMo/V cofactors in FeMo/V proteins using 51 crystallographic data sets of FeMo/V proteins from the Protein Data Bank. The calculations show molybdenum(III) to be present in MoFe₇S₉C(Cys)(HHis)[R-(H)homocit] (where H₄homocit is homocitric acid, HCys is cysteine and HHis is histidine) in FeMo cofactors, while vanadium(III) with a more reduced iron complement is obtained for FeV cofactors. Using an error analysis of the calculated valences, it was found that in FeMo cofactors Fe1, Fe6 and Fe7 can be unambiguously assigned as iron(III), while Fe2, Fe3, Fe4 and Fe5 show different degrees of mixed valences for the individual Fe atoms. For the FeV cofactors in PDB entry 5n6y, Fe4, Fe5 and Fe6 correspond to iron(II), iron(II) and iron(III), respectively, while Fe1, Fe2, Fe3 and Fe7 exhibit strongly mixed valences. Special situations such as CO-bound and selenium-substituted FeMo cofactors and O(N)H-bridged FeV cofactors are also discussed and suggest rearrangement of the electron configuration on the substitution of the bridging S atoms.

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.

A model for dinitrogen binding in the E4 state of nitrogenase

Molybdenum nitrogenase is one of the most intriguing metalloenzymes in nature, featuring an exotic iron-molybdenum-sulfur cofactor, FeMoco, whose mode of action remains elusive. In particular, the molecular and electronic structure of the N2-binding E4 state is not known. In this study we present theoretical QM/MM calculations of new structural models of the E4 state of molybdenum-dependent nitrogenase and compare to previously suggested models for this enigmatic redox state. We propose two models as possible candidates for the E4 state. Both models feature two hydrides on the FeMo cofactor, bridging atoms Fe2 and Fe6 with a terminal sulfhydryl group on either Fe2 or Fe6 (derived from the S2B bridge) and the change in coordination results in local lower-spin electronic structure at Fe2 and Fe6. These structures appear consistent with the bridging hydride proposal put forward from ENDOR studies and are calculated to be lower in energy than other proposed models for E4 at the TPSSh-QM/MM level of theory. We critically analyze the DFT method dependency in calculations of FeMoco that has resulted in strikingly different proposals for this state. Importantly, dinitrogen binds exothermically to either Fe2 or Fe6 in our models, contrary to others, an effect rationalized via the unique ligand field (from the hydrides) at the Fe with an empty coordination site. A low-spin Fe site is proposed as being important to N2 binding. Furthermore, the geometries of these states suggest a feasible reductive elimination step that could follow, as experiments indicate. Via this step, two electrons are released, reducing the cofactor to yield a distorted 4-coordinate Fe2 or Fe6 that partially activates N2. We speculate that stabilization of an N2-bound Fe(i) at Fe6 (not found for Fe2 model) via reductive elimination is a crucial part of N2 activation in nitrogenases, possibly aided by the apical heterometal ion (Mo or V). By using protons from the sulfhydryl group (to regenerate the sulfide bridge between Fe2 and Fe6) and the nearby homocitrate hydroxy group, we calculate a plausible route to yield a diazene intermediate. This is found to be more favorable with the Fe6-bound model than the Fe2-bound model; however, this protonation is uphill in energy, suggesting protonation of N2 might occur later in the catalytic cycle or via another mechanism.

Preliminary Assignment of Protonated and Deprotonated Homocitrates in Extracted FeMo-Cofactors by Comparisons with Molybdenum(IV) Lactates and Oxidovanadium Glycolates

A similar pair of protonated and deprotonated mononuclear oxidovanadium glycolates [VO(Hglyc)(phen)(H₂O)]Cl·2H₂O (1) and [VO(glyc)(bpy)(H₂O)] (2) and a mixed-(de)protonated oxidovanadium triglycolate (NH₄)₂[VO(Hglyc)₂(glyc)]·H₂O (3) were isolated and examined. The ≡C-O(H) (≡C-OH or ≡C-O) groups coordinated to vanadium were spectroscopically and structurally identified. The glycolate in 1 features a bidentate chelation through protonated α-hydroxy and α-carboxy groups, whereas the glycolate in 2 coordinates through deprotonated α-alkoxy and α-carboxy groups. The glycolates in 3 coordinate to vanadium through α-alkoxy or α-hydroxy and α-carboxy groups and thus have both protonated ≡C-OH and deprotonated ≡C-O bonds simultaneously. Structural investigations revealed that the longer protonated V-O_α-hydroxy bonds [2.234(2) Å and 2.244(2) Å] in 1 and 3 are close to those of FeV-cofactor (FeV-co) 2.17 Ź (FeMo-co 2.17 Ų), while deprotonated V-O_α-alkoxy bonds [2, 1.930(2); 3, 1.927(2) Å] were obviously shorter. This shows a similar elongated trend as the Mo-O distances in the previously reported deprotonated vs protonated molybdenum lactates (Wang, S. Y. et al. Dalton Trans. 2018, 47, 7412-7421) and these vanadium and molybdenum complexes have the same local V/Mo-homocitrate structures as those of FeV/Mo-cos of nitrogenases. The IR spectra of these oxidovanadium and the previously synthesized molybdenum complexes including different substituted ≡C-O(H) model compounds show red-shifts for ≡C-OH vs ≡C-O alternation, which further assign the two IR bands of extracted FeMo-co at 1084 and 1031 cm^-1 to ≡C-O and ≡C-OH vibrations, respectively. Although the structural data or IR spectra for some of the previously synthesized Mo/V complexes and extracted FeMo-co were measured earlier, this is the first time that the ≡C-O(H) coordinated peaks are assigned. The overall structural and IR results well suggest the coexistence of homocitrates coordinated with α-alkoxy (deprotonated) and α-hydroxy (protonated) groups in the extracted FeMo-co.

Critical computational analysis illuminates the reductive-elimination mechanism that activates nitrogenase for N2 reduction

Recent spectroscopic, kinetic, photophysical, and thermodynamic measurements show activation of nitrogenase for N₂ → 2NH₃ reduction involves the reductive elimination (re) of H₂ from two [Fe-H-Fe] bridging hydrides bound to the catalytic [7Fe-9S-Mo-C-homocitrate] FeMo-cofactor (FeMo-co). These studies rationalize the Lowe-Thorneley kinetic scheme's proposal of mechanistically obligatory formation of one H₂ for each N₂ reduced. They also provide an overall framework for understanding the mechanism of nitrogen fixation by nitrogenase. However, they directly pose fundamental questions addressed computationally here. We here report an extensive computational investigation of the structure and energetics of possible nitrogenase intermediates using structural models for the active site with a broad range in complexity, while evaluating a diverse set of density functional theory flavors. (i) This shows that to prevent spurious disruption of FeMo-co having accumulated 4[e ^-/H⁺] it is necessary to include: all residues (and water molecules) interacting directly with FeMo-co via specific H-bond interactions; nonspecific local electrostatic interactions; and steric confinement. (ii) These calculations indicate an important role of sulfide hemilability in the overall conversion of E ₀ to a diazene-level intermediate. (iii) Perhaps most importantly, they explain (iiia) how the enzyme mechanistically couples exothermic H₂ formation to endothermic cleavage of the N≡N triple bond in a nearly thermoneutral re/oxidative-addition equilibrium, (iiib) while preventing the "futile" generation of two H₂ without N₂ reduction: hydride re generates an H₂ complex, but H₂ is only lost when displaced by N₂, to form an end-on N₂ complex that proceeds to a diazene-level intermediate.

Reliable Estimation of Prediction Uncertainty for Physicochemical Property Models

One of the major challenges in computational science is to determine the uncertainty of a virtual measurement, that is the prediction of an observable based on calculations. As highly accurate first-principles calculations are in general unfeasible for most physical systems, one usually resorts to parameteric property models of observables, which require calibration by incorporating reference data. The resulting predictions and their uncertainties are sensitive to systematic errors such as inconsistent reference data, parametric model assumptions, or inadequate computational methods. Here, we discuss the calibration of property models in the light of bootstrapping, a sampling method that can be employed for identifying systematic errors and for reliable estimation of the prediction uncertainty. We apply bootstrapping to assess a linear property model linking the ⁵⁷Fe Mössbauer isomer shift to the contact electron density at the iron nucleus for a diverse set of 44 molecular iron compounds. The contact electron density is calculated with 12 density functionals across Jacob's ladder (PWLDA, BP86, BLYP, PW91, PBE, M06-L, TPSS, B3LYP, B3PW91, PBE0, M06, TPSSh). We provide systematic-error diagnostics and reliable, locally resolved uncertainties for isomer-shift predictions. Pure and hybrid density functionals yield average prediction uncertainties of 0.06-0.08 mm s^-1 and 0.04-0.05 mm s^-1, respectively, the latter being close to the average experimental uncertainty of 0.02 mm s^-1. Furthermore, we show that both model parameters and prediction uncertainty depend significantly on the composition and number of reference data points. Accordingly, we suggest that rankings of density functionals based on performance measures (e.g., the squared coefficient of correlation, r², or the root-mean-square error, RMSE) should not be inferred from a single data set. This study presents the first statistically rigorous calibration analysis for theoretical Mössbauer spectroscopy, which is of general applicability for physicochemical property models and not restricted to isomer-shift predictions. We provide the statistically meaningful reference data set MIS39 and a new calibration of the isomer shift based on the PBE0 functional.

Reactivity of hydride bridges in a high-spin [Fe3(μ-H)3]3+ cluster: reversible H2/CO exchange and Fe-H/B-F bond metathesis

The triiron trihydride complex Fe₃H₃L (1) [where L^3– is a tris(β-diketiminate)cyclophanate] reacts with CO and with BF₃·OEt₂ to afford (Fe^ICO)₂Fe^II(μ₃-H)L (2) and Fe₃F₃L (3), respectively.

The triiron trihydride complex Fe₃H₃L (1) [where L^3– is a tris(β-diketiminate)cyclophanate] reacts with CO and with BF₃·OEt₂ to afford (Fe^ICO)₂Fe^II(μ₃-H)L (2) and Fe₃F₃L (3), respectively. Variable-temperature and applied-field Mössbauer spectroscopy support the assignment of two high-spin (HS) iron(i) centers and one HS iron(ii) ion in 2. Preliminary studies support a CO-induced reductive elimination of H₂ from 1, rather than CO trapping a species from an equilibrium mixture. This complex reacts with H₂ to regenerate 1 under a dihydrogen atmosphere, which represents a rare example of reversible CO/H₂ exchange and the first to occur at high-spin metal centers, as well as the first example of a reversible multielectron redox reaction at a designed high-spin metal cluster. The formation of 3 proceeds through a previously unreported net fluoride-for-hydride substitution, and 3 is surprisingly chemically inert to Si–H bonds and points to an unexpectedly large difference between the Fe–F and Fe–H bonds in this high-spin system.

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.

X-ray Absorption and Emission Spectroscopic Studies of [L2Fe2S2](n) Model Complexes: Implications for the Experimental Evaluation of Redox States in Iron-Sulfur Clusters

Herein, a systematic study of [L2Fe2S2](n) model complexes (where L = bis(benzimidazolato) and n = 2-, 3-, 4-) has been carried out using iron and sulfur K-edge X-ray absorption (XAS) and iron Kβ and valence-to-core X-ray emission spectroscopies (XES). These data are used as a test set to evaluate the relative strengths and weaknesses of X-ray core level spectroscopies in assessing redox changes in iron-sulfur clusters. The results are correlated to density functional theory (DFT) calculations of the spectra in order to further support the quantitative information that can be extracted from the experimental data. It is demonstrated that due to canceling effects of covalency and spin state, the information that can be extracted from Fe Kβ XES mainlines is limited. However, a careful analysis of the Fe K-edge XAS data shows that localized valence vs delocalized valence species may be differentiated on the basis of the pre-edge and K-edge energies. These findings are then applied to existing literature Fe K-edge XAS data on the iron protein, P-cluster, and FeMoco sites of nitrogenase. The ability to assess the extent of delocalization in the iron protein vs the P-cluster is highlighted. In addition, possible charge states for FeMoco on the basis of Fe K-edge XAS data are discussed. This study provides an important reference for future X-ray spectroscopic studies of iron-sulfur clusters.

Nitrogenase FeMoco investigated by spatially resolved anomalous dispersion refinement

The [Mo:7Fe:9S:C] iron-molybdenum cofactor (FeMoco) of nitrogenase is the largest known metal cluster and catalyses the 6-electron reduction of dinitrogen to ammonium in biological nitrogen fixation. Only recently its atomic structure was clarified, while its reactivity and electronic structure remain under debate. Here we show that for its resting S=3/2 state the common iron oxidation state assignments must be reconsidered. By a spatially resolved refinement of the anomalous scattering contributions of the 7 Fe atoms of FeMoco, we conclude that three irons (Fe1/3/7) are more reduced than the other four (Fe2/4/5/6). Our data are in agreement with the recently revised oxidation state assignment for the molybdenum ion, providing the first spatially resolved picture of the resting-state electron distribution within FeMoco. This might provide the long-sought experimental basis for a generally accepted theoretical description of the cluster that is in line with available spectroscopic and functional data.

The [Mo:7Fe:9S:C] iron-molybdenum cofactor (FeMoco) of nitrogenase is a large metal cluster with an important role in biological nitrogen fixation. Here, the authors use spatially resolved refinement of the anomalous scattering contributions of the iron atoms to determine the resting-state electron distribution of FeMoco.

Combining a Nitrogenase Scaffold and a Synthetic Compound into an Artificial Enzyme

Nitrogenase catalyzes substrate reduction at its cofactor center ([(Cit)MoFe7S9C](n-); designated M-cluster). Here, we report the formation of an artificial, nitrogenase-mimicking enzyme upon insertion of a synthetic model complex ([Fe6S9(SEt)2](4-); designated Fe6(RHH)) into the catalytic component of nitrogenase (designated NifDK(apo)). Two Fe6(RHH) clusters were inserted into NifDK(apo), rendering the conformation of the resultant protein (designated NifDK(Fe)) similar to the one upon insertion of native M-clusters. NifDK(Fe) can work together with the reductase component of nitrogenase to reduce C2H2 in an ATP-dependent reaction. It can also act as an enzyme on its own in the presence of Eu(II)DTPA, displaying a strong activity in C2H2 reduction while demonstrating an ability to reduce CN(-) to C1-C3 hydrocarbons in an ATP-independent manner. The successful outcome of this work provides the proof of concept and underlying principles for continued search of novel enzymatic activities based on this approach.

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.

The discovery of Mo(III) in FeMoco: reuniting enzyme and model chemistry

Biological nitrogen fixation is enabled by molybdenum-dependent nitrogenase enzymes, which effect the reduction of dinitrogen to ammonia using an Fe₇MoS₉C active site, referred to as the iron molybdenum cofactor or FeMoco. In this mini-review, we summarize the current understanding of the molecular and electronic structure of FeMoco. The advances in our understanding of the active site structure are placed in context with the parallel evolution of synthetic model studies. The recent discovery of Mo(III) in the FeMoco active site is highlighted with an emphasis placed on the important role that model studies have played in this finding. In addition, the reactivities of synthetic models are discussed in terms of their relevance to the enzymatic system.

Valence to core X-ray emission spectroscopy

This Progress Report discusses the chemical sensitivity of Kβ valence to core X-ray emission spectroscopy (vtc-XES) and its applications for investigating 3d-transition-metal based materials. Vtc-XES can be used for ligand identification and for the characterization of the valence electronic levels. The technique provides information that is similar to valence band photoemission spectroscopy but the sample environment can be chosen freely and thus allows measurements in presence of gases and liquids and it can be applied for measurements under in situ/operando or extreme conditions. The theoretical basis of the technique is presented using a one-electron approach and the vtc-XES spectral features are interpreted using ground state density functional theory calculations. Some recent results obtained by vtc-XES in various scientific fields are discussed to demonstrate the potential and future applications of this technique. Resonant X-ray emission spectroscopy is briefly introduced with some applications for the study of 3d and 5d-transition-metal based systems.

A journey into the active center of nitrogenase

Nitrogenase catalyzes the reduction of N2 to NH3, a key step in the global nitrogen cycle. This article describes our journey toward the definition of a complete molecular structure of the active site of nitrogenase, with an emphasis on the discovery of the interstitial carbide and the radical SAM-dependent insertion of this atom into the active FeMo cofactor site of nitrogenase.

Identification of a spin-coupled Mo(iii) in the nitrogenase iron–molybdenum cofactor

The molybdenum atom in FeMoco is imperative to the high activity of the enzyme and has been proposed to be Mo(iv). We demonstrate that only Mo(iii) fits Mo HERFD XAS data, the first example of Mo(iii) in biology. Theoretical calculations further reveal an unusual spin-coupled Mo(iii).

The positions of inner hydroxide groups and aluminium ions in exfoliated kaolinite as indicators of the external chemical environment

New insights: a molecular cluster model was created for exfoliated kaolinite using coordination chemistry principles highlighting the remarkable structural differences relative to crystalline kaolinite.