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

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

Quantum mechanical/molecular mechanical investigation of the mechanism of C-H hydroxylation of camphor by cytochrome P450cam: theory supports a two-state rebound mechanism

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.     

Is the ruthenium analogue of compound I of cytochrome p450 an efficient oxidant? A theoretical investigation of the methane hydroxylation reaction

High-valent metal-oxo complexes catalyze C-H bond activation by oxygen insertion, with an efficiency that depends on the identity of the transition metal and its oxidation state. Our study uses density functional calculations and theoretical analysis to derive fundamental factors of catalytic activity, by comparison of a ruthenium-oxo catalyst with its iron-oxo analogue toward methane hydroxylation. The study focuses on the ruthenium analogue of the active species of the enzyme cytochrome P450, which is known to be among the most potent catalysts for C-H activation. The computed reaction pathways reveal one high-spin (HS) and two low-spin (LS) mechanisms, all nascent from the low-lying states of the ruthenium-oxo catalyst (Ogliaro, F.; de Visser, S. P.; Groves, J. T.; Shaik, S. Angew. Chem. Int. Ed. 2001, 40, 2874-2878). These mechanisms involve a bond activation phase, in which the transition states (TS's) appear as hydrogen abstraction species, followed by a C-O bond making phase, through a rebound of the methyl radical on the metal-hydroxo complex. However, while the HS mechanism has a significant rebound barrier, and hence a long lifetime of the radical intermediate, by contrast, the LS ones are effectively concerted with small barriers to rebound, if at all. Unlike the iron catalyst, the hydroxylation reaction for the ruthenium analogue is expected to follow largely a single-state reactivity on the LS surface, due to a very large rebound barrier of the HS process and to the more efficient spin crossover expected for ruthenium. As such, ruthenium-oxo catalysts (Groves, J. T.; Shalyaev, K.; Lee, J. In The Porphyrin Handbook; Biochemistry and Binding: Activation of Small Molecules, Vol. 4; Kadish, K. M., Smith, K. M., Guilard, R., Eds.; Academic Press: New York, 2000; pp 17-40) are expected to lead to more stereoselective hydroxylations compared with the corresponding iron-oxo reactions. It is reasoned that the ruthenium-oxo catalyst should have larger turnover numbers compared with the iron-oxo analogue, due to lesser production of suicidal side products that destroy the catalyst (Ortiz de Montellano, P. R.; Beilan, H. S.; Kunze, K. L.; Mico, B. A. J. Biol. Chem. 1981, 256, 4395-4399). The computations reveal also that the ruthenium complex is more electrophilic than its iron analogue, having lower hydrogen abstraction barriers. These reactivity features of the ruthenium-oxo system are analyzed and shown to originate from a key fundamental factor, namely, the strong 4d(Ru)-2p(O,N) overlaps, which produce high-lying pi(Ru-O), sigma(Ru-O), and sigma(Ru-N) orbitals and thereby to lead to a preference of ruthenium for higher-valent oxidation states with higher electrophilicity, for the effectively concerted LS hydroxylation mechanism, and for less suicidal complexes. As such, the ruthenium-oxo species is predicted to be a more robust catalyst than its iron-oxo analogue.     

细胞色素P450催化4-氯-N-环丙基-N-异丙基苯胺Ca-H羟基化反应机理的理论研究

采用密度泛函理论的B3LYP方法,研究了细胞色素P450催化4-氯-N-环丙基-N-异丙基苯胺Ca-H羟基化的反应机理,该反应包含环丙基的羟基化和异丙基的羟基化两个反应途径,且这两个反应路径都是包含氢原子传递的协同过程,二重态的能垒明显低于四重态,反应主要在二重态上进行.通过比较这两个反应路径中Ca-H羟基化反应的活化...     

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.     

Kinetic isotope effect is a sensitive probe of spin state reactivity in C-H hydroxylation of N,N-dimethylaniline by cytochrome P450

This paper presents DFT calculations of C-H hydroxylation of N,N-dimethylaniline by Compound I (Cpd I) of cytochrome P450. The reaction involves two processes nascent from the two spin states of Cpd I, the low-spin (LS) and high-spin (HS) states. The calculations demonstrate that the kinetic isotope effects (KIEs) of the two processes are very different, and only KIELS fits the experimental datum. As such, KIE can be a sensitive probe of spin state reactivity.     

Revisiting the mechanism of P450 enzymes with the radical clocks norcarane and spiro[2,5]octane

Norcarane (1) and spiro[2.5]octane (2) yield different product distributions depending on whether they are oxidized via concerted, radical, or cationic mechanisms. For this reason, these two probes were used to investigate the mechanisms of hydrocarbon hydroxylation by two mammalian and two bacterial cytochrome P450 enzymes. Products indicative of a radical intermediate with a lifetime ranging from 16 to 52 ps were detected during the oxidation of norcarane by P450(cam) (CYP101), P450(BM3) (CYP102), CYP2B1, and CYP2E1. Trace amounts of the cation rearrangement product were observed with norcarane for all but CYP2E1, while no cation or radical rearrangement products were observed for spiro[2.5]octane. The results for the oxidation of norcarane with a radical rearrangement rate of 2 x 10(8) s(-1) are consistent with the involvement of a two-state radical rebound mechanism, while for the slower (5 x 10(7) s(-1)) spiro[2,5]oct-4-yl radical rearrangement products were beyond detection. Taken together with earlier data for the hydroxylation of bicyclo[2.1.0]pentane, which also suggested a 50 ps radical lifetime, these three structurally similar and functionally simple substrates show a consistent pattern of rearrangement that supports a radical rebound mechanism for this set of cytochrome P450 enzymes.     

Theoretical study of the mechanism of valence tautomerism in cobalt complexes

Valence tautomerism is studied in the [Co(II-HS)(sq)(2)(bpy)]/[Co(III-LS)(sq)(cat)(bpy)] mononuclear cobalt complex by using DFT methods (HS, high spin; LS, low spin; cat, catecholate; sq, semiquinone; bpy, 2,2'-bipyridine). Calculations at the B3LYP* level of theory reproduce well the energy gap between the Co(II-HS) and Co(III-LS) forms giving an energy gap of 4.4 kcal/mol, which is comparable to the experimental value of 8.9 kcal/mol. Potential energy surfaces and crossing seams of the electronic states of the doublet, quartet, and sextet spin states are calculated along minimum energy paths connecting the energy minima corresponding to the different spin states. The calculated minimum energy crossing points (MECPs) are located at 8.8 kcal/mol in the doublet/sextet surfaces, at 10.2 kcal/mol in the doublet/quartet surfaces, and at 8.4 kcal/mol in the quartet/sextet surfaces relative to the doublet ground state. Considering the energy of the three spin states and the crossing points, the one-step relaxation mechanism between the Co(II-HS) and Co(III-LS) forms is the most probable. This research shows that mapping MECPs can be a useful strategy to analyze the potential energy surfaces of systems with complex deformation modes.     

Interplay between two spin states determines the hydroxylation catalyzed by P450 monooxygenase from Trichoderma brevicompactum

Tri11 (now renamed as tri22) encoded cytochrome P₄₅₀ monooxygenase in Trichoderma brevicompactum catalyzes the C-4 C-H hydroxylation of 12, 13-epoxytrichothec-9-ene (EPT) to produce trichodermol in the biosynthetic pathway of trichodermin/harzianum A. The density functional theory (DFT)-quantum mechanics (QM) approach is applied to elucidate the hydroxylation of EPT by using a model active species of P₄₅₀ (Cpd I). The QM calculations were performed on the active site complex, to find out transition-state structure, intermediate, and product complexes for the two spin states at different potential energy surfaces. The two state reactivity rebound-free product formation resulted from the interplay of two spin states (doublet and quartet).     

Kinetic isotope effects implicate two electrophilic oxidants in cytochrome p450-catalyzed hydroxylations

Intramolecular kinetic isotope effects (KIEs) were determined for cytochrome P450-catalyzed hydroxylation reactions of methyl-dideuterated trans-2-phenylcyclopropylmethane-d2 (1-d2), which gives two products from oxidation of the methyl group, trans-2-phenylcyclopropylmethanol (2) and 1-phenyl-3-buten-1ol (3). In oxidations of each enantiomer of 1-d2 with three P450 enzymes (CYP2B1, CYPDelta2E1, and CYPDelta2E1 T303A), the apparent intramolecular KIEs were different for products 2 and 3 in all cases and different for each enzyme-substrate combination. In oxidations of each enantiomer of undeuterated 1-d0 and trideuteriomethyl 1-d3 by CYP2B1 and CYPDelta2E1, the ratio of products 2/3 decreased for 1-d3 in comparison to 1-d0 in all cases. The results require multiple pathways for P450-catalyzed hydroxylation and are consistent with the "two-oxidants" model, where hydroxylation is effected by both the hydroperoxy-iron species and the iron-oxo species. The results are not consistent with predictions of the "two-states" model for P450-catalyzed hydroxylations, where oxidations occur from a low-spin state and a high-spin state of iron-oxo.     

Multistate reactivity in styrene epoxidation by compound I of cytochrome p450: mechanisms of products and side products formation

Density functional theoretical calculations are used to elucidate the epoxidation mechanism of styrene with a cytochrome P450 model Compound I, and the formation of side products. The reaction features multistate reactivity (MSR) with different spin states (doublet and quartet) and different electromeric situations having carbon radicals and cations, as well as iron(III) and iron(IV) oxidation states. The mechanisms involve state-specific product formation, as follows: a) The low-spin pathways lead to epoxide formation in effectively concerted mechanisms. b) The high-spin pathways have finite barriers for ring-closure and may have a sufficiently long lifetime to undergo rearrangement and lead to side products. c) The high-spin radical intermediate, (4)2(rad)-IV, has a ring closure barrier as small as the C--C rotation barrier. This intermediate will therefore lose stereochemistry and lead to a mixture of cis and trans epoxides. The barriers for the production of aldehyde and suicidal complexes are too high for this intermediate. d) The high-spin radical intermediate, (4)2(rad)-III, has a substantial ring closure barrier and may survive long enough time to lead to suicidal, phenacetaldehyde and 2-hydroxostyrene side products. e) The phenacetaldehyde and 2-hydroxostyrene products both originate from crossover from the (4)2(rad)-III radical intermediate to the cationic state, (4)2(cat,z(2) ). The process involves an N-protonated porphyrin intermediate that re-shuttles the proton back to the substrate to form either phenacetaldehyde or 2-hydroxostyrene products. This resembles the internally mediated NIH-shift observed during benzene hydroxylation.     

