计算机辅助药物设计(完整版).pdf

《计算机辅助药物设计(完整版).pdf》由会员分享，可在线阅读，更多相关《计算机辅助药物设计(完整版).pdf（11页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、 计算机辅助药物设计完整版 第1章概论 一、药物发现一般过程 新药的研究有三个决定阶段：机辅助药物设计的主要任务就是先导化合物的发现与优化。二、合理药物设计 1、合理药物设计（rational drug design）是依据与药物作用的靶点，即广义上的受体，如酶、受 体、离子通道、病毒、核酸、多糖等，寻找和设计合理的药物分子。通过对药物和受体的结构在分 子水平甚至电子水平的全面准确了解进行 作用方式及产生生理活性的机理的认识 子作为效应子作用于机体的靶点，考察其 效应及运动协调性等进行分子设计。2、方法分类（1）合理药物设计有基于靶点结构的 三维结构搜索 和全新药物设计 等方法。后者分为模板定

2、位 法、原子生长法、分子碎片法。（2）根据受体是否已知分为直接药物设计和间接药物设计。前者即通过结构测定已知受体或受 体-配体复合物的三维结构，根据受体的三维结构要求设计新药的结构。受体结构测定方法：同源模 建（知道氨基酸序列不知道空间结构时），X 射线衍射（可结晶并得到晶体时），多维核磁共振技术（在体液即在水溶液环境中）。后者通过一些配体的结构知识（SAR，计算机图形显示等）推测受体 的图像，提出假想受体，采用建立药效团模型或 3D-QSAR 和基于药效团模型的三维结构搜索等方 法，间接进行药物设计。三、计算化学 计算化学包括分子模型、计算方法、计算机辅助分子设计（计。计算方法基本上可分为两

3、大类：分子力学（采用经典的物理学定律只考虑分子的核而忽略外围 的电子）和量子力学（采用薛定谔方程考虑外围电子的影响，分为 从头计算方法和半经验方法）。常用的计算应用有：（1）单点能计算：根据模型中原子的空间位置给出相应原子坐标的势能；（2）几何优化：系统的修改原子坐标使原子的三维构象能量最小化；（3）性质计算：预测某些物理 化学性质，如电荷、偶极矩、生成热等；（4）构象搜索：寻找能量最低的构象；（5）分子动力学模 拟：模拟分子的构象变化。方法选择主要有三个标准：（1）模型大小；（2）可用的参数；（3）计算机资源 四、计算化学中的基本概念 1、坐标系统 分为笛卡尔坐标（三维空间坐标）和 内坐标（

4、Z 矩阵表示，参数为键长、键角、二面角数据）。前者适合于描述一系列的不同分子，多用于分子力学程序，有 3N 个坐标；后者常用于描述单分子 系统内各原子的相互关系，多用于量子力学程序，有 3N-6 个坐标。2、原子类型：用来标记原子属性。3、势能面 体系能量的变化被认为能量在一个多维的面上运动，这个面被称为势能面。坐标上能量的一阶 导数为零的点为定点（原子力为零，局部或全局最稳定）4、面积 Van der Waals 面积：原子以 van der Waals 为半径的球的简单堆积。先导化合物的发现，新药物的优化研究，临床与开发研究。计算 基于结构的药物设计 基于机理的药物设计。形状互补，性质互补

5、 和通过对靶点的结构、功能、与药物 CADD 通过内源性物质或外源性小分（包括氢键、疏水性、静电等），溶剂 CAMD）、化学数据库及有机合成设 分子面积：试探分子（常为半径 1.4?的水分子）在 Van der Waals 面积上滚动的面积（包括试 探球与分子的接触面积和分子空穴产生的悬空面积）。可接近面积：试探球在分子 van der Waals 表面滚动时试探球原点处所产生的面积。5、单位：键长多用？（埃，angstroms），键能多用 kcal/mol 表示。五、计算机辅助药物设计软件及限制 目前 CADD 存在的问题：蛋白质结构三维结构的真实性和可用性问题（细胞膜上的受体或跨膜 蛋白离

6、开原先环境，空间排列会发生很大变化，难以得到真实的三维空间结构；大量受体结构未知；很多受体只有一级结构，获得的三维结构有限）；受体-配体相互作用的方式问题；设计的分子能否 进行化学合成；药物体内转运、代谢和体内毒副作用问题等。第2章分子力学 分子力学 是基于原子间存在化学键、非键原子之间的范德华及静电相互作用这一经典理论，通 过分子几何、能量、振动光谱及其他物理性质的计算寻求分子的平衡构型及能量，确定有机分子的 结构、构象、能量及动力学模型。其计算 忽略了电子的贡献，只考虑核。计算较小仅与分子中原子 数目的平方成正比。一般分子动力学软件提供三种位能面采样算法：单点：只是对位能面上某一点计算，给

