張召寶， 侯 林， 崔清華， 潘 晴， 馬魯豫

（1.山東中醫(yī)藥大學(xué)藥學(xué)院，山東濟(jì)南250300；2.山東省醫(yī)學(xué)科學(xué)院，山東濟(jì)南250062）

中藥三萜皂苷合成通路的生物信息學(xué)分析

張召寶1， 侯 林1， 崔清華1， 潘 晴1， 馬魯豫2*

（1.山東中醫(yī)藥大學(xué)藥學(xué)院，山東濟(jì)南250300；2.山東省醫(yī)學(xué)科學(xué)院，山東濟(jì)南250062）

目的從系統(tǒng)進(jìn)化學(xué)角度對(duì)中藥三萜皂苷合成通路中關(guān)鍵酶進(jìn)行生物信息學(xué)分析。方法從NCBI數(shù)據(jù)庫下載已有17種中藥中的鯊烯合成酶、鯊烯環(huán)氧酶、2，3-氧化鯊烯環(huán)化酶的蛋白質(zhì)全長(zhǎng)序列，通過生物信息學(xué)軟件對(duì)其進(jìn)行理化性質(zhì)分析、保守域分析、蛋白質(zhì)三維結(jié)構(gòu)預(yù)測(cè)和系統(tǒng)進(jìn)化樹的構(gòu)建。結(jié)果序列分析表明，鯊烯合成酶、鯊烯環(huán)氧酶、2，3-氧化鯊烯環(huán)化酶的蛋白質(zhì)序列保守性均較高。其中鯊烯合成酶具有兩個(gè)高度保守的區(qū)域——富含天冬氨酸的鎂離子結(jié)合位點(diǎn)1和2——可作為認(rèn)定鯊烯合成酶的依據(jù)；鯊烯環(huán)氧酶和2，3-氧化鯊烯環(huán)化酶的蛋白三維結(jié)構(gòu)保守，無規(guī)則卷曲部分為高變區(qū)。結(jié)論上游合成在進(jìn)化中趨于保守和穩(wěn)定；三萜皂苷眾多類型的形成可能與下游細(xì)胞色素P450、糖基轉(zhuǎn)移酶、β-糖苷酶等的修飾有關(guān)。

三萜皂苷；鯊烯合成酶；鯊烯環(huán)氧酶；2，3-氧化鯊烯環(huán)化酶

三萜皂苷（triterpenoid saponins）是一類重要的植物次生代謝產(chǎn)物，廣泛分布于人參［1］、靈芝［2］、三七［3］、墨旱蓮［4］等中藥中，并作為其主要的活性成分，具有抗癌、抗氧化、消炎解毒等功效。

目前，人們對(duì)三萜皂苷的合成通路已有較為明確的認(rèn)識(shí)［5］。首先通過甲羥戊酸等途徑合成異戊二烯焦磷酸，然后經(jīng)過香葉基二磷酸合酶、法尼基焦磷酸合酶、鯊烯合成酶和鯊烯環(huán)氧酶的催化，形成2，3-氧化鯊烯，再通過2，3-氧化鯊烯環(huán)化酶環(huán)化形成三萜類骨架，最后通過細(xì)胞色素P450、糖基轉(zhuǎn)移酶、β-糖苷酶的修飾作用，形成復(fù)雜的糖苷化合物——三萜皂苷。

三萜皂苷主要從栽培藥材中提取，但存在種類少、產(chǎn)量低、受環(huán)境影響大、收獲周期長(zhǎng)等問題，難以滿足市場(chǎng)需求［6］。本研究系統(tǒng)分析了三萜皂苷合成通路中的關(guān)鍵酶，可增加對(duì)其進(jìn)化規(guī)律的認(rèn)識(shí)，發(fā)掘更多的相關(guān)藥用資源，為今后利用基因工程技術(shù)實(shí)現(xiàn)其大規(guī)模人工合成打下堅(jiān)實(shí)基礎(chǔ)。

1 資料來源

