Structure-based design of estrogen receptor-beta selective ligands

Enhancement of Female Rat Fertility via Ethanolic Extract from Nigella sativa L. (Black Cumin) Seeds Assessed via HPLC-ESI-MS/MS and Molecular Docking

The characteristic chemical composition of Nigella seeds is directly linked to their beneficial properties. This study aimed to investigate the phytochemical composition of Nigella sativa seeds using a 100% ethanolic extract using HPLC-ESI-MS/MS. Additionally, it explored the potential biological effects of the extract on female rat reproduction. Follicle Stimulating Hormone (FSH), Luteinizing Hormone (LH), Estrogen (E2), and Progesterone (P4) hormone levels were also assessed, along with the morphological and histological effects of the extract on ovarian, oviductal, and uterine tissues. Molecular docking was performed to understand the extract's activity and its role in regulating female reproduction by assessing its binding affinity to hormonal receptors. Twenty metabolites, including alkaloids, saponins, terpenes, flavonoids, phenolic acids, and fatty acids, were found in the ethanolic extract of N. sativa seeds through the HPLC-ESI-MS/MS study. The N. sativa seed extract exhibited strong estrogenic and LH-like activities (p < 0.05) with weak FSH-like activity. Furthermore, it increased the serum levels of LH (p < 0.05), P4 hormones (p < 0.001), and E2 (p < 0.0001). Molecular docking results displayed a strong interaction with Erβ, LH, GnRH, and P4 receptors, respectively. Based on these findings, N. sativa seeds demonstrated hormone-like activities, suggesting their potential as a treatment for improving female fertility.

Evaluation of structural and thermodynamic insight of ERβ with DPN and derivatives through MMGBSA/MMPBSA methods

Estrogen receptors (ERs) are nuclear factors that exist as two subtypes: ERα and ERβ. Among the different selective ERβ agonist ligands, the widely used ERβ-selective agonist DPN (diarylpropionitrile) is highlighted. Recent experimental and thermodynamic information between R-DPN and S-DPN enantiomers with ERβ is important for evaluating further the ability of MD simulations combined with end-point methods to reproduce experimental binding affinity and generate structural insight not provided through crystallographic data. In this research, starting from crystallographic data and experimental binding affinities, we explored the structural and thermodynamic basis of the molecular recognition of ERβ with DPN and derivatives through triplicate MD simulations combined with end-point methods. Conformational analysis showed some regions with the highest mobility linked to ligand association that, at the time, impacted the total protein fluctuation. Binding free energy (ΔG) analysis revealed that the Molecular Mechanics Generalized-Born Surface Area (MMGBSA) approach was able to reproduce the experimental tendency with a strong correlation (R = 0.778), whereas per-residue decomposition analysis revealed that all the systems interacted strongly with eight residues (L298, E305, L339, M340, L343, F356, H475, and L476). The comparison between theoretical studies using the MMGBSA approach with experimental results provides new insights for drug designing of new DPN derivatives.

Developmental toxicity and estrogenicity of glyphosate in zebrafish in vivo and in silico studies

As the most heavily used herbicide globally, glyphosate (GLY) has been detected in a variety of environments and has raised concerns about its ecological and health effects. There is debate as to whether GLY may disrupt the endocrine system. Here, we investigated the developmental toxicity of GLY in zebrafish based on deep learning-enabled morphometric analysis (DLMA). In addition, the estrogenic activity of GLY was assessed by endocrine disruption prediction, docking study and in vivo experiments. Results showed that exposure to environmental concentrations of GLY negatively impacted zebrafish development, causing yolk edema and pericardial edema. Endocrine disruption prediction suggested that GLY may target estrogen receptors (ER). Molecular docking analysis revealed binding of GLY to three zebrafish ER. In vivo zebrafish experiment, GLY enhanced the protein levels of ERα and the mRNA levels of cyp19a, HSD17b1, vtg1, vtg2, esr1, esr2a and esr2b. These results suggest that GLY may act as an endocrine disruptor by targeting ER, which warrants further attention for its potential toxicity to aquatic animals.

The identification of cianidanol as a selective estrogen receptor beta agonist and evaluation of its neuroprotective effects on Parkinson's disease models

AIM

The present study aims to identify selective estrogen receptor beta (ERβ) agonists and to evaluate the neuroprotective mechanism in Parkinson's disease (PD) models.

MAIN METHODS

In-silico studies were carried out using Maestro and GROMACS. Neuroprotective activity and apoptosis were evaluated using cytotoxicity assay and flow cytometry respectively. Gene expression studies were carried out by reverse transcription polymerase chain reaction. Motor and cognitive functions were assessed by actophotometer, rotarod, catalepsy, and elevated plus maze. The neuronal population in the substantia nigra and striatum of rats was assessed by hematoxylin and eosin staining.

KEY FINDINGS

Cianidanol was identified as a selective ERβ agonist through virtual screening. The cianidanol-ERβ complex is stable during the 200 ns simulation and was able to retain the interactions with key amino acid residues. Cianidanol (25 μM) prevents neuronal toxicity and apoptosis induced by rotenone in differentiated SH-SY5Y cells. Additionally, cianidanol (25 μM) increases the expression of ERβ, cathepsin D, and Nrf2 transcripts. The neuroprotective effects of cianidanol (25 μM) were reversed in the presence of a selective ERβ antagonist. In this study, we found that selective activation of ERβ could decrease the transcription of α-synuclein gene. Additionally, cianidanol (10, 20, 30 mg/kg, oral) improves the motor and cognitive deficit in rats induced by rotenone.

SIGNIFICANCE

Cianidanol shows neuroprotective action in PD models and has the potential to serve as a novel therapeutic agent for the treatment of PD.

In silico guided in vitro study of traditionally used medicinal plants reveal the alleviation of post-menopausal symptoms through ERβ binding and MAO-A inhibition

The slumping level of estrogen and serotonin in menopausal women is directly associated with the occurrence of menopausal symptoms where, estrogen receptor-β (ERβ) and monoamine oxidase-A (MAO-A) are directly involved. The present investigation aimed for validation of promising plants traditionally used to alleviate menopausal symptoms with ERβ mediated MAO-A inhibition potential through in silico disease-target network construction using Cytoscape plugins followed by molecular docking of phytomolecules through AutoDock vina. ADMET parameters of identified bioactive phytomolecules were analysed through swissADME and ProTox II. The efficacy of promising plant leads was further established through in vitro ERβ competitive binding, MAO-A inhibition, enzyme kinetics and free radical quenching assays. In silico analysis suggested glabrene (ΔG = -9.7 Kcal/mol) as most promising against ERβ in comparison to 17β-estradiol (ΔG = -11.4 Kcal/mol) whereas liquiritigenin (ΔG = -9.4 Kcal/mol) showed potential binding with MAO-A in comparison to standard harmine (ΔG = -8.8 Kcal/mol). In vitro analysis of promising plants segregated Glycyrrhiza glabra (IC50 = 0.052 ± 0.007 μg/ml) as most promising, followed by Hypericum perforatum (IC50 = 0.084 ± 0.01 μg/ml), Trifolium pratense (IC50 = 0.514 ± 0.01 μg/ml) and Rumex nepalensis (IC50 = 2.568 ± 0.11 μg/ml). The enzyme kinetics of promising plant leads showed reversible and competitive nature of inhibition against MAO-A. The potency of plant extracts in quenching free radicals was at par with ascorbic acid. The identified four potent medicinal plants with ERβ selective, MAO-A inhibitory and free radical quenching abilities could be used against menopausal symptoms however, finding needs to be validated further for menopausal symptoms in in vivo conditions for drug development.Communicated by Ramaswamy H. Sarma.

