Computational Study of Drugs Targeting Nuclear Receptors

In silico profiling of endocrine-disrupting potential of bisphenol analogues and their halogenated transformation products

Due to its endocrine-disrupting properties, bisphenol A (BPA) is being phased out from plastics, thermal paper and epoxy resins, and its replacements are being introduced into the market. Bisphenols are released into the environment, where they can undergo halogenation. Unlike BPA, the endocrine-disrupting potential of BPA analogues and their halogenated transformation products has not been extensively studied. The aim of this study was to evaluate the endocrine-disrupting potential of 18 BPA analogues and their halogenated derivatives by calculating affinities for 14 human nuclear receptors utilizing the Endocrine Disruptome and VirtualToxLab™ in silico tools. Our simulations identified AR, ERs, and GR as the most favorable targets of bisphenols and their derivatives. Several BPA analogues displayed a higher predicted potential for endocrine disruption than BPA. Our models highlighted BPZ and BPPH as the most hazardous in terms of predicted endocrine activities. Halogenation, in general, was predicted to increase the binding affinity of bisphenols for AR, ERβ, MR, GR, PPARγ, and TRβ. Notably, mono- or 2,2'-di-halogenated bisphenols exhibited the highest potential for endocrine disruption. In vitro corroboration of the obtained results should be the next milestone in evaluating the safety of BPA substitutes and their halogenated transformation products.

Effects of central nervous system drugs on androgen, estrogen α, glucocorticoid, and thyroid receptors

Some drugs that act on the central nervous system (CNS) are known to affect the endocrine system, although the mechanisms of endocrine toxicity are not well characterized to date. Such CNS drugs include antipsychotics, anticonvulsants, and antidepressants. In the present study, in-vitro firefly luciferase reporter-gene assays using the AR-EcoScreen assay using Chinese hamster ovary (CHO) cell line, hERα-HeLa9903, MDA-kb2, and GH3.TRE-Luc cell lines were used to determine the effects of nine CNS drugs on the androgen receptor, estrogen receptor α, glucocorticoid receptor, and thyroid hormone receptor, respectively. In the AR-EcoScreen assay using CHO cells, anti-androgenic activities were shown for carbamazepine (IC₅₀, 167 μM), clonazepam (IC₅₀, 26.7 μM), eslicarbazepine acetate (IC₅₀, 375 μM), fluoxetine (at 25 μM), lorazepam (IC₅₀, 16.4 μM), and sertraline (IC₅₀, 8.7 μM). In the hERα-HeLa-9903 cells, estrogen receptor α agonistic activities were shown for fluoxetine, paroxetine, and sertraline (at 10 μM and 25 μM), and in the GH3.TRE-Luc cells, the same three CNS drugs showed antithyroid activities (IC₅₀s, 11.6, 11.9, 2.7 μM, respectively). In the hERα-HeLa-9903 cells, estrogen receptor α antagonistic activities were shown for carbamazepine (IC₅₀, 114.3 μM), clonazepam (IC₅₀, 52.9 μM), and eslicarbazepine acetate (IC₅₀, 376.6 μM). When the CNS drugs were tested in the MDA-kb2 cells, none of them showed any activities toward glucocorticoid receptors. Little to no effects were seen toward any of these nuclear receptors for paliperidone and risperidone. The increased signal in the estrogen receptor α agonism assay seen for fluoxetine and paroxetine was confirmed to be mediated through estrogen receptor α. Additionally, we examined the interference of these CNS drugs with the firefly luciferase enzyme. These data elucidate the potential for adverse endocrine effects for some of these CNS drugs, which should therefore contribute to informed choice when prescribing them. However, long-term exposure to therapeutic concentrations of CNS drugs that have activities on the endocrine system should be explored further also in vivo.

In Silico Identification of Potential Thyroid Hormone System Disruptors among Chemicals in Human Serum and Chemicals with a High Exposure Index

