Email updates

Keep up to date with the latest news and content from BMC Systems Biology and BioMed Central.

Open Access Research article

Computational modelling elucidates the mechanism of ciliary regulation in health and disease

Nikolay V Kotov12, Declan G Bates3, Antonina N Gizatullina2, Bulat Gilaziev2, Rustem N Khairullin4, Michael ZQ Chen56, Ignat Drozdov7, Yoshinori Umezawa8, Christian Hundhausen8, Alexey Aleksandrov9, Xing-gang Yan10, Sarah K Spurgeon10, C Mark Smales1 and Najl V Valeyev1*

Author Affiliations

1 Centre for Molecular Processing, School of Biosciences, University of Kent, Canterbury, Kent CT2 7NJ, UK

2 Biophysics & Bionics Lab, Department of Physics, Kazan State University, Kazan 420008, Russia

3 Centre for Systems, Dynamics and Control, College of Engineering, Mathematics and Physical Sciences, University of Exeter, Harrison Building, North Park Road, Exeter EX4 4QF, UK

4 Inter-regional Diagnostic Centre, Karbisheva-12A, Kazan 420101, Russia

5 Department of Mechanical Engineering, The University of Hong Kong, Pokfulam Road, Hong Kong

6 School of Automation, Nanjing University of Science and Technology, 200 Xiao Ling Wei, Nanjing 210094, P. R. China

7 Centre for Bioinformatics, Department of Computer Science, School of Physical Sciences & Engineering, King's College London, Strand, London WC2R 2LS, UK

8 St. John's Institute of Dermatology, King's College London, 9th Floor Tower Wing, Guy's Hospital, Great Maze Pond, SE1 9RT London, UK

9 Laboratoire de Biochimie, CNRS UMR7654, Department of Biology, Ecole Polytechnique, 91128 Palaiseau, France

10 School of Engineering and Digital Arts, University of Kent, Canterbury, Kent CT2 7NT, UK

For all author emails, please log on.

BMC Systems Biology 2011, 5:143  doi:10.1186/1752-0509-5-143

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1752-0509/5/143


Received:13 June 2011
Accepted:15 September 2011
Published:15 September 2011

© 2011 Kotov et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Ciliary dysfunction leads to a number of human pathologies, including primary ciliary dyskinesia, nephronophthisis, situs inversus pathology or infertility. The mechanism of cilia beating regulation is complex and despite extensive experimental characterization remains poorly understood. We develop a detailed systems model for calcium, membrane potential and cyclic nucleotide-dependent ciliary motility regulation.

Results

The model describes the intimate relationship between calcium and potassium ionic concentrations inside and outside of cilia with membrane voltage and, for the first time, describes a novel type of ciliary excitability which plays the major role in ciliary movement regulation. Our model describes a mechanism that allows ciliary excitation to be robust over a wide physiological range of extracellular ionic concentrations. The model predicts the existence of several dynamic modes of ciliary regulation, such as the generation of intraciliary Ca2+ spike with amplitude proportional to the degree of membrane depolarization, the ability to maintain stable oscillations, monostable multivibrator regimes, all of which are initiated by variability in ionic concentrations that translate into altered membrane voltage.

Conclusions

Computational investigation of the model offers several new insights into the underlying molecular mechanisms of ciliary pathologies. According to our analysis, the reported dynamic regulatory modes can be a physiological reaction to alterations in the extracellular environment. However, modification of the dynamic modes, as a result of genetic mutations or environmental conditions, can cause a life threatening pathology.

Background

Cilia are cellular protrusions which have been conserved in a wide range of organisms ranging from protozoa to the digestive, reproductive and respiratory systems of vertebrates [1]. Mobile or immotile cilia exist on every cell of the human body [2] and the insufficiently recognised importance of the cilium compartment in human physiology has been recently highlighted [1,3]. Cilia are present on most eukaryotic cell surfaces with the exception of the cells of higher plants and fungi [4]. Ciliary motility is important for moving fluids and particles over epithelial surfaces, and for the cell motility of vertebrate sperm and unicellular organisms. The cilium contains a microtubule-based axoneme that extends from the cell surface into the extracellular space. The axoneme consists of nine peripheral microtubule doublets arranged around a central core that may or may not contain two central microtubules (9+2 or 9+0 axoneme, respectively). Cilia can be broadly classified as 9+2 motile cilia or 9+0 immotile sensory cilia, although there are examples of 9+2 sensory cilia and 9+0 motile cilia. In mammals, motile 9+2 cilia normally concentrate in large numbers on the cell surface, beat in an orchestrated wavelike fashion, and are involved in fluid and cell movement. In contrast to motile cilia, primary cilia project as single immotile organelles from the cell surface. Primary cilia are found on nearly all cell types in mammals [5] and many are highly adapted to serve specialized sensory functions. The 9+2 cilia usually have dynein arms that link the microtubule doublets and are motile, while most 9+0 cilia lack dynein arms and are non-motile. In total, eight different types of cilia has been identified to date [6]. In this study, we investigate the mechanism of movement regulation for the motile type of cilia.

Although each individual cilium represents a tiny hair-like protrusion of only 0.25 μm in diameter and approximately 5-7 μm in length, cilia covering human airways can propel mucus with trapped particles of length up to 1 mm at a speed of 0.5 mm/second [7]. Such efficiency can be achieved due to the coordination between cilia and stimulus-dependent regulation of the rate of cilia beat. Dysfunction of ciliary regulation gives rise to pathologic phenotypes that range from being organ specific to broadly pleiotropic [3]. A link between ciliary function and human disease was discovered when individuals suffering from syndromes with symptoms including respiratory infections, anosmia, male infertility and situs inversus, were shown to have defects in ciliary structure and function [6].

Microscopic organisms that possess motile cilia which are used exclusively for either locomotion or to simply move liquid over their surface include Paramecia, Karyorelictea, Tetrahymena, Vorticella and others. The human mucociliary machinery operates in at least two different modes, corresponding to a low and high rate of beating. It has been shown that the high rate mode is mediated by second messengers [8], including purinergic, adrenergic and cholinergic receptors [9-19]. This mode enables a rapid response, which can last a significant period of time, to various stimuli by drastically increasing the ciliary beat frequency (CBF). At the same time, several ciliary movement modes have been reported in a ciliate Paramecium caudatum [20]. The remarkable conservation of ciliary mechanisms [21-25] creates grounds for the speculation that there can more than two ciliary beating modes in human tissues. It is, therefore, reasonable to suggest, that some human diseases, associated with aberrant ciliary motility, can arise due to modifications in the beating mode. Clearly, the development of therapeutic strategies against ciliary-associated pathologies will require advanced understanding of ciliary beating regulation mechanisms.

The periodic beating of cilia is governed by the internal apparatus of the organelle [26]. Its core part, the axoneme, contains nine microtubule pairs encircling the central pair. The transition at the junction of the cellular body and the ciliary axoneme is demarcated by Y-shaped fibres, which extend from the microtubule outer doublets to the ciliary membrane. The transition area, in combination with the internal structure of the basal body, is thought to function as a filter for the cilium, regulating the molecules that can pass into or out of the cilium. Ciliary motility is accomplished by dynein motor activity in a phosphorylation-dependent manner, which allows the microtubule doublets to slide relative to one another [12]. The dynein phosphorylation that controls ciliary activity is regulated by the interplay of calcium (Ca2+) and cyclic nucleotide pathways. The beating pattern of cilia consists of a fast effective stroke and a slower recovery stroke. During the effective stroke cilia are in an almost upright position, generating force for mucus movement. During the recovery stroke, the cilia are recovering from the power strike to the original position by moving in the vicinity of the cell surface.

Current theories which attempt to explain the workings of the Ca2+-dependent CBF regulation mechanism are incomplete and highly controversial. Elevation of intraciliary Ca2+ is one of the major regulators of ciliary movement. Calcium influx regulates ciliary activity by increasing intraciliary Ca2+ only, while the cytosolic bulk remains at a low level. Separate ciliary compartmentalisation for Ca2+ allows prolonged activation of ciliary beating without damaging the cell through high Ca2+ concentrations. It is well known that calcium fluxes via calcium channels lead to changes in organisms' swimming behaviour [27-29]. In mucus-transporting cilia, Ca2+ mediates CBF increase [19,30-32]. It has also been shown that there are some differences in the Ca2+-dependent CBF regulation in single cell organisms and in humans [12]. Sustained CBF increase requires prolonged elevation of Ca2+ levels which can be lethal to the cell [33,34]. It has been suggested that Ca2+-dependent ciliary regulation takes place locally in the vicinity or within the ciliary compartment, almost independently from intracellular Ca2+ concentration [35]. Given that the gradient of free Ca2+ in the cytosol dissipates within 1-2 seconds [36], it appears more likely that cilia form their own compartment where Ca2+ is regulated by active Ca2+ transport in a similar fashion to the intracellular Ca2+ regulatory system. This hypothesis resolves the problem of maintaining physiological levels of intracellular Ca2+ concentration. A number of experimental studies have reported several controversial results relating to the Ca2+-dependent mechanism of cilia regulation. For example, it has been reported that spontaneous cilia beat does not require alterations in Ca2+ [31,35], while nucleotide-dependent CBF increase requires Ca2+ [8]. It has also been shown that uncoupling between Ca2+ and CBF can be achieved by inhibition of Ca2+-dependent protein calmodulin (CaM) or the cyclic nucleotide pathway [19,32,37].

These findings suggest that intraciliary Ca2+ does not regulate cilia beat in isolation, but instead does so as part of more complex signalling network. Although it was originally believed that Ca2+, cyclic adenosine monophosphate (cAMP) and guanosine monophosphate (cGMP) regulate ciliary beat independently, numerous reports now strongly indicates that all three pathways are tightly interconnected [12,19,32,38-42]. The cAMP-dependent protein kinase (PKA) phosphorylates dynein in the bases of cilia and thereby increases the forward swimming speed in Paramecium [43-45]. Similar effects have been reported for PKA-dependent phosphorylation of axonemal targets in mammalian respiratory cilia [46]. Several lines of evidence indicate that PKA and cGMP-dependent kinase (PKG) both phosphorylate specific axonemal targets in a cAMP and cGMP-dependent manner. A schematic diagram for the underlying biochemical machinery for cilia movement regulation is shown in Figure 1C. It is striking that, despite significant experimental characterisation of this system, there is still rather limited mechanistic understanding of how intraciliary Ca2+ and nucleotide interplay relates to CBF.

