生物信息學(xué)從零開始學(xué)(5)——蛋白質(zhì)-小分子分子對接詳細(xì)步驟(docking)

疊合兩個蛋白質(zhì)三維結(jié)構(gòu)

蛋白質(zhì)三級結(jié)構(gòu)的比對

蛋白質(zhì)三維結(jié)構(gòu)比對就是對蛋白質(zhì)三維空間結(jié)構(gòu)的相似性進(jìn)行比較样屠。

  • 可用于探索蛋白質(zhì)進(jìn)化及同源關(guān)系
  • 改進(jìn)序列比對的精度
  • 改進(jìn)蛋白質(zhì)結(jié)構(gòu)預(yù)測工具
  • 為蛋白質(zhì)結(jié)構(gòu)分類提供依據(jù)
  • 幫助了解蛋白質(zhì)功能

結(jié)構(gòu)比對的結(jié)果可以用很多種參數(shù)來衡量搔驼,最常用的是root mean squared deviations(RMSD)迄靠。如果兩個結(jié)構(gòu)的RMSD為0埃瞳筏,則結(jié)構(gòu)一致,可以完全重合慌盯;一般來說RMSD小于3埃時周霉,認(rèn)為兩個結(jié)構(gòu)相似。

superpose蛋白結(jié)構(gòu)疊合
http://superpose.wishartlab.com/
SPDVB疊合
https://spdbv.vital-it.ch/

載入需要疊合的兩個蛋白質(zhì)PDB文件亚皂。

自動疊合

也可以選擇性疊合

## 蛋白質(zhì)表面的分子性質(zhì)
  • 表面形狀
  • 表面電荷分布
  • 表面殘基可溶性

首先需要下載VMD插件-APBS俱箱,下載后放入VMD安裝目錄。

  • 讀取需要分析的蛋白質(zhì)PDB文件jinrizhisen.PDB
  • 為jinrizhisen.PDB創(chuàng)建psf文件孕讳,psf文件是另一種存儲蛋白質(zhì)結(jié)構(gòu)信息的文件匠楚,

可以看到生成了以下文件

到這里,準(zhǔn)備工作結(jié)束厂财,下面開始做電荷分布

此時需要關(guān)閉VMD再重新打開芋簿,依次載入剛創(chuàng)建的jinrizhisen_autopsf.pdb文件和jinrizhisen_autopsf.psf文件。

run結(jié)束后產(chǎn)生了一個文件夾

刪除原來的分子璃饱,需要重新依次載入四個文件夾与斤,分別為jinrizhisen_autopsf.pdb,jinrizhisen_autopsf.psf荚恶,jinrizhisen_autopsf.pqr撩穿,pot.dx

上圖中藍(lán)色代表正電荷,紅色代表負(fù)電荷谒撼。白色不帶電荷食寡。

換個姿勢

同樣的方法也可以做殘基可溶性分布圖

誰能猜到這個蛋白是哪一類型的蛋白呢?或者具有哪些特別的結(jié)構(gòu)和功能廓潜?

蛋白質(zhì)四級結(jié)構(gòu)的獲取

  • 數(shù)據(jù)庫獲取
https://dip.doe-mbi.ucla.edu/dip/Main.cgi
https://thebiogrid.org/
https://string-db.org/
  • 分子對接(docking)兩種蛋白質(zhì)-蛋白質(zhì)分子對接
  • Rigid Docking
  • Flexible Docking

ZDOCK:全能在線對接工具

https://zdock.umassmed.edu/

半天打不開…… 分析一下

https://www.ebi.ac.uk/msd-srv/prot_int/cgi-bin/piserver

由于沒有上一步結(jié)果抵皱,這里也就先不分析了善榛。蛋白質(zhì)-小分子對接

#需要安裝python2.5環(huán)境,所有安裝路徑不能出現(xiàn)中文
#然后安裝autodock
http://autodock.scripps.edu/
#再安裝myltools
http://mgltools.scripps.edu/downloads

