 |
 |
 |
| 24小時(shí)熱門版塊排行榜 |
- >論壇更新日志 (1101)
- >考研 (175)
- >導(dǎo)師招生 (106)
- >蟲(chóng)友互識(shí) (47)
- >考博 (26)
- >文獻(xiàn)求助 (24)
- >論文投稿 (22)
- >博后之家 (15)
- >碩博家園 (15)
- >文獻(xiàn)求助完結(jié)子版 (14)
- >找工作 (14)
- >基金申請(qǐng) (9)
- >論文道賀祈福 (8)
- >教師之家 (5)
- >職場(chǎng)人生 (4)
- >休閑灌水 (4)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 11004 | 回復(fù): 16 | ||||||||||||||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||||||||||||||
御劍江湖榮譽(yù)版主 (著名寫(xiě)手)
天池冶海
|
[交流]
【轉(zhuǎn)帖】用Amber+Gaussian做小分子力場(chǎng) 已有14人參與
|
|||||||||||||||
|
用Amber+Gaussian做小分子力場(chǎng)，不是我的研究領(lǐng)域，不過(guò)閱讀此篇文章，發(fā)現(xiàn)不錯(cuò)，轉(zhuǎn)載給此研究領(lǐng)域的蟲(chóng)子們。。。 第一步：生成小分子模板 蛋白質(zhì)中各氨基酸殘基的力參數(shù)是預(yù)先存在的，但是很多模擬過(guò)程會(huì)涉及配體分子，這些有機(jī)小分子有很高的多樣性，他們的力參數(shù)和靜電信息不可能預(yù)存在庫(kù)文件中，需要根據(jù)需要自己計(jì)算生成模板。amber中的antechamber 程序就是生成小分子模板的。生成模板要進(jìn)行量子化學(xué)計(jì)算，這一步可以由antechamber中附帶的mopac完成，也可以由gaussian完成，這里介紹用gaussian計(jì)算的過(guò)程。 建議在計(jì)算前用sybyl軟件將小分子預(yù)先優(yōu)化，不要用gaussian優(yōu)化，大基組從頭計(jì)算進(jìn)行幾何優(yōu)化花費(fèi)的計(jì)算時(shí)間太長(zhǎng)。gaussian計(jì)算的輸入文件可以用antechamber程序直接生成，生成后去掉其中關(guān)于幾何優(yōu)化的參數(shù)即可將小分子優(yōu)化后的結(jié)構(gòu)存儲(chǔ)為mol2各式，上傳到工作目錄，用antechamber程序生成gaussian輸入文件，命令如下： antechamber -i 49.mol2 -fi mol2 -o 49.in -fo gzmat 這樣可以生成49.in文件，下載到windows環(huán)境，運(yùn)行g(shù)aussian計(jì)算這個(gè)文件，如果發(fā)現(xiàn)計(jì)算時(shí)間過(guò)長(zhǎng)或者內(nèi)存不足計(jì)算中斷，可以修改文件選擇小一些的基組。 獲得輸出文件49.out之后將它上傳到工作目錄，再用antechamber生成模板，命令如下： antechamber -i 49.out -fi gout -o 49mod.mol2 -fo mol2 -c resp 運(yùn)行之后就會(huì)生成一個(gè)新的mol2文件，如果用看圖軟件打開(kāi)這個(gè)文件會(huì)發(fā)現(xiàn)，原子的顏色很怪異，這是因?yàn)閙ol2的原子名稱不是標(biāo)準(zhǔn)的原子名稱，看圖軟件無(wú)法識(shí)別。下面一步是檢查參數(shù)，因?yàn)榭赡軙?huì)有一些特殊的參數(shù)在gaff中不存在需要程序注入， 命令如下： parmchk -i 49mod.mol2 -f mol2 -o 49mod 這樣那些特殊的參數(shù)就存在49mod這個(gè)文件中了 第二步：處理蛋白質(zhì)文件 amber自帶的leap程序是處理蛋白質(zhì)文件的，他可以讀入PDB各式的蛋白質(zhì)文件，根據(jù)已有的力場(chǎng)模板為蛋白質(zhì)賦予鍵參數(shù)和靜電參數(shù)。 PDB格式的文件有時(shí)會(huì)帶有氫原子和孤對(duì)電子的信息，但是在這種格式下氫原子和孤對(duì)電子的命名不是標(biāo)準(zhǔn)命名，力場(chǎng)模板無(wú)法識(shí)別這種不標(biāo)準(zhǔn)的命名，因此需要將兩者的信息刪除 ATOM 12 1H ARG A 82 12.412 8.891 34.128 1.00 0.00 H 在PDB各式里面，氫原子的信息會(huì)在第13或者14列出現(xiàn)H字符，可以應(yīng)用grep命令刪除在13或者14列出現(xiàn)H的行命令如下： grep -v ‘^………….H’ 1t4j.pdb > x grep -v ‘^…………H’ x > 1t4j_noh.pdb 除了刪除氫和孤對(duì)電子，還應(yīng)該把文件中的結(jié)晶水、乙酸等分子刪除，這些分子的信息常常集中在文件的尾部，可以直接刪除。 處理過(guò)之后的蛋白質(zhì)文件，只包括各氨基酸殘基和小分子配體的重原子信息，模擬需要的氫原子和水分子將在leap中添加接下來(lái)需要進(jìn)一步整理蛋白質(zhì)文件，主要關(guān)注氨基酸的不同存在形態(tài)和小分子的原子名稱。半胱氨酸有兩種存在形態(tài)，有的是兩個(gè)半胱氨酸通過(guò)二硫鍵相連，有的則是自由的，兩種半胱氨酸在力場(chǎng)文件中是用不同的unit來(lái)表示的，這相當(dāng)于是兩個(gè)完全不同的氨基酸，需要手動(dòng)更改蛋白質(zhì)文件中半胱氨酸的名字，橋連的要用CYX，自由的用CYS 組氨酸有若干種質(zhì)子態(tài)，和半胱氨酸一樣，也需要查閱文獻(xiàn)確定這些質(zhì)子態(tài)，并更改殘基名稱 最后需要修改的是配體分子的原子名，這是工作量最大的事情，仔細(xì)觀察可以發(fā)現(xiàn)，一般PDB文件中配體的各個(gè)原子名稱，和我們上面通過(guò)antechamber 生成的49mod.mol2中原子名稱并不一致，這會(huì)造成leap無(wú)法識(shí)別這些原子，無(wú)法為這些原子賦予力參數(shù)和靜電參數(shù)， 因此需要手動(dòng)更改蛋白質(zhì)文件中配體分子的原子名稱。 進(jìn)行這一步可以同時(shí)用看圖軟件打開(kāi)49mod.mol2和蛋白質(zhì)文件，隱藏蛋白質(zhì)文件中除配體分子以外的所有結(jié)構(gòu)，旋轉(zhuǎn)兩個(gè)文件，讓他們姿態(tài)相近，以方便觀察，并且在各自均標(biāo)識(shí)原子名稱，然后用文字編輯軟件打開(kāi)蛋白質(zhì)文件，核對(duì)看圖軟件中兩個(gè)分子對(duì)應(yīng)的原子名稱，手動(dòng)逐一修改。 修改原子名稱需要極大的耐心和細(xì)心，如果發(fā)生錯(cuò)誤下一步的操作會(huì)無(wú)法繼續(xù)。我現(xiàn)在想到的也只有這個(gè)笨辦法，如果誰(shuí)還有別的好法子，歡迎告訴我。 現(xiàn)在文件的準(zhǔn)備工作都已經(jīng)完成，該開(kāi)始正式的模擬了 第三步：生成拓?fù)湮募妥鴺?biāo)文件 用amber進(jìn)行分子動(dòng)力學(xué)模擬需要坐標(biāo)和拓?fù)湮募�，坐�?biāo)文件記錄了各個(gè)質(zhì)點(diǎn)所座落的坐標(biāo)，拓?fù)湮募涗浟苏麄�(gè)體系各質(zhì)點(diǎn)之間的鏈接狀況、力參數(shù)電荷等信息。這兩個(gè)文件是由leap 程序生成的amber中有兩個(gè)leap程序，一個(gè)是純文字界面的tleap，一個(gè)是帶有圖形界面的Xleap。但是amber的圖形界面做得很差，用遠(yuǎn)程登錄使用圖形界面不僅麻煩而且慢，所以我比較偏愛(ài)使用tleap，兩個(gè)leap的命令是完全一樣的，其實(shí)用哪一個(gè)都無(wú)所謂。 啟動(dòng)tleap，在shell里輸入命令tleap，leap就啟動(dòng)了，shell里會(huì)顯示 -I: Adding /usr/local/amber8/antechamber-1.23/dat/leap/prep to search path. -I: Adding /usr/local/amber8/antechamber-1.23/dat/leap/lib to search path. -I: Adding /usr/local/amber8/antechamber-1.23/dat/leap/parm to search path. -I: Adding /usr/local/amber8/antechamber-1.23/dat/leap/cmd to search path. Welcome to LEaP! (no leaprc in search path) > 這個(gè)>是leap的提示符 下面要調(diào)入庫(kù)文件。amber是模擬生物分子的好手，主要就是依靠專門為蛋白質(zhì)多糖核酸量身訂做的amber力場(chǎng)，力場(chǎng)的所有參數(shù)都存 儲(chǔ)在庫(kù)文件里，所以打開(kāi)leap第一件事便是調(diào)入庫(kù)文件。 amber提供了很多種庫(kù)，這里我們只用到兩個(gè)庫(kù)，gaff和02庫(kù)，輸入命令： >source leaprc.gaff >source leaprc.ff02 之后兩個(gè)庫(kù)就調(diào)入進(jìn)來(lái)了 這時(shí)可以用list命令看看庫(kù)里都有什么： 這里面羅列的就是庫(kù)里面的unit，包括20種氨基酸、糖以及核酸還有一些常見(jiàn)離子的參數(shù) 下面一步是調(diào)入配體分子的模板，首先補(bǔ)全參數(shù)，輸入命令： >loadamberparams 49mod 然后讀入模板文件,輸入命令： >MOL = loadmol2 49mod.mol2 其中MOL是unit的名字，要保證這個(gè)名字和pdb文件中配體的殘基名完全一致，否則系統(tǒng)仍然無(wú)法識(shí)別pdb文件中的小分子 現(xiàn)在再輸入list命令會(huì)發(fā)現(xiàn)庫(kù)里面多了一個(gè)unit: 那個(gè)就是配體分子的模板 下面讀入pdb文件，輸入命令： >comp = loadpdb 1t4j_noh.pdb 如果輸入這個(gè)命令之后，屏幕上出現(xiàn)了大量的創(chuàng)建unit或者atom的信息，如下所示，則說(shuō)明上面一步的pdb文件處理沒(méi)有做好， 還需要重新處理，通常這種情況都發(fā)生在配體分子上面，有時(shí)則是因?yàn)榈鞍踪|(zhì)中存在特殊殘基。 Creating new UNIT for residue: FRJ sequence: 1 Created a new atom named: O36 within residue: .R Created a new atom named: S33 within residue: .R Created a new atom named: O35 within residue: .R Created a new atom named: N34 within residue: .R 如果屏幕僅僅顯示添加原子，這說(shuō)明輸入的PDB文件中缺失了部分重原子，leap根據(jù)模板自動(dòng)補(bǔ)全了這些缺失的原子，這種情況 不會(huì)影響進(jìn)一步的計(jì)算 Added missing heavy atom: .R.A Added missing heavy atom: .R.A Added missing heavy atom: .R.A Added missing heavy atom: .R.A 根據(jù)體系的具體情況，還可能要將成對(duì)的CYX殘基中的二硫鍵相連，有時(shí)候還會(huì)鏈接其他原子，比如將糖基鏈接在氨基酸殘基上 ，用bond命令可以完成，命令用法如下： >bond comp.35.SG comp.179.SG 其中comp是剛才讀入的分子名稱，35和179是殘基序號(hào)，SG是CYX殘基模板中硫原子的名稱，用comp.35.SG這樣的語(yǔ)法就可以定 位一個(gè)原子 如果希望進(jìn)行考慮溶劑效應(yīng)的動(dòng)力學(xué)模擬，可能還需要為體系加上水，加水有很多種命令，這里只列舉一個(gè)： >solvatebox comp TIP3PBOX 10.0 solvatebox命令是說(shuō)要加上一個(gè)方形的周期水箱，comp指要加水的分子，TIP3PBOX是選擇的水模板名稱，10.0是水箱子的半徑 有的體系總電荷不為0，為了模擬穩(wěn)定，需要加入抗衡離子，系統(tǒng)會(huì)自動(dòng)計(jì)算體系的靜電場(chǎng)分布，在合適的位置加上離子，命令如下： >addions comp Na+ 0 意思是用鈉離子把體系總電荷補(bǔ)平，還可以選擇其他庫(kù)里面有的離子。 完成這一步之后就可以輸出拓?fù)湮募妥鴺?biāo)文件了，但是為了方便起見(jiàn)，在運(yùn)行最后一步之前要先把leap里加工好的分子單獨(dú)存 成一個(gè)庫(kù)文件，以后可以直接調(diào)用這個(gè)庫(kù)文件，免得重復(fù)上面的操作： >saveoff comp 1taj.off 這樣就生成了一個(gè)off文件在那里，下面輸出拓?fù)湮募妥鴺?biāo)文件 >saveamberparm comp 1t4j.prmtop 1t4j.inpcrd Checking Unit. Building topology. Building atom parameters. Building bond parameters. Building angle parameters. Building proper torsion parameters. Building improper torsion parameters. total 1 improper torsion applied Building H-Bond parameters. Not Marking per-residue atom chain types. Marking per-residue atom chain types. (Residues lacking connect0/connect1 - these don’t have chain types marked: res total affected CMET 1 ) (no restraints) >quit 現(xiàn)在準(zhǔn)備好拓?fù)湮募妥鴺?biāo)文件，接下來(lái)就可以開(kāi)始能量?jī)?yōu)化和動(dòng)力學(xué)模擬了。如果愿意的話還可以用ambpdb這個(gè)命令生成一 個(gè)pdb文件，直觀地看一看生成了什么東西，命令如下： [snowyowls@localhost actualamber]$ ambpdb -p 1t4j.prmtop <1t4j.inpcrd> kankan.pdb | New format PARM file being parsed. | Version = 1.000 Date = 09/08/06 Time = 16:36:09 [snowyowls@localhost actualamber]$ 可以下載之后用看圖軟件欣賞，如果加了溶劑盒子的話，看的時(shí)候可要小心，會(huì)恨嚇人的樣子。 第四步：能量?jī)?yōu)化 用amber進(jìn)行分子動(dòng)力學(xué)模擬需要坐標(biāo)和拓?fù)湮募�，這在上一步已經(jīng)完成了，分別生成了1t4j.prmtop 和1t4j.inpcrd兩個(gè)文件 ，下面一步就是動(dòng)力學(xué)模擬之前的能量?jī)?yōu)化了。 由于我們進(jìn)行的起始結(jié)構(gòu)來(lái)自晶體結(jié)構(gòu)或者同源模建，所以在分子內(nèi)部存在著一定的張力，能量?jī)?yōu)化就是在真正的動(dòng)力學(xué)之前， 釋放這些張力，如果沒(méi)有這個(gè)步驟，在動(dòng)力學(xué)模擬開(kāi)始之后，整個(gè)分子可能會(huì)因此散架。 能量?jī)?yōu)化由sander模塊完成，運(yùn)行sander至少需要三個(gè)輸入文件，分別是分子的拓?fù)湮募�，坐�?biāo)文件以及sander的控制文件。 現(xiàn)在分子的拓?fù)湮募妥鴺?biāo)文件已經(jīng)完成，需要建立輸入文件，min_1.in Initial minimisation of our structures &cntrl imin=1, maxcyc=4000, ncyc=2000, cut=10, ntb=1,ntr=1, restraint_wt=0.5 restraintmask=’:1-283′ / 文件首行是說(shuō)明，說(shuō)明這項(xiàng)任務(wù)的基本情況； &cntrl與/之間的部分是模擬的參數(shù) 其中imin=1表示任務(wù)是能量?jī)?yōu)化，maxcyc=4000表示能量?jī)?yōu)化共進(jìn)行4000步，ncyc=2000表示在整個(gè)能量?jī)?yōu)化的4000步中， 前 2000步采用最陡下降法，在2000步之后轉(zhuǎn)換為共軛梯度法，如果模擬的時(shí)候不希望進(jìn)行方法的轉(zhuǎn)換，可以再加入另一 個(gè)關(guān)鍵詞NTMIN，如果NTMIN =0則全程使用共軛梯度法，NTMIN=2則全程使用最陡下降法，此外還有=3和=4的選項(xiàng)，分別是xmin法和lmod法 ，具體情況可以看手冊(cè)。 第二行的cut=10表示非鍵相互作用的截?cái)嘀�，單位是埃�?ntb=1表示使用周期邊界條件，這個(gè)選項(xiàng)要和前面生成的拓?fù)湮募鴺?biāo)文件相匹配， 如果前面加溶劑時(shí)候用的是盒子水，就設(shè)置ntb=1，如果加的是層水，那就應(yīng)該選擇ntb=0；ntr=1表示在能量?jī)?yōu)化的過(guò)程中要約束一些原子， 約束的是哪些原子呢？后面有。 第三行和第四行都是關(guān)于約束原子的信息，restraint_wt=0.5限定了約束的力常數(shù)，在這里約束原子就是把原子用一根彈簧拉在固定的 位置上，一旦原子偏離固定的位置，系統(tǒng)就會(huì)給他施加一個(gè)回復(fù)力，偏離的越遠(yuǎn)，回復(fù)力越大，回復(fù)力就是由這個(gè)力常數(shù)決定的，單位是 Kcal/(mol*A)。 restraintmask=’:1-283′表示約束的是從1到283號(hào)殘基，在這個(gè)分子中，1-283號(hào)殘基是蛋白質(zhì)上的氨基酸殘基，從284 號(hào)開(kāi)始，就都是水了，所以這個(gè)命令的意思就是，約束整個(gè)蛋白質(zhì)，自由地優(yōu)化溶劑分子。因?yàn)槿軇┓肿邮乔懊娴膖leap自動(dòng)給加上的， 所以一定要額外多關(guān)注一些。 下面運(yùn)行sander來(lái)執(zhí)行這個(gè)優(yōu)化： [snowyowls@localhost actualamber]$ sander -O -i min_1.in -p 1t4j.prmtop -c 1t4j.i npcrd -ref 1t4j.inpcrd -r 1t4j_min1.rst -o 1t4j_min1.out 命令中，-O表示覆蓋所有同名文件，-i min_1.in表示sander的控制文件是min_1.in，-p 1t4j.prmtop表示分子的拓?fù)湮募?br /> ，-c 1t4j.inpcrd表示坐標(biāo)文件，-ref 1t4j.inpcrd是參考坐標(biāo)文件，只有在控制文件中出現(xiàn)關(guān)鍵詞ntr=1的時(shí)候才需要給 sander指定-ref文件，這是約束原子的參考坐標(biāo)，- ref 1t4j.inpcrd就是說(shuō)以1t4j.inpcrd中的坐標(biāo)為準(zhǔn)進(jìn)行約束原子的優(yōu)化。 以上這四個(gè)是輸入文件。-r 1t4j_min1.rst表示經(jīng)過(guò)模擬之后新的原子坐標(biāo)會(huì)輸出到1t4j_min1.rst文件中，-o 1t4j_min1.out 則表示優(yōu)化過(guò)程中的相關(guān)信息都會(huì)寫(xiě)入到1t4j_min1.out文件中。 運(yùn)行起這個(gè)命令之后，等拿到 1t4j_min1.rst文件就意味著對(duì)溶劑的優(yōu)化已經(jīng)差不多了，顯然下面還需要對(duì)蛋白質(zhì)本身進(jìn)行優(yōu)化， 這個(gè)優(yōu)化還要分兩步進(jìn)行，控制文件分別是min_2.in 和min_3.in min_2.in Initial minimisation of our structures &cntrl imin=1, maxcyc=5000, ncyc=2500, cut=10, ntb=1,ntr=1, restraint_wt=0.5 restraintmask=’:1-283@CA,N,C’ / 在這里發(fā)生變化的是約束原子的范圍， ‘:1-283@CA,N,C’表示1-283號(hào)殘基中名叫CA，N和C的原子，這些原子實(shí)際上是蛋白質(zhì)主鏈上 的原子，也就是說(shuō)這一次的優(yōu)化是約束了蛋白質(zhì)主鏈上的原子之后，對(duì)溶劑和側(cè)鏈原子進(jìn)行自由優(yōu)化。 min_3.in Initial minimisation of our structures &cntrl imin=1, maxcyc=10000, ncyc=5000, cut=10, ntb=1, / 在這里不再進(jìn)行約束原子的優(yōu)化了，對(duì)整個(gè)分子進(jìn)行全原子優(yōu)化。 三步優(yōu)化的命令如下： [snowyowls@localhost actualamber]$ sander -O -i min_1.in -p 1t4j.prmtop -c 1t4j.inpcrd -ref 1t4j.inpcrd -r 1t4j_min1.rst -o 1t4j_min1.out [snowyowls@localhost actualamber]$ sander -O -i min_2.in -p 1t4j.prmtop -c 1t4j_min1.rst -ref 1t4j_min1.rst -r 1t4j_min2.rst -o 1t4j_min2.out [snowyowls@localhost actualamber]$ sander -O -i min_3.in -p 1t4j.prmtop -c 1t4j_min2.rst -r 1t4j_heat0.rst -o 1t4j_min3.out 接下來(lái)就是真正的分子動(dòng)力學(xué)模擬了…… |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 物理前沿 | xuexijisuan | 分子模擬要術(shù) | amber |
| 分子模擬 | 精華帖轉(zhuǎn)載 | 分子對(duì)接 | Amber |
| 計(jì)算化學(xué) | 搜集 | 申請(qǐng) | 精華 |
| 分子動(dòng)力學(xué)模擬 |
» 本帖已獲得的紅花（最新10朵）
wl507586» 猜你喜歡
一志愿大連理工大學(xué)材料求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料與化工調(diào)劑一志愿大連海事085600，349
已經(jīng)有6人回復(fù)
-
一志愿中海洋320化學(xué)工程與技術(shù)學(xué)碩求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料科學(xué)與工程調(diào)劑
已經(jīng)有11人回復(fù)
-
311求調(diào)劑一志愿合肥工業(yè)大學(xué)
已經(jīng)有4人回復(fù)
-
一志愿食品科學(xué)與工程083200求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料科學(xué)與工程求調(diào)劑
已經(jīng)有10人回復(fù)
-
一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑
已經(jīng)有11人回復(fù)
-
一志愿 南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩
已經(jīng)有4人回復(fù)
-
085600材料與化工調(diào)劑
已經(jīng)有15人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 怎么處理UNIT MISSING 的錯(cuò)誤？
已經(jīng)有4人回復(fù)
- 分子力場(chǎng)中非鍵能
已經(jīng)有4人回復(fù)
- 關(guān)于pdb庫(kù)使用問(wèn)題
已經(jīng)有8人回復(fù)
- 關(guān)于MS力場(chǎng)編輯
已經(jīng)有11人回復(fù)
- Amber 分子動(dòng)力學(xué) 求助！�。�！
已經(jīng)有10人回復(fù)
- 求計(jì)算一個(gè)簡(jiǎn)單小分子的尺寸-用Gaussian
已經(jīng)有4人回復(fù)
- Gromacs中添加CHARMM力場(chǎng)
已經(jīng)有41人回復(fù)
- 使用amber力場(chǎng)，一些參數(shù)無(wú)法獲得
已經(jīng)有5人回復(fù)
- 高斯產(chǎn)生小分子參數(shù)問(wèn)題
已經(jīng)有6人回復(fù)
- 關(guān)于用amber跑REMD計(jì)算的流程及輸入文件的設(shè)置
已經(jīng)有6人回復(fù)
- 如何將sybyl對(duì)接后的蛋白質(zhì)和配體保存成一個(gè)PDB文件
已經(jīng)有8人回復(fù)
- 怎么自己做力場(chǎng)-及氣體力場(chǎng)的選擇
已經(jīng)有6人回復(fù)
- 一個(gè)橋連H的問(wèn)題
已經(jīng)有6人回復(fù)
- ATP分子力場(chǎng)參數(shù)怎么做
已經(jīng)有4人回復(fù)
- GAUSSIAN計(jì)算語(yǔ)法錯(cuò)誤
已經(jīng)有7人回復(fù)
- 糖蛋白---分子動(dòng)力學(xué)模擬
已經(jīng)有4人回復(fù)
- pdb文件中加入氫原子
已經(jīng)有7人回復(fù)
- 一個(gè)模擬體系中能否用兩個(gè)不同的力場(chǎng)？
已經(jīng)有7人回復(fù)
- Gaussian模擬SERS
已經(jīng)有6人回復(fù)
- 【求助】Amber中水分子的問(wèn)題
已經(jīng)有4人回復(fù)
- 【討論】gaussian優(yōu)化過(guò)程中具體在執(zhí)行什么操作
已經(jīng)有3人回復(fù)
- 【求助】量子化學(xué)方法求amber力場(chǎng)參數(shù)
已經(jīng)有10人回復(fù)
- 【討論】Gaussian 09調(diào)用內(nèi)存的問(wèn)題
已經(jīng)有3人回復(fù)
4樓2011-04-10 23:10:46
xd200620940
銅蟲(chóng) (正式寫(xiě)手)
- 模擬EPI: 2
- 應(yīng)助: 0 (幼兒園)
- 金幣: 753.9
- 紅花: 5
- 帖子: 508
- 在線: 137.7小時(shí)
- 蟲(chóng)號(hào): 494350
- 注冊(cè): 2008-01-13
- 專業(yè): 化學(xué)計(jì)量學(xué)與化學(xué)信息學(xué)
★ ★ ★ ★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖交流
御劍江湖(金幣+3): 謝謝！ 2011-04-10 11:28:01
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖交流
御劍江湖(金幣+3): 謝謝！ 2011-04-10 11:28:01
|
最后需要修改的是配體分子的原子名，這是工作量最大的事情，仔細(xì)觀察可以發(fā)現(xiàn)，一般PDB文件中配體的各個(gè)原子名稱，和我們上面通過(guò)antechamber 生成的49mod.mol2中原子名稱并不一致，這會(huì)造成leap無(wú)法識(shí)別這些原子，無(wú)法為這些原子賦予力參數(shù)和靜電參數(shù)， 因此需要手動(dòng)更改蛋白質(zhì)文件中配體分子的原子名稱。 進(jìn)行這一步可以同時(shí)用看圖軟件打開(kāi)49mod.mol2和蛋白質(zhì)文件，隱藏蛋白質(zhì)文件中除配體分子以外的所有結(jié)構(gòu)，旋轉(zhuǎn)兩個(gè)文件，讓他們姿態(tài)相近，以方便觀察，并且在各自均標(biāo)識(shí)原子名稱，然后用文字編輯軟件打開(kāi)蛋白質(zhì)文件，核對(duì)看圖軟件中兩個(gè)分子對(duì)應(yīng)的原子名稱，手動(dòng)逐一修改。 這個(gè)是用antechamber做小分子力場(chǎng)最為困難的地方，相當(dāng)?shù)牟蝗诵曰? |
2樓2011-04-10 11:25:22
3樓2011-04-10 18:00:41
xd200620940
銅蟲(chóng) (正式寫(xiě)手)
- 模擬EPI: 2
- 應(yīng)助: 0 (幼兒園)
- 金幣: 753.9
- 紅花: 5
- 帖子: 508
- 在線: 137.7小時(shí)
- 蟲(chóng)號(hào): 494350
- 注冊(cè): 2008-01-13
- 專業(yè): 化學(xué)計(jì)量學(xué)與化學(xué)信息學(xué)
★ ★ ★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖交流
御劍江湖(金幣+2): 期待專帖詳細(xì)介紹 2011-04-12 12:15:25
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖交流
御劍江湖(金幣+2): 期待專帖詳細(xì)介紹 2011-04-12 12:15:25
|
因?yàn)橹爸挥肎ROMACS，對(duì)AMBER了解不多。最近GROMACS將AMBER力場(chǎng)包含了進(jìn)來(lái)，對(duì)于配體小分子的轉(zhuǎn)化在網(wǎng)站的教程中也給出了一個(gè)方案（但沒(méi)有實(shí)例）。 根據(jù)該教程提到的方案，可以采用ACPYPE + ANTECHAMBER得到GROMACS格式的配體分子topology和GRO，因此上周我對(duì)該問(wèn)題進(jìn)行了系統(tǒng)地摸索。得到幾點(diǎn)體會(huì)： 1.antechamber做小分子參數(shù)時(shí)轉(zhuǎn)化成功與否與原子類型無(wú)關(guān)； 2.antechamber轉(zhuǎn)化小分子參數(shù)時(shí)要求PDB文件是氫飽和的，即從PDB庫(kù)中直接得到的小分子需加氫，才能進(jìn)行參數(shù)計(jì)算； 3.若是采用繪圖軟件得到的PDB文件，即使氫已飽和，直接使用antechamber做參數(shù)也不會(huì)成功，卡在leap階段。經(jīng)過(guò)研讀教程，再比較pdb文件，找到了問(wèn)題的關(guān)鍵。也得出第1點(diǎn)的結(jié)論。 關(guān)于antechamber以及acpype具體的使用步驟，近兩天我將開(kāi)一個(gè)專帖詳細(xì)介紹，敬請(qǐng)蟲(chóng)友們關(guān)注，謝謝！ [ Last edited by xd200620940 on 2011-4-11 at 19:32 ] |
5樓2011-04-11 19:30:47
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[有機(jī)交流]
 甲基亞磺磺酸鈉和甲基磺酸酯反應(yīng)機(jī)理
 10+3
