Am J Physiol Gastrointest Liver Physiol 294: G989-G995, 2008.
First published February 14, 2008; doi:10.1152/ajpgi.00528.2007
0193-1857/08 $8.00
NEUROREGULATION AND MOTILITY
Quantitative cellular description of gastric slow wave activity
Alberto Corrias and
Martin L. Buist
Division of Bioengineering, National University of Singapore, Singapore
Submitted 13 November 2007
; accepted in final form 11 February 2008
 |
ABSTRACT
|
|---|
Interstitial cells of Cajal (ICC) are responsible for the spontaneous and omnipresent electrical activity in the stomach. A quantitative description of the intracellular processes whose coordinated activity is believed to generate electrical slow waves has been developed and is presented here. In line with recent experimental evidence, the model describes how the interplay between the mitochondria and the endoplasmic reticulum in cycling intracellular Ca2+ provides the primary regulatory signal for the initiation of the slow wave. The major ion channels that have been identified as influencing slow wave activity have been modeled according to data obtained from isolated ICC. The model has been validated by comparing the simulated profile of the slow waves with experimental recordings and shows good correspondence in terms of frequency, amplitude, and shape in both control and pharmacologically altered conditions.
interstitial cells of Cajal; mathematical model; gastric electrophysiology
WITHIN THE GASTROINTESTINAL tract, the stomach generates and propagates an omnipresent rhythmic electrical activity. This regular pattern of depolarizations is often referred to as the gastric slow wave. The interstitial cells of Cajal (ICC) are now widely regarded as the pacemaker cells that generate the slow waves that propagate to neighboring ICCs and smooth muscle (SM) cells.
To date, a consensus has not been reached regarding the cellular mechanisms through which ICC generate slow waves, and several theories exist. Some investigators have identified a Cl– conductance in ICC and suggest that ICC generate spontaneous Ca2+-dependent Cl– currents that initiate slow wave activity (15, 43). The attribution of the spontaneous currents generated by ICC to a Cl– conductance has been challenged by other authors on grounds of the pharmacological agent used to block it and the characteristics of the recorded current itself (36). Nevertheless the involvement of Cl– currents in ICC electrophysiology has also been proposed by other investigators (14), but their presence and function remain unclear. Many experimental observations seem to identify the coordinated activity of mitochondria and endoplasmic reticulum (ER) in Ca2+ handling as a fundamental process for the generation of pacemaker signals (17, 20, 44). Koh et al. (21) proposed that a nonselective cationic conductance (NSCC), regulated by the local intracellular Ca2+ concentration, provides the electrical signal that initiates slow wave activity (36). After the generation of the primary pacemaker signal, a variety of voltage-dependent ion channels, mainly conducting Na+, Ca2+, and K+ ions and embedded in the plasma membrane of the ICC, are activated. The interactions between depolarizing and hyperpolarizing currents result in the observed slow wave profile.
Computational electrophysiological models provide a unique tool to succinctly describe an intricate network of interactions and allow the user to investigate the contribution of each component to the overall observed cellular behavior. Over the past few years, several modeling efforts have been directed toward the gastrointestinal tract at different spatial scales. At the tissue and organ level, electrophysiological descriptions of the stomach (33) and intestine (24) have been developed with the aim of studying slow wave propagation. These models, however, do not employ biophysically based descriptions of the cells whose activity is responsible for the slow waves. Instead, they generally employ a set of polynomial equations that can replicate the shape of a slow wave but do not allow for more in-depth investigations. Similarly, a cable model was developed recently by Edwards and Hirst (8) in which multiple ICC and SM cells were resistively coupled to replicate gastric slow wave activity. Although good correlations were achieved, the electrical descriptions of the two cell types did not attempt to describe the details of the underlying cellular processes.
When compared with other electrically active cell types, relatively few models of single cells from the gastrointestinal tract have appeared in the literature. Few models of gastrointestinal SM cells exist (31, 38), and only the more recent ones (6) are consistent with the current theory that SM cells are not able to produce slow waves themselves but need to be activated by the ICC network. To the best of our knowledge, only one mathematical description of an ICC has been developed to date (46). The authors describe the behavior of an intestinal ICC and carefully characterize most of the intracellular players and ion channels that are believed to participate in the generation and development of the slow waves. Nevertheless, they do not take into account the role played by the mitochondria in the generation of the slow waves despite the fact that several experiments have shown the importance of mitochondrial activity in slow wave generation (17, 44). Moreover, they do not include some of the ion channels whose presence has been confirmed in ICC, such as the Ca2+-activated (BK) K+ channels (5), ether-a-go-go (ERG) K+ channels (30). and Na+ channels (39). No electrophysiological model of a gastric ICC has been developed so far.
We therefore present here a mathematical description of single cell gastric ICC electrophysiology. The model details the main intracellular processes that are currently believed to provide the pacemaker signal as well as the membrane ion channels that have been identified and functionally characterized in ICC.
 |
MATERIALS AND METHODS
|
|---|
Overview of the model.
The model is based on a classical Hodgkin-Huxley approach where the cell membrane is described as an equivalent electrical circuit consisting of a capacitance connected in parallel with a number of variable conductances representing the different pathways for the movement of charged ions. The time dependence of the membrane potential is governed by
 | (1) |
