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。
延伸资源
- 对照工具:207《ProLIF》;结构修复:215《PDBFixer》;
- 可视化:211《PyMOL 开源版》;相互作用概念见「结构与模拟」模块。