静电互补 Electrostatic Complementarity：为什么电荷分布影响结合

形状互补决定「能不能装进去」，静电互补决定「装进去之后能不能稳定结合」。它是特异性与选择性的主要来源——也是最容易在设计中被搞错的维度。

静电相互作用的类型

类型	强度（真空中）	距离依赖	特点
盐桥	强	1/r	需要电荷互补；有方向性要求
氢键	中	短程	方向性最强，特异性的主要来源
偶极-偶极	弱到中	1/r³	—
阳离子-π	中	—	与芳环的富电子面作用
π-π 堆积	弱到中	—	面对面或边对面
卤键	弱到中	—	高度方向性（σ 空穴）

氢键的方向性是特异性的关键：疏水接触几乎不挑剔几何，而氢键要求供体-氢-受体接近直线、且受体的孤对电子朝向合适。这个严格的几何要求，正是选择性的来源。

去溶剂化：静电作用的隐藏代价

# 关键认识：
#   在水中，静电相互作用的净贡献远小于真空中的计算值
#
# 原因：
#   带电或极性基团在水中已经与水形成了良好的相互作用
#   结合时必须先脱去这些水（去溶剂化，代价高）
#   然后才能与蛋白形成新的相互作用（收益）
#
#   净贡献 = 新形成的相互作用 - 去溶剂化代价
#
# 后果：
#   1) 埋藏的极性基团如果没有形成氢键，是净损失
#      → 「未配对的埋藏极性原子」是设计中的严重问题
#
#   2) 盐桥在暴露于溶剂的位置贡献很小
#      （水可以屏蔽电荷）
#      → 埋藏的盐桥贡献才大
#
#   3) 疏水基团的去溶剂化代价低
#      → 这是「疏水结合便宜」的物理原因
#      → 也是 logP 与活性正相关的原因（见 072）
#
# 设计含义：
#   引入极性基团时必须确保它能形成相互作用，
#   否则不如用疏水基团

静电互补的可视化与量化

# 1) 静电势面（ESP surface）
#    计算分子表面的静电势，用颜色显示
#    红色 = 负电（富电子），蓝色 = 正电（缺电子）
#
#    好的互补：配体的红色区域对应蛋白的蓝色区域
#
#    计算方法：
#      - 量化（Psi4、ORCA、xtb）算精确的 ESP
#      - 或用力场的部分电荷近似（快但粗糙）
#
#    在 PyMOL 中显示：
#      set surface_ramp_above_mode, 1
#      ramp_new esp_ramp, protein_map, [-5, 0, 5], [red, white, blue]
#      set surface_color, esp_ramp, protein

# 2) 静电互补性得分
#    在界面上逐点比较配体与蛋白的静电势
#    计算相关系数或互补性指数
#    → 负相关 = 好的互补（一正一负）

# 3) 用 RDKit 算部分电荷（快速近似）
from rdkit import Chem
from rdkit.Chem import AllChem

mol = Chem.MolFromSmiles("CC(=O)[O-]")
mol = Chem.AddHs(mol)
AllChem.ComputeGasteigerCharges(mol)
for atom in mol.GetAtoms():
    q = atom.GetDoubleProp("_GasteigerCharge")
    print(f"{atom.GetSymbol()}{atom.GetIdx()}: {q:+.3f}")

# Gasteiger 电荷是经验方法，快但粗糙
# 精确计算需要量化（如 RESP、AM1-BCC，见 202）

常见的设计错误

• 把未配对的极性基团埋进疏水口袋。这是最常见也最严重的错误。去溶剂化代价高但没有补偿，净效果是降低亲和力。检查方法：在共晶或对接结构中，找出所有埋藏但未形成氢键的极性原子；
• 忽略电荷的方向性：盐桥不只要求电荷相反，还要求几何合适；
• 过度依赖表面盐桥：暴露于溶剂的盐桥贡献很小；
• 忽略 pKa 的变化：基团在口袋中的 pKa 可能与在水中不同——埋在疏水环境中的羧酸可能不电离，附近有正电残基则可能促进电离。这会改变实际的电荷状态；
• 忽略水介导的相互作用：许多「看起来没有直接接触」的极性基团实际上通过水桥相互作用（见 108《PLIP》）。

利用静电获得选择性

# 静电是选择性设计的有力工具
#
# 策略：
#   1) 找出靶点与脱靶在结合位点的电荷差异
#      如：靶点有 Asp 而脱靶是 Asn（同样大小，电荷不同）
#      → 设计带正电的基团与 Asp 作盐桥
#      → 在脱靶上没有这个作用，甚至可能不利
#
#   2) 利用局部的静电环境差异
#      即使残基相同，周围环境不同也会改变局部静电势
#
#   3) 利用 pKa 差异
#      某些残基在不同蛋白中的 pKa 不同
#
# 优势：
#   静电相互作用的方向性与距离依赖强
#   → 提供比疏水作用更好的判别力
#
# 代价：
#   引入带电基团通常损害通透性（见 060）
#   → 需要在选择性与成药性间权衡

卤键：被低估的相互作用

• 机制：与碳相连的重卤素（Cl、Br、I）在 C–X 键的延长线上有一个正电区域（σ 空穴），可与蛋白的负电原子（羰基氧、氮）作用；
• 高度方向性：要求 C–X⋯Y 接近 180°，因此特异性好；
• 强度：I > Br > Cl >> F（氟不形成卤键——它电负性太强，没有 σ 空穴）；
• 设计价值：在合适的几何位置引入氯或溴，可以获得额外的结合能而不显著增加极性；
• 注意：多数常用的对接打分函数对卤键建模不足，可能低估其贡献。

延伸资源

• 形状互补：093《形状互补 Shape Complementarity》；口袋分析：091《蛋白口袋 Protein Pocket》；
• 相互作用分析：108《PLIP》、109《ProLIF》；LLE：072《脂溶性配体效率 LLE》。