24小时热门版块排行榜    

查看: 1785  |  回复: 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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 283求调剂 +6 小楼。 2026-03-12 9/450 2026-03-14 22:15 by cfxczx2015
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 307求调剂 +4 辛仁豆腐 2026-03-08 6/300 2026-03-14 04:36 by JourneyLucky
[考研] 070300化学求调剂 +6 扑风铃的猫 2026-03-08 11/550 2026-03-14 03:39 by JourneyLucky
[考研] 一志愿中国科学技术大学,086000,生物与医药,求调剂 +3 12341234789 2026-03-08 4/200 2026-03-14 03:14 by JourneyLucky
[考研] 一志愿天大化工专调剂总分338 +3 蔡大美女 2026-03-08 3/150 2026-03-14 03:12 by JourneyLucky
[考研] 一志愿安徽大学材料工程专硕313分,求调剂的学校 +8 Yu先生 2026-03-10 10/500 2026-03-14 01:04 by JourneyLucky
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 0856材料与化工309分求调剂 +6 ZyZy…… 2026-03-10 6/300 2026-03-14 00:38 by JourneyLucky
[考研] 321求调剂 +3 CUcat 2026-03-10 3/150 2026-03-14 00:25 by JourneyLucky
[考研] 279求调剂 +3 抓着星星的女孩 2026-03-10 3/150 2026-03-13 23:47 by userper
[考研] 0703化学调剂 +4 快乐的香蕉 2026-03-11 4/200 2026-03-13 22:41 by JourneyLucky
[考研] 333求调剂 +3 球球古力 2026-03-11 3/150 2026-03-13 21:27 by JourneyLucky
[考研] 0703一志愿211 285分求调剂 +4 ly3471z 2026-03-13 4/200 2026-03-13 13:00 by JourneyLucky
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-13 3/150 2026-03-13 10:43 by houyaoxu
[考研] 一志愿江南大学085701环境工程专硕总分287求调剂 +5 18266118446 2026-03-09 5/250 2026-03-11 16:51 by 2020015
[考研] 材料工程307,求调剂 +7 我要燃烧你的梦 2026-03-08 7/350 2026-03-11 09:09 by 剑诗杜康
[考研] 0856材料与化工353分求调剂 +11 NIFFFfff 2026-03-09 11/550 2026-03-10 18:36 by suyuanhai
[考研] 0703化学调剂 +3 三dd. 2026-03-10 3/150 2026-03-10 15:45 by peike
[考研] 一志愿山东大学,总分327,英语二79,有论文,有竞赛,已过四六级 +3 木木目目1 2026-03-09 3/150 2026-03-09 19:52 by yuningshan
信息提示
请填处理意见