An in silico pipeline prioritizes candidate viral and bacterial mimics of hybrid insulin peptide autoantigens in type 1 diabetes
تقدم هذه الدراسة مساراً حاسوبياً يعتمد على الهيكل أولاً، يحدد ويرتب عشرة ببتيدات مرشحة فيروسية وبكتيرية تحاكي المستضدات الذاتية لببتيد الأنسولين الهجين في مرض السكري من النوع الأول، وذلك عبر دمج فحص الأنماط، وتوقعات الارتباط بـ HLA، والنمذجة الهيكلية القائمة على AlphaFold 3، مما يسلط الضاح في النهاية على ثلاثة ببتيدات محددة لفيروس سارس-كوف-2 وفيروس رينو البشري 3 للتحقق التجريبي.
يُعد مرض السكري من النوع الأول حالة يقوم فيها جهاز المناعة في الجسم بمهاجمة الخلايا التي تنتج الأنسولين عن طريق الخطأ. لعقود من الزمن، عرف العلماء أن الوراثة تلعب دوراً رئيسياً في تحديد من يصاب بالمرض، وتحديداً جينات معينة تساعد الجهاز المناعي على التعرف على الغزاة الأجانب. ومع ذلك، ظل المحفز الدقيق الذي يبدأ هذا الهجوم لغزاً. تشير إحدى النظريات الرائدة إلى أن فيروساً أو بكتيريا قد يخدع الجهاز المناعي؛ ويحدث هذا عندما يشبه جزء من ميكروب ما جزءاً من الجسم بشكل كبير، بحيث تبدأ الخلايا المناعية، المدربة على صيد الميكروب، في مهاجمة أنسجة الجسم نفسها عن طريق الخطأ. وقد أضاف اكتشاف أحدث لمسة من التعقيد إلى هذه القصة: فأحياناً، يقوم الجسم بإنشاء بروتينات هجينة غريبة عن طريق لصق قطع من بروتينات مختلفة معاً. هذه "ببتيدات الأنسولين الهجينة" لا توجد في المخطط الجيني الأصلي للجسم، مما يعني أن الجهاز المناعي لم يتعلم أبداً تجاهلها خلال تدريبه المبكر. وإذا صادف وجود ميكروب يمتلك قطعة تشبه أحد هذه الهجن النادرة، فقد يتجاوز ضوابط السلامة في الجسم ويثير المرض.
قام باحث في معهد بناروييا للأبحنت (Benaroya Research Institute) ببناء مسار رقمي للبحث عن هذه الأشكال الشبيهة الخطيرة. وبدلاً من مجرد مقارنة التسلسل الحرفي للبروتينات، وهو أمر قد يكون مضللاً، يركز هذا النهج الجديد على الشكل ثلاثي الأبعاد للجزيئات. بدأ الباحث بخمسة أمثلة معروفة حيث تعرفت الخلية المناعية للمريض بنجاح على ببتيد أنسولين هجين. وباستخدام نمذجة حاسوبية متقدمة، قاموا برسم خرائط توضح بالضبط أي أجزاء من الببتيد تلامست مع الخلية المناعية والبروتين الذي يحمله. ثم قاموا بترجمة هذه الاتصالات الفيزيائية إلى مجموعة من القواعد التي تصف الخصائص الكيميائية المطلوبة في نقاط محددة. استُخدمت هذه القواعد لمسح المخططات الجينية لعشرة أنواع مختلفة من الفيروسات والبكتيريا المعروفة بارتباطها بالسكري. بحث الكمبيوتر عن أي قطع ميكروبية تتناسب مع شكل وملف تعريف الهجين الخطير، بغض النظر عما إذا كانت تسلسلاتها الحرفية متطابقة أم لا.
بدأ البحث بأكثر من ألف تطابق محتمل. ثم طبق الباحث سلسلة من المرشحات الصارمة لتضييق القائمة. أولاً، تنبأ الكمبيوتر بما إذا كانت هذه القطع ستلتصق بالفعل بمستقبل المناعة المحدد المرتبط بخطر الإصابة بالسكري. بعد ذلك، قام بنمذجة المعقد بأكمله لمعرفة ما إذا كان الشكل سيبقى متماسكاً ومستقراً. وأخيراً، قارن شكل كل مرشح بالببتيد الهجين الأصلي لضمان تشابههما هيكلياً. أدت هذه العملية إلى استبعاد معظم المرشحين، مما ترك مجموعة صغيرة اجتازت كل الاختبارات. ومن المثير للاهتمام، وجدت الدراسة أن مجرد النظر إلى مدى تشابه الهيكل العظمي (backbone shape) لم يفصل بين المرشحين الجيدين والسيئين؛ حيث بدت نماذج "الخداع" (decoy controls) المستخدمة في الاختبار مشابهة تماماً للمرشحين الحقيقيين. وهذا يشير إلى أن تشابه الشكل وحده ليس كافياً لإثبات وجود محاكٍ (mimic).
بعد استبعاد المرشحين القادمين من فيروسات خاصة بالحيوانات، تبقى لدى الباحث عشرة منافسين أقوياء للاختبار التجريبي. برز اثنان منها كأكثر المسارات واعدة: أحدهما قطعة من فيروس سارس-كوف-2 (SARS-CoV-2)، المسبب لمرض كوفيد-19، والآخر يأتي من فيروس نزلات البرد الشائع المعروف باسم فيروس رينو البشري 3 (human rhinovirus 3). أظهرت قطعة سارس-كوف-2 درجة ثقة عالية جداً في النماذج الحاسوبية وطابقت شكل ببتيد الهجين بشكل جيد. أما قطعة فيروس رينو، ورغم أنها كانت أقل ثقة بقليل من حيث الشكل، فقد أظهرت طاقة ربط كيميائي أقوى، مما يشير إلى أنها قد تلتصق بقوة شديدة بمستقبل المناعة. لم يعلن الباحث عن فائز نهائي واحد؛ بل يتم التعامل مع كليهما كأولويات متساوية. وتخضع هاتان القطعتان الآن للاختبار في المختبر لمعرفة ما إذا كانتا قادرتقتين فعلياً على تحفيز الخلايا المناعية التي تقود مرض السكري من النوع الأول.
توفر الدراسة قائمة ملموسة من الأهداف للعلماء للتحقيق فيها، لكنها تضع أيضاً حدوداً واضحة لما يمكن للنماذج الحاسوبية وما لا يمكنها فعله. تُظهر النتائج أنه بينما نجح المسار الرقمي في تحديد عشرة مرشحين محتملين، إلا أنه لم يستطع إثبات أن هذه المحاكيات أكثر شيوعاً في الطبيعة مما قد يوحي به الاحتمال العشوائي. وتستبعد الدراسة صراحةً فكرة أن مقياساً بسيطاً لتشابه الهيكل العظمي يمكن أن يميز بين المحاكي الحقيقي وبين ببتيد عشوائي. إن قائمة العشرة ببتيدات النهائية هي نقطة انطلاق للتجارب الواقعية، وليست إجابة نهائية. يسلط هذا العمل الضوء على أنه بينما يمكن للمحاكاة الحاسوبية توجيه البحث عن أسباب الأمراض المعقدة، فإن الإثبات النهائي يتطلب اختباراً فيزيائياً في المختبر. وتخضع حالياً المرشحتان الرائدتان، إحداهما من فيروس كورونا والأخرى من فيروس نزلات البرد الشائع، لهذه المرحلة التجريبية الحاسمة لتحديد ما إذا كانت هما بالفعل المفتاح لفك لغز ما يحفز مرض السكري من النوع الأول.
ملخص تقني: خط معالجة حاسوبي (in silico) يعطي الأولوية لمرشحات من المحاكيات الفيروسية والبكتيرية لبيبتينات المناعة الذاتية لـ "ببتين الهجين للأنسولين" في مرض السكري من النوع الأول
بيان المشكلة مرض السكري من النوع الأول (T1D) هو مرض مناعة ذاتية مدفوع بتدمير الخلايا بيتا البنكرياسية بواسطة الخلايا التائية (T-cells). وبينما ترتبط القابلية الوراثية بأنماط جينية من فئة HLA class II (تحديداً DR3-DQ2 و DR4-DQ8)، لا تزال المسببات البيئية التي تحفز المناعة الذاتية غير واضحة. وتوجد فجوة حرجة في فهم كيفية قيام العدوى الميكروبية بتحفيز السكري من النوع الأول عبر "المحاكاة الجزيئية" (molecular mimicry). وتعد بيبتينات الهجين للأنسولين (HIPs)—وهي وصلات اندماجية تتشكل عن طريق الانقسام البروتيني لبروتين الأنسولين الأولي وبروتينات جزيرة أخرى—من المستضدات ذاتية المناعة في مرض السكري من النوع الأول. ولأن هذه التسلسلات الوصلية لا يتم تشفيرها في الخط الجرثومي (germline)، فقد تكون غائبة أثناء عملية الانتقاء السلبي في الغدة الصدغية (thymic negative selection)، مما يسمح للخلايا التائية المتفاعلة مع الـ HIP بالهروب من التسامح المركزي. ويتمثل التحدي في تحديد البيبتينات الميكروبية التي تحاكي بنيوياً وصلات الـ HIP هذه، مما قد يؤدي إلى تهيئة خلايا تائية متصالبة التفاعل. إن عمليات البحث التقليدية القائمة على التسلسل (مثل BLAST) غير كافية لأن مستقبلات الخلايا التائية (TCRs) تتعرف على التضاريس السطحية ثلاثية الأبعاد بدلاً من تطابق التسلسل الخطي.
المنهجية تقدم الدراسة خط معالجة يعتمد على مبدأ "البنية أولاً" (structure-first) مصمم لتحديد المحاكيات الميكروبية لخمسة أنظمة مرجعية من (HIP–TCR–MHC)، مقيدة بـ HLA-DQ8 أو DQ8-trans. تسير سير العمل كما يلي:
إعادة بناء المراجع: تم إعادة بناء خمسة معقدات (TCR–HIP–MHC) باستخدام تسلسلات جينات V/J وتسلسلات CDR3 المنشورة. تم تجميع سلاسل TCR باستخدام Stitchr، مع تقليم المناطق الثابتة.
تعريف الفارماتوفور (Pharmacophore): باستخدام نماذج AlphaFold 3 للمعقدات المرجعية، تم تحديد البقايا (residues) الواقعة ضمن نطاق 3.0 أنجستروم من ذرات TCR وMHC. تم تعميم هذه الاتصالات إلى فئات فيزيائية كيميائية لإنشاء أنماط ScanProsite، مع التركيز على التكافؤ البنيوي بدلاً من تطابق التسلسل.
مسح البروتيوم: تم تطبيق الأنماط على قاعدة بيانات Swiss-Prot عبر عشر مجموعات تصنيفية مرتبطة بالسكري من النوع الأول (ست مجموعات فيروسية، وثلاث أنواع بكتيرية، ومجموعة واحدة غير ممثلة)، مما أسفر عن 1,092 ضربة نمطية (motif hits).
التصفية والنمذجة:
ارتباط MHC: تمت تصفية الضربات باستخدام NetMHCIIpan-4.2 (رتبة EL < 5%).
موثوقية المعقد: تم نمذجة البيبتينات الناجية باستخدام AlphaFold 3 في معقد مع الـ TCR والـ MHC ذوي الصلة. تم الاحتفاظ فقط بالمرشحات التي سجلت درجة (ipTM) للواجهة > 0.75، مما نتج عنه 55 مرشحاً.
نمذجة TCR: تم نمذجة الـ 55 مرشحاً بشكل أكبر باستخدام TCRModel2 مع خمس بذور عشوائية لتقييم التباين البنيوي.
بوابات تحديد الأولويات:
بوابة الـ HIP الأصلية (Native-HIP): تمت مقارنة المرشحين مع مرجع الـ HIP الأصلي باستخدام بوابة الثلاثة معايير: ipTM لـ (TCR–pMHC)، وinterface-pLDDT، وطاقة جيبس الحرة (ΔG). كان على المرشحين استيفاء أو تجاوز فترات الثقة بنسبة 95% للمرجع الأصلي لهذه المقاييس.
الارتباط البشري: تم استبعاد البيبتينات المشتقة من السلالات المرتبطة بالحيوانات (تحديداً Rotavirus A).
الضوابط والتحقق: استخدمت الدراسة ضوابط بروتيومية مشوشة (shuffled-proteome) محافظة على التركيب لاختبار إثراء عائد النمط، وتضمنت ضوابط سالبة (اصطناعية وفيروسية) لاختبار توزيعات جذر متوسط مربع الانحراف (RMSD) للهيكل العظمي.
النتائج الرئيسية
عائد المرشحين: قلل خط المعالجة من 1,092 ضربة نمط أولية إلى 55 مرشحاً بعد تصفية ارتباط HLA وثقة AlphaFold 3.
الضوابط السالبة:
عائد النمط: لم يظهر ضابط البروتيوم المشوش المطابق للتركيب أي إثراء إجمالي لضربات النمط (نسبة الملاحظ/العدمي 0.97)، مما يشير إلى أن الأنماط تعمل كفلاتر قابلية للمعالجة (tractability filters) بدلاً من إظهار إثراء متأصل في البروتيوم.
RMSD الهيكل العظمي: وقعت الضوابط السالبة ضمن توزيع RMSD الهيكل العظمي للمرشحين (المئويات 21.8، 61.8، و81.8). لم يستطع أي قطع واحد لـ RMSD فصل المرشحين عن الضوابم، مما يشير إلى أن هندسة هيكل Class II MHC المشتركة هي قيد بنيوي وليست مميزاً للمحاكاة.
تحديد الأولويات: احتفظت بوابة الـ HIP الأصلية ذات الثلاثة معايير بـ 13 مرشحاً. بعد استبعاد ثلاثة بيبتينات من فيروس روتالفيروس A المرتبط بالحيوانات، تم ترشيح عشرة بيبتينات للاختبار التجريبي.
أفضل المرشحين: ظهرت بيبتين كأولويات مشتركة للمتابعة التجريبية:
RRNVATLQAENVTG من SARS-CoV-2.
SDLLKSVDSEEVRE من Human Rhinovirus 3.
أظهر بيبتين SARS-CoV-2 متوسط ipTM أعلى وتغيرات اتجاهية مواتية عبر جميع المقايلات الثلاثة.
أظهر بيبتين Rhinovirus 3 طاقة واجهة (ΔG) أكثر ملاءمة.
تشير الورقة إلى أن تفوق أي منهما على الـ HIP الأصلي لم يتم إثباته إحصائياً، وأن ترتيبهما النسبي ليس نهائياً.
المتانة الإحصائية: أظهر 48 من أصل 55 مرشحاً درجات ipTM أعلى بكثير من ضوابئهم السالبة المحددة (اختبار Welch's t-test, p<0.05). ومع ذلك، يشير المؤلف إلى أن هذه المقارنة مشوبة بالارتباك لأن المرشحين تم تصفيتهم مسبقاً للحصول على ipTM عالٍ، بينما لم يتم تصفية الضوابط.
الأهمية والادعاءات تدعي الورقة أنها توفر خط معالجة كمي يعتمد على البنية يولد مجموعة محددة وقابلة للاختبار من عشرة بيبتينات ميكروبية للتحقق التجريبي كمحاكيات للمستضدات ذاتية المناعة في السكري من النوع الأول.
ومن المهم أن المؤلف يصيغ مساهمات الدراسة بتواضع كتعريف للمرشحات والحدود معاً:
الفرضيات القابلة للاختبار: يوفر خط المعالجة تركيبات (بيبتين–TCR–MHC) محددة (لا سيما مرشحي SARS-CoV-2 و Rhinovirus 3) للاختبار التجريبي في معهد بناروي للأبحاث.
الحدود المنهجية: توثق الدراسة أن:
RMSD الهيكل العظمي ليس فلتراً مستقلاً للمحاكاة في هذا السياق، حيث تتداخل الضوابط بشكل كبير مع المرشحين.
عمليات البحث عن الأنماط المحافظة على التركيب لا تنتج إثراءً إجمالياً فوق الخلفيات المشوشة، مما يشيع أن الأنماط الفيزيائية الكيميائية الخطية ضرورية ولكنها غير كافية لتحديد المحاكيات الحقيقية دون نمذجة بنيوية لاحقة.
الفصل بين المرشحين والضوابم في درجات ipTM مشوب بالارتباك بسبب خطوة التصفية المسبقة، ولا يشكل معيار تمييز تم التحقق منه.
تتجنب الورقة صراحةً الادعاء بأن هذه البيبتينات تسبب السكري من النوع الأول أو أن SARS-CoV-2 مرتبط بالتأكيد بالمرض عبر هذه الآلية. بدلاً من ذلك، فهي تقدم فرضية ضيقة وقابلة للاختبار: وهي أن بيبتينات فيروسية محددة يمكنها محاكاة مستضدات الـ HIP بنيوياً وطاقياً في سياق HLA-DQ8 ومستقبلات TCR المشتقة من المرضى. يعمل هذا العمل كجسر بين البيولوجيا البنيوية الحاسوبية وعلم المناعة التجريبي، حيث يعطي الأولوية للأهداف مع توثيق حدود المقاييس الحالية للحوسبة الحيوية بدقة.