Why Is Quinidine an Inhibitor of Cytochrome P450 2D6?

In vitro comparative analysis of metabolic capabilities and inhibitory profiles of selected CYP2D6 alleles on tramadol metabolism

Tramadol, the 41st most prescribed drug in the United States in 2021 is a prodrug activated by CYP2D6, which is highly polymorphic. Previous studies showed enzyme-inhibitor affinity varied between different CYP2D6 allelic variants with dextromethorphan and atomoxetine metabolism. However, no study has compared tramadol metabolism in different CYP2D6 alleles with different CYP2D6 inhibitors. We hypothesize that the inhibitory effects of CYP2D6 inhibitors on CYP2D6-mediated tramadol metabolism are inhibitor- and CYP2D6-allele-specific. We performed comparative analyses of CYP2D6*1, CYP2D6*2, CYP2D6*10, and CYP2D6*17 using recombinant enzymes to metabolize tramadol to O-desmethyltramadol, measured via UPLC-MS/MS. The Michaelis constant (Km) and maximum velocity (Vmax) for each CYP2D6 allele, and IC50 values for different inhibitors were determined by nonlinear regression analysis. Intrinsic clearance was calculated as Vmax/Km. The intrinsic clearance of tramadol was almost double for CYP2D6*2 (180%) but was much lower for CYP2D6*10 and *17 (20% and 10%, respectively) compared to CYP2D6*1. The inhibitor potencies (defined by Ki) for the various inhibitors for the CYP2D6*1 allele were quinidine > terbinafine > paroxetine ≈ duloxetine >>bupropion. CYP2D6*2 showed the next greatest inhibition, with Ki ratios compared to CYP2D6*1 ranging from 0.96 to 3.87. For each inhibitor tested, CYP2D6*10 and CYP2D6*17 were more resistant to inhibition than CYP2D6*1 or CYP2D6*2, with most Ki ratios in the 3-9 range. Three common CYP2D6 allelic variants showed different metabolic capacities toward tramadol and genotype-dependent inhibition compared to CYP2D6*1. Further studies are warranted to understand the clinical consequences of inhibitor and CYP2D6 genotype-dependent drug-drug interactions on tramadol bioactivation.

Physiologically-based pharmacokinetic modeling of quinidine to establish a CYP3A4, P-gp, and CYP2D6 drug-drug-gene interaction network

The antiarrhythmic agent quinidine is a potent inhibitor of cytochrome P450 (CYP) 2D6 and P-glycoprotein (P-gp) and is therefore recommended for use in clinical drug-drug interaction (DDI) studies. However, as quinidine is also a substrate of CYP3A4 and P-gp, it is susceptible to DDIs involving these proteins. Physiologically-based pharmacokinetic (PBPK) modeling can help to mechanistically assess the absorption, distribution, metabolism, and excretion processes of a drug and has proven its usefulness in predicting even complex interaction scenarios. The objectives of the presented work were to develop a PBPK model of quinidine and to integrate the model into a comprehensive drug-drug(-gene) interaction (DD(G)I) network with a diverse set of CYP3A4 and P-gp perpetrators as well as CYP2D6 and P-gp victims. The quinidine parent-metabolite model including 3-hydroxyquinidine was developed using pharmacokinetic profiles from clinical studies after intravenous and oral administration covering a broad dosing range (0.1-600 mg). The model covers efflux transport via P-gp and metabolic transformation to either 3-hydroxyquinidine or unspecified metabolites via CYP3A4. The 3-hydroxyquinidine model includes further metabolism by CYP3A4 as well as an unspecific hepatic clearance. Model performance was assessed graphically and quantitatively with greater than 90% of predicted pharmacokinetic parameters within two-fold of corresponding observed values. The model was successfully used to simulate various DD(G)I scenarios with greater than 90% of predicted DD(G)I pharmacokinetic parameter ratios within two-fold prediction success limits. The presented network will be provided to the research community and can be extended to include further perpetrators, victims, and targets, to support investigations of DD(G)Is.

Human Cytochrome P450 1, 2, 3 Families as Pharmacogenes with Emphases on Their Antimalarial and Antituberculosis Drugs and Prevalent African Alleles

Precision medicine gives individuals tailored medical treatment, with the genotype determining the therapeutic strategy, the appropriate dosage, and the likelihood of benefit or toxicity. Cytochrome P450 (CYP) enzyme families 1, 2, and 3 play a pivotal role in eliminating most drugs. Factors that affect CYP function and expression have a major impact on treatment outcomes. Therefore, polymorphisms of these enzymes result in alleles with diverse enzymatic activity and drug metabolism phenotypes. Africa has the highest CYP genetic diversity and also the highest burden of malaria and tuberculosis, and this review presents current general information on CYP enzymes together with variation data concerning antimalarial and antituberculosis drugs, while focusing on the first three CYP families. Afrocentric alleles such as CYP2A6*17, CYP2A6*23, CYP2A6*25, CYP2A6*28, CYP2B6*6, CYP2B6*18, CYP2C8*2, CYP2C9*5, CYP2C9*8, CYP2C9*9, CYP2C19*9, CYP2C19*13, CYP2C19*15, CYP2D6*2, CYP2D6*17, CYP2D6*29, and CYP3A4*15 are implicated in diverse metabolic phenotypes of different antimalarials such as artesunate, mefloquine, quinine, primaquine, and chloroquine. Moreover, CYP3A4, CYP1A1, CYP2C8, CYP2C18, CYP2C19, CYP2J2, and CYP1B1 are implicated in the metabolism of some second-line antituberculosis drugs such as bedaquiline and linezolid. Drug-drug interactions, induction/inhibition, and enzyme polymorphisms that influence the metabolism of antituberculosis, antimalarial, and other drugs, are explored. Moreover, a mapping of Afrocentric missense mutations to CYP structures and a documentation of their known effects provided structural insights, as understanding the mechanism of action of these enzymes and how the different alleles influence enzyme function is invaluable to the advancement of precision medicine.

