The valence bond way: reactivity patterns of cytochrome P450 enzymes and synthetic analogs

Machine learning-aided engineering of a cytochrome P450 for optimal bioconversion of lignin fragments

Using machine learning, molecular dynamics simulations, and density functional theory calculations we gain insight into the selectivity patterns of substrate activation by the cytochromes P450. In nature, the reactions catalyzed by the P450s lead to the biodegradation of xenobiotics, but recent work has shown that fungi utilize P450s for the activation of lignin fragments, such as monomer and dimer units. These fragments often are the building blocks of valuable materials, including drug molecules and fragrances, hence a highly selective biocatalyst that can produce these compounds in good yield with high selectivity would be an important step in biotechnology. In this work a detailed computational study is reported on two reaction channels of two P450 isozymes, namely the O-deethylation of guaethol by CYP255A and the O-demethylation versus aromatic hydroxylation of p-anisic acid by CYP199A4. The studies show that the second-coordination sphere plays a major role in substrate binding and positioning, heme access, and in the selectivity patterns. Moreover, the local environment affects the kinetics of the reaction through lowering or raising barrier heights. Furthermore, we predict a site-selective mutation for highly specific reaction channels for CYP199A4.

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.

High-valent nonheme Fe(IV)O/Ru(IV)O complexes catalyze C-H activation reactivity and hydrogen tunneling: a comparative DFT investigation

A comprehensive density functional theory investigation has been presented towards the comparison of the C-H activation reactivity between high-valent iron-oxo and ruthenium-oxo complexes. A total of four compounds, e.g., [Ru(IV)O(tpy-dcbpy)] (1), [Fe(IV)O(tpy-dcbpy)] (1'), [Ru(IV)O(TMCS)] (2), and [Fe(IV)O(TMCS)] (2'), have been considered for this investigation. The macrocyclic ligand framework tpy(dcbpy) implies tpy = 2,2':6',2''-terpyridine, dcbpy = 5,5'-dicarboxy-2,2'-bipyridine, and TMCS is TMC with an axially tethered -SCH2CH2 group. Compounds 1 and 2' are experimentally synthesized standard complexes with Ru and Fe, whereas compounds 1' and 2 were considered to keep the macrocycle intact when switching the central metal atom. Three reactants including benzyl alcohol, ethyl benzene, and dihydroanthracene were selected as substrates for C-H activation. It is noteworthy to mention that Fe(IV)O complexes exhibit higher reactivity than those of their Ru(IV)O counterparts. Furthermore, regardless of the central metal, the complex featuring a tpy-dcbpy macrocycle demonstrates higher reactivity than that of TMCS. Here, a thorough analysis of the reactivity-controlling characteristics-such as spin state, steric factor, distortion energy, energy of the electron acceptor orbital, and quantum mechanical tunneling-was conducted. Fe(IV)O exhibits the exchanged enhanced two-state-reactivity with the quintet reactive state, whereas Ru(IV)O has only a triplet reactive state. Both the distortion energy and acceptor orbital energy are low in the case of Fe(IV)O supporting its higher reactivity. All the investigated C-H activation processes involve a significant contribution from hydrogen tunneling, which is more pronounced in the case of Ru, although it cannot alter the reactivity pattern. Furthermore, it has also been found that, independent of the central metal, aliphatic hydroxylation is always preferable to aromatic hydroxylation. Overall, this work is successful in establishing and investigating the cause of enzymes' natural preference for Fe over Ru as a cofactor for C-H activation enzymes.

Influence of Equatorial Co-Ligands on the Reactivity of LFeIIIOIPh

