Structure-activity relationship study on the binding of PBDEs with thyroxine transport proteins

Molecular docking and three-dimensional quantitative structure-activity relationships (3D-QSAR) were used to develop models to predict binding affinity of polybrominated diphenyl ether (PBDE) compounds to the human transthyretin (TTR). Based on the molecular conformations derived from the molecular docking, predictive comparative molecular similarity indices analysis (CoMSIA) models were developed. The results of CoMSIA models were as follows: leave-one-out (LOO) cross-validated squared coefficient q² (LOO) = 0.827 (full model, for all 28 compounds); q² (LOO) = 0.752 (split model, for 22 compounds in the training set); leave-many-out (LMO) cross-validated squared coefficient q² (LMO, two groups) = 0.723 ± 0.100 (full model, for all 28 compounds); q² (LMO, five groups) = 0.795 ± 0.030 (full model, for all 28 compounds); and the predictive squared correlation coefficient r²(pred)  = 0.928 (for six compounds in the test set). The developed CoMSIA models can be used to infer the activities of compounds with similar structural characteristics. In addition, the interaction mechanism between hydroxylated polybrominated diphenyl ethers (HO-PBDEs) and the TTR was explored. Hydrogen bonding with amino acid residues Asp74, Ala29, and Asn27 may be an important determinant for HO-PBDEs binding to TTR. Among them, forming hydrogen bonds with amino acid residues Asp74 might exert a more important function.