where Vm (in mV) represents the transmembrane potential, Cm (25 pF) is the cell capacitance, and Iion (in pA) represents the sum of the ionic currents crossing the cell membrane.
Figure 1 presents a schematic overview of the model where the four main intracellular compartments [bulk cytoplasm, mitochondria, ER, and submembrane space (SS)] are displayed together with the ion channels and transport mechanisms included in the plasma membrane.

View larger version (38K):
[in this window]
[in a new window]
|
Fig. 1. Schematic view of the model. The intracellular space is divided in four compartments: cytoplasm, mitochondria, endoplasmic reticulum (ER), and submembrane space (SS). All membrane ion channels and transport mechanisms included in the model are depicted. Details of the Pacemaker Unit (PU) are given in Fig. 2.
|
|
The fraction of volume occupied by the mitochondria was estimated at 13% of the total intracellular space for the following reasons: first, it was noted by several authors that ICC display particularly abundant mitochondria (22) and it was therefore deemed reasonable to choose a value higher than estimates given for other cell types [e.g., 5% on pancreatic β-cells (9)]. Second, a value of 13% is compatible with confocal images obtained in the presence of MitoTracker Green FM on cultured murine jejunal ICC when compared with whole cell c-kit staining of the same cell (44).
The proportion of the intracellular volume occupied by the ER was set to 10%, similar to estimates made by other investigators on intestinal ICC (46), pancreatic β-cells (9), and cardiac myocytes (16). The SS is defined to be a subcellular compartment limited by the plasma membrane, mitochondria, and ER. It has been suggested that Ca2+ shuttling between mitochondria and ER within this enclosed space is responsible for the pacemaking capabilities of ICC (36). Quantification of the volume occupied by the SS is made difficult because of the lack of specific morphological information about its geometry. Sanders et al. (36) suggest that its volume should be "very small." Given a lack of quantitative information, a value 0.1% of the cell volume was chosen based on its ability to be able to yield appropriate physiological Ca2+ concentrations. The bulk cytoplasm, defined as the space available for ionic diffusion, is assumed to occupy 70% of the cellular volume. This value is slightly higher than the one used by Fall and Keizer (9; F&K) (50%) because of experimental observations on ICC that have shown a particularly "electron-dense" cytoplasm (22).
Pacemaker unit.
The intracellular mechanisms responsible for the pacemaker activity in gastrointestinal ICC have been the object of controversy. Although many intracellular players have been implicated in the generation of the slow waves, several experimental observations appear to agree on identifying the interplay between mitochondria and ER in Ca2+ handling as the main process involved in initiating the electrical activity in gastrointestinal ICC (36, 42, 44). Based on these observations, comprehensive theories have recently been proposed (36, 20). Although the interested reader may find their detailed description and extensive references in these published reviews, here we briefly summarize pacemaking theory as it applies to the present model (see Fig. 2).

View larger version (35K):
[in this window]
[in a new window]
|
Fig. 2. Schematic view of the PU. The Ca2+ concentration in the SS is the result of the interplay between mitochondrial Ca2+ handling proteins (uniporter and Na+/Ca2+ exchanger), processes involving the ER [uptake through the sarco(endo)plasmic reticulum Ca2+-ATPase (SERCA) pump and inositol trisphosphate (IP3)-mediated release], and a small leakage between the SS and the cytoplasm. Red arrows represents Ca2+ fluxes, whereas blue arrows represent monovalent cationic fluxes (Na+ or K+).
|
|
The initiating event is proposed to be a Ca2+ release from inositol trisphosphate (IP3)-operated stores in the ER into an intracellular subspace (the SS) formed by the ER, mitochondria, and plasma membrane. Upon the subsequent rise in Ca2+ concentration in this subspace, the Ca2+ uniporter on the mitochondrial membrane is gated open, and Ca2+ ions flow into the mitochondria down the steep electrochemical gradient. This is thought to remove a larger number of Ca2+ ions from the subspace than had previously entered from the ER, causing a temporary drop in the SS Ca2+ concentration. A recently identified Ca2+-inhibited NSCC (21) responds to the decrease in SS Ca2+ by opening the channel and allowing a cationic influx that depolarizes the plasma membrane. These depolarizations are referred to as unitary potentials and are thought to be the initiating electrical events that, in turn, activate voltage-gated Ca2+ channels and subsequently all of the other ion channels on the plasma membrane responsible for the observed slow wave profile.
A number of mathematical models describing intracellular Ca2+ dynamics have appeared in literature in the past few decades (29, 34, 37), each focusing on particular aspects of the ubiquitous role of Ca2+ ions in the cytoplasm. After an extensive search, the model developed by Fall and Keizer (9) (F&K) was chosen to represent the behavior of the intracellular Ca2+ handling proteins for the following reasons. First, it describes the activity of mitochondrial uniporter, mitochondrial Na+/Ca2+ exchanger, sarco(endo)plasmic reticulum Ca2+-ATPase pump, and IP3-mediated Ca2+ release from the ER, all of which appear to be included in the pacemaker mechanism of gastric ICC (36). Second, the frequency of the Ca2+ cycling predicted by the F&K model is close to the typical frequency of the generation of the slow wave in the stomach (3 cycles/min) (35). Finally, the F&K model explicitly describes the dynamics of the mitochondrial membrane potential, since this has been shown to play a crucial role in the pacemaker activity of ICC in the guinea pig stomach (44). Briefly, the F&K model describes the behavior of the mitochondrial membrane potential in terms of the main ionic transport proteins on the mitochondrial membrane, including the F0F1 ATPase proton pump, the respiration-driven proton pump, the adenine nucleotide translocator, and the aforementioned Ca2+-handling proteins.
The F&K model is itself largely based on previous descriptions of intracellular dynamics (25, 26, 27), and interested readers may find details and validations of these models in the referenced publications. The kinetics of the F&K model were left unaltered with the exception of the parameters describing the volume fractions of each subspace, which were modified to include the SS.
Membrane ion channels.
Each of the major membrane ionic conductances that have been experimentally described in ICC (shown in Fig. 1) were included in the Iion term of Eq. 1. Details of these conductances are given below.
Calcium channels.
At least two types of Ca2+ channel have been identified in ICC. One is a DHP-resistant conductance that was resolved with patch-clamp experiments on ICC isolated from the murine small intestine (19), whereas the other one has been molecularly classified as an L-type Ca2+ channel (5). We therefore include two equations, one for each Ca2+ channel conductance:
 | (2) |
 | (3) |
