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

Two highly selected mutations in the tandemly duplicated CYP6P4a and CYP6P4b genes drive pyrethroid resistance in Anopheles funestus in West Africa

BACKGROUND

Gaining a comprehensive understanding of the genetic mechanisms underlying insecticide resistance in malaria vectors is crucial for optimising the effectiveness of insecticide-based vector control methods and developing diagnostic tools for resistance management. Considering the heterogeneity of metabolic resistance in major malaria vectors, the implementation of tailored resistance management strategies is essential for successful vector control. Here, we provide evidence demonstrating that two highly selected mutations in CYP6P4a and CYP6P4b are driving pyrethroid insecticide resistance in the major malaria vector Anopheles funestus, in West Africa.

RESULTS

Continent-wide polymorphism survey revealed escalated signatures of directional selection of both genes between 2014 and 2021. In vitro insecticide metabolism assays with recombinant enzymes from both genes showed that mutant alleles under selection exhibit higher metabolic efficiency than their wild-type counterparts. Using the GAL4-UAS expression system, transgenic Drosophila flies overexpressing mutant alleles exhibited increased resistance to pyrethroids. These findings were consistent with in silico predictions which highlighted changes in enzyme active site architecture that enhance the affinity of mutant alleles for type I and II pyrethroids. Furthermore, we designed two DNA-based assays for the detection of CYP6P4a-M220I and CYP6P4b-D284E mutations, showing their current confinement to West Africa. Genotype/phenotype correlation analyses revealed that these markers are strongly associated with resistance to types I and II pyrethroids and combine to drastically reduce killing effects of pyrethroid bed nets.

CONCLUSIONS

Overall, this study demonstrated that CYP6P4a and CYP6P4b contribute to pyrethroid resistance in An. funestus and provided two additional insecticide resistance molecular diagnostic tools that would contribute to monitoring and better management of resistance.

How does multiple substrate binding lead to substrate inhibition of CYP2D6 metabolizing dextromethorphan? A theoretical study

CYP2D6 is one of the most important metalloenzymes involved in the biodegradation of many drug molecules in the human body. It has been found that multiple substrate binding can lead to substrate inhibition of CYP2D6 metabolizing dextromethorphan (DM), but the corresponding theoretical mechanism is rarely reported. Therefore, we chose DM as the probe and performed molecular dynamics simulations and quantum mechanical calculations on CYP2D6-DM systems to investigate the mechanism of how the multiple substrate binding leads to the substrate inhibition of CYP2D6 metabolizing substrates. According to our results, three gate residues (Arg221, Val374, and Phe483) for the catalytic pocket are determined. We also found that the multiple substrate binding can lead to substrate inhibition by reducing the stability of CYP2D6 binding DM and increasing the reactive activation energy of the rate-determining step. Our findings would help to understand the substrate inhibition of CYP2D6 metabolizing the DM and enrich the knowledge of the drug-drug interactions for the cytochrome P450 superfamily.

Sequential phase I metabolism of pyrethroids by duplicated CYP6P9 variants results in the loss of the terminal benzene moiety and determines resistance in the malaria mosquito Anopheles funestus

Pyrethroid resistance in Anopheles funestus is threatening the eradication of malaria. One of the major drivers of pyrethroid resistance in An. funestus are cytochrome P450 monooxygenases CYP6P9a and CYP6P9b, which are found upregulated in resistant An. funestus populations from Sub-Saharan Africa and are known to metabolise pyrethroids. Here, we have functionally expressed CYP6P9a and CYP6P9b variants and investigated their interactions with azole-fungicides and pyrethroids. Some azole fungicides such as prochloraz inhibited CYP6P9a and CYP6P9b at nanomolar concentrations, whereas pyrethroids were weak inhibitors (>100 μM). Amino acid sequence comparisons suggested that a valine to isoleucine substitution at position 310 in the active site cavity of CYP6P9a and CYP6P9b, respectively, might affect substrate binding and metabolism. We therefore swapped the residues by site directed mutagenesis to produce CYP6P9a^I310V and CYP6P9b^V310I. CYP6P9b^V310I produced stronger metabolic activity towards coumarin substrates and pyrethroids, particularly permethrin. The V310I mutation was previously also detected in a pyrethroid resistant field population of An. funestus in Benin. Additionally, we found the first metabolite of permethrin and deltamethrin after hydroxylation, 4'OH permethrin and 4'OH deltamethrin, were also suitable substrates for CYP6P9-variants, and were depleted by both enzymes to a higher extent than as their respective parent compounds (approximately 20% more active). Further, we found that both metabolites were toxic against An. funestus FANG (pyrethroid susceptible) but not towards FUMOZ-R (pyrethroid resistant) mosquitoes, the latter suggesting detoxification by overexpressed CYP6P9a and CYP6P9b. We confirmed by mass-spectrometric analysis that CYP6P9a and CYP6P9b are capable of cleaving phenoxybenzyl-ethers in type I pyrethroid permethrin and type II pyrethroid deltamethrin and that both enzymes preferentially metabolise trans-permethrin. This provides new insight into the metabolism of pyrethroids and a greater understanding of the molecular mechanisms of pyrethroid resistance in An. funestus.

The regioselectivity of the interaction between dextromethorphan and CYP2D6

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

