A multi-layer computational framework for structural prioritization of transcriptome-derived proteins

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

وهنا يأتي دور دراسة جديدة، تقدم طريقة منهجية للتنقل عبر هذا الجبل من البيانات. فقد طور باحثون من جامعة ABC الفيدرالية في البرازيل إطارًا حوسبيًا متعدد الطبقات مصممًا لغربلة آلاف البروتينات المرشحة وتحديد أكثرها واعدة لمزيد من الدراسة. وقد ركزوا على نبات "سيريوس جاماكارو" (Cereus jamacaru)، وهو صبار موطنه الأصلي منطقة "كواتينجا" في البرازيل، والمعروف بإمكانياته المضادة للالتهابات. وبينما ركزت الأبحاث السابقة على هذا النبات في الغالب على المركبات الكيميائية الصغيرة، أراد هذا الفريق استكشاف البروتينات الكامنة داخل شفرته الجينية. لم يكن هدفهم إثبات أن بروتينات النبات هذه تعالج الأمراض مباشرة، بل إنشاء طريقة موثوقة وخطوة بخطيبه لإيجاد مرشحين بروتينيين ذوي بنية مستقرة يمكنهم التفاعل مع إشارات الالتهاب البشرية. ومن خلال الجمع بين عدة تقنيات حوسبية مختلفة، نجحوا في تقليص قائمة تضم أكثر من 129,000 بروتين مترجم إلى تسعة مرشحين فقط عاليي الثقة، مما أظهر استراتيجية قابلة للتوسع يمكن تطبيقها على البيانات الجينية لأي كائن حي.

بدأت الرحلة بمجموعة بيانات ضخمة تحتوي على 129,052 تسلسلًا بروتينيًا مشتقًا من جذور وبشرة الصبار. ولإدراك هذا العدد الهائل، طبق الباحثون أولاً مرشحًا واسع النطاق يعتمد على الهياكل المعروفة. فقد قارنوا بروتينات الصبار بقاعدة بيانات للأشكال البروتينية المحددة تجريبيًا، مع الاحتفاظ فقط بتلك التي تشترك في أكثر من نصف تسلسلها مع هيكل معروف. قللت هذه الخطوة الأولية المجموعة إلى حوالي 14,800 مرشح. ثم استخدموا نظام تسجيل لتقييم جودة الأشكال المتوقعة، واستبعدوا تلك التي بدت غير مستقرة أو غير موثوقة. أدت هذه العملية إلى تقليص القائمة إلى حوالي 3,500 نموذج هيكلي عالي الجودة، مما خلق مجموعة يمكن إدارتها من المرشحين لإجراء تحليل أعمق.

بعد ذلك، وجه الفريق اهتمامهم إلى ثلاثة بروتينات بشرية محددة تُعرف باسم السيتوكينات: عامل نخر الورم ألفا (TNF-alpha)، وإنترلوكين-1 بيتا (IL-1 beta)، وإنترفيرون-ألفا (IFN-alpha). تعد هذه الجزيئات منظمات مركزية للاستجابة الالتهابية في الجسم وغالبًا ما تكون مرتبطة بحالات مثل التهاب المفاصل الروماتويدي. أراد الباحثون معرفة ما إذا كان أي من بروتينات الصبار يمكن أن يتفاعل فيزيائيًا مع هذه السيتوكينات البشرية. وللقيام بذلك، استخدموا استراتيجية "الرسو" (docking) ثنائية الخطوات. أولاً، استخدموا طريقة فحص سريعة تعتمد على الجسم الصلب (rigid-body) لاختبار مدى ملاءمة الـ 3,500 بروتين من الصبار مقابل السيتوكينات البشرية الثلاثة. حددت هذه الخطوة مجموعة أصغر من البروتينات التي أظهرت تطابقًا أوليًا واعدًا. ولتحسين هذه النتائج، استخدموا طريقة رسو ثانية أكثر مرونة تسمت بالسماح بحركات طفيفة في أشكال البروتين، محاكية الطريقة التي قد تعدل بها الجزيئات الحقيقية نفسها عند اجتماعها. ومن خلال البحث عن البروتينات التي قدمت أداءً جيدًا في كلتا الطريقتين، حدد الباحثون تسعة مرشحين أظهروا باستمرار توافقًا هيكليًا قويًا مع السيتوكينات البشرية.

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

ولاختبار ما إذا كانت هذه التفاعلات ستصمد في الظروف الواقعية، أجرى الباحثون عمليات محاكاة ديناميكية. وضعوا أكثر أربعة أزواج بروتينية واعدة في بيئة افتراضية مليئة بجزيئات الماء وحاكوا سلوكها بمرور الوقت. استغرقت عمليات المحاكاة هذه 100 نانو ثانية لكل منها، مع تكرارها ثلاث مرات لكل زوج لضمان إمكانية تكرار النتائج. كان الهدف هو معرفة ما إذا كانت البروتينات ستظل مقفلة معًا أم أنها ستنفصل مع حركتها واهتزازها. أظهرت النتائج أن التفاعلات ظلت مستقرة طوال عمليات المحاكاة؛ حيث حافظت البروتينات على شكلها ونقاط اتصالها، حتى أثناء حركتها في الماء الافتراضي. أكدت هذه الاستقرارية الديناميكية أن التفاعلات المتوقعة لم تكن مجرد لقطات ثابتة، بل مثلت معقدات قوية ومنطقية فيزيائيًا.

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

تكمن القيمة القصوى لهذا العمل في الإطار نفسه. لم يدّع الباحثون أنهم اكتشفوا دواءً جديدًا، بل طوروا طريقة موثوقة لإيجاد المرشحين. ومن خلال دمج طبقات متعددة من الأدلة — النمذجة الهيكلية، ونهج الرسو المزدوج، والتحليل التطوري، والمحاكاة الديناميكية — أنشأوا سير عمل يقلل من عدم اليقين في كل خطوة. يسمح هذا النهج للعلماء بالانتقال من مجموعات البيانات الضخمة وغير القابلة للإدارة إلى قائمة منسقة صغيرة من الأهداف عالية الثقة الجاهزة للاختبار التجريبي. وتخلص الدراسة إلى أن هذه الاستراتيجية قابلة للنقل، مما يعني إمكانية تطبيقها على أي "ترانسكريبتوم" (transcriptome) أو مجموعة بيانات بروتينية كبيرة لتحديد المرشحين ذوي البنية القوية. إنها تقدم مسارًا عمليًاًا لاستكشاف الموارد البيولوجية، وتحويل التدفق الهائل للبيانات الجينية إلى مجموعة مركزة من الفرضيات التي يمكن أن توجه أبحاث الكيمياء الحيوية المستقبلية وجهود هندسة البروتينات.