where IVDDR and IL-type represent the ionic currents flowing through the dihydropyridine-resistant and L-type channels, respectively. ECa is the Nernst potential for Ca2+ ions, whereas dVDDR and fVDDR represent activation and inactivation gating variables. Their steady-state equations have been taken directly from experimental values, and their time constants have been chosen to replicate the results from voltage-clamp experiments (19) (Table 1).
In the absence of direct measurements of L-type Ca2+ currents from isolated gastric ICC, a description of L-type Ca2+ channels in gastric SM was adopted (6). Again the dL-type and fL-type represent voltage-dependent activation and inactivation gating variables (Table 1) while fCa represents the Ca2+-dependent inactivation.
 | (4) |
where
[Ca]i is the variation of cytoplasmic concentration (in nM) from the value in resting conditions while hCa (201.4 nM) and sCa (13.1 nM) are the half-concentration and slope factor, respectively. The time constant for this inactivation was chosen to be 2 ms (6). The values of the maximal conductances GVDDR and GL-type (3 and 2 nS, respectively) were chosen to reproduce the voltage-clamp experiments of Kim et al. (19).
Potassium channels.
Voltage-dependent K+ channels encoded by the Kv1.1 gene have been identified in cultured murine ICC by immunohistochemical experiments using the double-labeling technique and functionally characterized with patch-clamp experiments on Xenopus ooctyes transfected with mRNA encoding canine Kv1.1 (12). We adopt the following equation
 | (5) |
Here Ikv11 is the current flowing through these voltage-dependent K+ channels. Gkv11 (6.3 nS) is the maximal conductance, EK is the Nernst potential for K+ ions, and dkv11 and fkv11 represent voltage-dependent gating variables whose parameters have been fitted to voltage-clamp experiments (12) (Table 1).
ERG channels have been proposed as important regulators of pacemaker activity in ICC (47). The characteristics of the ERG channels in ICC have been studied in detail in ICC by McKay et al. (30) from which the following equation was derived
 | (6) |
where IERG is the current flowing through ERG channels. dERG is a voltage-dependent gating variable, and GERG (2.5 nS) is the maximal conductance.
The presence of BK channels has been confirmed in ICC via immunocytochemical techniques in the murine jejunum and lower esophageal sphincter (5). A similar conductance had been previously identified in canine colonic ICC (23), and no major difference between BK conductance in ICC and surrounding myocytes was found in terms of their voltage and intracellular Ca2+ dependency (10). The BK channels in ICC were therefore assumed to be not substantially different from those in gastric SM cells, and the description of Corrias and Buist (6) was adopted.
 | (7) |
Here IBK is the ionic current flowing through BK channels, GBK (23 nS) represents the temperature-dependent maximal conductance, and PBK is the Ca2+ and voltage-dependent open probability that is given by
 | (8) |
where the values of KBK (–17), Caset (100 µM), and hBK (2) are taken from Carl et al. (4).
A background K+ conductance was also included in the model and was described as
 | (9) |
where IKb represents the total background potassium current with a maximal conductance GKb of 0.15 nS.
Sodium channels.
An SCN5A-encoded Na+ conductance has been found in human intestinal ICC by RT-PCR experiments (39). This current, believed to influence the rise phase of the slow waves, is tetrodotoxin resistant. The following equation was used to describe this current.
 | (10) |
Here INa represents the ionic current flowing through Na+ channel, ENa is the Nernst potential for Na+ ions, and dNa and fNa represent voltage-dependent gating variables whose parameters (Table 1) along with the maximal conductance GNa (20 nS) have been chosen to replicate voltage-clamp experiments (39).
NSCC.
NSCC are theorized to be the primary pacemaking conductance and allow ionic fluxes into the SS. The following equation was introduced to describe the NSCC channels
 | (11) |
