Department of Mathematics and Department of Biology, Duke University, Durham, USA.
Department of Mathematics, Tulane University, New Orleans, USA.
Bull Math Biol. 2021 Jan 16;83(3):21. doi: 10.1007/s11538-020-00847-3.
In developmental biology as well as in other biological systems, emerging structure and organization can be captured using time-series data of protein locations. In analyzing this time-dependent data, it is a common challenge not only to determine whether topological features emerge, but also to identify the timing of their formation. For instance, in most cells, actin filaments interact with myosin motor proteins and organize into polymer networks and higher-order structures. Ring channels are examples of such structures that maintain constant diameters over time and play key roles in processes such as cell division, development, and wound healing. Given the limitations in studying interactions of actin with myosin in vivo, we generate time-series data of protein polymer interactions in cells using complex agent-based models. Since the data has a filamentous structure, we propose sampling along the actin filaments and analyzing the topological structure of the resulting point cloud at each time. Building on existing tools from persistent homology, we develop a topological data analysis (TDA) method that assesses effective ring generation in this dynamic data. This method connects topological features through time in a path that corresponds to emergence of organization in the data. In this work, we also propose methods for assessing whether the topological features of interest are significant and thus whether they contribute to the formation of an emerging hole (ring channel) in the simulated protein interactions. In particular, we use the MEDYAN simulation platform to show that this technique can distinguish between the actin cytoskeleton organization resulting from distinct motor protein binding parameters.
在发育生物学以及其他生物系统中,可以使用蛋白质位置的时间序列数据来捕捉新兴的结构和组织。在分析这种时变数据时,不仅要确定拓扑特征是否出现,还要确定它们形成的时间,这是一个常见的挑战。例如,在大多数细胞中,肌动蛋白丝与肌球蛋白马达蛋白相互作用,并组织成聚合物网络和更高阶结构。环道就是这种结构的一个例子,它保持直径不变,在细胞分裂、发育和伤口愈合等过程中发挥关键作用。鉴于在体内研究肌动蛋白与肌球蛋白相互作用的局限性,我们使用复杂的基于代理的模型在细胞中生成蛋白质聚合物相互作用的时间序列数据。由于数据具有丝状结构,我们建议沿着肌动蛋白丝进行采样,并在每个时间点分析所得点云的拓扑结构。基于持久同调中的现有工具,我们开发了一种拓扑数据分析(TDA)方法,用于评估动态数据中有效环的生成。该方法通过与数据中组织出现相对应的路径,将拓扑特征与时序相关联。在这项工作中,我们还提出了评估感兴趣的拓扑特征是否显著的方法,从而评估它们是否有助于模拟蛋白质相互作用中新兴孔(环道)的形成。特别是,我们使用 MEDYAN 模拟平台来证明该技术可以区分不同马达蛋白结合参数导致的肌动蛋白细胞骨架组织。