黃曉峰+曹同成+陳波

摘要：通過對2013年諾貝爾化學(xué)獎的成果介紹，敘述了理論和計算化學(xué)學(xué)科的建立和發(fā)展歷程，指出了計算化學(xué)在新世紀(jì)對化學(xué)學(xué)科理論發(fā)展，以及指導(dǎo)化學(xué)實驗方面具有的重要意義。

關(guān)鍵詞：2013年諾貝爾化學(xué)獎；理論與計算化學(xué)；計算機輔助；模型化學(xué)

文章編號：1005–6629（2014）3–0011–04 中圖分類號：G633.8 文獻標(biāo)識碼：B

2013年的諾貝爾化學(xué)獎被授予了Martin Karplus、Michael Levitt以及Arieh Warshel三位美國科學(xué)家，以表彰他們在發(fā)展復(fù)雜化學(xué)體系多尺度模型方面所做出的杰出貢獻。我們知道，長久以來，化學(xué)學(xué)科的奠基和發(fā)展始終離不開化學(xué)家在實驗室中的辛勤勞動，但與此同時，隨著實踐知識的不斷豐富和完善，以及運算能力的突飛猛進，理論和計算化學(xué)有可能也應(yīng)當(dāng)在新世紀(jì)在化學(xué)學(xué)科的傳統(tǒng)領(lǐng)域發(fā)揮更大的作用。當(dāng)前，解開每個人生命背后的謎團也是人們的興趣所在。Karplus，Levitt和Warshel三位科學(xué)家將經(jīng)典力學(xué)模擬方法結(jié)合最新發(fā)展的量子物理計算方法，為建立和發(fā)展多尺度復(fù)雜模型的理論模擬研究做出了基礎(chǔ)性的貢獻。那么，到底什么是理論模擬方法？它有什么重要的科學(xué)意義？對我們又有什么啟迪？

1 理論與計算化學(xué)的建立和發(fā)展歷程

20世紀(jì)初量子力學(xué)的發(fā)現(xiàn)為科學(xué)家們打開了深層次研究分子和原子的大門。量子力學(xué)中著名的薛定諤方程以其優(yōu)美簡潔的形式描述了原子和分子的重要組成部分——電子的行為[1]。1927年，Walter Heitler以及Fritz London兩位科學(xué)家利用薛定諤方程解開了氫氣分子電子結(jié)構(gòu)[2]，理論化學(xué)從此悄然興起。隨后，價鍵理論[3]、Hartree-Fock理論[4]、分子軌道理論[5]等的建立極大地豐富了理論化學(xué)的內(nèi)容。從此，化學(xué)學(xué)科可以說與物理學(xué)一樣，開始了真正的兩條腿走路，而不再只是依靠實驗知識的獲取跛足而行。

早在20世紀(jì)50年代，科學(xué)家利用半經(jīng)驗的方法對原子軌道進行了計算。50至60年代期間，各種各樣基于現(xiàn)代量子理論的計算已經(jīng)被用來計算一些簡單分子的電子結(jié)構(gòu)和相互作用。20世紀(jì)70年代，例如Gaussian？、ATMOL？、IBMOL？等量子化學(xué)計算軟件的開發(fā)也擴充了計算化學(xué)的內(nèi)涵。

與此同時，新的化學(xué)合成與表征技術(shù)的開發(fā)使得越來越多新穎的分子被制造出來，人們不僅需要認識這些新分子，而且也需要借助一定手段來指導(dǎo)新分子的合成。在這樣的前提下，就需要借助計算機對分子進行模擬。

1990年，密度泛函理論（Density Functional Theory）的提出將理論和計算化學(xué)帶到了一個新紀(jì)元。和以往的方法相比，密度泛函理論解決了以往的分子模型中電子交換和相關(guān)作用的近似，由其得出的分子幾何結(jié)構(gòu)和電子結(jié)構(gòu)的預(yù)測與實驗數(shù)據(jù)吻合得非常好。直至目前，密度泛函理論依然是分子和化學(xué)反應(yīng)模擬中最重要也是最為常用的方法，兩位科學(xué)家Walter Kohn[6]和John Pople[7]因為分別發(fā)展了密度泛函理論以及將這種量子力學(xué)計算方法融入到計算化學(xué)中去而獲得了1998年的諾貝爾化學(xué)獎，這是諾貝爾化學(xué)獎第一次被授予理論和計算化學(xué)領(lǐng)域的科學(xué)家。獲獎?wù)咧坏腜ople也是著名量子化學(xué)計算軟件Gaussian[8]的開發(fā)者之一，該軟件在2009年又進行了一次更新，是當(dāng)今功能最完善、計算最有效、生命力最長的量子化學(xué)計算軟件。

