Role of the somersault rearrangement in the oxidation step for flavin monooxygenases (FMO). A comparison between FMO and conventional xenobiotic oxidation with hydroperoxides

Model quantum mechanical calculations presented for C-4a-flavin hydroperoxide (FlHOOH) at the B3LYP/6-311+G(d,p) level suggest a new mechanism for flavoprotein monooxygenase (FMO) oxidation involving a concerted homolytic O-O bond cleavage in concert with hydroxyl radical transfer from the flavin hydroperoxide rather than an S(N)2-like displacement by the substrate on the C-4a-hydroperoxide OOH group. Homolytic O-O bond cleavage in a somersault-like rearrangement of hydroperoxide C-4a-flavinhydroperoxide (1) (FLHO-OH → FLHO···HO) produces an internally hydrogen-bonded HO(?) radical intermediate with a classical activation barrier of 27.0 kcal/mol. Model hydroperoxide 1 is used to describe the transition state for the key oxidation step in the paradigm aromatic hydroxylase, p-hydroxybenzoate hydroxylase (PHBH). A comparison of the electron distribution in the transition structures for the PHBH hydroxylation of p-hydroxybenzoic acid (ΔE(?) = 23.0 kcal/mol) with that of oxidation of trimethylamine (ΔE(?) = 22.3 kcal/mol) and dimethyl sulfide (ΔE? = 14.1 kcal/mol) also suggests a mechanism involving a somersault mechanism in concert with transfer of an HO(?) radical to the nucleophilic heteroatom center with a hydrogen transfer back to the FLH-O residue after the barrier is crossed to produce the final product, FLH-OH. In each case the hydroxylation barrier was less than that of the O-O rearrangement barrier in the absence of a substrate supporting an overall concerted process. All three transition structures bear a resemblance to the TS for the comparable hydroxylation of isobutane (ΔE(?) = 29.2 kcal/mol) and for simple Fenton oxidation by aqueous iron(III) hydroperoxides. To our surprise the oxidation of N- and S-nucleophiles with conventional oxidants such as alkyl hydroperoxides and peracids also proceeds by HO(?) radical transfer in a manner quite similar to that for tricyclic hydroperoxide 1. Stabilization of the developing oxyradical produced by somersault rearrangement for concerted enzymatic oxidation with tricyclic hydroperoxide 1 results in a reduced overall activation barrier.     

What is the active species of cytochrome P450 during camphor hydroxylation? QM/MM studies of different electronic states of compound I and of reduced and oxidized iron-oxo intermediates

We have investigated C-H hydroxylation of camphor by Compound I (Cpd I) of cytochrome P450cam in different electronic states and by its one-electron reduced and oxidized forms, using QM/MM calculations in the native protein/solvent environment. Cpd I species with five unpaired electrons (pentaradicaloids) are ca. 12 kcal/mol higher in energy than the ground state Cpd I species with three unpaired electrons (triradicaloids). The H-abstraction transition states of pentaradicaloids lie ca. 21 (9) kcal/mol above the triradicaloid (pentaradicaloid) reactants. Hydroxylation via pentaradicaloids is thus facile provided that they can react before relaxing to the ground-state triradicaloids. Excited states of Cpd I with an Fe(V)-oxo moiety lie more than 20 kcal/mol above the triradicaloid ground state in single-point gas-phase calculations, but these electronic configurations are not stable upon including the point-charge protein environment which causes SCF convergence to the triradicaloid ground state. One-electron reduced species (Cpd II) show sluggish reactivity compared with Cpd I in agreement with experimental model studies. One-electron oxidized species are more reactive than Cpd I but seem too high in energy to be accessible. The barriers to hydrogen abstraction for the various forms of Cpd I are generally not affected much by the chosen protonation states of the Asp297 and His355 residues near the propionate side chains of the heme or by the appearance of radical character at Asp297, His355, or the propionates.     

Importance of cytochromes in cyclization reactions: Quantum chemical study on a model reaction of proguanil to cycloguanil

