第十章分子對(duì)接

1、第十章第十章 藥物設(shè)計(jì)的基本方法藥物設(shè)計(jì)的基本方法分子對(duì)接分子對(duì)接內(nèi)容內(nèi)容分子對(duì)接簡(jiǎn)介分子對(duì)接簡(jiǎn)介分子對(duì)接法的基本原理分子對(duì)接法的基本原理分子對(duì)接軟件分子對(duì)接軟件DockAutoDockCADD基于結(jié)構(gòu)方法分類10.1 引言引言 分子對(duì)接就是兩個(gè)或多個(gè)分子之間通過分子對(duì)接就是兩個(gè)或多個(gè)分子之間通過幾何匹配幾何匹配和和能量匹配能量匹配而相互識(shí)別的過程。而相互識(shí)別的過程。 分子對(duì)接在酶學(xué)研究以及藥物設(shè)計(jì)中具有分子對(duì)接在酶學(xué)研究以及藥物設(shè)計(jì)中具有十分重要的意義。十分重要的意義。 在酶激活劑、酶抑制劑與酶相互作用以及藥物分在酶激活劑、酶抑制劑與酶相互作用以及藥物分子產(chǎn)生藥理反應(yīng)的過程中，小分子（通常意

2、義上的子產(chǎn)生藥理反應(yīng)的過程中，小分子（通常意義上的LigandLigand）與靶酶（通常意義上的）與靶酶（通常意義上的ReceptorReceptor）相互互結(jié)合，）相互互結(jié)合，首先就需要兩個(gè)分子充分接近，采取合適的取向，使首先就需要兩個(gè)分子充分接近，采取合適的取向，使兩者在必要的部位相互契合，發(fā)生相互作用，繼而通兩者在必要的部位相互契合，發(fā)生相互作用，繼而通過適當(dāng)?shù)臉?gòu)象調(diào)整，得到一個(gè)穩(wěn)定的復(fù)合物構(gòu)象。過適當(dāng)?shù)臉?gòu)象調(diào)整，得到一個(gè)穩(wěn)定的復(fù)合物構(gòu)象。 通過分子對(duì)接確定復(fù)合物中兩個(gè)分子正確的相對(duì)通過分子對(duì)接確定復(fù)合物中兩個(gè)分子正確的相對(duì)位置和取向，研究?jī)蓚€(gè)分子的構(gòu)象，特別是底物構(gòu)象位置和取向，研究?jī)?/p>

3、個(gè)分子的構(gòu)象，特別是底物構(gòu)象在形成復(fù)合物過程中的變化，是確定酶激活劑、抑制在形成復(fù)合物過程中的變化，是確定酶激活劑、抑制劑作用機(jī)制以及藥物作用機(jī)制，設(shè)計(jì)新藥的基礎(chǔ)。劑作用機(jī)制以及藥物作用機(jī)制，設(shè)計(jì)新藥的基礎(chǔ)。 大腸桿菌乙酰輔酶A羧化酶肉桂醛 在藥物設(shè)計(jì)中，分子對(duì)接方法主要用來從小分子數(shù)據(jù)庫中在藥物設(shè)計(jì)中，分子對(duì)接方法主要用來從小分子數(shù)據(jù)庫中搜尋與受體生物大分子有較好親和力的小分子，進(jìn)行藥理搜尋與受體生物大分子有較好親和力的小分子，進(jìn)行藥理測(cè)試，從中發(fā)現(xiàn)新的先導(dǎo)化合物。測(cè)試，從中發(fā)現(xiàn)新的先導(dǎo)化合物。 優(yōu)點(diǎn)：小分子配體來自現(xiàn)有化合物數(shù)據(jù)庫或制藥公司自己優(yōu)點(diǎn)：小分子配體來自現(xiàn)有化合物數(shù)據(jù)庫或制藥公司

