Literature DB >> 21378983

Characterization of a selective inhibitor of the Parkinson's disease kinase LRRK2.

Xianming Deng1, Nicolas Dzamko, Alan Prescott, Paul Davies, Qingsong Liu, Qingkai Yang, Jiing-Dwan Lee, Matthew P Patricelli, Tyzoon K Nomanbhoy, Dario R Alessi, Nathanael S Gray.   

Abstract

Mutations in leucine-rich repeat kinase 2 (LRRK2) are strongly associated with late-onset autosomal dominant Parkinson's disease. We employed a new, parallel, compound-centric approach to identify a potent and selective LRRK2 inhibitor, LRRK2-IN-1, and demonstrated that inhibition of LRRK2 induces dephosphorylation of Ser910 and Ser935 and accumulation of LRRK2 within aggregate structures. LRRK2-IN-1 will serve as a versatile tool to pharmacologically interrogate LRRK2 biology and study its role in Parkinson's disease.

Entities:  

Mesh:

Substances:

Year:  2011        PMID: 21378983      PMCID: PMC3287420          DOI: 10.1038/nchembio.538

Source DB:  PubMed          Journal:  Nat Chem Biol        ISSN: 1552-4450            Impact factor:   15.040


Parkinson’s disease (PD) is the second most common neurodegenerative disease in the world, with a prevalence of more than 1% after the age of 65 years[1]. Mutations in the Leucine-rich repeat kinase-2 (LRRK2) gene were first linked to PD in 2004[2,3] and extensive genetic, cellular and biochemical studies suggest that it is a key driver of the disease. Among the genes associated with PD, such as parkin and PINK1, LRRK2 is unique because a missense mutation, G2019S, is frequently found in not only familial but also sporadic PD cases[4,5]. The most common and extensively characterized mutation, G2019S, enhances kinase activity, suggesting that small molecule LRRK2 kinase inhibitors may serve as potential therapeutic agents for the treatment of PD[6]. To date, only relatively nonspecific LRRK2 inhibitors have been reported such as staurosporine[7,8], K252A[7,8], sunitinib[7-9], H-1152[9], and GW5074[10]. Here we report the identification of a selective inhibitor of LRRK2 and demonstrate its ability both in vitro and in vivo to inhibit LRRK2. In an effort to discover small molecules that can selectively inhibit LRRK2 kinase activity, we applied a novel, highly parallel, compound-centric strategy[11]. We screened a 300-membered compound library designed to target the ATP-binding site against a panel of 442 diverse kinases using an in vitro ATP-site competition binding assay[12,13]. In contrast to the traditional, linear process of inhibitor discovery, this high-throughput kinase profiling enables rapid identification of compounds capable of selectively inhibiting LRRK2 and other kinases of interest. A series of 2-amino-5,11-dimethyl-5H-benzo[e]pyrimido-[5,4-b][1,4]diazepin-6(11H)-one analogs emerged as potential LRRK2 inhibitors with strong binding affinities against both wild-type LRRK2 and LRRK2 mutant G2019S. LRRK2-IN-1 was identified as a promising inhibitor following the synthesis of approximately 50 analogs guided by iterative rounds of synthesis and testing in biochemical and cellular assays (to be described elsewhere) (Fig. 1a and Supplementary Methods). LRRK2-IN-1 inhibits both wild-type and G2019S mutant LRRK2 kinase activity with IC50 values of 13 nM and 6 nM (Fig. 1b), respectively with 0.1 mM ATP in the assay. Consistent with LRRK2-IN-1 being an ATP competitive inhibitor, the IC50 for G2019S mutant was increased when ATP was increased (Supplementary Fig. 1). Previous work identified a LRRK2[A2016T] mutant that is normally active, but 32-fold less sensitive to H-1152 and 12-fold less sensitive to sunitinib[9]. Remarkably, LRRK2[A2016T] and LRRK2[A2016T+G2019S] mutant were ~400-fold more resistant to LRRK2-IN-1 (Fig. 1b). A molecular model of LRRK2-IN-1 binding to the ATP-site of LRRK2 demonstrates that the A2016T mutation likely introduces a disfavorable steric interaction with the anthranilic acid phenyl group of LRRK2-IN-1 (Supplementary Fig. 2).
Figure 1

Enzymatic activity of LRRK2-IN-1 and its selectivity

