A 175-psec molecular dynamics simulation of camphor-bound cytochrome P-450cam.

The structure and internal motions of cytochrome P-450cam, a monooxygenase heme enzyme with 414 amino acid residues, with camphor bound at the active site have been evaluated on the basis of a 175-psec molecular dynamics simulation carried out at 300 K. All hydrogen atoms were explicitly modeled, and 204 crystallographic waters were included in the simulation. Based on an analysis of the time course of the trajectory versus potential energy, root mean square deviation, radius of gyration, and hydrogen bonding, the simulation was judged to be stable and representative of the average experimental structure. The averaged structural properties of the enzyme were evaluated from the final 135 psec of the simulation. The average atomic displacement from the X-ray structure was 1.39 A for all heavy atoms and 1.17 A for just C-alpha atoms. The average root-mean-square (rms) fluctuations of all heavy atoms and backbone atoms were 0.42 and 0.37 A, respectively. The computed rms fluctuations were in reasonable agreement with the experimentally determined temperature factors. All 13 segments of alpha-helix and 5 segments of beta-sheet were well preserved with the exception of the N-terminal half of helix F which alternated between an alpha-helix and a 310-helix. In addition there were in general only small variations in the relative orientation of adjacent alpha-helices. The rms fluctuations of the backbone dihedral angles in the secondary structure elements were almost uniformly smaller, with the fluctuation in alpha-helices and beta-sheets, 31 and 10% less, respectively, than those in nonsecondary structure regions. The reported crystal structure contains kinks in both helices C and I. In the simulation, both of these regions showed high mobility and large deviations from their starting positions. Since the kink in the I helix is at the oxygen binding site, these motions may have mechanistic implications.     

Molecular Dynamics Simulations Indicate that F87W,T185F-Cytochrome P450cam May Reductively Dehalogenate 1,1,1-Trichloroethane

Abstract

Cytochrome P450cam is capable of reductively dehalogenating several chlorinated alkanes at low, but measurable, rates. In previous investigations of structure-function relationships in this enzyme using molecular dynamics simulations, we noticed that 1,1,1-trichloroethane (TCA) exhibits a very high degree of mobility in the active site due to its smaller molecular volume relative to the native substrate, camphor(1,2). Several amino acid sidechains lining the active site also exhibit significant dynamic fluctuations, possibly as a result of poor steric complementarity to TCA. Guided by these results, we modeled double (F87W, T185F) and triple (F87W, T185F, V295I) mutants of P450cam, which provide additional bulk in the active site and increase the frequency of heme-substrate collision. Molecular dynamics simulations (300 ps on each protein) indicate that these mutants do not significantly perturb the three-dimensional fold of the enzyme, or local structure in the region of the active site. Both mutants bind the substrate more stably near the heme than the wild-type. Interestingly, however, the bulkier triple mutant seems to actually inhibit heme-substrate interactions relative to the double mutant. Over the final 200 ps of simulation, TCA is within 1 Å of nonbonded contact with the heme 25% more often in the double mutant versus the wild-type. The triple mutant, on the other hand, binds TCA within 1 Å of the heme only 15% as often as the wild-type. These results indicate that the double mutant may reductively dehalogenate TCA, a property not observed for the native protein. Implications for other experimentally measurable parameters are discussed.     

Molecular dynamics simulations of phenylimidazole inhibitor complexes of cytochrome P450 cam

Molecular dynamics simulations have been performed on three phenylimidazole inhibitor complexes ofP450 _cam, utilizing the X-ray structures and the AMBER suite of programs. Compared to their corresponding optimized X-ray structures, very similar features were observed for the 1-phenylimidazole (1-PI) and 2-phenylimidazole (2-PI) complexes during a 100 ps MD simulation. The 1-PI inhibitor binds as a Type II complex with the imidazole nitrogen as a ligand of the heme iron. Analysis of the inhibitor-enzyme interctions during the MD simulations reveals that electrostatic interactions of the imidazole with the heme and van der Waals interactions of the phenyl ring with nearby hydrophobic residues are dominant. By contrast, 2-PI binds as a Type I inhibitor in the substrate binding pocket, but not as a ligand of the iron. The interactions of this inhibitor are qualitatively different from that of the Type II 1-PI, being mainly electrostatic/H-bonding interactions with a bound water and polar residues. Although the third compound, 4-PI, in common with 1-PI, also binds as a Type II inhibitor, with one nitrogen of the imidazole as a ligand to the iron, the MD average binding orientation deviates significantly from the X-ray structure. The most important changes observed include: (1) the rotation of the imidazole ring of this inhibitor by about 90° to enhance electrostatic interactions of the imidazole NH group with the carbonyl group of LEU244, and (2) the rotation of the carbonyl group of ASP251 to form a H-bond with VAL254. An analysis of the H-bonding network surrounding this substrate in the optimized crystal structure revealed that there is no H-bonding partner either for the free polar NH group in the imidazole ring of 4-phenylimidazole or for the polar carbonyl group of the nearby ASP251 residue. The deviation of the dynamically averaged inhibitor-enzyme structure of the 4-PI complex from the optimized crystal structure can therefore be rationalized as a consequence of the optimization of the electrostatic interactions among the polar groups.     

