中药三萜皂苷的生物信息学分析范文

《中成药杂志》2015年第六期

1方法

1.1蛋白质理化性质分析蛋白质全长序列的理化性质应用在线软件ExPASy分析。其中，利用ProtParam模块根据具体的氨基酸序列，进行氨基酸数目、相对分子质量、理论等电点和不稳定指数的测算，然后利用ProtScale模块进行亲水性/疏水性分析。

1.2结构域分析在NCBI的ConservedDomainDatabase和SMART数据库中进行蛋白质的保守域和功能域预测，并通过WebLogo［9］对其做进一步分析。利用基于同源建模的分析工具SWISS-MODEL预测蛋白质的三维结构。

1.3构建系统进化树利用MEGA5.0软件［11］中的ClustalW进行序列比对，然后利用Neighbor-Joining(NJ)法构建系统进化树(bootstrap1000)。

2结果与分析

2.1鲨烯合成酶鲨烯合成酶(squalenesyn-thase，EC2.5.1.21)是一种定位于内质网上的膜结合蛋白［12］。在类异戊二烯代谢途径中，它以法尼二磷酸为底物，催化合成30C产物角鲨烯。由于处于甾醇与三萜合成通路的分支点，决定着法尼基焦磷酸的流向，因此鲨烯合成酶发挥着重要的调节作用。已有研究表明，鲨烯合成酶的结构和反应机制在动植物中相当保守。在哺乳动物中，该成分大约包含416aa，相对分子质量约为47kDa;在模式植物拟南芥中，该成分有2个(NP_195190.1和NP_195191.2)，其外显子具有高度同源性［13］。本研究从NCBI数据库中共下载得到17种中药的24条蛋白质全长序列，长度在405～418aa之间，相对分子质量在45.89～47.89kDa之间，理论等电点在6.07～8.53之间，多数为不稳定蛋白质(不稳定指数＞40)。以人参中的三条鲨烯合成酶序列(BAD08242.1，ACV88718.1和ACZ71037.1)为例，利用ProtScale模块进行疏水性/亲水性分析，由于其序列高度保守，所展现出的亲水性/疏水性规律较为一致，见图1。另外由图可知，鲨烯合成酶的N端分布有亲水域(分值＜－0.5)，而C端则分布有一段较短的强疏水域(分值＞0.5)，这与之前有关鲨烯合成酶细胞定位的报道相一致。鲨烯合成酶通过短的C端跨膜区锚定于内质网上，而N端的催化域则伸向胞质溶胶中。在对其序列进行理化性质分析的基础上，进一步研究其保守域。首先采用WebLogo对Lee［15］等报道的6个保守域进行分析，结果发现在17种中药中，Ⅲ、Ⅳ、Ⅴ三个保守域高度保守，而Ⅰ、Ⅱ、Ⅵ三个保守域保守性较差，见图2A。由此可见，之前对这6个保守域的报道具有一定局限性。根据NCBI的ConservedDomainDatabase的预测结果比对发现，这17种中药中的鲨烯合成酶具有两个高度保守(100%)的区域———富含天冬氨酸的镁离子结合位点1和2，见图2B(结合位点1位于77～81aa，结合位点2位于213～217aa)，该区域与15C法尼基和20C香叶基异戊二烯基二磷酸的1'－1'位缩合反应相关，其保守性保证了鲨烯合成酶与底物特异性结合，并可作为今后鉴定中药中该成分的重要依据。以人参中的三条鲨烯合成酶序列(BAD08242.1，ACV88718.1和ACZ71037.1)为例，SWISS-MODEL建模的结果显示，人参中鲨烯合成酶的三维结构相对保守，以α螺旋为主，而β折叠片较少，无规则卷曲等结构均很保守。以模式物种拟南芥等植物中的鲨烯合成酶蛋白质全长构建进化树，见图3。由图可知，进化树首先分为两大支，其中上面一支主要包含较高等的植物类群，主要为双子叶植物;下面一支以紫杉Taxuschinensis(Pilger)Rehd．、马尾松Pinusmasso-niana为代表，多为较低等的植物类群。整体而言，鲨烯合成酶具有较高的保守性和较为普遍的基因倍增现象，例如人参、丹参、北柴胡等6种中药具有两个及两个以上鲨烯合成酶。由于其具有较为重要的生理生化功能，故认为适应性进化可能导致了普遍的基因倍增，从而产生了不同的亚型，之间可能存在的功能上的互补和增强，又是酶功能进化的一种体现。

