Skip to content

Effects of the Haemodynamic Stimulus on the Location of Carotid Plaques Based on a Patient-Specific Mechanobiological Plaque Atheroma Formation Model

Hernández-López, Patricia; Cilla, Myriam; Martínez, Miguel et al. · 2021 · Frontiers in Bioengineering and Biotechnology

Metadata

Authors: Hernández-López, Patricia; Cilla, Myriam; Martínez, Miguel; Peña, Estefanía

DOI: 10.3389/fbioe.2021.690685

Tags: #atherosclerosis #finite-element #continuum-mechanics #reaction-diffusion

摘要概述

本文针对颈动脉分叉处 atheroma 斑块的空间定位问题，提出并验证了一个由 haemodynamic stimulus 驱动的 patient-specific mechanobiological 斑块生长模型。作者将 Cilla et al. (2014) 的轴对称模型扩展到三维，耦合 Navier-Stokes 方程（管腔瞬态血流）、Darcy 定律与 Kedem-Katchalsky 三孔模型（跨内皮物质输运）以及 convection-diffusion-reaction 方程（LDL、monocytes、macrophages、foam cells、cytokines、CSMC、SSMC、collagen 等九类物质的质量守恒），并以内皮细胞 shape index (SI) 作为血流动力学刺激与通透性之间的桥梁。研究在 8 个由临床 CT/MRI 分割得到的 carotid bifurcation 几何（4 名男性 atherosclerosis 患者 + 1 名健康志愿者）上模拟 30 年斑块发展，比较 TAWSS、OSI 及一个新提出的 TAWSS-OSI 组合刺激对斑块位置与 stenosis ratio 的预测能力。核心结论：新刺激在多数几何中优于单一指标，14 个真实斑块中正确预测 11 个，stenosis ratio 最小误差 2.77%、最大误差 32.89%；但在少数几何中血流动力学无法预测斑块位置，提示其他生物学或遗传因素可能更为主导。适用人群限定于高胆固醇血症假设下的 carotid bifurcation 场景，最大限制在于无 fluid-structure interaction、无 patient-specific 生化参数且缺少斑块前真实几何。

临床问题与研究设计

  • PICO/PECO：P = 存在或不存在 carotid atheroma 斑块的患者（4 名男性 atherosclerosis 患者 + 1 名健康志愿者，§ 2.1, PDF p. 2-3）；I/E = 三种 haemodynamic stimulus（TAWSS、OSI、新组合刺激）作为 mechanobiological 模型的输入并驱动 30 年斑块生长（§ 2.6, PDF p. 9）；C = 三种刺激相互比较，并与真实临床影像中分割出的斑块位置和 stenosis ratio 对照（§ 3, PDF p. 11）；O = 斑块位置一致性（14 个真实斑块的命中数）与 stenosis ratio 误差（§ 2.7, Eq. 44, PDF p. 10）。无治疗干预，研究为 computational patient-specific modelling，不适用经典治疗对照。
  • 研究类型：other（patient-specific computational modelling study with clinical imaging validation；非 RCT/cohort/case-control/diagnostic/meta-analysis）。
  • 注册与方案：未报告。临床影像由 Hospital Clínico Universitario in Zaragoza, Spain 依据医院伦理指南提供（§ 2.1, PDF p. 2），但未给出伦理批号、注册号或预设分析方案。
  • 中心、时间与随访：单中心影像来源（Hospital Clínico Universitario in Zaragoza）；随访为计算模拟的 30 年虚拟时长（§ 2.1, PDF p. 3, 文献 Insull, 2009）；真实患者的影像采集时间窗与随访未报告。

