# MOE基于片段的药物设计（一）：Scaffold Replacement——如何在保留取代基空间关系的同时替换核心骨架？
在先导化合物优化过程中，我们经常会遇到这样一个问题：
一个分子的两端已经能够与蛋白形成不错的相互作用，但中间的核心骨架并不理想。我们希望保留原有的关键取代基，同时寻找一个新的核心结构，以探索新的化学空间。
这就是Scaffold Replacement（骨架替换）要解决的问题。
在 MOE 的 Fragment-Based Drug Design 工作流中，Scaffold Replacement 是利用已有蛋白–配体复合物的三维结构，寻找能够保持原有取代基空间排布的新骨架。
本文以 MOE 教程中的1KV2 蛋白–配体复合物为例，介绍 Scaffold Replacement 的基本思路和操作流程。
1. Scaffold Replacement 要解决什么问题？
可以先把一个配体抽象成：
其中：
Scaffold：分子的核心骨架；
Fragment：从骨架向外延伸的取代基。
假设Fragment已经分别占据了两个重要的结合区域，并形成了较好的蛋白–配体相互作用。
我们现在希望：
也就是说：
改变中间核心，但尽量不要破坏两边已有的空间排布。
这样做可能用于探索不同的分子骨架、改善化合物性质或进一步优化活性。
原理补充：它与骨架跃迁（Scaffold Hopping）有什么关系？
Scaffold Replacement 可以理解为一种典型的Scaffold Hopping（骨架跃迁）思路。
重点是：
新骨架能否在三维空间中把原来的关键取代基放在相近的位置和方向。
因此，一个二维结构看起来完全不同的骨架，只要能够维持合适的三维几何关系，也可能成为有价值的替换方案。
2. MOE 是怎么寻找新骨架的？
Scaffold Replacement 的重要概念是：
Bond Vector （键矢量）
删除 scaffold 后，保留下来的 R1 和 R2 与原骨架之间会留下两个连接方向。
这些连接方向可以粗略理解为：
MOE 会利用这些bond vectors构建一个三维查询。
随后在 fragment/linker database 中寻找：
有没有fragment，其连接点距离、方向和空间排布能够与这些 bond vectors 匹配？
因此 Scaffold Replacement 的核心是：
“哪个新骨架能够以合适的三维几何关系重新连接这些保留部分？”
这也是为什么这种方法必须基于三维结构来理解。
3. 为什么开始计算之前要先看原始结合模式？
骨架替换第一步应该是：
弄清楚原配体哪些部分值得保留。
首先载入：
$MOE/sample/mol/1kv2.pdb.gz然后依次进行：
1. QuickPrep；
2. SiteView；
3. 查看蛋白–配体相互作用；
4. 打开 Ligand Interactions；
5. 查看 Receptor Surface。
这些操作的目的是在回答两个问题：
哪些相互作用已经很好？
如果某个官能团形成了非常明确的氢键或其他重要相互作用，那么后续设计时应该谨慎破坏它。
结合口袋还有什么空间？
Receptor Surface 可以帮助观察配体周围是否存在尚未充分占据的区域，以及不同取代基目前处在怎样的空间环境中。
在1KV2的例子里，可以关注到配体与Glu71和Asp168附近的关键氢键相互作用。
这些相互作用后面还会再次出现。
4. 第一步：准备蛋白–配体复合物
打开 1KV2 后，使用：
MOE | RHS | QuickPrep
对体系进行准备。
QuickPrep 在这里主要服务于后续蛋白–配体结构优化，包括：
结构准备；
质子化；
处理部分水分子；
对不同区域设置 Fixed、Tethered 和 Free atoms；
为后续 MM minimization 准备体系。
对于本文的 Scaffold Replacement 来说，不需要逐项修改默认参数，按照教程使用默认设置即可。
为什么这一步不能省？
后面的新骨架可能需要在结合口袋中进行结构优化。
如果输入结构本身存在不合理的质子化、缺失氢原子或明显几何问题，那么后续得到的结果也会受到影响。
因此：
可靠的起始结构，是所有 structure-based fragment design 的基础。
5. 第二步：确认原配体的关键相互作用
使用：
MOE | RHS | Ligand | Ligand Interactions...
打开 Ligand Interactions 图。
这里区分：
哪些官能团负责主要相互作用；
哪部分是准备保留的；
哪部分属于准备替换的 scaffold。
对于 1KV2，可以关注到配体与Glu71和Asp168附近的关键氢键相互作用。
这一步实际上是在给后面的骨架替换划边界：
保留部分 —— 准备替换的部分 —— 保留部分6. 第三步：打开 Scaffold Replacement
进入：
MOE | Compute | Fragments | Scaffold Replacement...
打开Scaffold Replacement面板。
这个窗口后面还会反复出现，因为 MOE 的多个 Fragment-Based Design 功能都整合在类似的工作框架中。
窗口中比较重要的区域包括：
Ligand；
Receptor；
Operation；
Linker Database；
Descriptor Filter；
Model File；
Pharmacophore；
Output Database；
Minimize generated structures；
Scoring。
本文先使用最基础的 Scaffold Replacement，不加入 Pharmacophore Filter。
7. 第四步：告诉 MOE“我要换哪一部分”
打开配体的 2D 显示：
MOE | Footer | 2D
然后用鼠标框选中央 scaffold。
选中后，如果这个查询定义合理，MOE 会：
用绿色标出待替换 scaffold；
显示对应的 bond vectors。
教程中的查询显示：
这里最需要理解的是Exit Vectors。
Exit Vector 是什么？
可以把它理解成：
待删除部分与保留结构之间的“出口方向”。
这个例子有两个 Exit Vectors，意味着：
保留结构 A │ ↓ [待替换 Scaffold] ↑ │ 保留结构 B新的 scaffold 必须能够同时连接这两端。
因此，定义 scaffold 的过程实际上同时定义了：
MOE 后面搜索新 fragment 时要满足的三维几何条件。
8. Linker Database：新骨架从哪里来？
教程使用的是 MOE 默认数据库：
$MOE/opt/linker.mdb其中包含大量 linker 的三维构象。
这里非常容易产生一个误解：
MOE 是不是现场“生成”了一个全新的 scaffold？
在这个 Scaffold Replacement 示例中，更准确的理解是：
MOE 根据当前 bond vector query，在指定的 linker/fragment database 中搜索能够满足几何条件的候选结构。
因此搜索结果会受到 fragment database 本身的影响。
如果换成公司内部构建的 fragment/linker database，搜索到的化学空间也会随之改变。
这一点在之后的教程中还会专门介绍。
9. Descriptor Filter 为什么突然变红了？
教程默认的 Descriptor Filter 中包含：
Weight<500 TPSA[40,140] !reactive运行前点击：
Check Filters
MOE 会检查当前 ligand 与这些过滤条件是否兼容。
教程中的 1KV2 ligand 分子量为：
Weight = 528.68已经超过了：
Weight < 500所以这里出现了警告。
教程随后将条件修改为：
Weight<600 TPSA[40,140] !reactive再进行搜索。
这个细节很重要。
Check Filters 是在提醒：
当前设置的筛选规则可能会把与你起始配体处于相似大小范围的新设计直接过滤掉。
如果起始分子本身已经超过 500，而你仍强制要求所有新结构 Weight<500，那么搜索空间可能会被过度限制。
因此 Filter 参数不应该机械照抄。
应该先问：
我现在希望通过这个筛选条件解决什么问题？
10. 开始 Scaffold Replacement 搜索
设置完成后点击：
Search
MOE 会开始遍历 linker database 中的候选结构。
运行过程中可以看到：
已搜索的 linkers 数量；
已生成的 unique molecules 数量。
最终结果写入：
scaffrepl.mdb并在 MOE Database Viewer 中打开。
到这里发生的事情可以概括成：
原始配体 ↓ 指定待替换 scaffold ↓ 提取 exit / bond vectors ↓ 搜索 3D linker database ↓ 找到几何匹配的新 scaffold ↓ 重新连接保留的 R-groups ↓ 生成新的 analogs11. 怎么看 Scaffold Replacement 的结果？
在 Database Viewer 中：
File | Browse...
打开 Database Browser。
随后逐个浏览生成的结构，并与原始 1KV2 ligand 比较。
真正应该看的是：
① 中央 scaffold 是否真的发生改变？
这是最基本的问题。
② 两端原有取代基的空间方向是否基本保留？
Scaffold Replacement 的价值就在这里。
如果新 scaffold 导致两端 R-groups 完全偏离原来的结合区域，那么即使二维结构看起来合理，也不是我们想要的结果。
③ 有没有明显的蛋白–配体冲突？
几何上能够连接，并不意味着整个新分子在口袋中一定合理。
④ 原来的关键相互作用还在不在？
这是本教程第一次 Scaffold Replacement 最重要的结果。
12. 第一次计算暴露了一个关键问题
基础 Scaffold Replacement 的确生成了一批具有明显化学多样性的 analogs。
而且它们能够：
替换中央核心，同时保持 R-groups 大致相同的相对空间方向。
但问题也随之出现：
部分新结构失去了原配体与 Glu71、Asp168 之间的重要氢键。
这揭示了 Scaffold Replacement 一个非常重要的边界：
满足 bond vector 的几何关系，并不等于一定能够保留原来的蛋白–配体相互作用。
换句话说：
几何能够接上 ≠ 关键药效特征一定保留 ≠ 实验活性一定提高
这也是为什么教程下一步没有继续盲目增加搜索数量，而是引入：
Pharmacophore Filter
用药效团约束告诉 MOE：
新骨架不仅要接得上，还要在关键位置保留指定的氢键受体和氢键供体特征。
这会是下一篇文章的核心内容。
13. Scaffold Replacement 的完整思路
最后把这一篇压缩成一张流程：
14. 这项功能什么时候比较适合用？
情况一：核心骨架希望发生变化
例如已有活性分子，但希望探索不同的 central core。
情况二：两端取代基已经有比较明确的结合模式
这时我们希望：
“两边尽量不动，中间换一种连接方式。”
情况三：希望探索新的化学空间
Scaffold Replacement 可以产生与原分子二维结构不同，但三维取代基排列相近的 analogs。
15. 使用时需要注意什么？
① 起始 binding mode 必须尽可能可靠
如果原始 ligand pose 本身就不可靠，那么根据它定义的 bond vectors 也没有可靠基础。
② Bond vector 只是几何约束
它不能自动保证：
原有氢键一定保留；
新分子的亲和力一定更高；
新分子一定能合成；
ADME 性质一定改善。
③ 不要只看生成数量
“生成 100 个”不意味着比“生成 10 个”更好。
最终仍需要结合：
protein–ligand interactions；
pocket fit；
分子性质；
后续 minimization / scoring；
medicinal chemistry judgment
进行判断。
④ Score 只能用于辅助排序
后续即使加入 GBVI/WSA dG 等评分，也不能把评分值直接理解为实验结合自由能或实验活性。
计算结果更适合用来：
帮助缩小候选范围和提出下一轮实验设计假设。
小结
MOE 的 Scaffold Replacement 可以概括成一句话：
删除现有配体的一段核心骨架，根据保留结构留下的三维 bond vectors，在 fragment/linker database 中寻找能够维持原有取代基空间排布的新 scaffold。
它解决的是：
“在尽量保留已有三维结合信息的前提下，我还能用什么新的化学骨架连接这些关键结构？”
基础 Scaffold Replacement 主要解决的是几何匹配问题。
但 1KV2 的例子进一步告诉我们：
几何匹配并不能保证关键氢键等药效特征被保留。
因此下一篇将继续介绍：
Pharmacophore Filter：怎样让骨架替换不仅“接得上”，还尽量保留原来的关键蛋白–配体相互作用？