Revisiting the Mössbauer Isomer Shifts of the FeMoco Cluster of Nitrogenase and the Cofactor Charge

Divergent Interaction of (Iso)nitriles with a Linear Iron(I) Silylamide─A Combined Structural, Spectroscopic, and Computational Study

Nitriles and isonitriles are important σ-donor ligands in coordination chemistry. Isonitriles also function in low-valent complexes as π-acceptor ligands similar to CO. Herein we present the unusual behavior of the highly reducing, high-spin iron(I) complex [Fe(hmds)2]- toward these compound classes. Rare examples of side-on coordination of nitriles to the metal center are observed. Insights gained by 57Fe Mössbauer spectroscopy as well as DFT and CASSCF calculations give an interplay between the resonance structures of not only an iron(I) π-complex and an iron(III) metallacycle but also point to the importance of an iron(II) nitrile radical anion. For an aromatic isonitrile end-on coordination is observed, which is best described as an iron(I) complex with only minor unpaired spin transfer onto the isonitrile. For aliphatic isonitriles, the selective R-CN bond cleavage occurs and yields stoichiometric mixtures of alkyl iron(II) and cyanido iron(II) complexes. Attempts to isolate presumed (iso)nitrile radical anions void of 3d-metal coordination give for the reaction of an aromatic isonitrile with KC8 facile reductive coupling to the corresponding diamido acetylene.

Low-Coordinate Iron Hydride Chemistry at an N,N,C-Heteroscorpionate Platform

Locally high-spin iron hydrides are proposed to play a critical role as intermediates in iron-molybdenum cofactor (FeMoco)-catalyzed N2 fixation. Inspired by these biological systems, we report herein our initial investigations into low-coordinate iron hydride chemistry supported by our N,N,C-heteroscorpionate ligands. Those ligands with smaller steric profiles are unable to completely suppress the formation of a binuclear [Fe(μ2-H)]2 complex; however, the incorporation of more substantial steric bulk allows for the isolation of a rare example of a terminal, high-spin (S = 2) Fe2+ hydride. Fourier transform infrared spectroscopy suggests an unusually weak Fe-H bond in the latter molecule. Mössbauer spectroscopies, coupled with density functional theory calculations, highlights the substantial influence of the terminal hydride ligand on 57Fe isomer shift.

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.

Restricted Open-Shell Hartree-Fock Method for a General Configuration State Function Featuring Arbitrarily Complex Spin-Couplings

In this work, we present a general spin restricted open-shell Hartree-Fock (ROHF) implementation that is able to generate self-consistent field (SCF) wave functions for an arbitrary configuration state function (CSF). These CSFs can contain an arbitrary number of unpaired electrons in arbitrary spin-couplings. The resulting method is named CSF-ROHF. We demonstrate that starting from the ROHF energy expression, for example, the one given by Edwards and Zerner, it is possible to obtain the values of the ROHF vector-coupling coefficients by setting up an open-shell for each group of consecutive parallel-coupled spins dictated by the unique spin-coupling pattern of any given CSF. To achieve this important and nontrivial goal, we employ the machinery of the iterative configuration expansion configuration interaction (ICE-CI) method, which is able to tackle general CI problems on the basis of spin-adapted CSFs. This development allows for the efficient generation of SCF spin-eigenfunctions for systems with complex spin-coupling patterns, such as polymetallic chains and metal clusters, while maintaining SCF scaling with system size (quadratic or less, depending on the specific algorithm and approximations chosen).

Unusual S=1 Four-Coordinate Fe(IV) Complexes Supported by Bisamide Ligands: Syntheses, Characterization, and Electronic Structures

The catalytic relevance of Fe(IV) species in non-heme iron catalysis has motivated synthetic advances in well-defined five- and six-coordinate Fe(IV) complexes for a better understanding of their fundamental electronic structures and reactivities. Herein, we report the syntheses of FeDipp2 and FeMes2, a pair of unusual four-coordinate non-heme formally Fe(IV) complexes with S=1 ground states supported by strongly donating bisamide ligands. By combining spectroscopic characterization and computational modeling, we found that small variations in ligand aryl substituents resulted in substantial changes in both structures and bonding. This work highlights the strong donor capabilities and modularity of the bisamide ligand set. More broadly, it is a critical contribution to the utilization of ligand design to modulate molecular geometries and electronic structures of low-coordinate, high-valent iron complexes.

Machine Learning-Enhanced Quantum Chemistry-Assisted Refinement of the Active Site Structure of Metalloproteins