Previous biomimetic studies clearly proved that equatorial ligands significantly influence the redox potential and thus the stability/reactivity of biologically important oxoiron intermediates; however, no such studies were performed on FeIIIOIPh species. In this study, the influence of substituted pyridine co-ligands on the reactivity of iron(III)-iodosylbenzene adduct has been investigated in sulfoxidation and epoxidation reactions. Selective oxidation of thioanisole, cis-cyclooctene, and cis- and trans-stilbene in the presence of a catalytic amount of [FeII(PBI)3](OTf)2 with PhI(OAc)2 provide products in good to excellent yields through an FeIIIOIPh intermediate depending on the co-ligand (4R-Py) used. Several mechanistic studies were performed to gain more insight into the mechanism of oxygen atom transfer (OAT) reactions to support the reactive intermediate and investigate the effect of the equatorial co-ligands. Based on competitive experiments, including a linear free-energy relationship between the relative reaction rates (logkrel) and the σp (4R-Py) parameters, strong evidence has been observed for the electrophilic character of the reactive species. The presence of the [(PBI)2(4R-Py)FeIIIOIPh]3+ intermediates and the effect of the co-ligands was also supported by UV-visible measurements, including the color change from red to green and the hypsochromic shifts in the presence of co-ligands. This is another indication that the title iron(III)-iodosylbenzene adduct is able to oxygenate sulfides and alkenes before it is transformed into the oxoiron form by cleavage of the O-I bond.

Regio- and stereoselective steroid hydroxylation by CYP109A2 from Bacillus megaterium explored by X-ray crystallography and computational modeling

The P450 monooxygenase CYP109A2 from Bacillus megaterium DSM319 was previously found to convert vitamin D3 (VD3) to 25-hydroxyvitamin D3. Here, we show that this enzyme is also able to convert testosterone in a highly regio- and stereoselective manner to 16β-hydroxytestosterone. To reveal the structural determinants governing the regio- and stereoselective steroid hydroxylation reactions catalyzed by CYP109A2, two crystal structures of CYP109A2 were solved in similar closed conformations, one revealing a bound testosterone in the active site pocket, albeit at a nonproductive site away from the heme-iron. To examine whether the closed crystal structures nevertheless correspond to a reactive conformation of CYP109A2, docking and molecular dynamics (MD) simulations were performed with testosterone and vitamin D3 (VD3) present in the active site. These MD simulations were analyzed for catalytically productive conformations, the relative occurrences of which were in agreement with the experimentally determined stereoselectivities if the predicted stability of each carbon-hydrogen bond was taken into account. Overall, the first-time determination and analysis of the catalytically relevant 3D conformation of CYP109A2 will allow for future small molecule ligand screening in silico, as well as enabling site-directed mutagenesis toward improved enzymatic properties of this enzyme.

QM Calculations Revealed that Outer-Sphere Electron Transfer Boosted O-O Bond Cleavage in the Multiheme-Dependent Cytochrome bd Oxygen Reductase

The cytochrome bd oxygen reductase catalyzes the four-electron reduction of dioxygen to two water molecules. The structure of this enzyme reveals three heme molecules in the active site, which differs from that of heme-copper cytochrome c oxidase. The quantum chemical cluster approach was used to uncover the reaction mechanism of this intriguing metalloenzyme. The calculations suggested that a proton-coupled electron transfer reduction occurs first to generate a ferrous heme b₅₉₅. This is followed by the dioxygen binding at the heme d center coupled with an outer-sphere electron transfer from the ferrous heme b₅₉₅ to the dioxygen moiety, affording a ferric ion superoxide intermediate. A second proton-coupled electron transfer produces a heme d ferric hydroperoxide, which undergoes efficient O-O bond cleavage facilitated by an outer-sphere electron transfer from the ferrous heme b₅₉₅ to the O-O σ* orbital and an inner-sphere proton transfer from the heme d hydroxyl group to the leaving hydroxide. The synergistic benefits of the two types of hemes rationalize the highly efficient oxygen reduction repertoire for the multi-heme-dependent cytochrome bd oxygen reductase family.

Computation-Aided Engineering of Cytochrome P450 for the Production of Pravastatin

CYP105AS1 is a cytochrome P450 from Amycolatopsis orientalis that catalyzes monooxygenation of compactin to 6-epi-pravastatin. For fermentative production of the cholesterol-lowering drug pravastatin, the stereoselectivity of the enzyme needs to be inverted, which has been partially achieved by error-prone PCR mutagenesis and screening. In the current study, we report further optimization of the stereoselectivity by a computationally aided approach. Using the CoupledMoves protocol of Rosetta, a virtual library of mutants was designed to bind compactin in a pro-pravastatin orientation. By examining the frequency of occurrence of beneficial substitutions and rational inspection of their interactions, a small set of eight mutants was predicted to show the desired selectivity and these variants were tested experimentally. The best CYP105AS1 variant gave >99% stereoselective hydroxylation of compactin to pravastatin, with complete elimination of the unwanted 6-epi-pravastatin diastereomer. The enzyme-substrate complexes were also examined by ultrashort molecular dynamics simulations of 50 × 100 ps and 5 × 22 ns, which revealed that the frequency of occurrence of near-attack conformations agreed with the experimentally observed stereoselectivity. These results show that a combination of computational methods and rational inspection could improve CYP105AS1 stereoselectivity beyond what was obtained by directed evolution. Moreover, the work lays out a general in silico framework for specificity engineering of enzymes of known structure.

