Posts
学习与思考
-
Python-Matplotlib做gmx_MMPBSA计算结果展示
Python, Numpy, Matplotlib, Pandas, Molecular Simulation ·学习使用GROMACS已经很久了，但是一直停留在很初级的应用上面，对高级技巧并不了解。就连做生物体系几乎必做的MMPBSA分析都没有尝试过。乘着这个假期了解了一下MMPBSA，才知道这个计算方法不仅仅是用于蛋白质-配体体系的结合自由能计算，对于任何一个二聚体都是可以的。另外，Jerkwin老师已经发展出了比之前常用的GMXMMPBSA与g_mmpbsa两种脚本/程序更加简单易用的计算脚本——gmx_mmpbsa。gmx_mmpbsa不但没有GROMACS与APBS程序的版本要求，还可以一步运行获得所有结果，大大降低了MMPBSA的学习成本。
Jerkwin老师的博客中有详尽的使用说明以及与其他两种方法计算结果的比对，另外，gmx_mmpbsa脚本的最新版本也在可以他的github库中找到。
gmx_mmpbsa使用前需安装GRMOMACS与APBS，脚本会自动调用这两个程序进行计算。Ubuntu环境下APBS可以使用
sudo apt install apbs直接进行安装。在修改脚本内的变量内容时，如果GROMACS与APBS都已添加进了环境变量，则可简写为：gmx='gmx'以及apbs='apbs'。脚本运行过程中如果出现某些
awk函数未定义的错误，那么还需要安装一下gawk，使用sudo apt install gawk即可。计算完成后会生成一系列不同结果的文档，这里编写了一个python脚本来进行绘图，顺便复习了一下Pandas与Matplotlib的使用方法。其中有关饼状图的修饰来自于Lemonbit的知乎专栏文章。
# This script is design to plot the mmpbsa calculate results # from Jerkwin's gmx_mmpbsa script # Author: Lewisbase # Date: 2020.02.29 import numpy as np import pandas as pd...