从噪声测量中提取准确的细胞骨架肌动蛋白速度分布

Extraction of accurate cytoskeletal actin velocity distributions from noisy measurements.

作者信息Cayla M Miller, Elgin Korkmazhan, Alexander R Dunn
PMID35963858
发布时间2022-08-13
DOI10.1038/s41467-022-31583-y

实验完整度

高

包含体外细胞模型(HFF、HUVEC)、固定细胞校准、药物扰动(latrunculin A、blebbistatin)及温度变化,并通过贝叶斯拟合和蒙特卡洛模拟进行机制验证。

主要模型

人包皮成纤维细胞(HFF) 人脐静脉内皮细胞(HUVEC) 固定HFF细胞(用于定位误差校准)

重点核对

SiR-actin标记浓度20–80 nM,孵育9分钟 成像:TIRF显微镜,300 ms曝光,每2秒一帧,共60帧(120秒) 固定细胞定位误差拟合参数:σxy1=23 nm,σxy2=50 nm,f1=0.57 药物处理:200 nM latrunculin A或25 μM blebbistatin 温度条件:37 °C(对照)或26 °C

摘要

肌动蛋白细胞骨架的动态重塑对许多细胞过程至关重要。追踪单个肌动蛋白丝的运动原则上可以揭示这种复杂行为如何在分子水平上产生。然而,从这些测量中提取的信息往往受到低信噪比的限制。我们开发了一种贝叶斯统计方法,利用量化的定位不确定性,从单个肌动蛋白相关荧光团的轨迹中估计真实的、潜在的速度分布。我们发现,成纤维细胞和内皮细胞中丝状(F)肌动蛋白的运动用统计跳跃过程比用丝状物经历连续扩散运动的模型描述得更好。特别是,一个具有指数分布的跳跃长度和时间尺度的模型再现了细胞皮层、基于整合素的黏附和应力纤维中测量到的肌动蛋白丝速度分布,表明一个共同的物理模型可能潜在地描述多种细胞环境中的肌动蛋白丝动力学。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
如何从低信噪比的单粒子追踪数据中准确推断活细胞中F-肌动蛋白丝的真实速度分布,并区分连续扩散与不连续跳跃模型?
核心机制
F-肌动蛋白运动更符合统计跳跃过程(指数分布的跳跃距离和等待时间),而非连续扩散加漂移模型;加入方向反转概率的1D跳跃模型可用3个全局参数描述2–40秒时间尺度上的数据。
主要证据
在HFF和HUVEC中,通过TIRF显微镜和QFSM追踪SiR-actin标记的单个荧光团,利用固定细胞量化定位误差,贝叶斯推断真实位移分布;韦布尔分布拟合优于指数、高斯和瑞利分布;跳跃模型拟合优于扩散加漂移模型。
研究意义
提供了一种适用于低信噪比单粒子追踪的贝叶斯分析方法,可能用于其他粒子追踪测量和超分辨数据分析;跳跃模型表明尽管细胞骨架分子尺度复杂,其涌现物理性质可能用共同物理框架理解。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

细胞培养与荧光标记

建立可用于活细胞单分子追踪的F-肌动蛋白标记体系。

HFF和HUVEC细胞培养,使用SiR-actin以低浓度标记F-肌动蛋白,固定细胞同样处理用于定位误差校准。

2

活细胞TIRF成像与单粒子追踪

获取单个肌动蛋白相关荧光团的时空轨迹。

使用TIRF显微镜在细胞基底面成像SiR-actin puncta,每2秒一帧,共120秒;使用QFSM软件识别puncta和轨迹,去噪后亚像素定位。

3

定位误差量化

估计测量噪声,为贝叶斯推断提供先验。

在固定细胞中追踪SiR-actin,将位移分布拟合为两个非中心χ分布的混合(s=0),通过最大似然估计得到定位误差参数。

4

贝叶斯推断真实位移分布

从测量的位移分布中反卷积定位误差,得到真实位移分布。

使用贝叶斯定理,假设真实位移服从韦布尔分布、指数分布、折叠高斯分布或瑞利分布,通过最大似然拟合推断真实分布。

5

机制模型拟合与比较

检验跳跃过程模型与扩散加漂移模型对数据的描述能力。

比较扩散加漂移模型(Rice分布)与1D位置跳跃模型(指数跳跃距离和等待时间),并加入方向反转概率;使用蒙特卡洛模拟和迭代拟合。

6

扰动实验验证

检验跳跃模型参数对细胞骨架动力学扰动的响应。

在HFF中分别进行温度降低(26 °C)、latrunculin A处理(200 nM)和blebbistatin处理(25 μM),追踪F-actin并拟合跳跃模型参数。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
高糖DMEM培养基Gibco21063-029
胎牛血清Corning35-011-CV
MEM非必需氨基酸Gibco11140050
青霉素/链霉素Gibco15140122
EGM-2基础培养基LonzaCC-3156
EGM-2 MV Bullet试剂盒LonzaCC-4147
SiR-肌动蛋白Cytoskeleton Inc.CY-SC001
Prolong Live抗淬灭试剂InvitrogenP36975
latrunculin A----
blebbistatin----

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
活细胞成像与追踪
SiR-actin浓度(20–80 nM,批次间需调整以达到相似标记密度)、孵育时间9分钟、TIRF显微镜曝光300 ms、帧间隔2秒、总时长120秒
阅读提示:Methods: Experimental setup and imaging; Fig. 1图注
定位误差校准
固定细胞制备(4%多聚甲醛15分钟)、SiR-actin标记浓度、拟合模型(两个非中心χ分布混合)、拟合参数σxy1=23 nm,σxy2=50 nm,f1=0.57
阅读提示:Methods: Actin tracking analysis; Fig. 2c; Supplementary Note 1
贝叶斯推断真实分布
假设的真实分布函数形式(韦布尔分布)、时间间隔范围2–40秒、最大似然估计方法
阅读提示:Methods: Model fitting; Supplementary Note 2
跳跃模型拟合
模拟粒子数10,000、模拟时长至少40秒、方向反转概率参数、优化算法(遗传算法)
阅读提示:Methods: Model fitting; Fig. 4; Supplementary Note 3
扰动实验
latrunculin A浓度200 nM、blebbistatin浓度25 μM、温度26 °C vs 37 °C、处理时间(文中未明确说明)
阅读提示:Results: F-actin velocity distribution fits can estimate physical parameters of underlying motion; Supplementary Figs. 9–10