小菜蛾碱性磷酸酯酶受体表达与分子模拟

张 霄，胡晓丹，仲建锋，武爱华，徐重新，刘 媛，张存政，谢雅晶，刘贤金



小菜蛾碱性磷酸酯酶受体表达与分子模拟

张 霄，胡晓丹，仲建锋，武爱华，徐重新，刘 媛，张存政，谢雅晶，刘贤金

（江苏省农业科学院食品质量安全与检测研究所/江苏省食品质量安全重点实验室-省部共建国家重点实验室培育基地/农业部农产品质量安全 控制技术与标准重点实验室，南京 210014）

【目的】利用原核表达小菜蛾()中肠膜结合碱性磷酸酯酶（membrane-bound alkaline phosphatase, mALP）并经Ligand blot验证其具有与Cry1Ac毒素结合的能力；通过同源建模和分子对接研究Cry1Ac-mALP的结合模式，预测毒素和受体结合区域及关键氨基酸位点（热点残基），为了解毒素-受体互作机制及分子改造增强Cry毒素活性的研究打下基础。【方法】针对小菜蛾mALP全长设计引物，并以小菜蛾cDNA为模板扩增mALP基因，双酶切后用T4连接酶连接至pET-26b原核表达载体，将构建的pET-26b-mALP载体转化Trans1-T1克隆感受态，挑取克隆并提取质粒后进行PCR、双酶切和测序验证，将验证无误的重组质粒转化BL21（DE3）表达感受态细胞，进行诱导表达。将诱导表达后的mALP转至PVDF膜上，通过Western blot和Ligand blot分别验证mALP是否成功表达以及是否具有与Cry1Ac毒素结合的能力。对mALP进行同源建模、分子动力学模拟以及模型评价，获得的mALP最佳三维结构与Cry1Ac毒素利用PatchDOCK和FireDock程序进行分子对接试验，对确定的最佳毒素-受体复合物进行结合区域和结合氨基酸位点分析，并通过计算机辅助的丙氨酸突变扫描试验确定毒素和受体参与的关键氨基酸残基。【结果】 扩增出小菜蛾mALP基因并克隆至pET-26b原核表达载体，转化BL21（DE3）表达感受态后挑取阳性克隆提取质粒后进行PCR、双酶切和测序均显示构建正确。通过原核表达和Western blot验证成功表达了mALP蛋白，并经Ligand blot试验证实了原核表达的mALP具有和Cry1Ac毒素结合的能力。利用同源建模成功获得了mALP的三维结构，通过PatchDOCK和FireDock分子对接程序，获得毒素和受体的对接复合物，通过溶剂可及表面积变化计算和Ligplot分析，确定毒素结构域II和结构域III均参与了受体结合，并且毒素和受体均以疏水结合和氢键结合模式参与结合，最后通过热点残基预测发现Cry1Ac毒素和mALP中分别有3个氨基酸残基（376Asn、443Ser和486Ser）和4个氨基酸残基（452Arg、499Thr、502Tyr和513Tyr）是参与互作的关键氨基酸位点。【结论】经原核表达的小菜蛾mALP同样具有与Cry1Ac毒素结合的能力，并利用分子模拟技术预测了小菜蛾mALP三维结构及与Cry1Ac毒素结合模式。

碱性磷酸酯酶；原核表达；配体印迹；同源建模; 分子对接；热点残基预测

0 引言

【研究意义】苏云金芽孢杆菌（，Bt）可在其生命周期中产生一种专一性极强的晶体毒素（crystal toxin，Cry toxin），导致昆虫肠道细胞溶解或引发细胞凋亡程序，从而杀死昆虫[1-3]。Cry1Ac毒素作为应用最为广泛的一类微生物杀虫蛋白，在防治大量田间害虫中已取得巨大成效[4-5]。然而，Cry毒素的抗药性问题已成为世界级难题，尤其是小菜蛾（）被认为是田间抗药性最严重和最难防治的世界性害虫之一[6-8]。目前认为受体的变异或缺失导致毒素结合能力的改变是抗性产生的最主要原因[9-11]。因此研究毒素-受体结合模式对改善抗性问题显得尤为重要，分子模拟正是针对这一问题而提出的理性解决方法[12-13]。【前人研究进展】碱性磷酸酯酶（alkaline phosphatase，ALP）作为Cry毒素在害虫体内的靶标受体之一[1-3]，已有很多研究显示ALP受体与害虫抗性密切相关。Jurat-Fuentes等[14]比较分析了敏感性和抗性烟芽夜蛾（）幼虫中ALP与Cry1Ac的结合特征，推测ALP为使其产生抗性的重要因素；Jurat-Fuentes等[15]利用逆转录定量PCR（qRT-PCR）发现多种鳞翅目靶标昆虫中肠膜结合碱性磷酸酯酶（membrane- bound alkaline phosphatase，mALP）的表达量下调是导致其产生抗性的重要因素；Chen等[16]报道了敏感性和抗性棉铃虫ALP受体上影响与Cry毒素结合区域（toxin-binding region，TBR），首次发现棉铃虫ALP受体TBR区域对Cry毒素活性发挥的作用。小菜蛾是田间最早报道对Cry1Ac毒素产生抗性的昆虫，因而其抗性问题也一直是本领域的研究热点[17-19]。Yang等[20]分别对敏感性和抗性小菜蛾的碱性磷酸酯酶以及其他潜在受体分别进行表达和生物学研究，为小菜蛾抗性治理提供了指导；刘洁等[21]克隆小菜蛾抗性、敏感种群中的部分碱性磷酸酶基因，通过与敏感种群相比发现抗性种群中该基因片段发生基因突变，该结果对于进一步研究小菜蛾全基因结构、功能有重要意义；近期Guo等[22]报道并证实了小菜蛾mALP可被MAPK（丝裂原活化蛋白激酶）信号途径反式调控从而导致小菜蛾对Cry毒素的高抗性。在众多研究毒素-受体结合模式的方法中，利用分子模拟（同源建模和分子对接）研究，在已知毒素和受体可以发生结合的情况下，正确预测复合物的结合模式是一种极为有效的途径[23-24]，Shan等[25]利用同源建模构建了棉铃虫的ALP三维模型并通过生物信息学在线网站对ALP受体功能进行了理论研究。【本研究切入点】利用分子模拟研究Cry1Ac和小菜蛾mALP的结合模式还未见报道。【拟解决的关键问题】原核表达具有Cry1Ac毒素结合活性的mALP受体，用于体外进行与毒素模拟物的结合分析，并利用分子模拟研究毒素和受体互作模式以及参与互作的关键氨基酸位点，进一步了解Cry毒素与mALP的互作机制，为利用分子改造提升Cry毒素与受体结合能力以增强Cry毒素活性的研究提供新的研究思路和技术途径。

