應(yīng)用ab initio方法和ABEEMσπ模型研究含鋁金屬酶＊

楊忠志，寧方達(dá)

（遼寧師范大學(xué)化學(xué)化工學(xué)院，遼寧大連 116029）

應(yīng)用ab initio方法和ABEEMσπ模型研究含鋁金屬酶＊

楊忠志，寧方達(dá)

（遼寧師范大學(xué)化學(xué)化工學(xué)院，遼寧大連 116029）

應(yīng)用原子－鍵電負(fù)性均衡方法中的σπ模型（ABEEMσπ模型），通過大量量子化學(xué)計(jì)算，擬合確定了含鋁金屬酶體系的ABEEMσπ參數(shù).將這些參數(shù)應(yīng)用到含鋁金屬酶大分子體系的電荷分布及Fukui函數(shù)的計(jì)算，結(jié)果顯示，ABEEMσπ模型計(jì)算得到的電荷分布及Fukui函數(shù)與從頭算和實(shí)驗(yàn)結(jié)論均有很好的一致性.還進(jìn)一步計(jì)算分析了1L3R酶與絲氨酸結(jié)合前后的分子各區(qū)域的電荷分布，結(jié)果表明，Al3＋是1L3R酶的活性中心，根據(jù)結(jié)合后分子的Fukui函數(shù)可以得出絲氨酸會(huì)使1L3R酶的活性降低.另外，通過比較兩者結(jié)合前后Al3＋的廣義Fukui函數(shù)，證明了廣義Fukui函數(shù)可用于該體系分子間反應(yīng)活性的比較，同時(shí)也說明利用ABEEMσπ模型來預(yù)測(cè)含鋁金屬酶的抑制劑是可行的.

從頭算方法；ABEEMσπ模型；電荷分布；Fukui函數(shù)

分子的電荷分布對(duì)于解釋和預(yù)測(cè)其結(jié)構(gòu)和性質(zhì)非常重要［1］.含鋁金屬酶是受人們關(guān)注的重要生命物質(zhì)，因體系大，對(duì)它的研究較困難.筆者應(yīng)用ABEEMσπ模型計(jì)算了多個(gè)含鋁金屬酶的原子、化學(xué)鍵和孤對(duì)電子區(qū)域的電荷分布，這些分子的ABEEMσπ電荷分布與HF／STO-3G電荷分布的線性相關(guān)系數(shù)均在0.96以上，表明ABEEMσπ模型可以用于準(zhǔn)確地計(jì)算大分子的電荷分布.

Parr和Yang提出的Fukui函數(shù)，使得前線分子軌道（FMO）理論和軟硬酸堿（HSAB）理論聯(lián)系起來.Geerlings等［2］提出應(yīng)用Fukui函數(shù)研究分子內(nèi)的反應(yīng)活性.筆者的計(jì)算結(jié)果表明，Al3＋均是所研究的含鋁金屬酶的活性中心，這與實(shí)驗(yàn)結(jié)論相一致，說明ABEEMσπ模型可以用于準(zhǔn)確地確定含鋁金屬酶的活性中心，此方法可以向其他金屬酶推廣.

生物體內(nèi)的一些小分子可以降低酶的反應(yīng)活性.筆者計(jì)算了當(dāng)1L3R酶與絲氨酸結(jié)合后的Fukui函數(shù)和電荷變化，也計(jì)算了Al3＋的廣義Fukui函數(shù).發(fā)現(xiàn)兩者結(jié)合后，1L3R酶的活性變小，說明ABEEMσπ模型可應(yīng)用于預(yù)測(cè)含鋁金屬酶的抑制劑.

1 計(jì)算方法

1.1 從頭計(jì)算法

利用從頭計(jì)算方法，基組的選擇非常重要［3－4］.大的基組在計(jì)算能量和優(yōu)化結(jié)構(gòu)時(shí)準(zhǔn)確度高，但計(jì)算電荷分布時(shí)用小的基組更好，因?yàn)楫?dāng)較大的基組應(yīng)用于某原子時(shí)，一定程度上也會(huì)將與其相鄰的原子包含在內(nèi)，這樣就會(huì)對(duì)分布在該原子上的Mulliken電荷估計(jì)過高.Wilson M S等［5］認(rèn)為，采用大的基組計(jì)算分子的電荷分布時(shí)，對(duì)分子內(nèi)原子間的電荷轉(zhuǎn)移的計(jì)算值會(huì)過高.Derouane E G等［6］研究發(fā)現(xiàn)，6-21G基組比STO-3G基組計(jì)算的電荷值高.Jakalian A等［7］指出，6-31G＊基組的計(jì)算值比STO-3G基組的計(jì)算值高10%～15%.筆者利用STO-3G基組計(jì)算分子的電荷分布.

