نوع مقاله: مقاله پژوهشی

نویسندگان

1 گروه زیست شناسی، دانشکده علوم، دانشگاه فردوسی مشهد، مشهد، ایران

2 - گروه سلولی و مولکولی، پژوهشکده فناوری زیستی، دانشگاه فردوسی مشهد، مشهد، ایران 3- گروه زیست شناسی، دانشکده علوم، دانشگاه فردوسی مشهد، مشهد، ایران

3 گروه شیمی، دانشکده علوم، دانشگاه فردوسی مشهد، مشهد، ایران

چکیده

در خصوص ویژگی‌های ساختاری و عملکردی مجموعه سموم آفت کش Xpt اطلاعات کمی موجود است. آگاهی از دمین‌های عملکردی و اتصالی در جهت ایجاد سازه‌های ژنی جدید با حذف نواحی غیر‌ضروری یا قرار دادن دمین‌های عملکردی در مجاورت پروتئین‌های دیگر و تولید سموم چند‌‌توان حائز اهمیت است. در پژوهش حاضر محاسبات زیستی این مجموعه جهت آگاهی از ویژگی‌های نواحی مختلف آن انجام گردید. توالی ژنی و پروتئینی این مجموعه شامل زیر واحد‌های XptA1 ، XptA2 ، XptB1 و XptC1 از پایگاه‌های داده استحصال و بررسی ساختاری و عملکردی شامل خصوصیات فیزیکوشیمیایی، توپولوژی، دامنه میزبانی، ساختارهای ثانویه و سه‌بعدی توسط برنامه‌های رایج محاسباتی صورت پذیرفت. قابلیت ترشح مستقل از سیگنال پپتید، دمین‌های مرتبط با خاصیت حشره‌کشی و عملکرد سایر ‌دومین ها آشکار گردید. همگون‌یابی توالی نواحی عملکردی منجر به ردیابی آن ها در طیف وسیع از باکتری‌های گرم منفی شد. ساختار ثانویه نشان‌ دهنده ماهیت هلیکسی ناحیه عملکردی XptA1 و XptA2 و صفحه‌ای XptC1 است. قابلیت عملکردی مطلوب XptC1 نسبت به سایر اجزا در اتصال به اکتین نشان داده شد. نتایج حاصل از این تحقیق با آشکار‌سازی اطلاعات پایه‌ای نواحی مختلف مجموعه سموم Xpt امکان طراحی سازه‌های ژنی جدید مبتنی بر زیست شناسی سنتزی را فراهم می‌آورد

کلیدواژه‌ها

موضوعات

عنوان مقاله [English]

Assessment of structural and functional properties of Xpt complex protein to develop molecular approach in design of novel generation of pesticide

نویسندگان [English]

  • Mahsa Jalili-Manesh 1
  • Aliakbar Haddad-Mashadrizeh 2
  • Ali Makhdoumi 1
  • Mohammad Reaz Housaindokht 3

1 Department of Biology, Faculty of Science, Ferdowsi University of Mashhad, Mashhad, Iran

2 2- Cell and Molecular Biotechnology Research Group, Institute of Biotechnology, Ferdowsi University of Mashhad, Mashhad, Iran, 3- Department of Biology, Faculty of Science, Ferdowsi University of Mashhad, Mashhad, Iran

3 Department of Chemistry, Faculty of Science, Ferdowsi University of Mashhad, Mashhad, Iran

چکیده [English]

Little information is available about structural and functional characteristics of insecticidal Xpt protein complex. Awareness of the functional and binding domains is important to create new and more efficient gene constructs by eliminating unessential areas or put up functional domains adjacent to other proteins and produce new multi-functional toxins. Therefore, in this study the biological computing of this complex was done in order to acquiring data about the characteristics of its different areas. After extracting of this complex gene and protein sequences consist of XptA1, XptA2,XptB1 and XptC1 subunits from databases, structural and functional monitoring include physicochemical properties, topology, host range, secondary and three-dimensional structures was done by common softwares. The results reveal the secretion independent of signal peptides, domains related to insecticidal properties and other functions. Also functional domains BLAST leading to tracing broad spectrum of gram-negative bacteria with these sequences that up to now their ability in biocontrol not to be exist. Secondary structure indicates helix nature of functional area in XptA1 and XptA2 and strand of XptC1. Design and evaluation of three- dimensional structures report acceptable quality of all subunits except XptB1. In addition, the desirable functionality of XptC1 compared with other subunits in connection to actin was obtained. The results of this study reveal basic information of different areas of Xpt toxin complex that provides designing of new gene constructs based on synthetic biology.

کلیدواژه‌ها [English]

  • Xpt protein
  • Biocontrol
  • Biological computing

ارزیابی ویژگیهای ساختاری وعملکردی مجموعه پروتئینی Xpt به منظور توسعه‌ روشهای مولکولی در طراحی  نسل جدید آفت کشها

مهسا جلیلی منش1، علی اکبر حدادمشهدریزه1، 2*، علی مخدومی1 و محمدرضا حسین‌دخت3

1 مشهد، دانشگاه فردوسی مشهد، دانشکده علوم، گروه زیست شناسی

2 مشهد، دانشگاه فردوسی مشهد، پژوهشکده فناوری زیستی، گروه سلولی و مولکولی

3 مشهد، دانشگاه فردوسی مشهد، دانشکده علوم، گروه شیمی

تاریخ دریافت: 7/7/96                  تاریخ پذیرش: 16/11/96

چکیده

در خصوص ویژگیهای ساختاری و عملکردی مجموعه سموم آفت کش Xpt اطلاعات کمی موجود است. آگاهی از دمینهای عملکردی و اتصالی در جهت ایجاد سازه‌های ژنی جدید با حذف نواحی غیر‌ضروری یا قرار دادن دمینهای عملکردی در مجاورت پروتئینهای دیگر و تولید سموم چند‌‌توان حائز اهمیت است. در پژوهش حاضر محاسبات زیستی این مجموعه جهت آگاهی از ویژگیهای نواحی مختلف آن انجام گردید. توالی ژنی و پروتئینی این مجموعه شامل زیر واحد‌های XptA1، XptA2، XptB1 و XptC1 از پایگاههای داده استحصال و بررسی ساختاری و عملکردی شامل خصوصیات فیزیکوشیمیایی، توپولوژی، دامنه میزبانی، ساختارهای ثانویه و سه‌بعدی توسط برنامه‌های رایج محاسباتی صورت پذیرفت. قابلیت ترشح مستقل از سیگنال پپتید، دمینهای مرتبط با خاصیت حشره‌کشی و عملکرد سایر دمینها آشکار گردید. واکاوی توالی نواحی عملکردی منجر به آشکارسازی آنها در طیف وسیعی از باکتریهای گرم منفی شد. ساختار ثانویه نشان‌ دهنده ماهیت هلیکسی ناحیه عملکردی XptA1 و XptA2 و صفحه‌ای XptC1 است.  قابلیت عملکردی مطلوب XptC1 نسبت به سایر اجزاء در اتصال به اکتین نشان داده شد. به طورکلی، نتایج حاصل از این تحقیق با آشکار‌سازی اطلاعات پایه‌ای نواحی مختلف مجموعه سموم Xpt امکان طراحی سازه‌های ژنی جدید مبتنی بر زیست شناسی سنتزی را فراهم می‌آورد.

