التطبيقات البيولوجية لنظرية التشعب

توفر التطبيقات البيولوجية لنظرية التشعب إطارًا لفهم سلوك الشبكات البيولوجية التي تُنمذج كنظم ديناميكية . في سياق النظام البيولوجي، تصف نظرية التشعب كيف يمكن لتغيرات طفيفة في مُدخلات النظام أن تُحدث تشعبًا أو تغييرًا نوعيًا في سلوكه. غالبًا ما تكون القدرة على إحداث تغيير جذري في مخرجات النظام أساسية لوظيفة الكائن الحي، ولذلك فإن التشعبات شائعة في الشبكات البيولوجية، مثل مراحل دورة الخلية .

الشبكات البيولوجية والأنظمة الديناميكية

تنشأ الشبكات البيولوجية من التطور ، ولذلك فهي تتميز بمكونات أقل توحيدًا وتفاعلات أكثر تعقيدًا من الشبكات التي يصممها الإنسان، كالشبكات الكهربائية . على المستوى الخلوي، قد تشمل مكونات الشبكة مجموعة واسعة من البروتينات، يختلف الكثير منها بين الكائنات الحية. تحدث تفاعلات الشبكة عندما يؤثر بروتين واحد أو أكثر على وظيفة بروتين آخر من خلال النسخ ، أو الترجمة ، أو الانتقال ، أو الفسفرة ، أو آليات أخرى. تُفعّل هذه التفاعلات عمل البروتين المستهدف أو تثبطه بطريقة ما. بينما يبني الإنسان الشبكات مع مراعاة البساطة والنظام، تكتسب الشبكات البيولوجية التكرار والتعقيد عبر مسار التطور. لذلك، قد يكون من المستحيل التنبؤ بالسلوك الكمي للشبكة البيولوجية من خلال معرفة تنظيمها. وبالمثل، من المستحيل وصف تنظيمها من خلال سلوكها فقط، مع أن السلوك قد يشير إلى وجود أنماط شبكية معينة .

الشكل 1. مثال على شبكة بيولوجية بين الجينات والبروتينات التي تتحكم في الدخول في طور S.

مع ذلك، فإنه بمعرفة تفاعلات الشبكة ومجموعة من المعايير الخاصة بالبروتينات وتفاعلاتها (والتي تُستخلص عادةً من خلال البحوث التجريبية )، يصبح من الممكن في كثير من الأحيان بناء نموذج للشبكة كنظام ديناميكي . بشكل عام، بالنسبة لـن{\displaystyle n}بالنسبة للبروتينات، يأخذ النظام الديناميكي الشكل التالي [ 1 ] حيثx{\displaystyle x}عادةً ما يكون تركيز البروتين:

x˙1=دx1دت=و1(x1،...،xن){\displaystyle {\dot {x}}_{1}={\frac {dx_{1}}{dt}}=f_{1}(x_{1},\ldots ,x_{n})}
{\displaystyle \vdots }
x˙أنا=دxأنادت=وأنا(x1،...،xن){\displaystyle {\dot {x}}_{i}={\frac {dx_{i}}{dt}}=f_{i}(x_{1},\ldots ,x_{n})}
{\displaystyle \vdots }
x˙ن=دxندت=ون(x1،...،xن){\displaystyle {\dot {x}}_{n}={\frac {dx_{n}}{dt}}=f_{n}(x_{1},\ldots ,x_{n})}

غالباً ما تكون هذه الأنظمة صعبة الحل، لذا فإن نمذجة الشبكات كنظم ديناميكية خطية أسهل. لا تحتوي الأنظمة الخطية على أي نواتج بينx{\displaystyle x}يمكن حل المعادلات دائمًا. ولها الشكل التالي لجميعأنا{\displaystyle i}:

وأنا=أأنا1x1+أأنا2x2++أأنانxن{\displaystyle f_{i}=a_{i1}x_{1}+a_{i2}x_{2}+\cdots +a_{in}x_{n}\,}

لسوء الحظ، غالباً ما تكون الأنظمة البيولوجية غير خطية ، وبالتالي فهي تحتاج إلى نماذج غير خطية.

أنماط الإدخال/الإخراج

