## Services on Demand

## Article

## Indicators

- Cited by SciELO
- Access statistics

## Related links

- Cited by Google
- Similars in SciELO
- Similars in Google

## Share

## Computational & Applied Mathematics

*On-line version* ISSN 1807-0302

### Comput. Appl. Math. vol.31 no.2 São Carlos 2012

#### http://dx.doi.org/10.1590/S1807-03022012000200004

**The sodium pump controls the frequency of action-potential-induced calcium oscillations**

**Shivendra G. Tewari**^{*}

Systems Science and Informatics Unit, Indian Statistical Institute, 8^{th} Mile, Mysore Road, Bangalore 560059, India E-mail: tewarisg@gmail.com

**ABSTRACT**

Calcium plays a significant role in a number of cellular processes, like muscle contraction, gene expression, synaptic plasticity, signal transduction, but the significance of calcium oscillations (CaOs) is not yet completely understood in most of the cell types. It is a widely accepted fact that CaOs are a frequency encoded signal that allows a cell to use calcium as a second messenger while avoiding its toxic effects. These intracellular CaOs are primarily driven by some agonist-dependent pathways or fluctuations in membrane potential. The present mathematical model is of the latter type. The model incorporates expression for all major intracellular ionic species and membrane proteins. Especially, it integrates the coupling effect of sodium pump and Na^{+} / Ca^{2+} exchanger over CaOs. By varying sodium pump current, it is found that, sodium pump is a key player in modulating intracellular CaOs. The model predicts that the sodium pump can play a decisive role in regulating intercellular cell signaling process. The present study forms the basis for sodium pump controlled intercellular signaling process and requires further experimental verification.

**Mathematical subject classification:** 34M10, 92C20.

**Key words:** Na^{+}/Ca^{2+} exchanger, sodium pump, calcium oscillations, membrane potential.

**1 Introduction**

But till date the function of cytosolic Calcium oscillations (CaOs) has not been completely understood in most cell types. CaOs are known to play a key role ina number of mechanisms like activation of extracellular signal regulated kinase (ERK) [4, 5], the contraction of smooth muscle [6], increase in the frequency of synaptic currents [7] and maturation of *Xenopus laevis Oocyte* [8]. These CaOs are supposed to contain frequency-encoded signals that help in using Ca^{2+} as a second messenger while avoiding its high intracellular concentrations [9]. Also, in the process of signal transduction, intracellular Ca^{2+} behaves like a switch and decides whether a particular signal needs to be further propagated or not. The increase in intracellular concentration is facilitated by the opening of transmembrane Ca^{2+} channels which lead to the opening of channels at the intracellular stores. There are mainly Rynodine Receptors (RyRs) or Inositol Triphosphate Receptors (IP_{3}Rs) that are located at the membrane of endoplasmic reticulum (in neurons) or sarcoplasmic reticulum (in myocytes) which causes an efflux of Ca^{2+} from the intracellular stores. The release of Ca^{2+} through IP_{3}Rs is as a result of some agonist or neurotransmitter binding to its receptor which can cause via G-protein link to phospholipase C (PLC), the cleavage of phosphotidylinositol (4,5)-bisphosphate (PIP_{2}) to inositol triphosphate (IP_{3}) and diacylglycerol (DAG). This released IP_{3} is free to diffuse through the cytosol and binds with IP_{3}R and leading to the subsequent opening of these receptors and release of Ca^{2+} from the intracellular stores. CaOs can be classified into mainly two types:

1) that is induced by changing membrane potential as in the case of an action potential; and

2) that occur in the presence of voltage clamp. The latter part can be further categorized based on the fact that the oscillatory Ca

^{2+}flux is from RyRs or IP_{3}Rs but our focus is on the first type.

