Substrate selectivity of human cytochrome P450 2C9: importance of residues 476, 365, and 114 in recognition of diclofenac and sulfaphenazole and in mechanism-based inactivation by tienilic acid

Formation of potentially toxic metabolites of drugs in reactions catalyzed by human drug-metabolizing enzymes

Data are presented on the formation of potentially toxic metabolites of drugs that are substrates of human drug metabolizing enzymes. The tabular data lists the formation of potentially toxic/reactive products. The data were obtained from in vitro experiments and showed that the oxidative reactions predominate (with 96% of the total potential toxication reactions). Reductive reactions (e.g., reduction of nitro to amino group and reductive dehalogenation) participate to the extent of 4%. Of the enzymes, cytochrome P450 (P450, CYP) enzymes catalyzed 72% of the reactions, myeloperoxidase (MPO) 7%, flavin-containing monooxygenase (FMO) 3%, aldehyde oxidase (AOX) 4%, sulfotransferase (SULT) 5%, and a group of minor participating enzymes to the extent of 9%. Within the P450 Superfamily, P450 Subfamily 3A (P450 3A4 and 3A5) participates to the extent of 27% and the Subfamily 2C (P450 2C9 and P450 2C19) to the extent of 16%, together catalyzing 43% of the reactions, followed by P450 Subfamily 1A (P450 1A1 and P450 1A2) with 15%. The P450 2D6 enzyme participated in an extent of 8%, P450 2E1 in 10%, and P450 2B6 in 6% of the reactions. All other enzymes participate to the extent of 14%. The data show that, of the human enzymes analyzed, P450 enzymes were dominant in catalyzing potential toxication reactions of drugs and their metabolites, with the major role assigned to the P450 Subfamily 3A and significant participation of the P450 Subfamilies 2C and 1A, plus the 2D6, 2E1 and 2B6 enzymes contributing. Selected examples of drugs that are activated or proposed to form toxic species are discussed.

The molecular basis of dapsone activation of CYP2C9-catalyzed nonsteroidal anti-inflammatory drug oxidation

Positive heterotropic cooperativity, or "activation," results in an instantaneous increase in enzyme activity in the absence of an increase in protein expression. Thus, cytochrome P450 (CYP) enzyme activation presents as a potential drug-drug interaction mechanism. It has been demonstrated previously that dapsone activates the CYP2C9-catalyzed oxidation of a number of nonsteroidal anti-inflammatory drugs in vitro. Here, we conducted molecular dynamics simulations (MDS) together with enzyme kinetic investigations and site-directed mutagenesis to elucidate the molecular basis of the activation of CYP2C9-catalyzed S-flurbiprofen 4'-hydroxylation and S-naproxen O-demethylation by dapsone. Supplementation of incubations of recombinant CYP2C9 with dapsone increased the catalytic efficiency of flurbiprofen and naproxen oxidation by 2.3- and 16.5-fold, respectively. MDS demonstrated that activation arises predominantly from aromatic interactions between the substrate, dapsone, and the phenyl rings of Phe114 and Phe476 within a common binding domain of the CYP2C9 active site, rather than involvement of a distinct effector site. Mutagenesis of Phe114 and Phe476 abrogated flurbiprofen and naproxen oxidation, and MDS and kinetic studies with the CYP2C9 mutants further identified a pivotal role of Phe476 in dapsone activation. MDS additionally showed that aromatic stacking interactions between two molecules of naproxen are necessary for binding in a catalytically favorable orientation. In contrast to flurbiprofen and naproxen, dapsone did not activate the 4'-hydroxylation of diclofenac, suggesting that the CYP2C9 active site favors cooperative binding of nonsteroidal anti-inflammatory drugs with a planar or near-planar geometry. More generally, the work confirms the utility of MDS for investigating ligand binding in CYP enzymes.

Human Cytochrome P450 2C9 and Its Polymorphic Modifications: Electroanalysis, Catalytic Properties, and Approaches to the Regulation of Enzymatic Activity

The electrochemical properties of cytochrome P450 2C9 (CYP2C9) and polymorphic modifications P450 2C9*2 (CYP2C9*2) and P450 2C9*3 (CYP2C9*3) were studied. To analyze the comparative electrochemical and electrocatalytic activity, the enzymes were immobilized on electrodes modified with a membrane-like synthetic surfactant (didodecyldimethylammonium bromide (DDAB)). An adequate choice of the type of modified electrode was confirmed by cyclic voltammetry of cytochromes P450 under anaerobic conditions, demonstrating well-defined peaks of reduction and oxidation of the heme iron. The midpoint potential, Emid, of cytochrome P450 2C9 is −0.318 ± 0.01 V, and Emid = −0.324 ± 0.01 V, and Emid = −0.318 ± 0.03 V for allelic variant 2C9*2 and allelic variant 2C9*3, respectively. In the presence of substrate diclofenac under aerobic conditions, cytochrome P450 2C9 and its polymorphic modifications P450 2C9*2 and P450 2C9*3 exhibit catalytic properties. Stimulation of the metabolism of diclofenac by cytochrome P450 2C9 in the presence of antioxidant medications mexidol and taurine was shown.

Differential inhibition of human CYP2C8 and molecular docking interactions elicited by sorafenib and its major N-oxide metabolite

The tyrosine kinase inhibitor sorafenib (SOR) is being used increasingly in combination with other anticancer agents like paclitaxel, but this increases the potential for drug toxicity. SOR inhibits several human CYPs, including CYP2C8, which is a major enzyme in the elimination of oncology drugs like paclitaxel and imatinib. It has been reported that CYP2C8 inhibition by SOR in human liver microsomes is potentiated by NADPH-dependent biotransformation. This implicates a SOR metabolite in enhanced inhibition, although the identity of that metabolite is presently unclear. The present study evaluated the capacity of the major N-oxide metabolite of SOR (SNO) to inhibit CYP2C8-dependent paclitaxel 6α-hydroxylation. The IC₅₀ of SNO against CYP2C8 activity was found to be 3.7-fold lower than that for the parent drug (14 μM versus 51 μM). In molecular docking studies, both SOR and SNO interacted with active site residues in CYP2C8, but four additional major hydrogen and halogen bonding interactions were identified between SNO and amino acids in the B-B' loop region and helixes F' and I that comprise the catalytic region of the enzyme. In contrast, the binding of both SOR and SNO to active site residues in the closely related human CYP2C9 enzyme was similar, as were the IC₅₀s determined against CYP2C9-mediated losartan oxidation. These findings suggest that the active metabolite SNO could impair the elimination of coadministered drugs that are substrates for CYP2C8, and mediate toxic adverse events, perhaps in those individuals in whom SNO is formed extensively.

