Endocrine-disrupting chemicals (EDCs) are notable for their persistence, bioaccumulation, and associations with cancer. Human nuclear receptors (hNRs) are primary targets disrupted by these persistent EDCs, resulting in alterations to xenobiotic metabolism, lipid homeostasis, and endocrine function, which can lead to carcinogenic effects. Despite their hazardous effects, comprehensive studies on EDC interactions and their impacts on hNRs remain limited. Here, we profiled the interactions of persistent EDCs, including PFAS, plastic additives, bisphenols, polybrominated diphenyl ethers, and phthalates, with key hNRs such as PXR, CAR, PPARα, PPARγ, PPARδ, AR, and RORγt. Through controlled molecular docking simulations, we observed strong binding of the EDCs to these receptors. Further analysis showed that EDCs exhibit strong binding activity towards hNRs by preferentially interacting with hydrophobic amino acids, namely leucine, isoleucine, methionine, and phenylalanine. PFAS demonstrated the highest binding affinity, characterized by a combination of complementary hydrophobic interactions from their fluorinated carbon chains and polar interactions from their functional groups (e.g., carboxylate, sulfonate) across all receptors. Distinct polycyclic and hydrophobic trends, contributing to strong NR binding, were evident in non-PFAS and nonplastic EDCs. The hNR activity assay in HepG2 cells revealed agonistic effects of dicyclohexyl phthalate (DCHP) and di-2-ethylhexyl phthalate (DEHP) on most receptors, except for PPARα. The hNR transcription factor pathway assay in HepG2 cells demonstrated increased gene expression of VDRE and PXR, suggesting potential chronic effects on xenobiotic metabolism and calcium homeostasis. Overall, our findings demonstrate the need for further research into the endocrine disruption and carcinogenic effects of these persistent EDCs.

Per- and polyfluoroalkyl substances (PFAS) are environmentally persistent organofluorines linked to cancer, organ dysfunction, and other health problems. This study used quantitative structure–property relationship (QSPR) and quantitative structure–activity relationship (QSAR) modeling to examine the binding of PFAS to estrogen receptor alpha (ERα) and beta (ERβ). Molecular docking of 14,591 PFAS compounds was performed, and docking scores were used as a measure of receptor affinity. QSPR models were built for two datasets: the ERα and ERβ top binders (TBs), and a set of commonly exposed (CE) PFAS. These models quantified how chemical descriptors influence binding affinity. Across the models, higher density and electrophilicity indicated positive correlations with affinity, while surface tension indicated negative correlations. Electrostatic descriptors, including HOMO energy and positive Fukui index (F⁺ max), were part of the models but showed inconsistent trends. The CE QSPR models displayed correlations that conflicted with those of the TB models. Following QSPR analysis, 66 QSAR models were developed using a mix of top binders and experimental data. These models achieved strong performance, with R² values averaging 0.95 for training sets and 0.78 for test sets, that indicated reliable predictive ability. To improve generalizability, large-set QSAR models were created for each receptor. After outlier removal, these models reached R² values of 0.68–0.71, which supports their use in screening structurally diverse PFAS. Overall, QSPR and QSAR analyses reveal key chemical features that influence PFAS–ER binding. This predictive approach provides a scalable framework to assess the binding interactions of structurally diverse PFAS to ERs and other nuclear receptors. All the codes, data, and the GUI visualization of the results are freely available at sivaGU/QSPR-QSAR-Molecular-Visualization-Tool.