2-呋喃甲醛-4-羥基苯甲酰腙及其Cu配合物的晶體結(jié)構(gòu)及與CT-DNA的結(jié)合性能

劉向榮　孫秀超　楊再文　趙順省　楊水蘭　閆森

(西安科技大學(xué)化學(xué)與化工學(xué)院，西安710054)

2-呋喃甲醛-4-羥基苯甲酰腙及其Cu配合物的晶體結(jié)構(gòu)及與CT-DNA的結(jié)合性能

劉向榮＊孫秀超楊再文趙順省楊水蘭閆森

(西安科技大學(xué)化學(xué)與化工學(xué)院，西安710054)

通過2-呋喃甲醛與4-羥基苯甲酰肼縮合得到2-呋喃甲醛-4-羥基苯甲酰腙(H2L)，并以其為配體與Cu配位得到了配合物[Cu(HL)2]·2H2O(1)，采用元素分析和X射線單晶衍射對配體和配合物進行結(jié)構(gòu)表征，結(jié)果表明配體屬于正交晶系，Pna21空間群。配合物屬于單斜晶系，P21/c空間群。利用熱重實驗研究了H2L和1的熱穩(wěn)定性，并計算了它們主要熱分解過程的表觀活化能，發(fā)現(xiàn)H2L和1的熱穩(wěn)定性都較高。通過紫外吸收光譜研究了H2L和1與CT-DNA的相互作用方式，并利用微量熱計測量了其作用過程的熱效應(yīng)，結(jié)果表明H2L和1均以插入方式與CT-DNA結(jié)合，且作用過程放出的熱量1大于H2L，說明配合物與CTDNA的結(jié)合能力強于配體。

酰腙化合物；Cu配合物；晶體結(jié)構(gòu)；CT-DNA；微量熱

0　引言

酰腙化合物分子結(jié)構(gòu)中含有生物活性基團-CONHN=CH-，不僅具有優(yōu)良的抑菌、殺菌、抗腫瘤、除草性等生物藥理活性，而且分子中的O、N均可作為配位原子，具有較強的配位能力和多樣的配位方式[1-2]，尤其是酰腙的過渡金屬配合物，由于融合了有機和無機單元，在生物醫(yī)學(xué)、發(fā)光材料和催化領(lǐng)域應(yīng)用比較廣泛[3-5]。據(jù)文獻報道[6-8]，大部分具有生物活性的化合物易以插入方式與病毒或腫瘤的DNA作用，從而達到抗病毒、抗腫瘤的目的。因此，本文擬將活性基團呋喃引入酰腙分子中實現(xiàn)活性疊加，獲得更具抑、殺菌能力的新物質(zhì)，故以2-呋喃甲醛和4-羥基苯甲酰肼為原料縮合得到2-呋喃甲醛-4-羥基苯甲酰腙，并以其為配體制備Cu配合物。通過X射線單晶衍射和元素分析確定配體及配合物的晶體結(jié)構(gòu)，利用熱重實驗分析它們的熱分解過程，采用紫外吸收光譜探索配體及配合物與CTDNA的作用模式，并借助微量熱法測定其與CTDNA相互作用過程的熱效應(yīng)，以此評價配體和配合物與CT-DNA的結(jié)合能力。

1　實驗部分

1.1主要試劑與儀器

2-呋喃甲醛、4-羥基苯甲酰肼、甲醇均為分析純,CT-DNA(小牛胸腺DNA)購自于美國Sigma公司。

XT4-100B熔點儀；PE-2400-Ⅱ元素分析儀；BRUKER SMART APEXⅡCCD單晶衍射儀；METTLER-TOLEDO TG-DSC1 HT熱重分析儀(N2流量為100.00 mL·min-1)；TU-1900型紫外-可見分光光度計；SETARAM C80微量熱儀。

1.2配體H2L的合成及單晶培養(yǎng)

配體和配合物的合成路線見Scheme 1。

Scheme 1

稱取0.153 2 g(1 mmol)4-羥基苯甲酰肼溶于15.0 mL甲醇中，量取90 μL(1 mmol)2-呋喃甲醛，慢慢滴加到4-羥基苯甲酰肼中，65℃攪拌下水浴回流3 h，冷卻至室溫后過濾，靜置約1周后產(chǎn)生棕黃色塊狀晶體(可供測試用)，產(chǎn)率67.45%。m.p.248～249℃。元素分析按C12H10N2O3計算的理論值(%)：C 62.60，H 4.38，N 12.17；實驗值：C 62.61，H 4.43，N 11.81。

