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

Electronic isomerism in a heterometallic nickel-iron-sulfur cluster models substrate binding and cyanide inhibition of carbon monoxide dehydrogenase

The nickel-iron carbon monoxide dehydrogenase (CODH) enzyme uses a heterometallic nickel-iron-sulfur ([NiFe4S4]) cluster to catalyze the reversible interconversion of carbon dioxide (CO2) and carbon monoxide (CO). These reactions are essential for maintaining the global carbon cycle and offer a route towards sustainable greenhouse gas conversion but have not been successfully replicated in synthetic models, in part due to a poor understanding of the natural system. Though the general protein architecture of CODH is known, the electronic structure of the active site is not well-understood, and the mechanism of catalysis remains unresolved. To better understand the CODH enzyme, we have developed a protein-based model containing a heterometallic [NiFe3S4] cluster in the Pyrococcus furiosus (Pf) ferredoxin (Fd). This model binds small molecules such as carbon monoxide and cyanide, analogous to CODH. Multiple redox- and ligand-bound states of [NiFe3S4] Fd (NiFd) have been investigated using a suite of spectroscopic techniques, including resonance Raman, Ni and Fe K-edge X-ray absorption spectroscopy, and electron paramagnetic resonance, to resolve charge and spin delocalization across the cluster, site-specific electron density, and ligand activation. The facile movement of charge through the cluster highlights the fluidity of electron density within iron-sulfur clusters and suggests an electronic basis by which CN- inhibits the native system while the CO-bound state continues to elude isolation in CODH. The detailed characterization of isolable states that are accessible in our CODH model system provides valuable insight into unresolved enzymatic intermediates and offers design principles towards developing functional mimics of CODH.

Three-Coordinate Nickel and Metal-Metal Interactions in a Heterometallic Iron-Sulfur Cluster

Biological multielectron reactions often are performed by metalloenzymes with heterometallic sites, such as anaerobic carbon monoxide dehydrogenase (CODH), which has a nickel-iron-sulfide cubane with a possible three-coordinate nickel site. Here, we isolate the first synthetic iron-sulfur clusters having a nickel atom with only three donors, showing that this structural feature is feasible. These have a core with two tetrahedral irons, one octahedral tungsten, and a three-coordinate nickel connected by sulfide and thiolate bridges. Electron paramagnetic resonance (EPR), Mössbauer, and superconducting quantum interference device (SQUID) data are combined with density functional theory (DFT) computations to show how the electronic structure of the cluster arises from strong magnetic coupling between the Ni, Fe, and W sites. X-ray absorption spectroscopy, together with spectroscopically validated DFT analysis, suggests that the electronic structure can be described with a formal Ni1+ atom participating in a nonpolar Ni-W σ-bond. This metal-metal bond, which minimizes spin density at Ni1+, is conserved in two cluster oxidation states. Fe-W bonding is found in all clusters, in one case stabilizing a local non-Hund state at tungsten. Based on these results, we compare different M-M interactions and speculate that other heterometallic clusters, including metalloenzyme active sites, could likewise store redox equivalents and stabilize low-valent metal centers through metal-metal bonding.

Valence Delocalization and Metal-Metal Bonding in Carbon-Bridged Mixed-Valence Iron Complexes

The carbide ligand in the iron-molybdenum cofactor (FeMoco) in nitrogenase bridges iron atoms in different oxidation states, yet it is difficult to discern its ability to mediate magnetic exchange interactions due to the structural complexity of the cofactor. Here, we describe two mixed-valent diiron complexes with C-based ketenylidene bridging ligands, and compare the carbon bridges with the more familiar sulfur bridges. The ground state of the [Fe2 (μ-CCO)2 ]+ complex with two carbon bridges (4) is S=1 / 2 ${{ 1/2 }}$ , and it is valence delocalized on the Mössbauer timescale with a small thermal barrier for electron hopping that stems from the low Fe-C force constant. In contrast, one-electron reduction of the [Fe2 (μ-CCO)] complex with one carbon bridge (2) affords a mixed-valence species with a high-spin ground state (S=7 / 2 ${ 7/2 }$ ), and the Fe-Fe distance contracts by 1 Å. Spectroscopic, magnetic, and computational studies of the latter reveal an Fe-Fe bonding interaction that leads to complete valence delocalization. Analysis of near-IR intervalence charge transfer transitions in 5 indicates a very large double exchange constant (B) in the range of 780-965 cm-1 . These results show that carbon bridges are extremely effective at stabilizing valence delocalized ground states in mixed-valent iron dimers.

Isostructural bridging diferrous chalcogenide cores [FeII(μ-E)FeII] (E = O, S, Se, Te) with decreasing antiferromagnetic coupling down the chalcogenide series