7、出该构象下的系统能量和 点在位能曲线上陡度）几何优化：对单点位能面采样，寻找梯度为零的构象，局部最小 分子动力学：对势能面增加动能，导致分子系统按 Newton 定律运动，在势能低点运动加快。主要用于能量最小化和构象搜索。受体结构已知，该法计算药物与受体的结合能；受体未知通过已 知配体导出药效团模型。一、理论简介 分子力学基本思想是通过选择一套 势函数和从实验中得到的一套力常数，从给定的分子体系原 子的空间坐标的初值，用分子力场描述的体系总能量对于原子坐标的梯度，通过多次迭代的数值算 法来得到合理的分子体系的结构。分子的化学键具有一定的键长、键角，分子要调整它的几何形状（构象），必须使其键长值

8、和键 角值尽可能接近标准值，同时非键作用能处于最小的状态，由这些键长和键角调节构象，给出核位 置的最佳分布，即分子的平衡构型。分子力学优化只能是 局部优化，若为了找到全局能量最低构象，须将所有可能的初始构象分别 进行优化，最后进行比较确定分子体系的最优构象。二、分子力场 分子力学有能力处理大分子体系，它从经典力学的观点来描述分子中原子的拓扑结构，是通过 分子立场这个分子模拟的基石实现的。如果一个解析表达式能拟合位能面，则此解析表达式就成为分子力场，亦即 一个力场的确定就 是选择解析函数形式和确定参数。分子力学用几个典型 结构参数和作用力来描述结构的变化，由 分 子内相互作用（键伸缩，角弯曲，扭

9、转能，面外弯曲等）和 分子间相互作用（静电、氢键、vdW）构成。三、能量最小化 按是否采用能量的导数分为两类：非导数法（即单纯型法：以逐个改变原子的位置来寻找能量 最小值，找到的并不一定是局部最小值，主要用于调整分子的起始构象）和导数法。一阶能量导数 的方向指向能量最小化的点，梯度反映该点的陡度，有 最陡下降法（SD）、共轭梯度法、任意步长 逼近法；二阶能量导数预测何处能量梯度方向发生变化，有 牛顿-拉普森法。SD:梯度是进行搜索的方向，每次搜索之后旧的方向被新点处的梯度取代，适合优化最初段，梯度（反应能量下降方向上该 尤其是减少大量的非键相互作用非常有用，适用于大分子。共轭梯度法：不仅运用当

10、前梯度，也采 用先前的最小化历史来确定下一步，收敛比 SD 快，用于大分子。Newt on-Ra phso n法：原则上可一 步收敛，但存储导数的矩阵太大，不适用于大体系。四、常用的分子力场 以适合生物大分子的 Amber 和适合小分子的 MM2 为代表。MM3:MM 力场对静电的相互作用采用键的偶极方法，对于极性或电荷系统不能充分模拟，适 用于小非极性分子的结构和热动力学模拟。AMBER 力场：广泛用于蛋白质和核酸，不适合用于小分子。OPLS 力场：用于蛋白质和核酸，特别适用在液相系统中模拟物理性质。CHARMM 力场：适用于生物大分子，充分考虑溶剂和溶剂、溶质和溶质、溶剂和溶质之间的 相互

11、作用。BIO+力场：CHARMM 力场的补充，采用 CHARMM 立场的参数，结果与 CHARMM 一样。MMFF94 力场：运用凝聚态过程，适合大分子和小分子，且精确一致。Universal 力场：针对整个周期表的分子力学和动态模拟力场，有过渡元素时的最佳选择。COMP ASS 力场 五、分子动力学 以牛顿力学为基础，把每个原子看做符合牛顿运动定律的粒子，在一定时间内，连续几分牛顿 运动方程计算原子的位置和速度得出原子的运动轨道。分子动力学涉及 Newt on 运动方程的积分，需 要选取适当的时间步长，选取的时间步长和运动的频率有关。分子动力学模拟经过三个阶段：加热、模拟（包括平衡期和资料收

12、集期）、冷却。构建分子的时候是 0K（原子运动速度为 0），缓慢加热（使系统在每个步长里都接近平衡，在 较短时间达到模拟温度下的平衡）到模拟温度。对室温下的模拟，梯度应小于 3，是为了避免人为 的在高能区产生的局部力使分子在张力较大的地方发生断裂或扭曲。在模拟期的平衡很重要，可以 避免加热过程中引入人为因素。分子体系冷却可降低在较高温度时分子的张力，冷却过程也叫 退火，使分子从高能构象越过一定的能垒转向稳定的低能构象。分子动力学模拟溶剂的作用，可通过选择媒介的介电常数和周期边界条件模拟，水的介电常数 80.4。采用周期边界条件模拟是为了 消除刚性壁边界条件和自然边界条件的表面效应。去顶八面体