1.3配合物1的合成及單晶培養(yǎng)

稱取0.034 7 g(0.15 mmol)配體H2L溶于10.0 mL甲醇中，待其完全溶解后，慢慢向其滴加0.036 3 g(0.15 mmol)Cu(NO3)2·3H2O的甲醇溶液，65℃攪拌下水浴回流3h，冷卻至室溫后過濾，靜置一段時間后產(chǎn)生墨綠色塊狀晶體(可供測試用)。m.p.＞300℃。元素分析按CuC24H22N4O8計算的理論值(%)：C 51.61，H 3.94，N 10.04；實驗值：C 51.22，H 3.87，N 9.75。

1.4晶體結(jié)構(gòu)測試

選取尺寸分別為0.37 mm×0.29 mm×0.13 mm的H2L的單晶和0.35 mm×0.21 mm×0.12 mm的1的單晶置于Bruker APEX-ⅡCCD單晶衍射儀上，在296(2)K下用經(jīng)石墨單色器單色化的Mo Kα射線(λ=0.071 073 nm)，以φ-ω掃描方式，分別在2.68°～25.10°(H2L)和2.37°～25.10°(1)范圍內(nèi)收集單晶衍射數(shù)據(jù)，然后用SAINT程序進行還原。衍射強度數(shù)據(jù)經(jīng)SADABS程序作經(jīng)驗吸收校正后用SHELXS程序[9]由直接法解出晶體結(jié)構(gòu)，對非氫原子坐標及其各向異性熱參數(shù)用SHELXL程序[10]進行全矩陣最小二乘法修正，與其相連的氫原子由理論加氫法得到。晶體學(xué)數(shù)據(jù)詳見表1。

CCDC：1030242，H2L；1401955，1。

1.5熱重實驗

準確稱取(10.0±0.05)mg H2L或配合物1放入氧化鋁坩堝內(nèi)，升溫程序設(shè)置為室溫至800℃，升溫速率為5.00、10.00和15.00℃·min-1，在氮氣氣氛下進行熱重分析。

1.6紫外吸收光譜

將3 mL 0.01 mol·L-1的Tris-HCl緩沖溶液(pH=7.90)加入?yún)⒈缺壬笾校驑悠繁壬笾屑尤氲润w積現(xiàn)配制的1×10-4mol·L-1配體溶液，用微量進樣器分別往參比比色皿和樣品比色皿中加入相同體積(50 μL)100 mg·L-1的CT-DNA溶液，使配體與CT-DNA的濃度比值逐漸減小，連續(xù)加入5次，每次間隔約需5 min混合均勻，掃描波長范圍為250～500 nm；配合物采用同樣的測試方法。

1.7微量熱實驗

準確稱取2.0 mg的配體放入樣品池下層，分別向參比池和樣品池下層各加0.01 mol·L-1的Tris-HCl緩沖溶液(pH=7.90)2 mL，再分別向參比池和樣品池的上層各加1 mL 100 mg·L-1的CT-DNA溶液，25℃下，用微量熱儀測量配體與CT-DNA作用的熱量變化；用同樣的方法測量配合物與CT-DNA作用過程的熱量變化。

表1　H2L和1的晶體學(xué)數(shù)據(jù)Table 1Crystallographic data for H2L and 1

2　結(jié)果與討論

2.1晶體結(jié)構(gòu)

從表1可以看出，H2L的晶體屬于正交晶系，空間群為Pna21，其分子結(jié)構(gòu)見圖1，主要鍵長及鍵角見表2。從表2可以看出，所列出的鍵長和鍵角的數(shù)據(jù)都在正常范圍之內(nèi)，H2L中苯環(huán)與呋喃環(huán)的二面角為32.5°，表明兩個環(huán)不在同一平面上。N2-C8鍵長為0.128 2(3)nm，小于N1-C7(0.135 0(3)nm)的鍵長，說明N2和C8之間形成了雙鍵[11]。分子中O2-C7鍵長為0.123 6(3)nm，屬于典型的C=O雙鍵，證明配體是以酮羰基形式存在的。圖2為配體的氫鍵圖，從圖2可以看出配體分子間通過氫鍵N-H…O (N1-H1…O2i氫鍵：H1…O2i0.213 nm，N1…O2i0.291 4(2)nm，N1-H1…O2i151.1°)和氫鍵O-H…N (O1-H1A…N2ii氫鍵：H1A…N2ii0.216 nm，O1··· N2ii0.296 9(3)nm，O1-H1A…N2ii169.5°)的弱作用力拓展為如圖3所示波浪形堆積結(jié)構(gòu)。氫鍵鍵長和鍵角見表3。

