Docking and QSAR study on the binding interactions between polycyclic aromatic hydrocarbons and estrogen receptor

Abstract

Little is known about the estrogenic activities of polycyclic aromatic hydrocarbons (PAHs) and the underlying mechanisms on estrogenic activities are still unclear. Molecular docking and quantitative structure-activity relationship (QSAR) were used to understand the relationship between molecular structural features and estrogenic activity, and to predict the binding affinity of PAHs to estrogen receptor α (ERα). From molecular docking analysis, hydrogen bonding as well as hydrophobic and π interactions were found between PAHs and ERα. Based on the docking results, appropriate molecular structural parameters were adopted to develop a QSAR model. Five descriptors were included in the QSAR model, which indicated that the estrogenic activity was related to molecular size, van der Waals volumes, shape profiles, polarizabilities and electropological states were significant parameters explaining the estrogenicity. Comparatively, the developed QSAR model had good robustness, predictive ability and mechanistic interpretability. Moreover, the applicability domain of the model was described.