New aryl hydrocarbon receptor homology model targeted to improve docking reliability

The AhR Signaling Mechanism: A Structural Point of View

The Aryl hydrocarbon Receptor (AhR) is a well-known sensor of xenobiotics; moreover, it is considered a promising drug target as it is involved in the regulation of many patho-physiological processes. For these reasons the study of its ligand-activated transcription mechanism has stimulated several studies for over twenty years. In this review we highlight the key role of molecular structural information in understanding the different steps of the signaling mechanism. The architecture of the AhR cytosolic complex, encompassing the hsp90 chaperone protein and the XAP2 and p23 co-chaperones, has become available in the last year thanks to Cryo-EM experiments. The structure of the AhR ligand-binding (PAS-B) domain has remained elusive for a long time; it has been predicted by homology modelling, based on known PAS systems, and its ligand-bound forms were modelled through ligand molecular docking. Although very recently some structural information on this domain has become available, considerable efforts are still needed to determine the binding geometries of the AhR key ligands by experimental high-resolution studies. On the other hand, the dimeric structure of AhR with the ARNT protein, bound to the specific DNA responsive element, was partially determined by X-ray crystallography and it was completed by homology modelling. On the whole the current structural knowledge of the main protein complexes that form over the AhR mechanism opens the way to confirm and further investigate the main steps of the proposed ligand-activated transcription mechanism of the AhR.

Gene-Environment Interactions: My Unique Journey

I am deeply honored to be invited to write this scientific autobiography. As a physician-scientist, pediatrician, molecular biologist, and geneticist, I have authored/coauthored more than 600 publications in the fields of clinical medicine, biochemistry, biophysics, pharmacology, drug metabolism, toxicology, molecular biology, cancer, standardized gene nomenclature, developmental toxicology and teratogenesis, mouse genetics, human genetics, and evolutionary genomics. Looking back, I think my career can be divided into four distinct research areas, which I summarize mostly chronologically in this article: (a) discovery and characterization of the AHR/CYP1 axis, (b) pharmacogenomics and genetic prediction of response to drugs and other environmental toxicants, (c) standardized drug-metabolizing gene nomenclature based on evolutionary divergence, and (d) discovery and characterization of the SLC39A8 gene encoding the ZIP8 metal cation influx transporter. Collectively, all four topics embrace gene-environment interactions, hence the title of my autobiography.

Microparticle Phosphatidylserine Mediates Coagulation: Involvement in Tumor Progression and Metastasis

Tumor progression and cancer metastasis has been linked to the release of microparticles (MPs), which are shed upon cell activation or apoptosis and display parental cell antigens, phospholipids such as phosphatidylserine (PS), and nucleic acids on their external surfaces. In this review, we highlight the biogenesis of MPs as well as the pathophysiological processes of PS externalization and its involvement in coagulation activation. We review the available evidence, suggesting that coagulation factors (mainly tissue factor, thrombin, and fibrin) assist in multiple steps of tumor dissemination, including epithelial-mesenchymal transition, extracellular matrix remodeling, immune escape, and tumor angiogenesis to support the formation of the pre-metastatic niche. Platelets are not just bystander cells in circulation but are functional players in primary tumor growth and metastasis. Tumor-induced platelet aggregation protects circulating tumor cells (CTCs) from the blood flow shear forces and immune cell attack while also promoting the binding of CTCs to endothelial cells and extravasation, which activates tumor invasion and sustains metastasis. Finally, in terms of therapy, lactadherin can inhibit coagulation by competing effectively with coagulation factors for PS binding sites and may similarly delay tumor progression. Furthermore, we also investigate the therapeutic potential of coagulation factor inhibitors within the context of cancer treatment. The development of multiple therapies targeting platelet activation and platelet-tumor cell interactions may not only reduce the lethal consequences of thrombosis but also impede tumor growth and spread.

Is it really safe to replace decabromodiphenyl ether (BDE209) with decabromodiphenyl ethane (DBDPE)?: A perspective from hepatotoxicity

In this paper, the hepatocytotoxicity and aryl hydrocarbon receptor (AHR) activity of decabromodiphenyl ethane (DBDPE), decabromodiphenyl ether (BDE209) and other 18 analogues were evaluated in vitro using human normal liver cell L02. These dioxin-like compounds showed differential hepatocytotoxicity (EC₅₀  = 0.38-17.87 mg/L) and AHR activity (EROD activity = 4.53-46.35 U/μg). In silico study indicated the distance of π-π bonds between the benzene ring of compounds and residue Phe234 of AHR played a key role in the binding of AHR, and the substituents on the benzene ring also influenced the activity. Combining molecular biology and bioomics, the comprehensive investigations on the hepatotoxic mechanisms have demonstrated the AHR signaling pathway was the key mediation mechanism for the hepatotoxicity of DBDPE/BDE209. The cytochrome P450s (CYP2 family) mediated formation of reactive oxygenated intermediates might be the dominant toxic mechanism, which could produce oxidative stress or cause genotoxicity. Although the experimental toxicity of DBDPE was smaller relative to BDE209, the health risk of DBDPE may be much greater than we expected, due to the high potential to form a variety of dioxin-like intermediates by microbial oxidation of ethyl group. Therefore, whether it is really safe to replace BDE209 with DBDPE is a debatable question, and more ecotoxicological and health data are needed to clarify this issue.

In-depth analysis of the interactions of various aryl hydrocarbon receptor ligands from a computational perspective

Aryl hydrocarbon receptor (AhR) is a ligand-dependent transcription factor that acts as a machinery that controls the expression of many genes, including cytochrome P450 CYP1A1, CYP1A2 and CYP1B1. It plays a principal role in numerous biological and toxicological functions, making it a promising target for developing therapeutic agents. Several novel small molecules targeting the AhR signaling pathway are currently under investigation as antitumor agents. Some have already advanced into clinical trials in patients with various tumors. Activation of AhR by diverse chemicals either endogenous or exogenous is initiated by the binding of these ligands to the PAS-B domain, which modulates AhR functions. There is, however, limited information about how various ligands interact with the PAS-B domain for activating or inhibiting the AhR. To better understand the mode of action of AhR agonists/antagonists. The current work proposes a combination of several computational tools to build dynamical models for the PAS-B domain bound to different ligands in mouse and human. Our findings reveal the essential roles of specific PAS-B residues (e.g., S365, V381& Q383), which mediate the AhR ligand-binding process. Our results also explain how these residues regulate the promiscuity of AhR in accommodating various chemicals in its binding PAS-B ligand-binding pocket.

The Highly Potent AhR Agonist Picoberin Modulates Hh-Dependent Osteoblast Differentiation

Identification and analysis of small molecule bioactivity in target-agnostic cellular assays and monitoring changes in phenotype followed by identification of the biological target are a powerful approach for the identification of novel bioactive chemical matter in particular when the monitored phenotype is disease-related and physiologically relevant. Profiling methods that enable the unbiased analysis of compound-perturbed states can suggest mechanisms of action or even targets for bioactive small molecules and may yield novel insights into biology. Here we report the enantioselective synthesis of natural-product-inspired 8-oxotetrahydroprotoberberines and the identification of Picoberin, a low picomolar inhibitor of Hedgehog (Hh)-induced osteoblast differentiation. Global transcriptome and proteome profiling revealed the aryl hydrocarbon receptor (AhR) as the molecular target of this compound and identified a cross talk between Hh and AhR signaling during osteoblast differentiation.

Structural insight into the ligand binding mechanism of aryl hydrocarbon receptor

The aryl hydrocarbon receptor (AHR), a member of the basic helix-loop-helix (bHLH) Per-Arnt-Sim (PAS) family of transcription factors, plays important roles in regulating xenobiotic metabolism, cellular differentiation, stem cell maintenance, as well as immunity. More recently, AHR has gained significant interest as a drug target for the development of novel cancer immunotherapy drugs. Detailed understanding of AHR-ligand binding has been hampered for decades by the lack of a three-dimensional structure of the AHR PAS-B domain. Here, we present multiple crystal structures of the Drosophila AHR PAS-B domain, including its apo, ligand-bound, and AHR nuclear translocator (ARNT) PAS-B-bound forms. Together with biochemical and cellular assays, our data reveal structural features of the AHR PAS-B domain, provide insights into the mechanism of AHR ligand binding, and provide the structural basis for the future development of AHR-targeted therapeutics.