目前，專門刊登量子化學(xué)理論、模型化學(xué)和計算化學(xué)的學(xué)術(shù)期刊也紛紛涌現(xiàn)，如，美國化學(xué)會（American Chemistry Society）下已有Journal of Chemical Information and Modeling， Journal of Chemical Theory and Computation， Journal of Physical Chemistry A三本期刊出版，而著名學(xué)術(shù)出版集團Elsevier也有Journal of Molecular Graphics and Modeling， Journal of Molecular Modeling， International Journal of Quantum Chemistry和Computational and Theoretical Chemistry等專刊，國內(nèi)也有例如《物理化學(xué)學(xué)報》和《計算機及應(yīng)用化學(xué)》等期刊。

2 復(fù)雜化學(xué)體系多尺度模型的建模以及應(yīng)用

1976年，Michael Levitt和Arieh Warshel二人提出了酶催化生物化學(xué)反應(yīng)的通用理論研究方法[10]。這個方法將生物酶-底物間的復(fù)合物和溶劑作用一起考慮在整個體系之內(nèi)，并且用量子力學(xué)和經(jīng)典力學(xué)兩種方法探討了所有可能影響催化路徑的因素。其中，量子力學(xué)包含了酶-底物鍵的斷裂，底物與酶結(jié)合時電荷的重新分布；而經(jīng)典力學(xué)部分則考慮了酶和底物之間的立體作用能和靜電作用能。綜合考慮以上兩點，兩位作者以一種水解酶裂解糖苷鍵為實例，首次進行了水解酶-糖苷這個復(fù)雜化學(xué)體系多尺度模型的理論計算（圖1）。如今復(fù)雜化學(xué)體系的QM/MM方法已經(jīng)被廣泛應(yīng)用到酶-底物催化反應(yīng)，有機反應(yīng)以及DNA/RNA的相關(guān)研究中去。

那么，如何建立一個合理的多尺度復(fù)雜模型？科學(xué)家們和軟件工程師們通力合作開發(fā)出了各種功能強大的分子建模和可視化軟件。對于小分子的構(gòu)建，最為常用的為PerkinElmer公司下屬的劍橋軟件公司開發(fā)的ChemBioOffice？系列軟件，包括了ChemBioDraw？和ChemBio 3D？兩個模塊（圖2）。當(dāng)在軟件窗口的右側(cè)ChemDraw？面板畫出感興趣的分子后，左邊的窗口就會立即顯示出分子的3D模型。本軟件還包括了其他很多內(nèi)容，例如對分子進行簡單的幾何結(jié)構(gòu)優(yōu)化操作或者分子動力學(xué)計算，根據(jù)計算結(jié)果畫出分子的部分電荷、分子軌道等信息。

GaussView？是Gaussian公司開發(fā)的用于分子建模的軟件包，目前已經(jīng)更新到GaussView5.0b版本。此軟件包的功能類似于ChemBioOffice？，該軟件并不如ChemBioOffice？那樣還具有計算功能，而只是作為量子化學(xué)計算軟件Gaussian？的圖形輸入接口，圖3是利用GaussView？創(chuàng)建了聯(lián)苯分子，當(dāng)利用Gaussian？軟件對分子進行計算完畢之后，也能夠展示分子軌道的圖形。

以上兩種軟件不僅可以在各自的軟件內(nèi)部進行計算，而且ChemBioOffice？軟件還提供了Gaussian？計算軟件的接口。我們可以在ChemBioOffice？中構(gòu)建完小分子，并設(shè)置運行參數(shù)之后在Gaussian？中進行對應(yīng)的計算。

