高等学校化学学报 ›› 2025, Vol. 46 ›› Issue (4): 20240518.doi: 10.7503/cjcu20240518

• 物理化学 • 上一篇    下一篇

基于键偶极模型的多肽-水体系多体极化作用强度高效计算新方法

祝佳怡, 姜笑楠, 郑笑函, 郝强(), 王长生()   

  1. 辽宁师范大学化学化工学院,大连 110629
  • 收稿日期:2024-11-25 出版日期:2025-04-10 发布日期:2024-12-31
  • 通讯作者: 郝强 E-mail:qh@luun.edu.cn;chwangcs@luun.edu.cn
  • 作者简介:王长生, 男, 博士, 教授, 主要从事理论与计算化学研究. E-mail: chwangcs@luun.edu.cn
  • 基金资助:
    国家自然科学基金(21773102);辽宁省教育厅基础研究项目(LJKMZ20221411)

A Novel Efficient Method Based on the Bond Dipole Model for Computing Many-body Polarization Strength in Peptide-Water Systems

ZHU Jiayi, JIANG Xiaonan, ZHENG Xiaohan, HAO Qiang(), WANG Changsheng()   

  1. School of Chemistry and Chemical Engineering,Liaoning Normal University,Dalian 116029,China
  • Received:2024-11-25 Online:2025-04-10 Published:2024-12-31
  • Contact: HAO Qiang E-mail:qh@luun.edu.cn;chwangcs@luun.edu.cn
  • Supported by:
    the National Natural Science Foundation of China(21773102);the General Research Program of the Department of Education of Liaoning Province, China(LJKMZ20221411)

摘要:

为解决传统非极化力场无法准确描述三体及更高阶多体作用的问题, 并精确模拟液态水环境下蛋白质的结构和功能, 提出了一种可以快速准确计算多肽-水体系中的多体极化作用强度的新方法. 该方法将多肽分子中的N—H和C=O极性化学键及水分子的O—H键视为键偶极, 周围环境对化学键的极化作用使得该化学键上产生了诱导键偶极, 使用键偶极之间的相互作用来描述体系中的多体极化作用. 通过拟合模型体系中三体作用能随分子间距离变化的势能曲线确定所需参数. 将该方法应用于多肽-水体系中的三体作用强度计算, 并与高精度的MP2方法和AMOEBABIO18可极化力场方法进行了系统比较. 结果表明, 本文方法计算的6个多肽水团簇(共包含92290个三体作用)的总三体作用强度与MP2方法的计算结果具有极高的线性相关性(线性相关系数为0.9965, 均方根误差为7.29 kJ/mol), 准确度略优于AMOEBABIO18可极化力场 方法 (其与MP2方法的线性相关系数为0.9950, 均方根误差为10.74 kJ/mol). 此外, 本文方法的计算效率与 AMOEBABIO18可极化力场方法相比也有显著提高, 在涉及超过20000次三体相互作用的模拟中, 计算时间缩短了约50%. 本文提出的高效方法为大规模蛋白质-水体系的多体作用模拟提供了新思路, 在相关领域中具有重要应用潜力.

关键词: 多体极化作用, 化学键偶极, 诱导键偶极, 多肽-水体系, 三体作用能

Abstract:

To address the limitations of traditional non-polarizable force fields in describing three-body and higher-order many-body interactions, and to accurately simulate protein structure and function in liquid water environments, we present a novel method for the efficient and accurate calculation of many-body polarization strengths in peptide- water systems. The N—H and C=O polar bonds in peptide molecules and the O—H bonds in water were treated as bond dipoles. The polarization effect of the environment on the chemical bond causes an induced bond dipole, and interactions of bond dipoles were employed to describe the many-body polarization interaction. The required parameters were determined by fitting the three-body interaction energy curves of the model molecules with different interaction distance. The method and parameters were applied to calculate the three-body polarization strength in peptide-water systems and the results were compared with those of high-precision MP2 method and the AMOEBABIO18 polarizable force field method. It is showed that the results of our method have a high linear correlation with those of MP2 method(correlation coefficient=0.9965, RMSE=7.29 kJ/mol) in the calculation of three-body interaction strengths in six peptide-water clusters(with a total of 92290 three-body interactions), which is superior to AMOEBABIO18 polarizable force field method(correlation coefficient=0.9950, RMSE=10.74 kJ/mol with the MP2 method). Moreover, the computational efficiency is significantly improved, with calculation time reduced by approximately 50% compared to AMOEBABIO18 polarizable force field method in simulations involving clusters with over 20000 three-body interactions. This efficient method offers a new approach for large-scale simulations of many-body interactions in protein-water systems, demonstrating significant potential for applications in related fields.

Key words: Many-body polarization interaction, Chemical bond dipole, Induced bond dipole, Peptide?water system, Three-body interaction energy

中图分类号: 

TrendMD: