Nitrogenase FeMoco investigated by spatially resolved anomalous dispersion refinement

Composition Determination of Heterometallic Trinuclear Clusters via Anomalous X-ray and Neutron Diffraction

Anomalous X-ray diffraction (AXD) and neutron diffraction can be used to crystallographically distinguish between metals of similar electron density. Despite the use of AXD for structural characterization in mixed metal clusters, there are no benchmark studies evaluating the accuracy of AXD toward assessing elemental occupancy in molecules with comparisons with what is determined via neutron diffraction. We collected resonant diffraction data on several homo and heterometallic clusters and refined their anomalous scattering components to determine metal site occupancies. Theoretical resonant scattering terms for Fe0, Co0, and Zn0 were compared against experimental values, revealing theoretical values are ill-suited to serve as references for occupancy determination. The cluster featuring distinct cation and anion metal compositions [CoCp2*][(tbsL)Fe3(μ3-NAr)] was used to assess the accuracy of different f' references for occupancy determination (f'theoretical ± 15-17%; f'experimental ± 10%). This methodology was applied toward calculating the occupancy of three different clusters: (tbsL)Fe2Zn(py) (6), (tbsL)Fe2Zn(μ3-NAr)(py) (7), and [CoCp*2][(tbsL)Fe2Zn(μ3-NAr)] (8). The first two clusters maintain 100% Fe/Zn site isolation, whereas 8 showed metal mixing within the sites. The large crystal size of 8 enabled collection of neutron diffraction data which was compared against the results found with AXD. The ability of AXD to replicate the metal occupancies as determined by neutron diffraction supports the AXD occupancy methodology developed herein. Furthermore, the advantages innate to AXD (e.g., smaller crystal sizes, shorter collection times, and greater availability of synchrotron resources) versus neutron diffraction further support the need for its development as a standard technique.

Utilizing anomalous signals for element identification in macromolecular crystallography

AlphaFold2 has revolutionized structural biology by offering unparalleled accuracy in predicting protein structures. Traditional methods for determining protein structures, such as X-ray crystallography and cryo-electron microscopy, are often time-consuming and resource-intensive. AlphaFold2 provides models that are valuable for molecular replacement, aiding in model building and docking into electron density or potential maps. However, despite its capabilities, models from AlphaFold2 do not consistently match the accuracy of experimentally determined structures, need to be validated experimentally and currently miss some crucial information, such as post-translational modifications, ligands and bound ions. In this paper, the advantages are explored of collecting X-ray anomalous data to identify chemical elements, such as metal ions, which are key to understanding certain structures and functions of proteins. This is achieved through methods such as calculating anomalous difference Fourier maps or refining the imaginary component of the anomalous scattering factor f''. Anomalous data can serve as a valuable complement to the information provided by AlphaFold2 models and this is particularly significant in elucidating the roles of metal ions.

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.

Improvement of Single-Crystal Structures of Very Heavy Element Compounds by Refining Anomalous Dispersion Parameters

Refining the anomalous dispersion parameters of the four uranium compounds NaUF5, NaU2F9, Cs2(UO2)TiO4, and Cs2(UO2)Ti2O6 gave insights into the crystallographic model improvement of very heavy atoms. We found that the values for the dispersive and absorptive parts, f' and f″, closely followed the X-ray absorption spectra on their L3, L2, and L1 edges. The obtained values are sensitive to the chemical environment at each crystallographically independent position. An incorrect treatment of the anomalous dispersion correction can lead to a wrong crystallographic model. The above-mentioned, already published structures were improved by this process. General guidelines were given for the crystal structure determination of very heavy compounds. When using Mo Kα radiation with uranium compounds, the proximity of its energy to the uranium L-edges causes a noticeable effect.

Cu site differentiation in tetracopper(i) sulfide clusters enables biomimetic N2O reduction

Copper clusters feature prominently in both metalloenzymes and synthetic nanoclusters that mediate catalytic redox transformations of gaseous small molecules. Such reactions are critical to biological energy conversion and are expected to be crucial parts of renewable energy economies. However, the precise roles of individual metal atoms within clusters are difficult to elucidate, particularly for cluster systems that are dynamic under operating conditions. Here, we present a metal site-specific analysis of synthetic Cu4(μ4-S) clusters that mimic the Cu Z active site of the nitrous oxide reductase enzyme. Leveraging the ability to obtain structural snapshots of both inactive and active forms of the synthetic model system, we analyzed both states using resonant X-ray diffraction anomalous fine structure (DAFS), a technique that enables X-ray absorption profiles of individual metal sites within a cluster to be extracted independently. Using DAFS, we found that a change in cluster geometry between the inactive and active states is correlated to Cu site differentiation that is presumably required for efficient activation of N2O gas. More precisely, we hypothesize that the Cu δ+⋯Cu δ- pairs produced upon site differentiation are poised for N2O activation, as supported by computational modeling. These results provide an unprecedented level of detail on the roles of individual metal sites within the synthetic cluster system and how those roles interplay with cluster geometry to impact the reactivity function. We expect this fundamental knowledge to inform understanding of metal clusters in settings ranging from (bio)molecular to nanocluster to extended solid systems involved in energy conversion.

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.

