《PLOS Computational Biology》：A compartmentalized SCN model explains temperature resistance through interregional phase differences
编辑推荐：
温度周期能够使许多外周组织的昼夜节律发生同步化，然而哺乳动物的主昼夜节律起搏器——视交叉上核（suprachiasmatic nucleus, SCN）却对这些相同的输入信号基本不敏感。尽管中枢和外周时钟具有相似的分子回路，但其对温度刺激产生不同反应的潜在机制
温度周期能够使许多外周组织的昼夜节律发生同步化，然而哺乳动物的主昼夜节律起搏器——视交叉上核（suprachiasmatic nucleus, SCN）却对这些相同的输入信号基本不敏感。尽管中枢和外周时钟具有相似的分子回路，但其对温度刺激产生不同反应的潜在机制仍不清楚。为了解决这一问题，研究人员开发了一个整合了热休克蛋白（heat shock protein, HSP）通路的SCN昼夜节律分区模型，模拟了SCN网络对温度变化的响应。利用该模型，研究人员成功复现了关键的实验观察结果，包括热刺激能够使外周时钟的相位发生倒转，而SCN则保持其抵抗性。进一步的分析揭示，SCN亚区之间的内在相位异质性（phase heterogeneity）是这种同步化抵抗性的基础。具体而言，背侧与腹侧SCN之间的相位差导致两者在受到同时刺激时产生不同的相位反应，从而使得整体SCN相位的净偏移减弱。此外，背侧和腹侧区域被同时刺激的热输入模式可能进一步增强SCN对温度同步化的稳健性。这些结果表明，SCN对温度同步化的抵抗机制可能在机理上不同于光同步化，并强调了其空间结构在维持昼夜节律稳定性中的潜在重要性。
研究背景、现存问题与研究目的
昼夜节律是广泛存在于生物体中的基本计时机制，在哺乳动物中呈现层级组织。主起搏器位于下丘脑前部的视交叉上核（suprachiasmatic nucleus, SCN），而外周时钟存在于几乎所有组织和器官中。外周振荡器可受摄食时间表和温度周期等外界线索调节，SCN则通过神经体液途径协调外周时钟，以维持机体在不同环境条件下的昼夜节律一致性。SCN包含约20000个密集排列的神经元，通常分为腹外侧区和背内侧区（也称为核心区和外壳区），这两个亚区通过有序的细胞间耦合网络相连，支持SCN内一致的昼夜节律振荡。其中，腹侧SCN富含血管活性肠肽（vasoactive intestinal peptide, VIP），而背侧SCN则以精氨酸加压素（arginine vasopressin, AVP）表达为特征。
在细胞层面，昼夜节律由保守的转录翻译反馈环生成，其中BMAL1:CLOCK复合物激活Period（Per1, Per2）和Cryptochrome（Cry1, Cry2）基因的转录，随后积累的PER和CRY蛋白形成复合物并反馈抑制BMAL1:CLOCK活性。温度波动可通过热休克蛋白（heat shock protein, HSP）和热休克因子（heat shock factor 1, HSF1）介导的细胞反应重置昼夜节律。实验研究表明，尽管核心时钟机制和温度响应信号通路高度保守，SCN与外周组织对温度同步化的反应存在显著差异：外周组织如垂体和肺可被倒置的温度周期重置相位，而完整SCN的相位基本保持不变。更关键的是，当SCN的腹外侧和背内侧区域被物理分离后，SCN变得对温度诱导的相移敏感。这表明SCN对温度同步化的抵抗并非单个细胞振荡器的固有属性，而是耦合网络的涌现特性。然而，空间组织、区域间耦合和内部相位结构如何共同产生这种稳健性仍不清楚。为此，研究人员开发了一个分区数学模型，将温度响应信号与SCN的区域结构相结合，旨在阐明内部相位偏移如何塑造组织水平对热输入的响应。
关键技术方法
本研究采用的主要技术方法是数学建模与计算机模拟。研究人员首先建立了一个无空间异质性的温度调控昼夜节律模型，该模型将热休克蛋白信号级联与核心昼夜节律基因网络耦合，作为描述外周时钟的基础框架。在此基础上，他们扩展构建了一个分区SCN模型，将背侧和腹侧区域表示为独立但耦合的两个隔室，每个隔室包含相同的核心时钟网络和局部温度响应通路。模型通过VIP和AVP两种神经肽介导区域间通讯，其中VIP由腹侧产生并对两个区域均有作用，而AVP由背侧产生并仅作用于背侧。研究人员利用常微分方程组描述这些动力学过程，并通过改变初始条件系统性地调节背侧-腹侧相位差，以考察其对SCN整体热响应的影响。此外，他们还引入了一个简化的双振荡器庞加莱模型（Poincaré model）来验证核心机制的普适性，并通过拉丁超立方采样进行参数敏感性分析。
研究结果
模型捕捉到外周时钟和SCN不同的热响应
研究人员首先验证了模型能否复现外周时钟与SCN对温度刺激反应的实验差异。均匀模型（代表外周时钟）对急性热脉冲的响应显示，Per2 mRNA出现急性升高，其诱导幅度落在实验观测范围内。在模拟倒置温度周期（36°C/38.5°C）的实验中，均匀模型的Per2振荡最终演变为近180°的反相位状态，这与在肺和垂体组织中观察到的实验结果一致。相比之下，分区SCN模型在相同倒置温度周期下仅表现出微小的相位调制，未能发展出反相位的Per2节律，这与SCN抵抗温度同步化的实验证据相符。值得注意的是，SCN模型在受到单次热脉冲时仍表现出明显的急性Per2升高，表明其抵抗性并非源于缺乏热响应能力。
区域相位差异是SCN抵抗温度同步化的基础
为了探究内部相位异质性如何塑造集体热响应，研究人员分析了SCN对单次热脉冲的反应。当在CT6施加2小时热脉冲时，整体SCN的Per2仅延迟约1小时，但腹侧SCN显示出约4小时的相位提前，而背侧SCN则延迟约2.4小时。由于整体SCN信号是腹侧和背侧振荡的加权积分，这两个相反方向的反应在组织水平上部分抵消，导致净相移很小。类似现象在CT18施加热脉冲时也被观察到。虽然并非在所有刺激时间都出现相反方向反应，但比较腹侧、背侧和整体SCN的相位响应曲线（phase response curve, PRC）发现，整体SCN通常表现出相对于亚区反应的衰减相移。通过改变腹侧:背侧权重比和进行拉丁超立方参数扰动分析，研究人员确认这一结论具有稳健性，即只要两亚区间存在相位差，整体SCN的相移通常介于相应的腹侧和背侧相移之间。
SCN热抵抗性随背侧-腹侧相位差缩小而降低
研究人员通过解耦两个隔室并调整初始条件来系统性地改变背侧-腹侧相位差。结果显示，当相位差仅为1小时时，热脉冲后整体SCN的相位延迟约为2.6小时，接近各亚区的单独相移。随着相位差增大至4小时，整体SCN的净相移降至约1.9小时；当相位差增至8小时（接近自然状态下的8.6小时），净相移进一步减少至约1.1小时。系统量化表明，整体SCN的相移随区域间相位差的增加而单调递减，并在约8小时处达到最小值。在任何给定的刺激时间，更大的背侧-腹侧相位差都会产生更小的SCN整体相移，证实了区域间相位差是SCN热抵抗性的关键决定因素。
热输入的空间模式可能有助于SCN热抵抗性
研究人员利用简化庞加莱模型和分子分区模型，比较了双重刺激（两个振荡器同时受扰）与仅腹侧刺激对整体相位的影响。结果显示，在大约半个昼夜周期内，双重刺激产生的整体相移小于仅腹侧刺激，尽管前者输入更强。这是因为在这些时间段内，两个振荡器对相同扰动产生了相反方向的相移，从而降低了组合信号的净相移。即使在增强热刺激强度的情况下，双重刺激仍然产生最小的相位漂移，同时伴随较大的急性Per2升高，表明在区域相位异质性存在时，双重刺激可以在某些昼夜相位限制整体相位漂移的同时产生显著的急性反应。
讨论与结论
本研究通过分区SCN模型识别出群体水平的组织是热抵抗性的重要来源，超越了分子层面的不敏感性。背侧和腹侧SCN亚区之间的耦合及其固有的不同相位，通过区域间补偿削弱了净相移。研究还表明输入传递的空间模式也影响稳健性，在某些相位，同时激活两个亚区能在产生较大急性反应的同时诱导较少的净相位重置。这些发现强调SCN架构（包括内部相位结构和刺激模式）是热扰动下功能稳定性的重要贡献者。该机制区别于振幅驱动的易感性和耦合介导的稳定化，其主要源于同一热刺激在两个SCN亚区诱发不同的相移，进而在整体SCN相位中部分抵消。然而，模型预测的残余相移（超过1小时）仍大于实验中观察到的极小漂移（约0.3小时），提示体内可能存在更强的区域内耦合等额外稳定化过程。研究结论指出，SCN对温度同步化的抵抗机制可能在机理上不同于光同步化，并强调了其空间结构在维持昼夜节律稳定性中的潜在重要性。
