Prioritization of Androgenetic Alopecia-Associated Lactylation Sites through Multi-Omics Integration and the HybridKla Deep Learning Framework
تدمج هذه الدراسة بيانات الأوميكس المتعددة من أنسجة فروة الرأس مع إطار تعلم عميق هجين (HybridKla) لتحديد مواقع اللاكتلة المرتبطة بالصلع الوراثي، وتحديد ZBTB7B K206 وK212 كمرشحين رئيسيين لاستقصاء التنظيم ما بعد الترجمة في البيئة المناعية للجريبات.
غالبًا ما يُنظر إلى تساقط الشعر على أنه مشكلة تجميلية، ولكن بالنسبة للملايين من الأشخاص المتأثرين بالصلع الوراثي (androgenetic alopecia)، فإنه يمثل لغزًا بيولوجيًا معقدًا. هذه الحالة، المعروفة باسم الصلع النمطي لدى الذكور أو الإناث، تتضمن انكماش بصيلات الشعر بمرور الوقت حتى لا تعود قادرة على إنتاج شعر مرئي. وبينما درس العلماء لفترة طويلة الجينات والبروتينات المشاركة في هذه العملية، برز مؤخرًا طبقة جديدة من البيولوجيا لتلفت الأنظار. فالخلايا لا تعتمد فقط على التعليمات المكتوبة في حمضها النووي (DNA)؛ بل تستخدم أيضًا علامات كيميائية للارتباط ببروتيناتها، تعمل مثل المفاتيح التي تشغل وظائف معينة أو تطفئها. إحدى هذه العلامات، وتسمى "اللاكتلة" (lactylation)، تربط استقلاب طاقة الخلية بكيفية سلوكها. فعندما تنتج الخلية "اللاكتات"، وهي ناتج ثانوي لاستخدام الطاقة، يمكنها ربط هذا الجزيء بالبروتينات، مما قد يغير طريقة عمل تلك البروتينات. إن فهم هذه التغييرات الكيميائية الدقيقة قد يكشف عن سبب توقف نمو الشعر في مناطق معينة من فروة الرأس، مما يقدم طريقة جديدة للنظر إلى المرض.
قرر فريق من الباحثين في جامعة تونغجي في شنغهاي استقصاء هذا الاحتمال من خلال فحص فروة الرأس نفسها. حيث جمعوا عينات من أنسجة من مقدمة ومؤخرة رؤوس عدة رجال يعانون من تساقط الشعر. تمت مقارنة مقدمة الرأس، حيث يخف الشعر عادةً، بمؤخرة الرأس، حيث يظل الشعر كثيفًا عادةً. استخدم العلماء أدوات متقدمة لقراءة النشاط الجيني للخلايا الفردية، وقياس إجمالي كمية البروتينات الموجودة، والبحث تحديدًا عن تلك البروتينات التي تحمل علامة "اللاكتلة". لم يكن هدفهم إيجاد علاج فورًا، بل تحديد دلائل جزيئية محددة يمكن أن تفسر ما يحدث داخل البصيلة. لقد ركزوا على بروتين يسمى ZBTB7B، المعروف بدوره في تنظيم الخلايا المناعية، لأن دراسات سابقة أشارت إلى أن النشاط المناعي يلعب دورًا في تساقط الشعر.
وجد الباحثون تناقضًا صارخًا في بياناتهم. ففي الشعر الخفيف في مقدمة فروة الرأس، كان إجمالي كمية بروتين ZBTB7B أعلى في الواقع مما هو عليه في الشعر الكثيف في الخلف. ومع ذلك، عندما بحثوا تحديدًا عن علامات "اللاكتلة" المرتبطة بهذا البروتين، وجدوا اتجاهًا عكسيًا. فالمواقع المحددة على البروتين حيث ترتبط هذه العلامات عادةً كانت أقل شيوعًا بكثير في الشعر الخفيف. كان الأمر كما لو أن المصنع ينتج المزيد من الآلات، ولكن عددًا أقل من هذه الآلات يحمل الأداة المحددة التي قد تغير طريقة عملها. كان هذا النمط ثابتًا عبر مختلف المتبرعين، مما يشير إلى أنه إشارة بيولوجية حقيقية وليس خطأً عشوائيًا. وقد تم تحديد موقعين محددين على البروتين حيث ظهر هذا الاختلاف بشكل أوضح وهما K206 وK212.
ولفهم هذه النتائج، لجأ الفريق إلى نموذج حاسوبي مصمم للتنبؤ بالمواقع التي من المرجح أن تظهر فيها علامات "اللاكتلة" على البروتين. لقد قاموا بتكييف إطار عمل ذكاء اصطناعي قوي، صُمم في الأصل للتعرف على الأنماط في تسلسلات البروتين، لتحليل بياناتهم الخاصة. ومن خلال تدريب النموذج على آلاف الأمثلة المعروفة ثم اختباره على البروتينات من دراستهم، أكدوا أن المواقع K206 وK212 هي بالفعل مرشحات محتملة لعملية "اللاكتلة". أشار النموذج الحاسوبي، جنبًا إلى جنب مع الأدلة المادية من المختبر، إلى هذين الموقعين باعتبارهما أكثر الأماكن واعدة للبحث القادم. كما تحقق الباحثون من بنية البروتين باستخدام خريطة رقمية لشكلة، ملاحظين أن هذه المواقع تقع في منطقة مرنة بين جزأين وظيفيين رئيسيين من الجزيء، مما يجعلها مواقع معقولة للتنظيم.
لم تثبت الدراسة أن "اللاكتلة" تسبب تساقط الشعر، كما لم تثبت أن إصلاح هذه العلامات سيؤدي إلى إعادة نمو الشعر. بدلاً من ذلك، قدمت قائمة أولويات واضحة للتجارب المستقبلية. يشير العمل إلى أن دور الجهاز المناعي في تساقط الشعر قد يتم التحكم فيه من خلال هذه المفاتيح الكيميائية الصغيرة على بروتين محدد. ومن خلال تضييق التركيز من الجينوم أو البروتينوم بأكمله إلى مجرد موقعين محددين على بروتين واحد، قدم الباحثون للمجتمع العلمي نقطة انطلاق ملموسة. لقد أظهروا أن النظر في إجمالي كمية البروتين ليس كافيًا؛ بل يجب أيضًا النظر في كيفية تعديل ذلك البروتين كيميائيًا. يربط هذا النهج بين البيئة الواسعة لأنسجة فروة الرأس والتغيرات الجزيئية الدقيقة التي تحدث داخل الخلايا الفردية، مما يوفر خريطة أكثر تفصيلًا لفهم بيولوجيا تساقط الشعر.
ملخص تقني: تحديد أولويات مواقع اللاكتيلة المرتبطة بالصلع الأندروجيني من خلال التكامل متعدد الأوميات وإطار التعلم العميق HybridKla
بيان المشكلة يتميز الصلع الأندروجيني (AGA) بتصغير حجم الجريبات الشعرية، ومع ذلك، لا تزال الأساسات الجزيئية لهذه العملية داخل البيئة المجهرية للجريب غير مفهومة بشكل كامل. وبينما حددت الدراسات الحديثة للنسخ المتماثل أحادي الخلية والنسخ المكاني تغيرات في حالات الخلايا والتواصل بين الخلايا، فإن البيانات النسخية وحدها لا تلتقط وفرة البروتين أو حالات تعديل ما بعد الترجمة (PTM). وتحديداً، فإن مساهمة اللاكتيلة على الليسين (Kla) — وهي تعديل ما بعد ترجمي (PTM) يربط استقلاب اللاكتات بتنظيم الجينات — في الصلع الأندروجيني لا تزال غير مشخصة بشكل جيد. ويكمن التحدي في التمييز بين ما إذا كانت التغيرات في وفرة الببتيدات المعدلة تعكس تغييرات حقيقية في تنظيم المواقع أم أنها مجرد تحولات في مستويات البروتين الكلي، وفي تحديد أولويات مواقع اللاكتيلة المحددة للتحقيق الوظيفي وسط تعقيد عينات الأنسجة الكاملة.
المنهجية استخدمت الدراسة سير عمل متكامل متعدد الأوميات يجمع بين التوصيف التجريبي والتنبؤ بالتسلسل القائم على التعلم العميق:
توليد بيانات متعددة الأوميات:
النسخ المتماثل أحادي الخلية: تحليل ست عينات مزدوجة من فروة الرأس الأمامية (FH) والقذالية (OC) من ثلاثة متبرعين. تم استخراج مجموعة خلوية محددة من الخلايا التائية (العنقود C21) للتحليل شبه المجمع (pseudobulk) لتحديد السياق النسخي.
البروتيوميات واللاكتيلوميات: إجراء تحليل طيف الكتلة بالاستحواذ المستقل للبيانات (DIA) على عينات FH وOC مزدوجة من متبرعين اثنين. شمل ذلك توصيف البروتين الكلي وإثراء الأجسام المضادة لليسين-لاكتيلة لتوصيف اللاكتيلوم.
القياس الكمي: تم قياس وفرة البروتين الكلي والببتيدات الملاكتيلة، وتطبيعها، ومقارنتها بين المناطق (FH مقابل OC) لتحديد الاتجاهات المتضاربة حيث أظهر البروتين الكلي والببتيدات المعدلة اتجاهات متعاكسة.
إطار التعلم العميق (تكييف HybridKla):
بنية النموذج: كيفت الدراسة إطار HybridKla، الذي يدمج ESM2 (نموذج لغة بروتيني)، وشبكة الذاكرة الطويلة قصيرة المدى (LSTM)، والميزات الفيزيائية والكيميائية والتركيبية المصنوعة يدوياً.
مقارنة التمثيل: تم تقييم خمس استراتيجيات تمثيل لـ ESM2 (رمز CLS، وتضمين الليسين المركزي، ومتوسط القناع، والدمج) لتحديد المدخل الأمثل للتنبؤ بموقع اللاكتيلة.
دمج الدرجات: قام نموذج "5D MetaModel" بدمج تنبؤات ESM2-E4، وLSTM، وثلاثة نماذج ميزات مصنوعة يدوياً (AAINDEX، وCKSAAP، وOBC).
استراتيجية التكييف: لتحسين الحساسية للبيانات المختبرية المحددة، تم تجميد هيكل ESM2، مع إجراء ضبط دقيق لرأس التصنيف فقط باستخدام مجموعة فرعية من بيانات اللاكتيلوم الخاصة بالدراسة (تكييف البيانات المختبرية).
تحديد أولويات المرشحين:
تم تحديد أولويات المرشحين بناءً على ثلاثة خطوط متقاربة من الأدلة: (1) الاتجاهات الإقليمية المتضاربة في وفرة البروتين الكلي مقابل الببتيد المعدل في الأنسجة ذات الصلة بالصلع الأندروجيني؛ (2) درجات التنبؤ العالية من نموذج التعلم العميق المكيّف؛ و(3) احتمالات التمركز المرتفعة في طيف الكتلة.
تم تقييم السياق الهيكلي باستخدام نماذج AlphaFold لتصور مواضع البقايا وثقة الموضع المحلية (pLDDT).
النتائج الرئيسية
السياق النسخي: كشف تحليل شبه المجمع لمجموعة الخلايا التائية C21 أن وحدة "مرتبطة بالمساعد" (بما في ذلك ZBTB7B) أظهرت زيادة ثابتة، وإن لم تكن ذات دلالة إحصائية، في فروة الرأس الأمامية عبر جميع المتبرعين الثلاثة.
اتجاهات الوفرة المتضاربة: في أزواج المتبرعين الاثنين، كانت وفرة بروتين ZBTB7B الكلية أعلى في فروة الرأس الأمامية (نسبة FH/OC = 1.806)، بينما أظهرت الببتيدات الملاكتيلة المرتبطة بالمواقع K206 وK212 وفرة أقل في فروة الرأس الأمامية (نسب FH/OC بلغت 0.523 و0.587 على التوالي). سلط هذا الاتجاه المتضاد الضوء على ضرورة التحقيق على مستوى الموقع بدلاً من تحليل البروتين الكلي.
أداء النموذج:
حققت تمثيلية E4 (دمج الليسين المركزي ومتوسط القناع) أعلى متوسط لـ ROC-AUC (0.8222) بين متغيرات ESM2.
حقق نموذج 5D MetaModel قيمة ROC-AUC قدرها 0.8281 وPR-AUC قدرها 0.7984 على مجموعة البيانات المرجعية.
أدى تكييف البيانات المختبرية إلى زيادة كبيرة في متوسط الاستدعاء للمواقع الإيجابية للبروتينات غير المرئية سابقاً من 78.09% إلى 90.75%، رغم أن هذا حسن حساسية الكشف بدلاً من خصوصية التصنيف الإجمالية.
تحديد المرشحين: تم التنبؤ بكل من K206 وK212 كمواقع إيجابية بواسطة 5D MetaModel. أكدت سجلات طيف الكتلة وجود الببتيدات الملاكتيلة في هذه المواقع باحتمالات تمركز عالية (K206: 0.896–1.000؛ K212: 1.000 عبر الجولات). أشارت المحاكاة الهيكلية إلى أن هذه البقايا تقع في منطقة منخفضة الثقة بين نطاق BTB/POZ ونطاق C2H2 zinc-finger.
المساهمات الرئيسية
سير عمل تحديد أولويات المواقع: تؤسس الدراسة سير عمل يربط الملاحظات متعددة الأوميات على مستوى الأنسجة (النسخ المتماثل، البروتيوميات، اللاكتيلوم) مع تنبؤات تسلسل البقايا لتحديد أولويات مواقع تعديل ما بعد الترجمة (PTM) المرشحة.
تحديد ZBTB7B K206/K212: يحدد العمل ZBTB7B K206 وK212 كمواقع لاكتيلة مرشحة في الصلع الأندروجيني، مشيراً إلى ملفها الفريد من اتجاهات الوفرة المتضاربة بالنسبة للبروتين الكلي.
التكييف المنهجي: توضح الورقة فائدة تكييف نماذج لغة البروتين المدربة مسبقاً (HybridKla) عبر الضبط الدقيق لرأس التصنيف لتحسين الاستدعاء لمجموعات البروتينات المتبقية المستمدة من نفس المصدر المختبري.
التكامل السياقي: تربط الدراسة الدور الراسخ لـ ZBTB7B في تنظيم سلالة الخلايا التائية مع إمكانية التنظيم عبر تعديل ما بعد الترجمة بواسطة اللاكتيلة، مما يوفر مدخلاً جزيئياً لاستقصاء الحالات المناعية في البيئة المجهرية للجريب.
الأهمية والادعاءات يذكر المؤلفون أن المساهمة الرئيسية لهذا العمل ليست إنشاء آلية مرضية سببية، بل توفير أساس مركز لفحص تنظيم ما بعد الترجمة ضمن البيئة المجهرية المناعية للجريب. ومن خلال دمج البيانات متعددة الأوميات مع التعلم العميق، تنقل الدراسة أبحاث الصلع الأندروجيني من الأوصاف الواسعة لاختلافات التعبير نحو فرضيات قابلة للاختبار ومحددة المواقع فيما يتعلق بدوران البروتين، وتجميع المعقدات التنظيمية، وتنظيم الكروماتين.
تصيغ الورقة نتائجها بتواضع باعتبارها "تحديد أولويات المواقع المرشحة" و"رؤى جزيئية". كما تقر صراحة بالقيود، بما في ذلك العدد الصغير من المتبرعين البيولوجيين، وغياب الضوابط السليمة، وحقيقة أن بيانات فروة الرضاء الكاملة لا يمكنها تعيين الأصل الخلوي للتعديلات بشكل نهائي. وتخلص الدراسة إلى أن التضارب المحدد بين وفرة البروتين الكلي والببتيد الملاكتيل لـ ZBTB7B يستوجب مزيداً من التحقق المستهدف لتحديد ما إذا كانت هذه المواقع المحددة تؤثر على النشاط النسخي لـ ZBTB7B أو نمو الجريبات الشعرية.