Tuning the electronic structure of octahedral iron complexes [FeL(X)] (L = 1-alkyl-4,7-bis(4-tert-butyl-2-mercaptobenzyl)-1,4,7-triazacyclononane, X = Cl, CH(3)O, CN, NO). The S = 1/2 <==>3/2 Spin equilibrium of [FeL(Pr)(NO)]

Cyanide replaces substrate in obligate-ordered addition of nitric oxide to the non-heme mononuclear iron AvMDO active site

Thiol dioxygenases are a subset of non-heme mononuclear iron oxygenases that catalyze the O₂-dependent oxidation of thiol-bearing substrates to yield sulfinic acid products. Cysteine dioxygenase (CDO) and 3-mercaptopropionic acid (3MPA) dioxygenase (MDO) are the most extensively characterized members of this enzyme family. As with many non-heme mononuclear iron oxidase/oxygenases, CDO and MDO exhibit an obligate-ordered addition of organic substrate before dioxygen. As this substrate-gated O₂-reactivity extends to the oxygen-surrogate, nitric oxide (NO), EPR spectroscopy has long been used to interrogate the [substrate:NO:enzyme] ternary complex. In principle, these studies can be extrapolated to provide information about transient iron-oxo intermediates produced during catalytic turnover with dioxygen. In this work, we demonstrate that cyanide mimics the native thiol-substrate in ordered-addition experiments with MDO cloned from Azotobacter vinelandii (AvMDO). Following treatment of the catalytically active Fe(II)-AvMDO with excess cyanide, addition of NO yields a low-spin (S = 1/2) (CN/NO)-Fe-complex. Continuous wave and pulsed X-band EPR characterization of this complex produced in wild-type and H157N variant AvMDO reveal multiple nuclear hyperfine features diagnostic of interactions within the first- and outer-coordination sphere of the enzymatic Fe-site. Spectroscopically validated computational models indicate simultaneous coordination of two cyanide ligands replaces the bidentate (thiol and carboxylate) coordination of 3MPA allowing for NO-binding at the catalytically relevant O₂-binding site. This promiscuous substrate-gated reactivity of AvMDO with NO provides an instructive counterpoint to the high substrate-specificity exhibited by mammalian CDO for L-cysteine.

Spectroscopic analysis of the mammalian enzyme cysteine dioxygenase

l-Cysteine (Cys) is an essential building block for the synthesis of new proteins and serves as a precursor for several biologically important sulfur-containing molecules, such as coenzyme A, taurine, glutathione, and inorganic sulfate. However, organisms must tightly regulate the concentration of free Cys, as elevated levels of this semi-essential amino acid can be extremely harmful. The non-heme iron enzyme cysteine dioxygenase (CDO) serves to maintain the proper levels of Cys by catalyzing its oxidation to cysteine sulfinic acid. Crystal structures of resting and substrate-bound mammalian CDO revealed two surprising structural motifs in the first and second coordination spheres of the Fe center. The first is the existence of a neutral three histidine (3-His) facial triad that coordinates the Fe ion, as opposed to an anionic 2-His-1-carboxylate facial triad that is typically observed in mononuclear non-heme Fe(II) dioxygenases. The second unusual structural feature exhibited by mammalian CDO is the presence of a covalent crosslink between the sulfur of a Cys residue and an ortho-carbon of a tyrosine residue. Spectroscopic studies of CDO have provided invaluable insights into the roles that these unusual features play with regards to substrate Cys and co-substrate O2 binding and activation. In this chapter, we summarize results obtained from electronic absorption, electron paramagnetic resonance, magnetic circular dichroism, resonance Raman, and Mössbauer spectroscopic studies of mammalian CDO carried out in the last two decades. Pertinent results obtained from complementary computational studies are also briefly summarized.

A Density Matrix Renormalization Group Study of the Low-Lying Excited States of a Molybdenum Carbonyl-Nitrosyl Complex

A density matrix renormalization group-self consistent field (DMRG-SCF) study has been carried out to calculate the low-lying excited states of CpMo(CO)₂ NO, a molybdenum complex containing NO and CO ligands. In order to automatically select an appropriate active space, a novel procedure employing the maximum single-orbital entropy for several states has been introduced and shown to be efficient and easy-to-implement when several electronic states are simultaneously considered. The analysis of the resulting natural transition orbitals and charge-transfer numbers shows that the lowest five excited electronic states are excitation into metal-NO antibonding orbitals, which offer the possibility for nitric oxide (NO) photorelease after excitation with visible light. Higher excited states are metal-centered excitations with contributions of metal-CO antibonding orbitals, which may serve as a gateway for carbon monoxide (CO) delivery. Time-dependent density functional theory calculations done for comparison, show that the state characters agree remarkably well with those from DMRG-SCF, while excitation energies are 0.4-1.0 eV red-shifted with respect to the DMRG-SCF ones.

