Structural characterization of human cytochrome P450 2C19: active site differences between P450s 2C8, 2C9, and 2C19

Enzyme Kinetics and Molecular Docking Studies on Cytochrome 2B6, 2C19, 2E1, and 3A4 Activities by Sauchinone

Sauchinone, an active lignan isolated from the aerial parts of Saururus chinensis (Saururaceae), exhibits anti-inflammatory, anti-obesity, anti-hyperglycemic, and anti-hepatic steatosis effects. As herb–drug interaction (HDI) through cytochrome P450s (CYPs)-mediated metabolism limits clinical application of herbs and drugs in combination, this study sought to explore the enzyme kinetics of sauchinone towards CYP inhibition in in vitro human liver microsomes (HLMs) and in vivo mice studies and computational molecular docking analysis. In in vitro HLMs, sauchinone reversibly inhibited CYP2B6, 2C19, 2E1, and 3A4 activities in non-competitive modes, showing inhibition constant (K_i) values of 14.3, 16.8, 41.7, and 6.84 μM, respectively. Also, sauchinone time-dependently inhibited CYP2B6, 2E1 and 3A4 activities in vitro HLMs. Molecular docking study showed that sauchinone could be bound to a few key amino acid residues in the active site of CYP2B6, 2C19, 2E1, and 3A4. When sibutramine, clopidogrel, or chlorzoxazone was co-administered with sauchinone to mice, the systemic exposure of each drug was increased compared to that without sauchinone, because sauchinone reduced the metabolic clearance of each drug. In conclusion, when sauchinone was co-treated with drugs metabolized via CYP2B6, 2C19, 2E1, or 3A4, sauchinone–drug interactions occurred because sauchinone inhibited the CYP-mediated metabolic activities.

Identification of karanjin isolated from the Indian beech tree as a potent CYP1 enzyme inhibitor with cellular efficacy via screening of a natural product repository

CYP1A1 is thought to mediate carcinogenesis in oral, lung and epithelial cancers. In order to identify a CYP1A1 inhibitor from an edible plant, 394 natural products in the IIIM's natural product repository were screened, at 10 μM concentration, using CYP1A1-Sacchrosomes™ (i.e. microsomal enzyme isolated from recombinant baker's yeast). Twenty-seven natural products were identified that inhibited 40-97% of CYP1A1's 7-ethoxyresorufin-O-deethylase activity. The IC50 values of the 'hits', belonging to different chemical scaffolds, were determined. Their selectivity was studied against a panel of 8 CYP-Sacchrosomes™. In order to assess cellular efficacy, the 'hits' were screened for their capability to inhibit CYP enzymes expressed within live recombinant human embryonic kidney (HEK293) cells from plasmids encoding specific CYP genes (1A2, 1B1, 2C9, 2C19, 2D6, 3A4). Isopimpinellin (IN-475; IC50, 20 nM) and karanjin (IN-195; IC50, 30 nM) showed the most potent inhibition of CYP1A1 in human cells. Isopimpinellin is found in celery, parsnip, fruits and in the rind and pulp of limes whereas different parts of the Indian beech tree, which contain karanjin, have been used in traditional medicine. Both isopimpinellin and karanjin negate the cellular toxicity of CYP1A1-mediated benzo[a]pyrene. Molecular docking and molecular dynamic simulations with CYP isoforms rationalize the observed trends in the potency and selectivity of isopimpinellin and karanjin.

Molecular dynamics investigations of regioselectivity of anionic/aromatic substrates by a family of enzymes: a case study of diclofenac binding in CYP2C isoforms

The CYP2C subfamily is of particular importance in the metabolism of drugs, food toxins, and procarcinogens. Like other P450 subfamilies, 2C enzymes share a high sequence identity, but significantly contribute in different ways to hepatic capacity to metabolize drugs. They often metabolize the same substrate to more than one product with different catalytic sites. Because it is challenging to characterize experimentally, much still remains unknown about the reason for why the substrate regioselectivity of these closely related subfamily members is different. Here, we have investigated the structural features of CYP2C8, CYP2C9, and CYP2C19 bound with their shared substrate diclofenac to elucidate the underlying molecular mechanism for the substrate regioselectivity of CYP2C subfamily enzymes. The obtained results demonstrate how a sequence divergence for the active site residues causes heterogeneous variations in the secondary structures and in major tunnel selections, and further affects the shape and chemical properties of the substrate-binding site. Structural analysis and free energy calculations showed that the most important determinants of regioselectivity among the CYP2C isoforms are the geometrical features of the active sites, as well as the hydrogen bonds and the hydrophobic interactions, mainly presenting as the various locations of Arg108 and substitutions of Phe205 for Ile205 in CYP2C8. The MM-GB/SA calculations combined with PMF results accord well with the experimental KM values, bridging the gap between the theory and the experimentally observed results of binding affinity differences. The present study provides important insights into the structure-function relationships of CYP2C subfamily enzymes, the knowledge of ligand binding characteristics and key residue contributions could guide future experimental and computational work on the synthesis of drugs with better pharmacokinetic properties so that CYP interactions could be avoided.

