Exploiting a subtype-specific mitochondrial vulnerability for successful treatment of colorectal peritoneal metastases.

Peritoneal metastases (PMs) from colorectal cancer (CRC) respond poorly to treatment and are associated with unfavorable prognosis. For example, the addition of hyperthermic intraperitoneal chemotherapy (HIPEC) to cytoreductive surgery in resectable patients shows limited benefit, and novel treatments are urgently needed. The majority of CRC-PMs represent the CMS4 molecular subtype of CRC, and here we queried the vulnerabilities of this subtype in pharmacogenomic databases to identify novel therapies. This reveals the copper ionophore elesclomol (ES) as highly effective against CRC-PMs. ES exhibits rapid cytotoxicity against CMS4 cells by targeting mitochondria. We find that a markedly reduced mitochondrial content in CMS4 cells explains their vulnerability to ES. ES demonstrates efficacy in preclinical models of PMs, including CRC-PMs and ovarian cancer organoids, mouse models, and a HIPEC rat model of PMs. The above proposes ES as a promising candidate for the local treatment of CRC-PMs, with broader implications for other PM-prone cancers.

GPCRDB: information system for G protein-coupled receptors.

The GPCRDB is a Molecular Class-Specific Information System (MCSIS) that collects, combines, validates and disseminates large amounts of heterogeneous data on G protein-coupled receptors (GPCRs). The GPCRDB contains experimental data on sequences, ligand-binding constants, mutations and oligomers, as well as many different types of computationally derived data such as multiple sequence alignments and homology models. The GPCRDB provides access to the data via a number of different access methods. It offers visualization and analysis tools, and a number of query systems. The data is updated automatically on a monthly basis. The GPCRDB can be found online at http://www.gpcr.org/7tm/.

Discovery of Targets for Immune-Metabolic Antitumor Drugs Identifies Estrogen-Related Receptor Alpha.