(a) Chemical structure of the LRRK2 inhibitor LRRK2-IN-1. (b) GST-LRRK2 (1326-2517), GST-LRRK2 [G2019S] (1326-2527), GST-LRRK2 [A2016T] (1326-2517) and GST-LRRK2 [A2016T+G2019S] (1326-2517) were assayed using 20µM Nictide in the presence of 100 µM ATP with the indicated concentrations of LRRK2-IN-1. The results are presented as percentage of kinase activity relative to the DMSO treated control. Results are averages of duplicate reactions with similar results obtained in at least one other experiment. The IC50 values, in µM, were derived from the graphs. (c) The combined kinase profiling hits with biochemical activities. Kinase hits from Ambit profiling with a score of less than 10% of the DMSO control at a concentration of 10 µM, ActivX profiling with higher than 50% inhibition at 1 µM, and Dundee profiling with greater than 50% inhibition at 1 µM are listed. The biochemical IC50 values of hits were measured by Invitrogen Adapta®, Z´-Lyte®, Lantha-Screen® assays (http://www.invitrogen.com/site/us/en/home/Products-and-Services/Services/Screening-and-Profiling-Services/SelectScreen-Profiling-ervice/SelectScreen-Kinase-Profiling-Service.html) and Dundee radioactive-base enzyme assay[14]. *Assessment of PLK1 cellular activity of LRRK2-IN-1 (see Supplementary Results).

The kinase selectivity of LRRK2-IN-1 was comprehensively evaluated using three independent and complementary experimental approaches that in aggregate assessed the selectivity against >470 kinases. These approaches included kinase-binding assays against a panel of 442 distinct kinases using the KINOMEscan™ methodology (Ambit Biosciences, San Diego, CA)[12,13], standard radioactive-based enzymatic assays against a panel of 105 kinases (Dundee profiling)[14], and activity-based proteomics against 260 kinases using the KiNativ™ technology (ActivX Biosciences, San Diego, CA)[15]. This analysis revealed that LRRK2-IN-1 possessed a selective profile (see Supplementary Dataset), inhibiting only 12 kinases among the ambit 442 diverse kinases panel with a score of less than 10% of the DMSO control at a concentration of 10 µM. In the biochemical assays, LRRK2-IN-1 possessed an IC50 of 45 nM for inhibition of DCLK2 and exhibited an IC50 of greater than 1 µM when evaluated in biochemical assays for AURKB, CHEK2, MKNK2, MYLK, NUAK1, and PLK1 (Supplementary Table 1). Biochemical assays are not available for MAPK7, DCLK1, PLK4, RPS6KA2 (Kin.Dom.2-C-terminal), and RPS6KA6 (Kin.Dom.2-C-terminal), therefore the dissociation constants (Kds) were measured by KINOMEscan binding assay (Supplementary Table 1). The selectivity score (S (3 µM)) for LRRK2-IN-1 calculated by dividing the number of kinases found to bind with a dissociation constant (Kd < 3 µM) by the total number of kinases tested is 0.029 (13/442), which represents high selectivity (see Supplementary Dataset)[13]. Further exploration of the selectivity of LRRK2-IN-1 was performed by profiling directly against all kinases detectable in human peripheral blood mononuclear cells (hPBMC, 180 kinases) and mouse brain tissue (195 kinases) using KiNativ™. This approach confirmed inhibition of LRRK2, DCLK1 and MAPK7, and supported >1 µM IC50s for CHEK2, and MYLK (smMLCK) and RPS6KA2. LRRK2-IN-1 was confirmed to inhibit MAPK7 with an EC50 of 160 nM (Supplementary Fig. 3). We are unable to evaluate whether LRRK2-IN-1 inhibits LRRK1 as LRRK1 kinase activity cannot be measured under the same conditions employed for LRRK2 (Supplementary Fig. 4). Overall the selectivity of LRRK2-IN-1 is clearly sufficient for use as a tool compound for assessing the effects of LRRK2 inhibition in cellular systems. LRRK2 binds to isoforms of 14-3-3 proteins and this interaction depends on the phosphorylation of LRRK2 at Ser910 and Ser935[16]. Recent work suggested that inhibition of LRRK2 kinase activity stimulates an as yet uncharacterized regulatory feedback loop promoting the release of 14-3-3 and dephosphorylation of LRRK2 at Ser910 and Ser935[17]. Indeed, LRRK2-IN-1 induced a dose dependent inhibition of Ser910 and Ser935 phosphorylation accompanied by loss of 14-3-3 binding to both wild type LRRK2 (Fig. 2a) and LRRK2[G2019S] stably transfected into HEK293 cells(Fig. 2b). Significant dephosphorylation of Ser910/Ser935 was observed at 1-3 µM of LRRK2-IN-1 for wild type LRRK2 and at slightly lower doses for LRRK2[G2019S] (Fig. 2a & 2b). Crucially, LRRK2-IN-1 failed to induce dephosphorylation of Ser910/Ser935 in the drug resistant LRRK2[A2016T] (Fig. 2c) and LRRK2[A2016T+G2019S] (Fig. 2d) mutants, establishing that the effects of LRRK2-IN-1 on LRRK2 phosphorylation are indeed mediated via inhibition of LRRK2 kinase activity and not due to off target effects.
Figure 2

LRRK2-IN-1 inhibits LRRK2 in vivo

(a) HEK 293 cells stably expressing wild type GFP-LRRK2 were treated with DMSO or increasing concentrations of LRRK2-IN-1 for 90 min. Cell lysates were subjected to immunoblotting for detection of LRRK2 phosphorylated at S910 and S935 and for total GFP-LRRK2. 14-3-3 binding was measured by overlay assay following immunoprecipitation of LRRK2 with GFP-binder Sepharose. (b) As in (a) except HEK 293 cells stably expressing GFP-LRRK2 [G2019S] were utilized. (c) As in (a) except HEK 293 cells stably expressing drug resistant GFP-LRRK2 [A2016T] were utilized. (d) As in (a) except HEK 293 cells stably expressing drug resistant GFPLRRK2[ A2016T+G2019S] cells were utilized. (e) Endogenous LRRK2 was immunoprecipitated with anti-LRRK2 100-500 from EBV immortalized human lymphoblastoid cells from a control subject and a Parkinson’s disease patient homozygous for the LRRK2 [G2019S] mutation. (f) Following treatment of the cells with DMSO or the indicated concentrations of LRRK2-IN-1 for 90 min. Immunoprecipitates were subjected to immunoblot analysis with the indicated antibody as well as 14-3-3 overlay for western analysis. Immunoprecipitations were performed in duplicate and results are representative of at least two independent experiments. (g) As in (e) except that human SHSY5Y cells were utilized. (h) 100 mg/kg LRRK2-IN-1 was administered by intraperitoneal injection into wild type mice (n=3 per group). Tissues were collected and endogenous LRRK2 was immunoprecipitated and analysed for S910 and S935 phosphorylation as above. Quantitation of immunoblots along with original images are available in Supplementary Fig. 11

Near maximal effects of 1 µM LRRK2-IN-1 on dephosphorylation of Ser910/Ser935 and 14-3-3 binding to wild type LRRK2 and LRRK2[G2019S] were observed within 15 min and sustained for at least 2 h (Supplementary Fig. 5). Phosphorylation of Ser910/Ser935 and 14-3-3 binding was restored to near control levels within 2 hours of removal of LRRK2-IN-1 from the tissue culture medium (Supplementary Fig. 5). We observed that kinase inactive LRRK2 was phosphorylated at Ser910/Ser935 to a slightly lower level than wild type LRRK2 and treatment with LRRK2-IN-1 did not induce significant dephosphorylation of kinase-inactive LRRK2 (Supplementary Fig. 6). We next examined the effect of LRRK2-IN-1 on endogenous LRRK2 phosphorylation and 14-3-3 binding in human lymphoblastoid cells derived from a control individual (Fig. 2e) and a Parkinson’s disease patient homozygous for the LRRK2 G2019S mutation (Fig. 2f). We found that increasing doses of LRRK2-IN-1 led to similar dephosphorylation of endogenous LRRK2 at Ser910 and Ser935, accompanied by loss of 14-3-3 binding, as was observed in HEK 293 cells stably expressing wild-type LRRK2 (compare Fig. 2a and 2e) or LRRK2[G2019S] (compare Fig. 2b and 2f). Moreover, endogenous LRRK2[G2019S] was also more sensitive to LRRK2-IN-1 than wild type LRRK2 (compare Fig. 2e and 2f). We also found that LRRK2-IN-1 induced a similar dose dependent Ser910/Ser935 dephosphorylation and loss of 14-3-3 binding to endogenous LRRK2 in human-derived neuroblastoma SHSY5Y cells (Fig. 2g) as well as mouse Swiss 3T3 cells (Supplementary Fig. 7). LRRK2-IN-1 exhibits favourable pharmacokinetics in mice with a T1/2 of 4.5 hr, AUC of 14758 hr*ng/mL and %F of 49.3% (Supplementary Table 2). Following intraperitoneal injection of LRRK2-IN-1 into mice (100 mg/kg) we observed complete Ser910/Ser935 dephosphorylation of LRRK2 in the kidney at the 1 h and 2 h time points (Fig. 2h). In contrast, no dephosphorylation of LRRK2 Ser910 or Ser935 was observed in brain, indicating that LRRK2-IN-1 is likely unable to efficiently cross the blood brain barrier (Fig. 2h). Previous work suggested that dephosphorylation of Ser910/Ser935 and loss of 14-3-3-binding, promoted relocalization of LRRK2 to more aggregate and fibrillar-like structures[17] that might resemble inclusion bodies that certain mutant forms of LRRK2 have been reported to accumulate within[16,18,19]. To verify how LRRK2-IN-1 affected LRRK2 localization, we treated HEK293 cell lines stably expressing forms of full length GFP-LRRK2 with 1 µM LRRK2-IN-1 over a 2 h time course. Prior to addition of LRRK2-IN-1, mutant LRRK2[G2019S] (Fig. 3a) and wild type LRRK2 (Supplementary Fig. 8a) were diffusely localized in the cytosol consistent with previous work[17]. Following addition of LRRK2-IN-1, we observed a significant accumulation of LRRK2[G2019S] (Fig. 3a, Movie 1 and Supplementary Fig. 9) and wild type LRRK2 (Supplementary Fig. 8a and Movie 2) into aggregate structures. Importantly, LRRK2-IN-1 did not affect the localization of drug resistant LRRK2[A2016T] (Supplementary Fig. 8b and Movie 3) or LRRK2[A2016T+G2019S] (Supplementary Fig. 8c and Movie 4) mutants, that remained diffusely localized in the cytoplasm. Consistent with LRRK2 becoming re-phosphorylated on Ser910 and Ser935 following inhibitor removal (Supplementary Fig. 5), LRRK2[G2019S] mutant (Fig. 3a & Movie 2) and wild type LRRK2 (Supplementary Fig. 8a and Movie 1) returned to a largely diffuse cytoplasmic re-localization within 2 hours of LRRK2-IN-1 removal. LRRK1 does not bind 14-3-3 and the localization of LRRK1 does not change with LRRK2-IN-1 treatment (Supplementary Fig. 10).
Figure 3

LRRK2-IN-1 alters the cytoplasmic localisation of LRRK2

HEK 293 cells stably expressing GFP-LRRK2 [G2019S] were treated with 1 µM LRRK2-IN-1 for 2 h. Cells were then washed and fixed with 4% paraformaldehyde for fluorescence imaging. To visualize LRRK2 localization following removal of LRRK2-IN-1, inhibitor-containing media was removed 2 h after treatment. Cells were then washed twice and incubated with fresh media without inhibitor for 2 h before again being washed and then fixed for fluorescence imaging. Representative micrographs of three independent experiments are presented.

In conclusion, we report upon the discovery of LRRK2-IN-1, the first selective and potent inhibitor of LRRK2. We establish that LRRK2-IN-1 rapidly suppresses LRRK2 kinase activity in vivo leading to dephosphorylation of Ser910/Ser935, loss of 14-3-3 binding and accumulation of LRRK2 within aggregate fibrillar structures. LRRK2-IN-1 has the potential to transform our understanding of LRRK2 function, in the same way that other signal transduction inhibitors including PD98059[20], SB203580[21] and LY294002[22] have aided with elucidation of the roles of Mek, p38 and PI 3-kinase respectively. We recommend use of drug resistant LRRK2[A2016T] mutant for validating that observed effects of LRRK2-IN-1 are indeed mediated via inhibition of LRRK2 kinase activity. LRRK2-IN-1 provides a useful first-generation ‘tool’ compound and provides a starting point for developing a compound that could ultimately be used to address the potential therapeutic benefit of inhibiting LRRK2 in Parkinson’s patients harboring the LRRK2 G2019S mutation.
  21 in total

Review 1.  LRRK2 and Parkinson disease.

Authors:  Justus C Dächsel; Matthew J Farrer
Journal:  Arch Neurol       Date:  2010-05

2.  A specific inhibitor of phosphatidylinositol 3-kinase, 2-(4-morpholinyl)-8-phenyl-4H-1-benzopyran-4-one (LY294002).

Authors:  C J Vlahos; W F Matter; K Y Hui; R F Brown
Journal:  J Biol Chem       Date:  1994-02-18       Impact factor: 5.157

3.  Phenotype, genotype, and worldwide genetic penetrance of LRRK2-associated Parkinson's disease: a case-control study.

Authors:  Daniel G Healy; Mario Falchi; Sean S O'Sullivan; Vincenzo Bonifati; Alexandra Durr; Susan Bressman; Alexis Brice; Jan Aasly; Cyrus P Zabetian; Stefano Goldwurm; Joaquim J Ferreira; Eduardo Tolosa; Denise M Kay; Christine Klein; David R Williams; Connie Marras; Anthony E Lang; Zbigniew K Wszolek; Jose Berciano; Anthony H V Schapira; Timothy Lynch; Kailash P Bhatia; Thomas Gasser; Andrew J Lees; Nicholas W Wood
Journal:  Lancet Neurol       Date:  2008-06-06       Impact factor: 44.182

4.  14-3-3 binding to LRRK2 is disrupted by multiple Parkinson's disease-associated mutations and regulates cytoplasmic localization.

Authors:  R Jeremy Nichols; Nicolas Dzamko; Nicholas A Morrice; David G Campbell; Maria Deak; Alban Ordureau; Thomas Macartney; Youren Tong; Jie Shen; Alan R Prescott; Dario R Alessi
Journal:  Biochem J       Date:  2010-09-15       Impact factor: 3.857

5.  LRRK2 regulates autophagic activity and localizes to specific membrane microdomains in a novel human genomic reporter cellular model.

Authors:  Javier Alegre-Abarrategui; Helen Christian; Michele M P Lufino; Ruxandra Mutihac; Lara Lourenço Venda; Olaf Ansorge; Richard Wade-Martins
Journal:  Hum Mol Genet       Date:  2009-07-29       Impact factor: 6.150

6.  Substrate specificity and inhibitors of LRRK2, a protein kinase mutated in Parkinson's disease.

Authors:  R Jeremy Nichols; Nicolas Dzamko; Jessica E Hutti; Lewis C Cantley; Maria Deak; Jennifer Moran; Paul Bamborough; Alastair D Reith; Dario R Alessi
Journal:  Biochem J       Date:  2009-10-23       Impact factor: 3.857

Review 7.  Leucine-rich repeat kinase 2 (LRRK2): a key player in the pathogenesis of Parkinson's disease.

Authors:  Payal N Gandhi; Shu G Chen; Amy L Wilson-Delfosse
Journal:  J Neurosci Res       Date:  2009-05-01       Impact factor: 4.164

8.  Identification of compounds that inhibit the kinase activity of leucine-rich repeat kinase 2.

Authors:  Jason P Covy; Benoit I Giasson
Journal:  Biochem Biophys Res Commun       Date:  2008-11-21       Impact factor: 3.575

9.  Inhibition of LRRK2 kinase activity leads to dephosphorylation of Ser(910)/Ser(935), disruption of 14-3-3 binding and altered cytoplasmic localization.

Authors:  Nicolas Dzamko; Maria Deak; Faycal Hentati; Alastair D Reith; Alan R Prescott; Dario R Alessi; R Jeremy Nichols
Journal:  Biochem J       Date:  2010-09-15       Impact factor: 3.857

10.  Inhibitors of leucine-rich repeat kinase-2 protect against models of Parkinson's disease.

Authors:  Byoung Dae Lee; Joo-Ho Shin; Jackalina VanKampen; Leonard Petrucelli; Andrew B West; Han Seok Ko; Yun-Il Lee; Kathleen A Maguire-Zeiss; William J Bowers; Howard J Federoff; Valina L Dawson; Ted M Dawson
Journal:  Nat Med       Date:  2010-08-22       Impact factor: 53.440

View more
  179 in total

1.  The Upshot of LRRK2 Inhibition to Parkinson's Disease Paradigm.

Authors:  A R Esteves; M G-Fernandes; D Santos; C Januário; S M Cardoso
Journal:  Mol Neurobiol       Date:  2014-11-15       Impact factor: 5.590

2.  Comprehensive analysis of kinase inhibitor selectivity.

Authors:  Mindy I Davis; Jeremy P Hunt; Sanna Herrgard; Pietro Ciceri; Lisa M Wodicka; Gabriel Pallares; Michael Hocker; Daniel K Treiber; Patrick P Zarrinkar
Journal:  Nat Biotechnol       Date:  2011-10-30       Impact factor: 54.908

3.  How chemoproteomics can enable drug discovery and development.

Authors:  Raymond E Moellering; Benjamin F Cravatt
Journal:  Chem Biol       Date:  2012-01-27

4.  Leucine-rich repeat kinase 2 modulates neuroinflammation and neurotoxicity in models of human immunodeficiency virus 1-associated neurocognitive disorders.

Authors:  Jenna M Puccini; Daniel F Marker; Tim Fitzgerald; Justin Barbieri; Christopher S Kim; Patrick Miller-Rhodes; Shao-Ming Lu; Stephen Dewhurst; Harris A Gelbard
Journal:  J Neurosci       Date:  2015-04-01       Impact factor: 6.167

5.  Mitochondrial clearance and maturation of autophagosomes are compromised in LRRK2 G2019S familial Parkinson's disease patient fibroblasts.

Authors:  Joanna A Korecka; Ria Thomas; Dan P Christensen; Anthony J Hinrich; Eliza J Ferrari; Simon A Levy; Michelle L Hastings; Penelope J Hallett; Ole Isacson
Journal:  Hum Mol Genet       Date:  2019-10-01       Impact factor: 6.150

6.  The G2385R variant of leucine-rich repeat kinase 2 associated with Parkinson's disease is a partial loss-of-function mutation.

Authors:  Iakov N Rudenko; Alice Kaganovich; David N Hauser; Aleksandra Beylina; Ruth Chia; Jinhui Ding; Dragan Maric; Howard Jaffe; Mark R Cookson
Journal:  Biochem J       Date:  2012-08-15       Impact factor: 3.857

7.  LRRK2 secretion in exosomes is regulated by 14-3-3.

Authors:  Kyle B Fraser; Mark S Moehle; João P L Daher; Philip J Webber; Jeri Y Williams; Carrie A Stewart; Talene A Yacoubian; Rita M Cowell; Terje Dokland; Tong Ye; Dongquan Chen; Gene P Siegal; Robert A Galemmo; Elpida Tsika; Darren J Moore; David G Standaert; Kyoko Kojima; James A Mobley; Andrew B West
Journal:  Hum Mol Genet       Date:  2013-07-25       Impact factor: 6.150

8.  Chemoproteomics-based design of potent LRRK2-selective lead compounds that attenuate Parkinson's disease-related toxicity in human neurons.

Authors:  Nigel Ramsden; Jessica Perrin; Zhao Ren; Byoung Dae Lee; Nico Zinn; Valina L Dawson; Danny Tam; Michael Bova; Manja Lang; Gerard Drewes; Marcus Bantscheff; Frederique Bard; Ted M Dawson; Carsten Hopf
Journal:  ACS Chem Biol       Date:  2011-08-10       Impact factor: 5.100

Review 9.  In Vitro Modeling of Leucine-Rich Repeat Kinase 2 G2019S-Mediated Parkinson's Disease Pathology.

Authors:  Scott C Vermilyea; Marina E Emborg
Journal:  Stem Cells Dev       Date:  2018-03-29       Impact factor: 3.272

10.  LRRK2 regulates synaptogenesis and dopamine receptor activation through modulation of PKA activity.

Authors:  Loukia Parisiadou; Jia Yu; Carmelo Sgobio; Chengsong Xie; Guoxiang Liu; Lixin Sun; Xing-Long Gu; Xian Lin; Nicole A Crowley; David M Lovinger; Huaibin Cai
Journal:  Nat Neurosci       Date:  2014-01-26       Impact factor: 24.884

View more

北京卡尤迪生物科技股份有限公司 © 2022-2023.