Docking and QSAR study on the binding interactions between polycyclic aromatic hydrocarbons and estrogen receptor

Synthesis, biological evaluation, theoretical calculations, QSAR and molecular docking studies of novel arylaminonaphthols as potent antioxidants and BChE inhibitors

A completely green protocol was developed for the synthesis of a series of arylaminonaphthol derivatives in the presence of N-ethylethanolamine (NEEA) as a catalyst under ultrasonic irradiation and solventless conditions. The major assets of this methodology were the use of non-toxic organic medium, available catalyst, mild reaction condition, and good to excellent yield of desired products. All of the synthesized products were screened for their in vitro antioxidant activity using DPPH, ABTS, and Ferric-phenanthroline assays and it was found that most of them are potent antioxidant agents. Also, their butyrylcholinesterase inhibitory activity has been investigated in vitro. All tested compounds exhibited potential inhibitory activity toward BuChE when compared to standard reference drug galantamine, however, compounds 4r, 4u, 4 g and 4x gave higher butyrylcholinesterase inhibitory with IC50 values of 14.78 ± 0.65 µM, 16.18 ± 0.50 µM, 20.00 ± 0.50 µM, and 20.28 ± 0.08 µM respectively. On the other hand, we employed density functional theory (DFT), calculations to analyze molecular geometry and global reactivity descriptors, and MESP analysis to predict electrophilic and nucleophilic attacks. A quantitative structure-activity relationship (QSAR) investigation was conducted on the antioxidant and butyrylcholinesterase properties of 25 arylaminonaphthol derivatives, resulting in robust and satisfactory models. To evaluate their anti-Alzheimer's activity, compounds 4 g, 4q, 4r, 4u, and 4x underwent docking simulations at the active site of the acetylcholinesterase (AChE) and butyrylcholinesterase (BChE), revealing why these compounds displayed superior activity, consistent with the biological findings.

Prepubertal exposure to polycyclic aromatic hydrocarbons are associated with early pubertal development onset in boys: A longitudinal study

OBJECTIVE

To investigate the effects of polycyclic aromatic hydrocarbons(PAHs) on puberty in boys.

METHODS

695 subjects were selected from four primary schools in Chongqing, China. 675 urine samples from these boys were collected four PAH metabolites: 1-hydroxypyrene, 2-hydroxynaphthoic, 2-hydroxyfluorene, and 9-hydroxyphenanthrene. Pubertal development of 695 boys was assessed at follow-up visits starting in December 2015 and occurring every six months thereafter until now, data used in this article ending in June 2021. A total of 12 follow-up visits were performed. Cox proportional hazards regression models were used to analyze the relationship between PAH metabolite concentrations and indicators of pubertal timing.

RESULTS

The mean age at puberty onset of testicular volume, facial hair, pubic hair, first ejaculation, and axillary hair in boys was 11.66, 12.43, 12.51, 12.72 and 13.70 years, respectively. Cox proportional hazards regression models showed that boys with moderate level of 1-OHPyr exposure was associated with earlier testicular development (hazard ratio [HR] = 1.276, 95% confidence interval [CI]: 1.006-1.619), with moderate level of 2-OHNap were at higher risk of early testicular development (HR = 1.273, 95% CI: 1.002-1.617) and early axillary hair development (HR = 1.355, 95% CI: 1.040-1.764), with moderate level of 2-OHFlu was associated with earlier pubic hair development (HR = 1.256, 95% CI: 1.001-1.577), with high level of 9-OHPhe were at higher risk of early fisrt ejaculation (HR = 1.333, 95% CI: 1.005-1.767) and early facial hair development (HR = 1.393, 95% CI: 1.059-1.831).

CONCLUSION

Prepubertal exposure to PAHs may be associated with earlier pubertal development in boys.

Tolerance mechanism of rice (Oryza sativa L.) seedings towards polycyclic aromatic hydrocarbons toxicity: The activation of SPX-mediated signal transduction to maintain P homeostasis

