208

PLIP:自动识别蛋白-配体相互作用

PLIP 自动识别氢键、疏水、π 堆积、盐桥等相互作用并生成报告,快速理解结合模式。这篇给出命令行用法、输出解读与它与 ProLIF 的分工。

PLIP(Protein-Ligand Interaction Profiler)自动检测蛋白-配体复合物中的各类非共价相互作用,并直接生成人类可读的报告。它的定位是「拿到一个复合物结构,五秒钟看懂它是怎么结合的」——不需要写代码,特别适合快速分析和做汇报材料。

安装与使用

pip install plip

# 直接分析本地 PDB
plip -f complex.pdb -t --name report      # -t 输出文本报告

# 直接按 PDB ID 联网分析
plip -i 1AKE -t -y                        # -y 同时生成 PyMOL 会话文件

# 输出 XML 便于程序解析
plip -f complex.pdb -x --name result

也有网页版可直接上传结构。-y 生成的 PyMOL 会话文件里相互作用已经标注好,打开即可出图,做幻灯片很省事。

能识别哪些相互作用

类型 典型几何判据
氢键 供受体距离 < 4.1 Å,角度 > 100°
疏水接触 疏水原子间 < 4.0 Å
π 堆积 环心距 < 5.5 Å,分平行与 T 型
π-阳离子 环心与正电中心 < 6.0 Å
盐桥 带电中心间 < 5.5 Å
卤键 符合方向性判据的卤原子接触
水桥 经水分子中介的氢键
金属配位 金属离子与配位原子

水桥和金属配位是 PLIP 比较有价值的两块:很多结合位点的关键相互作用是通过保守水分子中介的,手工看结构容易漏掉,而这类信息对分子设计(要不要置换那个水)很关键。

报告怎么读

文本报告按相互作用类型分节,列出残基编号、类型、距离、角度等。读的时候关注三点:

  • 哪些残基参与最多:出现在多种相互作用里的残基通常是结合的关键锚点,改造分子时要保住与它们的接触。
  • 相互作用类型的分布:以疏水接触为主说明结合靠形状互补,选择性可能不足;有明确的氢键网络和盐桥则通常特异性更好。
  • 有没有水桥:如果关键相互作用依赖某个水分子,设计时可以考虑用分子的一个基团去置换它——这往往是活性提升的机会点,也可能是陷阱(置换保守水的熵代价可能很高)。

和 ProLIF 怎么分工

  • PLIP:单个(或少量)静态结构,要的是快速、可读的报告和现成的出图。命令行一条,无需写代码。
  • ProLIF(见 207《ProLIF》):批量姿势或整条 MD 轨迹,要的是可统计、可做特征的结构化指纹。
  • 实际组合:先用 PLIP 看懂共晶结构的结合模式,明确「关键相互作用有哪些」;再用 ProLIF 在对接结果或 MD 轨迹上批量统计这些相互作用是否被复现、是否稳定。两者一静一动,配合得很自然。

注意事项

  • 依赖输入结构的质量:缺氢原子会严重影响氢键判定。分析前建议先用 PDBFixer(见 215《PDBFixer》)补氢并确定质子化态。
  • 几何判据是经验阈值:临界值附近的相互作用(如 4.0 Å 的氢键)要结合化学常识判断,别把工具输出当绝对真理。
  • 只描述不评价:它告诉你存在什么相互作用,不告诉你结合有多强。亲和力仍需打分函数、自由能计算或实验回答。

上手提示

  • 一条命令看懂结合模式,还能直接出 PyMOL 图,做汇报很省事;
  • 分析前务必补氢并定好质子化态,否则氢键判定不可靠;
  • 水桥信息容易被手工分析漏掉,是它的独特价值;
  • 静态结构用 PLIP,轨迹与批量姿势用 ProLIF。

延伸资源