thumbnailFigure 1. Schematic diagram of the Ca2+-dependent regulation of motile cilium. (A) Schematic representation of the motile cilium movement trajectory for one complete beating cycle. The cycle is divided into two phases: effective stroke and recovery stroke shown by red and green arrows, respectively. The frequency and the direction of cilia movement is regulated in a highly complex Ca2+ and membrane voltage-mediated manner. (B) A typical cilium consists of an axoneme of nine doublet microtubules. The axoneme is surrounded by a specialized ciliary membrane that is separated from the cell membrane by a zone of transition fibres. This separation creates an intraciliary compartment where key regulatory events take place somewhat independently from the cell body. In particular, intraciliary Ca2+ concentration can significantly differ from intracellular levels. (C) The ciliary motion is regulated by intraciliary Ca2+ levels. The Ca2+ concentration depends on the interplay of ion channels and membrane potential, Vm. The intraciliary Ca2+ concentration is dependent on currents via Ca2+, <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M101','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M101">View MathML</a>, and K+ channels, <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M102','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M102">View MathML</a>, the system of active, <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M103','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M103">View MathML</a>, and passive, <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M104','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M104">View MathML</a>, ions removal, Ca2+ leakage current, <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M105','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M105">View MathML</a>, hyperpolarisation-activated currents <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M106','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M106">View MathML</a>, inward current, I0, and the cilium-to-cell body current, <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M107','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M107">View MathML</a>. The conductivities of the channels are modulated by membrane potential. The result of the cross-talk between membrane potential and a variety of channels is that intraciliary Ca2+ can shift between several dynamic modes. The steady-state and dynamic Ca2+ alterations regulate the intraciliary levels of cAMP and cGMP in a Ca2+-CaM-dependent manner via the AC, GC and PDE isoforms [65]. Cyclic nucleotides, in turn, define the degree of phosphorylation of dynein filaments in the bases of ciliary axoneme via PKA and PKG kinases. Phosphorylation of dynein filaments regulates the relative doublet microtubules shift and thereby translates to the overall ciliary movement.

Another major regulator of ciliary beating is the membrane potential. A number of studies have reported the voltage-dependent effects of ciliary beating. The ciliate Didinium Nasutum has been shown to respond both to hyper- and de-polarization of the membrane [47]. The transmembrane potential alterations were shown to be mediated via the potential-dependent Ca2+ channels [48]. Electrophysiological studies in Paramecium caudatum have revealed complex relationships between ciliary Ca2+ currents, intraciliary Ca2+ concentration and transmembrane potential in the regulation of ciliary motility [49-55].

A number of previous computational studies have analysed various aspects of cilia movement regulation. One earlier model assessed the degree of synchronization between small ciliary areas [56]. The effects of viscosity have been investigated in mucus propelling cilia in [57]. The authors found that increasing the viscosity not only decreases CBF, but also changes the degree of correlation and synchronization between cilia. The mechanical properties of cilia motion were studied in an attempt to understand the ciliary dynamics in [58]. The authors concluded that bending and twisting properties of the cilium can determine self-organized beating patterns. While these reports offer valuable insights into the regulatory mechanisms of cilia, a number of essential questions remain unresolved. For example, there has not been a detailed analysis of how individual Ca2+ currents influence intraciliary Ca2+ levels. It also remains unclear how Ca2+ modulates nucleotide levels and membrane potential, and how such regulation affects ciliary movement. None of these reports have elucidated the underlying mechanisms governing the interplay between intraciliary Ca2+ and nucleotide alterations and CBF.

In this study, we integrate the available experimental information on the molecular pathways that regulate intraciliary Ca2+ concentration into a comprehensive mathematical model. By applying systems analysis, we elucidate the mechanisms of intraciliary Ca2+ spike generation, analyse the properties of such spikes and demonstrate the conditions under which the Ca2+ surges can become repetitive. We carry out detailed investigations of the individual current contributions to the regulation of the intraciliary Ca2+ concentrations and elucidate both steady-state and dynamic responses of Ca2+ currents and intraciliary Ca2+ concentration dynamics in response to the altered transmembrane potential shift. The model allows detailed elucidation of transmembrane potential and intraciliary Ca2+ coupling.

We employ the proposed model in order to understand the underlying molecular mechanisms of the crosstalk between Ca2+, membrane potential and nucleotide pathways that regulate ciliary movement. The systems model allows detailed analysis of the individual current contributions to the intraciliary homeostatic Ca2+ levels. Furthermore, we establish specific regulatory mechanisms for Ca2+ and cyclic nucleotide-dependent cilia movement characteristics. Crucially, our model predicts the possibility of several ciliary beating modes and describes specific conditions that initiate them. Specifically, we describe intraciliary Ca2+ dynamic modes that regulate healthy and pathologic cilia beating. We use these findings in order to propose experimentally testable hypotheses for possible therapeutic interventions in human diseases associated with pathologic cilia motility.

Results

A new model for the interplay between Ca2+ and K+ currents and transmembrane potential alterations

A new model for the regulation of ciliary movement that combines multiple Ca2+ and K+ currents [59-62] and transmembrane potential has been developed. In this model, the intraciliary Ca2+ levels are modulated by Ca2+ currents through the channels of passive and active Ca2+ transport, the current from the cilium into the cell body, the Ca2+ leakage current, and depolarisation and hyperpolarisation-activated currents. Variable extracellular conditions have continuous impact on the transmembrane potential which is intertwined with transmembrane ion currents and intraciliary Ca2+ homeostasis.

The overall network that regulates ciliary movement is divided into several functional modules (Figure 1C). One module combines all Ca2+ and K+ currents that define intraciliary Ca2+ homeostasis and the transmembrane potential. One of the most essential intraciliary Ca2+ binding proteins, CaM [63,64], selectively regulates the activities of adenylate cyclase (AC), guanylate cyclase (GC) and phosphodiesterases (PDE), and thereby modulates the intraciliary levels of adenosine monophosphate (cAMP) and guanosine monophosphate (cGMP) in a Ca2+ dependent manner [65]. The cAMP- and cGMP-dependent kinases phosphorylate dynein proteins in the bases of cilia and thereby induce the mechanical cilia movement. The complete set of equations making up the proposed model is presented in the Methods section. Below we provide a number of new insights into the mechanism of cilia regulation via a detailed investigation of the properties of this model.

The mechanism of Ca2+-dependent inhibition of Ca2+ channels

A subset of intraciliary Ca2+ channels have been reported to operate in an intraciliary Ca2+ dependent manner and have been proposed as major regulators of ciliary beat [49-51]. It is established that Ca2+ current is not inhibited by the double pulse application of depolarization impulses under voltage clamp conditions in those situations when the first transmembrane potential shift is equal to the equilibrium Ca2+ potential (+120 mV) [66]. Further experimental evidence reveals that Ca2+ current inactivation kinetics are delayed when Ca2+ ions are partially replaced by Ba2+ ions [67-71]. Altogether these findings suggest that the channels are not inhibited directly by the depolarizing shift of transmembrane potential, but that instead their conductivity is dependent on the intraciliary Ca2+ concentration. Some decrease of the inward current amplitude (by approximately 25%) upon transmembrane potential shift into the Ca2+ equilibrium level can be explained by the fact that K+ currents can contribute to the overall current measurements. Here we consider the intraciliary Ca2+ concentration-dependent Ca2+ channel inhibition and employ the developed model to analyse two potential scenarios for the Ca2+ channel conductivity regulation. In one case, Ca2+ ions bind to the Ca2+ binding site on the channel and thereby inhibit the channel's conductivity by direct interaction. The other possibility is that the Ca2+ binding protein interacts with the Ca2+ ion first and then this complex binds to the channel and inhibits its conductivity. In both cases the conductivity dependence on transmembrane potential is assumed to be monotonic according to the experimental data [66].

Direct Ca2+-dependent Ca2+ channel conductivity inhibition

We first investigate a potential intraciliary Ca2+ regulatory mechanism via direct Ca2+ ion binding-dependent Ca2+ channel inhibition. The relationship between the Ca2+ channel conductivity and the transmembrane potential has been experimentally characterized by an early study in Paramecium species [66]. The equation (18) in the Methods section approximates the experimentally established dependence. In the case of direct Ca2+-dependent Ca2+ channel inhibition, one can show that the nullclines for non-dimensional Ca2+ concentration (<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M1">View MathML</a>) and for the number of open channels (<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M2">View MathML</a>) from the system of differential equations (20) intersect at one stable point for all values in the physiological range of model parameters. The numerical solutions of the coupled differential equations (20) allow us to obtain the solutions for how the steady-state Ca2+ levels depend on the transmembrane potential. The model predictions for the steady-state ciliary Ca2+ channel conductivity dependence on the transmembrane potential under the voltage clamp conditions are shown on Figure 2A. The extracellular conditions are subject to constant change both in the case of ciliates as well as for multicellular organisms. The modifications in the external environment continuously shift the transmembrane potential. In order to estimate how the transmembrane potential alterations affect the inward Ca2+ current we derived the dependence for the Ca2+ ion flow (equation (21) in the Methods section). The experimental data-based (Figure 2A) [66] model for the inward Ca2+ current dependence on membrane potential (Figure 2B) predicts a significant reduction of the inward Ca2+ current amplitude as a function of membrane depolarization.

thumbnailFigure 2. Static and dynamic intraciliary Ca2+ concentration levels under membrane potential fixed conditions. (A) The model predictions for the relationship between the membrane potential and intraciliary Ca2+ concentration. (B) The maximum amplitude of Ca2+ current generated in response to the transmembrane potential shift is shown as function of the applied shift. The model predictions are obtained by setting the parameter α to 2, 3, and 4 in equations 18-21. The parameter α reflects the steepness of Ca2+ current dependence on membrane voltage. The model-based analysis shown here unravels the bell-shaped dependence of intraciliary Ca2+ concentration on membrane potential and the inverse relationship between the Ca2+ current magnitudes and the degree of membrane depolarisation.

We next investigated the dynamic alterations of the intraciliary Ca2+ concentrations and the inward Ca2+ current in response to the transmembrane potential shifts. Equations (20) and (21) in the Methods section were used for quantitative estimations of the Ca2+ concentration and current responses to the normalised membrane potential shifts. The model predicts that the rapid depolarising alterations of the ciliary transmembrane potential leads to the generation of single Ca2+ spikes (Figure 3). Such responses can take place when Ca2+ channels are inhibited by Ca2+ ions and the channel's conductivity depends on the membrane potential in a monotonic manner (equation 18). The amplitude of those impulses depends on the steepness of the Ca2+ channels conductivity dependence on the membrane potential. The mechanism of the Ca2+ spike generation is mainly due to the delay of the Ca2+-induced inhibition with respect to the Ca2+ conductivity alteration characteristic times. The described mechanism of the Ca2+ spike generation will only work if the Ca2+ currents significantly alter the Ca2+ concentration in the intraciliary compartments.

thumbnailFigure 3. Systems model predictions for intraciliary Ca2+ concentration and current responses to the transmembrane potential shift. The intraciliary Ca2+ concentrations (A) and Ca2+ currents (B) are calculated according to the model with direct Ca2+-mediated ciliary Ca2+ channels inhibition in response to variable degree of transmembrane potential shift. The responses are colour coded according to the degree of applied transmembrane potential shift. The violet and red coloured lines represent dynamic responses obtained after the smallest and the largest depolarising shift of membrane potential, respectively. The non dimensional membrane potential values following voltage shift from initial ψ0 = -1.2 are shown in (C) and remain the same throughout the figure. The calculations suggest that ciliary membrane depolarisation induces an intraciliary Ca2+ spike over a wide physiological range of depolarising conditions (C, E and G), whereas the current generates a biphasic response (D, F and H). The calculations were carried out for three different levels of parameter α, which reflects the steepness of the Ca2+ channels conductivity dependence on membrane potential.

Indirect Ca2+ channel conductivity regulation

In the previous section we considered Ca2+-dependent Ca2+ channel regulation under the assumptions that Ca2+ channels have an intracellular Ca2+ binding site and Ca2+ ion binding closes the channels. However, several experimental studies have suggested that the conductivity of Ca2+ channels in cilia can also be regulated indirectly, via a Ca2+ binding protein. At present, there is no direct experimental evidence that explicitly favours either direct or indirect regulatory mechanism. We, therefore, investigated the second possibility for indirect Ca2+-dependent Ca2+ channels conductivity inhibition.

