| 2 | 1/1 | 返回列表 |
| 查看: 932 | 回復: 1 | ||
[求助]
namd模擬中重原子的固定
|
|
新人求助: 初學分子模擬,老師教namd,給的任務是跑完2YHC這個蛋白,可是這個蛋白晶體結構有殘缺,從102-115殘基都缺,老師又讓先學modeller補全,補好之后要用namd對補出的那段(102-115)進行優(yōu)化,分別有兩個配置文件NVT.conf 和NPT.conf,我截取了101-116進行優(yōu)化,先固定重原子(occupancy列重原子設置成1,氫原子設置0), 依次跑了NVT和NPT,之后放開重原子(配置文件里的一個選項設置成off),又依次跑了NPT和NVT,四步過程中101和116的所有原子(包括氫原子)的occupancy列我都設置是1(我想的是這樣模擬完成之后101和116的坐標都沒變很好回拼),這樣在最后的一步NVT模擬后我將得到的dcd文件轉化成pdb文件希望能拼回去可是發(fā)現(xiàn)101和116殘基的坐標已經變了,老師說要經過平移旋轉等操作接回去,現(xiàn)在我技術還沒掌握,不過我卻對模擬操作有疑問。 模擬操作第三步和第四步中要放開重原子,然后要把conf文件中原子約束行設置成off,這樣我就糾結了,那我把101和116的所有原子都固定有什么意義呢,原子約束行設置成off后就不識別pdb文件中occupancy列的數(shù)字了吧,namd不會理解我想固定101和116的想法,這樣得到的最終結果中101和116的坐標肯定會變。然后老師說可以通過平移旋轉等接回去,我看了下101和116的坐標,他們和原坐標不具備簡單的統(tǒng)一平移關系,如果硬要把坐標變回去可能會破壞優(yōu)化的結果,這樣優(yōu)化就沒有意義了,所以我懷疑優(yōu)化過程中約束重原子部分的操作有問題,另一方面就是可能關于回拼的技術平移旋轉部分有深意我理解不夠,希望得到大家指點。 本人namd新手,老師給個方向自己摸索,理解不對之處多多,前段時間也發(fā)過類似的帖子這次是想印證下自己想法,說不清楚所以寫的略多。 本人新蟲,金幣不多,希望有想法的前輩不嫌棄給個回復。謝謝啦。 |
木蟲 (正式寫手)
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 266求調劑 +11 | 陽陽哇塞 2026-03-27 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 0703一志愿9,初試成績:338,四六級已過,有科研經歷,求調劑! +3 | Zuhui0306 2026-03-25 | 3/150 |
|
|
[考研] 339求調劑 +4 | 烤麥芽 2026-03-27 | 5/250 |
|
|
[考研] 348求調劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
[考研] 0856材料化工調劑 總分330 +8 | zhubinhao 2026-03-27 | 8/400 |
|
|
[考研] 求調劑 +6 | 林之夕 2026-03-24 | 6/300 |
|
|
[考研] 求調劑,一志愿 南京航空航天大學大學 ,080500材料科學與工程學碩 +4 | @taotao 2026-03-26 | 5/250 |
|
|
[考研] 生物學 296 求調劑 +4 | 朵朵- 2026-03-26 | 6/300 |
|
|
[考研] 340求調劑 +3 | Amber00 2026-03-26 | 3/150 |
|
|
[考研] 085601求調劑總分293英一數(shù)二 +4 | 鋼鐵大炮 2026-03-24 | 4/200 |
|
|
[考研] 0856求調劑 +8 | zhn03 2026-03-25 | 9/450 |
|
|
[考研] 一志愿中南大學化學學碩0703總分337求調劑 +7 | niko- 2026-03-22 | 7/350 |
|
|
[考研] 各位老師您好:本人初試372分 +5 | jj涌77 2026-03-25 | 6/300 |
|
|
[考研] 340求調劑 +5 | 話梅糖111 2026-03-24 | 5/250 |
|
|
[考研] 306求0703調劑一志愿華中師范 +10 | 紙魚ly 2026-03-21 | 11/550 |
|
|
[考研] 335求調劑 +4 | yuyu宇 2026-03-23 | 5/250 |
|
|
[考研] 284求調劑 +3 | yanzhixue111 2026-03-23 | 6/300 |
|
|
[考研] 336化工調劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|
|
[考研] 308求調劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|
[考研] 0703化學調劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|