第3章药物设计的基本原理和方法4.ppt
2023/11/3,1,第五节定量构效关系方法简介(Quantitative Structure Activity Relationships),定量构效关系（Quantitative Structure Activity Relationships，QSAR）：,1868年，有人提出药物的生物活性与化合物的结构特征呈函数关系：,是一种新药设计的研究方法，可以作为先导化合物优化的一种手段。也是计算机辅助药物设计的一个重要内容,Af（C）,生物活性,化合物的结构特征,由于学科的限制，并没有成功地将此关系用于药物设计,2023/11/3,2,20世纪 60年代,发展建立了三种定量构效关系的研究方法,使定量构效关系研究变为可能，并在药物设计中发挥重要作用,定量构效关系:,并根据信息进一步,是指选择一定的数学模式,，应用药物分子的物理化学参数、,结构参数,和拓扑参数,及生物活性，,对药物分子的化学结构,与其生物活性间关系进行定量分析，,找出结构与活性间的,量变规律，,或得到构效关系的数学方程，,对药物的化学结构进行优化。,2023/11/3,3,分三种方法:,第一种方法:Hannsch分析法，1964年由美国Hansch和日本藤田建立,该法的特点是以热力学为基础，应用化合物的疏水性参数、电性参数和立体参数表达药物的结构特征，分析结构与生物活性的构效关系。,第二种方法 Free-Wilson法,用数学加和模型表达药物的结构特征，分析其定量构效关系,第三种方法 Kier分子连接性法,根据拓扑学原理将分子连接性指数作为化合物结构的参数,结构参数的物理意义明确，便于研究构效关系，设计并预测化合物的活性，应用最为广泛,2023/11/3,4,Hansch方法认为，药物能呈现生物活性，是药物小分子与生物大分子相互作用的结果，这种相互作用，与药物的各种热力学性质有关，而且这些热力学性质具有加和性，Hansch分析法又称线性自由能相关模型。,Hansch方程的基本通式是：,Log 1/C=a2+b+c+dEs+kLog l/C=a logP2+blogP+c+dEs+k,上述参数除电性参数外，大多数有加和性，便于计算,化合物的生物活性参数,化合物的结构特性参数,结构的疏水多数,结构的电性参数,结构立体特征参数,2023/11/3,5,Hansch方程有时还使用I作为指示变量，I是一种半定量的参数，只有1和0两个值。,I用来描述不能用连续变量表达的药物分子特征,一、疏水性参数（Lipophilicity Parameters）,分子疏水性参数IogP,即分子的脂水分配系数(partition coefficient)，表示分子的疏水性,化合物在有机相和水相中分配平衡时的量（摩尔）浓度Co和CW之比值,P=CO/CW,P值一般较大，常用IogP表示,当分子中有该取代基时I为1，当分子中没有该取代基时I为0。,2023/11/3,6,logP的测定:,用经典的摇瓶法测定，为最准确的数据,logP还能用HPLC的方法,对一些已知此合物，可以从Hansch的数据手册或计算机辅助药物设计工作站的数据库查到logP的有关数据。,目前最快速的方法是在计算机工作站上建构化合物的二维结构，通过分子动力学计算优化，得到三维优势构象，用CLOGP商业软件模块计算，可自动得到logP的数据,研究对象是同源化合物时，由于具有相同的基本母核，可用取代基疏水常数(Substituent hydrophobic constant)值作疏水参数。,更多的值可查Hansch的数据手册,对氢来说，H=0取代基疏水常数的优点是可直接查表得到,2023/11/3,7,值可直接作疏水参数，也可以通过加和的方法，计算化合物logP,Logp=logPH+X+FX,母体化合物的lo g P,母体上各取代基为值的总和,需进行校正的各种因素之和,二、电性参数（Electronic Parameters）,可用查表方法直接得到的电性参数有：、*等,Hammett常数，表示芳香族取代基的诱导和共轭效应之和,取代基的诱导效应参数，对值则是取代基的共轭效应参数,各种与电性性质有关的物化参数如仍极矩（）、解离常数值等也均是常用的电性参数,2023/11/3,8,红外、紫外、NMR、MS等光谱数据，也可用来作构效关系研究的电性参数,三、立体参数（Steric Parameters）,该类参数的种类也很多，主要是：,（1）取代基的taft Es参数，这是从化学反应导出的立体参数。Es=logKx-logKH,（2）摩尔折射率(MR，molar refractivity)，近似代表分子体积,（3）Van der Waals体积，可直接描述取代基以及化合物的体积大小,（4）STERIMOL多维立体参数（verloop参数），是研究药物与受体结合大有用的立体参数，它可直接描述取代基在三维空间上的立体信息。,2023/11/3,9,四、Hausch方法在药物设计中的应用（Application of Hansch Method in Drug Design）,在新药设计中应用Hansch方法的条件：,1 所有化合物应是同源物，结构上具有相同的基本母核，与 同一个受体作用具有相同的作用机理；,2 所设计化合物的物理化学性质差异要大,所选择参数之间不能有相关性，要有比较大的差异，并且生物活性数据的变化幅度应大于一个对数单位（即大于10倍），否则得不到足够的信息；,2023/11/3,10,Hansch方程除了研究定量构效关系外，还能用来解释药物作用机理，推测和描述可能的受体模型，研究除活性以外的其他药代动力学定量关系,Hansch方法的一般操作过程,分五个步骤：,第一，从先导化合物出发，设计并合成首批化合物；,第二，用可靠的定量方法测活性；,第三，确定及计算化合物及取代基的各种理化参数或常数；,第四，用计算机程序计算Hansch方程，求出一个或几个显 著相关的方程；,第五，用所得方程定任地设计第二批新的化合物，并预测活性,2023/11/3,11,1 Hansch方法是二维的QSAR研究方法，该方法只考虑了化合物与受体作用的位点，没有考虑化合物与受体(酶)的结合时构象的变化,Hansch方法的主要缺陷：,2 只用于描述药物分子的二维结构，不能定量地描述三维结构与生物活性间的关系。,3 只能优化先导化合物，不能发现先导化合物。,所有参数只能表达二维意义上的结构特点，不能研究与受体三维空间作用的情况，不能研究药物构象和构型对活性的影响，因此不能全面解释生物活性的本质。,2023/11/3,12,第六节 计算机辅助药物设计简介（Computer-Aided Drug Design,CADD）,三维定量构效关系（3D-QSAR）是以药物分子和受体分子的三维结构特征为基础，分析结构与生物活性间的定量关系,20世纪80年代随着计算机辅助分子模拟技术（computer aided molecular modeling）的发展，计算机辅助药物设计（Computer-Aided Drug Design，CADD）得以迅速发展。,CADD是利用计算机的快速计算功能、全方位的逻辑推理功能、一目了然的图形显示功能，将量子化学、分子力学、药物化学、生物学科、计算机图形学和信息科学等学科交叉融会结合，从药物分子的作用机理入手来药物设计，减少了盲目性，节省了大量的人力和物力。,2023/11/3,13,从药物的作用机理，从靶点出发，考虑药物与受体的作用过程，并要模拟药物在体内吸收、转运、代谢等过程，比基于结构的药物设计更合理，但目前该法不成熟，还没有成功的例子。,CADD有两类方法：,基于机理的药物设计(Mechanism Based Drug Design，MBDD),基于结构的药物设计（Structure Based Drug Design，SBDD）,研究方法较为成熟，广为应用，本节重点介绍该法。,又分为两类,直接药物设计,间接药物设计,又称全新药物设计,2023/11/3,14,一、直接药物设计（de novo Drug Design）,直接药物设计又称全新药物设计受体靶点的三维结构是已知,该法的设计思路,以生物大分子的三维结构为基础，根据受体受点的形状和性质，研究药物与受体的相互作用，设计新的药物,受体的三维结构可通过X-单晶衍射得到，或通过多维NMR间接得到,用计算机分子模拟技术研究受体与药物结合部位的性质，如静电场、疏水场、氢键作用多位点的信息，根据靶标分子与药物分子相结合的活性部位的几何形状和化学特征，运用数据库搜寻设计与受体作用部位的形状和理化性质相匹配的分子。该法既能设计新的先导化合物，也能优化先导化合物。,2023/11/3,15,受体是生物体的特异性大分子，药物小分子称为配体（ligand）。在产生药理作用时，配体首先要分布到受体部位，并与受体结合（binding）。受体与配体结合的部位（binding She）是计算机辅助药物设计的重点研究问题，药物与受体作用并不是整体作用，实际是涉及受体中一部分重要的氨基酸残基。,直接药物设计常用的方法,分子对接法（Docking）,活性位点分析法(Active Site Analysis，ASA),首先要建立拥有大量化合物的三维结构数据库，用 DOCK程序，把库中的分子逐一与靶标分子进行“对接”(Docking)，通过不断优化小分子化合物的结构，寻找小分子与靶标大分子作用的最佳构象，计算其相互作用能，并在三线数据库搜寻与受体受点区域中心相匹配的最佳分子。,首先要分析哪些原子和基团可以和受体活性部位有较好的相互作用，把这些原子和基团作为碎片，连接成分子,模板定位法,通过药效团定出一个药物与受体作用的模板,2023/11/3,16,通过生长、旋转等得到基本骨架，按照受体的腔穴，定出靶标边界，这是一级结构的生成。从有关数据库搜索与受体受点结合的原子或原子团，设计新的化合物,2023/11/3,17,二、间接药物设计（Indirect Drug Design）,间接药物设计法受体的三维结构并不清楚,以小分子的构效关系为基础，从一组小分子化合物的结构和生物活性数据出发，研究结构与活性关系的规律,该法的设计思路,二维定量构效关系是一种重要的间接药物设计方法,三维定量构效关系是研究药物的三维结构，以活性好的受体搞抗剂或激动剂为依据，通过计算分子的优势构象来代替药效构象，推测受体的活性部位形状及作用方式等信息，得到虚拟受体模型，通过模拟药物与受体的结合，设计新的配基分子。,20世纪90年代发展起来,2023/11/3,18,三维定量构效关系(3DQSAR)研究的方法有三种,分子形状分析法,其理论基础是认为药物分子有多种构象。每种构象视为一种形状，药物活性与这些形状对受体的活性部位（空穴）的匹配性有关。,距离几何学方法,药物与受体的相互作用是通过药物的活性基团与受体受点直接作用产生，其活性强弱与药物与受点的结合能大小有关，通过选择合理的受体结合点分布模型，建立药物活性基团和受体结合点类型相关的能量参数。,先构造所研究药物的三维分子模型，计算药物的优势构象，求出各原子间距离的矩阵(S)，用若干个矩阵表示整个分子的模型。根据活性最高的药物构象，设计一些“空”的结合位点，再根据活性最低的构象，设计一些“实”的结合位点，从距离矩阵的上下限出发，进行搜寻，利用结合能参数，寻找药物基团与受体结合点的模型,该方法的操作过程是，先计算各药物分子的优势构象，一般以活性最强者为模板，其他分子与其合理叠会，求出分子形状参数L0（重叠长度）、S0（重塑面积）、V0（重叠体积）等，用分子形状参数进行QSAR计算，分析受体活性部位形状，并根据这些参数设计新的化合物。,比较分子场分析法,SeeNextSlide,2023/11/3,19,可采用化合物周围的静电场、范德华力场、氢键场的空间分布作为化合物结构特征变量，用数学方法建立化合物的生物活性与化合物周围各力场空间分布之间关系的模型。,比较分子场分析法(Comparative molecular field analysis，CoMFA),当药物与同一受体产生相互作用时，主要通过静电作用、疏水作用和范德华力等非共价键作用。,认为,分子与受体问的这些力场有某种相似性，活性与力场的大小及方向相关。,可用分子力场模型描述受体间的相互作用。,2023/11/3,20,CoMFA方法的一般操作过程,（1）逐一建造各化合物的分子，进行结构优化，研究药物的优势构象，并计算各原子的电荷密度。,（2）按照合理的叠加规则，把优势构象重叠在一个能包容全部化合物的空间网格上（lattice）。图3-9右是在药物设计工作站上31个甾体化合物的叠合的显示。,对活性的影响,2023/11/3,21,(3)根据化合物分子与受体的作用方式，选择合适的探针(Probe)，在空间同格中移动。所选择的探针，取决于被考察的力场的特征。,用H2O分子做深针，可研究药物与受体的流水作用和氢键作用,用CH3探针代表范德华力力场，可研究药物与受体的立体作用,用H+探针体现静电力场，研究药物与受体的电性作用,（4）探针每移动一个步长，计算其在空间网格上与各原子的相互作用能量,（5）建立等高线（contour maps），如图3-9右所示。该图形用红、蓝、黄、绿4种不同颜色，表达化合物立体和电性性质。,目前的计算机辅助药物设计还主要是基于结构的药物设计，今后的目标是向着基于机理的药物设计方向发展,已有成功的例子，如青蒿素 p60-61,