Analysis of early intermediate states of the nitrogenase reaction by regularization of EPR spectra

Due to the complexity of the catalytic FeMo cofactor site in nitrogenases that mediates the reduction of molecular nitrogen to ammonium, mechanistic details of this reaction remain under debate. In this study, selenium- and sulfur-incorporated FeMo cofactors of the catalytic MoFe protein component from Azotobacter vinelandii are prepared under turnover conditions and investigated by using different EPR methods. Complex signal patterns are observed in the continuous wave EPR spectra of selenium-incorporated samples, which are analyzed by Tikhonov regularization, a method that has not yet been applied to high spin systems of transition metal cofactors, and by an already established grid-of-error approach. Both methods yield similar probability distributions that reveal the presence of at least four other species with different electronic structures in addition to the ground state E0. Two of these species were preliminary assigned to hydrogenated E2 states. In addition, advanced pulsed-EPR experiments are utilized to verify the incorporation of sulfur and selenium into the FeMo cofactor, and to assign hyperfine couplings of 33S and 77Se that directly couple to the FeMo cluster. With this analysis, we report selenium incorporation under turnover conditions as a straightforward approach to stabilize and analyze early intermediate states of the FeMo cofactor.

What triggers the coupling of proton transfer and electron transfer at the active site of nitrogenase?

In converting N2 to NH3 the enzyme nitrogenase utilises 8 electrons and 8 protons in the complete catalytic cycle. The source of the electrons is an Fe4S4 reductase protein (Fe-protein) which temporarily docks with the MoFe-protein that contains the catalytic active cofactor, FeMo-co, and an electron transfer cluster called the P cluster. The overall mechanism involves 8 repetitions of a cycle in which reduced Fe-protein docks with the MoFe-protein, one electron transfers to the P-cluster, and then to FeMo-co, followed by dissociation of the two proteins and re-reduction of the Fe-protein. Protons are supplied serially to FeMo-co by a Grotthuss proton translocation mechanism from the protein surface along a conserved chain of water molecules (a proton wire) that terminates near S atoms of the FeMo-co cluster [CFe7S9Mo(homocitrate)] where the multiple steps of the chemical conversions are effected. It is assumed that the chemical mechanisms use proton-coupled electron-transfer (PCET) and that H atoms (e- + H+) are involved in each of the hydrogenation steps. However there is neither evidence for, or mechanism proposed, for this coupling. Here I report calculations of cluster charge distribution upon electron addition, revealing that the added negative charge is on the S atoms of FeMo-co, which thereby become more basic, and able to trigger proton transfer from H3O+ waiting at the near end of the proton wire. This mechanism is supported by calculations of the dynamics of the proton transfer step, in which the barrier is reduced by ca. 3.5 kcal mol-1 and the product stabilised by ca. 7 kcal mol-1 upon electron addition. H tunneling is probable in this step. In nitrogenase it is electron transfer that triggers proton transfer.

A database overview of metal-coordination distances in metalloproteins

Metalloproteins are ubiquitous in all living organisms and take part in a very wide range of biological processes. For this reason, their experimental characterization is crucial to obtain improved knowledge of their structure and biological functions. The three-dimensional structure represents highly relevant information since it provides insight into the interaction between the metal ion(s) and the protein fold. Such interactions determine the chemical reactivity of the bound metal. The available PDB structures can contain errors due to experimental factors such as poor resolution and radiation damage. A lack of use of distance restraints during the refinement and validation process also impacts the structure quality. Here, the aim was to obtain a thorough overview of the distribution of the distances between metal ions and their donor atoms through the statistical analysis of a data set based on more than 115 000 metal-binding sites in proteins. This analysis not only produced reference data that can be used by experimentalists to support the structure-determination process, for example as refinement restraints, but also resulted in an improved insight into how protein coordination occurs for different metals and the nature of their binding interactions. In particular, the features of carboxylate coordination were inspected, which is the only type of interaction that is commonly present for nearly all metals.

Analyses of the electronic structures of FeFe-cofactors compared with those of FeMo- and FeV-cofactors and their P-clusters