|
kaobao456 2026-03-29 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 277跪求調(diào)劑 +8 | 1915668 2026-03-27 | 12/600 |
|
|
|
[考研] 085600 295分求調(diào)劑 +7 | W55j 2026-03-30 | 9/450 |
|
|
|
[考研] 310求調(diào)劑 +10 | 爭(zhēng)取九點(diǎn)睡 2026-03-30 | 10/500 |
|
|
|
[考研] 298求調(diào)劑 +3 | 什么是胖頭魚(yú) 2026-03-30 | 5/250 |
|
|
|
[基金申請(qǐng)] 面上5B能上會(huì)嗎？ +7 | redcom 2026-03-29 | 7/350 |
|
|
|
[考研] 282求調(diào)劑 +4 | wcq131415 2026-03-24 | 4/200 |
|
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +4 | 邱gl 2026-03-26 | 7/350 |
|
|
|
[考研] 275求調(diào)劑 +15 | Micky11223 2026-03-25 | 20/1000 |
|
|
|
[考研] 356求調(diào)劑 +4 | gysy?s?a 2026-03-28 | 4/200 |
|
|
|
[考研] 298求調(diào)劑 +4 | 種圣賜 2026-03-28 | 4/200 |
|
|
|
[考研] 調(diào)劑考研 +3 | 王杰一 2026-03-29 | 3/150 |
|
|
|
[考研] 壓國(guó)家一區(qū)線，求導(dǎo)師收留，有恩必謝！ +7 | 迷人的哈哈 2026-03-28 | 7/350 |
|
|
|
[考研] 材料277求調(diào)劑 +7 | min3 2026-03-24 | 7/350 |
|
|
|
[考研] 085602 307分 求調(diào)劑 +7 | 不知道叫什么！ 2026-03-26 | 7/350 |
|
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫(kù) 2026-03-26 | 5/250 |
|
|
|
[考研] 351求調(diào)劑 +4 | 麥克阿磊 2026-03-24 | 4/200 |
|
|
|
[考研]
材料調(diào)劑
5+4
|
想要一壺桃花水 2026-03-25 | 10/500 |
|
|
|
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4 | 鋼鐵大炮 2026-03-24 | 4/200 |
|
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|