gmx_mmpbsa脚本更新:清理整理输出

类别:    标签: gmx   阅读次数:   版权: (CC) BY-NC-SA

前段时间有人指出, 我的gmx_mmpbsa脚本的能量分解输出中的残基名称与原始文件中的不一样. 我查看了一下, 确实如此, 因为脚本中对所有残基进行了重新编号. 如果原始文件中的残基编号不是从1开始的顺序编号, 那自然与脚本所给的不一致. 这也不能算是个错误, 只是输出时的不同考虑. 不过考虑到一般是以原始文件的编号为准, 修正一下可能更方便使用.

借着修正这个问题的机会, 我顺便把所有的输出文件也重新清理或整理了一下, 以期更方便查看和使用. 这不是很重要, 却很繁琐. 因为这个脚本已经很长了, 重构起来都有点吃力, 接近快要失去控制的感觉, 生怕修改后的计算结果与以前的不一致. 不过, 好在完成了, 初步测试也没有发现结果错误.

以前的脚本说明文章暂时没有机会修正, 这里简要列出一些更新说明吧:

+==============================================================+
             dG = RTln(Ki) = -RTln(KA)
             Ki = 1/KA = exp(dG/RT) = IC50 = EC50
+==============================================================+
| M(mol/L) |  m/u/pM | dG(kcal/mol) | dG(kJ/mol) | Affinity    |
|==============================================================|
| 0.1      | 100  mM |    -1.364    |   -5.708   |             |
| 0.01     |  10  mM |    -2.728    |  -11.416   |             |
| 0.001    |   1  mM |    -4.093    |  -17.124   |             |
| 0.0001   | 100  uM |    -5.457    |  -22.832   | Weak        |
|--------------------------------------------------------------|
| 0.00001  |  10  uM |    -6.821    |  -28.540   |             |
| 1.0E-06  |   1  uM |    -8.185    |  -34.248   | Medium      |
|--------------------------------------------------------------|
| 1.0E-07  | 100  nM |    -9.550    |  -39.956   |             |
| 1.0E-08  |  10  nM |   -10.914    |  -45.664   |             |
| 1.0E-09  |   1  nM |   -12.278    |  -51.372   | Strong      |
|--------------------------------------------------------------|
| 1.0E-10  | 100  pM |   -13.642    |  -57.080   |             |
| 1.0E-11  |  10  pM |   -15.007    |  -62.788   |             |
| 1.0E-12  |   1  pM |   -16.371    |  -68.496   | Very strong |
+==============================================================+
◆本文地址: , 转载请注明◆
◆评论问题: https://jerkwin.herokuapp.com/category/3/博客, 欢迎留言◆


前一篇: 卌惑记
后一篇: 类似gnuplot的在线绘图工具

访问人次(2015年7月 9日起): | 最后更新: 2023-08-17 03:23:12 UTC | 版权所有 © 2008 - 2023 Jerkwin