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

Food-Borne Polycyclic Aromatic Hydrocarbons and Circadian Disruption

Circadian disruption has been found to increase the risk of metabolic diseases, brain disorders, and cancer. The aryl hydrocarbon receptor (AhR), responsible for xenobiotic metabolism, is known to be activated by certain environmental stimuli, including polycyclic aromatic hydrocarbons (PAHs). Exposure to these stimuli may lead to diseases related to circadian disruption, with AhR activation suggested as a leading cause. Both the aryl hydrocarbon receptor nuclear translocator (ARNT) and aryl hydrocarbon receptor nuclear translocator-like (BMAL1) are class II basic helix-loop-helix/Per-ARNT-SIM (bHLH-PAS) proteins. These proteins form heterodimers with stimulated class I bHLH-PAS proteins, including circadian locomotor output cycles kaput (CLOCK) and AhR. Due to their sequential similarity, the overactivation of AhR by toxicants, such as PAHs, may lead to the formation of heterodimers with BMAL1, potentially causing circadian disruption. Dysregulation of BMAL1 can affect a wide range of metabolic genes, emphasizing its crucial roles. However, this issue has not been adequately addressed. Previous studies have reported that the inhibitory effects of phytochemicals on AhR activation can ameliorate diseases induced by environmental toxicants. Additionally, some phytochemicals have shown preventive effects on circadian misalignment. Therefore, this Review aims to explore potential strategies to prevent circadian disruption induced by food-borne toxicants, such as benzo[a]pyrene; to generate new ideas for future studies; and to highlight the importance of investigating these preventive strategies.

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.

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.

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.

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.

Current Therapeutic Landscape and Safety Roadmap for Targeting the Aryl Hydrocarbon Receptor in Inflammatory Gastrointestinal Indications

Target modulation of the AhR for inflammatory gastrointestinal (GI) conditions holds great promise but also the potential for safety liabilities both within and beyond the GI tract. The ubiquitous expression of the AhR across mammalian tissues coupled with its role in diverse signaling pathways makes development of a “clean” AhR therapeutically challenging. Ligand promiscuity and diversity in context-specific AhR activation further complicates targeting the AhR for drug development due to limitations surrounding clinical translatability. Despite these concerns, several approaches to target the AhR have been explored such as small molecules, microbials, PROTACs, and oligonucleotide-based approaches. These various chemical modalities are not without safety liabilities and require unique de-risking strategies to parse out toxicities. Collectively, these programs can benefit from in silico and in vitro methodologies that investigate specific AhR pathway activation and have the potential to implement thresholding parameters to categorize AhR ligands as “high” or “low” risk for sustained AhR activation. Exploration into transcriptomic signatures for AhR safety assessment, incorporation of physiologically-relevant in vitro model systems, and investigation into chronic activation of the AhR by structurally diverse ligands will help address gaps in our understanding regarding AhR-dependent toxicities. Here, we review the role of the AhR within the GI tract, novel therapeutic modality approaches to target the AhR, key AhR-dependent safety liabilities, and relevant strategies that can be implemented to address drug safety concerns. Together, this review discusses the emerging therapeutic landscape of modalities targeting the AhR for inflammatory GI indications and offers a safety roadmap for AhR drug development.

The aryl hydrocarbon receptor: A diagnostic and therapeutic target in glioma

Glioblastoma multiforme (GBM) is a deadly disease; 5-year survival rates have shown little improvement over the past 30 years. In vivo positron emission tomography (PET) imaging is an important method of identifying potential diagnostic and therapeutic molecular targets non-invasively. The aryl hydrocarbon receptor (AhR) is a transcription factor that regulates multiple genes involved in immune response modulation and tumorigenesis. The AhR is an attractive potential drug target and studies have shown that its activation by small molecules can modulate innate and adaptive immunity beneficially and prevent AhR-mediated tumour promotion in several cancer types. In this review, we provide an overview of the role of the AhR in glioma tumorigenesis and highlight its potential as an emerging biomarker for glioma therapies targeting the tumour immune response and PET diagnostics.

Evolution of Hominin Detoxification: Neanderthal and Modern Human Ah Receptor Respond Similarly to TCDD

In studies of hominin adaptations to fire use, the role of the aryl hydrocarbon receptor (AHR) in the evolution of detoxification has been highlighted, including statements that the modern human AHR confers a significantly better capacity to deal with toxic smoke components than the Neanderthal AHR. To evaluate this, we compared the AHR-controlled induction of cytochrome P4501A1 (CYP1A1) mRNA in HeLa human cervix epithelial adenocarcinoma cells transfected with an Altai-Neanderthal or a modern human reference AHR expression construct, and exposed to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). We compared the complete AHR mRNA sequences including the untranslated regions (UTRs), maintaining the original codon usage. We observe no significant difference in CYP1A1 induction by TCDD between Neanderthal and modern human AHR, whereas a 150–1,000 times difference was previously reported in a study of the AHR coding region optimized for mammalian codon usage and expressed in rat cells. Our study exemplifies that expression in a homologous cellular background is of major importance to determine (ancient) protein activity. The Neanderthal and modern human dose–response curves almost coincide, except for a slightly higher extrapolated maximum for the Neanderthal AHR, possibly caused by a 5′-UTR G-variant known from modern humans (rs7796976). Our results are strongly at odds with a major role of the modern human AHR in the evolution of hominin detoxification of smoke components and consistent with our previous study based on 18 relevant genes in addition to AHR, which concluded that efficient detoxification alleles are more dominant in ancient hominins, chimpanzees, and gorillas than in modern humans.

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.

