7.2 基于结构的虚拟筛选:分子对接 本节摘要:分子对接把配体放进靶点结合口袋,搜索姿态并用打分函数估算结合强度,是基于结构虚拟筛选(SBVS)的核心动作。本节走完"受体准备、配体准备、对接搜索、打分重排、位姿审查"的全流程,标定打分函数的误差量级与人工审查的看点,并用"一万分子对接换三成命中率"的案例给出完整的战役量感。 上一节靠已有活性分子间接推理,本节直接研究犯罪现场:靶点口袋长什么样、什么形状的分子能进去、进去后能赚多少相互作用能量。2.3 节的构象手艺在这里兑现第一价值。 有了靶点照片之后 受体准备。
本节摘要:分子对接把配体放进靶点结合口袋,搜索姿态并用打分函数估算结合强度,是基于结构虚拟筛选(SBVS)的核心动作。本节走完"受体准备、配体准备、对接搜索、打分重排、位姿审查"的全流程,标定打分函数的误差量级与人工审查的看点,并用"一万分子对接换三成命中率"的案例给出完整的战役量感。
上一节靠已有活性分子间接推理,本节直接研究犯罪现场:靶点口袋长什么样、什么形状的分子能进去、进去后能赚多少相互作用能量。2.3 节的构象手艺在这里兑现第一价值。
受体准备。对接前的靶点处理有四件事:选结构(解析质量、分辨率、配体共晶信息)、补氢与质子化状态(组氨酸质子化、天冬酰胺侧链翻转这类细节能改变口袋氢键网络)、处理水分子(保守水介导氢键要留下,松散水清走)、定义口袋(共晶配体周围、已知位点残基或全表面盲扫)。配体准备则是 2.3 节流程的复用:三维构象系综、质子化状态、部分电荷。
对接搜索。算法在口袋里给配体找姿态:位置、取向、可旋转键的二面角一起变。工程实现各异——增量构建、遗传算法、蒙特卡洛——但目标一致:在打分函数的地形上找低点,同时保持化学合理性(键长键角不许变形)。打分函数估算姿态的结合自由能,三条技术路线:经验打分(对氢键、疏水接触等项加权拟合实验数据)、知识库打分(从统计对接触计数)、力场打分(物理能量项)。它们的共同软肋是熵与水的处理——结合熵、水分子置换这些大头,简化打分函数只能间接近似。
误差量级要背下来:主流打分函数的相关系数普遍只有零点几,能量排序误差按千卡计,预测结合能差一个数量级是常态。所以对接打分的正确用途是粗排序,不是绝对预测;虚拟筛选的落地姿势是"对接粗排、共识重排、人工审查、实验验证"四连。

背景。某靶点拿到两埃出头的共晶结构(口袋里有片段配体),但没有可用的活性分子集——LBVS 缺弹药。团队决定对一万五千个类药采购分子做对接筛选,实验预算同样只够验证二十个。
操作。流水线五段全跑,两个细节决定成败。配体端用 2.3 节流程生成构象(每分子上限五十个、RMSD 阈值零点五埃),质子化状态按口袋环境 pH 估算;对接后不做"唯分数论",而是三打分共识(两种经验打分加一种力场重打分),取三者排名的几何平均,再由两位化学家独立审查前六十名的位姿:口袋里已知的锚点相互作用(共晶配体的关键氢键)是否复现、疏水骨架是否填充主要空腔、有没有与保守水网络打架。
结果。二十个送验分子里三个显示微摩尔级活性,命中率约百分之十五;加上两个弱活性的,可进入后续优化的系列有两套不同骨架。作为对照,未做位姿审查、纯按单打分函数取前二十的平行名单,只有一条命中——一倍以上的差距,全在审查工序里。
解读。三个可迁移的判断。其一,位姿审查的价值高于打分微调:打分函数之间相差的是噪声,人工审查筛掉的是"分数高但相互作用图景不合理"的假阳性(典型案例:配体悬在口袋口,靠一个边缘接触蹭到高分)。其二,命中率三成到一成五的量感来自靶点质量:口袋清晰、有共晶锚点的靶点 SBVS 表现好;口袋柔性大或结构预测分辨率低的靶点,先把结构功课补足再筛。其三,对接与 LBVS 并不互斥:本例若有几个活性分子,完全可以对接打分与指纹相似两路并联取并集,命中率与骨架多样性通常都更好。还有一条常被跳过的功课值得补记:受体准备要留版本——选用的结构、补氢方案、水分子处理、口袋定义各留一笔记录(10.2 节的留痕条款),因为对接结果对受体准备的敏感性极强,同一配体库换一版受体准备重跑,排名抖动两三成是常态;没有版本记录的对接结果,复盘时无从对质。
变式。沿着结构路线向上走:对对接命中的骨架做分子动力学重排与自由能微扰(FEP 类方法),把千卡级误差压到零点几卡级,代价是每个分子成百上千倍算力——只配用在优化阶段的三五极候选上;口袋柔性大时用诱导契合对接(残基侧链允许移动);共价靶点用共价对接处理反应项;RNA 与蛋白-蛋白界面则是另一套打分体系——同一个"搜索加打分"的骨架,参数与物理量全换。
位姿审查要练的就是认人——高频假阳性各有相貌,认得多了扫一眼就有数。
其一,口袋口的"挂在门把手上"。配体只把一个小基团探进口袋边缘,主体悬在溶剂区,靠单点接触蹭到不错的分数。判别动作:看配体质心与口袋中心的偏移、算埋入溶剂可及表面的比例,悬空过半的直接出局。其二,水网络打架。配体的极性基团顶进了必须留存的保守水分子位置,而打分函数没给这个冲突记账。判别动作:叠加共晶水分子坐标,凡与保守水近距重叠又不形成氢键替代关系的,重点怀疑。其三,质子化状态穿错衣服。配体制备时按溶液酸碱度定的质子态,可能与口袋微环境不符——组氨酸侧链与碱性胺最常中招,静电项随之整体失真。判别动作:对命中名单里的酸碱基团逐一核对口袋内可形成氢键的残基取向。其四,互变异构选错形态。酰胺与亚氨醇、酮式与烯醇式,选错形态等于让配体换了一张脸进考场。判别动作:对含互变体系的核心骨架,两种形态各对接一次,取物理上更可信的姿态而非分高的那个。其五,构象畸变凑分。柔性配体为了多挣一个相互作用,把自己扭成了高应变姿态,strain energy 一算反而倒赔。判别动作:抽查命中姿态与低能构象的重原子均方根偏差,偏差离谱者丢弃——2.3 节的构象系综在此就是裁判基准。
五相相貌的共同防御是把审查动作写成检查单:埋入比例、水冲突、质子态、互变形态、构象应变,五项各打勾或标注,命中名单逐条过。审查从"看看顺不顺眼"升级为可复核的工序,这也是 10.2 节"淘汰留痕"条款在对接场景的落点——每个出局者都带着哪一项不合格的记录。