In vitro inhibitory effects of glucosamine, chondroitin and diacerein on human hepatic CYP2D6

OBJECTIVES

Glucosamine, chondroitin and diacerein are natural compounds commonly used in treating osteoarthritis. Their concomitant intake may trigger drug-natural product interactions. Cytochrome P450 (CYP) has been implicated in such interactions. Cytochrome P450 2D6 (CYP2D6) is a major hepatic CYP involved in metabolism of 25% of the clinical drugs. This study aimed to investigate the inhibitory effect of these antiarthritic compounds on CYP2D6.

METHODS

CYP2D6 was heterologously expressed in Escherichia coli. CYP2D6-antiarthritic compound interactions were studied using in vitro enzyme kinetics assay and molecular docking.

RESULTS

The high-performance liquid chromatography (HPLC)-based dextromethorphan O-demethylase assay was established as CYP2D6 marker. All glucosamines and chondroitins weakly inhibited CYP2D6 (IC₅₀ values >300 µM). Diacerein exhibited moderate inhibition with IC₅₀ and K _i values of 34.99 and 38.27 µM, respectively. Its major metabolite, rhein displayed stronger inhibition potencies (IC₅₀=26.22 μM and K _i =32.27 μM). Both compounds exhibited mixed-mode of inhibition. In silico molecular dockings further supported data from the in vitro study. From in vitro-in vivo extrapolation, rhein presented an area under the plasma concentration-time curve (AUC) ratio of 1.5, indicating low potential to cause in vivo inhibition.

CONCLUSIONS

Glucosamine, chondroitin and diacerein unlikely cause clinical interaction with the drug substrates of CYP2D6. Rhein, exhibits only low potential to cause in vivo inhibition.

In vitro inhibitory effects of glucosamine, chondroitin and diacerein on human hepatic CYP2D6

OBJECTIVES

Glucosamine, chondroitin and diacerein are natural compounds commonly used in treating osteoarthritis. Their concomitant intake may trigger drug-natural product interactions. Cytochrome P450 (CYP) has been implicated in such interactions. Cytochrome P450 2D6 (CYP2D6) is a major hepatic CYP involved in metabolism of 25% of the clinical drugs. This study aimed to investigate the inhibitory effect of these antiarthritic compounds on CYP2D6.

METHODS

CYP2D6 was heterologously expressed in Escherichia coli. CYP2D6-antiarthritic compound interactions were studied using in vitro enzyme kinetics assay and molecular docking.

RESULTS

The high-performance liquid chromatography (HPLC)-based dextromethorphan O-demethylase assay was established as CYP2D6 marker. All glucosamines and chondroitins weakly inhibited CYP2D6 (IC50 values >300 µM). Diacerein exhibited moderate inhibition with IC50 and K i values of 34.99 and 38.27 µM, respectively. Its major metabolite, rhein displayed stronger inhibition potencies (IC50=26.22 μM and K i =32.27 μM). Both compounds exhibited mixed-mode of inhibition. In silico molecular dockings further supported data from the in vitro study. From in vitro-in vivo extrapolation, rhein presented an area under the plasma concentration-time curve (AUC) ratio of 1.5, indicating low potential to cause in vivo inhibition.

CONCLUSIONS

Glucosamine, chondroitin and diacerein unlikely cause clinical interaction with the drug substrates of CYP2D6. Rhein, exhibits only low potential to cause in vivo inhibition.

Tryptophan-75 Is a Low-Energy Channel-Gating Residue that Facilitates Substrate Egress/Access in Cytochrome P450 2D6

CYP2D6 is a major drug metabolizing enzyme with a buried active site. Channels leading to the active site from various enzyme surfaces are believed to facilitate ligand egress and access to the active site. The present study used molecular dynamics (MD) and in vitro studies with CYP2D6*1 and a Trp75-to-Ala mutant to examine channel gating in CYP2D6 by Trp75. MD simulations measured energy landscapes of Trp75 conformations and simulated substrate passage within channel 2b using bufuralol as a model substrate. Trp75 alternated between multiple stable states that supported substrate transport along channel 2b with low-energy barriers between states (∼ -1 kcal/mol). Trp75 conformations were stabilized primarily by hydrogen bonding between Trp75 and Glu222, Asn226, Ala225, or Gln72. Energy barriers were low between Trp75 conformations, allowing Trp75 to easily move between various conformations over time and to function in both binding to and moving substrates in the 2b channel of CYP2D6. Michaelis-Menten kinetic studies completed with purified enzyme in a reconstituted system showed overall reduced enzyme efficiency for metabolism of bufuralol and dextromethorphan by the Trp75Ala mutant compared with CYP2D6*1. In stopped-flow measurements, k off for dextromethorphan was decreased in the absence of Trp75. Our results support a role for Trp75 in substrate shuttling to the active site of CYP2D6. SIGNIFICANCE STATEMENT: Using combined molecular dynamics and in vitro assays, this study shows for the first time a role for Trp75 as a channel entrance gating residue in the mechanism of substrate binding/unbinding in CYP2D6. Energy landscapes derived from molecular dynamics were used to quantitate the strength of gating, and kinetics assays showed the impact on enzyme efficiency and k off of a Trp75Ala mutation.

