| 9 | 1/1 | 返回列表 |
| 查看: 2189 | 回復(fù): 8 | ||||
[交流]
Autodock對(duì)接結(jié)果可信度討論
|
||||
|
因?yàn)樽罱性谑褂胊utodock,所以有幾個(gè)問(wèn)題想請(qǐng)教和大家交流請(qǐng)教以下。 我是用Autodock做環(huán)糊精和藥物的自動(dòng)對(duì)接,本來(lái)是想先用Autodock找出藥物和環(huán)糊精最穩(wěn)定的結(jié)合位置和構(gòu)象,然后把該構(gòu)象拿到Gaussian進(jìn)一步優(yōu)化。 一開(kāi)始,我把格子大小設(shè)為60*60*60,可是結(jié)果出來(lái)的重復(fù)性不是很好,最高的那一簇能量總是變化,當(dāng)然每簇之間相差倒也不大。 后來(lái),我干脆把格子大小設(shè)為40*40*40,然后又提高對(duì)接次數(shù),結(jié)果出來(lái)后最高的那一簇能量到是符合了,而且也是能量最低點(diǎn)?墒牵疫不滿意的是,兩次對(duì)接結(jié)果中,一簇高為81,另一簇高為53,這個(gè)偏差就有點(diǎn)大了吧。 我想請(qǐng)教的是: 1。其實(shí)我把格子設(shè)為40*40*40只是把格子大小為60*60*60的對(duì)接結(jié)果(能量最低點(diǎn))拿來(lái)然后調(diào)整格子大小至40*40*40,觀察并沒(méi)有基團(tuán)伸到格子外,就拿來(lái)用了。但我總覺(jué)得這樣做好象不是很科學(xué)。請(qǐng)問(wèn),有沒(méi)有什么方法可以觀察這個(gè)對(duì)接設(shè)置是否可行? 2。我兩次對(duì)接結(jié)果能量最低點(diǎn)確實(shí)是最高簇,可是簇的高度相差有點(diǎn)大,一個(gè)81,一個(gè)53。這樣的話,結(jié)果是否可信呢?如果不可信,又該怎么辦呢? 3。在對(duì)接結(jié)果中,總會(huì)提到RMSD,也就是均方根偏差,就是結(jié)果構(gòu)象和目標(biāo)構(gòu)象的的偏差統(tǒng)計(jì)。我的理解就是每簇中的所有對(duì)接構(gòu)像并不都是一樣的,只是剛好RMSD在這個(gè)范圍內(nèi),所以才放到同一簇。所以才會(huì)有“結(jié)果中選簇最高中的能量最低點(diǎn)這個(gè)說(shuō)法”?我的理解就是即使有能量更低的簇,可是分子在那個(gè)構(gòu)像附近出現(xiàn)的概率較低。不知道這個(gè)看法對(duì)不對(duì)。 |
分子模擬 |
» 搶金幣啦!回帖就可以得到:
+1/183
+2/116
+1/90
+1/72
+2/40
+1/38
+1/36
+1/14
+1/10
+1/8
+1/6
+1/5
+1/5
+1/4
+1/3
+1/3
+1/3
+1/3
+1/3
+1/3
木蟲(chóng)之王 (文學(xué)泰斗)
銀蟲(chóng) (正式寫(xiě)手)
|
1.關(guān)于盒子大小的設(shè)置,用Autodock作對(duì)接你首先得充分了解你的活性部位,然后把它們?nèi)及鼑涂梢粤;虼蠡蛐?yīng)該問(wèn)題不大。因?yàn)楸緛?lái)這個(gè)對(duì)接就是個(gè)參考,計(jì)算過(guò)程中的參數(shù)好多都是近似的。并且算是半柔性的。 2.如果你想更可靠,你可以作二次對(duì)接,甚至三次對(duì)接,不過(guò)花費(fèi)的時(shí)間多。 3.RMSD是可以自己設(shè)置的,自己設(shè)置的越小分簇可能越多,如果最低能量和最高簇是同一簇,那么結(jié)果就會(huì)更可靠。 玩兒了近三年了。自己理解的吧。具體看你看文獻(xiàn)自己理解了。 |
銀蟲(chóng) (正式寫(xiě)手)


| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 289求調(diào)劑 +8 | 新時(shí)代材料 2026-03-27 | 8/400 |
|
|
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4 | Ripcord上岸 2026-03-27 | 4/200 |
|
|
[考研] 292求調(diào)劑 +13 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 14/700 |
|
|
[考研] 339求調(diào)劑 +4 | 烤麥芽 2026-03-27 | 5/250 |
|
|
[考研] 一志愿北京化工大學(xué)材料與化工 264分各科過(guò)A區(qū)國(guó)家線 +10 | 哈哈157349 2026-03-21 | 10/500 |
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
[考研] 303求調(diào)劑 +7 | 安憶靈 2026-03-22 | 8/400 |
|
|
[論文投稿] Journal of Mechanical Science and Technology +3 | Russ_ss 2026-03-25 | 5/250 |
|
|
[考研] 317求調(diào)劑 +7 | 蛋黃咸肉粽 2026-03-26 | 7/350 |
|
|
[考研] 071000生物學(xué)求調(diào)劑,初試成績(jī)343 +6 | 小小甜面團(tuán) 2026-03-25 | 6/300 |
|
|
[考研] 321求調(diào)劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|
|
[考研] 086000生物與醫(yī)藥292求調(diào)劑 +6 | 小小陳小小 2026-03-22 | 9/450 |
|
|
[考研] 資源與環(huán)境 調(diào)劑申請(qǐng)(333分) +9 | holy J 2026-03-21 | 9/450 |
|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap- 2026-03-22 | 10/500 |
|
|
[考研] 一志愿吉林大學(xué)材料與化工303分求調(diào)劑 +4 | 為學(xué)666 2026-03-24 | 4/200 |
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +10 | 紙魚(yú)ly 2026-03-21 | 11/550 |
|
|
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3 | 開(kāi)開(kāi)心心沒(méi)煩惱 2026-03-23 | 4/200 |
|
|
[考研] 材料/農(nóng)業(yè)專業(yè),07/08開(kāi)頭均可,過(guò)線就行 +3 | 呵唔哦豁 2026-03-23 | 4/200 |
|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|