医学基于结构的药物设计专题课件.ppt

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