4、自己 的數(shù)據(jù)庫，合成路線比較成熟；的數(shù)據(jù)庫，合成路線比較成熟； 分子對(duì)接方法從整體上考慮受體與配體結(jié)合的效果分子對(duì)接方法從整體上考慮受體與配體結(jié)合的效果 加速先導(dǎo)化合物的發(fā)現(xiàn)進(jìn)程。加速先導(dǎo)化合物的發(fā)現(xiàn)進(jìn)程。10.2 分子對(duì)接法的基本原理分子對(duì)接法的基本原理 配體與受體的結(jié)合過程是一個(gè)很復(fù)雜的過程，涉及到配體和配體與受體的結(jié)合過程是一個(gè)很復(fù)雜的過程，涉及到配體和受體的去溶劑化、配體和受體（主要是活性位點(diǎn)處的殘基）受體的去溶劑化、配體和受體（主要是活性位點(diǎn)處的殘基）的構(gòu)象變化以及配體與受體之間的相互作用。的構(gòu)象變化以及配體與受體之間的相互作用。 配體與受體的結(jié)合強(qiáng)弱取決于結(jié)合過程的自由能變化：配體

5、與受體的結(jié)合強(qiáng)弱取決于結(jié)合過程的自由能變化：lnbindingdGRTK 結(jié)合自由能也可以表示成：結(jié)合自由能也可以表示成：bindingbindingbindingGHT S 大部分的分子對(duì)接法忽略了全部的熵效應(yīng)，而在焓效應(yīng)中也只考慮配體與受體的相互作用能 配體與受體的相互作用：配體與受體的相互作用：bondHticelectrostavdweractionEEEEint其中氫鍵采用的形式和力場(chǎng)有關(guān)配體與受體的相互作用：配體與受體的相互作用：分子對(duì)接計(jì)算是把配體分子放在受體活性位點(diǎn)分子對(duì)接計(jì)算是把配體分子放在受體活性位點(diǎn)的位置，然后按照幾何互補(bǔ)、能量互補(bǔ)、化學(xué)的位置，然后按照幾何互補(bǔ)、能量互補(bǔ)

6、、化學(xué)環(huán)境互補(bǔ)的原則來實(shí)時(shí)評(píng)價(jià)配體與受體相互作環(huán)境互補(bǔ)的原則來實(shí)時(shí)評(píng)價(jià)配體與受體相互作用的好壞，并找到兩個(gè)分子之間最佳的結(jié)合模用的好壞，并找到兩個(gè)分子之間最佳的結(jié)合模式式 。 分子對(duì)接最初思想起源于分子對(duì)接最初思想起源于Fisher E.的的“鎖和鑰鎖和鑰匙模型匙模型”，認(rèn)為，認(rèn)為“鎖鎖”和和“鑰匙鑰匙”的相識(shí)別的的相識(shí)別的首要條件是他們?cè)诳臻g形狀上要互相匹配首要條件是他們?cè)诳臻g形狀上要互相匹配 配體和受體分子之間的識(shí)別要比配體和受體分子之間的識(shí)別要比“鎖和鑰匙鎖和鑰匙”模型復(fù)雜的多：模型復(fù)雜的多： 1.1.配體和受體分子的構(gòu)象是變化的，而不是剛性的，配體和配體和受體分子的構(gòu)象是變化的，而不是

7、剛性的，配體和受體在對(duì)接過程中互相適應(yīng)對(duì)方，從而達(dá)到更完美的匹配；受體在對(duì)接過程中互相適應(yīng)對(duì)方，從而達(dá)到更完美的匹配；2.2.其次，分子對(duì)接不但要滿足空間形狀的匹配，還要滿足能其次，分子對(duì)接不但要滿足空間形狀的匹配，還要滿足能量的匹配。量的匹配。 配體和受體之間的通過底物分子與靶酶分子能否結(jié)合以及結(jié)合的強(qiáng)度最終是由形成此復(fù)合物過程的結(jié)合自由能變化Gbind所決定的。 互補(bǔ)性（互補(bǔ)性（complementarity）和預(yù)組織）和預(yù)組織（pre-organization）是決定分子對(duì)接過）是決定分子對(duì)接過程的兩個(gè)重要原則，前者決定識(shí)別過程的程的兩個(gè)重要原則，前者決定識(shí)別過程的選擇性，而后者決定識(shí)別