واژه های کلیدی: پروتئین Xpt، کنترل زیستی، محاسبات زیستی

* نویسنده مسئول، تلفن: 09155187002 ، پست الکترونیکی: a.haddad@um.ac.ir

مقدمه

 

اهمیت دو موضوع اساسی امنیت و ایمنی غذایی ضرورت روزافزون کنترل آفات کشاورزی به ویژه حشرات به عنوان تهدید عمده در این حوزه را باعث شده است (3، 17، 20 و 22). تاکنون روشهای متعددی برای نیل به این اهداف مورد استفاده قرار گرفته که در این میان کنترل زیستی با تکیه بر باکتری Bacillus thurengiensis و محصولات مشتق شده از آن یکی از گسترده‌ترین و کم مخاطره‌ترین روشها برای محیط زیست و بروز مقاومت در حشرات آفت قلمداد می‌شود (7، 12 و 24) . با این حال، امروزه به دلیل وجود و بروز برخی کاستیها (8 و 17) در این نوع حشره‌کش تلاشهایی با هدف شناسایی، تجاری‌سازی و استفاده از ذخایر جدیدی از باکتریها و سموم پروتئینی حشره‌کشی آنها در حال انجام است، که منجر به شناسایی چندین باکتری جدید دارای خاصیت حشره‌کشی شده است. از جمله این باکتریها می‌توان به گونه‌های  Xenorhabdus اشاره کرد. این باکتری گرم منفی، متعلق به خانواده انتروباکتریاسه، بی‌هوازی اختیاری و همزیست با نماتود‌های جنس Steinernema می‌باشد (6 و 13). روابط همزیستی این باکتری‌ با نماتود به موازات مزایایی نظیر محدوده وسیع‌تر طیف اثر، کارآیی بالا، عملکرد سریع و سرعت بالا در ایجاد مرگ (23)، عملکرد اختصاصی و نیز عدم نیاز به پوشش محافظ در حین کار نسبت به Bt (Bacillus thurengiensis) و سایر حشره‌کشها، به دلیل مواردی مانند وابستگی شدید و عدم توانایی برای زیست مستقل در آب و خاک محدودیت استفاده از این باکتریها به عنوان حشره‌کش را در پی داشته است (8). در این راستا پژوهشهایی با هدف شناسایی ژنها و آنزیمهای دارای خاصیت حشره‌کشی در این باکتری منجر به شناسایی یک مجموعه پروتئینی مهم به نام Xpt (Xenorhabdus protein toxin) شده است (19، 26 و 29). این مجموعه پروتئینی متشکل از چهار زیر واحد XptA1، XptA2، XptB1 و XptC1 می‌باشد و شباهت ساختاری و عملکردی بالایی نسبت به مجموعه پروتئینی Tc (Toxin complex) در باکتری Photorhabdus به سبب ماندگاری، مشاهده خواص حشره کشی ژنها پس از کلون در E. coli دارند (21). بر پایه بیان متفاوت کمی و کیفی ژنهای xpt اختصاصیت حشره‌کشی و دز کشندگی مختلف در لاروهای Lepidoptera مشاهده شده است (7، 17 و 27). در پژوهشهای اخیر ترکیبات فعال زیستی در این مجموعه‌های سمی با قابلیت اثر بر روی اسکلت سلولی اکتین شناسایی شده که به دلیل مکانیسم عمل و اندامک هدف احتمالی متفاوت با Bt توانسته آن را به گزینه مناسب نسل آینده حشره‌کشهای میکروبی بدل نماید (24 و 25).زیر واحد XptA1 چهار بخشی قرینه با ساختاری شبیه قفس یا بطری مانند است. هر کدام از بخشها دارای سه دمین و یک پیچ‌خوردگی طولی با انتهایی باریک‌تر نسبت به انتهای دیگر دارد. پیش‌بینی محتوی ساختار ثانویه پروتئین حضور مارپیچها و کویلها را در بخش عمده توالی و نسبت کمتر صفحات را نشان داده است (19). در خصوص ساختار ثانویه یا سه بعدی سایر زیر واحدها اطلاعاتی گزارش نشده است. با این حال به دلیل قدمت اندک شناسایی این باکتری و توکسین مربوطه همچنان اطلاعات اندکی در خصوص خواص این مجموعه پروتئینی و مکانیسم دقیق آن و گیرنده سطح سلولی آن در دست است. امروزه از راه‌کار‌های مؤثر جهت افزایش کارآیی پروتئینهای نوترکیب با قابلیت عملکرد اختصاصی و نیز افزایش عملکرد آنها، ادغام پروتئینهای عملگرا و یا دمینهای عملکردی مستتر در این پروتئینهاست. در این راستا، در حوزه مبارزه با آفات افزایش کارآیی با افزایش دامنه اثر از طریق قرار دادن دمین اتصالی مناسب در کنار دمینهای کاتالیتیکی و افزایش قدرت آفت کشی با کنار هم قرار دادن دمینهای کاتالیتیکی سموم مختلف می‌تواند مورد بهره‌برداری قرار گیرد که خود نیازمند دسترسی به  اطلاعاتی کامل در مورد پروتئینها می‌باشد (25). لذا با هدف فراهم آوردن و افزایش اطلاعات پایه مولکولی، شناسایی و تفکیک نواحی دارای عملکرد در این مجموعه پروتئینی، ردیابی میزبانهای جدید حامل پروتئینهای مشابه به عنوان گزینه‌های احتمالی جهت تولید سموم کاراتر، پیش‌بینی و تعیین ساختار دوم و سه بعدی به منظور پیش‌نیاز توسعه و قدرتمند‌سازی بیشتر نسل جدید حشره‌کشهای زیستی این تحقیق به آن پرداخته است.

مواد و روشها

واکاوی داده‌ها: با هدف شناسایی اعضای مجموعه پروتئینی سمی Xpt بررسی منابع و مقالات انجام شد که منجر به دستیابی چهار عضو XptA1، XptA2، XptB1 و XptC1 در این خانواده پروتئینی شد.

استحصال توالیهای نوکلئوتیدی، پروتئینی و ساختار سه بعدی: دستیابی به توالی نوکلئوتیدی و پروتئینی ژنها و پروتئینهای مورد نظر به ترتیب با شماره شناسایی AJ308438 و شماره‌های دستیابی CAC38401، CAC38404، CAC38402 و CAC38403 از پایگاه اطلاعات نوکلئوتیدی و پروتئینی NCBI به آدرس http://www.ncbi.nlm.nih.gov)) میسر شد.