1 材料与方法

试验于2015年12月至2016年7月在江苏省农业科学院食品质量安全与检测研究所完成。

1.1 材料

1.1.1 供试昆虫、表达载体及菌株 小菜蛾（ 4龄幼虫）由笔者实验室于室内饲养多年，期间未施任何杀虫剂；表达载体pET-26b（+）购于Novagen公司；Trans1-T1、BL21（DE3）感受态细胞购自北京全式金生物技术有限公司。

1.1.2 试验试剂 Trizol和Super script III反转录试剂盒购自Invitrogen公司；Fast Pfu DNA聚合酶购自近岸蛋白质科技有限公司；卡那霉素和IPTG购自Sigma公司；PCR纯化试剂盒购自Promega公司；胶回收试剂和质粒提取试剂盒购自Axygen公司；I、I内切酶、T4 DNA连接酶购自美国NEB公司；HRP标记Anti His-Tag鼠单抗、HRP-羊抗兔IgG购自GE Healthcare公司；蛋白Marker购自Thermo公司；增强型HRP-DAB底物显色试剂盒购自天根生化科技（北京）有限公司；PVDF膜购自BIOSHRP；其他试剂均为分析纯级试剂。

1.2 方法

1.2.1 小菜蛾mALP基因扩增 将小菜蛾4龄幼虫置冰上冷冻，从其尾部拉取中肠，用0.5%的生理盐水清洗干净，放入匀浆器中，利用Trizol法提取总RNA，用反转录试剂盒将RNA反转录为cDNA第一链（按说明书方法）。根据Genbank公布的小菜蛾mALP基因序列（登录号：GenBank KC841472.2）设计全长引物。上游引物mALP-F：5′-CATGCCTCTC GCGTGGCGCGCCAGGTATC-3′（下划线为I酶切位点，CATG：保护碱基，CC：防止移码所加碱基）；下游引物mALP-R：5′-AATTAATAAGC GTCTCAGATACG-3′（下划线为I酶切位点，AAT：保护碱基）。PCR扩增条件：95℃预变性2 min；94℃变性30 s，57℃退火30 s，72℃延伸1 min，共30个循环；最后72℃终延伸10 min，4℃保存。扩增产物经1.0%琼脂糖凝胶电泳后，PCR产物纯化试剂盒纯化目的片段。

1.2.2 重组质粒pET-26b-mALP的构建及鉴定 纯化后的mALP基因片段和提取的pET-26b表达载体，同时进行Ⅰ和Ⅰ限制性内切酶酶切，纯化回收酶切产物，经T4连接酶连接过夜，次日转化至Trans1- T1感受态细胞并涂布于含50 µg·mL-1卡那霉素的LB平板上进行阳性转化子筛选，挑取阳性克隆过夜培养后提取质粒进行PCR鉴定、双酶切鉴定及序列测定（由上海生工生物工程公司完成），并将测序结果提交至GenBank进行序列同源性比对。