8、過理的結(jié)合能力。選擇性，而后者決定識(shí)別過理的結(jié)合能力。互補(bǔ)性包括空間結(jié)構(gòu)的互補(bǔ)性和電學(xué)性質(zhì)互補(bǔ)性包括空間結(jié)構(gòu)的互補(bǔ)性和電學(xué)性質(zhì)的互補(bǔ)性。的互補(bǔ)性。 1958年年Koshland提出了分子識(shí)別過程中提出了分子識(shí)別過程中的誘導(dǎo)契合（的誘導(dǎo)契合（induced fit）概念，指出配）概念，指出配體與受體相互結(jié)合時(shí)，受體將采取一個(gè)能體與受體相互結(jié)合時(shí)，受體將采取一個(gè)能同底物達(dá)到最佳結(jié)合的構(gòu)象。同底物達(dá)到最佳結(jié)合的構(gòu)象。分子對(duì)接方法根據(jù)不同的簡(jiǎn)化程度可以大致分分子對(duì)接方法根據(jù)不同的簡(jiǎn)化程度可以大致分為以下三類：為以下三類： （1）剛性對(duì)接；研究體系的構(gòu)象不發(fā)生變化）剛性對(duì)接；研究體系的構(gòu)象不發(fā)生變化 （

9、2）半柔性對(duì)接；研究體系尤其是配體的構(gòu)象）半柔性對(duì)接；研究體系尤其是配體的構(gòu)象允許在一定的范圍內(nèi)變化允許在一定的范圍內(nèi)變化 （3）柔性對(duì)接；研究體系的構(gòu)象基本上是可以）柔性對(duì)接；研究體系的構(gòu)象基本上是可以自由變化的自由變化的 剛性對(duì)接適合考察比較大的體系，比如蛋白質(zhì)和蛋白質(zhì)以剛性對(duì)接適合考察比較大的體系，比如蛋白質(zhì)和蛋白質(zhì)以及蛋白質(zhì)和核酸之間的相互作用，他計(jì)算較為簡(jiǎn)單，原理及蛋白質(zhì)和核酸之間的相互作用，他計(jì)算較為簡(jiǎn)單，原理也相對(duì)簡(jiǎn)單，主要是考慮構(gòu)象之間的契合程度。也相對(duì)簡(jiǎn)單，主要是考慮構(gòu)象之間的契合程度。 半柔性對(duì)接適合于處理小分子和大分子之間的對(duì)接。在對(duì)半柔性對(duì)接適合于處理小分子和大分子之間

10、的對(duì)接。在對(duì)接過程中，小分子的構(gòu)象一般是可以變化的，但大分子則接過程中，小分子的構(gòu)象一般是可以變化的，但大分子則是剛性的。由于小分子相對(duì)較小，因此在一定程度考察柔是剛性的。由于小分子相對(duì)較小，因此在一定程度考察柔性的基礎(chǔ)上，還可以保持較高的計(jì)算效率，在藥物設(shè)計(jì)尤性的基礎(chǔ)上，還可以保持較高的計(jì)算效率，在藥物設(shè)計(jì)尤其在基于分子數(shù)據(jù)庫的虛擬篩選過程中，一般采用半柔性其在基于分子數(shù)據(jù)庫的虛擬篩選過程中，一般采用半柔性的分子對(duì)接方法。的分子對(duì)接方法。 柔性對(duì)接方法一般用于精確考察分子之間的識(shí)別情況，由柔性對(duì)接方法一般用于精確考察分子之間的識(shí)別情況，由于在計(jì)算過程中體系的構(gòu)象是可以變化的，因此柔性對(duì)接于在

