利用分子模拟解释时间分辨实验的贝叶斯方法（A Bayesian Approach to Interpret Time-Resolved Experiments Using Molecular Simulations）

《ACS Physical Chemistry Au》：A Bayesian Approach to Interpret Time-Resolved Experiments Using Molecular Simulations

　　

时间分辨实验可提供蛋白质折叠、配体结合等动态过程的独特见解，但这些实验难以在分子层面解读，因其信号通常为多种构象态的平均，且所涉及过程常超出大多数模拟方法的可及时空尺度。本文描述一种利用分子模拟建模并解读时间分辨实验的方法——trBME（time-resolved Bayesian/maximum entropy，时间分辨贝叶斯/最大熵法），该方法将系统动力学模型与可与实验比对的三维蛋白结构模型相结合。研究人员利用合成的时间分辨小角X射线散射（trSAXS, time-resolved Small-Angle X-ray Scattering）数据演示该方法效用，展示其如何从实验测量中提取关于蛋白质去折叠过程的详细信息。

《利用分子模拟与时间分辨贝叶斯/最大熵法(trBME)解读时间分辨实验——以蛋白质去折叠trSAXS研究为例》论文解读

一、研究背景与意义

时间分辨实验（如时间分辨小角X射线散射 trSAXS、温度跃迁实验等）可在快速扰动（温度、压力、溶液条件改变）下捕捉生物大分子（如蛋白质）非平衡动态过程，有望提供"分子电影"般的关键生物学过程信息。传统分子模拟虽可研究生物分子动力学，但受限于可及时间与空间尺度及力场精度。现有的实验–模拟整合方法多为平衡（静态）实验数据的构象系综重加权(reweighting)，无法处理非平衡时间分辨数据；而直接对非平衡过程进行大量独立模拟或沿一维自由能面做朗文(Langevin)模拟受限于时间尺度，若模拟与实验不符亦无清晰策略用模拟解读实验；最大熵路径原则（Maximum Caliber, MaxCal）虽有应用但因需偏置模拟、实验数据点稀疏时效果差、复杂构象反应时间尺度难达等而难以实用。因此，亟需一种能在后处理（a posteriori）阶段整合时间分辨实验与分子模拟、兼容增强采样、不依赖直接模拟实验时间尺度的方法。本文由研究人员发表于《ACS Physical Chemistry Au》，提出 trBME（time-resolved Bayesian/maximum entropy）框架以解决上述问题。

二、主要关键技术方法

研究人员采用牛血清白蛋白（BSA）粗粒化C_α结构基模型与基于结构的势函数，先运行大量无偏高温短程模拟生成合成 trSAXS 数据（加散射角依赖噪声）；再在同温度下用阱偏好的 Well-Tempered Metadynamics (WTMetaD) 沿天然接触分数 Q（fraction of native contacts）为集合变量(CV, Collective Variable)采样获取自由能面 U(Q)=ΔG(Q)；假定扩散系数 D 为常数，以初始分布 p(Q₀, t₀)（折叠 basin高斯分布）与所得自由能面数值求解 Smoluchowski（Fokker–Planck）方程得一系列时间依赖概率密度即动力学前验(dynamical priors)；对每个实验时间点，以前述动力学前验为基础，用贝叶斯最大熵（BME, Bayesian Maximum Entropy）目标泛函（含约化 χ_R²项与相对熵 S_rel正则项，超参数 θ 控制先验信任度）逐时间点优化构象权重得后验(posterior)系综，并与均匀先验、平衡(WTMetaD偏差校正)先验对照验证。

三、研究结果

Introduction

现有重加权法只适用于平衡静态实验数据，对非平衡时间分辨实验不适用；已有替代方法（多轨迹平均、一维自由能面Langevin模拟、MaxCal偏置模拟）存在时间尺度受限、需新模拟适应不同实验、偏置实现困难等问题。研究人员指出时间分辨实验各时间点具准静态瞬时性，可用平衡方法分别处理，只需引入描述系统时间演化的动力学模型；且实验本身可提供非平衡动力学时间尺度，模拟只需提供构象系综，由此引出 trBME 思路。

Methods — Bayesian Maximum Entropy

标准 BME 对模拟轨迹 {X_i}_i=1…N赋予先验权重 w⁰_i（无偏 MD 为 1/N，增强采样可估算），计算 M 个静态可观测量 O_j^calc的加权平均 ?O_j^calc?(w)=Σ_iw_iO_j^calc(X_i)。通过最小化泛函 T(w|w⁰)= (1/M) Σ_j=1^M((O_j^exp? O_j^calc(w))/σ_j)²+ θ Σ_i=1^Nw_iln(w_i/w⁰_i) 获最优后验权重 w*，其中第一项为与实验偏差的约化 χ_R²，第二项为相对熵（负 Kullback–Leibler 散度）作正则防止过拟合，θ 为超参数，有效重采样比例 ?_eff=exp(S_rel) 表征重加权后保留的原模拟帧占比。该框架只对已采样构象重加权，未采样重要构象则效果不佳。