Understanding the fine structural details of inhibitor binding at the active site of metalloenzymes can have a profound impact on the rational drug design targeted to this broad class of biomolecules. Structural techniques such as NMR, cryo-EM, and X-ray crystallography can provide bond lengths and angles, but the uncertainties in these measurements can be as large as the range of values that have been observed for these quantities in all the published structures. This uncertainty is far too large to allow for reliable calculations at the quantum chemical (QC) levels for developing precise structure-activity relationships or for improving the energetic considerations in protein-inhibitor studies. Therefore, the need arises to rely upon computational methods to refine the active site structures well beyond the resolution obtained with routine application of structural methods. In a recent paper, we have shown that it is possible to refine the active site of cobalt(II)-substituted MMP12, a metalloprotein that is a relevant drug target, by matching to the experimental pseudocontact shifts (PCS) those calculated using multireference ab initio QC methods. The computational cost of this methodology becomes a significant bottleneck when the starting structure is not sufficiently close to the final one, which is often the case with biomolecular structures. To tackle this problem, we have developed an approach based on a neural network (NN) and a support vector regression (SVR) and applied it to the refinement of the active site structure of oxalate-inhibited human carbonic anhydrase 2 (hCAII), another prototypical metalloprotein target. The refined structure gives a remarkably good agreement between the QC-calculated and the experimental PCS. This study not only contributes to the knowledge of CAII but also demonstrates the utility of combining machine learning (ML) algorithms with QC calculations, offering a promising avenue for investigating other drug targets and complex biological systems in general.

Trimethylsilyldiazomethane Disassembly at a Three-Fold Symmetric Iron Site

