應(yīng)用易錯PCR定向進(jìn)化甘油脫氫酶

李梓君,方柏山,楊仲麗,劉嘉

(華僑大學(xué)工業(yè)生物技術(shù)福建省高等學(xué)校重點(diǎn)實(shí)驗(yàn)室,福建泉州 362021)

應(yīng)用易錯PCR定向進(jìn)化甘油脫氫酶

李梓君,方柏山,楊仲麗,劉嘉

(華僑大學(xué)工業(yè)生物技術(shù)福建省高等學(xué)校重點(diǎn)實(shí)驗(yàn)室,福建泉州 362021)

經(jīng)過連續(xù)易錯聚合酶鏈?zhǔn)椒磻?yīng)(Erro r-Prone PCR)向克雷伯桿菌(K lebsiella pneumoniae)甘油脫氫酶基因中引入突變,并構(gòu)建突變文庫.依據(jù)突變體催化番紅O顯色的特性,采用瓊脂平板法和96微孔板法相結(jié)合的高通量篩選方法,成功獲得突變體gldA-74,其酶活力是親本重組酶的2.73倍.通過軟件對突變體gldA-74酶基因和親本重組酶基因進(jìn)行分析對比,發(fā)現(xiàn)突變體gldA-74酶基因發(fā)生了一處點(diǎn)突變(A 964T),該突變使一處氨基酸被取代.將親本重組酶和進(jìn)化酶的高效表達(dá)產(chǎn)物經(jīng)Ni柱純化后,酶學(xué)性質(zhì)的測定結(jié)果表明,甘油的Km值由0.58 mmol·L-1升高至0.63 mmol·L-1,而NAD的Km值由0.74 mmol·L-1升高至0.77 mmol·L-1,并且酶的最適p H值由原來的11.5升高至12.5,而最適溫度由65℃降至55℃.

甘油脫氫酶;克雷伯桿菌;定向進(jìn)化;易錯聚合酶鏈?zhǔn)椒磻?yīng)

甘油脫氫酶(Glycerol Dehydrogenase,GDH)主要催化甘油生成二羥基丙酮(Dihydroxyacetone, DHA).根據(jù)催化反應(yīng)中產(chǎn)物和電子受體的不同,GDH可分為3種類型[1].第1種為依賴NAD+的GDH[2](EC1.1.1.6),主要存在于各種細(xì)菌中,如Aerobacter aerogenes,Escherichia coli,Bacillus substilis和Cellu lom onas sp.等.第2種為依賴NADP+的 GDH(EC1.1.1.72和EC1.1.1.56),大多存在于霉菌和動物組織中[3].第3種為不依賴輔酶位于細(xì)胞膜上的 GDH,在醋酸菌中發(fā)現(xiàn)存在這類GDH[4].GDH的催化產(chǎn)物DHA因其分子中含3種官能團(tuán),化學(xué)性質(zhì)活潑,廣泛參與聚合、縮合等反應(yīng),所以被廣泛應(yīng)用于輕工、醫(yī)藥和食品等行業(yè)[5].由于應(yīng)用過程中常常需要酶在長時間內(nèi),以及在極端環(huán)境中保持較高活性,所以需要對天然的甘油脫氫酶進(jìn)行改造.酶的定向進(jìn)化是在體外模擬自然進(jìn)化的過程,不需事先了解酶的空間結(jié)構(gòu)和催化機(jī)制,通過人為地創(chuàng)造特殊的條件,模擬自然進(jìn)化機(jī)制(隨機(jī)突變、重組和自然選擇),在體外改造酶基因,并定向選擇出所需性質(zhì)的突變酶.易錯PCR以其操作簡便、有效等優(yōu)點(diǎn)成為目前應(yīng)用最為廣泛的進(jìn)化手段之一[6].本研究利用易錯PCR對克雷伯桿菌中 GDH進(jìn)行定向進(jìn)化,以期獲得酶活大幅提高,酶學(xué)性質(zhì)有所改善的突變酶.

1 實(shí)驗(yàn)材料和方法

1.1 菌株與質(zhì)粒

E.coliBL21(DE3)(p ET-32gldA)的構(gòu)建,以及表達(dá)載體p ET-32a(+)、菌株和質(zhì)粒的保存,均由華僑大學(xué)工業(yè)生物技術(shù)福建省高等學(xué)校重點(diǎn)實(shí)驗(yàn)室完成.

