محاذاة التسلسل
في المعلوماتية الحيوية ، يُعدّ محاذاة التسلسل طريقةً لترتيب تسلسلات الحمض النووي (DNA) أو الحمض النووي الريبوزي (RNA ) أو البروتين لتحديد مناطق التشابه التي قد تكون ناتجةً عن علاقات وظيفية أو بنيوية أو تطورية بين التسلسلات. [ 1 ] تُمثَّل التسلسلات المحاذية لبقايا النيوكليوتيدات أو الأحماض الأمينية عادةً كصفوف ضمن مصفوفة . تُضاف فجوات بين البقايا بحيث تتم محاذاة الأحرف المتطابقة أو المتشابهة في أعمدة متتالية. تُستخدم محاذاة التسلسلات أيضًا لتسلسلات غير بيولوجية، مثل حساب تكلفة المسافة بين سلاسل نصية في لغة طبيعية ، أو لعرض البيانات المالية.

تفسير
إذا اشتركت سلسلتان في محاذاة واحدة في سلف مشترك، يُمكن تفسير عدم التطابق على أنه طفرات نقطية ، والفجوات على أنها طفرات إدخال أو حذف (أي طفرات إدخال أو حذف) حدثت في إحدى السلالتين أو كلتيهما منذ تباعدهما. في محاذاة تسلسلات البروتينات، يُمكن تفسير درجة التشابه بين الأحماض الأمينية التي تشغل موقعًا معينًا في التسلسل على أنها مقياس تقريبي لمدى حفظ منطقة أو نمط تسلسلي معين بين السلالات. يشير غياب الاستبدالات، أو وجود استبدالات محافظة للغاية فقط (أي استبدال الأحماض الأمينية التي تمتلك سلاسلها الجانبية خصائص كيميائية حيوية متشابهة) في منطقة معينة من التسلسل، إلى أن هذه المنطقة ذات أهمية هيكلية أو وظيفية. على الرغم من أن قواعد النيوكليوتيدات في الحمض النووي DNA والحمض النووي RNA أكثر تشابهًا من الأحماض الأمينية، إلا أن حفظ أزواج القواعد قد يشير إلى دور وظيفي أو هيكلي مماثل.
أساليب المحاذاة
يمكن محاذاة التسلسلات القصيرة جدًا أو المتشابهة جدًا يدويًا. مع ذلك، تتطلب معظم المشكلات المهمة محاذاة تسلسلات طويلة أو شديدة التباين أو كثيرة جدًا، والتي لا يمكن محاذاتها بالجهد البشري وحده. وقد طُوّرت خوارزميات متنوعة لإنتاج محاذاة تسلسلات عالية الجودة، وأحيانًا لتعديل النتائج النهائية لتعكس أنماطًا يصعب تمثيلها خوارزميًا (خاصةً في حالة تسلسلات النيوكليوتيدات). تنقسم الأساليب الحسابية لمحاذاة التسلسلات عمومًا إلى فئتين: المحاذاة الشاملة والمحاذاة المحلية . تُعد المحاذاة الشاملة شكلًا من أشكال التحسين الشامل الذي "يُجبر" المحاذاة على تغطية الطول الكامل لجميع تسلسلات الاستعلام. في المقابل، تُحدد المحاذاة المحلية مناطق التشابه داخل التسلسلات الطويلة التي غالبًا ما تكون متباعدة على نطاق واسع. غالبًا ما تكون المحاذاة المحلية مفضلة، ولكن قد يكون حسابها أكثر صعوبة بسبب التحدي الإضافي المتمثل في تحديد مناطق التشابه. [ 4 ] طُبّقت مجموعة متنوعة من الخوارزميات الحسابية على مشكلة محاذاة التسلسلات. وتشمل هذه الطرق البطيئة ولكنها صحيحة رسميًا مثل البرمجة الديناميكية . وتشمل هذه أيضًا الخوارزميات الفعالة والاستدلالية أو الأساليب الاحتمالية المصممة للبحث في قواعد البيانات واسعة النطاق، والتي لا تضمن العثور على أفضل التطابقات.
التمثيلات
تُعرض عمليات المحاذاة عادةً بيانيًا ونصيًا. في معظم تمثيلات محاذاة التسلسلات، تُكتب التسلسلات في صفوف مُرتبة بحيث تظهر البقايا المُحاذية في أعمدة متتالية. في النصوص، تُشار إلى الأعمدة المُحاذية التي تحتوي على أحرف متطابقة أو متشابهة بنظام من رموز الحفظ. كما في الصورة أعلاه، يُستخدم رمز النجمة أو الخط العمودي للدلالة على التطابق بين عمودين؛ ومن الرموز الأخرى الأقل شيوعًا النقطتان الرأسيتان للاستبدالات المحافظة والنقطة للاستبدالات شبه المحافظة. تستخدم العديد من برامج عرض التسلسلات الألوان لعرض معلومات حول خصائص عناصر التسلسل الفردية؛ في تسلسلات الحمض النووي DNA والحمض النووي الريبوزي RNA، يُعادل هذا تخصيص لون لكل نيوكليوتيد. في محاذاة البروتينات، كما في الصورة أعلاه، يُستخدم اللون غالبًا للإشارة إلى خصائص الأحماض الأمينية للمساعدة في تقييم حفظ استبدال حمض أميني مُعين. بالنسبة للتسلسلات المتعددة، غالبًا ما يكون الصف الأخير في كل عمود هو التسلسل التوافقي المُحدد بواسطة المحاذاة. كما يتم تمثيل التسلسل المتفق عليه غالبًا في شكل رسومي باستخدام شعار تسلسلي حيث يتناسب حجم كل حرف من النيوكليوتيدات أو الأحماض الأمينية مع درجة حفظه. [ 5 ]
يمكن تخزين محاذاة التسلسلات في مجموعة واسعة من تنسيقات الملفات النصية، وقد طُوِّر العديد منها في الأصل بالتزامن مع برنامج أو تطبيق محاذاة مُحدد. تسمح معظم الأدوات المستندة إلى الويب بعدد محدود من تنسيقات الإدخال والإخراج، مثل تنسيق FASTA وتنسيق GenBank ، كما أن الإخراج ليس سهل التعديل. تتوفر العديد من برامج التحويل التي توفر واجهات رسومية و/أو واجهات سطر أوامر، مثل READSEQ [ 6 ] و EMBOSS . كما توجد العديد من حزم البرمجة التي توفر وظيفة التحويل هذه، مثل BioPython و BioRuby و BioPerl . تستخدم ملفات SAM/BAM تنسيق سلسلة CIGAR (تقرير محاذاة الفجوات المُدمج) لتمثيل محاذاة تسلسل مع مرجع من خلال ترميز سلسلة من الأحداث (مثل التطابق/عدم التطابق، والإضافات، والحذف). [ 7 ]
تنسيق السيجار
المرجع : GTCGTAGAATA القراءة : CACGTAG—TA CIGAR: 2S5M2D2M حيث: 2S = قصّان ناعمان (قد يكونان عدم تطابق، أو قراءة أطول من التسلسل المطابق) 5M = 5 تطابقات أو عدم تطابقات 2D = حذفان 2M = تطابقان أو عدم تطابقان
لم يميز تنسيق CIGAR الأصلي من برنامج المحاذاة exonerate بين حالات عدم التطابق أو التطابقات التي تحتوي على الحرف M.
تحدد وثيقة مواصفات SAMv1 رموز CIGAR الأحدث. في معظم الحالات، يُفضل استخدام الرمزين '=' و'X' للدلالة على التطابقات أو عدم التطابقات بدلاً من الرمز 'M' القديم، الذي يُعد غامضاً.
| رمز السيجار | عدد صحيح BAM | وصف | يستهلك الاستعلام | مرجع المستهلك |
|---|---|---|---|---|
| م | 0 | تطابق المحاذاة (يمكن أن يكون تطابقًا أو عدم تطابق في التسلسل) | نعم | نعم |
| أنا | 1 | إدراج في المرجع | نعم | لا |
| د | 2 | حذف من المرجع | لا | نعم |
| شمال | 3 | المنطقة المتجاوزة من المرجع | لا | نعم |
| S | 4 | القطع الناعم (التسلسلات المقطوعة الموجودة في التسلسل) | نعم | لا |
| ح | 5 | القطع الحاد (التسلسلات المقطوعة غير موجودة في التسلسل) | لا | لا |
| P | 6 | الحشو (حذف صامت من المرجع المحشو) | لا | لا |
| = | 7 | تطابق التسلسل | نعم | نعم |
| X | 8 | عدم تطابق التسلسل | نعم | نعم |
- يشير "يستهلك الاستعلام" و "يستهلك المرجع" إلى ما إذا كانت عملية CIGAR تتسبب في أن يتحرك المحاذاة على طول تسلسل الاستعلام وتسلسل المرجع على التوالي.
- لا يمكن أن يكون H موجودًا إلا كعملية أولى و/أو أخيرة.
- لا يمكن أن تحتوي S إلا على عمليات H بينها وبين نهايات سلسلة CIGAR.
- في عملية محاذاة الحمض النووي الريبوزي الرسول مع الجينوم، تمثل العملية N إنترونًا. أما في أنواع المحاذاة الأخرى، فإن تفسير N غير محدد.
- يجب أن يساوي مجموع أطوال عمليات M/I/S/=/X طول التسلسل SEQ
التوافقات العالمية والمحلية
تُعدّ عمليات المحاذاة الشاملة، التي تسعى إلى محاذاة كل حمض أميني في كل تسلسل، مفيدةً للغاية عندما تكون التسلسلات في مجموعة الاستعلام متشابهة ومتساوية الحجم تقريبًا. (لا يعني هذا أن عمليات المحاذاة الشاملة لا يمكن أن تبدأ أو تنتهي في فجوات). ومن تقنيات المحاذاة الشاملة العامة خوارزمية نيدلمان-وونش ، التي تعتمد على البرمجة الديناميكية. أما عمليات المحاذاة المحلية فهي أكثر فائدة للتسلسلات غير المتشابهة التي يُشتبه في احتوائها على مناطق تشابه أو أنماط تسلسلية متشابهة ضمن سياق تسلسلها الأوسع. وتُعدّ خوارزمية سميث-واترمان طريقة عامة للمحاذاة المحلية تعتمد على نفس مخطط البرمجة الديناميكية، ولكن مع خيارات إضافية للبدء والانتهاء في أي مكان. [ 4 ]
تبحث الطرق الهجينة، المعروفة بالطرق شبه العالمية أو "العالمية-المحلية" (اختصارًا لـ " العالمية -المحلية " )، عن أفضل محاذاة جزئية ممكنة لتسلسلين (بمعنى آخر، يتم محاذاة أحد طرفي التسلسل أو كليهما، وأحد طرفيه أو كليهما). ويكون هذا مفيدًا بشكل خاص عندما يتداخل الجزء السفلي من أحد التسلسلين مع الجزء العلوي من التسلسل الآخر. في هذه الحالة، لا تُعد المحاذاة العالمية أو المحلية مناسبة تمامًا: فالمحاذاة العالمية تحاول إجبار المحاذاة على الامتداد إلى ما وراء منطقة التداخل، بينما قد لا تغطي المحاذاة المحلية منطقة التداخل بالكامل. [ 8 ] وهناك حالة أخرى تكون فيها المحاذاة شبه العالمية مفيدة، وهي عندما يكون أحد التسلسلين قصيرًا (مثل تسلسل جين) والآخر طويلًا جدًا (مثل تسلسل كروموسوم). في هذه الحالة، يجب محاذاة التسلسل القصير عالميًا (بشكل كامل)، بينما يُراد فقط محاذاة محلية (جزئية) للتسلسل الطويل.
يُشكّل التوسع السريع في البيانات الجينية تحديًا لسرعة خوارزميات محاذاة تسلسل الحمض النووي الحالية. وتتطلب الحاجة الماسة إلى طريقة فعّالة ودقيقة لاكتشاف متغيرات الحمض النووي مناهج مبتكرة للمعالجة المتوازية في الوقت الفعلي. وقد اقتُرحت أساليب الحوسبة الضوئية كبدائل واعدة للتطبيقات الكهربائية الحالية، إلا أن مدى قابليتها للتطبيق لا يزال قيد الاختبار..
محاذاة ثنائية
تُستخدم طرق محاذاة التسلسلات الثنائية لإيجاد أفضل تطابق بين تسلسلين متطابقين (محلي أو شامل). لا يمكن استخدام المحاذاة الثنائية إلا بين تسلسلين في كل مرة، ولكنها فعّالة من حيث الحساب، وغالبًا ما تُستخدم في الطرق التي لا تتطلب دقة فائقة (مثل البحث في قاعدة بيانات عن تسلسلات ذات تشابه عالٍ مع تسلسل الاستعلام). الطرق الثلاث الرئيسية لإنتاج المحاذاة الثنائية هي: طريقة المصفوفة النقطية، والبرمجة الديناميكية، وطريقة الكلمات؛ [ 1 ] ومع ذلك، يمكن لتقنيات محاذاة التسلسلات المتعددة أيضًا محاذاة أزواج من التسلسلات. على الرغم من أن لكل طريقة نقاط قوتها وضعفها، إلا أن جميع الطرق الثنائية الثلاث تواجه صعوبة مع التسلسلات المتكررة بكثرة ذات المحتوى المعلوماتي المنخفض ، خاصةً عندما يختلف عدد التكرارات بين التسلسلين المراد محاذاتهما.
أقصى تطابق فريد
إحدى طرق قياس فائدة محاذاة ثنائية معينة هي " التطابق الفريد الأقصى " (MUM)، أو أطول تسلسل فرعي يظهر في كلا تسلسلي الاستعلام. تعكس تسلسلات MUM الأطول عادةً قرابة أوثق [ 9 ] في محاذاة التسلسلات المتعددة للجينومات في علم الأحياء الحاسوبي . يُعد تحديد MUMs وغيرها من نقاط الارتكاز المحتملة الخطوة الأولى في أنظمة المحاذاة الأكبر مثل MUMmer . نقاط الارتكاز هي المناطق بين جينومين حيث يكونان متشابهين للغاية. لفهم ماهية MUM، يمكننا تحليل كل كلمة في الاختصار. تعني كلمة "تطابق" أن السلسلة الفرعية تظهر في كلا التسلسلين المراد محاذاتهما. تعني كلمة "فريد" أن السلسلة الفرعية تظهر مرة واحدة فقط في كل تسلسل. وأخيرًا، تعني كلمة "أقصى" أن السلسلة الفرعية ليست جزءًا من سلسلة أخرى أطول تستوفي الشرطين السابقين. الفكرة وراء ذلك هي أن التسلسلات الطويلة التي تتطابق تمامًا وتظهر مرة واحدة فقط في كل جينوم هي بالتأكيد جزء من المحاذاة الشاملة.
وبشكل أدق:
"بالنظر إلى جينومين A و B، فإن السلسلة الفرعية ذات التطابق الفريد الأقصى (MUM) هي سلسلة فرعية مشتركة بين A و B بطول أطول من الحد الأدنى المحدد d (افتراضيًا d = 20) بحيث
- إنها قصوى، أي لا يمكن تمديدها من أي من طرفيها دون حدوث عدم تطابق؛ و
- إنها فريدة في كلا التسلسلين" [ 10 ]
طرق المصفوفة النقطية
يُعدّ أسلوب المصفوفة النقطية، الذي يُنتج ضمنيًا مجموعة من المحاذاة لمناطق التسلسل الفردية، أسلوبًا نوعيًا وبسيطًا من الناحية المفاهيمية، ولكنه يستغرق وقتًا طويلًا عند تحليله على نطاق واسع. في غياب التشويش، يُمكن بسهولة تحديد بعض سمات التسلسل بصريًا - مثل الإضافات، والحذف، والتكرارات، أو التكرارات المعكوسة - من مخطط المصفوفة النقطية. لإنشاء مخطط المصفوفة النقطية ، يُكتب التسلسلان على طول الصف العلوي والعمود الأيسر لمصفوفة ثنائية الأبعاد، وتُوضع نقطة عند أي نقطة تتطابق فيها الأحرف في الأعمدة المناسبة - وهذا هو مخطط التكرار النموذجي . تُغيّر بعض التطبيقات حجم النقطة أو شدتها اعتمادًا على درجة تشابه الحرفين، لاستيعاب الاستبدالات المحافظة. ستظهر مخططات النقاط للتسلسلات شديدة الترابط كخط واحد على طول القطر الرئيسي للمصفوفة .
تشمل مشاكل مخططات النقاط كتقنية لعرض المعلومات: التشويش، وعدم الوضوح، وعدم سهولة الاستخدام، وصعوبة استخراج إحصائيات ملخص التطابق ومواقع التطابق على التسلسلين. كما يوجد مساحة مهدرة كبيرة حيث تتكرر بيانات التطابق بشكل طبيعي عبر القطر، وتشغل المساحة الفارغة أو التشويش معظم مساحة المخطط، وأخيرًا، تقتصر مخططات النقاط على تسلسلين. لا تنطبق أي من هذه القيود على مخططات محاذاة Miropeats، ولكنها تعاني من عيوبها الخاصة.
يمكن استخدام مخططات النقاط لتقييم التكرار في تسلسل واحد. يُرسم التسلسل مقابل نفسه، وتظهر المناطق التي تتشارك في أوجه تشابه كبيرة كخطوط خارج القطر الرئيسي. يحدث هذا التأثير عندما يتكون البروتين من نطاقات هيكلية متعددة متشابهة .
البرمجة الديناميكية
يمكن تطبيق تقنية البرمجة الديناميكية لإنتاج محاذاة شاملة باستخدام خوارزمية نيدلمان-وونش ، ومحاذاة محلية باستخدام خوارزمية سميث-واترمان . في الاستخدام المعتاد، تستخدم محاذاة البروتينات مصفوفة استبدال لتعيين درجات لتطابق أو عدم تطابق الأحماض الأمينية، وعقوبة فجوة لمطابقة حمض أميني في تسلسل ما مع فجوة في تسلسل آخر. قد تستخدم محاذاة الحمض النووي DNA والحمض النووي الريبي RNA مصفوفة تسجيل، ولكن عمليًا غالبًا ما يتم تعيين درجة تطابق موجبة، ودرجة عدم تطابق سالبة، وعقوبة فجوة سالبة. (في البرمجة الديناميكية القياسية، تكون درجة كل موضع حمض أميني مستقلة عن هوية جيرانه، وبالتالي لا تُؤخذ تأثيرات تكديس القواعد في الاعتبار. ومع ذلك، من الممكن مراعاة هذه التأثيرات عن طريق تعديل الخوارزمية). من الامتدادات الشائعة لتكاليف الفجوة الخطية القياسية تكاليف الفجوة الأفينية. هنا، يتم تطبيق عقوبتين مختلفتين لفتح فجوة ولتمديدها. عادةً ما تكون الأولى أكبر بكثير من الثانية، على سبيل المثال -10 لفتح الفجوة و-2 لتمديدها. ينتج عن ذلك عدد أقل من الفجوات في المحاذاة، مع الحفاظ على بقايا الأحماض الأمينية والفجوات معًا، وهي سمات أكثر تمثيلًا للتسلسلات البيولوجية. تُطبّق خوارزمية غوتو تكاليف الفجوات الخطية باستخدام ثلاث مصفوفات. [ 11 ] [ 12 ]
يمكن أن تكون البرمجة الديناميكية مفيدة في محاذاة تسلسلات النيوكليوتيدات مع تسلسلات البروتينات، وهي مهمة معقدة بسبب الحاجة إلى مراعاة طفرات الإزاحة الإطارية (عادةً ما تكون إدخالات أو حذف). تُنتج طريقة البحث عن الإطار سلسلة من المحاذيات الثنائية العالمية أو المحلية بين تسلسل نيوكليوتيدات مُستعلم عنه ومجموعة بحث من تسلسلات البروتينات، أو العكس. إن قدرتها على تقييم إزاحات الإطار التي تُزاح بعدد عشوائي من النيوكليوتيدات تجعل هذه الطريقة مفيدة للتسلسلات التي تحتوي على أعداد كبيرة من الإدخالات والحذف، والتي قد يصعب محاذاتها باستخدام طرق استدلالية أكثر كفاءة. عمليًا، تتطلب هذه الطريقة قدرة حاسوبية كبيرة أو نظامًا مُصممًا خصيصًا للبرمجة الديناميكية. توفر مجموعتا BLAST و EMBOSS أدوات أساسية لإنشاء محاذيات مُترجمة (على الرغم من أن بعض هذه الأساليب تستفيد من الآثار الجانبية لقدرات البحث عن التسلسل في هذه الأدوات). تتوفر طرق أكثر عمومية من خلال برامج مفتوحة المصدر مثل GeneWise .
تضمن طريقة البرمجة الديناميكية إيجاد محاذاة مثالية بناءً على دالة تقييم محددة؛ إلا أن تحديد دالة تقييم جيدة غالبًا ما يكون مسألة تجريبية وليست نظرية. ورغم إمكانية توسيع نطاق البرمجة الديناميكية لتشمل أكثر من سلسلتين، إلا أنها بطيئة للغاية عند التعامل مع أعداد كبيرة من السلاسل أو السلاسل الطويلة جدًا.
أساليب معالجة النصوص
تُعرف طرق الكلمات، أو طرق k -tuple، بأنها طرق استدلالية لا تضمن إيجاد حل محاذاة مثالي، ولكنها أكثر كفاءة بشكل ملحوظ من البرمجة الديناميكية. تُعد هذه الطرق مفيدة بشكل خاص في عمليات البحث في قواعد البيانات واسعة النطاق، حيث يُفترض أن نسبة كبيرة من التسلسلات المرشحة لن يكون لها تطابق يُذكر مع تسلسل الاستعلام. تشتهر طرق الكلمات بتطبيقها في أدوات البحث في قواعد البيانات FASTA وعائلة BLAST . [ 1 ] تُحدد طرق الكلمات سلسلة من التسلسلات الفرعية القصيرة غير المتداخلة ("الكلمات") في تسلسل الاستعلام، والتي تُطابق بعد ذلك مع تسلسلات قاعدة البيانات المرشحة. يُطرح الموضع النسبي للكلمة في التسلسلين المُقارنين للحصول على إزاحة؛ وهذا يُشير إلى منطقة محاذاة إذا أنتجت كلمات متعددة ومختلفة نفس الإزاحة. فقط في حال اكتشاف هذه المنطقة، تُطبق هذه الطرق معايير محاذاة أكثر حساسية؛ وبالتالي، يتم التخلص من العديد من المقارنات غير الضرورية مع التسلسلات التي لا يوجد تشابه يُذكر.
في طريقة FASTA، يُحدد المستخدم قيمة k لاستخدامها كطول الكلمة للبحث في قاعدة البيانات. تكون هذه الطريقة أبطأ ولكنها أكثر حساسية عند قيم k المنخفضة ، وهي القيم المُفضلة أيضًا لعمليات البحث التي تتضمن تسلسل استعلام قصير جدًا. توفر عائلة BLAST من طرق البحث عددًا من الخوارزميات المُحسّنة لأنواع مُحددة من الاستعلامات، مثل البحث عن تطابقات تسلسلية متباعدة الصلة. طُوّرت BLAST لتوفير بديل أسرع لـ FASTA دون التضحية بالكثير من الدقة؛ فمثل FASTA، تستخدم BLAST بحثًا عن الكلمات بطول k ، ولكنها تُقيّم فقط تطابقات الكلمات الأكثر أهمية، بدلًا من كل تطابق كلمة كما تفعل FASTA. تستخدم معظم تطبيقات BLAST طول كلمة افتراضيًا ثابتًا مُحسّنًا لنوع الاستعلام وقاعدة البيانات، ولا يتغير إلا في ظروف خاصة، مثل البحث باستخدام تسلسلات استعلام متكررة أو قصيرة جدًا. يمكن العثور على التطبيقات عبر عدد من البوابات الإلكترونية، مثل EMBL FASTA و NCBI BLAST .
محاذاة التسلسل المتعدد