所有軟件安裝完畢后打開autodocktools

首先加工處理小分子ad_ligand.pdb

給小分子加H原子

給小分子加電荷

加工完成后指定為要對接的ligand

自動確定扭矩中心

確定對接過程中鍵的旋轉(zhuǎn)與否呻畸,綠色為可旋轉(zhuǎn)移盆。

加工完成,將加工結(jié)果保存為pdbqt文件伤为。

接下來處理蛋白質(zhì)分子 先把剛才載入的小分子刪除

載入蛋白質(zhì)分子

然后和處理小分子的過程一樣咒循,加H,加電荷

蛋白質(zhì)處理完成后绞愚,將其指定為要對接的蛋白質(zhì)

指定完成后叙甸,自動保存為pdbqt文件 載入要涉及到的原子

指定盒子的范圍

設(shè)定完成后保存

保存目前的設(shè)置

下面運行autogrid,為蛋白質(zhì)分子創(chuàng)建map

autogrid4?-p?test.gpf?-l?test.glg

autogrid4:?Successful?Completion.

運行完成后產(chǎn)生了一系列map文件

下面設(shè)置docking的參數(shù)文件,選定蛋白質(zhì)和小分子

選擇docking算法

輸出docking參數(shù)文件

開始docking

E:\Autodocktest>autodock4?-p?test.dpf?-l?test.dlg

________________________________________________________________________________

autodock4:?Successful?Completion?on?"DESKTOP-H3AAFU3"

________________________________________________________________________________

這就是最終的輸出結(jié)果文件

下面將結(jié)果顯示出來爽醋,首先同樣刪除之前載入的所有文件

再打開蛋白質(zhì)文件

現(xiàn)在小分子的位置是docking之前的位置

看看docking之后的位置

看一下前50個對接結(jié)果

下面看一下排名第一的對接結(jié)果的具體情況

可以看到小分子與蛋白質(zhì)上的哪些氨基酸產(chǎn)生作用蚁署,并列出了產(chǎn)生氫鍵的位置 保存排名第一的對接狀態(tài) 刪除除了小分子之外的全部分子

然后保存

看一下rank1.pdb

將其內(nèi)容拷貝至ad.protein.pdb,保存為result.pdb
此時result.pdb就是既包含蛋白質(zhì)也包含小分子的復(fù)合體結(jié)構(gòu)然后就可以用VMD打開result.pdb進(jìn)行可視化

分子對接(docking):虛擬篩選Virtual screening(VS)

虛擬篩選(virtual screening,VS)也稱計算機篩選蚂四,即在進(jìn)行生物活性篩選之前光戈,利用計算機上的分子對接軟件模擬目標(biāo)靶點與候選藥物之間的相互作用,計算兩者之間的親和力大小遂赠,以降低實際篩選化合物數(shù)目久妆,同時提高先導(dǎo)化合物發(fā)現(xiàn)效率□文溃化合物小分子庫ZINC

https://zinc.docking.org/
http://vina.scripps.edu/tutorial.html

分子對接(docking):反向?qū)?/h2>

靶標(biāo)數(shù)據(jù)庫

http://bioinfo-pharma.u-strasbg.fr/scPDB/

分子動力學(xué)模擬(molecular dynamic simulation)

涉及軟件:NAMD筷弦,CHARMM,DESMOND抑诸,GAUSS等

https://www.ks.uiuc.edu/Gallery/Movies/ChannelProteins/

Water Channels in Cell Membranes

How Bacteria Makes a Tail

Electrical Current through a Bacterial Toxin

注:整理自山東大學(xué)《生物信息學(xué)》