11、計(jì)算過程中體系的構(gòu)象是可以變化的，因此柔性對(duì)接在提高了對(duì)接準(zhǔn)確性的同時(shí)卻需要耗費(fèi)較長(zhǎng)的計(jì)算時(shí)間。在提高了對(duì)接準(zhǔn)確性的同時(shí)卻需要耗費(fèi)較長(zhǎng)的計(jì)算時(shí)間。分子對(duì)接的目的：找到底物分子和受體分子的分子對(duì)接的目的：找到底物分子和受體分子的最佳結(jié)合位置最佳結(jié)合位置找到最佳的結(jié)合位置評(píng)價(jià)對(duì)接分子之間的結(jié)合強(qiáng)度構(gòu)象搜索方法按照分子對(duì)接程序?qū)ふ遗潴w與受體結(jié)合構(gòu)按照分子對(duì)接程序?qū)ふ遗潴w與受體結(jié)合構(gòu)象方法的差異，可以分為三種：象方法的差異，可以分為三種： 1、局部?jī)?yōu)化法、局部?jī)?yōu)化法 2、深度搜尋法、深度搜尋法 3、全局優(yōu)化法、全局優(yōu)化法局部?jī)?yōu)化法：不對(duì)配體和受體進(jìn)行構(gòu)象搜尋，局部?jī)?yōu)化法：不對(duì)配體和受體進(jìn)行構(gòu)象搜尋，只

12、是對(duì)初始構(gòu)象進(jìn)行優(yōu)化，得到配體與受體結(jié)只是對(duì)初始構(gòu)象進(jìn)行優(yōu)化，得到配體與受體結(jié)合的一個(gè)局部最優(yōu)構(gòu)象；合的一個(gè)局部最優(yōu)構(gòu)象；深度搜尋法：采用深度優(yōu)先搜尋法或廣度優(yōu)先深度搜尋法：采用深度優(yōu)先搜尋法或廣度優(yōu)先搜尋法，通過有限的步驟，找到一個(gè)相對(duì)較好搜尋法，通過有限的步驟，找到一個(gè)相對(duì)較好的局部最優(yōu)結(jié)合構(gòu)象。的局部最優(yōu)結(jié)合構(gòu)象。全局優(yōu)化法：在進(jìn)行構(gòu)象搜尋時(shí)，利用模擬退全局優(yōu)化法：在進(jìn)行構(gòu)象搜尋時(shí)，利用模擬退火算法或遺傳算法，尋找配體與受體的全局最火算法或遺傳算法，尋找配體與受體的全局最優(yōu)結(jié)合構(gòu)象。優(yōu)結(jié)合構(gòu)象。對(duì)接類型對(duì)接類型對(duì)接方法類型對(duì)接方法類型對(duì)接方法對(duì)接方法Flexible-Ligand Doc

13、king柔性配體對(duì)接Systematic系統(tǒng)方法Conformational 構(gòu)象Fragmentation 片段生長(zhǎng)Database 數(shù)據(jù)庫Random/stochastic隨機(jī)方法Monte Carlo (MC) 蒙特卡羅Genetic algorithm (GA) 遺傳算法Tabu Search 禁忌搜索Simulation methods模擬方法Molecular dynamics (MD) 分子動(dòng)力學(xué)Energy minimization 能量最小化Flexible-Protein Docking柔性蛋白對(duì)接Molecular dynamics (MD) 分子模擬Monte Carl

14、o (MC) 蒙特卡羅Rotamer libraries 旋轉(zhuǎn)異構(gòu)體庫Protein-ensemble grids 蛋白集合柵格Soft-receptor modeling 軟受體建模表1： 常用的分子對(duì)接方法 10.3 分子對(duì)接軟件分子對(duì)接軟件表2：有代表性的分子對(duì)接軟件名稱名稱優(yōu)化方法優(yōu)化方法評(píng)價(jià)函數(shù)評(píng)價(jià)函數(shù)速度速度DOCK片斷生長(zhǎng)分子力場(chǎng)、表面匹配得分、化學(xué)環(huán)境匹配得分快AutoDock遺傳算法半經(jīng)驗(yàn)自由能評(píng)價(jià)函數(shù)一般ICM-Docking隨機(jī)全局優(yōu)化半經(jīng)驗(yàn)自由能評(píng)價(jià)函數(shù)快GOLD遺傳算法半經(jīng)驗(yàn)自由能評(píng)價(jià)函數(shù)快FlexX片斷生長(zhǎng)半經(jīng)驗(yàn)自由能評(píng)價(jià)函數(shù)快Affinity蒙特卡羅/分子力學(xué)/分

