Waves in mathematical models of intracellular calcium and other excitable systems
Oscillations in cytoplasmic calcium concentration are a crucial control mechanism in almost every cell type. Two important classes of oscillation are of particular interest: solitary and periodic waves. Both types of waves are commonly observed in physical experiments and found in mathematical mode...
Main Author: | |
---|---|
Other Authors: | |
Published: |
ResearchSpace@Auckland
2012
|
Online Access: | http://hdl.handle.net/2292/14482 |
id |
ndltd-AUCKLAND-oai-researchspace.auckland.ac.nz-2292-14482 |
---|---|
record_format |
oai_dc |
spelling |
ndltd-AUCKLAND-oai-researchspace.auckland.ac.nz-2292-144822012-07-03T11:37:30ZWaves in mathematical models of intracellular calcium and other excitable systemsZhang, WenjunOscillations in cytoplasmic calcium concentration are a crucial control mechanism in almost every cell type. Two important classes of oscillation are of particular interest: solitary and periodic waves. Both types of waves are commonly observed in physical experiments and found in mathematical models of calcium dynamics and other excitable systems. In this thesis, we try to understand these two classes of wave solutions. We first investigate wave solutions of the canonical excitable model, the FitzHugh-Nagumo (FHN) equations. We analyze the FHN equations using geometric singular perturbation theory and numerical integration, and find some new codimension-two organizing centres of the overall dynamics. Many analytical results about the FHN model in its classical form have already been established. We devise a transformation to change the form of the FHN equations we study into the classical form to make use of the results. This enables us to show how basic features of the bifurcation structure of the FHN equations arise from the singular limit. We then study waves of a representative calcium model. We analyze the dynamics of the calcium model in the singular limit, and show how homoclinic and Hopf bifurcations of the full system arise as perturbations of singular homoclinic and Hopf bifurcations. We compare the wave solutions in the FHN model and the calcium model, and show that the dynamics of the two models differ in some respects (most importantly, in the way in which diffusion enters the equations). We conclude that the FHN model should not uniformly be used as a prototypical model for calcium dynamics. Motivated by phenomena seen in the FHN and calcium models, we then investigate reduction techniques for excitable systems, including the quasi-steady state approximation and geometric singular perturbation theory, and show that criticality of Hopf bifurcations may be changed when applying these reduction methods to slow-fast biophysical systems. This suggests that great care should be taken when using reduction techniques such as these, to ensure that spurious conclusions about the dynamics of a model are not drawn from the dynamics of a reduced version of the model. Finally, we describe the class of numerical algorithms used to compute features of the detailed bifurcation sets for the FHN and calcium models, and show how these were used to locate a non-structurally stable heteroclinic connection between periodic orbits in a calcium model; this is the first time such a global bifurcation has been computed.ResearchSpace@AucklandKirk, VivienSneyd, James2012-03-15T22:26:33Z2012-03-15T22:26:33Z2011Thesishttp://hdl.handle.net/2292/14482PhD Thesis - University of AucklandUoA2255444Items in ResearchSpace are protected by copyright, with all rights reserved, unless otherwise indicated. Previously published items are made available in accordance with the copyright policy of the publisher.https://researchspace.auckland.ac.nz/docs/uoa-docs/rights.htmhttp://creativecommons.org/licenses/by-nc-sa/3.0/nz/Copyright: The author |
collection |
NDLTD |
sources |
NDLTD |
description |
Oscillations in cytoplasmic calcium concentration are a crucial control mechanism in almost every cell type. Two important classes of oscillation are of particular interest: solitary and periodic waves. Both types of waves are commonly observed in physical experiments and found in mathematical models of calcium dynamics and other excitable systems. In this thesis, we try to understand these two classes of wave solutions. We first investigate wave solutions of the canonical excitable model, the FitzHugh-Nagumo (FHN) equations. We analyze the FHN equations using geometric singular perturbation theory and numerical integration, and find some new codimension-two organizing centres of the overall dynamics. Many analytical results about the FHN model in its classical form have already been established. We devise a transformation to change the form of the FHN equations we study into the classical form to make use of the results. This enables us to show how basic features of the bifurcation structure of the FHN equations arise from the singular limit. We then study waves of a representative calcium model. We analyze the dynamics of the calcium model in the singular limit, and show how homoclinic and Hopf bifurcations of the full system arise as perturbations of singular homoclinic and Hopf bifurcations. We compare the wave solutions in the FHN model and the calcium model, and show that the dynamics of the two models differ in some respects (most importantly, in the way in which diffusion enters the equations). We conclude that the FHN model should not uniformly be used as a prototypical model for calcium dynamics. Motivated by phenomena seen in the FHN and calcium models, we then investigate reduction techniques for excitable systems, including the quasi-steady state approximation and geometric singular perturbation theory, and show that criticality of Hopf bifurcations may be changed when applying these reduction methods to slow-fast biophysical systems. This suggests that great care should be taken when using reduction techniques such as these, to ensure that spurious conclusions about the dynamics of a model are not drawn from the dynamics of a reduced version of the model. Finally, we describe the class of numerical algorithms used to compute features of the detailed bifurcation sets for the FHN and calcium models, and show how these were used to locate a non-structurally stable heteroclinic connection between periodic orbits in a calcium model; this is the first time such a global bifurcation has been computed. |
author2 |
Kirk, Vivien |
author_facet |
Kirk, Vivien Zhang, Wenjun |
author |
Zhang, Wenjun |
spellingShingle |
Zhang, Wenjun Waves in mathematical models of intracellular calcium and other excitable systems |
author_sort |
Zhang, Wenjun |
title |
Waves in mathematical models of intracellular calcium and other excitable systems |
title_short |
Waves in mathematical models of intracellular calcium and other excitable systems |
title_full |
Waves in mathematical models of intracellular calcium and other excitable systems |
title_fullStr |
Waves in mathematical models of intracellular calcium and other excitable systems |
title_full_unstemmed |
Waves in mathematical models of intracellular calcium and other excitable systems |
title_sort |
waves in mathematical models of intracellular calcium and other excitable systems |
publisher |
ResearchSpace@Auckland |
publishDate |
2012 |
url |
http://hdl.handle.net/2292/14482 |
work_keys_str_mv |
AT zhangwenjun wavesinmathematicalmodelsofintracellularcalciumandotherexcitablesystems |
_version_ |
1716390978516942848 |