شبیه سازی دینامیک مولکولی برهمکنش ساخارین با پروتئین p53 انسانی

زمینه و هدف: نقش شیرین کننده‌های مصنوعی در وقوع خطر سرطان در طی چند دهه گذشته به طور گسترده‌ای مورد بحث قرار گرفته است. لذا هدف از این مطالعه تعیین شبیه‌سازی دینامیک و مولکولی بر همکنش ساخارین(sa) با پروتئینp53 انسانی بود.روش بررسی: این مطالعه بیوانفورماتیک در سال 1402 انجام شد. بر همکنش sa و سدیم ساخارین (ssa) با پروموتر(راه انداز) ژن p53 انسانی(pp53g) قبلا در سال 1399 در دو بخش نظری و تجربی چاپ شده است، اما در مطالعه حاضر، قابلیت اتصال و جایگاه اتصال لیگاند sa به عنوان یک شیرین کننده‌ مصنوعی با پروتئین p53 انسانی(رسپتور) به عنوان سرکوبگر تومور به صورت نظری انجام شد. رزیجوهای اسیدآمینه‌های درگیر در برهمکنش، انرژی آزاد اتصال و ثابت اتصال تعیین شد. برای محاسبات برهمکنش مولکولی، از داکینگ مولکولی استفاده شد. اطلاعات دقیق‌تر در خصوص شیوه اتصال کمپلکس لیگاندرسپتور با شبیه‌سازی دینامیک مولکولی md)) به دست آمد. فایل ساختار و توپولوژی برای پروتئینp53 انسانی استخراج شده از پایگاه داده‌ پروتئین بر اساس میدان نیروی amber 99 با برنامه gromacs 5.3.1 ساخته ‌شد. برنامه acpype/antechamber با gaff برای ایجاد فایل ساختار و توپولوژی لیگاند در شبیه سازی md استفاده ‌شد. این میدان نیرو با میدان نیروی  amber 99سازگار است. زمان شبیه سازی در حلال صریح 50 نانوثانیه برای کمپلکس پروتئین p53 با sa بود. داده‌های جمع‌آوری شده با استفاده از نرم‌افزارهای مختلف و مقایسه با نتایج مقالات مرتبط تجزیه و تحلیل شدند.یافته‌ها: نتایج داکینگ مولکولی نشان داد که ترکیب sa به پروتئین p53 انسانی با انرژی اتصالی(پیوندی)، 4.55 -کیلوکالری بر مول و ثابت اتصال 462.18 میکرومولار اتصال یافت. یک پیوند پیوند هیدروژنی بین sa با اسیدآمینه leu137 تشکیل شد. تغییرات کنفورماسیونی حاصل از شبیه‌سازی md برای کمپلکس لیگاند پروتئین نشان داد که sa می‌تواند از طریق دو پیوند هیدروژنی به his179 و arg196 به عنوان اسید آمینه‌های کلیدی پروتئین p53 در ناحیه اتصال به dna متصل شود. sa هم‌چنین می‌تواند در مجاورت اسید آمینه‌های اسید آمینه‌های leu137، ala138،his179 ، asp184 و met237 از طریق پیوندهای هیدروفوب قرار گیرد. مقادیر نمودارهای انحراف جذر میانگین مربعات(rmsd)، نوسان جذر میانگین مربعات(rmsf)، شعاع چرخش(rg)، برای پروتئین p53 آزاد و در حضور لیگاند sa نشان از اتصال پایدار sa به پروتئین p53 دارد.نتیجه‌گیری: این مطالعه می‌تواند اطلاعات ارزشمندی را در مورد مکانیسم اتصال sa و پروتئین p53 انسانی به عنوان یک درشت مولکول در سطح مولکولی با جزییات ریز اتمی ارایه کند. نتایج این مطالعه می‌تواند در تعیین احتمال خطر بالقوه سرطان‌زایی این شیرین کننده با توجه به مصرف زیاد آن و طراحی و سنتز شیرین کننده‌های مصنوعی جدیدتر و ایمن‌تر مفید باشد. 

molecular dynamics simulation of the interaction of saccharin with human p53 protein

background & aim: the role of artificial sweeteners in the occurrence of cancer risk has been widely discussed during the last few decades. therefore, the aim of the present study was to determine the dynamic and molecular simulation of the interaction of saccharin(sa) with human p53 protein.methods: the present bioinformatics study was conducted in 2023. the interaction of sa and sodium saccharin (ssa) with the human p53 gene promoter (pp53g) has already been published in 2019 in two theoretical and experimental sections. but in the present study, the binding ability and the binding site of sa ligand as a synthetic sweetener with human p53 protein (receptor) as a tumor suppressor was theoretically performed. the amino acid residues involved in the interaction, energy free binding and binding constant were determined. molecular docking was used for molecular interaction calculations. more detailed information about the binding method of the ligand-receptor complex was obtained by molecular dynamics (md) simulation. the structure and topology file for the human p53 protein extracted from the protein database was created based on the amber 99 force field with the gromacs 5.3.1 program. acpype/antechamber program with general amber force field (gaff) was used to create structure file and ligand topology in md simulation. this force field was compatible with the amber 99 force field. the simulation time in explicit solvent was 50 ns for sa-p53 protein complex. the collected data were analyzed using different software and compared with the results of related articles.results: the results of molecular docking indicated that the sa compound was bound to human p53 protein with a binding energy of -4.55 kcal/mol and a binding constant of 462.18 μm. a hydrogen bond was formed between sa and amino acid leu137. the conformational changes resulting from md simulation for the ligand-protein complex showed that sa can bind to arg196 and his179 as key amino acids of p53 protein in the dna binding region through two hydrogen bonds. sa can also be placed in the adjacent of amino acids leu137, ala138, his179, asp184 and met237 through hydrophobic bonds. the values of the plots of root mean square deviation (rmsd), root mean square fluctuation (rmsf), radius of gyration (rg) for free p53 protein and in the presence of sa ligand show the stable binding of sa to p53 protein.conclusion: the present study could provide valuable information about the binding mechanism of sa to human p53 protein as a macromolecule at the molecular level with subatomic details. the results can be useful in determining the potential carcinogenic risk of this sweetener due to its high consumption and the design and synthesis of newer and safer artificial sweeteners.