Mechanism of Melatonin Metabolism by CYP1A1: What Determines the Bifurcation Pathways of Hydroxylation versus Deformylation?

Melatonin, a widely applied cosmetic active ingredient, has a variety of uses as a skin protector through antioxidant and anti-inflammatory functions as well as giving the body UV-induced defenses and immune system support. In the body, melatonin is synthesized from a tryptophan amino acid in a cascade of reactions, but as melatonin is toxic at high concentrations, it is metabolized in the human skin by the cytochrome P450 enzymes. The P450s are diverse heme-based mono-oxygenases that catalyze oxygen atom-transfer processes that trigger metabolism and detoxification reactions in the body. In the catalytic cycle of the P450s, a short-lived high-valent iron(IV)-oxo heme cation radical is formed that has been proposed to be the active oxidant. How and why it activates melatonin in the human body and what the origin of the product distributions is, are unknown. This encouraged us to do a detailed computational study on a typical human P450 isozyme, namely CYP1A1. We initially did a series of molecular dynamics simulations with substrate docked into several orientations. These simulations reveal a number of stable substrate-bound positions in the active site, which may lead to differences in substrate activation channels. Using tunneling analysis on the full protein structures, we show that two of the four binding conformations lead to open substrate-binding pockets. As a result, in these open pockets, the substrate is not tightly bound and can escape back into the solution. In the closed conformations, in contrast, the substrate is mainly oriented with the methoxy group pointing toward the heme, although under a different angle. We then created large quantum cluster models of the enzyme and focused on the chemical reaction mechanisms for melatonin activation, leading to competitive O-demethylation and C⁶-aromatic hydroxylation pathways. The calculations show that active site positioning determines the product distributions, but the bond that is activated is not necessarily closest to the heme in the enzyme-substrate complex. As such, the docking and molecular dynamics positioning of the substrate versus oxidant can give misleading predictions on product distributions. In particular, in quantum mechanics cluster model I, we observe that through a tight hydrogen bonding network, a preferential 6-hydroxylation of melatonin is obtained. However, O-demethylation becomes possible in alternative substrate-binding orientations that have the C⁶-aromatic ring position shielded. Finally, we investigated enzymatic and non-enzymatic O-demethylation processes and show that the hydrogen bonding network in the substrate-binding pocket can assist and perform this step prior to product release from the enzyme.

C-N Bond Forming Radical Rebound Is the Enantioselectivity-Determining Step in P411-Catalyzed Enantioselective C(sp3)-H Amination: A Combined Computational and Experimental Investigation

Engineered metalloenzymes represent promising catalysts for stereoselective C-H functionalization reactions. Recently, P450 enzymes have been evolved to allow for new-to-nature intramolecular C(sp3)-H amination reactions via a nitrene transfer mechanism, giving rise to diamine derivatives with excellent enantiocontrol. To shed light on the origin of enantioselectivity, a combined computational and experimental study was carried out. Hybrid quantum mechanics/molecular mechanics calculations were performed to investigate the activation energies and enantioselectivities of both the hydrogen atom transfer (HAT) and the subsequent C-N bond forming radical rebound steps. Contrary to previously hypothesized enantioinduction mechanisms, our calculations show that the radical rebound step is enantioselectivity-determining, whereas the preceding HAT step is only moderately stereoselective. Furthermore, the selectivity in the initial HAT is ablated by rapid conformational change of the radical intermediate prior to C-N bond formation. This finding is corroborated by our experimental study using a set of enantiomerically pure, monodeuterated substrates. Furthermore, classical and ab initio molecular dynamics simulations were carried out to investigate the conformational flexibility of the carbon-centered radical intermediate. This key radical species undergoes a facile conformational change in the enzyme active site from the pro-(R) to the pro-(S) configuration, whereas the radical rebound is slower due to the spin-state change and ring strain of the cyclization process, thereby allowing stereoablative C-N bond formation. Together, these studies revealed an underappreciated enantioinduction mechanism in biocatalytic C(sp3)-H functionalizations involving radical intermediates, opening up new avenues for the development of other challenging asymmetric C(sp3)-H functionalizations.