The reaction of equimolar trimethylsilyldiazomethyllithium (LiTMSD) with high spin (S = 2) PhB(AdIm)3FeCl (PhB(AdIm)3- = tris(3-adamantylimidazol-2-ylidene)phenylborate) affords the corresponding N-nitrilimido complex PhB(AdIm)3Fe-N═N═C(SiMe3). This complex can be converted to the thermodynamically more favorable C-isocyanoamido isomer PhB(AdIm)3Fe-C═N═N(SiMe3) by reaction with an additional equivalent of LiTMSD. While the iron(II) complexes are four-coordinate, the diazomethane is bound side-on in the iron(I) congener PhB(AdIm)3Fe(N,N'-κ2-N2C(H)Si(CH3)3). The latter complex adopts high spin (S = 3/2) ground state and features an unusually weak C-H bond. Photolysis of the iron(II) complexes induces N═N bond cleavage, with the iron(II) cyanide PhB(AdIm)3Fe-C≡N and iron(IV) nitride PhB(AdIm)3Fe≡N complexes being the major products of the reaction. The same products are obtained when the iron(I) complex is photolyzed or treated with a fluoride source. The trimethylsilyldiazomethane-derived ligand disassembly reactions are contrasted with those observed for related tris(carbene)amine complexes.

Bionic Design of Iron-Doped MoSe2 Catalyst for Efficient Nitrogen Fixation

The catalytic system of biological nitrogen fixation in nature primarily relies on the "FeMo cofactor" within nitrogenase enzymes. Inspired by this natural structure, we have designed a bionic inorganic counterpart, iron-doped MoSe2, for the efficient electroreduction of dinitrogen to ammonia. The introduced Fe dopant significantly enhances nitrogen fixation activity of MoSe2. Furthermore, we constructed a heterostructure catalyst, the Fe-MoSe2/Mo2C with more versatile Mo valence states. The heterostructured electrocatalyst achieves an ammonia production rate of 3.38 μg cm-2 h-1, and a Faradaic efficiency of 30.8 %, which is ~5 fold higher than that of pristine MoSe2. This study presents a novel approach for designing bionic nitrogen fixation electrocatalysts.

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.

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.

From the Iron Pentacarbonyl Cation to Heteroleptic η6-arene Carbonyls and bis-η6-arene Cations

Partial ligand substitution at the iron pentacarbonyl radical cation generates novel half-sandwich complexes of the type [Fe(η6-arene)(CO)2]⋅+ (arene=1,3,5-tri-tert-butylbenzene, 1,3,5-trimethylbenzene, benzene and fluorobenzene). Of those, the bulkier 1,3,5-tri-tert-butylbenzene (mes*) derivative [Fe(mes*)(CO)2]⋅+ was fully characterized by XRD analysis, IR, NMR, cw-EPR, Mössbauer spectroscopy and cyclic voltammetry as the [Al(ORF)4]- (RF=C(CF3)3) salt. Chemical electronation, i. e., the single electron reduction, with decamethylferrocene generates neutral [Fe(mes*)(CO)2], whereas further deelectronation under CO-pressure leads to a dicationic three-legged [Fe(mes*)(CO)3]2+ salt with [Al(ORF)4]- counterion. The full substitution of the carbonyl ligands in [Fe(CO)5]⋅+[Al(ORF)4]- mainly resulted in disproportionation reactions, giving solid Fe(0) and the dicationic bis-arene salts [Fe(η6-arene)2]2+([Al(ORF)4]-)2 (arene=1,3,5-trimethylbenzene, benzene and fluorobenzene). Only by employing the very large fluoride bridged anion [F-{Al(ORF)3}2]-, it was possible to isolate an open shell bis-arene cation salt [Fe(C6H6)2]⋅+[F-{Al(ORF)3}2]-. The highly reactive cation was characterized by XRD analysis, cw-EPR, Mössbauer spectroscopy and cyclic voltammetry. The disproportionation of [Fe(C6H6)2]⋅+ salts to give solid Fe(0) and [Fe(C6H6)2]2+ salts was analyzed by a suitable cycle, revealing that the thermodynamic driving force for the disproportionation is a function of the size of the anion used and the polarity of the solvent.

Correlating Structure with Spectroscopy in Ascorbate Peroxidase Compound II

Structural and spectroscopic investigations of compound II in ascorbate peroxidase (APX) have yielded conflicting conclusions regarding the protonation state of the crucial Fe(IV) intermediate. Neutron diffraction and crystallographic data support an iron(IV)-hydroxo formulation, whereas Mössbauer, X-ray absorption (XAS), and nuclear resonance vibrational spectroscopy (NRVS) studies appear consistent with an iron(IV)-oxo species. Here we examine APX with spectroscopy-oriented QM/MM calculations and extensive exploration of the conformational space for both possible formulations of compound II. We establish that irrespective of variations in the orientation of a vicinal arginine residue and potential reorganization of proximal water molecules and hydrogen bonding, the Fe-O distances for the oxo and hydroxo forms consistently fall within distinct, narrow, and nonoverlapping ranges. The accuracy of geometric parameters is validated by coupled-cluster calculations with the domain-based local pair natural orbital approach, DLPNO-CCSD(T). QM/MM calculations of spectroscopic properties are conducted for all structural variants, encompassing Mössbauer, optical, X-ray absorption, and X-ray emission spectroscopies and NRVS. All spectroscopic observations can be assigned uniquely to an Fe(IV)═O form. A terminal hydroxy group cannot be reconciled with the spectroscopic data. Under no conditions can the Fe(IV)═O distance be sufficiently elongated to approach the crystallographically reported Fe-O distance. The latter is consistent only with a hydroxo species, either Fe(IV) or Fe(III). Our findings strongly support the Fe(IV)═O formulation of APX-II and highlight unresolved discrepancies in the nature of samples used across different experimental studies.

Protein3D: Enabling analysis and extraction of metal-containing sites from the Protein Data Bank with molSimplify

Metalloenzymes catalyze a wide range of chemical transformations, with the active site residues playing a key role in modulating chemical reactivity and selectivity. Unlike smaller synthetic catalysts, a metalloenzyme active site is embedded in a larger protein, which makes interrogation of electronic properties and geometric features with quantum mechanical calculations challenging. Here we implement the ability to fetch crystallographic structures from the Protein Data Bank and analyze the metal binding sites in the program molSimplify. We show the usefulness of the newly created protein3D class to extract the local environment around non-heme iron enzymes containing a two histidine motif and prepare 372 structures for quantum mechanical calculations. Our implementation of protein3D serves to expand the range of systems molSimplify can be used to analyze and will enable high-throughput study of metal-containing active sites in proteins.

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.

H2 formation from the E2-E4 states of nitrogenase

Nitrogenase is the only enzyme that can cleave the strong triple bond in N2, making nitrogen available for biological lifeforms. The active site is a MoFe7S9C cluster (the FeMo cluster) that binds eight electrons and protons during one catalytic cycle, giving rise to eight intermediate states E0-E7. It is experimentally known that N2 binds to the E4 state and that H2 is a compulsory byproduct of the reaction. However, formation of H2 is also an unproductive side reaction that should be avoided, especially in the early steps of the reaction mechanism (E2 and E3). Here, we study the formation of H2 for various structural interpretations of the E2-E4 states using combined quantum mechanical and molecular mechanical (QM/MM) calculations and four different density-functional theory methods. We find large differences in the predictions of the different methods. B3LYP strongly favours protonation of the central carbide ion and H2 cannot form from such structures. On the other hand, with TPSS, r2SCAN and TPSSh, H2 formation is strongly exothermic for all structures and En and therefore need strict kinetic control to be avoided. For the E2 state, the kinetic barriers for the low-energy structures are high enough to avoid H2 formation. However, for both the E3 and E4 states, all three methods predict that the best structure has two hydride ions bridging the same pair of Fe ions (Fe2 and Fe6) and these two ions can combine to form H2 with an activation barrier of only 29-57 kJ mol-1, corresponding to rates of 7 × 102 to 5 × 107 s-1, i.e. much faster than the turnover rate of the enzyme (1-5 s-1). We have also studied H-atom movements within the FeMo cluster, showing that the various protonation states can quite freely be interconverted (activation barriers of 12-69 kJ mol-1).

Protonation of Homocitrate and the E1 State of Fe-Nitrogenase Studied by QM/MM Calculations

Nitrogenase is the only enzyme that can cleave the strong triple bond in N2, making nitrogen available for biological life. There are three isozymes of nitrogenase, differing in the composition of the active site, viz., Mo, V, and Fe-nitrogenase. Recently, the first crystal structure of Fe-nitrogenase was presented. We have performed the first combined quantum mechanical and molecular mechanical (QM/MM) study of Fe-nitrogenase. We show with QM/MM and quantum-refinement calculations that the homocitrate ligand is most likely protonated on the alcohol oxygen in the resting E0 state. The most stable broken-symmetry (BS) states are the same as for Mo-nitrogenase, i.e., the three Noodleman BS7-type states (with a surplus of β spin on the eighth Fe ion), which maximize the number of nearby antiferromagnetically coupled Fe-Fe pairs. For the E1 state, we find that protonation of the S2B μ2 belt sulfide ion is most favorable, 14-117 kJ/mol more stable than structures with a Fe-bound hydride ion (the best has a hydride ion on the Fe2 ion) calculated with four different density-functional theory methods. This is similar to what was found for Mo-nitrogenase, but it does not explain the recent EPR observation that the E1 state of Fe-nitrogenase should contain a photolyzable hydride ion. For the E1 state, many BS states are close in energy, and the preferred BS state differs depending on the position of the extra proton and which density functional is used.

Mössbauer and Nuclear Resonance Vibrational Spectroscopy Studies of Iron Species Involved in N-N Bond Cleavage

Diketiminate-supported iron complexes are capable of cleaving the strong triple bond of N2 to give a tetra-iron complex with two nitrides (Rodriguez et al., Science, 2011, 334, 780-783). The mechanism of this reaction has been difficult to determine, but a transient green species was observed during the reaction that corresponds to a potential intermediate. Here, we describe studies aiming to identify the characteristics of this intermediate, using a range of spectroscopic techniques, including Mössbauer spectroscopy, electronic absorption spectroscopy, Raman spectroscopy, nuclear magnetic resonance (NMR) spectroscopy, and nuclear resonance vibrational spectroscopy (NRVS) complemented by density functional theory (DFT) calculations. We successfully elucidated the nature of the starting iron(II) species and the bis(nitride) species in THF solution, and in each case, THF breaks up the multiiron species. Various observations on the green intermediate species indicate that it has one N2 per two Fe atoms, has THF associated with it, and has NRVS features indicative of bridging N2. Computational models with a formally diiron(0)-N2 core are most consistent with the accumulated data, and on this basis, a mechanism for N2 splitting is suggested. This work shows the power of combining NRVS, Mössbauer, NMR, and vibrational spectroscopies with computations for revealing the nature of transient iron species during N2 cleavage.

Iron-Catalyzed Intermolecular C-H Amination Assisted by an Isolated Iron-Imido Radical Intermediate

Here we report the use of a base metal complex [(tBu pyrpyrr2 )Fe(OEt2 )] (1-OEt2 ) (tBu pyrpyrr2 2- =3,5-tBu2 -bis(pyrrolyl)pyridine) as a catalyst for intermolecular amination of Csp3 -H bonds of 9,10-dihydroanthracene (2 a) using 2,4,6-trimethyl phenyl azide (3 a) as the nitrene source. The reaction is complete within one hour at 80 °C using as low as 2 mol % 1-OEt2 with control in selectivity for single C-H amination versus double C-H amination. Catalytic C-H amination reactions can be extended to other substrates such as cyclohexadiene and xanthene derivatives and can tolerate a variety of aryl azides having methyl groups in both ortho positions. Under stoichiometric conditions the imido radical species [(tBu pyrpyrr2 )Fe{=N(2,6-Me2 -4-tBu-C6 H2 )] (1-imido) can be isolated in 56 % yield, and spectroscopic, magnetometric, and computational studies confirmed it to be an S = 1 FeIV complex. Complex 1-imido reacts with 2 a to produce the ferrous aniline adduct [(tBu pyrpyrr2 )Fe{NH(2,6-Me2 -4-tBu-C6 H2 )(C14 H11 )}] (1-aniline) in 45 % yield. Lastly, it was found that complexes 1-imido and 1-aniline are both competent intermediates in catalytic intermolecular C-H amination.

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.

Investigating Metal-Metal Bond Polarization in a Heteroleptic Tris-Ylide Diiron System

This article describes the synthesis, characterization, and S-atom transfer reactivity of a series of C3v-symmetric diiron complexes. The iron centers in each complex are coordinated in distinct ligand environments, with one (FeN) bound in a pseudo-trigonal bipyramidal geometry by three phosphinimine nitrogens in the equatorial plane, a tertiary amine, and the second metal center (FeC). FeC is coordinated, in turn, by FeN, three ylidic carbons in a trigonal plane, and, in certain cases, by an axial oxygen donor. The three alkyl donors at FeC form through the reduction of the appended N═PMe3 arms of the monometallic parent complex. The complexes were studied crystallographically, spectroscopically (NMR, UV-vis, and Mössbauer), and computationally (DFT, CASSCF) and found to be high-spin throughout, with short Fe-Fe distances that belie weak orbital overlap between the two metals. Further, the redox nature of this series allowed for the determination that oxidation is localized to the FeC. S-atom transfer chemistry resulted in the formal insertion of a S atom into the Fe-Fe bond of the reduced diiron complex to form a mixture of Fe4S and Fe4S2 products.

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.

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

The HD Reaction of Nitrogenase: a Detailed Mechanism

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

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 P¹⁺/P^N redox potential is nearly the same independently of whether P¹⁺ 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 E₀ state is MoIIIFe3IIFe4III , in agreement with previous experimental studies and QM calculations. Moreover, the redox potentials of the first five E₀-E₄ 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 (E₄-E₈) states of the reaction mechanism have redox potential that are more positive (i.e., more exothermic) than that of the E₀/E₁ couple.

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.

Stability and bonding of carbon(0)-iron-N2 complexes relevant to nitrogenase co-factor: EDA-NOCV analyses

The factors/structural features which are responsible for the binding, activation and reduction of N₂ to NH₃ by FeMoco of nitrogenase have not been completely understood well. Several relevant model complexes by Holland et al. and Peters et al. have been synthesized, characterized and studied by theoretical calculations. For a matter of fact, those complexes are much different than real active N₂ -binding Fe-sites of FeMoco, which possesses a central C(4-) ion having an eight valence electrons as an μ₆ -bridge. Here, a series of [(S₃ C(0))Fe(II/I/0)-N₂ ]^n- complexes in different charged/spin states containing a coordinated σ- and π-donor C(0)-atom which possesses eight outer shell electrons [carbone, (Ph₃ P)₂ C(0); Ph₃ P→C(0)←PPh₃ ] and three S-donor sites (i.e. ^- S-Ar), have been studied by DFT, QTAIM, and EDA-NOCV calculations. The effect of the weak field ligand on Fe-centres and the subsequent N₂ -binding has been studied by EDA-NOCV analysis. The role of the oxidation state of Fe and N₂ -binding in different charged and spin states of the complex have been investigated by EDA-NOCV analyses. The intrinsic interaction energies of the Fe-N₂ bond are in the range from -42/-35 to -67 kcal/mol in their corresponding ground states. The S₃ C(0) donor set is argued here to be closer to the actual coordination environment of one of the six Fe-centres of nitrogenase. In comparison, the captivating model complexes reported by Holland et al. and Peter et al. possess a stronger π-acceptor C-ring (S₂ C_ring donor, π-C donor) and stronger donor set like CP₃ (σ-C donor) ligands, respectively.

Metal-Metal Bond Umpolung in Heterometallic Extended Metal Atom Chains

Understanding the fundamental properties governing metal–metal interactions is crucial to understanding the electronic structure and thereby applications of multimetallic systems in catalysis, material science, and magnetism. One such property that is relatively underexplored within multimetallic systems is metal–metal bond polarity, parameterized by the electronegativities (χ) of the metal atoms involved in the bond. In heterobimetallic systems, metal–metal bond polarity is a function of the donor–acceptor (Δχ) interactions of the two bonded metal atoms, with electropositive early transition metals acting as electron acceptors and electronegative late transition metals acting as electron donors. We show in this work, through the preparation and systematic study of a series of Mo₂M(dpa)₄(OTf)₂ (M = Cr, Mn, Fe, Co, and Ni; dpa = 2,2′-dipyridylamide; OTf = trifluoromethanesulfonate) heterometallic extended metal atom chain (HEMAC) complexes that this expected trend in χ can be reversed. Physical characterization via single-crystal X-ray diffraction, magnetometry, and spectroscopic methods as well as electronic structure calculations supports the presence of a σ symmetry 3c/3e^– bond that is delocalized across the entire metal-atom chain and forms the basis of the heterometallic Mo₂–M interaction. The delocalized 3c/3e^– interaction is discussed within the context of the analogous 3c/3e^– π bonding in the vinoxy radical, CH₂CHO. The vinoxy comparison establishes three predictions for the σ symmetry 3c/3e^– bond in HEMACS: (1) an umpolung effect that causes the Mo–M interactions to become more covalent as Δχ increases, (2) distortion of the σ bonding and non-bonding orbitals to emphasize Mo–M bonding and de-emphasize Mo–Mo bonding, and (3) an increase in Mo spin population with increasing Mo–M covalency. In agreement with these predictions, we find that the Mo₂···M covalency increases with increasing Δχ of the Mo and M atoms (Δχ_Mo–M increases as M = Cr < Mn < Fe < Co < Ni), an umpolung of the trend predicted in the absence of σ delocalization. We attribute the observed trend in covalency to the decreased energic differential (ΔE) between the heterometal orbital and the σ bonding molecular orbital of the Mo₂ quadruple bond, which serves as an energetically stable, “ligand”-like electron-pair donor to the heterometal ion acceptor. As M is changed from Cr to Ni, the σ bonding and nonbonding orbitals do indeed distort as anticipated, and the spin population of the outer Mo group is increased by at least a factor of 2. These findings provide a predictive framework for multimetallic compounds and advance the current understanding of the electronic structures of molecular heteromultimetallic systems, which can be extrapolated to applications in the context of mixed-metal surface catalysis and multimetallic proteins.

This work describes how use of a metal−metal quadruply bonded metalloligand can reverse expected trends in metal−metal bond polarity through the preparation and systematic study of a novel series of Mo₂M(dpa)₄(OTf)₂ (M = Cr, Mn, Fe, Co, and Ni) heterotrimetallic extended metal atom chain (HEMAC) complexes. These complexes feature a 3c/3e⁻ metal−metal bond that is delocalized across the entire metal atom chain and is compared to the 3c/3e⁻ π bonding in the vinoxyl radical.

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.

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

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 r²SCAN 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 r²SCAN 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 r²SCAN, 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.

Thermodynamically Favourable States in the Reaction of Nitrogenase without Dissociation of any Sulfide Ligand

We have used combined quantum mechanical and molecular mechanical (QM/MM) calculations to study the reaction mechanism of nitrogenase, assuming that none of the sulfide ligands dissociates. To avoid the problem that there is no consensus regarding the structure and protonation of the E₄ state, we start from a state where N₂ is bound to the cluster and is protonated to N₂H₂, after dissociation of H₂. We show that the reaction follows an alternating mechanism with HNNH (possibly protonated to HNNH₂) and H₂NNH₂ as intermediates and the two NH₃ products dissociate at the E₇ and E₈ levels. For all intermediates, coordination to Fe6 is preferred, but for the E₄ and E₈ intermediates, binding to Fe2 is competitive. For the E₄, E₅ and E₇ intermediates we find that the substrate may abstract a proton from the hydroxy group of the homocitrate ligand of the FeMo cluster, thereby forming HNNH₂, H₂NNH₂ and NH₃ intermediates. This may explain why homocitrate is a mandatory component of nitrogenase. All steps in the suggested reaction mechanism are thermodynamically favourable compared to protonation of the nearby His‐195 group and in all cases, protonation of the NE2 atom of the latter group is preferred.

Molybdenum-Sulfur Bond: Electronic Structure of Low-Lying States of MoS

The molybdenum-sulfur bond plays an important role in many processes such as nitrogen-fixation, and it is found as a building block in layered materials such as MoS₂, known for its various shapes and morphologies. Here, we present an accurate theoretical and experimental investigation of the chemical bonding and the electronic structure of 20 low-lying states of the MoS molecule. Multireference and coupled cluster methodologies, namely, MRCISD, MRCISD + Q, RCCSD(T), and RCCSD[T], were employed in conjunction with basis sets up to aug-cc-pwCV5Z-PP/aug-cc-pwCV5Z for the study of these states. We note the significance of including the inner 4s²4p⁶ electrons of Mo and 2s²2p⁶ of S in the correlated space to obtain accurate results. Experimentally, the predissociation threshold of MoS was measured using resonant two-photon ionization spectroscopy, allowing for a precise measurement of the bond dissociation energy. Our extrapolated computational D₀ value for the ground state is 3.936 eV, in excellent agreement with our experimental measurement of 3.932 ± 0.004 eV. The largest calculated adiabatic D₀ (5.74 eV) and the largest dipole moment (6.50 D) were found for the ⁵Σ⁺ state, where a triple bond is formed. Finally, the connection of the chemical bonding of the isolated MoS species to the relevant solid, MoS₂, is emphasized. The low-lying septet states of the diatomic molecule are involved in the material as a building block, explaining the stability and the variety of the shapes and morphologies of the material.

A Combined Spectroscopic and Computational Study on the Mechanism of Iron-Catalyzed Aminofunctionalization of Olefins Using Hydroxylamine Derived N-O Reagent as the "Amino" Source and "Oxidant"

Herein, we study the mechanism of iron-catalyzed direct synthesis of unprotected aminoethers from olefins by a hydroxyl amine derived reagent using a wide range of analytical and spectroscopic techniques (Mössbauer, Electron Paramagnetic Resonance, Ultra-Violet Visible Spectroscopy, X-ray Absorption, Nuclear Resonance Vibrational Spectroscopy, and resonance Raman) along with high-level quantum chemical calculations. The hydroxyl amine derived triflic acid salt acts as the “oxidant” as well as “amino” group donor. It activates the high-spin Fe(II) (S_t = 2) catalyst [Fe(acac)₂(H₂O)₂] (1) to generate a high-spin (S_t = 5/2) intermediate (Int I), which decays to a second intermediate (Int II) with S_t = 2. The analysis of spectroscopic and computational data leads to the formulation of Int I as [Fe(III)(acac)₂-N-acyloxy] (an alkyl-peroxo-Fe(III) analogue). Furthermore, Int II is formed by N–O bond homolysis. However, it does not generate a high-valent Fe(IV)(NH) species (a Fe(IV)(O) analogue), but instead a high-spin Fe(III) center which is strongly antiferromagnetically coupled (J = −524 cm^–1) to an iminyl radical, [Fe(III)(acac)₂-NH·], giving S_t = 2. Though Fe(NH) complexes as isoelectronic surrogates to Fe(O) functionalities are known, detection of a high-spin Fe(III)-N-acyloxy intermediate (Int I), which undergoes N–O bond cleavage to generate the active iron–nitrogen intermediate (Int II), is unprecedented. Relative to Fe(IV)(O) centers, Int II features a weak elongated Fe–N bond which, together with the unpaired electron density along the Fe–N bond vector, helps to rationalize its propensity for N-transfer reactions onto styrenyl olefins, resulting in the overall formation of aminoethers. This study thus demonstrates the potential of utilizing the iron-coordinated nitrogen-centered radicals as powerful reactive intermediates in catalysis.

Bonding and stability of dinitrogen-bonded donor base-stabilized Si(0)/Ge(0) species [(cAACMe-Si/Ge)2(N2)]: EDA-NOCV analysis

Recently, dinitrogen (N₂) binding and its activation have been achieved by non-metal compounds like intermediate cAAC-borylene as (cAAC)₂(B-Dur)₂(N₂) [cAAC = cyclic alkyl(amino) carbene; Dur = aryl group, 2,3,5,6-tetramethylphenyl; B-Dur = borylene]. It has attracted a lot of scientific attention from different research areas because of its future prospects as a potent species towards the metal free reduction of N₂ into ammonia (NH₃) under mild conditions. Two (cAAC)(B-Dur) units, each of which possesses six valence electrons around the B-centre, are shown to accept σ-donations from the N₂ ligand (B ← N₂). Two B-Dur further provide π-backdonations (B → N₂) to a central N₂ ligand to strengthen the B–N₂–B bond, providing maximum stability to the compound (cAAC)₂(B-Dur)₂(N₂) since the summation of each pair wise interaction accounted for the total stabilization energy of the molecule. (cAAC)(B-Dur) unit is isolobal to cAAC–E (E = Si, Ge) fragment. Herein, we report on the stability and bonding of cAAC–E bonded N₂-complex (cAAC–E)₂(N₂) (1–2; Si, Ge) by NBO, QTAIM and EDA-NOCV analyses (EDA-NOCV = energy decomposition analysis coupled with natural orbital for chemical valence; QTAIM = quantum theory of atoms in molecule). Our calculation suggested that syntheses of elusive (cAAC–E)₂(N₂) (1–2; Si, Ge) species may be possible with cAAC ligands having bulky substitutions adjacent to the C_cAAC atom by preventing the homo-dimerization of two (cAAC)(E) units which can lead to the formation of (cAAC–E)₂. The formation of E Created by potrace 1.16, written by Peter Selinger 2001-2019 ]]> E bond is thermodynamically more favourable (E = Si, Ge) over binding energy of N₂ inbetween two cAAC–E units.

