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.     

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.     

Proton-shuffle mechanism of O-O activation for formation of a high-valent oxo-iron species of bleomycin

Bleomycins (BLMs) can utilize H2O2 to cleave DNA in the presence of ferric ions. DFT calculations were used to study the mechanism of O-O bond cleavage in the low-spin FeIII-hydroperoxo complex of BLM. The following alternative hypotheses were investigated using realistic structural models: (a) heterolytic cleavage of the O-O bond, generating a Compound I (Cpd I) like intermediate, formally BLM-FeV=O; (b) homolytic O-O cleavage, leading to a BLM-FeIV=O species and an OH* radical; and (c) a direct O-O cleavage/H-abstraction mechanism by ABLM. The calculations showed that (a) is a facile and viable mechanism; it involves acid-base proton reshuffle mediated by the side-chain linkers of BLM, causing thereby heterolytic cleavage of the O-O bond and generation of Cpd I. Formation of Cpd I is found to involve a barrier of 13.3 kcal/mol, which is lower than the barriers in the alternative mechanisms (b and c) that possess respective barriers of 31 and 17 kcal/mol. The so-formed Cpd I species with a radical on the side-chain linker, methylvalerate (V), adjacent to the BLM-FeIV=O complex, resembles the formation of the active species of cytochrome c peroxidase in the Poulos-Kraut proton-shuffle mechanism in heme peroxidases (Poulos, T. L.; Kraut, J. J. Biol. Chem. 1980, 255, 8199-8205). Experimental data are discussed and shown to be in accord with this proposal. It suggests that the high-valence Cpd I species of BLM participates in the DNA cleavage. This is an alternative mechanistic hypothesis to the exclusive reactivity scenario based on ABLM (FeIII-OOH).     

How Is a Metabolic Intermediate Formed in the Mechanism‐Based Inactivation of Cytochrome P450 by Using 1,1‐Dimethylhydrazine: Hydrogen Abstraction or Nitrogen Oxidation?

A precise understanding of the mechanism‐based inactivation of cytochrome P450 enzymes (P450s) at the quantum mechanical level should allow more reliable predictions of drug–drug interactions than those currently available. Hydrazines are among the molecules that act as mechanism‐based inactivators to terminate the function of P450s, which are essential heme enzymes responsible for drug metabolism in the human body. Despite its importance, the mechanism explaining how a metabolic intermediate (MI) is formed from hydrazine is not fully understood. We used density functional theory (DFT) calculations to compare four possible mechanisms underlying the reaction between 1,1‐dimethylhydrazine (or unsymmetrical dimethylhydrazine, UDMH) and the reactive compound I (Cpd I) intermediate of P450. Our DFT calculations provided a clear view on how an aminonitrene‐type MI is formed from UDMH. In the most favorable pathway, hydrogen is spontaneously abstracted from the N2 atom of UDMH by Cpd I, followed by a second hydrogen abstraction from the N2 atom by Cpd II. Nitrogen oxidation of nitrogen atoms and hydrogen abstraction from the C? H bond of the methyl group were found to be less favorable than the hydrogen abstraction from the N? H bond. We also found that the reaction of protonated UDMH with Cpd I is rather sluggish. The aminonitrene‐type MI binds to the ferric heme more strongly than a water molecule. This is consistent with the notion that the catalytic cycle of P450 is impeded when such an MI is produced through the P450‐catalyzed reaction.     

Cytochrome P450 compound I in the plane wave pseudopotential framework: GGA electronic and geometric structure of thiolate‐ligated iron(IV)–oxo porphyrin

The cytochromes P450 constitute a ubiquitous family of metalloenzymes, catalyzing manifold reactions of biological and synthetic importance via a thiolate‐ligated iron‐oxo (IV) porphyrin radical species denoted compound I (Cpd I). Experimental investigations have implicated this intermediate in a broad spectrum of biophysically interesting phenomena, further augmenting the importance of a Cpd I model system. Ab initio molecular dynamics, including Car–Parrinello and path integral methods, conjoin electronic structure theory with finite temperature simulation, affording tools most valuable to approach such enzymes. These methods are typically driven by density functional theory (DFT) in a plane‐wave pseudopotential framework; however, existing studies of Cpd I have been restricted to localized Gaussian basis sets. The appropriate choice of density functional and pseudopotential for such simulations is accordingly not obvious. To remedy this situation, a systematic benchmarking of thiolate‐ligated Cpd I is performed using several generalized‐gradient approximation (GGA) functionals in the Martins–Troullier and Vanderbilt ultrasoft pseudopotential schemes. The resultant electronic and structural parameters are compared to localized–basis DFT calculations using GGA and hybrid density functionals. The merits and demerits of each scheme are presented in the context of reproducing existing experimental and theoretical results for Cpd I. © 2013 Wiley Periodicals, Inc.     

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.     

