Simulation methods for spatiotemporal models of biochemical signaling networks

WANDA STRYCHALSKI: Simulation Methods for Spatiotemporal Models of Biochemical Signaling Networks (Under the direction of David Adalsteinsson and Timothy C. Elston) 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

[1]  K. Jacobson,et al.  Transient anchorage of cross-linked glycosyl-phosphatidylinositol–anchored proteins depends on cholesterol, Src family kinases, caveolin, and phosphoinositides , 2006, The Journal of cell biology.

[2]  James A. Sethian,et al.  Level Set Methods and Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid , 2012 .

[3]  Baba C. Vemuri,et al.  A fast level set based algorithm for topology-independent shape modeling , 1996, Journal of Mathematical Imaging and Vision.

[4]  D. Odde,et al.  Potential for Control of Signaling Pathways via Cell Size and Shape , 2006, Current Biology.

[5]  D. Gillespie Exact Stochastic Simulation of Coupled Chemical Reactions , 1977 .

[6]  Shuji Ishihara,et al.  A Mass Conserved Reaction–Diffusion System Captures Properties of Cell Polarity , 2006, PLoS Comput. Biol..

[7]  B Bunow,et al.  Pattern formation by reaction-diffusion instabilities: application to morphogenesis in Drosophila. , 1980, Journal of theoretical biology.

[8]  P. Colella,et al.  A higher-order embedded boundary method for time-dependent simulation of hyperbolic conservation laws , 2000 .

[9]  H. Dohlman,et al.  G proteins and pheromone signaling. , 2002, Annual review of physiology.

[10]  B. Kholodenko Cell-signalling dynamics in time and space , 2006, Nature Reviews Molecular Cell Biology.

[11]  Andrew Harvey,et al.  Forecasting, Structural Time Series Models and the Kalman Filter. , 1991 .

[12]  B N Kholodenko,et al.  Diffusion control of protein phosphorylation in signal transduction pathways. , 2000, The Biochemical journal.

[13]  Giovanni Dietler,et al.  Gradient of rigidity in the lamellipodia of migrating cells revealed by atomic force microscopy. , 2005, Biophysical journal.

[14]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[15]  C. Peskin,et al.  Implicit second-order immersed boundary methods with boundary mass , 2008 .

[16]  Phillip Colella,et al.  A cartesian grid embedded boundary method for the heat equation and poisson's equation in three dimensions , 2004 .

[17]  Phillip Colella,et al.  A Cartesian grid embedded boundary method for hyperbolic conservation laws , 2006 .

[18]  James A. Sethian,et al.  The Fast Construction of Extension Velocities in Level Set Methods , 1999 .

[19]  Andrew B Goryachev,et al.  Dynamics of Cdc42 network embodies a Turing‐type mechanism of yeast cell polarity , 2008, FEBS letters.

[20]  Katherine C. Chen,et al.  Sniffers, buzzers, toggles and blinkers: dynamics of regulatory and signaling pathways in the cell. , 2003, Current opinion in cell biology.

[21]  M. Ikawa,et al.  Glycosylphosphatidylinositol-anchor-deficient mice: implications for clonal dominance of mutant cells in paroxysmal nocturnal hemoglobinuria. , 1996, Blood.

[22]  Youcef Saad,et al.  A Basic Tool Kit for Sparse Matrix Computations , 1990 .

[23]  K. Hahn,et al.  Spatiotemporal dynamics of RhoA activity in migrating cells , 2006, Nature.

[24]  David S. Stoffer,et al.  Time series analysis and its applications , 2000 .

[25]  N. Britton Essential Mathematical Biology , 2004 .

[26]  E. Aandahl,et al.  Localized effects of cAMP mediated by distinct routes of protein kinase A. , 2004, Physiological reviews.

[27]  Ravi Iyengar,et al.  Cell Shape and Negative Links in Regulatory Motifs Together Control Spatial Information Flow in Signaling Networks , 2008, Cell.

[28]  Mario Mellado,et al.  Role of the Pi3k Regulatory Subunit in the Control of Actin Organization and Cell Migration , 2000, The Journal of cell biology.

[29]  Samuel A. Isaacson,et al.  Incorporating Diffusion in Complex Geometries into Stochastic Chemical Kinetics Simulations , 2006, SIAM J. Sci. Comput..

[30]  Glazier,et al.  Simulation of the differential adhesion driven rearrangement of biological cells. , 1993, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[31]  W. Rappel,et al.  Directional sensing in eukaryotic chemotaxis: a balanced inactivation model. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[32]  B N Kholodenko,et al.  Spatial gradients of cellular phospho‐proteins , 1999, FEBS letters.

[33]  Matthias Peter,et al.  Nuclear sequestration of the exchange factor Cdc24 by Far1 regulates cell polarity during yeast mating , 2000, Nature Cell Biology.

[34]  C. Peskin The immersed boundary method , 2002, Acta Numerica.

[35]  Leah Edelstein-Keshet,et al.  Phosphoinositides and Rho proteins spatially regulate actin polymerization to initiate and maintain directed movement in a one-dimensional model of a motile cell. , 2007, Biophysical journal.

[36]  C. Lawson,et al.  ICAM-1 signaling in endothelial cells , 2009, Pharmacological reports : PR.

[37]  J. Sethian,et al.  Transport and diffusion of material quantities on propagating interfaces via level set methods , 2003 .

