- 1State Key Laboratory of Biotherapy and Cancer Center, Sichuan University and Collaborative Innovation Center for Biotherapy, West China Hospital, Sichuan University, Chengdu, China
- 2College of Pharmacy, National and Local Joint Engineering Research Center of Targeted and Innovative Therapeutics, Chongqing Key Laboratory of Kinase Modulators as Innovative Medicine, Chongqing University of Arts and Sciences, Chongqing, China
- 3Department of Clinical Pharmacy, Sichuan Cancer Hospital and Institute, Sichuan Cancer Center, School of Medicine, University of Electronic Science and Technology of China, Chengdu, China
- 4Research Laboratory of Emergency Medicine, Department of Emergency Medicine, West China Hospital, Sichuan University, Chengdu, China
Gastric cancer (GC) is one of the most malignant cancers and is estimated to be fifth in incidence ratio and the third leading cause of cancer death worldwide. Despite advances in GC treatment, poor prognosis and low survival rate necessitate the development of novel treatment options. Fibroblast growth factor receptors (FGFRs) have been suggested to be potential targets for GC treatment. In this study, we report a novel selective FGFR inhibitor, RK-019, with a pyrido [1, 2-a] pyrimidinone skeleton. In vitro, RK-019 showed excellent FGFR1-4 inhibitory activities and strong anti-proliferative effects against FGFR2-amplification (FGFR2-amp) GC cells, including SNU-16 and KATO III cells. Treatment with RK-019 suppressed phosphorylation of FGFR and its downstream pathway proteins, such as FRS2, PLCγ, AKT, and Erk, resulting in cell cycle arrest and induction of apoptosis. Furthermore, daily oral administration of RK-019 could attenuate tumor xenograft growth with no adverse effects. Here, we reported a novel specific FGFR inhibitor, RK-019, with potent anti-FGFR2-amp GC activity both in vitro and in vivo.
Introduction
Gastric cancer (GC) is one of the most malignant cancers and is estimated to be fifth in incidence ratio and the third leading cause of cancer death worldwide (Sung et al., 2021). Patients with early-stage GC are mostly asymptomatic, and most diagnoses are made at the medium or advanced stage (Takahashi et al., 2013). However, at these stages, chemotherapy, radiotherapy, and surgical treatment efficacies are very low (Van Cutsem et al., 2016). Hence, developing new therapeutic strategies with high efficacy and low toxicity for the treatment of advanced GC is crucial.
Fibroblast growth factor receptors (FGFR1-4) belong to the receptor tyrosine kinase family that bind to numerous fibroblast growth factor (FGF) members (Ornitz and Itoh, 2001; Beenken and Mohammadi, 2009; Belov and Mohammadi, 2013). Extracellular ligand FGFs bind to FGFRs and induce FGFR dimerization and then initiate downstream intracellular signaling cascade pathways, including those involving Ras-Raf-Erk, PI3K-AKT-mTOR, and PLCγ/Ca2+ (Sarabipour and Hristova, 2016). Moreover, the JAK2-STAT pathway can be activated by FGFRs in certain cellular contexts (Deo et al., 2002; Cerliani et al., 2011; Li et al., 2019). Furthermore, by regulating signaling cascade pathways, FGFRs mediate physiological processes, such as development progress, cell proliferation, differentiation, and angiogenesis (Eswarakumar et al., 2005; Fukumoto, 2008; Brooks et al., 2012). Anomalously activated FGFR signaling due to FGFR mutations, amplifications, and translocations is involved in tumorigenesis and progression of cancers, including breast cancer, prostate cancer, GC, urothelial cancer, and cholangiocarcinoma (Grose and Dickson, 2005; Katoh, 2010; Lamont et al., 2011; Jain et al., 2018).
Studies indicated that FGFR1 overexpression facilitates peritoneal diffusion via epithelial-to-mesenchymal transition (EMT) in GC (Shimizu et al., 2018). Overexpression of FGFR1 in GC tissue samples was correlated with EphA4 protein expression whose synergy promoted GC development (Oki et al., 2008). In addition, FGFR2 is considered to play an important role in GC, especially in GC with chromosomal instability (Network, 2014). FGFR2 overexpression was identified in 60% of patients with GC, and FGFR2-amplification (FGFR2-amp) was found to occur in approximately 2–15% of these patients (Matsumoto et al., 2012; Su et al., 2014; Han et al., 2015; Ahn et al., 2016; Jia et al., 2016; Kim et al., 2019b). FGFR2 expression level, amplification, and mutations are associated with drug resistance and prognosis (Matsumoto et al., 2012; Tokunaga et al., 2016; Hosoda et al., 2017). Furthermore, FGFR2-amp is an adverse prognostic factor in GC patients. Compared to FGFR2-unamplified GC, FGFR2-amp is significantly associated with lymph node metastasis and worse survival rate (Kim et al., 2019a). Besides, FGFR3 and FGFR4 expression levels and mutations are associated with poor prognosis and drug resistance, inhibition of cell proliferation signals, and induction of apoptosis (Jang et al., 2001; Ye et al., 2012; Ye et al., 2013; Piro et al., 2016). Evidence suggests that FGFR knockdown or inhibition selectively inhibits GC cell line growth (Katoh, 2010; Jang et al., 2017; Zhang et al., 2019). Hence, pharmacological targeting of FGF/FGFRs signaling pathways may be effective for the treatment of FGFRs-altered GC.
Based on biochemistry and cell-based screening from our laboratory’s bioactive compound library, we found a novel FGFR inhibitor RK-019. In this study, we evaluated the anti-neoplastic activity of RK-019 both in vitro and in vivo. RK-019 showed excellent inhibition and great selectivity against FGFR family kinases, strong anti-proliferative effects, and anti-metastasis on FGFR2-amp GC cells. Therefore, RK-019 might be an effective treatment agent against FGFR2-amp GC, making it a candidate drug for GC treatment.
Materials and methods
Materials
RK-019 was synthesized in our laboratory (State Key Laboratory of Biotherapy, Sichuan University, Chengdu, China). RK-019 was dissolved in dimethyl sulfoxide (DMSO) from Sigma (D8418, St Louis, Mo, United States) at a stock concentration of 40 mM and stored in a -80°C refrigerator. The working solution was diluted to 100 nM, 10 nM, 5 nM, 1 nM, and 0.1 nM by DMSO and stored in a -20°C refrigerator. 3-(4,5-Dimethyl-2-thiazolyl)-2,5-diphenyl-2-H-tetra-zolium bromide (MTT) was purchased from CSNpharm (CSN12440, Shanghai, China). A PE Annexin V Apoptosis Detection Kit was purchased from BD Biosciences (559763, Franklin, NJ, United States). Antibodies: FGFR2 (23328), p-FRS2 (3864), PLCγ (5690), p-PLCγ (14008), AKT (4691), p-AKT (4060), Erk (4695), p-Erk (4370), CDK 2 (18048), CDK 4 (12790), CDK 6 (13331), Cyclin D1 (55506), Cyclin E (4136), p27 (3686), Caspase 3 (9662), and Cleaved-caspase 3 (9661) were purchased from Cell Signaling Technology (Beverly, MA, United States). p-FGFR2T653/T654 (AF8210), JAK2 (AF6022), p-JAK2T1007 (AF3022), STAT3 (AF6294), and p-STAT3T705 (AF3293) were purchased from Affinity Biosciences (Changzhou, China). MMP-2 (CY7164) and MMP-9 (CY5205) were purchased from Abways (Shanghai, China). FRS2 (R26776) and β-actin (200068-8F10) are from Zen-bio (Chengdu, China).
Synthesis of RK-019
Synthesis route diagrammed in Figure 1C
Step i: (E)-5-(((5-bromopyridin-2-yl)imino)methyl)-2,2-dimethyl-1,3-dioxane-4,6-dione (2). A mixture of triethyl orthoformate (8.5 g, 80 mmol) and 2,2-dimethyl-1,3-dioxane-4,6-dione (11.5 g, 80 mmol) was heated at 60°C for 2 h. Then, a solution of 5-bromopyridin-2-amine (1, 13.8 g, 80 mmol) in EtOH (80 ml) was added slowly and the resulting reaction mixture was stirred at 60°C for additional 2 h. Upon completion of the reaction, the mixture was cooled to room temperature. The precipitate was filtered and the filter cake was washed with a small amount of EtOH and dried in a vacuum oven to afford 2 as white solid (16.6g, 64%). 1H NMR (400 MHz, CDCl3): δ 11.39–11.24 (m, 1 H), 9.38–9.28 (m, 1 H), 8.47 (d, J = 2.4 Hz, 1H), 7.86 (dd, J = 8.5, 2.4 Hz, 1H), 6.95 (d, J = 8.5 Hz, 1H), and 1.76 (s, 6H). ESI-MS: C12H11BrN2O4, MS (ESI) m/z 327.1 [M + H]+.
FIGURE 1. Development of RK-019. (A) Chemical structure of RK-019. (B) IC50 values of RK-019 against kinases FGFR 1–4. (C) Synthesis route of RK-019. Reagents and conditions: (i) Triethyl orthoformate, 2,2-dimethyl-1,3-dioxane-4,6-dione, EtOH, reflux; (ii) Diphenyl oxide, 220°C; (iii) NIS, DMF, 80°C; (iv) Boronic acids or boronate esters, Pd (dppf)2Cl2, K2CO3, dioxane, H2O, 100°C; (v) 3,5-Dimethoxyaniline, Pd2 (dba)3, 2,2′-bis(diphenylphosphino)-1,1′-dinaphthalene, Cs2CO3, toluene, 100°C; (vi) (a) 2-(tert-Butyldimethylsilyloxy)bromoethane, NaH, DMF, 5°C–r.t.; (b) TBAF, THF, r. t.; and (vii) (a) MsCl, Et3N, CH2Cl2, 0°C; (b) (R)-1-(aminomethyl)ethanol, CH3CN, 100°C. (D) Docking results of RK-019 on FGFR2 protein (PDB ID: 6AGX).
Step ii: 7-bromo-4H-pyrido [1,2-a]pyrimidin-4-one (3). Ph2O (200 ml) was heated to 220°C, then 2 (16.6 g, 51 mmol) was slowly added into the solution. The mixture was stirred at 220°C for 30 min. TLC detected that the reaction was completed. The mixture was cooled and purified by column chromatography with petroleum ether/ethyl acetate (4:1) to afford 3 (10.7g, 93%). 1H NMR (400 MHz, DMSO-d6): δ 9.03 (d, J = 2.2 Hz, 1H), 8.32 (d, J = 6.4 Hz, 1H), 8.08 (dd, J = 9.4, 2.2 Hz, 1H), 7.65 (d, J = 9.4 Hz, 1H), and 6.46 (d, J = 6.4 Hz, 1H). ESI-MS: C8H5BrN2O, MS (ESI) m/z 225.2 [M + H]+.
Step iii: 7-bromo-3-iodo-4H-pyrido [1,2-a]pyrimidin-4-one (4). To a stirred solution of 3 (10.2 g, 45 mmol) in DMF (50 ml) was added NIS (13.6g, 60 mmol), the mixture was stirred at 80°C for 5 h. Upon completion of the reaction, the reaction mixture was cooled to room temperature, and then added with H2O (50 ml) under stirring. The precipitate was collected by filtration, washed with water dried to a constant weight to afford 4 (15.0 g, 95%).1H NMR (400 MHz, DMSO-d6): δ 9.02 (d, J = 2.0 Hz, 1H), 8.75 (s, 1H), 8.13 (dd, J = 9.4, 2.0 Hz, 1H), 7.67 (d, J = 9.4 Hz, 1H). MS (ESI) m/z 351.0 [M + H]+.
Step iv: 7-bromo-3-(1-methyl-1H-pyrazol-4-yl)-4H-pyrido [1,2-a]pyrimidin-4-one (5). 4 (14.1 g, 40 mmol) and 1-methyl-4-(4,4,5,5-tetramethyl-1,3,2-dioxaborolan-2-yl)-1H-pyrazole (8.3 g, 40 mmol), Na2CO3 (8.5 g, 80 mmol) were dissolved in dioxane (100 ml) and H2O (25 ml). The suspension was degassed under nitrogen bubbling for 10 min before Pd (dppf)2Cl2 (2.9 g, 4 mmol) was added. The reaction mixture was heated to 100°C for 5 h, and then diluted with ethyl acetate. The solution was washed with water and brine successively, dried over anhydrous sodium sulfate, and concentrated in vacuum. The residue was purified via silica gel chromatography with petroleum ether/ethyl acetate (2:1) to afford 5 (8.1 g, 66%). 1H NMR (400 MHz, DMSO-d6): δ 8.79–8.68 (m, 2H), 8.32 (s, 1H), 8.10 (s, 1H), 7.78 (dd, J = 9.6, 2.6 Hz, 1H), 7.69 (d, J = 9.6 Hz, 1H).
Step v: 7-((3,5-dimethoxyphenyl)amino)-3-(1-methyl-1H-pyrazol-4-yl)-4H-pyrido [1,2-a]pyrimidin-4-one (6). A mixture of 5 (8.1 g, 26.3 mmol), 3,5-dimethoxyaniline (4.8 g, 31.6 mmol), Pd2 (dba)3 (2.4 g, 2.63 mmol), (±)-BINAP (2.5 g, 3.95 mmol), and cesium carbonate (12.9 g, 39.5 mmol) in anhydrous toluene (120 ml) was degassed with N2 for 10 min. The reaction was heated to 100°C overnight under N2, and then cooled to room temperature. The mixture was filtrated through a Celite and the filter cake was washed with DCM/MeOH (10/1). After concentration of the filtrate, the residue was purified by column chromatography with methyl allylchloride/methanol (20:1) to give 6 (7.1 g, 56%).1H NMR (400 MHz, DMSO-d6): δ 10.70 (s, 1H), 8.67 (d, J = 2.6 Hz, 1H), 8.48 (s, 1H), 8.23 (s, 1H), 7.89 (s, 1H), 7.36 (d, J = 9.6 Hz, 1H), 7.25 (dd, J = 9.6, 2.6 Hz, 1H), 6.28–6.16 (m, 3H), 3.86 (s, 3H), and 3.74 (s, 6H).
Step vi: 7-((3,5-dimethoxyphenyl) (2-hydroxyethyl)amino)-3-(1-methyl-1H-pyrazol-4-yl)-4H-pyrido [1,2-a]pyrimidin-4-one (7). To a solution of 6 (4.0 g, 10.6 mmol) in DMF (60 ml) was added sodium hydride (60% dispersion in mineral oil, 933 mg, 23.32 mmol) slowly at 0 °C under argon. After being stirred at this temperature for 30 min, (2-bromoethoxy) (tert-butyl)dimethylsilane (4.6 ml, 21.2 mmol) was added and then the reaction mixture was stirred at room temperature overnight. Following this time, the reaction mixture was diluted with water and then extracted with ethyl acetate twice. The combined extracts were washed with water and brine, dried over anhydrous Na2SO4, and concentrated under reduced pressure. The residue was dissolved in tetrahydrofuran (20 ml) and tetrabutylammonium fluoride (1 M in tetrahydrofuran, 21.2 ml, 21.2 mmol) was added. The reaction mixture was stirred at room temperature overnight. Following this time, the reaction mixture was concentrated under vacuum and the resulting residue was purified by column chromatography with methyl allylchloride/methanol (10:1) to give the title compound 7 as yellow solid (2.5g, 56%). 1H NMR (400 MHz, DMSO-d6): δ 8.77–8.71 (m, 2H), 8.38 (s, 1H), 8.12 (s, 1H), 7.68 (dd, J = 9.6, 2.6 Hz, 1H), 7.61 (d, J = 9.6 Hz, 1H), 6.43–6.32 (m, 3H), 4.96 (t, J = 5.0 Hz, 1H), 3.93–3.83 (m, 5H), 3.71 (s, 6H), and 3.68–3.61 (m, 2H). HRMS: calculated for C22H23N5O4 [(M + H)+], 422.1824; found 422.1821.
Step vii: (R)-7-((3,5-dimethoxyphenyl) (2-((1-hydroxypropan-2-yl)amino)ethyl)amino)-3-(1-methyl-1H-pyrazol-4-yl)-4H-pyrido [1,2-a]pyrimidin-4-one (RK-019). To a mixture of 7 (2.5 g, 5.9 mmol) and Et3N (2.1 ml, 14.75 mmol) in DCM (30 ml) was added methanesulfonyl chloride (0.917 ml, 11.8 mmol) dropwise at 0°C under argon. The reaction mixture was stirred at room temperature for 2 h. After completion (monitored by TLC), the reaction mixture was quenched with water, then extracted with DCM twice. The combined extracts were washed with water and brine successively, dried over anhydrous Na2SO4, and concentrated under vacuum. The residue was purified by column chromatography to give the intermediate (2.2 g). A mixture of the intermediate (2.2 g, 4.4 mmol) and (R)-1-(aminomethyl)ethanol (3.3 g, 44 mmol) in acetonitrile (44 ml) was heated at 100°C overnight. After completion (monitored by TLC), the reaction mixture was concentrated under vacuum and the residue was purified on TLC-preparative plates with methyl allylchloride/methanol (10:1) to afford the desired product RK-019 as light-yellow solid (366 mg, 18%), mp.144–150°C. 1H NMR (400 MHz, DMSO-d6): δ 8.73 (s, 1H), 8.66 (d, J = 2.5 Hz, 1H), 8.38 (s, 1H), 8.12 (s, 1H), 7.66 (dd, J = 9.6, 2.5 Hz, 1H), 7.61 (d, J = 9.6 Hz, 1H), 6.38 (d, J = 2.1 Hz, 2H), 6.31 (d, J = 2.1 Hz, 1H), 4.51–4.43 (m, 1H), 3.89 (s, 3H), 3.84 (t, J = 6.5 Hz, 2H), 3.72 (s, 6H), 3.29–3.20 (m, 3H), 2.86–2.75 (m, 2H), 2.63–2.57 (m, 1H), and 0.90 (d, J = 6.2 Hz, 3H). 13C NMR (101 MHz, DMSO-d6) δ 161.88, 149.53, 148.44, 147.50, 143.96, 140.08, 138.25, 135.93, 131.16, 129.47, 121.78, 120.85, 109.54, 104.32, 97.49, 55.76, 52.78, 48.55, 44.28, and 23.31. HRMS: calculated for C25H30N6O4 [(M + H) +], 479.2403; found 479.2405.
Molecular docking
The 3D structure of FGFR2 was downloaded from the PDB (http://www.rcsb.org/, PDB ID: FGFR1: 5EW8, FGFR2: 6AGX, FGFR3: 6LVM, and FGFR4: 6NVK). The protein was prepared with discovery studio 3.1. Molecule RK-019 was built with ChemBio3D and optimized at the molecular mechanical level. Then, RK-019 was docked to the binding site of JNJ42756493 by employing a protein-ligand docking program GOLD 2.5, respectively. Scoring function GOLDSCORE was used for exhaustive searching, solid body optimizing, and interaction scoring. The final results for molecular docking were visualized by using the PyMol program.
Kinase inhibition and selectivity assay
Median-inhibitory concentration (IC50) and inhibition ratio of kinases by RK-019 was assessed by Eurofins Discovery Services (Dundee, United Kingdom) using the ATP-site competition binding assay in vitro. Firstly, gradient concentration of RK-019 against FGFRs, including FGFR1, FGFR2, FGFR3, and FGFR4, were tested. IC50 value was calculate by the GraphPad Prism software (v8.4.3, GraphPad Software, California, United States). Then, to measure the selectivity prolife of RK-019, a panel of 422 recombinant human kinases was screened at a concentration of 1 µM RK-019. Kinase MAP was drawn by webtools from Cell Signaling Technology, Inc. (www.cellsignal.com).
Cell lines
SNU-16, KATO III, AZ 521, MGC 80–3, HGC 27, N87, AGS, NUGC-4, GT 39, MKN 45, BGC823, and GES-1 were purchased from BeNa Culture Collection (Beijing, China). HGC27, AGS, NUGC-4, and GT39 cell lines were cultured in DMEM (L110KJ, BasalMedia, Shanghai, China) with 10% of fetal bovine serum (900–108, Gemini, California, USA). SNU-16, AZ521, MGC 80–3, N87, MKN45, BGC823, and GES-1cell lines were cultured in RPMI-1640 medium (L210KJ, BasalMedia) with 10% of fetal bovine serum. KATO III cell was cultured in IMDM medium (L610KJ, BasalMedia) with 20% of fetal bovine serum. All the culture medium was contained with 1% antibiotic (30–002-CI, Corning, Corning, NY, United States). All the cell lines were maintained in 5% CO2 condition at 37 °C.
Cell viability assay
For suspension and mixed suspension adherent cells, the cells were seeded in 96-well plate with different cell density, 0.5×104 per well for 96 h treatment, 1×104 per well for 72 h treatment, 2×104 per well for 48 h treatment, and 5×104 per well for 24 h treatment. Different doses of RK-019 were treated as soon as the cell seeded in 96-well plates. For adherent cells, the cells were seeded in a 96-well plate with 3×103 per well and cultured overnight followed by administrated different doses of RK-019. After indicated time treatment, MTT solution at a final concentration of 0.5 mg/ml was added and incubated 2–4 h at 37°C. Formazan formed by the living cells was dissolved with DMSO and the absorbance was measured using a microplate spectrophotometer (MultiskanFC, ThermoFisher, Waltham, MA, United States) at 570 nm. IC50 value of RK-019 was calculated by the Graphpad Prism software.
Quantitative reverse transcription PCR (qRT-PCR)
The GC cell total RNAs were extracted by TRIzol (15596026, ThermoFisher). RNA was reversed transcribed into cDNA by using HiScript III RT SuperMix for qPCR (R323-01, Vazyme, Nanjing, China) according to the manufacturer’s protocol. qRT-PCR for FGFR2 gene expression was carried out by ChamQ universal SYBR qPCR Master Mix (Q711-02, Vazyme) on Bio-Rad CFX96 Realtime PCR system. The PCR primer for FGFR1: forward 5‘- GCTACAAGGTCCGTTATGC -3’ and reverse 5′- CAATGCAGGTGTAGTTGCC -3’, FGFR2: forward 5‘- GGTGGCTGAAAAACGGGAAG -3’ and reverse 5′- AGATGGGACCACACTTTCCATA -3’, FGFR3: forward 5‘- TGCGTCGTGGAGAACAAGTTT -3’ and reverse 5′- GCACGGTAACGTAGGGTGTG -3’, FGFR4: forward 5‘- CCATAGGGACCCCTCGAATAG -3’ and reverse 5′- CAGCGGAACTTGACGGTGT -3’, and ACTB: forward 5‘- CACCATTGGCAATGAGCGGTTC -3’ and reverse 5′- AGGTCTTTGCGGATGTCCACGT -3’. The gene expression analysis was conducted by the ΔΔCq method.
Cell counting assay
SNU-16 and KATO III were seeded in 6-well plated with 1×105 per well. RK-019 was treated as soon as the cells seeded in the 6-well plates. After indicated time, the images were captured using an inverted microscope. Then, the cells in the captured well were resuspended and live cell numbers were determined by the trypan blue (BL627A, Biosharp, Hefei, China) method.
Edu staining assay
For Edu staining assay, 1×105 cells were seeded in 24-well plates. RK-019 and DMSO were treated as soon as the cells seeded into the plates. After 24 h, the cells were incubated with the medium contained 10 µM Edu for 2 h. After that, the cells were collected and fixated to slide, and then stained with the Cell-Light EdU Apollo488 In Vitro Kit (C10310, Ribobio, Guangzhou, China). After staining, the samples were photographed and analyzed by the ImageJ software (v1.53a, Wayne Rasband, United States).
Cell cycle and apoptosis analysis
Cell cycle and apoptosis assay were both measured by flow cytometry (FCM). The data were analyzed by the NovoExpress software (v1.4.0, Agilent, Santa Clara, CA, United States). For cell cycle assay, the cells were treated by different concentrations of RK-019 for 24 h, then were harvested and washed twice with cold phosphate buffered solution (PBS), followed by fixed with 75% ethanol overnight. Before test, the cells were washed twice with cold PBS and stained with PI stanning solution using cell cycle detecting kit (KGA512, KeyGEN, Nanjing, China), and then measured by FCM. For apoptosis assay, the cells were harvested after treated with different concentrations of RK-019 for 24 h, and stained with a PE Annexin V Apoptosis Detection Kit according to the manufacturer’s protocol followed by detected using FCM.
Transwell migration and invasion assay
In the migration assay, 1×105 cells were resuspended in a serum-free medium, then added into the upper chamber, meanwhile the medium containing 10% FBS was added at the bottom with 100 nM RK-019. After 24 h incubation, the migration cells on the filters were fixed with 4% paraformldehyde and washed with PBS, then stained with crystal violet solution for 20 min. In the invasion assay, the upper surface of the transwell was coated with Matrigel (356234, Corning) for 30 min at 37°C until Matrigel was solidified. 1×105 cells were resuspended in a serum-free medium, then added into the upper chamber, meanwhile the medium containing 10% FBS was added at the bottom with 100 nM RK-019.After 24 h , the invasion cells on the filters were hatched and fixed, then stained with crystal violet.
Western blot analysis
The cells treated with RK-019 at indicated concentration for 24 h were harvested and lysed in 1× RIPA buffer (20e188, MilliporeSigma, Burlington, MA, United States), which contained protease inhibitor cocktail (B14001, Bimake, Houston, TX, United States) and phosphatase inhibitor cocktail (B15001, Bimake), for 30 min and equalized by concentration of the total protein before loading. Using sodium dodecyl sulphate - polyacrylamide gel electrophoresis (SDS-PAGE) to separate the protein sample and transfer to PVDF membranes (ISEQ00010, MilliporeSigma). The membranes were incubated with relevant primary antibody and corresponding HRP-labeled secondary antibody. Then, chemiluminescence was used to detect target bands. ImageJ software were used to make the gradation analysis.
In western blot, results analysis of relative expression quantification level, data were firstly calculated according to the following formular: Protein Relative expression data = Gt/Gi. Gt means the gradation level of target protein band and Gi means the gradation level of internal reference protein band. Then, the relative expression level can be calculated by the following formular: Protein Relative expression level = Dt/Dc ×100%. Dt means the protein relative expression data of target group and Dc means the protein relative expression data of control group.
In western blot, results analysis of relative phosphorylation quantification level, data were firstly calculated according to the following formular: Protein Relative phosphorylation data = Gp/Gup. Gp means the gradation level of target phosphorylation protein band and Gup means the gradation level of target total protein band. Then, the relative phosphorylation level can be calculated by the following formular: Protein Relative expression level = Pt/Pc ×100%. Pt means the relative phosphorylation data or target group and Pc means the relative phosphorylation data of control group.
In vivo PK assay
All animal experiments in this study have been approved by the Institutional Animal Care and Treatment Committee of Sichuan University in China and were carried out in accordance with the approved guidelines. Sprague-Dawley rats (weight 180–200 g, Beijing HJF bioscience, Beijing, China) were used in this experiment and maintained in a specific-pathogen-free (SPF) condition facility. RK-019 was dissolve in ethanol first, followed by mixed with Kolliphor EL, then filled with saline to the calculated volume. The formula was 12.5% ethanol, 12.5% Kolliphor EL, and 75% saline. RK-019 was administrated via i. v. 3 mg/kg or p. o. 30 mg/kg. After 5, 10, 15, 30, 60, 120, 240, 360, 600, and 1,440 min of administration, the blood samples were collected. Plasma fraction was obtained by centrifuged blood, then deproteinized with methanol containing an internal standard. The compound concentrations of the target compound in the supernatant were measured by liquid chromatography–mass spectrometry/mass spectrometry (LC-MS/MS).
Subcutaneous xenograft model
NOD/SCID mice (6-week-old, weight 18–20 g, Beijing HJF bioscience) were used in this experiment and maintained in a specific-pathogen-free (SPF) condition facility. SNU-16 Cells were resuspended in PBS mixed with Matrigel (1:1) and adjust the density to 1×107 per mice (100 μL). The cell suspensions were inoculated subcutaneously on the right flank of the mice. When the mean tumor volume reached approximately 200 mm3, the mice were divided into four groups (five mice each group): i. Vehicle group, ii. RK-019 15 mg/kg group, iii. RK-019 30 mg/kg group, and iv. RK-019 45 mg/kg group. RK-019 was dissolved in same formula as the in vivo PK study. RK-019 was administrated once per day by oral gavage for 21 days. Tumor volumes were measured three times per week. The tumor volume was measured by a caliper and calculated according to the following formula:
At the end of administration, the mice were executed euthanasia, then tumors, blood, and organs were extracted. The complete blood count (CBC), serum phosphorus analysis, and blood biochemical analysis were completed by West China Frontier Pharma Tech (Chengdu, China). The pathological section and IHC stanning was completed by Servicebio (Wuhan, China). The tumor growth inhibition (TGI) values were calculated with the following formula:
Tn and T0 represent average tumor volume before treatment and that of day n after treatment in the treatment group. Cn and C0 represent average tumor volume before treatment and that of day n after treatment in the vehicle group. In this study, n is 21st day.
Statistical analyses
All data were analyzed by the GraphPad Prism 8.4.3 (GraphPad Software, California, United States), and was shown as mean value ±SD or SEM, details will be illustrated in the figure legends. Dose-effect curve analysis (IC50 calculation) was performed by the Graphpad Prism software. First, concentration value was converted to log value, followed by non-linear regression analyze which is log(inhibitor) vs. normalized response—Variable slope, was performed to calculate the IC50 values. The statistically significant p values were calculated by student’s t-test or two-way ANOVA and were shown in graphics which were labeled as follows: *p < 0.05; **p < 0.01; ***p < 0.001.
Result
Development of RK-019
In this study, we identified a novel small molecule pan-FGFR inhibitor, RK-019, whose chemical structure was named as (R)-7-((3,5-dimethoxyphenyl) (2-((1-hydroxypropan-2-yl)amino)ethyl)amino)-3-(1-methyl-1H-pyrazol-4-yl)-4H-pyrido [1,2-a]pyrimidin-4-one (Figure 1A). The method of synthesis of RK-019 is described in detail in the materials and methods section (Figure 1C). Briefly, RK-019 was synthesized from 5-bromopyridin-2-amine 1) with imidization, cyclization, iodination, and Suzuki coupling of 4 with 1-methyl-4-pyrazole boronic acid pinacol ester, followed by Buchwald coupling to get 6. Finally, nucleophilic substitution with 6 and 2-isopropylaminoethylchloride and deprotection by hydrochloric acid produced RK-019. 1H-NMR, 13C-NMR, and HRMS spectrometry results were shown in Supplementary Figure S3.
Molecular docking was conducted to model the binding of RK-019 in the ATP pocket using a reported crystal structure of FGFR kinase domain (PDB ID: 6AGX). As shown in Figure 1D, RK-019 could fit well in the ATP pocket of FGFR2 and the N-methyl pyrazole group could extend into the solvent region. RK-019 maintained H-bond interactions with FGFR2 at the backbone NH of Ala567 in the hinge region. The methoxyphenyl motif occupied the hydrophobic region in the ATP pocket and its methoxy group formed a H-bond with the NH of Asp644. The 2-amino-propanol group can interact with the main carbonyl group of Asp644 and Asn631 to form two hydrogen bonds. Overall, RK-019 could fit well into the ATP pocket of FGFR2. Similarity, RK-019 could also fit well in FGFR1 (PDB ID: 5EW8), FGFR3 (PDB ID: 6LVM), and FGFR4 (PDB ID: 6NVK) ATP binding pocket (Supplementary Figure S1). Then, the ATP-based kinase activity assay provided by kinase profile service (Eurofins Discovery) revealed that RK-019 was sufficient to inhibit the FGFR family kinase activity, with the IC50 of 9.1 nM (FGFR1), 4.6 nM (FGFR2), 26.3 nM (FGFR3), and 40.7 nM (FGFR4) (Figure 1B).
Kinase selectivity profile of RK-019
To further investigate the kinase selectivity, the kinase inhibitory profile of RK-019 was determined against a diverse panel of 422 recombinant human kinases from Eurofins Discovery by the ATP-site competition binding assay at a concentration of 1 µM (Supplementary Table S1). The kinase inhibition results were depicted in a kinase MAP at a 35% activity cutoff of the DMSO control (Figure 2). Of these, only seven hits (FGFR1, FGFR2, FGFR3, FGFR4, FLT4, RET, and LYN) showed over 90% of inhibition. A further enzyme activity assay was performed to determine the IC50 values of these seven hits. As results shown, RK-019 showed a great selectivity to the FGFRs family kinase, which indicated that the molecule was a potent and selective pan-FGFR inhibitor.
FIGURE 2. Kinase selectivity profile of RK-019. Measurement of 422 kinases was performed at 1 μM of RK-019. Each kinase was measured once. Data were cutoff on 35% kinase activity compared to the DMSO group. The TREE spot image was mapped with the KinMap software tool provided by Cell Signaling Technology, Inc. (www.cellsignal.com). The percentage of control means remaining active kinase percentage.
Anti-proliferation effects of RK-019 against FGFR2-amp GC cell lines
Different GC cell lines were used to verify the biological activity of RK-019, including SNU-16, KATO III, AZ 521, MGC 80–3, HGC 27, N87, AGS, NUGC-4, GT 39, MKN 45, BGC823, and gastric epithelial cell line GES-1. The FGFR2 expression level is relevant to poor pathological features and prognostic, and anti-FGFR2 agents could provide potential effectiveness in GC treatment [22, 31, 32]. Firstly, we measured the FGFRs mRNA expression levels in these GC cell lines by qRT-PCR. By comparing the expression data, we have found that FGFR2-amp cell lines, SNU-16, and KATO III exhibited the highest mRNA expression level of FGFR2 (Figure 3A, Supplementary Figure S2). Then, MTT assay was performed to test the viability of GC cell lines following treatment with RK-019. As shown in Figure 3B, SNU-16 and KATO III cell lines were demonstrated to be the most sensitive cell lines to RK-019, with the IC50 values of SNU-16 and KATO III cell lines were 3.96 ± 4.4 nM and 5.45 ± 5.3 nM, respectively.
FIGURE 3. Anti-proliferative action of RK-019. (A) qRT-PCR analysis of the FGFR2 relative expression level on different types of gastric cancer cell lines. GES-1 was used as control. Data were shown in mean ± SEM. (B) IC50 of RK-019 measured by MTT assay on different types of gastric cancer cell lines. Data were shown in mean ± SD. (C) SNU-16 and KATO III cell lines were treated with different concentrations of RK-019 for 24, 48, 72 and 96 h, respectively. The inhibition ratio indexes by MTT assay were the percentage of cells contrast to DMSO-treated group at respect time. Data were shown in mean ± SD. (D) Edu staining was used to detect the cell proliferation phenotype with 100 nM RK-019 treatment for 24 h (E) SNU-16 and KATO III were treated with RK-019 for 9 days in different concentrations, cells were counted by trypan blue and took pictures every 2 day. The picture only shows the 100 nM RK-019 treatment group on Day 0 and Day 9. Data were shown in mean ± SD. All the data were compiled from three independently repeat experiments. The significance was determined by Student’s t-test, *p < 0.05; **p < 0.01; ***p < 0.001, for the designated treatment group vs. the DMSO group.
Furthermore, SNU-16 and KATO III cells were exposed to graded concentrations of RK-019 for 24, 48, 72, and 96 h to evaluate the dose- and time-dependent relationships. As shown in Figure 3C, RK-019 showed limited suppression effect after 24 h treatment and obvious suppression effect after 48 and 72 h, but did not significantly increase after 96 h. Moreover, Edu staining was used to verify the changes in cell proliferation. After 24 h of RK-019 administration, the number of Edu-positive stained cells significantly decreased in SNU-16 and KATO III (Figure 3D). The results demonstrated that the inhibition of cell proliferation by RK-019 occurred in a time- and dose-dependent manner. In addition, a prolonged drug administration experiment was performed. SNU-16 and KATO III were treated with 1, 5, 10, and 100 nM of RK-019 for 9 days. Live cell numbers were determined by trypan blue staining and photos were taken every 2 days to verify the RK-019 anti-proliferation activity. As shown in Figure 3E, SNU-16 cells were inhibited by treatment with 10 and 100 nM of RK-019, and KATO III cells was suppressed by treatment with 5, 10, and 100 nM of RK-019.
RK-019 inhibited FGFR2 phosphorylation and the downstream signal pathway in SNU-16 and KATO III cells
When FGFR is bound to its ligand, auto-phosphorylation at T653/T654 site occurs, which is important for the FGFR kinase activity (Zou et al., 2012). Activated FGFR can phosphorylate downstream proteins, such as FRS2 and PLCγ (Xu et al., 1998; Goetz and Mohammadi, 2013), and further activate PI3K-AKT-mTOR and Ras-Raf-Erk pathways (Cailliau et al., 2001; Kamata et al., 2002; Tsang and Dawid, 2004). FGFR signaling activation could induce cell proliferation and survival through the signaling cascade. Thus, western blot analysis was used to study the molecular mechanisms underlying the anti-tumor effects of RK-019 and phosphorylation levels of FGFR2, and downstream proteins were determined and analyzed in SNU-16 and KATO III cells.
As shown in Figures 4A and B, after treated with 0.1, 1, 5, 10, and 100 nM of RK-019 for 24 h in SNU-16 and KATO III cells, the protein expression level of phosphorylated FGFR2T653/T654 was decreased. Furthermore, the FRS2 and PLCγ phosphorylation levels were also significantly decreased, indicating that RK-019 could suppress the FGFR signaling activity by inhibiting its auto-phosphorylation ability. In addition, we observed that the AKT and Erk phosphorylation levels were also decreased (Figures 4C and D). The AKT and Erk phosphorylation levels reflect the activities of the PI3K-AKT-mTOR and Ras-Raf-Erk pathways, which play an important role in cell proliferation and survival (Lawlor and Alessi, 2001; Zhang and Liu, 2002). Overall, our results suggested that RK-019 could inhibit the auto-phosphorylation of FGFR2, affect the activities of downstream proteins, such as FRS2 and PLCγ, and deactivate the PI3K-AKT-mTOR and Ras-Raf-Erk pathways, resulting in cell proliferation inhibition.
FIGURE 4. RK-019 effectively inhibits the phosphorylation of FGFR2 and the downstream protein in SNU-16 and KATO III cells. SNU-16 and KATO III cells were treated with different concentrations of RK-019 for 24 h. Then, the cells were harvested and lysed with RIPA, then western blot analyses were performed. Western blot results of FGFR signal related proteins in SNU-16 (A) and KATO III (B). Western blot results of downstream of FGFR signal in SNU-16 (C) and KATO III (D). Each experiment was repeated for three times independently. ImageJ software was used to quantify the analysis of image profile. Quantification data are performed in mean value ±SD. The significance was determined by Student′s t-test, *p < 0.05; **p < 0.01; ***p < 0.001, for the designated treatment group vs. the DMSO group.
RK-019 induced G0/G1 phase arrest
The FCM analysis was performed to detect the cell cycle distribution. As shown in Figure 5A and B, RK-019 induced G0/G1 arrest after 24 h treatment of RK-019 in SNU-16 and KATO III cells. The G0/G1 proportion was raised from 24.75% to 26.81%–50.29% and 67.88% in SNU-16 and KATO III cells, respectively, after 100 nM RK-019 administration.
FIGURE 5. RK-019 induced G0/G1 phase arrest. SNU-16 and KATO III cells were treated with different concentrations of RK-019 for 24 h (A), (B) FCM analysis cell cycle distribution by PI staining. Cell cycle distribution was analyzed and quantified by the NovoExpress software. Values were shown by mean ± SD. (C), (D) Western blot of cell cycle related proteins and quantification analysis. ImageJ software was used to quantify the analysis of image profile. Quantification data were performed in mean value ±SD. Each experiment was repeated for three times independently. Significance was determined by Student’s t-test *p < 0.05; **p < 0.01; ***p < 0.001, for the designated treatment group vs. the DMSO group.
We also examined the cell cycle related protein expression levels in cells by western blot analysis. As shown in Figures 5C and D, the G0/G1 cell cycle protein expression levels, including CDK2, 4, and 6, and Cyclin D1 and E, were obviously decreased in KATO III and SNU-16 cells, respectively. Meanwhile, the protein expression level of cell cycle suppression factor p27 was both increased in SNU-16 and KATO III cells. The protein expression levels of CDK2, 4, 6, and Cyclin D1 and E are related with G0/G1 arrest, and p27 can mediate the inhibition of cyclin-CDK2 complex function and result in cell cycle arrest (Massagué, 2004). Therefore, we conclude that RK-019 might induce G0/G1 arrest by inhibited the expression levels of cyclin-CDK complex.
RK-019 induced apoptosis
FGFR signal activation in different types of cancer can lead to the apoptosis resistance (Acevedo et al., 2009). Here, we used the Annexin V-PE/7-AAD dual-labeling method to evaluate the apoptosis levels. SNU-16 and KATO III cells were treated with different concentrations of RK-019 for 72 h and analyzed by FCM. The population of apoptotic cells (early and late apoptotic cells) was increased significantly in SNU-16 (45.5%) and KATO III (32.7%) cells (Figures 6A and B), indicating that apoptosis was induced by RK-019. Further, western blot analysis showed that the amount of cleaved caspase-3 in SNU-16 and KATO III was significantly increased after administering RK-019. In addition, the presence of activated caspase-3 and the cleaved version of poly (ADP-ribose) polymerase (PARP) verified the induction of apoptosis (Figures 6C,D).
FIGURE 6. RK-019 induced cell apoptosis. SNU-16 and KATO III cells were treated with different concentrations of RK-019 for 72 h (A), (B) FCM analysis cell apoptosis by 7-AAD and Annexin V-PE duel staining and quantified FCM apoptosis analysis data. Values were shown the total number of Q2-2 (early apoptotic) and Q2-4 (late apoptotic) by mean ± SD. (C), (D) Western blot of apoptosis related proteins and quantification analysis, data were shown by mean ± SD. Each experiment was repeated for three times independently. Significance was determined by Student’s t-test *p < 0.05; **p < 0.01; ***p < 0.001, for the designated treatment group vs. the DMSO group.
RK-019 inhibits cell migration and invasion
Peritoneal metastasis is one of the leading causes of death in patients with GC (Yao et al., 2020) and FGFRs can phosphorylate STAT to promote cancer cell metastasis (Babina and Turner, 2017). As shown in transwell assays, RK-019 could significantly suppress the migration and invasion ability in SNU-16 and KATO III cells (Figures 7A,B). Furthermore, western blot analysis for migration and invasion-related protein confirmed the effect of RK-019 (Figure 7C). The proteins expression levels of MMP-2 and MMP-9 were reduced obviously after RK-019 treatment. The phosphorylation levels of JAK2 and STAT3, which are FGFRs-related downstream proteins, were decreased. As the STAT signaling pathway can regulate the expression levels of MMP-2 and MMP-9 directly (Xie et al., 2004; Jia et al., 2017), we speculated that RK-019 may suppress the migration and invasion of GC cells by blocking the JAK2-STAT signal axis.
FIGURE 7. RK-019 inhibit SNU-16 and KATO III migration and invasion. Migration and invasion assay, cells (1 × 105) were seeded in the top chamber of transwell with serum-free medium and treated with RK-019 100 nM for 24 h, then the cells were stained with crystal violet. (A) Photographs of migration assay. (B) Photographs of invasion assay. (C) Western blot of FGFR downstream protein, including JAK2, p-JAK2, STAT3 and p-STAT3, and migration and invasion proteins, MMP-2 and MMP-9. (D), (E) Quantification analysis by the ImageJ software of western blot results, data were performed by mean ± SD. Each experiment was repeated for three times independently. Significance was determined by Student’s t-test *p < 0.05; **p < 0.01; ***p < 0.001, for the designated treatment group vs. the DMSO group.
In vivo efficacy of RK-019
Pharmacokinetic profile of RK-019 was conducted in Sprague-Dawley rats (i.v.,3 mg/kg; oral, 30 mg/kg). The pharmacokinetic parameters are listed in Table 1. After a single oral administration, RK-019 showed slow absorption (Tmax = 2.67), the peak plasma concentration (Cmax) was 234.65 ng/ml, the area under the plasma concentration time curve (AUC) was 1,448.41 ng h/mL, the biological half-life (T1/2) was 2.83 h, and the oral relative bioavailability was 19.00%.
To investigate the anti-tumor effect of RK-019 in vivo, we established a xenograft model using SNU-16 cells. RK-019 was administered at doses of 15, 30, and 45 mg/kg once a day and after 21 days of treatment, the TGI values were 34.3, 85.9, and 83.5%, respectively (Figure 8A). Meanwhile, a significant decrease in tumor weight in the RK-019-treated group was observed compared to that in the vehicle group (Figure 8B).
FIGURE 8. Anti-tumor efficacy of RK-019 in the SNU-16 xenograft model. NOD/SCID mice bearing SNU-16 were orally treated with vehicle or RK-019 once a day for 21 days. Tumor volume and body weight were measured three times per week. After 21 days of administration, mice were euthanasia and tumor tissue were collected, and were immunohistochemically analyzed with anti-Ki-67, anti-cleaved caspase-3, and anti p-FGFR2 antibodies. (A) Tumor growth curve, data were performed by mean ± SEM. (n = 5) (B) Tumor weight, data were performed by mean ± SD. (n = 5) Significance was determined by two-way ANNOVA *p < 0.05; **p < 0.01; ***p < 0.001, for the designated treatment group vs. the Vehicle group. (C) Mice body weight curve, data were performed by mean ± SD. (n = 5) (D) IHC results of tumor tissue after treatment. (E) Western blot analysis of FGFR downstream protein and migration and invasion proteins, including JAK2, p-JAK2, STAT3, p-STAT3, MMP-2, and MMP-9. Each group was randomly selected three tumor tissues, followed by lysis and western blot analysis was performed. Results were quantified by the ImageJ software, and data was performed by mean ± SD. Significance was determined by Student’s t-test *p < 0.05; **p < 0.01; ***p < 0.001, for the designated treatment group vs. the vehicle group.
To validate the results of the in vivo assay, we conducted IHC analyses using Ki67 as a cell proliferation marker, cleaved caspase-3 as an apoptotic cell marker, and p-FGFR2 as a FGFR2 signaling activity marker. As shown in Figure 8D, the number of p-FGFR2-positive and Ki67-positive cells decreased, while the number of cleaved-caspase 3-positive cells increased significantly. Moreover, western blot analysis was performed to detecting the influence of RK-019 on JAK2-STAT signal axis in vivo (Figure 8D). Although the changes of relative phosphorylation level of JAK2 and STAT3 are not significant, the JAK2 and STAT3 phosphorylation level decreased remarkably in 30 and 45 mg/kg RK-019 treatment groups. The downstream proteins expression of this axis, MMP2 and MMP9, have also decreased after RK-019 treatment. These results suggested that RK-019 could block the FGFR signaling pathway in vivo, and result in proliferation arrest and cell apoptosis induction, corroborating the in vitro observations. In the meantime, RK-019 exhibit the potential ability, suppress the migration, and invasion of GC cells in vivo.
Meanwhile, preliminary toxicity of RK-019 was evaluated. During the treatment, no significant body weight change occurred (Figure 8C). Moreover, no pathological changes were observed in the harvested organs of the RK-019-treated group (Figure 9D), and hematopoietic toxicity was absent (Figures 9A–C) at the end of treatment, suggesting that the mice were tolerant to RK-019 treatment. Taken together, RK-019 inhibited tumor growth in the SNU-16 xenograft model by inhibiting cell proliferation and inducing cell apoptosis.
FIGURE 9. Safety profile of RK-019 in the xenograft mice model. After 21 days of RK-019 treatment, the mice were euthanasia, then blood and organs were extracted and analysis. Values in graphics were plotted as mean ± SD (n = 5). (A) Complete blood count data. (B) Blood biochemical analysis. (C) Serum phosphorus analysis. (D) H&E staining of pathological section in the heart, liver, spleen, lungs, and kidneys.
Discussion
FGFR signaling is an important regulatory pathway. However, FGFR signaling disorder is closely related to many diseases, such as cancer and fibrosis (Robertson et al., 2000). Abnormally activated FGFR signaling in cancer activates a series of signaling pathways, such as PI3K-AKT-mTOR, Ras-Raf-Erk, and JAK-STAT, which results in uncontrolled cell proliferation, metastasis, and avoidance of cell death (Gschwind et al., 2004). Preclinical studies have identified FGFRs as a potential therapeutic target in many types of cancers, including lung, breast, gastric, and hematologic cancers. Excitingly, FGFR inhibitors, including Erdafitinib (JNJ-42756493) and Infigratinib (BGJ-398), have been approved for urothelial carcinoma and cholangiocarcinoma treatment (Markham, 2019; Montazeri and Bellmunt, 2020; Botrus et al., 2021).
GC is a highly heterogeneous disease from morphological and molecular standpoints (Gullo et al., 2018). This results in heterogeneity of treatment effect and makes the GC treatment development more difficult than other types of cancers. Until now, therapeutic schedules for GC are very limited. Traditional chemotherapy drugs, such as anti-tumor platinum drugs, docetaxel, and 5-Fu, still are the most common drugs in gastric cancer treatment (Smyth et al., 2020). Compared to the chemotherapy, targeted therapy exhibited more advantages in efficacy and safety. However, targeted drugs for GC treatment are very limited. Studies have revealed FGFR2 as a potential target for GC treatment (Lengyel et al., 2022). Therefore, discovery targeted therapy for FGFR2 is very significant.
However, FGFRs inhibitors still are face many challenges. Firstly, although regulation of FGFR signaling cascades has been widely investigated, their unique function and drug resistance mechanisms remain unclear (Babina and Turner, 2017; Simons, 2021). Use of FGFRs inhibitors might cause a series of problems such as mutations conferring resistance to FGFR-targeting drugs and side-effects (Yue et al., 2021). On the other hand, due to the patient heterogeneity, FGFRs kinase inhibitors and anti-body drug conjugates had shown lower efficacy than expected, and drug toxicity, multiple GC treatment clinical trials targeted on FGFRs have not made any progress (Repetto et al., 2021; Lengyel et al., 2022). Hence, the discovery of novel FGFRs inhibitors is necessary.
In our present study, we demonstrated a novel small molecule RK-019 by screening our library of pharmacologically active compounds. In the following research, we found RK-019 exhibited great inhibitory capacity and selectivity against the FGFRs family kinases, the IC50 were 9.1 nM (FGFR1), 4.6 nM (FGFR2), 26.3 nM (FGFR3), and 40.7 nM (FGFR4), respectively. This indicated RK-019 is a pan-FGFR inhibitor with great kinase inhibitory activity.
Then, we tested the anti-proliferative effect of RK-019 on multiple GC cell lines and found that it could efficiently inhibit the proliferation of FGFR2-amp GC cell lines, SNU-16 and KATO III, with average IC50 values of 3.96 ± 4.4 nM and 5.45 ± 5.3 nM, respectively. These IC50 values were significantly different from that of the gastric epithelial cells GES-1 or other types of GC cell lines, indicating that RK-019 could inhibit FGFR2-amp GC cells. However, prolonged treatment of SNU-16 cells with RK-019 (5 and 10 nM) showed contrasting results, likely due to the activation of drug resistance pathways. Thus, in the future, we will focus on identifying the mechanism of drug resistance and developing anti-drug resistance strategies to enhance the anti-tumor effect.
Our results verified that RK-019 could inhibit FGFR2T653/T654 auto-phosphorylation, which is the critical modification for the FGFR2 kinase activity and could prevent the phosphorylation of downstream proteins, including FRS2 and PLCγ. Phosphorylated FRS2 and PLCγ play a central role in the FGFR pathway activity, triggering PI3K-AKT-mTOR and Ras-Raf-Erk pathways (Katoh and Katoh, 2006; Acevedo et al., 2009; Ahmad et al., 2012; Ornitz and Itoh, 2015). In our study, we observed a decrease in the phosphorylation of AKT and Erk, which are known regulators of cell proliferation and survival. Furthermore, FCM and western blot results revealed G0/G1 phase arrest, validating our speculation about the inhibitory effects on the FGFR signaling.
Interestingly, although cell cycle arrest occurred immediately after RK-019 treatment, apoptosis was delayed. Apoptosis occurred after 72 h treatment, indicating that apoptosis might not be directly triggered by the inhibition of FGFR signaling. In addition, we found that apoptosis induced by RK-019 was not mediated by mitochondria. Mitochondrial membrane potential was increased, while the expression levels of mitochondrial apoptosis-related proteins did not change as expected after RK-019 treatment (Data not shown). Thus, we presume that RK-019 induced apoptosis via non-classic pathways.
Furthermore, we found that RK-019 could inhibit cell migration and invasion in SNU-16 and KATO III cells. Peritoneal metastasis is a leading cause of death in patients with gastric carcinoma. However, there are limited treatment options and no targeted therapy or immunotherapy for gastric carcinoma (Yao et al., 2020). RK-019 could inhibit the FGFR2 activity and repress JAK2 phosphorylation, thus mediating the STAT3 activity and decreasing the MMP-2 and MMP-9 levels. These findings demonstrate the ability of RK-019 to control cancer cell metastasis, which is crucial in GC therapy.
Finally, we verified the anti-tumor activity in vivo using the SNU-16 xenograft model. The results showed that RK-019 (30 mg/kg, daily oral administration for 21 days) could remarkably suppress SNU-16 tumor growth, with an inhibitory ratio of 85.9%. Furthermore, IHC revealed that RK-019 could prevent FGFR2 phosphorylation, induce apoptosis, and inhibit cell proliferation in the tumor sections.
Taken together, we have reported here a novel pan-FGFR inhibitor, RK-019, which exhibited excellent anti-tumorigenic activity against FGFR2-amp GC in vitro and in vivo.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding authors.
Ethics statement
The animal study was reviewed and approved by the Institutional Animal Care and Treatment Committee of Sichuan University.
Author contributions
LY is the project administrator. JZ, KR, YZ, ZL, and LY contributed to conception and design of the study. JZ, KR, XL, LT, and QW performed the experiments in this article. JZ, KR, and XL organized and validated the experiment results. JR gave material and resource support. JZ wrote the first draft of the manuscript. YZ, ZL, and LY performed reviewing and editing. All authors contributed to manuscript revision and read and approved the submitted version.
Funding
This study was supported by the National S&T Major Special Project on Major New Drug Innovations of China (2018ZX09201018), the Science and Technology Program of Sichuan Province (2021YFH0145), and the Scientific Research Project of the Health Planning Committee of Sichuan (20PJ108).
Conflict of interest
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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2022.998199/full#supplementary-material
References
Acevedo, V. D., Ittmann, M., and Spencer, D. M. (2009). Paths of FGFR-driven tumorigenesis. Cell Cycle 8 (4), 580–588. doi:10.4161/cc.8.4.7657
Ahmad, I., Iwata, T., and Leung, H. Y. (2012). Mechanisms of FGFR-mediated carcinogenesis. Biochim. Biophys. Acta 1823 (4), 850–860. doi:10.1016/j.bbamcr.2012.01.004
Ahn, S., Lee, J., Hong, M., Kim, S. T., Park, S. H., Choi, M. G., et al. (2016). FGFR2 in gastric cancer: Protein overexpression predicts gene amplification and high H-index predicts poor survival. Mod. Pathol. 29 (9), 1095–1103. doi:10.1038/modpathol.2016.96
Babina, I. S., and Turner, N. C. (2017). Advances and challenges in targeting FGFR signalling in cancer. Nat. Rev. Cancer 17 (5), 318–332. doi:10.1038/nrc.2017.8
Beenken, A., and Mohammadi, M. (2009). The FGF family: Biology, pathophysiology and therapy. Nat. Rev. Drug Discov. 8 (3), 235–253. doi:10.1038/nrd2792
Belov, A. A., and Mohammadi, M. (2013). Molecular mechanisms of fibroblast growth factor signaling in physiology and pathology. Cold Spring Harb. Perspect. Biol. 5 (6), a015958. doi:10.1101/cshperspect.a015958
Botrus, G., Raman, P., Oliver, T., and Bekaii-Saab, T. (2021). Infigratinib (BGJ398): An investigational agent for the treatment of FGFR-altered intrahepatic cholangiocarcinoma. Expert Opin. Investig. Drugs 30 (4), 309–316. doi:10.1080/13543784.2021.1864320
Brooks, A. N., Kilgour, E., and Smith, P. D. (2012). Molecular pathways: Fibroblast growth factor signaling: A new therapeutic opportunity in cancer. Clin. Cancer Res. 18 (7), 1855–1862. doi:10.1158/1078-0432.CCR-11-0699
Cailliau, K., Browaeys-Poly, E., and Vilain, J. P. (2001). Fibroblast growth factors 1 and 2 differently activate MAP kinase in Xenopus oocytes expressing fibroblast growth factor receptors 1 and 4. Biochim. Biophys. Acta 1538 (2-3), 228–233. doi:10.1016/s0167-4889(01)00074-x
Cerliani, J. P., Guillardoy, T., Giulianelli, S., Vaque, J. P., Gutkind, J. S., Vanzulli, S. I., et al. (2011). Interaction between FGFR-2, STAT5, and progesterone receptors in breast cancer. Cancer Res. 71 (10), 3720–3731. doi:10.1158/0008-5472.CAN-10-3074
Deo, D. D., Axelrad, T. W., Robert, E. G., Marcheselli, V., Bazan, N. G., and Hunt, J. D. (2002). Phosphorylation of STAT-3 in response to basic fibroblast growth factor occurs through a mechanism involving platelet-activating factor, JAK-2, and src in human umbilical vein endothelial cells: Evidence for a dual kinase mechanism. J. Biol. Chem. 277 (24), 21237–21245. doi:10.1074/jbc.M110955200
Eswarakumar, V. P., Lax, I., and Schlessinger, J. (2005). Cellular signaling by fibroblast growth factor receptors. Cytokine Growth Factor Rev. 16 (2), 139–149. doi:10.1016/j.cytogfr.2005.01.001
Fukumoto, S. (2008). Actions and mode of actions of FGF19 subfamily members. Endocr. J. 55 (1), 23–31. doi:10.1507/endocrj.kr07e-002
Goetz, R., and Mohammadi, M. (2013). Exploring mechanisms of FGF signalling through the lens of structural biology. Nat. Rev. Mol. Cell Biol. 14 (3), 166–180. doi:10.1038/nrm3528
Grose, R., and Dickson, C. (2005). Fibroblast growth factor signaling in tumorigenesis. Cytokine Growth Factor Rev. 16 (2), 179–186. doi:10.1016/j.cytogfr.2005.01.003
Gschwind, A., Fischer, O. M., and Ullrich, A. (2004). The discovery of receptor tyrosine kinases: Targets for cancer therapy. Nat. Rev. Cancer 4 (5), 361–370. doi:10.1038/nrc1360
Gullo, I., Carneiro, F., Oliveira, C., and Almeida, G. M. (2018). Heterogeneity in gastric cancer: From pure morphology to molecular classifications. Pathobiology 85 (1-2), 50–63. doi:10.1159/000473881
Han, N., Kim, M. A., Lee, H. S., and Kim, W. H. (2015). Evaluation of fibroblast growth factor receptor 2 expression, heterogeneity and clinical significance in gastric cancer. Pathobiology 82 (6), 269–279. doi:10.1159/000441149
Hosoda, K., Yamashita, K., Ushiku, H., Ema, A., Moriya, H., Mieno, H., et al. (2017). Prognostic relevance of FGFR2 expression in stage II/III gastric cancer with curative resection and S-1 chemotherapy. Oncol. Lett. 15, 1853–1860. doi:10.3892/ol.2017.7515
Jain, A., Borad, M. J., Kelley, R. K., Wang, Y., Abdel-Wahab, R., Meric-Bernstam, F., et al. (2018). Cholangiocarcinoma with FGFR genetic aberrations: A unique clinical phenotype. JCO Precis. Oncol. 2, 1–12. doi:10.1200/PO.17.00080
Jang, J.-H., Shin, K.-H., and Park, J.-G. (2001). Mutations in fibroblast growth factor receptor 2 and fibroblast growth factor receptor 3 genes associated with human gastric and colorectal cancers. Cancer Res. 61 (9), 3541. https://aacrjournals.org/cancerres/article/61/9/3541/508688/Mutations-in-Fibroblast-Growth-Factor-Receptor-2
Jang, J., Kim, H. K., Bang, H., Kim, S. T., Kim, S. Y., Park, S. H., et al. (2017). Antitumor effect of AZD4547 in a fibroblast growth factor receptor 2-amplified gastric cancer patient-derived cell model. Transl. Oncol. 10 (4), 469–475. doi:10.1016/j.tranon.2017.03.001
Jia, Y.-X., Li, T.-F., Zhang, D.-D., Fan, Z.-M., Fan, H.-J., Yan, J., et al. (2016). The coexpression and prognostic significance of c-MET, fibroblast growth factor receptor 2, and human epidermal growth factor receptor 2 in resected gastric cancer: A retrospective study. Onco. Targets. Ther. 9, 5919–5929. doi:10.2147/OTT.S111778
Jia, Z.-H., Jia, Y., Guo, F.-J., Chen, J., Zhang, X.-W., and Cui, M.-H. (2017). Phosphorylation of STAT3 at Tyr705 regulates MMP-9 production in epithelial ovarian cancer. Plos One 12 (8), e0183622. doi:10.1371/journal.pone.0183622
Kamata, T., Hattori, Y., Hamada, H., Kizaki, M., Terada, M., and Ikeda, Y. (2002). Keratinocyte growth factor regulates proliferation and differentiation of hematopoietic cells expressing the receptor gene K-sam. Exp. Hematol. 30 (4), 297–305. doi:10.1016/s0301-472x(01)00800-1
Katoh, M. (2010). Genetic alterations of FGF receptors: An emerging field in clinical cancer diagnostics and therapeutics. Expert Rev. Anticancer Ther. 10 (9), 1375–1379. doi:10.1586/era.10.128
Katoh, M., and Katoh, M. (2006). FGF signaling network in the gastrointestinal tract (review). Int. J. Oncol. 29 (1), 163–168. doi:10.3892/ijo.29.1.163
Kim, H. S., Kim, J. H., Jang, H. J., Han, B., and Zang, D. Y. (2019b). Pathological and prognostic impacts of FGFR2 overexpression in gastric cancer: A meta-analysis. J. Cancer 10 (1), 20–27. doi:10.7150/jca.28204
Kim, H. S., Kim, J. H., and Jang, H. J. (2019a). Pathologic and prognostic impacts of FGFR2 amplification in gastric cancer: A meta-analysis and systemic review. J. Cancer 10 (11), 2560–2567. doi:10.7150/jca.29184
Lamont, F., Tomlinson, D., Cooper, P. A., Shnyder, S. D., Chester, J. D., and Knowles, M. (2011). Small molecule FGF receptor inhibitors block FGFR-dependent urothelial carcinoma growth in vitro and in vivo. Br. J. Cancer 104 (1), 75–82. doi:10.1038/sj.bjc.6606016
Lawlor, M. A., and Alessi, D. R. (2001). PKB/Akt: A key mediator of cell proliferation, survival and insulin responses? J. Cell Sci. 114 (16), 2903–2910. doi:10.1242/jcs.114.16.2903
Lengyel, C. G., Hussain, S., Seeber, A., Jamil Nidhamalddin, S., Trapani, D., Habeeb, B. S., et al. (2022). FGFR pathway inhibition in gastric cancer: The golden era of an old target? Life 12 (1), 81. doi:10.3390/life12010081
Li, P., Huang, T., Zou, Q., Liu, D., Wang, Y., Tan, X., et al. (2019). FGFR2 promotes expression of PD-L1 in colorectal cancer via the JAK/STAT3 signaling pathway. J. Immunol. 202 (10), 3065–3075. doi:10.4049/jimmunol.1801199
Markham, A. (2019). Erdafitinib: First global approval. Drugs 79 (9), 1017–1021. doi:10.1007/s40265-019-01142-9
Massagué, J. (2004). G1 cell-cycle control and cancer. Nature 432 (7015), 298–306. doi:10.1038/nature03094
Matsumoto, K., Arao, T., Hamaguchi, T., Shimada, Y., Kato, K., Oda, I., et al. (2012). FGFR2 gene amplification and clinicopathological features in gastric cancer. Br. J. Cancer 106 (4), 727–732. doi:10.1038/bjc.2011.603
Montazeri, K., and Bellmunt, J. (2020). Erdafitinib for the treatment of metastatic bladder cancer. Expert Rev. Clin. Pharmacol. 13 (1), 1–6. doi:10.1080/17512433.2020.1702025
Network, C. G. A. R. (2014). Comprehensive molecular characterization of gastric adenocarcinoma. Nature 513 (7517), 202–209. doi:10.1038/nature13480
Oki, M., Yamamoto, H., Taniguchi, H., Adachi, Y., Imai, K., and Shinomura, Y. (2008). Overexpression of the receptor tyrosine kinase EphA4 in human gastric cancers. World J. Gastroenterol. 14 (37), 5650–5656. doi:10.3748/wjg.14.5650
Ornitz, D. M., and Itoh, N. (2001). Fibroblast growth factors. Genome Biol. 2 (3),S3005. doi:10.1186/gb-2001-2-3-reviews3005
Ornitz, D. M., and Itoh, N. (2015). The fibroblast growth factor signaling pathway. Wiley Interdiscip. Rev. Dev. Biol. 4 (3), 215–266. doi:10.1002/wdev.176
Piro, G., Carbone, C., Cataldo, I., Di Nicolantonio, F., Giacopuzzi, S., Aprile, G., et al. (2016). An FGFR3 autocrine loop sustains acquired resistance to trastuzumab in gastric cancer patients. Clin. Cancer Res. 22 (24), 6164–6175. doi:10.1158/1078-0432.ccr-16-0178
Repetto, M., Crimini, E., Giugliano, F., Morganti, S., Belli, C., and Curigliano, G. (2021). Selective FGFR/FGF pathway inhibitors: Inhibition strategies, clinical activities, resistance mutations, and future directions. Expert Rev. Clin. Pharmacol. 14 (10), 1233–1252. doi:10.1080/17512433.2021.1947246
Robertson, S. C., Tynan, J. A., and Donoghue, D. J. (2000). RTK mutations and human syndromes: When good receptors turn bad. Trends Genet. 16 (6), 265–271. doi:10.1016/s0168-9525(00)02021-7
Sarabipour, S., and Hristova, K. (2016). Mechanism of FGF receptor dimerization and activation. Nat. Commun. 7 (1), 10262–10312. doi:10.1038/ncomms10262
Shimizu, D., Saito, T., Ito, S., Masuda, T., Kurashige, J., Kuroda, Y., et al. (2018). Overexpression of FGFR1 promotes peritoneal dissemination via epithelial-to-mesenchymal transition in gastric cancer. Cancer Genomics Proteomics 15 (4), 313–320. doi:10.21873/cgp.20089
Simons, M. (2021). Fibroblast growth factors: The keepers of endothelial normalcy. J. Clin. Invest. 131 (17), e152716. doi:10.1172/jci152716
Smyth, E. C., Nilsson, M., Grabsch, H. I., van Grieken, N. C., and Lordick, F. (2020). Gastric cancer. Lancet 396 (10251), 635–648. doi:10.1016/S0140-6736(20)31288-5
Su, X., Zhan, P., Gavine, P., Morgan, S., Womack, C., Ni, X., et al. (2014). FGFR2 amplification has prognostic significance in gastric cancer: Results from a large international multicentre study. Br. J. Cancer 110 (4), 967–975. doi:10.1038/bjc.2013.802
Sung, H., Ferlay, J., Siegel, R. L., Laversanne, M., Soerjomataram, I., Jemal, A., et al. (2021). Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. Ca. Cancer J. Clin. 71 (3), 209–249. doi:10.3322/caac.21660
Takahashi, T., Saikawa, Y., and Kitagawa, Y. (2013). Gastric cancer: Current status of diagnosis and treatment. Cancers 5 (4), 48–63. doi:10.3390/cancers5010048
Tokunaga, R., Imamura, Y., Nakamura, K., Ishimoto, T., Nakagawa, S., Miyake, K., et al. (2016). Fibroblast growth factor receptor 2 expression, but not its genetic amplification, is associated with tumor growth and worse survival in esophagogastric junction adenocarcinoma. Oncotarget 7 (15), 19748–19761. doi:10.18632/oncotarget.7782
Tsang, M., and Dawid, I. B. (2004). Promotion and attenuation of FGF signaling through the Ras-MAPK pathway. Sci. STKE 2004 (228), pe17. doi:10.1126/stke.2282004pe17
Van Cutsem, E., Sagaert, X., Topal, B., Haustermans, K., and Prenen, H. (2016). Gastric cancer. Lancet 388 (10060), 2654–2664. doi:10.1016/S0140-6736(16)30354-3
Xie, T.-X., Wei, D., Liu, M., Gao, A. C., Ali-Osman, F., Sawaya, R., et al. (2004). Stat3 activation regulates the expression of matrix metalloproteinase-2 and tumor invasion and metastasis. Oncogene 23 (20), 3550–3560. doi:10.1038/sj.onc.1207383
Xu, H., Lee, K. W., and Goldfarb, M. (1998). Novel recognition motif on fibroblast growth factor receptor mediates direct association and activation of SNT adapter proteins. J. Biol. Chem. 273 (29), 17987–17990. doi:10.1074/jbc.273.29.17987
Yao, X., Ajani, J. A., and Song, S. (2020). Molecular biology and immunology of gastric cancer peritoneal metastasis. Transl. Gastroenterol. Hepatol. 5, 57. doi:10.21037/tgh.2020.02.08
Ye, Y.-W., Hu, S., Shi, Y.-Q., Zhang, X.-F., Zhou, Y., Zhao, C.-L., et al. (2013). Combination of the FGFR4 inhibitor PD173074 and 5-fluorouracil reduces proliferation and promotes apoptosis in gastric cancer. Oncol. Rep. 30 (6), 2777–2784. doi:10.3892/or.2013.2796
Ye, Y. W., Zhang, X., Zhou, Y., Wu, J., Zhao, C., Yuan, L., et al. (2012). The correlations between the expression of FGFR4 protein and clinicopathological parameters as well as prognosis of gastric cancer patients. J. Surg. Oncol. 106 (7), 872–879. doi:10.1002/jso.23153
Yue, S., Li, Y., Chen, X., Wang, J., Li, M., Chen, Y., et al. (2021). FGFR-TKI resistance in cancer: Current status and perspectives. J. Hematol. Oncol. 14 (1), 23. doi:10.1186/s13045-021-01040-2
Zhang, J., Tang, P. M. K., Zhou, Y., Cheng, A. S. L., Yu, J., Kang, W., et al. (2019). Targeting the oncogenic FGF-FGFR Axis in gastric carcinogenesis. Cells 8 (6), E637. doi:10.3390/cells8060637
Zhang, W., and Liu, H. T. (2002). MAPK signal pathways in the regulation of cell proliferation in mammalian cells. Cell Res. 12 (1), 9–18. doi:10.1038/sj.cr.7290105
Zou, L., Cao, S., Kang, N., Huebert, R. C., and Shah, V. H. (2012). Fibronectin induces endothelial cell migration through β1 integrin and Src-dependent phosphorylation of fibroblast growth factor receptor-1 at tyrosines 653/654 and 766. J. Biol. Chem. 287 (10), 7190–7202. doi:10.1074/jbc.M111.304972
Glossary
GC Gastric cancer
FGFR Fibroblast growth factor receptor
FGF Fibroblast growth factor
EMT Epithelial-to-mesenchymal transition
MTT 3-(4,5-Dimethyl-2-thiazolyl)-2,5-diphenyl-2-H-tetra-zolium bromide
DMSO Dimethyl sulfoxide
qRT-PCR Quantitative reverse transcription PCR
FCM Flow cytometry Flow cytometry
PBS Phosphate buffered solution
SDS-PAGE Sodium dodecyl sulphate-polyacrylamide gel electrophoresis
PVDF Polyvinylidene difluoride
SPF Specific-pathogen free
LC-MS/MS Liquid chromatography–mass spectrometry/mass spectrometry
PK Pharmacokinetics
FCM Flow cytometry Flow cytometry
PI Propidium iodide
PARP Poly (ADP-ribose) polymerase
CBC Complete blood count
TGI Tumor growth inhibition
IC50 Median-inhibitory concentration
ATP Adenosine triphosphate
AUC Area under the plasma concentration time curve
IHC Immunohistochemistry
FRS2 Fibroblast growth factor receptor substrate 2
PLCγ Phosphoinositide phospholipase C gamma
Erk Extracellular signal-regulated kinase
JAK Janus kinase
STAT Signal transducer and activator of transcription
CDK Cyclin-dependent kinases
MMP Matrix metalloproteinases
PI3K Phosphatidylinositol 3 kinase.
Keywords: gastric cancer, receptor tyrosine kinase, antineoplastic agents, FGFRs, cell apoptosis
Citation: Zeng J, Ran K, Li X, Tao L, Wang Q, Ren J, Hu R, Zhu Y, Liu Z and Yu L (2022) A novel small molecule RK-019 inhibits FGFR2-amplification gastric cancer cell proliferation and induces apoptosis in vitro and in vivo. Front. Pharmacol. 13:998199. doi: 10.3389/fphar.2022.998199
Received: 19 July 2022; Accepted: 30 August 2022;
Published: 21 September 2022.
Edited by:
Chen Shi, Huazhong University of Science and Technology, ChinaReviewed by:
Xiaoyan Huang, Shaanxi Provincial People’s Hospital, ChinaDmitry A. Khochenkov, Russian Cancer Research Center NN Blokhin, Russia
Copyright © 2022 Zeng, Ran, Li, Tao, Wang, Ren, Hu, Zhu, Liu and Yu. 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: Luoting Yu, eXVsdW90QHNjdS5lZHUuY24=; Zhihao Liu, bGl1emhpaGFvQHNjdS5lZHUuY24=; Yongxia Zhu, emh1eW9uZ3hpYTA0MDlAMTYzLmNvbQ==
†These authors contributed equally to this work and share first authorship