24小时热门版块排行榜    

查看: 1771  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上项目还需要AI说明吗? +3 liyundong 2026-03-08 3/150 2026-03-09 22:30 by kingkocxr
[考研] A区一本交叉课题组,低分调剂,招收机械电子信息通信等交叉方向 +54 lisimayy 2026-03-04 67/3350 2026-03-09 19:44 by 李玲玉123
[考研] 0703化学求调剂 +4 很老实人 2026-03-09 4/200 2026-03-09 13:38 by macy2011
[考研] 337求调剂 +3 睡醒,。 2026-03-09 3/150 2026-03-09 10:02 by 求调剂zz
[考研] 求0856调剂 +5 小力气珂珂 2026-03-08 5/250 2026-03-09 06:42 by houyaoxu
[考研] 290 材料与化工求调剂 +7 Nebulala 2026-03-08 7/350 2026-03-09 05:21 by houyaoxu
[考研] 297求调剂 +3 胡达灵 2026-03-05 5/250 2026-03-08 23:13 by 清风月
[考研] 一志愿武理314求调剂 +4 ( ̄~ ̄;) 2026-03-08 5/250 2026-03-08 23:12 by L135790
[考研] 考研求调剂,工科,总分304 +4 矮子不想长高了 2026-03-08 7/350 2026-03-08 22:21 by 小木虫tim
[考研] 301求调剂 +11 朝天椒的双马尾 2026-03-05 11/550 2026-03-08 21:50 by 知足常乐的乐
[考研] 085701环境工程专业,初试305,均过国家A区线 +7 卡卡来了@ 2026-03-07 8/400 2026-03-08 06:40 by 刘兵
[考研] 308求调剂 +7 倘若起风了呢 2026-03-05 9/450 2026-03-06 22:59 by 倘若起风了呢
[考研] 306求调剂 +7 Bahati 2026-03-05 7/350 2026-03-06 22:11 by 星空星月
[考研] 275求调剂 +4 大爆炸难民 2026-03-06 5/250 2026-03-06 09:21 by guoweigw
[考研] 材料调剂 +5 ws 上岸鸭 2026-03-05 5/250 2026-03-05 20:06 by houyaoxu
[考研] 070300化学 280 一志愿太原理工 求调剂 +8 拾玖壹 2026-03-04 8/400 2026-03-05 08:52 by anny19840123
[考研] 一志愿985材料与化工 326分求调剂 +3 Hz795795 2026-03-04 3/150 2026-03-04 20:54 by wutongshun
[考研] 学硕材料275调剂 +9 路三三 2026-03-03 9/450 2026-03-04 17:02 by 梦天888
[考研] 0703化学 一志愿东北大学 279分求调剂 +4 聂乐平 2026-03-03 4/200 2026-03-04 16:31 by zhukairuo
[考研] 材料专硕346求调剂 +3 旺一下 2026-03-04 3/150 2026-03-04 16:26 by sslc1985
信息提示
请填处理意见