Simulating Biochemical Signaling Networks in Complex Moving Geometries

Signaling networks regulate cellular responses to environmental stimuli through cascades of protein interactions. 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 developed a numerical scheme consisting of a cut cell finite volume spatial discretization coupled with level set methods to simulate the resulting advection-reaction-diffusion system. We then apply the method to several biochemical reaction networks in changing geometries. We found that a Turing instability can develop exclusively by cell deformations that maintain constant area. For a Turing system with a geometry-dependent single or double peak solution, simulations in a dynamically changing geometry suggest that a single peak solution is the only stable one, independent of the oscillation frequency. The method is also applied to a model of a signaling network in a migrating fibroblast.

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

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

[3]  Ronald Fedkiw,et al.  Level set methods and dynamic implicit surfaces , 2002, Applied mathematical sciences.

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

[5]  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.

[6]  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.

[7]  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.

[8]  S. Osher,et al.  Uniformly High-Order Accurate Nonoscillatory Schemes. I , 1987 .

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

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

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

[12]  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.

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

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

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

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

[17]  Timothy C Elston,et al.  A Cut Cell Method for Simulating Spatial Models of Biochemical Reaction Networks in Arbitrary Geometries. , 2010, Communications in applied mathematics and computational science.

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

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

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

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

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

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

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

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

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

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

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

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

[30]  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.

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

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

[33]  Jeff Hasty,et al.  Yeast Dynamically Modify Their Environment to Achieve Better Mating Efficiency , 2011, Science Signaling.

[34]  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.

[35]  B. Alberts,et al.  Molecular Biology of the Cell (Fifth Edition) , 2008 .

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

[37]  O. C. Zienkiewicz,et al.  The Finite Element Method: Its Basis and Fundamentals , 2005 .

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

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

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

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

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

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

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

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