Scaling Quantum Optimisation Beyond Hardware Limits for Real-World Scientific Workloads: Genome Assembly on Current Quantum Hardware
تُثبت هذه الدراسة البحثية أن إطار تحسين التفكيك الذاتي للهاميلتوني (HADOF) يمكنه التغلب على قيود أجهزة الحوسبة الكمومية ذات الحجم المتوسط والمشوبة بالضجيج (NISQ) الحالية لتجميع جينوم بكتيريا الزائفة الزنجارية (*Pseudomonas aeruginosa*) المكون من 7.1 مليون زوج قاعدي بنجاح على أجهزة كمومية حقيقية، محققاً نسبة 99.348% من الجينوم، مما يثبت جدوى التحسين الكمومي القابل للتوسع لأعباء العمل العلمية واسعة النطاق.
المؤلفون الأصليون:G Sankar, N., Miliotis, G., Caton, S.
يحمل كل كائن حي دليل تعليمات فريداً خاصاً به، مكتوباً بشفرة مكونة من أربعة أحرف كيميائية هي: A وC وG وT. ولفهم كيفية عمل الكائن الحي، أو كيف يسبب المرض، أو كيف قد يقاوم الدواء، يجب على العلماء أولاً قراءة هذا الدليل بأكمله. ومع ذلك، لا تستطيع الآلات الحديثة قراءة الكتاب بأكمله دفعة واحدة؛ بل تقوم بدلاً من ذلك بتقطيع الحمض النووي (DNA) إلى ملايين القطع الصغيرة والمتداخلة، تماماً مثل تمزيق رواية إلى آلاف القطع الصغيرة. ويتمثل التحدي في تجميع الجينوم في أخذ هذه الشظايا المبعثرة وإعادة تجميعها في التسلسل الصحيح والمكتمل. إنها أحجية ضخمة حيث لا تكون القطع مبعثرة فحسب، بل غالباً ما تحتوي أيضاً على أخطاء وتكرارات مربكة تجعل من الصعب معرفة مكان كل قطعة. ويعد إتقان هذا الأمر أمراً حيوياً لتتبع الأمراض المعدية، وفهم السرطان، ومراقبة كيفية تطور البكتيريا للبقاء على قيد الحياة أمام المضادات الحيوية.
لعقود من الزمن، اعتمد العلماء على الحواسيب التقليدية لحل هذه الأحجية، باستخدام طرق مختصرة ذكية لتخمين الترتيب الأفضل للقطع. ولكن مع ازدياد حجم الجينومات وتعقيدها، تفشل هذه الطرق المختصرة أحياناً، مما يترك فجوات أو يخلق اتصالات خاطئة. ويتساءل مجال جديد من البحث عما إذا كانت الحواسيب الكمومية، التي تعمل وفق القوانق الفيزيائية الغريبة التي تحكم عالم الجسيمات دون الذرية، يمكنها حل هذه الألغاز بفعالية أكبر. والأمل هو أن تتمكن هذه الآلات من استكشاف ملايين الترتيبات المحتملة في وقت واحد للعثور على المسار الوحيد الأفضل. ومع ذلك، لا تزال الحواسيب الكمومية الحالية في مراحلها الأولى؛ فهي صغيرة، وهشة، وعرضة لارتكاب الأخطاء، مما يعني أنها لا تستطيع بعد التعامل مع مجموعات البيانات الضخمة المطلوبة في البيولوجيا الواقعية.
في دراسة حديثة، تناول الباحثون هذه المشكلة تحديداً من خلال محاولة تجميع جينوم بكتيريا Pseudomonas aeruginosa، وهي بكتيريا شائعة يمكن أن تسبب عدوى خطيرة، باستخدام حاسوب كمومي حقيقي. لم يحاول الفريق فرض الأحجية بأكملها على الآلة دفعة واحدة، وهو أمر كان سيكون مستحيلاً نظراً للحجم المحدود للحاسوب. بدلاً من ذلك، قاموا بتفكيك المشكلة إلى آلاف القطع الصغيرة القابلة للإدارة. واستخدموا طريقة جديدة تسمى "إطار عمل تحسين التفكيك التلقائي للهاملتوني" (Hamiltonian Auto Decomposition Optimisation Framework - HADOF)، والذي يعمل بمثابة منسق. يقوم هذا الإطار بتقسيم مهمة التجميع الضخمة إلى مشكلات فرعية صغيرة تناسب الشريحة الكمومية، وحلها واحدة تلو الأخرى، ثم دمج النتائج معاً. واختبر الباحثون هذا النهج على معالج كمومي حقيقي يحتوي على 133 بت كمومي، أو "كيوبت" (qubit)، باستخدام بيانات تسلسل الحمض النووي الفعلية لبكتيريا يبلغ حجم جينومها 7.1 مليون زوج من القواعد.
أظهرت النتائج أنه على الرغم من أن الحاسوب الكمومي لم يكن مثالياً بعد، إلا أنه استطاع إنتاج إجابة مفيدة بيولوجياً. فقد نجحت الطريقة المدعومة بالكم في إعادة بناء 99.348 بالمائة من جينوم البكتيريا، وهي نتيجة قريبة بشكل ملحوظ مما يمكن أن تحققه أفضل الحواسيب التقليدية. وكان التسلسل المعاد بناؤه نسخة شبه مثالية من الأصل، مع عدم وجود أقسام مكررة أو مفقودة تقريباً. ويعد هذا النجاح مهماً لأنه يثبت أن التحسين الكمومي يمكن أن يعمل على نطاق أكبر بكثير مما سبق، متجاوزاً النماذج الصغيرة والرمزية إلى كائن حقيقي ذي صلة سريرية. ووجد الباحثون أن الحاسوب الكمومي لم يجد ببساطة الدرجة الرياضية "الأفضل" في حساباته؛ بل إن الدرجة المثالية رياضياً أدت أحياناً إلى جينوم مكسور أو غير صحيح. وبدلاً من ذلك، كان مفتاح النجاح يكمن في النظر إلى بنية الحل — وتحديداً، عدد القطع التي تم الاحتفاظ بها في المسار النهائي — واستخدام ذلك لتحديد التجميع الصحيح من بين مجموعة واسعة من الاحتمالات.
كما سلطت هذه الدراسة الضوء على القيود الحالية للتكنولوجيا. فعندما قام الباحثون بتشغيل نفس المشكلة على محاكاة مثالية وخالية من الضجيج للحاسوب الكمومي، كانت النتائج أفضل قليلاً مما كانت عليه في الجهاز الحقيقي. فقد تسبب الجهاز الحقيقي في حدوث أخطاء جعلت الدرجات الرياضية الخام تبدو أسوأ والحلول أكثر تشرذماً. ومع ذلك، أظهر الباحثون أنه حتى مع وجود هذه الأخطاء، فإن الطريقة الأساسية كانت قوية بما يكفي لاستعادة جينوم عالي الجودة. واكتشفوا أن عدد القطع التي تم الاحتفاظ بها في التجميع النهائي كان مؤشراً قوياً على النجاح: فالحلول التي احتفظت بحد أدنى معين من القطع أدت دائماً تقريباً إلى جينوم كامل، بغض النظر عن الدرجة الخام التي أعطاها الحاسوب لهم. وهذا يشير إلى أنه في المستقبل، قد لا يحتاج العلماء إلى انتظار حواسيب كمومية مثالية للحصول على نتائج مفيدة؛ بل يمكنهم استخدام آلات حالية غير مثالية إذا توفرت لديهم الطريقة الصحيحة لتفسير البيانات.
يعمل هذا العمل كجسر بين الوعد النظري والتطبيق العملي. فهو يوضح أنه من خلال الجمع بين الحوسبة الكمومية والاستراتيجيات الكلاسيكية الذكية، يمكن للعلماء البدء في حل المشكلات البيولوجية التي كانت بعيدة المنال سابقاً. وقد أتاح الباحثون شفرتهم وبرمجياتهم وبياناتهم للجمهور، مما يسمح للآخرين باختبار هذه الطرق على كائنات مختلفة ومع آلات مختلفة. ورغم أن هذا ليس بديلاً للأدوات القياسية المستخدمة في المستشفيات والمختبرات، إلا أنه دليل واضح على أن التحسين الكمومي يمكنه التعامل مع التعقيد البيولوجي في العالم الحقيقي. ويتضمن المسار المستقبلي تحسين هذه الطرق للتعامل مع جينومات أكبر وحتى تصفية الأخطاء التي تنتجها الأجهزة الحالية حتماً. وفي الوقت الحالي، تقف هذه الدراسة كخطوة ملموسة نحو مستقبل قد تساعد فيه الحواسيب الكمومية في فك رموز التعليمات الأكثر تعقيداً للحياة، قطعة قطعة.
ملخص تقني: توسيع نطاق التحسين الكمي لما وراء حدود الأجهزة لأعباء العمل العلمية في العالم الحقيقي
بيان المشكلة يعد تجميع الجينوم، وخاصة التجميع من النوع "de novo" لبيانات التسلسل طويلة القراءة، تحديًا أساسيًا في علم الجينوم ذو تطبيقات حيوية في مراقبة الأمراض المعدية وجينوم السرطان. تتضمن المهمة إعادة بناء جينوم كامل من قراءات مجزأة، وهي عملية غالبًا ما تُنمذج كمشكلة تحسين توافقي واسعة النطاق على رسم بياني للتداخل (overlap graph). وبينما تعتمد الخوارزميات الكلاسيكية (مثل Unicycler) على تبسيط الرسم البياني الاستدلالي، إلا أنها تواجه قيودًا في القابلية للتوسع والوصول إلى الأمثلية العالمية، خاصة بالنسبة لمجموعات البيانات المعقدة والمتكررة والمشوشة.
تقدم الحوسبة الكمية مزايا محتملة للتنقل في مساحات البحث المعقدة هذه عبر خوارزميات مثل خوارزمية التقريب الكمي للتحسين (QAOA) والتلدين الكمي (Quantum Annealing). ومع ذلك، فإن النشر العملي مقيد بشدة بحدود أجهزة الحوسبة الكمية ذات المقياس المتوسط المشوبة بالضجيج (NISQ):
عدد الكيوبتات (Qubits): تتطلب رسوم الجينوم البيانية في الواقع آلاف المتغيرات، وهو ما يتجاوز بكثير أعداد الكيوبتات في الأجهزة الحالية (مثل 133 كيوبت في جهاز ibm_torino من IBM).
الاتصال والضجيج: تؤدي معدلات الخطأ العالية والاتصال المحدود بين الكيوبتات إلى دوائر عميقة تحتوي على الكثير من البوابات ثنائية الكيوبت، مما يؤدي إلى تراكم كبير للأخطاء والوصول إلى حلول دون المستوى الأمثل.
قابلية توسع الصياغة: غالبًا ما يؤدي ترميز مشكلات التجميع في صيغ التحسين الثنائي التربيعي غير المقيد (QUBO) إلى انفجار سريع في المتغيرات الثنائية، مما يجعل التضمين المباشر مستحيلاً.
التعامل مع القيود: يصعب فرض القيود البيولوجية الصارمة (مثل المسارات الصالحة، وعدم وجود تفرع) بشكل مثالي في صيغ QUBO غير المقيدة دون استخدام حدود جزائية ثقيلة قد تؤدي إلى تجميعات غير صالحة أو غير منطقية بيولوجيًا.
المنهجية يقترح المؤلفون خط معالجة لتجميع الجينوم بمساعدة كمية يدمج إطار عمل "التحسين عبر التفكيك التلقائي للهاميلتوني" (HADOF) مع صياغة QUBO لمشكلة التجميع. يستبدل سير العمل خطوة تحسين تخطيط الرسم البياني في المجمع الكلاسيكي Unicycler بنهج كمي.
البيانات والمعالجة المسبقة: تستخدم الدراسة بيانات تسلسل طويلة القراءة حقيقية من نوع Oxford Nanopole Technology (ONT) لجينوم Pseudomonas aeruginosa (حوالي 7.1 ميجابايت). تمت معالجة القراءات الخام باستخدام Minimap لتحديد التداخلات و Miniasm لبناء رسم بياني للسلاسل (string graph). هذا الرسم البياني، الذي يحتوي على 524 عقدة و2,313 حافة موجهة، يعمل كمدخل للعملية.
صياغة QUBO: بدلاً من صياغة مصفوفة التبديل (التي تتطلب N2 من المتغيرات)، يستخدم المؤلفون صياغة QUBO القائمة على الحواف. يقلل هذا من عدد المتغيرات إلى عدد الحواف (M=∣E∣). تتضمن دالة الهدف:
حد خطي سالب يكافئ تضمين الحافة.
حدود تربيعية جزائية لمنع التفرع (وجود عدة حواف داخلة أو خارجة من نفس الرأس). هذه الصياغة هي نموذج اتصال مخفف وليست ترميزًا صارمًا للمسار الهاميلتوني، مما يسمح بمزيد من المرونة على الأجهزة المشوبة بالضجيج.
تنفيذ HADOF: للتغلب على حدود الأجهزة، يتم تفكيك QUBO العالمي إلى مشكلات فرعية صغيرة (مجموعات فرعية من 5 متغيرات) باستخدام HADOF.
يقوم إطار العمل بحل هذه المشكلات الفرعية بشكل تكراري على أجهزة كمية (أو محاكيات) باستخدام QAOA (تحديدًا نسخة متروتريزية، QAOA-t).
يحافظ الإطار على متجه "اعتقاد" عالمي للاحتمالات الهامشية (P(xi)) المستمدة من التكرارات السابقة لتوفير السياق للمتغيرات غير النشطة في المشكلات الفرعية اللاحقة.
يسمح هذا النهج الاتحادي بحل مشكلة مكونة من 2,313 متغير باستخدام جهاز يحتوي على 133 كيوبت عبر تنفيذ 462 دائرة مكونة من خمسة كيوبتات ودائرة واحدة مكونة من ثلاثة كيوبتات.
المعالجة اللاحقة: نظرًا لأن الحلول الكمية تنتج توزيعات احتمالية للحلول بدلاً من مخرج حتمي واحد، فإن خطوة معالجة لاحقة كلاسيكية مطلوبة. يتم تطبيق خوارزمية البحث بالعمق أولاً (DFS) على حلول الرسم البياني التي تم أخذ عينات منها لحل الانتهاكات (مثل العقد ذات الدرجة > 1)، واستخراج أطول مسار خطي، وإجراء عملية التدوير (circularization). بعد ذلك، يتم صقل التسلسلات الناتية باستخدام Racon وتقييمها مقابل المرجع الأساسي Unicycler باستخدام QUAST.
المساهمات الرئيسية
تجميع كمي قابل للتوسع: يوضح العمل، حسب علم المؤلفين، أكبر رسم بياني لتجميع الجينوم تمت دراسته على أجهزة كمية حقيقية حتى الآن (7.1 ميجابايت)، متجاوزًا الدراسات السابقة التي اقتصرت على الجينومات الاصطناعية أو الفيروسية (أقل من 5 كيلوبايت).
تكامل HADOF: نجح في دمج HADOF لتجزئة مشكلة QUBO واسعة النطاق عبر مشكلات فرعية صغيرة، مما مكن من التنفيذ على أجهزة NISQ الحالية رغم قيود الكيوبت والاتصال.
التقييم التجريبي: تقدم الدراسة مقارنة صارمة بين المنهج الكمي (HADOF-QAOA)، والتلدين المحاكي (SA) الكلاسيكي، والمجمع الكلاسيكي المتطور Unicycler.
تحليل المقاييس: وضعت الدراسة إطارًا لتقييم العلاقة بين مقاييس تحسين QUBO المجردة (درجة الهدف) وجودة التجميع البيولوجي (كسر الجينوم، نسبة التكرار)، مما يسلط الضوء على أن تقليل درجة QUBO وحدها لا يضمن تجميعًا صالحًا بيولوجيًا.
النتائج قيمت الدراسة تجميع جينوم P. aeruginosa باستخدام ثلاثة حلول: التلدين المحاكي (SA)، و HADOF-QAOA على محاكي مثالي خالٍ من الضجيج، و HADOF-QAOA على جهاز ibm_torino الحقيقي.
جودة التجميع:
Unicycler (المرجع الأساسي): حقق تجميعًا دائريًا قدره 7,139,404 زوج قاعدة مع 100% من كسر الجينوم.
SA: أنتج تجميعات بنسبة 100% من كسر الجينوم ونسبة تكرار قدرها 1.0، مما يطابق تقريبًا المرجع الأساسي Unicycler.
HADOF-QAOA (المثالي): حقق 100% من كسر الجينوم مع نسبة تكرار 1.0، رغم أنه استخدم عقدًا أقل في المسار النهائي مقارنة بـ Unicycler.
HADOF-QAOA (الأجهزة الحقيقية): رغم الضجيج الكبير ودرجات QUBO الهدف الأعلى (بمتوسط -370 مقابل -992 لـ SA)، حقق أفضل حل تم أخذ عينة منه 99.348% من كسر الجينوم مع نسبة تكرار 1.0.
التحسين مقابل الجودة البيولوجية:
لوحظ انفصال بين درجة QUBO الخام وجودة التجميع البيولوجي. لم تضمن درجات QUBO المنخفضة كسور جينوم أعلى؛ حيث شملت الحلول ذات الدرجات المتشابهة طيفًا كاملاً من الاكتمال (من 15% إلى 100%).
الاحتفاظ بالعقد: لوحظ ارتباط خطي قوي بين عدد العقد التي تم الاحتفاظ بها بعد معالجة DFS وكسر الجينوم. الحلول التي احتفظت بـ ≥200 عقدة حققت باستمرار ما يقرب من 100% من كسور الجينوم.
الطول كمؤشر: أظهر الطول الإجمالي للجينوم المجمع علاقة خطية مع كسر الجينوم، حيث استقر عند الحجم المرجعي (7.1 ميجابايت). أظهرت التجميعات التي تجاوزت هذا الطول زيادة في نسب التكرار، مما يشير إلى تجميع مفرط ناتج عن التكرارات.
أداء الحلول (Solvers): أظهر SA أكثر توزيع للحلول قرب المرجع الأمثل. كانت نتائج الأجهزة الحقيقية أكثر تشتتًا، لكن استراتيجية التفكيك سمحت باستخراج حل واحد عالي الجودة على الأقل (99.35% كسر) من التوزيع المشوب بالضجيج.
الأهمية والادعاءات يزعم البحث أن هذا العمل يمثل خطوة كبيرة نحو المعلوماتية الحيوية الكمية العملية من خلال إثبات أن أطر التفكيك الهجينة يمكن أن تسد الفجوة بين التحسين الكمي النظري وقيود الأجهزة الحالية.
إثبات المفهوم: يثبت التجميع الناجح لجينوم 7.1 ميجابايت على أجهزة حقيقية أن عمليات إعادة البناء الصالحة بيولوجيًا ممكنة رغم قيود NISQ، بشرما تم استخدام استراتيجيات التفكيك والمعالجة اللاحقة.
رؤية منهجية: تسلط الدراسة الضوء على أنه بالنسبة لتجميع الجينوم، فإن الحل الرياضي "الأفضل" (أقل طاقة QUBO) ليس بالضريد هو الأكثر دقة بيولوجيًا. بدلاً من ذلك، تعمل الأنماط الهيكلية في توزيع الحلول (مثل أعداد الاحتفاظ بالعقد وأطوال المسارات) كمؤشرات أفضل لجودة التجميع.
الاتجاهات المستقبلية: يؤكد المؤلفون أنه بينما الأجهزة الحالية محدودة، فإن إصدار تنفيذ HADOF مفتوح المصدر وخط أنابيب التقييم يوفر أساسًا للأبحاث المستقبلية. ويقترحون أن التحسينات المستقبلية يجب أن تركز على تصفية المرشحين الخالية من المراجع بناءً على الأنماط الهيكلية الملاحظة، وإعادة صياغة عقوبات QUBO لتشمل الطبوغرافيا المحلية، وتطوير سياسات تجميع خاصة برسوم الجينوم البيانية.
يخلص العمل إلى أنه بينما لم يتحقق التفوق الكمي من حيث السرعة أو الأمثلية المطلقة بعد، فإن القدرة على توليد تجميعات عالية الجودة وصالحة بيولوجيًا على أجهزة حقيقية باستخدام أطر عمل هجينة تثبت صحة إمكانات الأساليب المدعومة كميًا لأعباء العمل العلمية المعقدة.