1.2 試劑

LB培養(yǎng)基.T4 DNA連接酶、dN TPs,限制酶、r Taq酶,購自大連TaKaRa公司.膠回收試劑盒、質(zhì)粒提取試劑盒均購自美國Omega公司,輔酶Ⅰ購自Amersco公司,其余試劑均為國產(chǎn)分析純.

1.3 定向進(jìn)化方法

1.3.1 甘油脫氫酶基因突變體文庫構(gòu)建 (1)引物.PrimerⅠ為5’-CGTCGGA TCCTACA TGCGCACTTA TTTGAG-3’);PrimerⅡ(5’-AA TGCTCGAGCGAA TTAACGCGCCAGCCAC-3’).

(2)易錯PCR反應(yīng)體系.20μL體系中含1×Taq DNA聚合酶緩沖液、0.2 mmol·L-1dA TP和dGTP,1 mmol·L-1dCTP和d TTP,6.0～8.5 mmol·L-1M g2+,0.3～0.8 mmol·L-1M n2+,1 mmol·L-1PrimerⅠ和PrimerⅡ,模板DNA為質(zhì)粒p ET-32gldA約50 ng.其中,M g2+和M n2+的濃度被調(diào)整,可以得到不同突變頻率的文庫.

(3)擴(kuò)增與鑒定.94℃,5 min;94℃,1 min;62℃,1.5 min;72℃,1.5 min;30個循環(huán);72℃,15 min.用質(zhì)量分?jǐn)?shù)為0.8%的瓊脂糖凝膠電泳鑒定易錯PCR產(chǎn)物,利用膠回收試劑盒切膠回收PCR產(chǎn)物,于-20℃保存.

(4)文庫構(gòu)建.將上述易錯PCR產(chǎn)物經(jīng)膠回收純化后,經(jīng)Bam H I-XhoⅠ雙酶切,與同樣雙酶切的載體p ET-32a(+)相連接,重組子引入E.coliBL 21(DE3),構(gòu)建突變體文庫.

1.3.2 突變體的篩選和鑒定 (1)篩選模型的建立[7].將菌體用無菌水稀釋,涂布于番紅O篩選平板(胰蛋白胨6.0 g·L-1,酵母浸膏 3.0 g·L-1,氯化鈉 10.0 g·L-1,甘油 18.4 g·L-1,番紅O 0.1 mmol·L-1,IPTG 1.0 mmol·L-1,氨芐青霉素75.0 mg·L-1,瓊脂粉10.0 g·L-1),觀察平板的單菌落的生長及變色情況.挑取顏色顯示均勻桃紅色的單菌落置于含75 mg·L-1氨芐青霉素的LB液體培養(yǎng)基中培養(yǎng),經(jīng)乳糖誘導(dǎo)后收集菌體進(jìn)行復(fù)篩鑒定.

(2)陽性轉(zhuǎn)化子的鑒定.按質(zhì)粒提取試劑盒方法抽提正向重組質(zhì)粒,經(jīng) HindⅢ,PstⅠ單雙酶切,用質(zhì)量分?jǐn)?shù)為0.8%瓊脂糖凝膠電泳鑒定酶切產(chǎn)物.

1.4 誘導(dǎo)表達(dá)

挑取陽性克隆單菌落接種于含有75 m g·L-1氨芐青霉素的LB培養(yǎng)基中,于30℃培養(yǎng)至吸光度值(D)約0.4時,加入乳糖至終質(zhì)量濃度為2 g·L-1,于30℃下誘導(dǎo)培養(yǎng)5 h后離心收集菌體.

1.5 進(jìn)化酶酶活力及菌體蛋白含量測定

酶活力測定具體方法參照文[8].GDH使甘油脫氫生成二羥基丙酮,同時輔酶 ⅠNAD+被還原為NADH.NADH在340 nm的紫外光下有最大吸收,根據(jù)吸光度的增加情況,用初速度法可測定 GDH的活力.取1.5 mL反應(yīng)液(含30 mmol·L-1(NH4)2SO4,0.2 mol·L-1甘油,2 mmol·L-1NAD,1 μmol·L-1Fe(NH4)2(SO4)2,0.1 mol·L-1碳酸鉀緩沖溶液,p H=12.0),于45℃下恒溫后,加入適量酶液啟動反應(yīng),在340 nm波長下連續(xù)測定酶反應(yīng)過程中的吸光值(每隔6 s記錄一個數(shù)據(jù),時間1 min),通過計(jì)算可得到GDH的酶活(z).酶活的定義:在上述條件下,每分鐘還原1μmol甘油的酶量為1個酶活力單位.蛋白含量按Bradford法測定[9],以牛血清白蛋白為標(biāo)準(zhǔn)蛋白質(zhì).

