Two-State Reactivity in Alkane Hydroxylation by Non-Heme Iron−Oxo Complexes

Proton-triggered chemoselective halogenation of aliphatic C-H bonds with nonheme FeIV-oxo complexes

Halogenation of aliphatic C-H bonds is a chemical transformation performed in nature by mononuclear nonheme iron dependent halogenases. The mechanism involves the formation of an iron(IV)-oxo-chloride species that abstracts the hydrogen atom from the reactive C-H bond to form a carbon-centered radical that selectively reacts with the bound chloride ligand, a process commonly referred to as halide rebound. The factors that determine the halide rebound, as opposed to the reaction with the incipient hydroxide ligand, are not clearly understood and examples of well-defined iron(IV)-oxo-halide compounds competent in C-H halogenation are scarce. In this work we have studied the reactivity of three well-defined iron(IV)-oxo complexes containing variants of the tetradentate 1-(2-pyridylmethyl)-1,4,7-triazacyclononane ligand (Pytacn). Interestingly, these compounds exhibit a change in their chemoselectivity towards the functionalization of C-H bonds under certain conditions: their reaction towards C-H bonds in the presence of a halide anionleads to exclusive oxygenation, while the addition of a superacid results in halogenation. Almost quantitative halogenation of ethylbenzene is observed when using the two systems with more sterically congested ligands and even the chlorination of strong C-H bonds such as those of cyclohexane is performed when a methyl group is present in the sixth position of the pyridine ring of the ligand. Mechanistic studies suggest that both reactions, oxygenation and halogenation, proceed through a common rate determining hydrogen atom transfer step and the presence of the acid dictates the fate of the resulting alkyl radical towards preferential halogenation over oxygenation.

Enhancing the high-spin reactivity in C-H bond activation by Iron (IV)-Oxo species: insights from paclitaxel hydroxylation by CYP2C8

Previous theoretical studies have revealed that high-spin states possess flatter potential energy surfaces than low-spin states in reactions involving iron(IV)-oxo species of cytochrome P450 enzymes (P450s), nonheme enzymes, or biomimetic complexes. Therefore, actively utilizing high-spin states to enhance challenging chemical transformations, such as C-H bond activation, represents an intriguing research avenue. However, the inherent instability of high-spin states relative to low-spin states in pre-reaction complexes often hinders their accessibility around the transition state, especially in heme systems with strong ligand fields. Counterintuitively, our investigation of the metabolic hydroxylation of paclitaxel by human CYP2C8 using a hybrid quantum mechanics and molecular mechanics (QM/MM) approach showed that the high-spin sextet state exhibits unusually high stability, when the reaction follows a secondary reaction pathway leading to 6β-hydroxypaclitaxel. We thoroughly analyzed the factors contributing to the enhanced stabilization of the high-spin state, and the knowledge obtained could be instrumental in designing competent biomimetic catalysts and biocatalysts for C-H bond activation.

Unraveling the Geometry-Driven C═C Epoxidation and C-H Hydroxylation Reactivity of Tetra-Coordinated Nonheme Iron(IV)-Oxo Complexes