Versatility of the electronic structure of compound I in catalase-peroxidases

Catalase-peroxidases (KatGs) are bifunctional heme proteins, belonging to the family of class I peroxidases, that are able to catalyze both catalatic and peroxidatic reactions within a peroxidase-like structure. We investigated the electronic structure of reaction intermediates of the catalytic cycle of KatGs by means of density functional theory (DFT) QM/MM calculations. The outcome was that the ionization state of the KatG-specific covalent adduct (Met264-Tyr238-Trp111) affects the radical character of compound I (Cpd I). Specifically, in the optimized structures, substantial radical character is observed on the proximal Trp330 when Tyr238 is protonated, whereas when Tyr238 is deprotonated the radical localizes on the Met+-Tyr(O-)-Trp adduct. These findings are not affected by protein thermal fluctuations, although details of the spin density distribution are affected by the geometry of the active site. Calculations provide structures in good agreement with the crystal structure of BpKatG Cpd I. They also provide an explanation for the experimental findings of the mobile and catalatic-specific residue Arg426 being 100% in conformation R in the X-ray structure of BpKatG treated with organic peroxides. The role of different Cpd I forms in the catalase and peroxidase reaction pathways is discussed.     

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.     

DFT Calculations Suggest a New Type of Self‐Protection and Self‐Inhibition Mechanism in the Mammalian Heme Enzyme Myeloperoxidase: Nucleophilic Addition of a Functional Water rather than One‐Electron Reduction

The mammalian heme enzyme myeloperoxidase (MPO) catalyzes the reaction of Cl^? to the antimicrobial‐effective molecule HOCl. During the catalytic cycle, a reactive intermediate “Compound I” (Cpd I) is generated. Cpd I has the ability to destroy the enzyme. Indeed, in the absence of any substrate, Cpd I decays with a half‐life of 100 ms to an intermediate called Compound II (Cpd II), which is typically the one‐electron reduced Cpd I. However, the nature of Cpd II, its spectroscopic properties, and the source of the additional electron are only poorly understood. On the basis of DFT and time‐dependent (TD)‐DFT quantum chemical calculations at the PBE0/6‐31G* level, we propose an extended mechanism involving a new intermediate, which allows MPO to protect itself from self‐oxidation or self‐destruction during the catalytic cycle. Because of its similarity in electronic structure to Cpd II, we named this intermediate Cpd II′. However, the suggested mechanism and our proposed functional structure of Cpd II′ are based on the hypothesis that the heme is reduced by charge separation caused by reaction with a water molecule, and not, as is normally assumed, by the transfer of an electron. In the course of this investigation, we found a second intermediate, the reduced enzyme, towards which the new mechanism is equally transferable. In analogy to Cpd II′, we named it Fe^II′. The proposed new intermediates Cpd II′ and Fe^II′ allow the experimental findings, which have been well documented in the literature for decades but not so far understood, to be explained for the first time. These encompass a) the spontaneous decay of Cpd I, b) the unusual (chlorin‐like) UV/Vis, circular dichroism (CD), and resonance Raman spectra, c) the inability of reduced MPO to bind CO, d) the fact that MPO‐Cpd II reduces SCN^? but not Cl^?, and e) the experimentally observed auto‐oxidation/auto‐reduction features of the enzyme. Our new mechanism is also transferable to cytochromes, and could well be viable for heme enzymes in general.     

DFT Calculations Suggest a New Type of Self-Protection and Self-Inhibition Mechanism in the Mammalian Heme Enzyme Myeloperoxidase: Nucleophilic Addition of a Functional Water rather than One-Electron Reduction