截止到2014年3月27日，從NCBI數(shù)據(jù)庫下載鯊烯合成酶蛋白質(zhì)全長(zhǎng)序列，共31條，包括人參Panax ginseng C.A.Mey.、三七Panax notoginseng（Burk.）F.H.Chen等17種中藥的24條序列；下載鯊烯環(huán)氧酶蛋白質(zhì)全長(zhǎng)序列，共19條，包含人參、蛇足石杉Huperzia serrata（Thunb.ex Murray）Trev.等5種中藥材的7條序列；下載2，3-氧化鯊烯環(huán)化酶蛋白質(zhì)全長(zhǎng)序列，共8條，包括人參、蒲公英Taraxacum mongolicum Hand.-Mazz.、總狀升麻Cimicifuga racemosa（Turcz.）Maxim.和百脈根Lotus corniculatus這4種中藥材的4條序列。相關(guān)序列的NCBI號(hào)在系統(tǒng)進(jìn)化樹中列出。

2 方法

2.1 蛋白質(zhì)理化性質(zhì)分析 蛋白質(zhì)全長(zhǎng)序列的理化性質(zhì)應(yīng)用在線軟件ExPASy（http：//www.expasy.org/）分析。其中，利用ProtParam模塊根據(jù)具體的氨基酸序列，進(jìn)行氨基酸數(shù)目、相對(duì)分子質(zhì)量、理論等電點(diǎn)和不穩(wěn)定指數(shù)的測(cè)算，然后利用ProtScale模塊進(jìn)行親水性/疏水性分析。

2.2 結(jié)構(gòu)域分析 在NCBI的Conserved Domain Database［7］和SMART［8］（http：//smart.embl-heidelberg.de/）數(shù)據(jù)庫中進(jìn)行蛋白質(zhì)的保守域和功能域預(yù)測(cè)，并通過WebLogo［9］對(duì)其做進(jìn)一步分析。利用基于同源建模的分析工具SWISS-MODEL［10］（http：//swissmodel.expasy.org/）預(yù)測(cè)蛋白質(zhì)的三維結(jié)構(gòu)。

2.3 構(gòu)建系統(tǒng)進(jìn)化樹 利用MEGA 5.0軟件［11］中的ClustalW進(jìn)行序列比對(duì)，然后利用Neighbor-Joining（NJ）法構(gòu)建系統(tǒng)進(jìn)化樹（bootstrap= 1000）。

3 結(jié)果與分析

3.1 鯊烯合成酶 鯊烯合成酶（squalene synthase，EC 2.5.1.21）是一種定位于內(nèi)質(zhì)網(wǎng)上的膜結(jié)合蛋白［12］。在類異戊二烯代謝途徑中，它以法尼二磷酸為底物，催化合成30C產(chǎn)物角鯊烯。由于處于甾醇與三萜合成通路的分支點(diǎn)，決定著法尼基焦磷酸的流向，因此鯊烯合成酶發(fā)揮著重要的調(diào)節(jié)作用。

已有研究表明，鯊烯合成酶的結(jié)構(gòu)和反應(yīng)機(jī)制在動(dòng)植物中相當(dāng)保守。見表1。

表1 藥用植物鯊烯合成酶蛋白質(zhì)序列的理化性質(zhì)Tab.1 Physical and chem ical properties of medicinal plants squalene synthase protein sequences

在哺乳動(dòng)物中，該成分大約包含416 aa，相對(duì)分子質(zhì)量約為47 kDa；在模式植物擬南芥中，該成分有2個(gè) （NP_195190.1和NP_195191.2），其外顯子具有高度同源性［13］。本研究從NCBI數(shù)據(jù)庫中共下載得到17種中藥的24條蛋白質(zhì)全長(zhǎng)序列，長(zhǎng)度在405～418 aa之間，相對(duì)分子質(zhì)量在45.89～47.89 kDa之間，理論等電點(diǎn)在6.07～8.53之間，多數(shù)為不穩(wěn)定蛋白質(zhì) （不穩(wěn)定指數(shù)＞40）。以人參中的三條鯊烯合成酶序列（BAD08242.1，ACV88718.1和ACZ71037.1）為例，利用ProtScale模塊進(jìn)行疏水性/親水性分析，由于其序列高度保守，所展現(xiàn)出的親水性/疏水性規(guī)律較為一致，見圖1。另外由圖可知，鯊烯合成酶的N端分布有親水域 （分值＜-0.5），而C端則分布有一段較短的強(qiáng)疏水域 （分值＞0.5），這與之前有關(guān)鯊烯合成酶細(xì)胞定位的報(bào)道相一致。鯊烯合成酶通過短的C端跨膜區(qū)錨定于內(nèi)質(zhì)網(wǎng)上，而N端的催化域則伸向胞質(zhì)溶膠中［14］。