13、常用于球形分子。分子力学模拟一般都应先对分子进行几何优化，在几何优化基础上再对分子进行动力学模拟。分子动力学的作用：（1）分子动力学主要用于 能量最小化，和分子力学的区别：分子力学不能越过一定的能垒，只 是局部优化；分子动力学模拟则是全局优化低能构象（越过能垒高度与模拟温度有关）。（2）通过分子动力学模拟，可再现分子的各种构象形式，用于推测药物与受体相作用的构象。（3）分子力学适合处理分子内张力（键，角，二面角）或 当考虑到极性溶剂（如体内）或溶剂效应时，使用分子动力学方法，力学或模拟退火法。Monte Carlo 方法 利用随机取样处理问题的方法称为 Monte Carlo 方法，它是一种通

14、过的采取随机数和概率统计进 行猜测来研究问题。分子动力学不能越过的能垒，Monte Carlo 构象搜寻可以是跳跃式的，其优点是 取样的构象恰当，对低能构象取样几率大。七、Lan gev in Langevin 模拟是随机动力学模拟，通过给各原子分配分子在溶剂中与其他分子相互碰撞并随着 在溶剂中运动产生的摩擦力的值进行模拟。Langevin 动力学模拟只是对分子施加一个力模拟碰撞后 的能量损失，所以不需要指明溶剂分子，多用于长链分子和聚合物。特别适用于研究溶剂中的大分 子。八、构象分析 模拟 vdw 力等分子处于非极性溶剂中，Monte carlo 方法,Langevin 动 八、描述分子结构

15、的三个层次：分子构造，分子构型，分子构象。构象搜寻采用适当的方法产生各 种不同的构象，并对这些构象进行能量最小化，比较这些构象并找出其中能量最低的构象。根据产 生构象的方法不同，可分为（1）系统搜寻法（系统地搜寻分子的争购构象空间，寻找势能面上的极 小点。最基本的搜索方法是 格点搜索，即在分子构象空间中以柔性键的旋转角度小间隔为变量逐渐 搜索。只适合处理小分子体系，也不适用于环状结构。）和（2）非系统搜寻两类。九、随机搜寻法（1）模拟退火方法 是分子动力学，Monte Carlo 和 Langevin 动力学在模拟时采取温度缓慢降低的方法。它首先使体 系升温，使分子体系有足够的能量，克服柔性分

16、子中存在的各种旋转能垒和顺反异构能垒，搜寻全 部构象空间，在构象空间中选出一些能量相对极小的构象，然后逐步降温，再进行分子动力学模拟，此时较高能垒已无法越过，在极小化后去除能量较高的构象，最后可以得到能量最小的优势构象。Monte Carlo 退火方法采用 Monte Carlo 方法的 Metropolis 采样算法。在某温度下，体系有起始 构象，构象发生微小的随机变化产生新的构象，相应的能量发生变化。如果能量差不大于零则接受 构象变化，新构象成为下一步的初始构象。如果能量大于零，选择一个随机数将其与原相比较，若 能量小于零则接受，新的构象成为下一次随机变化的起始点；否则拒绝变化，老的构象仍

17、是下一次 随机变化的起始点（过程中产生 Boltzmann 分布）。（与 Monte Carlo 的区别在于 Monte Carlo 模拟退 火方法中温度也是体系的变量。模拟退火方法可有效地寻找分子的优势构象，因为其取舍构象时不仅接受能量下降的变化，也 接收部分部分能量上升的变化，并且该法不依赖于初始构象。（2）高温淬火动力学 是高温分子动力学与能量最低化相结合来判断一系列构象的分布，往往还需要结合模拟退火。（3）遗传算法（主要用于全局优化）算法步骤：随机产生初始群体，群体中个体以二进制序列标记计算适应值（个体优劣的度 量和下一代存活概率）通过复制算子、杂交算子和变异算子产生新一代更具适应的群

18、体。对于构象搜寻，染色体二进制序列的值表示分子可旋转的扭角。适应值为能量的函数。（4）距离几何法 该方法核心仍为随机技术。特别适用于导出大量信息无法手工解析的蛋白质与核酸的结构。第3章量子力学 一、量子化学理论简介 应用量子力学原理处理化学问题，形成分子轨道理论、价键理论与配位场理论。分子轨道理论（又称 Hartee-Fock 理论）在物理模型上有三个基本近似：非相对近似是电子在 原子核附近运动而不被的原子核俘获，必须保持很高运动度。近似认为电子质量等于电子静止的质 量，即电子质量恒为1 个单位。Born-Oppenheimer 近似是在计算分子总能量（电子总能量与核排 斥能值和）时，把电子的