Targeting the aryl hydrocarbon receptor with a novel set of triarylmethanes

The aryl hydrocarbon receptor (AhR) is a chemical sensor upregulating the transcription of responsive genes associated with endocrine homeostasis, oxidative balance and diverse metabolic, immunological and inflammatory processes, which have raised the pharmacological interest on its modulation. Herein, a novel set of 32 unsymmetrical triarylmethane (TAM) class of structures has been synthesized, characterized and their AhR transcriptional activity evaluated using a cell-based assay. Eight of the assayed TAM compounds (14, 15, 18, 19, 21, 22, 25, 28) exhibited AhR agonism but none of them showed antagonist effects. TAMs bearing benzotrifluoride, naphthol or heteroaromatic (indole, quinoline or thiophene) rings seem to be prone to AhR activation unlike phenyl substituted or benzotriazole derivatives. A molecular docking analysis with the AhR ligand binding domain (LBD) showed similarities in the binding mode and in the interactions of the most potent TAM identified 4-(pyridin-2-yl (thiophen-2-yl)methyl)phenol (22) compared to the endogenous AhR agonist 5,11-dihydroindolo[3,2-b]carbazole-12-carbaldehyde (FICZ). Finally, in silico predictions of physicochemical and biopharmaceutical properties for the most potent agonistic compounds were performed and these exhibited acceptable druglikeness and good ADME profiles. To our knowledge, this is the first study assessing the AhR modulatory effects of unsymmetrical TAM class of compounds.

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).

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.

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.

Complex immune responses and molecular reactions to pathogens and disease in a desert reptile (Gopherus agassizii)

Immune function plays an important role in an animal's defense against infectious disease. In reptiles, immune responses may be complex and counterintuitive, and diagnostic tools used to identify infection, such as induced antibody responses are limited. Recent studies using gene transcription profiling in tortoises have proven useful in identifying immune responses to various intrinsic and extrinsic stressors. As part of a larger experiment with Mojave desert tortoises (Gopherus agassizii), we facilitated the transmission of the pathogenic bacteria, Mycoplasma agassizii (Myag), to naïve adults and measured innate and induced immune reactions over time. Specifically, we evaluated clinical condition, presence of Myag in the nasal/oral cavity, induced antibody responses specific to Myag, and measured molecular reactions (gene transcript profiles) in 15 captive tortoises classified as naïve, exposed, or infected and 14 wild tortoises for comparison. Myag was confirmed inside the nasal/oral cavity in exposed tortoises within 30-60 days of introduction to infected animals, yet we did not detect Myag specific induced antibody responses in these individuals until 420-595 days post exposure. Surprisingly, we found no overall differences in the gene transcript profiles between our experimental treatment groups throughout this study. This work highlights the complexities in assessing immune function and diagnosing pathogen related infections in tortoises and other reptiles.

Identification of the novel role of butyrate as AhR ligand in human intestinal epithelial cells

The ligand activated transcription factor, aryl hydrocarbon receptor (AhR) emerged as a critical regulator of immune and metabolic processes in the gastrointestinal tract. In the gut, a main source of AhR ligands derives from commensal bacteria. However, many of the reported microbiota-derived ligands have been restricted to indolyl metabolites. Here, by screening commensal bacteria supernatants on an AhR reporter system expressed in human intestinal epithelial cell line (IEC), we found that the short chain fatty acid (SCFA) butyrate induced AhR activity and the transcription of AhR-dependent genes in IECs. We showed that AhR ligand antagonists reduced the effects of butyrate on IEC suggesting that butyrate could act as a ligand of AhR, which was supported by the nuclear translocation of AhR induced by butyrate and in silico structural modelling. In conclusion, our findings suggest that (i) butyrate activates AhR pathway and AhR-dependent genes in human intestinal epithelial cell-lines (ii) butyrate is a potential ligand for AhR which is an original mechanism of gene regulation by SCFA.

Comparative analysis of interactions between aryl hydrocarbon receptor ligand binding domain with its ligands: a computational study

BACKGROUND

Aryl hydrocarbon receptor (AhR) ligands may act as potential carcinogens or anti-tumor agents. Understanding how some of the residues in AhR ligand binding domain (AhRLBD) modulate their interactions with ligands would be useful in assessing their divergent roles including toxic and beneficial effects. To this end, we have analysed the nature of AhRLBD interactions with 2,3,7,8-tetrachlorodibenzo-ρ-dioxin (TCDD), 6-formylindolo[3,2-b]carbazole (FICZ), indole-3-carbinol (I3C) and its degradation product, 3,3'-diindolylmethane (DIM), Resveratrol (RES) and its analogue, Piceatannol (PTL) using molecular modeling approach followed by molecular dynamic simulations.