1.2.3 pET-26b-mALP在BL21（DE3）中的诱导表达和Ligand blot试验鉴定将pET-26b-mALP重组质粒转化至BL21（DE3）。挑取阳性克隆，经测序及双酶切验证后接种于含50 µg·mL-1卡那霉素的LB培养基中过夜培养，以pET-26b空载作为对照。次日按1%比例转接于含0.5%葡萄糖和50 µg·mL-1卡那霉素的LB培养基中，37℃，250 r/min振荡培养至OD600约为0.6—0.8，加入终浓度为0.5 mmol·L-1IPTG，25℃诱导表达16 h，过夜诱导表达产物与变性上样缓冲液混匀煮沸后经SDS-PAGE电泳分离后，用湿转法将蛋白转移至PVDF膜上，5% MPBS（含5%脱脂奶粉的PBS）室温封闭过夜，用PBST（含0.5%吐温-20的PBS）洗膜3次每次10 min，加入HRP标记的Anti His-Tag鼠单抗（1﹕3 000）作为二抗，室温孵育2 h后洗膜，用增强型HRP-DAB底物显色试剂盒显色，检测目的蛋白的表达。Ligand blot前期步骤均与Western blot一致，并增加pET26b空载诱导产物对照和mALP转膜后不与Cry毒素孵育对照，排除非特异性和假阳性。过夜封闭洗膜后，加入Cry1Ac毒素孵育2 h后洗膜，加入亲和纯化后的兔抗Cry1Ac多抗（笔者实验室自己制备）孵育2 h后洗膜，加入HRP-羊抗兔IgG（1﹕3 000）孵育2 h后洗膜，增强型HRP-DAB底物显色试剂盒显色，检测目的蛋白是否具有与Cry1Ac毒素结合能力。

1.2.4 小菜蛾mALP同源建模及模型评价 将本研究获得的mALP基因翻译后的氨基酸序列作为目标蛋白，用BLAST中的blastp suite对PDB蛋白结构数据库（http://blast.ncbi.nlm.nih.gov/Blast.cgi）分别进行多序列同源性搜索，取同源性最高的晶体结构（PDB ID：1K7H）作为同源建模模板，同步将蛋白序列提交至I-TASSER进行建模。Cry1Ac在PDB数据库中已有晶体结构（PDB ID：4ARY），可直接下载其pdb文件。利用Swiss-Model（http://swissmodel. expasy.org/）和I-TASSER（http://zhanglab.ccmb.med.umich.edu/I- TASSER/），分别预测mALP的三维结构。将模型能量最小化后，采用OpenMM Zephyr程序对其进行分子动力学模拟。最终将获得的模型结构采用生物信息学在线网站http://services.mbi.ucla.edu/SAVES/中的PROCHECK、ERRAT以及Verify3D模块进行评估，检测其合理性，比较分析后选取最佳mALP三维模型进行后续研究。

1.2.5 小菜蛾mALP与Cry1Ac分子对接 将Cry1Ac和mALP的pdb文件利用PatchDOCK对接程序分析（http://bioinfo3d.cs.tau.ac.il/PatchDock/），得到200个对接复合物，进一步利用其FireDock的refine模块，获得综合打分排序的复合物，同样进行能量最小化和分子动力学模拟验证。验证后的对接复合物进行ASA计算，分析毒素和受体参与结合的区域；并进一步通过Ligplot分析毒素和受体结合区域中参与互作的氨基酸位点。

1.2.6 关键氨基酸位点预测 利用DrugscorePPI在线服务（http://cpclab.uni-duesseldorf.de/dsppi/main.php），提交毒素-受体对接复合物pdb文件，基于丙氨酸扫描突变，计算出突变前后结合自由能的变化值，当结合自由能的变化值＞2 kJ·mol-1时被认为是关键氨基酸位点，通过计算分别预测出毒素和受体参与结合的关键氨基酸位点。

2 结果

2.1 小菜蛾mALP基因扩增

提取小菜蛾4龄幼虫中肠总RNA，利用反转录试剂盒成功制备cDNA。通过针对mALP全长设计的引物，PCR扩增获得了大小约为1 655 bp的条带，与预期大小一致（图1）。

M：DNA Marker；1：mALP 基因PCR扩增产物 PCR product of mALP gene

2.2 重组质粒pET-26b-mALP的构建及鉴定

将纯化后的PCR产物连接至pET-26b载体中，经PCR验证条带大小与预期大小一致（约1 700 bp）（图2-A），经Ⅰ和Ⅰ酶切得到mALP片段（图2-B），并对其序列进行测定。序列分析结果表明该小菜蛾mALP基因编码区序列为1 632 bp，编码567个氨基酸，理论预测分子量约为62 kD。经NCBI网站Blastp搜索比对后发现与已公布的小菜蛾mALP序列一致性为99%（仅有6个氨基酸存在差异）（图2-C）。说明pET-26b-mALP载体构建成功。

2.3 pET-26b-mALP在BL21中的诱导表达和Ligand blot鉴定

将鉴定无误的阳性克隆在25℃，转速220 r/min，终浓度为0.5 mmol·L-1IPTG条件下，诱导培养16 h进行表达。经SDS-PAGE电泳发现，全菌液在预期大小处出现明显的特异性蛋白条带（62 kD）并经Western blot证实了成功表达出小菜蛾mALP蛋白（图3-A、3-B），Ligand blot试验则进一步证实了通过原核表达的mALP具有与Cry1Ac体外结合的能力（图3-C）。

2.4 小菜蛾mALP同源建模及模型评价

通过SWISS-MODEL和I-TASSER，分别构建了mALP的三维结构模型、能量最小化和分子动力学模拟后，利用在线模型评价方法对其三维结构进行评价比较分析，结果发现I-TASSER方法构建的mALP三维结构最佳（图4-A）。mALP的三维结构以及与模板的RMSD值为0.24（图4-B）；经分子动力学模拟，模型构像趋于稳定后（图4-C）进行模型评价，结果表明构建的模型各项参数指标（表1）均显示模型是合理可信的，可用于后续分子对接试验。