19、运动和核的运动分开处理，忽略贡献很小的电子态之间偶合项，即非绝热 项，故又称绝热近似。单电子近似则认为 电子波函数为 n 个电子所占据的轨道（单电子函数）的 乘积。自洽场：每一个电子的运动状态不仅受核力的作用，而且还要受其他 n-1 个电子所产生的势 场的作用，这样在描述与电子势场时，必须考虑受作用的那一个电子状态和它本身作为其他电子运 动时对势场的贡献要一致起来，即自洽起来。这样的势场，成为自洽场。1、Hartree-Fock 方法固有近似性，它忽略了局部的电子-电子效应，即忽略了电子相关。因此 常采取组态相互作用 弥补这一缺点。组态指电子在其可能达到的各轨道上的排布方式的描述。通过 用不同

20、组态的波函数进行混合来表示分子的波函数称为 组态相互作用。利用组态相互作用给出基态 与激发态的能量可以预测电子的吸收频率、紫外可见光谱，用 研究键的断裂。2、另一个考虑电子相关能的方法是 MIler-Plesset 微扰法（MP）（MP2 最常用）。RHF 的组态相互作用而不是用 UHF 3、密度泛函（DFT）方法计算电子相关能。传统的 ab initio HF SCF 方法难以考虑电子相关作 用，而MP 等方法难以处理即使是稍稍较大的体系。与波函数方法相比。DFT 使得量子力学方法可 以直接用于大分子的计算，比如生物大分子。二、从头计算法 从头计算法（ab initio）是解全电子体系的非相

21、对论的量子力学方程，全部严格计算分子积分，不做任何近似处理。体系的 Hamilton 算符显全部电子（包括内层电子）的贡献。求解 Hartree-Fock 方程是一个试探和迭代的过程，全部计算达到进一步迭代时轨道能量（或总能量）不再变化（超过 某阈值）为止，于是就称这些轨道与它们所产生的位能场自洽，并将全部过程称为自洽场（SCF）方法。Roothaan 把分子轨道（M0）表示为原子轨道（AO）的基函数线性组合（LCAO-MO）,即用 LCAO 逼近Hartree-Fock 轨道，导出一组代数函数，即 Hartree-Fock-Roothaan方程或 Roothaan 方程。原子轨道的集合称为

22、基组（basis set）。极小基组 是 Slater 型轨道（Slater type orbitals,STO）,每 个占据轨道只用一个指数相表示，形式为 STO-nG,n 表示每个原子轨道的 Gaussian 函数个数，适当 表示 Slater 轨道至少需要有 3 个 Gaussian 函数。用具有不同指数来表示每个占据轨道，用两个指数 时的基组称为双E基组（由一个收缩指数和一个发散指数线性组合给出计算结果）。若对内层电子采 用一个指数，而对价电子采用两个指数，即为双E混合基组（一般表示为 3-21G,3 代表有 3 个 Gauss 函数描述内层轨道，价电子也有 3 个 gauss 函数，

23、2 个收缩函数，1 个发散函数）。在双E基组再加上 极化函数，用*号表示增加一组 d 型函数描述 2p 轨道分子中的极化情况，用*表示表示除增加 d 型 函数外，又增加 P 型函数表示 1s 轨道在分子中的计划情况，称为 双E扩展基组。从头计算法只是说明它是全电子的、非经验的计算方法。量子化学从头计算方法可以获得相当 高的精度，甚至达到所谓的化学精度。三、简化的从头计算方法 1、价电子从头计算方法（VE-AB INITIO）2、浮动球高斯轨道方法（FSGO）3、分子碎片法（MF）通常用于研究有机大分子。4、模拟从头计算分子轨道方法（SAMO）5、基于固体物理的 Slater-Xa方法 四、半经

24、验方法 定义：引入可调参数，体系 Hamilton 算符只显价电子的贡献。一种基于 Hartree-Fock-Roothaan 方程，借用经验或半经验参数代替分子积分。可分为处理 n价电子和处理全部价电子两大类。每一 类中又分为单电子法和双电子法两种。前者忽略了电子间的相互作用，后者则考察了电子间的排斥 能。（1）n价电子处理（Huckel 近似）：只考察n价电子，对于共轭分子是十分成功的。其单电子 法为简单分子轨道法（HMO），双电子法为 PPP 法。（2）全价电子处理：对于非平面型分子，已不能将 n电子单独分离出来，必须将 b电子包括在 内的所有价电子都考察进去。单电子法：包括推广的 HM

25、O 法（EHMO）。EHMO 从是否要求电荷自恰与自洽的严格程度又 分为简单 EHMO,电荷迭代 EHMO（IT-EMHO）和电荷迭代与组态相关 EHMO（C-EHMO,即 MWH 法）。EHMO 不能预测键长，计算电荷密度较高，一定条件下计算键角比较准确。由于 EHMO 方法 计算量小，又能在研究同系物和分子内部电荷分布、键性质、分子轨道能级及其对称性，较多应用 于大分子和含重元素的化合物或原子簇。双电子法：包括各种忽略微分重叠方法（NDO）及修正的 NDO 方法。NDO 根据双电子积分 时忽略程度又分为全忽略微分重叠方法（CNDO），间略忽略微分重叠方法（INDO）和忽略双原子 微分重叠方