Methods — Dynamical Priors

时间分辨实验各时间点可视作局域平衡系综，需为每个时间点 t 提供时间依赖先验 w_t⁰。假设系统低维 CV 空间 {ξ_k} 中过阻尼 Langevin 动力学近似成立，对应条件概率密度 p(ξ,t|ξ₀,t₀) 满足 Smoluchowski 方程（Fokker–Planck 型）。由模拟得平均力势 U(ξ)（即沿 CV 的自由能面）及可选扩散系数 D(ξ)，本文进一步简化设 D(ξ)=D（常数），求解方程得从初态至平衡的一系列时间依赖概率密度，将之映射为模拟轨迹各帧权重即得动力学前验，并可依实验 χ_R²优选时间匹配解，解除 FP 方程起始时刻依赖。

Methods — trBME Pipeline

流程分三步：(1) 收集时间分辨实验数据，运行可描述过程且能给出沿选定 CV 自由能面的分子模拟（如 metadynamics、umbrella sampling、parallel tempering 等增强采样）；(2) 以模拟所得自由能面为 Smoluchowski 方程势，数值求解得时间依赖概率密度作动力学前验，依与实验吻合度选最佳匹配时间尺度；(3) 对各实验时间点独立用标准 BME 以对应动力学前验 w_t⁰重加权获时间分辨后验构象系综。

Methods — Temperature-Induced Unfolding of BSA

为验证 trBME，研究人员构建 BSA 热诱导去折叠的合成 trSAXS 数据集：用 C_α粗粒化结构基模型运行 1000 条 50 ns 无偏高温模拟获去折叠轨迹，计算各帧 SAXS 强度并系综平均加噪声得"实验" trSAXS 曲线；同条件下运行 500 ns WTMetaD（偏置 Q）获双阱（折叠态 Q≈0.9、中间态 Q≈0.65–0.83、去折叠态 Q≈0.5）自由能面供 trBME 输入。

Results and Discussion

对合成 trSAXS 数据应用 trBME：WTMetaD 给出多极小自由能面，解 Smoluchowski 方程得从折叠 basin向去折叠迁移并瞬态布居中件中间态的概率密度（动力学前验）。逐时间点 BME 重加权后，模拟 SAXS 曲线与合成实验高度吻合（平均 χ_R²=1.27±0.02，平均 ?_eff=0.913±0.004），表明动力学前验已接近真实且仅需微小幅值重加权。重加权后系综的回转半径 R_g、天然接触分数 Q、RMSD 分布均与生成合成数据的无偏 MD 真值相符。对比试验显示：用平衡 WTMetaD 先验（w_eq⁰）或均匀先验（w_uniform⁰=1/N）经 BME 虽可拟合 trSAXS（χ_R²分别 2.8±0.2 与 0.96±0.01）但构象分布明显偏离真值，所需重加权幅度更大（?_eff分别 0.31±0.04 与 0.77±0.02），折叠/中间/去折叠态布居随时间演化亦不如动力学前验准确（Jensen–Shannon 散度更大）。证明引入恰当动力学前严于时间分辨 BME 重加权对恢复正确非平衡构象演化至关重要。

四、结论（Conclusions 部分翻译）

研究人员描述了 trBME（时间分辨贝叶斯/最大熵）细化方法——一种利用分子模拟解读时间分辨实验数据的途径。trBME 基于两点：(i) 许多时间分辨实验在时间局部是静态的；(ii) 可利用实验获取非平衡动力学的时间尺度，从而用模拟解释此动力学背后的结构系综。这与以往常试图用模拟预测体系动力学的研究不同——研究人员主张实验通常较易获得整体动力学与时间尺度，故将重点放在用模拟为时间分辨实验提供结构阐释。研究人员用合成 trSAXS 实验演示 trBME 并在多种构象性质上与真值高度吻合，如能捕获不同构象态布居的时间演化。trBME 具普适性，可用于其他时间分辨实验及其他大分子的动态过程；虽本文用 metadynamics 捕获实验所探运动，trBME 亦兼容其他采样方法（如马尔可夫态模型 MSM 描述底层动力学），有望为理解生物功能相关的生物大分子动态运动提供新见解。