1.6 突變基因序列測定

對篩選到的陽性克隆進(jìn)行培養(yǎng)復(fù)篩,挑選在高p H值條件下仍有較高酶活力的突變體,提取質(zhì)粒并用自動序列儀進(jìn)行測序.測序由遼寧省大連寶生物公司完成.

1.7 進(jìn)化酶的純化

按照常規(guī)方法進(jìn)行鎳柱親和層析純化[8].SDS-PAGE(十二烷基磺酸鈉-聚丙烯酰胺膠體電泳):采用12%分離膠的不連續(xù)垂直平板電泳,考馬斯亮藍(lán)G-250染色.

2 實(shí)驗(yàn)結(jié)果

2.1 易錯聚合酶鏈?zhǔn)椒磻?yīng)

在高離子濃度M g2+或M n2+存在的條件下,Taq DNA聚合酶的保真性大大降低.在延伸的過程中,一些堿基被錯誤地引入到基因中,從而使基因發(fā)生突變.通過改變M g2+和M n2+的濃度,可以得到不同突變頻率的片段多樣性文庫.在其他成分恒定的條件下,選擇不同濃度的M g2+,M n2+進(jìn)行易錯PCR,結(jié)果如圖1所示.圖1(a)中,1～6表示M g2+濃度分別為6.0,6.5,7.0,7.5,8.0,8.5 mmol·L-1;圖1(b)中,1～6表示M n2+濃度分別為0.3,0.4,0.5,0.6,0.7,0.8 mmol·L-1.

從圖1可知,M gCl2濃度范圍在7.0～8.5 mmol·L-1間均適合該基因易錯PCR反應(yīng).為了確保易錯PCR的突變率控制在每個基因1～3個堿基突變,選擇M gCl2濃度7.0 mmol·L-1作為該基因易錯PCR反應(yīng)體系的最適濃度.在以上實(shí)驗(yàn)得到的最適M gCl2濃度的基礎(chǔ)上,進(jìn)一步考察不同M n2+濃度對反應(yīng)體系的影響.結(jié)果表明,M n2+最適濃度為0.5 mmol·L-1.

圖1 不同離子濃度下的易錯PCRFig.1 Erro r-p rone PCR conducted with different ion concentration

2.2 突變體文庫的構(gòu)建及突變株的篩選和鑒定

在確定的易錯PCR條件下擴(kuò)增甘油脫氫酶基因,經(jīng)連續(xù)易錯PCR,產(chǎn)物經(jīng)純化后用Bam HⅠ-XhoⅠ雙酶切,與經(jīng)過同樣雙酶切的表達(dá)載體p ET-32a(+)相連接.連接產(chǎn)物轉(zhuǎn)化至E.coliBL 21(DE3)感受態(tài)細(xì)胞中構(gòu)建突變文庫.研究中共獲得5 000余株的突變株.

2.2.1 突變株的初篩 將轉(zhuǎn)化后的E.coliBL 21(DE3)細(xì)胞涂布于番紅O篩選平板上,在37℃下培養(yǎng)約12 h,對突變體進(jìn)行篩選.甘油脫氫酶在催化甘油生成二羥基丙酮的同時,把 NAD+轉(zhuǎn)化成NADH.所生成的NADH會通過呼吸鏈?zhǔn)闺妱萆?產(chǎn)生強(qiáng)質(zhì)子流,增強(qiáng)了大腸桿菌細(xì)胞對番紅O的攝入,并使得番紅O在大腸桿菌內(nèi)聚集[7,10].通過番紅O篩選平板涂布,具有甘油脫氫酶活性的克隆顯示均勻的桃紅色,而無活性的克隆外圍則有一個透明圈,如圖2所示.

圖2 番紅O篩選平板初篩結(jié)果Fig.2 The screening result in safranin O plates