The mammalian heme enzyme myeloperoxidase (MPO) catalyzes the reaction of Cl(-) to the antimicrobial-effective molecule HOCl. During the catalytic cycle, a reactive intermediate "Compound?I" (Cpd?I) is generated. Cpd?I has the ability to destroy the enzyme. Indeed, in the absence of any substrate, Cpd?I decays with a half-life of 100?ms to an intermediate called Compound?II (Cpd?II), which is typically the one-electron reduced Cpd?I. However, the nature of Cpd?II, its spectroscopic properties, and the source of the additional electron are only poorly understood. On the basis of DFT and time-dependent (TD)-DFT quantum chemical calculations at the PBE0/6-31G* level, we propose an extended mechanism involving a new intermediate, which allows MPO to protect itself from self-oxidation or self-destruction during the catalytic cycle. Because of its similarity in electronic structure to Cpd?II, we named this intermediate Cpd?II'. However, the suggested mechanism and our proposed functional structure of Cpd?II' are based on the hypothesis that the heme is reduced by charge separation caused by reaction with a water molecule, and not, as is normally assumed, by the transfer of an electron. In the course of this investigation, we found a second intermediate, the reduced enzyme, towards which the new mechanism is equally transferable. In analogy to Cpd?II', we named it Fe(II') . The proposed new intermediates Cpd?II' and Fe(II') allow the experimental findings, which have been well documented in the literature for decades but not so far understood, to be explained for the first time. These encompass a)?the spontaneous decay of Cpd?I, b)?the unusual (chlorin-like) UV/Vis, circular dichroism (CD), and resonance Raman spectra, c)?the inability of reduced MPO to bind CO, d)?the fact that MPO-Cpd?II reduces SCN(-) but not Cl(-) , and e)?the experimentally observed auto-oxidation/auto-reduction features of the enzyme. Our new mechanism is also transferable to cytochromes, and could well be viable for heme enzymes in general.     

Compound I in heme thiolate enzymes: a comparative QM/MM study

This study directly compares the active species of heme enzymes, so-called Compound I (Cpd I), across the heme-thiolate enzyme family. Thus, sixty-four different Cpd I structures are calculated by hybrid quantum mechanical/molecular mechanical (QM/MM) methods using four different cysteine-ligated heme enzymes (P450(cam), the mutant P450(cam)-L358P, CPO and NOS) with varying QM region sizes in two multiplicities each. The overall result is that these Cpd I species are similar to each other with regard to many characteristic features. Hence, using the more stable CPO Cpd I as a model for P450 Cpd I in experiments should be a reasonable approach. However, systematic differences were also observed, and it is shown that NOS stands out in most comparisons. By analyzing the electrical field generated by the enzyme on the QM region, one can see that (a) the protein exerts a large influence and modifies all the Cpd I species compared with the gas-phase situation and (b) in NOS this field is approximately planar to the heme plane, whereas it is approximately perpendicular in the other enzymes, explaining the deviating results on NOS. The calculations on the P450(cam) mutant L358P show that the effects of removing the hydrogen bond between the heme sulfur and L358 are small at the Cpd I stage. Finally, Mossbauer parameters are calculated for the different Cpd I species, enabling future comparisons with experiments. These results are discussed in the broader context of recent findings of Cpd I species that exhibit large variations in the electronic structure due to the presence of the substrate.     

Compound I of nitric oxide synthase: the active site protonation state

A quantum mechanical/molecular mechanical (QM/MM) study of the formation of the elusive active species Compound I (Cpd I) of nitric oxide synthase (NOS) from the oxyferrous intermediate shows that two protons have to be provided to produce a reaction that is reasonably exothermic and that leads to the appearance of a radical on the tetrahydrobiopterin cofactor. Molecular dynamics and energy considerations show that a possible source of proton is the water H-bond chain formed from the surface to the active site, but that a water molecule by itself cannot be the source of the proton; an H3O+ species that is propagated along the chain is more likely. The QM/MM calculations demonstrate that Cpd I and H2O are formed from the ferric-hydrogen peroxide complex in a unique heterolytic O-O cleavage mechanism. The properties of the so-formed Cpd I are compared with those of the known species of chloroperoxidase, and the geometry and spin densities are found to be compatible. The M?ssbauer parameters are calculated and may serve as experimental probes in attempts to characterize NOS Cpd I.     

Electronic state of the molecular oxygen released by catalase

In catalases, the high redox intermediate known as compound I (Cpd I) is reduced back to the resting state by means of hydrogen peroxide in a 2-electron reaction [Cpd I (Por(*+)-Fe(IV)O) + H(2)O(2) --> Enz (Por-Fe(III)) + H(2)O + O(2)]. It has been proposed that this reaction takes place via proton transfer toward the distal His and hydride transfer toward the oxoferryl oxygen (H(+)/H(-) scheme) and some authors have related it to singlet oxygen generation. Here, we consider the possible reaction schemes and qualitatively analyze the electronic state of the species involved to show that the commonly used association of the H(+)/H(-) scheme with singlet oxygen production is not justified. The analysis is complemented with density functional theory (DFT) calculations for a gas-phase active site model of the reactants and products.     

Spectroelectrochemical and computational studies on the mechanism of hypoxia selectivity of copper radiopharmaceuticals