Dinitrogen (N₂) binding and its activation have been achieved by non-metal compounds like intermediate cAACborylene with the general formula of (cAAC)₂(B-Dur)₂(N₂) [cAAC = cyclic alkyl(amino)carbene; Dur = aryl group, 2,3,5,6-tetramethylphenyl; B-Dur = aryl-borylene].

Stabilization of intermediate spin states in mixed-valent diiron dichalcogenide complexes

The electronic structure and ground spin states, S, observed for mixed-valent iron–sulfur dimers (Fe^II-Fe^III) are typically determined by the Heisenberg exchange interaction, J, that couples the magnetic interaction of the two metal centres either ferromagnetically (J > 0, S = 9/2) or antiferromagnetically (J < 0, S = 1/2). In the case of antiferromagnetically coupled iron centres, stabilization of the high-spin S = 9/2 ground state is also feasible through a Heisenberg double-exchange interaction, B, which lifts the degeneracy of the Heisenberg spin states. This theorem also predicts intermediate spin states for mixed-valent dimers, but those have so far remained elusive. Herein, we describe the structural, electron paramagnetic resonance and Mössbauer spectroscopic, and magnetic characterization of a series of mixed-valent complexes featuring [Fe₂Q₂]⁺ (Q = S^2–, Se^2–, Te^2–), where the Se and Te complexes favour S = 3/2 spin states. The incorporation of heavier chalcogenides in this series reveals a delicate balance of antiferromagnetic coupling, Heisenberg double-exchange and vibronic coupling.

