تعداد نشریات | 43 |
تعداد شمارهها | 1,652 |
تعداد مقالات | 13,408 |
تعداد مشاهده مقاله | 30,253,287 |
تعداد دریافت فایل اصل مقاله | 12,089,822 |
بررسی ساختار و عملکرد پروتئین اسپایک کروناویروس سندرم حاد تنفسی ۲ (SARS-CoV-2) با داکینگ مولکولی | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
زیست شناسی میکروبی | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
دوره 12، شماره 46، تیر 1402، صفحه 27-50 اصل مقاله (1.64 M) | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
نوع مقاله: پژوهشی- فارسی | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
شناسه دیجیتال (DOI): 10.22108/bjm.2022.133412.1468 | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
نویسندگان | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
نوید محمدجانی1؛ مراحم آشنگرف* 1؛ فرشاد درویشی* 2 | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1گروه علوم زیستی، دانشکده علوم، دانشگاه کردستان، سنندج، ایران | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
2گروه میکروبیولوژی، دانشکده علوم زیستی، دانشگاه الزهرا، تهران، ایران | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
چکیده | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
مقدمه: بیش از دو سال از شیوع سندرم حاد تنفسی ویروس کرونا 2 (SARS-CoV-2) در ووهان چین میگذرد. SARS-CoV-2 یک پاندمی بیسابقه (COVID-19) در عصر مدرن ایجاد کرده و تمام ساختارهای علمی، اقتصادی و اجتماعی جهان را به چالش کشیده است. تحقیقات گستردهای در سراسر جهان برای مطالعه ساختار و عملکرد پروتئین اسپایک SARS-CoV-2 از ابتدای پاندمی انجام شدهاند. بنا بر مطالعات مختلف، نقش و اهمیت پروتئین اسپایک در ورود SARS-CoV-2 به سلولهای انسانی ثابت شده است. ازجمله عوامل مهم و مؤثر در کنترل همهگیری کووید-19، درک تغییرات ساختاری پروتئین اسپایک در واریانتهای ویروس و اثربخشی داروها و واکسنها در برابر واریانتهای جدید SARS-CoV-2 است. با توجه به اهمیت فراوان پروتئین اسپایک SARS-CoV-2، سعی شد در این مطالعه ساختار و تعامل آن با گیرندههای سلولی با استفاده از ابزار بیوانفورماتیک بررسی شود. علاوه بر این، اثر داروها و واکسنهای موجود بر وارینتهای مختلف SARS-CoV-2 بررسی شد. مواد و روشها: روشهای بیوانفورماتیک و داکینگ مولکولی برای مطالعه عملکرد اتصال پروتئین اسپایک به گیرندههای سلولی مختلف آن در بدن انسان استفاده شدند. برای انجام فرایند داکینگ مولکولی از وب سرور HDOCK و در مرحله بعد از وب سرور PDBsum برای بررسیهای پس از داکینگ و اطمینان از صحت دادههای بهدستآمده استفاده شد. همچنین از نرمافزار PyMOL برای مطالعه و تجسم جهشهای جایگزین و حذف در پروتئین اسپایک استفاده شد. نتایج: تفاوت در پروتئین اسپایک واریانتهای SARS-CoV-2 با تصویرسازیهای انجامگرفته، توسط نرمافزار PyMOL مقایسه شد. دادههای حاصل از جهشهای گسترده در پروتئین اسپایک واریانت اومیکرون نشاندهندة نوعی بازطراحی دوبارة پروتئین اسپایک بود. این جهشهای وسیع، تغییرات گستردهای در ساختار و عملکرد پروتئین اسپایک به همراه دارد. براساس دادههای حاصل از داکینگ مولکولی و مقایسه میزان سطح انرژی (امتیاز داکینگ) اتصال پروتئین اسپایک واریانتهای مختلف ویروس با گیرندههای مختلف نشاندهندة بالابودن تمایل اتصال - با سطح انرژی پایینتر و پایدارتر پروتئین اسپایک - با گیرنده KREMEN1 نسبت به گیرنده اصلی ACE2 است. میزان سطوح انرژی اتصال پروتئین اسپایک در واریانتهای مختلف ویروس با سایر گیرندهها غیر از گیرنده KREMEN1 تقریباً مشابه بود یا اختلاف اندکی با سطوح انرژی اتصال پروتئین اسپایک در واریانتهای مختلف ویروس با گیرنده اصلی ACE2 دارد. بحث و نتیجه گیری: تفاوت زیادی در تمایل اتصال پروتئین اسپایک واریانتهای نوظهور و مهم ویروس (اصلی، آلفا، دلتا و اومیکرون) با گیرنده ACE2 وجود ندارد؛ اما میزان تمایل اتصال به گیرنده KREMEN1، رو به افزایش و در واریانت دلتا به اوج خود رسیده است. سطح انرژی اتصال پروتئین اسپایک به گیرنده KREMEN1 بهطور محسوسی نسبت به سایر گیرندهها منفیتر بوده و این نشاندهندة ایجاد اتصالات و پایداری بیشتر است. با توجه به سطح انرژیهای دریافتی از وب سرور HDOCK، میتوان نتیجه گرفت اگرچه گیرنده ACE2 گیرنده اصلی پروتئین اسپایک نام گرفته است، همچنان گیرندههای دیگری نیز با پایداری خوبی میتوانند به پروتئین اسپایک متصل شوند. براساس اطلاعات موجود، تاکنون هیچ گزارشی مبنی بر انجام فرایند داکینگ بین پروتئین اسپایک با سایر گیرندههای غیراصلی پروتئین اسپایک (AXL، ASGR1 و KREMEN1) صورت نگرفته است؛ بنابراین، برای تأیید یا رد این یافتهها باید پژوهشهای جامعتری انجام شود. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
کلیدواژهها | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
SARS-CoV-2؛ اسپایک؛ گیرنده ACE2؛ جهش؛ واریانتهای جدید؛ بیوانفورماتیک | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
اصل مقاله | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
مقدمه. بیماری کروناویروس ۲۰۱۹[1] یا کووید-19[2] یک بیماری عفونی نوظهور است که در اواخر سال 2019 آغاز شد و بهصورت پاندمی جهان را درگیر کرد. این بیماری توسط یک ویروس جدید از خانواده کروناویروسها به نام کروناویروس سندرم حاد تنفسی ۲ (SARS-CoV-2) به وجود آمد (1). SARS-CoV-2 ازطریق پروتئین اسپایک موجود در سطح خود به سلولهای میزبان حمله میکند. گیرنده اصلی این پروتئین در سلولهای انسانی، آنزیم مبدل آنژیوتانسین 2 (ACE2)[3] است و یک عامل تعیینکننده برای انتقال بین گونهای است. پس از اتصال اسپایک ویروس به گیرنده، پروتئازهای خاصی وارد عمل میشوند و این پروتئین را به زیرواحدهای آن یعنی S1 و S2 تقسیم میکنند (2). همانطور که اشاره شد، آنزیم مبدل آنژیوتانسین 2 (ACE2)، گیرنده اصلی برای ورود به سلول است؛ اما سلولها در اندام هدف اصلی (ریهها)، سطوح پایینی از بیان ACE2 دارند؛ درحالیکه SARS-CoV-2 بهطور مؤثری میتواند ریهها را درگیر کند. این شواهد نشان میدهند گیرندههای جایگزین دیگری میتوانند در فرایند ورود ویروس بهعنوان گیرنده نقش داشته باشند. بنا بر مشاهدات آزمایشگاهی گیرنده پروتئینی تایروزین کینازی (AXL)[4]، نشانگر پروتئینی در چشم و بینی (KREMEN1)[5] و گیرنده ASGR1[6] نیز میتوانند در ورود ویروس به سلولهای میزبان نقش داشته باشد (1و2). تاکنون واریانتهای مختلفی از این ویروس شناسایی شدهاند که در همه آنها پروتئین اسپایک مهمترین نقش را در جهشهای مربوطه داشته است؛ زیرا این پروتئین در فرایند ورود ویروس به سلول میزبان نقش اصلی را دارد و به همین دلیل هرگونه جهش در آن، تغییرات گستردهای در عملکرد پروتئین ایجاد میکند که بر میزان سرایتپذیری، بار ویروسی و همچنین افزایش یا کاهش مرگومیر نقش دارد (3). SARS-CoV-2 متعلق به خانواده بتاکروناویروسها است که بزرگترین ژنوم RNA را دارد. این ژنوم، 29 پروتئین ساختاری، غیرساختاری و کمکی را کد میکند که در ورود ویروس به سلولهای میزبان، تکثیر ژنوم، رونویسی، سرهمبندی و انتشار ویروس شرکت میکنند (4). بنا بر بررسیهای انجامگرفته بخشی از ژنوم SARS-CoV-2 و [7]MERS-CoV (سویه خطرناک قبلی از این خانواده ویروسی که در سال 2012 باعث بیماری شبیه سارس[8] به نام مرس[9] یا سندرم تنفسی کشنده خاورمیانه[10] شد) یکسان است (5و6). با توجه به جهشهایی که در SARS-CoV-2 رخ میدهد، یکی از نکات کلیدی درک عمیق تمامی این جهشها و شناسایی زودهنگام جهشهای مهم برای معرفی آنها بهعنوان واریانتهای خطرناک است تا جامعه بشری قبل از اینکه یک جهش بتواند به شکل گستردهای منتشر شود، با اتخاذ تدابیر پیشگیرانه از این امر جلوگیری کند. همچنین درک علت تکامل ویروس کرونا به مدل فعلی یعنی SARS-CoV-2 هنوز بهطور دقیق شناسایی نشده و بدون شک انجام پژوهشهای ترکیبی بیوانفورماتیکی و آزمایشگاهی برای شناسایی علل این امر لازم است تا از تکرار پاندمیهای مشابه در آینده پیشگیری کند. همچنین ارزیابی اثربخشی تمامی واکسنهای رایج در جهان تاکنون تا حد خوبی انجام شده است؛ اما به تحقیقات بیشتری در این زمینه برای کشف دقیقتر میزان اثرگذاری واریانتهای جدید بر اثر بخشی واکسنها لازم است (3و7). هدف این پژوهش بررسی بیوانفورماتیکی ساختار و عملکرد پروتئین اسپایک SARS-CoV-2 برای شناسایی و پیشبینی عملکرد واریانتهای جدید بر اثربخشی واکسنها و داروها است. بدین منظور در این مطالعه واریانتهای مهم پروتئین اسپایک با استفاده از تکنیک داکینگ مولکولی با گیرندههای بالقوه آن بررسی شدهاند تا نحوه اتصال و میزان تمایل هریک برای اتصال به یکدیگر مشخص شوند.
مواد و روشها. .بررسی بیوانفورماتیکی ساختار پروتئین اسپایک SARS-CoV-2: در این مطالعه برای بررسی بیوانفورماتیکی ساختار پروتئین اسپایک SARS-CoV-2، ابتدا دادههای پروتئین اسپایک با فرمت pdb با IDهای 6VXX، 7XK3، 7FEM، 7BNM، 7LWU، 7VX1، 7SBK، 6VYB و 6VXX از پایگاه دادة پروتئینی PDB[11] دریافت شدند. در مرحله بعد هریک از این دادهها با نرمافزار PyMOL به شکل جداگانه از لحاظ توالی و ساختار سهبعدی بررسی شدند. .پیشپردازش[12] دادههای پروتئینی اسپایک واریانتهای مختلف ویروس و گیرندههای مختلف سلولی در بدن انسان: در ادامه باید نحوه اتصالات پروتئین اسپایک در واریانتهای با اهمیت SARS-CoV-2 بررسی شوند. برای این کار دادهها با فرمت pdb پروتئین اسپایک در ویروس پایه (ووهان)، واریانتهای آلفا، دلتا و اومیکرون و همچنین دادههای مربوط به گیرندههای ACE2، ASGR1، KREMEN1 و AXL از پایگاه دادة پروتئینی PDB دریافت شدند. مطرحبودن احتمال دخالت گیرندههای دیگر موجود در بدن انسان، غیر از گیرنده اصلی ACE2 در مطالعات مختلف، دلیل دریافت و بررسی دادههای پروتئینی گیرندههایی غیر از گیرنده اصلی یعنی ACE2 بود. قبل از بررسی نحوه اتصالات پروتئین اسپایک با گیرندههای سلولی و داکینگ مولکولی[13]، دادههای پروتئینی باید پیشپردازش شوند. برای انجام این فرایند، ابتدا دادههای HETATOM دادههای دریافتی از پایگاه PDB(با IDهای 7DK3،7VX4 ، 6YAU، 5U6B، 6SNW، 7EDI، 7W92 و 7TGW) با نرمافزار ++Notepad (نسخه 8.3.3) حذف شدند. سپس عمل پیشپردازش دادهها با نرمافزار UCSF Chimera (نسخه 1.16) ادامه یافت؛ بدین صورت که زنجیره یا ماده اضافی همراه با پروتئین اسپایک یا گیرنده سلولی گزارششده، همگی حذف شدند و پروتئین اسپایک یا گیرنده خالص به دست آمد. در گام بعدی، آمینواسیدهای غیراستاندارد حذف شدند و سپس با استفاده از گزینه Dock Prep نرمافزار تمامی دادهها بنا به پیشفرض نرمافزار، پیشپردازش و آماده انجام فرایند داکینگ مولکولی شدند. داکینگ مولکولی دادههای پروتئینی اسپایک واریانتهای مختلف ویروس با گیرندههای مختلف سلولی در بدن انسان: در این مرحله هدف بررسی تمایل اتصال کلی پروتئین اسپایک و گیرندههای سلولی بود؛ به همین دلیل از روش داکینگ سرتاسری[14] استفاده شد تا تمایلات اتصال تمامی بخشهای پروتئین اسپایک در واریانتهای مهم و گیرندههای سلولی با هم بررسی شوند. بدین منظور از وب سرور [15]HDOCK استفاده شد و 16 حالت مختلف (4 واریانت پروتئین اسپایک و 4 نوع گیرنده) بررسی شدند. بهدلیل اینکه ساختار باز پروتئین اسپایک در اتصال به گیرنده نقش دارد، تمامی فایلهای مربوط به پروتئین اسپایک، فرم باز این پروتئین را داشتند. بعد از انجام محاسبات توسط سرور و اطلاعات بهدستآمده، 10 مدل اصلی با اولویت مدل اول و پایدارتر تجزیه و تحلیل شدند. پردازش نهایی پس از داکینگ مولکولی[16]: در این مرحله، مدل یک و بهترین مدل مطرحشده توسط وب سرور HDOCK (از لحاظ سطح انرژی و RMSD[17]) در رابطه با 16 حالت، بررسی و برای انجام پردازش پس از داکینگ و تحلیل دادههای وب سرور HDOCK استفاده شد. RMSD، درجه انطباق دو جزء با هم Dock شده است و این مقدار عددی هرچه به صفر نزدیکتر باشد، یعنی دو جزء (برای مثال پروتئین اسپایک و گیرنده) بهتر روی هم منطبق شدهاند. برای انجام فرایند پس از داکینگ مولکولی از وب سرور PDBsum[18] استفاده شد. این مرحله درواقع برای اطمینان از درستی دادهها و فرایند داکینگ انجام شد. بحث و نتیجهگیری .ساختار و عملکرد پروتئین اسپایک SARS-CoV-2 در مطالعات بیوانفورماتیکی: با توجه به اهمیت پروتئین اسپایک، در مکانیسمهای ورود ویروس به سلول میزبان، تلاشهای زیادی برای مطالعه ساختاری این پروتئین انجام گرفته است؛ این مطالعات ساختاری را میتوان به دو دسته تقسیم کرد: 1- روشهای آزمایشگاهی و 2- روشهای بیوانفورماتیکی (4). میکروسکوپ الکترونی کرایو (Cryo-EM): یکی از مناسبترین و در عین حال بهترین روش برای شناسایی ساختار پروتئین اسپایک استفاده از میکروسکوپ الکترونی کرایو است که میتواند تا سطح اتمی، پروتئین اسپایک را شناسایی کند و در ارزیابی و پیشبینی نحوه عملکرد این پروتئین بسیار خوب عمل کرده است. در میکروگراف الکترونی پروتئین اسپایک به راحتی میتوان تریمرهای آزاد و کمپلکسهای تریمر را دید. با استفاده از میکروگراف، حالت بسته و باز دومین اتصال گیرنده (RBD)[19] تشخیص داده میشود که درنتیجه میتوان با بررسی پروتئین اسپایک در حالت پیش از ادغام و پس از ادغام به تغییر شکل فضایی طی ادغام غشایی پی برد. همچنین با استفاده از این میکروگراف، زیرواحدهای S1 و S2 پروتئین اسپایک به شکلی کاملاً واضح مشاهده میشوند (8و50). .آنالیز گلیکانی سایت ویژه اسپایک (Site-specific glycan analysis): گلیکانها در گلیکوپروتئینهای ویروسی نقش بیبدیلی در تاخوردگی پروتئین، پایداری، شناسایی توسط سیستم ایمنی و فرار از سیستم ایمنی دارند. آنالیزهای گلیکانی که برای ارزیابی میزان گلیکوزیلاسیون انجام شدند، گلیکوزیلاسیون را در تمام 22 توالی گلیکان متصل به N موجود در پروتئین اسپایک SARS-CoV-2 نشان دادند. براساس مطالعات، 98 درصد پروتئین اسپایک SARS-CoV-2 گلیکوزیله است؛ اما در موقعیتهای 603 و 657 پروتئین اسپایک گلیکانی وجود ندارد؛ درحالیکه سایر نقاط پروتئین اسپایک به شدت گلیکوزیلهاند. نکته جالب توجه این است که گلیکوزیلهبودن اسپایک در بعضی از نقاط دیگر آن با میکروسکوپ الکترونی کرایو مشاهده نشده است که درواقع نشان میدهد آنالیز گلیکانی میتواند تکمیلکنندة روش میکروسکوپ الکترونی باشد (8). همترازی توالی: ترازکردن و همارزکردن توالیهای پروتئین اسپایک SARS-CoV-2 و SARS-CoV برای مشخصکردن میزان مشابهت و تفاوتهای این دو در سطح توالی انجام شدهاند که اینگونه مطالعات میتوانند بینش بسیار دقیقی را برای تولید واکسنها و مبارزه با کووید-19 در اختیار ما قرار بدهند (3). تجزیه و تحلیل جهش: اطلاعات توالی در این مطالعات از پایگاه داده GISAID[20] (ابتکار جهانی برای به اشتراکگذاری تمامی دادههای آنفلوانزا پرندگان) به دست میآید. در ابتدا توالیهای نوکلئوتیدی گلیکوپروتئین اسپایک با استفاده از نرمافزارهای بیوانفورماتیکی تبدیل به توالیهای پروتئینی میشوند. سپس این توالیهای پروتئینی با واریانت ووهان، همردیف و نقاط جهشیافته تجزیه و تحلیل میشوند (26). مطالعات داکینگ مولکولی: این رویکرد برای پیشبینی اتصال پروتئین اسپایک با گیرنده سلولی انجام میشود. همچنین این دادهها میتوانند در پیشبینی و ارزیابی اولیه توانایی مواد مختلف بهعنوان دارو در ایجاد تداخل در اتصال پروتئین اسپایک به گیرنده استفاده شوند (51). آنالیز مبتنی بر ساختارهای اپیتوپهای ساختاری روی پروتئین اسپایک: در این روش به آنتیبادیهای گوناگون که به نواحی مختلف پروتئین اسپایک حمله میکنند، امتیازات مختلفی براساس آنالیزهای بیوانفورماتیکی داده میشود. در این رویکرد میزان دسترسی آنتیبادی برای پروتئین اسپایک در حالت باز و بسته، تعیین و به شکل نمودارهای heatmap تحلیلی به نمایش گذاشته میشود (43). مطالعه حاضر براساس روشهای بیوانفورماتیکی بهویژه مطالعات داکینگ مولکولی انجام شده است. همانطور که در شکل 1 مشخص است، توالی 1273 آمینواسیدی پروتئین اسپایک SARS-CoV-2 دارای زیرواحد S1 (آمینواسیدهای 14 تا 685) و زیرواحد S2 (آمینواسیدهای 686 تا 1273) است که با رنگهای متفاوت از هم تفکیک شدند (1و2). گفتنی است جایگاه برخی از ریشههای آمینواسیدی در داخل زیرواحدها از دادههای دریافتی از پایگاه دادة PDB حذف شده بود که این امر بهعلت 100 درصد نبودن بازده دستگاههای توالییاب است و به همین دلیل در این مطالعه، این نواحی نادیده گرفته شدند. علاوه بر این در شکل 1، ساختار باز و بسته پروتئین اسپایک و حالت متصل به گیرنده ACE2 نشان داده شده است. پروتئین اسپایک یک گلیکوپروتئین است که بخشهای پروتئینی توسط مولکولهای پلیساکاریدی پوشیده میشوند تا آنها را استتار کنند و به فرار ویروس از سیستم ایمنی میزبان منجر شوند. اسپایک نقش کلیدی در شناسایی گیرنده و فرایند ادغام غشایی SARS-CoV-2 بر عهده دارد. در این پروتئین سه پیشساز به هم میرسند تا یک تریمر را تشکیل دهند و از دو زیرواحد S1 و S2 تشکیل شده است (1و3). پروتئین اسپایک حاوی 22 سایت به شدت گلیکوزیله است، مطالعات نشان میدهند گلیکوزیلاسیون پروتئین اسپایک میتواند بر عفونتزایی SARS-CoV-2 مؤثر باشد و کاهش گلیکوزیلاسیون این پروتئین میتواند میزان عفونتزایی ویروس را کاهش دهد (7). بنا بر بررسیهای صورتگرفته با میکروسکوپ الکترونی کرایو زیرواحد S1 حول یک محور سهگانه میپیچید و زیرواحد S2 را میپوشاند (4). عناصر ساختاری پروتئین اسپایک شامل یک توالی سیگنال جداشده (SS)[21]، دومین پایانه N (NTD)[22]، دومین اتصال به گیرنده (RBD)، محل برش پروتئاز، فیوژن پپتید (FP)[23]، هپتاد تکرارشونده 1(HR1)[24]، هپتاد تکرارشونده 2 (HR2)[25] هستند (8). بهطور متوسط 30 تا 60 تریمر اسپایک از پوشش SARS-CoV-2 با فاصله حدود 15 نانومتری از یکدیگر بیرون آمدهاند. زیرواحد S1 حاوی دومین اتصال به گیرنده (RBD) است که مسئولیت اتصال به گیرنده آنزیم مبدل آنژیوتانسین 2 (ACE2) را بهعنوان گیرنده اصلی ویروس دارد و زیرواحد S2 میانجی همجوشی ویروس با غشای سلولی میزبان است. اسپایک با اندازه 180 تا 200 کیلودالتون و با طول 1273 آمینواسید، یک پایانه N خارج سلولی دارد که حاوی زیرواحد S1 (آمینواسیدهای 14 تا 685) مسئول اتصال به گیرنده و زیرواحد S2 (آمینواسیدهای 686 تا 1273) مسئول همجوشی با غشا است. در زیرواحد S1 یک دومین N ترمینال (آمینواسید 14 تا 305) و دومین اتصال به گیرنده یا RBD (آمینواسید 319 تا 541) وجود دارد (قسمت A شکل 1) (1و2). زیرواحدهای S1 و S2 اسپایک در فضای خارج SARS-CoV-2 قرار میگیرند؛ اما بخش گذرندة غشایی[26] آن درون غشا قرار دارد و یک دومین پایانه C نیز در فضای داخل ویروس قرار میگیرد (9). اسپایک به شکل عادی در حالت پیشادغام[27] قرار دارد؛ اما هنگامی که ویروس با سلول میزبان تعامل پیدا میکند، تغییر شکل فضایی ساختاری گستردهای در پروتئین اسپایک اتفاق میافتد که درنهایت به ترکیب ویروس با غشای سلولی میزبان منجر خواهد شد. دومین اتصال به گیرنده قبل از اتصال در حالت بسته قرار دارد؛ اما در هنگام اتصال شکل باز آن به ACE2 متصل میشود. اتصال به ACE2 باعث ایجاد تغییرات ساختاری میشود که تعاملات بین زیرواحدهای S1 و S2 را بیثبات میکند. هنگام عفونت ویروسی زیرواحدهای S1 و S2 اسپایک SARS-CoV-2 توسط پروتئازهای سلولی از هم جدا میشوند و بیشتر بهصورت تریمرهای هترودایمر S1/S2 تجمع مییابند و به این روش فعالسازی پروتئین صورت میگیرد؛ سرین پروتئازTMPRSS2 [28] فعالکنندة اصلی اسپایک است. دومین اصلی برای فرایند پروتئولیز، در سایت S2 قرار دارد که توسط TMPRSS2 یا سایر پروتئازها یک پایانه N آزاد برای فیوژن پپتید ایجاد میکند و باعث ایجاد تغییرات ساختاری برای در دسترس قرارگرفتن فیوژن پپتید برای غشای سلولی میشود. این قسمت درواقع محل شکست اصلی پروتئین اسپایک است. فیوژن پپتید یک سرین پروتئاز اولیه است که در بسیاری از سلولهای اپیتلیال بیان میشود؛ به همین دلیل یک هدف ضدویروسی ایدئال برای مهار ورود ویروس است. علاوه بر این پروتئاز اصلی، سرین اندوپروتئاز پروپروتئین کانورتاز 1[29]، تریپیسن[30]، ماتریپتاز (سرین پپتیداز غشای شبهتریپسین[31]) و کاتپسین[32] در فرایند پروتئولیز اسپایک طی ورود ویروس به سلول میزبان میتوانند نقش داشته باشند. این تنوع در گیرندهها نیز وجود دارد (1و2).
شکل 1 - ساختار پروتئین اسپایک SARS-CoV-2. A) ساختار خطی پروتئین اسپایک SARS-CoV-2 و محل برش فورین. B) پروتئین اسپایک در حالت بسته؛ در این حالت پروتئین نمیتواند به گیرنده ACE2 متصل شود، زیرواحد S1 (قرمز)، زیرواحد S2 (آبی) (PDB ID 6VXX). C) پروتئین اسپایک در حالت باز و آماده اتصال (PDB ID 7XK3). D) پروتئین اسپایک درحال اتصال به گیرنده ACE2 که گیرنده ACE2به رنگ سفید نشان داده شده است (PDB ID7FEM).
همانطور که اشاره شد هنگامی که زیر واحد S1 به ACE2 متصل میشود یک ادغام رخ میدهد؛ مارپیچ N ترمینال ACE2 با سطح بیرونی دومین متصلشونده به گیرنده، در زیرواحد S1 متصل میشود که این اتصال شامل 16 آمینواسید در دومین متصلشونده به گیرنده و 20 آمینواسید در ACE2 است که شبکهای متشکل از 14 پیوند هیدروژنی و یک پل نمکی را تشکیل میدهند (4). مطالعات نشان میدهند شبهویروسهای SARS-CoV-2 فاقد محل برش برای پروتئین اسپایک هستند؛ بنابراین، به نظر میرسد به دست آوردن محل برش فورین یکی از اولین رخدادهای سازگاری خانواده ویروسهای کرونا برای ایجاد عفونتهای شدید در انسان است (13). توالیهای دومین متصلشونده به گیرنده SARS-COV و SARS-CoV-2 درجه بالای تشابه 74 درصدی دارند؛ البته برخی از جایگزینیها در مناطق کلیدی اتصالی در دومین اتصال به گیرنده SARS-CoV-2 رخ داده است که به تماسهای اتمی بیشتر بین دو پروتئین منجر میشوند و بهطور بالقوه به میل اتصال بالاتر به ACE2 در SARS-CoV-2 منجر میشود (حدود 4 برابر). یک موتیف باقیمانده پلی بازیک در مرز بین زیرواحدهای S1 و S2 وجود دارد که آن را مستعد شکافتن توسط فورین / شبهفورین پروتئازها در طول بیوژنز و ورود به سلول میکند. پروتئازهای شبهفورین تقریباً در تمام بافتها و اندامهای مختلف توزیع میشوند؛ این امر توجیهکنندة علت انتقالپذیری بسیار کارآمد SARS-CoV-2 بین جمعیتهای انسانی است. همچنین تنوع پروتئازهای استفادهشده در ورود ویروس، علت دیگر این انتقالپذیری بالاتر و عفونت سیستماتیک چند عضوی در بیماران آلوده به SARS-CoV-2 است (2). دادههای حاصل از SDS-PAGE و وسترنبلات دربارة نحوه ورود خانواده کروناویروسها به سلول میزبان، شکست پروتئولیز اسپایک در هنگام ورود ویروس را تأیید میکنند و این موضوع را میتوان از باندهای تشکیلشده در حالتهای قبل و بعد از ورود ویروس مشاهده کرد (14و15). .ایجاد واریانتهای جدید SARS-CoV-2 با تکیه بر جهشهای پروتئین اسپایک: SARS-CoV-2 جزء معدود ویروسهای RNAدار است که مکانیسم تصحیح ژنوم[33] طبیعی دارد و این عملکرد را توسط nsp[34] با فعالیت َ3 به 5َ اگزونوکلئازی انجام میدهد (16). گزارشهای ژنومی نشان میدهند تاکنون هزاران جهش شناختهشده و ناشناخته در SARS-CoV-2 رخ دادهاند؛ اما به جهشهایی بیشتر توجه میشود که میتوانند عملکرد ورود ویروس به سلولهای میزبان را بهبود بخشند و همچنین بر اثرگذاری واکسنها تأثیر بگذارند (17). برای مدلسازی پروتئین اسپایک در واریانتهای جدید SARS-CoV-2، تمامی جهشهای صورتگرفته (اعم از جهشهای حذف و جانشینی) از واریانتهای مختلف که از مقالات مختلف استخراج شده بودند برای واریانتهای آلفا، بتا، دلتا و اومیکرون روی پروتئین اسپایک اعمال شدند که با رنگ مجزا از بقیه نواحی پروتئین در شکل 2 و 3 تفکیک شدند. درواقع این بررسیها که ترکیبی از تکیه بر دادههای توالی و ظاهر پروتئین اسپایک هستند، دید عمیقی برای ارزیابی جهشهای صورتگرفته در پروتئین اسپایک به ما میدهند. جهش مهم D614G(B.1): این جهش از نوع نقطهای بوده (18) که شامل جایگزینی آمینواسید گلایسین به جای آمینواسید آسپارتات در موقعیت 614 پروتئین اسپایک است. یافتهها نشان میدهند این تغییر میتواند برهمکنش پروتئین اسپایک و گیرنده ACE2 را افزایش دهد؛ البته این جهش با افزایش شدت عفونت ناشی از SARS-CoV-2 مرتبط نیست؛ اما میتواند میزان بار ویروسی واردشده به سلولهای انسانی را افزایش دهد. در ویروس پایه ووهان، وجود آسپارتات در موقعیت 614 باعث ایجاد پلهای نمکی بین این نقطه پروتئین اسپایک با K854 و T859آن میشود؛ درنتیجه به فشردهسازی این پروتئین منجر خواهد شد؛ اما با جایگزینی گلایسین در این موقعیت، پلهای نمکی دیگر وجود نخواهند داشت؛ بنابراین، تریمر پروتئین اسپایک با استحکام کمتری در کنار هم قرار خواهند گرفت. این ترکیب سستتر توجیهکنندة این واقعیت است که چرا واریانت حاوی D614G نسبت به نوع ووهان عفونیتر است (قسمت A شکل 2). درواقع این جهش موجب میشود پروتئین اسپایک برای تعامل با گیرنده ACE2 آزادتر باشد و راحتتر به گیرنده خود متصل شود. این جهش میزان انتقال در مقادیر کم ویروس را افزایش میدهد. D614G شکست فورین را کاهش میدهد و درنتیجه خطر ریزش زودرس S1 را کم میکند و پایداری حرارتی اسپایک را افزایش خواهد داد. این تغییرات درنهایت به افزایش بار ویروسی در دستگاه تنفسی فوقانی (بینی و حلق) و نه در ریهها منجر میشوند؛ به همین دلیل این جهش موجب تسهیل انتقال فرد به فرد میشود (3و8و18-20). مطالعات نشان دادند ویروسهای با پروتئینهای اسپایک حاوی G614 با کارایی بالاتری وارد سلولهای میزبان میشوند (21) و همچنین G614 علاوه بر بار ویروسی بالاتر با سن کمتر بیماران مرتبط است (22). این تغییرات ساختاری موجب شد در مدت بسیار کمی جهش D614G در غالب واریانتهای نوظهور کووید-19 نمایان شود و به نوعی تمامی واریانتهای جدید نوادگان این جهش نقطهای مهم بودند. ویروسشناسان بر این عقیدهاند که این جهش به تسریع فرایند تکامل SARS-CoV-2 منجر شده است و به همین دلیل جهش بسیار حائز اهمیتی است (19). ردیابی این واریانت جدید نشان داده است میزان بیان ژن ACE2 در جمعیتهای مختلف انسانی در ایجاد این جهش مؤثر بوده است؛ زیرا بنا بر بررسیهای انجامشده بیان این ژن در جمعیتهای اروپایی، آمریکای شمالی و آفریقایی بیشتر است و به همین دلیل این واریانت در این مناطق بیشتر مشاهده شده است؛ درحالیکه در چین به این شکل گسترده دیده نشد. درواقع این ایده در اینجا مطرح میشود که ممکن است سطح بیان ACE2 نیروی محرکهای برای جهش D614G بوده باشد (17). با توجه به اینکه بیشتر شرکتهای واکسنسازی پروتئین اسپایک را بهعنوان پروتئین هدف اصلی خود در توسعه واکسنها به کار بردهاند، این جهش در کاهش اثرگذاری واکسنها تأثیر مهمی داشته است (23). تکامل واریانتهای جدید: تغییرات انجامشده در پروتئین اسپایک تاکنون به افزایش انعطافپذیری ویروس منجر شده است. بنا بر ارزیابیهای انجامشده سرعت جهش در پروتئین اسپایک حدود سه برابر بیشتر از کل ژنوم SARS-CoV-2 است (شکل 2 و 3). براساس توصیه سازمان بهداشت جهانی، تمامی جهشهای نوظهور کشفشده توسط دانشمندان در پایگاه داده GISAID[35] (ابتکار جهانی برای به اشتراکگذاری تمامی دادههای آنفلوانزا پرندگان) قرار داده میشوند و درحال حاضر برای هر آمینواسید پروتئین اسپایک یک جهش در این پایگاه داده به ثبت رسیده است.
شکل 2- پروتئین اسپایک در واریانتهای مختلف SARS-CoV-2 (مناطق جهشیافته با رنگ زرد مشخص شدند). A) جهش D614G باعث بازشدن بیشتر پروتئین شده است (PDB ID 7BNM)، B) واریانت آلفا (B.1.1.7) (PDB ID 7LWU)، C) واریانت بتا (B.1.351) (PDB ID 7VX1)، D) واریانت دلتا (B.1.617.2) (PDB ID 7SBK) ، D) نقاط جهش واریانت اومیکرون در مقایسه با ویروس ووهان (PDB ID 6VYB)، G) پروتئین اسپایک اولیه و فاقد جهش؛ دومین اتصال به گیرنده RBD (زرد)، زیرواحد S1 پروتئین اسپایک (قرمز)، زیرواحد S2 پروتئین اسپایک (آبی) (PDB ID 6VXX).
شکل 3- تصویر از دید بالا پروتئین اسپایک در واریانتهای جدیدتر اومیکرون و دلتا. تعداد جهشها در اومیکرون بهخصوص در زیرواحد S1 که مسئول اتصال پروتئین به گیرنده را بر عهده دارد، به شکل قابل ملاحظهای بالاتر بوده است و به نوعی میتوان این را بازطراحی دوباره پروتئین اسپایک تعبیر کرد (PDB ID 6VYB) و (PDB ID 7SBK).
تنها تعداد کمی از این جهشها میتوانند مزیتهای فرار از سیستم ایمنی، واکسنها و آنتیبادیهای خنثیکننده را به ویروس بدهند (19)؛ بنابراین، بین واریانتهای جدید تنها جهشهایی که اهمیت بالینی و نگرانکننده داشته باشند، عنوان VOC[36] را به خود اختصاص خواهند داد. علاوه بر این، واریانتهای با اهمیت کمتر با عنوان VOI[37] شناخته میشوند (24و25). مشخص شده است فرکانس نواسانات مولکولی ارتباط نزدیکی با خواص عملکردی مولکولهای زیستی دارند. کاهش نوسانات مولکولی در ناحیه اتصال به گیرنده گلیکوپروتئین اسپایک جهشیافته را میتوان از دو جنبه ارزیابی کرد. نخست، کاهش حرکت مولکولی به کاهش میل ترکیبی ACE2 میتواند منجر شود و دوم، کاهش حرکت مولکولی میتواند کمپلکس اسپایک-ACE2 را قادر کند پس از اتصال به ساختاری پایدارتر و محکمتر تبدیل شود. درواقع دو حالت up-formation و down-formation در اینجا مطرح است؛ در حالت up-formation پایداری مولکولی و ساختاری محکم روی خواهد داد که به اثر منفی بر اتصال بین اسپایک و ACE2 منجر خواهد شد. درواقع این یافتهها میتوانند درک عمیقی را دربارة شیوه پیشبینی و ارزیابی عملکرد هر جهش جدید روی عملکرد پروتئین اسپایک به ما بدهد و میتواند در کنترل پاندمی کمک شایانی به جامعه بشری ارائه کند (26). همانطور که در جدول 1 مشاهده میشود واریانتهای اپسیلون، اتا، لوتا، کاپا و لامبدا روی عملکرد آنتیبادیهای خنثیکننده تأثیرکمتری دارند و به همین دلیل با عنوان VOI شناخته میشوند (13). داکینگ مولکولی دادههای پروتئینی اسپایک واریانتهای مختلف ویروس با گیرندههای مختلف سلولی در بدن انسان: شایان ذکر است علت استفاده از توالی پروتئین اسپایک چهار واریانت اصلی، آلفا، دلتا و اومیکرون اهمیت بالینی بیشتر (VOC) و تأثیر بیشتر این جهشها در تداخل در عملکرد آنتیبادیهای خنثیکنندة این پروتئینهای اسپایک به نسبت سایر پروتئینهای اسپایک واریانتهای دیگر کمتر مهم (VOI) بوده است. پس از انجام فرایند داکینگ مولکولی توسط وب سرور HDOCK مدل اول و محتملتر هرکدام از محاسبات با نرمافزار PyMOL ترسیم شد. همانطور که در جدول 2 نمایان است گیرنده مدنظر با رنگ سبز، زیرواحد S1 با رنگ قرمز و زیرواحد S2 با رنگ آبی مشخص شدند. سطح انرژی یا امتیاز داکینگ منفیتر نشاندهندة پایداری بالاتر پیشبینیهای انجامگرفته است؛ درنتیجه میتوان گفت این عدد با تمایل اتصال دو جزء به هم ارتباط دارد. مقایسه سطح انرژی دریافتی از هرکدام از محسابات نشاندهندة تمایل بالای جهشهای جدید SARS-CoV-2 برای اتصال به گیرنده KREMEN1 بود که این تمایل اتصال حتی از گیرنده اصلی ACE2 نیز بیشتر بود. این در حالی است که هیچگونه تفاوت معنیداری در سطح انرژی اتصال گیرنده ACE2 با واریانتهای جدید پروتئین اسپایک در مقایسه با سایر گیرندهها مشاهده نشد و درواقع تمایل اتصال تفاوت چشمگیری نداشت. همانطور که مشخص است بهعلت منفیتربودن قابل ملاحظه و محسوس گفتنی است تمایل اتصال پروتئین اسپایک واریانتهای نوظهور کووید-19 نسبت به گیرنده جایگزین و غیراصلی KREMEN1 درحال افزایش است؛ بهطوریکه هرچه از واریانت اصلی به سمت دلتا پیش میرویم، تمایل اتصال به KREMEN1 (از 27/308- تا 51/388- کیلوژول بر مول) افزایش یافته است. پایداری ساختار بهوجودآمده در واریانت اومیکرون با KREMEN1 با کاهش محسوسی در سطح انرژی (06/323- کیلوژول بر مول) همراه بوده است؛ درحالیکه تفاوت محسوسی در تمایل اتصال پروتئین اسپایک به گیرنده ACE2 مشاهده نشده است. گفتنی است با وجود اینکه ما در این مطالعه از تکنیک داکینگ سرتاسری استفاده کردیم، باز هم زیرواحد S1 پروتئین اسپایک محل اصلی اتصال بین پروتئین اسپایک و گیرندههای آن بوده است؛ درحالیکه زیرواحد S2 که در اصل در داخل پوشش ویروس قرار میگیرد و قاعدتاً نقشی در اتصال ندارد، بنا بر مطالعات داکینگ، حتی در صورتی که به شکل خالص و جدا از ساختمان ویروس نیز بررسی شود، نقشی در اتصال ویروس نمیتواند داشته باشد. علت این امر مناسببودن ریشههای موجود در زیرواحد S1 برای اتصال به گیرنده بوده است و این اتصالات با زیرواحد S1 کمترین انرژی ممکن را دارند و پایدارترند؛ بنابراین، محتملتر نیز هستند. این اتصالات در تصاویر جدول 2 نمایاناند. برای اطمینان از انجام درست و بیعیب فرایند داکینگ مولکولی از وب سرور PDBsum میتوان استفاده کرد. تجزیه و تحلیلهای آماری که این وب سرور در اختیار قرار داد بررسی شدند. میزان G-Factor و نمودار راماچاندران در اولویت بودند. در بررسیهای انجامشده G-Factor بسیار غیرطبیعی (کمتر از 0/1-) در مورد هیچکدام از اتصالات جدول 2 مشاهده نشد. در کل با آنالیز نمودارهای راماچاندران و G-Factorها به این نتیجه رسیدیم که فرایند داکینگ مولکولی به خوبی صورت گرفته است و خطای سیستماتیک و بزرگی در تحلیل نتایج مشاهده نشد. .داروها و واکسنهای استفادهشده برای درمان کووید-19: داروهای استفادهشده برای درمان کووید-19 به سه دسته تقسیمبندی میشوند: 1) آنتیبادیهای خنثیکننده، 2) مهارکنندههای فیوژن و 3) مهارکنندههای پروتئازها. علاوه بر این داروها، داروهای ضدویروس غیراختصاصی مانند IFN-alpha نوترکیب انسانی، رمدیسیویر (ضدابولا و عملکرد آن مهار RNA پلیمرازها است)، کلروکین (تولید سایتوکینهای التهابی را کاهش میدهد)، فاویپیراویر (داروی ضدآنفلوانزا)، لوپیناویر - ریتوناویر (مهارکنندة پروتئازها) و آرابیدول (داروی ضدآنفلوانزا) ازنظر بالینی برای درمان کووید-19 استفاده شدهاند (1و2). همچنین هدف قراردادن فعل و انفعالات آلوستریک آنزیمهای دخیل در مکانیسم ویروس یک استراتژی پذیرفتنی را برای مداخله درمانی ارائه دادند (41). گفتنی است آنتیبادیهای مونوکلونال بهدستآمده از اشخاص بهبودیافته از کووید-19 میتوانند به اپیتوپهای جدیدی ازجمله حفرههای بسیار داخلی پروتئین اسپایک حمله کنند (4). علاوه بر داروها، واکسنها نیز برای کنترل پاندمی کووید-19 استفاده شدند و پلتفرمهای مختلفی برای ساخت واکسنها به کار رفتند؛ ازجمله 1) ویروس کامل غیرفعال، 2) پروتئین زیرواحد ویروسی، 3) استفاده از حامل آدنوویروس و 4) استفاده از قطعات mRNA پوشیدهشده با لیپوزوم (13). دلیل این تنوع گسترده در پلتفرمهای ساخت واکسن تلاش برای کنترل پاندمی SARS-CoV-2 با کمک واکسنهای ایمن و مؤثر بود (42). پروتئین اسپایک با توجه به نقش مهمی که در فرایند ورود ویروس به سلولهای میزبان بازی میکند، هدف اصلی برای ساخت واکسن و آنتیبادیهای خنثیکننده و مهارکنندههای پروتئازی است (1و2و8). .ارزیابی اثر داروها بر واریانتهای SARS-CoV-2: اثرات جهش بر آنتیبادیها ارزیابی شدهاند. مطالعات نشان میدهند جهشهای دومین اتصال به گیرنده (RBD) پروتئین اسپایک میتوانند بر برهمکنشهای آنتیژن - آنتیبادیهای استفادهشده تأثیرگذار باشند (43)؛ برای مثال، بنا بر ارزیابیهای انجامگرفته در مورد واریانت اومیکرون و اعلام سازمان بهداشت جهانی، کورتیکواستروئیدها و مهارکنندههای گیرنده اینترلوکین 6 همچنان مؤثرند؛ اما دربارة سایر داروها باید اثرگذاری آنها بهطور مجدد ارزیابی شود (39). .ارزیابی اثر واکسنها بر واریانتهای SARS-CoV-2: با توجه به ظهور واریانتهای جدید SARS-CoV-2 و اینکه هنوز افراد واکسینهشده بخش کوچکی از جمعیت جهانی را تشکیل میدهند، نمیتوان با استناد به میزان مرگومیر و تعداد مبتلایان به ارزیابی دقیقی از میزان اثرگذاری واکسنها دست یافت؛ اما با استفاده از آزمایش میزان خنثیسازی سرم خون افراد واکسینهشده در برابر انواع واریانتهای نوظهور میتوان ارزیابی تقریباً مناسبی از میزان کارایی واکسنهای مختلف داشت (18و43و44). دو ابزار آزمایشگاهی که در ارزیابی عملکرد واکسنها در مقابل انواع جهشهای SARS-CoV-2 استفاده میشوند (21)، شامل شبهویروسها و ذرات شبهویروسی هستند. شبهویروسها هترولوگهایی هستند که با پروتئینهای ساختاری ایمن و غیرقابل تکثیر ویروس مدنظر پوشیده شدند؛ اما ذرات شبهویروسی همان پروتئینهای ساختاری ویروس مطالعهشده هستند و بنابراین، خواص مشابهتری نسبت به ویروسهایی دارند که از آنها مشتق شدند و میتوان از آنها برای مطالعه کل مجموعه ویروسها استفاده کرد (13). علاوه بر روشهای خنثیسازی سرم میتوان از تکنیک بیوانفورماتیکی ISM[38] برای تشخیص اثربخشی انواع واکسنهای توسعهدادهشده در برابر واریانتهای نوظهور کووید-19 استفاده کرد. این روش قبلاً هر ساله برای ساخت واکسنهای فصلی آنفلوانزا استفاده میشد. بیشترین کاربرد ISM برای شناسایی جهشهای مهمی است که میتوانند بر اثربخشی واکسنها علیه ویروس تأثیرگذار باشند (42). همانطور که از نتایج فرایند داکینگ مولکولی در جدول 2 مشخص است، واریانتهای مختلف کویید-19 در تمایل اتصال این ویروس به گیرندهها میتوانند نقش داشته باشند و تمایل اتصال از واریانت اصلی تا اومیکرون به گیرندههای مختلف ویروس در بدن انسان افزایش یافته است؛ این دادهها میتوانند تأییدکنندة عملکرد خنثیسازی واکسنهای مختلف علیه پروتئین اسپایک باشند که در جدول خنثی شدهاند. درواقع کاهش در عملکرد خنثیسازی واکسنها (جدول 3) به نوعی میتواند مربوط به افزایش تمایلات اتصال پروتئین اسپایک در واریانتهای VOC با گیرندههای آن باشد.
جدول 1- جهشهای با اهمیت پروتئین اسپایک در واریانتهای SARS-CoV-2
جدول 2- نتایج داکینگ مولکولی بین پروتئین اسپایک اصلی و واریانتهای آلفا، دلتا و اومیکرون SARS-CoV-2 با گیرندههای مختلف سلولی در بدن انسان به همراه سطح انرژی اتصال و درجه انطباق پروتئین اسپایک و گیرنده با هم (rmsd)
* سطح انرژی پایینتر و درنتیجه ایجاد ساختار پایدارتر ساختار متصل پروتئین اسپایک با KREMEN1 مشهود است و در واریانت دلتا به اوج خود رسیده است. * rmsd: درجه انطباق دو جزء به هم متصلشده که این مقدار عددی هرچه به صفر نزدیکتر باشد، یعنی دو جزء (برای مثال پروتئین اسپایک و گیرنده) بهتر روی هم منطبق شدهاند. جدول 3- اثربخشی واکسنهای معروف در جهان پس از 14 روز از تزریق دو یا سه دوز در برابر انواع واریانتهای مهم SARS-CoV-2 خلاصه شده است.
جمعبندی مطالعه حاضر با ارائه نکات مهم در زمینة پروتئین اسپایک SARS-CoV-2، گامی در جهت درک بهتر و مبارزه با کووید-19 و همچنین بیماریهای نوظهور ویروسی در آینده است. بدون شک انجام آزمایشات جدید برای به دست آوردن بینش بهتر در مورد SARS-CoV-2 لازم است و این امر برای درک اثر واریانتهای نوظهور و استراتژیهای جدید برای مبارزه با این ویروس بسیار با اهمیت است. درواقع تغییرات انجامشده در ساختار SARS-CoV-2، بهخصوص پروتئین اسپایک باید پیوسته، ارزیابی و با ویروس اولیه و اصلی (ووهان) مقایسه شود. همچنین باید تحقیقات بیشتری در زمینة مکانیسمهای ورود این ویروس به سلولهای میزبان صورت بگیرد تا مشخص شود چه گیرندهها و پروتئازهای جایگزین دیگری میتوانند در مکانیسم ورود ویروس به سلولهای میزبان نقش داشته باشند. بدون شک درگیری ریه با وجود اینکه آنها ازجمله بافتهاییاند که پایینترین میزان بیان گیرنده ACE2 را بین بافتهای دیگر دارند، شاهد بزرگی بر این ادعا است؛ امید است با پژوهشهایی که در آینده در این زمینه انجام میگیرند، به این سؤال اساسی پاسخ داده شود. همانطور که در این مطالعه به آن اشاره شد، بنا بر مطالعات داکینگ مولکولی انجامگرفته در این مقاله تفاوت زیادی در تمایل اتصال پروتئین اسپایک واریانتهای نوظهور و مهم ویروس (اصلی، آلفا، دلتا و اومیکرون) با گیرنده ACE2 وجود ندارد؛ اما میزان تمایل اتصال به گیرنده جایگزین KREMEN1 رو به افزایش بوده و در واریانت دلتا به اوج خود رسیده است و سطح انرژی اتصال پروتئین اسپایک به گیرنده KREMEN1 بهطور محسوسی به نسبت سایر گیرندهها منفی بوده و این نشاندهندة پایداری بیشتر و احتمال بیشتر ایجاد این اتصالات بوده است که این موارد در جدول 2 مشاهده میشوند. همچنین با توجه به سطح انرژیهای دریافتی از وب سرور HDOCK، میتوان به این نتیجه رسید که اگرچه گیرنده ACE2 گیرنده اصلی پروتئین اسپایک نام گرفته است، همچنان گیرندههای دیگری با پایداری خوب میتوانند به پروتئین اسپایک متصل شوند. براساس اطلاعات و دانش ما و بنا بر بررسیهای صورتگرفته، تاکنون هیچ گزارشی مبنی بر انجام فرایند داکینگ بین پروتئین اسپایک با سایر گیرندههای غیراصلی و جایگزین پروتئین اسپایک (AXL، ASGR1 و KREMEN1) انجام نگرفته است (تا آوریل 2022)؛ به همین دلیل باید برای تأیید یا رد این یافته پژوهشهای گستردهتر و بیشتری انجام شود. یکی از چالشهای مهم و اساسی، ارزیابی میزان اثرگذاری واکسنها بر واریانتهای نوظهور است که متأسفانه در این زمینه دادههای متناقض فراوانی وجود دارند؛ البته تحقیقات گسترده دانشمندان در جهان در زمینة میزان اثربخشی واقعی واکسنها بر واریانتهای نوظهور درحال توسعه است که در صورت نیاز به واکسنهای جدید برای واریانتهای جدید به کار خواهند رفت. واکسنهای استفادهشده در جهان در مقابل واریانتهای جدید و مهم کارایی لازم را دارند؛ اما تقریباً کارایی همگی آنها به مرور و با ظهور واریانتهای جدید درحال کاهش است. یکی دیگر از نکات اساسی، شناخت دقیق نقاط مستعد جهش در پروتئین اسپایک ویروس SARS-CoV-2 است که با توجه به این نقاط مستعد جهش، با ارزیابیهای بیوانفورماتیکی جهشهای مهمی که ممکن است در آینده در پروتئین اسپایک رخ دهد را میتوان پیشبینی کرد. درواقع دانشمندان علاوه بر ارزیابی و ردیابی مداوم تمامی واریانتهای جدید از نمونههای گزارششده در جهان، باید به دنبال پیشبینی انواع جهشهای جدید و خطرناک بودن یا نبودن آنها بروند تا بتوان حتی قبل از ظهور جهشهای جدید خود را برای مقابله با آنها آماده کرد و در صورت لزوم واکسنهای جدیدی را در این زمینه طراحی و توسعه داد.
[1]- Coronavirus disease 2019 [2]- COVID-19 [3]- Angiotensin-converting enzyme 2 (encoded by ACE2 gene) [4]- Tyrosine-protein kinase receptor UFO (encoded by the AXL gene) [5]- Kringle-containing protein marking the eye and nose protein1(Kremen protein 1 in humans is encoded by the KREMEN1 gene) [6]- Asialogylcoprotein receptor 1 (encoded by the ASGR1 gene) [7]- Middle East respiratory syndrome coronavirus (MERS-CoV) [8]- SARS [9]- MERS [10]- Middle East respiratory syndrome [11]- Protein Data Bank (PDB) [12]- Preprocessing [13]- Molecular docking [14]- Global Docking [15]- http://hdock.phys.hust.edu.cn/ [16]- Postprocessing [17]- Root-mean-square deviation (RMSD) [18]- https://www.ebi.ac.uk/thornton-srv/databases/cgi-bin/pdbsum/GetPage.pl?pdbcode=index.html [19]- Receptor-binding domain (RBD) [20]- Global Initiative on Sharing Avian Influenza Data (GISAID) [21]- Signal sequence [22]- N-terminal domain [23]- Fusion peptide [24]- heptad-repeat region 1 [25]- heptad-repeat region 2 [26]- Transmembrane [27]- Prefusion [28]- Transmembrane serine protease 2 [29]- Serin endoprotease proprotein convertase1 [30]- Trypsin [31]- Trypsin-like integral membrane serine peptidase [32]- Cathepsin [33]- Proofreading [34]- nonstructural protein (nsp) [35]- Global Initiative on Sharing Avian Influenza Data (GISAID) [36]- Variants of Concern (VOC) [37]- Variant of Interest (VOI) [38]- Informational spectrum method (ISM) | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
مراجع | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
آمار تعداد مشاهده مقاله: 2,628 تعداد دریافت فایل اصل مقاله: 305 |