Caroline Chopra Department of Applied Mathematics, University of Leeds, Leeds LS2 9JT England
The pioneering experimental work of Alexandre Fabiato in 1975 saw the birth of new opportunities for thousands of scientists in the area of intracellular calcium dynamics. In examining skinned cardiac muscle cell preparations Fabiato showed that the elevated intracellular calcium concentration preceding muscle contraction was largely due to the release of calcium from intracellular stores and modulated by intracellular calcium itself, hence introducing the important role of calcium as an intracellular signaler. (For calcium reviews see[2-3,7,11,13].)
This finding opposed the then current hypothesis that intracellular calcium concentration increased mainly due to a flow of ions through electrically activated channels in the cell membrane. These revelations offered a new lease of life to mathematicians and computational biologists striving to better understand the heart beat through modelling.
The 1980s saw the development of a whole host of dynamical systems ranging in scope from simple qualitative models to very complicated biophysical models, taking the form of highly nonlinear coupled ordinary differential equations, which, more or less, accurately described Fabiato's results. Linear stability analysis and bifurcation theory have helped further understanding of these systems, together with the use of other general results relating to the existence and stability of periodic orbits (see [1,4-5] for useful texts in this area). In the wake of technological advances, such as super computers, it has become possible to perform extensive numerical experiments on these models incorporating both time and space variations of the calcium transient. In particular models have been spatially extended to form reaction-diffusion systems, the analytical and numerical work performed on generic excitable systems (extensively reviewed in ) has been adapted to these systems providing useful qualitative results [8-9]. Such simulations have provided illuminating results which in turn has set a precedent for further experimental work.
The advent of experimental techniques such as intracellular fluorescent dyes and voltage clamping in the late 1980s meant not only that the behaviour of the calcium transient, in the absence of electrical activity, could be studied whilst retaining the membrane but also that a direct measurement of intracellular calcium concentration could be obtained based on the fluorescence of the inserted dye. These advances have provided mathematicians with further evidence on which to base new and more accurate models.
The 1990s saw the discovery of intracellular calcium spiral waves in the presence of an obstacle, such as the nucleus. Mathematicians working in this area are currently trying to obtain more information on these types of waves by relating the speed normal to a curved travelling wave front to its curvature, known as the eikonal equation, in an attempt to characterise this behaviour. It is of interest to know whether calcium dynamics follow the same speed-curvature result as other dynamical systems and to postulate as to why differences may occur. Also of current interest are the conditions required for the anchoring of wave fronts to obstacle's. Although little is known about these conditions the current belief is that they are intrinsically linked to the excitability of the system and anchoring can therefore be achieved, or eliminated, by tweaking the model parameters.
From plane waves to spiral waves, the heartbeat is an interdisciplinary mystery evoking many questions with the answers lying in further collaborative research between experts from many scientific backgrounds. Further details of current research in and relating to this area are to be found at http://www.cbiol.leeds.ac.uk/.