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

Dimerization Rules of Mammalian PAS Proteins

The PAS (PER, ARNT, SIM) protein family plays a vital role in mammalian biology and human disease. This analysis arose from an interest in the signaling mechanics by the Ah receptor (AHR) and the Ah receptor nuclear translocator (ARNT). After more than fifty years by studying this and related mammalian sensor systems, describing the role of PAS domains in signal transduction is still challenging. In this perspective, we attempt to interpret recent studies of mammalian PAS protein structure and consider how this new insight might explain how these domains are employed in human signal transduction with an eye towards developing strategies to target and engineer these molecules for a new generation of therapeutics. Our approach is to integrate our understanding of PAS protein history, cell biology, and molecular biology with recent structural discoveries to help explain the mechanics of mammalian PAS protein signaling. As a learning set, we focus on sequences and crystal structures of mammalian PAS protein dimers that can be visualized using readily available software.

Promotion of ABCG2 gene expression by neolignans from Piper longum L

We focused on Piper longum L., a herbal drug produced in Myanmar, which has a renoprotective effect. Thus, we attempted to isolate and identify compounds that enhance the expression of the ABCG2 gene from the aerial parts of the plant except for the fruit. Among the various P. longum extracts, we isolated and identified the components. Using Caco-2 cells, the hABCG2 mRNA expression-enhancing effects of the isolated compounds were compared with the positive reference compound (3-methylcholanthrene [3MC]) using real-time polymerase chain reaction. Six compounds were isolated and identified from the methanol extract of P. longum. Among the isolated compounds, licarin A and neopomatene had lower toxicity and higher hABCG2 mRNA expression-enhancing effects in Caco-2 cells. Suppression of hAhR expression by siRNA reduced the activity of licarin A and neopomatene, as well as the hAhR agonist 3MC, suggesting that these 2 compounds may act as hAhR agonists to promote hABCG2 expression.

Molecular dynamics simulations depict structural motions of the whole human aryl hydrocarbon receptor influencing its binding of ligands and HSP90

The aryl hydrocarbon receptor (AhR) has broad biological functions when its ligands activate it; the non-binding interactions with AhR have not been fully elucidated due to the absence of a complete tridimensional (3D) structure. Therefore, utilization of the whole 3D structure from Homo sapiens AhR by in silico studies will allow us to better study and analyze the binding mode of its full and partial agonists, and antagonists, as well as its interaction with the HSP90 chaperone. The 3D AhR structure was obtained from I-TASSER and subjected to molecular dynamics (MD) simulations to obtain different structural conformations and determine the most populated AhR conformer by clustering analyses. The AhR-3D structures selected from MD simulations and those from clustering analyses were used to achieve docking studies with some of its ligands and protein-protein docking with HSP90. Once the AhR-3D structure was built, its Ramachandran maps and energy showed a well-qualified 3D model. MD simulations showed that the per-Arnt-Sim homology (PAS) PAS A, PAS B, and Q domains underwent conformational changes, identifying the conformation when agonists were binding also, and HSP90 was binding near the PAS A, PAS B, and Q domains. However, when antagonists are binding, HSP90 does not bind near the PAS A, PAS B, and Q domains. These studies show that the complex agonist-AhR-HSP90 can be formed, but this complex is not formed when an antagonist is binding. Knowing the conformations when the ligands bind to AHR and the behavior of HSP90 allows for an understanding of its activity.Communicated by Ramaswamy H. Sarma.

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.

Key Amino Acids for Transferase Activity of GDSL Lipases