26、法（NDDO）。半经验方法中常用的计算方法 CNDO、INDO:CNDO 是最简单的自洽方法，用于计算敞开和封闭体系基态的电子性质，是计 算药物电荷分布的首选方法（注：但并非精度最高的方法）。MINDO3 是将 INDO 许多相互作用的 计算用参数代替，计算分子基态性质较为成功，主要用于有机大分子，特别适用于含硫化合物。ZINDO/1 将 INDO 扩展到过渡金属，用于计算含过渡金属分子的能量与几何优化，ZINDO/S 可用于预测紫外可见光谱，不适用于计算几何优化和分子动力学。基于忽略双原子重叠方法有 MNDO 和 AMI 法，MNDO 方法对 NDDO 方法所做的的修正主要 是采用电荷多极矩

27、作用表示双中心电子积分，其次是引入带参数的函数来表示几个重要积分，最后,它用光谱数据对单中心双电子积分参数化时部份地考虑了相关能。AM1 是对 MNDO 方法的改进，也是最精确的方法之一，用于含有第一周期和第二周期元素的 有机分子，不适用于过渡金属，计算同时含有氮和氧的分子结果好于 MNDO。PM3 与 AM1 方法一 样，只是参数和不一样，其对非键相互作用计算优于 AM1，主要用于有机分子和主族元素。五、药物设计中量子力学活性指标 无论是药物分子还是生物大分子，都是具有一定结构的化学物质，它们的化学性质都是有其外 周电子的基本结构特性所决定的。1、轨道能量以前沿轨道能量最为重要，它们是 HO

28、MO（最高占有轨道能）和 LUMO（最低空 轨到能）。HOMO 能可作为分子给电子能力的量度，而 LUMO 能则可作为分子接收电子能力的量度。EHOMO与分子的电离势相关，作为分子给电子能力的量度，其值越小，该轨道中的电子越稳定，分 子给电子能力越小，ELUMO与分子的电子亲和能直接相关，其值越小，电子进入该轨道后体系能量 降低越多，该分子接收电子的能力越强。EHOMO与ELUMO之差是非常重要的稳定性指标，其差越大 稳定性越好，在化学反应中的活性越差，该值近似等于分子的最低激发能，所以有时也将其作为化 学反应的最低活化能（但这种想法忽略了分子在激发态时的重组，常出现错误）。2、电荷密度 的大

29、小可以反映各原子发生反应的倾向性，电子密度越大的位置与亲电试剂的反应 性越大，而电子密度越小的位置则于亲和试剂的反应性越大。3、键级（bond order，Prs）即键的数目，表示两个相邻原子间成键的强度。4、离域能亦称共轭能，是指通过电子的共轭作用发生超共轭而使体系趋于稳定，活化能降低的 能量。离域能越大，反应则越易于进行。超离域度（Sr是以轨道能级系数的倒数加权的电子密度，是离域能的度量参数。5、原子自极化率（nr）其值越大则表明他在攻击物质面前越容易调节自己的电荷，因而活性也 较大，常成为反应活性点。6、前沿电子密度 fr也称福井函数，前沿轨道的电子密度 fr是比较分子内部不同位置的化学

30、反应 活性的量度。注：前沿轨道电荷只能用来描述同一分子内不同原子的活性。7、静电势指一个单位正电荷从无限远处带到某一点所要做的功，比单个原子的静电荷密度更反 映实际。QM/MM 量子力学与分子力学相结合 在研究受体与药物相互作用时，有时要涉及到化学键的生成和断裂，近年来发展了量子力学与 经验力学相结合的方法 一一 QM/MM 模型。方法中最重要的一步是根据研究的目的划分体系，通常 将体系化分为如下几个区域：量子力学区域，分子力学区域和边界区域。第4章二维定量构效关系 定量构效关系是研究一组化合物的活性或 毒性或药代性质与其结构之间、物理化学性质之间之 间或者拓扑结构 之间的相关关系，用数理统计

31、和数理模型加以表征的研究方法。其意义为：在 受体 结构未知 的情况下，揭示化合物的结构与活性的依赖关系，建立表征这种关系的数学模型，以 新化合物的活性，演绎受体与药物结合并呈活性时的立体环境及物理化学要求。QSAR 数据主要由两部分组成：活性数据又可称为应变量，在 QSAR 中，应变量活性参数通常以产生标准生物效应时药物的物 质的量剂量八、预测 或物质的量浓度的负对数（Iog1/C）表示。化学描述变量（即为自变量）包括：（1）理化参数描述符 疏水性参数 LogP，脂水分配系数。通过摇瓶法和 HPLC 法测定。电性参数 Hammett?常数：表示芳香取代基的诱导和共轭效应之和；？值：取代基诱导效