DEEPCYPs: A deep learning platform for enhanced cytochrome P450 activity prediction

Cytochrome P450 (CYP) is a superfamily of heme-containing oxidizing enzymes involved in the metabolism of a wide range of medicines, xenobiotics, and endogenous compounds. Five of the CYPs (1A2, 2C9, 2C19, 2D6, and 3A4) are responsible for metabolizing the vast majority of approved drugs. Adverse drug-drug interactions, many of which are mediated by CYPs, are one of the important causes for the premature termination of drug development and drug withdrawal from the market. In this work, we reported in silicon classification models to predict the inhibitory activity of molecules against these five CYP isoforms using our recently developed FP-GNN deep learning method. The evaluation results showed that, to the best of our knowledge, the multi-task FP-GNN model achieved the best predictive performance with the highest average AUC (0.905), F1 (0.779), BA (0.819), and MCC (0.647) values for the test sets, even compared to advanced machine learning, deep learning, and existing models. Y-scrambling testing confirmed that the results of the multi-task FP-GNN model were not attributed to chance correlation. Furthermore, the interpretability of the multi-task FP-GNN model enables the discovery of critical structural fragments associated with CYPs inhibition. Finally, an online webserver called DEEPCYPs and its local version software were created based on the optimal multi-task FP-GNN model to detect whether compounds bear potential inhibitory activity against CYPs, thereby promoting the prediction of drug-drug interactions in clinical practice and could be used to rule out inappropriate compounds in the early stages of drug discovery and/or identify new CYPs inhibitors.

Treatment of Erythroid Precursor Cells from β-Thalassemia Patients with Cinchona Alkaloids: Induction of Fetal Hemoglobin Production

β-thalassemias are among the most common inherited hemoglobinopathies worldwide and are the result of autosomal mutations in the gene encoding β-globin, causing an absence or low-level production of adult hemoglobin (HbA). Induction of fetal hemoglobin (HbF) is considered to be of key importance for the development of therapeutic protocols for β-thalassemia and novel HbF inducers need to be proposed for pre-clinical development. The main purpose on this study was to analyze Cinchona alkaloids (cinchonidine, quinidine and cinchonine) as natural HbF-inducing agents in human erythroid cells. The analytical methods employed were Reverse Transcription quantitative real-time PCR (RT-qPCR) (for quantification of γ-globin mRNA) and High Performance Liquid Chromatography (HPLC) (for analysis of the hemoglobin pattern). After an initial analysis using the K562 cell line as an experimental model system, showing induction of hemoglobin and γ-globin mRNA, we verified whether the two more active compounds, cinchonidine and quinidine, were able to induce HbF in erythroid progenitor cells isolated from β-thalassemia patients. The data obtained demonstrate that cinchonidine and quinidine are potent inducers of γ-globin mRNA and HbF in erythroid progenitor cells isolated from nine β-thalassemia patients. In addition, both compounds were found to synergize with the HbF inducer sirolimus for maximal production of HbF. The data obtained strongly indicate that these compounds deserve consideration in the development of pre-clinical approaches for therapeutic protocols of β-thalassemia.

Physiologically-based Pharmacokinetic Modeling to Assess the Impact of CYP2D6-Mediated Drug-Drug Interactions on Tramadol and O-Desmethyltramadol Exposures via Allosteric and Competitive Inhibition

Tramadol is an opioid medication used to treat moderately severe pain. Cytochrome P450 (CYP) 2D6 inhibition could be important for tramadol, as it decreases the formation of its pharmacologically active metabolite, O‐desmethyltramadol, potentially resulting in increased opioid use and misuse. The objective of this study was to evaluate the impact of allosteric and competitive CYP2D6 inhibition on tramadol and O‐desmethyltramadol pharmacokinetics using quinidine and metoprolol as prototypical perpetrator drugs. A physiologically based pharmacokinetic model for tramadol and O‐desmethyltramadol was developed and verified in PK‐Sim version 8 and linked to respective models of quinidine and metoprolol to evaluate the impact of allosteric and competitive CYP2D6 inhibition on tramadol and O‐desmethyltramadol exposure. Our results show that there is a differentiated impact of CYP2D6 inhibitors on tramadol and O‐desmethyltramadol based on their mechanisms of inhibition. Following allosteric inhibition by a single dose of quinidine, the exposure of both tramadol (51% increase) and O‐desmethyltramadol (52% decrease) was predicted to be significantly altered after concomitant administration of a single dose of tramadol. Following multiple‐dose administration of tramadol and a single‐dose or multiple‐dose administration of quinidine, the inhibitory effect of quinidine was predicted to be long (≈42 hours) and to alter exposure of tramadol and O‐desmethyltramadol by up to 60%, suggesting that coadministration of quinidine and tramadol should be avoided clinically. In comparison, there is no predicted significant impact of metoprolol on tramadol and O‐desmethyltramadol exposure. In fact, tramadol is predicted to act as a CYP2D6 perpetrator and increase metoprolol exposure, which may necessitate the need for dose separation.