M：DNA Marker。A：1：重组质粒PCR产物 PCR product amplified by pET-26b-mALP recombinant plasmid；B：1：pET-26b-mALP重组质粒 pET-26b-mALP recombinant plasmid；2：重组质粒双酶切产物 Recombinant plasmid digested production by NcoⅠand NotⅠ；C: 本研究mALP（Query）与公布的小菜蛾mALP（Sbjet）序列比对Amino acid sequence alignment between mALP in this study (Query) and mALP in GenBank (GenBank ID: KC841472.2) (Sbjet)

M：蛋白Marker Protein Marker。A、B： 1：pET-26b空载 Empty pET-26b vector；2：诱导表达mALP Expression of mALP after IPTG induction； C：1、2：pET-26b 空载和诱导表达mALP与Cry1Ac毒素和抗体孵育 Expression productions of empty pET-26b vector and mALP were incubated with Cry toxin and anti-Cry toxin antibody；3：诱导表达mALP与CBS缓冲液和抗体孵育 Expression of mALP were incubated with CBS buffer and anti-Cry toxin antibody

图4 mALP三维结构（A）、与模板（红色）结构比对（B）、分子动力学模拟-RMSD值（C）及拉氏图分析（D）

表1 模型评价

2.5 小菜蛾mALP与Cry1Ac分子对接

利用PatchDOCK对接程序分析，得到200个对接复合物，进一步利用其FireDock的refine模块，获得综合打分排序的复合物，同样经过能量最小化和分子动力学模拟，最终确定最佳的对接复合物（图5-A）。通过溶剂可及表面积变化值（solvent accessible surface areas，ASA）计算分析发现Cry1Ac毒素Domain II和Domain III的ASA值分别为823 Å2和622 Å2，对结合受体的贡献相对最重要，其中Loop环2的ASA值为586 Å2，受体中参与结合的氨基酸总的ASA变化值为1186 Å2（表2）。进一步利用Ligplot分析发现毒素和受体分别都有6个氨基酸残基参与互作，并且都有大量的氨基酸残基以疏水形式参与结合（图5-B）。

表2 溶剂可及表面积变化值计算

B：形成氢键氨基酸中红色氨基酸字母与数字代表Cry1Ac毒素，绿色氨基酸字母与数字代表mALP受体；形成疏水性氨基酸中蓝色氨基酸字母与数字代表Cry1Ac毒素，黑色氨基酸字母与数字代表mALP受体 In formation of the hydrogen bonds amino acids the red amino acid names and numbers represented as Cry1Ac toxin, the green amino acid names and numbers represented as mALP receptor, respectively. While in formation of the hydrophobic interaction the blue amino acid names and numbers represented as Cry1Ac toxin, the black amino acid names and numbers represented as mALP receptor

2.6 关键氨基酸位点预测

经过DrugscorePPI计算预测对接复合物结合区域中主要参与结合的各残基结合自由能变化值，结果如图6-A所示，毒素和受体结合位点中分别有3个氨基酸残基（376ASN、443SER和486SER）和4个氨基酸残基（452ARG、499THR、502TYR和513TYR）为参与互作的关键氨基酸位点（结合自由能的变化值＞2 kJ·mol-1）；这些关键氨基酸位点的分布情况如图6-B所示。

A：DrugscorePPI计算预测复合物中主要参与结合的各残基结合自由能变化DrugscorePPI calculated the change of binding free energy of each residue in the predicted complexes；B：粉红色骨架区Pink ribbon region：Cry1Ac毒素 Cry1Ac toxin；红色球形氨基酸Red spherical amino acids：Cry毒素关键氨基酸位点Key amino acids of Cry1Ac toxin；黄色骨架区Yellow ribbon region：mALP receptor；蓝色球形氨基酸Blue spherical amino acids：mALP受体关键氨基酸位点Key amino acids of mALP receptor

3 讨论

目前，利用真核表达体系（杆状病毒-昆虫细胞）已成功表达了大量具Bt Cry毒素结合活性的各类受体[26-28]，其中ALP主要依靠真核表达的糖基化位点并利用糖基锚定（GPI-anchored）位点导致细胞穿孔[29]。而原核表达（大肠杆菌）作为应用最为广泛和成功的系统，也被研究者们应用于Cry毒素受体的表达和功能活性鉴定[30-31]，Chen等利用pET系统成功表达并鉴定了烟芽夜蛾和棉铃虫的ALP具有和Cry毒素结合的能力[16,32]。本研究将小菜蛾mALP基因克隆至pET-26b载体进行原核表达和配体印迹试验，证实了经原核表达的mALP同样具有与Cry1Ac毒素体外结合的能力。在原核表达试验中，分别优化了IPTG浓度、诱导温度、转速和诱导时间，但在所有优化的条件下，发现mALP均会形成一定量的包涵体，且经尿素溶解蛋白变性后的mALP在配体印迹试验中同样具有Cry毒素结合活性（结果未显示），说明mALP参与毒素结合的关键氨基酸表位极有可能是线性表位而非空间构象，这与Fernandez等[33]通过分段表达埃及伊蚊（）ALP受体蛋白区域与Cry11Aa毒素在变形条件下进行配体印记试验发现ALP具有两个结合表位区域相一致。

