定量构效关系

**定量構效關係（QSAR）** 是一種藉助分子的理化性質參數或結構參數，以數學和統計學手段定量研究有機小分子與生物大分子相互作用、有機小分子在生物體內吸收、分布、代謝、排泄等生理相關性質的方法。這種方法廣泛應用於藥物、農藥、化學毒劑等生物活性分子的合理設計，在早期的藥物設計中，定量構效關係方法占據主導地位，1990 年代以來隨著計算機計算能力的提高和眾多生物大分子三維結構的準確測定，基於結構的藥物設計逐漸取代了定量構效關係在藥物設計領域的主導地位，但是 QSAR 在藥學研究中仍然發揮著非常重要的作用。

## 發展歷史

定量構效關係是在傳統構效關係的基礎上，結合物理化學中常用的經驗方程的數學方法出現的，其理論歷史可以追溯到 1868 年提出的 Crum-Brown 方程，該方程認為化合物的生理活性可以用化學結構的函數來表示，但是並未建立明確的函數模型。最早的可以實施的定量構效關係方法是美國波蒙拿學院的 Hansch 在 1962 年提出的 Hansch 方程。Hansch 方程脫胎於 1935 年英國物理化學家哈密頓提出的哈密頓方程以及改進的塔夫托方程。哈密頓方程是一個計算取代苯甲酸解離常數的經驗方程，這個方程將取代苯甲酸解離常數的對數值與取代基團的電性參數建立了線性關係，塔夫托方程是在哈密頓方程的基礎上改進形成的計算脂肪族酯類化合物水解反應速率常數的經驗方程，它將速率常數的對數與電性參數和立體參數建立了線性關係。

Hansch 方程在形式上與哈密頓方程和塔夫托方程非常接近，以生理活性物質的半數有效量作為活性參數，以分子的電性參數、立體參數和疏水參數作為線性回歸分析的變量，隨後，Hansch 和日本訪問學者藤田稔夫等人一道改進了 Hansch 方程的數學模型，引入了指示變量、拋物線模型和雙線性模型等修正，使得方程的預測能力有所提高。

幾乎在 Hansch 方法發表的同時，Free 等人發表了 Free-Wilson 方法，這種方法直接以分子結構作為變量對生理活性進行回歸分析。其在藥物化學中的應用範圍遠不如 Hansch 方法廣泛。Hansch 方法、Free-Wilson 方法等方法均是將分子作為一個整體考慮其性質，並不能細緻地反應分子的三維結構與生理活性之間的關係，因而又被稱作**二維定量構效關係**

二維定量構效關係出現之後，在藥物化學領域產生了很大影響，人們對構效關係的認識從傳統的定性水平上升到定量水平。定量的結構活性關係也在一定程度上揭示了藥物分子與生物大分子結合的模式。在 Hansch 方法的指導下，人們成功地設計了諾氟沙星等喹諾酮類抗菌藥。

由於二維定量不能精確描述分子三維結構與生理活性之間的關係，1980 年代前後人們開始探討基於分子構象的**三維定量構效關係**的可行性。1979 年，Crippen 提出「距離幾何學的 3D-QSAR」；1980 年 Hopfinger 等人提出「分子形狀分析方法」；1988 年 Cramer 等人提出了「比較分子場方法」（CoMFA）。比較分子場方法一經提出便席捲藥物設計領域，成為應用最廣泛的基於定量構效關係的藥物設計方法；1990 年代，又出現了在比較分子場方法基礎上改進的「比較分子相似性方法」以及在「距離幾何學的 3D-QSAR」基礎上發展的「虛擬受體方法」等新的三維定量構效關係方法，但是老牌的 CoMFA 依然是使用最廣泛的定量構效關係方法。

## 二維定量構效關係

二維定量構效關係方法是將分子整體的結構性質作為參數，對分子生理活性進行回歸分析，建立化學結構與生理活性相關性模型的一種藥物設計方法，常見的二維定量構效關係方法有 Hansch 方法、Free-wilson 方法、分子連接性方法等，最為著名和應用最廣泛的是 Hansch 方法

### 活性參數

活性參數是構成二維定量構效關係的要素之一，人們根據研究的體系選擇不同的活性參數，常見的活性參數有：半數有效量、半數有效濃度、半數抑菌濃度、半數致死量、最小抑菌濃度等，所有活性參數均必須採用物質的量作為計量單位，以便消除分子量的影響，從而真實地反應分子水平的生理活性。為了獲得較好的數學模型，活性參數在二維定量構效關係中一般取負對數後進行統計分析。

### 結構參數

結構參數是構成定量構效關係的另一大要素，常見的結構參數有：疏水參數、電性參數、立體參數、幾何參數、拓撲參數、理化性質參數以及純粹的結構參數等

