Diversity in mechanisms of substrate oxidation by cytochrome P450 2D6. Lack of an allosteric role of NADPH-cytochrome P450 reductase in catalytic regioselectivity

Roles of selected non-P450 human oxidoreductase enzymes in protective and toxic effects of chemicals: review and compilation of reactions

This is an overview of the metabolic reactions of drugs, natural products, physiological compounds, and other (general) chemicals catalyzed by flavin monooxygenase (FMO), monoamine oxidase (MAO), NAD(P)H quinone oxidoreductase (NQO), and molybdenum hydroxylase enzymes (aldehyde oxidase (AOX) and xanthine oxidoreductase (XOR)), including roles as substrates, inducers, and inhibitors of the enzymes. The metabolism and bioactivation of selected examples of each group (i.e., drugs, "general chemicals," natural products, and physiological compounds) are discussed. We identified a higher fraction of bioactivation reactions for FMO enzymes compared to other enzymes, predominately involving drugs and general chemicals. With MAO enzymes, physiological compounds predominate as substrates, and some products lead to unwanted side effects or illness. AOX and XOR enzymes are molybdenum hydroxylases that catalyze the oxidation of various heteroaromatic rings and aldehydes and the reduction of a number of different functional groups. While neither of these two enzymes contributes substantially to the metabolism of currently marketed drugs, AOX has become a frequently encountered route of metabolism among drug discovery programs in the past 10-15 years. XOR has even less of a role in the metabolism of clinical drugs and preclinical drug candidates than AOX, likely due to narrower substrate specificity.

The regioselectivity of the interaction between dextromethorphan and CYP2D6

CYP2D6 is an important enzyme of the cytochrome P450 superfamily, and catalyzes nearly 25% of the drugs sold in the market. For decades, the interactions and metabolism between CYP2D6 and substrates have been a hot topic. However, the key factors of the catalytic regioselectivity for CYP2D6 still remain controversial. Here, we construct four systems to explore the interaction between dextromethorphan (DM) and CYP2D6. A new binding mode of CYP2D6 is defined, and two key residues (residue Asp301 and residue Glu216) are discovered working simultaneously to stabilize the DM at the reactive site by forming water bridge hydrogen bonds when CYP2D6 binds DM. Our results also indicate that the substrate concentration could mediate the binding mode between the substrate and CYP2D6 by decreasing the volume of the catalytic pocket, which is not conducive to the O-demethylation of DM but benefits the N-demethylation of DM. These results could shed light on the process of CYP2D6 binding to the substrate, and help to better understand the regioselectivity of CYP2D6 catalyzing the substrates.

Oxygen Surrogate Systems for Supporting Human Drug-Metabolizing Cytochrome P450 Enzymes

Oxygen surrogates (OSs) have been used to support cytochrome P450 (P450) enzymes for diverse purposes in drug metabolism research, including reaction phenotyping, mechanistic and inhibition studies, studies of redox partner interactions, and to avoid the need for NADPH or a redox partner. They also have been used in engineering P450s for more cost-effective, NADPH-independent biocatalysis. However, despite their broad application, little is known of the preference of individual P450s for different OSs or the substrate dependence of OS-supported activity. Furthermore, the biocatalytic potential of OSs other than cumene hydroperoxide (CuOOH) and hydrogen peroxide (H₂O₂) is yet to be explored. Here, we investigated the ability of the major human drug-metabolizing P450s, namely CYP3A4, CYP2C9, CYP2C19, CYP2D6, and CYP1A2, to use the following OSs: H₂O₂, tert-butyl hydroperoxide (tert-BuOOH), CuOOH, (diacetoxyiodo)benzene, and bis(trifluoroacetoxy)iodobenzene. Overall, CuOOH and tert-BuOOH were found to be the most effective at supporting these P450s. However, the ability of P450s to be supported by OSs effectively was also found to be highly dependent on the substrate used. This suggests that the choice of OS should be tailored to both the P450 and the substrate under investigation, underscoring the need to employ screening methods that reflect the activity toward the substrate of interest to the end application. SIGNIFICANCE STATEMENT: Cytochrome P450 (P450) enzymes can be supported by different oxygen surrogates (OSs), avoiding the need for a redox partner and costly NADPH. However, few data exist comparing relative activity with different OSs and substrates. This study shows that the choice of OS used to support the major drug-metabolizing P450s influences their relative activity and regioselectivity in a substrate-specific fashion and provides a model for the more efficient use of P450s for metabolite biosynthesis.

The relationships between cytochromes P450 and H2O2: Production, reaction, and inhibition

In this review we address the relationship between cytochromes P450 (P450) and H2O2. This association can affect biology in three distinct ways. First, P450s produce H2O2 as a byproduct either during catalysis or when no substrate is present. This reaction, known as uncoupling, releases reactive oxygen species that may have implications in disease. Second, H2O2 is used as an oxygen-donating co-substrate in peroxygenase and peroxidase reactions catalyzed by P450s. This activity has proven to be important mainly in reactions involving prokaryotic P450s, and investigators have harnessed this reaction with the aim of adaptation for industrial use. Third, H2O2-dependent inhibition of human P450s has been studied in our laboratory, demonstrating heme destruction and also the inactivating oxidation of the heme-thiolate ligand to a sulfenic acid (-SOH). This reversible oxidative modification of P450s may have implications in the prevention of uncoupling and may give new insights into the oxidative regulation of these enzymes. Research has elucidated many of the chemical mechanisms involved in the relationship between P450 and H2O2, but the application to biology is difficult to evaluate. Further studies are needed reveal both the harmful and protective natures of reactive oxygen species in an organismal context.

Structure and function of the cytochrome P450 peroxygenase enzymes