Temperature-Dependent Reactivity of a Non-heme FeIII(OH)(SR) Complex: Relevance to Isopenicillin N Synthase

Non-heme iron complexes with cis-Fe^III(OH)(SAr/OAr) coordination were isolated and examined for their reactivity with a tertiary carbon radical. The sulfur-ligated complex shows a temperature dependence on ^•OH versus ArS^• transfer, whereas the oxygen-ligated complex does not. These results provide the first working model for C-S bond formation in isopenicillin N synthase and indicate that kinetic control may be a key factor in the selectivity of non-heme iron "rebound" processes.

57Fe Mössbauer parameters from domain based local pair-natural orbital coupled-cluster theory

We report on applications of the domain based local pair-natural orbital (PNO) coupled-cluster method within the singles and doubles approximation (DLPNO-CCSD) to the calculation of ⁵⁷Fe isomer shifts and quadrupole splittings in a small training set of iron complexes consisting of large molecular ligands and iron atoms in varying charge, spin, and oxidation states. The electron densities and electric field gradients needed for these calculations were obtained within the recently implemented analytic derivative scheme. A method for the direct treatment of scalar relativistic effects in the calculation of effective electron densities is described by using the first-order Douglas-Kroll-Hess Hamiltonian and a Gaussian charge distribution model for the nucleus. The performance of DLPNO-CCSD is compared with four modern-day density functionals, namely, RPBE, TPSS, B3LYP, and B2PLYP, as well as with the second-order Møller-Plesset perturbation theory. An excellent correlation between the calculated electron densities and the experimental isomer shifts is attained with the DLPNO-CCSD method. The correlation constant a obtained from the slope of the linear correlation plot is found to be ≈-0.31 a.u.³ mm s^-1, which agrees very well with the experimental calibration constant α = -0.31 ± 0.04 a.u.³ mm s^-1. This value of a is obtained consistently using both nonrelativistic and scalar relativistic DLPNO-CCSD electron densities. While the B3LYP and B2PLYP functionals achieve equally good correlation between theory and experiment, the correlation constant a is found to deviate from the experimental value. Similar trends are observed also for quadrupole splittings. The value of the nuclear quadrupole moment for ⁵⁷Fe is estimated to be 0.15 b at the DLPNO-CCSD level. This is consistent with previous results and is here supported by a higher level of theory. The DLPNO-CCSD results are found to be insensitive to the intrinsic approximations in the method, in particular the PNO occupation number truncation error, while the results obtained with density functional theory (DFT) are found to depend on the choice of the functional. In a statistical sense, i.e., on the basis of the linear regression analysis, however, the accuracies of the DFT and DLPNO-CCSD results can be considered comparable.

A reagent for heteroatom borylation, including iron mediated reductive deoxygenation of nitrate yielding a dinitrosyl complex

4,4'-Bipyridyl is shown to be a catalyst for transfer of pinacolboryl groups from (Bpin)₂ to nitrogen heterocycles and to Me₃SiN₃. Using stoichiometric (Bpin)₂(pyrazine) or (Bpin)₂(bipyridine) in an analogous manner, an aromatic nitro group is deoxygenated and subsequently borylated, and four-fold deoxygenation of (DIM)Fe(NO₃)₂(MeCN) to yield the dinitrosyl complex (DIM)Fe(NO)₂ is facile. The co-product O(Bpin)₂ is the quantitative fate of the removed oxo groups. With borylation of both nitrogen heterocycles and doubly deoxygenating two nitrates coordinated to a single metal center, broad spectrum methodology is demonstrated.

Outer-Sphere Tyrosine 159 within the 3-Mercaptopropionic Acid Dioxygenase S-H-Y Motif Gates Substrate-Coordination Denticity at the Non-Heme Iron Active Site

