量子-经典异构计算驱动的生物分子模拟研究进展:从千原子到万原子尺度的突破

薛定谔了么
2026-06-04 19:57:44
生命科学
量子信息
行业动态

本文介绍2026年5月克利夫兰诊所、日本理化学研究所与IBM团队的量子计算生物医药应用突破。研究依托量子中心超级计算框架,结合经典超算与156量子比特Heron r2处理器,采用新型EWF-TrimSQD混合算法,首次完成含12635个原子的蛋白质-配体复合物模拟,规模达半年前同类研究的40倍,关键步骤精度最高提升210倍,片段能量计算匹配耦合簇(CCSD)黄金标准。该成果突破传统量子硬件仅能处理百原子级体系的局限,将模拟对象拓展至真实生物场景下的蛋白复合物,验证了量子计算支撑药物早期发现中分子互作预测、结合能计算核心需求的可行性,标志着量子技术从原理验证阶段迈入实用科研工具阶段,为缩短药物研发周期、解析复杂生物分子机制提供了全新路径。




一、研究背景与问题提出 


药物研发中候选分子与蛋白质的结合机制解析是生命科学领域的核心难题,也是制约新药开发效率的关键瓶颈——当前单款药物研发周期普遍超过10年、成本高企,核心痛点在于传统计算方法随分子体系规模扩大呈指数级算力增长,难以对上万原子级别的蛋白质-配体复合物实现高精度从头算(ab initio)模拟。量子计算理论上可突破经典计算的指数墙,但长期以来受限于硬件规模与噪声问题,仅能完成百原子级的小分子模拟,尚未触及真实药物研发所需的生物体系尺度。


20265月,克利夫兰诊所、日本理化学研究所(RIKEN)与IBM联合团队在arXiv发布预印本成果,首次将量子硬件可模拟的具生物学意义的分子体系规模推升至12,635个原子,标志着量子计算正式逼近真实药物研发场景需求。



 


二、核心研究方法 


