沒食子酸及其衍生物抗氧化活性的理論研究

裴 玲，張海霞，張 巖

?

沒食子酸及其衍生物抗氧化活性的理論研究

*裴 玲1,2，張海霞1,2，張 巖1,2

（1.濱州學院化學工程系，山東，濱州 256603；2.濱州市材料化學重點實驗室，山東，濱州 256603）

采用量子化學密度泛函理論(DFT)中的B3LYP方法在6-311+G(d,p)水平下對沒食子酸(GA)及其衍生物沒食子酸甲酯(MG)、沒食子酸乙酯(EG)、沒食子酸丙酯(PG)、沒食子酸異丙酯(IPG)、沒食子酸丁酯(BG)和沒食子酸叔丁酯(TBG)進行了優化計算，從7個分子的幾何結構、氫原子的NBO電荷、羥基解離能、電離勢、HOMO、LUMO、ΔE(LUMO-HOMO)等方面分析了7種物質清除自由基的可能機理。結果表明，7種物質C(2)位的酚羥基為最大可能活性位點，容易發生抽氫反應，其次是C(1)位和C(3)位。抽氫反應表明，在非極性溶劑中其抗氧化能力大小為TBG > BG > IPG > PG > EG > MG > GA；電子轉移機理表明，極性溶劑中其抗氧化能力大小為TBG > IPG > BG ≈ PG ≈ EG > MG > GA。

沒食子酸；抗氧化性；解離能；密度泛函

多酚類物質是一類廣泛存在于植物體內的次生代謝物質的混合物，具有眾多的生物活性，尤其是在清除自由基能力方面表現了很好的抗氧化性[1-4]，成為近年來備受關注的具有廣泛開發前景的天然抗氧化劑。其抗氧化結構-活性相關性研究引起了人們研究的興趣，并取得了一些研究成果[5-8]。沒食子酸(GA)及其衍生物是植物界中分布較廣的一類多酚類物質，具有抗氧化性，作為常見的油溶性抗氧化劑，已被廣泛應用于油脂及含油食品、油炸食品、肉類食品、花粉食品、化妝品等。馬耀等人[9]采用了循環伏安掃描和電子自旋共振相結合的方法對沒食子酸(GA)及其衍生物沒食子酸甲酯(MG)、沒食子酸乙酯(EG)、沒食子酸丙酯(PG)、沒食子酸異丙酯(IPG)、沒食子酸丁酯(BG)和沒食子酸叔丁酯(TBG)的電子轉移能力、清除自由基的能力和臨床應用之間的關系進行了研究，結果表明GA類化合物能有效清除自由基，故而顯示抗氧化性。為了進一步深入了解此類物質的抗氧化機理，本研究用GA類化合物清除自由基能力的大小代表其抗氧化性的相對強弱，運用密度泛函方法計算GA類化合物的結構性質，從理論上對其抗氧化活性進行詳細的理論分析，以期能從理論的角度，深化GA類物質的抗氧化性機理研究，為能更好的開發利用天然植物源多酚類提供理論依據。

1 計算方法

構建GA及其衍生物的初始構象后，運用密度泛函B3LYP方法，在6-311＋G(d,p)基組的水平下，對GA類化合物及其脫氫自由基的幾何構型進行優化和能量計算，得到優勢構象，頻率分析表明無虛頻，為勢能面上的穩定駐點。所有的計算通過Gaussian 03程序完成[10]。

2 結果與討論

2.1 GA類化合物及自由基的優化

2.1.1 GA類化合物的優化

根據理論分析三個羥基成順式時，相鄰兩OH可形成弱氫鍵，此時化合物較穩定，但三個OH的取向可能會影響化合物的穩定性，為了全面考慮化合物的穩定性，我們優化了三個OH成順式取向的兩種情況，兩種情況下優化的GA類化合物幾何構型見圖1。不同構型下優化的化合物的能量、相對能量列于表1。由表1的數據可知，各個化合物構型I和構型II的能量相差2 kJ·mol-1左右，說明兩種構型都是其穩定結構。

圖1 GA及其衍生物的幾何構型

表1 GA及其衍生物的能量

2.1.2 GA類化合物自由基的優化

由于酚類抗氧劑在清除自由基過程中主要是通過抽氫反應產生較穩定的半醌自由基來終止鏈式反應，因此我們主要考慮生成半醌自由基的情況，考慮到形成自由基后是否存在氫鍵，優化了7種物質可能形成的半醌式自由基，以PG為例畫出形成自由基時有無氫鍵的結構，結果見圖2。各個物質自由基的能量列于表2中。

