Simulation methods for spatiotemporal models of biochemical signaling networks Public Deposited

Downloadable Content

Download PDF
Last Modified
  • March 21, 2019
  • Strychalski, Wanda
    • Affiliation: College of Arts and Sciences, Department of Mathematics
  • Cells use signaling networks consisting of multiple interacting proteins to respond to changes in their environment. In many situations, such as chemotaxis, spatial and temporal information must be transmitted through a signaling network. Recent computational studies have emphasized the importance of cellular geometry in signal transduction, but have been limited in their ability to accurately represent complex cell morphologies. We present a finite volume method that addresses this problem. Our method uses Cartesian cut cells in a differential algebraic formulation to handle the complex boundary dynamics encountered in biological systems. The method is second order in space and time. Several models of signaling systems are simulated in realistic cell morphologies obtained from live cell images. We then examine the effects of geometry on signal transduction. External signals can trigger cells to polarize and move in a specific direction. During migration, spatially localized activity of proteins is maintained. To investigate the effects of morphological changes on intracellular signaling, we present a numerical scheme consisting of a cut cell finite volume spatial discretization coupled with level set methods to simulate the resulting advection-reaction-diffusion equation. We then show that shape deformations drive a Turing-type system into an unstable regime. The method is also applied to a model of a signaling network in a migrating fibroblast. Determining the signaling mechanisms used by membrane proteins that interact with the cytoskeleton is important for understanding phenomena such as T-cell activation and viral infection. To investigate these interactions, recent experiments have tracked the movements of single lipids and glycosyl-phosphatidylinositol (GPI) anchored protein clusters tagged with 40 nm gold particles. These experiments reveal regions of transient confinement and transient anchorage of the particles. The distribution of transient anchorage release times exhibits a long tail. We developed a stochastic model of the system to explain the transient anchorage release times and the underlying biochemical reaction system.
Date of publication
Resource type
Rights statement
  • In Copyright
  • Adalsteinsson, David
  • Elston, Timothy
Degree granting institution
  • University of North Carolina at Chapel Hill
  • Open access

This work has no parents.