Despite extensive investigations of mixed-valence complexes, molecules with intermediate spin states have remained elusive. Now, selenium- and tellurium-bridged mixed-valent iron dimers have been prepared in which a balance of Heisenberg exchange and double-exchange coupling of the unpaired electron, combined with moderate vibronic contributions, stabilizes S = 3/2 ground spin states.

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.

Stabilization of a high-spin three-coordinate Fe(III) imidyl complex by radical delocalization

High-spin, late transition metal imido complexes have attracted significant interest due to their group transfer reactivity and catalytic C−H activation of organic substrates. Reaction of a new two-coordinate iron complex,...

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.

Estimations of Fe-N2 Intrinsic Interaction Energies of Iron-Sulfur/Nitrogen-Carbon Sites: A Deeper Bonding Insight by EDA-NOCV Analysis of a Model Complex of the Nitrogenase Cofactor

The MoFe₇S₉C^1- unit of the nitrogenase cofactor (FeMoco) attracts chemists and biochemists due to its unusual ability to bind aerial dinitrogen (N₂) at ambient condition and catalytically convert it into ammonia (NH₃). The mode of N₂ binding and its reaction pathways are yet not clear. An important conclusion has been made based on the very recent synthesis and isolation of model Fe(I/0)-complexes with sulfur-donor ligands under the cleavage of one Fe-S bond followed by binding of N₂ at the Fe(0) center. These complexes are structurally relevant to the nitrogenase cofactor (MoFe₇S₉C^1-). Herein, we report the EDA-NOCV analyses and NICS calculations of the dinitrogen-bonded dianionic complex Fe⁰-N₂ (1) (having a C_Ar ← Fe π-bond) and monoanionic complex Fe^I-N₂ (2) (having a C_Ar-Fe σ-bond) to provide a deeper insight into the Fe-N₂ interacting orbitals and corresponding pairwise interaction energies (EDA-NOCV = energy decomposition analysis coupled with natural orbital for chemical valence; NICS = nucleus-independent chemical shifts). The orbital interaction in the Fe-N₂ bond is significantly larger than Coulombic interactions, with major pairwise contributions coming from d(Fe) orbitals to the empty π* orbitals of N₂ (three Fe → N₂). ΔE _int values are in the range of -61 to -77 kcal mol^-1. Very interestingly, NICS calculations have been carried out for the fragments before and after binding of the N₂ molecule. The computed σ- and π-aromaticity values are attributed to the position of the Fe atoms, oxidation states of Fe centers, and Fe-C bond lengths of these two complexes.

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.

