 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3686 | 回復(fù): 8 | |||||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||||||
御劍江湖榮譽(yù)版主 (著名寫手)
天池冶海
|
[交流]
【轉(zhuǎn)帖】GROMACS分析程序討論，在gromacs中添加Amber力場 已有7人參與
|
||||||
|
gromacs的分析程序很多，不知到各位平時都用到哪些呢？現(xiàn)在遇到的困惑是學(xué)會了做MD，但不會深入分析結(jié)果。我先說說我自己的,比較常規(guī)，拋磚引玉： g_energy：能量提取和分析 g_dist：研究兩個基團(tuán)之間的距離波動 g_rms：研究某些基團(tuán)的位置偏離 g_rmsf：計算溫度因子，查看哪些殘基波動較強(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ì)的二級結(jié)構(gòu)變化 group務(wù)必選擇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)容分別是這兩個原子序號，運(yùn)行比如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，把兩個原子都寫進(jìn)那一個group下。用的時候選這個group就行了，結(jié)果和g_dist的第二列的內(nèi)容一樣。 gromacs的分析程序很多，不知到各位平時都用到哪些呢？現(xiàn)在遇到的困惑是學(xué)會了做MD，但不會深入分析結(jié)果。我先說說我自己的,比較常規(guī)，拋磚引玉： g_energy：能量提取和分析g_dist：研究兩個基團(tuán)之間的距離波動g_rms：研究某些基團(tuán)的位置偏離g_rmsf：計算溫度因子，查看哪些殘基波動較強(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ì)的二級結(jié)構(gòu)變化 group務(wù)必選擇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)容分別是這兩個原子序號，運(yùn)行比如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，把兩個原子都寫進(jìn)那一個group下。用的時候選這個group就行了，結(jié)果和g_dist的第二列的內(nèi)容一樣。 gromacs的分析程序很多，不知到各位平時都用到哪些呢？現(xiàn)在遇到的困惑是學(xué)會了做MD，但不會深入分析結(jié)果。我先說說我自己的,比較常規(guī)，拋磚引玉：g_energy：能量提取和分析g_dist：研究兩個基團(tuán)之間的距離波動g_rms：研究某些基團(tuán)的位置偏離g_rmsf：計算溫度因子，查看哪些殘基波動較強(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ì)的二級結(jié)構(gòu)變化group務(wù)必選擇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)容分別是這兩個原子序號，運(yùn)行比如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，把兩個原子都寫進(jìn)那一個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)帶文獻(xiàn)例子 2.解壓縮 3.把a(bǔ)minoacides.dat、vdwradii.dat拷到top文件夾下覆蓋，實(shí)際上一股腦考過去就完了。 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 ] |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 精華 | GROMACS | gromacs | 力場 |
» 猜你喜歡
一志愿大連理工大學(xué)材料求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料與化工調(diào)劑一志愿大連海事085600，349
已經(jīng)有6人回復(fù)
-
一志愿中海洋320化學(xué)工程與技術(shù)學(xué)碩求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料科學(xué)與工程調(diào)劑
已經(jīng)有11人回復(fù)
-
311求調(diào)劑一志愿合肥工業(yè)大學(xué)
已經(jīng)有4人回復(fù)
-
一志愿食品科學(xué)與工程083200求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料科學(xué)與工程求調(diào)劑
已經(jīng)有10人回復(fù)
-
一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑
已經(jīng)有11人回復(fù)
-
一志愿 南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩
已經(jīng)有4人回復(fù)
-
085600材料與化工調(diào)劑
已經(jīng)有15人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 求Gromacs中Glycam力場數(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過程，但是出錯阿，這個怎么辦呢？謝謝
已經(jīng)有11人回復(fù)
- 【求助】GROMACS分析二面角隨時間的變化
已經(jīng)有4人回復(fù)
- 【原創(chuàng)】在gromacs的gro中插入小分子
已經(jīng)有12人回復(fù)
- 【求助】Gromacs和NAMD比Amber快很多嗎？如何選擇？
已經(jīng)有7人回復(fù)
白玉浴血
木蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2513.6
- 散金: 102
- 紅花: 8
- 帖子: 838
- 在線: 243.2小時
- 蟲號: 2684358
- 注冊: 2013-09-27
- 性別: MM
- 專業(yè): 高分子組裝與超分子結(jié)構(gòu)
9樓2015-06-03 17:43:52
changle
木蟲 (小有名氣)
- 應(yīng)助: 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
木蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 1800.6
- 紅花: 1
- 帖子: 95
- 在線: 97.5小時
- 蟲號: 1203702
- 注冊: 2011-02-15
- 性別: GG
- 專業(yè): 生物大分子結(jié)構(gòu)與功能
4樓2011-05-12 09:45:10
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿 南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +4 | @taotao 2026-03-30 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 285求調(diào)劑 +5 | AZMK 2026-03-30 | 8/400 |
|
|
|
[考研] 材料求調(diào)劑 一志愿哈工大總分298分，前三科223分 +10 | dongfang59 2026-03-27 | 10/500 |
|
|
|
[考研] 317求調(diào)劑 +8 | 十閑wx 2026-03-24 | 8/400 |
|
|
|
[考研] 生物學(xué)學(xué)碩，一志愿湖南大學(xué)，初試成績338 +7 | YYYYYNNNNN 2026-03-26 | 9/450 |
|
|
|
[考研] 一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑 +25 | 稻妻小編 2026-03-26 | 25/1250 |
|
|
|
[考研] 各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩 +10 | Koxui 2026-03-28 | 10/500 |
|
|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +12 | 辰沐5211314 2026-03-26 | 12/600 |
|
|
|
[考研] 085600，專業(yè)課化工原理，320分求調(diào)劑 +5 | 大饞小子 2026-03-29 | 5/250 |
|
|
|
[考研] 085602 化學(xué)工程專碩 340分求調(diào)劑 +4 | qianbai11 2026-03-29 | 4/200 |
|
|
|
[考研] 材料與化工272求調(diào)劑 +21 | 阿斯蒂芬2004 2026-03-28 | 21/1050 |
|
|
|
[考研] 一志愿中南大學(xué)化學(xué)0703總分337求調(diào)劑 +6 | niko- 2026-03-27 | 6/300 |
|
|
|
[考研] 調(diào)劑310 +12 | 溫柔的晚安 2026-03-25 | 13/650 |
|
|
|
[考研] 322求調(diào)劑 +5 | 舊吢 2026-03-24 | 5/250 |
|
|
|
[考研] 286求調(diào)劑 +12 | PolarBear11 2026-03-26 | 12/600 |
|
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
|
[考研] 292求調(diào)劑 +4 | 求求了收下我吧?/a> 2026-03-26 | 4/200 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
|
[考研] 07化學(xué)303求調(diào)劑 +5 | 睿08 2026-03-25 | 5/250 |
|
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-25 | 3/150 |
|