Heterologous expression of cytochrome P450 2D6 mutants, electron transfer, and catalysis of bufuralol hydroxylation: the role of aspartate 301 in structural integrity

Cytochrome b 5 Binds Tightly to Several Human Cytochrome P450 Enzymes

Numerous studies have been reported in the past 50-plus years regarding the stimulatory role of cytochrome b ₅ (b ₅) in some, but not all, microsomal cytochrome P450 (P450) reactions with drugs and steroids. A missing element in most of these studies has been a sensitive and accurate measure of binding affinities of b ₅ with P450s. In the course of work with P450 17A1, we developed a fluorescent derivative of a human b ₅ site-directed mutant, Alexa 488-T70C-b ₅, that could be used in binding assays at sub-μM concentrations. Alexa 488-T70C-b ₅ bound to human P450s 1A2, 2B6, 2C8, 2C9, 2E1, 2S1, 4A11, 3A4, and 17A1, with estimated K _d values ranging from 2.5 to 61 nM. Only weak binding was detected with P450 2D6, and no fluorescence attenuation was observed with P450 2A6. All of the P450s that bound b ₅ have some reported activity stimulation except for P450 2S1. The affinity of P450 3A4 for b ₅ was decreased somewhat by the presence of a substrate or inhibitor. The fluorescence of a P450 3A4•Alexa 488-T70C-b ₅ complex was partially restored by titration with NADPH-P450 reductase (POR) (K _d,apparent 89 nM), suggesting the existence of a ternary P450 3A4-b ₅-POR complex, as observed previously with P450 17A1. Gel filtration evidence was also obtained for this ternary complex with P450 3A4. Overall, the results indicated that the affinity of b ₅ for many P450s is very high, and that ternary P450-b ₅-POR complexes are relevant in P450 3A4 reactions as opposed to a shuttle mechanism. SIGNIFICANCE STATEMENT: High-affinity binding of cytochrome b ₅ (b ₅) (K _d < 100 nM) was observed with many drug-metabolizing cytochrome P450 (P450) enzymes. There is some correlation of binding with reported stimulation, with several exceptions. Evidence is provided for a ternary P450 3A4-b ₅-NADPH-P450 reductase complex.

Deciphering Reaction Determinants of Altered-Activity CYP2D6 Variants by Well-Tempered Metadynamics Simulation and QM/MM Calculations

The xenobiotic metabolizing enzyme CYP2D6 is the P450 cytochrome family member with the highest rate of polymorphism. This causes changes in the enzyme activity and specificity, which can ultimately lead to adverse reactions during drug treatment. To avoid or lower CYP-related toxicity risks, prediction of the most likely positions within a molecule where a metabolic reaction might occur is paramount. In order to obtain accurate predictions, it is crucial to understand all phenomena within the active site of the enzyme that contribute to an efficient substrate recognition and the subsequent catalytic reaction together with their relative weight within the overall thermodynamic context. This study aims to define the weight of the driving forces upon the C-H bond activation within CYP2D6 wild-type and a clinically relevant allelic variant with increased activity (CYP2D6*53) featuring two amino acid mutations in close vicinity of the heme. First, we investigated the steric and electrostatic complementarity of the substrate bufuralol using well-tempered metadynamics simulations with the aim to obtain the free energy profiles for each site of metabolism (SoM) within the different active sites. Second, the stereoelectronic complementarity was determined for each SoM within the two different active-site environments. Relying on the well-tempered metadynamics simulation energy profiles of each SoM, we identified the binding mode that was closest to the preferred transition-state geometry for efficient C-H bond activation. The binding modes were then used as starting structures for the quantum mechanics/molecular mechanics calculations performed to quantify the corresponding activation barriers. Our results show the relevance of the steric component in orienting the SoM in an energetically accessible position toward the heme. However, the corresponding intrinsic reactivity and electronic complementarity within the active site must be accurately evaluated in order to obtain a meaningful reaction prediction, from which the predominant SoM can be determined. The F120I mutation lowered the activation barrier for the major site and one of the minor SoMs. However, it had an impact neither on the CYP2D6 enantioselectivity preference of the oxidation reaction nor on the stereoselectivity from the substrate point of view.

Human cytochrome P450 enzymes bind drugs and other substrates mainly through conformational-selection modes

Cytochrome P450 (P450) enzymes are major catalysts involved in the oxidations of most drugs, steroids, carcinogens, fat-soluble vitamins, and natural products. The binding of substrates to some of the 57 human P450s and other mammalian P450s is more complex than a two-state system and has been proposed to involve mechanisms such as multiple ligand occupancy, induced-fit, and conformational-selection. Here, we used kinetic analysis of binding with multiple concentrations of substrates and computational modeling of these data to discern possible binding modes of several human P450s. We observed that P450 2D6 binds its ligand rolapitant in a mechanism involving conformational-selection. P450 4A11 bound the substrate lauric acid via conformational-selection, as did P450 2C8 with palmitic acid. Binding of the steroid progesterone to P450 21A2 was also best described by a conformational-selection model. Hexyl isonicotinate binding to P450 2E1 could be described by either a conformational-selection or an induced-fit model. Simulation of the binding of the ligands midazolam, bromocriptine, testosterone, and ketoconazole to P450 3A4 was consistent with an induced-fit or a conformational-selection model, but the concentration dependence of binding rates for varying both P450 3A4 and midazolam concentrations revealed discordance in the parameters, indicative of conformational-selection. Binding of the P450s 2C8, 2D6, 3A4, 4A11, and 21A2 was best described by conformational-selection, and P450 2E1 appeared to fit either mode. These findings highlight the complexity of human P450-substrate interactions and that conformational-selection is a dominant feature of many of these interactions.

Rolapitant Is a Reversible Inhibitor of CYP2D6

Rolapitant [(Varubi), 5S,8S)-8-[[(1R)-1-[3,5 bis(trifluoromethyl phenyl]ethoxy]methyl]-8-phenyl-1,7-diazaspiro[4.5]decan-2-one] is a high-affinity NK1 receptor antagonist that was approved in September 2015 as a treatment for nausea and vomiting caused by chemotherapy. In vivo rolapitant moderately inhibits CYP2D6 for at least 7 days after one 180 mg dose. Due to the long inhibition time, we investigated rolapitant as a possible mechanism-based inactivator of CYP2D6. Rolapitant docked in the active site of CYP2D6 and displayed type I binding to CYP2D6 with a K _s value of 1.2 ± 0.4 µM. However, in NADPH-, time-, and concentration-dependent assays of CYP2D6 activity, no evidence for mechanism-based inactivation and no metabolites of rolapitant were observed. Stopped-flow binding studies yielded a k_on /k_off (K _d) value of 6.2 µM. The IC₅₀ value for rolapitant inhibition of CYP2D6 activity was 24 µM, suggesting that inhibition is not due to tight binding of rolapitant to CYP2D6. By Lineweaver-Burk analysis, rolapitant behaved as a mixed, reversible inhibitor. The K _i values of 20 and 34 µM were determined by Dixon analysis, with bufuralol and dextromethorphan as reporter substrates, respectively, and drug-drug interaction modeling did not predict the reported in vivo inhibition. The interaction of rolapitant with CYP2D6 was also examined in 1 microsecond molecular dynamics simulations. Rolapitant adopted multiple low-energy binding conformations near the active site, but at distances not consistent with metabolism. Given these findings, we do not see evidence that rolapitant is a mechanism-based inactivator. Moreover, the reversible inhibition of CYP2D6 by rolapitant may not fully account for the moderate inhibition described in vivo.

