可极化AMOEBA模型用于模拟Mg2+蛋白核酸复合体

科学领域:

• 计算化学是一种计算化学.
• 生物物理学的生物物理.
• 结构生物学是结构生物学.

背景情况:

• 分子力学 (MM) 模拟为涉及核酸辅因子和双价的酶机制提供了详细的见解.
• 在MM模拟中模拟双价,特别是Mg2+,仍然是一个重大挑战,现有的极化力场在相互作用能量中显示出很大的错误,并且无法再现Mg2+·蛋白·ATP复合物的实验结构.

研究的目的:

• 系统评估和批判性修订可偏化的AMOEBA模型,以改进对Mg2+·蛋白·ATP复合物的模拟.
• 提高使用Mg2+和核酸辅因子的酶MM模拟的预测性能.

主要方法:

• 修订了使用高场校正 (AMOEBABIO18-HFC) 的AMOEBA蛋白模型,以改善Mg2+蛋白相互作用.
• 纳入多体非绑定固定 (NB-fix) 校正,以进一步减少相互作用能量错误.
• 开发一种新的AMOEBA模型,用于氨酸三酸盐 (ATP) 的修订极化,范德瓦尔斯 (vdW) 和二面参数.
• 与vdW-inclusive密度函数理论 (DFT) 和合集群 (CCSD(T)) 计算对比相互作用能量.
• 使用改进的模型,对Mg2+·Kinase·ATP复合体进行分子动力学 (MD) 模拟.

主要成果:

• 修订后的AMOEBABIO18-HFC模型显著降低了Mg2+蛋白相互作用能量误差 (MAE从17到10千卡/mol).
• 结合多体NB-fix纠正,MAE进一步降低至6kcal/mol (误差<2%).
• 新的ATP模型准确地预测了实验性Mg2+-ATP结合的自由能量,并提供了Mg2+关联的见解.
• 使用改进模型进行的MD模拟显示,模拟和实验确定的Mg2+·Kinase·ATP复合体结构 (X射线晶体学) 之间的一致性更好.

结论:

• 精细的AMOEBA模型,包括AMOEBABIO18-HFC和多体NB-fix校正,大大提高了模拟生物系统中Mg2+相互作用的准确性.
• 新的ATP模型增强了Mg2+-ATP结合的模拟.
• 这些进步使Mg2+依赖酶的MD模拟更加可靠,从而使得更好的结构预测与实验数据保持一致.