The bone ecosystem facilitates multiple myeloma relapse and the evolution of heterogeneous drug resistant disease.

Multiple myeloma (MM) is an osteolytic malignancy that is incurable due to the emergence of treatment resistant disease. Defining how, when and where myeloma cell intrinsic and extrinsic bone microenvironmental mechanisms cause relapse is challenging with current biological approaches. Here, we report a biology-driven spatiotemporal hybrid agent-based model of the MM-bone microenvironment. Results indicate MM intrinsic mechanisms drive the evolution of treatment resistant disease but that the protective effects of bone microenvironment mediated drug resistance (EMDR) significantly enhances the probability and heterogeneity of resistant clones arising under treatment. Further, the model predicts that targeting of EMDR deepens therapy response by eliminating sensitive clones proximal to stroma and bone, a finding supported by in vivo studies. Altogether, our model allows for the study of MM clonal evolution over time in the bone microenvironment and will be beneficial for optimizing treatment efficacy so as to significantly delay disease relapse.

CK1δ and CK1ε Signaling Sustains Mitochondrial Metabolism and Cell Survival in Multiple Myeloma.

Multiple myeloma remains an incurable malignancy due to acquisition of intrinsic programs that drive therapy resistance. Here we report that casein kinase-1δ (CK1δ) and CK1ε are therapeutic targets in multiple myeloma that are necessary to sustain mitochondrial metabolism. Specifically, the dual CK1δ/CK1ε inhibitor SR-3029 had potent in vivo and ex vivo anti-multiple myeloma activity, including against primary multiple myeloma patient specimens. RNA sequencing (RNA-seq) and metabolic analyses revealed inhibiting CK1δ/CK1ε disables multiple myeloma metabolism by suppressing genes involved in oxidative phosphorylation (OxPhos), reducing citric acid cycle intermediates, and suppressing complexes I and IV of the electron transport chain. Finally, sensitivity of multiple myeloma patient specimens to SR-3029 correlated with elevated expression of mitochondrial genes, and RNA-seq from 687 multiple myeloma patient samples revealed that increased CSNK1D, CSNK1E, and OxPhos genes correlate with disease progression and inferior outcomes. Thus, increases in mitochondrial metabolism are a hallmark of multiple myeloma progression that can be disabled by targeting CK1δ/CK1ε. SIGNIFICANCE: CK1δ and CK1ε are attractive therapeutic targets in multiple myeloma whose expression increases with disease progression and connote poor outcomes, and that are necessary to sustain expression of genes directing OxPhos.

Glutathione levels are associated with methotrexate resistance in acute lymphoblastic leukemia cell lines.

Introduction: Methotrexate (MTX), a folic acid antagonist and nucleotide synthesis inhibitor, is a cornerstone drug used against acute lymphoblastic leukemia (ALL), but its mechanism of action and resistance continues to be unraveled even after decades of clinical use. Methods: To better understand the mechanisms of this drug, we accessed the intracellular metabolic content of 13 ALL cell lines treated with MTX by 1H-NMR, and correlated metabolome data with cell proliferation and gene expression. Further, we validated these findings by inhibiting the cellular antioxidant system of the cells in vitro and in vivo in the presence of MTX. Results: MTX altered the concentration of 31 out of 70 metabolites analyzed, suggesting inhibition of the glycine cleavage system, the pentose phosphate pathway, purine and pyrimidine synthesis, phospholipid metabolism, and bile acid uptake. We found that glutathione (GSH) levels were associated with MTX resistance in both treated and untreated cells, suggesting a new constitutive metabolic-based mechanism of resistance to the drug. Gene expression analyses showed that eight genes involved in GSH metabolism were correlated to GSH concentrations, 2 of which (gamma-glutamyltransferase 1 [GGT1] and thioredoxin reductase 3 [TXNRD3]) were also correlated to MTX resistance. Gene set enrichment analysis (GSEA) confirmed the association between GSH metabolism and MTX resistance. Pharmacological inhibition or stimulation of the main antioxidant systems of the cell, GSH and thioredoxin, confirmed their importance in MTX resistance. Arsenic trioxide (ATO), a thioredoxin inhibitor used against acute promyelocytic leukemia, potentiated MTX cytotoxicity in vitro in some of the ALL cell lines tested. Likewise, the ATO+MTX combination decreased tumor burden and extended the survival of NOD scid gamma (NSG) mice transplanted with patient-derived ALL xenograft, but only in one of four ALLs tested. Conclusion: Altogether, our results show that the cellular antioxidant defense systems contribute to leukemia resistance to MTX, and targeting these pathways, especially the thioredoxin antioxidant system, may be a promising strategy for resensitizing ALL to MTX.

CancerCellTracker: a brightfield time-lapse microscopy framework for cancer drug sensitivity estimation.