The model for Ca2+ ion interactions followed by the interactions with the Ca2+ channels (described in the Methods section) is developed in line with a previously suggested modelling methodology for Ca2+-CaM interactions [63,64]. The model predictions for the number of open channels as a function of intraciliary Ca2+ concentration are shown on Figure 4 as derived by equation (25) in the Methods section. Given that the exact nature of the Ca2+ binding protein acting as a mediator between Ca2+ ions and Ca2+ channels is not established, equation (25) has been solved for different physiologically possible ratios of total number of channels to the dissociation constant for the Ca2+- binding protein interactions with Ca2+ channels. According to the derived models, the comparison of direct versus indirect Ca2+ channels inhibition can be carried out by setting the non dimensional Ca2+ concentration to zero (u = 0 in equation (25)). In this case, the predictions of equation (25) for indirect Ca2+ channels conductivity inhibition almost coincide with the model for direct Ca2+-dependent inhibition. In both models most of the Ca2+ channels are open in the lower range of Ca2+ concentrations. However, the model for indirect Ca2+ channels inhibition predicts that the number of open Ca2+ channels would equal <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M3">View MathML</a> when intraciliary Ca2+ reaches high concentrations. In other words, in the case of the indirect mechanism of inhibition, high Ca2+ does not inhibit the channels completely. The number of remaining channels in the open state would depend on the concentration of the regulatory Ca2+ binding protein.

thumbnailFigure 4. The predicted probabilities for Ca2+ channels to be in the open state. The indirect Ca2+-mediated inhibition of Ca2+ channels via a Ca2+-sensor protein, Cac0, is shown for a range of Ca2+-binding protein concentrations. The higher levels of Ca2+-sensor contributes to a greater degree of inhibition, although the general mechanism remains the same.

We next analysed the dynamics of the ciliary Ca2+ channels inhibition in response to a change in intraciliary Ca2+ concentration and in the transmembrane potential. As mentioned earlier, the Ca2+ channel conductivity is not inhibited by the membrane potential, but rather has a monotonic dependence on the transmembrane potential difference as shown on Figure 2B. We found that the characteristic time <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M4">View MathML</a> of Ca2+ channel alterations in response to step changes in Ca2+ is inversely proportional to the total number of channels. The model predictions for intraciliary Ca2+ concentration and inward Ca2+ current in response to depolarising transmembrane potential changes from V0 to V1, are shown on Figure 5. According to the derived equation (30) in the Methods section, the inward Ca2+ current changes with a characteristic time τV, which can only be estimated by considering the K+ current contribution. In the first instance we only consider active and passive Ca2+ transport (equation (31) in Methods). We neglected by the kinetics for the channels of active transport due to the assumption that the kinetics of alterations of active transport are much faster that the characteristic alteration times of passive transport (equation (32) for the Ca2+ currents). The model predictions for the steady-state dependence of transmembrane potential on intraciliary Ca2+ concentration and the dynamic Ca2+ current amplitude dependence on the membrane potential are shown in Figure 6A and 6B, respectively.

thumbnailFigure 5. The responses of Ca2+ current and concentration to the potential shift via indirect Ca2+ inhibition. The intraciliary Ca2+ concentration (A) and Ca2+ current (B) responses to membrane depolarisation are computed under the assumption of indirect Ca2+ channel inhibition for different concentrations of Ca2+-sensor protein, Cac0. The model predicts that in response to fast membrane depolarisation, the ciliary coupled Ca2+ and membrane potential system generates an intraciliary Ca2+ concentration spike similar to the one observed under the direct Ca2+-inhibition assumptions. The non dimensional membrane potential values following voltage shift from the initial ψ0 = -1.2 are shown in (C) and are the same throughout the figure. The major difference between direct (Figure 3) and indirect mechanisms of inhibition occurs in the steady-state levels of Ca2+ concentration that takes place after the transitional process. The comparison of intraciliary Ca2+ concentration responses for increasing levels of the Ca2+ sensor protein (C, E and G) predicts that intraciliary Ca2+ concentrations following the voltage shift are inversely dependent on the Ca2+ sensor concentration, Cac0. Different concentrations of Ca2+ sensor protein do not affect the membrane depolarisation-induced Ca2+ current (D, F and H).

thumbnailFigure 6. The membrane potential effects on Ca2+ current under fixed voltage. (A) The relationship between membrane potential and intraciliary Ca2+ concentration is shown under voltage clamp conditions. (B) The inward current amplitude dependence as a function of a fast shift in the holding membrane potential. Calculations were performed for different concentrations of Ca2+ sensor protein, CaC0. The model predicts qualitatively similar responses for a relatively wide range of Cac0 concentrations.

The cilia's external environment is subject to constant change and can significantly affect the behavioural responses of ciliates and modulate ciliary beating in multicellular organisms. In order to account for the effects of extraciliary Ca2+ variations we estimated the amplitudes of intraciliary Ca2+ spike generation under different external Ca2+ concentrations and variable transmembrane potentials. Figure 7A shows that the increase of the Ca2+ concentration in the external solution increases the amplitude of the generated intraciliary Ca2+ spike. The amplitude of the spikes goes to zero when the membrane potential equals the equilibrium potential for Ca2+ ions. The increase of membrane potential decreases the amplitude of the Ca2+ current (Figure 7B). This effect is due to the increase of the steady-state intracilia Ca2+ level and Ca2+-dependent inhibition of the Ca2+ channels.

thumbnailFigure 7. The model predictions for Ca2+ current in response to membrane potential shift. The model predictions for Ca2+ current amplitude shown as a function of extracellular Ca2+ concentration and the initial membrane potential values before the voltage shift was applied. (A) The dependence of Ca2+ current amplitude on the transmembrane potential alterations is calculated for different extracellular Ca2+ concentrations. (B) The graph shows the Ca2+ current amplitude dependence on the initial values of the potential when the depolarising voltage shift was applied. The model shows that the ciliary system properties are reasonably conserved under significant variations of extracellular Ca2+ and over a wide range of membrane potential values.

We noted earlier that there is a Ca2+ current in the cilia which transfers ions from the cilia into the cellular compartments. This current can be described by equation (11) in Methods. The contribution of cilium-to-cell body current to the intraciliary Ca2+ concentration dynamics was evaluated experimentally in [72,73]. It was shown that under depolarized membrane potential conditions the contribution of this current is very small and the intraciliary Ca2+ is mainly pumped out of the cilia into the extracellular space by the active Ca2+ transport. According to other observations, Ca2+ current from cilia into the cellular compartment can be larger than the current generated by the active Ca2+ transport. In order to investigate the role and contribution of the cilia-to-cell compartment current, we introduced its contribution to the intraciliary Ca2+ concentration dynamics (equation (34)). We performed qualitative analysis of the Ca2+ concentration alterations in the cilia in the presence of the cilium-to-cell current and compared the Ca2+ dynamics with the case when this current was not present. We found that although the cilium-to-cell body current influences the intraciliary Ca2+ concentration levels, it does not change the dynamics qualitatively when the membrane potential is depolarized and fixed.

Our findings suggest that the cilium represents an excitable system with unique properties. The Ca2+-dependent inhibition of Ca2+ channels inhibition allows for the generation of single impulses of variable amplitude proportional to the degree of membrane depolarisation caused by variations in the external concentrations of ions. This system is able to generate a single spike despite unpredictable variations of ionic concentrations in the environment and is, therefore, very robust to alterations in the external conditions. Another interesting aspect of the ciliary excitation is the ability of the system to generate regulatory intraciliary Ca2+ impulses proportional to the degree of membrane depolarisation (Figures 3 and 5). This property can allow cells to sense and "automatically" respond to alterations in their environment.

The contribution of K+ currents

In the previous section, we analysed the dynamic properties of the intraciliary Ca2+ system under voltage clamp conditions. Several lines of evidence suggest that K+ currents contribute to the currents registered in cilia under voltage clamped conditions. The existence of K+ currents in cilia is supported by a number of experimental studies. The experimental data shows that the measured current is not equal to zero when the membrane potential equals the equilibrium membrane potential for Ca2+ ions. Instead, the current equals zero when membrane potential is about 10 mV while the equilibrium potential for Ca2+ ions equals 120 mV [66]. This observation suggests that both Ca2+ and K+ currents contribute to the overall current measured at early stages of current registration under voltage clamp, and therefore both currents need to be taken into the consideration in order to advance understanding of the mechanisms involved in ciliary regulation. At the same time, it has so far been impossible to register Ca2+ currents by inhibiting the K+ contribution. Various compounds can only partially block the K+ current when applied from inside of the membrane. Ciliary K+ currents have also been measured separately from Ca2+ currents.

In order to account for the contribution of K+ currents to the regulation of intraciliary Ca2+ concentration, we developed a model for the regulation of K+ currents by membrane potential (equation (35) in Methods). The dependence of K+ conductivity on the membrane potential is described by equation (36). Figure 8 shows the measured experimental values and the approximating curve calculated according to equation (36). The current on Figure 8 is normalized to 1 when membrane potential equals 0. In order to approximate the current, we only used the experimental values obtained under the membrane depolarized conditions, and this experimental data is approximated by equation (36). The predictions for intraciliary Ca2+ dynamics and Ca2+ impulse amplitude in response to the shift in membrane potential are shown in Figure 9A and 9B, respectively. A comparison of the membrane potential shift-induced intraciliary Ca2+ spike generation in the presence of the K+ currents (Figure 9A) with the model when K+ channels are not incorporated (Figure 3 and 5), suggests that while the K+ currents change the quantitative values of the Ca2+ concentration dynamics, the general shape of the Ca2+ response remains the same.

thumbnailFigure 8. Current-voltage characteristic for the steady-state K+ current. The model predictions for the K+ current as a function of membrane potential are compared with the experimentally measured voltage-current. The model elicits accurate agreement with the experimental data for the physiologically relevant range of membrane potential alterations.

thumbnailFigure 9. Ca2+ concentration dynamics response to transmembrane potential shift in the presence of K+ current. (A) The intraciliary Ca2+ dynamics is computed in response to transmembrane potential shift by incorporating the contribution from K+ channels to the overall current. Comparison of these results with Figures 3 and 5 suggests that while K+ currents introduce some quantitative changes, the overall qualitative characteristics of the response remain the same. The non dimensional membrane potential values following voltage shift are indicated for each predicted response. (B) The amplitude of the generated Ca2+ concentration response is shown as a function of membrane potential alteration. The amplitude of the current is approximately linearly inversely proportional to the magnitude of the applied membrane potential shift.

The transmembrane potential dynamics in the absence of voltage clamp

In the previous sections we investigated the mechanisms of the transmembrane potential shift-dependent Ca2+ spike generation under voltage clamp conditions. However, Ca2+ currents themselves can alter the membrane potential. Here we incorporate the membrane potential dependence on Ca2+ currents and investigate the membrane potential dynamics in the absence of voltage clamp (equations (40) and (41)). The non dimensional Ca2+ concentration and membrane potential are described by equation (42).

The predictions for the individual Ca2+ and K+ current responses to various membrane potential shifts are shown in Figure 10. One can note that the dynamics of the responses significantly differs between the ones obtained under voltage clamp conditions and the situation when the membrane potential is not fixed, but dependent on the currents (Figures 3, 5 and 9). Voltage current characteristics for Ca2+, K+ and full currents calculated using the current amplitudes from the currents dynamics are shown on Figure 10H. One can note that the current is an almost exponentially growing function of membrane potential.

