RSK inhibitor BI-D1870 inhibits acute myeloid leukemia cell proliferation by targeting mitotic exit

Discovery of potential RSK1 inhibitors for cancer therapy using virtual screening, molecular docking, molecular dynamics simulation, and MM/GBSA calculations

The p90 ribosomal protein S6 Kinase (RSK) family belongs to Ser/Thr protein kinases that includes four isoform RSK1-4 in mammals. The ribosomal protein S6 Kinase 1 (RSK1) is also known as ribosomal protein S6 kinase alpha-1 (RPS6KA1) is a special protein due to their two catalytic regions that is associated with abundantly various cancers and it is proposed as a drug target. Several RSK1 isoform inhibitors have been reported but none of them are used in clinical studies. Thus, we aimed to perform ligand pharmacophore mapping with the known inhibitor and structure-based virtual screening studies to determine potential candidates against RSK1-terminal kinase domains CTKD and NTKD. The studied compounds from the databases (ApexBio, ChEMBL, ChemDiv). The molecular docking study was performed with the resulted candidates by using CDOCKER and Glide/SP methods. The four candidates with the highest docking scores were used for further 100-ns molecular dynamics (MD) simulations and Molecular Mechanics Generalised Born and Surface Area (MM/GBSA) calculations. The root mean square deviation (RMSD) for protein complexes were found between 2 Å and 4 Å. Solvent accessible surface area (SASA), radius of gyration (Rg), and polar surface area (PSA) values were calculated for compounds. The binding free energies were calculated between -72.22 kcal/mol and -82.44 kcal/mol. The interaction diagrams showed that hydrogen bond, alkyl, and π-alkyl interactions were observed with specific residues such as Leu144, Lys94, Asp142 for RSK1-NTKD, and Cys532, Cys556, Lys447, Asn540 for RSK1-CTKD. The identified compounds may be potential inhibitor candidates of RSK1 following the preclinical studies.Communicated by Ramaswamy H. Sarma.

RSK1 is an exploitable dependency in myeloproliferative neoplasms and secondary acute myeloid leukemia

Myeloid malignancies are heterogenous disorders characterized by distinct molecular drivers but share convergence of oncogenic signaling pathways and propagation by ripe pro-inflammatory niches. Here, we establish a comprehensive transcriptional atlas across the spectrum of myeloproliferative neoplasms (MPN) and secondary acute myeloid leukemia (sAML) through RNA-sequencing of 158 primary samples encompassing CD34+ hematopoietic stem/progenitor cells and CD14+ monocytes. Supported by mass cytometry (CyTOF) profiling, we reveal aberrant networks of PI3K/AKT/mTOR signalling and NFκB-mediated hyper-inflammation. Combining ATAC-Seq, CUT&Tag, RNA-seq, and CyTOF, we demonstrate that targeting of ribosomal protein S6 kinase A1 (RSK1) suppresses NFκB activation and diminishes pro-inflammatory mediators including tumor necrosis factor (TNF) associated with MPN disease severity and transformation. We further evaluate a therapeutic approach utilizing a first-in-class RSK inhibitor, PMD-026, currently in Phase 2 development for breast cancer, for use in myeloid malignancies. Treatment with PMD-026 suppressed disease burden across seven syngeneic and patient-derived xenograft leukemia mouse models spanning the spectrum of driver and disease-modifying mutations. These findings uncover a therapeutic avenue for a conserved dependency across MPN and sAML.

RSK1 dependency in FLT3-ITD acute myeloid leukemia