Crystal Structures of Drug-Metabolizing CYPs

The complex enzyme kinetics displayed by drug-metabolizing cytochrome P450 enzymes (CYPs) (see Chapter 9 ) can, in part, be explained by an examination of their crystallographic protein structures. Fortunately, despite low sequence similarity between different families of drug-metabolizing CYPs, there exists a high degree of structural homology within the superfamily. This similarity in the protein fold allows for a direct comparison of the structural features of CYPs that contribute toward differences in substrate binding, heterotropic and homotropic cooperativity, and genetic variability in drug metabolism. In this chapter, we first provide an overview of the nomenclature and the role of structural features that are common in all CYPs. We then apply these definitions to understand the different substrate specificities and functions in the CYP3A, CYP2C, and CYP2D families of enzymes.

Roles of the variable P450 substrate recognition sites SRS1 and SRS6 in esfenvalerate metabolism by CYP6AE subfamily enzymes in Helicoverpa armigera

The cotton bollworm P450s of the clustered CYP6AE subfamily share high sequence identities but differ dramatically in their capacity to metabolize xenobiotics, especially esfenvalerate. Among them, CYP6AE17 has the highest sequence identity with CYP6AE18 but shows ~7-fold higher metabolic efficiency. CYP6AE11 is most active towards esfenvalerate but CYP6AE20 is inactive even though the enzymes share 54.8% sequence identity. Sequence analysis revealed the SRS1 (Substrate Recognition Site) and SRS6 between CYP6AE17 and CYP6AE18, and SRS1 between CYP6AE11 and CYP6AE20 are the most variable among all six SRSs. In order to identify the key factors that underlie the observed catalytic difference, we exchanged these SRS sequences between two pairs of P450s and studied the activity of the resulting hybrid mutants or chimeras. In vitro metabolism showed that the CYP6AE17/18 chimeras had 2- and 14-fold decreased activities and the CYP6AE18/17 chimeras had 6- and 10-fold increased activities to esfenvalerate. Meanwhile, after exchanging SRS1 with each other, the CYP6AE11/20 chimera folded incorrectly but the CYP6AE20/11 chimera gained moderate activity to esfenvalerate. Molecular modelling showed that amino acids variants within SRS1 or SRS6 change the shape and chemical environment of the active sites, which may affect the ligand-binding interactions. These results indicate that the protein structure variation resulting from the sequence diversity of SRSs promotes the evolution of insect chemical defense and contributes to the development of insect resistance to pesticides.

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.

The marmoset cytochrome P450 superfamily: Sequence/phylogenetic analyses, genomic structure, and catalytic function

The common marmoset (Callithrix jacchus) is a New World monkey that has attracted much attention as a potentially useful primate model for preclinical testing. A total of 36 marmoset cytochrome P450 (P450) isoforms in the P450 1-51 subfamilies have been identified and characterized by the application of genome analysis and molecular functional characterization. In this mini-review, we provide an overview of the genomic structures, sequence identities, and substrate selectivities of marmoset P450s compared with those of human P450s. Based on the sequence identity, phylogeny, and genomic organization of marmoset P450s, orthologous relationships were established between human and marmoset P450s. Twenty-four members of the marmoset P450 1A, 2A, 2B, 2C, 2D, 2E, 3A, 4A, and 4F subfamilies shared high degrees of homology in terms of cDNA (>89%) and amino acid sequences (>85%) with the corresponding human P450s; P450 2C76 was among the exceptions. Phylogenetic analysis using amino acid sequences revealed that marmoset P450s in the P450 1-51 families were located in the same clades as their human and macaque P450 homologs. This finding underlines the evolutionary closeness of marmoset P450s to their human and macaque homologs. Most marmoset P450 1-4 enzymes catalyzed the typical drug-metabolizing reactions of the corresponding human P450 homologs, except for some differences of P450 2A6 and 2B6. Consequently, it appears that the substrate specificities of enzymes in the P450 1-4 families are generally similar in marmosets and humans. The information presented here supports a better understanding of the functional characteristics of marmoset P450s and their similarities and differences with human P450s. It is hoped that this mini-review will facilitate the successful use of marmosets as primate models in drug metabolism and pharmacokinetic studies.

Biotransformation, Using Recombinant CYP450-Expressing Baker's Yeast Cells, Identifies a Novel CYP2D6.10A122V Variant Which Is a Superior Metabolizer of Codeine to Morphine Than the Wild-Type Enzyme