Proguanil, an anti‐malarial prodrug, undergoes cytochrome P450 catalyzed biotransformation to the pharmacologically active triazine metabolite (cycloguanil), which inhibits plasmodial dihydrofolate reductase. This cyclization is catalyzed by CYP2C19 and many anti‐malarial lead compounds are being designed and synthesized to exploit this pathway. Quantum chemical calculations were performed using the model species (Cpd I for active species of cytochrome and N4‐isopropyl‐N6‐methylbiguanide for proguanil) to elucidate the mechanism of the cyclization pathway. The overall reaction involves the loss of a water molecule, and is exothermic by approximately 55 kcal/mol, and involves a barrier of approximately 17 kcal/mol. The plausible reaction pathway involves the initial H‐radical abstraction from the isopropyl group by Cpd I, followed by two alternative paths‐ (i) oxygen rebound to provide hydroxyl derivative and (ii) loss of additional H‐radical to yield 1,3,5‐triazatriene, which undergoes cyclization. This study helped in understanding the role of the active species of cytochromes in this important cyclization reaction. © 2014 Wiley Periodicals, Inc.     

Chemical behavior of the biradicaloid (HO...ONO) singlet states of peroxynitrous acid. The oxidation of hydrocarbons, sulfides, and selenides

Various high levels of theory have been applied to the characterization of two higher lying biradicaloid metastable singlet states of peroxynitrous acid. A singlet minimum (cis-2) was located that had an elongated O-O distance (2.17 A) and was only 12.2 kcal/mol [UB3LYP/6-311+G(3df,2p)+ZPVE] higher in energy than its ground-state precursor. A trans-metastable singlet (trans-2) was 10.9 kcal/mol higher in energy than ground-state HO-ONO. CASSCF(12,10)/6-311+G(d,p) calculations predict the optimized geometries of these cis- and trans-metastable singlets to be close to those obtained with DFT. Optimization of cis- and trans-2 within the COSMO solvent model suggests that both exist as energy minima in polar media. Both cis- and trans-2 exist as hydrogen bonded complexes with several water molecules. These collective data suggest that solvated forms of cis-2.3H(2)O and trans-2.3H(2)O represent the elusive higher lying biradicaloid minima that were recently (J. Am. Chem. Soc. 2003, 125, 16204) advocated as the metastable forms of peroxynitrous acid (HOONO). The involvement of metastable trans-2 in the gas phase oxidation of methane and isobutane is firmly established to take place on the unrestricted [UB3LYP/6-311+G(d,p)] potential energy surface (PES) with classical activations barriers for the hydrogen abstraction step that are 15.7 and 5.9 kcal/mol lower than the corresponding activation energies for producing products methanol and tert-butyl alcohol formed on the restricted PES. The oxidation of dimethyl sulfide and dimethyl selenide, two-electron oxidations, proceeds by an S(N)2-like attack of the heteroatom lone pair on the O-O bond of ground-state peroxynitrous acid. No involvement of metastable forms of HO-ONO was discernible.     

Quantum mechanical/molecular mechanical study on the mechanisms of compound I formation in the catalytic cycle of chloroperoxidase: an overview on heme enzymes

The formation of Compound I (Cpd I), the active species of the enzyme chloroperoxidase (CPO), was studied using QM/MM calculation. Starting from the substrate complex with hydrogen peroxide, FeIII-HOOH, we examined two alternative mechanisms on the three lowest spin-state surfaces. The calculations showed that the preferred pathway involves heterolytic O-O cleavage that proceeds via the iron hydroperoxide species, i.e., Compound 0 (Cpd 0), on the doublet-state surface. This process is effectively concerted, with a barrier of 12.4 kcal/mol, and is catalyzed by protonation of the distal OH group of Cpd 0. By comparison, the path that involves a direct O-O cleavage from FeIII-HOOH is less favored. A proton coupled electron transfer (PCET) feature was found to play an important role in the mechanism nascent from Cpd 0. Initially, the O-O cleavage progresses in a homolytic sense, but as soon as the proton is transferred to the distal OH, it triggers an electron transfer from the heme-oxo moiety to form water and Cpd I. This study enables us to generalize the mechanisms of O-O activation, elucidated so far by QM/MM calculations, for other heme enzymes, e.g., cytochrome P450cam, horseradish peroxidase (HRP), nitric oxide synthase (NOS), and heme oxygenase (HO). Much like for CPO, in the cases of P450 and HRP, the PCET lowers the barrier below the purely homolytic cleavage alternative (in our case, the homolytic mechanism is calculated directly from FeIII-HOOH). By contrast, the absence of PCET in HO, along with the robust water cluster, prefers a homolytic cleavage mechanism.     