Human Family 1-4 cytochrome P450 enzymes involved in the metabolic activation of xenobiotic and physiological chemicals: an update

This is an overview of the metabolic activation of drugs, natural products, physiological compounds, and general chemicals by the catalytic activity of cytochrome P450 enzymes belonging to Families 1-4. The data were collected from > 5152 references. The total number of data entries of reactions catalyzed by P450s Families 1-4 was 7696 of which 1121 (~ 15%) were defined as bioactivation reactions of different degrees. The data were divided into groups of General Chemicals, Drugs, Natural Products, and Physiological Compounds, presented in tabular form. The metabolism and bioactivation of selected examples of each group are discussed. In most of the cases, the metabolites are directly toxic chemicals reacting with cell macromolecules, but in some cases the metabolites formed are not direct toxicants but participate as substrates in succeeding metabolic reactions (e.g., conjugation reactions), the products of which are final toxicants. We identified a high level of activation for three groups of compounds (General Chemicals, Drugs, and Natural Products) yielding activated metabolites and the generally low participation of Physiological Compounds in bioactivation reactions. In the group of General Chemicals, P450 enzymes 1A1, 1A2, and 1B1 dominate in the formation of activated metabolites. Drugs are mostly activated by the enzyme P450 3A4, and Natural Products by P450s 1A2, 2E1, and 3A4. Physiological Compounds showed no clearly dominant enzyme, but the highest numbers of activations are attributed to P450 1A, 1B1, and 3A enzymes. The results thus show, perhaps not surprisingly, that Physiological Compounds are infrequent substrates in bioactivation reactions catalyzed by P450 enzyme Families 1-4, with the exception of estrogens and arachidonic acid. The results thus provide information on the enzymes that activate specific groups of chemicals to toxic metabolites.

Case Study 9: Probe-Dependent Binding Explains Lack of CYP2C9 Inactivation by 1-Aminobenzotriazole (ABT)

The potential for new chemical entities to inhibit the major cytochrome P450 (CYP) isoforms is routinely evaluated to minimize the risk of developing drugs with drug-drug interaction liabilities. CYP inhibition assays are routinely performed in a high-throughput format to efficiently screen large numbers of compounds. In evaluating a time-saving assay using diclofenac as the CYP2C9 probe substrate, a discrepancy was observed in which minimal inhibition was detected using diclofenac whereas using (S)-warfarin resulted in potent inhibition, supporting the presence of dual-binding sites in the relatively large CYP2C9 active site cavity.These observations provided further insights into explaining the reported ineffective inactivation of CYP2C9 for the pan-CYP inactivator 1-aminobenzotriazole (ABT). Mechanistic reversible and time-dependent inhibition experiments revealed that the ineffective CYP2C9 inactivation by ABT was also probe-dependent, with utilization of (S)-warfarin as the probe substrate resulting in more potent CYP2C9 inhibition by ABT compared to diclofenac. Addition of (S)-warfarin to the reversible and time-dependent inhibition experiments between ABT and diclofenac resulted in an attenuation of the inhibitory effects of ABT on CYP2C9-mediated diclofenac metabolism. Molecular docking studies further confirmed that (S)-warfarin and diclofenac preferentially bind in different regions of the CYP2C9 active site, with (S)-warfarin occupying a distal "warfarin-binding pocket" and diclofenac occupying a binding site close to the active heme moiety. ABT preferentially binds in the distal warfarin-binding pocket, supporting that diclofenac is minimally deterred from access to the CYP2C9 active site in the presence of ABT, thus resulting in minimal inactivation. Simultaneously docking of (S)-warfarin and ABT revealed that (S)-warfarin outcompetes ABT for the distal binding site and results in the binding of ABT to the CYP2C9 active site, supporting the observations of potent inactivation of CYP2C9 when (S)-warfarin is the probe substrate.These results highlight that probe selection is crucial when evaluating CYP inhibition potential, and it is recommended that multiple probes be utilized for CYP2C9, similar to the approach routinely employed for CYP3A4. Further, utilization of ABT as a pan-inhibitor of CYP activity for investigational compounds, both in vitro and in vivo, should be accompanied with the understanding that residual CYP-mediated oxidative metabolism could potentially be observed for CYP2C9 substrates and should not necessarily be attributed to non-P450-mediated metabolism.

Exploring the biotransformation of N-(2-hydroxyphenyl)-2-propylpentanamide (an aryl valproic acid derivative) by CYP2C11, using in silico predictions and in vitro studies

OBJECTIVES

N-(2-hydroxyphenyl)-2-propylpentanamide (HO-AAVPA), a derivative of valproic acid (VPA), has been proposed as a potential anticancer agent due to its improved antiproliferative effects in some cancer cell lines. Although there is evidence that VPA is metabolized by cytochrome P450 2C11 rat isoform, HO-AAVPA CYP-mediated metabolism has not yet been fully explored. Therefore, in this work, the biotransformation of HO-AAVPA by CYP2C11 was investigated.

METHODS

Kinetic parameters and spectral interaction between HO-AAVPA and CYP were evaluated using rat liver microsomes. The participation of CYP2C11 in metabolism of HO-AAVPA was confirmed by cimetidine (CIM) inhibition assay. Docking and molecular dynamics simulations coupled to MMGBSA methods were used in theoretical study.

KEY FINDINGS

HO-AAVPA is metabolized by CYP enzymes (K_M  = 38.94 µm), yielding a hydroxylated metabolite according to its HPLC retention time (5.4 min) and MS analysis (252.2 m/z). In addition, CIM inhibition in rat liver microsomes (K_i  = 59.23 µm) confirmed that CYP2C11 is mainly involved in HO-AAVPA metabolism. Furthermore, HO-AAVPA interacts with CYP2C11 as a type I ligand. HO-AAVPA is stabilized at the CYP2C11 ligand recognition site through a map of interactions similar to other typical CYP2C11 substrates.

CONCLUSION

Therefore, rat liver CYP2C11 isoform is able to metabolize HO-AAVPA.

Formation of a Toxic Quinoneimine Metabolite from Diclofenac: A Quantum Chemical Study

BACKGROUND