CYP2D6 Allelic Variants *34, *17-2, *17-3, and *53 and a Thr309Ala Mutant Display Altered Kinetics and NADPH Coupling in Metabolism of Bufuralol and Dextromethorphan and Altered Susceptibility to Inactivation by SCH 66712

Metabolic phenotype can be affected by multiple factors, including allelic variation and interactions with inhibitors. Human CYP2D6 is responsible for approximately 20% of cytochrome P450-mediated drug metabolism but consists of more than 100 known variants; several variants are commonly found in the population, whereas others are quite rare. Four CYP2D6 allelic variants-three with a series of mutations distal to the active site (*34, *17-2, *17-3) and one ultra-metabolizer with mutations near the active site (*53), along with reference *1 and an active site mutant of *1 (Thr309Ala)-were expressed, purified, and studied for interactions with the typical substrates dextromethorphan and bufuralol and the inactivator SCH 66712. We found that *34, *17-2, and *17-3 displayed reduced enzyme activity and NADPH coupling while producing the same metabolites as *1, suggesting a possible role for Arg296 in NADPH coupling. A higher-activity variant, *53, displayed similar NADPH coupling to *1 but was less susceptible to inactivation by SCH 66712. The Thr309Ala mutant showed similar activity to that of *1 but with greatly reduced NADPH coupling. Overall, these results suggest that kinetic and metabolic analysis of individual CYP2D6 variants is required to understand their possible contributions to variable drug response and the complexity of personalized medicine.

Sulfenylation of Human Liver and Kidney Microsomal Cytochromes P450 and Other Drug-Metabolizing Enzymes as a Response to Redox Alteration

The lumen of the endoplasmic reticulum (ER) provides an oxidizing environment to aid in the formation of disulfide bonds, which is tightly regulated by both antioxidant proteins and small molecules. On the cytoplasmic side of the ER, cytochrome P450 (P450) proteins have been identified as a superfamily of enzymes that are important in the formation of endogenous chemicals as well as in the detoxication of xenobiotics. Our previous report described oxidative inhibition of P450 Family 4 enzymes via oxidation of the heme-thiolate cysteine to a sulfenic acid (-SOH) (Albertolle, M. E. et al. (2017) J. Biol. Chem. 292, 11230-11242). Further proteomic analyses of murine kidney and liver microsomes led to the finding that a number of other drug-metabolizing enzymes located in the ER are also redox-regulated in this manner. We expanded our analysis of sulfenylated enzymes to human liver and kidney microsomes. Evaluation of the sulfenylation, catalytic activity, and spectral properties of P450s 1A2, 2C8, 2D6, and 3A4 led to the identification of two classes of redox sensitivity in P450 enzymes: heme-thiolate-sensitive and thiol-insensitive. These findings provide evidence for a mammalian P450 regulatory mechanism, which may also be relevant to other drug-metabolizing enzymes. (Data are available via ProteomeXchange with identifier PXD007913.).

Heme-thiolate sulfenylation of human cytochrome P450 4A11 functions as a redox switch for catalytic inhibition

Cytochrome P450 (P450, CYP) 4A11 is a human fatty acid ω-hydroxylase that catalyzes the oxidation of arachidonic acid to the eicosanoid 20-hydroxyeicosatetraenoic acid (20-HETE), which plays important roles in regulating blood pressure regulation. Variants of P450 4A11 have been associated with high blood pressure and resistance to anti-hypertensive drugs, and 20-HETE has both pro- and antihypertensive properties relating to increased vasoconstriction and natriuresis, respectively. These physiological activities are likely influenced by the redox environment, but the mechanisms are unclear. Here, we found that reducing agents (e.g. dithiothreitol and tris(2-carboxyethyl)phosphine) strongly enhanced the catalytic activity of P450 4A11, but not of 10 other human P450s tested. Conversely, added H₂O₂ attenuated P450 4A11 catalytic activity. Catalytic roles of five of the potentially eight implicated Cys residues of P450 4A11 were eliminated by site-directed mutagenesis. Using an isotope-coded dimedone/iododimedone-labeling strategy and mass spectrometry of peptides, we demonstrated that the heme-thiolate cysteine (Cys-457) is selectively sulfenylated in an H₂O₂ concentration-dependent manner. This sulfenylation could be reversed by reducing agents, including dithiothreitol and dithionite. Of note, we observed heme ligand cysteine sulfenylation of P450 4A11 ex vivo in kidneys and livers derived from CYP4A11 transgenic mice. We also detected sulfenylation of murine P450 4a12 and 4b1 heme peptides in kidneys. To our knowledge, reversible oxidation of the heme thiolate has not previously been observed in P450s and may have relevance for 20-HETE-mediated functions.

The catalytic function of cytochrome P450 is entwined with its membrane-bound nature

Cytochrome P450, a family of monooxygenase enzymes, is organized as a catalytic metabolon, which requires enzymatic partners as well as environmental factors that tune its complex dynamic. P450 and its reducing counterparts—cytochrome P450-reductase and cytochrome b ₅—are membrane-bound proteins located in the cytosolic side of the endoplasmic reticulum. They are believed to dynamically associate to form functional complexes. Increasing experimental evidence signifies the role(s) played by both protein-protein and protein-lipid interactions in P450 catalytic function and efficiency. However, the biophysical challenges posed by their membrane-bound nature have severely limited high-resolution understanding of the molecular interfaces of these interactions. In this article, we provide an overview of the current knowledge on cytochrome P450, highlighting the environmental factors that are entwined with its metabolic function. Recent advances in structural biophysics are also discussed, setting up the bases for a new paradigm in the study of this important class of membrane-bound enzymes.

Aromatic hydroxylation of salicylic acid and aspirin by human cytochromes P450

Aspirin (acetylsalicylic acid) is a well-known and widely-used analgesic. It is rapidly deacetylated to salicylic acid, which forms two hippuric acids-salicyluric acid and gentisuric acid-and two glucuronides. The oxidation of aspirin and salicylic acid has been reported with human liver microsomes, but data on individual cytochromes P450 involved in oxidation is lacking. In this study we monitored oxidation of these compounds by human liver microsomes and cytochrome P450 (P450) using UPLC with fluorescence detection. Microsomal oxidation of salicylic acid was much faster than aspirin. The two oxidation products were 2,5-dihydroxybenzoic acid (gentisic acid, documented by its UV and mass spectrum) and 2,3-dihydroxybenzoic acid. Formation of neither product was inhibited by desferrioxamine, suggesting a lack of contribution of oxygen radicals under these conditions. Although more liphophilic, aspirin was oxidized less efficiently, primarily to the 2,5-dihydroxy product. Recombinant human P450s 2C8, 2C9, 2C19, 2D6, 2E1, and 3A4 all catalyzed the 5-hydroxylation of salicylic acid. Inhibitor studies with human liver microsomes indicated that all six of the previously mentioned P450s could contribute to both the 5- and 3-hydroxylation of salicylic acid and that P450s 2A6 and 2B6 have contributions to 5-hydroxylation. Inhibitor studies indicated that the major human P450 involved in both 3- and 5-hydroxylation of salicylic acid is P450 2E1.

Mechanistic basis of electron transfer to cytochromes p450 by natural redox partners and artificial donor constructs

