Does substrate oxidation determine the regioselectivity of cyclohexene and propene oxidation by cytochrome p450?

DFT and QM/MM computations of allylic C-H hydroxylation versus C=C epoxidation in cyclohexene and propene by Compound I of P450cam demonstrate that the relative barriers for the oxidative processes themselves are not good predictors of the observed selectivity. However, a kinetic expression previously developed (Kozuch, S.; Shaik, S. J. Am. Chem. Soc. 2006, 128, 3355) for catalytic cycles under steady-state conditions, predicts, in accord with experiment, that propene will undergo exclusive C=C epoxidation, while cyclohexene will undergo both reactions with a small preference for epoxidation. The model expression for the effective barrier of the cycle forms a general basis for understanding and predicting the selectivity of P450 isozymes.     

NR transfer reactivity of azo-compound I of P450. How does the nitrogen substituent tune the reactivity of the species toward C-H and C=C activation?

We studied electronic structures and reactivity patterns of azo-compound I species (RN-Cpd I) by comparison to O-Cpd I of, e.g., cytochrome P450. The study shows that the RN-Cpd I species are capable of C=C aziridination and C-H amidation, in a two-state mechanism similar to that of O-Cpd I. However, unlike O-Cpd I, here the nitrogen substituent (R) exerts a major impact on structure and reactivity. Thus, it is demonstrated that Fe=NR bonds of RN-Cpd I will generally be substantially longer than Fe=O bonds; electron-withdrawing R groups will generate a very long Fe=N bond, whereas electron-releasing R groups should have the opposite effect and hence a shorter Fe=N bond. The R substituent controls also the reactivity of RN-Cpd I toward C=C and C-H bonds by exerting steric and electronic effects. Our analysis shows that an electron-releasing substituent will lower the barriers for both bond activation reactions, since the electronic factor makes the reactions highly exothermic, while an electron-withdrawing one should raise both barriers. The steric bulk of the substituent is predicted to inhibit more strongly the aziridination reactions. It is predicted that electron-releasing substituents with small bulk will create powerful aziridination reagents, whereas electron-withdrawing substituents like MeSO(2) will prefer C-H bond activation with preference that increases with steric bulk. Finally, the study predicts (i) that the reactions of RN-Cpd I will be less stereospecific than those of O-Cpd I and (ii) that aziridination will be more stereoselective than amidation.     

A valence bond modeling of trends in hydrogen abstraction barriers and transition states of hydroxylation reactions catalyzed by cytochrome P450 enzymes

The paper outlines the fundamental factors that govern the mechanisms of alkane hydroxylation by cytochrome P450 and the corresponding barrier heights during the hydrogen abstraction and radical rebound steps of the process. This is done by a combination of density functional theory calculations for 11 alkanes and valence bond (VB) modeling of the results. The energy profiles and transition states for the various steps are reconstructed using VB diagrams (Shaik, S. S. J. Am. Chem. Soc. 1981, 103, 3692-3701. Shaik, S.; Shurki, A. Angew. Chem. Int. Ed. 1999, 38, 586-625.) and the DFT barriers are reproduced by the VB model from raw data based on C-H bond energies. The model explains a variety of other features of P450 hydroxylations: (a) the nature of the polar effect during hydrogen abstraction, (b) the difference between the activation mechanisms leading to the Fe(IV) vs the Fe(III) electromers, (c) the difference between the gas phase and the enzymatic reaction, and (d) the dependence of the rebound barrier on the spin state. The VB mechanism shows that the active species of the enzyme performs a complex reaction that involves multiple bond making and breakage mechanisms by utilizing an intermediate VB structure that cuts through the high barrier of the principal transformation between reactants and products, thereby mediating the process at a low energy cost. The correlations derived in this paper create order and organize the data for a process of a complex and important enzyme. This treatment can be generalized to the reactivity patterns of nonheme systems and synthetic iron-oxo porphyrin reagents.     

Trapping-mediated dissociative chemisorption of cycloalkanes on Ru(001) and Ir(111): influence of ring strain and molecular geometry on the activation of C-C and C-H bonds.

