پیشگویی جهش‌های پایدارکننده بالقوه پروتئین اینترفرون بتا a1 (سینووکس) به‌وسیله شبیه‌سازی دینامیک مولکولی

داروی ایرانی سینووکس (اینترفرون بتا) یک داروی مهم در جلوگیری از پیشرفت بیماری‌ام. اس است و افزایش پایداری آن اهمیت زیادی دارد. در این تحقیق با طراحی پنج جهش برای این پروتئین و انجام شبیه‌سازی دینامیک مولکولی و داکینگ پایدارکننده‌ترین جهش، مشخص شد. برای انجام این کار در ابتدا ساختار طبیعی اینترفرون بتا به درگاه اینترنتی پاپ میوزیک داده شد تا جهش‌های پایدارکننده پیشگویی شوند. سپس مدل چهار جهش مناسب‌تر پیشگویی شده توسط این سرور و یک جهش ترکیبی (جهش ترکیبی 2 و 3 یا جهش 5) از این پروتئین تهیه شد. سپس مدل هرکدام از این جهش‌ها و ساختار طبیعی به‌طور جداگانه در دماهای 300 و 350 و 400 و 450 و 500 کلوین به مدت 20 نانوثانیه به‌وسیله برنامه گرومکس 5 شبیه‌سازی دینامیک مولکولی شدند. (درمجموع 600 نانوثانیه). سپس پارامترهای ساختاری برای هرکدام از این ساختارها در هر دما در 10 نانوثانیه آخر شبیه‌سازی‌ها محاسبه شد و نمودار دما در برابر میانگین پارامترها رسم شدند و شیب تغییرات ساختاری در اثر افزایش دما برای همه جهش‌ها و ساختار طبیعی به دست آمد تا ساختارهای پایدارتر مشخص شوند. نتایج به‌دست‌آمده نشان داد که جهش شماره سه (g127e) در بین پنج جهش پایداری بیشتری در مقایسه با ساختار اینترفرون طبیعی دارد. بعلاوه قدرت اتصال جهش شماره سه به گیرنده اینترفرون نسبت به حالت وحشی و سایر جهش‌ها بیشتر بود. به نظر می‌رسد که افزایش هیدروفیلیستی اطراف باقی‌مانده 127 در جهش 3 در سطح پروتئین منجر به افزایش حلالیت و پایداری بیشتر آن در خارج از بدن شده است.

prediction of potential stabilizer mutation of interferon beta 1 a (cinovex) via molecular dynamics simulation

iranian drug cinnovex (interferon beta) is an important drug for preventing the progress of m.s. disease and increase of its stability is very important. in this research, we designed five mutations in the structure of interferon and performed molecular dynamics simulation and docking and the best stabilizer mutation was determined. at first, the wild type structure of beta interferon was subjected to pop music server for prediction of stabilizer mutations. then the first four proper score mutations were selected and mutated structures and a mix mutation (mut 2+ mut 3) (mut 5) were made. then, molecular dynamics simulation performed via gromacs 5 package on wild and mutated structures in five temperatures (300, 350, 400, 450, 500 k) separately for 20 ns (in total 600 ns). the simulation parameters were calculated in each temperature for all structures during the last 10 ns of simulations. the plots of temperature against averages of parameters were drawn and the regression slope of parameters was calculated to find the more stable structure. the results revealed that mutation 3 (g127e) is the best mutation compare to wild interferon. also, the binding power of mutation 3 with its receptor is more than wild type and other mutations. it seems that the increase of hydrophilicity of protein around position 127 in mutation 3 on interferon surface leads to increase the stability and solubility of it ex vivo.