《药物设计学第四讲CADD理论与方法.ppt》由会员分享，可在线阅读，更多相关《药物设计学第四讲CADD理论与方法.ppt（81页珍藏版）》请在三一办公上搜索。
1、1,1,1,第三章 基于结构的药物设计,1,2,2,基于结构的药物设计（计算机辅助药物设计）,第一节 计算机辅助药物设计的基本理论和方法第二节 基于配体的药物设计方法（定量构效关系方法，药效团模型）第三节 基于受体的药物设计方法（分子对接，虚拟筛选，全新药物设计）,1,3,3,第一节 计算机辅助药物设计的基本理论和方法,一、基于药物-靶标相互作用的药物设计二、类药性预测方法与初步虚拟筛选三、分子三维结构表示,1,4,4,思考题,类药性预测方法有哪些？各自的主要内容是什么？一般的分子力场包含哪几项？如何通过理论计算方法得到一个小分子的最低能量结构？,1,5,5,基于先导化合物和靶分子结合模式的先
2、导化合物的结构优化,已有药物-靶结合作用的增强增加新的药物-靶结合作用复杂分子的结构简化柔性结构的构象限制,一、基于药物-靶标相互作用的药物设计,1,6,6,1、已有药物-靶结合作用的增强,通过研究先导化合物与靶分子的作用模式可以发现，先导化合物的重要官能团可能不是出于最佳的位置，或者与靶分子的相互作用较弱。可以通过各种结构修饰方法来增强先导化合物与靶分子的相互作用，提高先导化合物对靶分子的亲和力和选择性。,1,7,7,结合基团的替换,变换结合基团是一种简单的增强药物-靶标相互作用的方法，通过分子对接等方法可以预测最佳的结合基团。在变换结合基团的过程中，需要考虑合成的可行性，一般避免对原合成路
3、线有大的改动。,1,8,8,可作用于和肾上腺素受体,仅作用于肾上腺素受体,烷基替换,烷基替换最为常见：甲基、乙基、丙基、丁基、异丙基、异丁基、叔丁基等的变换来考察链的长度和体积对结合的影响。对于有些能同时作用于多个同类受体的药物，可以通过烷基的改变，提高药物选择性。,1,9,9,抗真菌药,三唑环用氟代嘧啶环替代后，大大增强了与唑类抗真菌药物作用靶点羊毛甾醇14-去甲基化酶的相互结合作用。,芳基替换,芳香环和杂环的替换：通常可以改变环的大小和杂原子来增强药物与靶分子的结合相互作用。,1,10,10,结合基团位置的变化,在研究先导化合物与靶分子的作用模式时发现，有时结合基团并不处于一个最佳的位置，
4、通过调整结合基团的位置，可以增强先导化合物与靶分子的结合相互作用。,1,11,11,血管紧张素(ACE)转化酶,1,12,12,结合基团电子性质的变化,结合基团的电子性质在药物与靶分子的静电相互作用和氢键相互作用中起到非常重要的作用。,如：胺基的可质子化程度决定了与结合位点羧基负离子之间作用的强弱。芳香胺弱于脂肪胺。,1,13,13,对氢键的影响。如当羰基作为氢键受体时，其氧原子上的电子云密度越高，形成氢键的能力越强。,1,14,14,2、增加新的药物-靶结合作用,在基团选取过程中，结合分子模拟方法，通过理论计算预测新分子与靶分子的相互作用能，可大大提高设计的成功率。,在研究药物与靶分子的作用
5、模式时，经常可以发现已有的分子并不“完美”，并不能与全部的活性结合位点形成相互作用。,1,15,15,与阿片受体形成新的疏水相互作用。,与血管紧张素转化酶形成新的疏水相互作用。,1,16,16,3、复杂分子的结构简化,通过研究天然产物中得到的先导化合物与靶分子的相互作用模式，找出它们结构中对产生活性必须的部分，然后将不必要的部分删除，简化结构，提高先导化合物合成的可行性，而且有可能减少因多余基团带来的毒副作用。,神经化学信使胆囊收缩素拮抗剂,1,17,注 意,1、随着分子简化，分子体积变小，有可能与其他靶分子作用，因此每一步简化都要跟踪其活性。2、分子简化必须适度。过于简化的分子往往会导致活性
6、和选择性的降低，以及毒副作用的增多。,17,1,18,18,4、柔性结构的构象限制,有些先导化合物含有多个可旋转键（柔性键），分子构象数目多，则分子难于以正确的构象与靶分子结合，降低了结合亲和力。同时，同一个分子的不同构象可能与不同的靶分子结合，缺乏选择性，带来副作用。限制构象的方法：成环固定 引入双键，三键、酰胺、芳环等刚性官能团 立体位阻,1,19,19,引入双键,选择性地作用于多巴胺D3，避免与D2的副作用,具有更高的抗肿瘤活性,安定药,在分子的构象限制过程中，往往会使分子结构更为复杂，因此在构象限制过程中必须同时考虑合成的可行性。,1,20,二、类药性预测方法与初步虚拟筛选,类药性指化
7、合物与已知药物的相似性。具有类药性的化合物并不是药物，但是具有成为药物的可能，这一类化合物称为类药性分子或药物类似物分子。类药性是一个模糊的概念，在药物研发中，可以说类药性分子是高质量的先导化合物。先导化合物在化学上容易合成，并且具有ADME的性质。在药物研发过程中，测定活性是确认化合物作为先导物的直接指标，而许多特性在先导化合物的优化过程中无法测量。为了降低研发的风险，类药性被更多用于评估化合物可能存在的导致失败的特性。,20,1,21,多种类药性预测方法已经用于化合物合成前的类药性评价、类药性化合物的广泛收集、类药组合化学库的设计、高通量筛选前化合物的类药性预测等方面。对提高筛选化合物的质
8、量，降低人力物力财力消耗，提高先导化合物发现的几率、促进药物发展的进程具有重要意义。,21,1,22,（1）结构特征分析法 主要用于研究分子结构中是否具有药物分子常用的结构片段，虽然不能直接用来区分药物与非药物，但有助于定义药物分子的结构，对用于全新药物设计的结构片段库的建立和组合化学库的设计具有重要意义。,22,1、类药性预测方法,5100个药物,CMC（comprehensive medicinal chemistry database）,1179中不同骨架,32个片段可描述半数药物,42个片段可描述1/4药物,最常出现的取代基：非芳香杂环出现的频率是芳香杂环的两倍。最常出现的功能团是脂肪
9、叔胺、醇羟基、酰胺基。,1,23,（2）Linpinski五倍率法,23,1,24,五规则对一些化合物的预测,24,1,25,一些药物化学家还仿照里宾斯基五规则提出了用于筛选类先导化合物的五规则：,化合物的分子量在100至380道尔顿之间化合物中有尽可能少的氢键受体化合物中有尽可能少的芳香环化合物脂水分配系数的对数值在1到3之间化合物能够与靶受体相结合,25,1,26,（3）Ghose法和Oprea法,Ghose分析了CMC数据库中的6304个化合物，发现80%符合以下规则：-0.04AlgP 5.6（平均2.52）160相对分子质量480（平均357）40摩尔折射率130（平均97）20原子
10、总数70（平均48）Oprea研究了许多类药库和非类药库，发现三个性质对于区分类药分子和非类药分子非常重要：,26,1,27,（4）类药指数法,类药指数（drug-like index,DLI）是通过分析CMC数据库中4836个化合物的25个结构参数，构建了“类药”簇中心。计算化合物的DLI，需计算这25个参数值，然后与类药簇中心的分布进行比较，若化合物的每个参数值都与类药簇中心非常接近，则化合物的DLI就会很高。,27,1,28,2、初步虚拟筛选方法,一般可根据分子性质的计算值和分子中可能存在的反应活性子结构和毒性子结构来区分类药和非类药化合物。（1）五倍率法（2）排出含有不利官能团的分子,
11、28,1,29,29,1,30,常见的毒性骨架结构和毒性部位或类型,30,畸形,肿瘤,肺,血液,突变,血管,1,31,31,心脏,神经失控,皮肤,腺体,泌尿系统,肠胃,脑和中枢神经,肝,1,32,（3）与已知药物的骨架和侧链比较,32,11.80%,2.50%,3.80%,2.30%,2.30%,2.80%,2.30%,4.80%,CMC数据库中药物最常见骨架及出现概率,1,33,33,三、分子三维结构表示,1,34,34,1、分子结构表示；2、分子三维结构信息的获得；3、能量最小化和构象搜索。,1,35,35,1、分子结构的表示,C8H9NO3,1,36,36,线型模型Sticks mode
12、lWith Hydrogen Atoms,线型模型Sticks modelWithout Hydrogen Atoms,1,37,37,空间填充模型Space Filling(CPK),棒型模型 Tubes model,球棍模型 Balls and Sticks,线型点状模型Sticks and dots model,1,38,38,Surfaces(http:/sci.org/Science/Compchem/feature14.html),VdW surfaceConnolly molecular surfaceSolvent-accessible surface(Lee and Rica
13、hrds,Connolly)Excluded volume,VdW Surface,Connolly surface,Probe,Solvent-accessible Surface,1,39,39,Examples,VdW Surface,Solvent Accessible Surface,1,40,40,Tryptophan:Molecular Surface,1,41,41,Hierarchy Model of Protein Structure,The secondary structure is the basic structure unit,1,42,42,蛋白结构的表示,碳图
14、模型C plot model,缎带模型Ribbon model,卡通模型Cartoon model,1,43,43,2、分子三维结构信息的获得,X-ray CrystallographyNMR SpectroscopyComputer-generated 3D structures,构象分析和几何优化,分子力学量子力学,1,44,44,通过对衍射的位置、强度计算，读出原子坐标值，解析结构得到晶体空间结构,计算机分子模型技术可将上述数学数值和符号转化为高分辨率的分子。,（1）,1,45,45,化合物测定须得到单晶蛋白质结晶蛋白质-药物复合物结晶X射线晶体学测得的是晶体状态下的优势构象,1,46,
15、46,可测定溶液中分子结构优势构象，更好代表生物环境下的分子，还能研究动力学特征。可得到蛋白质或核酸与药物相互作用的构象。,（2）,1,47,47,量子力学：可精确计算分子和电子的性质和能量分子力学：可计算出分子的静态性质，如结构和分子能量（能量中不涉及电子能量项）分子动力学：可计算分子的动态性质，如分子的运动轨迹，也能计算出分子体系的热力学性质，如焓、熵和自由能。,（3）Computer-generated 3D structures,1,48,48,量子化学（quantum chemistry）,量子化学是一门以量子力学的基本原理和方法来研究化学问题的学科。它从微观角度对分子的电子结构、成
16、键特征和规律、各种光谱和波谱以及分子间相互作用进行研究，并藉此阐明物质的反应性以及结构和性能关系等等。量子化学学科的建立标志着化学摆脱了单一依靠实验的局面，将化学带入了一个新的时代。,量子力学研究微观粒子（电子、原子、分子）运动规律的理论。它用波函数描写粒子的运动状态，以薜定谔方程确定波函数的变化规律并对各物理量进行计算。,1,49,49,从头计算方法（ab initio）,以基本物理常数（Planck常量、电子质量和电量）以及元素的原子序数，不借助于任何经验参数，求解薜定谔方程计算结果精度高，可靠性大，但是计算量极大，消耗计算机时太多从头计算法的软件有Gaussian，SPARTAN等,1,
17、50,50,量子化学半经验计算法Semiempirical methods,为了回避从头计算方法的复杂运算，发展了若干种近似性更大的分子轨道理论计算方法。在物理模型上，它们都引入可调参数，体系基于Hartree-Fock-Roothaan方程，借用经验或半经验参数代替分子积分，统称为半经验(semi-empirical method)分子轨道法。半经验计算法可以用于计算较大的分子体系。,1,51,51,半经验计算法主要包括：1、完全不考虑双电子作用的单电子近似方法，如HMO、EHMO法等；2、通常所说的半经验量子化学计算方法，主要指的是建立在零微分重叠（Zero Differential Ov
18、erlap,ZDO）近似基础上的计算方法，该方法忽略积分值很小、且运算复杂的双电子积分。,如CNDO、INDO、NDDO、改进的MINDO、MNDO以及AM1和PM3等,1,52,52,分子力学又称分子力场，基于经典牛顿力学方程的一种计算分子的平衡结构和能量的方法，来研究分子体系的结构和性质：将原子看成大小不同的橡皮球，把化学键看成长度不等的弹簧，原子之间的化学键，具有标准的键长、键角等，分子内还存在非键相互作用。分子调整自身构象以尽量给出原子位置的最佳排布。由于计算量小，分子力学可研究包括成千上万个原子的分子体系，包括有机小分子、生物大分子。,分子力学方法,1,53,53,1,54,54,S
19、teric Energy,The difference between the actual energy of the system and the energy of a hypothetical system is that all structural parameters are at their equilibrium(minimum energy)values.,Steric Energy=Estretch+Ebend+Etorsion+Estretch-bend+Etorsion-stretch+EVdW+Ehydrogen+Eelectrostatic+Edipole,Est
20、retchStretch energy(over all bonds)EbendBending energy(over all angles)EtorsionTorsional(dihedral)energy(over all dihedral angles)Estretch-bend Stretch bend energyEtorsion-stretchTorsion stretch energyEVdWVan Der Waals energy(over all atom pairs 1,3)EhydrogenHydrogen bonding energyEelectrostaticElec
21、trostatic energy(over all charged atom pairs 1,3)EdipoleDipole energy,1,55,1,56,常用的分子力场,56,1,57,57,Angle Bending,AMBER:,MM2:,MM3:,1,58,Uses of“Steric Energy”,“Steric energy”has NO physical meaning,and it is defined differently in different programsTherefore it CAN NOT be used to compare structures c
22、alculated by different programsIts use is limited to comparing ISOMERIC structures having the SAME number and kinds of bonds(conformers,stereoisomers).,1,59,Successes of MM Calculations,Calculations are very fastGeometry optimizations of small to medium-size molecules can be accomplished on a pcConf
23、ormations of macromolecules(including biomacromolecules such as peptides and polysaccharides)can be calculated using workstations or parallel processing computers.,1,60,60,A General Force Field Calculation,InputAtom TypesConnectivity Starting GeometryCalculate molecular properties at initial geometr
24、yGeometry optimizationCalculate molecular properties at final geometryOutputMolecular structureMolecular energyDipole momentsetc.etc.etc.,1,61,61,分子力学计算 MM2 计算方法 半经验公式 MOPAC 量子力学计算 从头计算 Gaussian,计算方法,1,62,62,量子力学计算法计算结果的精度要比分子力学计算法高的多，但每一钟计算方法都有它独特的用法。.分子力学计算方法MM2 应用范围：体系中有几千个原子；有机物、低聚核苷酸、羧氨酸以及糖等大型分
25、子；气象（只能用分子力学计算法）。包含的可用技术：能量最低化；构象搜索；同一分子的个别点能量比较；分子动力学研究。,1,63,63,.量子计算方法（1）半经验方法-MOPAC 含有120个重原子并且原子总数在300以内；有机物、有机金属化合物、低聚物（包括缩氨酸、核苷、糖类）；气相或绝对溶剂环境。（2）从头计算法 体系的原子总数在30个以内；气相或绝对溶剂环境；基态、过渡态和激发态的研究（主要方法）；静电势；偶极矩；过渡态几何形状和能量；键分裂能。（3）从头计算法Gaussian程序 Chem3D中添加了Gaussian的菜单，但需要将程序安装后才能使用。,1,64,64,3、Energy M
26、inimization（能量优化）,Starting geometry,Easy to implement;useful for well defined structuresDepends strongly on starting geometry,+,1,65,65,分子力学的能量优化方法较通用，计算速度快，但在最小化中只允许体系能量下降，而不能上升，因而有着陷入局部最小值的缺陷。,1,66,66,1,67,67,Population of Minima,Global minimum,Most populated minimum,Active Structure,1,68,68,基于导数能
27、量最小化方法总结,Steepest Descent(最速下降法)Robust.Convergence guaranteed on quadratic surfaces.Slow to converge especially near the minimum.Conjugate gradient,Quasi NR(共轭梯度法)Converges in approximately N steps for N degrees of freedom.Convergence guaranteed on quadratic surfaces.Unstable away from the minimum(i
28、.e.,when surface not quadratic).Newton-Raphson(牛顿-拉普森法)Converges in one step on quadratic surfaces.Unstable away from the minimum.Computationally expensive.,1,69,69,Population of Minima,No minimization method can guarantee the location of the global energy minimum.,But,is this what we really want?,1
29、,70,70,Conformational Search Outline,1,71,71,构象搜索法,系统搜索法（systematic searching）非系统搜索法（non-systematic searching）蒙特卡洛方法（Monte Carlo method）模拟退火方法（Simulated Annealing Algorithm）遗传算法（Genetic Algorithm）,1,72,系统搜索法（systematic searching），系统地搜索整个构象空间，寻找能量最低点。该法在搜索条件的参数下不会遗漏一个构象。最基本的系统搜索法是网格搜索（grid search），也称
30、为树状搜索（tree search）。,72,1,73,非系统搜索法,73,蒙特卡洛方法（Monte Carlo method），是一种统计抽样方法，它的命名来自于摩纳哥的赌城。基本思想是在求解的空间中随机采样并计算目标函数，以在足够多的采样点中找到一个叫高质量的最优解为最终目标。此方法取样恰当的话，对低能量构象取样概率大，缺点是计算量极大。,1,74,非系统搜索法,74,模拟退火方法（Simulated Annealing Algorithm），在分子动力学计算时，以模拟固体退火的过程来求解组和优化，以获得分子的全局优势构象。计算时先使体系升温，在高温（通常500-600K）下进行分子动力学
31、模拟，使分子体系有足够的能量，克服柔性分子中存在的各种能垒，越过能垒后从高温慢慢降低温度，搜索全部构象空间，从中选出一些能量相对较小的构象，然后逐渐降温，再进行分子动力学模拟，此时较高的能垒已无法越过，在极小化后去除能量较高的构象，最后可以得到相应的能量最小的优势构象。,1,75,非系统搜索法,75,遗传算法（Genetic Algorithm），将达尔文进化论的自然选择和优化过程引申到科学计算的优化过程中，把待解决的问题化解为一个寻优问题，利用编码方式把寻优问题表达成一个二进制串，经选择、交配和变异三种操作，不断产生新群体，解决组和优化问题。,1,76,76,理论计算获得一个分子最低能量构象
32、的步骤,2D分子结构的产生；转换2D分子结构为3D空间结构；构象搜索；能量最小化。,1,77,77,ChemDrawISISDraw,2D分子结构的产生,1,78,78,绘制分子结构,若知道分子的系统命名，则可迅速得到其分子结构，但本身收录的分子是有限的。操作步骤：点击Structure菜单里的Convert name to structure，弹出对话框，输入分子的名称，就可得到分子结构。,1,79,79,绘制结构式,1.普通分子式的绘制绘制右边的三肽分子结构打开文档在工具栏中选择苯环工具在串口单击，绘制苯，用五元环工具绘制五元环,1,80,80,选择单双键工具分别绘制单键、双键用文本输入工具输入元素符号选定全部分子，在Structure菜单中选择Clean Up Structure，调整分子结构，这样得到完美的三肽分子结构。最后，给分子结构命名。用菜单中的“Convert Structure to name”命令。,1,81,81,Chem3DHyperchemSybylDSSchrodinger,转换2D分子结构为3D空间结构分子结构的能量最小化,