Data on toxic effects are at large missing the prevailing understanding of the risks of industrial chemicals. Thyroid hormone (TH) system disruption includes interferences of the life cycle of the thyroid hormones and may occur in various organs. In the current study, high-throughput screening data available for 14 putative molecular initiating events of adverse outcome pathways, related to disruption of the TH system, were used to develop 19 in silico models for identification of potential thyroid hormone system-disrupting chemicals. The conformal prediction framework with the underlying Random Forest was used as a wrapper for the models allowing for setting the desired confidence level and controlling the error rate of predictions. The trained models were then applied to two different databases: (i) an in-house database comprising xenobiotics identified in human blood and ii) currently used chemicals registered in the Swedish Product Register, which have been predicted to have a high exposure index to consumers. The application of these models showed that among currently used chemicals, fewer were overall predicted as active compared to chemicals identified in human blood. Chemicals of specific concern for TH disruption were identified from both databases based on their predicted activity.

Effects of tyrosine kinase inhibitors on androgen, estrogen α, glucocorticoid and thyroid receptors

Modern anticancer therapies favor a targeted approach. Tyrosine kinase inhibitors (TKIs) are drugs that target molecular pathways involved in various types of malignancies. Although TKIs are safe and well tolerated, they remain not completely selective; e.g., endocrine-mediated adverse events have been observed with their use. In the present study, the effects of seven TKIs were determined on the activities of androgen receptor, estrogen receptor α (ERα), glucocorticoid receptor and thyroid receptor in vitro using stably transfected cell lines expressing firefly luciferase reporter gene: AR-EcoScreen, hERα-HeLa9903, MDA-kb2, and GH3.TRE-Luc cells, respectively. Antiandrogenic activity was seen for erlotinib, estrogenic activity for imatinib, antiestrogenic activity for dasatinib, erlotinib, nilotinib, regorafenib and sorafenib, glucocorticoid activity for erlotinib and ibrutinib, antiglucocorticoid activity for regorafenib and sorafenib, and antithyroid activity for ibrutinib. Additionally, synergism was seen for 1-5 μM dasatinib and 500 nM hydrocortisone combination for glucocorticoid activity in MDA-kb2 cells. The estrogenic activity of imatinib was confirmed as mediated through ERα, and interference of the TKIs with the reporter gene assays was ruled out in a cell-lysate-based firefly luciferase enzyme inhibition assay. Imatinib in combination with 4-hydroxytamoxifen showed concentration-dependent effects on the metabolic activity of ERα-expressing AN3CA, MCF-7, and SKOV3 cells, and on cell proliferation and adhesion of MCF-7 cells. These findings contribute to the understanding of the endocrine effects of TKIs, in terms of toxicity and effectiveness, and define the need to further evaluate the endocrine disrupting activities of TKIs to safeguard human and environmental health.

The effect of endocrine disrupting chemicals on the vitronectin-receptor (integrin αvβ3)-mediated cell adhesion of human umbilical vein endothelial cells

Endocrine disrupting chemicals (EDCs) are associated with cancer development and progression due to their promotion of increased cell invasiveness and metastasis formation. However, the effects of EDCs on cell adhesion mediated through integrins have not been well studied to date. Their actions are implicated by binding sites for hormones on the vitronectin receptor (VTNR; or integrin α_vβ₃), which is involved in tumor angiogenesis and metastasis. VTNR-expressing human umbilical vein endothelial cells (HUVECs) were used to determine the effects of EDCs and endogenous hormones on cell adhesion to vitronectin-coated surfaces, and on VTNR activation. Cell adhesion was significantly increased for bisphenol A, triclocarban, and triclosan (10, 100 nM; p < 0.05), with similar trends for bisphenols AF and S (10, 100 nM; p > 0.05). No changes in cell adhesion were seen for 5α-dihydrotestosterone, 17β-estradiol, triiodothyronine, imatinib and paroxetine. These data indicate that EDC-mediated increases in HUVEC adhesion to vitronectin are not mediated through androgenic, estrogenic, or thyroid activities, nor through activation of VTNR. Although these effects of EDCs on HUVEC adhesion require further investigation of the underlying mechanism(s) of action to define their biological relevance, the low-dose effects and nonmonotonic responses revealed here define the need for further investigation of these EDCs.

