Follow us 登录

KRAS小分子抑制剂研究进展

约30%人类癌症与RAS突变有关。RAS家族有3种亚型,即KRAS,HRAS和NRAS。KRAS已成为人类癌症中最常出现突变的基因之一,在多种癌症如胰腺癌、肺癌和结肠癌等均发现KRAS蛋白被上调或被突变激活。多年来,针对KRAS突变抗癌药物的临床活性并不乐观。安进公司利用KRAS突变产生的半胱氨酸残基成功开发上市了共价抑制剂AMG510,使KRAS再次成为新药研发热点。本文以小分子与KRAS蛋白的相互作用模式为依据,将KRAS小分子抑制剂分为直接作用抑制剂和间接作用抑制剂,综述了各类抑制剂的代表结构、作用机制和生物活性研究进展,为KRAS小分子抑制剂抗肿瘤新药的开发提供思路和方向。



正文

癌基因RAS是首个被发现的人类癌基因,主要有 KRAS,HRAS 和 NRAS 3 种亚型。RAS突变引发的癌症约占人类癌症的30%,而KRAS基因突变占 RAS 基因突变总数的 85%(NRAS 占12%,HRAS占3%),KRAS也被誉为癌症治疗的“靶点之王”。


KRAS蛋白具有低疏水性,经过法尼基化修饰后可克服自身低疏水性,结合到细胞膜内侧。通常,KRAS在细胞内以无活性的KRAS-鸟苷二磷酸(guanosine diphosphate,GDP)形式存在,在受体蛋白酪氨酸激酶受到刺激信号刺激后,可招募下游鸟嘌呤核苷转换因子(guanine nucleotide exchangefactors,GEF)。在 GEF作用下,KRAS释放 GDP,并迅速与鸟苷三磷酸(guanosine triphosphate,GTP)结合,以活性形式 KRAS-GTP存在。KRAS蛋白自身具有 GTP 酶活性,在 GTP 酶激活蛋白(GTPase-activating proteins,GAP)调控下活性增强,可将结合的GTP水解为GDP,使KRAS-GTP转变为无活性的KRAS-GDP。KRAS正是通过这种非活性状态(与GDP结合)和活性状态(与GTP结合)之间的切换,发挥二元开关的作用。这个过程在正常细胞中虽非常缓慢,但大多数KRAS突变会抑制GTP水解,影响KRAS-GTP向KRAS-GDP转换,进而引起KRAS功能异常激活,引发癌症。


针对这一调控机制,设计开发直接作用于KRAS的抑制剂,进而抑制其功能的异常激活,即成为新药研发的重要思路。也可针对 KRAS细胞分布及其活性增强的调控,设计间接作用于KRAS的抑制剂。如抑制KRAS释放GDP的一种鸟苷释放蛋白激活剂 SOS1(son of sevenless 1,SOS1)抑制剂、影响 KRAS 分布的法尼基转移酶抑制剂和KRAS作用于下游效应物的含Src同源2结构域蛋白酪氨酸磷酸酶(Src homology 2 domain contain⁃ing protein tyrosine phosphatase,SHP2)抑制剂等,都是重要的间接作用于KRAS的抑制剂。


KRAS参与调控3条重要信号通路:迅速加速纤维肉瘤激酶(rapidly accelerated fibrosarcoma,Raf)/丝裂原活化蛋白激酶激酶(mitogen-activatedprotein kinase kinase,MEK)/细胞外信号调节激酶(extracellular signal-regulated kinase,ERK)、磷脂酰肌醇3-激酶(phosphatidylinositol-4,5-bispho⁃sphate 3-kinase,PI3K)-蛋白激酶B(protein kinaseB, Akt)-哺乳动物雷帕霉素靶蛋白(mammalian tar⁃get of Rapamycin,mTOR)和 Ral/鸟嘌呤核苷酸解离刺激因子(Ral/GDP dissociation stimulator,Ral/GDS)。因而也可以设计作用于KRAS下游相关效应物的间接抑制剂。本文对直接作用抑制剂和间接作用抑制剂的代表结构、分子机制和生物活性研究进展进行综述。


1 直接作用于KRAS的抑制剂