where INSCC is the ionic current, GNSCC is the maximal conductance (12.15 nS), and ENSCC is the reversal potential for NSCC channels that was calculated with the Goldman-Hodgkin-Katz equation with the ratio between Na+ and K+ permeability set to 1.056 in agreement with the notion that NSCC channels are substantially equally permeable to Na+ and K+. The value of GNSCC was chosen to produce pacemaker potentials of sufficient size to activate the conductances on the plasma membrane responsible for the upstroke of the slow wave under normal conditions. PopenNSCC is the Ca2+-dependent open probability whose Ca2+ dependence has been briefly illustrated previously (36). Here a Hill equation was used to describe the experimental data:
 | (12) |
where [Ca]PU is the Ca2+ concentration in the pacemaker unit, Kopen NSCC and hopen NSCC are the half-value and Hill coefficient, respectively, and their values (74.5 µM and –85, respectively) have been chosen to yield a physiological response with respect to the expected Ca2+ concentrations in the SS.
Chloride channels.
The presence and functions of Cl– in ICC has been the object of controversy (36). Zhu et al. (48) identified an inwardly rectifying whole cell current that the authors believe is carried by Cl– ions. This dependence was described by the following equation
 | (13) |
where ICl is the ionic current, ECl is the Nernst potential for Cl– ions, GCl (10.1 nS) is the maximal conductance, and dCl is a gating variable that takes into account the mechanisms of intracellular regulation of Cl– channels. Although not fully understood, some investigators have suggested that intracellular Ca2+ may play a role in regulating Cl– currents in ICC (43, 48). Intracellular Ca2+ is already proposed as a regulatory component of Cl– channels in SM cells in different organs with thresholds of activation >100 nM (3). Here the following Hill equation was used to model the steady-state effect of intracellular Ca2+ on the Cl– current.
 | (14) |
Here KCl-act represents the activation threshold and hCl represents the Hill coefficient. Their values (140 nM and 3, respectively) have been chosen in view of the intracellular Ca2+ concentrations experienced in ICC as well as taking into account activation thresholds measured in other cell types for this type of conductance (4). A time constant of 30 ms has been used to represent the kinetics of Cl– channels opening in the calculation of dCl.
Ca2+ extrusion mechanism.
At least two proteins have been identified in ICC for the purpose of Ca2+ extrusion from the cytoplasm, the plasmalemmal Ca2+ pump and the Na+/Ca2+ exchanger (5). Because quantitative data pertaining to the kinetics of these proteins was not available, a phenomenological description of these Ca2+ extrusion mechanisms was included to provide long-term intracellular homeostasis.
 | (15) |
Here ICa-EXT is the Ca2+ efflux (in mM/s), ICa-EXT-MAX (0.0885 mM/s) and KCa-EXT (298 nM) are the maximal flux and the half-concentration, respectively.
 |
RESULTS
|
|---|
Integration of Eq. 1 with the incorporation of the other model equations yields the ICC membrane potential as a function of time as is shown in Fig. 3.

View larger version (17K):
[in this window]
[in a new window]
|
Fig. 3. A: slow waves recorded from guinea pig antral interstitial cells of Cajal [ICC; adapted from Hirst and Edwards (13)]. B: simulated gastric slow waves obtained after 10 min of simulated slow wave activity from the presented model.
|
|
The generation of slow waves occurs spontaneously and steadily, even for long simulation times (>1 h). The frequency is 3.02 waves/min, in general agreement with experimental recordings at 37 degrees from the canine gastric antrum [2.8 waves/min (44)] and the guinea pig gastric antrum [3.6 waves/min (13)]. Three waves per minute is also the dominant frequency in the normal human stomach. The value of the resting potential is –67.6 mV, again, in agreement with experimental measurements from guinea pig gastric antrum [ranging from –65.7 mV (13) to –70 mV (44)]. The primary component and plateau phase can be clearly distinguished in the profile of a simulated slow wave (Fig. 4). The amplitude of the primary component is 43.63 mV, and the peak amplitude is 43.8 mV. These values were found to have an average of 44.5 and 45.9 mV, respectively, in a series of experiments on guinea pig antral ICC (13). The half-width and maximal upstroke rate of the simulated slow wave are 8.31 s and 0.5 V/s, respectively; both values are very similar to the ones obtained from guinea pig antral ICC [8.36 s and 0.4 V/s, respectively (13)].
To further ensure the validity of the model, the effect of several different experimental conditions was simulated to determine their effect on the behavior of the generated slow waves.
Involvement of mitochondrial activity in generation of slow waves.
The involvement of mitochondria in the pacemaker mechanism in ICC has been the subject of several experiments. An addition of CGP-37157, a potent inhibitor of the mitochondrial Na+/Ca2+ exchanger, caused the inhibition of slow wave activity in the murine intestine (17). This was simulated by altering the conductance of the mitochondrial Na+/Ca2+ exchanger, and under these conditions the model also ceased to generate slow waves. Similarly, experiments have been performed to determine the involvement of the mitochondrial membrane potential on the generation of slow waves in the guinea pig (2, 11) and canine (45) stomach as well as in murine jejunum (44). In all cases, the addition of mitochondrial uncouplers such as CCCP caused a cessation of slow wave generation. Such an addition was simulated by decreasing the magnitude of the mitochondrial membrane potential. Again the model correctly showed a termination of slow waves generation under such conditions.
Involvement of the ER on slow wave activity.
Malysz et al. (28) investigated the influence of IP3-dependent Ca2+ release on slow wave activity in the murine small intestine. They found that blockade of IP3-dependent Ca2+ release by Xestospongin C eventually abolished slow wave activity. Similar observations were made in similar tissues by Ward et al. (44). On the other hand, it was observed that stimulation of IP3-dependent Ca2+ release by generation of IP3 through a G protein-dependent pathway (triggered by stimulation of
1-adrenoreceptors) caused an increase in the frequency of slow wave activity coupled with a suppressed plateau phase in the murine small intestine (28). The model was tested under both of these conditions by conducting two separate simulations, one with the value of the IP3-mediated Ca2+ efflux from the ER set to zero (blockade) and the other with an increased concentration of IP3. Figure 5, B and C, shows the results of simulated slow waves in such conditions compared with control conditions in Fig. 5A.

