| 查看: 4535 | 回復(fù): 12 | ||||||
xueer小雯金蟲 (初入文壇)
|
[求助]
請(qǐng)教gromacs蛋白與小分子相互作用分析 已有3人參與
|
|
是這樣的,我先做了蛋白與小分子的平衡,而后將小分子從蛋白的活性口袋中拉伸出來。 現(xiàn)在主要想分析整個(gè)過程中拉力的變化,看看峰值什么的,這個(gè)不知道該用什么命令來執(zhí)行。 然后再分析分析氫鍵作用,這個(gè)我知道用g_bond。 然后還有疏水作用和水橋,不知道這兩個(gè)用哪個(gè)命令可以執(zhí)行。 請(qǐng)教大家,請(qǐng)不吝賜教! 多謝! |
gromacs | xuexijisuan | 藥物化學(xué) | gromacs |
|
是這樣的,我先做了蛋白與小分子的平衡,而后將小分子從蛋白的活性口袋中拉伸出來。 現(xiàn)在主要想分析整個(gè)過程中拉力的變化,看看峰值什么的, 這個(gè)可以算 PMF 然后再分析分析氫鍵作用,這個(gè)我知道用g_bond ,貌似不是g_bond, 是g_hbond |
|
本帖內(nèi)容被屏蔽 |
木蟲 (小有名氣)
| 樓上說的都是很常用的,如果你是新手只想簡(jiǎn)單的看力和相互作用能量變化的話建議使用g_traj提取力文件,g_energy提取相互作用能量,g_hbond分析氫鍵,g_saltbr分析鹽橋;PMF可以用g_wham,疏水作用和水橋好像是不能直接分析的,需要借助其他軟件,水橋你要是用g_hbond自己計(jì)算就麻煩了。當(dāng)然你還可以輔助一些比如距離g_dist等分析峰值,有些文章是用力、距離和能量共同確定峰值的。 |
金蟲 (初入文壇)
金蟲 (初入文壇)
木蟲 (小有名氣)
|
不好意思,更正一下,g_traj不適合用作拉伸力的提前,他提出來的多是和坐標(biāo)相關(guān)的。拉伸用到的力是隨時(shí)間變化的那種,你打開你的文件夾,查找有沒有一個(gè)pullf的文件,在你設(shè)置pulling的時(shí)候,這個(gè)文件是默認(rèn)輸出的,如果沒有,那么你用mdrun -rerun -pf提前下,rerun的速度回很快 |
金蟲 (初入文壇)
金蟲 (初入文壇)
木蟲 (小有名氣)
|
在mdp文件設(shè)置的時(shí)候一般都是用的umbrella。恕水平有限不能將傘形抽樣解釋的很清楚,拙見如下,歡迎高手指正。 傘形取樣是為了克服你模擬時(shí)間太短而導(dǎo)致的取樣不足而生的,它能服務(wù)于計(jì)算PMF。模擬基本都在ns級(jí),現(xiàn)在好像到了ms級(jí),但是相對(duì)于真實(shí)反應(yīng)還是太短,這么短時(shí)間,模擬無法完全取樣所有過程(特別是高能量狀態(tài),能壘障礙太多啊)。傘形取樣就是把反應(yīng)坐標(biāo)分成若干個(gè)window,然后分別計(jì)算window的概率分別,最后wham出PMF。PMF這東西,拉伸文章里到是很常用,但是你不用也無妨,除非審稿人非要你補(bǔ)這數(shù)據(jù).... |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 293求調(diào)劑 +3 | 濤濤Wjt 2026-03-22 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 289材料與化工(085600)B區(qū)求調(diào)劑 +3 | 這么名字咋樣 2026-03-22 | 4/200 |
|
|
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +3 | niko- 2026-03-22 | 3/150 |
|
|
[考研] 求調(diào)劑 +7 | Auroracx 2026-03-22 | 7/350 |
|
|
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3 | 夏夏夏小正 2026-03-17 | 5/250 |
|
|
[考研] 0805 316求調(diào)劑 +3 | 大雪深藏 2026-03-18 | 3/150 |
|
|
[考研] 299求調(diào)劑 +4 | 某某某某位 2026-03-21 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | 白QF 2026-03-21 | 3/150 |
|
|
[考研] 346求調(diào)劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +12 | 墨墨漠 2026-03-18 | 13/650 |
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊(duì) 2026-03-20 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|