Detailed chemical, spectroelectrochemical and computational studies have been used to investigate the mechanism of hypoxia selectivity of a range of copper radiopharmaceuticals. A revised mechanism involving a delicate balance between cellular uptake, intracellular reduction, reoxidation, protonation and ligand dissociation is proposed. This mechanism accounts for observed differences in the reported cellular uptake and washout of related copper bis(thiosemicarbazonato) complexes. Three copper and zinc complexes have been characterised by X-ray crystallography and the redox chemistry of a series of copper complexes has been investigated by using electronic absorption and EPR spectroelectrochemistry. Time-dependent density functional theory (TD-DFT) calculations have also been used to probe the electronic structures of intermediate species and assign the electronic absorption spectra. DFT calculations also show that one-electron oxidation is ligand-based, leading to the formation of cationic triplet species. In the absence of protons, metal-centred one-electron reduction gives the reduced anionic copper(I) species, [CuIATSM](-), and for the first time it is shown that molecular oxygen can reoxidise this anion to give the neutral, lipophilic parent complexes, which can wash out of cells. The electrochemistry is pH dependent and in the presence of stronger acids both chemical and electrochemical reduction leads to quantitative and rapid dissociation of copper(I) ions from the mono- or diprotonated complexes, [CuIATSMH] and [Cu(I)ATSMH2]+. In addition, a range of protonated intermediate species have been identified at lower acid concentrations. The one-electron reduction potential, rate of reoxidation of the copper(I) anionic species and ease of protonation are dependent on the structure of the ligand, which also governs their observed behaviour in vivo.     

Catalytic mechanism of cytochrome P450 for 5'-hydroxylation of nicotine: fundamental reaction pathways and stereoselectivity

A series of computational methods were used to study how cytochrome P450 2A6 (CYP2A6) interacts with (S)-(-)-nicotine, demonstrating that the dominant molecular species of (S)-(-)-nicotine in CYP2A6 active site exists in the free base state (with two conformations, SR(t) and SR(c)), despite the fact that the protonated state is dominant for the free ligand in solution. The computational results reveal that the dominant pathway of nicotine metabolism in CYP2A6 is through nicotine free base oxidation. Further, first-principles quantum mechanical/molecular mechanical free energy (QM/MM-FE) calculations were carried out to uncover the detailed reaction pathways for the CYP2A6-catalyzed nicotine 5'-hydroxylation reaction. In the determined CYP2A6-(S)-(-)-nicotine binding structures, the oxygen of Compound I (Cpd I) can abstract a hydrogen from either the trans-5'- or the cis-5'-position of (S)-(-)-nicotine. CYP2A6-catalyzed (S)-(-)-nicotine 5'-hydroxylation consists of two reaction steps, that is, the hydrogen transfer from the 5'-position of (S)-(-)-nicotine to the oxygen of Cpd I (the H-transfer step), followed by the recombination of the (S)-(-)-nicotine moiety with the iron-bound hydroxyl group to generate the 5'-hydroxynicotine product (the O-rebound step). The H-transfer step is rate-determining. The 5'-hydroxylation proceeds mainly with the stereoselective loss of the trans-5'-hydrogen, that is, the 5'-hydrogen trans to the pyridine ring. The calculated overall stereoselectivity of ～97% favoring the trans-5'-hydroxylation is close to the observed stereoselectivity of 89-94%. This is the first time it has been demonstrated that a CYP substrate exists dominantly in one protonation state (cationic species) in solution, but uses its less-favorable protonation state (neutral free base) to perform the enzymatic reaction.     

Combined Experimental and Theoretical Study on the Reactivity of Compounds I and II in Horseradish Peroxidase Biomimetics

