姚建華 李佳 徐雯麗 蔣舒仰 胡靜 李靖靖 王宇飛

1中國科學院能量調(diào)控材料重點實驗室，中國科學院上海有機化學研究所（上海 200032）

2鄭州工程技術(shù)學院 （河南鄭州 450044）

眾所周知，一個化合物的性質(zhì)與其化學結(jié)構(gòu)式關(guān)系密切。通常做的分子設(shè)計或化合物設(shè)計，本質(zhì)上是要設(shè)計一個具有某種特定功能/性能的化合物。人們關(guān)注的目標是化合物的功能/性能，而非化合物的幾何結(jié)構(gòu)。

國際純粹與應用化學聯(lián)合會(International Union of Pure and Applied Chemistry，（IUPAC） 是世界化學命名、術(shù)語（包括周期表中新元素的命名）、標準化測量方法、原子量和許多其他重要評估數(shù)據(jù)的權(quán)威機構(gòu)，并為全球化學界提出開發(fā)和維護創(chuàng)造共同語言的建議。在IUPAC出版的化學術(shù)語概略（Compendium of Chemical Terminology）[1]書中，分子設(shè)計的定義如下所述：所謂分子設(shè)計，即采用各種技術(shù)去發(fā)現(xiàn)具有潛在應用所需的特定性質(zhì)的新化學實體。計算機輔助分子設(shè)計的定義為：所謂計算機輔助分子設(shè)計，即采用各種計算機輔助技術(shù)，開展發(fā)現(xiàn)、設(shè)計和優(yōu)化具有特定結(jié)構(gòu)和性質(zhì)的化合物的工作。

本文將介紹用于分子設(shè)計的計算機輔助方法。

1 經(jīng)驗型分子設(shè)計

經(jīng)驗型分子設(shè)計是以設(shè)計者的經(jīng)驗及相關(guān)理論為依據(jù)，設(shè)計具有特定功能/性能的化合物，其工作流程如圖1所示。

圖1所示的經(jīng)驗型分子設(shè)計工作流程表明，工作人員根據(jù)靈感和經(jīng)驗，以及已有的研究報道，提出一個可能具有某種功能/特性，且未見報道的化合物結(jié)構(gòu)。采用合成方法獲得該化合物后，用對應的實驗方法，測試該化合物的相關(guān)功能/特性。如果實驗測試結(jié)果符合設(shè)計要求，那么該化合物的分子設(shè)計工作即完成。如果實驗測試結(jié)果不符合設(shè)計要求，設(shè)計者則根據(jù)經(jīng)驗和靈感，修改化合物的化學結(jié)構(gòu)，然后進行合成、實驗測試。如果實驗測試結(jié)果仍不符合設(shè)計要求，則重復前一過程，直到獲得符合設(shè)計要求的化合物。

圖1 經(jīng)驗型分子設(shè)計工作流程示意圖

2 計算機輔助分子設(shè)計

IUPAC關(guān)于計算機輔助分子設(shè)計的內(nèi)容包含三部分：化合物發(fā)現(xiàn)、設(shè)計和化合物結(jié)構(gòu)優(yōu)化。可以認為，所謂化合物發(fā)現(xiàn)是利用計算機輔助技術(shù)在自然界的動物、植物和礦物中發(fā)現(xiàn)具有特定功能/性能的化合物；所謂化合物設(shè)計是利用計算機輔助技術(shù)設(shè)計具有特定功能/性能的化合物；所謂化合物結(jié)構(gòu)優(yōu)化，即利用計算機輔助技術(shù)對那些具有特定功能/性能的化合物結(jié)構(gòu)的取代基團及其位置作適當調(diào)整，以起到提高化合物性能的作用。

計算機輔助分子設(shè)計的方法主要有三種：基于數(shù)據(jù)（D）、基于邏輯（L）和基于原理（P）（如圖 2所示）?；跀?shù)據(jù)即是利用數(shù)據(jù)庫系統(tǒng)，獲得相應的化合物功能、性質(zhì)和化學結(jié)構(gòu)等數(shù)據(jù)；基于邏輯即是利用功能/性質(zhì)與結(jié)構(gòu)之間的關(guān)系及規(guī)則，根據(jù)化合物的化學結(jié)構(gòu)預測它的可能功能/性質(zhì)；基于原理即是利用量化計算方法，根據(jù)化合物的化學結(jié)構(gòu)，計算化學結(jié)構(gòu)對應的相關(guān)物化參數(shù)，并根據(jù)參數(shù)，推測可能的功能。一般而言，人們關(guān)注的化合物功能/性質(zhì)取決于化合物的應用領(lǐng)域，它們的預測策略和方法因應用領(lǐng)域不同而不同。如，醫(yī)藥或農(nóng)藥領(lǐng)域，化合物的生物活性預測，需要同時關(guān)注化合物的化學結(jié)構(gòu)以及作用靶標信息；材料領(lǐng)域，材料的性能不僅與單種成分化合物的性質(zhì)相關(guān)，還與配方中所有化合物的性質(zhì)及它們之間的相互作用相關(guān)。