Pharmacogenomics of CYP2C9: Functional and Clinical Considerations

CYP2C9 is the most abundant CYP2C subfamily enzyme in human liver and the most important contributor from this subfamily to drug metabolism. Polymorphisms resulting in decreased enzyme activity are common in the CYP2C9 gene and this, combined with narrow therapeutic indices for several key drug substrates, results in some important issues relating to drug safety and efficacy. CYP2C9 substrate selectivity is detailed and, based on crystal structures for the enzyme, we describe how CYP2C9 catalyzes these reactions. Factors relevant to clinical response to CYP2C9 substrates including inhibition, induction and genetic polymorphism are discussed in detail. In particular, we consider the issue of ethnic variation in pattern and frequency of genetic polymorphisms and clinical implications. Warfarin is the most well studied CYP2C9 substrate; recent work on use of dosing algorithms that include CYP2C9 genotype to improve patient safety during initiation of warfarin dosing are reviewed and prospects for their clinical implementation considered. Finally, we discuss a novel approach to cataloging the functional capabilities of rare 'variants of uncertain significance', which are increasingly detected as more exome and genome sequencing of diverse populations is conducted.

Identification of novel natural inhibitors of Opisthorchis felineus cytochrome P450 using structure-based screening and molecular dynamic simulation

Opisthorchis felineus is the etiological agent of opisthorchiasis in humans. O. felineus cytochrome P450 (OfCYP450) is an important enzyme in the parasite xenobiotic metabolism. To identify the potential anti-opisthorchid compound, we conducted a structure-based virtual screening of natural compounds from the ZINC database (n = 1,65,869) against the OfCYP450. The ligands were screened against OfCYP450 in four sequential docking modes that resulted in 361 ligands having better docking score. These compounds were evaluated for Lipinski and ADMET prediction, and 10 compounds were found to fit well with re-docking studies. After refinement by docking and drug-likeness analyses, four potential inhibitors (ZINC2358298, ZINC8790946, ZINC70707116, and ZINC85878789) were identified. These ligands with reference compounds (itraconazole and fluconazole) were further subjected to molecular dynamics simulation (MDS) and binding energy analyses to compare the dynamic structure of protein after ligand binding and the stability of the OfCYP450 and bound complexes. The binding energy analyses were also calculated. The results suggested that the compounds had a negative binding energy with -259.41, -110.09, -188.25, -163.30, -202.10, and -158.79 kJ mol^-1 for itraconazole, fluconazole, and compounds with IDs ZINC2358298, ZINC8790946, ZINC70707116, and ZINC85878789, respectively. These lead compounds displayed significant pharmacological and structural properties to be drug candidates. On the basis of MDS results and binding energy analyses, we concluded that ZINC8790946, ZINC70707116, and ZINC85878789 have excellent potential to inhibit OfCYP450.

Interactions of omeprazole-based analogues with cytochrome P450 2C19: a computational study

Cytochrome P450 2C19 (CYP2C19) is one of 57 drug metabolizing enzymes in humans and is responsible for the metabolism of ∼7-10% of drugs in clinical use. Recently omeprazole-based analogues were reported to be the potent inhibitors of CYP2C19 and have the potential to be used as the tool compounds for studying the substrate selectivity of CYP2C19. However, the binding modes of these compounds with CYP2C19 remain to be elucidated. In this study, a combination of molecular docking, molecular dynamics (MD), and MM/GBSA calculations was employed to systematically investigate the interactions between these compounds and CYP2C19. The binding modes of these analogues were analyzed in detail. The results indicated that the inclusion of explicit active site water molecules could improve binding energy prediction when the water molecules formed a hydrogen bonding network between the ligand and protein. We also found that the effect of active site water molecules on binding free energy prediction was dependent on the ligand binding modes. Our results unravel the interactions of these omeprazole-based analogues with CYP2C19 and might be helpful for the future design of potent CYP2C19 inhibitors with improved metabolic properties.

Identification of Potent and Selective CYP1A1 Inhibitors via Combined Ligand and Structure-Based Virtual Screening and Their in Vitro Validation in Sacchrosomes and Live Human Cells

Target structure-guided virtual screening (VS) is a versatile, powerful, and inexpensive alternative to experimental high-throughput screening (HTS). To discover potent CYP1A1 enzyme inhibitors for cancer chemoprevention, a commercial library of 50 000 small molecules was utilized for VS guided by both ligand and structure-based strategies. For experimental validation, 300 ligands were proposed based on combined analysis of fitness scores from ligand based e-pharmacophore screening and docking score, prime MMGB/SA binding affinity and interaction pattern analysis from structure-based VS. These 300 compounds were screened, at 10 μM concentration, for in vitro inhibition of CYP1A1-Sacchrosomes (yeast-derived microsomal enzyme) in the ethoxyresorufin-O-de-ethylase assay. Thirty-two compounds displayed >50% inhibition of CYP1A1 enzyme activity at 10 μM. 2-Phenylimidazo-[1,2-a]quinoline (5121780, 119) was found to be the most potent with 97% inhibition. It also inhibited ∼95% activity of CYP1B1 and CYP1A2, the other two CYP1 enzymes. The compound 5121780 (119) showed high selectivity toward inhibition of CYP1 enzymes with respect to CYP2 and CYP3 enzymes (i.e., there was no detectable inhibition of CYP2D6/CYP2C9/CYP2C19 and CYP3A4 at 10 μM). It was further investigated in live CYP-expressing human cell system, which confirmed that compound 5121780 (119) potently inhibited CYP1A1, CYP1A2, CYP1B1 enzymes with IC₅₀ values of 269, 30, and 56 nM, respectively. Like in Sacchrosomes, inhibition of CYP2D6/CYP2C9/CYP2C19 and CYP3A4 enzymes, expressed within live human cells, could hardly be detected at 10 μM. The compound 119 rescued CYP1A1 overexpressing HEK293 cells from CYP1A1 mediated benzo[a]pyrene (B[a]P) toxicity and also overcame cisplatin resistance in CYP1B1 overexpressing HEK293 cells. Molecular dynamics simulations of 5121780 (119) with CYP1 enzymes was performed to understand the interaction pattern to CYP isoforms. Results indicate that VS can successfully be used to identify promising CYP1A1 inhibitors, which may have potential in the development of novel cancer chemo-preventive agents.

Molecular dynamics investigations of membrane-bound CYP2C19 polymorphisms reveal distinct mechanisms for peripheral variants by long-range effects on the enzymatic activity

The results reveal distinct mechanisms for enzymatic activity deficiencies upon two peripheral variants in CYP2C19.

Functioning of drug-metabolizing microsomal cytochrome P450s: In silico probing of proteins suggests that the distal heme 'active site' pocket plays a relatively 'passive role' in some enzyme-substrate interactions

PURPOSE

The currently held mechanistic understanding of microsomal cytochrome P450s (CYPs) seeks that diverse drug molecules bind within the deep-seated distal heme pocket and subsequently react at the heme centre. To explain a bevy of experimental observations and meta-analyses, we indulge a hypothesis that involves a "diffusible radical mediated" mechanism. This new hypothesis posits that many substrates could also bind at alternate loci on/within the enzyme and be reacted without the pertinent moiety accessing a bonding proximity to the purported catalytic Fe-O enzyme intermediate.

METHODS

Through blind and heme-distal pocket centered dockings of various substrates and non-substrates (drug molecules of diverse sizes, classes, topographies etc.) of microsomal CYPs, we explored the possibility of access of substrates via the distal channels, its binding energies, docking orientations, distance of reactive moieties (or molecule per se) to/from the heme centre, etc. We investigated specific cases like- (a) large drug molecules as substrates, (b) classical marker drug substrates, (c) class of drugs as substrates (Sartans, Statins etc.), (d) substrate preferences between related and unrelated CYPs, (e) man-made site-directed mutants' and naturally occurring mutants' reactivity and metabolic disposition, (f) drug-drug interactions, (g) overall affinities of drug substrate versus oxidized product, (h) meta-analysis of in silico versus experimental binding constants and reaction/residence times etc.

RESULTS

It was found that heme-centered dockings of the substrate/modulator drug molecules with the available CYP crystal structures gave poor docking geometries and distances from Fe-heme centre. In conjunction with several other arguments, the findings discount the relevance of erstwhile hypothesis in many CYP systems. Consequently, the newly proposed hypothesis is deemed a viable alternate, as it satisfies Occam's razor.

CONCLUSIONS

The new proposal affords expanded scope for explaining the mechanism, kinetics and overall phenomenology of CYP mediated drug metabolism. It is now understood that the heme-iron and the hydrophobic distal pocket of CYPs serve primarily to stabilize the reactive intermediate (diffusible radical) and the surface or crypts of the apoprotein bind to the xenobiotic substrate (and in some cases, the heme distal pocket could also serve the latter function). Thus, CYPs enhance reaction rates and selectivity/specificity via a hitherto unrecognized modality.

Pharmacogenomics of the cytochrome P450 2C family: impacts of amino acid variations on drug metabolism

