如何确定PDB结构结合口袋中最优氢键供体/受体对位置

完全可行，这套定位逻辑和你已经掌握的π-π堆积位点确定思路本质一致：都是基于蛋白残基官能团的固有几何特征做静态推导，不需要跑分子动力学这类重计算，得到的位点精度足够支撑初始配体对接的约束设置、或者结合模式预判。

• 第一步先筛选候选作用原子
  先锁定你之前划定的口袋周边3Å范围内的所有氨基酸，不要只看侧链，主链的酰胺供体、羰基受体是占比很高的强作用位点，必须纳入统计，按规则拆分出所有氢键供体、受体重原子：
  • 氢键供体：连接了氢原子的电负性重原子，包括Ser/Thr/Tyr的羟基氧、Lys/Arg的氨基/胍基氮、Asn/Gln的酰胺氮、His咪唑环上连氢的氮、所有氨基酸主链的酰胺氮、Cys的巯基硫
  • 氢键受体：带孤对电子的电负性重原子，包括Asp/Glu的羧基氧、Asn/Gln的酰胺羰基氧、Ser/Thr/Tyr的羟基氧、His咪唑环上未连氢的氮、所有氨基酸主链的羰基氧
• 第二步沿特征方向计算最优落点
  这一步和你之前找芳香环法向延伸定π-π位点的逻辑完全对应，每个供体/受体的最优作用方向是由官能团杂化方式固定的，不需要算力场参数：

  供体位点计算：沿供体重原子和其共价连接的氢原子的成键方向，朝向口袋空腔（远离残基本体的一侧）延伸。蛋白-配体氢键的最优重原子间距为2.7_{3.3Å，因此从氢原子位置出发，沿成键方向向外走1.7}2.3Å的坐标，就是适配配体上对应受体原子的最优落点。
  受体位点计算：沿受体原子的孤对电子伸展方向，朝向口袋空腔延伸。不同杂化类型的原子孤对方向是固定的：sp2杂化的羰基氧，两个孤对在羰基平面内，与C=O键呈120°夹角分布；sp3杂化的羟基氧，两个孤对在C-O-H键的垂直平面上，与O-H键呈约109°夹角；咪唑环上的受体氮，孤对沿环平面向外，与环上相邻两个共价键呈120°夹角。同样按2.7_{3.3Å的最优重原子间距折算，从受体原子出发沿孤对方向向内走1.7}2.3Å的坐标，就是适配配体上对应供体原子的最优落点。

• 第三步过滤无效位点、标记优先级
  算完初始落点后需要做一轮清洗，避免把蛋白内部的无效位置算进去：
  • 把相互距离小于1Å的重叠落点合并，这类被多个残基共同指向的位点是置信度最高的强结合区
  • 剔除落在蛋白原子范德华半径内的落点：如果落点距离任意蛋白非氢原子小于1.2Å，说明位置在蛋白结构内部，没有配体结合空间，直接删除
  • 给剩余位点打优先级分：带电残基（Asp/Glu/Lys/Arg）对应的位点 > 极性中性残基（Ser/Thr/Asn/Gln/Tyr/His）对应的位点 > 主链原子对应的位点；和你之前定位到的π-π堆积位点距离小于2Å的氢键位点可以额外提权，这类位置是常见的协同作用热点。

• 日常分析用PyMOL就能批量完成所有计算，不需要手动逐个测量：先把3Å范围内的口袋残基选为初始集合，拆分供体、受体原子子集，写个简单的Python脚本遍历原子、按上述几何规则计算方向、生成落点标记即可，核心打标记用cmd.pseudoatom命令就能实现，没有复杂逻辑。
• 不用强行卡完美键长键角：天然解析的蛋白-配体复合物结构中，氢键允许±0.5Å的键长偏差、±30°的键角偏差，你标记的最优落点周围0.8Å半径内的区域，都属于高概率的配体作用范围。

内容的提问来源于stack exchange，提问作者DrTchocky