研究人群与干预或暴露

  • 纳入与排除标准：纳入标准未以正式 clinical-trial 形式陈述；作者从 4 名男性 atherosclerosis 患者与 1 名健康志愿者的临床影像中分割出 8 个 carotid bifurcation 几何（左/右侧），编码 A–H（§ 2.1, Fig. 1, PDF p. 2-3）。性别、年龄、合并症、用药状态、血压与 LDL 水平等的逐患者基线特征未报告。排除标准未报告。
  • 样本量、基线特征与代表性：共 8 个几何（7 个病理几何 A–G + 1 个健康几何 H）；样本量小且为回顾性影像便利样本，无 a priori 样本量计算；所有病理患者统一假设高胆固醇血症 LDL = 6.98 mol/m³（≈270 mg/dL，§ 2.1, PDF p. 3, 文献 Goldstein and Brown, 1977）。壁厚固定为 CCA 0.7 mm、ICA/ECA 0.53 mm（§ 2.1, PDF p. 3, 文献 Sommer et al., 2010）。代表性有限，且仅男性，未含女性。
  • 干预、对照或暴露定义：暴露为三种 mechanical stimulus 计算的 SI 场，分别驱动 endothelial leaky junction 比例与 monocyte 募集（§ 2.6, Eq. 32-38, PDF p. 9-10）：(1) TAWSS（Eq. 30-32，阈值 2 Pa，文献 Zhao et al., 2002; Younis et al., 2004; Filipovic et al., 2013）；(2) OSI（Eq. 33-34，阈值 0.1910，由 TAWSS=2 Pa 对应 SI 反推）；(3) 新组合刺激 SI = k8·e^(k9·OSI) + k10·e^(k11·(TAWSS)²)（Eq. 35, Fig. 5, PDF p. 9，拟合参数见 Table 4）。对照为真实临床影像中分割的斑块（Fig. 1 与 Fig. 8 最后一列）。
  • 失访、交叉与依从性：不适用（计算模拟，无患者随访）；真实患者无 LDL 水平、血压、用药依从性数据（§ 4, PDF p. 15，作者明确列为局限）。

终点与统计方法

  • Primary endpoint：斑块位置一致性，定义为模型预测的斑块高 SI / 生长区域与真实影像分割斑块在几何上的吻合度（§ 3, Fig. 7-8, PDF p. 11-14）。14 个真实斑块中正确命中的数量为汇总指标（Abstract；§ 4, PDF p. 15）。
  • Secondary/safety endpoints：stenosis ratio (SR, %) 误差，SR = (1 − lumen area with plaque / lumen area without plaque) × 100（Eq. 44, § 2.7, PDF p. 10）；按几何与动脉段（CCA/ICA/ECA）分层报告于 Table 6（PDF p. 15）。safety endpoint 不适用（计算研究）。
  • 效应量与模型：效应量以 SR 绝对误差（%）表达；新刺激最大误差 32.89%、最小误差 2.77%；TAWSS 最大 45.74%、最小 3.58%（高估）；OSI 所有案例误差 >9.53%（低估）（§ 3, Table 6, PDF p. 11, 15）。无 HR/RR/OR/AUC，无置信区间与 p 值，因无重复抽样统计推断框架。
  • 缺失数据、多重比较、亚组与敏感性分析：多重比较与 formal 统计检验未报告。敏感性分析以参数变异讨论形式呈现：作者指出 LDL 相关参数对斑块生长影响最大（因 LDL 启动整个炎症过程），并用 patient A 校准后外推到其余几何（§ 2.7, PDF p. 11）；mesh sensitivity analysis 对管腔与动脉壁分别完成（§ 2.2, PDF p. 4）；未报告针对血压、LDL 水平或壁厚的系统性 uncertainty quantification。

核心结果与证据

主要发现 1：新组合刺激对斑块位置的总体预测优于 TAWSS 与 OSI

  • 临床结论：在 8 个 patient-specific carotid 几何中，新提出的 TAWSS-OSI 组合刺激在多数几何中比单独使用 TAWSS 或 OSI 更准确地预测斑块位置；14 个真实斑块中模型共正确预测 11 个。TAWSS 倾向于在 CCA 远端预测出非生理性过度生长，OSI 则系统性低估斑块面积与生长。
  • 证据定位：§ 3, Fig. 8, Table 6, PDF p. 11, 14-15；Abstract（"The model predicts correctly eleven of the fourteen plaques"）。
  • 效应量：位置命中 11/14（约 78.6%）。SR 误差：新刺激最小 2.77%、最大 32.89%；TAWSS 最小 3.58%、最大 45.74%（高估）；OSI 所有案例误差 >9.53%（低估）。逐几何 SR 数值见 Table 6（如 patient C CCA 真实 61.14%，新刺激预测 52.65%，差 8.49%；ICA 真实 32.91%，新刺激 22.71%，差 10.2%）。
  • 临床意义：新刺激在位置预测层面具有明确增量价值，提示将 TAWSS 与 OSI 联合作为 endothelial SI 的驱动变量比单一指标更贴近临床影像观察；但 SR 误差仍可达 32.89%，个体层面预测精度尚不足以直接用于单患者治疗决策。
  • 不确定性：无置信区间，因单次确定性模拟；参数源自不同物种与血管（非均为人 carotid），且 patient A 校准后外推引入系统性偏倚；3 个未命中斑块（patient E 全部、patient G 的 ICA）提示模型在这些位点失效。

