克雷伯氏菌ZD112氯氰菊酯降解酶基因的克隆与生物信息学分析.docx

克雷伯氏菌ZD112氯氰菊酯降解酶基因的克隆与生物信息学分析.docx

《克雷伯氏菌ZD112氯氰菊酯降解酶基因的克隆与生物信息学分析.docx》由会员分享，可在线阅读，更多相关《克雷伯氏菌ZD112氯氰菊酯降解酶基因的克隆与生物信息学分析.docx（5页珍藏版）》请在冰豆网上搜索。

克雷伯氏菌ZD112氯氰菊酯降解酶基因的克隆与生物信息学分析

克雷伯氏菌ZD112氯氰菊酯降解酶基因的克隆与生物信息学分析

（作者:

___________单位:

___________邮编:

___________）

【摘要】目的筛选新的菊酯农药降解酶基因并了解其特性，为解决食品和环境中的菊酯农药残留问题创造条件。

方法通过构建基因文库的方式筛选获得新的菊酯农药降解酶（EstP）基因，并对其进行一系列生物信息学分析。

结果获得了新的菊酯农药降解酶（EstP）基因。

EstP属于水解酶类，由638个氨基酸编码，预测分子量为73.4kDa，pI值为8.34，与其他蛋白相似性极低，不属于任何已知的蛋白超家族，仅与羧肽酶Taq（M32）的保守域具有一定相似性，推测Lys137,Glu116,Glu148为其必需氨基酸。

结论新的菊酯农药降解基因得到克隆，生物信息学分析为其定向进化提供了理论指导，为高效基因工程菌的构建打下了基础。

【关键词】菊酯农药降解酶；基因克隆；菊酯农药；生物信息学分析　　MolecularcloningandbioinformaticanalysisofanovelpyrethroidhydrolyzingesterasefromKlebsiellasp.StrainZD112　　WANGZhuoya,LIUYuhuan,LIHe　　1.GuangdongPharmaceuticalUniversity,Guangzhou,Guangdong510006;　　2.SUNYatsenUniversity,Guangzhou,Guangdong510275,China　　Abstract:

ObjectiveToobtainpyrethroidhydrolyzingesterase（EstP）geneandanalyzethecharacteristicsofit.MethodsBioinformaticsanalyzingtoolswereusedtoscreenaEstPgenefromKlebsiellasp.StrainZD112DNAlibraryandpredictthecharacterandfunctionoftheprotein.ResultsAnovelEstPwithhigheractivitywasfound.EstPwasaproteinof638aminoacidresidues,amolecularmassof73.4kDa,andpI8.34.Itshowednosimilaritieswithknownproteins,anddidnotbelongtoanyproteinsuperfamily.ItisexpectedtoberelatedtocarboxypeptidaseTaq（M32）metallopeptidase,andrevealedthatLys136,Glu114,Glu147wereimportantfortheesteraseactivity.ConclusionAnewEstPgenewasobtained,whichmightlayabasisforthefurtherstudyonthestructureanalysisandgenecloneofit.　　Keywords:

pyrethroidhydrolyzingesterase;molecularcloning;pyrethroidinsecticides;bioinformaticanalysis　　拟除虫菊酯为第三代杀虫剂，其长期广泛的使用，对土壤、大气和水源等造成了不同程度的污染［1］。

其对健康的危害也不容忽视，主要包括对免疫系统的胁迫，对淋巴结和脾脏的损伤、致癌以及环境激素效应等［2］。

寻找一种有效方法消除或降低食品和环境中的菊酯农药残留已成为当前迫切需要解决的问题。

　　　　在过去的研究工作中，本实验室获得了一株能以菊酯为唯一碳源和氮源生长的伯雷克氏菌株Klebsiellasp.ZD112，鉴于直接应用农药降解酶制剂比应用产酶的菌剂更有优势［3］，为了能够产业化生产菊酯农药降解酶，通过构建基因文库的方式获得了新的菊酯农药降解酶基因，并进行了一系列生物信息学分析，为高效基因工程菌的构建打下基础。

　　1材料与方法　　1.1培养基　　　　LB培养基与LA培养基按照分子克隆手册配制。

筛选培养基为加入100mg/L氨苄青霉素（ampicillin）和100μmol/L5bromo4chloro3indolylcaprylate（Xcaprylate）的LB培养基（Xcaprylate可被水解酶水解而形成蓝色的水解圈［4］）。