In Silico simulation of Cytochrome P450-Mediated metabolism of aromatic amines: A case study of N-Hydroxylation

Aromatic amines, the widely used raw materials in industry, cause long-term exposure to human bodies. They can be metabolized by cytochrome P450 enzymes to form active electrophilic compounds, which will potentially react with nucleophilic DNA to exert carcinogenesis. The short lifetime and versatility of the oxidant (a high-valent iron (IV)-oxo species, compound I) of P450 enzymes prompts us to use theoretical methods to investigate the metabolism of aromatic amines. In this work, the density functional theory (DFT) has been employed to simulate the hydroxylation metabolism through H-abstraction and to calculate the activation energy of this reaction for 28 aromatic amines. The results indicate that the steric effects, inductive effects and conjugative effects greatly contribute to the metabolism activity of the chemicals. The further correlation reveals that the dissociation energy of -NH₂ (BDE_N-H) can successfully predict the time-consuming calculated activation energy (R² for aromatic and heteroaromatic amines are 0.93 and 0.86, respectively), so BDE_N-H can be taken as a key parameter to characterize the relative stability of aromatic amines in P450 enzymes and further to quickly assess their potential toxicity. The validation results prove such relationship has good statistical performance (q_cv² for aromatic and heteroaromatic amines are 0.95 and 0.90, respectively) and can be used to other aromatic amines in the application domain, greatly reducing computational cost and providing useful support for experimental research.

The Many Chemists Who Could Have Proposed the Woodward-Hoffmann Rules (Including Roald Hoffmann) But Didn't: The Theoretical and Physical Chemists†

It is a reasonable question to ask, why, as of 1965 when the five Woodward-Hoffmann communication appeared, did no other physical chemist or chemical physicist or theoretical chemist discover the orbital symmetry rules for all pericyclic reactions? Two theoretical chemists - Luitzen Oosterhoff (in 1961) and Kenichi Fukui (in 1964) had discovered portions of the orbital symmetry rules; their stories appear in the papers immediately preceding this paper which is Paper 5 in a 27-paper series on the history of Woodward-Hoffmann rules. Concise yet telling stories of 19 other chemists who could have, might have, perhaps even should have discovered the Woodward-Hoffmann rules are presented with explanations as to why they did not do so. Social, political, and scientific explanations will summarize the analyses.

A framework for constructing linear free energy relationships to design molecular transition metal catalysts

A computational framework for ligand-driven design of transition metal complexes is presented in this work. We propose a general procedure for the construction of active site-specific linear free energy relationships (LFERs), which are inspired from Hammett and Taft correlations in organic chemistry and grounded in the activation strain model (ASM). Ligand effects are isolated and quantified in terms of their contribution to interaction and strain energy components of ASM. Scalar descriptors that are easily obtainable are then employed to construct the complete LFER. We successfully demonstrate proof-of-concept by constructing and applying an LFER to CH activation with enzyme-inspired [Cu2O2]2+ complexes. The key benefit of using ASM is a built-in compensation or error cancellation between LFER prediction of interaction and strain terms, resulting in accurate barrier predictions for 37 of the 47 catalysts examined in this study. The LFER is also transferable with respect to level of theory and flexible towards the choice of reference system. The absence of interaction-strain compensation or poor model performance for the remaining systems is a consequence of the approximate nature of the chosen interaction energy descriptor and LFER construction of the strain term, which focuses largely on trends in substrate and not catalyst strain.

CYP154C5 Regioselectivity in Steroid Hydroxylation Explored by Substrate Modifications and Protein Engineering*