Zinc-substituted cytochrome P450cam: characterization of protein conformers F420 and F450 by photoinduced electron transfer

Metal substitution of heme proteins is widely applied in the study of biologically relevant electron transfer (ET) reactions. It has been shown that many modified proteins remain in their native conformation and can provide useful insights into the molecular mechanism of electron transfer between the native protein and its substrates. We investigated ET reactions between zinc-substituted cytochrome P450(cam) and small organic compounds such as quinones and ferrocene, which are capable of accessing the protein's hydrophobic channel and binding close to the active site, like its native substrate, camphor. Following the substitution method developed by Gunsalus and co-workers [Wagner, G. C., et al. (1981) J. Biol. Chem. 256, 6262-6265], we have identified two dominant forms of the zinc-substituted protein, F450 and F420, that exhibit different photophysical and photochemical properties. The ET behavior of F420 suggests that hydrophobic redox-active ligands are able to penetrate the hydrophobic channel and place themselves in the direct vicinity of the Zn-porphyrin. In contrast, the slower ET quenching rates observed in the case of F450 indicate that the association is weak and occurs outside of the protein channel. Therefore, we conclude that F420 corresponds to the open structure of the native cytochrome P450(cam) while F450 has a closed or partially closed channel that is characteristic of the camphor-containing cytochrome P450(cam). The existence of two distinct conformers of Zn-bound P450(cam) is consistent with the findings of Goodin and co-workers [Lee, Y.-T., et al. (2010) Biochemistry 49, 3412-3419] and has significant consequences for future electron transfer studies on this popular metalloenzyme.     

Surface-modified mutants of cytochrome P450cam: enzymatic properties and electrochemistry.

We report the electrochemistry of genetic variants of the haem monooxygenase cytochrome P450cam. A surface cysteine-free mutant (abbreviated as SCF) was prepared in which the five surface cysteine residues Cys-58, Cys-85, Cys-136, Cys-148 and Cys-334 were changed to alanines. Four single surface cysteine mutants with an additional mutation, R72C, R112C, K344C or R364C, were also prepared. The haem spin-state equilibria, NADH turnover rates and camphor-hydroxylation properties, as well as the electrochemistry of these mutants are reported. The coupling of a redox-active label, N-ferrocenylmaleimide, to the single surface cysteine mutant SCF-K344C, and the electrochemistry of this modified mutant are also described.     

Mechanism of the plant cytochrome P450 for herbicide resistance: a modelling study

Plant cytochrome P450 is a key enzyme responsible for the herbicide resistance but the molecular basis of the mechanism is unclear. To understand this, four typical plant P450s and a widely resistant herbicide chlortoluron were analysed by carrying out homology modelling, molecular docking, molecular dynamics simulations and binding free energy analysis. Our results demonstrate that: (i) the putative hydrophobic residues located in the F-helix and polar residues in I-helix are critical in the herbicide resistance; (ii) the binding mode analysis and binding free energy calculation indicate that the distance between catalytic site of chlortoluron and heme of P450, as well as the binding affinity are key elements affecting the resistance for plants. In conclusion, this work provides a new insight into the interactions of plant P450s with herbicide from a molecular level, offering valuable information for the future design of novel effective herbicides which also escape from the P450 metabolism.     

Structural alterations of the heme environment of cytochrome P450cam and the Y96F mutant as deduced by resonance Raman spectroscopy