In silico prediction of nuclear receptor binding to polychlorinated dibenzofurans and its implication on endocrine disruption in humans and wildlife

Polychlorinated dibenzofurans (PCDFs) are known to cause endocrine disruption in humans and wildlife but the mechanisms underlying this disruption have not been adequately investigated. In this paper, the susceptibility of the endocrine system to disruption by PCDF congeners via nuclear receptor binding was studied using molecular docking simulation. Findings revealed that some PCDF congeners exhibit high probabilities of binding to androgen receptor in its agonistic and antagonistic conformations. In depth molecular docking analysis of the receptor-ligand complexes formed by PCDFs with androgen receptor in its agonistic and antagonistic conformations showed that, these complexes were stabilized by electrostatic, van der Waals, pi-effect and hydrophobic interactions. It was also observed that PCDF molecules mimic the modes of interaction observed in androgen-testosterone and androgen-bicalutamide complexes, utilizing between 65 and 83% of the amino acid residues used by the co-crystallized ligands for binding. This computational study suggests that some PCDF congeners may act as agonists and antagonists of androgen receptor in humans and wildlife via inapproprate binding to the receptor.

Incomplete thermal ablation-induced up-regulation of transcription factor nuclear receptor subfamily 2, group F, member 6 (NR2F6) contributes to the rapid progression of residual liver tumor in hepatoblastoma

Hepatoblastoma is a kind of extreme malignancy frequently diagnosed in children. Although surgical resection is considered as the first-line treatment for hepatoblastoma, a relatively large population of patients have lost the preferred opportunity for surgery. Administration of locoregional ablation enables local tumor control but with the deficiency of insufficient ablation, residual tumor, and rapid progression. In this study, we integrated 219 hepatoblastoma and 121 non-cancer liver tissues to evaluate the expression of NR2F6, from which a higher NR2F6 level was found in hepatoblastoma compared with non-cancer livers with a standard mean difference (SMD) of 1.04 (95% CI: 0.79, 1.29). The overexpression of NR2F6 also appeared to be an efficient indicator in distinguishing hepatoblastoma tissues from non-cancer liver tissues from the indication of a summarized AUC of 0.90, with a pooled sensitivity of 0.76 and a pooled specificity of 0.89. Interestingly, nude mouse xenografts provided direct evidence that overexpressed NR2F6 was also detected in residual tumor compared to untreated hepatoblastoma. Chromatin immunoprecipitation-binding data in HepG2 cells and transcriptome analysis of HepG2 xenografts were combined to identify target genes regulated by NR2F6. We finally selected 150 novel target genes of NR2F6 in residual tumor of incomplete ablation, and these genes appeared to be associated with the biological regulation of lipid metabolism-related pathway. Accordingly, targeting NR2F6 holds a therapeutic promise in treating residual recurrent hepatoblastoma after incomplete ablation.

A computational insight into endocrine disruption by polychlorinated biphenyls via non-covalent interactions with human nuclear receptors

Production of polychlorinated biphenyls (PCBs) was banned a long time ago because of their harmful health effects but humans continue to be exposed to residual PCBs in the environment. In this study, the susceptibility of human nuclear receptors to binding by PCBs was investigated using molecular docking simulation. Findings revealed that PCBs belonging to ortho-substituted, mono-ortho-substituted and non-ortho-substituted congeners could bind to agonistic conformations of androgen (AR), estrogen (ER α and ER β), glucocorticoid (GR) and thyroid hormone (TR α and TR β) receptors as well as antagonistic conformation of androgen receptor (AR an) but only ortho-substituted and mono-ortho-substituted PCBs could bind to estrogen receptors in their antagonistic conformations (ER α an and ER β an). Further molecular docking analyses showed that PCBs mimic the modes of interaction observed for the co-crystallized ligands in the crystal structures of the affected receptors, utilizing 81%, 83%, 78%, 60%, 75%, 60%, 86%, 100% and 75% of the amino acid residues utilized by the co-crystallized ligands for binding in AR, AR an, ER α, ER α an, ER β, ER β an, GR, TR α and TR β respectively. This computational study suggests that PCBs may cause endocrine disruption via formation of non-covalent interactions with androgen, estrogen, glucocorticoid and thyroid hormone receptors.

