分子对接怎么做：靶点分型、软件匹配与结果验证的实战框架

分子对接怎么做——这个问题在开始任何软件操作之前，需要先回答两个更前置的问题：靶点是什么类型？配体的化学空间有多大？靶点的结构特性（水解酶vs激酶vs GPCR vs核酸结合蛋白）和配体的复杂性（类药小分子vs多肽vs天然产物）决定了后续对接方法的选择。跳过靶点分型直接跑软件，是新手对接翻车的头号模式。

靶点分型：对接策略因靶点而异

激酶（Kinase） 的ATP结合位点是一个深埋的裂隙状腔体，铰链区（hinge）固定的1-3条关键氢键是激酶抑制剂结合的核心motif。对接激酶的要点：盒子中心必须精确定位在铰链区主链的O/N原子处——如果盒子偏移>2 Å，采样空间绕过了氢键最佳距离，top pose很可能缺少铰链区氢键，是个无效pose。激酶对接的成功率比较高（~70-80% redocking RMSD<2 Å），因为ATP腔体形状互补天然引导配体采样。

蛋白酶（Protease） 的催化位点以催化三联体/二联体（Ser-His-Asp或Cys-His）为中心，底物结合在延伸的裂隙中。对接蛋白酶时的关注点是催化残基的质子化状态——催化的His（pKa可被环境抬升到~7-8，在中性pH下处于中性态）、催化的Ser/Cys（需要去质子化才能亲核攻击）——这些残基的质子化状态直接决定对接中能否形成正确的催化构象。Vina不会自动处理这种非标准质子化，预处理阶段手动调节是关键。

核受体（Nuclear Receptor） 的配体结合域（LBD）是一个大型疏水腔体（~600-1000 ų），配体以疏水堆积和少数极性接触的方式结合。大型腔体对接的主要风险是搜索空间过大——配体在腔内的位置可以大范围滑动，打分函数难以区分真实的结合位点与”贴着疏水墙的非特异性结合”。核受体对接的高成功率技巧：利用共晶配体的结构坐标确定盒子的精确位置和尺寸，并在对接后对pose做cluster一致性检查。

金属蛋白酶（Metalloprotease） 的锌离子配位是对接的硬骨头——标准AD4/Vina对Zn²⁺的vdW参数化不够（默认半径范德华赋予值可能让配体与Zn的距离偏大0.5-1.0 Å），需要手动调整金属离子的vdW半径或使用专门参数化的金属离子力场。如果对接pose中配体的供电子原子与Zn的距离>2.5 Å，几乎可以肯定不是正确的配位结合模式。

软件匹配：不是越新越好

AutoDock Vina是目前最广泛使用的开源对接工具——速度快（一分钟内完成一个标准配体的对接）、精度可接受（~70% redocking成功率）、可编程性强（Python接口smina/vina）。对于初始虚拟筛选和hit识别阶段，Vina是首选。

AutoDock 4（AD4）相比Vina的优势在于可定制性强——用户可以编辑原子类型参数、范德华半径和去溶剂化参数。对于非标准原子类型（含硼、含硅药物分子）、卤素键等Vina未充分参数化的相互作用，AD4允许手动干预参数表，比Vina更能精准调校。

商业软件（Glide SP/XP、GOLD、MOE Dock）的核心优势在于打分函数的校准更精细（基于更大的专有训练集）和诱导契合流程的自动化程度——Schrödinger的IFD-MD可以自动对活性位点侧链做MD采样+对接迭代，对诱导契合效应明显的靶点（如HIV蛋白酶、BACE1）的pose预测精度优于Vina。

单配体对接→批量虚拟筛选→结果分析的标准管线

一个完整的对接项目通常管线化：准备靶点→单配体redocking验证参数→库准备（配体数据库，如ZINC、Enamine REAL）→批量对接→得分过滤（如top 10% per target）→pose聚类（1% cluster cutoff）→目视检查。每一步的过滤标准都需要根据靶点特性设定。

Redocking验证是质控的起点：将共晶配体从蛋白结构中提取出来→对接回蛋白→比较对接pose和晶体pose的RMSD（通常以重原子RMSD<2.0 Å为成功阈值）。如果redocking的RMSD>3 Å→对接参数设置有问题，回到盒子定位、配体质子化状态、exhaustiveness等参数逐一排查。redocking不通过就不该开始真正的虚拟筛选——不然筛出来的top hits可能跟靶点的结合模式完全不一样。

跑完对接后，对top hits的化学多样性做二次筛选——结构冗余度大的top cluster不能都带进下一轮，取化学空间分散的top结构（基于ECFP4或MACCS指纹做Tanimoto距离过滤，取距离>0.4的独立簇）。

分子对接怎么做——从靶点分型开始，到软件匹配合适、到redocking通过、再到结果质控，每一步都是在最小化”垃圾进垃圾出”的风险。

更多内容请访问 https://www.keyanxueshu.com/