Multiscale Analysis and Simulation of a Signaling Process With Surface Diffusion

We present and analyse a model for cell signalling processes in biological tissues. The model includes diffusion and nonlinear reactions on the cell surfaces, and both inter- and intracellular signalling. Using techniques from the theory of two-scale convergence as well the unfolding method, we show convergence of the solutions to the model to solutions of a two-scale macroscopic problem. We also present a two-scale bulk-surface finite element method for the approximation of the macroscopic model and we report on some benchmarking results as well as numerical simulations in a biologically relevant regime that illustrate the influence of cell-scale heterogeneities on the macroscopic concentrations.

[1]  J A Sherratt,et al.  Keratinocyte growth factor signalling: a mathematical model of dermal-epidermal interaction in epidermal wound healing. , 2000, Mathematical biosciences.

[2]  C. M. Elliott,et al.  Finite element analysis for a coupled bulk-surface partial differential equation , 2013 .

[3]  M. A. Alonso,et al.  The role of lipid rafts in signalling and membrane trafficking in T lymphocytes. , 2001, Journal of cell science.

[4]  Dagmar Iber,et al.  Kidney branching morphogenesis under the control of a ligand–receptor-based Turing mechanism , 2013, Physical biology.

[5]  Charles M. Elliott,et al.  A Coupled Ligand-Receptor Bulk-Surface System on a Moving Domain: Well Posedness, Regularity, and Convergence to Equilibrium , 2016, SIAM J. Math. Anal..

[6]  John A. Mackenzie,et al.  A computational method for the coupled solution of reaction–diffusion equations on evolving domains and manifolds: Application to a model of cell migration and chemotaxis , 2016, J. Comput. Phys..

[7]  Christoph Schwab,et al.  High-Dimensional Finite Elements for Elliptic Problems with Multiple Scales , 2005, Multiscale Modeling & simulation.

[8]  M. Ptashnyk,et al.  Unfolding-based corrector estimates for a reaction–diffusion system predicting concrete corrosion , 2011, 1106.3233.

[9]  Anna K. Marciniak-Czochra,et al.  Derivation of a Macroscopic Receptor-Based Model Using Homogenization Techniques , 2008, SIAM J. Math. Anal..

[10]  G. Galdi.  An Introduction to the Mathematical Theory of the Navier-Stokes Equations : Volume I: Linearised Steady Problems , 1994 .

[11]  O. Lakkis,et al.  RATE OF CONVERGENCE FOR A GALERKIN SCHEME APPROXIMATING A TWO-SCALE REACTION-DIFFUSION SYSTEM WITH NONLINEAR TRANSMISSION CONDITION (Nonlinear evolution equations and mathematical modeling) , 2010, 1002.3793.

[12]  Matthias Röger,et al.  Turing instabilities in a mathematical model for signaling networks , 2012, Journal of mathematical biology.

[13]  Erik S. Welf,et al.  Quantitative statistical description of integrin clusters in adherent cells. , 2009, IET systems biology.

[14]  Doina Cioranescu,et al.  Homogenization of Reticulated Structures , 1999 .

[15]  G. Allaire.  Homogenization and two-scale convergence , 1992 .

[16]  A. Marciniak-Czochra.  Receptor-Based Models with Diffusion-Driven Instability for Pattern Formation in Hydra , 2003 .

[17]  Maria Neuss-Radu,et al.  Some extensions of two-scale convergence , 1996 .

[18]  G. Dal Maso,et al.  An extension theorem from connected sets, and homogenization in general periodic domains , 1992 .

[19]  J. Groves,et al.  Signaling clusters in the cell membrane. , 2011, Current opinion in cell biology.

[20]  H. Brezis.  Functional Analysis, Sobolev Spaces and Partial Differential Equations , 2010 .

[21]  Kai Simons,et al.  Lipid rafts and signal transduction , 2000, Nature Reviews Molecular Cell Biology.

[22]  Doina Cioranescu,et al.  The Periodic Unfolding Method in Homogenization , 2008, SIAM J. Math. Anal..

[23]  J. Mackenzie,et al.  Local modulation of chemoattractant concentrations by single cells: dissection using a bulk-surface computational model , 2016, Interface Focus.

