| 8 | 1/1 | 返回列表 |
| 查看: 4207 | 回復(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ā)表 | |
|---|---|---|---|---|
|
[考研] 085600材料與化工調(diào)劑 324分 +8 | llllkkkhh 2026-03-18 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)(081702)315分求調(diào)劑 +10 | yangfz 2026-03-17 | 10/500 |
|
|
[考研] 材料專碩英一數(shù)二306 +4 | z1z2z3879 2026-03-18 | 4/200 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | pupcoco 2026-03-17 | 7/350 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
[考研] 268求調(diào)劑 +6 | 簡單點0 2026-03-17 | 6/300 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[考研] 290求調(diào)劑 +6 | 孔志浩 2026-03-12 | 11/550 |
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨; 2026-03-16 | 3/150 |
|
|
[考研] 機(jī)械專碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
[基金申請]
今年的國基金是打分制嗎?
50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請 +6 | zlingli 2026-03-13 | 6/300 |
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy! 2026-03-13 | 6/300 |
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|
|
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8 | 努力奮斗112 2026-03-12 | 9/450 |
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|