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

"Commandeuring" Xenobiotic Metabolism: Advances in Understanding Xenobiotic Metabolism

The understanding of how exogenous chemicals (xenobiotics) are metabolized, distributed, and eliminated is critical to determine the impact of the chemical and its metabolites to the (human) organism. This is part of the research and educational discipline ADMET (absorption, distribution, metabolism, elimination, and toxicity). Here, we review the work of Jan Commandeur and colleagues who have not only made a significant impact in understanding of phase I and phase II metabolism of several important compounds but also contributed greatly to the development of experimental techniques for the study of xenobiotic metabolism. Jan Commandeur’s work has covered a broad area of research, such as the development of online screening methodologies, the use of a combination of enzyme mutagenesis and molecular modeling for structure–activity relationship (SAR) studies, and the development of novel probe substrates. This work is the bedrock of current activities and brings the field closer to personalized (cohort-based) pharmacology, toxicology, and hazard/risk assessment.

In-Silico Modeling in Drug Metabolism and Interaction: Current Strategies of Lead Discovery

BACKGROUND

Drug metabolism is a complex mechanism of human body systems to detoxify foreign particles, chemicals, and drugs through bio alterations. It involves many biochemical reactions carried out by invivo enzyme systems present in the liver, kidney, intestine, lungs, and plasma. After drug administration, it crosses several biological membranes to reach into the target site for binding and produces the therapeutic response. After that, it may undergo detoxification and excretion to get rid of the biological systems. Most of the drugs and its metabolites are excreted through kidney via urination. Some drugs and their metabolites enter into intestinal mucosa and excrete through feces. Few of the drugs enter into hepatic circulation where they go into the intestinal tract. The drug leaves the liver via the bile duct and is excreted through feces. Therefore, the study of total methodology of drug biotransformation and interactions with various targets is costly.

METHODS

To minimize time and cost, in-silico algorithms have been utilized for lead-like drug discovery. Insilico modeling is the process where a computer model with a suitable algorithm is developed to perform a controlled experiment. It involves the combination of both in-vivo and in-vitro experimentation with virtual trials, eliminating the non-significant variables from a large number of variable parameters. Whereas, the major challenge for the experimenter is the selection and validation of the preferred model, as well as precise simulation in real physiological status.

RESULTS

The present review discussed the application of in-silico models to predict absorption, distribution, metabolism, and excretion (ADME) properties of drug molecules and also access the net rate of metabolism of a compound.

CONCLUSION

It helps with the identification of enzyme isoforms; which are likely to metabolize a compound, as well as the concentration dependence of metabolism and the identification of expected metabolites. In terms of drug-drug interactions (DDIs), models have been described for the inhibition of metabolism of one compound by another, and for the compound-dependent induction of drug-metabolizing enzymes.

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.

Ontogenesis of phase I hepatic drug metabolic enzymes in sheep

Cytochrome P450 (CYP) enzymes are important for the metabolism of many drugs. While there is information on their identity and ontogeny in humans and rodents, similar data in sheep are lacking. In the present study, cDNA sequences of several CYP enzymes (CYP2A6, CYP2C19, CYP2D6) were cloned by rapid amplification of cDNA ends. In adult, newborn and fetal sheep the mRNA and protein levels of these CYPs and the regulatory factor, hepatic nuclear factor 4α (HNF4α) were determined in liver samples using real-time PCR and western blotting. The effect of antenatal glucocorticoid on these enzymes was also studied by i.v. infusion of cortisol (0.45 mg h–1; 80 h) to another group of fetuses. The mRNA and protein levels of the CYPs and HNF4α were low or absent in the fetus, followed by increasing levels in the newborn and adult. Fetal cortisol administration significantly increased the mRNA and protein levels of CYP2D6. Moreover, the correlation observed between the CYP and HNF4α mRNA levels suggests a possible regulatory role for this transcription factor. The findings suggest that fetal and newborn lambs have a low ability to metabolise drugs that are substrates of these enzymes, and that this ability increases with advancing postnatal age, similar to the situation in humans.

IDSite: An accurate approach to predict P450-mediated drug metabolism

