关于生物节律的数学模型

关于生物节律的数学模型

关于生物节律的数学模型

摘要：为了对人体的生命节律进行有效的解释，通过对人体生物钟的分析，得到生命节律的变化的实质就是生物钟的变化，而人体生物钟是在多个基因及蛋白的共同作用下进行调节的，为了将生物钟很好地运用数学语言进行表达，本文采用系统动力学方程模型对生物钟进行解释，其中，整个生命体就是系统，由于节律有振幅和频率的变化，因此，本文对其两种情况分别建立Hopf分叉判别模型和劳思—赫尔维茨判据对其分析，同时本文还描述了当振幅和频率都发生变化的时候的处理方法，即将Hopf分叉判别模型和劳思—赫尔维茨判据共同作用即可，为了很好的检验本文所建立的模型的准确性及实际运用意义，本文通过对心肌细胞进行分析处理，得到了心脏在用内在和外在的影响的情况下的振荡图形，与实际振荡基本一致，从而验证了该模型的准确性及实际意义。

关键字：系统动力学、Hopf分叉、劳思—赫尔维茨判据、混沌

一、问题的重述

自然界的生命体或者人工智能系统中，可以发现有许多周期或者近似周期节律的现象存在。这些周期节律的现象会根据自身内部的反馈机制以及外界的刺激或者外界输入的变化而变化：(A)或者节律变快了，但是振幅变化不太大；(B)或者振幅变化很大，但是节律相对稳定；(C)亦或者两者都会有较大的变化。例如，正常人的心跳节律在有外界刺激时就会有(A)类现象的出现；正常人和病人的体温昼夜变化又会有(B) 类现象的出现。

（1）请设计若干数学模型，以便能够充分反映上述描述的这些节律现象，并解释模型产生这些现象的机制。

（2）考虑将设计的模型具体运用于某一个具有实际背景的系统中，并评价模型的优点和缺点。

二、模型假设

三、符号说明

及

，

=

i

K

(等

3,2,1

clo

)

bm a

25

...

表示Michaelis常数

i

)

V

(=

i

3,2,1

25

...

表示各部分的最大速率

i

定义：小写rev（或reverbα）、bma（或bmall）、clo（或clock）、1

per、cry、2

2

per、crie、1

cry表示相应基因，而大写的上述字母则表示相应蛋白CB表示CLOCK蛋白与BMALL蛋白的异二聚体

四、问题分析

生物节律现象的各种变化，根本上是由生物钟决定的。生物钟是指激发生物体内的生物节律并使之稳定维持的内部定时机制。生命科学家发现，生物钟是多种多样的，就人体而言，已发现一百多种。人体的生理指标，如体温、血压、脉搏；人的体力、情绪、智力和妇女的月经周期；体内的信号，如脑电波、心电波、经络电位、体电磁场的变化等等，都会随着昼夜变化作各种周期性变化。

生物钟主要是由生物体内的一种被称作时钟基因的细胞核或细胞质决定的。生物钟基因普遍存在于生物体内，起着调控机体生命活动和近日节律运转的作用，它们及其编码的蛋白质产物构成的自主调控的转录和翻译负反馈环是生物钟运转的分子机制。在生物体内，存在着不止一个时钟基因，同时还存在着大量对生物体的时钟基因有影响的基因，时钟基因之间相互耦合或彼此之间相互影响会同时对生物节律变化起作用。由于生物体内的生物钟存在对光或其他外界噪声、湿度以及温度等的敏感性，所以生物体的时钟基因与外界环境存在着一定的关系。总而言之，生物节律是系统作用的结果。从系统科学的观点考虑，它是一个负反馈控制的系统。因此，可以运用系统动力学的有关方法建立数学模型，以反映题目所描述的节律现象。（意义？？？）

对于一般的生化系统来说，都存在着一个稳定态。本文要判断一个系统是否处于稳定状态，必须得进行稳定性分析。通过对该生命节律模型的稳定性分析，本文可以得到该系统的稳定态标准。

系统的震荡性是生物节律发生变化的基础。要产生稳定性震荡，首先得满足两个条件：1.系统具有自我增加的元素，2系统是具有负反馈的部分

相关临界值则反映了生命节律的稳定态和周期震荡，通过求出相关临界值，可以判断生命节律变化处于何种状态。

五、模型建立

生物节律的起始点是一种“转录—翻译—抑转录”机制构成的反馈环。反馈环的正性因素通过PAS 结构域相互配对，结合于生物钟基因上游的启动子上，启动基因转录，进而翻译产生生物钟蛋白，参与生物钟现象。而生物钟蛋白同时作为负性因素，进人细胞核通过与正性因素结合而抑制生物钟基因启动子的活化，使mRNA 水平下降，最终蛋白自身减少，反馈抑制作用逐渐减弱，完成一个循环。各生物钟基因相互作用机制如下图所示：

生物钟基因转录因子主要反馈环示意图

注：-表示负反馈，+表示正反馈， 表示自身震荡

19

K 18

K

per2

crie

clock

per1

cry2

cry1

21

21,V K

PER1

16

K CRY2

25

25,V K

44,V K 33,V K

55,V K 66,V K 7

7,V K

PER2

17

K

22

22,V K

CRY1 24

24,V K

20

K

CRIE

23

23,V K

1

1,V K

14

14,V K

1313,V K CLOCK

BMALL

11

11,V K

12

12,V K 15

15,V K 10

K

REVBRB α

reverb α

2

2,V K

9

9,V K

bmall

8

8,V K