24小时热门版块排行榜    

查看: 1772  |  回复: 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的回帖

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的回帖
查看全部 3 个回答

zxhuizhang

木虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
我想问一下,关于linux系统的电脑配置问题。大概需要什么样子的
见此
2楼2021-02-02 12:12:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 295复试调剂 +4 简木ChuFront 2026-03-09 4/200 2026-03-09 21:49 by wangjihu
[考研] 085600材料与化工 326 求调剂 +3 热爱生活ing 2026-03-09 3/150 2026-03-09 20:12 by 雾散后相遇lc
[考研] 086000生物与医药319分求调剂 +4 Tolkien 2026-03-07 4/200 2026-03-09 13:49 by macy2011
[考研] 0703化学求调剂 +4 很老实人 2026-03-09 4/200 2026-03-09 13:38 by macy2011
[考博] 申博 +5 添菜了哈 2026-03-04 8/400 2026-03-09 10:52 by jamesqi
[考研] 中科大材料299求调剂 +10 DAIjiayo 2026-03-05 16/800 2026-03-09 10:45 by 斩魂滴兔子!
[考研] 0856材料与化工290求调剂 +7 Nebulala 2026-03-08 8/400 2026-03-09 05:25 by houyaoxu
[考研] 一志愿武理314求调剂 +4 ( ̄~ ̄;) 2026-03-08 5/250 2026-03-08 23:12 by L135790
[考研] 一志愿郑大071000分数282求调剂 +3 研研颜 2026-03-05 7/350 2026-03-07 22:33 by 帆船哥
[考研] 求调剂,不管什么专业,我是可塑造的人才一枚,希望遇到知己老师捞捞我 +4 13102137290 2026-03-06 5/250 2026-03-07 21:21 by lissomchan
[考研] 材料专硕323求调剂 +5 李白26 2026-03-07 5/250 2026-03-07 16:46 by njzyff
[考研] 求调剂 +4 呼呼?~+123456 2026-03-06 4/200 2026-03-06 23:11 by L135790
[考研] 求调剂 +5 danyyyy 2026-03-04 5/250 2026-03-06 08:16 by Leeding1356
[考研] 考研282分求调剂,接受跨专业 +4 刘淄博 2026-03-04 7/350 2026-03-05 22:56 by 刘淄博
[考研] 一志愿武理085601专硕347分求调剂 +4 啊欧欧欧 2026-03-04 5/250 2026-03-05 19:27 by Leding1356
[考研] 武理材料工程302调剂 +3 Doleres 2026-03-03 6/300 2026-03-04 15:54 by zhukairuo
[考研] 281求调剂 +3 我是小小葱葱 2026-03-03 5/250 2026-03-04 14:23 by kakakapanpan
[考研] 0703化学学硕 288分 +4 旅行中的紫葡萄 2026-03-03 4/200 2026-03-03 23:41 by ch1979
[考研] 306求调剂 +7 张张张张oo 2026-03-03 7/350 2026-03-03 21:35 by L135790
[考研] 0856材料工程,初试313调剂 +7 卖个关子吧 2026-03-03 7/350 2026-03-03 21:33 by L135790
信息提示
请填处理意见