BioSolveIT ReCore：骨架跃迁和核心替换工具

ReCore 是 BioSolveIT 的骨架跃迁工具。它解决一个具体的设计问题：保留分子中与蛋白形成关键相互作用的取代基，只把中间的连接核心换掉。这是先导优化中很常见但技术上不简单的需求。

骨架跃迁的动机

动机	说明
改善性质	原核心太脂溶、溶解度差、代谢不稳定
规避专利	进入未被覆盖的化学空间（法律判断仍需专利律师）
改善选择性	换核心可能改变与脱靶蛋白的匹配
简化合成	换成更易合成的核心
调整构象	改变取代基的相对朝向与刚性
解决毒性	原核心含结构警报或产生活性代谢物

几何约束的思路

ReCore 的做法与基于相似性的方法不同——它从三维几何出发：

• 1. 定义出口向量：在已知结合姿势中，标出核心与各取代基的连接点，记录这些连接键的位置与方向（向量）；
• 2. 搜索片段库：在预建的三维片段库中，寻找能提供相同出口向量几何的核心结构——即连接点的空间距离与角度都匹配；
• 3. 组装与评估：把新核心与原取代基组装，检查空间冲突与相互作用是否保留。

关键洞察是：真正重要的不是核心本身，而是它把取代基「摆在什么位置、什么朝向」。只要新核心能把关键取代基放到同样的空间位置，化学结构完全不同也可能保持活性。

成功的关键前提

• 必须知道正确的结合姿势：出口向量是从三维结构中读出来的。没有共晶结构或可靠的对接姿势，这套方法无从谈起。
• 要判断哪些取代基是关键的：与蛋白形成关键氢键、盐桥、占据关键疏水口袋的取代基必须保留位置；伸向溶剂的部分容差大。
• 核心本身可能也参与相互作用：如果原核心与蛋白有直接接触（如芳环的 π 堆积），换掉它风险大得多。
• 构象刚性影响大：新核心的柔性不同，会改变取代基的构象分布，即使静态几何匹配也可能活性下降。

开源实现思路

# 1. 用 RDKit 做基于 MCS 的核心替换
from rdkit import Chem
from rdkit.Chem import AllChem, rdRGroupDecomposition

# 分解出核心与 R 基团
core = Chem.MolFromSmarts("c1ccc2[nH]ccc2c1")   # 原核心
res, unmatched = rdRGroupDecomposition.RGroupDecompose([core], mols, asSmiles=True)
# res 中每个分子被拆成 Core + R1 + R2 + ...

# 2. 枚举新核心并重组
new_cores = ["c1ccc2occc2c1", "c1ccc2sccc2c1", "c1cnc2[nH]ccc2c1"]
for nc in new_cores:
    # 用带原子映射的 SMARTS 反应把 R 基团接到新核心上
    rxn = AllChem.ReactionFromSmarts(f"...")
    # 生成新分子

# 3. 用约束嵌入检查几何可行性
from rdkit.Chem import rdMolAlign
# 把新分子的关键取代基原子对齐到原分子的对应位置
# 计算 RMSD，判断几何是否匹配

# 更实用的开源路径：
#   - CReM：基于化学规则的分子生成/变换库，支持核心替换
#     pip install crem
#   - mmpdb：从数据中挖掘匹配分子对，找出真实发生过的核心替换
#   - SwissBioisostere / SwissSimilarity：在线等排与相似性检索
#   - REINVENT4 的 LinkInvent 模式（见 221）：给定两个片段设计连接部分

CReM 和 mmpdb 是最值得掌握的两个开源工具：前者做规则化的结构变换，后者从真实数据中挖掘「哪些替换实际发生过且保持了活性」——后者尤其有价值，因为它反映的是经验证的化学知识。

验证与风险

• 骨架跃迁的成功率不高：几何匹配只是必要条件，实际活性受构象熵、局部静电、去溶剂化等多重因素影响。做好「大部分尝试会失活」的心理准备。
• 务必做对接或 FEP 验证：候选新核心先用计算评估，再决定合成哪些。
• 合成可行性要提前确认：新核心可能极难合成或没有商业砌块。
• 专利判断交给律师：计算工具能提供化学思路，自由实施（FTO）的法律判断必须由专业人士完成（见 234《SureChEMBL》）。

延伸资源

• 同系列：303《BioSolveIT SeeSAR》、307《BioSolveIT infiniSee》；等排替换：281《Cresset Spark》；
• 骨架与 R 基团概念：040《Scaffold 骨架》、042《R-group Decomposition》。