Preferred binding orientations of phenacetin in CYP1A1 and CYP1A2 are associated with isoform-selective metabolism

Uncoupling of Cytochrome P450 2B6 and stimulation of reactive oxygen species production in pharmacogenomic alleles affected by interethnic variability

Cytochrome P450 mediated substrate metabolism is generally characterized by the formation of reactive intermediates. In vitro and in vivo reaction uncoupling, results in the accumulation and dissociation of reactive intermediates, leading to increased ROS formation. The susceptibility towards uncoupling and altered metabolic activity is partly modulated by pharmacogenomic alleles resulting in amino acid substitutions. A large variability in the prevalence of these alleles has been demonstrated in CYP2B6, with some being predominantly unique to African populations. The aim of this study is to characterize the uncoupling potential of recombinant CYP2B6*1, CYP2B6*6 and CYP2B6*34 metabolism of specific substrates. Therefore, functional effects of these alterations on enzyme activity were determined by quantification of bupropion, efavirenz and ketamine biotransformation using HPLC-MS/MS. Determination of H2O2 levels was performed by the AmplexRed/horseradish peroxidase assay. Our studies of the amino acid substitutions Q172H, K262R and R487S revealed an exclusive use of the peroxide shunt for the metabolism of bupropion and ketamine by CYP2B6*K262R. Ketamine was also identified as a trigger for the peroxide shunt in CYP2B6*1 and all variants. Concurrently, ketamine acted as an uncoupler for all enzymes. We further showed that the expressed CYP2B6*34 allele results in the highest H2O2 formation. We therefore conclude that the reaction uncoupling and peroxide shunt are directly linked and can be substrate specifically induced with K262R carriers being most likely to use the peroxide shunt and R487S carrier being most prone to reaction uncoupling. This elucidates the functional diversity of pharmacogenomics in drug metabolism and safety.

In vitro differences in toddalolactone metabolism in various species and its effect on cytochrome P450 expression

CONTEXT

Toddalolactone, the main component of Toddalia asiatica (L.) Lam. (Rutaceae), has anticancer, antihypertension, anti-inflammatory, and antifungal activities.

OBJECTIVE

This study investigated the metabolic characteristics of toddalolactone.

MATERIALS AND METHODS

Toddalolactone metabolic stabilities were investigated by incubating toddalolactone (20 μM) with liver microsomes from humans, rabbits, mice, rats, dogs, minipigs, and monkeys for 0, 30, 60, and 90 min. The CYP isoforms involved in toddalolactone metabolism were characterized based on chemical inhibition studies and screening assays. The effects of toddalolactone (0, 10, and 50 µM) on CYP1A1 and CYP3A5 protein expression were investigated by immunoblotting. After injecting toddalolactone (10 mg/kg), in vivo pharmacokinetic profiles using six Sprague-Dawley rats were investigated by taking 9-time points, including 0, 0.25, 0.5, 0.75, 1, 2, 4, 6 and 8 h.

RESULTS

Monkeys showed the greatest metabolic capacity in CYP-mediated and UGT-mediated reaction systems with short half-lives (T_1/2) of 245 and 66 min, respectively, while T_1/2 of humans in two reaction systems were 673 and 83 min, respectively. CYP1A1 and CYP3A5 were the major CYP isoforms involved in toddalolactone biotransformation. Induction of CYP1A1 protein expression by 50 μM toddalolactone was approximately 50% greater than that of the control (0 μM). Peak plasma concentration (C_max) for toddalolactone was 0.42 μg/mL, and T_max occurred at 0.25 h post-dosing. The elimination t_1/2 was 1.05 h, and the AUC_0-t was 0.46 μg/mL/h.

CONCLUSIONS

These findings demonstrated the significant species differences of toddalolactone metabolic profiles, which will promote appropriate species selection in further toddalolactone studies.

Candidate Proficiency Test Chemicals to Address Industrial Chemical Applicability Domains for in vitro Human Cytochrome P450 Enzyme Induction