2.2.2 突變株的復(fù)篩 將初篩獲得的克隆子接種于LB培養(yǎng)基(含75 mg·L-1氨芐青霉素)中,在30℃下培養(yǎng)至D約0.4時,加入乳糖至終質(zhì)量濃度為2 g·L-1;30℃誘導(dǎo)培養(yǎng)5 h后,超聲破碎菌體并離心收集上清液.采用96孔板結(jié)合酶標(biāo)儀,測定初篩獲得的菌株在不同p H值下的酶活,最終獲得最佳突變體gldA-74,其甘油脫氫酶活力為2.202 m kat·L-1.

相比由張婷婷等[8]構(gòu)建的親本重組酶(酶活為0.805 mkat·L-1),酶活提高2.73倍;而相比于陳宏文等[1]由出發(fā)菌克雷伯桿菌經(jīng)培養(yǎng)基和培養(yǎng)條件優(yōu)化后,所測得無細(xì)胞抽提液的酶活0.062 m kat· L-1,酶活提高35.7倍.

2.2.3 突變株酶切鑒定結(jié)果 將均勻桃紅色的克隆子在LB液體培養(yǎng)基中進(jìn)行培養(yǎng),然后抽提重組質(zhì)粒.重組質(zhì)粒用HindⅢ進(jìn)行單酶切鑒定時,應(yīng)出現(xiàn)約6 987 bp大小的一條帶;而用 HindⅢ/PstⅠ進(jìn)行雙酶切鑒定時,由于gldA-74上還存在兩個PstⅠ酶切位點(diǎn),應(yīng)得到兩條DNA條帶,其大小分別約為4 808,1 584 bp.

酶切鑒定結(jié)果表明連接正確,如圖3所示.圖3中:1,2分別為經(jīng) HindⅢ和經(jīng) HindⅢ/PstⅠ進(jìn)行酶切p ET-32gldA-74.

2.2.4 進(jìn)化酶與親本重組酶序列比較結(jié)果 突變后的基因委托上海生工生物工程技術(shù)服務(wù)有限公司進(jìn)行測序,將進(jìn)化酶與親本重組酶基因序列進(jìn)行對比.測序結(jié)果表明,利用易錯PCR的方法,成功地向甘油脫氫酶基因引入了突變,其中酶活力增高、最適p H值增大的gldA-74突變體酶基因有1個堿基發(fā)生了突變(即964位的A變?yōu)門).這處突變使得322位的I(異亮氨酸)突變?yōu)镕(苯丙氨酸).

2.3 甘油脫氫酶的分離與純化

選用的表達(dá)載體p ET-32a(+)中含有表達(dá)His6-Tagged的基因序列,進(jìn)化酶含有該融合標(biāo)簽,可以用鎳柱純化.另外,在該進(jìn)化酶的N端還含有一個Rrx-Tag TM,該融合標(biāo)簽有利于提高目的蛋白的可溶部分比例及活性.研究中,進(jìn)化酶經(jīng)鎳柱親和層析純化后可見單一條帶,而穿透液幾乎不含該條帶,基本上可以認(rèn)為得到比較純的蛋白.

當(dāng)以p ET-32a(+)為表達(dá)載體時,表達(dá)的目標(biāo)蛋白所含的融合基因的相對分子質(zhì)量約為20.4 ku,而研究中的 GDH的相對分子質(zhì)量約為34 ku[1].因此,突變質(zhì)粒所表達(dá)的融合蛋白的相對分子質(zhì)量應(yīng)為54 ku,如圖4所示.圖4中:1,2,3分別為純蛋白、穿透液和粗蛋白.

表1為純化結(jié)果.表1中:mt為總蛋白量;zt為總活力;ρ為蛋白質(zhì)量濃度;σ為比活力;η為回收率;n為純化倍數(shù).從表1可以看出,純化后進(jìn)化酶的回收率為24.1%,純化倍數(shù)是2.1.

圖3 表達(dá)載體酶切鑒定Fig.3 Enzyme digestion of exp ression vector

圖4 進(jìn)化酶的純化Fig.4 Purification of evolved emzyme

表1 進(jìn)化GDH的純化Tab.1 Purification of evolved GDH

2.4 甘油脫氫酶的酶學(xué)性質(zhì)分析

