Docking and QSAR comparative studies of polycyclic aromatic hydrocarbons and other procarcinogen interactions with cytochromes P450 1A1 and 1B1

Dibenzyl trisulfide binds to and competitively inhibits the cytochrome P450 1A1 active site without impacting the expression of the aryl hydrocarbon receptor

The toxicological manifestation of many pollutants relies upon their binding to the aryl hydrocarbon receptor (AHR), and it follows a cascade of reactions culminating in an elevated expression of cytochrome P450 (CYP) 1 enzymes. CYP1A1 and CYP1B1 are associated with enhanced carcinogenesis when chronically exposed to certain polyaromatic hydrocarbons, and their inhibition may lead to chemoprevention. We evaluated dibenzyl trisulfide (DTS), expressed in the ethnomedical plant, Petiveria alliacea, for such potential chemoprevention. Using recombinant human CYP1A1 and CYP1B1 bactosomes on a fluorogenic assay, we first demonstrated that DTS moderately inhibited both enzymes with half maximal inhibitory concentration (IC50) values of 1.3 ± 0.3 and 1.7 ± 0.3 μM, respectively. Against CYP1A1, DTS was a reversible, competitive inhibitor with an apparent inhibitory constant (Ki) of 4.55 ± 0.37 μM. In silico molecular modeling showed that DTS binds with an affinity of -39.8 kJ·mol-1, situated inside the binding pocket, approximately 4.3 Å away from the heme group, exhibiting interactions with phenylalanine residue 123 (Phe-123), Phe-224, and Phe-258. Lastly, zebrafish (Danio rerio) embryos were exposed to 0.08-0.8 μM DTS from 24 to 96 h post fertilization (hpf) with the in vivo ethoxyresorufin-O-deethylase (EROD) assay, and, at 96 hpf, DTS significantly suppressed EROD CYP1A activity in a dose-dependent manner, with up to 60% suppression in the highest 0.8 μM exposure group. DTS had no impact on gene transcription levels for cyp1a and aryl hydrocarbon receptor 2 (ahr2). In co-exposure experiments, DTS suppressed CYP1A activity induced by both B[a]P and PCB-126, although these reductions were not significant. Taken together, these results demonstrate that DTS is a direct, reversible, competitive inhibitor of the carcinogen-activating CYP1A enzyme, binding in the active site pocket close to the heme site, and shows potential in chemoprevention.

Urinary polycyclic aromatic hydrocarbons and sex hormones in children and adolescents: Evidence from NHANES

BACKGROUND

Evidences showed that polycyclic aromatic hydrocarbons (PAHs) do harm to human body. However, the association between PAHs and sex hormones in children and adolescents remains unclear.

OBJECTIVES

The study aims to investigate the associations between PAHs and sex hormones in the general children and adolescent population.

METHODS

967 participants aged 6-19 with complete data of PAHs exposure biomarkers, covariates and sex hormones [total testosterone (TT), estradiol (E2) and sex hormone binding globulin (SHBG)] were recruited from National Health and Nutrition Examination Survey (NHANES), 2013-2016. Free androgen index (FAI) was calculated with TT/SHBG. Multivariate linear regression models were performed in six subgroups (male children, male adolescents, male late adolescents, female children, female adolescents and female late adolescents) to estimate the associations between sex hormone alterations and PAHs exposure.

RESULTS

In male puberty adolescents, weighted multivariate linear regression indicated that negative trends for 2-Hydroxynaphthalene, 1-Hydroxyphenanthrene, 2&3-Hydroxyphenanthrene and E2 (2-Hydroxynaphthalene: β: -0.104, 95%CI: -0.180, -0.029, P < 0.01; 1-Hydroxyphenanthrene: β: -0.112, 95%CI: -0.206, -0.018, P = 0.019; 2&3-Hydroxyphenanthrene: β: -0.125, 95%CI: -0.232, -0.018, P = 0.022), while exposure to 2-Hydroxynaphthalene was related to TT reduction (β: -0.099, 95%CI: -0.177, -0.020, P = 0.014). Same pattern between 2&3-Hydroxyphenanthrene and E2 alteration (2&3-Hydroxyphenanthrene: β: -0.139, 95%CI: -0.236, -0.041, P < 0.01) was also observed in male late adolescents. In male children, we determined that 1-Hydroxyphenanthrene was negatively associated with SHBG (β: -0.121, 95%CI: -0.205, -0.037, P < 0.01), while the same patterns were observed in male puberty children. We did not observe any significant result in female subgroups. All these results above were determined to have q value < 0.05.