The electronic structures of FeFe-cofactors (FeFe-cos) in resting and turnover states, together with their PN clusters from iron-only nitrogenases, have been calculated using the bond valence method, and their crystallographic data were reported recently and deposited in the Protein Data Bank (PDB codes: 8BOQ and 8OIE). The calculated results have also been compared with those of their homologous Mo- and V-nitrogenases. For FeFe-cos in the resting state, Fe1/2/4/5/6/7/8 atoms are prone to Fe3+, while the Fe3 atom shows different degrees of mixed valences. The results support that the Fe8 atom at the terminal positions of FeFe-cos possesses the same oxidation states as the Mo3+/V3+ atoms of FeMo-/FeV-cos. In the turnover state, the overall oxidation state of FeFe-co is slightly reduced than those in the resting species, and its electronic configuration is rearranged after the substitution of S2B with OH, compatible with those found in CO-bound FeV-co. Moreover, the calculations give the formal oxidation states of 6Fe2+-2Fe3+ for the electronic structures of PN clusters in Fe-nitrogenases. By the comparison of Mo-, V- and Fe-nitrogenases, the overall oxidation levels of 7Fe atoms (Fe1-Fe7) for both FeFe- and FeMo-cos in resting states are found to be higher than that of FeV-co. For the PN clusters in MoFe-, VFe- and FeFe-proteins, they all exhibit a strong reductive character.

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.

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

Uncovering a CF3 Effect on X-ray Absorption Energies of [Cu(CF3 )4 ]- and Related Copper Compounds by Using Resonant Diffraction Anomalous Fine Structure (DAFS) Measurements

Understanding the electronic structures of high-valent metal complexes aids the advancement of metal-catalyzed cross coupling methodologies. A prototypical complex with formally high valency is [Cu(CF3 )4 ]- (1), which has a formal Cu(III) oxidation state but whose physical analysis has led some to a Cu(I) assignment in an inverted ligand field model. Recent examinations of 1 by X-ray spectroscopies have led previous authors to contradictory conclusions, motivating the re-examination of its X-ray absorption profile here by a complementary method, resonant diffraction anomalous fine structure (DAFS). From analysis of DAFS measurements for a series of seven mononuclear Cu complexes including 1, here it is shown that there is a systematic trifluoromethyl effect on X-ray absorption that blue shifts the resonant Cu K-edge energy by 2-3 eV per CF3 , completely accounting for observed changes in DAFS profiles between formally Cu(III) complexes like 1 and formally Cu(I) complexes like (Ph3 P)3 CuCF3 (3). Thus, in agreement with the inverted ligand field model, the data presented herein imply that 1 is best described as containing a Cu(I) ion with dn count approaching 10.

On the Shoulders of Giants-Reaching for Nitrogenase

Only a single enzyme system-nitrogenase-carries out the conversion of atmospheric N2 into bioavailable ammonium, an essential prerequisite for all organismic life. The reduction of this inert substrate at ambient conditions poses unique catalytic challenges that strain our mechanistic understanding even after decades of intense research. Structural biology has added its part to this greater tapestry, and in this review, I provide a personal (and highly biased) summary of the parts of the story to which I had the privilege to contribute. It focuses on the crystallographic analysis of the three isoforms of nitrogenases at high resolution and the binding of ligands and inhibitors to the active-site cofactors of the enzyme. In conjunction with the wealth of available biochemical, biophysical, and spectroscopic data on the protein, this has led us to a mechanistic hypothesis based on an elementary mechanism of repetitive hydride formation and insertion.

Catalysis and structure of nitrogenases

In providing bioavailable nitrogen as building blocks for all classes of biomacromolecules, biological nitrogen fixation is an essential process for all organismic life. Only a single enzyme, nitrogenase, performs this task at ambient conditions and with ATP as an energy source. The assembly of the complex iron-sulfur enzyme nitrogenase and its catalytic mechanism remains a matter of intense study. Recent progress in the structural analysis of the three known isoforms of nitrogenase-differentiated primarily by the heterometal in their active site cofactor-has revealed a degree of structural plasticity of these clusters that suggest two distinct binding sites for substrates and reaction intermediates. A mechanistic proposal based on this finding integrates most of the available experimental data. Furthermore, the first applications of high-resolution cryo-electron microscopy have highlighted further dynamic conformational changes. Structures obtained under turnover conditions support the proposed alternating half-site reactivity in the C2-symmetric nitrogenase complex.

Metal Site-Specific Electrostatic Field Effects on a Tricopper(I) Cluster Probed by Resonant Diffraction Anomalous Fine Structure (DAFS)