Halogen Transfer to Carbon Radicals by High-Valent Iron Chloride and Iron Fluoride Corroles

High-valent iron halide corroles were examined to determine their reactivity with carbon radicals and their ability to undergo radical rebound-like processes. Beginning with Fe(Cl)(ttppc) (1) (ttppc = 5,10,15-tris(2,4,6-triphenylphenyl)corrolato^3-), the new iron corroles Fe(OTf)(ttppc) (2), Fe(OTf)(ttppc)(AgOTf) (3), and Fe(F)(ttppc) (4) were synthesized. Complexes 3 and 4 are the first iron triflate and iron fluoride corroles to be structurally characterized by single crystal X-ray diffraction. The structure of 3 reveals an Ag^I-pyrrole (η²-π) interaction. The Fe(Cl)(ttppc) and Fe(F)(ttppc) complexes undergo halogen transfer to triarylmethyl radicals, and kinetic analysis of the reaction between (p-OMe-C₆H₄)₃C• and 1 gave k = 1.34(3) × 10³ M^-1 s^-1 at 23 °C and 2.2(2) M^-1 s^-1 at -60 °C, ΔH^⧧ = +9.8(3) kcal mol^-1, and ΔS^⧧ = -14(1) cal mol^-1 K^-1 through an Eyring analysis. Complex 4 is significantly more reactive, giving k = 1.16(6) × 10⁵ M^-1 s^-1 at 23 °C. The data point to a concerted mechanism and show the trend X = F^- > Cl^- > OH^- for Fe(X)(ttppc). This study provides mechanistic insights into halogen rebound for an iron porphyrinoid complex.

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.