CYP154C5 from Nocardia farcinica is a P450 monooxygenase able to hydroxylate a range of steroids with high regio- and stereoselectivity at the 16α-position. Using protein engineering and substrate modifications based on the crystal structure of CYP154C5, an altered regioselectivity of the enzyme in steroid hydroxylation had been achieved. Thus, conversion of progesterone by mutant CYP154C5 F92A resulted in formation of the corresponding 21-hydroxylated product 11-deoxycorticosterone in addition to 16α-hydroxylation. Using MD simulation, this altered regioselectivity appeared to result from an alternative binding mode of the steroid in the active site of mutant F92A. MD simulation further suggested that the entrance of water to the active site caused higher uncoupling in this mutant. Moreover, exclusive 15α-hydroxylation was observed for wild-type CYP154C5 in the conversion of 5α-androstan-3-one, lacking an oxy-functional group at C17. Overall, our data give valuable insight into the structure-function relationship of this cytochrome P450 monooxygenase for steroid hydroxylation.

Reactivity patterns of vanadium(iv/v)-oxo complexes with olefins in the presence of peroxides: a computational study

Vanadium porphyrin complexes are naturally occurring substances found in crude oil and have been shown to have medicinal properties as well. Little is known on their activities with substrates; therefore, we decided to perform a detailed density functional theory study on the properties and reactivities of vanadium(iv)- and vanadium(v)-oxo complexes with a TPPCl8 or 2,3,7,8,12,13,17,18-octachloro-meso-tetraphenylporphyrinato ligand system. In particular, we investigated the reactivity of [VV(O)(TPPCl8)]+ and [VIV(O)(TPPCl8)] with cyclohexene in the presence of H2O2 or HCO4-. The work shows that vanadium(iv)-oxo and vanadium(v)-oxo are sluggish oxidants by themselves and react with olefins slowly. However, in the presence of hydrogen peroxide, these metal-oxo species can be transformed into a side-on vanadium-peroxo complex, which reacts with substrates more efficiently. Particularly with anionic axial ligands, the side-on vanadium-peroxo and vanadium-oxo complexes produced epoxides from cyclohexene via small barrier heights. In addition to olefin epoxidation, we investigated aliphatic hydroxylation mechanisms by the same oxidants and some oxidants show efficient and viable cyclohexene hydroxylation mechanisms. The work implies that vanadium-oxo and vanadium-peroxo complexes can react with double bonds through epoxidation, and under certain conditions also undergo hydroxylation, but the overall reactivity is highly dependent on the equatorial ligand, the local environment and the presence or absence of anionic axial ligands.

Synthetic Fe/Cu Complexes: Toward Understanding Heme-Copper Oxidase Structure and Function

Heme-copper oxidases (HCOs) are terminal enzymes on the mitochondrial or bacterial respiratory electron transport chain, which utilize a unique heterobinuclear active site to catalyze the 4H+/4e- reduction of dioxygen to water. This process involves a proton-coupled electron transfer (PCET) from a tyrosine (phenolic) residue and additional redox events coupled to transmembrane proton pumping and ATP synthesis. Given that HCOs are large, complex, membrane-bound enzymes, bioinspired synthetic model chemistry is a promising approach to better understand heme-Cu-mediated dioxygen reduction, including the details of proton and electron movements. This review encompasses important aspects of heme-O2 and copper-O2 (bio)chemistries as they relate to the design and interpretation of small molecule model systems and provides perspectives from fundamental coordination chemistry, which can be applied to the understanding of HCO activity. We focus on recent advancements from studies of heme-Cu models, evaluating experimental and computational results, which highlight important fundamental structure-function relationships. Finally, we provide an outlook for future potential contributions from synthetic inorganic chemistry and discuss their implications with relevance to biological O2-reduction.

What drives the H-abstraction reaction in bio-mimetic oxoiron-bTAML complexes? A computational investigation

Monomeric iron-oxo units have been confirmed as intermediates involved in the C-H bond activation in various metallo-enzymes. Biomimetic oxoiron complexes of the biuret modified tetra-amido macrocyclic ligand (bTAML) have been demonstrated to oxidize a wide variety of unactivated C-H bonds. In the current work, density functional theory (DFT) has been employed to investigate the hydrogen abstraction (HAT) reactivity differences across a series of bTAML complexes. The cause for the differences in the HAT energy barriers has been found to be the relative changes in the energy of the frontier molecular orbitals (FMOs) induced by electronic perturbation.

How Do Ring Size and π-Donating Thiolate Ligands Affect Redox-Active, α-Imino-N-heterocycle Ligand Activation?

