Development of leads targeting ER-α in breast cancer: An in silico exploration from natural domain

The Role of Simulation Science in Public Health at the Agency for Toxic Substances and Disease Registry: An Overview and Analysis of the Last Decade

Environmental exposures are ubiquitous and play a significant, and sometimes understated, role in public health as they can lead to the development of various chronic and infectious diseases. In an ideal world, there would be sufficient experimental data to determine the health effects of exposure to priority environmental contaminants. However, this is not the case, as emerging chemicals are continuously added to this list, furthering the data gaps. Recently, simulation science has evolved and can provide appropriate solutions using a multitude of computational methods and tools. In its quest to protect communities across the country from environmental health threats, ATSDR employs a variety of simulation science tools such as Physiologically Based Pharmacokinetic (PBPK) modeling, Quantitative Structure-Activity Relationship (QSAR) modeling, and benchmark dose (BMD) modeling, among others. ATSDR's use of such tools has enabled the agency to evaluate exposures in a timely, efficient, and effective manner. ATSDR's work in simulation science has also had a notable impact beyond the agency, as evidenced by external researchers' widespread appraisal and adaptation of the agency's methodology. ATSDR continues to advance simulation science tools and their applications by collaborating with researchers within and outside the agency, including other federal/state agencies, NGOs, the private sector, and academia.

Estrogenic flavonoids and their molecular mechanisms of action

Flavonoids are a major group of phytoestrogens associated with physiological effects, and ecological and social impacts. Although the estrogenic activity of flavonoids was reported by researchers in the fields of medical, environmental and food studies, their molecular mechanisms of action have not been comprehensively reviewed. The estrogenic activity of the respective classes of flavonoids, anthocyanidins/anthocyanins, 2-arylbenzofurans/3-arylcoumarins/α-methyldeoxybenzoins, aurones/chalcones/dihydrochalcones, coumaronochromones, coumestans, flavans/flavan-3-ols/flavan-4-ols, flavanones/dihydroflavonols, flavones/flavonols, homoisoflavonoids, isoflavans, isoflavanones, isoflavenes, isoflavones, neoflavonoids, oligoflavonoids, pterocarpans/pterocarpenes, and rotenone/rotenoids, was summarized through a comprehensive literature search, and their structure-activity relationship, biological activities, signaling pathways, and applications were discussed. Although the respective classes of flavonoids contained at least one chemical mimicking estrogen, the mechanisms varied, such as those with estrogenic, anti-estrogenic, non-estrogenic, and biphasic activities, and additional activities through crosstalk/bypassing, which exert biological activities through cell signaling pathways. Such mechanistic variations of estrogen action are not limited to flavonoids and are observed among other broad categories of chemicals, thus this group of chemicals can be termed as the "estrogenome". This review article focuses on the connection of estrogen action mainly between the outer and the inner environments, which represent variations of chemicals and biological activities/signaling pathways, respectively, and form the basis to understand their applications. The applications of chemicals will markedly progress due to emerging technologies, such as artificial intelligence for precision medicine, which is also true of the study of the estrogenome including estrogenic flavonoids.

Molecular docking, QSAR, pharmacophore modeling, and dynamics studies of some chromone derivatives for the discovery of anti-breast cancer agents against hormone-dependent breast cancer

In search of new anti-breast cancer agents, the present study envisaged the design and synthesis of a series of benzopyran-chalcones. All the synthesized compounds were assayed for their in-vitro anticancer activity against ER + MCF-7 and triple-negative MDA-MB-231 breast cancer cell lines using SRB assay. The synthesized compounds were found active against ER + MCF-7 cell lines. Based on the in-vitro data, in-silico analysis was performed using hormone-dependent breast cancer targets such as hER-α and aromatase because the compounds showed activity against MCF-7 cells and none was active against MDA-MB-231. The in-silico results supported the in-vitro anticancer activity suggesting the affinity of compounds toward hormone-dependant breast cancer. Compounds 4A1 to 4A3 were found to be most cytotoxic to MCF-7 cells with IC50 values of 31.87, 22.95, and 20.34 μg/ml, respectively (Doxorubicin IC50: <10 μg/ml). In addition, they showed the interactions with the amino acid residues of a binding cavity of an hER-α. Furthermore, quantitative structure-activity relationship (QSAR) studies were performed to reveal the vital structural features required for anticancer activity against breast cancer. Molecular dynamic simulation studies of hER-α and 4A3 in comparison with the raloxifene complex ensure the appropriate refinement of compounds in the dynamic system. Additionally, a generated pharmacophore model explored the essential pharmacophoric features of the synthesized scaffolds with respect to clinically used drug molecules for optimal hormone-dependant anti-breast cancer activity.Communicated by Ramaswamy H. Sarma.

