| 5 | 1/1 | 返回列表 |
| 查看: 998 | 回復(fù): 4 | ||
亍寞金蟲 (小有名氣)
|
[求助]
求助VMD里面水的表示。。或者存儲rst7的格式。。。
|
|
想要取出來某原子半徑15A內(nèi)的所有水分子和蛋白,然后存儲成標(biāo)準(zhǔn)格式的rst7,也就是6F12.7,一行六個坐標(biāo)兩個原子的那個,但是直接存儲出來是一行十個坐標(biāo),后續(xù)處理不方便啊。。。這個有啥辦法可以存儲成6F12.7的么? 然后想著,如果不行從頭開始,于是想要用VMD的save coordinates來存pdb然后leap再處理, 但是water and within 15.0 of name C, 但是在邊緣區(qū)域的水分子不完整,只有個H或者只有個O,用leap試圖添加缺少的原子但是好像leap或者是咱。。還是不夠智能總是失敗。。。 編程實在乏力 ![]() ![]() ![]() 求助 |
|
VMD 可以的,多一條命令: set sel [atomselect top "water and name OW and within 15.0 of name C"] #這一步選出C 15埃米以內(nèi)的水分子中的氧原子 set waterID [$sel get resname] #這一步記下上一步選出的水分子的 resid set theWantedWater [atomselect top "water and resid $waterID"] # theWantedWater 就是你所需要的了,這樣就沒有殘缺的情況了 $theWantedWater writepdb theWantedWater.pdb #theWantedWater.pdb 就保存好了,繼續(xù)你剩下的工作吧 |
金蟲 (小有名氣)
|
感謝回復(fù),我之前沒用過VMD的命令行方式,不知道為什么最后存出來的pdb是空文件?我輸入文件里面水里面的氧,名字是O,于是我把第一條命令中的OW換成O,其他的沒有變化,然后執(zhí)行完了第二條命令之后,終端里面顯示了出來了WAT WAT WAT WAT。。。。感覺應(yīng)該是成功選中了,但是最后出來是空文件,不知怎解? 順便想問一下。。。您對VMD源碼熟悉么?我想直接去改http://www.ks.uiuc.edu/Research/ ... ugin_8c-source.html 這里面的8.3部分改成需要的12.7,但是以前沒改過vmd,找了半天也不知道在哪里改。。。 |
金蟲 (小有名氣)
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[有機(jī)交流]
20+3
|
FENGSHUJEI 2026-03-23 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 289求調(diào)劑 +7 | 碩星赴 2026-03-23 | 7/350 |
|
|
[考研] 08工學(xué)調(diào)劑 +10 | 用戶573181 2026-03-20 | 15/750 |
|
|
[考研] 招08考數(shù)學(xué) +7 | laoshidan 2026-03-20 | 16/800 |
|
|
[考研] 一志愿211 初試270分 求調(diào)劑 +5 | 谷雨上岸 2026-03-23 | 6/300 |
|
|
[考研] 材料專碩331求調(diào)劑 +4 | 鮮當(dāng)牛 2026-03-24 | 4/200 |
|
|
[考研] 327求調(diào)劑 +5 | prayer13 2026-03-23 | 5/250 |
|
|
[考研] 269求調(diào)劑 +4 | 我想讀研11 2026-03-23 | 4/200 |
|
|
[考研] 307求調(diào)劑 +11 | 冷笙123 2026-03-17 | 11/550 |
|
|
[考研] 尋找調(diào)劑 +4 | 倔強(qiáng)芒? 2026-03-21 | 4/200 |
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 4/200 |
|
|
[考研] 311求調(diào)劑 +3 | 勇敢的小吳 2026-03-20 | 3/150 |
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
[考研] 274求調(diào)劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
[考研] 321求調(diào)劑 +9 | 何潤采123 2026-03-18 | 11/550 |
|
|
[考研] 考研調(diào)劑求學(xué)校推薦 +3 | 伯樂29 2026-03-18 | 5/250 |
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|