High-level ab initio studies of hydrogen abstraction from prototype hydrocarbon systems

Symmetric and nonsymmetric hydrogen abstraction reactions are studied using state-of-the-art ab initio electronic structure methods. Second-order M?ller-Plesset perturbation theory (MP2) and the coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] methods with large correlation consistent basis sets (cc-pVXZ, where X = D,T,Q) are used in determining the transition-state geometries, activation barriers, and thermodynamic properties of several representative hydrogen abstraction reactions. The importance of basis set, electron correlation, and choice of zeroth-order reference wave function in the accurate prediction of activation barriers and reaction enthalpies are also investigated. The ethynyl radical (*CCH), which has a very high affinity for hydrogen atoms, is studied as a prototype hydrogen abstraction agent. Our high-level quantum mechanical computations indicate that hydrogen abstraction using the ethynyl radical has an activation energy of less than 3 kcal mol(-1) for hydrogens bonded to an sp(2) or sp(3) carbon. These low activation barriers further corroborate previous studies suggesting that ethynyl-type radicals would make good tooltips for abstracting hydrogens from diamondoid surfaces during mechanosynthesis. Modeling the diamond C(111) surface with isobutane and treating the ethynyl radical as a tooltip, hydrogen abstraction in this reaction is predicted to be barrierless.     

QM/MM modeling of the hydroxylation of the androstenedione substrate catalyzed by cytochrome P450 aromatase (CYP19A1)

CYP19A1 aromatase is a member of the Cytochrome P450 family of hemeproteins, and is the enzyme responsible for the final step of the androgens conversion into the corresponding estrogens, via a three‐step oxidative process. For this reason, the inhibition of this enzyme plays an important role in the treatment of hormone‐dependent breast cancer. The first catalytic subcycle, corresponding to the hydroxilation of androstenedione, has been proposed to occur through a first hydrogen abstraction and a subsequent oxygen rebound step. In present work, we have studied the mechanism of the first catalytic subcycle by means of hybrid quantum mechanics/molecular mechanics methods. The inclusion of the protein flexibility has been achieved by means of Free Energy Perturbation techniques, giving rise to a free energy of activation for the hydrogen abstraction step of 13.5 kcal/mol. The subsequent oxygen rebound step, characterized by a small free energy barrier (1.5 kcal/mol), leads to the hydroxylated products through a highly exergonic reaction. In addition, an analysis of the primary deuterium kinetic isotopic effects, calculated for the hydrogen abstraction step, reveals values (～10) overpassing the semiclassical limit for the C? H, indicating the presence of a substantial tunnel effect. Finally, a decomposition analysis of the interaction energy for the substrate and cofactor in the active site is also discussed. According to our results, the role of the enzymatic environment consists of a transition state stabilization by means of dispersive and polarization effects. © 2015 Wiley Periodicals, Inc.     

QM/MM study of mechanisms for compound I formation in the catalytic cycle of cytochrome P450cam

In the catalytic cycle of cytochrome P450cam, after molecular oxygen binds as a ligand to the heme iron atom to yield a ferrous dioxygen complex, there are fast proton transfers that lead to the formation of the active species, Compound I (Cpd I), which are not well understood because they occur so rapidly. In the present work, the conversion of the ferric hydroperoxo complex (Cpd 0) to Cpd I has been investigated by combined quantum-mechanical/molecular-mechanical (QM/MM) calculations. The residues Asp(251) and Glu(366) are considered as proton sources. In mechanism I, a proton is transported to the distal oxygen atom of the hydroperoxo group via a hydrogen bonding network to form protonated Cpd 0 (prot-Cpd0: FeOOH(2)), followed by heterolytic O-O bond cleavage that generates Cpd I and water. Although a local minimum is found for prot-Cpd0 in the Glu(366) channel, it is very high in energy (more than 20 kcal/mol above Cpd 0) and the barriers for its decay are only 3-4 kcal/mol (both toward Cpd 0 and Cpd I). In mechanism II, an initial O-O bond cleavage followed by a concomitant proton and electron transfer yields Cpd I and water. The rate-limiting step in mechanism II is O-O cleavage with a barrier of about 13-14 kcal/mol. According to the QM/MM calculations, the favored low-energy pathway to Cpd I is provided by mechanism II in the Asp(251) channel. Cpd 0 and Cpd I are of similar energies, with a slight preference for Cpd I.     