View larger version (17K):
[in this window]
[in a new window]
|
Fig. 5. A: simulated slow waves in control condition. B: simulated presence of Xestospongin C by setting to zero the Ca2+ release from the ER. C: simulated presence of 1-adrenoreceptor stimulation with consequent increase in IP3 level by 70 nM.
|
|
Here the lack of IP3-mediated Ca2+ release abolished slow waves (Fig. 5B), in agreement with the observed effects of Xestospongin C (28). Increasing the concentration of IP3 by 70 nM caused an increase in slow wave frequency to 4.0 cycles/min and a reduced plateau phase (half-width reduced from 8.2 to 5.8 s; Fig. 5C), in agreement with the experimental observations (28).
2-Aminoethoxydiphenyl borate (2-APB) is known to be a membrane-permeable agent capable of compromising Ca2+ handling by the ER. Recent reports suggest that it directly inhibits Ca2+ entry pathways in the ER (3, 7). Fukuta et al. (11) observed that addition of 2-APB caused a reversible blockade of slow wave activity in the guinea pig stomach (Fig. 6A). We simulated the presence of 2-APB by inhibiting the Ca2+ entry pathway in the ER. The results of these simulations showed a reversible cessation of slow wave activity (Fig. 6B) in agreement with the experimental data.

View larger version (40K):
[in this window]
[in a new window]
|
Fig. 6. A: effect of addition of 2-aminoethoxydiphenyl borate (2-APB) on ICC gastric slow waves from guinea pig adapted from Fukuta et al. (11). B: simulated slow waves in presence of 2-APB, a blocker of Ca2+ entry in ER.
|
|
 |
