Comparison of metabolic soft spot predictions of CYP3A4, CYP2C9 and CYP2D6 substrates using MetaSite and StarDrop

Metabolite identification study plays an important role in determining the sites of metabolic liability of new chemical entities (NCEs) in drug discovery for lead optimization. Here we compare the two predictive software, MetaSite and StarDrop, available for this purpose. They work very differently but are used to predict the site of oxidation by major human cytochrome P450 (CYP) isoforms. Neither software can predict non-CYP catalyzed metabolism nor the rates of metabolism. For the purpose of comparing the two software packages, we tested known probe substrate for these enzymes, which included 12 substrates of CYP3A4 and 18 substrates of CYP2C9 and CYP2D6 were analyzed by each software and the results were compared. It is possible that these known substrates were part of the training set but we are not aware of it. To assess the performance of each software we assigned a point system for each correct prediction. The total points assigned for each CYP isoform experimentally were compared as a percentage of the total points assigned theoretically for the first choice prediction for all substrates for each isoform. Our results show that MetaSite and StarDrop are similar in predicting the correct site of metabolism by CYP3A4 (78% vs 83%, respectively). StarDrop appears to do slightly better in predicting the correct site of metabolism by CYP2C9 and CYP2D6 metabolism (89% and 93%, respectively) compared to MetaSite (63% and 70%, respectively). The sites of metabolism (SOM) from 34 in-house NCEs incubated in human liver microsomes or human hepatocytes were also evaluated using two prediction software packages and the results showed comparable SOM predictions. What makes this comparison challenging is that the contribution of each isoform to the intrinsic clearance (Clint) is not known. Overall the software were comparable except for MetaSite performing better for CYP2D6 and that MetaSite has a liver model that is absent in StarDrop that predicted with 82% accuracy.     

Regioselectivity prediction of CYP1A2-mediated phase I metabolism

A kinetic, reactivity-binding model has been proposed to predict the regioselectivity of substrates meditated by the CYP1A2 enzyme, which is responsible for the metabolism of planar-conjugated compounds such as caffeine. This model consists of a docking simulation for binding energy and a semiempirical molecular orbital calculation for activation energy. Possible binding modes of CYP1A2 substrates were first examined using automated docking based on the crystal structure of CYP1A2, and binding energy was calculated. Then, activation energies for CYP1A2-mediated metabolism reactions were calculated using the semiempirical molecular orbital calculation, AM1. Finally, the metabolic probability obtained from two energy terms, binding and activation energies, was used for predicting the most probable metabolic site. This model predicted 8 out of 12 substrates accurately as the primary preferred site among all possible metabolic sites, and the other four substrates were predicted into the secondary preferred site. This method can be applied for qualitative prediction of drug metabolism mediated by CYP1A2 and other CYP450 family enzymes, helping to develop drugs efficiently.     

In silico modeling of p450 substrates, inhibitors, activators, and inducers

Cytochrome P450 enzymes are the predominant mediators of phase I metabolism of exogenous small molecules. As a result of their extensive role in metabolism of xenobiotics, drug compounds, and endogenous compounds, as well as their wide tissue distribution, significant drug discovery resources are spent to avoid interacting with this class of enzymes. Here we review historical and recent in silico modeling of 7 cytochrome P450 enzymes of particular interest, specifically CYP1A2, CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2D6, and CYP3A4. For each we provide a brief biological background including known inhibitors, substrates, and inducers, as well as details of computational modeling efforts and advances in structural biology. We also provide similar details for 3 nuclear receptors known to regulate gene expression of these enzyme families.     

Capillary electrophoretic investigation of the enantioselective metabolism of propafenone by human cytochrome P-450 SUPERSOMES: Evidence for atypical kinetics by CYP2D6 and CYP3A4

An enantioselective CE method was used to identify the ability of CYP450 enzymes and their stereoselectivity in catalyzing the transformation of propafenone (PPF) to 5-hydroxy-propafenone (5OH-PPF) and N-despropyl-propafenone (NOR-PPF). Using in vitro incubations with single CYP450 enzymes (SUPERSOMES), 5OH-PPF is shown to be selectively produced by CYP2D6 and N-dealkylation is demonstrated to be mediated by CYP2D6, CYP3A4, CYP1A2, and CYP1A1. For the elucidation of kinetic aspects of the metabolism with CYP2D6 and CYP3A4, incubations with individual PPF enantiomers and racemic PPF were investigated. With the exception of the dealkylation in presence of R-PPF only, which can be described by the Michaelis-Menten model, all CYP2D6-induced reactions were found to follow autoactivation kinetics. For CYP3A4, all NOR-PPF enantiomer formation rates as function of PPF enantiomer concentration were determined to follow substrate inhibition kinetics. The formation of NOR-PPF by the different enzymes is stereoselective and is reduced significantly when racemic PPF is incubated. Clearance values obtained for CYP3A4 dealkylation are stereoselective whereas those of CYP2D6 hydroxylation are not. This paper reports the first investigation of the PPF hydroxylation and dealkylation kinetics by the CYP2D6 enzyme and represents the first report in which enantioselective CE data provide the complete in vitro kinetics of metabolic steps of a drug.     