Internal tandem duplications (ITD) in fms-like tyrosine kinase 3 (FLT3) represent the most common genetic alteration in de novo acute myeloid leukemia (AML). Here, we identify ribosomal protein s6 kinase a1 (RSK1) as a core dependency in FLT3-ITD AML and unveil the existence of crucial bi-directional regulation. RSK1 perturbation resulted in marked apoptosis and abrogated phosphorylation of FLT3 and associated downstream signaling cascades in FLT3-ITD AML cell lines. Using cycloheximide, MG-132, and ubiquitination assays, we further demonstrate mechanistically that RSK1 regulates FLT3-ITD activity, and protein stability through deubiqutinase USP1, which we identify as a second dependency. Importantly, multivariate analysis revealed heightened expression of RPS6KA1 and USP1 to be associated with poor patient prognosis, and these effectors may serve as biomarkers predictive of patient survival and therapeutic response to FLT3-ITD inhibitors. Lastly, RSK1 inhibition utilizing a first-in-class RSK inhibitor, PMD-026, that is currently undergoing Phase 2 development for breast cancer, diminished leukemic disease burden in MV4-11 xenograft and syngeneic Flt3ITDTet2KO leukemia models. These findings illustrate an unconventional and promising therapeutic strategy targeting FLT3-ITD leukemia.

p90RSK pathway inhibition synergizes with cisplatin in TMEM16A overexpressing head and neck cancer

Head and neck squamous cell carcinoma (HNSCC) constitutes one of the most common types of human cancers and often metastasizes to lymph nodes. Platinum-based chemotherapeutic drugs are commonly used for treatment of a wide range of cancers, including HNSCC. Its mode of action relies on its ability to impede DNA repair mechanisms, inducing apoptosis in cancer cells. However, due to acquired resistance and toxic side-effects, researchers have been focusing on developing novel combinational therapeutic strategies to overcome cisplatin resistance. In the current study, we identified p90RSK, an ERK1/2 downstream target, as a key mediator and a targetable signaling node against cisplatin resistance. Our results strongly support the role of p90RSK in cisplatin resistance and identify the combination of p90RSK inhibitor, BI-D1870, with cisplatin as a novel therapeutic strategy to overcome cisplatin resistance. In addition, we have identified TMEM16A expression as a potential upstream regulator of p90RSK through the ERK pathway and a biomarker of response to p90RSK targeted therapy in the context of cisplatin resistance.

Therapeutic targeting of p90 ribosomal S6 kinase

The Serine/Threonine protein kinase family, p90 ribosomal S6 kinases (RSK) are downstream effectors of extracellular signal regulated kinase 1/2 (ERK1/2) and are activated in response to tyrosine kinase receptor or G-protein coupled receptor signaling. RSK contains two distinct kinase domains, an N-terminal kinase (NTKD) and a C-terminal kinase (CTKD). The sole function of the CTKD is to aid in the activation of the NTKD, which is responsible for substrate phosphorylation. RSK regulates various homeostatic processes including those involved in transcription, translation and ribosome biogenesis, proliferation and survival, cytoskeleton, nutrient sensing, excitation and inflammation. RSK also acts as a major negative regulator of ERK1/2 signaling. RSK is associated with numerous cancers and has been primarily studied in the context of transformation and metastasis. The development of specific RSK inhibitors as cancer therapeutics has lagged behind that of other members of the mitogen-activated protein kinase signaling pathway. Importantly, a pan-RSK inhibitor, PMD-026, is currently in phase I/1b clinical trials for metastatic breast cancer. However, there are four members of the RSK family, which have overlapping and distinct functions that can vary in a tissue specific manner. Thus, a problem for transitioning a RSK inhibitor to the clinic may be the necessity to develop isoform specific inhibitors, which will be challenging as the NTKDs are very similar to each other. CTKD inhibitors have limited use as therapeutics as they are not able to inhibit the activity of the NTKD but could be used in the development of proteolysis-targeting chimeras.

Identification and study of cuproptosis-related genes in prognostic model of multiple myeloma

BACKGROUND

Multiple myeloma (MM) is a highly heterogeneous disease. Cuproptosis is a novel mode of death that is closely associated with several diseases, such as hepatocellular carcinoma. However, its role in MM is unknown.

METHODS

