Search

When oxidants favour cleaving a strong C–H bond at the expense of weaker ones, which are otherwise inherently preferred due to their favourable reaction energy, reactivity factors such as the polarity match effect are often invoked.

Bifurcating reactions yield two different products emerging from one single transition state and are therefore archetypal examples of reactions that cannot be described within the framework of the traditional Eyring's transition state theory (TST). With the growing number and importance of these reactions in organic and biosynthetic chemistry, there is also an increasing demand for a theoretical tool that would allow for the accurate quantification of reaction outcome at low cost. Here, we introduce such an approach that fulfils these criteria, by evaluating bifurcation selectivity through the energy distribution within the reactive mode of the key transition state. The presented method yields an excellent agreement with experimentally reported product ratios and predicts the correct selectivity for 89% of nearly 50 various cases, covering pericyclic reactions, rearrangements, fragmentations and metal-catalyzed processes as well as a series of trifurcating reactions. With 71% of product ratios determined within the error of less than 20%, we also found that the methodology outperforms three other tested protocols introduced recently in the literature. Given its predictive power, the procedure makes reaction design feasible even in the presence of complex non-TST chemical steps., Reactive Mode Composition Factor (RMCF) analysis is a powerful tool to forecast the product distribution of bifurcating reactions through analysis of the kinetic energy distribution within the first transition state traversed by the reacting system.

Complexes with two or more magnetically coupled metal ions have attracted considerable attention as catalysts of many vital processes, single-molecule magnets, or spin-crossover compounds. Elucidation of their electronic structures is essential for understanding their catalytic and magnetic properties. Here, we provide an unprecedented insight into exchange-coupling mechanisms between the magnetic centers in six prototypical bis-μ-oxo bimetallic M2O2 complexes, including two biologically relevant models of non-heme iron enzymes. Employing multiconfigurational/multireference methods and related orbital entanglement analysis, we revealed the essential and counterintuitive role of predominantly unoccupied valence metal d orbitals in their strong antiferromagnetic coupling. We found that the participation of these orbitals is twofold. First, they enhance the superexchange between the singly occupied d orbitals. Second, they become substantially occupied and thus directly magnetically active, which we perceive as a new mechanism of the exchange interaction between the magnetic transition metal centers.

By employing the computational protocol for calculation of reduction potentials of the Fe4 S4 -containing species validated using a representative series of well-defined synthetic complexes, we focused on redox properties of two prototypical radical SAM enzymes to reveal how they transform SAM into the reactive 5'-deoxyadenosyl radical, and how they tune this radical for its proper biological function. We found the reduction potential of SAM is indeed elevated by 0.3-0.4 V upon coordination to Fe4 S4 , which was previously speculated in the literature. This makes a generation of 5'-deoxyadenosyl radical from SAM less endergonic (by ca. 7-9 kcal mol-1 ) and hence more feasible in both enzymes as compared to the identical process in water. Furthermore, our calculations indicate that the enzyme-bound 5'-deoxyadenosyl radical has a significantly lower reduction potential than in referential aqueous solution, which may help the enzymes to suppress potential side redox reactions and simultaneously elevate its proton-philic character, which may, in turn, promote the radical hydrogen-atom abstraction ability.

The α-ketoglutarate (αKG)-dependent oxygenases catalyze a diverse range of chemical reactions using a common high-spin FeIV═O intermediate that, in most reactions, abstract a hydrogen atom from the substrate. Previously, the FeIV═O intermediate in the αKG-dependent halogenase SyrB2 was characterized by nuclear resonance vibrational spectroscopy (NRVS) and density functional theory (DFT) calculations, which demonstrated that it has a trigonal-pyramidal geometry with the scissile C-H bond of the substrate calculated to be perpendicular to the Fe-O bond. Here, we have used NRVS and DFT calculations to show that the FeIV═O complex in taurine dioxygenase (TauD), the αKG-dependent hydroxylase in which this intermediate was first characterized, also has a trigonal bipyramidal geometry but with an aspartate residue replacing the equatorial halide of the SyrB2 intermediate. Computational analysis of hydrogen atom abstraction by square pyramidal, trigonal bipyramidal, and six-coordinate FeIV═O complexes in two different substrate orientations (one more along [σ channel] and another more perpendicular [π channel] to the Fe-O bond) reveals similar activation barriers. Thus, both substrate approaches to all three geometries are competent in hydrogen atom abstraction. The equivalence in reactivity between the two substrate orientations arises from compensation of the promotion energy (electronic excitation within the d manifold) required to access the π channel by the significantly larger oxyl character present in the pπ orbital oriented toward the substrate, which leads to an earlier transition state along the C-H coordinate.

Iron porphyrin carbenes (IPCs) are important reaction intermediates in engineered carbene transferase enzymes and homogeneous catalysis. However, discrepancies between theory and experiment complicate the understanding of IPC electronic structure. In the literature, this has been framed as whether the ground state is an open- vs closed-shell singlet (OSS vs CSS). Here we investigate the structurally dependent ground and excited spin-state energetics of a free carbene and its IPC analogs with variable trans axial ligands. In particular, for IPCs, multireference ab initio wave function methods are more consistent with experiment and predict a mixed singlet ground state that is dominated by the CSS (Fe(II) ← {:C(X)Y}0) configuration (i.e., electrophilic carbene) but that also has a small, non-negligible contribution from an Fe(III)–{C(X)Y}−• configuration (hole in d(xz), i.e., radical carbene). In the multireference approach, the “OSS-like” excited states are metal-to-ligand charge transfer (MLCT) in nature and are energetically well above the CSS-dominated ground state. The first, lowest energy of these “OSS-like” excited states is predicted to be heavily weighted toward the Fe(III)–{C(X)Y}−• (hole in d(yz)) configuration. As expected from exchange considerations, this state falls energetically above a triplet of the same configuration. Furthermore, potential energy surfaces (PESs) along the IPC Fe–C(carbene) bond elongation exhibit increasingly strong mixings between CSS/OSS characters, with the Fe(III)–{C(X)Y}−• configuration (hole in d(xz)) growing in weight in the ground state during bond elongation. The relative degree of electrophilic/radical carbene character along this structurally relevant PES can potentially play a role in reactivity and selectivity patterns in catalysis. Future studies on IPC reaction coordinates should evaluate contributions from ground and excited state multireference character.