32、应 之和；偶极矩（?？；解离常数 PKa;红外、紫外、NMR、MS 等光谱数据。立体参数 Taft 立体参数：Es=lg?（kx/kH）A Es 越小，表示体积越大，水解速度越慢。摩尔折 射：？?2?2-1?2+2?/?数值越大代表体积越大；Van der Waals 体积；STERIMOL 多 维立体参数。（2）分子片段描述符 描述符将分子中某一特征片断，如原子片断、环片断以及亚结构片断作为 描述符代码，是一种 拓扑学范畴的描述符。（3）分子连接性指数反映了分子中各原子排列状况、分支大小且与多种理化常数及生物活性相 关。可用质谱的质荷比 作为描述符。（4）其他参数 QSAR 的分析方法：线性

33、回归分析法和模式判别法。回归分析：是对一组数据进行最小二乘拟合处理并建立函数关系的过程。当有几种性质可能对 活性有贡献时，可用多元回归来处理。拟合函数的统计评价也是这种分析的一部分。常用的包括 Free-Wilson 方法和 Hansch 分析法等。线性关系的判别：因变量与自变量之间是否存在线性关系，可通过对回归系数进行 t 检验来断 定。进行 t 检验后，如果发现有的偏回归系数不显著，那么就要从回归方程中删除这些没有明显作 用的自变量，删除变量时不可同时将几个不显著的自变量一起去掉，应当先删去 t 值最小的一个变 量，重新计算回归方程，对新的回归方程的回归系数作检验，再删去新方程中一个不显著

34、的自变量，如此重复，直到回归方程中所有的自变量都显著为止。方程的显著性检验可用 复相关系数取值 R 断定。复相关系数 R 反映了因变量与所有自变量之间 回归关系密切的程度，检验 R 的显著性就是检验回归方程的显著性。因变量与自变量总体相关并不 意味因变量与每个自变量都显著相关，因此还需按上述方法对偏回归系数作显著性检验。虽然 R 是衡量总回归效果的重要标志，但 R 值的大小与回归方程中自变量的个数 n,及因变量 所取的观察值的个数 k 有关。当 n 相对于 k 不很大时，会获得较大的 R 值，即容易产生偶然相关。进行多元回归时要注意 n 与 k 的比例。一般认为 n 至少是 k 的 4 倍或

35、5 倍以上，就是说 1 个自变量 要求因变量有 45 个观察值对应，以消除偶然相关的影响。另外，需要注意，由自变量预测因变量 时，所用的数据不应超出建立回归方程时数据范围，不能随意外推。数据拟合后所得方程的好坏也可用 R2和 S 两种统计量来判断。R2是方程方差在数据方差中所 占的份数，R2=1 表示数据对方程完全适合，而 R2=0.50 表示数据中只有 50%的方差可用方程解释；补充统计量 S 是观察值与方程预测值的标准偏差。如果方程的 S 值较实际测定的标准偏差小，表明 数据拟合得较理想。回归分析通常被运用于结构非常近似的化合物。当化合物结构差别较大时，由于描述符类型过 多且难以产生相关性

36、，此时有 模式识别：以建立能区分活性种类的判别函数为目的。常用的方法：（1）聚类分析 预先将为数众多的不同取代基（称为样本或样品），按其化学结构参数（称作变量 或指标）进行分类，使化学结构参数相近者归为一类，不同者归为其它类别。对样本分类的方法为 Q 型聚类方法，对变量分类的方法为 R 型聚类分析法。（2）主成分分析合理地从 k 个主成分中挑选出少数几个彼此无相关性的主成分作代表，就可以 获取由原始变量提供的绝大部分信息。（3）非线性变换（4）因子分析 和主成分分析相当类似，利 用相关系数矩阵以少数几个互不相关的主因子来代表原始变量所提供的信息。（5）人工神经网络 QSAR 应用：（1）预测同

37、源物质的生物活性（2）避免合成过多的化合物；（3）更有目的地提高 化合物的选择性作用；（4）预测化合物的某些成药性质；（5）帮助了解药物的作用机制及描述受体 的图像。2D-QSAR 限制：使用前提（1）假定化合物的结构和生物活性之间存在一定的关系，也就是说，结构和活性之间存在函数关系（2）根据已知化合物结构一活性数据建立的函数，可以外推置新的化 合物（3）化合物的结构可用适当的结构描述符来表示（4）所有化合物的限速反应均相同（5）化合 物与酶或受体结合时，药物、酶或受体所引起的构象变化不予考虑（6）与受体亲和力有关的线性自 由能参数的具有加和性（7）药物在体内代谢的差异均予以忽略。定量构效关系