* **疏水參數**：藥物在體內吸收和分布的過程與其疏水性密切相關，因而疏水性是影響藥物生理活性的一個重要性質，在二維定量構效關係中採用的疏水參數最常見的是脂水分配係數，其定義為分子在正辛醇與水中分配的比例，對於分子母環上的取代基，脂水分配係數的對數值具有加和性，可以通過簡單的代數計算獲得某一取代結構的疏水參數。
* **電性參數**：二維定量構效關係中的電性參數直接繼承了哈密頓公式和塔夫托公式中的電性參數的定義，用以表徵取代基團對分子整體電子分配的影響，其數值對於取代基也具有加和性。
* **立體參數**：立體參數可以表徵分子內部由於各個基團相互作用對藥效構象產生的影響以及對藥物和生物大分子結合模式產生的影響，常用的立體參數有塔夫托立體參數、摩爾折射率、范德華半徑等。
* **幾何參數**：幾何參數是與分子構象相關的立體參數，因為這類參數常常在定量構效關係中占據一定地位，故而將其與立體參數分割考慮，常見的幾何參數有分子表面積、溶劑可及化表面積、分子體積、多維立體參數等
* **拓撲參數**：在分子連接性方法中使用的結構參數，拓撲參數根據分子的拓撲結構將各個原子編碼，用形成的代碼來表徵分子結構。
* **理化性質參數**：偶極矩、分子光譜數據、前線軌道能級、酸鹼解離常數等理化性質參數有時也用做結構參數參予定量構效關係研究
* **純粹的結構參數**：在 Free-Wilson 方法中，使用純粹的結構參數，這種參數以某一特定結構的分子為參考標準，依照結構母環上功能基團的有無對分子結構進行編碼，進行回歸分析，為每一個功能基團計算出回歸係數，從而獲得定量構效關係模型。

### 數學模型

二維定量構效關係中最常見的數學模型是線性回歸分析，Hansch 方程和 Free-Wilson 方法均採用回歸分析。

經典的 Hansch 方程形式為：

 $\lg \left({\frac {1}{C}}\right)=a\pi +b\sigma +cE_{s}+k$ 其中 $\pi $ 為分子的疏水參數，其與分子脂水分配係數 $P_{X}$ 的關係為： $\pi =\lg \left({\frac {P_{x}}{P_{H}}}\right)$ ， $\sigma $ 為哈密頓電性參數， $E_{s}$ 為塔夫托立體參數，其中 a，b，c，k 均為回歸係數。

日本學者藤田稔夫對經典的 Hansch 方程作出一定改進，用拋物線模型描述疏水性與活性的關係：

 $\lg \left({\frac {1}{C}}\right)=a\pi +b{\pi }^{2}+c\sigma +dE_{s}+k$ 這一模型擬合效果更好。

Hansch 方程進一步，以雙直線模型描述疏水性與活性的關係：

 $\lg \left({\frac {1}{C}}\right)=a\lg P-b\lg(\beta P+1)+D$ 其中的 P 為分子的脂水分配係數， $a,b,\beta $ 為回歸係數，D 代表方程的其他部分。雙直線模型的預測能力比拋物線模型進一步加強。

Free-Wilson 方法的方程形式為：

 $\lg \left({\frac {1}{C}}\right)=\sum _{i}\sum _{j}G_{ij}X_{ij}+\mu $ 其中 $X_{ij}$ 為結構參數，若結構母環中第 i 個位置有第 j 類取代基則結構參數取值為 1 否則為 0， $\mu $ 為參照分子的活性參數， $G_{ij}$ 為回歸係數。

除了回歸分析，遺傳算法、人工神經網絡、偏最小二乘分析、模式識別、單純形方法等統計分析方法也會應用於二維定量構效關係數學模型的建立

### 發展

目前，二維定量構效關係的研究集中在兩個方向：結構數據的改良和統計方法的優化。

傳統的二維定量構效關係使用的結構數據常僅能反應分子整體的性質，通過改良結構參數，使得二維結構參數能夠在一定程度上反應分子在三維空間內的伸展狀況，成為二維定量構效關係的一個發展方向。

引入新的統計方法，如遺傳算法、人工神經網絡、偏最小二乘回歸等，擴展二維定量構效關係能夠模擬的資料結構的範圍，提高 QSAR 模型的預測能力是 2D-QSAR 的主要發展方向。

## 三維定量構效關係

三維定量構效關係是引入了藥物分子三維結構信息進行定量構效關係研究的方法，這種方法間接地反映了藥物分子與大分子相互作用過程中兩者之間的非鍵相互作用特徵，相對於二維定量構效關係有更加明確的物理意義和更豐富的信息量，因而 1980 年代以來，三維定量構效關係逐漸取代了二維定量構效關係的地位，成為基於機理的合理藥物設計的主要方法之一。目前應用最廣泛的三維定量構效關係方法是 CoMFA 和 CoMSIA 即比較分子場方法和比較分子相似性方法，除了上述兩種方法，3D-QSAR 還有 DG 3D-QSAR、MSA、GERM 等眾多方法。