Mechanistic Biomarkers in Toxicology

Biomarkers are important parameters that are reliable, applicable, reproducible, and generally inexpensive. All biomarkers have a significant role in human health, especially mechanistic biomarkers, which are the most important for the prevention of toxic effects and diseases. They demonstrate the possibility of diagnosis, prognosis, recurrence, and spread of disease. Furthermore, they show the exposure levels to numerous chemical, biological, and physical agents. To date, the development and application of biomarkers require the knowledge of mechanisms underlying their production. Therefore, the present study focused on the possible mechanistic biomarkers.

Feature importance of machine learning prediction models shows structurally active part and important physicochemical features in drug design

The objective of this study was to obtain the indicators of physicochemical parameters and structurally active sites to design new chemical entities with desirable pharmacokinetic profiles by investigating the process by which machine learning prediction models arrive at their decisions, which are called explainable artificial intelligence. First, we developed the prediction models for metabolic stability, CYP inhibition, and P-gp and BCRP substrate recognition using 265 physicochemical parameters for designing the molecular structures. Four important parameters, including the well-known indicator h_logD, are common in some in vitro studies; as such, these can be used to optimize compounds simultaneously to address multiple pharmacokinetic concerns. Next, we developed machine learning models that had been programmed to show structurally active sites. Many types of machine learning models were developed using the results of in vitro metabolic stability study of around 30000 in-house compounds. The metabolic sites of in-house compounds predicted using some prediction models matched experimentally identified metabolically active sites, with a ratio of number of metabolic sites (predicted/actual) of over 90%. These models can be applied to several screening projects. These two approaches can be employed for obtaining lead compounds with desirable pharmacokinetic profiles efficiently.

Predicting drug metabolism and pharmacokinetics features of in-house compounds by a hybrid machine-learning model

We constructed machine learning-based pharmacokinetic prediction models with very high performance. The models were trained on 26138 and 16613 compounds involved in metabolic stability and cytochrome P450 inhibition, respectively. Because the compound features largely differed between the publicly available and in-house compounds, the models learned on the public data could not predict the in-house compounds, suggesting that outside of a certain applicability domain (AD), the prediction results are unreliable and can mislead the design of novel compounds. To exclude the uncertain prediction results, we constructed another machine learning model that determines whether the newly designed compound is inside or outside the AD. The AD was evaluated multi-dimensionally with some explanatory variables: The structural similarities and the probability obtained from the pharmacokinetic prediction model. The accuracy of predicting metabolic stability was 79.9% on the test set, increasing significantly to 93.6% after excluding the low-reliability compounds. The model properly classified the reliability of the compounds. After learning on the in-house compounds, the reliability model classified almost all (90%) of the public compounds as low reliability, indicating that the AD was properly determined by the model.

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.

Insights into the metabolic characteristics of aminopropanediol analogues of SYLs as S1P1 modulators: from structure to metabolism

SYL927 and SYL930, two aminopropanediol analogues, are novel Sphingosine-1-phosphate receptor 1 (S1P₁) modulators with higher selectivity and pharmacological activity compared with FTY720. Although the immunosuppressive activity of SYLs has been well demonstrated, information regarding the metabolic fates of the two chemicals is limited except for the CYP-catalyzed hydroxylation of SYL930. In this study, the biotransformation schemes of the two promising chemicals were investigated and compared using liver microsomes, S9 fractions and recombinant enzymes, and relevant molecular mechanism was primarily demonstrated by ligand-enzyme docking analysis (CDOCKER). As a result, the hydroxylation at alkyl chain on oxazole ring by the action of CYPs was found for both SYLs in vivo. The SULT-catalyzed sulfonation of the hydroxide was observed for SYL927 while the ADH/ALDH-catalyzed oxidation was only discovered for SYL930. The docking analysis suggested that specific non-covalent forces and/or bonding conformations of the hydroxides with biomacromolecules might be involved in the disparate metabolism of SYLs. Exploring the metabolic characteristics will help clarify the substance base for efficacy and safety of the two drugs. The uncovered structure-metabolism relationship in this study may provide an implication for the design and optimization for other S1P modulators.

New insights into quetiapine metabolism using molecular networking