Pharmacogenomics investigates DNA and RNA variations in the human genome related to drug responses. Cytochrome P450 (CYP) is a supergene family of drug-metabolizing enzymes responsible for the metabolism of approximately 90% of human drugs. Among the major CYP isoforms, the CYP2C subfamily is of clinical significance because it metabolizes approximately 20% of clinically administrated drugs and represents several variant alleles leading to adverse drug reactions or altering drug efficacy. Here, we review recent progress on understanding the interindividual variability of the CYP2C members and the functional and clinical impact on drug metabolism. We summarize current advances in the molecular modeling of CYP2C polymorphisms and discuss the structural bases and molecular mechanisms of amino acid variants of CYP2C members that affect drug metabolism.

Hepatic cytochromes P450: structural degrons and barcodes, posttranslational modifications and cellular adapters in the ERAD-endgame

The endoplasmic reticulum (ER)-anchored hepatic cytochromes P450 (P450s) are enzymes that metabolize endo- and xenobiotics i.e. drugs, carcinogens, toxins, natural and chemical products. These agents modulate liver P450 content through increased synthesis or reduction via inactivation and/or proteolytic degradation, resulting in clinically significant drug-drug interactions. P450 proteolytic degradation occurs via ER-associated degradation (ERAD) involving either of two distinct routes: Ubiquitin (Ub)-dependent 26S proteasomal degradation (ERAD/UPD) or autophagic lysosomal degradation (ERAD/ALD). CYP3A4, the major human liver/intestinal P450, and the fast-turnover CYP2E1 species are degraded via ERAD/UPD entailing multisite protein phosphorylation and subsequent ubiquitination by gp78 and CHIP E3 Ub-ligases. We are gaining insight into the nature of the structural determinants involved in CYP3A4 and CYP2E1 molecular recognition in ERAD/UPD [i.e. K48-linked polyUb chains and linear and/or "conformational" phosphodegrons consisting either of consecutive sequences on surface loops and/or disordered regions, or structurally-assembled surface clusters of negatively charged acidic (Asp/Glu) and phosphorylated (Ser/Thr) residues, within or vicinal to which, Lys-residues are targeted for ubiquitination]. Structural inspection of select human liver P450s reveals that such linear or conformational phosphodegrons may indeed be a common P450-ERAD/UPD feature. By contrast, although many P450s such as the slow-turnover CYP2E1 species and rat liver CYP2B1 and CYP2C11 are degraded via ERAD/ALD, little is known about the mechanism of their ALD-targeting. On the basis of our current knowledge of ALD-substrate targeting, we propose a tripartite conjunction of K63-linked Ub-chains, P450 structural "LIR" motifs and selective cellular "cargo receptors" as plausible P450-ALD determinants.

Biotransformation of 8:2 fluorotelomer alcohol by recombinant human cytochrome P450s, human liver microsomes and human liver cytosol

Biotransformation of 8:2 fluorotelomer alcohol (8:2 FTOH) can form potentially more toxic metabolites. However, the responsible cytochrome P450 (CYP) isoform(s) and phase II metabolism have not been studied in humans. Here, we characterized the in vitro metabolism of 8:2 FTOH by recombinant human CYPs, human liver microsomes, and human liver cytosol. The results showed that among the 11 isoforms investigated, CYP2C19 was the only enzyme capable of catalyzing 8:2 FTOH with Km and Vmax values of 18.8 μM and 8.52 pmol min(-1) pmol(-1) P450, respectively. The phase I metabolite was identified as 8:2 fluorotelomer aldehyde (8:2 FTAL). HLMs also catalyzed 8:2 FTOH transformation, with the Vmax and intrinsic clearance (CLint) values similar to those of CYP2C19 after the protein content is taken into account. Molecular docking showed that the hydroxyl group of 8:2 FTOH accesses the heme iron-oxo of CYP2C19 in an energetically favored orientation. 8:2 FTOH was also transformed by phase II enzymes to form O-glucuronide and O-sulfate conjugates. The CLint values follow the order of sulfation > oxidation > glucuronidation, suggesting that conjugation is the major metabolic pathway, which explains the low yield of perfluoroalkyl acids (PFCAs). These results provide new insight into fluorotelomer alcohol biotransformation and indirect human exposure to PFCAs.

Repurposing Resveratrol and Fluconazole To Modulate Human Cytochrome P450-Mediated Arachidonic Acid Metabolism