Plant tolerance to abiotic stress depends on fast molecular cascades involving stress perception, signal transduction, gene expression alterations, and metabolic rearrangement. This study sheds light on the tolerance mechanism of rice (Oryza sativa L.) towards the toxicity of the polycyclic aromatic hydrocarbons (PAHs), including phenanthrene (Phe), pyrene (Pyr), and benzo[a]pyrene (BaP). Results showed that three PAHs significantly activated the phosphoinositide signaling system involving the phosphorus (P) metabolism and homeostasis in rice roots. This activation increased phytic acid (IP6) levels to over 54.12% of the control (p < 0.05). Molecular docking verified that three PAHs occupied the IP6 binding site in SPX3, a negative regulatory factor of P homeostasis, where ARG229 interacted with PAHs via the van der Waals force. Moreover, the expression of gene encoding SPX3 was significantly downregulated 2.81-, 2.83-, and 2.18-fold under Phe, Pyr, and BaP stress, respectively, relative to the control. Conversely, the expression levels of the gene coding SDEL2 was significantly increased, promoting the degradation of SPX3. Ultimately, P absorption and nucleic acid synthesis were enhanced, alleviating the inhibition effect of PAHs on rice growth. Notably, Pyr demonstrated the strongest binding affinity for SPX3, confirming its critical interference with P homeostasis. These findings provide insight into the molecular mechanisms regulating plant responses to PAHs, and offer guidance for improving crop resistance against organic pollutants and protecting food security.

Association of the urinary polycyclic aromatic hydrocarbons with sex hormones stratified by menopausal status older than 20 years: a mixture analysis

We examined the relationships between exposure to polycyclic aromatic hydrocarbons (PAH) metabolites and sex hormones in pre- and postmenopausal women from the 2013-2016 National Health and Nutrition Examination Survey. The study comprised 648 premenopausal and 370 postmenopausal women (aged 20 years or older) with comprehensive data on PAH metabolites and sex steroid hormones. To evaluate the correlations between individual or mixture of the PAH metabolites and sex hormones stratified by menopausal status, we used linear regression and Bayesian kernel machine regression (BKMR). After controlling for confounders, 1-Hydroxynaphthalene (1-NAP) was inversely associated with total testosterone (TT), and 1-NAP, 3-Hydroxyfluorene (3-FLU), and 2-Hydroxyfluorene (2-FLU) were inversely associated with estradiol (E2). 3-FLU was positively associated with sex hormone-binding globulin (SHBG) and TT/E2, whereas 1-NAP and 2-FLU were inversely associated with free androgen index (FAI). In the BKMR analyses, chemical combination concentrations at or above the 55th percentile were inversely connected to E2, TT, and FAI values but positively correlated with SHBG when compared with the matching 50th percentile. In addition, we only found that mixed exposure to PAHs was positively associated with TT and SHBG in premenopausal women. Exposure to PAH metabolites, either alone or as a mixture, was negatively associated with E2, TT, FAI, and TT/E2 but positively associated with SHBG. These associations were stronger among postmenopausal women.

Antiandrogenic and Estrogenic Activity Evaluation of Oxygenated and Nitrated Polycyclic Aromatic Hydrocarbons Using Chemically Activated Luciferase Expression Assays

To establish the risk of the endocrine disrupting activity of polycyclic aromatic compounds, especially oxygenated and nitrated polycyclic aromatic hydrocarbons (oxy-PAHs and nitro-PAHs, respectively), antiandrogenic and estrogenic activities were determined using chemically activated luciferase expression (CALUX) assays with human osteoblast sarcoma cells. A total of 27 compounds including 9 oxy-PAHs (polycyclic aromatic ketones and quinones) and 8 nitro-PAHs was studied. The oxy-PAHs of 7H-benz[de]anthracen-7-one (BAO), 11H-benzo[a]fluoren-11-one (B[a]FO), 11H-benzo[b]fluoren-11-one (B[b]FO), and phenanthrenequinone (PhQ) exhibited significantly the potent inhibition of AR activation. All nitro-PAHs exhibited high antiandrogenic activities (especially high for 3-nitrofluoranthene (3-NFA) and 3-nitro-7H-benz[de]anthracen-7-one (3-NBAO)), and the AR inhibition was confirmed as noncompetitive for 3-NFA, 3-NBAO, and 1,3-dinitropyrene (1,3-DNPy). Antiandrogenic activity of 3-NFA demonstrated characteristically a U-shaped dose-response curve; however, the absence of fluorescence effect on the activity was confirmed. The prominent estrogenic activity dependent on dose-response curve was confirmed for 2 oxy-PAHs (i.e., B[a]FO and B[b]FO). Elucidating the role of AR and ER on the effects of polycyclic aromatic compounds (e.g., oxy- and nitro-PAHs) to endocrine dysfunctions in mammals and aquatic organisms remains a challenge.