يُعدّ محاذاة التسلسلات المتعددة امتدادًا للمحاذاة الثنائية، حيث يشمل أكثر من تسلسلين في آنٍ واحد. وتسعى طرق المحاذاة المتعددة إلى محاذاة جميع التسلسلات في مجموعة استعلام مُعطاة. تُستخدم المحاذاة المتعددة غالبًا لتحديد مناطق التسلسل المحفوظة عبر مجموعة من التسلسلات التي يُفترض أنها مرتبطة تطوريًا. ويمكن استخدام هذه الأنماط التسلسلية المحفوظة جنبًا إلى جنب مع المعلومات الهيكلية والآلية لتحديد المواقع النشطة التحفيزية للإنزيمات . كما تُستخدم المحاذاة للمساعدة في تحديد العلاقات التطورية من خلال بناء الأشجار التطورية . تُعدّ محاذاة التسلسلات المتعددة عملية حسابية معقدة، وتؤدي معظم صياغات هذه المشكلة إلى مسائل تحسين توافقي كاملة من نوع NP . [ 13 ] [ 14 ] ومع ذلك، فقد أدت فائدة هذه المحاذاة في المعلوماتية الحيوية إلى تطوير مجموعة متنوعة من الطرق المناسبة لمحاذاة ثلاثة تسلسلات أو أكثر.
البرمجة الديناميكية
تُعدّ تقنية البرمجة الديناميكية قابلة للتطبيق نظريًا على أي عدد من التسلسلات؛ إلا أنها، نظرًا لتكلفتها الحسابية العالية من حيث الوقت والذاكرة ، نادرًا ما تُستخدم لأكثر من ثلاثة أو أربعة تسلسلات في أبسط صورها. تتطلب هذه الطريقة إنشاء مكافئ ذي n بُعد لمصفوفة التسلسل المُشكّلة من تسلسلين، حيث n هو عدد التسلسلات في الاستعلام. تُطبّق البرمجة الديناميكية القياسية أولًا على جميع أزواج تسلسلات الاستعلام، ثم تُملأ "مساحة المحاذاة" من خلال النظر في التطابقات أو الفجوات المحتملة في المواضع الوسيطة، لتُنشئ في النهاية محاذاة بين كل محاذاة لتسلسلين. على الرغم من أن هذه التقنية مُكلفة حسابيًا، إلا أن ضمانها لحل أمثل عالمي يُعدّ مفيدًا في الحالات التي لا تتطلب سوى محاذاة عدد قليل من التسلسلات بدقة. وقد طُبّقت إحدى طرق تقليل المتطلبات الحسابية للبرمجة الديناميكية، والتي تعتمد على دالة الهدف "مجموع الأزواج"، في حزمة برامج MSA [ تم حذف الرابط ] . [ 15 ]
الأساليب التدريجية
تُنشئ الطرق التدريجية أو الهرمية أو الشجرية محاذاة تسلسلية متعددة من خلال محاذاة التسلسلات الأكثر تشابهًا أولًا، ثم إضافة التسلسلات أو المجموعات الأقل تشابهًا تدريجيًا إلى المحاذاة حتى يتم دمج مجموعة الاستعلام بأكملها في الحل. تعتمد الشجرة الأولية التي تصف ترابط التسلسلات على مقارنات ثنائية قد تتضمن طرق محاذاة ثنائية استدلالية مشابهة لـ FASTA . تعتمد نتائج المحاذاة التدريجية على اختيار التسلسلات "الأكثر ترابطًا"، وبالتالي قد تتأثر بعدم دقة المقارنات الثنائية الأولية. بالإضافة إلى ذلك، تُرجّح معظم طرق المحاذاة التدريجية المتعددة التسلسلات في مجموعة الاستعلام وفقًا لترابطها، مما يقلل من احتمالية اختيار تسلسلات أولية غير مناسبة، وبالتالي يُحسّن دقة المحاذاة.
تُستخدم العديد من صيغ تطبيق Clustal التدريجي [ 16 ] [ 17 ] [ 18 ] في محاذاة التسلسلات المتعددة، وبناء الأشجار التطورية، وكبيانات إدخال للتنبؤ ببنية البروتين . ويُعرف أحد المتغيرات الأبطأ ولكن الأكثر دقة من الطريقة التدريجية باسم T-Coffee . [ 19 ]
الأساليب التكرارية
تسعى الطرق التكرارية إلى التغلب على الاعتماد الكبير على دقة عمليات المحاذاة الثنائية الأولية، والتي تُعدّ نقطة ضعف الطرق التقدمية. تعمل هذه الطرق على تحسين دالة الهدف بناءً على طريقة مُختارة لتقييم المحاذاة، وذلك بتحديد محاذاة شاملة أولية ثم إعادة محاذاة مجموعات فرعية من التسلسل. بعد ذلك، تُحاذى المجموعات الفرعية المُعاد محاذاتها بدورها لإنتاج محاذاة التسلسل المتعددة للتكرار التالي. تُستعرض طرق مختلفة لاختيار مجموعات التسلسل الفرعية ودالة الهدف في المرجع [ 20 ] .
إيجاد الزخارف
تعتمد تقنية البحث عن الأنماط، المعروفة أيضًا بتحليل الملفات التعريفية، على إنشاء محاذاة تسلسلية متعددة شاملة، بهدف تحديد أنماط التسلسل المحفوظة القصيرة بين التسلسلات في مجموعة الاستعلام. ويتم ذلك عادةً بإنشاء محاذاة تسلسلية متعددة شاملة أولًا، ثم عزل المناطق المحفوظة بدرجة عالية واستخدامها لإنشاء مجموعة من مصفوفات الملفات التعريفية. تُرتّب مصفوفة الملف التعريفي لكل منطقة محفوظة على غرار مصفوفة التسجيل، ولكن تُستمدّ قيم التكرار لكل حمض أميني أو نيوكليوتيد في كل موضع من توزيع الأحرف في المنطقة المحفوظة، بدلًا من توزيع تجريبي أكثر عمومية. تُستخدم مصفوفات الملفات التعريفية بعد ذلك للبحث في تسلسلات أخرى عن وجود النمط الذي تُميّزه. في الحالات التي تحتوي فيها مجموعة البيانات الأصلية على عدد قليل من التسلسلات، أو على تسلسلات شديدة الصلة فقط، تُضاف قيم افتراضية لتطبيع توزيعات الأحرف المُمثلة في النمط.
تقنيات مستوحاة من علوم الحاسوب