MM transcriptomic and clinical data were obtained from UCSC Xena and gene expression omnibus (GEO) databases. Following MM samples were divided into different subtypes based on the cuproptosis genes, the differentially expressed genes (DEGs) among different subtypes, namely, candidate cuproptosis related genes were analyzed by univariate Cox and least absolute shrinkage and selection operator (LASSO) regression to construct a cuproptosis-related risk model. After the independent prognostic analysis was performed, a nomogram was constructed. Finally, Functional enrichment analysis and immune infiltration analysis were performed in the high- and low-risk groups, potential therapeutic agents were then predicted.

RESULTS

The 784 MM samples in UCSC Xena cohorts were divided into three different subtypes, and 4 out of 346 candidate cuproptosis related genes, namely CDKN2A, BCL3, KCNA3 and TTC14 were used to construct a risk model. Risk score was considered a reliable independent prognostic factor for MM patients. It was investigated that the pathway of cell cycle was significantly enriched in the high-risk group. In addition, immune score, ESTIMATE score and cytolytic activity were significantly different between different risk groups, as well as 13 immune cells such as memory B cells. Nine drugs were predicted in our study.

CONCLUSION

A cuproptosis-related prognostic model was constructed, which may have a potential guiding role in the treatment of MM.

The ribosomal protein S6 kinase alpha-1 (RPS6KA1) induces resistance to venetoclax/azacitidine in acute myeloid leukemia

Venetoclax/azacitidine combination therapy is effective in acute myeloid leukemia (AML) and tolerable for older, multimorbid patients. Despite promising response rates, many patients do not achieve sustained remission or are upfront refractory. Identification of resistance mechanisms and additional therapeutic targets represent unmet clinical needs. By using a genome-wide CRISPR/Cas9 library screen targeting 18,053 protein- coding genes in a human AML cell line, various genes conferring resistance to combined venetoclax/azacitidine treatment were identified. The ribosomal protein S6 kinase A1 (RPS6KA1) was among the most significantly depleted sgRNA-genes in venetoclax/azacitidine- treated AML cells. Addition of the RPS6KA1 inhibitor BI-D1870 to venetoclax/azacitidine decreased proliferation and colony forming potential compared to venetoclax/azacitidine alone. Furthermore, BI-D1870 was able to completely restore the sensitivity of OCI-AML2 cells with acquired resistance to venetoclax/azacitidine. Analysis of cell surface markers revealed that RPS6KA1 inhibition efficiently targeted monocytic blast subclones as a potential source of relapse upon venetoclax/azacitidine treatment. Taken together, our results suggest RPS6KA1 as mediator of resistance towards venetoclax/azacitidine and additional RPS6KA1 inhibition as strategy to prevent or overcome resistance.

P90 ribosomal S6 kinases: A bona fide target for novel targeted anticancer therapies?

The 90 kDa ribosomal S6 kinase (RSK) family of proteins is a group of highly conserved Ser/Thr kinases. They are downstream effectors of the Ras/ERK/MAPK signaling cascade. ERK1/2 activation directly results in the phosphorylation of RSKs, which further, through interaction with a variety of different downstream substrates, activate various signaling events. In this context, they have been shown to mediate diverse cellular processes like cell survival, growth, proliferation, EMT, invasion, and metastasis. Interestingly, increased expression of RSKs has also been demonstrated in various cancers, such as breast, prostate, and lung cancer. This review aims to present the most recent advances in the field of RSK signaling that have occurred, such as biological insights, function, and mechanisms associated with carcinogenesis. We additionally present and discuss the recent advances but also the limitations in the development of pharmacological inhibitors of RSKs, in the context of the use of these kinases as putative, more efficient targets for novel anticancer therapeutic approaches.

BI-D1870 Induces Mitotic Dysfunction and Apoptosis in Neuroblastoma by Regulating the PI3K-Akt-mTORC1 Signal Axis