酚類化合物清除自由基的活性，取決于環上羥基抽氫反應的難易和抽氫后生成自由基的穩定性。從每種分子自由基的能量來看，三個酚羥基各去掉一個H原子后形成的自由基，有兩個氫鍵的自由基(1-OH-a、2-OH-a、3-OH-a)最穩定，其次是含有一個氫鍵的自由基(1-OH-b、1-OH-c、2-OH-b、2-OH-c、3-OH-c)，最不穩定的是不含氫鍵的自由基(2-OH-d)。但含有一個氫鍵自由基時，OH的H與相鄰C=O雙鍵的O形成氫鍵的自由基比OH的H與相鄰OH的O形成氫鍵的自由基更穩定，如1-OH-b比1-OH-c穩定，3-OH-c比3-OH-b穩定。

表2 各個自由基的能量(單位：a.u.)

2.2 反應活性位點確定

2.2.1 GA類化合物及其自由基的結構參數

根據表1知，每種分子的結構I的能量較低，較穩定，表3列出了各化合物構型I的主要結構參數。酚類化合物清除自由基活性與酚羥基鍵強度成反比，O-H鍵長越長，鍵強度越低，則抗氧化劑活性越高。從表2可以看出，7種物質的RO(2)-H(2)> RO(1)-H(1)> RO(3)-H(3)，說明C(2)位的酚羥基鍵最弱，即7種物質的C(2)位的羥基氫原子可能更容易失去。

表3 GA及其衍生物的主要幾何結構參數(鍵長單位：nm)

鑒于各自由基的穩定性，后面討論半醌式自由基時主要以形成的雙氫鍵自由基作為基礎來討論。比較表2含有兩個氫鍵的自由基，發現各個物質的2-OH-a自由基能量最低，其次是1-OH-a，3-OH-a自由基的能量最高，據此我們可以初步判斷當C(2)上OH失去H原子后形成的自由基最穩定。

2.2.2 NBO 電荷布居分析

親電基團傾向于進攻負電荷多的原子，親核基團傾向于進攻正電荷多的原子。人體內自由基多是氧自由基，其上的氧原子帶有較大的負電荷。酚類化合物清除自由基的最主要機制是酚羥基與自由基發生抽氫反應。所以酚羥基H原子所帶正電荷越正，與自由基上帶負電荷的氧原子吸引作用越強，就越容易與自由基發生反應。因此帶較大正電荷的酚羥基H原子就是自由基反應的最大可能活性位點。表4列出了7種分子的酚羥基H的NBO電荷。從表4可以看出，7種分子的C(2)位的羥基H 原子所帶正電荷稍大，所以2位C 原子上的羥基氫原子與自由基作用稍大，2位C原子上的羥基為最大可能活性位點。

表4 羥基H原子上的NBO電荷布局數

2.2.3 O-H鍵解離焓與電離勢

目前普遍認為[11-12]在非極性溶劑中，抗氧化劑傾向于通過一步抽氫反應清除自由基，以鍵解離焓(BDE)為理論指標；在極性溶劑中，抗氧化劑清除自由基主要通過電子轉移反應進行，以電離勢(IP)為理論指標。BDE為理論指標對應抗氧化劑清除自由基的一步抽氫機理，BDE越低說明O-H 鍵越弱，當與自由基反應時，越容易斷裂，抽氫反應越容易進行，相應的抗氧化劑活性可能越強。IP為理論指標對應抗氧化劑清除自由基的電子轉移機理，IP越小越容易發生電子轉移，相應的抗氧化劑活性可能越強。鍵解離焓和電離勢可根據文獻[13]公式、來計算，其中是母體分子失去氫原子后產生的自由基的焓，是氫原子的焓，是母體分子的焓，是相應的自由基陽離子的生成焓。計算時設定溫度為標準溫度，校正因子v = 0.9804[14]。表5列出了計算得到的GA類化合物的不同羥基的脫氫解離焓和電離勢。

從表5的結果發現，當GA類化合物生成半醌自由基時，比較這3種斷鍵方式，發現C(2)上O-H鍵BDE最小，其次是C(1)上O-H和C(3)上O-H鍵的BDE，分別比2號位要高大約25 kJ·mol-1和27 kJ·mol-1，IP也有類似現象。由此可以判斷生成3種可能的半醌自由基由易到難依次為：2-OH、1-OH、3-OH自由基。此結論和表2中自由基能量分析的結果是一致的，2號位在容易發生反應，體現出較好的抗氧化性。主要原因在于2號位上發生抽氫反應以后，形成更加穩定的分子內氫鍵，從而使其產物的穩定性較好，容易發生反應。

