《Journal of Computational Science》：Boolean Modeling and AI Validation of the ATM/miR-34a–miR-200 Axis Regulating Survival–Apoptosis switch in Breast Cancer
编辑推荐：
共济失调毛细血管扩张突变（ATM）蛋白在阻碍癌症进展中发挥关键作用，传统上因其参与DNA损伤应答（DDR）而闻名。近期研究表明，ATM诱导Snail1和ZEB1的表达，这两者是上皮-间充质转化（EMT）和药物耐药的关键驱动因子。微小RNA-34a（miR-34
共济失调毛细血管扩张突变（ATM）蛋白在阻碍癌症进展中发挥关键作用，传统上因其参与DNA损伤应答（DDR）而闻名。近期研究表明，ATM诱导Snail1和ZEB1的表达，这两者是上皮-间充质转化（EMT）和药物耐药的关键驱动因子。微小RNA-34a（miR-34a）和微小RNA-200（miR-200）均为p53诱导的miRNA，它们直接靶向Snail1和ZEB1，从而抑制EMT和药物耐药。然而，ATM信号通路内的动态调控相互作用仍不清楚。为解决这一问题，研究人员构建了一个基于文献癌症研究的动态布尔模型，并使用来自PubMed的独立实验证据进行了验证。随后，利用AI驱动的文献挖掘平台PubTator3评估模型预测，以验证其与近期实验发现的预测一致性。功能获得和功能丧失扰动分析进一步揭示了ATM和p53诱导的miRNA在促进或抑制DDR和凋亡通路中的双重调控作用。值得注意的是，研究人员发现miR-34a与ZEB1之间存在一个双负反馈环路，该环路作为一个双稳态开关，调控EMT、干性、药物耐药和凋亡。总之，本研究展示了布尔网络建模与基于AI的文献评估相结合如何能够阐明癌症系统中的涌现调控行为。此处揭示的ATM/miR-34a–miR-200轴为乳腺癌中DDR调节和细胞命运动力学提供了机制性和治疗性见解。
研究背景与目的
共济失调毛细血管扩张突变（ATM）蛋白是一种丝氨酸/苏氨酸激酶，长期以来被认为在维持基因组稳定性和抑制肿瘤发生中发挥关键作用。ATM主要通过协调细胞对DNA双链断裂（DSBs）的应答，激活肿瘤抑制因子p53蛋白（p53）信号通路，该通路对于DNA修复、细胞周期检查点激活和凋亡诱导至关重要，从而维护基因组完整性。然而，近年来的证据表明，ATM还能通过直接磷酸化并稳定关键的EMT转录因子Snail1和ZEB1，将DNA损伤信号与上皮-间充质转化（EMT）诱导、干细胞性和药物耐药联系起来，这与其传统的抑癌角色相矛盾。尽管已知p53可诱导miR-34a和miR-200，这些miRNA又能靶向Snail1和ZEB1，但ATM信号通路内这种复杂的动态调控相互作用尚不清晰。因此，本研究的目的是探究ATM在DNA损伤应答（DDR）中是否与p53诱导的miRNA（miR-34a和miR-200）共同作为关键的分子开关，协调促存活（EMT/耐药）与促凋亡状态之间的平衡。
研究方法
研究人员首先利用PubMed和BIOGRID数据库系统梳理已知的调控关系，以此为基础构建了包含15个信号组分（包括两个miRNA）的动态布尔模型。模型构建使用了GINsim软件，并通过MaBoSS软件进行随机模拟以分析网络动态。为了验证模型预测，研究人员采用了AI驱动的文献挖掘工具PubTator3，对模型的预测结果进行独立的交叉验证。具体而言，研究人员设定了功能获得（GoF）和功能丧失（LoF）的扰动条件，模拟了不同基因过表达或敲低对细胞表型的影响，并将模型预测与已有的独立实验观察结果进行比对。
研究结果
3.1 动态模型与野生型场景分析
通过对野生型网络的稳态分析，研究人员发现了三种不同的稳态。第一种对应增殖状态，特征为Mdm2和AKT上调。第二种发生在DDR背景下，表现为ATM、EMT标志物（Snail1、ZEB1）和癌症干细胞标志物SOX2激活，同时凋亡被抑制。第三种则是由p53响应通路（包括miR-34a、miR-200、BAX和Caspase3）触发的凋亡性细胞死亡。MaBoSS模拟显示，在无输入（DDR）时，有50%的概率进入增殖和耐药状态；而当DDR激活后，则有39%的概率转向EMT和干细胞性，另有11%的概率走向凋亡和药物敏感性。
3.2 动态模型准确性评估
通过将模型的功能获得和功能丧失扰动结果与未用于模型构建的独立实验证据进行对比，研究人员验证了模型的准确性。例如，DDR中ATM过表达会诱导Snail1和ZEB1活性，促进EMT和耐药；而p53敲低则会增强EMT和耐药。相反，协同过表达miR-34a和miR-200则会完全转向凋亡和药物敏感性。所有模型预测均与已有实验结果一致。
3.3 ATM在DDR中的双重调控功能
进一步的扰动分析揭示了ATM在不同条件下的双重作用。ATM过表达联合p53敲低会强力激活EMT、干细胞性和耐药（100%）。ATM与p53同时过表达则产生双稳态，约53%的概率偏向促存活表型，47%的概率走向凋亡。ATM敲低联合p53过表达则完全转向凋亡。此外，ATM过表达联合Snail1和ZEB1敲低，或联合miR-34a和miR-200过表达，都能完全逆转至凋亡和药物敏感状态。这表明ATM的作用取决于它与p53、Snail1/ZEB1以及miRNA之间的平衡。
3.4 EMT、干细胞性与凋亡中的反馈环路探索
通过GINsim分析，研究人员在基因调控网络中鉴定出五个功能性回路，其中四个是已知的正反馈回路（miR-34a/Snail1、miR-200/ZEB1、miR-200/Snail1）和一个负反馈回路（p53/Mdm2）。第五个是预测的miR-34a/ZEB1正反馈回路。对该预测回路的扰动分析表明，过表达ZEB1同时敲低miR-34a会完全诱导EMT、干细胞性和耐药；反之，过表达miR-34a同时敲低ZEB1则会完全诱导凋亡和药物敏感性，证实了这个双负反馈环路可能作为一个双稳态开关，调控细胞在耐药与凋亡状态间的切换。
3.5 连接AI与系统生物学：挖掘文献中的实验证据
利用PubTator3进行的AI辅助文献挖掘结果显示，大多数模型预测（如ATM敲低+p53过表达促进凋亡；ATM过表达+p53敲低诱导EMT）在现有研究中得到强烈支持。然而，涉及ATM过表达联合miR-34a/miR-200过表达或联合ZEB1/Snail1敲低的两个预测场景，虽然单个组分的效应已被证实，但ATM在这些组合中的具体作用尚未被直接实验检验，因此被标记为“部分支持”。这指出了未来实验验证的关键方向。
讨论与结论
研究人员的模型揭示了ATM在DDR中充当一个双稳态开关，通过p53–miR-34a/miR-200–ZEB1/SNAIL1轴整合DNA修复控制与EMT可塑性。ATM与p53的平衡决定了细胞是走向促存活的EMT/耐药状态还是促凋亡状态。特别地，新预测的miR-34a/ZEB1双负反馈环路为理解可逆EMT和混合E/M状态提供了新的分子机制，有助于解释肿瘤异质性和治疗抵抗。AI辅助的文献验证增强了模型的可信度，但也指出需要直接的实验来验证ATM与miRNA或EMT转录因子协同作用的精确效果。
研究结论部分翻译如下：
总之，本研究重申了ATM在DDR中的多重作用，证明了其基于与Snail1和/或ZEB1的调控而产生的对比效应。具体来说，ATM与Snail1和/或ZEB1协同作用会促进药物耐药和EMT，同时抑制p53通路。相反，ATM过表达联合Snail1和ZEB1敲低则会强力激活p53、miR-34a和miR-200通路，导致显著的凋亡诱导。此外，一个涉及miR-34a和ZEB1的新型调控回路的鉴定，连同miR-34a和miR-200的关键作用，为理解DDR介导的细胞过程增添了新的维度。这些发现阐明了调控细胞对DNA损伤反应的复杂调控网络，凸显了ATM、miR-34a和miR-200作为癌症治疗中有前景的治疗靶点。