Metabolism is involved in both pharmacology and toxicology of most xenobiotics including drugs. Yet, visualization tools facilitating metabolism exploration are still underused, despite the availibility of pertinent bioinformatics solutions. Since molecular networking appears as a suitable tool to explore structurally related molecules, we aimed to investigate its interest in in vitro metabolism exploration. Quetiapine, a widely prescribed antipsychotic drug, undergoes well-described extensive metabolism, and is therefore an ideal candidate for such a proof of concept. Quetiapine was incubated in metabolically competent human liver cell models (HepaRG) for different times (0 h, 3 h, 8 h, 24 h) with or without cytochrom P450 (CYP) inhibitor (ketoconazole as CYP3A4/5 inhibitor and quinidine as CYP2D6 inhibitor), in order to study its metabolism kinetic and pathways. HepaRG culture supernatants were analyzed on an ultra-high performance liquid chromatography coupled with tandem mass spectrometry (LC-HRMS/MS). Molecular networking approach on LC-HRMS/MS data allowed to quickly visualize the quetiapine metabolism kinetics and determine the major metabolic pathways (CYP3A4/5 and/or CYP2D6) involved in metabolite formation. In addition, two unknown putative metabolites have been detected. In vitro metabolite findings were confirmed in blood sample from a patient treated with quetiapine. This is the first report using LC-HRMS/MS untargeted screening and molecular networking to explore in vitro drug metabolism. Our data provide new evidences of the interest of molecular networking in drug metabolism exploration and allow our in vitro model consistency assessment.

Effect of berberine on in vitro metabolism of sulfonylureas: A herb-drug interactions study

Patients with type 2 diabetes may co-ingest herbal and prescription medicines to control their blood sugar levels. Competitive binding of drug and herb may mutually affect their metabolism. This can alter the level of drug and its kinetics in the body, potentially causing toxicities or loss of efficacy. Understanding how the metabolism of sulfonylureas like glyburide and gliclazide can be affected by the presence of berberine and vice versa can provide valuable information on the possible risk of toxicities caused by co-ingestion of drugs.

METHODS

Berberine and sulfonylureas (glyburide and gliclazide) were co-incubated with rat liver microsomes in the presence of a NADPH-regenerating system. The metabolites of berberine and sulfonylureas were analysed using liquid chromatography with high-resolution mass spectrometry in the positive ion mode. The role of individual isozymes in the metabolism of berberine, glyburide and gliclazide was investigated by using specific inhibitors.

RESULTS

In vitro metabolism of berberine led to the formation of demethyleneberberine (B1a) and its isomer B1b through demethylenation. Berberrubine (B2a) and its isomer B2b were formed through demethylation. The isozymes CYP3A and CYP2D were found to be involved in the metabolism of berberine. In vitro metabolism of glyburide and gliclazide led to the formation of hydroxylated metabolites. The isozymes CYP3A and CYP2C were found to be involved in the metabolism of glyburide. Gliclazide was metabolised by CYP2C. In vitro co-incubation of glyburide or gliclazide with berberine showed that each drug's metabolism was compromised as they share a common isozyme. A strong negative linear correlation of glyburide or gliclazide metabolite levels and the concentration of berberine confirmed the effect of berberine on the metabolism of sulfonylureas.

CONCLUSIONS

The metabolism of sulfonylureas and berberine was affected when these compounds were co-incubated with each other. This may be attributable to competitive binding of the herb and drug to the catalytic sites of the same isozymes.

Molecular Dynamics Simulations Reveal Structural Differences among Allelic Variants of Membrane-Anchored Cytochrome P450 2D6

Cytochrome P450 2D6 (CYP2D6) is an enzyme that is involved in the metabolism of roughly 25% of all marketed drugs and therefore belongs to the most important enzymes in drug metabolism. CYP2D6 features a high degree of genetic polymorphism that can significantly affect the metabolic activity of an individual. In extreme cases, structural changes at the level of single amino acids can either increase its enzymatic activity abolishing the drug therapeutic effect or completely disable the enzyme and elevate drug plasma level potentially leading to adverse effects. In this study, starting from the crystal structure, we built a full-length membrane-anchored all-atom model of the wild-type CYP2D6 as well as five of its variants differing in the enzymatic activity. We validated our models with available experimental data and compared their structural properties with molecular dynamics simulations. The main focus of this study was to identify differences that could mechanistically explain the altered activity of the variants and improve our understanding of their functioning. We observed differences in the opening frequencies and minimal diameters of tunnels that connect the buried active site to the surrounding solvent environment. The variants CYP2D6*4 and CYP2D6*10 associated with missing or decreased activity showed less frequent opening of the tunnels compared to the wild-type. Both CYP2D6*10 and CYP2D6*17 showed a deprivation of an important ligand tunnel suggesting a feasible reason for their altered substrate specificity. Next, the altered fold at the N-terminal anchor region and the decreased active site volume caused by the amino acid mutations of the CYP2D6*4 variant offer an explanation for the absence of its metabolic activity. The mutations in CYP2D6*53 contributed to a significant enlargement of an important ligand tunnel and an extension of the active site cavity. This could explain the altered metabolic profile as well as the enhanced metabolic rates of this particular variant supporting its designation as a possible cause for the ultrarapid metabolizer phenotype. We believe these novel structural insights could advance the fields of personalized medicine and enzyme engineering. Furthermore, they could aid in guiding laboratory as well as computational experiments in the future.

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.

In vitro inhibitory mechanisms and molecular docking of 1'-S-1'-acetoxychavicol acetate on human cytochrome P450 enzymes

BACKGROUND

