定量構效關係

定量構效關係（quantitative structure-activity relationship，QSAR）又稱量化構效關係，是利用理論計算和統計分析方法研究系列化合物分子結構與生物活性之間所存在的定量關係，可作為定量預測未知物質性質的一種方法。

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

定量構效關係模型（QSAR models）是化學和生物科學與工程中使用的迴歸或分類模型。與其他迴歸模型一樣，QSAR 迴歸模型將一組預測變數（X）與響應變數（Y）的效力相關聯，而 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 方法

2.1 活性引數

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

2.2 結構引數

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

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

2.3 數學模型

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

經典的 Hansch 方程形式為：

其中為分子的疏水引數，其與分子脂水分配係數的關係為：，為哈密頓電性引數，為塔夫託立體引數，其中 a，b，c，k 均為迴歸係數。

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

這一模型擬合效果更好。

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

其中的 P 為分子的脂水分配係數，為迴歸係數，D 代表方程的其他部分。雙直線模型的預測能力比拋物線模型進一步加強。

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

其中為結構引數，若結構母環中第 i 個位置有第 j 類取代基則結構引數取值為 1 否則為 0，為參照分子的活性引數，為迴歸係數。

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

2.4 發展

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

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

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

三維定量構效關係

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

3.1 CoMFA&CoMSIA

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

作為探針可以探測靜電場等等，一些成熟的比較分子場程式可以提供數十種探針粒子供使用者選擇。

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

區域則提示正電性基團更有利於活性。除了直觀的圖形化結果，CoMFA 還能獲得迴歸方程，以定量描述分子場與活性的關係。

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

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

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

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

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

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

方法評價

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

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