AFFINIKEY AI / 数据探索
第一阶段 · 公开实验数据 · 2026.10

小分子结构,能否帮 EGFR 和 ALK 活性排个序?

我们整理公开实验中的小分子结构和 IC₅₀ 数据,先检查模型能否预测靶点蛋白的酶活性。以下是目前的探索结果。

内部研究摘要 · 回顾性评估
初步看到预测信号,但还不能证明它能稳定挑出有效新分子。

结果来自一批既有公开数据和一次固定拆分;还需要更严格的独立评估及实验验证。

EGFR 模型训练 + 验证2,657条实验测量记录
EGFR 固定测试集241种分子结构 · 285 条测量
ALK 骨架测试集84种分子结构
候选库已整理,尚未打分≈5 万每个靶点约 5 万种分子

这项研究具体在做什么?

先回答一个窄问题:只看小分子的二维结构,能否估计它在某种 EGFR 或 ALK 酶实验中的活性?

IC₅₀ 不一定是在测癌细胞生长

本次主数据是靶点酶活性实验。它测的是小分子抑制 EGFR 或 ALK 酶活性的浓度,不是癌细胞分裂速度,也不能直接代表抗癌疗效。

测试里的“241”是结构数

同一种分子可能有多条实验测量。因此 EGFR 测试集有 241 种不同结构、285 条测量记录;两者不是同一个计数。

EGFR:模型误差有多大?

按化学骨架拆分:测试集中的分子骨架没有出现在训练集。每种方法都在同一批测试数据上评估。

平均预测误差(MAE)

单位:pIC₅₀ · 越短越好 · 条形按 0–1.2 刻度
随机森林
0.560
DeepSeek V4 Pro
0.605
DeepSeek Flash
0.617
均值基线
1.145

怎么读:MAE 是预测与实验结果的平均差距,越低越好。随机森林误差小于简单均值基线;它与两个 DeepSeek 结果的置信区间有重叠,目前无法判断谁明确更好。

预测排序与实验排序有多一致?

Spearman ρ · 越接近 1 越一致
随机森林
0.836
DeepSeek Pro
0.816
DeepSeek Flash
0.807

三个数值相近。它们表示同一测试集上的排序关联,不等于候选在真实实验里会命中。

方法结构平衡 MAE ↓95% 区间排序相关 ρ ↑
传统模型 随机森林0.5600.491–0.6300.836
DeepSeek V4 Pro0.6050.525–0.6880.816
DeepSeek Flash0.6170.543–0.6950.807
均值基线1.145——

ρ(Spearman 排序相关)越接近 1,预测排序越接近实验排序。这里按 285 条测量记录计算;一分子的重复测量会分别计入。MAE 和 ρ 是不同角度的指标。

多加特征、换算法后,有提升吗?

EGFR 同一骨架拆分下,结果没有支持更复杂方案稳定胜出。

尝试测试集结果目前能说什么
增加分子结构特征MAE 0.565原始基线为 0.560;未改善。
搜索传统模型与参数126 组;最低测试 MAE 0.567仍高于原始随机森林的 0.560。
DeepSeek Pro / FlashMAE 0.605 / 0.617差异不确定;不能判定 Pro 明确更好。

特征组合的测试区间为 0.500–0.640,与基线区间 0.492–0.631 大幅重叠。这里是探索性结果,模型选择并非一次全新的独立确认。

数据覆盖到哪里?

目前先用 ChEMBL 实验标签做建模。其他来源用于补充核对,不是把所有数据库的数字简单相加。

EGFR2,942 条

IC₅₀ 测量 · 2,358 种结构
随机森林测试 MAE:0.560

ALK877 条

IC₅₀ 测量 · 801 种结构
随机森林测试 MAE:0.577

EGFR 候选库50,289

种候选结构已整理,尚未批量预测

ALK 候选库52,300

种候选结构已整理,尚未批量预测

“未注明变体”表示公开记录没有写清具体蛋白变化,不等于已确认的野生型或同一种蛋白构建体。

三个词,快速看懂

IC₅₀

在特定实验中,让测量到的作用降低一半所需的分子浓度。本次测的是酶活性。

pIC₅₀

把 IC₅₀ 换到对数尺度,方便建模;通常数值越高,IC₅₀ 越低。这里的误差不是 nM 浓度差。

骨架拆分

把结构骨架不同的分子放到测试集,检查模型面对新形状时表现如何。

结论边界与下一步

这是研究流程检查,不是药物发现结果。

尚未证明跨实验稳定

不同实验条件可能影响 IC₅₀。来源拆分下 EGFR 排序相关接近 0,assay 还需要人工审查。

没有细胞或临床结论

模型没有预测肿瘤缩小、患者获益或治疗效果;也还没有对候选库批量打分。

大模型有记忆风险

DeepSeek 可能在训练中见过公开分子数据;这次回顾性测试无法排除这种影响。

建议先复核数据

人工检查实验类型和蛋白版本,再锁定评估方案;最后用新实验验证候选排序是否有用。