Computational insights on agonist and antagonist mechanisms of estrogen receptor α induced by bisphenol A analogues

Structural analogues of bisphenol A (BPA) have become widely used as alternatives in BPA-free products. Most toxicological investigations have focused on the estrogenic activities of these analogues, which have been considered as potential environmental estrogens. However, recent studies revealed that certain BPA analogues could dramatically inhibit the proliferation of breast cancer cells, and exhibited strong anti-estrogenic effects compared with the antagonist 4-hydroxytamoxifen (OHT). Thus, we adopted computational models combining molecular dynamics simulations and binding free energy calculations to explore the underlying molecular basis of BPA analogues binding to estrogen receptor α (ERα). We also evaluated ligand-induced structural rearrangements of ERα at the atomic level. Conformational analyses showed that induced-fit H-bonding recognition by Thr347 was an important factor distinguishing antagonist from agonist BPA analogues. Moreover, antagonists of BPA analogues could indirectly induce the structural reposition of key helix 12 and produce an antagonistic conformation of ERα. Compared with OHT, the binding affinity of BPA analogues is stronger for antagonists than agonists. Taken together, we therefore propose computational indicators for screening of anti-estrogenic activities of BPA analogues, which may be beneficial for predicting the estrogenic or anti-estrogenic effects of BPA alternatives.