劉靖麗 權(quán)彥 閆浩

摘 ? ? ?要：基于密度泛函（DFT）方法研究甘草中黃酮類化合物的抗氧化活性，對甘草素、甘草苷、異甘草素和異甘草苷4種化合物的分子結(jié)構(gòu)及其相應(yīng)的自由基進行了理論計算。分別在氣相和溶劑中獲得了4種化合物及其相應(yīng)自由基的最優(yōu)幾何結(jié)構(gòu)、NBO電荷、BDE和IP等量子化學參數(shù)。計算結(jié)果顯示：4種化合物的抗氧化活性由大到小的順序為異甘草素、異甘草苷、甘草素、甘草苷，與實驗活性順序一致，并且4種化合物中C7位是自由基反應(yīng)的最大可能活性位點。

關(guān) ?鍵 ?詞：甘草;黃酮類化合物;抗氧化活性;密度泛函理論

中圖分類號：O641 ? ? ? 文獻標識碼： A ? ? ? 文章編號： 1671-0460（2020）10-2163-05

Abstract： In order to study the relationship between the structure and antioxidant activity of flavonoids from licorice， the density functional theory （DFT） B3LYP method with 6-31+G** basis sets was used to optimize the structures of liquiritigenin， liquiritin， isoliquiritigenin， isoliquiritoside and the corresponding free radicals. The properties of the stable geometries were discussed on basis of the molecular geometries， the natural bond order （NBO） charges， the bond dissociation enthalpies （BDE）， ionization potential （IP） and the frontier molecular orbitals theory. The results revealed that the sequence of antioxidant activity of these four flavonoid compounds was： isoliquiritigenin> isoliquiritoside>liquiritigenin>liquiritin， which was consistent with the experimental conclusions. The calculated BDE values for flavonoid compounds clearly demonstrated that the site of C7—OH was the most active site.

Key words： Licorice; Flavonoid compound; Antioxidation activity; DFT

甘草為豆科植物甘草（Glvcyrrhiza uralensis Fisch.）的根及根狀莖，屬多年生草本植物，是一種補益中草藥，具有補脾益氣，祛痰止咳，清熱解毒，調(diào)和諸藥等功效[1]。近年來，相關(guān)學者對甘草進行了大量的系統(tǒng)研究，其主要的化學成分為黃酮類、三萜皂苷類、香豆素類、生物堿及多糖等，這些物質(zhì)具有廣泛的藥理活性[2]。

通過對甘草黃酮類成分的藥理研究發(fā)現(xiàn)，甘草中黃酮類成分具有抗氧化、抗?jié)?、抗腫瘤、抗HIV等作用[3-4]，因而具有潛在的藥用開發(fā)價值。醫(yī)學研究表明[5]，活性氧自由基與一些疾病（如衰老、心血管疾病、炎癥、腫瘤等）的發(fā)生有關(guān)，因而開發(fā)天然抗氧化劑對預(yù)防這些疾病發(fā)揮著至關(guān)重要的作用。傅乃武等[6]研究了甘草黃酮類成分對活性氧的清除作用，實驗表明甘草黃酮類化合物具有顯著的抗氧化活性。隨后崔譽蓉[7]等測定了甘草中的4種黃酮類成分甘草素、異甘草素、甘草苷和異甘草苷（分子結(jié)構(gòu)見圖1）的體外抗氧化活性，并初步探討了4種化合物的結(jié)構(gòu)與抗氧化活性的關(guān)系。結(jié)果表明，4種化合物的抗氧化活性順序為異甘草素、異甘草苷、甘草素和甘草苷。課題組前期合成了異甘草素，研究了其體外抗氧化活性[8]。潘燕[9]等發(fā)現(xiàn)甘草水溶性總黃酮具有抗心肌缺血的活性，其機理是水溶性總黃酮可通過抑制自由基的生成，從而減少自由基對人體的危害，保護心肌受損。LEE[10]等通過實驗發(fā)現(xiàn)甘草的甲醇提取物可保護細胞，免受氧化傷害，因而具有抗氧化的能力。AMAROWICZ[11]等通過體外抗氧化實驗發(fā)現(xiàn)，甘草的乙醇提取物對DPPH·和HO·自由基有清除作用，但是不同溶劑提取物對化合物抗氧化活性的影響未知。