Aryl Hydrocarbon Receptor-Dependent and -Independent Pathways Mediate Curcumin Anti-Aging Effects

The aryl hydrocarbon receptor (AhR) is a ligand-activated transcription factor whose activity can be modulated by polyphenols, such as curcumin. AhR and curcumin have evolutionarily conserved effects on aging. Here, we investigated whether and how the AhR mediates the anti-aging effects of curcumin across species. Using a combination of in vivo, in vitro, and in silico analyses, we demonstrated that curcumin has AhR-dependent or -independent effects in a context-specific manner. We found that in Caenorhabditis elegans, AhR mediates curcumin-induced lifespan extension, most likely through a ligand-independent inhibitory mechanism related to its antioxidant activity. Curcumin also showed AhR-independent anti-aging activities, such as protection against aggregation-prone proteins and oxidative stress in C. elegans and promotion of the migratory capacity of human primary endothelial cells. These AhR-independent effects are largely mediated by the Nrf2/SKN-1 pathway.

Phytochemicals From Vicia faba Beans as Ligands of the Aryl Hydrocarbon Receptor to Regulate Autoimmune Diseases

Legumes are associated with gut health benefits, and increasing evidence indicates that their consumption reduces the risk of chronic diseases that include autoimmunity. Beans are rich sources of compounds with health-promoting effects, and recent metabolomic approaches have enabled the comprehensive characterization of the chemical composition of Vicia faba L. This article reviewed whether the phytocompounds in broad beans might modulate the aryl hydrocarbon receptor (AhR), which plays an essential role in autoantigen tolerance as a potential dietary strategy for autoimmune disease management. Therefore, thirty molecules present in Vicia faba of the chemical classes of flavonoids, chalcones, stilbenes, jasmonates, alkaloids, and amino acids, and either a human- or microbiome-derived product of biotransformation, retrieved from the literature or predicted in silico were evaluated by docking for affinity against the ligand-binding domain of AhR. Most analyzed compounds showed high affinity even after their metabolism which indicate that some AhR modulators remain active despite several steps in their biotransformation. Hence, our results suggest that in similitude with the gut metabolism of the tryptophan, phytocompounds mainly polyphenols also lead to metabolites that induce the AhR pathway. Furthermore, wyerone acid, wyerone epoxide, jasmonic acid, stizolamine, vicine, and convicine and their metabolite derivatives are reported for the first time as potential AhR ligands. Overall, chronic consumption of phytochemicals in Vicia faba L. and their gut biotransformation may protect against autoimmune disease pathogenesis by AhR modulation.

Review of in silico studies dedicated to the nuclear receptor family: Therapeutic prospects and toxicological concerns

Being in the center of both therapeutic and toxicological concerns, NRs are widely studied for drug discovery application but also to unravel the potential toxicity of environmental compounds such as pesticides, cosmetics or additives. High throughput screening campaigns (HTS) are largely used to detect compounds able to interact with this protein family for both therapeutic and toxicological purposes. These methods lead to a large amount of data requiring the use of computational approaches for a robust and correct analysis and interpretation. The output data can be used to build predictive models to forecast the behavior of new chemicals based on their in vitro activities. This atrticle is a review of the studies published in the last decade and dedicated to NR ligands in silico prediction for both therapeutic and toxicological purposes. Over 100 articles concerning 14 NR subfamilies were carefully read and analyzed in order to retrieve the most commonly used computational methods to develop predictive models, to retrieve the databases deployed in the model building process and to pinpoint some of the limitations they faced.

The aryl hydrocarbon receptor as a model PAS sensor

Proteins containing PER-ARNT-SIM (PAS) domains are commonly associated with environmental adaptation in a variety of organisms. The PAS domain is found in proteins throughout Archaea, Bacteria, and Eukarya and often binds small-molecules, supports protein-protein interactions, and transduces input signals to mediate an adaptive physiological response. Signaling events mediated by PAS sensors can occur through induced phosphorelays or genomic events that are often dependent upon PAS domain interactions. In this perspective, we briefly discuss the diversity of PAS domain containing proteins, with particular emphasis on the prototype member, the aryl hydrocarbon receptor (AHR). This ligand-activated transcription factor acts as a sensor of the chemical environment in humans and many chordates. We conclude with the idea that since mammalian PAS proteins often act through PAS-PAS dimers, undocumented interactions of this type may link biological processes that we currently think of as independent. To support this idea, we present a framework to guide future experiments aimed at fully elucidating the spectrum of PAS-PAS interactions with an eye towards understanding how they might influence environmental sensing in human and wildlife populations.

Identification of viable TCDD access pathways to human AhR PAS-B ligand binding domain

Unintentionally released in the environment as by-products of industrial activities, dioxins, exemplified by 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), represent a primary concern for human health. Exposure to these chemicals is known to produce a broad spectrum of adverse effects, including cancer. The main mechanism of action of TCDD in humans involves binding to the Aryl hydrocarbon Receptor (AhR). Although qualitatively established, TCDD capture by the AhR remains poorly characterized at the molecular level. Starting from a recently developed structural model of the human AhR PAS-B domain, in this work we attempt the identification of viable TCDD access pathways to the human AhR ligand binding domain by means of molecular dynamics. Based on the result of metadynamics simulations, we identify two main regions that may potentially serve as access paths for TCDD. For each path, we characterize the residues closely interacting with TCDD, thereby suggesting a possible mechanism for TCDD capture. Our results are reviewed and discussed in the light of the available information about Human AhR structure and functions.

Biochemical properties of human full-length aryl hydrocarbon receptor (AhR)

The aryl hydrocarbon receptor (AhR) is a very unstable protein. AhR binds to the molecular chaperone complex (HSP90-p23-XAP2) to maintain a stable structure in the cytoplasm. After binding to ligands, such as dioxin, AhR translocates from the cytoplasm to the nucleus with a molecular chaperone complex. The protein forms a heterodimer with Arnt after nuclear transfer, functions as a transcription factor by binding to a xenobiotic responsive element (XRE), and induces the cytochrome P450 1A1 (CYP1A1). Because of the unstable protein, expression of the full-length AhR in the E. coli expression system is very difficult. Many studies investigated AhR using AhR domains in vitro. We expressed and purified the human full-length AhR in E. coli expression system. Furthermore, specific antibodies were prepared. Purified full-length AhR could bind to ligand. In the presence of ligand, α-helix and random coil of AhR increased and β-sheet decreased on CD spectrum. Full-length AhR could bind to HSP90, XAP2 and p23 in the presence or absence of ligand. We now show the biochemical properties of full-length AhR.

Improvement in the screening performance of potential aryl hydrocarbon receptor ligands by using supervised machine learning

The aryl hydrocarbon receptor (AhR), which is a ligand-dependent transcription factor, plays a crucial role in the regulation of xenobiotic metabolism. There are a large number of artificial or natural molecules in the environment that can activate AhR. In this study, we developed a virtual screening procedure to identify potential ligands of AhR. One structure-based method and two ligand-based methods were used for the virtual screening procedure. The results showed that the precision rate (0.96) and recall rate (0.64) of our procedure were significantly higher than those of a procedure used in a previous study, which suggests that supervised machine learning techniques can greatly improve the performance of virtual screening. Moreover, a pesticide dataset including 777 frequently used pesticides was screened. Seventy-seven pesticides were identified as potential AhR ligands by all three screening methods, among which 12 have never been previously reported as AhR agonists. Two non-agonist AhR ligands and 14 of the 77 pesticides were randomly selected for testing by in vitro and in vivo assays. All 14 pesticides showed different degrees of AhR agonistic activity, and none of the two non-agonist AhR ligand pesticides showed AhR agonistic activity, which suggests that our procedure had good robustness. Four of the pesticides were reported as AhR agonists for the first time, suggesting that these pesticides may need further toxicity assessment. In general, our procedure is a rapid, powerful and computationally inexpensive tool for predicting chemicals with AhR agonistic activity, which could be useful for environmental risk prediction and management.