Considerable effort has been devoted to the development of first-row transition-metal catalysts containing redox-active imino-pyridine ligands that are capable of storing multiple reducing equivalents. This property allows abundant and inexpensive first-row transition metals, which favor sequential one-electron redox processes, to function as competent catalysts in the concerted two-electron reduction of substrates. Herein we report the syntheses and characterization of a series of iron complexes that contain both π-donating thiolate and π-accepting (α-imino)-N-heterocycle redox-active ligands, with progressively larger N-heterocycle rings (imidazole, pyridine, and quinoline). A cooperative interaction between these complementary redox-active ligands is shown to dictate the properties of these complexes. Unusually intense charge-transfer (CT) bands, and intraligand metrical parameters, reminiscent of a reduced (α-imino)-N-heterocycle ligand (L^•-), initially suggested that the electron-donating thiolate had reduced the N-heterocycle. Sulfur K-edge X-ray absorption spectroscopic (XAS) data, however, provides evidence for direct communication, via backbonding, between the thiolate sulfur and the formally orthogonal (α-imino)-N-heterocycle ligand π*-orbitals. DFT calculations provide evidence for extensive delocalization of bonds over the sulfur, iron, and (α-imino)-N-heterocycle, and TD-DFT shows that the intense optical CT bands involve transitions between a mixed Fe/S donor, and (α-imino)-N-heterocycle π*-acceptor orbital. The energies and intensities of the optical and S K-edge pre-edge XAS transitions are shown to correlate with N-heterocycle ring size, as do the redox potentials. When the thiolate is replaced with a thioether, or when the low-spin S = 0 Fe(II) is replaced with a high-spin S = 3/2 Co(II), the N-heterocycle ligand metrical parameters and electronic structure do not change relative to the neutral L⁰ ligand. With respect to the development of future catalysts containing redox-active ligands, the energy cost of storing reducing equivalents is shown to be lowest when a quinoline, as opposed to imidazole or pyridine, is incorporated into the ligand backbone of the corresponding Fe complex.

Influence of Ligand Architecture in Tuning Reaction Bifurcation Pathways for Chlorite Oxidation by Non-Heme Iron Complexes

Reaction bifurcation processes are often encountered in the oxidation of substrates by enzymes and generally lead to a mixture of products. One particular bifurcation process that is common in biology relates to electron transfer versus oxygen atom transfer by high-valent iron(IV)-oxo complexes, which nature uses for the oxidation of metabolites and drugs. In biomimicry and bioremediation, an important reaction relates to the detoxification of ClO_x^- in water, which can lead to a mixture of products through bifurcated reactions. Herein we report the first three water-soluble non-heme iron(II) complexes that can generate chlorine dioxide from chlorite at ambient temperature and physiological pH. These complexes are highly active oxygenation oxidants and convert ClO₂^- into either ClO₂ or ClO₃¯ via high-valent iron(IV)-oxo intermediates. We characterize the short-lived iron(IV)-oxo species and establish rate constants for the bifurcation mechanism leading to ClO₂ and ClO₃^- products. We show that the ligand architecture of the metal center plays a dominant role by lowering the reduction potential of the metal center. Our experiments are supported by computational modeling, and a predictive valence bond model highlights the various factors relating to the substrate and oxidant that determine the bifurcation pathway and explains the origins of the product distributions. Our combined kinetic, spectroscopic, and computational studies reveal the key components necessary for the future development of efficient chlorite oxidation catalysts.

A mononuclear non-heme-iron dioxygen-carrying protein?

The ability of mononuclear non-heme iron complexes to function as molecular oxygen transporters is investigated by density functional theory. The factors governing the efficiency of the reversible binding of dioxygen at the active site of the dinuclear non-heme iron enzyme hemerythrin, including antiferromagnetic coupling and the conversion of dioxygen to hydroperoxo by a proton coupled 2-electron transfer mechanism, are revisited and considered as possible tools in mononuclear non-heme complexes. Several mononuclear non-heme model complexes, including active sites of enzymes already known to interact with dioxgenic ligands, are constructed and the molecular oxygen transportation capabilities of these complexes are examined computationally. The high-spin nature of the ground state of these complexes implies an intrinsic kinetic lability of the oxy structures, as also evident from potential energy surface calculations towards iron-dioxygen cleavage. Proton affinities as calibrated with reference compounds showed that these complexes are highly unlikely to undergo protonation to form hydroperoxo-like adducts. Mixed superoxo descriptions of the dissociated dioxygenic ligands in all complexes add to the overall conclusion that these model structures are significantly disadvantaged in any attempt to be employed for molecular oxygen transportation.