Studies of multinuclear metal complexes are greatly enhanced by resonant diffraction measurements, which probe X-ray absorption profiles of crystallographically independent metal sites within a cluster. In particular, X-ray diffraction anomalous fine structure (DAFS) analysis provides data that can be interpreted akin to site-specific XANES, allowing for differences in metal K-edge resonances to be deconvoluted even for different metal sites within a homometallic system. Despite the prevalence of Cu-containing clusters in biology and energy science, DAFS has yet to be used to analyze multicopper complexes of any type until now. Here, we report an evaluation of trends using a series of strategically chosen Cu(I) and Cu(II) complexes to determine how energy dependencies of anomalous scattering factors are impacted by coordination geometry, ligand shell, cluster nuclearity, and oxidation state. This calibration data is used to analyze a formally tricopper(I) complex that was found by DAFS to be site-differentiated due to the unsymmetrical influence on different Cu sites of the electrostatic field from a proximal K+ cation.

Facing the phase problem

The marvel of X-ray crystallography is the beauty and precision of the atomic structures deduced from diffraction patterns. Since these patterns record only amplitudes, phases for the diffracted waves must also be evaluated for systematic structure determination. Thus, we have the phase problem as a central complication, both intellectually for the field and practically so for many analyses. Here, I discuss how we - myself, my laboratory and the diffraction community - have faced the phase problem, considering the evolution of methods for phase evaluation as structural biology developed to the present day. During the explosive growth of macromolecular crystallography, practice in diffraction analysis evolved from a universal reliance on isomorphous replacement to the eventual domination of anomalous diffraction for de novo structure determination. As the Protein Data Bank (PDB) grew and familial relationships among proteins became clear, molecular replacement overtook all other phasing methods; however, experimental phasing remained essential for molecules without obvious precedents, with multi- and single-wavelength anomalous diffraction (MAD and SAD) predominating. While the mathematics-based direct methods had proved to be inadequate for typical macromolecules, they returned to crack substantial selenium substructures in SAD analyses of selenomethionyl proteins. Native SAD, exploiting the intrinsic S and P atoms of biomolecules, has become routine. Selenomethionyl SAD and MAD were the mainstays of structural genomics efforts to populate the PDB with novel proteins. A recent dividend has been paid in the success of PDB-trained artificial intelligence approaches for protein structure prediction. Currently, molecular replacement with AlphaFold models often obviates the need for experimental phase evaluation. For multiple reasons, we are now unfazed by the phase problem. Cryo-EM analysis is an attractive alternative to crystallography for many applications faced by today's structural biologists. It simply finesses the phase problem; however, the principles and procedures of diffraction analysis remain pertinent and are adopted in single-particle cryo-EM studies of biomolecules.

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.

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.

The binding of reducible N2 in the reaction domain of nitrogenase

The binding of N₂ to FeMo-co, the catalytic site of the enzyme nitrogenase, is central to the conversion to NH₃, but also has a separate role in promoting the N₂-dependent HD reaction (D₂ + 2H⁺ + 2e^- → 2HD). The protein surrounding FeMo-co contains a clear channel for ingress of N₂, directly towards the exo-coordination position of Fe2, a position which is outside the catalytic reaction domain. This led to the hypothesis [I. Dance, Dalton Trans., 2022, 51, 12717] of 'promotional' N₂ bound at exo-Fe2, and a second 'reducible' N₂ bound in the reaction domain, specifically the endo-coordination position of Fe2 or Fe6. The range of possibilities for the binding of reducible N₂ in the presence of bound promotional N₂ is described here, using density functional simulations with a 486 atom model of the active site and surrounding protein. The pathway for ingress of the second N₂ through protein, past the first N₂ at exo-Fe2, and tumbling into the binding domain between Fe2 and Fe6, is described. The calculations explore 24 structures involving 6 different forms of hydrogenated FeMo-co, including structures with S2BH unhooked from Fe2 but tethered to Fe6. The calculations use the most probable electronic states. End-on (η¹) binding of N₂ at the endo position of either Fe2 or Fe6 is almost invariably exothermic, with binding potential energies ranging up to -18 kcal mol^-1. Many structures have binding energies in the range -6 to -14 kcal mol^-1. The relevant entropic penalty for N₂ binding from a diffusible position within the protein is estimated to be 4 kcal mol^-1, and so the binding free energies for reducible N₂ are suitably negative. N₂ binding at endo-Fe2 is stronger than at endo-Fe6 in three of the six structure categories. In many cases the reaction domain containing reducible N₂ is expanded. These results inform computational simulation of the subsequent steps in which surrounding H atoms transfer to reducible N₂.

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.

The HD Reaction of Nitrogenase: a Detailed Mechanism

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

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.

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.

Dinitrogen Binding and Activation: Bonding Analyses of Stable V(III/I)-N2-V(III/I) Complexes by the EDA-NOCV Method from the Perspective of Vanadium Nitrogenase