Accurate prediction of drug metabolism is crucial for drug design. Since a large majority of drugs metabolism involves P450 enzymes, we herein describe a computational approach, IDSite, to predict P450-mediated drug metabolism. To model induced-fit effects, IDSite samples the conformational space with flexible docking in Glide followed by two refinement stages using the Protein Local Optimization Program (PLOP). Sites of metabolism (SOMs) are predicted according to a physical-based score that evaluates the potential of atoms to react with the catalytic iron center. As a preliminary test, we present in this paper the prediction of hydroxylation and O-dealkylation sites mediated by CYP2D6 using two different models: a physical-based simulation model, and a modification of this model in which a small number of parameters are fit to a training set. Without fitting any parameters to experimental data, the Physical IDSite scoring recovers 83% of the experimental observations for 56 compounds with a very low false positive rate. With only 4 fitted parameters, the Fitted IDSite was trained with the subset of 36 compounds and successfully applied to the other 20 compounds, recovering 94% of the experimental observations with high sensitivity and specificity for both sets.

Using a homology model of cytochrome P450 2D6 to predict substrate site of metabolism

CYP2D6 is an important enzyme that is involved in first pass metabolism and is responsible for metabolizing ~25% of currently marketed drugs. A homology model of CYP2D6 was built using X-ray structures of ligand-bound CYP2C5 complexes as templates. This homology model was used in docking studies to rationalize and predict the site of metabolism of known CYP2D6 substrates. While the homology model was generally found to be in good agreement with the recently solved apo (ligand-free) X-ray structure of CYP2D6, significant differences between the structures were observed in the B' and F-G helical region. These structural differences are similar to those observed between ligand-free and ligand-bound structures of other CYPs and suggest that these conformational changes result from induced-fit adaptations upon ligand binding. By docking to the homology model using Glide, it was possible to identify the correct site of metabolism for a set of 16 CYP2D6 substrates 85% of the time when the 5 top scoring poses were examined. On the other hand, docking to the apo CYP2D6 X-ray structure led to a loss in accuracy in predicting the sites of metabolism for many of the CYP2D6 substrates considered in this study. These results demonstrate the importance of describing substrate-induced conformational changes that occur upon binding. The best results were obtained using Glide SP with van der Waals scaling set to 0.8 for both the receptor and ligand atoms. A discussion of putative binding modes that explain the distribution of metabolic sites for substrates, as well as a relationship between the number of metabolic sites and substrate size, are also presented. In addition, analysis of these binding modes enabled us to rationalize the typical hydroxylation and O-demethylation reactions catalyzed by CYP2D6 as well as the less common N-dealkylation.

A single amino acid residue, Ala 105, confers 16alpha-hydroxylase activity to human cytochrome P450 17alpha-hydroxylase/17,20 lyase

In adrenal steroidogenesis, CYP17 catalyses the 17alpha-hydroxylation of pregnenolone and progesterone and the subsequent 17,20-lyase reaction, yielding adrenal androgens. The enzyme exhibits distinctly different selectivities towards these substrates in various species. CYP17 has also been shown to exhibit 16alpha-hydroxylase activity towards progesterone in some species, with only human and chimp CYP17 catalysing the biosynthesis of substantial amounts of 16-OHprogesterone. The 16alpha-hydroxylase activity was investigated by introducing an Ala105Leu substitution into human CYP17. The converse mutation, Leu105Ala was introduced into the baboon, goat and pig enzymes. Wt human CYP17 converted approximately 30% progesterone to 16-OHprogesterone while the Ala105Leu mutant converted negligible amounts to 16-OHprogesterone ( approximately 9%), comparable to wt CYP17 of the other three species when expressed in COS-1 cells. The ratio of 17-hydroxylated products to 16-OHprogesterone of human CYP17 was 2.7 and that of the mutant human construct 10.5. Similar ratios were observed for human and goat CYP17 with the corresponding Ala or Leu residues. Although the Leu105Ala mutation of both baboon and pig CYP17 exhibited the same trend regarding the ratios, the rate of progesterone conversion was reduced. Coexpression with cytochrome b(5) significantly decreased the ratio of 17-hydroxylated products to 16-OHprogesterone in the Leu105 constructs, while effects were negligible with Ala at this position. Homology models show that Ala105 faces towards the active pocket in the predicted B'-C domain of CYP17. The smaller residue allows more flexibility of movement in the active pocket than Leu, presenting both the C16 and C17 of progesterone to the iron-oxy complex.

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.

