 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3685 | 回復: 8 | |||||||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||||||
御劍江湖榮譽版主 (著名寫手)
天池冶海
|
[交流]
【轉(zhuǎn)帖】GROMACS分析程序討論，在gromacs中添加Amber力場 已有7人參與
|
||||||
|
gromacs的分析程序很多，不知到各位平時都用到哪些呢？現(xiàn)在遇到的困惑是學會了做MD，但不會深入分析結(jié)果。我先說說我自己的,比較常規(guī)，拋磚引玉： g_energy：能量提取和分析 g_dist：研究兩個基團之間的距離波動 g_rms：研究某些基團的位置偏離 g_rmsf：計算溫度因子，查看哪些殘基波動較強 g_rmsf可取平均結(jié)構(gòu)，用-ox參數(shù)輸出 g_gyrate：分析蛋白質(zhì)的回旋半徑變化 蛋白質(zhì)的回旋半徑反應了蛋白質(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ì)的二級結(jié)構(gòu)變化 group務必選擇protein，否則…… do_dssp –s md.tpr –f md.xtc –o ss.xpm g_dist和g_bond：提取兩個原子間距離隨時間變化 提取兩個原子距離隨時間變化，用g_dist和g_bond都可以。 g_dist的index寫法是，在里面加入兩個group，內(nèi)容分別是這兩個原子序號，運行比如g_dist -f a-sol-md.xtc -s a-sol-md.tpr -n index.ndx，之后會提示選擇group，分別選擇那兩個group就行了，結(jié)果輸出在dist.xvg，第二列是距離，后面三列是x/y/z的距離。 如果用g_bond，則index里新加入一個group，把兩個原子都寫進那一個group下。用的時候選這個group就行了，結(jié)果和g_dist的第二列的內(nèi)容一樣。 gromacs的分析程序很多，不知到各位平時都用到哪些呢？現(xiàn)在遇到的困惑是學會了做MD，但不會深入分析結(jié)果。我先說說我自己的,比較常規(guī)，拋磚引玉： g_energy：能量提取和分析g_dist：研究兩個基團之間的距離波動g_rms：研究某些基團的位置偏離g_rmsf：計算溫度因子，查看哪些殘基波動較強 g_rmsf可取平均結(jié)構(gòu)，用-ox參數(shù)輸出g_gyrate：分析蛋白質(zhì)的回旋半徑變化蛋白質(zhì)的回旋半徑反應了蛋白質(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ì)的二級結(jié)構(gòu)變化 group務必選擇protein，否則…… do_dssp –s md.tpr –f md.xtc –o ss.xpmg_dist和g_bond：提取兩個原子間距離隨時間變化 提取兩個原子距離隨時間變化，用g_dist和g_bond都可以。g_dist的index寫法是，在里面加入兩個group，內(nèi)容分別是這兩個原子序號，運行比如g_dist -f a-sol-md.xtc -s a-sol-md.tpr -n index.ndx，之后會提示選擇group，分別選擇那兩個group就行了，結(jié)果輸出在dist.xvg，第二列是距離，后面三列是x/y/z的距離。如果用g_bond，則index里新加入一個group，把兩個原子都寫進那一個group下。用的時候選這個group就行了，結(jié)果和g_dist的第二列的內(nèi)容一樣。 gromacs的分析程序很多，不知到各位平時都用到哪些呢？現(xiàn)在遇到的困惑是學會了做MD，但不會深入分析結(jié)果。我先說說我自己的,比較常規(guī)，拋磚引玉：g_energy：能量提取和分析g_dist：研究兩個基團之間的距離波動g_rms：研究某些基團的位置偏離g_rmsf：計算溫度因子，查看哪些殘基波動較強g_rmsf可取平均結(jié)構(gòu)，用-ox參數(shù)輸出g_gyrate：分析蛋白質(zhì)的回旋半徑變化蛋白質(zhì)的回旋半徑反應了蛋白質(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ì)的二級結(jié)構(gòu)變化group務必選擇protein，否則……do_dssp –s md.tpr –f md.xtc –o ss.xpmg_dist和g_bond：提取兩個原子間距離隨時間變化提取兩個原子距離隨時間變化，用g_dist和g_bond都可以。g_dist的index寫法是，在里面加入兩個group，內(nèi)容分別是這兩個原子序號，運行比如g_dist -f a-sol-md.xtc -s a-sol-md.tpr -n index.ndx，之后會提示選擇group，分別選擇那兩個group就行了，結(jié)果輸出在dist.xvg，第二列是距離，后面三列是x/y/z的距離。如果用g_bond，則index里新加入一個group，把兩個原子都寫進那一個group下。用的時候選這個group就行了，結(jié)果和g_dist的第二列的內(nèi)容一樣。 在gromacs中添加Amber力場 （GROMACS4.5自帶Amber力場） Install FF: 所有信息來自http://folding.stanford.edu/ffamber/ 1.下載力場 ： 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)帶文獻例子 2.解壓縮 3.把aminoacides.dat、vdwradii.dat拷到top文件夾下覆蓋，實際上一股腦考過去就完了。 4.子文件夾下的東西都拷到外面，也就是top里面。 5.修改FF.dat,修改第一行的數(shù)字，加上新添的力場數(shù)。 6.為每個新的力場增加新行，如： 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 ] |
回復此樓» 收錄本帖的淘帖專輯推薦
 精華 | GROMACS | gromacs | 力場 |