An aryl hydrocarbon receptor from the caecilian Gymnopis multiplicata suggests low dioxin affinity in the ancestor of all three amphibian orders

The aryl hydrocarbon receptor (AHR) plays pleiotropic roles in the development and physiology of vertebrates in conjunction with xenobiotic and endogenous ligands. It is best known for mediating the toxic effects of dioxin-like pollutants such as 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). While most vertebrates possess at least one AHR that binds TCDD tightly, amphibian AHRs bind TCDD with very low affinity. Previous analyses of AHRs from Xenopus laevis (a frog; order Anura) and Ambystoma mexicanum (a salamander; order Caudata) identified three amino acid residues in the ligand-binding domain (LBD) that underlie low-affinity binding. In X. laevis AHR1β, these are A354, A370, and N325. Here we extend the analysis of amphibian AHRs to the caecilian Gymnopis multiplicata, representing the remaining extant amphibian order, Gymnophiona. G. multiplicata AHR groups with the monophyletic vertebrate AHR/AHR1 clade. The LBD includes all three signature residues of low TCDD affinity, and a structural homology model suggests that its architecture closely resembles those of other amphibians. In transactivation assays, the EC50 for reporter gene induction by TCDD was 17.17 nM, comparable to X. laevis AhR1β (26.23 nM) and Ambystoma AHR (34.09 nM) and dramatically higher than mouse AhR (0.13 nM), a trend generally reflected in direct measures of TCDD binding. These shared properties distinguish amphibian AHRs from the high-affinity proteins typical of both vertebrate groups that diverged earlier (teleost fish) and those that appeared more recently (other tetrapods). These findings suggest the hypothesis that AHRs with low TCDD affinity represent a characteristic that evolved in a common ancestor of all three extant amphibian groups.

Hypoxia and Oxygen-Sensing Signaling in Gene Regulation and Cancer Progression

Oxygen homeostasis regulation is the most fundamental cellular process for adjusting physiological oxygen variations, and its irregularity leads to various human diseases, including cancer. Hypoxia is closely associated with cancer development, and hypoxia/oxygen-sensing signaling plays critical roles in the modulation of cancer progression. The key molecules of the hypoxia/oxygen-sensing signaling include the transcriptional regulator hypoxia-inducible factor (HIF) which widely controls oxygen responsive genes, the central members of the 2-oxoglutarate (2-OG)-dependent dioxygenases, such as prolyl hydroxylase (PHD or EglN), and an E3 ubiquitin ligase component for HIF degeneration called von Hippel–Lindau (encoding protein pVHL). In this review, we summarize the current knowledge about the canonical hypoxia signaling, HIF transcription factors, and pVHL. In addition, the role of 2-OG-dependent enzymes, such as DNA/RNA-modifying enzymes, JmjC domain-containing enzymes, and prolyl hydroxylases, in gene regulation of cancer progression, is specifically reviewed. We also discuss the therapeutic advancement of targeting hypoxia and oxygen sensing pathways in cancer.

The Cellular and Molecular Determinants of Naphthoquinone-Dependent Activation of the Aryl Hydrocarbon Receptor

1,2-naphthoquinone (1,2-NQ) and 1,4-naphthoquinone (1,4-NQ) are clinically promising biologically active chemicals that have been shown to stimulate the aryl hydrocarbon receptor (AhR) signaling pathway, but whether they are direct or indirect ligands or activate the AhR in a ligand-independent manner is unknown. Given the structural diversity of AhR ligands, multiple mechanisms of AhR activation of gene expression, and species differences in AhR ligand binding and response, we examined the ability of 1,2-NQ and 1,4-NQ to bind to and activate the mouse and human AhRs using a series of in vitro AhR-specific bioassays and in silico modeling techniques. Both NQs induced AhR-dependent gene expression in mouse and human hepatoma cells, but were more potent and efficacious in human cells. 1,2-NQ and 1,4-NQ stimulated AhR transformation and DNA binding in vitro and was inhibited by AhR antagonists. Ligand binding analysis confirmed the ability of 1,2-NQ and 1,4-NQ to competitively bind to the AhR ligand binding cavity and the molecular determinants for interactions were predicted by molecular modeling methods. NQs were shown to bind distinctly differently from that of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and differences were also observed between species. Mutation of amino acid residues (F289, M334, and M342) involved in critical NQ:AhR binding interactions, decreased NQ- and AhR-dependent gene expression, consistent with a role for these residues in binding and activation of the AhR by NQs. These studies provide insights into the molecular mechanism of action of NQs and contribute to the development of emerging NQ-based therapeutics.

Marine Polyhydroxynaphthoquinone, Echinochrome A: Prevention of Atherosclerotic Inflammation and Probable Molecular Targets

The effect of low doses of echinochrome A (EchA), a natural polyhydroxy-1,4-naphthoquinone pigment from the sea urchin Scaphechinus mirabilis, has been studied in clinical trials, when it was used as an active substance of the drug Histochrome® and biologically active supplement Thymarin. Several parameters of lipid metabolism, antioxidant status, and the state of the immune system were analyzed in patients with cardiovascular diseases (CVD), including contaminating atherosclerosis. It has been shown that EchA effectively normalizes lipid metabolism, recovers antioxidant status and reduces atherosclerotic inflammation, regardless of the method of these preparations' administrations. Treatment of EchA has led to the stabilization of patients, improved function of the intracellular matrix and decreased epithelial dysfunction. The increased expression of surface human leukocyte antigen DR isotype (HLA-DR) receptors reflects the intensification of intercellular cooperation of immune cells, as well as an increase in the efficiency of processing and presentation of antigens, while the regulation of CD95 + expression levels suggests the stimulation of cell renewal processes. The immune system goes to a different level of functioning. Computer simulations suggest that EchA, with its aromatic structure of the naphthoquinone nucleus, may be a suitable ligand of the cytosolic aryl cell receptor, which affects the response of the immune system and causes the rapid expression of detoxification enzymes such as CYP and DT diaphorase, which play a protective role with CVD. Therefore, EchA possesses not only an antiradical effect and antioxidant activity, but is also a SOD3 mimetic, producing hydrogen peroxide and controlling the expression of cell enzymes through hypoxia-inducible factors (HIF), peroxisome proliferator-activated receptors (PPARs) and aryl hydrocarbon receptor (AhR).

Antimigraine Drug Avitriptan Is a Ligand and Agonist of Human Aryl Hydrocarbon Receptor That Induces CYP1A1 in Hepatic and Intestinal Cells

The efforts for therapeutic targeting of the aryl hydrocarbon receptor (AhR) have emerged in recent years. We investigated the effects of available antimigraine triptan drugs, having an indole core in their structure, on AhR signaling in human hepatic and intestinal cells. Activation of AhR in reporter gene assays was observed for Avitriptan and to a lesser extent for Donitriptan, while other triptans were very weak or no activators of AhR. Using competitive binding assay and by homology docking, we identified Avitriptan as a low-affinity ligand of AhR. Avitriptan triggered nuclear translocation of AhR and increased binding of AhR in CYP1A1 promotor DNA, as revealed by immune-fluorescence microscopy and chromatin immune-precipitation assay, respectively. Strong induction of CYP1A1 mRNA was achieved by Avitriptan in wild type but not in AhR-knockout, immortalized human hepatocytes, implying that induction of CYP1A1 is AhR-dependent. Increased levels of CYP1A1 mRNA by Avitriptan were observed in human colon carcinoma cells LS180 but not in primary cultures of human hepatocytes. Collectively, we show that Avitriptan is a weak ligand and activator of human AhR, which induces the expression of CYP1A1 in a cell-type specific manner. Our data warrant the potential off-label therapeutic application of Avitriptan as an AhR-agonist drug.

Crosstalk between hydroxytyrosol, a major olive oil phenol, and HIF-1 in MCF-7 breast cancer cells