Cytochrome P450 (P450) enzymes metabolize arachidonic acid (AA) to several biologically active epoxyeicosatrienoic acids (EETs) and hydroxyeicosatetraenoic acids (HETEs). Repurposing clinically-approved drugs could provide safe and readily available means to control EETs and HETEs levels in humans. Our aim was to determine how to significantly and selectively modulate P450-AA metabolism in humans by clinically-approved drugs. Liquid chromatography-mass spectrometry was used to determine the formation of 15 AA metabolites by human recombinant P450 enzymes, as well as human liver and kidney microsomes. CYP2C19 showed the highest EET-forming activity, while CYP1B1 and CYP2C8 showed the highest midchain HETE-forming activities. CYP1A1 and CYP4 showed the highest subterminal- and 20-HETE-forming activity, respectively. Resveratrol and fluconazole produced the most selective and significant modulation of hepatic P450-AA metabolism, comparable to investigational agents. Monte Carlo simulations showed that 90% of human population would experience a decrease by 6-22%, 16-39%, and 16-35% in 16-, 18-, and 20-HETE formation, respectively, after 2.5 g daily of resveratrol, and by 22-31% and 14-23% in 8,9- and 14,15-EET formation after 50 mg of fluconazole. In conclusion, clinically-approved drugs can provide selective and effective means to modulate P450-AA metabolism, comparable to investigational drugs. Resveratrol and fluconazole are good candidates to be repurposed as new P450-based treatments.

Effects of protein flexibility and active site water molecules on the prediction of sites of metabolism for cytochrome P450 2C19 substrates

Structure-based prediction of sites of metabolism (SOMs) mediated by cytochrome P450s (CYPs) is of great interest in drug discovery and development.

Toxicity Originating from Thiophene Containing Drugs: Exploring the Mechanism using Quantum Chemical Methods