» 猜你喜歡
學碩274求調(diào)劑
已經(jīng)有11人回復
-
一志愿鄭大材料工程290求調(diào)劑
已經(jīng)有12人回復
-
環(huán)境工程 085701，267求調(diào)劑
已經(jīng)有8人回復
-
083000學碩274求調(diào)劑
已經(jīng)有10人回復
-
一志愿085600中科院寧波所276分求調(diào)劑
已經(jīng)有4人回復
-
材料專業(yè)申博
已經(jīng)有3人回復
-
307求調(diào)劑
已經(jīng)有20人回復
-
285求調(diào)劑
已經(jīng)有8人回復
-
材料與化工（0856）304求B區(qū)調(diào)劑
已經(jīng)有8人回復
-
370求調(diào)劑
已經(jīng)有3人回復
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 求Gromacs中Glycam力場數(shù)據(jù)包。。非常感謝
已經(jīng)有5人回復
- Gromacs中加入離子
已經(jīng)有5人回復
- 【討論】gromacs電荷不為整數(shù)
已經(jīng)有6人回復
- 【討論】gromacs 并行 mpirun
已經(jīng)有4人回復
- 【求助】如何使用gromacs中的constraint
已經(jīng)有8人回復
- 【求助】這是我安裝gromacs過程，但是出錯阿，這個怎么辦呢？謝謝
已經(jīng)有11人回復
- 【求助】GROMACS分析二面角隨時間的變化
已經(jīng)有4人回復
- 【原創(chuàng)】在gromacs的gro中插入小分子
已經(jīng)有12人回復
- 【求助】Gromacs和NAMD比Amber快很多嗎？如何選擇？
已經(jīng)有7人回復
7樓2012-02-16 10:32:45
changle
木蟲 (小有名氣)
- 應助: 1 (幼兒園)
- 金幣: 1800.6
- 紅花: 1
- 帖子: 95
- 在線: 97.5小時
- 蟲號: 1203702
- 注冊: 2011-02-15
- 性別: GG
- 專業(yè): 生物大分子結(jié)構(gòu)與功能
2樓2011-04-08 17:21:39
3樓2011-05-11 12:07:11
changle
木蟲 (小有名氣)
- 應助: 1 (幼兒園)
- 金幣: 1800.6
- 紅花: 1
- 帖子: 95
- 在線: 97.5小時
- 蟲號: 1203702
- 注冊: 2011-02-15
- 性別: GG
- 專業(yè): 生物大分子結(jié)構(gòu)與功能
4樓2011-05-12 09:45:10
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 學碩274求調(diào)劑 +11 | Li李魚 2026-03-26 | 11/550 |
|
|---|---|---|---|---|
|
|
[考研] 福建理工大學材料學院先進合金團隊招收考研調(diào)劑學生 +3 | 大華金商都 2026-03-30 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | 圖鑒212 2026-03-30 | 3/150 |
|
|
|
[考研] 哈爾濱工業(yè)大學材料與化工專碩378求調(diào)劑 +3 | 塔比烏斯 2026-03-30 | 3/150 |
|
|
|
[考研] 26考研-291分-廈門大學（085601）-柔性電子學院材料工程專業(yè)求調(diào)劑 +5 | min3 2026-03-24 | 6/300 |
|
|
[考研] 348求調(diào)劑 +6 | 小懶蟲不懶了 2026-03-28 | 6/300 |
|
|
|
[考研] 375求調(diào)劑 +6 | 雨夏整夜 2026-03-29 | 6/300 |
|
|
|
[考研] 330分求調(diào)劑 +5 | qzenlc 2026-03-29 | 5/250 |
|
|
|
[考研] 數(shù)一英一271專碩（085401）求調(diào)劑，可跨 +7 | 前行必有光 2026-03-28 | 8/400 |
|
|
|
[考研] 一志愿太原理工安全工程300分，求調(diào)劑 +5 | 0857求調(diào)劑. 2026-03-24 | 6/300 |
|
|
|
[考研] 316求調(diào)劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
|
[考研] 266分，求材料冶金能源化工等調(diào)劑 +7 | 哇呼哼呼哼 2026-03-27 | 9/450 |
|
|
|
[考研] 藥學105500求調(diào)劑 +3 | Ssun。。 2026-03-28 | 3/150 |
|
|
|
[考研] 085405 考的11408求各位老師帶走 +3 | Qiu學ing 2026-03-28 | 3/150 |
|
|
|
[考研] 安徽大學專碩生物與醫(yī)藥專業(yè)(086000)324分，英語已過四六級，六級521，求調(diào)劑 +4 | 美味可樂雞翅 2026-03-26 | 4/200 |
|
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
|
[考研] 351求調(diào)劑 +4 | 麥克阿磊 2026-03-24 | 4/200 |
|
|
|
[考研] 一志愿天津大學339材料與化工求調(diào)劑 +3 | 江往賣魚 2026-03-26 | 3/150 |
|
|
|
[考研] 086003食品工程求調(diào)劑 +6 | 淼淼111 2026-03-24 | 6/300 |
|
|
|
[考研] 277分求調(diào)劑，跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|