圖1 鯊烯合成酶疏水性/親水性分析Fig.1 Hydroponic/hydrophilic analysis of squalene synthases in P.ginseng

在對(duì)其序列進(jìn)行理化性質(zhì)分析的基礎(chǔ)上，進(jìn)一步研究其保守域。首先采用WebLogo對(duì)Lee［15］等報(bào)道的6個(gè)保守域進(jìn)行分析，結(jié)果發(fā)現(xiàn)在17種中藥中，Ⅲ、Ⅳ、Ⅴ三個(gè)保守域高度保守，而Ⅰ、Ⅱ、Ⅵ三個(gè)保守域保守性較差，見圖2 A。由此可見，之前對(duì)這6個(gè)保守域的報(bào)道具有一定局限性。根據(jù)NCBI的Conserved Domain Database的預(yù)測(cè)結(jié)果比對(duì)發(fā)現(xiàn)，這17種中藥中的鯊烯合成酶具有兩個(gè)高度保守 （100%）的區(qū)域——富含天冬氨酸的鎂離子結(jié)合位點(diǎn)1和2，見圖2 B（結(jié)合位點(diǎn)1位于77～81 aa，結(jié)合位點(diǎn)2位于213～217 aa），該區(qū)域與15C法尼基和20C香葉基異戊二烯基二磷酸的1＇-1＇位縮合反應(yīng)相關(guān)，其保守性保證了鯊烯合成酶與底物特異性結(jié)合，并可作為今后鑒定中藥中該成分的重要依據(jù)。以人參中的三條鯊烯合成酶序列（BAD08242.1，ACV88718.1和ACZ71037.1）為例，SWISS-MODEL建模的結(jié)果顯示，人參中鯊烯合成酶的三維結(jié)構(gòu)相對(duì)保守，以α螺旋為主，而β折疊片較少，無規(guī)則卷曲等結(jié)構(gòu)均很保守（圖2 C）。

圖2 鯊烯合成酶保守域分析Fig.2 Conserved domain analysis of squalene synthases

以模式物種擬南芥等植物中的鯊烯合成酶蛋白質(zhì)全長(zhǎng)構(gòu)建進(jìn)化樹，見圖3。由圖可知，進(jìn)化樹首先分為兩大支，其中上面一支主要包含較高等的植物類群，主要為雙子葉植物；下面一支以紫杉Taxus chinensis（Pilger）Rehd.、馬尾松Pinusmassoniana為代表，多為較低等的植物類群。整體而言，鯊烯合成酶具有較高的保守性和較為普遍的基因倍增現(xiàn)象，例如人參、丹參、北柴胡等6種中藥具有兩個(gè)及兩個(gè)以上鯊烯合成酶。由于其具有較為重要的生理生化功能，故認(rèn)為適應(yīng)性進(jìn)化可能導(dǎo)致了普遍的基因倍增，從而產(chǎn)生了不同的亞型，之間可能存在的功能上的互補(bǔ)和增強(qiáng)，又是酶功能進(jìn)化的一種體現(xiàn)。

圖3 鯊烯合成酶系統(tǒng)進(jìn)化樹 （NJ法）Fig.3 Phylogenetic tree of squalene synthases

3.2 鯊烯環(huán)氧酶 鯊烯環(huán)氧酶（squalene epoxidase，EC 1.14.13.132）是一種細(xì)胞膜結(jié)合蛋白，催化三萜合成的第一步氧化反應(yīng)，是該合成通路的限速酶之一［16］。

從NCBI中共下載得到7條中藥鯊烯環(huán)氧酶的蛋白質(zhì)全長(zhǎng)序列，其長(zhǎng)度在537～561 aa之間，相對(duì)分子質(zhì)量在59.13～62.79 kDa之間，理論等電點(diǎn)在8.63～9.19之間，全部為不穩(wěn)定蛋白質(zhì) （不穩(wěn)定指數(shù)＞40）（見表2）。