在一個復(fù)雜化學(xué)體系中，往往還要涉及到生物大分子的構(gòu)建?，F(xiàn)在科學(xué)家們已經(jīng)構(gòu)建起了大分子結(jié)構(gòu)庫，最著名就是由美國布魯克海文（Brookhaven）國家實驗室建立的蛋白質(zhì)數(shù)據(jù)庫（Protein Data Bank，http：// www.rcsb.org）。庫內(nèi)包含了蛋白質(zhì)、多肽、DNA、RNA等95644個晶體結(jié)構(gòu)數(shù)據(jù)。我們可以通過下載數(shù)據(jù)來得到生物大分子的晶體結(jié)構(gòu)。

Accelrys公司開發(fā)的Discovery Studio Client？軟件能夠讀取從Protein Data Bank下載的pdb文件，如圖4展示的是Discovery Studio Client？的界面，展示了人體血清白蛋白和一種DNA的結(jié)構(gòu)。

此外，Discovery Studio Client？還具有將小分子和大分子組裝結(jié)合在一起的功能，如圖5分別是將一種長鏈的污染物分子結(jié)合到了脂肪酸結(jié)合酶和人體血清白蛋白中，這就完成了一個復(fù)雜化學(xué)體系的模型構(gòu)建。

VMD？軟件也是一種常用的可視化軟件，相對于Discovery Studio Client？，其功能更側(cè)重于動態(tài)展現(xiàn)動力學(xué)情況下分子的運動和形變情況。圖6則是VMD？軟件的界面以及其展示的人體血清白蛋白分子和DNA分子。

在分子建模完成之后，就可以對一個建立完成的化學(xué)體系進行理論的計算，預(yù)測這個復(fù)雜化學(xué)體系的物理化學(xué)性質(zhì)。對于一個多尺度模型的計算，計算方法的選擇也是多尺度的。首先，對需要模擬的化學(xué)反應(yīng)的區(qū)域要進行界定。在界定了這個區(qū)域之后，必須對這個區(qū)域內(nèi)的分子進行高精度的量子化學(xué)計算，模擬或預(yù)測該區(qū)域內(nèi)可能存在的化學(xué)鍵以及鍵的斷裂。在界定的反應(yīng)區(qū)域之外，由于不牽涉到化學(xué)反應(yīng)，所以不需要高精度的量子化學(xué)計算方法，而只需要相對簡單的半經(jīng)驗的計算方法或者更簡單的分子力學(xué)方法進行計算?？偠灾?，這就是復(fù)雜化學(xué)體系多尺度模型的計算，即QM/MM計算。涉及量子化學(xué)部分的QM計算，需要用到包含量子化學(xué)計算的軟件，例如最著名的Gaussian？，GAMESS？等。在這些軟件中，也可以采用ONIOM方法[12]進行計算。

3 復(fù)雜化學(xué)體系多尺度模型建立的科學(xué)意義及其展望

結(jié)合理論以及計算化學(xué)發(fā)展本身的歷程來看，復(fù)雜化學(xué)體系多尺度模型具有十分重要的科學(xué)意義。首先，此模型的建立使我們從簡單分子的化學(xué)反應(yīng)進入到了生物大分子體系的理論計算研究。利用理論計算這個強有力的工具，生命科學(xué)的奧秘將很快被解開，人們對生命科學(xué)背后的化學(xué)機制的認識將會上升到分子層面，對帶動化學(xué)，乃至生命科學(xué)學(xué)科具有舉足輕重的作用。其次，多尺度模型的建立也能夠促進理論和計算化學(xué)本身的發(fā)展，豐富理論和計算化學(xué)本身的內(nèi)涵，并且隨著研究體系的進一步復(fù)雜化，將在現(xiàn)有的多尺度基礎(chǔ)上提出新的超尺度模型的可能。

