药物设计的基本原理和方法幻灯片.ppt

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