Cytochromes P450 (P450s) are hemoproteins catalyzing oxidative biotransformation of a vast array of natural and xenobiotic compounds. Reducing equivalents required for dioxygen cleavage and substrate hydroxylation originate from different redox partners including diflavin reductases, flavodoxins, ferredoxins and phthalate dioxygenase reductase (PDR)-type proteins. Accordingly, circumstantial analysis of structural and physicochemical features governing donor-acceptor recognition and electron transfer poses an intriguing challenge. Thus, conformational flexibility reflected by togging between closed and open states of solvent exposed patches on the redox components was shown to be instrumental to steered electron transmission. Here, the membrane-interactive tails of the P450 enzymes and donor proteins were recognized to be crucial to proper orientation toward each other of surface sites on the redox modules steering functional coupling. Also, mobile electron shuttling may come into play. While charge-pairing mechanisms are of primary importance in attraction and complexation of the redox partners, hydrophobic and van der Waals cohesion forces play a minor role in docking events. Due to catalytic plasticity of P450 enzymes, there is considerable promise in biotechnological applications. Here, deeper insight into the mechanistic basis of the redox machinery will permit optimization of redox processes via directed evolution and DNA shuffling. Thus, creation of hybrid systems by fusion of the modified heme domain of P450s with proteinaceous electron carriers helps obviate the tedious reconstitution procedure and induces novel activities. Also, P450-based amperometric biosensors may open new vistas in pharmaceutical and clinical implementation and environmental monitoring.

Metoclopramide is metabolized by CYP2D6 and is a reversible inhibitor, but not inactivator, of CYP2D6

1. Metoclopramide is a widely used clinical drug in a variety of medical settings with rare acute dystonic events reported. The aim of this study was to assess a previous report of inactivation of CYP2D6 by metoclopramide, to determine the contribution of various CYPs to metoclopramide metabolism, and to identify the mono-oxygenated products of metoclopramide metabolism. 2. Metoclopramide interacted with CYP2D6 with Type I binding and a Ks value of 9.56 ± 1.09 µM. CYP2D6 was the major metabolizer of metoclopramide and the two major products were N-deethylation of the diethyl amine and N-hydroxylation on the phenyl ring amine. CYPs 1A2, 2C9, 2C19, and 3A4 also metabolized metoclopramide. 3. While reversible inhibition of CYP2D6 was noted, CYP2D6 inactivation by metoclopramide was not observed under conditions of varying concentration or varying time using Supersomes(TM) or pooled human liver microsomes. 4. The major metabolites of metoclopramide were N-hydroxylation and N-deethylation formed most efficiently by CYP2D6 but also formed by all CYPs examined. Also, while metoclopramide is metabolized primarily by CYP2D6, it is not a mechanism-based inactivator of CYP2D6 in vitro.

Functional influence of human CYP2D6 allelic variations: P34S, E418K, S486T, and R296C

CYP2D6 is responsible for the oxidative metabolism of 20-25 % of clinical drugs and its genetic polymorphisms can significantly influence the drug metabolism. In this study, we analyzed the functional activities of four nonsynonymous single nucleotide polymorphisms from CYP2D6*52 allele, which were recently found, and one found frequently in CYP2D6 alleles. Recombinant variant enzymes of E418K, S486T, and R296C were successfully expressed in Escherichia coli and purified. However, a CYP holoenzyme spectrum of P34S variant was not detected in E. coli whole cell level. Structural analysis indicated that P34S mutation seemed to perturb a highly conserved proline-rich N-terminus of CYP2D6. Steady state kinetic analyses showed the significant reductions of enzymatic activities in E418K and R296C variants. In the case of bufuralol 1'-hydroxylation, a novel mutant, E418K, showed 32 % decrease in catalytic efficiency (k cat/K m) mainly due to the decrease of k cat value. R296C showed much greater reduction in the catalytic efficiency (9 % of wild-type) due to both of a decrease of k cat value and an increase of K m value. In the case of dextromethorphan O-demethylation, E418K showed both of a decrease of k cat value and an increase K m value to result in ~43 % reduction of catalytic efficiency. A highly decreased catalytic efficiency (~6 % of wild-type) in the mutant of R296C also was observed mainly due to the dramatic change of k cat value of dextromethorphan O-demethylation. These results suggested that individuals carrying these allelic variants are likely to have the altered metabolic abilities of many clinical drugs therefore, these polymorphisms of CYP2D6 should be much concerned for reliable drug treatment.

Metabolism of 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine by mitochondrion-targeted cytochrome P450 2D6: implications in Parkinson disease

1-Methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP) is a neurotoxic side product formed in the chemical synthesis of desmethylprodine opioid analgesic, which induces Parkinson disease. Monoamine oxidase B, present in the mitochondrial outer membrane of glial cells, catalyzes the oxidation of MPTP to the toxic 1-methyl-4-phenylpyridinium ion (MPP(+)), which then targets the dopaminergic neurons causing neuronal death. Here, we demonstrate that mitochondrion-targeted human cytochrome P450 2D6 (CYP2D6), supported by mitochondrial adrenodoxin and adrenodoxin reductase, can efficiently catalyze the metabolism of MPTP to MPP(+), as shown with purified enzymes and also in cells expressing mitochondrial CYP2D6. Neuro-2A cells stably expressing predominantly mitochondrion-targeted CYP2D6 were more sensitive to MPTP-mediated mitochondrial respiratory dysfunction and complex I inhibition than cells expressing predominantly endoplasmic reticulum-targeted CYP2D6. Mitochondrial CYP2D6 expressing Neuro-2A cells produced higher levels of reactive oxygen species and showed abnormal mitochondrial structures. MPTP treatment also induced mitochondrial translocation of an autophagic marker, Parkin, and a mitochondrial fission marker, Drp1, in differentiated neurons expressing mitochondrial CYP2D6. MPTP-mediated toxicity in primary dopaminergic neurons was attenuated by CYP2D6 inhibitor, quinidine, and also partly by monoamine oxidase B inhibitors deprenyl and pargyline. These studies show for the first time that dopaminergic neurons expressing mitochondrial CYP2D6 are fully capable of activating the pro-neurotoxin MPTP and inducing neuronal damage, which is effectively prevented by the CYP2D6 inhibitor quinidine.

Crystal structure of human cytochrome P450 2D6 with prinomastat bound

Human cytochrome P450 2D6 contributes to the metabolism of >15% of drugs used in clinical practice. This study determined the structure of P450 2D6 complexed with a substrate and potent inhibitor, prinomastat, to 2.85 Å resolution by x-ray crystallography. Prinomastat binding is well defined by electron density maps with its pyridyl nitrogen bound to the heme iron. The structure of ligand-bound P450 2D6 differs significantly from the ligand-free structure reported for the P450 2D6 Met-374 variant (Protein Data Bank code 2F9Q). Superposition of the structures reveals significant differences for β sheet 1, helices A, F, F', G", G, and H as well as the helix B-C loop. The structure of the ligand complex exhibits a closed active site cavity that conforms closely to the shape of prinomastat. The closure of the open cavity seen for the 2F9Q structure reflects a change in the direction and pitch of helix F and introduction of a turn at Gly-218, which is followed by a well defined helix F' that was not observed in the 2F9Q structure. These differences reflect considerable structural flexibility that is likely to contribute to the catalytic versatility of P450 2D6, and this new structure provides an alternative model for in silico studies of substrate interactions with P450 2D6.

Molecular dynamics analysis reveals structural insights into mechanism of nicotine N-demethylation catalyzed by tobacco cytochrome P450 mono-oxygenase

CYP82E4, a cytochrome P450 monooxygenase, has nicotine N-demethylase (NND) activity, which mediates the bioconversion of nicotine into nornicotine in senescing tobacco leaves. Nornicotine is a precursor of the carcinogen, tobacco-specific nitrosamine. CYP82E3 is an ortholog of CYP82E4 with 95% sequence identity, but it lacks NND activity. A recent site-directed mutagenesis study revealed that a single amino acid substitution, i.e., cysteine to tryptophan at the 330 position in the middle of protein, restores the NND activity of CYP82E3 entirely. However, the same amino acid change caused the loss of the NND activity of CYP82E4. To determine the mechanism of the functional turnover of the two molecules, four 3D structures, i.e., the two molecules and their corresponding cys–trp mutants were modeled. The resulting structures exhibited that the mutation site is far from the active site, which suggests that no direct interaction occurs between the two sites. Simulation studies in different biological scenarios revealed that the mutation introduces a conformation drift with the largest change at the F-G loop. The dynamics trajectories analysis using principal component analysis and covariance analysis suggests that the single amino acid change causes the opening and closing of the transfer channels of the substrates, products, and water by altering the motion of the F-G and B-C loops. The motion of helix I is also correlated with the motion of both the F-G loop and the B-C loop and; the single amino acid mutation resulted in the curvature of helix I. These results suggest that the single amino acid mutation outside the active site region may have indirectly mediated the flexibility of the F-G and B-C loops through helix I, causing a functional turnover of the P450 monooxygenase.

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).

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.

Anandamide oxidation by wild-type and polymorphically expressed CYP2B6 and CYP2D6

Anandamide is an arachidonic acid-derived endogenous cannabinoid that regulates normal physiological functions and pathophysiological responses within the central nervous system and in the periphery. Several cytochrome P450 (P450) isoforms metabolize anandamide to form hydroxylated and epoxygenated products. Human CYP2B6 and CYP2D6, which are expressed heterogeneously throughout the brain, exhibit clinically significant polymorphisms and are regulated by external factors, such as alcohol and smoking. Oxidative metabolism of anandamide by these two P450s may have important functional consequences for endocannabinoid system signaling. In this study, we investigated the metabolism of anandamide by wild-type CYP2B6 (2B6.1) and CYP2D6 (2D6.1) and by their common polymorphic mutants 2B6.4, 2B6.6, 2B6.9, and 2D6.34. Major differences in anandamide metabolism by the two isoforms and their mutants were found in vitro with respect to the formation of 20-hydroxyeicosatetraenoic acid ethanolamide (20-HETE-EA) and 14,15-epoxyeicosatetraenoic acid ethanolamide (14,15-EET-EA). Pharmacological studies showed that both 20-HETE-EA and 14,15-EET-EA bind to the rat brain cannabinoid CB1 receptor with lower affinities relative to that of anandamide. In addition, both products are degraded more rapidly than anandamide in rat brain homogenates. Their degradation occurs via different mechanisms involving either fatty acid amide hydrolase (FAAH), the major anandamide-degrading enzyme, or epoxide hydrolase (EH). Thus, the current findings provide potential new insights into the actions of inhibitors FAAH and EH, which are being developed as novel therapeutic agents, as well as a better understanding of the interactions between the cytochrome P450 monooxygenases and the endocannabinoid system.

Nitric oxide blocks cellular heme insertion into a broad range of heme proteins