Drugs that kill tumors through multiple mechanisms have the potential for broad clinical benefits. Here, we first developed an in silico multiomics approach (BipotentR) to find cancer cell-specific regulators that simultaneously modulate tumor immunity and another oncogenic pathway and then used it to identify 38 candidate immune-metabolic regulators. We show the tumor activities of these regulators stratify patients with melanoma by their response to anti-PD-1 using machine learning and deep neural approaches, which improve the predictive power of current biomarkers. The topmost identified regulator, ESRRA, is activated in immunotherapy-resistant tumors. Its inhibition killed tumors by suppressing energy metabolism and activating two immune mechanisms: (i) cytokine induction, causing proinflammatory macrophage polarization, and (ii) antigen-presentation stimulation, recruiting CD8+ T cells into tumors. We also demonstrate a wide utility of BipotentR by applying it to angiogenesis and growth suppressor evasion pathways. BipotentR (http://bipotentr.dfci.harvard.edu/) provides a resource for evaluating patient response and discovering drug targets that act simultaneously through multiple mechanisms. SIGNIFICANCE: BipotentR presents resources for evaluating patient response and identifying targets for drugs that can kill tumors through multiple mechanisms concurrently. Inhibition of the topmost candidate target killed tumors by suppressing energy metabolism and effects on two immune mechanisms. This article is highlighted in the In This Issue feature, p. 517.

Comparative analysis of pharmacophore screening tools.

The pharmacophore concept is of central importance in computer-aided drug design (CADD) mainly because of its successful application in medicinal chemistry and, in particular, high-throughput virtual screening (HTVS). The simplicity of the pharmacophore definition enables the complexity of molecular interactions between ligand and receptor to be reduced to a handful set of features. With many pharmacophore screening softwares available, it is of the utmost interest to explore the behavior of these tools when applied to different biological systems. In this work, we present a comparative analysis of eight pharmacophore screening algorithms (Catalyst, Unity, LigandScout, Phase, Pharao, MOE, Pharmer, and POT) for their use in typical HTVS campaigns against four different biological targets by using default settings. The results herein presented show how the performance of each pharmacophore screening tool might be specifically related to factors such as the characteristics of the binding pocket, the use of specific pharmacophore features, and the use of these techniques in specific steps/contexts of the drug discovery pipeline. Algorithms with rmsd-based scoring functions are able to predict more compound poses correctly as overlay-based scoring functions. However, the ratio of correctly predicted compound poses versus incorrectly predicted poses is better for overlay-based scoring functions that also ensure better performances in compound library enrichments. While the ensemble of these observations can be used to choose the most appropriate class of algorithm for specific virtual screening projects, we remarked that pharmacophore algorithms are often equally good, and in this respect, we also analyzed how pharmacophore algorithms can be combined together in order to increase the success of hit compound identification. This study provides a valuable benchmark set for further developments in the field of pharmacophore search algorithms, e.g., by using pose predictions and compound library enrichment criteria.

Estrogen-related receptor alpha drives mitochondrial biogenesis and resistance to neoadjuvant chemoradiation in esophageal cancer.

Neoadjuvant chemoradiotherapy (nCRT) improves outcomes in resectable esophageal adenocarcinoma (EAC), but acquired resistance precludes long-term efficacy. Here, we delineate these resistance mechanisms. RNA sequencing on matched patient samples obtained pre-and post-neoadjuvant treatment reveal that oxidative phosphorylation was the most upregulated of all biological programs following nCRT. Analysis of patient-derived models confirms that mitochondrial content and oxygen consumption strongly increase in response to nCRT and that ionizing radiation is the causative agent. Bioinformatics identifies estrogen-related receptor alpha (ESRRA) as the transcription factor responsible for reprogramming, and overexpression and silencing of ESRRA functionally confirm that its downstream metabolic rewiring contributes to resistance. Pharmacological inhibition of ESRRA successfully sensitizes EAC organoids and patient-derived xenografts to radiation. In conclusion, we report a profound metabolic rewiring following chemoradiation and demonstrate that its inhibition resensitizes EAC cells to radiation. These findings hold broader relevance for other cancer types treated with radiation as well.

ss-TEA: Entropy based identification of receptor specific ligand binding residues from a multiple sequence alignment of class A GPCRs.

BACKGROUND: G-protein coupled receptors (GPCRs) are involved in many different physiological processes and their function can be modulated by small molecules which bind in the transmembrane (TM) domain. Because of their structural and sequence conservation, the TM domains are often used in bioinformatics approaches to first create a multiple sequence alignment (MSA) and subsequently identify ligand binding positions. So far methods have been developed to predict the common ligand binding residue positions for class A GPCRs. RESULTS: Here we present 1) ss-TEA, a method to identify specific ligand binding residue positions for any receptor, predicated on high quality sequence information. 2) The largest MSA of class A non olfactory GPCRs in the public domain consisting of 13324 sequences covering most of the species homologues of the human set of GPCRs. A set of ligand binding residue positions extracted from literature of 10 different receptors shows that our method has the best ligand binding residue prediction for 9 of these 10 receptors compared to another state-of-the-art method. CONCLUSIONS: The combination of the large multi species alignment and the newly introduced residue selection method ss-TEA can be used to rapidly identify subfamily specific ligand binding residues. This approach can aid the design of site directed mutagenesis experiments, explain receptor function and improve modelling. The method is also available online via GPCRDB at http://www.gpcr.org/7tm/.

Olfactomedin-4 regulation by estrogen in the human endometrium requires epidermal growth factor signaling.

Olfactomedin-4 (OLFM-4) is an extracellular matrix protein that is highly expressed in human endometrium. We have examined the regulation and function of OLFM-4 in normal endometrium and in cases of endometriosis and endometrial cancer. OLFM-4 expression levels are highest in proliferative-phase endometrium, and 17ß-estradiol up-regulates OLFM-4 mRNA in endometrial explant cultures. Using the luciferase reporter under control of the OLFM-4 promoter, it was shown that both 17ß-estradiol and OH-tamoxifen induce luciferase activity, and epidermal growth factor receptor-1 is required for this estrogenic response. In turn, EGF activates the OLFM-4 promoter, and estrogen receptor-α is needed for the complete EGF response. The cellular functions of OLFM-4 were examined by its expression in OLFM-4-negative HEK-293 cells, which resulted in decreased vimentin expression and cell adherence as well as increased apoptosis resistance. In cases of endometriosis and endometrial cancer, OLFM-4 expression correlated with the presence of epidermal growth factor receptor-1 and estrogen receptor-α (or estrogen signaling). An increase of OLFM-4 mRNA was observed in the endometrium of endometriosis patients. No change in OLFM-4 expression levels were observed in patients with endometrial cancer relative with controls. In conclusion, cross-talk between estrogen and EGF signaling regulates OLFM-4 expression. The role of OLFM-4 in endometrial tissue remodeling before the secretory phase and during the predisposition and early events in endometriosis can be postulated but requires additional investigation.

