1
|
Buenzli PR, Simpson MJ. Curvature dependences of wave propagation in reaction–diffusion models. Proc Math Phys Eng Sci 2022. [DOI: 10.1098/rspa.2022.0582] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
Reaction–diffusion waves in multiple spatial dimensions advance at a rate that strongly depends on the curvature of the wavefronts. These waves have important applications in many physical, ecological and biological systems. In this work, we analyse curvature dependences of travelling fronts in a single reaction–diffusion equation with general reaction term. We derive an exact, non-perturbative curvature dependence of the speed of travelling fronts that arises from transverse diffusion occurring parallel to the wavefront. Inward-propagating waves are characterized by three phases: an establishment phase dominated by initial and boundary conditions, a travelling-wave-like phase in which normal velocity matches standard results from singular perturbation theory and a dip-filling phase where the collision and interaction of fronts create additional curvature dependences to their progression rate. We analyse these behaviours and additional curvature dependences using a combination of asymptotic analyses and numerical simulations.
Collapse
Affiliation(s)
- Pascal R. Buenzli
- School of Mathematical Sciences, Queensland University of Technology (QUT), Brisbane, Queensland 4000, Australia
| | - Matthew J. Simpson
- School of Mathematical Sciences, Queensland University of Technology (QUT), Brisbane, Queensland 4000, Australia
| |
Collapse
|
2
|
Kuze M, Kitahata H, Nakata S. Traveling waves propagating through coupled microbeads in the Belousov-Zhabotinsky reaction. Phys Chem Chem Phys 2021; 23:24175-24179. [PMID: 34673865 DOI: 10.1039/d1cp03916d] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Spatio-temporal patterns, namely global oscillations (GO) and traveling waves (TW), were investigated in spherical microbeads loaded with a catalyst for the Belousov-Zhabotinsky (BZ) reaction onto the surface (2D-loaded) or the entire volume of the bead (3D-loaded). GO and TW selectively appeared in the 2D- and 3D-loaded beads, respectively, placed on a polyethylene terephthalate (PET) sheet in the catalyst-free BZ solution. We examined two types of coupling of the two beads: 2D-3D and 3D-3D couplings. In both cases, synchronization occurred when the minimum distance between the two beads, l, was shorter than the threshold. Herein, we reported not only temporal information, that is, phase difference, but also spatial information, that is, the directions of the TW propagating through the coupled BZ beads. In the synchronization for the 2D-3D coupling, TW in the 3D-loaded bead were initiated from the point near the 2D-loaded bead as a pacemaker and propagated in the opposite direction. By contrast, the directions of the TW in the 3D-loaded bead changed depending on l in the synchronization for the 3D-3D coupling. These experimental results can be quantitatively reproduced by numerical calculations based on the diffusion dynamics of an activator of the BZ reaction. Our results suggest that the features of spatio-temporal wave propagation are indicative of the configuration of the oscillators.
Collapse
Affiliation(s)
- Masakazu Kuze
- Graduate School of Integrated Sciences for Life, Hiroshima University, 1-3-1 Kagamiyama, Higashi-Hiroshima 739-8526, Japan.
| | - Hiroyuki Kitahata
- Graduate School of Science, Chiba University, Yayoi-cho 1-33, Inage-ku, Chiba 263-8522, Japan
| | - Satoshi Nakata
- Graduate School of Integrated Sciences for Life, Hiroshima University, 1-3-1 Kagamiyama, Higashi-Hiroshima 739-8526, Japan.
| |
Collapse
|
3
|
McGuire MK, Fuller CA, Lindner JF, Manz N. Geographic tongue as a reaction-diffusion system. CHAOS (WOODBURY, N.Y.) 2021; 31:033118. [PMID: 33810753 DOI: 10.1063/5.0020906] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/05/2020] [Accepted: 11/04/2020] [Indexed: 06/12/2023]
Abstract
Geographic tongue or benign migratory glossitis is a condition of an unknown cause characterized by chronic lesions that slowly migrate across the surface of the tongue. The condition's characteristic wavefronts suggest that it can be modeled as a reaction-diffusion system. Here, we present a model for geographic tongue pattern evolution using reaction-diffusion equations applied to portions of spheroids and paraboloids that approximate a tongue shape. We demonstrate that the observed patterns of geographic tongue lesions can be explained by propagating reaction-diffusion waves on these variably curved surfaces.
Collapse
Affiliation(s)
| | - Chase A Fuller
- Department of Physics, The College of Wooster, Wooster, Ohio 44691, USA
| | - John F Lindner
- Department of Physics, The College of Wooster, Wooster, Ohio 44691, USA
| | - Niklas Manz
- Department of Physics, The College of Wooster, Wooster, Ohio 44691, USA
| |
Collapse
|
4
|
Kuze M, Hiranishi Y, Okamoto Y, Shioi A, Nakata S. Coupling of Two Microbeads Exhibiting Different Features of Oscillations in the Belousov-Zhabotinsky Reaction. CHEM LETT 2019. [DOI: 10.1246/cl.190292] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Affiliation(s)
- Masakazu Kuze
- Graduate School of Integrated Science for Life, Hiroshima University, 1-3-1 Kagamiyama, Higashi-Hiroshima, Hiroshima 739-8526, Japan
| | - Yuri Hiranishi
- Faculty of Science, Hiroshima University, 1-3-1 Kagamiyama, Higashi-Hiroshima, Hiroshima 739-8526, Japan
| | - Yasunao Okamoto
- Department of Chemical Engineering and Materials Science, Doshisha University, Kyoto 610-0321, Japan
| | - Akihisa Shioi
- Department of Chemical Engineering and Materials Science, Doshisha University, Kyoto 610-0321, Japan
| | - Satoshi Nakata
- Graduate School of Integrated Science for Life, Hiroshima University, 1-3-1 Kagamiyama, Higashi-Hiroshima, Hiroshima 739-8526, Japan
| |
Collapse
|
5
|
Marino F, Giacomelli G. Excitable Wave Patterns in Temporal Systems with Two Long Delays and their Observation in a Semiconductor Laser Experiment. PHYSICAL REVIEW LETTERS 2019; 122:174102. [PMID: 31107096 DOI: 10.1103/physrevlett.122.174102] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/21/2018] [Indexed: 06/09/2023]
Abstract
Excitable waves arise in many spatially extended systems of either a biological, chemical, or physical nature due to the interplay between local reaction and diffusion processes. Here we demonstrate that similar phenomena are encoded in the time dynamics of an excitable system with two, hierarchically long delays. The transition from 1D localized structures to curved wave segments is experimentally observed in an excitable semiconductor laser with two feedback loops and reproduced by numerical simulations of a prototypical model. While closely related to those found in 2D excitable media, wave patterns in delayed systems exhibit unobserved features originating from causality related constraints. An appropriate dynamical representation of the data uncovers these phenomena and permits us to interpret them as the result of an effective 2D advection-reaction-diffusion process.
Collapse
Affiliation(s)
- Francesco Marino
- CNR-Istituto Nazionale di Ottica, largo E. Fermi 6, I-50125 Firenze, Italy
| | - Giovanni Giacomelli
- CNR-Istituto dei Sistemi Complessi, via Madonna del Piano 10, I-50019 Sesto Fiorentino, Italy
| |
Collapse
|
6
|
Markovic VM, Bánsági T, McKenzie D, Mai A, Pojman JA, Taylor AF. Influence of reaction-induced convection on quorum sensing in enzyme-loaded agarose beads. CHAOS (WOODBURY, N.Y.) 2019; 29:033130. [PMID: 30927847 DOI: 10.1063/1.5089295] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/18/2019] [Accepted: 03/06/2019] [Indexed: 06/09/2023]
Abstract
In theory, groups of enzyme-loaded particles producing an acid or base may show complex behavior including dynamical quorum sensing, the appearance of synchronized oscillations above a critical number or density of particles. Here, experiments were performed with the enzyme urease loaded into mm-sized agarose beads and placed in a solution of urea, resulting in an increase in pH. This behavior was found to be dependent upon the number of beads present in the array; however, reaction-induced convection occurred and plumes of high pH developed that extended to the walls of the reactor. The convection resulted in the motion of the mm-sized particles and conversion of the solution to high pH. Simulations in a simple model of the beads demonstrated the suppression of dynamical quorum sensing in the presence of flow.
Collapse
Affiliation(s)
- Vladimir M Markovic
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16, 11158 Belgrade, Serbia
| | - Tamás Bánsági
- Department of Chemical and Biological Engineering, The University of Sheffield, Mappin Street, Sheffield S1 3JD, United Kingdom
| | - Dennel McKenzie
- Department of Chemistry, Louisiana State University, 232 Choppin Hall, Baton Rouge, Louisiana 70803, USA
| | - Anthony Mai
- Department of Chemistry, Louisiana State University, 232 Choppin Hall, Baton Rouge, Louisiana 70803, USA
| | - John A Pojman
- Department of Chemistry, Louisiana State University, 232 Choppin Hall, Baton Rouge, Louisiana 70803, USA
| | - Annette F Taylor
- Department of Chemical and Biological Engineering, The University of Sheffield, Mappin Street, Sheffield S1 3JD, United Kingdom
| |
Collapse
|
7
|
Nevins TD, Kelley DH. Front tracking velocimetry in advection-reaction-diffusion systems. CHAOS (WOODBURY, N.Y.) 2018; 28:043122. [PMID: 31906630 DOI: 10.1063/1.5020055] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
In advection-reaction-diffusion systems, the spreading of a reactive scalar can be significantly influenced by the flow field in which it grows. In systems with sharp boundaries between reacted and unreacted regions, motion of the reaction fronts that lie at those boundaries can quantify spreading. Here, we present an algorithm for measuring the velocity of reaction fronts in the presence of flow, expanding previous work on tracking reaction fronts without flow. The algorithm provides localized measurements of front speed and can distinguish its two components: one from chemical dynamics and another from the underlying flow. We validate that the algorithm returns the expected front velocity components in two simulations and then show that in complex experimental flows, the measured front velocity maps fronts from one time step to the next self-consistently. Finally, we observe a variation of the chemical speed with flow speed in a variety of experiments with different time scales and length scales.
Collapse
Affiliation(s)
- Thomas D Nevins
- Department of Physics and Astronomy, University of Rochester, Rochester, New York 14627, USA
| | - Douglas H Kelley
- Department of Mechanical Engineering, University of Rochester, Rochester, New York 14627, USA
| |
Collapse
|
8
|
Nevins TD, Kelley DH. Front tracking for quantifying advection-reaction-diffusion. CHAOS (WOODBURY, N.Y.) 2017; 27:043105. [PMID: 28456164 DOI: 10.1063/1.4979668] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
We present an algorithm for measuring the speed and thickness of reaction fronts, and from those quantities, the diffusivity and the reaction rate of the active chemical species. This front-tracking algorithm provides local measurements suitable for statistics and requires only a sequence of concentration fields. Though our eventual goal is front tracking in advection-reaction-diffusion, here we demonstrate the algorithm in reaction-diffusion. We test the algorithm with validation data in which front speed and thickness are prescribed, as well as simulation results in which diffusivity and reaction rate are prescribed. In all tests, measurements closely match true values. We apply the algorithm to laboratory experiments using the Belousov-Zhabotinsky reaction, producing speed, diffusivity, and reaction rate measurements that are statistically more robust than in prior studies. Finally, we use thickness measurements to quantify the concentration profile of chemical waves in the reaction.
Collapse
Affiliation(s)
- Thomas D Nevins
- Department of Physics and Astronomy, University of Rochester, Rochester, New York 14627, USA
| | - Douglas H Kelley
- Department of Mechanical Engineering, University of Rochester, Rochester, New York 14627, USA
| |
Collapse
|
9
|
Kneer F, Obermayer K, Dahlem MA. Analyzing critical propagation in a reaction-diffusion-advection model using unstable slow waves. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2015; 38:95. [PMID: 25704900 DOI: 10.1140/epje/i2015-15010-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/05/2014] [Revised: 11/18/2014] [Accepted: 01/21/2015] [Indexed: 06/04/2023]
Abstract
The effect of advection on the propagation and in particular on the critical minimal speed of traveling waves in a reaction-diffusion model is studied. Previous theoretical studies estimated this effect on the velocity of stable fast waves and predicted the existence of a critical advection strength below which propagating waves are not supported anymore. In this paper, an analytical expression for the advection-velocity relation of the unstable slow wave is derived. In addition, the critical advection strength is calculated taking into account the unstable slow wave solution. We also analyze a two-variable reaction-diffusion-advection model numerically in a wide parameter range. Due to the new control parameter (advection) we can find stable wave propagation in the otherwise non-excitable parameter regime, if the advection strength exceeds a critical value. Comparing theoretical predictions to numerical results, we find that they are in good agreement. Theory provides an explanation for the observed behaviour.
Collapse
Affiliation(s)
- Frederike Kneer
- Department of Software Engineering and Theoretical Computer Science, Technische Universität Berlin, Ernst-Reuter-Platz 7, D-10587, Berlin, Germany,
| | | | | |
Collapse
|
10
|
Tissue voltage discordance during tachycardia versus sinus rhythm: implications for catheter ablation. Heart Rhythm 2013; 10:800-4. [PMID: 23434619 DOI: 10.1016/j.hrthm.2013.02.020] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 12/11/2012] [Indexed: 11/23/2022]
Abstract
BACKGROUND Electroanatomic mapping systems are an important tool to identify cardiac chamber voltage and assess channels of slow conduction. OBJECTIVE To assess the correlation between electroanatomic mapping voltage maps obtained during macroreentrant tachycardia compared to sinus rhythm (SR) with a contact mapping system. METHODS We retrospectively evaluated patients with atrial flutter (AFL) referred for radiofrequency ablation with electroanatomic voltage maps obtained during AFL and SR. The atrium was divided into predetermined segments. Overall atrial and segmental peak-to-peak bipolar voltages in AFL and SR were assessed. To directly compare a region within the same patient, tissue voltage differences during AFL and SR were assessed on the basis of mean voltage difference. RESULTS Sixteen patients (87% men) had available voltage maps. Eighty-one percent had typical cavotricuspid isthmus-dependent right AFL. A mean of 441.7±153.9 vs 398.1±125.4 total points (P = .22) were sampled during AFL and SR, with a mean of 99.5±58.9 vs 91.2±60.4 points (P = .45) sampled per region. Overall right atrial mean voltage was significantly higher during AFL than SR (0.554±0.092mV vs 0.473±0.079mV; P≤.001), with the lateral wall (0.707±0.120mV vs 0.573±0.097mV; P = .0004) and the cavotricuspid isthmus (0.559±0.100mV vs 0.356±0.066mV; P<.0001) also showing higher mean voltage during AFL. When compared within an individual patient, 19% (14 of 75) of the patient regions had a>0.5mV mean voltage difference and 40% (30 of 75) had a>0.25mV mean voltage difference. CONCLUSIONS These data suggest that voltage maps performed during macroreentrant atrial arrhythmias often vary significantly from maps obtained during SR.
Collapse
|
11
|
Ertl G. Oscillatory kinetics and spatio-temporal self-organization in reactions at solid surfaces. Science 2010; 254:1750-5. [PMID: 17829239 DOI: 10.1126/science.254.5039.1750] [Citation(s) in RCA: 349] [Impact Index Per Article: 24.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2022]
Abstract
Chemical reactions far from equilibrium on solid surfaces may exhibit typical phenomena of nonlinear dynamics, as exemplified by the catalytic oxidation of carbon monoxide on a platinum(110) single-crystal surface. Depending on the external parameters (temperature and partial pressures of the reactants), the temporal variation of the reaction rate may become oscillatory or even chaotic. In a parallel way, the concentration distributions of the adsorbed species on the surface form spatio-temporal patterns including propagating and standing waves, rotating spirals, as well as irregular and rapidly changing structures denoted "chemical turbulence."
Collapse
|
12
|
Foerster P, Müller SC, Hess B. Critical size and curvature of wave formation in an excitable chemical medium. Proc Natl Acad Sci U S A 2010; 86:6831-4. [PMID: 16594068 PMCID: PMC297943 DOI: 10.1073/pnas.86.18.6831] [Citation(s) in RCA: 60] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The critical radius for the outward propagation of waves in an excitable solution of the Belousov-Zhabotinskii reaction was experimentally analyzed and found to be approximately 20 mum, being in a range predicted by theory. Thus, the wave initiation depends on the critical radius in an all-or-none fashion. For waves having high positive curvature of wave fronts, a linear relationship between the curvature and their normal velocity was established, allowing computation of a diffusion coefficient of 1.9 x 10(-5) cm(2)/s for the autocatalytic species, which agrees well with results previously obtained for negatively curved wave fronts. The analysis of the dispersion of wave velocity yielded the decrease of wave velocity for small initiation periods as predicted theoretically.
Collapse
Affiliation(s)
- P Foerster
- Max-Planck-Institut für Ernährungsphysiologie, Rheinlanddamm 201, D-4600 Dortmund 1, Federal Republic of Germany
| | | | | |
Collapse
|
13
|
|
14
|
Bishop KJM, Fiałkowski M, Grzybowski BA. Micropatterning Chemical Oscillations: Waves, Autofocusing, and Symmetry Breaking. J Am Chem Soc 2005; 127:15943-8. [PMID: 16277538 DOI: 10.1021/ja054851o] [Citation(s) in RCA: 17] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Arrays of chemical oscillators are micropatterned by Wet Stamping. The technique is used to demonstrate that chemical waves can be initiated and controlled in oscillatory systems and that such waves can give rise to phenomena not observed in excitable media. Interoscillator coupling and synchronization, kinetic autofocusing, and twist-symmetry breaking are a consequence of the dependence of the oscillation phase on the local concentrations of reagents and on systems' geometry. Conditions under which a generic oscillatory system would exhibit such behaviors are determined.
Collapse
Affiliation(s)
- Kyle J M Bishop
- Department of Chemical and Biological Engineering and Northwestern Institute for Complexity, Northwestern University, 2145 Sheridan Road, Illinois 60208, USA
| | | | | |
Collapse
|
15
|
Kuklik P, Zebrowski JJ. Reentry wave formation in excitable media with stochastically generated inhomogeneities. CHAOS (WOODBURY, N.Y.) 2005; 15:33301. [PMID: 16252987 DOI: 10.1063/1.1947427] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/05/2023]
Abstract
Clinical research shows that the frequency of arrhythmia events depends on the number and area of the border zones of infarct scars. We investigate the possibility that arrhythmia is initiated by reentry waves generated by the inhomogeneity of conduction velocity at the border zone. The interaction of a plane wave with a spatially extended inhomogeneity is simulated in the FitzHugh- Nagumo model. The inhomogeneity is introduced into the model by modifying the spatial dependence of the diffusion coefficient in a stochastic manner. This results in a rich variety of spatial distributions of conductivity. A plane wave propagating through such a system may break up on the regions with low conductivity and produce numerous spiral waves. The frequency of reentry wave formation is studied as a function of the parameters of the inhomogeneity generation algorithm. Three main scenarios of reentry wave formation were found: unidirectional block, main wave-wavelet collision, and wave break up during collision, on a region in which a conduction velocity gradient occurs. These scenarios are likely candidates for the mechanisms of arrhythmia initiation in a damaged tissue, e.g., the border zone of an infarct scar.
Collapse
Affiliation(s)
- Paweł Kuklik
- Faculty of Physics and Center of Excellence for Complex Systems Research at Warsaw University of Technology, ul. Koszykowa 75, Warsaw, Poland.
| | | |
Collapse
|
16
|
Marts B, Martinez K, Lin AL. Front dynamics in an oscillatory bistable Belousov-Zhabotinsky chemical reaction. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2004; 70:056223. [PMID: 15600744 DOI: 10.1103/physreve.70.056223] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/30/2004] [Indexed: 05/24/2023]
Abstract
We observe breathing front dynamics which select three distinct types of bistable patterns in the 2:1 resonance regime of the periodically forced oscillatory Belousov-Zhabotinsky reaction. We measure the curvature-driven shrinking of a circular domain R approximately t(1/2) at forcing frequencies below a specific value, and show that the fast time scale front oscillations (breathing) drive this slow time scale shrinking. Above a specific frequency, we observe fronts of higher curvature grow instead of shrink and labyrinth patterns form. Just below the transition frequency is a relatively narrow range of frequencies where the curvature-driven coarsening is balanced by a competing front interaction, which leads to a pattern of localized structures. The length scale of the localized structure and labyrinth patterns is set by the front interactions.
Collapse
Affiliation(s)
- Bradley Marts
- Center for Nonlinear and Complex Systems and Department of Physics, Duke University, Durham, North Carolina 27708, USA
| | | | | |
Collapse
|
17
|
Kirner T, Ackermann J. Pattern selection in a cooperative biochemical in vitro amplification system: the role of parasites. J Theor Biol 2003; 224:539-49. [PMID: 12957125 DOI: 10.1016/s0022-5193(03)00202-9] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Previously, numerical simulations have shown that evolving systems can be stabilized against emerging parasites by pattern formation in spatially extended flow reactors. Hence, it can be argued that pattern formation is a prerequisite for any experimental investigation of the biochemical evolution of cooperative function. Here, we study a model of an experimental biochemical system for the cooperative in vitro amplification of DNA strands and show that emerging parasites can induce a complex pattern formation even when no pattern formation occurs without parasites. In an adiabatic approximation where the cooperative amplification reaction is assumed to adapt fast to slowly emerging parasites, the parasite concentration itself acts as a Steuer parameter for the selection of various complex patterns. Without such an adiabatic approximation only transient patterns emerge. As any species can grow for very low concentrations, the parasite is able to infect the entire reactor and the system is finally diluted out. In the experimental biochemical system, however, the species are individual molecules and the growth of spatially separated, non-infected regions becomes feasible. Hence a cutoff threshold for the minimal concentration is applied. In these simulations the otherwise lethal infection by parasites induces the formation of spatiotemporal spirals, and this spatial structure help the host and parasitoid species to survive together. These theoretical results describe an inherent property of cooperative reactions and have an important impact on experimental investigations on the molecular evolution and complex function in spatially extended reactors. Since the formation of the complex pattern is restricted either to a rather large cutoff value or a special choice of the kinetic parameters, we, however, conclude that the persistence of evolving cooperative amplification is not possible in a simple reaction-diffusion reactor. Experimental set-ups with patchy environments, e.g. biomolecular amplification in coupled microstructured flow chambers or in microemulsion, are eligible candidates for the observation of such a self-organized pattern selection.
Collapse
Affiliation(s)
- T Kirner
- Fraunhofer Gesellschaft, Biomolecular Information Processing, Schloss Birlinghoven, D-53754 Sankt Augustin, Germany
| | | |
Collapse
|
18
|
Bub G, Shrier A. Propagation through heterogeneous substrates in simple excitable media models. CHAOS (WOODBURY, N.Y.) 2002; 12:747-753. [PMID: 12779603 DOI: 10.1063/1.1502481] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/24/2023]
Abstract
The interaction of waves and obstacles is simulated by adding heterogeneities to a FitzHugh-Nagumo model and a cellular automata model. The cellular automata model is formulated to account for heterogeneities by modelling the interaction between current sources and current sinks. In both models, wave fronts propagate if the size of the heterogeneities is small, and block if the size of the heterogeneities is large. For intermediate values, wave fronts break up into numerous spiral waves. The theoretical models give insights concerning spiral wave formation in heterogeneous excitable media. (c) 2002 American Institute of Physics.
Collapse
Affiliation(s)
- Gil Bub
- Center for Nonlinear Dynamics in Physiology and Medicine, Department of Physiology, McIntyer Medical Sciences Building, McGill University, Montreal, Quebec, Canada H3G 1Y6
| | | |
Collapse
|
19
|
Wussling MH, Krannich K, Drygalla V, Podhaisky H. Calcium waves in agarose gel with cell organelles: implications of the velocity curvature relationship. Biophys J 2001; 80:2658-66. [PMID: 11371442 PMCID: PMC1301453 DOI: 10.1016/s0006-3495(01)76235-2] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022] Open
Abstract
Calcium oscillations and waves have been observed not only in several types of living cells but also in less complex systems of isolated cell organelles. Here we report the determination of apparent Ca2+ diffusion coefficients in a novel excitable medium of agarose gel with homogeneously distributed vesicles of skeletal sarcoplasmic reticulum. Spatiotemporal calcium patterns were visualized by confocal laser scanning fluorescence microscopy. To obtain characteristic parameters of the velocity curvature relationship, namely, apparent diffusion coefficient, velocity of plane calcium waves, and critical radius, positively and negatively curved wave fronts were analyzed. It is demonstrated that gel-immobilized cell organelles reveal features of an excitable medium. Apparent Ca2+ diffusion coefficients of the in vitro system, both in the absence or in the presence of mitochondria, were found to be higher than in cardiac myocytes and lower than in unbuffered agarose gel. Plane calcium waves propagated markedly slower in the in vitro system than in rat cardiac myocytes. Whereas mitochondria significantly reduced the apparent Ca2+ diffusion coefficient of the in vitro system, propagation velocity and critical size of calcium waves were found to be nearly unchanged. These results suggest that calcium wave propagation depends on the kinetics of calcium release rather than on diffusion.
Collapse
Affiliation(s)
- M H Wussling
- Julius Bernstein Institute of Physiology, Martin Luther University Halle-Wittenberg, D-06097 Halle, Germany
| | | | | | | |
Collapse
|
20
|
Danse PW, Garratt CJ, Mast F, Allessie MA. Preferential depression of conduction around a pivot point in rabbit ventricular myocardium by potassium and flecainide. J Cardiovasc Electrophysiol 2000; 11:262-73. [PMID: 10749349 DOI: 10.1111/j.1540-8167.2000.tb01795.x] [Citation(s) in RCA: 17] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
INTRODUCTION During reentrant arrhythmias, the circulating wavefront often makes a sharp turn around a functional or anatomic barrier. We tested the hypothesis that lowering the safety factor for conduction by high K+ or flecainide preferentially depresses conduction of sharply turning wavefronts. METHODS AND RESULTS In 16 Langendorff-perfused rabbit hearts, a thin layer of anisotropic ventricular myocardium was made using a cryoprocedure. In this layer, a linear radiofrequency lesion was made parallel to the fiber orientation. The tip of the lesion was extended by a short incision. U-turning wavefronts were initiated by pacing at one side of the lesion. A mapping electrode (240 electrodes, resolution 350 to 700 microm) was used to measure conduction times and velocity of planar waves (longitudinal and transverse) and U-turning wavefronts. The safety factor for conduction was lowered by high potassium (8, 10, and 12 mmol/L) and flecainide (1 and 2 mg/L). On average, high potassium and flecainide increased the conduction times of U-turning wavefronts 1.6 times more than longitudinal or transverse planar wavefronts (P < 0.01). At a critical lowering of the excitatory current, functional conduction block occurred at the pivot point, which forced the wavefront to make a longer U-turn. In these cases, the total U-turn conduction time increased from 27+/-9 msec to 75+/-37 msec. About 40% of this delay was caused by a shift of the pivot point and consequent lengthening of the returning pathway. CONCLUSION Lowering the amount of excitatory current by potassium or flecainide preferentially impairs U-turn conduction. The occurrence of long delays and conduction block at pivot points may explain the mode of action of Class I drugs.
Collapse
Affiliation(s)
- P W Danse
- Department of Physiology, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University, The Netherlands
| | | | | | | |
Collapse
|
21
|
Shoji S, Adaniya H, Numano F, Hiraoka M. Factors determining clockwise and counterclockwise conduction patterns in atrial reentrant tachycardias: a rabbit model of atrial flutter. J Cardiovasc Electrophysiol 2000; 11:311-23. [PMID: 10749354 DOI: 10.1111/j.1540-8167.2000.tb01800.x] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Abstract
INTRODUCTION In the development of atrial flutter due to reentry, the crista terminalis is supposed to pose a conduction barrier, but the role of its longitudinal conduction in determining the propagation pattern of the reentrant impulse is not known. In rabbit right atrial preparations, we induced reentrant atrial tachycardias and examined the effects of transverse section of the crista terminalis on the development and conduction patterns of arrhythmias. METHODS AND RESULTS Right atrial preparations from 12 albino rabbits were placed endocardial surface down in a chamber with an array of 48 bipolar electrodes to draw activation maps. A single premature stimulus was delivered to induce tachycardias at the free wall. In the control, five instances of tachycardia per preparation were induced and another five were induced after cutting the crista terminalis. In the control, the mean duration of tachycardia was 127.1+/-25.2 seconds. The tachycardia was counterclockwise in 39 of 60 instances, clockwise in 12, and undetermined in 4 defined as "atypical." After transverse section of the crista terminalis, the duration was prolonged to 372.6+/-30.4 seconds, but the conduction patterns were not changed. In the free wall, counterclockwise reentry had a broader wavefront and faster conduction than clockwise reentry. CONCLUSION Longitudinal conduction block at the crista terminalis contributed to maintenance of reentrant atrial tachycardias, but had no influence on their propagation patterns. Clockwise and counterclockwise rotation of impulses in reentrant tachycardias had different paths and velocities of the wavefront in the free wall of the right atrium.
Collapse
Affiliation(s)
- S Shoji
- Department of Cardiovascular Diseases, Medical Research Institute, Tokyo Medical and Dental University, Japan
| | | | | | | |
Collapse
|
22
|
Cabo C, Pertsov AM, Davidenko JM, Jalife J. Electrical turbulence as a result of the critical curvature for propagation in cardiac tissue. CHAOS (WOODBURY, N.Y.) 1998; 8:116-126. [PMID: 12779715 DOI: 10.1063/1.166292] [Citation(s) in RCA: 17] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/24/2023]
Abstract
In cardiac tissue, the propagation of electrical excitation waves is dependent on the active properties of the cell membrane (ionic channels) and the passive electrical properties of cardiac tissue (passive membrane properties, distribution of gap junctions, and cell shapes). Initiation of cardiac arrhythmias is usually associated with heterogeneities in the active and/or passive properties of cardiac tissue. However, as a result of the effect of wave front geometry (curvature) on propagation of cardiac waves, inexcitable anatomical obstacles, like veins and arteries, may cause the formation of self-sustained vortices and uncontrolled high-frequency excitation in normal homogeneous myocardium. (c) 1998 American Institute of Physics.
Collapse
Affiliation(s)
- Candido Cabo
- Department of Pharmacology, College of Physicians and Surgeons of Columbia University, New York, New York 10032
| | | | | | | |
Collapse
|
23
|
Wussling MH, Scheufler K, Schmerling S, Drygalla V. Velocity-curvature relationship of colliding spherical calcium waves in rat cardiac myocytes. Biophys J 1997; 73:1232-42. [PMID: 9284291 PMCID: PMC1181023 DOI: 10.1016/s0006-3495(97)78156-6] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Abstract
Colliding spherical calcium waves in enzymatically isolated rat cardiac myocytes develop new wavefronts propagating perpendicular to the original direction. When investigated by confocal laser scanning microscopy (CLSM), using the fluorescent Ca2+ indicator fluo-3 AM, "cusp"-like structures become visible that are favorably approximated by double parabolae. The time-dependent position of the vertices is used to determine propagation velocity and negative curvature of the wavefront in the region of collision. It is evident that negatively curved waves propagate faster than positively curved, single waves. Considering two perfectly equal expanding circular waves, we demonstrated that the collision of calcium waves is due to an autocatalytic process (calcium-induced calcium release), and not to a simple phenomenon of interference. Following the spatiotemporal organization in simpler chemical systems maintained under conditions far from the thermodynamic equilibrium (Belousov-Zhabotinskii reaction), the dependence of the normal velocity on the curvature of the spreading wavefront is given by a linear relation. The so-called velocity-curvature relationship makes clear that the velocity is enhanced by curvature toward the direction of forward propagation and decreased by curvature away from the direction of forward propagation (with an influence of the diffusion coefficient). Experimentally obtained velocity data of both negatively and positively curved calcium waves were approximated by orthogonal weighted regression. The negative slope of the straight line resulted in an effective diffusion coefficient of 1.2 x 10(-4) mm2/s. From the so-called critical radius, which must be exceeded to initiate a traveling calcium wave, a critical volume (with enhanced [Ca2+]i) of approximately 12 microm3 was calculated. This is almost identical to the volume that is occupied by a single calcium spark.
Collapse
Affiliation(s)
- M H Wussling
- Julius Bernstein Institute of Physiology, Martin Luther University, Halle-Wittenberg, Halle, Germany.
| | | | | | | |
Collapse
|
24
|
Milton RA, Scott SK. Wave selection and lateral instabilities in mixed autocatalytic systems. J Chem Phys 1997. [DOI: 10.1063/1.473627] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
25
|
Diffraction of spherical waves by a toroidal obstacle: eikonal approach to excitable reaction—diffusion systems. Proc Math Phys Eng Sci 1996. [DOI: 10.1098/rspa.1996.0147] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
|
26
|
Ševčíková H, Schreiber I, Marek M. Dynamics of Oxidation Belousov−Zhabotinsky Waves in an Electric Field. ACTA ACUST UNITED AC 1996. [DOI: 10.1021/jp961789w] [Citation(s) in RCA: 34] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Hana Ševčíková
- Department of Chemical Engineering, Prague Institute of Chemical Technology, Technická 5, 166 28 Prague 6, Czech Republic
| | - Igor Schreiber
- Department of Chemical Engineering, Prague Institute of Chemical Technology, Technická 5, 166 28 Prague 6, Czech Republic
| | - Miloš Marek
- Department of Chemical Engineering, Prague Institute of Chemical Technology, Technická 5, 166 28 Prague 6, Czech Republic
| |
Collapse
|
27
|
Brazhnik PK, Tyson JJ. Nonspiral excitation waves beyond the eikonal approximation. PHYSICAL REVIEW. E, STATISTICAL PHYSICS, PLASMAS, FLUIDS, AND RELATED INTERDISCIPLINARY TOPICS 1996; 54:4338-4346. [PMID: 9965583 DOI: 10.1103/physreve.54.4338] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/12/2023]
|
28
|
Affiliation(s)
- B M Sager
- Rowland Institute for Science, Cambridge, Massachusetts 02142, USA
| |
Collapse
|
29
|
Hwang SC, Halpin-Healy T. Chemical wave refraction phenomena. PHYSICAL REVIEW. E, STATISTICAL PHYSICS, PLASMAS, FLUIDS, AND RELATED INTERDISCIPLINARY TOPICS 1996; 54:3009-3012. [PMID: 9965423 DOI: 10.1103/physreve.54.3009] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
|
30
|
Affiliation(s)
- K. Agladze
- National Institute of Materials and Chemical Research, 1-1 Higashi, Tsukuba, Ibaraki 305, Japan
| | - R. R. Aliev
- National Institute of Materials and Chemical Research, 1-1 Higashi, Tsukuba, Ibaraki 305, Japan
| | - T. Yamaguchi
- National Institute of Materials and Chemical Research, 1-1 Higashi, Tsukuba, Ibaraki 305, Japan
| | - K. Yoshikawa
- Graduate School for Human Informatics, Nagoya University, Nagoya, 464-01, Japan
| |
Collapse
|
31
|
Brazhnik PK, Tyson JJ. Propagation of waves through a line of discontinuity in two-dimensional excitable media: Refraction and reflection of autowaves. PHYSICAL REVIEW. E, STATISTICAL PHYSICS, PLASMAS, FLUIDS, AND RELATED INTERDISCIPLINARY TOPICS 1996; 54:1958-1968. [PMID: 9965279 DOI: 10.1103/physreve.54.1958] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
|
32
|
Steinbock O, Kettunen P. Chemical clocks on the basis of rotating waves. Measuring irrational numbers from period ratios. Chem Phys Lett 1996. [DOI: 10.1016/0009-2614(96)00099-1] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
33
|
Cabo C, Pertsov AM, Davidenko JM, Baxter WT, Gray RA, Jalife J. Vortex shedding as a precursor of turbulent electrical activity in cardiac muscle. Biophys J 1996; 70:1105-11. [PMID: 8785270 PMCID: PMC1225040 DOI: 10.1016/s0006-3495(96)79691-1] [Citation(s) in RCA: 120] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023] Open
Abstract
In cardiac tissue, during partial blockade of the membrane sodium channels, or at high frequencies of excitation, inexcitable obstacles with sharp edges may destabilize the propagation of electrical excitation waves, causing the formation of self-sustained vortices and turbulent cardiac electrical activity. The formation of such vortices, which visually resembles vortex shedding in hydrodynamic turbulent flows, was observed in sheep epicardial tissue using voltage-sensitive dyes in combination with video-imaging techniques. Vortex shedding is a potential mechanism leading to the spontaneous initiation of uncontrolled high-frequency excitation of the heart.
Collapse
Affiliation(s)
- C Cabo
- Department of Pharmacology, SUNY Health Science Center at Syracuse 13210, USA
| | | | | | | | | | | |
Collapse
|
34
|
Ševčíková H, Kosek J, Marek M. Splitting of 2D Waves of Excitation in a Direct Current Electric Field. ACTA ACUST UNITED AC 1996. [DOI: 10.1021/jp952004+] [Citation(s) in RCA: 24] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
35
|
Epstein IR, Showalter K. Nonlinear Chemical Dynamics: Oscillations, Patterns, and Chaos. ACTA ACUST UNITED AC 1996. [DOI: 10.1021/jp953547m] [Citation(s) in RCA: 391] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
36
|
Ong JJ, Cha YM, Kriett JM, Boyce K, Feld GK, Chen PS. The relation between atrial fibrillation wavefront characteristics and accessory pathway conduction. J Clin Invest 1995; 96:2284-96. [PMID: 7593615 PMCID: PMC185879 DOI: 10.1172/jci118284] [Citation(s) in RCA: 21] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Abstract
Although the source-sink relationship for impulse propagation in cardiac tissues has been demonstrated in vitro, there has been no verification of this hypothesis in humans. Accordingly, eight patients undergoing surgical division of their accessory pathways were studied. A 56-channel (7 x 8) bipolar plaque electrode array was placed over the atrioventricular groove on the accessory pathway and atrial fibrillation electrically induced. 10 episodes of QRS transition from consecutively preexcited to nonpreexcited complexes were analyzed. This showed that consecutively preexcited QRS complexes were always associated with uniform large atrial wavefronts. Immediately prior to QRS transition, four general types of changes were observed: (a) premature invasion by secondary wavefronts creating local conduction block (n = 5); (b) wavefront collision leading to wavefront curvature (n = 2); (c) transition from a uniform large atrial wavefront to multiple fractionated small wavefronts (n = 1); and (d) uniform atrial wavefronts "marching" into the accessory pathway refractory period (n = 2). We conclude that local atrial wavefront characteristics are important factors influencing impulse propagation through the accessory pathway. The findings that local wavefront collision, curvature, or fractionation often precede loss of accessory pathway conduction support the notion that source-sink relationship is an important determinant of the safety factor for impulse propagation in the human heart.
Collapse
Affiliation(s)
- J J Ong
- Department of Medicine, Cedars-Sinai Medical Center, Los Angeles, California 90048, USA
| | | | | | | | | | | |
Collapse
|
37
|
Knisley SB, Hill BC. Effects of bipolar point and line stimulation in anisotropic rabbit epicardium: assessment of the critical radius of curvature for longitudinal block. IEEE Trans Biomed Eng 1995; 42:957-66. [PMID: 8582725 DOI: 10.1109/10.464369] [Citation(s) in RCA: 47] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
Excitation front shape and velocity were studied in anisotropic perfused rabbit epicardium stained with potentiometric fluorescent dye. In the combined results from all experiments, convex excitation fronts produced by stimulation with a single electrode propagated longitudinally 13.3% slower than flat excitation fronts produced by stimulation with a line of electrodes. For transverse propagation, the two stimulation methods produced similar flat excitation fronts and velocities. The critical excitation front radius of curvature for longitudinal block (Rcr), calculated from excitable media theory, was 92 microns in control hearts. In hearts exposed to diacetyl monoxime (20 mmol/L), which decreases inward sodium current, Rcr was 175 microns. The slower longitudinal propagation velocity of convex fronts versus flat fronts and the theoretically predicted critical radius of curvature may be important for propagation and block of ectopic depolarizations in the heart.
Collapse
Affiliation(s)
- S B Knisley
- Division of Cardiovascular Disease, School of Medicine, University of Alabama, Birmingham 35294-0019, USA
| | | |
Collapse
|
38
|
Yoneyama M. Fluorescence microscopic studies of spatial structures in Ru(bpy)32+ monolayers. Colloids Surf A Physicochem Eng Asp 1995. [DOI: 10.1016/0927-7757(95)03249-d] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
39
|
Graham MD, Bär M, Kevrekidis IG, Asakura K, Lauterbach J, Rotermund H, Ertl G. Catalysis on microstructured surfaces: Pattern formation during CO oxidation in complex Pt domains. PHYSICAL REVIEW. E, STATISTICAL PHYSICS, PLASMAS, FLUIDS, AND RELATED INTERDISCIPLINARY TOPICS 1995; 52:76-93. [PMID: 9963407 DOI: 10.1103/physreve.52.76] [Citation(s) in RCA: 29] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
|
40
|
Lee KJ, Swinney HL. Lamellar structures and self-replicating spots in a reaction-diffusion system. PHYSICAL REVIEW. E, STATISTICAL PHYSICS, PLASMAS, FLUIDS, AND RELATED INTERDISCIPLINARY TOPICS 1995; 51:1899-1915. [PMID: 9962849 DOI: 10.1103/physreve.51.1899] [Citation(s) in RCA: 39] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
|
41
|
Cabo C, Pertsov AM, Baxter WT, Davidenko JM, Gray RA, Jalife J. Wave-front curvature as a cause of slow conduction and block in isolated cardiac muscle. Circ Res 1994; 75:1014-28. [PMID: 7525101 DOI: 10.1161/01.res.75.6.1014] [Citation(s) in RCA: 178] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
We have investigated the role of wave-front curvature on propagation by following the wave front that was diffracted through a narrow isthmus created in a two-dimensional ionic model (Luo-Rudy) of ventricular muscle and in a thin (0.5-mm) sheet of sheep ventricular epicardial muscle. The electrical activity in the experimental preparations was imaged by using a high-resolution video camera that monitored the changes in fluorescence of the potentiometric dye di-4-ANEPPS on the surface of the tissue. Isthmuses were created both parallel and perpendicular to the fiber orientation. In both numerical and biological experiments, when a planar wave front reached the isthmus, it was diffracted to an elliptical wave front whose pronounced curvature was very similar to that of a wave front initiated by point stimulation. In addition, the velocity of propagation was reduced in relation to that of the original planar wave. Furthermore, as shown by the numerical results, wave-front curvature changed as a function of the distance from the isthmus. Such changes in local curvature were accompanied by corresponding changes in velocity of propagation. In the model, the critical isthmus width was 200 microns for longitudinal propagation and 600 microns for transverse propagation of a single planar wave initiated proximal to the isthmus. In the experiments, propagation depended on the width of the isthmus for a fixed stimulation frequency. Propagation through an isthmus of fixed width was rate dependent both along and across fibers. Thus, the critical isthmus width for propagation was estimated in both directions for different frequencies of stimulation. In the longitudinal direction, for cycle lengths between 200 and 500 milliseconds, the critical width was < 1 mm; for 150 milliseconds, it was estimated to be between 1.3 and 2 mm; and for the maximum frequency of stimulation (117 +/- 15 milliseconds), it was > 2.5 mm. In the transverse direction, critical width was between 1.78 and 2.32 mm for a basic cycle length of 200 milliseconds. It increased to values between 2.46 and 3.53 mm for a basic cycle length of 150 milliseconds. The overall results demonstrate that the curvature of the wave front plays an important role in propagation in two-dimensional cardiac muscle and that changes in curvature may cause slow conduction or block.
Collapse
Affiliation(s)
- C Cabo
- Department of Pharmacology, SUNY Health Science Center, Syracuse 13210
| | | | | | | | | | | |
Collapse
|
42
|
Steinbock O, Zykov VS, Müller SC. Wave propagation in an excitable medium along a line of a velocity jump. PHYSICAL REVIEW. E, STATISTICAL PHYSICS, PLASMAS, FLUIDS, AND RELATED INTERDISCIPLINARY TOPICS 1993; 48:3295-3298. [PMID: 9960985 DOI: 10.1103/physreve.48.3295] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
|
43
|
Steinbock O, Siegert F, Müller SC, Weijer CJ. Three-dimensional waves of excitation during Dictyostelium morphogenesis. Proc Natl Acad Sci U S A 1993; 90:7332-5. [PMID: 8394018 PMCID: PMC47131 DOI: 10.1073/pnas.90.15.7332] [Citation(s) in RCA: 71] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023] Open
Abstract
Cells in Dictyostelium slugs follow well-defined patterns of motion. We found that the chemotactic cell response is controlled by a scroll wave of messenger concentration in the highly excitable prestalk zone of the slug that decays in the less-excitable prespore region into planar wave fronts. This phenomenon is investigated by numerical solutions of partial differential equations that couple local nonlinear kinetics and diffusive transport of the chemotactic signal. In the interface of both regions a complex twisted scroll wave is formed that reduces the wave frequency in the prespore zone. The spatio-temporal dynamics of waves and filaments are followed over 33 periods of rotation. These results yield an explanation of collective self-organized cell motion in a multicellular organism.
Collapse
Affiliation(s)
- O Steinbock
- Max-Planck-Institut für molekulare Physiologie, Dortmund, Federal Republic of Germany
| | | | | | | |
Collapse
|
44
|
Nettesheim S, von Oertzen A, Rotermund HH, Ertl G. Reaction diffusion patterns in the catalytic CO‐oxidation on Pt(110): Front propagation and spiral waves. J Chem Phys 1993. [DOI: 10.1063/1.464323] [Citation(s) in RCA: 238] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
45
|
Siegert F, Weijer CJ. Three-dimensional scroll waves organize Dictyostelium slugs. Proc Natl Acad Sci U S A 1992; 89:6433-7. [PMID: 1631140 PMCID: PMC49515 DOI: 10.1073/pnas.89.14.6433] [Citation(s) in RCA: 162] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
To test the hypothesis that periodic signals and chemotaxis direct later morphogenesis in Dictyostelium discoideum, we investigated cell behavior and cell movement in slugs. Trails of neutral red-stained prestalk and anterior-like cells were recorded by high-resolution digital image processing. Neutral red-stained anterior-like cells in the prespore zone of slugs move straight forward in the direction of slug migration and, furthermore, show coherent periodic cell movement. In contrast, cells in the prestalk zone move along completely different trajectories. Prestalk cells move perpendicular to the direction of slug migration; that is, they rotate around the tip. The cell movement data show that the chemotactic signal in the slug propagates as a three-dimensional scroll wave in the prestalk zone and as a planar wave in the prespore zone. The different behavior of prestalk and prespore cells is most likely caused by a difference in the oscillatory properties of the two cell types. We provide evidence that the slug stage of Dictyostelium behaves like an excitable system and that a (twisted) scroll wave organizes slug formation and migration.
Collapse
Affiliation(s)
- F Siegert
- Zoologisches Institut, Universität München, Federal Republic of Germany
| | | |
Collapse
|
46
|
Nagy‐Ungvarai Z, Ungvarai J, Müller SC, Hess B. The role of curvature and pulse width for transition to unstable wave fronts in the Belousov–Zhabotinsky reaction. J Chem Phys 1992. [DOI: 10.1063/1.463279] [Citation(s) in RCA: 17] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
47
|
Abstract
The drift of spiral waves toward regions of higher light intensity was observed experimentally in the ruthenium-catalyzed Belousov-Zhabotinsky reaction. A light gradient can thus be used to manipulate optical information in new computational systems based on photochemical media. The drift of a gradient that is rotationally invariant in space is three to four times as fast as that of a translationally invariant gradient. Simulations based on the use of a cellular automaton, which is made isotropic by a semirandom distribution of cells, are in agreement with the experimental results.
Collapse
|
48
|
Abstract
We derive an "eikonal-curvature" equation to describe the propagation of action potential wavefronts in myocardium. This equation is used to study the effects of fiber orientation on propagation in the myocardial wall. There are significant computational advantages to the use of an eikonal-curvature equation over a full ionic model of action potential spread. With this model, it is shown that the experimentally observed misalignment of spreading action potential "ellipses" from fiber orientation in level myocardial surfaces is adequately explained by the rotation of fiber orientation through the myocardial wall. Additionally, it is shown that apparently high propagation velocities on the epicardial and endocardial surfaces are the result of propagation into the midwall region and acceleration along midwall fibers before reemergence at an outer surface at a time preceding what could be accomplished with propagation along the surface alone.
Collapse
Affiliation(s)
- J P Keener
- Department of Mathematics, University of Utah, Salt Lake City 84112
| |
Collapse
|
49
|
Lechleiter J, Girard S, Peralta E, Clapham D. Spiral calcium wave propagation and annihilation in Xenopus laevis oocytes. Science 1991; 252:123-6. [PMID: 2011747 DOI: 10.1126/science.2011747] [Citation(s) in RCA: 466] [Impact Index Per Article: 14.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Intracellular calcium (Ca2+) is a ubiquitous second messenger. Information is encoded in the magnitude, frequency, and spatial organization of changes in the concentration of cytosolic free Ca2+. Regenerative spiral waves of release of free Ca2+ were observed by confocal microscopy in Xenopus laevis oocytes expressing muscarinic acetylcholine receptor subtypes. This pattern of Ca2+ activity is characteristic of an intracellular milieu that behaves as a regenerative excitable medium. The minimal critical radius for propagation of focal Ca2+ waves (10.4 micrometers) and the effective diffusion constant for the excitation signal (2.3 x 10(-6) square centimeters per second) were estimated from measurements of velocity and curvature of circular wavefronts expanding from foci. By modeling Ca2+ release with cellular automata, the absolute refractory period for Ca2+ stores (4.7 seconds) was determined. Other phenomena expected of an excitable medium, such as wave propagation of undiminished amplitude and annihilation of colliding wavefronts, were observed.
Collapse
Affiliation(s)
- J Lechleiter
- Department of Pharmacology, Mayo Foundation, Rochester, MN 55905
| | | | | | | |
Collapse
|
50
|
Propagation of Chemical Waves in Discrete Excitable Media: Anisotropic and Isotropic Wave Fronts. NONLINEAR WAVE PROCESSES IN EXCITABLE MEDIA 1991. [DOI: 10.1007/978-1-4899-3683-7_43] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/11/2023]
|