The compound, 1'-S-1'-acetoxychavicol acetate (ACA), isolated from the rhizomes of a Malaysian ethno-medicinal plant, Alpinia conchigera Griff. (Zingiberaceae), was previously shown to have potential in vivo antitumour activities. In the development of a new drug entity, potential interactions of the compound with the cytochrome P450 superfamily metabolizing enzymes need to be ascertain.

PURPOSE

The concomitant use of therapeutic drugs may cause potential drug-drug interactions by decreasing or increasing plasma levels of the administered drugs, leading to a suboptimal clinical efficacy or a higher risk of toxicity. Thus, evaluating the inhibitory potential of a new chemical entity, and to clarify the mechanism of inhibition and kinetics in the various CYP enzymes is an important step to predict drug-drug interactions.

STUDY DESIGN

This study was designed to assess the potential inhibitory effects of Alpinia conchigera Griff. rhizomes extract and its active constituent, ACA, on nine c-DNA expressed human cytochrome P450s (CYPs) enzymes using fluorescent CYP inhibition assay.

METHODS/RESULTS

The half maximal inhibitory concentration (IC₅₀) of Alpinia conchigera Griff. rhizomes extract and ACA was determined for CYP1A2, CYP2A6, CYP2B6, CYP2C8, CYP2C19, CYP2D6, CYP2E1, CYP3A4 and CYP3A5. A. conchigera extract only moderately inhibits on CYP3A4 (IC₅₀ = 6.76 ± 1.88µg/ml) whereas ACA moderately inhibits the activities of CYP1A2 (IC₅₀ = 4.50 ± 0.10µM), CYP2D6 (IC₅₀ = 7.50 ± 0.17µM) and CYP3A4 (IC₅₀ = 9.50 ± 0.57µM) while other isoenzymes are weakly inhibited. In addition, mechanism-based inhibition studies reveal that CYP1A2 and CYP3A4 exhibited non-mechanism based inhibition whereas CYP2D6 showed mechanism-based inhibition. Lineweaver-Burk plots depict that ACA competitively inhibited both CYP1A2 and CYP3A4, with a K_i values of 2.36 ± 0.03 µM and 5.55 ± 0.06µM, respectively, and mixed inhibition towards CYP2D6 with a K_i value of 4.50 ± 0.08µM. Further, molecular docking studies show that ACA is bound to a few key amino acid residues in the active sites of CYP1A2 and CYP3A4, while one amino residue of CYP2D6 through predominantly Pi-Pi interactions.

CONCLUSION

Overall, ACA may demonstrate drug-drug interactions when co-administered with other therapeutic drugs that are metabolized by CYP1A2, CYP2D6 or CYP3A4 enzymes. Further in vivo studies, however, are needed to evaluate the clinical significance of these interactions.

Characterization of oxycodone in vitro metabolism by human cytochromes P450 and UDP-glucuronosyltransferases

The hepatic metabolism of oxycodone by cytochromes P450 (CYP) and the UDP-glucuronosyltransferases (UGT), the main metabolic enzymes of phase I and phase II, respectively, was assessed in vitro. The N-demethylation by CYP3A4/5 and the O-demethylation by CYP2D6 in human liver microsomes (HLM) followed Michaelis-Menten kinetics, with intrinsic clearances of 1.46μL/min/mg and 0.35μL/min/mg, respectively. Although noroxycodone and oxymorphone mainly contribute to the elimination of oxycodone, the simulated total in vivo clearance using in vitro phase I metabolism was underestimated. For the first time, metabolism of oxycodone by UGT was deeply investigated using HLM, recombinant enzymes and selective inhibitors. Oxycodone-glucuronide was mainly produced by UGT2B7 (K_m=762±153μM, V_max=344±20 peak area/min/mg) and to a lesser extent by UGT2B4 (K_m=2454±497μM, V_max=201±19 peak area/min/mg). Finally, the kinetics of the drug-drug interactions were assessed using two CYP and UGT cocktail approaches. Incubations of HLM with phase I and phase II drug probes showed that oxycodone mainly decreased the in vitro activities of CYP2D6, CYP3A4/5, UGT1A3, UGT1A6 and UGT2B subfamily with an important impact on UGT2B7.

An evaluation of the CYP2D6 and CYP3A4 inhibition potential of metoprolol metabolites and their contribution to drug-drug and drug-herb interaction by LC-ESI/MS/MS

The aim of the present study was to evaluate the contribution of metabolites to drug-drug interaction and drug-herb interaction using the inhibition of CYP2D6 and CYP3A4 by metoprolol (MET) and its metabolites. The peak concentrations of unbound plasma concentration of MET, α-hydroxy metoprolol (HM), O-desmethyl metoprolol (ODM) and N-desisopropyl metoprolol (DIM) were 90.37 ± 2.69, 33.32 ± 1.92, 16.93 ± 1.70 and 7.96 ± 0.94 ng/mL, respectively. The metabolites identified, HM and ODM, had a ratio of metabolic area under the concentration-time curve (AUC) to parent AUC of ≥0.25 when either total or unbound concentration of metabolite was considered. In vitro CYP2D6 and CYP3A4 inhibition by MET, HM and ODM study revealed that MET, HM and ODM were not inhibitors of CYP3A4-catalyzed midazolam metabolism and CYP2D6-catalyzed dextromethorphan metabolism. However, DIM only met the criteria of >10% of the total drug related material and <25% of the parent using unbound concentrations. If CYP inhibition testing is solely based on metabolite exposure, DIM metabolite would probably not be considered. However, the present study has demonstrated that DIM contributes significantly to in vitro drug-drug interaction. Copyright © 2016 John Wiley & Sons, Ltd.

