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

Estrogenic disruption effects and formation mechanisms of transformation products during photolysis of preservative parabens

The ubiquitous presence of emerging contaminants (ECs) in the environment and their associated adverse effects has raised concerns about their potential risks. The increased toxicity observed during the environmental transformation of ECs is often linked to the formation of their transformation products (TPs). However, comprehension of their formation mechanisms and contribution to the increased toxicity remains an unresolved challenge. To address this gap, by combining quantum chemical and molecular simulations with photochemical experiments in water, this study investigated the formation of TPs and their molecular interactions related to estrogenic effect using the photochemical degradation of benzylparaben (BZP) preservative as a representative example. A non-targeted analysis was carried out and three previously unknown TPs were identified during the transformation of BZP. Noteworthy, two of these novel TPs, namely oligomers BZP-o-phenol and BZP-m-phenol, exhibited higher estrogenic activities compared to the parent BZP. Their IC50 values of 0.26 and 0.50 μM, respectively, were found to be lower than that of the parent BZP (6.42 μM). The binding free energies (ΔGbind) of BZP-o-phenol and BZP-m-phenol (-29.71 to -23.28 kcal·mol-1) were lower than that of the parent BZP (-20.86 kcal·mol-1), confirming their stronger binding affinities toward the estrogen receptor (ER) α-ligand binding domain. Subsequent analysis unveiled that these hydrophobic residues contributed most favorably to ER binding, with van der Waals interactions playing a significant role. In-depth examination of the formation mechanisms indicated that these toxic TPs primarily originated from the successive cleavage of ester bonds (OCH2C6H5 and COO group), followed by their combination with BZP*. This study provides valuable insight into the mechanisms underlying the formation of toxic TPs and their binding interactions causing the endocrine-disrupting effects. It offers a crucial framework for elucidating the toxicological patterns of ECs with similar structures.

Screening of the Antagonistic Activity of Potential Bisphenol A Alternatives toward the Androgen Receptor Using Machine Learning and Molecular Dynamics Simulation

Over the past few decades, extensive research has indicated that exposure to bisphenol A (BPA) increases the health risks in humans. Toxicological studies have demonstrated that BPA can bind to the androgen receptor (AR), resulting in endocrine-disrupting effects. In recent investigations, many alternatives to BPA have been detected in various environmental media as major pollutants. However, related experimental evaluations of BPA alternatives have not been systematically implemented for the assessment of chemical safety and the effects of structural characteristics on the antagonistic activity of the AR. To promote the green development of BPA alternatives, high-throughput toxicological screening is fundamental for prioritizing chemical tests. Therefore, we proposed a hybrid deep learning architecture that combines molecular descriptors and molecular graphs to predict AR antagonistic activity. Compared to previous models, this hybrid architecture can extract substantial chemical information from various molecular representations to improve the model's generalization ability for BPA alternatives. Our predictions suggest that lignin-derivable bisguaiacols, as alternatives to BPA, are likely to be nonantagonist for AR compared to bisphenol analogues. Additionally, molecular dynamics (MD) simulations identified the dihydrotestosterone-bound pocket, rather than the surface, as the major binding site of bisphenol analogues. The conformational changes of key helix H12 from an agonistic to an antagonistic conformation can be evaluated qualitatively by accelerated MD simulations to explain the underlying mechanism. Overall, our computational study is helpful for toxicological screening of BPA alternatives and the design of environmentally friendly BPA alternatives.

Antagonistic mechanisms of bisphenol analogues on the estrogen receptor α in zebrafish embryos: Experimental and computational studies