Mechanism-based inactivation of CYP2B1 and its F-helix mutant by two tert-butyl acetylenic compounds: covalent modification of prosthetic heme versus apoprotein

The mechanism-based inactivation of cytochrome CYP2B1 [wild type (WT)] and its Thr205 to Ala mutant (T205A) by tert-butylphenylacetylene (BPA) and tert-butyl 1-methyl-2-propynyl ether (BMP) in the reconstituted system was investigated. The inactivation of WT by BPA exhibited a k(inact)/K(I) value of 1343 min(-1)mM(-1) and a partition ratio of 1. The inactivation of WT by BMP exhibited a k(inact)/K(I) value of 33 min(-1)mM(-1) and a partition ratio of 10. Liquid chromatography/tandem mass spectrometry analysis (LC/MS/MS) of the WT revealed 1) inactivation by BPA resulted in the formation of a protein adduct with a mass increase equivalent to the mass of BPA plus one oxygen atom, and 2) inactivation by BMP resulted in the formation of multiple heme adducts that all exhibited a mass increase equivalent to BMP plus one oxygen atom. LC/MS/MS analysis indicated the formation of glutathione (GSH) conjugates by the reaction of GSH with the ethynyl moiety of BMP or BPA with the oxygen being added to the internal or terminal carbon. For the inactivation of T205A by BPA and BMP, the k(inact)/K(I) values were suppressed by 100- and 4-fold, respectively, and the partition ratios were increased 9- and 3.5-fold, respectively. Only one major heme adduct was detected following the inactivation of the T205A by BMP. These results show that the Thr205 in the F-helix plays an important role in the efficiency of the mechanism-based inactivation of CYP2B1 by BPA and BMP. Homology modeling and substrate docking studies were presented to facilitate the interpretation of the experimental results.

Understanding CYP2D6 interactions

Owing to the polymorphic nature of CYP2D6, clinically significant issues can arise when drugs rely on that enzyme either for clearance, or metabolism to an active metabolite. Available screening methods to determine if the compound is likely to cause drug-drug interactions, or is likely to be a victim of inhibition of CYP2D6 by other compounds will be described. Computational models and examples will be given on strategies to design out the CYP2D6 liabilities for both heme-binding compounds and non-heme-binding compounds.

AmineDB: Large scale docking of amines with CYP2D6 and scoring for druglike properties—towards defining the scope of the chemical defense against foreign amines in humans

Organic amines are prevalent in nature and in drugs, especially the psychotherapeutic agents, and a major defense against potentially toxic amines is metabolism by CYP2D6. In order to understand better the constraints on the broad specificity of CYP2D6, 4207 amines were docked into the binding site of this enzyme. Docking poses were found predominantly with the positively charged amino groups closest to Asp301, with aromatic rings close to Phe120 and sometimes extending as far as Phe483. Organic amines that bind best to CYP2D6 tend to have larger molecular weights and logP values. Organic amines that score highly as being druglike, based on a Bayesian model constructed using a 5223-drug training set, are least likely to bind to CYP2D6. This correlation suggests that the set of known drugs, which have been largely designed or selected to avoid high affinity CYP binding, partially encodes the binding site preferences (or rather anti-preferences) of CYP2D6. Finally, in order to benchmark our docking and druglike scoring procedures, an analysis of psychotherapeutic agents is presented. All of these data, including the 4207 AM1-optimized ligand structures in proper ionization states, docking poses and scores, Druglike Scores and Lipinski properties, can be viewed from an online database, the AmineDB.

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.

State-of-the-art tools for computational site of metabolism predictions: comparative analysis, mechanistical insights, and future applications

In drug design, it is crucial to have reliable information on how a chemical entity behaves in the presence of metabolizing enzymes. This requires substantial experimental efforts. Consequently, being able to predict the likely site/s of metabolism in any compound, synthesized or virtual, would be highly beneficial and time efficient. In this work, six different methodologies for predictions of the site of metabolism (SOM) have been compared and validated using structurally diverse data sets of drug-like molecules with well-established metabolic pattern in CYP3A4, CYP2C9, or both. Three of the methods predict the SOM based on the ligand's chemical structure, two additional methods use structural information of the enzymes, and the sixth method combines structure and ligand similarity and reactivity. The SOM is correctly predicted in 50 to 90% of the cases, depending on method and enzyme, which is an encouraging rate. We also discuss the underlying mechanisms of cytochrome P450 metabolism in the light of the results from this comparison.

