第4章：格点参数设置——Grid Box大小与位置、格点分辨率、格点类型
格点参数，说白了就是你在受体周围画一个“探测范围”。
做过AutoDock对接的朋友都知道，Grid Box设大了，算得慢还容易出假阳性；设小了，配体可能跑出框外。我记得第一次做虚拟筛选时，用的还是AutoDock 4.2，Box设得太小，跑了三天筛了几万个化合物，结果全是无效对接——配体压根没放进盒子里。
所以这一章，我们把格点参数的三个核心维度掰开揉碎讲清楚：位置与大小、分辨率、格点类型。你掌握了这些，后面的特征构建才有根基。
4.1 Grid Box大小与位置——你画圈，配体往里站
Grid Box的中心，通常指向受体结合位点。我习惯先用Pymol或Discovery Studio打开受体结构，手动挑出一两个关键残基——比如激酶铰链区的Met、Asp——把它们的坐标取个均值作为Box中心。
快速定位Box中心的命令（以AutoDock Tools为例）
# 在ADT中点击 Grid -> Grid Options
# 手动输入Center坐标（x,y,z）
# 或者选中关键残基后：
Grid -> Set Center -> From Atom SelectionBox大小怎么定？黄金法则：Box边长 = 配体最长维度 + 5～8 Å。
| 场景 | Box边长建议 | 为什么这么定 |
| 已知结合位点+小配体 | 15～18 Å | 让配体有足够空间旋转，但不过度探索无关区域 |
| 已知位点+中等配体（MW 300-500） | 20～25 Å | 标准虚拟筛选常用大小 |
| 盲对接（未知位点） | 全蛋白包裹（30～60 Å） | 覆盖整个受体，但计算量陡增 |
| 大分子复合物对接 | 40～60 Å | 必须包含全部界面残基 |
我个人比较习惯的做法是：先用较大的Box跑一轮粗筛，看看配体倾向结合在哪里。找到热点区域后，缩小Box，精细重跑。这样既省时间，结果又可信。
遇到过一次很麻烦的情况：某次做虚拟筛选，Box中心偏移了2 Å，结果筛出来的“高活性”化合物全是假结合。后来才发现，配体的柔性端在打分时“伸”到了一个不可能的构象里。那之后我每次都会检查Box是否完全包住了配体的所有可能位置。
⚠️ 易犯错误：Box设得太大
不要以为Box越大越保险。Box大了，空区域里的噪声打分反而会拉低真实结合构象。我自己的经验是：Box最好只比配体最长尺寸多6 Å左右。多一分冗余，少一分浪费。
4.2 格点分辨率——精细还是粗糙？看场景
分辨率，说白了就是把Grid Box切成多少个小立方体。默认值通常是0.375 Å（对应大约 1/2 碳碳键长）。
为什么是这个数字？嗯，这是经过大量测试的经验值。太粗了（比如0.5 Å以上），对配体微小位姿变化不敏感；太细了（0.2 Å以下），格点数量爆炸，算到天荒地老。
| 分辨率 | 格点数量（Box 20ų） | 适用场景 |
| 0.250 Å | ~512,000 | 最终精修、研究关键残基侧链构象 |
| 0.375 Å | ~152,000 | 标准对接、中等规模筛选 |
| 0.500 Å | ~64,000 | 大规模初筛、阳性化合物粗选 |
一个实用的参数组合（我自己常用的）
# 虚拟筛选首轮（粗筛）
resolution: 0.500
box_size: 22.5 (边长)
grid_types: electrostatic, vdw
# 命中化合物精修
resolution: 0.250
box_size: 15.0 (边长)
grid_types: electrostatic, vdw, hbond我曾经有一回偷懒，一个项目全用0.5 Å去跑，最后Top10的化合物做活性测试，九个是假的。后来用0.25 Å重跑，结果才稳定下来。所以说，分辨率不只是消耗资源，它直接决定了你的筛选命中率有多高。
不过话说回来，也不是越细越好。你想想看，如果你的晶体结构分辨率本身只有3.0 Å，那格点弄到0.2 Å就有点自欺欺人了——初始结构里就有0.3 Å以上的误差。
4.3 格点类型：静电、vdW、氢键——三足鼎立
AutoDock系列在计算格点时，默认会算三种势场：
- 静电势（Electrostatic）：基于库伦公式计算。需要受体原子带部分电荷。这部分对极性配体、带电配体非常敏感。
- 范德华势（vdW）：用的是Lennard-Jones 12-6势能。说白了就是常规的吸引力+泡利排斥。这是分子对接中的基础项，几乎每个对接都少不了它。
- 氢键势（Hydrogen Bond）：额外补充的方向性势能。AutoDock等程序把它单独拎出来，因为氢键在药物-靶标结合中太重要了。
我遇到过的最大的坑：格点类型选少了。有个项目里我只开了静电和vdW，结果跑出来的最优构象里的氢键全是歪的。打开AutoDock生成的 .glg 日志才发现，氢键格点根本没参与打分。
# 检查格点日志中是否包含氢键
# 打开对接输出 .dlg 文件，搜索 "grid":
Grid Map: receptor_elec.map (静电 ✓)
Grid Map: receptor_desolv.map (去溶剂 ✓)
Grid Map: receptor_hb.map (氢键 ✓) -> 确认这一行存在关于格点类型的几点个人建议
- 做常规小分子对接：静电 + vdW + 氢键 全开。不要省。
- 做高度疏水的配体（如脂溶性分子）：可以考虑关掉静电，只开vdW和氢键，但不算常用。
- 做金属蛋白体系：额外注意金属配位键的势能项。AutoDock 4默认没有专门处理金属配位的格点类型，这时候就需要调整参数文件了。
还有一个细节：水分子和溶剂效应。部分格点程序（如AutoDock 4）会额外生成desolvation map（去溶剂化格点）。这是跟氢键相关的，通常在格点类型里也算“隐藏选项”。你跑对接时，如果发现配体里有很多极性基团，一定要保证desolvation map被正确计算，否则溶剂效应会被低估。
# 一个典型的AutoGrid参数文件（.gpf）示例
# 通常会包含以下map行：
receptor_a.map # 原子类型A（vdW）
receptor_c.map # 原子类型C（vdW）
receptor_n.map # 原子类型N（vdW）
receptor_o.map # 原子类型O（vdW）
receptor_s.map # 原子类型S（vdW）
receptor_h.map # 氢键供体/受体
receptor_elec.map # 静电
receptor_desolv.map # 去溶剂化你可能会问，为什么氢键格点单独列出来而不是合并到vdW里？原因很简单：氢键的方向性极强（大约10°～30°的偏差就会显著影响能量）。如果混在vdW里，原子类型区分会不精细，方向性势能丢失。所以做到现在，几乎所有主流对接软件都把氢键单独作为一类格点。
4.4 参数协同配置“三合一”清单
最后，我把Box大小、分辨率、格点类型这三点放在一起，整理成一个快速参考清单。你可以把它当成一个检查表，每次配置参数时过一遍。
| 维度 | 检查要点 | 常见盲区 |
| Box大小 | 是否包住了配体的所有可能构象？ | 对柔性配体，Box要额外大1-2 Å |
| Box位置 | 中心点是否对准了结合位点核心残基？ | 盲目用默认原点（0,0,0）是大忌 |
| 格点分辨率 | 晶体结构分辨率 vs 格点分辨率是否匹配？ | 分辨率过高会放大结构噪声 |
| 格点类型（静电） | 带电配体是否开启了静电格点？ | 忽略静电会导致电荷互补缺失 |
| 格点类型（vdW） | 是否覆盖了所有配体原子类型？ | 缺原子类型map，该原子得分为0 |
| 格点类型（氢键） | 极性基团是否依赖氢键？ | 关闭氢键会导致“假结合”高排位 |
我记得有一次帮一个实验室的同学看对接问题，所有参数看着都对，就是结果不对。我排查了半天，最后发现Grid Box 只覆盖了半个结合位点——因为蛋白晶体在包盒子时遇到对称单元边界，自动截断了。解决办法很简单：手动把Box沿Z轴方向再放大2 Å。
所以啊，你虽然可以信任参数配置工具，但验证同样重要。跑完AutoGrid之后，我建议你用可视化软件打开格点地图和蛋白，肉眼检查一下：Box是不是真的包住了你想让配体探索的区域？格点地图边界有没有奇怪的空白？这一步花不了5分钟，但能省下后面几天跑假数据的痛苦。
⚠️ 最后提醒一句
格点参数不是一次性设好就一劳永逸的。不同项目、不同配体，都需要重新评估。你的系统是激酶、蛋白酶还是GPCR？配体是带正电还是中性？有没有糖基化位点或金属离子影响格点计算？这些因素会直接影响你最终的参数选择。
嗯，格点参数这块，大致的框架就这些了。配合好Box大小、分辨率和格点类型，你的对接就跑在了一个正确的起跑线上。