　　1.2菌株　　大肠杆菌DH5α［supE44,ΔlacU169（φ80lacZΔM15）,hsdR17,recA1,endA1,gyrA96,thi1,relA1］由本实验室保存，菊酯农药降解酶生产菌株Klebsiellasp.ZD112菌株由本实验室从污染土壤分离。

pUC18质粒购自TaKaRa公司。

　　1.3酶和试剂　　各种限制性内切酶、DNAmarker、XGal、IPTG购自TaKaRa公司,T4连接酶、去磷酸化酶购自MBI公司。

其他生物化学试剂均为国产分析纯。

　　1.4基因文库的构建　　参照分子克隆手册提取Klebsiellasp.ZD112基因组DNA和pUC18质粒。

ZD112基因组DNA经HindⅢ不完全酶切后电泳，按照3SDNAGelPurificationkit试剂盒说明书割胶回收2～10kb片断。

pUC18经HindⅢ酶切后去磷酸化后并电泳，割胶回收。

将酶切并回收的DNA片断与去磷酸化的克隆载体按5∶1的比例16℃连接过夜，Ca转入大肠杆菌DH5α，并涂布LA平板，进行蓝白斑筛选。

　　1.5目标亚克隆的筛选　　挑取白色的转化子到含酯酶底物（Xcaprylate）的筛选平板上培养过夜，挑取有蓝色水解圈形成的转化子保存并接种于LA液体培养基中振荡培养。

提取质粒酶切验证，同时破碎菌体提取粗酶液,测定其对氯氰菊酯的降解率。

　　1.6亚克隆测序与分析　　测序由上海博亚生物公司完成。

对插入片断通过ORFFinder（http:

//www.ncbi.nlm.nih.gov/gorf.html）寻找读码框，并采用BLASTn（http:

//www.ncbi.nlm.nih.gov/BLAST/）进行同源性分析。

　　1.7目的基因生物信息学分析　　分子量和等电点计算采用DNAstar5.0软件计算；蛋白保守域查询采用rpsBlast（http:

//www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi），数据库选择为CDDv2.06-11530PSMMs，鉴于EstP与其他蛋白的同源性非常低，E值选择时适当放宽，设为1；酶与非酶分析采用基于蛋白氨基酸序列性质的ArchaeaFun1.0Server（http:

//genome.cbs.dtu.dk/services/ArchaeaFun）；蛋白超家族分析采用基于隐马尔科夫（HiddenMarkovModel,HMM）算法的Superfamily1.61（http:

//supfam.mrclmb.cam.ac.uk/SUPERFAMILY/），参数选择：

examplesequence:

Aminoacidsequence，Scoringoptions:

Local/localmodel/sequencescoring，BlastPrefilterPvalue10e10；蛋白亚细胞定位采用Loctree服务器（http:

//cubic.bioc.columbia.edu/services/loctree/）；EstP的多序列比对使用ClustalX软件，相似序列由rpsBlast获得，蛋白势能矩阵选择为Blosumseries。

　　2结果与分析　　2.1基因文库的构建　　通过蓝白斑筛选，共获得白色克隆约3万个。

ZD112基因组DNA不完全酶切结果见图1，质粒pUC18的HindⅢ酶切图见图2。

　　2.2目标亚克隆的筛选与鉴定　　挑取1500个阳性克隆至到含底物（Xcaprylate）筛选培养基。

获得4个能在筛选培养基产生蓝色水解圈的单菌落，挑取后加入Amp含量为100mg/L的LB液体培养基培养,甘油保存。

其粗酶液对氯氰菊酯的降解率见表1，ZD112粗酶液作对照［5］。

　　挑选降解率高的阳性克隆子LA5进行酶切（图3），从图中可知，L5A连接上了大约3.6kb左右的片段,送测序公司测序（图4）。

用NCBI上的工具ORFfinder分析LA5的克隆片段,发现其含有2个较大的开放读码框（ORF）,一个长度为1914bp（ORF1，登录号:

AY995176）,另外一个长度为756bp（ORF2，登录号:

AY954696）。

其中ORF1仅与2个细菌预测蛋白有一定相似性（hypotheticalproteinSC171：

88%相似；ebA4602：

36%相似），并含有-10,-35以及RBS（ribosomalbindingsite）序列（图5），因此该ORF可能是编码菊酯降解酶（EstP）的新基因。

表1各阳性克隆与ZD112的降解率（略）　　2.3生物信息学分析结果　　　　EstP由637个氨基酸组成，经DNAstar5.0计算出其分子量为73.4kDa，pI值为8.34。

