4. 基于片段的骨架跃迁设计

基于片段的骨架跃迁(Fragment-based Scaffold Hopping, FBSH)与整分子相似性搜索的思路不同:它先把已知活性分子拆成若干“有意义的片段”,再在片段层面寻找可替换的骨架单元,最后把片段重新组装成完整分子。这条路径的价值在于,它允许骨架发生较大变化,同时保留与靶点关键相互作用相关的片段特征。本节围绕三个环节展开:片段库如何构建、片段如何连接与生长、片段命中如何转化为可合成的跃迁候选。

一、片段库构建:从“拆什么”到“存什么”

片段库不是简单地把分子切成小块。构建时首先要确定切割规则。常用的策略包括:按可旋转键切割、按环系与连接链切割、按药效团特征切割,以及按逆合成可及性切割。以可旋转键切割为例,一个含 3 个可旋转键的分子可被切成 4 个片段,但并非所有片段都值得保留。通常需要过滤掉分子量小于 80 Da 的碎片、不含任何环或特征官能团的烷基链,以及无法独立合成的片段。

片段库的存储字段应至少包含:片段 SMILES、分子量、环数、氢键供体/受体数、可旋转键数、logP 贡献值、连接点类型(如单键连接、双键连接、芳环连接)、以及该片段在已知活性分子中出现的频次。连接点类型尤其重要,因为它决定了后续片段能否以正确的几何方向拼接。例如,一个通过芳环对位连接的片段与通过间位连接的片段,在空间取向上完全不同,不能混用。

片段库构建的实用原则:优先保留含环片段和含极性基团的片段,因为它们是决定骨架形状和靶点结合的主要因素;烷基链片段可只保留长度分布统计,不必逐条存储。

二、片段连接与生长:从片段到完整分子的组装逻辑

片段连接与生长是 FBSH 的核心操作。连接指把两个或多个片段通过化学键拼接成完整分子;生长指在片段命中基础上逐步添加原子或小基团,直到达到目标分子量或活性阈值。实际操作中,连接与生长往往交替进行。

连接策略可分为三类。第一类是“头尾连接”,即两个片段各有一个连接点,直接成键。第二类是“桥接连接”,两个片段之间插入一个连接链或连接环,如亚甲基、酰胺键、三唑环等。第三类是“融合连接”,两个片段共享一个环或一个原子,形成稠环体系。生长策略则包括:在片段上添加取代基、延长连接链、引入新的环系、以及进行生物电子等排替换。

判断连接与生长是否合理,需要同时考虑几何和能量两个维度。几何上,连接后的分子应能保持关键片段与靶点残基的相互作用方向;能量上,连接后的构象能不应显著高于片段单独存在时的能量之和。一个常用的检查方法是:对连接后的分子做快速构象搜索,观察关键片段是否仍能落在与原始活性分子相似的结合口袋位置。

基于片段的骨架跃迁设计流程 片段库构建 切割规则 过滤条件 连接点类型 频次统计 输出:片段 SMILES + 属性表 片段连接与生长 头尾连接 桥接连接 融合连接 取代基生长 输出:候选分子集合 跃迁候选转化 几何检查 能量评估 合成可及性 活性预测 输出:优先合成列表 三类连接方式示意 头尾连接 片段A 片段B 直接成键,无插入原子 桥接连接 片段A 连接链 片段B 插入酰胺/三唑/亚甲基等 融合连接 片段A 片段B 共享环/原子,形成稠环 迭代优化:根据评估结果返回调整片段库或连接策略

三、从片段命中到跃迁候选的转化路径

片段命中通常来自两种途径:一是通过片段筛选实验(如 NMR、X 射线晶体学、SPR)得到与靶点弱结合的片段;二是通过计算片段对接或片段相似性搜索得到候选片段。无论哪种途径,片段命中只是起点,真正的跃迁候选需要经过多轮转化。

转化路径可以概括为“四步筛选”。第一步,片段命中确认:检查片段是否确实与靶点结合,结合位点是否与原始活性分子重叠。第二步,片段连接方案生成:根据片段在靶点中的相对位置,枚举可能的连接方式,包括直接连接、桥接连接和融合连接。第三步,候选分子生成与初步评估:对每个连接方案生成完整分子,计算分子量、logP、氢键供体/受体数、可旋转键数等类药性指标,并做快速构象搜索。第四步,优先排序与合成决策:综合几何匹配度、能量合理性、合成可及性和活性预测值,选出优先合成列表。

以一个具体示例说明。假设原始活性分子含有一个苯并咪唑骨架,与靶点激酶铰链区形成两个氢键。片段筛选得到一个含氨基吡啶的片段,晶体结构显示该片段与铰链区形成类似的氢键,但位置略有偏移。此时可以尝试用氨基吡啶片段替换苯并咪唑片段,并通过一个亚甲基桥接连接到分子的其余部分。生成候选分子后,检查其构象:氨基吡啶的氮原子与铰链区残基的距离是否在 2.8–3.2 Å 范围内,苯环平面是否与原始分子保持相似的取向。若满足,则进入合成可及性评估;若不满足,则调整连接链长度或连接点位置,重新生成候选。

转化路径中的关键判断依据:几何匹配度优先于能量评分。因为骨架跃迁的核心目标是保持关键相互作用,几何方向错误时,即使能量评分再低,也难以获得真正的活性分子。

四、实操中的参数与判断标准

在片段库构建阶段,建议设置以下参数:片段分子量范围 80–250 Da;环数 1–3;氢键供体数 ≤3;氢键受体数 ≤6;可旋转键数 ≤4。这些参数并非绝对,可根据靶点类型调整。例如,针对蛋白-蛋白相互作用靶点,可适当放宽分子量上限至 300 Da。

在连接与生长阶段,连接链长度通常控制在 1–4 个原子。过短可能导致片段间几何冲突,过长则增加构象熵损失和合成难度。桥接连接中,酰胺键和三唑环是常用选择,因为它们在药物化学中合成成熟且代谢稳定。

在候选转化阶段,优先排序可参考以下权重:几何匹配度 40%、能量合理性 20%、合成可及性 20%、活性预测 20%。几何匹配度可通过关键片段与靶点残基的距离偏差、角度偏差来量化。能量合理性可用分子力学快速评估。合成可及性可参考 SAscore 或逆合成路线长度。活性预测可使用基于片段的 QSAR 模型或对接打分。

阶段 关键输入 核心操作 输出与判断
片段库构建 已知活性分子、切割规则 切割、过滤、属性标注 片段 SMILES + 连接点类型 + 频次
片段连接与生长 片段库、靶点结构 头尾/桥接/融合连接、取代基生长 候选分子集合,检查几何与能量
跃迁候选转化 候选分子集合 四步筛选、优先排序 优先合成列表,含评分与理由

最后需要强调,基于片段的骨架跃迁不是一次性计算,而是一个迭代过程。片段库可以随着新数据不断扩充,连接策略可以根据评估结果调整,候选分子也可以在多轮循环中逐步优化。把每一轮的输入、操作和判断依据记录下来,才能形成可复现、可积累的设计经验。