Computation Sheds Insight into Iron Porphyrin Carbenes' Electronic Structure, Formation, and N-H Insertion Reactivity

Iron porphyrin carbenes constitute a new frontier of species with considerable synthetic potential. Exquisitely engineered myoglobin and cytochrome P450 enzymes can generate these complexes and facilitate the transformations they mediate. The current work harnesses density functional theoretical methods to provide insight into the electronic structure, formation, and N-H insertion reactivity of an iron porphyrin carbene, [Fe(Por)(SCH3)(CHCO2Et)](-), a model of a complex believed to exist in an experimentally studied artificial metalloenzyme. The ground state electronic structure of the terminal form of this complex is an open-shell singlet, with two antiferromagnetically coupled electrons residing on the iron center and carbene ligand. As we shall reveal, the bonding properties of [Fe(Por)(SCH3)(CHCO2Et)](-) are remarkably analogous to those of ferric heme superoxide complexes. The carbene forms by dinitrogen loss from ethyl diazoacetate. This reaction occurs preferentially through an open-shell singlet transition state: iron donates electron density to weaken the C-N bond undergoing cleavage. Once formed, the iron porphyrin carbene accomplishes N-H insertion via nucleophilic attack. The resulting ylide then rearranges, using an internal carbonyl base, to form an enol that leads to the product. The findings rationalize experimentally observed reactivity trends reported in artificial metalloenzymes employing iron porphyrin carbenes. Furthermore, these results suggest a possible expansion of enzymatic substrate scope, to include aliphatic amines. Thus, this work, among the first several computational explorations of these species, contributes insights and predictions to the surging interest in iron porphyrin carbenes and their synthetic potential.

Tuning reactivity and selectivity in hydrogen atom transfer from aliphatic C-H bonds to alkoxyl radicals: role of structural and medium effects

Hydrogen atom transfer (HAT) is a fundamental reaction that takes part in a wide variety of chemical and biological processes, with relevant examples that include the action of antioxidants, damage to biomolecules and polymers, and enzymatic and biomimetic reactions. Moreover, great attention is currently devoted to the selective functionalization of unactivated aliphatic C-H bonds, where HAT based procedures have been shown to play an important role. In this Account, we describe the results of our recent studies on the role of structural and medium effects on HAT from aliphatic C-H bonds to the cumyloxyl radical (CumO(•)). Quantitative information on the reactivity and selectivity patterns observed in these reactions has been obtained by time-resolved kinetic studies, providing a deeper understanding of the factors that govern HAT from carbon and leading to the definition of useful guidelines for the activation or deactivation of aliphatic C-H bonds toward HAT. In keeping with the electrophilic character of alkoxyl radicals, polar effects can play an important role in the reactions of CumO(•). Electron-rich C-H bonds are activated whereas those that are α to electron withdrawing groups are deactivated toward HAT, with these effects being able to override the thermodynamic preference for HAT from the weakest C-H bond. Stereoelectronic effects can also influence the reactivity of the C-H bonds of ethers, amines, and amides. HAT is most rapid when these bonds can be eclipsed with a lone pair on an adjacent heteroatom or with the π-system of an amide functionality, thus allowing for optimal orbital overlap. In HAT from cyclohexane derivatives, tertiary axial C-H bond deactivation and tertiary equatorial C-H bond activation have been observed. These effects have been explained on the basis of an increase in torsional strain or a release in 1,3-diaxial strain in the HAT transition states, with kH(eq)/kH(ax) ratios that have been shown to exceed one order of magnitude. Medium effects on HAT from aliphatic C-H bonds to CumO(•) have been also investigated. With basic substrates, from large to very large decreases in kH have been measured with increasing solvent hydrogen bond donor (HBD) ability or after addition of protic acids or alkali and alkaline earth metal ions, with kinetic effects that exceed 2 orders of magnitude in the reactions of tertiary alkylamines and alkanamides. Solvent hydrogen bonding, protonation, and metal ion binding increase the electron deficiency and the strength of the C-H bonds of these substrates deactivating these bonds toward HAT, with the extent of this deactivation being modulated by varying the nature of the substrate, solvent, protic acid, and metal ion. These results indicate that through these interactions careful control over the HAT reactivity of basic substrates toward CumO(•) and other electrophilic radicals can be achieved, suggesting moreover that these effects can be exploited in an orthogonal fashion for selective C-H bond functionalization of substrates bearing different basic functionalities.