15、子動(dòng)力學(xué)分子力場(chǎng)慢ZDock&RDock幾何匹配/分子動(dòng)力學(xué)CAPRI* /分子力場(chǎng)慢FlexiDock遺傳算法分子力場(chǎng)慢eHiTS系統(tǒng)搜索半經(jīng)驗(yàn)自由能評(píng)價(jià)函數(shù)快Hex幾何匹配CAPRI*快分子對(duì)接軟件最常用的分子對(duì)接軟件被引次數(shù)統(tǒng)計(jì) 10.4 DOCK 基本情況基本情況 第一個(gè)分子對(duì)接程序。由加利福尼亞州立大學(xué)舊金山分校第一個(gè)分子對(duì)接程序。由加利福尼亞州立大學(xué)舊金山分校Kuntz小組于小組于1982年開發(fā)。目前版本為年開發(fā)。目前版本為DOCK6.610.4.1 基本原理基本原理10.4.1.1 形狀匹配形狀匹配受體與配體的原子之間的距離大于其范德受體與配體的原子之間的距離大于其范德華

16、半徑之和；華半徑之和；配體中能形成氫鍵的原子，都有一個(gè)受體配體中能形成氫鍵的原子，都有一個(gè)受體的氮原子或氧原子在其的氮原子或氧原子在其0.35 nm范圍之內(nèi)范圍之內(nèi)所有的配體原子都在受體的結(jié)合位點(diǎn)內(nèi)所有的配體原子都在受體的結(jié)合位點(diǎn)內(nèi)Dock利用利用球形集合球形集合來表示受體活性位點(diǎn)以來表示受體活性位點(diǎn)以及配體的形狀。及配體的形狀。對(duì)于受體，程序產(chǎn)生一個(gè)球集，使之填充受體分子表面的口袋。如果配體與受體匹配較好的話，配體的球集應(yīng)該能較好地填充受體的球集。10.4.1.2 匹配算法匹配算法 大分子配體的對(duì)接時(shí)間計(jì)算時(shí)間較長(zhǎng)；在大分子配體的對(duì)接時(shí)間計(jì)算時(shí)間較長(zhǎng)；在復(fù)合物結(jié)構(gòu)中，由于配體的大部分表面不復(fù)

17、合物結(jié)構(gòu)中，由于配體的大部分表面不會(huì)參與相互作用。因此，將表面分成不同會(huì)參與相互作用。因此，將表面分成不同的部分，每一部分分別進(jìn)行匹配，這一過的部分，每一部分分別進(jìn)行匹配，這一過程稱為分簇（程稱為分簇（cluster）.也就是說采取了也就是說采取了“分而治之分而治之”的匹配算法。的匹配算法。10.4.1.3 基于經(jīng)驗(yàn)勢(shì)能函數(shù)的評(píng)價(jià)函數(shù) Dock中引入經(jīng)驗(yàn)勢(shì)能函數(shù)做為評(píng)價(jià)函數(shù)。 優(yōu)點(diǎn)：預(yù)先對(duì)受體原子的作用進(jìn)行加和，并存入 網(wǎng)格文件中，計(jì)算能力增大，對(duì)接結(jié)果精確。 Dock采用分子力學(xué)方法，并基于Amber力場(chǎng)參數(shù)。力場(chǎng)中的非鍵作用只包含靜電作用和范德力場(chǎng)中的非鍵作用只包含靜電作用和范德華作用，氫

18、鍵作為靜電作用處理。華作用，氫鍵作為靜電作用處理。10.4.2 柔性對(duì)接柔性對(duì)接Dock4.0之前的版本沒有考慮配體的柔性。之前的版本沒有考慮配體的柔性。 之后的版本，增加了配體和受體的柔性，之后的版本，增加了配體和受體的柔性，大大改善對(duì)接結(jié)果。大大改善對(duì)接結(jié)果?？紤]配體分子的柔性確定剛性片段柔性搜尋1.確定剛性片段確定剛性片段 柔性分子可以看成由剛性片段組成；柔性分子可以看成由剛性片段組成； 剛性片段中的原子由不能旋轉(zhuǎn)的鍵連接；剛性片段中的原子由不能旋轉(zhuǎn)的鍵連接； 相鄰的剛性片段由可以旋轉(zhuǎn)的鍵連接。相鄰的剛性片段由可以旋轉(zhuǎn)的鍵連接。 環(huán)的處理確定可旋轉(zhuǎn)的鍵2. 柔性搜尋柔性搜尋 Dock有兩