The electronic structure and reactivity of tetra-coordinated nonheme iron(IV)-oxo complexes have remained unexplored for years. The recent synthesis of a closed-shell iron(IV)-oxo complex [(quinisox)FeIV(O)]+ (1) has set up a platform to understand how such complexes compare with the celebrated open-shell iron-oxo chemistry. Herein, using density functional theory and ab initio calculations, we present an in-depth electronic structure investigation of the C═C epoxidation [oxygen atom transfer (OAT)] and C-H hydroxylation [hydrogen atom transfer (HAT)] reactivity of 1. Using a solvent-coordinated geometry of 1 (1') and other potential tetra-coordinated iron(IV)-oxo complexes bearing rigid ligands (2 and 3), we established the geometric origin of spin-state energetics and reactivity of 1. Complex 1 featuring a strong Fe-O bond exhibits OAT and HAT reactivity in its quintet state. The lowest quintet OAT pathway has a lower barrier by ∼4 kcal/mol than the quintet HAT pathway, corroborating the experimentally observed gas-phase OAT reactivity preference. A conventional HAT reactivity preference for 2 and a comparable OAT and HAT reactivity for 3 are observed. This further supports the geometry-driven reactivity preference for 1. Noncovalent interaction analyses reveal a pronounced π-π interaction between the substrate and ligand in the OAT transition state, rationalizing the origin of the observed reactivity preference for 1.

The spin-forbidden transition in iron(IV)-oxo catalysts relevant to two-state reactivity

Quintet oxoiron(IV) intermediates are often invoked in nonheme iron enzymes capable of performing selective oxidation, while most well-characterized synthetic model oxoiron(IV) complexes have a triplet ground state. These differing spin states lead to the proposal of a two-state reactivity model, where the complexes cross from the triplet to an excited quintet state. However, the energy of this quintet state has never been measured experimentally. Here, magnetic circular dichroism is used to assign the singlet and triplet excited states in a series of triplet oxoiron(IV) complexes. These transition energies are used to determine the energies of the quintet state via constrained fitting of 2p3d resonant inelastic x-ray scattering. This allowed for a direct correlation between the quintet energies and substrate C─H oxidation rates.

An Active Site Tyr Residue Guides the Regioselectivity of Lysine Hydroxylation by Nonheme Iron Lysine-4-hydroxylase Enzymes through Proton-Coupled Electron Transfer

Lysine dioxygenase (KDO) is an important enzyme in human physiology involved in bioprocesses that trigger collagen cross-linking and blood pressure control. There are several KDOs in nature; however, little is known about the factors that govern the regio- and stereoselectivity of these enzymes. To understand how KDOs can selectively hydroxylate their substrate, we did a comprehensive computational study into the mechanisms and features of 4-lysine dioxygenase. In particular, we selected a snapshot from the MD simulation on KDO5 and created large QM cluster models (A, B, and C) containing 297, 312, and 407 atoms, respectively. The largest model predicts regioselectivity that matches experimental observation with rate-determining hydrogen atom abstraction from the C4-H position, followed by fast OH rebound to form 4-hydroxylysine products. The calculations show that in model C, the dipole moment is positioned along the C4-H bond of the substrate and, therefore, the electrostatic and electric field perturbations of the protein assist the enzyme in creating C4-H hydroxylation selectivity. Furthermore, an active site Tyr233 residue is identified that reacts through proton-coupled electron transfer akin to the axial Trp residue in cytochrome c peroxidase. Thus, upon formation of the iron(IV)-oxo species in the catalytic cycle, the Tyr233 phenol loses a proton to the nearby Asp179 residue, while at the same time, an electron is transferred to the iron to create an iron(III)-oxo active species. This charged tyrosyl residue directs the dipole moment along the C4-H bond of the substrate and guides the selectivity to the C4-hydroxylation of the substrate.

C-H Bond Oxidation by MnIV-Oxo Complexes: Hydrogen-Atom Tunneling and Multistate Reactivity

The reactivity of six MnIV-oxo complexes in C-H bond oxidation has been examined using a combination of kinetic experiments and computational methods. Variable-temperature studies of the oxidation of 9,10-dihydroanthracene (DHA) and ethylbenzene by these MnIV-oxo complexes yielded activation parameters suitable for evaluating electronic structure computations. Complementary kinetic experiments of the oxidation of deuterated DHA provided evidence for hydrogen-atom tunneling in C-H bond oxidation for all MnIV-oxo complexes. These results are in accordance with the Bell model, where tunneling occurs near the top of the transition-state barrier. Density functional theory (DFT) and DLPNO-CCSD(T1) computations were performed for three of the six MnIV-oxo complexes to probe a previously predicted multistate reactivity model. The DFT computations predicted a thermal crossing from the 4B1 ground state to a 4E state along the C-H bond oxidation reaction coordinate. DLPNO-CCSD(T1) calculations further confirm that the 4E transition state offers a lower energy barrier, reinforcing the multistate reactivity model for these complexes. We discuss how this multistate model can be reconciled with recent computations that revealed that the kinetics of C-H bond oxidation by this set of MnIV-oxo complexes can be well-predicted on the basis of the thermodynamic driving force for these reactions.

Computational Study Into the Oxidative Ring-Closure Mechanism During the Biosynthesis of Deoxypodophyllotoxin

The nonheme iron dioxygenase deoxypodophyllotoxin synthase performs an oxidative ring-closure reaction as part of natural product synthesis in plants. How the enzyme enables the oxidative ring-closure reaction of (-)-yatein and avoids substrate hydroxylation remains unknown. To gain insight into the reaction mechanism and understand the details of the pathways leading to products and by-products we performed a comprehensive computational study. The work shows that substrate is bound tightly into the substrate binding pocket with the C7'-H bond closest to the iron(IV)-oxo species. The reaction proceeds through a radical mechanism starting with hydrogen atom abstraction from the C7'-H position followed by ring-closure and a final hydrogen transfer to form iron(II)-water and deoxypodophyllotoxin. Alternative mechanisms including substrate hydroxylation and an electron transfer pathway were explored but found to be higher in energy. The mechanism is guided by electrostatic perturbations of charged residues in the second-coordination sphere that prevent alternative pathways.

DFT and TDDFT exploration on electronic transitions and bonding aspect of DPA and PTDC ligated transition metal complexes

CONTEXT

In this study, we have investigated the structure, reactivity, bonding, and electronic transitions of DPA and PDTC along with their Ni-Zn complexes using DFT/TD-DFT methods. The energy gap between the frontier orbitals was computed to understand the reactivity pattern of the ligands and metal complexes. From the energies of FMO's, the global reactivity descriptors such as electron affinity, ionization potential, hardness (η), softness (S), chemical potential (μ), electronegativity (χ), and electrophilicity index (ω) have been calculated. The complexes show a strong NLO properties due to easily polarization as indicated by the narrow HOMO-LUMO gap. The polarizability and hyperpolarizabilities of the complexes indicate that they are good candidates for NLO materials. Molecular electrostatic potential (MEP) maps identified electrophilic and nucleophilic sites on the surfaces of the complexes. TDDFT and NBO analyses provided insights into electronic transitions, bonding, and stabilizing interactions within the studied complexes. DPA and PDTC exhibited larger HOMO-LUMO gaps and more negative electrostatic potentials compared to their metal complexes suggesting the higher reactivity. Ligands (DPA and PDTC) had absorption spectra in the range of 250 nm to 285 nm while their complexes spanned 250 nm to 870 nm. These bands offer valuable information on electronic transitions, charge transfer and optical behavior. This work enhances our understanding of the electronic structure and optical properties of these complexes.

METHODS

Gaussian16 program was used for the optimization of all the compounds. B3LYP functional in combination with basis sets, such as LanL2DZ for Zn, Ni and Cu while 6-311G** for other atoms like C, H, O, N, and S was used. Natural bond orbital (NBO) analysis is carried out to find out how the filled orbital of one sub-system interacts with the empty orbital of another sub-system. The ORCA software is used for computing spectral features along with the zeroth order regular approximation method (ZORA) to observe its relativistic effects. TD-DFT study is carried out to calculate the excitation energy by using B3LYP functional.

How Do Differences in Electronic Structure Affect the Use of Vanadium Intermediates as Mimics in Nonheme Iron Hydroxylases?

We study active-site models of nonheme iron hydroxylases and their vanadium-based mimics using density functional theory to determine if vanadyl is a faithful structural mimic. We identify crucial structural and energetic differences between ferryl and vanadyl isomers owing to the differences in their ground electronic states, i.e., high spin (HS) for Fe and low spin (LS) for V. For the succinate cofactor bound to the ferryl intermediate, we predict facile interconversion between monodentate and bidentate coordination isomers for ferryl species but difficult rearrangement for vanadyl mimics. We study isomerization of the oxo intermediate between axial and equatorial positions and find the ferryl potential energy surface to be characterized by a large barrier of ca. 10 kcal/mol that is completely absent for the vanadyl mimic. This analysis reveals even starker contrasts between Fe and V in hydroxylases than those observed for this metal substitution in nonheme halogenases. Analysis of the relative bond strengths of coordinating carboxylate ligands for Fe and V reveals that all of the ligands show stronger binding to V than Fe owing to the LS ground state of V in contrast to the HS ground state of Fe, highlighting the limitations of vanadyl mimics of native nonheme iron hydroxylases.

How Axial Coordination Regulates the Electronic Structure and C-H Amination Reactivity of Fe-Porphyrin-Nitrene?

Detailed electronic structure and its correlation with the intramolecular C-H amination reactivity of Fe-porphyrin-nitrene intermediates bearing different "axial" coordination have been investigated using multiconfigurational complete active space self-consistent field (CASSCF), N-electron valence perturbation theory (NEVPT2), and hybrid density functional theory (DFT-B3LYP) calculations. Three types of "axial" coordination, -OMe/-O(H)Me (1-Sul/2-Sul), -SMe/-S(H)Me (3-Sul/4-Sul), and -NMeIm (MeIm = 3-methyl-imidazole) (5-Sul) mimicking serine, cysteine, and histidine, respectively, along with no axial coordination (6-Sul) have been considered to decipher how the "axial" coordination of different strengths regulates the electronic integrity of the Fe-N core and nitrene-transfer reactivity of Fe-porphyrin-nitrene intermediates. CASSCF-based natural orbitals reveal two distinct classes of electronic structures: Fe-nitrenes (1-Sul and 3-Sul) with relatively stronger axial coordination (-OMe and -SMe) display "imidyl" nature and those (2-Sul, 4-Sul, and 6-Sul) with weaker axial coordination (-O(H)Me, -S(H)Me and no axial coordination) exhibit "imido-like" character. A borderline between the two classes is also observed with NMeIm axial coordination (5-Sul). Axial coordination of different strengths not only regulates the electronic structure but also modulates the Fe-3d orbital energies, as revealed through the d-d transition energies obtained by CASSCF/NEVPT2 calculations. The relatively lower energy of Fe-3dz2 orbital allows easy access to low-lying high-spin quintet states in the cases of weaker "axial" coordination (2-Sul, 4-Sul, and 6-Sul), and the associated hydrogen atom transfer (HAT) reactivity appears to involve two-state triplet-quintet reactivity through minimum energy crossing point (3,5MECP) between the spin states. In stark contrast, Fe-nitrenes with relatively stronger "axial" coordination (1-Sul and 3-Sul) undergo triplet-only HAT reactivity. Overall, this in-depth electronic structure investigation and HAT reactivity evaluation reveal that the weaker axial coordination in Fe-porphyrin-nitrene complexes (2-Sul, 4-Sul, and 6-Sul) can promote more efficient C-H oxidation through the quintet spin state.

The Nature of the Chemical Bonds of High-Valent Transition-Metal Oxo (M=O) and Peroxo (MOO) Compounds: A Historical Perspective of the Metal Oxyl-Radical Character by the Classical to Quantum Computations

This review article describes a historical perspective of elucidation of the nature of the chemical bonds of the high-valent transition metal oxo (M=O) and peroxo (M-O-O) compounds in chemistry and biology. The basic concepts and theoretical backgrounds of the broken-symmetry (BS) method are revisited to explain orbital symmetry conservation and orbital symmetry breaking for the theoretical characterization of four different mechanisms of chemical reactions. Beyond BS methods using the natural orbitals (UNO) of the BS solutions, such as UNO CI (CC), are also revisited for the elucidation of the scope and applicability of the BS methods. Several chemical indices have been derived as the conceptual bridges between the BS and beyond BS methods. The BS molecular orbital models have been employed to explain the metal oxyl-radical character of the M=O and M-O-O bonds, which respond to their radical reactivity. The isolobal and isospin analogy between carbonyl oxide R2C-O-O and metal peroxide LFe-O-O has been applied to understand and explain the chameleonic chemical reactivity of these compounds. The isolobal and isospin analogy among Fe=O, O=O, and O have also provided the triplet atomic oxygen (3O) model for non-heme Fe(IV)=O species with strong radical reactivity. The chameleonic reactivity of the compounds I (Cpd I) and II (Cpd II) is also explained by this analogy. The early proposals obtained by these theoretical models have been examined based on recent computational results by hybrid DFT (UHDFT), DLPNO CCSD(T0), CASPT2, and UNO CI (CC) methods and quantum computing (QC).

Probing the Origins of Puzzling Reactivity in Fe/Mn-Oxo/Hydroxo Species toward C-H Bonds: A DFT and Ab Initio Perspective

Activation of C-H bonds using an earth-abundant metal catalyst is one of the top challenges of chemistry, where high-valent Mn/Fe-oxo(hydroxo) biomimic species play an important role. There are several open questions related to the comparative oxidative abilities of these species, and a unifying concept that could accommodate various factors influencing reactivity is lacking. To shed light on these open questions, here, we have used a combination of density functional theory (DFT) (B3LYP-D3/def2-TZVP) and ab initio (CASSCF/NEVPT2) calculations to study a series of high-valent metal-oxo species [Mn+H3buea(O/OH)] (M = Mn and Fe, n = II to V; H3buea = tris[(N'-tert-butylureaylato)-N-ethylene)]aminato towards the activation of dihydroanthracene (DHA). The H-bonding network in the ligand architecture influences the ground state-excited state gap and brings several excited states of the same spin multiplicity closer in energy, which triggers reactivity via one of those excited states, reducing the kinetic barriers for the C-H bond activation and rationalizing several puzzling reactivity trends observed in various high-valent Mn/Fe-oxo(hydroxo) species.

Two state reactivity (TSR) and hydrogen tunneling reaction kinetics measured in the Co+ mediated decomposition of CH3CHO

The electronic structure of a transition metal atom allows it to act as a catalytic active site by providing lower energy alternative pathways in chemical transformations. We have identified and kinetically characterized three such pathways in the title reaction. One is an adiabatic pathway that occurs on a single potential energy surface described within the Born-Oppenheimer approximation. A second pathway opens microseconds into the reaction as a portion of the reacting population competitively transitions from triplet to singlet multiplicity to circumvent energetic barriers on the triplet surface. These pathways are single- and two-state reactive (SSR and TSR) where the Co+ cation mediates an oxidative addition/reductive elimination sequence of the CH3CHO molecule. The third observed reaction pathway is the aldehyde hydrogen tunneling through an Eyring barrier to form high-spin products. First-order rate constants for the adiabatic and nonadiabatic energy lowered pathways, and the hydrogen tunneling pathway, are each measured using the single photon initiated dissociative rearrangement reaction (SPIDRR) experimental technique. We believe that this is the first experimental study where such disparate dynamic features (SSR, TSR, and H-tunneling) are disentangled in a system's chemistry, attributing specific rate constant values to each effect and quantifying the various competitions. Moreover, multi-reference CASSCF/CASPT2 calculations indicate that structures with covalent Co-H bonds are present exclusively along the excited singlet surface. This phenomenon significantly reduces these structures' energy relative to their triplet counterparts, thus enabling the surface crossing and spin inversion that cause the observed two-state reactivity.

Two-State Hydrogen Atom Transfer Reactivity of Unsymmetric [Cu2(O)(NO)]2+ Complexes

We report the temperature-dependent spin switching of dicopper oxo nitrosyl [Cu2(O)(NO)]2+ complexes and their influence on hydrogen atom transfer (HAT) reactivity. Electron paramagnetic resonance (EPR) and Evans method analysis suggest that [Cu2(O)(NO)]2+ complexes transition from the S = 1/2 to the S = 3/2 state around ca. 202 K. At low temperatures (198 K) where S = 3/2 dominates, a strong correlation between the rate of HAT (kHAT) and the population of the S = 1/2 state was identified (R2 = 0.988), suggesting that the HAT by [Cu2(O)(NO)]2+ complexes proceeds by the S = 1/2 isomer. Installation of functional groups that introduce an unsymmetric secondary coordination environment accelerates the HAT rates through perturbation of the spin equilibria. Given the often unsymmetric coordination sphere of bimetallic active sites in natural proteins, we anticipate that similar strategies could be employed by metalloenzymes to control HAT reactions.

Single-Atom MnN5 Catalytic Sites Enable Efficient Peroxymonosulfate Activation by Forming Highly Reactive Mn(IV)-Oxo Species

Four-nitrogen-coordinated transitional metal (MN₄) configurations in single-atom catalysts (SACs) are broadly recognized as the most efficient active sites in peroxymonosulfate (PMS)-based advanced oxidation processes. However, SACs with a coordination number higher than four are rarely explored, which represents a fundamental missed opportunity for coordination chemistry to boost PMS activation and degradation of recalcitrant organic pollutants. We experimentally and theoretically demonstrate here that five-nitrogen-coordinated Mn (MnN₅) sites more effectively activate PMS than MnN₄ sites, by facilitating the cleavage of the O-O bond into high-valent Mn(IV)-oxo species with nearly 100% selectivity. The high activity of MnN₅ was discerned to be due to the formation of higher-spin-state N₅Mn(IV)═O species, which enable efficient two-electron transfer from organics to Mn sites through a lower-energy-barrier pathway. Overall, this work demonstrates the importance of high coordination numbers in SACs for efficient PMS activation and informs the design of next-generation environmental catalysts.

Comparing the reactivity of an oxoiron(IV) cation radical and its oxoiron(V) tautomer towards C-H bonds

An oxoiron(IV) cation radical is generated upon two-electron oxidation of an iron(III) complex bearing an electron-rich methoxy substituted bTAML framework and thoroughly characterized via multiple spectroscopic techniques and density functional theory (DFT). Reactivity studies demonstrate faster rates for oxidation of strong aliphatic sp³ C-H bonds than for its corresponding oxoiron(V) valence tautomer.

The role of the intermediate triplet state in iron-catalyzed multi-state C-H activation

Efficient activation and functionalization of the C-H bond under mild conditions are of a great interest in chemical synthesis. We investigate the previously proposed spin-accelerated activation of the C(sp²)-H bond by a Fe(II)-based catalyst to clarify the role of the intermediate triplet state in the reaction mechanism. High-level electronic structure calculations on a small model of a catalytic system utilizing the coupled cluster with the single, double, and perturbative triple excitations [CCSD(T)] are used to select the density functional for the full-size model. Our analysis indicates that the previously proposed two-state quintet-singlet reaction pathway is unlikely to be efficient due to a very weak spin-orbit coupling between these two spin states. We propose a more favorable multi-state quintet-triplet-singlet reaction pathway and discuss the importance of the intermediate triplet state. This triplet state facilitates a spin-accelerated reaction mechanism by strongly coupling to both quintet and singlet states. Our calculations show that the C-H bond activation through the proposed quintet-triplet-singlet reaction pathway is more thermodynamically favorable than the single-state quintet and two-state singlet-quintet mechanisms.

What Drives Radical Halogenation versus Hydroxylation in Mononuclear Nonheme Iron Complexes? A Combined Experimental and Computational Study

Nonheme iron halogenases are unique enzymes in nature that selectively activate an aliphatic C–H bond of a substrate to convert it into C–X (X = Cl/Br, but not F/I). It is proposed that they generate an Fe^III(OH)(X) intermediate in their catalytic cycle. The analogous Fe^III(OH) intermediate in nonheme iron hydroxylases transfers OH^• to give alcohol product, whereas the halogenases transfer X^• to the carbon radical substrate. There remains significant debate regarding what factors control their remarkable selectivity of the halogenases. The reactivity of the complexes Fe^III(BNPA^Ph2O)(OH)(X) (X = Cl, Br) with a secondary carbon radical (R^•) is described. It is found that X^• transfer occurs with a secondary carbon radical, as opposed to OH^• transfer with tertiary radicals. Comprehensive computational studies involving density functional theory were carried out to examine the possible origins of this selectivity. The calculations reproduce the experimental findings, which indicate that halogen transfer is not observed for the tertiary radicals because of a nonproductive equilibrium that results from the endergonic nature of these reactions, despite a potentially lower reaction barrier for the halogenation pathway. In contrast, halogen transfer is favored for secondary carbon radicals, for which the halogenated product complex is thermodynamically more stable than the reactant complex. These results are rationalized by considering the relative strengths of the C–X bonds that are formed for tertiary versus secondary carbon centers. The computational analysis also shows that the reaction barrier for halogen transfer is significantly dependent on secondary coordination sphere effects, including steric and H-bonding interactions.

Electrostatic Perturbations in the Substrate-Binding Pocket of Taurine/α-Ketoglutarate Dioxygenase Determine its Selectivity

Taurine/α-ketoglutarate dioxygenase is an important enzyme that takes part in the cysteine catabolism process in the human body and selectively hydroxylates taurine at the C1 -position. Recent computational studies showed that in the gas-phase the C2 -H bond of taurine is substantially weaker than the C1 -H bond, yet no evidence exists of 2-hydroxytaurine products. To this end, a detailed computational study on the selectivity patterns in TauD was performed. The calculations show that the second-coordination sphere and the protonation states of residues play a major role in guiding the enzyme to the right selectivity. Specifically, a single proton on an active site histidine residue can change the regioselectivity of the reaction through its electrostatic perturbations in the active site and effectively changes the C1 -H and C2 -H bond strengths of taurine. This is further emphasized by many polar and hydrogen bonding interactions of the protein cage in TauD with the substrate and the oxidant that weaken the pro-R C1 -H bond and triggers a chemoselective reaction process. The large cluster models reproduce the experimental free energy of activation excellently.

Measurement of time dependent product branching ratios indicates two-state reactivity in metal mediated chemical reactions

For several decades, the influence of Two State Reactivity (TSR) has been implicated in a host of reactions, but has lacked a stand-alone, definitive experimental kinetic signature identifying its occurrence. Here, we demonstrate that the measurement of a temporally dependent product branching ratio is indicative of spin inversion and is a kinetic signature of TSR. This is caused by products exiting different hypersurfaces with different rates and relative exothermicities. The composite measurement of product intensities with the same mass but with different multiplicities yield biexponential temporal dependences with the sampled product ratio changing in time. These measurements are made using the single photon initiated dissociative rearrangement reaction (SPIDRR) technique which identifies TSR but further determines the kinetic parameters for reaction along the original ground electronic surface in competition with spin inversion and its consequent TSR.

A radical mechanism for C–H bond cross-coupling and N2 activation catalysed by β-diketiminate iron complexes

Density functional theory calculations and electronic structure analyses reveal a radical mechanism with spin-crossovers for C–H bond cross-coupling and N2 activation catalysed by β-diketiminate iron complexes.

Can the isonitrile biosynthesis enzyme ScoE assist with the biosynthesis of isonitrile groups in drug molecules? A computational study

Computational studies show that the isonitrile synthesizing enzyme ScoE can catalyse the conversion of γ-Gly substituents in substrates to isonitrile. This enables efficient isonitrile substitution into target molecules such as axisonitrile-1.

Large Isotope Effects in Organometallic Chemistry

The kinetic isotope effect (KIE) is key to understanding reaction mechanisms in many areas of chemistry and chemical biology, including organometallic chemistry. This ratio of rate constants, k_H /k_D , typically falls between 1-7. However, KIEs up to 105 have been reported, and can even be so large that reactivity with deuterium is unobserved. We collect here examples of large KIEs across organometallic chemistry, in catalytic and stoichiometric reactions, along with their mechanistic interpretations. Large KIEs occur in proton transfer reactions such as protonation of organometallic complexes and clusters, protonolysis of metal-carbon bonds, and dihydrogen reactivity. C-H activation reactions with large KIEs occur with late and early transition metals, photogenerated intermediates, and abstraction by metal-oxo complexes. We categorize the mechanistic interpretations of large KIEs into the following three types: (a) proton tunneling, (b) compound effects from multiple steps, and (c) semi-classical effects on a single step. This comprehensive collection of large KIEs in organometallics provides context for future mechanistic interpretation.

Structure and Functional Differences of Cysteine and 3-Mercaptopropionate Dioxygenases: A Computational Study

Thiol dioxygenases are important enzymes for human health; they are involved in the detoxification and catabolism of toxic thiol-containing natural products such as cysteine. As such, these enzymes have relevance to the development of Alzheimer's and Parkinson's diseases in the brain. Recent crystal structure coordinates of cysteine and 3-mercaptopropionate dioxygenase (CDO and MDO) showed major differences in the second-coordination spheres of the two enzymes. To understand the difference in activity between these two analogous enzymes, we created large, active-site cluster models. We show that CDO and MDO have different iron(III)-superoxo-bound structures due to differences in ligand coordination. Furthermore, our studies show that the differences in the second-coordination sphere and particularly the position of a positively charged Arg residue results in changes in substrate positioning, mobility and enzymatic turnover. Furthermore, the substrate scope of MDO is explored with cysteinate and 2-mercaptosuccinic acid and their reactivity is predicted.

pH Changes That Induce an Axial Ligand Effect on Nonheme Iron(IV) Oxo Complexes with an Appended Aminopropyl Functionality

Nonheme iron enzymes often utilize a high-valent iron(IV) oxo species for the biosynthesis of natural products, but their high reactivity often precludes structural and functional studies of these complexes. In this work, a combined experimental and computational study is presented on a biomimetic nonheme iron(IV) oxo complex bearing an aminopyridine macrocyclic ligand and its reactivity toward olefin epoxidation upon changes in the identity and coordination ability of the axial ligand. Herein, we show a dramatic effect of the pH on the oxygen-atom-transfer (OAT) reaction with substrates. In particular, these changes have occurred because of protonation of the axial-bound pendant amine group, where its coordination to iron is replaced by a solvent molecule or anionic ligand. This axial ligand effect influences the catalysis, and we observe enhanced cyclooctene epoxidation yields and turnover numbers in the presence of the unbound protonated pendant amine group. Density functional theory studies were performed to support the experiments and highlight that replacement of the pendant amine with a neutral or anionic ligand dramatically lowers the rate-determining barriers of cyclooctene epoxidation. The computational work further establishes that the change in OAT is due to electrostatic interactions of the pendant amine cation that favorably affect the barrier heights.

Thebaine is Selectively Demethylated by Thebaine 6-O-Demethylase and Codeine-3-O-demethylase at Distinct Binding Sites: A Computational Study

Two homologous 2-oxoglutarate-dependent (ODD) nonheme enzymes thebaine 6-O-demethylase (T6ODM) and codeine-3-O-demethylase (CODM), are involved in the morphine biosynthesis pathway from thebaine, catalyzing the O-demethylation reaction with precise regioselectivity at C6 and C3 positions of thebaine respectively. We investigated the origin of the regioselectivity of these enzymes by combined molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) calculations and found that Thebaine binds at the two distinct sites of T6ODM and CODM, which determines the regioselectivity of the enzymes. A remarkable oxo rotation is observed in the decarboxylation process. Starting from the closed pentacoordinate configuration, the C-terminal lid adopts an open conformation in the octahedral Fe(IV) = O complex to facilitate the subsequent demethylation. Phe241 and Phe311 stabilize the substrate in the binding pocket, while Arg219 acts as a gatekeeper residue to stabilize the substrate. Our results unravel the regioselectivity in 2-OG dependent nonheme enzymes and may shed light for exploring the substrate scope of these enzymes and developing novel biotechnology for morphine biosynthesis.

Density Functional Theory Study into the Reaction Mechanism of Isonitrile Biosynthesis by the Nonheme Iron Enzyme ScoE

The nonheme iron enzyme ScoE catalyzes the biosynthesis of an isonitrile substituent in a peptide chain. To understand details of the reaction mechanism we created a large active site cluster model of 212 atoms that contains substrate, the active oxidant and the first- and second-coordination sphere of the protein and solvent. Several possible reaction mechanisms were tested and it is shown that isonitrile can only be formed through two consecutive catalytic cycles that both use one molecule of dioxygen and α-ketoglutarate. In both cycles the active species is an iron(IV)-oxo species that in the first reaction cycle reacts through two consecutive hydrogen atom abstraction steps: first from the N–H group and thereafter from the C–H group to desaturate the NH-CH2 bond. The alternative ordering of hydrogen atom abstraction steps was also tested but found to be higher in energy. Moreover, the electronic configurations along that pathway implicate an initial hydride transfer followed by proton transfer. We highlight an active site Lys residue that is shown to donate charge in the transition states and influences the relative barrier heights and bifurcation pathways. A second catalytic cycle of the reaction of iron(IV)-oxo with desaturated substrate starts with hydrogen atom abstraction followed by decarboxylation to give isonitrile directly. The catalytic cycle is completed with a proton transfer to iron(II)-hydroxo to generate the iron(II)-water resting state. The work is compared with experimental observation and previous computational studies on this system and put in a larger perspective of nonheme iron chemistry.

Activation of C-H bonds by a nonheme iron(IV)-oxo complex: mechanistic evidence through a coupled EDXAS/UV-Vis multivariate analysis

The understanding of reactive processes involving organic substrates is crucial to chemical knowledge and requires multidisciplinary efforts for its advancement. Herein, we apply a combined multivariate, statistical and theoretical analysis of coupled time-resolved X-ray absorption (XAS)/UV-Vis data to obtain detailed mechanistic information for on the C-H bond activation of 9,10-dihydroanthracene (DHA) and diphenylmethane (Ph2CH2) by the nonheme FeIV-oxo complex [N4Py·FeIV(O)]2+ (N4Py = N,N-bis(2-pyridylmethyl)-N-bis(2-pyridyl)methylamine) in CH3CN at room temperature. Within this approach, we determine the number of key chemical species present in the reaction mixtures and derive spectral and concentration profiles for the reaction intermediates. From the quantitative analysis of the XAS spectra the transient intermediate species are structurally determined. As a result, it is suggested that, while DHA is oxidized by [N4Py·FeIV(O)]2+ with a hydrogen atom transfer-electron transfer (HAT-ET) mechanism, Ph2CH2 is oxidized by the nonheme iron-oxo complex through a HAT-radical dissociation pathway. In the latter process, we prove that the intermediate FeIII complex [N4Py·FeIII(OH)]2+ is not able to oxidize the diphenylmethyl radical and we provide its structural characterization in solution. The employed combined experimental and theoretical strategy is promising for the spectroscopic characterization of transient intermediates as well as for the mechanistic investigation of redox chemical transformations on the second to millisecond time scales.

A DFT study on the C-H oxidation reactivity of Fe(iv)-oxo species with N4/N5 ligands derived from l-proline

The hydroxylation of hexane by two FeIVO complexes bearing a pentadentate ligand (N5, Pro3Py) and a tetradentate ligand (N4, Pro2PyBn) derived from l-proline was studied by DFT calculations. Theoretical results predict that both FeIVO complexes hold triplet ground states. The hydrogen atom abstraction (HAA) processes by both FeIVO species proceed through a two-state reactivity, thus indicating that HAA occurs via a low-barrier quintet surface. Beyond the conventional rebound step, the dissociation path is also calculated and is found to potentially occur after HAA.

A mechanistic study of the manganese porphyrin-catalyzed C–H isocyanation reaction

The favourable radical rebound pathway is NCO-rebound from the Mn(TMP)(NCO)₂ complex due to the stronger trans effect of the axial ligand NCO and the electron-donating aryl substituents on the porphyrin ligand.

Computational studies of DNA base repair mechanisms by nonheme iron dioxygenases: selective epoxidation and hydroxylation pathways

DNA base repair mechanisms of alkylated DNA bases is an important reaction in chemical biology and particularly in the human body. It is typically catalyzed by an α-ketoglutarate-dependent nonheme iron dioxygenase named the AlkB repair enzyme. In this work we report a detailed computational study into the structure and reactivity of AlkB repair enzymes with alkylated DNA bases. In particular, we investigate the aliphatic hydroxylation and C[double bond, length as m-dash]C epoxidation mechanisms of alkylated DNA bases by a high-valent iron(iv)-oxo intermediate. Our computational studies use quantum mechanics/molecular mechanics methods on full enzymatic structures as well as cluster models on active site systems. The work shows that the iron(iv)-oxo species is rapidly formed after dioxygen binding to an iron(ii) center and passes a bicyclic ring structure as intermediate. Subsequent cluster models explore the mechanism of substrate hydroxylation and epoxidation of alkylated DNA bases. The work shows low energy barriers for substrate activation and consequently energetically feasible pathways are predicted. Overall, the work shows that a high-valent iron(iv)-oxo species can efficiently dealkylate alkylated DNA bases and return them into their original form.

Extending nudged elastic band method to reaction pathways involving multiple spin states

There are diverse reactions including spin-state crossing, especially the reactions catalyzed by transition metal compounds. To figure out the mechanisms of such reactions, the discussion of minimum energy intersystem crossing (MEISC) points cannot be avoided. These points may be the bottleneck of the reaction or inversely accelerate the reactions by providing a better pathway. It is of great importance to reveal their role in the reactions by computationally locating the position of the MEISC points together with the reaction pathway. However, providing a proper initial guess for the structure of the MEISC point is not as easy as that of the transition state. In this work, we extended the nudged elastic band (NEB) method for multiple spin systems, which is named the multiple spin-state NEB method, and it is successfully applied to find the MEISC points while optimizing the reaction pathway. For more precisely locating the MEISC point, a revised approach was adopted. Meanwhile, our examples also suggest that special attention should be paid to the criterion to define an image optimized as the MEISC point.

Spin-inversion mechanisms in O2 binding to a model heme compound: A perspective from nonadiabatic wave packet calculations

Spin-inversion dynamics in O₂ binding to a model heme complex, which consisted of Fe(II)-porphyrin and imidazole, were studied using nonadiabatic wave packet dynamics calculations. We considered three active nuclear degrees of freedom in the dynamics, including the motions along the Fe-O distance, Fe-O-O angle, and Fe out-of-plane distance. Spin-free potential energy surfaces for the singlet, triplet, quintet, and septet states were developed using density functional theory calculations, and spin-orbit coupling elements were obtained from CASSCF-level electronic structure calculations. The spin-inversion mainly occurred between the singlet state and one of the triplet states due to large spin-orbit couplings and the contributions of other states were extremely small. The present quantum dynamics calculations suggested that the narrow crossing region model plays a dominant role in the O₂ binding dynamics. In addition, the one-dimensional Landau-Zener model underestimated the nonadiabatic transition probability.

Iron and manganese oxo complexes, oxo wall and beyond

High-valent metal-oxo species with multiply-bonded M-O groups have been proposed as key intermediates in many biological and abiological catalytic oxidation reactions. These intermediates are implicated as active oxidants in alkane hydroxylation, olefin epoxidation and other oxidation reactions. For example, [Fe^ivO(porphyrinato^•-)]⁺ cofactors bearing π-radical porphyrinato^•- ligands oxidize organic substrates in cytochrome P450 enzymes, which are common to many life forms. Likewise, high-valent Mn-oxo species are active for H₂O oxidation in photosystem II. The chemistry of these native reactive species has inspired chemists to prepare highly oxidized transition-metal complexes as functional mimics. Although many synthetic Fe-O and Mn-O complexes now exist, the analogous oxo complexes of the late transition metals (groups 9-11) are rare. Indeed, late-transition-metal-oxo complexes of tetragonal (fourfold) symmetry should be electronically unstable, a rule commonly referred to as the 'oxo wall'. A few late metal-oxos have been prepared by targeting other symmetries or unusual spin states. These complexes have been studied using spectroscopic and theoretical methods. This Review describes mononuclear non-haem Fe-O and Mn-O species, the nature of the oxo wall and recent advances in the preparation of oxo complexes of Co, Ni and Cu beyond the oxo wall.