1.2 ABEEMσπ模型

由電負(fù)性均衡原理和密度泛函理論，Yang等提出并發(fā)展了原子與鍵電負(fù)性均衡方法（ABEEM）的σπ模型（ABEEMσπ模型）.在該模型中，雙鍵被劃分為1個(gè)σ鍵區(qū)域和4個(gè)π鍵區(qū)域，σ鍵區(qū)域的電荷中心選擇在成鍵原子的共價(jià)半徑之比處，π鍵區(qū)域的電荷中心選擇在垂直于雙鍵所在的平面并位于雙鍵原子上下兩側(cè).

根據(jù)電負(fù)性均衡原理，分子中各區(qū)域的有效電負(fù)性與分子的電負(fù)性相等.這對(duì)于任意包含有i個(gè)單鍵原子、j個(gè)雙鍵原子、k個(gè)單鍵、l個(gè)雙鍵中的σ鍵區(qū)域、4l個(gè)π鍵區(qū)域和m個(gè)孤對(duì)電子區(qū)域的分子，就會(huì)同時(shí)有（i＋j＋k＋l＋4l＋m）個(gè)方程.再聯(lián)立分子的總電荷守恒方程，則能夠得到分子的電負(fù)性和分子中每一個(gè)原子、鍵和孤對(duì)電子的電荷分布.

2 結(jié)果與討論

2.1 含鋁金屬酶電荷分布的計(jì)算

計(jì)算大分子的電荷分布，首先要選擇模型分子，選用HF／STO-3G方法計(jì)算電荷，然后利用最小二乘法，運(yùn)用相關(guān)程序，擬合Al3＋等的參數(shù)，即價(jià)態(tài)電負(fù)性和價(jià)態(tài)硬度.該程序中，基點(diǎn)的選擇非常關(guān)鍵.基點(diǎn)是指程序中所需要用到的初始電荷.筆者在Yang等的研究基礎(chǔ)上擬合并確定了含鋁金屬酶體系的參數(shù)，部分基點(diǎn)和ABEEMσπ參數(shù)如表1所示.

表1 含鋁金屬酶體系的部分基點(diǎn)和ABEEMσπ參數(shù)

研究發(fā)現(xiàn)，ABEEMσπ模型計(jì)算的模型分子的電荷分布可以很好地模擬從頭算的結(jié)果，因此文中所擬合的ABEEMσπ參數(shù)是合理的.ABEEMσπ模型中的原子、化學(xué)鍵和孤對(duì)電子類型是根據(jù)其周圍的化學(xué)環(huán)境恰當(dāng)定義的，因此它們的參數(shù)是可轉(zhuǎn)移的.運(yùn)用所擬合的參數(shù)，自編程序，計(jì)算PDB代號(hào)為1BS1，1L3R，3UKD的含鋁金屬酶的電荷分布［8－10］.1BS1的分子式為C26H48AlF3N9O16P2，共含有105個(gè)原子.1L3R的分子式為C26H57AlF3N8O31P5，共含有131個(gè)原子.3UKD的分子式為C25H43AlF3N12O20P3，共含有107個(gè)原子.對(duì)于這3種含鋁金屬酶，ABEEMσπ模型計(jì)算的電荷分布與HF／STO-3G方法計(jì)算的電荷分布的線性相關(guān)如圖1所示.

圖1 ABEEMσπ模型與HF／STO-3G方法計(jì)算的含鋁金屬酶的電荷分布的線性相關(guān)

圖1表明，利用ABEEMσπ模型計(jì)算的3個(gè)大分子回歸到原子上的電荷分布，與HF／STO-3G方法計(jì)算的電荷分布的線性相關(guān)方程的斜率和線性相關(guān)系數(shù)均接近于1，截距均接近于0.這說明ABEEMσπ模型可以用于準(zhǔn)確計(jì)算大分子的電荷分布.

2.2 含鋁金屬酶活性中心的研究