Evidence for the functional role of residues in the B'-C loop of baboon cytochrome P450 side-chain cleavage (CYP11A1) obtained by site-directed mutagenesis, kinetic analysis and homology modelling

To gain further insight into the structure/function relationship of cytochrome P450 side-chain cleavage (CYP11A1), this enzyme was investigated in the Cape baboon (Papio ursinus). Four constructs were cloned and characterised in non-steroidogenic mammalian COS-1 cells. Wild type recombinant baboon CYP11A1 cDNA yielded a K(m) value of 1.6 microM for 25-hydroxycholesterol. The single amino acid substitutions, I98Q and I98K resulted in a 1.7- and 2.8-fold increases in K(m) values, respectively. Conversely, the introduction of the mutation, K103A, resulted in a 1.8-fold decrease in K(m). A homology model of CYP11A1, based on the crystal structures of CYP102 and CYP2C5, revealed that residues 98 and 103 lie within the B'-C loop and contribute to the spatial orientation and structural integrity of this domain. Based on these results we propose a topological model of the CYP11A1 active pocket, which is supported by substrate docking analysis and kinetic studies.

Catalytic site prediction and virtual screening of cytochrome P450 2D6 substrates by consideration of water and rescoring in automated docking

Automated docking strategies successfully applied to binding mode predictions of ligands in Cyt P450 crystal structures in an earlier study (de Graaf et al. J. Med. Chem. 2005, 7, 2308-2318) were used for the catalytic site prediction (CSP) of 65 substrates in a CYP2D6 homology model. The consideration of water molecules at predicted positions in the active site and the rescoring of pooled docking poses from four different docking programs (AutoDock, FlexX, GOLD-Goldscore, and GOLD-Chemscore) with the SCORE scoring function enabled the successful prediction of experimentally reported sites of catalysis of more than 80% of the substrates. Three docking algorithms (FlexX, GOLD-Goldscore, and GOLD-Chemscore) were subsequently used in combination with six scoring functions (Chemscore, DOCK, FlexX, GOLD, PMF, and SCORE) to assess the ability of docking-based virtual screening methods to prioritize known CYP2D6 substrates seeded into a drug-like chemical database (in the absence and presence of active-site water molecules). Finally, the optimal docking strategy in terms of virtual screening accuracy, GOLD-Chemscore with the consideration of active-site water (60% of known substrates recovered in the top 5% of the ranked drug-like database), was verified experimentally; it was successfully used to identify high-affinity CYP2D6 ligands among a larger proprietary database.

IN SILICO PREDICTION OF DRUG BINDING TO CYP2D6: IDENTIFICATION OF A NEW METABOLITE OF METOCLOPRAMIDE

Patients with cancer often take many different classes of drugs to treat the effects of their malignancy and the side effects of treatment, as well as their comorbidities. The potential for drug-drug interactions that may affect the efficacy of anticancer treatment is high, and a major source of such interactions is competition for the drug-metabolizing enzymes, cytochromes P450 (P450s). We have examined a series of 20 drugs commonly prescribed to cancer patients to look for potential interactions via CYP2D6. We used a homology model of CYP2D6, together with molecular docking techniques, to perform an in silico screen for binding to CYP2D6. Experimental IC50 values were determined for these compounds and compared with the model predictions to reveal a correlation with a regression coefficient of r2= 0.61. Importantly, the docked conformation of the commonly prescribed antiemetic metoclopramide predicted a new site of metabolism that was further investigated through in vitro analysis with recombinant CYP2D6. An aromatic N-hydroxy metabolite of metoclopramide, consistent with predictions from our modeling studies, was identified by high-performance liquid chromatography/mass spectrometry. This metabolite was found to represent a major product of metabolism in human liver microsomes, and CYP2D6 was identified as the main P450 isoform responsible for catalyzing its formation. In view of the prevalence of interindividual variation in the CYP2D6 genotype and phenotype, we suggest that those experiencing adverse reactions with metoclopramide, e.g., extrapyramidal syndrome, are likely to have a particular CYP2D6 genotype/phenotype. This warrants further investigation.

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.

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.