Molecular dynamics and density functional studies on the metabolic selectivity of antipsychotic thioridazine by cytochrome P450 2D6: Connection with crystallographic and metabolic results

CYP2D6, a cytochrome P450 isoform, significantly contributes to the metabolism of many clinically important drugs. Thioridazine (THD) is one of the phenothiazine-type antipsychotics, which exhibit dopamine D2 antagonistic activity. THD shows characteristic metabolic profiles compared to other phenothiazine-type antipsychotics such as chlorpromazine. The sulfur atom attached to the phenothiazine ring is preferentially oxidized mainly by CYP2D6, that is, the 2-sulfoxide is a major metabolite, and interestingly this metabolite shows more potent activity against dopamine D2 receptors than THD. On the other hand, the formation of this metabolite causes many serious problems for its clinical use. Wójcikowski et al. (Drug Metab. Dispos. 2006, 34, 471) reported a kinetic study of THD formed by CYP2D6. Recently, Wang et al. (J. Biol. Chem. 2012, 287, 10834 and J. Biol. Chem. 2015, 290, 5092) revealed the crystallographic structure of THD with CYP2D6. In the current study, the binding and reaction mechanisms at the atomic and electronic levels were computationally examined based on the assumption as to whether or not the different crystallographic binding poses correspond to the different metabolites. The binding and oxidative reaction steps in the whole metabolic process were investigated using molecular dynamics and density functional theory calculations, respectively. The current study demonstrated the essential importance of the orientation of the substrate in the reaction center of CYP2D6 for the metabolic reaction.

Guanfu base A, an antiarrhythmic alkaloid of Aconitum coreanum, Is a CYP2D6 inhibitor of human, monkey, and dog isoforms

Guanfu base A (GFA) is a novel heterocyclic antiarrhythmic drug isolated from Aconitum coreanum (Lèvl.) rapaics and is currently in a phase IV clinical trial in China. However, no study has investigated the influence of GFA on cytochrome P450 (P450) drug metabolism. We characterized the potency and specificity of GFA CYP2D inhibition based on dextromethorphan O-demethylation, a CYP2D6 probe substrate of activity in human, mouse, rat, dog, and monkey liver microsomes. In addition, (+)-bufuralol 1'-hydroxylation was used as a CYP2D6 probe for the recombinant form (rCYP2D6), 2D1 (rCYP2D1), and 2D2 (rCYP2D2) activities. Results show that GFA is a potent noncompetitive inhibitor of CYP2D6, with inhibition constant Ki = 1.20 ± 0.33 μM in human liver microsomes (HLMs) and Ki = 0.37 ± 0.16 μM for the human recombinant form (rCYP2D6). GFA is also a potent competitive inhibitor of CYP2D in monkey (Ki = 0.38 ± 0.12 μM) and dog (Ki = 2.4 ± 1.3 μM) microsomes. However, GFA has no inhibitory activity on mouse or rat CYP2Ds. GFA did not exhibit any inhibition activity on human recombinant CYP1A2, 2A6, 2C8, 2C19, 3A4, or 3A5, but showed slight inhibition of 2B6 and 2E1. Preincubation of HLMs and rCYP2D6 resulted in the inactivation of the enzyme, which was attenuated by GFA or quinidine. Beagle dogs treated intravenously with dextromethorphan (2 mg/ml) after pretreatment with GFA injection showed reduced CYP2D metabolic activity, with the Cmax of dextrorphan being one-third that of the saline-treated group and area under the plasma concentration-time curve half that of the saline-treated group. This study suggests that GFA is a specific CYP2D6 inhibitor that might play a role in CYP2D6 medicated drug-drug interaction.

Quantum mechanical modeling: a tool for the understanding of enzyme reactions

Most enzyme reactions involve formation and cleavage of covalent bonds, while electrostatic effects, as well as dynamics of the active site and surrounding protein regions, may also be crucial. Accordingly, special computational methods are needed to provide an adequate description, which combine quantum mechanics for the reactive region with molecular mechanics and molecular dynamics describing the environment and dynamic effects, respectively. In this review we intend to give an overview to non-specialists on various enzyme models as well as established computational methods and describe applications to some specific cases. For the treatment of various enzyme mechanisms, special approaches are often needed to obtain results, which adequately refer to experimental data. As a result of the spectacular progress in the last two decades, most enzyme reactions can be quite precisely treated by various computational methods.

Enzyme kinetic and molecular docking studies for the inhibitions of miltirone on major human cytochrome P450 isozymes