Multi-state epoxidation of ethene by cytochrome P450: a quantum chemical study.

The epoxidation of ethene by a model for Compound I of cytochrome P450, studied by the use of density functional B3LYP calculations, involves two-state reactivity (TSR) with multiple electromer species, hence "multi-state epoxidation". The reaction is found to proceed in stepwise and effectively concerted manners. Several reactive states are involved; the reactant is an (oxo)iron(IV) porphyrin cation radical complex with two closely lying spin states (quartet and doublet), both of which react with ethene to form intermediate complexes with a covalent C-O bond and a carbon-centered radical (radical intermediates). The radical intermediates exist in two electromers that differ in the oxidation state of iron; Por(+)(*)Fe(III)OCH(2)CH(2)(*) and PorFe(IV)OCH(2)CH(2)(*) (Por = porphyrin). These radical intermediates exist in both the doublet- and quartet spin states. The quartet spin intermediates have substantial barriers for transformation to the quartet spin PorFe(III)-epoxide complex (2.3 kcal mol(-)(1) for PorFe(IV)OCH(2)CH(2)(*) and 7.2 kcal mol(-)(1) for Por(+)(*)Fe(III)OCH(2)CH(2)(*)). In contrast, the doublet spin radicals collapse to the corresponding PorFe(III)-epoxide complex with virtually no barriers. Consequently, the lifetimes of the radical intermediates are much longer on the quartet- than on the doublet spin surface. The loss of isomeric identity in the epoxide and rearrangements to other products arise therefore mostly, if not only, from the quartet process, while the doublet state epoxidation is effectively concerted (Scheme 7). Experimental trends are discussed in the light of the computed mechanistic scheme, and a comparison is made with closely related mechanistic schemes deduced from experiment.     

A theoretical study on the mechanism of camphor hydroxylation by compound I of cytochrome p450

Mechanistic and energetic aspects for the conversion of camphor to 5-exo-hydroxycamphor by the compound I iron-oxo species of cytochrome P450 are discussed from B3LYP DFT calculations. This reaction occurs in a two-step manner along the lines that the oxygen rebound mechanism suggests. The activation energy for the first transition state of the H atom abstraction at the C5 atom of camphor is computed to be more than 20 kcal/mol. This H atom abstraction is the rate-determining step in this hydroxylation reaction, leading to a reaction intermediate that involves a carbon radical species and the iron-hydroxo species. The second transition state of the rebound step that connects the reaction intermediate and the product alcohol complex lies a few kcal/mol below that for the H atom abstraction on the doublet and quartet potential energy surfaces. This energetic feature allows the virtually barrierless recombination in both spin states, being consistent with experimentally observed high stereoselectivity and brief lifetimes of the reaction intermediate. The overall energetic profile of the catalytic mechanism of camphor hydroxylation particularly with respect to why the high activation energy for the H atom abstraction is accessible under physiological conditions is also considered and calculated. According to a proton source model involving Thr252, Asp251, and two solvent water molecules (Biochemistry 1998, 37, 9211), the energetics for the conversion of the iron-peroxo species to compound I is studied. A significant energy over 50 kcal/mol is released in the course of this dioxygen activation process. The energy released in this chemical process is an important driving force in alkane hydroxylation by cytochrome P450. This energy is used for the access to the high activation energy for the H atom abstraction.     

Gauging the relative oxidative powers of compound I, ferric-hydroperoxide, and the ferric-hydrogen peroxide species of cytochrome P450 toward C-H hydroxylation of a radical clock substrate

