Skip to content

Karkhaneh Yousefi, Ali Akbar, Petit, Claudie, Ben Hassine, Amira et al. · 2023 · J. Mech. Behav. Biomed. Mater.

#continuum-mechanics #finite-element #vascular-G&R

??? info "Metadata" DOI: 10.1016/j.jmbbm.2023.105990

📄 Quick Note

摘要概述

本文提出一种基于连续介质力学的新型计算模型，用以描述主动脉平滑肌细胞 (aortic SMCs) 通过 acto-myosin 收缩实现 stiffness sensing 的机制。核心创新是将每根 stress fiber 的 myosin 激活等效为"热应变" (eigenstrain αΔT)，并引入双段应力-应变关系（线性段 + softening 段，后者刻画 actin/myosin 丝滑动），仅用 5 个本构参数即可刻画单根纤维的力学行为。作者用 Abaqus 建立含随机数量/排布 stress fiber、Weibull 分布激活水平的细胞群体有限元模型，并将预测与两种 SMC谱系 (healthy ASMC_1、aneurysmal ASMC_2) 的 traction force microscopy 实验对照，发现模型不仅准确预测了 traction 随基底刚度的两段式响应（threshold ≈ 12 kPa），还能再现细胞间统计分散。模型进一步预测基底刚度通过 cytoskeleton 直接传递至核膜与细胞核，引起核形变，提示力学-基因表达耦合路径。

关键图表

  • Figure 1（stress fiber 应力-应变示意）：纤维先以斜率 E_SF 线性加载至 σ_max，超过该弹性限后进入 softening 段（斜率 E_SF'），用以表征 actin/myosin 丝在最大收缩后的滑动效应——是整篇文章的核心本构假设。
  • Figure 5/7（ASMC_1 / ASMC_2 的 traction vs. 基底刚度）：模型预测与 TFM 实验点高度吻合，尤其带 softening 的版本在大刚度区正确再现了 traction 不再单调上升甚至下降的"optimal stiffness"现象；去掉 softening 则在高刚度区明显偏离。
  • Figure 9（核膜力学性质敏感性）：核膜刚度变化对细胞总体 traction 影响很小，但核内与核膜 von Mises 应力呈反向变化——揭示 cytoskeletal 力可经 perinuclear actin cap 直接影响核力学环境，潜在影响基因表达。

与我的关联

该工作把细胞尺度的 stiffness sensing 简化为可在标准 FE 框架内实现的本构关系（eigenstrain + softening），与 vascular G&R 多尺度建模天然衔接——其 stress fiber 本构和"细胞核应力传递"结论可直接作为血管壁力学-生物学耦合的下游模块，值得在 SMC 主动张力建模中借鉴。

🔬 Deep Note

摘要概述

Karkhaneh Yousefi 等（Mines Saint-Etienne, SAINBIOSE INSERM）提出了一个基于连续介质力学的主动脉平滑肌细胞（SMC）刚度感知模型。核心创新是用热应变的类比（eigenstrain αΔT）建模 stress fiber 的 acto-myosin 收缩：α 为收缩系数，ΔT 为 myosin 激活水平（0-1，满足 Weibull 分布）。每根 stress fiber 满足双段应力-应变关系（弹性段 + 软化段，后者描述 actin/myosin 滑移）。模型用 Abaqus FE 实现，单细胞含随机数量的 stress fiber（均值 25），在 4/8/12/25 kPa 四种基底刚度上模拟。结果表明模型不仅准确预测了牵引力随基底刚度的变化（12 kPa 达峰值，即"最优刚度"），还成功复现了细胞间牵引力的统计变异。最后将模型用于计算核膜与细胞核应力，显示基底刚度引起的骨架力变化直接导致核形变，可能改变基因表达。

研究问题与动机

主动脉 SMC 通过 acto-myosin 收缩感知环境刚度，但连续介质力学建模这一过程的工作极少——既有 motor-clutch 模型参数过多，且未在 FE 框架内考虑细胞几何与 stress fiber 排列。作者要回答：能否用最少的参数（5 个本构参数）在连续介质力学框架下预测 SMC 的刚度感知行为，并同时解释细胞间统计变异？