Resolution of Low-Energy States in Spin-Exchange Transition-Metal Clusters: Case Study of Singlet States in [Fe(III)4S4] Cubanes

Polynuclear transition-metal (PNTM) clusters owe their catalytic activity to numerous energetically low-lying spin states and stable oxidation states. The characterization of their electronic structure represents one of the greatest challenges of modern chemistry. We propose a theoretical framework that enables the resolution of targeted electronic states with ease and apply it to two [Fe(III)₄S₄] cubanes. Through direct access to their many-body wave functions, we identify important correlation mechanisms and their interplay with the geometrical distortions observed in these clusters, which are core properties in understanding their catalytic activity. The simulated magnetic coupling constants predicted by our strategy allow us to make qualitative connections between spin interactions and geometrical distortions, demonstrating its predictive power. Moreover, despite its simplicity, the strategy provides magnetic coupling constants in good agreement with the available experimental ones. The complexes are intrinsically frustrated anti-ferromagnets, and the obtained spin structures together with the geometrical distortions represent two possible ways to release spin frustration (spin-driven Jahn–Teller distortion). Our paradigm provides a simple, yet rigorous, route to uncover the electronic structure of PNTM clusters and may be applied to a wide variety of such clusters.

Structure-Spectroscopy Correlations for Intermediate Q of Soluble Methane Monooxygenase: Insights from QM/MM Calculations