圖1　H2L的分子結(jié)構(gòu)Fig.1　Molecular structure of H2L

配合物1屬于單斜晶系，P21/c空間群，其分子結(jié)構(gòu)見圖4，主要鍵長及鍵角見表2，從表2可以看出所列出的鍵長和鍵角的數(shù)據(jù)都在正常范圍之內(nèi)。1的晶體單胞結(jié)構(gòu)中包含1個Cu2+、2個配體負離子，并含有結(jié)晶水分子。1屬于雙齒配位，每個配體負離子上的亞胺N和羰基O都參與配位，形成四配位的單核結(jié)構(gòu)，Cu2+所在環(huán)的內(nèi)角求和是539.9°，說明Cu2+與配位原子組成了閉合的五元環(huán)。O2-Cu1-O2i和N2-Cu1-N2i鍵角均為180.0°，表明Cu2+周圍的2個五元環(huán)幾乎處于同一平面。配位后分子中的O2-C7鍵長變長，N1-C7鍵長變短，且配合物中N1-C7鍵長為0.133 2(9)nm，比吡啶環(huán)中N-C雙鍵(0.135 2 nm)短，說明配體是以烯醇形式配位的。圍繞Cu2+的Cu1-O2鍵長為0.191 2(5)nm，Cu1-N2鍵長為0.192 5(6)nm，與類似結(jié)構(gòu)的銅配合物的鍵長基本一致[12-13]。

表2H2L和1的主要鍵長(nm)及鍵角(°)Table 2Selected bond lengths(nm)and angles(°)for H2L and 1

表3　H2L的氫鍵鍵長(nm)及鍵角(°)Table 3Hydrogen bond distances(nm)and angles(°)of H2L

圖2　H2L的氫鍵圖Fig.2　View of the hydrogen bonds of H2L

圖3　H2L的晶胞堆積圖Fig.3　Crystal packing of H2L

如圖5所示，配合物1中苯環(huán)上羥基的O1、H1A和結(jié)晶水中的O4形成氫鍵O1-H1A…O4ii(0.272 4(9)nm，151.9°)，結(jié)晶水中的O4、H4B、H4C和N1、O1形成氫鍵O4-H4C…N1iii(0.301 3(10)nm，150(7)°)和O4-H4B…O1iv(0.287 1(10)nm，157(7)°)，單體分子通過以上3種弱氫鍵作用相互連接成如圖6所示規(guī)則交叉排列的三維超分子結(jié)構(gòu)。氫鍵參數(shù)詳見表4。

表4　配合物1的氫鍵鍵長(nm)及鍵角(°)Table 4Hydrogen bond distances(nm)and angles(°)of complex 1

圖4　配合物1的分子結(jié)構(gòu)Fig.4Molecular structure of complex 1

圖5　配合物1的氫鍵圖Fig.5Hydrogen bond diagram of complex 1

圖6　配合物1的晶胞堆積圖Fig.6Crystal packing diagram of complex 1

圖7　H2L在5.00℃·min-1下的TG-DTG曲線Fig.7TG-DTG curves of H2L at the heating rate of 5.00℃·min-1

2.2熱重分析

2.2.1配體H2L熱重分析

配體H2L在5.00℃·min-1升溫速率下的TGDTG曲線如圖7所示，在升溫速率分別為5.00、10.00和15.00℃·min-1時的DTG曲線見圖8。可以看出配體在3種升溫速率下的熱分解過程相似，都只有一個階段，分解峰分別出現(xiàn)在291.50、303.55和310.89℃，表明配體在291.50℃以下可以穩(wěn)定存在。升溫速率為5.00℃·min-1時實驗失重率為70.67%，計算值70.90%(與OH-C6H4-CO-NH-N=C-的百分含量相同)，因此可推測熱分解過程中C8-C9鍵斷裂。