COMPARISON OF METHODS FOR THE PREDICTION OF THE METABOLIC SITES FOR CYP3A4-MEDIATED METABOLIC REACTIONS

Predictions of the metabolic sites for new chemical entities, synthesized or only virtual, are important in the early phase of drug discovery to guide chemistry efforts in the synthesis of new compounds with reduced metabolic liability. This information can now be obtained from in silico predictions, and therefore, a thorough and unbiased evaluation of the computational techniques available is needed. Several computational methods to predict the metabolic hot spots are emerging. In this study, metabolite identification using MetaSite and a docking methodology, GLUE, were compared. Moreover, the published CYP3A4 crystal structure and computed CYP3A4 homology models were compared for their usefulness in predicting metabolic sites. A total of 227 known CYP3A4 substrates reported to have one or more metabolites adding up to 325 metabolic pathways were analyzed. Distance-based fingerprints and four-point pharmacophore derived from GRID molecular interaction fields were used to characterize the substrate and protein in MetaSite and the docking methodology, respectively. The CYP3A4 crystal structure and homology model with the reactivity factor enabled achieved a similar prediction success (78%) using the MetaSite method. The docking method had a relatively lower prediction success (approximately 57% for the homology model), although it still may provide useful insights for interactions between ligand and protein, especially for uncommon reactions. The MetaSite methodology is automated, rapid, and has relatively accurate predictions compared with the docking methodology used in this study.

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.

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.

IN SILICO AND IN VITRO SCREENING FOR INHIBITION OF CYTOCHROME P450 CYP3A4 BY COMEDICATIONS COMMONLY USED BY PATIENTS WITH CANCER

Cytochrome P450 3A4 (CYP3A4) is the major enzyme responsible for phase I drug metabolism of many anticancer agents. It is also a major route for metabolism of many drugs used by patients to treat the symptoms caused by cancer and its treatment as well as their other illnesses, for example, cardiovascular disease. To assess the ability to inhibit CYP3A4 of drugs most commonly used by our patients during cancer therapy, we have made in silico predictions based on the crystal structures of CYP3A4. From this set of 33 common comedicated drugs, 10 were predicted to be inhibitors of CYP3A4, with the antidiarrheal drug loperamide predicted to be the most potent. There was significant correlation (r(2) = 0.75-0.66) between predicted affinity and our measured IC(50) values, and loperamide was confirmed as a potent inhibitor (IC(50) of 0.050 +/- 0.006 microM). Active site docking studies predicted an orientation of loperamide consistent with formation of the major (N-demethylated) metabolite, where it interacts with the phenylalanine cluster and Arg-212 and Glu-374; experimental evidence for the latter interaction comes from the approximately 12-fold increase in K(M) for loperamide observed for the Glu-374-Gln mutant. The commonly prescribed drugs loperamide, amitriptyline, diltiazem, domperidone, lansoprazole, omeprazole, and simvastatin were identified by our in silico and in vitro screens as relatively potent inhibitors of CYP3A4 that have the potential to interact with cytotoxic agents to cause adverse effects, highlighting the likelihood of drug-drug interactions affecting chemotherapy treatment.

Chemical proteomic tool for ligand mapping of CYP antitargets: an NMR-compatible 3D QSAR descriptor in the Heme-Based Coordinate System

Chemical proteomic strategies strive to probe and understand protein-ligand interactions across gene families. One gene family of particular interest in drug and xenobiotic metabolism are the cytochromes P450 (CYPs), the topic of this article. Although numerous tools exist to probe affinity of CYP-ligand interactions, fewer exist for the rapid experimental characterization of the structural nature of these interactions. As a complement to recent advances in X-ray crystallography, NMR methods are being developed that allow for fairly high throughput characterization of protein-ligand interactions. One especially promising NMR approach involves the use of paramagnetic induced relaxation effects to measure distances of ligand atoms from the heme iron in CYP enzymes. Distances obtained from these T(1) relaxation measurements can be used as a direct source of 1-dimensional structural information or to restrain a ligand docking to generate a 3-dimensional data set. To facilitate such studies, we introduce the concept of the Heme-Based Coordinate System and present how it can be used in combination with NMR T(1) relaxation data to derive 3D QSAR descriptors directly or in combination with in silico docking. These descriptors should have application in defining the binding preferences of CYP binding sites using 3D QSAR models. They are especially well-suited for the biasing of fragment assembly and combinatorial chemistry drug design strategies, to avoid fragment or reagent combinations with enhanced affinity for CYP antitargets.

