24小时热门版块排行榜    

查看: 1773  |  回复: 2

MolAICal

新虫 (初入文坛)

[交流] 使用MolAICal基于NAMD模拟结果计算小分子和蛋白MM/GBSA的教程 已有1人参与

使用MolAICal基于NAMD模拟结果计算小分子和蛋白MM/GBSA的教程


更多教程(含英文教程)请见如下:
MolAICal官方主页:https://molaical.github.io
MolAICal 文章介绍:https://doi.org/10.1093/bib/bbaa161
MolAICal中文博客:https://molaical.github.io/cntutorial.html
MolAICal blogspot:https://qblab.blogspot.com


1. 简介
在本教程中介绍了基于NAMD的分子动力学模拟结果,使用MolAICal计算小分子和Mpro蛋白受体MM/GBSA的方法。本教程只是一个简单演示。为了节省运行及存储空间,本教程仅选择了Mpro复合物分子动力学模拟的25帧用于计算。

2.工具
2.1. 所需软件下载地址
1) MolAICal : https://molaical.github.io
2) NAMD: https://www.ks.uiuc.edu/Research/namd/

2.2. 操作示例文件
所有用到的操作教程文件均可在下面的网站下载:
https://github.com/MolAICal/tutorials/tree/master/004-MMGBSA

3. 操作流程
转到以下目录:
#> cd 004-MMGBSA

3.1. 提取蛋白与配体复合物的轨迹文件
CODE:
#> vmd -dispdev text -psf "mpro.psf" -e stripDCD.vmd -args protein,or,resname,LIG "mpro.dcd"  "complex" mpro.psf mpro.pdb

-args: 其用法类似VMD软件中的“atomselect”命令,比如"atomselect top protein or resname LIG",此处逗号"," 代表空格。其中脚本文件stripDCD.vmd可以在本教程材料或MolAICal软件的“scripts”目录里面找到。

执行上述命令后生成complex.psf, complex.pdb 和 complex.dcd文件。将“GBIS” 和“sasa”参数设置为on。打开并按照下文内容修改“complex.conf”文件:
--------------------------------------------------------------------------
structure          complex.psf
coordinates        complex.pdb
outputName         complex

paraTypeCharmm      on
parameters          par_all36_prot.prm
parameters          par_all36_cgenff.prm
parameters          ligand.str
parameters          toppar_water_ions.str

coorfile open dcd complex.dcd
--------------------------------------------------------------------------
本教程中命令在CPU上运行。你可以选择GPU进行运算。在Linux系统下运行NAMD命令,如下:
CODE:
#> namd2 +p3 complex.conf >& complex.log &

其中符号“&”代表程序在Linux系统中进行后台运行,如果你使用的是Windows操作系统,请不要用“&”,例如,命令换成这样:
CODE:
#> namd2 +p3 complex.conf > complex.log

3.2. 仅提取蛋白的轨迹文件:
CODE:
#> vmd -dispdev text -psf "mpro.psf" -e stripDCD.vmd -args protein "mpro.dcd" "protein" mpro.psf mpro.pdb

上述命令会生成 protein.psf, protein.pdb和protein.dcd。打开“protein.conf”,参考 “complex.conf”修改相关参数。

本教程中命令在CPU上运行。你可以选择GPU进行运算。在Linux系统下运行NAMD命令,如下:
CODE:
#> namd2 +p3 protein.conf >& protein.log &

3.3. 仅提取配体的轨迹文件:
CODE:
#> vmd -dispdev text -psf "mpro.psf" -e stripDCD.vmd -args resname,LIG "mpro.dcd" "ligand" mpro.psf mpro.pdb

上述命令会生成 ligand.psf, ligand.pdb和ligand.dcd。打开“ligand.conf”,参考 “complex.conf”修改相关参数。

本教程中命令在CPU上运行。你可以选择GPU进行运算。在Linux系统下运行NAMD命令,如下:
CODE:
#> namd2 +p3 ligand.conf >& ligand.log &   

