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.     

How does product isotope effect prove the operation of a two-state "rebound" mechanism in C-H hydroxylation by cytochrome P450?

C-H hydroxylation is a fundamental process. In Nature it is catalyzed by the enzyme cytochrome P450, in a still-debated mechanism that poses a major intellectual challenge for both experiment and theory; currently, the opinions keep swaying between the original single-state rebound mechanism, a two-oxidant mechanism (where ferric peroxide participates as a second oxidant, in addition to the primary active species, the high-valent iron-oxo species), and two-state reactivity (TSR) mechanism (where two spin states are involved). Recent product isotope effect (PIE) measurements for the trans-2-phenyl-methyl cyclopropane probe (1), led Newcomb and co-workers (Newcomb, M.; Aebisher, D.; Shen, R.; Esala, R.; Chandrasena, P.; Hollenberg, P.; Coon, M. J. J. Am. Chem. Soc. 2003, 125, 6064-6065) to rule out TSR in favor of the two-oxidant scenario, since the direction of the PIE was at odds with the one predicted from calculations on methane hydroxylation. The present report describes a density functional theoretical study of C-H hydroxylation of the Newcomb probe, 1, leading to rearranged (3) and unrearranged (2) products. Our study shows that the reaction occurs via TSR in which the high-spin pathway gives dominant rearranged products, whereas the low-spin pathway favors unrearranged products. The calculated PIE(2/3) values based on TSR are found to be in excellent agreement with the experimental data of Newcomb and co-workers. This match between experiment and theory makes a strong case that the reaction occurs via TSR mechanism.     

Do spin state changes matter in organometallic chemistry? A computational study

Spin changes occur often in organometallic chemistry, and their effect on kinetics is not well understood. We report computations on the singlet and triplet potential energy surfaces of several processes of this type and show that the topology of the individual surfaces, as well as of the crossing regions between them, can be used to rationalize the observed reactivity in all cases. In particular, the slow addition of dihydrogen to W[N(CH(2)CH(2)NSiMe(3))(3)]H (Schrock, R. R.; Shih, K. Y.; Dobbs, D. A.; Davis, W. M. J. Am. Chem. Soc. 1995, 117, 6609) is shown to be a "spin-blocked" reaction with a high barrier due to the crossing between reactant triplet and product singlet surfaces. In contrast, addition of CO to TpCo(CO) (Detrich, J. L.; Reinaud, O. M.; Rheingold, A. L.; Theopold, K. H. J. Am. Chem. Soc. 1995, 117, 11745) is fast because the triplet and singlet surfaces cross at low energy. Particular care is taken to use DFT methods which are in adequate agreement with experimental and high-level computational energetics for these systems.     

A two-state reactivity rationale for counterintuitive axial ligand effects on the C-H activation reactivity of nonheme FeIV=O oxidants

This paper addresses the observation of counterintuitive reactivity patterns of iron-oxo reagents, TMC(L)FeO(2+,1+); L=CH(3)CN, CF(3)CO(2) (-), N(3) (-), and SR(-), in O-transfer to phosphines versus H-abstraction from, for example, 1,4-cyclohexadiene. Experiments show that O-transfer reactivity correlates with the electrophilicity of the oxidant, but H-abstraction reactivity follows an opposite trend. DFT/B3 LYP calculations reveal that two-state reactivity (TSR) serves as a compelling rationale for these trends, whereby all reactions involve two adjacent spin-states of the iron(IV)-oxo species, triplet and quintet. The ground state triplet surface has high barriers, whereas the excited state quintet surface features lower ones. The barriers, on any single surface, are found to increase as the electrophilicity of TMC(L)FeO(2+,1+) decreases. Thus, the counterintuitive behavior of the H-abstraction reactions cannot be explained by considering the reactivity of only a single spin state but can be rationalized by a TSR model in which the reactions proceed on the two surfaces. Two TSR models are outlined: one is traditional involving a variable transmission coefficient for crossover from triplet to quintet, followed by quintet-state reactions; the other considers the net barrier as a blend of the triplet and quintet barriers. The blending coefficient (x), which estimates the triplet participation, increases as the quintet-triplet energy gap of the TMC(L)FeO(2+,1+) reagent increases, in the following order of L: CH(3)CN > CF(3)CO(2) (-) > N(3) (-) > SR(-). The calculated barriers predict the dichotomic experimental trends and the counterintuitive behavior of the H-abstraction series. The TSR approaches make a variety of testable predictions.     