thumbnailFigure 10. The dynamics of current alterations in the absence of voltage clamp. The Ca2+ and K+ current dynamics were calculated in response to the membrane potential shift in the absence of voltage holding conditions (A-G). The membrane potential was depolarised from the normalised value V0 = -30 mV to -12.5 mV (A) and 25 mV (G), which represent the smallest and the largest change shift, respectively. The currents elicited by the membrane potential alteration within the range are shown in (B)-(F). <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M108','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M108">View MathML</a>, <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M109','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M109">View MathML</a> and IFull represent calcium, potassium and full currents, respectively. The comparison with the Ca2+ current responses to membrane potential shift obtained under the voltage clamped conditions (Figure 3, 5 and 9) reveals significant differences in the dynamics. (H) Voltage current characteristic is calculated using the current amplitudes from the currents dynamics for the full (IFull) and Ca2+ (<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M110','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M110">View MathML</a>) currents. The steady-state voltage current relationship is calculated according to the stationary current values, (IPS).

The monotonic dependence of Ca2+ current on transmembrane potential, and simultaneous Ca2+-dependent inhibition of Ca2+ channels, represents a classical problem of two interconnected variables: intraciliary Ca2+ and membrane potential. In this system, increasing Ca2+ current with transmembrane potential depolarisation represents a positive feedback loop mechanism, whereas the intraciliary Ca2+ concentration-dependent Ca2+ channels inhibition represents a negative feedback loop. We, therefore, sought to investigate the range of potential dynamical properties of the ciliary system emerging from the coupling of Ca2+ current and membrane potential described by equations (42).

Figure 11 shows the Ca2+ current and membrane potential dynamics and the phase diagrams for an increasing range of inward current. We found that the inward current into the cilium can modify the dynamic properties of the Ca2+-membrane potential system. In all cases, the null cline <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M5">View MathML</a> represents the N-shaped curve. In the physiological range of non-dimensional membrane potential V (50 mV, -0.025 mV) and intraciliary Ca2+ (u) from 0.04 to 40 μM, the null cline <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M1">View MathML</a> shows a monotonic growth.

thumbnailFigure 11. The model predictions for coupled intraciliary Ca2+ concentration and membrane potential alterations. The systems model for intraciliary Ca2+ concentrations, Ca2+ and K+ currents coupled with membrane potential responds to applied inward current in a dynamically diverse fashion. The dynamic mode of the ciliary system is defined by the intersection of nullclines shown on phase diagrams for three representative values of inward current, I0. A. A small inward current, I0, initiates the generation of a small pulse of intraciliary Ca2+ (B) and membrane potential (C). D. A larger inward current can shift the system into the mode of sustained oscillations of both intraciliary Ca2+ (E) and membrane potential (F). G. Significantly higher values of inward current shift the intersection of the null clines to a higher region with a stable solution. This leads to a switch type of response to a higher steady-state level of Ca2+ (H) and membrane potential (I). The model describes a novel mechanism of ciliary excitability and predicts that the ciliary system can generate either a single impulse, generate sustainable oscillations, or operate as a switch between lower and higher Ca2+ and membrane potential levels.

One can clearly see that there is significantly different response for different values of the inward current. When the influx of the ions is relatively small, the <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M1">View MathML</a> null cline intersects the <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M6">View MathML</a> = 0 null cline in the left descending area (Figure 11A); such a null cline crossing results in a stable solution. In this case the system responds by the generation of a single impulse of both intraciliary Ca2+ concentration and the membrane potential followed by a return to homeostatic levels (Figure 11A, B and 11C). Further increasing the current causes the null cline <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M1">View MathML</a> to intersect with the null cline <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M5">View MathML</a> in the middle region of the ascending area, leading to an unstable solution with a limit cycle formed around the area that represents the oscillations. (Figure 11D, E and 11F). However, further increase of the current causes the <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M1">View MathML</a> null cline to intersect with the null cline <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M5">View MathML</a> in the right descending area, resulting in a stable solution with a slight increase of the homeostatic Ca2+ and membrane potential levels (Figure 11G, H and 11I). The key conclusion from this analysis is that the external ionic conditions can initiate essentially different dynamic properties of the system regulating ciliary movement. One of the key factors that affect the ciliary beat cycle is the level of intraciliary Ca2+. Our findings suggest that in response to the external conditions, there are several possibilities for intraciliary Ca2+ upregulation. The system can generate a single spike (Figure 11B) of variable amplitude (data not shown), permanently increase Ca2+ in a dynamic fashion and maintain the high intraciliary levels (Figure 11E), or operate in a monostable multivibrator mode (cilia can generate a Ca2+ spike in response to any alteration of membrane potential) (Figure 11H). These three possibilities can be associated with the different modes of ciliary beat observed in human cilia as well as in various ciliates.

The dynamic properties of excitable systems with two interdependent variables are reasonably well understood at a theoretical level. In the present case, Ca2+ and membrane potential represent the slow and fast variables, respectively. This study, therefore, establishes that the dynamic properties of ciliary systems, where the Ca2+ and K+ channel conductivities represent monotonic function of membrane potential and the Ca2+ channels conductivity inversely depends on intraciliary Ca2+ concentration, are comparable with the properties of excitable systems based on the "N-shape" dependence of the Na2+ channel conductivity on membrane potential [74]. At the same time, it is essential to note that the mechanism of excitation described in motile cilia is different from the "classical" one described in most excitable cells and systems that involve IP3 Ca2+ channels [75,76].

The membrane hyperpolarisation-dependent currents modulate the excitatory properties of the ciliary system

The ciliary transmembrane potential can shift in two directions. In the previous section we investigated the intraciliary Ca2+ responses caused by membrane depolarisation. Here we assess the implications of the membrane hyperpolarisation which has been shown to activate the current from cilia into the cell body [77,78]. We introduced the corresponding term into our model for the Ca2+ ions movement via the membrane as a function of the corresponding membrane potential shift (equation 43). By assuming the potential independent mechanism for Ca2+ and K+ ion expulsion, the system of intraciliary Ca2+ and membrane potential is derived as shown in equation (46) in the Methods section.

The model predictions for intraciliary Ca2+ and membrane potential dynamics in the absence of the hyperpolarisation-induced current are shown in Figure 11D, E and 11F. In this case the system remains in the mode of steady oscillations. We found that such oscillatory mode can be significantly modulated by the hyperpolarisation-induced and the external conditions-dependent current. Figure 12A shows the phase diagram in the presence of the hyperpolarising current. In contrast to the situation described on Figure 11 when the hyperpolarising current was absent, the ciliary system generates a single spike in response to alteration of external ionic concentrations over the whole physiological range of the inward current (Figure 12B and 12C). Under certain combinations of the hyperpolarisation-induced currents and the Ca2+-dependent K+ currents the phase diagram modifies in a manner so that the system responds by generating a single spike in response to alteration of the inward current of any magnitude (Figure 13). This indicates that the ciliary system under membrane hyperpolarising conditions can become a monostable multivibrator.

thumbnailFigure 12. The "monostable multivibrator" mode occurs in response to membrane hyperpolarisation. The model predictions are shown for the ciliary response to membrane hyperpolarisation. The phase diagram analysis (A) suggests that the nullclines for Ca2+ and membrane potential V always intersect at a single point that correspond to a steady-state solution in the whole physiological range of possible inward current values, I0. The model shows that this type of system always generates a Ca2+ impulse in response to any alteration of inward current, I0, within the physiological range (B).

thumbnailFigure 13. The hyperpolarisation-mediated Ca2+ spike is proportional to inward current. The effects of hyperpolarisation-mediated currents are described as a function of inward current magnitude. The computational calculations predict that in the presence of hyperpolarisation induced currents, the ciliary system always generates Ca2+ impulses proportional to the applied inward current in both positive (A) and negative (B) ranges of values. The comparison of phase diagrams in the presence (C) and absence (D) of the hyperpolarisation-induced current vCah suggests that the Ca2+ spike generation in response to inward current alteration is due to the hyperpolarisation-mediated currents. In the absence of hyperpolarisation the solution becomes unstable with a limit cycle and the system undergoes sustained oscillations (D).

The role of cilia-to body Ca2+ current under membrane hyperpolarisation

Intraciliary Ca2+ concentration has been experimentally estimated to be approximately one order of magnitude higher in comparison with the intracellular levels. The Ca2+ current generated by the ion flow from the ciliary compartment into the cell has been reported by a number of groups [28,52,53], but the role this current plays in the regulation of the ciliary beat remains unclear. In order to address this question, we introduced the term for this current into the equation that describes the intraciliary Ca2+ concentration (equation (47) in the Methods section). In the absence of direct measurements of the dependence of conductivity on membrane potential we set the conductivity to increase in response to hyperpolarisation of membrane potential. A summary of the model responses for all channel conductivities as a function of membrane potential is shown in Figure 14. By numerically solving the coupled equations for intraciliary Ca2+ and membrane potential with the cilia to the cell body contribution, we found that this current does not qualitatively change the general dynamic properties of the system.

thumbnailFigure 14. Key ciliary ion channels conductivity dependence on membrane potential. The overall ciliary excitability properties are due to the unique combination of ion channel conductivity dependence on membrane potential. The normalised ion channel conductivities are shown for Ca2+ currents v(V), K+ currents <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M111','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M111">View MathML</a>, hyperpolarisation-induced currents vh(V), and the cilia-to-cell currents vt(V). The hyperpolarisation and cilia-to-cell currents are activated by membrane hyperpolarisation, whereas the activation of Ca2+ and K+ channels takes place under depolarised conditions.

Despite the lack of a noticeable contribution to the ciliary dynamic properties, this current requires a special consideration. Experimental studies have clearly demonstrated that intraciliary Ca2+ is significantly higher than intracellular Ca2+ concentration. At the same time, if the conductivity of protein structures governing the Ca2+ ions movement from cilia to the body is high, most of the intraciliary ions would move from cilia into the cell body in a very short time. A simple calculation suggests that if Ca2+ could freely flow from cilia into the body, the intraciliary concentration would become equal to the intracellular Ca2+ concentration in less than 100 μs due to the difference in the volumes of the cell body and intraciliary compartments. Experimental measurements in ciliates show that the hyperpolarisation-induced backwards movements can last longer than 100 μseconds. It is also known that the avoidance reaction that requires long term elevation of intraciliary Ca2+ concentration can be observed in hyperpolarizing solutions. During all this time the intraciliary Ca2+ concentration can be several orders of magnitude higher than the intraciliary concentration. In this study, we have demonstrated that the steady-state Ca2+ current under the depolarized membrane potential conditions can only be reduced by the Ca2+-dependent inhibition of Ca2+ channels. All these observations suggest that the Ca2+ removal from cilia to the cell body occurs in a membrane potential dependent manner.

The mechanism of Ca2+ and cyclic nucleotide-dependent CBF regulation

In addition to intraciliary Ca2+ and K+ potassium levels being coupled with the membrane potential modulation, cyclic nucleotides contribute to the regulation of one of the major ciliary beat parameters, frequency. Intraciliary Ca2+ levels activate a variety of adenylate cyclases (AC) and phosphodiesterases (PDE) that produce and hydrolyse cyclic nucleotides, respectively, and thereby modulate the intraciliary cAMP and cGMP levels. At the same time, cAMP and cGMP-dependent kinases phosphorylate dynein arms [45] in the bases of cilia and thereby induce the ciliary movement [79].

In a previous work we showed that the ciliary beat frequency can have a "double" bell shape dependence on Ca2+ concentration [65] due to the differential regulation of adenylate and guanylate cyclase isoforms in a Ca2+-calmodulin (CaM) dependent manner [63,64]. One bell-shape was due to the cAMP production by a combination of AC and PDE, whereas another one was mediated by cGMP production and degradation. Our results proposed an explanation for seemingly conflicting experimental evidence suggesting that CBF can both decrease and increase with increasing Ca2+ concentration. Here we extend our previous analysis and describe the conditions when one or the other peak can be significantly reduced or even disappear. Figure 15A shows the model predictions for the ciliary beat frequency in comparison with the experimental data [65,80]. Recent studies demonstrated that hormones and pharmacological agents can regulate both function and structure of cilia by interfering with the cyclic nucleotide signalling pathways [81,82]. These findings support the possibility for the development of novel therapeutic strategies for ciliary pathologies by modulating the dynamic mode of beating via ciliary membrane receptors [83].

thumbnailFigure 15. Predictions for ciliary beat frequency modulation. (A) The comparison of experimental data and model predictions for the CBF dependence on Ca2+ concentrations. The experimental data for CBF [80] is shown as circles whereas the model prediction is represented as continuous blue line [65]. (B) GC activity depends on a number of intracellular mediators. The model describes the correlation of the right "peak" with the GC activity levels. (C) The model predicts that extracellular signals via G-protein pathway can significantly modulate the left "peak" of the plot of CBF dependence on Ca2+ concentration. In both cases (B and C) the diminishing amplitude is shown as green and red lines. The model describes the distinct contribution of both internal and external signals to CBF modulation.

Figure 15B and 15C demonstrate that the "amplitude" of each peak can be significantly diminished if the activity of the AC or GC, respectively is modulated by a temporary or permanent, internal or external signal. Under such a scenario, CBF can only increase or decrease if it happens to be on one slope of the bell-shaped dependence. Therefore, according to our analysis, different organisms with the same underlying ciliary regulatory system can achieve all possible CBF regulatory modes as a function of Ca2+ concentration: the reverse bell-shaped dependence, if the "peak" values shown on Figure 15A occur at the lower and higher limits of the physiological range for Ca2+ concentration, the bell shape dependence that can be either cAMP and cGMP dependent, and either monotonic increase or decrease if the physiological range of Ca2+ concentrations occur at one of the slopes. Our model, therefore, describes the core Ca2+-dependent regulatory mechanisms of cilia beat, but also provides an explanation for the differences observed between cilia in different single cell organisms as well as tissue specific differences. It also unravels the mechanism for how various stimuli modulate the rate of CBF by signalling via Ca2+- and G-protein mediated pathways.

Discussion

We develop a new computational model for Ca2+ and membrane potential-dependent ciliary regulation that explains how different ciliary beating regimes are regulated. The model describes a novel mechanism of excitability based on the membrane potential-dependence of Ca2+ currents (Figure 2) and simultaneous intraciliary Ca2+-concentration mediated inhibition of Ca2+ channels (Figure 4). Our analysis shows that motile cilia constitute an excitable system with a novel mechanism of excitability. The ciliary system is able to generate a Ca2+ spike in response to a wide range of transmembrane depolarisation (Figure 3, 5 and 9). The major difference in the ciliary excitation described here, with respect to classical excitation mechanisms, is that ciliary excitability is robust to a wide range of ionic variations in the environment.

The excitability mechanism of cells in evolutionary advanced organisms is based on a combination of the N-shaped dependence of the quick inward cationic current on the transmembrane potential and slow alterations of the K+ conductivity [84-87]. The ciliary voltage-current characteristic (Figure 10H) suggests several functional dynamic modes of operation: i) single impulse generation, ii) oscillator, iii) trigger (Figure 11), all initiated by membrane depolarisation. At the same time, the hyperpolarisation-induced Ca2+ currents switch the system into the mode of a monostable multivibrator, when cilia can generate a Ca2+ spike in response to any alteration of membrane potential. The dynamics of such a system depends on the transmembrane potential. In other words, any alterations in the transmembrane potential (for example, initiated by variations of the external ion concentrations) switch functional performance of the system or make it non-excitable.