4. 用MolAICal计算MM/GBSA
CODE:
#> molaical.exe -mmgbsa -c complex.log -r protein.log -l ligand.log

输出结果中给出下文所示的结合自由能△G:
--------------------------------------------------------------------------
delta E(internal): -4.0000007572871255E-6
delta E(electrostatic) + deltaG(sol): 7.702936000001536
delta E(VDW) + deltaG(sol): -44.43611599999989
delta G binding: -36.73318399999911
--------------------------------------------------------------------------
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zxhuizhang

木虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
我想问一下,关于linux系统的电脑配置问题。大概需要什么样子的
见此
2楼2021-02-02 12:12:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

MolAICal

新虫 (初入文坛)

引用回帖:
2楼: originally posted by zxhuizhang at 2021-02-02 12:12:09
我想问一下,关于linux系统的电脑配置问题。大概需要什么样子的

如果已经跑完md, 仅仅用molaical算, 一般普通家用的linux系统就行。。当然用服务器更好了
3楼2021-04-27 12:26:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 MolAICal 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 环境调剂 +3 晓看天暮看云 2026-03-09 3/150 2026-03-10 00:34 by wdwdoo
[考研] 接受26届调剂生 +20 猪猪猪毛 2026-03-06 20/1000 2026-03-10 00:32 by Dynamic_
[考研] 294 英二数二物化 求调剂 +6 米饭团不好吃 2026-03-09 6/300 2026-03-09 23:55 by barlinike
[考研] 288求调剂 +9 王晓阳- 2026-03-09 11/550 2026-03-09 22:43 by NaiHe95273
[考研] 求调剂,数一英一274分 +4 小菲会努力 2026-03-08 4/200 2026-03-09 12:40 by 一定上岸哟_
[考研] 一志愿武理314求调剂 +4 ( ̄~ ̄;) 2026-03-08 5/250 2026-03-08 23:12 by L135790
[考研] 0701-322 求调剂 +3 jiliuxian 2026-03-06 8/400 2026-03-08 19:31 by jiliuxian
[考研] 招调剂 +20 帆船哥 2026-03-04 23/1150 2026-03-08 16:25 by taycat
[考研] 化学工程与技术316求调剂 +7 Sixuan wang 2026-03-06 7/350 2026-03-08 11:23 by 我的船我的海
[考研] 化学工程求调剂 +12 化工人999 2026-03-04 12/600 2026-03-08 11:17 by 我的船我的海
[考研] 070300化学求调剂292分 +3 打烊eee 2026-03-07 3/150 2026-03-07 20:04 by houyaoxu
[考研] 诚求调剂,323分 有专利/科研/实习经历 +6 橙子cyx 2026-03-06 8/400 2026-03-07 19:06 by 新篇章DFSS
[考研] 2026调剂】考试A区0703化学类323分 诚求接收 +3 卷柏卷柏 2026-03-05 4/200 2026-03-07 17:00 by yuandd_2001
[考研] 085600材料与化工 292分 求调剂 +6 程晴之 2026-03-06 6/300 2026-03-07 09:22 by 斩魂滴兔子!
[考研] 287求调剂 +3 看看我. 2026-03-05 6/300 2026-03-06 10:40 by Iveryant
[考研] 316求调剂 +3 林小星发大财 2026-03-05 5/250 2026-03-06 09:05 by Iveryant
[考研] 085600,一志愿郑州大学,280分求调剂 +7 Wuqi725 2026-03-05 7/350 2026-03-06 08:12 by 杨杨杨紫
[考研] 274环境工程求调剂 +6 扶柳盈江 2026-03-05 6/300 2026-03-05 13:16 by 梦天888
[考研] 347分电子信息专硕求调剂 +3 咕嘟咕嘟1 2026-03-03 3/150 2026-03-04 12:32 by 热情沙漠
[考研] 0703化学306调剂 +4 26要上岸 2026-03-03 4/200 2026-03-03 20:08 by yoyohj0427
信息提示
请填处理意见