| 5 | 1/1 | 返回列表 |
| 查看: 1769 | 回復(fù): 2 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
[交流]
使用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的分子動力學(xué)模擬結(jié)果,使用MolAICal計算小分子和Mpro蛋白受體MM/GBSA的方法。本教程只是一個簡單演示。為了節(jié)省運(yùn)行及存儲空間,本教程僅選擇了Mpro復(fù)合物分子動力學(xué)模擬的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上運(yùn)行。你可以選擇GPU進(jìn)行運(yùn)算。在Linux系統(tǒng)下運(yùn)行NAMD命令,如下: 其中符號“&”代表程序在Linux系統(tǒng)中進(jìn)行后臺運(yùn)行,如果你使用的是Windows操作系統(tǒng),請不要用“&”,例如,命令換成這樣: 3.2. 僅提取蛋白的軌跡文件: 上述命令會生成 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. 僅提取配體的軌跡文件: 上述命令會生成 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計算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 -------------------------------------------------------------------------- |
木蟲 (正式寫手)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 理學(xué)287分求調(diào)劑 +3 | 看看我. 2026-03-05 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 085701環(huán)境工程317分求調(diào)劑 +9 | 6汆尼9 2026-03-07 | 9/450 |
|
|
[考研] 085701環(huán)境工程專業(yè),初試305,均過國家A區(qū)線 +7 | 卡卡來了@ 2026-03-07 | 8/400 |
|
|
[考研] 一志愿蘇大材料工程專碩293求調(diào)劑 +4 | 鋼鐵大炮 2026-03-04 | 4/200 |
|
|
[考研] 一志愿鄭大071000分?jǐn)?shù)282求調(diào)劑 +3 | 研研顏 2026-03-05 | 7/350 |
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | G212 2026-03-03 | 6/300 |
|
|
[考研] 070300化學(xué)求調(diào)劑292分 +3 | 打烊eee 2026-03-07 | 3/150 |
|
|
[考研] 2026調(diào)劑】考試A區(qū)0703化學(xué)類323分 誠求接收 +3 | 卷柏卷柏 2026-03-05 | 4/200 |
|
|
[考研] 材料調(diào)劑 +10 | ounce. 2026-03-04 | 12/600 |
|
|
[考研] 一志愿華中科技大學(xué),化學(xué)專業(yè)344分,求調(diào)劑 +7 | 邢xing1 2026-03-02 | 7/350 |
|
|
[考研] 085600,一志愿鄭州大學(xué),280分求調(diào)劑 +7 | Wuqi725 2026-03-05 | 7/350 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-04 | 3/150 |
|
|
[考研] 材料328求調(diào)劑 +10 | 一個蘿卜02 2026-03-03 | 10/500 |
|
|
[考研] 求調(diào)劑院校 +6 | 云朵452 2026-03-02 | 12/600 |
|
|
[基金申請] 沒有青基直接申請面上,感覺自己瘋了 +6 | kevin63t 2026-03-02 | 12/600 |
|
|
[考研] 266求調(diào)劑 +7 | 哇塞王帥 2026-03-03 | 7/350 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 將未有期 2026-03-04 | 3/150 |
|
|
[基金申請] 請問大家,研究風(fēng)險與應(yīng)對措施那里, 大家都怎么寫呢 ? +3 | cauasen 2026-03-02 | 3/150 |
|
|
[考研] 26考研報考西工大材料308分求調(diào)劑 +4 | weizhong123 2026-03-01 | 5/250 |
|
|
[考研] 求調(diào)劑 +11 | yunziaaaaa 2026-03-01 | 13/650 |
|