Cheminformatics Driven Development of Novel Therapies for Drug Resistant Prostate Cancer

Exploring the Role of Chemoinformatics in Accelerating Drug Discovery: A Computational Approach

Cheminformatics and its role in drug discovery is expected to be the privileged approach in handling large number of chemical datasets. This approach contributes toward the pharmaceutical development and assessment of chemical compounds at a faster rate efficiently. Additionally, as technological advancement impacts research, cheminformatics is being used more and more in the field of health science. This chapter describes the concepts of cheminformatics along with its involvement in drug discovery with a case study.

Cheminformatic Analysis and Machine Learning Modeling to Investigate Androgen Receptor Antagonists to Combat Prostate Cancer

Prostate cancer (PCa) is a major leading cause of mortality of cancer among males. There have been numerous studies to develop antagonists against androgen receptor (AR), a crucial therapeutic target for PCa. This study is a systematic cheminformatic analysis and machine learning modeling to study the chemical space, scaffolds, structure-activity relationship, and landscape of human AR antagonists. There are 1678 molecules as final data sets. Chemical space visualization by physicochemical property visualization has demonstrated that molecules from the potent/active class generally have a mildly smaller molecular weight (MW), octanol-water partition coefficient (log P), number of hydrogen-bond acceptors (nHA), number of rotatable bonds (nRot), and topological polar surface area (TPSA) than molecules from intermediate/inactive class. The chemical space visualization in the principal component analysis (PCA) plot shows significant overlapping distributions between potent/active class molecules and intermediate/inactive class molecules; potent/active class molecules are intensively distributed, while intermediate/inactive class molecules are widely and sparsely distributed. Murcko scaffold analysis has shown low scaffold diversity in general, and scaffold diversity of potent/active class molecules is even lower than intermediate/inactive class molecules, indicating the necessity for developing molecules with novel scaffolds. Furthermore, scaffold visualization has identified 16 representative Murcko scaffolds. Among them, scaffolds 1, 2, 3, 4, 7, 8, 10, 11, 15, and 16 are highly favorable scaffolds due to their high scaffold enrichment factor values. Based on scaffold analysis, their local structure-activity relationships (SARs) were investigated and summarized. In addition, the global SAR landscape was explored by quantitative structure-activity relationship (QSAR) modelings and structure-activity landscape visualization. A QSAR classification model incorporating all of the 1678 molecules stands out as the best model from a total of 12 candidate models for AR antagonists (built on PubChem fingerprint, extra trees algorithm, accuracy for training set: 0.935, 10-fold cross-validation set: 0.735 and test set: 0.756). Deeper insights into the structure-activity landscape highlighted a total of seven significant activity cliff (AC) generators (ChEMBL molecule IDs: 160257, 418198, 4082265, 348918, 390728, 4080698, and 6530), which provide valuable SAR information for medicinal chemistry. The findings in this study provide new insights and guidelines for hit identification and lead optimization for the development of novel AR antagonists.

Allosteric binding on nuclear receptors: Insights on screening of non-competitive endocrine-disrupting chemicals

Endocrine-disrupting chemicals (EDCs) can compete with endogenous hormones and bind to the orthosteric site of nuclear receptors (NRs), affecting normal endocrine system function and causing severe symptoms. Recently, a series of pharmaceuticals and personal care products (PPCPs) have been discovered to bind to the allosteric sites of NRs and induce similar effects. However, it remains unclear how diverse EDCs work in this new way. Therefore, we have systematically summarized the allosteric sites and underlying mechanisms based on existing studies, mainly regarding drugs belonging to the PPCP class. Advanced methods, classified as structural biology, biochemistry and computational simulation, together with their advantages and hurdles for allosteric site recognition and mechanism insight have also been described. Furthermore, we have highlighted two available strategies for virtual screening of numerous EDCs, relying on the structural features of allosteric sites and lead compounds, respectively. We aim to provide reliable theoretical and technical support for a broader view of various allosteric interactions between EDCs and NRs, and to drive high-throughput and accurate screening of potential EDCs with non-competitive effects.

Interleukin-10 and Small Molecule SHIP1 Allosteric Regulators Trigger Anti-inflammatory Effects through SHIP1/STAT3 Complexes

The anti-inflammatory actions of interleukin-10 (IL10) are thought to be mediated primarily by the STAT3 transcription factor, but pro-inflammatory cytokines such as interleukin-6 (IL6) also act through STAT3. We now report that IL10, but not IL6 signaling, induces formation of a complex between STAT3 and the inositol polyphosphate-5-phosphatase SHIP1 in macrophages. Both SHIP1 and STAT3 translocate to the nucleus in macrophages. Remarkably, sesquiterpenes of the Pelorol family, which we previously described as allosteric activators of SHIP1 phosphatase activity, could induce SHIP1/STAT3 complex formation in cells and mimic the anti-inflammatory action of IL10 in a mouse model of colitis. Using crystallography and docking studies we identified a drug-binding pocket in SHIP1. Our studies reveal new mechanisms of action for both STAT3 and SHIP1 and provide a rationale for use of allosteric SHIP1-activating compounds, which mimic the beneficial anti-inflammatory actions of IL10.