2.4.1 酶動力學(xué)參數(shù)比較 親本重組酶與進(jìn)化酶都以1/V和1/S作Lineweaver Burk雙倒數(shù)圖,結(jié)果如圖5所示.由此計(jì)算,得到親本重組酶的Km(甘油)為0.58 mmol·L-1,Km(NAD)為0.74 mmol· L-1,Vmax(甘油)為0.75 mmol·(L·min)-1,Vmax(NAD)為0.96 mmol·(L·min)-1;進(jìn)化酶的Km(甘油)為0.63 mmol·L-1,Km(NAD)為0.77 mmol·L-1,Vmax(甘油)為3.69 mmol·(L·m in)-1,Vmax(NAD)為4.75 mmol·(L·min)-1.

2.4.2 最適溫度比較 親本重組酶與進(jìn)化酶的最適反應(yīng)溫度,如圖6所示.從圖6可見,親本重組酶的最適溫度為65℃,說明GDH能在較高的溫度下進(jìn)行反應(yīng).但是,再繼續(xù)升高溫度,酶活力迅速下降;在70℃時,殘留的相對酶活力僅有20%.由圖6又可知,進(jìn)化酶的最適溫度為55℃,比親本重組酶的最適溫度降低了10℃.在50～60℃范圍內(nèi),進(jìn)化酶的相對酶活均維持在80%以上.

2.4.3 最適p H值的比較 在45℃,不同p H值的緩沖液中,檢測親本重組酶和進(jìn)化酶活力,結(jié)果如圖7所示.由圖7可以看出,甘油脫氫酶的酶活隨著反應(yīng)液p H值的變化而變化.親本重組酶在p H值為11.0時,甘油脫氫酶的酶活最大;而在p H值為10～12時,GDH的相對活性能保持在70%以上.

圖5 酶動力學(xué)參數(shù)的Line weaver Burk雙倒數(shù)圖Fig.5 Enzyme kinetic parameters of the Line weaver Burk double-reciprocal graph

由圖7還可以看出,進(jìn)化酶在p H值為12.5時酶活最大,比親本重組酶的最適p H值高.在p H值為10～13時,進(jìn)化酶的相對活性能保持在70%以上.中性環(huán)境酶活損失很大,在p H值為7的條件下,酶活僅為20%左右.進(jìn)化酶催化反應(yīng)耐強(qiáng)堿,說明 GDH存有較大的工業(yè)應(yīng)用價值.

圖7 pH值對酶活力的影響Fig.7 the effect of reaction pH on the specific activity of recombinant GDH and evolved GDH

圖6 溫度對酶活力的影響Fig.6 Effect of reaction temperature on the activity of recombinant GDH and evolved GDH

3 討論

易錯PCR操作簡單,容易引入點(diǎn)突變.其突變頻率高于傳統(tǒng)的物理和化學(xué)誘變,是目前定向進(jìn)化領(lǐng)域最為常用的隨機(jī)突變手段.因此,采用易錯PCR對GDH進(jìn)行定向進(jìn)化.

首先運(yùn)用番紅O篩選平板對5 000余株突變體進(jìn)行初篩;然后,再利用96孔板結(jié)合酶標(biāo)儀測定經(jīng)初篩所得的菌株在不同p H值條件下的酶活,最終成功獲得突變體gldA-74.該突變體有一個有益突變(即I322F),使得進(jìn)化酶的酶活是親本重組酶的2.73倍.

將親本重組酶和進(jìn)化酶的高效表達(dá)產(chǎn)物經(jīng)鎳柱純化后,其酶學(xué)性質(zhì)測定表明:Km(甘油)由0.58 mmol·L-1升高至0.63 mmol·L-1,Km(NAD)由0.74 mmo l·L-1升高至0.77 mmol·L-1,并且酶的最適p H值也由11.5升高至12.5,最適溫度由65℃降至55℃.綜上所述,最終獲得的突變體gldA-74與親本重組酶相比,具有較高的催化活力和偏堿的最適p H值,更適合在偏堿的環(huán)境中應(yīng)用.

今后的研究方向有以下幾個方面:(1)應(yīng)用生物信息學(xué)工具,對定向進(jìn)化后的突變基因庫、蛋白質(zhì)序列、3D結(jié)構(gòu)、活性與功能的變化等做比對分析;(2)從根本上揭示突變酶的催化機(jī)理和折疊機(jī)制、蛋白質(zhì)一級序列與三級結(jié)構(gòu)的關(guān)系,以及結(jié)構(gòu)與功能的關(guān)系等,為探索酶法同時合成DHA和1,3-PD的新途徑提供理論依據(jù).