表2 藥用植物鯊烯環(huán)氧酶蛋白質(zhì)序列的理化性質(zhì)Tab.2 Physical and chem ical properties of medicinal p lants squalene epoxidase protein sequences

以三七和人參中的鯊烯環(huán)氧酶序列為例，利用ProtScale模塊進(jìn)行疏水性/親水性分析，可以看出，三七和人參中的兩條鯊烯環(huán)氧酶蛋白質(zhì)全長(zhǎng)序列出現(xiàn)了較為明顯的分化；三七中的ABE73759.1與人參的BAA24448.1具有較為一致的親水性/疏水性規(guī)律；三七中的 AFV92748.1與人參的ACJ24907.2具有較為一致的親水性/疏水性規(guī)律；這說明這兩種亞型之間的分化早于三七和人參之間的分化，并在隨后的進(jìn)化中形成了自己獨(dú)特的親水性/疏水性規(guī)律 （圖4）。

圖4 鯊烯環(huán)氧酶疏水性/親水性分析Fig.4 Hydroponic/hydrophilic analysis of squalene epoxidases in P.notoginseng and P.ginseng

對(duì)其進(jìn)行功能域分析，見圖5 A。由圖可知，鯊烯環(huán)氧酶的C端含有一個(gè)高度保守的跨膜結(jié)構(gòu)域（5～27 aa），然后是1個(gè)重疊的FAD結(jié)合域（NADPH結(jié)合域，60～408 aa）和鯊烯環(huán)氧酶結(jié)構(gòu)域 （pfam08491，211～484 aa）。藥用植物中的鯊烯環(huán)氧酶蛋白質(zhì)序列呈現(xiàn)出較高的同源性，以模式植物擬南芥和中藥代表物種人參為例，通過SWISS-MODEL進(jìn)行三維結(jié)構(gòu)的預(yù)測(cè)，見圖5 B。由圖可知，鯊烯環(huán)氧酶的α螺旋和β折疊片等結(jié)構(gòu)保守，無規(guī)則卷曲部分為可變區(qū)，其中連接β折疊片者較為高變。

以模式物種擬南芥等植物中的鯊烯環(huán)氧酶蛋白質(zhì)全長(zhǎng)構(gòu)建進(jìn)化樹，見圖6。由圖可知，進(jìn)化地位較為低等的蕨類植物——蛇足石杉作為外群，位于進(jìn)化樹的底部；五加科的人參、三七、刺五加等植物聚為一支，顯示出較近的親緣關(guān)系。該進(jìn)化樹較好地符合物種之間的進(jìn)化規(guī)律，說明鯊烯環(huán)氧酶伴隨物種的進(jìn)化進(jìn)程趨于保守，只在特定類群中出現(xiàn)了基因倍增的現(xiàn)象，如五加科，這可能與該科物種特定的生境與生理習(xí)性相關(guān)，例如，某種環(huán)境壓力迫使其特定通路的完善或者特定代謝產(chǎn)物的富集。

圖5 鯊烯環(huán)氧酶保守域及三維結(jié)構(gòu)分析Fig.5 Conserved domains and 3D structure analysis of squalene epoxidases

圖6 鯊烯環(huán)氧酶系統(tǒng)進(jìn)化樹 （NJ法）Fig.6 Phylogenetic tree of squalene epoxidases

3.3 2，3-氧化鯊烯環(huán)化酶 氧化鯊烯環(huán)化酶（oxidosqualene cyclases，EC 5.4.99.7）是一大類氧化鯊烯環(huán)化酶的統(tǒng)稱，能催化一系列復(fù)雜的環(huán)化和重排反應(yīng)。在植物中，由于它可催化三萜皂苷合成中的第一步環(huán)化反應(yīng)，因而被認(rèn)為是合成通路中的關(guān)鍵酶。

表3 藥用植物氧化鯊烯環(huán)化酶蛋白質(zhì)序列的理化性質(zhì)Tab.3 Physical and chem ical properties of medicinal p lants oxidosqualene cyclase p rotein sequences