通过比较Cry毒素和靶标昆虫潜在受体（敏感性和抗性）在结合模式下参与互作氨基酸的变化，对了解毒素-受体作用的分子机理和利用蛋白质工程方法提高Cry毒素杀虫效果具有重要意义[34-35]。但由于Cry毒素种类繁多以及其在不同昆虫中的潜在受体不断被发现且分子量普遍偏大，而目前在蛋白质数据库（protein data bank，PDB）中仅有十几种Cry毒素晶体结构和双翅目昆虫疟蚊的氨基肽酶受体APN晶体结构被解析[36-37]。所以利用同源建模和分子对接可为研究毒素-受体结合模式提供强有力的支撑[38-39]。虽然配体印迹试验是在受体蛋白变形条件下进行的，但结果证实了ALP结合Cry毒素是不依赖其空间构象的，在非变性条件下同样具备和Cry毒素结合的能力[29]，所以可以利用预测的小菜蛾ALP三维结构与Cry毒素进行分子模拟研究。利用Swiss-Model和I-TASSER分别进行了mALP的三维建模，经过能量最小化和分子动力学分析，以及模型评价分析，最终确定了I-TASSER构建的mALP三维模型最佳，可一步通过分子对接分析其结合表位特征。

通过PatchDOCK和FireDock程序以及分子动力学模拟方法，对Cry1Ac-mALP进行分子对接和对接复合物评价，获得参数最优的对接复合物，经溶剂可及表面积变化和Ligplot氢键及疏水性分析，发现Cry1Ac毒素的结构域II和结构域III均参与受体结合，结构域II中的Loop2和Loop3区域对受体结合的贡献相对更大；Loop3区域中有2个氨基酸残基与受体形成氢键，而Loop2区域中有多达7个氨基酸与受体可形成疏水结合作用，这对形成稳定和牢靠的复合物至关重要，且与目前报道的Cry毒素Loop2区域主要参与结合ALP的研究相一致[40-41]。mALP有6个氨基酸与毒素形成氢键结合，更多的氨基酸则是通过疏水作用与受体结合。经过计算机辅助丙氨酸突变扫描，预测了毒素和受体中参与互作的关键氨基酸残基-“热点”（Hot-Spots），它们对蛋白质间的结合自由能（∆GBind）有着显著的贡献，是毒素与受体相互作用主要依赖的少数几个氨基酸残基[42]。本研究中预测出的这些热点残基可以在接下来的试验中进一步通过生物分子互作试验（例如Biacore）验证，并且mALP的热点残基对于今后认知和定位受体的毒素结合区域以及在比较敏感型和抗性昆虫受体氨基酸变异上有很大帮助，同时Cry毒素中的热点残基也是今后利用现代分子生物学手段改造毒素用于提高毒素活力和改善抗性的首要选择。

4 结论

利用原核表达系统，成功制备了能够结合Cry1Ac毒素的小菜蛾mALP；预测出小菜蛾mALP的三维结构及与毒素的结合模式，并在此基础上通过计算机辅助丙氨酸突变扫描获得毒素-受体结合的关键氨基酸位点，为进一步阐明Cry毒素与mALP的互作机制及分子改造增强Cry毒素活性的研究打下基础。

References

[1] 喻子牛, 孙明, 刘子铎, 戴经元, 陈亚华, 喻凌, 罗曦霞. 苏云金芽孢杆菌的分类及生物活性蛋白基因. 中国生物防治, 1996, 12(2): 85-89.

Yu Z N, Sun M, Liu Z D, Dai J Y, Chen Y H, Yu L, Luo X X. The classification ofand their biological active protein genes., 1996, 12(2): 85-89. (in Chinese)

[2] Bravo A, Likitvivatanavong S, Gill S S, Soberón M.: a story of a successful bioinsecticide., 2011, 41(7): 423-431.

[3] Vachon V, Laprade R, Schwartz J L. Current models of the mode of action ofinsecticidal crystal proteins: a critical review., 2012, 111(1): 1-12.

[4] Fabrick J A, Tabashnik B E. Binding oftoxin Cry1Ac to multiple sites of cadherin in pink bollworm., 2007, 37(2): 97-106.

[5] 潘家荣, 乔艳红, 张维, 林敏, 张杰. Bt晶体蛋白CrylAc放射免疫检测技术研究. 核农学报, 2006, 20(6): 544-547.

Pan J R, Qiao Y H, Zhang W, Lin M, Zhang J. Study on radioimmunoassary of Bt Cry1Ac., 2006, 20(6): 544-547. (in Chinese)

[6] Lee D W, Choi J Y, Kim W T, Je Y H, Song J T, Chung B K, Boo K S, Koh Y H. Mutations of acetylcholinesterase1 contribute to prothiofos-resistance in(L.)., 2007, 353: 591-597.

[7] 李怡萍, 梁革梅, 仵均祥, 陈豪, 马康生, 吴孔明, 郭予元. 苏云金芽孢杆菌杀虫机理及害虫对其抗性机制的研究进展. 西北农林科技大学学报(自然科学版), 2010, 38(9): 118-128.