圖2 三種計算機輔助分子設(shè)計方法

圖2 所示的In house數(shù)據(jù)庫系統(tǒng)，一般為自建的數(shù)據(jù)庫系統(tǒng)。通常，In house數(shù)據(jù)庫系統(tǒng)由兩部分組成：數(shù)據(jù)庫管理系統(tǒng)和數(shù)據(jù)庫。數(shù)據(jù)庫管理系統(tǒng)是用于管理和查詢數(shù)據(jù)庫中的數(shù)據(jù)；數(shù)據(jù)庫中包含的結(jié)構(gòu)化的數(shù)據(jù)，可以是本單位自用的專用實驗與計算研究數(shù)據(jù)，也可以是專門從文獻中采集的化合物及其特定性質(zhì)數(shù)據(jù)。

圖2所示的基于邏輯的性質(zhì)預測系統(tǒng)，具有根據(jù)化合物的化學結(jié)構(gòu)預測其相關(guān)性質(zhì)的功能。而基于原理的方法，通常是利用程序計算一些參數(shù)，根據(jù)這些參數(shù)預測相關(guān)的性質(zhì)。

計算機輔助的分子設(shè)計的流程如圖3所示。

圖3 計算機輔助分子設(shè)計的流程

圖3 顯示的流程表明，計算機輔助分子設(shè)計過程中，先利用數(shù)據(jù)庫系統(tǒng)查詢憑經(jīng)驗和靈感設(shè)計出的化合物是否已存在（即基于數(shù)據(jù)的方法）。如果數(shù)據(jù)庫中已收錄了對應的化合物，則需要根據(jù)預測模型、經(jīng)驗和靈感修改化合物結(jié)構(gòu)；如果該化合物尚未被收錄到數(shù)據(jù)庫中，則需要利用結(jié)構(gòu)與性質(zhì)之間的關(guān)系及規(guī)則，預測該化合物的性質(zhì)；如果預測結(jié)果符合設(shè)計要求，則進入合成階段；如果不符合要求，設(shè)計者則需要根據(jù)預測模型、經(jīng)驗和靈感優(yōu)化化合物結(jié)構(gòu)，利用基于數(shù)據(jù)的方法，判斷該化合物是否被數(shù)據(jù)庫收錄；利用基于邏輯或原理的方法預測該化合物性質(zhì)，并判斷是否符合設(shè)計要求；重復這個過程，直到獲得預測結(jié)果符合設(shè)計要求，然后進入合成階段。之后的工作流程與經(jīng)驗型的基本相同。

自1946年計算機問世，化學數(shù)據(jù)管理系統(tǒng)就是化學家們期盼的文獻和數(shù)據(jù)查詢工具。早在20世紀50年代，美國國家標準化辦公室數(shù)據(jù)處理部門的工作人員，在《科學（Science）》雜志發(fā)表了題為“用數(shù)字化計算機查找化學記錄”的文章[2]。在這篇文章中，作者介紹了化學數(shù)據(jù)、化合物結(jié)構(gòu)計算機處理的一些方法和策略，以及當時已建立的SEAC［NBS（National Bureau of Standards）Electronic Automatic Computer］系統(tǒng)。

由于計算機系統(tǒng)的不斷更新，對應的化學數(shù)據(jù)管理系統(tǒng)一級數(shù)據(jù)庫結(jié)構(gòu)也隨之不斷更新升級。從單機版的化學數(shù)據(jù)庫系統(tǒng)，逐步轉(zhuǎn)化為“Client/Server”，現(xiàn)在已更新為“基于網(wǎng)絡(luò)”的系統(tǒng)?！盎诰W(wǎng)絡(luò)”系統(tǒng)的優(yōu)點是：數(shù)據(jù)及其管理系統(tǒng)的維護及時且方便。目前，在化學及相關(guān)領(lǐng)域常用的商業(yè)化合物數(shù)據(jù)庫系統(tǒng)有：美國化學文摘社的SciFinder[3]，BIOVIA公司的化合物毒性數(shù)據(jù)庫系統(tǒng)、化合物代謝數(shù)據(jù)庫系統(tǒng)、可用化學品目錄數(shù)據(jù)庫系統(tǒng)[4]，愛斯唯爾的Reaxys化學數(shù)據(jù)庫系統(tǒng)等等[5]。In house或非商業(yè)化學數(shù)據(jù)庫系統(tǒng)一般不公開，這里只介紹一下中國科學院上海有機化學研究所姚建華團隊建立的一些非商業(yè)專用的化學數(shù)據(jù)庫系統(tǒng)，如化合物疏水常數(shù)數(shù)據(jù)庫系統(tǒng)、農(nóng)用化學品數(shù)據(jù)庫、食品添加劑信息管理系統(tǒng)和危險化學品信息查詢系統(tǒng)等等[6]。