The Gly-Asp-Ser-Leu (GDSL) motif of esterase/lipase family proteins (GELPs) generally exhibit esterase activity, whereas transferase activity is markedly preferred in several GELPs, including the Tanacetum cinerariifolium GDSL lipase TciGLIP, which is responsible for the biosynthesis of the natural insecticide, pyrethrin I. This transferase activity is due to the substrate affinity regulated by the protein structure and these features are expected to be conserved in transferase activity-exhibiting GELPs (tr-GELPs). In this study, we identified two amino acid residues, [N/R]208 and D484, in GELP sequence alignments as candidate key residues for the transferase activity of tr-GELPs by two-entropy analysis. Molecular phylogenetic analysis demonstrated that each tr-GELP is located in the clusters for non-tr-GELPs, and most GELPs conserve at least one of the two residues. These results suggest that the two conserved residues are required for the acquisition of transferase activity in the GELP family. Furthermore, substrate docking analyses using ColabFold-generated structure models of both natives and each of the two amino acids-mutated TciGLIPs also revealed numerous docking models for the proper access of substrates to the active site, indicating crucial roles of these residues of TciGLIP in its transferase activity. This is the first report on essential residues in tr-GELPs for the transferase activity.

Molecular and structural basis of interactions of vitamin D3 hydroxyderivatives with aryl hydrocarbon receptor (AhR): An integrated experimental and computational study

To better understand the molecular and structural basis underlying the interaction of vitamin D3 hydroxyderivatives with AhR, molecular simulation was used to probe the binding of 1,20(OH)₂D3, 1,25(OH)₂D3, 20,23(OH)₂D3 and 20(OH)D3 to AhR. qPCR showed that vitamin D3 derivatives stimulate expression of cyp1A1 and cyp1B1 genes that are downstream targets of AhR signaling. These secosteroids stimulated the translocation of the AhR to the nucleus, as measured by flow cytometry and western blotting. Molecular dynamics simulations were used to model the binding of vitamin D3 derivatives to AhR to examine their influence on the structure, conformation and dynamics of the AhR ligand binding domain (LBD). Binding thermodynamics, conformation, secondary structure, dynamical motion and electrostatic potential of AhR were analyzed. The molecular docking scores and binding free energy were all favorable for the binding of D3 derivatives to the AhR. These established ligands and the D3 derivatives are predicted to have different patterns of hydrogen bond formation with the AhR, and varied residue conformational fluctuations and dynamical motion for the LBD. These changes could alter the shape, size and electrostatic potential distribution of the ligand binding pocket, contributing to the different binding affinities of AhR for the natural ligands and D3 derivatives.

Elucidating the aryl hydrocarbon receptor antagonism from a chemical-structural perspective

The aryl hydrocarbon receptor (AhR) plays an important role in several biological processes such as reproduction, immunity and homoeostasis. However, little is known on the chemical-structural and physicochemical features that influence the activity of AhR antagonistic modulators. In the present report, in vitro AhR antagonistic activity evaluations, based on a chemical-activated luciferase gene expression (AhR-CALUX) bioassay, and an extensive literature review were performed with the aim of constructing a structurally diverse database of contaminants and potentially toxic chemicals. Subsequently, QSAR models based on Linear Discriminant Analysis and Logistic Regression, as well as two toxicophoric hypotheses were proposed to model the AhR antagonistic activity of the built dataset. The QSAR models were rigorously validated yielding satisfactory performance for all classification parameters. Likewise, the toxicophoric hypotheses were validated using a diverse set of 350 decoys, demonstrating adequate robustness and predictive power. Chemical interpretations of both the QSAR and toxicophoric models suggested that hydrophobic constraints, the presence of aromatic rings and electron-acceptor moieties are critical for the AhR antagonism. Therefore, it is hoped that the deductions obtained in the present study will contribute to elucidate further on the structural and physicochemical factors influencing the AhR antagonistic activity of chemical compounds.

Computational prediction of cytochrome P450 inhibition and induction