[38]  Jung-Chi Liao,et al.  Extending the absorbing boundary method to fit dwell-time distributions of molecular motors with complex kinetic pathways , 2007, Proceedings of the National Academy of Sciences.

[39]  H. Meinhardt,et al.  A theory of biological pattern formation , 1972, Kybernetik.

[40]  Carol S. Woodward,et al.  Enabling New Flexibility in the SUNDIALS Suite of Nonlinear and Differential/Algebraic Equation Solvers , 2020, ACM Trans. Math. Softw..

[41]  Alex Mogilner,et al.  Multiscale Two-Dimensional Modeling of a Motile Simple-Shaped Cell , 2005, Multiscale Model. Simul..

[42]  David Trebotich,et al.  Embedded boundary algorithms and software for partial differential equations , 2008 .

[43]  L. Rajendran,et al.  Lipid rafts and membrane dynamics , 2005, Journal of Cell Science.

[44]  Jaap A. Kaandorp,et al.  Computational methods for diffusion-influenced biochemical reactions , 2007, Bioinform..

[45]  A. M. Turing,et al.  The chemical basis of morphogenesis , 1952, Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences.

[46]  Peter McCorquodale,et al.  A Cartesian grid embedded boundary method for the heat equation on irregular domains , 2001 .

[47]  O. Petersen,et al.  Transformation of local Ca2+ spikes to global Ca2+ transients: the combinatorial roles of multiple Ca2+ releasing messengers , 2002, The EMBO journal.

[48]  Carolyn R. Bertozzi,et al.  The Glycosylphosphatidylinositol Anchor: A Complex Membrane-Anchoring Structure for Proteins† , 2008, Biochemistry.

[49]  S. Su,et al.  The critical role of Toll-like receptor signaling pathways in the induction and progression of autoimmune diseases. , 2009, Current molecular medicine.

[50]  J. Tyson,et al.  Numerical analysis of a comprehensive model of M-phase control in Xenopus oocyte extracts and intact embryos. , 1993, Journal of cell science.

[51]  John J. Tyson,et al.  Hysteresis drives cell-cycle transitions in Xenopus laevis egg extracts , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[52]  Richard A. Davis,et al.  Time Series: Theory and Methods (2nd ed.). , 1992 .

[53]  Alex M. Andrew,et al.  Level Set Methods and Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid Mechanics, Computer Vision, and Materials Science (2nd edition) , 2000 .

[54]  E. Rouy,et al.  A viscosity solutions approach to shape-from-shading , 1992 .

[55]  E R Kandel,et al.  Spatially resolved dynamics of cAMP and protein kinase A subunits in Aplysia sensory neurons. , 1993, Science.

[56]  P. Colella,et al.  A Cartesian Grid Embedded Boundary Method for Poisson's Equation on Irregular Domains , 1998 .

[57]  Randall J. LeVeque,et al.  Finite difference methods for ordinary and partial differential equations - steady-state and time-dependent problems , 2007 .

[58]  Zhilin Li,et al.  The Immersed Interface Method: Numerical Solutions of PDEs Involving Interfaces and Irregular Domains (Frontiers in Applied Mathematics) , 2006 .

[59]  K. Hahn,et al.  Activation of Endogenous Cdc42 Visualized in Living Cells , 2004, Science.

[60]  H. Pelham,et al.  Slow Diffusion of Proteins in the Yeast Plasma Membrane Allows Polarity to Be Maintained by Endocytic Cycling , 2003, Current Biology.

[61]  John B. Moore,et al.  Hidden Markov Models: Estimation and Control , 1994 .

[62]  K. Hahn,et al.  Localized Rac activation dynamics visualized in living cells. , 2000, Science.

[63]  Krister Wennerberg,et al.  Rho and Rac Take Center Stage , 2004, Cell.

[64]  P. Colella,et al.  Embedded boundary grid generation using the divergence theorem, implicit functions, and constructive solid geometry , 2008 .

[65]  Ronald Fedkiw,et al.  The immersed interface method. Numerical solutions of PDEs involving interfaces and irregular domains , 2007, Math. Comput..

[66]  S. Osher,et al.  Uniformly high order accurate essentially non-oscillatory schemes, 111 , 1987 .

[67]  Linda R. Petzold,et al.  Using Krylov Methods in the Solution of Large-Scale Differential-Algebraic Systems , 1994, SIAM J. Sci. Comput..

[68]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .

[69]  Linda R. Petzold,et al.  Numerical solution of initial-value problems in differential-algebraic equations , 1996, Classics in applied mathematics.

[70]  James C. Schaff,et al.  Analysis of nonlinear dynamics on arbitrary geometries with the Virtual Cell. , 2001, Chaos.

[71]  Alexandra Jilkine,et al.  Polarization and Movement of Keratocytes: A Multiscale Modelling Approach , 2006, Bulletin of mathematical biology.

[72]  Attila Csikász-Nagy,et al.  Analysis of a generic model of eukaryotic cell-cycle regulation. , 2006, Biophysical journal.

[73]  Jonathan Richard Shewchuk,et al.  Triangle: Engineering a 2D Quality Mesh Generator and Delaunay Triangulator , 1996, WACG.

[74]  Jason M Haugh,et al.  Spatial analysis of 3' phosphoinositide signaling in living fibroblasts, III: influence of cell morphology and morphological Polarity. , 2005, Biophysical journal.

[75]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .