《Frontiers in Bioinformatics》：CUL3-mediated NRF2 dysregulation in TP53-driven non-small cell lung cancer: neddylation-associated structural perturbations
编辑推荐：
CUL3编码Cullin-3，是Cullin-RING连接酶3（CRL3）的支架，通过KEAP1依赖的泛素化介导NRF2的蛋白酶体降解。尽管KEAP1或NFE2L2突变等经典NRF2改变在非小细胞肺癌（NSCLC）中已被充分表征，CUL3突变的作用仍定义不足。
CUL3编码Cullin-3，是Cullin-RING连接酶3（CRL3）的支架，通过KEAP1依赖的泛素化介导NRF2的蛋白酶体降解。尽管KEAP1或NFE2L2突变等经典NRF2改变在非小细胞肺癌（NSCLC）中已被充分表征，CUL3突变的作用仍定义不足。研究人员假设CUL3错义变异在TP53改变的NSCLC中呈条件性富集，EGFR癌基因亚群作为明确的疾病背景而非协同需求。利用两层Fisher精确检验和含EGFR×TP53交互项的逻辑回归，对多队列NSCLC数据集（n = 47,862；其中10,228例具致癌EGFR突变并确认CUL3覆盖）进行分析，并在TCGA LUAD/LUSC中进行转录组学验证。CUL3突变在TP53突变EGFR背景中显示条件性富集（OR = 2.92，p = 4.68 × 10?29）；EGFR×TP53交互不显著（OR = 0.86，p = 0.547）证实富集由TP53缺失而非EGFR协同驱动。CUL3突变肿瘤在三个独立特征中显示NRF2激活升高，伴增殖信号增加和药物代谢改变。对失活和NEDD8化活性状态下野生型、L321F、G402E和M711I CUL3:KEAP1复合物的全原子分子动力学模拟揭示变异特异性构象行为。L321F使NEDD8化复合物进入构象受限状态，伴骨架柔性降低、回转半径减小和自由能最低点偏移；G402E显示渐进性失稳，M711I与野生型相当。自由能景观分析证实L321F + NEDD8相对野生型构象能量最低点偏移，表明尽管复合物形成完整但催化柔性受损。这些发现将CUL3定义为EGFR驱动NSCLC中TP53相关的NRF2调节修饰因子，通过NEDD8化活性CRL3复合物的变构刚性化发挥作用。
研究背景方面，NSCLC是基因组分层疾病，EGFR激活突变定义主要亚群，但EGFR驱动肿瘤并非生物学均一，TP53等共变异常见并与基因组不稳定、应激反应扭曲和临床结局差异相关。CUL3是Cullin-RING E3泛素连接酶（CRL3）核心支架，参与蛋白周转、细胞周期、蛋白稳态和应激适应，并通过与KEAP1形成复合物介导NRF2泛素化降解；KEAP1与NFE2L2改变在NSCLC中已被广泛描述，而CUL3突变低频且系统意义不清。现有研究多聚焦经典NRF2轴改变，尚不清楚CUL3突变是否随机分布、是否独立于经典NRF2轴，以及其结构后果；NEDD8化活性CRL3构象瞬态且难以常规实验解析。研究人员因此整合多队列基因组、转录组和全原子分子动力学，检验CUL3突变在TP53改变背景中的条件性富集及其对NRF2调节和CRL3活性构象的影响。研究结论将CUL3定义为EGFR驱动NSCLC中TP53相关的NRF2调节情境依赖修饰因子，并提出其通过NEDD8化活性CRL3复合物变构刚性化发挥作用的机制模型，为CUL3突变、TP53改变腺癌的NRF2通路或neddylation导向干预及风险分层提供候选方向。
关键技术方法方面，研究人员整合TCGA、MSK-IMPACT、AACR GENIE v19.0与CPTAC等NSCLC队列，患者层面去重后共47,862例，主要分析采用AACR GENIE v19.0，并以v20.0作大队列验证；以Fisher精确检验和含EGFR×TP53交互项的逻辑回归评估条件性富集。转录组分析使用TCGA LUAD/LUSC，计算多个NRF2特征评分、单样本基因集富集分析（ssGSEA）免疫特征、药物代谢基因及neddylation通路共表达。结构分析以PDB 5NLB为模板，结合AlphaFold2建模CUL3:KEAP1复合物，在失活与NEDD8化状态开展全原子分子动力学模拟（GROMACS，CHARMM36m，200 ns），并进行主成分分析（PCA）与自由能景观（FEL）分析。
研究结果中，3.1 多队列EGFR驱动NSCLC中CUL3突变景观：通过五个独立NSCLC测序库汇总，47,862例中10,228例为EGFR致癌突变且具CUL3覆盖，100例CUL3突变，三阳性（EGFR + TP53 + CUL3）排除混杂后为61例；各队列CUL3突变率0.30%—1.72%，合并0.98%，提示CUL3突变在EGFR驱动NSCLC中稳定低频而非明显技术伪影。3.2 CUL3突变的条件性富集：两层Fisher精确检验显示全NSCLC队列中TP53与CUL3共现OR = 2.81，EGFR致癌亚组中OR = 2.51，均显著，表明CUL3突变与TP53突变非随机共现。3.3 EGFR与TP53的非加性交互：逻辑回归显示TP53主效应显著（OR = 2.92，p = 4.68 × 10
?29），EGFR主效应未达显著，EGFR×TP53交互不显著（OR = 0.86，p = 0.547）；驱动分层及v20.0验证显示驱动阴性、TP53定义肿瘤中富集最强，LUAD中关联保留，支持富集由TP53改变驱动而非EGFR协同。3.4 NRF2轴混杂评估：三阳性组中KEAP1、STK11和NFE2L2共突变率分别为22.1%、20.6%和10.3%，六种序贯排除场景后CUL3富集模式仍保留，提示观察信号并非仅由经典NRF2轴病变驱动。3.5 三阳性肿瘤TMB升高：可得TMB的三阳性样本中位TMB为32.9 mut/Mb，高于EGFR + TP53和EGFR-only组，但仅9例且存在 ascertainment 偏倚，研究者将其作为描述性、假设生成观察。3.6 转移富集作为情境严重度指标：活检来源分析未显示三阳性组转移比例显著高于对照，三阳性组介于EGFR-only与EGFR + TP53组之间；组织学分层显示LUAD中CUL3–TP53关联仍显著。3.7 CUL3突变肿瘤中NRF2通路激活：三个独立NRF2特征评分在CUL3突变肿瘤中均显著升高，支持CUL3改变与NRF2转录程序激活相关。3.8 增殖与免疫微环境变化：ssGSEA显示CUL3突变肿瘤增殖特征显著升高，干扰素反应呈方向性降低并接近显著，其他免疫细胞特征方向性降低但多未达FDR显著。3.9 CUL3突变肿瘤中药物代谢基因表达改变：UGT1A1和UGT1A9显著升高，CYP3A5和UGT2B7显著降低，提示向结合型II相清除偏移的转录重编程，但仅为转录关联。3.10 Neddylation通路与NRF2活性耦合：NAE1、UBA3、UBE2F、RBX1、KEAP1和NFE2L2与NRF2激活评分正相关，UBA3在CUL3突变组升高，支持NRF2高输出肿瘤更多参与neddylation机制并支持对NEDD8化活性CRL3建模。3.11 CUL3突变对CUL3–KEAP1支架的结构影响：L321F位于埋藏螺旋区，G402E位于溶剂暴露甘氨酸环，M711I邻近K712 NEDD8化位点；SIFT、PolyPhen-2和ESM-2评分提示多个变异具功能影响，结构推断CUL3变异可能通过变构调节而非直接破坏界面发挥作用。3.12 NEDD8化CUL3:KEAP1复合物构建与轨迹分析：构建K712–G76异肽键并完成200 ns模拟；L321F + NEDD8显示CUL3与KEAP1骨架RMSD、RMSF、Rg降低及界面氢键和CUL3:NEDD8接触增加，G402E + NEDD8呈渐进性骨架失稳，M711I + NEDD8接近野生型，提示L321F导致NEDD8化复合物刚性化而G402E导致失稳。3.13 构象景观分析：PCA与FEL显示野生型、G402E和M711I在NEDD8化后采样空间扩展，而L321F + NEDD8局限为单一窄盆地且能量最低点相对野生型偏移，提示L321F限制构象采样并损害催化柔性。
讨论与结论部分，原文将基因组、转录组和结构分析汇聚为CUL3情境依赖修饰因子模型：CUL3错义变异以L321F为例，通过NEDD8化活性CRL3复合物变构刚性化影响NRF2调节；基因组富集由TP53改变而非EGFR–TP53协同驱动，并对NRF2轴混杂排除稳健；转录组支持NRF2通路激活、增殖信号、免疫微环境重塑和药物代谢改变；分子动力学提示L321F使NEDD8化复合物构象受限、能量最低点偏移和采样幅度降低。由于各系统为单条轨迹，变异特异性构象差异系跨四种CUL3背景比较解释，正式统计置信仍需独立重复。Neddylation共表达与结构发现支持CUL3突变为该亚组NRF2轴失调的非经典但功能相关途径。作为整合计算研究，这些结果作为待实验检验的假设提出，需在CUL3变异模型中通过NRF2报告、CRL3:KEAP1共免疫沉淀和体外neddylation动力学试验验证；将CUL3变异识别为TP53相关、NRF2激活事件，提名CUL3突变、TP53改变腺癌作为NRF2通路或neddylation导向干预及NRF2风险分层的候选情境。研究结论翻译如下：本研究整合多队列基因组分析、转录组分析和全原子分子动力学模拟，将CUL3表征为EGFR驱动非小细胞肺癌中NRF2通路调节的情境依赖修饰因子。三个分析阶段的汇聚发现支持一机制模型：以L321F为代表的CUL3错义变异通过NEDD8化活性CRL3复合物的变构刚性化发挥作用。基因组分析确立CUL3突变富集由TP53改变而非EGFR–TP53协同驱动，EGFR×TP53交互不显著，并在EGFR致癌疾病背景中观察，对连续NRF2轴混杂排除稳健，在六种敏感性场景中保留。转录组分析确认NRF2通路激活升高作为功能后果，并伴增殖信号、免疫微环境重塑和药物代谢基因表达改变，符合NRF2驱动重编程。分子动力学模拟揭示L321F替换使NEDD8化复合物进入构象受限状态，伴能量最低点偏移和构象采样幅度降低，表明尽管复合物完整性保留但催化柔性受损。由于每个系统以单条轨迹模拟，变异特异性构象差异在四种CUL3背景间作比较解释，正式统计置信需独立重复。Neddylation通路共表达数据与结构发现支持一模型：CUL3突变代表该亚组中NRF2轴失调的非经典但功能相关途径。作为整合计算研究，这些基因组、转录组和结构结果作为待实验检验的假设提出，需在CUL3变异模型中通过NRF2报告、CRL3:KEAP1共免疫沉淀和体外neddylation动力学试验验证。将CUL3变异识别为TP53相关、NRF2激活事件，将提名CUL3突变、TP53改变的腺癌作为NRF2通路或neddylation导向干预及NRF2风险分层的候选情境，为未来转化研究方向。