Although the insertion of heme into proteins enables their function in bioenergetics, metabolism, and signaling, the mechanisms and regulation of this process are not fully understood. We developed a means to study cellular heme insertion into apo-protein targets over a 3-h period and then investigated how nitric oxide (NO) released from a chemical donor (NOC-18) might influence heme (protoporphyrin IX) insertion into seven targets that present a range of protein structures, heme ligation states, and functions (three NO synthases, two cytochrome P450's, catalase, and hemoglobin). NO blocked cellular heme insertion into all seven apo-protein targets. The inhibition occurred at relatively low (nM/min) fluxes of NO, was reversible, and did not involve changes in intracellular heme levels, activation of guanylate cyclase, or inhibition of mitochondrial ATP production. These aspects and the range of protein targets suggest that NO can act as a global inhibitor of heme insertion, possibly by inhibiting a common step in the process.

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.

Human liver mitochondrial cytochrome P450 2D6--individual variations and implications in drug metabolism

Constitutively expressed human cytochrome P450 2D6 (CYP2D6; EC 1.14.14.1) is responsible for the metabolism of approximately 25% of drugs in common clinical use. It is widely accepted that CYP2D6 is localized in the endoplasmic reticulum of cells; however, we have identified this enzyme in the mitochondria of human liver samples and found that extensive inter-individual variability exists with respect to the level of the mitochondrial enzyme. Metabolic assays using 7-methoxy-4-aminomethylcoumarin as a substrate show that the human liver mitochondrial enzyme is capable of oxidizing this substrate and that the catalytic activity is supported by mitochondrial electron transfer proteins. In the present study, we show that CYP2D6 contains an N-terminal chimeric signal that mediates its bimodal targeting to the endoplasmic reticulum and mitochondria. In vitro mitochondrial import studies using both N-terminal deletions and point mutations suggest that the mitochondrial targeting signal is localized between residues 23-33 and that the positively-charged residues at positions 24, 25, 26, 28 and 32 are required for mitochondrial targeting. The importance of the positively-charged residues was confirmed by transient transfection of a CYP2D6 mitochondrial targeting signal mutant in COS-7 cells. Both the mitochondria and the microsomes from a CYP2D6 stable expression cell line contain the enzyme and both fractions exhibit bufuralol 1'-hydroxylation activity, which is completely inhibited by CYP2D6 inhibitory antibody. Overall, these results suggest that the targeting of CYP2D6 to mitochondria could be an important physiological process that has significance in xenobiotic metabolism.

Versatile capacity of shuffled cytochrome P450s for dye production

DNA family shuffling is a relatively new method of directed evolution used to create novel enzymes in order to improve their existing properties or to develop new features. This method of evolution in vitro has one basic requirement: a high similarity of initial parental sequences. Cytochrome P450 enzymes are relatively well conserved in their amino acid sequences. Members of the same family can have more than 40% of sequence identity at the protein level and are therefore good candidates for DNA family shuffling. These xenobiotic-metabolising enzymes have an ability to metabolise a wide range of chemicals and produce a variety of products including blue pigments such as indigo. By applying the specifically designed DNA family shuffling approach, catalytic properties of cytochrome P450 enzymes were further extended in the chimeric progeny to include a new range of blue colour formations. This mini-review evokes the possibility of exploiting directed evolution of cytochrome P450s and the novel enzymes created by DNA family shuffling for the production of new dyes.

The endocannabinoid anandamide is a substrate for the human polymorphic cytochrome P450 2D6

Members of the cytochrome P450 (P450) family of drug-metabolizing enzymes are present in the human brain, and they may have important roles in the oxidation of endogenous substrates. The polymorphic CYP2D6 is one of the major brain P450 isoforms and has been implicated in neurodegeneration, psychosis, schizophrenia, and personality traits. The objective of this study was to determine whether the endocannabinoid arachidonoylethanolamide (anandamide) is a substrate for CYP2D6. Anandamide is the endogenous ligand to the cannabinoid receptor CB1, which is also activated by the main psychoactive component in marijuana. Signaling via the CB1 receptor alters sensory and motor function, cognition, and emotion. Recombinant CYP2D6 converted anandamide to 20-hydroxyeicosatetraenoic acid ethanolamide and 5,6-, 8,9-, 11,12-, and 14,15-epoxyeicosatrienoic acid ethanolamides (EET-EAs) with low micromolar K(m) values. CYP2D6 further metabolized the epoxides of anandamide to form novel dioxygenated derivatives. Human brain microsomal and mitochondrial preparations metabolized anandamide to form hydroxylated and epoxygenated products, respectively. An inhibitory antibody against CYP2D6 significantly decreased the mitochondrial formation of the EET-EAs. To our knowledge, anandamide and its epoxides are the first eicosanoid-like molecules to be identified as CYP2D6 substrates. Our study suggests that anandamide may be a physiological substrate for brain mitochondrial CYP2D6, implicating this polymorphic enzyme as a potential component of the endocannabinoid system in the brain. This study also offers support to the hypothesis that neuropsychiatric phenotype differences among individuals with genetic variations in CYP2D6 could be ascribable to interactions of this enzyme with endogenous substrates.

Analysis of CYP2D6 substrate interactions by computational methods

Cytochrome P450 CYP2D6 is involved in the oxidation of well over 150 drugs and, in general, those which contain a basic nitrogen atom in the molecule. To clarify how the residues of CYP2D6 are utilized for orientating a wide range of its specific substrates and distinguishing them from a variety of other organic compounds, docking studies by AutoDock and molecular dynamics (MD) simulations were conducted. Specific ligands were docked to both the homology model and crystal structures optimally to estimate the site of reaction on the ligand molecule and the binding energy for the complex, which were generally in good agreement with the experimental data. MD simulation for the CYP2D6-propranolol complex was then carried out to reveal the amino acid residues interacting with the substrate at the active site. Phe-120, Glu-216, Asp-301, and Phe-483 are identified as the substrate-binding residues in agreement with previously reported site-directed mutagenesis data and the crystal structure reported recently (PDB code: 2F9Q). As well as these residues, our theoretical prediction suggests that Phe-219 and Glu-222 are also important residues for mediating oxidation of substrates, especially propranolol.

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.

Effects of N-terminal modification of recombinant human cytochrome P450 1A2 on catalytic activity

1. The high-level expression of mammalian cytochrome P450 in bacteria usually requires modification of the amino-terminal region of the enzyme. The effect of altering amino acids in the N-terminus of human recombinant CYP1A2 on its catalytic activity was investigated herein. 2. Rates of 7-ethoxyresorufin O-deethylation by CYP1A2a (a form made by altering the amino acids LLL of CYP1A2 to RER at positions 3-5) in reconstituted systems were significantly low compared with those of other CYP1A2 N-terminal variants at a low ratio of cytochrome P450 to NADPH-cytochrome P450 reductase, but not at higher reductase concentrations. 3. CYP1A2a-dependent ethoxyresorufin O-deethylase activity in a cumene hydroperoxide-supported system was approximately 2-fold higher than other CYP1A2 N-terminal variants. 4. Our results suggest that modification of three N-terminal amino acids in CYP1A2 alters the interaction between CYP1A2 and the reductase in reconstituted phospholipid vesicles and in the bicistronic membranes.

Functional interaction of nitrogenous organic bases with cytochrome P450: A critical assessment and update of substrate features and predicted key active-site elements steering the access, binding, and orientation of amines

The widespread use of nitrogenous organic bases as environmental chemicals, food additives, and clinically important drugs necessitates precise knowledge about the molecular principles governing biotransformation of this category of substrates. In this regard, analysis of the topological background of complex formation between amines and P450s, acting as major catalysts in C- and N-oxidative attack, is of paramount importance. Thus, progress in collaborative investigations, combining physico-chemical techniques with chemical-modification as well as genetic engineering experiments, enables substantiation of hypothetical work resulting from the design of pharmacophores or homology modelling of P450s. Based on a general, CYP2D6-related construct, the majority of prospective amine-docking residues was found to cluster near the distal heme face in the six known SRSs, made up by the highly variant helices B', F and G as well as the N-terminal portion of helix C and certain beta-structures. Most of the contact sites examined show a frequency of conservation < 20%, hinting at the requirement of some degree of conformational versatility, while a limited number of amino acids exhibiting a higher level of conservation reside close to the heme core. Some key determinants may have a dual role in amine binding and/or maintenance of protein integrity. Importantly, a series of non-SRS elements are likely to be operative via long-range effects. While hydrophobic mechanisms appear to dominate orientation of the nitrogenous compounds toward the iron-oxene species, polar residues seem to foster binding events through H-bonding or salt-bridge formation. Careful uncovering of structure-function relationships in amine-enzyme association together with recently developed unsupervised machine learning approaches will be helpful in both tailoring of novel amine-type drugs and early elimination of potentially toxic or mutagenic candidates. Also, chimeragenesis might serve in the construction of more efficient P450s for activation of amine drugs and/or bioremediation.

Crystal Structure of Human Cytochrome P450 2D6

Cytochrome P450 2D6 is a heme-containing enzyme that is responsible for the metabolism of at least 20% of known drugs. Substrates of 2D6 typically contain a basic nitrogen and a planar aromatic ring. The crystal structure of human 2D6 has been solved and refined to 3.0A resolution. The structure shows the characteristic P450 fold as seen in other members of the family, with the lengths and orientations of the individual secondary structural elements being very similar to those seen in 2C9. There are, however, several important differences, the most notable involving the F helix, the F-G loop, the B'helix, beta sheet 4, and part of beta sheet 1, all of which are situated on the distal face of the protein. The 2D6 structure has a well defined active site cavity above the heme group, containing many important residues that have been implicated in substrate recognition and binding, including Asp-301, Glu-216, Phe-483, and Phe-120. The crystal structure helps to explain how Asp-301, Glu-216, and Phe-483 can act as substrate binding residues and suggests that the role of Phe-120 is to control the orientation of the aromatic ring found in most substrates with respect to the heme. The structure has been compared with published homology models and has been used to explain much of the reported site-directed mutagenesis data and help understand the metabolism of several compounds.

Structural motif-based homology modeling of CYP27A1 and site-directed mutational analyses affecting vitamin D hydroxylation

Human CYP27A1 is a mitochondrial cytochrome P450, which is principally found in the liver and plays important roles in the biological activation of vitamin D(3) and in the biosynthesis of bile acids. We have applied a systematic analysis of hydrogen bonding patterns in 11 prokaryotic and mammalian CYP crystal structures to construct a homology-based model of CYP27A1. Docking of vitamin D(3) structures into the active site of this model identified potential substrate contact residues in the F-helix, the beta-3 sheet, and the beta-5 sheet. Site-directed mutagenesis and expression in COS-1 cells confirmed that these positions affect enzymatic activity, in some cases shifting metabolism of 1alpha-hydroxyvitamin D(3) to favor 25- or 27-hydroxylation. The results suggest that conserved hydrophobic residues in the beta-5 hairpin help define the shape of the substrate binding cavity and that this structure interacts with Phe-248 in the F-helix. Mutations directed toward the beta-3a strand suggested a possible heme-binding interaction centered on Asn-403 and a structural role for substrate contact residues Thr-402 and Ser-404.

Interactions of mammalian cytochrome P450, NADPH-cytochrome P450 reductase, and cytochrome b(5) enzymes

An immobilized system was developed to detect interactions of human cytochromes P450 (P450) with the accessory proteins NADPH-P450 reductase and cytochrome b(5) (b(5)) using an enzyme-linked affinity approach. Purified enzymes were first bound to wells of a polystyrene plate, and biotinylated partner enzymes were added and bound. A streptavidin-peroxidase complex was added, and protein-protein binding was monitored by measuring peroxidase activity of the bound biotinylated proteins. In a model study, we examined protein-protein interactions of Pseudomonas putida putidaredoxin (Pdx) and putidaredoxin reductase (PdR). A linear relationship (r(2)=0.96) was observed for binding of PdR-biotin to immobilized Pdx compared with binding of Pdx-biotin to immobilized PdR (the estimated K(d) value for the Pdx.PdR complex was 0.054muM). Human P450 2A6 interacted strongly with NADPH-P450 reductase; the K(d) values (with the reductase) ranged between 0.005 and 0.1muM for P450s 2C19, 2D6, and 3A4. Relatively weak interaction was found between holo-b(5) or apo-b(5) (devoid of heme) with NADPH-P450 reductase. Among the rat, rabbit, and human P450 1A2 enzymes, the rat enzyme showed the tightest interaction with b(5), although no increases in 7-ethoxyresorufin O-deethylation activities were observed with any of the P450 1A2 enzymes. Human P450s 2A6, 2D6, 2E1, and 3A4 interacted well with b(5), with P450 3A4 yielding the lowest K(d) values followed by P450s 2A6 and 2D6. No appreciable increases in interaction between human P450s with b(5) or NADPH-P450 reductase were observed when typical substrates for the P450s were included. We also found that NADPH-P450 reductase did not cause changes in the P450.substrate K(d) values estimated from substrate-induced UV-visible spectral changes with rabbit P450 1A2 or human P450 2A6, 2D6, or 3A4. Collectively, the results show direct and tight interactions between P450 enzymes and the accessory proteins NADPH-P450 reductase and b(5), with different affinities, and that ligand binding to mammalian P450s did not lead to increased interaction between P450s and the reductase.

Assessment of arginine 97 and lysine 72 as determinants of substrate specificity in cytochrome P450 2C9 (CYP2C9)

CYP2C9 is distinguished by a preference for substrates bearing a negative charge at physiological pH. Previous studies have suggested that CYP2C9 residues R97 and K72 may play roles in determining preference for anionic substrates by interaction at the active site or in the access channel. The aim of the present study was to assess the role of these two residues in determining substrate selectivity. R97 and K72 were substituted with negative, uncharged polar and hydrophobic residues using a degenerate polymerase chain reaction-directed strategy. Wild-type and mutant enzymes were expressed in bicistronic format with human cytochrome P450 reductase in Escherichia coli. Mutation of R97 led to a loss of holoenzyme expression for R97A, R97V, R97L, R97T, and R97E mutants. Low levels of hemoprotein were detected for R97Q, R97K, R97I, and R97P mutants. Significant apoenzyme was observed, suggesting that heme insertion or protein stability was compromised in R97 mutants. These observations are consistent with a structural role for R97 in addition to any role in substrate binding. By contrast, all K72 mutants examined (K72E, K72Q, K72V, and K72L) could be expressed as hemoprotein at levels comparable to wild-type. Type I binding spectra were obtained with wild-type and K72 mutants using diclofenac and ibuprofen. Mutation of K72 had little or no effect on the interaction with these substrates, arguing against a critical role in determining substrate specificity. Thus, neither residue appears to play a role in determining substrate specificity, but a structural role for R97 can be proposed consistent with recently published crystallographic data for CYP2C9 and CYP2C5.

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.

Non-specific inhibition of human cytochrome P450-catalyzed reactions by hemin

Hemin, a stable form of heme, is known to have an antimutagenic effect. Inhibitory effects of hemin on the cytochrome P450 (CYP)-catalyzed reactions of human liver microsomes and reconstituted systems containing purified CYP and NADPH-cytochrome P450 reductase (NPR) were seen. Hemin non-specifically inhibited all of the microsomal CYP activities examined. Hemin also inhibited 7-ethoxyresorufin O-deethylation, 3-[2-(N,N-diethyl-N-methylammonium)ethyl]-7-methoxy-4-methylcoumarin O-demethylation, and testosterone 6beta-hydroxylation catalyzed by purified CYPs 1A2, 2D6, and 3A4, with IC50 values of 27, 19, and 2.4 microM, respectively. Hemin also inhibited reduction of cytochrome c and ferricyanide by NPR, as much as 47%. Spectrally detectable CYP was destroyed in human liver microsomes and in a reconstituted system in the presence of hemin and an NADPH-generating system. We propose that the antimutagenic effect of hemin might be due to inhibition of CYP and NPR enzymes involved in the bioactivation of mutagens.

Phe120 contributes to the regiospecificity of cytochrome P450 2D6: mutation leads to the formation of a novel dextromethorphan metabolite

Although the residues that determine the preference of CYP2D6 (cytochrome P450 2D6) for compounds containing a basic nitrogen are well characterized, the contribution of other active site residues to substrate binding and orientation is less well understood. Our structural model of CYP2D6 identifies the aromatic residue Phe120 as a likely major feature of the active site. To examine the role of Phe120, mutants of CYP2D6 in which this residue has been substituted by alanine, leucine, tyrosine, serine, histidine, tryptophan or methionine residues have been prepared in bacterial membranes co-expressing human cytochrome NADPH cytochrome P450 oxidoreductase. The mutants have been characterized using the prototypical bufuralol 1' hydroxylase and dextromethorphan O- and N-demethylase activities of CYP2D6. Larger effects on K(m) values are observed for dextromethorphan O-demethylation than for bufuralol 1' hydroxylation, indicating that the Phe120 side chain is more important in dextromethorphan than in bufuralol binding. A role for this side chain in determining the regiospecificity of substrate oxidation was indicated by changes in the relative rates of O- and N-demethylation of dextromethorphan and, notably, by the formation of 7-hydroxy dextromethrophan, a novel dextromethorphan metabolite, in mutants in which it had been substituted. Computational studies of dextromethorphan binding to the active site of the Phe120-->Ala mutant were carried out to throw light on the way in which the removal of this side chain leads to different modes of ligand binding.

Mutagenesis and molecular dynamics suggest structural and functional roles for residues in the N-terminal portion of the cytochrome P450 2B1 I helix

To investigate their potential roles in ligand access, binding, and subsequent metabolism, residues in the N-terminal portion of the cytochrome P450 2B1 I helix were mutated to alanine and phenylalanine. Of the 18 mutants from E286 to S294 only 7 yielded holoprotein in an Escherichia coli expression system. Substitutions at positions 289, 290, 292, and 294 caused >/= 2-fold changes in kcat and/or Km for two or more of the 2B1 substrates examined, testosterone, 7-ethoxy-4-trifluoromethylcoumarin, 7-benzyloxyresorufin, and benzphetamine. I290 substitutions had the largest effects on steady-state parameters for three substrates and increased benzphetamine affinity. Steered molecular dynamics simulations of testosterone egress along the I helix identified hydrophobic interactions with I290, L293, and S294 and water bridges to E286 and S294. Sensitivity of holoprotein formation to substitution and effects on substrate binding and metabolism suggest structural and functional roles for residues in the N-terminus of the cytochrome P450 2B1 I helix.

Functional interaction of cytochrome P450 with its redox partners: a critical assessment and update of the topology of predicted contact regions

The problem of donor-acceptor recognition has been the most important and intriguing one in the area of P450 research. The present review outlines the topological background of electron-transfer complex formation, showing that the progress in collaborative investigations, combining physical techniques with chemical-modification and immunolocalization studies as well as site-directed mutagenesis experiments, has increasingly enabled the substantiation of hypothetical work resulting from homology modelling of P450s. Circumstantial analysis reveals the contact regions for redox proteins to cluster on the proximal face of P450s, constituting parts of the highly conserved, heme-binding core fold. However, more variable structural components located in the periphery of the hemoprotein molecules also participate in donor docking. The cross-reactivity of electron carriers, purified from pro- and eukaryotic sources, with a diversity of P450 species points at a possible evolutionary conservation of common anchoring domains. While electrostatic mechanisms appear to dominate orientation toward each other of the redox partners to generate pre-collisional encounter complexes, hydrophobic forces are likely to foster electron transfer events by through-bonding or pi-stacking interactions. Moreover, electron-tunneling pathways seem to be operative as well. The availability of new P450 crystal structures together with improved validation strategies will undoubtedly permit the production of increasingly satisfactory three-dimensional donor-acceptor models serving to better understand the molecular principles governing functional association of the redox proteins.

Influence of lipophilicity on the interactions of N-alkyl-4-phenyl-1,2,3,6-tetrahydropyridines and their positively charged N-alkyl-4-phenylpyridinium metabolites with cytochrome P450 2D6

The relationship between lipophilicity and CYP2D6 affinity of cyclic tertiary (N-alkyl-4-phenyl-1,2,3,6-tetrahydropyridines) and quaternary (N-alkyl-4-phenylpyridinium) amines was examined. The 1,2,3,6-tetrahydropyridine scaffold was chosen due to its common occurrence in the structures of CYP2D6 ligands such as the Parkinsonian neurotoxin 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP) and the dehydrated haloperidol metabolite N-[4-(4-fluorophenyl)-4-oxobutyl]-4-(4-chlorophenyl)-1,2,3,6-tetrahydropyridine (HPTP). Likewise, the pyridinium framework is found in and 4-(4-chlorophenyl)-1-[4-(4-fluorophenyl)-4-oxobutyl]pyridinium and N-methyl-4-phenylpyridinium (MPP(+)), the positively charged metabolites of MPTP and haloperidol. The lack of CYP2D6 inhibition by MPTP and its pyridinium metabolite MPP(+) was due to their hydrophilic nature since higher N-alkyl homologs revealed substantial increases in inhibitory potency against recombinant CYP2D6-mediated bufuralol-1'-hydroxylation. The reasonable correlation between lipophilicity and CYP2D6 inhibition by pyridiniums and 1,2,3,6-tetrahydropyridines was only limited to straight chain N-alkyl analogs, since certain N-alkylaryl analogs of lower lipophilicity were better CYP2D6 inhibitors. CYP2D6 substrate properties of straight chain N-alkyltetrahydropyridines were also governed by lipophilicity, and N-heptyl-4-phenyl-1,2,3,6-tetrahydropyridine was the optimal substrate (K(mapp) = 0.63 microM). Metabolism studies indicated that the N-heptyl analog underwent monohydroxylation on the aromatic ring and on the N-heptyl group suggesting that 1,2,3,6-tetrahydropyridines can bind in more than one conformation in the CYP2D6 active site. Increased lipophilicity of haloperidol metabolites did not correlate with inhibitory potency since the more lipophilic HPTP metabolite was less potent as an inhibitor than reduced-haloperidol and reduced-HPTP. Furthermore, HPTP and reduced-HPTP, of comparable lipophilicity to the N-heptyltetrahydropyridine analog were inactive as CYP2D6 substrates. This observation suggests that steric constraints rather than lipophilicity are responsible for the lack of CYP2D6 substrate properties of cyclic tertiary amines tethered to bulky N-substituents. This phenomenon appears to be a common theme among several cyclic tertiary amine-containing anti-depressants and should be taken into consideration when designing central nervous system agents devoid of CYP2D6 substrate properties.

