Discovery of estrogen receptor modulators: a review of virtual screening and SAR efforts

Estrogen Receptor Alpha Binders for Hormone-Dependent Forms of Breast Cancer: e-QSAR and Molecular Docking Supported by X-ray Resolved Structures

Cancer, a life-disturbing and lethal disease with a high global impact, causes significant economic, social, and health challenges. Breast cancer refers to the abnormal growth of cells originating from breast tissues. Hormone-dependent forms of breast cancer, such as those influenced by estrogen, prompt the exploration of estrogen receptors as targets for potential therapeutic interventions. In this study, we conducted e-QSAR molecular docking and molecular dynamics analyses on a diverse set of inhibitors targeting estrogen receptor alpha (ER-α). The e-QSAR model is based on a genetic algorithm combined with multilinear regression analysis. The newly developed model possesses a balance between predictive accuracy and mechanistic insights adhering to the OECD guidelines. The e-QSAR model pointed out that sp2-hybridized carbon and nitrogen atoms are important atoms governing binding profiles. In addition, a specific combination of H-bond donors and acceptors with carbon, nitrogen, and ring sulfur atoms also plays a crucial role. The results are supported by molecular docking, MD simulations, and X-ray-resolved structures. The novel results could be useful for future drug development for ER-α.

The Effect of 42-Day Exposure to a Low Deoxynivalenol Dose on the Immunohistochemical Expression of Intestinal ERs and the Activation of CYP1A1 and GSTP1 Genes in the Large Intestine of Pre-pubertal Gilts

Deoxynivalenol (DON) is a mycotoxin that contaminates various plant materials. Exposure to DON can disrupt hormonal homeostasis, decrease body weight gains and modulate the immune system in pigs. It can also cause diarrhea, vomiting, leukocytosis, hemorrhaging or even death. Prolonged exposure to low doses of DON can have serious health implications in mammals. This is the first in vivo study to show that per os administration of low DON doses probably contributes to specific dysfunctions in steroidogenesis processes by inducing the immunohistochemical expression of estrogen receptors alpha (ERα) in the entire gastrointestinal tract in strongly stained cells (3 points) and estrogen receptors beta (ERβ), but only in both investigated segments of the duodenum in pre-pubertal gilts. Therefore, the aim of this study was to determine whether a NOAEL dose of DON (12 μg DON/kg BW) administered per os over a period of 42 days induces changes in the immunohistochemical expression of ER in different intestinal segments and the transcriptional activation of CYP1A1 and GSTP1 genes in the large intestine of pre-pubertal gilts. This is the first report to demonstrate the expression of ER, in particular ERβ, with the associated consequences. The expression of ER was accompanied by considerable variations in the activation of CYP1A1 and GSTP1 genes, but it supported the maintenance of a stable consensus between the degree of mycotoxin exposure and the detoxifying effect in pre-pubertal gilts.

Glycosides and Glycoconjugates of the Diterpenoid Isosteviol with a 1,2,3-Triazolyl Moiety: Synthesis and Cytotoxicity Evaluation

Several glycoconjugates of the diterpenoid isosteviol (16-oxo-ent-beyeran-19-oic acid) with a 1,2,3-triazolyl moiety were synthesized, and their cytotoxicity was evaluated against some human cancer and normal cell lines. Most of the synthesized compounds demonstrated weak inhibitory activities against the M-HeLa and MCF-7 human cancer cell lines. Three lead compounds, 54, 56 and 57, exhibited high selective cytotoxic activity against M-HeLa cells (IC₅₀ = 1.7-1.9 μM) that corresponded to the activity of the anticancer drug doxorubicin (IC₅₀ = 3.0 μM). Moreover, the lead compounds were not cytotoxic with respect to a Chang liver human normal cell line (IC₅₀ > 100 μM), whereas doxorubicin was cytotoxic to this cell line (IC₅₀ = 3.0 μM). It was found that cytotoxic activity of the lead compounds is due to induction of apoptosis proceeding along the mitochondrial pathway. The present findings suggest that 1,2,3-triazolyl-ring-containing glycoconjugates of isosteviol are a promising scaffold for the design of novel anticancer agents.

Insight into the structural stability of coumestrol with human estrogen receptor α and β subtypes: a combined approach involving docking and molecular dynamics simulation studies

Structural stability of coumestrol with human estrogen receptor α and β subtypes.

Versatility or promiscuity: the estrogen receptors, control of ligand selectivity and an update on subtype selective ligands

The estrogen receptors (ERs) are a group of versatile receptors. They regulate an enormity of processes starting in early life and continuing through sexual reproduction, development, and end of life. This review provides a background and structural perspective for the ERs as part of the nuclear receptor superfamily and discusses the ER versatility and promiscuity. The wide repertoire of ER actions is mediated mostly through ligand-activated transcription factors and many DNA response elements in most tissues and organs. Their versatility, however, comes with the drawback of promiscuous interactions with structurally diverse exogenous chemicals with potential for a wide range of adverse health outcomes. Even when interacting with endogenous hormones, ER actions can have adverse effects in disease progression. Finally, how nature controls ER specificity and how the subtle differences in receptor subtypes are exploited in pharmaceutical design to achieve binding specificity and subtype selectivity for desired biological response are discussed. The intent of this review is to complement the large body of literature with emphasis on most recent developments in selective ER ligands.

Alternative (non-animal) methods for cosmetics testing: current status and future prospects-2010

The 7th amendment to the EU Cosmetics Directive prohibits to put animal-tested cosmetics on the market in Europe after 2013. In that context, the European Commission invited stakeholder bodies (industry, non-governmental organisations, EU Member States, and the Commission's Scientific Committee on Consumer Safety) to identify scientific experts in five toxicological areas, i.e. toxicokinetics, repeated dose toxicity, carcinogenicity, skin sensitisation, and reproductive toxicity for which the Directive foresees that the 2013 deadline could be further extended in case alternative and validated methods would not be available in time. The selected experts were asked to analyse the status and prospects of alternative methods and to provide a scientifically sound estimate of the time necessary to achieve full replacement of animal testing. In summary, the experts confirmed that it will take at least another 7-9 years for the replacement of the current in vivo animal tests used for the safety assessment of cosmetic ingredients for skin sensitisation. However, the experts were also of the opinion that alternative methods may be able to give hazard information, i.e. to differentiate between sensitisers and non-sensitisers, ahead of 2017. This would, however, not provide the complete picture of what is a safe exposure because the relative potency of a sensitiser would not be known. For toxicokinetics, the timeframe was 5-7 years to develop the models still lacking to predict lung absorption and renal/biliary excretion, and even longer to integrate the methods to fully replace the animal toxicokinetic models. For the systemic toxicological endpoints of repeated dose toxicity, carcinogenicity and reproductive toxicity, the time horizon for full replacement could not be estimated.

Studies of tricyclic piperazine/piperidine furnished molecules as novel integrin αvβ3/αIIbβ3 dual antagonists using 3D-QSAR and molecular docking

The development of injectable integrin α(v)β(3)/α(IIb)β(3) dual antagonists attracts much attention of research for treating of acute ischemic diseases in recent years. In this work, based on a dataset composed of 102 tricyclic piperazine/piperidine furnished dual α(v)β(3) and α(IIb)β(3) antagonists, a variety of in silico modeling approaches including the comparative molecular field analysis (CoMFA), comparative similarity indices analysis (CoMSIA), and molecular docking were applied to reveal the requisite 3D structural features impacting the biological activities. Our statistical results show that the ligand-based 3D-QSAR models for both the α(v)β(3) and α(IIb)β(3) studies exhibited satisfactory internal and external predictability, i.e., for the CoMFA models, results of Q(2)=0.48, R(ncv)(2)=0.87, R(pred)(2)=0.71 for α(v)β(3) and Q(2)=0.50, R(ncv)(2)=0.85, R(pred)(2)=0.72 for α(IIb)β(3) analysis were obtained, and for the CoMSIA ones, the outcomes of Q(2)=0.55, R(ncv)(2)=0.90, R(pred)(2)=0.72 for α(v)β(3) and Q(2)=0.52, R(ncv)(2)=0.88, R(pred)(2)=0.74 for α(IIb)β(3) were achieved respectively. In addition, through a comparison between 3D-QSAR contour maps and docking results, it is revealed that that the most crucial interactions occurring between the tricyclic piperazine/piperidine derivatives and α(v)β(3)/α(IIb)β(3) receptor ligand binding pocket are H-bonding, and the key amino acids impacting the interactions are Arg214, Asn215, Ser123, and Lys253 for α(v)β(3), but Arg214, Asn215, Ser123 and Tyr190 for α(IIb)β(3) receptors, respectively. Halogen-containing groups at position 15 and 16, benzene sulfonamide substituent at position 23, and the replacement of piperazine with 4-aminopiperidine of ring B may increase the α(v)β(3)/α(IIb)β(3) antagonistic activity. The potencies for antagonists to inhibit isolated α(v)β(3) and α(IIb)β(3) are linear correlated, indicating that similar interaction mechanisms may exist for the series of molecules. To our best knowledge this is the first report on 3D-QSAR modeling of these dual α(v)β(3)/α(IIb)β(3) antagonists. The results obtained should provide information for better understanding of the mechanism of antagonism and thus be helpful in design of novel potent dual α(v)β(3)/α(IIb)β(3) antagonists.

3D-QSAR and molecular docking studies on derivatives of MK-0457, GSK1070916 and SNS-314 as inhibitors against Aurora B kinase

Development of anticancer drugs targeting Aurora B, an important member of the serine/threonine kinases family, has been extensively focused on in recent years. In this work, by applying an integrated computational method, including comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA), homology modeling and molecular docking, we investigated the structural determinants of Aurora B inhibitors based on three different series of derivatives of 108 molecules. The resultant optimum 3D-QSAR models exhibited (q² = 0.605, r²_pred = 0.826), (q² = 0.52, r²_pred = 0.798) and (q² = 0.582, r²_pred = 0.971) for MK-0457, GSK1070916 and SNS-314 classes, respectively, and the 3D contour maps generated from these models were analyzed individually. The contour map analysis for the MK-0457 model revealed the relative importance of steric and electrostatic effects for Aurora B inhibition, whereas, the electronegative groups with hydrogen bond donating capacity showed a great impact on the inhibitory activity for the derivatives of GSK1070916. Additionally, the predictive model of the SNS-314 class revealed the great importance of hydrophobic favorable contour, since hydrophobic favorable substituents added to this region bind to a deep and narrow hydrophobic pocket composed of residues that are hydrophobic in nature and thus enhanced the inhibitory activity. Moreover, based on the docking study, a further comparison of the binding modes was accomplished to identify a set of critical residues that play a key role in stabilizing the drug-target interactions. Overall, the high level of consistency between the 3D contour maps and the topographical features of binding sites led to our identification of several key structural requirements for more potency inhibitors. Taken together, the results will serve as a basis for future drug development of inhibitors against Aurora B kinase for various tumors.

