Comprehensive assessment of estrogenic activities of parabens by in silico approach and in vitro assays.

Parabens are ubiquitous pollutants in the environment and humans due to their wide applications in food, pharmaceuticals, and personal care products. Although the estrogenic activity of some parabens has been confirmed, the underlying mechanisms and the structure-estrogenic activity relationship are still largely unclear. Here, we systematically used in silico and in vitro approaches to investigate the estrogenic potency of typical parabens, including methyl-, ethyl-, propyl-, iso-propyl-, butyl-, iso-butyl- and benzyl-paraben. Molecular dynamics simulations and binding free energy calculations were combined to investigate the atomic-level mechanism of paraben binding to estrogen receptors (ERs). Computational analysis showed that ER were the targets of tested parabens and kept a stable agonist conformation. The calculated total binding free energies suggested that van der Waals interactions were the major driving forces for paraben-ER interaction and correlated with the structure of paraben side chains. In in vitro assays, paraben with an aromatic side chain, benzyl-paraben, showed the strongest estrogenic activity at 0.01 μM and the EC at 0.796 ± 0.307 μM, on par with levels commonly detected in human organs. Among tested parabens with an alkyl side chain, the estrogenicity increased as the side chain length increased from 1 to 4, but no significant difference appeared between parabens with isomeric alkyl side chains (propyl- vs isopropyl and butyl- vs iso-butylparaben). The estrogenic activity of parabens was significantly related to the calculated binding energies (R = 0.94, p = 0.0012), depending on the side chains of parabens. Our findings provide a significant mechanism for parabens to disrupt estrogenic function and considerations for structural optimization from the perspective of environmental protection.