Mitra, Amit

Abstract Multiple myeloma (MM) is the second-most common hematological malignancy in the United States. Myeloma tumors have significant genetic heterogeneity that impacts disease progression and promotes resistance to treatment. RAS is a key oncogenic driver that activates the mitogen-activated protein kinase (MAPK) and the phosphatidylinositol 3-kinase (PI3K)/AKT pathways..Mutations in genes that code for RAS isozymes, KRAS, HRAS, or NRAS, lead to increased levels of activated GTP-RAS and are associated with the most aggressive forms of cancer. Mutations in the RAS family of genes are most common (40-50%) in myeloma. We and others have reported that RAS mutations may not be disease-initiating but rather arise during disease progression and relapse or as subclonal events that impart selective growth advantage and emerging drug resistance. In addition, non-mutant RAS can also be hyperactivated from upstream signaling (like receptor tyrosine kinases) or loss-of-function mutations in GTPase-activating proteins (GAPs). Given the importance of RAS mutations, allele-specific RAS inhibitors have long been considered an attractive strategy for myeloma treatment. However, multiple studies have reported various mechanisms of resistance to allele-specific RAS inhibitors as well as isoform-specific pan-RAS inhibitors, including newly acquired activating RAS mutations, RAS amplification, unchecked activation of co-expressed wild-type RAS (WT-RAS) isozymes, or acquired bypass mechanisms like new oncogenic fusions. Since MM remains incurable, we sought to improve upon currently available treatments by developing small-molecule inhibitors that preferentially inhibit myeloma cells with activated RAS isozymes. We have created NRasQ61X and KRasG12X mutations by CRISPR-Cas directed edits in a single endogenous RAS allele in human myeloma cell lines (HMCLs) that initially had two wild-type alleles. We found that these mutant-NRAS and mutant-KRAS edited cells have increased proliferation compared to cells with WT-RAS, as well as acquired resistance to several commonly used therapeutics. In this study, we used a triangulation method that combines our drug combination candidate prediction algorithm called ‘secDrug’ and the ‘Direct to Drug’ screening platform that evaluated 84 FDA-approved oncology drugs and emerging therapeutics in these RAS mutant myeloma lines. We identified BIRC5 as a potential target for RAS-mutant myeloma. BIRC5 expression is a key mechanism for the evasion of apoptosis. Notably, oncogenic RAS has been shown to promote MYC stabilization via upregulation of BIRC5. Earlier studies suggested that BIRC5 inhibitors are potent in targeting drug-resistant subclones in hematological malignancies. We performed in vitro cytotoxicity studies and cell-based assays using an in-house panel of HMCLs (n>50) representing WT- and mutant-RAS, as well as innate and acquired response/resistance to the standard of care drugs, proteasome inhibitors (PIs) and immunomodulatory drugs (IMiDs), to confirm the potency of novel selective BIRC5 inhibitors, as single agent (IC50 range 0.22-15.30 nM; Mean 4.431±1.041 nM) and in combination with PIs and IMiDs (Chou-Talalay's Combination Index or CI <0.3, indicating strong synergism). Further, we used bulk and single-cell whole-transcriptomics, siRNA-mediated gene knockdown studies, and functional validation to characterize the treatment-induced genes and molecular pathways underpinning the mechanism of action and drug synergy at the subclone level. Finally, we compared the efficacy of these BIRC5 inhibitors against several RAS inhibitors, including cyclorasin B4-27, a bicyclic peptidyl pan-RAS inhibitor, RMC-6236, a multi-selective noncovalent RAS(on) inhibitor, and ADT-007, a novel, highly potent, and selective pan-RAS inhibitor that we have synthesized, and observed comparable effects on cell death and apoptosis. Ex vivo analysis in bone marrow-derived CD138+ cells from myeloma patients using high-throughput Mass Cytometry (CyTOF or Cytometry by Time-Of-Flight) confirmed the dose-dependent reduction of MCL1, Ki-67, CCND1, c-KIT, and MYC. We will present these results that demonstrate the clinical potential of BIRC5 inhibitors as novel candidates in curbing progression and drug resistance in RAS-mutant myeloma.

Mutations in genes that code for RAS isozymes (N- and K-) are the most common (40-50%) somatic variations in multiple myeloma (MM). While ~46% of patients with newly diagnosed MM carry a RAS/Raf mutation, >64% of patients with refractory/relapsed MM carry these mutations. RAS activation can also occur from upstream signaling through receptor tyrosine kinases (RTK) or loss-of-function mutations in GTPase-activating proteins (GAPs), which are prevalent in MM. Although covalent inhibitors targeting KRAS alleles are in various stages of clinical development, multiple mechanisms of resistance have been reported. We discovered and developed a novel, highly potent, and selective pan-RAS inhibitor, ADT-007, capable of suppressing MAPK/AKT signaling by blocking GTP activation of RAS. Here, we show that ADT-007 is effective against MM, particularly those MM cell lines harboring KRAS and NRAS mutations that are representative of relapsed/refractory myeloma (RRMM). For these studies, we compiled a panel of over 50 human myeloma cell lines (HMCLs) with broad biological and genetic heterogeneity that represent innate drug-response/resistance (refractory disease), as well as 10 pairs of parental and clonally-derived proteasome inhibitor (PI)- and immunomodulatory drug (IMD)- resistant HMCLs as models of acquired resistance (relapse). Using CRISPR-cas-directed gene editing, we also have generated clonal NRAS and KRAS mutants of drug-sensitive myeloma lines that initially had two wild-type alleles. All cell lines have been authenticated using the GenePrint 24 System and tested for mycoplasma contamination. First, we performed predictive analysis using single-cell RNA sequencing (scRNAseq; 10X Genomics) on all MM cell lines and demonstrated the potential efficacy of ADT-007 based on the subclonal expression of the target genes/pathways. Using in vitro chemosensitivity assays, ADT-007 was highly potent and effective in killing chemo-resistant and RAS-mutant myeloma cells as a single agent with IC50 values ranging from 0.76 to 12 nM with synergistic benefits if combined with PIs and IMDs. Combination indices (CI) were found to be consistently below 0.9, indicating a synergistic drug interaction. (Drug synergy was determined by CalcuSyn software based on Chou-Talalay's CI method and the isobologram algorithm). Caspase 3/7 assay and flow cytometry (Annexin V-FITC + propidium iodide) confirmed apoptosis induction by ADT-007. Furthermore, the in vitro efficacy of ADT-007 was significantly better (nM vs. µM IC50) than cyclorasin B4-27, a bicyclic peptidyl that has been reported as a potent pan-RAS inhibitor against several other cancers, and RMC-6236, a multi-selective noncovalent RAS(on) inhibitor which is in clinical trial for patients with advanced solid tumors harboring G12X, G13X, and Q61X mutations. Pre-vs-post-treatment bulk-RNA sequencing followed by pathway analysis (using Ingenuity Pathway Analysis and Gene Set Enrichment Analysis) showed that ADT-007 suppressed KRAS and NRAS in a dose-dependent manner. Immunoblotting analysis confirmed the dose-dependent reduction of protein markers implicated in cell survival and chemotherapeutic resistance in myeloma, including c-myc, Ki-67, NFkB, IRF4, IKZF1 (Ikaros) and IKZF3 (Aiolos). Finally, we performed pre-vs-post-treatment single-cell transcriptomics analysis, which revealed a distinct shift in clusters representing myeloma single-cell subpopulations following ADT-007 treatment. High-dimensional mass cytometry (cytometry by time of flight or CyTOF) analysis using heavy metal-conjugated proteins in RRMM patient bone marrow-derived primary myeloma cells is currently underway to confirm the cytotoxic effect of ADT-007 ex vivo. ADT-1004 is an orally bioavailable prodrug of ADT-007 designed to improve water solubility and metabolic stability. ADT-1004 has shown strong antitumor activity in highly aggressive and clinically relevant mouse models of pancreatic cancer at dosages, is well tolerated, and is being developed by the NCI Stepping Stones program. Next, we plan to evaluate the antitumor activity of ADT-1004 in mouse xenograft models of human myeloma using immunocompromised mice. Our research will lay the groundwork for future development efforts needed to advance ADT-1004 to clinical trials involving patients with RRMM.

Abstract Multiple myeloma (MM) is the most common blood cancer in African American (AA) individuals. AAs have an increased risk of developing Monoclonal Gammopathy of Undetermined Significance (MGUS) and MM compared to non-Hispanic White (NHW) individuals. While several studies have shown a pivotal role for the immune tumor microenvironment (iTME) in MM pathogenesis, the composition of the iTME in AAs with MM is under investigated. Viably frozen whole bone marrow (BM) aspirates from 128 individuals (53 AA and 75 NHW) with either MM (n=65), MGUS (n=28), smoldering MM (SMM, n=13), or healthy controls (n=22) were evaluated by single cell RNA-sequencing (scRNA-seq) or mass cytometry (CyTOF). We characterized 34 (14 AA and 20 NHW) by scRNA-seq and 94 (39 AA and 55 NHW) by CyTOF. CyTOF samples were stained with Maxpar Direct Immune Profiling Assay Kit (Fluidigm), acquired on a Helios mass cytometer and analyzed in Cytobank. Data were analyzed by logistic equivalent regression model with p<0.05 for significance. ScRNA-seq samples were processed as in Yao, 2022. Data were normalized using scvi-tools and clustered by Uniform manifold approximation and projection (UMAP). Differentially expressed markers between clusters were identified using Model-based Analysis of Single-cell Transcriptomics (MAST) in Seurat R package at FDR < 0.05 and log2 fold change > 0.55. ScRNA-seq of the CD45+ non-tumor lymphocytes revealed 17 immune populations (HSC, erythroblast, MDSC, GMP, CD1c+ DC, CD4+ and CD8+ T cell, CD14+ M1 macrophage, CD16+monocyte, monocyte, pDC, naïve B, pre-B, pro-B, memory B, NK/MAIT). CD3E/D+cells were clustered into naïve/central memory (CCR7+, LEF1+, SELL+), memory (GZMK+), effector (GZMB+, GZMH+, FGFBP2+), g/d (TRDC+), MAIT (SLC4A10+, KLRB1+) and IFN CD8+ T cells (IFIT1+, IFIT2+,IFIT3+). Disease progression was associated with increases in median CD4+ and CD8+ GZMK+ T cells, memory B cells and a reduction in all monocyte clusters. Within the diseased cohort in comparison to NHW patients, AA patients had increased CD3+ T cells (CD4 and CD8+) (AA 57.5% vs. NHW 44.3%) and a reduction in overall monocytes (AA 10.5% vs. 26.5%). These results were validated by CyTOF, which confirmed an increase in CD3+ T cells (AA 53.7% vs. NHW 35.7%, p<0.001) including increases in CD8+ (AA 7.0% vs. NHW 3.2%, p<0.05) and CD4+ terminal effector T cells (AA 5.7% vs. NHW 2.7%, p<0.05) and a reduction in monocytes (AA 18.3% vs. NHW 26.6%) among AAs. After adjusting for age and sex, AAs retained an increase in T cells (OR 1.59 (95%CI 1.2-2.01, p=0.00023)) and a reduction in monocytes (OR 0.88 (95%CI 0.78-1, p=0.045)). Despite the increase in CD8+ T cells in AA patients, we identified reduced GZMK+ T cells and reductions in exhaustion markers (TIGIT, LAG3, TOX) in CD8+ T cells in comparison to NHWs. Nevertheless, the distribution of the immune populations was similar between healthy AA and NHW individuals. Our multiomics approach to characterize the iTME identified significant differences among AA compared to NHW patients with MM. Citation Format: Suganti Shivaram, Hongwei Tang, Huihuang Yan, Shulan Tian, Guyu Qin, Emilie Anderson, Neeraj Sharma, Abiola Bolarinwa, Cinthya Zepeda-Mendoza, Stacey Lehman, Felicia Gomez, Amit Mitra, Surendra Dasari, Taxiarchis Kourelis, Shaji Kumar, Linda Baughn. Characterization of the Immune Tumor Microenvironment in a Diverse Cohort of Patients with Multiple Myeloma [abstract]. In: Proceedings of the Blood Cancer Discovery Symposium; 2024 Mar 4-6; Boston, MA. Philadelphia (PA): AACR; Blood Cancer Discov 2024;5(2_Suppl):Abstract nr P06.

Multiple myeloma is the second-most common hematological malignancy in the US that remains a challenging disease to cure despite recent advancements in treatment strategies. Drug resistance is a major cause of concern in myeloma chemotherapy. We have created an optimization-regularization-based computational prediction method called secDrug that identified several novel secondary drug (‘secDrug‘) combinations against drug-resistant myeloma (relapsed/refractory multiple myeloma or RRMM). Top among our predicted secDrugs includes the survivin inhibitor, YM155. Survivin is a member of the inhibitor of apoptosis protein (IAP) family that has dual roles in mitosis and apoptosis. YM155 inhibits survivin at the transcriptional level by inhibiting the survivin gene promoter. Several studies have demonstrated the efficacy of YM155 in solid tumors and heme malignancies, including myeloma. However, the precise mechanism of action is not yet known. MM cell lines express high levels of IAPs. We observed high in vitro cytotoxicity of YM155 as a single agent (IC 50 range 0.41 to 12.85nM) as well as in combination with the primary drugs like proteasome inhibitors and Immunomodulatory drugs (IMiDs) against >15 human myeloma cell lines representing innate/refractory and acquired/emerging (relapse) resistance. Next, we explored the molecular mechanisms underlying the efficacy of YM155 using pre-vs-post-treatment tumor whole-transcriptome analysis (RNAseq). The top differentially regulated genes include DDIT3 (CHOP), ATF3, TXNDC5 SEC24D, and DHX15 ( Figure 1A). Ingenuity Pathway Analysis (IPA) predicted that, in addition to regulating apoptosis, survivin inhibition in myeloma significantly upregulated unfolded protein response and downregulated vehicle trafficking and genes associated with immune cell function. Further, we performed single-cell transcriptomics (scRNAseq; 10X Genomics) analysis to explore treatment-induced changes to the subclonal architecture. Our scRNAseq results revealed a distinct shift in clusters representing myeloma single-cell subpopulations following YM155 treatment ( Figure 1B). For example, Cluster 1 (53.11%), characterized by high expression of cell survival markers Ki67 and MYC, was pre-dominant in untreated myeloma cells, while Cluster 4 (52.19%) was the primary subclone in post-treatment tumors. Cluster 3 was shared between the untreated and treated groups. Detailed distinctions between each single-cell cluster will be presented. Furthermore, our scRNAseq data alluded to the benefits of combining survivin-inhibitors and Hsp90-inhibitors based on the expression of drug target genes. Hsp90 inhibitors indirectly target survivin by disrupting the interaction between survivin and Hsp90, resulting in a destabilized survivin protein that ultimately results in survivin degradation. In fact, we have earlier shown the effectiveness of Hsp90 inhibitors (17-AAG, AICAR, and CCT018159) as potent secDrugs against RRMM. Finally, high-dimensional mass cytometry or CyTOF analysis using heavy metal-conjugated proteins in patient bone marrow-derived primary myeloma cells ( ex vivo model system) showed that YM155 induced dose-dependent reduction of several protein markers implicated in cell survival and chemotherapeutic resistance in multiple myeloma, including myc, Ki-67, CD28, CD81, and CD274 (PD-L1). Previous studies have shown that CD28 and CD81 are high-risk immunophenotypes associated with significantly worse prognosis in patients with myeloma. PD-L1 (Programmed death-ligand 1) is a cell-surface glycoprotein expressed in a large number of malignancies that is involved in the termination of immune response and immune evasion. An earlier study has shown that YM155 enhances daratumumab-mediated cellular lysis of multiple myeloma cells. Our RNAseq and CyTOF results thus suggest a potential basis of synergy between YM155 with immunotherapies approved for myeloma, which we will test further. Our approach thus integrates in silico prediction with single-cell multi-ome analysis to identify molecular mechanisms potentially underlying subclonal response to novel combination therapy candidates (secDrugs) for the treatment of RRMM.

Myeloma is well known to have significant genetic heterogeneity among patients' tumors, affecting both progression as well as therapeutic responses. Moreover, subclonal diversity can lead to selective outgrowths that are influenced by new mutations, as well as the bone marrow microenvironment and therapeutic exposures. The ras family of genes (N -and K-) are the most common mutations (40-50% ) in myeloma. Notably, ras mutations occur among a variety of other co-existing chromosomal abnormalities. Ras mutations were found to be more commonly associated with cyclin D1 high expression (e.g., resulting from 11;14 translocations). We and others have also shown that rasmutations may not be disease-initiating but may arise during disease progression and relapse; or are subclonal events that impart selective growth advantage and emerging drug resistance. There is evidence that activating mutations (codons 12,13, 61) in K- ras may have different clinical impact than activating mutations in N- ras yet there are conflicting clinical outcomes that may be due to therapies analyzed or differences in other co-existing genetic mutations. Previous studies using transfection of genes with targeted mutations show random integration and often in multiple copies, thus, limiting interpretations. We chose to create N ras61 and K ras12 mutations by CRISPR-cas directed edits in a single endogenous ras allele in a myeloma cell line (FLAM76) that initially had two wild type alleles (and an 11;14 translocation). We then compared the impact of these edits on growth, signaling, gene expression, and drug responses. We find both the mutN ras and mutK ras edited cells have increased proliferation compared to wild type. However, we find that the N ras mutation increases MAPK signaling and downstream gene activation to a greater extent than the K rasmutation. Based on RNAseq analysis, a MAPK gene activation score was higher in mutN ras edited cells than mutK ras edited cells. Further, we see distinct gene expression patterns between mutN ras versus mutK ras edited cells (See heat map Fig) impacting distinct pathways. Notably, mutNras cells have much higher NOXA expression, increasing sensitivity to apoptosis-inducing drugs (e.g., proteasome inhibitors); and Kras induces anti-apoptotic Bcl-xL, reducing sensitivity. We have conducted a high throughput drug screen with 84 drugs, and while some drugs are active against both mutN ras and mutK ras edits (Proteasome inhibitors), we have identified drugs with differential response to K versus N ras mutations. For example, mutK ras FLAM76 cells are more resistant to proteasome inhibitors; and particularly more resistant to HSP inhibitors and HDAC inhibitors compared to wild type or mutN ras cells. In contrast, mutK ras cells are more sensitive to PikFyve and IL-12 inhibitors than mutNras or wild type cells. Gene expression analysis shows that mutKras cells result in more downregulation of apoptotic genes and higher upregulation of cell survival genes, as well as ER stress pathway genes. Our results suggest that treating myelomas containing K versus N ras activating mutations may benefit from understanding distinct pathway activations and targeting selective drug therapies. We are currently testing new pan-Ras inhibitors to determine if there are distinct responses between our N- and K- ras edited cells, and we are expanding testing to a wider panel of myeloma cell lines containing ras mutations in the context of other genetic abnormalities. Details of gene expression distinctions, pathways, and drug responses that distinguish mutN ras and mutK ras in myeloma will be presented.

Abstract Multiple myeloma (MM) is the second-most common hematological malignancy in the US. MM is an incurable, age-dependent plasma cell neoplasm with a 5-year survival rate of less than 50%. Extensive inter-individual variation in response to standard-of-care drugs like proteasome inhibitors (PIs) and immunomodulatory drugs (IMiDs), drug resistance, and dose-limiting toxicities are critical problems for the treatment of MM. Clinical success in anti-myeloma treatment, therefore, warrants continuous development of novel combination therapy strategies with the explicit goal to improve the therapeutic efficacy by concomitantly targeting multiple signaling pathways. Previously, we have reported the development of an in-house computational pipeline called secDrug that applies greedy algorithm-based set-covering computational optimization method followed by a regularization technique to predict secondary drugs that can be repurposed as novel synergistic partners of standard-of-care drugs for the management of refractory/ resistant MM. Top among these secondary drugs (secDrugs) were the HSP90 inhibitor 17-AAG. In this study, we used 17-AAG as a proof of principle to establish a pipeline that integrates our in silico predictions with in vitro and ex vivo validation as well as multi-omics technologies to identify, validate, and characterize therapeutic agents that could be used either alone or in combination with standard-of-care drugs for the treatment of R/R MM patients (Figure 1). To screen and validate our in silico prediction results, we performed in vitro cytotoxicity assays using 17-AAG on a panel of human myeloma cell lines (HMCLs; in vitro model systems) that captures a wide range of biological and genetic heterogeneity representing the complexities encountered in clinical settings. These cell lines include HMCLs representing innate sensitive/resistance, >10 pairs of parental and clonally-derived PI- and IMiD-resistant pairs (P vs VR or LenR; representing acquired/emerging resistance/relapse), NRAS mutants which leads to the constitutive activation of oncogenic Ras signaling, and CRISPR-edited HSP90 knockdown cell line. Our results showed that 17-AAG has high synergistic activity in combination with PI in inducing apoptosis even in innate and acquired PI-resistant HMCLs and significantly reduces the effective dose of PI required to achieve IC 50 (Chou-Talalay's Dose Reduction Index or DRI 7±1.4). Moreover, 17-AAG+IMiD showed synergistic cell killing activity in clonally-derived IMiD resistant HMCL. Further, 17-AAG induced cell death was comparable with Hsp90 knockdown as evident from the cytotoxicity assay using PI and 17-AAG in combination in RPMI8226-wild type and RPMI-HSP90AA1 knocked down cell line. Notably, 17-AAG was strikingly effective against the NRAS-mutant cell line indicating an additional niche (NRas mutant myeloma) where 17-AAG could be most effective. Next, we performed RNA sequencing to elucidate the molecular mechanisms behind 17-AAG drug action, drug synergy, 17-AAG-induced cell death. Our gene expression profiling (GEP) followed by Ingenuity Pathway Analysis (IPA) analysis revealed protein ubiquitination, aryl hydrocarbon receptor signalling pathway as the top canonical pathways. 17-AAG induced apoptosis via mitochondrial mediated pathway in myeloma. 17-AAG exerts its cytotoxic effect by activating intrinsic pathway of apoptosis which we further confirmed through the increase in reactive oxygen species generation and decrease in mitochondrial membrane potential. 17-AAG was also effective in reducing the expression of hallmarks of MM such as p65/NF-kB, IRF4, c-Myc. Finally, we performed mass cytometry (CyTOF; Cytometry by time of flight) on primary bone-marrow cells (PMCs) from myeloma patients for further validation of proteomic signatures at the single-cell level. CyTOF analysis confirmed 17-AAG-induced cell death and key changes in MM-specific proteomic markers. 17-AAG treated PMCs showed elevated cleaved caspase levels and down-regulation of IRF4 and phospho-STAT3. GEP and CyTOF results were confirmed using immunoblotting assays. Together, our study demonstrates a unique pipeline for drug repositioning that has the potential to revolutionize clinical decision-making by minimizing the number of drugs required for discovering successful combination chemotherapy regimens against drug-resistant myeloma. Figure 1 Figure 1. Disclosures Kumar: BMS: Consultancy, Research Funding; Abbvie: Consultancy, Membership on an entity's Board of Directors or advisory committees, Research Funding; Amgen: Consultancy, Research Funding; Takeda: Consultancy, Membership on an entity's Board of Directors or advisory committees, Research Funding; Tenebio: Research Funding; Beigene: Consultancy; Oncopeptides: Consultancy; Antengene: Consultancy, Honoraria; Carsgen: Research Funding; Janssen: Consultancy, Membership on an entity's Board of Directors or advisory committees, Research Funding; KITE: Consultancy, Membership on an entity's Board of Directors or advisory committees, Research Funding; Merck: Research Funding; Astra-Zeneca: Consultancy, Membership on an entity's Board of Directors or advisory committees, Research Funding; Novartis: Research Funding; Roche-Genentech: Consultancy, Research Funding; Bluebird Bio: Consultancy; Celgene: Membership on an entity's Board of Directors or advisory committees, Research Funding; Adaptive: Membership on an entity's Board of Directors or advisory committees, Research Funding; Sanofi: Research Funding.

Abstract Multiple myeloma (MM) is characterized by significant genetic diversity at subclonal levels that likely plays a defining role in the heterogeneity of tumor progression, clinical aggressiveness and drug sensitivity. Such heterogeneity is a driving factor in the evolution of MM, from founder clones through outgrowth of subclonal fractions. DNA Sequencing studies on MM samples have indeed demonstrated such heterogeneity in subclonal architecture at diagnosis based on recurrent mutations in pathologically relevant genes that may ultimately to lead to relapse. However, no study so far has reported a predictive gene expression signature that can identify, distinguish and quantify drug sensitive and drug-resistant subpopulations within a bulk population of myeloma cells. In recent years, our laboratory has successfully developed a gene expression profile (GEP)-based signature that could not only distinguish drug response of MM cell lines, but also was effective in stratifying patient outcomes when applied to GEP profiles from MM clinical trials using proteasome inhibitors (PI) as chemotherapeutic agents. Further, we noted myeloma cell lines that responded to the drug often contained residual sub-population of cells that did not respond, and likely were selectively propagated during drug treatment in vitro, and in patients. In this study, we performed targeted qRT-PCR analysis of single cells using a gene panel that included PI sensitivity genes and gene signatures that could discriminate between low and high-risk myeloma followed by intensive bioinformatics and statistical analysis for the classification and prediction of PI response in individual cells within bulk multiple myeloma tumors. Fluidigm's C1 Single-Cell Auto Prep System was used to perform automated single-cell capture, processing and cDNA synthesis on 576 pre-treatment cells from 12 cell lines representing a wide range of PI-sensitivity and 370 cells from 7 patient samples undergoing PI treatment followed by targeted gene expression profiling of single cells using automated, high-throughput on-chip qRT-PCR analysis using 96.96 Dynamic Array IFCs on the BioMark HD System. Probability of resistance for each individual cell was predicted using a pipeline that employed the machine learning methods Random Forest, Support Vector Machine (radial and sigmoidal), LASSO and kNN (k Nearest Neighbor) for making single-cell GEP data-driven predictions/ decisions. The weighted probabilities from each of the algorithms were used to quantify resistance of each individual cell and plotted using Ensemble forecasting algorithm. Using our drug response GEP signature at the single cell level, we could successfully identify distinct subpopulations of tumor cells that were predicted to be sensitive or resistant to PIs. Subsequently, we developed a R Statistical analysis package (http://cran.r-project.org), SCATTome (Single Cell Analysis of Targeted Transcriptome), that can restructure data obtained from Fluidigm qPCR analysis run, filter missing data, perform scaling of filtered data, build classification models and successfully predict drug response of individual cells and classify each cell's probability of response based on the targeted transcriptome. We will present the program output as graphical displays of single cell response probabilities. This package provides a novel classification method that has the potential to predict subclonal response to a variety of therapeutic agents. Disclosures Kumar: Skyline: Consultancy, Honoraria; BMS: Consultancy; Onyx: Consultancy, Research Funding; Sanofi: Consultancy, Research Funding; Janssen: Consultancy, Research Funding; Novartis: Research Funding; Takeda: Consultancy, Research Funding; Celgene: Consultancy, Research Funding.