激素受体AmEcR的化合物模拟筛选及应用

埃及基础与应用科学杂志4（2017）288完整文章Ariadne merioneecdysone receptor（AmEcR）蛋白：一种比较激动剂和拮抗剂化合物的计算机模拟Al-Sundaravadivelana，Easwaran Murugeshb，Mathew Preethyc，Prabu Sivaprasathc，a印度泰米尔纳德邦塞勒姆636 011，帕德玛瓦尼女子艺术和科学学院动物学系b印度泰米尔纳德邦哥印拜陀641 046 Bharathiar大学生物信息学系c印度泰米尔纳德邦哥印拜陀641 021 Eachanari Post Karpagam大学生命科学学院生物技术系阿提奇莱因福奥文章历史记录：2017年7月2日收到2017年10月8日收到修订版，2017年2017年10月18日在线提供保留字：阿里阿德涅·梅里奥内20羟基蜕皮激素（20E）薛定谔A B S T R A C T蜕皮激素信号转导在昆虫变态过程中起着重要作用，20-羟基蜕皮激素（20-hydroxyecdysone，20 E）与蜕 皮激 素受 体 配体 结 合域 （ EcR-LBD） 组 成的 核受 体 结合 ， 触发 昆 虫的 发 育转 变。 对 阿里 阿德 涅（ Ariadne merione ） 蜕 皮 激 素 受 体 （ AmEcR ） cDNA 进 行 了 扩 增 和 部 分 序 列 测 定 。 利 用 BIOEDIT（v7.2.5）软件预测AmEcR蛋白的理论分子量（MW）为21.192 kDa，等电点（pI）为9.31，脂肪指数（PI）为101.739。鉴定了A. merione与Precis coenia基因的相似性最大。本研究利用配体-受体工程技术，筛选具有拮抗作用的化合物，辅助配制昆虫特异性杀虫剂。 使用来自Zinc数据库的5554个分子进行了使用Schrödinger maestro和虚拟筛选建模的AmEcR的EcR蛋白质3D结构，其中ZINC 20031812显示了最高的滑移评分-6.257，并且根据文献选择了乙恶唑并显示了最佳滑移评分-6.671。用分子动力学模拟方法对拮抗剂和激动剂（20E）进行了激动剂和拮抗剂的均方根偏差（RMSD）值表明结合在水中是稳定的，其距离范围为2.3至2.6 μ m，1.8至1.5 μ m。2.3 1.9 ~ 2.3 μ m，时间尺度变化为1 ps。由于Etoxazole和ZINC 20031812是拮抗剂，因此计算它们比20E更稳定。©2017曼苏拉大学。由爱思唯尔公司制作和主持这是一篇开放获取的文章，CC BY-NC-ND许可证（http://creativecommons.org/licenses/by-nc-nd/4.0/）。1. 介绍半变态和全变态昆虫发育的一个共同特征是类固醇激素20羟基蜕皮激素（20 E）的周期性脉冲决定了每个发育过渡。蜕皮激素受体复合物是昆虫蜕皮过程中蜕皮激素诱导的生理和形态变化的关键元件，蜕皮激素如20-羟基蜕皮激素（20 E）及其类似物与蜕皮激素受体的配体结合结构域结合[1]。虽然20 E在全变态昆虫中的分子作用已被广泛研究，但在半变态昆虫中的数据很少[2]。激素20 E是亲脂性的，其合成并释放到血淋巴中，并进入细胞负责变态。蜕皮激素受体基因是一种核受体，具有两个主要结构域，称为配体结合结构域（LBD）和DNA结合结构域（DBD），*通讯作者。电子邮件地址：sivaprasathibt@gmail.com（P. Sivaprasath）。在细胞质中处于非活性状态。20 E与蜕皮激素受体基因的LBD结合，并激活DBD。蜕皮激素受体蛋白（ecdysone receptor protein，EcR）的基本结构由5个结构域组成，A/B转录激活结构域、C-DNA结合结构域、D铰链区和E配体结合结构域[3]。蜕皮过程由核受体超家族中的许多转录因子启动这一结果表明，激素途径中的几个晚期基因上调，有助于介导蜕皮过程[4]。这些核受体由于其异二聚体性质而具有穿过核孔的能力，并结合到DNA的特定位点，其经历主动复制、染色体重塑并最终导致翻译。蜕皮期间角质层的降解和新角质层的形成是蜕皮激素的主要信号功能之一，而N-乙酰基-β-D-葡萄糖胺（几丁质）是昆虫角质层的主要成分[5]。鳞翅目昆虫中的一些是严重的农业害虫，通过落叶破坏作物本研究是基于对特定发育基因及其编码蛋白的生物信息学分析而提出的假设。https://doi.org/10.1016/j.ejbas.2017.10.0022314- 808 X/©2017曼苏拉大学。由爱思唯尔公司制作和主持这是一篇基于CC BY-NC-ND许可证的开放获取文章（http://creativecommons.org/licenses/by-nc-nd/4.0/）。可在ScienceDirect上获得目录列表埃及基础与应用科学杂志杂志主页：www.elsevier.com/locate/ejbasC. Sundaravadivelan等人/埃及基础与应用科学杂志4（2017）288289××特异性发育蛋白（昆虫特异性）的靶向可能有助于控制目标害虫。筛选出的化合物模拟20羟基蜕皮激素，精确干扰受体功能，有助于控制农业害虫的大量数量这可以通过使用有效的生物信息学工具来设计与蜕皮激素受体蛋白对接的完美分子来实现2. 材料和方法2.1. RNA分离从四龄幼虫中提取总RNA。使用TRIzol试剂（Ambion ®，LifeTechnology，USA），将所述细胞培养在merione（在蜕皮间阶段）中。一个个体（=100 mg）A. 使用无菌特氟隆手动均化器，用1ml TRIzol试剂充分研磨四龄的merione。将含有匀浆的试管在室温下孵育5 min，以使核完全解离为了制备α蛋白复合物，加入200 μl预冷氯仿，剧烈混合15 s，并将悬浮液在室温下孵育2-3 min，然后离心（7500× g，4 ℃，15 min）。将上清液（~ 500 μ l）转移至无菌管中，之后加入500 μ l 75%异丙醇并在室温下孵育10 min，然后离心（7500 g，4 ℃，15 min）。 弃去上清液，然后将Iml的75%乙醇加入到沉淀中，涡旋并离心（6800g，4 °C，5min）。将上清液完全弃去，最后将沉淀风干，再悬浮于20 μl DEPC水中，并储存在-80 °C以供进一步使用。2.2. 聚合酶链反应从A.根据制造商的方案，使用第一链cDNA合成试剂盒（Roche，Germany）逆转录merione幼虫。两个基因特异性引物，正向引物（50-AGATGACCATCCTCACCGTG-3 0 ） 和 反 向 引 物 （ 5 0 -AGATGACCATCCTCACCGTG -30）Fig. 1. 编码阿里阿德涅梅里翁蜕皮激素受体蛋白的部分核苷酸序列（1290C. Sundaravadivelan等人/埃及基础与应用科学杂志4（2017）288≤≤≤≤≤2.3. DNA测序和序列分析在Eurofins Scientific Ltd.（Bangalore，India）使用Sanger双脱氧 技 术 对 纯 化 和 提 取 的 样 品 进 行 测 序 。 采 用 BIOEDIT 软 件（ver.7.2.5）对获得的序列进行分析，并对AmEcR核苷酸及其编码蛋白进行比对。2.4. 同源建模和结构验证表1图二. AmEcR蛋白的3D建模结构。使用Schrödinger maestro（版本9.3）（Schrödinger Inc.）预测AmEcR翻译蛋白质序列的从头算建模使用Prime-SP（版本3.1）（标准精密度），它有助于比较建模，包括比对、构建结构、折叠识别和分子力学-广义Born模型，并增加了疏水溶剂可及表面积计算。使用ProCheck验证模型结构的范围和结构相似性针对AmEcR建模结构中存在的氨基酸总数分析甘氨酸和脯氨酸的残基使用Ramachan-dran图对该AmEcR模型化结构进行细化和验证Ramachandran图中AmEcR的残留百分比。立体化学稳定性EcR受体残基相似性（%）最受欢迎区域15894其他允许区域106在慷慨允许区域的00不允许区域00非甘氨酸和非脯氨酸残基168末端残基数（不包括Gly和Pro）2甘氨酸残基的数量（显示为三角形）7脯氨酸残基7残基总数和相似184100.00图三. 模拟AmEcR蛋白的Ramachandran图显示在比对程序中与模板结构的根据其它鳞翅目昆虫EcRs基因序列的保守区，设计了5 0 -ACGTCCAGATCTCCTCGAG - 3 0聚体。将退火温度固定在57 °C持续35个循环。通过QIAquick凝胶提取试剂盒（Qiagen）纯化PCR产物。2.5. 网格生成使用Glide程序通过点击受体网格生成面板设置受体网格生成。在网格生成过程中，使用了默认参数，没有引入特殊约束。采用默认网格大小，所有活跃的网站。2.6. 配体回收根据PubChem数据库检索的文献调查[6，7]，选择啶虫脒、色虫酰肼、二苯甲酰肼、乙螨唑、唑螨酯、甲氧虫酰肼、吡丙醚和虫酰肼进行研究。另一方面，通过设置从ZINC数据库下载结构数据文件（SDF基于Tice规则的化合物性质。Tice规则指出，潜在的杀虫化合物应具有：（1）分子量小于500 g/mol，（2）氢供体的数目为3，（3）氢供体的数目为1，（4）氢供体的数目为1，（5）氢供体的数目为1，（6）氢供体的数目为的氢键受体12，（4）log P分配系数（亲油性）5，和（5）可旋转键数12[7，8]。下载了一个29兆字节的SDF文件，其中包含5554个化合物。 从PubChem数据库中检索20羟基蜕皮激素（20E），并作为对照，比较研究2.7. 虚拟筛选与对接在薛定谔模块中使用来自Zinc数据库的5554种化合物的SDF文件进行虚拟筛选。 此SDF文件已导入系统。开始虚拟筛选，然后进行500次迭代以分析配体结构和与受体结合的能力。基于虚拟筛选，得分最高的单一配体用于对接研究。虚拟筛选工作流面板设置LigPrep（2.6版）、Qik-Prop（3.5版）和Glide（5.8版）配体对接的输入文件，并按顺序将其提交给选定的主机。滑翔结果检查的重点是视觉而不是数字评估。第一组练习使用项目表来显示标准精密（SP）滑翔对接工作的结果，检查单个配体姿势及其与输入受体结构的接触。第二组练习使用Glideexpress precision（XP）可视化工具C. Sundaravadivelan等人/埃及基础与应用科学杂志4（2017）288291表2基于Tice规则选择激动剂和拮抗剂分子进行对接研究配体pHxlogP拓扑极表面积（K2）氢键供体氢键受体净电荷分子量（g/mol）可旋转键20E（对照）7.01.3660.2370480.6425乙螨6.23.5653.3050359.405ZINC200318127.02.9756150393.5125见图4。所选激动剂和拮抗剂的2D结构，A）20-羟基蜕皮激素，B）乙恶唑和C）ZINC 20031812。表3AmEcR模拟蛋白质和配体相互作用谱。配体名称交互部分粘合长度（mm）Glide评分氨基酸位置氨基酸部分配体部分20EALA 54THR 2TYR 64HHHOoo2.022.372.08-8.248乙螨ZINC20031812160岁ARG 39PHE 53TYR 64THR 2160岁TYR 64OHCHHOHHOHOOHO1.81二点三三，二点一三2.701.961.751.971.74-6.671-6.257面板显示Glide XP评分函数中有助于配体结合的术语信息2.8. 蛋白质-配体复合物使用Desmond（3.1版）（D.E.Shaw Research，New York），其具有内置的优化的势配体模拟（OPLS-2005），其为计算生物分子结构之间的结合能提供了强有力的框架。本研究采用分子动力学方法对AmEcR在1000 ps平衡后，为100个样品生成1000 ps的轨迹。均方根偏差（RMSD）用于测量两个结构中相同类型原子之间的标量距离。在目前的研究中，我们使用C（alpha）原子之间的RMSD来测量蛋白质同源物之间的拟合。在这个计算中，我们使用重原子的RMSD来比较时间上的结构和原始结构（在时间= 0ps时）之间的空间偏差。通常情况下，重原子的RMSD在分子模拟时间的纳秒内不应改变超过3微秒。3. 结果和讨论3.1. DNA测序和序列分析在随后的基因清洗和PCR产物的测序中，使用BIOEDIT（版本7.2.5）的核酸含量显示，胞 嘧 啶 占 32.91% ， 其 次 是 鸟 嘌 呤 、 腺 嘌 呤 和 胸 腺 嘧 啶 ， 分 别 占26.58%、21.70%和18.81%。根据这些核苷酸数目，EcR的部分编码序列区（CDS）的长度为553 bp，并将AmEcR核苷酸与其编码蛋白进行了比 对 （ 图 1 ） 。 AmEcR 基 因 序 列 在 NCBI 数 据 库 中 提 交 （ 登 录 号KJ652504）。3.2. 同源建模和结构验证对AmEcR序列进行与具有每个原子构型的结构构象的蛋白质数据库的同源性建模。Schrödinger（ver.9.3）模块选择烟芽夜蛾（PDB登录代码该AmEcR与模板蛋白的同源性为89.31%。使用ProCheck验证AmEcR的建模结构。理想情况下，该结构显示约86%的残基在其核心区域中，这与比对程序中的模板结构的残基相似性约94%。核心区域中残留物的百分比是立体化学质量的更好指南（表1;图2）。 3）。3.3. 筛选化合物已知的20E配体参与变态过程中与受体结合。随后，筛选出5554种化合物，其中两种化合物，即Etoxazole和ZINC 20031812显示出对接相容性，因此用于进一步的研究（表2，图4）。这些分子292C. Sundaravadivelan等人/埃及基础与应用科学杂志4（2017）288--指定生物相关的质子化状态，并注释了诸如分子量、计算的Log P和可 旋 转 键 数 的 性 质 。 该 数 据 库 可 以 免 费 下 载（ http://zinc.docking.org ） ， 有 几 种 常 见 的 文 件 格 式 ， 包 括SMILES，mol2，3D SDF和DOCK flexibase格式[9]。虚拟筛选有助于鉴定与已知EcR配体相似的新型非甾体配体[10]。3.4. 对接和蛋白质-配体相互作用谱20 E与受体模型的对接显示，配体分子可以以与其他类固醇激素-受体复合物类似的方式与受体相互作用[11]。在本研究中，计算机对接显示蛋白质中的结合结构域和特定口袋中配体的相互作用如表3所示。 第64位的酪氨酸是与所有配体共同结合的氨基酸，这表明这些拮抗剂结合在与真正的20E结合的相同口袋中。 在分析了几种晶体结构后，发现DNA结合结构域（DBD）和配体结合结构域（LBD）在昆虫EcR中高度保守[12]。一种方法是通过使用的配体结合结构域的蜕皮激素受体（EcR）的蛾Heliotis virescens以及虚拟分子库的类似物的已知的二酰基肼（DAH）型蜕皮激素激动剂的晶体结构。通过将DAH与EcR的结合口袋对接，然后进行对接构象的CoMFA（比较分子场分析）和CoMSIA（比较分子相似性指数分析），可以绘制迄今为止未探索的受体腔区域[13]。3.5. 模拟AmEcR与20EAmEcR与20 E的分子对接结果表明，配体在5个位点上通过配体间的相 互 作 用 与 受 体 ALA54 （ H ） 、 THR2 （ H ） 、 TYR64 （ H ） 、ARG160（O）、ARG39（H）发生C（碳）、H（氢）和O（氧）形成碳、氢和氧键，键距为2.02、2.37、2.08、1.81、2.33和2.13埃（图5A）。从Glide评分8.248可以注意到AmEcR与20E的更高相互作用（表3）。氢键是配体与蜕皮激素受体结合的重要物理化学性质之一，在模拟的配体-受体复合物中手动计数配体分子与受体之间可能的氢键在变态过程中，LBD发挥主要作用，充当受体二聚化、配体识别和辅因子相互作用[14]。蜕皮激素受体和配体结合模型试图解释20 E和二苯甲酰肼如何与配体结合口袋相互作用，同源性模型复合物提供了新的见解，可用于合理设计新的环境安全杀虫剂[113.6. 模拟AmEcR与乙恶唑AmEcR与乙恶唑的相互作用表明，该配体在AmEcR的3个位点上相互作用，具有残留原子类型PHE 53（C）、TYR 64（H）和THR 2（H），结合距离为2.70、1.96和1.75 nm（图5B）。Glide评分6.671表明AmEcR与乙恶唑的相互作用最高（表3）。LBD含有配体结合口袋（LBP），其结合ecdys-图五、AmEcR与测试配体的相互作用曲线，A）AmEcR与20羟基蜕皮激素（20 E），B）AmEcR与乙恶唑和C）AmEcR与ZINC 20031812。第64位的酪氨酸是与所有配体共同结合的氨基酸，这表明Etoxazole和ZINC 20031812（拮抗剂）结合在20 E结合的相同口袋中C. Sundaravadivelan等人/埃及基础与应用科学杂志4（2017）288293-图第六章AmEcR-20 E复合物的分子动力学，A）RMSD图B）势能图和C）AmEcR和激动剂20 E的1- AmEcR和20 E复合物的研究表明，AmEcR在2.3和2.6 μ m的距离变化范围内是稳定的。该配合物的势能偏差为-26100 ~-26500E。类固醇以及某些非甾体类EcR激动剂，如基于DAH的杀虫剂[16]。3.7. 模拟AmEcR与ZINC 20031812将ZINC 20031812与AmEcR模型对接显示，配体在具有残余原子类型THR 160（O）、THR 64（H）的AmEcR的2个位点处相互作用，并且其结合距离为1.97和1.74 π（图5C）。Glide评分为6.257，表明AmEcR与ZINC 20031812存在潜在相互作用（表3）。3.8. AmEcR激动剂和拮抗剂复合物的分子模拟在分子模拟中，对应于原子的粒子的位置和速度根据经典物理定律演化[17]。本文通过计算RMSD值，给出了AmEcR-20 E、AmEcR-乙恶唑和AmEcR-ZINC 20031812复合物的叠加构象。RMSD测量指定帧的选定原子与参考帧相比的位移的平均变化，参考帧用于轨迹中的所有帧。在蛋白质RMSD中，蛋白质框架最初将在参考框架骨架上对齐，其中RMSD基于原子选择来计算这在整个模拟过程中提供了对蛋白质结构构象的清晰了解进一步的RMSD分析还将说明平衡模拟（如果有的话），其在模拟结束时围绕热平均结构的波动。在模拟过程中，对于小的球状蛋白质，从1到3 μ m的顺序变化变化远大于上述指定值表明模拟过程中蛋白质的构象变化较大。应注意，如果模拟收敛，则RMSD值稳定在固定值附近。在模拟结束时蛋白质的RMSD的增加或减少模式指示系统未平衡。配体RMSD表示配体相对于蛋白质及其结合口袋的稳定性将蛋白质-配体复合物的RMSD在参照物的蛋白质骨架上对齐，随后测量配体重原子的RMSD这表明，如果该值显著大于蛋白质的RMSD，则可能配体已扩散远离其初始结合位点[18]。3.9. AmEcR-20 E复合物根据AmEcR与20 E复合物的重叠构象计算出RMSD，表明AmEcR在2.3 ~ 2.6 μ m294C. Sundaravadivelan等人/埃及基础与应用科学杂志4（2017）288----图第七章AmEcR-乙恶唑复合物的分子模拟，A）RMSD图B）势能图和C）AmEcR和拮抗剂乙恶唑的1-AmEcR与乙恶唑复合物的重叠构象显示在1.8和2.3 nm之间的水溶液中结构稳定，势能偏差范围为-26000 E至-26400 E，较低的势能证实了结合稳定性。（图6A）。对AmEcR和20 E复合物的分析表明，其势能为模拟的，并且对于观察到的复合物结构具有低势能，范围为26100 E至26500 E（图6B）。势能的降低表明分子间键合的复杂性增加。从1到100的轨迹被叠加并检查复杂结构的运动。结构排列相对于水环境表现出更多的复杂性，因为环境是疏水性的，并且配体在与AmEcR的内部键合中具有更大的影响（图6C）。蛋白质-配体结合自由能的准确预测是计算药物设计的主要目标。使用液体模拟优化势（OPLS-2005）力场，并获得与实验溶剂化自由能的高度相关性和大多数官能团的低平均无符号误差[19]。3.10. AmEcR-拮抗剂复合物AmEcR和乙氧唑复合物的叠加构象结构表明，AmEcR在1.8和2.3μ m之间的距离在水中是稳定的（图1）。 7A）和ZINC 20031812显示从1.9到2.3 μ m的变化（图7A）。 8 A）。Etoxazole的作用模式揭示了几丁质生物合成抑制剂活性。在秋粘虫草地贪夜蛾中观察到蜕皮缺陷[20]。模拟的AmEcR和两种拮抗剂复合物具有范围从26000E至26400E的相同的电位偏差，较低的势能证实了结合稳定性（图7 B，图7B）。 8 B）。从1到100的轨迹被叠加并检查复杂结构的运动。结构对齐-由于水环境是疏水的，配体在与AmEcR的内部键合中具有更大的作用，因此相对于水环境，AmEcR表现出更复杂的性质（图7 C，图7C）。 8 C）。测试了乙螨唑对蔬菜主要害虫甜菜夜蛾（Spodopteraexigua）、小菜蛾（Plutella xylostella）、豆蚜（Aphis craccivora）和胭脂红蜘蛛（Tetranychus cinnabari nus）的几丁质合成抑制剂活性，并记录了其有效LC50[21]。比较激动剂和拮抗剂的RMSD值，表明它们在水中稳定，2.3 到2.6 μ m，1.8到2.3 μ m和1.9到2.3 μ m，在时间尺度上的变化为1 ps。由于Etoxazole和ZINC 20031812是拮抗剂，计算它们比20E更稳定在目前的工作中，我们提出了一个目标选择性潜在的化合物，可用于控制特定的昆虫害虫，而不会对环境中的其他非靶标生物造成在体外或体内研究之前，计算机模拟分析能够提供有关筛选、可能的作用模式和选择控制剂的更多信息。一种新的毛虫控制剂，命名为tebufenozide，其对非靶标生物的危害非常小，并且对生态系统安全它可用于广泛的农业害虫，并有效地取代对环境有毒的广谱杀虫剂[22]。与甾体和非甾体配体的复合物中的EcR结构的比较揭示了根本不同且仅部分重叠的配体结合口袋，其不能通过分子建模和对接研究预测[23]。非甾体化合物如RH 5849（1，2-二苯甲酰基-1-叔丁基肼）在烟草角虫Manduca sexta的所有幼虫发育阶段引起过早开始蜕皮。RH 5849的活性是其30至>670倍C. Sundaravadivelan等人/埃及基础与应用科学杂志4（2017）288295见图8。AmEcR-ZINC20031812复合物的分子模拟AmEcR和ZINC 20031812的叠加结构-26000E至-26400E。就像20-羟基蜕皮激素蜕皮激素一样[24]。通过将20E和二苯甲酰肼对接至蜕皮激素受体，天然和合成分子的可能的新型叠加使得能够设计环境友好的杀昆虫剂[15]。过去对节肢动物的控制主要依赖于化学杀虫剂，这些化学杀虫剂污染环境，并且对非目标物种也有影响。化学杀虫剂不仅消耗我们食物的营养价值，而且还会积累。研究不断发现食物中的农药残留，食用它会导致无数的疾病。世界卫生组织（世卫组织）估计，每年有300万起农药中毒病例，多达22万人死亡，主要发生在发展中国家，儿童更容易受到化学农药的伤害。现在，有更多的数据库提供了数千种植物基化合物。这些生物活性化合物可能是有效的控制严重的农业害虫，而且它们是生态友好的。这些生物活性化合物可能具有变构调节特性，它可以靶向蛋白质的其他潜在结合位点，从而有效地改变蛋白质的构象结构，抑制天然20-羟基蜕皮激素在其专用口袋中的结合。柠檬苦素类化合物是属于四降三萜类化合物组的化合物，通常，属于该组的化合物对昆虫表现出广泛的生物活性，如杀虫、拒食和生长调节活性[25]。 在Sf 9（昆虫）细胞系上分析了印楝内酯（一种柠檬酸类化合物）的细胞毒性作用，而9.8μ M的浓度对Sf细胞系作用迅速，并诱导了印楝内酯的细胞毒性作用。细胞膜的破坏[26]。此外，从伞形花（十字花科）种子中分离的葫芦素B和D已被证明在防止黑腹果蝇BII永久细胞系中20 E诱导的形态学变化方面具有拮抗活性[27]。许多鳞翅目昆虫的幼虫是农业上的主要害虫。主要害虫科为夜蛾科、螟蛾科和卷蛾科。阿里阿德涅merione属于鳞翅目和家庭若虫。虽然它不被认为是一种严重的害虫，但计算机模拟研究有助于保持关于目标蛋白质、化合物的数据，并有助于设计昆虫特异性生物农药。4. 结论为了从生防角度了解拮抗剂和拮抗剂分子的作用机制，对AmEcR由于使用生物信息学工具设计的受体特异性拮抗剂，生物防治方法将更加精确，并且这个概念将是昆虫特异性的，它有助于保护有益昆虫和环境。确认作者感谢管理层，卡尔帕加姆大学，高等教育卡尔帕加姆学院，哥印拜陀提供本研究期间所需的所有设施296C. Sundaravadivelan等人/埃及基础与应用科学杂志4（2017）288引用[1] Jayachandran B，Hussain M，Asgari S. miR-14对棉铃虫蜕皮激素受体的调控及其与杆状病毒感染的潜在联系JInvertebr Pathol 2013;114：151-7.[2] Mané-Padrós D，Borràs-Castells F，Belles X，Martín D.核受体HR4在蜕皮激素触发的半代谢昆虫德国小蠊（Blattella germanica）发育的基因级联反应中起重要作用。Mol CellEndocrinol 2012;348：322-30.[3] Thummel CS.从胚胎发生到变态：果蝇核受体超家族成员的调节和功能。Cell1995;83：871-7.[4] 郑万华，杨东东，王继新，宋圣强，李佳，赵晓芳。Hsc 70与超气门结合导致棉铃虫20-羟基蜕皮激素应答基因上调分子细胞内分泌学2010;315：282-91。[5] Merzendorfer H，Zimoch L.昆虫的几丁质代谢：几丁质酶和几丁质酶的结构、功能和调控。第206章. 4393-12[6] Smagghe G，Decombel L，Tirry L.四种蜕皮激素激动剂在多重抗性棉叶虫体内吸收、氧化及结合毒性的意义。昆虫生物化学生理学2001;46：127-9.[7] TiceCM. 选择合适的化合物进行筛选：Lipinski的药物5规则适用于农用化学品吗？害虫管理科学2001;57：3-16.[8] 吴伟杰，王晓刚，王晓刚.非印楝素类柠檬苦素与棉铃虫蜕皮激素受体的电子对接研究Med Chem Res2015;24：2621-31.[9] Irwin JJ，Shoichet BK. ZINC- 用 于 虚 拟 筛 选 的 市 售 化 合 物 的 免 费 数 据 库 。 JChem Inf Model 2005;45：177-82.[10] 张文忠，张文忠，张文忠.用蛋白质模拟系统评价蜕皮甾体在配体-受体相互作用中的氢键。Bioorg Med Chem2011;17：5868-73.[11] [10]张文，张文. 蜕皮激素激动剂与受体的结合模式：比较建模和对接研究。J MolModel 2003;9：58-65.[12] Iwema T，Chaumot A，Studer RA，Rechavi MR，Billas IML，Moras D，et al.昆虫中USP和蜕皮激素受体EcR之间异二聚化界面的结构和进化创新。Mol BiolEvol2009;26：753-68.[13] Bordas B，Belai I，Lopata A，Szanto Z.使用3D分子场解释评分函数。昆虫蜕皮激素受体二酰肼结合口袋的定位。J Chem Inf Model 2007; 47：176 - 85.[14] Nakagawa Y ， Henrich VC. 节 肢 动 物 核 受 体 及 其 在 蜕 皮 中 的 作 用 。 FEBS J2009;276：6128-57。[15] [10]杨文军，李文军. 20-羟基蜕皮激素和二苯甲酰肼结合的新模型：一同源性建模和对接方法。 蛋白质科学2000;9：1073-84.[16] Tohidi EsfahaniD，Lawrence MC，Graham LD，Hannan GN，Simpson AM，Hill RJ. 异翅类稻绿象甲蜕皮激素受体的亚型：蛋白质表征，RH5992杀虫剂结合和同源性建模。害虫管理科学2011;67：1457-67.[17] Bowers KJ，Chow E，Xu H，Dror RO ，Eastwood MP，Gregersen BA et al.Scalablealgorithmsformoleculardynamicssimulationsoncommodityclusters ， ProceedingsoftheACM/IEEEconferenceonsupercomputing（SC06），Tampa，Florida; 2006. p. 11比17[18] 薛定谔发布，Maestro-Desmond互操作性工具，版本3.6，薛定谔。D.E. ShawResearch，New York; 2013.[19] ShivakumarD，Williams J，Wu Y，Damm W，Shelley J，Sherman W. 用分子动力学自由能微扰和OPLS力场预测绝对溶剂化自由能。J Chem Theory Comput2010;6：1509-19.[20] Nauen R，Smagghe G. Etoxazole的作用方式。害虫管理科学2004;62：379-82.[21] 李英，杨宁，魏X，凌Y，杨X，王Q。乙螨唑对我国蔬菜害虫和螨类的药效评价。昆虫科学杂志2014;14：1-14.[22] 丁虫酰肼：一种具有异常高靶标选择性的新型毛虫防治剂。绿色化学合成与工艺。见：ACS系列研讨会，宾夕法尼亚州; 2000年。p. 八比十七[23] 杨文，李文.蜕皮激素受体配体结合口袋的结构适应性。Lett Nat 2003;426：91-6.[24] Wing KD，Slawecki RA，Carlson GR. RH-5849，非甾体蜕皮激素激动剂：对鳞翅目幼虫的影响。科学1988;241：470-2.[25] Roy A，Saraf S. Limestone：植物界重要生物活性三萜类化合物概述。《生物医药通报》 2006;29：191-201。[26] Cohen E，Quistad GB，Casida JE.印楝杀虫剂中印楝内酯、环氧氮杂二酮及其它柠檬苦素类化合物的细胞毒性。 生命科学1996;58：1075-81.[27] Dinan L，Whiting P，Girault JP，Lafont R，Dhadialla TS，Cress DE，等. 葫芦素是作用于蜕皮激素受体的昆虫类固醇激素拮抗剂。生物化学杂志1997;327：643-50.