[1]陳宏文,吳雅紅,吳振華,等.克雷伯桿菌甘油脫氫酶的分離純化及性質(zhì)[J].無錫輕工大學(xué)學(xué)報(bào),2005,24(1):1-5.

[2]RUXHEIN IKOV 1 S N,BURKE1 J,SEDELN IKOVA 1 S,et al.Glycerol dehydrogenase structure,specificity, and mechanism of a familyⅢpolyol dehydrogenase[J].Structure,2001,9(9):789-802.

[3]YAMADA H,NAGAO A,N ISH ISE H,et al.Fo rmation of glycerol dehydrogenase by microorganism s[J].Agric Biol Chem,1982,46(9):2325-2331.

[4]AM EYAMA M,SH INAGAWA E,MA T SUSH ITA K,et al.Solubilization,purification and properties of membrane-bound glycerol dehydrogenase fromGluconobacter industrius[J].Agric Biol Chem,1985,49(4):1001-1010.

[5]KAWASH IMA K,ITOH H,CHGA TEJ.Nonenzymatic browning reactions of dihydroxyacetone with amino acids o r their esters[J].Agri Biol Chem,1980,44(7):1595-1599.

[6]MOORE J C,JIN H M,KUCHNER O,et al.Strategies for thein vitroevolution of protein function:Enzyme evolution by random recombination of imp roved sequences[J].Journal of Molecular Biology,1997,272(3):336-347.

[7]LOO V B,HARALD J,SPELBERGL.Directed evolution of epoxide hydrolase fromA.radiobactertoward higher enantioselectivity by erro r-p rone PCR and DNA shuffling[J].Chemistry&Biology,2004,11(7):981-990.

[8]張婷婷,方柏山,王耿,等.克雷伯桿菌甘油脫氫酶基因的克隆表達(dá)與純化[J].生物工程學(xué)報(bào),2008,24(3):495-499.

[9]汪家政,汪明.蛋白質(zhì)技術(shù)手冊[M].北京:科學(xué)出版社,2000:42-47.

[10]方柏山,鄭媛媛.酶體外定向進(jìn)化(Ⅰ):突變基因文庫構(gòu)建技術(shù)及其新進(jìn)展[J].華僑大學(xué)學(xué)報(bào):自然科學(xué)版, 2004,25(4):337-342.

(責(zé)任編輯:黃仲一英文審校:陳國華)

Directed Evolution of Glycerol Dehydrogenase by Error Prone PCR

L IZi-jun,FANGBai-shan, YANG Zhong-li,L IU Jia
(Key Laboratory of Industrial Biotechnology Fujian Province,Huaqiao University,Quanzhou 362021,China)

The GDH enzyme gene fromKlebsiella pneumoniaewas repeatedly amp lified by two sequential error-prone polymerase chain reaction(PCR),and the gene library was built.The mutant can catalyze safranin O to show colors,and on the base of this property,the gldA-74 mutant was obtained by high throughput screening method using active agar plates in combination with 96-wellmicrop lates.The activity of the gldA-74 mutant was showed 2.73 fold of that of the parent recombinase.Gene analysis of the gldA-74mutant indicated that the mutant enzyme had a pointmutation(A 964T) which resulted in an amino acid being rep laced.Then,the highly exp ressed productions of recombinase GDH and evolved GDH were purified by Ni-NTA and the enzymatic properties were determined.The results showed that theKmof GDH increased from 0.58 mmol·L-1to 0.63 mmol·L-1,and theKmof NAD increased as well,from 0.71 mmol·L-1to 0.77 mmol·L-1.The op timum pH of mutant enzyme also increased from 11.5 to 12.5 while the op timum temperature decreased from 65℃to 55℃.

glycerol dehydrogenase;Klebsiella pneumoniae;directed evolution;error-prone polymerase chain reaction

Q 554+.903

A

1000-5013(2010)06-0661-06

2008-12-28

方柏山(1957-),男,教授,主要從事生物化工的研究.E-mail:fangbs@hqu.edu.cn.

國家重點(diǎn)基礎(chǔ)研究發(fā)展(973)計(jì)劃項(xiàng)目(2006AA 020103);國家自然科學(xué)基金資助項(xiàng)目(20676048)