Resonance Raman spectroscopy at 2.5cm(-1) resolution was used to probe differences in wild-type and Y96F mutant P450cam (CYP101), both with and without bound camphor or styrene substrates. In the substrate-free state, the spin state equilibrium is shifted from 6-coordinate low spin (6CLS) toward more 5-coordinate high spin (5CHS) when tyrosine-96 in the substrate pocket is replaced by phenylalanine. About 25% of substrate-free Y96F mutant is 5CHS as opposed to 8% for substrate-free wild-type P450cam. Spin equilibrium constants calculated from Raman intensities indicate that the driving force for electron transfer from putidaredoxin, the natural redox partner of P450cam, is significantly smaller on styrene binding than for camphor binding. Spectral differences suggest that there is a tilt in camphor toward the pyrrole III ring on Y96F mutation. This finding is consistent with the altered product distribution found for camphor hydroxylation by the Y96F mutant relative to the single enantiomer produced by the wild-type enzyme.     

Crystallographic study on the dioxygen complex of wild-type and mutant cytochrome P450cam. Implications for the dioxygen activation mechanism

Two key amino acids, Thr252 and Asp251, are known to be important for dioxygen activation by cytochrome P450cam. We have solved crystal structures of a critical intermediate, the ferrous dioxygen complex (Fe(II)-O2), of the wild-type P450cam and its mutants, D251N and T252A. The wild-type dioxygen complex structure is very much the same as reported previously (Schlichting, I., Berendzen, J., Chu, K., Stock, A. M., Maves, S. A., Benson, D. E., Sweet, R. M., Ringe, D., Petsko, G. A., and Sligar, S. G. (2000) Science 287, 1615-1622) with the exception of higher occupancy and a more ordered structure of the iron-linked dioxygen and two "catalytic" water molecules that form part of a proton relay system to the iron-linked dioxygen. Due to of the altered conformation of the I helix groove these two waters are missing in the D251N dioxygen complex which explains its lower catalytic activity and slower proton transfer to the dioxygen ligand. Similarly, the T252A mutation was expected to disrupt the active site solvent structure leading to hydrogen peroxide formation rather than substrate hydroxylation. Unexpectedly, however, the two "catalytic" waters are retained in the T252A mutant. Based on these findings, we propose that the Thr(252) accepts a hydrogen bond from the hydroperoxy (Fe(III)-OOH) intermediate that promotes the second protonation on the distal oxygen atom, leading to O-O bond cleavage and compound I formation.     

Leu254 residue and calcium ions as new structural determinants of carboxypeptidase T substrate specificity

New determinants of Thermoactinomyces vulgaris carboxypeptidase T (CPT) substrate specificity--structural calcium ions and Leu254 residue--were found by means of steady-state kinetics and site-directed mutagenesis. The removal of calcium ions shifted the selectivity profile of hydrolysis of tripeptide substrates with C-terminal Leu, Glu, and Arg from 64/1.7/1 to 162/1.3/1. Substitution of the hydrophobic Leu254 in CPT for polar Asn did not change hydrolysis efficiency of substrates with C-terminal Leu and Arg, but resulted in more than 28-fold decrease in activity towards the substrate with C-terminal Glu. It is shown that the His68 residue is not a structural determinant of CPT specificity.     

Smoluchowski dynamics of the vnd/NK‐2 homeodomain from Drosophila melanogaster: First‐order mode‐coupling approximation

This work is the first in a series devoted to applying mode coupling diffusion theory to the derivation of local dynamics properties of proteins in solution. The first‐order mode‐coupling approximation, or optimized Rouse–Zimm local dynamics (ORZLD), is applied here to derive the rotational dynamics of the bonds and compare the calculated with the experimental nmr ¹⁵N spin–lattice relaxation time behavior of the vnd/NK‐2 homeodomain from Drosophila melanogaster. The starting point for the calculations is the experimental three‐dimensional structure of the homeodomain determined by multidimensional nmr spectroscopy. The results of the computations are compared with experimentally measured ¹⁵N spin–lattice relaxation times T₁, at 34.5 and 60.8 MHz, to check the first‐order approximation. To estimate the relative importance of internal and overall rotation, both rigid and fluctuating dynamic models are examined, with fluctuations evaluated using molecular dynamics (MD) simulations. The correlation times for the fundamental bond vector time correlation function and for the second‐order bond orientational TCF are obtained as a function of the residue number for vnd/NK‐2. The stability of the corresponding local dynamics pattern for the fluctuating structure as a function of the length of the MD trajectory is presented. Diffusive dynamics, which is essentially free of model parameters even at first order in the mode‐coupling diffusion approach, confirm that local dynamics of proteins can be described in terms of rotational diffusion of a fluctuating quasi‐rigid structure. The comparison with the nmr data shows that the first‐order mode coupling diffusion approximation accounts for the correct order of magnitude of the results and of important qualitative aspects of the data sensitive to conformational changes. Indications are obtained from this study to efficiently extend the theory to higher order in the mode‐coupling expansion. These results demonstrate the promise of the mode‐coupling approach, where the local dynamics of proteins is described in terms of rotational diffusion of a fluctuating quasi‐rigid structure, to analyze nmr spin–lattice relaxation behavior. © 1999 John Wiley & Sons, Inc. Biopoly 49: 235–254, 1999     

