نموذج NK

نموذج NK هو نموذج رياضي وصفه مخترعه الرئيسي ستيوارت كوفمان بأنه سطح لياقة "قابل للتعديل" . وتجسد عبارة "قابل للتعديل" فكرة إمكانية تعديل كل من الحجم الإجمالي لسطح اللياقة وعدد "التلال والوديان" المحلية فيه من خلال تغيير معياريه.شمال{\displaystyle N}وك{\displaystyle K}، معشمال{\displaystyle N}كونه طول سلسلة من التطور وك{\displaystyle K}تحديد مستوى وعورة المناظر الطبيعية.

لقد وُظِّف نموذج NK في مجالاتٍ عديدة، منها الدراسات النظرية لعلم الأحياء التطوري ، وعلم المناعة ، والتحسين ، والتطور التكنولوجي ، وعلم الفريق، [ 1 ] والأنظمة المعقدة . كما اعتُمد النموذج في نظرية التنظيم ، حيث يُستخدم لوصف كيفية بحث كيانٍ ما في بيئةٍ ما من خلال التلاعب بخصائصه المختلفة. فعلى سبيل المثال، يمكن أن يكون الكيان منظمة ، وتمثل التلال والوديان الربح (أو تغيراته)، ويتطلب التحرك في البيئة قراراتٍ تنظيمية (مثل إضافة خطوط إنتاج أو تغيير الهيكل التنظيمي)، والتي تتفاعل فيما بينها وتؤثر على الربح بطريقةٍ معقدة. [ 2 ]