Thiol dioxygenases are non-heme mononuclear iron enzymes that catalyze the O₂-dependent oxidation of free thiols (-SH) to produce the corresponding sulfinic acid (-SO₂^-). Regardless of the phylogenic domain, the active site for this enzyme class is typically comprised of two major features: (1) a mononuclear ferrous iron coordinated by three protein-derived histidines and (2) a conserved sequence of outer Fe-coordination-sphere amino acids (Ser-His-Tyr) spatially adjacent to the iron site (∼3 Å). Here, we utilize a promiscuous 3-mercaptopropionic acid dioxygenase cloned from Azotobacter vinelandii (Av MDO) to explore the function of the conserved S-H-Y motif. This enzyme exhibits activity with 3-mercaptopropionic acid (3mpa), l-cysteine (cys), as well as several other thiol-bearing substrates, thus making it an ideal system to study the influence of residues within the highly conserved S-H-Y motif (H157 and Y159) on substrate specificity and reactivity. The pK_a values for these residues were determined by pH-dependent steady-state kinetics, and their assignments verified by comparison to H157N and Y159F variants. Complementary electron paramagnetic resonance and Mössbauer studies demonstrate a network of hydrogen bonds connecting H157-Y159 and Fe-bound ligands within the enzymatic Fe site. Crucially, these experiments suggest that the hydroxyl group of Y159 hydrogen bonds to Fe-bound NO and, by extension, Fe-bound oxygen during native catalysis. This interaction alters both the NO binding affinity and rhombicity of the 3mpa-bound iron-nitrosyl site. In addition, Fe coordination of cys is switched from thiolate only to bidentate (thiolate/amine) for the Y159F variant, indicating that perturbations within the S-H-Y proton relay network also influence cys Fe binding denticity.

An efficient, step-economical strategy for the design of functional metalloproteins

The bottom-up design and construction of functional metalloproteins remains a formidable task in biomolecular design. Although numerous strategies have been used to create new metalloproteins, pre-existing knowledge of the tertiary and quaternary protein structure is often required to generate suitable platforms for robust metal coordination and activity. Here we report an alternative and easily implemented approach (metal active sites by covalent tethering or MASCoT) in which folded protein building blocks are linked by a single disulfide bond to create diverse metal coordination environments within evolutionarily naive protein-protein interfaces. Metalloproteins generated using this strategy uniformly bind a wide array of first-row transition metal ions (MnII, FeII, CoII, NiII, CuII, ZnII and vanadyl) with physiologically relevant thermodynamic affinities (dissociation constants ranging from 700 nM for MnII to 50 fM for CuII). MASCoT readily affords coordinatively unsaturated metal centres-including a penta-His-coordinated non-haem Fe site-and well-defined binding pockets that can accommodate modifications and enable coordination of exogenous ligands such as nitric oxide to the interfacial metal centre.

Non-heme High-Spin {FeNO}6-8 Complexes: One Ligand Platform Can Do It All

Heme and non-heme iron-nitrosyl complexes are important intermediates in biology. While there are numerous examples of low-spin heme iron-nitrosyl complexes in different oxidation states, much less is known about high-spin (hs) non-heme iron-nitrosyls in oxidation states other than the formally ferrous NO adducts ({FeNO}⁷ in the Enemark-Feltham notation). In this study, we present a complete series of hs-{FeNO}^6-8 complexes using the TMG₃tren coligand. Redox transformations from the hs-{FeNO}⁷ complex [Fe(TMG₃tren)(NO)]²⁺ to its {FeNO}⁶ and {FeNO}⁸ analogs do not alter the coordination environment of the iron center, allowing for detailed comparisons between these species. Here, we present new MCD, NRVS, XANES/EXAFS, and Mössbauer data, demonstrating that these redox transformations are metal based, which allows us to access hs-Fe(II)-NO^-, Fe(III)-NO^-, and Fe(IV)-NO^- complexes. Vibrational data, analyzed by NCA, directly quantify changes in Fe-NO bonding along this series. Optical data allow for the identification of a "spectator" charge-transfer transition that, together with Mössbauer and XAS data, directly monitors the electronic changes of the Fe center. Using EXAFS, we are also able to provide structural data for all complexes. The magnetic properties of the complexes are further analyzed (from magnetic Mössbauer). The properties of our hs-{FeNO}^6-8 complexes are then contrasted to corresponding, low-spin iron-nitrosyl complexes where redox transformations are generally NO centered. The hs-{FeNO}⁸ complex can further be protonated by weak acids, and the product of this reaction is characterized. Taken together, these results provide unprecedented insight into the properties of biologically relevant non-heme iron-nitrosyl complexes in three relevant oxidation states.

Density Functional Calculations for Prediction of 57Fe Mössbauer Isomer Shifts and Quadrupole Splittings in β-Diketiminate Complexes

The relative ease of Mössbauer spectroscopy and of density functional theory (DFT) calculations encourages the use of Mössbauer parameters as a validation method for calculations, and the use of calculations as a double check on crystallographic structures. A number of studies have proposed correlations between the computationally determined electron density at the iron nucleus and the observed isomer shift, but deviations from these correlations in low-valent iron β-diketiminate complexes encouraged us to determine a new correlation for these compounds. The use of B3LYP/def2-TZVP in the ORCA platform provides an excellent balance of accuracy and speed. We provide here not only this new correlation and a clear guide to its use but also a systematic analysis of the limitations of this approach. We also highlight the impact of crystallographic inaccuracies, DFT model truncation, and spin states, with intent to assist experimentalists to use Mössbauer spectroscopy and calculations together.