此外，作為一門交叉學(xué)科，理論和計算化學(xué)的發(fā)展也勢必會帶動其他相關(guān)學(xué)科的進一步發(fā)展。90年代開始，納米學(xué)科蓬勃發(fā)展，各種新材料如雨后春筍般出現(xiàn)，得益于理論化學(xué)中平面波和贗勢（即將離子實的內(nèi)部勢能用假想的勢能取代真實的勢能，但在求解波動方程時，不改變能量本征值和離子實之間區(qū)域的波函數(shù)）的發(fā)展，對具有周期性結(jié)構(gòu)的晶體材料性質(zhì)的模擬和預(yù)測也成為可能。目前，已經(jīng)有Material Studio？、VASP？等多種模擬軟件。在藥物合成方面，計算機輔助藥物合成的概念已經(jīng)深入人心（Computer-aided Drug Design）。顧名思義，計算機輔助藥物設(shè)計利用計算化學(xué)這個強有力的工具來發(fā)現(xiàn)或者研究具有生物活性的藥物分子的行為，其最基本的目標(biāo)就是通過計算化學(xué)來預(yù)測一個分子與靶生物分子是否會結(jié)合，并且其結(jié)合能力有多強，能夠?qū)崿F(xiàn)這一功能的軟件則包括了GOLD？、SYBYL？等等。

可以說，理論和計算化學(xué)已經(jīng)成為輔助化學(xué)家們探索世界的重要工具，也成為了指引科學(xué)家探索未知世界的新羅盤。

參考文獻：

[1] Schrodinger E. An Undulatory Theory of the Mechanics of Atoms and Molecules. Phys. Rev.， 1926， 28， 1049～1070.

[2] Heitler， W. & London， F. Wechselwirkung Neutraler Atome und hom？opolare Bindung nach der Quantenmechanik. Zeitschrift fur Physic 1927， （44）： 455～472.

[3] Pauling， L. Electronic Structure of the Benzene Molecule. Nature， 1987： 325， 396.

[4] Levine， I. N. Quantum Chemistry （4th edition）， Englewood Cliffes， New Jersey： Prentice Hall.

[5] Mulliken， R. S. Electronic Structures of Polyatomic Molecules and Valence. II. General Considerations. Phys. Rev. 1932， （41）： 49～71.

[6] Hohenberg P.； Kohn， W. Inhomogeneous Electron Gas. Phys. Rev. 1964， （136）： B864～B871.

[7] Pople， J. A. Molecular Association in Liquids： II. A Theory of the Structure of Water. Proc. Royal Soc. A， 1951： 205， 163.

[8] Gaussian 09， Revision D.01， M. J. Frisch， G. W. Trucks， H. B. Schlegel， G. E. Scuseria， M. A. Robb， J. R. Cheeseman， G. Scalmani， V. Barone， B. Mennucci， G. A. Petersson， H. Nakatsuji， M. Caricato， X. Li， H. P. Hratchian， A. F. Izmaylov， J. Bloino， G. Zheng， J. L. Sonnenberg， M. Hada， M. Ehara， K. Toyota， R. Fukuda， J. Hasegawa， M. Ishida， T. Nakajima， Y. Honda， O. Kitao， H. Nakai， T. Vreven， J. A. Montgomery， Jr.， J. E. Peralta， F. Ogliaro， M. Bearpark， J. J. Heyd， E. Brothers， K. N. Kudin， V. N. Staroverov， R. Kobayashi， J. Normand， K. Raghavachari， A. Rendell， J. C. Burant， S. S. Iyengar， J. Tomasi， M. Cossi， N. Rega， J. M. Millam， M. Klene， J. E. Knox， J. B. Cross， V. Bakken， C. Adamo， J. Jaramillo， R. Gomperts， R. E. Stratmann， O. Yazyev， A. J. Austin， R. Cammi， C. Pomelli， J. W. Ochterski， R. L. Martin， K. Morokuma， V. G. Zakrzewski， G. A. Voth， P. Salvador， J. J. Dannenberg， S. Dapprich， A. D. Daniels， ？. Farkas， J. B. Foresman， J. V. Ortiz， J. Cioslowski， and D. J. Fox， Gaussian， Inc.， Wallingford CT， 2009.