隨著計算化學的發(fā)展，密度泛函（DFT）方法被用來研究分子的抗氧活性。WRIGHT等最先采用DFT方法計算了多酚類抗氧化劑的量子化學參數(shù)酚羥基解離能（BDE）和絕熱電離勢（IP）[12]。本課題組前期運用密度泛函理論方法研究了燈盞花乙素及其苷元[13]和黃芪中黃酮類化合物的抗氧化活 ?性[14]，同時模擬了溶劑對抗氧化活性的影響。雖然張鑫[15]等采用DFT方法研究了甘草素、異甘草素、甘草苷和異甘草苷的抗氧化活性，通過量子化學參數(shù)初步分析和比較了4個化合物的抗氧化活性大小，但是他們只考慮了抽氫反應(yīng)機理，而沒有考慮電子轉(zhuǎn)移機理，也沒有考慮溶劑效應(yīng)對化合物抗氧化機理的影響。而在中草藥有效成分黃酮的提取和分離過程中，甲醇和乙醇是常用的溶劑，可能會對有效成分的活性產(chǎn)生影響，進而影響抗氧劑的抗氧化機理。為此，本文選取甘草中的4種黃酮類化合物甘草素、異甘草素、甘草苷和異甘草苷為研究對象，比較直接的抽氫反應(yīng)機理和電子轉(zhuǎn)移機理，同時考慮溶劑效應(yīng)對抗氧化活性的影響，從原子水平的角度深入揭示4種甘草黃酮類化合物的結(jié)構(gòu)與抗氧化活性關(guān)系，以期對中藥甘草的開發(fā)利用和生物活性研究提供理論上的幫助。

1 ?計算方法

采用DFT（B3LYP）方法以6-31+G**為基組對4種黃酮類化合物及其自由基的分子結(jié)構(gòu)進行全優(yōu)化，對于閉殼層的分子選用限制性的計算，而開殼層的自由基用非限制性的計算。分子的幾何平衡構(gòu)型通過頻率計算進行驗證，沒有虛頻，說明得到的是能量極小值。最后對化合物及其相應(yīng)的脫氫自由基和陽離子自由基進行熱力學和NBO電荷計算，基于4種化合物酚羥基的O—H鍵長、NBO電荷、鍵解離能（BDE）、電離勢值（IP）和前線分子軌道等量子化學參數(shù)，分析4種甘草黃酮類化合物的抗氧化活性和結(jié)構(gòu)之間的關(guān)系及規(guī)律。為了評估溶劑對化合物抗氧化活性的影響，以甲醇、乙醇和水為溶劑，采用PCM模型[16]對化合物及其自由基在水溶液中的結(jié)構(gòu)、頻率和NBO進行了計算。所有計算均由Gaussian 16軟件包完成[17]。

2 ?結(jié)果與討論

2.1 ?酚羥基的鍵長

對于黃酮類化合物來說，它們是通過斷裂 ? ?O—H鍵而起抗氧化作用的，故而根據(jù)分子中O—H的數(shù)目和鍵長可以初步推斷其抗氧化活性的強弱。通過對4種化合物的結(jié)構(gòu)進行分析，發(fā)現(xiàn)異甘草素有3個酚羥基，異甘草苷和甘草素有2個酚羥基，而甘草苷只有1個酚羥基。根據(jù)化合物所具有的酚羥基的數(shù)目，推測異甘草素的清除自由基活性最強，甘草苷的自由基清除活性最弱。根據(jù)酚羥基鍵長分析，計算結(jié)果列于表1。

