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.     

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.     

Anilinic N‐Oxides Support Cytochrome P450‐Mediated N‐Dealkylation through Hydrogen‐Atom Transfer

The mechanism of N‐dealkylation mediated by cytochrome P450 (P450) has long been studied and argued as either a single electron transfer (SET) or a hydrogen atom transfer (HAT) from the amine to the oxidant of the P450, the reputed iron–oxene. In our study, tertiary anilinic N‐oxides were used as oxygen surrogates to directly generate a P450‐mediated oxidant that is capable of N‐dealkylating the dimethylaniline derived from oxygen donation. These surrogates were employed to probe the generated reactive oxygen species and the subsequent mechanism of N‐dealkylation to distinguish between the HAT and SET mechanisms. In addition to the expected N‐demethylation of the product aniline, 2,3,4,5,6‐pentafluoro‐N,N‐dimethylaniline N‐oxide (PFDMAO) was found to be capable of N‐dealkylating both N,N‐dimethylaniline (DMA) and N‐cyclopropyl‐N‐methylaniline (CPMA). Rate comparisons of the N‐demethylation of DMA supported by PFDMAO show a 27‐fold faster rate than when supported by N,N‐dimethylaniline N‐oxide (DMAO). Whereas intermolecular kinetic isotope effects were masked, intramolecular measurements showed values reflective of those seen previously in DMAO‐ and the native NADPH/O₂‐supported systems (2.33 and 2.8 for the N‐demethylation of PFDMA and DMA from the PFDMAO system, respectively). PFDMAO‐supported N‐dealkylation of CPMA led to the ring‐intact product N‐cyclopropylaniline (CPA), similar to that seen with the native system. The formation of CPA argues against a SET mechanism in favor of a P450‐like HAT mechanism. We suggest that the similarity of KIEs, in addition to the formation of the ring‐intact CPA, argues for a similar mechanism of Compound I (Cpd I) formation followed by HAT for N‐dealkylation by the native and N‐oxide‐supported systems and demonstrate the ability of the N‐oxide‐generated oxidant to act as an accurate mimic of the native P450 oxidant.     

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.     

Theoretical study of N-demethylation of substituted N,N-dimethylanilines by cytochrome P450: the mechanistic significance of kinetic isotope effect profiles

The mechanism of N-demethylation of N,N-dimethylanilines (DMAs) by cytochrome P450, a highly debated topic in mechanistic bioinorganic chemistry (Karki, S. B.; Dinnocenczo, J. P.; Jones, J. P.; Korzekwa, K. R. J. Am. Chem. Soc. 1995, 117, 3657), is studied here using DFT calculations of the reactions of the active species of the enzyme, Compound I (Cpd I), with four para-(H, Cl, CN, NO2) substituted DMAs. The calculations resolve mechanistic controversies, offer a consistent mechanistic view, and reveal the following features: (a) the reaction pathways involve C-H hydroxylation by Cpd I followed by a nonenzymatic carbinolamine decomposition. (b) C-H hydroxylation is initiated by a hydrogen atom transfer (HAT) step that possesses a "polar" character. As such, the HAT energy barriers correlate with the energy level of the HOMO of the DMAs. (c) The series exhibits a switch from spin-selective reactivity for DMA and p-Cl-DMA to two-state reactivity, with low- and high-spin states, for p-CN-DMA and p-NO2-DMA. (d) The computed kinetic isotope effect profiles (KIEPs) for these scenarios match the experimentally determined KIEPs. Theory further shows that the KIEs and TS structures vary in a manner predicted by the Melander-Westheimer postulate: as the substituent becomes more electron withdrawing, the TS is shifted to a later position along the H-transfer coordinate and the corresponding KIEs increases. (e) The generated carbinolaniline can readily dissociate from the heme and decomposes in a nonenzymatic environment, which involves water assisted proton shift.     

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.     

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

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