[9] CHARMM： The Biomolecular Simulation Program， J. Comp. Chem. 30， 1545～1615 （2009）， by B. R. Brooks， C. L. Brooks III， A. D. Mackerell， L. Nilsson， R. J. Petrella， B. Roux， Y. Won， G. Archontis， C. Bartels， S. Boresch A. Caflisch， L. Caves， Q. Cui， A. R. Dinner， M. Feig， S. Fischer， J. Gao， M. Hodoscek， W. Im， K. Kuczera， T. Lazaridis， J. Ma， V. Ovchinnikov， E. Paci， R. W. Pastor， C. B. Post， J. Z. Pu， M. Schaefer， B. Tidor， R. M. Venable， H. L. Woodcock， X. Wu， W. Yang， D. M. York， and M. Karplus.

[10][11] Warshel， A.， Levitt， M. Theoretical Studies of Enzymic Reactions： Dielectric， Electrostatic and Stabilization of the Carbonium Ion in the Reaction of Lysozyme. J. Mol. Biol.， 1976： 103， 227～249.

[12] Svensson， M.， Humber， S.， Froese， R. D. J.， Matsubara， T.， Sieber， S.， Morokuma， K. Oniom. A Multilayered Intergrated MO+MM Method for Geometry Optimizations and Single Point Energy Predictions. A Test for Diels-Alder Reactions and Pt（P（t-Bu）3）2+H2 Oxidative Addition. J. Phys. Chem.， 1996， （100）： 19357～19363.

[5] Mulliken， R. S. Electronic Structures of Polyatomic Molecules and Valence. II. General Considerations. Phys. Rev. 1932， （41）： 49～71.

[6] Hohenberg P.； Kohn， W. Inhomogeneous Electron Gas. Phys. Rev. 1964， （136）： B864～B871.

[7] Pople， J. A. Molecular Association in Liquids： II. A Theory of the Structure of Water. Proc. Royal Soc. A， 1951： 205， 163.

[8] Gaussian 09， Revision D.01， M. J. Frisch， G. W. Trucks， H. B. Schlegel， G. E. Scuseria， M. A. Robb， J. R. Cheeseman， G. Scalmani， V. Barone， B. Mennucci， G. A. Petersson， H. Nakatsuji， M. Caricato， X. Li， H. P. Hratchian， A. F. Izmaylov， J. Bloino， G. Zheng， J. L. Sonnenberg， M. Hada， M. Ehara， K. Toyota， R. Fukuda， J. Hasegawa， M. Ishida， T. Nakajima， Y. Honda， O. Kitao， H. Nakai， T. Vreven， J. A. Montgomery， Jr.， J. E. Peralta， F. Ogliaro， M. Bearpark， J. J. Heyd， E. Brothers， K. N. Kudin， V. N. Staroverov， R. Kobayashi， J. Normand， K. Raghavachari， A. Rendell， J. C. Burant， S. S. Iyengar， J. Tomasi， M. Cossi， N. Rega， J. M. Millam， M. Klene， J. E. Knox， J. B. Cross， V. Bakken， C. Adamo， J. Jaramillo， R. Gomperts， R. E. Stratmann， O. Yazyev， A. J. Austin， R. Cammi， C. Pomelli， J. W. Ochterski， R. L. Martin， K. Morokuma， V. G. Zakrzewski， G. A. Voth， P. Salvador， J. J. Dannenberg， S. Dapprich， A. D. Daniels， ？. Farkas， J. B. Foresman， J. V. Ortiz， J. Cioslowski， and D. J. Fox， Gaussian， Inc.， Wallingford CT， 2009.

[9] CHARMM： The Biomolecular Simulation Program， J. Comp. Chem. 30， 1545～1615 （2009）， by B. R. Brooks， C. L. Brooks III， A. D. Mackerell， L. Nilsson， R. J. Petrella， B. Roux， Y. Won， G. Archontis， C. Bartels， S. Boresch A. Caflisch， L. Caves， Q. Cui， A. R. Dinner， M. Feig， S. Fischer， J. Gao， M. Hodoscek， W. Im， K. Kuczera， T. Lazaridis， J. Ma， V. Ovchinnikov， E. Paci， R. W. Pastor， C. B. Post， J. Z. Pu， M. Schaefer， B. Tidor， R. M. Venable， H. L. Woodcock， X. Wu， W. Yang， D. M. York， and M. Karplus.