Identification of novel potential cyclooxygenase-2 inhibitors using ligand- and structure-based virtual screening approaches

Cyclooxygenase 2 (COX-2) is a well-established target for the design of anti-inflammatory intermediates. Celecoxib was selected as a template molecule to perform ligand-based virtual screening, i.e. to search for structures with similarity in shape and electrostatic potential, with a gradual increase in accuracy through the combined fitting of several steps using eight commercial databases. The molecules ZINC408709 and ZINC2090319 reproduced values within the limits established in an initial study of absorption and distribution in the body. No alert was fired for possible toxic groups when these molecules were subjected to toxicity prediction. Molecular docking results with these compounds showed a higher binding affinity in comparison to rofecoxib for the COX-2 target. Additionally, ZINC408709 and ZINC2090319 were predicted to be potentially biologically active. In in silico prediction of endocrine disruption potential, it was established that the molecule ZINC2090319 binds strongly to the target related to cardiovascular risk in a desirable way as a non-steroidal antagonist and ZINC408709 binds strongly to the target that is associated with the treatment of inflammatory pathologies and similar to celecoxib. Metabolites generated from these compounds are less likely to have side effects. Simulations were used to evaluate the interaction of compounds with COX-1 and COX-2 during 200 ns. Despite the differences, ZINC408709 molecule showed better stability for COX-2 during molecular dynamics simulation. In the calculations of free energy MM/PBSA, the molecule ZINC408709 ΔGbind value has a higher affinity to celecoxib and rofecoxib COX-2. This demonstrates that the selected substances can be considered as promising COX-2 inhibitors. Communicated by Ramaswamy H. Sarma.

Triclocarban, Triclosan, Bromochlorophene, Chlorophene, and Climbazole Effects on Nuclear Receptors: An in Silico and in Vitro Study

BACKGROUND

Endocrine-disrupting chemicals can interfere with hormonal homeostasis and have adverse effects for both humans and the environment. Their identification is increasingly difficult due to lack of adequate toxicological tests. This difficulty is particularly problematic for cosmetic ingredients, because in vivo testing is now banned completely in the European Union.

OBJECTIVES

The aim was to identify candidate preservatives as endocrine disruptors by in silico methods and to confirm endocrine receptors' activities through nuclear receptors in vitro.

METHODS

We screened preservatives listed in Annex V in the European Union Regulation on cosmetic products to predict their binding to nuclear receptors using the Endocrine Disruptome and VirtualToxLab™ version 5.8 in silico tools. Five candidate preservatives were further evaluated for androgen receptor (AR), estrogen receptor (ER α ), glucocorticoid receptor (GR), and thyroid receptor (TR) agonist and antagonist activities in cell-based luciferase reporter assays in vitro in AR-EcoScreen, hER α -HeLa- 9903 , MDA-kb2, and GH3.TRE-Luc cell lines. Additionally, assays to test for false positives were used (nonspecific luciferase gene induction and luciferase inhibition).

RESULTS

Triclocarban had agonist activity on AR and ER α at 1 μ M and antagonist activity on GR at 5 μ M and TR at 1 μ M . Triclosan showed antagonist effects on AR, ER α , GR at 10 μ M and TR at 5 μ M , and bromochlorophene at 1 μ M (AR and TR) and at 10 μ M (ER α and GR). AR antagonist activity of chlorophene was observed [inhibitory concentration at 50% (IC50) IC 50 = 2.4 μ M ], as for its substantial ER α agonist at > 5 μ M and TR antagonist activity at 10 μ M . Climbazole showed AR antagonist (IC 50 = 13.6 μ M ), ER α agonist at > 10 μ M , and TR antagonist activity at 10 μ M .

DISCUSSION

These data support the concerns of regulatory authorities about the endocrine-disrupting potential of preservatives. These data also define the need to further determine their effects on the endocrine system and the need to reassess the risks they pose to human health and the environment. https://doi.org/10.1289/EHP6596.