Olive oil intake has been linked with a lower incidence of breast cancer. Hypoxic microenvironment in solid tumors, such as breast cancer, is known to play a crucial role in cancer progression and in the failure of anticancer treatments. HIF-1 is the foremost effector in hypoxic response, and given that hydroxytyrosol (HT) is one of the main bioactive compounds in olive oil, in this study we deepen into its modulatory role on HIF-1. Our results in MCF-7 breast cancer cells demonstrate that HT decreases HIF-1α protein, probably by downregulating oxidative stress and by inhibiting the PI3K/Akt/mTOR pathway. Strikingly, the expression of HIF-1 target genes does not show a parallel decrease. Particularly, adrenomedullin and vascular endothelial growth factor are up-regulated by high concentrations of HT even in HIF-1α silenced cells, pointing to HIF-1-independent mechanisms of regulation. In fact, we show, by in silico modelling and transcriptional analysis, that high doses of HT may act as an agonist of the aryl hydrocarbon receptor favoring the induction of these angiogenic genes. In conclusion, we suggest that the effect of HT in a hypoxic environment is largely affected by its concentration and involves both HIF-1 dependent and independent mechanisms.

Transitional States in Ligand-Dependent Transformation of the Aryl Hydrocarbon Receptor into Its DNA-Binding Form

The aryl hydrocarbon receptor (AhR) is a ligand-activated transcription factor that mediates the biological and toxicological effects of an AhR lacking the entire PASB structurally diverse chemicals, including halogenated aromatic hydrocarbons. Ligand-dependent transformation of the AhR into its DNA binding form involves a ligand-dependent conformational change, heat shock protein 90 (hsp90), dissociation from the AhR complex and AhR dimerization with the AhR nuclear translocator (ARNT) protein. The mechanism of AhR transformation was examined using mutational approaches and stabilization of the AhR:hsp90 complex with sodium molybdate. Insertion of a single mutation (F281A) in the hsp90-binding region of the AhR resulted in its constitutive (ligand-independent) transformation/DNA binding in vitro. Mutations of AhR residues within the Arg-Cys-rich region (R212A, R217A, R219A) and Asp371 (D371A) impaired AhR transformation without a significant effect on ligand binding. Stabilization of AhR:hsp90 binding with sodium molybdate decreased transformation/DNA binding of the wild type AhR but had no effect on constitutively active AhR mutants. Interestingly, transformation of the AhR in the presence of molybdate allowed detection of an intermediate transformation ternary complex containing hsp90, AhR, and ARNT. These results are consistent with a stepwise transformation mechanism in which binding of ARNT to the liganded AhR:hsp90 complex results in a progressive displacement of hsp90 and conversion of the AhR into its high affinity DNA binding form. The available molecular insights into the signaling mechanism of other Per-ARNT-Sim (PAS) domains and structural information on hsp90 association with other client proteins are consistent with the proposed transformation mechanism of the AhR.

Assessment of Risks of Dioxins for Aryl Hydrocarbon Receptor-Mediated Effects in Polar Bear (Ursus maritimus) by in Vitro and in Silico Approaches

Polar bear (Ursus maritimus) populations accumulate dioxins and related compounds (DRCs) at levels that are of health concern. The toxicities of DRCs are primarily mediated via aryl hydrocarbon receptor (AHR) signaling pathway. To evaluate the sensitivity and responses to DRCs in polar bears, we assessed the activation potencies of polar bear-specific AHR (pbAHR) by DRCs through in vitro and in silico approaches. In vitro assays showed that the pbAHR was as sensitive to DRCs as C3H/lpr mouse AHR, which is well-known to be highly sensitive to DRCs. Comparison of pbAHR transactivation potencies indicated that TCDF, 2,3,4,7,8-PeCDF, and BaP exhibited high induction equivalency factors (IEFs). Considering the accumulation levels of DRCs in polar bears, PCB126 was found to be the most active inducer of pbAHR. The in vitro transactivation potencies of ligands of pbAHR showed a significant relationship with in silico ligand docking energies in a pbAHR homology model. The protein ligand interaction fingerprint (PLIF) analysis showed different interaction patterns depending on the ligands. Several amino acids which are highly conserved among mammals may be involved in species-specific responses via backbone interactions with neighboring amino acid residues which are specific to pbAHR. We document high susceptibility of polar bears to DRCs, through a mechanistic approach, for the first time.

Regulation of Aryl Hydrocarbon Receptor Signaling Pathway and Dioxin Toxicity by Novel Agonists and Antagonists

Dioxins, mostly through activation of aryl hydrocarbon receptor (AhR), are potent toxic substances widely distributed in the environment, while moderated suppression of AhR also exhibits anti-tumor effects. Therefore, the proper modulation of AhR activity may counteract AhR-mediated toxicities and certain diseases. In this investigation, we identified several novel AhR moderate agonists and antagonists using chemical biology approaches. The mechanisms and mode of interactions with AhR by these hits were also revealed using both experimental and computational studies. The newly identified AhR moderate agonists and antagonists were predicted to bind to AhR and modulate AhR signaling. The structure-activity relationships of moderate agonists and antagonists and their unique binding features with AhR have created a solid framework for further optimization of the next generation of AhR modulators.

Modulation of lymphocyte-mediated tissue repair by rational design of heterocyclic aryl hydrocarbon receptor agonists

Aryl hydrocarbon receptor (AHR) is an essential regulator of gut immunity and a promising therapeutic target for inflammatory bowel disease (IBD). Current AHR agonists are inadequate for clinical translation due to low activity, inadequate pharmacokinetics, or toxicity. We synthesized a structurally diverse library and used integrated computational and experimental studies to discover mechanisms governing ligand-receptor interaction and to design potent drug leads PY109 and PY108, which display physiochemical drug-likeness properties, desirable pharmacokinetic profiles, and low toxicity. In a murine model of dextran sulfate sodium-induced colitis, orally administered compounds increase interleukin-22 (IL-22) production and accelerate mucosal healing by modulating mucosal adaptive and innate lymphoid cells. AHR and IL-22 pathway induction was confirmed using RNA sequencing and characterization of the lymphocyte protein-protein interaction network. Significant induction of IL-22 was also observed using human T cells from patients with IBD. Our findings support rationally designed AHR agonists for IBD therapy.

The aryl hydrocarbon receptor (AhR) as a breast cancer drug target

Breast cancer is the most common cancer in women, with more than 1.7 million diagnoses worldwide per annum. Metastatic breast cancer remains incurable, and the presence of triple-negative phenotypes makes targeted treatment impossible. The aryl hydrocarbon receptor (AhR), most commonly associated with the metabolism of xenobiotic ligands, has emerged as a promising biological target for the treatment of this deadly disease. Ligands for the AhR can be classed as exogenous or endogenous and may have agonistic or antagonistic activity. It has been well reported that agonistic ligands may have potent and selective growth inhibition activity in a number of oncogenic cell lines, and one (aminoflavone) has progressed to phase I clinical trials for breast cancer sufferers. In this study, we examine the current state of the literature in this area and elucidate the promising advances that are being made in hijacking the cytosolic-to-nuclear pathway of the AhR for the possible future treatment of breast cancer.

A Biomimetic, One-Step Transformation of Simple Indolic Compounds to Malassezia-Related Alkaloids with High AhR Potency and Efficacy

