No abstract

A two-state rebound mechanism of alkane hydroxylation by a model active species of the enzyme cytochrome P450 is studied using density functional theoretic calculations. Theory corroborates Groves's rebound mechanism (Groves, J. T. J. Chem. Educ. 1985, 62, 928), with a key difference, namely that in the two-state rebound the reactivity and product distribution result from the interplay of two reactive states of the active ferryl-oxene (Por+•FeO) species of the enzyme:  one state is low-spin (doublet) and the other high-spin (quartet). Transition-state structures, intermediates, and product complexes are identified for the two states. The bond activation in either one of the two states involves a hydrogen abstraction-like transition structure. However, while in the high-spin state there forms a radical that has a significant barrier for rebound, in the low-spin state the rebound is virtually barrierless. Even though one cannot ignore incursion of a small amount of radicals in the low-spin state, it is clear that the radical has a significant lifetime mainly on the high-spin surface. The results are used to gain insight into puzzling experimental data which emerge from studies of ultrafast radical clocks (e.g., Toy, P. H.; Newcomb, M.; Hollenberg, P. F., J. Am. Chem. Soc. 1998, 120, 7719), vis à vis the nature the transition state, deduced from kinetic isotope effect measurements (Manchester, J. I.; Dinnocenzo, J. P.; Higgins, L. A.; Jones, J. P. J. Am. Chem. Soc. 1997, 119, 5069) and stereochemical scrambling patterns (Groves, J. T.; McClusky, G. A.; White, R. E.; Coon, M. J. Biochem. Biophys. Res. Commun. 1978, 81, 154). Understanding the electronic structure of the various species leads to a predictive structure−reactivity picture, based on the two-state reactivity scenario (Shaik, S.; Filatov, M.; Schröder, D.; Schwarz, H. Chem. Eur. J. 1998, 4, 193). The model makes it possible to predict the dependence of the relative rates of the two states, and of the corresponding steps as a function of the nature of the alkane, the resulting alkyl radical, and the binding capability of the thiolate proximal ligand of the active species.

The primary oxidant of cytochrome P450 enzymes, Compound I, is hard to detect experimentally; in the case of cytochrome P450(cam), this intermediate does not accumulate in solution during the catalytic cycle even at temperatures as low as 200 K (ref 4). Theory can play an important role in characterizing such elusive species. We present here combined quantum mechanical/molecular mechanical (QM/MM) calculations of Compound I of cytochrome P450(cam) in the full enzyme environment as well as density functional studies of the isolated QM region. The calculations assign the ground state of the species, quantify the effect of polarization and hydrogen bonding on its properties, and show that the protein environment and its specific hydrogen bonding to the cysteinate ligand are crucial for sustaining the Fe-S bond and for preventing the full oxidation of the sulfur.

The stereospecific cytochrome P450-catalyzed hydroxylation of the C(5)-H((5-exo)) bond in camphor has been studied theoretically by a combined quantum mechanical/molecular mechanical (QM/MM) approach. Density functional theory is employed to treat the electronic structure of the active site (40-100 atoms), while the protein and solvent environment (ca. 24,000 atoms) is described by the CHARMM force field. The calculated energy profile of the hydrogen-abstraction oxygen-rebound mechanism indicates that the reaction takes place in two spin states (doublet and quartet), as has been suggested earlier on the basis of calculations on simpler models ("two-state reactivity"). While the reaction on the doublet potential energy surface is nonsynchronous, yet effectively concerted, the quartet pathway is truly stepwise, including formation of a distinct intermediate substrate radical and a hydroxo-iron complex. Comparative calculations in the gas phase demonstrate the effect of the protein environment on the geometry and relative stability of intermediates (in terms of spin states and redox electromers) through steric constraints and electronic polarization.

