Theoretical investigation of AhR binding property with relevant structural requirements for AhR-mediated toxicity of polybrominated diphenyl ethers.

Polybrominated diphenyl ethers (PBDEs) are more frequently suspected with the induction of toxicity via signal transduction pathway of cytosolic aryl hydrocarbon receptor (AhR), the initial binding to which is assumed to be an essential prerequisite during the ligand-dependent activation. However, the AhR binding property and associated toxicity of PBDEs is yet to be clearly known for lacking insights into the structural requirements at molecular level. To understand the AhR binding property of PBDEs, the ligand binding domain (LBD) of AhR was simulatively developed on homologous protein after basic validation of geometrical rationality and the binding interaction profile was visually described using molecular docking approach. For AhR binding, the offset or edge-on π-π stackings with aromatic motifs including Phe289, Phe345 and His285 were shown to be structurally required whereas the electrostatic attraction validated for AhR binding to dioxins might be less effective for 2,2',3,4,4'-pentabromodiphenyl ether (BDE-85). Besides the demands of less steric hindrance from alanines and weak formulation of hydrogen bonds, the dispersion force through large contact and polarization of S-π electrons seemed to be impactful when BDE-85 were closer to Cys327, Met334 or Met342. With theoretical computation of AhR binding energies, the more significant correlativity with bioassays was derived especially for the lowly/moderately brominated congeners, and could be used to predict the AhR binding affinity on certain degree. The informative results would thus not only help well understand the molecular basis of AhR-mediated toxicity but give an approach for accelerative evaluation of AhR binding and toxicity of PBDEs.