本研究對(duì)其中的2，3-氧化鯊烯環(huán)化酶進(jìn)行深入的分析，發(fā)現(xiàn)它在真菌、動(dòng)物和植物中均較為保守。從NCBI中共下載得到4條中藥的2，3-氧化鯊烯環(huán)化酶全長(zhǎng)蛋白序列，包括人參、百脈根、蒲公英、總狀升麻，長(zhǎng)度在757～780 aa之間，相對(duì)分子質(zhì)量在86.18～88.83 kDa之間，理論等電點(diǎn)在5.89～8.07之間，全部為不穩(wěn)定蛋白質(zhì) （不穩(wěn)定指數(shù)＞40），見表3。利用ProtScale模塊對(duì)其進(jìn)行疏水性/親水性分析，可見4條中藥的2，3-氧化鯊烯環(huán)化酶序列具有較為一致的親水性/疏水性整體規(guī)律，但親水區(qū)域和疏水區(qū)域的強(qiáng)弱程度略有不同，見圖7。

圖7 2，3-氧化鯊烯環(huán)化酶疏水性/親水性分析Fig.7 Hydroponic/hydrophilic analysis of 2，3-oxidosqualene cyclases

通過SWISS-MODEL進(jìn)行三維結(jié)構(gòu)的預(yù)測(cè)，見圖8。由圖可知，2，3-氧化鯊烯環(huán)化酶結(jié)構(gòu)穩(wěn)定，與模式植物擬南芥相比，無規(guī)則卷曲部分存在一定程度變化，而且其多分布于蛋白質(zhì)表層。但由于序列有限，構(gòu)建的系統(tǒng)進(jìn)化樹并不能很好地反應(yīng)該酶的系統(tǒng)進(jìn)化規(guī)律，見圖9。

圖8 2，3-氧化鯊烯環(huán)化酶三維結(jié)構(gòu)分析Fig.8 3D structure analysis of 2，3-oxidosqualene cyclases

圖9 2，3-氧化鯊烯環(huán)化酶系統(tǒng)進(jìn)化樹 （NJ法）Fig.9 Phylogenetic tree of 2，3-oxidosqualene cyclases

4 討論

中藥中所含的三萜皂苷成分類型眾多，而且結(jié)構(gòu)復(fù)雜。通過對(duì)其合成通路中關(guān)鍵酶的系統(tǒng)分析發(fā)現(xiàn)，鯊烯合成酶、鯊烯環(huán)氧酶、2，3-氧化鯊烯環(huán)化酶均較為保守，說明雖然三萜皂苷的合成是一個(gè)高度復(fù)雜的動(dòng)態(tài)過程，但其苷元碳環(huán)骨架構(gòu)成等上游步驟在進(jìn)化中趨于保守和穩(wěn)定，而且該成分眾多類型的形成可能與下游細(xì)胞色素P450、糖基轉(zhuǎn)移酶、β-糖苷酶等的修飾有關(guān)。另外，不同類型細(xì)胞色素P450催化的羥化反應(yīng)和糖基轉(zhuǎn)移酶催化的糖基反應(yīng)可能是物種特異性三萜皂苷形成的原因。但由于細(xì)胞色素P450和糖基轉(zhuǎn)移酶均為超家族［17-19］，序列的同源性低，專一性高，大多尚未驗(yàn)證在三萜皂苷通路中的酶活性，難以進(jìn)行系統(tǒng)進(jìn)化分析，是目前通路水平三萜皂苷研究的難點(diǎn)［20］。除了應(yīng)用傳統(tǒng)的基因克隆方法外，還可以采用高通量測(cè)序的手段，對(duì)三萜皂苷合成通路及關(guān)鍵酶進(jìn)行全面闡釋。例如在人參的轉(zhuǎn)錄組分析中，共發(fā)現(xiàn)了133個(gè)細(xì)胞色素P450和235個(gè)糖基轉(zhuǎn)移酶，其中有6個(gè)UDP-糖基轉(zhuǎn)移酶可能直接參與人參皂苷的合成［20］；在西洋參的轉(zhuǎn)錄組中發(fā)現(xiàn)了150個(gè)細(xì)胞色素P450和235個(gè)糖基轉(zhuǎn)移酶，其中有1個(gè)細(xì)胞色素P450和4個(gè)UDP-糖基轉(zhuǎn)移酶可能直接參與皂苷的合成［21］。