38、的应用注意以下几个方面：（1）只能应用于作用机制相同的化合物，作用机制不 同的化合物难于应用。一般认为，结构相近的同源物，其在体内的作用机制是相同的；（2）只能预 测同源物的生物活性，对于非同源物由于作用机制不同不能预测；（3）在预测生物活性方面并不都 是成功的。因为所得到的定量关系式还不能完全解释化合物与受体或酶间的作用情况，同时现有的 理化参数还不足以真正描述化合物生物活性的本质；（4）一组同源物质的生物活性变化幅度若小于 一个对数单位时，往往难于得到满意的相关结果，这可以由相关系数 与回归分析的化合物数目与所得到的相关式中参数项数目之比至少为 若有偏离较大的化合物需要除掉时，除掉的数据点

39、不能多于原化合物总数的 的因素而导致机会相关。对于被除去的数据点应给以合理的解释；（6）定量构效关系的研究不能代 替药物设计的所用工作，也不能发现结构全新的先导化合物，对于一个全新的化合物也无法预测他 的生物活性，这一方面可以通过三维定量构效关系进行。非极性化合物与水混合时会形成互不相溶的两相，即非极性分子有离开水相进入非极性相的趋 势，即疏水性，非极性溶质与水溶剂的相互作用则称为 疏水效应。疏水常数的定义：化合物在脂水两相中平衡浓度之比，通常将分子在 水-正辛醇体系中的分配系 数 logP 作为疏水性的度量，选择水-正辛醇体系是由于该体系 与生物体系相似。LogP 的测定方法：（1）摇瓶法。

40、摇动一个装有两种互不相溶溶剂和一种溶质的烧瓶，待达到平衡后，分析溶质在 其中一相或两相中的浓度。缺点：速度慢费时、繁琐、易受溶质的稳定性或纯度的影响。（2）薄层色谱法（原理是 log P 与色谱保留指数 RM 具有线性关系）；（3）反相高效液相色谱法；（4）估算方法碎片加合性法（原理是把分子划分为基本片段，每种特定的基本片段具有特 定的贡献值，整个分子的 logP 值是其所含的所有片段贡献的加和。由于该法假定分子的疏水常数具 有加和性，因而若偏离该假定则就有较大误差。）基于分子性质的计算方法（注意：含羧基等极性 基团的分子在疏水过程中包含电离、水合和形成离子对等现象，其 logP 的计算偏差较

41、大。）用 20 个已知 LOGP 的甾体化合物库名为建立只有 5 个参数的回归方程，预测化合物 LOGP 的 取舍的原则（有可能以实验题考核）：（1）若交叉相关系数 0.9,说明两参数是高度相关的，也就是说回归方程中只需其中一个就 可（2）若删除与 logP 相关系数W 0.05 的参数导致 F 值下降，则不应删除（3）对两个参数其它条件 相似可删除与 logP相关性小的参数（4）删除过程中应注意复相关系数 R2 和 Fisher 值，以判断参数 的取舍。Hansh 回归方程（即多元线性回归方程，又称二维定量构效关系）提出了 同源物的生物活性与 各种取代基的理化参数之间的依赖关系，采用最小二乘

42、法经多重回归，用与自由能相关的参数方程 表示也就是 线性自由能相关方法。Log1/C=-kin 2+k2n+k3+k4Es+ks,其中 C 为化合物产生指定 生物效应的物质的量浓度、n为疏水参数、b为电性参数、Es 为立体参数、k1,k2,k3,k4 代表权重即 r 的计算方法造成的；（5）参 5:1，在进行回归分析计算中，10%，以避免掺入人为 各因素贡献大小的系数。1、疏水参数n是被取代的化合物的脂水分配系数（正辛醇-水）（logP）与未取代的母体化合物 的脂水分配系数之差。-n 2表示药物的亲脂性和生物活性呈倒抛物线关系。随着脂溶性增加，到达 作用部位的概率也逐渐增加，达高峰后，又由于亲

43、脂性过大，则药物难以在水相中运转，而在脂肪 组织中含量增加，使到达作用部位的概率又下降，所以药物需要适当的疏水性，活性最大时的疏水 性为最适n值（logP0）。（参见下面的 HANSCH 抛物线模型）2、电性参数指取代基的电性效应对分子反应性能的影响。Hammett 取代常数（）反映芳环间 位或对位上侧链取代基的影响（共轭效应与诱导效应之和），b为正值表示为吸电子基；负值表示为推 电子基。Taft 取代常数（b*是反映脂肪族化合物中未构成共轭体系的取代基团极性效应的量度。3、立体参数：反映所代表基团的大小及它对配体-受体位相互作用接近 的影响。Taft Es 值为 用酸性介质中的水解速率常数反

44、映立体因素的影响。氢的 Es 值为零，基团越大，Es值越负。摩尔 折射（MR）反应是液体的物质的量体积，若方程中 MR 项为负值，即表示空间位阻为重要影响因 素，如果为正值，说明由诱导极化产生的色散力为主要影响因素。van der Waals 参数，用 van der Waals 体积（Vs）和半径（rv）来表示基团的实际大小，van der Waals 半径和 Taft Es值用作立体参数 的主要缺点是具有这些常数的基团数目很有限。最小立体差异（MSD）参数：反映一个化学式与另 一个化学式相互重叠在一起时，以其不重合部分的差异为基础。Kubinyi 双线性模型：该模型与 Hansch 的抛物