The FeVco cofactor of nitrogenase (VFe₇S₈(CO₃)C) is an alternative in the molybdenum (Mo)-deficient free soil living azotobacter vinelandii. The rate of N₂ reduction to NH₃ by FeVco is a few times higher than that by FeMoco (MoFe₇S₉C) at low temperature. It provides a N source in the form of ammonium ions to the soil. This biochemical NH₃ synthesis is an alternative to the industrial energy-demanding production of NH₃ by the Haber-Bosch process. The role of vanadium has not been clearly understood yet, which has led chemists to come up with several stable V-N₂ complexes which have been isolated and characterized in the laboratory over the past three decades. Herein, we report the EDA-NOCV analyses of dinitrogen-bonded stable complexes V(III/I)-N₂ (1-4) to provide deeper insights into the fundamental bonding aspects of V-N₂ bond, showing the interacting orbitals and corresponding pairwise orbital interaction energies (ΔE _orb(n)). The computed intrinsic interaction energy (ΔE _int) of V-N₂-V bonds is significantly higher than those of the previously reported Fe-N₂-Fe bonds. Covalent interaction energy (ΔE _orb) is more than double the electrostatic interaction energy (ΔE _elstat) of V-N₂-V bonds. ΔE _int values of V-N₂-V bonds are in the range of -172 to -204 kcal/mol. The V → N₂ ← V π-backdonation is four times stronger than V ← N₂ → V σ-donation. V-N₂ bonds are much more covalent in nature than Fe-N₂ bonds.

Refinement of anomalous dispersion correction parameters in single-crystal structure determinations

Correcting for anomalous dispersion is part of any refinement of an X-ray dif-fraction crystal structure determination. The procedure takes the inelastic scattering in the diffraction experiment into account. This X-ray absorption effect is specific to each chemical compound and is particularly sensitive to radiation energies in the region of the absorption edges of the elements in the compound. Therefore, the widely used tabulated values for these corrections can only be approximations as they are based on calculations for isolated atoms. Features of the unique spatial and electronic environment that are directly related to the anomalous dispersion are ignored, although these can be observed spectroscopically. This significantly affects the fit between the crystallographic model and the measured intensities when the excitation wavelength in an X-ray diffraction experiment is close to an element's absorption edge. Herein, we report on synchrotron multi-wavelength single-crystal X-ray diffraction, as well as X-ray absorption spectroscopy experiments which we performed on the mol-ecular compound Mo(CO)6 at energies around the molybdenum K edge. The dispersive (f') and absorptive (f'') terms of the anomalous dispersion can be refined as independent parameters in the full-matrix least-squares refinement. This procedure has been implemented as a new feature in the well-established OLEX2 software suite. These refined parameters are in good agreement with the independently recorded X-ray absorption spectrum. The resulting crystallographic models show significant improvement compared to those employing tabulated values.

Calculating the chemical mechanism of nitrogenase: new working hypotheses

The enzyme nitrogenase converts N₂ to NH₃ with stoichiometry N₂ + 8H⁺ + 8e^- → 2NH₃ + H₂. The mechanism is chemically complex with multiple steps that must be consistent with much accumulated experimental information, including exchange of H₂ and N₂ and the N₂-dependent hydrogenation of D₂ to HD. Previous investigations have developed a collection of working hypotheses that guide ongoing density functional investigations of mechanistic steps and sequences. These include (i) hypotheses about the serial provision of protons and their conversion to H atoms bonded to S and Fe atoms of the FeMo-co catalytic site, (ii) the migration of H atoms over the surface of FeMo-co, (iii) the roles of His195, (iv) identification of three protein channels, one for the ingress of N₂, a separate pathway for the passage of exogenous H₂ (D₂) and product H₂ (HD), and a hydrophilic pathway for egress of product NH₃. Two additional working hypotheses are described in this paper. N₂ passing along the N₂ channel approaches and binds end-on to the exo coordination position of Fe2, with favourable energetics when FeMo-co is pre-hydrogenated. This exo-Fe2-N₂ is apparently not reduced but has a promotional role by expanding the reaction zone. A second N₂ can enter via the N₂ ingress channel and bind at the endo-Fe6 position, where it is surrounded by H atom donors suitable for the N₂ → NH₃ conversion. It is proposed that this endo-Fe6 position is also the binding site for H₂ (generated or exogenous), accounting for the competitive inhibition of N₂ reduction by H₂. The HD reaction occurs at the endo-Fe6 site, promoted by N₂ at the exo-Fe2 site. The second hypothesis concerns the most stable electronic states of FeMo-co with ligands bound at Fe2 and Fe6, and provides a protocol for management of electronic states in mechanism calculations.

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

Gases like H2, N2, CO2, 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 N2, CO2, and CO and the production of H2 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 N2 fixation by nitrogenase and H2 production by hydrogenase as well as CO2 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.