Density functional calculations were performed in response to the controversies regarding the identity of the oxidant species in cytochrome P450. The calculations were used to gauge the relative C-H hydroxylation reactivity of three potential oxidant species of the enzyme, the high-valent oxo-iron species Compound I (Cpd I), the ferric hydroperoxide Compound 0 (Cpd 0), and the ferric-hydrogen peroxide complex Fe(H(2)O(2)). The results for the hydroxylation of a radical probe substrate, 1, show the following trends: (a) Cpd I is the most reactive species; in its presence the other two reagents will be silent. (b) In the absence of Cpd I, substrate oxidation by Cpd 0 and Fe(H(2)O(2)) will take place via a stepwise mechanism that involves initial O-O homolysis followed by H-abstraction from 1. (c) Cpd 0 will undergo mostly porphyrin hydroxylation and only approximately 15% of substrate oxidation producing mostly the rearranged alcohol, 3 (Scheme 2). (d) Fe(H(2)O(2)) will generate mostly free hydrogen peroxide (uncoupling). A small fraction will perform substrate oxidation and lead mostly to 3. Reactivity probes for these reagents are kinetic isotope effect (KIE) and the product ratio of unrearranged to rearranged alcohols, [2/3]. Thus, for substrate oxidation by Cpd 0 or Fe(H(2)O(2)) KIE will be small, approximately 2, while Cpd I will have large KIE values. Typically both Cpd 0 and Fe(H(2)O(2)) will lead to a [2/3] ratio < 1, while Cpd I will lead to ratios > 1. In addition, the product isotope effect (KIE(2)/KIE(3) not equal 1) is expected from the reactivity of Cpd I.     

Active species of horseradish peroxidase (HRP) and cytochrome P450: two electronic chameleons

The active site of HRP Compound I (Cpd I) is modeled using hybrid density functional theory (UB3LYP). The effects of neighboring amino acids and of environmental polarity are included. The low-lying states have porphyrin radical cationic species (Por(*)(+)). However, since the Por(*)(+) species is a very good electron acceptor, other species, which can be either the ligand or side chain amino acid residues, may participate in electron donation to the Por(*)(+) moiety, thereby making Cpd I behave like a chemical chameleon. Thus, this behavior that was noted before for Cpd I of P450 is apparently much more wide ranging than initially appreciated. Since chemical chameleonic behavior property was found to be expressed not only in the properties of Cpd I itself, but also in its reactivity, the roots of this phenomenon are generalized. A comparative discussion of Cpd I species follows for the enzymes HRP, CcP, APX, CAT (catalase), and P450.     

The directive of the protein: how does cytochrome P450 select the mechanism of dopamine formation?

Dopamine can be generated from tyramine via arene hydroxylation catalyzed by a cytochrome P450 enzyme (CYP2D6). Our quantum mechanical/molecular mechanical (QM/MM) results reveal the decisive impact of the protein in selecting the 'best' reaction mechanism. Instead of the traditional Meisenheimer-complex mechanism, the study reveals a mechanism involving an initial hydrogen atom transfer from the phenolic hydroxyl group of the tyramine to the iron-oxo of the compound I (Cpd I), followed by a ring-π radical rebound that eventually leads to dopamine by keto-enol rearrangement. This mechanism is not viable in the gas phase since the O-H bond activation by Cpd I is endothermic and the process does not form a stable intermediate. By contrast, the in-protein reaction has a low barrier and is exothermic. It is shown that the local electric field of the protein environment serves as a template that stabilizes the intermediate of the H-abstraction step and thereby mediates the catalysis of dopamine formation at a lower energy cost. Furthermore, it is shown that external electric fields can either catalyze or inhibit the process depending on their directionality.     

A theoretical study of reactivity and regioselectivity in the hydroxylation of adamantane by ferrate(VI)

