跳到主要内容
返回时间线
arXiv来源发表:

SO3LR Tauri 把力场扩到85种元素:一套权重同时跑蛋白、核酸、膜与含铂DNA

核心概要

SO3LR Tauri 把此前覆盖8种元素、400万结构的 SO3LR 力场扩展到85种元素、约1.48亿结构,并改写了短程排斥、静电与色散三个物理模块。它保留快速的 SO3krates 消息传递主干,在耦合簇参考上复现核碱基配对、蛋白–配体相互作用与糖构象排序,并在无体系专属拟合的动力学中给出水与有机溶剂密度、离子水合和多种生物分子结构。模型分 S/M/L 三档发布,单张 GPU 上可扩展到百万原子量级。

Source-provided article image: SO3LR Tauri: Generally Applicable Force Fields for Reactive (Bio)Molecular Simulations Covering the Periodic Table
arXiv

深度剖析

一个力场现在同时覆盖蛋白质、核酸、脂质、糖、水与有机溶剂、离子、含金属生物分子和反应路径,不再为每一类体系单独拟合参数。 此前的 SO3LR 只在8种元素、400万结构上训练,单卡可算20万原子;经典经验力场则按水、离子、金属中心等分别参数化,把几类生物分子放在一起模拟时需要拼接并非共同设计的参数。 在耦合簇参考的 SAPT10k 上,SO3LR-T-M 的相互作用能 MAE 为0.45 kcal/mol,是 SO3LR 的0.90 kcal/mol 的一半;在蛋白–配体活性位点集 PLA15 上为3.81 kcal/mol,相对误差2.1%。

模型把学习到的局部描述与三个显式物理项相加,长程静电和色散用粒子网格 Ewald 求和,因此在周期性大体系中不必截断。 SO3LR 的短程排斥是训练中拟合的 ZBL 型屏蔽函数,色散则用单一学习到的 Hirshfeld 体积比统一缩放自由原子系数,这种单标量缩放在整个周期表上过于受限。 短程排斥换成由专门 DFT 计算预表化的、按元素对区分的 NLH 势,网络改为逐原子预测极化率比与 C6 比;在300 K、1 atm 下四种独立重复的 NPT 模拟中,水与15种有机溶剂的密度平均相对误差为6.3%,OPLS4 为2.4%。

同一套预训练权重可以跑含铂的 DNA 交联和未修饰的 B-DNA 十二聚体,金属配位由势函数本身维持。 经典力场通常要为每种配位环境单独推导铂参数,并用固定配位几何的 improper 扭转项,配位数和电荷分布都无法随环境响应。 在顺铂 1,2-d(GpG) 交联的四个重复轨迹中,Pt(II) 在每一帧都保持四配位,没有配体交换,也未施加任何键合约束或 improper 扭转;Pt–N 键长与晶体学平均值的偏差在原文给出的阈值以内。

启示与展望

已展示的范围是静态基准、三类反应过渡态搜索,以及水、有机溶剂、盐溶液、肽 pKa、蛋白、糖蛋白、POPC 膜、铂交联 DNA 与 B-DNA 十二聚体的动力学模拟,全部在无体系专属拟合的前提下完成。作者同时说明,反应性研究中只有脯氨酸催化的羟醛缩合过渡态是盲搜得到的,安息香缩合与钴胺素甲基转移是从已发表结构出发精修的。

作者列出仍需推进的方向:静态与动态液体性质的预测精度、反应势垒、化学与构象空间覆盖,以及描述生物分子过程所需的采样与推理速度。文中还指出,有机液体密度存在系统性高估,部分可能来自核量子效应;水的自扩散约为实验值的一半,第一氧–氧峰仍偏高;肽 pKa 的 MAE 为0.98个单位,伯胺类基团被低估;B-DNA 的糖环构象在 O4′-endo 附近存在过量布居。