In Silico Analysis of Pathogenic Missense Single Nucleotide Polymorphisms in the Bovine IGF-1 Gene

تخيل أن جينوم البقرة عبارة عن مكتبة ضخمة وقديمة تحتوي على كتيبات التعليمات لبناء وتشغيل كائن حي. داخل هذه المكتبة، تمثل الحروف A وC وG وT الأبجدية، وعندما تترتب بشكل صحيح، فإنها تكتب وصفات البروتينات — تلك الآلات الجزيئية الصغيرة التي تبقي الحيوان حياً، وتنمو، وتنتج الحليب. ولكن في بعض الأحيان، يتم استبدال حرف واحد في الوصفة بآخر، تماماً مثل خطأ مطبعي في كتاب طهي. في عالم الوراثة، يسمى هذا "تعدد أشكال النوكليوتيدات المفردة" (Single Nucleotide Polymorphism) أو SNP. معظم هذه الأخطاء غير ضارة، مثل تغيير كلمة "ملح" إلى "ملح" (بخطأ إملائي) في وصفة حساء؛ حيث يظل طعم الحساء جيداً. ومع ذلك، فإن بعض الأخطاء تكون كارثية، حيث تحول "السكر" إلى "ملح" وتفسد الطبق.

واحدة من أهم الوصفات في مكتبة البقرة هي لوقود هرمون يسمى IGF-1. فكر في IGF-1 كأنه رئيس عمال في موقع بناء، يصرخ بالأوامر للعمال لبناء العضلات، وتنمية العظام، وتطوير الغدد الثديية اللازمة لإنتاج الحليب. إذا كانت تعليمات رئيس العمال مشوشة بسبب خطأ مطبعي في الوصفة، فقد يتوقف موقع البناء، أو قد يكون المبنى مائلاً. بالنسبة لمربي الماشية، يعني وجود رئيس عمال معطل نقصاً في الحليب وعدد العجول، وهو أمر سيء للاقتصاد. لطالما عرف العلماء وجود هذه الأخطاء المطبعية الجينية، لكن العثور على الأخطاء المحددة التي تعطل الآلة بالفعل من بين ملايين الحروف غير الضارة يشبه البحث عن ترس واحد مكسور في ساعة ضخمة. وهنا تأتي قوة التحليل "في السيليكو" (in silico) أو التحليل الحاسوبي. فبدلاً من اختبار الأبقار في الحظيرة، يستخدم الباحثون أجهزة كمبيوتر فائقة لمحاكاة آثار هذه الأخطاء المطبعية، حيث يعملون كالمحققين الرقميين الذين يمكنهم التنبؤ بالأخطاء الجينية التي ستسبب أكبر قدر من المشاكل قبل حدوثها في الواقع.

في هذه القصة التحقيقية الرقمية، وضع فريق من الباحثين من جامعة لورستان في إيران خطة للعثور على أسوأ الأخطاء المطبعية في جين IGF-1 البقري. بدأوا بجمع قائمة تضم 66 نوعاً مختلفاً من "النيوكليوتيدات المفردة الخاطئة في المعنى" (missense SNPs) — وهي أخطاء مطبعية محددة تغير بالفعل حمضاً أمينياً واحداً (لبنة بناء البروتين) إلى حمض آخر مختلف، مما قد يغير شكل "رئيس العمال" IGF-1. وباستخدام مجموعة من الأدوات الحاسوبية المتطورة، عملوا كمرشح (فلتر)، يغربلون الضجيج للعثور على إشارات الخطر.

أولاً، مرروا القائمة عبر أداة تسمى PredictSNP، والتي تعمل كلوحة من الحكام الخبراء. ومن بين 66 مشتبهاً به، تم تحديد 18 منها على أنها "ضارة" (deleterious)، مما يعني أنها من المرجح أن تسبب ضرراً. ثم ركز الباحثون على هذه الـ 18 ليروا مدى توافقها مع التاريخ التطوري للبروتين. وباست استخدام أداة تسمى ConSurf، اكتشفوا أن هذه المواضع السيئة تقع في مناطق "محفوظة للغاية" (highly conserved) — وهي أجزاء من البروتين ظلت كما هي تماماً لملايين السنين لأنها حاسمة للبقاء. الأمر يشبه اكتشاف أن الخطأ المطبعي حدث في عجلة القيادة للسيارة؛ فأنت تعلم أن هذا مكان سيء للعبث به.

