[论文解读] Bayesian Estimation of the Hydroxyl Radical Diffusion Coefficient at Low Temperature and High Pressure from Atomistic Molecular Dynamics

Carter T. Butts, Rachel W. Martin|arXiv (Cornell University)|Aug 25, 2021

Spectroscopy and Quantum Chemical Studies参考文献 60被引用 7

一句话总结

本研究采用分层贝叶斯推断框架，利用TIP4P/2005水模型的原子尺度分子动力学（MD）模拟，在低温高压条件下估算水溶液中羟基自由基（·OH）的扩散系数。作者通过优化Lennard-Jones参数，开发出适用于CHARMM的·OH力场，使其在298 K和1 atm下重现实验测得的扩散系数0.23 Ų/ps，并基于贝叶斯后验推断构建了预测D·OH随温度和压力变化的对数-对数多项式模型，同时通过Yeh和Hummer的解析校正方法修正有限尺寸效应。

ABSTRACT

The hydroxyl radical is the primary reactive oxygen species produced by the radiolysis of water, and is a significant source of radiation damage to living organisms. Mobility of the hydroxyl radical at low temperatures and/or high pressures is hence a potentially important factor in determining the challenges facing psychrophilic and/or barophilic organisms in high-radiation environments (e.g., ice-interface or undersea environments in which radiative heating is a potential heat and energy source). Here, we estimate the diffusion coefficient for the hydroxyl radical in aqueous solution, using a hierarchical Bayesian model based on atomistic molecular dynamics trajectories in TIP4P/2005 water over a range of temperatures and pressures.

研究动机与目标

• 在实验测量不可行的低温高压条件下，估算水溶液中羟基自由基（·OH）的扩散系数。
• 开发一种适用于CHARMM的·OH力场，使其在TIP4P/2005水中准确重现298 K和1 atm下已知的0.23 Ų/ps扩散系数。
• 利用Yeh和Hummer的解析校正方法，修正分子动力学模拟中的有限尺寸效应。
• 基于贝叶斯后验推断，构建D·OH随温度和压力变化的预测对数-对数多项式模型。
• 通过量化极端条件下·OH的扩散行为，评估其在冰冻或深海环境中极端嗜极生物辐射损伤中的作用。

提出的方法

• 采用分层贝叶斯模型，从原子尺度MD轨迹中推断·OH扩散系数，整合多种模拟条件下的不确定性。
• 通过优化氧原子和氢原子的Lennard-Jones势阱深度（ε）与半径（rmin/2）参数，参数化基于CHARMM的·OH力场，使其在298 K和1 atm下匹配实验测得的D·OH = 0.23 Ų/ps。
• 采用准随机Halton序列对250组参数集进行采样，以探索εO、εH、rminO/2、rminH/2的四维参数空间。
• 每次模拟采用1 ns的生产运行，使用NpT系综，通过朗之万动力学控制温度，采用朗之万-Nosé-Hoover活塞控制压力，模拟介质为TIP4P/2005水。
• 应用Bullerjahn等人提出的基于协方差的估计器计算DMD，随后利用Yeh和Hummer的公式进行无限尺寸校正：D∞ = DMD + (ζkBT)/(2πηL)。
• 参数选择基于氢键作用能，优先选择O–H供体与O–O受体几何构型之间能量关系正确失稳的配置，以确保物理合理性。

实验结果

研究问题

• RQ1 在缺乏实验数据的低温高压条件下，水溶液中羟基自由基的扩散系数是多少？
• RQ2 如何参数化一种物理上准确的、适用于CHARMM的·OH力场，以重现常温常压下的已知扩散系数？
• RQ3 MD模拟中的有限尺寸效应在多大程度上会偏差估计的扩散系数？如何进行校正？
• RQ4 ·OH扩散系数如何随温度和压力变化？能否从模拟数据中推导出预测模型？
• RQ5 ·OH的迁移性在嗜冷菌和嗜压菌中辐射诱导生物大分子损伤过程中起什么作用？

主要发现

• 最优·OH力场参数为εO = -0.30 kcal/mol，rminO/2 = 1.5 Å，εH = -0.15 kcal/mol，rminH/2 = 1.0 Å，该参数组合在298 K和1 atm下准确重现实验测得的D·OH = 0.23 Ų/ps。
• Yeh和Hummer的校正使估计的扩散系数平均降低约10%，且在较小的模拟盒子中校正幅度最大。
• 最终的D·OH后验分布最适于用对数-对数多项式模型描述：log10(D·OH) = a + b·log10(T) + c·log10(P) + d·log10(T)·log10(P)，其中系数通过贝叶斯推断获得。
• 模型预测D·OH随温度降低和压力升高显著减小，与溶剂黏度增加和分子迁移率降低一致。
• 在D·OH估计中绝对相对误差（ARE）最小的参数集强烈偏好在1.3–1.8 Å范围内O–O排斥作用，表明空间位阻在准确模拟扩散中至关重要。
• 氢键分析表明，所选参数集正确地使O–O面对面构型相对于O–H供体构型能量失稳，确保了合理的溶剂化行为。