The conversion of adamantane to adamantanols mediated by ferrate (FeO(4)(2)(-)), monoprotonated ferrate (HFeO(4)(-)), and diprotonated ferrate (H(2)FeO(4)) is discussed with the hybrid B3LYP density functional theory (DFT) method. Diprotonated ferrate is the best mediator for the activation of the C-H bonds of adamantane via two reaction pathways, in which 1-adamantanol is formed by the abstraction of a tertiary hydrogen atom (3 degrees ) and 2-adamantanol by the abstraction of a secondary hydrogen atom (2 degrees ). Each reaction pathway is initiated by a C-H bond cleavage via an H-atom abstraction that leads to a radical intermediate, followed by a C-O bond formation via an oxygen rebound step to lead to an adamantanol complex. The activation energies for the C-H cleavage step are 6.9 kcal/mol in the 1-adamantanol pathway and 8.4 kcal/mol in the 2-adamantanol pathway, respectively, at the B3LYP/6-311++G level of theory, whereas those of the second reaction step corresponding to the rebound step are relatively small. Thus, the rate-determining step in the two pathways is the C-H bond dissociation step, which is relevant to the regioselectivity for adamantane hydroxylation. The relative rate constant (3 degrees )/(2 degrees ) for the competing H-atom abstraction reactions is calculated to be 9.30 at 75 degrees C, which is fully consistent with an experimental value of 10.1.     

Reaction mechanism between carbonyl oxide and hydroxyl radical: a theoretical study

The reaction mechanism of carbonyl oxide with hydroxyl radical was investigated by using CASSCF, B3LYP, QCISD, CASPT2, and CCSD(T) theoretical approaches with the 6-311+G(d,p), 6-311+G(2df, 2p), and aug-cc-pVTZ basis sets. This reaction involves the formation of H2CO + HO2 radical in a process that is computed to be exothermic by 57 kcal/mol. However, the reaction mechanism is very complex and begins with the formation of a pre-reactive hydrogen-bonded complex and follows by the addition of HO radical to the carbon atom of H2COO, forming the intermediate peroxy-radical H2C(OO)OH before producing formaldehyde and hydroperoxy radical. Our calculations predict that both the pre-reactive hydrogen-bonded complex and the transition state of the addition process lie energetically below the enthalpy of the separate reactants (DeltaH(298K) = -6.1 and -2.5 kcal/mol, respectively) and the formation of the H2C(OO)OH adduct is exothermic by about 74 kcal/mol. Beyond this addition process, further reaction mechanisms have also been investigated, which involve the abstraction of a hydrogen of carbonyl oxide by HO radical, but the computed activation barriers suggest that they will not contribute to the gas-phase reaction of H2COO + HO.     

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.     

Catalytic mechanism of dopamine beta-monooxygenase mediated by Cu(III)-oxo

Mechanisms of dopamine hydroxylation by the Cu(II)-superoxo species and the Cu(III)-oxo species of dopamine beta-monooxygenase (DBM) are discussed using QM/MM calculations for a whole-enzyme model of 4700 atoms. A calculated activation barrier for the hydrogen-atom abstraction by the Cu(II)-superoxo species is 23.1 kcal/mol, while that of the Cu(III)-oxo, which can be viewed as Cu(II)-O*, is 5.4 kcal/mol. Energies of the optimized radical intermediate in the superoxo- and oxo-mediated pathways are 18.4 and -14.2 kcal/mol, relative to the corresponding reactant complexes, respectively. These results demonstrate that the Cu(III)-oxo species can better mediate dopamine hydroxylation in the protein environment of DBM. The side chains of three amino acid residues (His415, His417, and Met490) coordinate to the Cu(B) atom, one of the copper sites in the catalytic core that plays a role for the catalytic function. The hydrogen-bonding network between dopamine and the three amino acid residues (Glu268, Glu369, and Tyr494) plays an essential role in substrate binding and the stereospecific hydroxylation of dopamine to norepinephrine. The dopamine hydroxylation by the Cu(III)-oxo species is a downhill and lower-barrier process toward the product direction with the aid of the protein environment of DBM. This enzyme is likely to use the high reactivity of the Cu(III)-oxo species to activate the benzylic C-H bond of dopamine; the enzymatic reaction can be explained by the so-called oxygen rebound mechanism.     

Quantum mechanical/molecular mechanical study of mechanisms of heme degradation by the enzyme heme oxygenase: the strategic function of the water cluster