In silico prediction of estrogen receptor subtype binding affinity and selectivity using statistical methods and molecular docking with 2-arylnaphthalenes and 2-arylquinolines

Over the years development of selective estrogen receptor (ER) ligands has been of great concern to researchers involved in the chemistry and pharmacology of anticancer drugs, resulting in numerous synthesized selective ER subtype inhibitors. In this work, a data set of 82 ER ligands with ERα and ERβ inhibitory activities was built, and quantitative structure-activity relationship (QSAR) methods based on the two linear (multiple linear regression, MLR, partial least squares regression, PLSR) and a nonlinear statistical method (Bayesian regularized neural network, BRNN) were applied to investigate the potential relationship of molecular structural features related to the activity and selectivity of these ligands. For ERα and ERβ, the performances of the MLR and PLSR models are superior to the BRNN model, giving more reasonable statistical properties (ERα: for MLR, R_tr² = 0.72, Q_te² = 0.63; for PLSR, R_tr² = 0.92, Q_te² = 0.84. ERβ: for MLR, R_tr² = 0.75, Q_te² = 0.75; for PLSR, R_tr² = 0.98, Q_te² = 0.80). The MLR method is also more powerful than other two methods for generating the subtype selectivity models, resulting in R_tr² = 0.74 and Q_te² = 0.80. In addition, the molecular docking method was also used to explore the possible binding modes of the ligands and a relationship between the 3D-binding modes and the 2D-molecular structural features of ligands was further explored. The results show that the binding affinity strength for both ERα and ERβ is more correlated with the atom fragment type, polarity, electronegativites and hydrophobicity. The substitutent in position 8 of the naphthalene or the quinoline plane and the space orientation of these two planes contribute the most to the subtype selectivity on the basis of similar hydrogen bond interactions between binding ligands and both ER subtypes. The QSAR models built together with the docking procedure should be of great advantage for screening and designing ER ligands with improved affinity and subtype selectivity property.