Cytochrome P450 (CYP) enzymes play a key role in the metabolism of both xenobiotics and endogenous chemicals, and the activity of some CYP isoforms are susceptible to induction and/or inhibition by certain chemicals. As CYP induction/inhibition can bring about significant alterations in the level of in vivo exposure to CYP substrates and metabolites, CYP induction/inhibition data is needed for regulatory chemical toxicity hazard assessment. On the basis of available human in vivo pharmaceutical data, a draft Organisation for Economic Co-operation and Development Test Guideline (TG) for an in vitro CYP HepaRG test method that is capable of detecting the induction of four human CYPs (CYP1A1/1A2, 2B6, and 3A4), has been developed and validated for a set of pharmaceutical proficiency chemicals. However to support TG adoption, further validation data was requested to demonstrate the ability of the test method to also accurately detect CYP induction mediated by industrial and pesticidal chemicals, together with an indication on regulatory uses of the test method. As part of "GOLIATH", a European Union Horizon-2020 funded research project on metabolic disrupting chemical testing approaches, work is underway to generate supplemental validated data for an additional set of chemicals with sufficient diversity to allow for the approval of the guideline. Here we report on the process of proficiency chemical selection based on a targeted literature review, the selection criteria and considerations required for acceptance of proficiency chemical selection for OECD TG development (i.e. structural diversity, range of activity, relevant chemical sectors, global restrictions etc). The following 13 proposed proficiency chemicals were reviewed and selected as a suitable set for use in the additional validation experiments: tebuconazole, benfuracarb, atrazine, cypermethrin, chlorpyrifos, perfluorooctanoic acid, bisphenol A, N,N-diethyl-m-toluamide, benzo-[a]-pyrene, fludioxonil, malathion, triclosan, and caffeine. Illustrations of applications of the test method in relation to endocrine disruption and non-genotoxic carcinogenicity are provided.

Deep Neural Network-Assisted Drug Recommendation Systems for Identifying Potential Drug-Target Interactions

In silico methods to identify novel drug-target interactions (DTIs) have gained significant importance over conventional techniques owing to their labor-intensive and low-throughput nature. Here, we present a machine learning-based multiclass classification workflow that segregates interactions between active, inactive, and intermediate drug-target pairs. Drug molecules, protein sequences, and molecular descriptors were transformed into machine-interpretable embeddings to extract critical features from standard datasets. Tools such as CHEMBL web resource, iFeature, and an in-house developed deep neural network-assisted drug recommendation (dNNDR)-featx were employed for data retrieval and processing. The models were trained with large-scale DTI datasets, which reported an improvement in performance over baseline methods. External validation results showed that models based on att-biLSTM and gCNN could help predict novel DTIs. When tested with a completely different dataset, the proposed models significantly outperformed competing methods. The validity of novel interactions predicted by dNNDR was backed by experimental and computational evidence in the literature. The proposed methodology could elucidate critical features that govern the relationship between a drug and its target.

The relationships between cytochromes P450 and H2O2: Production, reaction, and inhibition

In this review we address the relationship between cytochromes P450 (P450) and H2O2. This association can affect biology in three distinct ways. First, P450s produce H2O2 as a byproduct either during catalysis or when no substrate is present. This reaction, known as uncoupling, releases reactive oxygen species that may have implications in disease. Second, H2O2 is used as an oxygen-donating co-substrate in peroxygenase and peroxidase reactions catalyzed by P450s. This activity has proven to be important mainly in reactions involving prokaryotic P450s, and investigators have harnessed this reaction with the aim of adaptation for industrial use. Third, H2O2-dependent inhibition of human P450s has been studied in our laboratory, demonstrating heme destruction and also the inactivating oxidation of the heme-thiolate ligand to a sulfenic acid (-SOH). This reversible oxidative modification of P450s may have implications in the prevention of uncoupling and may give new insights into the oxidative regulation of these enzymes. Research has elucidated many of the chemical mechanisms involved in the relationship between P450 and H2O2, but the application to biology is difficult to evaluate. Further studies are needed reveal both the harmful and protective natures of reactive oxygen species in an organismal context.