Preparation, characterization, and substrate metabolism of gold-immobilized cytochrome P450 2C9

The cytochrome P450 enzymes represent an important class of heme-containing enzymes. There is considerable interest in immobilizing these enzymes on a surface so that interactions between a single enzyme and other species can be studied with respect to electron transfer, homodimer or heterodimer interactions, or for construction of biological-based chips for standardizing cytochrome P450 metabolism or for high-throughput screening of pharmaceutical agents. Previous studies have generally immobilized P450 enzymes in a matrix or on a surface. Here, we have attached CYP2C9 to gold substrates such that the resulting construct maintains the ability to bind and metabolize substrates in the presence of NADPH and cytochrome P450 reductase. The activity of these chips is directly dependent upon the linkers used to attach CYP2C9 and to the presence of key molecules in the active site during enzyme attachment. A novel method to detect substrate-enzyme binding, namely, superconducting quantum interference device (SQUID) magnetometry, was used to monitor the binding of substrates. Most significantly, conditions that allow measurable CYP2C9 metabolism to occur have been developed.     

Characterization of the CYP3A4 active site by homology modeling

Human microsomal cytochrome P450s participate in drug metabolism and detoxification. Among them, CYP3A4 is the most important isoform for drug-drug interactions. To gain a better understanding of the active site, a homology model of CYP3A4 was constructed based on the crystallographic coordinates of mammalian CYP2C5. The putative active site is much larger than that of CYP2C5 and is divided into three parts (i.e. a proximal and two distal sites from the heme). Most residues reported to be important for ligand-binding are located in the active site of the model. Moreover, some inhibitors (paclitaxel etc.) docked into the model have complementary shapes to the pocket. Pharmacophore docking of 14 substrates was also performed using Ph4Dock of MOE. Calculated interaction energies showed a moderate correlation with the logarithm of apparent K(m) values. These results suggest that this model is reliable enough to be used in the design of compounds for removing undesirable CYP3A4 inhibition.     

Molecular mechanisms of cytochrome p450 induction: potential for drug-drug interactions

The human liver cytochromes P450 (CYP P450s) are superfamily of hemoproteins responsible for catalyzing the oxidative metabolism of drugs and xenobiotics entering human body. Drug-drug/xenobiotic interactions are a major cause of therapeutic failures and adverse events. The concomitant administration of inducers with other drugs that are metabolized by CYP450 can result in their altered metabolism in the gastrointestinal tract and/ or liver. The clinical importance of such interactions includes auto induction leading to suboptimal or failed treatment. It is a major concern for the drug companies while developing new drugs. The present understanding of the mechanisms of induction of CYP P450s enzymes and their regulation has made considerable progress during last few years. However there are still gaps in our understanding on molecular aspects of CYP enzymes. Therefore, it remains the subject of intense scientific research to ascertain their in vivo function, and also better understand how the expression of CYP enzymes is regulated at the molecular level. This review analyzes and presents recent findings and concepts on xenosensors and their target genes. Emphasis is given to the molecular mechanisms and signaling pathways of CYP P450 mediated induction by xenobiotics and their potential for drug-drug interactions.     

基于液相色谱-质谱联用的肝微粒体中药物代谢酶定量方法研究进展

肝微粒体中的药物代谢酶是参与药物代谢的主要酶类,不管是合成药物,还是组成复杂的中药,主要经过肝脏中药物代谢酶的代谢转化成易于排泄的化合物。由于药物代谢酶的表达水平具有较大的个体差异性,因此,药物代谢酶的准确定量对药物的药理、药物相互作用和临床应用具有重要意义。本文综述了近十年来药物代谢酶定量方法的研究进展。     

Kinetic Characteristics of Curcumin and Germacrone in Rat and Human Liver Microsomes: Involvement of CYP Enzymes

Curcumin and germacrone, natural products present in the Zingiberaceae family of plants, have several biological properties. Among these properties, the anti-NSCLC cancer action is noteworthy. In this paper, kinetics of the two compounds in rat liver microsomes (RLMs), human liver microsomes (HLMs), and cytochrome P450 (CYP) enzymes (CYP3A4, 1A2, 2E1, and 2C19) in an NADPH-generating system in vitro were evaluated by UP-HPLC–MS/MS (ultrahigh-pressure liquid chromatography–tandem mass spectrometry). The contents of four cytochrome P450 (CYP) enzymes, adjusting by the compounds were detected using Western blotting in vitro and in vivo. The t_1/2 of curcumin was 22.35 min in RLMs and 173.28 min in HLMs, while 18.02 and 16.37 min were gained for germacrone. The V_max of curcumin in RLMs was about 4-fold in HLMs, meanwhile, the V_max of germacrone in RLMs was similar to that of HLMs. The single enzyme t1/2 of curcumin was 38.51 min in CYP3A4, 301.4 min in 1A2, 69.31 min in 2E1, 63.01 min in 2C19; besides, as to the same enzymes, t1/2 of germacrone was 36.48 min, 86.64 min, 69.31 min, and 57.76 min. The dynamic curves were obtained by reasonable experimental design and the metabolism of curcumin and germacrone were selected in RLMs/HLMs. The selectivities in the two liver microsomes differed in degradation performance. These results meant that we should pay more attention to drugs in clinical medication–drug and drug–enzyme interactions.     

