香豆素–嘧啶结构先前分子对接排名与MDA-MB-231/MCF-10A细胞活性数据的回顾性基准评估

《Biology》：Retrospective Benchmarking of Previously Modeled Coumarin–Pyrimidine Structures Against MDA-MB-231/MCF-10A Cellular Activity

　　

本研究回顾性评估了先前已发表的香豆素–嘧啶化合物库的靶点层面分子建模排名与14个具有明确结构归属的实验评估化合物的现有细胞数据集之间的一致性。将实验室报道的多浓度MTT检测所得的MDA-MB-231和MCF-10A IC50值与已发表的AutoDock Vina评分进行基于秩次的分析比较。化合物5e在实验室报道数据中具有最低的MDA-MB-231 IC50（12.4 ± 0.2 μM；n = 2）以及最大的描述性选择性指数（7.75）；该比值仅针对MDA-MB-231/MCF-10A这一细胞模型比较，不解释为普遍肿瘤选择性的证据。对接排名与细胞效力之间仅呈现弱相关性，且经Benjamini–Hochberg校正后，没有任何对接–表型关联保持统计学显著。任何基础二维描述符均未显示具有说服力的校正后关联。完整的1H/13C NMR及经验证的HRMS文档、定量纯度测量数据以及完整全系列MTT原始导出数据均不可获取；因此，对结构依赖性和效力相关结论均采取审慎解读。总体而言，该分析表明在MDA-MB-231/MCF-10A数据集中，靶点层面对接排名与全细胞活性之间的跨终点可转移性有限，且未建立经过验证的构效关系、泛TNBC活性、凋亡诱导、BCL-2机制或转化疗效。

论文解读

一、研究背景与问题

乳腺癌仍是全球癌症发病和死亡的主要原因之一，涵盖生物学上截然不同的亚型。三阴性乳腺癌（triple-negative breast cancer, TNBC）因缺乏雌激素受体、孕激素受体及HER2过表达/扩增，缺乏内分泌治疗和HER2靶向治疗的选择，且常表现出侵袭性强、早期复发和治疗耐药等特征，临床上极具挑战性。这些特征持续推动针对生存、增殖、DNA修复和凋亡等互补通路的小分子药物开发。

杂合药物化学策略可将不同靶点识别模式相关的药效团元素结合起来。香豆素衍生物作为抗癌骨架被广泛探索，嘧啶基序常见于激酶导向及其他肿瘤学小分子药物中。两者组合为探索多靶点药理学以及结构变化如何同时影响靶点层面分子识别和表型细胞反应提供了化学上可行的平台。

研究人员先前的工作中，一个包含64个成员的香豆素–嘧啶化合物库（含羟基和巯基嘧啶系列）仅通过分子对接、200 ns分子动力学模拟、MM-GBSA和MM-PBSA分析、药效团建模及计算机ADMET分析，对PARP-1、EGFR、HER2和BCL-2四个靶点进行了纯计算机评估，并明确提出了后续合成和实验评估的建议。本研究中实验研究的化合物属于同一先前建模的化学空间，因此先前的计算结果属于先前预测而非本工作的新结果。

现有细胞数据的出现创造了一个不同于简单重复原始建模流程的机会：可将计算排名与观察到的细胞活性进行回顾性检验。此类比较具有信息价值，因为对接分数近似于简化模型下靶点–配体相互作用的有利程度，而全细胞细胞毒性整合了通透性、细胞内暴露、靶点表达、代谢稳定性、多靶点药理学及下游细胞反应。因此，靶点特异性分数与表型效力之间的不完全一致具有科学信息价值，而非本质上的矛盾。

二、研究内容与结论概述

本研究对14个化合物（5a–5n）在MDA-MB-231细胞和非致瘤性MCF-10A细胞中的可用实验随访数据进行了回顾性检验。其具体贡献不在于“对接分数与细胞毒性代表不同终点”这一一般性观察，而在于对先前明确前瞻性发表的多靶点建模研究，针对相应化学系列后续生成的实验数据集进行化合物层面的回顾性审计。这种结构匹配的比较识别了个体计算优先级在细胞测试后被保留或取代的位置，为使用实验观察之前发表的预测进行模型评估提供了透明的案例研究。

研究的主要结论是：在可用数据集中，化合物5e具有最低的MDA-MB-231 IC₅₀（12.4 ± 0.2 μM）和最大的描述性选择性指数（7.75）；先前发表的对接排名未能直接重现全细胞效力排名；所有对接–表型关联经多重性校正后均不显著；二维描述符分析也未发现具有说服力的关联。由于化学表征文档和生物学重复不完整，研究聚焦于描述性跨终点基准评估，不声称结构新颖性、经过验证的构效关系、跨TNBC疾病类别的活性或已确认的分子机制。

三、主要技术方法

研究人员采用的主要技术方法包括：回顾性基准评估设计（将先前发表的AutoDock Vina 1.1.2对接评分作为先验预测，与实验室报道的MDA-MB-231和MCF-10A IC₅₀值进行Spearman秩相关和Kendall τb分析）；Benjamini–Hochberg错误发现率多重性校正；留一法敏感性分析；基于RDKit 2025.09.4的二维理化描述符分析（包括分子量、Crippen cLogP、拓扑极性表面积、氢键供受体计数、可旋转键数、sp3碳分数和芳香环数）；以及描述性选择性指数计算。细胞来源方面，MDA-MB-231人TNBC细胞获自NCCS（印度浦那），MCF-10A非致瘤性乳腺上皮细胞获自ATCC。统计分析在Python 3.13.5中执行，使用SciPy 1.17.0和statsmodels 0.14.6。

