Department of Mathematics and Biology, Duke University, Durham, USA.
School of Mathematical and Statistical Sciences, Arizona State University, Tempe, USA.
Bull Math Biol. 2020 Sep 16;82(10):126. doi: 10.1007/s11538-020-00797-w.
In many biological systems, the movement of individual agents is characterized having multiple qualitatively distinct behaviors that arise from a variety of biophysical states. For example, in cells the movement of vesicles, organelles, and other intracellular cargo is affected by their binding to and unbinding from cytoskeletal filaments such as microtubules through molecular motor proteins. A typical goal of theoretical or numerical analysis of models of such systems is to investigate effective transport properties and their dependence on model parameters. While the effective velocity of particles undergoing switching diffusion dynamics is often easily characterized in terms of the long-time fraction of time that particles spend in each state, the calculation of the effective diffusivity is more complicated because it cannot be expressed simply in terms of a statistical average of the particle transport state at one moment of time. However, it is common that these systems are regenerative, in the sense that they can be decomposed into independent cycles marked by returns to a base state. Using decompositions of this kind, we calculate effective transport properties by computing the moments of the dynamics within each cycle and then applying renewal reward theory. This method provides a useful alternative large-time analysis to direct homogenization for linear advection-reaction-diffusion partial differential equation models. Moreover, it applies to a general class of semi-Markov processes and certain stochastic differential equations that arise in models of intracellular transport. Applications of the proposed renewal reward framework are illustrated for several case studies such as mRNA transport in developing oocytes and processive cargo movement by teams of molecular motor proteins.
在许多生物系统中,个体代理的运动表现出具有多种不同的定性行为,这些行为源于各种生物物理状态。例如,在细胞中,囊泡、细胞器和其他细胞内货物的运动受到它们与细胞骨架丝(如微管)结合和解离的影响,这是通过分子马达蛋白实现的。对这些系统的理论或数值模型进行分析的典型目标是研究有效传输特性及其对模型参数的依赖性。虽然经历切换扩散动力学的粒子的有效速度通常可以很容易地根据粒子在每个状态下花费的长时间分数来描述,但有效扩散率的计算更加复杂,因为它不能简单地用粒子在一个时刻的传输状态的统计平均值来表示。然而,这些系统通常是再生的,也就是说它们可以分解为独立的循环,这些循环以返回到基本状态为标志。通过使用这种分解,我们通过计算每个循环内动力学的矩并应用更新奖励理论来计算有效传输特性。这种方法为线性对流-反应-扩散偏微分方程模型的直接均匀化提供了一种有用的大时间分析替代方法。此外,它适用于一类一般的半马尔可夫过程和某些在细胞内运输模型中出现的随机微分方程。提出的更新奖励框架的应用在几个案例研究中得到了说明,例如发育中的卵母细胞中的 mRNA 运输和分子马达蛋白团队的连续货物运动。