CYP2D6, a cytochrome P450 (CYP) enzyme, metabolizes codeine to morphine. Within the human body, 0-15% of codeine undergoes O-demethylation by CYP2D6 to form morphine, a far stronger analgesic than codeine. Genetic polymorphisms in wild-type CYP2D6 (CYP2D6-wt) are known to cause poor-to-extensive metabolism of codeine and other CYP2D6 substrates. We have established a platform technology that allows stable expression of human CYP genes from chromosomal loci of baker's yeast cells. Four CYP2D6 alleles, (i) chemically synthesized CYP2D6.1, (ii) chemically synthesized CYP2D6-wt, (iii) chemically synthesized CYP2D6.10, and (iv) a novel CYP2D6.10 variant CYP2D6-C (i.e., CYP2D6.10^A122V) isolated from a liver cDNA library, were cloned for chromosomal integration in yeast cells. When expressed in yeast, CYP2D6.10 enzyme shows weak activity compared with CYP2D6-wt and CYP2D6.1 which have moderate activity, as reported earlier. Surprisingly, however, the CYP2D6-C enzyme is far more active than CYP2D6.10. More surprisingly, although CYP2D6.10 is a known low metabolizer of codeine, yeast cells expressing CYP2D6-C transform >70% of codeine to morphine, which is more than twice that of cells expressing the extensive metabolizers, CYP2D6.1, and CYP2D6-wt. The latter two enzymes predominantly catalyze formation of codeine's N-demethylation product, norcodeine, with >55% yield. Molecular modeling studies explain the specificity of CYP2D6-C for O-demethylation, validating observed experimental results. The yeast-based CYP2D6 expression systems, described here, could find generic use in CYP2D6-mediated drug metabolism and also in high-yield chemical reactions that allow the formation of regio-specific dealkylation products.

Molecular Determinants of Substrate Affinity and Enzyme Activity of a Cytochrome P450BM3 Variant

Cytochrome P450_BM3 catalyzes the hydroxylation and/or epoxidation of fatty acids, fatty amides, and alcohols. Protein engineering has produced P450_BM3 variants capable of accepting drug molecules normally metabolized by human P450 enzymes. The enhanced substrate promiscuity has been attributed to the greater flexibility of the lid of the substrate channel. However, it is not well understood how structurally different and highly polar drug molecules can stably bind in the active site nor how the activity and coupling efficiency of the enzyme may be affected by the lack of enzyme-substrate complementarity. To address these important aspects of non-native small molecule binding, this study investigated the binding of drug molecules with different size, charge, polar surface area, and human P450 affinity on the promiscuous R47L/F87V/L188Q/E267V/F81I pentuple mutant of P450_BM3. Binding free energy data and energy decomposition analysis showed that pentuple mutant P450_BM3 stably binds (i.e., negative ΔG_b°) a broad range of substrate and inhibitor types because dispersion interactions with active site residues overcome unfavorable repulsive and electrostatic effects. Molecular dynamics simulations revealed that 1) acidic substrates tend to disrupt the heme propionate A-K69 salt bridge, which may reduce heme oxidizing ability, and 2) the lack of complementarity leads to high substrate mobility and water density in the active site, which may lead to uncoupling. These factors must be considered in future developments of P450_BM3 as a biocatalyst in the large-scale production of drug metabolites.

Microsecond MD simulations of human CYP2D6 wild-type and five allelic variants reveal mechanistic insights on the function

Characterization of cytochrome P450 2D6 (CYP2D6) and the impact of the major identified allelic variants on the activity of one of the most dominating drug-metabolising enzymes is essential to increase drug safety and avoid adverse reactions. Microsecond molecular dynamics simulations have been performed to capture the dynamic signatures of this complex enzyme and five allelic variants with diverse enzymatic activity. In addition to the apo simulations, three substrates (bufuralol, veliparib and tamoxifen) and two inhibitors (prinomastat and quinidine) were included to explore their influence on the structure and dynamical features of the enzyme. Our results indicate that the altered enzyme activity can be attributed to changes in the hydrogen bonding network within the active site, and local structural differences in flexibility, position and shape of the binding pocket. In particular, the increased (CYP2D6*53) or the decreased (CYP2D6*17) activity seems to be related to a change in dynamics of mainly the BC loop due to a modified hydrogen bonding network around this region. In addition, the smallest active site volume was found for CYP2D6*4 (no activity). CYP2D6*2 (normal activity) showed no major differences in dynamic behaviour compared to the wild-type.

Celecoxib is a substrate of CYP2D6: Impact on celecoxib metabolism in individuals with CYP2C9*3 variants

Celecoxib was characterized as a substrate of human cytochrome P450 (CYP) 2D6 in vitro. In recombinant CYP2D6, celecoxib hydroxylation showed atypical substrate inhibition kinetics with apparent K_m, K_i, and V_max of 67.2 μM, 12.6 μM, and 1.33 μM/min, respectively. In human liver microsomes (HLMs), a concentration-dependent inhibition of celecoxib hydroxylation by quinidine was observed after CYP2C9 and CYP3A4 were inhibited. In individual HLMs with variable CYP2D6 activities, a significant correlation was observed between celecoxib hydroxylation and CYP2D6-selective dextromethorphan O-demethylation when CYP2C9 and CYP3A4 activities were suppressed (r = 0.97, P < 0.0001). Molecular modeling showed two predominant docking modes of celecoxib with CYP2D6, resulting in either a substrate or an inhibitor. A second allosteric binding antechamber, which stabilized the inhibition mode, was revealed. Modeling results were consistent with the observed substrate inhibition kinetics. Using HLMs from individual donors, the relative contribution of CYP2D6 to celecoxib metabolism was found to be highly variable and dependent on CYP2C9 genotypes, ranging from no contribution in extensive metabolizers with CYP2C9*1*1 genotype to approximately 30% in slow metabolizers with allelic variants CYP2C9*1*3 and CYP2C9*3*3. These results demonstrate that celecoxib may become a potential victim of CYP2D6-associated drug-drug interactions, particularly in individuals with reduced CYP2C9 activity.