Why Is Quinidine an Inhibitor of Cytochrome P450 2D6?

We have previously shown that Phe(120), Glu(216), and Asp(301) in the active site of cytochrome P450 2D6 (CYP2D6) play a key role in substrate recognition by this important drug-metabolizing enzyme (Paine, M. J., McLaughlin, L. A., Flanagan, J. U., Kemp, C. A., Sutcliffe, M. J., Roberts, G. C., and Wolf, C. R. (2003) J. Biol. Chem. 278, 4021-4027 and Flanagan, J. U., Maréchal, J.-D., Ward, R., Kemp, C. A., McLaughlin, L. A., Sutcliffe, M. J., Roberts, G. C., Paine, M. J., and Wolf, C. R. (2004) Biochem. J. 380, 353-360). We have now examined the effect of mutations of these residues on interactions of the enzyme with the prototypical CYP2D6 inhibitor, quinidine. Abolition of the negative charge at either or both residues 216 and 301 decreased quinidine inhibition of bufuralol 1'-hydroxylation and dextromethorphan O-demethylation by at least 100-fold. The apparent dissociation constants (K(d)) for quinidine binding to the wild-type enzyme and the E216D and D301E mutants were 0.25-0.50 microm. The amide substitution of Glu(216) or Asp(301) resulted in 30-64-fold increases in the K(d) for quinidine. The double mutant E216Q/D301Q showed the largest decrease in quinidine affinity, with a K(d) of 65 microm. Alanine substitution of Phe(120), Phe(481),or Phe(483) had only a minor effect on the inhibition of bufuralol 1'-hydroxylation and dextromethorphan O-demethylation and on binding. In contrast to the wild-type enzyme, a number of the mutants studied were found to be able to metabolize quinidine. E216F produced O-demethylated quinidine, and F120A and E216Q/D301Q produced both O-demethylated quinidine and 3-hydroxyquinidine metabolites. Homology modeling and molecular docking were used to predict the modes of quinidine binding to the wild-type and mutant enzymes; these were able to rationalize the experimental observations.

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.

Synthesis and evaluation of imidazolylmethylenetetrahydronaphthalenes and imidazolylmethyleneindanes: potent inhibitors of aldosterone synthase

Elevated plasma aldosterone levels play a detrimental role in certain forms of congestive heart failure and myocardial fibrosis. We proposed aldosterone synthase (CYP11B2) as a novel target for the treatment of these diseases. In this study, the synthesis and biological evaluation of substituted E- and Z-imidazolylmethylenetetrahydronaphthalenes and E- and Z-imidazolylmethyleneindanes (compounds 1a,b-9a,b) is described. The compounds were prepared by a Wittig-like reaction. They were tested for activity using bovine CYP11B and human CYP11B2 expressed in fission yeast and V79 MZh cells. Selectivity was determined toward human CYP11B1, CYP19, and CYP17. Especially in the case of CYP11B1 (steroid 11beta-hydroxylase), selectivity is a crucial issue, since sequence homology between this enzyme and the target enzyme is very high (93%). On the basis of the X-ray structure of human CYP2C9, a protein model of CYP11B2 was developed and docking experiments with the title compounds were performed. The biological results revealed highly potent inhibitors of CYP11B2 (IC(50) = 4-93 nM). The Z-isomers usually were more active than the corresponding E-isomers. Different inhibitory profiles could be observed: rather selective inhibitors of CYP11B1, dual inhibitors of both enzymes, and rather selective inhibitors of CYP11B2. The chloro derivative 8b was found to be a highly potent CYP11B2 inhibitor (IC(50) = 4 nM) showing a 5-fold selectivity for CYP11B1 (IC(50) = 20 nM). This compound could be an interesting lead for further optimization as a therapeutic agent. It also could be used as well as the CYP11B1 selective compounds as a pharmacological tool.

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.

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.

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.

Validation of Model of Cytochrome P450 2D6: An in Silico Tool for Predicting Metabolism and Inhibition