MOTIVATION: Time-lapse microscopy is a powerful technique that relies on images of live cells cultured ex vivo that are captured at regular intervals of time to describe and quantify their behavior under certain experimental conditions. This imaging method has great potential in advancing the field of precision oncology by quantifying the response of cancer cells to various therapies and identifying the most efficacious treatment for a given patient. Digital image processing algorithms developed so far require high-resolution images involving very few cells originating from homogeneous cell line populations. We propose a novel framework that tracks cancer cells to capture their behavior and quantify cell viability to inform clinical decisions in a high-throughput manner. RESULTS: The brightfield microscopy images a large number of patient-derived cells in an ex vivo reconstruction of the tumor microenvironment treated with 31 drugs for up to 6 days. We developed a robust and user-friendly pipeline CancerCellTracker that detects cells in co-culture, tracks these cells across time and identifies cell death events using changes in cell attributes. We validated our computational pipeline by comparing the timing of cell death estimates by CancerCellTracker from brightfield images and a fluorescent channel featuring ethidium homodimer. We benchmarked our results using a state-of-the-art algorithm implemented in ImageJ and previously published in the literature. We highlighted CancerCellTracker's efficiency in estimating the percentage of live cells in the presence of bone marrow stromal cells. AVAILABILITY AND IMPLEMENTATION: https://github.com/compbiolabucf/CancerCellTracker. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Plasma cell dependence on histone/protein deacetylase 11 reveals a therapeutic target in multiple myeloma.

The clinical utility of histone/protein deacetylase (HDAC) inhibitors in combinatorial regimens with proteasome inhibitors for patients with relapsed and refractory multiple myeloma (MM) is often limited by excessive toxicity due to HDAC inhibitor promiscuity with multiple HDACs. Therefore, more selective inhibition minimizing off-target toxicity may increase the clinical effectiveness of HDAC inhibitors. We demonstrated that plasma cell development and survival are dependent upon HDAC11, suggesting this enzyme is a promising therapeutic target in MM. Mice lacking HDAC11 exhibited markedly decreased plasma cell numbers. Accordingly, in vitro plasma cell differentiation was arrested in B cells lacking functional HDAC11. Mechanistically, we showed that HDAC11 is involved in the deacetylation of IRF4 at lysine103. Further, targeting HDAC11 led to IRF4 hyperacetylation, resulting in impaired IRF4 nuclear localization and target promoter binding. Importantly, transient HDAC11 knockdown or treatment with elevenostat, an HDAC11-selective inhibitor, induced cell death in MM cell lines. Elevenostat produced similar anti-MM activity in vivo, improving survival among mice inoculated with 5TGM1 MM cells. Elevenostat demonstrated nanomolar ex vivo activity in 34 MM patient specimens and synergistic activity when combined with bortezomib. Collectively, our data indicated that HDAC11 regulates an essential pathway in plasma cell biology establishing its potential as an emerging theraputic vulnerability in MM.

IAP and HDAC inhibitors interact synergistically in myeloma cells through noncanonical NF-κB- and caspase-8-dependent mechanisms.

Interactions between the inhibitor of apoptosis protein antagonist LCL161 and the histone deacetylase inhibitor panobinostat (LBH589) were examined in human multiple myeloma (MM) cells. LCL161 and panobinostat interacted synergistically to induce apoptosis in diverse MM cell lines, including those resistant to bortezomib (PS-R). Similar interactions were observed with other histone deacetylase inhibitors (MS-275) or inhibitors of apoptosis protein antagonists (birinapant). These events were associated with downregulation of the noncanonical (but not the canonical) NF-κB pathway and activation of the extrinsic, caspase-8-related apoptotic cascade. Coexposure of MM cells to LCL161/LBH589 induced TRAF3 upregulation and led to TRAF2 and NIK downregulation, diminished expression of BCL-XL, and induction of γH2A.X. Ectopic expression of TRAF2, NIK, or BCL-XL, or short hairpin RNA TRAF3 knock-down, significantly reduced LCL161/LBH589 lethality, as did ectopic expression of dominant-negative FADD. Stromal/microenvironmental factors failed to diminish LCL161/LBH589-induced cell death. The LCL161/LBH589 regimen significantly increased cell killing in primary CD138+ cells (N = 31) and was particularly effective in diminishing the primitive progenitor cell-enriched CD138-/19+/20+/27+ population (N = 23) but was nontoxic to normal CD34+ cells. Finally, combined LCL161/LBH589 treatment significantly increased survival compared with single-agent treatment in an immunocompetent 5TGM1 murine MM model. Together, these findings argue that LCL161 interacts synergistically with LBH589 in MM cells through a process involving inactivation of the noncanonical NF-κB pathway and activation of the extrinsic apoptotic pathway, upregulation of TRAF3, and downregulation of TRAF2/BCL-XL. Notably, this regimen overcomes various forms of resistance, is active against primary MM cells, and displays significant in vivo activity. This strategy warrants further consideration in MM.

A pharmacodynamic model of clinical synergy in multiple myeloma.