Heme degradation by heme oxygenase (HO) enzymes is important in maintaining iron homeostasis and prevention of oxidative stress, etc. In response to mechanistic uncertainties, we performed quantum mechanical/molecular mechanical investigations of the heme hydroxylation by HO, in the native route and with the oxygen surrogate donor H2O2. It is demonstrated that H2O2 cannot be deprotonated to yield Fe(III)OOH, and hence the surrogate reaction starts from the FeHOOH complex. The calculations show that, when starting from either Fe(III)OOH or Fe(III)HOOH, the fully concerted mechanism involving O-O bond breakage and O-C(meso) bond formation is highly disfavored. The low-energy mechanism involves a nonsynchronous, effectively concerted pathway, in which the active species undergoes first O-O bond homolysis followed by a barrier-free (small with Fe(III)HOOH) hydroxyl radical attack on the meso position of the porphyrin. During the reaction of Fe(III)HOOH, formation of the Por+*FeIV=O species, compound I, competes with heme hydroxylation, thereby reducing the efficiency of the surrogate route. All these conclusions are in accord with experimental findings (Chu, G. C.; Katakura, K.; Zhang, X.; Yoshida, T.; Ikeda-Saito, M. J. Biol. Chem. 1999, 274, 21319). The study highlights the role of the water cluster in the distal pocket in creating "function" for the enzyme; this cluster affects the O-O cleavage and the O-Cmeso formation, but more so it is responsible for the orientation of the hydroxyl radical and for the observed alpha-meso regioselectivity of hydroxylation (Ortiz de Montellano, P. R. Acc. Chem. Res. 1998, 31, 543). Differences/similarities with P450 and HRP are discussed.     

Two states and two more in the mechanisms of hydroxylation and epoxidation by cytochrome P450

Past studies have shown that oxidation reactions by P450 Compound I (Cpd I) can be described by two competing quartet and doublet spin states, which possess three unpaired electrons, hence tri-radicals. One electron excitation from the delta orbital to sigma* xy generates two states that possess five unpaired electrons, so-called penta-radicals, in sextet and quartet situations, and which were shown by theory to lie only approximately 12-14 kcal/mol higher in energy than the tri-radical ground states (ref 7). The present study focuses on the C-H hydroxylation and C=C epoxidation of propene by these penta-radical states. It is shown that the initial energy differences, between the penta-radical and tri-radical states, diminish along the reaction pathway, due to the favorable and cumulative exchange stabilization of the more open-shell species. Furthermore, theory suggests that hydrogen bonding to the thiolate ligand, and general polarity of the environment, reduce these gaps further, thereby making the penta-radical states accessible to ground-state reactivity. The interconversion between the tri-radical and penta-radical states along the reaction coordinate will depend on the dynamics of spin-flips and energy barriers between the states. Especially interesting should be the region of the reaction intermediates; for both epoxidation and hydroxylation, this region is typified by a dense manifold of spin states and electromeric states (that differ by the oxidation state of iron), such that the total reactivity would be expected to reflect the interplay of these states, giving rise to multistate reactivity.     

A density functional study of O-O bond cleavage for a biomimetic non-heme iron complex demonstrating an Fe(v)-intermediate

Density functional theory using the B3LYP hybrid functional has been employed to investigate the reactivity of Fe(TPA) complexes (TPA = tris(2-pyridylmethyl)amine), which are known to catalyze stereospecific hydrocarbon oxidation when H(2)O(2) is used as oxidant. The reaction pathway leading to O-O bond heterolysis in the active catalytic species Fe(III)(TPA)-OOH has been explored, and it is shown that a high-valent iron-oxo intermediate is formed, where an Fe(V) oxidation state is attained, in agreement with previous suggestions based on experiments. In contrast to the analogous intermediate [(Por.)Fe(IV)=O](+1) in P450, the TPA ligand is not oxidized, and the electrons are extracted almost exclusively from the mononuclear iron center. The corresponding homolytic O-O bond cleavage, yielding the two oxidants Fe(IV)=O and the OH. radical, has also been considered, and it is shown that this pathway is inaccessible in the hydrocarbon oxidation reaction with Fe(TPA) and hydrogen peroxide. Investigations have also been performed for the O-O cleavage in the Fe(III)(TPA)-alkylperoxide species. In this case, the barrier for O-O homolysis is found to be slightly lower, leading to loss of stereospecificity and supporting the experimental conclusion that this is the preferred pathway for alkylperoxide oxidants. The difference between hydroperoxide and alkylperoxide as oxidant derives from the higher O-O bond strength for hydrogen peroxide (by 8.0 kcal/mol).