惯例跑了一组 mRNA-LNP 冻融实验-20/-80/-196°C 三温度 × 0/1/3/5 次循环DLS 复测粒径和 PDI想看包保护剂蔗糖/海藻糖到底护住了没有。这份笔记用 Pandas 做完整的数据整理scipy 从 DSC 热流曲线提取 Tg’/Tm 相变点最后用 matplotlib 画出粒径漂移 vs 冻融循环的稳定性曲线。依赖pandas、numpy、scipy、matplotlib。一、为什么 LNP 怕冻mRNA-LNP 在 2–8°C 缓慢降解水解、氧化长期储存依赖冷冻。但冷冻本身是个破坏过程冰晶挤压冻结时水结冰排出溶质LNP 被浓缩进非冻结相颗粒间距离从微米级压缩到纳米级聚集乃至融合机械应力冰晶生长的物理挤压直接破坏颗粒结构pH 迁移磷酸盐缓冲液冻结时磷酸氢二钠先结晶残余液相 pH 可漂移 3 个单位枸橼酸和组氨酸缓冲液冻融 pH 稳定得多——这就是冻干制剂几乎全用组氨酸的原因相变脂质核心在低温下从液晶态向凝胶态转变Tm脂肪链排直、晶格收缩把内水相的 RNA 挤出来。所以稳定性数据要同时看两件事热学上有没有危险的相变点DSC 的 Tg’ 和 Tm颗粒上冻融几次后粒径分布漂了多少DLS 复测的粒径漂移曲线。二、数据结构一张长表实验设计 3×4 因子温度 × 循环数 冻干保护剂 2 水平5%/10% 蔗糖 阳性对照无保护剂。DLS 每条件 3 独立瓶 × 3 读数。存成一张 CSVtemp_c,cycles,sucrose_pct,rep,reading,d_nm,pdi -20,0,5,1,1,83.2,0.112 -20,0,5,1,2,84.1,0.115 ...读入和描述统计importpandasaspdimportnumpyasnp dfpd.read_csv(freetz_data.csv)print(df.groupby([temp_c,cycles,sucrose_pct]).agg(d_mean(d_nm,mean),d_sd(d_nm,std),pdi_mean(pdi,mean),n(d_nm,size)))先做正态性检查再谈平均冻融后的粒径分布经常出现主体未变 少量 1–3 µm 大颗粒的双峰。此时 Z-average 被大颗粒拖高强度加权对大颗粒极其敏感d⁶ 加权均值既不代表主体也不代表尾巴。两个习惯报 Z-average 的同时报 PDI 和数量分布峰值对漂移率做 log 变换后再比较对数空间下的漂移近似线性。df[drift_pct]df.groupby([temp_c,sucrose_pct,rep])[d_nm]\.transform(lambdas:(s/s.iloc[0]-1)*100)注意这里按 rep独立瓶分组后取第 0 次循环读数为基线不要用全表 cycles0 的均值当基线——瓶间差异会污染循环效应。三、DSC 热流曲线提取 Tg’/Tm差示扫描量热DSC以固定速率降温/升温热流曲线上的吸热/放热台阶对应相变。Tg’最大冻结浓缩液的玻璃化转变是冻结相从橡胶态进入玻璃态的温度——冻到 Tg’ 以下分子扩散近乎停止这就是安全储存温度的下界。Tm 是脂质凝胶-液晶相变。3.1 一步法一阶导数峰值对升温曲线求一阶导数Tg’ 表现为导数曲线上的台阶Tm 表现为峰。用 scipy 的 savgol_filter 平滑后找峰importnumpyasnpfromscipy.signalimportsavgol_filter,find_peaks# dsc.txt: 两列 温度(°C) 热流(mW)升温段 10 K/minT,hfnp.loadtxt(dsc.txt,delimiter,,unpackTrue)# 平滑 一阶导dhfsavgol_filter(hf,window_length51,polyorder3,deriv1,delta1.0)# Tg: 导数台阶 → 对导数再求导二阶导极值d2savgol_filter(hf,window_length101,polyorder3,deriv2,delta1.0)tg_candidates,_find_peaks(np.abs(d2),heightnp.percentile(np.abs(d2),95))print(Tg 候选:,T[tg_candidates])# Tm: 导数峰tm_candidates,_find_peaks(dhf,height0)print(Tm 候选:,T[tm_candidates])一步法快但 Tg’ 的台阶在原始曲线很微弱尤其冻干样品含糖量高时实际方法学里更稳的做法是 DSC 常用的两遍扫描第一遍升温消除热历史第二遍才是分析曲线。3.2 更稳的两点法半高法对 Tm 峰取峰两侧斜率最大点作切线半高交点定熔点。工程上够用的近似峰起点温度 0.5×(峰结束-起点)。deftm_half_height(T,hf,i_peak,window40):峰位 i_peak 处的半高温度近似lomax(0,i_peak-window);himin(len(T),i_peakwindow)basenp.median(hf[lo:hi])halfbase0.5*(hf[i_peak]-base)# 升温侧最后一个低于半高的位置i_halflonp.where(hf[lo:i_peak]half)[0][-1]returnT[i_half]这段近似代码按需取用严格切线法建议直接用仪器软件本文重点在数据流不在切线几何。四、冻融循环粒径漂移曲线把每个条件下、每次冻融后相对于 cycle 0 的粒径漂移串成曲线importmatplotlib.pyplotasplt g(df.groupby([temp_c,cycles,sucrose_pct])[drift_pct].agg([mean,std]).reset_index())fig,axplt.subplots(figsize(7,4.5))colors{-20:#1f77b4,-80:#ff7f0e,-196:#2ca02c}fortemp,subing.groupby(temp_c):s5sub[sub.sucrose_pct5].sort_values(cycles)s10sub[sub.sucrose_pct10].sort_values(cycles)ax.errorbar(s5.cycles,s5[mean],yerrs5[std],markero,labelf{temp}°C, 5% 蔗糖,colorcolors[temp])ax.errorbar(s10.cycles,s10[mean],yerrs10[std],markers,ls--,labelf{temp}°C, 10% 蔗糖,colorcolors[temp])ax.axhline(0,colorgray,lw0.8)ax.set_xlabel(冻融循环数);ax.set_ylabel(粒径漂移 (%))ax.set_title(LNP 冻融稳定性粒径漂移 vs 循环数)ax.legend(fontsize8)plt.tight_layout();plt.show()看图说话的三个判据斜率漂移 vs 循环数近线性 → 应力是累积性的冰晶挤压主导第 1 次循环跳变后平台 → 稳定核不稳定亚群结构分装的瓶间异质性大放行风险高;温度效应-196°C液氮漂移反而常比 -80 小——降温速率越快冰晶越小、停留玻璃化时间短保护剂剂量5%→10% 蔗糖把漂移从 15% 压到 5% 以内比换缓冲液划算得多。五、统计检验到底稳没稳三条判据都是肉眼看图。要不要说5 次循环后粒径无显著漂移得做检验。对数空间漂移、按瓶为独立单位fromscipyimportstats subdf[(df.temp_c-80)(df.sucrose_pct5)]c0np.log(sub[sub.cycles0].groupby(rep)[d_nm].mean())c5np.log(sub[sub.cycles5].groupby(rep)[d_nm].mean())t_stat,pstats.ttest_rel(c0,c5)print(f配对 t 检验t{t_stat:.2f}, p{p:.3f})注意单位是瓶不是读数——同一瓶 3 个读数是技术重复当独立样本会假性显著p 值虚低一个数量级。生物重复才配得上统计推断这是冻融这种每瓶一个命运实验最容易犯的错。五、灌装材料的吸附冻融之外的隐藏变量冻融实验做漂亮了还有一个隐藏变量容器吸附。LNP 对硅胶注射器内壁的吸附损失可到 5–15%加了吐温80或 BSA 的缓冲体系能压下去。冻融研究若在玻璃小瓶做、临床产品在预灌封注射器两者吸附行为不同结论不能直接平移。方法学上至少加一组注射器室温放置 4h的对照。六、小结冻融稳定性 DSC 相变点Tg’/Tm安全温度窗 DLS 循环漂移颗粒完整性两条线都要有数据Z-average 是强度加权对少量大颗粒极敏感双峰分布下均值不可单独引用技术重复 ≠ 生物重复统计单位永远是独立瓶容器吸附是冻融之外常被漏掉的变量预灌封注射器需单独考察。参考Dong Y, et al. Lipid nanoparticles for mRNA delivery. Nature ProtocolsmRNA-LNP 制备、冻存与表征章节《中国药典》通则 0632 热分析法DSC 方法学口径冻干化学药物注射剂药学研究技术原则缓冲液选择与冻干工艺依据DSC 仪器厂商应用笔记生物制剂 Tg’ 测定实例 LNP 冻融稳定性数据分析：Pandas 整理 DSC 相变点 + 冻融循环粒径漂移可视化
惯例跑了一组 mRNA-LNP 冻融实验-20/-80/-196°C 三温度 × 0/1/3/5 次循环DLS 复测粒径和 PDI想看包保护剂蔗糖/海藻糖到底护住了没有。这份笔记用 Pandas 做完整的数据整理scipy 从 DSC 热流曲线提取 Tg’/Tm 相变点最后用 matplotlib 画出粒径漂移 vs 冻融循环的稳定性曲线。依赖pandas、numpy、scipy、matplotlib。一、为什么 LNP 怕冻mRNA-LNP 在 2–8°C 缓慢降解水解、氧化长期储存依赖冷冻。但冷冻本身是个破坏过程冰晶挤压冻结时水结冰排出溶质LNP 被浓缩进非冻结相颗粒间距离从微米级压缩到纳米级聚集乃至融合机械应力冰晶生长的物理挤压直接破坏颗粒结构pH 迁移磷酸盐缓冲液冻结时磷酸氢二钠先结晶残余液相 pH 可漂移 3 个单位枸橼酸和组氨酸缓冲液冻融 pH 稳定得多——这就是冻干制剂几乎全用组氨酸的原因相变脂质核心在低温下从液晶态向凝胶态转变Tm脂肪链排直、晶格收缩把内水相的 RNA 挤出来。所以稳定性数据要同时看两件事热学上有没有危险的相变点DSC 的 Tg’ 和 Tm颗粒上冻融几次后粒径分布漂了多少DLS 复测的粒径漂移曲线。二、数据结构一张长表实验设计 3×4 因子温度 × 循环数 冻干保护剂 2 水平5%/10% 蔗糖 阳性对照无保护剂。DLS 每条件 3 独立瓶 × 3 读数。存成一张 CSVtemp_c,cycles,sucrose_pct,rep,reading,d_nm,pdi -20,0,5,1,1,83.2,0.112 -20,0,5,1,2,84.1,0.115 ...读入和描述统计importpandasaspdimportnumpyasnp dfpd.read_csv(freetz_data.csv)print(df.groupby([temp_c,cycles,sucrose_pct]).agg(d_mean(d_nm,mean),d_sd(d_nm,std),pdi_mean(pdi,mean),n(d_nm,size)))先做正态性检查再谈平均冻融后的粒径分布经常出现主体未变 少量 1–3 µm 大颗粒的双峰。此时 Z-average 被大颗粒拖高强度加权对大颗粒极其敏感d⁶ 加权均值既不代表主体也不代表尾巴。两个习惯报 Z-average 的同时报 PDI 和数量分布峰值对漂移率做 log 变换后再比较对数空间下的漂移近似线性。df[drift_pct]df.groupby([temp_c,sucrose_pct,rep])[d_nm]\.transform(lambdas:(s/s.iloc[0]-1)*100)注意这里按 rep独立瓶分组后取第 0 次循环读数为基线不要用全表 cycles0 的均值当基线——瓶间差异会污染循环效应。三、DSC 热流曲线提取 Tg’/Tm差示扫描量热DSC以固定速率降温/升温热流曲线上的吸热/放热台阶对应相变。Tg’最大冻结浓缩液的玻璃化转变是冻结相从橡胶态进入玻璃态的温度——冻到 Tg’ 以下分子扩散近乎停止这就是安全储存温度的下界。Tm 是脂质凝胶-液晶相变。3.1 一步法一阶导数峰值对升温曲线求一阶导数Tg’ 表现为导数曲线上的台阶Tm 表现为峰。用 scipy 的 savgol_filter 平滑后找峰importnumpyasnpfromscipy.signalimportsavgol_filter,find_peaks# dsc.txt: 两列 温度(°C) 热流(mW)升温段 10 K/minT,hfnp.loadtxt(dsc.txt,delimiter,,unpackTrue)# 平滑 一阶导dhfsavgol_filter(hf,window_length51,polyorder3,deriv1,delta1.0)# Tg: 导数台阶 → 对导数再求导二阶导极值d2savgol_filter(hf,window_length101,polyorder3,deriv2,delta1.0)tg_candidates,_find_peaks(np.abs(d2),heightnp.percentile(np.abs(d2),95))print(Tg 候选:,T[tg_candidates])# Tm: 导数峰tm_candidates,_find_peaks(dhf,height0)print(Tm 候选:,T[tm_candidates])一步法快但 Tg’ 的台阶在原始曲线很微弱尤其冻干样品含糖量高时实际方法学里更稳的做法是 DSC 常用的两遍扫描第一遍升温消除热历史第二遍才是分析曲线。3.2 更稳的两点法半高法对 Tm 峰取峰两侧斜率最大点作切线半高交点定熔点。工程上够用的近似峰起点温度 0.5×(峰结束-起点)。deftm_half_height(T,hf,i_peak,window40):峰位 i_peak 处的半高温度近似lomax(0,i_peak-window);himin(len(T),i_peakwindow)basenp.median(hf[lo:hi])halfbase0.5*(hf[i_peak]-base)# 升温侧最后一个低于半高的位置i_halflonp.where(hf[lo:i_peak]half)[0][-1]returnT[i_half]这段近似代码按需取用严格切线法建议直接用仪器软件本文重点在数据流不在切线几何。四、冻融循环粒径漂移曲线把每个条件下、每次冻融后相对于 cycle 0 的粒径漂移串成曲线importmatplotlib.pyplotasplt g(df.groupby([temp_c,cycles,sucrose_pct])[drift_pct].agg([mean,std]).reset_index())fig,axplt.subplots(figsize(7,4.5))colors{-20:#1f77b4,-80:#ff7f0e,-196:#2ca02c}fortemp,subing.groupby(temp_c):s5sub[sub.sucrose_pct5].sort_values(cycles)s10sub[sub.sucrose_pct10].sort_values(cycles)ax.errorbar(s5.cycles,s5[mean],yerrs5[std],markero,labelf{temp}°C, 5% 蔗糖,colorcolors[temp])ax.errorbar(s10.cycles,s10[mean],yerrs10[std],markers,ls--,labelf{temp}°C, 10% 蔗糖,colorcolors[temp])ax.axhline(0,colorgray,lw0.8)ax.set_xlabel(冻融循环数);ax.set_ylabel(粒径漂移 (%))ax.set_title(LNP 冻融稳定性粒径漂移 vs 循环数)ax.legend(fontsize8)plt.tight_layout();plt.show()看图说话的三个判据斜率漂移 vs 循环数近线性 → 应力是累积性的冰晶挤压主导第 1 次循环跳变后平台 → 稳定核不稳定亚群结构分装的瓶间异质性大放行风险高;温度效应-196°C液氮漂移反而常比 -80 小——降温速率越快冰晶越小、停留玻璃化时间短保护剂剂量5%→10% 蔗糖把漂移从 15% 压到 5% 以内比换缓冲液划算得多。五、统计检验到底稳没稳三条判据都是肉眼看图。要不要说5 次循环后粒径无显著漂移得做检验。对数空间漂移、按瓶为独立单位fromscipyimportstats subdf[(df.temp_c-80)(df.sucrose_pct5)]c0np.log(sub[sub.cycles0].groupby(rep)[d_nm].mean())c5np.log(sub[sub.cycles5].groupby(rep)[d_nm].mean())t_stat,pstats.ttest_rel(c0,c5)print(f配对 t 检验t{t_stat:.2f}, p{p:.3f})注意单位是瓶不是读数——同一瓶 3 个读数是技术重复当独立样本会假性显著p 值虚低一个数量级。生物重复才配得上统计推断这是冻融这种每瓶一个命运实验最容易犯的错。五、灌装材料的吸附冻融之外的隐藏变量冻融实验做漂亮了还有一个隐藏变量容器吸附。LNP 对硅胶注射器内壁的吸附损失可到 5–15%加了吐温80或 BSA 的缓冲体系能压下去。冻融研究若在玻璃小瓶做、临床产品在预灌封注射器两者吸附行为不同结论不能直接平移。方法学上至少加一组注射器室温放置 4h的对照。六、小结冻融稳定性 DSC 相变点Tg’/Tm安全温度窗 DLS 循环漂移颗粒完整性两条线都要有数据Z-average 是强度加权对少量大颗粒极敏感双峰分布下均值不可单独引用技术重复 ≠ 生物重复统计单位永远是独立瓶容器吸附是冻融之外常被漏掉的变量预灌封注射器需单独考察。参考Dong Y, et al. Lipid nanoparticles for mRNA delivery. Nature ProtocolsmRNA-LNP 制备、冻存与表征章节《中国药典》通则 0632 热分析法DSC 方法学口径冻干化学药物注射剂药学研究技术原则缓冲液选择与冻干工艺依据DSC 仪器厂商应用笔记生物制剂 Tg’ 测定实例