Assessing the potential of pharmaceuticals and their transformation products to cause mutagenic effects: Implications for gene expression profiling

The selection and prioritization of pharmaceuticals and their transformation products for evaluating effects on the environment and human health is a challenging task. One common approach is based on compounds (e.g., mixture composition, concentrations), and another on biology (e.g., relevant endpoint, biological organizational level). Both of these approaches often resemble a Lernaean Hydra-they can create more questions than answers. The present study embraces this complexity, providing an integrated approach toward assessing the potential effects of transformation products of pharmaceuticals by means of mutagenicity, estrogenicity, and differences in the gene expression profiles. Mutagenicity using the tk kinase assay was applied to assess a list of 11 priority pharmaceuticals, namely, atenolol, azithromycin, carbamazepine, diclofenac, ibuprofen, erythromycin, metoprolol, ofloxacin, propranolol, sulfamethoxazole, and trimethoprim. The most mutagenic compounds were found to be β-blockers. In parallel, the photolabile pharmaceuticals were assessed for their mixture effects on mutagenicity (tk assay), estrogenicity (T47D- KBluc assay), and gene expression (microarrays). Interestingly, the mixtures were mutagenic at the µg/L level, indicating a synergistic effect. None of the photolysed mixtures were statistically significantly estrogenic. Gene expression profiling revealed effects related mainly to certain pathways, those of the p53 gene, mitogen-activated protein kinase, alanine, aspartate, and glutamate metabolism, and translation-related (spliceosome). Fourteen phototransformation products are proposed based on the m/z values found through ultra-performance liquid chromatography-tandem mass spectrometry analysis. The transformation routes of the photolysed mixtures indicate a strong similarity with those obtained for each pharmaceutical separately. This finding reinforces the view that transformation products are to be expected in naturally occurring mixtures. Environ Toxicol Chem 2016;35:2753-2764. © 2016 SETAC.

Recent Structural Insights into Cytochrome P450 Function

Cytochrome P450 (P450) enzymes are important in the metabolism of drugs, steroids, fat-soluble vitamins, carcinogens, pesticides, and many other types of chemicals. Their catalytic activities are important issues in areas such as drug-drug interactions and endocrine function. During the past 30 years, structures of P450s have been very helpful in understanding function, particularly the mammalian P450 structures available in the past 15 years. We review recent activity in this area, focusing on the past 2 years (2014-2015). Structural work with microbial P450s includes studies related to the biosynthesis of natural products and the use of parasitic and fungal P450 structures as targets for drug discovery. Studies on mammalian P450s include the utilization of information about 'drug-metabolizing' P450s to improve drug development and also to understand the molecular bases of endocrine dysfunction.

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.

Marmoset cytochrome P450 2D8 in livers and small intestines metabolizes typical human P450 2D6 substrates, metoprolol, bufuralol and dextromethorphan

1. Although the New World non-human primate, the common marmoset (Callithrix jacchus), is a potentially useful animal model, comprehensive understanding of drug metabolizing enzymes is insufficient. 2. A cDNA encoding a novel cytochrome P450 (P450) 2D8 was identified in marmosets. The amino acid sequence deduced from P450 2D8 cDNA showed a high sequence identity (83-86%) with other primate P450 2Ds. Phylogenetic analysis showed that marmoset P450 2D8 was closely clustered with human P450 2D6, unlike P450 2Ds of miniature pig, dog, rabbit, guinea pig, mouse or rat. 3. Marmoset P450 2D8 mRNA was predominantly expressed in the liver and small intestine among the tissues types analyzed, whereas marmoset P450 2D6 mRNA was expressed predominantly in the liver where P450 2D protein was detected by immunoblotting. 4. By metabolic assays using marmoset P450 2D8 protein heterologously expressed in Escherichia coli, although P450 2D8 exhibits lower catalytic efficiency compared to marmoset and human P450 2D6 enzymes, P450 2D8 mediated O-demethylations of metoprolol and dextromethorphan and bufuralol 1'-hydroxylation. 5. These results suggest that marmoset P450 2D8 (also expressed in the extrahepatic tissues) has potential roles in drug metabolism in a similar manner to those of human and marmoset P450 2D6.

Double site saturation mutagenesis of the human cytochrome P450 2D6 results in regioselective steroid hydroxylation

The human cytochrome P450 2D6 (CYP2D6) is one of the major human drug metabolizing enzymes and acts preferably on substrates containing a basic nitrogen atom. Testosterone - just as other steroids - is an atypical substrate and only poorly metabolized by CYP2D6. The present study intended to investigate the influence of the two active site residues 216 and 483 on the capability of CYP2D6 to hydroxylate steroids such as for example testosterone. All 400 possible combinatorial mutations at these two positions have been generated and expressed individually in Pichia pastoris. Employing whole-cell biotransformations coupled with HPLC-MS analysis the testosterone hydroxylase activity and regioselectivity of every single CYP2D6 variant was determined. Covering the whole sequence space, CYP2D6 variants with improved activity and so far unknown regio-preference in testosterone hydroxylation were identified. Most intriguingly and in contrast to previous literature reports about mutein F483I, the mutation F483G led to preferred hydroxylation at the 2β-position, while the slow formation of 6β-hydroxytestosterone, the main product of wild-type CYP2D6, was further reduced. Two point mutations have already been sufficient to convert CYP2D6 into a steroid hydroxylase with the highest ever reported testosterone hydroxylation rate for this enzyme, which is of the same order of magnitude as for the conversion of the standard substrate bufuralol by wild-type CYP2D6. Furthermore, this study is also an example for efficient human CYP engineering in P. pastoris for biocatalytic applications and to study so far unknown pharmacokinetic effects of individual and combined mutations in these key enzymes of the human drug metabolism.