A full understanding of the catalytic action of non-heme iron (NHFe) and non-heme diiron (NHFe2) enzymes is still beyond the grasp of contemporary computational and experimental techniques. Many of these enzymes exhibit fascinating chemo-, regio-, and stereoselectivity, in spite of employing highly reactive intermediates which are necessary for activations of most stable chemical bonds. Herein, we study in detail one intriguing representative of the NHFe2 family of enzymes: soluble Δ9 desaturase (Δ9D), which desaturates rather than performing the thermodynamically favorable hydroxylation of substrate. Its catalytic mechanism has been explored in great detail by using QM(DFT)/MM and multireference wave function methods. Starting from the spectroscopically observed 1,2-μ-peroxo diferric P intermediate, the proton-electron uptake by P is the favored mechanism for catalytic activation, since it allows a significant reduction of the barrier of the initial (and rate-determining) H-atom abstraction from the stearoyl substrate as compared to the "proton-only activated" pathway. Also, we ruled out that a Q-like intermediate (high-valent diamond-core bis-μ-oxo-[FeIV]2 unit) is involved in the reaction mechanism. Our mechanistic picture is consistent with the experimental data available for Δ9D and satisfies fairly stringent conditions required by Nature: the chemo-, stereo-, and regioselectivity of the desaturation of stearic acid. Finally, the mechanisms evaluated are placed into a broader context of NHFe2 chemistry, provided by an amino acid sequence analysis through the families of the NHFe2 enzymes. Our study thus represents an important contribution toward understanding the catalytic action of the NHFe2 enzymes and may inspire further work in NHFe(2) biomimetic chemistry.

Hydrogen atom abstraction (HAA) is central to life and its importance in synthetic chemistry continues to grow. Enzymes rely on HAA to trigger life-sustaining reaction cascades, and greener synthetic routes are attainable by in situ capture of the carbon-centered radicals generated by HAA. Despite the potential of HAA for the diversification of molecular complexity and the late-stage functionalization of bioactive compounds, readily applicable and reliable models translating experimentally or computationally accessible thermodynamic quantities into relative free energy barriers are missing. In this work, we discovered a complete thermodynamic basis for the description of HAA reactivity, which consists of three components. Besides, the traditional linear free energy relationship and the recently introduced factor of asynchronicity (Srnec et al, PNAS 2018, 115, E10287-E10294), we present the third thermodynamic component of H-atom abstraction reactions – factor of frustration that arises from the dissimilarity of the species competing over a hydrogen atom in their overall ability to acquire electron and proton. Incorporating these non-classical descriptors into a Marcus-type model, the approach herein presented allows nearly quantitative prediction of relative barriers in six sets of metal-oxo-mediated HAA reactions, outperforming existing methods even in a stringent test with >200 computational HAA reactions.