We have measured the initial probabilities of dissociative chemisorption of perhydrido and perdeutero cycloalkane isotopomers on the hexagonally close-packed Ru(001) and Ir(111) single-crystalline surfaces for surface temperatures between 250 and 1100 K. Kinetic parameters (activation barrier and preexponential factor) describing the initial, rate-limiting C-H or C-C bond cleavage reactions were quantified for each cycloalkane isotopomer on each surface. Determination of the dominant initial reaction mechanism as either initial C-C or C-H bond cleavage was judged by the presence or absence of a kinetic isotope effect between the activation barriers for each cycloalkane isotopomer pair, and also by comparison with other relevant alkane activation barriers. On the Ir(111) surface, the dissociative chemisorption of cyclobutane, cyclopentane, and cyclohexane occurs via two different reaction pathways: initial C-C bond cleavage dominates on Ir(111) at high temperature (T > approximately 600 K), while at low temperature (T < approximately 400 K), initial C-H bond cleavage dominates. On the Ru(001) surface, dissociative chemisorption of cyclopentane occurs via initial C-C bond cleavage over the entire temperature range studied, whereas dissociative chemisorption of both cyclohexane and cyclooctane occurs via initial C-H bond cleavage. Comparison of the cycloalkane C-C bond activation barriers measured here with those reported previously in the literature qualitatively suggests that the difference in ring-strain energies between the initial state and the transition state for ring-opening C-C bond cleavage effectively lowers or raises the activation barrier for dissociative chemisorption via C-C bond cleavage, depending on whether the transition state is less or more strained than the initial state. Moreover, steric arguments and metal-carbon bond strength arguments have been evoked to explain the observed trend of decreasing C-H bond activation barrier with decreasing cycloalkane ring size.     

Decomposition of methylamine on Mo(100) surface: A DFT study

The initial decomposition of methylamine on Mo(100) surface has been investigated by self-consistent (GGA-PW91) density functional theory combined with periodic slab model. The adsorption energies of possible species and the activation energies for possible elementary reactions involved are obtained in the present work. Our results indicate that the barriers decreased with the order of C-N>N-H>C-H. In addition, metastable adsorption of the abstracted hydrogen atom on the hollow site in the final state is also considered for the N-H and C-H bond breaking. For the C-H bond cleavage, the reaction barrier that the abstracted hydrogen located on the hollow site in the final state is lower than that on the bridge site. However, for the N H bond breaking, the barriers are alike for the abstracted hydrogen on both hollow and bridge sites in the final state.     

CX_2(X=H,F,Cl)与甲基异丙基醚C-H键插入反应的理论研究

用从头计算方法在MP2 /6 31G(d)水平上研究了CX2 (X =H ,F ,Cl)与甲基异丙基醚的C -H键插入反应。CCl2 与甲基异丙基醚两个不同的α C的C -H键插入势垒分别为 117.2kJ/mol (甲基 )和 2 0 .6kJ/mol (异丙基 )。CF2 与异丙基α C的C -H键上插入势垒为 12 0 .0kJ/mol,在插入甲基上C -H键时会引起C -O键的断裂。CH2 的插入反应则不需要势垒。对CX2 与二甲醚、甲乙醚、甲基异丙基醚、甲基苄基醚上各种不同的C -H键插入势垒进行了比较 ,甲基和苯基都促使其毗邻的C -H键更容易被CX2 所插入     

P450 fingerprinting method for rapid discovery of terpene hydroxylating P450 catalysts with diversified regioselectivity

Engineered P450 enzymes constitute attractive catalysts for the selective oxidation of unactivated C-H bonds in complex molecules. A current bottleneck in the use of P450 catalysis for chemical synthesis is the time and effort required to identify the P450 variant(s) with the desired level of activity and selectivity. In this report, we describe a method to map the active site configuration of engineered P450 variants in high throughput using a set of semisynthetic chromogenic probes. Through analysis of the resulting 'fingerprints', reliable predictions can be made regarding the reactivity of these enzymes toward complex substrates structurally related to the fingerprint probes. In addition, fingerprint analysis offers a convenient and time-effective means to assess the regioselectivity properties of the fingerprinted P450s. The described approach can represent a valuable tool to expedite the discovery of P450 oxidation catalysts for the functionalization of relevant natural products such as members of the terpene family.     

Is the ruthenium analogue of compound I of cytochrome p450 an efficient oxidant? A theoretical investigation of the methane hydroxylation reaction