Reliable Estimation of Prediction Uncertainty for Physicochemical Property Models

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

Structural and Spectroscopic Characterization of a High-Spin {FeNO}(6) Complex with an Iron(IV)-NO(-) Electronic Structure

Although the interaction of low-spin ferric complexes with nitric oxide has been well studied, examples of stable high-spin ferric nitrosyls (such as those that could be expected to form at typical non-heme iron sites in biology) are extremely rare. Using the TMG3 tren co-ligand, we have prepared a high-spin ferric NO adduct ({FeNO}(6) complex) via electrochemical or chemical oxidation of the corresponding high-spin ferrous NO {FeNO}(7) complex. The {FeNO}(6) compound is characterized by UV/Visible and IR spectroelectrochemistry, Mössbauer and NMR spectroscopy, X-ray crystallography, and DFT calculations. The data show that its electronic structure is best described as a high-spin iron(IV) center bound to a triplet NO(-) ligand with a very covalent iron-NO bond. This finding demonstrates that this high-spin iron nitrosyl compound undergoes iron-centered redox chemistry, leading to fundamentally different properties than corresponding low-spin compounds, which undergo NO-centered redox transformations.

Photoinitiated Reactivity of a Thiolate-Ligated, Spin-Crossover Nonheme {FeNO}(7) Complex with Dioxygen

The nonheme iron complex, [Fe(NO)(N3PyS)]BF4, is a rare example of an {FeNO}(7) species that exhibits spin-crossover behavior. The comparison of X-ray crystallographic studies at low and high temperatures and variable-temperature magnetic susceptibility measurements show that a low-spin S = 1/2 ground state is populated at 0-150 K, while both low-spin S = 1/2 and high-spin S = 3/2 states are populated at T > 150 K. These results explain the observation of two N-O vibrational modes at 1737 and 1649 cm(-1) in CD3CN for [Fe(NO)(N3PyS)]BF4 at room temperature. This {FeNO}(7) complex reacts with dioxygen upon photoirradiation with visible light in acetonitrile to generate a thiolate-ligated, nonheme iron(III)-nitro complex, [Fe(III)(NO2)(N3PyS)](+), which was characterized by EPR, FTIR, UV-vis, and CSI-MS. Isotope labeling studies, coupled with FTIR and CSI-MS, show that one O atom from O2 is incorporated in the Fe(III)-NO2 product. The O2 reactivity of [Fe(NO)(N3PyS)]BF4 in methanol is dramatically different from CH3CN, leading exclusively to sulfur-based oxidation, as opposed to NO· oxidation. A mechanism is proposed for the NO· oxidation reaction that involves formation of both Fe(III)-superoxo and Fe(III)-peroxynitrite intermediates and takes into account the experimental observations. The stability of the Fe(III)-nitrite complex is limited, and decay of [Fe(III)(NO2)(N3PyS)](+) leads to {FeNO}(7) species and sulfur oxygenated products. This work demonstrates that a single mononuclear, thiolate-ligated nonheme {FeNO}(7) complex can exhibit reactivity related to both nitric oxide dioxygenase (NOD) and nitrite reductase (NiR) activity. The presence of the thiolate donor is critical to both pathways, and mechanistic insights into these biologically relevant processes are presented.

Reactivity of a Fe(III)-Bound Methoxide Supported with a Tris(thiolato)phosphine Ligand: Activation of C-Cl Bond in CH2Cl2 by Nucleophilic Attack of a Fe(III)-OCH3 Moiety

Two mononuclear nonheme Fe(III) complexes, [PPh4][Fe(III)(PS3″)(OCH3)] (1) and [PPh4][Fe(III)(PS3″)(Cl)] (2), supported by a tris(benzenethiolato)phosphine derivative PS3″ (PS3″ = P(C6H3-3-Me3Si-2-S)3(3-)) have been synthesized and characterized. The structures resolved from X-ray crystallography show that Fe(III) centers in both complexes adopt distorted trigonal-bipyramidal geometry with a methoxide or a chloride binding in the axial position. The magnetic data for both are consistent with intermediate-spin Fe(III) centers with a C3 symmetry (S = 3/2 ground state). The bound methoxide in 1 is labile and can be replaced by a CH3CN molecule. The forming Fe(III)-CH3CN species can be further reduced by cobaltcene quantitatively to a stable Fe(II)-CH3CN complex, [Fe(PS3″)(CH3CN)](-). One-electron oxidation of 2 by ferrocenium gave a Fe(IV) analogue, [Fe(IV)(PS3″)(Cl)]. Importantly, the Fe(III)-OCH3 moiety in complex 1 acts as a strong nucleophile that activates the C-Cl bond in CH2Cl2, leading to the formation of complex 2 quantitatively. Complex 1 also reacts with other electrophiles, benzyl chloride and benzyl bromide, to generate Fe(III)-X species (X = Cl or Br). The reactions were investigated and monitored by UV-vis-NIR, NMR, and ESI-MS spectroscopies.

The "Gln-Type" Thiol Dioxygenase from Azotobacter vinelandii is a 3-Mercaptopropionic Acid Dioxygenase

Cysteine dioxygenase (CDO) is a non-heme iron enzyme that catalyzes the O2-dependent oxidation of l-cysteine to produce cysteinesulfinic acid. Bacterial CDOs have been subdivided as either "Arg-type" or "Gln-type" on the basis of the identity of conserved active site residues. To date, "Gln-type" enzymes remain largely uncharacterized. It was recently noted that the "Gln-type" enzymes are more homologous with another thiol dioxygenase [3-mercaptopropionate dioxygenase (MDO)] identified in Variovorax paradoxus, suggesting that enzymes of the "Gln-type" subclass are in fact MDOs. In this work, a putative "Gln-type" thiol dioxygenase from Azotobacter vinelandii (Av) was purified to homogeneity and characterized. Steady-state assays were performed using three substrates [3-mercaptopropionic acid (3mpa), l-cysteine (cys), and cysteamine (ca)]. Despite comparable maximal velocities, the "Gln-type" Av enzyme exhibited a specificity for 3mpa (kcat/KM = 72000 M(-1) s(-1)) nearly 2 orders of magnitude greater than those for cys (110 M(-1) s(-1)) and ca (11 M(-1) s(-1)). Supporting X-band electron paramagnetic resonance (EPR) studies were performed using nitric oxide (NO) as a surrogate for O2 binding to confirm obligate-ordered addition of substrate prior to NO. Stoichimetric addition of NO to solutions of 3mpa-bound enzyme quantitatively yields an iron-nitrosyl species (Av ES-NO) with EPR features consistent with a mononuclear (S = (3)/2) {FeNO}(7) site. Conversely, two distinct substrate-bound conformations were observed in Av ES-NO samples prepared with cys and ca, suggesting heterogeneous binding within the enzymatic active site. Analytical EPR simulations are provided to establish the relative binding affinity for each substrate (3map > cys > ca). Both kinetic and spectroscopic results presented here are consistent with 3mpa being the preferred substrate for this enzyme.

Orbital entanglement and CASSCF analysis of the Ru-NO bond in a Ruthenium nitrosyl complex

Multiconfigurational wavefunction analysis and entanglement measures based on von Neumann entropy shed light on the electronic structure of a Ru nitrosyl complex, in particular on the Ru–NO bond.

Complete active space self-consistent field (CASSCF) wavefunctions and an orbital entanglement analysis obtained from a density-matrix renormalisation group (DMRG) calculation are used to understand the electronic structure, and, in particular, the Ru–NO bond of a Ru nitrosyl complex. Based on the configurations and orbital occupation numbers obtained for the CASSCF wavefunction and on the orbital entropy measurements evaluated for the DMRG wavefunction, we unravel electron correlation effects in the Ru coordination sphere of the complex. It is shown that Ru–NO π bonds show static and dynamic correlation, while other Ru–ligand bonds feature predominantly dynamic correlation. The presence of static correlation requires the use of multiconfigurational methods to describe the Ru–NO bond. Subsequently, the CASSCF wavefunction is analysed in terms of configuration state functions based on localised orbitals. The analysis of the wavefunctions in the electronic singlet ground state and the first triplet state provides a picture of the Ru–NO moiety beyond the standard representation based on formal oxidation states. A distinct description of the Ru and NO fragments is advocated. The electron configuration of Ru is an equally weighted superposition of Ru^II and Ru^III configurations, with the Ru^III configuration originating from charge donation mostly from Cl ligands. However, and contrary to what is typically assumed, the electronic configuration of the NO ligand is best described as electroneutral.

Spectroscopic and computational study of a nonheme iron nitrosyl center in a biosynthetic model of nitric oxide reductase

A major barrier to understanding the mechanism of nitric oxide reductases (NORs) is the lack of a selective probe of NO binding to the nonheme FeB center. By replacing the heme in a biosynthetic model of NORs, which structurally and functionally mimics NORs, with isostructural ZnPP, the electronic structure and functional properties of the FeB nitrosyl complex was probed. This approach allowed observation of the first S=3/2 nonheme {FeNO}(7) complex in a protein-based model system of NOR. Detailed spectroscopic and computational studies show that the electronic state of the {FeNO}(7) complex is best described as a high spin ferrous iron (S=2) antiferromagnetically coupled to an NO radical (S=1/2) [Fe(2+)-NO(.)]. The radical nature of the FeB -bound NO would facilitate N-N bond formation by radical coupling with the heme-bound NO. This finding, therefore, supports the proposed trans mechanism of NO reduction by NORs.