Video Abstract

•

Loss of normal interleukin-10 (IL10) function results in inflammatory diseases

•

IL10 or SHIP1 agonists induce formation of SHIP1/STAT3 complexes

•

SHIP1 Y190 phosphorylation is required for SHIP1/STAT3 complex formation

•

SHIP1 agonists mimic IL10 anti-inflammatory action in a mouse model of colitis

Computer-Aided Ligand Discovery for Estrogen Receptor Alpha

Breast cancer (BCa) is one of the most predominantly diagnosed cancers in women. Notably, 70% of BCa diagnoses are Estrogen Receptor α positive (ERα+) making it a critical therapeutic target. With that, the two subtypes of ER, ERα and ERβ, have contrasting effects on BCa cells. While ERα promotes cancerous activities, ERβ isoform exhibits inhibitory effects on the same. ER-directed small molecule drug discovery for BCa has provided the FDA approved drugs tamoxifen, toremifene, raloxifene and fulvestrant that all bind to the estrogen binding site of the receptor. These ER-directed inhibitors are non-selective in nature and may eventually induce resistance in BCa cells as well as increase the risk of endometrial cancer development. Thus, there is an urgent need to develop novel drugs with alternative ERα targeting mechanisms that can overcome the limitations of conventional anti-ERα therapies. Several functional sites on ERα, such as Activation Function-2 (AF2), DNA binding domain (DBD), and F-domain, have been recently considered as potential targets in the context of drug research and discovery. In this review, we summarize methods of computer-aided drug design (CADD) that have been employed to analyze and explore potential targetable sites on ERα, discuss recent advancement of ERα inhibitor development, and highlight the potential opportunities and challenges of future ERα-directed drug discovery.

Feature vector clustering molecular pairs in computer simulations

A clustering framework is introduced to analyze the microscopic structural organization of molecular pairs in liquids and solutions. A molecular pair is represented by a representative vector (RV). To obtain RV, intermolecular atom distances in the pair are extracted from simulation trajectory as components of the key feature vector (KFV). A specific scheme is then suggested to transform KFV to RV by removing the influence of permutational molecular symmetry on the KFV as the predicted clusters should be independent of possible permutations of identical atoms in the pair. After RVs of pairs are obtained, a clustering analysis technique is finally used to classify all the RVs of molecular pairs into the clusters. The framework is applied to analyze trajectory from molecular dynamics simulations of an ionic liquid (trihexyltetradecylphosphonium bis(oxalato)borate ([P_6,6,6,14 ][BOB])). The molecular pairs are successfully categorized into physically meaningful clusters, and their effectiveness is evaluated by computing the product moment correlation coefficient (PMCC). (Willett, Winterman, and Bawden, J. Chem. Inf. Comput. Sci. 1986, 26, 109-118; Downs, Willett, and Fisanick, J. Chem. Inf. Comput. Sci. 1994, 34, 1094-1102) It is observed that representative configurations of two clusters are related to two energy local minimum structures optimized by density functional theory (DFT) calculation, respectively. Several widely used clustering analysis techniques of both nonhierarchical (k-means) and hierarchical clustering algorithms are also evaluated and compared with each other. The proposed KFV technique efficiently reveals local molecular pair structures in the simulated complex liquid. It is a method, which is highly useful for liquids and solutions in particular with strong intermolecular interactions. © 2019 Wiley Periodicals, Inc.

Molecules targeting the androgen receptor (AR) signaling axis beyond the AR-Ligand binding domain

Prostate cancer (PCa) is the second most common cause of cancer-related mortality in men in the United States. The androgen receptor (AR) and the physiological pathways it regulates are central to the initiation and progression of PCa. As a member of the nuclear steroid receptor family, it is a transcription factor with three distinct functional domains (ligand-binding domain [LBD], DNA-binding domain [DBD], and transactivation domain [TAD]) in its structure. All clinically approved drugs for PCa ultimately target the AR-LBD. Clinically active drugs that target the DBD and TAD have not yet been developed due to multiple factors. Despite these limitations, the last several years have seen a rise in the discovery of molecules that could successfully target these domains. This review aims to present and comprehensively discuss such molecules that affect AR signaling through direct or indirect interactions with the AR-TAD or the DBD. The compounds discussed here include hairpin polyamides, niclosamide, marine sponge-derived small molecules (eg, EPI compounds), mahanine, VPC compounds, JN compounds, and bromodomain and extraterminal domain inhibitors. We highlight the significant in vitro and in vivo data found for each compound and the apparent limitations and/or potential for further development of these agents as PCa therapies.