The cytochromes P450 (P450s or CYPs) constitute a large heme enzyme superfamily, members of which catalyze the oxidative transformation of a wide range of organic substrates, and whose functions are crucial to xenobiotic metabolism and steroid transformation in humans and other organisms. The P450 peroxygenases are a subgroup of the P450s that have evolved in microbes to catalyze the oxidative metabolism of fatty acids, using hydrogen peroxide as an oxidant rather than NAD(P)H-driven redox partner systems typical of the vast majority of other characterized P450 enzymes. Early members of the peroxygenase (CYP152) family were shown to catalyze hydroxylation at the α and β carbons of medium-to-long-chain fatty acids. However, more recent studies on other CYP152 family P450s revealed the ability to oxidatively decarboxylate fatty acids, generating terminal alkenes with potential applications as drop-in biofuels. Other research has revealed their capacity to decarboxylate and to desaturate hydroxylated fatty acids to form novel products. Structural data have revealed a common active site motif for the binding of the substrate carboxylate group in the peroxygenases, and mechanistic and transient kinetic analyses have demonstrated the formation of reactive iron-oxo species (compounds I and II) that are ultimately responsible for hydroxylation and decarboxylation of fatty acids, respectively. This short review will focus on the biochemical properties of the P450 peroxygenases and on their biotechnological applications with respect to production of volatile alkenes as biofuels, as well as other fine chemicals.

Detoxification of 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP) by cytochrome P450 enzymes: A theoretical investigation