四、研究结果

评估化合物系列及可用分析记录。 实验评估系列包含14个化合物，归属为羟基嘧啶–香豆素结构，差异在于芳基或杂芳基取代基。终产物为黄至橙色结晶固体，分离产率62–79%，熔点范围225–272 °C。13个文档化归属可明确匹配到先前建模库中的条目。现有正离子ESI-MS PDF和部分¹³C记录作为支持性分析记录保留，但不构成完整的结构表征包或定量纯度评估。

实验室报道的细胞毒性和描述性选择性。 外包实验室确认所有化合物5a–5n均经多浓度MTT测量评估。在可用数据集中，化合物5e具有最低的MDA-MB-231 IC₅₀（12.4 ± 0.2 μM），其次为5m（13.8 ± 0.3 μM）和5h（14.8 ± 0.2 μM）。大多数其他化合物IC₅₀值在约20–25 μM范围。5e在MCF-10A细胞中报道IC₅₀为96.1 ± 2.2 μM，对应描述性SI约7.75。5m和5h的SI分别约为6.77和5.72。该排名保持描述性，因每个IC₅₀仅基于两次独立实验，且完整化学表征和定量纯度数据不可用。

先前分子建模排名的跨终点基准评估。 主要回顾性基准包含13个可明确对应到先前库的化合物。PARP-1（Spearman ρ = 0.017，p = 0.956）、EGFR（ρ = 0.042，p = 0.893）和HER2（ρ = 0.081，p = 0.793）的靶点特异性对接评分与MDA-MB-231效力几乎无秩一致性。BCL-2对接与IC₅₀呈中等但非确证性关联（ρ = 0.492，p = 0.087）。四靶点平均秩与实验效力秩弱对应（ρ = 0.204，p = 0.505）。化合物5e实验效力排名第一但平均对接秩为4.25；5m实验排名第二且平均对接秩最优（3.50）；5h实验排名第三但平均对接秩较差（7.00）。九项主要探索性对接–表型比较经Benjamini–Hochberg校正后均不显著。BCL-2对接与SI的关联产生最小校正q值（ρ = ?0.680，p = 0.011，q = 0.095），未达预设q < 0.05阈值。留一法分析中该关联方向持续保留（ρ范围 ?0.769至?0.605），但后验检查不构成确证性证据。

二维理化描述符背景分析。 对全部14个结构计算了八种基础二维描述符。经16项描述符–表型测试的Benjamini–Hochberg校正后，没有任何描述符与MDA-MB-231 IC₅₀或SI显示具有说服力的关联。最大名义趋势为MW与SI之间的负相关（ρ = ?0.506，p = 0.065），仍不显著。这些探索性描述符分析不用于赋予5e排名以结构原因或推断构效关系。

五、讨论总结与结论翻译

讨论部分总结。 本研究最宜解读为先前发表的计算机发现项目与后续实验室报道细胞数据之间的回顾性比较。先前工作为64个香豆素–嘧啶结构建立了跨四个乳腺癌相关靶点的分子建模景观。回顾性基准明确区分了靶点层面建模与全细胞表型：PARP-1、EGFR和HER2的先前对接排名与报道的MDA-MB-231效力几乎无对应关系，BCL-2仅显示中等非确证性关系。这一结果与“对接排名直接预测细胞毒效力”的简单主张不一致，但这并非意外——对接分数近似特定蛋白口袋中的分子识别，而MTT细胞毒性整合了通透性、溶解度、细胞内暴露、代谢、靶点丰度、脱靶药理学和下游网络反应。该研究的具体贡献在于其时间性和结构层面的设计：靶点层面预测在细胞观察可用之前已发表，分析将先验预测映射到后续测试的相应化学系列，实现化合物层面的优先级保留或取代审计。研究同时指出，BCL-2对接与SI的探索性关联虽方向稳健但未达多重性校正显著性阈值，不赋予机制优先级。研究还讨论了TNBC治疗格局的演变、乳腺癌异质性对单一细胞系比较的局限、MTT终点的方法学限制、化学表征不完整的约束，以及未来研究需补充完整结构表征、定量纯度、正交活力验证和额外TNBC模型等证据。

结论翻译。 本工作对先前发表的靶点层面分子建模排名与相应化学系列的后续MDA-MB-231/MCF-10A细胞数据集进行了化合物层面的回顾性审计。在实验室报道数据中，化合物5e具有最低的报道MDA-MB-231 IC₅₀和最大的描述性SI，该SI仅针对此细胞系比较。先前对接排名未直接转移至细胞效力排名，且经多重性校正后无任何对接–表型关联保持统计学显著。由于完整结构表征、定量纯度数据、完整全系列MTT原始数据集、正交活力验证、充分生物学重复及额外TNBC模型均不可用，本研究不建立身份依赖性构效关系、普遍肿瘤选择性、跨TNBC疾病类别的活性、凋亡诱导、BCL-2机制或转化疗效。主要推论是方法学性的：前瞻性靶点层面分子建模排名与全细胞表型排名是不同终点，需要直接实验交叉验证。次级Annexin V/PI和RT-qPCR记录仅保留在补充信息中用于透明性，不构成本结论。化合物5e因此仅被描述为可用数据集中最活跃的成员，并作为未来重新合成、完整分析表征和独立前瞻性验证的优先对象。