- 1Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, China
- 2Department of Medicinal Chemistry, School of Pharmacy, Fudan University, Shanghai, China
Inhibition of aberrant Hedgehog (Hh) pathway had been proved to be a promising therapeutic intervention in cancers like basal cell carcinoma (BCC), medulloblastoma (MB), and so on. Two drugs (Vismodegib, Sonidegib) were approved to treat BCC and more inhibitors are in clinical investigation. However, the adverse effects and drug resistance restricted the use of Hh inhibitors. In the present study, 61 synthesized compounds containing central backbone of phthalazine or dimethylpyridazine were screened as candidates of new Hh signaling inhibitors by performing dual luciferase reporter assay. Among the compounds, L-4 exhibited an IC50 value of 2.33 nM in the Shh-Light II assay. L-4 strongly inhibited the Hh pathway in vitro and blocked the Hh pathway by antagonizing the smoothened receptor (Smo). Remarkably, L-4 could significantly suppress the Hh pathway activity provoked by Smo mutant (D473H) which showed strong resistant properties to existing drugs such as Vismodegib. Orally administered L-4 exhibited prominent dose-dependent anti-tumor efficacy in vivo in Ptch+/-; p53-/- MB allograft model. Furthermore, L-4 showed good tolerance in acute toxicity test using ICR mice. These evidences indicated that L-4 was a potent, well-tolerated, orally active inhibitor of Hedgehog pathway, and might be a promising candidate in development of Hh-targeted anti-cancer drugs.
Introduction
The Hedgehog pathway is evolutionarily conserved and plays vital roles in the early development of the mammals by controlling cell proliferation and differentiation (Pomeroy et al., 2002; Jiang and Hui, 2008). It was found in a state of suppression in most adult tissues and only involved in the maintenance and repair of some tissues (Rimkus et al., 2016). The main components of the Hh pathway in mammals include the following parts: three Hh ligands (Sonic hedgehog, Indian hedgehog, Desert hedgehog) which can stimulate the Hh pathway, PTCH which represses the activity of Smo, Smo which acts as a positive regulatory protein, and the GLI family (Ruiz and Altaba, 1997). PTCH and Smo are transmembrane proteins (Roudaut et al., 2011). In the absence of Hh ligands, PTCH inhibits Smo activity by repressing its trafficking and localization to the cilia (Cheng and Bishop, 2002; Goetz and Anderson, 2010). Once the Hh ligands bind to PTCH, PTCH would leave the cilia and Smo would be able to accumulate and activate in the cilia (Sharpe et al., 2015; Zhang et al., 2017). Activated Smo subsequently starts a downstream signaling leading to the activation of transcription factors, GLI. Then, the activated GLIs translocate into the nucleus to induce the expression of various context-specific genes which regulate cellular differentiation, proliferation and survival (Scales and de Sauvage, 2009; Stecca and Ruiz, 2010).
Aberrant activation of Hh pathway signaling leads to abnormal development of the tissues and is linked to many malignant tumors like BCC, MB, and so on (Epstein, 2008; Rimkus et al., 2016; Tan et al., 2018). MB, known as a primitive neuroectodermal tumor, is the most common malignant pediatric brain tumor (Romer et al., 2004). It happens in the cerebellum and accounts for about 20% of all pediatric brain tumors (Ellison et al., 2003). Previous studies had shown that Hh pathway was closely related to MB. Mutant PTCH1 was reported to cause the human developmental disorder Gorlin syndrome which increased a notable risk of advanced BCC and MB. This indicated that Hh pathway contributed to the formation of MB (Hahn et al., 1996). Actually, further researches showed abnormal expression of PTCH1, SMO, and SUFU in a large percentage of spontaneous MB (Epstein, 2008; Kool et al., 2008). Up to now, MB are termed four subgroups: WNT, SHH, group 3, and group 4 (Taylor et al., 2012). Among them, Hh pathway subtype accounts for 30% of MB (Kool et al., 2012; Korshunov et al., 2012). Thus, research and development of inhibitors targeting the Hh pathway might be a promising strategy for MB therapy (Rimkus et al., 2016).
Cyclopamine is the first Hh inhibitor to be found which advances our research of the Hh pathway as a therapeutic target substantially (Incardona et al., 1998; Tremblay et al., 2009). It is an allosteroid alkaloid extracted from the veratrine (Incardona et al., 1998). Cyclopamine could directly combine with the transmembrane domain of Smo proteins and prevent the conformational change and activation of Smo (Chen et al., 2016). However, the poor potency and oral solubility limit its application in clinical therapy (Tremblay et al., 2009). In recent years, great effort has been taken predominantly on targeting Smo. So far, a variety of Smo inhibitors have been reported. Among them, Vismodegib (GDC-0449) (Robarge et al., 2009; Kumar et al., 2017) and Sonidegib (NVP-LDE225) (Pan et al., 2010; Jalili et al., 2013) have been approved by FDA for the treatment of BCC. However, Vismodegib and Sonidegib showed undesirable adverse effect and drug resistance subsequently in clinic therapy (Scales and de Sauvage, 2009). According to the research, there are 15–33% untreated BCC patients harbor Smo mutations and the percentage increases to 69–77% in resistant tumors (Dong et al., 2018). The MB relapsed after a few months of Vismodegib treatment and it was found that the mutation of Smo D473H gave rise to the drug resistance (Yauch et al., 2009). Thus, the binding ability of Hh inhibitors to D473H is important to overcome the drug resistance.
NVP-LEQ506 (Peukert et al., 2013) and Taladegib (LY-2940680) (Bendell et al., 2018) are two representative Hh inhibitors in clinic trial. NVP-LEQ506 was developed by Novartis and LY-2940680 was developed by Lilly United States and Ignyta. LY-2940680 and Anta XV (lead compound of NVP-LEQ506) are both benzylphthalazine derivatives which were reported to be potent Hh pathway inhibitors with low nanomolar affinity for Smo. Based on these previous studies, we designed and synthesized a series of small molecular compounds containing central backbone of phthalazine or dimethylpyridazine as candidates of new Hh signaling inhibitors. The piperazine linker of lead compound Anta XV was replaced by piperidin-4-amine moiety and the benzyl group was optimized with different aromatic rings to generate compounds series A. We applied the similar strategy to replace the 4-methylamino-piperidine linker of LY-2940680 with pyrrolidin-3-amine moiety and the 1-methyl-1H-pyrazol group with benzene ring to generate compounds series L. Furthermore, in some compounds, the bicyclic phthalazine core was replaced with dimethylpyridazine to mimic the pharmacophore of NVP-LEQ506. After screening for inhibitive effects on Hh signaling pathway of these compounds, L-4 was found to be the most potent compound. Herein, evaluation of anti-tumor effects of L-4 in vitro and in vivo was conducted.
Materials and Methods
Reagents
Vismodegib, Taladegib, and SAG were purchased from Selleck (TX, United States). NVP-LEQ506 was bought from MedChemExpress (NJ, United States). L-4 was synthesized in-house. For in vitro experiments, compounds were dissolved in DMSO (Beyotime, Shanghai, China) and stock solutions were stored at -20°C. For in vivo experiments, Vismodegib, Taladegib and L-4 were formulated in 0.5% methylcellulose (Aladdin, Shanghai, China).
Cell Lines and Cell Culture
Shh-Light II cell line was a gift of Professor. Philip Beachy in Stanford University and cells were maintained in DMEM (Hyclone, UT, United States) containing 10% fetal bovine serum (FBS) (Biological Industries, Israel), zeocin (Invitrogen, CA, United States) 0.15 mg/mL and G418 (Invitrogen) 0.4 mg/mL at 37°C in a 5% CO2 atmosphere. HEK293 human epithelial kidney cell was purchased from Cell Bank of China Science Academy (Shanghai, China) and cells were maintained in DMEM containing 10% FBS, streptomycin 100 μg/mL and penicillin 100 U/mL as recommended at 37°C in a 5% CO2 atmosphere.
Animals
All animal (purchased from Shanghai SLAC Laboratory Animal Co., Ltd., Shanghai, China) care and experimental protocols for this study complied with the Chinese regulations and the Guidelines for the Care and Use of Laboratory Animals drawn up by the National Institutes of Health (USA) and were approved by the Institutional Animal Care and Use Committee of the East China Normal University. ICR mice and BALB/c nude mice (6–8 weeks old) were purchased from the Shanghai SLAC Laboratory Animal CO., LTD., Mice were maintained on a 12:12 light-dark cycle in a temperature-controlled (21∼23°C) and SPF conditional room. Mice were provided with standard rodent chow and water ad libitum. All animals were acclimatized for a week before beginning the study.
Gli-Luc Reporter Assay
Shh-Light II cells were seeded in 96-well plants at a density of 4000 cells/200 μL. After 48 h of various compounds treatment as indicated, the luciferase activity was examined with Dual-Luciferase Assay Reporter System (Promega, Madison, WI, United States) following the manufacturer’s instructions. The firefly luciferase values were normalized to Renilla values. GraphPad Prism software was used to generate dose–response curves and IC50 values with five parameters logarithmic non-linear regression fitting.
Medulloblastoma Primary Cells Culture and Brdu Assay
Ptch+/-; p53-/- mice were gained from hybridization of Ptch+/- mice and p53-/- mice. p53-/- mice were kindly provided by Professor Wei Zhang (East China Normal University, Shanghai, China) and Ptch+/- mice were bought from Model Animal Research Center of Nanjing University (Nanjing, China). The primary spontaneously MB in Ptch+/-; p53-/- mice were collected and allografted into BALB/c nude mice. After development in the nude mice, tumors were harvested. The cells were maintained in Neurobasal A medium (Invitrogen) containing B-27 supplement (Invitrogen), EGF 20 ng/ml (Invitrogen), bFGF 20 ng/ml (Invitrogen), non-essential amino acids (Invitrogen), N-acetylcysteine 60 μg/ml (Sigma-aldrich, St. Louis, MO, United States).
Medulloblastoma cells were seeded in 96-well plates at a density of 20000 cells/200 μL. After treated with various concentrations of L-4 for 36 h, the Brdu assay was conducted using the Brdu Cell Proliferation Kit (Abcam, Cambridge, United Kingdom) according to manufacturer’s instructions.
Microscopy Analysis of Fluorescent BODIPY-Cyclopamine Competition Assay
Human wild type Smo (hSmo) expression vector is a gift from Professor Ke Yu (Fudan University, Shanghai, China). Microscopy analysis of fluorescent BODIPY-cyclopamine competition assay was performed as following. In brief, the HEK293 cells were seeded in 6-well plates and transfected with wild type hSmo expression vector. Two days after transfection, the cells were seeded in 24-well plates and incubated with compounds. After incubation for 10 h, the cells were washed with PBS twice, fixed with 4% paraformaldehyde in dark for 10 min, washed with PBS again and incubated with Triton X-100 solution for 10 min. The cells were then stained with DAPI (Sigma-aldrich) and photographed using a fluorescence microscope.
Flow Cytometry Analysis of Fluorescent BODIPY-Cyclopamine Competition Assay
Fluorescent BODIPY-cyclopamine competition assay was adopted as following. In brief, HEK293 cells were seeded in 6-well plates and transfected with wild type hSmo expression vector. Two days later, cells were seeded in 24-well plates and treated with the compounds. After incubation for 10 h, the cells were harvest and tested by flow cytometry. The green fluorescence of BODIPY-cyclopamine bound to Smo was examined by a flow cytometer (Guava EasyCyte 6HT-2 L, Merck Millipore, Billerica, MA, United States). At least 5000 cells were analyzed for each data point. FACS data were analyzed with InCyte software.
Antagonizing Drug-Resistant Assay
Drug-resistant D473H Smo expression vector is a gift from Professor Ke Yu (Fudan University, Shanghai, China). Shh-Light II cells were transfected with wild type and mutant D473H Smo expression vectors. Then the cells were seeded in 96-well plates and treated with Vismodegib and L-4 at various concentrations. After incubation for 48h, a dual luciferase reporter assay (Promega, Madison, WI, United States) was performed following the manufacturer’s instructions. 293T cells transfected with mutant Smo D473H were treated with Vismodegib and L-4 for 48 h, then cells were harvested and subjected to WB analysis.
Molecular Modeling
Molecular docking was carried out by using Glide module in Schrödinger 2013 software package. The crystallographic structure of the human Smo 7TM receptor was retrieved from the RCSB protein data bank (ID: 4JKV). The SWISS-MODEL online server was employed to construct the mutant D473H by replace the specific residue Asp473. The WT and D473H structures were minimized by using AMBER16 software and then were prepared by using the Protein Preparation Wizard module. Two ligands were selected in the docking study: L4 and Vismodegib. Before docking, they have been processed by using the Ligprep module in Schrödinger 2013. These two molecules were docked into the binding pocket of the WT and D473H protein structure. The best docked poses were select for further analysis.
Western Blot Analysis
The proteins of whole cell lysates were prepared in RIPA buffer containing protease and phosphatase inhibitors (Thermo Fisher Scientific, Waltham, MA, United States), and quantified by the BCA method. Immunoblotting analysis of Gli1 and Ptch1 were performed as previously described (Dorward et al., 2016). The expression of β-actin was used as loading control. Antibodies against Smo and β-actin were purchased from Santa Cruz Biotechnology (CA, United States). Antibody against Ptch1 was purchased from Abcam. Antibody against Histone H3 was purchased from CST (Danvers, MA, United States).
Cell fractionation was performed as following. The collected cells were lysed by cytoplasmic lysis buffer [50 mM Tris-HCl (PH 7.4), 140 mM NaCl, 1.5 mM MgCl2, and 0.5% NP-40 (1.06 g/ml)] for 20 min. Then the lysate was centrifuged for 10 min at 3000 rpm at 4°C and the supernatant was collected as cytoplasmic lysate. The pellet was washed once by cytoplasmic lysis buffer and dissolved by RIPR buffer (Thermo Fisher Scientific, Waltham, MA, United States) for 20 min. Then the lysate was centrifuged for 10 min at 12000 rpm at 4°C and the supernatant was collected as nuclear lysate.
Medulloblastoma Allograft Model
The primary intracranial MB spontaneously developed in Ptch+/-; p53-/- mice were harvested and subcutaneously allografted into athymic nude mice. The tumors in nude mice were further collected and cut into 1 mm3 fragments. Then the fragments were inoculated subcutaneously into the right flank of other nude mice. When the tumor volume reached 150–200 mm3, the mice were randomized grouped into control and treatment groups to receive treatment accordingly. L-4 was orally administered at 5, 10, and 20 mg/kg once a day for 13 days. Taladegib was used as a positive control and orally administered at 20 mg/kg once a day. The tumor growth was recorded with the measurement of length (L) and width (W) by caliper every other day and calculated as tumor volume (V) = L × W2/2. Meanwhile, the body weights of mice were recorded. The TGI was calculated as TGI% = [1-(mean tumor volume of the treated group on the last day-mean tumor volume of the treated group on day 0)/(mean tumor volume of the control group on the last day-mean tumor volume of the control group on day 0)] × 100%. The tumor regression of individual mouse was defined as the tumor volume at the end was less than the tumor volume when treatment was initiated.
Acute Toxicity Test
According to guideline 425 Organisation for Economic Co-operation Development [OECD], 2008, L-4 was administered in a single dose of 2000 mg/kg by oral in five female ICR mice (n = 5). A negative control was established by administering a 0.5% CMCNa solution to a female mouse (n = 1). Body weight of the animals were recorded shortly before the administration of the L-4 and then daily, during the treatment period. After administration, animals were observed individually during the first 30 min and then daily for 14 days. Observations included mortality and changes in skin and fur, eyes and mucous membranes, and also respiratory, circulatory, autonomic and central nervous systems, and somatomotor activity and behavior pattern. Attention should be directed to observations of tremors, convulsions, salivation, diarrhea, lethargy, sleep, and coma.
Statistical Analysis
The statistical significance of differences between groups was evaluated by the unpaired Student’s t test and indicated with ∗∗P < 0.01, ∗P < 0.05. All statistical tests were two sided.
Results
Screening Hh Signaling Inhibitory Activities of the Synthesized Compounds
To evaluate the inhibitory effects of all the compounds, dual luciferase reporter assays were conducted using Shh-Light II cells, which were NIH-3T3 cells stably fused with a Gli-responsive firefly luciferase and Renilla-luciferase (Xin et al., 2014). As shown in Table 1, 61 compounds were screened and 14 compounds showed good inhibitory potency. Then, we selected 5 representative compounds to further check the in vitro concentration-inhibition IC50 values. NVP-LEQ506 and LY-2940680 were used as positive controls. The IC50 values ranged from 1 to 50 nM (Table 2). L-4 (Figure 1) showed the best inhibitory potency with an IC50 of 2.33 nM.
Table 1. Inhibitory rate of compounds to the Hh pathway tested by dual luciferase reporter assays in Shh-Light II cells.
Table 2. IC50 value of compounds to the Hh pathway tested by dual-luciferase reporter assays in Shh-Light II cells.
L-4 Inhibits the Proliferation of Medulloblastoma Cells With Activation of Hh Pathway Signaling
The potency of L-4 in inhibiting the proliferation of cancer cells related to activation of Hh pathway was tested. The Ptch+/-; p53-/- MB mouse model, a generally accepted Hh-driven mouse MB model which had been used to evaluate the anti-cancer effect of Smo inhibitors (Wetmore et al., 2001), was established. Spontaneous MB cells were then isolated from nude mice for testing the antiproliferative effect of compound L-4 using Brdu assay. Vismodegib and LY2940680 were used as positive controls. The results showed that L-4 dose-dependently inhibited the proliferation of MB cells in a similarity to that of Vismodegib and LY2940680 (Figure 2A). The IC50 of L-4 was 31.36 nM and the IC50 of Vismodegib and LY2940680 were 37.65 and 35.5 nM, respectively. Also, Western bolt results showed that L-4 dose-dependently inhibited the expression of Gli1 and Ptch1 in MB cells (Figure 2B).
Figure 2. L-4 suppressed the proliferation and Hh pathway of Ptch+/-; p53-/- medulloblastoma (MB) cells in vitro. (A) L-4 inhibited the proliferation of Ptch+/-; p53-/- MB cells assessed by Brdu analysis. MB cells were exposed to various concentrations of L-4 for 36 h. Vismodegib and LY-2940680 were used as positive controls. (B) L-4 inhibited the expression of Gli1 and Ptch1 in primary Ptch+/-; p53-/- MB cells. MB cells were treated with various concentration of L-4 after 48 h and collected for lysis. Protein extracts were probed with specific antibodies against Gli1, Ptch1, and Actin (as a loading control). Cells isolated from normal mouse brain was marked as “Normal Brain.”
L-4 Binds With Smo and Restrains the Expression of Gli1 and Ptch1
To check whether the inhibitory potency of L-4 against the Hh pathway is due to targeting Smo, direct interaction of L-4 with Smo was analyzed by using BODIPY-cyclopamine, a kind of cyclopamine derivative with fluorescent label (Solinas et al., 2012; Sharpe et al., 2015). Firstly, HEK293 cells were transfected with the wild type Smo expression vector and the effect of transfection was examined by western blot (Figure 3A). Then, fluorescent BODIPY-cyclopamine competition assay was performed. In the presence of 10 μM L-4, cellular fluorescence of BODIPY-cyclopamine was not observed by the fluorescence microscopy analysis (Figure 3B). It was similar to the positive control, cyclopamine (Schaefer et al., 2013). As shown in Figures 3C,D, results of quantitative analysis using flow cytometry suggested that there was a dose-dependent replacement of BODIPY-cyclopamine by L-4. Taken together, our data demonstrated that L-4 was able to displace the BODIPY-cyclopamine from Smo and inhibited the Hh signaling pathway by targeting Smo. Moreover, results of western blot (Figure 3E) suggested that L-4 dose-dependently reduced the expression of Gli1 and Ptch1 in Shh-Light II cells.
Figure 3. L-4 inhibited the Hh pathway by binding to Smo. (A) Western blot analysis of Smo expression in HEK293 cells. HEK293 cells were transfected with wild type hSmo expression vector and the Smo protein expression level was examined by western blot. Actin was used as a loading control. (B) BODIPI-cyclopamine competition analysis was examined by fluorescent microscope. Photographs were representatives from three independent experiments. (C,D) BODIPI-cyclopamine competition analysis was examined by flow cytometry. The results were expressed as mean ± SEM of three independent experiments (n = 3). (E) L-4 inhibited the expression of Gli1 and Ptch1 in Shh-Light II cells. Shh-Light II cells were treated with 100 nM SAG and various concentrations of L-4 for 48 h were harvested and subjected to WB analysis. (F) L-4 decreased the Gli1 expression both in the cytoplasm and the nucleus in Shh-Light II cells. The cells were treated with 100 nM SAG and various concentrations of L-4 for 48 h. After treatment, cell fractionation was performed and lysates were subjected to WB analysis.
It had been shown that Gli1 translocated from the cytoplasm to the nucleus to induce a series of downstream genes expression (Sharpe et al., 2015). Thus, the shuttle of Gli1 in intercellular compartments after treated with L-4 was observed. As showed in Figure 3F, SAG induced the shuttle of Gli1 to the nucleus in Shh-Light II cells. While, after treatment with different concentrations of L-4 for 48 h, a strong decrease of Gli1 both in the cytoplasm and the nucleus could be observed. The results indicated that the inhibition of Smo by L-4 not only had a negative effect on translocation of Gli1, but also inhibited the Gli1 expression.
L-4 Inhibits Drug-Resistant Mutant Smo
Whether L-4 could inhibit the aberrant Hh pathway provoked by mutant Smo D473H was checked. Mutant Smo D473H is an important reason for giving rise to the resistance to the current existing drugs (Dijkgraaf et al., 2011; Infante et al., 2016). Compared to the conspicuous inhibitory effect on the Shh-Light II cells stimulated by Smo wild type expression vectors, L-4 showed similar inhibitory effects on the Shh-light II cells with overexpression of mutant Smo D473H (Figure 4A). The IC50 values of L-4 in inhibiting wild type and D473H were 23.46 and 24.55 nM, respectively. On the contrary, the inhibiting effects of Vismodegib on Hh pathway provoked by overexpression of D473H were reduced (Figure 4B). In the present study, the IC50 values of Vismodegib in inhibiting wild type and D473H were 21.63 and 326.8 nM, respectively. To further confirm these results, we transfected mutant Smo D473H into 293T cells and there was a corresponding activation in Gli1 and PTCH1. Although these proteins were largely resistant to Vismodegib, they were dose dependently inhibited by L-4 (Figure 4C). To sum up, we indicate that L-4 could effectively target the mutant Smo D473H.
Figure 4. L-4 inhibited the Gli-luciferase activity stimulated by wild type Smo and mutant Smo D473H. Light II cells transfected with wild type Smo or Smo D473H plasmids were exposed to L-4 for 48 h and then were harvested for dual luciferase reporter assays. The results were expressed as the mean ± SEM. (A) L-4 dose-dependently inhibited the Gli luciferase activity proved by wild type Smo and mutant Smo D473H. (B) Vismodegib failed to inhibited the Gli luciferase activity initiated by mutant Smo D473H. (C) L-4 dose-dependently inhibited the expression of Gli1 and Ptch1 in 293T cells transfected with mutant Smo D473H. Transfected cells were treated with Vismodegib and L-4 for 48 h. Cells were harvested and subjected to WB analysis.
The Binding Mode of L-4 With Wild Type and Mutant Smo
The molecular docking simulation was used to study the binding modes of L-4 and Vismodegib with wild type and mutated D473H Smo. The predicted binding mode showed that the interaction between L-4 and wt-Smo was much stronger. As shown in Figure 5A, the inhibitor L-4 forms three critical hydrogen bonds with Smo protein: one hydrogen bond from the amide oxygen atom to Asn219 residue, two hydrogen bonds between the dimethylpyridazine core and Arg400. Furthermore, the phenyl ring and pyridazine ring also interacted with the benzene ring of Phe484 and Tyr394 via π-π stacking. It suggested the interaction between L4 and wt-Smo was much stronger. These interactions have hardly changed after the D473H mutation (Figure 5B), the distance of Hydrogen bond is changed from 2.9, 3.4, and 3.3 Å to 3.0, 3.1, and 3.4 Å, respectively. The amide oxygen atom of Vismodegib also formed one hydrogen bond with Arg400 (2.9 Å) in wt-Smo, which disappeared after D473H mutation (5.2 Å). In fact, the D473H mutation (relatively larger side chain) could change the structure of the binding pocket, especially the position of Arg400. There is only one hydrogen bond between Vismodegib and the wt-Smo, which could be easily disrupted by D473H mutation and leaded to drug resistance (Figures 5C,D). The Computational predicted binding modes were in good agreement with the experiment results.
Figure 5. Computational predicted binding mode of L-4 and Vismodegib with wild type and mutated D473H Smo. (A) L-4-Smo (WT), (B) L-4-SMO (D473H), (C) Vismodegib -Smo (WT), and (D) Vismodegib-Smo (D473H). All pictures were generated by using Pymol.
L-4 Inhibits the Hh Signaling Pathway Activity in vivo
Subcutaneous Ptch+/-; p53-/- MB allograft model was used to assess the in vivo anti-tumor efficacy of L-4. L-4 was orally administered at a dose of 5, 10, or 20 mg/kg once a day. As showed in Figures 6A,C,D, L-4 showed dose-dependent anti-tumor activity after 13 days of oral administration. Even at the lowest dose level of 5 mg/kg, tumor growth was significantly inhibited by L-4 compared to the vehicle control. TGI were 53.3, 84.3, and 93.3% under the dose of 5, 10, and 20 mg/kg, respectively. Body weights did not change significantly in mice treated with L4 at all doses (Figure 6B). Four hours after the last dose of vehicle, L-4 and LY-2940680 treatment, tumors were harvested and Gli1 and Pthc1 expression levels were checked as pharmacodynamic markers of Hh pathway. The results showed that L-4 could dose-dependently decrease the expression of Gli1 and Ptch1 in tumors (Figure 6E). Herein, these studies suggested that L-4 successfully induced TGI in vivo as it worked in vitro by inhibiting Hh pathway.
Figure 6. Anti-tumor activity of L-4 in Ptch+/-; p53-/- MB allograft model. Mice bearing in MB were orally administered with LY-2940680 and L-4 once a day for 13 days. (A) Tumor volume for indicated days was showed as mean ± SEM (n = 5). (B) Mean body weights were presented as mean ± SEM (n = 5). (C) The pictures of tumors. (D) Tumor weight. (E) L-4 inhibited the Gli1 and PTCH1 expression in the tumors. After 4 h of last dosage, tumor samples were collected for examining the expression of Gli1 and PTCH1 protein expression. (F) Body weight for the ICR mice in the acute toxicity test. Number 0 was the negative control mouse orally administered with 0.5% CMCNa solution. Number 1–5 were the mice orally administered with a single dose of 2000 mg/kg L-4, separately. ∗p < 0.05, ∗∗p < 0.01 vs. control.
L-4 Shows No Toxicity in 14-Day Test of ICR Mice
L-4 was run through acute toxicity test as part of the preclinical safety assay. In 14-day test using ICR mice, mice receiving a single dose of 2000 mg/kg L-4 did not exhibit any toxic symptom (Figure 6F). Briefly, it could be predicted that the LD50 of the L-4 might be higher than 2000 mg/kg.
Discussion
The abnormal expression of Hh pathway, which has a key role in the development and maintenance of many organs and tissues, has been linked to tumorigenesis. It had been proved that aberrant in this pathway can lead to BCC and MB. Recent researches demonstrated that Hh pathway could also be linked to many other cancers, including pancreatic (Rajurkar et al., 2012; Miyazaki et al., 2016; Kumar et al., 2017), colon (Park et al., 2011), gastric (Tang et al., 2018), lung (Zhu et al., 2017), breast (Song et al., 2016), and prostate cancers (Chen et al., 2011), leukemias and multiple myeloma (Blotta et al., 2012). As a result, development of novel Hh inhibitors is of great meaning and urgency. Here, we designed and synthesized a series of small molecular compounds containing central backbone of phthalazine or dimethylpyridazine as candidates of new Hh signaling inhibitors. After evaluation with dual luciferase reporter assays using Shh-Light II cells, among the all compounds, L-4 showed the best inhibitory potency with an IC50 of 2.33 nM.
In vitro, evaluating by Brdu assay, L-4 showed obvious inhibitory effect on the proliferation of primary MB cells with nanomolar IC50 just as Vismodegib and Taladegib. In addition, western blot analysis of primary MB cells treated with L-4 revealed that the expression of Hh pathway were repressed. These results confirmed that L-4 had inhibition effects on proliferation of MB cells by targeting Hh pathway.
BODIPI-cyclopamine is the derivative of cyclopamine with a fluorescent tag. Those inhibitors which share the same binding domain with cyclopamine can competitively prevent BODIPI-cyclopamine from binding to Smo protein. Correspondingly, the reduce of fluorescence signal of BODIPI-cyclopamine would be detected in the microscopy and flow cytometry analysis. In this occasion, we are able to confirm the binding domain between the inhibitor and Smo. To check whether L-4 could directly bind to Smo and to find the domain of binding between L-4 and Smo, BODIPI-cyclopamine competition assay was conducted in the present study. Both the results of fluorescence microscope analysis and the results of flow cytometry analysis showed the concentration-dependent competitive relationship between L-4 and BODIPI-cyclopamine. The results suggested that L-4 could directly bind to the transmembrane domain of Smo and occupied the same pocket as cyclopamine. As previous reported, in activating Hh pathway, Gli1 shuttles from the cell cytoplasm to nucleus to start the downstream signal in nucleus (Scales and de Sauvage, 2009). Under the treatment of L-4, there was an obvious decrease of Gli1 in both cytoplasm and nucleus. The reduce of Gli1 in the nucleus might be resulted from block of the Hh pathway by L-4. Under blocking of the Hh pathway, more Gli1 protein was restrained by the SuFu, and thus few Gli1 could go through the nucleus. Since Gli1 is also a Hh target gene and significantly influenced by activation of Hh pathway, there might be a negative feedback after L-4 binds to the Smo. Less Gli1 shuttles into nucleus might lead to reduction in the gene expression of Gli1 which resulted in the decrease of expression of Gli1 protein. Therefore, the level of Gli1 in cytoplasm also decreased.
We further checked the ability of L-4 in binding with mutant Smo D473H. In the antagonizing drug-resistant assay, L-4 exhibited equivalent potency in inhibiting the Gli-luciferase activity, and downstream biomarkers provoked by Smo wild type and Smo D473H. On the contrary, Vismodegib showed ameliorated ability in inhibiting the Gli-luciferase activity provoked by Smo D473H. In 293T cells, L-4 blocked the downstream Hh pathway biomarkers driven by mutant Smo D473H while Vismodegib failed. The results suggested possible better therapeutic effects of L-4 than Vismodegib in clinic because many patients suffered from drug resistance after treated by Vismodegib. Take molecular modeling into consideration, strong hydrogen bonds between L-4 and Smo protein suggested that the binding between L-4 and Smo might be stable. Particularly, results of calculational chemistry also showed that L-4 could equally bind to the wild type and mutant Smo, which supported results of our drug-resistant study. Vismodegib showed much weaker binding ability to mutant Smo. The docking scores of Vismodegib to wild type Smo and mutant Smo were -8.62 and -7.39, respectively. The inhibition ability of Vismodegib in antagonizing mutant D473H is more than ten folds weaker than that of wild type. The results were consistent with the IC50 values of Vismodegib in inhibiting wild type and mutant D473H transfected in the Shh-Light II cells, which were 21.6 and 326.8 nM, respectively. The mutation (relatively larger side chain) of Smo can change the structure of the binding pocket, especially the position of Arg400. There is only one hydrogen bond between Vismodegib and the Smo which can be easily disrupted by the mutation, while the interaction between L-4 and Smo is much stronger, thus their binding is not greatly affected by the mutation.
Of note, our tumor profiling data indicated that orally administration of L-4 showed remarkable dose-dependent anti-cancer effects in the Ptch+/-; p53-/- MB allograft model without inducing loss of body weight. At the dose of 20 mg/kg, tumors in all five mice almost stopped growing in the 13-day treatment. The expression levels of Gli1 and Ptch1 in tumors reduced along with augment of the dose of L-4. It suggested that L-4 inhibited MB tumor growth by inhibiting Hh pathway and Gli1 and Ptch1 might be important pharmacodynamic markers. Besides, according to the results of acute toxicity test, L-4 was well tolerated in adult mice and no obvious side effects were observed.
In all, L-4 exhibited excellent drug-like properties. L-4 could bind to both wild type and mutant D473H Smo and then inhibited the escapement of Gli from SUFU protein complex and translocation of Gli into the nucleus which would induce the expression of various context-specific genes related to regulation of cellular differentiation, proliferation, and survival (Figure 7). It had considerable anti-tumor efficacy in vitro and in vivo. It was highly soluble and with limited toxicity. Nowadays, many preclinical and clinical studies are conducted to develop novel Hh inhibitors or combine Hh inhibitors with other clinical medicine in hope to cure cancers with overexpressed Hh pathway (Riobo et al., 2006; Wang et al., 2012; Amakye et al., 2013). Although previous studied Hh inhibitors failed in clinical trials in all other tumors except BCC and MB, more effort should be taken to research and develop new Hh inhibitors. Since L-4 could bind with wild type and mutant Smo in the aberrant Hh pathway, it is reasonable that L-4 could be a promising candidate as new anti-cancer agent and create more opportunities in cancer therapy.
Figure 7. L-4 inhibited the Hedgehog pathway by targeting Smo. L-4 bound with wild type and mutant D473H Smo to prevent Gli from being released from SUFU protein complex and then inhibited the translocation of Gli into the nucleus to induce the expression of various context-specific genes which regulate cellular differentiation, proliferation, and survival.
Author Contributions
XZ, YF, WZ, and XD conceived and designed the study. MZ, HW, and CW performed the experiments. TZ performed the molecular docking analysis. MZ wrote the manuscript. XZ, WZ, and XD supervised the whole experimental work and revised the manuscript. All authors read and approved the final manuscript.
Funding
This work was supported by National Natural Science Foundation of China (81872496, 21572037, 21572263, and 21572067), the Science and Technology Commission of Shanghai Municipality (16DZ2280100), and pilot project of the Chinese Academy of Sciences (XDA12020356).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Abbreviations
BCA, bicinchoninic acid; BCC, basal cell carcinoma; CMCNa, carboxymethyl cellulose sodium; GLI, glioma associated oncogene; Hh, hedgehog; ICR mice, institute of cancer research mice; MB, medulloblastoma; PTCH, patched; SHH, sonic hedgehog; Smo, smoothened; SPF, specific pathogen-free; TGI, tumor growth inhibition; WNT, wingless.
References
Amakye, D., Jagani, Z., and Dorsch, M. (2013). Unraveling the therapeutic potential of the Hedgehog pathway in cancer. Nat. Med. 19, 1410–1422. doi: 10.1038/nm.3389
Bendell, J., Andre, V., Ho, A., Kudchadkar, R., Migden, M., Infante, J., et al. (2018). Phase I Study of LY2940680, a Smo antagonist, in patients with advanced cancer including treatment-naive and previously treated basal cell carcinoma. Clin. Cancer Res. 24, 2082–2091. doi: 10.1158/1078-0432.CCR-17-0723
Blotta, S., Jakubikova, J., Calimeri, T., Roccaro, A. M., Amodio, N., Azab, A. K., et al. (2012). Canonical and noncanonical Hedgehog pathway in the pathogenesis of multiple myeloma. Blood 120, 5002–5013. doi: 10.1182/blood-2011-07-368142
Chen, B., Trang, V., Lee, A., Williams, N. S., Wilson, A. N., Epstein, E. H., et al. (2016). Posaconazole, a second-generation triazole antifungal drug, inhibits the hedgehog signaling pathway and progression of basal cell carcinoma. Mol. Cancer Ther. 15, 866–876. doi: 10.1158/1535-7163.MCT-15-0729-T
Chen, M., Carkner, R., and Buttyan, R. (2011). The hedgehog/Gli signaling paradigm in prostate cancer. Expert Rev. Endocrinol. Metab. 6, 453–467. doi: 10.1586/EEM.11.2
Cheng, S. Y., and Bishop, J. M. (2002). Suppressor of Fused represses Gli-mediated transcription by recruiting the SAP18-mSin3 corepressor complex. Proc. Natl. Acad. Sci. U.S.A. 99, 5442–5447. doi: 10.1073/pnas.082096999
Dijkgraaf, G. J., Alicke, B., Weinmann, L., Januario, T., West, K., Modrusan, Z., et al. (2011). Small molecule inhibition of GDC-0449 refractory smoothened mutants and downstream mechanisms of drug resistance. Cancer Res. 71, 435–444. doi: 10.1158/0008-5472.CAN-10-2876
Dong, X., Wang, C., Chen, Z., and Zhao, W. (2018). Overcoming the resistance mechanisms of Smoothened inhibitors. Drug Discov. Today 23, 704–710. doi: 10.1016/j.drudis.2018.01.012
Dorward, H. S., Du, A., Bruhn, M. A., Wrin, J., Pei, J. V., Evdokiou, A., et al. (2016). Pharmacological blockade of aquaporin-1 water channel by AqB013 restricts migration and invasiveness of colon cancer cells and prevents endothelial tube formation in vitro. J. Exp. Clin. Cancer Res. 35:36. doi: 10.1186/s13046-016-0310-6
Ellison, D. W., Clifford, S. C., Gajjar, A., and Gilbertson, R. J. (2003). What’s new in neuro-oncology? Recent advances in medulloblastoma. Eur. J. Paediatr. Neurol. 7, 53–66. doi: 10.1016/s1090-3798(03)00014-x
Epstein, E. H. (2008). Basal cell carcinomas: attack of the hedgehog. Nat. Rev. Cancer 8, 743–754. doi: 10.1038/nrc2503
Goetz, S. C., and Anderson, K. V. (2010). The primary cilium: a signalling centre during vertebrate development. Nat. Rev. Genet. 11, 331–344. doi: 10.1038/nrg2774
Hahn, H., Wicking, C., Zaphiropoulous, P. G., Gailani, M. R., Shanley, S., Chidambaram, A., et al. (1996). Mutations of the human homolog of Drosophila patched in the nevoid basal cell carcinoma syndrome. Cell 85, 841–851. doi: 10.1016/S0092-8674(00)81268-4
Incardona, J. P., Gaffield, W., Kapur, R. P., and Roelink, H. (1998). The teratogenic veratrum alkaloid cyclopamine inhibits Sonic hedgehog signal transduction. Development 125, 3553–3562.
Infante, P., Alfonsi, R., Ingallina, C., Quaglio, D., Ghirga, F., D’Acquarica, I., et al. (2016). Inhibition of Hedgehog-dependent tumors and cancer stem cells by a newly identified naturally occurring chemotype. Cell Death Dis. 7:e2376. doi: 10.1038/cddis.2016.195
Jalili, A., Mertz, K. D., Romanov, J., Wagner, C., Kalthoff, F., Stuetz, A., et al. (2013). NVP-LDE225, a potent and selective SMOOTHENED antagonist reduces melanoma growth in vitro and in vivo. PLoS One 8:e69064. doi: 10.1371/journal.pone.0069064
Jiang, J., and Hui, C. C. (2008). Hedgehog signaling in development and cancer. Dev. Cell 15, 801–812. doi: 10.1016/j.devcel.2008.11.010
Kool, M., Korshunov, A., Remke, M., Jones, D. T., Schlanstein, M., Northcott, P. A., et al. (2012). Molecular subgroups of medulloblastoma: an international meta-analysis of transcriptome, genetic aberrations, and clinical data of WNT, SHH, Group 3, and Group 4 medulloblastomas. Acta Neuropathol. 123, 473–484. doi: 10.1007/s00401-012-0958-8
Kool, M., Koster, J., Bunt, J., Hasselt, N. E., Lakeman, A., van Sluis, P., et al. (2008). Integrated genomics identifies five medulloblastoma subtypes with distinct genetic profiles, pathway signatures and clinicopathological features. PLoS One 3:e3088. doi: 10.1371/journal.pone.0003088
Korshunov, A., Remke, M., Kool, M., Hielscher, T., Northcott, P. A., Williamson, D., et al. (2012). Biological and clinical heterogeneity of MYCN-amplified medulloblastoma. Acta Neuropathol. 123, 515–527. doi: 10.1007/s00401-011-0918-8
Kumar, V., Chaudhary, A. K., Dong, Y., Zhong, H. A., Mondal, G., Lin, F., et al. (2017). Design, synthesis and biological evaluation of novel hedgehog inhibitors for treating pancreatic cancer. Sci. Rep. 7:1665. doi: 10.1038/s41598-017-01942-7
Miyazaki, Y., Matsubara, S., Ding, Q., Tsukasa, K., Yoshimitsu, M., Kosai, K., et al. (2016). Efficient elimination of pancreatic cancer stem cells by hedgehog/GLI inhibitor GANT61 in combination with mTOR inhibition. Mol. Cancer 15:49. doi: 10.1186/s12943-016-0534-2
Organisation for Economic Co-operation Development [OECD] (2008). Guideline 425: Acute Oral Toxicity-Up and Down Procedure. Available at: https://www.oecd-ilibrary.org/environment/test-no-425-acute-oral-toxicity-up-and-down-procedure_9789264071049-en
Pan, S., Wu, X., Jiang, J., Gao, W., Wan, Y., Cheng, D., et al. (2010). Discovery of NVP-LDE225, a potent and selective smoothened antagonist. ACS Med. Chem. Lett. 1, 130–134. doi: 10.1021/ml1000307
Park, K. S., Martelotto, L. G., Peifer, M., Sos, M. L., Karnezis, A. N., Mahjoub, M. R., et al. (2011). A crucial requirement for Hedgehog signaling in small cell lung cancer. Nat. Med. 17, 1504–1508. doi: 10.1038/nm.2473
Peukert, S., He, F., Dai, M., Zhang, R., Sun, Y., Miller-Moslin, K., et al. (2013). Discovery of NVP-LEQ506, a second-generation inhibitor of smoothened. ChemMedChem. 8, 1261–1265. doi: 10.1002/cmdc.201300217
Pomeroy, S. L., Tamayo, P., Gaasenbeek, M., Sturla, L. M., Angelo, M., McLaughlin, M. E., et al. (2002). Prediction of central nervous system embryonal tumour outcome based on gene expression. Nature 415, 436–442. doi: 10.1038/415436a
Rajurkar, M., De Jesus-Monge, W. E., Driscoll, D. R., Appleman, V. A., Huang, H., Cotton, J. L., et al. (2012). The activity of Gli transcription factors is essential for Kras-induced pancreatic tumorigenesis. Proc. Natl. Acad. Sci. U.S.A. 109, E1038–E1047. doi: 10.1073/pnas.1114168109
Rimkus, T. K., Carpenter, R. L., Qasem, S., Chan, M., and Lo, H. W. (2016). Targeting the sonic hedgehog signaling pathway: review of smoothened and GLI inhibitors. Cancers 8:E22. doi: 10.3390/cancers8020022
Riobo, N. A., Lu, K., Ai, X. B., Haines, G. M., and Emerson, C. P. (2006). Phosphoinositide 3-kinase and Akt are essential for Sonic Hedgehog signaling. Proc. Natl. Acad. Sci. U.S.A. 103, 4505–4510. doi: 10.1073/pnas.0504337103
Robarge, K. D., Brunton, S. A., Castanedo, G. M., Cui, Y., Dina, M. S., Goldsmith, R., et al. (2009). GDC-0449-a potent inhibitor of the hedgehog pathway. Bioorg. Med. Chem. Lett. 19, 5576–5581. doi: 10.1016/j.bmcl.2009.08.049
Romer, J. T., Kimura, H., Magdaleno, S., Sasai, K., Fuller, C., Baines, H., et al. (2004). Suppression of the Shh pathway using a small molecule inhibitor eliminates medulloblastoma in Ptc1(+/-)p53(-/-) mice. Cancer Cell 6, 229–240. doi: 10.1016/j.ccr.2004.08.019
Roudaut, H., Traiffort, E., Gorojankina, T., Vincent, L., Faure, H., Schoenfelder, A., et al. (2011). Identification and mechanism of action of the acylguanidine MRT-83, a novel potent Smoothened antagonist. Mol. Pharmacol. 79, 453–460. doi: 10.1124/mol.110.069708
Ruiz, I., and Altaba, A. (1997). Catching a Gli-mpse of Hedgehog. Cell 90, 193–196. doi: 10.1016/S0092-8674(00)80325-6
Scales, S. J., and de Sauvage, F. J. (2009). Mechanisms of Hedgehog pathway activation in cancer and implications for therapy. Trends Pharmacol. Sci. 30, 303–312. doi: 10.1016/j.tips.2009.03.007
Schaefer, G. I., Perez, J. R., Duvall, J. R., Stanton, B. Z., Shamji, A. F., and Schreiber, S. L. (2013). Discovery of small-molecule modulators of the sonic hedgehog pathway. J. Am. Chem. Soc. 135, 9675–9680. doi: 10.1021/ja400034k
Sharpe, H. J., Wang, W., Hannoush, R. N., and de Sauvage, F. J. (2015). Regulation of the oncoprotein Smoothened by small molecules. Nat. Chem. Biol. 11, 246–255. doi: 10.1038/nchembio.1776
Solinas, A., Faure, H., Roudaut, H., Traiffort, E., Schoenfelder, A., Mann, A., et al. (2012). Acylthiourea, acylurea, and acylguanidine derivatives with potent hedgehog inhibiting activity. J. Med. Chem. 55, 1559–1571. doi: 10.1021/jm2013369
Song, L., Wang, W., Liu, D., Zhao, Y., He, J., Wang, X., et al. (2016). Targeting of sonic hedgehog-Gli signaling: a potential therapeutic target for patients with breast cancer. Oncol. Lett. 12, 1027–1033. doi: 10.3892/ol.2016.4722
Stecca, B., and Ruiz, I. A. A. (2010). Context-dependent regulation of the GLI code in cancer by Hedgehog and non-Hedgehog signals. J. Mol. Cell Biol. 2, 84–95. doi: 10.1093/jmcb/mjp052
Tan, I. L., Wojcinski, A., Rallapalli, H., Lao, Z., Sanghrajka, R. M., Stephen, D., et al. (2018). Lateral cerebellum is preferentially sensitive to high sonic hedgehog signaling and medulloblastoma formation. Proc. Natl. Acad. Sci. U.S.A 115, 3392–3397. doi: 10.1073/pnas.1717815115
Tang, C. T., Lin, X. L., Wu, S., Liang, Q., Yang, L., Gao, Y. J., et al. (2018). NOX4-driven ROS formation regulates proliferation and apoptosis of gastric cancer cells through the GLI1 pathway. Cell Signal. 46, 52–63. doi: 10.1016/j.cellsig.2018.02.007
Taylor, M. D., Northcott, P. A., Korshunov, A., Remke, M., Cho, Y. J., Clifford, S. C., et al. (2012). Molecular subgroups of medulloblastoma: the current consensus. Acta Neuropathol. 123, 465–472. doi: 10.1007/s00401-011-0922-z
Tremblay, M. R., Lescarbeau, A., Grogan, M. J., Tan, E., Lin, G., Austad, B. C., et al. (2009). Discovery of a potent and orally active hedgehog pathway antagonist (IPI-926). J. Med. Chem. 52, 4400–4418. doi: 10.1021/jm900305z
Wang, Y., Ding, Q., Yen, C. J., Xia, W., Izzo, J. G., Lang, J. Y., et al. (2012). The crosstalk of mTOR/S6K1 and Hedgehog pathways. Cancer Cell 21, 374–387. doi: 10.1016/j.ccr.2011.12.028
Wetmore, C., Eberhart, D. E., and Curran, T. (2001). Loss of p53 but not ARF accelerates medulloblastoma in mice heterozygous for patched. Cancer Res. 61, 513–516.
Xin, M., Zhang, L., Tang, F., Tu, C., Wen, J., Zhao, X., et al. (2014). Design, synthesis, and evaluation of pyrrolo[2,1-f][1,2,4]triazine derivatives as novel hedgehog signaling pathway inhibitors. Bioorg. Med. Chem. 22, 1429–1440. doi: 10.1016/j.bmc.2013.12.055
Yauch, R. L., Dijkgraaf, G. J., Alicke, B., Januario, T., Ahn, C. P., Holcomb, T., et al. (2009). Smoothened mutation confers resistance to a Hedgehog pathway inhibitor in medulloblastoma. Science 326, 572–574. doi: 10.1126/science.1179386
Zhang, X., Zhao, F., Wu, Y., Yang, J., Han, G. W., Zhao, S., et al. (2017). Crystal structure of a multi-domain human smoothened receptor in complex with a super stabilizing ligand. Nat. Commun. 8:15383. doi: 10.1038/ncomms15383
Keywords: L-4, hedgehog, Smo, medulloblastoma, D473H
Citation: Zhu M, Wang H, Wang C, Fang Y, Zhu T, Zhao W, Dong X and Zhang X (2019) L-4, a Well-Tolerated and Orally Active Inhibitor of Hedgehog Pathway, Exhibited Potent Anti-tumor Effects Against Medulloblastoma in vitro and in vivo. Front. Pharmacol. 10:89. doi: 10.3389/fphar.2019.00089
Received: 05 September 2018; Accepted: 23 January 2019;
Published: 21 February 2019.
Edited by:
Carmen Alvarez-Lorenzo, University of Santiago de Compostela, SpainReviewed by:
Violaine See, University of Liverpool, United KingdomEugenio Gaudio, Istituto Oncologico della Svizzera Italiana, Switzerland
Copyright © 2019 Zhu, Wang, Wang, Fang, Zhu, Zhao, Dong and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xiongwen Zhang, eHd6aGFuZ0BzYXQuZWNudS5lZHUuY24= Xiaochun Dong, eGNkb25nQGZ1ZGFuLmVkdS5jbg==
†These authors have contributed equally to this work