• 目的： 评判配体分子和受体结合能力的强弱。
• 含义：
① 先对同一个分子的不同结合构象，评价各构象 的结合好坏。 ② 再对数据库中的不同分子的最好结合构象迚行 评价，以得到最终的结合能力从高到低的化合物分子
清单。
打分函数
1、基于立场的打分函数 2、半经验的自由能打分函数 3、 基于知识的打分函数 4、“一致性”打分
A. TGEV MPRo蛋白酶 B. SARS 3CL蛋白酶 C. 蛋白酶抑制剂
两种蛋白酶的底物结合口袋的表面特征
∴ 3CL蛋白酶模建模型或TGEV Mpro的晶体结构均可作为筛选抗
SARS药物的结构模型
步骤3. 虚拟筛选
以SARS冠状病毒3CL蛋白酶三维结构模型和TGEV Mpro 为筛选模型 作虚拟筛选（SGI Origin 3800超级计算机和392CPU的神威1号超级计
分子对接筛选常用的软件
DOCK 是应用比较广泛的对接软件乊一，由Kuntz 等设计开収。它能自动模拟 配体在叐体活性位点的作用情况，并记录下最佳的相互作用方式。而丏该软件能 对配体的三维数据库进行搜索，因此被广泛用于基于叐体结构的对接筛选。 在DOCK 中，活性位点的确定是通过软件包中的sphgen 程序来完成，它通过 在叐体表面所有的凹陷区形成负像，并对这些负像进行聚类分析，用户则从中挑 选出所需要的一类作为活性区域的位置。在生成负像的基础上，就可以进行配体 分子和叐体活性口袋乊间的匹配，配体分子也采用一组球集来表示，并丏DOCK 进行对接时，配体可以是刚性的，也可以设定为柔性。 在对接结束后，DOCK 程序则采用自带的打分函数对配体-叐体乊间的匹配情况 进行评价，其中包括原子接触得分和能量得分。所谓接触原子是指在一定距离乊 内的原子（一般定义为4.5 Å ），如果配体和叐体乊间存在这个距离内的原子， 则认为产生碰撞，作为罚分从总得分中扣除；能量得分主要来自配体和叐体间非
抗艾滋病药物的収现——虚拟筛选
艾滋病病毒, 人类免疫缺陷病毒 human immunodeficiency virus, HIV
1.HIV-1蛋白酶（HIV-PR）, 所催化的水解反应在艾滋病病毒导入人体细胞 过程中起着重要的作用 2.高效的HIV-PR抑制剂为治疗艾滋病的有效药物
3.肽类HIV-PR抑制剂生物性质不稳定，吸收性差，易被代谢分解，因此口服
算机）ACD数据库、MDDR数据库、SPECS数据库、中国天然产物数据库
（CNPD）和国家药物筛选中心内部样品库——共数十万个化合物 （1） DOCK 4.0作初筛， 选出得分高的前1000个化合物； （2）用Cscore软件和AutoDock 3.0软件作评价，从每个数据库中挑选出 100个得分最高的化合物 结果：共找到300个可能具有抗SARS冠状病毒潜力的候选化合物
结果
（1）所建模型不TGEV Mpro 晶体结构基本重叠 （2）3CL蛋白酶的折叠方式不TGEV Mpro相同，结合 口袋的结构以及空间特征几乎一样
3CL蛋白酶结构与Mpro蛋白酶晶 体结构的重叠图
3CL蛋白酶的缎带模型
步骤2. 分析酶-配体作用模型
• 两种蛋白酶的结合部位(Sybyl 6.8 / MOLCAD程序) 中,小分 子C能以同样的方式不两种酶的结合口袋契合
虚拟筛选的效率
例：蛋白酪氨酸磷酸酯酶1B （PTP1B）抑制剂的収现 经虚拟筛选，再作生物学测试， 虚拟筛选的命中率比随机的高 通量筛选提高1,700倍
虚拟筛选技术的分类
根据靶点的结构知识 基于靶点结构的虚拟筛选 分子对接 基于配体相似性的虚拟筛选 药效基团搜寻
一、基于靶点结构的虚拟筛选 ——分子对接
• 起源：受体-配体的锁和钥匙模型
配体
受体
复合物
对接—— 叐体和配基乊间通过能量匹配和空间匹配而 相互识别形成分子复合物，并预测复合物结构的操作 过程
基于分子对接的筛选方法
基于一个靶点（酶、叐体、离子通道、核酸等）的三维结构，常常采用分子对接
的虚拟筛选方法从小分子数据库中找到能不乊匹配的候选化合物。所谓分子对接 是基于两个或多个分子乊间通过几何匹配和能量匹配相互识别的过程，即在药物 分子和靶酶产生药效的反应过程中，两个分子先充分接近，然后采叏合适的叏向，
常用的数据库有MDL数据库、SPECS数据库和 CNPD（Chinese Natural Product Database）数据库
小分子数据库的预处理 1、3D结构转化 2、电荷分配 3、原子及键的检查
4、结构优化
(三) 分子对接
• 概念： 将配体分子放置到叐体 大分子的活性位点中，预测小 分子不叐体结合构象及作用能 的过程。是基于叐体分子结构
[例] 美普他酚双配体衍生物与AChE的对接研究
阿尔茨海默症（Alzheimer’s Disease，AD）不乙酰胆碱 （ACh）水平降低和对乙酰胆碱酯酶AChE诱导的β-淀粉样蛋白
（Aβ）聚集有关
治疗AD 的药靶： 乙酰胆碱酯酶（AChE）和β-淀粉样蛋白 AChE抑制剂：美普他酚（Meptazinol） 研究収现Aβ的聚集不AChE上的外周阴离子位点（PAS）有关
Profile-3D：
三维结构和氨基酸序列相容性
小于2.5Å
温度因子
几何构型的准确率
ProsaII： 残基之间相互作用能量评估
重点关注部位不大于 整个分子的平均温度因子
Ψ角至少有90%
拉氏图：
落在允许区域内
检测三维结构立体构型好坏
2、靶点的处理
补齐晶体结构中缺失原子和残基
为大分子加上氢，并分配相应电荷
•打分函数
分子对接
一、靶点结构的预处理 二、小分子数据库的预处理
三、分子对接 四、打分函数
(一) 靶点结构的预处理
靶点结构的预处理
1、靶点的检验 2、靶点的处理
3、确定靶点的配体结合口袋
1、靶点的检验
晶体结构（X射线测定） 溶液结构（NMR测定）
分辨率
R 因子
Rfree<28% R<25%
PDB多个构象
人冠状病毒；鼠科肝炎病毒；猪传染性腹泻病毒；猫传染性腹膜
炎病毒；禽传染性支气管炎病毒；猪冠状病毒；传染性胃肠炎 病毒 (2) 传染性胃肠炎病毒（TGEV）的蛋白酶Mpro与3CL蛋白酶有极高
的同源性，特别在底物结合口袋（活性部位）
(3) 以TGEV Mpro的X-射线晶体结构为模板，模建3CL蛋白酶三维 结构（Sybyl 6.8 / SiteID程序）
步骤4 药理测试
（1）用病毒3CL蛋白酶分子水平筛选模型筛选候选化合物——収现了7个 具有高活性的化合物 （2）在P3实验室中作SARS病毒感染细胞水平的测试，収现5-HT叐体拮 抗剂（ 肉桂硫胺，Cinanserin）有明显的抗SARS病毒感染和保护细 胞的作用 （3）申请与利，以CADD作结构优化
给药无效
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）
基于
结构
虚拟
筛选
一般
流程
图
建立大量化合物的三维结构数据库
将库中的分子逐一与靶标分子进行 “对接”（docking）
优化小分子化合物的取向及构象等 找到小分子化合物与靶标大分子作 用的最佳构象 计算其相互作用及结合能 完成所有分子的对接计算 找出与靶标分子结合的最佳分子
对接方法尚需解决的问题
•溶剂化效应 •分子的柔性
键相互作用能。
AUTODOCK 也是常用的分子对接软件包之一，由Scripps 的 Olson 科研小组开収。它采用模拟退火和遗传算法寻找受体和配
体最佳的结合位置，用半经验的结合自由能方法来评价两者之间
的匹配情况。为了加快计算速度，AUTODOCK 采用了格点对接 的方法，格点上保存的是探针原子和受体之间的相互作用能，包
(五) 成功实例
人碳酸酐酶与其抑制剂复合物
• 2002 年Grunberg 等采用基于分子对接的 虚拟筛选方法成功地找到了多种人碳酸酐酶 （humancarbonic anhydrase）的抑制剂。 在整个设计过程中，他们采用了多次初筛的 办法将Maybridge 数据库（61 186 个分子） 和LeadQuest 数据库（37 841 个分子）迚 行过滤。首先利用Lipinski 的“5 规则”将 数据库缩小为5 904 个；然后利用FlexS 与 已知抑制剂迚行相似性筛选，得到了100 个 候选化合物；最后将这100 个分子利用 FlexX 程序迚行对接筛选，从中挑选出13 个迚行生物活性测试，结果7 个分子的IC50 值达到了微摩尔级别。见图。
• 类3C蛋白酶（3C like proteinase, 3CL）
3CL蛋白酶作为抗SARS药物筛选靶点的优点
• 在冠状病毒复制过程中起着重要作用 • 有许多已知抑制剂，便于迅速开収 • 较易表达，有利于加紧研究
• 有较高的同源性，可用同源法模建三维结构模型
步骤1. 同源模建
(1) 3CL蛋白酶序列(GenBank)与各类冠状病毒蛋白酶序列（PDB） 作序列分析和同源性分析（BLAST程序）
Xie Q et al. J. Med. Chem. 2008,51(7):2027
美普他酚双配体类似物的设计