نوع مقاله : مروری
نویسندگان
1 گروه بیوتکنولوژی، دانشکده علوم پایه، واحد شهرکرد، دانشگاه آزاداسلامی، شهرکرد، ایران
2 گروه زیست شناسی، دانشکده علوم پایه، واحد رشت، دانشگاه آزاد اسلامی، رشت، ایران
چکیده
کلیدواژهها
Introduction
On December 31, 2019, an unknown coronavirus disease was reported from China that cause acute respiratory syndrome which was named COVID-19. This virus was highly contagious and turned into a pandemic in just 3 months. It seems that the high tendency to evolve and the similarity to other viruses of the beta-coronavirus family are the main factors in the susceptibility of this virus to mutation. In all genomes of this virus, the tendency of the virus to spread are varied. However, they are common in human-to-human transmission. The virus can infect one person quickly but not infect others. It can be said that this virus is selective in choosing its host, but the reports have indicated that it has a high rate of 90% in accepting a new host. The angiotensin converting enzyme-2 (ACE-2) receptor is the receptor of SARS-COV-2 to enter human cells. Since the expression level of ACE2 in men living in East Asian countries is higher than in other countries, this population is more susceptible to the virus. In this review study, firs the molecular and cellular characteristics, proliferation, and replication of SARS-COV-2 are investigated. Then, by reviewing bioinformatics studies conducted regarding the analysis of conserved domain and homology and molecular docking, the biological roles of some specific proteins SARS-COV-2 are investigated.
Methods
In this study, 60 related articles were reviewed. Then, the information extracted from them were categorized including the types of SARS-COV-2 genomes in different parts of the world and the identified proteins and structural components of these viruses. In this regard, we first reviewed the studies conducted on the genome structure of the SARS-COV-2 and the genetic affinity of SARS-COV-2 to the known coronaviruses. Then, the function of structural and non-structural proteins of coronaviruses was reviewed. Finally, the bioinformatics studies on open reading frames (ORFs) of the virus and the possible origin of SARS-COV-2 were reviewed to know whether it is possible that this virus was created in laboratory or not. This section subdivided into two sections: (a) Bioinformatic study regarding that ORF8 and glycoprotein inhibit Heme metabolism by binding to porphyrin, and (b) bioinformatic study of the genome structure of SARS-COV-2 by examining the mutation in the receptor-binding domain of SARS-COV-2 and the polybasic cleavage site and O-glycan linkage.
Results
Analysis of SARS-COV-2 sequences showed that this virus had a typical genomic structure of beta-coronaviruses, including: Bat SARS-like-ZC45, Bat SARS-like-ZXC21, SARS-CoV and MERS-CoV. According to the phylogenetic tree of coronaviruses (Figure 1), SARS-COV-2 had a high genetic similarity to Bat SARS-like-ZC45 and Bat SARS-like-ZXC21, while SARS-CoV had the lowest similarity to them.
In this study, it was found that ORF8 and surface glycoprotein can bind to porphyrin, and ORF1ab, ORF10 and ORF3a can attack the heme part of hemoglobin to dissociate iron and form porphyrin. This attack causes a decrease in the hemoglobin to carry oxygen and carbon dioxide. Lung cells become inflamed due to the inability to exchange carbon dioxide and oxygen, which leads to large ground-glass opacities in CT scan images.
The three-dimensional image of the structures of three viral proteins (surface glycoprotein, envelope protein, and nucleocapsid phosphoprotein) was downloaded from the PDB database. The structural model indicated the inability of the surface glycoprotein to bind to the porphyrin in Heme, while the envelope protein could bind to porphyrin. For the phosphoprotein nucleocapsid, the results indicated the ability of this protein to bind to porphyrin with the highest possible binding energy (5206.53 kcal/mol). Figure 2 shows the docking results, a two-dimensional view of the binding sections of the envelope protein and nucleocapsid phosphoprotein with human porphyrin.
Figure 3 illustrated the two-dimensional view of the binding sections of viral non-structural protein ORF8 with human porphyrin.
Discussion
Based on the structural studies performed on SARS-COV-2, it seems that this virus can bind to the human cell receptor ACE-2. The S protein of SARS-COV-2 has a functional polybasic cleavage site called “Furin” at the S1/S2 boundary, which is formed by 12 nucleotides, resulting in three O-glycan binding sites around this polybasic cleavage site. Since the COVID-19 pandemic is a global problem, predictive methods to search for the protected domains of the virus are very important. Achieving the structure of protein molecules including ORF8 and nucleocapsid phosphoproteins (surface proteins) using bioinformatics study and analyzing the possibility and the method of their binding to human porphyrin in the heme part of hemoglobin is needed. According to the review of the bioinformatics studies, chloroquine can prevent ORF1ab, ORF10 and ORF3a from attacking hemoglobin to form porphyrin, and avoid the binding of ORF8 and surface glycoproteins to porphyrin to a certain extent, which effectively relieves acute respiratory symptoms. The results of bioinformatics studies are of high value in preventing the spread of COVID-19, developing drugs and vaccines, and clinical practice.
Ethical Considerations
Compliance with ethical guidelines
This research was written in accordance with the principles of research ethics approved by Islamic Azad University, Shahrekord Branch.
Funding
This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.
Authors' contributions
Conceptualization: Tohid Piri-Gharaghie; Research and review: Sheida Beiranvand; Editing and finalization written by: Sameh Haji Mohammadi; Amirhossein Ghadiri
Conflicts of interest
The authors declared no conflict of interest.
Acknowledgements
The authors would like to thank the staff members of the Biotechnology Research Center of the Islamic Azad University of Shahrekord Branch in Iran for their help and support.
مقدمه
در 31 دسامبر 2019، یک ویروس کرونا ناشناخته، از چین بهعنوان عامل ایجاد عفونتهای حاد تنفسی گزارش شد. این ویروس، قدرت سرایت بسیار بالایی داشت، بهطوریکه تنها در مدت 3 ماه، به یک همهگیری بزگ جهانی تبدیل شد [1]. ویروسکرونایی که به تازگی کشف شده است، از نظر ژنتیکی با دیگر کروناویروسها مثل سارس و مرس قرابت ژنتیکی دارد. بهطور سادهتر، میتوان گفت که این ویروس، از خانواده بتاکروناویروسها است. تعدادی از ژنومهای مهم شناساییشده از این گروه جدید، ایزوله و جدا شده است که در جدول شماره 1 انواع ژنوم مهم ویروس کرونا که تا امروز گزارش شده، نشان داده شده است [2]. به نظر میرسد تمایل بالا برای به تکامل رسیدن و شباهت داشتن به دیگر ویروسهای خانواده، بهعنوان عوامل اصلی مستعد بودن این ویروس برای جهش باشد که به دنبال آن، تمایل ویروس به انجام تغییرات در اثر جهش، باعث بروز و شکلگیری این ویروس جدید شده است [3].
در تمامی این ژنومهای یافتشده، تمایل ویروس به سرایت از یک انسان به انسان دیگر متغیر بوده است. البته مواردی نیز وجود داشته است که افراد ناقل، ویروس را به دیگران منتقل نکردهاند؛ در مواردی هم انتقال فوقالعاده بالای عفونت از فردی به فرد دیگر مشاهده شده است. بنابراین میتوان گفت ویروس از نظر انتقال از انسان به انسان، در هر دو حالت افراط دارد. این بدان معنی است که این ویروس میتواند یک فرد را به سرعت آلوده کند، اما فرد دیگری را مبتلا نکند. ازاینرو شاید بتوان گفت که این ویروس، در انتخاب میزبان خود، گزینشی عمل میکند.گزارشها حاکی از آن است که ویروس، نرخ بالای90 درصد را در پذیرش میزبان جدید از خود نشان میدهد و این بهمعنای همهگیری آن است [3]. یک توافق کلی در مورد گیرنده اختصاصی آنژیوناسین2 بهعنوان گیرنده ویروس کرونا جهت ورود به سلولهای انسانی وجود دارد. ازآنجاییکه میزان بیان آنژیوناسین2 در مردان شرقی، بیشتر از دیگر جمعیتهای انسانی است، ازاینرو این جمعیت برای ورود ویروس به بدن بسیار مستعد هستند [4].
ساختار ژنوم ویروس کرونا جدید
ظهور ویروس کرونای جدید به نامSARS-COV-2 که باعث شیوع پنومونی ذاتالریه غیرمعمول در دنیا شده است، چالش بزرگی را پیش روی محققان قرار داد. چن و همکاران باتوجهبه جهش این ویروس و جدید بودن آن، به تشریح ساختار ژنوم این ویروس پرداختند و مشخص کردند SARS-COV-2 مربوط به زیرخانواده کروناویرینا از خانواده کروناویریده از زیر رده نیدوویرالس است و این زیرخانواده شامل 4 جنس: گاماکروناویروس، بتاکروناویروس، آلفاکروناویروس و دلتا کروناویروس است [5].
ژنوم ویروس کرونا RNA تک رشته سنس تقریباً24 تا 34 کیلو بازی با ساختار سر 5’ cap و دم 3’ پلی A است. RNA ژنومی، در فرایند ترجمه مستقیماً برای ترجمه پلی پروتئین 1a/1ab مورد استفاده قرار میگیرد که پروتئینهای غیرساختاری را برای تشکیل کمپلکس همانندسازی-رونویسی در یک ویزیکول ۲غشایی کد میکند [6]. به دنبال آن، مجموعهای از RNAهای سابژنومیک توسط کمپلکس همانندسازی_رونویسی به شکل درونزا به روش رونویسی ناپیوسته سنتز میشوند. این mRNAهای سابژنومیک توالیهای 5’ پیشرو و انتهاهای3’ مشترک دارند. خاتمه رونویسی و متعاقباً ایجاد یک RNA پیشرو در توالی نظارتی رونویسی که در بین ORFها قرار دارد، رخ میدهد. این RNAهای سابژنومیک رشته منفی بهعنوان الگو برای تولید mRNAهای سابژنومیک عمل میکنند. ژنوم و زیرژنوم SARS-COV-2 حداقل 6 نوع بازخوانی یا ORF دارد [5، 6]. اولین ORF که ORF1a/b است، حدود دوسوم طول ژنوم بوده که 16 پروتئین غیرساختاری یا nsp1-16 را کد میکند. یک جهش نقطهای بین ORF1a و ORF1b منجر به تولید دو پلیپپتید pp1a و pp1ab میشود. این پلیپپتیدها به روش ویروسی توسط پروتئازهای شبه کیموتریپسین یا پروتئازهای اصلی و یک یا دو عدد پروتئاز شبهپاپائین در nsps16 پردازش میشوند.
ORFهای دیگر، در یکسوم ژنوم که در نزدیکی انتهای 3’ قرار گرفته است، حداقل 4 پروتئین ساختاری اصلی را کدگذاری میکند که شامل پروتئینهای اسپایک، غشا (M)، پوشش (E)، و پروتئینهای نوکلئوکپسید است. علاوهبراین، 4 پروتئین ساختاری اصلی، پروتئینهای ساختاری خاص و همچنین پروتئینهای جانبی دیگر مانند پروتئین HE، پروتئین 3a/b، پروتئین 4a/b و دیگر پروتئینها را کد میکنند [6]. تمام پروتئینهای ساختاری و جانبی از RNAهای سابژنومیک SARS-COV-2 ترجمه شده است [5]. نرخ جهش در همانندسازی ویروسهای RNAدار بسیار بیشتر از ویروسهای DNAدار است و ژنوم ویروسهای RNAدار معمولاً کمتر از 10K نوکلئوتید است. اندازه ژنوم این ویروس 30Kb بوده که در بین ویروسهای RNAدار دارای بزرگترین ژنوم هستند. اندازه این ژنوم، تقریباً 2 برابر بزرگتر از سایر ویروسهای RNAدار است و دومین ویروس RNAدار بزرگ ازنظر اندازه ژنوم به شمار میرود.
در ویروسهای RNAدار، نگهداری ژنومی با این سایز بزرگ، نیازمند ویژگیهایی خاص است که در ویروس کروناهای کمپلکس با چندین RNA آنزیمهایی پردازشکننده مانند 3’ و 5’ اگزوریبونوکلوئاز از nsp14 را بیان میکند که مسئول این عملکرد است. آنزیم 3’ و 5’ اگزوریبونوکلوئاز در بین تمام ویروسهای RNAدار مختص کروناویروسها بوده که بهعنوان بخشی از قسمت تصحیحکننده از RTC عمل میکند [6]. تجزیهوتحلیل توالیهای SARS-COV-2 توسط چن و همکاران، نشان داد که این ویروس دارای یک ساختار ژنومی معمولی از ویروس کروناهای متعلق به بتاکروناویروسها است که از Bat-SARS-like(SL)-ZC45، Bat-SL ZXC21، SARS-CoV و MERS-CoV تشکیل شده است. مطابق با درخت فیلوژنی کروناویروسها، SARS-COV-2 شباهت بسیار زیاد ژنتیکی به Bat-SL-CoV ZC45 وBat-SL-CoV ZXC21 دارد و در بین این نمونهها٬ کمترین شباهت را به SARS-CoV دارد. تصویر شماره 1 درخت فیلوژنی٬ قرابت ژنی این نمونهها را بهوضوح نشان میدهد [5، 6].
عملکرد پروتئینهای ساختاری و غیرساختاری ویروسهای کرونا با تأکید بر SARS-COV-2
اکثر پروتئینهای غیرساختاری که در همانندسازی کرونا ویروسها نقش کلیدی دارند، گزارش شدهاند. بااینحال، عملکرد برخی از nspsها هنوز ناشناخته هستند یا کاملاً مشخص نشدهاند [7]. چن و همکاران بهطور خلاصه، عملکرد پروتئینهای غیرساختاری کروناویروسها را توضیح دادهاند که در جدول شماره 2 نمایش داده شده است. علاوهبراین، 4 پروتئین ساختاری برای مونتاژ و ایجاد ویریون و عفونت ویروسهای کرونا ضروری است [۵]. هموتریمرهایی از S پروتئینها باعث افزایش پروتئینهای اسپایک در سطح ذرات ویروس میشود و این بهعنوان کلید اتصال ویروس به گیرندههای سلول میزبان عمل میکند. پروتئین M 3 دامین انتقال غشایی دارد که باعث شکلگیری ویریونها و خمیدگی غشایی میشود و به نوکلئوکپسید متصل میشود. پروتئین E در مونتاژ و رهاسازی ویروس نقش دارد و برای پاتوژنز ویروس لازم است. پروتئین N شامل 2 دامین است که هر دوی آنها با مکانیسمهای مختلف به ژنوم RNA ویروس متصل میشوند [۸].
همچنین گزارش شده که پروتئین N میتواند با اتصال به پروتئین غیرساختاری 3 و کمپلکس٬ همانندسازی_رونویسی را فعال کند و به بستهبندی ژنوم درون ویرون کمک کند. N همچنین یک آنتاگونیست اینترفرون است که سرکوبگر کدشده ویروسی یا از RNA اینترفرونها است که تکثیر و همانندسازی ویروسی را تسریع میبخشد. به احتمال زیاد، این پروتئینهای ساختاری در SARS-COV-2 به اتصال به رسپتورهای آنژیوناسین2 کمک میکنند و باعث تسهیل این اتصال میشوند. درنتیجه به ورود ویروس به سلول میزبان کمک میکنند و باعث بالا رفتن نرخ آلودگی و انتقال میشوند [۵].
جیهووا و همکاران در مطالعهای چنین عنوان کردند کهSARS-COV-2 نسبت به سارس تمایل بیشتری به آنژیوناسین2 دارد که این موضوع، بیانگر میزان بیشتر آلودگی و درصد بالاتر انتقالSARS-COV-2 نسبت به سارس است. همچنین پروتئینهای اسپایک٬SARS-COV-2 نسبت به سارس انرژی کمتری دارد که نشان از پایداری بیشتر SARS-COV-2 دارد و ممکن است توانایی زنده مانی در دماهای بالاتری نسبت به سارس داشته باشد. همچنین این امر میتواند بیانگر منشأ احتمالی SARS-COV-2 از خفاش باشد، زیرا خفاشها دمای بدن بالاتری نسبت به انسان دارند. علاوهبراین٬ RBD پروتئین SARS-COV-2 انعطافپذیری بسیار بالاتری بهویژه در نزدیکی سایت اتصال، نسبت به سارس دارد. این نشان میدهد که SARS-COV-2 برای اتصال به آنژیوناسین2 باید بر آنتروپی بالا غلبه کند [7]. بنابراین نسبت به سارس٬ در دمای بالا فعالیت کمتری دارد و با افزایش دما، بسیار کمتر از سارس باعث مرگومیر میشود. برخلاف سارس که پس از سال 2003 از بین رفت، ممکن است SARS-COV-2 بهعلت مطابقت آن با دمای بالاتر از دمای بدن انسان، توانایی زنده ماندن بالاتری در تابستان _که اکثر ویروس ها فعال نیستند_داشته باشد که این امر بهعلت انعطافپذیری پروتئین RBD ویروس SARS-COV-2 است. به همین دلیل SARS-COV-2 حتی در تابستان نیز عفونی است [7]. اما شواهد نشان میدهد که افزایش درجه حرارت بالاتر از 60 درجه باعث نابودی این ویروس میشود [۵، ۹].
بررسی مطالعات بیوانفورماتیکی
در اینجا مروری بر مطالعات بیوانفورماتیک در موردORFهای ویروس و منشأ احتمالی SARS-COV-2 می شود که آیا امکان آزمایشگاهی بودن این ویروس وجود دارد یا خیر و بدین منظور٬ سناریوهای مختلف ارائهشده مورد بحث قرار میگیرد. تجزیهوتحلیل بهدستآمده به لزوم مطالعات بیشتر درخصوص آزمایشگاهی بودن SARS-COV-2 اشاره میکند.
بررسی بیوانفورماتیکی ORF8 و سطح گلیکوپروتئین مهارکننده متابولیسم هِم توسط اتصال به پورفیرین
باتوجهبه محدودیتهای موجود در تجهیزات آزمایشگاهی، بسیاری از نقشهای پروتئینی SARS-COV-2 ازجمله ORF8 هنوز نامشخص است. بنابراین در صحنه همهگیری اخیر، پیشبینی نقش بیولوژیک پروتئینهای ویروسی از طریق بیوانفورماتیک از اهمیت علمی بالایی برخوردار است [۹]. ونژانگ لیو و همکاران نیز در مطالعهای بیوانفورماتیکی به پیشبینی نقش ORF8 این ویروس پرداختند [۱۰].
اکثر آزمایشهای برگرفته از بیماران دارای کووید-۱۹ نشان میدهد که تعداد نوتروفیلها در بیماران کاهش مییابد و مقادیر شاخص فرریتین سرم، سرعت رسوب اریتروسیت، پروتئین واکنشی C، آلبومین و دهیدروژن از بیماران بهطور معناداری افزایش مییابد [۴]. این، نشاندهنده کاهش میزان هموگلوبین در بیماران و افزایش میزان هِم است. با این اوصاف، بدن بسیاری از یونهای آهن مضر را انباشته خواهد کرد که باعث التهاب در بدن میشود و پروتئین واکنشی c و آلبومین را افزایش خواهد داد [۱۰]. سلولها به وجود التهاب واکنش نشان میدهند و حجم زیادی از فرریتین سرم را تولید میکنند تا یونهای آهن آزاد را جذب کنند و به این ترتیب آسیب ناشی از آن را کاهش دهند [۱۰].
هموگلوبین شامل 4 زیرواحد: 2α و 2β است و هر زیرواحد دارای یک قسمت هم متصل به آهن است. هِم یکی از اجزای مهم هموگلوبین و یک پورفیرین آهن است. ساختار بدون آهن آن، پورفیرین نامیده میشود. وقتی آهن دو ظرفیتی است، هموگلوبین ممکن است دیاکسیدکربن را آزاد کند و اتمهای اکسیژن را در سلولهای آلوِئولار جذب کند و باعث اکسیدهشدن آهن و ۳ ظرفیتیشدن آن شود [۱۱]. وقتی هموگلوبین در بدن از طریق جریان خون دردسترس سلولها قرار میگیرد، اتمهای اکسیژن را آزاد میکند. سپس دیاکسیدکربن را جذب میکند و آهن ۲ظرفیتی میشود. هیچ دارویی در حال حاضر برای درمان بیماری کووید-19 وجود ندارد.
بااینحال چند داروی قدیمی میتوانند برخی از توابع ویروس را مهار کنند [۱۱]. برای مثال٬ کلروکین فسفات، تأثیر مشخصی بر ذاتالریه کووید-۱۹ دارد. این دارو، یکی از داروهای ضد مالاریا است. اثرات درمانی کلروکین فسفات بر ذاتالریه کووید-19 حاکی از آن است که ذاتالریه این ویروس٬ ممکن است ارتباط نزدیکی با متابولیسم هموگلوبین غیرعادی در انسان داشته باشد [۱۱، ۱۲]. کلروکین، یک داروی رایج برای درمان پورفیریا است. بنابراین، این عقیده وجود دارد که ترکیب پروتئینهای ویروسی و پورفیرینها منجر به یک سری واکنشهای پاتولوژیکال انسان مانند کاهش هموگلوبین میشود [۱۲].
گروه ونژانگ لیو با استفاده از توالیهای پروتئینی موجود در پایگاه داده مرکز ملی اطلاعات زیستفناوری، تمام پروتئینهای ویروس کرونای جدید، پروتئینهای متصلشونده به هِم و هموکسیداز و توالی پروتئینی را برای تجزیهوتحلیل مورد استفاده قرار دادند. همچنین از تمامی پروتئینهای ویروس کرونای جدید چینی برای ساخت سازههای ۳بعدی با استفاده از مدلسازی همولوژی استفاده شده است [۱۰]. در مطالعه این گروه، مجموعهای از تحلیل بیوانفورماتیکی بر پایه توالی پروتئین زیستی منتشر شد که نتایج در تصویر شماره ۲ نشان داده شده است. در این تصویر، دامینهای حفاظتشده پروتئینهای ویروسی توسط سرور آنلاینMeMe تحلیل شده است. دامینهای حفاظتشده برای پیشبینی تفاوتهای عملکرد پروتئین ویروسی و انسانی مورد استفاده قرار گرفتند. همچنین ساختار ۳بعدی پروتئین ویروسی با استفاده از همولوژی مدلینگ سوئیس مدل ساخته شد [۱۰]. چون طول توالی بیش از 5000 نوکلئوتید بود، ابزار مدلینگ همولوژی Discovery-Studio 2016 انتخاب شد. در انتها با استفاده از تکنولوژی داکینگ مولکولی (LibDock tool) از Discovery-Studio 2016 اتصال گیرنده لیگاند پروتئین ویروسی با پروتئین هِم انسان شبیهسازی شد. روند این کار٬ مبتنی بر اصول تکاملی است. اگر چه مشخصه توالی بیولوژیکی شکلهای زیستی پیشرفته و ویروسی متفاوت است، اما مولکولهای با ساختار مشابه همیشه میتوانند به نقشهای زیستی مشابه دست یابند [۱۰].
MeMe یک وبسایت آنلاین است که بسیاری از ابزارهای پیشبینی و موتیفها را با هم ادغام میکند. موتیف یک دامین حفاظتشده از یک توالی کوچک در یک پروتئین است. مدلهای مبتنی بر موتیف میتوانند قابلیت اطمینان تحلیل فیلوژنیکی را ارزیابی کنند. پس از باز کردن ابزار آنلاین MeMe، توالیهای پروتئینی موردنظر٬ در یک فایل متنی ادغام میشوند و فرمت فایل بهصورت FASTA باقی میماند. سپس تعداد موتیفهایی که قرار است پیدا شوند، انتخاب میشوند و با کلیک بر دکمه go نتایج لازم به دست میآید [۱۰]. در پایان تحلیل دامینهای حفاظتشده پس از کلیک بر روی لینک نمایش داده میشود.
سوئیس مدل٬ یک سرور مدلسازی کاملاً خودکار برای ساخت مدلینگ همولوژی از ساختار پروتئین است. پس از ورود به نرمافزار و وارد کردن توالی روی Search Template کلیک میشود و پس از تکمیل جستوجو٬ الگو برای مدلسازی انتخاب میشود [۱۰، ۱۱]. سپس جستوجوی الگو با کلیک بر Build Model انجام میشود و یک مدل الگو بهطور خودکار انتخاب میشود. نرمافزار٬ چندین الگو را جستوجو میکند و مدلهای مختلف را ارائه میدهد. سپس مدل در قالب PDB دانلود میشود و در VMD شبیهسازی میشود. سپس داکینگ مولکولی طبق مراحل به پیش میٰرود. بدین گونه که ابتدا یک مدل لیگاند آمادهسازی میشود. یک فایل لیگاند از قبیل هِم باز میشود و روی Prepare Ligands در زیر منوی Dock Ligands از Receptor-Ligand interactions کلیک میشود تا یک مدل لیگاند هِم برای اتصالبه وجود آید.
در ابتدا اتم آهن از هِم حذف میشود. سپس با کلیک بر Prepare Ligands مدل پورفیرین تولید میشود. با استفاده از 0 XT باز شده باید Prepare Ligands را دوباره برای به دست آوردن مدل لیگاند کلروکین کلیک کرد [۱۰]. سپس برای آمادهسازی مدل گیرنده پروتئینی، باید فایلهای PDB بهدستآمده از همولوژی مدلینگ را باز کرد و روی Prepare Protein در زیر منوی Dock Ligands ازReceptor-Ligand interactions کلیک کرد تا یک مدل گیرنده پروتئین جهت داکینگ استفاده شود. سپس باید پارامترهای داکینگ را برای به دست آوردن داکینگ تنظیم کرد. برای این کار در ادامه generated protein receptor model انتخاب میشود و از منوی Define and Edit Binding Site در زیر منوی Receptor-Ligand interactions روی receptor Cavities کلیک میشود و با این کار٬ دایره قرمز روی نمودار مدل گیرنده پروتئین ظاهر میشود. پس از راستکلیک کردن٬ میتوان شعاع دایره قرمز را تغییر داد.
در منوی Receptor-Ligand interactions و زیرمنوی Dock Ligands، گزینه Dock Ligand (LibDock)، انتخاب میشود و از قسمت pop-up box تمام لیگاندهای جدید پیشبینیشده بهعنوان گیرنده، newly established receptor model-ALL انتخاب میشود و سایتهای دایره بهعنوان مختصات دایره انتخاب میشو. درنهایت برای شروع اتصال روی Run کلیک میشود.سپس پس از تکمیل اتصال٬ با محاسبه انرژی پیوندی و انتخاب بهترین انرژی اتصال بهعنوان لیگاند٬، مکانهای زیادی به نمایش گذاشته میشود [۱۰] که بدین منظور با باز کردن Docked view و کلیک روی دکمه Caculate Binding Energies از Dock Ligands منوی Receptor-Ligand interactions درpop-up box گیرنده بهعنوان مقدار پیشفرض و لیگاند بهعنوان docked model-ALL انتخاب میشود٬ سپس محاسبه انرژی اتصال آغاز میشود.
درنهایت با مقایسه انرژی اتصال بزرگترین انرژی اتصال انتخاب میشود. هرچه پایداری کمپلکس بیشتر باشد، انرژی اتصال بیشتر خواهد بود. سپس باید بخش مشترک را نمایش داد. پس از تنظیم استایل نمایشی از منطقه اتصال، روی دکمه 2D Map از زیرمنویView Interaction موجود در Receptor-Ligand interactions در pop up the view of binding section کلیک میشود. این شکل بهعنوان فایل تصویری قابل ذخیره شدن است [۱۰]، اما بهعنوان نتایج حاصل از این کار میتوان چنین گفت که در انسانها هموگلوبین را میتوان به 2 بخش گلوبین و هِم تجزیه کرد. قسمت هِم از پورفیرین و یون آهن تشکیل میشود و یون آهن در میان پورفیرین قرار دارد [۱۲، ۱۳]. همچنین در آب٬ نامحلول است و با پروتئینهای اتصالی هم ترکیب میشود تا یک کمپلکس را تشکیل دهد و به کبد منتقل شود.
پورفیرین بهصورت بیلیروبین تغییر شکل مییابد و از مجرای صفراوی خارج میشود و آهن موجود در مولکول میتواند توسط بدن مورد مصرف قرار گیرد [۱۳]. اگر پروتئینهای ویروس بتوانند به قسمت پورفیرین از هِم متصل شوند باید توانایی اتصال مشابهی را برای Heme-binding protein انسانی داشته باشند، یعنی پروتئینهای ویروسی و Heme-binding protein انسانی باید دامینهای حفاظتشده مشابهی داشته باشند [۱۳].
روشهای گفتهشده بالا برای پروتئینهای مشابه ویروسی انجام شد و 3 پروتئین ویروسی گلیکوپروتئین سطحی، پروتئین پوششی و نوکلئوکپسید فسفوپروتئین مورد استفاده قرار گرفت و دامینهای 3 پروتئین p کوچکی داشتند و متفاوت بودند که نشان میدهد توانایی پروتئین ساختاری برای اتصال به پورفیرین کمی متفاوت است. علاوه بر این٬ گلیگوپروتئین غشایی نمیتواند به پورفیرین متصل شود. سپس سرور آنلاین سوئیس مدل سطح گلیکوپروتئین ها را برای ایجاد یک ساختار 3 بعدی مدلسازی کرد و دو نوع فایل بر اساس مدل پروتئین اسپایک و نمونه E2 انتخاب شد [۱۰]. سپس فایل ساختاری 3D از بانک اطلاعات پروتئین دانلود شد. مدل ساختاری، حاکی از عدم توانایی گلیکوپروتیئن سطحی به قسمت پورفیرین از هِم بود٬ اما همین اقدامات برای پروتئین پوششی، حاکی از توانایی اتصال پروتئین پوششی به پورفیرین است. تصویر شماره 3، نمای 2 بعدی مقطع اتصال است که در آن تعداد 18 آمینواسید از پروتئینهای پوششی ویروس با پورفیرین انسانی تعامل دارد. اما همین اقدامات برای نوکلئوکپسید فسفوپروتئین نیز انجام شد که نتایج حاکی از توانایی اتصال نوکلئوکپسید فسفوپروتئین به پورفیرین با بالاترین انرژی اتصال ممکن یعنی 5206.53 kcal/mol است.
تصویر شماره 4 نمای داکینگ را نشان میدهد که نمای 2بعدی مقطع اتصال از نوکلئوکپسید فسفوپروتئین ویروس با پورفیرین انسانی است [۱۰-۱۳].
اگر پروتئینهای غیرساختاری بتوانند به قسمت پورفیرین هِم متصل شوند، باید توانایی اتصال مشابهی با پروتئین اتصالی هِم انسانی داشته باشند. سپس سرور آنلاین MEME برای جستوجوی دامینهای حفاظتشده بین پروتئینهای غیرساختاری ویروس و Heme binding protein انسانی مورد استفاده قرار میگیرد. برای مطالعه خصوصیات اتصال پروتئین ORF8 به هِم همان مراحل تجزیهوتحلیل همانند پروتئینهای ساختاری انجام میشود. چندین نوع شکل ساختاری متصل کردن پروتئین ORF8 به پورفیرین نشان از داشتن بالاترین سطح انرژی یعنی 25/859 kcal/mol بود. ضمن اینکه نتایج داکینگ مولکولی نیز تاGییدکننده اتصال مناسب پروتئین ORF8 به پورفیرین و اشتراک آنها در 18 آمینو اسید دارد. تصویر شماره 3 نمای 2بعدی مقطع اتصال از پروتئینهای غیرساختاری ORF8 ویروس با پورفیرین انسانی است [۱۰].
ازآنجاکه همهگیری SARS-COV-2 بهعنوان یک معضل جهانی مطرح است، روشهای پیشبینی برای جستوجوی دامینهای حفاظتشده بسیار حائز اهمیت است. دستیابی به ساختار مولکولهای پروتئین ازجمله ORF8 و نوکلئوکپسید فسفوپروتئینها (پروتئینهای سطحی) با استفاده از روشهای مدلسازی بیوانفورماتیکی و تجزیهوتحلیل امکان و نحوه اتصال آن به پورفرین انسانی از قسمت هِم در هموگلوبین بسیار حائز اهمیت است.
مطالعات ونژانگ لیو و همکاران نشان داد ORF8 و گلیکوپروتئینهای سطحی ویروس میتوانند بهصورت یک کمپلکس با پورفرین انسانی ترکیب شوند [۱۰]. همچنین این گروه نشان دادند که بهترتیب پروتئینهای ORF10، ORF1ab، ORF3a به قسمت هِم انسانی حمله میکنند تا آهن را برای تشکیل پورفیرین جدا کنند. پس از جدا شدن آهن از هِم و تشکیل پورفیرین دو پروتئین ORF8 و گلیکوپروتئینهای سطحی ویروس SARS-COV-2 با پورفیرین یک کمپلکس قوی تشکیل میدهند که مانع از باند شدن آهن به پورفیرین و در نهایت هموگلوبین میشوند که این امر باعث کاهش شدید اکسیژنرسانی به اندامها میشود [۱۰]. ازآنجاکه 2 کمپلکس پورفیرین تولیدشده در بدن انسان، مسیرهای آنابولیک هِم را مهار میکند، میتواند باعث ایجاد طیف وسیعی از بیماریهای عفونی شود. با این یافته میتوان چنین تصور کرد که کلوروکین میتواند از پروتئینهای ORF10، ORF1ab، ORF3a برای حمله به هِم جلوگیری کند و مانع از ایجاد پورفیرین شود. به همین دلیل بین دو پروتئین ORF8 و گلیکوپروتئین سطحی کروناویروس، با پورفرین بدن شخص بیمار، کمپلکس تشکیل می شود [۱۰، ۱۱].
بررسی بیوانفورماتیکی ساختار ژنوم SARS-COV-2
مقایسه ژنوم آلفا و بتاکروناویروس با هم دو نکته قابل توجه در مورد SARS-COV-2 را مشخص میکند: در ابتدا بر اساس مطالعات ساختاری که روی SARS-COV-2 انجام شده است، به نظر میرسد که این ویروس برای اتصال به گیرنده سلول انسانی آنژیوناسین2 بهینه شده است. درثانی، پروتئین S از SARS-COV-2 یک محل کلیواژ عملکردی polybasic به نام فارین در مرز S1-S2 دارد که از طریق وارد شدن 12 نوکلئوتید به وجود آمده است که منجر به ایجاد سه سایت پیوندی O-گلیکان در اطراف این سایت شده است [۱۴].
موتاسیون در دامین اتصال به گیرنده SARS-COV-2
دامین اتصال به گیرنده در پروتئین S یک بخش بسیار متغیر در ژنوم SARS-COV-2 است. یافتهها نشان میدهند که 6 آمینواسید RBD برای اتصال به گیرندههای آنژیوناسین2 و تأمین دامنه میزبان از ویروسهای شبهسارس بسیار حیاتی هستند. مطابق با تحلیلهای انجامشده که در مجله پزشکی طبیعی چاپ شد، این 6 آمینواسید در سارس بهترتیب Y442، L472، N479، D480، T487 و Y4911 است که دقیقاً با آمینواسیدهای L455، F486، Q493، S494، N501 و Y505 از SARS-COV-2 ازنظر محل قرارگیری مطابقت دارد [۱۴]. اما 5 مورد از این 6 رزیدو بین SARS-COV-2 و سارس متفاوت است (تصویر شماره 5).
براساس مطالعات ساختاری و آزمایشات بیوشیمیایی SARS-COV-2 به نظر میرسد که قسمت RBD از SARS-COV-2 توانایی اتصال با همولوژی بالا به گیرندههای آنژیوناسین2 از انسان، سمور، گربه و دیگر گونهها دارد. در حالی که تجزیهوتحلیل بالا نشان میدهد SARS-COV-2 ممکن است میل ترکیبی بالایی با آنژیوناسین2 انسانی داشته باشد، اما تجزیهوتحلیل محاسباتی پیشبینی میکند که این برهمکنشها ایدهآل نیست و توالی RBD با توالی سارس که یک توالی بهینه است، تفاوت دارد. در اقع این توالی در سارس نسبت به SARS-COV-2 تعاملات بهینهتری دارد [۱۴]. بنابراین کارایی بالای اتصال پروتئینS از SARS-COV-2 به آنژیوناسین2 انسانی، به احتمال زیاد٬ ناشی از انتخاب طبیعی روی یک انسان یا موجوداتی است که گیرنده آنژیوناسین2 شبهانسانی که اجازه اتصال بهینه را به ویروس میدهند. این٬ یک مدرک کاملاً معتبر بر آزمایشگاهی نبودن ویروس SARS-COV-2 است، اما تأئید آن نیازمند مطالعات موشکافانه بیشتر است [۱۴].
محل کلیواژ polybasic فارین و پیوند O-گلیکان
دومین قسمت قابلتوجه از SARS-COV-2 یک محل کلیواژ polybasic (RRAR) در نقطه اتصال S1-S2 یعنی 2 زیر واحد S است. این کار٬ امکان تقسیم مؤثر توسط فارین و پروتئازهای دیگر را فراهم میکند و در تعیین دامنه آلودگی ویروس و دامنه میزبانی نقش دارد [۱۵]. علاوهبراین یک پرولین نیز در این محل از SARS-COV-2 وجود دارد. بنابراین توالی واردشده بهصورت PRRA است. تغییرات ایجادشده توسط پرولین به نظر میرسد که منجر به اضافه شدن پیوند -Oگلیکانها به S673، T678 و S686 شود که محل تقسیم و کلیواژ و انحصاری SARS-COV-2 است [۱۶]. محلهای کلیواژ polybasic در اصل در دودمان B از بتاکروناویروسها وجود ندارند. اگرچه در دیگر کروناویروسهای انسانی نظیر HKU1 از دودمان A این سایتها و پیوند O-گلیکان وجود دارد [۱۵].
باتوجهبه سطح تنوع ژنتیکی بالا در اسپایک، این احتمال وجود دارد که ویروسSARS-COV-2 با جایگاههای کاتالیتیک کامل در گونههای دیگر نیز کشف شود. همچنین عملکرد محل کلیواژ polybasic که برای تشخیص انتقال و بیماریزایی ویروس در مدلهای حیوانی تعیینکننده است، مشخص نیست [۱۴]. آزمایشهایی با ویروس سارس نشان داد که الحاق یک محل کلیواژ فارین در نقطه اتصال S1-S2، فیوژ سلولی را بدون تأثیر بر ورود ویروس افزایش میدهد. در ویروس آنفولانزای مرغی، تکثیر و انتقال سریع در جمعیتهای بهشدت متراکم مرغ برای به دست آوردن محل کلیواژ polybasic در پروتئین هماگلوتینین انتخاب میشود [۱۷، ۱۸] که عملکردی مشابه با پروتئین S ویروس کرونا دارد.
دریافت محل کلیواژ polybasic به روش الحاق یا رونویسی، ویروس آنفولانزای مرغی را بهشدت بیماریزا میکند [۱۸]. این موضوع همچنین پس از ورود محل کلیواژ polybasic توسط پروتئین هماگلوتینین در محیط کشت سلولی و سپس در حیوانات مشاهده شده است. عملکرد پیوند O- گلیکان مشخص نیست، اما این پیوندها میتوانند یک دامین شبهموسین ایجاد کنند که اپیتوپها را محافظت میکند و یا میتوانند کلیدی برای رزیدوهای پروتئین S از SARS-COV-2 باشد [۱۹]. بسیاری از ویروسها از دامینهای شبهموسین بهعنوان سپرهای دفاعی استفاده میکنند. اگرچه به نظر میرسد که پیوند O-گلیکوزیلاسیون قوی است، اما مطالعات تجربی برای تعیین اینکه آیا این سایت ها در SARS-COV-2 استفاده میشود، مورد نیاز است. در رابطه با امکان دستکاری آزمایشگاهیSARS-COV-2 همانطور که در سطور قبل اشاره شد قسمت RBD از SARS-COV-2 متفاوتتراز آنچه که از قبل بوده، برای اتصال به آنژیوناسین2 انسانی بهینه شده است [۲۰]. علاوهبراین اگر دستکاری ژنتیکی انجام شده باشد، یکی از چند سیستم ژنتیکی موجود برای بتاکروناویروسها احتمالاً باید مورد استفاده قرار گرفته باشد. با این حال دادههای ژنتیکی نشان میدهد که SARS-COV-2 از هیچیک از ویروسهای استفادهشده قبلی مشتق نشده است [۲۱]، اما دو سناریو برای احتمال وجود SARS-COV-2 وجود دارد:
مورد اول: انتخاب طبیعی در یک میزبان حیوانی قبل از انتقال Zoonotic؛
مورد دوم: انتخاب طبیعی در انسانها پس از انتقال Zoonotic است [۱۴]. نتایج این مطالعات به لزوم مطالعات بیوانفورماتیکی بیشتر در خصوص سناریوی آزمایشگاهی بودن ویروس دلالت دارد.
نتیجهگیری
در دو دهه اخیر دو ویروس کرونا از جنس بتاکروناویروسها به نامهای سارس و مرس باعث وحشت عمومی مداوم شدند و به مهمترین وقایع بهداشت عمومی تبدیل شدند تا اینکه ویروس کرونای جدید 2019 با پنومونی وسیعی از ووهان چین گسترش یافت [۵۴]. در این مطالعه به جمعآوری مجموعه اطلاعات حاصل از پژوهشهای انجام شده در خصوص کووید-19 پرداخته شد. بررسی ویژگیهای ویروس نشان داد همانندسازی SARS-COV-2 توسط اتصال پروتئین S به گیرنده سطح سلول شروع میشود. SARS-COV-2 آنزیم سطح سلولی آنژیوناسین۲ را بهعنوان گیرنده اختصاصی جهت ورود به سلول استفاده میکند. ازآنجاییکه میزان بیان این آنزیم در جوامع آسیایشرقی بیشتر از سایر جوامع است، لذا شیوع آن در این جوامع بیشتر از سایرین بوده است [۵۵]. باتوجهبه خصوصیات منحصربهفرد ماده ژنتیکی کروناویروسها و همچنین مکانیسم تکثیر آنها، ظهور ویروس کرونای جدید برای ویروسشناسان دور از انتظار نبوده است. ازجمله این ویژگیهای خاص میتوان به میزان بالای نوترکیبی ژنتیکی این ویروس اشاره کرد که هرچند سال یکبار منجر به پدیدار شدن سویههای جدید و ناشناخته میشود. دلیل این مقدار بالای نوترکیبی را میتوان به طول بسیار بلند ماده ژنتیکی این ویروس در مقایسه با سایر ویروسها، تکثیر پیچیده و میزان بالای خطای آنزیم همانندساز این ویروس و درنهایت، به دامنه میزبانی وسیع این ویروسها در انسانها و حیوانات مختلف نسبت داد [۵۶]. از زمان شیوع بیماری همهگیر کووید-19، استفاده از بیوانفورماتیک برای تجزیهوتحلیل نقش پروتئینهای کروناویروس جدید (مانند ORF8 و گلیکوپروتئین سطحی) اهمیت علمی بالای خود را یادآور شد [۵۷].
در این مطالعه، چگونگی بررسی بیوافورماتیکی و روشهای پیشبینی دامنه برای جستوجوی دامنههای حفاظتشده بیان شدند. مطالعات بیوانفورماتیکی مختلف درخصوص ساختار مولکولهای پروتئینی مانند ORF8 و گلیکوپروتئین سطحی با استفاده از روشهای مدلسازی همولوژی و تکنولوژی داکینگ مولکولی برای تجزیهوتحلیل بخش اتصال پروتئین ویروسی به پورفرین بررسی شد [۵۸]. جمعبندی نتایج این مطالعات نشان میدهد که ORF8 و گلیکوپروتئین سطحی میتوانند بهترتیب با هم ترکیب شوند تا یک کمپلکس را شکل دهند. همچنین، ORFهای دیگر میتوانند روی زنجیره ۱- بتا در هموگلوبین حمله کنند تا آهن را تجزیه کرده و porphyrin را شکل دهند. این حمله منجر به کاهش میزان هموگلوبین در دسترس برای حمل اکسیژن و دیاکسید کربن میشود.
سلولهای ریه بهعلت ناتوانی در مبادله کربن دیاکسید و اکسیژن به شدت التهاب دارند که در نهایت منجر به تصاویر ground مانند در سیتیاسکن میشوند. با این یافتههای بیوانفورماتیکی و تجزیهوتحلیل بیشتر میتوان فهمید که کلروکین ممکن است از حمله به هموگلوبین برای تشکیل porphyrin جلوگیری کند و مانع از اتصال ORF8 و سطحی به porphyrins تا حد معینی شود و بهطور مؤثر علایم حاد تنفسی را تسکین دهد. در حال حاضر جنگ سختی بین ماده ژنتیکی و هوش ویروسی با عقل و هوش انسانی در جریان است و برای پیروزی در این نبرد، علاوه بر شناخت کافی از این ویروس مهاجم و مسری، نیاز به اتخاذ تصمیمات کنترلی صحیح و بموقع از سوی کشورها و اهتمام به انجام اقدامات محافظتی و بهداشتی شخصی است. باتوجهبه همهگیری کنونی، اعتقاد بر این است که نتایج تحقیقات بیوانفورماتیکی از ارزش بالایی در پیشگیری از شیوع ذاتالریه ویروس کرونا، توسعه داروها و واکسنها و درمان بالینی برخوردارند.
ملاحظات اخلاقی
پیروی از اصول اخلاق پژوهش
این پژوهش مطابق با اصول اخلاق پژوهشی مصوب دانشگاه آزاد اسلامی واحد شهرکرد، نگارش شده است.
حامی مالی
این تحقیق هیچگونه کمک مالی خاصی از سازمان های تأمین مالی در بخشهای عمومی، تجاری یا غیرانتفاعی دریافت نکرد.
مشارکت نویسندگان
مفهومسازی: توحید پیری قراقیه؛ تحقیق و بررسی: شیدا بیرانوند؛ ویراستاری و نهاییسازی نوشته: سامه حاجی محمدی و امیرحسین قدیری.
تعارض منافع
بنابر اظهار نویسندگان، در این مقاله هیچگونه تعارض منافعی را ندارند.
تشکر و قدردانی
نویسندگان از کمک و حمایت کارکنان مرکز تحقیقات بیوتکنولوژی دانشگاه آزاد اسلامی واحد شهرکرد در ایران تشکر میکنند.
References