بررسی ساختار و عملکرد پروتئین اسپایک کروناویروس سندرم حاد تنفسی 2 (sars-cov-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) صورت نگرفته است؛ بنابراین، برای تایید یا رد این یافته‌ها باید پژوهش‌های جامع‌تری انجام شود.

investigating the structure and function of sars-cov-2 spike protein by molecular docking

introduction: it has been more than two years since the outbreak of severe acute respiratory syndrome coronavirus 2 (sars-cov-2) in wuhan, china. the sars-cov-2 has created an unprecedented pandemic (covid-19) in the modern era and challenged all the scientific, economic, and social structures of the world. extensive worldwide research has been conducted to study the structure and function of the sars-cov-2 spike protein since the beginning of the pandemic. the role and importance of spike protein in the introduction of sars-cov-2 into human cells has been proven. among the important and effective factors in controlling the covid-19 pandemic is understanding the structural changes of spike protein in different variants of the virus and the effectiveness of drugs and vaccines against new sars-cov-2 variants. due to the great importance of spike protein sars-cov-2, it was tried to study its structure and interaction with cellular receptors using bioinformatics tools in this study. furthermore, the effects of different sars-cov-2 variants on existing drugs and vaccines were reviewed.materials and methods: the bioinformatics methods and molecular docking was used to study the binding function of spike protein to its various cellular receptors in the human body. hdock web server was used to perform the molecular docking process and in the next step, the pdbsum web server was used for post-docking checks and to ensure the accuracy of the obtained data. pymol software was also used to study and visualize replacement and deletion mutations in the spike protein.results: extensive mutations in the spike protein in the omicron variant indicate a redesign in the spike protein. these massive mutations cause extensive changes in the structure and function of the spike protein. the data obtained from molecular docking and comparing the energy level (docking score) of spike protein binding of different variants of the virus with different cellular receptors indicate a high tendency of the spike protein to bind, with lower energy levels and more stable, to kremen1 receptor than the main ace2 receptor.  the level of spike protein binding energy levels of different virus variants was almost similar to other receptors (except the kremen1 receptor) or slightly different from the spike protein binding energy levels with the main ace2 receptor.discussion and conclusion: there is not much difference in the tendency of spike protein binding of emerging and important variants of the virus (main, alpha, delta, and omicron) to ace2 receptor, but the tendency to bind to kremen1 receptor is increasing and has reached its peak in the delta variant. the energy level of the spike protein binding to the kremen1 receptor is significantly negative compared to other receptors, indicating greater binding and stability. also, considering the level of energy received from the hdock web server, it can be concluded that although the ace2 receptor is the main receptor for the spike protein, other receptors can still bind to the spike protein with good stability. based on the available information, so far there have been no reports of docking between spike protein and other non-main spike protein receptors (axl, asgr1, kremen1). therefore, more comprehensive research is needed to confirm or refute these findings.