بعد ذلك، قام الفريق بمحاكاة كيفية تأثير هذه الطفرات على الاستقرار الفيزيائي للبروتين. استخدموا موازين رقمية (I-Mutant2.0 و MUpro) لمعرفة ما إذا كانت الطفرات تجعل البروتين مهتزاً أو عرضة للتفكك. وأظهرت النت النتائج أن معظم الطفرات الـ 18 جعلت البروتين أقل استقراراً، مثل استبدال عارضة فولاذية متينة بعارضة بلاستيكية واهنة. كما تحققوا من مرونة البروتين باستخدام أداة تسمى MEDUSA، ووجدوا أن هذه الطفرات عالقة في أجزاء صلبة من الهيكل، مما يجعل من الصعب أكثر على البروتين القيام بعمله.

تقلص التحقيق أكثر باستخدام MutPred2، وهي أداة تتنبأ بكيفية تعطيل الطفرة للآليات الجزيئية للبروتين. حددت هذه الخطوة سبعة "أشرار خارقين" كانوا خطيرين بشكل خاص. وأخيراً، استخدم الباحثون Project HOPE للنظر في الشكل ثلاثي الأبعاد للبروتين. قارنت هذه الأداة الحجم، والشحنة، والقدرة على كره الماء (hydrophobicity) للبروتينات المتحولة مقابل البروتينات الطبيعية. ووجدوا أن الطفرة L105R كانت الأكثر تدميراً على الإطلاق. في هذه الطفرة، تم استبدال حمض أميني متعادل ومتوسط الحجم (لوسين) بآخر كبير وذو شحنة موجبة (أرجينين). وقد حاكت الأبحاث أن هذا التغيير سيؤدي إلى حدوث "نتوءات" في بنية البروتين لأن القطعة الجديدة كبيرة جداً، كما أن الشحنة الموجبة ستتنافر مع أجزاء أخرى من الجزيء، مما يؤدي إلى تعطيل التفاعلات الكارهة للماء الدقيقة التي تربط البروتين معاً.

ولرؤية التأثير الواقعي لطفرة L105R، قام الفريق بإجراء محاكاة "الارتباط البروتيني-البروتيني" (protein-protein docking). لقد قاموا بنمذجة كيف يقوم بروتين IGF-1 الطبيعي والنسخة المتحورة منه بـ "المصافحة" مع شريكهما، وهو مستقبل IGF-1. في المحاكاة، شكل البروتين الطبيعي مصافحة قوية مكونة من 25 رابطة هيدروجينية. ومع ذلك، استطاع بروتين L105R المتحور تحقيق 21 رابطة هيدروجينية فقط. ويشير هذا النقص في الروابط إلى أن البروتين المتحور سيعاني للارتباط بفعالية بالمستقبل الخاص به، مما قد يؤدي إلى قطع خط الاتصال اللازم للنمو وإنتاج الحليب.

تخلص الورقة البحثية إلى أنه بينما تستند هذه النتائج إلى محاكاة حاسوبية ولم تُختبر بعد في أبقار حية، فإن طفرة L105R تبرز كأكثر المرشحين احتمالاً لتعطيل وظيفة بروتين IGF-1 البقري. ويقترح المؤلفون أن هذه الطفرة المحددة يمكن أن تعمل كعلامة جزيئية لبرامج التربية، مما يساعد المزارعين على اختيار الأبقار الأقل عرضة لحمل هذا "الخطأ المطبعي" الجيني والأكثر عرضة للإنتاجية. لا تدعي الدراسة أنها حلت مشكلة إنتاج الألبان، لكنها تقدم مشتبهاً به ذا احتمالية عالية في حالة الخلل الوظيفي الجيني، مما يمهد الطريق للاختبارات الواقعية المستقبلية.