It was originally believed that Ca2+, cAMP and cGMP each represent an independent pathway of ciliary regulation, however, there is by now a significant amount of evidence that strongly suggests that all three pathways are intimately interconnected [88]. It is well established that cAMP and cGMP are synthesized by AC isoforms and hydrolysed by PDEs in a Ca2+-CaM-dependent manner. In this work we describe the mechanism of the cross talk between the three circuits and explain how CBF can be modulated via extra- and intraciliary pathways (Figure 15).

Conclusions

Therapeutic applications of systems model for intraciliary Ca2+ regulation

Our detailed analysis of the effects of several Ca2+ and potassium currents and membrane potential on intraciliary Ca2+ levels offers a new way of interpreting ciliary motility associated pathologies. There are two main Ca2+-mediated parameters that govern the motile function of cilia: the direction and the frequency of beat. Our model shows that there can be several dynamic regimes of intraciliary Ca2+ alterations during which the intraciliary Ca2+ concentration can be either at low or high levels, temporarily or for a significant period of time (Figure 16A). Experimental evidence suggests that high Ca2+ reverses the direction of cilia strike (Figure 16B), and modifies the frequency in a highly nonlinear manner (Figure 15) via synthesis and hydrolysis of cyclic nucleotides (Figure 15). In a previous work [89], we showed that genetic mutations that alter the dynamic properties of a system in a permanent manner can lead to disease. In this study, we propose a conceptually similar mechanism for the pathologies associated with ciliary motility that can take place in some human diseases [6]. Our study suggests that while CBF regulation in a multicellular organism can be modulated by a number of intracellular and extracellular factors (Figure 15), genetic mutations that directly or indirectly, affect the Ca2+-mediated CBF dependence can dramatically impair the essential processes such as clearing function in airways, male fertility, or the determination of the left-right axis during development (Figure 16C and 16D). The treatment of the pathologies associated with this mechanism would rely on the restoration of the original Ca2+-dependent CBF dependence.

thumbnailFigure 16. Systems model for ciliary excitation offers new avenues for mechanistic interpretation of ciliary pathologies. The systems model for intraciliary Ca2+ regulation offers new strategies for interpretation of experimental data and development of pharmaceutical interventions for ciliary motility-associated pathologies. The ciliary system is predicted to maintain either low (A) or high (B) levels of intraciliary Ca2+. The extracellular conditions can shift the functional modes of ciliary activity and cause a temporal, repetitive or a long term Ca2+ increase which causes cilia to reverse the direction of beat. The long term reversal of the direction of beat can explain the mechanism of the situs inversus disease, which is a congenital condition in which the major organs are mirrored with respect to their normal positions. Intraciliary Ca2+ levels modulate ciliary beat frequency [65] via either an external signal through the G-protein mediated pathways or by parametric regulation of GC activity. Such alterations can represent a physiological response to external and intracellular signals, but can also occur as a result of genetic mutations. According to our model, the latter case represents a potential pathology and in either case of permanent Ca2+-dependent CBF alteration (C or D) requires the development of therapeutic strategies to rescue the mutation-mediated alteration of the system.

Future perspective

At present there is limited understanding of the underlying biological mechanisms that govern ciliary motility. This study describes the modes of intraciliary Ca2+ dynamics in a highly detailed fashion. It shows the conditions that switch the system between the modes of Ca2+ spike generation, oscillatory dynamics and a trigger. The interdependent influences of Ca2+ and K+ currents, transmembrane potential and cyclic nucleotides modulate the ciliary beat frequency and the direction of beat in a highly nonlinear manner. The further development of mathematical models of this system is still required to represent ciliary movements as a function of Ca2+ concentration and obtain the detailed understanding of ciliary motility which will be crucial for the development of new treatments for human diseases. While the core protein regulatory machinery involved in ciliary motility is very likely to be conserved, some variations in response to increased Ca2+ between single cell ciliates and mammalian cilia have been reported [90]. We would argue that those differences are not due to the change in the mechanisms of Ca2+-dependent regulation but are rather caused by variations in the parameters of the regulatory circuits. The further investigation of single cell ciliates may allow a greater degree of characterisation of ciliary movement mechanisms, because in these systems alterations of ciliary motility translate into movement trajectories which can be easily observed.

Methods

Model Description

Figure 1 provides a schematic outline of the network regulating intraciliary Ca2+ concentration that is considered in our model. Intraciliary Ca2+ concentration is regulated by the currents of passive and active Ca2+ transport, as well as by Ca2+ leak into the extracellular space and into the cell body.

A basic mathematical model for intraciliary Ca2+ concentration and its relationship to transmembrane potential was proposed for the first time in [91]. A large number of recent experimental findings now allow the formulation of a more advanced model that includes the crucial aspects of the molecular mechanisms governing cilia movement. Below we describe the complete model for intraciliary Ca2+ regulation developed in this study.

The dynamics of intraciliary Ca2+ alteration are given by:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M7">View MathML</a>

(1)

where VR - is the cilium volume, SR - is the cilium surface area, and <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M8">View MathML</a> and <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M9">View MathML</a>-are the Ca2+ currents through the channels of passive and active Ca2+ transport, respectively. <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M10">View MathML</a> is the current from the cilium into the cell body. <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M11">View MathML</a>-is the Ca2+ leakage current. <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M12">View MathML</a> is the function that encounters Ca2+ binding to and release from CaM, the main Ca2+ binding protein in cilia, z = 2 is the Ca2+ ions charge, and F is the Faraday constant.

The dynamics of Ca2+ concentration alterations within the cilium are defined by the individual contribution of the Ca2+ currents. The current via the channels of passive Ca2+ transport is given by:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M13">View MathML</a>

(2)

where <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M14">View MathML</a> is the density of functionally active Ca2+ channels. <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M15">View MathML</a> is the time averaged current via a single Ca2+ channel. Ca2+ channel activity is modulated externally by a number of metabolic pathways. Therefore, we only consider the pool of functionally active Ca2+ channels. The time-averaged current via a single Ca2+ channel is given by:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M16">View MathML</a>

(3)

where <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M17">View MathML</a> is the conductivity of a single channel in the state i (in the most general state Ca2+ channels can have a number of states with different degrees of conductivity), <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M18">View MathML</a> is the Ca2+ potential in the equilibrium, Vm is the transmembrane potential of the cilia membrane.

The Ca2+ leakage current is given by:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M19">View MathML</a>

(4)

Assuming that the active Ca2+ transport system extrudes one Ca2+ ion per cycle, the current generated by the plasma membrane Ca2+ pump is given by:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M20">View MathML</a>

(5)

where <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M21">View MathML</a> is the density of the plasma membrane Ca2+ pump protein complexes bound to one Ca2+ ion, and iA is the time averaged Ca2+ current via a single Ca2+ channel. By introducing further assumptions that all the channels of active Ca2+ transport are saturated by ATP and that all bound Ca2+ molecules are released into the extracellular space, this current would be defined by the dynamics of the active transport channels bound to a Ca2+ ion:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M22">View MathML</a>