Bisphenol A (BPA) can disturb the estrogen receptor α (ERα)-mediated signaling pathway, which results in endocrine-disrupting effects and reproductive toxicity. Most BPA analogues as alternatives were evidenced to generate estrogenic activity as agonists or partial agonists of ERα. Recent studies indicated that certain BPA analogues, such as bisphenol M (BPM), bisphenol P (BPP), and bisphenol FL (BPFL), exhibited strong anti-estrogenic effects comparable with the typical antagonist 4-hydroxytamoxifen. However, conflicting findings were also observed for the compounds in different in vitro assays, and whether these BPA analogues can elicit an in vivo effect on ERα at environmentally relevant concentrations remains unknown. The underlying structural basis of estrogenic/anti-estrogenic activity should be further elucidated at the atomic level. To address these issues, we combined zebrafish-based in vivo and in silico methods to assess the effects of the compounds on ERα. The results show that the expressions of ERα-mediated downstream related genes in zebrafish embryos decreased after exposed to the compounds. Further molecular dynamics simulations were used to probe the antagonistic mechanisms of the compounds on ERα. The key H-bonding interactions were identified as important ligand recognition by ERα in the analysis of binding modes and binding free energy calculations. In summary, the current study provides preliminary in vivo evidence of fish species for the anti-estrogenic activity of certain BPA analogues.

Comparison of the mechanisms of estrogen disrupting effects between triphenyl phosphate (TPhP) and tris(1,3-dichloro-2-propyl) phosphate (TDCIPP)

As the typical aryl-organophosphate flame retardants (OPFRs), triphenyl phosphate (TPhP) and tris(1,3-dichloro-2-propyl) phosphate (TDCIPP) were reported to be estrogen disruptors. However, estrogen receptor α (ERα) binding experiments could not explain their biological effects. In this study, their action on ERα, G protein-coupled estrogen receptor (GPER) and the synthesis of 17β-estradiol (E2) were investigated using in vitro assays and molecular docking. The results showed that TPhP acted as an ERα agonist and recruited steroid receptor co-activator 1 (SRC1) and 3 (SRC3), which was found for the first time. Unlike TPhP, TDCIPP acted as an ERα antagonist. However, both TPhP and TDCIPP activated the estrogen pathway by GPER in SKBR3 cells which were lack of ERα. Although molecular docking results revealed that both TPhP and TDCIPP could dock into ERα and GPER, their substituent groups and combination mode might affect the receptor activation. In addition, by using estrogen biosynthesis assay in H295R cells, both of TPhP and TDCIPP were found to promote E2 synthesis and E2/T ratio involving their different alteration on levels of progesterone, testosterone and estrone, and expression of various key genes. Our data proposed estrogen-disrupting mechanism frameworks of TPhP and TDCIPP. Moreover, our results will contribute to future construction of adverse outcome pathway (AOP) framework of endocrine disruptors.

Structural basis for molecular recognition of G protein-coupled estrogen receptor by selected bisphenols

Complicated ligand-dependent signaling pathways of bisphenol A (BPA) and its analogues involve not only intranuclear estrogen receptor but also membrane receptor G protein-coupled estrogen receptor (GPER). However, the structural basis for molecular recognition of GPER by the environmental chemicals remains unknown. To reveal the structural dependence of GPER recognition by bisphenols, a systematic molecular dynamics simulation study was performed for selected bisphenols with different electron hybrid orbitals and substituents on their C atoms connecting two phenol rings. BPA was used as a control, bisphenol C(BPC) as an example for a connecting C with sp² hybrid orbitals to provide more ligand rigidity, bisphenol E(BPE) and bisphenol F(BPF) for decreased steric hindrance and hydrophobicity around the connecting C, and bisphenol B(BPB) and bisphenol AF(BPAF) for increased hydrophobicity and steric hindrance. All the tested bisphenols can bind with GPER at its classic orthosteric site to obtain GPER-ligand complexes, while van der Waals interactions and direct inter-molecular electrostatic energies provide the driving forces for ligand binding. Bulky substituents and structural rigidity of the connecting C dramatically impair hydrogen bonding between GPER and the bisphenols, which results in decreased contribution of both favorable intermolecular hydrogen bonds and unfavorable polar solvation effect to complex stability of BPB and BPC since decreased number of key residues is expected. Increase in substituent lipophilicity enhances the van der Waals interactions and favorable non-polar solvation effect. The six bisphenols of high structural similarity shared two key recognition residues, Leu137^TM3 and Trp272^TM6, the latter of which was in the highly conserved CWxP motif of TM6 and has been reported as key residue for G protein-coupled receptor activation. Based on the obtained knowledge, GPER affinity and relevant toxicity of BPA alternatives can be easily predicted, and the calculated binding free energies are consistent with the available experimental observations.

