陳海龍,田耀旗,李丹,金征宇*
1(江南大學(xué),食品科學(xué)與技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室,江蘇 無錫,214122)2(江南大學(xué) 食品學(xué)院,江蘇 無錫,214122) 3(紹興市食品藥品檢驗(yàn)研究院,浙江 紹興,312000)
脂肪酶交聯(lián)聚集體的制備及其催化合成月桂酸淀粉酯的研究
陳海龍1,2,田耀旗1,2,李丹1,3,金征宇1,2*
1(江南大學(xué),食品科學(xué)與技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室,江蘇 無錫,214122)2(江南大學(xué) 食品學(xué)院,江蘇 無錫,214122) 3(紹興市食品藥品檢驗(yàn)研究院,浙江 紹興,312000)
采用交聯(lián)酶聚集體(cross-linked enzyme aggregates, CLEAs)技術(shù)制備了脂肪酶CLEAs,優(yōu)化了制備條件:以硫酸銨為沉淀劑,以14 mmol/L戊二醛為交聯(lián)劑,交聯(lián)時(shí)間1.5 h,所得脂肪酶CLEAs的相對(duì)酶活力達(dá)到41.89%;與游離酶相比,脂肪酶CLEAs的熱穩(wěn)定性顯著增強(qiáng),重復(fù)使用性顯著提高,重復(fù)使用10次后相對(duì)酶活力達(dá)到50.26%。以該脂肪酶CLEAs作為催化劑,以玉米淀粉和月桂酸為原料,進(jìn)行了酯化反應(yīng),得到月桂酸淀粉酯的最大取代度為0.151。此外,利用紅外光譜(Fourier transform infrared spectoscopy,F(xiàn)T-IR)對(duì)月桂酸淀粉酯進(jìn)行了分子表征,發(fā)現(xiàn)在1 713 cm-1處出現(xiàn)了羰基CO的特征吸收峰,說明月桂酸與淀粉在脂肪酶CLEAs的催化作用下合成了月桂酸淀粉酯。同時(shí),熱重分析(TGA)表明,相比于原淀粉,酯化產(chǎn)物的熱穩(wěn)定有較低趨勢(shì)。
脂肪酶;交聯(lián)酶聚集體;酯化;月桂酸淀粉酯
淀粉是由多葡萄糖單元組成的天然高分子碳水化合物,具有較強(qiáng)的親水性,限制了其工業(yè)化應(yīng)用,通常需要通過酯化、醚化、氧化等方法進(jìn)行改性。酯化淀粉是通過淀粉結(jié)構(gòu)中的羥基與酸發(fā)生酯化反應(yīng)后得到的一類變性淀粉。相比于傳統(tǒng)的化學(xué)合成法,脂肪酶法合成淀粉酯具有條件溫和、綠色無毒、產(chǎn)物易于分離純化等優(yōu)點(diǎn),在食品、醫(yī)藥、化工等領(lǐng)域有著廣泛的潛在應(yīng)用。作為一種新型的無載體固定化酶技術(shù),CLEAs[1]技術(shù)無需其他載體,對(duì)酶的純度要求低,操作簡(jiǎn)單,成本較低。目前,關(guān)于脂肪酶法合成淀粉酯的研究主要集中于不同底物和反應(yīng)體系[2-3]、制備工藝的優(yōu)化[4]以及反應(yīng)動(dòng)力學(xué)分析[5]等方面,淀粉酯通常是由商業(yè)化固定脂肪酶催化合成的,而尚未有利用脂肪酶CLEAs催化合成淀粉酯的報(bào)道。
本研究對(duì)脂肪酶CLEAs的制備工藝進(jìn)行了優(yōu)化,并研究了其熱穩(wěn)定性和重復(fù)使用性;其次在脂肪酶CLEAs的催化作用下,以玉米淀粉和月桂酸為原料,考察了不同條件對(duì)產(chǎn)物取代度的影響,同時(shí)對(duì)淀粉酯產(chǎn)物進(jìn)行了分子表征,以期為淀粉酯的酶法高效合成提供新的途徑。
1.1 試劑與儀器
普通玉米淀粉,杭州普羅星淀粉有限公司;月桂酸、無水乙醇、戊二醛、乙二醇二甲醚、叔丁醇、丙酮、異丙醇、DMSO、(NH4)2SO4等試劑均為分析純,上海國(guó)藥集團(tuán);棕櫚酸對(duì)硝基苯酯(p-Nitrophenyl palmitate,p-NPP),上海阿拉丁有限公司;豬胰腺脂肪酶,Type II,美國(guó)Sigma-Aldrich公司。
傅立葉變換紅外光譜儀(FT-IR,5DXC),美國(guó)Nicolet公司;熱重分析儀(TGA,TGA/SDTA851e),瑞士Mettler-Toledo公司;高速冷凍離心機(jī)(BIOFUGEPRIMOR),美國(guó)Thermo Electron有限公司;雙光束紫外可見分光光度計(jì)(TU-1901),北京普析通用儀器有限責(zé)任公司。
1.2 實(shí)驗(yàn)方法
1.2.1 交聯(lián)酶聚集體的制備
參考CAO[6]等人的方法并略做改進(jìn)。取9 mL的沉淀劑于離心管中,逐滴加入1 mL的脂肪酶液(pH 7.5,0.05 mol/L),25 ℃下水浴振蕩沉淀1.0 h,然后加入戊二醛溶液,在25 ℃下繼續(xù)振蕩交聯(lián),將混合物在4 ℃、10 000 r/min條件下離心5 min,去掉上清液,用磷酸緩沖液清洗3次以除去殘留的戊二醛及未交聯(lián)的脂肪酶。將得到的沉淀重懸于緩沖液或冷凍干燥后,于4 ℃冰箱中冷藏備用。
(1)
1.2.2 脂肪酶水解活性的測(cè)定
參考VORDERWüLBECKE[7]等人的方法并略做改進(jìn)。將90 mgp-NPP溶解于30 mL異丙醇中,與450 mL磷酸緩沖液(pH 7.5,0.05 mol/L)混合均勻,然后取1 mL上述混合溶液,預(yù)熱5 min后加入1 mL稀釋酶液,37 ℃下恒溫水浴反應(yīng)5 min,立即加入2 mL乙醇以終止反應(yīng),于410 nm處測(cè)定吸光度,平行測(cè)定3次。同時(shí)以1 mL磷酸緩沖液作為空白對(duì)照。在一定條件下,每分鐘分解底物p-NPP產(chǎn)生1 μmol對(duì)硝基苯酚所需的酶量定義為1個(gè)酶活力單位。
1.2.3 熱穩(wěn)定性與重復(fù)使用性的測(cè)定
取一定量的游離酶和固定化酶分別置于30、40、50、60、70 ℃的水浴鍋中溫育2 h,參照1.3.2的方法測(cè)定殘余酶活,并定義初始酶活為相對(duì)酶活100%。
取一定量的固定化酶加入底物溶液中,于37 ℃條件下反應(yīng)5 min,過濾出酶并用磷酸緩沖液清洗,然后重新加入新的底物溶液,參照1.2.2的方法測(cè)定酶活。重復(fù)上述操作10次,并測(cè)定每次反應(yīng)結(jié)束之后的殘余酶活。將第一次反應(yīng)的酶活力定義為100%。
1.2.4 月桂酸淀粉酯的合成
參考ZARSKI[8]等人的方法并略做改進(jìn)。普通玉米淀粉于105 ℃烘箱中過夜干燥后,在二甲基亞砜中加入5%(w/v)的淀粉,90 ℃條件下充分?jǐn)嚢? h使淀粉溶解,然后將溶液冷卻至反應(yīng)溫度,加入月桂酸、脂肪酶交聯(lián)聚集體20%(w/w),水浴振蕩反應(yīng)4 h。反應(yīng)結(jié)束后,離心除去交聯(lián)酶聚集體,在上清液中攪拌加入無水乙醇,將淀粉酯沉淀析出,并用乙醇洗滌3~4次,再用丙酮洗滌分散,沉淀在60 ℃烘箱中過夜干燥,研磨粉碎獲得淀粉酯。
1.2.5 月桂酸淀粉酯取代度的測(cè)定
參考MILADINOV[9]等人的方法并略做改進(jìn)。準(zhǔn)確稱取0.5 g樣品置于150 mL錐形瓶中,加入5 mL乙醇和25 mL去離子水,以及25 mL 0.1 mol/L的NaOH溶液,在室溫下振蕩60 min,然后加入2~3滴酚酞指示劑,用0.1 mol/L的標(biāo)準(zhǔn)HCl溶液滴定至粉紅色剛好消失,記錄消耗的HCl體積,平行測(cè)定3次。同時(shí)以原淀粉作為空白試驗(yàn)。計(jì)算公式如下:
(2)
式中:DS,取代度,定義為每個(gè)D-吡喃葡萄糖殘基中的羥基被取代的平均數(shù)量;162,葡萄糖酐單元的分子量;183,月桂酸?;姆肿恿?;c,標(biāo)準(zhǔn)HCl的摩爾濃度,mol/L;m,干樣品的質(zhì)量,g;V0,滴定空白液用去的標(biāo)準(zhǔn)鹽酸的體積,mL;V1,滴定樣品消耗的標(biāo)準(zhǔn)鹽酸的體積,mL。
1.2.6 傅里葉轉(zhuǎn)換紅外光譜(FT-IR)
采用KBr壓片法制樣,分別將原淀粉與月桂酸淀粉酯粉末與干燥的KBr按合適的比例混勻,在紅外燈下使用瑪瑙研缽將其研磨成細(xì)粉,轉(zhuǎn)移至壓片機(jī)模具中,使用壓片機(jī)將樣品壓制成均勻薄片。以空氣作為空白,將薄片放入傅里葉變換紅外光譜儀中掃描。紅外光譜儀掃描波數(shù)范圍為4 000~400 cm-1,掃描次數(shù)32次,分辨率為4 cm-1。
1.2.7 熱重分析(TGA)
分別稱取3~4 mg的原淀粉與月桂酸淀粉酯于坩堝中,置于熱重分析儀系統(tǒng)中。設(shè)置氮?dú)饬髁繛?0 mL/min,升溫范圍為50~600 ℃,升溫速率為10 ℃/min。
1.2.8 數(shù)據(jù)分析
所有數(shù)據(jù)均為至少3次實(shí)驗(yàn)的平均值,誤差均為對(duì)應(yīng)的標(biāo)準(zhǔn)偏差;利用Origin 8.5 軟件對(duì)數(shù)據(jù)進(jìn)行圖表繪制和統(tǒng)計(jì)分析。
2.1 固定化酶的制備與性質(zhì)表征
2.1.1 固定化酶的制備
CLEAs的制備過程主要包括:先利用蛋白沉淀劑將溶解狀態(tài)的酶進(jìn)行物理沉淀形成聚集體,再利用交聯(lián)劑與酶發(fā)生交聯(lián)反應(yīng)生成不溶的CLEAs。其中,沉淀劑種類、交聯(lián)時(shí)間、戊二醛濃度和酶濃度對(duì)CLEAs的最終酶活影響較大。
沉淀過程是CLEAs制備步驟中關(guān)鍵的第一步。常用的蛋白沉淀劑主要包括中性鹽、有機(jī)溶劑和非離子型聚合物3類。沉淀劑對(duì)不同過程中酶活力的影響如表1所示。首先,加入沉淀劑使酶發(fā)生沉淀聚集,再將沉淀后的酶聚集體重新溶解回游離酶狀態(tài),考察此時(shí)復(fù)溶聚集體的相對(duì)酶活力。異丙醇做沉淀劑時(shí)復(fù)溶酶活高達(dá)101.37%,甚至出現(xiàn)了超酶活現(xiàn)象[10]。然后,又考察了沉淀劑對(duì)CLEAs最終酶活力的影響。以硫酸銨做沉淀劑時(shí),所得到的CLEAs的最終活性卻是最高的,達(dá)到41.89%。沉淀步驟的最佳沉淀劑并不是制備CLEAs所需的最佳沉淀劑,其原因可能是:一、沉淀步驟形成的酶聚集體可能處于一種不利于反應(yīng)的狀態(tài),但是當(dāng)它們復(fù)溶于緩沖液后卻會(huì)表現(xiàn)出正常的活性;二、當(dāng)酶聚集體被交聯(lián)后,會(huì)以這種不利于反應(yīng)的狀態(tài)保留下來,并且表現(xiàn)出很低的活性[11]。因此,選取硫酸銨作為沉淀劑。
表1 不同沉淀劑對(duì)復(fù)溶聚集體及CLEAs相對(duì)酶活力的影響
交聯(lián)時(shí)間是CLEAs制備過程中的重要影響因素之一。由圖1(a)可知,制備脂肪酶CLEAs的最佳交聯(lián)時(shí)間為1.5 h。當(dāng)交聯(lián)時(shí)間小于1.5 h時(shí),時(shí)間的延長(zhǎng)有利于酶分子的固定,生成的CLEAs酶活也逐漸上升;超過1.5 h后,由于CLEAs的過度交聯(lián),導(dǎo)致其相對(duì)酶活力開始下降。
戊二醛是制備CLEAs的過程中最常用的交聯(lián)劑。戊二醛對(duì)酶的交聯(lián)主要是通過其兩端的醛基和酶分子表面的α-賴氨酸發(fā)生共價(jià)結(jié)合,得到不溶的交聯(lián)酶聚集體。戊二醛濃度過低時(shí),交聯(lián)不充分,得到的交聯(lián)酶聚集體量較少,容易造成酶固定量的損失;但是當(dāng)濃度過高時(shí),酶分子間及分子內(nèi)易發(fā)生交聯(lián),戊二醛甚至?xí)c酶活性中心結(jié)合,造成酶活性的下降。如圖1(b)所示,戊二醛濃度為14 mmol/L時(shí)獲得的固定化酶相對(duì)酶活力最高。
酶濃度主要影響著CLEAs的顆粒大小。如圖1(c)所示,隨著酶濃度的增大,CLEAs的相對(duì)酶活力也逐漸升高;當(dāng)酶濃度高于30 mg/mL時(shí),CLEAs的活力開始下降,這可能是因?yàn)樾纬傻拿妇奂w顆粒逐漸增大,造成擴(kuò)散限制,阻礙了底物分子與酶活性中心的接觸。因此,選用30 mg/ mL的酶液來制備CLEAs。
圖1 交聯(lián)時(shí)間(a)、戊二醛濃度(b)、酶濃度(c)對(duì)CLEAs相對(duì)酶活力的影響Fig.1 Effects of the cross-linking time(a), concentration of glutaraldhyde(b) and enzyme concentration (c) on relative activity of CLEAs
2.1.2 固定化酶的熱穩(wěn)定性及重復(fù)使用性
2.1.2.1 固定化酶的熱穩(wěn)定性
實(shí)驗(yàn)考察了(a)游離脂肪酶和(b)CLEAs在不同溫度下的熱穩(wěn)定性,結(jié)果如圖2所示。隨著溫度的升高,游離酶和CLEAs的活力均呈現(xiàn)下降趨勢(shì),但是CLEAs的熱穩(wěn)定性遠(yuǎn)高于游離酶。在70 ℃,保溫2 h的條件下,游離酶只保留了21.40%的活性,但是CLEAs仍保留54.88%的活性。這是由于酶分子之間發(fā)生交聯(lián)之后,分子之間連接更緊密,酶失活需要克服更多的非共價(jià)及共價(jià)作用力,同時(shí)酶分子內(nèi)部基團(tuán)的熱振動(dòng)減小,酶分子發(fā)生熱伸展變性的可能性變小,因此CLEAs的熱穩(wěn)定性增強(qiáng),抗熱變性能提高。
圖2 游離酶(a)和CLEAs (b)的熱穩(wěn)定曲線Fig.2 Thermal stabilities of free enzyme (a) and CLEAs (b)
2.1.2.2 固定化酶的重復(fù)使用性
CLEAs的重復(fù)使用性如圖3所示。隨著使用次數(shù)的增多,CLEAs的酶活力逐漸降低,使用前5次其酶活力降低較快,但是之后,其酶活力損失很小。這可能是因?yàn)椋篊LEAs在制備的過程中,部分酶分子間并未發(fā)生交聯(lián),而只是吸附在其表面,在重復(fù)回收的操作過程中可能會(huì)發(fā)生酶的泄露及部分活性損失;使用10次后,固定化酶不再泄露并且保持穩(wěn)定,CLEAs仍保持50.26%的相對(duì)酶活力。
圖3 CLEAs的重復(fù)使用性Fig.3 Reusability of CLEAs
2.2 淀粉酯的制備與結(jié)構(gòu)表征
2.2.1 月桂酸淀粉酯的制備
取代度,是指淀粉的每個(gè)葡萄糖單元上的活性羥基被取代的物質(zhì)的量。由于淀粉分子中每個(gè)脫水葡萄糖單元上只有3個(gè)羥基能被取代,所以取代度的理論值最大是3。表2是在不同反應(yīng)溫度和月桂酸/ AGU摩爾比的條件下,得到的月桂酸淀粉酯的取代度。隨著溫度的升高,產(chǎn)物取代度不斷減小,其主要原因是因?yàn)闇囟雀哂?0 ℃時(shí),脂肪酶的穩(wěn)定性降低,酶發(fā)生了熱失活;隨著月桂酸/AGU摩爾比的不斷增大,月桂酸相對(duì)于淀粉是過量的,月桂酸的濃度變大,反應(yīng)平衡也就向著正方向移動(dòng),因此產(chǎn)物的取代度不斷增大。本實(shí)驗(yàn)結(jié)果表明,在溫度為60 ℃、底物摩爾比為3∶1時(shí),月桂酸淀粉酯的取代度最高為0.151。
表2 反應(yīng)溫度及月桂酸/AGU的摩爾比對(duì)月桂酸淀粉酯取代度的影響
2.2.2 月桂酸淀粉酯的結(jié)構(gòu)表征
2.2.2.1 紅外光譜分析
圖4為(a)原淀粉和(b)月桂酸淀粉酯的紅外光譜圖。圖4(a)玉米原淀粉在3 419 cm-1處有一寬峰,為淀粉葡萄糖殘基(AGU)上O—H的伸縮振動(dòng);2 930 cm-1出現(xiàn)的吸收峰為C—H的伸縮振動(dòng);1 157 cm-1、1 081 cm-1、1 017 cm-1為C—O的伸縮振動(dòng)吸收峰。與圖4(a)原淀粉相比,圖4(b)月桂酸淀粉酯在1 713 cm-1附近出現(xiàn)了一個(gè)新的吸收峰,這是羰基的特征吸收峰,表明淀粉AGU上—OH與月桂酸之間發(fā)生了酯化反應(yīng)。同時(shí),月桂酸淀粉酯在3 419 cm-1左右的吸收峰強(qiáng)度變小,這是淀粉AGU上—OH被酯化形成的CO部分取代的結(jié)果。
圖4 原淀粉(a)和月桂酸淀粉酯(b)的紅外光譜Fig.4 FT-IR spectra: native starch (a) and starch laurate (b)
2.2.2.2 熱重分析
原淀粉和月桂酸淀粉酯的TG及DTG曲線如圖5所示。
圖5 原淀粉(a)和月桂酸淀粉酯(b)的TG及DTG曲線Fig.5 TG and DTG curves: native starch (a); starchlaurate (b)
原淀粉與月桂酸淀粉酯的第1個(gè)熱失重過程發(fā)生于100 ℃左右,這是淀粉中結(jié)晶水蒸發(fā)的結(jié)果;而由于淀粉鏈裂解而產(chǎn)生的第2個(gè)失重過程則發(fā)生于270 ℃左右,且月桂酸淀粉酯的起始裂解溫度低于原淀粉。原淀粉與月桂酸淀粉酯的最大裂解溫度分別為309、320 ℃。通過比較原淀粉和月桂酸淀粉酯的裂解溫度,可以發(fā)現(xiàn)酯化反應(yīng)降低了產(chǎn)物的熱裂解溫度,導(dǎo)致其熱穩(wěn)定性下降。這也進(jìn)一步證實(shí)了淀粉與月桂酸成功發(fā)生了酯化反應(yīng),生成了月桂酸淀粉酯。
(1)以硫酸銨為沉淀劑,以14 mmol/L戊二醛為交聯(lián)劑,交聯(lián)1.5 h,得到的CLEAs相對(duì)酶活力達(dá)到41.89%;與游離酶相比,重復(fù)使用10次后相對(duì)酶活力50.26%,具有良好的重復(fù)使用性。
(2)以普通玉米淀粉和月桂酸為原料,利用脂肪酶CLEAs為催化劑合成了月桂酸淀粉酯,結(jié)果表明,在60 ℃和底物摩爾比為3∶1時(shí),月桂酸淀粉酯的取代度達(dá)到0.151。
(3)通過FT-IR、TGA對(duì)產(chǎn)物進(jìn)行結(jié)構(gòu)表征,證明月桂酸與淀粉發(fā)生了酯化反應(yīng),生成了月桂酸淀粉酯。
[1] CAO Liu-qiu, VAN R F, SHELDON R A. Cross-linked enzyme aggregates: a simple and effective method for the immobilization of penicillin acylase[J]. Organic Letters, 2000, 2(10): 1 361-1 364.
[2] RAJAN A, ABRAHAM T E. Enzymatic modification of cassava starch by bacterial lipase[J]. Bioprocess and biosystems engineering, 2006, 29(1): 65-71.
[3] LIN Ri-hui, LI He, LONG Han, et al. Structure and characteristics of lipase-catalyzed rosin acid starch[J]. Food Hydrocolloids, 2015, 43: 352-359.
[4] HORCHANI H, CHABOUNI M, GARGOURI Y, et al. Solvent-free lipase-catalyzed synthesis of long-chain starch esters using microwave heating: Optimization by response surface methodology[J]. Carbohydrate Polymers, 2010, 79(2): 466-474.
[5] WANG Yan, XIN Jia-ying, SHI Jia, et al. A kinetic study of starch palmitate synthesis by immobilized lipase-catalyzed esterification in solvent free system[J]. Journal of Molecular Catalysis B: Enzymatic, 2014, 101: 73-79.
[7] VORDERWüLBECKE T, KIESLICH K, ERDMANN H. Comparison of lipases by different assays[J]. Enzyme and Microbial Technology, 1992, 14(8): 631-639.
[8] ZARSKI A, PTAK S, SIEMION P, et al. Esterification of potato starch by a biocatalysed reaction in an ionic liquid[J]. Carbohydrate polymers, 2016, 137: 657-663.
[9] MILADINOV V D, HANNA M A. Physical and molecular properties of starch acetates extruded with water and ethanol[J]. Industrial & Engineering Chemistry Research, 1999, 38(10): 3 892-3 897.
[10] AYTAR B S, BAKIR U. Preparation of cross-linked tyrosinase aggregates[J]. Process Biochemistry, 2008, 43(2): 125-131.
[11] KARTAL F, JANSSEN M H, HOLLMANN F, et al. Improved esterification activity ofCandidarugosalipase in organic solvent by immobilization as cross-linked enzyme aggregates (CLEAs)[J]. Journal of Molecular Catalysis B: Enzymatic, 2011, 71(3): 85-89.
Preparation of cross linked lipase aggregates for synthesis of starch laurate
CHEN Hai-long1,2, TIAN Yao-qi1,2, LI Dan1,3, JIN Zheng-yu1,2*
1 (State Key Laboratory of Food Science and Technology,Jiangnan University, Wuxi 214122, China)2 (School of Food Science and Technology,Jiangnan University, Wuxi 214122, China)3(Shaoxing Testing Institute of Food and Drug, Shaoxing 312000, China)
In this work, lipase-cross linked aggregates (lipase-CLEAs) were prepared and the effects of the precipitation and cross-linking on CLEAs activity were investigated. The optimum conditions for immobilization process were obtained as follows: ammonium sulfate was used as precipitant,concentration of cross-linker and enzyme respectively was 14 mmol/L and 30 mg/mL,cross-linking time was 1.5 h. The maximum activity recovery was 41.89% for the lipase-CLEAs. Compared with free enzyme, the thermal stability of CLEAs was largely enhanced and the reusability was significantly improved to 50.26% of their initial activity after 10 recycles. Then maize starch was esterified with laurate by using lipase-CLEAs as catalyst. The effects of the reaction temperature and the molar ratio of lauric acid/anhydroglucose unit (AGU) in starch on the degree of substitution (DS) of starch laurate were studied. The product with the highest DS of 0.151 was prepared under the given conditions. Fourier transform infrared (FTIR) analysis showed that a new peak at 1 713 cm-1for CO stretching in the ester group was observed, which confirmed the formation of starch laurate. Furthermore, decrease of thermal stability of the starch laurate was indicated by thermal gravimetric analysis (TGA).
lipase; Cross-linked enzyme aggregates; esterification; starch laurate
10.13995/j.cnki.11-1802/ts.201702004
碩士研究生(金征宇教授為通訊作者,E-mail:Jinlab2008@yahoo.com)。
國(guó)家自然科學(xué)基金項(xiàng)目(31230057)
2016-09-18,改回日期:2016-11-01