Structure-Based Drug Design for Cytochrome P450 Family 1 Inhibitors

Cytochromes P450 are a class of metalloproteins which are responsible for electron transfer in a wide spectrum of reactions including metabolic biotransformation of endogenous and exogenous substrates. The superfamily of cytochromes P450 consists of families and subfamilies which are characterized by a specific structure and substrate specificity. Cytochromes P450 family 1 (CYP1s) play a distinctive role in the metabolism of drugs and chemical procarcinogens. In recent decades, these hemoproteins have been intensively studied with the use of computational methods which have been recently developed remarkably to be used in the process of drug design by the virtual screening of compounds in order to find agents with desired properties. Moreover, the molecular modeling of proteins and ligand docking to their active sites provide an insight into the mechanism of enzyme action and enable us to predict the sites of drug metabolism. The review presents the current status of knowledge about the use of the computational approach in studies of ligand-enzyme interactions for CYP1s. Research on the metabolism of substrates and inhibitors of CYP1s and on the selectivity of their action is particularly valuable from the viewpoint of cancer chemoprevention, chemotherapy, and drug-drug interactions.

Increased Phenacetin Oxidation upon the L382V Substitution in Cytochrome P450 1A2 is Associated with Altered Substrate Binding Orientation

Leucine382 of cytochrome P450 1A2 (CYP1A2) plays an important role in binding and O-dealkylation of phenacetin, with the L382V mutation increasing substrate oxidation (Huang and Szklarz, 2010, Drug Metab. Dispos. 38:1039–1045). This was attributed to altered substrate binding orientation, but no direct experimental evidence had been available. Therefore, in the current studies, we employed nuclear magnetic resonance (NMR) longitudinal (T₁) relaxation measurements to investigate phenacetin binding orientations within the active site of CYP1A2 wild type (WT) and mutants. Paramagnetic relaxation time (T_1P) for each proton of phenacetin was calculated from the T₁ value obtained from the enzymes in ferric and ferrous-CO state in the presence of phenacetin, and used to model the orientation of phenacetin in the active site. All aromatic protons of phenacetin were nearly equidistant from the heme iron (6.34–8.03 Å). In contrast, the distance between the proton of the –OCH₂– group, which is abstracted during phenacetin oxidation, and the heme iron, was much shorter in the L382V (5.93 Å) and L382V/N312L (5.96 Å) mutants compared to the N312L mutant (7.84 Å) and the wild type enzyme (6.55 Å), consistent with modeling results. These studies provide direct evidence for the molecular mechanism underlying increased oxidation of phenacetin upon the L382V mutation.

CYP1A1-Mediated Intramolecular Rearrangement of Aminoazepane in GDC-0339

GDC-0339 is a novel small molecule pan-Pim kinase inhibitor that was discovered as a potential treatment of multiple myeloma. During the in vitro and in vivo metabolite profiling of GDC-0339, a metabolite was detected that had the same elemental composition as the parent but was distinct with respect to its chromatographic separation and mass spectrometric fragmentation pattern. High resolution tandem mass spectrometry data indicated the metabolite was modified at the aminoazepane moiety. The structure was solved by nuclear magnetic resonance analysis of the isolated metabolite and further confirmed by comparing it to a synthetic standard. These results indicated that the metabolite was formed by an intramolecular amine replacement reaction with the primary amine forming a new attachment to pyrazole without any change in stereochemistry. In vitro experiments showed cytochrome P450s catalyzed the reaction and demonstrated high isoform selectivity by CYP1A1. Results from kinetic experiments showed that the CYP1A1-mediated rearrangement of GDC-0339 was an efficient reaction with apparent turnover number (kcat) and Michaelis constant (Km) of 8.4 minutes-1 and 0.6 μM, respectively. The binding of GDC-0339 to the cytochrome P450 active site was examined by characterizing the direct inhibition of CYP1A1-mediated phenacetin O-deethylation, and GDC-0339 was a potent competitive inhibitor with Ki of 0.9 μM. This high affinity binding was unexpected given a narrow active site for CYP1A1 and GDC-0339 does not conform structurally to known CYP1A1 substrates, which are mostly polyaromatic planar molecules. Further, we explored some of the structural requirements for the rearrangement reaction and identified several analogs to GDC-0339 that undergo this biotransformation.

Characterization of differences in substrate specificity among CYP1A1, CYP1A2 and CYP1B1: an integrated approach employing molecular docking and molecular dynamics simulations

Recent trends in new drug discovery of anticancer drugs have made oncologists more aware of the fact that the new drug discovery must target the developing mechanism of tumorigenesis to improve the therapeutic efficacy of antineoplastic drugs. The drugs designed are expected to have high affinity towards the novel targets selectively. Current research highlights overexpression of CYP450s, particularly cytochrome P450 1A1 (CYP1A1), in tumour cells, representing a novel target for anticancer therapy. However, the CYP1 family is identified as posing significant problems in selectivity of anticancer molecules towards CYP1A1. Three members have been identified in the human CYP1 family: CYP1A1, CYP1A2 and CYP1B1. Although sequences of the three isoform have high sequence identity, they have distinct substrate specificities. The understanding of macromolecular features that govern substrate specificity is required to understand the interplay between the protein function and dynamics, design novel antitumour compounds that could be specifically metabolized by only CYP1A1 to mediate their antitumour activity and elucidate the reasons for differences in substrate specificity profile among the three proteins. In the present study, we employed a combination of computational methodologies: molecular docking and molecular dynamics simulations. We utilized eight substrates for elucidating the difference in substrate specificity of the three isoforms. Lastly, we conclude that the substrate specificity of a particular substrate depends upon the type of the active site residues, the dynamic motions in the protein structure upon ligand binding and the physico-chemical characteristics of a particular ligand. Copyright © 2016 John Wiley & Sons, Ltd.

Comparative proteomics among cytochrome p450 family 1 for differential substrate specificity

Apart from playing key roles in drug metabolism and adverse drug-drug interactions, CYPs are potential drug targets to treat a variety of diseases. The intervention of over expression of P450 1A1 (CYP1A1) in tumor cells is identified as a novel strategy for anticancer therapy. We investigated three isoforms of CYP1 family (CYP1A1, CYP1A2, and CYP1B1) for their substrate specificity. The understanding of macromolecular features that govern substrate specificity is required to understand the interplay between the protein function and dynamics. This can help in design of new antitumor molecule specifically metabolized by CYP1A1 to mediate their antitumor activity. In the present study, we carried out the comparative protein structure analysis of the three isoforms. Sequence alignment, root mean square deviation (RMSD) analysis, B-factor analysis was performed to give a better understanding of the macromolecular features involved in substrate specificity and to understand the interplay between protein dynamics and functions which will have important implications on rational design of anticancer drugs. We identified the differences in amino acid residues among the three isoforms of CYP1 family, which may account for differential substrate specificity. Six putative substrate recognition sequences are characterized along with the regions they form in the protein structure. Further the RMSD and B-factor analysis provides the information about the identified residues having the maximum RMSD and B-factor deviations.

Paramagnetic nuclear magnetic resonance relaxation and molecular mechanics studies of the chloroperoxidase-indole complex: insights into the mechanism of chloroperoxidase-catalyzed regioselective oxidation of indole