DISCUSSION
|
|---|
Although the first hints of autonomous ICC electrical activity became apparent since the end of nineteenth century, it was only about two decades ago that ICC were definitively established as pacemaker cells in the gastrointestinal region (1, 23). Since then, ICC have received an ever-increasing amount of attention by many research groups, since several motility dysfunctions such as gastroparesis and gastrointestinal stromal tumors have been related to abnormal ICC behavior (40, 41). In spite of this, the cellular mechanisms that underlie slow wave generation are not yet fully understood. Here the hypothesis that a Ca2+-inhibited NSCC is the primary pacemaker conductance that has been adopted. This operates under the primary influence of the Ca2+ concentration in a space delimited by ER, mitochondria, and plasma membrane, as proposed by Sanders et al. (36). There is also debate regarding the possible molecular candidates for the NSCC. Earlier studies seemed to point at the involvement of the classic transient receptor potential (TRPC4); however, the observation that TRPC4–/– mice displayed normal slow wave activity prompted further investigations (18). Although a relationship between a melastatin-type transient receptor potential (TRPM7) and pacemaker currents in ICC has been suggested (18), it has not been widely accepted, and the molecular identity of the NSCC remains unknown. Although a definitive determination of the molecular identity of the NSCC would allow a better understanding of ICC pathophysiology, this would not have a direct impact on the validity of this model, since the properties of the channel and not its molecular identity are of primary importance here.
There exist other theories that differ from the NSCC hypothesis. Of particular note is the hypothesis that it is a Cl– conductance and not a NSCC that is the primary pacemaker conductance (48). Should the NSCC theory be displaced in the future, this model, because of its modular nature, can easily be updated by replacing the current primary pacemaker conductance (NSCC) with an equation describing the kinetics of the alternate mechanism.
Despite this ongoing controversy about the nature of the primary pacemaker conductance, there seem to be agreement on the role played by Ca2+ and its handling by the ER and mitochondria in regulating ion channels and, consequently, in the generation of slow waves. The model was therefore developed starting with a validated model of mitochondrial and ER activity (9). To this, descriptions of the main membrane ion channels that have been identified in ICC were coupled. Where possible, data were adopted from isolated gastric ICC. Where direct measurements of currents from gastric ICC were not available, data from intestinal ICC or from gastric SMC were used. The equations that were chosen to describe the ion channels were validated by simulated voltage-clamp experiments, and the results were compared with published experimental data under the same conditions wherever possible.
Some debate also exists with regard to the ionic conductances that are involved in the plateau phase. In particular, the form of the depolarizing current that counterbalances the activation of K+ conductances during the plateau has not been resolved. It has been suggested that more than one conductance may participate and that L-type Ca2+ channels may have a role because of their incomplete inactivation during voltage-clamp experiments (32). The same investigators deemed it unlikely that L-type Ca2+ channels were the sole depolarizing factor during the plateau phase and observed that manipulation of Cl– and Na+ concentrations led to modifications of the slow wave plateau. The involvement of Cl– (20) and possibly other Na+-dependent membrane ionic mechanisms was therefore hypothesized. A Cl– conductance was included in the model description because of reports of its presence and influence in ICC activity (15, 48), in particular in the regenerative component of the slow wave (14). Here the Cl– conductance does not constitute a primary pacemaker conductance but contributes to the plateau phase of the slow wave. Other investigators have suggested that the plateau of the slow wave might arise from the summation of many pacemaker potentials entrained successively by an initial event (36). The velocity of slow wave conduction would, however, suggest that in this case either not all of the pacemaker units are activated by the initial upstroke of the slow wave as it moves through the tissue or that some pacemaker units recover quickly enough to be discharged a second time during the plateau phase. As is the case in other cell types, here the SS represents the summation of individual subcellular spaces, but it is clear that in the future further investigation is needed in this area.
The four K+ conductances included in the model represent those that have been well quantified through patch-clamp experiments and are all believed to provide significant contributions to the slow wave. It should be noted, however, that this is not an exhaustive list of those K+ channels that have been found in ICC. In particular, ATP-dependent K+ channels and intermediate and small Ca2+-activated K+ channels may also be included in the model description once comprehensive quantitative characterizations become available.
The model can be thought of as two tightly coupled components, the internal calcium cycling mechanisms that initiate pacemaker activity and the components on the plasma membrane that form the slow wave in response to this activity. As demonstrated in Fig. 5C, spontaneous rhythmicity can be maintained over a range of pacing frequencies by an appropriate manipulation of the appropriate calcium cycling parameters. The shape of the resulting slow wave is primarily determined by the parameters of the plasma membrane components, and, in general, these may be adjusted to alter the shape of the slow wave without a cessation of activity.
In conclusion, a quantitative model of slow wave generation by a gastric ICC has been developed. Included are descriptions of the main ion channels and intracellular processes that are believed to contribute to the generation of the slow waves. Such models represent a situation whereby the underlying biology is adequately represented and the user has complete control over all of the cellular parameters. This control allows predictive investigations to be performed in relation to determining the likely cause of a certain pathology, or the likely consequences of a certain intervention or mutation. An example of this is whether an ion channel (genetic) mutation could result in clinically relevant anomalies. Experimentally, the behavior of a channel is usually studied in isolation, whereas a model can allow a determination of how the other wild-type channels react to a given situation and in addition can predict changes in whole cell behavior.
The development of this model is seen as an initial step, and further efforts are needed to establish computational models as an important tool for investigating gastric pathophysiology. In particular, the inclusion of existing single cell models of ICC and SM cells (6) in multicellular tissue-level models may allow important insights into the spatial propagation of slow waves in health and disease. A sample implementation of the model is available from the authors upon request.
 |
FOOTNOTES
|
|---|
Address for reprint requests and other correspondence: M. L. Buist, Division of Bioengineering, National Univ. of Singapore, 9 Engineering Dr. 1, Singapore 117576 (e-mail: biebml{at}nus.edu.sg)
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
 |