KRAS是拥有 GTP酶活性的鸟嘌呤核苷酸结合蛋白,主要由1个GTP酶催化结构域和紧随其后的高度可变区域组成(图1) ,其直接靶向抑制剂主要包括:① 直接作用于GTP酶活性位点的鸟苷类似物;② 直接作用于开关Ⅱ口袋的抑制剂;③ 其他方式直接作用于KRAS的抑制剂;④ 基于蛋白降解 靶 向 嵌 合 体(proteolysis targeting chimeras,PROTAC)技术靶向诱导KRAS突变蛋白降解。性位点的竞争性可逆抑制剂的开发长期未能取得突破,KRAS一度被认为是“不可成药靶点”。


image.png

Hunter等和Chen等报道的SML系列化合物能与KRAS的GTP酶活性位点发生不可逆共价结合,很好地解决了竞争性可逆抑制剂对 GTP酶活性位点结合能力弱的问题。SML系列化合物是一种GDP类似物,含有一个从β磷酸盐延伸出来的亲电弹头,经过迈克尔反应加到Cys12上,形成一个稳定的硫醚键,研究发现,SML 阻止了 KRAS与下游效应物Raf的结合。但该类化合物作为竞争性不可逆抑制剂在RAS不同亚基之间的选择性仍然是化合物结构设计的挑战。该类抑制剂的代表化合物SML-8-73-1(1,图2)在1 mmol·L -1 GDP/GTP环境中,与 KRAS G12C 的结合率超过 95%,抑制肺癌细胞 H358 增殖的 IC 50 值为26.6 μmol·L -1 。

image.png

1.1 作用于GTP酶活性位点


KRAS突变主要发生在Gly12,Gly13和Gln61残基上,而这3个氨基酸处于KRAS蛋白发挥GTP酶活性的催化位点,突变后对KRAS蛋白的GTP酶活性影响显著,因此找到可与KRAS的GTP酶活性位点结合的小分子抑制剂,对抑制KRAS异常激活意义重大。早期的 KRAS抑制剂开发思路即是寻找能竞争性结合KRAS的GTP酶活性位点的鸟苷类似物,但由于RAS各亚型中GTP酶结构域的序列和结构几乎相同,各亚型之间很难做到高选择性;且KRAS与GDP/GTP的亲和力常数达pmol水平,而细胞中 GTP浓度在 mmol水平,小分子药物很难对KRAS进行竞争性结合;同时KRAS蛋白表面缺少有利于小分子结合的空腔,使基于GTP酶活Boumelha等于2021年在美国癌症研究学会上报道了化合物RM-032,该分子可与 KRA SG12C和亲环蛋白 A 共价结合,抑制 KRAS G12C肿 瘤 细胞中 RAS通路信号传导和细胞增殖。RM-032具有口服利用度高、体外对细胞增殖抑制明显的特点,但其详细结构未见报道。


1. 2 作用于开关Ⅱ口袋


1.2.1 作用于KRAS G12C突变型

2013年,Janes等首次报道了利用小分子共价结合KRAS G12C突变体的可行性。对于KRASG12C型的突变,突变产生的半胱氨酸(C12)残基与相邻的组氨酸95(H95)残基产生了一个易于靶向的狭窄口袋开关Ⅱ口袋,可提供稳定的药物-蛋白质相互作用位点。小分子通过不可逆共价结合于此“口袋”,可将KRAS G12C锁定在非活性状态,从而影响肿瘤细胞的增殖扩散,而不破坏未突变KRAS 的正常功能。由于开关Ⅱ口袋的发现和不可逆共价结合技术的发展,使得直接靶向KRAS抑制剂的研发重新成为热点,直接靶向KRAS的小分子不可逆共价抑制剂不断被报道。


安进公司研发的AMG510(2,图3)是第一个进入临床研究的KRAS共价抑制剂,在30 min内可缓慢地将KRAS-TGP转换为KRAS-GDP。AMG510分子与 KRAS G12C共结晶(图 4)显示,其母核喹唑啉酮占据 KRAS 开关Ⅱ口袋,丙烯酰胺部分与C12形成共价键,吡啶环的取代基(异丙基)与Y96,H95和Q99紧密接触,很大程度上填充了H95残基的侧链旋转产生的口袋。此时,KRAS蛋白GTP酶活性位点无法再容纳3个磷酸基,与GTP的结合力大大降低,同时也会阻碍GEF催化GTP替换GDP的过程,将KRAS G12C锁定在结合GDP的非活性状态,阻止肿瘤细胞增殖扩散.

image.png


image.png

image.png