Role of conserved Asp293 of cytochrome P450 2C9 in substrate recognition and catalytic activity

Human cytochrome P450 2C9 (CYP2C9) is important in the metabolism of non-steroidal anti-inflammatory compounds such as diclofenac, the antidiabetic agent tolbutamide and other clinically important drugs, many of which are weakly acidic. Multiple sequence alignment of CYPs identified CYP2C9 Asp(293) as corresponding to Asp(301) of CYP2D6, which has been suggested to play a role in the binding of basic substrates to the latter enzyme. Replacement of Asp(293) with Ala (D293A) decreased activity by more than 90%, and led to an approx. 3- to 10-fold increase in K (m) values for the three test substrates tolbutamide, dextromethorphan and diclofenac. Conservative replacement of the carboxyl side chain in a Glu (D293E) mutant produced no significant changes in K (m) values and slight increases in k (cat) values. Changes in regiospecificity were observed for both the Ala and Glu substitutions; low levels of both dextromethorphan O- and N-demethylation were observed in the D293A mutant, whereas increased preference for O-demethylation was observed for the D293E mutant. Expression of constructs coding for Asn (D293N) and Gln (D293Q) substitutions failed to form a P450 correctly. Our analysis suggests a structural role for the carboxyl side chain of Asp(293) in CYP2C9 substrate binding and catalysis. The conservation of an Asp residue in other CYP families in a position equivalent to Asp(293) indicates a common mechanism for maintaining the active-site architecture.