利用本研究所編程序，可計(jì)算得到含鋁金屬酶各區(qū)域的Fukui函數(shù).所研究分子的骨架結(jié)構(gòu)及標(biāo)號(hào)如圖2所示.經(jīng)程序計(jì)算所得的其活性區(qū)域的Fukui函數(shù)值列于表2—4.從以表2—4可以看出，所研究的3個(gè)酶中，Al3＋的Fukui函數(shù)值均最大，因此這些酶的活性中心均為Al3＋，這與實(shí)驗(yàn)上的結(jié)論相一致［8－10］.由此可見，ABEEMσπ模型可以準(zhǔn)確地確定含鋁金屬酶的活性中心.

圖2 含鋁金屬酶的骨架結(jié)構(gòu)及標(biāo)號(hào)

表2 1BS1的Fukui函數(shù)

表3 1L3R的Fukui函數(shù)

表4 3UKD的Fukui函數(shù)

2.3 含鋁金屬酶的抑制劑的研究

筆者考察了絲氨酸1L3R酶活性的抑制作用.絲氨酸與PDB代號(hào)為1L3R的酶結(jié)合后的骨架結(jié)構(gòu)如圖3所示，表5給出了兩者結(jié)合前后的電荷分布.由表5可知，與絲氨酸結(jié)合后，相比較于兩者結(jié)合之前，Al3＋的電荷減小了0.122e，變化值大于其他位點(diǎn)的電荷變化，所以，絲氨酸與1L3R酶結(jié)合時(shí)Al3＋的電荷所受的影響最大，可見，Al3＋確實(shí)是1L3R酶的活性中心.

表6給出1L3R酶分子與絲氨酸結(jié)合后的Fukui函數(shù).由表6可以看出，與絲氨酸結(jié)合后，Al3＋的Fukui函數(shù)值不再是最大的，也就是說，1L3R酶的活性中心Al3＋的活性降低.因此，ABEEMσπ模型可以用于預(yù)測(cè)含鋁金屬酶的抑制劑.

近來，楊忠志等［11］提出應(yīng)用廣義活性指標(biāo)作為研究分子間的反應(yīng)活性指標(biāo)，并指出，廣義Fukui函數(shù)越大，則反應(yīng)活性也就越大.筆者改進(jìn)程序，經(jīng)計(jì)算可得，1L3R酶與絲氨酸結(jié)合前后，Al3＋的廣義Fukui函數(shù)分別為1.325 0和0.677 5，可見，結(jié)合絲氨酸以后，1L3R酶的活性明顯變小，也驗(yàn)證了廣義活性指標(biāo)可作為研究分子間的反應(yīng)活性指標(biāo).

圖3 1L3R和絲氨酸結(jié)合后的骨架結(jié)構(gòu)

表5 1L3R和絲氨酸結(jié)合前后的電荷分布

表6 1L3R和絲氨酸結(jié)合后的Fukui函數(shù)

3 結(jié)語(yǔ)

應(yīng)用ABEEMσπ模型計(jì)算了多個(gè)含鋁金屬酶的電荷分布，與HF／STO-3G的電荷分布結(jié)果均有很好的一致性，而且該方法更省時(shí).應(yīng)用ABEEMσπ模型計(jì)算的Fukui函數(shù)表明，Al3＋均是所研究的含鋁金屬酶的活性中心，這與實(shí)驗(yàn)結(jié)論相一致，說明ABEEMσπ模型可用于準(zhǔn)確地確定含鋁金屬酶的活性中心.此外，計(jì)算了當(dāng)PDB代號(hào)為1L3R的含鋁金屬酶與絲氨酸結(jié)合后的Fukui函數(shù)和電荷變化，結(jié)果表明，兩者結(jié)合后1L3R酶的活性變小，同時(shí)，驗(yàn)證了廣義Fukui函數(shù)可用于比較分子間的反應(yīng)活性，這表明ABEEMσπ模型也適用于預(yù)測(cè)含鋁金屬酶的抑制劑.該研究為進(jìn)一步作含鋁金屬酶分子的動(dòng)力學(xué)模擬研究奠定基礎(chǔ)，也為ABEEMσπ模型應(yīng)用于其他金屬酶開了先例.