阿达格拉西布(adagrasib,MRTX849)是另一种可口服的 KRAS G12C 小分子共价抑制剂(3,图3),具有与AMG510相同的作用机制。该化合物的上市申请已于2022年2月15日向美国食品药品管理局(FDA)提交,用于治疗KRAS G12C突变且此前接受过至少1次全身治疗的非小细胞肺癌患者。


国内外其他公司基于开关Ⅱ口袋作用机制设计开发的KRAS G12C共价小分子抑制剂化学结构见图 3,其体外细胞增殖抑制活性见表 1。其中贝达、信达、先声、劲方药业和首药控股的研制的该类化合物细胞增殖抑制活性IC 50 值达几十nmol水平,其他公司的小分子抑制剂紧随其后,均表现出良好的抑制活性。


1.2.2 作用于KRAS G12D突变型

KRAS G12D突变也是人类恶性肿瘤中最主要的变异之一,研究可抑制KRAS G12D的小分子抑制剂对治疗恶性肿瘤具有极高的价值。清华大学张永辉团队报道的针对KRAS G12D突变的抑制剂TH-Z816(18,图5),具有与KRAS G12C抑制剂截然不同的结构和作用机制。对 TH-Z816 与KRAS G12D 的共晶结构研究发现,KRAS G12D化合物结合口袋也由开关Ⅱ口袋(红色)、P-环(蓝绿色)和α2-螺旋(绿色)等蛋白结构组成(图6A)。TH-Z816的萘基部分深深嵌入开关Ⅱ口袋中,与残基Val9,Met72,Phe78,Gln99,Ile100和Val103形成疏水相互作用,同时还与残基 His95,Glu62,Gly60 和 Asp12 之间也有 4 对极性相互作用(图6B)。突变体中Asp12的羧酸基在生理条件下脱去质子后带负电荷,与碱性的、带正电荷的哌嗪基团形成正、负电荷相互吸引的盐桥。这种正、负电荷之间的强盐桥相互作用类似于建筑结构中的榫卯结构,使蛋白与小分子药物像榫和卯一样牢固连接(图 6C)。这种结合也“挤压”了 KRAS G12D 的GTP酶活性位点,使其结合GTP的能力大大降低,从而抑制KRAS因突变产生的过度活化,这一发现为KRAS G12D抑制剂的开发提供了很好的思路。相较针对KRAS G12C突变开发的不可逆共价抑制剂,盐桥在提供牢固的结合力的同时也使这种结合具有可逆性。因此,以这种作用机制开发的抑制剂属于变构位点的可逆抑制剂。

image.png

image.png

在张永辉团队之前,最早报道的KRAS G12D可 逆 抑 制 剂 是 Mirati Therapeutics 公 司 研 发 的MRTX1133(19,图 5),其具有与 TH-Z816 相同的哌嗪结构,并与其作用机制类似。该化合物对KRAS G12D 的 抑 制 具 有 选 择 性 和 可 逆 性 ,对KRAS 野生型细胞无影响。在胰腺癌和结直肠癌的动物荷瘤模型上均显示出对 KRAS依赖性信号的持续剂量依赖性抑制,且在 3 mg·kg -1 水平即可有效抑制 G12D 突变型移植瘤生长,但至今未见其进入临床研究的报道。我国恒瑞医药披露,其针对 KRAS G12D 开发的抑制剂 HRS-4642 注射液于 2022 年 7 月 26 日获批临床,无疑成为国内 KRAS G12D 抑制剂药物研发的先行者,但未见其化合物结构的报道。同时,国内加科思公司的化合物 20(图 5)和南方医科大学陈建军团队的化合物 KD8(21,图 5)也是含哌嗪结构的抑制剂。



1.3 其他方式直接作用于KRAS


Kessler等开发的BI-2852(22,图7A),结合在 KRAS上的开关Ⅰ口袋和开关Ⅱ口袋间的特殊口袋中。其与 KRAS的活性或非活性形式结合的亲和力达nmol水平。晶体结构研究发现,BI-2852与KRAS的活性形式KRAS-GTP结合时,开关Ⅰ口袋(图7B浅绿色)口袋和开关Ⅱ口袋(图7B浅红色)口袋相较于与非活性形式KRAS-GDP结合时均显著减小。通过这种特定的结合方式,不但阻断KRAS-GDP 与 SOS1 之间的相互作用,也抑制KRAS-GTP与其下游效应物Raf和PI3K之间的相互作用,这导致 GDP 到 GTP 的交换被抑制,GAP催化的GTP到GDP的交换也被抑制,从而可很好地抑制NCI-H358肿瘤细胞增殖 。BI-2852最终虽未能成药,但证明了开关Ⅰ口袋和开关Ⅱ口袋确实是可能的结合位点,给抑制KRAS药物的开发提供更多可能。