Nonheme Fe(IV) Oxo Complexes of Two New Pentadentate Ligands and Their Hydrogen-Atom and Oxygen-Atom Transfer Reactions

Two new pentadentate {N5} donor ligands based on the N4Py (N4Py = N,N-bis(2-pyridylmethyl)-N-bis(2-pyridyl)methylamine) framework have been synthesized, viz. [N-(1-methyl-2-benzimidazolyl)methyl-N-(2-pyridyl)methyl-N-(bis-2-pyridyl methyl)amine] (L(1)) and [N-bis(1-methyl-2-benzimidazolyl)methyl-N-(bis-2-pyridylmethyl)amine] (L(2)), where one or two pyridyl arms of N4Py have been replaced by corresponding (N-methyl)benzimidazolyl-containing arms. The complexes [Fe(II)(CH3CN)(L)](2+) (L = L(1) (1); L(2) (2)) were synthesized, and reaction of these ferrous complexes with iodosylbenzene led to the formation of the ferryl complexes [Fe(IV)(O)(L)](2+) (L = L(1) (3); L(2) (4)), which were characterized by UV-vis spectroscopy, high resolution mass spectrometry, and Mössbauer spectroscopy. Complexes 3 and 4 are relatively stable with half-lives at room temperature of 40 h (L = L(1)) and 2.5 h (L = L(2)). The redox potentials of 1 and 2, as well as the visible spectra of 3 and 4, indicate that the ligand field weakens as ligand pyridyl substituents are progressively substituted by (N-methyl)benzimidazolyl moieties. The reactivities of 3 and 4 in hydrogen-atom transfer (HAT) and oxygen-atom transfer (OAT) reactions show that both complexes exhibit enhanced reactivities when compared to the analogous N4Py complex ([Fe(IV)(O)(N4Py)](2+)), and that the normalized HAT rates increase by approximately 1 order of magnitude for each replacement of a pyridyl moiety; i.e., [Fe(IV)(O)(L(2))](2+) exhibits the highest rates. The second-order HAT rate constants can be directly related to the substrate C-H bond dissociation energies. Computational modeling of the HAT reactions indicates that the reaction proceeds via a high spin transition state.

Magnetic circular dichroism and computational study of mononuclear and dinuclear iron(IV) complexes

High-valent iron(IV)-oxo species are key intermediates in the catalytic cycles of a range of O2-activating iron enzymes. This work presents a detailed study of the electronic structures of mononuclear ([FeIV(O)(L)(NCMe)]2+, 1, L = tris(3,5-dimethyl-4-methoxylpyridyl-2-methyl)amine) and dinuclear ([(L)FeIV(O)(μ-O)FeIV(OH)(L)]3+, 2) iron(IV) complexes using absorption (ABS), magnetic circular dichroism (MCD) spectroscopy and wave-function-based quantum chemical calculations. For complex 1, the experimental MCD spectra at 2-10 K are dominated by a broad positive C-term band between 12000 and 18000 cm-1. As the temperature increases up to ~20 K, this feature is gradually replaced by a derivative-shaped signal. The computed MCD spectra are in excellent agreement with experiment, which reproduce not only the excitation energies and the MCD signs of key transitions but also their temperature-dependent intensity variations. To further corroborate the assignments suggested by the calculations, the individual MCD sign for each transition is independently determined from the corresponding electron donating and accepting orbitals. Thus, unambiguous assignments can be made for the observed transitions in 1. The ABS/MCD data of complex 2 exhibit ten features that are assigned as ligand-field transitions or oxo- or hydroxo-to-metal charge transfer bands, based on MCD/ABS intensity ratios, calculated excitation energies, polarizations, and MCD signs. In comparison with complex 1, the electronic structure of the FeIV=O site is not significantly perturbed by the binding to another iron(IV) center. This may explain the experimental finding that complexes 1 and 2 have similar reactivities toward C-H bond activation and O-atom transfer.