第3章药物设计的基本原理和方法(4)

2019/12/211第五节定量构效关系方法简介(QuantitativeStructure–ActivityRelationships)定量构效关系（QuantitativeStructure–ActivityRelationships，QSAR）：1868年，有人提出药物的生物活性与化合物的结构特征呈函数关系：是一种新药设计的研究方法，可以作为先导化合物优化的一种手段。也是计算机辅助药物设计的一个重要内容A＝f（C）生物活性化合物的结构特征由于学科的限制，并没有成功地将此关系用于药物设计2019/12/21220世纪60年代发展建立了三种定量构效关系的研究方法使定量构效关系研究变为可能，并在药物设计中发挥重要作用定量构效关系:并根据信息进一步是指选择一定的数学模式，应用药物分子的物理化学参数、结构参数和拓扑参数及生物活性，对药物分子的化学结构与其生物活性间关系进行定量分析，找出结构与活性间的量变规律，或得到构效关系的数学方程，对药物的化学结构进行优化。2019/12/213分三种方法:第一种方法:Hannsch分析法，1964年由美国Hansch和日本藤田建立该法的特点是以热力学为基础，应用化合物的疏水性参数、电性参数和立体参数表达药物的结构特征，分析结构与生物活性的构效关系。第二种方法Free-Wilson法用数学加和模型表达药物的结构特征，分析其定量构效关系第三种方法Kier分子连接性法根据拓扑学原理将分子连接性指数作为化合物结构的参数结构参数的物理意义明确，便于研究构效关系，设计并预测化合物的活性，应用最为广泛2019/12/214Hansch方法认为，药物能呈现生物活性，是药物小分子与生物大分子相互作用的结果，这种相互作用，与药物的各种热力学性质有关，而且这些热力学性质具有加和性，Hansch分析法又称线性自由能相关模型。Hansch方程的基本通式是：Log1/C=a2+b+c+dEs+···+kLogl/C=alogP2+blogP+c+dEs+···+k上述参数除电性参数外，大多数有加和性，便于计算化合物的生物活性参数化合物的结构特性参数结构的疏水多数结构的电性参数结构立体特征参数2019/12/215Hansch方程有时还使用I作为指示变量，I是一种半定量的参数，只有1和0两个值。I用来描述不能用连续变量表达的药物分子特征一、疏水性参数（LipophilicityParameters）分子疏水性参数IogP,即分子的脂水分配系数(partitioncoefficient)，表示分子的疏水性化合物在有机相和水相中分配平衡时的量（摩尔）浓度Co和CW之比值,P=CO/CWP值一般较大，常用IogP表示当分子中有该取代基时I为1，当分子中没有该取代基时I为0。2019/12/216logP的测定:用经典的摇瓶法测定，为最准确的数据logP还能用HPLC的方法对一些已知此合物，可以从Hansch的数据手册或计算机辅助药物设计工作站的数据库查到logP的有关数据。目前最快速的方法是在计算机工作站上建构化合物的二维结构，通过分子动力学计算优化，得到三维优势构象，用CLOGP商业软件模块计算，可自动得到logP的数据研究对象是同源化合物时，由于具有相同的基本母核，可用取代基疏水常数(Substituenthydrophobicconstant)π值作疏水参数。更多的π值可查Hansch的数据手册取代基BrCIFNO2HOHNH2CH3C2H5C6H5π值0.860.710.14-0.280.00-0.67-1.230.561.021.96对氢来说，πH=0取代基疏水常数的优点是可直接查表得到2019/12/217π值可直接作疏水参数，也可以通过加和的方法，计算化合物logPLogp=logPH+ΣπX+ΣFX母体化合物的logP母体上各取代基为值的总和需进行校正的各种因素之和二、电性参数（ElectronicParameters）可用查表方法直接得到的电性参数有：、、、*等Hammett常数，表示芳香族取代基的诱导和共轭效应之和取代基的诱导效应参数，对值则是取代基的共轭效应参数各种与电性性质有关的物化参数如仍极矩（μ）、解离常数值等也均是常用的电性参数2019/12/218红外、紫外、NMR、MS等光谱数据，也可用来作构效关系研究的电性参数三、立体参数（StericParameters）该类参数的种类也很多，主要是：（1）取代基的taftEs参数，这是从化学反应导出的立体参数。Es=logKx-logKH（2）摩尔折射率(MR，molarrefractivity)，近似代表分子体积（3）VanderWaals体积，可直接描述取代基以及化合物的体积大小（4）STERIMOL多维立体参数（verloop参数），是研究药物与受体结合大有用的立体参数，它可直接描述取代基在三维空间上的立体信息。2019/12/219四、Hausch方法在药物设计中的应用（ApplicationofHanschMethodinDrugDesign）在新药设计中应用Hansch方法的条件：1所有化合物应是同源物，结构上具有相同的基本母核，与同一个受体作用具有相同的作用机理；2所设计化合物的物理化学性质差异要大所选择参数之间不能有相关性，要有比较大的差异，并且生物活性数据的变化幅度应大于一个对数单位（即大于10倍），否则得不到足够的信息；2019/12/2110Hansch方程除了研究定量构效关系外，还能用来解释药物作用机理，推测和描述可能的受体模型，研究除活性以外的其他药代动力学定量关系Hansch方法的一般操作过程分五个步骤：第一，从先导化合物出发，设计并合成首批化合物；第二，用可靠的定量方法测活性；第三，确定及计算化合物及取代基的各种理化参数或常数；第四，用计算机程序计算Hansch方程，求出一个或几个显著相关的方程；第五，用所得方程定任地设计第二批新的化合物，并预测活性2019/12/21111Hansch方法是二维的QSAR研究方法，该方法只考虑了化合物与受体作用的位点，没有考虑化合物与受体(酶)的结合时构象的变化Hansch方法的主要缺陷：2只用于描述药物分子的二维结构，不能定量地描述三维结构与生物活性间的关系。3只能优化先导化合物，不能发现先导化合物。所有参数只能表达二维意义上的结构特点，不能研究与受体三维空间作用的情况，不能研究药物构象和构型对活性的影响，因此不能全面解释生物活性的本质。2019/12/2112第六节计算机辅助药物设计简介（Computer-AidedDrugDesign,CADD）三维定量构效关系（3D-QSAR）是以药物分子和受体分子的三维结构特征为基础，分析结构与生物活性间的定量关系20世纪80年代随着计算机辅助分子模拟技术（computeraidedmolecularmodeling）的发展，计算机辅助药物设计（Computer-AidedDrugDesign，CADD）得以迅速发展。CADD是利用计算机的快速计算功能、全方位的逻辑推理功能、一目了然的图形显示功能，将量子化学、分子力学、药物化学、生物学科、计算机图形学和信息科学等学科交叉融会结合，从药物分子的作用机理入手来药物设计，减少了盲目性，节省了大量的人力和物力。2019/12/2113从药物的作用机理，从靶点出发，考虑药物与受体的作用过程，并要模拟药物在体内吸收、转运、代谢等过程，比基于结构的药物设计更合理，但目前该法不成熟，还没有成功的例子。CADD有两类方法：基于机理的药物设计(MechanismBasedDrugDesign，MBDD)基于结构的药物设计（StructureBasedDrugDesign，SBDD）研究方法较为成熟，广为应用，本节重点介绍该法。又分为两类直接药物设计间接药物设计又称全新药物设计2019/12/2114一、直接药物设计（denovoDrugDesign）直接药物设计又称全新药物设计受体靶点的三维结构是已知该法的设计思路以生物大分子的三维结构为基础，根据受体受点的形状和性质，研究药物与受体的相互作用，设计新的药物受体的三维结构可通过X-单晶衍射得到，或通过多维NMR间接得到用计算机分子模拟技术研究受体与药物结合部位的性质，如静电场、疏水场、氢键作用多位点的信息，根据靶标分子与药物分子相结合的活性部位的几何形状和化学特征，运用数据库搜寻设计与受体作用部位的形状和理化性质相匹配的分子。该法既能设计新的先导化合物，也能优化先导化合物。2019/12/2115受体是生物体的特异性大分子，药物小分子称为配体（ligand）。在产生药理作用时，配体首先要分布到受体部位，并与受体结合（binding）。受体与配体结合的部位（bindingShe）是计算机辅助药物设计的重点研究问题，药物与受体作用并不是整体作用，实际是涉及受体中一部分重要的氨基酸残基。直接药物设计常用的方法分子对接法（Docking）活性位点分析法(ActiveSiteAnalysis，ASA)首先要建立拥有大量化合物的三维结构数据库，用DOCK程序，把库中的分子逐一与靶标分子进行“对接”(Docking)，通过不断优化小分子化合物的结构，寻找小分子与靶标大分子作用的最佳构象，计算其相互作用能，并在三线数据库搜寻与受体受点区域中心相匹配的最佳分子。首先要分析哪些原子和基团可以和受体活性部位有较好的相互作用，把这些原子和基团作为碎片，连接成分子模板定位法通过药效团定出一个药物与受体作用的模板2019/12/2116通过生长、旋转等得到基本骨架，按照受体的腔穴，定出靶标边界，这是一级结构的生成。从有关数据库搜索与受体受点结合的原子或原子团，设计新的化合物2019/12/2117二、间接药物设计（IndirectDrugDesign）间接药物设计法受体的三维结构并不清楚以小分子的构效关系为基础，从一组小分子化合物的结构和生物活性数据出发，研究结构与活性关系的规律该法的设计思路二维定量构效关系是一种重要的间接药物设计方法三维定量构效关系是研究药物的三维结构，以活性好的受体搞抗剂或激动剂为依据，通过计算分子的优势构象来代替药效构象，推测受体的活性部位形状及作用方式等信息，得到虚拟受体模型，通过模拟药物与受体的结合，设计新的配基分子。20世纪90年代发展起来2019/12/2118三维定量构效关系(3DQSAR)研究的方法有三种分子形状分析法其理论基础是认为药物分子有多种构象。每种构象视为一种形状，药物活性与这些形状对受体的活性部位（空穴）的匹配性有关。距离几何学方法药物与受体的相互作用是通过药物的活性基团与受体受点直接作用产生，其活性强弱与药物与受点的结合能大小有关，通过选择合理的受体结合点分布模型，建立药物活性基团和受体结合点类型相关的能量参数。先构造所研究药物的三维分子模型，计算药物的优势构象，求出各原子间距离的矩阵(S)，用若干个矩阵表示整个分子的模型。根据活性最高的药物构象，设计一些“空”的结合位点，再根据活性最低的构象，设计一些“实”的结合位点，从距离矩阵的上下限出发，进行搜寻，利用结合能参数，寻找药物基团与受体结合点的模型该方法的操作过程是，先计算各药物分子的优势构象，一般以活性最强者为模板，其他分子与其合理叠会，求出分子形状参数L0（重叠长度）、S0（重塑面积）、V0（重叠体积）等，用分子形状参数进行QSAR计算，分析受体活性部位形状，并根据这些参数设计新的化合物。比较分子场分析法SeeNextSlide2019/12/2119可采用化合物周围的静电场、范德华力场、氢键场的空间分布作为化合物结构特征变量，用数学方法建立化合物的生物活性与化合物周围各力场空间分布之间关系的模型。比较分子场分析法(Comparativemolecularfieldanalysis，CoMFA)当药物与同一受体产生相互作用时，主要通过静电作用、疏水作用和范德华力等非共价键作用。认为分子与受体问的这些力场有某种相似性，活性与力场的大小及方向相关。可用分子力