image.png


1.4 KRAS-PROTAC抑制剂


PROTAC技术是利用小分子化合物诱导靶向蛋白降解的技术。PROTAC分子是一种双功能小分子,一端结合靶蛋白的配体,另一端结合E3泛素连接酶的配体,结构中 2 个配体之间通过连接子(linker)相连。其在体内可将靶蛋白和E3酶拉近,使靶蛋白被打上泛素标签,然后通过泛素 — —蛋白酶体途径降解靶蛋白。Bond等报道了LC-2(23,图 8),它是第一个能够降解内源性 KRASG12C的 PROTAC 分子,LC-2 用阿达格拉西布弹头与 KRAS G12C 共价结合,并招募 E3 连接酶VHL,从而诱导快速而持续的 KRAS G12C 降解,导致纯合和杂合 KRAS G12C 细胞系中 MAPK 信号通路的抑制;在 NCI-H2030 细胞中,LC-2 很好地 诱 导 了 内 源 性 KRAS G12C 降 解 ,IC 50 值 为(0.59±0.20)μmol·L -1 。PROTAC 技术虽起步较晚,但因其在药物开发上的独特优势而备受关注,加之KRAS抑制剂研究日益广泛,该类降解剂的报道也不断涌现。Li等报道的化合物KP-14(24,图8)与LC-2具有类似的结构,其对NCI-H358细胞增殖的IC 50 值为17.41 μmol·L -1 。


2 间接作用于KRAS的抑制剂

KRAS在胞内信号传递通路中处于关键位置,通过阻断 KRAS上、下游信号转导,可开发出多种间接抑制KRAS的小分子化合物。主要包括:影响KRAS 在细胞中正常分布小分子抑制剂,抑制KRAS 释放 GDP 的 SOS1 抑制剂,抑制 KRAS 脱磷酸化的 SHP2 抑制剂和抑制KRAS调控的下游信号通路的抑制剂等。



2.1 影响KRAS在细胞中正常分布的抑制剂


2.1.1 法尼基转移酶抑制剂

法尼基化是存在于真核细胞中的异戊二烯化修饰形式,它可以修饰蛋白并使其锚定到胞膜上。实体瘤增殖与转移中常见法尼基转移酶的参与,且与转移酶的活性有关。通过抑制法尼基转移酶的活性,可以阻断 KRAS蛋白法尼基修饰,切断其与 细 胞 膜 结 合 路 径 ,从 而 抑 制 肿 瘤 细 胞 的 增殖。替吡法尼(tipifarnib,R115777)(25,图 9)是选择性抑制法尼基转移酶的小分子抑制剂,基于乳腺癌和肾癌的适应证的开发处于临床Ⅱ期试验阶段,基于胰腺癌适应证的开发处于Ⅲ期临床试验阶段。


2.1.2 PDE-δ δ靶向降解剂

KRAS为了到达质膜,其法尼基化的高变区与PDE-δ的穿梭因子结合。Cheng等报道了第一类 PDE-δ 降解物化合物 26(图 9)。该化合物与PDE-δ结合,并有效地诱导PDE-δ降解,最终诱导KRAS突变体SW480和HCT116细胞的凋亡,并在体内、外对KRAS突变的结直肠癌细胞均表现出明显的抗增殖能力,对 SW480 细胞增殖的 IC 50 值为(8.0±2.7)μmol·L -1 ,对HCT116细胞增殖的IC 50 值为(8.0±3.6)μmol·L -1 。



2.2 SOS1抑制剂


SOS1蛋白是一种在细胞中广泛表达的调控蛋白,作为信号通路中的关键蛋白,其在细胞内许多信号转导通路中起重要调控作用。SOS1是GEF和RAS的激活剂,可在变构位点结合KRAS-GTP,从而增强GEF功能,构成正反馈调节,其缺失或其GEF 功能特异性基因失活已被证明会降低携带KRAS突变的肿瘤细胞存活率。研发有效且选择性SOS1 抑制剂将与 MEK 抑制剂协同作用,从而在KRAS驱动癌症中产生强大且持续的通路阻断和良好的抗肿瘤功效,成为KRAS小分子抑制剂的研究方向。