High-valent metal-oxo complexes catalyze C-H bond activation by oxygen insertion, with an efficiency that depends on the identity of the transition metal and its oxidation state. Our study uses density functional calculations and theoretical analysis to derive fundamental factors of catalytic activity, by comparison of a ruthenium-oxo catalyst with its iron-oxo analogue toward methane hydroxylation. The study focuses on the ruthenium analogue of the active species of the enzyme cytochrome P450, which is known to be among the most potent catalysts for C-H activation. The computed reaction pathways reveal one high-spin (HS) and two low-spin (LS) mechanisms, all nascent from the low-lying states of the ruthenium-oxo catalyst (Ogliaro, F.; de Visser, S. P.; Groves, J. T.; Shaik, S. Angew. Chem. Int. Ed. 2001, 40, 2874-2878). These mechanisms involve a bond activation phase, in which the transition states (TS's) appear as hydrogen abstraction species, followed by a C-O bond making phase, through a rebound of the methyl radical on the metal-hydroxo complex. However, while the HS mechanism has a significant rebound barrier, and hence a long lifetime of the radical intermediate, by contrast, the LS ones are effectively concerted with small barriers to rebound, if at all. Unlike the iron catalyst, the hydroxylation reaction for the ruthenium analogue is expected to follow largely a single-state reactivity on the LS surface, due to a very large rebound barrier of the HS process and to the more efficient spin crossover expected for ruthenium. As such, ruthenium-oxo catalysts (Groves, J. T.; Shalyaev, K.; Lee, J. In The Porphyrin Handbook; Biochemistry and Binding: Activation of Small Molecules, Vol. 4; Kadish, K. M., Smith, K. M., Guilard, R., Eds.; Academic Press: New York, 2000; pp 17-40) are expected to lead to more stereoselective hydroxylations compared with the corresponding iron-oxo reactions. It is reasoned that the ruthenium-oxo catalyst should have larger turnover numbers compared with the iron-oxo analogue, due to lesser production of suicidal side products that destroy the catalyst (Ortiz de Montellano, P. R.; Beilan, H. S.; Kunze, K. L.; Mico, B. A. J. Biol. Chem. 1981, 256, 4395-4399). The computations reveal also that the ruthenium complex is more electrophilic than its iron analogue, having lower hydrogen abstraction barriers. These reactivity features of the ruthenium-oxo system are analyzed and shown to originate from a key fundamental factor, namely, the strong 4d(Ru)-2p(O,N) overlaps, which produce high-lying pi(Ru-O), sigma(Ru-O), and sigma(Ru-N) orbitals and thereby to lead to a preference of ruthenium for higher-valent oxidation states with higher electrophilicity, for the effectively concerted LS hydroxylation mechanism, and for less suicidal complexes. As such, the ruthenium-oxo species is predicted to be a more robust catalyst than its iron-oxo analogue.     

Diversity of P450 enzymes in the biosynthesis of natural products

Covering: 1985 to 2012Diverse oxygenation patterns of natural products generated by secondary metabolic pathways in microorganisms and plants are largely achieved through the tailoring reactions catalysed by cytochrome P450 enzymes (P450s). P450s are a large family of oxidative hemoproteins found in all life forms from prokaryotes to humans. Understanding the reactivity and selectivity of these fascinating C-H bond-activating catalysts will advance their use in generating valuable pharmaceuticals and products for medicine, agriculture and industry. A major strength of this P450 group is its set of established enzyme-substrate relationships, the source of the most detailed knowledge on how P450 enzymes work. Engineering microbial-derived P450 enzymes to accommodate alternative substrates and add new functions continues to be an important near- and long-term practical goal driving the structural characterization of these molecules. Understanding the natural evolution of P450 structure-function should accelerate metabolic engineering and directed evolutionary approaches to enhance diversification of natural product structures and other biosynthetic applications.     

Theoretical study of N-demethylation of substituted N,N-dimethylanilines by cytochrome P450: the mechanistic significance of kinetic isotope effect profiles

The mechanism of N-demethylation of N,N-dimethylanilines (DMAs) by cytochrome P450, a highly debated topic in mechanistic bioinorganic chemistry (Karki, S. B.; Dinnocenczo, J. P.; Jones, J. P.; Korzekwa, K. R. J. Am. Chem. Soc. 1995, 117, 3657), is studied here using DFT calculations of the reactions of the active species of the enzyme, Compound I (Cpd I), with four para-(H, Cl, CN, NO2) substituted DMAs. The calculations resolve mechanistic controversies, offer a consistent mechanistic view, and reveal the following features: (a) the reaction pathways involve C-H hydroxylation by Cpd I followed by a nonenzymatic carbinolamine decomposition. (b) C-H hydroxylation is initiated by a hydrogen atom transfer (HAT) step that possesses a "polar" character. As such, the HAT energy barriers correlate with the energy level of the HOMO of the DMAs. (c) The series exhibits a switch from spin-selective reactivity for DMA and p-Cl-DMA to two-state reactivity, with low- and high-spin states, for p-CN-DMA and p-NO2-DMA. (d) The computed kinetic isotope effect profiles (KIEPs) for these scenarios match the experimentally determined KIEPs. Theory further shows that the KIEs and TS structures vary in a manner predicted by the Melander-Westheimer postulate: as the substituent becomes more electron withdrawing, the TS is shifted to a later position along the H-transfer coordinate and the corresponding KIEs increases. (e) The generated carbinolaniline can readily dissociate from the heme and decomposes in a nonenzymatic environment, which involves water assisted proton shift.     

Ab initio reaction path analysis of benzene hydrogenation to cyclohexane on Pt(111)

First-principles density functional theory calculations were performed to obtain detailed insight into the mechanism of benzene hydrogenation over Pt(111). The results indicate that benzene hydrogenation follows a Horiuti-Polanyi scheme which involves the consecutive addition of hydrogen adatoms. A first-principles-based reaction path analysis indicates the presence of a dominant reaction path. Hydrogenation occurs preferentially in the meta position of a methylene group. Cyclohexadiene and cyclohexene are expected to be at best minor products, since they are not formed along the dominant reaction path. The only product that can desorb is cyclohexane. Along the dominant reaction path, two categories of activation energies are found: lower barriers at approximately 75 kJ/mol for the first three hydrogenation steps, and higher barriers of approximately 88 kJ/mol for steps four and six, where hydrogen can only add in the ortho position of two methylene groups. The highest barrier at 104 kJ/mol is calculated for the fifth hydrogenation step, which may potentially be the rate-determining step. The high barrier for this step is likely the result of a rather strong C-H...Pt interaction in the adsorbed reactant state (1,2,3,5-tetrahydrobenzene) which increases the barrier by approximately 15 kJ/mol. Benzene and hydrogen are thought to be the most-abundant reaction intermediates.     

The barriers to methyl group rotation in s-cis- and s-trans-methyl nitrite

Ab initio Hartree-Fock calculations at the 4–21 level reproduce the structures of s-cis- and s-trans-methyl nitrite as well as the surprisingly large difference in the methyl group rotational barriers of the two conformers. Furthermore, the computations reveal the significant structural changes accompanying internal rotation. These structural changes, together with the localized orbitais and population analyses, provide information for an understanding of the unusual barriers. The extraordinarily low barrier of the s-trans conformer results from stabilization of the higher-energy eclipsed form by an attractive interaction involving the in-plane methyl hydrogen atom and the lone pair of the nitrogen atom. The high barrier of the s-cis form is due to additional stabilization of the staggered conformation by a hyperconjugative interaction of the out-of-plane C-H groups with the NO double bond.     

Effect of the axial ligand on substrate sulfoxidation mediated by iron(IV)-oxo porphyrin cation radical oxidants

Cytochromes P450 catalyze a range of different oxygen-transfer processes including aliphatic and aromatic hydroxylation, epoxidation, and sulfoxidation reactions. Herein, we have investigated substrate sulfoxidation mediated by models of P450 enzymes as well as by biomimetic oxidants using density functional-theory methods and we have rationalized the sulfoxidation reaction barriers and rate constants. We carried out two sets of calculations: first, we calculated the sulfoxidation by an iron(IV)-oxo porphyrin cation radical oxidant [Fe(IV)=O(Por(+.))SH] that mimics the active site of cytochrome P450 enzymes with a range of different substrates, and second, we studied one substrate (dimethyl sulfide) with a selection of different iron(IV)-oxo porphyrin cation radical oxidants [Fe(IV)=O(Por(+.))L] with varying axial ligands L. The study presented herein shows that the barrier height for substrate sulfoxidation correlates linearly with the ionization potential of the substrate, thus reflecting the electron-transfer processes in the rate-determining step of the reaction. Furthermore, the axial ligand of the oxidant influences the pK(a) value of the iron(IV)-oxo group, and, as a consequence, the bond dissociation energy (BDE(OH) value correlates with the barrier height for the reverse sulfoxidation reaction. These studies have generalized substrate-sulfoxidation reactions and have shown how they fundamentally compare with substrate hydroxylation and epoxidation reactions.     

The final catalytic step of cytochrome p450 aromatase: a density functional theory study

B3LYP density functional theory calculations are used to unravel the mysterious third step of aromatase catalysis. The feasibility of mechanisms in which the reduced ferrous dioxygen intermediate mediates androgen aromatization is explored and determined to be unlikely. However, proton-assisted homolysis of the peroxo hemiacetal intermediate to produce P450 compound I and the C19 gem-diol likely proceeds with a low energetic barrier. Mechanisms for the aromatization and deformylation sequence which are initiated by 1beta-hydrogen atom abstraction by P450 compound I are considered. 1beta-Hydrogen atom abstraction from substrates in the presence of the 2,3-enol encounters strikingly low barriers (5.3-7.8 kcal/mol), whereas barriers for this same process rise to 17.0-27.1 kcal/mol in the keto tautomer. Transition states for 1beta-hydrogen atom abstraction from enolized substrates in the presence of the 19-gem-diol decayed directly to the experimentally observed products. If the C19 aldehyde remains unhydrated, aromatization occurs with concomitant decarbonylation and therefore does not support dehydration of the C19 aldehyde prior to the final catalytic step. On the doublet surface, the transition state connects to a potentially labile 1(10) dehydrogenated product, which may undergo rapid aromatization, as well as formic acid. Ab initio molecular dynamics confirmed that the 1beta-hydrogen atom abstraction and deformylation or decarbonylation occur in a nonsynchronous, coordinated manner. These calculations support a dehydrogenase behavior of aromatase in the final catalytic step, which can be summarized by 1beta-hydrogen atom abstraction followed by gem-diol deprotonation.     

Can ferric-superoxide act as a potential oxidant in P450(cam)? QM/MM investigation of hydroxylation, epoxidation, and sulfoxidation

In view of recent reports of high reactivity of ferric-superoxide species in heme and nonheme systems (Morokuma et al. J. Am. Chem. Soc. 2010, 132, 11993-12005; Que et al. Inorg. Chem. 2010, 49, 3618-3628; Nam et al. J. Am. Chem. Soc. 2010, 132, 5958-5959; J. Am. Chem. Soc. 2010, 132, 10668-10670), we use herein combined quantum mechanics/molecular mechanics (QM/MM) methods to explore the potential reactivity of P450(cam) ferric-superoxide toward hydroxylation, epoxidation, and sulfoxidation. The calculations demonstrate that P450 ferric-superoxide is a sluggish oxidant compared with the high-valent oxoiron porphyrin cation-radical species. As such, unlike heme enzymes with a histidine axial ligand, the P450 superoxo species does not function as an oxidant in P450(cam). The origin of this different behavior of the superoxo species of P450 vis-a?-vis other heme enzymes like tryptophan 2, 3-dioxygenase (TDO) is traced to the ability of the latter superoxo species to make a stronger FeOO-X (X = H,C) bond and to stabilize the corresponding bond-activation transition states by resonance with charge-transfer configurations. By contrast, the negatively charged thiolate ligand in the P450 superoxo species minimizes the mixing of charge transfer configurations in the transition state and raises the reaction barrier. However, as we demonstrate, an external electric field oriented along the Fe-O axis with a direction pointing from Fe toward O will quench Cpd I formation by slowing the reduction of ferric-superoxide and will simultaneously lower the barriers for oxidation by the latter species, thereby enabling observation of superoxo chemistry in P450. Other options for nascent superoxo reactivity in P450 are discussed.     

Predictable stereoselective and chemoselective hydroxylations and epoxidations with P450 3A4

Enantioselective hydroxylation of one specific methylene in the presence of many similar groups is debatably the most challenging chemical transformation. Although chemists have recently made progress toward the hydroxylation of inactivated C-H bonds, enzymes such as P450s (CYPs) remain unsurpassed in specificity and scope. The substrate promiscuity of many P450s is desirable for synthetic applications; however, the inability to predict the products of these enzymatic reactions is impeding advancement. We demonstrate here the utility of a chemical auxiliary to control the selectivity of CYP3A4 reactions. When linked to substrates, inexpensive, achiral theobromine directs the reaction to produce hydroxylation or epoxidation at the fourth carbon from the auxiliary with pro-R facial selectivity. This strategy provides a versatile yet controllable system for regio-, chemo-, and stereoselective oxidations at inactivated C-H bonds and demonstrates the utility of chemical auxiliaries to mediate the activity of highly promiscuous enzymes.     

On the mechanism of (PCP)Ir-catalyzed acceptorless dehydrogenation of alkanes: a combined computational and experimental study

Pincer complexes of the type ((R)PCP)IrH(2), where ((R)PCP)Ir is [eta(3)-2,6-(R(2)PCH(2))(2)C(6)H(3)]Ir, are the most effective catalysts reported to date for the "acceptorless" dehydrogenation of alkanes to yield alkenes and free H(2). We calculate (DFT/B3LYP) that associative (A) reactions of ((Me)PCP)IrH(2) with model linear (propane, n-PrH) and cyclic (cyclohexane, CyH) alkanes may proceed via classical Ir(V) and nonclassical Ir(III)(eta(2)-H(2)) intermediates. A dissociative (D) pathway proceeds via initial loss of H(2), followed by C-H addition to ((Me)PCP)Ir. Although a slightly higher energy barrier (DeltaE(+ +)) is computed for the D pathway, the calculated free-energy barrier (DeltaG(+ +)) for the D pathway is significantly lower than that of the A pathway. Under standard thermodynamic conditions (STP), C-H addition via the D pathway has DeltaG(o)(+ +) = 36.3 kcal/mol for CyH (35.1 kcal/mol for n-PrH). However, acceptorless dehydrogenation of alkanes is thermodynamically impossible at STP. At conditions under which acceptorless dehydrogenation is thermodynamically possible (for example, T = 150 degrees C and P(H)2 = 1.0 x 10(-7) atm), DeltaG(+ +) for C-H addition to ((Me)PCP)Ir (plus a molecule of free H(2)) is very low (17.5 kcal/mol for CyH, 16.7 kcal/mol for n-PrH). Under these conditions, the rate-determining step for the D pathway is the loss of H(2) from ((Me)PCP)IrH(2) with DeltaG(D)(+ +) approximately DeltaH(D)(+ +) = 27.2 kcal/mol. For CyH, the calculated DeltaG(o)(+ +) for C-H addition to ((Me)PCP)IrH(2) on the A pathway is 35.2 kcal/mol (32.7 kcal/mol for n-PrH). At catalytic conditions, the calculated free energies of C-H addition are 31.3 and 33.7 kcal/mol for CyH and n-PrH addition, respectively. Elimination of H(2) from the resulting "seven-coordinate" Ir-species must proceed with an activation enthalpy at least as large as the enthalpy change of the elimination step itself (DeltaH approximately 11-13 kcal/mol), and with a small entropy of activation. The free energy of activation for H(2) elimination (DeltaG(A)(+ +)) is hence found to be greater than ca. 36 kcal/mol for both CyH and n-PrH under catalytic conditions. The overall free-energy barrier of the A pathway is calculated to be higher than that of the D pathway by ca. 9 kcal/mol. Reversible C-H(D) addition to ((R)PCP)IrH(2) is predicted to lead to H/D exchange, because the barriers for hydride scrambling are extremely low in the "seven-coordinate" polyhydrides. In agreement with calculation, H/D exchange is observed experimentally for several deuteriohydrocarbons with the following order of rates: C(6)D(6) > mesitylene-d(12) > n-decane-d(22) > cyclohexane-d(12). Because H/D exchange in cyclohexane-d(12) solution is not observed even after 1 week at 180 degrees C, we estimate that the experimental barrier to cyclohexane C-D addition is greater than 36.4 kcal/mol. This value is considerably greater than the experimental barrier for the full catalytic dehydrogenation cycle for cycloalkanes (ca. 31 kcal/mol). Thus, the experimental evidence, in agreement with calculation, strongly indicates that the A pathway is not kinetically viable as a segment of the "acceptorless" dehydrogenation cycle.     

Methane-to-methanol oxidation by the hydrated iron(IV) oxo species in aqueous solution: a combined DFT and car-parrinello molecular dynamics study

Previously, we have shown that the ferryl ion ([FeIVO]2+) is easily produced from Fenton's reagent (i.e., a mixture of Fe2+ ions and H2O2 in aqueous solution), using DFT and Car-Parrinello MD calculations. To verify that the ferryl ion can indeed act as the active species in oxidation reactions with Fenton's reagent, we study in the present paper the reactivity of the ferryl ion toward an organic substrate, in particular the oxidation of methane to methanol. In the first part of this paper, we perform static DFT calculations on the reaction of CH4 with the [(H2O)5FeIVO]2+ complex in vacuo that show a strong prevalence of the oxygen-rebound mechanism over the methane coordination mechanism. This is in agreement with the static DFT results for methane oxidation by biocatalysts MMO and P450, but not with those for methane oxidation by bare metal-oxo ions, where the methane coordination mechanism prevails. The highest energy barrier in the oxygen-rebound mechanism is only 3 kcal/mol in vacuo, whereas in the methane coordination mechanism the highest barrier is 23 kcal/mol. Overall the oxidation reaction energy is downhill by 47 kcal/mol. We conclude that the ferryl ion can indeed act as the oxidative intermediate in the Fenton oxidation of organic species. In the second part of this paper, we perform a preliminary assessment of solvent effects on the oxidation by the ferryl ion in aqueous solution using the method of constrained (first principles) molecular dynamics. The free energy barrier of the H-abstraction reaction from methane by the ferryl ion (i.e., the first step in the rebound mechanism) in aqueous solution is, with 22 kcal/mol in solution, significantly higher than in vacuo. Given the fact that methane has a relatively strong C-H bond (ca. 10 kcal/mol stronger than the C-H bonds in the more typical Fenton's reagent substrates), we infer that for many organic substrates oxidation with the ferryl ion as an active intermediate may be a perfectly viable route.     

Substitution of hydrogen by deuterium changes the regioselectivity of ethylbenzene hydroxylation by an oxo-iron-porphyrin catalyst

Heme oxo-iron complexes are powerful oxygenation catalysts of environmentally benign hydroxylation processes. We have performed density functional theoretic calculations on a model system, that is, an oxo-iron-porphyrin (Por) complex [(Fe=O)Cl(Por)], and studied its reactivity toward a realistic substrate, namely, ethylbenzene. The calculations showed that the dominant reaction process in the gas phase is benzyl hydroxylation leading to 1-phenylethanol, with an energetic barrier of 9.1 kcal mol(-1), while the competing para-phenyl hydroxylation has a barrier 3.0 kcal mol(-1) higher in energy. This benzyl hydroxylation barrier is the lowest C-H hydroxylation barrier we have obtained so far for oxo-iron-porphyrin complexes. Due to electronic differences between the intermediates in the phenyl and benzyl hydroxylation processes, the phenyl hydroxylation process is considerably stabilised over the benzyl hydroxylation mechanism in environments with a large dielectric constant. In addition, we calculated kinetic isotope effects of the substitution of one or more hydrogen atoms of ethylbenzene by deuterium atoms and studied its effect on the reaction barriers. Thus, in a medium with a large dielectric constant, a regioselectivity change occurs between [H(10)]ethylbenzene and [D(10)]ethylbenzene whereby the deuterated species gives phenol products whereas the hydrogenated species gives mainly 1-phenylethanol products. This remarkable metabolic switching was analysed and found to occur due to 1) differences in strength between a C-H versus a C-D bond and 2) stabilisation of cationic intermediates in a medium with a large dielectric constant. We have compared our calculations with experimental work on synthetic oxo-iron-porphyrin catalysts as well as with enzyme-reactivity studies.