编辑推荐:
研究人员研究了受限泊肃叶流动(Poiseuille flow)中二维多组分囊泡的动力学,聚焦于膜组分、力学异质性与流体动力学约束之间的相互作用。膜组分采用相场模型(phase-field model)描述,局部弯曲刚度(bending rigidity)依赖于
研究人员研究了受限泊肃叶流动(Poiseuille flow)中二维多组分囊泡的动力学,聚焦于膜组分、力学异质性与流体动力学约束之间的相互作用。膜组分采用相场模型(phase-field model)描述,局部弯曲刚度(bending rigidity)依赖于相组分,从而将相分离与膜变形耦合。通过拉格朗日乘子(Lagrange multiplier)施加局部膜不可伸展性,并在内外黏度匹配条件下采用边界积分方法(boundary integral method)求解囊泡-流体耦合系统。在固定的五倍弯曲刚度对比下,研究人员考察了平均相组分、通道宽度和初始囊泡取向的影响。与通常弛豫至通道中心线附近对称子弹形构型的均匀囊泡不同,多组分囊泡表现出显著的对称性破缺和多种动力学状态,包括伴随膜坦克履带运动(tank-treading)的不对称稳态形状,以及依赖于受限程度的呼吸(breathing)和摆动(wagging)运动。力学异质性还会在瞬态变形过程中促进局部高曲率区域的形成,且较软的膜结构域优先占据曲率较大的区域。这些结果表明,与组分相关的弯曲刚度能够定性地改变压力驱动受限流动中的囊泡动力学,并凸显了膜异质性在窄通道内囊泡样细胞输运与变形中的重要性。
多组分囊泡在受限微通道中的流动是数学生物学中的经典动边界问题,并与镰状细胞病(sickle cell disease, SCD)的微循环功能障碍密切相关。SCD患者红细胞内血红蛋白S异常聚合导致细胞密度、形态和变形能力出现显著异质性,进而在直径5–15 μm的毛细血管中诱发异常变形与运动,增加流动阻力与血管黏附风险。已有研究广泛采用边界积分方法模拟均匀囊泡在泊肃叶流中的形态与迁移,揭示子弹形、降落伞形、拖鞋形等多种稳态;相场模型则被用于描述膜上脂相分离、出芽等形态过程。然而,多组分囊泡在直通道压力驱动流中的动力学仍缺乏系统刻画,特别是膜弯曲刚度随组分变化与流体动力受限之间的耦合机制尚不清楚。为此,研究人员建立了受限泊肃叶流中多组分囊泡的完全耦合数学模型,并采用稳定数值格式考察膜力学异质性对囊泡动力学的影响。
研究人员开展了以下研究。在建模方面,采用相场变量φ描述膜局部组分,以双阱势能的线能项描述相分离;膜弯曲能采用含组分依赖弯曲刚度B(φ)的Helfrich能量,并取B(φ)随φ线性插值,其中较软红色相为2×10
?19 J,较硬蓝色相为5倍参考值。局部不可伸展性通过拉格朗日乘子施加。流体问题用边界积分法求解周期性通道内单位胞元中的斯托克斯方程,利用近场镜像与代理源(proxy sources)构造远场修正,并通过扩展线性方程组同时施加通道壁无滑移与周期压力降条件。相场演化采用Cahn-Hilliard型方程,数值上结合Fourier伪谱方法与小尺度分解,用二阶指数时间差分推进。所有模拟使用512个膜离散点。该模型忽略了细胞骨架剪切弹性、胞内血红蛋白聚合和三维几何,作为隔离弯曲刚度效应的一种理想化模型。研究无样本队列来源。
研究结果如下。
3.1 相场演化基准测试:在固定细长囊泡几何、仅演化膜组分的基准算例中,φ趋向两个纯相,界面宽度较窄,并形成四段相边界;较硬的相优先占据低曲率区域,较软相趋向高曲率区域,表明弯曲刚度与曲率的耦合可以驱动相的空间分选。
3.2 通道中的均匀囊泡:与参考弯曲刚度相比,五倍弯曲刚度改变瞬态变形和长期构型。参考刚度下三种初始取向均弛豫至中心线对称子弹形;高刚度下水平初始取向演化为前缘膨大的哑铃形,而倾斜和垂直取向仍为子弹形,说明均匀囊泡的终态形状总体上保持对称性。
3.3 通道中的多组分囊泡:在固定五倍弯曲刚度对比下,系统考察平均相浓度、初始取向和通道宽度的影响。
3.3.1 两相囊泡:当硬相为少数相时,无论初始取向如何,囊泡最终收敛到同一不对称子弹形,刚性蓝色相位于前缘较平区域;瞬态变形中产生局部高曲率尖点。当两种相比例接近时,不同初始取向分别导致水平泪滴形或不同子弹形等长期状态,表明存在多个动力学状态;平均相浓度较高时也出现带明显尾尖的泪滴形或不对称形状,且相分布沿膜持续漂移。
3.3.2 不同通道宽度:对均匀囊泡,缩窄通道增强壁面约束,形成壁面引导构型;通道变宽后恢复子弹形。对多组分囊泡,较宽通道中出现伴随缓慢坦克履带运动的不对称子弹形;通道变窄时,大幅形状重排被抑制,坦克履带运动仍在但构型更受限。当平均相浓度增加时,中等受限程度出现呼吸运动,形状、质心和相分布呈周期振荡;最宽通道出现摆动运动,取向和横向位移周期性变化。这些结果表明受限程度与相分数共同控制动力学模式转变。
讨论:该研究的主要贡献在于揭示膜组分-曲率-流体动力耦合作为窄通道内囊泡形态与运动的控制机制。与均匀囊泡相比,异质膜通过相分布改变局部弯曲刚度,产生对称性破缺、曲率局部化和耦合相再分布,从而出现不对称构型及呼吸/摆动运动。这些结果对红细胞和SCD的关联是定性的,模型未量化水动力学阻力、细胞捕获或血管闭塞;多细胞聚集导致血管阻塞仍是需未来研究的假说。已有微流控和计算研究支持异常红细胞力学与微循环输运受损之间的定性联系,但尚不能证明本文观察到的高曲率形态本身足以增加阻力和捕获。模型局限包括内外黏度相等、二维几何、局部不可伸展膜和线性弯曲刚度插值;未来应引入黏度对比、三维几何、细胞骨架弹性、膜张力调控及真实微血管网络。
研究结论:与组分相关的弯曲刚度能够定性地改变压力驱动受限流动中的囊泡动力学;膜力学异质性在窄通道内囊泡样细胞的输运与变形中具有重要作用。