Diclofenac is a non-steroidal antiinflammatory drug. It is predominantly metabolized by CYP2C9. 4'-hydroxydiclofenac and its quinoneimine are the metabolites of diclofenac. However, few numbers of serious cases of idiosyncratic hepatotoxicity due to diclofenac metabolism were reported. The formation of the quinoneimine metabolite was found to be responsible for this idiosyncratic toxicity. Quinoneimine is an over-oxidized metabolite of diclofenac.

METHOD

In this work, computational studies were conducted to detail the formation of a quinoneimine metabolite from diclofenac. Further, the idiosyncratic toxicity of quinoneimine due to its reactivity was also investigated by quantum chemical analysis.

RESULTS & CONCLUSION

The results demonstrate the possibility of formation of quinoneimine metabolite due to various factors that are involved in the metabolism of diclofenac. The present study may provide the structural in-sights during the drug development processes to avoid the metabolism directed idiosyncratic toxicity.

Insights into molecular mechanisms of drug metabolism dysfunction of human CYP2C9*30

Cytochrome P450 2C9 (CYP2C9) metabolizes about 15% of clinically administrated drugs. The allelic variant CYP2C9*30 (A477T) is associated to diminished response to the antihypertensive effects of the prodrug losartan and affected metabolism of other drugs. Here, we investigated molecular mechanisms involved in the functional consequences of this amino-acid substitution. Molecular dynamics (MD) simulations performed for the active species of the enzyme (heme in the Compound I state), in the apo or substrate-bound state, and binding energy analyses gave insights into altered protein structure and dynamics involved in the defective drug metabolism of human CYP2C9.30. Our data revealed an increased rigidity of the key Substrate Recognition Sites SRS1 and SRS5 and shifting of the β turn 4 of SRS6 toward the helix F in CYP2C9.30. Channel and binding substrate dynamics analyses showed altered substrate channel access and active site accommodation. These conformational and dynamic changes are believed to be involved in the governing mechanism of the reduced catalytic activity. An ensemble of representative conformations of the WT and A477T mutant properly accommodating drug substrates were identified, those structures can be used for prediction of new CYP2C9 and CYP2C9.30 substrates and drug-drug interactions.

Pharmacogenomics of CYP2C9: Functional and Clinical Considerations

CYP2C9 is the most abundant CYP2C subfamily enzyme in human liver and the most important contributor from this subfamily to drug metabolism. Polymorphisms resulting in decreased enzyme activity are common in the CYP2C9 gene and this, combined with narrow therapeutic indices for several key drug substrates, results in some important issues relating to drug safety and efficacy. CYP2C9 substrate selectivity is detailed and, based on crystal structures for the enzyme, we describe how CYP2C9 catalyzes these reactions. Factors relevant to clinical response to CYP2C9 substrates including inhibition, induction and genetic polymorphism are discussed in detail. In particular, we consider the issue of ethnic variation in pattern and frequency of genetic polymorphisms and clinical implications. Warfarin is the most well studied CYP2C9 substrate; recent work on use of dosing algorithms that include CYP2C9 genotype to improve patient safety during initiation of warfarin dosing are reviewed and prospects for their clinical implementation considered. Finally, we discuss a novel approach to cataloging the functional capabilities of rare 'variants of uncertain significance', which are increasingly detected as more exome and genome sequencing of diverse populations is conducted.

Structural Basis of Single-Nucleotide Polymorphisms in Cytochrome P450 2C9

Single-nucleotide polymorphisms in drug-metabolizing cytochrome P450 (CYP) enzymes are important contributors to interindividual differences in drug metabolism leading to adverse drug reactions. Despite their extensive characterization and importance in pharmacogenetics of clinical drugs, the structural basis of CYP polymorphisms has remained scant. Here we report the crystal structures of human CYP2C9 and its polymorphic variants, *3 (I359L) and *30 (A477T), with an antihypertensive drug losartan. The structures show distinct interaction and occupation of losartan in the active site, the access channel, and the peripheral binding site. The I359L substitution located far from the active site remarkably altered the residue side chains near the active site and the access channel, whereas the T477 substitution illustrated hydrogen-bonding interaction with the reoriented side chain of Q214. The results yield structural insights into the reduced catalytic activity of the CYP2C9 variants and have important implications for understanding genetic polymorphisms in CYP-mediated drug metabolism.

Single-Protein Tracking Reveals That NADPH Mediates the Insertion of Cytochrome P450 Reductase into a Biomimetic of the Endoplasmic Reticulum

Cytochrome P450 reductase (CPR) is the redox partner for most human cytochrome P450 enzymes. It is also believed that CPR is an integral membrane protein exclusively. Herein, we report that, contrary to this belief, CPR can exist as a peripheral membrane protein in the absence of NADPH and will transition to an integral membrane protein in the presence of stoichiometric amounts of NADPH or greater. All experiments were performed in a solid-supported cushioned lipid bilayer that closely matched the chemical composition of the human endoplasmic reticulum and served as an ER biomimetic. The phase characteristics and fluidity of the ER biomimetic was characterized with fluorescence micrographs and temperature-dependent fluorescence recovery after photobleaching. The interactions of CPR with the ER biomimetic were directly observed by tracking single CPR molecules using time-lapse single-molecule fluorescence imaging and subsequent analysis of tracks. These studies revealed dramatic changes in diffusion coefficient and the degree of partitioning of CPR as a function of NADPH concentration.

Toxicity Originating from Thiophene Containing Drugs: Exploring the Mechanism using Quantum Chemical Methods