Rationalization of stereospecific binding of propranolol to cytochrome P450 2D6 by free energy calculations

Abstract

Cytochrome P450 2D6 is a major drug-metabolising enzyme with a wide substrate range. A single-point mutation introduced in this enzyme induces stereoselective binding of R and S-propranolol whereas the wild type has no preference. The system has previously been studied both experimentally and computationally (de Graaf et al. in Eur Biophys J 36:589–599, 2007a). The in silico study reported hysteresis and significant deviations from closure of thermodynamic cycles, probably because of lack of sampling. Here, we focus on the effect of prolonged simulation time and enhanced sampling methods, such as Hamiltonian replica exchange, to reduce these problems and to improve the precision of free energy calculations. Finally we rationalize the results at a molecular level and compare data with experimental findings and previously estimated free energies.

Graphical Abstract

Propranolol (PPD) binds stereospecifically to cytochrome P450 2D6, if the mutation F483A is introduced, by reducing binding affinity to R-PPD. This observation has previously been studied by molecular dynamics, by use of several different thermodynamic cycles. Previous computational results suffered from internal inconsistencies, probably because of insufficient sampling. Here we use prolonged simulation time and Hamiltonian replica exchange to overcome the inconsistencies, and give a possible rationalization of the phenomenon at a molecular level.

Steroid biotransformations in biphasic systems with Yarrowia lipolytica expressing human liver cytochrome P450 genes

BACKGROUND

Yarrowia lipolytica efficiently metabolizes and assimilates hydrophobic compounds such as n-alkanes and fatty acids. Efficient substrate uptake is enabled by naturally secreted emulsifiers and a modified cell surface hydrophobicity and protrusions formed by this yeast. We were examining the potential of recombinant Y. lipolytica as a biocatalyst for the oxidation of hardly soluble hydrophobic steroids. Furthermore, two-liquid biphasic culture systems were evaluated to increase substrate availability. While cells, together with water soluble nutrients, are maintained in the aqueous phase, substrates and most of the products are contained in a second water-immiscible organic solvent phase.

RESULTS

For the first time we have co-expressed the human cytochromes P450 2D6 and 3A4 genes in Y. lipolytica together with human cytochrome P450 reductase (hCPR) or Y. lipolytica cytochrome P450 reductase (YlCPR). These whole-cell biocatalysts were used for the conversion of poorly soluble steroids in biphasic systems.Employing a biphasic system with the organic solvent and Y. lipolytica carbon source ethyl oleate for the whole-cell bioconversion of progesterone, the initial specific hydroxylation rate in a 1.5 L stirred tank bioreactor was further increased 2-fold. Furthermore, the product formation was significantly prolonged as compared to the aqueous system. Co-expression of the human CPR gene led to a 4-10-fold higher specific activity, compared to the co-overexpression of the native Y. lipolytica CPR gene. Multicopy transformants showed a 50-70-fold increase of activity as compared to single copy strains.

CONCLUSIONS

Alkane-assimilating yeast Y. lipolytica, coupled with the described expression strategies, demonstrated its high potential for biotransformations of hydrophobic substrates in two-liquid biphasic systems. Especially organic solvents which can be efficiently taken up and/or metabolized by the cell might enable more efficient bioconversion as compared to aqueous systems and even enable simple, continuous or at least high yield long time processes.

Functional characterization of a first avian cytochrome P450 of the CYP2D subfamily (CYP2D49)

The CYP2D family members are instrumental in the metabolism of 20-25% of commonly prescribed drugs. Although many CYP2D isoforms have been well characterized in other animal models, research concerning the chicken CYP2Ds is limited. In this study, a cDNA encoding a novel CYP2D enzyme (CYP2D49) was cloned from the chicken liver for the first time. The CYP2D49 cDNA contained an open reading frame of 502 amino acids that shared 52%-57% identities with other CYP2Ds. The gene structure and neighboring genes of CYP2D49 are conserved and similar to those of human CYP2D6. Additionally, similar to human CYP2D6, CYP2D49 is un-inducible in the liver and expressed predominantly in the liver, kidney and small intestine, with detectable levels in several other tissues. Metabolic assays of the CYP2D49 protein heterologously expressed in E. coli and Hela cells indicated that CYP2D49 metabolized the human CYP2D6 substrate, bufuralol, but not debrisoquine. Moreover, quinidine, a potent inhibitor of human CYP2D6, only inhibited the bufuralol 1'-hydroxylation activity of CYP2D49 to a negligible degree. All these results indicated that CYP2D49 had functional characteristics similar to those of human CYP2D6 but measurably differed in the debrisoquine 4'-hydroxylation and quinidine inhibitory profile. Further structure-function investigations that employed site-directed mutagenesis and circular dichroism spectroscopy identified the importance of Val-126, Glu-222, Asp-306, Phe-486 and Phe-488 in keeping the enzymatic activity of CYP2D49 toward bufuralol as well as the importance of Asp-306, Phe-486 and Phe-488 in maintaining the conformation of CYP2D49 protein. The current study is only the first step in characterizing the metabolic mechanism of CYP2D49; further studies are still required.