نسخة مبكرة من النموذج، والتي لم تأخذ في الاعتبار سوى الأكثر سلاسة (ك=0{\displaystyle K=0}) والأكثر صلابة (ك=شمال-1{\displaystyle K=N-1}تم تقديم مفهوم المناظر الطبيعية في كتاب كوفمان وليفين (1987). [ 3 ] أما النموذج كما هو معروف حاليًا فقد ظهر لأول مرة في كتاب كوفمان وواينبرغر (1989). [ 4 ]

أحد أسباب اكتساب هذا النموذج اهتمامًا واسعًا في مجال التحسين هو كونه حالةً بسيطةً للغاية لما يُعرف بمسألة NP-كاملة [ 5 ] ، ما يعني صعوبة إيجاد الحلول المثلى الشاملة. وقد أُثبت مؤخرًا أن نموذج NK لـ K > 1 هو أيضًا مسألة PLS-كاملة [ 6 ] ، ما يعني صعوبة إيجاد حتى الحلول المثلى المحلية للياقة البدنية. ولهذا الأمر تبعاتٌ على دراسة التطور المفتوح .

مثال نموذجي: لياقة البلازميد

البلازميد عبارة عن حلقة صغيرة من الحمض النووي (DNA) داخل خلايا معينة ، ويمكنه التكاثر بشكل مستقل عن الخلايا المضيفة. لنفترض أننا نرغب في دراسة مدى كفاءة البلازميدات.

لتبسيط الأمر، نمثل البلازميد كحلقة من N جينًا محتملاً، مرتبة دائمًا بنفس الترتيب، ولكل منها حالتان محتملتان (نشط أو غير نشط، من النوع X أو من النوع Y، إلخ...). ثم يُنمذج البلازميد بسلسلة ثنائية طولها N ، وبالتالي تكون دالة اللياقة هيF:{0،1}شمالR{\displaystyle F:\{0,1\}^{N}\to \mathbb {R} }.

أبسط نموذج هو الذي لا تتفاعل فيه الجينات مع بعضها البعض، وبالتالي نحصل علىF(S1S2Sشمال)=و1(S1)+و2(S2)++وشمال(Sشمال){\displaystyle F(S_{1}S_{2}\cdots S_{N})=f_{1}(S_{1})+f_{2}(S_{2})+\cdots +f_{N}(S_{N})}حيث كلوأنا(Sأنا){\displaystyle f_{i}(S_{i})}يشير إلى مساهمة الجين في اللياقةSأنا{\displaystyle S_{i}}في الموقعأنا{\displaystyle i}.

لنمذجة التفاعل الجيني ، نُدخل عاملاً آخر K ، وهو عدد الجينات الأخرى التي يتفاعل معها جينٌ ما. من المعقول افتراض أنه على البلازميد، يتفاعل جينان إذا كانا متجاورين، مما يُعطيF(S1S2Sشمال)=و1(S1،S2)+و2(S2،S3)++وشمال-1(Sشمال-1،Sشمال)+وشمال(Sشمال،S1){\displaystyle F(S_{1}S_{2}\cdots S_{N})=f_{1}(S_{1},S_{2})+f_{2}(S_{2},S_{3})+\cdots +f_{N-1}(S_{N-1},S_{N})+f_{N}(S_{N},S_{1})}على سبيل المثال، عندما تكون K = 1 و N = 5 ،

F(00101)=و1(0،0)+و2(0،1)+و3(1،0)+و4(0،1)+و5(1،0){\displaystyle F(00101)=f_{1}(0,0)+f_{2}(0,1)+f_{3}(1,0)+f_{4}(0,1)+f_{5}(1,0)}

يقوم نموذج NK بتعميم هذا من خلال السماح بقيم K و N المحدودة التعسفية، بالإضافة إلى السماح بتعريف تعسفي لتجاور الجينات (لا تقع الجينات بالضرورة على دائرة أو قطعة مستقيمة).

التعريف الرياضي

يُعرّف نموذج NK فضاءً طوريًا توافقيًا ، يتكون من كل سلسلة (مختارة من أبجدية معينة) بطولشمال{\displaystyle N}لكل سلسلة نصية في فضاء البحث هذا، يتم تحديد قيمة عددية (تُسمى الملاءمة ). إذا تم تحديد مقياس مسافة بين السلاسل النصية، فإن البنية الناتجة تُسمى مشهدًا .

تُحدد قيم اللياقة وفقًا للنسخة المحددة من النموذج، ولكن السمة الرئيسية لنموذج NK هي أن لياقة سلسلة معينةS{\displaystyle S}هو مجموع المساهمات من كل موضعوأنا(S){\displaystyle f_{i}(S)}في السلسلة:

F(S)=أناو~أنا(S)،{\displaystyle F(S)=\sum _{i}{\tilde {f}}_{i}(S),}

وتعتمد مساهمة كل موقع بشكل عام على حالته وحالةك{\displaystyle K}مواقع أخرى:

و~أنا(S)=وأنا(Sأنا،Sكأنا1،...،Sكأناك)،{\displaystyle {\tilde {f}}_{i}(S)=f_{i}(S_{i},S_{k_{i1}},\dots ,S_{k_{iK}}),}

أينكأناج{\displaystyle k_{ij}}هو فهرسج{\displaystyle j}الجار رقم 1 للموقعأنا{\displaystyle i}.

وبالتالي، فإن دالة اللياقةوأنا{\displaystyle f_{i}}هي عبارة عن دالة تربط بين سلاسل طولها K  +  1 وقيم عددية، والتي يطلق عليها واينبرغر في أعماله اللاحقة اسم "مساهمات اللياقة". وغالبًا ما يتم اختيار مساهمات اللياقة هذه عشوائيًا من توزيع احتمالي محدد .

تصوير ثنائي الأبعاد لمشهد لياقة NK. تمثل الأسهم مسارات الطفرات المختلفة التي يمكن أن يتبعها السكان أثناء تطورهم على مشهد اللياقة.

مثال: نماذج الزجاج الدوار

يُكتب نموذج إيزينغ أحادي البعد للزجاج المغزلي عادةً على النحو التاليح=-أنا=1شمالجأنا،أنا+1SأناSأنا+1-μأنا=1شمالحأناSأنا{\displaystyle H=-\sum _{i=1}^{N}J_{i,i+1}S_{i}S_{i+1}-\mu \sum _{i=1}^{N}h_{i}S_{i}}أينح{\displaystyle H}هو الهاميلتوني، والذي يمكن اعتباره طاقة.

يمكننا إعادة صياغتها كحالة خاصة من نموذج NK مع K=1 :ح=F(S)=أنا،جوأنا،ج(Sأنا،Sج){\displaystyle H=F(S)=\sum _{i,j}f_{i,j}(S_{i},S_{j})}من خلال تحديدوأنا(Sأنا،Sأنا+1)=-جأنا،أنا+1SأناSأنا+1-μحأناSأنا{\displaystyle f_{i}(S_{i},S_{i+1})=-J_{i,i+1}S_{i}S_{i+1}-\mu h_{i}S_{i}}بشكل عام، نموذج إيزينغ ذو الأبعاد m على شبكة مربعة{1،2،...،ن}م{\displaystyle \{1,2,...,n\}^{m}}هو نموذج NK معشمال=نم،ك=م{\displaystyle N=n^{m},K=m}.

بما أن K يقيس تقريبًا "وعورة" مشهد اللياقة البدنية (انظر أدناه)، فإننا نرى أنه مع زيادة بُعد نموذج Ising، تزداد وعورته أيضًا.

متىμ=0{\displaystyle \mu =0}هذا هو نموذج إدواردز-أندرسون، وهو قابل للحل تمامًا.

يعمم نموذج شيرينغتون-كيركباتريك نموذج إيزينغ من خلال السماح لجميع أزواج اللفات الممكنة بالتفاعل (بدلاً من استخدام رسم بياني شبكي ، استخدم الرسم البياني الكامل )، وبالتالي فهو أيضًا نموذج NK معك=شمال-1{\displaystyle K=N-1}.

بجعل جميع التسلسلات الفرعية الممكنة للدوران تتفاعل، بدلاً من الأزواج فقط، نحصل على نموذج المدى اللانهائي، وهو أيضاً نموذج NK معك=شمال-1{\displaystyle K=N-1}.

طوبولوجيا قابلة للتعديل

توضيح لطوبولوجيا قابلة للتعديل في نموذج NK. تمثل العقد سلاسل ثنائية منفصلة، ​​وتربط الحواف السلاسل بمسافة هامينغ تساوي واحدًا بالضبط. (يسار) N = 5، K = 0. (وسط) N = 5، K = 1. (يمين) N = 5، K = 2. يشير لون العقدة إلى لياقتها، حيث تدل القيم الأكثر احمرارًا على لياقة أعلى. تم اختيار تضمين المكعب الفائق بحيث تكون أعلى قيمة للياقة في المركز. لاحظ أن سطح البيانات عند K = 0 يبدو أكثر سلاسة من الحالات ذات قيم K الأعلى.

تتحكم قيمة K في درجة التفاعل الجيني (Epistasis) في نموذج NK، أو مدى تأثير المواقع الجينية الأخرى على مساهمة موقع جيني معين في اللياقة. عندما تكون K = 0، تكون لياقة سلسلة معينة عبارة عن مجموع بسيط لمساهمات المواقع الجينية الفردية: بالنسبة لدوال اللياقة غير التافهة، يوجد حل أمثل عالمي يسهل تحديده (مجموعة الأصفار إذا كانت f (0) > f (1)، أو مجموعة الآحاد إذا كانت f (1) > f (0)). أما عندما تكون K غير صفرية ، فإن لياقة السلسلة هي مجموع لياقات السلاسل الفرعية، والتي قد تتفاعل معًا لتُعيق النظام (انظر كيفية تحقيق اللياقة المثلى في المثال أعلاه). وبالتالي، فإن زيادة K تزيد من وعورة سطح اللياقة.

تنويعات مع مساحات محايدة

لا يدعم نموذج NK الأساسي ظاهرة الفضاء المحايد ، أي مجموعات الجينومات المرتبطة بطفرات مفردة لها نفس قيمة اللياقة. وقد اقتُرح تعديلان لإدراج هذا التركيب المهم بيولوجيًا . يُدخل نموذج NKP مُعاملًا.P{\displaystyle P}نسبةP{\displaystyle P}التابع2ك{\displaystyle 2^{K}}تُضبط مساهمات اللياقة على الصفر، بحيث تكون مساهمات العديد من الأنماط الجينية متداخلة . يُدخل نموذج NKQ مُعاملًاسؤال{\displaystyle Q}ويفرض تجزئة على قيم مساهمة اللياقة البدنية المحتملة بحيث تأخذ كل مساهمة إحدى القيم التالية:سؤال{\displaystyle Q}القيم المحتملة، مما يُدخل مرة أخرى انحلالًا في مساهمات بعض الأنماط الجينية . يتوافق نموذج NK الأساسي معP=0{\displaystyle P=0}وسؤال={\displaystyle Q=\infty }الحالات في ظل هذه المعايير.

النتائج المعروفة

في عام 1991، نشر واينبرغر تحليلاً مفصلاً [ 7 ] للقضية التي1<<كشمال{\displaystyle 1<<k\leq N}ويتم اختيار مساهمات اللياقة عشوائيًا. وقد تبين لاحقًا أن تقديره التحليلي لعدد الحلول المثلى المحلية معيب . ومع ذلك، تدعم التجارب العددية المضمنة في تحليل واينبرغر نتيجته التحليلية القائلة بأن اللياقة المتوقعة لسلسلة ما تتوزع توزيعًا طبيعيًا بمتوسط ​​يقارب

μ+σ2ln(ك+1)ك+1{\displaystyle \mu +\sigma {\sqrt {{2\ln(k+1)} \over {k+1}}}}

وتباين يقارب

(ك+1)σ2شمال[ك+1+2(ك+2)ln(ك+1)]{\displaystyle {{(k+1)\sigma ^{2}} \over {N[k+1+2(k+2)\ln(k+1)]}}}.

التطبيقات

لقد وجد نموذج NK استخدامًا في العديد من المجالات، بما في ذلك دراسة الزجاج المغزلي ، وحل المشكلات الجماعية ، [ 8 ] والتفاعل الجيني والتعددية في علم الأحياء التطوري ، والتحسين التوافقي .

مراجع

  1. بوروماند، أمين؛ سمولدينو، بول إي. (2023). "يمكن لتحيز التفوق والتشويش في التواصل أن يعززا حل المشكلات الجماعي" . مجلة المجتمعات الاصطناعية والمحاكاة الاجتماعية . 26 (3). doi : 10.18564/jasss.5154 .
  2. ليفينثال، د.أ. (1997). "التكيف في المناظر الطبيعية الوعرة". علوم الإدارة . 43 (7): 934-950 . doi : 10.1287/mnsc.43.7.934 .
  3. كوفمان، س.؛ ليفين، س. (1987). "نحو نظرية عامة للمشي التكيفي على المناظر الطبيعية الوعرة". مجلة البيولوجيا النظرية . 128 (1): 11-45 . Bibcode : 1987JThBi.128...11K . doi : 10.1016/s0022-5193(87)80029-2 . PMID 3431131 . 
  4. كوفمان، س.؛ واينبرغر، إ. (1989). "نموذج الخلايا القاتلة الطبيعية للمناظر الطبيعية الوعرة للياقة البدنية وتطبيقه على نضج الاستجابة المناعية". مجلة البيولوجيا النظرية . 141 (2): 211-245 . Bibcode : 1989JThBi.141..211K . doi : 10.1016/s0022-5193(89)80019-0 . PMID 2632988 . 
  5. Weinberger, E. (1996), “NP-completeness of Kauffman's Nk model, a Tuneably Rugged Fitness Landscape”, Santa Fe Institute Working Paper, 96-02-003.
  6. كازناتشيف، أرتيم (2019). "التعقيد الحسابي كقيد نهائي على التطور" . علم الوراثة . 212 (1): 245-265 . doi : 10.1534/genetics.119.302000 . PMC 6499524. PMID 30833289 .  
  7. واينبرغر، إدوارد (15 نوفمبر 1991). "الخصائص المحلية لنموذج كوفمان Nk: سطح طاقة قابل للتعديل". مجلة Physical Review A. 10. 44 (10): 6399–6413 . Bibcode : 1991PhRvA..44.6399W . doi : 10.1103/physreva.44.6399 . PMID 9905770 . 
  8. بوروماند، أ. وسمولدينو، ب.إ، 2021. العمل الجاد، والمخاطرة، والتنوع في نموذج لحل المشكلات الجماعية. مجلة المجتمعات الاصطناعية والمحاكاة الاجتماعية، 24(4).