Integration of in silico methods to determine endocrine-disrupting tobacco pollutants binding potency with steroidogenic genes: comprehensive QSAR modeling and ensemble docking strategies

A myriad of tobacco-associated chemicals may have possibilities to developmental/reproductive axis and endocrine-disruption impacts. Mostly they breach the biotransformation of cholesterol in mitochondria by interfering with steroidogenic pathway genes, prompting to adverse effects in steroid biosynthesis; however, studies are scanty. The quantitative structure-activity relationship (QSAR) modeling and comparative docking strategies were used to understand structural features of dataset compounds that influence developmental/reproductive toxicity and estrogen and androgen receptor-binding abilities, and to predict binding levels of toxicants with steroidogenic acute regulatory protein (StAR) and cholesterol side-chain cleavage enzyme (CYP11A1) active sites. Developed QSAR models presented good robustness and predictive ability that were determined from the applicability domain and, clustering and classification of chemicals by performing self-organizing maps. Accordingly, the exorbitant amount of polycyclic aromatic hydrocarbons (PAHs) and a limited number of other chemicals including N-nitrosamines and nicotine was represented as potential developmental/reproductive toxicants as well as estrogen and androgen receptor binders. From the docking analysis, hydrogen bonding, nonpolar, atomic π-stacking, and π-cation interactions were found between PAHs (bay and fjord structural pockets) and functional hotspot residues of StAR and CYP11A1, which strengthened the subtle structural changes at domains. These govern barrier effects to cholesterol binding and/or locking cholesterol to complicate its ejection from the Ω1 loop of StAR, and further mitigates steroid biosynthesis through cholesterol by CYP11A1; therefore, they are presumably considered as block-cluster mechanisms. These outcomes are significant to be hopeful to estimate developmental/reproductive toxicity and endocrine-disruption activities of other environmental pollutants, and could be useful for further assessment to discover binding mechanisms of PAHs with other steroidogenesis pathway genes.

The role of estrogen receptor β in fine particulate matter (PM2.5) organic extract-induced pulmonary inflammation in female and male mice

Fine particulate matter organic extract (Po) was reported to promote inflammation in the lung. Sex differences were reported in many inflammatory diseases. In this study, we investigated the effects of Po exposure on pulmonary inflammatory response and evaluated the role of sex in this process. While mice were exposed to 100 µg/m³ Po for 12 weeks by an inhalation exposure system, the lung histopathological analysis shown obvious inflammation, the cell numbers in bronchoalveolar lavage fluid (BALF) were significantly increased, and most inflammatory cytokines in BALF were upregulated. The results of factorial analysis of variance shown that there was an interaction between sex and Po exposure in the inflammatory cell numbers and the levels of tumor necrosis factor-α (TNF-α), interleukin-5 (IL-5), and growth-related oncogene/keratinocyte chemoattractant (GRO/KC). Notably, these changes and interactions were diminished while Po-exposed mice were administered with the estrogen receptor β (ERβ) antagonist. We speculated that sex might affect the levels of inflammatory indicators in BALF of Po-exposed mice and female mice were more prone to inflammation while exposed to Po. Moreover, ERβ was involved in these processes. To our knowledge, this is the first investigation about the role of sex in Po-induced adverse effects.

Critical target identification and human health risk ranking of metal ions based on mechanism-driven modeling

