| 查看: 4071 | 回复: 7 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
请问一下大家有没有分析蛋白质二级结构含量变化的软件或者方法？
|
||
| 我用amber跑完的轨迹，想看看蛋白质二级结构是怎么变的。。 |
回复此楼» 猜你喜欢
论文终于录用啦！满足毕业条件了
已经有13人回复
-
2025年遐想
已经有5人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
-
求个博导看看
已经有18人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 怎样进行蛋白质高级结构的比对分析
已经有7人回复
- 蛋白质结构域的分析
已经有4人回复
- 蛋白质立体结构3D图 软件
已经有4人回复
- 对蛋白质的结构进行模拟和性质分析一般采用仕么软件
已经有4人回复
- 蛋白质二级结构红外测定
已经有15人回复
- 蛋白质二级结构、三级结构的预测
已经有8人回复
- 关于预测蛋白质分子量、二级结构
已经有6人回复
- 蛋白质二级结构预测
已经有9人回复
- 关于红外分析蛋白质二级结构的文章--投稿求助
已经有5人回复
- 急求：得到蛋白质后，想对其进行二级结构的检测，求方法和机构
已经有9人回复
- 蛋白质二级结构分析
已经有12人回复
- 用protean分析蛋白质二级结构
已经有4人回复
- 【求助/交流】请教各位：如何分析蛋白质的结构？
已经有4人回复
- 求一个认可度较高蛋白质二级结构预测软件
已经有1人回复
- 【求助】清楚的蛋白质二级结构图
已经有6人回复
- 【求助】现在还有没有可用的可画蛋白质拓扑结构的软件或网站？
已经有4人回复
xuxm03
金虫 (小有名气)
- 应助: 24 (小学生)
- 金币: 1053.6
- 散金: 5
- 红花: 2
- 帖子: 276
- 在线: 461.2小时
- 虫号: 757207
- 注册: 2009-04-25
- 专业: 生物化学
|
#vmd默认使用stride计算蛋白二级结构，这里给出vmd提出dssp计算的二级结构并输出到一个文件中用于作图，文件输出格式类似Amber的输出格式，其中X代表不规则卷曲或者loop，将以下代码保存为test.tcl可以直接在终端下运行以下命令： #vmd test.psf test.dcd -dispdev text -e test.tcl #如果是amber的拓扑文件和轨迹则为 #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 } |
7楼2013-08-01 18:10:33