A rational basis for the axial ligand effect in C-H oxidation by [MnO(porphyrin)(X)]+ (X = H2O, OH-, O2-) from a DFT study

Oxyl radical character in the MnO group of the title system is shown from a density functional theory study to be essential for efficient C-H cleavage, which is a key step in C-H oxidation. Since oxyl species have elongated Mn-O bonds relative to the more usual oxo species of type MnO, the normal expectation would be that high trans-influence ligands X should facilitate oxyl character by elongating the Mn-O bond and thus enhance both oxyl character and reactivity. Contrary to this expectation, but in line with the experimental data (Jin, N.; Ibrahim, M.; Spiro, T. G.; Groves, J. T. J. Am. Chem. Soc. 2007, 129, 12416), we find that reactivity increases along the series X = O(2-) < OH(-) < H2O for the following reasons. The ground-state singlet (S) is unreactive for all X, and only the higher-energy triplet (T) and quintet (Q) states have the oxyl character needed for reactivity, but the higher trans-influence X ligands are also shown to increase the S/T and S/Q gaps, thus making attainment of the needed T and Q states harder. The latter effect is dominant, and high trans-influence X ligands thus disfavor reaction. The higher reactivity in the presence of acid noted by Groves and co-workers is thus rationalized by the preference for having X = H2O over OH(-) or O(2-).     

Mode selective photodissociation dynamics in V+(OCO)

The electrostatic V+(OCO) complex has a vibrationally resolved photodissociation spectrum in the visible. Photodissociation produces V+ + CO2 (nonreactive pathway) and VO+ +CO (reactive pathway). Production of VO+ is energetically favored, but spin forbidden. One-photon dissociation studies confirm mode selectivity observed by Lessen et al. [J. Chem. Phys. 95, 1414 (1991)]: excitation of one quantum of rocking motion enhances VO+ production by >30%. Branching ratio measurements in one-photon dissociation are extended to higher energy. The effect of OCO antisymmetric stretch vibrations on reactivity is investigated using vibrationally mediated photodissociation, in which the OCO antisymmetric stretch is excited at 2390.9 cm(-1). Vibrationally excited molecules are then dissociated in the visible. Seven vibronic bands are investigated, involving the antisymmetric stretch alone and in combination with the CO2 bend, the V+(OCO) stretch and rock. Exciting the antisymmetric stretch leads to a approximately 15% increase in the reactive VO+ channel, compared to other states at similar energy. Combination bands involving the antisymmetric stretch all show slightly higher reactivity. Electronic structure calculations were performed to characterize the dissociation pathways and excited electronic states of V+(OCO). The geometries of reactants, products, and transition states and relative energies of quintet and triplet states were determined using hybrid density functional theory; energies were also calculated using the coupled cluster with single, double and perturbative triple excitations method. In addition, time-dependent density functional theory calculations were performed to predict the excited electronic states of quintet and triplet V+(OCO). Spin-orbit coupling of quintet states to triplet states was calculated and used to compute intersystem crossing rates, which reproduce many of the observed mode selective trends. The V+--OCO stretch and OCO antisymmetric stretch appear to enhance reactivity by increasing the intersystem crossing rate.     

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.     

Theoretical Investigation for Two-state Reactivity of CO_2 Hydrogenation Catalyzed by Ru in Gas Phase