Iron compounds containing a bridging oxo or sulfido moiety are ubiquitous in biological systems, but substitution with the heavier chalcogenides selenium and tellurium, however, is much rarer, with only a few examples reported to date. Here we show that treatment of the ferrous starting material [(^tBupyrpyrr₂)Fe(OEt₂)] (1-OEt₂) (^tBupyrpyrr₂ = 3,5-^tBu₂-bis(pyrrolyl)pyridine) with phosphine chalcogenide reagents E = PR₃ results in the neutral phosphine chalcogenide adduct series [(^tBupyrpyrr₂)Fe(EPR₃)] (E = O, S, Se; R = Ph; E = Te; R = ^tBu) (1-E) without any electron transfer, whereas treatment of the anionic starting material [K]₂[(^tBupyrpyrr₂)Fe₂(μ-N₂)] (2-N₂) with the appropriate chalcogenide transfer source yields cleanly the isostructural ferrous bridging mono-chalcogenide ate complexes [K]₂[(^tBupyrpyrr₂)Fe₂(μ-E)] (2-E) (E = O, S, Se, and Te) having significant deviation in the Fe-E-Fe bridge from linear in the case of E = O to more acute for the heaviest chalcogenide. All bridging chalcogenide complexes were analyzed using a variety of spectroscopic techniques, including ¹H NMR, UV-Vis electronic absorbtion, and ⁵⁷Fe Mössbauer. The spin-state and degree of communication between the two ferrous ions were probed via SQUID magnetometry, where it was found that all iron centers were high-spin (S = 2) Fe^II, with magnetic exchange coupling between the Fe^II ions. Magnetic studies established that antiferromagnetic coupling between the ferrous ions decreases as the identity of the chalcogen is tuned from O to the heaviest congener Te.

Exploiting Molecular Symmetry to Quantitatively Map the Excited-State Landscape of Iron-Sulfur Clusters

Cuboidal [Fe₄S₄] clusters are ubiquitous cofactors in biological redox chemistry. In the [Fe₄S₄]¹⁺ state, pairwise spin coupling gives rise to six arrangements of the Fe valences ("valence isomers") among the four Fe centers. Because of the magnetic complexity of these systems, it has been challenging to understand how a protein's active site dictates both the arrangement of the valences in the ground state as well as the population of excited-state valence isomers. Here, we show that the ground-state valence isomer landscape can be simplified from a six-level system in an asymmetric protein environment to a two-level system by studying the problem in synthetic [Fe₄S₄]¹⁺ clusters with solution C_3v symmetry. This simplification allows for the energy differences between valence isomers to be quantified (in some cases with a resolution of <0.1 kcal/mol) by simultaneously fitting the VT NMR and solution magnetic moment data. Using this fitting protocol, we map the excited-state landscape for a range of clusters of the form [(SIMes)₃Fe₄S₄-X/L]ⁿ, (SIMes = 1,3-dimesityl-imidazol-4,5-dihydro-2-ylidene; n = 0 for anionic, X-type ligands and n = +1 for neutral, L-type ligands) and find that a single ligand substitution can alter the relative ground-state energies of valence isomers by at least 10³ cm^-1. On this basis, we suggest that one result of "non-canonical" amino acid ligation in Fe-S proteins is the redistribution of the valence electrons in the manifold of thermally populated excited states.

Structural insight into halide-coordinated [Fe4S4XnY4-n]2- clusters (X, Y = Cl, Br, I) by XRD and Mössbauer spectroscopy

Iron sulphur halide clusters [Fe₄S₄Br₄]^2- and [Fe₄S₄X₂Y₂]^2- (X, Y = Cl, Br, I) were obtained in excellent yields (77 to 78%) and purity from [Fe(CO)₅], elemental sulphur, I₂ and benzyltrimethylammonium (BTMA⁺) iodide, bromide and chloride. Single crystals of (BTMA)₂[Fe₄S₄Br₄] (1), (BTMA)₂[Fe₄S₄Br₂Cl₂] (2), (BTMA)₂[Fe₄S₄Cl₂I₂] (3), and (BTMA)₂[Fe₄S₄Br₂I₂] (4) were isostructural to the previously reported (BTMA)₂[Fe₄S₄I₄] (5) (monoclinic, Cc). Instead of the chloride cubane cluster [Fe₄S₄Cl₄]^2-, we found the prismane-shaped cluster (BTMA)₃[Fe₆S₆Cl₆] (6) (P1̄). ⁵⁷Fe Mössbauer spectroscopy indicates complete delocalisation with Fe^2.5+ oxidation states for all iron atoms. Magnetic measurements showed small χ_MT values at 298 K ranging from 1.12 to 1.54 cm³ K mol^-1, indicating the dominant antiferromagnetic exchange interactions. With decreasing temperature, the χ_MT values decreased to reach a plateau at around 100 K. From about 20 K, the values drop significantly. Fitting the data in the Heisenberg-Dirac-van Vleck (HDvV) as well as the Heisenberg Double Exchange (HDE) formalism confirmed the delocalisation and antiferromagnetic coupling assumed from Mössbauer spectroscopy.