方法

核心本构——stress fiber 的"热应变类比"： - stress fiber 视为 truss 元素，myosin 激活产生 eigenstrain αΔT（类比热收缩 ΔT=温差，α=线膨胀系数） - ΔT ∈ [0,1] 满足 Weibull 分布：f(ΔT) = (k/λ)(ΔT/λ)^(k-1) e^(-(ΔT/λ)^k)，k=形状参数，λ=尺度参数（由中位 ΔT̅ 决定） - 双段应力-应变： - 弹性段（ε ≤ ε_max）：σ = E_SF (αΔT − ε_a)，E_SF = 50 MPa - 软化段（ε > ε_max）：σ = E'_SF (αΔT − ε_a) + (E_SF − E'_SF) ε_max，E'_SF < 0 描述 actin/myosin 滑移 - ε_max = ξ αΔT（弹性极限与激活水平成正比，ξ 为线性系数 0<ξ<1）

FE 模型： - 半纺锤形细胞 250×30 µm²，含细胞膜、细胞质、stress fiber（~25 根，均匀分布 N~U(25, 6.05²)）、核膜、细胞核 - 一半 stress fiber 连接两端 focal adhesion（跨核），另一半连接核膜（perinuclear actin cap，LINC 复合体） - 细胞膜/核膜：neo-Hookean，c10=600 kPa - 细胞质/核：neo-Hookean，c10=100 Pa - 基底：线弹性，E=4/8/12/25 kPa，ν=0.45