RESULTS

Results showed that each of the AhR ligands, TCDD, FICZ, I3C, DIM, RES and PTL affect the local and global conformations of AhRLBD.

CONCLUSION

The data presented in this study provide a structural understanding of AhR with its ligands and set the basis for its functions in several pathways and their related diseases.

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.

Ligand-induced perturbation of the HIF-2α:ARNT dimer dynamics

Hypoxia inducible factors (HIFs) are transcription factors belonging to the basic helix−loop−helix PER-ARNT-SIM (bHLH-PAS) protein family with a role in sensing oxygen levels in the cell. Under hypoxia, the HIF-α degradation pathway is blocked and dimerization with the aryl hydrocarbon receptor nuclear translocator (ARNT) makes HIF-α transcriptionally active. Due to the common hypoxic environment of tumors, inhibition of this mechanism by destabilization of HIF-α:ARNT dimerization has been proposed as a promising therapeutic strategy. Following the discovery of a druggable cavity within the PAS-B domain of HIF-2α, research efforts have been directed to identify artificial ligands that can impair heterodimerization. Although the crystallographic structures of the HIF-2α:ARNT complex have elucidated the dimer architecture and the 0X3-inhibitor placement within the HIF-2α PAS-B, unveiling the inhibition mechanism requires investigation of how ligand-induced perturbations could dynamically propagate through the structure and affect dimerization. To this end, we compared evolutionary features, intrinsic dynamics and energetic properties of the dimerization interfaces of HIF-2α:ARNT in both the apo and holo forms. Residue conservation analysis highlighted inter-domain connecting elements that have a role in dimerization. Analysis of domain contributions to the dimerization energy demonstrated the importance of bHLH and PAS-A of both partners and of HIF-2α PAS-B domain in dimer stabilization. Among quaternary structure oscillations revealed by Molecular Dynamics simulations, the hinge-bending motion of the ARNT PAS-B domain around the flexible PAS-A/PAS-B linker supports a general model for ARNT dimerization in different heterodimers. Comparison of the HIF-2α:ARNT dynamics in the apo and 0X3-bound forms indicated a model of inhibition where the HIF-2α-PAS-B interfaces are destabilised as a result of water-bridged ligand-protein interactions and these local effects allosterically propagate to perturb the correlated motions of the domains and inter-domain communication. These findings will guide the design of improved inhibitors to contrast cell survival in tumor masses.

A low oxygen condition, called hypoxia, often occurs in tumor masses and generally correlates with worse prognosis. Cells in a tumor react to low oxygen levels with a metabolism modification induced by the activation of hypoxia inducible factors (HIFs) through dimerization with a partner protein and binding to a DNA target. Disrupting this protein-protein interaction could be a potential therapeutic strategy, but directly interfering with dimer formation can be troublesome because of the difficulty to design drugs that bind to protein interfaces. However, ligands that bind internal protein cavities can indirectly perturb the interfaces reducing dimers stability. Albeit protein crystallography had offered a detailed static picture of a HIF dimer bound to candidate inhibitors, it is not able to describe either the perturbation caused by binding or the molecular mechanism of dimer destabilization. Here we exploit molecular dynamics to identify the crucial interfaces in the HIF dimer stabilization and, by comparing the results obtained in the bound and unbound forms, we reveal the mechanism of ligand inhibition at atomic detail. All these findings will guide toward the design of improved dimerization inhibitors, to contrast cell survival in tumor masses.

Structural modeling of the AhR:ARNT complex in the bHLH-PASA-PASB region elucidates the key determinants of dimerization

Elucidation of the dimerization process of the aryl hydrocarbon receptor (AhR) with the AhR nuclear translocator (ARNT) is crucial for understanding the mechanisms underlying the functional activity of AhR, including mediation of the toxicity of environmental contaminants. In this work, for the first time a structural model of the AhR:ARNT dimer encompassing the entire bHLH-PASA-PASB domain region is proposed. It is developed by using a template-based modeling approach, relying on the recently available crystallographic structures of two dimers of homologous systems in the bHLH-PAS family of proteins: the CLOCK:BMAL1 and the HIF2α:ARNT heterodimers. The structural and energetic characteristics of the modeled AhR:ARNT protein-protein interface are determined by evaluating the variations in solvent accessible surface area, the total binding free energy and the per-residue free energy contributions obtained by the MM-GBSA method and the Energy Decomposition Analysis. The analyses of the intricate network of inter-domain interactions at the dimerization interfaces provide insights into the key determinants of dimerization. These are confirmed by comparison of the computational findings with the available experimental mutagenesis and functional analysis data. The results presented here on the AhR:ARNT dimer structure and interactions provide a framework to start analyzing the mechanism of AhR transformation into its functional DNA binding form.