REFERENCES
|
|---|
- Barajas-Lopez C, Berezin I, Daniel EE, Huizinga JD. Pacemaker activity recorded in interstitial cells of Cajal of the gastrointestinal tract. Am J Physiol Cell Physiol 257: C830–C835, 1989.[Abstract/Free Full Text]
- Beckett EA, Bayguinov YR, Sanders KM, Ward SM, Hirst GD. Properties of unitary potentials generated by intramuscular interstitial cells of Cajal in the murine and guinea-pig gastric fundus. J Physiol 559: 259–269, 2004.[Abstract/Free Full Text]
- Bootman MD, Collins TJ, Mackenzie L, Roderick HL, Berridge MJ, Peppiatt CM. 2-Aminoethoxydiphenyl borate (2-APB) is a reliable blocker of store-operated Ca2+ entry but an inconsistent inhibitor of InsP3-induced Ca2+ release. FASEB J 16: 1145–1150, 2002.[Abstract/Free Full Text]
- Carl A, Lee HK, Sanders KM. Regulation of ion channels in smooth muscles by calcium. Am J Physiol Cell Physiol 271: C9–C34, 1996.[Abstract/Free Full Text]
- Cho WJ, Daniel EE. Proteins of interstitial cells of Cajal and intestinal smooth muscle, colocalized with caveolin-1. Am J Physiol Gastrointest Liver Physiol 288: G571–G585, 2005.[Abstract/Free Full Text]
- Corrias A, Buist ML. A quantitative model of gastric smooth muscle cellular activation. Ann Biomed Eng 35: 1595–1607, 2007.[CrossRef][Web of Science][Medline]
- Dobrydneva Y, Blackmore P. 2-Aminoethoxydiphenyl borate directly inhibits store-operated calcium entry channels in human platelets. Mol Pharmacol 60: 541–552, 2001.[Abstract/Free Full Text]
- Edwards FR, Hirst GD. An electrical description of the generation of slow waves in the antrum of the guinea-pig. J Physiol 564: 213–232, 2005.[Abstract/Free Full Text]
- Fall CP, Keizer JE. Mitochondrial modulation of intracellular Ca(2+) signaling. J Theor Biol 210: 151–165, 2001.[CrossRef][Web of Science][Medline]
- Farrugia G. Ionic conductances in gastrointestinal smooth muscles and interstitial cells of Cajal. Annu Rev Physiol 61: 45–84, 1999.[CrossRef][Web of Science][Medline]
- Fukuta H, Kito Y, Suzuki H. Spontaneous electrical activity and associated changes in calcium concentration in guinea-pig gastric smooth muscle. J Physiol 540: 249–260, 2002.[Abstract/Free Full Text]
- Hatton WJ, Mason HS, Carl A, Doherty P, Latten MJ, Kenyon JL, Sanders KM, Horowitz B. Functional and molecular expression of a voltage-dependent K(+) channel (Kv1.1) in interstitial cells of Cajal. J Physiol 533: 315–327, 2001.[Abstract/Free Full Text]
- Hirst GD, Edwards FR. Generation of slow waves in the antral region of guinea-pig stomach-a stochastic process. J Physiol 535: 165–180, 2001.[Abstract/Free Full Text]
- Hirst GD, Bramich NJ, Teramoto N, Suzuki H, Edwards FR. Regenerative component of slow waves in the guinea-pig gastric antrum involves a delayed increase in [Ca(2+)](i) and Cl(–) channels. J Physiol 540: 907–919, 2002.[Abstract/Free Full Text]
- Huizinga JD, Zhu Y, Ye J, Molleman A. High-conductance chloride channels generate pacemaker currents in interstitial cells of Cajal. Gastroenterology 123: 1627–1636, 2002.[CrossRef][Web of Science][Medline]
- Jafri MS, Gillo B. A membrane potential model with counterions for cytosolic calcium oscillations. Cell Calcium 16: 9–19, 1994.[CrossRef][Web of Science][Medline]
- Kim BJ, Jun JY, So I, Kim KW. Involvement of mitochondrial Na+–Ca2+ exchange in intestinal pacemaking activity. World J Gastroenterol 12: 796–799, 2006.[Medline]
- Kim BJ, So I, Kim KW. The relationship of TRP channels to the pacemaker activity of interstitial cells of Cajal in the gastrointestinal tract. J Smooth Muscle Res 42: 1–7, 2006.[CrossRef][Medline]
- Kim YC, Koh SD, Sanders KM. Voltage-dependent inward currents of interstitial cells of Cajal from murine colon and small intestine. J Physiol 541: 797–810, 2002.[Abstract/Free Full Text]
- Kito Y, Suzuki H. Electrophysiological properties of gastric pacemaker potentials. J Smooth Muscle Res 39: 163–173, 2003.[CrossRef][Medline]
- Koh SD, Jun JY, Kim TW, Sanders KM. A Ca(2+)-inhibited non-selective cation conductance contributes to pacemaker currents in mouse interstitial cell of Cajal. J Physiol 540: 803–814, 2002.[Abstract/Free Full Text]
- Komuro T. Comparative morphology of interstitial cells of Cajal: ultrastructural characterization. Microsc Res Tech 47: 267–285, 1999.[CrossRef][Web of Science][Medline]
- Langton P, Ward SM, Carl A, Norell MA, Sanders KM. Spontaneous electrical activity of interstitial cells of Cajal isolated from canine proximal colon. Proc Natl Acad Sci USA 86: 7280–7284, 1989.[Abstract/Free Full Text]
- Lin AS, Buist ML, Smith NP, Pullan AJ. Modelling slow wave activity in the small intestine. J Theor Biol 242: 356–362, 2006.[CrossRef][Web of Science][Medline]
- Magnus G, Keizer J. Minimal model of beta-cell mitochondrial Ca2+ handling. Am J Physiol Cell Physiol 273: C717–C733, 1997.[Abstract/Free Full Text]
- Magnus G, Keizer J. Model of beta-cell mitochondrial calcium handling and electrical activity. I. Cytoplasmic variables. Am J Physiol Cell Physiol 274: C1158–C1173, 1998.[Abstract/Free Full Text]
- Magnus G, Keizer J. Model of beta-cell mitochondrial calcium handling and electrical activity. II. Mitochondrial variables. Am J Physiol Cell Physiol 274: C1174–C1184, 1998.[Abstract/Free Full Text]
- Malysz J, Donnelly G, Huizinga JD. Regulation of slow wave frequency by IP3-sensitive calcium release in the murine small intestine. Am J Physiol Gastrointest Liver Physiol 280: G439–G448, 2001.[Abstract/Free Full Text]
- Marhl M, Haberichter T, Brumen M, Heinrich R. Complex calcium oscillations and the role of mitochondria and cytosolic proteins. Biosystems 57: 75–86, 2000.[CrossRef][Web of Science][Medline]
- McKay CM, Ye J, Huizinga JD. Characterization of depolarization-evoked ERG K currents in interstitial cells of Cajal. Neurogastroenterol Motil 18: 324–333, 2006.[CrossRef][Web of Science][Medline]
- Miftakhov RN, Abdusheva GR, Christensen J. Numerical simulation of motility patterns of the small bowel. I. Formulation of a mathematical model. J Theor Biol 197: 89–112, 1999.[CrossRef][Web of Science][Medline]
- Molleman J. Ion Channels involved in gastrointestinal action potential generation. In: Pacemaker Activity and Intracellular Communication. Boca Raton, FL: CRC, 1995, p. 223–235.
- Pullan A, Cheng L, Yassi R, Buist M. Modelling gastrointestinal bioelectric activity. Prog Biophys Mol Biol 85: 523–550, 2004.[CrossRef][Web of Science][Medline]
- Roux E, Noble PJ, Noble D, Marhl M. Modelling of calcium handling in airway myocytes. Prog Biophys Mol Biol 90: 64–87, 2006.[CrossRef][Web of Science][Medline]
- Sanders KM, Publicover N. Electrophysiology of the gastric musculature. In: Handbook of Physiology. New York: Academic, 1995, p. 187–216.
- Sanders KM, Koh SD, Ward SM. Interstitial cells of cajal as pacemakers in the gastrointestinal tract. Annu Rev Physiol 68: 307–343, 2006.[CrossRef][Web of Science][Medline]
- Schuster S, Marhl M, Hofer T. Modelling of simple and complex calcium oscillations. From single-cell responses to intercellular signalling. Eur J Biochem 269: 1333–1355, 2002.[Web of Science][Medline]
- Skinner FK, Ward CA, Bardakjian BL. Pump and exchanger mechanisms in a model of smooth muscle. Biophys Chem 45: 253–272, 1993.[CrossRef][Web of Science][Medline]
- Strege PR, Ou Y, Sha L, Rich A, Gibbons SJ, Szurszewski JH, Sarr MG, Farrugia G. Sodium current in human intestinal interstitial cells of Cajal. Am J Physiol Gastrointest Liver Physiol 285: G1111–G1121, 2003.[Abstract/Free Full Text]
- Streutker CJ, Huizinga JD, Driman DK, Riddell RH. Interstitial cells of Cajal in health and disease. Part I: normal ICC structure and function with associated motility disorders. Histopathology 50: 176–189, 2007.[CrossRef][Web of Science][Medline]
- Streutker CJ, Huizinga JD, Driman DK, Riddell RH. Interstitial cells of Cajal in health and disease. Part II: ICC and gastrointestinal stromal tumours. Histopathology 50: 190–202, 2007.[CrossRef][Web of Science][Medline]
- Takaki M. Gut pacemaker cells: the interstitial cells of Cajal (ICC). J Smooth Muscle Res 39: 137–161, 2003.[CrossRef][Medline]
- Tokutomi N, Maeda H, Tokutomi Y, Sato D, Sugita M, Nishikawa S, Nishikawa S, Nakao J, Imamura T, Nishi K. Rhythmic Cl- current and physiological roles of the intestinal c-kit-positive cells. Pflugers Arch 431: 169–177, 1995.[CrossRef][Web of Science][Medline]
- Ward SM, Ordog T, Koh SD, Baker SA, Jun JY, Amberg G, Monaghan K, Sanders KM. Pacemaking in interstitial cells of Cajal depends upon calcium handling by endoplasmic reticulum and mitochondria. J Physiol 525: 355–361, 2000.[Abstract/Free Full Text]
- Ward SM, Dixon RE, Faoite A, Sanders KM. Voltage-dependent calcium entry underlies propagation of slow waves in canine gastric antrum. J Physiol 561: 793–810, 2004.[Abstract/Free Full Text]
- Youm JB, Kim N, Han J, Kim E, Joo H, Leem CH, Goto G, Noma A, Earm YE. A mathematical model of pacemaker activity recorded from mouse small intestine. Philos Transact A Math Phys Eng Sci 364: 1135–1154, 2006.[CrossRef][Medline]
- Zhu Y, Golden CM, Ye J, Wang XY, Akbarali HI, Huizinga JD. ERG K+ currents regulate pacemaker activity in ICC. Am J Physiol Gastrointest Liver Physiol 285: G1249–G1258, 2003.[Abstract/Free Full Text]
- Zhu Y, Mucci A, Huizinga JD. Inwardly rectifying chloride channel activity in intestinal pacemaker cells. Am J Physiol Gastrointest Liver Physiol 288: G809–G821, 2005.[Abstract/Free Full Text]
This article has been cited by other articles:

|
 |

|
 |
 
D. P. Nickerson and M. L. Buist
A physiome standards-based model publication paradigm
Phil Trans R Soc A,
May 28, 2009;
367(1895):
1823 - 1844.
[Abstract]
[Full Text]
[PDF]
|
 |
|
Copyright © 2008 by the American Physiological Society.