单一终点的相对活力随检测时间变化
西班牙马德里理工大学(Universidad Politécnica de Madrid)等团队在 Pharmaceutics 研究化疗反应的定量方式。研究以噻唑蓝(thiazolyl blue tetrazolium bromide,MTT)代谢活性检测为输入,比较常规半数抑制浓度(half-maximal inhibitory concentration,IC50)和根据时间序列计算的有效增长率。
模型使用 5 个结直肠癌细胞系和 1 个乳腺癌细胞系:SW480、SW620、DLD1、HCT116、HT29 与 MCF7。论文材料包括奥沙利铂和顺铂,主要图表展示奥沙利铂浓度相关结果。
检测在接种时以及处理后 24、48、72 小时分别进行,每个细胞系包含 3 次独立实验,每次至少 3 个技术重复。检测记录 546 nm 吸光度,常规相对活力以处理组与同时间未处理组的信号比值表示。
不同终点获得的活力曲线和 IC50 明显变化。同一细胞系在相同终点使用不同拟合工具,也产生不同数值;常规软件之间的差别还涉及误差估计方式。
先检查时间序列是否符合指数增长
团队将 MTT 吸光度按接种时信号归一化,再对其对数与时间拟合直线。斜率作为有效增长率:正值表示代谢信号增加,接近零表示净变化很小,负值表示信号随时间下降。
未处理培养物需要先呈现近似指数增长,才能在这一时间窗口使用恒定增长率。HCT116 对照及不同药物浓度下的半对数图支持该拟合,论文同时把空间、营养和高密度等长期约束列为影响因素。
研究以 1000 次自助重采样(bootstrap)估计斜率分布,随后对各药物浓度下的增长率拟合带偏移的指数关系。第二层进行 1000 次蒙特卡洛拟合,把增长率估计的不确定性传递到浓度指标,输出均值和相应区间。
MTT 测量反映代谢活性,本文据此估计培养物的净增长。药物可能同时改变细胞代谢、分裂和存活,这套计算没有直接逐个细胞计数,也没有独立分离出生与死亡过程。

窄屏可横向滑动查看完整图表。
两个指标分别对应停止增长和增长减半
研究提出 ICr0,表示拟合增长率达到零所需的药物浓度;ICrmed 则表示增长率降到未处理组一半所需浓度。这两个定义直接指定增长状态,常规 IC50 指定的是某时点处理组相对对照的活力。
对 HCT116,奥沙利铂浓度升高后,有效增长率从对照约 0.30 天⁻¹逐渐降至负值。表中的 6.25 μg/mL 已接近零,较高浓度下信号进一步减少。
表中 1.6 μg/mL 对应约 0.14 天⁻¹,50 μg/mL 对应约 −0.44 天⁻¹。每个估计同时报告误差,接近零的单个浓度点与用整条曲线计算的零增长浓度分别呈现。
原文表 3 中,HCT116 的 ICr0 约为 8.3 μg/mL,ICrmed 约为 2.9 μg/mL;研究的新方法 IC50 约为 5.5 μg/mL。各数值分别指定零净增长、增长减半和相对活力减半。
当细胞在测试浓度内持续增长时,达到零增长所需浓度可能难以确定。作者因此建议,在这类反应中优先考虑 ICrmed,以减少为估计零增长阈值而向测试范围之外延伸。
同一数据的拟合方式也会改变药效排序
HCT116 的奥沙利铂 IC50,在 GraphPad Prism 中约为 5.37 μg/mL,另一在线计算器给出 8.22 μg/mL,新方法给出约 5.5 μg/mL。SW480 对应结果约为 7.89、13.52 和 5.8 μg/mL。
这些比较使用同一研究的数据,显示曲线形式和估计方式可以影响细胞间的相对敏感性。新方法同时给出参数分布,而在线计算器在本次比较中没有提供相同的误差描述。
HT29 的 ICr0 与 ICrmed 分别约为 37.0 和 12.0 μg/mL,高于 HCT116 的对应值。在这套培养条件和药物测试中,较高浓度阈值对应较难减慢增长的细胞。
时间独立的定义仍依赖模型和输入条件
研究还将 4 个时间点减少为接种时和 72 小时两个点,多数估计接近完整时间序列。HCT116 的 ICr0 均约为 8.3 μg/mL,HT29 则由约 37.0 变为 34.4 μg/mL,区间也随输入变化。
少量时间点可以计算近似值,但完整序列更便于检查是否符合指数增长。本文所称时间独立,是在短期恒定增长率和相应浓度拟合成立时,由增长状态定义指标。
这项研究提供了人源细胞药效分析的具体补充:保留初始信号,检查时间关系,再比较浓度下的净增长。6 个细胞系的化疗实验支持其可行性,更多培养形式和不同作用机制仍需分别检查拟合条件。