2.2鲨烯环氧酶鲨烯环氧酶(squaleneepoxi-dase，EC1.14.13.132)是一种细胞膜结合蛋白，催化三萜合成的第一步氧化反应，是该合成通路的限速酶之一。以三七和人参中的鲨烯环氧酶序列为例，利用ProtScale模块进行疏水性/亲水性分析，可以看出，三七和人参中的两条鲨烯环氧酶蛋白质全长序列出现了较为明显的分化;三七中的ABE73759.1与人参的BAA24448.1具有较为一致的亲水性/疏水性规律;三七中的AFV92748.1与人参的ACJ24907.2具有较为一致的亲水性/疏水性规律;这说明这两种亚型之间的分化早于三七和人参之间的分化，并在随后的进化中形成了自己独特的亲水性/疏水性规律。对其进行功能域分析，见图5A。由图可知，鲨烯环氧酶的C端含有一个高度保守的跨膜结构域(5～27aa)，然后是1个重叠的FAD结合域(NADPH结合域，60～408aa)和鲨烯环氧酶结构域(pfam08491，211～484aa)。药用植物中的鲨烯环氧酶蛋白质序列呈现出较高的同源性，以模式植物拟南芥和中药代表物种人参为例，通过SWISS-MODEL进行三维结构的预测，见图5B。由图可知，鲨烯环氧酶的α螺旋和β折叠片等结构保守，无规则卷曲部分为可变区，其中连接β折叠片者较为高变。以模式物种拟南芥等植物中的鲨烯环氧酶蛋白质全长构建进化树，见图6。由图可知，进化地位较为低等的蕨类植物———蛇足石杉作为外群，位于进化树的底部;五加科的人参、三七、刺五加等植物聚为一支，显示出较近的亲缘关系。该进化树较好地符合物种之间的进化规律，说明鲨烯环氧酶伴随物种的进化进程趋于保守，只在特定类群中出现了基因倍增的现象，如五加科，这可能与该科物种特定的生境与生理习性相关，例如，某种环境压力迫使其特定通路的完善或者特定代谢产物的富集。

2.32，3-氧化鲨烯环化酶氧化鲨烯环化酶(ox-idosqualenecyclases，EC5.4.99.7)是一大类氧化鲨烯环化酶的统称，能催化一系列复杂的环化和重排反应。在植物中，由于它可催化三萜皂苷合成中的第一步环化反应，因而被认为是合成通路中的关键酶。本研究对其中的2，3-氧化鲨烯环化酶进行深入的分析，发现它在真菌、动物和植物中均较为保守。从NCBI中共下载得到4条中药的2，3-氧化鲨烯环化酶全长蛋白序列，包括人参、百脉根、蒲公英、总状升麻，长度在757～780aa之间，相对分子质量在86.18～88.83kDa之间，理论等电点在5.89～8.07之间，全部为不稳定蛋白质(不稳定指数＞40)，见表3。利用ProtScale模块对其进行疏水性/亲水性分析，可见4条中药的2，3-氧化鲨烯环化酶序列具有较为一致的亲水性/疏水性整体规律，但亲水区域和疏水区域的强弱程度略有不同，见图7。通过SWISS-MODEL进行三维结构的预测，见图8。由图可知，2，3-氧化鲨烯环化酶结构稳定，与模式植物拟南芥相比，无规则卷曲部分存在一定程度变化，而且其多分布于蛋白质表层。但由于序列有限，构建的系统进化树并不能很好地反应该酶的系统进化规律，见图9。

3讨论

中药中所含的三萜皂苷成分类型众多，而且结构复杂。通过对其合成通路中关键酶的系统分析发现，鲨烯合成酶、鲨烯环氧酶、2，3-氧化鲨烯环化酶均较为保守，说明虽然三萜皂苷的合成是一个高度复杂的动态过程，但其苷元碳环骨架构成等上游步骤在进化中趋于保守和稳定，而且该成分众多类型的形成可能与下游细胞色素P450、糖基转移酶、β-糖苷酶等的修饰有关。另外，不同类型细胞色素P450催化的羟化反应和糖基转移酶催化的糖基反应可能是物种特异性三萜皂苷形成的原因。但由于细胞色素P450和糖基转移酶均为超家族［17-19］，序列的同源性低，专一性高，大多尚未验证在三萜皂苷通路中的酶活性，难以进行系统进化分析，是目前通路水平三萜皂苷研究的难点［20］。除了应用传统的基因克隆方法外，还可以采用高通量测序的手段，对三萜皂苷合成通路及关键酶进行全面阐释。例如在人参的转录组分析中，共发现了133个细胞色素P450和235个糖基转移酶，其中有6个UDP-糖基转移酶可能直接参与人参皂苷的合成［20］;在西洋参的转录组中发现了150个细胞色素P450和235个糖基转移酶，其中有1个细胞色素P450和4个UDP-糖基转移酶可能直接参与皂苷的合成［21］。虽然很多三萜苷元碳环骨架上复杂的修饰反应尚不清楚，但是通过对通路上游关键酶的分析，已经具备了人工合成三萜苷元碳环骨架的理论基础，增加了对三萜皂苷合成通路进化规律的认识。鲨烯合成酶所具有的2个高度保守区域可作为今后认定该类酶的重要依据，并发掘更多的三萜皂苷药用资源。下一步研究的重点将是考察通路中不同基因的表达模式［23］，并且结合转录组和代谢组，高通量地筛选药用植物中的细胞色素P450和UDP-糖基转移酶，系统地比较不同修饰体的活性差异，探究其修饰规律，推进人工合成三萜皂苷的进程。

作者：张召宝 侯林 崔清华 潘晴 马鲁豫 单位：山东中医药大学药学院 山东省医学科学院