版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

藍(lán)紫色玻璃心

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1434.5
帖子: 66
在線: 38.3小時
蟲號: 973642
注冊: 2010-03-17
專業(yè): 化學(xué)計量學(xué)與化學(xué)信息學(xué)

[求助] VMD分析殘基RMSD

看網(wǎng)上說用VMD自帶的rmsd-fullthrottle.tcl文件就可以對指定殘基計算RMSD，但是這個文件不是很明白（對程序一點都不明白

），想問一下有沒有用這個程序計算RMSD的？我對這個文件里面應(yīng)該修改什么？或者能不能給個具體步驟？一點頭緒都沒有！

# This script contains a procedure called rmsd_residue_over_time that calculates the average RMSD for each residue in a selection over all frames in a trajectory. The procedure is called as:
# rmsd_residue_over_time mol sel_resid
#where mol is the molecule in VMD and sel_resid is a list of the residue numbers in that selection.

#You can use the procedure for any residue or list of residues. Here, as an example, we will make a selection for all residues in the protein. Note that this will take a long time to calculate:

set sel_resid [[atomselect top "protein and alpha"] get resid]

#The procedure is presented below.  It also sets the B value to the value calculated, so you can color the protein by RMSD.  The call for the procedure is at the end of the file.

proc rmsd_residue_over_time {{mol top} res} {

# use frame 0 for the reference
set reference [atomselect 1 "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]

foreach r $res {
set rmsd($r) 0
}

#loop over all frames in the trajectory
for {set frame 1} {$frame < $num_steps} {incr frame} {
puts "Calculating rmsd for frame $frame ..."
# get the correct frame
$compare 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 "chain U and resid $r and noh" frame 0]
      set comp [atomselect $mol "chain U 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)"
      set res_b [atomselect $mol "resid $r"]
         $res_b set beta $rmsd($r)
         $res_b delete
      set ave [expr $ave + $rmsd($r)]

}
set ave [expr $ave/[llength $res]]
puts " Average rmsd per residue: $ave"

}

#Call the procedure

rmsd_residue_over_time top $sel_resid

回復(fù)此樓

» 猜你喜歡

085405 考的11408求各位老師帶走已經(jīng)有3人回復(fù)
339求調(diào)劑已經(jīng)有6人回復(fù)
材料與化工272求調(diào)劑已經(jīng)有4人回復(fù)
085701環(huán)境工程，267求調(diào)劑已經(jīng)有15人回復(fù)
材料求調(diào)劑一志愿哈工大324 已經(jīng)有9人回復(fù)
292求調(diào)劑已經(jīng)有15人回復(fù)
304求調(diào)劑已經(jīng)有5人回復(fù)
材料與化工考研調(diào)劑已經(jīng)有17人回復(fù)
321求調(diào)劑已經(jīng)有5人回復(fù)
329求調(diào)劑已經(jīng)有5人回復(fù)

1樓 2012-03-16 09:46:38

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

智能機器人

Robot (super robot)

我們都愛小木蟲

找到一些相關(guān)的精華帖子，希望有用哦~

鄰硝基芐負(fù)離子已經(jīng)有5人回復(fù)
3-甲基-2噻吩甲醛的2位上溴問題已經(jīng)有15人回復(fù)
轉(zhuǎn)基因植物安全嗎？已經(jīng)有12人回復(fù)
【求助】關(guān)于BOC保護(hù)胺基的反應(yīng)后處理問題�。�！已經(jīng)有10人回復(fù)
青基參與人有無年齡限制已經(jīng)有4人回復(fù)
含有芐溴的硝基苯的還原問題已經(jīng)有10人回復(fù)
急問：青基，參考文獻(xiàn)要實行實引制嗎？即文中提到的文獻(xiàn)與后面的參考文獻(xiàn)一一對應(yīng)嗎？已經(jīng)有9人回復(fù)
青基排版格式的-括號內(nèi)的內(nèi)容刪除了沒？急啊已經(jīng)有4人回復(fù)
液相的基線一直波動是怎么回事已經(jīng)有9人回復(fù)
醛轉(zhuǎn)化成氰基能一步到位嗎已經(jīng)有3人回復(fù)
聚合物以及聚合物基復(fù)合材料斷裂韌性GIC或者KIC測試已經(jīng)有6人回復(fù)
轉(zhuǎn)基因植株的RT-PCR結(jié)果分析已經(jīng)有8人回復(fù)
國青基申請人介紹，沒有留參加人的介紹的框啊已經(jīng)有4人回復(fù)
【求助】最終RMSD結(jié)果終于分析出來了，但是很不好。已經(jīng)有5人回復(fù)
【求助】DLPOLY計算出來的結(jié)果，是否要編程才能分析？已經(jīng)有15人回復(fù)
【求助】（請教！�。┤绾卧赩MD中對結(jié)構(gòu)圖形中的特定殘基進(jìn)行標(biāo)記？？已經(jīng)有9人回復(fù)

點擊這里搜索更多相關(guān)資源

科研從小木蟲開始，人人為我，我為人人

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主藍(lán)紫色玻璃心的主題更新