على الرغم من التعقيد والتنوع الكبيرين المحتملين للشبكات البيولوجية، فإن سلوك الشبكة من الدرجة الأولى يُعمم على أحد أربعة أنماط محتملة للمدخلات والمخرجات: النمط الزائدي أو نمط ميكايليس-مينتين ، والنمط فائق الحساسية ، والنمط ثنائي الاستقرار ، والنمط ثنائي الاستقرار غير العكوس (وهو نمط ثنائي الاستقرار يتطلب مدخلات سالبة، وبالتالي مستحيلة بيولوجيًا، للعودة من حالة المخرجات العالية). يمكن الاطلاع على أمثلة لكل نمط في سياقات بيولوجية على صفحاتها الخاصة.

تُظهر الشبكات فائقة الحساسية، والشبكات ثنائية الاستقرار، والشبكات ثنائية الاستقرار غير القابلة للعكس، جميعها تغييراً نوعياً في سلوك الشبكة حول قيم معينة للمعلمات  - وهذه هي نقاط التشعب الخاصة بها.

التشعبات الأساسية في وجود الخطأ

الشكل 2. مخطط طور تشعب العقدة السرجية، حيث يتم تغيير معامل التحكم (يُشار إليه بـ ε بدلاً من r ، ولكنه مكافئ وظيفيًا). عندما تنخفض قيمة ε، تتقارب النقاط الثابتة وتتلاشى؛ وعندما تزداد قيمة ε، تظهر النقاط الثابتة. يُرمز إلى معدل تغير dx/dt بالرمز v.

يمكن فهم الأنظمة الديناميكية غير الخطية بسهولة أكبر من خلال نظام مثال أحادي البعد حيث يعتمد التغير في كمية معينة x (مثل تركيز البروتين) على نفسها فقط:

x˙=دxدت=و(x){\displaystyle {\dot {x}}={\frac {dx}{dt}}=f(x)\,}

بدلاً من حل النظام تحليلياً، وهو ما قد يكون صعباً أو مستحيلاً بالنسبة للعديد من الدوال، غالباً ما يكون من الأسرع والأكثر إفادة اتباع نهج هندسي ورسم مخطط طوري . المخطط الطوري هو رسم تخطيطي نوعي لسلوك المعادلة التفاضلية، يُظهر حلول التوازن أو النقاط الثابتة وحقل المتجهات على خط الأعداد الحقيقية.

تصف التشعبات تغيرات في استقرار أو وجود نقاط ثابتة مع تغير أحد معايير التحكم في النظام. ولتبسيط شرح التشعب في نظام ديناميكي، تخيل جسمًا متوازنًا على عارضة رأسية. يمكن اعتبار كتلة الجسم معيار التحكم ( r )، وانحراف العارضة عن المحور الرأسي المتغير الديناميكي (x) . مع ازدياد قيمة r ، يبقى x مستقرًا نسبيًا. ولكن عندما تصل الكتلة إلى نقطة معينة  - نقطة التشعب  - تنحني العارضة فجأة، في اتجاه يعتمد على عيوب طفيفة في التركيب. هذا مثال على تشعب الشوكة. أدت التغيرات في معيار التحكم في النهاية إلى تغيير السلوك النوعي للنظام.

تشعب العقدة السرجية

للحصول على مثال أكثر دقة، انظر إلى النظام الديناميكي الموضح في الشكل 2، والموصوف بالمعادلة التالية:

x˙=-x2+ر{\displaystyle {\dot {x}}=-x^{2}+r}

حيث يُمثل r مرة أخرى مُعامل التحكم (المُشار إليه بـ ε في الشكل 2). تُمثل النقاط الثابتة للنظام بنقاط تقاطع منحنى مخطط الطور مع المحور السيني. يُمكن تحديد استقرار نقطة ثابتة مُعينة من خلال اتجاه التدفق على المحور السيني؛ على سبيل المثال، في الشكل 2، النقطة الخضراء غير مستقرة (تدفق مُتباعد)، والنقطة الحمراء مستقرة (تدفق مُتقارب). في البداية، عندما تكون قيمة r أكبر من صفر، يمتلك النظام نقطة ثابتة مستقرة ونقطة ثابتة غير مستقرة. مع انخفاض قيمة تتحرك النقاط الثابتة معًا، وتصطدم لفترة وجيزة لتُشكل نقطة ثابتة شبه مستقرة عند r = 0، ثم تختفي عندما تكون r < 0.

