BioSolveIT HYDE：打分函数与结合能解释

HYDE（HYdrogen bond and DEhydration）是 BioSolveIT 的打分函数，与多数经验打分函数不同，它建立在一个明确的物理图景上：配体结合的净能量 = 形成的相互作用（主要是氢键）− 脱去水化壳层的代价。这个清晰的物理框架让它的分数可以分解到每个原子，从而变得可解释。

物理基础

配体从水中进入蛋白口袋，要经历两个能量过程：

• 脱水（去溶剂化）：配体和口袋原本都被水分子包围，结合时这些水必须被排开。这是要付出能量代价的，极性原子的代价尤其高（它们与水形成了良好的氢键）。
• 形成新相互作用：配体与蛋白之间形成氢键、疏水接触等，释放能量。
• 净效果：两者之差决定结合强度。

这个框架解释了一个重要的设计原则：把极性基团放进疏水口袋是双重损失——付出了高昂的脱水代价，却没有得到补偿的氢键。反之，疏水基团进入疏水口袋，脱水代价低（疏水原子与水的相互作用本来就弱），是「便宜」的结合。

可解释性带来什么

常规打分	HYDE 式分解
「这个分子得分 −9.2」	「这个羟基脱水代价高但没形成氢键，扣了 1.5」
不知道该改哪里	直接指向具体的改造位点
只能比较分子间	可比较同一分子内各部分
黑箱	可与化学直觉对照验证

在 SeeSAR（见 303《BioSolveIT SeeSAR》）中，这种分解被可视化成原子上的彩色球——绿色表示该原子对结合有利，红色表示不利。药化看一眼就知道该动哪里。

局限

• 仍是经验估算：脱水项与氢键项的参数是拟合出来的，不是从第一性原理算的。绝对值不可当实测亲和力用。
• 对熵的处理简化：构象熵、水分子释放的熵贡献处理得比较粗糙，而这些在真实结合中很重要。
• 不处理特殊相互作用：卤键、π-阳离子、金属配位等的建模不如氢键充分。
• 受体刚性假设：不考虑诱导契合。
• 依赖正确的姿势：姿势错了，分解分析全错。
• 共同的边界：与所有经验打分函数一样，与实测活性的相关性有限，用于指导设计方向比用于预测排名更可靠。

用开源做能量分解

# MM/GBSA 逐残基分解是最常用的开源路径
# gmx_MMPBSA（基于 AmberTools，支持 GROMACS/OpenMM 轨迹）

mamba install -c conda-forge gmx_mmpbsa

# mmpbsa.in 配置
&general
  startframe=100, endframe=500, interval=5,
  verbose=2,
/
&gb
  igb=8, saltcon=0.150,
/
&decomp
  idecomp=2,                    # 逐残基分解
  dec_verbose=3,
/

# 运行后得到每个残基对结合能的贡献，
# 可以看出哪些残基是主要贡献者、哪些是不利的

# 更轻量：用 ProLIF（207）做相互作用层面的分解
# 虽然不给能量数值，但能给出
# 「哪个残基、什么类型的相互作用、占有率多少」
# 对设计的指导价值往往不亚于能量分解，且更稳健

设计中的实用原则

无论用什么工具，这些由脱水-氢键框架推出的原则都成立：

• 埋藏的极性原子必须形成氢键，否则是净损失；
• 疏水接触是「便宜」的结合能，但过度增加疏水性会损害溶解度与选择性（配体高效性的意义正在于此）；
• 置换结合位点的水分子可能带来显著增益，但要看那个水是否「高能」（与蛋白氢键不佳的水容易被置换且有利）；
• 不要只加基团：每加一个原子都有代价，看每重原子的效率（LE）而非绝对分数。

延伸资源

• 配套工具：303《BioSolveIT SeeSAR》、305《BioSolveIT FlexX》；打分函数概念：095《Docking Score》；
• 相互作用分析：207《ProLIF》、208《PLIP》；自由能：201《OpenFE》。