Characterizing the Intermediates Compound I and II in the Cytochrome P450 Catalytic Cycle with Nonlinear X‐ray Spectroscopy: A Simulation Study

Cytochrome P450 enzymes are an important family of biocatalysts that oxidize chemically inert C?H bonds. There are many unresolved questions regarding the catalytic reaction intermediates, in particular P450 Compound I (Cpd‐I) and II (Cpd‐II). By using simple molecular models, we simulate various X‐ray spectroscopy signals, including X‐ray absorption near‐edge structure (XANES), resonant inelastic X‐ray scattering (RIXS), and stimulated X‐ray Raman spectroscopy (SXRS) of the low‐ and high‐spin states of Cpd‐I and II. Characteristic peak patterns are presented and connected to the corresponding electronic structures. These X‐ray spectroscopy techniques are complementary to more conventional infrared and optical spectroscopy and they help to elucidate the evolving electronic structures of transient species along the reaction path.     

New features in the catalytic cycle of cytochrome P450 during the formation of compound I from compound 0

Density functional theory (DFT) is applied to the dark section of the catalytic cycle of the enzyme cytochrome P450, namely, the formation of the active species, Compound I (Cpd I), from the ferric-hydroperoxide species (Cpd 0) by a protonation-assisted mechanism. The chosen 96-atom model includes the key functionalities deduced from experiment: Asp(251), Thr(252), Glu(366), and the water channels that relay the protons. The DFT model calculations show that (a) Cpd I is not formed spontaneously from Cpd 0 by direct protonation, nor is the process very exothermic. The process is virtually thermoneutral and involves a significant barrier such that formation of Cpd I is not facile on this route. (b) Along the protonation pathway, there exists an intermediate, a protonated Cpd 0, which is a potent oxidant since it is a ferric complex of water oxide. Preliminary quantum mechanical/molecular mechanical calculations confirm that Cpd 0 and Cpd I are of similar energy for the chosen model and that protonated Cpd 0 may exist as an unstable intermediate. The paper also addresses the essential role of Thr(252) as a hydrogen-bond acceptor (in accord with mutation studies of the OH group to OMe).     

Direct Comparison of the Reactivity of Model Complexes for Compounds 0, I,and II in Oxygenation,Hydrogen‐Abstraction,and Hydride‐Transfer Processes

The iron(III) meso‐tetramesitylporphyrin complex is a good biomimetic to study the catalytic reactions of cytochrome P450. All of the three most discussed reactive intermediates concerning P450 catalysis (namely, Cpd 0, Cpd I, and Cpd II) can be selectively produced, identified, and stabilized for many minutes in solution at low temperature by choosing appropriate reaction conditions. In this way, their reactivity against various substrates was determined by utilizing low‐temperature rapid‐scan UV/Vis spectroscopy. Since all reactive intermediates are derived from a single model complex, the results of these kinetic measurements provide for the first time a full comparability of the determined rate constants for the three intermediates. The rate constants reveal a significant dependence of the reactivity on the type of reaction (e.g., oxygenation, hydrogen abstraction, or hydride transfer), which closely correlates with the chemical nature of Cpds 0, I, and II. The detailed knowledge of the reactivity of these intermediates provides a valuable tool to evaluate their particular role in biological systems.     

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.     

A Compound I Mimic Reveals the Transient Active Species of a Cytochrome P450 Enzyme: Insight into the Stereoselectivity of P450-Catalysed Oxidations

Catching the structure of cytochrome P450 enzymes in flagrante is crucial for the development of P450 biocatalysts, as most structures collected are found trapped in a precatalytic conformation. At the heart of P450 catalysis lies Cpd I, a short-lived, highly reactive intermediate, whose recalcitrant nature has thwarted most attempts at capturing catalytically relevant poses of P450s. We report the crystal structure of P450BM3 mimicking the state in the precise moment preceding epoxidation, which is in perfect agreement with the experimentally observed stereoselectivity. This structure was attained by incorporation of the stable Cpd I mimic oxomolybdenum mesoporphyrin IX into P450BM3 in the presence of styrene. The orientation of styrene to the Mo-oxo species in the crystal structures sheds light onto the dynamics involved in the rotation of styrene to present its vinyl group to Cpd I. This method serves as a powerful tool for predicting and modelling the stereoselectivity of P450 reactions.     

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.     

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.     