Theoretical 57Fe Mössbauer spectroscopy: isomer shifts of [Fe]-hydrogenase intermediates

A computational protocol for ⁵⁷Fe isomer shifts, based on the relativistic eXact 2-Component Hamiltonian (X2C), is applied to discriminate between proposed intermediates of [Fe]-hydrogenase. Detailed analysis reveals that the difference in isomer shifts between two intermediates is due to an overlap effect.

Nitrite activation to nitric oxide via one-fold protonation of iron(II)-O,O-nitrito complex: relevance to the nitrite reductase activity of deoxyhemoglobin and deoxyhemerythrin

The reversible transformations [(Bim)3Fe(κ(2)-O2N)][BF4] (3) <-> [(Bim)3Fe(NO)(κ(1)-ONO)][BF4]2 (4) were demonstrated and characterized. Transformation of O,O-nitrito-containing complex 3 into [(Bim)3Fe(μ-O)(μ-OAc)Fe(Bim)3](3+) (5) along with the release of NO and H2O triggered by 1 equiv of AcOH implicates that nitrite-to-nitric oxide conversion occurs, in contrast to two protons needed to trigger nitrite reduction producing NO observed in the protonation of [Fe(II)-nitro] complexes.

Preparation of non-heme {FeNO}7 models of cysteine dioxygenase: sulfur versus nitrogen ligation and photorelease of nitric oxide

We present the synthesis and spectroscopic characterization of [Fe(NO)(N3PyS)]BF4 (3), the first structural and electronic model of NO-bound cysteine dioxygenase. The nearly isostructural all-N-donor analogue [Fe(NO)(N4Py)](BF4)2 (4) was also prepared, and comparisons of 3 and 4 provide insight regarding the influence of S vs N ligation in {FeNO}(7) species. One key difference occurs upon photoirradiation, which causes the fully reversible release of NO from 3, but not from 4.

Spectroscopic and computational characterization of the NO adduct of substrate-bound Fe(II) cysteine dioxygenase: insights into the mechanism of O2 activation

Cysteine dioxygenase (CDO) is a mononuclear nonheme iron(II)-dependent enzyme critical for maintaining appropriate cysteine (Cys) and taurine levels in eukaryotic systems. Because CDO possesses both an unusual 3-His facial ligation sphere to the iron center and a rare Cys-Tyr cross-link near the active site, the mechanism by which it converts Cys and molecular oxygen to cysteine sulfinic acid is of broad interest. However, as of yet, direct experimental support for any of the proposed mechanisms is still lacking. In this study, we have used NO as a substrate analogue for O2 to prepare a species that mimics the geometric and electronic structures of an early reaction intermediate. The resultant unusual S = (1)/2 {FeNO}(7) species was characterized by magnetic circular dichroism, electron paramagnetic resonance, and electronic absorption spectroscopies as well as computational methods including density functional theory and semiempirical calculations. The NO adducts of Cys- and selenocysteine (Sec)-bound Fe(II)CDO exhibit virtually identical electronic properties; yet, CDO is unable to oxidize Sec. To explore the differences in reactivity between Cys- and Sec-bound CDO, the geometries and energies of viable O2-bound intermediates were evaluated computationally, and it was found that a low-energy quintet-spin intermediate on the Cys reaction pathway adopts a different geometry for the Sec-bound adduct. The absence of a low-energy O2 adduct for Sec-bound CDO is consistent with our experimental data and may explain why Sec is not oxidized by CDO.

A step beyond the Feltham-Enemark notation: spectroscopic and correlated ab initio computational support for an antiferromagnetically coupled M(II)-(NO)- description of Tp*M(NO) (M = Co, Ni)

Multiple spectroscopic and computational methods were used to characterize the ground-state electronic structure of the novel {CoNO}(9) species Tp*Co(NO) (Tp* = hydro-tris(3,5-Me(2)-pyrazolyl)borate). The metric parameters about the metal center and the pre-edge region of the Co K-edge X-ray absorption spectrum were reproduced by density functional theory (DFT), providing a qualitative description of the Co-NO bonding interaction as a Co(II) (S(Co) = 3/2) metal center, antiferromagnetically coupled to a triplet NO(-) anion (S(NO) = 1), an interpretation of the electronic structure that was validated by ab initio multireference methods (CASSCF/MRCI). Electron paramagnetic resonance (EPR) spectroscopy revealed significant g-anisotropy in the S = ½ ground state, but the linear-response DFT performed poorly at calculating the g-values. Instead, CASSCF/MRCI computational studies in conjunction with quasi-degenerate perturbation theory with respect to spin-orbit coupling were required for obtaining accurate modeling of the molecular g-tensor. The computational portion of this work was extended to the diamagnetic Ni analogue of the Co complex, Tp*Ni(NO), which was found to consist of a Ni(II) (S(Ni) = 1) metal center antiferromagnetically coupled to an S(NO) = 1 NO(-). The similarity between the Co and Ni complexes contrasts with the previously studied Cu analogues, for which a Cu(I) bound to NO(0) formulation has been described. This discrepancy will be discussed along with a comparison of the DFT and ab initio computational methods for their ability to predict various spectroscopic and molecular features.