[24]  J. Simon.  Compact sets in the spaceLp(O,T; B) , 1986 .

[25]  Giovanni P. Galdi,et al.  An Introduction to the Mathematical Theory of the Navier-Stokes Equations: Steady-State Problems , 2011 .

[26]  Charles M. Elliott,et al.  Coupled Bulk-Surface Free Boundary Problems Arising from a Mathematical Model of Receptor-Ligand Dynamics , 2015, SIAM J. Math. Anal..

[27]  Roland Eils,et al.  Integration of Activating and Inhibitory Receptor Signaling by Regulated Phosphorylation of Vav1 in Immune Cells , 2011, Science Signaling.

[28]  C. Venkataraman,et al.  Stability analysis and simulations of coupled bulk-surface reaction–diffusion systems , 2015, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[29]  Dagmar Iber,et al.  An interplay of geometry and signaling enables robust lung branching morphogenesis , 2014, Development.

[30]  J. Gálvez,et al.  Mathematical modelling and computational study of two-dimensional and three-dimensional dynamics of receptor–ligand interactions in signalling response mechanisms , 2013, Journal of Mathematical Biology.

[31]  Doina Cioranescu,et al.  The Periodic Unfolding Method in Domains with Holes , 2012, SIAM J. Math. Anal..

[32]  H. Wiley,et al.  Relationship between epidermal growth factor receptor occupancy and mitogenic response. Quantitative analysis using a steady state model system. , 1984, The Journal of biological chemistry.

[33]  Malte A. Peter,et al.  Diffusion on Surfaces and the Boundary Periodic Unfolding Operator with an Application to Carcinogenesis in Human Cells , 2014, SIAM J. Math. Anal..

[34]  Peter Wall,et al.  Two-scale convergence , 2002 .

[35]  K. Gaus,et al.  Domain-specific lipid distribution in macrophage plasma membranes Published, JLR Papers in Press, May 1, 2005. DOI 10.1194/jlr.M500103-JLR200 , 2005, Journal of Lipid Research.

[36]  C. Marshall,et al.  Regulation of cell signalling by uPAR , 2010, Nature Reviews Molecular Cell Biology.

[37]  D A Lauffenburger,et al.  Analysis of intracellular receptor/ligand sorting. Calculation of mean surface and bulk diffusion times within a sphere. , 1986, Biophysical journal.

[38]  Willi Jäger,et al.  Diffusion, convection, adsorption, and reaction of chemicals in porous media , 1991 .

[39]  K. Burridge,et al.  The on-off relationship of Rho and Rac during integrin-mediated adhesion and cell migration , 2014, Small GTPases.

[40]  Steven J. Ruuth,et al.  Simple computation of reaction–diffusion processes on point clouds , 2013, Proceedings of the National Academy of Sciences.

[41]  S. Clouse.  Brassinosteroid Signal Transduction: From Receptor Kinase Activation to Transcriptional Networks Regulating Plant Development , 2011, Plant Cell.

[42]  C. M. Elliott,et al.  Modelling cell motility and chemotaxis with evolving surface finite elements , 2012, Journal of The Royal Society Interface.

[43]  Dagmar Iber,et al.  Feedback, receptor clustering, and receptor restriction to single cells yield large Turing spaces for ligand-receptor-based Turing models. , 2014, Physical review. E, Statistical, nonlinear, and soft matter physics.

[44]  Omar Lakkis,et al.  Implicit-Explicit Timestepping with Finite Element Approximation of Reaction-Diffusion Systems on Evolving Domains , 2011, SIAM J. Numer. Anal..

[45]  J. Stock,et al.  Signal Transduction: Receptor Clusters as Information Processing Arrays , 2002, Current Biology.

[46]  Jussi Timonen,et al.  Protein Diffusion in Mammalian Cell Cytoplasm , 2011, PloS one.

[47]  G. Nguetseng.  A general convergence result for a functional related to the theory of homogenization , 1989 .

[48]  V. Thomée.  Galerkin Finite Element Methods for Parabolic Problems (Springer Series in Computational Mathematics) , 2010 .

[49]  O. A. Ladyzhenskai︠a︡,et al.  Linear and Quasi-linear Equations of Parabolic Type , 1995 .