干细胞无标记表征中电穿孔与阻抗谱混合建模研究

《Scientific Reports》：Hybrid modeling of electroporation and impedance spectroscopy for label free characterization of stem cells

　　

在再生医学与细胞治疗中，对干细胞分化状态进行无标记、非破坏性表征仍是一项重要目标。现有计算框架常将电穿孔（electroporation）置于组织尺度或简化单细胞几何中处理，较少有研究将时域电穿孔可观测量与微流控平台上测得的扫频阻抗特征相联系。本研究提出一种修正的混合解析–数值与实验框架，用于比较未分化人源间充质干细胞（human mesenchymal stem cells, hMSCs）与成骨定向 hMSCs。数值模型采用细胞类型特异性参数：未分化 hMSC 模型代表半径 RU=10 μm、胞质电导率 σi,U=0.32 S m-1、膜电容 Cm,U=1×10-2 F m-2、特征电穿孔电压 Uep,U=0.258 V；成骨 hMSC 模型 RO=13 μm、σi,O=0.24 S m-1、Cm,O=8×10-3 F m-2、Uep,O=0.32 V。两模型置于同一微流控电极环境，受 1–5 kV cm-1、上升时间 1 ns 的电场脉冲激励。被动 Schwan RC 时间常数分别为 4.69×10-7 s 与 5.96×10-7 s；阈值后上升时间更短，约 1×10-7–2×10-7 s。5 kV cm-1 下被动极化跨膜电位分别约 7.5 V 与 9.75 V。对两群体行 1 kHz–1 MHz 扫频阻抗谱（impedance spectroscopy），在 1、5、10、15、20、25 V 外加电压下提取低频阻抗、串联电阻、电抗谷深、相角与电压依赖阻抗降。实验显示成骨 hMSCs 基线阻抗更高（1 V、1 kHz 下 |Z|=11029 Ω vs 7863 Ω），未分化 hMSCs 在 25 V 下高压阻抗降更强（约 92.4% vs 86.9%）。以 10.4 μm 与 24.9 μm 聚苯乙烯微球作无细胞尺寸标准。联合结果定义细胞类型特征空间 Fcombined={τcharge, Vm, N(t), rp(t), σm(t), fc, Δ|Z|fc,0, Rs,1kHz, |Xs,pk|} 供未来无标记分类研究使用。本工作应理解为匹配的建模与阻抗分析框架；确定性生物学分类、直接孔成像、活性验证及训练分类器性能不在本研究范围。

研究背景方面，电穿孔作为外电场诱导脂质双分子层瞬时通透的技术，在生物医学中具基石地位，其过程涵盖膜电容充电、孔核化与孔扩张重封三阶段， efficacy 受场强、脉宽、细胞生物物理属性及介质电导率共同调控。纳米秒脉冲电场（nsPEFs）可扰动质膜与细胞内膜且限热。在干细胞领域，nsPEFs 还被报道可调节分化潜能。现有计算研究多孤立处理单细胞膜模型或组织尺度电穿孔，阻抗法虽已用于干细胞分化监测，但三条鸿沟依然存在：其一，尚无框架在微流控几何中同时解析跨膜电位、孔密度、孔半径与动态电导率四大量；其二，孔半径动力学与阻抗谱解释工作流鲜有显式关联；其三，超越 V_m 与 N 的量化特征不足以区分同脉冲下不同分化态干细胞。该研究即针对上述缺口，在《Scientific Reports》发表，建立匹配建模与阻抗分析框架。