Huge amounts of metals have been released into environment due to various anthropogenic activities, such as smelting and processing of metals and subsequent application in construction, automobiles, batteries, optoelectronic devices, and so on, resulting in widespread detection in environmental media. However, some metal ions are considered as "Environmental health hazards", leading to serious human health concerns through affecting critical targets. Hence, it is necessary to quickly and effectively recognize the key target of metal ions in living organisms. Fortunately, the development of high-throughput analysis and in silico approaches offer a promising tool for target identification. In this study, the key oncogenic target (tumor suppressor protein, p53) was screened by network analysis based on the comparative toxicogenomics database (CTD). Some metal ions could bind to p53 core domain, impair its function and induce the development of cancer risk, but its mechanisms were still unclear. Therefore, a quantitative structure-activity relationship (QSAR) model was constructed to characterize the binding constants (K_a) between DNA binding domain of p53 (p53 DBD) and nine metal ions (Mg²⁺, Ca²⁺, Cu²⁺, Zn²⁺, Co²⁺, Ni²⁺, Mn²⁺, Fe³⁺ and Ba²⁺). It had good robustness and predictive ability, which could be used to predict the K_a values of other six metal ions (Li⁺, Ag⁺, Cs⁺, Cd²⁺, Hg²⁺ and Pb²⁺) within application domain. The results showed strong binding affinity between Cd²⁺/Hg²⁺/Pb²⁺ and p53 DBD. Subsequent mechanism analyses revealed that first hydrolysis constant (|logK_OH|) and polarization force (Z²/r) were key metal ion-characteristic parameters. The metal ions with weak hydrolysis constants and strong polarization forces could readily interact with N-containing histidine and S-containing cysteine of p53 DBD, which resulted in high K_a values. This study identified p53 as potential target for metal ions, revealed the key characteristics affecting the actions and provide a basic understanding of metal ions-p53 DBD interaction.

Aromatic hydrocarbon receptors in mitochondrial biogenesis and function

Mitochondria are ubiquitous membrane-bound organelles that not only play a key role in maintaining cellular energy homeostasis and metabolism but also in signaling and apoptosis. Aryl hydrocarbons receptors (AhRs) are ligand-activated transcription factors that recognize a wide variety of xenobiotics, including polyaromatic hydrocarbons and dioxins, and activate diverse detoxification pathways. These receptors are also activated by natural dietary compounds and endogenous metabolites. In addition, AhRs can modulate the expression of a diverse array of genes related to mitochondrial biogenesis and function. The aim of the present review is to analyze scientific data available on the AhR signaling pathway and its interaction with the intracellular signaling pathways involved in mitochondrial functions, especially those related to cell cycle progression and apoptosis. Various evidence have reported the crosstalk between the AhR signaling pathway and the nuclear factor κB (NF-κB), tyrosine kinase receptor signaling and mitogen-activated protein kinases (MAPKs). The AhR signaling pathway seems to promote cell cycle progression in the absence of exogenous ligands, whereas the presence of exogenous ligands induces cell cycle arrest. However, its effects on apoptosis are controversial since activation or overexpression of AhR has been observed to induce or inhibit apoptosis depending on the cell type. Regarding the mitochondria, although activation by endogenous ligands is related to mitochondrial dysfunction, the effects of endogenous ligands are not well understood but point towards antiapoptotic effects and inducers of mitochondrial biogenesis.

Assessing the chemical-induced estrogenicity using in silico and in vitro methods

Multiple substances are considered endocrine disrupting chemicals (EDCs). However, there is a significant gap in the early prioritization of EDC's effects. In this work, in silico and in vitro methods were used to model estrogenicity. Two Quantitative Structure-Activity Relationship (QSAR) models based on Logistic Regression and REPTree algorithms were built using a large and diverse database of estrogen receptor (ESR) agonism. A 10-fold external validation demonstrated their robustness and predictive capacity. Mechanistic interpretations of the molecular descriptors (C-026, nArOH,PW5, B06[Br-Br]) used for modelling suggested that the heteroatomic fragments, aromatic hydroxyls, and bromines, and the relative bond accessibility areas of molecules, are structural determinants in estrogenicity. As validation of the QSARs, ESR transactivity of thirteen persistent organic pollutants (POPs) and suspected EDCs was tested in vitro using the MMV-Luc cell line. A good correspondence between predictions and experimental bioassays demonstrated the value of the QSARs for prioritization of ESR agonist compounds.

New insights into the mechanism of hepatocyte apoptosis induced by typical organophosphate ester: An integrated in vitro and in silico approach