(6)

where <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M23">View MathML</a> is the density of the active Ca2+ transport channels, <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M24">View MathML</a> and <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M25">View MathML</a> are the association and dissociation constants for the Ca2+ ion interaction with the active Ca2+ transport channels, respectively. <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M24">View MathML</a> is the constant that defines the Ca2+ ion transition from the bound state into the Ca2+ channel. By introducing new non-dimensional variables: <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M26">View MathML</a>, (where nm is the dissociation constant of the protein regulating the passive Ca2+ transport channels), equation (6) takes the following form:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M27">View MathML</a>

(7)

The steady state-current through these active transport channels is given by:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M28">View MathML</a>

(8)

where <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M29','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M29">View MathML</a>.

The model in [91] represented the leakage current from cilium into the cell body by the following expression:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M30','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M30">View MathML</a>

(9)

where ϑ is the effective diffusion constant, and <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M31','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M31">View MathML</a>and <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M32','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M32">View MathML</a> are the intraciliary and intracellular Ca2+ concentrations, respectively. The Nernst equations allow a more accurate modelling of the leakage current from cilium into the cell body as follows:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M33','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M33">View MathML</a>

(10)

where gt(Vm) is the overall conductivity of the cilium base area, Vrt is the difference of the potential between cell body and cilia, <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M34','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M34">View MathML</a>is the intraciliary Ca2+ concentration, and <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M35','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M35">View MathML</a> is the intracellular Ca2+ concentration. Equation (10) can then be represented in the following non dimensional form:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M36">View MathML</a>

(11)

where <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M37','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M37">View MathML</a>.

In the following sections we derive the models and analyse the individual contributions of the different types of Ca2+ currents to the intraciliary Ca2+ homeostasis.

Model for intraciliary Ca2+-dependent Ca2+ channel conductivity inhibition

In this model we assume that Ca2+ channels located within cilia have a Ca2+ binding site on the intracellular site of the channel. According to such a model, Ca2+ ion binding to that site mediates the channel's transition into the closed state with no conductivity. We further assume that the characteristic time for the transition from the conductive to the non conductive states is much smaller that the characteristic Ca2+ alteration times. In that case, the dynamics of the Ca2+ channels transition into the closed state in response to Ca2+ increase is given by:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M38','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M38">View MathML</a>

(12)

where N is the number of channels in the open state, N0 is the total number of channels, np and nm are the association and dissociation constants for the interaction of Ca2+ ions with the Ca2+ channels, respectively. The steady-state solution of equation (12) is given by:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M39','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M39">View MathML</a>

(13)

where <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M40','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M40">View MathML</a>. In the non dimensional form this solution is given by:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M41">View MathML</a>

(14)

where <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M42">View MathML</a>.

When u = 0, n = 1 and for u = ∞ n = 0, the current via these channels equals:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M43','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M43">View MathML</a>

(15)

The alteration of the overall conductivity of the Ca2+ channels in the Paramecium cilia in response to the shift of transmembrane potential from V0 to V1 under the voltage clamped condition is given by:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M44">View MathML</a>

(16)

where τp is the characteristic time of the transmembrane potential alteration from V0 to V1.

The equation for the conductivity alterations (16) can be represented as follows:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M45','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M45">View MathML</a>

(17)

where <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M46','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M46">View MathML</a>, v(ψ, t) is the Ca2+ channel conductivity dependence on the transmembrane potential and on time.

An early study investigated the Ca2+ channels' conductivity dependence on transmembrane potential ψ in Paramecia [66] and approximated it by the following equation:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M47','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M47">View MathML</a>

(18)

where α, d and λ are the parameter values that allow the best representation of the available experimental data. In this model, the steepness of the dependence of the conductivity on membrane potential is represented by the parameter α.

For simplicity we assume that the elementary intraciliary volume with inward Ca2+ current does not contain any Ca2+ binding proteins and Ca2+ is extruded into the extracellular space by the active Ca2+ transport only. The Ca2+ binding to the Ca2+ channels is assumed to occur much faster than the characteristic times of Ca2+ channel inhibition. Under such assumptions, the intraciliary Ca2+ dynamics in response to the transmembrane potential shift under the voltage clamp is given by:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M48','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M48">View MathML</a>

(19)

In non dimensional form, equation (19) takes the following form:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M49','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M49">View MathML</a>

(20)

where <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M50','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M50">View MathML</a>

The initial conditions for the equations (20) when η = 0 and ψ = ψ0 are equal to u0 (ψ0) and n(u0), which can be obtained by numerical solution of these equations. The non dimensional Ca2+ current is equal to:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M51','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M51">View MathML</a>

(21)

Indirect Ca2+ channel conductivity regulation

The second model of intraciliary calcium regulation assumes the Ca2+ binding protein interacts with the Ca2+ channel and thereby causes the Ca2+ channels to close. By considering that the Ca2+ ion-Ca2+ binding protein occurs on a faster time scale than the Ca2+ binding protein-Ca2+ channels interaction, the steady-state solution for the Ca2+ binding protein in complex with a Ca2+ ion is given by:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M52','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M52">View MathML</a>

(22)

where [CaC0] is the total concentration of the Ca2+ binding protein and <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M53','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M53">View MathML</a> is the equilibrium dissociation constant.

The Ca2+ binding protein interaction with the Ca2+ channels is given by:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M54','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M54">View MathML</a>

(23)

where C is the number of channels in the open conductive state, and C0 is the total number of channels. The steady-state solution of equation (23) is given by:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M55','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M55">View MathML</a>

(24)

where <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M56','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M56">View MathML</a> is the equilibrium dissociation constant for the Ca2+ binding protein-Ca2+ channels interaction. By combing equations (22) and (24) one obtains the dependence of the open Ca2+ channels on the Ca2+ concentration:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M57','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M57">View MathML</a>

(25)

where <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M58','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M58">View MathML</a>,

Equation (23) can then be represented in the following non dimensional form:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M59','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M59">View MathML</a>

(26)

where <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M60','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M60">View MathML</a>,.

The solution of equation (26) that describes the dynamics of Ca2+ channels in the open conductive state in response to a Ca2+ shift from u0 to u1 Ca2+ level is given by:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M61','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M61">View MathML</a>

(27)

where <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M62','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M62">View MathML</a>,.

The solution (27) suggests that the number of open channels would change exponentially in response to a Ca2+ surge. The characteristic time for such an exponential change is given by:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M63','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M63">View MathML</a>

(28)

For cases when cac0 > > 1 and u > > 1, the characteristic time approximately equals <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M64','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M64">View MathML</a>

Under the assumption that the alteration of the transmembrane potential difference influences the Ca2+ channel conductivity, the channel conductivity as a function of Ca2+ concentration can be represented as follows:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M65','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M65">View MathML</a>

(29)

When membrane potential changes from ψ0 to ψ1, and the alteration of the intraciliary Ca2+ concentration is delayed, the conductivity changes from one value to another exponentially with the characteristic time τV:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M66','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M66">View MathML</a>

(30)

Here we consider a simplified scenario when only passive and active Ca2+ transport channels are present. For such a model, the alteration of Ca2+ concentration in a cilium is given by:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M67','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M67">View MathML</a>

(31)

By substituting formulas for the passive and active Ca2+ currents into equation (31), one obtains:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M68','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M68">View MathML</a>

(32)

In this equation, we include the kinetics for the active Ca2+ channels due to the assumption that the dynamics of currents via the active Ca2+ channels is much faster than the dynamics of currents through the passive Ca2+ transport.

The non dimensional representation of Equations (32) and (23) is given by:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M69','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M69">View MathML</a>

(33)

where <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M70','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M70">View MathML</a>

By substituting the non dimensional representation of the Ca2+ current (21) into (33) the equation for the intraciliary Ca2+ dynamics (33) take the following form:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M71','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M71">View MathML</a>

(34)

where <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M72','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M72">View MathML</a>.

Potassium current

Here we consider the mechanism of K+ current activation by membrane potential. The full K+ current in a cilium is given by:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M73','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M73">View MathML</a>

(35)

where <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M74','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M74">View MathML</a> is the number of open K+ channels, <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M75','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M75">View MathML</a> is the maximal conductivity, and <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M76','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M76">View MathML</a> is the equilibrium K+ potential.

We assume that the K+ channel conductance is dependent on the membrane potential (the opposite of the case of Ca2+ channel conductance). The K+ channel conductance dependence on the membrane potential is given by:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M77','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M77">View MathML</a>

(36)

We further assume that in response to the fast membrane potential alteration from ψ0 to ψ1, the K+ conductivity is changing exponentially with a characteristic time <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M78','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M78">View MathML</a>. In that case the conductivity dynamics over time are given by:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M79','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M79">View MathML</a>

(37)

The parameters for the K+ channels' conductivity dependence on the membrane potential can be estimated from experimental measurements of the fast K+ current (10 msec) as a function of membrane potential. In such a short time K+ current does not reach the maximum value and the Ca2+ current is almost equal to zero. Under such conditions, the equation for the K+ current is given by:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M80','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M80">View MathML</a>

(38)

The full current measured under the voltage clamp conditions is then given by:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M81','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M81">View MathML</a>

(39)

where <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M82','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M82">View MathML</a>.

The transmembrane potential dynamics

All the models described above were specifically developed under the voltage clamp conditions. However, Ca2+ as well as K+ currents can alter the membrane potential. In order to account for this effect, we considered the membrane potential dynamics in the absence of voltage clamp. Having established that the main contributors to the registered currents are the Ca2+ and K+ currents, we analyse the membrane potential dynamics as a function of Ca2+ and K+ current contributions:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M83','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M83">View MathML</a>

(40)

In the non dimensional form the equation for the membrane potential dynamics is given by:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M84','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M84">View MathML</a>

(41)

where <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M85','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M85">View MathML</a>, and <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M86','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M86">View MathML</a> are the steady-state Ca2+ and K+ channel conductivities, respectively, <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M87','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M87">View MathML</a>,.

In this equation, we take into consideration the contribution of the independent parts of the Ca2+ and K+ currents, conductivities <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M88','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M88">View MathML</a> and <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M86','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M86">View MathML</a>, respectively, that do not depend on membrane potential or Ca2+ concentration. We also use the fact that Ca2+ and K+ channel conductivities change on a much faster time scale in comparison with membrane potential or Ca2+ concentration. Such an approximation allows us to employ the steady-state solutions for the channel conductivities as a function of membrane potential and Ca2+ concentration. The described assumptions and considerations lead to the following system of two nonlinear coupled equations for intraciliary Ca2+ concentration and membrane potential:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M89','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M89">View MathML</a>

(42)

where I0 is the non dimensional inward current.

Currents activated by the membrane hyperpolarisation

The membrane hyperpolarisation current is given by:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M90','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M90">View MathML</a>

(43)

where <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M91','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M91">View MathML</a> is the number of open, membrane hyperpolarisation-activated Ca2+ channels located on the cell body, <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M92','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M92">View MathML</a>is the conductivity of the channel, and <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M93','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M93">View MathML</a> is the equilibrium potential for Ca2+ ions.

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M94','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M94">View MathML</a>

(44)

where <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M95','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M95">View MathML</a>

By combining the contribution from different types of currents, one can derive the dynamics of the membrane potential alteration as follows:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M96','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M96">View MathML</a>

(45)