For the exploration of the intrinsic reactivity of two key active species in the catalytic cycle of horseradish peroxidase (HRP), Compound I (HRP‐I) and Compound II (HRP‐II), we generated in situ [Fe^IV?O(TMP^+.)(2‐MeIm)]⁺ and [Fe^IV?O(TMP)(2‐MeIm)]⁰ (TMP=5,10,15,20‐tetramesitylporphyrin; 2‐MeIm=2‐methylimidazole) as biomimetics for HRP‐I and HRP‐II, respectively. Their catalytic activities in epoxidation, hydrogen abstraction, and heteroatom oxidation reactions were studied in acetonitrile at ?15 °C by utilizing rapid‐scan UV/Vis spectroscopy. Comparison of the second‐order rate constants measured for the direct reactions of the HRP‐I and HRP‐II mimics with the selected substrates clearly confirmed the outstanding oxidizing capability of the HRP‐I mimic, which is significantly higher than that of HRP‐II. The experimental study was supported by computational modeling (DFT calculations) of the oxidation mechanism of the selected substrates with the involvement of quartet and doublet HRP‐I mimics (^2,4Cpd I) and the closed‐shell triplet spin HRP‐II model (³Cpd II) as oxidizing species. The significantly lower activation barriers calculated for the oxidation systems involving ^2,4Cpd I than those found for ³Cpd II are in line with the much higher oxidizing efficiency of the HRP‐I mimic proven in the experimental part of the study. In addition, the DFT calculations show that all three reaction types catalyzed by HRP‐I occur on the doublet spin surface in an effectively concerted manner, whereas these reactions may proceed in a stepwise mechanism with the HRP‐II mimic as oxidant. However, the high desaturation or oxygen rebound barriers during C?H bond activation processes by the HRP‐II mimic predict a sufficient lifetime for the substrate radical formed through hydrogen abstraction. Thus, the theoretical calculations suggest that the dissociation of the substrate radical may be a more favorable pathway than desaturation or oxygen rebound processes. Importantly, depending on the electronic nature of the oxidizing species, that is, ^2,4Cpd I or ³Cpd II, an interesting region‐selective conversion phenomenon between sulfoxidation and H‐atom abstraction was revealed in the course of the oxidation reaction of dimethylsulfide. The combined experimental and theoretical study on the elucidation of the intrinsic reactivity patterns of the HRP‐I and HRP‐II mimics provides a valuable tool for evaluating the particular role of the HRP active species in biological systems.     

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.     

On the functional role of a water molecule in clade 3 catalases: a proposal for the mechanism by which NADPH prevents the formation of compound II

X-ray structures of the 13 different monofunctional heme catalases published to date were scrutinized in order to gain insight in the mechanism by which NADPH in Clade 3 catalases may protect the reactive ferryloxo intermediate Compound I (Cpd I; por (*+)Fe (IV)O) against deactivation to the catalytically inactive intermediate Compound II (Cpd II; porFe (IV)O). Striking similarities in the molecular network of the protein subunits encompassing the heme center and the surface-bound NADPH were found for all of the Clade 3 catalases. Unique features in this region are the presence of a water molecule (W1) adjacent to the 4-vinyl group of heme and a serine residue or a second water molecule hydrogen-bonded to both W1 and the carbonyl group of a threonine-proline linkage, with the proline in van der Waals contact with the dihydronicotinamide group of NADPH. A mechanism is proposed in which a hydroxyl anion released from W1 undergoes reversible nucleophilic addition to the terminal carbon of the 4-vinyl group of Cpd I, thereby producing a neutral porphyrin pi-radical ferryloxo (HO-por (*)Fe (IV)O) species of reduced reactivity. This structure is suggested to be the elusive Cpd II' intermediate proposed in previous studies. An accompanying proton-shifting process along the hydrogen-bonded network is believed to facilitate the NADPH-mediated reduction of Cpd I to ferricatalase and to serve as a funnel for electron transfer from NADPH to the heme center to restore the catalase Fe (III) resting state. The proposed reaction paths were fully supported as chemically reasonable and energetically feasible by means of density functional theory calculations at the (U)B3LYP/6-31G* level. A particularly attractive feature of the present mechanism is that the previously discussed formation of protein-derived radicals is avoided.     

Quantum and molecular mechanical study of the first proton transfer in the catalytic cycle of cytochrome P450cam and its mutant D251N

In the catalytic cycle of cytochrome P450cam, the hydroperoxo intermediate (Cpd 0) is formed by proton transfer from a reduced oxyheme complex (S5). This process is drastically slowed down when Asp251 is mutated to Asn (D251N). We report quantum mechanical/molecular mechanical (QM/MM) calculations that address this proton delivery in the doublet state through a hydrogen-bond network in the Asp251 channel, both for the wild-type enzyme and the D251N mutant, using four different active-site models. For the wild-type, we find a facile concerted mechanism for proton transfer from protonated Asp251 via Wat901 and Thr252 to the FeOO moiety, with a barrier of about 1 kcal/mol and a high exothermicity of more than 20 kcal/mol. In the D251N mutant with a neutral Asn251 residue, the proton transfer is almost thermoneutral or slightly exothermic in the three models considered. It is still very facile when the Asn251 residue adopts a conformation analogous to Asp251 in the wild-type enzyme, but the barrier increases significantly when the Asn251 side chain flips (as indicated by classical molecular dynamics simulations). This flip disrupts the hydrogen-bond network and hence the proton-transfer pathway, which causes a longer lifetime of S5 in the D251N mutant (consistent with experimental observations). The entry of an additional water molecule into the active site of D251N with flipped Asn251 regenerates the hydrogen-bond network and provides a viable mechanism for proton delivery in the mutant, with a moderate barrier of about 7 kcal/mol.