In silico simulations and molecular descriptors to predict in vitro transactivation potencies of Baikal seal estrogen receptors by environmental contaminants

Baikal seals (Pusa sibirica) are vulnerable to high levels of organic pollutants. Here, we evaluated the transactivation potencies of bisphenols (BPs) and hydroxylated polychlorinated biphenyls (OH-PCBs) via the Baikal seal estrogen receptor α and β (bsERα and bsERβ) using in vitro and in silico approaches. In vitro reporter gene assays showed that most BPs and OH-PCBs exhibited estrogenic activity with bsER sub-type-specific potency. Among the BPs tested, bisphenol AF showed the lowest EC50 for both bsERs. 4'-OH-CB50 and 4'-OH-CB30 showed the lowest EC50 among OH-PCBs tested for bsERα and bsERβ, respectively. 4-((4-Isopropoxyphenyl)-sulfonyl)phenol, 4'-OH-CB72, and 4'-OH-CB121 showed weak bsERα-specific transactivation. Only 4-OH-CB107 did not affect both bsERs. In silico docking simulations revealed the binding affinities of these chemicals to bsERs and partially explained the in vitro results. Using the in silico simulations and molecular descriptors as explanatory variables and the in vitro results as objective variables, the quantitative structure-activity relationship (QSAR) models constructed for classification and regression accurately separated bsER-active compounds from non-active compounds and predicted the in vitro bsERα- and bsERβ-transactivation potencies, respectively. The QSAR models also suggested that chemical polarity, van der Waals surface area, bridging atom structure, position of the phenolic-OH group, and ligand interactions with key residues of the ligand binding pocket are critical variables to account for the bsER transactivation potency of the test compounds. We also succeeded in constructing computational models for predicting in vitro transactivation potencies of mouse ERs in the same manner, demonstrating the applicability of our approach independent of species-specific responses.

Molecular Docking Study of Isoxazole Indole Derivatives (B2A2 Series) as Promising Selective Estrogen Receptor Modulators & Anticancer Drugs

