| 8 | 1/1 | 返回列表 |
| 查看: 4211 | 回復(fù): 7 | |||
[求助]
請問一下大家有沒有分析蛋白質(zhì)二級結(jié)構(gòu)含量變化的軟件或者方法?
|
| 我用amber跑完的軌跡,想看看蛋白質(zhì)二級結(jié)構(gòu)是怎么變的。。 |
| 樓主可以考慮DSSP. 這里是鏈接http://swift.cmbi.ru.nl/gv/dssp/ |

新蟲 (初入文壇)
|
DSSP可以輸入pdb文件,輸出每個氨基酸的二級結(jié)構(gòu)。 所以先導(dǎo)出pdb,然后dssp處理,然后直接分析輸出結(jié)果就行了。 最后可以畫以下幾種圖: 1. 橫坐標(biāo)時間,縱坐標(biāo)是殘基號,用不同顏色表示不同二級結(jié)構(gòu) 2. 每個時間點各種二級結(jié)構(gòu)所占的比例,或者占初始此二級結(jié)構(gòu)的比例 |
金蟲 (小有名氣)
|
#vmd默認(rèn)使用stride計算蛋白二級結(jié)構(gòu),這里給出vmd提出dssp計算的二級結(jié)構(gòu)并輸出到一個文件中用于作圖,文件輸出格式類似Amber的輸出格式,其中X代表不規(guī)則卷曲或者loop,將以下代碼保存為test.tcl可以直接在終端下運(yùn)行以下命令: #vmd test.psf test.dcd -dispdev text -e test.tcl #如果是amber的拓?fù)湮募蛙壽E則為 #vmd -parm7 test.prmtop -crdbox test.mdcrd -dispdev text -e test.tcl #輸出的文件為sec4traj.txt 1 # dealing with data in DSSP file 2 proc myproc1 {a b framei} { 3 # creat a empty list 4 set newsecdata 5 set linenumber 0 6 foreach a $b { 7 incr linenumber 8 if {[string equal [string index $a 2] #]} { 9 puts "The line number is $linenumber" 10 set secdata [lrange $b $linenumber end] 11 # puts $secdata 12 foreach secinfo $secdata { 13 set sectype [string index $secinfo 16] 14 if {[string is space $sectype]} then { 15 lappend newsecdata X 16 } else { 17 lappend newsecdata $sectype 18 } 19 } 20 } 21 } 22 set mynewsecdata [join $newsecdata \t] 23 set secdataframe [linsert $mynewsecdata 0 $framei] 24 puts $secdataframe 25 set sectraj [open sec4traj.txt a] 26 puts $sectraj $secdataframe 27 close $sectraj 28 } 29 30 # generate each frame pdb file and analyze secondary structure with dssp 31 set n [molinfo top get numframes] 32 puts "===Analysis number of frames: $n" 33 set myprotein [atomselect top protein] 34 #write resids 35 set myresids [lsort -integer -unique [$myprotein get resid]] 36 37 set sectraj [open sec4traj.txt w] 38 puts $sectraj $myresids 39 close $sectraj 40 41 #analysis of frames 42 for { set i 0 } { $i < $n } { incr i } { 43 $myprotein frame $i 44 $myprotein update 45 $myprotein writepdb myptrotein_$i.pdb 46 exec mkdssp -i myptrotein_$i.pdb -o myptrotein_$i.dssp 47 48 set dsspfile [open myptrotein_$i.dssp r] 49 set dsspline [read $dsspfile ] 50 set dsspdata [split $dsspline \n] 51 close $dsspfile 52 myproc1 splitdata $dsspdata $i 53 54 55 file delete myptrotein_$i.pdb 56 file delete myptrotein_$i.dssp 57 } |
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 328求調(diào)劑,英語六級551,有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-16 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 一志愿天大材料與化工(085600)總分338 +5 | 蔡大美女 2026-03-13 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑 +8 | 妮妮ninicgb 2026-03-15 | 12/600 |
|
|
[考研] 085410人工智能專碩317求調(diào)劑(0854都可以) +3 | xbxudjdn 2026-03-18 | 3/150 |
|
|
[考研] 266求調(diào)劑 +5 | 陽陽哇塞 2026-03-14 | 9/450 |
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
[考研] 0854可跨調(diào)劑,一作一項核心論文五項專利,省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
[考研] 288求調(diào)劑,一志愿華南理工大學(xué)071005 +4 | ioodiiij 2026-03-17 | 4/200 |
|
|
[考研] 環(huán)境工程調(diào)劑 +8 | 大可digkids 2026-03-16 | 8/400 |
|
|
[考研] 材料與化工求調(diào)劑 +6 | 為學(xué)666 2026-03-16 | 6/300 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
[考研] 一志愿華中師范071000,325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy! 2026-03-13 | 6/300 |
|
|
[考研] 學(xué)碩285求調(diào)劑 +13 | Wisjxn 2026-03-12 | 46/2300 |
|
|
[考研] 招收0805(材料)調(diào)劑 +3 | 18595523086 2026-03-13 | 3/150 |
|
|
[考研] [0860]321分求調(diào)劑,ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|
|
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9 | 2026調(diào)劑俠 2026-03-12 | 9/450 |
|