分子对接实验室汇报

分子对接 定义 u虚拟筛选（virtual screening, VS) 也称分子对接，即在进行生物活性筛选之前，在计算机上对化合物分子进行预筛选，以降低实际筛选化合物数目，同时提高先导化合物发现效率。 化 合 物 数 据 库虚 拟 的 实 际 存 在 的优 势不 消 耗 样 品 ， 降 低 筛 选 成 本考 虑 化 合 物 分 子 的 药 动 学 性 质 和 毒 性 ， 增 加 筛 选 的 内 涵 虚拟筛选的效率 例：蛋白酪氨酸磷酸酯酶1B（PTP1B）抑制剂的发现 经虚拟筛选，再作生物学测试，虚拟筛选的命中率比随机的高通量筛选提高1,700倍 虚拟筛选技术的分类根 据 靶 点 的 结 构 知 识基 于 靶 点 结 构 的 虚 拟 筛 选基 于 配 体 相 似 性 的 虚 拟 筛 选 分 子 对 接药 效 基 团 搜 寻 一、基于靶点结构的虚拟筛选 分子对接 起源：受体-配体的锁和钥匙模型 配 体 受 体 复 合 物 对接 受体和配基之间通过能量匹配和空间匹配而相互识别形成分子复合物，并预测复合物结构的操作过程 基于分子对接的筛选方法 基于一个靶点（酶、受体、离子通道、核酸等）的三维结构，常常采用分子对接的虚拟筛选方法从小分子数据库中找到能与之匹配的候选化合物。所谓分子对接是基于两个或多个分子之间通过几何匹配和能量匹配相互识别的过程，即在药物分子和靶酶产生药效的反应过程中，两个分子先充分接近，然后采取合适的取向，使两者在必要的部位相互契合、发生作用，进而通过构象调整形成稳定的复合物。 在通常的分子对接中，小分子的构象是柔性的，如果筛选一个分子数目过多的数据库（如ACD-SC），那么整个虚拟筛选过程将非常耗时。因此，通常可先设定一些条件，比如Lipinski 的“5 倍律经验规则”等一系列类药性条件，先对该库进行过滤，从而快速缩小三维数据库的规模。另外，如果针对某个靶点，已经 获得相关抑制剂的结构，则可采用分子形状匹配的方法（如FlexS），对数据库进行初筛，保留其中与已知抑制剂形状相似的分子。经过这些初步筛选之后，再采用基于分子对接的虚拟筛选从数据库中找出可能与靶点相互匹配的有机小分子。 分子对接筛选常用的软件 DOCK 是应用比较广泛的对接软件之一，由Kuntz 等设计开发。它能自动模拟配体在受体活性位点的作用情况，并记录下最佳的相互作用方式。而且该软件能对配体的三维数据库进行搜索，因此被广泛用于基于受体结构的对接筛选。 在DOCK 中，活性位点的确定是通过软件包中的sphgen 程序来完成，它通过在受体表面所有的凹陷区形成负像，并对这些负像进行聚类分析，用户则从中挑选出所需要的一类作为活性区域的位置。在生成负像的基础上，就可以进行配体分子和受体活性口袋之间的匹配，配体分子也采用一组球集来表示，并且DOCK 进行对接时，配体可以是刚性的，也可以设定为柔性。 在对接结束后，DOCK 程序则采用自带的打分函数对配体-受体之间的匹配情况 进行评价，其中包括原子接触得分和能量得分。所谓接触原子是指在一定距离之内的原子（一般定义为4.5 ），如果配体和受体之间存在这个距离内的原子，则认为产生碰撞，作为罚分从总得分中扣除；能量得分主要来自配体和受体间非键相互作用能。 AUTODOCK 也是常用的分子对接软件包之一，由Scripps 的Olson 科研小组开发。它采用模拟退火和遗传算法寻找受体和配体最佳的结合位置，用半经验的结合自由能方法来评价两者之间的匹配情况。为了加快计算速度，AUTODOCK 采用了格点对接的方法，格点上保存的是探针原子和受体之间的相互作用能，包括了范德华相互作用能、静电作用能和氢键相互作用能等。在AUTODOCK3.0 以前，只能实现单个配体和受体分子之间的对接，但目前的版本（AUTODOCK4.0）已经能实现对数据库对接筛选的功能。 除此之外，还有一些常用的对接筛选程序，比如FlexX、GOLD、Affinity、Glide 等，它们都各自开发出一套相应的对接筛选策略、打分函数，使得对接筛选的应用越来越广泛。 基于结构虚拟筛选一般流程图 建 立 大 量 化 合 物 的 三 维 结 构 数 据 库将 库 中 的 分 子 逐 一 与 靶 标 分 子 进 行“ 对 接 ” （ docking）优 化 小 分 子 化 合 物 的 取 向 及 构 象 等找 到 小 分 子 化 合 物 与 靶 标 大 分 子 作用 的 最 佳 构 象计 算 其 相 互 作 用 及 结 合 能 完 成 所 有 分 子 的 对 接 计 算找 出 与 靶 标 分 子 结 合 的 最 佳 分 子 对 接 方 法 尚 需 解 决 的 问 题分 子 的 柔 性溶 剂 化 效 应打 分 函 数 (一) 靶点结构的预处理靶 点 结 构 的 预 处 理 1、 靶 点 的 检 验2、 靶 点 的 处 理3、 确 定 靶 点 的 配 体 结 合 口 袋 1、 靶 点 的 检 验晶 体 结 构 （ X射 线 测 定 ） 溶 液 结 构 （ NMR测 定 ）分 辨 率 R 因 子温 度 因 子 几 何 构 型 的 准 确 率 小 于 2.5 Rfree28%R25%重 点 关 注 部 位 不 大 于整 个 分 子 的 平 均 温 度 因 子 角 至 少 有 90%落 在 允 许 区 域 内 PDB多 个 构 象Profile-3D：三 维 结 构 和 氨 基 酸 序 列 相 容 性ProsaII：残 基 之 间 相 互 作 用 能 量 评 估拉 氏 图 ：检 测 三 维 结 构 立 体 构 型 好 坏 2、 靶 点 的 处 理补 齐 晶 体 结 构 中 缺 失 原 子 和 残 基为 大 分 子 加 上 氢 ， 并 分 配 相 应 电 荷确 定 带 电 残 基 的 质 子 化 状 态去 除 不 重 要 的 小 分 子 3、 确 定 靶 点 的 配 体 结 合 口 袋靶 点 分 子 为 复 合 物靶 点 分 子 不 是 复 合 物实 验 数 据 未 知 为 配 体 为 中 心 的 5-7区 域 内 的 氨 基 酸 残 基根 据 同 源 蛋 白 或定 点 变 突 的 数 据 确 定 利 用 软 件 搜 寻 小 分 子 数 据 库 的 来 源 商 用 化 合 物 数 据 库公 司 或 研 究 机 构 自 有 数 据 库设 计 的 虚 拟 化 合 物 库 常 用 的 数 据 库 有 MDL数 据 库 、 SPECS数 据 库 和CNPD（ Chinese Natural Product Database） 数 据 库(二) 小分子数据库的预处理 (三) 分子对接 概念： 将配体分子放置到受体大分子的活性位点中，预测小分子与受体结合构象及作用能的过程。是基于受体分子结构虚拟筛选的核心。 目的：从小分子数据库中发现合适的化合物作为受体大分子的配体。从整体上考虑配体与受体结合的效果。 关键环节：配体结合构象的优化系 统 搜 索 随 机 搜 索确 定 性 搜 索 片 段 生 长 法构 象 搜 索 法构 象 库 方 法 蒙 特 卡 罗 法模 拟 退 火遗 传 算 法禁 忌 搜 索 法 分 子 动 力 学 模 拟 目的： 评判配体分子和受体结合能力的强弱。 含义： 先对同一个分子的不同结合构象，评价各构象的结合好坏。 再对数据库中的不同分子的最好结合构象进行评价，以得到最终的结合能力从高到低的化合物分子清单。 (四) 打分函数 1、 基 于 力 场 的 打 分 函 数 2、 半 经 验 的 自 由 能 打 分 函 数 3、 基 于 知 识 的 打 分 函 数 4、 “ 一 致 性 ” 打 分 (五) 部分实例 人碳酸酐酶与其抑制剂复合物 2002 年Grunberg 等采用基于分子对接的虚拟筛选方法成功地找到了多种人碳酸酐酶（humancarbonic anhydrase）的抑制剂。在整个设计过程中，他们采用了多次初筛的办法将Maybridge 数据库（61 186 个分子）和LeadQuest 数据库（37 841 个分子）进行过滤。首先利用Lipinski 的“5 规则”将数据库缩小为5 904 个；然后利用FlexS 与已知抑制剂进行相似性筛选，得到了100 个候选化合物；最后将这100 个分子利用 FlexX 程序进行对接筛选，从中挑选出13 个进行生物活性测试，结果7 个分子的IC50 值达到了微摩尔级别。见图。 抗艾滋病药物的发现 虚拟筛选 艾 滋 病 病 毒 , 人 类 免 疫 缺 陷 病 毒human immunodeficiency virus, HIV 1.HIV-1蛋白酶（HIV-PR）, 所催化的水解反应在艾滋病病毒导入人体细胞过程中起着重要的作用2.高效的HIV-PR抑制剂为治疗艾滋病的有效药物3.肽类HIV-PR抑制剂生物性质不稳定，吸收性差，易被代谢分解，因此口服给药无效 1. X-射线晶体结构2. 搜寻数据库3. 生物测试: 高选择性, 高活性( Ki = 0.1 nM) 抗 SARS冠 状 病 毒 药 物 的 设 计基 于 SARS-CoV 3CL蛋 白 酶 的 虚 拟 筛 选严 重 急 性 呼 吸 道 综 合 征 病 原 体 SARS冠 状 病 毒 SARS-Cov 感 染 宿 主 细 胞 起 重 要 作 用 的 结 构 蛋 白 E蛋白（small envelope protein，小信封蛋白） S蛋白（spike glycoprotein，刺突糖蛋白） M蛋白（membrane protein，膜蛋白） N蛋白（nucleocapsid protein，核衣壳蛋白） 多聚酶（polymerase） 类3C蛋白酶（3C like proteinase, 3CL） 3CL蛋 白 酶 作 为 抗 SARS药 物 筛 选 靶 点 的 优 点 在冠状病毒复制过程中起着重要作用 有许多已知抑制剂，便于迅速开发 较易表达，有利于加紧研究 有较高的同源性，可用同源法模建三维结构模型 步骤1. 同源模建(1) 3CL蛋白酶序列(GenBank)与各类冠状病毒蛋白酶序列（PDB）作序列分析和同源性分析（BLAST程序） 人冠状病毒；鼠科肝炎病毒；猪传染性腹泻病毒；猫传染性腹膜炎病毒；禽传染性支气管炎病毒；猪冠状病毒；传染性胃肠炎病毒(2) 传染性胃肠炎病毒（TGEV）的蛋白酶M pro与3CL蛋白酶有极高的同源性，特别在底物结合口袋（活性部位）(3) 以TGEV Mpro的X-射线晶体结构为模板，模建3CL蛋白酶三维结构（Sybyl 6.8 / SiteID程序） 结果（1）所建模型与TGEV Mpro 晶体结构基本重叠（2）3CL蛋白酶的折叠方式与TGEV Mpro相同，结合口袋的结构以及空间特征几乎一样 3CL蛋 白 酶 结 构 与 Mpro蛋 白 酶 晶体 结 构 的 重 叠 图 3CL蛋 白 酶 的 缎 带 模 型 步骤2. 分析酶-配体作用模型 两种蛋白酶的结合部位(Sybyl 6.8 / MOLCAD程序) 中,小分子C能以同样的方式与两种酶的结合口袋契合A. TGEV M PRo蛋 白 酶B. SARS 3CL蛋 白 酶C. 蛋 白 酶 抑 制 剂两 种 蛋 白 酶 的 底 物 结 合 口 袋 的 表 面 特 征 3CL蛋白酶模建模型或TGEV Mpro的晶体结构均可作为筛选抗SARS药物的结构模型 步骤3. 虚拟筛选以SARS冠状病毒3CL蛋白酶三维结构模型和TGEV Mpro 为筛选模型作虚拟筛选（SGI Origin 3800超级计算机和392CPU的神威1号超级计算机）ACD数据库、MDDR数据库、SPECS数据库、中国天然产物数据库（CNPD）和国家药物筛选中心内部样品库 共数十万个化合物（1） DOCK 4.0作初筛， 选出得分高的前1000个化合物；（2）用Cscore软件和AutoDock 3.0软件作评价，从每个数据库中挑选出 100个得分最高的化合物结果：共找到300个可能具有抗SARS冠状病毒潜力的候选化合物 步骤4 药理测试（1）用病毒3CL蛋白酶分子水平筛选模型筛选候选化合物 发现了7个具有高活性的化合物（2）在P3实验室中作SARS病毒感染细胞水平的测试，发现5-HT受体拮抗剂（ 肉桂硫胺，Cinanserin）有明显的抗SARS病毒感染和保护细胞的作用（3）申请专利，以CADD作结构优化 例 美 普 他 酚 双 配 体 衍 生 物 与 AChE的 对 接 研 究 阿尔茨海默症（Alzheimers Disease，AD）与乙酰胆碱（ACh）水平降低和对乙酰胆碱酯酶AChE诱导的-淀粉样蛋白（A）聚集有关 治疗AD 的药靶： 乙酰胆碱酯酶（AChE）和-淀粉样蛋白 AChE抑制剂：美普他酚（Meptazinol） 研究发现A的聚集与AChE上的外周阴离子位点（PAS）有关 Xie Q et al. J. Med. Chem. 2008,51(7):2027 美 普 他 酚 双 配 体 类 似 物 的 设 计 在美普他酚分子中通过引入不同长度连接链，希望使另一个配体能与PAS相互作用合 成 n=2-12的 多 个 美 普 他 酚 双 配 体 类 似 物 ， 发 现 n=9的 AChE抑 制 活 性 最 强 （ IC50 3.9nM） , 比 美 普 他 酚 高 2万 倍 ， 同 时 对 AChE诱 导 的 A聚 集 具 有 明显 的 抑 制 作 用 （ IC50 79M） , 具 双 重 作 用 . 以 分 子 对 接 证 明 作 用 模 式 小鼠乙酰胆碱酯酶（mAChE） PDB数据库中的琥珀酰胆碱与mAChE的复合物（PDB编号：2HA2） 双配体分子（n=9）的三维结构 以CORINA软件生成，将N原子质子化，再用Tripos力场进行分子力学优化 以GOLD 3.0进行分子对接，最后用SYBYL /CScore一致性打分和半经验自由能评价函数XSCORE打分，评价出最优结合构象 对 接 模 型 ： 证 明 了 AChE抑 制 作 用 和 A 聚 集 抑 制 作 用 双配基配体与电鳐TcAChE复合物的单晶X衍射证明：结合方式与对接相似