Multireference ab initio quantum mechanics/molecular mechanics study on intermediates in the catalytic cycle of cytochrome P450(cam)

We have investigated the elusive reactive species of cytochrome P450(cam) (Compound I), the hydroxo complex formed during camphor hydroxylation, and the ferric hydroperoxo complex (Compound 0) by combined quantum mechanical/molecular mechanical (QM/MM) calculations, employing both density functional theory (DFT) and correlated ab initio methods. The first two intermediates appear multiconfigurational in character, especially in the doublet state and less so in the quartet state. DFT(B3LYP)/MM calculations reproduce the relative energies from correlated ab initio QM/MM treatments quite well, except for the splitting of the lowest A(1u)-A(2u) radical states. The inclusion of dynamic correlation is crucial for the proper ab initio treatment of these intermediates.     

Hydroxylation by the hydroperoxy-iron species in cytochrome P450 enzymes

Intramolecular and intermolecular kinetic isotope effects (KIEs) were determined for hydroxylation of the enantiomers of trans-2-(p-trifluoromethylphenyl)cyclopropylmethane (1) by hepatic cytochrome P450 enzymes, P450s 2B1, Delta2B4, Delta2B4 T302A, Delta2E1, and Delta2E1 T303A. Two products from oxidation of the methyl group were obtained, unrearranged trans-2-(p-trifluoromethylphenyl)cyclopropylmethanol (2) and rearranged 1-(p-trifluoromethylphenyl)but-3-en-1-ol (3). In intramolecular KIE studies with dideuteriomethyl substrates (1-d(2)) and in intermolecular KIE studies with mixtures of undeuterated (1-d(0)) and trideuteriomethyl (1-d(3)) substrates, the apparent KIE for product 2 was consistently larger than the apparent KIE for product 3 by a factor of ca. 1.2. Large intramolecular KIEs found with 1-d(2) (k(H)/k(D) = 9-11 at 10 degrees C) were shown not to be complicated by tunneling effects by variable temperature studies with two P450 enzymes. The results require two independent isotope-sensitive processes in the overall hydroxylation reactions that are either competitive or sequential. Intermolecular KIEs were partially masked in all cases and largely masked for some P450s. The intra- and intermolecular KIE results were combined to determine the relative rate constants for the unmasking and hydroxylation reactions, and a qualitative correlation was found for the unmasking reaction and release of hydrogen peroxide from four of the P450 enzymes in the absence of substrate. The results are consistent with the two-oxidants model for P450 (Vaz, A. D. N.; McGinnity, D. F.; Coon, M. J. Proc. Natl. Acad. Sci. U.S.A. 1998, 95, 3555), which postulates that a hydroperoxy-iron species (or a protonated analogue of this species) is a viable electrophilic oxidant in addition to the consensus oxidant, iron-oxo.     