### CoMFA&CoMSIA

CoMFA 和 CoMISA 是應用最廣泛的合理藥物設計方法之一，這種方法認為，藥物分子與受體間的相互作用取決於化合物周圍分子場的差別，以定量化的分子場參數作為變量，對藥物活性進行回歸分析便可以反應藥物與生物大分子之間的相互作用模式進而有選擇地設計新藥。

比較分子場方法將具有相同結構母環的分子在空間中疊合，使其空間取向儘量一致，然後用一個探針粒子在分子周圍的空間中遊走，計算探針粒子與分子之間的相互作用，並記錄下空間不同坐標中相互作用的能量值，從而獲得分子場數據。不同的探針粒子可以探測分子周圍不同性質的分子場，甲烷分子作為探針可以探測立體場，水分子作為探針可以探測疏水場，氫離子作為探針可以探測靜電場等等，一些成熟的比較分子場程序可以提供數十種探針粒子供用戶選擇。

探針粒子探測得到的大量分子場信息作為自變量參與對分子生理活性數據的回歸分析，由於分子場信息數據量很大，屬於高維化學數據，因而在回歸分析過程中必須採取數據降維措施，最常用的方式是偏最小二乘回歸，此外主成分分析也用於數據的分析。

統計分析的結果可以圖形化地輸出在分子表面，用以提示研究者如何有選擇地對先導化合物進行結構改造。右圖為一 CoMFIA 計算的結果輸出，圖中藍色區域若以負電性基團取代則會提高藥物的活性，紅色區域則提示正電性基團更有利於活性。除了直觀的圖形化結果，CoMFA 還能獲得回歸方程，以定量描述分子場與活性的關係。

CoMSIA 是對 CoMFA 方法的改進，他改變了探針粒子與藥物分子相互作用能量的計算公式，從而獲得更好的分子場參數。

### 其他三維定量構效關係方法

除了比較分子場方法，三維定量構效關係還有距離幾何學三位定量構效關係（DG 3D-QSAR）、分子形狀分析（MSA）、虛擬受體等方法（FR）等

距離幾何學三維定量構效關係嚴格來講是一種介於二維和三維之間的 QSAR 方法。這種方法將藥物分子劃分為若干功能區塊定義藥物分子活性位點，計算低能構象時各個活性位點之間的距離，形成距離矩陣；同時定義受體分子的結合位點，獲得結合位點的距離矩陣，通過活性位點和結合位點的匹配為每個分子生成結構參數，對生理活性數據進行統計分析。

分子形狀分析認為藥物分子的藥效構象是決定藥物活性的關鍵，比較作用機理相同的藥物分子的形狀，以各分子間重疊體積等數據作為結構參數進行統計分析獲得構效關係模型。

虛擬受體方法是 DG 3D-QSAR 和 CoMFA 方法的延伸與發展，其基本思路是採用多種探針粒子在藥物分子周圍建立一個虛擬的受體環境，以此研究不同藥物分子之間活性與結構的相關性。其原理較之 CoMFA 方法更加合理，是目前定量構效關係研究的熱點之一。

## 方法評價

定量構效關係研究是人類最早的合理藥物設計方法之一，具有計算量小，預測能力好等優點。在受體結構未知的情況下，定量構效關係方法是最準確和有效地進行藥物設計的方法，根據 QSAR 計算結果的指導藥物化學家可以更有目的性地對生理活性物質進行結構改造。在 1980 年代計算機技術爆炸式發展之前，QSAR 是應用最廣泛也幾乎是唯一的合理藥物設計手段。

但是 QSAR 方法不能明確給出回歸方程的物理意義以及藥物 - 受體間的作用模式，物理意義模糊是對 QSAR 方法最主要的置疑之一。另外在定量構效關係研究中大量使用了實驗數據和統計分析方法，因而 QSAR 方法的預測能力很大程度上受到試驗數據精度的限制，同時時常要面對「統計方法欺詐」的置疑。

## 參見

* 化學相似性
* 構效關係
* 分子對接
* 藥物設計

## 外部連結

* History of QSAR（頁面存檔備份，存於網際網路檔案館）
* QSAR 與藥物設計
* 第十一屆定量構效關係國際研討會
* QSAR（頁面存檔備份，存於網際網路檔案館）

## 參考書目

* 李仁利 1998 年藥物的構效關係北京大學藥學院講義
* 陳凱先等 2000 年計算機輔助藥物設計 —— 原理、方法及應用上海科學技術出版社 ISBN 7-5323-5551-9
* 徐筱傑等 2004 年計算機輔助藥物分子設計化學工業出版社 ISBN 7-5025-5520-X