Apoptosis is one of the typical features of liver diseases, therefore molecular targets of hepatic apoptosis and regulatory mechanisms need to be further investigated. The caspases play important functions in the execution of apoptosis and many studies have focused on classical caspase-dependent cell death pathways. However, other types of cell death pathways (such as mitochondrial poly (ADP-ribose) polymerase-1 (PARP1) pathway) are suggested to be also as important as the caspase-mediated pathways in reflection of early toxic effects in hepatocytes, which requires additional research. In this work, an approach integrated in silico and in vitro was used to investigate the underlying toxicological mechanisms of hepatocyte apoptosis through the PARP1 dependent cell death pathway induced by triphenyl phosphate (TPP). Docking view showed that TPP could interact with helix αJ to affect the activation of PARP1 as a molecular initial event. In vitro assays suggested some biochemical events downstream of PARP1 activation, such as mitochondrial injury, apoptosis inducing factor (AIF) release, reactive oxygen species (ROS) production, and DNA damage. Moreover, the apoptosis was alleviated when cells were pretreated with PJ34 hydrochloride (PARP1 inhibitor), suggesting the mitochondrial PARP1 dependent pathway played a pivotal role in L02 cells apoptosis. This study indicated that PARP1 was an important molecular target in this process. And it also helped to understand the mechanism of hepatocytes apoptosis, early hepatic toxicity, and even liver diseases.

Urinary polycyclic aromatic hydrocarbons in relation to anthropometric measures and pubertal development in a cohort of Northern California girls

Supplemental Digital Content is available in the text.

Polycyclic aromatic hydrocarbons (PAHs) are a class of ubiquitous, environmental chemicals that may have endocrine disrupting capabilities. We investigated whether childhood exposure to PAHs was associated with adiposity and pubertal timing in a longitudinal study of 404 girls enrolled in the Northern California site of the Breast Cancer and the Environment Research Program cohort.

The critical factors affecting typical organophosphate flame retardants to mimetic biomembrane: An integrated in vitro and in silico study

Organophosphate flame retardants (OPFRs) have been reported to induce cytotoxicity in a structure-dependent manner. The toxic effects may be due to the damage of biomembrane integrity and/or the interference of membrane signal pathway. In this study, the damages of fifteen typical OPFRs (chlorinated phosphates, alkyl phosphates, aryl phosphates, and alkoxy phosphates) to mimetic biomembrane were determined by the electrochemical impedance spectroscopy (EIS). The molecular structure descriptors that characterized the action mechanisms were screened by stepwise regression. The six molecular descriptors (MATS7e, DLS_05, Mor19m, Mor22v, Mor12v and MATS8m) were screened to study the actions between OPFRs and mimetic biomembrane. A quantitative structure-activity relationship (QSAR) model was developed by the partial least squares (PLS) method. Statistical results indicated that the QSAR model had good robustness and mechanism interpretability. The distribution of atomic electronegativities (MATS7e) and atomic masses in three dimensional spaces (Mor19m) were the key factors influencing the actions between OPFRs and simulated biofilms. The compounds with strong electron-withdrawing property could invade the inner layer of membrane and destroy its integrity. High levels of steric hindrance could impair the damage capacity caused by electronegativity. Moreover, drug-like index (DLS_05), spatial structures of particle (Mor22v, Mor12v) and atomic masses (MATS8m) also affected the actions. The results revealed the mechanism of the actions of OPFRs with simulated biofilms and elucidated the key structural characteristics affecting the actions of OPFRs, which could provide theoretical basis for ecological risk assessment of OPFRs.

Urinary biomarkers of polycyclic aromatic hydrocarbons in pre- and peri-pubertal girls in Northern California: Predictors of exposure and temporal variability

BACKGROUND

Polycyclic aromatic hydrocarbons (PAHs), a class of chemicals produced as combustion by-products, have been associated with endocrine disruption. To understand exposure in children, who have been less studied than adults, we examined PAH metabolite concentrations by demographic characteristics, potential sources of exposure, and variability over time, in a cohort study of pre- and peri-pubertal girls in Northern California.

METHODS

Urinary concentrations of ten PAH metabolites and cotinine were quantified in 431 girls age 6-8 years at baseline. Characteristics obtained from parental interview, physical exam, and linked traffic data were examined as predictors of PAH metabolite concentrations using multivariable linear regression. A subset of girls (n = 100) had repeat measures of PAH metabolites in the second and fourth years of the study. We calculated the intraclass correlation coefficient (ICC), Spearman correlation coefficients, and how well the quartile ranking by a single measurement represented the four-year average PAH biomarker concentration.

RESULTS

Eight PAH metabolites were detected in ≥ 95% of the girls. The most consistent predictors of PAH biomarker concentrations were cotinine concentration, grilled food consumption, and region of residence, with some variation by demographics and season. After adjustment, select PAH metabolite concentrations were higher for Hispanic and Asian girls, and lower among black girls; 2-naphthol concentrations were higher in girls from lower income households. Other than 1-naphthol, there was modest reproducibility over time (ICCs between 0.18 and 0.49) and the concentration from a single spot sample was able to reliably rank exposure into quartiles consistent with the multi-year average.

