Spiking at the edge: Excitability at interfaces in reaction–diffusion systems

Edited by Jacques Prost, Institut Curie, Paris, France; received May 12, 2023; accepted October 25, 2023 by Editorial Board Member Paul Chaikin
January 12, 2024
121 (3) e2307996120

Significance

Spiking is a general phenomenon that is crucial in the firing of neurons, beating of hearts, and spread of diseases. In homogeneous media, spiking arises from a local competition between amplifying and suppressing forces. But most real-world systems are far from homogeneous. Here, we demonstrate that inhomogeneities such as interfaces and boundaries (that spatially segregate these two forces) can promote spiking, even if the system does not spike when these forces are evenly mixed. We mathematically derive a spiking phase diagram in terms of interfacial diffusion and amplification strength. Our findings apply to chemical reactions, predator–prey dynamics, and recent electrophysiology experiments, in which localized action potentials were observed at the interface of distinct, nonspiking bioelectric tissues.

Abstract

Excitable media, ranging from bioelectric tissues and chemical oscillators to forest fires and competing populations, are nonlinear, spatially extended systems capable of spiking. Most investigations of excitable media consider situations where the amplifying and suppressing forces necessary for spiking coexist at every point in space. In this case, spikes arise due to local bistabilities, which require a fine-tuned ratio between local amplification and suppression strengths. But, in nature and engineered systems, these forces can be segregated in space, forming structures like interfaces and boundaries. Here, we show how boundaries can generate and protect spiking when the reacting components can spread out: Even arbitrarily weak diffusion can cause spiking at the edge between two non-excitable media. This edge spiking arises due to a global bistability, which can occur even if amplification and suppression strengths do not allow spiking when mixed. We analytically derive a spiking phase diagram that depends on two parameters: i) the ratio between the system size and the characteristic diffusive length-scale and ii) the ratio between the amplification and suppression strengths. Our analysis explains recent experimental observations of action potentials at the interface between two non-excitable bioelectric tissues. Beyond electrophysiology, we highlight how edge spiking emerges in predator–prey dynamics and in oscillating chemical reactions. Our findings provide a theoretical blueprint for a class of interfacial excitations in reaction–diffusion systems, with potential implications for spatially controlled chemical reactions, nonlinear waveguides and neuromorphic computation, as well as spiking instabilities, such as cardiac arrhythmias, that naturally occur in heterogeneous biological media.

Continue Reading

Data, Materials, and Software Availability