根據(jù)分子價鍵理論：分子的鍵長越短，鍵能越大，鍵越不容易斷裂，化學反應(yīng)也越不活潑，抗氧化活性也越低;反之，鍵長越長，抗氧化活性越強。因此通過鍵長大小可以對反應(yīng)活性做初步的判斷。由表1鍵長數(shù)據(jù)分析，對于甘草素來說，C7位酚羥基鍵長大于C4′位酚羥基的鍵長，可知C7位可能是最大活性位點。而甘草苷則是在甘草素的C4′位發(fā)生了糖基化，只有一個C7位的酚羥基。對于異甘草素來說，酚羥基的鍵長由小到大順序為： 9—OH、 ?4'—OH 、7—OH，預(yù)測C7位是可能的最大活性位點，更容易清除自由基。而異甘草苷是在異甘草素的C4′位發(fā)生了糖基化，2個酚羥基的鍵長順序為9—OH鍵長小于7—OH鍵長。由以上分析可知，4種化合物的C7位均為最大活性位點，因此從 ? ? ?C7—OH的鍵長判斷化合物的抗氧化能力由大到小的順序為異甘草素、異甘草苷、甘草素、甘草苷。

2.2 ?NBO（Natural Bond Orbital） 電荷分析

親電基團容易進攻帶負電荷多的原子，親核基團容易進攻帶正電荷多的原子。人體內(nèi)的自由基主要為帶負電荷的氧自由基，所以酚羥基H原子就成為抗氧化反應(yīng)中的活性位點，酚羥基H所帶的正電荷越大，與氧自由基的相互吸引作用越強，化合物的抗氧化活性越強。計算所得的酚羥基H原子的NBO電荷列于表2。

由表2數(shù)據(jù)可看出4種化合物C7位羥基氫上所帶電荷都稍偏大，所以判斷C7位是自由基反應(yīng)的最大可能活性位點。4種化合物C7位羥基H所帶的正電荷由大到小順序為異甘草素、異甘草苷、甘草素、甘草苷，由此推測抗氧化活性由大到小順序為：異甘草素、異甘草苷、甘草素、甘草苷，與前面的鍵長分析結(jié)果一致。

2.3 ?酚羥基的鍵解離焓和電離勢分析

抗氧化劑的抗氧化機理可以根據(jù)化合物酚羥基的鍵解離焓（BDE）和電離勢（IP）來進行分析。抽氫反應(yīng)機理是由最低的BDE值決定，單電子轉(zhuǎn)移反應(yīng)機理是由最低的IP值決定。為此，分別在B3LYP/6-31+G**水平上計算了BDE和IP值。模擬水提和醇提對化合物抗氧化活性的影響[9-11]，表3列出了不同溶劑中的BDE，表4列出了IP值。

為了考察4種黃酮化合物中各個位置酚羥基的抗氧化能力，分別對每個位置的酚羥基失去H后的自由基進行計算，通過比較各個位置的BDE值，來判斷化合物的最大反應(yīng)活性位點。而張鑫[15]等只對失去C7位酚羥基H后的自由基進行了計算。由表3中4種化合物的BDE值可見，對于甘草素來說，不管是在氣相中還是在溶液中，7—OH的BDE值都是最低的，這就表明C7—OH是最大反應(yīng)活性位點;分析異甘草素和異甘草苷的BDE值可知， ? ?C7—OH的BDE值也是最低的;這些結(jié)果與前面的鍵長和電荷分析結(jié)果一致。在氣相和溶劑中，4種化合物的BDE值由大到小順序為：甘草苷、甘草素、異甘草苷、異甘草素，這就表明其抗氧化活性由大到小順序應(yīng)該為：異甘草素、異甘草苷、甘草素、甘草苷。對于甘草素和異甘草素來說，C4′—OH的BDE值略高于C7—OH，說明C4'—OH對化合物的抗氧化活性有貢獻，這可能也是導致成苷后抗氧化活性降低的原因之一。