Residues glutamate 216 and aspartate 301 are key determinants of substrate specificity and product regioselectivity in cytochrome P450 2D6

Cytochrome P450 2D6 (CYP2D6) metabolizes a wide range of therapeutic drugs. CYP2D6 substrates typically contain a basic nitrogen atom, and the active-site residue Asp-301 has been implicated in substrate recognition through electrostatic interactions. Our recent computational models point to a predominantly structural role for Asp-301 in loop positioning (Kirton, S. B., Kemp, C. A., Tomkinson, N. P., St.-Gallay, S., and Sutcliffe, M. J. (2002) Proteins 49, 216-231) and suggest a second acidic residue, Glu-216, as a key determinant in the binding of basic substrates. We have evaluated the role of Glu-216 in substrate recognition, along with Asp-301, by site-directed mutagenesis. Reversal of the Glu-216 charge to Lys or substitution with neutral residues (Gln, Phe, or Leu) greatly decreased the affinity (K(m) values increased 10-100-fold) for the classical basic nitrogen-containing substrates bufuralol and dextromethorphan. Altered binding was also manifested in significant differences in regiospecificity with respect to dextromethorphan, producing enzymes with no preference for N-demethylation versus O-demethylation (E216K and E216F). Neutralization of Asp-301 to Gln and Asn had similarly profound effects on substrate binding and regioselectivity. Intriguingly, removal of the negative charge from either 216 or 301 produced enzymes (E216A, E216K, and D301Q) with elevated levels (50-75-fold) of catalytic activity toward diclofenac, a carboxylate-containing CYP2C9 substrate that lacks a basic nitrogen atom. Activity was increased still further (>1000-fold) upon neutralization of both residues (E216Q/D301Q). The kinetic parameters for diclofenac (K(m) 108 microm, k(cat) 5 min(-1)) along with nifedipine (K(m) 28 microm, k(cat) 2 min(-1)) and tolbutamide (K(m) 315 microm, k(cat) 1 min(-1)), which are not normally substrates for CYP2D6, were within an order of magnitude of those observed with CYP3A4 or CYP2C9. Neutralizing both Glu-216 and Asp-301 thus effectively alters substrate recognition illustrating the central role of the negative charges provided by both residues in defining the specificity of CYP2D6 toward substrates containing a basic nitrogen.