采用Kissinger和Ozawa公式[14]計算配體熱分解過程表觀活化能，計算公式分別見(1)和(2)式。

圖8　H2L在3種升溫速率下的DTG曲線Fig.8DTG curves of H2L at three heating rates

其中，Tp為分解峰的溫度，A為指前因子，Ea為表觀活化能，R為氣體摩爾常數(shù)，β為升溫速率，G(α)為積分機理函數(shù)。利用公式(1)和(2)可求得表觀活化能Ea及指前因子A，計算結(jié)果列于表5。

2.2.2配合物1熱重分析

配合物1的熱重分析曲線如圖9～10所示，配合物在3種升溫速率下的DTG曲線很相似。從圖9可見配合物的熱分解過程主要分為2個階段：第1階段分解的實驗失重率為7.97%，該階段是配合物中2分子結(jié)晶水的失去，第2階段對應(yīng)的實驗失重率為60.16%，最大失重率對應(yīng)的熱分解溫度配合物與配體接近，均大于290℃，說明配合物和配體的熱穩(wěn)定性都較高。利用公式(1)和(2)計算的配合物熱分解過程的表觀活化能也列于表5。從表5可以看出，配體和配合物的主要熱分解階段的表觀活化能分別為：145.8 kJ·mol-1和264.1 kJ·mol-1。

圖9　配合物1在5.00℃·min-1下的TG-DTG曲線Fig.9TG-DTG curves of complex 1 at the heating rateof 5.00℃·min-1

圖10　配合物1在3種升溫速率下的DTG曲線Fig.10DTG curves of complex 1 at three heating rates

表5H2L和1在3種升溫速率下的熱分解動力學(xué)參數(shù)Table 5Kinetic parameters of thermal decomposition for H2L and 1

2.3CT-DNA結(jié)合性能

2.3.1H2L和1與CT-DNA作用的紫外吸收光譜

分析

配體和配合物與CT-DNA相互作用的紫外吸收光譜圖分別見圖11和圖12，從圖11可以看出，隨CT-DNA濃度的增加，配體體系的吸光度值逐漸減小(0.42→0.38)，在326 nm處發(fā)生了減色效應(yīng)和微弱的紅移現(xiàn)象，可認為配體與CT-DNA發(fā)生插入作用，這是由于CT-DNA堿基對π電子軌道與配體中π*電子空軌道發(fā)生耦合，從而導(dǎo)致π-π*躍遷能減小[15-19]，產(chǎn)生紅移現(xiàn)象，此外，配體插入CT-DNA后，長鏈DNA的“屏蔽”使π平面接受的光通量減少，也會導(dǎo)致減色效應(yīng)[20-21]。圖12中配合物體系的吸光度也隨CT-DNA濃度的增加而逐漸減小(0.22→0.19，0.15→0.13)，分別在297 nm和355 nm處發(fā)生了減色效應(yīng)和微弱的紅移現(xiàn)象，表明配合物也是以插入作用與CT-DNA結(jié)合，配合物與CTDNA作用時，具有芳香性的配體會部分嵌入到DNA的相鄰堿基對之間，這一插入作用是藥物分子與DNA結(jié)合的一個重要模式，而且金屬插入劑是這種結(jié)合模式極有用的探針[22]。

圖11　H2L與CT-DNA作用的紫外光譜圖Fig.11UV-Vis spectra of H2L when interacting with CT-DNA

圖12　配合物1與CT-DNA作用的紫外光譜圖Fig.12UV-Vis spectra of 1 when interacting with CT-DNA

配體和配合物與CT-DNA的結(jié)合常數(shù)Kb可以通過收集在指定波長下的滴定曲線的吸光度并利用方程(3)[23]求得，其中εa，εf和εb分別表示任意CTDNA濃度下Tris-HCl緩沖溶液的摩爾消光系數(shù)、配體或配合物的摩爾消光系數(shù)和配體或配合物與CTDNA完全鍵合后的摩爾消光系數(shù)[24]。以CDNA/(εa-εf)對CCT-DNA作圖，配體和配合物的CDNA/(εa-εf)-CDNA關(guān)系圖分別見圖13和圖14，斜率與截距的比值即為樣品與CT-DNA的結(jié)合常數(shù)Kb。通過計算，配體和配合物的結(jié)合常數(shù)分別為：1.98×106L·mol-1和3.68×107L·mol-1，表明配體和配合物與CT-DNA都有較強的插入作用[23]，并且配合物與CT-DNA的結(jié)合能力大于配體。