考查溶劑效應(yīng)對抗氧化機理的影響，由表3可見，隨著溶劑極性的增加，BDE值增大，化合物抗氧化活性減弱，說明溶劑的極性會影響抗氧化劑的抗氧化機理。LITWINIENKO[19]等研究苯酚與自由基在各種溶劑中的反應(yīng)機理發(fā)現(xiàn)，抗氧化機理與溶劑的體積、介電常數(shù)和溶劑化能力有關(guān)。他們指出在介電常數(shù)小的溶劑中，抽氫反應(yīng)機理是主要的。因此，可以得出在非極性溶劑中，4種甘草黃酮類化合物的抗氧化機理以抽氫反應(yīng)機理為主。

由表4中數(shù)據(jù)可知，在氣相中甘草苷的電離勢IP值最大，為205.96 kcal·mol-1，異甘草素的電離勢IP值最小，為162.39 kcal·mol-1。在溶劑中4種化合物的IP值由小到大順序為：異甘草素、異甘草苷、甘草素、甘草苷，由IP值判斷的抗氧化活性由大到小順序亦為：異甘草素、異甘草苷、甘草素、甘草苷，與BDE的分析結(jié)果一致。由此可知，不管是BDE的分析還是IP分析，4種化合物的抗氧化活性由大到小順序為：異甘草素、異甘草苷、甘草素、甘草苷，計算所得的結(jié)論與崔譽蓉[7]等的實驗結(jié)果及張鑫[14]等在氣相中的計算結(jié)果一致。

另外，通過IP值的分析可以看出，在溶劑中，隨著溶劑極性的增加，IP值依次降低，化合物的抗氧化活性增強，說明極性溶劑有利于單電子轉(zhuǎn)移反應(yīng)。IP值與分子的供電子能力有著直接的關(guān)系，而極性溶劑有利于電子的離域，單電子轉(zhuǎn)移反應(yīng)更容易進行[13-14]。這一結(jié)論與LITWINIENKO[19]等得出的溶劑對抗氧化機理影響的結(jié)論一致，他們指出在介電常數(shù)高的溶劑中，例如在水、甲醇和乙醇中，抗氧化劑以單電子轉(zhuǎn)移反應(yīng)為主。

2.4 ?前線分子軌道能級分析

根據(jù)前線分子軌道理論，最高占據(jù)軌道（HOMO）和最低空軌道（LUMO）對分子活性的影響非比尋常。最高占據(jù)軌道HOMO所反映的是分子給出電子的能力，HOMO軌道能量值越大，該軌道上的電子就越活潑，越容易被激發(fā)，反應(yīng)就越容易進行;LUMO軌道能量值越低，則該軌道越容易接受電 ?子[19]。HOMO能級相對越高，HOMO和LUMO之間的能級差（?E）值越小，化學反應(yīng)就越容易發(fā)生，化合物的活性就越強。4種化合物在氣相和溶劑中?E的列于表5。

根據(jù)表5中的數(shù)據(jù)分析，不管是在氣相還是在溶劑中，異甘草素的能級差ΔE最低，說明異甘草素的抗氧化活性最大，然后是異甘草苷、甘草素，能級差最大的是甘草苷。所以根據(jù)前線分子軌道理論可判斷，4種甘草黃酮類成分的抗氧化活性由大到小順序為：異甘草素、異甘草苷、甘草素、甘草苷，這個結(jié)果與BDE和IP的分析結(jié)果一致，與實驗結(jié)果[7]吻合。從能級差分析可得出，溶劑極性增加，能級差減小，化合物的抗氧化活性增強。因此，在極性溶劑中化合物以單電子轉(zhuǎn)移反應(yīng)為主，而在非極性溶劑中，則以抽氫反應(yīng)為主。

3 ?結(jié) 論

