STAT1–p53协同调控模型：通过调节p21与PUMA编程应激条件下的细胞命运

《Journal of Theoretical Biology》：A Model of STAT1–p53 Coordination: Regulating p21 and PUMA to Program Cell Fate under Stress

　　

信号通路在应激诱导的细胞命运决定中发挥精密调控作用，其中p53的稳定性通过p53–Mdm2负反馈的衰减实现，是有效响应的核心机制。信号转导与转录激活因子1（STAT1）既能抑制Mdm2转录，又能共激活p53靶基因（p21和PUMA），但这两重功能尚未在统一框架下被定量描述。研究人员构建了一个机制模型，将生物学校准的应激尺度S与STAT1活性Σ耦合至经典的p53–Mdm2核心模块及下游p21与PUMA模块。关键过程不再由固定参数描述，而是显式依赖于S的函数，使模型能够随应激水平自适应调整，无需人工重新拟合参数。这一形式确保振荡、阻尼或平台型响应可仅通过改变S自然出现。模型经校准再现了DNA双链断裂（DSB）类损伤（如电离辐射及放射模拟剂、拓扑异构酶抑制剂）下典型的p53脉动/振荡模式。STAT1作为转录共激活因子调控启动子强度。分岔分析揭示了（S, Σ）平面上的二维霍普夫区域；Σ升高将该区域推向更高S值并逐渐收窄，最终在高Σ下抑制振荡。在具有时变S与Σ的通用细胞模拟中，结合两阶段决策规则（早期瞬态检测+稳态读取），可将p53动态映射至不同命运：持续的中等振荡对应细胞周期阻滞，中等阻尼响应对应衰老，强阻尼的高平台对应凋亡。该模型在不重新拟合动力学参数的情况下复现了MCF-7、HCT116、U2OS等不同细胞系的命运分布，并凸显细胞类型对p21与PUMA的敏感性差异。此框架将STAT1定义为应激响应p53网络的可调放大器与振荡抑制因子，提供了基于动力学的可检验命运调控预测。

研究背景方面，细胞在面临多种应激源（如紫外线、γ射线、病毒感染）时会启动复杂的信号网络，其中肿瘤抑制蛋白p53处于核心地位。基础状态下，p53受其主要调控因子Mdm2的负反馈抑制，维持在低水平；一旦DNA受损，p53稳定性提高并诱导Mdm2表达，形成振荡性调控特征。已有研究表明，STAT1除参与细胞因子信号外，还能稳定p53并增强其靶基因表达，但缺乏将其整合入p53–Mdm2动力学模型的定量研究。这种缺失限制了对不同应激与STAT1活性组合如何共同塑造下游效应蛋白p21与PUMA的动态及最终细胞命运的理解。

研究人员在比利时布鲁塞尔自由大学的理论时间生物学实验室构建了包含应激尺度S与STAT1活性Σ输入的机制模型，将STAT1的双重作用（抑制Mdm2转录与共激活p21/PUMA）纳入经典p53–Mdm2模块。技术方法上，模型采用应力依赖的参数函数替代固定参数，实现跨应激水平的自适应；引入两阶段命运判定规则（早期瞬态阈值检测与60–72小时稳态读取），并在（S, Σ）平面上进行分岔分析，识别振荡区域及其变化规律；此外，在单一通用细胞及多个细胞系（MCF-7、HCT116、U2OS）中验证模型，无需重新拟合生化参数。

研究结果分为以下部分：

1. 1.
   模型构建：将应激S与STAT1活性Σ作为外部输入，纳入p53–Mdm2核心模块，并通过Hill函数等形式将S嵌入关键过程，使模型能自动产生振荡、阻尼或平台响应。
2. 2.
   稳定性与分岔：分岔分析揭示二维霍普夫区域，Σ升高导致振荡区域右移并缩小，高Σ下振荡消失。
3. 3.
   细胞命运判定：基于p21与PUMA的阈值及停留时间，将p53动态映射为阻滞、衰老或凋亡三类命运。
4. 4.
   从单细胞到原型细胞系：同一判定规则在不同细胞系中复现文献报道的命运分布，揭示p21与PUMA敏感性的细胞类型差异。
5. 5.
   讨论与结论：该应力感知模型成功整合STAT1的调控轴，再现DNA损伤响应的典型动态序列；在无需参数重调的情况下预测不同S与Σ组合下的命运走向，并提出可检验的实验预测。

讨论部分指出，本研究弥补了p53–Mdm2模型未纳入STAT1调控的空白，并展示了应激与细胞因子信号联合调控的细胞命运决策机制。结论强调，STAT1既是p53网络的放大器，也是振荡抑制因子，其作用取决于应激强度与STAT1活性的组合，这为理解应激响应及开发靶向干预策略提供了理论基础。论文发表于《Journal of Theoretical Biology》。