数理解析研究所講究録984 巻 1997 年 16-3016周期解の発生と概日リズムのモデル
鈴木

保夫 (YASUO SUZUKI)

大阪府立美木多高等学校

大阪府堺市城山台 4–1–1
大阪府立大学工学部数理工学科にて、 研修中

ABSTRACT. Mathematical models,which are based on a understanding of biochemical process and biological knowl\’ege,can provide nonintuitive insights into
the circadian rhythmo and can suggest new avenues for experimentation.A simple
mathematical approach is developed to explore the dynamics of the Drosophila cir-protein andcadian components,that are per-RN
PER-TIM complex.The model simulations are compared to data on the generation of circadian rhythms,entrainment to LD cycles,phase response curves (PRCs),lossof rhythmicity in constant light, per mutations and temperature compensation.These
results indicate that the model is sufficient to exhibit the Drosophila circadian$\mathrm{A},PER-pr\mathrm{o}te\mathrm{i}\mathrm{n},tim-\mathrm{R}\mathrm{N}A,TIM$

clock.

1. 序

1.1. 生物のいろいろな周期現象
生物は、様々な周期を持っている。 その周期の発生原因によって、 大きく 2 つに分けることが出来る。 第 1 は、 生物個体が外界の周期現象 (例えば昼夜とか潮汐など) に、 同調して周期が発生する場合である。 第 2 は、 生物個体自身の中から自律的に周期が、発生している場合である。 これらは、 恒常条件に隔離した条件のもとで、 周期の発生を調べれば、 区別することが出来る。

