Comprehensive binding analysis of polybrominated diphenyl ethers and aryl hydrocarbon receptor via an integrated...

Polybrominated diphenyl ethers (PBDEs) are often suspected to activate the signal transduction pathway of aryl hydrocarbon receptor (AhR), a ligand-activated transcription factor, for the induction of toxicity. Hence, the binding property of PBDEs with AhR is assumed to be associated with the ligand-dependent activation of AhR that may introduce many drug-metabolizing enzymes of genes encoding. However, the binding mechanism and the structural effect of PBDEs on their binding properties of AhR still need to be unraveled for toxicology research. A comprehensive study of the PBDEs-AhR binding mechanism was investigated using an integrated molecular modeling approach with two-dimensional quantitative structure-activity relationships (2D-QSAR), three-dimensional QSAR (3D-QSAR), and molecular docking simulation. Molecular docking revealed the differences in binding domains among 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD)-AhR complex and two PBDE-AhR complexes. A 2D-QSAR model was developed to analyze the overall structural effects of PBDEs on the binding affinity of AhR. It provided an insight into major physico-chemical properties by multiple linear regression based on genetic algorithm with reasonable results. The 3D-QSAR modeling discovered the detailed interaction features of binding sites, configurations and interaction fields of AhR with different PBDE ligands. This study demonstrated that the descriptors of Smin69 and MoRSEC15 were related to electronic properties and had a great effect on the relative binding affinities. The position of Br substitutions exhibited a significant influence on the interactions between AhR and PBDEs, including halogen interaction, π-S interaction, π-π stacking interaction, and hydrophobic effect. This integrated molecular modeling approach provided a comprehensive analysis of the structural effects of PBDEs on their binding properties with AhR at molecular level.