关于电穿孔模型及其跨膜分子转运机制需重新审视的学术报告

一、 论文基本信息 本文由斯洛文尼亚卢布尔雅那大学电气工程学院的 Maria Scuderi、Janja Dermol-Černe、Lea Rems*（通讯作者）与荷兰代尔夫特理工大学化学工程系的 Clarissa Amaral da Silva、Aswin Muralidharan、Pouyan E. Boukany 共同完成，于2022年7月25日在线发表于爱思唯尔（Elsevier）旗下的学术期刊《Bioelectrochemistry》（第147卷，文章号108216）。这是一篇开放获取文章，采用CC BY-NC-ND 4.0许可协议。

二、 学术背景与研究目标 本研究属于生物电化学与生物物理学交叉领域，聚焦于电穿孔（Electroporation，亦称电渗透）技术的核心理论建模问题。电穿孔是一种利用高强度脉冲电场暂时增加细胞膜通透性的技术，已广泛应用于基因电转染、电化学疗法、药物输送及食品加工等多个领域。为了优化电脉冲参数、预测分子递送效率并深入理解其物理机制，科学家们开发了多种多物理场、多尺度的数学模型，用于在单细胞水平上模拟电穿孔过程及随后的分子跨膜转运。

然而，一个引人深思的现象是：尽管这些模型在理论描述上存在显著差异（例如，对膜通透性增加的机制假设不同），但它们各自在与其原始开发数据对比时，往往都报告了与实验的良好一致性。这引发了一个根本性问题：这些不同的模型不可能同时完全正确。那么，是电穿孔分子转运过程对模型细节不敏感，还是模型开发和验证方式存在普遍问题？

基于此，本研究设定了明确目标：1）批判性地审视现有的、允许进行时空分辨数值模拟的机理性单细胞电穿孔模型；2）通过将这些模型与一系列广泛的、使用不同脉冲参数和小分子的实验测量数据进行对比，评估其中任何一个模型是否足够“普适”，能够普遍描述不同实验条件下的跨膜转运现象；3）揭示现有模型的局限性，并为未来模型的验证和改进提供具体的指导方针和建议。

三、 研究流程与方法详述 本研究主要采用理论分析、数值模拟与实验数据对比验证相结合的研究范式，并未进行新的湿实验，而是对现有模型和已发表的实验数据集进行系统性评估。其工作流程可详细分解如下：

第一步：文献调研与模型筛选。 研究团队首先在Web of Science平台上进行了系统的文献检索，筛选出截至2021年底发表的、描述动物单细胞水平电穿孔及小分子转运的数学模型。他们设定了严格的入选标准：1）模型需描述电穿孔介导的小分子/离子摄取；2）能进行时空分辨的数值模拟；3）必须是机理性模型（即基于膜结构变化的分子机制）；4）已与实验数据（定性或定量）进行过比较。最终，从23篇候选论文中筛选出8个完全符合标准的模型。

第二步：模型分类与代表模型选取。 根据对膜通透性增加的理论描述方式，他们将现有模型分为三类：(i) 孔态模型（Pore states models）：基于动力学方案，描述膜在完整态与多个多孔态之间的转变。(ii) 孔分布模型（Pore distribution models）：基于Fokker-Planck或Smoluchowski方程，用分布函数描述不同半径的孔在电场下的动态演化。(iii) 其他现象学模型（Phenomenological models）：使用强经验性函数描述通透性变化，缺乏明确的分子机制基础。由于现象学模型严重依赖于拟合特定数据集，本研究不予进一步测试。从前两类中，他们选取了三个具有代表性且被广泛引用或作为后续研究基础的机理模型进行深入测试： 1. MT2010模型（Miklavčič and Towhidi, 2010）：孔态模型的代表，引入了“记忆态”来描述长期的扩散和内吞样转运。 2. S2011模型（Smith, 2011）：孔分布模型的重要分支，详细考虑了分子尺寸、电荷、在孔内的阻碍效应（hindrance）及分配效应（partitioning）对转运的影响。 3. LL2011模型（Li and Lin, 2011）：另一重要的孔分布模型，基于Krassowska和Filev的框架，主要关注电泳和扩散转运。