Metabolic‐intermediate complex formation with cytochrome P450: Theoretical studies in elucidating the reaction pathway for the generation of reactive nitroso intermediate

Mechanism‐based inhibition (MBI) of cytochrome P450 (CYP) can lead to drug–drug interactions and often to toxicity. Some aliphatic and aromatic amines can undergo biotransformation reactions to form reactive metabolites such as nitrosoalkanes, leading to MBI of CYPs. It has been proposed that the nitrosoalkanes coordinate with the heme iron, forming metabolic‐intermediate complex (MIC), resulting in the quasi‐irreversible inhibition of CYPs. Limited mechanistic details regarding the formation of reactive nitroso intermediate and its coordination with heme‐iron have been reported. A quantum chemical analysis was performed to elucidate potential reaction pathways for the generation of nitroso intermediate and the formation of MIC. Elucidation of the energy profile along the reaction path, identification of three‐dimensional structures of reactive intermediates and transition states, as well as charge and spin density analyses, were performed using the density functional B3LYP method. The study was performed using Cpd I [iron (IV‐oxo] heme porphine with SH^? as the axial ligand) to represent the catalytic domain of CYP, simulating the biotransformation process. Three pathways: (i) N‐oxidation followed by proton shuttle, (ii) N‐oxidation followed by 1,2‐H shift, and (iii) H‐abstraction followed by rebound mechanism, were studied. It was observed that the proton shuttle pathway was more favorable over the whole reaction leading to reactive nitroso intermediate. This study revealed that the MIC formation from a primary amine is a favorable exothermic process, involving eight different steps and preferably takes place on the doublet spin surface of Cpd I . The rate‐determining step was identified to be the first N‐oxidation of primary amine. © 2012 Wiley Periodicals, Inc.     

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

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

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

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

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.     

QM/MM study of the active species of the human cytochrome P450 3A4, and the influence thereof of the multiple substrate binding

Cytochrome P450 3A4 is involved in the metabolism of 50% of all swallowed drugs. The enzyme functions by means of a high-valent iron-oxo species, called compound I (Cpd I), which is formed after entrance of the substrate to the active site. We explored the features of Cpd I using hybrid quantum mechanical/molecular mechanical calculations on various models that are either substrate-free or containing one and two molecules of diazepam as a substrate. M?ssbauer parameters of Cpd I were computed. Our major finding shows that without the substrate, Cpd I tends to elongate its Fe-S bond, localize the radical on the sulfur, and form hydrogen bonds with A305 and T309, which may hypothetically lead to Cpd I consumption by H-abstraction. However, the positioning of diazepam close to Cpd I, as enforced by the effector molecule, was found to strengthen the NH...S interactions of the conserved I443 and G444 residues with the proximal cysteinate ligand. These interactions are known to stabilize the Fe-S bond, and as such, the presence of the substrate leads to a shorter Fe-S bond and it prevents the localization of the radical on the sulfur. This diazepam-Cpd I stabilization was manifested in the 1W0E conformer. The effector substrate did not influence Cpd I directly but rather by positioning the active substrate close to Cpd I, thus displacing the hydrogen bonds with A305 and T309, and thereby giving preference to substrate oxidation. It is hypothesized that these effects on Cpd I, promoted by the restrained substrate, may be behind the special metabolic behavior observed in cases of multiple substrate binding (also called cooperative binding). This restraint constitutes a mechanism whereby substrates stabilize Cpd I sufficiently long to affect monooxygenation by P450s at the expense of Cpd I destruction by the protein residues.