Hydrogen-deuterium exchange mass spectrometry for investigation of backbone dynamics of oxidized and reduced cytochrome P450cam

Backbone dynamics of the camphor monoxygenase cytochrome P450(cam) (CYP101) as a function of oxidation/ligation state of the heme iron were investigated via hydrogen/deuterium exchange (H/D exchange) as monitored by mass spectrometry. Main chain amide NH hydrogens can exchange readily with solvent and the rate of this exchange depends upon, among other things, dynamic fluctuations in local structural elements. A fluxional region of the polypeptide will exchange more quickly with solvent than one that is more constrained. In most regions of the enzyme, exchange rates were similar between oxidized high-spin camphor-bound and reduced camphor- and CO-bound CYP101 (CYP-S and CYP-S-CO, respectively). However, in regions of the protein that have previously been implicated in substrate access by structural and molecular dynamics investigations, the reduced enzyme shows significantly slower exchange rates than the oxidized CYP-S. This observation corresponds to increased flexibility of the oxidized enzyme relative to the reduced form. Structural features previously found to be perturbed in CYP-S-CO upon binding of the biologically relevant effector and reductant putidaredoxin (Pdx) as determined by nuclear magnetic resonance are also more protected from exchange in the reduced state. To our knowledge, this study represents the first experimental investigation of backbone dynamics within the P450 family using this methodology.     

Cytochrome b5 increases the rate of product formation by cytochrome P450 2B4 and competes with cytochrome P450 reductase for a binding site on cytochrome P450 2B4

The kinetics of product formation by cytochrome P450 2B4 were compared in the presence of cytochrome b(5) (cyt b(5)) and NADPH-cyt P450 reductase (CPR) under conditions in which cytochrome P450 (cyt P450) underwent a single catalytic cycle with two substrates, benzphetamine and cyclohexane. At a cyt P450:cyt b(5) molar ratio of 1:1 under single turnover conditions, cyt P450 2B4 catalyzes the oxidation of the substrates, benzphetamine and cyclohexane, with rate constants of 18 +/- 2 and 29 +/- 4.5 s(-1), respectively. Approximately 500 pmol of norbenzphetamine and 58 pmol of cyclohexanol were formed per nmol of cyt P450. In marked contrast, at a cyt P450:CPR molar ratio of 1:1, cyt P450 2B4 catalyzes the oxidation of benzphetamine congruent with100-fold (k = 0.15 +/- 0.05 s(-1)) and cyclohexane congruent with10-fold (k = 2.5 +/- 0.35 s(-1)) more slowly. Four hundred picomoles of norbenzphetamine and 21 pmol of cyclohexanol were formed per nmol of cyt P450. In the presence of equimolar concentrations of cyt P450, cyt b(5), and CPR, product formation is biphasic and occurs with fast and slow rate constants characteristic of catalysis by cyt b(5) and CPR. Increasing the concentration of cyt b(5) enhanced the amount of product formed by cyt b(5) while decreasing the amount of product generated by CPR. Under steady-state conditions at all cyt b(5):cyt P450 molar ratios examined, cyt b(5) inhibits the rate of NADPH consumption. Nevertheless, at low cyt b(5):cyt P450 molar ratios 相似文献   

Coupling and uncoupling mechanisms in the methoxythreonine mutant of cytochrome P450cam: a quantum mechanical/molecular mechanical study