image.png

BI-3406(27,图9)是一种高效、高选择性、高口服生物利用度的小分子SOS1抑制剂。该抑制剂可和 SOS1 催化结构域结合,从而阻断 SOS1 与KRAS之间相互作用。BI-3406可抑制肿瘤生长,在MIA-PaCa-2 荷瘤小鼠中,每日 2 次给予 BI-340650 mg·kg -1 与曲美替尼(trametinib)0.125 mg·kg -1 ,可在长达10 h内持续抑制肿瘤细胞的增殖;该化合物代谢约24 h后,MIA-PaCa-2肿瘤的ERK磷酸化水平恢复到基础水平。由于 SOS1 在所有肿瘤类型中均表达,其抑制剂 BI-3406 可广泛应用于KRAS突变引起的肿瘤治疗。



2.3 SHP2抑制剂


SHP2磷酸酶由PTPN11基因编码,参与多种细胞致癌信号级联反应,是第一个报道的致癌磷酸酶。SHP2可直接使KRAS脱磷酸化,增强其与效应蛋白RAF的结合,从而激活下游MEK/ERK信号通路。TNO155(28,图9)是诺华制药公司通过优化吡嗪类变构开发的SHP2抑制剂,其对SHP2的IC 50 为0.011 mmol·L -1[60] ,TNO155可抑制RAS/RAF/ERK 信号通路,并且在小鼠和大鼠中的口服生物利用度分别为 78% 和 86% 。JAB-3068(29,图9)是由加科思制药公司开发的SHP2的小分子抑制剂,可干扰 MAPK 信号通路,抑制表达SHP2 的肿瘤细胞增殖。Nichols 等报道了SHP2 的小分子抑制剂 RMC-4550(30,图 9),其IC 50 为0.583 nmol·L -1 ,可通过抑制KRAS蛋白的活性阻断RAS/ERK信号通路和肿瘤生长。SHP099(31,图9)也是一种有效的SHP2抑制剂。SHP099的极性官能团和离子官能团暴露,显著提高了其选择性和口服生物利用度,其对 ERK磷酸化的IC 50 为0.07 mmol·L -1。



2.4 抑制RAS/RAF/MEK/ERK信号通路的抑制剂


在抗癌药物研究信号通路筛选中,癌症细胞生长以及侵袭和转移等过程均有 RAS/RAF/MEK/ERK通路参与。由于该信号通路负责细胞增殖,因此通过研究小分子抑制剂将期阻断成为重要的研究方向。在RAS/RAF/MEK/ERK信号通路中,MAPK表达异常是判断癌症发生的标志,目前阻断KRAS突变肿瘤的策略主要集中在抑制下游MAPK信号转导。索拉非尼(sorafenib)(32,图 10)是可口服的、可进行多靶点治疗的KRAS小分子抑制剂。通过作用于RAF/MEK/ERK/信号通路,其可直接抑制信号转导,从而阻止肿瘤细胞增殖。可比替尼(cobimetinib,GDC - 0973)(33,图 10)是MEK/Akt双通路抑制剂,该化合物的杂氮环丁烷及酰胺键增加了细胞代谢稳定性,可通过抑制KRAS G13D和B-RAF G464V突变从而抑制肿瘤发生。



2.5 抑制PI3K/Akt/mTOR信号通路的抑制剂


布帕利西(buparlisib,BKM120)(34,图10)是一种广谱 PI3K 抑制剂,作用于 p110α,p110β,p110δ 和 p110γ 时,IC 50 值分别为 52,166,116 和262 nmol·L -1 。BEZ-235(35,图10)是咪唑并喹啉衍生物,是 PI3K 和 mTOR 的小分子抑制剂 。利用BEZ-235对p0α,p0γ,p0δ,p0β和mTOR作用时,IC 50 值分别为4,5,7,75 和20.7 nmol·L -1 。另外研究显示,BEZ-235在临床上可抑制晚期乳腺癌细胞增长,对实体肿瘤也具有治疗作用。BEZ-235抑制PI3K和mTOR突变,阻断PI3K/Akt/mTOR信号通路,抑制肿瘤细胞抑制肿瘤细胞的恶性增殖。临床证明,该药可维持糖代谢稳定,增强肿瘤细胞对细胞毒性药物如长春新碱等的敏感性。