CONCLUSIONS

These results confirm diet and environmental tobacco smoke exposure as the main sources of PAHs. Controlling for these sources, differences in concentrations still existed by race for specific PAH metabolites and by income for 2-naphthol. The modest temporal variability implies adequate exposure assignment using concentrations from a single sample to define a multi-year exposure timeframe for epidemiologic exposure-response studies.

Rigorous 3-dimensional spectral data activity relationship approach modeling strategy for ToxCast estrogen receptor data classification, validation, and feature extraction

The estrogenic potential (expressed as a score composite of 18 high throughput screening bioassays) of 1528 compounds from the ToxCast database was modeled by a 3-dimensional spectral data activity relationship approach (3D-SDAR). Due to a lack of ¹⁷ O nuclear magnetic resonance (NMR) simulation software, the most informative carbon-carbon 3D-SDAR fingerprints were augmented with indicator variables representing oxygen atoms from carbonyl and carboxamide, ester, sulfonyl, nitro, aliphatic hydroxyl, and phenolic hydroxyl groups. To evaluate the true predictive performance of the authors' model the United States Environmental Protection Agency provided them with a blind test set consisting of 2008 compounds. Of these, 543 had available literature data-their binding affinity served to estimate the external classification accuracy of the developed model: predictive accuracy of 0.62, sensitivity of 0.71, and specificity of 0.53 were obtained. Compared with alternative modeling techniques, the authors' model displayed very little reduction in performance between the modeling and the prediction set. A 3D-SDAR mapping technique allowed identification of structural features essential for estrogenicity: 1) the presence of a phenolic OH group or cyclohexenone, 2) a second aromatic or phenolic ring at a distance of 6 Å to 8 Å from the oxygen of the first phenol ring, 3) the presence of a methyl group approximately 6 Å away from the centroid of a phenol ring, and 4) a carbonyl group in close proximity (∼4 Å measured to the centroid) to 1 of the phenol rings. Environ Toxicol Chem 2017;36:823-830. Published 2016 Wiley Periodicals Inc. on behalf of SETAC. This article is a US government work and, as such, is in the public domain in the United States of America.

Environmental PAH exposure and male idiopathic infertility: a review on early life exposures and adult diagnosis

The male reproductive system is acutely and uniquely sensitive to a variety of toxicities, including those induced by environmental pollutants throughout the lifespan. Early life hormonal and morphological development results in several especially sensitive critical windows of toxicity risk associated with lifelong decreased reproductive health and fitness. Male factor infertility can account for over 40% of infertility in couples seeking treatment, and 44% of infertile men are diagnosed with idiopathic male infertility. Human environmental exposures are poorly understood due to limited available data. The latency between maternal and in utero exposure and a diagnosis in adulthood complicates the correlation between environmental exposures and infertility. The results from this review include recommendations for more and region specific monitoring of polycyclic aromatic hydrocarbon (PAH) exposure, longitudinal and clinical cohort considerations of exposure normalization, gene-environment interactions, in utero exposure studies, and controlled mechanistic animal experiments. Additionally, it is recommended that detailed semen analysis and male fertility data be included as endpoints in environmental exposure cohort studies due to the sensitivity of the male reproductive system to environmental pollutants, including PAHs.

In silico prediction and in vivo analysis of antiestrogenic potential of 2-isopropylthioxanthone (2-ITX) in juvenile goldfish (Carassius auratus)