表5 GA及其衍生物穩定自由基的OH鍵解離能和電離勢(單位：kJ·mol-1)

2.3 抗氧化活性比較

2.3.1 一步抽氫機理

當GA類化合物形成半醌自由基時，各種物質的BDE大小順序為：TBG < BG < IPG < PG < EG < MG < GA，由此可以判定在非極性溶劑中7種物質的抗氧化活性為TBG > BG > IPG > PG > EG > MG > GA。影響BDE的一個很重要的因素是分子內的氫鍵效應，分子內氫鍵可以穩定化合物，也可以穩定自由基，但對后者的穩定作用更大。另一個因素是電子效應，給電子取代基降低O-H的BDE，吸電子取代基升高O-H的BDE，7種基團(H, CH3, C2H5, C3H7, CH(CH3)2,C4H9, C(CH3)3)的給電子能力依次增強，因此其降低O-H的BDE依次增強，結果導致7種物質的BDE大小為TBG < BG < IPG < PG < EG < MG < GA。

3.3.2 電子轉移機理

（1）電離勢分析

根據表5中的數據可以看出，在形成半醌自由基時，各物質的IP大小順序為：TBG < IPG < BG < PG < EG < MG < GA，由此可以判斷在極性溶劑中7種分子的抗氧化性為TBG > IPG > BG > PG > EG > MG > GA。

（2）前線分子軌道及軌道能級分析

前線分子軌道理論認為分子在反應過程中優先起作用的是前線分子軌道。它對于分子的化學性質特別是其反應活性有著決定性作用。一般地講，HOMO對電子束縛較為松弛，具有電子給予體的性質，表現為EHOMO越大，其上電子受核吸引力越小，電子越不穩定，越易失去?？寡趸瘷C理中參數IP與最高占據軌道(HOMO)的能級密切相關。在以IP為指標的單電子轉移機理中，從母體化合物中轉移一個電子的難易程度，直接決定了化合物抗氧化能力的強弱。而對母體化合物而言，其最高占據軌道(HOMO)能級越高，軌道中的電子受到原子核的束縛越松弛，就越容易發生質子伴隨的電子轉移反應，分子的抗氧化性就越好。7個分子的HOMO數值見表6。從表6中可看到，7種分子的HOMO能級順序為TBG > IPG > BG ≈ PG ≈ EG > MG > GA，由此可以判斷在極性溶劑中，抗氧化性TBG優于IPG，BG、PG、EG三物質的抗氧化性相差不大，都比MG、GA的抗氧化性強。另外，對比7種分子的前線軌道能量和其去掉H原子以后的自由基前線軌道能量，可以看出失去H原子后，導致其電子最高占據軌道能級降低，從而使其失電子能力降低，即還原性降低。從理論上得出7種分子分子具有還原性的特征是由于它容易失去H 自由基。

表6 GA及其衍生物的最高占據軌道能 (a.u.)

通過電離勢和前線分子軌道分析知，在極性溶劑中抗氧化性TBG優于IPG，BG、PG、EG三物質的抗氧化性相差不大，都比MG、GA強。區別于抽氫反應機理，IPG的抗氧化性優于BG。

3 結論

采用B3LYP/6-311＋G(d,p)方法對7種GA類化合物抗氧化活性進行了研究，計算結果表明，7種物質在生成半醌自由基時2-OH活性最高，羥基上的氫最容易被含氧自由基提取，最有可能參與自由基的清除，這主要是由于GA類化合物C(2)位相鄰羥基之間的弱氫鍵相互作用，這種額外的共軛穩定性以及協同作用增強了2-羥基的活性，這跟馬耀等[9]的實驗研究結果是一致的。同時比較7種物質的抗氧化活性發現，在極性溶劑中7種化合物的抗氧化活性：TBG > BG > IPG > PG > EG > MG > GA，非極性溶劑中抗氧化性為TBG > IPG > BG ≈ PG ≈ EG > MG > GA。GA類化合物的抗氧化活性不同，主要受苯環上羰基取代基性質的影響。

[1] Sun Y M,Zhang H Y,Chen D Z,et al. Theoretical elucidation on the antioxidant mechanism of curcumin: a DFT study[J]. Organic letters, 2002, 4(17): 2909-2911.

[2] Lu L,Zhu S,Zhang H, et al. Improvement of antioxidative activity of resveratrol by elongating conjugated chain: A DFT theoretical study[J]. Computational and Theoretical Chemistry, 2013, 1019(1):39-47.