CONCLUSION

PAHs exposure was associated with the alterations of sex hormones in male adolescents and children. Considering the cross-sectional study design, further large-scale epidemiological study is necessary.

Human CYP1A1 inhibition by flavonoids

Cytochrome P4501A1 (CYP1A1) is involved in the metabolism of several genotoxic/carcinogenic environmental xenobiotics including polycyclic aromatic hydrocarbons (PAHs) like benzo[a]pyrene. Several authors had proposed CYP1A inhibition as a plausible strategy for cancer chemoprevention. Using ethoxyresorufin O-deethylase activity (EROD), we tested the inhibitory properties of nine flavonoids: quercetin, miricetin, luteolin, fisetin, morin, kaempferol, 5-hydroxyflavone (5-HF), 3-hydroxyflavone (3-HF), and flavone (F) against human recombinant CYP1A1. The last three compounds exerted the highest inhibitory effect with IC₅₀ values of 0.07, 0.10 and 0.08 μM respectively; the more hydroxyl-groups were present, the lower the potency of inhibition was. Biochemical characterization leads to the conclusion that flavone and its hydroxy derivatives are mixed-type inhibitors. In silico studies have shown that, Phe224 and other aromatic residues in the human CYP1A1 active site play an important role in flavonoid-CYP interaction, through a π/π stacking between the aminoacid and the flavonoid C-ring. Outside the active site, the three flavonoids bind preferentially between A and K helices of the enzyme. Results from the Ames test using human S9 fraction revealed that none of the three compounds was mutagenic. We can consider 5-HF, 3-HF, and F as potential chemopreventive agents against genotoxic damage caused by metabolites resulting from CYP1A1 activity.

Copper-Catalyzed Synthesis, Bio-Evaluation, and in Silico Studies of 2-Aryl-N-alkylbenzimidazoles as Neuroprotective Agents

2-aryl-N-alkylbenzimidazole derivatives synthesized by CuI/PPh3 promoted direct coupling of N-alkylbenzimidazoles with aryl bromides. In vitro neurotoxicities of 20 compounds were evaluated, and the neuroprotective abilities of low-neurotoxic compounds (3b, 3g, 3h, 3i, 3j, 3k, 3o, 3q, 3s and 3t) were investigated against toxicity induced by 1-methyl-4-phenylpyridinium ion (MPP+) in SH-SY5Y neuronal cells. In silico studies revealed that compound 3g could have molecule docking with the following proteins: the bone morphogenetic protein receptor type 1B (BMPR1B), human cytochrome P450 1B1(CYP1B1), Metabotropic glutamate receptor 7 (GRM7), histone deacetylase 6 (HDAC6), 5-hydroxytryptamine receptor 5A (HTR5A), human topoisomerase II beta (TOP2B). A molecular docking simulation of model compound 3g and model protein CYP1B1 has been shown.

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.

Review of Ligand Specificity Factors for CYP1A Subfamily Enzymes from Molecular Modeling Studies Reported to-Date

The cytochrome P450 (CYP) family 1A enzymes, CYP1A1 and CYP1A2, are two of the most important enzymes implicated in the metabolism of endogenous and exogenous compounds through oxidation. These enzymes are also known to metabolize environmental procarcinogens into carcinogenic species, leading to the advent of several types of cancer. The development of selective inhibitors for these P450 enzymes, mitigating procarcinogenic oxidative effects, has been the focus of many studies in recent years. CYP1A1 is mainly found in extrahepatic tissues while CYP1A2 is the major CYP enzyme in human liver. Many molecules have been found to be metabolized by both of these enzymes, with varying rates and/or positions of oxidation. A complete understanding of the factors that govern the specificity and potency for the two CYP 1A enzymes is critical to the development of effective inhibitors. Computational molecular modeling tools have been used by several research groups to decipher the specificity and potency factors of the CYP1A1 and CYP1A2 substrates. In this review, we perform a thorough analysis of the computational studies that are ligand-based and protein-ligand complex-based to catalog the various factors that govern the specificity/potency toward these two enzymes.

