分子对接

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,分子对接,Molecular docking,1,计算机辅助药物设计的方法学,基于配体的药物设计方法,基于受体的药物设计方法,基于结构的药物设计,2,基于受体的药物设计方法,通过受体的特征以及受体和药物分子之间的相互作用方式来进行药物设计的方法。,“有的放矢”,主要方法,分子对接法,从头设计,3,分子对接,1.,概念,2.,原理,3.,一般过程,4.,分类,5.,代表性软件,6.DOCK,软件,7.AUTODOCK,软件,4,1,分子对接的概念,受体和药物分子之间通过空间匹配和能量匹配而相互识别形成分子复合物，并预测复合物结构的操作过程。,整体上考虑配体与受体的结合效果，较好地避免局部作用较好、整体结合欠佳的情况。,5,2,分子对接的原理,理论基础：,“锁和钥匙模型”,“诱导契合模型”,重要原则：,互补性：决定识别过程的选择性,预组织性：决定识别过程的结合能力,6,分子对接的最初思想起源于,Fisher E,提出的“锁和钥匙模型”。即受体与配体的相互识别首要条件是空间结构的匹配,配体 受体 复合物,受体配体的锁和钥匙模型,7,Oh boy! What a perfect match,这类方法首先要建立大量化合物（例如几十至上百万个化合物）的三维结构数据库，然后将库中的分子逐一与靶标分子进行“对接”（,docking,），通过不断优化小分子化合物的位置（取向）以及分子内部柔性键的二面角（构象），寻找小分子化合物与靶标大分子作用的最佳构象，计算其相互作用及结合能。在库中所有分子均完成了对接计算之后，即可从中找出与靶标分子结合的最佳分子（前,50,名或前,100,名）,8,3,分子对接的基本原理,药物与受体的结合强度取决于结合的自由能变化,G,结合,=,H,结合,-,T,S,结合,= -,RT,ln,K,i,大部分的分子对接法忽略了全部的熵效应，而在焓效应也只考虑配体与受体的相互作用能，即：,E,interaction,= E,vdw,+ E,electrostatic,+ E,h-bond,9,3,分子对接的一般过程,找出空穴，定出表面,调整受体位点或药物的构象,计算对接时受体药物相互作用能量,进行分子动力学模拟,复合物的全局最优结合构象,10,11,12,13,14,15,4,分子对接的分类,刚性对接：研究体系的构象不发生变化。,半柔性对接：研究体系尤其是配体的构象允许在一定范围内变化。,柔性对接：研究体系的构象是可以自由变化的。,16,3,分子对接的基本方法,（一） 刚性的分子对接方法,这种方法是最初的分子对接的方法，在对接中，小分子和蛋白质两种都保持刚性。,（,1,）基于最大团搜索的方法,（,Clique-Search Based Approaches,）,对接两个刚性分子可以理解为分子在空间的匹配问题，这种匹配可以是一种形状上的互补或相互作用。如氢键受体与氢键给体的互补。搜索在三维空间中有效的条件下的最大匹配,17,受体的活性位点,配体,有效匹配的距离图集,受体配体的示意图，字母代表特征部分如氢键等，相应的有效匹配的图集如右，三个环性顶点组织的三角形为这个图集的一个最大团,(clique),18,Dock,对接程序中刚性对接的算法就是基于这种思想,Dock,利用球集来表示受体活性位点和配体的形状,19,一系列的球集填充在受体活性位点的表面，这些球集代表能被配体占据的体积。配体可以用球集表示或者用自己的原子表示，在,Dock,程序中，四个有效匹配的对应点被考虑，先考虑配体中第一个球集与活性位点的球集的匹配，第二个点则满足,d ,，其中,d,为第二个匹配点中配体和受体的球心与第一个点球心的距离，第三个点又必需满足与前两个球心的距离限制，以上过程一直进行到找不到更多匹配点为止。,20,（,2,）基于几何哈希技术“,geometric hashing”,的方法,第一部分中，几何哈希表从被对接的一个配体或一系列配体中构建 。哈希矩阵含有配体名字和能调整配体在空间方向的参考框架。,第二部分即识别阶段，蛋白质的特征用来识别哈希矩阵，每一次匹配表示蛋白质的特征与哈希矩阵中已定义好方位的配体相匹配，具有大量匹配信息的哈希矩阵代表着具有几个吻合特征的配体和方位,21,(3,）基于,pose clustering,的方法,这种方法与几何哈希的方法相类似，也是一种基于模式识别的方法。,在,LUDI,模型中，如图所示，对每一个作用基团，定义作用中心和作用表面。受体的作用表面近似地用离散的点表示，和对应的配体的中心目标点相匹配。,三个氢键受体的作用表面,Pose clustering,算法中的作用点,22,（二）柔性对接的方法,（,1,）构象的系综方法,Flexibase,用来储存小分子库中每个分子的一系列不同构象，用距离几何和能量最小化的方法产生构象，每个分子根据,rmsd,的差异选择,25,个系列构象。每个构象采用,FLOG,刚性对接的方法进行对接。,23,（,2,）片段的方法,片断的方法是处理小分子柔性的最通用的方法，配体分割成一些小的片断，这些片断可以认为是刚性构象或一个小的构象系综。一般，有两种方法来处理,:,第一种方法是把一个片段放入受体的作用位点，然后加上余下的片段，这种方法称为连续构建 “,incremental construction”.,第二种方法把所有或一部分片段独立地放入受体的作用位点，再重新连接至到构成一个完整的配体分子，这种策略称为“放置,&,加” “,place & join”,24,第一个连续构建的算法是,Kuntz,发展的,Dock,程序，首先，一个单独的锚碎片通过手工选择对接进受体的活性结合位点，并且考虑了氢键的作用。其次这个锚的优势位置主要包含有有大量匹配氢键对，高打分值，和低相似性。接着，一种回溯的算法（,backtracking algorithm,）用来搜索整个配体在结合位点的非重叠放置空间，在当前位置加上一个碎片后，优化的方法用来减少立体的张力和改善氢键的几何构型。最后的位置通过过滤，优化和基于力场的方法来打分评价。,FlexX,也是一个基于连续构建算法的对接程序,25,（,3,） 遗传算法和进化规划,遗传算法开始应用到分子对接技术，其特点为 ：,第一步，一个称为染色体的线性表示符能够描述构型的所有自由度，找到这个染色体描述符是算法中最困难的一步。第二步，确定一个一个类似如打分函数的目标函数。,著名的,GOLD,软件包括了这种算法,26,（,4,）基于分子模拟的方法,模拟退火的方法,，,Autodock,程序就采用了这种方法,分子动力学的方法,Monte Carlo,模拟,，一种统计力学的方法，这种算法中最重要的两部分是自由度的描述和能量的评价，合适的自由度描述可以避免较高能量的构象，用键角、扭曲角等内座标来描述配体的柔性比用笛卡儿空间的三维座标描述要强，同样，能量的评价也是最耗时，这一步时间必须足够的长。,27,分子对接的主要问题,如何找到最佳的结合位置,遗传算法,模拟退火,如何评价对接分子之间的结合强度,非键作用能,基于分子表面的溶剂化计算,半经验的自由能计算,28,5,代表性对接软件,名称 构象搜索方法 结合评价方法 速度,Flex X (Sybyl),片段生长法半经验自由能快,LigandFit(Cerius2),蒙地卡罗模拟 半经验自由能快,Glide (,薛定谔软件,),系统搜索 半经验自由能一般,Gold,遗传算法半经验自由能快,Affinity,（,InsightII),蒙地卡罗,/MM/MD,分子力场慢,AutoDock,遗传算法半经验自由能一般,Dock,片段生长法分子力场 快,ICM-Dock,随机全局优化半经验自由能快,Fred (openeye),系统搜索半经验自由能快,29,Flex X,对接软件,商业软件模块（,Sybyl),对接过程,确定核心结构片段,采用形态聚类算法，放置核心结构,采用树形搜索算法，其他部分片段依次生长连接在核心结构上,30,LigandFit,对接软件,商业软件模块（,Cerius2),对接过程,发现活性位点,进行配体分子的构象搜索,进行分子对接,计算评分,31,Dock,对接软件,学术用户免费软件,http:/docking.org/,对接过程,准备蛋白质和配体分子,计算,MS,表面积,进行活性位点特征描述,建立评分格点,对接计算,评分函数,基于,AMBER,力场的分子间能量评分（范德华静电）、接触评分,32,AutoDock,对接软件,免费软件,http:/autodock.scripps.edu/,只能进行一个分子与蛋白质的对接计算，不能进行数据库对接，没有平行化功能。,不需要预先知道活性位点。,33,DOCK,软件,自动地模拟配体分子在受体活性位点地作用情况，并把理论预测最佳地相互作用方式记录下来。,自动搜索配体的,三维结构数据库。,34,Dock,软件,步骤,1.,分子的准备工作,2.,活性位点的确定,3.,格点对接,4.,柔性对接,35,分子的准备工作,在,Chimera,软件中进行,加氢原子,加电荷,得到的文件：,1,）,rec_charged.mol2,2,）,rec_noH.pdb,3,）,lig_charged.mol2,36,活性位点的确定,填充受体分子表面,口袋或凹槽的球集,每个负像代表一个,可能的活性位点,聚类分析,37,负像个数的简化原则,只考虑半径最小的负像,仅保留夹角小于,90,度的负像（位于口袋浅表）,与负像相切的表面点必须间隔,4,个残基,去掉半径大于,5,埃的负像,38,程序命令,生产分子表面点,dms input_file -n -w 1.4,o output_file,例如：,dms rec_noH.pdb-n -w 1.4,o rec.ms,得到文件：,rec.ms,39,不含氢原子的受体文件,形成负像文件名,程序命令,产生负像文件,sphgen,默认的参数文件：,INSPH,默认输入文件为,rec.mc,默认得到负像文件为,rec.sph,默认得到一个计算过程描述文件：,OUTSPH,40,41,格点计算,依据受体表面的这些结合点与配体分子的距离匹配原则，将配体分子投映到受体分子表面,保留仅是与受体有个的特征值,命令,showbox box.in,输出文件：,rec_box.pdb,42,43,44,45,46,47,AutoDock,对接软件,准备蛋白质和配体,蛋白质：加电荷和溶剂化，存为,pdbqs,格式,配体：加电荷，定义搜索的二面角和根片段,格点对接：保留探针原子和受体之间的相互作用能,对接计算：遗传算法,评价函数：半经验的自由能计算,范德华相互作用,氢键相互作用,静电相互作用,溶剂化作用,48,8,分子对接计算的注意点,小分子问题,起始构象对对接结果有一定影响,对接时应以代谢物的结构进行,对分子进行加电荷和加氢处理,蛋白质问题如何选择合理的蛋白质活性位点,对接问题,搜索结合模式的正确性、对接的效率、评分的正确性,采用多个软件进行评价，减少结合模式搜索误差,定量指标，需要结合分子动力学进一步评价,49,评价：打分函数,每一个对接的算术都会采用平衡了时效和精确度的简单自由能预测方法，现在的打分函数主要包括三种：基于经验的回归参数的方法；基于分子力场的方法和基于知识的方法、基于知识的打分函数 。,50,（,1,）基于力场的方法,只考虑热焓对能量的贡献，不考虑熵的影响，一般情况下，采用标准力场的非键作用能如真空静电和范德华作用能用作打分函数，如,DOCK,程序中采用,AMBER,的能量函数：,51,（,2,） 基于经验的打分函数,基于经验的打分函数用多元回归的方法拟合各种物理参数对结合自由能的贡献，如,FlexX,程序中采用下列函数，所采用的方程包括，配体旋转键的个数、氢键、离子键，疏水和芳香环的,堆积作用，以及亲水作用。这种方法能快速直接地估算结合自由能，,52,（,3,）基于知识的打分函数,最初应用于蛋白质结构预测，打分函数用统计力学的方法得自蛋白质配体的复合物结构，结合自由能用函数为分子间距离的平均能的加和来计算。基于知识的打分函数是一种比较有前途的方法,53,虚拟筛选的具体流程,54,包括,4,个步骤：受体模型的建立；小分子库的产生；计算机筛选和命中化合物的后处理。,第一步，受体模型的建立：,蛋白质结构的准备是虚拟筛选的重要一步。虚拟筛选的蛋白靶标的结构可以从,PDB,库,(,http:/www.rcsb.org/pdb/index.html,),中直接下载使用,也可以通过和家族中同源蛋白的序列、结构信息比较，同源模建而得,1,）大分子结构获取,55,2,）接着是,结合位点,的描述，选择合适的配体结合口袋对分子对接至关重要,一种是直接从配体受体复合物结构中抽出；,选择口袋有两种方式：,如果没有复合物结构，则需要根据生物功能如结合、突变等实验信息来手动选择结合部位,56,第二步，建立小分子数据库,二维结构用结构转换程序如,CORINA,、,CONCORD,实现三维结构的转化。,建好的三维结构加氢加电荷后，便可以用于对接程序。,57,第三步，对接和打分，这一步是虚拟筛选的核心步骤。,对接操作就是把每个小分子放到受体蛋白的配体结合位点，优化配体构像和位置，使之与受体有最佳的结合作用，给最佳结合构象打分，对所有化合物根据打分排序，然后从化合物库中挑出打分最高的小分子。,58,最后一步是命中化合物的后处理,通过计算分子的类药性质,ADME/T (,吸收,absorption,、器官分布,distribution,、体内代谢,metabolism,、排泄,excretion,和毒性,toxicity),性质的估算，排除那些不具有类药性质的分子。,可以利用一些经验规则如“五规则” 等，快速排除那些不适合进一步药物开发的分子。,59,通过以上四步处理，大部分分子从化合物库中剔除，形成一个合理大小的化合物库，仅对这些适合成药的化合物或购买、或合成、或分离得到，然后再进行实际的生物测试。,60,对接方法尚需解决的问题：,分子的柔性,溶剂化效应,打分函数,61,分子对接的应用,靶,标,靶标分类,靶标结构,小分子库大小,所用方法,抑制剂活性,M,实验数据,AmpC,-lactamse,Hydrolase,X-ray,200k,NWU DOCK,26,X-ray,复合物,BCR,ABL,Kinase,X-ray,200k,DOCK,25,细胞的抑制活性实验,Anthrax EF,Adenylyl cyclase,X-ray,200k,NWU DOCK,20,酶动力学实验,IMPDH,Dehydrogrnase,X-ray,3500k,FlexX,30,酶动力学实验,Casein kinase II,Kinase,Homology,400k,DOCK,0.08,抑制活性、构效关系,K,通道,Ion channel,Homology,50k,DOCK,10,细胞的抑制活性实验,Thyroid homone receptor,Nuclear receptor,Homology,250k,ICM,0,75,抑制活性实验,CDK2,Kinase,X-ray,50k,LIDAEUS,2,X-ray,复合物,TGF,RK,Kinase,X-ray,200k,Catalyst,0,005,X-ray,复合物,cyclophilin,Immunophilin,X-ray,Unity/FlexX,6,细胞的抑制活性实验,tRNA guanine transglycoslase,X-ray,800k,Unity/FlexX,0.25,酶动力学实验,PfDHFR,Reductase,Homology,230k,Catalyst/DOCK,0.9,酶动力学实验,-Amylase,Hydrolase,X-ray,200k,Unity/FlexX,NMR,SPR,层析,62,（一） 配体对,CDK2,，,CDK4,激酶选择性的研究,细胞周期蛋白依赖性蛋白激酶（,cyclin-dependent kinases, CDKs,）是一类,Ser/Thr,蛋白激酶，直接参与调控细胞分裂周期，顾名思义，,CDK,只在周期蛋白,Cyclin,活化下才能工作,.,现在已经发现了,13,种,CDK,蛋白激酶和,25,种周期蛋白,Cyclin,，但它们具体的生物功能还不十分明了。其中,CDK2,、,CDK5,、,CDK6,已经解析出晶体结构，但,CDKs,的一级序列存在高的同源性，它们之间一级序列的高同源性暗示了三维结构的相似性,CDK1,、,CDK2,、,CDK4,、,CDK6,是基于,CDKs,结构化学治疗癌症的主要靶标。至今为止，文献报告的,CDKs,抑制剂有多种,63,尽管化合物结构各异，但所有,CDKs,的小分子抑制剂有一些共同的特点：,（,1,）一般小的分子量 （, autoMS protein.pdb ,78,（四） 用,sphgen,产生表征口袋形状的球集，得到,protein.sph,文件，该文件即为球集文件。用,showsphere,程序由,protein.sph,生成,sphere.pdb,，用,Sybyl,或,Weblab,等程序观看球集的分布。,79,（五）用,showbox,产生一表征网格范围的文件,box.pdb,80,（六）计算打分用的网格文件，首先建立,grid4,的输入文件,grid.in,，所需文件还有,protein.mol2,和,box.pdb,，然后用以下命名生成网格文件：, grid4 i grid.in,81,（七）建立,dock,的输入文件,dock.in,，所需文件还有,protein.sph,、,database.mol2,和网格文件，然后进行,dock,：, dock4 i dock.in,82,（八）,Dock,完毕后得到结果文件,result.info,和,result.mol2,，以及保留恢复文件,result.rst,。,（九）中途非正常退出，可以用如下命令继续：,dock4 i dock.in -r,83,84,85,86,87,88,89,90,