The Thr252 residue plays a vital role in the catalytic cycle of cytochrome P450cam during the formation of the active species (Compound I) from its precursor (Compound 0). We investigate the effect of replacing Thr252 by methoxythreonine (MeO-Thr) on this protonation reaction (coupling) and on the competing formation of the ferric resting state and H₂O₂ (uncoupling) by combined quantum mechanical/molecular mechanical (QM/MM) methods. For each reaction, two possible mechanisms are studied, and for each of these the residues Asp251 and Glu366 are considered as proton sources. The computed QM/MM barriers indicate that uncoupling is unfavorable in the case of the Thr252MeO-Thr mutant, whereas there are two energetically feasible proton transfer pathways for coupling. The corresponding rate-limiting barriers for the formation of Compound I are higher in the mutant than in the wild-type enzyme. These findings are consistent with the experimental observations that the Thr252MeO-Thr mutant forms the alcohol product exclusively (via Compound I), but at lower reaction rates compared with the wild-type enzyme.     

Cloning and characterization of the rat cytochrome P450 4F5 (CYP4F5) gene

The analysis of a non-redundant set of human proteins, for which both the crystallographic structures and the corresponding gene sequences are available, show that bases at third codon position are non-uniformly distributed along the coding sequences. Significant compositional differences are found by comparing the gene regions corresponding to the different secondary structures of the proteins. Inter-and intra-structure differences were most pronounced in the GC-richest genes. These results are not compatible with any proposed hypotheses based on a neutral process of formation/maintenance of the high GC₃ levels of the genes localized in the GC-richest isochores of the human genome.     

Compressibility and uncoupling of cytochrome P450cam: high pressure FTIR and activity studies

The effect of the hydrostatic pressure on the CO ligand stretch vibration in cytochrome P450cam-CO bound with various substrates is studied by FTIR. The vibration frequency is linearily shifted to lower values with increasing pressure. The slope of the shift gives the isothermal compressibility of the heme pocket and is found to be related to the high-spin state content in an opposite direction to that previously observed from the pressure-induced shift of the Soret band. This opposite behaviour is explained by the dual effect of heme pocket water molecules both on the CO ligand and on electrostatic potentials produced by the protein at the distal side. The latter effect disturbs ligand-distal side contacts which are needed for a specific proton transfer in oxygen activation when dioxygen is the ligand. Their loss results in uncoupled H(2)O(2) formation.     

Active-site residues move independently from the rest of the protein in a 200 ns molecular dynamics simulation of cytochrome P450 CYP119

The conformational dynamics of cytochrome P450 enzymes are critical to their catalytic activity. In this study, the correlated motion between residues in a 200 ns molecular dynamics trajectory of the thermophilic CYP119 was analyzed to parse out conformational relationships. Residues that are structurally related, for example residues within a helix, generally have highly correlated motion. In addition, clusters of non-adjacent residues that show correlated motion (“hot spots”) are seen in various regions, including at the base of the F and G helices that make up the most dynamic region of the enzyme. A modified k-means algorithm that clusters residues based on their correlated motion indicates that functionally related residues are in the same cluster (e.g., the catalytic threonines and the heme). Tightly coupled clusters form a solvent-exposed “shell” around the enzyme, whereas less coupling between clusters is seen in regions that are critical to ligand interactions, redox partner interactions, and catalysis. Most notably, we find that residues in the active site move independently from the rest of the enzyme, effectively insulating the catalytic machinery from other regions of the protein.     

Heterotropic and homotropic cooperativity by a drug-metabolising mutant of cytochrome P450 BM3

Recently, we described a triple mutant of the bacterial cytochrome P450 BM3 as the first mutant with affinity for drug-like compounds. In this paper, we show that this mutant, but not wild-type BM3, is able to metabolise testosterone and several drug-like molecules such as amodiaquine, dextromethorphan, acetaminophen, and 3,4-methylenedioxymethylamphetamine that are known substrates of human P450s. Interestingly, the metabolism of 3,4-methylenedioxymethylamphetamine and acetaminophen could be stimulated up to 70-fold by the addition of caffeine, a known activator of rat P450 3A2. With testosterone metabolism, homotropic cooperativity was observed. This shows that heterotropic and homotropic cooperativity, known to occur in the P450 3A family, can also take place in BM3. BM3 therefore can be used as a model system to study atypical kinetics in mammalian P450s. Second, this study shows that BM3 can be engineered to a drug-metabolising enzyme, making it a promising candidate to use as biocatalyst in drug discovery and synthesis.     

Structural and spectroscopic analysis of the F393H mutant of flavocytochrome P450 BM3.

