layout: post title: 二维红外谱图的模拟计算及gmx_ir脚本更新 categories:
实验常用的谱图一般都有对应的二维谱, 如二维核磁, 二维红外等. 基于分子动力学模拟计算二维谱图并不容易, 涉及多方面的知识, 我涉足不多, 不敢妄议. 这里只将自己目前了解到的些许和想法随意放在这里, 或有所裨益, 或徒增误导, 自便.
如果我们聚焦在红外谱图的分子动力学模拟计算上, 那么可以试试最近出的一个工具AIM, 论文见 Kim E. Van Adrichem, Thomas L. C. Jansen; AIM: A Mapping Program for Infrared Spectroscopy of Proteins; J. Chem. Theory Comput. 18(5): 3089-3098, 2022. 可以算一维和二维的红外谱图, 而且给出了GROMACS和NAMD示例. 但用起来还是有点麻烦的, 所以我也就没有测试.
另一种获得二维谱图的方法是基于通用的相关光谱方法, 大致属于一种分析方法, 可用于任意的光谱, 任意的外界扰动, 具体理论和谱图解读可以参考:
可用的程序比较简单, 2DShige. 如果我们以时间作为参数, 计算不同时间段的红外, 或以温度为参数, 计算不同温度下的红外, 或者在非平衡模拟下计算不同时间段的红外, 都可以分析相应的二维相关图.
这里我以TIP3P柔性水模型来测试一下. 计算方法类似以前的设置, NPT平衡模拟, 步长2 fs, 总时间100 ps, 保存每步的坐标和速度. 模拟完成后使用gmx_ir.bsh脚本计算每10 ps的红外谱图, 共10组, 0-10, 10-20, 20-30, ⋯, 90-100.
因为是平衡模拟, 所以这10组谱图并没有太大差距.
将10组数据转换成csv格式(波数取到5000即可), 使用2DShige对这10组谱图进行二维相关分析
手动调整Zmax的值, 各种谱图如下:
说实话, 我没能力解读这些谱图所揭示的信息. 或许是因为平衡模拟的原因? 大致看来, 4000附近的OH振动与2000附近的剪式振动是正相关的, 而与1000附近的振动是负相关的? 或许, 用更复杂的分子, 用不同温度下的谱图进行分析能得到更容易解读的信息.
这种相关分析方法可用于任意的数据分析, 鉴于MD模拟可以提供很多时间相关的数据, 或许可以将其用于更多的MD性质分析, 如残基运动分析, 相互作用分析等. 感兴趣的可以试试.
gmx_ir脚本更新-b, -e选项选择不同的时间段, 类似GROMACS的相应选项.