To unravel the mechanism of chloroperoxidase (CPO)-catalyzed regioselective oxidation of indole, we studied the structure of the CPO-indole complex using nuclear magnetic resonance (NMR) relaxation measurements and computational techniques. The dissociation constant (KD) of the CPO-indole complex was calculated to be approximately 21 mM. The distances (r) between protons of indole and the heme iron calculated via NMR relaxation measurements and molecular docking revealed that the pyrrole ring of indole is oriented toward the heme with its 2-H pointing directly at the heme iron. Both KD and r values are independent of pH in the range of 3.0-6.5. The stability and structure of the CPO-indole complex are also independent of the concentration of chloride or iodide ion. Molecular docking suggests the formation of a hydrogen bond between the NH group of indole and the carboxyl O of Glu 183 in the binding of indole to CPO. Simulated annealing of the CPO-indole complex using r values from NMR experiments as distance restraints reveals that the van der Waals interactions were much stronger than the Coulomb interactions in the binding of indole to CPO, indicating that the association of indole with CPO is primarily governed by hydrophobic rather than electrostatic interactions. This work provides the first experimental and theoretical evidence of the long-sought mechanism that leads to the "unexpected" regioselectivity of the CPO-catalyzed oxidation of indole. The structure of the CPO-indole complex will serve as a lighthouse in guiding the design of CPO mutants with tailor-made activities for biotechnological applications.

Human cytochrome P450 1A1 structure and utility in understanding drug and xenobiotic metabolism

Cytochrome P450 (CYP) 1A1 is an extrahepatic monooxygenase involved in the metabolism of endogenous substrates and drugs, as well as the activation of certain toxins and environmental pollutants. CYP1A1 is particularly well known for its ability to biotransform polycyclic aromatic hydrocarbons, such as benzo[a]pyrene in tobacco smoke, into carcinogens. CYP1A1 possesses functional similarities and differences with human CYP1A2 and CYP1B1 enzymes, but the structural basis for this has been unclear. We determined a 2.6 Å structure of human CYP1A1 with the inhibitor α-naphthoflavone. α-Naphthoflavone binds within an enclosed active site, with the planar benzochromen-4-one core packed flat against the I helix that composes one wall of the active site, and the 2-phenyl substituent oriented toward the catalytic heme iron. Comparisons with previously determined structures of the related cytochrome P450 1A2 and 1B1 enzymes reveal distinct features among the active sites that may underlie the functional variability of these enzymes. Finally, docking studies probed the ability of CYP1A structures to assist in understanding their known in vitro interactions with several typical substrates and inhibitors.

Quantitative synthesis of the association between the cytochrome P450 1A1 Ile462Val polymorphism and prostate cancer risk

The association between the cytochrome P450 1A1 (CYP1A1) Ile462Val polymorphism and prostate cancer risk remains inconclusive owing to the conflicting findings from previous studies. To get a more precise estimate of the possible association, we performed the present meta-analysis. We searched the PUBMED, EMBASE, and Wanfang databases for the studies which met the inclusion criteria. The pooled odds ratio (OR) with corresponding 95 % confidence interval (95 % CI) was used to estimate the association between CYP1A1 Ile462Val polymorphism and prostate cancer risk. A total of 13 studies with 2,350 cases and 2,992 controls were included in the meta-analysis. The results indicated that there was an obvious association between CYP1A1 Ile462Val polymorphism and increased risk of prostate cancer (for Val versus Ile: OR=1.27, 95 % CI 1.13-1.43, P<0.001; for ValVal versus IleIle: OR=1.51, 95 % CI 1.14-2.01, P=0.004; for ValVal + ValIle versus IleIle: OR=1.31, 95 % CI 1.14-1.51, P<0.001; for ValVal versus IleIle + ValIle: OR=1.38, 95 % CI 1.05-1.81, P=0.020). Subgroup analyses by ethnicity suggested that CYP1A1 Ile462Val polymorphism was associated with prostate cancer risk in Asians but not in Caucasians. This meta-analysis suggests that there is an association between CYP1A1 Ile462Val polymorphism and increased risk of prostate cancer. More studies with large sample are needed to further assess the association in Caucasians.