Ligand-Based Site of Metabolism Prediction for Cytochrome P450 2D6

A ligand-based method based on the SMARTCyp approach that predicts the sites of cytochrome P450 2D6-mediated metabolism of druglike molecules has been developed. The method uses only two descriptors besides the reactivity from SMARTCyp: the distance to a protonated nitrogen atom and the distance to the end of the molecule. Hence, the site of metabolism is predicted directly from the 2D structure of a molecule, without requiring calculation of electronic properties or generation of 3D structures. Testing on an independent test set gives an area under the curve value of 0.94, and a site of metabolism is found among the top two ranked atoms for 91% of the compounds.

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.

Prediction of CYP2D6 drug interactions from in vitro data: evidence for substrate-dependent inhibition

Predicting the magnitude of potential drug-drug interactions is important for underwriting patient safety in the clinical setting. Substrate-dependent inhibition of cytochrome P450 enzymes may confound extrapolation of in vitro results to the in vivo situation. However, the potential for substrate-dependent inhibition with CYP2D6 has not been well characterized. The inhibition profiles of 20 known inhibitors of CYP2D6 were characterized in vitro against four clinically relevant CYP2D6 substrates (desipramine, dextromethorphan, metoprolol, and thioridazine) and bufuralol. Dextromethorphan exhibited the highest sensitivity to in vitro inhibition, whereas metoprolol was the least sensitive. In addition, when metoprolol was the substrate, inhibitors with structurally constrained amino moieties (clozapine, debrisoquine, harmine, quinidine, and yohimbine) exhibited at least a 5-fold decrease in inhibition potency when results were compared with those for dextromethorphan. Atypical inhibition kinetics were observed for these and other inhibitor-substrate pairings. In silico docking studies suggested that interactions with Glu216 and an adjacent hydrophobic binding pocket may influence substrate sensitivity and inhibition potency for CYP2D6. The in vivo sensitivities of the clinically relevant CYP2D6 substrates desipramine, dextromethorphan, and metoprolol were determined on the basis of literature drug-drug interaction (DDI) outcomes. Similar to the in vitro results, dextromethorphan exhibited the highest sensitivity to CYP2D6 inhibition in vivo. Finally, the magnitude of in vivo CYP2D6 DDIs caused by quinidine was predicted using desipramine, dextromethorphan, and metoprolol. Comparisons of the predictions with literature results indicated that the marked decrease in inhibition potency observed for the metoprolol-quinidine interaction in vitro translated to the in vivo situation.

Does compound I vary significantly between isoforms of cytochrome P450?

The cytochrome P450 (CYP) enzymes are important in many areas, including pharmaceutical development. Subtle changes in the electronic structure of the active species, Compound I, have been postulated previously to account partly for the experimentally observed differences in reactivity between isoforms. Current predictive models of CYP metabolism typically assume an identical Compound I in all isoforms. Here we present a method to calculate the electronic structure and to estimate the Fe-O bond enthalpy of Compound I, and apply it to several human and bacterial CYP isoforms. Conformational flexibility is accounted for by sampling large numbers of structures from molecular dynamics simulations, which are subsequently optimized with density functional theory (B3LYP) based quantum mechanics/molecular mechanics. The observed differences in Compound I between human isoforms are small: They are generally smaller than the spread of values obtained for the same isoform starting from different initial structures. Hence, it is unlikely that the variation in activity between human isoforms is due to differences in the electronic structure of Compound I. A larger difference in electronic structure is observed between the human isoforms and P450(cam) and may be explained by the slightly different hydrogen-bonding environment surrounding the cysteinyl sulfur. The presence of substrate in the active site of all isoforms studied appears to cause a slight decrease in the Fe-O bond enthalpy, apparently due to displacement of water out of the active site, suggesting that Compound I is less stable in the presence of substrate.

Understanding the determinants of selectivity in drug metabolism through modeling of dextromethorphan oxidation by cytochrome P450

Cytochrome P450 enzymes play key roles in the metabolism of the majority of drugs. Improved models for prediction of likely metabolites will contribute to drug development. In this work, two possible metabolic routes (aromatic carbon oxidation and O-demethylation) of dextromethorphan are compared using molecular dynamics (MD) simulations and density functional theory (DFT). The DFT results on a small active site model suggest that both reactions might occur competitively. Docking and MD studies of dextromethorphan in the active site of P450 2D6 show that the dextromethorphan is located close to heme oxygen in a geometry apparently consistent with competitive metabolism. In contrast, calculations of the reaction path in a large protein model [using a hybrid quantum mechanical-molecular mechanics (QM/MM) method] show a very strong preference for O-demethylation, in accordance with experimental results. The aromatic carbon oxidation reaction is predicted to have a high activation energy, due to the active site preventing formation of a favorable transition-state structure. Hence, the QM/MM calculations demonstrate a crucial role of many active site residues in determining reactivity of dextromethorphan in P450 2D6. Beyond substrate binding orientation and reactivity of Compound I, successful metabolite predictions must take into account the detailed mechanism of oxidation in the protein. These results demonstrate the potential of QM/MM methods to investigate specificity in drug metabolism.

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