Influence of thiolate ligands on reductive N-O bond activation. Probing the O2(-) binding site of a biomimetic superoxide reductase analogue and examining the proton-dependent reduction of nitrite

Nitric oxide (NO) is frequently used to probe the substrate-binding site of "spectroscopically silent" non-heme Fe(2+) sites of metalloenzymes, such as superoxide reductase (SOR). Herein we use NO to probe the superoxide binding site of our thiolate-ligated biomimetic SOR model [Fe(II)(S(Me(2))N(4)(tren))](+) (1). Like NO-bound trans-cysteinate-ligated SOR (SOR-NO), the rhombic S = 3/2 EPR signal of NO-bound cis-thiolate-ligated [Fe(S(Me(2))N(4)(tren)(NO)](+) (2; g = 4.44, 3.54, 1.97), the isotopically sensitive ν(NO)(ν((15)NO)) stretching frequency (1685(1640) cm(-1)), and the 0.05 Å decrease in Fe-S bond length are shown to be consistent with the oxidative addition of NO to Fe(II) to afford an Fe(III)-NO(-) {FeNO}(7) species containing high-spin (S = 5/2) Fe(III) antiferromagnetically coupled to NO(-) (S = 1). The cis versus trans positioning of the thiolate does not appear to influence these properties. Although it has yet to be crystallographically characterized, SOR-NO is presumed to possess a bent Fe-NO similar to that of 2 (Fe-N-O = 151.7(4)°). The N-O bond is shown to be more activated in 2 relative to N- and O-ligated {FeNO}(7) complexes, and this is attributed to the electron-donating properties of the thiolate ligand. Hydrogen-bonding to the cysteinate sulfur attenuates N-O bond activation in SOR, as shown by its higher ν(NO) frequency (1721 cm(-1)). In contrast, the ν(O-O) frequency of the SOR peroxo intermediate and its analogues is not affected by H-bonds to the cysteinate sulfur or other factors influencing the Fe-SR bond strength; these only influence the ν(Fe-O) frequency. Reactions between 1 and NO(2)(-) are shown to result in the proton-dependent heterolytic cleavage of an N-O bond. The mechanism of this reaction is proposed to involve both Fe(II)-NO(2)(-) and {FeNO}(6) intermediates similar to those implicated in the mechanism of NiR-promoted NO(2)(-) reduction.

S K-edge X-ray absorption spectroscopy and density functional theory studies of high and low spin {FeNO}7 thiolate complexes: exchange stabilization of electron delocalization in {FeNO}7 and {FeO2}8

S K-edge X-ray absorption spectroscopy (XAS) is a direct experimental probe of metal ion electronic structure as the pre-edge energy reflects its oxidation state, and the energy splitting pattern of the pre-edge transitions reflects its spin state. The combination of sulfur K-edge XAS and density functional theory (DFT) calculations indicates that the electronic structures of {FeNO}(7) (S = 3/2) (S(Me2)N4(tren)Fe(NO), complex I) and {FeNO}(7) (S = 1/2) ((bme-daco)Fe(NO), complex II) are Fe(III)(S = 5/2)-NO(-)(S = 1) and Fe(III)(S = 3/2)-NO(-)(S = 1), respectively. When an axial ligand is computationally added to complex II, the electronic structure becomes Fe(II)(S = 0)-NO•(S = 1/2). These studies demonstrate how the ligand field of the Fe center defines its spin state and thus changes the electron exchange, an important factor in determining the electron distribution over {FeNO}(7) and {FeO2}(8) sites.

Experimental and theoretical Mössbauer study of an extended family of [Fe8(μ4-O)4(μ-4-R-px)12X4] clusters

Six [Fe(8)(μ(4)-O)(4)(μ-4-R-pyrazolato)(12)X(4)] complexes containing an identical Fe(8)(μ(4)-O)(4) core have been structurally characterized and studied by Mössbauer spectroscopy. In each case, an inner μ(4)-O bridged Fe(III) cubane core is surrounded by four trigonal bipyramidal iron centers, the two distinct sites occurring in a 1:1 ratio. The Mössbauer spectrum of each of the clusters consists of two quadrupole doublets, which, with one exception (X = NCS, R = H), overlap to give three absorption lines. The systematic variation of X and R causes significant changes in the Mössbauer spectra. A comparison with values for the same clusters computed using density functional theory allows us to establish an unequivocal assignment of these peaks in terms of a nested model for the overlapping doublets. The changes in Mössbauer parameters (both experimental and computed) for the 1-electron reduced species [Fe(8)(μ(4)-O)(4)(μ-4-Cl-pyrazolato)(12)Cl(4)](-) are consistent with a redox event that is localized within the cubane core.

Spin coupling in Roussin's red and black salts

Although DFT calculations have provided a first-order electronic-structural description for Roussin's red and black salts, a detailed study of spin coupling in these species has yet to be reported. Such an analysis is presented here for the first time, based on broken-symmetry density functional theory (DFT, chiefly OLYP/STO-TZP) calculations. Both the Noodleman and Yamaguchi formulas were used to evaluate the Heisenberg coupling constants (J). Three nitrosylated binuclear clusters were studied: [Fe(2)(NO)(2)(Et-HPTB)(O(2)CPh)](2+) (1; Et-HPTB=N,N,N',N'-tetrakis-(N-ethyl-2-benzimidazolylmethyl)-2-hydroxy-1,3-diaminopropane), [Fe(NO)(2){Fe(NO)(NS(3))}-S,S'] (2), and Roussin's red salt anion [Fe(2)(NO)(4)(μ-S)(2)](2-) (3). Although the Heisenberg J for 1 is small (≈10(2) cm(-1)), 2 and 3 exhibit J values that are at least an order of magnitude higher (≈10(3) cm(-1)), where the J values refer to the following Heisenberg spin Hamiltonian: ℋ=JS(A)⋅S(B). For Roussin's black salt anion, [Fe(4)(NO)(7)(μ(3)-S)(3)](-) (4), the Heisenberg spin Hamiltonian describing spin coupling between the {FeNO}(7) unit (S(A)=3/2) and the three {Fe(NO)(2)}(9) units (S(B)=S(C)=S(D)=1/2) in [Fe(4)(NO)(7)(μ(3)-S)(3)](-) was assumed to have the form: ℋ=J(12)(S(A)⋅S(B)+S(A)⋅S(C)+S(A)⋅S(D))+J(22)(S(B)⋅S(C)+S(B)⋅S(D)+S(C)⋅S(D)), in which J(12) corresponds to the interaction between the apical iron and a basal iron, and J(22) refers to that between any two basal iron centers. Although the basal-basal coupling constant J(22) was found to be small (≈10(2) cm(-1)), the apical-basal coupling constant J(12) is some forty times higher (≈4000 cm(-1)). Thus, the nitrosylated iron-sulfur clusters feature some exceptionally high J values relative to the non-nitrosylated {2Fe2S} and {4Fe4S} clusters. An analysis of spin-dependent bonding energies shed light on this curious feature. In essence, the energy difference between the high-spin (i.e., ferromagnetically coupled iron sites) and low-spin (i.e., maximum spin coupling) states of Roussin's salts are indeed rather similar to those of analogous non-nitrosylated iron-sulfur clusters. However, the individual Fe(NO)(x) (x=1, 2) site spins are lower in the nitrosylated systems, resulting in a smaller denominator in both the Noodleman and Yamaguchi formulas for J, which in turn translates into the very high J values.