Mechanistic Insights into the Regio- and Stereoselectivities of Testosterone and Dihydrotestosterone Hydroxylation Catalyzed by CYP3A4 and CYP19A1

The hydroxylation of nonreactive C−H bonds can be easily catalyzed by a variety of metalloenzymes, especially cytochrome P450s (P450s). The mechanism of P450 mediated hydroxylation has been intensively studied, both experimentally and theoretically. However, understanding the regio- and stereoselectivities of substrates hydroxylated by P450s remains a great challenge. Herein, we use a multi-scale modeling approach to investigate the selectivity of testosterone (TES) and dihydrotestosterone (DHT) hydroxylation catalyzed by two important P450s, CYP3A4 and CYP19A1. For CYP3A4, two distinct binding modes for TES/DHT were predicted by dockings and molecular dynamics simulations, in which the experimentally identified sites of metabolism of TES/DHT can access to the catalytic center. The regio- and stereoselectivities of TES/DHT hydroxylation were further evaluated by quantum mechanical and ONIOM calculations. For CYP19A1, we found that sites 1β, 2β and 19 can access the catalytic center, with the intrinsic reactivity 2β>1β>19. However, our ONIOM calculations indicate that the hydroxylation is favored at site 19 for both TES and DHT, which is consistent with the experiments and reflects the importance of the catalytic environment in determining the selectivity. Our study unravels the mechanism underlying the selectivity of TES/DHT hydroxylation mediated by CYP3A4 and CYP19A1 and is helpful for understanding the selectivity of other substrates that are hydroxylated by P450s.     

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.     

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

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