في هذه الحالة، ولأن سلوك النظام يتغير بشكل ملحوظ عندما تكون قيمة معامل التحكم r تساوي صفرًا، فإن الصفر يمثل نقطة تشعب . ومن خلال تتبع مواقع النقاط الثابتة في الشكل 2 مع تغير قيمة r ، يمكن إنشاء مخطط التشعب الموضح في الشكل 3.

تُعدّ أنواع أخرى من التشعبات مهمة أيضًا في الأنظمة الديناميكية، لكن تشعب العقدة السرجية يميل إلى أن يكون الأكثر أهمية في علم الأحياء. تخضع الأنظمة البيولوجية الحقيقية لتغيرات عشوائية طفيفة تُدخل حدود خطأ في المعادلات الديناميكية، وهذا عادةً ما يؤدي إلى تشعبات أكثر تعقيدًا تُبسط إلى عقد سرجية ونقاط ثابتة منفصلة. نناقش أدناه مثالين على التشعبات "غير الكاملة" التي قد تظهر في علم الأحياء. تجدر الإشارة إلى أن العقدة السرجية نفسها، في وجود الخطأ، تنتقل ببساطة في المستوى xr ، دون أي تغيير في السلوك النوعي؛ ويمكن إثبات ذلك باستخدام التحليل نفسه الموضح أدناه.

الشكل 4. تشعبات عابرة غير مضطربة (باللون الأسود) وغير كاملة (باللون الأحمر)، متراكبة. يُشار إلى u و p بالرمزين x و r على التوالي في بقية المقال. وكما في السابق، الخطوط المتصلة تمثل التشعبات المستقرة، والخطوط المتقطعة تمثل التشعبات غير المستقرة.

تشعب غير كامل فوق الحرج

يُعد التشعب العابر أحد أنواع التشعبات البسيطة الشائعة ، ويُعطى بواسطة

دxدت=رx-x2{\displaystyle {dx \over dt}=rx-x^{2}} ومخطط التشعب في الشكل 4 (المنحنيات السوداء). تظهر مخططات الأطوار في الشكل 5. عند تتبع نقاط التقاطع مع المحور السيني في مخطط الطور مع تغير قيمة r ، نجد مسارين ثابتين يتقاطعان عند نقطة الأصل؛ وهذه هي نقطة التشعب (بشكل بديهي، عندما يتغير عدد نقاط التقاطع مع المحور السيني في مخطط الطور). النقطة الثابتة اليسرى غير مستقرة دائمًا، بينما اليمنى مستقرة.

الشكل 5. مخططات طور التشعب فوق الحرج المثالي. تُحدد النقاط الثابتة على المحور السيني، كل مسار بلون مختلف. يشير اتجاه الأسهم إلى اتجاه حركتها مع ازدياد قيمة r . النقطة الحمراء غير مستقرة، والنقطة الزرقاء غير مستقرة. النقطة السوداء عند نقطة الأصل مستقرة عندما r < 0، وغير مستقرة عندما r > 0.

والآن، لنفترض إضافة حد خطأ h ، حيث 0 < h << 1. أي،

دxدت=رx-x2-ح{\displaystyle {dx \over dt}=rx-x^{2}-h}

يُزيح حد الخطأ جميع مخططات الطور عموديًا، نحو الأسفل إذا كانت قيمة h موجبة. في النصف الأيسر من الشكل 6 ( x < 0)، تُمثل النقاط الثابتة السوداء والحمراء والخضراء نقاطًا شبه مستقرة وغير مستقرة ومستقرة على التوالي. وينعكس هذا في النقاط الأرجوانية والسوداء والزرقاء في النصف الأيمن ( x > 0). وبالتالي، يتصرف كل نصف من هذه النصفين كنقطة تشعب عقدة سرجية؛ بمعنى آخر، يمكن تقريب التشعب العابر غير الكامل بتشعبين عقدة سرجية عند الاقتراب من النقاط الحرجة، كما هو واضح في المنحنيات الحمراء في الشكل 4.