Frontiers in metalloprotein crystallography and cryogenic electron microscopy

Metalloproteins comprise at least a third of all proteins that utilize redox properties of transition metals on their own or as parts of cofactors. The development of third generation storage ring sources and X-ray free-electron lasers with femtosecond pulses in the first decade of the 21st century has transformed metalloprotein crystallography. In the past decade, cryogenic-electron microscopy single-particle analysis, which does not require crystallization of biological samples has been extensively utilized, particularly for membrane-bound metalloprotein systems. Here, we explore recent frontiers in metalloprotein crystallography and cryogenic electron microscopy, organized for convenience under three metalloprotein-centered biological cycles, focusing on contributions from each technique, their synergy and the ability to preserve metals' redox states when subjected to a particular probe.

Analysis of the Geometric and Electronic Structure of Spin-Coupled Iron-Sulfur Dimers with Broken-Symmetry DFT: Implications for FeMoco

The open-shell electronic structure of iron-sulfur clusters presents considerable challenges to quantum chemistry, with the complex iron-molybdenum cofactor (FeMoco) of nitrogenase representing perhaps the ultimate challenge for either wavefunction or density functional theory. While broken-symmetry density functional theory has seen some success in describing the electronic structure of such cofactors, there is a large exchange-correlation functional dependence in calculations that is not fully understood. In this work, we present a geometric benchmarking test set, FeMoD11, of synthetic spin-coupled Fe-Fe and Mo-Fe dimers, with relevance to the molecular and electronic structure of the Mo-nitrogenase FeMo cofactor. The reference data consists of high-resolution crystal structures of metal dimer compounds in different oxidation states. Multiple density functionals are tested on their ability to reproduce the local geometry, specifically the Fe-Fe/Mo-Fe distance, for both antiferromagnetically coupled and ferromagnetically coupled dimers via the broken-symmetry approach. The metal-metal distance is revealed not only to be highly sensitive to the amount of exact exchange in the functional but also to the specific exchange and correlation functionals. For the antiferromagnetically coupled dimers, the calculated metal-metal distance correlates well with the covalency of the bridging metal-ligand bonds, as revealed via the corresponding orbital analysis, Hirshfeld S/Fe charges, and Fe-S Mayer bond order. Superexchange via bridging ligands is expected to be the dominant interaction in these dimers, and our results suggest that functionals that predict accurate Fe-Fe and Mo-Fe distances describe the overall metal-ligand covalency more accurately and in turn the superexchange of these systems. The best performing density functionals of the 16 tested for the FeMoD11 test set are revealed to be either the nonhybrid functionals r2SCAN and B97-D3 or hybrid functionals with 10-15% exact exchange: TPSSh and B3LYP*. These same four functionals are furthermore found to reproduce the high-resolution X-ray structure of FeMoco well according to quantum mechanics/molecular mechanics (QM/MM) calculations. Almost all nonhybrid functionals systematically underestimate Fe-Fe and Mo-Fe distances (with r2SCAN and B97-D3 being the sole exceptions), while hybrid functionals with >15% exact exchange (including range-separated hybrid functionals) overestimate them. The results overall suggest r2SCAN, B97-D3, TPSSh, and B3LYP* as accurate density functionals for describing the electronic structure of iron-sulfur clusters in general, including the complex FeMoco cluster of nitrogenase.

Determining the oxidation state of elements by X-ray crystallography

Protein-mediated redox reactions play a critical role in many biological processes and often occur at centres that contain metal ions as cofactors. In order to understand the exact mechanisms behind these reactions it is important to not only characterize the three-dimensional structures of these proteins and their cofactors, but also to identify the oxidation states of the cofactors involved and to correlate this knowledge with structural information. The only suitable approach for this based on crystallographic measurements is spatially resolved anomalous dispersion (SpReAD) refinement, a method that has been used previously to determine the redox states of metals in iron–sulfur cluster-containing proteins. In this article, the feasibility of this approach for small, non-iron–sulfur redox centres is demonstrated by employing SpReAD analysis to characterize Sulfolobus tokodaii sulerythrin, a ruberythrin-like protein that contains a binuclear metal centre. Differences in oxidation states between the individual iron ions of the binuclear metal centre are revealed in sulerythrin crystals treated with H2O2. Furthermore, data collection at high X-ray doses leads to photoreduction of this metal centre, showing that careful control of the total absorbed dose is a prerequisite for successfully determining the oxidation state through SpReAD analysis.

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.

Complete genome of the thermophilic purple sulfur Bacterium Thermochromatium tepidum compared to Allochromatium vinosum and other Chromatiaceae

The complete genome sequence of the thermophilic purple sulfur bacterium Thermochromatium tepidum strain MC^T (DSM 3771^T) is described and contrasted with that of its mesophilic relative Allochromatium vinosum strain D (DSM 180^T) and other Chromatiaceae. The Tch. tepidum genome is a single circular chromosome of 2,958,290 base pairs with no plasmids and is substantially smaller than the genome of Alc. vinosum. The Tch. tepidum genome encodes two forms of RuBisCO and contains nifHDK and several other genes encoding a molybdenum nitrogenase but lacks a gene encoding a protein that assembles the Fe-S cluster required to form a functional nitrogenase molybdenum-iron cofactor, leaving the phototroph phenotypically Nif^-. Tch. tepidum contains genes necessary for oxidizing sulfide to sulfate as photosynthetic electron donor but is genetically unequipped to either oxidize thiosulfate as an electron donor or carry out assimilative sulfate reduction, both of which are physiological hallmarks of Alc. vinosum. Also unlike Alc. vinosum, Tch. tepidum is obligately phototrophic and unable to grow chemotrophically in darkness by respiration. Several genes present in the Alc. vinosum genome that are absent from the genome of Tch. tepidum likely contribute to the major physiological differences observed between these related purple sulfur bacteria that inhabit distinct ecological niches.

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.

Revealing redox isomerism in trichromium imides by anomalous diffraction

In polynuclear biological active sites, multiple electrons are needed for turnover, and the distribution of these electrons among the metal sites is affected by the structure of the active site. However, the study of the interplay between structure and redox distribution is difficult not only in biological systems but also in synthetic polynuclear clusters since most redox changes produce only one thermodynamically stable product. Here, the unusual chemistry of a sterically hindered trichromium complex allowed us to probe the relationship between structural and redox isomerism. Two structurally isomeric trichromium imides were isolated: asymmetric terminal imide (^tbsL)Cr₃(NDipp) and symmetric, μ³-bridging imide (^tbsL)Cr₃(μ³–NBn) ((^tbsL)⁶⁻ = (1,3,5-C₆H₉(NC₆H₄-o-NSi^tBuMe₂)₃)⁶⁻). Along with the homovalent isocyanide adduct (^tbsL)Cr₃(CNBn) and the bisimide (^tbsL)Cr₃(μ³–NPh)(NPh), both imide isomers were examined by multiple-wavelength anomalous diffraction (MAD) to determine the redox load distribution by the free refinement of atomic scattering factors. Despite their compositional similarities, the bridging imide shows uniform oxidation of all three Cr sites while the terminal imide shows oxidation at only two Cr sites. Further oxidation from the bridging imide to the bisimide is only borne at the Cr site bound to the second, terminal imido fragment. Thus, depending on the structural motifs present in each [Cr₃] complex, MAD revealed complete localization of oxidation, partial localization, and complete delocalization, all supported by the same hexadentate ligand scaffold.

Application of high-resolution Multiwavelength Anomalous Diffraction (MAD) allows the assignment of localized, partly delocalized, and fully delocalized oxidation in a series of trichromium imide isomers.

Partial synthetic models of FeMoco with sulfide and carbyne ligands: Effect of interstitial atom in nitrogenase active site

