| 查看: 2025 | 回复: 6 | ||||||
[交流]
【讨论】大家看看我的RMSD图，给些指点，谢谢了 已有3人参与
|
|
跑完平衡动力学，在求RMSD for individual residues时发现输出的RMSD图有些诡异，我用的是这个程序： # % $Id: residue_rmsd.tcl,v 1.4 2006/03/06 23:56:46 timisgro Exp $ proc rmsd_residue_over_time {{mol top} res} { # use frame 0 for the reference set reference [atomselect $mol "protein" frame 0] # the frame being compared set compare [atomselect $mol "protein"] #make a selection with all atoms set all [atomselect top all] #get the number of frames set num_steps [molinfo $mol get numframes] #open file for writing set fil [open residue_rmsd.dat w] foreach r $res { set rmsd($r) 0 } #loop over all frames in the trajectory for {set frame 0} {$frame < $num_steps} {incr frame} { puts "Calculating rmsd for frame $frame ..." # get the correct frame $compare frame $frame $all frame $frame # compute the transformation set trans_mat [measure fit $compare $reference] # do the alignment $all move $trans_mat # compute the RMSD #loop through all residues foreach r $res { set ref [atomselect $mol "protein and resid $r and noh" frame 0] set comp [atomselect $mol "protein and resid $r and noh" frame $frame] set rmsd($r) [expr $rmsd($r) + [measure rmsd $comp $ref]] $comp delete $ref delete } } set ave 0 foreach r $res { set rmsd($r) [expr $rmsd($r)/$num_steps] # print the RMSD puts "RMSD of residue $r is $rmsd($r)" puts $fil " $r \t $rmsd($r)" set res_b [atomselect $mol "resid $r"] $res_b set user $rmsd($r) $res_b delete set ave [expr $ave + $rmsd($r)] } set ave [expr $ave/[llength $res]] puts " Average rmsd per residue: $ave" close $fil } 大家看看程序的问题吗?还有，之前的任何过程都没出现任何的error。先谢谢各位的指点了。  |
回复此楼» 收录本帖的淘帖专辑推荐
 分子动力学 | MDs-Gromacs | gromacs |
» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【口语听力技巧26】和美国人学说话：那些老师教你的错误英语
已经有14人回复
- 【讨论】大家来看看我的这个黄鸣龙反应
已经有0人回复
- 【讨论】“生死”一个月
已经有2人回复
- 【讨论】各位大侠，看看我的计划如何？
已经有20人回复
- 【交流讨论】今早外导的回信让我很迷惑。。。我该往下准备么。。
已经有23人回复
- 【交流讨论】请看我的邀请函
已经有4人回复
- 【交流讨论】请大家看看我的邀请函有没有问题，还有外导写了好多字，超过一页了，会不
已经有2人回复
- 【讨论】电面之后
已经有55人回复
- 【讨论】折磨我一个月的羟醛缩合反应
已经有5人回复
- 【讨论】讨论下这人制作的“design expert 用于响应面分析”视频
已经有7人回复
- 【讨论】假期里，我在坚持学习啊，哈哈，只玩了一天。
已经有2人回复
- 【资源】施一公谈学生如何阅读和写作论文
已经有52人回复
- 【原创】BJ J-1 7.28 1000水果~~顺便问一下有童鞋知道去美国要打什么疫苗么？
已经有7人回复
- 【交流讨论】CSC 我的担心与苦恼
已经有41人回复
|
是不是我跑平衡动力学时，有些参数设计的不对啊，这是我平衡动力学的参数： ## JOB DESCRIPTION ## ############################################################# # Minimization and Equilibration of # Ubiquitin in a Water Box ############################################################# ## ADJUSTABLE PARAMETERS ## ############################################################# structure ../common/1kim_wb.psf coordinates ../common/1kim_wb.pdb set temperature 310 set outputname 1kim_wb_eq firsttimestep 0 ############################################################# ## SIMULATION PARAMETERS ## ############################################################# # Input paraTypeCharmm on parameters ../common/par_all27_prot_lipid.inp temperature $temperature # Force-Field Parameters exclude scaled1-4 1-4scaling 1.0 cutoff 20. switching on switchdist 18. pairlistdist 22. # Integrator Parameters timestep 2.0 ;# 2fs/step rigidBonds all ;# needed for 2fs steps nonbondedFreq 1 fullElectFrequency 2 stepspercycle 10 # Constant Temperature Control langevin on ;# do langevin dynamics langevinDamping 5 ;# damping coefficient (gamma) of 5/ps langevinTemp $temperature langevinHydrogen off ;# don't couple langevin bath to hydrogens # Periodic Boundary Conditions cellBasisVector1 79. 0. 0. cellBasisVector2 0. 72. 0. cellBasisVector3 0. 0. 72. cellOrigin 0. 0. 0. wrapAll on # PME (for full-system periodic electrostatics) PME yes PMEGridSizeX 82 PMEGridSizeY 75 PMEGridSizeZ 75 # Constant Pressure Control (variable volume) useGroupPressure yes ;# needed for rigidBonds useFlexibleCell no useConstantArea no langevinPiston on langevinPistonTarget 1.01325 ;# in bar -> 1 atm langevinPistonPeriod 100. langevinPistonDecay 50. langevinPistonTemp $temperature # Output outputName $outputname restartfreq 500 ;# 500steps = every 1ps dcdfreq 250 xstFreq 250 outputEnergies 100 outputPressure 100 ############################################################# ## EXTRA PARAMETERS ## ############################################################# ############################################################# ## EXECUTION SCRIPT ## ############################################################# # Minimization minimize 5000 reinitvels $temperature run 20000 ;# 5ps |
2楼2010-12-15 15:38:50
|
这个是主链的随时RMSD，使用的脚本是： set outfile [open rmsd.dat w]; set nf [molinfo top get numframes] set frame0 [atomselect top "protein and backbone and noh" frame 0] set sel [atomselect top "protein and backbone and noh"] # rmsd calculation loop for {set i 1 } {$i < $nf } { incr i } { $sel frame $i $sel move [measure fit $sel $frame0] puts $outfile "$i [measure rmsd $sel $frame0]" } close $outfile  |
3楼2010-12-15 18:04:55
4楼2010-12-15 18:26:56
5楼2010-12-15 18:41:25
6楼2010-12-15 22:33:10
7楼2010-12-15 23:12:39
30