Malassezia furfur isolates from diseased skin preferentially biosynthesize compounds which are among the most active known aryl-hydrocarbon receptor (AhR) inducers, such as indirubin, tryptanthrin, indolo[3,2-b]carbazole, and 6-formylindolo[3,2-b]carbazole. In our effort to study their production from Malassezia spp., we investigated the role of indole-3-carbaldehyde (I3A), the most abundant metabolite of Malassezia when grown on tryptophan agar, as a possible starting material for the biosynthesis of the alkaloids. Treatment of I3A with H₂O₂ and use of catalysts like diphenyldiselenide resulted in the simultaneous one-step transformation of I3A to indirubin and tryptanthrin in good yields. The same reaction was first applied on simple indole and then on substituted indoles and indole-3-carbaldehydes, leading to a series of mono- and bisubstituted indirubins and tryptanthrins bearing halogens, alkyl, or carbomethoxy groups. Afterward, they were evaluated for their AhR agonist activity in recombinant human and mouse hepatoma cell lines containing a stably transfected AhR-response luciferase reporter gene. Among them, 3,9-dibromotryptanthrin was found to be equipotent to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) as an AhR agonist, and 3-bromotryptanthrin was 10-times more potent than TCDD in the human HG2L7.5c1 cell line. In contrast, 3,9-dibromotryptanthrin and 3-bromotryptanthrin were ∼4000 and >10,000 times less potent than TCDD in the mouse H1L7.5c3 cell line, respectively, demonstrating that they are species-specific AhR agonists. Involvement of the AhR in the action of 3-bromotryptanthrin was confirmed by the ability of the AhR antagonists CH223191 and SR1 to inhibit 3-bromotryptanthrin-dependent reporter gene induction in human HG2L7.5c1 cells. In conclusion, I3A can be the starting material used by Malassezia for the production of both indirubin and tryptanthrin through an oxidation mechanism, and modification of these compounds can produce some highly potent, efficacious and species-selective AhR agonists.

Modeling the binding of diverse ligands within the Ah receptor ligand binding domain

The Ah receptor (AhR) is a ligand-dependent transcription factor belonging to the basic helix-loop-helix Per-Arnt-Sim (bHLH-PAS) superfamily. Binding to and activation of the AhR by a variety of chemicals results in the induction of expression of diverse genes and production of a broad spectrum of biological and toxic effects. The AhR also plays important roles in several physiological responses, which has led it to become a novel target for the development of therapeutic drugs. Differences in the interactions of various ligands within the AhR ligand binding domain (LBD) may contribute to differential modulation of AhR functionality. We combined computational and experimental analyses to investigate the binding modes of a group of chemicals representative of major classes of AhR ligands. On the basis of a novel computational approach for molecular docking to the homology model of the AhR LBD that includes the receptor flexibility, we predicted specific residues within the AhR binding cavity that play a critical role in binding of three distinct groups of chemicals. The prediction was validated by site-directed mutagenesis and evaluation of the relative ligand binding affinities for the mutant AhRs. These results provide an avenue for understanding ligand modulation of the AhR functionality and for rational drug design.

The aryl hydrocarbon receptor 2 potentially mediates cytochrome P450 1A induction in the jungle crow (Corvus macrorhynchos)

To understand the role of aryl hydrocarbon receptor (AHR) isoforms in avian species, we investigated the functional characteristics of two AHR isoforms (designated as jcAHR1 and jcAHR2) of the jungle crow (Corvus macrorhynchos). Two amino acid residues corresponding to Ile³²⁴ and Ser³⁸⁰ (high sensitive type) in chicken AHR1 that are known to determine dioxin sensitivity were Ile³²⁵ and Ala³⁸¹ (moderate sensitive type) in jcAHR1 and Val³⁰⁶ and Ala³⁶² (low sensitive type) in jcAHR2. The quantitative comparison of the two jcAHR mRNA expression levels in a Tokyo jungle crow population showed that jcAHR2 accounted for 92.4% in the liver, while jcAHR1 accounted for only 7.6%. Both in vitro-expressed jcAHR1 and jcAHR2 proteins exhibited a specific binding to [³H]-labeled 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). Transactivation potencies for jcAHR1 and jcAHR2 in in vitro reporter gene assays were measured in jcAHR-expressed cells exposed to 16 dioxins and related compounds (DRCs). Both jcAHR1 and jcAHR2 were activated in a congener- and an isoform-specific manner. EC₅₀ value of TCDD for jcAHR2 (0.61 nM) was six-fold higher than that for jcAHR1 (0.098 nM), but jcAHR2 had higher transactivation efficacy than jcAHR1 in terms of the magnitude of response. The high transactivation efficacy of jcAHR2 in DRCs is in contrast to that of AHR2s in other avian species with low transactivation efficacy. Molecular docking simulations of TCDD with in silico jcAHR1 and jcAHR2 homology models showed that the two sensitivity-decisive amino acids indirectly controlled TCDD-binding modes through their surrounding amino acids. Deletion assays of jcAHR2 revealed that 736-805 amino acid residues in the C-terminal region were critical for its transactivation. We suggest that jcAHR2 plays a critical role in regulating the AHR signaling pathway, at least in its highly expressed organs.

Comparative In Vitro and In Silico Analysis of the Selectivity of Indirubin as a Human Ah Receptor Agonist

The aryl hydrocarbon receptor (AhR) is a ligand-dependent transcription factor that modulates gene expression following its binding and activation by structurally diverse chemicals. Species differences in AhR functionality have been observed, with the mouse AhR (mAhR) and human AhR (hAhR) exhibiting significant differences in ligand binding, coactivator recruitment, gene expression and response. While the AhR agonist indirubin (IR) is a more potent activator of hAhR-dependent gene expression than the prototypical ligand 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), it is a significantly less potent activator of the mAhR. DNA binding analysis confirmed the greater potency/efficacy of IR in stimulating transformation/DNA binding of the hAhR in vitro and domain-swapping experiments demonstrated that the enhanced response to IR was primarily due to the hAhR ligand binding domain (LBD). Site-directed mutagenesis and functional analysis studies revealed that mutation of H326 and A349 in the mAhR LBD to the corresponding residues in the hAhR LBD significantly increased the potency of IR. Since these mutations had no significant effect on ligand binding, these residues likely contribute to an enhanced efficiency of transformation/DNA binding by IR-bound hAhR. Molecular docking to mAhR LBD homology models further elucidated the different roles of the A375V mutation in TCDD and IR binding, as revealed by [³H]TCDD competitive binding results. These results demonstrate the differential binding of structurally diverse ligands within the LBD of a given AhR and confirm that amino acid differences within the LBD of AhRs contribute to significant species differences in ligand response.

In vitro toxicity and in silico docking analysis of two novel selective AH-receptor modulators

The mediator of dioxin toxicity, aryl hydrocarbon receptor (AHR), has also important physiological functions. Selective AHR modulators (SAHRMs) share some effects of dioxins, except for their marked toxicity. We recently characterised toxicologically two novel SAHRMs, prodrugs IMA-08401 and IMA-07101 in rats, demonstrating that they are far less deleterious than the most toxic AHR-agonist, TCDD. Here, we analysed the in vitro toxicity and in silico AHR binding of the respective active, deacetylated metabolites, IMA-06201 (N-ethyl-N-phenyl-5-chloro-1,2-dihydro-4-hydroxy-1-methyl-2-oxo-quinoline-3-carboxamide) and IMA-06504 (N-(4-trifluoromethylphenyl)-1,2-dihydro-4-hydroxy-5-methoxy-1-methyl-2-oxo-quinoline-3-carboxamide). In H4IIE rat hepatoma cells, IMA-06201 and IMA-06504 induced CYP1A1 with comparable potencies and efficacies to those of TCDD. They had little effect on cell viability as assessed by LDH leakage and MTT reduction assays, and were not mutagenic in the Ames test, but IMA-06504 elicited a maximally 2.7-fold increase in micronuclei. Molecular docking simulations showed that similar to TCDD, they occupy the central region of AHR ligand binding cavity. Hence, while showing low to negligible in vitro toxicity, these novel SAHRMs bind to the AHR qualitatively in a similar fashion to TCDD, and appear comparably powerful AHR agonists. Combined with our earlier results demonstrating that they seem considerably less toxic in vivo than TCDD, these compounds are thus highly interesting new SAHRMs.

Methylindoles and Methoxyindoles are Agonists and Antagonists of Human Aryl Hydrocarbon Receptor