Nitrogen-fixing organisms perform dinitrogen reduction to ammonia at an Fe-M (M = Mo, Fe, or V) cofactor (FeMco) of nitrogenase. FeMco displays eight metal centers bridged by sulfides and a carbide having the MFe₇S₈C cluster composition. The role of the carbide ligand, a unique motif in protein active sites, remains poorly understood. Toward addressing how the carbon bridge affects the physical and chemical properties of the cluster, we isolated synthetic models of subsite MFe₃S₃C displaying sulfides and a chelating carbyne ligand. We developed synthetic protocols for structurally related clusters, [Tp*M'Fe₃S₃X]^n-, where M' = Mo or W, the bridging ligand X = CR, N, NR, S, and Tp* = Tris(3,5-dimethyl-1-pyrazolyl)hydroborate, to study the effects of the identity of the heterometal and the bridging X group on structure and electrochemistry. While the nature of M' results in minor changes, the chelating, μ₃-bridging carbyne has a large impact on reduction potentials, being up to 1 V more reducing compared to nonchelating N and S analogs.

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.

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.

Synthesis and characterisation of a very low-coordinate diferrous [2Fe-2S]0 unit

Here we present the synthesis of a unique diferrous [2Fe-2S]⁰ complex with only three-coordinate iron ions via reduction of a four-coordinate diferric [2Fe-2S]²⁺ complex with concomitant ligand loss. The obtained compounds were thoroughly examined for their properties (e.g. by ⁵⁷Fe Mössbauer spectroscopy and magnetic susceptibility measurements). Facile cleavage of the [2Fe-2S] rhombus, commonly seen as rather stable, by CS₂ is also shown.

Bayesian phase difference estimation: a general quantum algorithm for the direct calculation of energy gaps

Quantum computers can perform full configuration interaction (full-CI) calculations by utilising the quantum phase estimation (QPE) algorithms including Bayesian phase estimation (BPE) and iterative quantum phase estimation (IQPE). In these quantum algorithms, the time evolution of wave functions for atoms and molecules is simulated conditionally with an ancillary qubit as the control, which make implementation to real quantum devices difficult. Also, most of the problems in chemistry discuss energy differences between two electronic states rather than total energies themselves, and thus direct calculations of energy gaps are promising for future applications of quantum computers to real chemistry problems. In the race of finding efficient quantum algorithms to solve quantum chemistry problems, we test a Bayesian phase difference estimation (BPDE) algorithm, which is a general algorithm to calculate the difference of two eigenphases of unitary operators in the several cases of the direct calculations of energy gaps between two electronic states on quantum computers, including vertical ionisation energies, singlet-triplet energy gaps, and vertical excitation energies. In the BPDE algorithm, state preparation is carried out conditionally on the ancillary qubit, and the time evolution of the wave functions in superposition of two electronic states are executed unconditionally. Based on our test, we conclude that BPDE is capable of computing the energy gap with an accuracy similar to BPE without controlled-time evolution simulations and with the smaller number of iterations in Bayesian optimisations.

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.

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.

Two ligand-binding sites in CO-reducing V nitrogenase reveal a general mechanistic principle

Besides its role in biological nitrogen fixation, vanadium-containing nitrogenase also reduces carbon monoxide (CO) to hydrocarbons, in analogy to the industrial Fischer-Tropsch process. The protein yields 93% of ethylene (C₂H₄), implying a C-C coupling step that mandates the simultaneous binding of two CO at the active site FeV cofactor. Spectroscopic data indicated multiple CO binding events, but structural analyses of Mo and V nitrogenase only confirmed a single site. Here, we report the structure of a two CO-bound state of V nitrogenase at 1.05 Å resolution, with one μ-bridging and one terminal CO molecule. This additional, specific ligand binding site suggests a mechanistic route for CO reduction and hydrocarbon formation, as well as a second access pathway for protons required during the reaction. Moreover, carbonyls are strong-field ligands that are chemically similar to mechanistically relevant hydrides that may be formed and used in a fully analogous fashion.

A Molecular Low-Coordinate [Fe-S-Fe] Unit in Three Oxidation States

A [Fe-S-Fe] subunit with a single sulfide bridging two low-coordinate iron ions is the supposed active site of the iron-molybdenum co-factor (FeMoco) of nitrogenase. Here we report a dinuclear monosulfido bridged diiron(II) complex with a similar complex geometry that can be oxidized stepwise to diiron(II/III) and diiron(III/III) complexes while retaining the [Fe-S-Fe] core. The series of complexes has been characterized crystallographically, and electronic structures have been studied using, inter alia, 57 Fe Mössbauer spectroscopy and SQUID magnetometry. Further, cleavage of the [Fe-S-Fe] unit by CS2 is presented.

The electronic structure of FeV-cofactor in vanadium-dependent nitrogenase

The electronic structure of the active-site metal cofactor (FeV-cofactor) of resting-state V-dependent nitrogenase has been an open question, with earlier studies indicating that it exhibits a broad S = 3/2 EPR signal (Kramers state) having g values of ∼4.3 and 3.8, along with suggestions that it contains metal-ions with valencies [1V3+, 3Fe3+, 4Fe2+]. In the present work, genetic, biochemical, and spectroscopic approaches were combined to reveal that the EPR signals previously assigned to FeV-cofactor do not correlate with active VFe-protein, and thus cannot arise from the resting-state of catalytically relevant FeV-cofactor. It, instead, appears resting-state FeV-cofactor is either diamagnetic, S = 0, or non-Kramers, integer-spin (S = 1, 2 etc.). When VFe-protein is freeze-trapped during high-flux turnover with its natural electron-donating partner Fe protein, conditions which populate reduced states of the FeV-cofactor, a new rhombic S = 1/2 EPR signal from such a reduced state is observed, with g = [2.18, 2.12, 2.09] and showing well-defined 51V (I = 7/2) hyperfine splitting, a iso = 110 MHz. These findings indicate a different assignment for the electronic structure of the resting state of FeV-cofactor: S = 0 (or integer-spin non-Kramers state) with metal-ion valencies, [1V3+, 4Fe3+, 3Fe2+]. Our findings suggest that the V3+ does not change valency throughout the catalytic cycle.

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.