تم تطبيق مجموعة متنوعة من خوارزميات التحسين العامة الشائعة الاستخدام في علوم الحاسوب على مشكلة محاذاة التسلسلات المتعددة. استُخدمت نماذج ماركوف المخفية لإنتاج درجات احتمالية لمجموعة من محاذاة التسلسلات المتعددة المحتملة لمجموعة استعلام معينة؛ على الرغم من أن الطرق المبكرة القائمة على نماذج ماركوف المخفية لم تُحقق أداءً مُرضيًا، إلا أن التطبيقات اللاحقة وجدتها فعالة بشكل خاص في اكتشاف التسلسلات ذات الصلة البعيدة لأنها أقل عرضة للتشويش الناتج عن الاستبدالات المحافظة أو شبه المحافظة. [ 21 ] كما استُخدمت الخوارزميات الجينية ومحاكاة التلدين في تحسين درجات محاذاة التسلسلات المتعددة وفقًا لدالة تسجيل مثل طريقة مجموع الأزواج. يمكن الاطلاع على مزيد من التفاصيل وحزم البرامج في المقالة الرئيسية " محاذاة التسلسلات المتعددة" .
تم تطبيق تحويل بوروز-ويلر بنجاح على محاذاة القراءات القصيرة السريعة في أدوات شائعة مثل Bowtie و BWA . انظر فهرس FM .
المحاذاة الهيكلية
تستخدم عمليات المحاذاة البنيوية، التي عادةً ما تكون خاصة بتسلسلات البروتينات وأحيانًا الحمض النووي الريبي (RNA)، معلوماتٍ حول البنية الثانوية والثالثية لجزيء البروتين أو الحمض النووي الريبي للمساعدة في محاذاة التسلسلات. يمكن استخدام هذه الطرق لتسلسلين أو أكثر، وتُنتج عادةً محاذاةً محلية؛ ومع ذلك، نظرًا لاعتمادها على توافر المعلومات البنيوية، لا يمكن استخدامها إلا مع التسلسلات التي تُعرف بنيتها المقابلة (عادةً من خلال علم البلورات بالأشعة السينية أو مطيافية الرنين النووي المغناطيسي ). ولأن بنية كلٍ من البروتين والحمض النووي الريبي أكثر حفظًا تطوريًا من التسلسل، [ 22 ] فإن عمليات المحاذاة البنيوية قد تكون أكثر موثوقية بين التسلسلات المتباعدة جدًا والتي تباعدت بشكلٍ كبير بحيث لا يمكن لمقارنة التسلسلات الكشف عن تشابهها بشكلٍ موثوق.
تُستخدم المحاذاة البنيوية كمعيار ذهبي في تقييم المحاذاة المستخدمة في التنبؤ ببنية البروتين بناءً على التشابه [ 23 ] ، لأنها تُحاذي صراحةً مناطق تسلسل البروتين المتشابهة بنيويًا، بدلًا من الاعتماد كليًا على معلومات التسلسل. مع ذلك، من الواضح أنه لا يمكن استخدام المحاذاة البنيوية في التنبؤ بالبنية، لأن أحد التسلسلات على الأقل في مجموعة الاستعلامات هو التسلسل المستهدف المراد نمذجته، والذي لا تُعرف بنيته. وقد ثبت أنه عند معرفة المحاذاة البنيوية بين تسلسل مستهدف وتسلسل نموذجي، يمكن إنتاج نماذج عالية الدقة لتسلسل البروتين المستهدف؛ إلا أن العقبة الرئيسية في التنبؤ بالبنية بناءً على التشابه تكمن في إنتاج محاذاة دقيقة بنيويًا بالاعتماد على معلومات التسلسل فقط. [ 23 ]
دالي
طريقة DALI، أو محاذاة مصفوفة المسافة ، هي طريقة تعتمد على التجزئة لإنشاء محاذاة هيكلية بناءً على أنماط تشابه التلامس بين الببتيدات السداسية المتتالية في تسلسلات الاستعلام. [ 24 ] يمكنها توليد محاذاة ثنائية أو متعددة وتحديد الجيران الهيكليين لتسلسل الاستعلام في بنك بيانات البروتين (PDB). وقد استُخدمت لإنشاء قاعدة بيانات محاذاة FSSP الهيكلية (تصنيف الطيات بناءً على محاذاة بنية البروتينات، أو عائلات البروتينات المتشابهة هيكليًا). يمكن الوصول إلى خادم DALI عبر الإنترنت من خلال DALI ، بينما تقع قاعدة بيانات FSSP في قاعدة بيانات DALI .
SSAP
برنامج SSAP (برنامج محاذاة البنية التسلسلية) هو طريقة تعتمد على البرمجة الديناميكية لمحاذاة البنية، وتستخدم متجهات الذرات في فضاء البنية كنقاط مقارنة. وقد تم توسيعه منذ وصفه الأصلي ليشمل المحاذاة المتعددة والثنائية، [ 25 ] واستُخدم في بناء قاعدة بيانات CATH (الفئة، البنية، الطوبولوجيا، التماثل) الهرمية لتصنيف طيات البروتين. [ 26 ] يمكن الوصول إلى قاعدة بيانات CATH من خلال CATH Protein Structure Classification .
الامتداد التوافقي
تُنشئ طريقة التمديد التوافقي للمحاذاة الهيكلية محاذاة هيكلية ثنائية باستخدام الهندسة المحلية لمحاذاة أجزاء قصيرة من البروتينين قيد التحليل، ثم تُجمّع هذه الأجزاء في محاذاة أكبر. [ 27 ] استنادًا إلى مقاييس مثل متوسط الجذر التربيعي للمسافة الصلبة ، ومسافات الأحماض الأمينية، والبنية الثانوية المحلية، وخصائص البيئة المحيطة مثل كراهية الماء بين الأحماض الأمينية المجاورة ، تُنشأ محاذاة محلية تُسمى "أزواج الأجزاء المُحاذية" وتُستخدم لبناء مصفوفة تشابه تُمثل جميع المحاذاة الهيكلية الممكنة ضمن معايير قطع مُحددة مسبقًا. بعد ذلك، يُتتبع مسار من حالة بنية بروتينية إلى أخرى عبر المصفوفة عن طريق تمديد المحاذاة المتنامية جزءًا تلو الآخر. يُحدد المسار الأمثل محاذاة التمديد التوافقي. يوجد خادم ويب يُنفذ هذه الطريقة ويُوفر قاعدة بيانات للمحاذاة الثنائية للهياكل في بنك بيانات البروتين على موقع التمديد التوافقي .
التحليل الوراثي
يُعدّ علم الوراثة العرقي ومحاذاة التسلسل مجالين وثيقي الصلة نظرًا لضرورة تقييم القرابة بين التسلسلات. [ 28 ] يستخدم علم الوراثة العرقي محاذاة التسلسلات على نطاق واسع في بناء وتفسير الأشجار الوراثية العرقية ، والتي تُستخدم لتصنيف العلاقات التطورية بين الجينات المتماثلة الموجودة في جينومات الأنواع المتباينة. ترتبط درجة اختلاف التسلسلات في مجموعة الاستعلام نوعيًا بالمسافة التطورية بينها. وبشكل عام، تشير نسبة التطابق العالية في التسلسلات إلى أن التسلسلات المعنية تشترك في سلف مشترك حديث نسبيًا ، بينما تشير نسبة التطابق المنخفضة إلى أن التباعد أقدم. هذا التقريب، الذي يعكس فرضية " الساعة الجزيئية " التي تفترض إمكانية استخدام معدل ثابت تقريبًا للتغير التطوري لاستقراء الوقت المنقضي منذ تباعد جينين لأول مرة (أي وقت الاندماج )، يفترض أن تأثيرات الطفرة والانتقاء ثابتة عبر سلالات التسلسل. لذا، لا يأخذ هذا النموذج في الحسبان الاختلافات المحتملة بين الكائنات الحية أو الأنواع في معدلات إصلاح الحمض النووي ، أو الحفاظ الوظيفي المحتمل لمناطق محددة في التسلسل. (في حالة تسلسلات النيوكليوتيدات، تتجاهل فرضية الساعة الجزيئية في أبسط صورها أيضًا الاختلاف في معدلات القبول بين الطفرات الصامتة التي لا تُغير معنى كودون معين ، والطفرات الأخرى التي تؤدي إلى دمج حمض أميني مختلف في البروتين). تسمح الطرق الأكثر دقة إحصائيًا بتفاوت معدل التطور على كل فرع من فروع الشجرة التطورية، مما يُنتج تقديرات أفضل لأوقات اندماج الجينات.
تُنتج تقنيات المحاذاة المتعددة التدريجية شجرة تطورية بالضرورة، لأنها تُدمج التسلسلات في المحاذاة المتنامية وفقًا لترتيب القرابة. أما التقنيات الأخرى التي تُجمّع محاذاة التسلسلات المتعددة والأشجار التطورية، فتقوم أولًا بتقييم الأشجار وفرزها، ثم تحسب محاذاة التسلسلات المتعددة من الشجرة ذات أعلى درجة. وتعتمد الطرق الشائعة الاستخدام لبناء الأشجار التطورية بشكل أساسي على الاستدلال ، لأن مشكلة اختيار الشجرة المثلى، مثل مشكلة اختيار محاذاة التسلسلات المتعددة المثلى، تُعدّ من المسائل الصعبة حسابيًا (NP-hard) . [ 29 ]
تقييم الأهمية
تُعدّ محاذاة التسلسلات مفيدة في المعلوماتية الحيوية لتحديد تشابه التسلسلات، وإنتاج الأشجار التطورية، وتطوير نماذج التماثل لبنية البروتينات. مع ذلك، فإن الأهمية البيولوجية لمحاذاة التسلسلات ليست واضحة دائمًا. يُفترض غالبًا أن المحاذاة تعكس درجة من التغير التطوري بين التسلسلات المنحدرة من سلف مشترك؛ ومع ذلك، من الممكن نظريًا أن يحدث تطور تقاربي ينتج عنه تشابه ظاهري بين بروتينات غير مرتبطة تطوريًا ولكنها تؤدي وظائف متشابهة ولها بنى متشابهة.
في عمليات البحث في قواعد البيانات، مثل BLAST، تُستخدم الأساليب الإحصائية لتحديد احتمالية ظهور محاذاة معينة بين تسلسلات أو مناطق تسلسلية بالصدفة، وذلك بناءً على حجم قاعدة البيانات وتكوينها. وتختلف هذه الاحتمالية اختلافًا كبيرًا تبعًا لنطاق البحث. فعلى وجه الخصوص، تزداد احتمالية العثور على محاذاة معينة بالصدفة إذا كانت قاعدة البيانات تتكون فقط من تسلسلات من نفس الكائن الحي الذي ينتمي إليه تسلسل الاستعلام. كما أن التسلسلات المتكررة في قاعدة البيانات أو الاستعلام قد تُشوه نتائج البحث وتقييم الدلالة الإحصائية؛ لذا يقوم BLAST تلقائيًا بتصفية هذه التسلسلات المتكررة في الاستعلام لتجنب النتائج الظاهرية التي هي في الواقع نتائج غير دقيقة إحصائيًا.
تتوفر في الأدبيات العلمية طرق لتقدير الدلالة الإحصائية لمحاذاة التسلسلات ذات الفجوات. [ 28 ] [ 30 ] [ 31 ] [ 32 ] [ 33 ] [ 34 ] [ 35 ] [ 36 ]
تقييم المصداقية
تشير الدلالة الإحصائية إلى احتمالية ظهور محاذاة بجودة معينة بالصدفة، لكنها لا تُبين مدى تفوق محاذاة معينة على المحاذاة البديلة لنفس التسلسلات. تُشير مقاييس مصداقية المحاذاة إلى مدى تشابه أفضل المحاذاة المُسجلة لزوج معين من التسلسلات. تتوفر في الأدبيات العلمية طرق لتقدير مصداقية المحاذاة للتسلسلات ذات الفجوات. [ 37 ]
دوال التسجيل
يُعدّ اختيار دالة تقييم تعكس الملاحظات البيولوجية أو الإحصائية حول التسلسلات المعروفة أمرًا بالغ الأهمية لإنتاج محاذاة جيدة. غالبًا ما تُحاذى تسلسلات البروتين باستخدام مصفوفات الاستبدال التي تعكس احتمالات استبدال الأحرف. تُشفّر سلسلة من المصفوفات تُسمى مصفوفات PAM (مصفوفات الطفرات المقبولة نقطيًا، والتي عرّفتها مارغريت دايوف في الأصل ، ويُشار إليها أحيانًا باسم "مصفوفات دايوف") بشكل صريح تقريبات تطورية تتعلق بمعدلات واحتمالات طفرات الأحماض الأمينية المحددة. وهناك سلسلة أخرى شائعة من مصفوفات التقييم، تُعرف باسم BLOSUM (مصفوفة استبدال الكتل)، تُشفّر احتمالات الاستبدال المستمدة تجريبيًا. تُستخدم متغيرات من كلا النوعين من المصفوفات للكشف عن التسلسلات ذات مستويات التباين المختلفة، مما يسمح لمستخدمي BLAST أو FASTA بتقييد عمليات البحث إلى التطابقات الأكثر ترابطًا أو توسيع نطاقها للكشف عن التسلسلات الأكثر تباينًا. تُراعي عقوبات الفجوات إدخال فجوة - في النموذج التطوري، طفرة إدخال أو حذف - في كل من تسلسلات النيوكليوتيدات والبروتينات، وبالتالي ينبغي أن تتناسب قيم العقوبات مع المعدل المتوقع لمثل هذه الطفرات. وعليه، تعتمد جودة عمليات المحاذاة الناتجة على جودة دالة التقييم.
قد يكون من المفيد والمثمر للغاية تجربة نفس عملية المحاذاة عدة مرات مع اختيار قيم مختلفة لمصفوفة التقييم و/أو جزاءات الفجوات، ثم مقارنة النتائج. ويمكن غالباً تحديد المناطق التي يكون فيها الحل ضعيفاً أو غير فريد من خلال ملاحظة أي مناطق المحاذاة تكون ثابتة في مواجهة التغيرات في معايير المحاذاة.
استخدامات بيولوجية أخرى
يمكن محاذاة الحمض النووي الريبوزي المُتسلسل، مثل علامات التسلسل المُعبر عنها وجزيئات mRNA كاملة الطول، مع جينوم مُتسلسل لتحديد مواقع الجينات والحصول على معلومات حول التضفير البديل [ 38 ] وتعديل الحمض النووي الريبوزي [ 39 ] . تُعد محاذاة التسلسل جزءًا من تجميع الجينوم ، حيث تُحاذى التسلسلات للعثور على مناطق التداخل لتكوين قطع متجاورة (أجزاء طويلة من التسلسل) [ 40 ] . ومن الاستخدامات الأخرى تحليل تعدد أشكال النوكليوتيدات المفردة (SNP )، حيث تُحاذى التسلسلات من أفراد مختلفين للعثور على أزواج القواعد المفردة التي غالبًا ما تختلف في مجتمع ما [ 41 ] .
الاستخدامات غير البيولوجية
وجدت أساليب محاذاة التسلسلات البيولوجية تطبيقات في مجالات أخرى، أبرزها معالجة اللغات الطبيعية والعلوم الاجتماعية ، حيث يُشار عادةً إلى خوارزمية نيدلمان-وونش باسم المطابقة المثلى . [ 42 ] وقد استعارت التقنيات التي تُولّد مجموعة العناصر التي تُختار منها الكلمات في خوارزميات توليد اللغة الطبيعية تقنيات محاذاة التسلسلات المتعددة من المعلوماتية الحيوية لإنتاج نسخ لغوية من البراهين الرياضية المُولّدة حاسوبيًا . [ 43 ] وفي مجال اللغويات التاريخية والمقارنة ، استُخدمت محاذاة التسلسلات لأتمتة جزئية للطريقة المقارنة التي يستخدمها اللغويون تقليديًا لإعادة بناء اللغات. [ 44 ] كما طبّقت أبحاث الأعمال والتسويق تقنيات محاذاة التسلسلات المتعددة في تحليل سلاسل المشتريات عبر الزمن. [ 45 ]
برمجة
تتوفر قائمة أكثر شمولاً بالبرامج المتاحة، مصنفة حسب الخوارزمية ونوع المحاذاة، على موقع برامج محاذاة التسلسل . تشمل الأدوات البرمجية الشائعة المستخدمة في مهام محاذاة التسلسل العامة ClustalW2 [ 46 ] وT-coffee [ 47 ] للمحاذاة، وBLAST [ 48 ] وFASTA3x [ 49 ] للبحث في قواعد البيانات. كما تتوفر أدوات تجارية مثل DNASTAR وLasergene و Geneious و PatternHunter . تُدرج الأدوات المصنفة على أنها تُجري محاذاة التسلسل في سجل bio.tools .
يمكن مقارنة خوارزميات وبرامج المحاذاة مباشرةً باستخدام مجموعة معيارية من محاذاة التسلسلات المتعددة المرجعية ، تُعرف باسم BAliBASE. [ 50 ] تتكون مجموعة البيانات من محاذاة هيكلية، والتي يمكن اعتبارها معيارًا تُقارن به الطرق القائمة على التسلسل فقط. جُدول الأداء النسبي للعديد من طرق المحاذاة الشائعة في مشاكل المحاذاة المتكررة، ونُشرت نتائج مختارة على الإنترنت في BAliBASE. [ 51 ] [ 52 ] يمكن حساب قائمة شاملة بنتائج BAliBASE للعديد من أدوات المحاذاة المختلفة (12 أداة حاليًا) ضمن منصة STRAP لتحليل البروتينات. [ 53 ]
انظر أيضاً
مراجع
- 1 2 3 ماونت، د.م. (2004). المعلوماتية الحيوية: تحليل التسلسل والجينوم ( الطبعة الثانية). مطبعة مختبر كولد سبرينغ هاربور: كولد سبرينغ هاربور، نيويورك. ISBN 978-0-87969-608-5.
- ↑ "أسئلة وأجوبة حول رموز Clustal" . Clustal . مؤرشف من الأصل بتاريخ 24 أكتوبر 2016. تم الاطلاع عليه بتاريخ 8 ديسمبر 2014 .
- ↑ نغ بي سي؛ هينيكوف إس (مايو 2001). "التنبؤ باستبدالات الأحماض الأمينية الضارة" . أبحاث الجينوم . 11 (5): 863-74 . doi : 10.1101/gr.176601 . PMC 311071. PMID 11337480 .
- 1 2 بوليانوفسكي، ف.و.؛ رويتبرغ، م.أ.؛ تومانيان، ف.ج. (2011). "تحليل مقارن لجودة خوارزمية شاملة وخوارزمية محلية لمحاذاة سلسلتين" . خوارزميات البيولوجيا الجزيئية . 6 (1): 25. doi : 10.1186/ 1748-7188-6-25 . PMC 3223492. PMID 22032267. S2CID 2658261 .
- ↑ شنايدر تي دي؛ ستيفنز آر إم (1990). "شعارات التسلسل: طريقة جديدة لعرض التسلسلات المتفق عليها" . مجلة أبحاث الأحماض النووية . 18 (20): 6097-6100 . doi : 10.1093 / nar/18.20.6097 . PMC 332411. PMID 2172928 .
- ↑ READSEQ
- ↑ "مواصفات تنسيق محاذاة التسلسل/الخريطة" (PDF) .
- ↑ برودنو م؛ مالدي س؛ بولياكوف أ؛ دو س ب؛ كورون أ؛ دوبتشاك إ؛ باتزوغلو س (2003). "المحاذاة العالمية المحلية: اكتشاف عمليات إعادة الترتيب أثناء المحاذاة". المعلوماتية الحيوية . 19. ملحق 1 (90001): i54–62. doi : 10.1093/bioinformatics/btg1005 . PMID 12855437 .
- ↑ ديلشر، أ. ل.؛ كاسيف، س.؛ فليشمان، ر. د.؛ بيترسون، ج.؛ وايت، أ.؛ سالزبرغ، س. ل. (1999). "محاذاة الجينومات الكاملة" . مجلة أبحاث الأحماض النووية . 27 (11): 2369-2376 . doi : 10.1093 / nar/30.11.2478 . PMC 148804. PMID 10325427 .
- ↑ وينغ-كين، سونغ (2010). الخوارزميات في المعلوماتية الحيوية: مقدمة عملية ( الطبعة الأولى). بوكا راتون: تشابمان آند هول/سي آر سي برس. رقم ISBN 978-1-4200-7033-0.
- ↑ غوتو، أوسامو (15 ديسمبر 1982). "خوارزمية محسّنة لمطابقة التسلسلات البيولوجية" . مجلة البيولوجيا الجزيئية . 162 (3): 705-708 . Bibcode : 1982JMBio.162..705G . doi : 10.1016/0022-2836(82)90398-9 . ISSN 0022-2836 . PMID 7166760 .
- ↑ غوتو، أوسامو (1 يناير 1999). "محاذاة التسلسلات المتعددة: الخوارزميات والتطبيقات" . التقدم في الفيزياء الحيوية . 36 : 159-206 . doi : 10.1016/S0065-227X(99)80007-0 . ISSN 0065-227X . PMID 10463075 .
- ↑ وانغ ل؛ جيانغ ت. (1994). "حول تعقيد محاذاة التسلسلات المتعددة". مجلة علم الأحياء الحاسوبي . 1 (4): 337-348 . Bibcode : 1994JCoB....1..337W . CiteSeerX 10.1.1.408.894 . doi : 10.1089/cmb.1994.1.337 . PMID 8790475 .
- ↑ إلياس، إسحاق (2006). "حل معضلة المحاذاة المتعددة". مجلة علم الأحياء الحاسوبي . 13 (7): 1323-1339 . CiteSeerX 10.1.1.6.256 . doi : 10.1089/cmb.2006.13.1323 . PMID 17037961 .
- ↑ ليبمان دي جيه؛ ألتشول إس إف؛ كيسيسي أوغلو جيه دي (1989). "أداة لمحاذاة التسلسلات المتعددة" . وقائع الأكاديمية الوطنية للعلوم في الولايات المتحدة الأمريكية . 86 ( 12): 4412-4415 . رمز Bibcode : 1989PNAS...86.4412L . doi : 10.1073/pnas.86.12.4412 . PMC 287279. PMID 2734293 .
- ↑ هيغينز دي جي ، شارب بي إم (1988). "CLUSTAL: حزمة برمجية لإجراء محاذاة تسلسل متعددة على حاسوب صغير". جين . 73 (1): 237-244 . doi : 10.1016/0378-1119(88)90330-7 . PMID 3243435 .
- ↑ تومسون، جيه دي؛ هيغينز، دي جي ؛ جيبسون، تي جيه. (1994). "CLUSTAL W: تحسين حساسية محاذاة التسلسل المتعدد التدريجي من خلال ترجيح التسلسل، وعقوبات الفجوات الخاصة بالموقع، واختيار مصفوفة الترجيح" . مجلة أبحاث الأحماض النووية . 22 (22): 4673-4680 . doi : 10.1093/nar/22.22.4673 . PMC 308517. PMID 7984417 .
- ↑ تشينا ر؛ سوغاوارا هـ؛ كويكي ت؛ لوبيز ر؛ جيبسون ت ج؛ هيغينز د ج؛ طومسون ج د. (2003). " محاذاة التسلسلات المتعددة باستخدام سلسلة برامج Clustal" . مجلة أبحاث الأحماض النووية . 31 (13): 3497-3500 . doi : 10.1093/nar/gkg500 . PMC 168907. PMID 12824352 .
- ↑ نوتردام سي؛ هيغينز دي جي ؛ هيرينغا جيه. (2000). "تي-كوفي: طريقة جديدة لمحاذاة تسلسل متعددة سريعة ودقيقة". مجلة علم الأحياء الجزيئي . 302 (1): 205-217 . doi : 10.1006/jmbi.2000.4042 . PMID 10964570. S2CID 10189971 .
- ↑ هيروساوا م؛ توتوكي ي؛ هوشيدا م؛ إيشيكاوا م. (1995). "دراسة شاملة حول الخوارزميات التكرارية لمحاذاة التسلسلات المتعددة". تطبيقات الحاسوب في العلوم البيولوجية . 11 (1): 13-18 . doi : 10.1093/bioinformatics/11.1.13 . PMID 7796270 .
- ↑ كاربلوس ك؛ باريت س؛ هيوجي ر. (1998). "نماذج ماركوف المخفية للكشف عن التشابهات البروتينية البعيدة" . المعلوماتية الحيوية . 14 (10): 846-856 . CiteSeerX 10.1.1.57.2762 . doi : 10.1093/bioinformatics/14.10.846 . PMID 9927713 .
- ↑ تشوثيا سي؛ ليسك إيه إم . (أبريل 1986). "العلاقة بين تباين التسلسل والبنية في البروتينات" . مجلة EMBO . 5 (4): 823-826 . doi : 10.1002/j.1460-2075.1986.tb04288.x . PMC 1166865. PMID 3709526 .
- 1 2 Zhang Y; Skolnick J. (2005). "يمكن حل مشكلة التنبؤ ببنية البروتين باستخدام مكتبة PDB الحالية" . وقائع الأكاديمية الوطنية للعلوم في الولايات المتحدة الأمريكية . 102 ( 4): 1029-34 . Bibcode : 2005PNAS..102.1029Z . doi : 10.1073/pnas.0407152101 . PMC 545829. PMID 15653774 .
- ↑ هولم إل؛ ساندر سي (1996). " رسم خريطة عالم البروتين". مجلة ساينس . 273 (5275): 595-603 . Bibcode : 1996Sci...273..595H . doi : 10.1126/science.273.5275.595 . PMID 8662544. S2CID 7509134 .
- ↑ تايلور دبليو آر؛ فلوريس تي بي؛ أورينغو سي إيه. (1994). "محاذاة بنية البروتين المتعددة" . علوم البروتين . 3 (10): 1858-1870 . doi : 10.1002/pro.5560031025 . PMC 2142613. PMID 7849601 .
- ↑ أورينغو، سي. أ.؛ ميتشي، أ. د.؛ جونز، س.؛ جونز، د. ت.؛ سوينديلز، م. ب.؛ ثورنتون، ج. م. (1997). "CATH - تصنيف هرمي لبنى نطاقات البروتين" . مجلة Structure . 5 (8): 1093-108 . doi : 10.1016/S0969-2126(97)00260-8 . PMID 9309224 .
- ↑ شينديالوف، آي. إن.؛ بورن، بي. إي. (1998). "محاذاة بنية البروتين عن طريق التمديد التوافقي التدريجي للمسار الأمثل". هندسة البروتين . 11 (9): 739-47 . doi : 10.1093/protein/11.9.739 . PMID 9796821 .
- 1 2 أورتيه ب؛ باستيان أ (2010). "من أين يأتي شكل توزيع درجة المحاذاة؟" . المعلوماتية الحيوية التطورية . 6 EBO.S5875: 159-187 . doi : 10.4137/EBO.S5875 . PMC 3023300. PMID 21258650 .
- ^ فيلسنشتاين ج. (2004). استنتاج السلالات . سيناور أسوشيتس: سندرلاند، MA. رقم ISBN 978-0-87893-177-4.
- ↑ ألتشول إس إف؛ غيش دبليو (1996). "إحصائيات المحاذاة المحلية". أساليب الحاسوب لتحليل تسلسل الجزيئات الكبيرة . أساليب في علم الإنزيمات. المجلد 266. الصفحات 460-480 . doi : 10.1016/S0076-6879(96)66029-7 . ISBN 978-0-12-182167-8PMID 8743700
- ↑ هارتمان، أ. ك. (2002). "أخذ عينات من الأحداث النادرة: إحصائيات محاذاة التسلسلات المحلية". مجلة الفيزياء E. 65 ( 5) 056102. arXiv : cond-mat/0108201 . Bibcode : 2002PhRvE..65e6102H . doi : 10.1103 /PhysRevE.65.056102 . PMID 12059642. S2CID 193085 .
- ↑ نيوبيرغ، ل . أ. (2008). "أهمية محاذاة التسلسلات ذات الفجوات" . مجلة علم الأحياء الحاسوبي . 15 (9): 1187-1194 . doi : 10.1089/cmb.2008.0125 . PMC 2737730. PMID 18973434 .
- ↑ إيدي إس آر؛ روست، بوركهارد (2008). روست، بوركهارد (محرر). "نموذج احتمالي لمحاذاة التسلسل المحلي يُبسط تقدير الدلالة الإحصائية" . مجلة PLOS للحوسبة البيولوجية . 4 (5) e1000069. رمز Bibcode : 2008PLSCB...4E0069E . doi : 10.1371/journal.pcbi.1000069 . PMC 2396288. PMID 18516236. S2CID 15640896 .
- ↑ باستيان أو؛ أود جي سي؛ روي إس؛ مارشال إي (2004). "أساسيات المحاذاة الثنائية الآلية الضخمة لتسلسلات البروتين: الأهمية النظرية لإحصائيات قيمة Z" . المعلوماتية الحيوية . 20 (4): 534-537 . CiteSeerX 10.1.1.602.6979 . doi : 10.1093/bioinformatics/btg440 . PMID 14990449 .
- ↑ أغراوال أ؛ هوانغ إكس (2011). "الأهمية الإحصائية الثنائية لمحاذاة التسلسل المحلي باستخدام مصفوفات الاستبدال الخاصة بالتسلسل والموقع". معاملات IEEE/ACM في علم الأحياء الحاسوبي والمعلوماتية الحيوية . 8 (1): 194-205 . Bibcode : 2011ITCBB...8..194A . doi : 10.1109/TCBB.2009.69 . PMID 21071807. S2CID 6559731 .
- ↑ أغراوال أ؛ بريندل ف ب؛ هوانغ إكس (2008). "الأهمية الإحصائية الثنائية والتحديد التجريبي لعقوبات فتح الفجوات الفعالة لمحاذاة التسلسل الموضعي للبروتين" . المجلة الدولية لعلم الأحياء الحاسوبي وتصميم الأدوية . 1 (4): 347-367 . doi : 10.1504/IJCBDD.2008.022207 . PMID 20063463 .
{{cite journal}}: CS1 maint: deprecated archiveal service ( link ) - ↑ نيوبيرغ إل إيه؛ لورانس سي إي (2009). " الحساب الدقيق للتوزيعات على الأعداد الصحيحة، مع تطبيق على محاذاة التسلسل" . مجلة علم الأحياء الحاسوبي . 16 (1): 1-18 . doi : 10.1089/cmb.2008.0137 . PMC 2858568. PMID 19119992 .
- ↑ كيم ن؛ لي سي (2008). "الكشف عن التضفير البديل باستخدام المعلوماتية الحيوية". المعلوماتية الحيوية . مناهج في البيولوجيا الجزيئية. المجلد 452. الصفحات 179-197 . doi : 10.1007/978-1-60327-159-2_9 . ISBN 978-1-58829-707-5PMID 18566765
- ↑ لي جيه بي، ليفانون إي واي، يون جيه كيه، وآخرون (مايو 2009). "تحديد مواقع تعديل الحمض النووي الريبوزي البشري على مستوى الجينوم من خلال التقاط الحمض النووي وتسلسله بالتوازي". مجلة ساينس . 324 (5931): 1210-1213 . رمز Bibcode : 2009Sci...324.1210L . doi : 10.1126/science.1170995 . PMID 19478186. S2CID 31148824 .
- ↑ بلازويتش ج، بريجا م، فيجليروفيتش م، وآخرون (يونيو 2009). "تجميع الجينوم الكامل من مخرجات تسلسل 454 عبر مفهوم رسم بياني معدل للحمض النووي". مجلة الكيمياء الحيوية الحاسوبية . 33 (3): 224-230 . doi : 10.1016/j.compbiolchem.2009.04.005 . PMID 19477687 .
- ↑ دوران سي؛ أبلبي إن؛ فاردي إم؛ إيميلفورت إم؛ إدواردز دي؛ باتلي جيه (مايو 2009). "اكتشاف تعدد أشكال النوكليوتيدات المفردة في الشعير باستخدام autoSNPdb" . مجلة التكنولوجيا الحيوية النباتية. 7 ( 4): 326-333 . Bibcode : 2009PBioJ...7..326D . doi : 10.1111/j.1467-7652.2009.00407.x . PMID 19386041 .
- ↑ أبوت، أ.؛ تساي، أ. (2000). "تحليل التسلسل وطرق المطابقة المثلى في علم الاجتماع: مراجعة وتوقعات". مناهج وبحوث علم الاجتماع . 29 (1): 3-33 . doi : 10.1177/0049124100029001001 . S2CID 121097811 .
- ↑ بارزيلاي ر؛ لي ل. (2002). "تحديد اختيار المفردات باستخدام محاذاة التسلسلات المتعددة" (ملف PDF) . وقائع مؤتمر ACL-02 حول الأساليب التجريبية في معالجة اللغة الطبيعية - EMNLP '02 . المجلد 10. الصفحات 164-171 . arXiv : cs/0205065 . Bibcode : 2002cs........5065B . doi : 10.3115/1118693.1118715 . S2CID 7521453 .
- ↑ كوندراك، غريغورز (2002). خوارزميات لإعادة بناء اللغة (ملف PDF) (أطروحة). جامعة تورنتو. مؤرشف من الأصل (ملف PDF) في 17 ديسمبر 2008. تم الاطلاع عليه في 21 يناير 2007 .
- ↑ برينزي أ.؛ د. فان دن بويل (2006). "دمج المعلومات التسلسلية في نماذج التصنيف التقليدية باستخدام نموذج SAM حساس للعنصر/الموقع" . نظم دعم القرار . 42 (2): 508-526 . doi : 10.1016/j.dss.2005.02.004 .انظر أيضًا ورقة برينزي وفان دن بويل البحثية : برينزي، أ؛ فان دن بويل، د (2007). "التنبؤ بتسلسل اقتناء الأجهزة المنزلية: ماركوف/ماركوف للتمييز وتحليل البقاء لنمذجة المعلومات التسلسلية في نماذج NPTB" . نظم دعم القرار . 44 (1): 28-45 . doi : 10.1016/j.dss.2007.02.008 .
- ↑ EMBL-EBI. "ClustalW2 < محاذاة التسلسل المتعدد < EMBL-EBI" . www.EBI.ac.uk. تم الاطلاع عليه بتاريخ 12 يونيو 2017 .
- ↑ تي-كوفي
- ↑ "BLAST: أداة البحث الأساسية للمحاذاة المحلية" . blast.ncbi.nlm.NIH.gov . تم الاطلاع عليه بتاريخ 12 يونيو 2017 .
- ↑ "خادم UVA FASTA" . fasta.bioch.Virginia.edu . تم الاطلاع عليه بتاريخ 12 يونيو 2017 .
- ↑ تومسون، جيه دي؛ بليونياك، إف؛ بوتش، أو (1999). "BAliBASE: قاعدة بيانات مرجعية للمحاذاة لتقييم برامج المحاذاة المتعددة" . المعلوماتية الحيوية . 15 (1): 87-88 . doi : 10.1093/bioinformatics/15.1.87 . PMID 10068696 .
- ↑ BAliBASE
- ↑ تومسون، جيه دي؛ بلونياك ، إف؛ بوتش، أو. (1999). "مقارنة شاملة لبرامج محاذاة التسلسلات المتعددة" . مجلة أبحاث الأحماض النووية . 27 (13): 2682-2690 . doi : 10.1093/nar/27.13.2682 . PMC 148477. PMID 10373585 .
- ↑ "محاذاة التسلسل المتعدد: Strap" . 3d-alignment.eu . تم الاطلاع عليه بتاريخ 12 يونيو 2017 .
روابط خارجية
الوسائط المتعلقة بمحاذاة التسلسل على ويكيميديا كومنز
- خوارزميات المعلوماتية الحيوية
- علم الوراثة الحاسوبي
- خوارزميات محاذاة التسلسل
- علم الأحياء التطوري النمائي
- خوارزميات على السلاسل النصية


