Search

Supplementary Methods

Pheochromocytomas and paragangliomas (PPGLs) provide some of the clearest genetic evidence for the critical role of metabolism in the tumorigenesis process. Approximately 40% of PPGLs are caused by driver germline mutations in 16 known susceptibility genes, and approximately half of these genes encode members of the tricarboxylic acid (TCA) cycle. Taking as a starting point the involvement of the TCA cycle in PPGL development, we aimed to identify unreported mutations that occurred in genes involved in this key metabolic pathway and that could explain the phenotypes of additional individuals who lack mutations in known susceptibility genes. To accomplish this, we applied a targeted sequencing of 37 TCA-cycle-related genes to DNA from 104 PPGL-affected individuals with no mutations in the major known predisposing genes. We also performed omics-based analyses, TCA-related metabolite determination, and (13)C(5)-glutamate labeling assays. We identified five germline variants affecting DLST in eight unrelated individuals (∼7%); all except one were diagnosed with multiple PPGLs. A recurrent variant, c.1121G>A (p.Gly374Glu), found in four of the eight individuals triggered accumulation of 2-hydroxyglutarate, both in tumors and in a heterologous cell-based assay designed to functionally evaluate DLST variants. p.Gly374Glu-DLST tumors exhibited loss of heterozygosity, and their methylation and expression profiles are similar to those of EPAS1-mutated PPGLs; this similarity suggests a link between DLST disruption and pseudohypoxia. Moreover, we found positive DLST immunostaining exclusively in tumors carrying TCA-cycle or EPAS1 mutations. In summary, this study reveals DLST as a PPGL-susceptibility gene and further strengthens the relevance of the TCA cycle in PPGL development.

Aberrant metabolism of cancer cells is well appreciated, but the identification of cancer subsets with specific metabolic vulnerabilities remains challenging. We conducted a chemical biology screen and identified a subset of neuroendocrine tumors displaying a striking pattern of sensitivity to inhibition of the cholesterol biosynthetic pathway enzyme squalene epoxidase (SQLE). Using a variety of orthogonal approaches, we demonstrate that sensitivity to SQLE inhibition results not from cholesterol biosynthesis pathway inhibition, but rather surprisingly from the specific and toxic accumulation of the SQLE substrate, squalene. These findings highlight SQLE as a potential therapeutic target in a subset of neuroendocrine tumors, particularly small cell lung cancers., Cancer cells are metabolically adaptable and the identification of specific vulnerabilities is challenging. Here the authors identify a subset of neuroendocrine cell lines exquisitely sensitive to inhibition of SQLE, an enzyme in the cholesterol biosynthetic pathway, due to the toxic accumulation of pathway intermediate squalene.

Somatic point mutations at a key arginine residue (R132) within the active site of the metabolic enzyme isocitrate dehydrogenase 1 (IDH1) confer a novel gain of function in cancer cells, resulting in the production of d-2-hydroxyglutarate (2-HG), an oncometabolite. Elevated 2-HG levels are implicated in epigenetic alterations and impaired cellular differentiation. IDH1 mutations have been described in an array of hematologic malignancies and solid tumors. Here, we report the discovery of AG-120 (ivosidenib), an inhibitor of the IDH1 mutant enzyme that exhibits profound 2-HG lowering in tumor models and the ability to effect differentiation of primary patient AML samples ex vivo. Preliminary data from phase 1 clinical trials enrolling patients with cancers harboring an IDH1 mutation indicate that AG-120 has an acceptable safety profile and clinical activity.

Somatic gain-of-function mutations in isocitrate dehydrogenases (IDH) 1 and 2 are found in multiple hematologic and solid tumors, leading to accumulation of the oncometabolite (R)-2-hydroxyglutarate (2HG). 2HG competitively inhibits α-ketoglutarate–dependent dioxygenases, including histone demethylases and methylcytosine dioxygenases of the TET family, causing epigenetic dysregulation and a block in cellular differentiation. In vitro studies have provided proof of concept for mutant IDH inhibition as a therapeutic approach. We report the discovery and characterization of AG-221, an orally available, selective, potent inhibitor of the mutant IDH2 enzyme. AG-221 suppressed 2HG production and induced cellular differentiation in primary human IDH2 mutation–positive acute myeloid leukemia (AML) cells ex vivo and in xenograft mouse models. AG-221 also provided a statistically significant survival benefit in an aggressive IDH2R140Q-mutant AML xenograft mouse model. These findings supported initiation of the ongoing clinical trials of AG-221 in patients with IDH2 mutation–positive advanced hematologic malignancies. Significance: Mutations in IDH1/2 are identified in approximately 20% of patients with AML and contribute to leukemia via a block in hematopoietic cell differentiation. We have shown that the targeted inhibitor AG-221 suppresses the mutant IDH2 enzyme in multiple preclinical models and induces differentiation of malignant blasts, supporting its clinical development. Cancer Discov; 7(5); 478–93. ©2017 AACR. See related commentary by Thomas and Majeti, p. 459. See related article by Shih et al., p. 494. This article is highlighted in the In This Issue feature, p. 443