参数标定：TFM 实验，两个细胞系（ASMC_1=商业 30 岁女性；ASMC_2=ATAA 72 岁男性原代），每系 ~160 细胞 × 4 刚度。Weibull 参数 (k, ΔT̅) 由标定化牵引力分布拟合；(α, ξ, E'_SF) 由 FE 拟合 4/8/12 kPa（25 kPa 用于估计软化）。

5 个本构参数：E_SF（固定 50 MPa）、ΔT（Weibull 分布）、α（收缩系数）、ξ（弹性极限系数）、E'_SF（软化系数）。

结果

标定的参数值（Table 1）： | 参数 | ASMC_1 | ASMC_2 | |---|---|---| | k | 1（指数分布） | 2（Rayleigh 分布） | | ΔT̅ | 0.132 | 0.363 | | α | 0.026 | 0.031 | | ξ | 0.68 | 0.68 | | E'_SF (MPa) | -17.4 | -80.2 |

关键发现： 1. 12 kPa 为"最优刚度"——牵引力达峰值；>12 kPa 后下降（softening 段生效，actin/myosin 滑移） 2. 软化行为是关键——无软化的模型（ξ=1）在高刚度基底上过估牵引力；含软化的模型与实验高度一致 3. 统计变异成功复现——虚拟 300 细胞群体（随机 stress fiber 数 + 随机 ΔT）的牵引力分布与实测分布吻合 4. 核力学耦合——基底刚度 → 骨架力变化 → 核膜/核形变 → 潜在基因表达改变（mechanotransduction 到核） 5. ASMC_2（ATAA 来源）更高激活（ΔT̅=0.363 vs 0.132）但更低软化系数（E'_SF=-80.2 vs -17.4）——提示动脉瘤 SMC 的 acto-myosin 动力学改变

讨论

作者立场：热应变类比是最简洁的 SMC 刚度感知连续介质力学表达——无需 motor-clutch 模型的数十个参数，5 个本构参数即可预测单细胞 + 群体统计。

局限性： - 2D 培养基底，非 3D ECM 环境（作者承认需扩展到 3D） - 单细胞模型，未考虑细胞-细胞相互作用 - ΔT 的 Weibull 分布是唯象拟合，无分子机制解释 - 核形变→基因表达的链条是推测，未实验验证 - 线弹性基底，未考虑 ECM 非线性力学

个人批注

与我 SMC G&R 工作的对接： 1. SMC 收缩力作为 G&R 的驱动信号——constrained mixture 框架中 SMC 成分的 deposition stretch 和 production rate 需要力学输入信号。本文的 αΔT eigenstrain 提供了一个从 ECM 刚度 → SMC 收缩力 → mechanotransduction 的连续介质力学桥梁。可考虑将 αΔT 作为 G&R 模型中 SMC active stress 的输入。 2. "最优刚度" 12 kPa——这是 SMC 力学稳态的 target 值的候选。Taylor & Humphrey 2009 说"壁厚倾向于将环向应力维持 target value"，本文提供了 SMC 感知刚度的微观基础——G&R 模型可设 homeostatic stiffness ≈ 12 kPa 区间。 3. 软化行为与 G&R 的连接——当 ECM 刚度超过 12 kPa，SMC 牵引力反而下降（actin/myosin 滑移）。这意味着 G&R 模型中 SMC 的 active tone 不能简单随刚度单调递增——需要分段：弹性段（<12 kPa）+ 软化段（>12 kPa）。 4. 核力学耦合——mechanotransduction 到基因表达的物理路径。这对 G&R 模型中"SMC 表型切换"（Bonetti 2021 的 KLF-4 通路）提供了力学入口：刚度 → 核形变 → 基因表达 → 表型切换。

可复用点： - eigenstrain αΔT 概念可直接在 FEniCS 中实现为 stress fiber 的 active strain 分量 - Weibull 分布的 ΔT 可作为 SMC 群体异质性的模型 - 5 参数本构（E_SF, ΔT, α, ξ, E'_SF）足够简洁，可嵌入 G&R 框架

疑问： - 12 kPa "最优刚度"在体内对应什么血管/什么病理状态？体内血管刚度通常更高（~100 kPa 量级），这个差异是 2D 培养基底的局限吗 - ATAA 来源 ASMC_2 的更高 ΔT̅ 是疾病原因还是结果 - 如何将单细胞模型扩展到血管壁中的 SMC 群体（考虑细胞间 ECM 介导的力学耦合）

与上下文的关系

本文在引用网络中的位置： - 建立在：团队前期 TFM 实验（Petit 2022 [33], AFM stiffness mapping [61]）、motor-clutch 模型（Bangasser 2013 [45], Isomursu 2022 [46]） - 被后续引用：作为 SMC 刚度感知的简洁连续介质力学模型，被血管生物力学和 mechanotransduction 领域引用

与同类工作的对比： - vs motor-clutch 模型：参数更少（5 vs 数十个），但牺牲了分子细节；适合嵌入组织级 G&R 模型 - vs 2009-VascularBiomech-Taylor（本队列）：Taylor&Humphrey 给宏观 G&R 框架（constrained mixture + FSG），本文给 SMC 力学感知的微观模型——二者可嵌套：本文的 αΔT 作为 Taylor&Humphrey 框架中 SMC active tone 的子模型 - vs 2021-NitroRedox-Atherosclerosis-Bonetti（本队列）：Bonetti 给 NO → RhoA → 表型切换的分子通路，本文给 RhoA 下游的 acto-myosin 力学——二者互补：分子通路 → 力学执行器

与我队列中其他论文的关系： - 2009-VascularBiomech-Taylor：宏观 G&R 框架 vs 本文的微观 SMC 模型——可嵌套 - 2021-NitroRedox-Atherosclerosis-Bonetti：NO→RhoA 分子通路 vs 本文的 RhoA→acto-myosin 力学输出 - 2016-MorphoelasticWound-Bowden：morphoelastic growth 模型可借鉴本文的 eigenstrain 概念描述 active contraction - 2023-Cardiac-Elastin-Shi：Shi 给 ECM 结构（elastin/collagen 组合），本文给 SMC 如何感知 ECM 刚度——二者构成"ECM 结构 → 力学环境 → SMC 感知"的链条