Roles of human liver cytochrome P450 3A4 and 1A2 enzymes in the oxidation of myristicin

The aim of this work was to identify the form(s) of human liver cytochrome P450 (CYP) involved in the hepatic transformation of myristicin to its major metabolite, 5-allyl-1-methoxy-2,3-dihydroxybenzene. When microsomes prepared from different human liver samples were compared, the activity of 5-allyl-1-methoxy-2,3-dihydroxybenzene formation was well correlated (r(2)=0.87) with nifedipine oxidation (a marker of CYP3A4). With a microsomal sample having high CYP3A4 activity, microsomal oxidation of myristicin to the major metabolite (5-allyl-1-methoxy-2,3-dihydroxybenzene) was markedly inhibited by gestodene and ketoconazole, selective inhibitors of CYP3A enzymes, but not by any of several other P450 inhibitors. Antibodies raised against CYPs 3A4 and 1A2 could also inhibit the oxidation of myristicin, but antibodies recognizing other CYPs had no effect. The oxidation of myristicin to 5-allyl-1-methoxy-2,3-dihydroxybenzene was catalyzed by purified bacterial recombinant CYPs 3A4 and 1A2. These results provide evidence that CYP3A4 (and possibly other CYP3A enzymes) and CYP1A2 play roles in the formation of the major metabolite, 5-allyl-1-methoxy-2,3-dihydroxybenzene.

Homology modeling of rat and human cytochrome P450 2D (CYP2D) isoforms and computational rationalization of experimental ligand-binding specificities

The ligand-binding characteristics of rat and human CYP2D isoforms, i.e., rat CYP2D1-4 and human CYP2D6, were investigated by measuring IC(50) values of 11 known CYP2D6 ligands using 7-methoxy-4-(aminomethyl)coumarin (MAMC) as substrate. Like CYP2D6, all rat CYP2D isozymes catalyzed the O-demethylation of MAMC with K(m) and V(max) values ranging between 78 and 145 microM and 0.048 and 1.122 min(-1), respectively. To rationalize observed differences in the experimentally determined IC(50) values, homology models of the CYP2D isoforms were constructed. A homology model of CYP2D6 was generated on the basis of crystallized rabbit CYP2C5 and was validated on its ability to reproduce binding orientations corresponding to metabolic profiles of the substrates and to remain stable during unrestrained molecular dynamics simulations at 300 K. Twenty-two active site residues, sharing up to 59% sequence identity, were identified in the CYP2D binding pockets and included CYP2D6 residues Phe120, Glu216, and Asp301. Electrostatic potential calculations displayed large differences in the negative charge of the CYP2D active sites, which was consistent with observed differences in absolute IC(50) values. MD studies on the binding mode of sparteine, quinidine, and quinine in CYP2D2 and CYP2D6 furthermore concurred well with experimentally determined IC(50) values and metabolic profiles. The current study thus provides new insights into differences in the active site topology of the investigated CYP2D isoforms.

Impact of incorporating the 2C5 crystal structure into comparative models of cytochrome P450 2D6

Cytochrome P450 2D6 (CYP2D6) metabolizes approximately one third of the drugs in current clinical use. To gain insight into its structure and function, we have produced four different sets of comparative models of 2D6: one based on the structures of P450s from four different microorganisms (P450 terp, P450 eryF, P450 cam, and P450 BM3), another on the only mammalian P450 (2C5) structure available, and the other two based on alternative amino acid sequence alignments of 2D6 with all five of these structures. Principal component analysis suggests that inclusion of the 2C5 crystal structure has a profound effect on the modeling process, altering the general topology of the active site, and that the models produced differ significantly from all of the templates. The four models of 2D6 were also used in conjunction with molecular docking to produce complexes with the substrates codeine and 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP); this identified Glu 216 [in the F-helix; substrate recognition site (SRS) 2] as a key determinant in the binding of the basic moiety of the substrate. Our studies suggest that both Asp 301 and Glu 216 are required for metabolism of basic substrates. Furthermore, they suggest that Asp 301 (I-helix, SRS-4), a residue thought from mutagenesis studies to bind directly to the basic moiety of substrates, may play a key role in positioning the B'-C loop (SRS-1) and that the loss of activity on mutating Asp 301 may therefore be the result of an indirect effect (movement of the B'-C loop) on replacing this residue.

Rate-limiting steps in cytochrome P450 catalysis

Cytochrome P450 (P450) reactions are of interest because of their relevance to the oxidative metabolism of drugs, steroids, carcinogens, and other chemicals. One of the considerations about functional characterization is which steps of the catalytic cycle are rate-limiting. Detailed analysis indicates that several different steps can be rate-limiting with individual P450 reactions. N-Dealkylation of para-substituted N,N-dimethylanilines is a function of the electron withdrawing/donating properties of the substituent and the oxidation-reduction potential of the substrate, supporting a role in rate-limiting electron transfer from substrate to the high valent P450. In the oxidations of ethanol and acetaldehyde by human P450 2E1, a step following product formation must be the slow step (but not product release per se). Several oxidations catalyzed by human P450s 1A2 and 2D6 show slow C-H bond breaking, and apparent high-valent iron complexes accumulate in the reaction steady-state. Kinetic simulations were used to test the suitability of potential schemes and to probe the effects of changes in individual reaction steps.

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 mechanisms of substrate oxidation by cytochrome P450 2D6. Lack of an allosteric role of NADPH-cytochrome P450 reductase in catalytic regioselectivity

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. Differences in the regioselectivity of oxidation products formed in systems containing NADPH-P450 reductase/NADPH and the model oxidant cumene hydroperoxide have been proposed by others to be due to an allosteric influence of the reductase on P450 2D6 (Modi, S., Gilham, D. E., Sutcliffe, M. J., Lian, L.-Y., Primrose, W. U., Wolf, C. R., and Roberts, G. C. K. (1997) Biochemistry 36, 4461-4470). We examined the differences in the formation of oxidation products of N-methyl-4-phenyl-1,2,5,6-tetrahydropyridine, metoprolol, and bufuralol between reductase-, cumene hydroperoxide-, and iodosylbenzene-supported systems. Catalytic regioselectivity was not influenced by the presence of the reductase in any of the systems supported by model oxidants, ruling out allosteric influences. The presence of the reductase had little effect on the affinity of P450 2D6 for any of these three substrates. The addition of the reaction remnants of the model oxidants (cumyl alcohol and iodobenzene) to the reductase-supported system did not affect reaction patterns, arguing against steric influences of these products on catalytic regioselectivity. Label from H(2)18O was quantitatively incorporated into 1'-hydroxybufuralol in the iodosylbenzene- but not in the reductase- or cumene hydroperoxide-supported reactions. We conclude that the P450 systems utilizing NADPH-P450 reductase, cumene hydroperoxide, and iodosylbenzene use similar but distinct chemical mechanisms. These differences are the basis for the variable product distributions, not an allosteric influence of the reductase.