SummaryAlthough aberrant metabolism in tumors has been well described, the identification of cancer subsets with particular metabolic vulnerabilities has remained challenging. Here, we conducted an siRNA screen focusing on enzymes involved in the tricarboxylic acid (TCA) cycle and uncovered a striking range of cancer cell dependencies on OGDH, the E1 subunit of the alpha-ketoglutarate dehydrogenase complex. Using an integrative metabolomics approach, we identified differential aspartate utilization, via the malate-aspartate shuttle, as a predictor of whether OGDH is required for proliferation in 3D culture assays and for the growth of xenograft tumors. These findings highlight an anaplerotic role of aspartate and, more broadly, suggest that differential nutrient utilization patterns can identify subsets of cancers with distinct metabolic dependencies for potential pharmacological intervention.

Small molecule modulators of GPR88 activity (agonists, antagonists, or modulators) are of interest as potential agents for the treatment of a variety of psychiatric disorders including schizophrenia. A series of phenylglycinol and phenylamine analogs have been prepared and evaluated for their GPR88 agonist activity and pharmacokinetic (PK) properties.

Cancer cells must carefully regulate their metabolism to maintain growth and division under varying nutrient and oxygen levels. Compelling data support the investigation of numerous enzymes as therapeutic targets to exploit metabolic vulnerabilities common to several cancer types. We discuss the rationale for developing such drugs and review three targets with central roles in metabolic pathways crucial for cancer cell growth: pyruvate kinase muscle isozyme splice variant 2 (PKM2) in glycolysis, glutaminase in glutaminolysis, and mutations in isocitrate dehydrogenase 1 and 2 isozymes (IDH1/2) in the tricarboxylic acid cycle. These targets exemplify the drugging approach to cancer metabolism, with allosteric modulation being the common theme. The first glutaminase and mutant IDH1/2 inhibitors have entered clinical testing, and early data are promising. Cancer metabolism provides a wealth of novel targets, and targeting allosteric sites promises to yield selective drugs with the potential to transform clinical outcomes across many cancer types.

AG-221 or enasidenib is a first-in-class selective inhibitor of mutated isocitrate dehydrogenase 2 (IDH2) with early demonstrated clinical efficacy in acute myeloid leukemia as a single agent, yet with persistence of mutant IDH2 clones. Two papers in this issue of Cancer Discovery provide further insight into the biological activity of AG-221 in promoting differentiation of IDH2 mutant cells and reversing aberrant DNA methylation over time, and demonstrating pre-clinical activity in combination with a targeted FLT3 kinase inhibitor to eliminate IDH2 mutant clones.

SummaryProliferating tumor cells use aerobic glycolysis to support their high metabolic demands. Paradoxically, increased glycolysis is often accompanied by expression of the lower activity PKM2 isoform, effectively constraining lower glycolysis. Here, we report the discovery of PKM2 activators with a unique allosteric binding mode. Characterization of how these compounds impact cancer cells revealed an unanticipated link between glucose and amino acid metabolism. PKM2 activation resulted in a metabolic rewiring of cancer cells manifested by a profound dependency on the nonessential amino acid serine for continued cell proliferation. Induction of serine auxotrophy by PKM2 activation was accompanied by reduced carbon flow into the serine biosynthetic pathway and increased expression of high affinity serine transporters. These data support the hypothesis that PKM2 expression confers metabolic flexibility to cancer cells that allows adaptation to nutrient stress.