圖13　H2L的CDNA/(εa-εf)-CDNA關(guān)系圖Fig.13Plot of CDNA/(εa-εf)vs CDNAfor H2L

圖14　配合物1的CDNA/(εa-εf)-CDNA關(guān)系圖Fig.14Plot of CDNA/(εa-εf)vs CDNAfor complex 1

2.3.2H2L和1與CT-DNA作用的微量熱分析

圖15和圖16分別是配體及配合物與CT-DNA相互作用的微量熱實驗熱譜圖。從圖15可以看出配體與CT-DNA混合后迅速放熱，在11.20 min時出現(xiàn)一個最大放熱峰，35.10 min時結(jié)束，通過計算峰面積得出配體與CT-DNA作用的焓變值為-4.56 kJ·mol-1。從圖16可以看出配合物與CT-DNA混合后發(fā)生了與配體類似的放熱過程，最大放熱峰出現(xiàn)在12.36 min，在49.37 min時結(jié)束，配合物與CTDNA作用的焓變值為-22.53 kJ·mol-1，從焓變值可知配合物與CT-DNA有更強的作用。由于配體和配合物與CT-DNA作用焓變的絕對值都小于30 kJ· mol-1，因此可認為配體和配合物與CT-DNA沒有發(fā)生化學(xué)反應(yīng)，不存在化學(xué)鍵[25-26]，這也印證了2.3.1節(jié)中的結(jié)果，即配體和配合物與CT-DNA均以插入方式相互作用。

圖15　H2L與CT-DNA作用熱譜圖Fig.15Thermogenic curve of H2L when interacting with CT-DNA

圖16　配合物1與CT-DNA作用熱譜圖Fig.16Thermogenic curve of 1 when interacting with CT-DNA

3　結(jié)論

以2-呋喃甲醛和4-羥基苯甲酰肼為原料制備了配體2-呋喃甲醛-4-羥基苯甲酰腙(H2L)，并進一步制得其Cu配合物[Cu(HL)2]·2H2O(1)，利用元素分析和X射線單晶衍射確定了其結(jié)構(gòu)。熱重實驗結(jié)果說明H2L和1的熱穩(wěn)定性都較高。通過紫外吸收光譜法和微量熱法研究了H2L和1與CT-DNA的相互作用模式和強弱，結(jié)果表明H2L和1與CT-DNA皆以插入模式結(jié)合，作用過程放熱，1與CT-DNA的結(jié)合能力強于H2L。

[1]YANG Rui(楊銳),HE Shui-Yang(何水樣),WU Wang-Ting (武望婷),et al.Acta Chim.Sinica(化學(xué)學(xué)報),2004,62(20): 2040-2044

[2]Monfort M,Resino J,Ribas J,et al.Inorg.Chem.,2000,3(9): 2572-2576

[3]REN Hai-Xian(任海仙),TANG Jing(唐靜),WEI Tai-Bao(魏太保),et al.Chinese J.Inorg.Chem.(無機化學(xué)學(xué)報),2007, 23(11):1907-1911

[4]Schuetz S A,Silvernail C M,Incarvito C D.Inorg.Chem., 2004,4(3):6203-6207

[5]Adsule S,Barve V,Ahmed F,et al.J.Med.Chem.,2006,49 (24):7242-7246

[6]Liu S,Liu F,Yu X,et al.Bioorg.Med.Chem.,2006,14: 1425-1430

[7]LI Ling(李玲),DONG Tong-Yi(董同義),LI Xiu-Lu(李修錄), et al.Acta Pharm.Sin.(藥學(xué)學(xué)報),1994,29:128-131

[8]Cao S L,Feng Y P,Jiang Y Y,et al.Bioorg.Med.Chem. Lett.,2005,15:1915-1917

[9]Sheldrick G M.SHELXS-97,Program for the Solution of Crystal Structures,University of G?ttingen,Germany,1997.

[10]Sheldrick G M.SHELXL-97,Program for the Refinement of Crystal Structures,University of G?ttingen,Germany,1997.

[11]LONG De-Qing(龍德清),CHEN Sheng-Sheng(陳勝勝), CHEN Hui(陳輝).Chinese J.Synth.Chem.(合成化學(xué)), 2006,14(1):69-71