自律的な周期の例として、神経インパルス、昆虫のはねの振動数心臓の拍動、歩) ホルモンの周
行や動作のリズム、 呼吸のリズム、 概日リズム (circadian期、 個体数の変動、 生理的な寿命などがあげられる。 自律的な周期現象の発生原因$\mathrm{r}\mathrm{h}\mathrm{y}\mathrm{t}\mathrm{h}\mathrm{m}\mathrm{s}$$\text{、}$は、 昆虫のはねの振動数や歩行や動作のリズムのように力学的な固有振動数に関係する場合や神経インパルス、 心臓の拍動、 呼吸のリズム、概日リズム、 ホルモンの周期のように、 自己調節系から周期が出現する場合もある。 自己調節系の周期現象

は、一般に自励系の微分方程式の周期解で表現される。以下では、ショウジョウバエの概日リズムの数理モデルについて、数学的な準備およびモデルの ミユレーショ
ンと事実との適合性について述べる。 なお、 これらの議論の多くは、 Runge-Kutta$\sqrt[\backslash ]{}$

17

4 次近似による数値解析によって求めた結果に基づいている。

2. 数理モデルを作成する上での視点
数理モデルを作る場合、 いろいろなアプローチがあるけれども、 今回の場合は、現在までの研究によってわかっていることをもとにしてモデルを作ろうとしている。
そのとき、分かっていない事柄を数理モデルの上でどう表現するか、 そして観測された量の極大の時刻間のずれの原因を、 モデル式の中でどう考えるかが、 モデル作
成上のポイントとなっている。 まず、隔日リズムを作っている化学物質の量の変化を微分方程式で表現する場合の検討をしてみよう。

2.1. 時間遅れを含まないモデル
1. 自励系、 1 元、 1 階、 時間遅れを含まない常微分方程式
最も簡単な例から考えることにする。 下図のように、 が I の速度で増加し、$S_{1}$

速度で減少している。 そして、

$S_{1}$の量に応じた negative feedback を I が受けて量を調節している。 negative feedback$\text{の式を}\frac{\beta}{1+\alpha S_{1}^{\rho}}$とする。 このとき$S_{1}$$k_{1}$の$S_{1}$のの量を表現するモデル式は、 次のようになる。 この式は、周期解を持たないことが知られている。

$\dot{S}_{1}=\frac{I\beta}{1+\alpha S_{1}^{\rho}}-k1S_{1}$2. 自励系、 多元、 1 階、 時間遅れを含まない常微分方程式
化学物質が逐次反応している場合を考える。反応系の最後の Sn の量に応じた negative
feedback を が受けて調節している。 このときの量を表現するモデ
ル式は、 次のようになる。 この式は、 周期解を持つ場合がある [1]。
$I$$\dot{S}_{1}=\frac{I\beta}{1+\alpha S_{n}^{\rho}}-k_{1}s_{1}$....... .. . . . . . . . .. ..,$S_{1},$ $S_{2},$,$\dot{S}_{2}=k_{1}S1^{-}k2S2$,$\dot{S}_{n}=k_{n-1}Sn-1^{-}kns_{n}$.

$\cdots,$$S_{n}$

18

22. 時間遅れを導入したモデル

1. 自励系、 1 元、 1 階、 時間遅れを含む常微分方程式
だけ遅れて、 negative feedback が働く場合を考える。 そのとき、 negative
時間$\tau$

feedback の式を$\frac{\beta}{1+.\alpha S_{1}(t-\mathcal{T})^{\rho}}$とするとよい。 このとき、 周期解を持つ場合があ.

る [2]。

$\beta$$1+\alpha s_{1}(t-\mathcal{T})^{\rho}$

En$\epsilon$$I\ovalbox{\tt\small REJECT}_{s_{\mathit{1}}}\mathrm{f}\mathrm{e}\mathrm{e}\mathrm{d}kJ\mathrm{b}\mathrm{a}\mathrm{c}\mathrm{k}t\mathrm{s}\tau$$\dot{S}_{1}(t)=\frac{I\beta}{1+\alpha S_{1}(i-\mathcal{T})^{\rho}}-k1S1(t)$2. 自励系、 多元、 1 階、 時間遅れを含む常微分方程式
物質や波の移動には、 必ず時間がかかるので、 同–の空間に存在しないかぎり 「自然現象は、 time delay を含む」 と考えるのが自然である。 移動に要する時間は、化学反応では短くて無視される場合が多いけれども、無視できない現象もあるだろう。
もちろん、 この式は周期解をもつ場合がある。

,$\dot{S}_{1}(t)=\frac{I\beta}{1+\alpha S_{n}(t-\mathcal{T}_{n})^{\rho}}-k_{1}S_{1}(t)$.....................

$*\cdots..-\ldots\ldots\ldots..\cdot.$,$\dot{S}_{2}(i)=k_{1}s_{1}(t-\tau 1)-k2S_{2}(\mathrm{t})$,$\dot{S}_{n}(t)--k_{n-}1s_{n-}1(t-\mathcal{T}_{n-}1)-knsn(t)$.

23. 周期解の発生と 2 つの time delay
自励系、 1 階の常微分方程式について、 1 元、時間遅れのない常微分方程式は周期解を持たないけれども、 –元、 時間遅れを持つ常微分方程式および多元の常微分方
程式は周期解を持つ場合がある。 2 階以上の場合は、 1 階の多元の方程式に変換することができる。 このことから、 周期解の発生には、 時間遅れの要素が関係していることがうかがわれる。
微分方程式に delay parts として、直接表現されている time delay を explicit timedelay、式の形や係数の値によって生じる極値の時刻の遅れをうみだす原因となる過去の影響を、 implicit time delay と呼ぶことにする。 たとえば、

19$x’(t)= \frac{\beta}{1+\alpha y(t)^{\rho}}-_{X}(t)$初期値$=0,t_{0}=0$.

$y’(i)=X(t)-y(t)$ .

のとき、 定数変化法により$x’(t)=. \frac{\beta}{1+\alpha(e^{-t}\int_{0^{t}}e\mathit{8}x(_{S})ds)^{\rho}}-x(t)$ここで、積分の項が implicit time delay を表す。 この式は、周期解を持っ場合がある [1]。 したがって、 自励系、 1 階、 1 元の常微分方程式であっても、現在の自己が過去の自己の影響 (積分の項、 あるいな delay parts) を受けるとき、周期解が生じる場

合があると考えられる。言いかえれば、「自励系の常微分方程式が、周期解を発生していれば、広義の time delay を含む閉ループが存在する」 と考えられる。

さて、モデル作成上では、物質の移動時間による遅れは、式の中に直接 delay partsとして、 explicit time delay によって直接表現されるべきものであり、 反応時間による遅れは、 式の形や係数の値によって表現されるべきものである。 explicit timedelay と implicit time delay の働きにより、極値をとる時刻に遅れが生じる。たとえば、 次のような系を、微分方程式で表現すると (各ステップは、 時間 の遅れをもち、 速さ で変化する)$\tau$$\mathrm{k}$入力$5+ \mathit{5}sin(\frac{2\pi}{24}t)\frac{\backslash }{k\prime}$$s_{\mathit{1}}$$arrow k\tau$$s_{\mathit{2}}$$arrow k\tau$$\dot{S}_{2}(t)=k(S_{1}(t-\tau)-S_{2}(t))$,$\dot{S}_{8}(t)=k(S7(t-\mathcal{T})-S8(t))$.

$. \dot{S}_{1}..(t.)..=.k..(.5..+.5.s.i.n..(..\frac{2\pi}{24}.l.)...-..S_{1}..(t.).\mathrm{I}..’,$となる。 このとき、 極値をとる時刻の遅れを、解曲線から求めると下の表のようになる。 この単純な例では、 $\tau=0$ から $\tau=1$ になると遅れは、 それぞれする。 また、$\mathrm{k}$が 2 倍になると遅れは、 ほぼ

緬確か $y,$$\lambda$$\frac{1}{2}$$\tau$だけ増大になる。

賎勿 [m\urcorner ‘屋勲

implicit time delay の代わりに、 delay parts を使って、 explicit time delay によって表現することはできないけれども、未知の調節や反応による implicit time delayに相当する値を、 delay parts に加えることによって、 同–の周期をもつ式を作ることはできる。 したがって、未知の調節や反応を表現する式の代わりに、 delay parts

20を変化させて同–の周期を出現させることにより、 black box を含んでいる現象を暫定的に表現できる。 したがって、 black box をふくむモデルの delay paxts は、未知の調笛や厘応に上る imnlicit time delav と分子の移動時間による
time$\mathrm{e}\mathrm{X}\mathrm{D}\mathrm{l}\mathrm{i}\mathrm{C}\mathrm{i}\mathrm{t}$モテルを作成し、 適当な係数を与えて$\backslash \nearrow\backslash$

ミユレーションを行い、 既知の実駅結呆との比較検討をしよう。

3. 概日リズムの数理モデル
3.1. 概日リズムの数理モデル
概日リズムは、 単細胞生物からヒトまでの様々な系統の生物にみられる。 1 個の細胞内でも概日リズムが働くことから、 細胞間の相互作用は主要なことではないと考えられる [ 。また、概日リズムが多くの真核生物にみられ、原核生物にはあまり$3|$みられないことは、真核生物のみが持っている核膜や細胞内膜系との関係がうかがわれた。 しかし、藍藻の–種である Synechococcus sp(原核生物である) は、概日リズムと温度補償性を確実に持っていることがわかっている [ 。けれども、 まだ真核$4|$生物の時計と原核生物の時計とが、等しいかどうかわかっていない。 ここでは、ショウジョウバエの概日リズムのモデルを考えるが、原核生物や他の生物の概容リズムについても示唆することがあるだろう。

21

ショウジョウバエの概日リズムに関する最近の研究成果は、何がリズムの中心的役割をはたしているのかに迫っていると思われる。 1996 年に別々の 4 チームによって、 よく似た内容のリズムの位相変位 (light-entrainment mechanism) のモデルが提案されている [5, 6, 7, 。これらの位相変位メカニズムの中には、 リズム発生のメ
カニズムが含まれていると考えられる。 生命現象は、 分子の言葉で表現されるもの
であることは言うまでもないことであるが、 分子の反応順序や周期的な分子の濃度変化の記述のみでは、 リズム発生のメカニズムの予想は立てられても、 実際に 24時間の周期をもつリズムを起こすことを示せないであろう。 リズムの発生は、 DNA
転写調節に関する negative feedback によって引き起こされることがすでに提案されている $[9, 10]$ 。しかし、 フィードバックから必ず、 リズムが発生するとは限らない
ので、 リズムが発生することを示すためには、分子生物学の成果に基づく数理モデ
ルが作られることが必要である。 また、そのことによってはじめて、周期の突然変$8|$異や位相変位や温度補償性についての詳細な検討をなすことができると考えられる。
生物時計を構成している遺伝子から作られる Per-RNA と PER-protein は,、 DNA
転写調節に negative feedback loop をもつことによって、 24 時間の周期をもつと考per-RNA は核内でつくられ、 RNA processing control,RNAtransport control,translation control という様々な調節 $[15, 16]$ を受けたのち、核外にでて PER-protein を合成する。このとき、 per-RNA の極大の時刻と PER-proteinの極大の時刻との間には、約 6 時間 (少なくとも 4 時間) の遅れをもっことがしられている $[3, 5]$ 。もう –つの生物時計を構成している遺伝子から作られる tim-RNA も、同様にして核外にでて、 $TIM$-protein を合成する。 つぎに、 PER-protein や TIMprotein は、 単独では核膜を通過できず、 PER-TIM complex になって核膜を通過し、核内に入って DNA 転写調節に negative feedback として働いて、 RNA を減少させる。 PER-TIM complex の核膜の通過には、 調節があることがしられている[5, 6, 7, 8, 11, 12]。さらに、 これらのタンパク質はリン酸化や他の過程によって活えられている $[9, 10]$。

性を調節されているらしい $[13, 14]$ 。また、 $TIM$-protein が光によって減少すること、そして光による $TIM$-protein の減少が位相の変化をもたらしていると考えられている [5, 6, 7, 8]。さらに、 $TIM$-protein と PER-protein はよく似た位相をもつ量的変化をすることがわかっている [

$5|$

。

これらの知見をもとにして、Per-RNA.

-RNA.

$PER- \mathrm{p}\mathrm{r}\mathrm{o}\mathrm{t}\mathrm{e}\mathrm{i}\mathrm{n}_{\text{、}}tim$$TIM- \mathrm{p}\mathrm{r}\mathrm{o}\mathrm{t}\mathrm{e}\mathrm{i}\mathrm{n}_{\text{、}}$PER-TIM complex の各分子の量 (同じ意味で濃度あるいは活性) の変化を数理モデルで表現しよう。 そのとき、未知の調節や反応によって起こされる time delay お
よび、分子の移動時間による time delay を考慮して、時間遅れを含む微分方程式によって表現することにする。 DNA 転写調節の詳細は、知られていないので、 negativeを用いる。 この式を含む方程式には、周期解を
feedback を表す式として、$\underline{\beta_{ }}$$1+\alpha x^{\rho}$持つ場合があることが知られている $[1, 2]$ 。概日リズムのような恒常条件での自律的22$\dot{S}_{1}(t)=\frac{I_{1}\sqrt 1}{1+\alpha_{1}S5(t-\mathcal{T}_{3})^{\rho_{1}}}-k_{1}S_{1}(t)$$\dot{S}_{2}(t)=I_{27}kS1(t-\tau_{1})-k2S_{2}(t)-k6s2(t)S4(t)$$\dot{S}_{3}(t)=\frac{I_{3}\ }{1+\alpha_{2}S5(t-\mathcal{T}3)\rho 2}-k_{3}s_{3(t})$$\dot{S}_{4}(t)=I4k_{8}S_{3}(t-\tau 2)-k_{4}S_{4}(t)-k_{6}s2(t)s_{4}(t)$$\dot{S}_{5}(t)=k_{6}S_{2}(t)S_{4}(t)-k_{5}s_{5}(t)$Figl. ショウジョウバエの概日リズムの数理モデル
ドットは、 時間微分 ; は時刻 ; $S(t-\tau)$ は、 時刻 $(t-\tau)$ の
$t$

$S$

の値という意味 ;はそれぞれ per-RNA,PER-protein,tim-RNA,TIM-protein,PERは、 per-RNA が、$TIM$ complex の量 (同じ意味で濃度、 あるいは活性) を表す。
$S_{1},$ $S_{2},$$s3,$ $S4,$ $s_{5}$$\tau_{1}$核から細胞質へ未知の反応や調節を経て移動する時間を表現する time
は、 tim-RNA が、核から細胞質へ未知の反応や調節を経て移動する時間を表現する time乃は、 PER-TIM complex が、 細胞質から核への移動するのに使$\mathrm{d}\mathrm{e}\mathrm{l}\mathrm{a}\mathrm{y}_{\text{、}}\tau_{2}$$\mathrm{d}\mathrm{e}\mathrm{l}\mathrm{a}\mathrm{y}_{\text{、}}$

への転
から
は
われる時間を表現する time delay を表している。
濃度に比例して、 タンパク質が合成されるときの反応
は、
写速度、$I_{1},I_{3}$$\mathrm{R}\mathrm{N}\mathrm{A}$$I_{2}k_{7},$ $I_{4}k_{8}$

速度、$k_{1},$$k_{2},$$\mathrm{R}\mathrm{N}\mathrm{A}$$\mathrm{D}\mathrm{N}\mathrm{A}$$k3,$ $k_{4},$ $k_{5}$は、 基質が別のものになる反応速度、$k_{6}$

は、 $TIM$ と PER から PER-TIM complex ができるときの反応速度を表している。 周期解を持つとき、

周= 関数 ($\tau_{1},$$\tau_{2,3}\mathcal{T},$Il, $I3,$$I_{2}k7,$$I_{4}k_{8},$ $k_{1},$$k2,$ $k3,$ $k4,$ $k5,$ $k6,$ $\alpha 1,$$\alpha_{2},$$\beta 1,$oe,$\rho 1,$$\rho 2$) である。

23

なリズムの発生は、数学でいえば、初期関数が–定のとき、 自励系の微分方程式が、周期解を発生することに相当する。 周期解の発生や極大値の時刻の遅れには、 上記
2 つの time delay (2.3 をみよ) の和の大きさが強く関係していることが、数値解析,
によって予測されている。

以上のことより、ショウジョウバエの概日リズムの数理モデルを提案する (Figl)。
は、$S_{5}$

が移動に要する時間を考慮して、乃時間前の$\text{を受けて}\frac{I_{1}\beta_{1}}{1+\alpha_{1}s_{5}(t-\tau 3)\rho_{1}}$の速度で増加し、$k_{1}$$S_{5}$$S_{1}$の量に対して negative feedbackの速度で減少する。

$S_{\mathit{2}}$は、 タンパク

合成開始までに要する時間を考慮して、 \tau 1 時間前の の量に比例したの速度で増加し、 の速度で減少する。 は、 が移動に要する時間を考慮して、$S_{1}$$k_{2}$

時間前の層し、$\ovalbox{\tt\small REJECT}$$k_{3}$$S_{5}$

は、$k_{6}$$S_{5}$の量に対して negative feedbackの速度で減少する。

時間前の$S_{5}$$S_{3}$$S_{3}$$S_{\mathit{4}}$$S_{2}$と$S_{\mathit{4}}$$\text{を受けて}\frac{I_{3}\ }{1+\alpha_{2}S_{5}(t-\tau_{\mathit{3}})\rho 2}$$\tau_{3}$の速度で増は、 タンパク合成開始までに要する時間を考慮して、) の速度で増加し、 の速度で減少する。
より合成されて増加し、 の速度で減少する。

の量に比例したの速度で$I_{2}k_{7}S_{1}(t-\mathcal{T}_{1})$$I_{\mathit{4}}k_{\mathit{8}}S_{3(-\tau}t2$$k_{3}$$k_{5}$

3.2. 周期を決定する要因
この数理モデルの周期は、$\beta_{2},p_{1},\rho_{2}$$\tau_{1},$$\tau_{2,3,}\mathcal{T}$I,$I_{3,2}Ik_{\mathit{7}},$$I_{4}k\mathit{8},$$k1,$ $k_{2},$ $k3,$$k_{\mathit{4}},$$k5,$ $k6,$ $\alpha_{1},\alpha_{2},\beta 1$,の関数となっている。 しかし、数学的には、 この周期を求めることは困難であり、 まだ求められていない。 また、 モデル式が周期解を持つことは、 数学的には証明されていないが、 数値解析によって予想される。 この数理モデルの周期が変化するためには、 フィードバックの式が変化する場合
は、 2 種類の意味を含む time delay が変化する場合度が変化する場合$(\alpha_{1}, \alpha_{\mathit{2}}, \beta_{1}, \alpha, p1,p2)_{\text{、}}$$(_{\mathcal{T}_{1},\tau_{2}}$,$\tau_{3}^{)_{\text{、}}}$

あるい

あるいは、 反応速のいずれかである。 周期解を持つ範囲で、 \tau 1,\tau 2,\tau 3 の増加は、 周期の増加をもたらす。 \alpha 1, \alpha 2 の増加は、 周期の$(I_{1}, I_{3}, I_{2}k_{7}, I_{\mathit{4}}k\mathit{8}, k1, k_{2}, k3, k_{4}, k_{5}, k6)$

増加をもたらす場合が多い。
,12k7,I4k8,kl,k2,k3,k4,k5 の増加は、 周期の減少をもたらす場合が多い。 したがって、概日リズムの周期が変化した突然変異であ$\beta_{1},\beta_{2},I_{1},I_{3}$る $Per^{0S},per,per^{L},time\iota esS$ は、係数値のいずれかが、 野生種と異なる値となる変化が起こったと考えられる。 さらに、 一般に温度の上昇に伴って反応速度が増加する
ので、 そのとき周期を–定にする調節がないと周期は短くなる。 概日リズムは、 温

度変化にかかわらずに周期を–定にするという温度補償性を持っている。 したがって、 温度の上昇したときにはリズムが長くなるように、 上記の 3 点のうちのいずれかを調節していると考えられる。
と ; RNA の移動時間、 RNA プロセッ
また、 DNA 転写調節には、シング調節、 RNA 輸送調節と翻訳の調節には、 と ; -RNA の分解調節には、と; タンパク質の活性調節には、と; タンパク質の移動時間と輸送調節; 以上のように、各係数の値を変化させることによって、様々な調節をシには、$\alpha_{1},$$\alpha_{2},$$\beta_{1},$$\beta_{2},$$\rho 1$$\rho_{2}$$\tau_{1}$$k_{\mathit{3}}$$k_{\mathit{2}},$$\tau_{3}$

ミ $=$ レーションすることができる。

$k_{\mathit{4}}$$k_{5}$$\tau_{2}$$\mathrm{m}$$k_{1}$

24

4. モデルのシミ$=$

レーション結果と検討

41. 24 時間周期の発生
PER-protein の量は、 PER-antigen antibody を用いて測定されている。 PERantigen antibody は、 単独の PER-protei ばかりでなく PER-TIM complex の中の PER-protein にも反応している [6] ので、 これまでに報告されている PERprotein の量は、 PER-protein と PER-TIM complex を合計した値と考えられる。
したがって、 per-RNA の peak の時刻と PER-protein と PER-TIM complex の和の peak の時刻との差が約 6 時間であることを満たすようにとし、 周期. .
が 24 時間になるように他の係数の値を定める (Fi $g2$ ) 。この結果から、核膜を隔てた分子の移動および未知の反応や調節 [11, 12, 15, 16] によって作られる大きな timeを含む、 DNA 転写調節の negative feedback によって、 24 時間の周
delay
$n$

$\tau_{1}=\tau_{2}=4.0$$(\tau_{1}, \mathcal{T}_{2}, \tau_{\mathit{3}})$

期を引き起こすことができることが示された (Fig2AB)。いうまでもなく周期が 24
時間になる係数の値の組は無数にあるが、 進化を通じて適した値の組が、 自然選択された結果、実際の係数の値の組を持つ概日リズムが作られたと考えられる。 また、
レーションの per-gene の feedback の働きと $tim$-gene の feedback の働
このシミ
きは異なっているので、 単独の PER-protein と単独の $TIM$-protein の peak の時$\supset-$

刻は異なっている。

42.

$\mathrm{L}\mathrm{D}$(Light-Dark)cycle への同調

光によって、 $TIM$-protein が減少することは、$k_{4},$$k_{5}$が、 増大することによって

シミ $=$ レーションできる。 恒暗条件では、 $TIM$-protein は、$\mathrm{Z}\mathrm{T}20$($\mathrm{Z}e\mathrm{i}\mathrm{t}\mathrm{g}e\mathrm{b}\mathrm{e}\mathrm{r}$time)のとき、 1 時間で 15 倍に増加し、 $\text{光_{の}照射^{によ}って減少して}\frac{1}{6}$ に漸近することが知られている [5]。

このことから、 S4 の増加速度を 1 、減少速度を$\dot{S}_{\mathit{4}}(t)=1-as4(t)$. より

したがって、 $S_{4}(0)=1$ とすると、$S_{4}(1)=(1- \frac{1}{a})e-a+\frac{1}{a}=1.5$$-..—\mathrm{A}\cdot\backslash --$

’

.

$-$

光の照射により、

1

$a$

とすると、

$S_{\mathit{4}}(t)=e_{0}e-at+\underline{1}$1 時間で 15 倍になる恒暗条件の時は、これより、 $a=.0.4$

6 に漸近するときは、$\backslash ---\tau-,$$-\mathrm{Y}$.

-L.

$\tau$.

$-’\backslash$

$-$

1111

これより、

$\lim_{larrow\infty \mathit{4}}S(t)=-a\mathrm{A}=\overline{6}\perp$$a=6$

モデルでは、 $k_{5}=0.3$ としたので、 光の効果は、

PER-TIM-comPlex ので、$\dot{S}_{2}(t)$$k_{4},$$k_{5}$

を 20 倍にすることおよび-protein がなくなると、 単独の PER-protein が残るの$\dot{T}IM$式の右辺に、 $(20-1)\cross 0.3S_{5}(t)$ を加えることによってシミ$\supset-$

レーション

することができる。

その結果、各成分は、 L12D1\nu cyc1e(明期 12 時間、暗期 12 時間) に同調する。 さ25らに、 per-RNA,PER-protein,tim-RNA,TIM-protein,PER-TIM complex の振幅は、 恒暗条件より大きくなり、 かつ極大値と極小値の比は、 大きくなる$(\mathrm{F}\mathrm{i}\mathrm{g}2\mathrm{B}\mathrm{C})$。

-cycle では、直航条件と比較して、 リズムが、 より強く明確におこることを意味している。 また、 単独の PER-protein の極大の時刻は、 ZT1 と) 。同様に、他のZT12 となると予想される (Ficycle においても、各成分
は、 LD–cycle に同調し、 また、 単独の PER-protein の極大の時刻は、 光を消した
時刻と光をつけ始めてから約 1 時間後の時刻に極大があらわれる場合が多い。 ただし、光の効果を 20 倍から順に小さくするにしたがって、つまり明期の光を弱くする
このことは、$\mathrm{L}12\mathrm{D}12$$\mathrm{L}\mathrm{D}$$g2\mathrm{C}$にしたがって同調しない状態が出現してくる。

43. 光の影響による位相の変化 (PRCs)
恒暗条件の下で、 circadian time ごとに、 一定量の光を与えたとき、 その後の周期が進んだり遅れたりすることが知られている [17]。弱い光の効果と強い光の効果
を、 $TIM$-protein の減少速度をそれぞれ 4 倍、 あるいは 20 倍することによって、シミ $=$ レーションすることができる。 その結果、 光照射を与える時刻によって、 その後の位相が進んだり、 遅れたりする
る 2 番目の$(\mathrm{F}\mathrm{i}\mathrm{g}3\mathrm{A}\mathrm{B})$。次に、刺激終了後 13 時間以降にあの極大の時刻によって、位相のずれを調べ、位相反応曲線 (PRCs) を求めた。 このとき、 弱い光では 1 型、 強い光では、 型になった。実$S_{5}$$0$$(\mathrm{F}\mathrm{i}\mathrm{g}3\mathrm{C}\mathrm{D})[17]$際のショウジョウバエの羽化や活動量の PRCs に大変よく似た曲線を描くことがわかった。 また、 このシミ $=$ レーションでは、刺激終了時刻が、$\mathrm{t}=86$, 光の効果を 62倍にすると、 リズムが弱くなり、 回復するのに 100 時間ほどかかる消去点に近い値となる。 適当な位相に型を起こす刺激より弱い適当な大きさの刺激を与えると、リズムが消滅することが知られている [17]。以上のことからこのモデルは、 光による位相の変化を反映していることがわかる。
$0$44. 恒明条件でのリズムの消失
野生種の PER-protein の量は、 恒暗条件では変動せずリズムを持たないことが知られている [18]。恒明条件での各成分の変化を、 光の効果を$\mathrm{F}\mathrm{i}g2\mathrm{C}$と同様に変化させてシミ $=$ レーションした。 そのとき、光の効果を 20 倍にすると、時計の各成分
は、 リズムを持たなかった$(\mathrm{F}\mathrm{i}\mathrm{g}4\mathrm{A})$。次に光の効果をしだいに弱くするにしたがって、周期と振幅がどう変化するかを調べたより大きくなると、振幅が$0$$(\mathrm{F}\mathrm{i}\mathrm{g}4\mathrm{B})$。そのとき、光の効果が、 87 倍となりリズムが失われることが分かった。 また、光の効果が、 87 倍より小さくなるにしたがって振幅が大となり、周期も大となる。 このことは、一定の強さの以上の光のもとで、 リズムが消滅することを示している。 また、 一定の強さ以下の光のもとでは、 恒暗条件と比較して、 振幅が小さいので、 弱いリズムが発現すること、 および 短い周期となることが予想される。

26

A

Fig2 シ ュレーションの結果
Figl の式に次の係数の値をもちいる。
$\approx\sim$$\tau_{1}=\tau_{2}=4.\mathrm{o},$$k_{4}=0.20,$$A:S_{5}- t$$\mathrm{B}:S_{1},$$\tau_{3}=1.0,$$k_{\mathit{5}}=0.30,$$k_{6}=0.20,$$I_{\mathit{3}}=2.0,$$\alpha 1=1.0$$I_{2}k_{7}=I_{4}k_{8}=1.\mathrm{o},$$k_{1}=0.38,$, \alpha 2=10,\beta 1=[ら $=1.0,$の初期関数 $=0$ ,$k_{2}=0.38,$$k_{3}=0.30$,$\rho_{1}=\rho_{2}=3.0$

恒暗条件。 24 時間周期が発生する。
曲線恒暗剣牛。tim-RNA の\Phi f[の時刻を CT12 (circadian time)

曲線$S2,$ $S_{3},$$I_{1}=1.0,$$S_{1}(t),$ $S_{2}(t),$ $s_{3(}t),$ $S\mathit{4}(t),$ $s_{5}(t)$$S_{\mathit{4}},$$S_{\mathit{5}}-.t$

ときめる。

, ,
時刻は、

$\mathrm{C}:S_{1}$$S_{2}$$S_{3},$曲線。 (明期):D(暗期) $=12:12$ 。単独の PER-protein の極大のと ZT12 となる。 $24n\leq t<12+24n$ のとき、 明期、 光の効果とし$S_{\mathit{4}},$$\mathrm{Z}\mathrm{T}1$$S_{5}-t$$\mathrm{L}$て、 $k_{\mathit{4}}=20\cross 0.2,k_{5}=20\cross 0.3;\dot{S}_{2}(\iota)=S_{1}(t-4)-\mathrm{o}.38s2(t)-\mathrm{o}.2s_{2}(\iota)S\mathit{4}(\iota)+(20-$$1)\cross 0.3S\mathit{5}(t)$にする。 $12+24n\leq t<24+24n$ のとき、 暗期、 Fig2 の係数値を使う。 $(n=0,1,2\cdots\cdot)$

27

Fig3. 位相反応曲線 (PRCs)
曲線位相が遅れる。実線は、恒暗条件の時で、$t<80$ に光を与えている。光の効果は、 $78\leq t<80$ のとき$\mathrm{A}:S_{5}- t$のとき$\mathrm{F}\mathrm{i}\mathrm{g}2\mathrm{A}$$k_{\mathit{4}}=4\cross 0.2$と等しい。点線は、 $79\leq$$k_{\mathit{5}}=4\mathrm{X}0.3_{\text{、}}79\leq t<80$と $\dot{S}_{2}(t)=S1(t-4)-0.\mathrm{s}8S2(t)-0.2s2(t)S4(t)+(4-1)\cross 0.3s_{5}(t)$ .

その他の時刻では、 Fig2 の値をとる。

$k_{5}$

を 2 時間変化させることは、 PER-TIM

complex が、核内に入るときの time delay 1 時間を考慮すると、実際には、光を 1
時間あてたことに相当する。
曲線 位相が進む場合。実線は、恒暗条件の時で、 Fiと等しい。点線は、$\mathrm{B}:S_{5}- t$$g2\mathrm{A}$$89\leq t<90$ に光を与えている。

C:弱い光をあてたときの$\mathrm{P}\mathrm{R}\mathrm{C}_{0}$タイプ 1 となる。 光の効果を、 4 倍としたとき。 刺激終了後 13 時間以降にある 2 番目の の極大の時刻を比較することによって位相のずれを調べ、
を描いた。
の場合に
は,t
は,t
あたる。位相の遅れは負、 進みは正であらわす。 Circadian time の から 12 は主
観的昼をあらわし、 CT12 から CT24 は、 主観的夜をあらわす。
タイプ O となる。光の効果を、 20 倍としたとき。他
D:強い光を当てたときの$S_{\mathit{5}}$$\mathrm{P}\mathrm{R}\mathrm{C}$$\mathrm{F}\mathrm{i}\mathrm{g}3\mathrm{A}$$=80,\mathrm{C}\mathrm{T}\mathrm{l}\mathrm{l}.\mathrm{F}\mathrm{i}\mathrm{g}3\mathrm{B}$$=90,\mathrm{C}\mathrm{T}21$$0$$\mathrm{P}\mathrm{R}\mathrm{C}_{\mathrm{o}}$は、$\mathrm{C}$

と等しい。

28

Fig4:恒明条件における光の効果と周期や振幅との関係。
の曲線 光の効果を、 20 倍にする。
, 恒強光下での
4 $-\mathit{0}.38s_{2}(t)-0.2S_{2}(t)S\mathit{4}(t)+(20-1)\mathrm{X}0.3S5(t)$の曲線 光の効果は、 4 倍にする。
, 恒弱光下での$\mathrm{A}:a$$\dot{S}_{2}(t)=S_{1}(t$ -$S_{5^{-}}t$$k_{4}=20\cross 0.2_{\text{、}}k_{5}=20\cross \mathit{0}.3_{\text{、}}$$S_{5^{-}}t$$k_{4}=4\cross 0.2_{\text{、}}k_{5}=4\cross 0.3_{\text{、}}$

$)$

$b$

。

$\dot{S}_{2}(t)=S_{1}(t - 4)-$ $0.38S_{2}(t)-0.2s2(t)s_{4}(t)+(4-1)\cross 0.3S_{\mathit{5}}(t)$ 。
$c$, 恒暗条件での$S_{5^{-}}t$の曲線条件は Fi$g2\mathrm{A}$と同–。

B:光の効果と周期や振幅との関係。 光の効果が、 87 より大きいとき、振幅が、 に$0$なる。 そのとき、 リズムが消滅する。

45. per 突然変異と温度補償性について
は、完全な per-protein が作らにすることによりシミ レー
れないことが分かっているので、 I1=0、あるいは島ションできる。 このときは、安定な解が得られる。 周期が短く変化する突然変異である pers や周期が長く変化する突然変異である Per については、 32 で述べた様々な可能性が考えられ、何が原因であるのか確定しがたい。 また、 この温度補償性についても同様に様々な可能性や調節を考えることができる。すでに行われた実験結果の中に、実験者の意図にかかわらずに、 これらの疑問に答えうる内容が含まれているかもしれない。 また、数理モデルを使うことによって、 どのようなことがわかれば疑問に答えられるかについて詳細に検討することができるだろう。
この数理モデルは、 $per^{0},per,$ per L, timeless の各 mutation は、 何が変化した結果であるのか、 また、 温度補償性はどのように実現されているのかという疑問や実験から生じた数々の疑問にたいして、 実験を実施する上での指針を与えることができ
概日リズムの周期が消滅する突然変異である$per^{0}$$=\mathit{0}$$\epsilon$

るであろう。

$=$

29

5. まとめモデルを表現する微分方程式に関して

1. 微分方程式に delay parts として、 直接表現されている time delay を explicit
time delay. 式の形や係数の値によって生じる極値の時刻の遅れをうみだす原因となる過去の影響を、 implicit time delay と呼ぶことにする。
2. 自励系の常微分方程式が、周期解を発生していれば、広義の time delay を含む
閉ループが存在すると考えられる。
3. ここで考えているような単純な反応を表現する微分方程式では、未知の調節や反応 (black box) が生み出す implicit time delay と explicit time delay の和に相当する値を、 explicit time delay に加えることによって、 同–の周期をもつ式を作ることができる。 また、解曲線の極値の時刻の遅れは、 implicit time delay と explicit

time delay に強くに依存し、 これらの和が大きくなると遅れは大きくなる。
4. 得られた概日リズムの数理モデルは、シミ レーションによって、 24 時間周期と$=$なる定数があることはわかったが、これを数学的に証明することが今後の課題となる。
ショウジョウバエの概日リズムの数理モデルに関して

1. 分子生物学的研究の結果にもとづいて、 DNA の転写調節の negative feedback
loop に explicit time delay を含む概日リズムの数理モデルを作成した。 このモデル
は、 24 時間の周期を発生させることができる.

2. この数理モデルは、

cycle に同調する。 そのとき、恒暗条件より強いリと ZT12 に極大値をもつと予ズムとなっている。 また、単独の PER-protein は想される。
3. この数理モデルは、 実際の PRCs と同様に、 光の強さに応じて 2 種類の
PRCsを描いて、 位相が変化する。
4. この数理モデルは、 強い恒明条件のもとではリズムの消失がおこる。 また、 弱い恒明条件のもとでは、 周期が短くてかつ弱いリズムになると予想される。
5. この数理モデルは、係数の値を変えることによって、周期を変化させることが$\mathrm{L}12\mathrm{D}12$$\mathrm{Z}\mathrm{T}1$できる。 したがって、 周期の突然変異である

timeless をシミレーションすることができる。 また、 温度補償性をシミ . レーションすることもで$per^{0_{\text{、}}}per_{\text{、}^{}s}per^{L}\text{、}$$\iota$

$=$

きる。

6.

すでに発表されている実験結果を基にして、数理モデルから per 突然変異の原
因や温度補償性を求めること、 あるいは何を調べればよいかを示すことが今後の課
題となる。

30

参考文献
1. C. Walter, The Occurrence and the Significance of Limit Cycle
in Controlled Biochemical Systems,J.Theor.Biol.,27, 25&272(1970).
2. C. M. Mackey and L. Glass, Oscillation and Chaos inControl
Systems,SCIENCE, 197,287-289(1977).
3. H. Zeng,P. E. Hardin,M. Rosbash, Constitutiveof the Drosophila period
inhibits period
EMBO Joumal, $13,3590- 3598(1994)$ .
4. T. Kondo et , Circadian Clock Mutants of Cyanobac teria,S CIENCE,266,1233-1236(1994).
5. M. Hunter-Ensor,A. Ousley,A. Sehgal, Regulation of Drosophila Protein Timeless Suggests a
Mechanism for Resetting the Circadian Clock by Light ].
6. H. Zeng,Z. Qian,M. Myers,M. Rosbash, A light-entrainment mechanism for the Drvsophila
circadian clock,NATURE,380, 129135 (1996).
7. M. P. Myers,K. Wager-Smith,A. Rothenfluh-Hilfiker,M. W. Young, Light-Induced Degradation of Timeless and Entrainment of the Drosophila Circadian Clock,SCIENCE,271, 17361740(1996).
8. C. Lee,V. Parikh,T. Itsukaichi,K. Bae,I. Edery, Resetting the Drvsophila Clock by Photic Rqulation of PER and a PER-TIM Complez,SCIEN CE,271,1740-1744(1996).
9. P. E. Hardin,J. C. Hall,M. Rosbash, Feedback of thegene Product on circadian
cycling of its messenger $RNA$ levels,NATURE,343,53&540(1990).
10. B. D. Aronson,K. A. Johnson,J. J. Loros,J. C. Dunlap, Negative Feedback Defining a Circadianof the Clock Gene ffequency,SCIENCE,263, 157&1584(1994).
11. L. B. Vosshall,J. L. Price,A. Sehgal,L. Saez,M. W. Young, Block in Nuclear Localization ofProtein by a Second Clock Mutation, timeless,SCIENCE,263,1606-1609(1994).
12. K. D. Curtin,Z. J. Huang,M. Rosbash,
Regulated Nuclear Entry of the $m_{so_{P^{hila}}}$period Protein Contributes to the Circadian Clock,Neuron,14,365-372(1995).
13. I. Edery,L. J. Zwiebel,M. E. Dembinska,M. Rosbash, Temporal phosphorylation of the
Drosophila period.
14. Z. J. Huang,K. D. Curtin,M. Rosbash, PER
Interactions and Temperature Compensation of a Circadian Clock in Drvsophila,SCIEN,116&1172(1995).
\’ei15. B. Alberts et al., Molecular Biology of the Cell Gaxland Publishing,p401-47,$Beh\beta\dot{m}or$$Phy_{\dot{\Re}}ologiCal$$over_{\mathrm{W}^{r}n}es\dot{\alpha}O$$prot\dot{\alpha}n$$mRNAcycling,\mathrm{T}\mathrm{h}\mathrm{e}$$\mathrm{a}1$$\mathrm{C}\mathrm{e}$$1,84,677-\mathit{6}\mathit{8}\mathit{5}(\mathit{1}\mathit{9}\mathit{9}\mathit{6})$$DrosoPhi\iota_{a}P^{e}\dot{n}od$$ol_{\mathit{0}\mathbb{C}}k.\cdot Autoregu\iota_{a}tion$$p\alpha\dot{r}od$$T\alpha np_{om}lly$$p\mathrm{r}ot\dot{\alpha}n,\mathrm{p}\mathrm{r}\mathrm{o}\mathrm{c}.\mathrm{N}\mathrm{a}\mathrm{t}\mathrm{l}.\mathrm{A}\mathrm{c}\mathrm{a}\mathrm{d}.\mathrm{s}\mathrm{C}\mathrm{i}.\mathrm{U}\mathrm{S}\mathrm{A},\mathit{9}1,226\alpha \mathit{2}\mathit{2}64(1\mathit{9}\mathit{9}4)$$Prot\dot{\alpha}n$$C\mathrm{E},\mathit{2}67$$1\mathit{9}\mathit{9}4,(\mathrm{t}\mathrm{h}\mathrm{i}\mathrm{r}\mathrm{d}$tion).

16. J. E. Darnell Jr, Variety in the level of gene control in eukaryotic cdls,NATURE,2$\mathit{9}7,36\triangleright$

371(1982).

17. A. T. Winfree, Integrated View of Resetting a Circadian Clock,J.Theor.Biol.,28,327-374(1970).
18. J. L. Price,M. E. Dembinska,M. W. Young,M. Rosbash, SuPPression of Period
abundance and circadian cycling by the Drosophila clock mutation timdess,The EMBO
Joumal, 14,4044-4049(1995).
$\mathrm{p}rot_{\dot{\theta}n}$謝辞

数理モデルの数理解析と数値解析について、 大阪府立大学工学部数理工学科の原惟行教授、宮崎倫子博士をはじめとする研究室の人々の暖かい指導と有益な議論に感謝する。また、大阪府教育委員会より、研彦の機会を与えられたことに感謝する。