Introduction: Neuroblastoma (NB) is one of the most common extracranial solid malignant tumors in children. The 5-year survival rate of high-risk or refractory NB is less than 50%. Therefore, developing new effective therapeutics for NB remains an urgent challenge. Materials and Methods: Based on the NB dataset TARGET-NBL in the TCGA database, the prognosis-related genes were analyzed using univariate cox regression (p < 0.01). The protein network interaction of prognostic genes was analyzed using STRING to obtain 150 hub genes with HR > 1 and 150 hub genes with HR < 1. The Connectivity Map database was used to predict a therapeutic drug: BI-D1870, a ribosomal S6 kinase inhibitor. The inhibitory effect of BI-D1870 on NB was investigated through in vivo and in vitro experiments, and its inhibitory mechanism was explored. Results: Both the in vivo and in vitro experiments showed that BI-D1870 could inhibit tumor proliferation and induce tumor apoptosis. Furthermore, we proved that BI-D1870 caused G2/M phase arrest and mitosis damage in cells. RNA-seq of cells showed that BI-D1870 may inhibit the growth of NB by inhibiting the PI3K-Akt-mTOR axis. Western blot and immunofluorescence testing showed that BI-D1870 inhibited the PI3K-Akt-mTORC1 signal pathway to regulate the phosphorylation of RPS6 and 4E BP1 proteins, inhibit protein translation, and inhibit microtubule formation, thus preventing mitotic proliferation and inducing apoptosis. Conclusions: This study provides strong support that BI-D1870 may be a potential adjuvant therapy for NB.

RSK inhibitors as potential anticancer agents: Discovery, optimization, and challenges

Ribosomal S6 kinase (RSK) family is a group of serine/threonine kinases, including four isoforms (RSK1/2/3/4). As a downstream effector of the Ras-mitogen-activated protein kinase (Ras-MAPK) pathway, RSK participates in many physiological activities such as cell growth, proliferation, and migration, and is intimately involved in tumor occurrence and development. As a result, it is recognized as a potential target for anti-cancer and anti-resistance therapies. There have been several RSK inhibitors discovered or designed in recent decades, but only two have entered clinical trials. Low specificity, low selectivity, and poor pharmacokinetic properties in vivo limit their clinical translation. Published studies performed structure optimization by increasing interaction with RSK, avoiding hydrolysis of pharmacophores, eliminating chirality, adapting to binding site shape, and becoming prodrugs. Besides enhancing efficacy, the focus of further design will move towards selectivity since there are functional differences among RSK isoforms. This review summarized the types of cancers associated with RSK, along with the structural characteristics and optimization process of the reported RSK inhibitors. Furthermore, we addressed the importance of RSK inhibitors' selectivity and discussed future drug development directions. This review is expected to shed light on the emergence of RSK inhibitors with high potency, specificity, and selectivity.

CDC20: a novel therapeutic target in cancer

Cell division cycle protein 20 (Cdc20) is a member of the cell cyclin family. In the early stage of mitosis, it activates the anaphase-promoting complex (APC) and forms the E3 ubiquitin ligase complex APC^Cdc20, which destroys key regulators of the cell cycle and promotes mitosis. Cdc20 serves as a target for the spindle checkpoint, ensuring proper chromosome segregation. As an oncoprotein, Cdc20 is highly expressed in a variety of malignant tumors, and Cdc20 overexpression is associated with poor prognosis of these tumors. This review aims to dissect the tumorigenic role of Cdc20 in human malignancies and its targeting strategies.

DUSP6 mediates resistance to JAK2 inhibition and drives leukemic progression

Myeloproliferative neoplasms (MPNs) exhibit a propensity for transformation to secondary acute myeloid leukemia (sAML), for which the underlying mechanisms remain poorly understood, resulting in limited treatment options and dismal clinical outcomes. Here, we performed single-cell RNA sequencing on serial MPN and sAML patient stem and progenitor cells, identifying aberrantly increased expression of DUSP6 underlying disease transformation. Pharmacologic dual-specificity phosphatase (DUSP)6 targeting led to inhibition of S6 and Janus kinase (JAK)-signal transducer and activator of transcription (STAT) signaling while also reducing inflammatory cytokine production. DUSP6 perturbation further inhibited ribosomal S6 kinase (RSK)1, which we identified as a second indispensable candidate associated with poor clinical outcome. Ectopic expression of DUSP6 mediated JAK2-inhibitor resistance and exacerbated disease severity in patient-derived xenograft (PDX) models. Contrastingly, DUSP6 inhibition potently suppressed disease development across Jak2V617F and MPLW515L MPN mouse models and sAML PDXs without inducing toxicity in healthy controls. These findings underscore DUSP6 in driving disease transformation and highlight the DUSP6-RSK1 axis as a vulnerable, druggable pathway in myeloid malignancies.