第三步：模型实施与对比实验选择。 研究团队独立重新实现了这三个选定的模型（代码已在GitHub开源）。为了进行公正的测试，他们特意选择了这些模型在其原始论文中未曾对比过的实验数据集，以检验其外推和普适能力。选用的关键实验数据包括： - Puc等人（2003）测量的DC-3F成纤维细胞在单次100 μs或1 ms矩形脉冲后，荧光染料荧光黄（Lucifer Yellow）的摄取量。 - Canatella等人（2001）测量的DU 145前列腺癌细胞在单次指数衰减脉冲后，染料钙黄绿素（Calcein）的摄取量。 - Sozer等人（2018）测量的U937淋巴瘤细胞在单次220 μs脉冲或10次6 ns脉冲后，钙黄绿素和碘化丙啶（Propidium Iodide）摄取的动力学过程（时间进程）及空间分布（阳极侧与阴极侧的不对称性）。 - Gabriel和Teissié（1999）测量的CHO细胞中钙离子（与Fluo-3结合）摄取的二维荧光图像。

第四步：模型与实验的定量/定性对比分析。 这是本研究的核心环节。团队将每个模型在对应实验条件下的模拟结果（如最终细胞内浓度、摄取动力学曲线、分子在细胞两侧的分布）与上述实验测量值进行细致比较。他们调整了模型参数（如细胞半径、介质电导率、脉冲波形）以匹配实验条件，但刻意避免为了拟合新数据而重新优化模型的核心参数（如孔扩散系数、孔形成速率等），旨在检验模型的“即用”普适性。

第五步：结果分析与局限性探讨。 基于对比结果，作者深入分析了每个模型在解释不同实验数据时暴露出的矛盾与不足，并追溯这些不足背后可能的理论假设缺陷。

四、 主要研究结果详述 研究结果通过几个精心设计的对比案例，清晰地揭示了现有主流电穿孔模型的局限性：

1. 定性验证不足以确立模型的有效性：以MT2010模型为例。 该模型最初仅与Kotnik等人（2003）的荧光黄摄取实验进行了定性比较（如预测矩形脉冲比三角波、正弦波脉冲导致更多摄取），并显示一致。然而，当本研究将其与Puc等人（2003）的定量数据对比时（图2），发现尽管模型能定性地重现“摄取量随电场增强而增加”的趋势，但其预测的细胞内荧光黄浓度比实验测量值高出近一个数量级。这表明，仅凭定性吻合远不能证明模型的定量预测能力或机制的正确性。

2. 固定时间点的定量验证无法反映转运动力学的正确性：以S2011模型为例。 S2011模型在其原始论文中，通过调整多个参数，与Puc等人（荧光黄）和Canatella等人（钙黄绿素）在脉冲后某一固定时间点测得的分子摄取总量数据达成了极佳的定量吻合（图3）。这曾是该模型被认为强大的重要依据。然而，当将其应用于Sozer等人测量的动力学过程数据时，问题显现。对于220 μs脉冲（图4）： - 动力学失配：模型严重高估了脉冲期间的电泳转运，模拟显示脉冲施加时（t=5s）细胞内浓度瞬间骤升，而实验数据中未见此剧烈跳跃。 - 空间不对称性矛盾：对于带正电的碘化丙啶，实验显示其在阴极侧的摄取多于阳极侧，这暗示扩散为主导机制且阴极侧通透性更高。而原始S2011模型（假设电泳主导）预测了相反的趋势（阳极侧更多）。 - 参数调整的随意性：作者演示了如何通过“合理”但任意的调整（如忽略电泳、增大最小孔半径）使模型结果与实验吻合（图4c）。这恰恰说明复杂模型通过调整多个参数可以拟合多种数据，但这并不能证明模型机制的正确性，反而揭示了其潜在的“过拟合”风险和参数不确定性问题。 对于6 ns脉冲序列（图5），S2011模型同样出现定量偏差（预测的浓度变化比实验值高3-10倍），并且未能正确预测多次脉冲相对于单次脉冲的累积效应（模型预测相似，而实验显示多次脉冲效果强得多）。

