 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2496 | 回復(fù): 14 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
aizazadi525鐵桿木蟲(chóng) (正式寫手)
|
[求助]
能量?jī)?yōu)化之后還是大于1000 已有2人參與
|
||
|
小蟲(chóng)剛接觸gromacs，現(xiàn)在在對(duì)一個(gè)自己對(duì)接的蛋白復(fù)合物進(jìn)行動(dòng)力學(xué)模擬，在能量?jī)?yōu)化過(guò)程中，使用如下.mdp文件 ; Title of run ; Parameters describing what to do, when to stop and what to save integrator = steep ; Algorithm (steep = steepest descent minimization) emtol = 1000.0 ; Stop minimization when the maximum force < 10.0 kJ/mol emstep = 0.01 ; Energy step size nsteps = 50000 ; Maximum number of (minimization) steps to perform energygrps = system ; Which energy group(s) to write to disk ; Parameters describing how to find the neighbors of each atom and how to calculate the interactions nstlist = 1 ; Frequency to update the neighbor list and long range forces ns_type = grid ; Method to determine neighbor list (simple, grid) rlist = 1.0 ; Cut-off for making neighbor list (short range forces) coulombtype = PME ; Treatment of long range electrostatic interactions rcoulomb = 1.0 ; long range electrostatic cut-off rvdw = 1.0 ; long range Van der Waals cut-off pbc = xyz ; Periodic Boundary Conditions (yes/no) 但是運(yùn)行總是很快就結(jié)束了，返回結(jié)果說(shuō)： Stepsize too small, or no change in energy. Converged to machine precision, but not to the requested precision Fmax< 1000 試著調(diào)了幾個(gè)參數(shù)，如把nsteps調(diào)整為100000，或把emstep調(diào)整為0.02或0.03，但還是出現(xiàn)同樣的結(jié)果。想著不報(bào)錯(cuò)就繼續(xù)運(yùn)行，但是在NVT平衡的時(shí)候，就會(huì)出現(xiàn)LINCE WARNING,出現(xiàn)太多程序就自動(dòng)停止了。 有沒(méi)有哪位大神知道一下，這種情況該如何處理。謝謝了�。� 2015-11-06_172833.png 2015-11-06_172803.png 2015-11-06_172600.png |
回復(fù)此樓» 猜你喜歡
招收博士1-2人
已經(jīng)有4人回復(fù)
-
08工科 320總分 求調(diào)劑
已經(jīng)有7人回復(fù)
-
生物學(xué)一志愿985，分?jǐn)?shù)349求調(diào)劑
已經(jīng)有7人回復(fù)
-
0856材料專碩353求調(diào)劑
已經(jīng)有4人回復(fù)
-
山東省面上項(xiàng)目限額評(píng)審
已經(jīng)有4人回復(fù)
-
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑
已經(jīng)有5人回復(fù)
-
生物學(xué)071000 329分求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿華中科技大學(xué)071000，求調(diào)劑
已經(jīng)有4人回復(fù)
-
生物學(xué)調(diào)劑
已經(jīng)有4人回復(fù)
-
求調(diào)劑院校信息
已經(jīng)有4人回復(fù)
aizazadi525
鐵桿木蟲(chóng) (正式寫手)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 5372.7
- 散金: 5
- 紅花: 7
- 帖子: 625
- 在線: 377.2小時(shí)
- 蟲(chóng)號(hào): 1146537
- 注冊(cè): 2010-11-14
- 性別: GG
- 專業(yè): 高分子合成化學(xué)
|
你好！我的蛋白是直接從PDB數(shù)據(jù)庫(kù)中下載的，使用Chimera軟件修復(fù)了氨基酸殘基，里邊的配體是我們自己分離得到的。我用sybyl軟件構(gòu)建了配體的三維結(jié)構(gòu)、并進(jìn)行了能量?jī)?yōu)化、構(gòu)象優(yōu)化。隨后使用了Surflex dock模塊進(jìn)行了對(duì)接，在對(duì)接前對(duì)蛋白進(jìn)行了加氫、加電荷及一些優(yōu)化。對(duì)對(duì)接結(jié)果另存為PDB文件，然后進(jìn)行MD模擬。這就是我的試驗(yàn)流程。我自己在網(wǎng)上也百度了一些解決方法，如2樓說(shuō)的修改原子的坐標(biāo)、修改integrator參數(shù)為cg、修改emstep為0.00001等方法，輸出結(jié)果還是stepsize too small。你說(shuō)的這個(gè)優(yōu)化蛋白，可不可以給以個(gè)具體的方法？謝謝了！ |
9樓2015-11-10 13:33:48
angjzy
木蟲(chóng) (小有名氣)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 1745.2
- 散金: 1491
- 紅花: 5
- 帖子: 251
- 在線: 172.3小時(shí)
- 蟲(chóng)號(hào): 647898
- 注冊(cè): 2008-11-06
- 專業(yè): 藥物材料
2樓2015-11-06 21:11:52
awaken2013
銅蟲(chóng) (小有名氣)
- 模擬EPI: 2
- 應(yīng)助: 26 (小學(xué)生)
- 金幣: 137.9
- 紅花: 15
- 帖子: 243
- 在線: 149.6小時(shí)
- 蟲(chóng)號(hào): 2854550
- 注冊(cè): 2013-12-07
- 性別: GG
- 專業(yè): 生物大分子結(jié)構(gòu)與功能
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
aizazadi525: 金幣+10, ★★★很有幫助 2015-11-09 11:07:28
感謝參與，應(yīng)助指數(shù) +1
aizazadi525: 金幣+10, ★★★很有幫助 2015-11-09 11:07:28
|
只能說(shuō)你自己對(duì)接的這個(gè)構(gòu)象，無(wú)法被能量最小化到1000以內(nèi)。（報(bào)錯(cuò)寫的很清楚） 樓上說(shuō)的方法可能可行。 若不行，還是仔細(xì)要看看你對(duì)接的結(jié)果 |
» 本帖已獲得的紅花（最新10朵）
aizazadi5253樓2015-11-09 02:10:58
aizazadi525
鐵桿木蟲(chóng) (正式寫手)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 5372.7
- 散金: 5
- 紅花: 7
- 帖子: 625
- 在線: 377.2小時(shí)
- 蟲(chóng)號(hào): 1146537
- 注冊(cè): 2010-11-14
- 性別: GG
- 專業(yè): 高分子合成化學(xué)
4樓2015-11-09 15:32:34
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 286求調(diào)劑 +10 | Faune 2026-03-21 | 10/500 |
|
|---|---|---|---|---|
|
|
[考研] 280求調(diào)劑 +11 | 咕嚕曉曉 2026-03-18 | 12/600 |
|
|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
[考研] 初試 317 +7 | 半拉月丙 2026-03-20 | 7/350 |
|
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +5 | 張vvvv 2026-03-15 | 7/350 |
|
|
|
[考研] 一志愿南大，0703化學(xué)，分?jǐn)?shù)336，求調(diào)劑 +3 | 收到VS 2026-03-21 | 3/150 |
|
|
|
[考研] 297求調(diào)劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
|
[考研] 材料 271求調(diào)劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
|
[考研] 22 350 本科985求調(diào)劑，求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|
|
[考研] 材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué) +4 | 順順順mr 2026-03-18 | 5/250 |
|
|
|
[考研] 機(jī)械專碩299求調(diào)劑至材料 +3 | kkcoco25 2026-03-16 | 4/200 |
|
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +12 | yangfz 2026-03-17 | 12/600 |
|
|
|
[考研] 307求調(diào)劑 +10 | 冷笙123 2026-03-17 | 10/500 |
|
|
|
[考研] 材料專業(yè)求調(diào)劑 +6 | hanamiko 2026-03-18 | 6/300 |
|
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
|
[考研] 一志愿福大288有機(jī)化學(xué)，求調(diào)劑 +3 | 小木蟲(chóng)200408204 2026-03-18 | 3/150 |
|
|
|
[考研] 0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專利，省、國(guó)級(jí)證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|