- Primary research
- Open Access
Harmine, a dual-specificity tyrosine phosphorylation-regulated kinase (DYRK) inhibitor induces caspase-mediated apoptosis in neuroblastoma
Cancer Cell International volume 18, Article number: 82 (2018)
Neuroblastoma (NB) is an early childhood malignancy that arises from the developing sympathetic nervous system. Harmine is a tricyclic β-carboline alkaloid isolated from the harmal plant that exhibits both cytostatic and cytotoxic effects. Harmine is capable of blocking the activities of dual-specificity tyrosine phosphorylation-regulated kinase (DYRK) family proteins and mitogen-activated protein kinase. These kinases promote proliferation and inhibit apoptosis.
Four human NB cell lines were used to study the effects of harmine treatment: SKNBE and KELLY (MYCN-amplified) as well as SKNAS and SKNFI (MYCN non-amplified). The anti-cancer properties of harmine were examined by RealTime-Glo MT cell viability assays, caspase activity assays, PARP cleavage using Western blot analysis, and flow cytometry-based Annexin V detection. A molecular interaction model of harmine bound to the DYRK2 family kinase was generated by computational docking using X-ray structures. NB tumors from human patients were profiled for DYRK mRNA expression patterns and clinical correlations using the R2 platform.
The IC50 values for harmine after 72 h treatment were 169.6, 170.8, and 791.7 μM for SKNBE, KELLY, and SKNFI, respectively. Exposure of these NB cell lines to 100 μM of harmine resulted in caspase-3/7 and caspase-9 activation as well as caspase-mediated PARP cleavage and Annexin V-positive stained cells, as early as 24 h after treatment, clearly suggesting apoptosis induction, especially in MYCN-amplified cell lines. Elevated DYRK2 mRNA levels correlated with poor prognosis in a large cohort of NB tumors.
Harmine is a known inhibitor of DYRK family kinases. It can induce apoptosis in NB cell lines, which led us to investigate the clinical correlations of DYRK family gene expression in NB tumors. The patient results support our hypothesis that DYRK inhibition by harmine and the subsequent triggering of caspase-mediated apoptosis might present a novel approach to NB therapy.
Neuroblastoma (NB) is an early childhood malignancy that arises from the developing sympathetic nervous system, resulting in aggressive tumor formation in the sympathetic ganglia and/or the adrenal glands . Amplification of the MYCN gene, leading to over-expression of the MYCN protein, is the most prevalent NB genetic aberration. It is found in ~ 20% of NB, predominantly in high-stage tumors and has been linked to high risk disease, and poor patient prognosis [1, 2]. Of those high risk patients that respond initially to chemotherapy, the majority will succumb to the disease after a relapse into a chemotherapy-resistant state . In addition to MYCN gene amplifications, mutations in genes encoding the mitogen-activated protein kinase (MAPK) pathway, and in the ALK gene have been identified as main drivers in the majority of NB .
The harmal plant (Peganum harmala L., family Zygophyllaceae), also called Syrian rue, is a perennial shrub native to the eastern Mediterranean region. Various parts of the plant have long been used in traditional folk medicine . Harmine is a tricyclic β-carboline alkaloid isolated from harmal seeds and acts as a monoamine oxidase A (MAO-A) inhibitor [3,4,5]. Intake of alkaloids from the harmal plant can have anti-depressive, analgesic, and anti-bacterial pharmacological effects [6,7,8].
Harmine has been shown to induce apoptosis and inhibit cell proliferation, migration, and invasion in a dose-dependent manner in various human cancer cell lines including C334, CCD18LU, HeLa, HL-60, K562, SW480, BGC-823, and SGC-790 . While lower concentrations of harmine typically induced cytostasis, higher concentrations are associated with cytotoxicity .
The dual-specificity tyrosine phosphorylation-regulated kinase (DYRK) family proteins are related to the MAPK family. However, the activating tyrosine phosphorylation of DYRK family kinases is not catalyzed by upstream kinases but occurs through autophosphorylation . There is growing interest in the role of DYRK family kinases in cancer, as they can act as regulators of protein stability during the cell cycle and regulate the activity of the proteasome [11, 12]. Remarkably, harmine inhibits all DYRK family members (DYRK1A, DYRK1B, DYRK2, and DYRK4), with the highest affinity for DYRK1A [13, 14].
In this study, we identified that harmine induces apoptotic cell death in NB cells, generated a molecular interaction model for harmine bound to DYRK2, and showed that DYRK2 mRNA expression patterns in a large cohort of human NB tumors suggest the involvement of DYRK2 in NB tumorigenesis. Together, our results offer a potential new route of NB therapy.
Harmine (Fig. 1), LDN-192960, and INDY were purchased from Cayman Chemical. The compound solids were stored at − 20 °C. Stock solutions were prepared by dissolving harmine (100 mM), LDN-192960 (40 mM) and INDY (40 mM) into sterile DMSO (VWR). Stock concentrations were filtered prior to being added to the cell cultures.
The human NB cell lines KELLY (also known as N206; #92110411 from Sigma), SK-N-AS (called SKNAS in this study; CRL-213 from ATCC), SK-N-BE (clone SKNBE(2)C, called SKNBE in this study; CRL-2268 from ATCC), and SK-N-FI (called SKNFI in this study; from the Children’s Oncology Group) were cultured in RPMI 1640 (VWR), supplemented with 10% heat-inactivated fetal bovine serum (Invitrogen), Penicillin (100 IU/mL) and Streptomycin (100 μg/mL) (30-002-CI, Corning). All cells were purchased from their respective suppliers within the last 2 years. The cells were maintained at 37 °C in a humidified atmosphere containing 5% CO2.
Cell viability assay
Cell viability and IC50 was determined using the RealTime-Glo MT cell viability assay (G9712, Promega). This reagent allows continuous measurement of cell viability in the same well. Cells were plated at a density of 4000 cells/well into white-walled, opaque assay plates. After the plated cells had been given 24 h to adhere, they were treated with harmine concentrations ranging from 0 to 1 mM. The MT Cell Viability Substrate and NanoLuc Enzyme were equilibrated to 37 °C, 2× RealTime-Glo reagent was prepared, and an equal volume was added to each well. For time zero measurements, cells were incubated with reagent for 20 min at 37 °C, and luminescence was measured on a Biotek Synergy microplate reader. Luminescence was measured at 24, 48, and 72 h after addition of harmine.
The colorimetric SRB assay was used to measure cytotoxicity as previously described, following the treatment with DYRK2 inhibitor LDN-192960 or DYRK1A/B inhibitor INDY [15,16,17]. Briefly, NB cells were plated in transparent flat 96-well plates and allowed to attach overnight. At the initiation of each experiment (t = 0) and after drug treatments, cells were fixed with 10% TCA at 4 °C for 1 h, washed with deionized water, and dried at room temperature. Cells were then stained with 100 μl of 0.4% SRB in 1% acetic acid for 20 min at room temperature, rinsed five times with 1% acetic acid and allowed to dry at room temperature. One hundred µl of 10 mM Tris–HCl pH 7.0 was added to each well, shaken for 10 min at room temperature and read at 540 nm using a Molecular Devices Flexstation 3 microplate reader.
Caspase activity assays
The quantification of the relative caspase activity in harmine-treated cells was carried out using the Caspase-Glo 3/7 and Caspase-Glo 9 Assay kits (G8091 and G8211, respectively, from Promega). Cell lines were plated at a density of 16,000 cells/well into white-walled, opaque 96-well plates. Twenty-four hours after plating, cells were treated with 0, 25, 50 and 100 µM harmine for 24 h. Caspase Glo reagents were added to the cells and incubated at room temperature. Luminescence was measured using a Biotek Synergy microplate reader every 20 min for 3 h.
Western blot analysis
Whole cell lysates were prepared using radioimmunoprecipitation assay (RIPA) buffer, (20 mM Tris–HCl [pH 7.5], 135 mM NaCl, 2 mM EDTA, 0.1% (w/v) sodium lauryl sulfate, 10% (v/v) glycerol, 0.5% (w/v) sodium deoxycholate, and 1% (v/v) Triton X-100). The RIPA buffer was supplemented with cOmplete™ Protease Inhibitor Cocktail (Roche), and 0.27 mM Na3VO4 and 20 mM NaF as phosphatase inhibitors. Protein concentration was determined using the Bradford dye reagent protein assay (Bio-Rad). Equal amounts of protein were resolved using 12% SDS-PAGE, and transferred to 0.45 µM polyvinylidene difluoride Immobilon-P membrane (Millipore). The transferred proteins were incubated with the following primary antibodies: PARP rabbit polyclonal antibody (#9542) from Cell Signaling Technologies at 1:1000 dilution and GAPDH mouse monoclonal antibody (#SC-47724) from Santa Cruz Biotechnology at 1:1000 dilution. After incubating the blots with the primary antibody at 4 °C for at least 8 h, they were washed and incubated with secondary antibodies at 1:10,000 dilution for 1 h at room temperature. The secondary antibodies used were LiCor IRDye® 680RD Goat anti-Mouse IgG, (925-68070), and LiCor IRDye® 680RD Goat anti-Rabbit IgG, (925-68710).
Annexin V detection
To quantify the percentage (%) of apoptosis induction in harmine-treated NB cells, the Annexin V Detection Kit APC (88-8007, Invitrogen) assay was used. Samples were prepared according to the manufacturer’s protocol and analyzed by flow cytometry. In brief, the cells were plated overnight and then treated with 0 μM or 100 μM harmine for 24 h at 37 °C. Cells were collected by centrifugation and stained with Annexin V APC/DAPI staining solution at room temperature for 20 min in the dark. Next, cells were collected by centrifugation suspended in PBS (pH 7.4) and analyzed immediately using flow cytometry Excitation/emission for Annexin APC and DAPI were 633/700 and 350/450 nm respectively.
Molecular docking of DYRK2 and harmine
The molecular docking model for DYRK2 in complex with harmine was constructed by identifying the target structures from the Protein Data Bank (PDB) . The crystal structure of human DYRK1A (PDB ID: 3ANR) in complex with harmine was previously published [13, 18]. To predict if DYRK2 also forms a complex with harmine, we generated a molecular docking model with coordinates from an existing crystal structure of DYRK2 in complex with Leucettine (PDB:4AZF) using the online docking web server SwissDock (http://www.swissdock.ch/) [18,19,20,21]. Chain A of 4AZF was isolated from the original crystal structure, and the Leucettine L41 ligand was removed. The modified structure was then entered into SwissDock, along with the chemical structure of harmine (ZINC ID: 27646846) taken from the ZINC molecule database (http://zinc.docking.org/) [19, 20]. The conformational ΔG was calculated and deemed best fit by the parameters set by SwissDock. The resulting molecular interaction model was visualized using visual molecular dynamics (VMD) .
NCBI BLAST sequence alignment
Alignment of the DYRK family sequences was carried out using the Blastp tool provided by the National Center for Biotechnology Information (NCBI) [23, 24]. The amino acid sequences for human DYRK1A and human DYRK2 used for the alignment were taken from the Basic Local Alignment Search Tool (BLAST) database, (Accession No. Q13627 and AAH06375, respectively) (https://blast.ncbi.nlm.nih.gov/Blast.cgi) [18, 24]. The query cover reported herein was given as part of the output of the Blastp program, as was the percentage of identity and the corresponding E-value.
NB public mRNA expression dataset analysis
For analysis of DYRK gene expression in human NB patients, the largest public NB cohort for which genome-wide tumor RNA-sequencing has been performed, SEQC-498, (n = 498; GSE62564) was analyzed using the R2 genomics analysis and visualization platform developed in the Department of Oncogenomics at the Academic Medical Center—University of Amsterdam (http://r2.amc.nl). Expression data for the datasets were retrieved from the public Gene Expression Omnibus (GEO) dataset on the NCBI website (http://www.ncbi.nlm.nih.gov/geo/) and analyzed as previously described .
The statistical significance for the cell viability measurements (Fig. 2b) was calculated using the GraphPad Prism 7 package (https://www.graphpad.com/) and an unpaired student’s t-test assuming the null hypothesis was performed. The change in caspase activity was quantified by calculating the average fold change, after normalizing experimental samples to their respective controls. The significance of both caspase activation and Annexin V apoptotic measurements (Figs. 3, 4) was calculated using an unpaired student’s t-test, assuming the null hypothesis. DYRK gene NB tumor mRNA expression correlation with survival probability (Fig. 7) was evaluated by Kaplan–Meier analysis using the log-rank test as described . To determine the optimal value of gene expression to set as cutoff value, all tumor samples were first sorted according to gene mRNA expression and subsequently divided into two groups. Analyses were performed on groups separated by median or average tumor mRNA expression values. DYRK mRNA expression correlation with tumor MYCN gene amplification was determined using the non-parametric, rank-based Kruskal–Wallis test. For all tests, a P-value < 0.05 was considered to be statistically significant.
Harmine induces dose- and time-dependent NB cell death
To study the effect of harmine on the morphology of NB cells, four NB cell lines were exposed to 100 μM harmine. As early as 24 h after treatment, the plated cells began to show signs of morphological changes associated with apoptosis (Additional file 1: Fig. S1). In comparison to control cells, treated cells were smaller in size, rounded in shape, and more detached from the plate surface. After 72 h, the majority of treated cells had become spherical in shape and had detached (Fig. 2a). To determine if harmine induces dose- and time-dependent cell death, the four NB cell lines were treated with increasing drug concentrations (0, 6, 12.5, 25, 50, and 100 μM) and cell viability was measured after 24, 48 or 72 h of treatment (Fig. 2b). Increased treatment length with harmine caused decreased IC50 in all four cell lines. The IC50 values for SKNBE, KELLY, and SKNFI after a 72 h treatment of harmine were 169.6 ± 0.10, 170.8 ± 0.10, and 791.7 ± 0.77 μM, respectively. The IC50 value of SKNAS after 72 h could not be calculated (for all IC50 values including 24 and 48 h time points; see Additional file 1: Table S1). The IC50 decrease was consistent with the morphological changes observed in the treated cells (Fig. 2a). Moreover, the IC50 values suggest that harmine is more toxic to MYCN-amplified NB cell lines (SKNBE and KELLY), than to NB cell lines with a normal MYCN gene copy number (SKNAS and SKNFI).
Harmine activates caspase-3/7 and caspase-9 in NB cells
The dose- and time-dependent cell death and associated morphological changes in response to harmine treatment prompted us to determine if harmine induces caspases, known to be activated during apoptosis. NB cells were treated with 0, 50, and 100 μM of harmine for 24 h, after which caspase activity was measured using the Promega Caspase-Glo 3/7 and Caspase-Glo 9 Assay kits. Caspase-3/7 activity increased with increasing harmine concentration (Fig. 3a) and this was significant for all four cell lines at the highest concentration (100 μM). The average fold changes for the caspase-3/7 activity were: SKNBE = 4.32 (P < 1.0 × 10−4), KELLY = 8.26 (P < 1.0 × 10−4), SKNAS = 6.98 (P < 1.0 × 10−4), and SKNFI = 6.77 (P < 1.0 × 10−4). Under identical cell treatment conditions, a significant increase in caspase-9 activation was observed in all four cell lines (Fig. 3b). The average fold changes for caspase-9 activity were: SKNBE = 2.51 (P < 1.0 × 10−4), KELLY = 2.96 (P < 1.0 × 10−4), SKNAS = 4.18 (P < 1.0 × 10−4), and SKNFI = 3.86 (P < 1.0 × 10−4). The results suggest that harmine triggers apoptotic cell death in NB cells.
Harmine induces progressive apoptosis in NB cells
To confirm that caspase activation leads to progressive apoptosis, harmine-treated NB cells were analyzed for caspase-mediated PARP cleavage. As shown in Fig. 4, cleaved PARP appeared in whole cell lysates of harmine-treated SKNBE and KELLY cells, but not SKNAS or SKNFI cells. These results confirm that harmine induces apoptosis, with a profound effect on MYCN-amplified NB cells (SKNBE and KELLY); both cell lines representing the most aggressive sub-types of NB tumors.
To further confirm actual apoptosis in the harmine-treated NB cells, the amount of apoptotic cells in four NB cell lines was quantified using Annexin V staining and flow cytometry (Fig. 5a). The total amount of apoptotic cells in the SKNBE cells significantly increased from 12.5% ± 0.60 (control) to 36.8% ± 14.5 in cells that had been treated with harmine (100 μM) (Fig. 5b). Similarly, harmine significantly increased the total amount of apoptotic cells in KELLY cells from 11.3% ± 2.31 (control) to 45.3% ± 6.96 in the presence of harmine (100 μM). The change in apoptotic cells in SKNAS cells was less dramatic and increased from 6.20% ± 2.24 (control) to 16.6% ± 5.31 (100 μM) and there was no change in the percentage of apoptotic cells in SKNFI cells (Fig. 5b). The results further confirm our observation that harmine appears more toxic to MYCN-amplified NB cells (SKNBE, KELLY) than to NB cells with a normal MYCN gene copy number (SKNAS, SKNFI).
Molecular docking model predicts interaction of DYRK2 with harmine
Based on literature supporting the potential interaction between harmine and DYRK2, we performed a docking simulation using SwissDock [19, 20]. To accomplish this, we utilized the crystal structures of human DYRK1A in complex with harmine (PDB ID: 3ANR) (Fig. 6a) and human DYRK2 in complex with Leucettine (PDB:4AZF) [13, 21]. Chain A of human DYRK2 was isolated and was entered into the SwissDock portal with harmine as the ligand, after removing Leucettine. The structure of harmine was provided by the ZINC molecular database (ZINC ID: 27646846). The most favorable conformation of the DYRK2/harmine complex had a ΔG value of − 6.89 kJ mol−1 (Fig. 6a). The figure displays the highest affinity binding between human DYRK2 and harmine. The crystal structure of chain A of DYRK1A in complex with harmine is provided for comparison (Fig. 6b). The two proteins were also compared by performing a sequence alignment between human DYRK1A (Accession No. Q13627) and human DYRK2 (Accession No. AAH06375) using the NCBI Blastp online tool (Fig. 6c) . The query identity was 41%, with an E-value of 4 × 10−92.
DYRK family genes are expressed in NB tumors and DYRK1B/2/3 predict poor patient outcome
Our in vitro data in NB cell lines were supported by crystal structure-derived computational docking models and suggest that targeting DYRK family kinases through harmine treatment may be a potential new route of therapy for patients with NB. This notion prompted us to investigate DYRK family kinases in human NB tumors. To accomplish this, we analyzed the mRNA expression of the DYRK gene family in SEQC-498, the largest RNASeq dataset on human NB samples in the public domain. Analyses were performed using the R2 website (see Materials and Methods). We first determined whether DYRK genes are expressed in human NB tumors (Fig. 7a). All five DYRK genes show robust expression, especially DYRK1A and DYRK2.
We next determined whether high DYRK gene expression is beneficial to NB tumorigenesis and translates into poor patient outcome (Fig. 7b–f). Kaplan–Meier graphs representing the overall survival prognosis of NB patients in the SEQC-498 cohort based on grouping of the patients according to median DYRK gene tumor mRNA expression showed that high-level tumor mRNA expression of DYRK1B, DYRK2, and DYRK3, but not of DYRK1A or DYRK4, is significantly predictive of poor patient outcome (Fig. 7b–f). Although DYRK1A and DYRK2 are the most highly expressed, only DYRK2 is prognostic of survival outcome, suggesting that in NB the deleterious effect of harmine treatment acts mostly through DYRK2. Harmine toxicity tests and apoptosis analyses on NB cell lines suggest that MYCN-amplified cell lines are more sensitive to harmine than NB cell lines with two MYCN gene copies. We were interested whether this was reflected in the correlation between DYRK gene mRNA expression and MYCN tumor amplification in human NB tumors. We found that DYRK2 and DYRK3 mRNA expression in NB tumors was significantly higher in samples with MYCN gene amplification (Fig. 8). In contrast, DYRK1B and DYRK4 tumor expression was not correlated to MYCN amplification status, and DYRK1A was actually lower in tumors with MYCN amplification. These results might indicate that the most aggressive, MYCN-amplified NB tumors that have the highest DYRK2 and DYRK3 expression could be more sensitive to harmine treatment, similar to the results we found in NB cell lines. These expression profiles along with the harmine-DYRK2 interaction model would suggest that the cytotoxic effect of harmine in NB is primarily through DYRK2 inhibition.
NB is a rare but lethal childhood tumor. Patients often present with advanced disease with survival chances below 50%, in spite of aggressive, multimodal therapy . Almost all survivors of high-stage NB suffer from late therapy effects that severely impact quality of life. Current treatment might have reached a therapeutic plateau, clarifying an urgent need for more specific, more effective, novel therapies. The use of gene-targeting treatments is still very rare, and no current therapies target MYCN, the most important NB oncogene, that is almost invariantly over-expressed in advanced NB.
We investigated the use of harmine, a tricyclic β-carboline alkaloid with known, dose-dependent in vitro cytotoxicity for cell lines from different cancer types [9,10,11]. Harmine has been shown to inhibit MAO-A activity, increase BDNF levels, and inhibit human topoisomerase I [6, 27]. It is also a potent inhibitor of the DYRK family that interferes with neurite formation and DYRK1A plays a role in neurodegenerative disorders . Although harmine displayed the highest in vitro potency against recombinant DYRK1A, it also inhibited DYRK1B and DYRK2, and, at much lower potency, DYRK4 . Recently, harmine was found to synergize with the anthracycline doxorubicin in the breast cancer cell line MCF-7  and harmine-mediated inhibition of DYRK1A destabilized epidermal growth factor receptor (EGFR) in aggressive glioblastomas and reduced EGFR-dependent glioblastoma growth . Although harmine clearly binds to and inhibits DYRK proteins, it cannot be entirely ruled out that its anticancer effects at least in part were influenced by the inhibition of other potential harmine targets such as MAO-A or topoisomerase.
We observed strong cell death in harmine-treated NB cell lines that acting in part through the activation of caspase-3/7 and caspase-9. PARP cleavage and Annexin V staining analyses showed that apoptosis progressed into late stages. Further increasing the dose range of harmine resulted in similar effects (Additional file 1: Fig. S2). Excitingly, late apoptosis was significantly more prevalent in SKNBE and KELLY, two NB cell lines with MYCN amplification and concomitant MYCN oncoprotein expression that represent high-stage NB, than in two NB cell lines, SKNAS and SKNFI, with normal MYCN copy number and without MYCN over-expression. It has been previously shown that MYCN has the ability to promote apoptosis and specific targeted therapy combinations have been tested to exploit the apoptosis-primed state of MYCN-amplified NB cells. Moreover, caspases have been shown to contribute to cell cycle regulation independent of apoptosis, and caspase inhibitors prevented cell proliferation through induction of cell cycle arrest at mitotic phase [30, 31]. Our observations in vitro were further supported by mRNA expression data we generated in human NB tumors. We found that DYRK2 and DYRK3, both harmine targets, are most highly expressed in NB patients with poor prognosis, and show MYCN-correlated expression. Although the expression of DYRK1A was also high in NB, it was not of prognostic value. These results strongly support our in vitro findings.
To verify the importance of DYRK2 in NB, we treated four NB cell lines SKNBE, KELLY, SKNAS, and SKNFI with LDN-192960, a DYRK2/haspin inhibitor and INDY, a DYRK1A/B inhibitor [13, 32]. Strikingly, inhibition of DYRK2 with LDN-192960 induced significant cytotoxic effects whereas inhibition of DYRK1A/B with INDY had no effect under identical conditions (Fig. 9), thus underlining the importance of DYRK2 in NB.
The childhood tumor NB is responsible for 16% of childhood cancer deaths. Therapy for high-stage disease that is often accompanied by tumor MYCN over-expression, is not sufficiently effective, and causes frequent severe late effects. The use of the natural product harmine that targets DYRK family proteins to treat NB cell lines in vitro showed extensive cytotoxicity, especially in cell lines with MYCN gene amplification and over-expression. These experimental data combined with our in silico modeling and epidemiologic data strongly suggest a role for DYRK2 in NB tumorigenesis and the potential for harmine as an effective treatment option for MYCN-amplified NB.
dual-specificity tyrosine phosphorylation-regulated kinase
Becker W. Emerging role of DYRK family protein kinases as regulators of protein stability in cell cycle control. Cell Cycle. 2012;11:3389–94.
Huang M, Weiss WA. Neuroblastoma and MYCN. Cold Spring Harb Perspect Med. 2013;3:a014415.
Moloudizargari M, Mikaili P, Aghajanshakeri S, Asghari MH, Shayegh J. Pharmacological and therapeutic effects of Peganum harmala and its main alkaloids. Pharmacogn Rev. 2013;7:199–212.
Oliveira CDR, Moreira CQ, de Sá LRM, de Spinosa HS, Yonamine M. Maternal and developmental toxicity of ayahuasca in Wistar rats. Birth Defects Res B Dev Reprod Toxicol. 2010;89:207–12.
Zhang L, Zhang F, Zhang W, Chen L, Gao N, Men Y, et al. Harmine suppresses homologous recombination repair and inhibits proliferation of hepatoma cells. Cancer Biol Ther. 2015;16:1585–92.
Fortunato JJ, Réus GZ, Kirsch TR, Stringari RB, Stertz L, Kapczinski F, et al. Acute harmine administration induces antidepressive-like effects and increases BDNF levels in the rat hippocampus. Prog Neuropsychopharmacol Biol Psychiatry. 2009;33:1425–30.
Farouk L, Laroubi A, Aboufatima R, Benharref A, Chait A. Evaluation of the analgesic effect of alkaloid extract of Peganum harmala L.: possible mechanisms involved. J Ethnopharmacol. 2008;115:449–54.
Nenaah G. Antibacterial and antifungal activities of (beta)-carboline alkaloids of Peganum harmala (L) seeds and their combination effects. Fitoterapia. 2010;81:779–82.
Matthay KK, Maris JM, Schleiermacher G, Nakagawara A, Mackall CL, Diller L, et al. Neuroblastoma. Nat Rev Dis Primers. 2016;2:16078.
Himpel S, Panzer P, Eirmbter K, Czajkowska H, Sayed M, Packman LC, et al. Identification of the autophosphorylation sites and characterization of their effects in the protein kinase DYRK1A. Biochem J. 2001;359(Pt 3):497–505.
Maddika S, Chen J. Protein kinase DYRK2 is a scaffold that facilitates assembly of an E3 ligase. Nat Cell Biol. 2009;11:409–19.
VerPlank JJS, Goldberg AL. Regulating protein breakdown through proteasome phosphorylation. Biochem J. 2017;474:3355–71.
Ogawa Y, Nonaka Y, Goto T, Ohnishi E, Hiramatsu T, Kii I, et al. Development of a novel selective inhibitor of the Down syndrome-related kinase Dyrk1A. Nat Commun. 2010;1:86.
Göckler N, Jofre G, Papadopoulos C, Soppa U, Tejedor FJ, Becker W. Harmine specifically inhibits protein kinase DYRK1A and interferes with neurite formation. FEBS J. 2009;276:6324–37.
Orellana EA, Kasinski AL. Sulforhodamine B (SRB) assay in cell culture to investigate cell proliferation. Bio Protoc. 2016;6:e1984.
Skehan P, Storeng R, Scudiero D, Monks A, McMahon J, Vistica D, et al. New colorimetric cytotoxicity assay for anticancer-drug screening. J Natl Cancer Inst. 1990;82:1107–12.
Vichai V, Kirtikara K. Sulforhodamine B colorimetric assay for cytotoxicity screening. Nat Protoc. 2006;1:1112–6.
Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, et al. The protein data bank. Nucleic Acids Res. 2000;28:235–42.
Grosdidier A, Zoete V, Michielin O. SwissDock, a protein-small molecule docking web service based on EADock DSS. Nucleic Acids Res. 2011;39(Web Server issue):W270–7.
Grosdidier A, Zoete V, Michielin O. Fast docking using the CHARMM force field with EADock DSS. J Comput Chem. 2011;32:2149–59.
Tahtouh T, Elkins JM, Filippakopoulos P, Soundararajan M, Burgy G, Durieu E, et al. Selectivity, cocrystal structures, and neuroprotective properties of leucettines, a family of protein kinase inhibitors derived from the marine sponge alkaloid leucettamine B. J Med Chem. 2012;55:9312–30.
Humphrey W, Dalke A, Schulten K. VMD: visual molecular dynamics. J Mol Graph. 1996;14:27.
NCBI Resource Coordinators. Database resources of the national center for biotechnology information. Nucleic Acids Res. 2017;45:D12–7.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.
Revet I, Huizenga G, Chan A, Koster J, Volckmann R, van Sluis P, et al. The MSX1 homeobox transcription factor is a downstream target of PHOX2B and activates the Delta-Notch pathway in neuroblastoma. Exp Cell Res. 2008;314:707–19.
Bewick V, Cheek L, Ball J. Statistics review 12: survival analysis. Crit Care. 2004;8:389–94.
Jahaniani F, Ebrahimi SA, Rahbar-Roshandel N, Mahmoudian M. Xanthomicrol is the main cytotoxic component of Dracocephalum kotschyii and a potential anti-cancer agent. Phytochemistry. 2005;66:1581–92.
Atteya R, Ashour ME, Ibrahim EE, Farag MA, El-Khamisy SF. Chemical screening identifies the β-carboline alkaloid harmine to be synergistically lethal with doxorubicin. Mech Ageing Dev. 2017;161:141–8.
Pozo N, Zahonero C, Fernández P, Liñares JM, Ayuso A, Hagiwara M, et al. Inhibition of DYRK1A destabilizes EGFR and reduces EGFR-dependent glioblastoma growth. J Clin Invest. 2013;123:2475–87.
Ham J, Costa C, Sano R, Lochmann TL, Sennott EM, Patel NU, et al. Exploitation of the apoptosis-primed state of MYCN-amplified neuroblastoma to develop a potent and specific targeted therapy combination. Cancer Cell. 2016;29:159–72.
Hashimoto T, Kikkawa U, Kamada S. Contribution of caspase(s) to the cell cycle regulation at mitotic phase. PLoS ONE. 2011;6:e18449.
Cuny GD, Robin M, Ulyanova NP, Patnaik D, Pique V, Casano G, et al. Structure-activity relationship study of acridine analogs as haspin and DYRK2 kinase inhibitors. Bioorg Med Chem Lett. 2010;20:3491–4.
ASB conceived the project, designed and supervised the study, and assembled the manuscript. KU performed all experiments under direct supervision of CS. DG generated the Kaplan–Meier diagrams and gene expression array data. KU and DG wrote the corresponding experimental sections and ASB assembled the manuscript and finalized all figures. All authors read and approved the final manuscript.
The authors gratefully acknowledge Dr. Rachael Sheridan, manager of the Flow Cytometry Core at the Van Andel Research Institute, for performing flow cytometry for the Annexin V APC analysis, and Dr. Jan Koster (Department of Oncogenomics, Academic Medical Center—University of Amsterdam) for providing access to and help with the R2 platform.
The authors declare that they have no competing interests.
Availability of data and materials
All expression values and other details for SEQC-498 can be obtained through its GSE identifier on the NCBI GEO website.
Consent for publication
Ethics approval and consent to participate
Approval for the use of human tissue and data was obtained from the ‘Medisch Ethische Commissie (MEC) van het AMC (Amsterdam)’, the local research and ethics committee. In addition, all analyses were in compliance with the ‘Declaration of Helsinki for Medical Research involving Human Subjects’ (http://www.wma.net/en/30publications/10policies/b3/index.html).
The study was supported by internal Michigan State University-discretional funds awarded to André S. Bachmann.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.