主要发现 2：健康几何验证模型稳定性——无 atheroprone 刺激时斑块不生长

  • 临床结论：在健康志愿者几何（H）上，模型未产生显著斑块生长，OSI 刺激下生长可忽略，验证了模型在无 atheroprone haemodynamic 环境下的稳定性与"无假阳性生长"特性。
  • 证据定位：§ 3, Fig. 6, PDF p. 11-12；Abstract（"a healthy geometry has been computed to check that plaque is not developed with the model in this case"）。
  • 效应量：未报告数值化生长量；以定性"negligible growth"描述（OSI 刺激下）。
  • 临床意义：模型满足"健康场景不产生斑块"的最低必要条件，支持其作为风险定位工具的概念合理性；但单一健康样本不足以证明模型在更大健康人群中的特异性。
  • 不确定性：仅 1 个健康几何，无量化阈值判定"negligible"；未报告健康几何中 SI 场与病理几何的定量对比。

主要发现 3：部分斑块无法由 haemodynamics 预测，提示生物学/遗传因素主导

  • 临床结论：在 patient E 中，三种刺激均无法预测斑块位置；patient E 与 G 的 ICA 斑块亦未被任何刺激预测。作者推断这些位点的 haemodynamic disturbance 不显著，疾病触发可能源于其他系统性或遗传因素（如 external lesion、intimal layer pathological weakness）。
  • 证据定位：§ 3, Fig. 7-8, PDF p. 11-14；§ 4 Discussion, PDF p. 13（"for patients E and G, the stenosis appears in a zone where no haemodynamical disturbance is observed"）。
  • 效应量：未命中斑块数 3/14（约 21.4%）。Table 6 显示 patient E CCA 真实 SR 19.61%，新刺激预测 15.84%；patient G ICA 真实 45.81%，新刺激预测 12.91%。
  • 临床意义：这是对模型适用边界的关键负结果——它界定了 haemodynamic-driven mechanobiological 模型的解释上限，提示在约 1/5 斑块中需引入非血流动力学机制（遗传、局部血管壁结构异常）。
  • 不确定性：作者明确为假设性推断（"could be due to other systemic or genetic conditions"），未提供患者层面的遗传或组织学证据；未量化"haemodynamical disturbance 不显著"的判据。

主要发现 4：ICA 斑块预测失败可能与 haemodynamics-growth 解耦有关——更新几何后改善

  • 临床结论：作者假设 ICA 斑块预测失败源于模型中血流与斑块生长的 uncoupled 处理（血流仅在初始计算，不随斑块生长更新）。对 patient B 进行 15 年生长后更新几何再模拟 15 年，ICA 斑块预测与真实影像吻合度改善，支持 fluid-structure interaction (FSI) 耦合的必要性。
  • 证据定位：§ 4 Discussion, Fig. 9, PDF p. 13（"we computed the growth of B geometry until 15 years, and then we actualized the geometry"）。
  • 效应量：定性改善（"the plaque in the ICA better matches the real plaque"）；未给出更新前后 SR 数值对比。作者说明完整 3D FSI 因计算成本过高未实现。
  • 临床意义：指明了模型改进方向——coupled haemodynamics-growth 框架可能显著提升 ICA 等下游斑块的预测能力，对未来 patient-specific 模型设计有直接方法学价值。
  • 不确定性：仅 1 个 patient (B) 的初步验证；未量化改善幅度；未与完整 FSI 参考解对比。

