分子對接
分子對接
分子對接是通過受體的特徵以及受體和藥物分子之間的相互作用方式來進行藥物設計的方法。主要研究分子間(如配體和受體)相互作用,並預測其結合模式和親合力的一種理論模擬方法。近年來,分子對接方法已成為計算機輔助藥物研究領域的一項重要技術。
分子對接這一想法的歷史可以追溯到19世紀提出的受體學說,Fisher提出的受體學說認為,藥物與體內的蛋白質大分子即受體會發生類似鑰匙與鎖的識別關係,這種識別關係主要依賴兩者的空間匹配。隨著受體學說的發展,人們對生理活性分子與生物分子的相互作用有了更加深刻的認識,從基於空間匹配的剛性模型逐漸發展成為基於空間匹配和能量匹配的柔性模型。模型的優化使通過計算模擬分子間相互作用的設想更容易實現。另一方面,計算機和計算科學的迅速發展又使得人們能夠處理大量數據,這兩方面的因素共同促成了分子對接方法的出現。
早期的分子對接方法用分子力學方法或者量子化學方法計算小分子之間分子識別,在一些分子模擬軟體包中也含有分子對接的模塊。但是由於演演算法和計算機處理能力的限制,早期的對接方法較難處理含有大分子的分子對接過程。
1995年由Accelrys公司開發的計算化學軟體Affinity上市,這是第一個可以進行有大分子參與的商業化分子對接軟體,此後,商業化和免費的分子對接軟體層出不窮。現在應用中的分子對接軟體涵蓋了剛性對接、半柔性對接、柔性對接等各種對接方法,在能量優化方面則使用了人工神經網路、遺傳演演算法、模擬退火、禁忌搜索、局部搜索等各種方法,目前的分子對接方法是研究小分子與大分子相互作用模式、生物大分子間識別、分子自組裝、超分子結構等課題的常用方法之一。
分子對接方法的兩大課題是分子之間的空間識別和能量識別。空間匹配是分子間發生相互作用的基礎,能量匹配是分子間保持穩定結合的基礎。對於幾何匹配的計算,通常採用格點計算、片斷生長等方法,能量計算則使用模擬退火、遺傳演演算法等方法。
各種分子對接方法對體系均有一定的簡化,根據簡化的程度和方式,可以將分子對接方法分為三類。
剛性對接:剛性對接方法在計算過程中,參與對接的分子構像不發生變化,僅改變分子的空間位置與姿態,剛性對接方法的簡化程度最高,計算量相對較小,適合於處理大分子之間的對接。
半柔性對接:半柔性對接方法允許對接過程中小分子構像發生一定程度的變化,但通常會固定大分子的構像,另外小分子構像的調整也可能受到一定程度的限制,如固定某些非關鍵部位的鍵長、鍵角等,半柔性對接方法兼顧計算量與模型的預測能力,是應用比較廣泛的對接方法之一。
柔性對接:柔性對接方法在對接過程中允許研究體系的構像發生自由變化,由於變數隨著體系的原子數呈幾何級數增長,因此柔性對接方法的計算量非常大,消耗計算機時很多,適合精確考察分子間識別情況。
Dock是目前應用最廣泛的分子對接軟體之一,由Kuntz課題組開發。Dock應用半柔性對接方法,固定小分子的鍵長和鍵角,將小分子配體拆分成若干剛性片斷,根據受體表面的幾何性質,將小分子的剛性片斷重新組合,進行構像搜索。在能量計算方面,Dock考慮了靜電相互作用、范德華力等非鍵相互作用,在進行構像搜索的過程中搜索體系勢能面。最終軟體以能量評分和原子接觸罰分之和作為對接結果的評價依據。
Autodock是另外一個應用廣泛的分子對接程序,由Olson科研組開發。AutoDock應用半柔性對接方法,允許小分子的構像發生變化,以結合自由能作為評價對接結果的依據。自從AutoDock3.0版本以後,對能量的優化採用拉馬克遺傳演演算法(LGA),LGA將遺傳演演算法與局部搜索方法相結合,以遺傳演演算法迅速搜索勢能面,用局部搜索方法對勢能面進行精細的優化。
FlexX是德國國家信息技術研究中心生物信息學演演算法和科學計算研究室開發的分子對接軟體,目前已經作為分子設計軟體包BioSolveIT LeadIT的一個模塊實現商業化。FlexX使用碎片生長的方法尋找最佳構像,根據對接自由能的數值選擇最佳構像。FlexX程序對接速度快效率高,可以用於小分子資料庫的虛擬篩選。