Previous studies have shown that major tanshinones isolated from Danshen (Salvia miltiorrhiza) inhibited human and rat CYP450 enzymes-mediated metabolism of model probe substrates, with potential in causing herb-drug interactions. Miltirone, another abietane type-diterpene quinone isolated from Danshen, has been reported for its anti-oxidative, anxiolytic and anti-cancer effects. The aim of this study was to study the effect of miltirone on the metabolism of model probe substrates of CYP1A2, 2C9, 2D6 and 3A4 in pooled human liver microsomes. Miltirone showed moderate inhibition on CYP1A2 (IC(50)=1.73 μM) and CYP2C9 (IC(50)=8.61 μM), and weak inhibition on CYP2D6 (IC(50)=30.20 μM) and CYP3A4 (IC(50)=33.88 μM). Enzyme kinetic studies showed that miltirone competitively inhibited CYP2C9 (K(i)=1.48 μM), and displayed mixed type inhibitions on CYP1A2, CYP2D6 and CYP3A4 with K(i) values of 3.17 μM, 24.25 μM and 35.09 μM, respectively. Molecular docking study further confirmed the ligand-binding conformations of miltirone in the active sites of these human CYP450 isoforms, and provided some information on structure-activity relationships for the CYPs inhibition by tanshinones. Taken together, CYPs inhibitions of miltirone were weaker than dihydrotanshinone, but stronger than cryptotanshinone, tanshinone I and tanshinone IIA.

Molecular docking and enzyme kinetic studies of dihydrotanshinone on metabolism of a model CYP2D6 probe substrate in human liver microsomes

The effects of Danshen and its active components (tanshinone I, tanshinone IIA, dihydrotanshinone and cryptotanshinone) on CYP2D6 activity was investigated by measuring the metabolism of a model CYP2D6 probe substrate, dextromethorphan to dextrorphan in human pooled liver microsomes. The ethanolic extract of crude Danshen (6.25-100 μg/ml) decreased dextromethorphan O-demethylation in vitro (IC(50)=23.3 μg/ml) and the water extract of crude Danshen (0.0625-1 mg/ml) showed no inhibition. A commercially available Danshen pill (31.25-500 μg/ml) also decreased CYP2D6 activity (IC(50)=265.8 μg/ml). Among the tanshinones, only dihydrotanshinone significantly inhibited CYP2D6 activity (IC(50)=35.4 μM), compared to quinidine, a specific CYP2D6 inhibitor (IC(50)=0.9 μM). Crytotanshinone, tanshinone I and tanshinone IIA produced weak inhibition, with IC(20) of 40.8 μM, 16.5 μM and 61.4 μM, respectively. Water soluble components such as salvianolic acid B and danshensu did not affect CYP2D6-mediated metabolism. Enzyme kinetics studies showed that inhibition of CYP2D6 activity by the ethanolic extract of crude Danshen and dihydrotanshinone was concentration-dependent, with K(i) values of 4.23 μg/ml and 2.53 μM, respectively, compared to quinidine, K(i)=0.41 μM. Molecular docking study confirmed that dihydrotanshinone and tanshinone I interacted with the Phe120 amino acid residue in the active cavity of CYP2D6 through Pi-Pi interaction, but did not interact with Glu216 and Asp301, the key residues for substrate binding. The logarithm of free binding energy of dihydrotanshinone (-7.6 kcal/mol) to Phe120 was comparable to quinidine (-7.0 kcal/mol) but greater than tanshinone I (-5.4 kcal/mol), indicating dihydrotanshinone has similar affinity to quinidine in binding to the catalytic site on CYP2D6.

Pinpointing P450s associated with pyrethroid metabolism in the dengue vector, Aedes aegypti: developing new tools to combat insecticide resistance

BACKGROUND

Pyrethroids are increasingly used to block the transmission of diseases spread by Aedes aegypti such as dengue and yellow fever. However, insecticide resistance poses a serious threat, thus there is an urgent need to identify the genes and proteins associated with pyrethroid resistance in order to produce effective counter measures. In Ae. aegypti, overexpression of P450s such as the CYP9J32 gene have been linked with pyrethroid resistance. Our aim was to confirm the role of CYP9J32 and other P450s in insecticide metabolism in order to identify potential diagnostic resistance markers.

METHODOLOGY/PRINCIPAL FINDINGS