Cytochrome P450 (CYP) enzymes play an important role in the phase I metabolism of many xenobiotics. Most drug-drug interactions (DDIs) associated with CYP are caused by either CYP inhibition or induction. The early detection of potential DDIs is highly desirable in the pharmaceutical industry because DDIs can cause serious adverse events, which can lead to poor patient health and drug development failures. Recently, many computational studies predicting CYP inhibition and induction have been reported. The current computational modeling approaches for CYP metabolism are classified as ligand- and structure-based; various techniques, such as quantitative structure-activity relationships, machine learning, docking, and molecular dynamic simulation, are involved in both the approaches. Recently, combining these two approaches have resulted in improvements in the prediction accuracy of DDIs. In this review, we present important, recent developments in the computational prediction of the inhibition of four clinically crucial CYP isoforms (CYP1A2, 2C9, 2D6, and 3A4) and three nuclear receptors (aryl hydrocarbon receptor, constitutive androstane receptor, and pregnane X receptor) involved in the induction of CYP1A2, 2B6, and 3A4, respectively.

Druggability assessment of mammalian Per-Arnt-Sim [PAS] domains using computational approaches

Per-Arnt-Sim (PAS) domains are key regions that occur in different regulatory proteins from all kingdoms of life. PAS domains show a remarkably conserved structural scaffold, despite a highly variable primary sequence. In this study we have attempted to address some of the gaps in knowledge regarding the druggability of PAS-A domains, differences in structure and dynamics within the PAS domain family and how this affects the druggability potential, as well as give insight into the druggability of steroid receptor coactivators and putative binding modes of the NCOA1. Investigations were performed through a range of computational methods including molecular docking studies, atomistic molecular dynamics simulations, and hotspot mapping. Atomistic molecular dynamics simulations show that the function of the AhR PAS-B domain is regulated by the dynamics of the highly conserved tyrosine Y322 residue, which acts as a "gatekeeper" controlling the access to the binding cavity and finely tuning the binding affinity. Furthermore, the transition between the partially unfolded and helical conformation of the loop1 segment within PAS-B domains was shown to be essential for the generation of "druggable" sites, especially for the NCOA1 PAS-B domain. Finally, our simulations indicated the undruggability of PAS-A domains, caused by the inherent characteristics of their putative binding sites. In conclusion, this work emphasises the role of intrinsic dynamics in tuning the druggability of PAS-B domains and shows that PAS-B domains of steroid receptor coactivators, such as NCOA1, can be targeted by small molecule ligands, which highlights the potential of developing new therapeutics designed to target these coactivators using structure-based approaches.

Bioassay system for the detection of aryl hydrocarbon receptor agonists in waterborne pesticides using zebrafish cyp1a1 promoter-luciferase recombinant hepatic cells

Currently used pesticides are widely distributed in surface water. Most of them are harmful to aquatic animals. In the present study, a novel recombinant hepatic cell line was developed for detecting the activity of aryl hydrocarbon receptor (AhR) agonists, and their potential bio-toxicity to aquatic animals was evaluated. First, HepG2 cells were stably transfected with a luciferase reporter plasmid containing a zebrafish (Danio rerio) cyp1a1 promoter with twelve copies of dioxin-responsive elements. The minimal detection limit was a 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) concentration of 0.3 ng/L (equivalent to 0.93 pM). The AhR agonist activity of thirty-seven pesticides was screened by measuring luciferase expression in the recombinant cells. Four pesticides (ipconazole, mepanipyrim, quinoxyfen and pencycuron) significantly induced luciferase expression. Additionally, cells treated with these four pesticides increased the cyp1a1 mRNA expression level, which further confirmed that they are AhR agonists. In conclusion, for the first time, we developed recombinant hepatic cells with a luciferase reporter plasmid containing a modified zebrafish cyp1a1 promoter. The cyp1a1 promoter luciferase reporter cells can be used as a sensitive and specific model to screen AhR agonists among currently used pesticides in water, which could be beneficial for risk assessment in aquatic environments.

Mechanisms and therapeutic prospects of polyphenols as modulators of the aryl hydrocarbon receptor

Polyphenolic AhR modulators displayed concentration-, XRE-, gene-, species- and cell-specific agonistic/antagonistic activity.

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.