Molecular modeling studies of benzothiophene-containing derivatives as promising selective estrogen receptor downregulators: a combination of 3D-QSAR, molecular docking and molecular dynamics simulations

Selective estrogen receptor downregulators (SERDs) for the treatment of positive breast cancer can act both as estrogen alpha receptor (ERα) antagonists and degraders. In this study, the optimal antagonist models (CoMFA-A, q² = 0.660, r² = 0.996; CoMSIA-A, q² = 0.728, r² = 0.992) and degrader models (CoMFA-D, q² = 0.850, r² = 0.996; CoMSIA-D, q² = 0.719, r² = 0.995) of a series of potent benzothiophene-containing SERDs were constructed to explore the three-dimensional quantitative structure-activity relationship. Internal and external validation indicated that all models exhibited good applicability, high predictive ability and robustness. Contour maps revealed the relationships between the essential structural features and antagonistic and degradation activities. Additionally, molecular docking, molecular dynamics and free energy calculation studies were further performed to investigate the detailed binding mode. Results indicated that several key residues, ARG394, GLU353, PHE404 and ILE424, were crucial for the stability of the ligand binding domain. The hydrophobic, electrostatic and Van der Waals interactions played significant effect on the binding affinity. Finally, ten novel compounds were designed based on above findings, where the predicted activity of compound D8 was equivalent to that of the compound LSZ102. 3D-QSAR, ADMET and bioavailability predictions indicated that all designed compounds with good predicted activity, good physicochemical and bioavailability could be potential candidates for SERDs. These results and combinations of computational methods provided guidance for the rational drug design of novel potential SERDs.Communicated by Ramaswamy H. Sarma.

Elucidating the aryl hydrocarbon receptor antagonism from a chemical-structural perspective

The aryl hydrocarbon receptor (AhR) plays an important role in several biological processes such as reproduction, immunity and homoeostasis. However, little is known on the chemical-structural and physicochemical features that influence the activity of AhR antagonistic modulators. In the present report, in vitro AhR antagonistic activity evaluations, based on a chemical-activated luciferase gene expression (AhR-CALUX) bioassay, and an extensive literature review were performed with the aim of constructing a structurally diverse database of contaminants and potentially toxic chemicals. Subsequently, QSAR models based on Linear Discriminant Analysis and Logistic Regression, as well as two toxicophoric hypotheses were proposed to model the AhR antagonistic activity of the built dataset. The QSAR models were rigorously validated yielding satisfactory performance for all classification parameters. Likewise, the toxicophoric hypotheses were validated using a diverse set of 350 decoys, demonstrating adequate robustness and predictive power. Chemical interpretations of both the QSAR and toxicophoric models suggested that hydrophobic constraints, the presence of aromatic rings and electron-acceptor moieties are critical for the AhR antagonism. Therefore, it is hoped that the deductions obtained in the present study will contribute to elucidate further on the structural and physicochemical factors influencing the AhR antagonistic activity of chemical compounds.

Discovery of novel natural compound inhibitors targeting estrogen receptor α by an integrated virtual screening strategy

Estrogen receptor (ER) is a nuclear hormone receptor and plays an important role in mediating the cellular effects of estrogen. ER can be classified into two receptors: estrogen receptor alpha (ERα) and beta (ERβ), and the former is expressed in 50~80% of breast tumors and has been extensively investigated in breast cancer for decades. Excessive exposure to estrogen can obviously stimulate the growth of breast cancers primarily mediated by ERα, and thus anti-estrogen therapies by small molecules are of concern to clinicians and pharmaceutical industry in the treatment of ERα-positive breast cancers. Although a series of estrogen receptor modulators have been developed, these drugs can lead to resistance and side effects. Therefore, the development of small molecule inhibitors with high target specificity has been intensified. In this pursuit, an integrated computer-aided virtual screening technique, including molecular docking and pharmacophore model screening, was used to screen traditional Chinese medicine (TCM) databases. The compounds with high docking scores and fit values were subjected to ADME (adsorption, distribution, metabolism, excretion) and toxicity prediction, and ten hits were identified as potential inhibitors targeting ERα. Molecular docking was used to investigate the binding modes between ERα and three most potent hits, and molecular dynamic simulations were chosen to explore the stability of these complexes. The rank of the predicted binding free energies evaluated by MM/GBSA is consistent with the docking score. These novel scaffolds discovered in the present study can be used as critical starting point in the drug discovery process for treating ERα-positive breast cancer. Graphical abstract .