安全性、负结果与亚组

  • 安全性：不适用（computational modelling study，无患者安全性终点）。
  • 负结果：(1) 14 个真实斑块中 3 个未被任何 haemodynamic stimulus 预测（patient E 全部、patient G 的 ICA），见 § 3 与 Fig. 7-8；(2) OSI 在所有案例中低估 SR（误差 >9.53%），TAWSS 在 CCA 预测出非生理性斑块（§ 3, Table 6, PDF p. 11, 15）；(3) 新刺激虽总体最优，但在 patient A 的 ICA（真实 70.45%，预测 25.30%）、patient G 的 ICA（真实 45.81%，预测 12.91%）等位点仍存在大幅低估。
  • 亚组：按几何（A–H）与动脉段（CCA/ICA/ECA）分层的 SR 结果见 Table 6；未按性别（全为男性）、年龄、LDL 水平、血压或用药状态进行亚组分析，相关患者层面协变量未报告。

关键图表

  • 图表定位：Fig. 1（真实临床影像几何 A–H，斑块以箭头标示，PDF p. 3）；Fig. 2（计算流程，PDF p. 4）；Fig. 5（SI-TAWSS-OSI 拟合曲面，PDF p. 10）；Fig. 6（健康几何 SI 与斑块生长，PDF p. 12）；Fig. 7（病理几何 A–G 的 SI 分布，PDF p. 13）；Fig. 8（病理几何斑块生长预测 vs. 真实，PDF p. 14）；Fig. 9（patient B 的 15 年更新几何实验，PDF p. 15）；Table 6（逐几何 SR 对比，PDF p. 15）。
  • 核心数字：Table 6 是定量核心——新刺激 SR 误差范围 2.77–32.89%，TAWSS 3.58–45.74%（高估），OSI 全部 >9.53%（低估）；patient C CCA 真实 61.14% vs. 新刺激 52.65%；patient A ICA 真实 70.45% vs. 新刺激 25.30%。
  • 阅读陷阱：(1) Fig. 7 仅展示 SI 分布，非斑块本身，需与 Fig. 8 生长结果对照阅读；(2) Table 6 的 "Location" 列按动脉段（CCA/ICA）分行，同一几何可有多行，阅读时需成对对照真实值；(3) Fig. 9 的"更新几何"非完整 FSI，仅为半程重启的近似，易被误读为已实现耦合；(4) SR 定义基于面积而非直径，与临床常用 NASCET 直径法不直接可比。

偏倚、局限与外部有效性

  • 选择、测量、混杂和报告偏倚：(1) 选择偏倚——仅 4 名男性患者 + 1 名健康志愿者，无女性、无样本量计算、无正式纳入/排除标准；(2) 测量偏倚——真实斑块由 thresholding segmentation 获取，小侧支被人为剔除以简化计算（§ 2.1, PDF p. 3），且作者无真实斑块前几何，只能"数字移除"斑块反推健康几何，引入重构误差；(3) 混杂偏倚——所有患者统一假设 LDL=270 mg/dL 高胆固醇血症，但真实患者的 LDL、血压、用药、病程均未知（§ 4, PDF p. 15）；(4) 报告偏倚——参数多源自不同物种与血管的文献（非均为人 carotid），且 patient A 校准后外推至其余几何，存在 overfitting 风险，作者未提供 hold-out 验证。
  • 因果解释：部分成立。模型展示了 haemodynamic stimulus → endothelial SI → leaky junction → LDL/monocyte 跨壁输运 → 炎症级联 → 斑块生长的机制链，但为 computational association 而非因果干预证明；patient E/G 的失败案例本身反证了纯 haemodynamic 因果链的不完备性。
  • 适用患者：高胆固醇血症假设下的男性 carotid bifurcation atherosclerosis 患者，用于斑块位置风险定位的计算辅助。
  • 不适用患者：女性、非 carotid 血管床、非高胆固醇血症表型、已存在复杂斑块形态或遗传性血管病变的患者；亦不适用于需要个体化治疗决策或急性事件预测的临床场景。
  • 与指南或标准治疗的关系：本次未核实（文章未引用或对比任何现行 atherosclerosis 临床指南，如 ESC 或 AHA/ACC 指南）。