Initial Stages of Spontaneous Binding of Folate-Based Vectors to Folate Receptor-α Observed by Unbiased Molecular Dynamics

Active targeting is a prospective strategy for controlled drug delivery to malignant tumor tissues. One of the approaches relies on recognition of a bioactive ligand by a receptor expressed abundantly on the surface of cancer cell membranes. A promising ligand-receptor pair is folic acid (or its dianionic form, folate) combined with the folate receptor-α (FRα). A number of targeting drug delivery systems based on folate have been suggested, but the mechanism of binding of the ligand or its derivatives to the receptor is not fully known at the molecular level. The current study summarizes the results from unbiased all-atom molecular dynamics simulations at physiological conditions describing the binding of two forms of folate and four of its synthetically available derivatives to FRα. The models (ca. 185,000 atoms) contain one receptor molecule, embedded in the outer leaflet of a lipid bilayer, and one ligand, all immersed in saline. The bilayer represents a human cancer cell membrane and consists of 370 asymmetrically distributed lipid molecules from 35 types. The ability of the vector molecules to bind to the receptor, the position of binding, and the interactions between them are analyzed. Spontaneous binding on the nanosecond scale is observed for all molecules, but its time, position, and persistence depend strongly on the ligand. Only folate, 5-methyltetrahydrofolate, and raltitrexed bind selectively at the active site of the receptor. Two binding poses are observed, one of them (realized by raltitrexed) corresponding qualitatively to that reported for the crystallographic structure of the complex folate-FRα. Pemetrexed adsorbs nonspecifically on the protein surface, while methotrexate and pteroyl ornithine couple much less to the receptor. The molecular simulations reproduce qualitatively correctly the relative binding affinity measured experimentally for five of the ligands. Analysis of the interactions between the ligands and FRα shows that in order to accomplish specific binding to the active site, a combination of hydrogen bonding, π-stacking, and van der Waals and Coulomb attraction should be feasible simultaneously for the vector molecule. The reported results demonstrate that it is possible to observe receptor-ligand binding without applying bias by representing the local environment as close as possible and contain important molecular-level guidelines for the design of folate-based systems for targeted delivery of anticancer drugs.

Identification of active and inactive agonists/antagonists of estrogen receptor based on Tox21 10K compound library: Binomial analysis and structure alert

Endocrine disrupting chemicals can mimic, block, or interfere with hormones in organisms and subsequently affect their development and reproduction, which has raised significant public concern over the past several decades. To investigate (quantitative) structure-activity relationship, 8280 compounds were compiled from the Tox21 10K compound library. The results show that 50% activity concentrations of agonists are poorly related to that of antagonists because many compounds have considerably different activity concentrations between the agonists and antagonists. Analysis on the chemical classes based on mode of action (MOA) reveals that estrogen receptor (ER) is not the main target site in the acute toxicity to aquatic organisms. Binomial analysis of active and inactive ER agonists/antagonists reveals that ER activity of compounds is dominated by octanol/water partition coefficient and excess molar refraction. The binomial equation developed from the two descriptors can classify well active and inactive ER chemicals with an overall prediction accuracy of 73%. The classification equation developed from the molecular descriptors indicates that estrogens react with the receptor through hydrophobic and π-n electron interactions. At the same time, molecular ionization, polarity, and hydrogen bonding ability can also affect the chemical ER activity. A decision tree developed from chemical structures and their applications reveals that many hormones, proton pump inhibitors, PAHs, progestin, insecticides, fungicides, steroid and chemotherapy medications are active ER agonists/antagonists. On the other hand, many monocyclic/nonaromatic chain compounds and herbicides are inactive ER compounds. The decision tree and binomial equation developed here are valuable tools to predict active and inactive ER compounds.