机械测试、在体光学相干弹性成像与个性化有限元建模的整合用于预测角膜交联术的几何学 outcomes

《Annals of Biomedical Engineering》：Integration of Mechanical Testing, In Vivo Optical Coherence Elastography and Personalized Finite Element Modeling to Predict Geometrical Outcomes of Corneal Cross-Linking

　　

目的：角膜交联术（CXL）诱导角膜产生机械学与几何学改变，这些改变通常在术前规划中被忽视。研究人员开发并校准了一个患者特异性有限元模型（FEM），用于预测CXL引起的地形图改变，并将其应用于三例圆锥角膜（KC）患者作为概念验证。方法：为校准该模型，研究人员对五枚人供体角膜在CXL前后进行了纳米压痕（NI）及离体光学相干弹性成像（OCE）膨胀测试。纳米压痕结果用于调整黏超弹性参数，而离体OCE轴向应变用于验证。个性化角膜模型基于三例KC患者的地形图构建，受影响区域的区域刚度根据在体压力调制OCE测量的轴向应变进行调整。随后将模拟CXL结果与6个月临床结果进行比较。结果：CXL使纤维相关机械参数增加16倍，并将黏弹性时间常数降低为原来的三分之一。在体OCE测量显示KC区域机械学弱化平均达57%。与6个月随访临床地形图相比，模型预测的CXL诱导曲率改变分别为：三例患者分别为?1.5 D vs. ?1.76 D、?1.65 D vs. ?1.91 D、以及?1.76 D vs. ?1.57 D。结论：通过将FEM与在体角膜机械学特性表征相结合，可以预测患者特异性地形图改变，这有助于改善CXL治疗的规划。

该研究发表于《Annals of Biomedical Engineering》，聚焦于圆锥角膜（KC）这一进行性眼病的治疗优化问题。KC以角膜变薄、前凸呈锥形为特征，伴随局部生物力学退化，表现为胶原纤维排列紊乱及弹性模量降低。尽管角膜交联术（CXL）通过紫外线A照射联合核黄素光敏剂建立交联以加固角膜、延缓KC进展，但CXL诱导的屈光改变尚未被纳入术前规划，其疗效与机制仍有待阐明。现有有限元模型（FEM）多忽略组织黏弹性，且缺乏独立数据集校准及临床验证。为此，研究人员整合机械测试与逆有限元建模（iFEM），建立可预测CXL地形图改变的个性化模型，并在三例临床病例中验证其适用性。

研究采用的主要关键技术方法包括：纳米压痕（NI）测试获取黏超弹性参数，离体光学相干弹性成像（OCE）膨胀测试验证模型，以及基于Python的自动化网格生成与FEBio有限元求解。临床样本方面，五枚人供体角膜（75±5.4岁）由LIONS Cornea Bank Baden-Württemberg at the Eye Center, University of Freiburg提供；三例KC患者数据集来自Canton of Zurich伦理委员会批准的研究（2021-02275）。在体OCE采用Heidelberg Engineering公司ANTERION设备进行压力调制测量，个性化模型依据患者术前Pentacam地形图构建，CXL模拟参照ELZA-PACE定制方案。

研究结果部分包含以下方面：

"Ex Vivo Experiments"：CCT在CXL前后分别为419±60 μm和484±33 μm（p=0.125）。NI显示CXL区域E_HZ从99.96±15.8 kPa升至141.0±4.7 kPa（p<0.01），C_IT从19.02±2.0%降至16.48±1.3%（p<0.01）。离体OCE显示CXL后前300 μm应变幅度降低（?5.92±0.8‰ vs 1.1±2.6‰，p=0.06）。

"Mechanical Parameter Identification"：iFEM优化后，CXL使纤维相关参数k₁增强16倍，g₁增加175%，τ₁降低288%。OCE膨胀测试模拟中，FEM轴向应变（?5.94‰ vs ?1.0‰）落入实验数据范围内，验证了模型准确性。

"Clinical Validation of the Model"：在体OCE显示三例患者KC区域相对其余角膜变形高2.15?2.72倍，平均机械弱化57.1%。k₁据此分别降低至原值的45.4%、46.5%和36.7%。模型预测与6个月临床随访的ΔKmaxZonalMean3高度吻合（?1.5 D vs ?1.76 D、?1.65 D vs ?1.91 D、?1.76 D vs ?1.57 D），尤其在3 mm光学区内一致性最佳。

讨论部分指出，该FEM同时考量了组织的压缩与拉伸行为、黏弹性特征及深度依赖性刚度分布，并首次整合在体OCE数据实现个性化机械属性赋值。与既往研究比较：Nohava等报告的弹性模量增幅较低可能源于固定方式差异；Roy和Dupps报道的200?300%杨氏模量增加对应更大曲率改变，与本研究?1.75 D平均降幅相符；Wollensak等的350%增加源于更高应力水平下的测试。研究也承认局限性：iFEM基于平均几何模型且未纳入零应力几何；供体角膜年龄偏大、无KC病变，可能影响参数外推；实验用自制遮光罩可能存在散射光污染；仅三例患者属概念验证，缺乏统计效力与IOP个体化数据；未获取术后在体OCE数据以评估生物重塑效应。尽管存在上述局限，该框架完全基于开源软件，单次运行约三小时，为临床转化奠定基础。未来可结合物理信息神经网络加速计算，并整合偏振敏感OCT获取个体化纤维分布。

结论：通过整合压缩与膨胀测试以全面表征角膜的时间依赖性与非线性力学行为，研究人员开发了新型患者特异性CXL模型。该模型依据在体测量的几何与力学属性进行得最优化制，其预测结果与CXL术后六个月地形图矫正高度吻合，表明所提出的有限元模型能够改善定制化CXL手术的术前规划。