Previous studies have shown both anti-estrogenic and anti-androgenic activities of 2-isopropylthioxanthone (2-ITX), a well known food contaminant, in in vitro assays. However, no data are available on the anti-estrogenic potentials and risks of 2-ITX in aquatic organisms. This work evaluated the potential endocrine disrupting effects of 2-ITX at the level of estrogen receptor (ER) signaling cascade using juvenile goldfish (Carassius auratus) as model. Firstly, we investigated the ligand binding efficiency of 2-ITX to the ligand binding domains (LBD) of goldfish ER subtypes using a molecular docking approach. Secondly, we assessed the effects of 2-ITX on E2-induced hepatic expression of ERα1, ERβ1, ERβ2, and vitellogenin (VTG) in vivo. Crosstalk between ER-VTG and aryl hydrocarbon receptor 2 (AhR2)-cytochrome P4501A (CYP1A) was also investigated. Fish were injected with increasing doses of 2-ITX ranging from 2 to 10µg/g BW, and results were compared to the effect of tamoxifen, a well-known ER modulator. We observed that compared to ERβ, the interaction potentials of 2-ITX to goldfish ERα1 LBD was more stable in the inactive receptor conformation. The in silico docking simulation analysis also revealed that 2-ITX acted as agonist for the goldfish AhR2 LBDs suggesting the ability of this compound to activate the cross-talk between the ERα- and AhR-signaling pathways. In vivo experiments confirm in silico simulation predictions demonstrating that 2-ITX reduced the estrogenicity of E2 at both transcriptional and post-transcriptional levels, indicating a clear anti-estrogenic effect. Co-exposure of E2 and 2-ITX also resulted in a significant decrease of CYP1A gene expression with respect to 2-ITX alone. Results from these studies collectively revealed that the antiestrogenic property of 2-ITX can be ascribed to a combination of effects on multiple signaling pathways suggesting the potential for this environmental contaminant to affect the hormonal control of reproductive processes in fish.

Use of computational modeling approaches in studying the binding interactions of compounds with human estrogen receptors

Estrogens have a whole host of physiological functions in many human organs and systems, including the reproductive, cardiovascular, and central nervous systems. Many naturally-occurring compounds with estrogenic or antiestrogenic activity are present in our environment and food sources. Synthetic estrogens and antiestrogens are also important therapeutic agents. At the molecular level, estrogen receptors (ERs) mediate most of the well-known actions of estrogens. Given recent advances in computational modeling tools, it is now highly practical to use these tools to study the interaction of human ERs with various types of ligands. There are two common categories of modeling techniques: one is the quantitative structure activity relationship (QSAR) analysis, which uses the structural information of the interacting ligands to predict the binding site properties of a macromolecule, and the other one is molecular docking-based computational analysis, which uses the 3-dimensional structural information of both the ligands and the receptor to predict the binding interaction. In this review, we discuss recent results that employed these and other related computational modeling approaches to characterize the binding interaction of various estrogens and antiestrogens with the human ERs. These examples clearly demonstrate that the computational modeling approaches, when used in combination with other experimental methods, are powerful tools that can precisely predict the binding interaction of various estrogenic ligands and their derivatives with the human ERs.

The interaction between organic phosphate ester and p53: an integrated experimental and in silico approach

Concerns have been raised in regards to the environmental impact of the more used organophosphate flame retardants (OPFRs). In this study, to better understand the relationship between molecular structural features of OPFRs and binding affinity for the tumor suppressor p53, an integrated experimental and in silico approach was used. From docking analysis, hydrogen bonding and hydrophobic interactions were found to be the dominant interactions, which implied the binding affinities of the compounds. The binding constants of 5 OPFRs were determined by surface plasmon resonance technology (SPR). Based on the observed interactions, appropriate molecular structural parameters were adopted to develop a quantitative structure-activity relationship (QSAR) model. The developed QSAR model had good robustness, predictive ability and mechanism interpretability. The interactions between the OPFRs and p53 (Ebinding) and the partition ability of the OPFRs into the bio-phase are main factors governing the binding affinities.

Endocrine disruption of the epigenome: a breast cancer link

The heritable component of breast cancer accounts for only a small proportion of total incidences. Environmental and lifestyle factors are therefore considered to among the major influencing components increasing breast cancer risk. Endocrine-disrupting chemicals (EDCs) are ubiquitous in the environment. The estrogenic property of EDCs has thus shown many associations between ongoing exposures and the development of endocrine-related diseases, including breast cancer. The environment consists of a heterogenous population of EDCs and despite many identified modes of action, including that of altering the epigenome, drawing definitive correlations regarding breast cancer has been a point of much discussion. In this review, we describe in detail well-characterized EDCs and their actions in the environment, their ability to disrupt mammary gland formation in animal and human experimental models and their associations with exposure and breast cancer risk. We also highlight the susceptibility of early-life exposure to each EDC to mediate epigenetic alterations, and where possible describe how these epigenome changes influence breast cancer risk.