雖然很多三萜苷元碳環(huán)骨架上復(fù)雜的修飾反應(yīng)尚不清楚，但是通過對(duì)通路上游關(guān)鍵酶的分析，已經(jīng)具備了人工合成三萜苷元碳環(huán)骨架的理論基礎(chǔ)，增加了對(duì)三萜皂苷合成通路進(jìn)化規(guī)律的認(rèn)識(shí)。鯊烯合成酶所具有的2個(gè)高度保守區(qū)域可作為今后認(rèn)定該類酶的重要依據(jù)，并發(fā)掘更多的三萜皂苷藥用資源。下一步研究的重點(diǎn)將是考察通路中不同基因的表達(dá)模式［23］，并且結(jié)合轉(zhuǎn)錄組和代謝組，高通量地篩選藥用植物中的細(xì)胞色素P450和UDP-糖基轉(zhuǎn)移酶，系統(tǒng)地比較不同修飾體的活性差異，探究其修飾規(guī)律，推進(jìn)人工合成三萜皂苷的進(jìn)程。

［1］明乾良，韓 婷，黃 芳，等.人參皂苷生物合成途徑及其相關(guān)酶的研究進(jìn)展［J］.中草藥，2010，41（11）：1913-1917.

［2］邢增濤，郁瓊花，張勁松，等.不同品種靈芝中三萜類化合物的比較研究［J］.中藥材，2004，27（8）：575-576.

［3］吳耀生，朱 華，李 珅，等.三七鯊烯合酶基因在三七根，莖，蘆頭中的轉(zhuǎn)錄表達(dá)與三萜皂苷合成［J］.中國生物化學(xué)與分子生物學(xué)報(bào)，2007，23（12）：1000-1005.

［4］程 敏，胡正海.墨旱蓮的生物學(xué)和化學(xué)成分研究進(jìn)展［J］.中草藥，2010，41（12）：2116-2118.

［5］Choi DW，Jung JD，Ha Y I，et al.Analysis of transcripts in methyl jasmonate-treated ginseng hairy roots to identify genes involved in the biosynthesis of ginsenosides and other secondary metabolites［J］.Plant Cell Rep，2005，23（8）：557-566.

［6］焦曉林，高微微.環(huán)境因子對(duì)藥用植物三萜皂苷合成影響的研究進(jìn)展［J］.中草藥，2011，42（2）：398-402.

［7］Marchler-Bauer A，Lu S，Anderson JB，et al.CDD：a conserved domain database for the functional annotation of proteins［J］.Nucleic Acids Res，2011，39：D225-229.

［8］Letunic I，Doerks T，Bork P.SMART 7：recent updates to the protein domain annotation resource［J］.Nucleic Acids Res，2012，40：D302-305.

［9］Crooks GE，Hon G，Chandonia JM，et al.WebLogo：a sequencelogo generator［J］.Genome Res，2004，14（6）：1188-1190.

［10］Kiefer F，Arnold K，KünzliM，et al.The SWISS-MODEL Repository and associated resources［J］.Nucleic Acids Res，2009，37：D387-392.

［11］Tamura K，Peterson D，Peterson N，et al.MEGA5：molecular evolutionary genetics analysis usingmaximum likelihood，evolutionary distance，andmaximum parsimonymethods［J］.Mol Biol Evol，2011，28（10）：2731-2739.

［12］張風(fēng)俠，梁新華，王 俊.植物三萜皂苷生物合成及關(guān)鍵酶鯊烯合酶的研究進(jìn)展［J］.農(nóng)業(yè)科學(xué)研究，2009，30（3）：64-68.

［13］Nakashima T，Inoue T，Oka A，et al.Cloning，expression，and characterization of cDNAs encoding Arabidopsis thaliana squalene synthase［J］.Proc Natl Acad Sci U S A，1995，92（6）：2328-2332.

［14］Tansey TR，Shechter I.Structure and regulation ofmammalian squalene synthase［J］.Biochim Biophys Acta，2000，1529（1-3）：49-62.

［15］Lee JH，Yoon Y H，Kim H Y，etal.Cloning and expression of squalene synthase cDNA from hot pepper（Capsicum annuum L.）［J］.Mol Cells，2002，13（3）：436-443.

