药物设计学第五章药物发现的虚拟筛选方法

1、第五章第五章 药物发现的虚拟筛选方法药物发现的虚拟筛选方法化学信息学化学信息学生物信息学生物信息学虚拟筛选虚拟筛选先导物先导物化学生物学化学生物学候选药物候选药物靶点大分子靶点大分子第一节第一节 概述概述虚拟筛选：针对重要疾病特定靶标生物大分子的三维结构或虚拟筛选：针对重要疾病特定靶标生物大分子的三维结构或定量构效关系定量构效关系(qsar)(qsar)模型，从现有小分子数据库中，搜寻与模型，从现有小分子数据库中，搜寻与靶标生物大分子结合或符合靶标生物大分子结合或符合qsar qsar 模型的化合物，进行实验模型的化合物，进行实验筛选研究。筛选研究。虚拟筛选流程虚拟筛选流程基于分子对接虚拟筛选

2、流程基于分子对接虚拟筛选流程sarssars爆发，从现有药物中快速筛选，老药，安全，直接进临爆发，从现有药物中快速筛选，老药，安全，直接进临床。床。sarssars冠状病毒冠状病毒(sars-cov)(sars-cov)蛋白水解酶蛋白水解酶(mpro)(mpro)的同源蛋白模拟的同源蛋白模拟结构；结构；虚拟筛选现有药物库虚拟筛选现有药物库mdl/cmc(http:/ mpro sars-cov mpro 的抑的抑制剂；制剂；(a) (a) 肉桂硫胺的化学结构；肉桂硫胺的化学结构；(b) (b) 肉桂硫胺与肉桂硫胺与sars-cov mprosars-cov mpro活性位点的作用方式活性位点的

3、作用方式（绿线：（绿线：-nh-nh 或或-ch-ch 氢键相互作用；红线：氢键相互作用；红线：-ch-ch o o 氢键相互作用；蓝线：氢键相互作用；蓝线：c145 c145 的巯基的巯基(-sh)(-sh)与肉桂硫胺的与肉桂硫胺的c=c c=c 双键作用双键作用酶水平抑制活性酶水平抑制活性ic50 ic50 为为5.0m5.0m，能抑制病毒颗粒，能抑制病毒颗粒的复制；的复制；细胞水平抑制细胞水平抑制sars-cov sars-cov 的的活性为活性为31m31m；是冠状病毒；是冠状病毒mpro mpro 的专一性抑制剂。的专一性抑制剂。第二节第二节 化学信息处理化学信息处理一、化学信息的表

4、示方法一、化学信息的表示方法1 1、一维结构表示、一维结构表示smilessmiles编码的立体化学信息表示编码的立体化学信息表示第二节第二节 化学信息处理化学信息处理一、化学信息的表示方法一、化学信息的表示方法2 2、二维结构表示、二维结构表示用图表示苯丙氨酸的结构用图表示苯丙氨酸的结构（1 1）图论基础和图的矩阵表示）图论基础和图的矩阵表示第二节第二节 化学信息处理化学信息处理一、化学信息的表示方法一、化学信息的表示方法2 2、二维结构表示、二维结构表示乙醛的邻接矩阵表示和矩阵简化步骤乙醛的邻接矩阵表示和矩阵简化步骤（2 2）化合物结构的矩阵表示）化合物结构的矩阵表示两个原子间有键相两个原

5、子间有键相连，矩阵元素为连，矩阵元素为1 1第二节第二节 化学信息处理化学信息处理一、化学信息的表示方法一、化学信息的表示方法3 3、三维结构表示三维结构表示甲烷立体结构的笛卡儿坐标系统甲烷立体结构的笛卡儿坐标系统原子坐标随坐标原点、原子坐标随坐标原点、分子的摆放位置而变分子的摆放位置而变1 1，2-2-二氯乙烷的内坐标系统二氯乙烷的内坐标系统描述分子中原子描述分子中原子和键的相对位置和键的相对位置第二节第二节 化学信息处理化学信息处理一、化学信息的表示方法一、化学信息的表示方法常见的分子存储格式常见的分子存储格式（1 1）molmol格式文件格式文件（2 2）sdfsdf格式文件格式文件第二

6、节第二节 化学信息处理化学信息处理二、化合物数据库的生成二、化合物数据库的生成isisisis信息管理系统信息管理系统化学信息管理系统（化学信息管理系统（isis/hostisis/host）化学信息生成和管理软件（化学信息生成和管理软件（isis/baseisis/base）化学绘图软件（化学绘图软件（isis/drawisis/draw）第二节第二节 化学信息处理化学信息处理二、化合物数据库的生成二、化合物数据库的生成1 1、化合物数据库的生成、化合物数据库的生成2 2、三维结构数据库、三维结构数据库19871987年出现年出现concordconcord程序，第一个能把二维连接表转换成三

7、维坐程序，第一个能把二维连接表转换成三维坐标标, ,摆脱了三维结构信息的来源限制摆脱了三维结构信息的来源限制, , 使任何拥有二维分子数据使任何拥有二维分子数据库的实验室和公司都可以建立自己的三维结构数据库。之后产库的实验室和公司都可以建立自己的三维结构数据库。之后产生生chem-xchem-x、cobracobra、corinacorina和和molgeomolgeo等一些三维结构转换程序。等一些三维结构转换程序。利用利用concordconcord程序得到的三维结构数据库有：程序得到的三维结构数据库有：fcd-3dfcd-3d、mddr-3dmddr-3d、cmc-3dcmc-3d、cas

8、rfcasrf、cast-3dcast-3d、pomona-92cpomona-92c目前，各大制药公司和分子设计公司建立的三维结构数据库有：目前，各大制药公司和分子设计公司建立的三维结构数据库有：acd-3dacd-3d、sybyl-3dbsybyl-3db、chem-3dbschem-3dbs、nci-3dnci-3d。剑桥结构数据库（剑桥结构数据库（csdcsd）：小分子晶体结构）：小分子晶体结构布鲁克海文国家实验室蛋白质数据库（布鲁克海文国家实验室蛋白质数据库（pdbpdb）：大分子晶体结构）：大分子晶体结构蛋白质结构的蛋白质结构的4 4个层次：个层次：一级结构：组成蛋白质的氨基酸序列

9、；一级结构：组成蛋白质的氨基酸序列；二级结构：骨架原子间的相互作用形成的局部结构，二级结构：骨架原子间的相互作用形成的局部结构，如如alphaalpha螺旋，螺旋，betabeta片层和片层和looploop区等；区等；三级结构：二级结构在更大范围内堆积形成的空间三级结构：二级结构在更大范围内堆积形成的空间结构；结构；四级结构：主要描述不同亚基之间的相互作用。四级结构：主要描述不同亚基之间的相互作用。 第三节第三节 生物信息处理生物信息处理一、序列分析一、序列分析1 1、单个序列分析、单个序列分析酪氨酸磷酸酶的蛋白质序列酪氨酸磷酸酶的蛋白质序列酪氨酸磷酸酶的二酪氨酸磷酸酶的二级结构，其中级结构

10、，其中h h 代代表螺旋，表螺旋，e e 代表折代表折叠，叠，b b表示表示桥，桥，g g表示表示310310螺旋，螺旋，i i表表示示螺旋，螺旋，t t表示氢表示氢键转角，键转角，s s代表转向，代表转向，图图7.17.1（c c）显示的）显示的是该。是该。 一、序列分析一、序列分析1 1、单个序列分析、单个序列分析2 2、双重序列比较、双重序列比较序列对比序列对比3 3、多重序列比较、多重序列比较序列对比序列对比鼠和小龙虾的胰蛋白酶鼠和小龙虾的胰蛋白酶mouseivggynceensvpyqvslns-gyhfcggslineqwvvsaghcyk-sriqvcrayfishivggtda

11、vlgefpyqlsfqetflgfsfhfcgasiynenyaitaghcvygddyenpsglqi*mouserlgehnievlegneqfinaakiirhpqydrktlnndimliklssravinarvstislptacrayfishvageldmsvnegseqtitvskiilhenfdydlldndisllklsgsltfnnnvapialpaq|-s-s-|mouseppatgtkclisgwgntassgadypdelqcldapvlsqakceasypg-kitsnmfcvgflecrayfishghtatgnvivtgwg-ttseggntpdvlqkvtvp

12、lvsdaecrddygadeifdsmicagvpe+*|-s-s-|mouseggkdscqgdsggpvvcng-qlqgvvswgdgcaqknkpgvytkvynyvkwikntiaancrayfishggkdscqgdsggplaasdtgstylagivswgygcarpgypgvytevsyhvdwikanav-数据库数据库pdb5&/swiss-prot#n;pir#owlembltremblgenbank/y3s-i0:?9(v%6s:nprosite3r,b/_7p7i-k/uswiss-model!f)&f,#sw8)fswiss-3dimagedsspfsspscop

13、cath2z5o*;a6e-y&_pfam说明说明1j1e#f.q8u5蛋白质三维结构*s!i%h7*s7l)r6k蛋白质序列数据库7a,.d6,h8q蛋白质序列数据库非冗余蛋白质序列核酸序列数据库!embl的翻译数据库核酸序列数据库蛋白质功能位点:o*d4a-d/tt;p6l65m,s从序列模建结构/s2k,&n,o1|三维结构图示蛋白质二级结构参数已知空间结构的蛋白质家族蛋白质分类数据库蛋白质分类数据库;e7nq#$h&r3?132r(b蛋白质家族和结构域y.t1s5bh!;d.i%a网址链接网址链接/pdb/spr

14、ot//6u;x/s8p2m4e;c*_7w)dhttp:/www.bioinf.man.ac.uk/dbbrowser/owl/57:c/q0e7lhttp:/www.embl-heidelberg.de//sprot/,l&e$a8g(p2i0c+whttp:/www.ncbi./genbank//prosite/9d41a8?8v;?/swissmod/swiss-model.html$3j+y#x%phttp:

15、//sw3d/http:/www.cmbi.kun.nl/gv/dssp/*6dd&stxhttp:/www.ebi.ac.uk/dali/fssp/fssp.htmlhttp:/scop.mrc-lmb.cam.ac.uk/scop/&r9s&e/v&y:lhttp:/www.biochem.ucl.ac.uk/bsm/cath//二、蛋白质折叠与三维结构预测二、蛋白质折叠与三维结构预测从头（从头（ab initioab initio）预测法）预测法穿针引线穿针引线(threading)(threading)法法同源模建（同源

16、模建（homology modelinghomology modeling）法）法蛋白质的空间结构信息蕴涵于氨基酸序列中。蛋白质的空间结构信息蕴涵于氨基酸序列中。蛋白质预测工作的理论基础：系统的稳定状态通蛋白质预测工作的理论基础：系统的稳定状态通常是能量最小的状态常是能量最小的状态酪氨酸磷酸酶的折叠结构酪氨酸磷酸酶的折叠结构1.1.寻找同源模板：分辨率，有否晶体寻找同源模板：分辨率，有否晶体x-x-射线射线衍射结构，衍射结构，blast scoreblast scoreg g蛋白偶联受体（蛋白偶联受体（gpcrgpcr）：牛视紫质（）：牛视紫质（pdbpdb）2.2.用用t-coffeet-c

17、offee计算计算, ,确定目标和模板之间一确定目标和模板之间一级序列氨基酸残基间的对应位置级序列氨基酸残基间的对应位置第四节第四节 虚拟筛选虚拟筛选三维结构数据库搜寻三维结构数据库搜寻 相应的三维结构搜索软件促进了三维结构搜索的发展相应的三维结构搜索软件促进了三维结构搜索的发展, ,使三维结构搜索方法真正走向实用。使三维结构搜索方法真正走向实用。输入药效团结构，在三维结构数据库中搜索含有该药效团输入药效团结构，在三维结构数据库中搜索含有该药效团的分子；的分子；搜索到的分子有可能具有药效团所对应的某种生物活性，搜索到的分子有可能具有药效团所对应的某种生物活性，它可以在实验筛选之前为分子的生物活

18、性作可能性预报；它可以在实验筛选之前为分子的生物活性作可能性预报；当搜索结果中有一些分子，与已知活性的分子结构不同，当搜索结果中有一些分子，与已知活性的分子结构不同，意味着可能发现先导物。意味着可能发现先导物。数据库的分子是已经存在的化合物，合成方法已知，甚至数据库的分子是已经存在的化合物，合成方法已知，甚至已成为商品，可直接购买较快地进入生物测试阶段，提高新已成为商品，可直接购买较快地进入生物测试阶段，提高新药开发的效率。但化合物结构全部来自药开发的效率。但化合物结构全部来自x x射线晶体衍射和中射线晶体衍射和中子衍射实验，三维结构搜索发展缓慢。子衍射实验，三维结构搜索发展缓慢。 基于配体的

19、数据库搜寻方法：基于配体的数据库搜寻方法： 根据药效基团根据药效基团（pharmacophorepharmacophore）模型进行三维模型进行三维结构数据库搜寻。结构数据库搜寻。hivproteaseinhibitorshh4.1-6.1 4.4-6.4 1.8-3.8 建立一系列活性分子的药效构象；建立一系列活性分子的药效构象；抽提出共有的药效基团；抽提出共有的药效基团；在现有的数据库中寻找符合药效基在现有的数据库中寻找符合药效基团模型的化合物。团模型的化合物。比较著名的软件有：比较著名的软件有：catalystcatalyst（应用普遍）（应用普遍）unityunity分子对接法（分子对

20、接法（dockingdocking）：）： 将小分子配体对接到受体的活性位点，搜寻其合理的取将小分子配体对接到受体的活性位点，搜寻其合理的取向和构象，使得配体与受体的形状和相互作用的匹配最佳。向和构象，使得配体与受体的形状和相互作用的匹配最佳。 在药物设计中，分子对接方法主要用来从化合物数据库在药物设计中，分子对接方法主要用来从化合物数据库中搜寻与受体生物大分子有较好亲和力的小分子，从而发现中搜寻与受体生物大分子有较好亲和力的小分子，从而发现全新的先导化合物。分子对接由于从整体上考虑配体与受体全新的先导化合物。分子对接由于从整体上考虑配体与受体的结合效果，所以能较好地避免其他方法中容易出现的局部的结合效果，所以能较好地避免其他方法中容易出现的局部作用较好，整体结合欠佳的情况。作用较好，整体结合欠佳的情况。目前具代表性的分子对接软件：目前具代表性的分子对接软件： dockdock、f1exx f1exx 和和 goldgo