The determination of the diiron core intermediate structures involved in the catalytic cycle of soluble methane monooxygenase (sMMO), the enzyme that selectively catalyzes the conversion of methane to methanol, has been a subject of intense interest within the bioinorganic scientific community. Particularly, the specific geometry and electronic structure of the intermediate that precedes methane binding, known as intermediate Q (or MMOH_Q), has been debated for over 30 years. Some reported studies support a bis-μ-oxo-bridged Fe(IV)₂O₂ closed-core conformation Fe(IV)₂O₂ core, whereas others favor an open-core geometry, with a longer Fe–Fe distance. The lack of consensus calls for a thorough re-examination and reinterpretation of the spectroscopic data available on the MMOH_Q intermediate. Herein, we report extensive simulations based on a hybrid quantum mechanics/molecular mechanics approach (QM/MM) approach that takes into account the complete enzyme to explore possible conformations for intermediates MMOH_ox and MMOH_Q of the sMMOH catalytic cycle. High-level quantum chemical approaches are used to correlate specific structural motifs with geometric parameters for comparison with crystallographic and EXAFS data, as well as with spectroscopic data from Mössbauer spectroscopy, Fe K-edge high-energy resolution X-ray absorption spectroscopy (HERFD XAS), and resonance Raman ¹⁶O–¹⁸O difference spectroscopy. The results provide strong support for an open-core-type configuration in MMOH_Q, with the most likely topology involving mono-oxo-bridged Fe ions and alternate terminal Fe-oxo and Fe-hydroxo groups that interact via intramolecular hydrogen bonding. The implications of an open-core intermediate Q on the reaction mechanism of sMMO are discussed.