[论文解读] Particle-based multiscale modeling of intracellular calcium dynamics

Ulrich Dobramysl, Sten Rüdiger+1 more|arXiv (Cornell University)|Apr 1, 2015

Gene Regulatory Network Analysis参考文献 24被引用 6

一句话总结

本文提出了一种基于粒子的多尺度模型，将随机钙通道门控与单个Ca²⁺离子的布朗运动相结合，用于模拟细胞内钙脉冲。研究发现，只有当通道抑制的时间尺度足够长时，钙脉冲才会出现，通过平均场分析和离子追踪确定了脉冲生成的关键区域。

ABSTRACT

Intracellular calcium is regulated by the release of Ca$^{2+}$ ions from the endoplasmic reticulum via inisitol-4,5-triphosphate receptor channels. The resulting dynamics are highly diverse, lead to local calcium "puffs" as well as global waves propagating through cells. Local fluctuations in the number of calcium ions play a crucial role in the onset of these features. Previous modeling studies have focused on stochastic channel models coupled to deterministic diffusion of ions, thereby neglecting local fluctuations of the ion number. Tracking of individual ions is computationally difficult due to the scale separation in the Ca$^{2+}$ concentration when channels are in the open or closed states. In this paper, a spatial multiscale model for investigating of the dynamics of puffs is presented. It couples Brownian motion (diffusion) of ions with a stochastic channel gating model. The model is used to analyze calcium puff statistics. Concentration time traces as well as channel state information are studied. We identify the regime in which puffs can be found and develop a mean-field theory to extract the boundary of this regime. Puffs are only possible when the time scale of channel inhibition is sufficiently large. Implications for the understanding of puff generation and termination are discussed.

研究动机与目标

• 解决先前模型忽略钙信号中局部离子数波动的局限性。
• 开发一种空间分辨的基于粒子的模型，以同时捕捉随机通道行为和离子扩散。
• 分析钙脉冲的统计特性，包括浓度时间轨迹和通道状态动力学。
• 确定钙脉冲能够出现和终止的条件。
• 推导平均场理论，以预测脉冲生成区域的边界。

提出的方法

• 该模型将单个Ca²⁺离子的布朗运动与肌醇-1,4,5-三磷酸受体（IP3R）通道的随机门控模型相耦合。
• 使用适当的阻尼和噪声项，通过朗之万动力学模拟离子扩散。
• 基于实验数据导出的转换速率，随机更新通道状态（开放、关闭、失活）。
• 该模型同时追踪离子浓度随时间的变化以及通道状态的时序演化。
• 推导出平均场近似，以解析估计脉冲生成区域的边界。
• 模拟生成浓度轨迹和通道状态序列，用于脉冲特征的统计分析。

实验结果

研究问题

• RQ1 在存在随机离子数波动和通道门控的情况下，钙脉冲在何种条件下会出现？
• RQ2 通道抑制的时间尺度如何影响钙脉冲的稳定性和持续时间？
• RQ3 离子扩散、通道开放动力学与局部钙脉冲形成之间存在何种关系？
• RQ4 平均场理论如何预测参数空间中脉冲生成区域的边界？
• RQ5 局部离子数波动在脉冲启动和终止过程中起什么作用？

主要发现

• 只有当通道抑制的时间尺度足够长时，脉冲才可能发生，表明脉冲形成存在关键的动力学要求。
• 通过整合单个离子动力学与随机通道门控，该模型成功捕捉了真实的脉冲统计特性。
• 浓度时间轨迹和通道状态序列表明，脉冲的起始和终止由离子扩散与通道失活动力学之间的相互作用所决定。
• 平均场理论成功预测了脉冲生成区域的边界，验证了模拟结果。
• 局部离子数波动显著影响脉冲动力学，而这类效应在确定性扩散模型中被忽略。
• 该模型揭示，脉冲生成不仅依赖于通道开放，还依赖于离子释放与通道失活之间的时序相关性。