All data and code supporting this article are available in Zenodo (DOI: https://www.doi.org/10.5281/zenodo.10426295) (62).

Acknowledgments

We would like to thank William Bialek, Ned Wingreen, Joshua W. Shaevitz, David J. Schwab, James P. Sethna, Itai Cohen, Eric R. Dufresne, Erwin Frey, Fridtjof Brauns, Pankaj Mehta, Mary Silber, Brent Dorion, Massimo Vergassola, Suraj Shankar, Noah P. Mitchell, Danny S. Seara, Michel Fruchart, and Tali Khain for useful discussions. C.S. acknowledges support from the Bloomenthal Fellowship and the NSF Graduate Research Fellowship under Grant No. 1746045. V.V. acknowledges support from the Simons Foundation, the Army Research Office under grants W911NF-22-2-0109 and W911NF-23-1-0212, and the University of Chicago Materials Research Science and Engineering Center, which is funded by the NSF under Award No. DMR-2011854. This research was also supported from the National Science Foundation under grant DMR-2118415 and through the Center for Living Systems (grant no. 2317138). A.E.C. acknowledges support from a Vannevar Bush Faculty Fellowship grant N00014-18-1-2859, and NSF Quantum Sensing for Biophysics and Bioengineering (QuBBE) Quantum Leap Challenge Institute (QLCI) grant OMA-2121044. H.O. acknowledges support from an European Molecular Biology Organization (EMBO) Fellowship ALTF 543-2020.

Author contributions

C.S., H.O., A.E.C., and V.V. designed research; C.S., H.O., A.E.C., and V.V. performed research; C.S. and H.O. analyzed data; and C.S., H.O., A.E.C., and V.V. wrote the paper.

Competing interests

The authors declare no competing interest.

Supporting Information

Appendix 01 (PDF)
Movie S1.
Experimental observation of spiking at a bioelectric interface: Two tissues of human embryonic kidney (HEK293) cells were genetically modified to express either sodium (NaV1.5) or potassium (Kir2.1) channels. Neither tissue alone is able to spike. However, upon stimulation by a laser, an action potential is observed to propagate along their interface, as revealed by a voltage sensitive red dye. Adapted from Ref. (8).
Movie S2.
2D nonlinear waveguides built from excitable interfaces: A network of excitable interfaces forms a two-dimensional nonlinear waveguide capable of performing computations. Simulations of the binary half adder in Fig. 5c-d are shown. See SI §5 for details of the numerics.

References

1
A. T. Winfree, Electrical turbulence in three-dimensional heart muscle. Science 266, 1003–1006 (1994).
2
D. M. Bers, Cardiac excitation-contraction coupling. Nature 415, 198–205 (2002).
3
K. H. W. J. ten Tusscher, D. Noble, P. J. Noble, A. V. Panfilov, A model for human ventricular tissue. Am. J. Physiol.-Heart Circul. Physiol. 286, H1573–H1589 (2004).
4
H. Cheng, W. J. Lederer, M. B. Cannell, Calcium sparks: Elementary events underlying excitation-contraction coupling in heart muscle. Science 262, 740–744 (1993).
5
M. Stern, Theory of excitation-contraction coupling in cardiac muscle. Biophys. J. 63, 497–517 (1992).
6
F. X. Witkowski et al., Spatiotemporal evolution of ventricular fibrillation. Nature 392, 78–82 (1998).
7
F. Rieke, D. Warland, R. De Ruyter van Steveninck, W. Bialek, Spikes: Exploring the Neural Code (Bradford Book/MIT Press, 1997).
8
B. Drossel, F. Schwabl, Self-organized critical forest-fire model. Phys. Rev. Lett. 69, 1629–1632 (1992).
9
R. M. Anderson, R. M. May, Population biology of infectious diseases: Part i. Nature 280, 361–367 (1979).
10
J. D. Murray, E. A. Stanley, D. L. Brown, On the spatial spread of rabies among foxes. Proc. R. Soc. London. Ser. B. Biol. Sci. 229, 111–150 (1986).
11
R. M. Anderson, H. C. Jackson, R. M. May, A. M. Smith, Population dynamics of fox rabies in Europe. Nature 289, 765–771 (1981).
12
P. Rohani, D. J. D. Earn, B. T. Grenfell, Opposite patterns of synchrony in sympatric disease metapopulations. Science 286, 968–971 (1999).
13
S. Kondo, T. Miura, Reaction-diffusion model as a framework for understanding biological pattern formation. Science 329, 1616–1620 (2010).
14
J. A. Bourret, R. G. Lincoln, B. H. Carpenter, Fungal endogenous rhythms expressed by spiral figures. Science 166, 763–764 (1969).
15
M. Loose, E. Fischer-Friedrich, J. Ries, K. Kruse, P. Schwille, Spatial regulators for bacterial cell division self-organize into surface waves in vitro. Science 320, 789–792 (2008).
16
N. Tompkins et al., Testing Turing’s theory of morphogenesis in chemical cells. Proc. Natl. Acad. Sci. U.S.A. 111, 4397–4402 (2014).
17
H. H. Rotermund, W. Engel, M. Kordesch, G. Ertl, Imaging of spatio-temporal pattern evolution during carbon monoxide oxidation on platinum. Nature 343, 355–357 (1990).
18
O. Steinbock, P. Kettunen, K. Showalter, Anisotropy and spiral organizing centers in patterned excitable media. Science 269, 1857–1860 (1995).
19
M. Vinson, S. Mironov, S. Mulvey, A. Pertsov, Control of spatial orientation and lifetime of scroll rings in excitable media. Nature 386, 477–480 (1997).
20
Y. Fuseya, H. Katsuno, K. Behnia, A. Kapitulnik, Nanoscale Turing patterns in a bismuth monolayer. Nat. Phys. 17, 1031–1036 (2021).
21
T. H. Tan et al., Topological turbulence in the membrane of a living cell. Nat. Phys. 16, 657–662 (2020).
22
A. T. Winfree, Persistent tangled vortex rings in generic excitable media. Nature 371, 233–236 (1994).
23
W. Van Saarloos, V. Vitelli, Z. Zeravcic, Soft Matter: Concepts, Phenomena and Applications (Princeton University Press, 2023).
24
A. M. Turing, The chemical basis of morphogenesis. Philos. Trans. R. Soc. London. Ser. B, Biol. Sci. 237, 37–72 (1952).
25
A. Nakamasu, G. Takahashi, A. Kanbe, S. Kondo, Interactions between zebrafish pigment cells responsible for the generation of Turing patterns. Proc. Natl. Acad. Sci. U.S.A. 106, 8429–8434 (2009).
26
S. Kondo, R. Asai, A reaction-diffusion wave on the skin of the marine angelfish Pomacanthus. Nature 376, 765–768 (1995).
27
S. A. Newman, H. L. Frisch, Dynamics of skeletal pattern formation in developing chick limb. Science 205, 662–668 (1979).
28
N. P. Mitchell et al., Visceral organ morphogenesis via calcium-patterned muscle constrictions. eLife 11, e77355 (2022).
29
M. C. Wigbers et al., A hierarchy of protein patterns robustly decodes cell shape information. Nat. Phys. 17, 578–584 (2021).
30
S. Di Talia, M. Vergassola, Waves in embryonic development. Annu. Rev. Biophys. 51, 327–353. 35119944 (2022).
31
M. Vergassola, V. E. Deneke, S. D. Talia, Mitotic waves in the early embryogenesis of Drosophila: Bistability traded for speed. Proc. Natl. Acad. Sci. U.S.A. 115, E2165–E2174 (2018).
32
J. Lechleiter, S. Girard, E. Peralta, D. Clapham, Spiral calcium wave propagation and annihilation in Xenopus laevis oocytes. Science 252, 123–126 (1991).
33
J. B. Michaux, F. B. Robin, W. M. McFadden, E. M. Munro, Excitable RhoA dynamics drive pulsed contractions in the early C. elegans embryo. J. Cell Biol. 217, 4230–4252 (2018).
34
J. B. Chang, J. E. Ferrell Jr., Mitotic trigger waves and the spatial coordination of the Xenopus cell cycle. Nature 500, 603–607 (2013).
35
J. M. Davidenko, A. V. Pertsov, R. Salomonsz, W. Baxter, J. Jalife, Stationary and drifting spiral waves of excitation in isolated cardiac muscle. Nature 355, 349–351 (1992).
36
R. A. Gray, A. M. Pertsov, J. Jalife, Spatial and temporal organization during cardiac fibrillation. Nature 392, 75–78 (1998).
37
F. H. Fenton, E. M. Cherry, H. M. Hastings, S. J. Evans, Multiple mechanisms of spiral wave breakup in a model of cardiac electrical activity. Chaos: Interdiscip. J. Nonlinear Sci. 12, 852–892 (2002).
38
J. Halatek, E. Frey, Rethinking pattern formation in reaction-diffusion systems. Nat. Phys. 14, 507–514 (2018).
39
M. C. Cross, P. C. Hohenberg, Pattern formation outside of equilibrium. Rev. Mod. Phys. 65, 851–1112 (1993).
40
M. Kim et al., Controlling chemical turbulence by global delayed feedback: Pattern formation in catalytic CO oxidation on Pt(110). Science 292, 1357–1360 (2001).
41
F. Brauns, J. Halatek, E. Frey, Phase-space geometry of mass-conserving reaction-diffusion dynamics. Phys. Rev. X 10, 041036 (2020).
42
S. Alonso, F. Sagués, A. S. Mikhailov, Taming Winfree turbulence of scroll waves in excitable media. Science 299, 1722–1725 (2003).
43
H. M. McNamara et al., Bioelectrical domain walls in homogeneous tissues. Nat. Phys. 16, 357–364 (2020).
44
T. Eckstein, E. Vidal-Henriquez, A. Gholami, Experimental observation of boundary-driven oscillations in a reaction-diffusion-advection system. Soft Matter 16, 4243–4255 (2020).
45
E. Vidal-Henriquez, V. Zykov, E. Bodenschatz, A. Gholami, Convective instability and boundary driven oscillations in a reaction-diffusion-advection model. Chaos: Interdiscip. J. Nonlinear Sci. 27, 103110 (2017).
46
W. M. Ni, J. Wei, On the location and profile of spike-layer solutions to singularly perturbed semilinear Dirichlet problems. Commun. Pure Appl. Math. 48, 731–768 (1995).
47
G. Bub, A. Shrier, L. Glass, Spiral wave generation in heterogeneous excitable media. Phys. Rev. Lett. 88, 058101 (2002).
48
Z. F. Mainen, T. J. Sejnowski, Influence of dendritic structure on firing pattern in model neocortical neurons. Nature 382, 363–366 (1996).
49
M. C. Wigbers, F. Brauns, T. Hermann, E. Frey, Pattern localization to a domain edge. Phys. Rev. E 101, 022414 (2020).
50
F. Brauns, J. Halatek, E. Frey, Diffusive coupling of two well-mixed compartments elucidates elementary principles of protein-based pattern formation. Phys. Rev. Res. 3, 013258 (2021).
51
G. Bub, A. Shrier, L. Glass, Spiral wave generation in heterogeneous excitable media. Phys. Rev. Lett. 88, 058101 (2002).
52
K. Agladze, J. P. Keener, S. C. Müller, A. Panfilov, Rotating spiral waves created by geometry. Science 264, 1746–1748 (1994).
53
M. F. Staddon, E. M. Munro, S. Banerjee, Pulsatile contractions and pattern formation in excitable actomyosin cortex. PLoS Comput. Biol. 18, 1–21 (2022).
54
S. Luther et al., Low-energy control of electrical turbulence in the heart. Nature 475, 235–239 (2011).
55
A. Murugan, S. Vaikuntanathan, Topologically protected modes in non-equilibrium stochastic systems. Nat. Commun. 8, 13881 (2017).
56
C. L. Kane, T. C. Lubensky, Topological boundary modes in isostatic lattices. Nat. Phys. 10, 39–45 (2014).
57
M. Z. Hasan, C. L. Kane, Colloquium: Topological insulators. Rev. Mod. Phys. 82, 3045–3067 (2010).
58
S. Shankar, A. Souslov, M. J. Bowick, M. C. Marchetti, V. Vitelli, Topological active matter. Nat. Rev. Phys. 4, 380–398 (2022).
59
B. G. ge Chen, N. Upadhyaya, V. Vitelli, Nonlinear conduction via solitons in a topological mechanical insulator. Proc. Natl. Acad. Sci. U.S.A. 111, 13004–13009 (2014).
60
X. Mao, T. C. Lubensky, Maxwell lattices and topological mechanics. Ann. Rev. Condens. Matter Phys. 9, 413–433 (2018).
61
S. D. Huber, Topological mechanics. Nat. Phys. 12, 621–623 (2016).
62
C. Scheibner, H. Ori, A. E. Cohen, V. Vitelli, Code and data for “Spiking at the edge: Excitability at interfaces in reaction-diffusion systems”. Zenodo. https://doi.org/10.5281/zenodo.10426295. Deposited 26 December 2023.
63
H. Ori et al., Observation of topological action potentials in engineered tissues. Nat. Phys. 19, 290–296 (2023).
64
R. FitzHugh, Impulses and physiological states in theoretical models of nerve membrane. Biophys. J. 1, 445–466 (1961).
65
C. Xu et al., Computational design of transmembrane pores. Nature 585, 129–134 (2020).
66
J. Payandeh, T. Scheuer, N. Zheng, W. A. Catterall, The crystal structure of a voltage-gated sodium channel. Nature 475, 353–358 (2011).
67
C. C. Conley, J. A. Smoller, Algebraic and Topological Invariants for Reaction-Diffusion Equations, J. M. Ball, Eds. (Springer, Netherlands, Dordrecht, 1983), pp. 3–24.
68
K. Mischaikow, M. Mrozek, “Conley index” in Handbook of Dynamical Systems, B. Fiedler, Ed. (Elsevier Science, 2002), vol. 2, pp. 393–460.
69
J. J. Tyson, The Belousov-Zhabotinskii Reaction (Springer-Verlag, Berlin, New York, 1976).
70
M. A. Budroni, L. Lemaigre, D. M. Escala, A. P. Muñuzuri, A. De Wit, Spatially localized chemical patterns around an a + b → oscillator front J. Phys. Chem. A 120, 851–860 (2016).
71
B. Dúzs, P. De Kepper, I. Szalai, Turing patterns and waves in closed two-layer gel reactors. ACS Omega 4, 3213–3219 (2019).
72
S. N. Semenov et al., Autocatalytic, bistable, oscillatory networks of biologically relevant organic reactions. Nature 537, 656–660 (2016).
73
R. Yoshida, Self-oscillating gels driven by the Belousov-Zhabotinsky reaction as novel smart materials. Adv. Mater. 22, 3463–3483 (2010).
74
G. Rabai, K. Kustin, I. R. Epstein, A systematically designed pH oscillator: The hydrogen peroxide-sulfite-ferrocyanide reaction in a continuous-flow stirred tank reactor. J. Am. Chem. Soc. 111, 3870–3874 (1989).
75
A. Testa et al., Sustained enzymatic activity and flow in crowded protein droplets. Nat. Commun. 12, 6293 (2021).
76
A. Adamatzky, B. De Lacy Costello, T. Asai, Reaction-Diffusion Computers (Elsevier Science, 2005).
77
J. Holley, I. Jahan, B. De Lacy Costello, L. Bull, A. Adamatzky, Logical and arithmetic circuits in Belousov-Zhabotinsky encapsulated disks. Phys. Rev. E 84, 056110 (2011).
78
Á. Tóth, K. Showalter, Logic gates in excitable media. J. Chem. Phys. 103, 2058–2066 (1995).
79
O. Steinbock, Á. Tóth, K. Showalter, Navigating complex labyrinths: Optimal paths from chemical waves. Science 267, 868–871 (1995).
80
H. M. McNamara et al., Geometry-dependent arrhythmias in electrically excitable tissues. Cell Syst. 7, 359–370.e6 (2018).
81
E. M. Cherry, J. R. Ehrlich, S. Nattel, F. H. Fenton, Pulmonary vein reentry-properties and size matter: Insights from a computational analysis. HeartRythm 4, 1553–1562 (2007).

Information & Authors

Information

Published in

The cover image for PNAS Vol.121; No.3
Proceedings of the National Academy of Sciences
Vol. 121 | No. 3
January 16, 2024
PubMed: 38215183

Classifications

Data, Materials, and Software Availability

All data and code supporting this article are available in Zenodo (DOI: https://www.doi.org/10.5281/zenodo.10426295) (62).

Submission history

Received: May 12, 2023
Accepted: October 25, 2023
Published online: January 12, 2024
Published in issue: January 16, 2024

Keywords

  1. reaction–diffusion equations
  2. dynamical systems
  3. electrophysiology
  4. morphogenesis
  5. topological phenomena

Acknowledgments

We would like to thank William Bialek, Ned Wingreen, Joshua W. Shaevitz, David J. Schwab, James P. Sethna, Itai Cohen, Eric R. Dufresne, Erwin Frey, Fridtjof Brauns, Pankaj Mehta, Mary Silber, Brent Dorion, Massimo Vergassola, Suraj Shankar, Noah P. Mitchell, Danny S. Seara, Michel Fruchart, and Tali Khain for useful discussions. C.S. acknowledges support from the Bloomenthal Fellowship and the NSF Graduate Research Fellowship under Grant No. 1746045. V.V. acknowledges support from the Simons Foundation, the Army Research Office under grants W911NF-22-2-0109 and W911NF-23-1-0212, and the University of Chicago Materials Research Science and Engineering Center, which is funded by the NSF under Award No. DMR-2011854. This research was also supported from the National Science Foundation under grant DMR-2118415 and through the Center for Living Systems (grant no. 2317138). A.E.C. acknowledges support from a Vannevar Bush Faculty Fellowship grant N00014-18-1-2859, and NSF Quantum Sensing for Biophysics and Bioengineering (QuBBE) Quantum Leap Challenge Institute (QLCI) grant OMA-2121044. H.O. acknowledges support from an European Molecular Biology Organization (EMBO) Fellowship ALTF 543-2020.
Author contributions
C.S., H.O., A.E.C., and V.V. designed research; C.S., H.O., A.E.C., and V.V. performed research; C.S. and H.O. analyzed data; and C.S., H.O., A.E.C., and V.V. wrote the paper.
Competing interests
The authors declare no competing interest.

Notes

This article is a PNAS Direct Submission. J.P. is a guest editor invited by the Editorial Board.
*
Notice that Eqs. 3 and 4 do not contain advective transport, which has also been shown to give rise oscillations near Dirichlet boundaries, for example, in models of and experiments on Dictyostelium discoideum (44, 45).
For simplicity, in this example, we are using the same hopping rate ϵ within the forest as between the forest and desert. This distinction becomes irrelevant in the continuum limit (large ϵ and large N) because this subextensive heterogeneity is absorbed into the Dirichlet boundary condition at an edge or into the continuity requirements across an interface.

Authors

Affiliations

Colin Scheibner
Department of Physics and The James Franck Institute, The University of Chicago, Chicago, IL 60637
Kadanoff Center for Theoretical Physics, The University of Chicago, Chicago, IL 60637
Hillel Ori
Department of Chemistry and Chemical Biology, Harvard University, Cambridge, MA 02138
Adam E. Cohen
Department of Chemistry and Chemical Biology, Harvard University, Cambridge, MA 02138
Department of Physics, Harvard University, Cambridge, MA 02138
Department of Physics and The James Franck Institute, The University of Chicago, Chicago, IL 60637
Kadanoff Center for Theoretical Physics, The University of Chicago, Chicago, IL 60637
Institute for Biophysical Dynamics, The University of Chicago, Chicago, IL 60637

Notes

1
To whom correspondence may be addressed. Email: [email protected].

Metrics & Citations

Metrics

Note: The article usage is presented with a three- to four-day delay and will update daily once available. Due to ths delay, usage data will not appear immediately following publication. Citation information is sourced from Crossref Cited-by service.


Altmetrics

Citations

Export the article citation data by selecting a format from the list below and clicking Export.

Cited by

    Loading...

    View Options

    View options

    PDF format

    Download this article as a PDF file

    DOWNLOAD PDF

    Login options

    Check if you have access through your login credentials or your institution to get full access on this article.

    Personal login Institutional Login

    Recommend to a librarian

    Recommend PNAS to a Librarian

    Purchase options

    Purchase this article to access the full text.

    Single Article Purchase

    Spiking at the edge: Excitability at interfaces in reaction–diffusion systems
    Proceedings of the National Academy of Sciences
    • Vol. 121
    • No. 3

    Figures

    Tables

    Media

    Share

    Share

    Share article link

    Share on social media