| 第670行: |
第670行: |
| | | | |
| | | | |
| − | 生物中的基因调控网络也可以用CTMC来描述,考虑一个两节点的基因调控网络,其中基因A有自发激活和降解,并激活调节子R;调节子R有自发降解,并抑制基因A。
| + | 在系统生物学中,基因调控网络(Gene Regulatory Network, GRN)描述了基因、蛋白质和其他分子之间复杂的相互作用,共同控制基因的表达水平。为了定量分析这类网络的动态行为,特别是涉及随机性的情况(如基因表达的噪声),常用CTMC作为数学建模的框架。考虑一个简化的两节点基因调控网络模型,包含以下两个组件。基因A:具有自发激活(Spontaneous Activation)和基础失活(Basal Inactivation)的能力。当基因A处于激活状态时,它能促进转录调节子R(Transcriptional Regulator R)的表达或激活。调节子R :具有基础失活(Basal Inactivation)的能力。当调节子R处于激活状态时,它会抑制基因A的表达。这个网络形成了一个负反馈回路:A激活R,而R抑制A。 |
| | + | [[文件:基因调控网络.png|缩略图|居中|两节点基因调控网络]] |
| | + | 根据基因A和调节子R的状态组合,整个系统可以处在一下四种离散状态之一:<math>S=\{A^-R^-,A^-R^+,A^+R^-,A^+R^+\}</math>。并对各状态转移参数进行建模,<math>k_1</math>是A的自发激活速率,<math>\gamma_1</math>是A的基础失活速率,<math>\gamma_2</math>是R介导的A抑制速率,<math>\lambda</math>是A激活R的速率,<math>\mu</math>是R的基础失活速率。下图展示了该基因调控网络四种状态之间可能的转移路径及其对应的速率参数: |
| | + | [[文件:基因调控网络状态转移示意图.png|缩略图|居中|基因调控网络状态转移示意图]] |
| | + | 基于上述状态定义和转移速率,可以构建描述该CTMC模型的转移速率矩阵Q。矩阵的行和列按顺序对应状态<math>\{A^-R^-,A^-R^+,A^+R^-,A^+R^+\}</math> |
| | + | |
| | + | <math>Q=\begin{pmatrix} |
| | + | -k_1 & 0 & k_1 & 0 \\ |
| | + | \mu & -\mu & 0 & 0 \\ |
| | + | \gamma_1 & 0 & -(\gamma_1+\mu) & \lambda \\ |
| | + | 0 & \gamma_2 & \mu & -(\gamma_2+\mu) |
| | + | \end{pmatrix}</math> |
| | + | |
| | + | 这个简单的两节点模型虽然结构简洁,通过调整参数的相对大小,能够展现出双稳态或者震荡的行为。同样的建模思路可以推广到更复杂的调控网络中,为研究基因表达噪声、细胞命运决策(如双稳态)、生物节律(如震荡)等生物学现象提供了强大的定量分析工具。 |
| | | | |
| | ===转移过程和演化性质=== | | ===转移过程和演化性质=== |