The coupled system of differential equations for the Ca2+ and membrane potential dynamics can be further formulated as:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M97','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M97">View MathML</a>

(46)

where vh(ψ) is the Ca2+ current contribution, activated by membrane depolarization, and <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M98','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M98">View MathML</a> is the Ca2+-dependent K+ current contribution.

Cilia-to body Ca2+ current

The incorporation of the Ca2+ current from cilia to cell into the equations for the intraciliary Ca2+ and membrane potential system leads to an additional term being added to the equations (46):

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M99','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M99">View MathML</a>

(47)

In the general case the conductivity of the protein complexes governing the Ca2+ current from cilia to the cell body depends on the membrane potential:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/143/mathml/M100','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/143/mathml/M100">View MathML</a>

(48)

All parameter values used in the above equations are given in Table 1. The relationship between dimensional and non-dimensional quantities for Ca2+ concentration and membrane potential are given in Table 2.

Table 1. Parameter values employed in the systems model for the ciliary excitation

Table 2. The relationship between dimensional and non-dimensional quantities for Ca2+ concentration and membrane potential

Authors' contributions

NVK, ANG, BG, MZQC, ID, CH and AA developed and implemented the project under the supervision of DGB, RNK, YU and NVV. XY, SKS and CMS contributed to the analysis of the model. All authors contributed to the writing of the final manuscript. All authors read and approved the final manuscript.

Acknowledgements

This work was supported by the Strategic Research Development Award from Faculty of Sciences, University of Kent (NVV) and the Russian Fund for Basic Research (NVK).