[12]WU Qiong-Jie(吳瓊潔),CHEN Xiao-Hua(陳小華),CAI Bi-Qiong(蔡碧瓊),et al.Chinese J.Inorg.Chem.(無機化學(xué)學(xué)報),2012,28(1):201-206

[13]ZHENG Chang-Zheng(鄭長征),WANG Yan-Yan(王燕燕), WANG Liang(王亮).Chinese J.Inorg.Chem.(無機化學(xué)學(xué)報),2014,30(7):1474-1480

[14]HU Rong-Zu(胡榮祖),GAO Sheng-Li(高勝利),ZHAO Feng -Qi(趙鳳起),et al.Thermal Analysis Kinetics.2nd(熱分析動力學(xué).2版).Beijing:Science Press,2008:79-120

[15]GAO Ting(高婷),ZHANG Wan-Ju(張萬舉),WANG Fang (王芳),et al.Guangzhou Chem.Ind.(廣州化工),2010,38(5): 23-25

[16]Wolfe A,Shimer G H,Meehan T.Biochem.,1987,26:6392-6396

[17]Pyle A M,Rehmann J P,Meshoyrer R,et al.J.Am.Chem. Soc.,1989,111:3051-3058

[18]Raman N,Sobha S.Inorg.Chem.Commun.,2012,17:120-123

[19]Eberhard J,Stoll I,Brockhinke R,et al.CrystEngComm, 2013,15:4225-4248

[20]Rajput C,Rutkaite R,Swanson L,et al.Chem.Eur.J., 2006,12:4611-4619

[21]Mabeswari P U,Palaniandavar M.J.Inorg.Biochem.,2004, 98:219-230

[22]Carter M T,Rodriguez M,Bard A J.J.Inorg.Biochem., 2008,102:330-341

[23]ZHU Li(朱莉),PENG Bin(彭斌),LING You(凌友),et al. Acta Chim.Sinica(化學(xué)學(xué)報),2008,66(24):2705-2711

[24]Shi S,Geng X T,Zhao J,et al.Biochimie,2010,92:370-377

[25]Solimani R.Biochim.Biophys.Acta,1997,1336(2):281-294

[26]Peter J,Ingermar W.J.Biochem.Biopphys.Methods,1997, 35(103):103-114

2-Furancarbaldehyde-4-hydroxy-benzoylhydrazone and Its CuComplex: Crystal Structures and Binding Ability with CT-DNA

LIU Xiang-Rong＊SUN Xiu-ChaoYANG Zai-Wen ZHAO Shun-ShengYANG Shui-LanYAN Sen
(College of Chemistry and Chemical Engineering,Xi′an University of Science and Technology,Xi′an 710054,China)

2-furancarbaldehyde-4-hydroxy-benzoylhydrazone(H2L)was synthesized from 4-hydroxybenzoylhydrazine and 2-furancarbaldehyde,and its Cucoordination polymer[Cu(HL)2]·2H2O(1)also was prepared.H2L and 1 were characterized by elemental analysis and single crystal X-ray diffraction analysis.The crystal structural analyses show that the crystal of H2L belongs to orthorhombic system space group Pna21,and the crystal of 1 crystallizes in monoclinic system space group P21/c.Thermal gravity analyses were used to investigate the thermal stabilities of H2L and 1,and their apparent activation energy of the decompositions were also calculated.It was found that both of the complex and the ligand possessed higher thermal stability.The binding modes of H2L and 1 with CT-DNA were studied by UV-Vis absorption,and their interaction enthalpies were measured by microcalorimetry.The results showed that the interaction of H2L and 1 with CT-DNA belonged to partial intercalation mode and 1 presented stronger interaction with CT-DNA than H2L.CCDC:1030242,H2L;1401955,1.

hydrazone compounds;Cucomplex;crystal structure;CT-DNA;microcalorimetry

O614.121

A

1001-4861(2016)02-0250-09

10.11862/CJIC.2016.037

2015-08-19。收修改稿日期：2015-11-22。

國家自然科學(xué)基金項目(No.21073139，21301139，21103135)和陜西省自然科學(xué)基礎(chǔ)研究計劃項目(No.2015JQ2043)資助。

＊通信聯(lián)系人。E-mail：liuxiangrongxk@163.com；會員登記號：S06N9491M1008。