1、第四章 基于结构的药物设计Structure-based drug designSBDD.第一节 概述1.1 1.1 受体理论和药物受体相互作用受体理论和药物受体相互作用1.2 1.2 药物结构与生物活性的关系药物结构与生物活性的关系1.3 1.3 基于结构的药物设计基于结构的药物设计.1.1 1.1 受体理论和药物受体理论和药物-受体相互作用受体相互作用受体:受体:能识别和结合生物活性物质,并产生 生物效应的结构。生物活性物质生物活性物质:内源性活性调节物:维持机体机能的基本生理机制。外源性药物:干预机体生理生化作用。.外源性药物外源性药物受体激动剂(受体激动剂(agonistagonist
2、):药物与内源性物质一样,产生相似的生物作用。受体拮抗剂(受体拮抗剂(antagonistantagonist):药物与受体结合后阻碍了内源性物质与受体结合,而导致生物作用的抑制。.1.1.1 1.1.1 药物药物-受体相互作用受体相互作用空间形状契合空间形状契合 形状互补形状互补诱导契合诱导契合*通过各自的变构作用通过各自的变构作用形状互补形状互补相互作用力契合相互作用力契合 性质互补性质互补.A.A.空间形状契合空间形状契合:受体生物大分子的结构具有不对称性,药物受体生物大分子的结构具有不对称性,药物和受体的形状互补才能结合起效。和受体的形状互补才能结合起效。肾上腺素能激动肾上腺素能激动-
3、和和-肾上腺素受体肾上腺素受体而异丙肾上腺素仅能激动而异丙肾上腺素仅能激动-肾上腺素受体肾上腺素受体肾上腺素肾上腺素异丙肾上腺素异丙肾上腺素.立体因素的影响立体因素的影响 药物多种对应异构体中只有一个能与受体发生药物多种对应异构体中只有一个能与受体发生特异性结合。特异性结合。-OH的立体结构对活性影响显著,的立体结构对活性影响显著,-碳碳为为R构型的肾上构型的肾上腺素左旋体是腺素左旋体是-碳碳为为S构型的右旋体的构型的右旋体的12倍。倍。.一个对映体完全占据受体结合部位一个对映体完全占据受体结合部位 另一个对映体仅能部分匹配另一个对映体仅能部分匹配.B.诱导契合诱导契合(induced-fit
4、)受体与配体相互作用时,受体结合部位受体与配体相互作用时,受体结合部位受药物的诱导而发生构象的变化;柔性的药受药物的诱导而发生构象的变化;柔性的药物分子也会发生构象变化,产生可逆的、互物分子也会发生构象变化,产生可逆的、互补性的契合,即为诱导契合。补性的契合,即为诱导契合。.C.C.相互作用力契合相互作用力契合化学反应(即成键)、分子间作用力(即缔合)化学反应(即成键)、分子间作用力(即缔合)药物药物-受体相互作用力的类型和性质:受体相互作用力的类型和性质:共价键共价键离子键离子键氢键氢键疏水键疏水键离子离子偶极键偶极键范德华力范德华力偶极偶极偶极偶极诱导偶极诱导偶极.有机磷杀虫剂可以将乙酰胆
5、碱酶中的丝氨酸残基的羟基酯化,使其失去水解乙酰胆碱的活性。共价键共价键(covalent bond):两个或多个原子共同使用它们的外层电子,达到电子饱和的状态,由此组成比较稳定的化学结构。有机磷类脱磷酰基水解作用其中PX:代表有机磷酸酯杀虫剂;X:代表脱离基;PE:磷酰化酶,反应的三步为:(1)形成可逆性复合体(在酶的阴离子部位)。(2)磷酰化反应。有机磷是通过酰化反应产生对AchE的抑制作用。(3)酶的复活。.范德华力范德华力(van der Waals force,VDW)分为三种作用力:色散力、诱导力和取向力色散力、诱导力和取向力色散力色散力:所有分子或原子间都存在。分子由于电子的运由于
6、电子的运动产生的动产生的的瞬时偶极间的作用力。分子的变形性越大(一般分子量愈大分子量愈大,变形性愈大),色散力越大。分子的电离势越低(分子内所含的电子数愈多电子数愈多),色散力越大。.取向力取向力(orientation force)取向力发生在极性分子极性分子之间。当极性分子互相接近时,分子的永久偶极之间同极相斥、异极相吸,使分子在空间按一定取向排列吸引,而处于较稳定的状态。这种永久偶极间的吸引力称为取向力。取向力与分子的偶极矩平方成正比,即分子的极性越分子的极性越大大,取向力越大取向力越大。.诱导力诱导力(induction force):在极性分子非极性分子之间,极性分子极性分子之间,都
7、存在诱导力。极性分子与非极性分子接近时,极性分子的永久偶极产生的电场使非极性分子极化产生诱导偶极。永久偶极与诱导偶极间的吸引力称为诱导力。诱导力同样存在于极性分子之间,对极性分子来说,诱导力是附加的取向力。.疏水键(疏水键(hydrophobic bond):):疏水键的本质跟两个基团的疏水性有关,并非两个集团之间的共价或非共价作用,就是两个集团的某区域都对水有排斥作用所以将两个基团连接在一起。.氢键氢键(hydrogen bond):氢原子与负电性杂原子共价结合,与另一具有未共用电子对的杂原子形成一种弱的静电引力。.氢键氢键.氢键氢键.离子键(离子键(ionic bond,ion-ion b
8、ond)两个具有相反电荷的离子间的静电引力。两个具有相反电荷的离子间的静电引力。.偶极偶极-偶极作用(偶极作用(dipole-dipole interaction)在药物和受体分子中,在药物和受体分子中,当碳原子和其他电负当碳原子和其他电负性较大的原子,如性较大的原子,如N、O、S、卤素等成键时,、卤素等成键时,由于电负性较大原子由于电负性较大原子的诱导作用使得电荷的诱导作用使得电荷分布不均匀,导致电分布不均匀,导致电子的不对称分布,产子的不对称分布,产生生电偶极电偶极。极性分子本身存在的极性分子本身存在的偶极,是偶极,是固有偶极固有偶极。.诱导偶极作用(诱导偶极作用(induced dipo
9、le interaction)极性分子极性分子对非极性分子或极性分子诱导产生诱导产生的偶极。.离子离子-诱导偶极作用(诱导偶极作用(ion-induced dipole interaction)药物分子药物分子的偶极偶极受到来自生物大分子的离子离子吸引,而产生相互作用,这种相互作用对稳定药物受体复合物起到重要作用,但是这种离子-偶极的作用比离子产生的静电作用要弱得多。离子-偶极、偶极-偶极相互作用的例子通常见于羰基羰基类化合物。.1.2 药物的结构与生物活性的关系药物的结构与生物活性的关系(SAR)结构性质分子体积分子体积电荷分布电荷分布立体结构立体结构.药药 物物药效基团药效基团药动基团药动
10、基团毒性基团毒性基团.1.2.1 药效基团药效基团(pharmacophore或或biophore)指一系列生物活性分子所指一系列生物活性分子所共有共有的、对药物活性起的、对药物活性起决定作用决定作用的的分子片段及其三维空间位置分子片段及其三维空间位置结构特征。结构特征。能被受体所能被受体所识别识别、与受体受点、与受体受点结合结合起关键作用。起关键作用。当与受体受点结合后,产生特定的当与受体受点结合后,产生特定的生理活性。生理活性。具有类似结构的化合物往往有着具有类似结构的化合物往往有着相近的药理作用。相近的药理作用。.沙可来新沙可来新 嘧啶苯芥嘧啶苯芥双双-氯乙氨基氯乙氨基.阿片类止痛药的药
11、效基团阿片类止痛药的药效基团.1.2.2 毒性基团毒性基团(toxicophore)在一些药物中,如果药效基团所产生的生物效应在一些药物中,如果药效基团所产生的生物效应为为毒性反应毒性反应,该基团为毒性基团。,该基团为毒性基团。毒性基团往往存在于对病原体(微生物或癌细胞)毒性基团往往存在于对病原体(微生物或癌细胞)的化学治疗药物中,毒性的的化学治疗药物中,毒性的选择性越好,则药物选择性越好,则药物越安全。越安全。其他类药物应其他类药物应避免避免有毒性基团或潜在的毒性基团有毒性基团或潜在的毒性基团(能在体内经转化生成毒性基团)的存在。(能在体内经转化生成毒性基团)的存在。.1.2.3 药动基团(
12、kinetophore)药物中参与体内药物的吸收、分布、代谢和排泄药物中参与体内药物的吸收、分布、代谢和排泄(ADME)过程的基团,它本身不具有显著的生物活性,)过程的基团,它本身不具有显著的生物活性,决定着药物的药动学性质。决定着药物的药动学性质。药动基团通常是模拟自然界存在的物质,比如氨基酸、磷药动基团通常是模拟自然界存在的物质,比如氨基酸、磷酸基、糖基等生物代谢基本物质,经化学键与药效基团结酸基、糖基等生物代谢基本物质,经化学键与药效基团结合,使药物分子具易被合,使药物分子具易被转运转运的性质。它可改变药物在体内的性质。它可改变药物在体内的转运,或使作用的转运,或使作用定位定位化。化。沙
13、可来新沙可来新 嘧啶苯芥嘧啶苯芥双双-氯乙氨基氯乙氨基.1.2.4 药物构象药物构象优势构象:优势构象:在有机化合物分子中,碳-碳单键的自由旋转,引起结合在碳原子上的原子或基团的相对位置发生改变,产生若干种不同的空间排布方式,称为该分子的构象。其中,能量最低的构象最稳定,称为优势构象。活性构象活性构象:药物为了与受体的结合部位结合而发生的构象重组,形成能量较高的特定构象模式,即活性构象活性构象。药效构象药效构象:药物分子与受体结合时所采取的实际构象实际构象,并不一定采取它的优势构象,实际构象为药效构象。药效构象。.药效构象药效构象(pharmacophoric conformation).1.
14、3 基于结构的药物设计基于结构的药物设计间接药物设计间接药物设计 直接药物设计直接药物设计 .第二节第二节 计算机辅助药物设计计算机辅助药物设计Computer-Aided Drug Design.2.1 计算机辅助药物设计计算机辅助药物设计 利用计算机的计算、逻辑判断、图形显示等功能进行合理药物设计。特点:特点:由经典的定量构效关系发展为由经典的定量构效关系发展为三维定量构效关系。三维定量构效关系。设计中浩繁的计算、数据的存储和处理、显示、预测设计中浩繁的计算、数据的存储和处理、显示、预测等,均由计算机来完成。等,均由计算机来完成。将药物将药物-受体作用受体作用可视化。可视化。.2.2 2.
15、2 计算机辅助药物设计所依赖的理论和技术计算机辅助药物设计所依赖的理论和技术计算机辅助计算机辅助病因和靶点的确认病因和靶点的确认 化合物分析制备化合物分析制备 分子模建分子模建 .2.3 2.3 计算机辅助药物设计的意义计算机辅助药物设计的意义指导有目的地开发新药,减少盲目性和偶然性。加快研制新药速度,节省人力、物力和财力。为研究者提供理论思维形象化的表达,直观设计,理解和解释实验结果。局限:只是辅助性工具,仍需研究者的经验判断和指导。.2.4.1 三维结构的理论计算方法理论计算方法理论计算方法量子化学量子化学分子力学分子力学分子动力学分子动力学2.4 CADD2.4 CADD的基本理论和技术
16、的基本理论和技术.量子化学量子化学(quantum chemistry)1913年年Bohr提出:原子中的电子提出:原子中的电子只能处于包含基态在内的定态上只能处于包含基态在内的定态上,电子在两个定态之间跃迁而改,电子在两个定态之间跃迁而改变它的能量,同时辐射出一定波变它的能量,同时辐射出一定波长的光,光的波长取决于定态之长的光,光的波长取决于定态之间的能量差间的能量差量子论量子论量子力学研究微观粒子(电子、原子、分研究微观粒子(电子、原子、分子)运动规律的理论。它用波函数描写粒子的运动状态,以薜定谔方程确定波函数的变化规律并对各物理量进行计算。.量子化学量子化学以量子力学的基本原理和方法来研
17、究化学以量子力学的基本原理和方法来研究化学问题的学科。它从微观角度对分子的问题的学科。它从微观角度对分子的电子结构、电子结构、成键特征和规律、各种光谱和波谱以及分子间相成键特征和规律、各种光谱和波谱以及分子间相互作用互作用进行研究并藉此阐明物质的特性以及结构进行研究并藉此阐明物质的特性以及结构和性能的关系等。和性能的关系等。计算方法:计算方法:从头计算法从头计算法 半经验计算法半经验计算法量子化学可计算出分子的各种参数:量子化学可计算出分子的各种参数:分子结构、分子结构、电子结构、系统总能量、各个轨道的分子信息。电子结构、系统总能量、各个轨道的分子信息。.从头计算法从头计算法(ab initi
18、o methods)以基本物理常数以及元素的原子序数,不借助于任何经验参数,求解薜定谔方程(Schrdinger Equation)计算结果精度高,可靠性大,但是计算量极大,消耗计算机时太多 从头计算法的软件有Gaussian,SPARTAN等等.半经验计算法半经验计算法(semiempirical methods)采用实验值拟合的经验参数,大大提高计算采用实验值拟合的经验参数,大大提高计算速度,计算精度较差。速度,计算精度较差。半经验计算法的软件有:半经验计算法的软件有:MOPAC和和AMPAC等。等。.量子化学计算的优缺点量子化学计算的优缺点可计算出分子的理化参数、分子结构的可计算出分子的
19、理化参数、分子结构的电子电子分布分布和易与亲电试剂或亲核试剂和易与亲电试剂或亲核试剂反应的部反应的部位。位。只适用于计算分子量较小分子,计算时间只适用于计算分子量较小分子,计算时间长长.分子力学分子力学(molecular mechanics)基于经典牛顿力学方程的一种计算分子的平衡结构和能量的方法。在分子中原子间存在化学键,具有标准的键长、键角等,分子内还存在非键相互作用。分子调整自身构象以尽量给出原子位置的最佳排布。由于计算量小,分子力学可研究包括成千上万个原子的分子体系,包括有机小分子、生物大分子。力场总能量,即分子总能量=键合作用键合作用+非键作用非键作用.键合能的构成:.力场的总能量
20、,即分子总能量Etotal:Etotal=Ec+Eb+Et+Ev+Eh+Ee+EdEc 键的伸缩能(compression energy)Eb 键角的弯曲能(bending energy)Et 键的二面角扭转能(torsional energy)Ev 范德华作用能(van de Waals energy)Eh 氢键作用能(hydrogen bonding energy)Ee 静电作用能(electrostatic energy)Ed 偶极作用能(dipole energy).分子力学计算的优缺点分子力学计算的优缺点 可计算出分子的总能量可计算出分子的总能量计算速度快计算速度快 计算优势构象时易
21、于陷于局部最小的能量计算优势构象时易于陷于局部最小的能量.分子动力学分子动力学(molecular dynamics)原子在某一时刻由于运动而其坐标发生原子在某一时刻由于运动而其坐标发生变化。以原子的牛顿运动方程计算每一原变化。以原子的牛顿运动方程计算每一原子的位置、作用力和加速度,按照分子瞬子的位置、作用力和加速度,按照分子瞬时运动状态,模拟分子运动的过程。时运动状态,模拟分子运动的过程。计算体系的计算体系的热力学量热力学量和和其他宏观性质其他宏观性质。.分子处于系列构象的互变平衡之中。分子起功能作用时往往伴随着构象的运动。分子动力学是研究分子构象及其它性质随时间变化的重要工具。计算机模拟的
22、计算机模拟的柔性分子低能构象柔性分子低能构象.构象空间搜寻方法构象空间搜寻方法系统搜寻法系统搜寻法非系统搜寻法非系统搜寻法蒙特卡洛方法蒙特卡洛方法模拟退火法模拟退火法分子动力学法分子动力学法2.4.2 2.4.2 构象空间搜寻方法构象空间搜寻方法.体系能量的变化在一个多维的面上运动,这个面称体系能量的变化在一个多维的面上运动,这个面称为势能面为势能面(energy surface).系统搜寻法系统搜寻法系统地搜寻整个构象空间,寻找能量最低系统地搜寻整个构象空间,寻找能量最低点。在搜寻条件内不会遗漏一个构象。点。在搜寻条件内不会遗漏一个构象。最基本的系统搜寻法是网格搜寻(最基本的系统搜寻法是网格
23、搜寻(gridsearch),也称树状搜寻(),也称树状搜寻(tree search)。)。.非系统搜寻法非系统搜寻法分子动力学法分子动力学法 搜寻构象空间,势能的波动对应着分子搜寻构象空间,势能的波动对应着分子构象的变化,当总能量出现最小值时,求得低能构象。构象的变化,当总能量出现最小值时,求得低能构象。蒙特卡洛方法蒙特卡洛方法(Monte Carlo Method)一种统计抽一种统计抽样方法,在求解的空间中随机采样并计算目标函数,以在样方法,在求解的空间中随机采样并计算目标函数,以在足够多的采样点中找到一个较高质量的最优解作为最终解。足够多的采样点中找到一个较高质量的最优解作为最终解。模拟
24、退火法模拟退火法(simulated annealing algorithm)先升先升温,高温下分子动力学模拟,分子体系有足够的能量,搜温,高温下分子动力学模拟,分子体系有足够的能量,搜寻全部构象空间。在构象空间中选出一些能量相对极小的寻全部构象空间。在构象空间中选出一些能量相对极小的构象,然后逐渐降温,再进行分子动力学模拟,在极小化构象,然后逐渐降温,再进行分子动力学模拟,在极小化后去除能量较高的构象,最后可以得到全局优势构象。后去除能量较高的构象,最后可以得到全局优势构象。.2.4.3 2.4.3 分子三维结构测定分子三维结构测定X射线晶体学射线晶体学核磁共振核磁共振药物分子结构药物分子结
25、构生物大分子结构生物大分子结构生物大分子生物大分子-配体复合物结构配体复合物结构.X射线晶体学射线晶体学(X-ray crystallography)通过对衍射的位置、强度计算,读出原子坐标值,解通过对衍射的位置、强度计算,读出原子坐标值,解析结构得到晶体空间结构析结构得到晶体空间结构,计算机分子模型技术可将上述计算机分子模型技术可将上述数学数值和符号转化为高分辨率的分子。数学数值和符号转化为高分辨率的分子。.化合物测定须得到单晶化合物测定须得到单晶蛋白质结晶蛋白质结晶蛋白质蛋白质-药物复合物结晶药物复合物结晶X射线晶体学测得的是晶体状态下的优势构象射线晶体学测得的是晶体状态下的优势构象.核磁
26、共振技术(NMR)可测定溶液中分子结构优势构象,更好代表生物环境下的可测定溶液中分子结构优势构象,更好代表生物环境下的分子,还能研究动力学特征。分子,还能研究动力学特征。可得到蛋白质或核酸与药物相互作用的构象。可得到蛋白质或核酸与药物相互作用的构象。.2.4.4 2.4.4 分子模拟分子模拟(m molecular olecular m modelingodeling,MMMM)利用计算机图形学进行分子模拟的技术利用计算机图形学进行分子模拟的技术(computer、aided molecular modeling,CAMM)为三维结构研究的有效手段利用计算机利用计算机来构造、显示、分析和贮存复
27、杂的分子模来构造、显示、分析和贮存复杂的分子模型,提供直观的分子立体形象观察分子间型,提供直观的分子立体形象观察分子间的相互作用。的相互作用。.表示分子结构的实体模型表示分子结构的实体模型.分子结构的计算机分子模型分子结构的计算机分子模型特点特点以各种不同的方式显示以各种不同的方式显示在屏幕上随心移动分子在屏幕上随心移动分子 围绕假定的轴旋转围绕假定的轴旋转 翻转并移向某个部位翻转并移向某个部位 可使特定的键转动可使特定的键转动在屏幕上随心操作分子在屏幕上随心操作分子 剪裁、连接、装配、组合、建造、叠合剪裁、连接、装配、组合、建造、叠合计算分子的性质计算分子的性质 原子间距、角度、分子体积、表
28、面积、分子形状原子间距、角度、分子体积、表面积、分子形状 分子的电子特性、氢键、供体分子的电子特性、氢键、供体/受体或带电基团的性质受体或带电基团的性质显示分子间的相互作用显示分子间的相互作用.线型模型计算机分子模型的表示方式计算机分子模型的表示方式.棒型模型棒型模型.球棍模型球棍模型.线型点状模型线型点状模型.CPK模型模型.生物大分子结构的表示方式绳代表肽链无规则卷曲圆柱代表a螺旋缎带中的箭头代表折叠的走向小球代表肽结合的金属原子宽带代表蛋白质分子的 主链.分子网格表面 表示体积.分子表面 表示电性.受体-配体结合模型隐藏大部分受体结构,仅显示结合部位和配体隐藏大部分受体结构,仅显示结合部
29、位和配体.受体-配体结合模型配体与受体以不同形配体与受体以不同形式表示式表示测定配体原子和基团测定配体原子和基团与结合部位邻近原子与结合部位邻近原子和基团的距离,确定和基团的距离,确定相互作用相互作用结合力结合力作用匹配,距离相符作用匹配,距离相符.动态显示受体-配体结合模型.2.4.5 2.4.5 计算机硬件和软件计算机硬件和软件 计算机硬件和工作站 微机 小型计算机图形工作站超级计算机.小型计算机(minicomputer)规模介于微机与中型计算机之间的计算机规模介于微机与中型计算机之间的计算机品牌:品牌:VEX(DEC).计算机图形工作站(graphic workstation)可用于图
30、形等方面具有较佳配备的高档计算机。运算能力可用于图形等方面具有较佳配备的高档计算机。运算能力和显示功能较强和显示功能较强特点:高速特点:高速CPU、高分辨率的彩色显示器、大容量内存、高分辨率的彩色显示器、大容量内存、图形功能强、类似图形功能强、类似UNIX的操作系统的操作系统品牌:品牌:Silicon Graphics Incorporated(SGI),SunGraphics,HP-PA,DEC Alpha,Apollo,Dell.超级计算机(supercomputer)速度快、功能强、规模大速度快、功能强、规模大 品牌:品牌:SGI,曙,曙.软件和专家系统软件和专家系统计算机辅助药物设计软
31、件计算机辅助药物设计软件大多在研究过程中产生的程序的基础上,由商业性公司进大多在研究过程中产生的程序的基础上,由商业性公司进一步作模块化加工,成为实用性软件。一步作模块化加工,成为实用性软件。商业化软件都具有漂亮的界面和使用方便的菜单,方便用商业化软件都具有漂亮的界面和使用方便的菜单,方便用户使用。一般有多种可供选购的功能模块户使用。一般有多种可供选购的功能模块专家系统(专家系统(expert system)应用人工智能技术,模拟学科专家处理专业问题的应用人工智能技术,模拟学科专家处理专业问题的思维、推理方法,作出专家水平的判断和决策,用计算思维、推理方法,作出专家水平的判断和决策,用计算机代
32、替或部分代替专家工作的软件系统。机代替或部分代替专家工作的软件系统。.例:分子设计系列软件Sybyl/Base基本平台基本平台Dynamics分子动力学分子动力学QSAR with CoMFA3D QSARAdvanced Computation3D-构象分析构象分析Biopolymer生物大分子模建生物大分子模建Composer同源蛋白预测同源蛋白预测FlexX对接对接.SYBYL.SYBYL.第三节第三节 基于靶点结构的药物设计基于靶点结构的药物设计-直接法直接法.基于靶点结构的药物设计流程:基于靶点结构的药物设计流程:.主要内容主要内容.3.1 3.1 全新药物设计全新药物设计:根据靶点
33、结构直接造出性状和性质互补的新根据靶点结构直接造出性状和性质互补的新的配体分子的三维结构。由于它能提取结构全新的配体分子的三维结构。由于它能提取结构全新的具有启发性的先导物,故称为全新或从头设计。的具有启发性的先导物,故称为全新或从头设计。全新药物设计的基本流程全新药物设计的基本流程.定义受体活性位点定义受体活性位点1.受体与配体的复合物分析受体与配体的复合物分析 获得受体空腔,分析受体配体结合部位的氨获得受体空腔,分析受体配体结合部位的氨基酸残基,推出可能的药效基团。基酸残基,推出可能的药效基团。2.活性位点分析法活性位点分析法 不知靶点和配体的结合位点,用一些简单的分不知靶点和配体的结合位
34、点,用一些简单的分子或碎片作为探针(水、苯环),找到他们的结子或碎片作为探针(水、苯环),找到他们的结合位置。合位置。.从二氢叶酸还原酶从二氢叶酸还原酶(dihydroflate reductase,DHFR)-甲氨蝶呤甲氨蝶呤(Methotrexate,MTX)复合物获取受体空腔复合物获取受体空腔1、观察空腔表面的氨基酸残基与药物间的作用;、观察空腔表面的氨基酸残基与药物间的作用;2、受体、受体-药物复合物中剥离药物分子。药物复合物中剥离药物分子。DHFR-MTX复合物复合物DHFR 受体空腔图受体空腔图受体与配体的复合物分析受体与配体的复合物分析.活性位点分析法活性位点分析法分析表明受体的
35、活性部位具有分析表明受体的活性部位具有3个重要特征:疏水区、氢键个重要特征:疏水区、氢键接受体区、氢键供体区接受体区、氢键供体区.3.2 生长法生长法 自动模板定位法(自动模板定位法(molecular template-molecular template-directed)directed)原子生长法原子生长法(atom build)(atom build)碎片生长法碎片生长法(fragment build)(fragment build).模板库,与受体性状互补的模版放在活性位点。从有关数据库搜索与受体受点结合的原子或原子团自动模板定位法:自动模板定位法:在受体的结合位点用模版构建出一在
36、受体的结合位点用模版构建出一个形状互补的三维分子骨架。个形状互补的三维分子骨架。.原子生长法原子生长法(atom build)(atom build):在受体的表面逐个增加原子,在受体的表面逐个增加原子,形成与结合位点性状和性质互补的分子形成与结合位点性状和性质互补的分子。起始点种类:起始点种类:1.已知配体或配体的一部分;已知配体或配体的一部分;2.受体结合位点表面形成较强结合的原子受体结合位点表面形成较强结合的原子(锚原子、种子原子锚原子、种子原子)。.碎片生长法碎片生长法(fragment build)(fragment build):每个碎片由单一每个碎片由单一官能团(如羰基、羟基或苯
37、环等)构成,以碎片逐个官能团(如羰基、羟基或苯环等)构成,以碎片逐个生长来构建成一个分子。生长来构建成一个分子。.3.3 连接法连接法 碎片连接法碎片连接法(linked-fragment approach)位点连接法位点连接法(site point connection).碎片连接法:碎片连接法:根据受体结合位点区域性质的不同,分根据受体结合位点区域性质的不同,分为若干能容纳一个碎片的子区域,然后搜寻碎片库进行对接。为若干能容纳一个碎片的子区域,然后搜寻碎片库进行对接。再搜寻连接子库,将各区域的碎片连接起来,形成完整的分再搜寻连接子库,将各区域的碎片连接起来,形成完整的分子。子。搜寻碎片搜寻
38、碎片搜寻连接子搜寻连接子.氢键位点疏水位点在活性位点上匹配碎片在活性位点上匹配碎片在相邻位点上匹配碎片;在相邻位点上匹配碎片;与放好的碎片连接与放好的碎片连接位点连接法:位点连接法:以与靶点活性位点匹配的小分子碎片为以与靶点活性位点匹配的小分子碎片为构建模块,通过与相邻活性位点匹配的碎片进行生长和连接,构建模块,通过与相邻活性位点匹配的碎片进行生长和连接,得到全新的配体分子。得到全新的配体分子。.优点优点:配体结构可能是全新的,不受现有知识的约束,也不受人的思维束缚。该方法既可用于先导化合物的发现,也可用于对先导化合物的结构优化。缺点缺点:可能只是理论分子,仅对接有效,且可能不易合成。全新药物
39、设计的优缺点全新药物设计的优缺点.第四节第四节 基于配体结构的药物设计基于配体结构的药物设计间接法间接法.基于配体结构的药物设计流程:基于配体结构的药物设计流程:.主要内容主要内容.4.1 4.1 药效基团模型法药效基团模型法 对一系列活性化合物进行对一系列活性化合物进行3D-QSAR研研究,并结合构象分析,得出该类药物的药究,并结合构象分析,得出该类药物的药效基团模型。可在受体三维结构未知的情效基团模型。可在受体三维结构未知的情况下,构造出虚拟的药效基团模型,再由况下,构造出虚拟的药效基团模型,再由此出发,进行药物设计此出发,进行药物设计。.药效基团单元药效基团单元pharmacophori
40、c elements 组成药效基团所共有的组成药效基团所共有的一组原子一组原子(团团)。药效基团单元中的原子或官能团可通过药效基团单元中的原子或官能团可通过氢键、氢键、静电作用、范德华作用静电作用、范德华作用和和疏水键疏水键等与受体中受等与受体中受点发生结合。点发生结合。药效基团中的药效基团单元往往用药效基团中的药效基团单元往往用点点来定义,来定义,若干个点构成平面或空间的结构。若干个点构成平面或空间的结构。药效基团是药效基团是药效基团单元药效基团单元的集合的集合.结合部位认定的自动分析方法:利用探针分析活性位点的结构,即用探针原子或基团探测能与靶点产生较好相互作用的部位 主要软件:Sybyl
41、/Receptor;Disco/DiscoTech;Accelrys Catalyst “负片”映射:基于靶点的药效基团识别基于靶点的药效基团识别.基于配体的药效基团识别基于配体的药效基团识别 分子叠合分子叠合.简化方法:最强活性化合物的优势构象简化方法:最强活性化合物的优势构象.简化方法:刚性药物模板简化方法:刚性药物模板右旋氯筒箭毒碱N2受体拮抗剂受体拮抗剂.2.活性类似物法活性类似物法 (Active Analogue Approach,AAA)活性类似物活性类似物共同活性体积共同活性体积(药效基团)(药效基团).例例 活性类似物法(活性类似物法(AAA)确定药效基团)确定药效基团1.测
42、得甲氧氯普胺、可卡因和贝美司琼等共测得甲氧氯普胺、可卡因和贝美司琼等共9个药物的个药物的5-HT3受体活性。受体活性。2.SYBYL计算得到最低能量构象式,重叠得到重叠图。计算得到最低能量构象式,重叠得到重叠图。.3.受体映射,受体映射,包络面为虚拟包络面为虚拟受体结合构象受体结合构象4.以以AAA法推测法推测药效基团单元药效基团单元,获得药效基团获得药效基团模型模型.5.根据药效基团模型设计新分子根据药效基团模型设计新分子.4.2 4.2 基于配体相似性的虚拟筛选基于配体相似性的虚拟筛选 以三维药效基团模型为搜寻标准,搜寻小分子三维结构数据库,得到新配体分子。药效基团识别搜寻化合物库.4.3
43、 2D-QSAR方法方法Two-Dimensional Quantitative Structure-Activity Relationship定量构效关系定量构效关系QSAR 揭示一组化合物的生物揭示一组化合物的生物活性与其分子结构特征之间的相互关系,以数学活性与其分子结构特征之间的相互关系,以数学模型表达和概括出量变规律,以此设计新的化合模型表达和概括出量变规律,以此设计新的化合物。物。.线性自由能相关方法(linear free energy related method),即,即Hansch分析法分析法(Hansch Analysis)log1/C=-k1+k2+k3+k4 Es+ksC:产生指定生物效应的克分子浓度产生指定生物效应的克分子浓度:疏水性参数:疏水性参数:电子效应参数:电子效应参数Es:立体参数:立体参数k1、k2、k3和和k4:权重系数:权重系数ks:常数项常数项.4.4 3D-QSAR方法方法Three-Dimensional Quantitative Structure-Activity Relationships 拟合出药物的理化和三维结构参数与药拟合出药物的理化和三维结构参数与药效的定量关系,以此关系设计新的化合物。效的定量关系,以此关系设计新的化合物。.总结.