Jafri et al. [6] showed CaOs for changing membrane potential of endoplasmic reticulum (ER), Atri et al. [10] showed CaOs for Ca^{2+} flux through the IP_{3}Rs in Xenopus laevis Oocyte and determined the intermediate range of IP_{3} for CaOs. Wagner and Keizer [11] showed the effect of rapid buffering over CaOs. Later on, Kusters et al. [12] proposed an integrated model which combines excitable membrane with an IP_{3} mediated Ca^{2+} oscillator for normal rat kidney (NRK) fibroblast. Recently, Silva et al. [13] proposed a mathematical model for endothelial cells which incorporated nearly all the important biophysical parameters but is unable to exhibit CaOs. Thus, we can say that in all the investigations on CaOs carried out by researchers so far, none of the investigators have tried to incorporate the effect of changing cytosolic Na^{+} and K^{+} ions over CaOs. Thus, in this article, we have proposed a mathematical model governing CaOs for changing membrane potential in the absence of fluxes from the intracellular stores. Holmgren et al. [14] determined the three distinct steps of Na^{+} ion release from the sodium pump with the help of high speed voltage jumps. Here, we have tried to incorporate the impact of these distinct steps of sodium pump over cytosolic CaOs, in case of an action potential. Thus, we have incorporated L-type Ca^{2+} channel, Na^{+} channel, K^{+} channel, Plasma-Membrane (PM) Ca^{2+} ATPase, Na^{+} / Ca^{2+} exchanger (NCX), Na^{+} / K^{+} ATPase (Sodium pump), inward rectifier potassium channel (K_{ir}), Ca^{2+} dependent intermediate potassium channel (IK_{Ca}, Ca^{2+} dependent small potassium channel (SK_{Ca}) and dynamic membrane potential. The gating mechanism of the trasmembrane channels emulate the gating mechanism of the famous Hodgkin and Huxley model [15]. Further, the pumps and proteins are modeled to have realistic gating mechanism such that they are in agreement with the biological facts. The proposed mathematical model leads to a system of non-linear ordinary differential equations. We have used Euler's method for the simulation of the proposed model for which a MATLAB script has been written.

**2 Mathematical formulation**

Our cell model assumes that cell is cylindrical in shape. The diameter of cell is assumed to be 20 *µ*m and its length to be 100 *µ*m (see Fig. 1). Its specific membrane capacitance is taken to be 1 *µ*F/cm^{2}. As mentioned in literature [16, 17, 18], we have taken actual surface area (A_{cap}) of cell to be larger than its geometrical surface area. Further, the channels, proteins and pumps are supposed to be homogenously distributed throughout the membrane. The formulation of the proposed mathematical model comprises of different components which are elaborated in the following subsections:

*2.1 The Ca ^{2+} channel*

The L-type Ca^{2+} channel is supposed to have a permeability ratio of 36000:1:18 for Ca^{2+}, Na^{+} and K^{+} ions, respectively [19]. The voltage gated Ca^{2+} channel was modeled using non-linear Goldman-Hodgkin-Katz (GHK) current equation [9, 20], which can be stated as,

where, *S* is any of the ions, [*S*]_{i} and [*S*]_{o} are the intracellular and extracellular concentration of*S*ion respectively (in mM), *P _{S}* is permeability (in cm/s) of

*S*ion,

*z*is its valence,

_{S}*γ*

_{i},

*γ*

_{o}are the activity coefficient (a.c.) of the

*S*ion,

*F*is Faradays constant (in Coulombs/moles),

*V*is membrane potential (in volts),

_{m}*R*is real gas constant (in J/K moles) and

*T*is absolute temperature (in K). The total current through the L-type Ca

^{2+}channel is

*I*

_{Ca,t}= I_{Ca,L}+I_{Na,L}+I_{K,L}.Further, equation (1) is converted into fluxes (in mM/second), before being used as an expression for individual ionic concentrations, by using Faraday's constant, Volume of the cytosol (V_{cyt}) and using the fact that 1 L = 10^{-3} m^{3},

where, *S* is any of the ion and *L* signifies *L*-type Ca^{2+} channel. In equation (2) there is a negative sign against it because by convention inward current is taken to be negative.

*2.2 PM Ca ^{2+} ATPase*

PM Ca^{2+} ATPase (PMCA) is a P-type ATPase. The energy required to extrude Ca^{2+} out of the cytosol is met by ATP. The kinetics of this pump follows the enzyme-substrate formulism and hence using *Michaelis Menten* type kinetics [21], one can formulate the net efflux of Ca^{2+} ions out of the cytosol by,

where, *Î*_{pump} is the pump current given by the following equation,

where, [Ca^{2+}] is the intracellular Ca^{2+} concentraton (in mM), H,pump is the Hill's coefficient for PMCA, is the Ca^{2+} concentration at which the maximum pump current is halved (in mM).

*2.3 Na ^{+} / Ca^{2+} exchanger*

This protein is known to play an important role in excitation-contraction coupling in cardiac myocytes [22]. In neurons this protein helps in the extrusion of cytosolic Ca^{2+} concentration and hence helps in the modulation of neurotransmitter release [23]. It is known that *the cardiac type* 3Na^{+} / 1Ca^{2+} exchanger is dominant in brain [23]. Thus, we have used the same exchange for our model. We know that the amount of energy required to extrude an ion against its concentration gradient is given by [21, 20, 9],

where, *S* is the extruded ion. Introducing energy barrier, *η*, and using the fact that

_{Ca}= 3 Δ

_{Na}

we can write NCX current equation with an allosteric dependence over [Ca^{2+}] as [13],

where, [*Ca*^{2+}]_{o} is the extracellular Ca^{2+} concentration (in mM), [*Na*^{+}] is the intracellular Na^{+} concentration (in mM), [*Na*^{+}]_{o} is the extracellular Na^{+} concentration (in mM), is the Ca^{2+} concentration at which *I _{NCX}* is halved,

*H, N C X*is the Hill's coefficient of NCX,

*d*is constant for saturability of

_{NCX}*I*,

_{NCX}*g*is the conductance of

_{NCX}*NCX*(in nS).

*2.4 Na ^{+} / K^{+} ATPase*

Na^{+} / K^{+} ATPase (NaK) is also a P-type ATPase which is also known as the *sodium pump* and is a 147 kDa membrane protein [24]. It is known for the extrusion of Na^{+} ions at an expense of some ATP and inflow of K^{+} ions. Its formulation is based on the steps given by Holmgren et al. [14], where heused high speed voltage jumps to determine three distinct steps of Na^{+} ions deocclusion from the pump. The current through the sodium pump has the following form,

here, *I _{NaK}* is the scaling factor of NaK current (in

*µ*A/cm

^{2}),

*k*(in ms) is the forward (deocclusion) rate constant,

_{f}*k*(in ms) is the backward (occlusion) rate constant,

_{b}*K*

_{0.5}(0) is half activating [

*Na*

^{+}]

_{o}concentration at 0 mV,

*H, N a K*is the Hill's coefficient for half activating NaK current, λ is the fraction of electrical field dropped along the access channel and τ

_{NaK}(in ms) is someconstant.

*2.5 Cytosolic Ca ^{2+} buffers*

It is assumed that a single buffer specie is present inside the cytosol and follows the following bi-molecular reaction

and which can be formulated in terms of the following differential equations,

If we also assume that there are no sources and no sinks present for buffer. Then letting *B _{T}* represent the total buffer concentration, the above equation can be written in the reduced form as,

where, *k*^{+} is the buffer association rate, *k*^{-} is the buffer dissociation rate, [*C a B*] represents bound buffer concentration.

*2.6 Na ^{+} and K^{+} channels*

To generate action potentials Na^{+} current and K^{+} channels are taken as modeled by Hodgkin and Huxley [15]. The transmembrane current due to Na^{+} and K^{+} channels has been modeled using the linear current-voltage relationship derived with the help of *Ohm's Law*,

where,* S *is either Na^{+} or K^{+} ion, *g _{S}* is conductance of the given ion,

*V*is the reversal potential of the given ion determined by

_{S}*Nernst Equilibrium Potential*equation (or simply Nernst equation),

where, 1000 is used to convert volts into milli-volts. All other symbols have their usual meanings. Here, and in all other instances, individual ionic reversal potential has been determined using Nernst equation at each integration step during runtime.

*2.7 Ca ^{2+} activated small and intermediate K^{+} channel*

The current through SK_{Ca} is modeled using a linear current voltage relation as follows,

Here, *g _{SKCa}* is SK

_{Ca}channel conductance per unit area (in mS/cm

^{2}),

*P*is its Ca

_{o,SKCa}^{2+}dependent open probability,

*V*is the reversal potential of

_{K}*K*

^{+}ions. Similarly, Ca

^{2+}activated intermediate K

^{+}current is modeled using a linear current voltage relation as follows,

Here, *g _{IKCa}* is

*I K*channel conductance per unit area (in mS/cm

_{Ca}^{2}),

*P*is its Ca

_{o, I KCa}^{2+}dependent open probability,

*V*is the reversal potential of K

_{K}^{+}ions.

*2.8 Inward rectifier K ^{+} channel*

Inward rectifier K^{+} current is known to contribute to resting membrane potential. It is also modeled using the linear current voltage relationship,

where,

and *g _{Kir}* is the K

_{ir}conductance (in nS),

*V*is the reversal potential for K

_{K}^{+}given by Nernst equation, [

*K*

^{+}]

_{o}is the extracellular K

^{+}ion concentration. Here, g

_{Kir}is converted into mS/cm

^{2}using A

_{cap}before being used in the equation governing membrane potential.

*2.9 Ca ^{2+} and Na^{+} leak currents*

To balance the net effect of *I _{NCX}* and

*I*there is supposed to be a Ca

_{pump}^{2+}leak current given by,

*I*(

_{Ca,b}= g_{Ca,b}*V*)

_{m}— V_{Ca}where, *g _{Ca,b}* is the Ca

^{2+}leak conductance per unit area (in mS/cm

^{2}),

*V*is the reversal potential (in mV) for Ca

_{Ca}^{2+}given by Nernst equation. Similarly, we can formulate Na

^{+}leak current to balance the net effect of

*I*and

_{NCX}*Î*,

_{NaK}*I*(

_{Na,b}= g_{Na,b}*V*)

_{m}— V_{Na}where, *g _{Na,b}* is the Na

^{+}leak conductance per unit area (in mS/cm

^{2}),

*V*is the reversal potential (in mV) for Na

_{Na}^{+}given by Nernst equation. The current due to all other ions is considered as leak and is incorporated as,

where, *g _{L}* is leak conductance (in mS/cm

^{2}) and

*V*is leak reversal potential assumed to be constant.

_{L}*2.10 Membrane potential*

Like the formulation of Hodgkin and Huxley [15] we have divided the total membrane current into capacitive current and ionic currents. Thus, for capacitive current we have,

where, *I _{app}* is the applied membrane current density (in

*µ*A/cm

^{2}),

*V*is the membrane potential (in mV),

_{m}*C*is the specific membrane capacitance (

_{m}*µ*F/cm

^{2}), I

_{i}accounts for all the transmembrane currents discussed earlier and

*t*is time (in ms). The gating mechanism of the transmembrane currents follows Hodgkin and Huxley [15]. Combining equation (1)-(13) we can write the mathematical model governing CaOs with relevance to Na

^{+}and K

^{+}ions, as in the case of an action potential as,

In equation (14),

*α*

_{i},

*β*

_{i}(

*i = m, n, h, m*)

_{c}, h_{c}are rate constants which vary with membrane potential but not with time (*ms*^{-1}) and *m, n,h, m _{c}, h_{c}* are dimensionless gating variables with values lying between 0 and 1. In this model we assume that fluxes from IP

_{3}Rs are absent. This can be achieved by blocking IP

_{3}R channel by using an IP

_{3}R antagonist like heparin [25]. This assumption has been taken to exclude the effect of intracellular stores over CaOs. The initial condition of the system is,

All ionic concentrations are in the units of mM. The ordinary differential equations governing the gating variables (*m, n, h, m _{c}, h_{c}*) are

Here, *m* and (1 — * m*) are representing *on* and *off* state of the variable *m*, respectively. Variables *n, h, m _{c}* and

*h*also follow likewise. The mathematical expressions of the voltage-dependent rate constants in equation (14) are as follows,

_{c} For the solution of equations (14)-(16), we have used Euler's method and written a script in MATLAB that has been simulated on an AMD Turion 64 × 2 machine with 1.6 GHz processing speed and 2.5 GB memory. The time taken per simulation is ~9 sec when simulating for 30 ms using 4000 time steps i.e. Δ*t* = 0.0075 ms. The numerical results obtained are used to study the effect of varying transmembrane currents over CaOs which are discussed in the following section.

**3 Results and discussion**

In all the figures, it is assumed that cytosolic Ca^{2+} is buffered with 50 *µ*M Ethylene Glycol Tetraacetic Acid (EGTA). The standard biophysical parameters used for simulation of the model are listed in Table 1-4 unless stated along with the figures. Since our main objective is to study CaOs, we have shown results that are pertinent to CaOs only.

In Figure 2 we observe the effect of an impulse of 10 *µ*A/cm^{2} over membrane potential. Such an effect has been studied in great detail by Hodgkin and Huxley [15], Luo and Rudy [16, 17, 18] thus, we need not give much emphasis over it here. In Figure 3, we have shown different current densities. All these current densities result in the shown action potential. The rest of the results shown are relevant to CaOs and has been studied in detail in the following figures.

In Figure 4 we have shown Ca^{2+} oscillation and buffered Ca^{2+} curves with respect to standard parameters listed in Table 1. Initially it was assumed that 2.3 *µ*M of Ca^{2+} is buffered. Further it is apparent from Figure 4 that the amplitude of first Ca^{2+} spike is more than the second spike. It is because of the slow dissociation constant of EGTA which results in higher buffered Ca^{2+} and lower cytosolic Ca^{2+} concentration. Comparing Figure 4 with Figure 2, it is clear that both of them are positively correlated. As membrane potential rises the Ca^{2+} concentration rises and when membrane potential drops Ca^{2+} concentration drops.

In Figure 5 the effect of increasing total buffer concentration of EGTA isshown. The results are shown for *B _{T}* = 50

*µ*M (dark line) and

*B*= 200

_{T}*µ*M (broken line). As expected increasing buffer concentration results in lower amplitude of Ca

^{2+}oscillation, which is also evident from Figure 5.

In Figure 6 we observe the effect of increasing NCX conductance. Our simulation is in support of the biological fact that at negative potentials NCX works in reverse direction i.e. outflow of Ca^{2+} ions and inflow of 3Na^{+} ions. As increasing NCX conductance results in higher amplitude of Ca^{2+} oscillation while there is no change resting Ca^{2+} concentration at more positive membrane potentials as we have used a leak to neutralize the effect of NCX and pump currents.

In Figure 7 we observe the results for which the mathematical model was proposed. It is widely believed that cells encode information in the frequency of Ca^{2+} oscillations rather than its amplitude. There are a number of authors who have shown different roles of this ubiquitous sodium pump [26, 27]. Matchkov et al. [27] also experimentally demonstrated that sodium pump plays a significant role in regulating CaOs via regulation of cytosolic Na^{+} ions. Similar philosophy is suggested by our present simulations. We increased pumping rate, *I _{NaK} *= 1.5

*µ*A/cm

^{2}(dark line), 3

*µ*A/cm

^{2}(broken line), 4.5

*µ*A/cm

^{2}(dotted line), of sodium pump and observed changes in CaOs. It is seen that an increase in I

_{NaK}results in an increase in the period of Ca

^{2+}oscillations. The changes were quite apparent and are reflected from Figure 7.

In Figures 8-10, the effect of different extracellular concentrations of Ca^{2+}, K^{+} and Na^{+} are shown over CaOs. The findings of Figure 8 are quite obvious but should be mentioned to show accordance with the biological facts. The curves are shown for [Ca^{2+}]_{o} = 1.8 mM (dark line) and [Ca^{2+}]_{o} = 1 mM (broken line). It is apparent from Figure 8 that lowering [Ca^{2+}]_{o} results in lower amplitude of Ca^{2+} oscillation which is because of a corresponding decrease in Ca^{2+} gradient.

In Figure 9, the findings are worth mentioning as an increase in [K^{+}]_{o} concentration resulted in an increase in frequency of Ca^{2+} oscillation. Although these findings are also pretty obvious as increasing [K^{+}]_{o} leads to an increase in reversal potential of K^{+} ions and hence increases the frequency of action potential which in turn increases the frequency of CaOs. In Figure 10, the effect of decreasing [*Na*^{+}]_{o} concentration is shown. The curves are shown for [*Na*^{+}]_{o} = 145 mM (dark line) and [*Na*^{+}]_{o} = 140 mM (broken line), it is apparent from figure that decreasing [*Na*^{+}]_{o} results in an increase in amplitude of CaOs. As in the previous case changing [*Na*^{+}]_{o} concentration changes the reversal potential of Na^{+} ion. But the observed change in *V _{Na}* is minimal and obviously does not affect amplitude of action potential. The reason behind the increase in amplitude and latency in Ca

^{2+}oscillation is the change in Na

^{+}ion gradient. This gradient regulates the pumping rate of NCX exchanger decreasing the gradient means decreasing the pumping rate of NCX exchanger. Hence, affecting the net extrusion of Ca

^{2+}ions via NCX exchanger; resulting in an increase in amplitude and latency of CaOs.

The results obtained in this paper are new and are subject to CaOs. The intent behind the present study was to investigate the effect of Na^{+} / K^{+} ATPase over Ca^{2+} oscillation influenced by the experimental results obtained by Matchkov et al. [27]. The results obtained by our simulations are quite convincing with biological facts. The obtained results also confirmed the hypothesis of Matchkov et al. [27] that interaction between NCX and Na^{+} / K^{+} ATPase modulates intercellular communication. It was observed that increasing NaK current decreases the frequency of CaOs. The results obtained by previous investigators regarding CaOs have been mainly concerned with membrane potential, inositol triphosphate (IP_{3}) or ryanodine receptor [6, 28, 29, 10, 11, 12, 13]. None of the earlier investigators gave much emphasis over this interaction of NCX and sodium pump which in turn effects CaOs. Thus, in this article, we have looked into and demonstrated a novel mechanism which modulates frequency of Ca^{2+} oscillation. Here, we have proposed a mathematical model which can be used for problems related to similar cell processes. The results obtained in this paper give new and useful insight for neurologists to look into the paradigm of CaOs at a different perspective. Also, the results obtained are relevant to biomedical scientists for developing protocols for diagnosis and treatment of neurological disorders.

**Acknowledgments. ** The author acknowledges fruitful discussions with Dr. Ronald J. Clarke, School of Chemistry, The University of Sydney, Australia for giving useful insights over the kinetics of Na^{+} / K^{+} ATPase.

**References**

[1] A.C. Charles, C.C.G. Naus, D. Zhu, G.M. Kidder, E.R. Dirksen and M.J. Sanderson, *Intercellular Calcium Signaling via Gap Junctions in Glioma Cells.* The Journal of Cell Biology, **118** (1992), 195-201. [ Links ]

[2] A. Peskoff and G.A. Langer, *Calcium Concentration and Movement in the Ventricular Cardiac Cell during an Excitation-Contraction Cycle*. Biophys. J., **74** (1998), 153-174. [ Links ]

[3] J. Shuai, J.E. Pearson and I. Parker, *Modeling Ca ^{2+} Feedback on a Single Inositol 1,4,5-Trisphosphate Receptor and Its Modulation by Ca^{2+} Buffers*. Biophys. J.,

**95**(2008), 3738-3752. [ Links ]

[4] O. Melien, L.S. Nilssen, O.F. Dajani, K.L. Sand, J-G Iversen, D.L Sandnes and T. Christoffersen, *Ca ^{2+}-mediated activation of ERK in hepatocytes by norepinephrine and prostaglandin F2 role of calmodulin and src kinases*. BMC Cell Biol.,

**3**(2002). [ Links ]

[5] C.J. Dixon, J.F. Hall, T.E. Webb and M.R. Boarder, *Regulation of Rat Hepatocyte Function by P2Y Receptors: Focus on Control of Glycogen Phosphorylase and Cyclic AMP by 2-Methylthioadenosine 5-Diphosphate*. The Journal of Pharmacology and Experimental Therapeutics, **311** (2004), 334-341. [ Links ]

[6] M.S. Jafri, S.P. Vajda, S. Pasik and B. Gillo, *A membrane model for cytosolic calcium oscillations: A study using Xenopus oocytes*. Biophys. J., **63** (1992), 235-246. [ Links ]

[7] T.A. Fiacco and K.D. McCarthy, *Intracellular Astrocyte Calcium Waves In SituIncrease the Frequency of Spontaneous AMPA Receptor Currents in CA1 Pyramidal Neurons*. The Journal of Neuroscience, **24** (2004), 722-732. [ Links ]

[8] L. Sun, R. Hodeify, S. Haun, A. Charlesworth, A.M. MacNicol, S. Ponnappan, U. Ponnappan, C. Prigent and K. Machaca, *Ca ^{2+} Homeostasis Regulates Xenopus Oocyte Maturation*. Biology of Reproduction,

**78**(2008), 726-735. [ Links ]

[9] J. Keener and J. Sneyd, *Mathematical Physiology*. Springer, **8** (1998). [ Links ]

[10] A. Atri, J. Amundson, D. Clapham and J. Sneyd, *A Single-Pool Model for Intracellular Calcium Oscillations and Waves in the Xenopus laevis Oocyte*. Biophys. J., **65** (1993), 1727-1739. [ Links ]

[11] J. Wagner and J. Keizer, *Effects of Rapid Buffers on Ca ^{2+} Diffusion and Ca^{2+} Oscillations*. Biophys. J.,

**67**(1994), 447-456. [ Links ]

[12] J.M.A.M. Kusters, M.M. Dernison, W.P.M. van Meerwijk, D.L. Ypey, A.P.R. Theuvenet and C.C.A.M. Gielen, *Stabilizing Role of Calcium Store-Dependent Plasma Membrane Calcium Channels in Action-Potential Firing and Intracellular Calcium Oscillations*. Biophys. J., **89** (2005), 3741-3756. [ Links ]

[13] H.S. Silva, A. Kapela and N.M. Tsoukias, *A mathematical model of plasma membrane electrophysiology and calcium dynamics in vascular endothelial cells*. Am. J. Physiol. Cell. Physiol., **293** (2007), C277-C293. [ Links ]

[14] M. Holmgren, J. Wagg, F. Bezanilla, R.F. Rakowski, P. De. Weer and D.C. Gadsby, *Three distinct and sequential steps in the release of sodium ions by the Na+ / K+ ATPase*. Nature, **403** (2000), 898-901. [ Links ]

[15] A.L. Hodgkin and A.F. Huxley, *A Quantitative Description of Membrane Current and its Application to Conduction and Excitation in Nerve*. J. Physiol., **117** (1952), 500-544. [ Links ]

[16] C.H. Luo and Y. Rudy, *A model of the ventricular cardiac action potential. Depolarization, repolarization, and their interaction*. Circ. Res., **68** (1991), 1501-1526. [ Links ]

[17] C.H. Luo and Y. Rudy, *A dynamic model of the cardiac ventricular action potential. I. Simulations of ionic currents and concentration changes*. Circ. Res., **74** (1994), 1071-1096. [ Links ]

[18] C.H. Luo and Y. Rudy, *A dynamic model of the cardiac ventricular action potential. II. After depolarizations, triggered activity, and potentiation*. Circ. Res., **74** (1994), 1097-1113. [ Links ]

[19] T.R. Shannon, F. Wang, J. Puglisi, C. Weber and D.M. Bers, *A Mathematical Treatment of Integrated Ca ^{2+} Dynamics Within the Ventricular Myocyte*. Biophys. J.,

**87**(2004), 3351-3371. [ Links ]

[20] G.L. Fain, *Moleculer and cellular physiology of neurons*. Harvard University Press (1999). [ Links ]

[21] D.L. Nelson and M.M. Cox, *Lehninger Principles of Biochemistry*. W.H. Freeman (2005). [ Links ]

[22] Y. Fujioka, K. Hiroe and S. Matsuoka, *Regulation kinetics of Na+-Ca ^{2+} exchange current in guinea-pig ventricular myocytes*. J. Physiol.,

**529**(2000), 611-623. [ Links ]

[23] M.P. Blaustein and W.J. Lederer, *Sodium / Calcium exchange: its physiological implications*. Physiol. Rev., **79** (1999), 763-854. [ Links ]

[24] R.J. Clarke and D. J. Kane, *Two Gears of Pumping by the Sodium Pump*. Biophys. J., **93** (2007), 4187-4196. [ Links ]

[25] L.Y. Bourguignon, N. Iida, L. Sobrin and G.J. Bourguignon, *Identification of an IP _{3} receptor in endothelial cells*. J. Cell. Physiol.,

**159**(1994), 29-34. [ Links ]

[26] A. Miyakawa-Naito, P. Uhlen, M. Lal, O. Aizman, K. Mikoshiba, H. Brismar, S. Zelenin and A. Aperia, *Cell Signaling Microdomain with Na,K-ATPase and Inositol 1,4,5-Trisphosphate Receptor Generates Calcium Oscillations*. J. Bio. Chem., **278** (2003), 50355-50361.

[27] V.V. Matchkov, H. Gustafsson, A. Rahman, D.M. Boedtkjer, S. Gorintin, A.K. Hansen, E.V. Bouzinova, H.A. Praetorius, C. Aalkjaer and H. Nilsson, *Interaction Between Na+/K+-Pump and Na+/Ca ^{2+}-Exchanger Modulates Intercellular Communication*. Circ. Res.,

**100**(2007), 1026-1035. [ Links ]

[28] G.W. De Young and J. Keizer, *A single-pool inositol 1,4,5-trisphosphate-receptor-based model for agonist-stimulated oscillations in Ca ^{2+} concentration.* Proc.Natl. Acad. Sci. USA,

**89**(1992), 9895-9899. [ Links ]

[29] J. Sneyd, S. Girard and D. Clapham, *Calcium wave propagation by calcium-induced calcium release: an unusual excitable system.* Bull. Math. Biol., **55** (1993), 315-344. [ Links ]

#CAM-365/11.

Received: 02/V/11.

Accepted: 11/IX/11.

*Present address: Biotechnology & Bioengineering Center, and Department of Physiology, Medical College of Wisconsin, 8701 Watertown Plank Road, Milwaukee, WI 53226, USA.