In the preceding paper in this issue [Ost, T. W. B., Miles, C. S., Munro, A. W., Murdoch, J., Reid, G. A., and Chapman, S. K. (2001) Biochemistry 40, 13421-13429], we have established that the primary role of the phylogenetically conserved phenylalanine in flavocytochrome P450 BM3 (F393) is to control the thermodynamic properties of the heme iron, so as to optimize electron-transfer both to the iron (from the flavin redox partner) and onto molecular oxygen. In this paper, we report a detailed study of the F393H mutant enzyme, designed to probe the structural, spectroscopic, and metabolic profile of the enzyme in an attempt to identify the factors responsible for causing the changes. The heme domain structure of the F393H mutant has been solved to 2.0 A resolution and demonstrates that the histidine replaces the phenylalanine in almost exactly the same conformation. A solvent water molecule is hydrogen bonded to the histidine, but there appears to be little other gross alteration in the environment of the heme. The F393H mutant displays an identical ferric EPR spectrum to wild-type, implying that the degree of splitting of the iron d orbitals is unaffected by the substitution, however, the overall energy of the d-orbitals have changed relative to each other. Magnetic CD studies show that the near-IR transition, diagnostic of heme ligation state, is red-shifted by 40 nm in F393H relative to wild-type P450 BM3, probably reflecting alteration in the strength of the iron-cysteinate bond. Studies of the catalytic turnover of fatty acid (myristate) confirms NADPH oxidation is tightly coupled to fatty acid oxidation in F393H, with a product profile very similar to wild-type. The results indicate that gross conformational changes do not account for the perturbations in the electronic features of the P450 BM3 heme system and that the structural environment on the proximal side of the P450 heme must be conformationally conserved in order to optimize catalytic function.     

Molecular dynamics simulation of the substitution of serine for the conserved glycine in the G-loop in the cdc28-srm yeast mutant using the crystal lattice of human CDK2 kinase

Two-nanosecond molecular dynamics simulations of the crystal lattice of an active complex of pT160-CDK2 kinase/cyclin A/ATP-Mg²⁺/substrate were performed. The simulations showed that the structures of the wild-type CDK2 complex and the mutant CDK2 complex involving the substitution G16S-CDK2 corresponding to the yeast substitution G20S-CDC28 differ noticeably and the differences between the structural conformations are most pronounced in the regions that play a key role in the kinase functioning. The results of the computer calculations were used to consider the structural elements that may affect the kinase activity, the regulatory phosphorylation, and the binding of protein kinase with cyclins and substrates.     

On the identity and reactivity patterns of the "second oxidant" of the T252A mutant of cytochrome P450cam in the oxidation of 5-methylenenylcamphor

Density functional calculations show that in the absence of Compound I, the primary oxidant species of P450, the precursor species, Compound 0 (FeOOH), can effect double bond activation of 5-methylenylcamphor (1). The mechanism is initiated by homolytic cleavage of the O–O bond and formation of an OH radical bound to the Compound II species by hydrogen bonding interactions. Subsequently, the so-formed OH radical can either activate the double bond of 1 or attack the meso position of the heme en route to heme degradation. The calculations show that double bond activation is preferred over attack on the heme. Past the double bond activation, the intermediate can either lead to epoxidation or to a glycol formation. The glycol formation is predicted to be preferred, although in the P450_cam pocket the competition may be closer. Therefore, in the absence of Compound I, Compound 0 will be capable of epoxidizing double bonds. Previous studies [E. Derat, D. Kumar, H. Hirao, S. Shaik, J. Am. Chem. Soc. 128 (2006) 473–484] showed that in the case of a substrate that can undergo only C–H activation, the bound OH prefers heme degradation over hydrogen abstraction. Since the epoxidation barrier for Compound I is much smaller than that of Compound 0 (12.8 vs. 18.9 kcal/mol), when Compound I is present in the cycle, Compound 0 will be silent. As such, our mechanism explains lucidly why T252A P450_cam can epoxidize olefins like 5-methylenylcamphor but is ineffective in camphor hydroxylation [S. Jin, T.M. Makris, T. A. Bryson, S.G. Sligar, J.H. Dawson, J. Am. Chem. Soc. 125 (2003) 3406–3407]. Our calculations show that the glycol formation is a marker reaction of Compound 0 with 5-methylenylcamphor. If this product can be found in T252A P450_cam or in similar mutants of other P450 isozymes, this will constitute a more definitive proof for the action of Cpd 0 in P450 enzymes.