image.png


2.6 抑制Ral/GDS/Ral信号通路的抑制剂


鸟嘌呤核苷酸交换因子,具有直接和特异性催化GDP/GTP交换能力,从而促进其活性信号状态,减少Ral A/B与鸟苷核苷酸亲和能力,有利于结合GTP。Ral/GDS/Ral信号通路负责细胞因子释放,能促进Ral蛋白GDP/GTP转换。Ral由Ral A和Ral B亚型组成,Ral A对肿瘤生长起重要作用,而Ral B在肿瘤转移和侵袭中起重要作用。通过抑制 Ral/GDS/Ral 通路可抑制 Ral 下游效应信号通路。研究发现,二氢吡咯环衍生物P61A6(结构未报道)可有效抑制非小细胞肺癌异种移植模型中肿瘤细胞增殖,P61A6对H358,H23和H1507细胞的IC 50 值为5~15 μmol·L -1。


3 结语与展望

KRAS蛋白具有独特的结构特点,由于表面相对光滑且口袋腔较浅,与药物分子结合存在一定难度,一度被认为是“不可成药”的靶点。近年来共价抑制剂技术的发展以及对KRAS结构认识的深入,特别是 AMG510 发现以来,直接或间接靶向抑制KRAS突变的小分子抑制剂相继出现,使KRAS成为抗癌药物热门研究靶点。国内外已有多家企业聚焦基于KRAS靶点的创新药开发,申报临床的有十余种,前景值得期待。其中再鼎医药从 MIrati公司引进的阿达格拉西布进度较快,信达生物的 GF-105、加科思JAB-21822、益方生物D-1553、贝达药业BPI-421286和恒瑞医药HRS-4642等均已进入临床研究。先声、济民可信的化合物均处在非临床药学研究,细胞增殖抑制活性均达到nmol水平,其他公司的小分子抑制剂也紧随其后,均表现出良好抑制活性。目前,KRAS小分子抑制剂的研究重心主要是针对 KRAS G12C。随着研究的深入,针对KRAS G12D与KRAS/SOS1相互作用的抑制剂及KRAS降解剂也将不断涌现。同时,目前研究结果显示,KRAS靶向药物与其他免疫药物或靶向药物联用具有更突出效果,如阿达格拉西布与PD-1抗体潘利珠单抗(pembrolizumab)联用治疗携带KRAS G12C突变的非小细胞肺癌患者,Ⅰ期临床试验结果显示疾病控制率达到100%。此外,PROTAC技术凭借其独特优势将在针对KRAS靶点的药物开发中发挥出巨大潜力。


因KRAS突变导致的肿瘤种类较多,从适应证等方面来看,KRAS是十分值得深入广泛研究的癌症治疗靶点,但目前国内针对该靶点的小分子药物研究同质化严重。就已公布的小分子化合物看,结构差异不大,适应证集中,作用效果相当,这可能导致最终的新药申报、临床选用和市场占有出现十分激烈的竞争。目前 KRAS G12C 不可逆抑制剂AMG510已上市,阿达格拉西布已递交上市申请,二者的IC 50 值在40~100 nmol·L -1 ,临床用药剂量较大,其中后者临床有效剂量为600 mg,AMG510>960 mg,且AMG510已出现耐药性。因此,发现活性更高、生物利用度更好的新化合物仍具有重要意义。提高化合物的耐药性也将是 KRAS抑制剂未来发展的挑战和机遇。




作者:焦龙华1,2 ,田从魁 1,2 ,王胜超 1,3 ,王永广 1,4 ,郜红利 1,2 ,史伯安 1,2


机构:

1.湖北民族大学武陵山天然药物研究所,湖北 恩施 445000

2. 风湿性疾病发生与干预湖北省重点实验室,湖北 恩施 445000

3. 宜昌东阳光长江药业股份有限公司,湖北 宜昌 443000

4. 苏州施安鼎泰生物医药技术有限公司,江苏 苏州 215000

来源:中国药理学与毒理学杂志2023年3月第37卷第3期




Contact Us

Feel free to call us on
025-85998075

Our Email

Drop us a line anytime at
sales@popebiotech.com,
and we’ll get back soon.

Our Address

Come visit us at 南京市江宁区科学园乾德路5号