Novel methylindoles were identified as endobiotic and xenobiotic ligands of the human aryl hydrocarbon receptor (AhR). We examined the effects of 22 methylated and methoxylated indoles on the transcriptional activity of AhRs. Employing reporter gene assays in AZ-AHR transgenic cells, we determined full agonist, partial agonist, or antagonist activities of tested compounds, having substantially variable EC₅₀, IC₅₀, and relative efficacies. The most effective agonists (E_MAX relative to 5 nM dioxin) of the AhR were 4-Me-indole (134%), 6-Me-indole (91%), and 7-MeO-indole (80%), respectively. The most effective antagonists of the AhR included 3-Me-indole (IC₅₀; 19 μM), 2,3-diMe-indole (IC₅₀; 11 μM), and 2,3,7-triMe-indole (IC₅₀; 12 μM). Reverse transcription polymerase chain reaction analyses of CYP1A1 mRNA in LS180 cells confirmed the data from gene reporter assays. The compound leads, 4-Me-indole and 7-MeO-indole, induced substantial nuclear translocation of the AhR and enriched binding of the AhR to the CYP1A1 promoter, as observed using fluorescent immunohistochemistry and chromatin immunoprecipitation assays, respectively. Molecular modeling and docking studies suggest the agonists and antagonists likely share the same binding pocket but have unique binding modes that code for their affinity. Binding pocket analysis further revealed that 4-methylindole and 7-methoxyindole can simultaneously bind to the pocket and produce synergistic interactions. Together, these data show a dependence on subtle and specific chemical indole structures as AhR modulators and furthermore underscore the importance of complete evaluation of indole compounds as nuclear receptor ligands.

Identification of potential aryl hydrocarbon receptor ligands by virtual screening of industrial chemicals

We have developed a virtual screening procedure to identify potential ligands to the aryl hydrocarbon receptor (AhR) among a set of industrial chemicals. AhR is a key target for dioxin-like compounds, which is related to these compounds' potential to induce cancer and a wide range of endocrine and immune system-related effects. The virtual screening procedure included an initial filtration aiming at identifying chemicals with structural similarities to 66 known AhR binders, followed by 3 enrichment methods run in parallel. These include two ligand-based methods (structural fingerprints and nearest neighbor analysis) and one structure-based method using an AhR homology model. A set of 6445 commonly used industrial chemicals was processed, and each step identified unique potential ligands. Seven compounds were identified by all three enrichment methods, and these compounds included known activators and suppressors of AhR. Only approximately 0.7% (41 compounds) of the studied industrial compounds was identified as potential AhR ligands and among these, 28 compounds have to our knowledge not been tested for AhR-mediated effects or have been screened with low purity. We suggest assessment of AhR-related activities of these compounds and in particular 2-chlorotrityl chloride, 3-p-hydroxyanilino-carbazole, and 3-(2-chloro-4-nitrophenyl)-5-(1,1-dimethylethyl)-1,3,4-oxadiazol-2(3H)-one.

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.

Molecular modeling of the AhR structure and interactions can shed light on ligand-dependent activation and transformation mechanisms

Molecular modeling has given important contributions to elucidation of the main stages in the AhR signal transduction pathway. Despite the lack of experimentally determined structures of the AhR functional domains, information derived from homologous systems has been exploited for modeling their structure and interactions. Homology models of the AhR PASB domain have provided information on the binding cavity and contributed to elucidate species-specific differences in ligand binding. Molecular Docking simulations of the ligand binding process have given insights into differences in binding of diverse agonists, antagonists, and selective AhR modulators, and their application to virtual screening of large databases of compounds have allowed identification of novel AhR ligands. Recently available structural information on protein-protein and protein-DNA complexes of other bHLH-PAS systems has opened the way for modeling the AhR:ARNT dimer structure and investigating the mechanisms of AhR transformation and DNA binding. Future research directions should include simulation of the protein dynamics to obtain a more reliable description of intermolecular interactions involved in signal transmission.

Docking Methodologies and Recent Advances

Docking, a molecular modelling method, has wide applications in identification and optimization in modern drug discovery. This chapter addresses the recent advances in the docking methodologies like fragment docking, covalent docking, inverse docking, post processing, hybrid techniques, homology modeling etc. and its protocol like searching and scoring functions. Advances in scoring functions for e.g. consensus scoring, quantum mechanics methods, clustering and entropy based methods, fingerprinting, etc. are used to overcome the limitations of the commonly used force-field, empirical and knowledge based scoring functions. It will cover crucial necessities and different algorithms of docking and scoring. Further different aspects like protein flexibility, ligand sampling and flexibility, and the performance of scoring function will be discussed.

Editor's Highlight: Microbial-Derived 1,4-Dihydroxy-2-naphthoic Acid and Related Compounds as Aryl Hydrocarbon Receptor Agonists/Antagonists: Structure-Activity Relationships and Receptor Modeling

1,4-Dihydroxy-2-naphthoic acid (1,4-DHNA) is a bacterial-derived metabolite that binds the aryl hydrocarbon receptor (AhR) and exhibits anti-inflammatory activity in the gut. The structure-dependent AhR activity of hydroxyl/carboxy-substituted naphthoic acids (NAs) was determined in young adult mouse colonic (YAMC) cells and human Caco2 colon cancer cells using CYP1A1/CYP1B1 mRNAs as Ah-responsive genes. Compounds used in this study include 1,4-, 3,5-, and 3,7-DHNA, 1,4-dimethoxy-2-naphthoic acid (1,4-DMNA), 1- and 4-hydroxy-2-naphthoic acid (1-HNA, 4-HNA), 1- and 2-naphthoic acid (1-NA, 2-NA), and 1- and 2-naphthol (1-NOH, 2-NOH). 1,4-DHNA was the most potent compound among hydroxyl/carboxy naphthalene derivatives, and the fold induction response for CYP1A1 and CYP1B1 was similar to that observed for 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) in YAMC and Caco2 cells. 1- and 4-HNA were less potent than 1,4-DHNA but induced maximal (TCDD-like) response for CYP1B1 (both cell lines) and CYP1A1 (Caco2 cells). With the exception of 1- and 2-NA, all compounds significantly induced Cyp1b1 in YAMC cells and these responses were not observed in AhR-deficient YAMC cells generated using CRISPR/Cas9 technology. In addition, we also observed that 1- and 2-NOH (and 1,4-DHNA) were weak AhR agonists, and 1- and 2-NOH also exhibited partial AhR antagonist activity. Structure-activity relationship studies for CYP1A1 but not CYP1B1 were similar in both cell lines, and CYP1A1 induction required one or both 1,4-dihydroxy substituents and activity was significantly enhanced by the 2-carboxyl group. We also used computational analysis to show that 1,4-DHNA and TCDD share similar interactions within the AhR binding pocket and differ primarily due to the negatively charged group of 1,4-DHNA.

In silico prediction and in vivo analysis of antiestrogenic potential of 2-isopropylthioxanthone (2-ITX) in juvenile goldfish (Carassius auratus)

Previous studies have shown both anti-estrogenic and anti-androgenic activities of 2-isopropylthioxanthone (2-ITX), a well known food contaminant, in in vitro assays. However, no data are available on the anti-estrogenic potentials and risks of 2-ITX in aquatic organisms. This work evaluated the potential endocrine disrupting effects of 2-ITX at the level of estrogen receptor (ER) signaling cascade using juvenile goldfish (Carassius auratus) as model. Firstly, we investigated the ligand binding efficiency of 2-ITX to the ligand binding domains (LBD) of goldfish ER subtypes using a molecular docking approach. Secondly, we assessed the effects of 2-ITX on E2-induced hepatic expression of ERα1, ERβ1, ERβ2, and vitellogenin (VTG) in vivo. Crosstalk between ER-VTG and aryl hydrocarbon receptor 2 (AhR2)-cytochrome P4501A (CYP1A) was also investigated. Fish were injected with increasing doses of 2-ITX ranging from 2 to 10µg/g BW, and results were compared to the effect of tamoxifen, a well-known ER modulator. We observed that compared to ERβ, the interaction potentials of 2-ITX to goldfish ERα1 LBD was more stable in the inactive receptor conformation. The in silico docking simulation analysis also revealed that 2-ITX acted as agonist for the goldfish AhR2 LBDs suggesting the ability of this compound to activate the cross-talk between the ERα- and AhR-signaling pathways. In vivo experiments confirm in silico simulation predictions demonstrating that 2-ITX reduced the estrogenicity of E2 at both transcriptional and post-transcriptional levels, indicating a clear anti-estrogenic effect. Co-exposure of E2 and 2-ITX also resulted in a significant decrease of CYP1A gene expression with respect to 2-ITX alone. Results from these studies collectively revealed that the antiestrogenic property of 2-ITX can be ascribed to a combination of effects on multiple signaling pathways suggesting the potential for this environmental contaminant to affect the hormonal control of reproductive processes in fish.

