1. 引言:当AI开始“画靶”,我们该如何审视“射箭”?
最近几年,AI驱动的分子生成领域可谓风生水起。从早期的基于规则和片段拼接,到如今基于深度生成模型(如VAE、GAN、流模型、扩散模型)的“从头设计”,我们似乎已经站在了药物发现的“奇点”边缘。打开任何一篇顶刊论文,你都能看到令人惊叹的指标:生成分子的新颖性高达99%,与靶标蛋白的结合亲和力预测值(如docking score)远超已知活性分子,ADMET(吸收、分布、代谢、排泄、毒性)性质预测一片飘绿。标题中提到的“基于靶标的从头分子生成”,正是这个领域的皇冠明珠——给定一个靶标蛋白的三维结构或序列,AI模型就能像一位不知疲倦的超级化学家,源源不断地“构想”出全新的、理论上能与之完美结合的候选药物分子。
这听起来无比美好,不是吗?但作为一名在这个交叉领域摸爬滚打了多年的从业者,我越来越频繁地听到一个灵魂拷问,它也正是这篇博文标题所隐喻的:我们引以为傲的“理性设计”,究竟是在进行一场严谨的科学探索,还是在玩一场“先射箭,后画靶”的游戏?
所谓“射箭画靶”,指的是先有了一个预设的、漂亮的“结果”(比如一个生成模型产出的、各项计算指标都完美的分子),然后我们再去为这个结果寻找合理性,甚至调整“靶心”(评价标准)来确保箭矢命中。在分子生成中,这可能表现为:我们精心优化模型,使其在某个特定的、有缺陷的评估体系下(比如单一的对接打分函数)产出“高分”分子,然后宣称我们实现了“理性设计”。然而,这些在虚拟世界中光芒万丈的分子,一旦进入湿实验验证,往往黯然失色,成功率远低于预期。
这篇博文,我想和你一起,暂时抛开那些令人眼花缭乱的模型架构和SOTA指标,沉下心来,从第一性原理和实践经验出发,系统地审视“基于靶标的从头分子生成”这一范式。我们将深入探讨其背后的技术逻辑、评估体系的潜在陷阱、以及从“虚拟分子”到“真实药物”之间那道看似很近、实则遥远的鸿沟。这不是要否定这项技术的价值,恰恰相反,是为了让它走得更稳、更远。
2. “理性设计”的技术基石:我们究竟在优化什么?
要理解“射箭画靶”的风险,首先得明白当前主流技术是如何“射箭”的。基于靶标的从头分子生成,其技术栈通常可以分解为三个核心环节:靶标表征、分子生成模型、以及驱动生成的优化目标(即“靶心”)。问题往往就埋藏在这第三个环节。
2.1 靶标表征:从静态结构到动态画像
绝大多数生成模型依赖于靶标蛋白的静态三维结构,通常来自蛋白质数据库(PDB)中的晶体结构或AlphaFold2等工具的预测结构。我们将这个结构转化为模型可读的输入,常见方法包括:
- 网格化(Grid-based):将结合口袋划分为三维网格,每个网格点用特征(如原子类型、电荷、疏水性)填充,形成一种3D图像。
- 图表示(Graph-based):将蛋白质视为图,残基或原子是节点,相互作用是边,并附上丰富的化学和物理特征。
- 表面描述符:计算结合口袋的几何形状、静电势、疏水斑块等。
注意:这里埋下了第一个隐患。蛋白质是动态的,存在构象变化(构象异构)。一个静态的“快照”能否代表其在生理状态下的真实形态?尤其是对于某些柔性较大的靶点,依赖单一静态结构进行设计,无异于“刻舟求剑”。
2.2 分子生成模型:从随机噪声到候选分子
这是AI大显身手的舞台。主流模型都在学习一个从随机噪声或潜在空间,到满足特定条件的分子结构的映射。近年来,扩散模型(Diffusion Model)因其强大的生成质量和稳定性,已成为该领域的首选。
- 前向过程(加噪):将一个真实的分子(或它的3D坐标与特征)逐步添加高斯噪声,最终变成纯随机噪声。
- 反向过程(去噪):模型学习从噪声中逐步恢复出分子结构。关键在于,这个去噪过程可以被“条件”所引导。
- 条件控制:我们将靶标口袋的特征作为条件输入到去噪模型中。这样,在生成过程的每一步,模型都在问自己:“在当前噪声状态下,如何向既符合化学规律、又能与旁边这个靶标口袋相匹配的分子结构去演变?”
这个过程本质上是在化学空间和几何空间中进行的一次极其复杂的、受控的采样。
2.3 优化目标:“靶心”是如何被定义的?
这是所有问题的焦点。我们如何告诉模型什么是“好”分子?通常,这个“靶心”由一系列计算预测指标构成,它们被组合成一个复合的奖励函数(Reward)或条件,用于训练或微调生成模型。
一个典型的“靶心”可能包括:
- 结合亲和力(Docking Score):通过分子对接软件(如AutoDock Vina, Glide)快速估算生成分子与靶标的结合强度。这是最核心、最常用的指标。
- 类药性(Drug-likeness):如Lipinski五规则(Ro5)、定量估算(QED)分数,确保分子具有成为口服药物的基本潜质。
- 合成可行性(Synthesizability):如SA Score(合成可及性分数),评估分子在实验室中被化学合成的难易程度。
- 新颖性(Novelty):与已知分子数据库(如ChEMBL)的相似度,鼓励模型发现全新骨架。
- 其他性质:如预测的溶解度、代谢稳定性、毒性等。
“射箭画靶”的苗头就此出现:我们整个生成系统的“理性”,完全依赖于这些计算预测指标的准确性和完整性。如果这些指标本身存在偏差,或者它们与最终的真实目标(体外/体内活性、成药性)关联性不强,那么无论我们的模型多么强大,优化过程多么精巧,都只是在完美地命中一个“错误的靶心”。
3. “画靶”的诱惑:评估体系中的陷阱与循环论证
让我们深入看看,在追求高分的过程中,我们是如何不知不觉地“画靶”的。
3.1 对接打分的“暴政”与过拟合
分子对接是虚拟筛选的基石,但其局限性在生成任务中被急剧放大。
- 打分函数的局限性:对接打分函数(如Vina, Glide SP)是经验性或基于力场的快速估算工具,并非精确的自由能计算。它们对某些相互作用(如氢键、π-π堆积)的权重可能失真,且无法很好地处理溶剂化效应、蛋白质柔性等复杂因素。
- 生成模型对打分函数的过拟合:这是一个致命问题。如果生成模型在训练或强化学习过程中,唯一或主要的优化目标就是降低某个特定对接程序的打分(比如Vina score),那么模型会迅速学会“欺骗”这个打分函数。它可能生成一些在Vina计算框架下能产生异常高分(低能量)的分子,但这些“优势”可能源于打分函数自身的漏洞或简化假设,而非真实的生物物理相互作用。
- 案例:模型可能学会生成一些具有特殊构象或官能团的分子,这些官能团恰好与打分函数中某个项的权重强烈耦合,从而获得高分,但这些基团在现实中可能不稳定、有毒性或根本不适合合成。
这就好比学生不是通过学习知识来取得好成绩,而是通过深入研究某位出题老师的个人偏好和历年真题的固定套路,来“优化”自己的答案。最终考试成绩(对接打分)很高,但真实能力(结合活性)未必强。
3.2 数据偏差与“内卷”式生成
生成模型的训练数据通常来源于已知的活性分子-靶标对(如PDBbind, BindingDB)。这些数据存在天然偏差:
- 化学空间覆盖不全:数据库中的分子多是已被合成和测试的,它们只代表了浩瀚化学空间中极小、且有偏的一部分(更偏向于“可合成”“已报道”)。
- 靶标偏好性:某些热门靶点(如激酶、GPCR)数据极多,而许多重要但难啃的靶点数据稀少。
在这种有偏的数据上训练,模型会倾向于生成与训练集分子“相似”的新分子,尽管它们可能具有新颖的SMILES字符串,但核心骨架和相互作用模式可能并未跳出原有窠臼。为了追求“新颖性”指标,模型又可能被推向生成一些极其怪异、合成路线未知的分子。这种在“模仿已知”和“追求怪异”之间的摇摆,并不是真正的理性探索。
3.3 多目标优化的“纸面平衡”
为了获得“完美”分子,我们会设置多目标优化:高结合力、高类药性、高合成可行性、高新颖性。在计算中,我们可以通过加权求和或帕累托前沿等方法,找到一系列“最优折衷”的分子。
然而,这种平衡是“纸面上”的。每个指标都是一个预测模型,都有误差。当这些带有误差的预测值被组合时,误差会被叠加甚至放大。一个在各项预测指标上都“均衡优秀”的分子,可能在实际中因为某一项预测的微小偏差(比如毒性被低估)而彻底失败。我们优化的是一个由不完美代理指标构成的“完美幻象”。
4. 从“虚拟靶心”到“真实命中”:缺失的验证环节
“射箭画靶”最直接的后果,就是虚拟世界与真实世界的脱节。一个在计算评估中全面胜出的分子,需要经历怎样的考验才能证明其价值?
4.1 湿实验验证:照妖镜与第一道坎
计算生成的分子必须进行化学合成与生物活性测试,这是无可替代的验证。
- 合成可行性陷阱:SA Score等指标仅是基于反应规则的粗略估计。一个SA Score尚可的分子,其合成路线可能长达十几步,总产率极低,成本高昂,根本不具备开发价值。我经历过多次,计算化学家兴高采烈地送来一个“完美”分子,有机合成同事看了直摇头,表示关键步骤没有已知方法,或需要用到极其昂贵、不稳定的试剂。
- 生物活性测试的复杂性:
- 测定方法:常见的测试包括酶活性抑制实验、细胞增殖抑制实验等。但计算预测的结合亲和力(dG)与实验测得的IC50/Ki值并非线性关系,还受到测试体系、底物浓度等多种因素影响。
- 假阳性:分子可能因在溶液中聚集、与测定体系中的组分发生非特异性结合、或自身发生化学降解而导致假阳性信号。
- 选择性:对目标靶点的高活性,可能伴随着对相关同源蛋白的交叉反应,导致脱靶毒性,这在计算阶段很难充分评估。
4.2 早期ADMET评估:被忽略的“隐形靶心”
即使一个分子在简单的生化水平上有活性,它要成为药物,还必须跨越ADMET这座大山。计算预测在这方面同样力有不逮:
- 渗透性与代谢:细胞渗透性(Caco-2, PAMPA预测)、代谢稳定性(肝微粒体稳定性预测)是早期淘汰率极高的环节。许多分子因为透膜性差或首过代谢太快而失败。当前的生成模型很少将这类复杂、噪声更大的预测模型深度整合到优化循环中。
- 毒性:预测药物毒性是计算化学中最具挑战性的任务之一。常见的致突变性(Ames)、心脏毒性(hERG)预测模型特异性虽有一定提升,但远未达到可靠的程度。依赖这些预测来“过滤”分子,风险极高。
真正的“理性设计”,应该从一开始就将这些关键的成药性因素作为“硬约束”或“核心优化目标”之一,而不是事后过滤的筛子。但目前的技术框架,往往更擅长优化一个清晰、易计算但可能失真的目标(如对接打分),而非一系列模糊、嘈杂但更真实的目标。
5. 走向更坚实的理性:技术反思与未来路径
认识到问题所在,是为了寻求解决之道。要避免“射箭画靶”,我们需要在技术路径和评估哲学上进行革新。
5.1 构建更接近真实的评估体系
- 采用更精确的计算方法:在关键环节,用更耗时但更准确的方法替代快速预测。例如,对于最终筛选出的少量顶级虚拟分子,使用分子动力学模拟(MD)来评估结合模式的稳定性、计算结合自由能(如MM/PBSA, MM/GBSA),这比静态对接能提供更丰富的动态信息和更可靠的能量估计。
- 集成多构象受体:不要只用一个静态晶体结构。可以采集来自MD模拟的多个代表性构象,或直接使用AlphaFold2预测的多条状态结构,让生成模型面对一个“靶标集合”,从而设计出对靶标微小构象变化不敏感的、更具鲁棒性的分子。
- 发展并审慎使用多参数优化(MPO):设计更智能的MPO策略,不仅考虑指标的加权和,更考虑指标的可靠性权重和阶段优先级。在早期生成阶段,可以更注重新颖性和合成可行性;在后期优化阶段,则需引入更可靠的物化性质和初步毒性预测。
5.2 改进生成模型本身
- 物理信息嵌入:将基本的物理原理和力场项作为模型训练的约束或正则化项,而不仅仅是事后的过滤条件。例如,在扩散模型的训练中,引入分子力场计算的能量项作为辅助损失,引导模型生成构象能更低的分子。
- 融合实验反馈的主动学习循环:这是打破虚拟循环的关键。建立“计算生成 -> 合成与测试 -> 数据反馈 -> 模型更新”的闭环。即使初期湿实验失败的数据(哪些分子没有活性、哪些难以合成),对于模型来说也是极其宝贵的负样本,可以用于修正模型对有偏数据集的依赖,使其学习到更接近真实世界的化学与生物学边界。
- 分阶段、分层次的生成策略:模仿药物化学家的思维,不要试图一步到位生成“完美”分子。可以先让模型生成满足核心药效团和初步结合模式的分子骨架(Scaffold),再通过后续的优化模型(如基于图的修饰模型)来逐步添加侧链、优化性质,每一步都有更聚焦、更可靠的评估。
5.3 调整社区的评价标准与文化
最后,也是最难的一点,是改变我们评价工作的“靶心”。
- 从“SOTA指标”到“实验验证”:一篇论文的价值,不应仅仅体现在它生成的分子在某个标准测试集上的对接打分比前人高了多少,或者新颖性达到了99%。社区应更看重工作是否包含了切实的合成与活性验证,哪怕只有几个分子。即使验证结果不理想,对失败原因的分析也具有极高的价值。
- 公开与共享:鼓励研究者公开其生成分子的详细列表、计算评估结果,以及最重要的——后续的实验测试数据(无论成功与否)。这能帮助整个社区更准确地评估不同生成方法和评估流程的真实有效性,避免大家在各自封闭的“画靶”体系中内卷。
- 拥抱不确定性:在报告中,不仅展示生成的“最佳”分子,也应坦诚地说明整个生成过程的不确定性。例如,展示分子在多次生成中的出现频率、其对打分函数微小变化的敏感性、以及不同评估指标之间的冲突情况。一个对扰动不敏感、各项指标一致的分子,通常比一个各项指标顶尖但极其脆弱的分子更可靠。
在我个人的项目经历中,曾耗费数月优化一个生成模型,使其在内部测试集上的对接打分和类药性分数达到极致。我们满怀信心地将排名前10的分子交给合作方合成,结果仅有1个被成功合成出来,且该分子在最初的细胞活性测试中毫无信号。这次挫折让我们彻底反思,将大量精力投入到了如何将合成路线预测(基于逆合成分析工具)和初步的细胞毒性预测更紧密地耦合到生成循环的早期阶段。虽然这增加了计算开销,也使得生成分子的“纸面数据”不再那么惊艳,但后续批次分子的合成成功率和初步活性筛选的“命中率”有了实实在在的提升。
生成式AI为药物发现带来了前所未有的可能性,但它不是点石成金的魔术。“基于靶标的从头分子生成”要想从一项炫技的研究,成长为一种可靠的药物设计工具,我们必须时刻保持警惕,避免陷入“射箭画靶”的自娱自乐。这要求我们不仅追求模型的“智能”,更要追求评估体系的“真实”和科学验证的“严谨”。理性设计的灵魂,不在于算法的复杂,而在于我们对化学、生物学复杂性的敬畏,以及将虚拟与现实反复对照、迭代修正的耐心。这条路很长,但每一步都值得踏稳。