本研究采用量子中心超级计算(Quantum-Centric Supercomputing异构框架,结合算法创新与跨设施算力调度,具体技术路径如下:



  1. 任务拆分策略:依托经典超级计算机(日本富岳Miyabi-G)采用量子嵌入法,将完整蛋白质-配体复合物拆解为多个可独立计算的量子化学片段,规避全体系直接模拟的算力爆炸问题。

  2. 量子-经典分工执行:电子结构计算等量子优势任务分配至部署于克利夫兰诊所、RIKENIBM 156量子比特Heron r2处理器,单次最多调用94个量子比特,累计运行9200余个量子线路、耗时超100小时;经典超算同步负责片段构型采样、子空间对角化优化与最终结果整合,分布式线性代数内核并行效率达72.5%

  3. 算法优化:团队提出新型混合算法EWF-TrimSQD,通过局域子空间截断大幅降低量子资源开销,相较此前同类方法计算效率显著提升。



 图1 研究的系统(A)采用的总体算法方案(B)片段的表示(C)以及本工作中提出的低缩放片段构建的细节(D)


 


三、关键实验结果 


本次研究选取两种典型作用模式的蛋白质-配体复合物作为测试对象:分散作用主导的T4溶菌酶体系(11,608原子)、静电作用主导的胰蛋白酶体系(12,635原子),均置于水溶液环境中模拟,更贴近生理真实场景,核心指标突破包括:



四、研究价值与领域影响


应用层面


该成果验证了量子计算支撑药物早期发现两大核心需求的可行性:一是大规模生物体系的原子级动力学模拟,二是高精度结合自由能计算,未来可进一步扩展至酶催化机制解析、难成药靶点相互作用预测等场景,有望从源头缩短药物研发周期、降低试错成本。


学科层面


此次突破标志着量子计算的发展范式发生根本转变:行业评估指标从单一的量子比特数、门保真度、错误率,转向可解决实际科学问题的规模与价值,量子硬件正式从原理验证工具升级为可产出具科研意义结果的实用化装置。


五、研究基础与后续方向


该工作是三方长期合作的延续:此前团队已在《Science Advances》发表硫化铁分子电子态模拟成果,并完成首个20氨基酸全量子中心模拟体系(Trp-笼,303原子)。后续研究将聚焦于进一步优化量子纠错方案、提升异构工作流的自动化程度,目标实现更大规模生物分子体系的动态模拟,最终将量子计算整合为常规药物研发算力栈的一部分。




图2 实验:大基组下的化学计算




参考文献


[1] Merz K M, Shajan A, Kaliakin D, et al. Crossing the 12,000-atom barrier with heterogeneous quantum-classical supercomputing: quantum chemistry of protein-ligand complexes[J]. arXiv preprint arXiv:2605.0XXXX, 2026.


[2] IBM News Room. Cleveland Clinic, RIKEN, and IBM Model a 12,635-Atom Protein, the Largest Kno


65
1
0
0
关于作者
相关文章
  • npj Comp. Mater. | 统一扩散框架加速无机晶体材料设计:DiffCry ...
    传统无机晶体生成扩散模型多采用晶格、原子种类与坐标分步扩散架构,存在参数量庞大、采样效率偏 ...
    了解详情 
  • Cell: 当 Ising 模型遇见 Perturb-seq:用统计物理反推百万细胞 ...
    本文介绍加州理工学院与洛克菲勒大学团队发表于《Cell》的研究——D-SPIN(单细胞扰动 ...
    了解详情 
领取成功
本月5个550bit真机配额已发放给您,配额将在2个月后到期,请及时使用哦~
活动中心
联系我们
二维码
返回顶部
返回
活动中心

完成任务,轻松获取真机配额

×
每日必做
新手任务
长期任务
其他任务
快速回复 返回顶部 返回列表
玻色有奖小调研
填写问卷,将免费赠送您1个1000bit真机配额
(单选) 您是从哪个渠道得知我们的?*
您是从哪个社交媒体得知我们的?*
您是通过哪个学校的校园宣讲得知我们的呢?
取消

提交成功

真机配额已发放到您的账户,可前往【云平台】查看

量子AI开发者认证

考核目标

开发者能够成功搭建Kaiwu-PyTorch-Plugin项目基础环境,并成功运行示例代码,根据示例提示,输出指定的值并填写至相应的输入框中。

通过奖励

5个一年效期的1000量子比特真机配额

专属「量子AI开发者」社区认证标识

开发者权益

每月固定权益:5个550量子比特真机配额
前往考核

第一步

按照README提示成功安装Kaiwu-PyTorch-Plugin库环境依赖
前往GitHub

第二步

运行 community-assessment 分支下的 run_rbm.py 代码示例

第三步

理解示例代码,手动打印并填写如下数值:

正相采样的状态

负相采样的状态

正相的能量值

负相的能量值

*

提交答案

开发者权益

每月固定权益:5个550量子比特的真机配额

恭喜您完成考核

您将获得量子AI开发者认证标识及考核奖励

1000 bit*5

配额

Quantum AI Developer Certification

Assessment Objectives

Developers should successfully set up the basic environment for the Kaiwu-PyTorch-Plugin project, run the QBM-VAE sample code, and calculate the correct FID value based on the random seed value provided by the system.

Pass Rewards

10 quotas for 550-qubit real quantum machines with a one-year validity period

Exclusive "Quantum AI Developer" Community Certification Badge

Developer Benefits

Fixed Monthly Benefits: 5 quotas for 550-qubit real quantum machines
Proceed to Assessment

Step 1

Install the environment dependencies for the Kaiwu-PyTorch-Plugin library according to the README instructions
Go to GitHub

Step 2

Replace the Seed Value

Your seed value is

Step 3

Enter the FID Value You Calculated

*

Submit Answer

Developer Benefits

Fixed Monthly Benefits: 5 quotas of 550-qubit real machines

Congratulations on Completing the Assessment

You will receive the Quantum AI Developer Certification Badge and Assessment Rewards

550bit*10

Quotas