References

  1. Ainsworth C: Cilia: tails of the unexpected.

    Nature 2007, 448:638-641. PubMed Abstract | Publisher Full Text OpenURL

  2. Dawe HR, Farr H, Gull K: Centriole/basal body morphogenesis and migration during ciliogenesis in animal cells.

    J Cell Sci 2007, 120:7-15. PubMed Abstract | Publisher Full Text OpenURL

  3. Nigg EA, Raff JW: Centrioles, centrosomes, and cilia in health and disease.

    Cell 2009, 139:663-678. PubMed Abstract | Publisher Full Text OpenURL

  4. Mitchell DR, Nakatsugawa M: Bend propagation drives central pair rotation in Chlamydomonas reinhardtii flagella.

    J Cell Biol 2004, 166:709-715. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  5. Wheatley DN, Wang AM, Strugnell GE: Expression of primary cilia in mammalian cells.

    Cell Biol Int 1996, 20:73-81. PubMed Abstract | Publisher Full Text OpenURL

  6. Afzelius BA: Cilia-related diseases.

    J Pathol 2004, 204:470-477. PubMed Abstract | Publisher Full Text OpenURL

  7. King M, Gilboa A, Meyer FA, Silberberg A: On the transport of mucus and its rheologic simulants in ciliated systems.

    Am Rev Respir Dis 1974, 110:740-745. PubMed Abstract OpenURL

  8. Ma W, Silberberg SD, Priel Z: Distinct axonemal processes underlie spontaneous and stimulated airway ciliary activity.

    J Gen Physiol 2002, 120:875-885. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. Gheber L, Priel Z: Metachronal activity of cultured mucociliary epithelium under normal and stimulated conditions.

    Cell Motil Cytoskeleton 1994, 28:333-345. PubMed Abstract | Publisher Full Text OpenURL

  10. Mason SJ, Paradiso AM, Boucher RC: Regulation of transepithelial ion transport and intracellular calcium by extracellular ATP in human normal and cystic fibrosis airway epithelium.

    Br J Pharmacol 1991, 103:1649-1656. PubMed Abstract | PubMed Central Full Text OpenURL

  11. Ovadyahu D, Eshel D, Priel Z: Intensification of ciliary motility by extracellular ATP.

    Biorheology 1988, 25:489-501. PubMed Abstract OpenURL

  12. Salathe M: Regulation of mammalian ciliary beating.

    Annu Rev Physiol 2007, 69:401-422. PubMed Abstract | Publisher Full Text OpenURL

  13. Salathe M, Lipson EJ, Ivonnet PI, Bookman RJ: Muscarinic signaling in ciliated tracheal epithelial cells: dual effects on Ca2+ and ciliary beating.

    Am J Physiol 1997, 272:L301-310. PubMed Abstract | Publisher Full Text OpenURL

  14. Sanderson MJ, Dirksen ER: Mechanosensitive and beta-adrenergic control of the ciliary beat frequency of mammalian respiratory tract cells in culture.

    Am Rev Respir Dis 1989, 139:432-440. PubMed Abstract | Publisher Full Text OpenURL

  15. Verdugo P: Ca2+-dependent hormonal stimulation of ciliary activity.

    Nature 1980, 283:764-765. PubMed Abstract | Publisher Full Text OpenURL

  16. Villalon M, Hinds TR, Verdugo P: Stimulus-response coupling in mammalian ciliated cells. Demonstration of two mechanisms of control for cytosolic [Ca2+].

    Biophys J 1989, 56:1255-1258. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Weiss T, Gheber L, Shoshan-Barmatz V, Priel Z: Possible mechanism of ciliary stimulation by extracellular ATP: involvement of calcium-dependent potassium channels and exogenous Ca2+.

    J Membr Biol 1992, 127:185-193. PubMed Abstract OpenURL

  18. Wong LB, Yeates DB: Luminal purinergic regulatory mechanisms of tracheal ciliary beat frequency.

    Am J Respir Cell Mol Biol 1992, 7:447-454. PubMed Abstract | Publisher Full Text OpenURL

  19. Zagoory O, Braiman A, Gheber L, Priel Z: Role of calcium and calmodulin in ciliary stimulation induced by acetylcholine.

    Am J Physiol Cell Physiol 2001, 280:C100-109. PubMed Abstract | Publisher Full Text OpenURL

  20. Valeyev NV, Bates DG, Umezawa Y, Gizatullina AN, Kotov NV: Systems biology of cell behavior.

    Methods Mol Biol 2010, 662:79-95. PubMed Abstract | Publisher Full Text OpenURL

  21. Bettencourt-Dias M, Glover DM: Centrosome biogenesis and function: centrosomics brings new understanding.

    Nat Rev Mol Cell Biol 2007, 8:451-463. PubMed Abstract | Publisher Full Text OpenURL

  22. Branche C, Kohl L, Toutirais G, Buisson J, Cosson J, Bastin P: Conserved and specific functions of axoneme components in trypanosome motility.

    J Cell Sci 2006, 119:3443-3455. PubMed Abstract | Publisher Full Text OpenURL

  23. Carvalho-Santos Z, Machado P, Branco P, Tavares-Cadete F, Rodrigues-Martins A, Pereira-Leal JB, Bettencourt-Dias M: Stepwise evolution of the centriole-assembly pathway.

    J Cell Sci 2010, 123:1414-1426. PubMed Abstract | Publisher Full Text OpenURL

  24. Hildebrandt F, Otto E: Cilia and centrosomes: a unifying pathogenic concept for cystic kidney disease?

    Nat Rev Genet 2005, 6:928-940. PubMed Abstract OpenURL

  25. Piasecki BP, Burghoorn J, Swoboda P: Regulatory Factor × (RFX)-mediated transcriptional rewiring of ciliary genes in animals.

    Proc Natl Acad Sci USA 2010, 107:12969-12974. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  26. Satir P, Christensen ST: Structure and function of mammalian cilia.

    Histochem Cell Biol 2008, 129:687-693. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Eckert R, Naitoh Y, Friedman K: Sensory mechanisms in Paramecium. I. Two components of the electric response to mechanical stimulation of the anterior surface.

    J Exp Biol 1972, 56:683-694. PubMed Abstract | Publisher Full Text OpenURL

  28. Murakami A, Eckert R: Cilia: activation coupled to mechanical stimulation by calcium influx.

    Science 1972, 175:1375-1377. PubMed Abstract | Publisher Full Text OpenURL

  29. Naitoh Y, Eckert R, Friedman K: A regenerative calcium response in Paramecium.

    J Exp Biol 1972, 56:667-681. PubMed Abstract | Publisher Full Text OpenURL

  30. Braiman A, Priel Z: Intracellular stores maintain stable cytosolic Ca(2+) gradients in epithelial cells by active Ca(2+) redistribution.

    Cell Calcium 2001, 30:361-371. PubMed Abstract | Publisher Full Text OpenURL

  31. Salathe M, Bookman RJ: Mode of Ca2+ action on ciliary beat frequency in single ovine airway epithelial cells.

    J Physiol 1999, 520(Pt 3):851-865. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  32. Zagoory O, Braiman A, Priel Z: The mechanism of ciliary stimulation by acetylcholine: roles of calcium, PKA, and PKG.

    J Gen Physiol 2002, 119:329-339. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  33. Dong Z, Saikumar P, Weinberg JM, Venkatachalam MA: Calcium in cell injury and death.

    Annu Rev Pathol 2006, 1:405-434. PubMed Abstract | Publisher Full Text OpenURL

  34. Orrenius S, Zhivotovsky B, Nicotera P: Regulation of cell death: the calcium-apoptosis link.

    Nat Rev Mol Cell Biol 2003, 4:552-565. PubMed Abstract | Publisher Full Text OpenURL

  35. Korngreen A, Priel Z: Purinergic stimulation of rabbit ciliated airway epithelia: control by multiple calcium sources.

    J Physiol 1996, 497(Pt 1):53-66. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  36. Braiman A, Gold'Shtein V, Priel Z: Feasibility of a sustained steep Ca(2+)Gradient in the cytosol of electrically non-excitable cells.

    J Theor Biol 2000, 206:115-130. PubMed Abstract | Publisher Full Text OpenURL

  37. Uzlaner N, Priel Z: Interplay between the NO pathway and elevated [Ca2+]i enhances ciliary activity in rabbit trachea.

    J Physiol 1999, 516(Pt 1):179-190. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  38. Geary CA, Davis CW, Paradiso AM, Boucher RC: Role of CNP in human airways: cGMP-mediated stimulation of ciliary beat frequency.

    Am J Physiol 1995, 268:L1021-1028. PubMed Abstract | Publisher Full Text OpenURL

  39. Schmid A, Bai G, Schmid N, Zaccolo M, Ostrowski LE, Conner GE, Fregien N, Salathe M: Real-time analysis of cAMP-mediated regulation of ciliary motility in single primary human airway epithelial cells.

    J Cell Sci 2006, 119:4176-4186. PubMed Abstract | Publisher Full Text OpenURL

  40. Wyatt TA, Forget MA, Adams JM, Sisson JH: Both cAMP and cGMP are required for maximal ciliary beat stimulation in a cell-free model of bovine ciliary axonemes.

    Am J Physiol Lung Cell Mol Physiol 2005, 288:L546-551. PubMed Abstract | Publisher Full Text OpenURL

  41. Yang B, Schlosser RJ, McCaffrey TV: Dual signal transduction mechanisms modulate ciliary beat frequency in upper airway epithelium.

    Am J Physiol 1996, 270:L745-751. PubMed Abstract | Publisher Full Text OpenURL

  42. Zhang L, Sanderson MJ: The role of cGMP in the regulation of rabbit airway ciliary beat frequency.

    J Physiol 2003, 551:765-776. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  43. Hamasaki T, Barkalow K, Richmond J, Satir P: cAMP-stimulated phosphorylation of an axonemal polypeptide that copurifies with the 22S dynein arm regulates microtubule translocation velocity and swimming speed in Paramecium.

    Proc Natl Acad Sci USA 1991, 88:7918-7922. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  44. Hamasaki T, Barkalow K, Satir P: Regulation of ciliary beat frequency by a dynein light chain.

    Cell Motil Cytoskeleton 1995, 32:121-124. PubMed Abstract | Publisher Full Text OpenURL

  45. Satir P, Barkalow K, Hamasaki T: Ciliary beat frequency is controlled by a dynein light chain phosphorylation.

    Biophys J 1995, 68:222S. PubMed Abstract | PubMed Central Full Text OpenURL

  46. Salathe M, Pratt MM, Wanner A: Cyclic AMP-dependent phosphorylation of a 26 kD axonemal protein in ovine cilia isolated from small tissue pieces.

    Am J Respir Cell Mol Biol 1993, 9:306-314. PubMed Abstract | Publisher Full Text OpenURL

  47. Pernberg J, Machemer H: Voltage-dependence of ciliary activity in the ciliate Didinium nasutum.

    J Exp Biol 1995, 198:2537-2545. PubMed Abstract | Publisher Full Text OpenURL

  48. Nakaoka Y, Imaji T, Hara M, Hashimoto N: Spontaneous fluctuation of the resting membrane potential in Paramecium: amplification caused by intracellular Ca2+.

    J Exp Biol 2009, 212:270-276. PubMed Abstract | Publisher Full Text OpenURL

  49. Brehm P, Eckert R: An electrophysiological study of the regulation of ciliary beating frequency in Paramecium.

    J Physiol 1978, 283:557-568. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  50. Iwadate Y, Nakaoka Y: Calcium regulates independently ciliary beat and cell contraction in Paramecium cells.

    Cell Calcium 2008, 44:169-179. PubMed Abstract | Publisher Full Text OpenURL

  51. Naitoh Y: Reversal response elicited in nonbeating cilia of paramecium by membrane depolarizatin.

    Science 1966, 154:660-662. PubMed Abstract | Publisher Full Text OpenURL

  52. Naitoh Y, Eckert R: Ionic mechanisms controlling behavioral responses of paramecium to mechanical stimulation.

    Science 1969, 164:963-965. PubMed Abstract | Publisher Full Text OpenURL

  53. Schultz JE, Guo Y, Kleefeld G, Volkel H: Hyperpolarization- and depolarization-activated Ca2+ currents in Paramecium trigger behavioral changes and cGMP formation independently.

    J Membr Biol 1997, 156:251-259. PubMed Abstract | Publisher Full Text OpenURL

  54. Oami K, Takahashi M: K+-induced Ca2+ conductance responsible for the prolonged backward swimming in K+-agitated mutant of Paramecium caudatum.

    J Membr Biol 2003, 195:85-92. PubMed Abstract | Publisher Full Text OpenURL

  55. Oami K, Takahashi M: Identification of the Ca2+ conductance responsible for K+-induced backward swimming in Paramecium caudatum.

    J Membr Biol 2002, 190:159-165. PubMed Abstract | Publisher Full Text OpenURL

  56. Gheber L, Priel Z: Synchronization between beating cilia.

    Biophys J 1989, 55:183-191. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  57. Gheber L, Korngreen A, Priel Z: Effect of viscosity on metachrony in mucus propelling cilia.

    Cell Motil Cytoskeleton 1998, 39:9-20. PubMed Abstract | Publisher Full Text OpenURL

  58. Hilfinger A, Julicher F: The chirality of ciliary beats.

    Phys Biol 2008, 5:016003. PubMed Abstract | Publisher Full Text OpenURL

  59. Castillo K, Bacigalupo J, Wolff D: Ca2+-dependent K+ channels from rat olfactory cilia characterized in planar lipid bilayers.

    FEBS Lett 2005, 579:1675-1682. PubMed Abstract | Publisher Full Text OpenURL

  60. Doughty MJ, Dryl S: Control of ciliary activity in Paramecium: an analysis of chemosensory transduction in a eukaryotic unicellular organism.

    Prog Neurobiol 1981, 16:1-115. PubMed Abstract | Publisher Full Text OpenURL

  61. Iwadate Y, Suzaki T: Ciliary reorientation is evoked by a rise in calcium level over the entire cilium.

    Cell Motil Cytoskeleton 2004, 57:197-206. PubMed Abstract | Publisher Full Text OpenURL

  62. Kawai F: Ca2+-activated K+ currents regulate odor adaptation by modulating spike encoding of olfactory receptor cells.

    Biophys J 2002, 82:2005-2015. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  63. Valeyev NV, Bates DG, Heslop-Harrison P, Postlethwaite I, Kotov NV: Elucidating the mechanisms of cooperative calcium-calmodulin interactions: a structural systems biology approach.

    BMC Syst Biol 2008, 2:48. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  64. Valeyev NV, Heslop-Harrison P, Postlethwaite I, Kotov NV, Bates DG: Multiple calcium binding sites make calmodulin multifunctional.

    Mol Biosyst 2008, 4:66-73. PubMed Abstract | Publisher Full Text OpenURL

  65. Valeyev NV, Heslop-Harrison P, Postlethwaite I, Gizatullina AN, Kotov NV, Bates DG: Crosstalk between G-protein and Ca2+ pathways switches intracellular cAMP levels.

    Mol Biosyst 2009, 5:43-51. PubMed Abstract | Publisher Full Text OpenURL

  66. Eckert R, Brehm P: Ionic mechanisms of excitation in Paramecium.

    Annu Rev Biophys Bioeng 1979, 8:353-383. PubMed Abstract | Publisher Full Text OpenURL

  67. Bannai H, Yoshimura M, Takahashi K, Shingyoji C: Calcium regulation of microtubule sliding in reactivated sea urchin sperm flagella.

    J Cell Sci 2000, 113(Pt 5):831-839. PubMed Abstract | Publisher Full Text OpenURL

  68. Hayashi S, Shingyoji C: Bending-induced switching of dynein activity in elastase-treated axonemes of sea urchin sperm--roles of Ca2+ and ADP.

    Cell Motil Cytoskeleton 2009, 66:292-301. PubMed Abstract | Publisher Full Text OpenURL

  69. Nakano I, Kobayashi T, Yoshimura M, Shingyoji C: Central-pair-linked regulation of microtubule sliding by calcium in flagellar axonemes.

    J Cell Sci 2003, 116:1627-1636. PubMed Abstract | Publisher Full Text OpenURL

  70. Tamm S: Ca/Ba/Sr-induced conformational changes of ciliary axonemes.

    Cell Motil Cytoskeleton 1990, 17:187-196. PubMed Abstract | Publisher Full Text OpenURL

  71. Tamm SL: Control of reactivation and microtubule sliding by calcium, strontium, and barium in detergent-extracted macrocilia of Beroe.

    Cell Motil Cytoskeleton 1989, 12:104-112. PubMed Abstract | Publisher Full Text OpenURL

  72. Doughty MJ: Control of ciliary activity in Paramecium--II. Modification of K+-induced ciliary reversal by cholinergic ligands and quaternary ammonium compounds.

    Comp Biochem Physiol C 1978, 61C:375-384. OpenURL

  73. Doughty MJ: Control of ciliary activity in Paramecium--I. Modification of K+-induced ciliary reversal by temperature and ruthenium red.

    Comp Biochem Physiol C 1978, 61C:369-373. OpenURL

  74. Rudy Y, Silva JR: Computational biology in the study of cardiac ion channels and cell electrophysiology.

    Q Rev Biophys 2006, 39:57-116. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  75. Valeyev NV, Downing AK, Skorinkin AI, Campbell ID, Kotov NV: A calcium dependent de-adhesion mechanism regulates the direction and rate of cell migration: a mathematical model.

    In Silico Biol 2006, 6:545-572. PubMed Abstract OpenURL

  76. Bezprozvanny I, Watras J, Ehrlich BE: Bell-shaped calcium-response curves of Ins(1, 4, 5)P3- and calcium-gated channels from endoplasmic reticulum of cerebellum.

    Nature 1991, 351:751-754. PubMed Abstract | Publisher Full Text OpenURL

  77. Frolenkov GI: Regulation of electromotility in the cochlear outer hair cell.

    J Physiol 2006, 576:43-48. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  78. Mogami Y, Baba SA: Super-helix model: a physiological model for gravitaxis of Paramecium.

    Adv Space Res 1998, 21:1291-1300. PubMed Abstract | Publisher Full Text OpenURL

  79. Satir P, Barkalow K, Hamasaki T: The control of ciliary beat frequency.

    Trends Cell Biol 1993, 3:409-412. PubMed Abstract | Publisher Full Text OpenURL

  80. Nakaoka Y, Tanaka H, Oosawa F: Ca2+-dependent regulation of beat frequency of cilia in Paramecium.

    J Cell Sci 1984, 65:223-231. PubMed Abstract | Publisher Full Text OpenURL

  81. Abdul-Majeed S, Nauli SM: Dopamine receptor type 5 in the primary cilia has dual chemo- and mechano-sensory roles.

    Hypertension 2011, 58:325-331. PubMed Abstract | Publisher Full Text OpenURL

  82. Abdul-Majeed S, Moloney BC, Nauli SM: Mechanisms regulating cilia growth and cilia function in endothelial cells.

    Cell Mol Life Sci 2011. OpenURL

  83. Nauli SM, Haymour HS, WA A, Lo ST, Nauli AM: Primary Cilia are Mechanosensory Organelles in Vestibular Tissues. In Mechanosensitivity and Mechanotransduction. Mechanosensitivity in Cells and Tissues. Volume 4. Edited by Kamkin A, Kiseleva I. Dordrecht Heidelberg London New York: Springer; 2011::317-350. OpenURL

  84. Brady AJ, Tan ST: The ionic dependence of cardiac excitability and contractility.

    J Gen Physiol 1966, 49:781-791. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  85. Einwachter HM, Haas HG, Kern R: Membrane current and contraction in frog atrial fibres.

    J Physiol 1972, 227:141-171. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  86. Kobatake Y, Tasaki I, Watanabe A: Phase transition in membrane with reference to nerve excitation.

    Adv Biophys 1971, 2:1-31. PubMed Abstract OpenURL

  87. Sah R, Ramirez RJ, Oudit GY, Gidrewicz D, Trivieri MG, Zobel C, Backx PH: Regulation of cardiac excitation-contraction coupling by action potential repolarization: role of the transient outward potassium current (I(to)).

    J Physiol 2003, 546:5-18. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  88. Braiman A, Priel Z: Efficient mucociliary transport relies on efficient regulation of ciliary beating.

    Respir Physiol Neurobiol 2008, 163:202-207. PubMed Abstract | Publisher Full Text OpenURL

  89. Valeyev NV, Hundhausen C, Umezawa Y, Kotov NV, Williams G, Clop A, Ainali C, Ouzounis C, Tsoka S, Nestle FO: A systems model for immune cell interactions unravels the mechanism of inflammation in human skin.

    PLoS Comput Biol 2010, 6:e1001024. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  90. Satir P, Christensen ST: Overview of structure and function of mammalian cilia.

    Annu Rev Physiol 2007, 69:377-400. PubMed Abstract | Publisher Full Text OpenURL

  91. Hook C, Hildebrand E: Excitation of paramecium.

    Journal of Mathematical Biology 1979, 8:197-214. Publisher Full Text OpenURL