Two types of detoxification routes, N-demethylation to form 4-phenyl-1,2,3,6-tetrahydropyridine (PTP) and aromatic hydroxylation to generate 4-(4'-hydroxyphenyl)-1-methyl-1,2,3,6-tetrahydropyridine (MPTP-OH), for 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP) mediated by Compound I (Cpd I) of cytochrome P450 are investigated theoretically using hybrid density functional calculations. Quantum chemical results reveal that for the N-demethylation, the initial C-H bond activation is achieved via a hydrogen atom transfer (HAT) mechanism. This is followed by a subsequent O-rebound to yield the carbinolamine intermediate. Due to the nature of pericyclic reaction, the generated carbinolamine decomposes in a non-enzymatic aqueous environment with the assistance of water molecules, forming amine and hydrated formaldehyde. For the aromatic hydroxylation, an initial addition of Cpd I to the substrate occurs mainly through a side-on approach with a subsequent proton shuttle to form the phenol product. A comparison of the energy barriers for both routes indicates that the N-demethylation (7.5/5.7kcal/mol for the quartet/doublet state in solvent) is thermodynamically more favorable than the aromatic hydroxylation process (14.9/14.8kcal/mol for the quartet/doublet state in solvent). This trend is in good agreement with the experimental product distribution, viz., the N-demethylation product PTP is more than the aromatic hydroxylation product MPTP-OH. Taken together, these observations not only enrich our knowledge on the mechanistic details of the N-dealkylation and the aromatic hydroxylation by P450s, but also provide certain insights into the metabolism of other analogous toxins.

Integrated structure- and ligand-based in silico approach to predict inhibition of cytochrome P450 2D6

MOTIVATION

Cytochrome P450 (CYP) is a superfamily of enzymes responsible for the metabolism of drugs, xenobiotics and endogenous compounds. CYP2D6 metabolizes about 30% of drugs and predicting potential CYP2D6 inhibition is important in early-stage drug discovery.

RESULTS

We developed an original in silico approach for the prediction of CYP2D6 inhibition combining the knowledge of the protein structure and its dynamic behavior in response to the binding of various ligands and machine learning modeling. This approach includes structural information for CYP2D6 based on the available crystal structures and molecular dynamic simulations (MD) that we performed to take into account conformational changes of the binding site. We performed modeling using three learning algorithms--support vector machine, RandomForest and NaiveBayesian--and we constructed combined models based on topological information of known CYP2D6 inhibitors and predicted binding energies computed by docking on both X-ray and MD protein conformations. In addition, we identified three MD-derived structures that are capable all together to better discriminate inhibitors and non-inhibitors compared with individual CYP2D6 conformations, thus ensuring complementary ligand profiles. Inhibition models based on classical molecular descriptors and predicted binding energies were able to predict CYP2D6 inhibition with an accuracy of 78% on the training set and 75% on the external validation set.

Toward reduction in animal sacrifice for drugs: molecular modeling of Macaca fascicularis P450 2C20 for virtual screening of Homo sapiens P450 2C8 substrates

Macaca fascicularis P450 2C20 shares 92% identity with human cytochrome P450 2C8, which is involved in the metabolism of more than 8% of all prescribed drugs. To date, only paclitaxel and amodiaquine, two substrate markers of the human P450 2C8, have been experimentally confirmed as M. fascicularis P450 2C20 drugs. To bridge the lack of information on the ligands recognized by M. fascicularis P450 2C20, in this study, a three-dimensional homology model of this enzyme was generated on the basis of the available crystal structure of the human homologue P450 2C8 using YASARA. The results indicated that 90.0%, 9.0%, 0.5%, and 0.5% of the residues of the P450 2C20 model were located in the most favorable, allowed, generously allowed, and disallowed regions, respectively. The root-mean-square deviation of the C-alpha superposition of the M. fascicularis P450 2C20 model with the Homo sapiens P450 2C8 was 0.074 Å, indicating a very high similarity of the two structures. Subsequently, the 2C20 model was used for in silico screening of 58 known P450 2C8 substrates and 62 inhibitors. These were also docked in the active site of the crystal structure of the human P450 2C8. The affinity of each compound for the active site of both cytochromes proved to be very similar, meaning that the few key residues that are mutated in the active site of the M. fascicularis P450 do not prevent the P450 2C20 from recognizing the same substrates as the human P450 2C8.

Metabolite profile resulting from the activation/inactivation of 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine and 2-methyltetrahydro-β-carboline by oxidative enzymes

Metabolic enzymes are involved in the activation/deactivation of the 1-methyl-4-phenyl-1,2,3,6-tetrahydropyiridine (MPTP) neurotoxin and its naturally occurring analogs 2-methyltetrahydro-β-carbolines. The metabolic profile and biotransformation of these protoxins by three enzymes, monoamine oxidase (MAO), cytochrome P450, and heme peroxidases (myeloperoxidase and lactoperoxidase), were investigated and compared. The metabolite profile differed among the enzymes investigated. MAO and heme peroxidases activated these substances to toxic pyridinium and β-carbolinium species. MAO catalyzed the oxidation of MPTP to 1-methyl-4-phenyl-2,3-dihydropyridinium cation (MPDP(+)), whereas heme peroxidases catalyzed the oxidation of MPDP(+) to 1-methyl-4-phenylpyridinium (MPP(+)) and of 2-methyltetrahydro-β-carboline to 2-methyl-3,4-dihydro-β-carbolinium cation (2-Me-3,4-DH β C(+)). These substances were inactivated by cytochrome P450 2D6 through N-demethylation and aromatic hydroxylation (MPTP) and aromatic hydroxylation (2-methyltetrahydro-β-carboline). In conclusion, the toxicological effects of these protoxins might result from a balance between the rate of their activation to toxic products (i.e., N-methylpyridinium-MPP(+) and MPDP(+)- and N-methyl--β--carbolinium- βC(+)-) by MAO and heme peroxidases and the rate of inactivation (i.e., N-demethylation, aromatic hydroxylation) by cytochrome P450 2D6.

Substituted imidazole of 5-fluoro-2-[4-[(2-phenyl-1H-imidazol-5-yl)methyl]-1-piperazinyl]pyrimidine Inactivates cytochrome P450 2D6 by protein adduction

5-Fluoro-2-[4-[(2-phenyl-1H-imidazol-5-yl)methyl]-1-piperazinyl]pyrimidine (SCH 66712) is a potent mechanism-based inactivator of human cytochrome P450 2D6 that displays type I binding spectra with a K(s) of 0.39 ± 0.10 μM. The partition ratio is ~3, indicating potent inactivation that addition of exogenous nucleophiles does not prevent. Within 15 min of incubation with SCH 66712 and NADPH, ∼90% of CYP2D6 activity is lost with only ~20% loss in ability to bind CO and ~25% loss of native heme over the same time. The stoichiometry of binding to the protein was 1.2:1. SDS-polyacrylamide gel electrophoresis with Western blotting and autoradiography analyses of CYP2D6 after incubations with radiolabeled SCH 66712 further support the presence of a protein adduct. Metabolites of SCH 66712 detected by mass spectrometry indicate that the phenyl group on the imidazole ring of SCH 66712 is one site of oxidation by CYP2D6 and could lead to methylene quinone formation. Three other metabolites were also observed. For understanding the metabolic pathway that leads to CYP2D6 inactivation, metabolism studies with CYP2C9 and CYP2C19 were performed because neither of these enzymes is significantly inhibited by SCH 66712. The metabolites formed by CYP2C9 and CYP2C19 are the same as those seen with CYP2D6, although in different abundance. Modeling studies with CYP2D6 revealed potential roles of various active site residues in the oxidation of SCH 66712 and inactivation of CYP2D6 and showed that the phenyl group of SCH 66712 is positioned at 2.2 Å from the heme iron.

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.

Polymorphism of human cytochrome P450 enzymes and its clinical impact

Pharmacogenetics is the study of how interindividual variations in the DNA sequence of specific genes affect drug response. This article highlights current pharmacogenetic knowledge on important human drug-metabolizing cytochrome P450s (CYPs) to understand the large interindividual variability in drug clearance and responses in clinical practice. The human CYP superfamily contains 57 functional genes and 58 pseudogenes, with members of the 1, 2, and 3 families playing an important role in the metabolism of therapeutic drugs, other xenobiotics, and some endogenous compounds. Polymorphisms in the CYP family may have had the most impact on the fate of therapeutic drugs. CYP2D6, 2C19, and 2C9 polymorphisms account for the most frequent variations in phase I metabolism of drugs, since almost 80% of drugs in use today are metabolized by these enzymes. Approximately 5-14% of Caucasians, 0-5% Africans, and 0-1% of Asians lack CYP2D6 activity, and these individuals are known as poor metabolizers. CYP2C9 is another clinically significant enzyme that demonstrates multiple genetic variants with a potentially functional impact on the efficacy and adverse effects of drugs that are mainly eliminated by this enzyme. Studies into the CYP2C9 polymorphism have highlighted the importance of the CYP2C9*2 and *3 alleles. Extensive polymorphism also occurs in other CYP genes, such as CYP1A1, 2A6, 2A13, 2C8, 3A4, and 3A5. Since several of these CYPs (e.g., CYP1A1 and 1A2) play a role in the bioactivation of many procarcinogens, polymorphisms of these enzymes may contribute to the variable susceptibility to carcinogenesis. The distribution of the common variant alleles of CYP genes varies among different ethnic populations. Pharmacogenetics has the potential to achieve optimal quality use of medicines, and to improve the efficacy and safety of both prospective and currently available drugs. Further studies are warranted to explore the gene-dose, gene-concentration, and gene-response relationships for these important drug-metabolizing CYPs.

Mammalian cytochrome P450 enzymes catalyze the phenol-coupling step in endogenous morphine biosynthesis

A cytochrome P450 (P450) enzyme in porcine liver that catalyzed the phenol-coupling reaction of the substrate (R)-reticuline to salutaridine was previously purified to homogeneity (Amann, T., Roos, P. H., Huh, H., and Zenk, M. H. (1995) Heterocycles 40, 425-440). This reaction was found to be catalyzed by human P450s 2D6 and 3A4 in the presence of (R)-reticuline and NADPH to yield not a single product, but rather (-)-isoboldine, (-)-corytuberine, (+)-pallidine, and salutaridine, the para-ortho coupled established precursor of morphine in the poppy plant and most likely also in mammals. (S)-Reticuline, a substrate of both P450 enzymes, yielded the phenol-coupled alkaloids (+)-isoboldine, (+)-corytuberine, (-)-pallidine, and sinoacutine; none of these serve as a morphine precursor. Catalytic efficiencies were similar for P450 2D6 and P450 3A4 in the presence of cytochrome b(5) with (R)-reticuline as substrate. The mechanism of phenol coupling is not yet established; however, we favor a single cycle of iron oxidation to yield salutaridine and the three other alkaloids from (R)-reticuline. The total yield of salutaridine formed can supply the 10 nm concentration of morphine found in human neuroblastoma cell cultures and in brain tissues of mice.

Insights into drug metabolism by cytochromes P450 from modelling studies of CYP2D6-drug interactions

The cytochromes P450 (CYPs) comprise a vast superfamily of enzymes found in virtually all life forms. In mammals, xenobiotic metabolizing CYPs provide crucial protection from the effects of exposure to a wide variety of chemicals, including environmental toxins and therapeutic drugs. Ideally, the information on the possible metabolism by CYPs required during drug development would be obtained from crystal structures of all the CYPs of interest. For some years only crystal structures of distantly related bacterial CYPs were available and homology modelling techniques were used to bridge the gap and produce structural models of human CYPs, and thereby obtain useful functional information. A significant step forward in the reliability of these models came seven years ago with the first crystal structure of a mammalian CYP, rabbit CYP2C5, followed by the structures of six human enzymes, CYP1A2, CYP2A6, CYP2C8, CYP2C9, CYP2D6 and CYP3A4, and a second rabbit enzyme, CYP2B4. In this review we describe as a case study the evolution of a CYP2D6 model, leading to the validation of the model as an in silico tool for predicting binding and metabolism. This work has led directly to the successful design of CYP2D6 mutants with novel activity-including creating a testosterone hydroxylase, converting quinidine from inhibitor to substrate, creating a diclofenac hydroxylase and creating a dextromethorphan O-demethylase. Our modelling-derived hypothesis-driven integrated interdisciplinary studies have given key insight into the molecular determinants of CYP2D6 and other important drug metabolizing enzymes.

Binding of bufuralol, dextromethorphan, and 3,4-methylenedioxymethylamphetamine to wild-type and F120A mutant cytochrome P450 2D6 studied by resonance Raman spectroscopy

Cytochrome P450 2D6 (CYP2D6) is one of the most important drug-metabolizing enzymes in humans. Resonance Raman data, reported for the first time for CYP2D6, show that the CYP2D6 heme is found to be in a six-coordinated low-spin state in the absence of substrates, and it is perturbed to different extents by bufuralol, dextromethorphan, and 3,4-methylenedioxymethylamphetamine (MDMA). Dextromethorphan and MDMA induce in CYP2D6 a significant amount of five-coordinated high-spin heme species and reduce the polarity of its heme-pocket, whereas bufuralol does not. Spectra of the F120A mutant CYP2D6 suggest that Phe120 is involved in substrate-binding of dextromethorphan and MDMA, being responsible for the spectral differences observed between these two compounds and bufuralol. These differences could be explained postulating a different substrate mobility for each compound in the CYP2D6 active site, consistently with the role previously suggested for Phe120 in binding dextromethorphan and MDMA.

Topological role of cytochrome P450 2D6 active site residues

Recent reports have identified Phe120, Asp301, Thr309, and Glu216 as important residues in cytochrome P450 2D6 (CYP2D6) substrate binding and catalysis. Complementary homology models have located these amino acids within the binding pocket of CYP2D6 and in the present study we have used aryldiazenes to test these models and gain further insight in the role these amino acids have in maintaining the integrity of the active site cavity. When Phe120 was replaced to alanine, there was a significant increase in probe migration to pyrrole nitrogens C and D, in agreement with homology models which have located the phenyl side-chain of Phe120 above these two pyrrole rings. No changes in topology were observed with the D301Q mutant, supporting claims that in this mutant the electrostatic interactions with the B/C-loop are largely maintained and the loop retains its native orientation. The T309V mutation resulted in significant topological alteration suggesting that, in addition to its potential role in dioxygen activation, Thr309 plays an important structural role within the active site crevice. Replacement of Ile106 with Glu, engineered to cause electrostatic repulsion with Glu216, had a profound topological effect in the higher region within the active site cavity and impaired the catalytic activity towards CYP2D6 probe substrates.

Role of the conserved threonine 309 in mechanism of oxidation by cytochrome P450 2D6

Based on sequence alignments and homology modeling, threonine 309 in cytochrome P450 2D6 (CYP2D6) is proposed to be the conserved I-helix threonine, which is supposed to be involved in dioxygen activation by CYPs. The T309V mutant of CYP2D6 displayed a strong shift from O-dealkylation to N-dealkylation reactions in oxidation of dextromethorphan and 3,4-methylenedioxymethylamphetamine. This may be explained by an elevated ratio of hydroperoxo-iron to oxenoid-iron of the oxygenating species. In consistence, using cumene hydroperoxide, which directly forms the oxenoid-iron, the T309V mutant again selectively catalyzed the O-dealkylation reactions. The changed ratio of oxygenating species can also explain the decreased activity and changed regioselectivity that were observed in 7-methoxy-4-(aminomethyl)-coumarin and bufuralol oxidation, respectively, by the T309V mutant. Interestingly, the T309V mutant always showed a significantly increased, up to 75-fold, higher activity compared to that of the wild-type when using cumene hydroperoxide. These results indicate that T309 in CYP2D6 is involved in maintaining the balance of multiple oxygenating species and thus influences substrate and regioselectivity.

The role of phenylalanine 483 in cytochrome P450 2D6 is strongly substrate dependent

The polymorphic cytochrome P450 2D6 (CYP2D6) is involved in the metabolism of 30% of the drugs currently prescribed, and is thus clinically relevant. Typical CYP2D6 substrates generally contain a basic nitrogen atom and an aromatic moiety adjacent to the site of metabolism. Recently, we demonstrated the importance of active site residue F120 in substrate binding and catalysis in CYP2D6. On the basis of protein homology models, it is claimed that another active site phenylalanine, F483, may also play an important role in the interaction with the aromatic moiety of CYP2D6 substrates. Experimental data to support this hypothesis, however, is not yet available. In fact, in the only study performed, mutation of F483 to isoleucine or tryptophan did not affect the 1'-hydroxylation of bufuralol at all [Smith G, Modi S, Pillai I, Lian LY, Sutcliffe MJ, Pritchard MP, et al., Determinants of the substrate specificity of human cytochrome P-450 CYP2D6: design and construction of a mutant with testosterone hydroxylase activity. Biochem J 1998;331:783-92]. In the present study, the role of F483 in ligand binding and metabolism by CYP2D6 was examined experimentally using site-directed mutagenesis. Replacement of F483 by alanine resulted in a 30-fold lower V(max) for bufuralol 1'-hydroxylation, while the K(m) was hardly affected. The V(max) for 3,4-methylenedioxy-methylamphetamine O-demethylenation on the other hand decreased only two-fold, whereas the effect on the K(m) was much larger. For dextromethorphan, in addition to dextrorphan (O-demethylation) and 3-methoxymorphinan (N-demethylation), two other metabolites were formed that could not be detected for the wild-type. The substrate 7-methoxy-4-(aminomethyl)-coumarin was not metabolised at all by CYP2D6[F483A], a phenomenon that was reported also for CYP2D6[F120A]. The presented data show that next to F120, residue F483 plays a very important role in the metabolism of typical CYP2D6 substrates. The influence of F483 on metabolism was found to be strongly substrate-dependent.

NON-MICHAELIS-MENTEN KINETICS IN CYTOCHROME P450-CATALYZED REACTIONS

The cytochrome P450 monooxygenases (CYPs) are the dominant enzyme system responsible for xenobiotic detoxification and drug metabolism. Several CYP isoforms exhibit non-Michaelis-Menten, or “atypical,” steady state kinetic patterns. The allosteric kinetics confound prediction of drug metabolism and drug-drug interactions, and they challenge the theoretical paradigms of allosterism. Both homotropic and heterotropic ligand effects are now widely documented. It is becoming apparent that multiple ligands can simultaneously bind within the active sites of individual CYPs, and the kinetic parameters change with ligand occupancy. In fact, the functional effect of any specific ligand as an activator or inhibitor can be substrate dependent. Divergent approaches, including kinetic modeling and X-ray crystallography, are providing new information about how multiple ligand binding yields complex CYP kinetics.

Metabolic Regio- and Stereoselectivity of Cytochrome P450 2D6 towards 3,4-Methylenedioxy-N-alkylamphetamines: in Silico Predictions and Experimental Validation

A series of 3,4-methylenedioxy-N-alkylamphetamines (MDAAs) were automatically docked and subjected to molecular dynamics (MD) simulations in a cytochrome P450 2D6 (CYP2D6) protein model. The predicted substrate binding orientations, sites of oxidation, and relative reactivities were compared to the experimental data of wild-type and Phe120Ala mutant CYP2D6. Automated docking results were not sufficient to accurately rationalize experimental binding orientations of 3,4-methylenedioxy-N-methylamphetamine (MDMA) in the two enzymes as measured with spin lattice relaxation NMR. Nevertheless, the docking results could be used as starting structures for MD simulations. Predicted binding orientations of MDMA and sites of oxidation of the MDAAs derived from MD simulations matched well with the experimental data. It appeared the experimental results were best described in MD simulations considering the nitrogen atoms of the MDAAs in neutral form. Differences in regioselectivity and stereoselectivity in the oxidative metabolism of the MDAAs by the Phe120Ala mutant CYP2D6 were correctly predicted, and the effects of the Phe120Ala mutation could be rationalized as well.

P(450)/NADPH/O(2)- and P(450)/PhIO-catalyzed N-dealkylations are mechanistically distinct

A high-valent iron-oxo species analogous to the compound I of peroxidases has been thought to be the activated oxygen species in P450-catalyzed reactions. Spectroscopic characterization of the catalytically competent iron-oxo species in iodosobenzene (PhIO)-supported model reactions and parallels between these model reactions and PhIO- and NADPH/O2-supported P450 reactions have been taken as strong evidence for this proposal. To support this proposal, subtle differences observed in regio- and chemoselectivities, isotope effects, and source of oxygen, etc., between NADPH/O2- and PhIO-supported P450 reactions have been generally attributed to reasons other than the mechanistic differences between the two systems. In the present study, we have used a series of sensitive mechanistic probes, 4-chloro-N-cyclopropyl-N-alkylanilines, to compare and contrast the chemistries of the NADPH/O2- and PhIO-supported purified CYP2B1 N-dealkylation reactions. Herein we present the first experimental evidence to demonstrate that the NADPH/O2- and PhIO-supported P450 N-dealkylations are mechanistically distinct and, thus, the P450/PhIO system may not be a good mechanistic model for P450/NADPH/O2-catalyzed N-dealkylations.

Progress in cytochrome P450 active site modeling

Models capable of predicting the possible involvement of cytochromes P450 in the metabolism of drugs or drug candidates are important tools in drug discovery and development. Ideally, functional information would be obtained from crystal structures of all the cytochromes P450 of interest. Initially, only crystal structures of distantly related bacterial cytochromes P450 were available-comparative modeling techniques were used to bridge the gap and produce structural models of human cytochromes P450, and thereby obtain some useful functional information. A significant step forward in the reliability of these models came four years ago with the first crystal structure of a mammalian cytochrome P450, rabbit CYP2C5, followed by the structures of two human enzymes, CYP2C8 and CYP2C9, and a second rabbit enzyme, CYP2B4. The evolution of a CYP2D6 model, leading to the validation of the model as an in silico tool for predicting binding and metabolism, is presented as a case study.

Influence of phenylalanine 120 on cytochrome P450 2D6 catalytic selectivity and regiospecificity: crucial role in 7-methoxy-4-(aminomethyl)-coumarin metabolism

The polymorphic human debrisoquine hydroxylase, cytochrome P450 2D6 (CYP2D6), is one of the most important phase I drug metabolising enzymes. It is responsible for metabolising a large number of compounds that mostly share similarity in having a basic N-atom and an aromatic moiety. In homology modelling studies, it has been suggested that in fixation of this aromatic moiety, there may be an important role for phenylalanine 120 (Phe(120)). In this study, the role of Phe(120) in ligand binding and catalysis was experimentally examined by mutating it into an alanine. Strikingly, this substitution led to a completely abolished 7-methoxy-4-(aminomethyl)-coumarin (MAMC) O-demethylating activity of CYP2D6. On the other hand, bufuralol metabolism was hardly affected (K(m) of 1-hydroxylation mutant: 1.2 microM, wild-type: 2.9 microM, 4-hydroxylation mutant: 1.5 microM, and wild-type: 3.2 microM) and neither was affected dextromethorphan O-demethylation (K(m) mutant: 1.2 microM, wild-type: 2 microM, k(cat) mutant: 4.5 min(-1), and wild-type: 3.3 min(-1)). However, the Phe(120)Ala mutant also formed 3-hydroxymorphinan, the double demethylated form of dextromethorphan, which was not detected using wild-type CYP2D6. 3,4-Methylenedioxymethamphetamine (MDMA) was demethylenated by both mutant and wild-type CYP2D6 to 3,4-dihydroxymethamphetamine (3,4-OH-MA K(m) of mutant: 55 microM and wild-type: 2 microM). In addition, the mutant formed two additional metabolites; 3,4-methylenedioxyamphetamine (MDA) and N-hydroxy-3,4-methylenedioxymethamphetamine (N-OH-MDMA). Inhibition experiments of dextromethorphan O-demethylation showed a decreased affinity of the Phe(120)Ala mutant for quinidine (IC(50) mutant: 240 nM and wild-type, 40 nM), while IC(50)s for quinine were equal (1 microM). These data indicate the importance of Phe(120) in the selectivity and regiospecificity in substrate binding and catalysis by CYP2D6.

Formation of a novel reversible cytochrome P450 spectral intermediate: role of threonine 303 in P450 2E1 inactivation

tert-Butyl acetylene (tBA) is a mechanism-based inactivator of cytochromes P450 2E1 and 2E1 T303A; however, the inactivation of the T303A mutant could be reversed by overnight dialysis. The inactivation of P450 2E1 T303A, but not the wild-type 2E1 enzyme, by tBA resulted in the formation of a novel reversible acetylene-iron spectral intermediate with an absorption maximum at 485 nm. The formation of this intermediate required oxygen and could be monitored spectrally with time. Although the alternate oxidants tert-butyl hydroperoxide (tBHP) and cumene hydroperoxide (CHP) supported the inactivation of wild-type P450 2E1 by tBA in a reductase- and NADPH-free system, only tBHP supported the inactivation of the 2E1 T303A mutant. The losses in enzymatic activity occurred concomitantly with losses in the native P450 heme, which were accompanied by the formation of tBA-adducted heme products. The inactivations supported by tBHP and CHP were completely irreversible with overnight dialysis. Spectral binding constants (K(s)) for the binding of tBA to the 2E1 P450s together with models of the enzymes with the acetylenic inactivator bound in the active site suggest that the T303A mutation results in increased hydrophobic interactions between tBA and nearby P450 residues, leading to a higher binding affinity for the acetylene compound in the mutant enzyme. Together, these data support a role for the highly conserved T303 residue in proton delivery to the active site of P450 2E1 and in the inactivation of the 2E1 P450s by small acetylenic compounds.

Rate-limiting steps in oxidations catalyzed by rabbit cytochrome P450 1A2

Several issues regarding the rate-limiting nature of individual reaction steps in catalysis by rabbit liver cytochrome P450 (P450) 1A2 were addressed using anisoles and other substrates. Substrate binding is very fast (k > 10(6) M(-1) s(-1)). Product release is not rate-limiting, as shown by the absence of bursts, placing rate-limiting steps at or before product formation. We had previously shown that the first 1-electron reduction step is fast (k > 700 min(-1)), even in the absence of ligand [Guengerich, F. P., and Johnson, W. W. (1997) Biochemistry 36, 14741-147500]. O(2) binding to ferrous P450 is fast (k >/= 10(6) M(-1) s(-1)). The decay of the P450 Fe(2+)-substrate-O(2) complex was slow in the absence of NADPH-P450 reductase, with a first-order rate constant of 14 min(-1) at 25 degrees C. During the decay, product was formed (from the substrate methacetin) in 61% theoretical yield, although this reaction requires electron transfer among P450 molecules and may not be related to normal turnover. Steady-state spectra suggest that one or more iron-oxygen complexes accumulate, representing entities between the Fe(2+)-O(2) complex and putative FeO(3+) entity. Kinetic isotope effect experiments were done with several substrates, mainly anisoles. Apparent intrinsic deuterium isotope effects as high as 15 were measured. In all cases, the C-H bond-breaking step is at least partially rate-limiting. The isotope effects were not strongly attenuated in noncompetitive or competitive experiments, consistent with relatively rapid P450-substrate exchange, except with the active enzyme Fe-O complex. Kinetic simulations with the available data (i) are consistent with the view that C-H bond breaking is a major rate-limiting step, (ii) demonstrate that increasing the rate of this step will affect k(cat), K(m), and kinetic hydrogen isotope effects but will only increase catalytic efficiency to a certain degree, (iii) indicate that increasing ground-state binding can increase catalytic efficiency but not k(cat), and (iv) suggest that nonproductive binding modes and abortive reduction of O(2) are factors that attenuate catalytic efficiency.

Binding of low molecular mass compounds to proteins studied by liquid chromatographic techniques

The newest achievements in the application of miscellaneous liquid chromatographic techniques such as size-exclusion, ion-exchange and reversed-phase high-performance liquid chromatography, and thin-layer chromatography for the elucidation of the various aspects of the binding of ligands to proteins are compiled and briefly discussed. Examples of employment in pharmaceutical and clinical chemistry, drug design, enzyme kinetic studies and environmental protection are presented.

Role of glutamic acid 216 in cytochrome P450 2D6 substrate binding and catalysis

Human cytochrome P450 (P450) 2D6 is an important enzyme involved in the metabolism of drugs, many of which are amines or contain other basic nitrogen atoms. Asp301 has generally been considered to be involved in electrostatic docking with the basic substrates, on the basis of previous modeling studies and site-directed mutagenesis. Substitution of Glu216 with a residue other than Asp strongly attenuated the binding of quinidine, bufuralol, and several other P450 2D6 ligands. Catalytic activity with the substrates bufuralol and 4-methoxyphenethylamine was strongly inhibited by neutral or basic mutations at Glu216 (>95%), to the same extent as the substitution of Asn at Asp301. Unlike the Asp301 mutants, the Gln216 mutant (E216Q) retained 40% enzyme efficiency with the substrate spirosulfonamide, devoid of basic nitrogen, suggesting that the substitutions at Glu216 affect binding of amine substrates more than other catalytic steps. Attempts to induce catalytic specificity toward new substrates by substitutions at Asp301 and Glu216 were unsuccessful. Collectively, the results provide evidence for electrostatic interaction of amine substrates with Glu216, and we propose that both of these acidic residues plus at least another residue(s) is (are) involved in binding the repertoire of P450 2D6 ligands.

Oxidation of methoxyphenethylamines by cytochrome P450 2D6. Analysis of rate-limiting steps

Cytochrome P450 (P450) 2D6 is involved in the oxidation of a large fraction ( approximately 30%) of drugs used by humans and also catalyzes the O-demethylation of the model substrates 3- and 4-methoxyphenethylamine followed by subsequent ring hydroxylation to dopamine. Burst kinetics were not observed; rate-limiting step(s) must occur prior to product formation. Rates of reduction of ferric P450 2D6 were stimulated by 3- or 4-methoxyphenethylamine or the inhibitor quinidine; reduction is not the most rate-limiting step. The non-competitive intramolecular deuterium isotope effect, an estimate of the intrinsic isotope effect, for 4-methoxyphenethylamine O-demethylation was 9.6. Intermolecular non-competitive deuterium isotope effects of 3.1-3.8 were measured for k(cat) and k(cat)/K(m) for both O-demethylation reactions, implicating at least partially rate-limiting C-H bond breaking. Simulation of steady-state kinetic data yielded a catalytic mechanism dominated by the rates of (i) Fe(2+)O(2)(-) protonation (plus O-O bond scission) and (ii) C-H bond breaking, consistent with the appearance of the spectral intermediates in the steady state, attributed to iron-oxygen complexes. However, all the rates of individual steps (or rates of combined steps) are considerably higher than k(cat), and the contributions of several steps must be considered in understanding rates of the P450 2D6 reactions.

Diversity in the oxidation of substrates by cytochrome P450 2D6: lack of an obligatory role of aspartate 301-substrate electrostatic bonding

Cytochrome P450 (P450) 2D6 was first identified as the polymorphic human debrisoquine hydroxylase and subsequently shown to catalyze the oxidation of a variety of drugs containing a basic nitrogen. Residue Asp301 has been characterized as being involved in electrostatic interactions with substrates on the basis of homology modeling and site-directed mutagenesis experiments [Ellis, S. W., Hayhurst, G. P., Smith, G., Lightfoot, T., Wong, M. M. S., Simula, A. P., Ackland, M. J., Sternberg, M. J. E., Lennard, M. S., Tucker, G. T., and Wolf, C. R. (1995) J. Biol. Chem. 270, 29055-29058]. However, pharmacophore models based on the role of Asp301 in substrate binding are compromised by reports of catalytic activity toward substrates devoid of a basic nitrogen, which have generally been ignored. We characterized a high-affinity ligand for P450 2D6, also devoid of a basic nitrogen atom, spirosulfonamide [4-[3-(4-fluorophenyl)-2-oxo-1-oxaspiro[4.4]non-3-en-4-yl]benzenesulfonamide], with K(s) 1.6 microM. Spirosulfonamide is a substrate for P450 2D6 (k(cat) 6.5 min(-)(1) for the formation of a syn spiromethylene carbinol, K(m) 7 microM). Mutation of Asp301 to neutral residues (Asn, Ser, Gly) did not substantially affect the binding of spirosulfonamide (K(s) 2.5-3.5 microM). However, the hydroxylation of spirosulfonamide was attenuated in these mutants to the same extent (90%) as for the classic nitrogenous substrate bufuralol, and the effect of the D301N substitution was manifested on k(cat) but not K(m). Analogues of spirosulfonamide were also evaluated as ligands and substrates. Analogues in which the sulfonamide moiety was modified to an amide, thioamide, methyl sulfone, or hydrogen were ligands with K(s) values of 1.7-32 microM. All were substrates, and the methyl sulfone analogue was oxidized to the syn spiromethylene carbinol analogue of the major spirosulfonamide product. The D301N mutation produced varying changes in the oxidation patterns of the spirosulfonamide analogues. The peptidometic ritonavir and the steroids progesterone and testosterone had been reported to be substrates for P450 2D6, but the affinities (K(s)) were unknown; these were estimated to be 1.2, 1.5, and 15 microM, respectively (cf. 6 microM for the classic substrate bufuralol). The results are consistent with a role of Asp301 other than electrostatic interaction with a positively charged ligand. H-Bonding or electrostatic interactions probably enhance binding of some substrates, but our results show that it is not required for all substrates and explain why predictive models fail to recognize the proclivity for many substrates, especially those containing no basic nitrogen.

Catalytic specificity of CYP2D isoforms in rat and human

In rats, six cytochrome P450 (P450) 2D isoforms have been genetically identified. Nonetheless, there is little evidence of catalytic properties of each CYP2D isoform. In this study, using recombinant CYP2D isoforms (rat CYP2D1, CYP2D2, CYP2D3, and CYP2D4 and human CYP2D6) or hepatic microsomes, we investigated the catalytic specificity toward bufuralol, debrisoquine, and propranolol, which are frequently used as CYP2D substrates. Bufuralol was oxidized to three metabolites by rat and human hepatic microsomes. 1'-Hydroxybufuralol was the major metabolite. 1'2'-Ethenylbufuralol, one of the others, was identified as a novel metabolite. The formation of 1'-hydroxybufuralol and 1'2'-ethenylbufuralol in hepatic microsomes was inhibited by anti-CYP2D antibody, suggesting that these metabolites were formed by CYP2D isoforms. All rat and human recombinant CYP2D isoforms possessed activity for the 1'-hydroxylation of bufuralol, indicating that this catalytic property was common to all CYP2D isoforms. However, the 1'2'-ethenylation of bufuralol was catalyzed only by rat CYP2D4 and human CYP2D6. Debrisoquine was oxidized to two metabolites, 3-hydroxydebrisoquine, and 4-hydroxydebrisoquine, by hepatic microsomes. Recombinant CYP2D2 and CYP2D6 had very high levels of activity for the 4-hydroxylation of debrisoquine with low K(m) values. Only CYP2D1 had a higher level of 3-hydroxylation than 4-hydroxylation activity. Propranolol 4-hydroxylation was catalyzed by CYP2D2, CYP2D4, and CYP2D6. The 7-hydroxylation of propranolol was catalyzed only by CYP2D2. In conclusion, in rats, bufuralol 1'2'-ethenylation activity was specific to CYP2D4 and debrisoquine 4-hydroxylation and propranolol 7-hydroxylation activities were specific to CYP2D2. These catalytic activities are useful as a probe for rat CYP2D isoforms.

4-Hydroxylation of debrisoquine by human CYP1A1 and its inhibition by quinidine and quinine

A panel of 15 recombinant cytochromes P450 expressed in human B-lymphoblastoid cells was used to study debrisoquine 4-hydroxylation. Both CYP2D6 and CYP1A1 carried out the reaction. The apparent K(m) (micromolar) and V(max) (picomoles per minute per picomole of P450) for CYP2D6 were 12.1 and 18.2 and for CYP1A1 were 23.1 and 15.2, respectively. CYP1A1 debrisoquine 4-hydroxylase was inhibited by the CYP1A1 inhibitor alpha-naphthoflavone and the CYP1A1 substrate 7-ethoxyresorufin. Additionally and surprisingly, this reaction was also inhibited by quinidine and quinine, with respective IC(50) values of 1.38 +/- 0.10 and 3.31 +/- 0.14 microM, compared with those for CYP2D6 debrisoquine 4-hydroxylase of 0.018 +/- 0.05 and 3.75 +/- 2.07 microM, respectively. Anti-CYP1A1 monoclonal antibody (mAb) 1-7-1 abolished CYP1A1 debrisoquine hydroxylase and anti-CYP2D6 mAb 50-1-3 eradicated CYP2D6 debrisoquine 4-hydroxylase. Three further CYP2D6-specific reactions were tested: dextromethorphan O-demethylation, bufuralol 1'-hydroxylation, and sparteine dehydrogenation. The CYP2D6 specificity, judged by the CYP2D6/CYP1A1 activity ratios was 18.5, 7.0, 6.0, and 1.6 for dextromethorphan, bufuralol, sparteine, and debrisoquine, respectively. Thus, debrisoquine is not a specific CYP2D6 substrate and quinidine is not a specific CYP2D6 inhibitor. These findings have significant implications for the conduct of in vitro drug metabolism inhibition studies and underscore the fallacy of "specific chemical inhibitors" of a supergene family of enzymes that have overlapping substrate specificities. The use of highly specific mAbs in such studies is mandated. It is unclear as yet whether these findings have implications for the relationship between CYP2D6 genotype and in vivo debrisoquine 4-hydroxylase activity.