The aryl hydrocarbon receptor (AHR) mediates toxic responses to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and other dioxin-like cornpounds (DLCs). Avian species possess multiple AHR isoforms (AHR1, AHR1,6, 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 ho/o human HIF-2 alpha (PDB entry 3H7VV). Molecular docking simulations of TCDD and other DLCs with avian AHRls 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 silica simulations of AHRs and ligands could be used to predict isoform-, ligand-, and species-specific interactions.