There has been much interest in the development of a predictive model of cytochrome P450 2D6 particularly because this enzyme is involved in the oxidation of at least 50 drugs. Previously we have described the combined use of homology modeling and molecular docking to correctly position a range of substrates in the CYP2D6 active site with the known sites of metabolism above the heme. Here, our approach identifies correctly the site of metabolism of the atypical (no basic nitrogen) cytochrome P450 2D6 substrate, spirosulfonamide. The same method is used to screen a small compound database for cytochrome P450 2D6 inhibition. A database containing 33 compounds from the National Cancer Institute database was docked into our cytochrome P450 2D6 homology model using the program GOLDv2.0. Experimental IC50 values for the 33 compounds were determined; comparison with the corresponding docked scores revealed a correlation with a regression coefficient of r2 = 0.61 (q2 = 0.59). The method was able to discriminate between tight and weak binding compounds and correctly identified several novel inhibitors. The results therefore suggest that our approach, which combines homology modeling with molecular docking, has produced a useful predictive in silico tool for cytochrome P450 2D6 inhibition, which is best used as one filter in a multifilter database screen.

The functional role of threonine-205 in the mechanism-based inactivation of P450 2B1 by two ethynyl substrates: the importance of the F helix in catalysis

We have previously demonstrated that substituting Val for Thr-205 in P450 2B1 abolishes the 16beta-hydroxylation of testosterone and markedly decreases the ability of 2-ethnylnaphthalene (2EN) and 17alpha-ethynylestradiol (17EE) to inactivate P450 2B1. The role of Thr-205 has been further investigated by measuring the kinetics of the mechanism-based inactivation of the 7-ethoxy-(trifluoromethyl)coumarin deethylation activity of 2B1 by 2EN and 17EE in wild-type (WT) and mutant P450s. In general, the kinetics of the inactivation of the Ser and Ala mutants was not significantly altered compared with WT. In contrast, the efficiency of the inactivation of the Val mutant decreased by approximately 6- and approximately 30-fold for 2EN and 17EE, respectively. High-pressure liquid chromatography (HPLC) analysis and SDS gel electrophoresis demonstrated the covalent binding of radiolabeled 2EN- and 17EE-reactive intermediates to the WT apoprotein, but not the Val mutant. The Val mutant was able to metabolize 2EN to 2-naphthylacetic acid, except the initial rate was slower than the WT. HPLC analysis of the 17EE incubation mixtures revealed three major metabolites and showed a correlation between the efficiency of inactivation and the generation of one of the major metabolites (C). Metabolite C was generated by the WT, Ser mutant, and Ala mutant. Metabolite C may be formed by the oxidation of the ethynyl group, and this reactive intermediate contributes to the inactivation of P450 2B1 by 17EE. The site-specific mutation of one residue, Thr-205 to Val, is sufficient to alter the profile of products formed during 17EE metabolism, such that very low levels of metabolite C are formed and inactivation is essentially abolished.

Differential Roles of Arg97, Asp293, and Arg108 in Enzyme Stability and Substrate Specificity of CYP2C9

CYP2C9 metabolizes a wide range of drugs, many of which are negatively charged at physiological pH. Therefore, it has been thought that complementarily charged amino acid(s) are critically involved in substrate binding. Previous studies have implicated arginine residues at positions 97, 105, and 108 and aspartate at position 293 in the normal catalytic function of the enzyme. To elucidate the role of these amino acids in the substrate specificity of CYP2C9, a series of mutants were constructed and analyzed for functional activity, thermal stability, and ligand binding. Charge-modifying mutations at positions 97, 105, and 293 decreased catalytic activity toward diclofenac, (S)-warfarin, and pyrene in a substrate-independent manner with Arg105 the least, and Arg97 the most, sensitive amino acids in this regard. Decreases in functional activity paralleled thermal instability of the mutants, suggesting that loss of function reflects more generalized structural changes rather than the absence of a specifically charged amino acid at these three positions. The R108H mutant was inactive toward all three substrates because of unexpected nitrogen ligation to the heme. Conversely, the R108F mutant exhibited substrate-dependent catalytic behavior, with almost complete loss of activity toward (S)-warfarin and diclofenac, but preservation of pyrene metabolism. In addition, the R108F mutation abrogated the Type I difference spectra induced by flurbiprofen and benzbromarone, obligate anions at physiological pH. These data identify critical roles for Arg97 and Asp293 in the structural stability of the enzyme and demonstrate a selective role for Arg108 in the binding and metabolism of negatively charged substrates of CYP2C9.