Snooker: a structure-based pharmacophore generation tool applied to class A GPCRs.

G-protein coupled receptors (GPCRs) are important drug targets for various diseases and of major interest to pharmaceutical companies. The function of individual members of this protein family can be modulated by the binding of small molecules at the extracellular side of the structurally conserved transmembrane (TM) domain. Here, we present Snooker, a structure-based approach to generate pharmacophore hypotheses for compounds binding to this extracellular side of the TM domain. Snooker does not require knowledge of ligands, is therefore suitable for apo-proteins, and can be applied to all receptors of the GPCR protein family. The method comprises the construction of a homology model of the TM domains and prioritization of residues on the probability of being ligand binding. Subsequently, protein properties are converted to ligand space, and pharmacophore features are generated at positions where protein ligand interactions are likely. Using this semiautomated knowledge-driven bioinformatics approach we have created pharmacophore hypotheses for 15 different GPCRs from several different subfamilies. For the beta-2-adrenergic receptor we show that ligand poses predicted by Snooker pharmacophore hypotheses reproduce literature supported binding modes for ∼75% of compounds fulfilling pharmacophore constraints. All 15 pharmacophore hypotheses represent interactions with essential residues for ligand binding as observed in mutagenesis experiments and compound selections based on these hypotheses are shown to be target specific. For 8 out of 15 targets enrichment factors above 10-fold are observed in the top 0.5% ranked compounds in a virtual screen. Additionally, prospectively predicted ligand binding poses in the human dopamine D3 receptor based on Snooker pharmacophores were ranked among the best models in the community wide GPCR dock 2010.

A prospective cross-screening study on G-protein-coupled receptors: lessons learned in virtual compound library design.

We present the systematic prospective evaluation of a protein-based and a ligand-based virtual screening platform against a set of three G-protein-coupled receptors (GPCRs): the ß-2 adrenoreceptor (ADRB2), the adenosine A(2A) receptor (AA2AR), and the sphingosine 1-phosphate receptor (S1PR1). Novel bioactive compounds were identified using a consensus scoring procedure combining ligand-based (frequent substructure ranking) and structure-based (Snooker) tools, and all 900 selected compounds were screened against all three receptors. A striking number of ligands showed affinity/activity for GPCRs other than the intended target, which could be partly attributed to the fuzziness and overlap of protein-based pharmacophore models. Surprisingly, the phosphodiesterase 5 (PDE5) inhibitor sildenafil was found to possess submicromolar affinity for AA2AR. Overall, this is one of the first published prospective chemogenomics studies that demonstrate the identification of novel cross-pharmacology between unrelated protein targets. The lessons learned from this study can be used to guide future virtual ligand design efforts.

Discovery and optimization of 1-(4-(pyridin-2-yl)benzyl)imidazolidine-2,4-dione derivatives as a novel class of selective cannabinoid CB2 receptor agonists.

Here, we report the identification and optimization of 1-(4-(pyridin-2-yl)benzyl)imidazolidine-2,4-dione derivatives as a novel chemotype with selective cannabinoid CB2 receptor agonist activity. 1 is a potent and selective cannabinoid CB2 receptor agonist (hCB2 pEC(50) = 8.6). The compound was found to be metabolically unstable, which resulted in low oral bioavailability in rat (F(po) = 4%) and possessed off-target activity at the hERG ion channel (pK(i) = 5.5). Systematic modification of physicochemical properties, such as lipophilicity and basicity, was used to optimize the pharmacokinetic profile and hERG affinity of this novel class of cannabinoid CB2 receptor agonists. This led to the identification of 44 as a potent, selective, and orally bioavailable cannabinoid CB2 receptor agonist (hCB2 pEC(50) = 8.0; hERG pK(i) < 4; F(po) = 100%), which was active in a rat spinal nerve ligation model of neuropathic pain.