2. 分子对接原理：锁钥模型与诱导契合
讲分子对接，绕不开两个经典模型——锁钥模型和诱导契合。很多人以为这只是理论课上的概念，跟实际操作关系不大。但我跟你说，理解清楚这两个东西，直接决定了你对接结果靠不靠谱。
我自己刚入行时，拿着一个刚性受体做了几百次对接，结果不管怎么调参数，对接姿势都怪怪的。后来才发现，问题出在——我把受体当成了“死”的，而它其实是“活”的。
2.1 锁钥模型：最朴素的匹配逻辑
1894年，Emil Fischer提出了锁钥模型。大意就是：受体是锁，配体是钥匙。钥匙形状对得上锁孔，才能插进去、转得动。
这个模型很直观，也支撑了早期的分子对接算法。在那个年代，计算资源金贵得不行，把受体当成刚性结构处理，是唯一可行的路。
- 受体刚性：对接过程中，受体原子坐标固定不动。
- 配体柔性：只允许配体分子发生构象变化。
- 打分依赖形状互补：对接分数主要由空间匹配度和静电作用决定。
但问题来了。真实情况哪有那么理想？蛋白质不是铁疙瘩，它是有弹性的。配体一过来，口袋周围的侧链甚至主链都可能挪位置。锁钥模型在这一点上，明显跟现实脱节。
2.2 诱导契合模型：更真实的结合过程
1958年，Daniel Koshland提出了诱导契合模型。核心思想很简单：受体和配体都是可变的。配体靠近时，受体的结合口袋会发生形变，两者相互适应，最终达到稳定的结合状态。
说白了，不是钥匙找锁，而是钥匙和锁孔一起变形，直到咬合。我当年做一个激酶项目，有个配体用刚性对接死活对不上，换成柔性对接后，发现口袋入口处一个酪氨酸侧链向外翻转了——那个姿势，正好让配体插进去。
- 初始接触：配体以较低亲和力靠近受体。
- 构象调整：受体侧链、loop区乃至整体结构发生位移。
- 稳定结合：双方找到能量最低的构象组合，形成复合物。
2.3 两种模型在AutoDock中的体现
AutoDock系列工具中，两种模型都有对应的实现方式。我们来拆开看。
| 模型 | AutoDock对应模式 | 速度 | 适用场景 |
|---|---|---|---|
| 锁钥模型 | 刚性对接（受体固定） | 快 | 大规模虚拟筛选、已知刚性口袋 |
| 诱导契合 | 柔性对接（指定侧链可旋转） | 慢 | 重要位点验证、难对接分子 |
在AutoDock 4中，如果你想做柔性对接，需要用AutoGrid准备柔性残基的格点图，再指定哪些侧链的二面角可以旋转。操作上比刚性对接多两步，但效果差异很大。
下面是我常用的柔性残基设置命令，写在一个配置文件中：
flexible_residues
residue_number 1
side_chain_atoms CB CG CD1 CD2
rotatable_bonds CB-CG CG-CD1 CG-CD2
2.4 从理论到实践：选择策略
你在实际项目中到底用哪个模型？我给一个简单判断流程：
- 如果你在跑预筛选：用锁钥模型（刚性对接）。库里有几十万分子，不可能每个都做柔性。
- 如果你拿到一个候选分子：换成诱导契合（柔性对接），确认结合姿势是否合理。
- 如果已知口袋有柔性：直接上柔性对接。我见过有些GPCR体系，口袋入口残基位移超过2Å，刚性对接直接报错。
- 如果打分结果跟实验对不上：检查是不是忽略了诱导契合效应。有时候调整一个侧链，打分能差好几倍。
2.5 两种模型的核心区别
再浓缩一下，你记住这张表就够了：
| 维度 | 锁钥模型 | 诱导契合 |
|---|---|---|
| 受体状态 | 完全刚性 | 部分柔性 |
| 配体自由度 | 构象可旋转 | 构象可旋转 |
| 计算量 | 低 | 高 |
| 真实度 | 偏低（忽略蛋白变形） | 较高（逼近真实结合） |
| 典型工具 | AutoDock刚性模式 | AutoDock柔性残基模式 |
嗯，说到底，没有哪个模型是万能的。我自己的做法是：把两种模型当成工具箱里的两把扳手，看情况换着用。遇到需要精细调参的项目，我甚至会在同一次对接中混用——先刚性跑全局，再柔性跑局部优化。
最后，给你一张图，直观感受两种模型在对接流程中的位置和区别。