Numerous reports have described the differential metabolism of cancer cells as compared to their normal counterparts. However, only relatively few metabolic genes with cancer-specific mutations have been reported and the identification of cancer subsets with particular metabolic vulnerabilities remains a challenge. To explore potential cancer-specific dependencies, we conducted a chemical biology screen utilizing a collection of small molecule inhibitors targeting diverse metabolic pathways in a large panel of cancer cell lines. A subset of neuroendocrine tumors, particularly small cell lung cancers (SCLC), displayed a striking dependence on squalene epoxygenase, SQLE, an enzyme in the cholesterol biosynthetic pathway. To develop further confidence in these findings, we have determined the first three-dimensional SQLE structure and further advanced a pharmacological toolbox for SQLE. Using these tools, we showed that the observed effects are on target and that the patterns of cellular sensitivity observed in vitro display excellent translation to in vivo xenografts studies. Interestingly, using a variety of orthogonal approaches, we demonstrated that SQLE sensitivity appears not to be related to overall inhibition of the cholesterol pathway but rather to specific and toxic accumulation of the SQLE substrate, squalene. Collectively, these findings highlight the utility of chemical biology screens and identify SQLE as a potential therapeutic target in a subset of neuroendocrine tumors, particularly SCLC. Citation Format: Christopher Mahoney, David Pirman, Victor Chubukov, Taryn Sleger, Anil Padyana, Stefan Gross, Sebastian Hayes, Zi Peng Fan, Gabrielle McDonald, Yu Chen, Joshua Murtie, Giovanni Cianchetta, Raj Nagaraja, Rohini Narayanaswamy, Sung Choe, Stuart Murray, Shengfang Jin, Scott Biller, Thomas Roddy, Gromoslaw A. Smolen. A chemical biology screen identifies a unique vulnerability of neuroendocrine cancer cells to SQLE inhibition [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 3504.

Trytophan Hydroxylase Type I (TPH1), most abundantly expressed in the gastrointestinal tract, initiates the synthesis of serotonin by catalyzing hydroxylation of tryptophan in the presence of biopterin and oxygen. We have previously described three series of novel, periphery-specific TPH1 inhibitors that selectively deplete serotonin in the gastrointestinal tract. We have now determined co-crystal structures of TPH1 with three of these inhibitors at high resolution. Analysis of the structural data showed that each of the three inhibitors fills the tryptophan binding pocket of TPH1 without reaching into the binding site of the cofactor pterin, and induces major conformational changes of the enzyme. The enzyme-inhibitor complexes assume a compact conformation that is similar to the one in tryptophan complex. Kinetic analysis showed that all three inhibitors are competitive versus the substrate tryptophan, consistent with the structural data that the compounds occupy the tryptophan binding site. On the other hand, all three inhibitors appear to be uncompetitive versus the cofactor 6-methyltetrahydropterin, which is not only consistent with the structural data but also indicate that the hydroxylation reaction follows an ordered binding mechanism in which a productive complex is formed only if tryptophan binds only after pterin, similar to the kinetic mechanisms of tyrosine and phenylalanine hydroxylase.

Tryptophan hydroxylase (TPH) is a key enzyme in the synthesis of serotonin. As a neurotransmitter, serotonin plays important physiological roles both peripherally and centrally. Here we describe the discovery of substituted triazines as a novel class of tryptophan hydroxylase inhibitors. This class of TPH inhibitors can selectively reduce serotonin levels in murine intestine after oral administration without affecting levels in the brain. These TPH inhibitors may provide novel treatments for gastrointestinal disorders associated with dysregulation of the serotonergic system, such as chemotherapy-induced emesis and irritable bowel syndrome.

We report here a general method for the prediction of hERG potassium channel blockers using computational models generated from correlation analyses of a large dataset and pharmacophore-based GRIND descriptors. These 3D-QSAR models are compared favorably with other traditional and chemometric based HQSAR methods.

Trying to understand the complex interactions that substrates and inhibitors have with the efflux transporter P-glycoprotein has been the subject of various publications. In this work, we have confined our study to substrates by picking a diverse set of 129 compounds based on the efflux ratios from Caco-2 permeability measurements. These compounds were then evaluated for P-glycoprotein inhibition using a calcein-AM assay. The subsequent data was used in a 3D-QSAR analysis using GRIND pharmacophore-based and physicochemical descriptors. Pharmacophore-based descriptors produced a much more robust model than the one obtained from physicochemical-based descriptors. This supports the process proposed by Seelig and co-workers previously published whereby the substrate enters the membrane as the first step and is then recognized by P-glycoprotein in a second step. The strong correlation, highlighted by PLS statistical analysis, between pharmacophoric descriptors and inhibition values suggests that substrate interaction, with perhaps the mouth of the protein or another binding site, plays a key role in the efflux process, yielding a model in which diffusion across the membrane is less important than substrate-protein interaction. One pharmacophore emerged from the analysis of the model. We pose that the recognition elements, at least determined by the molecules used in this study, are two hydrophobic groups 16.5 A apart and two hydrogen-bond-acceptor groups 11.5 A apart and that the dimensions of the molecule also plays a role in its recognition as a substrate.

The first allosteric, type III inhibitor of LIM-kinase 2 (LIMK2) is reported. A series of molecules that feature both an N-phenylsulfonamide and tertiary amide were not only very potent at LIMK2 but also were extremely selective against a panel of other kinases. Enzymatic kinetic studies showed these molecules to be noncompetitive with ATP, suggesting allosteric inhibition. X-ray crystallography confirmed that these sulfonamides are a rare example of a type III kinase inhibitor that binds away from the highly conserved hinge region and instead resides in the hydrophobic pocket formed in the DFG-out conformation of the kinase, thus accounting for the high level of selectivity observed.

Since the introduction of NMR prediction software, medicinal chemists have imagined submitting their compounds to corporate compound registration systems that would ultimately display a simplified pass/fail result. We initially implemented such a system based on HPLC and liquid chromatography mass spectrometry (LCMS) data that is embedded within our industry standard sample submission and registration process. By using gradient-heteronuclear single quantum coherence (HSQC) experiments, we have extended this concept to NMR data through a comparison of experimentally acquired data against predicted 1H and 13C NMR data. Integration of our compound registration system with our analytical instruments now provides our chemists unattended and automated NMR verification for collections of submitted compounds. The benefits achieved from automated processing and interpretation of results produced enhanced confidence in our compound library and released the chemists from the tedium of manipulating large amounts of data. This allows scientists to focus more of their attention to the drug discovery process. Copyright © 2008 John Wiley & Sons, Ltd.

Ligand-based virtual screening approaches were applied to search for new chemotype KCOs activating Kir6.2/SUR1 KATP channels. A total of 65 208 commercially available compounds, extracted from the ZINC archive, served as database for screening. In a first step, pharmacokinetic filtering via VolSurf reduced the initial database to 1913 compounds. Afterward, six molecules were selected as templates for similarity searches: similarity scores, obtained toward these templates, were calculated with the GRIND, FLAP, and TOPP approaches, which differently encode structural information into potential pharmacophores. In this way, we obtained 32 hit candidates, 16 via GRIND and eight each via FLAP and TOPP. For biological testing of the hit candidates, their effects on membrane potentials in HEK 293 cells expressing Kir6.2/SUR1 were studied. GRIND, FLAP, and TOPP all yielded hits, but no method top-ranked all the actives. Thus, parallel application of different approaches probably improves hit detection.

An extended VolSurf approach, that additionally includes SHAPE descriptors, was applied to a dataset of 55 quinolones. Bactericidal activity was measured at Bayer AG, Germany, for Gram-negative (Escherichia coli and Pseudomonas aeruginosa) and Gram-positive bacteria (Staphylococcus aureus and Enterococcus faecalis). Chemometric analysis was first approached via a classical VolSurf approach. The following descriptors were found most important: bactericidal activity particularly increases with high values of the best volume (BV11(OH2)) and the minimum energy (Emin1(OH2)) of the water probe, high values of the integy moment (ID(DRY)) of the lipophilic probe, and high values of the hydrophilic region (W(O)) of the hydrogen bond acceptor probe. Best volume (BV31(OH2)) of the water probe and best volume (BV12(DRY)) and lipophilic regions (D(DRY)) of the lipophilic probe as well as H-bonding capacity derived with the CO probe (HB(O)) are inversely related to activity. PLS analysis yields a five-component model with an r(2) of 0.83 and a q(2) of 0.43 after variable selection via fractional factorial design (FFD). Chemometric modeling could be improved by including newly derived SHAPE descriptors, which were merged with the VolSurf descriptors and subjected to PLS analysis. The global model of this extended VolSurf approach is optimal with two components and exhibits a significantly improved statistical quality; a marginally reduced r(2) (0.75 versus 0.83) is more than compensated by a highly improved predictivity with a q(2) of 0.63 versus 0.43. To prove model quality, external prediction of seven test set quinolones was performed. The precise prediction of all test set molecules nicely demonstrates the robustness and statistical significance of the obtained chemometric model using the extended VolSurf approach.