Drug metabolism of thiophene containing substrates by cytochrome P450s (CYP450) leads to toxic side effects, for example, nephrotoxicity (suprofen, ticlopidine), hepatotoxicity (tienilic acid), thrombotic thrombocytopenic purpura (clopidogrel), and aplastic anemia (ticlopidine). The origin of toxicity in these cases has been attributed to two different CYP450 mediated metabolic reactions: S-oxidation and epoxidation. In this work, the molecular level details of the bioinorganic chemistry associated with the generation of these competitive reactions are reported. Density functional theory was utilized (i) to explore the molecular mechanism for S-oxidation and epoxidation using the radical cationic center Cpd I [(iron(IV)-oxo-heme porphine system with SH(-) as the axial ligand, to mimic CYP450s] as the model oxidant, (ii) to establish the 3D structures of the reactants, transition states, and products on both the metabolic pathways, and (iii) to examine the potential energy (PE) profile for both the pathways to determine the energetically preferred toxic metabolite formation. The energy barrier required for S-oxidation was observed to be 14.75 kcal/mol as compared to that of the epoxidation reaction (13.23 kcal/mol) on the doublet PE surface of Cpd I. The formation of the epoxide metabolite was found to be highly exothermic (-23.24 kcal/mol), as compared to S-oxidation (-8.08 kcal/mol). Hence, on a relative scale the epoxidation process was observed to be thermodynamically and kinetically more favorable. The energy profiles associated with the reactions of the S-oxide and epoxide toxic metabolites were also explored. This study helps in understanding the CYP450-catalyzed toxic reactions of drugs containing the thiophene ring at the atomic level.

Mechanistic studies of the cationic binding pocket of CYP2C9 in vitro and in silico: metabolism of nonionizable analogs of tienilic acid

Tienilic acid (TA) is selectively oxidized at the C-5 position of the thiophene ring by the human liver enzyme cytochrome P450 2C9 (CYP2C9). This oxidation is mediated by the proximal positioning of the thiophene over the heme iron, which is proposed to be coordinated by an interaction of the TA carboxylic acid to a cationic binding pocket in the enzyme active site. In this study, we investigated how chemical modification of TA influences the bioactivation by CYP2C9. For this investigation, nine analogs of TA were chosen with substitutions on either side of the molecule. We tested three parameters, including CYP2C9 inhibition, metabolic profiling, and in silico docking. Of the 10 compounds tested, only two (TA and a noncarboxyl analog) resulted in competitive and time-dependent inhibition of CYP2C9. Metabolic profiling revealed a trend in which substitution of the carboxylate with nonionizable functional groups resulted in metabolic switching from oxidation of the aromatic ring to dealkylation reactions at the opposite side of the structure. The in silico modeling predicted an opposite binding orientation to that of TA for many analogs, including the 3-thenoyl regio-isomer analog, which contradicts previous models. Together these data show that disrupting interactions with the cationic binding pocket of CYP2C9 will impact the sites of metabolism and inhibition of the enzyme.

Bioactivation potential of thiophene-containing drugs

Thiophene is a five-membered, sulfur-containing heteroaromatic ring commonly used as a building block in drugs. It is considered to be a structural alert, as its metabolism can lead to the formation of reactive metabolites. Thiophene S-oxides and thiophene epoxides are highly reactive electrophilic thiophene metabolites whose formation is cytochrome P450-dependent. These reactive thiophene-based metabolites are quite often responsible for drug-induced hepatotoxicity. Tienilic acid is an example of a thiophene-based drug that was withdrawn from the market after only a few months of use, due to severe cases of immune hepatitis. However, inclusion of the thiophene moiety in drugs does not necessarily result in toxic effects. The presence of other, less toxic metabolic pathways, as well as an effective detoxification system in our body, protects us from the bioactivation potential of the thiophene ring. Thus, the presence of a structural alert itself is insufficient to predict a compound's toxicity. The question therefore arises as to which factors significantly influence the toxicity of thiophene-containing drugs. There is no easy way to answer this question. However, the findings presented here indicate that, for a number of reasons, daily dose and alternative metabolic pathways are important factors when predicting toxicity and will therefore be discussed together with examples.

Structural and energetic analysis to provide insight residues of CYP2C9, 2C11 and 2E1 involved in valproic acid dehydrogenation selectivity

Docking and molecular dynamics (MD) simulation have been two computational techniques used to gain insight about the substrate orientation within protein active sites, allowing to identify potential residues involved in the binding and catalytic mechanisms. In this study, both methods were combined to predict the regioselectivity in the binding mode of valproic acid (VPA) on three cytochrome P-450 (CYP) isoforms CYP2C9, CYP2C11, and CYP2E1, which are involved in the biotransformation of VPA yielding reactive hepatotoxic intermediate 2-n-propyl-4-pentenoic acid (4nVPA). There are experimental data about hydrogen atom abstraction of the C4-position of VPA to yield 4nVPA, however, there are not structural evidence about the binding mode of VPA and 4nVPA on CYPs. Therefore, the complexes between these CYP isoforms and VPA or 4nVPA were studied to explore their differences in binding and energetic stabilization. Docking results showed that VPA and 4nVPA are coupled into CYPs binding site in a similar conformation, but it does not explain the VPA hydrogen atom abstraction. On the other hand, MD simulations showed a set of energetic states that reorient VPA at the first ns, then making it susceptible to a dehydrogenation reaction. For 4nVPA, multiple binding modes were observed in which the different states could favor either undergo other reaction mechanism or ligand expulsion from the binding site. Otherwise, the energetic and entropic contribution point out a similar behavior for the three CYP complexes, showing as expected a more energetically favorable binding free energy for the complexes between CYPs and VPA than with 4nVPA.

High warfarin sensitivity in carriers of CYP2C9*35 is determined by the impaired interaction with P450 oxidoreductase

Cytochrome P450 2C9 (CYP2C9) metabolizes many clinically important drugs including warfarin and diclofenac. We have recently reported a new allelic variant, CYP2C9*35, found in a warfarin hypersensitive patient with Arg125Leu and Arg144Cys mutations. Here, we have investigated the molecular basis for the functional consequences of these polymorphic changes. CYP2C9.1 and CYP2C9-Arg144Cys expressed in human embryonic kidney 293 cells effectively metabolized both S-warfarin and diclofenac in NADPH-dependent reactions, whereas CYP2C9-Arg125Leu or CYP2C9.35 were catalytically silent. However, when NADPH was replaced by a direct electron donor to CYPs, cumene hydroperoxide, hereby bypassing the CYP oxidoreductase (POR), all variant enzymes were active, indicating unproductive interactions between CYP2C9.35 and POR. In silico analysis revealed a decrease of the electrostatic potential of CYP2C9-Arg125Leu-POR interacting surface and the loss of stabilizing salt bridges between these proteins. In conclusion, our data strongly suggest that the Arg125Leu mutation in CYP2C9.35 prevents CYP2C9-POR interactions resulting in the absence of NADPH-dependent CYP2C9-catalyzed activity in vivo, thus influencing the warfarin sensitivity in the carriers of this allele.

The dynamics of camphor in the cytochrome P450 CYP101D2

The recent crystal structures of CYP101D2, a cytochrome P450 protein from the oligotrophic bacterium Novosphingobium aromaticivorans DSM12444 revealed that both the native (substrate-free) and camphor-soaked forms have open conformations. Furthermore, two other potential camphor-binding sites were also identified from electron densities in the camphor-soaked structure, one being located in the access channel and the other in a cavity on the surface near the F-helix side of the F-G loop termed the substrate recognition site. These latter sites may be key intermediate positions on the pathway for substrate access to or product egress from the active site. Here, we show via the use of unbiased atomistic molecular dynamics simulations that despite the open conformation of the native and camphor-bound crystal structures, the underlying dynamics of CYP101D2 appear to be very similar to other CYP proteins. Simulations of the native structure demonstrated that the protein is capable of sampling many different conformational substates. At the same time, simulations with the camphor positioned at various locations within the access channel or recognition site show that movement towards the active site or towards bulk solvent can readily occur on a short timescale, thus confirming many previously reported in silico studies using steered molecular dynamics. The simulations also demonstrate how the fluctuations of an aromatic gate appear to control access to the active site. Finally, comparison of camphor-bound simulations with the native simulations suggests that the fluctuations can be of similar level and thus are more representative of the conformational selection model rather than induced fit.

[Metabolism of xenobiotics: beneficial and adverse effects]

The systems developed by living organisms for the metabolism of xenobiotics play a key role in the adaptation of living species to their chemical environment. Recent data about mammalian cytochrome P450 structures have led to a better understanding of the molecular basis for the adaptability of these enzymes to xenobiotics exhibiting highly variable structures. The action of these enzymes on xenobiotics leads to other beneficial effects such as the bioactivation of some drugs, but also to adverse effects with the formation of aggressive metabolites for the cell that are responsible for the appearance of many toxicities.

Distal effect of amino acid substitutions in CYP2C9 polymorphic variants causes differences in interatomic interactions against (S)-warfarin

Cytochrome P450 2C9 (CYP2C9) is crucial in excretion of commonly prescribed drugs. However, changes in metabolic activity caused by CYP2C9 polymorphisms inevitably result in adverse drug effects. CYP2C9*2 and *3 are prevalent in Caucasian populations whereas CYP2C9*13 is remarkable in Asian populations. Single amino acid substitutions caused by these mutations are located outside catalytic cavity but affect kinetic activities of mutants compared to wild-type enzyme. To relate distal effects of these mutations and defective drug metabolisms, simulations of CYP2C9 binding to anti-coagulant (S)-warfarin were performed as a system model. Representative (S)-warfarin-bound forms of wild-type and mutants were sorted and assessed through knowledge-based scoring function. Interatomic interactions towards (S)-warfarin were predicted to be less favorable in mutant structures in correlation with larger distance between hydroxylation site of (S)-warfarin and reactive oxyferryl heme than wild-type structure. Using computational approach could delineate complication of CYP polymorphism in management of drug therapy.

MOLECULAR DOCKING STUDY OF THE AFFINITY OF CYP2C9 AND CYP2D6 FOR IMRECOXIB

With the aid of the automatic molecular docking, the affinity of CYP2C9 and CYP2D6 for imrecoxib was studied by InsightII/Affinity program. The results indicate that CYP2C9–imrecoxib complex has higher stability and stronger affinity because CYP2C9 has more favorable interaction energy (-62.72 kcal/mol) and higher Ludi score (610) with imrecoxib than CYP2D6 (-50.22 kcal/mol and 551) and this is consistent with the results of the kinetic experiments by Li et al. By analyzing the theoretical results combined with the experimental ones, we suggest that the affinity difference is caused by the difference of the structure between CYP2C9 and CYP2D6, and the most important residues for enzyme–substrate complexes, such as Phe476, Asn204, Phe100, Leu366 and Arg108 of CYP2C9 and Phe120, Glu216, and Phe483 of CYP2D6 were also identified.

Probing small-molecule binding to cytochrome P450 2D6 and 2C9: An in silico protocol for generating toxicity alerts

Drug metabolism, toxicity, and their interaction profiles are major issues in the drug-discovery and lead-optimization processes. The cytochromes P450 (CYPs) 2D6 and 2C9 are enzymes involved in the oxidative metabolism of a majority of marketed drugs. Therefore, the prediction of the binding affinity towards CYP2D6 and CYP2C9 would be beneficial for identifying cytochrome-mediated adverse effects triggered by drugs or chemicals (e.g., toxic reactions, drug-drug, and food-drug interactions). By identifying the binding mode by using pharmacophore prealignment, automated flexible docking, and by quantifying the binding affinity by multidimensional QSAR (mQSAR), we validated a model family of 56 compounds (46 training, 10 test) and 85 compounds (68 training, 17 test) for CYP2D6 and CYP2C9, respectively. The correlation with the experimental data (cross-validated r²=0.811 for CYP2D6 and 0.687 for CYP2C9) suggests that our approach is suited for predicting the binding affinity of compounds towards CYP2D6 and CYP2C9. The models were challenged by Y-scrambling and by testing an external dataset of binding compounds (15 compounds for CYP2D6 and 40 for CYP2C9). To assess the probability of false-positive predictions, datasets of nonbinders (64 compounds for CYP2D6 and 56 for CYP2C9) were tested by using the same protocol. The two validated mQSAR models were subsequently added to the VirtualToxLab (VTL, http://www.virtualtoxlab.org).

New insights into the structural characteristics and functional relevance of the human cytochrome P450 2D6 enzyme

To date, the crystal structures of at least 12 human CYPs (1A2, 2A6, 2A13, 2C8, 2C9, 2D6, 2E1, 2R1, 3A4, 7A1, 8A1, and 46A1) have been determined. CYP2D6 accounts for only a small percentage of all hepatic CYPs (< 2%), but it metabolizes approximately 25% of clinically used drugs with significant polymorphisms. CYP2D6 also metabolizes procarcinogens and neurotoxins, such as 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine, 1,2,3,4-tetrahydroquinoline, and indolealkylamines. Moreover, the enzyme utilizes hydroxytryptamines and neurosteroids as endogenous substrates. Typical CYP2D6 substrates are usually lipophilic bases with an aromatic ring and a nitrogen atom, which can be protonated at physiological pH. Substrate binding is generally followed by oxidation (5-7 A) from the proposed nitrogen-Asp301 interaction. A number of homology models have been constructed to explore the structural features of CYP2D6, while antibody studies also provide useful structural information. Site-directed mutagenesis studies have demonstrated that Glu216, Asp301, Phe120, Phe481, and Phe483 play important roles in determining the binding of ligands to CYP2D6. The structure of human CYP2D6 has been recently determined and shows the characteristic CYP fold observed for other members of the CYP superfamily. The lengths and orientations of the individual secondary structural elements in the CYP2D6 structure are similar to those seen in other human CYP2 members, such as CYP2C9 and 2C8. The 2D6 structure has a well-defined active-site cavity located above the heme group with a volume of approximately 540 A(3), which is larger than equivalent cavities in CYP2A6 (260 A(3)), 1A2 (375 A(3)), and 2E1 (190 A(3)), but smaller than those in CYP3A4 (1385 A(3)) and 2C8 (1438 A(3)). Further studies are required to delineate the molecular mechanisms involved in CYP2D6 ligand interactions and their implications for drug development and clinical practice.

Testosterone-metabolizing capacity and characteristics of adrenal microsomes in human fetus in vitro

We describe the metabolite spectrum of testosterone (T) in human fetal adrenal in vitro, and verify possible roles of CYP3A and 2C isoforms of human fetal adrenal in T metabolism which respond to T metabolism in liver. Isolation and identification of T and its metabolites were carried out by isocratic high-performance liquid chromatography. CYP isoforms involved in T metabolism were examined by an inhibition study. Four metabolites formed by fetal adrenal microsomes were androstenedione (A), 16alpha-hydroxytestosterone (16alpha-HT), 16beta-hydroxytestosterone (16beta-HT) and 6beta-hydroxytestosterone (6beta-HT), which made up 61%, 25%, 12% and 2%, respectively, of the total. However, fetal liver microsomes produced six metabolites, which were, from high to low: A, 2alpha-HT, 6beta-HT, 2beta-HT, 16beta-HT and 16alpha-HT. Generation of A, 16alpha-HT and 16beta-HT in the adrenals was much greater than that in the liver. Erythromycin inhibited the formation of 6beta-HT and 16beta-HT with maximal inhibition of 76% and 47%. Sulfaphenazole and omeprazone had no inhibitive effect on the formations of T metabolites. The predominant metabolites of T in fetal adrenal, A and 16alpha-HT, were not decreased by these three inhibitors. These results suggest that the human fetal adrenal has greater T-metabolizing ability and a different metabolizing pathway from adult and fetal liver. Some other CYPs but not CYP3A7, CYP2C9 and CYP2C19 might play the most important part in the metabolism of T in human fetal adrenal although CYP3A7 is responsible for the formation of 6beta-HT and partly participates in the formation of 16beta-HT.

Molecular modeling of human cytochrome P450 2W1 and its interactions with substrates

The human cytochrome P450 2W1 (CYP2W1) was categorized into the so-called "orphan" CYPs because of its unknown enzymatic function. However, recent studies showed that the recombinant CYP2W1 exhibited broad catalytic activity towards several chemicals. Furthermore, this enzyme was selectively expressed in some forms of cancers, whereas a very low expression was found in human normal issues. These render CYP2W1 as a potential drug target for cancer therapy. At present, however, little information is available on the active site topology and the substrate binding modes of CYP2W1. In this study, the three-dimensional model of CYP2W1 was constructed using the homology modeling method. Two known substrates, benzphetamine and indole, were then docked into the active site, and refined by molecular dynamics simulations. The interaction energy between the substrates and the enzyme was calculated and analyzed by using the MM-GBSA method. The results indicated that the constructed CYP2W1 model can account for the regioselectivity of this enzyme towards the known substrates and van der Waals interactions were the driving force for the substrate binding. Several key residues were identified to be responsible for the binding of indole and benzphetamine with CYP2W1. These findings provide useful information for the detailed characterization of the biological roles of CYP2W1 and structure-based drug design of this enzyme.

A comparative analysis of binding sites between mouse CYP2C38 and CYP2C39 based on homology modeling, molecular dynamics simulation and docking studies

Mouse CYP2C38 and CYP2C39 are two closely related enzymes with 91.8% sequence identity. But they exhibit different substrate binding features. In this study, three-dimensional models of CYP2C38 and CYP2C39 were constructed using X-ray crystal structure of human CYP2C8 as the template based on homology modeling methods and molecular dynamics simulations. Tolbutamide as the common substrate of CYP2C38 and CYP2C39 was docked into them and positioned in their active sites with different orientation. All-trans retinoic acid (atRA) is a specific substrate for CYP2C39 and not catalyzed by CYP2C38. By comparison of active site architectures between CYP2C38 and CYP2C39, the possible reasons affecting their substrate binding were proposed. In addition, Arg241, Glu300, Leu366 and Leu476 are identified as critical residue for substrates binding.

Identification of selectivity-determining residues in cytochrome P450 monooxygenases: a systematic analysis of the substrate recognition site 5

The large and diverse family of cytochrome P450 monooxygenases (CYPs) was systematically analyzed to identify selectivity- and specificity-determining residues in the substrate recognition site 5, which is located in close vicinity to the heme center. A positively charged heme-interacting residue was identified in the structures of 29 monooxygenases and in 97.7% of the 6379 CYP sequences investigated here. This heme-interacting residue restricts the conformation of the substrate recognition site 5 and is preferentially located at position 10 or 11 after the conserved ExxR motif (in 94.4% of the sequences), in 3.3% of the sequences at position 9 or 12. As a result, a classification by the position of the heme-interacting residue allows to predict residues that are closest to the heme center and restrict its accessibility. In 98.4% of all CYP sequences a preferentially hydrophobic residue is located at position 5 after the ExxR motif that is predicted to point close to the heme center. Replacing this residue by hydrophobic residues of different size has been shown to change substrate specificity and regioselectivity for CYPs of different superfamilies. Twenty-seven percent of all CYPs are predicted to contain a second selectivity-determining residue at position 9 after the ExxR motif that can be identified by the pattern EXXR-X(7)-{P}-x-P-[HKR].

CYP2C9 amino acid residues influencing phenytoin turnover and metabolite regio- and stereochemistry

Phenytoin has been an effective anticonvulsant agent for over 60 years, although its clinical use is complicated by nonlinear pharmacokinetics, a narrow therapeutic index, and metabolically based drug-drug interactions. Although it is well established that CYP2C9 is the major cytochrome P450 enzyme controlling metabolic elimination of phenytoin through its oxidative conversion to (S)-5-(4-hydroxyphenyl)-5-phenylhydantoin (p-HPPH), nothing is known about the amino acid binding determinants within the CYP2C9 active site that promote metabolism and maintain the tight stereocontrol of hydroxy metabolite formation. This knowledge gap was addressed here through the construction of nine active site mutants at amino acid positions Phe100, Arg108, Phe114, Leu208, and Phe476 and in vitro analysis of the steady-state kinetics and stereochemistry of p-HPPH formation. The F100L and F114W mutants exhibited 4- to 5-fold increases in catalytic efficiency, whereas the F100W, F114L, F476L, and F476W mutants lost >90% of their phenytoin hydroxylation capacity. This pattern of effects differs substantially from that found previously for (S)-warfarin and (S)-flurbiprofen metabolism, suggesting that these three ligands bind within discrete locations in the CYP2C9 active site. Only the F114L, F476L, and L208V mutants altered phenytoin's orientation during catalytic turnover. The L208V mutant also uniquely demonstrated enhanced 6-hydroxylation of (S)-warfarin. These latter data provide the first experimental evidence for a role of the F-G loop region in dictating the catalytic orientation of substrates within the CYP2C9 active site.

Functional analysis of phenylalanine residues in the active site of cytochrome P450 2C9

The two published crystal structures of cytochrome P450 2C9, complexed with ( S)-warfarin or flurbiprofen, implicate a cluster of three active site phenylalanine residues (F100, F114, F476) in ligand binding. However, these three residues appear to interact differently with these two ligands based on the static crystal structures. To elucidate the importance of CYP2C9's active site phenylalanines on substrate binding, orientation, and catalytic turnover, a series of leucine and tryptophan mutants were constructed and their interactions with ( S)-warfarin and ( S)-flurbiprofen examined. The F100-->L mutation had minor effects on substrate binding and metabolism of each substrate. In contrast, the F114L and F476L mutants exhibited substantially reduced ( S)-warfarin metabolism and altered hydroxy metabolite profiles but only modestly decreased nonsteroidal antiinflammatory drug (NSAID) turnover while maintaining product regioselectivity. The F114-->W and F476-->W mutations also had opposing effects on ( S)-warfarin versus NSAID turnover. Notably, the F476W mutant increased the efficiency of ( S)-warfarin metabolism 5-fold, yet decreased the efficiency of ( S)-flurbiprofen turnover 20-fold. (1)H NMR T 1 relaxation studies suggested a slightly closer positioning of ( S)-warfarin to the heme in the F476W mutant relative to the wild-type enzyme, and stoichiometry studies indicated enhanced coupling of reducing equivalents to product formation for ( S)-warfarin, again in contrast to effects observed with ( S)-flurbiprofen. These data demonstrate that F114 and F476, but not F100, influence ( S)-warfarin's catalytic orientation. Differential interactions of F476 mutants with the two substrates suggest that their catalytically productive binding modes are not superimposable.

Mechanism-based inactivation of cytochrome P450 2C9 by tienilic acid and (+/-)-suprofen: a comparison of kinetics and probe substrate selection

In vitro experiments were conducted to compare k(inact), K(I) and inactivation efficiency (k(inact)/K(I)) of cytochrome P450 (P450) 2C9 by tienilic acid and (+/-)-suprofen using (S)-flurbiprofen, diclofenac, and (S)-warfarin as reporter substrates. Although the inactivation of P450 2C9 by tienilic acid when (S)-flurbiprofen and diclofenac were used as substrates was similar (efficiency of approximately 9 ml/min/micromol), the inactivation kinetics were characterized by a sigmoidal profile. (+/-)-Suprofen inactivation of (S)-flurbiprofen and diclofenac hydroxylation was also described by a sigmoidal profile, although inactivation was markedly less efficient (approximately 1 ml/min/micromol). In contrast, inactivation of P450 2C9-mediated (S)-warfarin 7-hydroxylation by tienilic acid and (+/-)-suprofen was best fit to a hyperbolic equation, where inactivation efficiency was moderately higher (10 ml/min/micromol) and approximately 3-fold higher (3 ml/min/micromol), respectively, relative to that of the other probe substrates, which argues for careful consideration of reporter substrate when mechanism-based inactivation of P450 2C9 is assessed in vitro. Further investigations into the increased inactivation seen with tienilic acid relative to that with (+/-)-suprofen revealed that tienilic acid is a higher affinity substrate with a spectral binding affinity constant (K(s)) of 2 microM and an in vitro half-life of 5 min compared with a K(s) of 21 microM and a 50 min in vitro half-life for (+/-)-suprofen. Lastly, a close analog of tienilic acid with the carboxylate functionality replaced by an oxirane ring was devoid of inactivation properties, which suggests that an ionic binding interaction with a positively charged residue in the P450 2C9 active site is critical for recognition and mechanism-based inactivation by these close structural analogs.

Re-engineering of CYP2C9 to probe acid-base substrate selectivity

A common feature of many CYP2C9 ligands is their weak acidity. As revealed by crystallography, the structural basis for this behavior involves a charge-pairing interaction between an anionic moiety on the substrate and an active site R108 residue. In the present study we attempted to re-engineer CYP2C9 to better accept basic ligands by charge reversal at this key residue. We expressed and purified the R108E and R108E/D293N mutants and compared their ability with that of native CYP2C9 to interact with (S)-warfarin, diclofenac, pyrene, propranolol, and ibuprofen amine. As expected, the R108E mutant maintained all the native enzyme's pyrene 1-hydroxylation activity, but catalytic activity toward diclofenac and (S)-warfarin was abrogated. In contrast, the double mutant displayed much less selectivity in its behavior toward these control ligands. Neither of the mutants displayed significant enhancement of propranolol metabolism, and all three preparations exhibited a type II (inhibitor) rather than type I (substrate) spectrum with ibuprofen amine, although binding became progressively weaker with the single and double mutants. Collectively, these data underscore the importance of the amino acid at position 108 in the acid substrate selectivity of CYP2C9, highlight the accommodating nature of the CYP2C9 active site, and provide a cautionary note regarding facile re-engineering of these complex cytochrome P450 active sites.

Defining targets for investigating the pharmacogenomics of adverse drug reactions to antifungal agents

Adverse drug reactions (ADRs) associated with antifungal therapy are major problems in patients with invasive fungal infections. Whether by clinical history or patterns of genetic variation, the identification of patients at risk for ADRs should result in improved outcomes while minimizing deleterious side effects. A major contributing factor to ADRs with antifungal agents relates to drug distribution, metabolism and excretion. Genetic variation in key genes can alter the structure and expression of genes and gene products (e.g., proteins). Thus far, the effort has focused on identifying polymorphisms with either empirical or predicted in silico functional consequences; the best candidate genes encode phase I and II drug-metabolizing enzymes (e.g., CYP2C19 and N-acetyltransferase), plasma proteins (albumin and lipoproteins) and drug transporters (P-glycoprotein and multidrug resistance proteins), which can affect the disposition of antifungal agents, eventually leading to dose-dependent (type A) toxicity. Less is known regarding the key genes that interact with antifungal agents, resulting in idiosyncratic (type B) ADRs. The possible role of certain gene products and genetic polymorphisms in the toxicities of antifungal agents are discussed in this review. The preliminary data address the following: low-density lipoproteins and cholesteryl ester transfer protein in amphotericin B renal toxicity; toll-like receptor 1 and 2 in amphotericin B infusion-related ADRs; phosphodiesterase 6 in voriconazole visual adverse events; flavin-containing monooxygenase, glutathione transferases and multidrug resistance proteins 1 and 2 in ketoconazole and terbinafine hepatotoxicity; CYP enzymes and P-glycoprotein in drug interactions between azoles and coadministered medications; multidrug resistance proteins 8 and 9 on 5-flucytosine bone marrow toxicity; and mast cell activation in caspofungin histamine release. This will focus on high-priority candidate genes, which could provide a starting point for molecular studies to elucidate the potential mechanisms for understanding toxicity associated with antifungal drugs as well as identifying candidate genes for large population prospective genetic association studies.

Improvement of cyclophosphamide activation by CYP2B6 mutants: from in silico to ex vivo

Cyclophosphamide (CPA) is a chemotherapeutic agent that is primarily activated in the liver by cytochrome P4502B6 (CYP2B6) and then transported to the tumor via blood flow. To prevent deleterious secondary effects, P450-based gene-directed enzyme prodrug therapy (GDEPT) consists of expressing CYP2B6 in tumor cells before CPA treatment. Given the relatively low affinity of CYP2B6 for CPA, the aim of our work was to modify CYP2B6 to increase its catalytic efficiency (V(max)/K(m)) to metabolize CPA into 4'-OH CPA. A molecular model of CYP2B6 was built, and four residues in close contact with the substrate were subjected to mutagenesis. Canine CYP2B11 exhibiting a particularly low K(m) to CPA, the amino acids exclusively present in the CYP2B11 substrate recognition sequences were substituted in human CYP2B6. All mutants (n = 26) were expressed in Saccharomyces cerevisiae and their enzymatic constants (K(m), V(max)) evaluated using CPA as substrate. Five mutants exhibited a 2- to 3-fold higher catalytic efficiency than wild-type CYP2B6. A double mutant, comprising the two most effective mutations, showed a 4-fold increase in K(m)/V(max). Molecular dynamic simulations of several mutants were found to be consistent with the observed modifications in catalytic efficiency. Finally, expression of the CYP2B6 114V/477W double mutant, contrary to wt CYP2B6, allowed switching of a resistant human head and neck cancer cell line (A-253) into a sensitive cell line toward CPA. Thus, we were able to obtain a new efficient CYP2B6 mutant able to metabolize CPA, an important step in the GDEPT strategy for human cancer treatment.

Stereoselective hexobarbital 3'-hydroxylation by CYP2C19 expressed in yeast cells and the roles of amino acid residues at positions 300 and 476

We examined the enzymatic function of recombinant CYP2C19 in enantiomeric hexobarbital (HB) 3'-hydroxylation, and searched the roles of amino acid residues, such as Phe-100, Phe-114, Asp-293, Glu-300, and Phe-476 of CYP2C19 in the stereoselective HB 3'-hydroxylation, using a yeast cell expression system and site-directed mutagenesis method. CYP2C19 wild-type exerted substrate enantioselectivity of (R)-HB>>(S)-HB and metabolite diastereoselectivity of 3'(R)<3'(S) in 3'-hydroxylation of HB enantiomers. The substitution of Asp-293 by alanine failed to yield an observable peak at 450 nm in its reduced carbon monoxide-difference spectrum. CYP2C19-E300A and CYP2C19-E300V with alanine and valine, respectively, in place of Glu-300 exerted total HB 3'-hydroxylation activities of 45 and 108%, respectively, that of the wild-type. Interestingly, these two mutants showed substrate enantioselectivity of (R)-HB<(S)-HB, which is opposite to that of the wild-type, while metabolite diasteroselectivity remained unchanged. The replacement of Phe-476 by alanine increased total HB 3'-hydroxylation activity to approximately 3-fold that of the wild-type. Particularly, 3'(S)-OH-(S)-HB-forming activity elevated to 7-fold that of the wild-type, resulting in the reversal of the substrate enantioselectivity. In contrast, the substitution of phenylalanine at positions 100 and 114 by alanine did not produce a remarkable change in the total activity or the substrate enantioselectivity. These results indicate that Glu-300 and Phe-476 are important in stereoselective oxidation of HB enantiomers by CYP2C19.

Selective, competitive and mechanism-based inhibitors of human cytochrome P450 2J2

Twenty five derivatives of the drugs terfenadine and ebastine have been designed, synthesized and evaluated as inhibitors of recombinant human CYP2J2. Compound 14, which has an imidazole substituent, is a good non-competitive inhibitor of CYP2J2 (IC(50)=400nM). It is not selective towards CYP2J2 as it also efficiently inhibits the other main vascular CYPs, such as CYP2B6, 2C8, 2C9 and 3A4; however, it could be an interesting tool to inhibit all these vascular CYPs. Compounds 4, 5 and 13, which have a propyl, allyl and benzo-1,3-dioxole terminal group, respectively, are selective CYP2J2 inhibitors. Compound 4 is a high-affinity, competitive inhibitor and alternative substrate of CYP2J2 (K(i)=160+/-50nM). Compounds 5 and 13 are efficient mechanism-based inhibitors of CYP2J2 (k(inact)/K(i) values approximately 3000Lmol(-1)s(-1)). Inactivation of CYP2J2 by 13 is due to the formation of a stable iron-carbene bond which occurs upon CYP2J2-catalyzed oxidation of 13 with a partition ratio of 18+/-3. These new selective inhibitors should be interesting tools to study the biological roles of CYP2J2.