 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 9 | 1/1 | 返回列表 |
| 查看: 3684 | 回復(fù): 8 | |||||||
御劍江湖榮譽(yù)版主 (著名寫手)
天池冶海
|
[交流]
【轉(zhuǎn)帖】GROMACS分析程序討論，在gromacs中添加Amber力場(chǎng) 已有7人參與
|
|
gromacs的分析程序很多，不知到各位平時(shí)都用到哪些呢？現(xiàn)在遇到的困惑是學(xué)會(huì)了做MD，但不會(huì)深入分析結(jié)果。我先說說我自己的,比較常規(guī)，拋磚引玉： g_energy：能量提取和分析 g_dist：研究?jī)蓚�(gè)基團(tuán)之間的距離波動(dòng) g_rms：研究某些基團(tuán)的位置偏離 g_rmsf：計(jì)算溫度因子，查看哪些殘基波動(dòng)較強(qiáng) g_rmsf可取平均結(jié)構(gòu)，用-ox參數(shù)輸出 g_gyrate：分析蛋白質(zhì)的回旋半徑變化 蛋白質(zhì)的回旋半徑反應(yīng)了蛋白質(zhì)分子的體積和形狀。同一體系的回旋半徑越大，說明體系發(fā)生了膨脹。 g_gyrate –f md.trr –s md.tpr –o gyrate.xvg g_sas：分析溶劑的可及表面積（SASA） 它是描述蛋白質(zhì)疏水性的重要手段，氨基酸殘基的疏水性是影響蛋白質(zhì)折疊的重要物理作用。 g_sas -f md.trr -s md.tpr -o area.xvg -or resarea.xvg -oa atomarea.xvg do_dssp：分析蛋白質(zhì)的二級(jí)結(jié)構(gòu)變化 group務(wù)必選擇protein，否則…… do_dssp –s md.tpr –f md.xtc –o ss.xpm g_dist和g_bond：提取兩個(gè)原子間距離隨時(shí)間變化 提取兩個(gè)原子距離隨時(shí)間變化，用g_dist和g_bond都可以。 g_dist的index寫法是，在里面加入兩個(gè)group，內(nèi)容分別是這兩個(gè)原子序號(hào)，運(yùn)行比如g_dist -f a-sol-md.xtc -s a-sol-md.tpr -n index.ndx，之后會(huì)提示選擇group，分別選擇那兩個(gè)group就行了，結(jié)果輸出在dist.xvg，第二列是距離，后面三列是x/y/z的距離。 如果用g_bond，則index里新加入一個(gè)group，把兩個(gè)原子都寫進(jìn)那一個(gè)group下。用的時(shí)候選這個(gè)group就行了，結(jié)果和g_dist的第二列的內(nèi)容一樣。 gromacs的分析程序很多，不知到各位平時(shí)都用到哪些呢？現(xiàn)在遇到的困惑是學(xué)會(huì)了做MD，但不會(huì)深入分析結(jié)果。我先說說我自己的,比較常規(guī)，拋磚引玉： g_energy：能量提取和分析g_dist：研究?jī)蓚�(gè)基團(tuán)之間的距離波動(dòng)g_rms：研究某些基團(tuán)的位置偏離g_rmsf：計(jì)算溫度因子，查看哪些殘基波動(dòng)較強(qiáng) g_rmsf可取平均結(jié)構(gòu)，用-ox參數(shù)輸出g_gyrate：分析蛋白質(zhì)的回旋半徑變化蛋白質(zhì)的回旋半徑反應(yīng)了蛋白質(zhì)分子的體積和形狀。同一體系的回旋半徑越大，說明體系發(fā)生了膨脹。 g_gyrate –f md.trr –s md.tpr –o gyrate.xvg g_sas：分析溶劑的可及表面積（SASA）它是描述蛋白質(zhì)疏水性的重要手段，氨基酸殘基的疏水性是影響蛋白質(zhì)折疊的重要物理作用。 g_sas -f md.trr -s md.tpr -o area.xvg -or resarea.xvg -oa atomarea.xvgdo_dssp：分析蛋白質(zhì)的二級(jí)結(jié)構(gòu)變化 group務(wù)必選擇protein，否則…… do_dssp –s md.tpr –f md.xtc –o ss.xpmg_dist和g_bond：提取兩個(gè)原子間距離隨時(shí)間變化 提取兩個(gè)原子距離隨時(shí)間變化，用g_dist和g_bond都可以。g_dist的index寫法是，在里面加入兩個(gè)group，內(nèi)容分別是這兩個(gè)原子序號(hào)，運(yùn)行比如g_dist -f a-sol-md.xtc -s a-sol-md.tpr -n index.ndx，之后會(huì)提示選擇group，分別選擇那兩個(gè)group就行了，結(jié)果輸出在dist.xvg，第二列是距離，后面三列是x/y/z的距離。如果用g_bond，則index里新加入一個(gè)group，把兩個(gè)原子都寫進(jìn)那一個(gè)group下。用的時(shí)候選這個(gè)group就行了，結(jié)果和g_dist的第二列的內(nèi)容一樣。 gromacs的分析程序很多，不知到各位平時(shí)都用到哪些呢？現(xiàn)在遇到的困惑是學(xué)會(huì)了做MD，但不會(huì)深入分析結(jié)果。我先說說我自己的,比較常規(guī)，拋磚引玉：g_energy：能量提取和分析g_dist：研究?jī)蓚�(gè)基團(tuán)之間的距離波動(dòng)g_rms：研究某些基團(tuán)的位置偏離g_rmsf：計(jì)算溫度因子，查看哪些殘基波動(dòng)較強(qiáng)g_rmsf可取平均結(jié)構(gòu)，用-ox參數(shù)輸出g_gyrate：分析蛋白質(zhì)的回旋半徑變化蛋白質(zhì)的回旋半徑反應(yīng)了蛋白質(zhì)分子的體積和形狀。同一體系的回旋半徑越大，說明體系發(fā)生了膨脹。g_gyrate –f md.trr –s md.tpr –o gyrate.xvgg_sas：分析溶劑的可及表面積（SASA）它是描述蛋白質(zhì)疏水性的重要手段，氨基酸殘基的疏水性是影響蛋白質(zhì)折疊的重要物理作用。g_sas -f md.trr -s md.tpr -o area.xvg -or resarea.xvg -oa atomarea.xvgdo_dssp：分析蛋白質(zhì)的二級(jí)結(jié)構(gòu)變化group務(wù)必選擇protein，否則……do_dssp –s md.tpr –f md.xtc –o ss.xpmg_dist和g_bond：提取兩個(gè)原子間距離隨時(shí)間變化提取兩個(gè)原子距離隨時(shí)間變化，用g_dist和g_bond都可以。g_dist的index寫法是，在里面加入兩個(gè)group，內(nèi)容分別是這兩個(gè)原子序號(hào)，運(yùn)行比如g_dist -f a-sol-md.xtc -s a-sol-md.tpr -n index.ndx，之后會(huì)提示選擇group，分別選擇那兩個(gè)group就行了，結(jié)果輸出在dist.xvg，第二列是距離，后面三列是x/y/z的距離。如果用g_bond，則index里新加入一個(gè)group，把兩個(gè)原子都寫進(jìn)那一個(gè)group下。用的時(shí)候選這個(gè)group就行了，結(jié)果和g_dist的第二列的內(nèi)容一樣。 在gromacs中添加Amber力場(chǎng) （GROMACS4.5自帶Amber力場(chǎng)） Install FF: 所有信息來自http://folding.stanford.edu/ffamber/ 1.下載力場(chǎng) ： http://folding.stanford.edu/ffamber/ffamber_v3.3.1.tar.gz (490 kB) http://folding.stanford.edu/ffamber/ffamber_v3.3.1-doc.tar.gz (6.54 MB)帶文獻(xiàn)例子 2.解壓縮 3.把a(bǔ)minoacides.dat、vdwradii.dat拷到top文件夾下覆蓋，實(shí)際上一股腦考過去就完了。 4.子文件夾下的東西都拷到外面，也就是top里面。 5.修改FF.dat,修改第一行的數(shù)字，加上新添的力場(chǎng)數(shù)。 6.為每個(gè)新的力場(chǎng)增加新行，如： ffamber94 AMBER94 Cornell protein/nucleic forcefield ffamber99 AMBER99 Wang protein/nucleic acid forcefield ffamber99p AMBER99p protein/nucleic forcefield ffamber03 AMBER03 Duan protein/nucleic forcefield [ Last edited by 御劍江湖 on 2011-4-9 at 09:48 ] |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 精華 | GROMACS | gromacs | 力場(chǎng) |
» 猜你喜歡
一志愿085600中科院寧波所276分求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料專業(yè)申博
已經(jīng)有3人回復(fù)
-
307求調(diào)劑
已經(jīng)有20人回復(fù)
-
285求調(diào)劑
已經(jīng)有8人回復(fù)
-
材料與化工（0856）304求B區(qū)調(diào)劑
已經(jīng)有8人回復(fù)
-
370求調(diào)劑
已經(jīng)有3人回復(fù)
-
070300化學(xué)279求調(diào)劑
已經(jīng)有5人回復(fù)
-
福建理工大學(xué)材料學(xué)院先進(jìn)合金團(tuán)隊(duì)招收考研調(diào)劑學(xué)生
已經(jīng)有4人回復(fù)
-
英一數(shù)一408，總分284，二戰(zhàn)真誠(chéng)求調(diào)劑
已經(jīng)有5人回復(fù)
-
083000環(huán)境科學(xué)與工程調(diào)劑，總分281
已經(jīng)有4人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 求Gromacs中Glycam力場(chǎng)數(shù)據(jù)包。。非常感謝
已經(jīng)有5人回復(fù)
- Gromacs中加入離子
已經(jīng)有5人回復(fù)
- 【討論】gromacs電荷不為整數(shù)
已經(jīng)有6人回復(fù)
- 【討論】gromacs 并行 mpirun
已經(jīng)有4人回復(fù)
- 【求助】如何使用gromacs中的constraint
已經(jīng)有8人回復(fù)
- 【求助】這是我安裝gromacs過程，但是出錯(cuò)阿，這個(gè)怎么辦呢？謝謝
已經(jīng)有11人回復(fù)
- 【求助】GROMACS分析二面角隨時(shí)間的變化
已經(jīng)有4人回復(fù)
- 【原創(chuàng)】在gromacs的gro中插入小分子
已經(jīng)有12人回復(fù)
- 【求助】Gromacs和NAMD比Amber快很多嗎？如何選擇？
已經(jīng)有7人回復(fù)
changle
木蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 1800.6
- 紅花: 1
- 帖子: 95
- 在線: 97.5小時(shí)
- 蟲號(hào): 1203702
- 注冊(cè): 2011-02-15
- 性別: GG
- 專業(yè): 生物大分子結(jié)構(gòu)與功能
2樓2011-04-08 17:21:39
3樓2011-05-11 12:07:11
changle
木蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 1800.6
- 紅花: 1
- 帖子: 95
- 在線: 97.5小時(shí)
- 蟲號(hào): 1203702
- 注冊(cè): 2011-02-15
- 性別: GG
- 專業(yè): 生物大分子結(jié)構(gòu)與功能
4樓2011-05-12 09:45:10
陳應(yīng)廣
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 296.1
- 散金: 1
- 帖子: 13
- 在線: 10.6小時(shí)
- 蟲號(hào): 1456893
- 注冊(cè): 2011-10-23
- 專業(yè): 蛋白質(zhì)組學(xué)
5樓2011-12-21 09:39:10
6樓2011-12-21 11:23:33
向日葵88
銀蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 319
- 帖子: 6
- 在線: 8.6小時(shí)
- 蟲號(hào): 1340083
- 注冊(cè): 2011-07-07
- 專業(yè): 有機(jī)化工
7樓2012-02-16 10:32:45
fanc232
新蟲 (小有名氣)
- 應(yīng)助: 8 (幼兒園)
- 金幣: 1238.4
- 散金: 206
- 紅花: 3
- 帖子: 253
- 在線: 108.7小時(shí)
- 蟲號(hào): 1046781
- 注冊(cè): 2010-06-24
- 性別: MM
- 專業(yè): 細(xì)胞生物學(xué)研究中的新方法
8樓2015-03-19 16:23:32
白玉浴血
木蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2513.6
- 散金: 102
- 紅花: 8
- 帖子: 838
- 在線: 243.2小時(shí)
- 蟲號(hào): 2684358
- 注冊(cè): 2013-09-27
- 性別: MM
- 專業(yè): 高分子組裝與超分子結(jié)構(gòu)
9樓2015-06-03 17:43:52
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 御劍江湖 的主題更新
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿085600中科院寧波所276分求調(diào)劑 +3 | 材料學(xué)257求調(diào)劑 2026-03-28 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 277跪求調(diào)劑 +8 | 1915668 2026-03-27 | 12/600 |
|
|
|
[考研] 279求調(diào)劑 +12 | j的立方 2026-03-29 | 12/600 |
|
|
[考研] 322求調(diào)劑 +10 | 宋明欣 2026-03-27 | 10/500 |
|
|
|
[考研] 26考研-291分-廈門大學(xué)（085601）-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +5 | min3 2026-03-24 | 6/300 |
|
|
|
[考研] 328求調(diào)劑 +8 | 嗯滴的基本都 2026-03-27 | 8/400 |
|
|
|
[考研] 求調(diào)劑323材料與化工 +10 | 1124361 2026-03-24 | 10/500 |
|
|
|
[考研] 262求調(diào)劑 +4 | ZZ..000 2026-03-30 | 4/200 |
|
|
|
[考研] 324求調(diào)劑 +9 | hanamiko 2026-03-26 | 11/550 |
|
|
|
[碩博家園] 求調(diào)劑 有機(jī)化學(xué)考研356分 +10 | Nadiums 2026-03-25 | 11/550 |
|
|
|
[考研] 070300一志愿211，312分求調(diào)劑院校 +4 | 小黃鴨寶 2026-03-30 | 4/200 |
|
|
|
[考研] 生物技術(shù)與工程 +7 | 1294608413 2026-03-25 | 8/400 |
|
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +4 | 邱gl 2026-03-26 | 7/350 |
|
|
|
[考研] 考研調(diào)劑 +7 | 小蠟新筆 2026-03-29 | 7/350 |
|
|
|
[考研] 11408軟件工程求調(diào)劑 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +7 | 爭(zhēng)取九點(diǎn)睡 2026-03-28 | 8/400 |
|
|
|
[考研] 本科新能源科學(xué)與工程，一志愿華理能動(dòng)285求調(diào)劑 +7 | AZMK 2026-03-28 | 11/550 |
|
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
|
[考研] 279 分 求調(diào)劑 +4 | 睡個(gè)好覺_16 2026-03-24 | 4/200 |
|
|
|
[考研] 334分 一志愿武理-080500 材料求調(diào)劑 +4 | 李李不服輸 2026-03-25 | 4/200 |
|