關(guān)于化合物性質(zhì)預測的工作，在20世紀40年代就有人提出化合物生物活性與它的化學結(jié)構(gòu)相關(guān)的概念[7]。目前，化合物性質(zhì)預測的方法主要有以下三大類：（1）以分子描述符[8]為自變量，性質(zhì)為因變量的方程式作為預測依據(jù)；（2）以實驗得到的規(guī)律作為預測依據(jù)的專家型推理；（3）以分子結(jié)構(gòu)信息與性質(zhì)的對應關(guān)系為預測依據(jù)。

第一類預測方法大多通過已有的軟件計算多種的分子描述符，分析這些描述符與性質(zhì)的相關(guān)性，建立對應的方程式，并用于預測化合物性質(zhì)。這種方法的優(yōu)點：軟件開發(fā)的工作量不大；不足之處：必須依靠計算分子描述符的軟件。

第二類預測方法通過收集已報道的實驗總結(jié)得出的結(jié)構(gòu)與性質(zhì)關(guān)系，比較被預測化合物結(jié)構(gòu)與規(guī)則相符度，預測化合物的性質(zhì)。這種方法的優(yōu)點：有效地利用了實驗結(jié)果；不足之處：實驗規(guī)則數(shù)量的有限性。

第三類預測方法通過收集已報道的實驗結(jié)果，分析化學結(jié)構(gòu)與性質(zhì)之間的關(guān)系，建立化學結(jié)構(gòu)與性質(zhì)的關(guān)系模型，并用于預測化合物性質(zhì)。這種方法的優(yōu)點：有效地利用了實驗結(jié)果及其結(jié)果中隱含的規(guī)則；不足之處：軟件開發(fā)工作量極大。

分子描述符主要包含以下10大類：拓撲類（topological）、 幾 何 類 （geometrical）、 電 子 類（electronic）、物化類 （physicochemical）、指紋類（fingerprints）、官能團類（functional groups）、特性類（properties）、電荷類（charge-related）、半經(jīng)驗類（semi-empirical）和熱力學類（thermodynamical）。描述符計算軟件各自定義描述符類型及其計算方法。

目前，現(xiàn)有的分子描述符的計算軟件主要包括幾類[9]：ADAPT，ADMET Predictor，ADRIANA.Code，CODESSA，DRAGON，GRID，ISIDA，MARVIN Beans，MOE，MOLCONN-Z，MOLGEN-QSPR，PowerMV，HyperChem[10]和Gaussian[11]等等。其中部分軟件還具有描述符與性質(zhì)的關(guān)系數(shù)據(jù)分析功能，如CODESSA和DRAGON等。

第二類預測方法的軟件有DEREK（Deductive Estimation of Risk from Existing Knowledge）[12]，該軟件可用于預測化合物的毒性。

第三類預測方法的軟件有PASS[13]、CISOCPSMT、CISOC-PSAT 和 CISOC-pKa 等[14]。

在計算機輔助分子設(shè)計工作中，常用的預測軟件有：PASS，Sybyl（原 Tripos公司的產(chǎn)品），CoMFA（Comparative Molecular Field Analysis）[15]，CoMSIA（Molecular Similarity Indices in a Comparative Analysis）[16]，DEREK，MultiCASE[17]，CISOC-PSCT，CISOC-PSMT，CISOC-PSAT，CISOC-logP，CISOC-pKa，CISOC-logS[14]，TOPKAT，Discovery Studio 和Material Studio[18]等等。其中，PASS，Sybyl 和Discovery Studio適用于藥物或農(nóng)藥研發(fā)中化合物生物活性的預測；Material Studio適用于材料研發(fā)中化合物物化參數(shù)和材料性能的預測；DEREK，MultiCASE，CISOC-PSCT，CISOC-PSMT，CISOCPSAT和TOPKAT可用于預測化合物的毒性；CISOC-logP，CISOC-pKa和CISOC-logS分別預測化合物的疏水常數(shù)、酸堿解離常數(shù)和水溶解度。

與經(jīng)驗型的分子設(shè)計相比，計算機輔助分子設(shè)計有效提高了分子設(shè)計的效率和精準性，降低了實驗的盲目性和化合物設(shè)計過程中產(chǎn)生的不符合設(shè)計要求的化合物數(shù)量，降低了實驗廢棄物處理工作量。