3. 某些实验对模型验证缺乏区分度：以LL2011和S2011模型对比为例。 LL2011模型最初因其能定性重现Gabriel和Teissié的钙离子摄取二维图像而获得验证（图6a-c）。然而，本研究发现，尽管LL2011和S2011模型在孔动力学和分子通量描述上存在差异，但它们都能很好地模拟该实验的钙离子空间分布模式。这意味着该实验无法区分这两个理论不同的模型。而当用LL2011模型去模拟Puc等人的荧光黄定量数据时，它却严重高估了摄取量（图6d）。这说明，仅依靠单一类型（尤其是定性）的实验数据来验证模型是远远不够的，容易导致对模型普适性的错误信心。

4. 模型根本局限源于对电穿孔分子机制的简化假设。 作者指出，当前几乎所有机理模型都基于一个核心假设：膜通透性的增加完全归因于脂质孔的形成与关闭。然而，越来越多的证据表明这一假设可能过于简化： - 长寿命孔的能量困境：为解释实验中持续数分钟的物质转运，模型需假设脂质孔关闭存在极高的能量壁垒（>10 kT），这与纯脂质体系的自由能计算和实验观测不符。 - 其他潜在机制：实验和模拟表明，电脉冲还可能引起脂质过氧化损伤、膜蛋白变性、以及细胞骨架破坏等，这些都可能形成性质不同于纯脂质孔的、更持久的膜缺陷或孔洞。Pavlin等人早先也提出了“短寿命孔”和“长寿命孔”的区分。这些机制可能共同作用，但现有模型均未纳入考虑。

五、 研究结论与价值 本研究得出核心结论：目前没有一个现有的单细胞电穿孔模型能够普适地描述所有不同脉冲参数和小分子下的跨膜转运实验数据。 即使通过参数调整能使某个模型匹配特定数据集，其内在的物理机制描述也可能是不完整或不准确的。模型在应用于其原始验证条件之外时需格外谨慎，必须辅以实验验证。

本文的科学价值与应用意义在于： 1. 批判性评估与警示：首次系统性地对主流电穿孔模型进行了“交叉检验”，揭示了它们在普适性上的严重局限，挑战了领域内可能存在的对某些模型的过度信任。 2. 指明了模型发展的未来方向：明确指出，仅仅对现有模型进行参数优化或方程修补可能走入死胡同。根本出路在于更深入地理解电穿孔的分子机制。必须超越“单一脂质孔”的范式，考虑多种可能的膜缺陷类型（如氧化损伤孔、蛋白质介导的孔）及其协同作用。这需要结合分子动力学模拟、增强采样方法以及对包含复杂脂质、蛋白质和细胞骨架的模型膜系统的实验研究。 3. 为实验验证提供了具体指南：为确保未来开发的模型能得到有效验证，作者提出了设计验证实验的详细建议，强调实验应提供：定量的分子摄取/泄漏数据、摄取的时间动力学曲线、显示不对称性的二维空间分布信息，并且要覆盖宽范围的脉冲参数（时长、幅度、数量）和不同电荷的分子（至少一种阴离子和一种阳离子）。同时，实验报告必须详尽，包括细胞尺寸形状、密度、介质电导率、脉冲波形、分子特性、温度等关键参数。

六、 研究亮点 1. 研究视角新颖：不是提出新模型，而是对现有主流模型进行“第三方”的、系统性的比较评估，发现了模型宣称的“成功”与其实验普适性之间的巨大差距。 2. 论证逻辑严密：通过将模型与其未训练过的实验数据进行对比，有力揭示了模型的局限性。特别是利用Sozer等人的动力学和空间分辨数据进行检验，比传统的终点法数据更能暴露模型在机制描述上的问题。 3. 深入机制批判：不仅指出了模型与数据的不符，更进一步追溯到当前理论框架的核心假设缺陷（即仅考虑脂质孔），并提出了结合多尺度研究探索多元机制的未来路径。 4. 实践指导性强：不仅指出问题，还为建设性地改进模型发展和实验验证提供了清晰、可操作的具体建议和数据共享倡导（如开源代码），对领域发展具有直接的推动作用。

七、 其他有价值的内容 作者在文中强调了数据与模型共享的重要性。他们指出，许多已发表模型的描述不够详细或存在印刷错误，阻碍了重现和使用。为此，他们将自己在本研究中实现的所有模型代码公开在GitHub仓库中，为领域的透明化和可重复性研究树立了榜样。这一行动本身也是对科学共同体的一种重要贡献。