Li Y P, Liang G M, Wu J X, Chen H, Ma K S, Wu K M, Guo Y Y. Progress in insecticidal mechanism of Bt and resistance mechanism of pest insect to Bt., 2010, 38(9): 118-128. (in Chinese)

[8] 张丽阳, 刘承兰. 昆虫抗药性机制及抗性治理研究进展. 环境昆虫学报, 2016, 38(3): 640-647.

Zhang L Y, Liu C L. Research progress on mechanism of insect resistance to insecticides and its management., 2016, 38(3): 640-647. (in Chinese)

[9] Hecke D G, Gahan L J, Baxter S W, Zhao J Z, Shelton A M, Gould F, Tabashnik B E. The diversity of Bt resistance genes in species of Lepidoptera., 2007, 95: 192-197.

[10] 谭声江, 陈晓峰, 李典谟. 昆虫对Bt毒素的抗性机理研究进展. 昆虫知识, 2001, 38(1): 12-17.

Tan S J, Chen X F, Li D M. Research progresses on mechanism of insect resistance totoxin., 2001, 38(1): 12-17. (in Chinese)

[11] Pigott C R, Ellar D J. Role of receptors incrystal toxin activity., 2007, 71(2): 255-281.

[12] 张长胜, 来鲁华. 蛋白质相互作用预测、设计与调控. 物理化学学报, 2012, 28(10): 2363-2380.

Zhang C S, Lai L H. Protein-protein interaction: Prediction, design, and modulation.,2012, 28(10): 2363-2380. (in Chinese)

[13] Rosenfeld L, Heyne M, Shifman J M, Papo N. Protein engineering by combined computational andevolution approaches., 2016, 41(5): 421-433.

[14] Jurat-Fuentes J L, Adang M J. Characterization of a Cry1Ac- receptor alkaline phosphatase in susceptible and resistantlarvae., 2004, 271: 3127-3135.

[15] Jurat-Fuentes J L, Karumbaiah L, Jakka S R K, Ning C, Liu C, Wu K, Jackson J, Gould F, Blanco C, Portilla M, Perera O, Adang M. Reduced levels of membrane-bound alkaline phosphatase are common to Lepidopteran strains resistant to Cry toxins from, 2011, 6(3): e17606.

[16] Chen W, Liu C, Xiao Y, Zhang D, Zhang Y, Li X, Tabashnik B E, Wu K. A toxin-binding alkaline phosphatase fragment synergizes Bt toxin Cry1Ac against susceptible and resistant., 2015, 10(4): e0126288.

[17] 杨峰山, 吴青君, 徐宝云, 曹丽波, 朱国仁, 张友军. 小菜蛾对Bt毒素Cry1Ac和Bt制剂抗性的选育及其抗性种群的生物学适应性. 昆虫学报, 2006, 49(1): 64-69.

Yang F S, Wu Q J, Xu B Y, Cao L B, Zhu G R, Zhang Y J. Resistance selection ofby Cry1Ac toxin and Bt formulation and biological fitness of the resistant populations., 2006, 49(1): 64-69. (in Chinese)

[18] Guo Z J, Kang S, Zhu X, Wu Q J, Wang S L, Xie W, Zhang Y J. The midgut cadherin-like gene is not associated with resistance totoxin Cry1Ac in(L.)., 2015, 126: 21-30.

[19] 张涛, 张丽丽, 魏纪珍, 肖玉涛, 梁革梅, 周瑞阳. Cry1Ac抗、感棉铃虫碱性磷酸酯酶(ALP1)的表达量比较. 中国农业科学, 2013, 46(17): 3580-3586.

Zhang T, Zhang L L, Wei J Z, Xiao Y T, Liang G M, Zhou R Y. The expression level of alkaline phosphatase (ALP1) in Cry1Ac-resistant and susceptible cotton bollwormHübner., 2013, 46(17): 3580-3586. (in Chinese)

[20] Yang Z X, Wu Q J, Wang S L, Chang X, Wang J H, Guo Z J, Lei Y Y, Xu B Y, Zhang Y J. Expression of cadherin, aminopeptidase N and alkaline phosphatase genes in Cry1Ac- susceptible and Cry1Ac-resistant strains of(L.).,2012, 136: 539-548.

[21] 刘洁, 文礼章, 王少丽, 吴青君, 张友军. 小菜蛾碱性磷酸酶基因cDNA片段的克隆及其序列分析. 昆虫知识, 2010, 47(2): 270-274.

Liu J, Wen L Z, Wang S L, Wu Q J, Zhang Y J. Cloning and sequence analysis of alkaline phosphatase gene in the diamondback moth,., 2010, 47(2): 270-274. (in Chinese)

[22] Guo Z, Kang S, Chen D F, Wu Q J, Wang S L, Xie W, Zhu X, Baxter S W, Zhou X G, Juerat-Fuentes J L, Zhang Y J. MAPK signaling pathway alters expression of midgut ALP and ABCC genes and causes resistance toCry1Ac toxin in diamondback moth., 2015, 11(4): e1005124.