[10][11] Warshel， A.， Levitt， M. Theoretical Studies of Enzymic Reactions： Dielectric， Electrostatic and Stabilization of the Carbonium Ion in the Reaction of Lysozyme. J. Mol. Biol.， 1976： 103， 227～249.

[12] Svensson， M.， Humber， S.， Froese， R. D. J.， Matsubara， T.， Sieber， S.， Morokuma， K. Oniom. A Multilayered Intergrated MO+MM Method for Geometry Optimizations and Single Point Energy Predictions. A Test for Diels-Alder Reactions and Pt（P（t-Bu）3）2+H2 Oxidative Addition. J. Phys. Chem.， 1996， （100）： 19357～19363.

[5] Mulliken， R. S. Electronic Structures of Polyatomic Molecules and Valence. II. General Considerations. Phys. Rev. 1932， （41）： 49～71.

[6] Hohenberg P.； Kohn， W. Inhomogeneous Electron Gas. Phys. Rev. 1964， （136）： B864～B871.

[7] Pople， J. A. Molecular Association in Liquids： II. A Theory of the Structure of Water. Proc. Royal Soc. A， 1951： 205， 163.

[8] Gaussian 09， Revision D.01， M. J. Frisch， G. W. Trucks， H. B. Schlegel， G. E. Scuseria， M. A. Robb， J. R. Cheeseman， G. Scalmani， V. Barone， B. Mennucci， G. A. Petersson， H. Nakatsuji， M. Caricato， X. Li， H. P. Hratchian， A. F. Izmaylov， J. Bloino， G. Zheng， J. L. Sonnenberg， M. Hada， M. Ehara， K. Toyota， R. Fukuda， J. Hasegawa， M. Ishida， T. Nakajima， Y. Honda， O. Kitao， H. Nakai， T. Vreven， J. A. Montgomery， Jr.， J. E. Peralta， F. Ogliaro， M. Bearpark， J. J. Heyd， E. Brothers， K. N. Kudin， V. N. Staroverov， R. Kobayashi， J. Normand， K. Raghavachari， A. Rendell， J. C. Burant， S. S. Iyengar， J. Tomasi， M. Cossi， N. Rega， J. M. Millam， M. Klene， J. E. Knox， J. B. Cross， V. Bakken， C. Adamo， J. Jaramillo， R. Gomperts， R. E. Stratmann， O. Yazyev， A. J. Austin， R. Cammi， C. Pomelli， J. W. Ochterski， R. L. Martin， K. Morokuma， V. G. Zakrzewski， G. A. Voth， P. Salvador， J. J. Dannenberg， S. Dapprich， A. D. Daniels， ？. Farkas， J. B. Foresman， J. V. Ortiz， J. Cioslowski， and D. J. Fox， Gaussian， Inc.， Wallingford CT， 2009.

[9] CHARMM： The Biomolecular Simulation Program， J. Comp. Chem. 30， 1545～1615 （2009）， by B. R. Brooks， C. L. Brooks III， A. D. Mackerell， L. Nilsson， R. J. Petrella， B. Roux， Y. Won， G. Archontis， C. Bartels， S. Boresch A. Caflisch， L. Caves， Q. Cui， A. R. Dinner， M. Feig， S. Fischer， J. Gao， M. Hodoscek， W. Im， K. Kuczera， T. Lazaridis， J. Ma， V. Ovchinnikov， E. Paci， R. W. Pastor， C. B. Post， J. Z. Pu， M. Schaefer， B. Tidor， R. M. Venable， H. L. Woodcock， X. Wu， W. Yang， D. M. York， and M. Karplus.

[10][11] Warshel， A.， Levitt， M. Theoretical Studies of Enzymic Reactions： Dielectric， Electrostatic and Stabilization of the Carbonium Ion in the Reaction of Lysozyme. J. Mol. Biol.， 1976： 103， 227～249.

[12] Svensson， M.， Humber， S.， Froese， R. D. J.， Matsubara， T.， Sieber， S.， Morokuma， K. Oniom. A Multilayered Intergrated MO+MM Method for Geometry Optimizations and Single Point Energy Predictions. A Test for Diels-Alder Reactions and Pt（P（t-Bu）3）2+H2 Oxidative Addition. J. Phys. Chem.， 1996， （100）： 19357～19363.