تحليل الاستقرار الخطي

الشكل 6. مخططات طور التشعب فوق الحرج غير الكامل. تظهر خمس قيم لـ r ، مُعطاة بالنسبة للنقطتين الحرجتين. لاحظ أن قيمة التقاطع مع المحور y هي نفسها h ، أو مقدار عدم الكمال. النقاط الخضراء والزرقاء مستقرة، بينما النقاط الخضراء والحمراء والبنفسجية غير مستقرة. تشير النقاط السوداء إلى نقاط ثابتة شبه مستقرة.

إلى جانب مراقبة التدفق في مخططات الطور، من الممكن أيضًا إثبات استقرار النقاط الثابتة المختلفة باستخدام تحليل الاستقرار الخطي . أولًا، يتم إيجاد النقاط الثابتة في مخطط الطور عن طريق جعل معادلة التشعب تساوي صفرًا.

دxدت=و(x)=رx-(x)2-ح0=رx*-(x*)2-ح{\displaystyle {\begin{aligned}{dx \over dt}=f(x)&=rx-(x)^{2}-h\\0&=rx^{*}-(x^{*})^{2}-h\end{aligned}}}

استخدام الصيغة التربيعية لإيجاد النقاط الثابتة x* :

x*=-ر±ر2-4(-1)ح2(-1)=ر±ر2+4ح2ر±ر22{\displaystyle {\begin{aligned}x^{*}&={-r\pm {\sqrt {r^{2}-4(-1)h}} \over 2(-1)}\\&={r\pm {\sqrt {r^{2}+4h}} \over 2}\\&\approx {r\pm {\sqrt {r^{2}}} \over 2}\end{aligned}}}

حيث تم في الخطوة الأخيرة استخدام التقريب 4h << r2 ، وهو أمر معقول لدراسة النقاط الثابتة التي تتجاوز نقطة التشعب، مثل المنحنيات الزرقاء الفاتحة والخضراء في الشكل 6. وبتبسيط الأمر أكثر،