High-throughput screening of inhibitory potential of nine cytochrome P450 enzymes in vitro using liquid chromatography/tandem mass spectrometry

The early detection of potential drug-drug interactions is an important issue of drug discovery that has led to the development of high-throughput screening (HTS) methods for potential drug interactions. We developed a HTS method for potential interactions of inhibitory drugs for nine human P450 enzymes using cocktail incubation and tandem mass spectrometry in vitro. This new method involves incubation of two cocktail doses and single cassette analysis. The two cocktail doses in vitro were developed to minimize solvent effects and mutual drug interactions among substrates: cocktail A was composed of phenacetin for CYP1A2, coumarin for CYP2A6, paclitaxel for CYP2C8, S-mephenytoin for CYP2C19, dextromethorphan for CYP2D6, and midazolam for CYP3A4; and cocktail B was composed of three substrates including bupropion for CYP2B6, tolbutamide for CYP2C9, and chlorzoxazone for CYP2E1. In the incubation study of these cocktails, the reaction mixtures were pooled and simultaneously analyzed using liquid chromatography/tandem mass spectrometry employing a fast gradient. The method was validated by comparing the inhibition data obtained from the incubation of each individual probe substrate alone with data from the new method. The IC50 value of each inhibitor in the cocktail agreed well with that of the individual probe drug as well as with values previously reported in the literature. As a HTS method for potential interactions of the inhibition of these nine P450 enzymes, this new method will be useful in the drug discovery process and for the mechanistic understanding of drug interactions.     

Cytochrome P450 reaction phenotyping and inhibition and induction studies of pinostrobin in human liver microsomes and hepatocytes

Pinostrobin (PI, 5‐hydroxy‐7‐methoxyflavanone) is a natural flavonoid known for its rich pharmacological activities. The objective of this study was to identify the human liver cytochrome P450 (CYP450) isoenzymes involved in the metabolism of PI. A single hydoxylated metabolite was obtained from PI after an incubation with pooled human liver microsomes (HLMs). The relative contributions of different CYP450s were evaluated using CYP450‐selective inhibitors in HLMs and recombinant human CYP450 enzymes, and the results revealed the major involvement of CYP1A2, CYP2C9 and CYP2E1 in PI metabolism. We also evaluated the ability of PI to inhibit and induce human cytochrome P450 enzymes in vitro . High‐performance liquid chromatography and liquid chromatography–tandem mass spectrometry analytical techniques were used to estimate the enzymatic activities of seven drug‐metabolizing CYP450 isozymes in vitro . In HLMs, PI did not inhibit CYP 1A2, CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2D6 or CYP3A4 (IC₅₀ > 100 μm ). In the induction studies, PI had minimal effects on CYP1A2, CYP2B6and CYP3A4 activity. Based on these results, PI would not be expected to cause clinically significant CYP450 inhibition or induction.     

Homology modeling, binding site identification and docking in flavone hydroxylase CYP105P2 in Streptomyces peucetius ATCC 27952

Homology models of cytochrome P450 105P2 (CYP105P2) were constructed using four P450 structures, CYP105A1, CYP105, CYP165B3 and CYP107L1, as templates for the model building. Using Accelrys Discovery Studio 2.1 software, the lowest energy CYP105P2 model was then assessed for stereochemical quality and side-chain environment. Further active site optimization of the CYP105P2 model built using these templates was performed by molecular dynamics to generate the final CYP105P2 model. The substrates, flavone, flavanone, quercetin and naringenin, were docked into the model. The model-flavone complex was used to validate the active site architecture, and structurally and functionally important residues were identified by subsequent characterization of the secondary structure.     

Novel hierarchical classification and visualization method for multiobjective optimization of drug properties: application to structure-activity relationship analysis of cytochrome P450 metabolism

In the lead optimization process, medicinal chemists must consider various chemical properties of active compounds, including ADME/Tox properties, and find the best compromise among these. This study presents a novel data mining method for multiobjective optimization of chemical properties, which consists of the hierarchical classification and visualization of multidimensional data. A hierarchical classification tree model is generated by an extension of recursive partitioning that utilizes averaged information gains for multiple objective variables as a quality-of-split criterion. All the hierarchically structured data objects are represented using a large-scale data visualization technique. The technique is an extension of HeiankyoView, which displays data objects as colored icons and group nodes as rectangular borders. Each icon is divided into subregions with different colors, so that it can present multidimensional data according to brightness of the colors. The proposed method was applied to the structure-activity relationship analysis for cytochrome P450 (CYP) substrates. The substrate specificity of six CYP isoforms was successfully delineated: e.g., CYP2C9 substrates are anionic compounds, while CYP2D6 substrates are cationic; and CYP2E1 substrates are smaller compounds, while CYP3A4 substrates are larger compounds.