The prediction of Fe Mössbauer parameters by the density functional theory: a benchmark study

We report the performance of eight density functionals (B3LYP, BPW91, OLYP, O3LYP, M06, M06-2X, PBE, and SVWN5) in two Gaussian basis sets (Wachters and Partridge-1 on iron atoms; cc-pVDZ on the rest of atoms) for the prediction of the isomer shift (IS) and the quadrupole splitting (QS) parameters of Mössbauer spectroscopy. Two sources of geometry (density functional theory-optimized and X-ray) are used. Our data set consists of 31 iron-containing compounds (35 signals), the Mössbauer spectra of which were determined at liquid helium temperature and where the X-ray geometries are known. Our results indicate that the larger and uncontracted Partridge-1 basis set produces slightly more accurate linear correlations of electronic density used for the prediction of IS and noticeably more accurate results for the QS parameter. We confirm and discuss the earlier observation of Noodleman and co-workers that different oxidation states of iron produce different IS calibration lines. The B3LYP and O3LYP functionals have the lowest errors for either IS or QS. BPW91, OLYP, PBE, and M06 have a mixed success whereas SVWN5 and M06-2X demonstrate the worst performance. Finally, our calibrations and conclusions regarding the best functional to compute the Mössbauer characteristics are applied to candidate structures for the peroxo and Q intermediates of the enzyme methane monooxygenase hydroxylase (MMOH), and compared to experimental data in the literature.