TOPOLOGICAL CHANGES IN THE CYP3A4 ACTIVE SITE PROBED WITH PHENYLDIAZENE: EFFECT OF INTERACTION WITH NADPH-CYTOCHROME P450 REDUCTASE AND CYTOCHROMEB5AND OF SITE-DIRECTED MUTAGENESIS

The active site topology of heterologously expressed CYP3A4 purified from an Escherichia coli expression system was examined using phenyldiazene. Incubation of CYP3A4 with phenyldiazene and subsequent oxidation yielded all four potential N-phenylprotoporphyrin IX regioisomers derived from attack on an available nitrogen atom in pyrrole rings B, A, C, or D (N(B):N(A):N(C):N(D) = 6:73:7: 13). Further study using 28 active site mutants showed that substitution of residues closer to the heme, Ala-305, Thr-309, or Ala-370, with a larger residue caused the most drastic changes in regioisomer formation, which reflected the location of each amino acid residue replaced in a CYP3A4 homology model. Previous studies have suggested a conformational change in CYP3A4 upon binding of NADPH-cytochrome P450 reductase (CPR) or cytochrome b(5) (b(5)). Therefore, regioisomer formation was also compared in the absence of redox partners and in the presence of CPR, b(5), or both. Formation of all four regioisomers in CYP3A4 wild type, particularly the minor ones, was reduced in the presence of b(5). CPR also greatly decreased the three minor isomers but increased the major isomer significantly. The presence of b(5) and CPR restored minor isomer formation and suppressed the enhancement of N(A) formation caused by CPR alone. Interestingly, the effects of the redox partners differed among representative active site mutants. In particular, the increase in N(C) upon substitution of Ala-370 with Phe was significantly reversed in the presence of redox partners, strongly suggesting that a conformational change occurs around pyrrole ring C due to protein-protein interactions between CYP3A4 and CPR or b(5).

Contribution of individual cytochrome P450 isozymes to the O-demethylation of the psychotropic beta-carboline alkaloids harmaline and harmine

The psychotropic beta-carboline alkaloids, showing high affinity for 5-hydroxytryptamine, dopamine, benzodiazepine, and imidazoline receptors and the stimulation of locus coeruleus neurons, are formed endogenously from tryptophan-derived indolealkylamines through the Pictet-Spengler condensation with aldehydes in both plants and mammals. Cytochromes P450 1A1 (18.5), 1A2 (20), and 2D6 (100) catalyzed the O-demethylation of harmaline, and CYP1A1 (98.5), CYP1A2 (35), CYP2C9 (16), CYP2C19 (30), and CYP2D6 (115) catalyzed that of harmine (relative activities). The dehydrogenation/aromatization of harmaline to harmine was not carried out by aromatase (CYP19), CYP1A2, CYP2C9, CYP2D6, CYP3A4, pooled recombinant cytochromes P450, or human liver microsomes (HLMs). Kinetic parameters were calculated for the O-demethylations mediated by each isozyme and by pooled HLMs. K(cat) (min(-1)) and K(m) Awake M) values for harmaline were: CYP1A1, 10.8 and 11.8; CYP1A2, 12.3 and 13.3; CYP2C9, 5.3 and 175; CYP2C19, 10.3 and 160; and CYP2D6, 39.9 and 1.4. Values for harmine were: CYP1A1, 45.2 and 52.2; CYP1A2, 9.2 and 14.7; CYP2C9, 11.9 and 117; CYP2C19, 21.4 and 121; and CYP2D6, 29.7 and 7.4. Inhibition studies using monoclonal antibodies confirmed that CYP1A2 and CYP2D6 were the major isozymes contributing to both harmaline (20% and 50%, respectively) and harmine (20% and 30%) O-demethylations in pooled HLMs. The turnover numbers for CYP2D6 are among the highest ever reported for a CYP2D6 substrate. Finally, CYP2D6-transgenic mice were found to have increased harmaline and harmine O-demethylase activities as compared with wild-type mice. These findings suggest a role for polymorphic CYP2D6 in the pharmacology and toxicology of harmine and harmaline.

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.