Molecular docking and QSAR studies: noncovalent interaction between acephate analogous and the receptor site of human acetylcholinesterase

Twelve new compounds of acephate (Ace) analogues were synthesized and characterized by (31)P, (13)C, and (1)H NMR and IR spectroscopy. The probable insecticide potential of these compounds as well as 23 previously prepared molecules with a general skeleton of RC(O)-NH-P(O)X1X2 was predicted by PASS software. Docking analysis showed that hydrophobic interaction and hydrogen bonding were created between the functional groups of Ace derivatives and the receptor sites of acetylcholinesterase. PCA-QSAR indicated that the electronic descriptors are dominated in comparison with the structural descriptors. The experimental-QSAR (R(2) = 0.903 and VIF < 2.997) and DFT-QSAR (R(2) = 0.990 and VIF ≤ 10) models clarified that the net charge of functional groups contributes an important function in an inhibition mechanism. Validity and integrity of this model were confirmed by the LOO cross-validation method with q(2) = 0.940 and low residuals between the training and testing sets. The correlation matrix of DFT-QSAR model confirmed the molecular docking results.

Phthalates efficiently bind to human peroxisome proliferator activated receptor and retinoid X receptor α, β, γ subtypes: an in silico approach

This exhaustive in silico study looks into the molecular interactions of phthalates and their metabolites with human peroxisome proliferator-activated receptor (hPPAR) and retinoid X receptor (hRXR) α, β and γ subtypes--the nuclear receptor proteins function as transcription factors by regulating the expression of downstream genes. Apart from the much discussed plasticizer bisphenol A, we examined the binding affinities of 15 common diphthalates and their monophthalates, natural (linoleic acid, conjugated linoleic acid) and synthetic (bezafibrate, pioglitazone, GW 50156) ligands with hPPARs. In addition to these phthalates, specific natural (retinoic and phytanic acids) and synthetic (bexarotene, rosiglitazone) ligands were examined with hRXRs. The Maestro, Schrödinger Suite 2012 was used for the molecular docking study. In general, natural ligands of hPPAR showed less binding efficiencies than phthalic acid esters and drugs. The diphthalate di-iso-decyl phthalate showed the highest G score (-9.99) with hPPAR (γ), while its monophthalate (mono-iso-decyl phthalate) showed a comparatively less G score (-9.56). Though the PPAR modulator GW 50156 showed strong affinity with all hPPAR subtypes, its highest G score (-12.43) was with hPPARβ. Hazardous di(2-ethylhexyl)phthalate generally showed a greater preference to hRXRs than hPPARs, but its highest G score (-10.87) was with hRXRα; while its monophthalate (Mono(2-ethylhexyl)phthalate) showed a lesser G score (-8.59). The drug bexarotene showed the highest G score (-13.32) with hRXRβ. Moreover, bisphenol A showed more affinity towards hRXR. Briefly, this study gives an overview on the preference of phthalic acid esters, natural and synthetic ligands on to hPPAR and hRXR subtypes, which would lead to further in vitro mechanistic as well as in vivo preclinical and clinical studies.

Assessment of hydroxylated metabolites of polychlorinated biphenyls as potential xenoestrogens: a QSAR comparative analysis∗

Alternative methods, including quantitative structure-activity relationships (QSAR), are being used increasingly when appropriate data for toxicity evaluation of chemicals are not available. Approximately 40 mono-hydroxylated polychlorinated biphenyls (OH-PCBs) have been identified in humans. They represent a health and environmental concern because some of them have been shown to have agonist or antagonist interactions with human hormone receptors. This could lead to modulation of steroid hormone receptor pathways and endocrine system disruption. We performed QSAR analyses using available estrogenic activity (human estrogen receptor ER alpha) data for 71 OH-PCBs. The modelling was performed using multiple molecular descriptors including electronic, molecular, constitutional, topological, and geometrical endpoints. Multiple linear regressions and recursive partitioning were used to best fit descriptors. The results show that the position of the hydroxyl substitution, polarizability, and meta adjacent un-substituted carbon pairs at the phenolic ring contribute towards greater estrogenic activity for these chemicals. These comparative QSAR models may be used for predictive toxicity, and identification of health consequences of PCB metabolites that lack empirical data. Such information will help prioritize such molecules for additional testing, guide future basic laboratory research studies, and help the health/risk assessment community understand the complex nature of chemical mixtures.