AMBER分子动力学模拟之结果分析(MMGB/PBSA)-- HIV蛋白酶-抑制剂复合物(4)
AMBER分子动力学模拟之结果分析(MMGB/PBSA)-- HIV蛋白酶-抑制剂复合物(4)
结合自由能计算
我们首先计算焙变,用到的是pbsa和gbsa方法。我们需要一下文件
三个top文件,pro.prmtop
lig.prmtop
com.prmtop
;输入文件MM_GBSA.in
;将要进行运算的轨迹文件md2.crd
;执行文件run.sh
;
MMGBSA
vim MM_GBSA.in
## 所用轨迹的参数
&general
startframe = 1, ## 帧的开始
endframe = 400, ## 帧的结束
interval=40, ## 帧间隔
use_sander =1, #
netcdf=1, ##轨迹是压缩格式吗
keep_files=0, ## 是否保存中间文件
debug_printlevel = 0,
verbose = 1,
entropy = 0 ## 是否计算熵
/
### 计算GBSA的参数
&gb
igb = 2, # GB类型
saltcon = 0, # 带电吗
ifant = 0, #
molsurf = 0,
surften = 0.005,
surfoff = 0
/
## 是否做残基分解
&decomp
idecomp=1,
dec_verbose=0
/
## 是否用nmode计算熵变
#&nmode
#nmstartframe =1,
#nmendframe =10,
#nminterval =1,
#dielc = 1,
#maxcyc = 500000,
#drms = 0.001,
#nmode_igb =0
#/
vim run.sh
单核计算
python MMPBSA.py -O -i MM_GBSA.in -o MM_GBSA.dat -eo MM_GBSA.csv -do MM_GBSA_DECOMP.dat -deo MM_GBSA_DECOMP.csv -cp com.top -rp pro.top -lp ../top/lig.top -y ../md/md2.crd> MM_PBGBSA.log并行计算
mpirun -np2 MMPBSA.py.MPI -O -i MM_GBSA.in -o MM_GBSA.dat -eo MM_GBSA.csv -do MM_GBSA_DECOMP.dat -deo MM_GBSA_DECOMP.csv -cp com.top -rp pro.top -lp ../top/lig.top -y ../md/md2.crd> MM_PBGBSA.log
参数说明
-np
2 2个core并行
-i
MM_GBSA.in input文件
-o
MM_GBSA.dat 结果文件
-eo
MM_GBSA.csv 详细结果文件
-do
MM_GBSA_DECOMP.dat 残基分解的结果文件(总的)
-deo
MM_GBSA_DECOMP.csv 残基分解的详细结果文件(每一帧)
-cp
com.top 复合物的 top
-rp
pro.top 文件受体的 top
-lp
lig.top 文件配体的 top
-y
…/md2/md2.crd 文件轨迹文件
MMPBSA
## 所用轨迹的参数
&general
startframe = 1, ## 帧的开始
endframe = 400, ## 帧的结束
interval=40, ## 帧间隔
use_sander =1, #
netcdf=1, ##轨迹是压缩格式吗
keep_files=0, ## 是否保存中间文件
debug_printlevel = 0,
verbose = 1,
entropy = 0 ## 是否计算熵
/
&pb
indi=1,
exdi=80.0
inp=1
cavity_offset=0.92,
scale=2.0,
istrng=0.1,
linit=1000,
prbrad=1.4,
radiopt=0
/&decomp
idecomp=1,
dec_verbose=0
/## 是否用nmode计算熵变
#&nmode
#nmstartframe =1,
#nmendframe =10,
#nminterval =1,
#dielc = 1,
#maxcyc = 500000,
#drms = 0.001,
#nmode_igb =0
#/
vim run.sh
MMPBSA.py -O -i MM_PBSA.in -o MM_PBSA.dat -eo MM_PBSA.csv -do MM_PBSA_DECOMP.dat -deo MM_PBSA_DECOMP.csv -cp com.top -rp pro.top -lp ../top/lig.top -y ../md/md2.crd> MM_PBGBSA.logmpirun -np2 MMPBSA.py.MPI -O -i MM_PBSA.in -o MM_GBSA.dat -eo MM_PBSA.csv -do MM_PBSA_DECOMP.dat -deo MM_PBSA_DECOMP.csv -cp com.top -rp pro.top -lp ../top/lig.top -y ../md/md2.crd> MM_PBGBSA.log
结果分析
运行结果MMGBSA.dat
,详细的每一帧的结果在MMGBSA.csv
残基分解的结果MM_GBSA_DECOMPdat
,详细的每一帧的结果在MM_GBSA_DECOMP.csv
AMBER分子动力学模拟之结果分析(MMGB/PBSA)-- HIV蛋白酶-抑制剂复合物(4)相关推荐
- AMBER分子动力学模拟之结果分析(构象分析)-- HIV蛋白酶-抑制剂复合物(6)
AMBER分子动力学模拟之结果分析(构象分析)-- HIV蛋白酶-抑制剂复合物(6) RMSD RMSF b-facto计算 RMSD RMSD measures the deviation of a ...
- AMBER分子动力学模拟之结果分析(突变型的能量计算,丙氨酸扫描)-- HIV蛋白酶-抑制剂复合物(5)
AMBER分子动力学模拟之结果分析(突变型的能量计算,丙氨酸扫描)-- HIV蛋白酶-抑制剂复合物(5) 丙氨酸扫描 在带电残基上引入一个或几个丙氨酸,观察这些改变对蛋白功能的影响.置换成丙氨酸,去除 ...
- AMBER分子动力学模拟之结果分析(最低能量结果)-- HIV蛋白酶-抑制剂复合物(3)
AMBER分子动力学模拟之结果分析(最低能量结果)-- HIV蛋白酶-抑制剂复合物(3) 在analysis目录下 解析.out文件 下载process_mdout.perl 脚本 perl proc ...
- AMBER分子动力学模拟之TOP准备-- HIV蛋白酶-抑制剂复合物(1)
AMBER分子动力学模拟之TOP准备-- HIV蛋白酶-抑制剂复合物(1) 我们以HIV蛋白酶-抑制剂复合物为例子,跑Amber动力学模拟 下载1phv 从PBD下载文件:https://www.rc ...
- CADD分子对接、薛定谔分子对接、AMBER分子动力学能量优化与分析、AIDD人工智能(机器学习与深度学习)药物发现
分子动力学模拟是分子模拟中最接近实验条件的模拟方法,能够从原子层面给出体系的微观演变过程,直观的展示实验现象发生的机理与规律,促使我们的研究向着更高效.更经济.更有预见性的方向发展.分子动力学可以解决 ...
- amber分子动力学模拟干货总结
AMBER软件安装运行.模型构建及力场文件生成.能量优化.溶剂模型.MM/PBSA方法结合自由能计算.另配合专业案例及经典文章案例复现 分子动力学入门理论 教学目标:了解本方向内容.理论基础.研究意义 ...
- gromacs manual_GROMACS蛋白配体分子动力学模拟结果分析简要笔记
0. 引言 本文以前文(https://zhuanlan.zhihu.com/p/149862369)为基础,对蛋白配体复合物分子模拟体系的结果进行一系列的粗浅分析,本文记述了简要的分析方法. 1 M ...
- Amber小分子-蛋白复合体分子动力学模拟
Amber小分子-蛋白复合体分子动力学模拟 以前经常用GROMACS进行分子动力学模拟,后来试了一下Amber后发现,在我当前配置的GPU资源上,果然还是Amber更快一些,GROMACS太吃CPU资 ...
- 分子动力学模拟手把手教你
如果你是AMBER的新用户 或对一般的分子动力学模拟毫无了解 可通过此教程入门. 教程简介 这个教程旨在介绍如何使用Amber进行分子动力学模拟,并假设您以前没有使用过Amber. 它专门为想要了解如 ...
最新文章
- ?类型 ,?类型(值)
- FastReport.net 使用记录
- 从excel表中生成批量SQL,将数据录入到数据库中
- 【OS学习笔记】九 实模式:从汇编的角度理解栈结构
- 算法(21)-leetcode-剑指offer5
- 红橙Darren视频笔记 IOC注解框架 自己写个注解框架
- jsp网页上实现计算圆面积小程序
- C语言目录文件操作补充
- [基于子串搜索的方法] BNDM算法
- tensorflow.python.framework.tensor_shape.is_fully_defined()
- 计算机在线给手机杀毒,手机在线杀毒,教您手机怎么杀毒
- 吉首大学期末计算机考试,吉首大学微机原理期末考试试卷.doc
- Bzoj 2563: 阿狸和桃子的游戏 题解
- 微信小游戏开发新手教程14-整合到一起,做出你的小游戏
- 快速下载/上传google drive文件的方法
- 【计算机网络】第三章 数据链路层的共同问题、点对点协议PPP
- SAP770系统FI模块配置(给科目表分配公司代码)
- DPDK 编译安装(meson ninja)
- 无极性电容的定义及应用
- Qt之QIODevice笔记