［16］Phillips D R，Rasbery JM，Bartel B，et al.Biosynthetic diversity in plant triterpene cyclization［J］.Curr Opin Plant Biol，2006，9（3）：305-314.

［17］Werck-Reichhart D，F(xiàn)eyereisen R.Cytochromes P450：a success story［J］.Genome Biol，2000，1（6）：1-9.

［18］Ross J，LiY，Lim E K，et al.Higher plantglycosyltransferases［J］.Genome Biol，2001，2（2）：1-6.

［19］徐潔森，魏 建，陶韻文，等.植物細(xì)胞色素P450在三萜皂苷生物合成中的功能研究進(jìn)展［J］.中草藥，2012，43（8）：1635-1640.

［20］鄭夏生，羅秀秀，徐 暉，等.崗梅轉(zhuǎn)錄組及其烏索烷型三萜皂苷生物合成相關(guān)酶基因的發(fā)掘［J］.世界科學(xué)技術(shù)-中醫(yī)藥現(xiàn)代化，2014，16（7）：1505-1512.

［21］Chen S，Luo H，Li Y，et al.454 EST analysis detects genes putatively involved in ginsenoside biosynthesis in Panax ginseng［J］.Plant Cell Rep，2011，30（9）：1593-1601.

［22］Sun C，Li Y，Wu Q，et al.De novo sequencing and analysis of the American ginseng root transcriptome using a GSFLX Titanium platform to discover putative genes involved in ginsenoside biosynthesis［J］.BMCGenomics，2010，11：262-273.

［23］牛云云，朱孝軒，羅紅梅，等.三萜皂苷合成生物學(xué)元件的初步開發(fā)：三七鯊烯環(huán)氧酶編碼基因克隆及表達(dá)模式分析［J］.藥學(xué)學(xué)報(bào)，2013，48（2）：211-218.

Bioinformatics analysis of triterpenoid saponins biosynthesis pathway in herbalmedicines

ZHANG Zhao-bao1， HOU Lin1， CUIQing-hua1， PAN Qing1， MA Lu-yu2*
（1.College of Pharmacy，Shandong University of Traditional ChineseMedicine，Jinan 250300，China；2.Shandong Academy ofMedical Sciences，Jinan 250062，China）

AIMTo conductan informative analysis of key enzyme in biosynthesis pathway of triterpenoid saponins from herbalmedicines in the perspective of phylogenetics.METHODSFull protein sequences of squalene synthase，squalene epoxidase and 2，3-oxidosqualene cyclase，belonging to seventeen kinds of herbs，were downloaded from NCBIdatabase.Afterwards，analysis of physicaland chemical propertieswas conducted，conserved domain analysis，3D structures and Neighbor-Joining（NJ）phylogenetic analysis were carried out by bioinformatic softwares.RESULTSIt is shown that all the three enzymes were highly conservative.For squalene synthases，there were two completely conserved domains，Mg2+binding site with substrates，which were aspartic-rich and could be used as sequence tags.For squalene epoxidases and 2，3-oxidosqualene cyclases，their3D structureswere also highly conservative with few flexible random coils.CONCLUSIONThe upstream syntheses are highly conserved in evolution，the formation ofmany types of triperpenoid saponins may be associated with modification of downstream cytochrome P450，glycosyltranferases and beta-glucosidases.

triterpenoid saponin；squalene synthase；squalene epoxidase；2，3-oxidosqualene cyclase

R284.1

：A

：1001-1528（2015）06-1255-07

10.3969/j.issn.1001-1528.2015.06.021

2014-08-16

國家自然科學(xué)基金項(xiàng)目 （81273394）；“十二五”重大新藥創(chuàng)制國家科技重大專項(xiàng) （2013ZX09103002-005）

張召寶（1988—），男，碩士，研究方向?yàn)橹兴帉W(xué)。Tel：（0531）82919863，E-mail：zhangzhaobao2010@163.com

*通信作者：馬魯豫 （1972—），男，博士，副研究員，研究方向?yàn)橹兴幮滤庨_發(fā)。Tel：（0531）82919863，E-mail：maluyuyky@ 126.com