Structural-functional adaptations of porcine CYP1A1 to metabolize polychlorinated dibenzo-p-dioxins

Polychlorinated dibenzo-p-dioxins (PCDDs) are widespread by-products of human industrial activity. They accumulate in tissues of animals and humans, exerting numerous adverse effects on different systems. In living organisms, dioxins are metabolized by enzymes of the cytochrome P450 family, including CYP1A1. Particular dioxin congeners differ in their toxicity level and ability to undergo biodegradation. Since the molecular mechanisms underlying dioxin susceptibility or resistance to biodegradation are unknown, in the present study the molecular interactions between five selected dioxins and porcine CYP1A1 protein were investigated. It was found that the ability of a dioxin to undergo CYP1A1-mediated degradation is associated mainly with the number and position of chlorine atoms in the dioxin molecule. Among all examined congeners, 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) demonstrated the highest affinity to CYP1A1 and, at the same time, the greatest distance to the active site of the enzyme. Interestingly, in contrast to other dioxins, the binding of the TCDD molecule to the porcine CYP1A1 active site resulted in a rapid and continuous closure of substrate channels. All the information may help to explain the extended half-life of TCDD in living organisms as well as its high toxicity.

Drug Metabolism in Preclinical Drug Development: A Survey of the Discovery Process, Toxicology, and Computational Tools

BACKGROUND

While establishing efficacy in translational models and humans through clinically-relevant endpoints for disease is of great interest, assessing the potential toxicity of a putative therapeutic drug is critical. Toxicological assessments in the pre-clinical discovery phase help to avoid future failure in the clinical phases of drug development. Many in vitro assays exist to aid in modular toxicological assessment, such as hepatotoxicity and genotoxicity. While these methods have provided tremendous insight into human toxicity by investigational new drugs, they are expensive, require substantial resources, and do not account for pharmacogenomics as well as critical ADME properties. Computational tools can fill this niche in toxicology if in silico models are accurate in relating drug molecular properties to toxicological endpoints as well as reliable in predicting important drug-target interactions that mediate known adverse events or adverse outcome pathways (AOPs).

METHODS

We undertook an unstructured search of multiple bibliographic databases for peer-reviewed literature regarding computational methods in predictive toxicology for in silico drug discovery. As this review paper is meant to serve as a survey of available methods for the interested reader, no focused criteria were applied. Literature chosen was based on the writers' expertise and intent in communicating important aspects of in silico toxicology to the interested reader.

CONCLUSION

This review provides a purview of computational methods of pre-clinical toxicologic assessments for novel small molecule drugs that may be of use for novice and experienced investigators as well as academic and commercial drug discovery entities.

Associations between CYP1A1 rs1048943 A > G and rs4646903 T > C genetic variations and colorectal cancer risk: Proof from 26 case-control studies

Cytochrome P450 1A1 (CYP1A1) enzyme is one of the most important metabolizing enzymes responsible for the metabolism of numerous xenobiotics. Numerous individual case-control studies have investigated the associations between the CYP1A1 rs1048943 A > G and rs4646903 T > C genetic variations and colorectal cancer (CRC) risk, but the conclusions were controversial. To obtain a scientific conclusion, we performed a meta-analysis based on a total of 26 publications, including 20 studies with 8665 cases and 9953 controls on rs1048943 A > G and 19 studies with 6416 cases and 7551 controls on rs4646903 T > C, respectively. The pooled analysis indicated that rs1048943 A > G was associated with an increased risk of CRC (G vs. A: OR = 1.28, 95% CI = 1.08-1.52; GG vs. AA: OR = 1.54, 95% CI = 1.25-1.91; GA vs. AA: OR = 1.26, 95% CI = 1.00-1.60; GG/GA vs. AA: OR = 1.31, 95% CI = 1.05-1.64; GG vs. GA/AA. OR = 1.56, 95% CI = 1.26-1.91). Stratification analysis showed the association between rs1048943 A > G and CRC risk was more obvious in studies with the population-based (PB) design or high quality score. The association between rs4646903 T > C and CRC risk did not reach statistical significance in the pooled analysis as well as stratification analysis. This meta-analysis demonstrated CYP1A1 rs1048943 A > G may increase the susceptibility to CRC instead of rs4646903 T > C. This conclusion suggested CYP1A1 may contribute to the pathogenesis of CRC.