ساختارهای سه بعدی با کد 1vw1 جهت زیر واحد XptA1 و XptA2، 4igl و 4o9x به ترتیب برای زیر واحدهای XptB1 و XptC1 به عنوان ساختارهای الگو جهت مدل‌سازی سه بعدی از پایگاه RCSB (http://www.rcsb.org/pdb/home/home.do)   استحصال شدند.

از سویی دیگر، توالی پروتئینی اکتین حشره Plutella xylostalla (یکی از مخرب‌ترین حشرات آفت راسته Lepidoptera در حوزه کشاورزی و به نمایندگی از سایر آفات این گروه) به عنوان اندامک سلولی هدف احتمالی این کمپلکس پروتئینی سمی با شماره شناسایی AEP25396 از پایگاه داده NCBI استحصال و مشابه‌ترین ساختار سه‌بعدی آن به صورت مونومری و بدون هیچ‌گونه لیگاندی از پایگاه Swissmodel به آدرس (http://www.swissmodel.expasy.org)  تحت فرمت pdb و با کد 3b5u.1.N ذخیره شد.

بررسی ویژگیهای ساختاری و عملکردی توالیهای انتخابی: بررسی توالی نوکلئوتیدی از لحاظ طول و محتوی GC به ترتیب توسط پایگاه اطلاعات نوکلئوتید به آدرس (http://www.ncbi.nlm.nih.gov/nuccore/?term) و Endmemo  (http://www.endmemo.com/bio) صورت پذیرفت. بررسی ویژگیهای ساختاری توالیهای پروتئینی انتخابی با استفاده از پایگاه پروتئینی  NCBI (http://www.ncbi.nlm.nih.gov/protein/term)، برنامه Protparam موجود در پایگاه Expasy به آدرس (http://web.expasy.org/protparam) جهت بررسی ویژگیهای فیزیکوشیمیایی، برنامه TMHMM (http://www.cbs.dtu.dk/services/TMHMM-2.0) با هدف دست‌یابی به توپولوژی ساختاری زیر واحد‌ها، برنامهSignal peptide prediction با آدرس (http://sigpep.services.came.sbg.ac.at) با هدف آشکارسازی توالی‌ پپتید‌های ترشحی و برنامه NetNGlyc، NetPhos و Sulfinator با آدرسهای (http://www.cbs.dtu.dk/services/NetNGlyc)، (http://www.cbs.dtu.dk/services/NetPhos) و (http://web.expasy.org/Sulfinator) در جهت بررسی وجود تغییرات پس از ترجمه در توالیها صورت پذیرفت. آشکارسازی دمینهای عملکردی پروتئینهای این خانواده در سه پایگاه ConservedDomain ، Motifscan و InterProScan 5 به ترتیب با آدرسهای (http://www.ncbi.nlm.nih.gov/cdd/)، (http://myhits.isb-sib.ch/cgi-bin/motif_scan) و (http://www.ebi.ac.uk/interpro/sequence-search) به صورت مجزا انجام گرفت.

بررسی دامنه میزبانی و قرابت‌یابی توالیها: برنامه تحت شبکه Protein Blast با آدرس (http:blast.ncbi.nlm.nih.gov/Blast) مبتنی بر ماتریکس PAM250 جهت آشکارسازی توالیهای پروتئینی مشابه استفاده شد. نتایج مناسب بر اساس معیار‌هایی نظیر ارزشE، درصد همپوشانی و همسانی انتخاب، خوشه‌بندی نتایج براساس الگوریتم Neighbor joining برنامه MEGA6 انجام و در نهایت قرابت توالیها نسبت به یکدیگر تحت درخت فیلوژنیک مشخص شد.

بررسی ساختار ثانویه توالیهای پروتئینی: به منظور پیش‌بینی و جمع‌آوری اطلاعات در خصوص وجود، میزان و موقعیت قرارگیری اجزای ساختار ثانویه مانند هلیکسها و صفحات و همچنین حضور یا عدم حضور ساختار Coiled-coil به ترتیب از برنامه‌های تحت شبکه Psipred با روشهای Psipred و Dompered، Paircoil2 با تنظیمات پیش‌فرض و COILs تحت ماتریکس MTIDK به آدرسهای (http://bioinf.cs.ucl.ac.uk/psipred)، (http://paircoil2.csail.mit.edu/paircoil2.html) و (http://www.ch.embnet.org/software/COILS) استفاده شد.

تعیین ساختار سه‌بعدی توالیهای پروتئینی و ارزیابی کیفیت مدلهای طراحی شده: ساختار 3D زیر واحدهای مورد بررسی توسط نرم‌افزار MODELLER9.15 تعیین شدند. به این منظور، انتخاب ساختار الگو توسط روش همگون‌یابی در پایگاه داده Pdb صورت پذیرفت، سپس بهترین نتیجه از لحاظ همسانی و همپوشانی انتخاب و از پایگاه PDB به آدرس (http://www.rcsb.org/pdb/home/home.do) ساختار سه بعدی به فرمت Pdb ذخیره شد. مدل‌سازی بر پایه همردیف‌سازی بین توالی الگو و هدف انجام، انرژی مدلهای طراحی شده براساس معیار dope score توسط نرم‌افزار محاسبه شدند و بهترین مدل براساس پایین‌ترین dope score انتخاب و سپس بهینه‌سازی ساختار‌های طراحی شده توسط نرم‌افزار MOE صورت پذیرفت. اعتبار‌سنجی مدلها براساس ارزیابی شیمی فضایی (28) با استفاده از نرم‌افزار تحت وب ProSA (https://prosa.services.came.sbg.ac.at/prosa.php) (30)و نمودار راماچاندران برنامه MOE به صورت مجزا انجام شد.

ارزیابی میل اتصال پروتئینهای انتخابی به پروتئین اکتین: بر‌همکنش مولکولی و میل اتصال پروتئینهای سمی انتخابی به هدف درون سلولی احتمالی اکتین بر پایه اصول اشکال مکمل توسط برنامه تحت وب Patchdock (11) به آدرس (http://bioinfo3d.cs.tau.ac.il/PatchDock) انجام شد. نتایج براساس انرژی واکنش‌پذیری اتمی(ACE) از میان 20 راه حل اول با بیشترین اعتبار مورد بررسی قرار گرفتند و مناسب‌ترین نتایج گزینش شدند. آشکارسازی صحت اتصال در موقعیت مناسب با استفاده از برنامه Pymol1 صورت پذیرفت و توسط تصاویر حاصل از این نرم‌ افزار نمایش داده شد.

نتایج

ویژگیهای ساختاری و عملکردی توالیهای مورد بررسی: بررسی ساختاری توالیهای موردنظر منجربه آشکارسازی اطلاعاتی در ارتباط با طول ژن و پروتئین، درصد GC ژنها و نیز تغییرات پس از ترجمه‌ای متفاوت این پروتئینها شد (جدول 1).

 

جدول 1- ویژگیهای ساختاری و تغییرات پس از ترجمه توالیهای ژنی و پروتئینی

ردیف

نام زیرواحد

طول ژن

درصدGC

طول پروتئین(aa)

سیگنال پپتید

گلیکوزیلاسیون/فسفریلاسیون/سولفوریلاسیون

1

XptA1

7572

88/40

2523

+

+/+/+

2

XptA2

7617

16/49

2538

-

+/+/+

3

XptB1

3045

14/50

1014

-

+/+/-

4

XptC1

4206

38/50

1401

-

+/+/+

 


مرور نتایج حاصله طویل و بزرگ بودن ژنها و پروتئینها و تشابه در محتوی GC ژنها را نسبت به یکدیگر آشکار نمود. نتایج ارزیابی تغییرات پس از ترجمه شامل گلیکوزیلاسیون، فسفوریلاسیون، سولفوریلاسیون هر زیر واحد به صورت مجزا منجر به آشکارسازی این مدیفیکاسیونها در تمامی زیر واحدها شد (جدول 2).

ویژگیهای فیزیکوشیمیایی توالیهای پروتئینی: ارزیابی در این حوزه منجر به آشکار‌سازی وزن مولکولی بالای پروتئینها در محدوده 110 تا 287 کیلودالتونی، بار الکتریکی 1- تا 74-، pH ایزوالکتریک (pI) 5 تا 7، هیدروفوبیسیته 46/0- تا 276/0-، نیمه عمر بیش از 10 ساعت و شاخص ناپایداری حدود 30 تا 40 و شاخص آلیفاتیک ( شاخص آلیفاتیک حجم نسبی زنجیره آلیفاتیک در اسید‌های آمینه آلانین، والین، لوسین و ایزولوسین یک پروتئین می‌باشد که اثری مثبت در افزایش مقاومت حرارتی پروتئینهای کروی دارا می‌باشد.)  75 تا 90 شده است (جدول2). همان‌طور که در این جدول نشان داده شده زیرواحد XptA1 بیشترین وزن مولکولی، شاخص آلیفاتیک، بار منفی و شاخص ناپایداری بالای 40 را به خود اختصاص داده که ناپایداری آن را نسبت به سایر زیرواحدها نشان می‌دهد.


جدول 2- مشخصه‌های فیزیکوشیمیایی توالیهای پروتئینی زیرواحد‌ها

ردیف

نام پروتئین

وزن مولکولی (دالتون)

بار الکتریکی

pI

هیدروفوبیسیتی

نیمه عمر (ساعت)

شاخص ناپایداری

شاخص الیفاتیک

1

XptA1

9/287002

74-

23/5

319/0-

>10

02/40

03/90

2

XptA2

6/283997

47-

19/5

276/0-

>10

28/36

40/86

3

XptB1

4/110257

1-

05/7

336/0-

>10

52/30

36/82

4

XptC1

4/157736

40-

17/5

463/0-

>10

74/39

18/75

 


ویژگیهای توپولوژی توالیهای مورد بررسی: بررسی و پیش‌بینی موقعیت سلولی توالیهای مورد‌ نظر موقعیت خارج سلولی و ترشحی بودن زیرواحد‌های XptA1، XptA2 و XptC1 را آشکار کرد. با این وجود، زیر واحد XptB1 نتایج متفاوتی را نشان داد، به گونه‌ای که 2 دمین گذرنده از غشاء در موقعیت اسید‌های آمینه 620 تا 646 و 709 تا 1008 و همچنین نواحی داخل سلولی از محدوده اسید‌آمینه 642 تا 726 و 736 تا 1014( انتهای پروتئین) در آن آشکار شد. از سویی ارزیابی ماهیت ترشحی بودن کلیه زیرواحدها منجر به آشکارسازی تنها یک توالی احتمالی موردنظر از اسیدآمینه 71 تا 87 با درصد همسانی 52 و امتیاز 6/21 در زیرواحد XptA1 شد. بررسیها در این حوزه در سایر نقاط توالی XptA1 و سایر زیرواحدها نتیجه‌ای در پی نداشت و سیگنال پپتید احتمالی دیگری ردیابی نشد.

بررسی عملکردی توالیهای پروتئینی: این بررسیها منجر به آشکارسازی چندین دمین عملکردی با خاصیت حشره‌کشی و با طولهای مختلف در زیرواحد‌های XptA1، XptA2 و XptC1 و نیز دمینهایی با عملکردهای غیرمرتبط با خاصیت حشره‌کشی در تمامی زیرواحد‌ها و چندین دمین با عملکرد ناشناخته در زیرواحد XptB1 شد (جدول 3).

 

جدول 3-ویژگیهای عملکردی توالی پروتئینی زیرواحدهای مورد بررسی

ردیف

نام پروتئین

نام دمین

موقعیت

طول

عملکرد

1

XptA1

Insecticidal_toxin/plasmid_vir

(Salmonella virulence plasmid 28.1kDa A protein)

398-1

398

تهاجم موفق به سلول میزبان

Hydrogenase expression/synthesis hypA family

2087-2074

14

متالوچپرون (ورود اتم نیکل به هیدروژنازها)

2

XptA2

Insecticidal_toxin/plasmid_vir

(Salmonella virulence plasmid 28.1kDa A protein)

291-8

284

تهاجم موفق به سلول میزبان

Cu-oxidase_4
Multi-copper polyphenol oxidoreductase laccase

1335-1315

21

تجزیه ترکیبات فنلی و آروماتیک و برخی ترکیبات غیر فنلی

Flu_NS2
Influenza non-structural protein (NS2)

2249-2234

16

واسطه خروج هسته ای ریبونوکلئوپروتئینی ویروسی

3

XptB1

Rhs_assc_core super family

679-652

28

عملکرد نامشخص

DUF456

826-725

102

عملکرد نامشخص

RhsA

677-7

671

عملکرد نامشخص

FlgD
Flagellar hook capping protein

26-7

20

پروتئین درپوش قلاب فلاژلی(اسمبل قلاب)

Borrelia_rep
Borrelia repeat protein

197-180

18

فاقد عملکرد خاص

Rick_17kDa_Anti
Rickettsia 17 kDa surface antigen

903-860

44

آنتی ژن سطحی ریکتزیا

4

XptC1

SpvB

Salmonella virulence plasmid 65kDa B protein

331-28

304

ریبوزیله کردن اکتین(اثر بر اسکلت سلولی)

TcdB_toxin_midN

835-655

181

مرتبط با توالی تکراری Rhs(نقش احتمالی در اتصال به لیگاند)

TcdB_toxin_midC

1031-884

148

ناحیه میانی انتهای C سم حشره کش TcdB

VCBS repeat domain

535-474

62

نقش احتمالی در چسبندگی

FG-GAP repeat

493-466

28

توالی تکراری

544-516

29

 


پراکنش دامنه میزبانی زیرواحد‌ها و تعیین قرابت توالیهای پروتئینی: نتایج همگون‌یابی توالیهای پروتئینی در هر زیرواحد منجر به ردیابی تعداد زیادی همولوگ در محدوده همسانی و همپوشانی‌ مختلف (همسانی و همپوشانی 40 تا 100 درصد در زیر واحد XptA1، همسانی 34 تا 100 درصد و همپوشانی بالای 90 درصد برای زیر واحد XptA2، همسانی 47 تا 100 درصد  و همپوشانی 65 درصد به بالا در XptB1، همسانی 43 تا 100 درصد و همپوشانی 96 تا 100 درصد در زیر واحد XptC1 ) و ارزش E صفر در تمامی زیر واحد‌ها شد. توالیهای مشابه یافت شده در تمامی زیرواحدها منحصراً متعلق به باکتریهای گرم منفی بوده و هیچ موردی از آنها در باکتریهای گرم مثبت آشکار نشد. همچنین نتایج مبین حوزه وسیع‌تر پراکنش میزبانی توالیهای مورد بررسی نسبت به گزارشهایی که تاکنون ارائه شده‌اند می‌باشد، به گونه‌ای که چندین جنس و گونه جدید دارای توالیهای همولوگ با همسانی و همپوشانی مطلوب و ارزش E مناسب در باکتریهای Burkholderia pseudomallei، Trabulsiella odontotermitis، Erwinia toletana،Candidatus Paraburkholderia kirkii، Morganella morganii، Kosakonia sacchari،Cedecea neteri، Enterobacter cloacae، Shewanella baltica، Kosakonia sacchari SP1،Vibrio campbellii و Mortierella elongata بودند. قرابت‌یابی توالیهای آشکارشده در تمامی زیرواحد‌ها با همسانی بالای 50 درصد، منجر به قرارگیری گونه‌های مختلفی از Xenorhabdus و Erwinia و Photorhabdus با توالیهای مورد نظر در یک خوشه شد (شکل 1).

ساختار دوم توالیهای پروتئینی: محتوی ساختار دوم زیرواحدهایXptA1 و XptA2 منجربه آشکارسازی عمده ساختار آنها به صورت هلیکس با توزیع و تراکم بیشتر در یک سوم ابتدایی و انتهایی توالی شد.

 

 

شکل 1- قرابت توالیهای پروتئینی بررسی شده نسبت به یکدیگر و همولوگ‌های یافت شده (علائم ( بولتها )بیانگر توالیهای مورد بررسی است و شاخه‌های رنگی مبین نزدیک‌ترین همگونها به پروتئینهای بررسی شده است (آبی: زیرواحد XptB1- سبز : زیرواحد XptC1- صورتی : XptA2و قرمز : XptA1).

 

در همین راستا در زیر واحد‌های XptA1 و XptA2 به ترتیب تنها 9 و 7 درصد توالیها به صفحات بتا اختصاص یافت که بیشترین تراکم این ساختار در قسمت میانی در دو زیرواحد بود. بررسی ساختار دوم زیر واحد‌های XptB1 و XptC1 نتایج متفاوتی را آشکار نمود، به‌طوری‌که مقادیر صفحات بتا در زیرواحد XptB1 نسبت به هلیکس بیشتر بوده (5/10 درصد) که تراکم آن در نیمه ابتدایی بیشتر می‌باشد. همچنین مقادیر پیش‌بینی شده برای ساختار هلیکس زیرواحد XptC1 بسیار اندک در قیاس با سایر زیرواحد‌ها بود (5/5 درصد) که در نیمه انتهایی توالی مقدار این ساختار بیشتر است. با این حال توزیع صفحات بتا در این زیرواحد نسبت به سایرین یکنواخت‌تر پیش‌بینی شده است. از سویی دیگر، نواحی پیش‌بینی شده حضور ساختار Coiled-coil در زیرواحد‌ XptA1 در انتهای C با موقعیت 2240-2196 و موقعیتهای 2136-2103، 2256-2211 و 2273-2257 برای XptA2 ردیابی شد.

مدل‌سازی سه بعدی و ارزیابی کیفیت ساختار‌های سه بعدی: مقادیر Dope score برای زیرواحدهای XptA1، XptA2، XptB1 و XptC1 به ترتیب 81250/271933-، 71875/269717-، 06600/64722- و 87500/142070- حاصل شد. ارزیابی کیفیت کلی مدلها توسط نرم‌افزار تحت وب ProSA مقادیر امتیاز Z برای زیرواحد‌های A1، A2، B1 و C1 را به ترتیب 15/17-، 16/18-، 44/2- و 13/12- گزارش کرد. این مقادیر برای زیر واحد‌های A1، A2 و B1 علی رغم منفی بودن خارج از محدوه حالت طبیعی به دست آمده توسط روش کریستالوگرافی و NMR مشاهده شدند. منحصراً در زیر واحد C1 این مقدار در محدوده طبیعی و مطلوب واقع شده است ( شکل 2).

 

شکل2- نمودار انرژی کلی ساختار‌های طراحی شده (نقطه مشکی در هر نمودار امتیاز به دست آمده و جایگاه قرارگیری در محدوده ساختار طبیعی و غیرطبیعی برای زیرواحد موردنظر است.)

در نمودارهای انرژی محلی (شکل 3) انرژی ریشه‌های آمینواسید اغلب نواحی در طول توالی پروتئینی در زیرواحد A1 و A2 مقادیر زیر صفر و منفی را به خود اختصاص داده‌اند و تنها تعداد کمی پیک با انرژی بالاتر از صفر در قسمت میانی پروتئینها مشاهده شده‌ است. در خصوص زیر واحد B1 حدود 90 درصد شاهد مقادیر انرژی بالاتر از صفر بوده و در نمودار زیر واحد C1 میزان انرژی بیشتر در محدوده منفی و مطلوب قرار دارد (شکل 3) . تطبیق این نمودارها با ساختار سه بعدی نمایش داده شده توسط نرم‌افزار Jsmol در شکل 3 ریشه‌های دارای انرژی بالای صفر را به رنگ قرمز نشان داده در حالی که اسید‌های آمینه با کمترین میزان انرژی و پایدارترین حالت به رنگ آبی نشان داده شده‌اند.

بررسی نمودار‌های راماچاندران به دست آمده توسط نرم‌افزار MOE نشان دهنده قرارگیری 42، 37، 35 و 26 اسید‌آمینه خارج از محدوده مجاز به ترتیب در پروتئینهای XptA1، XptA2، XptB1و XptC1است، که موقعیت و نوع این اسیدهای آمینه تحت جدول 4 تنظیم شده است.

ارزیابی میل اتصالی توکسین‌های انتخابی به اکتین: نتایج و بررسی میل اتصال و محل بر‌همکنش اکتین به پروتئینهای طراحی شده مبین وجود چندین ناحیه بر‌همکنشی دارای امتیاز‌ها و انرژی واکنش‌پذیری اتمی مختلف و قابل قبول می‌باشند. از این رو مطلوب‌ترین نتیجه از بین 20 راه‌حل اولیه در زیر‌واحدهای XptA1 ، XptA2، XptB1 و XptC1 به ترتیب دارایACE و امتیاز 12/200- و 14848، 16/214-و 17112، 53/407- و 23202 و 24/527- و 17784 حاصل شد (شکل 4).  همان‌طور که در نمودار مشهود است زیر واحد XptC1 محکمترین اتصال را با اکتین برقرار نموده است. بر‌همکنش و اتصال اکتین با زیر واحد‌ها در شکل 5 نشان داده شده است.

 

جدول 4- موقعیت و نوع اسیدهای آمینه غیر مجاز بر اساس نمودار راماچاندران(موارد قرمز اسید‌های آمینه غیر‌مجاز در موقعیت دمین عملکردی است.)

ردیف

1

2

3

4

5

6

7

8

9

XptA1

ASN216

LEU217

SER299

GLU352

GLU365

ASN375

LYS378

ASN420

PRO421

SER422

ASP537

ASN621

PHE730

LYS837

MET899

LYS995

ALA1126

GLU1170

THR1383

ASN1394

PHE1398

ASP1440

CYS1482

TYR1548

ILE1550

TYR1565

PHE1612

LYS1655

ASN1670

ILE1671

ASP1672

SER1675

ASP1744

LYS1756

LYS1759

MET1777

GLY1828

ASP1843

PRO1864

ALA2318

THR2361

ASN2464

 

 

 

 

XptA2

SER3

LYS326

TYR340

ASP341

ASN362

THR424

LYS626

ALA794

ASP800

GLU821

GLN863

SER907

GLU957

MET1055

LYS1154

ASP1312

SER1387

SER1406

ASP1407

TYR1421

LEU1462

LYS1478

PHE1494

SER1496

SER1522

SER1556

ASN1562

ASP1668

GLU1671

ASP1759

PHE1887

PRO1969

ALA2039

THR2332

SER2376

ASN2489

ASP2527

 

 

 

 

 

 

 

 

 

XptB1

ASN17

GLN35

THR175

LYS198

ASP199

SER201

SER221

PHE222

ALA278

HIS287

GLN302

TYR327

ASN350

VAL354

GLN366

THR397

ASP406

SER418

GLY458

THR489

LYS539

ILE545

ALA550

LYS619

LEU624

ARG720

THR724

ALA743

GLU775

SER806

GLU860

ARG861

 

 

 

 

 

XptC1

TYR90

THR120

ALA184

ASP219

ASP260

GLU293

ARG491

SER497

GLU501

ASP556

HIS559

LEU560

LYS563

HIS605

ALA638

ASP671

THR711

SER813

ASP989

ASP1022

LYS1138

THR1139

GLU1149

ASN1291

SER1334

HIS1395

 

 

شکل 3- نمودار انرژی محلی و ساختار سه بعدی توسط نرم افزار Jsmol (تصاویر به ترتیب از چپ به راست : XptA1، XptA2، XptB1 و XptC1)

 

شکل 4- مقایسه ارزش انرژی تماس اتمی پروتئینهای طراحی شده با اکتین براساس  ACE

 

شکل 5- بر‌همکنش و اتصال زیر واحد‌ها با مولکول اکتین ( پروتئینهای سمی با رنگ آبی و اکتین با رنگ سبز نشان داده شده است.)


بحث و نتیجه گیری

کنترل جمعیت حشرات آفت همواره برای بشر با هدف تأمین امنیت غذایی امری مهم تلقی شده است، از این رو حشره‌کشهای سنتزی به کار گرفته شدند (8). با این حال به دلیل اثرات مخرب این ترکیبات، امروزه روش کنترل زیستی به خصوص باکتریها و محصولات تولید شده توسط آنها به ویژه Bt جایگزین مناسب این گروه محسوب می‌شود (8و29). با این وجود، گزارش به روز مقاومت علیه این سموم تلاشها را در جهت شناسایی باکتریها و سموم جدید برای نسل بعدی حشره‌کشهای میکروبی برانگیخت، که منجر به شناسایی باکتریهای همزیست نماتود از جمله فتورابدوس و زنورابدوس شد. سموم Tc و Xpt این باکتریها از طریق سازوکار احتمالی متفاوت با Bt منجر به مرگ سریع لارو حشره می‌شود (17و18). در این راستا، مجموعه پروتئین سمی Xpt به دلیل کلون موفق  و حفظ عملکرد پس از کلون شدن نسبت به همتای خود توانسته توجهات را بیشتر به خود جلب کند (27). با این حال تاکنون پژوهشهای مولکولی اندکی در خصوص این مجموعه پروتئینی که به عنوان کاندید مناسب در کنترل زیستی محسوب می‌شود انجام شده است. بنابراین شناخت دقیق، بررسی ویژگیهای ساختاری و عملکردی و دستیابی به ساختار سه بعدی و جایگاههای فعال این مجموعه پروتئینی سمی طبق روشهای in silico ضمن آنکه می‌تواند راه‌کاری در جهت به کارگیری بهتر، هدفمندتر، سریع تر و کم‌هزینه‌تر این سموم در کنترل زیستی فراهم آورد، می‌تواند به ایجاد و توسعه باکتریهای تراریخت قدرتمندتر و همچنین گیاهان خود محافظ که بتوانند از سازوکارهای مختلف قابلیت حشره‌کشی داشته باشند منجر شود، که در این پژوهش در دستور کار قرار گرفت (1). مطالعه منابع مجموعه پروتئینی دارای 4 زیر واحد اصلی XptA1 ، XptA2 ، XptB1 و XptC1 را برای القا و ایجاد مرگ در حشرات معرفی کرده است (27). همچنین محتوی GC آنها در حد میانی و بسیار نزدیک به همتا خود یعنی توالی ژنی پروتئین Tc بود (8) (جدول 1). به تبع اندازه بزرگ ژن، وزن مولکولی بالای زیرواحدهای پروتئینی نیز مورد انتظار بود. پروتئین XptA1 به میزان جزئی دارای شاخص ناپایداری بالاتر از حد معمول می‌باشد. با این حال به دلیل دارا بودن بالاترین شاخص آلیفاتیک در بین سایر زیرواحدها می‌توان آن را پایدار در شرایط محیطی و دارای مقاومت دمایی مناسب در نظر گرفت. در خصوص سایر زیرواحدها این ویژگیها در محدوه مجاز و مناسب قرار دارد (جدول 2). pH ایزوالکتریک گزارش شده در جدول 2 مبین قرارگیری این معیار در محدوه اسیدی و خنثی در تمامی زیرواحدها می‌باشد، از این رو با نتایج پژوهش قبلی در مورد XptA1 مبنی بر پایداری ساختار ثانویه در محدوده pH خنثی و قلیایی (19) کمی متفاوت به نظر می‌رسد که می‌توان آن را از طریق تفاوت در توالیهای مورد بررسی در دو تحقیق و همچنین تفاوت در محدوده pH در مورد پایداری ساختاری پروتئین و بیشینه عملکرد و فعالیت پروتئین در pH خاص توجیه کرد. نتایج توپولوژی حاصل کاملاً تأیید کننده ماهیت خارج سلولی و غیر تجمعی این مجموعه پروتئینی در سلول مشابه گزارشهای قبل می‌باشد (4). با این حال به دلیل آشکار نشدن توالیهای سیگنال پپتید در تمامی زیرواحد‌ها به جز XptA1، دو فرضیه مطرح می‌شود که این پروتئینها همانند برخی از پروتئینهای تاکنون شناسایی شده در باکتریهای گرم منفی از مسیر و مکانیسم مستقل و غیر وابسته به سیگنال پپتید می‌توانند به بیرون از سلول ترشح شوند (9 و 16)، و یا اینکه توالی سیگنال پپتید در این زیرواحد‌ها وجود دارد اما تاکنون در پایگاههای داده سیگنال پپتید مشابه این توالیها شناسایی وثبت نشده است. وقوع تغییرات پس از ترجمه در تمامی زیرواحدهای پروتئینی (جدول 1) این اعضا را برای تولید و توسعه موجودات تراریخت یوکاریوتی به ویژه گیاهان‌ تراریخت و خود محافظ در جهت بهره‌وری بیشتر (21) بسیار مناسب و موفق گردانده است، که همین ویژگی می‌تواند مزیت، امتیاز مثبت و ارجحیت این کمپلکس پروتئینی سمی در مقابل دیگر سموم حشره‌کش کاندید شده به حساب آید. نتایج بررسی عملکردی، آشکار‌کننده چندین دمین مرتبط با فعالیت کشندگی سلول در تمامی زیرواحد‌ها به جز XptB1 است که کاملاً موافق با یافته‌‌های قبلی در مورد خاصیت حشره‌کشی زیرواحد‌ها به صورت مجزا و تلفیقی بر روی چندین گونه از حشرات آفت Lepidoptera با قدرت‌ و سرعت متفاوت می‌باشد (21 و 27). عملکرد مؤثر بر اسکلت سلولی در پروتئین Tc باکتری فتورابدوس (18) در این بررسی در زیرواحد XptC1 ردیابی شد. بخشی از توالی در پروتئینهای مورد بررسی فاقد نقش عملکردی مشخص و دمین شناخته شده هستند که دلیل این موضوع را می‌توان عدم اطلاع از نحوه عمل دقیق این مجموعه بر روی سلول هدف، ناشناخته بودن گیرنده و همچنین وجود برخی دمینها و موتیفهایی با توالی و عمل ناشناخته برشمرد. لذا به سبب پاسخ به سوالات و ابهامات موجود در خصوص مکانیسم اثر و گیرنده سلولی این اعضا پروتئینی مطالعات آزمایشگاهی و بیوانفورماتیکی بیشتر و پیشرفته‌تری در این زمینه احساس می‌شود. همگون‌یابی مشابه با روش شکری و همکاران انجام گردید و پراکنش میزبانی آشکار شده در این پژوهش علاوه بر تأیید حضور توالیهای همگون در باکتریهای مورد انتظار و معرفی شده در سایر تحقیقات مانند گونه‌های فتورابدوس و سراشیا با میزان تشابه‌ بسیار نزدیک به مقالات گزارش شده (2، 5، 10 و15) ، چندین جنس و گونه جدید که تاکنون وجود توالیهای مشابه و بنابراین خاصیت القای مرگ سلولی و احتمالاً حشره‌کشی در آنها ردیابی نشده بود را آشکار نمود. از این رو، این یافته می‌تواند در جهت دسترسی و استفاده از منابع و ذخایر وسیع‌تر پروتئینی برای سنتز حشره‌کشهای میکروبی جدید و مقابله با بروز مقاومت در آینده راه‌گشا باشد (17). در نتایج حاصل از تحلیل ساختار دوم آشکار‌ شد که در زیرواحد‌های XptA1 و XptA2 دمین عملکردی با ویژگی احتمالی حشره‌کشی بیشتر از هلیکسها تشکیل شده است ( جدول 3). در همین راستا، پایداری نسبت به دما و مواد شیمیایی با افزایش میزان این نوع ساختارها در پروتئین نمایش داده شده است (14)، لذا دمینهای آشکارشده در این پروتئینها در صورت استفاده مجزا در ساختارهای هیبریدی ضمن آنکه می‌توانند دارای اثر مورد نظر در کنترل زیستی ‌باشند، دارای قابلیت ماندگاری در شرایط محیطی خواهند بود. در زیر واحد XptC1 نتایج عکس مشاهده شد و دمین مؤثر بر اکتین و دمین احتمالی اتصالی بیشتر دارای ساختار صفحات بوده که از لحاظ مقاومت و پایداری در مقایسه با موارد بالا می‌تواند متفاوت باشد. بررسی وجود ساختار‌های Coiled-coil، تکرار‌های 7 تایی با اولین و چهارمین اسید‌آمینه که اغلب هیدروفوب هستند و پنجمین و هفتمین که غالباً باردار یا قطبی هستند (19)، در زیرواحد  XptA1 آشکارکننده یک توالی از این ساختار، مشابه با نتیجه گزارش شده قبلی در این حوزه است (19). در همین راستا، در خصوص زیرواحد XptA2 سه ناحیه با این ساختار پیش‌بینی شد. این نواحی احتمالاً برهمکنش Coiled-coil را با زیرواحدهای مجاور شکل می‌دهند، با این وجود تا کنون در این مورد گزارشی ارائه نشده است. لازم به ذکر است که گفته شود این ساختار گزینه ای مناسب برای ایجاد بر‌همکنش بین زیرواحد‌ها با یکدیگر هستند (23)، از این رو می‌توان نتیجه گرفت که زیرواحد‌ها به صورت مجموعه  با یکدیگر همکاری می‌کنند مشابه آنچه تاکنون انتظار میرفت (26). ارزیابی ساختار‌های سه بعدی طراحی شده در ProSA به ویژه در زمینه نمودار‌های انرژی کل (شکل 2) نشان داد که مقادیر امتیاز همه موارد منفی و مناسب بوده ولی به دلیل قرار نگرفتن در محدوده مجاز نمودار احتمالاً کیفیت ساختار ‌در زیرواحد‌های XptA1 ، XptA2 و XptB1 در حد عالی نیستند که می‌تواند به دلیل محدودیت تعداد و همسانی در توالیهای الگو یافت شده و همچنین محدودیتهای نرم‌افزاری و ارزیابی باشد، به دلیل اینکه تاکنون دمین دارای عملکرد حشره‌کشی در XptB1 آشکار نشده است. بنابراین از مدل فضایی آن نیز می‌توان صرف‌نظر کرد. ارزیابی دو زیرواحد اول در خصوص قرار گرفتن خارج از محدوده مجاز نمودار نتایج مطلوبی نداشت با اینکه معیار‌های مدل‌سازی این‌دو مانند dope score و GA341 score در حد بسیار عالی حاصل شد که احتمالاً به دلیل محدودیت پایگاه proSA در سنجش پروتئینهایی با اندازه بزرگتر از 1000 اسیدآمینه می‌باشد. با این حال نمودار‌های انرژی محلی به دست آمده و ساختار سه بعدی نمایش داده شده (شکل 4) نشان دهنده پایین بودن سطح انرژی و بنابراین پایدار بودن ساختار است. بررسی نمودار راماچاندران حاکی از قرارگرفتن تعداد نسبتاً زیادی از اسید‌های آمینه در منطقه غیرمجاز است که البته تعداد کمی از این تعداد در موقعیت دمین عملکردی قرار دارند (جدول 4)، که البته باید در پژوهشهای تکمیلی بررسی شوند که آیا این اسیدهای آمینه تاثیری در عملکرد دمین مربوطه دارند و همچنین شاید بتوان با اتخاذ روشهای دیگر و پیشرفته‌تر مدل‌سازی، بهینه سازی و ارزیابی و به خصوص روش تجربی و آزمایشگاهی به نتیجه قطعی و مطلوب رسید. نتایج مشاهده شده در ساختار سه بعدی زیرواحد XptA1 ساختار قفس مانند و حفره میانی را مشابه با ساختار XptA1 بررسی شده در مقاله Lee و همکاران نشان می‌دهد که می‌تواند تا حدودی اثبات کننده صحت مدل‌سازی مجازی باشد (19). به علاوه شباهت توالی مشاهده شده بینXptA1 و XptA2 که قبلاً عنوان شدند (21) منجر به ایجاد تشابه در ساختار سه بعدی این دو نیز شد (شکل 3 و 5). مطالعات و تحلیلهای میل اتصالی با اکتین آشکار کننده وجود جایگاههای اتصالی در هر چهار زیرواحد شد که نشان می‌دهد تمامی این اجزاء دارای فعالیت و اثر بر اسکلت سلولی اکتین هستند. با این وجود میل اتصالی XptA1 نسبت به سایر اجزاء پایین‌تر حاصل شده که می تواند اثبات کننده فعالیت حشره‌کشی پایین این زیرواحد به تنهایی و نیاز آن به زیرواحد‌های دیگر برای القای سمیت کامل در بدن حشره داشته باشد (21). همچنین بیشترین میل اتصال به پروتئین اکتین در زیر واحد XptC1 حاصل شد که با نتایج حاصله از جدول 3 در ارتباط با وجود دمین عملکردی با وظیفه اثر بر اکتین کاملاً موافق می‌باشد.

1- جواهری مقدم، م.، آریاپور، ح.، دهنوخلجی، ا، 1394. طراحی ترکیبات کرومنی جدید با فعالیت ضد سرطانی و بررسی چگونگی میانکنش آنها با توبولین به روش داکینگ مولکولی. مجله زیست شناسی ایران، ج 28، ش 2. صفحه 190-178

2- شکری، ا.، نصیری، ن.، نعمت زاده، ق، 1393. همسانه سازی مولکولی ترادف رمزکننده ژن H+- ATPase غشای پلاسمایی در گراس شورزی Aeluropus littoralis. مجله زیست شناسی ایران، ج 27، ش 3. صفحه 388-377

3- Ansari, m.s., Ahmad, s., Ahmad, n., Ahmad, t. & Hasan, f. (2013).Microbial insecticides: food security and human health. In Management of Microbial Resources in the Environment, 341-360: Springer.

4- Bowen, d. (2000). Novel insecticidal toxins from nematode-symbiotic bacteria. Cellular and Molecular Life Sciences CMLS 57(5): 828-833.

5- Bowen, d. j. & Ensign, j. c. (1995). Characterization of a high molecular weight insecticidal protein complex produced by the entomopathogenic bacterium Photorhabdus luminescens.

6- Brown, s. e., Cao, a. t., Hines, e. r., Akhurst, r. j. & East, p. d. (2004). A novel secreted protein toxin from the insect pathogenic bacterium Xenorhabdus nematophila. Journal of Biological Chemistry 279(15): 14595-14601.

7- Castagnola, a. & Stock, s. p. (2014). Common virulence factors and tissue targets of entomopathogenic bacteria for biological control of lepidopteran pests. Insects 5(1): 139-166.

8- Chattopadhyay, a., Bhatnagar, n. & Bhatnagar, r. (2004). Bacterial insecticidal toxins. Critical reviews in microbiology 30(1): 33-54.

9- Delepelaire, p. & Wandersman, c. (1990). Protein secretion in gram-negative bacteria. The extracellular metalloproteaseB from Erwinia chrysanthemi contains a C-terminal secretion signal analogous to that of Escherichia coli alpha-hemolysin. Journal of Biological Chemistry 265(28): 17118-17125.

10- Dowling, a. & Waterfield, n.r. (2007). Insecticidal toxins from Photorhabdus bacteria and their potential use in agriculture. Toxicon 49(4): 436-451.

11- Duhovny, d., Nussinov, r. & Wolfson, h.j. (2002).Efficient unbound docking of rigid molecules. In International Workshop on Algorithms in Bioinformatics, 185-200: Springer.

12- Federici, b. a. (2005). Insecticidal bacteria: an overwhelming success for invertebrate pathology. Journal of invertebrate pathology 89(1): 30-38.

13- Forst, s., Dowds, b., Boemare, n. & Stackebrandt, e. (1997). Xenorhabdus and Photorhabdus spp.: bugs that kill bugs. Annual Reviews in Microbiology 51(1): 47-72.

14- Guo, j., Harn, n., Robbins, a., Dougherty, r. & Middaugh, c.r. (2006). Stability of helix-rich proteins at high concentrations. Biochemistry 45(28): 8686-8696.

15- Hurst, m.r., Glare, t.r., Jackson, t.a. & Ronson, c.w. (2000). Plasmid-located pathogenicity determinants of Serratia entomophila, the causal agent of amber disease of grass grub, show similarity to the insecticidal toxins of Photorhabdus luminescens. Journal of Bacteriology 182(18): 5127-5138.

16- Koronakis, v. & Hughes, c. (1993).Bacterial signal peptide-independent protein export: HlyB-directed secretion of hemolysin. In Seminars in cell biology, Vol. 4, 7-15: Elsevier.

17- Kumari, p., Kant, s., Zaman, s., Mahapatro, g.k., Banerjee, n. & Sarin, n. b. (2014). A novel insecticidal GroEL protein from Xenorhabdus nematophila confers insect resistance in tobacco. Transgenic research 23(1): 99-107.

18- Lang, a. e., Schmidt, g., Schlosser, a., Hey, t. d., Larrinua, i. m., Sheets, j.j., Mannherz, h. g. & Aktories, k. (2010). Photorhabdus luminescens toxins ADP-ribosylateactin and RhoA to force actin clustering. Science 327(5969): 1139-1142.

19- Lee, s. c., Stoilova-Mcphie, s., Baxter, l., Fülöp, v., Henderson, j., Rodger, a., Roper, d. i., Scott, d. j., Smith, c. j. & Morgan, j. a. w. (2007). Structural characterization of the insecticidal toxin XptA1, reveals a 1.15 MDa tetramer with a cage-like structure. Journal of molecular biology 366(5): 1558-1568.

20- Montesinos, e. (2003). Development, registration and commercialization of microbial pesticides for plant protection. International Microbiology 6(4): 245-252.

21- Morgan, j. a. w., Sergeant, m., Ellis, d., Ousley, m. & Jarrett, p. (2001). Sequence Analysis of Insecticidal Genes from Xenorhabdus nematophilus PMFI296. Applied and environmental microbiology 67(5): 2062-2069.

22- Pérez-García, a., Romero, d. & De Vicente, a. (2011). Plant protection and growth stimulation by microorganisms: biotechnological applications of Bacilli in agriculture. Current Opinion in Biotechnology 22(2): 187-193.

23- Poinar Jr, g. o. (1986). Entomophagous nematodes. Fortschritte der Zoologie (Germany, FR.)

24- Ruiu, l., Satta, a. & Floris, i. (2013). Emerging entomopathogenic bacteria for insect pest management. Bull. Insectology 66: 181-186.

25- Sainsbury, f., Benchabane, m., Goulet, m.-c. & Michaud, d. (2012). Multimodal protein constructs for herbivore insect control. Toxins 4(6): 455-475.

26- Sergeant, m., Baxter, l., Jarrett, p., Shaw, e., Ousley, m., Winstanley, c. & Morgan, j. a. w. (2006). Identification, typing, and insecticidal activity of Xenorhabdus isolates from entomopathogenic nematodes in United Kingdom soil and characterization of the xpt toxin loci. Applied and environmental microbiology 72(9): 5895-5907.

27- Sergeant, m., Jarrett, p., Ousley, m. & Morgan, j. a. w. (2003). Interactions of insecticidal toxin gene products from Xenorhabdus nematophilus PMFI296. Applied and environmental microbiology 69(6): 3344-3349.

28- Sudan, a. k. & Vakhlu, j. (2015). Isolation and in silico characterization of novel esterase gene with β-lactamase fold isolated from metagenome of north western Himalayas. 3 Biotech 5(4): 553-559.

29- Wang, q. y., Nangong, z. y., Yang, j., Song, p., Wang, y., Cui, l. & Cui, l. (2012). Toxic activity of a protein complex purified from Xenorhabdus nematophila HB310 to Plutella xylostella larvae. Insect Science 19(3): 329-336.

30- Wiederstein, m. & Sippl, m. j. (2007). ProSA-web: interactive web service for the recognition of errors in three-dimensional structures of proteins. Nucleic acids research 35(suppl 2): W407-W410.