| 2 | 1/1 | 返回列表 |
| 查看: 3766 | 回復(fù): 1 | ||
gn02530640銀蟲(chóng) (小有名氣)
|
[求助]
關(guān)于Material Studio 的pdb輸出檔 已有1人參與
|
|
先說(shuō)一下,我想問(wèn)的問(wèn)題其實(shí)跟分子模擬其實(shí)沒(méi)什么關(guān)系… 情況是這樣的。我個(gè)人目前使用NAMD作模擬動(dòng)力學(xué)計(jì)算,我的做法是先在Material Studio 或者是Gaussian 下創(chuàng)建我所需要的單體分子/聚合物以及作前處理(能量?jī)?yōu)化,計(jì)算電荷,etc……),甚至是用Amorphous Tools把分子放進(jìn)盒子后,再輸出成NAMD所需要的pdb文件?傊鞒淌沁@樣: 用 Gaussian 創(chuàng)建分子并且做結(jié)構(gòu)優(yōu)化和電荷計(jì)算后,輸出成pdb文件 -> 用 MS 讀取該 pdb文件,用 Amorphous Tools 把分子放進(jìn)盒子里 -> 將得到的 xtd 文件輸出成pdb文件,用于NAMD (這個(gè)過(guò)程可能不完全正確甚至有點(diǎn)麻煩,但至少我目前的確是這樣做的) 問(wèn)題來(lái)了,我從MS轉(zhuǎn)換輸出最后的pdb文件,分子的殘基序號(hào)都是相同的,舉例來(lái)說(shuō): HETATM 1 OH2 HOH W 1 HETATM 2 H1 HOH W 1 HETATM 3 H2 HOH W 1 TER 4 HETATM 5 OH2 HOH W 1 HETATM 6 H1 HOH W 1 HETATM 7 H2 HOH W 1 TER 8 HETATM 9 OH2 HOH W 1 HETATM 10 H1 HOH W 1 HETATM 11 H2 HOH W 1 TER 12 可以發(fā)現(xiàn)3個(gè)水分子的殘基序號(hào)都是1,這樣的pdb文件在產(chǎn)生psf文件的過(guò)程中肯定會(huì)出錯(cuò)。相反的,如果是在MS下直接畫(huà)好水分子并且輸出成pdb文件的話(huà),就不會(huì)有這種情況發(fā)生: ATOM 1 O1 MOL 1 ATOM 2 H2 MOL 1 ATOM 3 H3 MOL 1 ATOM 4 O4 MOL 2 ATOM 5 H5 MOL 2 ATOM 6 H6 MOL 2 ATOM 7 O7 MOL 3 ATOM 8 H8 MOL 3 ATOM 9 H9 MOL 3 TER 10 回到主題,我目前的解決方式是用文本編輯器開(kāi)啟pdb文件后手動(dòng)修改殘基序號(hào),但效率太低…只有幾個(gè)分子的話(huà)還好處理,如果是幾千個(gè)分子…… 因此我的問(wèn)題很簡(jiǎn)單:有沒(méi)有什么辦法讓MS在輸出pdb文件時(shí)自動(dòng)修改殘基序號(hào)? 或者是有沒(méi)有什么辦法可以快速修改殘基序號(hào)(比方說(shuō)寫(xiě)腳本?) 有點(diǎn)太啰嗦了,總之想請(qǐng)各位蟲(chóng)友指點(diǎn)一下迷津,感激不盡! |
鐵桿木蟲(chóng) (著名寫(xiě)手)
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0856調(diào)劑,是學(xué)校就去 +7 | sllhht 2026-03-19 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 一志愿武漢理工材料工程專(zhuān)碩調(diào)劑 +6 | Doleres 2026-03-19 | 6/300 |
|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +8 | 墨墨漠 2026-03-18 | 8/400 |
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過(guò)A區(qū)線 +8 | 步川酷紫123 2026-03-13 | 8/400 |
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 6/300 |
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 10/500 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
[考研] 287求調(diào)劑 +3 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
[考研] 0703化學(xué) 305求調(diào)劑 +4 | FY_yy 2026-03-14 | 4/200 |
|
|
[考研] 0854可跨調(diào)劑,一作一項(xiàng)核心論文五項(xiàng)專(zhuān)利,省、國(guó)級(jí)證書(shū)40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
[考博] 環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 5/250 |
|
|
[考研] 0703化學(xué)求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
[考研] 材料工程專(zhuān)碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼?~+123456 2026-03-14 | 3/150 |
|