Drug metabolism of thiophene containing substrates by cytochrome P450s (CYP450) leads to toxic side effects, for example, nephrotoxicity (suprofen, ticlopidine), hepatotoxicity (tienilic acid), thrombotic thrombocytopenic purpura (clopidogrel), and aplastic anemia (ticlopidine). The origin of toxicity in these cases has been attributed to two different CYP450 mediated metabolic reactions: S-oxidation and epoxidation. In this work, the molecular level details of the bioinorganic chemistry associated with the generation of these competitive reactions are reported. Density functional theory was utilized (i) to explore the molecular mechanism for S-oxidation and epoxidation using the radical cationic center Cpd I [(iron(IV)-oxo-heme porphine system with SH(-) as the axial ligand, to mimic CYP450s] as the model oxidant, (ii) to establish the 3D structures of the reactants, transition states, and products on both the metabolic pathways, and (iii) to examine the potential energy (PE) profile for both the pathways to determine the energetically preferred toxic metabolite formation. The energy barrier required for S-oxidation was observed to be 14.75 kcal/mol as compared to that of the epoxidation reaction (13.23 kcal/mol) on the doublet PE surface of Cpd I. The formation of the epoxide metabolite was found to be highly exothermic (-23.24 kcal/mol), as compared to S-oxidation (-8.08 kcal/mol). Hence, on a relative scale the epoxidation process was observed to be thermodynamically and kinetically more favorable. The energy profiles associated with the reactions of the S-oxide and epoxide toxic metabolites were also explored. This study helps in understanding the CYP450-catalyzed toxic reactions of drugs containing the thiophene ring at the atomic level.

Identification and in silico prediction of metabolites of the model compound, tebufenozide by human CYP3A4 and CYP2C19

The metabolites of tebufenozide, a model compound, formed by the yeast-expressed human CYP3A4 and CYP2C19 were identified to clarify the substrate recognition mechanism of the human cytochrome P450 (CYP) isozymes. We then determined whether tebufenozide metabolites may be predicted in silico. Hydrogen abstraction energies were calculated with the density functional theory method B3LYP/6-31G(∗). A docking simulation was performed using FRED software. Several alkyl sites of tebufenozide were hydroxylated by CYP3A4 whereas only one site was modified by CYP2C19. The accessibility of each site of tebufenozide to the reaction center of CYP enzymes and the susceptibility of each hydrogen atom for metabolism by CYP enzymes were evaluated by a docking simulation and hydrogen abstraction energy estimation, respectively.

Time-Dependent Inhibition of CYP2C19 by Isoquinoline Alkaloids: In Vitro and In Silico Analysis

The cytochrome P450 2C19 (CYP2C19) enzyme plays an important role in the metabolism of many commonly used drugs. Relatively little is known about CYP2C19 inhibitors, including compounds of natural origin, which could inhibit CYP2C19, potentially causing clinically relevant metabolism-based drug interactions. We evaluated a series (N = 49) of structurally related plant isoquinoline alkaloids for their abilities to interact with CYP2C19 enzyme using in vitro and in silico methods. We examined several common active alkaloids found in herbal products such as apomorphine, berberine, noscapine, and papaverine, as well as the previously identified mechanism-based inactivators bulbocapnine, canadine, and protopine. The IC50 values of the alkaloids ranged from 0.11 to 210 µM, and 42 of the alkaloids were confirmed to be time-dependent inhibitors of CYP2C19. Molecular docking and three-dimensional quantitative structure-activity relationship analysis revealed key interactions of the potent inhibitors with the enzyme active site. We constructed a comparative molecular field analysis model that was able to predict the inhibitory potency of a series of independent test molecules. This study revealed that many of these isoquinoline alkaloids do have the potential to cause clinically relevant drug interactions. These results highlight the need for studying more profoundly the potential interactions between drugs and herbal products. When further refined, in silico methods can be useful in the high-throughput prediction of P450 inhibitory potential of pharmaceutical compounds.

In vitro functional analysis of 24 novel CYP2C19 variants recently found in the Chinese Han population

1. CYP2C19 is a highly polymorphic enzyme responsible for the metabolism of a wide range of clinical drugs. Alterations to the CYP2C19 gene contribute to the variability of CYP2C19 enzyme activity, which causes pharmacokinetics and drug efficacies to vary and adverse drug reactions to occur in different persons. Recently, we identified 24 novel CYP2C19 allelic variants in the Chinese Han population. The purpose of present study is to assess the impact of these newly found nucleotide mutations on the enzymatic activity of the CYP2C19 protein. 2. Dual-expression vectors were constructed and transiently transfected into 293FT cells. Forty-eight hours after transfection, cells were re-suspended and incubated with two typical probe substrates, omeprazole and S-mephenytoin, to determine the activities of each variant relative to the wild-type protein. 3. Immunoblotting results showed that the protein expression levels of the CYP2C19 variants were diverse. Enzymatic ability analysis showed that the variant 35FS exhibited no functional activity, and most of the other variants showed significantly decreased metabolic activities toward both omeprazole and S-mephenytoin compared with wild-type. 4. These findings greatly enrich the knowledge of biological effects of these newly found CYP2C19 mutations and aid the application of this knowledge to future individualized drug therapy in clinic.

Prediction of cytochrome P450 mediated metabolism

Cytochrome P450 enzymes (CYPs) form one of the most important enzyme families involved in the metabolism of xenobiotics. CYPs comprise many isoforms, which catalyze a wide variety of reactions, and potentially, a large number of different metabolites can be formed. However, it is often hard to rationalize what metabolites these enzymes generate. In recent years, many different in silico approaches have been developed to predict binding or regioselective product formation for the different CYP isoforms. These comprise ligand-based methods that are trained on experimental CYP data and structure-based methods that consider how the substrate is oriented in the active site or/and how reactive the part of the substrate that is accessible to the heme group is. We will review key aspects for various approaches that are available to predict binding and site of metabolism (SOM), what outcome can be expected from the predictions, and how they could potentially be improved.

Modeling of interactions between xenobiotics and cytochrome P450 (CYP) enzymes

The adverse effects to humans and environment of only few chemicals are well known. Absorption, distribution, metabolism, and excretion (ADME) are the steps of pharmaco/toxicokinetics that determine the internal dose of chemicals to which the organism is exposed. Of all the xenobiotic-metabolizing enzymes, the cytochrome P450 (CYP) enzymes are the most important due to their abundance and versatility. Reactions catalyzed by CYPs usually turn xenobiotics to harmless and excretable metabolites, but sometimes an innocuous xenobiotic is transformed into a toxic metabolite. Data on ADME and toxicity properties of compounds are increasingly generated using in vitro and modeling (in silico) tools. Both physics-based and empirical modeling approaches are used. Numerous ligand-based and target-based as well as combined modeling methods have been employed to evaluate determinants of CYP ligand binding as well as predicting sites of metabolism and inhibition characteristics of test molecules. In silico prediction of CYP–ligand interactions have made crucial contributions in understanding (1) determinants of CYP ligand binding recognition and affinity; (2) prediction of likely metabolites from substrates; (3) prediction of inhibitors and their inhibition potency. Truly predictive models of toxic outcomes cannot be created without incorporating metabolic characteristics; in silico methods help producing such information and filling gaps in experimentally derived data. Currently modeling methods are not mature enough to replace standard in vitro and in vivo approaches, but they are already used as an important component in risk assessment of drugs and other chemicals.

Drug metabolism by cytochrome p450 enzymes: what distinguishes the pathways leading to substrate hydroxylation over desaturation?

Cytochrome P450 enzymes are highly versatile biological catalysts in our body that react with a broad range of substrates. Key functions in the liver include the metabolism of drugs and xenobiotics. One particular metabolic pathway that is poorly understood relates to the P450 activation of aliphatic groups leading to either hydroxylation or desaturation pathways. A DFT and QM/MM study has been carried out on the factors that determine the regioselectivity of aliphatic hydroxylation over desaturation of compounds by P450 isozymes. The calculations establish multistate reactivity patterns, whereby the product distributions differ on each of the spin-state surfaces; hence spin-selective product formation was found. The electronic and thermochemical factors that determine the bifurcation pathways were analysed and a model that predicts the regioselectivity of aliphatic hydroxylation over desaturation pathways was established from valence bond and molecular orbital theories. Thus, the difference in energy of the OH versus the OC bond formed and the π-conjugation energy determines the degree of desaturation products. In addition, environmental effects of the substrate binding pocket that affect the regioselectivities were identified. These studies imply that bioengineering P450 isozymes for desaturation reactions will have to include modifications in the substrate binding pocket to restrict the hydroxylation rebound reaction.

Clinical oxidative stress during leprosy multidrug therapy: impact of dapsone oxidation

This study aims to assess the oxidative stress in leprosy patients under multidrug therapy (MDT; dapsone, clofazimine and rifampicin), evaluating the nitric oxide (NO) concentration, catalase (CAT) and superoxide dismutase (SOD) activities, glutathione (GSH) levels, total antioxidant capacity, lipid peroxidation, and methemoglobin formation. For this, we analyzed 23 leprosy patients and 20 healthy individuals from the Amazon region, Brazil, aged between 20 and 45 years. Blood sampling enabled the evaluation of leprosy patients prior to starting multidrug therapy (called MDT 0) and until the third month of multidrug therapy (MDT 3). With regard to dapsone (DDS) plasma levels, we showed that there was no statistical difference in drug plasma levels between multibacillary (0.518±0.029 µg/mL) and paucibacillary (0.662±0.123 µg/mL) patients. The methemoglobin levels and numbers of Heinz bodies were significantly enhanced after the third MDT-supervised dose, but this treatment did not significantly change the lipid peroxidation and NO levels in these leprosy patients. In addition, CAT activity was significantly reduced in MDT-treated leprosy patients, while GSH content was increased in these patients. However, SOD and Trolox equivalent antioxidant capacity levels were similar in patients with and without treatment. These data suggest that MDT can reduce the activity of some antioxidant enzyme and influence ROS accumulation, which may induce hematological changes, such as methemoglobinemia in patients with leprosy. We also explored some redox mechanisms associated with DDS and its main oxidative metabolite DDS-NHOH and we explored the possible binding of DDS to the active site of CYP2C19 with the aid of molecular modeling software.

The role of drug metabolizing enzymes in clearance

INTRODUCTION

Metabolism is one of the most important clearance pathways representing the major clearance route of 75% drugs. The four most common drug metabolizing enzymes (DME) that contribute significantly to elimination pathways of new chemical entities are cytochrome P450s, UDP-glucuronosyltransferases, aldehyde oxidase and sulfotransferases. Accurate prediction of human in vivo clearance by these enzymes, using both in vitro and in vivo tools, is critical for the success of drug candidates in human translation.

AREAS COVERED

Important recent advances of key DME are reviewed and highlighted in the following areas: major isoforms, tissue distribution, generic polymorphism, substrate specificity, species differences, mechanism of catalysis, in vitro-in vivo extrapolation and the importance of using optimal assay conditions and relevant animal models.

EXPERT OPINION

Understanding the clearance mechanism of a compound is the first step toward successful prediction of human clearance. It is critical to apply appropriate in vitro and in vivo methodologies and physiologically based models in human translation. While high-confidence prediction for P450-mediated clearance has been achieved, the accuracy of human clearance prediction is significantly lower for other enzyme classes. More accurate predictive methods and models are being developed to address these challenges.

Use of density functional theory in drug metabolism studies

INTRODUCTION

The cytochrome P450 enzymes (CYPs) metabolize many drug compounds. They catalyze a wide variety of reactions, and potentially, a large number of different metabolites can be generated. Density functional theory (DFT) has, over the past decade, been shown to be a powerful tool to rationalize and predict the possible metabolites generated by the CYPs as well as other drug-metabolizing enzymes.

AREAS COVERED

We review applications of DFT on reactions performed by the CYPs and other drug-metabolizing enzymes able to perform oxidation reactions, with an emphasis on predicting which metabolites are produced. We also cover calculations of binding energies for complexes in which the ligands interact directly with the heme iron atom.

EXPERT OPINION

DFT is a useful tool for prediction of the site of metabolism. The use of small models of the enzymes work surprisingly well for most CYP isoforms. This is probably due to the fact that the binding of the substrates is not the major determinant. When binding of the substrate plays a significant role, the well-known issue of determining the free energy of binding is the challenge. How approaches taking the protein environment into account, like docking, MD and QM/MM, can be used are discussed.

Identification of key licorice constituents which interact with cytochrome P450: evaluation by LC/MS/MS cocktail assay and metabolic profiling

Licorice has been shown to affect the activities of several cytochrome P450 enzymes. This study aims to identify the key constituents in licorice which may affect these activities. Bioactivity assay was combined with metabolic profiling to identify these compounds in several complex licorice extracts. Firstly, the inhibition potencies of 40 pure licorice compounds were tested using an liquid chromatography/tandem mass spectrometry cocktail method. Significant inhibitors of human P450 isozymes 1A2, 2C9, 2C19, 2D6, and 3A4 were then selected for examination of their structural features by molecular docking to determine their molecular interaction with several P450 isozymes. Based on the present in vitro inhibition findings, along with our previous in vivo metabolic studies and the prevalence of individual compounds in licorice extract, we identified several licorice constituents, viz., liquiritigenin, isoliquiritigenin, together with seven isoprenylated flavonoids and arylcoumarins, which could be key components responsible for the herb-drug interaction between cytochrome P450 and licorice. In addition, hydrophilic flavonoid glycosides and saponins may be converted into these P450 inhibitors in vivo. These studies represent a comprehensive examination of the potential effects of licorice components on the metabolic activities of P450 enzymes.

Correlating structure and function of drug-metabolizing enzymes: progress and ongoing challenges

This report summarizes a symposium sponsored by the American Society for Pharmacology and Experimental Therapeutics at Experimental Biology held April 20-24 in Boston, MA. Presentations discussed the status of cytochrome P450 (P450) knowledge, emphasizing advances and challenges in relating structure with function and in applying this information to drug design. First, at least one structure of most major human drug-metabolizing P450 enzymes is known. However, the flexibility of these active sites can limit the predictive value of one structure for other ligands. A second limitation is our coarse-grain understanding of P450 interactions with membranes, other P450 enzymes, NADPH-cytochrome P450 reductase, and cytochrome b5. Recent work has examined differential P450 interactions with reductase in mixed P450 systems and P450:P450 complexes in reconstituted systems and cells, suggesting another level of functional control. In addition, protein nuclear magnetic resonance is a new approach to probe these protein/protein interactions, identifying interacting b5 and P450 surfaces, showing that b5 and reductase binding are mutually exclusive, and demonstrating ligand modulation of CYP17A1/b5 interactions. One desired outcome is the application of such information to control drug metabolism and/or design selective P450 inhibitors. A final presentation highlighted development of a CYP3A4 inhibitor that slows clearance of human immunodeficiency virus drugs otherwise rapidly metabolized by CYP3A4. Although understanding P450 structure/function relationships is an ongoing challenge, translational advances will benefit from continued integration of existing and new biophysical approaches.

A structural snapshot of CYP2B4 in complex with paroxetine provides insights into ligand binding and clusters of conformational states

An X-ray crystal structure of CYP2B4 in complex with the drug paroxetine [(3S,4R)-3-[(2H-1,3-benzodioxol-5-yloxy)methyl]-4-(4-fluorophenyl)piperidine] was solved at 2.14 Å resolution. The structure revealed a conformation intermediate to that of the recently solved complex with amlodipine and that of the more compact complex with 4-(4-chlorophenyl)imidazole in terms of the placement of the F-G cassette. Moreover, comparison of the new structure with 15 previously solved structures of CYP2B4 revealed some new insights into the determinants of active-site size and shape. The 2B4-paroxetine structure is nearly superimposable on a previously solved closed structure in a ligand-free state. Despite the overall conformational similarity among multiple closed structures, the active-site cavity volume of the paroxetine complex is enlarged. Further analysis of the accessible space and binding pocket near the heme reveals a new subchamber that resulted from the movement of secondary structural elements and rearrangements of active-site side chains. Overall, the results from the comparison of all 16 structures of CYP2B4 demonstrate a cluster of protein conformations that were observed in the presence or absence of various ligands.

Structural diversity of eukaryotic membrane cytochrome p450s

X-ray crystal structures are available for 29 eukaryotic microsomal, chloroplast, or mitochondrial cytochrome P450s, including two non-monooxygenase P450s. These structures provide a basis for understanding structure-function relations that underlie their distinct catalytic activities. Moreover, structural plasticity has been characterized for individual P450s that aids in understanding substrate binding in P450s that mediate drug clearance.

The contribution of atom accessibility to site of metabolism models for cytochromes P450

Three different types of atom accessibility descriptors are investigated in relation to site of metabolism predictions. To enable the integration of local accessibility we have constructed 2DSASA, a method for the calculation of the atomic solvent accessible surface area that is independent of 3D coordinates. The method was implemented in the SMARTCyp site of metabolism prediction models and improved the results by up to 4 percentage points for nine cytochrome P450 isoforms. The final models are made available at http://www.farma.ku.dk/smartcyp.