[23] 李红亮, 张林雅, 庄树林, 倪翠侠, 韩宝瑜, 商晗武. 中华蜜蜂普通气味结合蛋白ASP2的气味结合功能模式分析. 中国农业科学, 2013, 46(1): 154-161.

Li H L, Zhang L Y, Zhuang S L, Ni C X, Han B Y, Shang H W. Interpretation of odorant binding function and mode of general odorant binding protein ASP2 in Chinese honeybee ()., 2013, 46(1): 154-161. (in Chinese)

[24] 庄绪静, 尹姣, 李克斌, 曹雅忠. 华北大黑鳃金龟气味结合蛋白HoblOBP2的生物信息学分析. 植物保护, 2013, 39(1): 50-55.

Zhuang X J, Yin J, Li K B, Cao Y Z. Bioinformatics analysis of the odorant-bindig protein HoblOBP2 in olfactory sensilla of the scarab bettle.,2013, 39(1): 50-55. (in Chinese)

[25] Shan S P, Xia L Q, Ding X Z, Zhang Y M, Hu S B, Sun Y J, Yu, Z Q, Han L Z. Homology modeling of Cry1Ac toxin-binding alkaline phosphatase receptor fromand its functional interpretation., 2011, 29: 427432.

[26] Zhang X, Candas M, Griko N B, Rose-Young L, Jr Bulla L A. Cytotoxicity ofCry1Ab toxin depends on specific binding of the toxin to the cadherin receptor BT-R1 expressed in insect cells., 2005, 12: 1407-1416.

[27] 王俊华, 周小毛, 吴青君, 王少丽, 谢文, 陈得峰, 徐宝云, 张友军. 小菜蛾中肠氨肽酶N2在昆虫细胞中的表达. 农药学学报, 2012, 14(2): 151-157.

Wang J H, Zhou X M, Wu Q J, Wang S L, Xie W, Chen D F, Xu B Y, Zhang Y J. Expression of aminopeptidase N (APN2) frommidgut in Sf9 cells., 2012, 14(2): 151-157. (in Chinese)

[28] Zhang X B, Griko N B, Corona S K, Jr Bulla L A. Enhanced exocytosis of the receptor BT-R1 induced by the Cry1Ab toxin ofdirectly correlates to the execution of cell death., 2008, 149(4): 581-588.

[29] Ning C, Wua, K, Liu C, Gao Y, Jurat-Fuentes J L, Gao X. Characterization of a Cry1Ac toxin-binding alkaline phosphatase in the midgut from(Hübner) larvae., 2010, 56: 666-672.

[30] 李红亮, 聂文敏, 高其康, 程家安. 中华蜜蜂气味结合蛋白ASP2 cDNA的克隆及原核表达. 中国农业科学, 2008, 41(3): 933-938.

Li H L, Nie W M, Gao Q K, Cheng J A. Cloning of cDNA encoding odorant binding protein ASP2 in working bee’s antenna ofand its prokaryotic expression., 2008, 41(3): 933-938. (in Chinese)

[31] 林慧岩, 周子珊, 束长龙, 高继国, 张杰. APN1的原核表达及其与Cry1Ac蛋白体外结合. 植物保护, 2015, 41(4): 23-28.

Lin H Y, Zhou Z S, Shu C L, Gao J G, Zhang J. Prokaryotic expression of aminopeptidase N1 fromand binding analysis with Cry1Ac toxin., 2015, 41(4): 23-28. (in Chinese)

[32] Perera O P, Willis J D, Adang M J, Jurat-Fuentes J L. Cloning and characterization of the Cry1Ac-binding alkaline phosphatase (HvALP) from., 2009, 39: 294-302.

[33] Fernandez L E, Martinez-Anaya C, Lira E, Chen J, Evans A, Hernández-Martínez S, Lanz-Mendoza H, Bravo A, Gill S S, Soberón M. Cloning and epitope mapping of Cry11Aa-binding sites in the Cry11Aa-receptor alkaline phosphatase from., 2009, 48: 8899-8907.

[34] 张丽丽, 梁革梅, 曹广春, 高希武, 郭予元. 增强Bt Cry毒素杀虫作用的重要途径: 增效因子的利用及晶体蛋白的遗传改良. 环境昆虫学报, 2010, 32(3): 525-531.

Zhang L L, Liang G M, Cao G C, Gao X W, Guo Y Y. Approaches for enhancing the insecticidal activity ofCry toxins: application of synergistic factors and genetic improvement of crystal protein., 2010, 32(3): 525-531. (in Chinese)

[35] Zhang S P, Cheng H M, Gao Y L, Wang G R, Liang G M, Wu K M. Mutation of an aminopeptidase N gene is associated withresistance toCry1Ac toxin., 2009, 39: 421-429.

[36] 束长龙, 张风娇, 黄颖, 李艳秋, 张杰. Bt杀虫基因研究现状与趋势. 中国科学: 生命科学, 2016, 46(5): 548-555.

Shu C L, Zhang F J, Huang Y, Li Y Q, Zhang J. Current status and research trends of Bt insecticidal gene.,2016, 46(5): 548-555. (in Chinese)

[37] Atkinson S C, Armistead J S, Mathias D K, Sandeu M M, Tao D Y, Borhani-Dizaji N, Tarimo B B, Morlais I, Dinglasan R R, Borg N A. Structural analysis ofmidgut aminopeptidase N reveals a novel malaria transmission- blocking vaccine B-cell epitope., 2015, 22(7): 532-539.