The rate of elementary reactions usually rises with increasing temperature. In rare cases, however, a slowdown is observed instead. One example of this is hydrogen-atom abstraction from the iron(II...

We report and analyze chemoselectivity in the gas phase reactions of cycloalkenes (cyclohexene, cycloheptene, cis-cyclooctene, 1,4-cyclohexadiene) with a non-heme iron(IV)-oxo complex [(PyTACN)Fe(O)(Cl)]+, which models the active species in iron-dependent halogenases. Unlike in the halogenases, we did not observe any chlorination of the substrate. However, we observed two other reaction pathways: allylic hydrogen atom transfer (HAT) and alkene epoxidation. The HAT is clearly preferred in the case of 1,4-cyclohexadiene, both pathways have comparable reaction rates in reaction with cyclohexene, and epoxidation is strongly favored in reactions with cycloheptene and cis-cyclooctene. This preference for epoxidation differs from the reactivity of iron(IV)-oxo complexes in the condensed phase, where HAT usually prevails. To understand the observed selectivity, we analyze effects of the substrate, spin state, and solvation. Our DFT and CASPT2 calculations suggest that all the reactions occur on the quintet potential energy surface. The DFT-calculated energies of the transition states for the epoxidation and hydroxylation pathways explain the observed chemoselectivity. The SMD implicit solvation model predicts the relative increase of the epoxidation barriers with solvent polarity, which explains the clear preference of HAT in the condensed phase. Electronic supplementary material The online version of this article (10.1007/s13361-019-02251-1) contains supplementary material, which is available to authorized users.

Recently, the analysis of single-orbital entropy and mutual information has been introduced as a tool for the investigation of contributions to the exchange (J) coupling between open-shell metal ions [Stein et al. J. Phys. Chem. Lett. 2019, 10, 6762-6770]. Here, we show that this analysis may lead to an incorrect interpretation of the J-coupling mechanism. Instead, we propose an orbital-entanglement analysis that is based on the two-electron density and that provides a coherent picture of the contributing exchange pathways, which seems fully consistent with the available J values. For this purpose, we used a prototypical bis-μ-oxo binuclear manganese complex ([Mn2O2(NH3)8]4+) and demonstrated that its antiferromagnetism (J < 0), calculated by using the active space composed of all valence pO and dMn orbitals, correlates well with the largest elements in the differential low-spin vs high-spin entanglement map. These elements correspond to interactions between the pairs of dMn orbitals mediated by the oxo-bridging out-of-plane p orbitals, representing the π superexchange pathway. We also show that the reduction of active space to manifold of the singly occupied magnetic orbitals does not lead to discrepancy between the calculated J values and entanglement maps. This contrasts to analysis of mutual information, which suggests the "direct" dMn-dMn interactions to play a dominant role for the J coupling, irrespective of the size of active space as well as of the antiferromagnetism expected. The failure is attributed to the large contribution of spin entanglement contained in the mutual information of the low-spin state, which may be regarded as the origin of the different complexity of its wave function and electron density.

The α-ketoglutarate (αKG)-dependent oxygenases catalyze a diverse range of chemical reactions using a common high-spin Fe(IV)=O intermediate and, in most reactions, abstracts a hydrogen atom from the substrate. Previously, the Fe(IV)=O intermediate in the αKG-dependent halogenase SyrB2 was characterized by nuclear resonance vibrational spectroscopy (NRVS) and density functional theory (DFT) calculations, which demonstrated that it has a trigonal-pyramidal geometry with the scissile C-H bond of the substrate calculated to be perpendicular to the Fe-O bond. Here, we have used NRVS and DFT calculations to show that the Fe(IV)=O complex in taurine dioxygenase (TauD), the αKG-dependent hydroxylase in which this intermediate was first characterized, also has a trigonal bipyramidal geometry, but with an aspartate residue replacing the equatorial halide of the SyrB2 intermediate. Computational analysis of hydrogen atom abstraction by square pyramidal, trigonal bipyramidal, and six-coordinate Fe(IV)=O complexes in two different substrate orientations (one more along [σ channel] and another more perpendicular [π channel] to the Fe-O bond) reveals similar activation barriers. Thus, both substrate approaches to all three geometries are competent in hydrogen atom abstraction. The equivalence in reactivity between the two substrate orientations arises from compensation of the promotion energy (electronic excitation within the d manifold) required to access the π channel by the significantly larger oxyl character present in the pπ orbital oriented towards the substrate, which leads to an earlier transition state along the C-H coordinate.

By employing the computational protocol for calculation of reduction potentials of the Fe

The selective functionalization of C-H bonds is one of the Grails of synthetic chemistry. In this work, we demonstrate that the selectivity toward fast hydroxylation or radical diffusion (known as the OH-rebound and dissociation mechanisms) following H-atom abstraction (HAA) from a substrate C-H bond by high-valent iron-oxo oxidants is already encoded in the HAA step when the post-HAA barriers are much lower than the preceding one. By applying the reactive mode composition factor (RMCF) analysis, which quantifies the kinetic energy distribution (KED) at the reactive mode (RM) of transition states, we show that reactions following the OH-rebound coordinate concentrate the RM kinetic energy on the motion of the reacting oxygen atom and the nascent substrate radical, whereas reactions following the dissociation channel localize most of their kinetic energy in H-atom motion. These motion signatures serve to predict the post-HAA selectivity, and since KED is affected by the free energy of reaction and asynchronicity (factor η) of HAA, we show that bimolecular HAA reactions in solution that are electron transfer-driven and highly exergonic have the lowest fraction of KED on the transferred H-atom and the highest chance to follow rebound hydroxylation. Finally, the RMCF analysis predicts that the H/D primary kinetic isotope effect can serve as a probe for these mechanisms, as confirmed in virtually all reported examples in the literature.

Despite the growing number of reported FeIVO complexes, an unambiguous experimental characterization of their redox properties, such as one-electron reduction potentials, remains a challenging task. To this aim, we describe an efficient and straightforward theoretical protocol for accurate calculations of redox potentials and calibrate the protocol on a set of diverse 37 mononuclear nonheme iron (NHFe) redox couples. It is shown that the methodology, further applied to a set of 10 FeIVO species, not only serves for near-quantitative predictions of reduction potentials, but also is an elegant tool for interpretation of the experimental electrochemical data. The general need for such a computational methodology is illustrated on the prototypical example of the (N4Py)FeIVO complex, whose electrochemistry has been studied for many years and still raises many questions.

We report the generation of iron(V) nitride complexes, which are targets of biomimetic chemistry. Temperature-dependent ion spectroscopy shows that this reaction is governed by the spin state population of their iron(III) azide precursors and can be tuned by temperature. The complex [(MePy2TACN)Fe(N3)]2+ exists as a mixture of sextet and doublet spin states at 300 K, whereas only the doublet state is populated at 3 K. Photofragmentation of the sextet state complex leads to the reduction of the iron centre. The doublet state complex photodissociates to the desired iron(V) nitride complex. To generalize these findings, we show results for complexes with cyclam-based ligands.

Binuclear non-heme iron enzymes activate O2 for diverse chemistries that include oxygenation of organic substrates and hydrogen atom abstraction. This process often involves the formation of peroxo-bridged biferric intermediates, only some of which can perform electrophilic reactions. To elucidate the geometric and electronic structural requirements to activate peroxo reactivity, the active peroxo intermediate in 4-aminobenzoate N-oxygenase (AurF) has been characterized spectroscopically and computationally. A magnetic circular dichroism study of reduced AurF shows that its electronic and geometric structures are poised to react rapidly with O2. Nuclear resonance vibrational spectroscopic definition of the peroxo intermediate formed in this reaction shows that the active intermediate has a protonated peroxo bridge. Density functional theory computations on the structure established here show that the protonation activates peroxide for electrophilic/single-electron-transfer reactivity. This activation of peroxide by protonation is likely also relevant to the reactive peroxo intermediates in other binuclear non-heme iron enzymes.

Iron porphyrin carbenes (IPCs) are important reaction intermediates in engineered carbene transferase enzymes and homogeneous catalysis. However, discrepancies between theory and experiment complicate the understanding of IPC electronic structure (i.e., open- vs. closed-shell singlet (OSS vs. CSS)). Here we investigate the structurally dependent ground and excited spin state energetics of a free carbene and its IPC analogs. Only multireference ab initio wave function methods are consistent with experiment and predict a CSS ground state (Fe(II)←{:C(X)Y}0), contrary to density functional theory (DFT). The OSS is a high-lying metal-to-ligand charge transfer (MLCT) excited state that is sensitive to the nature of the axial ligand. Furthermore, potential energy surfaces (PESs) along the Fe–C bond elongation coordinate exhibit strong mixings between CSS/OSS characters, which can be an important feature for describing reaction mechanisms. Future studies on IPC reaction coordinates should evaluate contributions from ground and excited state multireference character.

A simple method for the evaluation of the kinetic energy distribution within the reactive mode of a transition state (TS), denoted as the Reactive Mode Composition Factor (RMCF), is presented. It allows one to directly map the barrier properties onto the atomic-motion components of the reaction coordinate at the TS, which has potential to shed light onto some mechanistic features of a chemical process. To demonstrate the applicability of RMCF to reactivity, we link the kinetic energy distribution within a reactive mode with the asynchronicity (η) in C-H bond activation, as they both evolve in a series of coupled proton-electron transfer (CPET) reactions between FeIVO oxidants and 1,4-cyclohexadiene. RMCF shows how the earliness or lateness of a process manifests as a redistribution of kinetic energy in the reactive mode as a function of the free energy of reaction (ΔG0) and η. Finally, the title analysis can be applied to predict H-atom tunneling contributions and kinetic isotope effects in a set of reactions, yielding a transparent rationalization based on the kinetic energy distributions in the reactive mode.

Aquest mateix article està publicat a l'edició alemanya d''Angewandte Chemie' (ISSN 0044-8249, EISSN 1521-3757), 2019, vol.131, núm. 28, p. 9721-9726. DOI https://doi.org/10.1002/ange.201904546 Terminal oxo complexes of late transition metals are frequently proposed reactive intermediates. However, they are scarcely known beyond Group 8. Using mass spectrometry, we prepared and characterized two such complexes: [(N4Py)CoIII(O)]+ (1) and [(N4Py)CoIV(O)]2+ (2). Infrared photodissociation spectroscopy revealed that the Co−O bond in 1 is rather strong, in accordance with its lack of chemical reactivity. On the contrary, 2 has a very weak Co−O bond characterized by a stretching frequency of ≤659 cm−1. Accordingly, 2 can abstract hydrogen atoms from non-activated secondary alkanes. Previously, this reactivity has only been observed in the gas phase for small, coordinatively unsaturated metal complexes. Multireference ab-initio calculations suggest that 2, formally a cobalt(IV)-oxo complex, is best described as cobalt(III)-oxyl. Our results provide important data on changes to metal-oxo bonding behind the oxo wall and show that cobalt-oxo complexes are promising targets for developing highly active C−H oxidation catalysts The project was funded by the European Research Council (ERC CoG No. 682275), the Czech Ministry of Education, Youth and Sports (LTAUSA17026), the COST action ECOSTBio, MINECO of Spain (CTQ2015‐70795‐P), the Catalan DIUE of the Generalitat de Catalunya (2017SGR01378, a BFI PhD grant to M.R., and an ICREA‐Academia award), and the Grant Agency of the Czech Republic (Grant No. 18‐13093S)

Low temperature magnetic circular dichroism (LT MCD) spectroscopy in combination with quantum-chemical calculations are used to define the electronic structure associated with the geometric structure of the Fe(IV)═O intermediate in SyrB2 that was previously determined by nuclear resonance vibrational spectroscopy. These studies elucidate key frontier molecular orbitals (FMOs) and their contribution to H atom abstraction reactivity. The VT MCD spectra of the enzymatic S = 2 Fe(IV)═O intermediate with Br(-) ligation contain information-rich features that largely parallel the corresponding spectra of the S = 2 model complex (TMG3tren)Fe(IV)═O (Srnec, M.; Wong, S. D.; England, J; Que, L; Solomon, E. I. Proc. Natl. Acad. Sci. USA 2012, 109, 14326-14331). However, quantitative differences are observed that correlate with π-anisotropy and oxo donor strength that perturb FMOs and affect reactivity. Due to π-anisotropy, the Fe(IV)═O active site exhibits enhanced reactivity in the direction of the substrate cavity that proceeds through a π-channel that is controlled by perpendicular orientation of the substrate C-H bond relative to the halide-Fe(IV)═O plane. Also, the increased intrinsic reactivity of the SyrB2 intermediate relative to the ferryl model complex is correlated to a higher oxyl character of the Fe(IV)═O at the transition states resulting from the weaker ligand field of the halogenase.

Hydrogen atom abstraction (HAA) reactions are cornerstones of chemistry. Various (metallo)enzymes performing the HAA catalysis evolved in nature and inspired the rational development of multiple synthetic catalysts. Still, the factors determining their catalytic efficiency are not fully understood. Herein, we define the simple thermodynamic factor η by employing two thermodynamic cycles: one for an oxidant (catalyst), along with its reduced, protonated, and hydrogenated form; and one for the substrate, along with its oxidized, deprotonated, and dehydrogenated form. It is demonstrated that η reflects the propensity of the substrate and catalyst for (a)synchronicity in concerted H + /e − transfers. As such, it significantly contributes to the activation energies of the HAA reactions, in addition to a classical thermodynamic (Bell–Evans–Polanyi) effect. In an attempt to understand the physicochemical interpretation of η, we discovered an elegant link between η and reorganization energy λ from Marcus theory. We discovered computationally that for a homologous set of HAA reactions, λ reaches its maximum for the lowest |η|, which then corresponds to the most synchronous HAA mechanism. This immediately implies that among HAA processes with the same reaction free energy, Δ G 0 , the highest barrier (≡Δ G ≠ ) is expected for the most synchronous proton-coupled electron (i.e., hydrogen) transfer. As proof of concept, redox and acidobasic properties of nonheme Fe IV O complexes are correlated with activation free energies for HAA from C−H and O−H bonds. We believe that the reported findings may represent a powerful concept in designing new HAA catalysts.

The Rieske dioxygenases are a major subclass of mononuclear nonheme iron enzymes that play an important role in bioremediation. Recently, a high-spin FeIII–(hydro)-peroxy intermediate (BZDOp) has been trapped in the peroxide shunt reaction of benzoate 1,2-dioxygenase. Defining the structure of this intermediate is essential to understanding the reactivity of these enzymes. Nuclear resonance vibrational spectroscopy (NRVS) is a recently developed synchrotron technique that is ideal for obtaining vibrational, and thus structural, information on Fe sites, as it gives complete information on all vibrational normal modes containing Fe displacement. In this study, we present NRVS data on BZDOp and assign its structure using these data coupled to experimentally calibrated density functional theory calculations. From this NRVS structure, we define the mechanism for the peroxide shunt reaction. The relevance of the peroxide shunt to the native FeII/O2 reaction is evaluated. For the native FeII/O2 reaction, an FeIII–superoxo intermediate is found to react directly with substrate. This process, while uphill thermodynamically, is found to be driven by the highly favorable thermodynamics of proton-coupled electron transfer with an electron provided by the Rieske [2Fe-2S] center at a later step in the reaction. These results offer important insight into the relative reactivities of FeIII–superoxo and FeIII–hydroperoxo species in nonheme Fe biochemistry.

The CuO+ core is a central motif of reactive intermediates in copper-catalysed oxidations occurring in nature. The high reactivity of CuO+ stems from a weak bonding between the atoms, which cannot be described by a simple classical model. To obtain the correct picture, we have investigated the acetonitrile-ligated CuO+ ion using neon-tagging photodissociation spectroscopy at 5 K. The spectra feature complex vibronic absorption progressions in NIR and visible regions. Employing Franck-Condon analyses, we derived low-lying triplet potential energy surfaces that were further correlated with multireference calculations. This provided insight into the ground and low-lying excited electronic states of the CuO+ unit and elucidated how these states are perturbed by the change in ligation. Thus, we show that the bare CuO+ ion has prevailingly a copper(I)-biradical oxygen character. Increasing the number of ligands coordinated to copper changes the CuO+ character towards the copper(II)-oxyl radical structure.

We report for the first time infrared spectra of three non-heme pseudo-octahedral iron(V) nitride complexes with assigned Fe-N stretching vibrations. The intensities of the Fe-N bands in two of the complexes are extremely weak. Their detection was enabled by the high resolution and sensitivity of the experiments performed at 3 K for isolated complexes in the gas phase. Multireference CASPT2 calculations revealed that the Fe-N bond in the ground doublet state is influenced by two low-lying excited doublet states. In particular, configuration interaction between the ground and the second excited state leads to avoided crossing of their potential energy surfaces along the Fe-N coordinate, which thus affects the ground-state Fe-N stretching frequency and intensity. Therefore, DFT calculated Fe-N stretching frequency strongly depends on the amount of Hartree-Fock exchange potential. As a result, by tuning the amount of Hartree-Fock exchange potential in the B3LYP functional, it was possible to obtain theoretical spectra perfectly consistent with the experimental data. The theory shows that the intensity of the Fe-N stretching vibration can almost vanish due to strong coupling with other stretching modes of the ligands.

The ability of an FeIV=O intermediate in SyrB2 to perform chlorination vs. hydroxylation was computationally evaluated for different substrates that had been studied experimentally. The π-trajectory for H-atom abstraction (FeIV=O oriented perpendicular to the C—H bond of substrate) was found to lead to the S = 2 five-coordinate HO—FeIII—Cl complex with the C• of the substrate, π-oriented relative to both the Cl− and OH− ligands. From this ferric intermediate, hydroxylation is thermodynamically favored, but chlorination is intrinsically more reactive due to the energy splitting between two key redox-active dπ* frontier molecular orbitals (FMOs). The splitting is determined by the differential ligand field effect of Cl− vs. OH− on the Fe center. This makes chlorination effectively competitive with hydroxylation. Chlorination vs. hydroxylation selectivity is then determined by the orientation of the substrate with respect to the HO—Fe—Cl plane that controls either the Cl− or OH− to rebound depending on the relative π-overlap with the substrate C radical. The differential contribution of the two FMOs to chlorination vs. hydroxylation selectivity in SyrB2 is related to a reaction mechanism that involves two asynchronous transfers: electron transfer from the substrate radical to the iron center followed by late ligand (Cl− or OH−) transfer to the substrate.

The results of density matrix renormalization group complete active space self-consistent field (DMRG-CASSCF) and second-order perturbation theory (DMRG-CASPT2) calculations are presented on various structural alternatives for the O-O and first C-H activating step of the catalytic cycle of the binuclear nonheme iron enzyme Δ(9) desaturase. This enzyme is capable of inserting a double bond into an alkyl chain by double hydrogen (H) atom abstraction using molecular O2. The reaction step studied here is presumably associated with the highest activation barrier along the full pathway; therefore, its quantitative assessment is of key importance to the understanding of the catalysis. The DMRG approach allows unprecedentedly large active spaces for the explicit correlation of electrons in the large part of the chemically important valence space, which is apparently conditio sine qua non for obtaining well-converged reaction energetics. The derived reaction mechanism involves protonation of the previously characterized 1,2-μ peroxy Fe(III)Fe(III) (P) intermediate to a 1,1-μ hydroperoxy species, which abstracts an H atom from the C10 site of the substrate. An Fe(IV)-oxo unit is generated concomitantly, supposedly capable of the second H atom abstraction from C9. In addition, several popular DFT functionals were compared to the computed DMRG-CASPT2 data. Notably, many of these show a preference for heterolytic C-H cleavage, erroneously predicting substrate hydroxylation. This study shows that, despite its limitations, DMRG-CASPT2 is a significant methodological advancement toward the accurate computational treatment of complex bioinorganic systems, such as those with the highly open-shell diiron active sites.

The non-heme ferryl active sites are of significant interest for their application in biomedical and green catalysis. These sites have been shown to have an S = 1 or S = 2 ground spin state; the latter is functional in biology. Low-temperature magnetic circular dichroism (LT MCD) spectroscopy probes the nature of the excited states in these species including ligand-field (LF) states that are otherwise difficult to study by other spectroscopies. In particular, the temperature dependences of MCD features enable their unambiguous assignment and thus determination of the low-lying excited states in two prototypical S = 1 and S = 2 NHFe(IV)[double bond, length as m-dash]O complexes. Furthermore, some MCD bands exhibit vibronic structures that allow mapping of excited-state interactions and their effects on the potential energy surfaces (PESs). For the S = 2 species, there is also an unusual spectral feature in both near-infrared absorption and MCD spectra - Fano antiresonance (dip in Abs) and Fano resonance (sharp peak in MCD) that indicates the weak spin-orbit coupling of an S = 1 state with the S = 2 LF state. These experimental data are correlated with quantum-chemical calculations that are further extended to analyze the low-lying electronic states and the evolution of their multiconfigurational characters along the Fe-O PESs. These investigations show that the lowest-energy states develop oxyl Fe(III) character at distances that are relevant to the transition state (TS) for H-atom abstraction and define the frontier molecular orbitals that participate in the reactivity of S = 1 vs. S = 2 non-heme Fe(IV)[double bond, length as m-dash]O active sites. The S = 1 species has only one available channel that requires the C-H bond of a substrate to approach perpendicular to the Fe-oxo bond (the π channel). In contrast, there are three channels (one σ and two π) available for the S = 2 non-heme Fe(IV)[double bond, length as m-dash]O system allowing C-H substrate approach both along and perpendicular to the Fe-oxo bond that have important implications for enzymatic selectivity.

Mononuclear non-haem iron (NHFe) enzymes catalyse a broad range of oxidative reactions, including halogenation, hydroxylation, ring closure, desaturation and aromatic ring cleavage reactions. They are involved in a number of biological processes, including phenylalanine metabolism, the production of neurotransmitters, the hypoxic response and the biosynthesis of secondary metabolites. The reactive intermediate in the catalytic cycles of these enzymes is a high-spin S = 2 Fe(IV)=O species, which has been trapped for a number of NHFe enzymes, including the halogenase SyrB2 (syringomycin biosynthesis enzyme 2). Computational studies aimed at understanding the reactivity of this Fe(IV)=O intermediate are limited in applicability owing to the paucity of experimental knowledge about its geometric and electronic structure. Synchrotron-based nuclear resonance vibrational spectroscopy (NRVS) is a sensitive and effective method that defines the dependence of the vibrational modes involving Fe on the nature of the Fe(IV)=O active site. Here we present NRVS structural characterization of the reactive Fe(IV)=O intermediate of a NHFe enzyme, namely the halogenase SyrB2 from the bacterium Pseudomonas syringae pv. syringae. This intermediate reacts via an initial hydrogen-atom abstraction step, performing subsequent halogenation of the native substrate or hydroxylation of non-native substrates. A correlation of the experimental NRVS data to electronic structure calculations indicates that the substrate directs the orientation of the Fe(IV)=O intermediate, presenting specific frontier molecular orbitals that can activate either selective halogenation or hydroxylation.

In this minireview, we provide an account of the current state-of-the-art developments in the area of mono- and binuclear non-heme enzymes (NHFe and NHFe2) and the smaller NHFe(2) synthetic models, mostly from a theoretical and computational perspective. The sheer complexity, and at the same time the beauty, of the NHFe(2) world represents a challenge for experimental as well as theoretical methods. We emphasize that the concerted progress on both theoretical and experimental side is a conditio sine qua non for future understanding, exploration and utilization of the NHFe(2) systems. After briefly discussing the current challenges and advances in the computational methodology, we review the recent spectroscopic and computational studies of NHFe(2) enzymatic and inorganic systems and highlight the correlations between various experimental data (spectroscopic, kinetic, thermodynamic, electrochemical) and computations. Throughout, we attempt to keep in mind the most fascinating and attractive phenomenon in the NHFe(2) chemistry, which is the fact that despite the strong oxidative power of many reactive intermediates, the NHFe(2) enzymes perform catalysis with high selectivity. We conclude with our personal viewpoint and hope that further developments in quantum chemistry and especially in the field of multireference wave function methods are needed to have a solid theoretical basis for the NHFe(2) studies, mostly by providing benchmarking and calibration of the computationally efficient and easy-to-use DFT methods.

Large-scale quantum and molecular mechanical methods (QM/MM) and QM calculations were carried out on the soluble Δ(9) desaturase (Δ(9)D) to investigate various structural models of the spectroscopically defined peroxodiferric (P) intermediate. This allowed us to formulate a consistent mechanistic picture for the initial stages of the reaction mechanism of Δ(9)D, an important diferrous nonheme iron enzyme that cleaves the C-H bonds in alkane chains resulting in the highly specific insertion of double bonds. The methods (density functional theory (DFT), time-dependent DFT (TD-DFT), QM(DFT)/MM, and TD-DFT with electrostatic embedding) were benchmarked by demonstrating that the known spectroscopic effects and structural perturbation caused by substrate binding to diferrous Δ(9)D can be qualitatively reproduced. We show that structural models whose spectroscopic (absorption, circular dichroism (CD), vibrational and Mössbauer) characteristics correlate best with experimental data for the P intermediate correspond to the μ-1,2-O(2)(2-) binding mode. Coordination of Glu196 to one of the iron centers (Fe(B)) is demonstrated to be flexible, with the monodentate binding providing better agreement with spectroscopic data, and the bidentate structure being slightly favored energetically (1-10 kJ mol(-1)). Further possible structures, containing an additional proton or water molecule are also evaluated in connection with the possible activation of the P intermediate. Specifically, we suggest that protonation of the peroxide moiety, possibly preceded by water binding in the Fe(A) coordination sphere, could be responsible for the conversion of the P intermediate in Δ(9)D into a form capable of hydrogen abstraction. Finally, results are compared with recent findings on the related ribonucleotide reductase and toluene/methane monooxygenase enzymes.

We used two theoretical methods to estimate reduction potentials and acidity constants in Mn superoxide dismutase (MnSOD), namely combined quantum mechanical and molecular mechanics (QM/MM) thermodynamic cycle perturbation (QTCP) and the QM/MM-PBSA approach. In the latter, QM/MM energies are combined with continuum solvation energies calculated by solving the Poisson-Boltzmann equation (PB) or by the generalised Born approach (GB) and non-polar solvation energies calculated from the solvent-exposed surface area. We show that using the QTCP method, we can obtain accurate and precise estimates of the proton-coupled reduction potential for MnSOD, 0.30±0.01 V, which compares favourably with experimental estimates of 0.26-0.40 V. However, the calculated potentials depend strongly on the DFT functional used: The B3LYP functional gives 0.6 V more positive potentials than the PBE functional. The QM/MM-PBSA approach leads to somewhat too high reduction potentials for the coupled reaction and the results depend on the solvation model used. For reactions involving a change in the net charge of the metal site, the corresponding results differ by up to 1.3 V or 24 pK(a) units, rendering the QM/MM-PBSA method useless to determine absolute potentials. However, it may still be useful to estimate relative shifts, although the QTCP method is expected to be more accurate.

Reduction potentials of several M(2+/3+) (M = Ru, Os) octahedral complexes, namely, [M(H2O)6](2+/3+), [MCl6](4-/3-), [M(NH3)6](2+/3+), [M(en)3](2+/3+) [M(bipy)3](2+/3+), and [M(CN)6](4-/3-), were calculated using the CASSCF/CASPT2/CASSI and MRCI methods including spin-orbit coupling (SOC) by means of first-order quasi-degenerate perturbation theory. It was shown that the effect of SOC accounts for a systematic shift of approximately -70 mV in the reduction potentials of the studied ruthenium (II/III) complexes and an approximately -300 mV shift for the osmium(II/III) complexes. SOC splits the sixfold-degenerate (2)T(2g) ground electronic state (in ideal octahedral symmetry) of the M(3+) ions into the E((5/2)g) Kramers doublet and G((3/2)g) quartet, which were calculated to split by 1354-1573 cm(-1) in the Ru(3+) complexes and 4155-5061 cm(-1) in the Os(3+) complexes. It was demonstrated that this splitting represents the main contribution to the stabilization of the M(3+) ground state with respect to the closed-shell (1)A(1g) ground state in M(2+) systems. Moreover, it was shown that the accuracy of the calculated reduction potentials depends on the calculated solvation energies of both the oxidized and reduced forms. For smaller ligands, it involves explicit inclusion of the second solvation sphere into the calculations, whereas implicit solvation models yield results of sufficient accuracy for complexes with larger ligands. In such cases (e.g., [M(bipy)3](2+/3+) and its derivatives), very good agreement between the calculated (SOC-corrected) values of the reduction potentials and the available experimental values was obtained. These results led us to the conclusion that especially for Os(2+/3+) complexes, inclusion of SOC is necessary to avoid systematic errors of approximately 300 mV in the calculated reduction potentials.

The influence of geometrical changes on the spin multiplicity of the ground states of the octahedral ruthenium(II/III) and osmium(II/III) complexes is investigated using the TD-DFT and MRCI methods. On the example of the [RuCl6]4- complex, we show that only after the optimisation of the molecular geometry in a solvent (using a polarised continuum model), which shortens the M-L bond lengths by ~0.2 Å, is the correct order of spin states obtained (i.e. a singlet is correctly predicted to be the ground state). On the contrary, in terms of the in vacuo optimised geometries of this negatively charged species, both the DFT and MRCI calculations predict a quintet ground state. This finding is further analysed by calculating the low- and high-spin potential energy curves corresponding to an elongation of the M-L distance, which makes it possible to predict the critical point at which the crossing of the two spin states occurs. Finally, it is complemented by the TD-DFT calculations of the lowest excited states in each spin multiplicity for a series of prototypical ligands. It is demonstrated that the calculated results correlate well with the known strengths of the ligand field. The two findings presented in this work are a small contribution to our understanding of the electronic structure and properties of the octahedral ruthenium(II/III) and osmium(II/III) complexes, which are relevant both in biomolecular and material sciences.

In the monomeric form, the title systems assume unexpected optical properties, and as oligomers they can serve as candidates for molecular devices. In bulk they are attractive in the area of material science. A broad set of quantum chemical methods ranging from density functional theory by the Hartree–Fock method and the Moller–Plesset perturbation theory to the coupled-cluster method, in connection with nonrelativistic Dunning's basis sets as well as with relativistic SDD basis sets were used. Electronic spectra were analyzed by means of time-dependent DFT, with the symmetry-adapted cluster configuration interaction and with the internally contracted multireference configuration interaction. The Douglas–Kroll–Hess quasirelativistic Hamiltonian served as a basis for the estimation of the role of relativistic effects. All the diatomics representing 25 elements of a matrix formed by the combination of group IIa (Be, ..., Ba) with group VIa (O, ..., Po) atoms were characterized by calculated bond length, valence vibration, dipole moment, and the Hartree–Fock frontier orbital energies. Calculated characteristics were confronted with experimental data. The geometry of the oligomers (from dimers through hexamers) was investigated systematically, and the nature of the located stationary points on the respective potential energy surfaces was established. Special attention was paid to the four lightest elements of the group IIa–VIa matrix, i.e. BeO, BeS, MgO, and MgS. As for electronic spectra, the systems of the first column (BeO, ..., BaO), the first row (BeO, ..., BePo), and the systems of the main diagonal (BeO, ..., BaPo) were calculated in the form of monomers and dimers. Analogous calculations were performed for a few group Ia–VIIa, IIIa–Va, IVa–IVa, and IIb–VIa systems. The group IIa–VIa and IIb–VIa diatomics exhibit electronic transitions in the visible region of the spectrum, and the longest wavelength bands are located in the near-infrared region. MgO represents an extreme with the first band appearing at about 3500 cm–1, which strictly speaking makes the use of the Born–Oppenheimer approximation questionable. Whenever experimental transitions are available, the agreement between the calculated and observed band positions is good. Passing from monomers to oligomers is always associated with a significant hypsochromic shift in the first transition.(© Wiley-VCH Verlag GmbH & Co. KGaA, 69451 Weinheim, Germany, 2007)

Ab initio calculations of structural stability and properties of selected IIIa-Va monomers (BN, BP, AlN, AlP), and dimers and trimers thereof are presented. It was demonstrated that multireference methods are needed for an appropriate description of the electronic structure of the diatomics under study. The energy differences between the quasi-degenerated states of monomers were determined: in the range of 220-257 cm(-1) for the X3Pi --> A1Sigma+ in BN, 277-311 cm(-1) for the X3Pi --> A3Sigma(-) in AlN, and 344-768 cm(-1) for the X3Sigma(-) --> A3Pi in the AlP molecule. Many stable isomers of dimers and trimers were optimized and the structural patterns of the species were established. Although the (BN)(i) isomers prefer linear and cyclic geometries with a metal-nonmetal sequence, the structure of the (AlN)(i) species is determined by strong N--N interactions; the geometric patterns of the (BP)(i) and (AlP)(i) systems are more involved. The electronic structure and physical and chemical properties of the most stable dimers and trimers were calculated. As to electronic spectra, absorption bands in the visible and infrared region were found. The calculated low-lying excited states, electron distribution and partially biradical character of structures along with the variety of stable isomers indicate high chemical reactivity of the studied IIIa-Va microclusters.

A robust and efficient theoretical approach for calculation of the reduction potentials of charged species in aqueous solution is presented. Within this approach, the reduction potential of a charged complex (with a charge |n| ≥ 2) is probed by means of the reduction potential of its neutralized (protonated/deprotonated) cognate, employing one or several H-atom addition/abstraction thermodynamic cycles. This includes a separation of one-electron reduction from protonation/deprotonation through the temperature dependence. The accuracy of the method has been assessed for the set of 15 transition-metal complexes that are considered as highly challenging systems for computational electrochemistry. Unlike the standard computational protocol(s), the presented approach yields results that are in excellent agreement with experimental electrochemical data. Last but not least, the applicability and limitations of the approach are thoroughly discussed.

A quantum chemical study of title systems and some related species (CaO, SrO, BaO and BaPo) was performed by B3LYP, H.F., MP2, TD DFT, and SAC-CI procedures. Structural features of isomers, vibrational frequencies, population analysis, ionization potential and in particular, electronic spectra. Electronic spectra of several monomers exhibit in accordance with experiments, the first electronic transitions in the near IR region (SAC-CI and TD DFT techniques). The passage to oligomers is accompanied by a very significant hypsochromic shift of the longest wavelength bands, increase of ionization potential and decrease of electron affinity.

Electronic spectra of conjugated polyynes [H(C≡C)nTH, nT = 1 to 7] and cumulenes [H2(C)nCH2, nC = 1 to 9] were calculated by means of the time-dependent density functional theory (TD DFT) and, for a group of selected molecules, also by the symmetry-adapted cluster configuration interaction method (SAC-CI). A comparison was made between calculated and published experimental spectral data. It turned out that the TD DFT (with B3LYP) was a reliable tool for calculation of band positions in the longest-wavelength region. In the case of cumulenes, except allene, only data for alkyl and phenyl derivatives are available. The characteristics of cumulenes split into two distinctly separate classes, planar (with an even number of carbon atoms, D2h) and non-planar (with an odd number of carbon atoms, D2d). Special attention was paid to the influence of substituents of various types on the position of the first intensive bands of polyynes and polyenes. Plotting wavenumbers of these bands against the reciprocal number of the CC triple bonds (polyynes) or the CC double bonds (polyenes) lead to (partial) rectification and made the extrapolation to infinite chains easier. While the extrapolated value for infinite polyynes (and derivatives like dimethyl or dicyano) amounts to about 16 000 cm-1, for anions and radical anions derived from the parent polyynes, the extrapolated value tends to significantly smaller values or possibly to zero. The situation with polyenes and their derivatives and heteroanalogues is similar; small deviations from zero with chains without significant bond alternation seem likely.

The results of density matrix renormalization group complete active space self-consistent field (DMRG-CASSCF) and second-order perturbation theory (DMRG-CASPT2) calculations are presented on various structural alternatives for the O-O and first C-H activating step of the catalytic cycle of the binuclear nonheme iron enzyme Δ(9) desaturase. This enzyme is capable of inserting a double bond into an alkyl chain by double hydrogen (H) atom abstraction using molecular O2. The reaction step studied here is presumably associated with the highest activation barrier along the full pathway; therefore, its quantitative assessment is of key importance to the understanding of the catalysis. The DMRG approach allows unprecedentedly large active spaces for the explicit correlation of electrons in the large part of the chemically important valence space, which is apparently conditio sine qua non for obtaining well-converged reaction energetics. The derived reaction mechanism involves protonation of the previously characterized 1,2-μ peroxy Fe(III)Fe(III) (P) intermediate to a 1,1-μ hydroperoxy species, which abstracts an H atom from the C10 site of the substrate. An Fe(IV)-oxo unit is generated concomitantly, supposedly capable of the second H atom abstraction from C9. In addition, several popular DFT functionals were compared to the computed DMRG-CASPT2 data. Notably, many of these show a preference for heterolytic C-H cleavage, erroneously predicting substrate hydroxylation. This study shows that, despite its limitations, DMRG-CASPT2 is a significant methodological advancement toward the accurate computational treatment of complex bioinorganic systems, such as those with the highly open-shell diiron active sites.

Mononuclear non-heme Fe (NHFe) enzymes play key roles in DNA repair, the biosynthesis of antibiotics, the response to hypoxia, cancer therapy, and many other biological processes. These enzymes catalyze a diverse range of oxidation reactions, including hydroxylation, halogenation, ring closure, desaturation, and electrophilic aromatic substitution (EAS). Most of these enzymes use an Fe(II) site to activate dioxygen, but traditional spectroscopic methods have not allowed researchers to insightfully probe these ferrous active sites. We have developed a methodology that provides detailed geometric and electronic structure insights into these NHFe(II) active sites. Using these data, we have defined a general mechanistic strategy that many of these enzymes use: they control O2 activation (and limit autoxidation and self-hydroxylation) by allowing Fe(II) coordination unsaturation only in the presence of cosubstrates. Depending on the type of enzyme, O2 activation either involves a 2e(-) reduced Fe(III)-OOH intermediate or a 4e(-) reduced Fe(IV)═O intermediate. Nuclear resonance vibrational spectroscopy (NRVS) has provided the geometric structure of these intermediates, and magnetic circular dichroism (MCD) has defined the frontier molecular orbitals (FMOs), the electronic structure that controls reactivity. This Account emphasizes that experimental spectroscopy is critical in evaluating the results of electronic structure calculations. Therefore these data are a key mechanistic bridge between structure and reactivity. For the Fe(III)-OOH intermediates, the anticancer drug activated bleomycin (BLM) acts as the non-heme Fe analog of compound 0 in heme (e.g., P450) chemistry. However BLM shows different reactivity: the low-spin (LS) Fe(III)-OOH can directly abstract a H atom from DNA. The LS and high-spin (HS) Fe(III)-OOHs have fundamentally different transition states. The LS transition state goes through a hydroxyl radical, but the HS transition state is activated for EAS without O-O cleavage. This activation is important in one class of NHFe enzymes that utilizes a HS Fe(III)-OOH intermediate in dioxygenation. For Fe(IV)═O intermediates, the LS form has a π-type FMO activated for attack perpendicular to the Fe-O bond. However, the HS form (present in the NHFe enzymes) has a π FMO activated perpendicular to the Fe-O bond and a σ FMO positioned along the Fe-O bond. For the NHFe enzymes, the presence of π and σ FMOs enables enzymatic control in determining the type of reactivity: EAS or H-atom extraction for one substrate with different enzymes and halogenation or hydroxylation for one enzyme with different substrates.

S = 2 Fe IV ═O species are key intermediates in the catalysis of most nonheme iron enzymes. This article presents detailed spectroscopic and high-level computational studies on a structurally-defined S = 2 Fe IV ═O species that define its frontier molecular orbitals, which allow its high reactivity. Importantly, there are both π- and σ-channels for reaction, and both are highly reactive because they develop dominant oxyl character at the transition state. These π- and σ-channels have different orientation dependences defining how the same substrate can undergo different reactions (H-atom abstraction vs. electrophilic aromatic attack) with Fe IV ═O sites in different enzymes, and how different substrates can undergo different reactions (hydroxylation vs. halogenation) with an Fe IV ═O species in the same enzyme.