CN102750459B - 一种以人乳头瘤病毒hpv的dna解旋酶e1为靶点的抗hpv病毒药物的虚拟筛选方法 - Google Patents
一种以人乳头瘤病毒hpv的dna解旋酶e1为靶点的抗hpv病毒药物的虚拟筛选方法 Download PDFInfo
- Publication number
- CN102750459B CN102750459B CN201210195233.2A CN201210195233A CN102750459B CN 102750459 B CN102750459 B CN 102750459B CN 201210195233 A CN201210195233 A CN 201210195233A CN 102750459 B CN102750459 B CN 102750459B
- Authority
- CN
- China
- Prior art keywords
- hpv
- docking
- active
- screening
- smaller ligand
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 239000003814 drug Substances 0.000 title claims abstract description 34
- 238000000034 method Methods 0.000 title claims abstract description 27
- 241000700605 Viruses Species 0.000 title claims abstract description 26
- 238000003041 virtual screening Methods 0.000 title claims abstract description 26
- 229940079593 drug Drugs 0.000 title claims abstract description 25
- 108090000133 DNA helicases Proteins 0.000 title claims abstract description 11
- 102000003844 DNA helicases Human genes 0.000 title claims abstract description 11
- 208000015181 infectious disease Diseases 0.000 title claims abstract description 6
- 230000002458 infectious effect Effects 0.000 title claims abstract description 6
- 239000003446 ligand Substances 0.000 claims abstract description 43
- 238000003032 molecular docking Methods 0.000 claims abstract description 42
- 238000012216 screening Methods 0.000 claims abstract description 29
- 101710185720 Putative ethidium bromide resistance protein Proteins 0.000 claims abstract description 4
- 208000019065 cervical carcinoma Diseases 0.000 claims description 10
- 230000003993 interaction Effects 0.000 claims description 8
- 125000001493 tyrosinyl group Chemical group [H]OC1=C([H])C([H])=C(C([H])=C1[H])C([H])([H])C([H])(N([H])[H])C(*)=O 0.000 claims description 8
- 206010059313 Anogenital warts Diseases 0.000 claims description 6
- 231100000041 toxicology testing Toxicity 0.000 claims description 6
- 239000013078 crystal Substances 0.000 claims description 4
- 238000011156 evaluation Methods 0.000 claims description 4
- 238000004513 sizing Methods 0.000 claims description 3
- 230000002265 prevention Effects 0.000 claims description 2
- DWNBOPVKNPVNQG-LURJTMIESA-N (2s)-4-hydroxy-2-(propylamino)butanoic acid Chemical group CCCN[C@H](C(O)=O)CCO DWNBOPVKNPVNQG-LURJTMIESA-N 0.000 claims 1
- 238000002560 therapeutic procedure Methods 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 19
- 150000001875 compounds Chemical class 0.000 abstract description 16
- 238000011160 research Methods 0.000 abstract description 9
- 206010008342 Cervix carcinoma Diseases 0.000 abstract description 5
- 208000006105 Uterine Cervical Neoplasms Diseases 0.000 abstract description 5
- 201000010881 cervical cancer Diseases 0.000 abstract description 5
- 241001465754 Metazoa Species 0.000 abstract description 3
- 125000000539 amino acid group Chemical group 0.000 abstract description 3
- 239000002547 new drug Substances 0.000 abstract description 3
- 208000022361 Human papillomavirus infectious disease Diseases 0.000 description 41
- QTBSBXVTEAMEQO-UHFFFAOYSA-N Acetic acid Chemical compound CC(O)=O QTBSBXVTEAMEQO-UHFFFAOYSA-N 0.000 description 18
- 239000013641 positive control Substances 0.000 description 13
- 102000004190 Enzymes Human genes 0.000 description 9
- 108090000790 Enzymes Proteins 0.000 description 9
- 108020004414 DNA Proteins 0.000 description 8
- -1 phenylacetic acid compound Chemical class 0.000 description 7
- 230000029812 viral genome replication Effects 0.000 description 7
- 208000000260 Warts Diseases 0.000 description 6
- 235000011054 acetic acid Nutrition 0.000 description 6
- 229960000583 acetic acid Drugs 0.000 description 6
- 102000004169 proteins and genes Human genes 0.000 description 6
- 108090000623 proteins and genes Proteins 0.000 description 6
- 201000010153 skin papilloma Diseases 0.000 description 6
- BMSPUAFKNKGVHS-QHCPKHFHSA-N 2-[4-[3-[2-[[(1s)-2,3-dihydro-1h-inden-1-yl]amino]-2-oxoethyl]phenyl]phenyl]sulfonylacetic acid Chemical compound C1=CC(S(=O)(=O)CC(=O)O)=CC=C1C1=CC=CC(CC(=O)N[C@@H]2C3=CC=CC=C3CC2)=C1 BMSPUAFKNKGVHS-QHCPKHFHSA-N 0.000 description 5
- 108060004795 Methyltransferase Proteins 0.000 description 5
- 238000007621 cluster analysis Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 230000000840 anti-viral effect Effects 0.000 description 4
- 238000007877 drug screening Methods 0.000 description 4
- 230000002255 enzymatic effect Effects 0.000 description 4
- 150000002611 lead compounds Chemical class 0.000 description 4
- 229930014626 natural product Natural products 0.000 description 4
- 235000018102 proteins Nutrition 0.000 description 4
- 231100000419 toxicity Toxicity 0.000 description 4
- 230000001988 toxicity Effects 0.000 description 4
- VYNUATGQEAAPAQ-UHFFFAOYSA-N 2-sulfonylacetic acid Chemical compound OC(=O)C=S(=O)=O VYNUATGQEAAPAQ-UHFFFAOYSA-N 0.000 description 3
- 206010028980 Neoplasm Diseases 0.000 description 3
- 241001631646 Papillomaviridae Species 0.000 description 3
- 201000011510 cancer Diseases 0.000 description 3
- 201000010099 disease Diseases 0.000 description 3
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 3
- 229940000406 drug candidate Drugs 0.000 description 3
- 239000000284 extract Substances 0.000 description 3
- 238000000338 in vitro Methods 0.000 description 3
- 238000002360 preparation method Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000010076 replication Effects 0.000 description 3
- 239000000725 suspension Substances 0.000 description 3
- 230000035502 ADME Effects 0.000 description 2
- 230000004543 DNA replication Effects 0.000 description 2
- 241000701806 Human papillomavirus Species 0.000 description 2
- 101000775744 Human papillomavirus type 18 Replication protein E1 Proteins 0.000 description 2
- 208000009608 Papillomavirus Infections Diseases 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 238000012790 confirmation Methods 0.000 description 2
- ZUOUZKKEUPVFJK-UHFFFAOYSA-N diphenyl Chemical compound C1=CC=CC=C1C1=CC=CC=C1 ZUOUZKKEUPVFJK-UHFFFAOYSA-N 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 229910052739 hydrogen Inorganic materials 0.000 description 2
- 239000001257 hydrogen Substances 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 239000002504 physiological saline solution Substances 0.000 description 2
- 238000000547 structure data Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 229940121819 ATPase inhibitor Drugs 0.000 description 1
- 101710153593 Albumin A Proteins 0.000 description 1
- 102000009027 Albumins Human genes 0.000 description 1
- 108010088751 Albumins Proteins 0.000 description 1
- 239000004475 Arginine Substances 0.000 description 1
- 101100001675 Emericella variicolor andJ gene Proteins 0.000 description 1
- NYHBQMYGNKIUIF-UUOKFMHZSA-N Guanosine Chemical compound C1=NC=2C(=O)NC(N)=NC=2N1[C@@H]1O[C@H](CO)[C@@H](O)[C@H]1O NYHBQMYGNKIUIF-UUOKFMHZSA-N 0.000 description 1
- 241000341655 Human papillomavirus type 16 Species 0.000 description 1
- OUYCCCASQSFEME-QMMMGPOBSA-N L-tyrosine Chemical compound OC(=O)[C@@H](N)CC1=CC=C(O)C=C1 OUYCCCASQSFEME-QMMMGPOBSA-N 0.000 description 1
- 108010075285 Nucleoside-Triphosphatase Proteins 0.000 description 1
- 102000008021 Nucleoside-Triphosphatase Human genes 0.000 description 1
- 108700026244 Open Reading Frames Proteins 0.000 description 1
- 108020005091 Replication Origin Proteins 0.000 description 1
- 108020005202 Viral DNA Proteins 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 238000007259 addition reaction Methods 0.000 description 1
- 239000000362 adenosine triphosphatase inhibitor Substances 0.000 description 1
- 230000003281 allosteric effect Effects 0.000 description 1
- 125000003275 alpha amino acid group Chemical group 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- ODKSFYDXXFIFQN-UHFFFAOYSA-N arginine Natural products OC(=O)C(N)CCCNC(N)=N ODKSFYDXXFIFQN-UHFFFAOYSA-N 0.000 description 1
- 125000004429 atom Chemical group 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- WLJVXDMOQOGPHL-UHFFFAOYSA-N benzyl-alpha-carboxylic acid Natural products OC(=O)CC1=CC=CC=C1 WLJVXDMOQOGPHL-UHFFFAOYSA-N 0.000 description 1
- 230000000975 bioactive effect Effects 0.000 description 1
- 235000010290 biphenyl Nutrition 0.000 description 1
- 239000004305 biphenyl Substances 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000006555 catalytic reaction Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000000795 chinese herbal drug Substances 0.000 description 1
- 229940125782 compound 2 Drugs 0.000 description 1
- 229940126214 compound 3 Drugs 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 210000005002 female reproductive tract Anatomy 0.000 description 1
- 235000013305 food Nutrition 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 208000021145 human papilloma virus infection Diseases 0.000 description 1
- 125000004435 hydrogen atom Chemical class [H]* 0.000 description 1
- 230000007062 hydrolysis Effects 0.000 description 1
- 238000006460 hydrolysis reaction Methods 0.000 description 1
- 239000003112 inhibitor Substances 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 150000002632 lipids Chemical class 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000010534 mechanism of action Effects 0.000 description 1
- 230000004630 mental health Effects 0.000 description 1
- 239000000178 monomer Substances 0.000 description 1
- 230000007524 negative regulation of DNA replication Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 150000002894 organic compounds Chemical class 0.000 description 1
- 230000001717 pathogenic effect Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 229960003424 phenylacetic acid Drugs 0.000 description 1
- 239000003279 phenylacetic acid Substances 0.000 description 1
- 230000003334 potential effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000003753 real-time PCR Methods 0.000 description 1
- 230000003362 replicative effect Effects 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 239000002342 ribonucleoside Substances 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000002904 solvent Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006277 sulfonation reaction Methods 0.000 description 1
- 229940043263 traditional drug Drugs 0.000 description 1
- OUYCCCASQSFEME-UHFFFAOYSA-N tyrosine Natural products OC(=O)C(N)CC1=CC=C(O)C=C1 OUYCCCASQSFEME-UHFFFAOYSA-N 0.000 description 1
- 241000701447 unidentified baculovirus Species 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 230000009385 viral infection Effects 0.000 description 1
- 210000003905 vulva Anatomy 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Landscapes
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
- Investigating Or Analysing Biological Materials (AREA)
Abstract
本发明涉及一种以人乳头瘤病毒HPV的DNA解旋酶E1为靶点的抗HPV病毒药物的虚拟筛选方法。包括以下步骤：1）确定E1蛋白结构的PDB结构数据；2）采用分子对接软件，依据E1的活性氨基酸残基位点确定活性中心，设定活性口袋；3）对小分子配体进行筛选；4）建立对接用小分子配体数据库；5）根据设定的活性口袋，利用分子对接软件，将对接用小分子配体数据库中的小分子配体与活性口袋一一做对接；6）根据对接结果的排序，初步确定具有抗宫颈癌效果的先导药物。采用本发明可在短时间内获得活性化合物的线索，然后再利用动物水平筛选或者高通量细胞平台筛选具有抗宫颈癌药物，大大提高了速度和效率，缩短新药研究周期。
Description
技术领域
本发明涉及药物筛选方法，尤其涉及以HPVE1为靶点、利用分子对接进行的新型抗HPV病毒药物的虚拟筛选方法。
背景技术
宫颈癌是严重危害妇女健康和生命的疾病，每年约有50万新发病例，25万人死亡，是引起妇女死亡的第二大原因，其中大部分在发展中国家，占80%左右。宫颈癌是我国女性生殖道恶性肿瘤中发病率最高的一类肿瘤。据统计，我国每年新发宫颈癌病例约13.5万人，其中约5万人因此丧生，使广大妇女的身心健康和生命安全受到严重威胁。1977年德国学者ZurHausen等从宫颈癌标本中发现了人乳头状病毒（humanpapillomavirus，HPV）DNA，推测HPV感染与宫颈癌发生有关(ZurH.atel.,CancerRes,1976,36:530.)。进一步的分子流行病学证据表明，HPV与人类多种癌症，特别是与宫颈癌有着密切的关系。由国际癌症研究机构（IARC）组织在22个国家进行了针对引起侵袭性宫颈癌（ICC）的HPV类型的调查。组织学确诊为ICC的1000个病例中有99.7%发现HPVDNA阳性，其中HPV的类型主要为HPV16（53%）和HPV18（15%）。由于HPV与多种癌症有关，对人类危害极大，预计抗HPV病毒药物将成为今后的研发和市场热点，受到更多关注。
传统的药物筛选无论是动物水平筛选还是高通量细胞平台筛选，都存在着天然药物及食物成分复杂，从确认提取物有潜在功效到确定有效成分，到最终分离出单体成分需要花费大量的人力、时间和精力，效率极低。因此，近年来虚拟筛选方法发展和在新药研究中的应用受到了广泛的重视。于是建立一套高效、有效、速效抗HPV病毒药物和防治宫颈癌药物筛选体系成为当务之急。
发明内容
本发明的目的是建立一种新颖，快速，高效的一种以人乳头瘤病毒HPV的DNA解旋酶E1为靶点的抗HPV病毒药物的虚拟筛选方法，以丰富现有的筛选方法。
本发明利用先进的计算机虚拟筛选技术，结合HPV药物体外细胞模型筛选方法，旨在发现一批能有效结合HPVE1蛋白（HPV基因组编码的蛋白中唯一的酶类蛋白，具有ATP依赖的DNA解旋酶活性）的活性部位从而抑制其ATP依赖的解旋酶活性的药物，从而阻断HPV病毒的复制进而达到预防以及治疗宫颈癌以及其他HPV病毒感染所致疾病的目的。
本发明，利用计算机进行药物虚拟筛选最关键的是基础理论研究的可靠性以及相关靶点的选择。所有的乳头瘤病毒都含有一个小的环状的双链只能编码八个特异性蛋白的DNA（如图1所示）。众多的高度保守的蛋白中，只有E1解旋酶具有酶活性。E1蛋白约68kDa，在HPV阳性细胞中低水平表达，其主要的生化特性是具有ATP酶和3’→5’解旋酶活性，能够识别HPV的复制起点。E1蛋白识别并结合到富含AT的HPV复制起点序列，在E2蛋白的配合下，解除HPV的DNA超螺旋，让复制复合物进入并启动DNA复制。因为E1是病毒复制和发病机制的根本，所以在抗HPV病毒药物的研究中，E1是具有较大吸引力的靶点。这一点在白尾野兔乳头瘤病毒模型中已经得到验证(Wu,X.,W.Xiao,andJ.L.Brandsma.1994;J.Virol.68:6097–6102)。
E1是乳头瘤病毒的复制解旋酶，与E2协同作用后结合到复制起始点。两者的结合不仅仅影响了两个蛋白结合到特殊DNA元件的能力，也包括两者在E2N端激活区和E1ATP酶/解旋酶功能域的关联作用。最近报道了一系列的结合到E2激活区的和阻止与E1关联的HPVDNA复制的抑制剂。E1-E2的复杂的结合能推进E1分子复合物聚集到起点从而形成六聚物，当ATP结合时发生应激反应。这些六聚物能催化E1的激活型从而在复制叉前解旋病毒DNA。作为众多解旋酶中的一例，E1的DNA解旋活性能被ATP的水解产生动力。最近的通过杆状病毒群表达模型做出了一些具有酶活性的特异性的高纯度的HPV6-11型E1蛋白。这些研究表明了HPV6-11中的E1蛋白对于ATP有相似的KM值，结合到E2上的半个E1的C端组成了核苷三磷酸酶和解旋活性。更有趣的是，E2通过提高其ATP的KM值七倍进而阻碍E1的核苷三磷酸酶活性。但是，也发现ATP能通过E1-E2之间蛋白相互作用而影响到两者结合到起点的能力。这一结果表明ATP不仅仅能激活E1的解旋酶活性，也能够作为一个在引起DNA复制的过程中帮助分割E1-E2关联作用进而在起始位点形成六聚物的一个分子感应器。进一步的研究发现，位于589位的精氨酸指（ArginineFinger）在E1发挥ATP酶活性上具有重要作用，而且R589在HPV6,11,16,18中高度保守。另外，现在已经发现，BiphenylsulfonaceticAcid的一系列衍生物可以通过作用于HPV6E1蛋白486位上的酪氨酸（Y486）而变构性地发挥抑制其ATP酶活性进而抑制其解旋酶活性的作用，而且该类化合物对于在E1蛋白具有高度保守的Y486的HPV18也能发挥显著的ATP酶抑制活性（HPV6的E1蛋白保守残基Y486在HPV18为Y492），提示Y486及其附近区域可能是ATP酶抑制剂发挥生物活性的活性中心（ATP结合位点为Lys482，HPV18E1为Lys490）。在了解以上的作用机理后，我们运用分子对接方法，以HPV18E1与E2蛋白复合物的六聚体的三维立体结构数据（PDBID:1TUE）为基准，分别以ATP酶功能域中高度保守的活性残基R589以及Y492为筛选的活性中心，从天然产物数据库和中草药数据库中筛选出具有E1蛋白ATP酶抑制活性，从而抑制HPV病毒复制的候选药物，并且以HPV感染HaCAT人永生化角质细胞为模型确认候选药物的抗HPV复制的作用，产生了一系列新型抗HPV病毒先导化合物，预期在防治宫颈癌以及其他HPV感染疾病方面发挥重要作用。
基于上述机理，本发明的构思如下：本发明中，选用国际天然产物数据库(DunkelM,atel.,NucleicAcidsRes2006;34(Databaseissue):D678-D683)及中国天然产物数据库，首先利用internet蛋白质三级结构数据库中下载HPV18E1ATP酶的蛋白晶体三级结构数据，该三级结构中含有HPV18E1蛋白的ATP酶功能域。以此为基础，利用国际上应用最为广泛的分子对接软件Autodock4.2软件，以E1蛋白的Y492为活性中心构建活性口袋，通过虚拟筛选和精确对接两种方法，从几万种天然产物中筛选出与该活性口袋对接较好的分子，并与联苯磺化乙酸类化合物作为对接分子的阳性对照，确定出一批具有抗HPV活性的先导化合物，再以人永生化上表皮细胞株HaCaT细胞的细胞平台的高通量筛选进行进一步验证。筛选结果初步证实，本方法丰富了抗HPV病毒药物筛选技术，具有快速、高效的特点。
本发明提供的一种以人乳头瘤病毒HPV的DNA解旋酶E1为靶点的抗HPV病毒药物的虚拟筛选方法，包括以下步骤：
1)确定E1蛋白结构的PDB结构数据；
2)采用分子对接软件，依据E1的活性氨基酸残基位点确定活性中心，设定活性口袋；
3)对小分子配体进行筛选，筛选项目为毒性筛选、类药性筛选；
4)建立对接用小分子配体数据库；
5)根据设定的活性口袋，利用分子对接软件，将对接用小分子配体数据库中的小分子配体与活性口袋一一做对接；
6)根据对接结果的排序，小分子配体与活性中心的相互作用方式，进行综合评价，初步确定具有抗HPV病毒效果的先导药物。
上述的虚拟筛选方法：步骤3）置于步骤1）至步骤6）任何一个步骤之前。即：可以对尚未进行分子对接或者已经对接后的小分子数据库中的小分子配体进行毒性、类药性情况进行筛选，而不受排序限制。
上述虚拟筛选方法，所述的E1蛋白的分子对接活性中心，可以是E1蛋白做为DNA解旋酶的酶活性中心。
上述虚拟筛选方法，所述的E1蛋白的分子对接活性中心，可以是E1蛋白做为DNA解旋酶的变构中心。
上述虚拟筛选方法，所述的E1蛋白的分子对接活性中心，可以是E1蛋白做为DNA解旋酶的影响酶活性的其他氨基酸残基。
本发明的有益效果是：本发明通过药物的虚拟筛选，可以在短时间内获得活性化合物的线索，将研究目标从几百万个化合物集中到几百个化合物，而且，虚拟筛选的阳性率（5%～30%）远远高于高通量实验筛选（0.01%～0.1%）。然后再从虚拟筛选后的先导化合物中，利用动物水平筛选或者高通量细胞平台，筛选具有抗宫颈癌药物，因此大大提高了筛选化合物的速度和效率，缩短新药研究的周期。
附图说明
图1HPV病毒的基因组结构和开放阅读框架（以HPV11为例）。
图2HPV18E1的蛋白质一级结构序列。
图3E1的活性功能域的蛋白质三级结构模式图。
图4小分子配体对接程序中使用的SHELL脚本。
图5A化合物2-(4-phenylphenyl)sulfonylacetat的结构式。
图5B化合物2({4'-[2-oxo-2-(phenylamino)ethyl]biphenyl-4-yl}sulfonyl)aceticacid的结构式。
图5C化合物3({3'-[(2-methylbenzyl)carbamoyl]biphenyl-4-yl}sulfonyl)aceticacid的结构式。
图6A阳性对照分子2-(4-phenylphenyl)sulfonylacetate与E1活性中心的成簇分析及相互作用图。
图6B阳性对照分子2({4'-[2-oxo-2-(phenylamino)ethyl]Biphenyl-4-yl}sulfonyl)aceticacid与E1活性中心的成簇分析及相互作用图。
图6C阳性对照分子3({3'-[(2-methylbenzyl)carbamoyl]biphenyl-4-yl}sulfonyl)aceticacid与E1活性中心的成簇分析及相互作用图。
图7A[(3'-{2-[(1S)-2,3-dihydro-1H-inden-1-ylamino]-2-oxoethyl}biphenyl-4-yl)sulfonyl]aceticacid的结构式。
图7B[(3'-{2-[(1S)-2,3-dihydro-1H-inden-1-ylamino]-2-oxoethyl}biphenyl-4-yl)sulfonyl]aceticacid与E1活性中心的成簇分析及相互作用图。
图8[(3'-{2-[(1S)-2,3-dihydro-1H-inden-1-ylamino]-2-oxoethyl}biphenyl-4-yl)sulfonyl]aceticacid对HPV-DNA复制抑制效果检测。
具体实施方式
以下实施例将有助于本领域的普通技术人员进一步理解本发明，但不以任何形式限制本发明。
实施例1以HPV18E1蛋白的酪氨酸残基Y492为活性位点，进行抗宫颈癌药物的虚拟筛选。
（一）虚拟筛选步骤如下：
1、ProteinDataBank中HPV18E1蛋白三维晶体结构的获得
利用ProteinDataBank（PDB）数据库，从中将我们实验所需的HPV18E1蛋白三维晶体结构（图3）下载下来，并获得E1蛋白的氨基酸序列（FASTA格式，图2），从而进行下步运算。
2、采用Autodock分子对接软件，依据E1的活性酪氨酸残基位点Y492确定活性中心，设定活性口袋
本步骤中，首先利用已知的可与E1的活性酪氨酸残基位点Y492发生加成反应的分子做阳性对照分子，依据其在以活性位点Y492为中心构建的活性口袋中与口袋相互作用的结合能及簇的分布情况为标准，确定E1蛋白中以Y492为中心的活性口袋的大小。
2.1利用Autodock分子对接软件进行分子对接计算的程序设置
主要步骤如下：
1)将步骤1中，E1活性功能区域的蛋白质三级结构在MGLTools中进行处理。主要包括去除水分子，加非极性氢，加上Gasteiger电荷。最后保存为1TUEforMD.pdbqt。
2)计算格点能量(AutoGrid)，在Grid>SetMapTypes>Directly中把要进行格点能量计算的原子类型输入进去。运行AutoGrid，完成之后会生成一系列文件，主要为每个原子的map文件。
3)为每个小分子配体准备对接所需的文件夹，这个过程要使用SHELL脚本实现（图4）。先制作一个小分子配体的列表文件ligands.list，然后按照图4所示的SHELL脚本的命令依次输入命令运行脚本a.csh；完成上面的操作后就会为每个小分子配体生成一个文件夹，该文件夹包含对接所需的所有文件：配体、受体、AutoGrid产生的map文件、dpf文件等。下一步即可以对接。
4)循环运行AutoDock：循环运行AutoDock，对接每一个小分子配体，并提取结合能等对接结果并对对接结果排序。这个脚本可以同时运行多个，互不干扰，因此可以完全利用计算机资源。
5)提取对接结果：最后得到summary.sort：虚拟筛选至此做完。
2.2以E1的活性酪氨酸残基位点Y492为活性中心，设定活性口袋
以E1蛋白的Y492为活性中心，以三种已知与Y492活性中心发生相互作用的化合物小分子，(1)2-(4-phenylphenyl)sulfonylacetate；（2）2({4'-[2-oxo-2-(phenylamino)ethyl]biphenyl-4-yl}sulfonyl)aceticacid；(3)3({3'-[(2-methylbenzyl)carbamoyl]biphenyl-4-yl}sulfonyl)aceticacid；为阳性对照分子，三个阳性对照分子的结构式分别如图5A、图5B和图5C所示。调整活性口袋的大小，进行分子对接，从而确定活性口袋的大小。确定活性口袋大小原则为综合考察评价使阳性对照分子与E1蛋白的相互作用的结合能和构象簇的集中情况。
对接采用计算机广泛应用的Lamarckian遗传算法（LGA），每个配体产生50个对接构象，ga_pop_size设为300ga_num_evals设为250000，其余参数采用程序默认值。对接时小分子配体可以自由转动，蛋白质分子则被固定。计算结束后以RMSD值作为聚类分析的依据。
根据反复验证，我们最终确定了活性口袋的设置条件。具体如下：以E1的活性酪氨酸残基位点Y492为活性中心的活性口袋：中心坐标为x为28.807、y为32.339、z为125.869，网格大小值为64×50×44，步长（spacing）为0.375nm。
下面是在确定的这一活性口袋下，三个阳性对照分子在活性口袋中与活性中心及周围基团相互作用的分析。
作用图截取簇中构象数大且结合能也较高的进行分析，结果如图6所示。
通过观察这三个参照分子的对接结果可得出如下结论：bindingenergy值较高，大概处于-7左右，结合能较强，配体分子能较好的与受体分子结合。通过以上阳性对照，确定了上述以E1蛋白中Y492为活性中心的活性中心及活性口袋设置的合理性，可以进行下一步的小分子配体筛选。
3、对小分子配体进行筛选，筛选项目为毒性筛选、类药性筛选
在正式对接前对小分子配体数据库中的小分子配体进行筛选。筛选项目为毒性筛选、类药性筛选。
小分子配体数据库选定免费的国际天然产物数据库（http://bioinformatics.charite.de/supernatural/），该数据库中含有45917种目前国际上已经报道过的并且是容易从公司购买得到的天然有机化合物的结构数据。
3.1毒性筛选：运用化合物毒性预测软件Osiris对小分子配体数据库中的小分子配体进行毒性预测，淘汰毒性预测结果显示阳性的分子。
3.2类药性筛选：运用ADME（药代动力学）预测软件或相关软件的ADME预测功能计算小分子配体数据库中各个小分子配体的类药性指数，指标包括多项，比如：分子量、氢键受体、氢键给体、可旋转键数、脂水分配系数、水溶性指数、分子极性表面积、分子体积、基于经验的药物类似性指数、化学焓变指数（分子稳定性指数）。依据各指数反映的分子类药性优劣，再给予各指数以优劣评分，指数最优者为1，指数最劣者为0。将每个分子的10项指数优劣性评分加和，得出总类药性评分。
4、建立对接用小分子配体数据库
综合步骤3中的毒性筛选结果，类药性评估结果的小分子配体的搜索结果，建立适用于本筛选目的对接用小分子配体数据库。
5、根据设定的活性口袋，利用分子对接软件，将对接用小分子配体数据库中的小分子配体与活性口袋一一做对接
根据上述步骤确定的活性口袋及对接用小分子配体数据库，运行2.1中设置的分子对接程序。将小分子配体与活性口袋一一做对接，确定先导药物。分为大批量筛选与精确筛选两个步骤。
大批量筛选：设置Autodock中的参数ga_num_evals=2500000及ga_run=10，运行2.1中设置的分子对接程序，进行一一对接。
当然也可依据构建的或选定的其它小分子配体数据库的大小进行ga_num_evals及ga_run的设置。
精确筛选：将Autodock软件中的ga_num_evals和ga_run参数设置到最高，即ga_num_evals=250000，ga_run=50，运行2.1中设置的分子对接程序，进行精确筛选。
6、根据对接结果的排序，初步确定具有抗HPV病毒复制效果的先导药物
经过第5步的精确筛选，获得对接后小分子配体数据库，根据对接结果的排序，初步确定具有抗HPV病毒复制效果的先导药物。
当然，为了进一步缩小先导药物的个数还可以采用：综合评价对接结果的排序，小分子配体药效团含有情况以及小分子配体与活性中心的相互作用方式，来确定具有抗HPV病毒复制效果的先导药物。
经过以上筛选，初步确定与E1蛋白Y492活性中心结合的结合能在-7左右的具有较强结合能力的先导化合物200余个，其中包含与阳性对照分子具有相同母核的化合物，这为进一步的抗HPV病毒复制药物活性筛选极大地缩小了范围。
（二）利用细胞筛选平台筛选具有抗HPV病毒复制效果的先导药物。
通过上述的虚拟筛选，从免费的天然产物数据中筛选出一批与E1蛋白的Y492为中心构建的活性口袋相互作用较好的候选化合物。然后根据从临床分离HPV6病毒并培养，加入药物后，进一步验证候选化合物的抗病毒效果。
我们从虚拟筛选后的候选化合物中，挑选结合能在-7左右和作用方式与阳性对照药类似的候选药物，[(3'-{2-[(1S)-2,3-dihydro-1H-inden-1-ylamino]-2-oxoethyl}biphenyl-4-yl)sulfonyl]aceticacid，分子结构式如图7A所示，其与E1蛋白的Y492活性中心相互作用情况如附图7B所示。选用该化合物进行如下体外抗病毒实验。
人乳头瘤病毒感染是引起尖锐湿疣(CA)的主要病因，而且主要是由HPV6/11型感染所致。从临床皮肤科采集来自于尖锐湿疣患者外阴部位的疣体标本(经鉴定为HPV6/11型)，集中到一只离心管中，称重30mg，用玻璃匀浆器匀浆，再加入5ml生理盐水，配成6mg/ml的离体CA混悬液，备用。吸取25μl样品液和25μlCA混悬液至115ml离心管中，振摇均匀，使样品与疣体能充分作用，置37℃水浴箱培养，另设生理盐水阳性对照、各提取溶媒空白对照。共设置四组，在样品处理疣体3、5、7天后，各取一组进行以下试验。从温浴箱中取出经样品处理的CA混悬液，提取DNA，作为模板进行FQ-PCR（定量PCR）扩增，以检验DNA含量。结果如附图8所示，在样品处理5天后，候选物40μM时，可以显著抑制HPV-DNA的复制，处理7天时PCR检测结果呈阴性。结论：确定[(3'-{2-[(1S)-2,3-dihydro-1H-inden-1-ylamino]-2-oxoethyl}biphenyl-4-yl)sulfonyl]aceticacid为抗HPV病毒药物候选物，可进一步用于动物实验。
通过虚拟筛选出的其它候选化合物，都可以按照上述方法，继续进行体外抗病毒实验，从而确定其抗病毒效果。
实施例2一种以E1蛋白的R589为活性位点构建活性口袋，进行的抗HPV病毒药物的虚拟筛选。
方法与实施例1相同，不同点在于：设置活性口袋时：以R589的中心坐标为活性口袋的中心，最终确定活性口袋的设置条件如下：以E1的活性酪氨酸残基位点R589为活性中心的活性口袋：中心坐标为x为36.410，y为25.423，z为106.056，网格大小值为40×40×40，步长（spacing）为0.375nm。
Claims (3)
1.一种以人乳头瘤病毒HPV的DNA解旋酶E1为靶点的抗HPV病毒药物的虚拟筛选方法，其特征在于，包括以下步骤：
1）确定E1蛋白结构的PDB晶体结构数据；
2）采用Autodock分子对接软件，依据E1的活性氨酸残基位点确定分子对接的活性中心，设定活性口袋：当依据E1的活性酪氨酸残基位点R589确定分子对接的活性中心，设定活性口袋时，活性口袋的设置条件如下：以E1的活性酪氨酸残基位点R589为活性中心的活性口袋：中心坐标为x为36.410，y为25.423，z为106.056，网格大小值为40×40×40，步长为0.375nm；
3）对小分子配体进行筛选，筛选项目为毒性筛选、类药性筛选；
4）建立对接用小分子配体数据库；
5）根据设定的活性口袋，利用分子对接软件，将对接用小分子配体数据库中的小分子配体与活性口袋一一做对接；
6）根据对接结果的排序，小分子配体与活性中心的相互作用方式，进行综合评价，初步确定具有抗HPV病毒效果的先导药物。
2.按照权利要求1所述虚拟筛选方法，其特征在于：步骤3）置于步骤1）至步骤6）任何一个步骤之前。
3.按照权利要求1或2所述的虚拟筛选方法筛选出的抗HPV病毒药物在制备预防和治疗宫颈癌和尖锐湿疣药物中的应用。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210195233.2A CN102750459B (zh) | 2012-06-14 | 2012-06-14 | 一种以人乳头瘤病毒hpv的dna解旋酶e1为靶点的抗hpv病毒药物的虚拟筛选方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210195233.2A CN102750459B (zh) | 2012-06-14 | 2012-06-14 | 一种以人乳头瘤病毒hpv的dna解旋酶e1为靶点的抗hpv病毒药物的虚拟筛选方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102750459A CN102750459A (zh) | 2012-10-24 |
| CN102750459B true CN102750459B (zh) | 2016-04-13 |
Family
ID=47030635
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210195233.2A Expired - Fee Related CN102750459B (zh) | 2012-06-14 | 2012-06-14 | 一种以人乳头瘤病毒hpv的dna解旋酶e1为靶点的抗hpv病毒药物的虚拟筛选方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102750459B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104504301A (zh) * | 2014-12-31 | 2015-04-08 | 贵州大学 | 一种以烟草花叶病毒rna解旋酶为靶标的抗病毒药物虚拟筛选方法 |
| CN105608341B (zh) * | 2015-12-19 | 2018-01-19 | 中国人民解放军国防科学技术大学 | 基于动态数据划分和冲突消解的高通量药物虚拟筛选方法 |
| CN106909785B (zh) * | 2017-02-25 | 2019-08-09 | 辽宁大学 | 一种以24-脱氢胆固醇还原酶(dhcr24)为靶点进行降胆固醇药物的虚拟筛选方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101916330A (zh) * | 2010-08-06 | 2010-12-15 | 辽宁大学 | 一种以Keap1为靶点的新型防癌抗癌药物的虚拟筛选方法 |
| CN102222176A (zh) * | 2011-06-01 | 2011-10-19 | 山东大学 | 一种快速发现以Bcl-2蛋白为靶点的先导化合物的方法 |
| CN102270281A (zh) * | 2011-06-01 | 2011-12-07 | 山东大学 | 一种以cdk1为靶点的抗癌药物的筛选方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100839600B1 (ko) * | 2006-12-04 | 2008-06-19 | 한국전자통신연구원 | 특이 결합부위 자동추출을 이용한 리간드 검색 장치 및 그방법 |
-
2012
- 2012-06-14 CN CN201210195233.2A patent/CN102750459B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101916330A (zh) * | 2010-08-06 | 2010-12-15 | 辽宁大学 | 一种以Keap1为靶点的新型防癌抗癌药物的虚拟筛选方法 |
| CN102222176A (zh) * | 2011-06-01 | 2011-10-19 | 山东大学 | 一种快速发现以Bcl-2蛋白为靶点的先导化合物的方法 |
| CN102270281A (zh) * | 2011-06-01 | 2011-12-07 | 山东大学 | 一种以cdk1为靶点的抗癌药物的筛选方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102750459A (zh) | 2012-10-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zhavoronkov et al. | Potential 2019-nCoV 3C-like protease inhibitors designed using generative deep learning approaches | |
| Chowdhury et al. | A computational approach to design potential siRNA molecules as a prospective tool for silencing nucleocapsid phosphoprotein and surface glycoprotein gene of SARS-CoV-2 | |
| Foglierini et al. | HCMV envelope glycoprotein diversity demystified | |
| CN102750459B (zh) | 一种以人乳头瘤病毒hpv的dna解旋酶e1为靶点的抗hpv病毒药物的虚拟筛选方法 | |
| Deng et al. | The epidemic origin and molecular properties of B′: a founder strain of the HIV-1 transmission in Asia | |
| CN104834831B (zh) | 一种基于三维定量构效关系模型的一致性模型构建方法 | |
| CN101916330A (zh) | 一种以Keap1为靶点的新型防癌抗癌药物的虚拟筛选方法 | |
| Wan et al. | Spatial-and valence-matched neutralizing DNA nanostructure blocks wild-type SARS-CoV-2 and omicron variant infection | |
| Wu et al. | The inhibition of Mpro, the primary protease of COVID-19, by Poria cocos and its active compounds: a network pharmacology and molecular docking study | |
| Sattari et al. | Repositioning therapeutics for COVID-19: virtual screening of the potent synthetic and natural compounds as SARS-CoV-2 3CLpro inhibitors | |
| Gupta et al. | In silico accelerated identification of structurally conserved CD8+ and CD4+ T-cell epitopes in high-risk HPV types | |
| El Mazouri et al. | Genetic diversity and evolutionary dynamics of the Omicron variant of SARS-CoV-2 in Morocco | |
| Lin et al. | Precise genotyping and recombination detection of Enterovirus | |
| Zhou et al. | Molecular epidemiological characteristics of SARS-CoV-2 in imported cases from 2021 to 2022 in Zhejiang Province, China | |
| Alizon et al. | Towards a multi-level and a multi-disciplinary approach to DNA oncovirus virulence | |
| Frank | “Just in Time”: the role of cryo-electron microscopy in combating recent pandemics | |
| CN110379462A (zh) | 一种基于Illumina技术组装中华金腰叶绿体基因组序列的方法 | |
| Sahlan et al. | In Silico Identification of Propolis Compounds Potential as COVID-19 Drug Candidates Against SARS-CoV-2 Spike Protein | |
| Ejeh et al. | Structure-based design, drug-likeness, and pharmacokinetic studies of novel substituted pyrimidine derivatives as potent HCV NS3/A4 protease inhibitors | |
| CN102178668B (zh) | 亚油酸在制备抗hpv药物中的应用 | |
| Islam et al. | High-throughput virtual screening of marine algae metabolites as high-affinity inhibitors of ISKNV major capsid protein: An analysis of in-silico models and DFT calculation to find novel drug molecules for fighting infectious spleen and kidney necrosis virus (ISKNV) | |
| Joshi et al. | Structure-based screening of novel lichen compounds for SARS coronavirus main protease (Mpro) and angiotensin-converting enzyme 2 (ACE2) inhibitory potentials as multi-target inhibitors of COVID-19 | |
| Tallungan et al. | Analysis of Inhibitory Potential of Bioactive Compounds from Langusei (Ficus minhassae Tesym. & De Vr.) against SARS-CoV-2 using an in silico Approach | |
| CN105803109A (zh) | 一种利用等温核酸扩增技术快速检测黄热病毒的方法 | |
| 王迪 | Screening of differential microRNA and prediction of target genes between human immunodeficiency virus-Burkitt lymphoma and diffuse large B cell lymphoma |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20211227 Address after: 210001 zone k146, floor 5, Yindu Jinchuang Plaza, No. 2, Shuiximen street, Qinhuai District, Nanjing City, Jiangsu Province Patentee after: Nanjing Yitian Biotechnology Co.,Ltd. Address before: 110136 58 Shenbei New Area Road South, Shenyang, Liaoning. Patentee before: LIAONING University |
|
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160413 |