Phospho-proteomics reveals that RSK signaling is required for proliferation of natural killer cells stimulated with IL-2 or IL-15

Natural killer (NK) cells are cytotoxic lymphocytes that play a critical role in the innate immune system. Although cytokine signaling is crucial for the development, expansion, and cytotoxicity of NK cells, the signaling pathways stimulated by cytokines are not well understood. Here, we sought to compare the early signaling dynamics induced by the cytokines interleukin (IL)-2 and IL-15 using liquid chromatography-mass spectrometry (LC-MS)-based phospho-proteomics. Following stimulation of the immortalized NK cell line NK-92 with IL-2 or IL-15 for 5, 10, 15, or 30 min, we identified 8,692 phospho-peptides from 3,023 proteins. Comparing the kinetic profiles of 3,619 fully quantified phospho-peptides, we found that IL-2 and IL-15 induced highly similar signaling in NK-92 cells. Among the IL-2/IL-15-regulated phospho-peptides were both well-known signaling events like the JAK/STAT pathway and novel signaling events with potential functional significance including LCP1 pSer5, STMN1 pSer25, CHEK1 pSer286, STIM1 pSer608, and VDAC1 pSer104. Using bioinformatic approaches, we sought to identify kinases regulated by IL-2/IL-15 stimulation and found that the p90 ribosomal S6 kinase (p90RSK) family was activated by both cytokines. Using pharmacological inhibitors, we then discovered that RSK signaling is required for IL-2 and IL-15-induced proliferation in NK-92 cells. Taken together, our analysis represents the first phospho-proteomic characterization of cytokine signaling in NK cells and increases our understanding of how cytokine signaling regulates NK cell function.

RSK Inhibition Induces Apoptosis by Downregulating Protein Synthesis in a Variety of Acute Myeloid Leukemia Cell Lines

Fms-like tyrosine kinase 3 (FLT3) and isocitrate dehydrogenase 1/2 (IDH1/2) mutations drive malignancy in acute myeloid leukemia (AML), which accounts for approximately 40% of AML cases. Treatment with FLT3 or IDH1/2 inhibitors is used for such patients; however, it is not considered for most patients with AML who lack mutations on the respective genes. In this study, p90 ribosomal S6 kinase (RSK) was found to serve as a new therapeutic target in various AMLs with or without FLT3 mutations. BI-D1870, a potent inhibitor of RSK, significantly suppressed the proliferation of AML cell lines, among which three encoded wild-type FLT3 and three contained FLT3 driver mutations, compared with chronic myeloid leukemia K562 cells or other adherent cancer cells. BI-D1870 inhibited protein synthesis by dephosphorylating the p70 S6 kinase and eukaryotic initiation factor 4E-binding protein 1 in all AML cells except KG-1a cells. Meanwhile, the expression of microtubule-associated protein light chain 3B-I and -II increased in KG-1a cells treated with BI-D1870. BI-D1870 induced caspase-dependent apoptosis in all AML cells, including KG-1a cells. We next investigated the synergistic effect of BI-D1870 with cytarabine, a traditional anticancer drug used in AML. Synergistic effects of BI-D1870 and cytarabine were not observed in any of the cell lines. The findings suggested that BI-D1870 alone exerts an adequate antiproliferative effect on AML with or without FLT3 mutations and serves as a novel AML therapeutic agent.

RSK Isoforms in Acute Myeloid Leukemia