本文通過密度泛函理論（DFT），采用B3LYP/6-31+G**方法，對中藥甘草中的4種黃酮類化合物甘草素、異甘草素、甘草苷和異甘草苷的分子結(jié)構(gòu)與抗氧化活性的關(guān)系進行了理論研究。在抽氫反應(yīng)機理中，根據(jù)4種化合物的BDE值分析，4種化合物在氣相中的抗氧化活性由大到小順序為：異甘草素、異甘草苷、甘草素、甘草苷。而本文通過比較C4′—和C7—位酚羥基的BDE值，確定 ? C4'—OH對化合物的抗氧化活性有貢獻，但是 ? ?C7—位酚羥基仍然是最大活性位點。另外，本文還考察了單電子轉(zhuǎn)移反應(yīng)機理，通過IP值分析得出的抗氧化活性由大到小順序也為：異甘草素、異甘草苷、甘草素、甘草苷，進一步說明異甘草素的抗氧化活性是最高的。同時為了模擬真實的實驗環(huán)境，以水、甲醇和乙醇為溶劑，考察溶劑效應(yīng)對化合物抗氧化活性的影響，得出在抽氫反應(yīng)中，隨著溶劑極性的增加，抗氧化活性減弱;而在單電子轉(zhuǎn)移反應(yīng)中，隨著溶劑極性的增加，抗氧化活性增強，與之前的研究結(jié)論是一致的[13]。因此，在極性溶劑中化合物以單電子轉(zhuǎn)移反應(yīng)為主，而在非極性溶劑中，則以抽氫反應(yīng)為主。

參考文獻：

[1]ASL M N， HOSSEINZADEH H. Review of pharmacological effects of Glycyrrhiza sp. and its bioactive compounds[J]. Phytotherapy Research， 2008， 22（6）： 709-724.

[2]ZHANG Q Y， MIN Y， XIE P S， et al. Chemical analysis of the Chinese herbal medicine Gan-Cao （licorice）[J]. Journal Chromatography A， 2009， 1216（11）： 1954-1969.

[3]WANG Z Y， NIXON D W， Licorice and Cancer[J]. Nutrition and Cancer， 2001， 39（1）： 1-11.

[4]TAKASHIMA M， HORIE M， SHICHIRI Y， et al. Assessment of antioxidant capacity for scavenging free radicals in vitro： a rational basis and practical application[J]. Free Radical Biology and Medicine， 2012， 52（7）： 1242-1252.

[5]GOLUBEV A， HANSON A D， GLADYSHEV V N. A Tale of Two Concepts： Harmonizing the Free Radical and Antagonistic Pleiotropy Theories of Aging[J]. Antioxidant Redox Signaling， 2018， 29（10）： 1003-1017.

[6]傅乃武，劉朝陽，張如意，等. 甘草黃酮體抗促癌、抗致突和抗氧化作用的研究[J]. 天然產(chǎn)物研究與開發(fā)，1995，7（4）：29-34.

[7]崔譽榮，陳鵬，劉軍花，等. 4種甘草黃酮類化合物抗氧化構(gòu)效關(guān)系研究[J]. 時珍國醫(yī)國藥，2010，21（12）：3041-3043.

[8]張光輝，張拴，龍旭，等.異甘草素合成工藝優(yōu)化及其清除自由基活性研究[J].當代化工，2017，46（9）：1749-1751.

[9]潘燕. 甘草水溶性總黃酮抗心肌缺血作用的研究[J]. 遼寧中醫(yī)雜志，2004，31（2）：173.

[10]LEE S E， HWANG H J， HA J S， et al. Screening of medicinal plant extracts for antioxidant activity[J]. Life Science， 2003， 73（2）： 167-179.

[11]AMAROWICZ R， PEGG R B， RAHIMI-MOGHADDAM P， et al. Free-radical scavenging capacity and antioxidant activity of selected-plant species from the canadian prairies[J]. Food Chemistry， 2004， 84（4）： 551-562.