核酸小分子对接

核酸小分子对接，是每一个想算清楚”这个小分子和 DNA/RNA 到底怎么结合”的人，踩坑最多的对接方向。我刚做 G-四链体配体筛选时，拿蛋白对接那套默认参数直接塞进去，结果 pose 满天飞还带一堆 clash，我居然还当真拿去写结合模式。后来才懂，核酸小分子对接难的从来不是把配体放进沟里，而是核酸那个强电荷环境和金属离子介导，普通对接力场根本没法原样照搬。

一、背景：为什么核酸对接比蛋白难核酸小分子对接

核酸是高度带电的大沟小沟体系，小分子（尤其是阳离子配体）和磷酸骨架的静电主导了结合，这和蛋白口袋以疏水/氢键为主的图景完全不同。做抗肿瘤的人筛 G-四链体配体，做抗菌的人筛 RNA 调控子，做抗病毒的人筛 RNA 元件抑制剂。可以说，核酸小分子对接是你把”配体怎么卡进沟里、靠什么稳定”说清楚的关键。很多新手直接套蛋白对接流程，等到别人问”你的 Mg²⁺ 在哪、离子介导你怎么处理的”，才发现核酸体系最关键的稳定因素被整个漏掉了。更现实的是，核酸柔性大，刚性受体假设下对接出的 pose 常常连物理接触都不对。

二、核心原理：从强静电到离子介导

核酸对接的核心难点在静电：骨架磷酸带负电，配体若带正电，长程库仑主导结合姿态，而经典对接的加和式打分对长程静电和去溶剂化描述粗糙。金属离子（尤其 Mg²⁺、K⁺）往往桥接配体和碱基，是结合模式的决定因素，但默认对接不处理离子介导。关键认知是，核酸-配体相互作用是”静电主导 + 沟槽形状匹配 + 离子桥接”的复合，单纯形状互补打分会严重失真。还有一个点：RNA 比 DNA 柔性更大、碱基可翻转，小沟大沟的几何随序列剧烈变化，受体刚性假设下的对接姿态可信度低，需要半柔性甚至柔性受体。

三、关键技术要点：电荷、离子与柔性受体

电荷处理是第一关：核酸力场（如 AMBER 的 OL15 核酸力场）要正确加原子电荷，配体也要用匹配残基的电荷模型，别用蛋白对接那套忽略核酸特异的电荷。离子必须显式建模——Mg²⁺/K⁺ 桥接要用显式离子或在对接中约束离子位置，很多 G-四链体配体就是靠 K⁺ 桥稳定。受体柔性上，我习惯先做含离子的短 MD 采样核酸构象，再对采样的若干受体构象分别对接，取一致的 pose。我早年算一个 RNA 适配体配体，刚性受体对接出一个漂亮 pose，含离子 MD 一跑发现配体实际被离子顶到了沟外——刚性假设直接把结合模式搞反。另一个常被忽略的点：溶剂化，核酸沟里的有序水对结合很关键，隐式溶剂下会丢水介导的相互作用。

四、实操流程：从核酸到可信 pose核酸小分子对接

第一步，准备核酸结构（加离子、中和、用核酸力场最小化）；第二步，含离子的短 MD 采样，取若干代表性受体构象；第三步，准备配体（正确质子化、带匹配电荷）；第四步，对每个受体构象做对接（用支持核酸的力场/参数，或加离子约束）；第五步，对 top pose 做含离子的短时 MD 验证，确认结合模式稳定而非对接假象。我习惯把”多个受体构象一致 + MD 验证稳定”作为可信 pose 的标准，单构象对接出的结果我只当线索。对 G-四链体还要专门确认拓扑和离子配位，对双链沟槽要结合小沟/大沟的几何差异。

核酸对接我踩过的最大坑，是漏掉离子介导。蛋白对接那套刚性受体直接搬过来，pose 满天飞还带 clash，看着还挺美，含离子 MD 一跑才发现配体被 K⁺ 顶到了沟外。现在核酸对接必带离子做 MD 验证，没有离子介导验证的 pose 我不认，因为那往往只是力场没算离子时的假象。

配体质子化在核酸环境里也要重新确认。沟里偏负电，配体在结合态的质子化态可能和溶液不同，用错电荷对接姿态直接偏。我先确认结合态电荷再对接，这一步省不得。核酸不是带电的蛋白，它的电荷环境和离子共同导演了结合，照搬蛋白流程一定会翻车。

小沟大沟的几何差异我也一定先确认。同一个配体在小沟和大沟的结合模式可能完全不同，盲对接容易把配体塞进几何不匹配的沟。我交付前会比对结合位点的沟槽宽度和电荷，确认 pose 落在物理合理的沟里，而不是只信打分最高的那一个。

五、常见问题

Q1：pose 满天飞还 clash？刚性受体假设错，核酸要半柔性+离子。 Q2：结合模式不稳定？没做含离子 MD 验证，补短时 MD 确认。 Q3：电荷对不上？用了蛋白力场电荷，换核酸力场重算。 Q4：离子介导漏了？没显式建模 Mg²⁺/K⁺，加离子约束再对接。 Q5：RNA 结果飘？柔性大，多构象采样后对接取一致。 Q6：去溶剂化描述差？换更准的打分或 MM-GBSA 精修。 Q7：小沟大沟搞混？结合位点几何要先确认，别盲对接。

讲一个具体的坑：有次算一个双链小沟结合配体，蛋白对接式的刚性受体给出配体平躺在小沟里、还带两个氢键，看着挺美。结果含 K⁺ 的 MD 一跑，K⁺ 把配体顶去了大沟外侧，原本的氢键全断——那个漂亮 pose 是力场没算离子时的假象。从那以后我固定一条规矩：核酸对接必带离子做 MD 验证，没有离子介导验证的 pose 我不认。另一个常被忽略的点：配体质子化，核酸环境偏负电，配体在沟里的质子化态可能和溶液不同，我一定先确认结合态电荷再对接。

六、复盘总结

回过头看，核酸小分子对接最容易被当成”蛋白对接换个受体”的照搬活，但它真正的难点在强电荷、离子介导和受体柔性：这三样少处理一样，pose 就可能是假象。我现在的习惯是核酸力场就位、离子显式建模、多构象对接加 MD 验证，三步少一步都不敢把结合模式写进结论。核酸不是带电的蛋白，它是静电和离子共同导演的结合剧场。把电荷、离子、柔性三件事钉死，对接出的 pose 才敢拿去解释调控机理。