19、種柔性搜尋方法： 一種是錨優(yōu)先搜尋（anchor-first search） 一種是同時(shí)搜尋（simultaneous search） a. 錨優(yōu)先搜尋錨優(yōu)先搜尋對(duì)接錨片段加一個(gè)剛性片段，先內(nèi)后外 先大后小3）構(gòu)象搜尋，得到NcNt個(gè)構(gòu)象4）根據(jù)構(gòu)象能量差異按照一定 標(biāo)準(zhǔn)，選取Nc個(gè)構(gòu)象5）回到第2步，不斷加錨，直到 所有的錨片段都對(duì)接過，跳過 5，得到對(duì)接構(gòu)象。b. 同時(shí)搜尋 先進(jìn)行構(gòu)象搜尋，然后將搜尋得到的構(gòu)先進(jìn)行構(gòu)象搜尋，然后將搜尋得到的構(gòu)象分別對(duì)接到受體的結(jié)合位點(diǎn)中。象分別對(duì)接到受體的結(jié)合位點(diǎn)中。 根據(jù)參數(shù)設(shè)置不同，可以進(jìn)行系統(tǒng)搜尋根據(jù)參數(shù)設(shè)置不同，可以進(jìn)行系統(tǒng)搜尋或隨機(jī)搜尋。或隨機(jī)搜

20、尋。10.5 AutoDock AutoDock是是The Scripps Research Institute的的OIson科研小組使用科研小組使用C語語言開發(fā)的分子對(duì)接軟件包，目前最言開發(fā)的分子對(duì)接軟件包，目前最新的版本為新的版本為4.2。 AutoDock是一個(gè)軟件包，其中主要是一個(gè)軟件包，其中主要包含包含AutoGrid和和AutoDock兩個(gè)程序。兩個(gè)程序。其中其中AutoGrid主要負(fù)責(zé)格點(diǎn)中相關(guān)主要負(fù)責(zé)格點(diǎn)中相關(guān)能量的計(jì)算，而能量的計(jì)算，而AutoDock則負(fù)責(zé)構(gòu)則負(fù)責(zé)構(gòu)象搜索及評(píng)價(jià)。象搜索及評(píng)價(jià)。 10.5.1 評(píng)價(jià)函數(shù)評(píng)價(jià)函數(shù)A. 3.0以前的版本采用經(jīng)驗(yàn)勢(shì)函數(shù)為評(píng)價(jià)函數(shù)以前

21、的版本采用經(jīng)驗(yàn)勢(shì)函數(shù)為評(píng)價(jià)函數(shù)B. 3.0版本開始采用經(jīng)驗(yàn)結(jié)合自由能為評(píng)價(jià)函數(shù)版本開始采用經(jīng)驗(yàn)結(jié)合自由能為評(píng)價(jià)函數(shù)A. 經(jīng)驗(yàn)勢(shì)能函數(shù)經(jīng)驗(yàn)勢(shì)能函數(shù) AutoDock的經(jīng)驗(yàn)勢(shì)能函數(shù)包括三項(xiàng)：靜電的經(jīng)驗(yàn)勢(shì)能函數(shù)包括三項(xiàng)：靜電作用、范德華作用和氫鍵作用。作用、范德華作用和氫鍵作用。 力場(chǎng)采用力場(chǎng)采用Amber力場(chǎng)力場(chǎng) 根據(jù)熱力學(xué)定律，離解過程的自由能變化根據(jù)熱力學(xué)定律，離解過程的自由能變化和焓的變化：和焓的變化： 實(shí)驗(yàn)測(cè)定的自由能與離解常數(shù)：實(shí)驗(yàn)測(cè)定的自由能與離解常數(shù)：GHT S 0lndGGRTK B. 經(jīng)驗(yàn)結(jié)合自由能經(jīng)驗(yàn)結(jié)合自由能 當(dāng)體系達(dá)到平衡時(shí)，自由能等于當(dāng)體系達(dá)到平衡時(shí)，自由能等于00lnd