Ribosomal S6 Kinases (RSKs) are a group of serine/threonine kinases that function downstream of the Ras/Raf/MEK/ERK signaling pathway. Four RSK isoforms are directly activated by ERK1/2 in response to extracellular stimuli including growth factors, hormones, and chemokines. RSKs phosphorylate many cytosolic and nuclear targets resulting in the regulation of diverse cellular processes such as cell proliferation, survival, and motility. In hematological malignancies such as acute myeloid leukemia (AML), RSK isoforms are highly expressed and aberrantly activated resulting in poor outcomes and resistance to chemotherapy. Therefore, understanding RSK function in leukemia could lead to promising therapeutic strategies. This review summarizes the current information on human RSK isoforms and discusses their potential roles in the pathogenesis of AML and mechanism of pharmacological inhibitors.

The landscape of gene co-expression modules correlating with prognostic genetic abnormalities in AML

Background

The heterogenous cytogenetic and molecular variations were harbored by AML patients, some of which are related with AML pathogenesis and clinical outcomes. We aimed to uncover the intrinsic expression profiles correlating with prognostic genetic abnormalities by WGCNA.

Methods

We downloaded the clinical and expression dataset from BeatAML, TCGA and GEO database. Using R (version 4.0.2) and ‘WGCNA’ package, the co-expression modules correlating with the ELN2017 prognostic markers were identified (R² ≥ 0.4, p < 0.01). ORA detected the enriched pathways for the key co-expression modules. The patients in TCGA cohort were randomly assigned into the training set (50%) and testing set (50%). The LASSO penalized regression analysis was employed to build the prediction model, fitting OS to the expression level of hub genes by ‘glmnet’ package. Then the testing and 2 independent validation sets (GSE12417 and GSE37642) were used to validate the diagnostic utility and accuracy of the model.

Results

A total of 37 gene co-expression modules and 973 hub genes were identified for the BeatAML cohort. We found that 3 modules were significantly correlated with genetic markers (the ‘lightyellow’ module for NPM1 mutation, the ‘saddlebrown’ module for RUNX1 mutation, the ‘lightgreen’ module for TP53 mutation). ORA revealed that the ‘lightyellow’ module was mainly enriched in DNA-binding transcription factor activity and activation of HOX genes. The ‘saddlebrown’ module was enriched in immune response process. And the ‘lightgreen’ module was predominantly enriched in mitosis cell cycle process. The LASSO- regression analysis identified 6 genes (NFKB2, NEK9, HOXA7, APRC5L, FAM30A and LOC105371592) with non-zero coefficients. The risk score generated from the 6-gene model, was associated with ELN2017 risk stratification, relapsed disease, and prior MDS history. The 5-year AUC for the model was 0.822 and 0.824 in the training and testing sets, respectively. Moreover, the diagnostic utility of the model was robust when it was employed in 2 validation sets (5-year AUC 0.743–0.79).

Conclusions

We established the co-expression network signature correlated with the ELN2017 recommended prognostic genetic abnormalities in AML. The 6-gene prediction model for AML survival was developed and validated by multiple datasets.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-021-02914-2.

N-Substituted pyrrolopyrimidines and purines as p90 ribosomal S6 protein kinase-2 (RSK2) inhibitors

The RSK2 kinase is the downstream effector of the Ras/Raf/MEK/ERK pathway, that is often aberrantly activated in acute myeloid leukemia (AML). Recently, we reported a structure-activity study for BI-D1870, the pan-RSK inhibitor, and identified pteridinones that inhibited cellular RSK2 activity that did not result in concomitant cytotoxicity. In the current study, we developed a series of pyrrolopyrimidines and purines to replace the pteridinone ring of BI-D1870, with a range of N-substituents that extend to the substrate binding site to probe complementary interactions, while retaining the 2,6-difluorophenol-4-amino group to maintain interactions with the hinge domain and the DFG motif. Several compounds inhibited cellular RSK2 activity, and we identified compounds that uncoupled cellular RSK2 inhibition from potent cytotoxicity in the MOLM-13 AML cell line. These N-substituted probes have revealed an opportunity to further examine substituents that extend from the ATP- to the substrate-binding site may confer improved RSK potency and selectivity.