本文使用 文章同步助手 同步

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末烂琴,一起剝皮案震驚了整個濱河市,隨后出現(xiàn)的幾起案子蜕乡,更是在濱河造成了極大的恐慌奸绷,老刑警劉巖,帶你破解...
    沈念sama閱讀 206,013評論 6 481
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件层玲,死亡現(xiàn)場離奇詭異号醉,居然都是意外死亡,警方通過查閱死者的電腦和手機辛块,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,205評論 2 382
  • 文/潘曉璐 我一進(jìn)店門畔派,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人润绵,你說我怎么就攤上這事线椰。” “怎么了尘盼?”我有些...
    開封第一講書人閱讀 152,370評論 0 342
  • 文/不壞的土叔 我叫張陵憨愉,是天一觀的道長。 經(jīng)常有香客問我莱衩,道長,這世上最難降的妖魔是什么笨蚁? 我笑而不...
    開封第一講書人閱讀 55,168評論 1 278
  • 正文 為了忘掉前任趟庄,我火速辦了婚禮括细,結(jié)果婚禮上,老公的妹妹穿的比我還像新娘戚啥。我一直安慰自己奋单,他們只是感情好,可當(dāng)我...
    茶點故事閱讀 64,153評論 5 371
  • 文/花漫 我一把揭開白布猫十。 她就那樣靜靜地躺著览濒,像睡著了一般。 火紅的嫁衣襯著肌膚如雪拖云。 梳的紋絲不亂的頭發(fā)上贷笛,一...
    開封第一講書人閱讀 48,954評論 1 283
  • 那天,我揣著相機與錄音宙项,去河邊找鬼乏苦。 笑死,一個胖子當(dāng)著我的面吹牛尤筐,可吹牛的內(nèi)容都是我干的汇荐。 我是一名探鬼主播,決...
    沈念sama閱讀 38,271評論 3 399
  • 文/蒼蘭香墨 我猛地睜開眼盆繁,長吁一口氣:“原來是場噩夢啊……” “哼掀淘!你這毒婦竟也來了?” 一聲冷哼從身側(cè)響起改基,我...
    開封第一講書人閱讀 36,916評論 0 259
  • 序言:老撾萬榮一對情侶失蹤繁疤,失蹤者是張志新(化名)和其女友劉穎,沒想到半個月后秕狰,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體稠腊,經(jīng)...
    沈念sama閱讀 43,382評論 1 300
  • 正文 獨居荒郊野嶺守林人離奇死亡,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 35,877評論 2 323
  • 正文 我和宋清朗相戀三年鸣哀,在試婚紗的時候發(fā)現(xiàn)自己被綠了架忌。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點故事閱讀 37,989評論 1 333
  • 序言:一個原本活蹦亂跳的男人離奇死亡我衬,死狀恐怖叹放,靈堂內(nèi)的尸體忽然破棺而出饰恕,到底是詐尸還是另有隱情,我是刑警寧澤井仰,帶...
    沈念sama閱讀 33,624評論 4 322
  • 正文 年R本政府宣布埋嵌,位于F島的核電站,受9級特大地震影響俱恶,放射性物質(zhì)發(fā)生泄漏雹嗦。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點故事閱讀 39,209評論 3 307
  • 文/蒙蒙 一合是、第九天 我趴在偏房一處隱蔽的房頂上張望了罪。 院中可真熱鬧,春花似錦聪全、人聲如沸泊藕。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,199評論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽娃圆。三九已至,卻和暖如春鹤竭,著一層夾襖步出監(jiān)牢的瞬間踊餐,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 31,418評論 1 260
  • 我被黑心中介騙來泰國打工臀稚, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留吝岭,地道東北人。 一個月前我還...
    沈念sama閱讀 45,401評論 2 352
  • 正文 我出身青樓吧寺,卻偏偏與公主長得像窜管,于是被迫代替她去往敵國和親。 傳聞我的和親對象是個殘疾皇子稚机,可洞房花燭夜當(dāng)晚...
    茶點故事閱讀 42,700評論 2 345

推薦閱讀更多精彩內(nèi)容