[论文解读] Modeling network dynamics: the lac operon, a case study

José M. G. Vilar, Călin C. Guet+1 more|ArXiv.org|Nov 17, 2004

Gene Regulatory Network AnalysisBiochemistry, Genetics and Molecular Biology参考文献 16被引用 17

一句话总结

本文提出了一种大肠杆菌乳糖操纵子的多尺度模型，整合了分子、细胞和种群水平的动力学，以模拟基因调控。通过结合随机基因表达、代谢反馈和种群异质性，该模型再现了关键的实验行为，如双稳态、诱导阈值和种群水平的双峰性，展示了系统级建模在理解基因网络方面的强大能力。

ABSTRACT

We use the lac operon in Escherichia coli as a prototype system to illustrate the current state, applicability, and limitations of modeling the dynamics of cellular networks. We integrate three different levels of description -molecular, cellular, and that of cell population- into a single model, which seems to capture many experimental aspects of the system.

研究动机与目标

• 开发一个统一分子、细胞和种群水平动力学的综合模型，以描述乳糖操纵子系统。
• 研究随机基因表达和代谢反馈如何塑造细菌种群的表型异质性。
• 利用单一集成模型框架再现实验中观察到的关键行为，如双稳态和诱导阈值。
• 评估当前建模方法在捕捉复杂网络动力学方面的局限性和适用性。
• 为一个广泛研究的生物网络中的基因调控提供系统级理解。

提出的方法

• 该模型整合了分子水平上的随机基因表达模拟，包括转录、翻译和乳糖阻遏蛋白结合。
• 通过乳糖代谢和cAMP-CRP激活引入代谢反馈，将基因表达与细胞能量状态联系起来。
• 通过追踪单个细胞随时间的变化来模拟细胞水平的动力学，考虑细胞生长、分裂和随机蛋白表达。
• 种群水平的行为通过模拟数千个单个细胞而自然涌现，从而允许对表型分布和双峰性进行分析。
• 该模型采用混合方法，对代谢通量使用微分方程，对随机基因表达使用Gillespie算法。
• 参数化基于文献中的实验数据，确保模型的生物学相关性和预测能力。

实验结果

研究问题

• RQ1 乳糖操纵子中的分子水平相互作用如何导致种群水平的异质性？
• RQ2 代谢反馈在塑造诱导阈值和双稳态行为中起什么作用？
• RQ3 一个单一的多尺度模型能否在不同生物尺度上再现关键的实验观察结果？
• RQ4 随机基因表达如何促进同源细菌种群中的表型多样性？
• RQ5 当前建模框架在捕捉基因网络完整动力学方面存在哪些局限性？

主要发现

• 该模型成功再现了实验中观察到的双稳态诱导行为，表现出非诱导态和诱导态之间尖锐的转变。
• 种群水平的模拟显示了基因表达的双峰分布，与实验的流式细胞术数据一致。
• 该模型预测，通过cAMP-CRP介导的代谢反馈显著调节诱导阈值和响应动力学。
• 基因表达中的随机性即使在同源种群中也会导致显著的表型异质性。
• 该模型揭示，分子噪声和代谢反馈对于捕捉实验行为的完整范围都至关重要。
• 将多个尺度整合到单一框架中，证明了多尺度建模在系统生物学中的可行性与必要性。