We have expressed CYP9J32 in Escherichia coli and show that the enzyme can metabolize the pyrethroids permethrin and deltamethrin. In addition, three other Ae. aegypti P450s (CYP9J24, CYP9J26, CYP9J28) were found capable of pyrethroid metabolism, albeit with lower activity. Both Ae. aegypti and Anopheles gambiae P450s (CYP's 6M2, 6Z2, 6P3) were screened against fluorogenic and luminescent substrates to identify potential diagnostic probes for P450 activity. Luciferin-PPXE was preferentially metabolised by the three major pyrethroid metabolisers (CYP9J32, CYP6M2 and CYP6P3), identifying a potential diagnostic substrate for these P450s.

CONCLUSIONS/SIGNIFICANCE

P450s have been identified with the potential to confer pyrethroid resistance in Ae.aegypti. It is recommended that over expression of these enzymes should be monitored as indicators of resistance where pyrethroids are used.

Bioactivities of berberine metabolites after transformation through CYP450 isoenzymes

Background

Berberine (BBR) is a drug with multiple effects on cellular energy metabolism. The present study explored answers to the question of which CYP450 (Cytochrome P450) isoenzymes execute the phase-I transformation for BBR, and what are the bioactivities of its metabolites on energy pathways.

Methods

BBR metabolites were detected using LC-MS/MS. Computer-assistant docking technology as well as bioassays with recombinant CYP450s were employed to identify CYP450 isoenzymes responsible for BBR phase-I transformation. Bioactivities of BBR metabolites in liver cells were examined with real time RT-PCR and kinase phosphorylation assay.

Results

In rat experiments, 4 major metabolites of BBR, berberrubine (M1), thalifendine (M2), demethyleneberberine (M3) and jatrorrhizine (M4) were identified in rat's livers using LC-MS/MS (liquid chromatography-tandem mass spectrometry). In the cell-free transformation reactions, M2 and M3 were detectable after incubating BBR with rCYP450s or human liver microsomes; however, M1 and M4 were below detective level. CYP2D6 and CYP1A2 played a major role in transforming BBR into M2; CYP2D6, CYP1A2 and CYP3A4 were for M3 production. The hepatocyte culture showed that BBR was active in enhancing the expression of insulin receptor (InsR) and low-density-lipoprotein receptor (LDLR) mRNA, as well as in activating AMP-activated protein kinase (AMPK). BBR's metabolites, M1-M4, remained to be active in up-regulating InsR expression with a potency reduced by 50-70%; LDLR mRNA was increased only by M1 or M2 (but not M3 and M4) with an activity level 35% or 26% of that of BBR, respectively. Similarly, AMPK-α phosphorylation was enhanced by M1 and M2 only, with a degree less than that of BBR.

Conclusions

Four major BBR metabolites (M1-M4) were identified after phase-I transformation in rat liver. Cell-free reactions showed that CYP2D6, CYP1A2 and CYP3A4 seemed to be the dominant CYP450 isoenzymes transforming BBR into its metabolites M2 and M3. BBR's metabolites remained to be active on BBR's targets (InsR, LDLR, and AMPK) but with reduced potency.

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

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

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

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.

Polymorphism of human cytochrome P450 enzymes and its clinical impact

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

Computational prediction of drug binding and rationalisation of selectivity towards cytochromes P450

BACKGROUND

Early in-vitro consideration of metabolism and inhibition of cytochrome P450 has proven its merits over the last 15 years. Simultaneously, many computational drug-design methods have been developed, and are being applied to study the interactions between drug candidates and cytochrome P450 enzymes (P450s).

OBJECTIVE

This review discusses the recent advances of these methods and the implications that are specific for P450s.

METHODS

Mainly focusing on the prediction of binding affinity and ligand selectivity, we outline the applicability of the different methods to answer specific questions. Special emphasis is put on the different levels of theory that are being used in recent computational descriptions of ligand-P450 interactions.

CONCLUSION

P450s offer an additional challenge for computational methods, considering the ambiguities of the catalytic cycle and the significant flexibility of the active site. Different computational methods display different limitations, which is crucial to take into account when choosing the method appropriate to each application.

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.

CYP53A15 of Cochliobolus lunatus, a target for natural antifungal compounds

A novel cytochrome P450, CYP53A15, was identified in the pathogenic filamentous ascomycete Cochliobolus lunatus. The protein, classified into the CYP53 family, was capable of para hydroxylation of benzoate. Benzoate is a key intermediate in the metabolism of aromatic compounds in fungi and yet basically toxic to the organism. To guide functional analyses, protein structure was predicted by homology modeling. Since many naturally occurring antifungal phenolic compounds are structurally similar to CYP53A15 substrates, we tested their putative binding into the active site of CYP53A15. Some of these compounds inhibited CYP53A15. Increased antifungal activity was observed when tested in the presence of benzoate. Some results suggest that CYP53A15 O-demethylation activity is important in detoxification of other antifungal substances. With the design of potent inhibitors, CYP53 enzymes could serve as alternative antifungal drug targets.

Functional expression and comparative characterization of nine murine cytochromes P450 by fluorescent inhibition screening

The increasing number of transgenic or gene knockout mouse models generated for use in drug metabolism studies has meant that a greater understanding of the function and substrate specificities of murine cytochromes P450 (P450s) has become essential, particularly with the recent advances in "humanized" mouse models. In this study, we have heterologously expressed nine murine P450s--Cyp1a1, Cyp1a2, Cyp1b1, Cyp2a4, Cyp2b20, Cyp2c29, Cyp2d22, Cyp2e1, and Cyp3a11--individually with human P450 oxidoreductase to generate functional monooxygenase systems in Escherichia coli. We have identified a suitable fluorogenic probe for each P450 and determined the apparent kinetic parameters. These probes have enabled the screening of a panel of 31 test compounds classified as "drugs," "natural compounds," "endogenous compounds," and "pesticides" by measurement of IC(50), thus allowing the comparison of binding affinities. Human P450s CYP2C9, CYP2D6, and CYP3A4 were also included in the study to enable direct comparisons to be made with the mouse enzymes. Although there were general similarities between human and mouse P450s, perhaps the most significant finding in this study was the observation that, despite 77% amino acid identity, Cyp2d22 and CYP2D6 were remarkably dissimilar in a range of enzymatic properties, with potentially serious implications for pharmacokinetic studies using CYP2D substrates. The data presented in this study provide a solid foundation with which to assess the degree of similarity (or difference) between mouse and human P450s involved in xenobiotic metabolism and can be used as a basis for further studies.