在GROMACS中使用FMM静电学模拟恒定pH值. (B) GPU 加速汉密尔顿式插值

科学领域:

• 计算化学和生物物理学
• 分子动力学模拟的模拟.
• 生物分子建模模型

背景情况:

• 生物分子中可定位组的质子状态极大地影响其结构和功能.
• 标准的分子动力学 (MD) 模拟通常忽略了质子化状态和生物分子构造之间的动态合.
• 在恒定的pH值下模拟系统需要专门的方法,这些方法在计算上昂贵.

研究的目的:

• 引入一种严格而高效的方法来执行恒定pH分子动力学模拟.
• 克服现有的哈密尔顿插值 (HI) λ-动力学方法的计算局限性,特别是粒子网Ewald (PME).
• 为了能够准确有效地模拟大型生物分子系统中的质子化动态.

主要方法:

• 在GROMACS模拟套件中开发多极加速哈密尔顿插值 (MAHI) 方法.
• 使用预计算的电荷尺度哈密尔顿数与快速多极方法 (FMM) 进行有效的力计算.
• 将MAHI与电荷插值 (QI) 和标准HI方法进行比较.

主要成果:

• 与基于FMM的标准MD相比,MAHI实现了精确的恒定pH模拟,计算开销是可以忽略不计的.
• 哈密尔顿插值 (HI) 显示了更频繁的质子化状态转换,导致更好的采样和精确的电荷插值 (QI).
• 在一个百万原子系统上,模拟了多达512个可定位位点,运行时间的增加不到20%.

结论:

• MAHI为大规模恒定pH分子动力学模拟提供了一个计算可行的方法.
• 该方法可将恒定pH模拟无整合到标准力场工作流中.
• 这一进步使得对生物宏分子的pH依赖性行为进行更准确的研究成为可能.