[38] Tajnea S, Sanam R, Gundla R, Gandhi N S, Mancera R L, Boddupally D, Vudem D R, Khareedu V R. Molecular modeling ofCry1Ac (DI-DII)-ASAL (lectin)-fusion protein and its interaction with aminopeptidase N (APN) receptor of., 2012, 33: 61-76.

[39] Ahamd A, Javed M R, Rao A Q, Khan M A, Ahad A, Din S U, Shahid A A, Husnain T.determination of insecticidal potential offusion protein against Lepidopterantargets using molecular docking., 2015, 6: Article 1081.

[40] Shao E, Lin L, Chen C, Chen H, Zhuang H H, Wu S Q, Shao L, Guan X, Huang Z P. Loop replacements with gut-binding peptides in Cry1Ab domain II enhanced toxicity against the brown planthopper,(Stål)., 2016, 6(12): e84022.

[41] Bravo A, Gómeza I, Porta H, García-Gómez B I, Rodriguez-Almazan C, Pardo L, Soberón M. Evolution ofCry toxins insecticidal activity., 2012, 6(1): 17-26.

[42] Cukuroglu E, Engin H B, Gursoy A, Keskin O. Hot spots in protein-protein interfaces: Towards drug discovery., 2014, 116(2/3): 165-173.

（责任编辑 岳梅）

Expression and molecular simulation of alkaline phosphatase receptor of

ZHANG Xiao, HU Xiao-dan, ZHONG Jian-feng, WU Ai-hua, XU Chong-xin, LIU Yuan, ZHANG Cun-zheng, XIE Ya-jing, LIU Xian-jin

(Institute of Food Quality Safety and Detection Research, Jiangsu Academy of Agricultural Sciences/Key Laboratory of Food Quality and Jiangsu Province-State Key Laboratory Breeding Base/Key Laboratory of Control Technology and Standard for Agro-product Safety and Quality, Ministry of Agriculture, Nanjing 210014)

【Objective】The objective of this study is to confirm the binding ability of membrane-bound alkaline phosphatase (mALP) ofwith Cry1Ac toxin using prokaryotic expression and Ligand blot, and to predict toxin-receptor binding region and key amino acid binding sites (hot-spots) employed by homology modeling and docking study of Cry1Ac-mALP binding mode. It will provide a basis for the study of toxin-receptor interaction mechanism and molecular modification to enhance the activity of Cry toxin. 【Method】The mALP offull-length primers were designed and amplified by PCR. The restricted products of mALP and pET-26b (+) were ligated by T4 DNA ligase after the dual-enzyme digestion procedures. The recombinant pET-26b-mALP vectors were transferred into the Trans1-T1 phage resistant chemically competent cells, then picked clones were analyzed by PCR amplification, dual-enzyme digestion and sanger sequencing. The positive recombinant vectors (anchoring the corrected mALP gene) were transferred into theBL21 (DE3) competent cells for prokaryotic expression. The inducible expression products of mALP were transformed onto PVDF membrane. Preparation of mALP and binding activity of Cry1Ac with malp were verified through Western blot and Ligand blot, respectively. Three-dimensional structure of mALP was predicted by homology modeling, molecular dynamics simulation and model evaluation. The toxin-receptor docking complexes were generated by using the PatchDocK and FireDock web-servers with molecular dynamics simulations. The toxin-receptor complex was analyzed to determine the interaction region and the amino acid binding sites, key amino acid residues involved in Cry toxin and ALP receptor by computer-aided alanine mutation scanning tests. 【Result】mALP gene was successfully amplified, followed with the prokaryotic expression of mALP receptor protein. Binding of Cry1Ac toxin with prepared mALP protein was verified. The three-dimensional structure of mALP was successfully obtained by homology modeling, then the Cry toxin-ALP complex was determined. By the changed solvent accessible surface areas calculation and Ligplot analysis, the results showed that the domain II and domain III of Cry toxin were involved in binding to receptor, and Cry toxin and ALP were interacted mainly depending on hydrophobic and hydrogen bonding patterns. Finally, through the computer-aided alanine mutation scanning hot residues analysis, there were three key amino acid residues (376Asn, 443Ser and 486Ser) from Cry toxin and four key amino acid residues (452Arg, 499Thr, 502Tyr and 513Tyr) from ALP were participated in the interaction of toxin-receptor complex, respectively. 【Conclusion】It can be determined that the mALP receptor also has the ability to bind Cry1Ac toxin by prokaryotic expression, the three-dimensional structure of mALP was predicted and the toxin-receptor binding model was studied using molecular simulation.

alkaline phosphatase (ALP); prokaryotic expression; Ligand blot; homology modeling; molecular docking; hot-spots prediction

2016-09-05；接受日期：2016-10-08

国家自然科学基金（31401813、31630061）、江苏省自然基金-青年基金（BK20140744）、江苏省食品质量安全重点实验室-省部共建国家重点实验室培育基地自主研究课题（3201613）

张霄，E-mail：zxwin2008@126.com。通信作者刘贤金，E-mail：jaasliu@jaas.ac.cn