【药靶作用力分析】氢键 —— 药物设计中最熟悉的陌生人

经典定义我们都熟：氢原子跟一个电负性强的原子（供体，D）形成共价键之后，又跟另一个电负性原子或者π电子体系（受体，A）产生非共价吸引。最典型的D/A组合就是N–H···O=C和O–H···O。

但在药物-蛋白界面，我们需要用更精细的几何与能量语言来理解它：

在分子动力学模拟和晶体结构解析里，氢键通常要同时满足两个条件才算数：

重原子D到A的距离 < 3.5 Å（部分软件采用3.0–3.2 Å更严格标准）

D–H···A 夹角 > 120°（越接近180°越理想）

参考《Journal of Medicinal Chemistry》上大量结合自由能分解的工作，蛋白-配体界面的氢键大致能分成三类：

例如酰胺–酰胺之间，单个氢键对结合自由能的贡献一般在0.5–1.5 kcal/mol。看似很弱，但当配体形成3–5个这样的氢键，就能贡献10倍以上的亲和力差异。

当供体或受体带有形式电荷时（如羧酸根–胍基），单个氢键的贡献可以飙升至3–6 kcal/mol，甚至达到共价键强度的十分之一。盐桥就是一种极端的电荷辅助氢键，我们会在后面专门讲。

C–H···O、C–H···π 等，强度通常 < 0.5 kcal/mol，但在疏水口袋中大量累积后，对结合构象的预组织和选择性不可忽视。

这是药物设计里最容易低估的一笔“隐性成本”。结合之前，配体的极性基团和蛋白的极性残基都泡在水里，靠氢键被水分子溶剂化。结合时，这些水分子要被排开，形成蛋白-配体氢键释放出的焓，首先得填上去溶剂化消耗的能量。很多氢键算下来净能量甚至“不赚不赔”，但它们能精准地锚定结合构象 —— 所以千万别简单粗暴地认为“氢键越多越好”。

在实际药物发现项目中，氢键并不是“画出来”的，而是计算和统计的结果。

最基础的是通过PDB晶体结构或对接pose，用几何判据来数氢键。PyMOL、VMD、MOE、Schrödinger Maestro这些软件都有内置工具。比如VMD的hbonds插件，默认用3.0 Å的距离和20°的角度偏离（也就是D–H···A > 150°）做标准。你还可以在模拟轨迹里统计某个氢键的占据率 —— 这个数字对判断它到底重不重要，非常关键。

静态结构只是一张快照。一个氢键稳不稳定，要看它在100 ns甚至更长的MD模拟里，占有率有没有超过80%。典型的分析管线是：对MD轨迹的每一帧识别氢键 → 统计“氢键生存率” → 绘制氢键热图。有些晶体结构里看到的氢键，在溶液里占有率可能不到30%，那很可能只是结晶堆积的假象，不是真正驱动结合的力。

在计算化学这边，SAPT（对称匹配微扰理论）和DFT能量分解可以把分子间相互作用拆成静电、交换排斥、色散和诱导项。氢键主要落在静电项上，但电荷转移的贡献经常被忽略。近年的ALMO-EDA等方法显示，某些强氢键里电荷转移能能达到2–5 kcal/mol，已经接近弱共价作用的范畴了。

一键生成配体-蛋白相互作用指纹，氢键会用不同颜色标出是主链还是侧链、中性还是电荷型。

免费网页工具，也有Python库，上传PDB就能自动识别氢键、卤键、π堆叠等。

对新手友好，做高质量的二维相互作用图，发论文很实用。

氢键对结合稳定性的贡献，最直观地体现在“突变-亲和力”实验里。如果蛋白侧链上一个关键氢键残基发生突变，配体亲和力掉个10–100倍，这类残基就被称为“热点”。

比如HIV蛋白酶抑制剂里，跟催化天冬氨酸（Asp25/25'）形成的氢键网络就是保守的稳定化核心，任何一个抑制剂要是少了跟Asp25的氢键，哪怕疏水部分填得再好，活性也会断崖式下跌。

所有他汀（比如阿托伐他汀）都带一个3,5-二羟基庚酸片段，羧基跟酶活性位点的Lys735形成电荷辅助氢键，旁边的羟基再跟Ser684的主链羰基形成一连串氢键。这个氢键网格提供了结合的核心焓驱动力，保证整个分子的长疏水尾部能精准地插进底物口袋。

同一家族激酶的ATP结合口袋长得非常像，想设计高选择性抑制剂，最有力的策略之一就是利用非保守残基形成的独特氢键。

COX-1和COX-2的活性位点只有一个关键残基的差别：COX-1的523位是异亮氨酸，COX-2是缬氨酸。这个微小变化在口袋侧面多出一个亲水小囊。塞来昔布的磺酰胺基团恰好伸进这个囊里，跟Arg513和His90形成特有的氢键网络，而这些残基在COX-1里根本够不着。正是这组选择性氢键，让塞来昔布对COX-2的选择性达到COX-1的375倍，从而避免了传统非甾体抗炎药常见的胃肠道副作用。