蛋白亚细胞定位分析结果显示，EstP属于细胞质内蛋白。

酶与非酶的分析结果表明，EstP是酶，并且属于水解酶类（图6）。

　　　进一步分析发现，EstP不属于任何已知的蛋白超家族，与EstP具同源性的已知蛋白也非常少。

rpsBlast分析（图7）表明，EstP仅有一段很短的序列（71～149氨基酸区段）与羧肽酶Taq（M32）的保守域有最高的相似性，这段序列还与预测膜蛋白COGO0012有一定的相似性。

通过对已知结构蛋白Taq（M32）的三级结构（图8）的分析发现，其与EstP具有同源性的序列呈V字型，且Taq（M32）的活性位点正位于这一区段内，由此，推测EstP的活性中心可能位于其71～149氨基酸区段内。

　　3讨论　　　　作为目前最有发展前途的杀虫剂，拟除虫菊酯尚未找到新型的替代化合物，在未来十几年内仍将广泛使用。

采用微生物农药降解酶解决菊酯农药残留污染问题是当前环境科学研究的热点,也是一个较新的研究领域。

由于农药降解酶在野生菌株中含量太低，生产成本高昂且难以大量生产，如何产业化地生产农药降解酶成为农药降解酶应用必须解决一个关键性问题。

农药降解酶基因的克隆、表达以及定向进化可以构建出降解谱广、降解能力强的工程菌，为微生物降解农药开辟了新途径。

　　　　在之前的研究中，本实验室从农药厂的污泥中筛选出能降解菊酯农药的伯雷克氏菌株Klebsiellasp.StrainZD112。

通过构建基因文库的方式，获得了一个具有菊酯降解活性的克隆子，对测序结果进行ORF分析，该ORF1在NCBI上比对的结果都是一些细菌的预测蛋白，如：

它与Salmonellaentericasubsp.entericaserovarCholeraesuis的预测蛋白hypotheticalproteinSC171（登录号:

AAS76447.1）有88%一致性；与Azoarcussp.EbN1的预测蛋白ebA4602（登录号:

YP159633.1）有36%的一致性。

他们的功能都还没有确定,目前为止尚未见有菊酯类农药降解基因的报道。

因此本研究克隆出来的基因有可能是一个新的菊酯农药降解基因。

　　　　对该基因的生物信息学分析显示，EstP为水解酶且为胞内蛋白，这与前期实验结论完全吻合［5］。

根据其与羧肽酶Taq（M32）保守域所表现的较高相似性，推测其氨基酸序列的71～149为其活性中心所在，必需氨基酸为Glu116,Lys137，Glu148，这一结论的得出为EstP的定向进化提供了理论指导。

然而，EstP与其他已知蛋白序列相似性极低，且不属于任何蛋白超级家族，这增加了对其研究的难度。

　　　　对这样一个全新水解酶的进一步研究，一方面，可深入了解菊酯降解酶的作用机制［6］，具有深刻的理论意义；另一方面，高效工程菌的开发有助于菊酯农药残留降解，在保护生态环境、突破绿色壁垒、维护人民健康等方面具有很大的应用潜力。

【参考文献】　　［1］MIYAMOTOJ,KANEKOH,TSUJIR,etal.Pyrethroids,nervepoisons:

howtheirriskstohumanhealthshouldbeassessed［J］.ToxicolLett,1995,83:

933-940.　　［2］PARKEK,KIMJH,GEESJ,etal.Determinationofpyrethroidresiduesinagriculturalproductsbyanenzymelinkedimmunosorbentassay［J］.JAgricFoodChem,2004,52:

5572-5576.　　［3］虞云龙，樊德方，陈鹤鑫．农药微生物降解的研究现状与发展策略［J］．环境科学进展，1996，4：

28-34．　　［4］CHOIYJ,MIGUEZCB,LEEBH.CharacterizationandheterologousgeneexpressionofanovelesterasefromLactobacilluscaseiCL96［J］.ApplEnvironMicrobiol,2004,70:

3213-3221.　　［5］梁卫驱,刘玉焕，李荷.氯氰菊酯降解菌的分离鉴定及其降解特性研究［J］.广东药学院学报,2007,23

（2）：

199.　　［6］皮埃尔·巴迪尔，索恩·布鲁纳克．生物信息学－机械学习方法［M］．北京：

中信出版社，2003：

211-219．