x*ر±ر2={0،نقطة ثابتة أقرب إلى الأصلر،نقطة ثابتة أبعد عن نقطة الأصل{\displaystyle {\begin{aligned}x^{*}&\approx {r\pm r \over 2}\\&={\begin{cases}0,&{\text{نقطة ثابتة أقرب إلى الأصل}}\\r,&{\text{نقطة ثابتة أبعد عن الأصل}}\end{cases}}\end{aligned}}}

بعد ذلك، حدد ما إذا كان منحنى صورة الطور يتزايد أو يتناقص عند النقاط الثابتة، وهو ما يمكن تقييمه عن طريق إدخال x * في المشتق الأول لمعادلة التشعب.

و(x)=ر-2xو(0)=ر={>0،لو ر>0غير مستقر (أرجواني)<0،لو ر<0مستقر (أخضر)و(ر)=-ر={<0،لو ر>0ثابت (أزرق)>0،لو ر<0غير مستقر (أحمر){\displaystyle {\begin{aligned}f'(x)&=r-2x\\f'(0)&=r={\begin{cases}>0,&{\text{إذا كان }}r>0\rightarrow {\text{غير مستقر (أرجواني)}}\\<0,&{\text{إذا كان }}r<0\rightarrow {\text{مستقر (أخضر)}}\end{cases}}\\f'(r)&=-r={\begin{cases}<0,&{\text{إذا كان }}r>0\rightarrow {\text{مستقر (أزرق)}}\\>0,&{\text{إذا كان }}r<0\rightarrow {\text{غير مستقر (أحمر)}}\end{cases}}\end{aligned}}}

تتعقد النتائج بسبب إمكانية أن تكون قيمة r موجبة أو سالبة؛ ومع ذلك، تبقى الاستنتاجات كما هي فيما يتعلق باستقرار كل نقطة ثابتة. وهذا ليس بالأمر المفاجئ، إذ أن المشتقة الأولى تحتوي على نفس المعلومات التي يوفرها تحليل تدفق مخطط الطور. وتتوافق الألوان في الحل أعلاه مع الأسهم في الشكل 6.

تشعب غير كامل على شكل شوكة

يُعد مثال العارضة المنحنية المذكور سابقًا مثالًا على تشعب الشوكة (أو ربما يُطلق عليه بشكل أدق "التشعب الثلاثي"). يظهر تشعب الشوكة "المثالي" على يسار الشكل 7، كما هو موضح في

دxدت=رx-x3{\displaystyle {dx \over dt}=rx-x^{3}}

عند r = 0 يحدث التشعب، ويمثله النقطة السوداء في نقطة الأصل في الشكل 8. مع ازدياد قيمة r بعد الصفر، تنقسم النقطة السوداء إلى ثلاثة مسارات: النقطة الزرقاء الثابتة المستقرة التي تتحرك يمينًا، والنقطة الحمراء الثابتة التي تتحرك يسارًا، ونقطة ثالثة غير مستقرة تبقى عند نقطة الأصل. الخطان الأزرق والأحمر متصلان في الشكل 7 (يسارًا)، بينما يمثل الخط المنقط على طول المحور السيني الموجب المسار الأسود غير المستقر.

كما في السابق، ضع في اعتبارك حد الخطأ h ، حيث 0 < h << 1، أي

الشكل 8. مخططات طور التشعب المثالي. تُحدد النقاط الثابتة على المحور السيني، كل مسار بلون مختلف. يشير اتجاه الأسهم إلى اتجاه حركتها مع ازدياد قيمة r . النقطتان الحمراء والزرقاء مستقرتان، وهناك نقطة ثابتة ثالثة غير مستقرة عند نقطة الأصل، مُشار إليها بالنقطة السوداء. عند r = r crit = 0، تُشير النقطة السوداء أيضًا إلى نقطة شبه مستقرة تظهر وتنقسم إلى المسارات الثلاثة الأخرى مع ازدياد قيمة r.

دxدت=رx-x3+ح{\displaystyle {dx \over dt}=rx-x^{3}+h}

مرة أخرى، تُنقل مخططات الطور لأعلى بمقدار ضئيل للغاية، كما هو موضح في الشكل 9. يُؤدي تتبع نقاط التقاطع مع المحور السيني في مخطط الطور مع تغير قيمة r إلى تحديد النقاط الثابتة، والتي تُعيد إنتاج النتيجة النوعية من الشكل 7 (على اليمين). وبشكل أكثر تحديدًا، تُقابل النقطة الثابتة الزرقاء في الشكل 9 المسار العلوي في الشكل 7 (على اليمين)؛ والنقطة الثابتة الخضراء هي المسار المنقط؛ والنقطة الثابتة الحمراء هي المسار السفلي. وبالتالي، في الحالة غير المثالية ( h ≠ 0)، يتبسط تشعب الشوكة إلى نقطة ثابتة مستقرة واحدة مقترنة بتشعب عقدة سرجية.

يمكن إجراء تحليل الاستقرار الخطي هنا أيضًا، ولكن باستخدام الحل المعمم لمعادلة تكعيبية بدلًا من معادلة تربيعية. وتكون العملية نفسها: 1) مساواة المعادلة التفاضلية بالصفر وإيجاد الصيغة التحليلية للنقاط الثابتة x* ، 2) تعويض كل x* في المشتقة الأولىو(x)=ددxدxدت{\displaystyle f'(x)={d \over dx}{dx \over dt}}ثم 3) تقييم الاستقرار بناءً على ما إذاو(x*){\displaystyle f'(x^{*})}إيجابي أم سلبي.

الشكل 9. مخططات طورية لتشعب الشوكة غير الكاملة. تظهر أربع قيم مختلفة لـ r نسبةً إلى r crit . لاحظ أن قيمة التقاطع مع المحور y هي نفسها h ، أو مقدار عدم الكمال. النقطتان الحمراء والزرقاء مستقرتان، بينما النقطة الخضراء (التي كانت مخفية سابقًا عند نقطة الأصل) غير مستقرة. كما في الشكل 5، تشير النقطة السوداء إلى نقطة شبه مستقرة تظهر وتنقسم إلى النقطتين الحمراء والخضراء مع ازدياد قيمة r.

تعدد الاستقرار

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

أمثلة محددة في علم الأحياء

تتحكم الشبكات ذات التشعبات الديناميكية في العديد من التحولات المهمة في دورة الخلية . تعمل تحولات G1/S و G2/M والطور الاستوائي-الطور الانفصالي كمفاتيح كيميائية حيوية في دورة الخلية . على سبيل المثال، يتم تحفيز مستخلصات بيض الضفدع الأفريقي (Xenopus laevis) للدخول والخروج من الانقسام المتساوي بشكل لا رجعة فيه عن طريق التغذية الراجعة الإيجابية في فسفرة Cdc2، وهو كيناز معتمد على السيكلين . [ 3 ]

في علم بيئة السكان ، يمكن أن تُظهر ديناميكيات شبكات التفاعلات الغذائية تشعبات هوبف . على سبيل المثال، في نظام مائي يتكون من منتج أولي ، ومورد معدني، وحيوان عاشب، وجد الباحثون أن أنماط التوازن، والدوران، وانقراض السكان يمكن وصفها نوعيًا بنموذج غير خطي بسيط مع تشعب هوبف. [ 4 ]

يمكن قياس استهلاك الجلاكتوز في الخميرة المتبرعمة (S. cerevisiae) من خلال التعبير عن بروتين GFP المُحفَّز بواسطة مُحفِّز GAL كدالة لتغير تركيزات الجلاكتوز. ويُظهر النظام تبديلاً ثنائي الاستقرار بين حالتي التحفيز وعدم التحفيز. [ 5 ]

وبالمثل، يُظهر استهلاك اللاكتوز في بكتيريا الإشريكية القولونية كدالة لتركيز ثيو-ميثيل جالاكتوزيد (نظير اللاكتوز) المقاس بواسطة محفز لاك المُعبِّر عن البروتين الفلوري الأخضر، استقرارًا ثنائيًا وتخلفًا (الشكل 10، على اليسار واليمين على التوالي). [ 6 ]

انظر أيضاً

مراجع

  1. ستروغاتز إس إتش (1994)، الديناميكا غير الخطية والفوضى، دار نشر بيرسيوس بوكس
  2. ديفيد أنجيلي، جيمس إي. فيريل الابن، وإدواردو د. سونتاج. الكشف عن تعدد الاستقرار، والتشعبات، والتخلف في فئة كبيرة من أنظمة التغذية الراجعة الإيجابية البيولوجية. وقائع الأكاديمية الوطنية للعلوم، 17 فبراير 2004، المجلد 101، العدد 7، الصفحات 1822-1827
  3. شا، وي؛ مور، جوناثان؛ تشين، كاثرين؛ لاساليتا، أنطونيو د.؛ يي، تشونغ سيون؛ تايسون، جون ج.؛ سيبل، جيل س. (4 فبراير 2003). "التخلف الخلوي يحفز تحولات دورة الخلية في مستخلصات بيض الضفدع الأفريقي المخلبي" . وقائع الأكاديمية الوطنية للعلوم في الولايات المتحدة الأمريكية . 100 (3): 975-980 . doi : 10.1073/pnas.0235349100 . ISSN 0027-8424 . PMC 298711. PMID 12509509 .   
  4. غريغور ف. فوسمان، ستيفن ب. إلنر، كايل و. شيرتزر، ونيلسون ج. هيرستون الابن. عبور تشعب هوبف في نظام مفترس-فريسة حي. مجلة ساينس . 17 نوفمبر 2000: 290 (5495)، 1358-1360. doi : 10.1126/science.290.5495.1358
  5. سونغ سي، فينيكس إتش، عابدي في، سكوت إم، إنغالز بي بي، وآخرون. 2010. تقدير بنية التشعب العشوائي للشبكات الخلوية. PLoS Comput Biol 6(3): e1000699. doi : 10.1371/journal.pcbi.1000699
  6. ^ أرطغرل م. أوزبوداك، موكوند ثاتاي، هان إن ليم، بوريس آي شريمان وألكسندر فان أوديناردن. تعدد الاستقرار في شبكة استخدام اللاكتوز في الإشريكية القولونية. طبيعة . 19 فبراير 2004؛ 427(6976):737–40