引申：激酶抑制剂的“铰链区氢键”策略

几乎所有ATP竞争性激酶抑制剂都必须跟激酶铰链区的主链形成1–3个氢键，模拟腺嘌呤环的识别方式。这部分氢键保证了基础亲和力，而选择性则来自铰链区以外、跟周围特异残基形成的额外氢键。 比如伊马替尼，除了铰链区氢键，它的苯酰胺基团还能跟Abl激酶独有的Glu286形成关键氢键，再配合疏水作用，从而在茫茫激酶组里精准地抑制Bcr-Abl。

氢键不仅贡献亲和力，某些情况下它直接决定了药物是激动剂还是拮抗剂。

全反式维甲酸（ATRA）与RAR结合时，其羧基与受体形成氢键网络，诱导螺旋H12发生“关盖”构象，招募共激活因子。而当某些拮抗剂结合时，虽然同样占据配体口袋，但其庞大的侧链阻碍了H12与配体间某个关键氢键的形成，导致H12采取不同构象，转而招募共抑制因子。一个氢键的有无，决定了转录机器的走向。

讲了这么多原理，我们看几个从靶标发现到临床药物的真实优化路径。

早期HIV蛋白酶抑制剂大多是拟肽分子，带着好几个酰胺键，能跟催化中心Asp25/Asp25′形成稳定的氢键，活性很高，但口服生物利用度差，还容易被降解。

后来研究人员用羟基乙胺这类过渡态模拟结构去替代易水解的酰胺键，在保留关键氢键的同时大幅提升了代谢稳定性。从沙奎那韦到安普那韦的优化过程中，引入四氢呋喃氨基甲酸酯片段，形成了更稳定的外围氢键网络，同时增加了分子刚性，最终改善了口服吸收和药代性质。

这个案例告诉我们，保持关键氢键的同时优化分子骨架和结合模式，完全有可能同步提升活性和DMPK性质。

EGFR发生T790M突变后，第一代EGFR抑制剂很容易失效，所以克服耐药就成了第三代药物设计的关键。奥希替尼没有死磕突变位点附近的相互作用，而是保留了跟EGFR铰链区Met793主链形成的关键氢键，同时引入丙烯酰胺基团，跟Cys797形成不可逆共价键，这样即使在突变体里也能保持较高的结合能力和抑制活性。

这个案例提醒我们，当关键结合环境发生变化时，与其固守原有氢键，不如重新构建整体结合模式，从而有效识别耐药突变体。

❶ 不要迷恋“氢键数”—— 氢键多不等于药好。每增加一个极性氢键供体或受体，都要认真评估它带来的去溶剂化惩罚和透膜性损失。

❷ 优先瞄准主链氢键：侧链容易突变产生耐药，而跟主链酰胺形成的氢键在靶标突变中通常更保守，是长效抑制剂的可靠根基。

❸ 用“氢键占据率”代替“晶体里看到几个氢键”来做决策。在微秒级的全原子MD模拟里，只有占据率超过80%的氢键才可能成为真正的结合热点，值得投入合成资源去优化。

❹ 巧用水分子介导的氢键：如果蛋白口袋某个区域总是“住”着一个稳定的水分子，不妨考虑让配体跟这个水分子形成氢键桥，而不是非得把它排走。这种策略有时候能带来意想不到的熵增收益。

❺ 氘代的隐秘效应值得留意：近年研究表明，碳–氘键比碳–氢键更短，振动频率更低，在极少数体系里，氘代能微调C–H···O这类弱氢键的强度，从而影响选择性。虽然还属于前沿，但已经有好几篇《Science》报道过它在药物优化中的应用了。

梳理下来，氢键的真实分量不取决于数量，而在于“位置”和“代价”—— 前者决定能否锚定关键残基，后者涉及去溶剂化惩罚与构象熵的取舍。

他汀、COX-2抑制剂和奥希替尼等案例反复印证：氢键设计贵在精准利用，而非盲目堆砌。然而，许多在晶体中看似完美的氢键，在液相环境中可能占据率极低，甚至只是结晶堆积的假象。因此，要真正将氢键转化为可控的设计工具，必须依赖动力学模拟来验证其稳定性，并结合能量分解明确其净贡献。

而“动态构象验证、结合能定量、残基级贡献解析”这些需求恰是当前药物研发中从“静态对接”迈向“动态精准”的关键瓶颈。为此，我们向您推荐：

点击文末阅读原文，

获取谱度众合专属分子动力学模拟相关实验方案~

告别氢键这个最熟悉的角色，下一期我们把目光投向一个长期被低估的选手——卤键。它的本质是σ-hole驱动的方向性静电作用，在靶点选择性上有自己独特的优势。你会发现，一个溴原子或碘原子在骨架上的位置挪动哪怕几埃，最终的活性数据就可能差出三个数量级。下篇我们接着聊。

参考资料