从中药等天然产物混合物备选库中筛选sars冠状病毒主蛋白酶抑制剂的新方法

文档序号:6127149阅读:645来源:国知局
专利名称:从中药等天然产物混合物备选库中筛选sars冠状病毒主蛋白酶抑制剂的新方法
技术领域
本发明涉及药物的筛选方法,更具体而言,涉及以SARS冠状病毒主蛋白酶的晶体结构为基础,从中药等天然产物混合物备选库中筛选SARS冠状病毒主蛋白酶的抑制剂的新方法及基于此方法发现的SARS冠状病毒主蛋白酶的抑制剂。
背景技术
SARS(Severe Acute Respiratory Sydrome)为引起严重急性呼吸道综合症(非典型性肺炎)的简称,其病原体为冠状病毒属的一种病毒(Peiris J.et al.Lancet,2003,3611319-1325)。冠状病毒是正链RNA病毒。冠状病毒属隶属于冠状病毒科。在目前已知的正链RNA病毒中,它们的基因组是最大的(Siddell,S.G.等人Coronaviruses,toroviruses,and arteriviruses.in Topley & Wilson′s Microbiology andMicrobia Infections,10th edition,Vol.Virology(eds.Mahy,B.W.J.& ter Meulen,V.)823-856(Hodder Arnold,London,2005)。这个属含有约26个种;按照它们的自然宿主、基因序列以及血清型关系,这个属又可以被分为三个组群(group)其中第一组群有TGEV,Porcine Transmissible Gastroenteritis Virus(猪传染性胃肠炎病毒)等;第二组群有SARS冠状病毒等;第三组群有AIBV,Avian InfectiousBronchitis Virus(禽传染性支气管炎病毒)等(Spaan,W.J.M.and Cavanagh,D.Coronaviridae.in Virus taxonomy,VIIIth Report of the ICTV.945-62(Elsevier-Academic Press.,London,2004)。
SARS冠状病毒2/3到3/4的基因组编码了两个复制酶多蛋白(replicasepolyproteins)pp1a和pp1ab(Ziebuhr,Snijder et al.2000;Thiel,Herold et al.2001;Anand,Yang et al.2005;Ziebuhr 2005),它们只有在病毒编码的蛋白酶切割成独立亚基之后才能使病毒完成正常转录、复制功能(Ziebuhr,Snijder et al.2000;Anand,Yang et al.2005;Bartlam,Yang et al.2005)。SARS冠状病毒主要蛋白酶(mainprotease,简称主蛋白酶)在这个过程中起主要作用。如果能够抑制SARS冠状病毒主蛋白酶的水解作用,那么将会有效地抵御SARS冠状病毒对人体的侵染。因此,SARS冠状病毒的主蛋白酶是抗SARS药物筛选的一个主要靶标。
天然产物又称次级代谢产物(Secondary Metabolites),具有结构的多样性和生物活性的多样性。天然产物及其衍生物在以往的疾病治疗中发挥了无可限量的作用,也是当今药物研发过程中最具潜力的资源之一(Newman DJ,GraggGM,Snader KM.Nat Prod Rep,2000,17(3)215-234;Lee K-H.J NatProd.2004,67(2)273-283)。现在对中药等传统药物、海洋生物以及微生物代谢过程中的天然产物研究方兴未艾,每年都会有大量结构新颖的化合物被发现提供出来,这些新颖化合物是合成方法所无法实现的药物和药物先导化合物的重要源泉,在新药和先导化合物的发现中起着重要作用。对于天然产物尤其是来源于中药等的天然产物的研究是很有必要的,因为中药在我国已有数千年应用的历史,是一个伟大的药物小分子化合物的宝库,因此从中药等天然产物中进行重要病毒或重要疾病相关靶蛋白抑制剂的筛选是很有必要的。
陈曙霞的CN1600308A(2005.03.30)中公开了从众多的中药中筛选出唯一一味低毒有效成分——槐果碱单体(sophocarpine,简称SC),其用于抗SARS病毒是有效果。但是其实质上并不是对中药整体的筛选,而是根据基于SARS病毒与柯萨奇B病毒存在的共同点,而槐果碱单体对于柯萨奇B病毒有治疗效果而推导得出的,其中并没有阐述该单体的作用机理。
徐筱杰等的CN1602853A(2005.04.06)中公开了一种抑制SARS冠状病毒感染的中药有效成分及其生物活性测定方法,具体涉及从植物五倍子中分离出的抑制SARS冠状病毒侵染宿主细胞的有效成分,该成分针对的是冠状病毒的关键结构域S2蛋白。黄建平的CN1803168A(2006.07.19)公开了一种治疗严重急性呼吸综合征(SARS)的中药,该中药来源于临床治疗ARDS有效验方,依据上述SARS病因病机及治则加减组方,由大青叶(干品)15-30g,生桅子5-15g,金银花15-30g,炒僵蚕10-15g,玄参10-15g,牡丹皮10-15g,紫丹参15-30g,西洋参5-10g,生甘草15-30g,桑白皮10-30g,葶苈子15-30g,汉防己5-10g,生大黄5-10g(后下)共十三味中药材组成,具有清热平喘,解毒化浊,凉血消斑,活血化瘀,利水祛湿,泻肺通腑,益气养阴之功效,主要用于SARS喘憋期治疗,入汤剂口服,一日两剂,一次煎取药液600毫升,每隔1-2小时喝一次。该文献中并没有对各种中药的作用机理做出任何阐述。
目前还没有以SARS冠状病毒主蛋白酶的晶体结构为基础,从中药等天然药物系统(如中药等天然产物混合物备选库)中来筛选SARS冠状病毒主蛋白酶抑制剂的技术。本发明提供了一种从来源于中药等天然产物混合物备选样品中筛选靶向与SARS冠状病毒主蛋白酶结合,且结合能力最好的小分子抑制剂的方法,并从原子水平阐释了其对SARS冠状病毒主蛋白酶的抑制机理。

发明内容
本发明提供了一种从天然产物混合物备选库中筛选SARS冠状病毒主蛋白酶抑制剂的方法,此方法包括步骤A测定天然产物混合物备选库中的备选样品对SARS冠状病毒主蛋白酶的抑制活性,获得对SARS冠状病毒主蛋白酶有抑制活性的样品,这些样品称为“有抑制活性的样品”,将SARS冠状病毒主蛋白酶的晶体与“有抑制活性的样品”接触,得到SARS冠状病毒主蛋白酶与小分子抑制剂复合体的晶体。SARS冠状病毒主蛋白酶的晶体与“有抑制活性的样品”接触时,根据需要,可分别与单一的“有抑制活性的样品”接触,也可与混合的“有抑制活性的样品”接触。
在得到SARS冠状病毒主蛋白酶与小分子抑制剂复合体的晶体后,再进行步骤B收集并分析步骤A中得到的复合体的晶体的X射线衍射数据,获得与SARS冠状病毒主蛋白酶结合的小分子抑制剂的电子密度。
获得与SARS冠状病毒主蛋白酶结合的小分子抑制剂的电子密度后,再进行步骤C以步骤B中获得的与SARS冠状病毒的主蛋白酶结合的小分子抑制剂的电子密度为基础,鉴定与SARS冠状病毒主蛋白酶结合的小分子抑制剂,获得小分子抑制剂与SARS冠状病毒主蛋白酶复合体的精细三维结构。
在步骤C之后,再进行步骤D测定小分子抑制剂对SARS冠状病毒主蛋白酶的抑制活性。
本发明中用来筛选SARS冠状病毒主蛋白酶抑制剂的天然产物混合物备选库,可以是以任何天然产物为原料,以任何方法制备得到的备选样品。
本发明优选用如下初级原料制备用于筛选SARS冠状病毒主蛋白酶抑制剂的天然产物混合物备选库包括传统药材、按中药处方规定的药材种类和用量组成的混合药材等中药,还未进行药用开发的动植物品种、海洋生物、微生物等等。
由于中药等天然产物粗提物所含的成分是相当复杂的,因此必须对天然产物的粗提物进行优化,目的是去除对SARS冠状病毒主蛋白酶体外活性测定和晶体浸泡实验中的干扰因素如蛋白、多糖等,保留成药性较好的小分子并使之尽量集中,提高相对浓度,以制备适用于以SARS冠状病毒主蛋白酶晶体结构为基础进行抑制剂筛选的备选库。本发明提供了可以实现这一目的制备备选库的技术方案,其基本原理为对初级原料提取后,提取物先使用不同极性的溶剂如乙酸乙酯进行萃取、萃取后的水相再用大孔树脂分离的方法进行处理,这样既能保证去除蛋白、多糖等干扰因素,又能利用两种不同分离机制的差别,使小分子集中到极性不同的段位,由此得到含有适中分子量的并能够用于以SARS冠状病毒主蛋白酶晶体结构为基础进行抑制剂筛选的天然产物混合物备选库。
本发明提供的制备用于筛选SARS冠状病毒主蛋白酶抑制剂天然产物混合物备选库的方法,具体包括下述步骤(a)、取一定量的初级原料,用乙醇溶液回流提取,合并提取液,浓缩得浸膏;(b)、称取一定量步骤(a)中所得的浸膏,用双蒸水悬浮,用有机溶剂萃取,分离有机溶剂相和水相;合并水相,备用;合并每次有机溶剂萃取物,蒸除其中的有机溶剂,得到干膏状有机溶剂相样品,备用;(c)、蒸除步骤(b)中萃取后的水相部分的水中溶解的少量有机溶剂,用大孔树脂柱层析分离,依次用水和乙醇溶液为流动相洗脱;弃去水洗部分,蒸除乙醇溶液洗脱部分的溶剂,得干膏状样品备用;
(d)、上述步骤(b)中制得的干膏状有机溶剂相样品和步骤(c)中制得的干膏状样品即为源自一种初级原料的2个备选样品。
本发明提供的制备天然产物混合物备选库的优选方法为步骤(a)中乙醇溶液与初级原料的用量质量比为6∶1-12∶1;所用乙醇溶液浓度为70%-95%;步骤(b)中浸膏与悬浮用双蒸水的质量比为1∶5-1∶15,萃取用有机溶剂为乙酸乙酯或氯仿,其体积与悬浮用双蒸水的体积比为6∶1-10∶1;步骤(c)中所用乙醇溶液浓度为50%。
上述制备方法中,优选的初级原料为中药,更优选的初级原料为鱼腥草,野菊花,大青叶,芦根,金银花,连翘,黄芩,知母,四季青,穿心莲,鸭跖草,栀子,射干,山豆根,马勃,玄参,地骨皮、川藏香茶菜、鄂西香茶菜。
此外,本发明提供的从天然产物混合物备选库中筛选SARS冠状病毒主蛋白酶抑制剂的方法,更优选在步骤A中SARS冠状病毒主蛋白酶晶体与“有抑制活性的样品”接触时,使用的方法为晶体浸泡或共结晶的晶体浸泡方法。
在一个实施方案中,本发明获得了SARS冠状病毒主蛋白酶-川藏香茶菜丙素复合体的晶体。
在另一个实施方案中,本发明筛选得到了SARS冠状病毒主蛋白酶的小分子抑制剂川藏香茶菜丙素。
在另一个实施方案,本发明提供了川藏香茶菜可以用于治疗或预防SARS冠状病毒、传染性胃肠炎病毒或者禽传染性支气管炎病毒等冠状病毒感染。
在另一个实施方案中,本发明提供的川藏香茶菜丙素能够有效抑制冠状病毒主蛋白酶活性,可以用于治疗或者预防SARS冠状病毒、传染性胃肠炎病毒或者禽传染性支气管炎病毒等冠状病毒感染。
在另一个实施方案中,本发明提供了能够有效抑制SARS冠状病毒主蛋白酶活性的川藏香茶菜乙酸乙酯萃取样品、鄂西香茶菜乙酸乙酯萃取样品、黄芩乙酸乙酯萃取样品、鱼腥草乙酸乙酯萃取样品、山豆根乙酸乙酯萃取样品、栀子乙酸乙酯萃取样品、知母乙酸乙酯萃取样品7种中药样品及其制备方法。
本发明选取SARS冠状病毒中晶体结构已知的主蛋白酶,对我国传统中药等天然产物混合物样品进行活性成分筛选,是从一种新的角度去研究中药学、天然产物学和生物学。其实质内容可包括两方面第一,靶蛋白的表达、纯化及结晶,这属于生物大分子的研究范畴;第二,使中药等天然产物混合物样品中非孤立存在的天然产物小分子自由地与靶蛋白——SARS冠状病毒主蛋白酶进行竞争性的结合,这属于化学特别是天然产物化学研究的范畴。
本发明的抑制剂筛选的方法以SARS冠状病毒主蛋白酶晶体结构为基础,从中药等天然产物混合物备选库中进行筛选,能够快速、准确的寻找对SARS冠状病毒主蛋白酶的特异性的抑制剂,提高了筛选的准确性。本发明的方法相对简单、成本较低,但筛选效率高、快捷、准确,并可以从混合物样品中,得到结合能力最好的抑制剂小分子。
具体而言,与现有方法相比,本发明的抑制剂的筛选方法的优点如下(1)首先,该方法省去了繁琐的天然产物分离纯化步骤,大大减少工作量,而且本方法的中药等天然产物混合物备选库的中的样品适宜于进行SARS冠状病毒主蛋白酶的体外活性测定和X射线晶体学测试,极大地提高了筛选效率和命中率。
(2)使用SARS冠状病毒主蛋白酶的晶体作为诱饵,浸泡有活性的天然产物混合物样品,可以把对SARS冠状病毒主蛋白酶有抑制作用的小分子抑制剂从混合物样品中分离出来,获得SARS冠状病毒主蛋白酶与小分子抑制剂结合的复合体的精细三维结构,并从原子水平上阐释该小分子抑制剂对SARS冠状病毒主蛋白酶的抑制机理。。
(3)使用本发明方法,通过使用SARS冠状病毒主蛋白酶浸泡天然产物混合物样品,可以从天然产物混合物样品中找到与SARS冠状病毒主蛋白酶结合能力最好的小分子(竞争结合的原理),优化选择了小分子抑制剂。


图1本发明的野菊花,芦根,金银花,连翘,四季青,鸭跖草,栀子,射干,马勃,玄参,地骨皮的水相50%乙醇洗脱样品对SARS冠状病毒主蛋白酶的抑制活性曲线图。样品浓度100μg/mL图2本发明的鱼腥草,大青叶,黄芩,知母,穿心莲,山豆根的水相50%乙醇洗脱样品对SARS冠状病毒主蛋白酶的抑制活性曲线图。样品浓度100μg/mL图3本发明的金银花,黄芩,穿心莲,四季青,鸭跖草,马勃,玄参,芦根,大青叶,地骨皮的乙酸乙酯萃取样品对SARS冠状病毒主蛋白酶的抑制活性曲线图。样品浓度100μg/mL图4本发明的鱼腥草,野菊花,连翘,知母,栀子,射干,山豆根的乙酸乙酯萃取样品对SARS冠状病毒主蛋白酶的抑制活性曲线图。样品浓度100μg/mL图5本发明的鱼腥草、山豆根、栀子、知母、黄芩的乙酸乙酯萃取样品对SARS冠状病毒主蛋白酶的抑制活性曲线图。样品浓度10μg/mL图6本发明的鄂西香茶菜乙酸乙酯萃取样品、水相50%乙醇洗脱样品对SARS冠状病毒主蛋白酶的抑制活性曲线图。样品浓度100μg/mL图7本发明的鄂西香茶菜乙酸乙酯萃取样品、水相50%乙醇洗脱样品对SARS冠状病毒主蛋白酶的抑制活性曲线图。样品浓度10μg/mL图8本发明的川藏香茶菜乙酸乙酯萃取样品、水相50%乙醇洗脱样品对SARS冠状病毒主蛋白酶的抑制活性曲线图。样品浓度100μg/mL图9本发明的川藏香茶菜乙酸乙酯萃取样品、水相50%乙醇洗脱样品对SARS冠状病毒主蛋白酶的抑制活性曲线图。样品浓度10μg/mL图10本发明的SARS冠状病毒主蛋白酶活性位点附近电子密度对比图。a为SARS冠状病毒主蛋白酶母体晶体活性位点的电子密度显示;b为SARS冠状病毒主蛋白酶与川藏香茶菜乙酸乙酯萃取样品浸泡后的晶体活性位点的电子密度显示;c为SARS冠状病毒主蛋白酶与川藏香茶菜丙素浸泡后的晶体活性位点的电子密度显示。图中所示的电子密度,蓝色为2fofc电子密度,为1.0sigma;青色和绿色为1fofc差值电子密度,为2.5sigma图11本发明的川藏香茶菜乙酸乙酯萃取物样品中的存在的小分子的电子密度与SARS冠状病毒主蛋白酶活性位点Cys145残基电子密度形成共价键的显示图。蛋白质分子用白色飘带图表示,活性位点的Cys145残基用球棍模型表示;蓝色电子密度表示Cys145的电子密度,红色电子密度表示未知小分子的电子密度。
图12本发明的SARS冠状病毒主蛋白酶与川藏香茶菜提取物乙酸乙酯萃取样品作用前后的质谱图。蓝色表示SARS冠状病毒主蛋白酶本身的质谱峰;红色表示SARS冠状病毒与川藏香茶菜乙酸乙酯萃取样品作用后的质谱峰。
图13 20μM的川藏香茶菜丙素对SARS冠状病毒主蛋白酶的抑制活性曲线14 2μM的川藏香茶菜丙素对SARS冠状病毒主蛋白酶的抑制活性曲线15 20μM的川藏香茶菜丙素对TGEV主蛋白酶的抑制活性曲线16 2μM的川藏香茶菜丙素对TGEV主蛋白酶的抑制活性曲线17 100μM的川藏香茶菜丙素对AIBV主蛋白酶的抑制活性曲线18 10μM的川藏香茶菜丙素对AIBV主蛋白酶的抑制活性曲线图为了更加详细地解释本发明,下面将结合附图给出本发明的具体实施例。
具体实施例方式
实施例11.SARS冠状病毒主蛋白酶的表达、纯化和晶体生长根据文献(Yang H etal.Proc Natl Acad Sci.2003 Nov;100(23)13190-5)进行SARS冠状病毒主蛋白酶(氨基酸序列见序列表部分)的表达、纯化和晶体生长,得到衍射分辨率较高(高于2.6)的SARS冠状病毒主蛋白酶的晶体。具体地a.SARS冠状病毒主蛋白酶的表达载体构建具体步骤包括(1.1)利用北京华大基因中心提供的编号为BJ01的SARS病毒毒株的cDNA文库,用PCR技术进行体外扩增。
正向引物5′-CGGGATCCAGTGGTTTTAGGAAAATG-3′反向引物5′-CCGCTCGAGTCATTGGAAGGTAACACCAGA-3′(1.2)经PCR扩增的基因片段经BamHI和XhoI双酶酶切后,用琼脂糖凝胶电泳回收大小为1kb左右的片段。
(1.3)将回收片段与T载体进行连接,然后用连接产物转化大肠肝菌Escherichia coli DH5α感受态细胞,并涂在LB平板(含100mg/L氨苄青霉素)上,培养过夜。
(1.4)从平板上挑取多个单克隆,分别接种于装有约5mL的LB的试管(该LB溶液中加入氨卞青霉素,使其终浓度为100mg/L)中,培养过夜。然后用质粒提取试剂盒(博大泰克公司B型质粒小量快速提取试剂盒)提取质粒,并用BamHI和XhoI酶切,然后用琼脂糖凝胶回收大小约为1kb左右的目标基因片段。
(1.5)将目标载体pGEX-4T-1(购自Pharmacia公司)用BamHI和XhoI酶切,然后用琼脂糖凝胶回收酶切的片段。
(1.6)将(1.4)和(1.5)获得的片段连接(将酶切回收后的目标基因片断和目标载体片段按照摩尔数3∶1-6∶1的比率混合,按照Takara DNA Ligation的要求在16℃反应30分钟-18小时),转化大肠肝菌Escherichia coli DH5α感受态细胞,涂在LB平板(含100mg/L氨苄青霉素)上培养过夜。将筛选到的阳性克隆,用于鉴定和测序。
测序结果表明,SARS冠状病毒的主蛋白酶的编码基因已被正确地克隆到pGEX-4T-1载体中。
b.SARS冠状病毒主蛋白酶的表达与纯化包括步骤(2.1)将上述步骤a中得到的含有编码SARS冠状病毒主蛋白酶基因的pGEX-4T-1载体转化Escherichia coli BL21(DE3)的菌株,并用LB平板(含100mg/L氨苄青霉素)筛选阳性克隆。
(2.2)在(2.1)中所述的LB平板上挑取阳性克隆(在含有氨苄青霉素的LB平板上生长出来的单克隆),培养过夜,然后转入1L的LB培养基(含100mg/L氨苄青霉素),当OD600达到0.6-0.8时,加入1mM左右的IPTG,在16℃培养12小时左右;(2.3)5000-8000rpm离心10-15min收集细胞,然后冰浴超声破菌20-30分钟;破菌液13000rpm-15000rpm离心20-40min后收集上清液。
(2.4)将上清液加入PBS预平衡的GST亲和层析柱(GE公司)中,用PBS淋洗20-30个柱床体积去除杂蛋白。最后加入2ml 0.1mg/ml左右的人类鼻病毒3C蛋白酶,在4℃酶切12-20小时,之后收集SARS冠状病毒主蛋白酶。
(2.5)将(2.4)获得的蛋白再用Mono Q(GE公司)阴离子交换层析进行纯化c.SARS冠状病毒的主蛋白酶的晶体生长SARS冠状病毒的主蛋白酶的晶体生长是采用悬滴气相扩散法,在289-297K(K是温度单位开尔文,摄氏温度16℃~24℃)下,用组织培养板进行的。用Hampton Research公司的晶体生长试剂盒(Crystal Screening Kit I & II,PEG/IONKit,PEG 6000 Grid Kit和(NH4)2SO4Kit)对初始结晶条件进行筛选。在PEG6000 Grid Kit的C3(20%PEG 6000,0.1M MES(pH 6.0))条件下的液滴中有红色球状沉淀,并在沉淀的边缘出现微晶。在此结晶条件下,我们尝试了不同pH值的缓冲液(如pH值从5.5到6.5之间进行)、不同聚合度的PEG(如从PEG 3350到PEG 10000之间进行)、不同蛋白浓度(如从2mg/ml到30mg/ml之间进行)和各种添加剂(Hampton Research公司地添加剂试剂盒)等来优化晶体。最后,在1.8~12%的PEG 6000,3%二甲基亚砜(DMSO),1mM DTT,100mM MES缓冲液(pH 5.0~PH6.0)以及蛋白浓度为5~12mg/ml时获得了衍射分辨率较高(高于2.6)的晶体。
2.中药等天然产物混合物备选库的制备源自川藏香茶菜(I.pharicus(Prain)Murata)的天然产物混合物备选样品的制备(1)制备川藏香茶菜提取物浸膏取500克干燥粉碎的川藏香茶菜药材,用4500毫升95%乙醇回流提取,提取3次,每次2小时,合并提取液,浓缩得浸膏;(2)制得川藏香茶菜乙酸乙酯萃取样品称取此浸膏10.5克,用85毫升双蒸水悬浮,倾入1000毫升分液漏斗中,用乙酸乙酯萃取。每次萃取用乙酸乙酯850毫升,充分振摇,静置6小时分层。共萃取6次,合并各次乙酸乙酯萃取物。萃取后得到乙酸乙酯相和水相。乙酸乙酯萃取物用EYELA N1001型旋转蒸发仪蒸除溶剂,得干膏状川藏香茶菜乙酸乙酯萃取样品,备用;(3)制得川藏香茶菜水相大孔树脂50%乙醇洗脱样品萃取后的水相部分用EYELA N1001型旋转蒸发仪蒸除水中溶解的少量乙酸乙酯,用HP-20型大孔树脂柱层析分离(大孔树脂用170毫升,玻璃柱规格30×300mm),依次用水、50%乙醇为流动相洗脱。水洗脱8次,每次洗脱250毫升,弃去水洗部分。50%乙醇洗脱5次,每次洗脱300毫升。将50%乙醇洗脱部分使用EYELA N1001型旋转蒸发仪蒸除溶剂,得干膏状川藏香茶菜水相大孔树脂50%乙醇洗脱样品,备用;(4)获得源自川藏香茶菜的天然产物混合物备选样品上述步骤(2)制得的川藏香茶菜乙酸乙酯萃取样品、步骤(3)制得的川藏香茶菜水相大孔树脂50%乙醇洗脱样品组成了源自初级原料川藏香茶菜的2个备选样品(按照上述方法制备的每一个备选样品都含有多种天然产物小分子化合物,因此称为天然产物混合物备选样品)。
源自其它初级原料的天然产物混合物备选样品的制备由于中药中的清热药往往有抗菌、抗病毒作用,并且非典期间国家中医药管理局公布的三个处方中也大多含此类药(比如处方一药物组成为鱼腥草15克、野菊花6克、茵陈15克、佩兰10克、草果3克。其功效清热解毒、利湿化浊,其中鱼腥草、野菊花为清热药;处方二药物组成为蒲公英15克、金莲花6克、大青叶10克、葛根10克、苏叶6克。其功效清热解毒、散风透邪,其中蒲公英、金莲花、大青叶为清热药;处方三药物组成为芦根15克、金银花10克、连翘10克、薄荷6克、生甘草5克。其功效清热解表、疏风透邪,其中芦根、金银花、连翘为清热药)。故先从这些药入手。查中药学专著(黄兆胜主编,人民卫生出版社,2002年8月第1版)常用的此类药有44种,其中17种归肺经,为鱼腥草,野菊花,大青叶,芦根,金银花,连翘,黄芩,知母,四季青,穿心莲,鸭跖草,栀子,射干,山豆根,马勃,玄参,地骨皮(中国科学院云南昆明植物所提供),列为用于进行SARS冠状病毒主蛋白酶抑制剂筛选的中药等天然产物混合物备选库的初级原料。
另有鄂西香茶菜(中国科学院云南昆明植物所提供)是我国民间用来治疗多种炎症的有效药物,且最近从该属植物中分离出多种被证明具有抗病毒、抗肿瘤等生理活性的二萜类成分,故也列为用于进行SARS冠状病毒主蛋白酶抑制剂筛选的中药等天然产物混合物备选库的初级原料。
源自以上18种初级原料的36个天然产物混合物备选样品按照上述源自川藏香茶菜(I.pharicus(Prain)Murata)的天然产物混合物备选样品的制备方法制备。
按照上述方法制备的38个天然产物混合物备选样品组成了用于进行SARS冠状病毒主蛋白酶抑制剂筛选的中药等天然产物混合物备选库。
3.天然产物混合物备选样品对SARS冠状病毒主蛋白酶抑制活性测定SARS冠状病毒主蛋白酶的活性测定是使用荧光底物MCA-AVLQSGFR-Lys(Dnp)-Lys-NH2(纯度大于95%,上海吉尔生化有限公司)来完成的。该荧光底物的氨基酸序列来源于SARS冠状病毒主蛋白酶的N端自剪切序列。
用于荧光强度测定的仪器为Fluoraskan Ascent荧光仪(ThermoLabsystems,Helsinki,Finland),激发光和发射光的波长分别为320nm和405nm。
按照本发明制备的用于进行SARS冠状病毒主蛋白酶抑制剂筛选的中药等天然产物混合物备选库中的样品称为天然产物混合物备选样品。将步骤2中制备的38个干膏状的天然产物混合物备选样品溶于二甲基亚砜(DMSO)中使其终浓度为50mg/mL,分别制得38个天然产物混合物备选样品的DMSO溶解物。
在缓冲溶液(50mM Tris-HCl(pH 7.3),1mM EDTA(含或不含DTT))中加入SARS冠状病毒主蛋白酶(终浓度0.5μM),加入天然产物混合物备选样品的DMSO溶解物(如川藏香茶菜乙酸乙酯萃取样品的DMSO溶解物)使其终浓度为100μg/mL,底物浓度为20μM,298K放置10分钟后,迅速加入荧光标记底物(MCA-AVLQ↓SGFRL(DNP)L-NH2,终浓度20μM)。激发波长和发射波长分别为320nm和405nm,温度保持298K,每2秒钟记录一次荧光读数。
对照不加入备选样品,其余条件相同。结果示于图1-图9中。
图1-图4以及图6、图8为38个天然产物混合物备选样品对SARS冠状病毒主蛋白酶抑制活性的测定结果,其中每个备选样品的浓度为100μg/mL。
由以上图1-图9可以看出,其中大部分备选样品均对SARS冠状病毒主蛋白酶有抑制活性,其中抑制活性较好的有7个,分别是黄芩、鱼腥草、山豆根、栀子、知母、川藏香茶菜、鄂西香茶菜的乙酸乙酯萃取物样品。
图5、图7和图9是以上7个抑制活性较好的黄芩、鱼腥草、山豆根、栀子、知母、川藏香茶菜、鄂西香茶菜的乙酸乙酯萃取样品的浓度稀释至10μg/mL时的对SARS冠状病毒主蛋白酶抑制活性测定结果,可以发现,将黄芩、鱼腥草、山豆根、栀子、知母、川藏香茶菜、鄂西香茶菜的乙酸乙酯萃取物浓度为10μg/mL时,对SARS冠状病毒主蛋白酶仍有较强的抑制活性。
这就说明在这7种中药的乙酸乙酯萃取样品中可能含有某种成分能够有效的抑制SARS冠状病毒主蛋白酶的活性,也说明这7种中药有可能成为治疗SARS冠状病毒的有效药物。为了了解它们对SARS冠状病毒主蛋白酶抑制作用的确切的分子机理,我们可以选择这7种中药的乙酸乙酯萃取样品进行SARS冠状病毒主蛋白酶的晶体浸泡或共结晶实验。
4、SARS冠状病毒主蛋白酶晶体浸泡川藏香茶菜乙酸乙酯萃取样品由于用于SARS冠状病毒主蛋白酶抑制剂筛选的中药等天然产物混合物备选库中的每一个样品都含有许多种天然产物小分子(因此称为天然产物混合物备选样品),若采用共结晶的方法,对SARS冠状病毒主蛋白酶结晶体系的影响较大,无法获得衍射质量好、数量充足的SARS冠状病毒主蛋白酶的晶体,甚至可能无法生长晶体,因此主要采用晶体浸泡的方法来获得SARS冠状病毒主蛋白酶与小分子抑制剂复合体的晶体。
为了尽可能的在保证晶体质量的基础上,提高浸泡时川藏香茶菜乙酸乙酯样品中小分子化合物的浓度,可以采用如下两种方法制备SARS冠状病毒主蛋白酶的晶体浸泡液。
方法一、将川藏香茶菜乙酸乙酯萃取样品的DMSO溶解物(按照步骤3制备)10倍稀释至SARS冠状病毒主蛋白酶晶体生长池液(1.8~12%的PEG 6000,3%二甲基亚砜(DMSO),1mM DTT,100mM MES缓冲液(pH 5.0~PH6.0))中,13000-15000rpm高速离心,之后取上清液,得浸泡液1,备用。
方法二、将干膏状的川藏香茶菜乙酸乙酯萃取样品(步骤2制备)直接溶解于SARS冠状病毒主蛋白酶晶体生长池液(1.8~12%的PEG 6000,3%二甲基亚砜(DMSO),1mM DTT,100mM MES缓冲液(pH 5.0~PH6.0))中,13000-15000rpm高速离心,之后取上清液,得浸泡液2,备用。
为了保证晶体浸泡时,不至于使川藏香茶菜乙酸乙酯萃取样品中的小分子化合物浓度过高影响晶体的衍射质量,因此将浸泡液1和2分别稀释至原液浓度的10倍,备用。
浸泡时,利用尼龙晶体环等工具,将已生长好的SARS冠状病毒主蛋白酶的晶体从结晶池液中取出,分别加入浸泡液1和2中以及相应的10倍稀释液中,浸泡2-48小时。
5.晶体衍射数据的收集和分析在步骤4的基础之上,使用具备较好衍射能力(分辨率优选大于2.6)的晶体,进行X射线衍射数据的收集和分析。
进行数据收集时,首先使用Hampton Research公司的尼龙晶体环,从浸泡液1和浸泡液2中获取浸泡一定时间(通常为2-48小时)的晶体;并迅速使用Rigaku公司的冷却系统或Oxford Cyrosystem公司的冷却系统,将晶体在以上所述的两种冷冻系统产生的低温氮气流中冷冻至-150℃--180℃;使X射线通过晶体,使用旋进法收集X射线衍射数据。
选择衍射质量好的晶体(分辨率最好大于2.6)进行数据收集、处理,检测是否有小分子的结合。数据收集的统计如下表1表1SARS冠状病毒主蛋白酶母体、SARS冠状病毒主蛋白酶与川藏香茶菜乙酸乙酯萃取物样品浸泡后的复合体数据收集统计

通过观察相同角度的SARS冠状病毒主蛋白酶母体和复合体的活性位点电子密度图(图10)(蓝色是蛋白质本身的电子密度),可以发现,在复合体的活性位点能够观察到清晰的非蛋白质本身的电子密度(青色表示),提示可能有小分子化合物结合在蛋白质分子上。图10a中显示的是SARS冠状病毒主蛋白酶母体结构的活性位点附近的电子密度,图10b为SARS冠状病毒主蛋白酶与川藏香茶菜乙酸乙酯萃取样品浸泡后晶体活性位点附近的电子密度。
仔细观察该处的电子密度(如图10所示),我们可以发现存在于川藏香茶菜乙酸乙酯萃取物样品中的未知小分子的电子密度与SARS冠状病毒主蛋白酶活性位点处Cys145的电子密度紧密结合在一起,提示我们这应该是小分子与Cys145的γ硫原子形成了共价键,从而抑制了SARS冠状病毒主蛋白酶的活性(如图11所示)。
6.结合小分子的分析与鉴定我们将SARS冠状病毒主蛋白酶(浓度为12mg/ml的SARS冠状病毒主蛋白酶水溶液)与川藏香茶菜乙酸乙酯萃取样品按照质量浓度比1∶1-1∶4的比例混合,并在4℃-16℃下作用6-18小时。将作用后的样品,13,000-15,000rpm高速离心,弃去沉淀后备用于质谱分析。质谱采用MALDI(matrix-assistedlaser-desorption ionization)谱对样品进行分析,结果参见图12。
从质谱的结果我们可以发现,在加入川藏香茶菜乙酸乙酯萃取物样品后,分子量出现了明显的偏移,这主要是因为蛋白质分子与川藏香茶菜乙酸乙酯萃取物样品中的小分子结合的原因。通过质谱结果可以计算(用较高的质谱峰值减去较低的质谱峰)发现,结合的小分子分子量约为390Da。
通过分析已有的川藏香茶菜成分,我们发现,在该分子量水平上,主要存在如下所示的化合物
我们可以发现在川藏香茶菜的主要成分中,在该分子量范围内,主要成分以二萜类化合物为主,包括川藏香茶甲、乙、丙、丁、戊、己素(分别对应上面PseurataA、PseurataB、PseurataC、PseurataD、PseurataE、PseurataF),Isodomedin和Dihydropseurata)等。又由于结合在SARS冠状病毒主蛋白酶活性口袋处的小分子,和Cys145的γ硫原子形成了共价键,推知该小分子为Michael受体类小分子,具有α,β共轭不饱和双键。
基于分子量和化学环境的分析,得出结论结合的小分子为川藏香茶菜丙素或者Isodomedin。
为了进一步验证,我们从中国科学院云南昆明植物所获得了川藏香茶菜丙素及Isodomedin的纯化合物,按照以上所述的方法制备川藏香茶菜丙素以及Isodomedin与SARS冠状病毒主蛋白酶复合体的晶体,并进行了数据收集(数据统计结果示于表2)。通过观察电子密度,我们只在川藏香茶菜丙素浸泡过的SARS冠状病毒主蛋白酶的晶体发现了相似的电子密度,如图10.c所示。因此可以确定,从川藏香茶菜乙酸乙酯萃取样品中发现的与SARS冠状病毒主蛋白酶相结合的抑制剂小分子应为川藏香茶菜丙素。
在川藏香茶菜乙酸乙酯萃取物样品和川藏香茶菜丙素浸泡过的SARS冠状病毒主蛋白酶的晶体活性位点,未知电子密度还有所区别,主要是因为不同晶体的数据质量不同,而且在川藏香茶菜乙酸乙酯萃取物样品中川藏香茶菜丙素的相对浓度比纯品肯定要低,因此可能占有率也有所下降。
因此,可以确定结合的小分子即为川藏香茶菜丙素。
表2SARS冠状病毒主蛋白酶母体、SARS冠状病毒主蛋白酶与川藏香茶菜丙素的复合体数据收集统计

由此,我们获得了SARS冠状病毒主蛋白酶与川藏香茶菜丙素复合体的精细三维结构。
这里,我们可以看到本发明的新方法与传统的抑制剂筛选方法相比,能够对筛选效率有极大的提高。按照我们的方法,(1)仅需要筛选一次就可以从川藏香茶菜乙酸乙酯萃取物中快速发现能够与SARS冠状病毒主蛋白酶结合的小分子,(2)该小分子是川藏香茶菜乙酸乙酯萃取物中与SARS冠状病毒主蛋白酶的结合能力最好的小分子,(3)无需在抑制剂筛选时制备众多的纯品化合物;而如果利用传统方法,不仅需要在进行筛选前分离提取制备出至少8种化合物,而且需要至少进行8次的活性测定、并进行至少8次晶体浸泡实验等,才有可能发现结合能力最好的小分子化合物。
7、川藏香茶菜丙素对SARS冠状病毒的主蛋白酶抑制活性的测定我们使用纯品化合物小分子(川藏香茶菜丙素),测定了其对SARS冠状病毒主蛋白酶的抑制活性。川藏香茶菜丙素对SARS冠状病毒主蛋白酶抑制活性的测定按照如下步骤进行在缓冲溶液(50mM Tris-HCl(pH 7.3),1mM EDTA(含或不含DTT))中加入SARS冠状病毒主蛋白酶(终浓度0.5μM),川藏香茶菜丙素(终浓度为20μM),底物浓度为20μM,298K放置10分钟后,迅速加入荧光标记底物(MCA-AVLQ↓SGFRL(DNP)L-NH2,终浓度20μM)。激发波长和发射波长分别为320nm和405nm,温度保持298K,每2秒钟记录一次荧光读数(结果示于图13)。然后改变川藏香茶菜丙素的浓度,在2μM测定其对SARS冠状病毒主蛋白酶的抑制活性(结果示于图14)。对照不加入备选样品,其余条件相同。活性测定结果见图13-图14。从图13-图14中可以发现,20μM的川藏香茶菜丙素对SARS冠状病毒主蛋白酶具有较强的抑制活性,当其浓度为2μM时对SARS冠状病毒主蛋白酶仍然有抑制活性。这就说明中药川藏香茶菜里面含有的川藏香茶菜丙素能够有效的抑制SARS冠状病毒主蛋白酶的活性,也说明川藏香茶菜丙素有可能成为治疗SARS冠状病毒的有效药物。
8、川藏香茶菜丙素对猪传染性胃肠炎病毒(TGEV)的主蛋白酶抑制活性的测定在缓冲液(20mM Tris-HCl pH 7.0,1mM DTT)中,加入所述蛋白酶(0.5μM),川藏香茶菜丙素(20μM),298K放置10分钟后,迅速加入荧光标记底物(MCA-AVLQSGFRL(DNP)L-NH2,20μM)。激发波长和发射波长分别为320nm和405nm,温度保持298K,每2秒钟记录一次荧光读数。然后改变抑制剂浓度,在2μM下测定抑制活性。对照不加入抑制剂,其余条件相同。结果示于图15-图16中。TGEV的表达纯化参考文献Conservation of substratespecificities among coronavirus main proteases.J.Gen.Virol.2002;83(Pt 3)595-9。
9.川藏香茶菜丙素对禽传染性支气管炎病毒(AIBV)的主蛋白酶抑制活性的测定在缓冲液(20mM Tris-HCl pH 7.0,1mM DTT)中,加入所述蛋白酶(1μM),川藏香茶菜丙素(100μM),298K放置10分钟后,迅速加入荧光标记底物(MCA-AVLQSGFRL(DNP)L-NH2,20μM)。激发波长和发射波长分别为320nm和405nm,温度保持298K,每2秒钟记录一次荧光读数。然后改变抑制剂浓度,在10μM下测定抑制活性。对照不加入抑制剂,其余条件相同。结果示于图17-图18中。
AIBV的表达纯化参考文献Preliminary crystallographic analysis of avianinfectious bronchitis virus main protease.Acta.Cryst.F.2007;63(Pt1)24-6。
由图13-图18可以看出,川藏香茶菜丙素对SARS、TGEV、AIBV的主蛋白酶都有抑制活性,其中对SARS冠状病毒主蛋白酶的抑制活性最强。由于TGEV属于冠状病毒属的第一组群,SARS冠状病毒属于冠状病毒属的第二组群,AIBV属于冠状病毒属的第三组群。因此,可以推断川藏香茶菜丙素对整个冠状病毒家族的主蛋白酶都有抑制活性。
本文中所涉及的各种实验用品(包括但不限于化学试剂、生物制品、细胞、生物体、仪器等)之中,对于那些特殊的或不易获得的,文中均已注明了制造商、参考文献或详细的制备方法;未经特别说明的,均为常规实验用品,在本申请申请日之前,可以通过各种方式(例如购买、自行制备等)很方便地获得。
应当理解,在不偏离本发明的精神和范围的情况下,本领域的普通技术人员可以在形式和细节上对其做出各种改变和改进,而这些均被认为落入本发明的保护范围中。
序列表<110>南开大学、清华大学、中国科学院生物物理研究所<120>从中药等天然产物混合物备选库中筛选SARS冠状病毒主蛋白酶抑制剂的新方法<160>2<210>1<211>306<212>prt<213>SARS冠状病毒(SARS Coronavirus)<400>1Ser Gly Phe Arg Lys Met Ala Phe Pro Ser Gly Lys Val Glu Gly Cys1 5 10 15Met Val Gln Val Thr Cys Gly Thr Thr Thr Leu Asn Gly Leu Trp Leu20 25 30Asp Asp Thr Val Tyr Cys Pro Arg His Val Ile Cys Thr Ala Glu Asp35 40 45Met Leu Asn Pro Asn Tyr Glu Asp Leu Leu Ile Arg Lys Ser Asn His50 55 60Ser Phe Leu Val Gln Ala Gly Asn Val Gln Leu Arg Val Ile Gly His65 70 75 80Ser Met Gln Asn Cys Leu Leu Arg Leu Lys Val Asp Thr Ser Asn Pro85 90 95Lys Thr Pro Lys Tyr Lys Phe Val Arg Ile Gln Pro Gly Gln Thr Phe100 105 110Ser Val Leu Ala Cys Tyr Asn Gly Ser Pro Ser Gly Val Tyr Gln Cys
115 120 125Ala Met Arg Pro Asn His Thr Ile Lys Gly Ser Phe Leu Asn Gly Ser130 135 140Cys Gly Ser Val Gly Phe Asn Ile Asp Tyr Asp Cys Val Ser Phe Cys145 150 155 160Tyr Met His His Met Glu Leu Pro Thr Gly Val His Ala Gly Thr Asp165 170 175Leu Glu Gly Lys Phe Tyr Gly Pro Phe Val Asp Arg Gln Thr Ala Gln180 185 190Ala Ala Gly Thr Asp Thr Thr Ile Thr Leu Asn Val Leu Ala Trp Leu195 200 205Tyr Ala Ala Val Ile Asn Gly Asp Arg Trp Phe Leu Asn Arg Phe Thr210 215 220Thr Thr Leu Asn Asp Phe Asn Leu Val Ala Met Lys Tyr Asn Tyr Glu225 230 235240Pro Leu Thr Gln Asp His Val Asp Ile Leu Gly Pro Leu Ser Ala Gln245250255Thr Gly Ile Ala Val Leu Asp Met Cys Ala Ala Leu Lys Glu Leu Leu260 265270Gln Asn Gly Met Asn Gly Arg Thr Ile Leu Gly Ser Thr Ile Leu Glu275 280 285Asp Glu Phe Thr Pro Phe Asp Val Val Arg Gln Cys Ser Gly Val Thr290 295300Phe Gln30权利要求
1.从天然产物混合物备选库中筛选SARS冠状病毒主蛋白酶抑制剂的方法,其特征在于包括步骤A测定天然产物混合物备选库中的备选样品对SARS冠状病毒主蛋白酶的抑制活性,获得对SARS冠状病毒主蛋白酶有抑制活性的“有抑制活性的样品”;将SARS冠状病毒主蛋白酶的晶体与“有抑制活性的样品”接触,得到SARS冠状病毒主蛋白酶与小分子抑制剂复合体的晶体。
2.根据权利要求1所述的方法,其特征在于步骤A之后还包括步骤B收集并分析步骤A中得到的复合体的晶体的X射线衍射数据,获得与SARS冠状病毒主蛋白酶结合的小分子抑制剂的电子密度。
3.根据权利要求2所述的方法,其特征在于步骤B之后还包括步骤C以步骤B中获得的与SARS冠状病毒的主蛋白酶结合的小分子抑制剂的电子密度为基础,鉴定与SARS冠状病毒主蛋白酶结合的小分子抑制剂,获得小分子抑制剂与SARS冠状病毒主蛋白酶复合体的精细三维结构。
4.根据权利要求3所述的方法,其特征在于在步骤C之后还包括步骤D测定小分子抑制剂对SARS冠状病毒主蛋白酶的抑制活性。
5.根据权利要求4所述的方法,其特征在于其中所述的天然产物混合物备选库中的备选样品,是用如下方法制备(a)、取一定量的初级原料,用乙醇溶液回流提取,合并提取液,浓缩得浸膏;(b)、称取一定量步骤(a)中所得的浸膏,用双蒸水悬浮,用有机溶剂萃取,分离有机溶剂相和水相;合并水相,备用;合并每次有机溶剂萃取物,蒸除其中的有机溶剂,得到干膏状有机溶剂相样品,备用;(c)、蒸除步骤(b)中萃取后的水相部分的水中溶解的少量有机溶剂,用大孔树脂柱层析分离,依次用水和乙醇溶液为流动相洗脱;弃去水洗部分,蒸除乙醇溶液洗脱部分的溶剂,得干膏状样品备用;(d)、上述步骤(b)中制得的干膏状有机溶剂相样品和步骤(c)中制得的干膏状样品即为源自一种初级原料的2个备选样品。
6.根据权利要求5所述的方法,其特征在于步骤(a)中乙醇溶液与初级原料的用量质量比为6∶1-12∶1;所用乙醇溶液浓度为70%-95%;步骤(b)中浸膏与悬浮用双蒸水的质量比为1∶5-1∶15,萃取用有机溶剂为乙酸乙酯或氯仿,其体积与悬浮用双蒸水的体积比为6∶1-10∶1;步骤(c)中所用乙醇溶液浓度为50%。
7.根据权利要求1-6中任一项权利要求所述的方法,其特征在于在步骤A中SARS冠状病毒主蛋白酶晶体与“有抑制活性的样品”接触时,使用的方法为晶体浸泡或共结晶的晶体浸泡方法。
8.根据权利要求7中所述的方法,其特征在于在步骤A中SARS冠状病毒主蛋白酶的晶体与“有抑制活性的样品”接触时,可分别与单一的“有抑制活性的样品”接触,也可与混合的“有抑制活性的样品”的接触。
9.根据权利要求8中任一项所述的方法,其中所述天然产物混合物备选库中的备选样品源自中药。
10.根据权利要求9中任一项所述的方法,其中所述天然产物混合物备选库中的备选样品源自下述初级原料鱼腥草,野菊花,大青叶,芦根,金银花,连翘,黄芩,知母,四季青,穿心莲,鸭跖草,栀子,射干,山豆根,马勃,玄参,地骨皮、川藏香茶菜、鄂西香茶菜。
11.根据权利要求1-10中任一项的方法获得的复合体的晶体,其是SARS冠状病毒主蛋白酶-川藏香茶菜丙素复合体的晶体。
12.根据权利要求1-10中任一项所述的方法,其中所述筛选出的小分子抑制剂是川藏香茶菜丙素。
13.川藏香茶菜在制备治疗冠状病毒感染的药物中的用途。
14.川藏香茶菜丙素在制备治疗或者预防冠状病毒感染的药物中的用途。
15.根据权利要求13或14所述的用途,其中所述的冠状病毒为SARS冠状病毒、猪传染性胃肠炎病毒或者禽传染性支气管炎病毒。
全文摘要
本发明提供了从天然产物混合物备选库中筛选SARS冠状病毒主蛋白酶抑制剂的方法将SARS冠状病毒主蛋白酶的晶体与“有抑制活性的样品”接触,获得与SARS冠状病毒主蛋白酶结合的小分子抑制剂的电子密度;以与SARS冠状病毒主蛋白酶结合的小分子抑制剂的电子密度为基础,鉴定小分子抑制剂,从而得到SARS冠状病毒主蛋白酶与小分子抑制剂复合体的精细三维结构。本发明据此方法筛选得到了SARS冠状病毒主蛋白酶的小分子抑制剂川藏香茶菜丙素。川藏香茶菜丙素能够有效抑制冠状病毒主蛋白酶活性,可以治疗或者预防SARS冠状病毒、传染性胃肠炎病毒或者禽传染性支气管炎病毒等冠状病毒感染。
文档编号G01N33/573GK101063234SQ200710065119
公开日2007年10月31日 申请日期2007年4月4日 优先权日2007年4月4日
发明者饶子和, 孙玉娜, 娄智勇, 向科辉 申请人:南开大学, 清华大学, 中国科学院生物物理研究所
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1