Magnetic exchange and valence delocalization in a mixed valence [Fe2+Fe3+Te2]+ complex: insights from theory and interpretations of magnetic and spectroscopic data

A mixed valence binuclear Fe^2.5+-Fe^2.5+ (Robin-Day Class III) transition metal complex, [Fe^2.5+μTe₂Fe^2.5+]^1-, composed of two FeN₂Te₂ pseudo-tetrahedral units with μ-bridging Te^2- ligands was reported to exist in an unprecedented S = 3/2 ground state (Nature Chemistry, https://doi.org/10.1038/s41557-021-00853-5). For this and the homologous complexes containing Se^2- and S^2-, the Anderson-Hasegawa double exchange spin-Hamiltonian was broadly used to interpret the corresponding structural, spectroscopic and magnetic data. First principles multireference ab initio calculations are used here to simulate magnetic and spectroscopic EPR data; analysis of the results affords a rationale for the stabilization of the S = 3/2 ground state of the Fe₂ pair. Complete Active Space Self-Consistent Field (CASSCF) calculations and dynamical correlation accounted for by means of N-Electron Valence Perturbation Theory to Second Order (NEVPT2) reproduce well the g-factors determined from simulations of X-band EPR spectra. A crucial technical tool to achieve these results is: (i) use of a localized orbital formulation of the many-particle problem at the scalar-relativistic CASSCF step; (ii) choice of state averaging over states of a given spin (at the CASCI/NEVPT2 step); and (iii) accounting for spin-orbit coupling within the non-relativistic Born-Oppenheimer (BO) many-particle basis using Quasi-Degenerate Perturbation Theory (QDPT). The inclusion of the S = 5/2 spin manifold reproduced the observed increase in the magnetic susceptibility (χT) in the high temperature range (T > 100 K), which is explained by thermal population of the S = 5/2 excited state at energy 160 cm^-1 above the S = 3/2 ground state. Theoretical values of χT from experimentally reported data points in the temperature range from 3 to 30 K were further computed and analyzed using a model which takes spin-phonon coupling into account. The model considerations and the computational protocols of this study are generally applicable to any Class I/II mixed valence dimer. The work can potentially stimulate further experimental and theoretical work on bi- and oligonuclear transition metal complexes of importance to bioinorganic chemistry and life sciences.

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.

Selenium Valence-to-Core X-ray Emission Spectroscopy and Kβ HERFD X-ray Absorption Spectroscopy as Complementary Probes of Chemical and Electronic Structure

Selenium X-ray absorption spectroscopy (XAS) has found widespread use in investigations of Se-containing materials, geochemical processes, and biologically active sites. In contrast to sulfur Kβ X-ray emission spectroscopy (XES), which has been found to contain electronic and structural information complementary to S XAS, Se Kβ XES remains comparatively underexplored. Herein, we present the first Se Valence-to-Core (VtC) XES studies of reduced Se-containing compounds and FeSe dimers. Se VtC XES is found to be sensitive to changes in covalent Se bonding interactions (Se–Se/Se–C/Se–H bonding) while being relatively insensitive to changes in Fe oxidation states as selenide bridges in FeSe dimers ([Fe₂Se₂]²⁺ vs [Fe₂Se₂]⁺). In contrast, Se Kβ HERFD XAS is demonstrated to be quite sensitive to changes in the Fe oxidation state with Se Kβ HERFD XAS demonstrating experimental resolution equivalent to Kα HERFD XAS. Additionally, computational studies reveal both Se VtC XES and XAS to be sensitive to selenium protonation in FeSe complexes.

Selenium is a trace element that plays vital roles in biological and geochemical cycles, energy storage, photovoltaics, and nanomaterials. Herein, selenium Valence-to-Core X-ray emission spectroscopy is explored as a new method of probing the chemical and electronic structure in selenium-containing compounds, demonstrating sensitivity to selenium bonding interactions. When paired with high-resolution Se X-ray absorption spectroscopy (HERFD XAS), these two methods have the potential to reveal greater insight into protonation and redox changes of Se-substituted FeS clusters.