45、线模型的不同点是曲线的上升及下降部分均为直 线，仅在最适 logP0 附近为抛物线，且两直线的斜率根据数据的分布情况而变化，而不是像抛物线 那样上升和下降两部分为对称的弧线，因此对数据可更好地拟合。HANSCH-FUJITA 法的注意：（1）所有化合物必须是同源物，即母环相同（2）与受体具有相 同的作用机理（3）需要化合物具有较大的物理化学性质的差异（4）所选择的参数之间不能有相关 性（5）化合物的活性差异至少要大于 10 倍（6）建模化合物的数量至少是回归分析选用参数的 HANSCH-FUJITA 法的限制：（1）只考虑化合物与受体的作用位点，而不考虑化合物化合物与 受体结合时的构象变化（2

46、）所有参数只能表达二维意义上的结构特征，不能表征三维特征；不能定 量的解释三维结构与生物活性之间的关系（3）只能用于先导物的优化，不能用于先导物的发现。第五章基于受体的药物设计（直接设计法）分子对接（Molecular Docking），是将小分子化合物（配体）识别、适配和结合于生物大分子（受 体）的虚拟操作，旨在发现与受体结合部位相匹配的小分子化合物。是从虚拟库中筛选苗头和先导 化合物的重要方法。受体和药物分子之间通过空间匹配和能量匹配而相互识别形成分子复合物，并 预测复合物结构的操作过程。分子对接的用途：（1）研究分子间的结合模式（2）预测分子间的结合能力（3）药物设计（先导 化合物的发现

47、或者改造）（4）蛋白质工程 分子对接的原理：药物与受体分子的结合强度取决于结合的自由能变化（G结合），结合的自由 能变化与药物-受体复合物的解离常数的对数成线性关系。G结合=-2.303 RT lgKi，即 G结合=-5.85 lgKi 基于该公式，由分子对接的能量值即可预测配体与受体间的亲和力。5 倍。结合自由能的变化与系统的焓和熵的关系为：G结合=H结合-T S结合在对接的过程中，通常忽 略熵效应的变化，而焓变通常只考虑配体与受体结合部位的相互作用能（E相互作用）。E相互作用包括氢 键作用（EHB）,静电作用（Eelec）和范德华作用（Evdw）。在一定意义上，氢键和静电作用表征了配 体-

48、受体结合部位的化学互补性，而范德华作用表征了相互间形状互补的疏水相互作用。分子对接的一般过程：（1）确定靶标，找出结合位点，定出表面；（2）制备并调整受体位点以 及配体的构象；（3）生成配体能量较低的多构象；（4）计算对接时所有构象-受体相互作用能量；（5）进行分子动力学模拟；（6）复合物的全局最优结合构象。分子对接的分类：（1）刚性对接：研究体系的构象不发生变化（精度不高，但速度快，适合于 针对大容量数据库的快速虚拟高通量筛选。）；（2）半柔性对接（适合于阐述配体-受体结合模式，以 及针对小容量数据库的虚拟筛选。）（3）柔性对接（度很高，速度慢，适合于揭示配体-受体间的动 态结合过程以及深入

49、的作用机制分析，不能用于高通量筛选。）分子对接的核心步骤：（1）如何找到最佳的结合位置：配体的可能构象；配体在活性位点的可能取向 受体与配体的结合是在三维空间中实现的，受体的结合部位和配体的分子构象决定了分子对接 的结果。受体结构由于结构复杂，计算冗长，一般不计其构象变化。而针对配体分子，需要进行构 象搜索，确定其低能态构象群，这是对接的最关键的步骤。分子对接中构象搜寻的主要方法：蒙特卡罗法，遗传算法，模拟退火法，系统搜索和片段生长 法。分子对接方法的分类，主要是依据构象搜寻的方法进行的。（2）如何评价对接分子之间的结合强度：非键作用能基于分子表面的溶剂化计算半经验 的自由能计算 可以利用打分

50、函数 第六章基于配体的药物设计（间接设计法）计算机辅助药物设计（Computer Aided Drug Design,CADD）,是药物化学的一个分支，是以计 算机作为操作界面和辅助手段，利用计算化学、分子图形学、统计学和数据库等技术，研究药物和 受体的相互作用，发现和设计新的生物活性分子的方法，在先导化合物的发现和优化合成中起重要 作用。根据受体结构是否已知，可将 CADD 分为（1）基于受体的药物设计（直接设计法）包含分 子对接同源模建分子动力学碎片分子设计（2）基于配体的药物设计（间接设计法）包含药 效团定量构效关系数据库构建 基于配体的药物设计：在不知道受体三维结构的情况下，从研究一系