数据、代码与可复现性

  • 数据：Data Availability Statement 声明"raw data supporting the conclusions of this article will be made available by the authors, without undue reservation"（§ Data Availability, PDF p. 16），但未提供公开仓库链接或 DOI；几何、参数与影像数据本次未核实是否实际可获取。
  • 代码：未提供。文章使用商业软件 COMSOL Multiphysics（COMSOL AB, Burlington, MA, USA）与 Materialise Mimics、Rhinoceros，未报告自定义脚本仓库或版本号（§ 2.1-2.2, PDF p. 3-4）。Supplementary Material 链接为 Frontiers 文章页（https://www.frontiersin.org/articles/10.3389/fbioe.2021.690685/full#supplementary-material），本次未核实其内容。
  • 复现判断：部分可复现。方程、边界条件、参数表（Table 1-5）与计算流程（Fig. 2, 四步法：transient flow → steady plasma → transient substances 30 年 → steady growth）在正文中披露较完整，理论上可在 COMSOL 中重建；但缺少公开代码、原始几何文件与 patient-specific 参数，且参数校准依赖 patient A 的未公开数据，独立复现仍受阻。

个人批注

对当前课题的意义

本文是将 WSS 衍生指标（TAWSS/OSI）与组织级 convection-diffusion-reaction 斑块生长模型耦合的代表性 patient-specific 建模范例，对我理解"机械刺激如何驱动血管壁重塑与斑块定位"的建模框架有直接参考价值。其将 endothelial shape index 作为血流动力学与物质通透性之间的桥梁变量，提供了一条可操作的耦合路径；同时 patient E/G 的负结果与 Fig. 9 的更新几何实验，清晰界定了 uncoupled haemodynamics-growth 框架的失效边界，对我当前关注的患者特异性建模误差来源与耦合策略选择有警示意义。

可复用的设计或统计方法

(1) 四步法计算流程（transient Navier-Stokes → steady Darcy/Kedem-Katchalsky → transient CDR 30 年 → steady growth）可作为多尺度血管建模的工作流模板；(2) 三孔模型（normal/leaky junction/vesicular）与 leaky junction 比例由 SI 驱动的参数化方式可直接复用；(3) SI 作为 TAWSS 与 OSI 的组合函数（Eq. 35, 指数和形式）提供了一种将多指标融合为单一通透性驱动变量的简洁参数化；(4) 以 stenosis ratio（Eq. 44, 面积法）作为可验证终点的量化策略；(5) 健康几何作为阴性对照验证模型特异性的做法值得借鉴。

疑问与后续动作

  • 疑问 1：新组合刺激 SI = k8·e^(k9·OSI) + k10·e^(k11·(TAWSS)²) 的拟合参数 k8–k11（Table 4）依赖 Sáez et al. (2015) 的 pseudo-experimental 数据，该拟合在外推到 carotid 之外的血管床时是否仍稳定？
  • 疑问 2：patient A 的参数校准是否构成 overfitting——若换用 patient B 或 C 校准再外推至 A，命中率与 SR 误差会如何变化？
  • 疑问 3：Fig. 9 的 15 年更新几何实验能否推广为一种"分段耦合"的廉价替代，以规避完整 3D FSI 的计算成本？
  • 后续动作：(1) 检索本文引用 Cilla et al. (2014) 的后续论文，确认是否有作者团队或第三方实现完整 FSI 耦合的版本；(2) 在本地论文队列中比对 Olgac et al. (2011) 的三孔模型实现细节，核对参数一致性；(3) 考虑在新建模工作中复用 Eq. 35 的 SI 组合形式并做敏感性分析。

与上下文的关系

本文建立在

本文是 Cilla et al. (2014) "Mathematical modelling of atheroma plaque formation and development in coronary arteries"（J. R. Soc. Interface 11:20130866）的直接三维化与患者特异化扩展，作者明确声明"evolution of the one previously proposed by the authors"（Abstract, § 4, PDF p. 11）。SI-TAWSS 关系基于 Levesque et al. (1986) 的狗主动脉实验数据（Eq. 32, Fig. 4A）；SI-OSI 与组合刺激基于 Sáez et al. (2015) 的理论模型（Eq. 34-35, Fig. 4B, Fig. 5）；三孔模型与跨内皮 LDL 输运基于 Olgac et al. (2008) 与 Weinbaum et al. (1985)（Eq. 5-10, Appendix A-B）；斑块生长的运动学基于 Garikipati et al. (2004) 的 continuum growth 框架（Eq. 39-41）。血流边界条件来自 Malvè et al. (2014)。

已核实的后续引用

本次未检索。

同类临床证据对比