Enhanced electron transfer and lauric acid hydroxylation by site-directed mutagenesis of CYP119

CYP119, a cytochrome P450 from a thermophilic organism for which a crystal structure is available, is shown here to hydroxylate lauric acid in a reaction supported by putidaredoxin and putidaredoxin reductase. This fatty acid hydroxylation activity is increased 15-fold by T214V and D77R mutations. The T214V mutation increases the rate by facilitating substrate binding and enhancing the associated spin state change, whereas the D77R mutation improves binding of the heterologous redox partner putidaredoxin to CYP119 and the rate of electron transfer from it to the heme group. A sequence alignment with P450(cam) can, therefore, be used to identify a part of the binding site for putidaredoxin on an unrelated P450 enzyme. This information can be used to engineer by mutagenesis an improved complementarity of the protein-protein interface that results in improved electron transfer from putidaredoxin to the P450 enzyme. As a result, the catalytic activity of the thermo- and barostable CYP119 has been incorporated into a catalytic system that hydroxylates fatty acids.     

The CYPome of Sorangium cellulosum So ce56 and identification of CYP109D1 as a new fatty acid hydroxylase

The first systematic study of the complete cytochrome P450 complement (CYPome) of Sorangium cellulosum So ce56, which is a producer of important secondary metabolites and has the largest bacterial?genome sequenced to date, is presented. We describe the bioinformatic analysis of the So ce56 cytochrome P450 complement consisting of 21 putative P450 genes. Because fatty acids play a pivotal role during the complex life cycle of myxobacteria, we focused our studies on the characterization of fatty acid hydroxylases. Three novel potential fatty acid hydroxylases (CYP109D1, CYP264A1, and CYP266A1) were used for detailed characterization. One of them, CYP109D1 was able to perform subterminal hydroxylation of saturated fatty acids with the support of two autologous and one heterologous electron transfer system(s). The kinetic parameters for the product hydroxylation were derived.     

CYP505E3: A Novel Self-Sufficient ω-7 In-Chain Hydroxylase

The self-sufficient cytochrome P450 monooxygenase CYP505E3 from Aspergillus terreus catalyzes the regioselective in-chain hydroxylation of alkanes, fatty alcohols, and fatty acids at the ω-7 position. It is the first reported P450 to give regioselective in-chain ω-7 hydroxylation of C10–C16 n-alkanes, thereby enabling the one step biocatalytic synthesis of rare alcohols such as 5-dodecanol and 7-tetradecanol. It shows more than 70 % regioselectivity for the eighth carbon from one methyl terminus, and displays remarkably high activity towards decane (TTN≈8000) and dodecane (TTN≈2000). CYP505E3 can be used to synthesize the high-value flavour compound δ-dodecalactone via two routes: 1) conversion of dodecanoic acid into 5-hydroxydodecanoic acid (24 % regioselectivity), which at low pH lactonises to δ-dodecalactone, and 2) conversion of 1-dodecanol into 1,5-dodecanediol (55 % regioselectivity), which can be converted into δ-dodecalactone by horse liver alcohol dehydrogenase.