[3] Moo-Huchin V M, Moo-Huchin M I, Estrada-León R J, et al. Antioxidant compounds, antioxidant activity and phenolic content in peel from three tropical fruits from Yucatan, Mexico[J]. Food chemistry, 2015, 166: 17-22.

[4] Deepha V, Praveena R, Sivakumar R, et al. Experimental and theoretical investigations on the antioxidant activity of isoorientin from Crotalaria globosa[J]. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, 2014, 121: 737-745.

[5] 李新芳,王學業,劉萬強,等. 酚類抗氧化劑結構—活性的ab initio研究[J]. 計算機與應用化學, 2005, 22(4):287-290.

[6] 張立偉,陳世榮,楊頻. 肉桂酸類抗氧化劑結構—活性關系研究[J]. 結構化學, 2003, 22(3), 341-345.

[7] 謝湖均,雷群芳,方文軍. 槲皮素抗氧化活性的密度泛函理論研究[J]. 化學學報, 2010, 68(15):1467-1472.

[8] Hou C Y. Theoretical study of antioxidative ability and antioxidative mechanism of norathyriol in solution[J]. Computational and Theoretical Chemistry, 2014, 1028:87-91.

[9] 馬耀,張建中,等. 沒食子酸及其衍生物抗氧化機制的研究[J]. 科學通報, 1993, 38(15):1411-1414.

[10] Frisch M J,Trucks G W,Schlegel H B,et al. Gaussian03-Version c01; Gaussian Inc[J]. Wallingford, CT, 2004.

[11] Zhang H Y, Wang L F, Sun Y M. Why B-ring is the active center for genistein to scavenge peroxyl radical: a DFT study[J]. Bioorganic & medicinal chemistry letters, 2003, 13(5): 909-911.

[12] Wang L F, Zhang H Y. A theoretical investigation on DPPH radical-scavenging mechanism of edaravone[J]. Bioorganic & medicinal chemistry letters, 2003, 13(21): 3789-3792.

[13] 張軍,陳德展. 橄欖油中兩種酚類化合物的抗氧化性理論研究[J]. 山東師范大學學報:自然科學版, 2005, 20(4):36-37.

[14] 張紅雨,王蘭芬. 茶多酚清除過氧自由基分子機理的量子化學計算研究[J]. 淄博學院學報:自然科學與工程版, 2002, 4(3):5-10.

Theory Study on the Antioxidation Activity of Gallic Acid and its Derivatives

*PEI Ling1,2, ZHANG Hai-xia1,2, ZHANG Yan1,2

(1.Department of Chemical Engineering, Binzhou University, Binzhou, Shandong 256603, China;2. Binzhou Key Laboratory of Material Chemistry, Binzhou, Shandong 256603, China)

The geometrical structures of gallic acid and its derivatives including methyl gallate (MG), ethyl gallate (EG), propyl gallate (PG), isopropyl gallate (IPG), butyl gallate (BG) and tertbutyl gallate (TBG) were calculated by using the Density Functional Theory (DFT) B3LYP method at the level of 6-311+G (d, p). The possible mechanism of scavenging free radicals of the seven compounds were discussed in detail based on the analysis of molecular structures, NBO charges on hydrogen atoms, the dissociation energy and ionization potential of O-H bonds, the energies of the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO), the energy gap between HOMO and LUMO, and etc. The calculated results showed that C2-OH should be the most active site for eliminating free radicals, next as the C1-OH and C3-OH. Pumping hydrogen reaction indicated that the order of the antioxidative activities of seven compounds in non-polar solvent listed as: TBG > BG > IPG > PG > EG > MG > GA. Electron transfer mechanism indicated that the order in polar solvent listed as: TBG > IPG > BG ≈ PG ≈ EG > MG > GA.

gallic acid; antioxidation; dissociation enthalpy; DFT

1674-8085(2015)04-0027-07

O 641

A

10.3969/j.issn.1674-8085.2015.04.006

2015-04-04；修改日期：2015-06-01

山東省濱州市科技發展計劃項目（2014ZC0306），濱州學院科研基金項目（BZXYG1210）

*裴 玲(1980-)，女，山東鄒平人，講師，碩士，主要從事量子化學理論計算研究(E-mail:peiling1201@163.com);

張海霞(1988-)，女，山東無棣人，助教，碩士，主要從事膠體與界面化學研究(Email:angle00521@163.com);

張 巖(1983-)，女，山東棗莊人，講師，主要從事天然有機物的提取研究(E-mail:yanzhang19@163.com).