Brain chemistry: how does P450 catalyze the O-demethylation reaction of 5-methoxytryptamine to yield serotonin?

Density functional theory has been applied to elucidate the mechanism of the O-demethylation reaction that generates serotonin from 5-methoxytryptamine (5-MT); a process that is efficiently catalyzed by P450 CYP2D6. Two substrates, the neutral 5-MT and the protonated 5-MTH(+), were used to probe the reactivity of CYP2D6 compound I. Notably, the H-abstraction process is found to be slightly more facile for 5-MT. However, our DFT augmented by docking results show that the amino acid Glu216 in the active site holds the NH(3)(+) tail of the 5-MTH(+) substrate in an upright conformation and thereby controls the regioselectivity of the bond activation. Thus, the substrate protonation serves an important function in maximizing the yield of serotonin. This finding is in accord with experimental conclusions that 5-MTH(+) serves as the substrate for the CYP2D6 enzyme. The study further shows that the H-abstraction follows two-state reactivity (TSR), whereas the rebound path may involve more states due to the appearance of both Fe(IV) and Fe(III) electromers during the reaction of 5-MTH(+).

Deletion of microsomal cytochrome b5 profoundly affects hepatic and extrahepatic drug metabolism

We demonstrated recently that cytochrome b(5) plays an important in vivo role in hepatic cytochrome P450 (P450) function [J Biol Chem 283:31385-31393, 2008]. We have now generated a model in which cytochrome b(5) has been deleted in all tissues [cytochrome b(5) complete null (BCN)], which surprisingly results in a viable mouse despite the putative in vivo roles of this protein in lipid and steroid hormone metabolism and the reduction of methemoglobin. In contrast to the liver-specific deletion, complete deletion of cytochrome b(5) leads to a neonatal increase in the expression of many hepatic P450s at both the protein and mRNA level. In extrahepatic tissues, some changes in P450 expression were also observed that were isoform-dependent. In vitro cytochrome P450 activities in liver, kidney, lung, and small intestine of BCN mice were determined for a range of model substrates and probe drugs; a profound reduction in the metabolism of some substrates, particularly in lung, kidney, and small intestine, was observed. In vivo, the metabolism of metoprolol was significantly altered in BCN mice, in contrast to the previous finding in the liver-specific cytochrome b(5) deletion, suggesting that extrahepatic cytochrome b(5) plays a significant role in its disposition. Testicular Cyp17 hydroxylase and lyase activities were also significantly reduced by cytochrome b(5) deletion, leading to significantly lower levels of testicular testosterone. The BCN mouse provides an additional model system with which to further investigate the functions of cytochrome b(5), particularly in extrahepatic tissues.

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.

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.

Metabolism ofN-substituted 7-methoxy-4-(aminomethyl) -coumarins by cytochrome P450 2D6 mutants and the indication of additional substrate interaction points

Previous studies have shown the critical roles residues F120 and F483 play in the oxidative metabolism of 7-methoxy-4-(aminomethyl)-coumarin (MAMC) by cytochrome P450 2D6 (CYP2D6). In the present study, a series of N-alkyl-7-methoxy-4-(aminomethyl)-coumarins (MAMC analogues) were used as substrates for the F120A and F483A mutants in order to probe the CYP2D6 active site. The F120A and F483A mutants of CYP2D6 displayed significant activity towards the MAMC analogues. Automated docking studies of the MAMC analogues in a CYP2D6 homology model suggested a distal hydrophobic active site binding cleft for the substrate N-alkyl chains, consisting of the residues L213 and V308.

The mechanism causing the difference in kinetic properties between rat CYP2D4 and human CYP2D6 in the oxidation of dextromethorphan and bufuralol

The capacity to oxidize bufuralol (BF) and dextromethorphan (DEX) was compared kinetically between human CYP2D6 and four rat CYP2D (CYP2D1, -2D2, -2D3 and -2D4) isoenzymes in a yeast cell expression system. In BF 1''-hydroxylation and DEX O-demethylation, only CYP2D4 showed hook-shaped Eadie-Hofstee plots, the other four CYP2D enzymes exhibiting linear plots. In DEX N-demethylation, rat CYP2D2 did not show any detectable activity under the conditions used, whereas the other four enzymes yielded linear Eadie-Hofstee plots. To elucidate the mechanisms causing the nonlinear kinetics, four CYP2D4 mutants, CYP2D4-F109I, -V123F, -L216F and -A486F, were prepared. CYP2D4-V123F, -L216F and -A486F yielded linear or linear-like Eadie-Hofstee plots for BF 1''-hydroxylation, whereas only CYP2D4-A486F exhibited linear plots for DEX O-demethylation. The substitution of Phe-109 by isoleucine did not have any effect on the oxidative capacity of CYP2D4 for either BF or DEX. These results suggest that the introduction of phenylalanine in the active-site cavity of CYP2D4 simplifies complicated interactions between the substrates and the amino acid residues, but the mechanisms causing the simplification differ between BF and DEX.

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.

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

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