| 3 | 1/1 | 返回列表 |
| 查看: 1768 | 回復(fù): 2 | |||
[交流]
使用MolAICal基于NAMD模擬結(jié)果計算小分子和蛋白MM/GBSA的教程 已有1人參與
|
|
使用MolAICal基于NAMD模擬結(jié)果計算小分子和蛋白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的分子動力學模擬結(jié)果,使用MolAICal計算小分子和Mpro蛋白受體MM/GBSA的方法。本教程只是一個簡單演示。為了節(jié)省運行及存儲空間,本教程僅選擇了Mpro復(fù)合物分子動力學模擬的25幀用于計算。 2.工具 2.1. 所需軟件下載地址 1) MolAICal : https://molaical.github.io 2) NAMD: https://www.ks.uiuc.edu/Research/namd/ 2.2. 操作示例文件 所有用到的操作教程文件均可在下面的網(wǎng)站下載: https://github.com/MolAICal/tutorials/tree/master/004-MMGBSA 3. 操作流程 轉(zhuǎn)到以下目錄: #> cd 004-MMGBSA 3.1. 提取蛋白與配體復(fù)合物的軌跡文件 -args: 其用法類似VMD軟件中的“atomselect”命令,比如"atomselect top protein or resname LIG",此處逗號"," 代表空格。其中腳本文件stripDCD.vmd可以在本教程材料或MolAICal軟件的“scripts”目錄里面找到。 執(zhí)行上述命令后生成complex.psf, complex.pdb 和 complex.dcd文件。將“GBIS” 和“sasa”參數(shù)設(shè)置為on。打開并按照下文內(nèi)容修改“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系統(tǒng)下運行NAMD命令,如下: 其中符號“&”代表程序在Linux系統(tǒng)中進行后臺運行,如果你使用的是Windows操作系統(tǒng),請不要用“&”,例如,命令換成這樣: 3.2. 僅提取蛋白的軌跡文件: 上述命令會生成 protein.psf, protein.pdb和protein.dcd。打開“protein.conf”,參考 “complex.conf”修改相關(guān)參數(shù)。 本教程中命令在CPU上運行。你可以選擇GPU進行運算。在Linux系統(tǒng)下運行NAMD命令,如下: 3.3. 僅提取配體的軌跡文件: 上述命令會生成 ligand.psf, ligand.pdb和ligand.dcd。打開“l(fā)igand.conf”,參考 “complex.conf”修改相關(guān)參數(shù)。 本教程中命令在CPU上運行。你可以選擇GPU進行運算。在Linux系統(tǒng)下運行NAMD命令,如下: 4. 用MolAICal計算MM/GBSA 輸出結(jié)果中給出下文所示的結(jié)合自由能△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 -------------------------------------------------------------------------- |
木蟲 (正式寫手)

| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 346分材料求調(diào)劑 +5 | snow_反季節(jié)版 2026-03-07 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 中科大材料299求調(diào)劑 +9 | DAIjiayo 2026-03-05 | 14/700 |
|
|
[考研] 化學290求調(diào)劑 +3 | 兩顆 西柚 2026-03-07 | 4/200 |
|
|
[考研] 一志愿中國石油大學(華東) 本科齊魯工業(yè)大學 求調(diào)劑 +3 | snw石 2026-03-07 | 3/150 |
|
|
[考研] 化工284求調(diào)劑 +6 | 小黑想變白 2026-03-07 | 6/300 |
|
|
[考研] 2026調(diào)劑】考試A區(qū)0703化學類323分 誠求接收 +3 | 卷柏卷柏 2026-03-05 | 4/200 |
|
|
[考研] 材料專碩調(diào)劑 +5 | 慕辰123 2026-03-05 | 8/400 |
|
|
[考研] 一志愿山東大學105500藥學專碩,總分302求調(diào)劑 +4 | 五維天空 2026-03-04 | 9/450 |
|
|
[考研] 求調(diào)劑 +4 | 呼呼?~+123456 2026-03-05 | 5/250 |
|
|
[考研] 306求調(diào)劑 +7 | Bahati 2026-03-05 | 7/350 |
|
|
[考研] 328求調(diào)劑 +6 | vuzhdkfjkx 2026-03-04 | 7/350 |
|
|
[考研] 304求調(diào)劑 +4 | 曼殊2266 2026-03-05 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | 泡了個椒 2026-03-04 | 4/200 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-04 | 3/150 |
|
|
[考研] 304分材料專碩求調(diào)劑 +11 | qiuzhigril 2026-03-03 | 14/700 |
|
|
[考研] 材料學碩080500復(fù)試調(diào)劑294 +3 | 四葉zjz 2026-03-04 | 3/150 |
|
|
[考研] 能動297求調(diào)劑,本科川大 +4 | 邵11 2026-03-04 | 4/200 |
|
|
[基金申請] 請問大家,研究風險與應(yīng)對措施那里, 大家都怎么寫呢 ? +3 | cauasen 2026-03-02 | 3/150 |
|
|
[考研] 085700資環(huán)求調(diào)劑,初始279,六級已過,英語能力強 +3 | 085700資環(huán)調(diào)劑 2026-03-03 | 4/200 |
|
|
[考研] 材料工程求調(diào)劑 +3 | 1431251 2026-03-03 | 3/150 |
|