In vitro and in silico evaluation of transactivation potencies of avian AHR1 and AHR2 by endogenous ligands: Implications for the physiological role of avian AHR2

Aryl hydrocarbon receptor (AHR) is well conserved from invertebrates to vertebrates, and it mediates the toxic effects of exogenous ligands, including dioxins. Recent studies reported that AHRs activated by endogenous ligands play critical roles in mammalian physiological homeostasis. Avian species possess at least two AHR isoforms (AHR1 and AHR2), which exhibit species- and isoform-specific transactivation potencies to exogenous ligands, whereas mammals possess a single AHR. To delineate the profiles and roles of endogenous ligands for avian AHR isoforms, we investigated in vitro transactivation potencies of avian AHRs (AHR1 and AHR2 from the jungle crow, Corvus macrorhynchos; common cormorant, Phalacrocorax carbo; and black-footed albatross, Phoebastria nigripes) treated with the endogenous tryptophan metabolites 6-formylindolo [3,2-b] carbazole (FICZ), l-kynurenine (l-Kyn), kynurenic acid (KYNA), and indoxyl sulfate (IS). Furthermore, we analyzed the binding mode of these ligands to each avian AHR isoform by in silico docking simulations. The EC50 of FICZ (0.009-0.032nM) was similar regardless of the species or isoform of AHR. The estimated in silico binding mode of FICZ to AHRs was well conserved in both isoforms. The transactivation potencies of avian AHRs to other tryptophan metabolites were 10(5)-10(7) fold lower than those for FICZ, and EC50 values varied in a species- and isoform-specific manner. This was consistent with poor conservation of the binding mode of l-Kyn, KYNA, and IS predicted in in silico docking simulations. Our results suggest that in avian species, FICZ is the most potent endogenous AHR ligand, and that AHR1 and AHR2 are physiologically functional.

Deciphering Dimerization Modes of PAS Domains: Computational and Experimental Analyses of the AhR:ARNT Complex Reveal New Insights Into the Mechanisms of AhR Transformation

The Aryl hydrocarbon Receptor (AhR) is a transcription factor that mediates the biochemical response to xenobiotics and the toxic effects of a number of environmental contaminants, including dioxins. Recently, endogenous regulatory roles for the AhR in normal physiology and development have also been reported, thus extending the interest in understanding its molecular mechanisms of activation. Since dimerization with the AhR Nuclear Translocator (ARNT) protein, occurring through the Helix-Loop-Helix (HLH) and PER-ARNT-SIM (PAS) domains, is needed to convert the AhR into its transcriptionally active form, deciphering the AhR:ARNT dimerization mode would provide insights into the mechanisms of AhR transformation. Here we present homology models of the murine AhR:ARNT PAS domain dimer developed using recently available X-ray structures of other bHLH-PAS protein dimers. Due to the different reciprocal orientation and interaction surfaces in the different template dimers, two alternative models were developed for both the PAS-A and PAS-B dimers and they were characterized by combining a number of computational evaluations. Both well-established hot spot prediction methods and new approaches to analyze individual residue and residue-pairwise contributions to the MM-GBSA binding free energies were adopted to predict residues critical for dimer stabilization. On this basis, a mutagenesis strategy for both the murine AhR and ARNT proteins was designed and ligand-dependent DNA binding ability of the AhR:ARNT heterodimer mutants was evaluated. While functional analysis disfavored the HIF2α:ARNT heterodimer-based PAS-B model, most mutants derived from the CLOCK:BMAL1-based AhR:ARNT dimer models of both the PAS-A and the PAS-B dramatically decreased the levels of DNA binding, suggesting this latter model as the most suitable for describing AhR:ARNT dimerization. These novel results open new research directions focused at elucidating basic molecular mechanisms underlying the functional activity of the AhR.

Computational modeling combined with experimental validation may give insight into structural and functional properties of protein systems. The basic Helix-Loop-Helix PER-ARNT-SIM (bHLH-PAS) proteins show conserved functional domains despite the broad range of functions exerted by the different systems. Within this protein family, the Aryl hydrocarbon Receptor (AhR) is known to mediate the toxic effects of a number of environmental contaminants, including dioxins and dioxin-like chemicals, and it also exerts other biochemical and physiological effects. Despite the absence of experimentally determined structures, theoretical models of the AhR PAS domains developed on the basis of homologous systems have allowed understanding of some aspects of the molecular mechanisms underlying its function. In this work we present alternative structural models of the transcriptionally active complex of AhR with the AhR Nuclear Translocator (ARNT) protein. Computational characterization of the modeled protein-protein interaction interfaces guided the design of mutagenesis experiments, and evaluation of the DNA binding ability of the resulting AhR:ARNT dimer mutants allowed validation of the models and selection of the most reliable one. These findings open new research directions for understanding the molecular mechanisms underlying the functional activity of the AhR.

The tertiary structures of porcine AhR and ARNT proteins and molecular interactions within the TCDD/AhR/ARNT complex

The aryl hydrocarbon receptor (AhR) is a ligand-dependent transcription factor that can be activated by structurally diverse synthetic and natural chemicals, including toxic environmental contaminant 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). In the present study, homology models of the porcine AhR-ligand binding domain (LBD) and the porcine aryl hydrocarbon receptor nuclear translocator-ligand binding domain (ARNT-LBD) were created on the basis of structures of closely related respective proteins i.e., human Hif-2α and ARNT. Molecular docking of TCDD to the porcine AhR-LBD model revealed high binding affinity (-8.8kcal/mol) between TCDD and the receptor. Moreover, formation of the TCDD/AhR-LBD complex was confirmed experimentally with the use of electrophoretic mobility shift assay (EMSA). It was found that TCDD (10nM, 2h of incubation) not only bound to the AhR in the porcine granulosa cells but also activated the receptor. The current study provides a framework for examining the key events involved in the ligand-dependent activation of the AhR.

Access Path to the Ligand Binding Pocket May Play a Role in Xenobiotics Selection by AhR

Understanding of multidrug binding at the atomic level would facilitate drug design and strategies to modulate drug metabolism, including drug transport, oxidation, and conjugation. Therefore we explored the mechanism of promiscuous binding of small molecules by studying the ligand binding domain, the PAS-B domain of the aryl hydrocarbon receptor (AhR). Because of the low sequence identities of PAS domains to be used for homology modeling, structural features of the widely employed HIF-2α and a more recent suitable template, CLOCK were compared. These structures were used to build AhR PAS-B homology models. We performed molecular dynamics simulations to characterize dynamic properties of the PAS-B domain and the generated conformational ensembles were employed in in silico docking. In order to understand structural and ligand binding features we compared the stability and dynamics of the promiscuous AhR PAS-B to other PAS domains exhibiting specific interactions or no ligand binding function. Our exhaustive in silico binding studies, in which we dock a wide spectrum of ligand molecules to the conformational ensembles, suggest that ligand specificity and selection may be determined not only by the PAS-B domain itself, but also by other parts of AhR and its protein interacting partners. We propose that ligand binding pocket and access channels leading to the pocket play equally important roles in discrimination of endogenous molecules and xenobiotics.

Measured and predicted affinities of binding and relative potencies to activate the AhR of PAHs and their alkylated analogues