Rational engineering of the fungal P450 monooxygenase CYP5136A3 to improve its oxidizing activity toward polycyclic aromatic hydrocarbons

A promising polycyclic aromatic hydrocarbon-oxidizing P450 CYP5136A3 from Phanerochaete chrysosporium was rationally engineered to enhance its catalytic activity. The residues W129 and L324 found to be critical in substrate recognition were transformed by single (L324F) and double (W129L/L324G, W129L/L324F, W129A/L324G, W129F/L324G and W129F/L324F) mutations, and the engineered enzyme forms were expressed in Pichia pastoris. L324F and W129F/L324F mutations enhanced the oxidation activity toward pyrene and phenanthrene. L324F also altered the regio-selectivity favoring C position 4 over 9 for hydroxylation of phenanthrene. This is the first instance of engineering a eukaryotic P450 for enhanced oxidation of these fused-ring hydrocarbons.

Insights on cytochrome p450 enzymes and inhibitors obtained through QSAR studies

The cytochrome P450 (CYP) superfamily of heme enzymes play an important role in the metabolism of a large number of endogenous and exogenous compounds, including most of the drugs currently on the market. Inhibitors of CYP enzymes have important roles in the treatment of several disease conditions such as numerous cancers and fungal infections in addition to their critical role in drug-drug interactions. Structure activity relationships (SAR), and three-dimensional quantitative structure activity relationships (3D-QSAR) represent important tools in understanding the interactions of the inhibitors with the active sites of the CYP enzymes. A comprehensive account of the QSAR studies on the major human CYPs 1A1, 1A2, 1B1, 2A6, 2B6, 2C9, 2C19, 2D6, 2E1, 3A4 and a few other CYPs are detailed in this review which will provide us with an insight into the individual/common characteristics of the active sites of these enzymes and the enzyme-inhibitor interactions.

Cytochrome P450 1A1-mediated anticancer drug discovery: in silico findings

INTRODUCTION

Target-specific drugs may offer fewer side/adverse effects in comparison with other anticancer agents and thus save normal healthy cells to a greater extent. The selective overexpression of cytochrome P450 1A1 (CYP1A1) in tumor cells induces the metabolism of benzothiazole and aminoflavone compounds to their reactive species, which are responsible for DNA adduct formation and cell death. This review encompasses the novelty of CYP1A1 as an anticancer drug target and explores the possible in silico strategies that would be applicable in the discovery and development of future antitumor compounds.

AREAS COVERED

This review highlights the various ligand-based and target-based in silico methodologies that were efficiently used in exploration of CYP1A1 as a novel antitumor target. These methodologies include electronic structure analysis, CoMFA studies, homology modeling, molecular docking, molecular dynamics analysis, pharmacophore mapping and quantitative structure activity relationship (QSAR) studies. It also focuses on the various approaches used in the development of the lysyl amide prodrug of 5F-203 (NSC710305) and dimethanesulfonate salt of 5-aminoflavone (NSC710464) as clinical candidates from their less potent analogues.

EXPERT OPINION

Selective overexpression of CYP1A1 in cancer cells offers tumor-specific drug design to ameliorate the current adverse effects associated with existing antitumor agents. Medicinal chemistry and in vitro driven approaches, in combination with knowledge-based drug design and by using the currently available tools of in silico methodologies, would certainly make it possible to design and develop novel anticancer compounds targeting CYP1A1.