1. 引言
由于地球的昼夜交替,各类生物体逐渐建立了相应的生理机能 [1]。这种在外界自然昼夜循环下,生物体为了保持与外界光暗循环相同周期而建立的生物节律,称为近日节律(一般在22小时到28小时之间)。近日节律的研究具有非常重要的生物学意义,在多个领域都有应用。这是复杂系统理论在生物体中的一个具体体现 [2] [3] [4] [5]。哺乳动物的近日节律受到体内位于下丘脑视交叉上核(SCN)的内源性时钟控制 [6] [7]。
SCN由大量神经元组成,分为左右两个核,在结构和功能上可划分为腹外侧(VL)和背内侧(DM) [8]。腹外侧能够直接接受来自视网膜的光信号,背内侧接受来自腹外侧的信号 [9]。这些神经元在没有周期光暗循环情况下能够周期地运动,称为自振荡 [10]。自振荡是SCN能够输出自振荡周期信号保持近日节律的根本原因 [11]。其核心机制是基因转录翻译组成的负反馈环 [12] [13]。SCN的神经元之间存在着耦合,不同的神经元之间自振荡周期不同,实验发现神经元的自振荡周期服从正太分布,在22小时到28小时之间 [14] [15]。神经元之间通过各种神经递质进行耦合,同步后能使整个SCN神经网络达到统一的周期 [16] [17]。若神经元受到影响,SCN未达到统一周期,腹外侧与背内侧周期不同步,称腹外侧与背内侧分裂。
SCN能够根据外界环境以及体内各种因素的变化控制生物的近日节律,维持生物体的正常运转 [18]。SCN与外界光暗循环周期保持同步的作用方式,称为导引 [19]。SCN能够导引的外界光暗循环周期范围称为导引范围。实验 [20] 发现光照强度能够增加SCN对外界光暗循环的导引能力。
由于光照和神经递质对SCN保持近日节律周期有重要影响 [21] [22],Gonze [14] 等人综合考虑光照和神经递质耦合,模拟单个神经元振子的自振荡,建立了描述近日节律的Goodwin模型。通过研究发现SCN的自振荡周期与耦合强度的大小成正比。顾长贵 [23] 在Goodwin模型的基础上,建立周期光暗循环下的近日节律模型,考虑耦合强度(服从正态分)的离散程度以及腹外侧神经元振子数占SCN总神经元振子数的比值对SCN的导引范围的影响,讨论了整个SCN网络的集体行为。发现在周期26小时的外界周期光暗循环下,耦合强度的离散能够显著增加SCN的导引能力,但在外界周期为22小时下,SCN的导引能力几乎不增加。
在顾长贵 [23] 的模型中,光照波形为方波。即在周期光暗循环中,假定有光照时其强度取恒定值,无光照时光照强度为0。但是实际生活中,一天中的光照非恒定值,而是随时间变化的。同时考虑到地理位置和季节的不同,实际光照强度会不同,耦合强度的离散程度以及腹外侧神经元振子数占SCN总神经元振子数的比值对SCN导引都有不同影响,因此有必要对在模拟实际光照条件下SCN的近日节律进行研究。
2. 模型描述
周期光暗循环下的近日节律模型 [21] 如下:
(1)
(2)
(3)
(4)
其中
、
、
、
和
分别是第i个神经元振子的mRNA、蛋白、抑制子、神经递质和标度因子,N为总振子数,
是光照项,n是Hill系数。时钟基因mRNA
产生蛋白
,然后激活转录抑制子
,组成单个神经元振子的负反馈环。神经递质
由mRNA
感生,不同神经元振子通过神经递质
的平均场F进行耦合。
(满足正态分布)表示每个神经元振子对平均场F的耦合强度,其平均值
,标准偏差
代表离散程度。当
时,耦合全同。按照文献 [23] 使标度因子
满足均值为1,标准偏差为0.05的正态分布,可以使每个神经元振子的自振荡周期固定在24小时,能够更好地讨论耦合强度的离散程度 [24] [25] 在模拟实际光照下对SCN导引范围的影响。
由于SCN的背内侧和腹外侧对光照信号的感应程度不同,实验 [26] 发现在同一核内,腹外侧能耦合背内侧,但背内侧对腹外侧的影响非常小。用
代表腹外侧神经元振子数占总神经元振子数的比例,定义标号为1到pN (取整数)的神经元振子属于腹外侧,其余神经元振子属于背内侧。其中腹外侧神经递质平均场 [23] 为
(5)
背内侧神经递质平均场 [23] 为
(6)
因此神经递质的平均场
(7)
相同步序参量可以描述腹外侧和背内侧之间的周期同步程度,其定义如下 [27]:
(8)
代表对时间求平均,其中
和
分别代表从平均场
和
的时间序列中用希尔伯特转换提取出来的相位信息。
时代表腹外侧和背内侧的相位完全不一致,
时代表腹外侧和背内侧的相位完全一致。
周期光暗循环下的近日节律模型中光照波形是方波形式,光照项为 [23]
(9)
考虑到实际光照强度在一天内随时间的变化情况,将光照项改为模拟实际光照,即
(10)
其中t为外界时间,
为外界光暗循环的周期,
为光照强度(
), 代表外界时间t对外界光暗循环的周期
取模,振子总数N = 100,其他模型参数同文献 [25]:
采用四阶Runge-Kutta算法进行计算,步长0.1 h,舍弃前2000步暂态,保留最后的4000步,计算腹外侧和背内侧的平均场时间序列的周期。
3. 结果分析
3.1. 腹外侧振子数占总振子数的比值p对SCN周期的影响
为方便起见,假设耦合强度全同。首先讨论腹外侧振子数占总振子数的比值p对腹外侧和背内侧导引到外界周期的影响。在模拟实际光照下,SCN周期变化情况如图1所示。图1(a)中,
,腹外侧和背内侧周期与外界周期同步程度很高,腹外侧与背内侧周期峰值接近;图1(b)中,
,腹外侧与背内侧的周期都与外界光暗循环的周期一致,腹外侧与背内侧周期峰值相同,且相位同步。因此在相同的光照强度下,增加腹外侧振子数可以增加SCN对外界光暗循环的导引能力,能够帮助SCN与外界光暗循环周期保持一致。在其他条件相同时,在周期光暗循环下的腹外侧振子数占总振子数的比值p对腹外侧和背内侧导引到外界周期的情况如图2所示。从图1(a)与图2(a)可以看出,当腹外侧振子数占总振子数的比值p较小时,模拟实际光照下的腹外侧和背内侧的周期更为同步;从图1(b)与图2(b)可以看出,当腹外侧振子数占总振子数的比值p较大时,模拟实际光照下的腹外侧与背内侧与外界周期保持一致的程度更高。这说明模拟实际光照能够增加SCN对外界的导引能力。因此,在实际光照下腹外侧神经元振子数目对SCN的导引有重要影响。
为了更清楚地了解模拟实际光照下,不同的 值对SCN与外界周期关系的影响,用p-Kf平面内的色相图来描述平均场的周期变化。图3显示在模拟实际光照下,耦合强度全同时,腹外侧和背内侧的周期在p-Kf平面内的变化。图3(a)、图3(b)表示在外界周期22小时时,腹外侧(背内侧)与外界环境达到统一周期(深蓝色)。图3(c)、图3(d)表示在外界周期26小时时,腹外侧(背内侧)与外界环境达到周期一致(红棕色)。在外界周期为22小时和26小时下,在光照强度大于一个阈值后,腹外侧周期与外界周期能够保持同步,不受p值的变化影响;但背内侧若要与外界周期保持一致则同时受Kf和p值的影响;在相同的光照下,背内侧与外界周期保持周期同步所需的p值更大。因此,在模拟实际光照下,外界周期为26小时比周期为22小时的腹外侧与背内侧更易发生分裂。
Figure 1. The time series of the mean field of the ventrolateral part (VL) and dorsomedial part (DVLM) under the simulated actual illumination with a period of 22 hours
图1. 在周期为22小时的模拟实际光照条件下,腹外侧(VL)和背内侧(DM)的平均场的时间序列图
Figure 2. The time series of mean field of the ventrolateral part (VL) and dorsomedial part (DM) under the periodic light-dark cycle with a period of 22 hours [23]
图2. 在周期为22小时的周期光暗循环条件下,腹外侧(VL)和背内侧(DM)的平均场的时间序列图 [23]
在周期光暗循环下,腹外侧和背内侧的周期在p-Kf平面内的变化如图4所示。与图4相比,在模拟实际光照下,代表腹外侧(背内侧)与外界周期一致的深蓝色(红棕色)区域面积均有增加,这表明模拟实际光照能帮助SCN导引到外界周期;在模拟实际光照下比在周期光暗循环下,光照强度的阈值和p值更小,这说明在模拟实际光照下腹外侧能够更快与外界周期达到同步。
3.2. 耦合强度的离散对SCN周期的影响
考虑耦合强度的离散程度对SCN周期的影响,用色相图描述在模拟实际光照下SCN的周期变化情况。图5显示在周期为22小时的模拟际光照下,耦合强度离散时的腹外侧和背内侧周期在p-Kf平面内的变化。在周期为22小时的周期光暗循环下,耦合强度离散时的SCN的周期在p-Kf平面内的变化如图6所示。
Figure 3. The ventrolateral part (VL) and dorsomedial part (DM) periodic hue diagrams in the p-Kf plane under the simulated actual illumination
图3. 模拟实际光照下,腹外侧(VL)和背内侧(DM)周期在p-Kf平面内的色相图
Figure 4. Periodic hue diagrams of ventrolateral part (VL) and dorsomedial part (DM) in the p-Kf plane under the periodic light-dark cycle [23]
图4. 周期光暗循环下,腹外侧(VL)和背内侧(DM)周期在p-Kf平面内的色相图 [23]
Figure 5. Under the simulated actual illumination with a period of 26 hours, the periodic hue diagrams of the ventrolateral part (VL) and dorsomedial part (DM) in the p-Kf plane
图5. 周期22小时模拟实际光照下,腹外侧(VL)和背内侧(DM)周期在p-Kf平面内的色相图
Figure 6. Under the periodic light-dark cycle with a period of 22 hours, the ventrolateral part (VL) and dorsomedial part (DM) periodic hue diagram in the p-Kf plane [23]
图6. 周期22小时周期光暗循环下,腹外侧(VL)和背内侧(DM)周期在p-Kf平面内的色相图 [23]
图5(a)、图5(b)表示在标准偏差
时,腹外侧(背内侧)与外界周期达到周期一致(深蓝色)。图5(c)~(d)表明在标准偏差
时,腹外侧(背内侧)与外界周期达到统一周期(深蓝色)。可以看出,在外界周期为22小时的模拟实际光照下,随着耦合强度离散程度的增加,代表腹外侧(背内侧)与外界周期达到同步的深蓝色区域面积显著增加,说明耦合强度离散程度的增加能帮助SCN更快地导引到外界周期。与图6相比,在外界周期为22小时的模拟实际光照下,代表腹外侧(背内侧)与外界周期达到一致的深蓝色区域面积均略有增加,说明模拟实际光照能够增加SCN导引到外界周期的能力。
在周期为26小时的模拟实际光照下,图7显示在耦合强度离散时的SCN周期在p-Kf平面内的变化。图7(a)、图7(b)表示在标准偏差
时,腹外侧(背内侧)与外界周期达到同步(红棕色)。图7(c)、图7(d)表明在标准偏差
时,腹外侧(背内侧)与外界周期达到统一(红棕色)。可以看出,在外界周期为26小时的模拟实际光照下,随着耦合强度离散程度的增加,代表腹外侧与外界周期一致的红棕色区域面积无明显变化,代表背内侧与外界周期同步的红棕色区域面积显著增加。这说明在周期为26小时的模拟实际光照下,增加耦合强度的离散能够显著增加SCN导引到外界周期的能力。
Figure 7. Under the simulated actual illumination with a period of 26 hours, the ventrolateral part (VL) and dorsomedial part (DM) periodic hue diagrams in the p-Kf plane
图7. 在周期为26小时的模拟实际光照下,腹外侧(VL)和背内侧(DM)的周期在p-Kf平面内的色相图
Figure 8. Underthe periodic light-dark cycle with a period of 26 hours, the periodic hue diagrams of the ventrolateral part (VL) and dorsomedial part (DM) in the p-Kf plane [23]
图8. 周期为26小时的周期光暗循环下,腹外侧(VL)和背内侧(DM)的周期在p-Kf平面内的色相图 [23]
在周期为26小时的模拟实际光照下,图8显示在耦合强度离散时的SCN周期在p-Kf平面内的变化。与图8对比,在模拟实际光照下,代表腹外侧(背内侧)与外界周期一致的红棕色区域面积略增加,说明模拟实际光照能增加SCN导引到外界周期的能力。
3.3. 相同步序参量
相同步序参量能够衡量网络的同步程度。我们进一步在p-Kf平面内讨论在实际光照下耦合强度的离散程度对相同步序参量的影响。
图9显示在模拟实际光照下相同步序参量在p-Kf平面内的变化。可以看出,在外界周期为22小时下,耦合强度的离散程度的增大使序参量 值略有减小,说明耦合强度离散程度的增大能增加腹外侧与背内侧的分裂。在外界周期为26小时的模拟实际光照下,耦合强度的离散程度的增大使序参量R值显著增加,即SCN与外界环境保持周期同步的能力增强。
在周期光暗循环下,相同步序参量在p-Kf平面内的变化如图10所示,腹外侧与背内侧在22小时的周期光暗循环下更易分裂。可以看出,在模拟实际光照下,序参量R值在外界周期为22小时下整体略微增大,腹外侧与背内侧的分裂减少;在外界周期为26小时下,序参量R值整体显著增大,SCN与外界环境更易保持周期同步。即在模拟实际光照的近日节律模型中相同步程度更高,这也说明我们改进的近日节律模型精度更高。
Figure 9. The variation diagram of phase synchronization order parameters in the p-Kf plane under the simulated actual illumination
图9. 模拟实际光照下,相同步序参量在p-Kf平面内的变化图
Figure 10. The variation diagram of phase synchronization order parameters in the p-Kf plane under the periodic light-dark cycle [23]
图10. 周期光暗循环下,相同步序参量在p-Kf平面内的变化图 [23]
4. 结论
本文首先提出更加符合实际的近日节律模型,然后讨论在不同光照强度下腹外侧神经元振子数占总神经元振子数的比值以及神经递质耦合强度的离散程度对生物钟的导引的影响。数值图像结果表明本文构建的近日节律模型比已有模型精度更高。说的更具体一些,和在周期光暗循环下一样,在模拟实际光照下的光照强度同样对SCN的导引起决定性作用,腹外侧所占神经元振子比值与神经递质耦合强度的离散程度对SCN的导引影响也较小;但在模拟实际光照下的背内侧与腹外侧的分裂现象更少,SCN导引到外界周期的能力更强。