A series of 7 compounds with isoxazole - indole - γ-resorcylic acid scaffold, segregated into B2 & A2 series, wherein, B2 comprises Compounds: 13, 14, 15 & 16 and A2 comprises Compounds: 10, 11 & 12, on the basis of the variable substituents at the indole, resorcinol and isoxazole end of the scaffold as in Figure: 1, were designed and docked with human estrogen receptor: 1ERRα. The Binding affinity (BA) and the interacting amino acids compared with reference selective estrogen receptor modulators (SERM's) such as Raloxifene, Estradiol, Bazedoxifene, Bisphenol, Genistein, Daidzein, Ormiloxifene, Tamoxifen, 6-hydroxy-naphthalen-2yl-benzo(D)-isoxazol-6-ol(1) using PyRx software and their ADME properties predicted with SWISS ADME online tool. Significant similarities and minor differences in the binding pattern between the key interacting aminoacids such as Arg 394, Glu 353, Asp 351, Leu 346, Leu 525, Trp 383, Phe 404, Ala 350, Leu 387, Met 421 responsible for ER agonist/antagonist affinity found in the binding cavity of a 1 Errα -Bazedoxifene/1 Errα -raloxifene/1 Errα -estradiol docked complex AND 1 Errα -isoxazole-indole- resorcinol docked complex indicate their promising potential to serve as potent ER agonists in bone or ER antagonists against breast cancer and other cancer diseases. The Compounds with highest BA is of the order: BA (A1series)>B1series>/<BA(A2 series)>/=BA (B2 series) exceptions: compounds: 4, 5 of B1 series & compound:13 of B2 series with identical and least BA values.BA(6)=BA(8)>BA(7)>BA(2)>BA(9)=BA(1)>BA(12)>BA(10)=BA(15)=BA(11)=BA(3)>BA(14)=BA(16)>BA(4)=BA(5)=BA(13).

Anti-Candida, docking studies, and in vitro metabolism-mediated cytotoxicity evaluation of Eugenol derivatives

The high morbidity and mortality rates of Candida infections, especially among immunocompromised patients, are related to the increased resistance rate of these species and the limited therapeutic arsenal. In this context, we evaluated the anti-Candida potential and the cytotoxic profile of eugenol derivatives. Anti-Candida activity was evaluated on C. albicans and C. parapsilosis strains by minimum inhibitory concentration (MIC), scanning electron microscopy (SEM), and molecular docking calculations at the site of the enzyme lanosterol-14-α-demethylase active site, responsible for ergosterol formation. The cytotoxic profile was evaluated in HepG2 cells, in the presence and absence of the metabolizing system (S9 system). The results indicated compounds 1b and 1d as the most active ones. The compounds have anti-Candida activity against both strains with MIC ranging from 50 to 100 μg ml^-1 . SEM analyses of 1b and 1d indicated changes in the envelope architecture of both C. albicans and C. parapsilosis like the ones of eugenol and fluconazole, respectively. Docking results of the evaluated compounds indicated a similar binding pattern of fluconazole and posaconazole at the lanosterol-14-α-demethylase binding site. In the presence of the S9 system, compound 1b showed the same cytotoxicity profile as fluconazole (1.08 times) and compound 1d had 1.23 times increase in cytotoxicity. Eugenol and other evaluated compounds showed a significant increase in cytotoxicity. Our results suggest compound 1b as a promising starting point candidate to be used in the design of new anti-Candida agent prototypes.

Development of a Novel 18F-Labeled Probe for PET Imaging of Estrogen Receptor β

Estrogen receptor beta (ERβ) is an important ER subtype that plays crucial roles in many physiological and pathological disorders. Herein, we developed the probe [¹⁸F]PVBO for in vivo ERβ targeted PET imaging and obtained promising results. The nonradioactive PVBO showed a 12.5-fold stronger binding affinity to ERβ than to ERα in vitro. In vitro assays revealed the specific uptake of [¹⁸F]PVBO by DU145 cells. The uptake of [¹⁸F]PVBO by DU145 xenografts increased during the 120 min dynamic scanning, with a maximum uptake of 2.80 ± 0.30% ID/g. Based on time activity curves (TACs), the injection of [¹⁸F]PVBO with unlabeled PVBO or ERB-041 resulted in a significant signal reduction with the tumor/muscle (T/M) ratio <1 at 30, 60, 75, and 120 min post-injection (p < 0.05). [¹⁸F]PVBO demonstrates the feasibility of noninvasively imaging ERβ-positive tumors by small-animal PET and provides a new strategy for visualizing ERβ in vivo.

Evaluating the correlation of binding affinities between isothermal titration calorimetry and fragment molecular orbital method of estrogen receptor beta with diarylpropionitrile (DPN) or DPN derivatives

Estrogen receptors (ERs) are ligand-activated transcription factors, with two subtypes ERα and ERβ. The endogenous ligand of ERs is the common 17β-estradiol, and the ligand-binding pocket of ERα and ERβ is very similar. Nevertheless, some ERβ-selective agonist ligands have been reported. DPN (diarylpropionitrile) is a widely used ERβ-selective agonist; however, the structure of the ERβ-DPN complex has not been solved. Therefore, the bound-state conformation of DPN and its enantioselectivity remain unresolved. In this report, we present the structures of the complexes of ERβ with DPN or its derivatives that include a chlorine atom by the X-ray crystallography. Additionally, we measured the binding affinity between ERβ and DPN or derivatives by isothermal titration calorimetry (ITC) and estimated the binding affinity by fragment molecular orbital (FMO) calculations. We also examined the correlation between the ITC data and results from the FMO calculations. FMO calculations showed that S-DPN interacts strongly with three amino acids (Glu305, Phe356, and His475) of ERβ, and ITC measurements confirmed that the chlorine atom of the DPN derivatives enhances binding affinity. The enthalpy change by ITC correlated strongly with the interaction energy (total IFIEs; inter-fragment interaction energies) calculated by FMO (R = 0.870). We propose that FMO calculations are a valuable approach for enhancing enthalpy contributions in drug design, and its scope of applications includes halogen atoms such as chlorine. This study is the first quantitative comparison of thermodynamic parameters obtained from ITC measurements and FMO calculations, providing new insights for future precise drug design.

Insilico Docking Study of Isoxazole Indole Linked Resorcinol Derivatives as Promising Selective Estrogen Receptor Modulators & Anticancer Drugs

A series of 9 compounds with isoxazole-indole-γ-resorcylic acid scaffold, segregated into B1 & A1 series, wherein, B1 comprises compounds:1,3,4,5, & 9 and A1comprises compounds: 2,6,7, & 8 , on the basis of variable substituents at the indole , resorcinol and isoxazole end of the scaffold as in Fig. 1, were designed and docked with human estrogen receptor:1ERRα. The binding affinity (BA) and the interacting amino acids compared with reference selective estrogen receptor modulators (SERMs) such as Raloxifene, Estradiol, Bazedoxifene, Bisphenol, Genistein, Daidzein, Ormiloxifene,Tamoxifen,6-hydroxy-naphthalen-2yl-benzo(D)-isoxazol-6-ol(1)(WAY-397) using PyRx software and their ADME properties predicted with SWISS ADME online tool. Significant similarities and minor differences in the binding pattern between the key interacting aminoacids such as Arg 394,Glu 353, Asp 351, Leu 346, Leu 525, Trp 383,Phe 404 ,Ala 350, Leu 387, Met 421 responsible for ER agonist/antagonist activity found in the binding cavity of a 1 Errα -Bazedoxifene/1 Errα -raloxifene/1 Errα -estradiol docked complex AND 1 Errα -isoxazole-indole- resorcinol docked complex indicate their promising potential to serve as potent ER agonists in bone or ER antagonists against breast cancer and other cancer diseases. The Compounds with Highest BA is of the order: BA (A1series)>B1 series & BA(6)=BA(8)>BA(7)>BA(2)>BA(9)=BA(1)>BA(3)>BA(4)=BA(5).

MolHyb: A Web Server for Structure-Based Drug Design by Molecular Hybridization

Molecular hybridization is a widely used ligand design method in drug discovery. In this study, we present MolHyb, a web server for structure-based ligand design by molecular hybridization. The input of MolHyb is a protein file and a seed compound file. MolHyb tries to generate novel ligands through hybridizing the seed compound with helper compounds that bind to the same protein target or similar proteins. To facilitate the job of getting helper compounds, we compiled a modeled protein-ligand structure database as an extension to crystal structures in the PDB database by placing the bioactive compounds in ChEMBL into their corresponding 3D protein binding pocket properly. MolHyb works by searching for helper compounds from the protein-ligand structure database and migrating chemical moieties from helper compounds to the seed compound efficiently. Hybridization is performed at both cyclic and acyclic bonds. The users can also input their own helper compounds to MolHyb. We hope that MolHyb will be a useful tool for rational drug design. MolHyb is freely available at http://molhyb.xundrug.cn/.

Human Estrogen Receptor Alpha Antagonists, Part 3: 3-D Pharmacophore and 3-D QSAR Guided Brefeldin A Hit-to-Lead Optimization toward New Breast Cancer Suppressants

The estrogen receptor α (ERα) is an important biological target mediating 17β-estradiol driven breast cancer (BC) development. Aiming to develop innovative drugs against BC, either wild-type or mutated ligand-ERα complexes were used as source data to build structure-based 3-D pharmacophore and 3-D QSAR models, afterward used as tools for the virtual screening of National Cancer Institute datasets and hit-to-lead optimization. The procedure identified Brefeldin A (BFA) as hit, then structurally optimized toward twelve new derivatives whose anticancer activity was confirmed both in vitro and in vivo. Compounds as SERMs showed picomolar to low nanomolar potencies against ERα and were then investigated as antiproliferative agents against BC cell lines, as stimulators of p53 expression, as well as BC cell cycle arrest agents. Most active leads were finally profiled upon administration to female Wistar rats with pre-induced BC, after which 3DPQ-12, 3DPQ-3, 3DPQ-9, 3DPQ-4, 3DPQ-2, and 3DPQ-1 represent potential candidates for BC therapy.

Naturally occurring plant-based anticancerous candidates as prospective ABCG2 inhibitors: an in silico drug discovery study

Abstract

ATP-binding cassette transporter G2 (ABCG2) is an efflux transporter related to the clinical multidrug resistance (MDR) phenomenon. Identifying ABCG2 inhibitors could help discover extraordinary curative strategies for carcinoma remediation. Hitherto, there is no medication drug inhibiting ABCG2 transporter, notwithstanding that a considerable number of drugs have been submitted to clinical-trial and investigational phases. In the search for unprecedented chemical compounds that could inhibit the ABCG2 transporter, an in silico screening was conducted on the Naturally Occurring Plant-based Anticancer Compound-Activity-Target (NPACT) database containing 1574 compounds. Inhibitor-ABCG2 binding affinities were estimated based on molecular docking and molecular minimization (MM) calculations and compared to a co-crystallized inhibitor (BWQ) acting as a reference inhibitor. Molecular dynamics (MD) simulations pursued by molecular mechanics-generalized Born surface area (MM-GBSA) binding energy estimations were further executed for compounds with MM-GBSA//MM binding energies lower than BWQ (calc. − 60.5 kcal/mol). NPACT00968 and NPACT01545 demonstrated auspicious inhibitory activities according to binding affinities (ΔGbinding) over the 100 ns MD simulations that were nearly one and a half folds compared to BWQ (− 100.4, − 94.7, and − 62.9 kcal/mol, respectively). Throughout the 100 ns MD simulations, structural and energetical analyses unveiled outstanding stability of the ABCG2 transporter when bound with NPACT00968 and NPACT01545. In silico calculations hold a promise for those two inhibitors as drug candidates of ABCG2 transporter and emphasize that further in vitro and in vivo experiments are guaranteed.

Graphical abstract

Graphene oxide enhanced the endocrine disrupting effects of bisphenol A in adult male zebrafish: Integrated deep learning and metabolomics studies

In our previous studies, it was found that graphene oxide (GO) reduced the endocrine disruption of bisphenol A (BPA) in zebrafish embryo and larvae, but through different mechanisms. In this study, adult male zebrafish were selected to further understand the interactions between GO and BPA considering that adult zebrafish have different uptake pathways and metabolism from embryo and larvae. BPA was predicted to bind with the estrogen receptor α (ERα) with a probability of 98.1% by training a directed-message passing deep neural network model, and was confirmed by molecular docking analysis. The results were in accordance with the significantly increased vitellogenin (VTG) and estradiol (E2) levels, while decreased testosterone (T) and follicle-stimulating hormone (FSH) levels in the adult male zebrafish after 7 d exposure to 500 μg/L BPA. Compared to BPA single exposure group, the presence of GO led to significantly lower T and FSH levels and fewer spermatozoa, indicating that GO enhanced the endocrine disruption effects of BPA in the adult zebrafish. Metabolomics analysis revealed that 5 μg/L BPA could elicit changes in the metabolome, and the responses were correlated with BPA concentrations. Metabolic pathway analysis revealed more disturbance was caused by the mixture of GO and BPA compared to BPA alone, including three additional pathways and stronger perturbations on carbohydrate, lipid, and amino acid metabolism, fortifying that GO exaggerated the toxic effects of BPA. This was opposite to the depression effect observed in zebrafish embryo and larvae, magnifying that the joint effects of exposure to nanomaterials and endocrine disrupting chemicals are relevant to the life stages of organisms.

Human Estrogen Receptor α Antagonists. Part 1: 3-D QSAR-Driven Rational Design of Innovative Coumarin-Related Antiestrogens as Breast Cancer Suppressants through Structure-Based and Ligand-Based Studies

The estrogen receptor α (ERα) represents a 17β-estradiol-inducible transcriptional regulator that initiates the RNA polymerase II-dependent transcriptional machinery, pointed for breast cancer (BC) development via either genomic direct or genomic indirect (i.e., tethered) pathway. To develop innovative ligands, structure-based (SB) three-dimensional (3-D) quantitative structure-activity relationship (QSAR) studies have been undertaken from structural data taken from partial agonists, mixed agonists/antagonists (selective estrogen receptor modulators (SERMs)), and full antagonists (selective ERα downregulators (SERDs)) correlated with either wild-type or mutated ERα receptors. SB and ligand-based (LB) alignments allow us to rule out guidelines for the SB/LB alignment of untested compounds. 3-D QSAR models for ERα ligands, coupled with SB/LB alignment, were revealed to be useful tools to dissect the chemical determinants for ERα-based anticancer activity as well as to predict their potency. The herein developed protocol procedure was verified through the design and potency prediction of 12 new coumarin-based SERMs, namely, 3DQ-1a to 3DQ-1e, that upon synthesis turned to be potent ERα antagonists by means of either in vitro or in vivo assays (described in the second part of this study).

FragRep: A Web Server for Structure-Based Drug Design by Fragment Replacement

The design of efficient computational tools for structure-guided ligand design is essential for the drug discovery process. We hereby present FragRep, a new web server for structure-based ligand design by fragment replacement. The input is a protein and a ligand structure, either from protein data bank or from molecular docking. Users can choose specific substructures they want to modify. The server tries to find suitable fragments that not only meet the geometric requirements of the remaining part of the ligand but also fit well with local protein environments. FragRep is a powerful computational tool for the rapid generation of ligand design ideas; either in scaffold hopping or bioisosteric replacing. The FragRep Server is freely available to researchers and can be accessed at http://xundrug.cn/fragrep.

Discovery of a multi-target compound for estrogen receptor-positive (ER+) breast cancer: Involvement of aromatase and ERs

Despite intense research, breast cancer remains the leading cause of cancer-related death in women worldwide, being estrogen receptor-positive (ER⁺) the most common subtype. Nowadays, aromatase inhibitors (AIs), the selective estrogen receptor modulator (SERM) tamoxifen and the selective estrogen receptor down-regulator (SERD) fulvestrant are used as therapeutic options for ER⁺ breast cancer, since they interfere directly with the production of estrogens and with the activation of estrogen-dependent signaling pathways. Despite the success of these treatments, the occurrence of resistance limits their clinical efficacy, demanding the development of novel therapies. Recently, multi-target compounds emerged as promising therapeutic strategies for ER⁺ breast cancer, as they can potentially modulate several important targets simultaneously. In line with this, in this work, the anti-cancer properties and multi-target action of 1,1-Bis(4-hydroxyphenyl)-2-phenylbut-1-ene, tamoxifen bisphenol (1,1-BHPE), were evaluated in an ER⁺ breast cancer cell model (MCF-7aro cells). Molecular docking analysis predicted that 1,1-BHPE was able to bind to aromatase, ERα and ERβ. In vitro studies showed that, although it did not present anti-aromatase activity, 1,1-BHPE reduced aromatase protein levels and interfered with ERα and ERβ signaling pathways, acting as an ERα antagonist and inducing ERβ up-regulation. Through these mechanisms, 1,1-BHPE was able to impair breast cancer growth and induce apoptosis. This represents an important therapeutic advantage because the main players responsible for estrogen production and signaling are modulated by a single compound. To the best of our knowledge, this is the first study describing the anti-cancer properties of 1,1-BHPE as a multi-target compound specific for ER⁺ breast cancer.

Synthesis and evaluation of 17α-triazolyl and 9α-cyano derivatives of estradiol

A variety of 17α-triazolyl and 9α-cyano derivatives of estradiol were prepared and evaluated for binding to human ERβ in both a TR-FRET assay, as well as ERβ and ERα agonism in cell-based functional assays. 9α-Cyanoestradiol (5) was nearly equipotent as estradiol as an agonist for both ERβ and ERα. The potency of the 17α-triazolylestradiol analogs is considerably more variable and depends on the nature of the 4-substituent of the triazole ring. While rigid protein docking simulations exhibited significant steric clashing, induced fit docking providing more protein flexibility revealed that the triazole linker of analogs 2d and 2e extends outside of the traditional ligand binding domain with the benzene ring located in the loop connecting helix 11 to helix 12.

Indole Derivative Interacts with Estrogen Receptor Beta and Inhibits Human Ovarian Cancer Cell Growth

Ovarian cancer remains the leading cause of mortality among gynecological tumors. Estrogen receptor beta (ERβ) expression has been suggested to act as a tumor suppressor in epithelial ovarian cancer by reducing both tumor growth and metastasis. ERβ expression abnormalities represent a critical step in the development and progression of ovarian cancer: for these reasons, its re-expression by genetic engineering, as well as the use of targeted ERβ therapies, still constitute an important therapeutic approach. 3-{[2-chloro-1-(4-chlorobenzyl)-5-methoxy-6-methyl-1H-indol-3-yl]methylene}-5-hydroxy-6-methyl-1,3-dihydro-2H-indol-2-one, referred to here as compound 3, has been shown to have cytostatic as well cytotoxic effects on various hormone-dependent cancer cell lines. However, the mechanism of its anti-carcinogenic activity is not well understood. Here, we offer a possible explanation of such an effect in the human ovarian cancer cell line IGROV1. Chromatin binding protein assay and liquid chromatography mass spectrometry were exploited to localize and quantify compound 3 in cells. Molecular docking was used to prove compound 3 binding to ERβ. Mass spectrometry-based approaches were used to analyze histone post-translational modifications. Finally, gene expression analyses revealed a set of genes regulated by the ERβ/3 complex, namely CCND1, MYC, CDKN2A, and ESR2, providing possible molecular mechanisms that underline the observed antiproliferative effects.

How technical progress reshaped behavioral neuroendocrinology during the last 50 years… and some methodological remarks

The first issue of Hormones and Behavior was published 50 years ago in 1969, a time when most of the techniques we currently use in Behavioral Endocrinology were not available. Researchers have during the last 5 decades developed techniques that allow measuring hormones in small volumes of biological samples, identify the sites where steroids act in the brain to activate sexual behavior, characterize and quantify gene expression correlated with behavior expression, modify this expression in a specific manner, and manipulate the activity of selected neuronal populations by chemogenetic and optogenetic techniques. This technical progress has considerably transformed the field and has been very beneficial for our understanding of the endocrine controls of behavior in general, but it did also come with some caveats. The facilitation of scientific investigations came with some relaxation of methodological exigency. Some critical controls are no longer performed on a regular basis and complex techniques supplied as ready to use kits are implemented without precise knowledge of their limitations. We present here a selective review of the most important of these new techniques, their potential problems and how they changed our view of the hormonal control of behavior. Fortunately, the scientific endeavor is a self-correcting process. The problems have been identified and corrections have been proposed. The next decades will obviously be filled with exciting discoveries in behavioral neuroendocrinology.

Construction of benzofuranone library via a metal-free, one-pot intermolecular condensation, and their application as efficient estrogen receptor β modulators

Facile synthesis of benzofuranone was achieved through a metal-free, one-pot intermolecular condensation between α-hydroxy aryl ketones and resorcinol derivatives. A library of 20 compounds with moderate to good overall yields was prepared. These compounds showed strong binding toward estrogen receptors along with good selectivity for ERβ (>190-fold over ERα). Anti-proliferative activity on DU-145, U-87, and MCF-7 cells gave inhibition IC₅₀ values in the low μM range, which suggested the promising potential therapeutic applications of these new classes of benzofuranones.

Stronger estrogenic and antiandrogenic effects on zebrafish larvae displayed by 6:2 polyfluoroalkyl phosphate diester than the 8:2 congener at environmentally relevant concentrations

Polyfluoroalkyl phosphate esters (PAPs) are one kind of emerging polyfluoroalkyl substances in the environment. However, their in vivo toxicities are largely unknown, especially at environmental relevant concentrations. To fill this gap, zebrafish embryos were exposed to 6:2 or 8:2 diPAP at environmentally relevant concentrations (0.5, 5, 50 ng/L) for 7 d. 6:2 and 8:2 diPAPs upregulated the mRNA and protein levels of aromatase in the exposed larvae, and elevated estradiol (E2) and vitellogenin (VTG) levels, but reduced testosterone (T) and 11-ketotestosterone (11-KT) levels, demonstrating estrogenic and antiandrogenic effects. Among the three ER subtypes, ERβ2 displayed the highest in vivo mRNA expression and the lowest in silico binding energies, suggesting that it was the main target ER subtype responsible for the estrogenic effect. Molecular simulation results indicated that diPAPs and E2 could bind to one common residue, arginine (Arg) 87, in the binding pocket of ERβ2, inducing similar estrogenic disruption mechanisms as E2. Both compounds could form hydrophobic interaction with glutamic acid (Glu) 12 and tryptophan (Trp) 80 and two hydrogen bonds with Arg81 of androgen receptor (AR) ligand-binding domains (LBDs) in antagonistic mode, resulting in a reduced level of AR upon exposure. The in silico binding energies of 6:2 diPAP with both ER and AR were lower than 8:2 diPAP, explaining the observed greater in vivo estrogenic and antiandrogenic activities of 6:2 diPAP. This study provided the first line of evidences that diPAPs could display adverse effects on the endocrine functions of fish species.

LARMD: integration of bioinformatic resources to profile ligand-driven protein dynamics with a case on the activation of estrogen receptor

Protein dynamics is central to all biological processes, including signal transduction, cellular regulation and biological catalysis. Among them, in-depth exploration of ligand-driven protein dynamics contributes to an optimal understanding of protein function, which is particularly relevant to drug discovery. Hence, a wide range of computational tools have been designed to investigate the important dynamic information in proteins. However, performing and analyzing protein dynamics is still challenging due to the complicated operation steps, giving rise to great difficulty, especially for nonexperts. Moreover, there is a lack of web protocol to provide online facility to investigate and visualize ligand-driven protein dynamics. To this end, in this study, we integrated several bioinformatic tools to develop a protocol, named Ligand and Receptor Molecular Dynamics (LARMD, http://chemyang.ccnu.edu.cn/ccb/server/LARMD/ and http://agroda.gzu.edu.cn:9999/ccb/server/LARMD/), for profiling ligand-driven protein dynamics. To be specific, estrogen receptor (ER) was used as a case to reveal ERβ-selective mechanism, which plays a vital role in the treatment of inflammatory diseases and many types of cancers in clinical practice. Two different residues (Ile373/Met421 and Met336/Leu384) in the pocket of ERβ/ERα were the significant determinants for selectivity, especially Met336 of ERβ. The helix H8, helix H11 and H7-H8 loop influenced the migration of selective agonist (WAY-244). These computational results were consistent with the experimental results. Therefore, LARMD provides a user-friendly online protocol to study the dynamic property of protein and to design new ligand or site-directed mutagenesis.

Rational approaches of drug design for the development of selective estrogen receptor modulators (SERMs), implicated in breast cancer

Drug discovery and development have gained momentum due to the rational drug design by engaging computational tools and bioinformatics methodologies. Bioisosteric replacements and hybrid molecular approaches are the other inventive processes, used by medicinal chemists for the desired modifications of leads for clinical drug candidates. SERMs, ought to produce inhibitory activity in breast, uterus and agonist activity in other tissues, are beneficial for estrogen-like actions. ER subtypes α and β are hormone dependent modulators of intracellular signaling and gene expression, and development of ER selective ligands could be an effective approach for treatment of breast cancer. This report has critically investigated the possible designing considerations of SERMs, their in silico interactions, and potent pharmacophore generation approaches viz. indole, restricted benzothiophene [3, 2-b] indole, carborane, xanthendione, combretastatin A-4, organometallic heterocycles, OBHS-SAHA hybrids, benzopyranones, tetrahydroisoquinolines, Dig G derivatives and their specifications in drug design and development, to rationally improve the understanding in drug discovery. This also includes various strategies for the development of dual inhibitors for the management of antiestrogenic resistance.

How To Design Selective Ligands for Highly Conserved Binding Sites: A Case Study Using N-Myristoyltransferases as a Model System

A model system of two related enzymes with conserved binding sites, namely N-myristoyltransferase from two different organisms, was studied to decipher the driving forces that lead to selective inhibition in such cases. Using a combination of computational and experimental tools, two different selectivity-determining features were identified. For some ligands, a change in side-chain flexibility appears to be responsible for selective inhibition. Remarkably, this was observed for residues orienting their side chains away from the ligands. For other ligands, selectivity is caused by interfering with a water molecule that binds more strongly to the off-target than to the target. On the basis of this finding, a virtual screen for selective compounds was conducted, resulting in three hit compounds with the desired selectivity profile. This study delivers a guideline on how to assess selectivity-determining features in proteins with conserved binding sites and to translate this knowledge into the design of selective inhibitors.

Discovery of novel natural flavonoids as potent antiviral candidates against hepatitis C virus NS5B polymerase

The non-structural 5B (NS5B) polymerase of hepatitis C virus (HCV) is an attractive target for antiviral intervention. Quercetagetin (Que) is a natural flavonoid, which has been exhibited to have anti-HCV property through inhibition of RNA binding to NS5B. The last few decades have witnessed a growing interest in the extraction of natural flavonoids with a plethora of different biological activities. Considering the high therapeutic potential of Que, the aim of this study is to explore wide structure entities with potent activity using Que as a prototype. A virtual screen protocol involving docking and molecular dynamics has been performed to examine the potency of forty-three natural flavonoids which recently extracted from plants for inhibition of NS5B. During two screening stages, two compounds 24 and 41 were identified to have more favorable binding affinity to NS5B as compared to Que. The comparative analysis showed that there is a significant difference in the binding free energy of Que and 41 (ΔΔG_bind = -11.17 kcal/mol). It was revealed that van der Waals (vdW) interaction drives the binding process of both 24 and 41 and plays an important role in increasing their activities relative to Que. PHE162 serves as a crucial residue in both the NS5B-24 and NS5B-41 systems, contributing the most vdW energy by π-π interaction, suggesting that aromatic interactions are critical for the binding of 24 and 41 to NS5B. Moreover, hydrogen bond analysis indicates that the hydrogen bonds formed by LYS98, THR137, ASP164 and ARG168, can play important roles in the increased binding affinity of 41 to NS5B relative to Que. The findings of this study will provide useful structure-activity relationship (SAR) guidelines for the design of novel inhibitors with improved/enhanced therapeutic activities in the treatment of hepatitis C.

A synthesis strategy for tetracyclic terpenoids leads to agonists of ERβ

Natural product and natural product-like molecules continue to be important for the development of pharmaceutical agents, as molecules in this class play a vital role in the pipeline for new therapeutics. Among these, tetracyclic terpenoids are privileged, with >100 being FDA-approved drugs. Despite this significant pharmaceutical success, there remain considerable limitations to broad medicinal exploitation of the class due to lingering scientific challenges associated with compound availability. Here, we report a concise asymmetric route to forging natural and unnatural (enantiomeric) C19 and C20 tetracyclic terpenoid skeletons suitable to drive medicinal exploration. While efforts have been focused on establishing the chemical science, early investigations reveal that the emerging chemical technology can deliver compositions of matter that are potent and selective agonists of the estrogen receptor beta, and that are selectively cytotoxic in two different glioblastoma cell lines (U251 and U87).

Many natural-product like drugs have a tetracyclic terpenoid core. Here, the authors developed a synthesis of triterpene-like tetracyclic systems, and apply this method to the preparation of a number of enantiomeric compounds, two of which are very selective ligands for estrogen receptor beta

Identification of 10-dehydrooxyglycyuralin E as a selective human estrogen receptor alpha partial agonist

Selective estrogen receptor modulators (SERMs) act as either agonist or antagonist of estrogen receptor (ER) in a tissue selective manner and have been used in several diseases such as breast cancer, postmenopausal syndrome, osteoporosis, and cardiovascular diseases. However, current SERMs may also increase the risk of serious side effects and trigger drug resistance. Herein, a screening program, that was designed to search for novel SERMs, resulted in the identification of a series of 2-arylbenzofuran-containing compounds that are ligands for ERα, when applying the Gaussia-luciferase reporter assay. One of these compounds, 10-dehydrooxyglycyuralin E (T9) was chemically synthesized. T9 showed anti-estrogenic/proliferative activity in ERα-positive breast cancer cells. Pretreatment of T9 prevented the mRNA expression of GREB1, which is an estrogen response gene. Furthermore, by an in silico docking simulation study we demonstrated that T9 showed interactions directly to ERα. Taken together, these results demonstrated that T9 is a candidate of SERMs and a useful seed compound for the foundation of the selective activity of SERMs.

Receptor-Receptor Interactions as a Widespread Phenomenon: Novel Targets for Drug Development?

The discovery of receptor-receptor interactions (RRI) has expanded our understanding of the role that G protein-coupled receptors (GPCRs) play in intercellular communication. The finding that GPCRs can operate as receptor complexes, and not only as monomers, suggests that several different incoming signals could already be integrated at the plasma membrane level via direct allosteric interactions between the protomers that form the complex. Most research in this field has focused on neuronal populations and has led to the identification of a large number of RRI. However, RRI have been seen to occur not only in neurons but also in astrocytes and, outside the central nervous system, in cells of the cardiovascular and endocrine systems and in cancer cells. Furthermore, RRI involving the formation of macromolecular complexes are not limited to GPCRs, being also observed in other families of receptors. Thus, RRI appear as a widespread phenomenon and oligomerization as a common mechanism for receptor function and regulation. The discovery of these macromolecular assemblies may well have a major impact on pharmacology. Indeed, the formation of receptor complexes significantly broadens the spectrum of mechanisms available to receptors for recognition and signaling, which may be implemented through modulation of the binding sites of the adjacent protomers and of their signal transduction features. In this context, the possible appearance of novel allosteric sites in the receptor complex structure may be of particular relevance. Thus, the existence of RRI offers the possibility of new therapeutic approaches, and novel pharmacological strategies for disease treatment have already been proposed. Several challenges, however, remain. These include the accurate characterization of the role that the receptor complexes identified so far play in pathological conditions and the development of ligands specific to given receptor complexes, in order to efficiently exploit the pharmacological properties of these complexes.

Noteworthy effect of slight variation in aliphatic chain length of trisubstituted imidazole inhibitors against epidermal growth factor receptor L858R/T790M/C797S mutant in cancer therapy

11h is a very potent inhibitor against epidermal growth factor receptor triple mutation L858R/T790M/C797S (EGFR^TM ) with 13-fold stronger potency than the FDA-approved osimertinib. Recently, two new EGFR^TM inhibitors, 11d and 11e, were reported which revealed 2.8- and 2.3-fold stronger potency than 11h, respectively. 11h, 11d, and 11e have the same structures but differ only in their aliphatic chain length. However, the exact effects of differential aliphatic chain length on the inhibitory potencies of these compounds require further elaboration at the atomistic level, hence the objective of this report. Various computational tools were employed for this purpose. From our findings, it was revealed that van der Waals (vdW) interactions modulate the binding mechanisms of these inhibitors and play the most important role in the differential inhibitory activities of 11d, 11h, and 11e. The strong hydrogen bond formation between the aliphatic chain of 11d and key residue ARG841 was recognized as the reason for its higher activity and inhibitory potency relative to 11h and 11e. Moreover, the extension of the N-terminal loop into the active site for vdW interaction with the phenyl group of 11e and carbon-hydrogen bond formed between the aliphatic chain of 11e and LEU718 engendered a higher activity of 11e than 11h.

Insight Analysis of Promiscuous Estrogen Receptor α-Ligand Binding by a Novel Machine Learning Scheme

Estrogen receptor α (ERα) plays a significant role in occurrence of breast cancer and may cause various adverse side-effects when ERα is an off-target protein. A theoretical model was derived to predict the binding affinity of ERα using the pharmacophore ensemble/support vector machine (PhE/SVM) scheme to consider the promiscuous characteristic of ERα. The estimations by PhE/SVM were discovered to be in good agreement with the observed values for those training molecules ( n = 31, r² = 0.80, q_CV² = 0.77, RMSE = 0.57, s = 0.58), test molecules ( n = 179, q² = 0.91-0.96, RMSE = 0.33, s = 0.26) and outliers ( n = 15, q² = 0.80-0.86, RMSE = 0.56, s = 0.49). When subjected to various statistical validations, the PhE/SVM model consistently fulfilled the strictest criteria. A mock test also asserted its predictivity. When compared with crystal structures, the calculated results are consistent with the reported ERα-ligand co-complex structure, and the plasticity nature of ERα is also disclosed. Consequently, this precise, fast, and robust model can be adopted to predict ERα-ligand binding affinities and to design safer non-ERα-targeted pharmaceuticals in the process of drug discovery and development.

A-C Estrogens as Potent and Selective Estrogen Receptor-Beta Agonists (SERBAs) to Enhance Memory Consolidation under Low-Estrogen Conditions

Estrogen receptor-beta (ERβ) is a drug target for memory consolidation in postmenopausal women. Herein is reported a series of potent and selective ERβ agonists (SERBAs) with in vivo efficacy that are A-C estrogens, lacking the B and D estrogen rings. The most potent and selective A-C estrogen is selective for activating ER relative to seven other nuclear hormone receptors, with a surprising 750-fold selectivity for the β over α isoform and with EC₅₀s of 20-30 nM in cell-based and direct binding assays. Comparison of potency in different assays suggests that the ER isoform selectivity is related to the compound's ability to drive the productive conformational change needed to activate transcription. The compound also shows in vivo efficacy after microinfusion into the dorsal hippocampus and after intraperitoneal injection (0.5 mg/kg) or oral gavage (0.5 mg/kg). This simple yet novel A-C estrogen is selective, brain penetrant, and facilitates memory consolidation.

Identification by Molecular Docking ofHomoisoflavones from Leopoldia comosa as Ligands of Estrogen Receptors

The physiological responses to estrogen hormones are mediated within specific tissues by at least two distinct receptors, ERα and ERβ. Several natural and synthetic molecules show activity by interacting with these proteins. In particular, a number of vegetal compounds known as phytoestrogens shows estrogenic or anti-estrogenic activity. The majority of these compounds belongs to the isoflavones family and the most representative one, genistein, shows anti-proliferative effects on various hormone-sensitive cancer cells, including breast, ovarian and prostate cancer. In this work we describe the identification of structurally related homoisoflavones isolated from Leopoldia comosa (L.) Parl. (L. comosa), a perennial bulbous plant, potentially useful as hormonal substitutes or complements in cancer treatments. Two of these compounds have been selected as potential ligands of estrogen receptors (ERs) and the interaction with both isoforms of estrogen receptors have been investigated through molecular docking on their crystallographic structures. The results provide evidence of the binding of these compounds to the target receptors and their interactions with key residues of the active sites of the two proteins, and thus they could represent suitable leads for the development of novel tools for the dissection of ER signaling and the development of new pharmacological treatments in hormone-sensitive cancers.

Identification of d-seco modified steroid derivatives with affinity for estrogen receptor α and β isoforms using a non-transcriptional fluorescent cell assay in yeast

Synthesis and biological evaluation of steroidal derivatives with anticancer properties is an active area of drug discovery. Here we measured the relative affinities of d-seco modified steroidal derivatives for estrogen receptor α, estrogen receptor β or androgen receptor ligand binding domains using an optimized non-transcriptional fluorescent cell assay in yeast. Ligand binding domains of steroid receptors were expressed in-frame with yellow fluorescent protein in the yeast Saccharomyces cerevisiae. Addition of known steroid ligands to yeast expressing the appropriate cognate receptor results in increased fluorescence intensity, enabling estimation of receptor binding affinities in a dose-response and time-dependent manner. Relative binding affinities of d-seco modified steroidal derivatives 1-4 were then evaluated using this yeast system by live cell fluorimetry and fluorescence microscopy, coupled with in vitro cytotoxicity and in silico molecular docking studies. d-Seco estratriene derivative 2displayed strong affinity for both estrogen receptor α and β ligand binding domains and negligible affinity for the androgen receptor ligand binding domain. Compound 2 also showed moderate cytotoxicity against estrogen receptor positive MCF-7 breast adenocarcinoma cells. In addition to identification of new ligands for steroid receptors, this assay could also be used to filter out compounds with potential for off-target interactions with steroid receptors during the early stages of compound screening.

Synthesis and structure-activity relationships of 1-benzylindane derivatives as selective agonists for estrogen receptor beta

The estrogen receptor beta (ERβ) selective agonist is considered a promising candidate for the treatment of estrogen deficiency symptoms in ERβ-expressing tissues, without the risk of breast cancer, and multiple classes of compounds have been reported as ERβ selective agonists. Among them, 6-6 bicyclic ring-containing structures (e.g., isoflavone phytoestrogens) are regarded as one of the cyclized analogues of isobutestrol 5b, and suggest that other cyclized scaffolds comprising 5-6 bicyclic rings could also act as selective ERβ ligands. In this study, we evaluated the selective ERβ agonistic activity of 1-(4-hydroxybenzyl)indan-5-ol 7a and studied structure-activity relationship (SAR) of its derivatives. Some functional groups improved the properties of 7a; introduction of a nitrile group on the indane-1-position resulted in higher selectivity for ERβ (12a), and further substitution with a fluoro or a methyl group to the pendant phenyl ring was also preferable (12b, d, and e). Subsequent chiral resolution of 12a identified that R-12a has a superior profile over S-12a. This is comparable to diarylpropionitrile (DPN) 5c, one of the promising selective ERβ agonists and indicates that this indane-based scaffold has the potential to provide better ERβ agonistic probes.

Treatment with an orally bioavailable prodrug of 17β-estradiol alleviates hot flushes without hormonal effects in the periphery

Estrogen deprivation has a profound effect on the female brain. One of the most obvious examples of this condition is hot flushes. Although estrogens relieve these typical climacteric symptoms, many women do not want to take them owing to unwanted side-effects impacting, for example, the uterus, breast and blood. Therefore, there is a need for developing safer estrogen therapies. We show here that treatment with 10β,17β-dihydroxyestra-1,4-dien-3-one (DHED), a novel brain-targeting bioprecursor prodrug of the main human estrogen, 17β-estradiol, alleviates hot flushes in rat models of thermoregulatory dysfunction of the brain. Oral administration of DHED elicits a significant reduction of tail skin temperature (TST) rise representing hot flushes in the morphine-dependent ovariectomized rat model and results in the restoration of estrogen deprivation-induced loss of diurnal rhythm in TST. These beneficial effects occur without detrimental peripheral hormonal exposure; thus, the treatment avoids potentially harmful stimulation of estrogen-sensitive peripheral organs, including the uterus and the anterior pituitary, or the proliferation of MCF-7a breast cancer cell xenografts. Our promising preclinical assessments warrant further considerations of DHED for the development of a brain-selective 17β-estradiol therapy to relieve hot flushes without undesirable peripheral side-effects.

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.

Estrogen Receptor β Activation Rapidly Modulates Male Sexual Motivation through the Transactivation of Metabotropic Glutamate Receptor 1a

In addition to the transcriptional activity of their liganded nuclear receptors, estrogens, such as estradiol (E2), modulate cell functions, and consequently physiology and behavior, within minutes through membrane-initiated events. The membrane-associated receptors (mERs) underlying the acute effects of estrogens on behavior have mostly been documented in females where active estrogens are thought to be of ovarian origin. We determined here, by acute intracerebroventricular injections of specific agonists and antagonists, the type(s) of mERs that modulate rapid effects of brain-derived estrogens on sexual motivation in male Japanese quail. Brain aromatase blockade acutely inhibited sexual motivation. Diarylpropionitrile (DPN), an estrogen receptor β (ERβ)-specific agonist, and to a lesser extent 17α-estradiol, possibly acting through ER-X, prevented this effect. In contrast, drugs targeting ERα (PPT and MPP), GPR30 (G1 and G15), and the Gq-mER (STX) did not affect sexual motivation. The mGluR1a antagonist LY367385 significantly inhibited sexual motivation but mGluR2/3 and mGluR5 antagonists were ineffective. LY367385 also blocked the behavioral restoration induced by E2 or DPN, providing functional evidence that ERβ interacts with metabotropic glutamate receptor 1a (mGluR1a) signaling to acutely regulate male sexual motivation. Together these results show that ERβ plays a key role in sexual behavior regulation and the recently uncovered cooperation between mERs and mGluRs is functional in males where it mediates the acute effects of estrogens produced centrally in response to social stimuli. The presence of an ER-mGluR interaction in birds suggests that this mechanism emerged relatively early in vertebrate history and is well conserved. Significance statement: The membrane-associated receptors underlying the acute effects of estrogens on behavior have mostly been documented in females, where active estrogens are thought to be of ovarian origin. Using acute intracerebroventricular injections of specific agonists and antagonists following blockade of brain aromatase, we show here that brain-derived estrogens acutely facilitate male sexual motivation through the activation of estrogen receptor β interacting with the metabotropic glutamate receptor 1a. This behavioral effect occurring within minutes provides a mechanistic explanation of how an estrogen receptor not intrinsically coupled to intracellular effectors can signal from the membrane to govern behavior in a very rapid fashion. It suggests that different subtypes of estrogen receptors could regulate the motivation versus performance aspects of behavior.

ERα36--Another piece of the estrogen puzzle

Although the nuclear action of estrogen receptors (ER) is a well-known fact, evidence supporting membrane estrogen receptors is steadily accumulating. New ER variants of unrecognized function have been discovered. ERα is a product of the ESR1 gene. It serves not only as a template for the full-length 66kDa protein, but also for smaller isoforms which exist as independent receptors. The recently discovered ERα36 (36kDa), consisting of 310 amino acids of total 595 ERα66 protein residues, is an example of that group. The transcription initiation site is identified in the first intron of the ESR1 gene. C-Terminal 27 amino acids are encoded by previously unknown exon 9. The presence of this unique C-terminal sequence creates an opportunity for the production of selective antibodies. ERα36 has been shown to have a high affinity to the cell membrane and as much as 90% of the protein can be bound with it. Post-translational palmitoylation is suspected to play a crucial role in ERα36 anchoring to the cell membrane. In silico analysis suggests the existence of a potential transmembrane domain in ERα36. ERα36 was found in most cells of animals at various ages, but its exact physiological function remains to be fully elucidated. It seems that cells traditionally considered as being deprived of ER are able to respond to hormonal stimulation via the ERα36 receptor. Moreover, ERα36 displays unique pharmacological properties and its action may be behind antiestrogen resistance. The use of ERα36 in cancer diagnosis gives rise to great expectations.

A molecular docking study of phytochemical estrogen mimics from dietary herbal supplements

PURPOSE

The purpose of this study is to use a molecular docking approach to identify potential estrogen mimics or anti-estrogens in phytochemicals found in popular dietary herbal supplements.

METHODS

In this study, 568 phytochemicals found in 17 of the most popular herbal supplements sold in the United States were built and docked with two isoforms of the estrogen receptor, ERα and ERβ (a total of 27 different protein crystal structures).

RESULTS

The docking results revealed six strongly docking compounds in Echinacea, three from milk thistle (Silybum marianum), three from Gingko biloba, one from Sambucus nigra, none from maca (Lepidium meyenii), five from chaste tree (Vitex agnus-castus), two from fenugreek (Trigonella foenum-graecum), and two from Rhodiola rosea. Notably, of the most popular herbal supplements for women, there were numerous compounds that docked strongly with the estrogen receptor: Licorice (Glycyrrhiza glabra) had a total of 26 compounds strongly docking to the estrogen receptor, 15 with wild yam (Dioscorea villosa), 11 from black cohosh (Actaea racemosa), eight from muira puama (Ptychopetalum olacoides or P. uncinatum), eight from red clover (Trifolium pratense), three from damiana (Turnera aphrodisiaca or T. diffusa), and three from dong quai (Angelica sinensis). Of possible concern were the compounds from men's herbal supplements that exhibited strong docking to the estrogen receptor: Gingko biloba had three compounds, gotu kola (Centella asiatica) had two, muira puama (Ptychopetalum olacoides or P. uncinatum) had eight, and Tribulus terrestris had six compounds.

CONCLUSIONS

This molecular docking study has revealed that almost all popular herbal supplements contain phytochemical components that may bind to the human estrogen receptor and exhibit selective estrogen receptor modulation. As such, these herbal supplements may cause unwanted side effects related to estrogenic activity.

Computational polypharmacology analysis of the heat shock protein 90 interactome

The design of a single drug molecule that is able to simultaneously and specifically interact with multiple biological targets is gaining major consideration in drug discovery. However, the rational design of drugs with a desired polypharmacology profile is still a challenging task, especially when these targets are distantly related or unrelated. In this work, we present a computational approach aimed at the identification of suitable target combinations for multitarget drug design within an ensemble of biologically relevant proteins. The target selection relies on the analysis of activity annotations present in molecular databases and on ligand-based virtual screening. A few target combinations were also inspected with structure-based methods to demonstrate that the identified dual-activity compounds are able to bind target combinations characterized by remote binding site similarities. Our approach was applied to the heat shock protein 90 (Hsp90) interactome, which contains several targets of key importance in cancer. Promising target combinations were identified, providing a basis for the computational design of compounds with dual activity. The approach may be used on any ensemble of proteins of interest for which known inhibitors are available.