BACKGROUND: Multiagent therapies, due to their ability to delay or overcome resistance, are a hallmark of treatment in multiple myeloma (MM). The growing number of therapeutic options in MM requires high-throughput combination screening tools to better allocate treatment, and facilitate personalized therapy. METHODS: A second-order drug response model was employed to fit patient-specific ex vivo responses of 203 MM patients to single-agent models. A novel pharmacodynamic model, developed to account for two-way combination effects, was tested with 130 two-drug combinations. We have demonstrated that this model is sufficiently parameterized by single-agent and fixed-ratio combination responses, by validating model estimates with ex vivo combination responses for different concentration ratios, using a checkerboard assay. This new model reconciles ex vivo observations from both Loewe and BLISS synergy models, by accounting for the dimension of time, as opposed to focusing on arbitrary time-points or drug effect. Clinical outcomes of patients were simulated by coupling patient-specific drug combination models with pharmacokinetic data. FINDINGS: Combination screening showed 1 in 5 combinations (21.43% by LD50, 18.42% by AUC) were synergistic ex vivo with statistical significance (P < 0.05), but clinical synergy was predicted for only 1 in 10 combinations (8.69%), which was attributed to the role of pharmacokinetics and dosing schedules. INTERPRETATION: The proposed framework can inform clinical decisions from ex vivo observations, thus providing a path toward personalized therapy using combination regimens. FUNDING: This research was funded by the H. Lee Moffitt Cancer Center Physical Sciences in Oncology (PSOC) Grant (1U54CA193489-01A1) and by H. Lee Moffitt Cancer Center's Team Science Grant. This work has been supported in part by the PSOC Pilot Project Award (5U54CA193489-04), the Translational Research Core Facility at the H. Lee Moffitt Cancer Center & Research Institute, an NCI-designated Comprehensive Cancer Center (P30-CA076292), the Pentecost Family Foundation, and Miles for Moffitt Foundation.

BCL2 Amplicon Loss and Transcriptional Remodeling Drives ABT-199 Resistance in B Cell Lymphoma Models.

Drug-tolerant "persister" tumor cells underlie emergence of drug-resistant clones and contribute to relapse and disease progression. Here we report that resistance to the BCL-2 targeting drug ABT-199 in models of mantle cell lymphoma and double-hit lymphoma evolves from outgrowth of persister clones displaying loss of 18q21 amplicons that harbor BCL2. Further, persister status is generated via adaptive super-enhancer remodeling that reprograms transcription and offers opportunities for overcoming ABT-199 resistance. Notably, pharmacoproteomic and pharmacogenomic screens revealed that persisters are vulnerable to inhibition of the transcriptional machinery and especially to inhibition of cyclin-dependent kinase 7 (CDK7), which is essential for the transcriptional reprogramming that drives and sustains ABT-199 resistance. Thus, transcription-targeting agents offer new approaches to disable drug resistance in B-cell lymphomas.

PLK1 stabilizes a MYC-dependent kinase network in aggressive B cell lymphomas.

Concordant activation of MYC and BCL-2 oncoproteins in double-hit lymphoma (DHL) results in aggressive disease that is refractory to treatment. By integrating activity-based proteomic profiling and drug screens, polo-like kinase-1 (PLK1) was identified as an essential regulator of the MYC-dependent kinome in DHL. Notably, PLK1 was expressed at high levels in DHL, correlated with MYC expression, and connoted poor outcome. Further, PLK1 signaling augmented MYC protein stability, and in turn, MYC directly induced PLK1 transcription, establishing a feed-forward MYC-PLK1 circuit in DHL. Finally, inhibition of PLK1 triggered degradation of MYC and of the antiapoptotic protein MCL-1, and PLK1 inhibitors showed synergy with BCL-2 antagonists in blocking DHL cell growth, survival, and tumorigenicity, supporting clinical targeting of PLK1 in DHL.

An Ex Vivo Platform for the Prediction of Clinical Response in Multiple Myeloma.

Multiple myeloma remains treatable but incurable. Despite a growing armamentarium of effective agents, choice of therapy, especially in relapse, still relies almost exclusively on clinical acumen. We have developed a system, Ex vivo Mathematical Myeloma Advisor (EMMA), consisting of patient-specific mathematical models parameterized by an ex vivo assay that reverse engineers the intensity and heterogeneity of chemosensitivity of primary cells from multiple myeloma patients, allowing us to predict clinical response to up to 31 drugs within 5 days after bone marrow biopsy. From a cohort of 52 multiple myeloma patients, EMMA correctly classified 96% as responders/nonresponders and correctly classified 79% according to International Myeloma Working Group stratification of level of response. We also observed a significant correlation between predicted and actual tumor burden measurements (Pearson r = 0.5658, P < 0.0001). Preliminary estimates indicate that, among the patients enrolled in this study, 60% were treated with at least one ineffective agent from their therapy combination regimen, whereas 30% would have responded better if treated with another available drug or combination. Two in silico clinical trials with experimental agents ricolinostat and venetoclax, in a cohort of 19 multiple myeloma patient samples, yielded consistent results with recent phase I/II trials, suggesting that EMMA is a feasible platform for estimating clinical efficacy of drugs and inclusion criteria screening. This unique platform, specifically designed to predict therapeutic response in multiple myeloma patients within a clinically actionable time frame, has shown high predictive accuracy in patients treated with combinations of different classes of drugs. The accuracy, reproducibility, short turnaround time, and high-throughput potential of this platform demonstrate EMMA's promise as a decision support system for therapeutic management of multiple myeloma. Cancer Res; 77(12); 3336-51. ©2017 AACR.