22、GRTK配體配體-受體結(jié)合的熱力學(xué)受體結(jié)合的熱力學(xué)配體配體-受體結(jié)合的熱力學(xué)過程受體結(jié)合的熱力學(xué)過程 藥物與受體存在與體液環(huán)藥物與受體存在與體液環(huán)境中，與周圍的水分子存境中，與周圍的水分子存在溶劑化作用。在溶劑化作用。對(duì)接過程中的去溶劑化效應(yīng) Autodock的經(jīng)驗(yàn)結(jié)合自由能評(píng)價(jià)函數(shù)：的經(jīng)驗(yàn)結(jié)合自由能評(píng)價(jià)函數(shù)：,1261210,()( )()i ji ji jI Jvdwhbondi ji jijijijijACDBGGGE trrrr ,( )ijelectortorsoli jijijq qGG NGr r AutoDock在早期版本中使用的是模擬退火算法（Simulated Anneal

23、ing Algorithm）來尋找配體與受體最佳的結(jié)合位置狀態(tài)，而從3.0版本開始使用一種改良的遺傳算法，即拉馬克遺傳算法（Lamarckian Genetic Algorithm，LGA）。 測(cè)試結(jié)果表明，LGA比傳統(tǒng)的遺傳算法和模擬退火具有更高的效率。在LGA方法中，作者把遺傳算法和局部搜索（Local search）結(jié)合在一起，遺傳算法用于全局搜索，而局部搜索用于能量?jī)?yōu)化。 10.5.2 柔性搜尋方法柔性搜尋方法a.模擬退火算法模擬退火算法 指定結(jié)合位點(diǎn)隨機(jī)地差生配體 的一個(gè)構(gòu)象隨機(jī)改變配體的位置、取向， 產(chǎn)生了一個(gè)新的構(gòu)象，并比較 兩者與受體的結(jié)合能 如果新的構(gòu)象能低于原來的，則 接受

24、新構(gòu)象如果能量高于原來的，則降低溫 度，在新溫度下重復(fù)步驟24. 在同一初始條件下反復(fù)進(jìn)行多次b. 遺傳算法遺傳算法模仿生物的遺傳和進(jìn)化，對(duì)多極值優(yōu)化問模仿生物的遺傳和進(jìn)化，對(duì)多極值優(yōu)化問題的解群進(jìn)行選擇、雜交和突變操作，使題的解群進(jìn)行選擇、雜交和突變操作，使得解群不斷優(yōu)化，直至收斂于全局最優(yōu)解得解群不斷優(yōu)化，直至收斂于全局最優(yōu)解或近似全局最優(yōu)解?；蚪迫肿顑?yōu)解。 AutoDock優(yōu)點(diǎn)優(yōu)點(diǎn)1. 同時(shí)考慮配體的柔性和蛋白質(zhì)結(jié)合部位同時(shí)考慮配體的柔性和蛋白質(zhì)結(jié)合部位側(cè)鏈的柔性側(cè)鏈的柔性2. 網(wǎng)格計(jì)算，加快計(jì)算網(wǎng)格計(jì)算，加快計(jì)算網(wǎng)格計(jì)算和網(wǎng)格計(jì)算和dock的區(qū)別：的區(qū)別：Dock中，格點(diǎn)上保存的不

