1 Department of Chemistry and Biochemistry, Concordia University, Montréal, Québec H4B 1R6, Canada. samira.norouzi@mail.concordia.ca.
2 Centre for Research in Molecular Modelling (CERMM), Concordia University, Montréal, Québec H4B 1R6, Canada.
3 Science College, Concordia University, Montréal, Québec H4B 1R6, Canada.
4 Department of Chemistry, John Abbott College, 21 275 Lakeshore Road, Sainte-Anne-de-Bellevue, Québec H9X 3L9, Canada.
5 School of Environmental Science and Engineering, Shanghai Jiao Tong University, Shanghai 200240, China.
6 Department of Chemical and Materials Engineering and Department of Physics, Concordia University, Montréal, Québec H4B 1R6, Canada.
7 School of Environmental Sciences, University of Guelph, Guelph, Ontario N1G 2W1, Canada.

To explore the use of molecular docking as a high throughput in silico screening tool for identifying chemicals of environmental health concern, we conducted a case study to assess endocrine disruption effects due to targeting of nuclear receptors (NRs) by chemicals with backbone structures like bisphenols, but with varied functional groups. The molecular docking analysis elucidates how functional groups of the chemicals, such as NH₂, Cl, and OCH₃, influence their interaction with the human estrogen receptor alpha (hERa), a key player in endocrine regulation. Through comparative docking analysis, we examined how bisphenol analogs interact with three distinct conformations of hERa: the apo structure and two structures with bound agonist and antagonist ligands. Water molecules within the protein and surrounding the ligand binding domain (LBD) were found to have little impact on the affinity of compounds binding to the receptor across various conformations. This can be attributed to the hydrophobic nature of the ligand-binding pocket, which consists mainly of hydrophobic amino acid residues and binding sites. In the assessment of bisphenol analogs compared to well established endocrine disrupting chemicals (EDCs), it was observed that these analogs exhibit characteristics commonly associated with endocrine disruptors. While compounds like BPA and BPF exhibited partial agonist activity, stimulating hERa activity to varying degrees, other compounds displayed non-agonist behavior, suggesting a different mode of interaction with the receptor. Further analysis revealed the significance of specific functional groups, such as hydroxyl or amine groups, on the aromatic ring of these compounds in modulating their binding affinity to hERa. Within the ligand binding site of hERa, amino acid residues Glu353, Arg394, and His524 have the capacity to form hydrogen bonds with hydroxyl or amine groups. Protonation or deprotonation of these groups can further alter their binding affinity, thereby influencing their interaction with estrogen receptors and subsequent estrogenic effects. Via this case study, we demonstrate the potential and provide best practices of using molecular docking as a new approach methodology (NAM) for chemical assessments and regulations.