Molecular docking, molecular dynamics simulation, and structure-based 3D-QSAR studies on the aryl hydrocarbon receptor agonistic activity of hydroxylated polychlorinated biphenyls

Abstract

The binding interactions between hydroxylated polychlorinated biphenyls (HO-PCBs) and the aryl hydrocarbon receptor (AhR) are suspected of causing toxic effects. To understand the binding mode between HO-PCBs and AhR, and to explore the structural characteristics that influence the AhR agonistic activities of HO-PCBs, the combination of molecular docking, three-dimensional quantitative structure–activity relationship (3D-QSAR), and molecular dynamics (MD) simulations was performed. Using molecular docking, the HO-PCBs were docked into the binding pocket of AhR, which was generated by homology modeling. Comparative molecular similarity index analysis (CoMSIA) models were subsequently developed from three different alignment rules. The optimum 3D-QSAR model showed good predictive ability (q2=0.583, R2=0.913) and good mechanism interpretability. The statistical reliability of the CoMSIA model was also validated. In addition, molecular docking and MD simulations were applied to explore the binding modes between the ligands and AhR. The results obtained from this study may lead to a better understanding of the interaction mechanism between HO-PCBs and AhR.