本文与本地队列中多篇 WSS/OSI 驱动的 atherosclerosis 建模与综述工作可比较：(1) 2013-Does-Low-Oscillatory-Peiffer 系统评估了低 TAWSS 与高 OSI 作为斑块定位指标的因果证据强度，本文的"OSI 阈值 0.1910"与"TAWSS 阈值 2 Pa"设定可与之对照；(2) 2011-Computed-High-Concentrations-Olgac 提供了本文直接引用的三孔模型与 LDL 跨壁输运计算框架，本文的 Eq. 5-10 与 Appendix A-B 与其同源；(3) 2013-Shear-Stress-Plaque-Wentzel 从临床角度综述 WSS 与斑块进展的关系，可作为本文机制预测的临床映射参照；(4) 2016-BiomechForces-Athero-Brown 与 2005-Role-Shear-Stress-Slager 提供了 biomechanical forces 驱动 atherogenesis 的更广生物学背景；(5) 2021-Plaque-Modeling-Cai 综述了斑块建模三代演进，本文属于其所述"第三代"（微环境组分演化耦合到形态变化）的实例；(6) 2017-VascularAdaptation-Garbey 提供了 agent-based model 替代路径，作者在 § 4 明确对比了 continuum model（本文）与 ABM 在计算成本与细胞随机行为建模上的互补优劣。本文相对于上述工作的差异在于：首次在 patient-specific 3D carotid 上将 TAWSS+OSI 组合刺激耦合到完整的九物质 CDR 生长模型，并报告逐几何 SR 误差。

与本地论文队列的关系

快速判断

  • 一句话结论：本文提出一个由血流动力学刺激驱动的 mechanobiological 动脉粥样硬化斑块生长模型，并在 8 个基于临床影像分割的 patient-specific 颈动脉分叉几何上模拟了 30 年的斑块发展过程。
  • 阅读范围：仅 metadata、first page、abstract、最多 3 个 figure/table captions 和 conclusion；未通读正文
  • 处理决定：保留参考
  • 决定理由：该文与模型、计算框架与血管或细胞系统的机制研究相关；受限材料足以保留研究线索，但不足以支持全文级方法或稳健性判断。

摘要概述

本文提出一个由血流动力学刺激驱动的 mechanobiological 动脉粥样硬化斑块生长模型，并在 8 个基于临床影像分割的 patient-specific 颈动脉分叉几何上模拟了 30 年的斑块发展过程。模型耦合了 Navier-Stokes 方程（管腔内瞬态血流）、Darcy 定律与 Kedem-Katchalsky 方程（内皮物质跨膜输运，三孔模型）以及 convection-diffusion-reaction 方程（LDL、foam cell 等物质的质量守恒），并通过内皮细胞形状指数调控物质通透性。作者比较了 TAWSS、OSI 及一个新提出的组合刺激（基于 TAWSS 与 OSI 拟合内皮细胞 shape index）对斑块位置与尺寸的预测能力，发现新刺激在多数几何中优于单一指标，14 个真实斑块中正确预测 11 个，狭窄度误差最小 2.77%、最大 32.89%。但在少数几何中血流动力学无法预测斑块位置，提示生物学或遗传因素可能更为主导。

关键证据

  • 证据 1（定位）：Conclusion
  • 展示或报告：但在少数几何中血流动力学无法预测斑块位置，提示生物学或遗传因素可能更为主导。
  • 支持的结论：但在少数几何中血流动力学无法预测斑块位置，提示生物学或遗传因素可能更为主导。
  • 注意事项：仅据受限 quick source；方法细节、完整定量结果与稳健性分析本次未核实。

局限与未核实项

  • 最大局限：当前仅完成 quick triage，不能据此确认正文中的全部实验、模型或统计假设。
  • 未核实项：完整方法、样本或参数设置、敏感性分析、局限讨论及数据与代码可用性。

与我的关联

  • 可复用点：文中的模型或计算框架可作为相关建模工作的候选参考；具体方程、参数与适用条件本次未核实。
  • 关联的当前问题：该工作是将血流动力学 WSS 衍生指标与基于反应-扩散-对流方程的组织级斑块生长模型相耦合的代表性 patient-specific 建模范例，对我理解机械刺激如何驱动血管壁重塑与斑块定位的建模框架有直接参考价值。
  • 下一步动作：保留参考；仅在当前问题需要其完整方法或定量证据时升级为 deep note。