量子加速药物虚拟筛选
前情提要:第 13 篇《应用落地:密码、材料、药物、AI》把"药物"列为量子四大应用之一,只说了"长期看有望加速筛分子";第 30 篇《量子模拟:算分子,最近的落脚点》则深入讲了怎么用量子算单个分子的基态能量。本篇把镜头拉到更大的战场——一个装了几百万分子的"化合物库",看量子能在"挑分子"这一步改进什么。(记住:30 算的是单分子,本篇算的是一整个库。)
新药研发为什么这么贵、这么慢?一个被低估的原因是:在真正动手合成之前,要先从海量分子里挑出值得试的那几个。这一步一旦做砸,后面几亿美元的投入可能全打水漂。
一、药物发现的第一步:从海量库里挑分子
化工厂商手里通常有上百万、甚至上千万个已知分子(称作"化合物库"),公开的 ChEMBL、BindingDB 等数据库已收录上千万条分子与靶点的活性记录。要找能结合某个疾病靶点的候选药,第一步不是设计,而是海选:把每个分子"塞"进靶蛋白的结合口袋,算它卡得紧不紧——这叫分子对接;再用自由能计算给结合强度打分。谁分高,谁进下一轮实验。这个"海选"本身就是一场算力赌博:一个中等规模的虚拟库动辄上百万分子,逐一认真算一遍自由能,可能要烧掉数月机时。
二、库太大,算不动
库大到这种量级,逐一对每个分子做高精度自由能计算,算力根本扛不住。所以经典流程只能先粗筛(快但糙)再精算(准但慢):用便宜的经验规则先砍掉 99%,只对剩下 1% 认真算。
代价很现实——真正好的分子,可能在第一轮粗筛里就被粗糙规则误杀,永远没机会进精算。因为经验规则"用近似换速度",对结构新颖、不在训练分布里的分子,打分往往会系统性偏倚,于是"库太大"直接造成了漏网之鱼。
三、让量子先帮你"排个序"
量子怎么帮?我们第 38 篇的 QAOA、第 28 篇的变分框架,恰好擅长"在超大组合空间里快速找到高分区"。把"从库里挑高潜力候选"建模成组合优化或采样问题,用量子在几百万分子里快速把最有戏的一小撮顶上来,而不是平均地、昂贵地精算每一个。这和第 38 篇物流优化是同一类"在巨大空间里找高分区"的问题,只是目标从"最短路径"换成了"最高结合潜力"。对药企来说,这一步省下的不是几个 CPU 小时,而是可能少走几年临床弯路,这也是为什么 QAOA 这类组合优化算法在药物信息学圈被反复研究。
注意分工:单个分子的精确能量,仍交给第 30 篇那种量子模拟(VQE)去算;本篇的量子,只负责"先排序、先筛选"。
四、一条库到候选的加速链
把"经典逐一大库、量子先顶高潜力"的链路画出来:
五、加速的是排序,不是精度
这里必须泼冷水、也划清边界:量子在药物里加速的是"排序和筛选",不是"把单个分子算得更准"——后者是 30 的活。今天 NISQ 能做的是混合式的"候选排序":量子先把高潜力分子顶上来,经典再精算确认。
至于"用量子把一千万分子的库全自动筛完",那要等纠错成熟。已有团队用变分量子算法在很小的分子子集上演示过"候选排序",虽然离真实药库还远,但验证了路径可行。所以务实的判断是:量子在药物发现里的第一个价值,是少漏掉好候选、少做无用功,而不是一夜之间算出新药。这也解释了为什么产业界对量子药物更务实的预期是"辅助排序工具",而非"自动发现机器"。
费曼自检
合上页面,用三句话复述今天的主线。参考版本:
1. 药物发现先要从百万级化合物库里挑候选,经典靠"分子对接 + 自由能"逐一打分。 2. 库太大只能先粗筛再精算,好分子可能在粗筛被误杀;量子把"筛选排序"建模成组合优化/采样,先顶高潜力候选。 3. 量子加速的是排序而非单分子精度(那是 30 的活),NISQ 先做混合候选排序,全库自动筛选还远。
两道自测题(点开看答案):
<details> <summary>自测一:本文讲的"量子加速药物"和第 30 篇有什么边界?</summary>
答案:第 30 篇是用量子模拟算单个分子的基态能量(求"这个分子多好");本篇是用量子在百万级化合物库里做筛选排序(求"哪些分子值得先看")。一个是单点精度,一个是库级排序,分工不同。 </details>
<details> <summary>自测二:为什么量子不能"直接把整个药库算完"?</summary>
答案:因为对库里每个分子做高精度自由能计算本身极贵,NISQ 机器既装不下、也扛不住这种规模;更现实的是量子先顶出少数高潜力候选,再由经典精算确认,全库自动化筛选要等容错机。 </details>
一件动手小事:找一个公开的分子对接入门数据集(如 BindingDB 的小子集或教程示例),思考如果把"从库里挑前 1%"建模成组合优化,第 38 篇的 QAOA 会在流程的哪一步介入。
下一篇,我们把同样的"量子补经典短板"思路,从药物转到材料:电池电解质和催化剂,该怎么选。下次见。