关键技术方法上，研究人员构建球形单壳（single-shell）细胞模型，膜作薄边界分布阻抗条件，耦合拉普拉斯准静态电位方程与 Smoluchowski 型孔密度方程及 Krassowska–Filev 孔半径动力学常微分方程，于 COMSOL Multiphysics 用 BDF 隐式时间步求解，网格 45210 单元、膜边界 0.05 μm。样本为文献代表性参数化的未分化 hMSCs（R_U=10 μm）与成骨定向 hMSCs（R_O=13 μm），实验端在 PDMS 微流控通道共面电极（间隙 200 μm）中以 PBS 悬浮细胞，施加 1 ns 上升阶跃脉冲（1–25 V 对应约 50–1250 V cm^-1），以 10.4 μm 与 24.9 μm 聚苯乙烯微球校准，随后 1 kHz–1 MHz 扫频阻抗谱提取 |Z|、R_s、X_s。仿真以 1–5 kV cm^-1 场强扫描作机理参扫。

研究结果部分，跨膜电位分布小节指出，未分化模型被动 Schwan RC 时间 4.69×10^-7 s、5 kV cm^-1 被动极电位 7.5 V，成骨模型 5.96×10^-7 s、9.75 V；过阈后两模型极电位钳位于约 ±1.4 V，孔介导电导限制进一步充电。孔密度演化小节显示，最高孔密度位于极化两极，未分化极密度约 3.6–3.7×10¹⁴ m^-2 略高于成骨 3.4–3.5×10¹⁴ m^-2，但成骨中间/边缘段平台更高。孔半径表征小节中，早期窗未分化可见至 9.2 nm、成骨至 7.9 nm，均随位置与场强非单调扩展。膜电导率动力学小节表明，σ_m 峰值量级 10^-4 S m^-1，沿膜周长呈双峰极区分布。膜重封行为小节指出，脉冲止后 V_m、N、σ_m 渐进衰减，成骨初期电导恢复略快，留残量差异，模型预测可逆重封但无活体染色验证。定量模型–数据比较小节将 σ_m(t) 增升与实验中 25 V 下 |Z| 降 92.4%（未分化）和 86.9%（成骨）物理关联，敏感性分析显示 τ_charge 对中 C_m、σ_e 中等敏感，钳位电压与 σ_m,max 扰动下变异小。

阻抗实验表征小节报告，1 V、1 kHz 下成骨 |Z|=11029 Ω 高于未分化 7863 Ω，R_s 相应 10804 Ω vs 7687 Ω，|X_s,pk| 1921 Ω vs 1609 Ω；≥15 V 出现对比反转，未分化低压阻抗更低（15 V 时 3409 Ω vs 7563 Ω），25 V 时未分化降至 599 Ω、成骨 1448 Ω，Δ|Z|_1kHz 达 92.4% vs 86.9%；β 色散角频随电压上移，基线交叉频 f_c,0 未分化 5.6 kHz、成骨 5.8 kHz。讨论部分确认时域 σ_m(t) 增升与频域低频 |Z|、R_s、|X_s 坍塌同源，亚阈特征由尺寸与容抗标定主导，超阈轨迹反映分化相关膜组成差异；九维特征空间 F_combined 统合两域供后续机器学习分类用。

结论部分翻译：本研究提出混合解析–数值框架解释 hMSC 微流控实验中电穿孔动力学与扫频阻抗谱，模型解析两分化态的 V_m、N、r_p、σ_m，阻抗测量提供电压依赖膜通透的频率域佐证。修订解读区分被动充电标度与过阈电压钳制，过阈后孔介导电导限压于约 1.3–1.4 V。孔密度、孔半径与 σ_m 作机制描述符而非已验证生物标志物。阻抗实验显示 |Z|、R_s、|X_s 与相角随电压降减，与模型 σ_m 增升同径；未分化 25 V 阻抗降 92.4% 强于成骨 86.9%，微球校准同几何尺寸对比。F_combined={τ_charge, V_clamp, N_max, r_p(t), σ_m,max, f_c, Δ|Z|_fc,0, R_s,1kHz, |X_s,pk|} 作无标记表征候选描述符。本研究不含训练分类器、交叉验证、孔直接成像或活性 assay，应作未来验证基础而非完结分类法。局限性含单壳球近似、忽略核与细胞骨架、孔参数由动力学推演非直视、仿真与实验电场非全等故属机制一致而非逐变量定量验证。