| 3 | 1/1 | 返回列表 |
| 查看: 1766 | 回復(fù): 2 | |||
[交流]
使用MolAICal基于NAMD模擬結(jié)果計(jì)算小分子和蛋白MM/GBSA的教程 已有1人參與
|
|
使用MolAICal基于NAMD模擬結(jié)果計(jì)算小分子和蛋白MM/GBSA的教程 更多教程(含英文教程)請(qǐng)見如下: MolAICal官方主頁(yè):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. 簡(jiǎn)介 在本教程中介紹了基于NAMD的分子動(dòng)力學(xué)模擬結(jié)果,使用MolAICal計(jì)算小分子和Mpro蛋白受體MM/GBSA的方法。本教程只是一個(gè)簡(jiǎn)單演示。為了節(jié)省運(yùn)行及存儲(chǔ)空間,本教程僅選擇了Mpro復(fù)合物分子動(dòng)力學(xué)模擬的25幀用于計(jì)算。 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",此處逗號(hào)"," 代表空格。其中腳本文件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上運(yùn)行。你可以選擇GPU進(jìn)行運(yùn)算。在Linux系統(tǒng)下運(yùn)行NAMD命令,如下: 其中符號(hào)“&”代表程序在Linux系統(tǒng)中進(jìn)行后臺(tái)運(yùn)行,如果你使用的是Windows操作系統(tǒng),請(qǐng)不要用“&”,例如,命令換成這樣: 3.2. 僅提取蛋白的軌跡文件: 上述命令會(huì)生成 protein.psf, protein.pdb和protein.dcd。打開“protein.conf”,參考 “complex.conf”修改相關(guān)參數(shù)。 本教程中命令在CPU上運(yùn)行。你可以選擇GPU進(jìn)行運(yùn)算。在Linux系統(tǒng)下運(yùn)行NAMD命令,如下: 3.3. 僅提取配體的軌跡文件: 上述命令會(huì)生成 ligand.psf, ligand.pdb和ligand.dcd。打開“l(fā)igand.conf”,參考 “complex.conf”修改相關(guān)參數(shù)。 本教程中命令在CPU上運(yùn)行。你可以選擇GPU進(jìn)行運(yùn)算。在Linux系統(tǒng)下運(yùn)行NAMD命令,如下: 4. 用MolAICal計(jì)算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ā)表 | |
|---|---|---|---|---|
|
[考研] 材料科學(xué)(0805)338 求調(diào)劑 +3 | xiaokang3286 2026-03-07 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 一志愿武漢理工085601,初試301分,請(qǐng)問能調(diào)劑到湖北嗎 +3 | 肖yang 2026-03-06 | 3/150 |
|
|
[考研] 081700學(xué)碩一志愿北京化工大學(xué)數(shù)二英一過六級(jí)有競(jìng)賽求調(diào)劑 +4 | galaxary 2026-03-07 | 6/300 |
|
|
[考研] 346分材料求調(diào)劑 +3 | snow_反季節(jié)版 2026-03-07 | 3/150 |
|
|
[考研] 我嘗試在小木蟲上招計(jì)算機(jī)專業(yè)的調(diào)劑學(xué)生 +4 | littlehu66 2026-03-04 | 6/300 |
|
|
[考研] 考研一志愿長(zhǎng)安大學(xué)材料與化工309分請(qǐng)求調(diào)劑 +3 | dtdxzxx 2026-03-06 | 5/250 |
|
|
[考研] 085602高分子方向求調(diào)劑 +7 | tlgudy 2026-03-04 | 7/350 |
|
|
[考研] 材料085601一志愿哈工大317 +4 | 壓迫感行 2026-03-04 | 4/200 |
|
|
[考研] 求材料調(diào)劑 +4 | berdmond 2026-03-05 | 4/200 |
|
|
[考研]
材料學(xué)學(xué)碩308分/本科東北大學(xué)/一志愿西北工業(yè)大學(xué)/
5+3
|
蘇堯幺幺 2026-03-03 | 8/400 |
|
|
[考研] 環(huán)境工程專碩307 一志愿211 四六級(jí)已過 求調(diào)劑 +5 | ccc! 2026-03-03 | 6/300 |
|
|
[考研] 0703化學(xué) 學(xué)碩 理工科均可 不區(qū)分研究方向 總分279求調(diào)劑 +7 | 1一11 2026-03-03 | 7/350 |
|
|
[考研] 環(huán)境調(diào)劑 +7 | 柒槿levana 2026-03-01 | 7/350 |
|
|
[考研] 282求調(diào)劑 +5 | 2103240126 2026-03-02 | 8/400 |
|
|
[考研] 322,求調(diào)劑 +3 | 菜菜愛玩 2026-03-04 | 3/150 |
|
|
[考研] 281求調(diào)劑 +3 | Y?l?h 2026-03-04 | 3/150 |
|
|
[考研] 331求調(diào)劑 +3 | zzZ&zZ 2026-03-03 | 3/150 |
|
|
[考研] 085600求調(diào)劑 +4 | LRZZZZZZ 2026-03-02 | 6/300 |
|
|
[考研] 材料工程269求調(diào)劑 +7 | 白刺玫 2026-03-02 | 7/350 |
|
|
[考研] 0856材料調(diào)劑 +5 | 沿岸有貝殼OUC 2026-03-02 | 5/250 |
|