Iron(III)-hydroperoxo, [Por(CysS)Fe(III)-OOH](-), a key species in the catalytic cycle of cytochrome P450, was recently identified by EPR/ENDOR spectroscopies (Davydov, R.; Makris, T. M.; Kofman, V.; Werst, D. E.; Sligar, S. G.; Hoffman, B. M. J. Am. Chem. Soc. 2001, 123, 1403-1415). It constitutes the last station of the preparative steps of the enzyme before oxidation of an organic compound and is implicated as the second oxidant capable of olefin epoxidation (Vaz, A. D. N.; McGinnity, D. F.; Coon, M. J. Proc. Natl. Acad. Sci. U.S.A. 1998, 95, 3555-3560), in addition to the penultimate active species, Compound I (Groves, J. T.; Han, Y.-Z. In Cytochrome P450: Structure, Mechanism and Biochemistry, 2nd ed.; Ortiz de Montellano, P. R., Ed.; Plenum Press: New York, 1995; pp 3-48). In response, we present a density functional study of a model species and its ethylene epoxidation pathways. The study characterizes a variety of properties of iron(III)-hydroperoxo, such as the O-O bonding, the Fe-S bonding, Fe-O and Fe-S stretching frequencies, its electron attachment, and ionization energies. Wherever possible these properties are compared with those of Compound I. The proton affinities for protonation on the proximal and distal oxygen atoms of iron(III)-hydroperoxo, and the effect of the thiolate ligand thereof, are determined. In accordance with previous results (Harris, D. L.; Loew, G. H. J. Am. Chem. Soc. 1998, 120, 8941-8948), iron(III)-hydroperoxo is a strong base (as compared with water), and its distal protonation leads to a barrier-free formation of Compound I. The origins of this barrier-free process are discussed using a valence bond approach. It is shown that the presence of the thiolate is essential for this process, in line with the "push effect" deduced by experimentalists (Sono, M.; Roach, M. P.; Coulter, E. D.; Dawson, J. H. Chem. Rev. 1996, 96, 2841-2887). Finally, four epoxidation pathways of iron(III)-hydroxperoxo are located, in which the species transfers oxygen to ethylene either from the proximal or from the distal sites, in both concerted and stepwise manners. The barriers for the four mechanisms are 37-53 kcal/mol, in comparison with 14 kcal/mol for epoxidation by Compound I. It is therefore concluded that iron(III)-hydroperoxo, as such, cannot be a second oxidant, in line with its significant basicity and poor electron-accepting capability. Possible versions of a second oxidant are discussed.

No abstract

The communication presents DFT calculations of 10 different C-H hydroxylation barriers by the active species of the enzyme cytochrome P450. The work demonstrates the existence of an excellent barrier-bond energy correlation. The so-obtained equation of the straight line is demonstrated to be useful for predicting barriers of related C-H activation processes, as well as for assessing barrier heights within the protein environment. This facility is demonstrated be estimating the barrier of camphor hydroxylation by P450cam.

C−H hydroxylation by the enzyme cytochrome P450 is one of Nature’s important and most ubiquitous processes. There is strong evidence that the mechanism proceeds by initial hydrogen abstraction, from the alkane, by the high‐valent iron−oxo species of the enzyme, followed by a rebound of the alkyl radical to form the ferric−alcohol product complex (the Groves “rebound” mechanism). Nevertheless, the “rebound” mechanism is still controversial due to the ultrashort radical lifetimes deduced from radical‐clock experiments. This review describes the main elements of the controversy and its updated resolution by theory, with an emphasis on the controversial rebound step. The theoretically derived model for alkane hydroxylation is found to involve two‐state reactivity (TSR). In TSR, radicals are produced on two different spin‐state surfaces, and thereafter they react differently; on the low‐spin surface, the rebound proceeds with no product rearrangement and the lifetime of the radicals is ultrashort (or zero), while on the high‐spin surface the barrier for rebound is substantial and the lifetime of the radical is sufficiently long that rearrangement may compete with product formation by rebound. A new valence‐bond model is developed to model the rebound barrier on the high‐spin surface and conceptualize its dependence on the nature of the alkyl radical. The possible intermediacy of carbocationic intermediates alongside radicals is discussed. It is shown that the TSR scenario provides a satisfactory rationale for the controversial findings in the field, and makes verifiable predictions. One of the predictions of TSR is that the ratio of unrearranged to rearranged products, [U/R], will be subject to an intrinsic isotope effect that is substrate‐dependent. This prediction and its possible verification by experiment are discussed. (© Wiley‐VCH Verlag GmbH & Co. KGaA, 69451 Weinheim, Germany, 2004)