25、是能量，而是中，格點(diǎn)上保存的不是能量，而是受體上原子作用加和；受體上原子作用加和；Autodock中，格點(diǎn)上保存的是探針原子中，格點(diǎn)上保存的是探針原子和受體之間的相互作用能。和受體之間的相互作用能。AutoDock格點(diǎn)作用能計(jì)算格點(diǎn)作用能計(jì)算通過通過AutoDock程序包中的程序包中的Autogrid程序程序完成，格點(diǎn)計(jì)算包括三部分：范德華作用、完成，格點(diǎn)計(jì)算包括三部分：范德華作用、氫鍵作用和庫倫作用氫鍵作用和庫倫作用對(duì)于范德華相互作用的計(jì)算，對(duì)于范德華相互作用的計(jì)算， 每個(gè)格點(diǎn)上每個(gè)格點(diǎn)上保存的范德華能量的值的數(shù)目與要對(duì)接的保存的范德華能量的值的數(shù)目與要對(duì)接的配體上的原子類型的數(shù)目一樣。配體

26、上的原子類型的數(shù)目一樣。如果一個(gè)配件中含有如果一個(gè)配件中含有C、O和和H三種原子類三種原子類型，那么在每個(gè)格點(diǎn)上就需要用三個(gè)探針型，那么在每個(gè)格點(diǎn)上就需要用三個(gè)探針原子來計(jì)算探針原子與受體之間的范德華原子來計(jì)算探針原子與受體之間的范德華相互作用值。當(dāng)配體和受體進(jìn)行分子對(duì)接相互作用值。當(dāng)配體和受體進(jìn)行分子對(duì)接時(shí)，配體中某個(gè)原子和受體之間的相互作時(shí)，配體中某個(gè)原子和受體之間的相互作用能通過周圍用能通過周圍8個(gè)格點(diǎn)上的這種原子類型個(gè)格點(diǎn)上的這種原子類型為探針的格點(diǎn)值用內(nèi)插法得到。為探針的格點(diǎn)值用內(nèi)插法得到。 靜電相互作用的計(jì)算采用了一個(gè)靜電勢(shì)格靜電相互作用的計(jì)算采用了一個(gè)靜電勢(shì)格點(diǎn)，在格點(diǎn)上儲(chǔ)存受體

27、分子的靜電勢(shì)。當(dāng)點(diǎn)，在格點(diǎn)上儲(chǔ)存受體分子的靜電勢(shì)。當(dāng)配體和受體分子對(duì)接時(shí)，某個(gè)原子和受體配體和受體分子對(duì)接時(shí)，某個(gè)原子和受體之間的靜電相互作用能通過周圍格點(diǎn)上靜之間的靜電相互作用能通過周圍格點(diǎn)上靜電勢(shì)以及原子上的部分電荷就可以計(jì)算得電勢(shì)以及原子上的部分電荷就可以計(jì)算得到。到。 計(jì)算氫鍵相互作用時(shí)，格點(diǎn)的處理和范德計(jì)算氫鍵相互作用時(shí)，格點(diǎn)的處理和范德華相互作用有點(diǎn)類似，每個(gè)格點(diǎn)上需要保華相互作用有點(diǎn)類似，每個(gè)格點(diǎn)上需要保存配體分子中所有氫鍵給體與氫鍵受體之存配體分子中所有氫鍵給體與氫鍵受體之間的相互作用能量，而且這些能量都是在間的相互作用能量，而且這些能量都是在氫鍵在最佳情況下的氫鍵能量值。氫鍵在最佳情況下的氫鍵能量值。 AutoDock格點(diǎn)對(duì)接示意圖 用圍繞受體活性位點(diǎn)的氨基酸殘基形成一個(gè)范圍更大的Box 用不同類型的原子作為探針進(jìn)行掃描，計(jì)算格點(diǎn)能量 AutoDock程序?qū)ε潴w在Box范圍內(nèi)進(jìn)行構(gòu)象搜索，根據(jù)能量評(píng)分對(duì)結(jié)果進(jìn)行排序。 分子對(duì)接在藥物設(shè)計(jì)中的應(yīng)用分子對(duì)接在藥物設(shè)計(jì)中的應(yīng)用一一、HIV-1蛋白酶Kuntz等利用Dock程序研究HIV-1蛋白酶體系，根據(jù)分子相似性對(duì)劍橋晶體數(shù)據(jù)庫進(jìn)行搜尋。然后選擇那些可以與活性位點(diǎn)處的天冬氨酸作用的帶有NH或者OH且易于合成的化合物。DesJ