Polycyclic aromatic hydrocarbons (PAHs) and their alkylated forms are important components of crude oil. Both groups of PAHs have been reported to cause dioxin-like responses, mediated by aryl hydrocarbon receptor (AhR). Thus, characterization of binding affinity to the AhR of unsubstituted or alkylated PAHs is important to understand the toxicological consequences of oil contamination on ecosystems. We investigated the potencies of major PAHs of crude oil, e.g., chrysene, phenanthrene and dibenzothiophene, and their alkylated forms (n=17) to upregulate expression of AhR-mediated processes by use of the H4IIE-luc transactivation bioassay. In addition, molecular descriptors of different AhR activation potencies among PAHs were investigated by use of computational molecular docking models. Based on responses of the H4IIE-luc in vitro assay, it was shown that potencies of PAHs were determined by alkylation in addition to the number and conformation of rings. Potencies of AhR-mediated processes were generally greater when a chrysene group was substituted, especially in 1-methyl-chrysene. Significant negative correlations were observed between the in vitro dioxin-like potency measured in H4IIE-luc cells and the binding distance estimated from the in silico modeling. The difference in relative potency for AhR activation observed among PAHs and their alkylated forms could be explained by differences among binding distances in the ligand binding domain of the AhR caused by alkylation. The docking model developed in the present study may have utility in predicting risks of environmental contaminants of which toxicities are mediated by AhR binding.

Naturally occurring marine brominated indoles are aryl hydrocarbon receptor ligands/agonists

The aryl hydrocarbon receptor (AhR) is a ligand-dependent transcription factor that mediates the toxic and biological effects of structurally diverse chemicals, including the environmental contaminant 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). As part of a larger effort to identify the full spectrum of chemicals that can bind to and activate the AhR, we have examined the ability of several naturally occurring marine-derived brominated indoles and brominated (methylthio)indoles (collectively referred to as brominated indoles) to bind to the AhR and stimulate AhR-dependent gene expression. Incubation of mouse, rat, and guinea pig recombinant cell lines containing a stably transfected AhR-responsive luciferase reporter gene with eight brominated indoles revealed that all compounds stimulated luciferase reporter gene activity, although some species-specific differences were observed. All compounds induced significantly more luciferase activity when incubated with cells for 4 h as compared to 24 h, demonstrating that these compounds are transient activators of the AhR signaling pathway. Three of the brominated indoles induced CYP1A1 mRNA in human HepG2 cells in vitro and Cyp1a mRNA in zebrafish embryos in vivo. The identification of the brominated indoles as direct ligands and activators/agonists of the AhR was confirmed by their ability to compete with [(3)H]TCDD for binding to the AhR and to stimulate AhR transformation and DNA binding in vitro. Taken together, these results indicate that marine-derived brominated indoles are members of a new class of naturally occurring AhR agonists.

In silico analysis of the interaction of avian aryl hydrocarbon receptors and dioxins to decipher isoform-, ligand-, and species-specific activations

The aryl hydrocarbon receptor (AHR) mediates toxic responses to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and other dioxin-like compounds (DLCs). Avian species possess multiple AHR isoforms (AHR1, AHR1β, and AHR2) that exhibit species- and isoform-specific responses to ligands. To account for the ligand preference in terms of the structural features of avian AHRs, we generated in silico homology models of the ligand-binding domain of avian AHRs based on holo human HIF-2α (PDB entry 3H7W ). Molecular docking simulations of TCDD and other DLCs with avian AHR1s and AHR2s using ASEDock indicated that the interaction energy increased with the number of substituted chlorine atoms in congeners, supporting AHR transactivation potencies and World Health Organization TCDD toxic equivalency factors of congeners. The potential interaction energies of an endogenous AHR ligand, 6-formylindolo [3,2-b] carbazole (FICZ) to avian AHRs were lower than those of TCDD, which was supported by a greater potency of FICZ for in vitro AHR-mediated transactivation than TCDD. The molecular dynamics simulation revealed that mean square displacements in Ile324 and Ser380 of TCDD-bound AHR1 of the chicken, the most sensitive species to TCDD, were smaller than those in other avian AHR1s, suggesting that the dynamic stability of these amino acid residues contribute to TCDD preference. For avian AHR2, the corresponding residues (Val/Ser or Val/Ala type) were not responsible for differential TCDD sensitivity. Application of the three-dimensional reference interaction site model showed that the stabilization of TCDD binding to avian AHRs may be due to the solvation effect depending on the characteristics of two amino acids corresponding to Ile324 and Ser380 in chicken AHR1. This study demonstrates that in silico simulations of AHRs and ligands could be used to predict isoform-, ligand-, and species-specific interactions.

A Structural Switch between Agonist and Antagonist Bound Conformations for a Ligand-Optimized Model of the Human Aryl Hydrocarbon Receptor Ligand Binding Domain

The aryl hydrocarbon receptor (AHR) is a ligand-activated transcription factor that regulates the expression of a diverse group of genes. Exogenous AHR ligands include the environmental contaminant 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), which is a potent agonist, and the synthetic AHR antagonist N-2-(1H-indol-3yl)ethyl)-9-isopropyl-2-(5-methylpyridin-3-yl)-9H-purin-6-amine (GNF351). As no experimentally determined structure of the ligand binding domain exists, homology models have been utilized for virtual ligand screening (VLS) to search for novel ligands. Here, we have developed an “agonist-optimized” homology model of the human AHR ligand binding domain, and this model aided in the discovery of two human AHR agonists by VLS. In addition, we performed molecular dynamics simulations of an agonist TCDD-bound and antagonist GNF351-bound version of this model in order to gain insights into the mechanics of the AHR ligand-binding pocket. These simulations identified residues 307–329 as a flexible segment of the AHR ligand pocket that adopts discrete conformations upon agonist or antagonist binding. This flexible segment of the AHR may act as a structural switch that determines the agonist or antagonist activity of a given AHR ligand.

In silico identification of an aryl hydrocarbon receptor antagonist with biological activity in vitro and in vivo

The aryl hydrocarbon receptor (AHR) is critically involved in several physiologic processes, including cancer progression and multiple immune system activities. We, and others, have hypothesized that AHR modulators represent an important new class of targeted therapeutics. Here, ligand shape-based virtual modeling techniques were used to identify novel AHR ligands on the basis of previously identified chemotypes. Four structurally unique compounds were identified. One lead compound, 2-((2-(5-bromofuran-2-yl)-4-oxo-4H-chromen-3-yl)oxy)acetamide (CB7993113), was further tested for its ability to block three AHR-dependent biologic activities: triple-negative breast cancer cell invasion or migration in vitro and AHR ligand-induced bone marrow toxicity in vivo. CB7993113 directly bound both murine and human AHR and inhibited polycyclic aromatic hydrocarbon (PAH)- and TCDD-induced reporter activity by 75% and 90% respectively. A novel homology model, comprehensive agonist and inhibitor titration experiments, and AHR localization studies were consistent with competitive antagonism and blockade of nuclear translocation as the primary mechanism of action. CB7993113 (IC50 3.3 × 10(-7) M) effectively reduced invasion of human breast cancer cells in three-dimensional cultures and blocked tumor cell migration in two-dimensional cultures without significantly affecting cell viability or proliferation. Finally, CB7993113 effectively inhibited the bone marrow ablative effects of 7,12-dimethylbenz[a]anthracene in vivo, demonstrating drug absorption and tissue distribution leading to pharmacological efficacy. These experiments suggest that AHR antagonists such as CB7993113 may represent a new class of targeted therapeutics for immunomodulation and/or cancer therapy.

AhR sensing of bacterial pigments regulates antibacterial defence

The aryl hydrocarbon receptor (AhR) is a highly conserved ligand-dependent transcription factor that senses environmental toxins and endogenous ligands, thereby inducing detoxifying enzymes and modulating immune cell differentiation and responses. We hypothesized that AhR evolved to sense not only environmental pollutants but also microbial insults. We characterized bacterial pigmented virulence factors, namely the phenazines from Pseudomonas aeruginosa and the naphthoquinone phthiocol from Mycobacterium tuberculosis, as ligands of AhR. Upon ligand binding, AhR activation leads to virulence factor degradation and regulated cytokine and chemokine production. The relevance of AhR to host defence is underlined by heightened susceptibility of AhR-deficient mice to both P. aeruginosa and M. tuberculosis. Thus, we demonstrate that AhR senses distinct bacterial virulence factors and controls antibacterial responses, supporting a previously unidentified role for AhR as an intracellular pattern recognition receptor, and identify bacterial pigments as a new class of pathogen-associated molecular patterns.