Gas-phase CO_2 catalyzed activation hydrogenation by Ru atoms was studied with density functional theory. Based on the structure optimization of different potential energy surfaces,there are two crossing points between singlet and triplet potential energy surfaces and there is a crossing point between quintet and triplet potential energy surfaces in the whole catalytic cycle. Spin transition probabilities in the vicinity of the intersections have been calculated by the Landau-Zener model theory. There are three minimum energy crossing points(MECPs) with strong spin-orbital coupling effect and higher spin transition probability,and all spin inversion occurred in s orbital and different d orbitals of ruthenium,indicating this is a typical two-state reactivity(TSR) reaction. Finally,the lowest energy reaction path is ensured.     

Conical intersection of the ground and first excited states of water: energies and reduced density matrices from the anti-Hermitian contracted Schrödinger equation

A conical intersection between the ground and first-excited states of water is computed through the direct calculation of two-electron reduced density matrices (2-RDMs) from solutions of the anti-Hermitian contracted Schr?dinger equation (ACSE). This study is an extension of a previous study in which the ACSE was used to compute the energies around a conical intersection in the triplet excited states of methylene [Snyder, J. W., Jr.; Rothman, A. E.; Foley, J. J.; Mazziotti, D. A. J. Chem. Phys. 2010, 132, 154109]. We compute absolute energies of the 1(1)A' and 2(1)A' states of water (H(2)O) and the location of the conical intersection. The ACSE energies are compared to those from ab initio wave function methods. To treat multireference correlation, we seed the ACSE with an initial 2-RDM from a multiconfiguration self-consistent field (MCSCF) calculation. Unlike the situation for methylene, the two states in the vicinity of the conical intersection of water both have the same spatial symmetry. Hence, the study demonstrates the ability of the ACSE to resolve states of the same spatial symmetry that are nearly degenerate in energy. The 2-RDMs from the ACSE nearly satisfy necessary N-representability conditions. Comparison of the results from double-ζ and augmented double-ζ basis sets demonstrates the importance of augmented (or diffuse) functions for determining the location of the conical intersection.     

Theoretical study of the mechanism of oxoiron(IV) formation from H2O2 and a nonheme iron(II) complex: O-O cleavage involving proton-coupled electron transfer

It has recently been shown that the nonheme oxoiron(IV) species supported by the 1,4,8,11-tetramethyl-1,4,8,11-tetraazacyclotetradecane ligand (TMC) can be generated in near-quantitative yield by reacting [Fe(II)(TMC)(OTf)(2)] with a stoichiometric amount of H(2)O(2) in CH(3)CN in the presence of 2,6-lutidine (Li, F.; England, J.; Que, L., Jr. J. Am. Chem. Soc. 2010, 132, 2134-2135). This finding has major implications for O-O bond cleavage events in both Fenton chemistry and nonheme iron enzymes. To understand the mechanism of this process, especially the intimate details of the O-O bond cleavage step, a series of density functional theory (DFT) calculations and analyses have been carried out. Two distinct reaction paths (A and B) were identified. Path A consists of two principal steps: (1) coordination of H(2)O(2) to Fe(II) and (2) a combination of partial homolytic O-O bond cleavage and proton-coupled electron transfer (PCET). The latter combination renders the rate-limiting O-O cleavage effectively a heterolytic process. Path B proceeds via a simultaneous homolytic O-O bond cleavage of H(2)O(2) and Fe-O bond formation. This is followed by H abstraction from the resultant Fe(III)-OH species by an ?OH radical. Calculations suggest that path B is plausible in the absence of base. However, once 2,6-lutidine is added to the reacting system, the reaction barrier is lowered and more importantly the mechanistic path switches to path A, where 2,6-lutidine plays an essential role as an acid-base catalyst in a manner similar to how the distal histidine or glutamate residue assists in compound I formation in heme peroxidases. The reaction was found to proceed predominantly on the quintet spin state surface, and a transition to the triplet state, the experimentally known ground state for the TMC-oxoiron(IV) species, occurs in the last stage of the oxoiron(IV) formation process.     

(Salen)Mn(III)-catalyzed epoxidation reaction as a multichannel process with different spin states. Electronic tuning of asymmetric catalysis: a theoretical study

The (salen)Mn(III)-catalyzed epoxidation reaction mechanism has been investigated using density functional theory (DFT). There is considerable interest in and controversy over the mechanism of this reaction. The results of experimental studies have offered some support for three different reaction mechanisms: concerted, stepwise radical, and metallooxetane mediated. In this paper, a theoretical examination of the reaction suggests a novel mechanism that describes the reaction as a multichannel process combining both concerted and stepwise radical pathways. The competing channels have different spin states: the singlet, the triplet, and the quintet. The singlet reaction pathway corresponds to a concerted mechanism and leads exclusively to a cis epoxide product. In contrast, the triplet and quintet reactions follow a stepwise mechanism and lead to a product mixture of cis and trans epoxides. We show that the experimentally observed dependence of isomer product ratios on electronic effects connected with the substitution of the catalyst ligands is due to changing the relative position and, hence, the relative activities of the channels with different cis-trans yields. Because the results and conclusions of the present work dramatically differ from the results and conclusion of the recent DFT theoretical investigation (Linde, C.; Akermark, B; Norrby, P.-O.; Svensson, M. J. Am. Chem. Soc. 1999, 121, 5083.), we studied possible sources for the deep contradictions between the two works. The choice of the DFT functional and a model has been shown to be crucial for accurate results. Using high level ab initio calculations (coupled cluster-CCSD(T)), we show that the computational procedure employed in this study generates significantly more reliable numerical results. It is also shown that the smaller cationic model without a chlorine ligand that was used by Linde et al. is too oversimplified with respect to our larger neutral model. For this reason, using the cationic model led to a qualitatively wrong quintet reaction profile that played a key role in theoretical postulates in the earlier work.     

Extension of the renormalized coupled-cluster methods exploiting left eigenstates of the similarity-transformed hamiltonian to open-shell systems: a benchmark study

The recently formulated completely renormalized coupled-cluster method with singles, doubles, and noniterative triples, exploiting the biorthogonal form of the method of moments of coupled-cluster equations (Piecuch, P.; W?och, M. J. Chem. Phys. 2005, 123, 224105; Piecuch, P.; W?och, M.; Gour, J. R.; Kinal, A. Chem. Phys. Lett. 2006, 418, 467), termed CR-CC(2,3), is extended to open-shell systems. Test calculations for bond breaking in the OH radical and the F2+ ion and singlet-triplet gaps in the CH2, HHeH, and (HFH)- biradical systems indicate that the CR-CC(2,3) approach employing the restricted open-shell Hartree--Fock (ROHF) reference is significantly more accurate than the widely used CCSD(T) method and other noniterative triples coupled-cluster approximations without making the calculations substantially more expensive. A few molecular examples, including the activation energies of the C2H4 + H --> C2H5 forward and reverse reactions and the triplet states of the CH2 and H2Si2O2 biradicals, are used to show that the dependence of the ROHF-based CR-CC(2,3) energies on the method of canonicalization of the ROHF orbitals is, for all practical purposes, negligible.     

What factors influence the ratio of C-H hydroxylation versus C=C epoxidation by a nonheme cytochrome P450 biomimetic?

Density functional calculations on a nonheme biomimetic (Fe=O(TMCS+) have been performed and its catalytic properties versus propene investigated. Our studies show that this catalyst is able to chemoselectively hydroxylate C=H bonds even in the presence of C=C double bonds. This phenomenon has been analyzed and found to occur due to Pauli repusions between protons on the TMCS ligand with protons attached to the approaching substrate. The geometries of the rate determining transition states indicate that the steric hindrance is larger in the epoxidation transition states than in the hydroxylation ones with much shorter distances; hence the hydroxylation pathway is favored over the epoxidation. Although, the reactant experiences close lying triplet and quintet spin states, the dominant reaction mechanism takes place on the quintet spin state surface; i.e., Fe=O(TMCS)+ reacts via single-state reactivity. Our calculations show that this spin state selectivity is the result of geometric orientation of the transition state structures, whereby the triplet ones are destabilized by electrostatic repulsions between the substrate and the ligand while the quintet spin transition states are aligned along the ideal axis. The reactivity patterns and geometries are compared with oxoiron species of dioxygenase and monoxygenase enzymes. Thus, Fe=O(TMCS)+ shows some similarities with P450 enzyme reactivity: it chemoselectively hydroxylates C=H bonds even in the presence of a C=C double bond and therefore is an acceptable P450 biomimetic. However, the absolute barriers of substrate oxidation by Fe=O(TMCS)+ are higher than the ones obtained with heme enzymes, but the chemoselectivity is lesser affected by external perturbations such as hydrogen bonding of a methanol molecule toward the thiolate sulfur or a dielectric constant. This is the first oxoiron complex whereby we calculated a chemoselective hydroxylation over epoxidation in the gas phase.     

Kinetic analysis of the reaction between (V2O5)n=1,2+ and ethylene

A systematic experimental and theoretical investigation of the influence of reactant energy on the reactivity of (V(2)O(5))(n)=1,2(+) clusters with ethylene (Justes, D. R.; Mitri?, R.; Moore, N. A.; Bonaci?-Koutecky, V.; Castleman, A. W., Jr. J. Am. Chem. Soc., 2003, 125, 6289) provided evidence of the rate controlling steps in the reaction. Herein, we present further experimental and theoretical evidence for our recently proposed radical cation mechanism for oxygen atom transfer from (V(2)O(5))(n)=1,2(+) clusters to ethylene. In particular the results of ab initio molecular dynamics simulations are found to further support the radical cation mechanism. Experimental reaction cross sections at the zero pressure limit and rate coefficients show that the energy dependence of the reaction cross section is in accord with the Langevin formula. Evidence is presented that ion-molecule association is the rate determining step, whereas subsequent hydrogen transfer and formation of acetaldehyde proceed without significant barriers. We propose a kinetic model for the reaction cross section that fully accounts for the experimental findings. The model offers the prospect of elucidating the details of the general reaction mechanisms through a study of the energy dependence of the reaction cross sections.     

Reactivity of Compound II: Electronic Structure Analysis of Methane Hydroxylation by Oxoiron(IV) Porphyrin Complexes

The methane hydroxylation reaction by a Compound II (Cpd II) mimic PorFe(IV)=O and its hydrosulfide-ligated derivative [Por(SH)Fe(IV)=O](-) is investigated by density functional theory (DFT) calculations on the ground triplet and excited quintet spin-state surfaces. On each spin surface both the σ- and π-channels are explored. H-abstraction is invariably the rate-determining step. In the case of PorFe(IV)=O the H-abstraction reaction can proceed either through the classic π-channel or through the nonclassical σ-channel on the triplet surface, but only through the classic σ-mechanism on the quintet surface. The barrier on the quintet σ-pathway is much lower than on the triplet channels so the quintet surface cuts through the triplet surfaces and a two state reactivity (TSR) mechanism with crossover from the triplet to the quintet surface becomes a plausible scenario for C-H bond activation by PorFe(IV)=O. In the case of the hydrosulfide-ligated complex the H-abstraction follows a π-mechanism on the triplet surface: the σ* is too high in energy to make a σ-attack of the substrate favorable. The σ- and π-channels are both feasible on the quintet surface. As the quintet surface lies above the triplet surface in the entrance channel of the oxidative process and is highly destabilized on both the σ- and π-pathways, the reaction can only proceed on the triplet surface. Insights into the electron transfer process accompanying the H-abstraction reaction are achieved through a detailed electronic structure analysis of the transition state species and the reactant complexes en route to the transition state. It is found that the electron transfer from the substrate σ(CH) into the acceptor orbital of the catalyst, the Fe-O σ* or π*, occurs through a rather complex mechanism that is initiated by a two-orbital four-electron interaction between the σ(CH) and the low-lying, oxygen-rich Fe-O σ-bonding and/or Fe-O π-bonding orbitals of the catalyst.     

A proton-shuttle mechanism mediated by the porphyrin in benzene hydroxylation by cytochrome p450 enzymes

Benzene hydroxylation is a fundamental process in chemical catalysis. In nature, this reaction is catalyzed by the enzyme cytochrome P450 via oxygen transfer in a still debated mechanism of considerable complexity. The paper uses hybrid density functional calculations to elucidate the mechanisms by which benzene is converted to phenol, benzene oxide, and ketone, by the active species of the enzyme, the high-valent iron-oxo porphyrin species. The effects of the protein polarity and hydrogen-bonding donation to the active species are mimicked, as before (Ogliaro, F.; Cohen, S.; de Visser, S. P.; Shaik, S. J. Am. Chem. Soc. 2000, 122, 12892-12893). It is verified that the reaction does not proceed either by hydrogen abstraction or by initial electron transfer (Ortiz de Montellano, P. R. In Cytochrome P450: Structure, Mechanism and Biochemistry, 2nd ed.; Ortiz de Montellano, P. R., Ed.; Plenum Press: New York, 1995; Chapter 8, pp 245-303). In accord with the latest experimental conclusions, the theoretical calculations show that the reactivity is an interplay of electrophilic and radicalar pathways, which involve an initial attack on the pi-system of the benzene to produce sigma-complexes (Korzekwa, K. R.; Swinney, D. C.; Trager, W. T. Biochemistry 1989, 28, 9019-9027). The dominant reaction channel is electrophilic and proceeds via the cationic sigma-complex,( 2)3, that involves an internal ion pair made from a cationic benzene moiety and an anionic iron porphyrin. The minor channel proceeds by intermediacy of the radical sigma-complex, (2)2, in which the benzene moiety is radicalar and the iron-porphyrin moiety is neutral. Ring closure in these intermediates produces the benzene oxide product ((2)4), which does not rearrange to phenol ((2)7) or cyclohexenone ((2)6). While such a rearrangement can occur post-enzymatically under physiological conditions by acid catalysis, the computations reveal a novel mechanism whereby the active species of the enzyme catalyzes directly the production of phenol and cyclohexenone. This enzymatic mechanism involves proton shuttles mediated by the porphyrin ring through the N-protonated intermediate, (2)5, which relays the proton either to the oxygen atom to form phenol ((2)7) or to the ortho-carbon atom to produce cyclohexenone product ((2)6). The formation of the phenol via this proton-shuttle mechanism will be competitive with the nonenzymatic conversion of benzene oxide to phenol by external acid catalysis. With the assumption that (2)5 is not fully thermalized, this novel mechanism would account also for the observation that there is a partial skeletal retention of the original hydrogen of the activated C-H bond, due to migration of the hydrogen from the site of hydroxylation to the adjacent carbon (so-called "NIH shift" (Jerina, D. M.; Daly, J. W. Science 1974, 185, 573-582)). Thus, in general, the computationally discovered mechanism of a porphyrin proton shuttle suggests thatthere is an enzymatic pathway that converts benzene directly to a phenol and ketone, in addition to nonenzymatic production of these species by conversion of arene oxide to phenol and ketone. The potential generality of protonated porphyrin intermediates in P450 chemistry is discussed in the light of the H/D exchange observed during some olefin epoxidation reactions (Groves, J. T.; Avaria-Neisser, G. E.; Fish, K. M.; Imachi, M.; Kuczkowski, R. J. Am. Chem. Soc. 1986, 108, 3837-3838) and the general observation of heme alkylation products (Kunze, K. L.; Mangold, B. L. K.; Wheeler, C.; Beilan, H. S.; Ortiz de Montellano, P. R. J. Biol. Chem. 1983, 258, 4202-4207). The competition, similarities, and differences between benzene oxidation viz. olefin epoxidation and alkanyl C-H hydroxylation are discussed, and comparison is made with relevant experimental and computational data. The dominance of low-spin reactivity in benzene hydroxylation viz. two-state reactivity (Shaik, S.; de Visser, S. P.; Ogliaro, F.; Schwarz, H.; Schr?der, D. Curr. Opin. Chem. Biol. 2002, 6, 556-567) in olefin epoxidation and alkane hydroxylation is traced to the loss of benzene resonance energy during the bond activation step.     

Can a single oxidant with two spin states masquerade as two different oxidants? A study of the sulfoxidation mechanism by cytochrome p450

Density functional calculations were performed on the sulfoxidation reaction by a model compound I (Cpd I) of cytochrome P450. By contrast to previous alkane hydroxylation studies, which exhibit a dominant low-spin (LS) pathway, the sulfoxidation follows a dominant high-spin (HS) reaction. Thus, competing hydroxylation and sulfoxidation processes as observed for instance by Jones et al. (Volz, T. J.; Rock, D. A.; Jones, J. P. J. Am. Chem. Soc. 2002, 124, 9724) are the result of a two-state reactivity scenario, whereby the hydroxylation originates from the LS pathway and the sulfoxidation from the HS pathway. In this manner, two spin states of a single oxidant (Cpd I) can be disguised as two different oxidants. The calculations rule out the possibility that a second oxidant (the ferric peroxide, Cpd 0 species) interferes in the observed results of Jones et al.     

The high-valent iron-oxo species of polyoxometalate, if it can be made, will be a highly potent catalyst for C-H hydroxylation and double-bond epoxidation

This study uses density functional theory (DFT) calculations to explore the reactivity of the putative high-valent iron-oxo reagent of the iron-substituted polyoxometalate (POM-FeO4-), derived from the Keggin species, PW12O40(3-). It is shown that POM-FeO4- is in principle capable of C-H hydroxylation and C=C epoxidation and that it should be a powerful oxidant, even more so than the Compound I species of cytochrome P450. The calculations indicate that in a solvent, the barriers, and especially those for epoxidation, become sufficiently small that one may expect an extremely fast reaction. An experimental investigation (by R.N. and A.M.K.) shows, however, that the formation of POM-FeO4- using the oxygen donor, F5PhI-O, leads to a persistent adduct, POM-FeO-I-PhF5(4-), which does not decompose to POM-FeO4- + F5Ph-I at the working temperature and exhibits sluggish reactivity, in accord with previous experimental results (Hill, C. L.; Brown, R. B., Jr. J. Am. Chem. Soc. 1986, 108, 536 and Mansuy, D.; Bartoli, J.-F.; Battioni, P.; Lyon, D. K.; Finke, R. G. J. Am. Chem. Soc. 1991, 113, 7222). Subsequent calculations indeed reveal that the gas-phase binding energy of F5PhI to POM-FeO4- is high (ca. 20 kcal/mol) compared to the corresponding binding energy of propene (ca. 2-3 kcal/mol). As such, the POM-FeO-I-PhF5(4-) complex is expected to be persistent toward the displacement of F5PhI by a substrate like propene, leading thereby to sluggish oxidative reactivity. According to theory, overcoming this technical difficulty may turn out to be very rewarding. The question is, can POM-FeO4- be made?     

Stability, reactivity, and aromaticity of compounds of a multivalent superatom

In this article, we analyze the stability, reactivity, and possible aromatic behavior of two recently reported clusters (Reveles, J. U.; Khanna, S. N.; Roach, P. J.; Castleman, A. W., Jr. Proc. Natl. Acad. Sci. 2006, 103, 18405), viz., Al(7)C(-) and Al(7)O(-) in the light of the principles of the maximum hardness and minimum electrophilicity as well as the nucleus-independent chemical shift values. Stability of these clusters in the context of addition/removal of an electron or an Al atom is now clearly understood.