［1］ YANG Z Z，WANG C S.Atom-Bond Electronegativity Equalization Method.1.Calculation of the Charge Distribution in Large Molecules［J］.J.Phys.Chem.A，1997，101：6 315－6 321.

［2］ GEERLINGS P，PROFT F D，LANGENAEKER W.Conceptual Density Functional Theory［J］.Chem.Rev.，2003，103：1 793－1 873.

［3］ KANG Y K，SCHERAGA H A.An Efficient Method for Calculating Atomic Charges of Peptides and Proteins from E-lectronic Populations［J］.J.Phys.Chem.B，2008，112：5 470－5 478.

［4］ REED A E，WEINSTOCK R B，WEINHOLD F.Natural Population Analysis［J］.J.Chem.Phys.，1985，83：735－746.

［5］ WILSON M S，ICHIKAWA S.Comparison Between the Geometric and Harmonic Mean Electronegativity Equilibration Techniques［J］.J.Phys.Chem.，1989，93：3 087－3 089.

［6］ DEROUANE E G，F(xiàn)RIPIAT J G，BALLMOOS R V.Quantum Mechanical Calculations on Molecular Sieves.2.Model Cluster Investigation of Sillcoaluminophosphates［J］.J.Phys.Chem.，1990，94：1 687－1 692.

［7］ JAKALIAN A，BUSH B L，JACK D B，et al.Fast，Efficient Generation of High-Quality Atomic Charges.AM1-BCC Model：Ⅰ.Method［J］.J.Comput.Chem.，2000，21（2）：132－146.

［8］ K?CK H，SANDMARK J，GIBSON K J，et al.Crystal Structure of Two Quaternary Complexes of Dethiobiotin Synthetase，Enzyme-MgADP-AlF3－Diaminopelargonic Acid and Enzyme-MgADP-Dethiobiotin-Phosphate；Implications for Catalysis［J］.Protein Science，1998，7：2 560－2 566.

［9］ MADHUSUDAN，AKAMINE P，XUONG N H，et al.Crystal Structure of a Transition State Mimic of the Catalytic Subunit of cAMP-Dependent Protein Kinase［J］.Nature Structural Biology，2002，9（4）：273－277.

［10］ SCHLICHTING I，REINSTEIN J.Structures of Active Conformations of UMP Kinase from Dictyostelium Discoideum Suggest Phosphoryl Transfer is Associative［J］.Biochemistry，1997，36：9 290－9 296.

［11］ 楊忠志，徐珍珍.應(yīng)用Fukui函數(shù)探討雙苯基－取代的自由基閉環(huán)反應(yīng)的區(qū)位選擇性［J］.遼寧師范大學(xué)學(xué)報(bào)：自然科學(xué)版，2012，35（2）：193－196.

（責(zé)任編輯 向陽(yáng)潔）

Study on Aluminum Metalloenzymes by ab initio Method and ABEEMσπModel

YANG Zhong-zhi，NING Fang-da
（School of Chemistry and Chemical Engineering，Liaoning Normal University，Dalian 116029，Liaoning China）

By applying the atom-bond electronegativity equalizationσπmodel（ABEEMσπmodel），a large number of quantum chemistry calculations were performed to determine the ABEEMσπparameters of aluminum metalloenzymes.Then these parameters were employed to study the charge distributions and Fukui function of aluminum metalloenzymes.Calculated results obtained by ABEEMσπmodel are in good agreement with those by the ab initio method and experimental conclusions.Further analysis of the charge distributions between 1L3Rand 1L3R-serine indicates that Al3＋is the active center of 1L3R，and serine would reduce the activity of 1L3Raccording to Fukui function.In addition，comparing the generalized Fukui function of Al3＋in 1L3Rwith that in 1L3R-serine，generalized Fukui function is proved to be appropriate in this system，and ABEEMσπmodel is feasible to predict the inhibitors of aluminum metalloenzymes.

ab initio method；ABEEMσπmodel；charge distributions；Fukui function

O641

A

10.3969／j.issn.1007－2985.2013.05.018

1007－2985（2013）05－0074－05

2013－04－19

國(guó)家自然科學(xué)基金資助項(xiàng)目（21133005）

楊忠志（1940－），男，吉林舒蘭人，遼寧師范大學(xué)化學(xué)化工學(xué)院教授，博士，博士生導(dǎo)師，主要從事理論與計(jì)算化學(xué)及其應(yīng)用等研究.