Modelling Cl- homeostasis and volume regulation of the cardiac cell

Gap junctions amplify spatial variations in cell volume in proliferating tumor spheroids

Sustained proliferation is a significant driver of cancer progression. Cell-cycle advancement is coupled with cell size, but it remains unclear how multiple cells interact to control their volume in 3D clusters. In this study, we propose a mechano-osmotic model to investigate the evolution of volume dynamics within multicellular systems. Volume control depends on an interplay between multiple cellular constituents, including gap junctions, mechanosensitive ion channels, energy-consuming ion pumps, and the actomyosin cortex, that coordinate to manipulate cellular osmolarity. In connected cells, we show that mechanical loading leads to the emergence of osmotic pressure gradients between cells with consequent increases in cellular ion concentrations driving swelling. We identify how gap junctions can amplify spatial variations in cell volume within multicellular spheroids and, further, describe how the process depends on proliferation-induced solid stress. Our model may provide new insight into the role of gap junctions in breast cancer progression.

Balance of Na+, K+, and Cl- Unidirectional Fluxes in Normal and Apoptotic U937 Cells Computed With All Main Types of Cotransporters

Fluxes of monovalent ions through the multiple pathways of the plasma membrane are highly interdependent, and their assessment by direct measurement is difficult or even impossible. Computation of the entire flux balance helps to identify partial flows and study the functional expression of individual transporters. Our previous computation of unidirectional fluxes in real cells ignored the ubiquitous cotransporters NKCC and KCC. Here, we present an analysis of the entire balance of unidirectional Na⁺, K⁺, and Cl^- fluxes through the plasma membrane in human lymphoid U937 cells, taking into account not only the Na/K pump and electroconductive channels but all major types of cotransporters NC, NKCC, and KCC. Our calculations use flux equations based on the fundamental principles of macroscopic electroneutrality of the system, water balance, and the generally accepted thermodynamic dependence of ion fluxes on the driving force, and they do not depend on hypotheses about the molecular structure of the channel and transporters. A complete list of the major inward and outward Na⁺, K⁺, and Cl^- fluxes is obtained for human lymphoid U937 cells at rest and during changes in the ion and water balance for the first 4 h of staurosporine-induced apoptosis. It is shown how the problem of the inevitable multiplicity of solutions to the flux equations, which arises with an increase in the number of ion pathways, can be solved in real cases by analyzing the ratio of ouabain-sensitive and ouabain-resistant parts of K⁺ (Rb⁺) influx (OSOR) and using additional experimental data on the effects of specific inhibitors. It is found that dynamics of changes in the membrane channels and transporters underlying apoptotic changes in the content of ions and water in cells, calculated without taking into account the KCC and NKCC cotransporters, differs only in details from that calculated for cells with KCC and NKCC. The developed approach to the assessment of unidirectional fluxes may be useful for understanding functional expression of ion channels and transporters in other cells under various conditions. Attached software allows reproduction of all calculated data under presented conditions and to study the effects of the condition variation.

A Tool for Computation of Changes in Na+, K+, Cl- Channels and Transporters Due to Apoptosis by Data on Cell Ion and Water Content Alteration

Monovalent ions are involved in a vast array of cellular processes. Their movement across the cell membrane is regulated by numerous channels and transporters. Identification of the pathways responsible for redistribution of ions and cell water in living cells is hampered by their strong interdependence. This difficulty can be overcome by computational analysis of the whole cell flux balance. Our previous computational studies were concerned with monovalent ion fluxes in cells under the conditions of balanced ion distribution or during transition processes after stopping the Na+/K+ pump. Here we analyze a more complex case-redistribution of ions during cell apoptosis when the parameters keep changing during the process. New experimental data for staurosporine-induced apoptosis of human lymphoma cells U937 have been obtained: the time course of changes in cellular K+, Na+, Cl-, and water content, as well as Rb+ fluxes as a marker of the Na/K pump activity. Using a newly developed computational tool, we found that alteration of ion and water balance was associated with a 55% decrease in the Na+/K+-ATPase rate coefficient over a 4-h period, with a time-dependent increase in potassium channel permeability, and a decrease in sodium channel permeability. The early decrease in [Cl-]i and cell volume were associated with an ~5-fold increase in chloride channel permeability. The developed approach and the presented executable file can be used to identify the channels and transporters responsible for alterations of cell ion and water balance not only during apoptosis but in other physiological scenarios.

The logic of ionic homeostasis: Cations are for voltage, but not for volume

Neuronal activity is associated with transmembrane ionic redistribution, which can lead to an osmotic imbalance. Accordingly, activity-dependent changes of the membrane potential are sometimes accompanied by changes in intracellular and/or extracellular volume. Experimental data that include distributions of ions and volume during neuronal activity are rare and rather inconsistent partly due to the technical difficulty of performing such measurements. However, progress in understanding the interrelations among ions, voltage and volume has been achieved recently by computational modelling, particularly “charge-difference” modelling. In this work a charge-difference computational model was used for further understanding of the specific roles for cations and anions. Our simulations show that without anion conductances the transmembrane movements of cations are always osmotically balanced, regardless of the stoichiometry of the pump or the ratio of Na⁺ and K⁺ conductances. Yet any changes in cation conductance or pump activity are associated with changes of the membrane potential, even when a hypothetically electroneutral pump is used in calculations and K⁺ and Na⁺ conductances are equal. On the other hand, when a Cl^- conductance is present, the only way to keep the Cl^-equilibrium potential in accordance with the changed membrane potential is to adjust cell volume. Importantly, this voltage-evoked Cl^--dependent volume change does not affect intracellular cation concentrations or the amount of energy that is necessary to support the system. Taking other factors into consideration (i.e. the presence of internal impermeant poly-anions, the activity of cation-Cl^- cotransporters, and the buildup of intra- and extracellular osmolytes, both charged and electroneutral) adds complexity, but does not change the main principles.

We have developed software that calculates membrane potential and cell volume that result from redistribution of principal ions (K⁺, Na⁺, and Cl^-) during normal cellular activity and experimental manipulations. Calculations in the model are done by an iterative charge-difference method that makes few assumptions about governing equations. Most of the features that were considered to be important for volume and voltage regulation were incorporated in the model, including the unique capability to perform calculations with different values of transmembrane water permeability. We have used the program to reexamine interactions between ionic fluxes, membrane potential, and cell volume and found that there was a previously unappreciated difference in the way that the distribution of cations and anions affect the cell. Na⁺ and K⁺, which are distributed unevenly across the membrane by the Na⁺/K⁺-ATPase, are primarily responsible for the membrane potential, but, contrary to popular belief, do not directly participate in volume regulation. On the other hand, the Cl^- conductance determines the extent of volume changes, because Cl^- has to follow the changes of membrane potential, which inevitably leads to changes in cell volume. The software is available to download and use for other investigations.

Quantitative Model for Ion Transport and Cytoplasm Conductivity of Chinese Hamster Ovary Cells

In mammalian cells cytoplasm ion concentrations and hence cytoplasm conductivity is an important indicator of their physiological state. Changes in the cytoplasm conductivity has been associated with physiological changes such as progression of cancer and apoptosis. In this work, a model that predicts the effects of physiological changes in ion transport on the cytoplasm conductivity of Chinese hamster ovary (CHO) cells is demonstrated. We determined CHO-specific model parameters, Na⁺/K⁺ ATPase pumps and ion channels densities, using a flux assay approach. The obtained sodium (P_Na), potassium (P_K) and chloride (P_Cl) permeability and Na⁺/K⁺ ATPase pump density were estimated to be 5.6 × 10^-8 cm/s, 5.6 × 10^-8 cm/s, 3.2 × 10^-7 cm/s and 2.56 × 10^-11 mol/cm², respectively. The model was tested by comparing the model predictions with the experimentally determined temporal changes in the cytoplasm conductivity of Na⁺/K⁺ ATPase pump inhibited CHO cells. Cells' Na⁺/K⁺ ATPase pumps were inhibited using 5 mM Ouabain and the temporal behavior of their cytoplasm conductivity was measured using dielectrophoresis cytometry. The measured results are in close agreement with the model-calculated values. This model will provide insight on the effects of processes such as apoptosis or external media ion concentration on the cytoplasm conductivity of mammalian cells.

Mathematical model of the post-ablation enhancement zone as a tissue-level oedematic response

PURPOSE

A hyperdense rim is commonly observed at the periphery of ablation zones during post-ablation imaging (e.g. ultrasound) in tumours. A mathematical model has been developed here to investigate the occurrence of this enhanced rim, caused by the ablated cells, giving an indication of the location of the final ablation region.

MATERIALS AND METHODS

The enhanced rim has been assumed here to be due to a tissue-level oedematic response of viable cells, which necessitated coupling multiple modelling elements in a spatially distributed system: thermal cell death, tissue-state dependent ion concentration dynamics, ion transport in the extracellular space, and osmotic cell volume regulation.

RESULTS

In response to the imposed temperature function, an ablation zone was predicted, distinguishing the tissue state between 'dead' and 'alive'. A disturbance in intracellular/extracellular ion concentrations was induced due to ion redistribution, which acted as an osmotic stress and contributed to significant cell swelling in a thin rim at the periphery of the ablation zone. It was also found that the rim size only changed slightly with varying lesion size, in response to different temperature profiles.

CONCLUSIONS

The study presents a novel mathematical model to understand the enhanced rim surrounding the ablation zone by assuming tissue-level cell oedema as the primary potential cause. The model links the direct response to thermal injury to an observable secondary response, which could be of clinical value in that the location of this bright ring could potentially be used for more accurate determination of the extent of the ablation zone.

Unidirectional Flux Balance of Monovalent Ions in Cells with Na/Na and Li/Na Exchange: Experimental and Computational Studies on Lymphoid U937 Cells

Monovalent ion traffic across the cell membrane occurs via various pathways. Evaluation of individual fluxes in whole cell is hampered by their strong interdependence. This difficulty can be overcome by computational analysis of the whole cell flux balance. However, the previous computational studies disregarded ion movement of the self-exchange type. We have taken this exchange into account. The developed software allows determination of unidirectional fluxes of all monovalent ions via the major pathways both under the balanced state and during transient processes. We show how the problem of finding the rate coefficients can be solved by measurement of monovalent ion concentrations and some of the fluxes. Interdependence of fluxes due to the mandatory conditions of electroneutrality and osmotic balance and due to specific effects can be discriminated, enabling one to identify specific changes in ion transfer machinery under varied conditions. To test the effectiveness of the developed approach we made use of the fact that Li/Na exchange is known to be an analogue of the coupled Na/Na exchange. Thus, we compared the predicted and experimental data obtained on U937 cells under varied Li+ concentrations and following inhibition of the sodium pump with ouabain. We found that the coupled Na/Na exchange in U937 cells comprises a significant portion of the entire Na+ turnover. The data showed that the loading of the sodium pump by Li/Na exchange involved in the secondary active Li+ transport at 1-10 mM external Li+ is small. This result may be extrapolated to similar Li+ and Na+ flux relationships in erythrocytes and other cells in patients treated with Li+ in therapeutic doses. The developed computational approach is applicable for studying various cells and can be useful in education for demonstrating the effects of individual transporters and channels on ion gradients, cell water content and membrane potential.

Mechanisms of astrocytic K(+) clearance and swelling under high extracellular K(+) concentrations

In response to the elevation of extracellular K(+) concentration ([K(+)]out), astrocytes clear excessive K(+) to maintain conditions necessary for neural activity. K(+) clearance in astrocytes occurs via two processes: K(+) uptake and K(+) spatial buffering. High [K(+)]out also induces swelling in astrocytes, leading to edema and cell death in the brain. Despite the importance of astrocytic K(+) clearance and swelling, the underlying mechanisms remain unclear. Here, we report results from a simulation analysis of astrocytic K(+) clearance and swelling. Astrocyte models were constructed by incorporating various mechanisms such as intra/extracellular ion concentrations of Na(+), K(+), and Cl(-), cell volume, and models of Na,K-ATPase, Na-K-Cl cotransporter (NKCC), K-Cl cotransporter, inwardly-rectifying K(+) (KIR) channel, passive Cl(-) current, and aquaporin channel. The simulated response of astrocyte models under the uniform distribution of high [K(+)]out revealed significant contributions of NKCC and Na,K-ATPase to increases of intracellular K(+) and Cl(-) concentrations, and swelling. Moreover, we found that, under the non-uniform distribution of high [K(+)]out, KIR channels localized at synaptic clefts absorbed excess K(+) by depolarizing the equivalent potential of K(+) (E K) above membrane potential, while K(+) released through perivascular KIR channels was enhanced by hyperpolarizing E K and depolarizing membrane potential. Further analysis of simulated drug effects revealed that astrocyte swelling was modulated by blocking each of the ion channels and transporters. Our simulation analysis revealed controversial mechanisms of astrocytic K(+) clearance and swelling resulting from complex interactions among ion channels and transporters.

Computational modeling reveals dendritic origins of GABA(A)-mediated excitation in CA1 pyramidal neurons

GABA is the key inhibitory neurotransmitter in the adult central nervous system, but in some circumstances can lead to a paradoxical excitation that has been causally implicated in diverse pathologies from endocrine stress responses to diseases of excitability including neuropathic pain and temporal lobe epilepsy. We undertook a computational modeling approach to determine plausible ionic mechanisms of GABA(A)-dependent excitation in isolated post-synaptic CA1 hippocampal neurons because it may constitute a trigger for pathological synchronous epileptiform discharge. In particular, the interplay intracellular chloride accumulation via the GABA(A) receptor and extracellular potassium accumulation via the K/Cl co-transporter KCC2 in promoting GABA(A)-mediated excitation is complex. Experimentally it is difficult to determine the ionic mechanisms of depolarizing current since potassium transients are challenging to isolate pharmacologically and much GABA signaling occurs in small, difficult to measure, dendritic compartments. To address this problem and determine plausible ionic mechanisms of GABA(A)-mediated excitation, we built a detailed biophysically realistic model of the CA1 pyramidal neuron that includes processes critical for ion homeostasis. Our results suggest that in dendritic compartments, but not in the somatic compartments, chloride buildup is sufficient to cause dramatic depolarization of the GABA(A) reversal potential and dominating bicarbonate currents that provide a substantial current source to drive whole-cell depolarization. The model simulations predict that extracellular K(+) transients can augment GABA(A)-mediated excitation, but not cause it. Our model also suggests the potential for GABA(A)-mediated excitation to promote network synchrony depending on interneuron synapse location - excitatory positive-feedback can occur when interneurons synapse onto distal dendritic compartments, while interneurons projecting to the perisomatic region will cause inhibition.

A mean-field model of ventricular muscle tissue

A theoretical model of the cross-linking topology of ventricular muscle tissue is developed. Using parameter estimation the terms of the theoretical model are estimated for normal and pathological conditions. The model represents the anisotropic structure of the tissue, reproduces published experimental data and characterizes the role of different tissue components in the observed macroscopic behavior. Changes in the material parameters are consistent with expected structural changes and the model is extended to reproduce force-Calcium relationships. Model results are invoked to argue that semisoft behavior and the material axis anisotropy arise from the constraints on the extracellular matrix cross-linking topology.

Steady-state solutions of cell volume in a cardiac myocyte model elaborated for membrane excitation, ion homeostasis and Ca2+ dynamics

The cell volume continuously changes in response to varying physiological conditions, and mechanisms underlying volume regulation have been investigated in both experimental and theoretical studies. Here, general formulations concerning cell volume change are presented in the context of developing a comprehensive cell model which takes Ca(2+) dynamics into account. Explicit formulas for charge conservation and steady-state volumes of the cytosol and endoplasmic reticulum (ER) are derived in terms of membrane potential, amount of ions, Ca(2+)-bound buffer molecules, and initial cellular conditions. The formulations were applied to a ventricular myocyte model which has plasma-membrane Ca(2+) currents with dynamic gating mechanisms, Ca(2+)-buffering reactions with diffusive and non-diffusive buffer proteins, and Ca(2+) uptake into or release from the sarcoplasmic reticulum (SR) accompanied by compensatory cationic or anionic currents through the SR membrane. Time-dependent volume changes in cardiac myocytes induced by varying extracellular osmolarity or by action potential generation were successfully simulated by the novel formulations. Through application of bifurcation analysis, the existence and uniqueness of steady-state solutions of the cell volume were validated, and contributions of individual ion channels and transporters to the steady-state volume were systematically analyzed. The new formulas are consistent with previous fundamental theory derived from simple models of minimum compositions. The new formulations may be useful for examination of the relationship between cell function and volume change in other cell types.

Regulatory volume decrease and P receptor signaling in fish cells: mechanisms, physiology, and modeling approaches

For animal cell plasma membranes, the permeability of water is much higher than that of ions and other solutes, and exposure to hyposmotic conditions almost invariably causes rapid water influx and cell swelling. In this situation, cells deploy regulatory mechanisms to preserve membrane integrity and avoid lysis. The phenomenon of regulatory volume decrease, the partial or full restoration of cell volume following cell swelling, is well-studied in mammals, with uncountable investigations yielding details on the signaling network and the effector mechanisms involved in the process. In comparison, cells from other vertebrates and from invertebrates received little attention, despite of the fact that e.g. fish cells could present rewarding model systems given the diversity in ecology and lifestyle of this animal group that may be reflected by an equal diversity of physiological adaptive mechanisms, including those related to cell volume regulation. In this review, we therefore present an overview on the most relevant aspects known on hypotonic volume regulation presently known in fish, summarizing transporters and signaling pathways described so far, and then focus on an aspect we have particularly studied over the past years using fish cell models, i.e. the role of extracellular nucleotides in mediating cell volume recovery of swollen cells. We, furthermore, present diverse modeling approaches developed on the basis of data derived from studies with fish and other models and discuss their potential use for gaining insight into the theoretical framework of volume regulation.

A physiome standards-based model publication paradigm

In this era of widespread broadband Internet penetration and powerful Web browsers on most desktops, a shift in the publication paradigm for physiome-style models is envisaged. No longer will model authors simply submit an essentially textural description of the development and behaviour of their model. Rather, they will submit a complete working implementation of the model encoded and annotated according to the various standards adopted by the physiome project, accompanied by a traditional human-readable summary of the key scientific goals and outcomes of the work. While the final published, peer-reviewed article will look little different to the reader, in this new paradigm, both reviewers and readers will be able to interact with, use and extend the models in ways that are not currently possible. Here, we review recent developments that are laying the foundations for this new model publication paradigm. Initial developments have focused on the publication of mathematical models of cellular electrophysiology, using technology based on a CellML- or Systems Biology Markup Language (SBML)-encoded implementation of the mathematical models. Here, we review the current state of the art and what needs to be done before such a model publication becomes commonplace.

Measuring and modeling chloride-hydroxyl exchange in the Guinea-pig ventricular myocyte

Protons are powerful modulators of cardiac function. Their intracellular concentration is regulated by sarcolemmal ion transporters that export or import H+-ions (or their ionic equivalent: HCO3-, OH-). One such transporter, which imports H+-equivalents, is a putative Cl-/OH- exchanger (CHE). A strong candidate for CHE is SLC26A6 protein, a product of the SLC26A gene family of anion transporters, which has been detected in murine heart. SLC26A6 protein is suggested to be an electrogenic 1Cl-/2OH-(2HCO3-) exchanger. Unfortunately, there is insufficient characterization of cardiac CHE against which the properties of heterologously expressed SLC26A6 can be matched. We therefore investigated the proton, Cl-, and voltage dependence of CHE activity in guinea-pig ventricular myocytes, using voltage-clamp, intracellular pH fluorescence, and mathematical modeling techniques. We find that CHE activity is tightly regulated by intracellular and extracellular pH, is voltage-insensitive over a wide range (+/-80 mV), and displays substrate dependence suggestive of electroneutral 1Cl-/1OH- exchange. These properties exclude electrogenic SLC26A6 as sole contributor to CHE. Either the SLC26A6 product in heart is electroneutral, or CHE comprises at least two transporters with oppositely balanced voltage sensitivity. Alternatively, CHE may comprise an H+-Cl- coinflux system, which cannot be distinguished kinetically from an exchanger. Irrespective of ionic mechanism, CHE's pH sensitivity helps to define resting intracellular pH, and hence basal function in the heart.

Role of Ca2+ Transporters and Channels in the Cardiac Cell Volume Regulation

Na+/K+ pump is one of key mechanisms to maintain cell volume. When it is inhibited, cells are at risk of swelling. However, in guinea pig ventricular myocytes, the cell area as an index of cell volume was almost constant during 90 min Na+/K+ pump blockade with 40 microM ouabain despite the marked membrane depolarization. In this study, involvements of Ca2+ transporters and channels in the cardiac cell volume regulation were proposed by conducting the computer simulation in parallel with the experimental validation.

Ionic mechanisms of cardiac cell swelling induced by blocking Na+/K+ pump as revealed by experiments and simulation

Although the Na(+)/K(+) pump is one of the key mechanisms responsible for maintaining cell volume, we have observed experimentally that cell volume remained almost constant during 90 min exposure of guinea pig ventricular myocytes to ouabain. Simulation of this finding using a comprehensive cardiac cell model (Kyoto model incorporating Cl(-) and water fluxes) predicted roles for the plasma membrane Ca(2+)-ATPase (PMCA) and Na(+)/Ca(2+) exchanger, in addition to low membrane permeabilities for Na(+) and Cl(-), in maintaining cell volume. PMCA might help maintain the [Ca(2+)] gradient across the membrane though compromised, and thereby promote reverse Na(+)/Ca(2+) exchange stimulated by the increased [Na(+)](i) as well as the membrane depolarization. Na(+) extrusion via Na(+)/Ca(2+) exchange delayed cell swelling during Na(+)/K(+) pump block. Supporting these model predictions, we observed ventricular cell swelling after blocking Na(+)/Ca(2+) exchange with KB-R7943 or SEA0400 in the presence of ouabain. When Cl(-) conductance via the cystic fibrosis transmembrane conductance regulator (CFTR) was activated with isoproterenol during the ouabain treatment, cells showed an initial shrinkage to 94.2 +/- 0.5%, followed by a marked swelling 52.0 +/- 4.9 min after drug application. Concomitantly with the onset of swelling, a rapid jump of membrane potential was observed. These experimental observations could be reproduced well by the model simulations. Namely, the Cl(-) efflux via CFTR accompanied by a concomitant cation efflux caused the initial volume decrease. Then, the gradual membrane depolarization induced by the Na(+)/K(+) pump block activated the window current of the L-type Ca(2+) current, which increased [Ca(2+)](i). Finally, the activation of Ca(2+)-dependent cation conductance induced the jump of membrane potential, and the rapid accumulation of intracellular Na(+) accompanied by the Cl(-) influx via CFTR, resulting in the cell swelling. The pivotal role of L-type Ca(2+) channels predicted in the simulation was demonstrated in experiments, where blocking Ca(2+) channels resulted in a much delayed cell swelling.

Modelling transcriptional feedback loops: the role of Gro/TLE1 in Hes1 oscillations

The transcriptional repressor Hes1, a basic helix-loop-helix family protein, periodically changes its expression in the presomitic mesoderm. Its periodic pattern of expression is retained in a number of cultured murine cell lines. In this paper, we introduce an extended mathematical model for Hes1 oscillatory expression that includes regulation of Hes1 transcription by Drosophila Groucho (Gro) or its vertebrate counterpart, the transducine-like enhancer of split/Groucho-related gene product 1 (TLE1). Gro/TLE1 is a necessary corepressor required by a number of DNA-binding transcriptional repressors, including Hes1. Models of direct repression via Hes1 typically display an expression overshoot after transcription initiation which is not seen in the experimental data. However, numerical simulation and theoretical predictions of our model show that the cofactor Gro/TLE1 reduces the overshoot and is thus necessary for a rapid and finely tuned response of Hes1 to activation signals. Further, from detailed linear stability and numerical bifurcation analysis and simulations, we conclude that the cooperativity coefficient (h) for Hes1 self-repression should be large (i.e. h>or=4). Finally, we introduce the characteristic turnaround duration, and show that for our model the duration of the repression loop is between 40 and 60 min.

Mathematical models in physiology

Computational modelling of biological processes and systems has witnessed a remarkable development in recent years. The search-term (modelling OR modeling) yields over 58000 entries in PubMed, with more than 34000 since the year 2000: thus, almost two-thirds of papers appeared in the last 5-6 years, compared to only about one-third in the preceding 5-6 decades. The development is fuelled both by the continuously improving tools and techniques available for bio-mathematical modelling and by the increasing demand in quantitative assessment of element inter-relations in complex biological systems. This has given rise to a worldwide public domain effort to build a computational framework that provides a comprehensive theoretical representation of integrated biological function-the Physiome. The current and next issues of this journal are devoted to a small sub-set of this initiative and address biocomputation and modelling in physiology, illustrating the breadth and depth of experimental data-based model development in biological research from sub-cellular events to whole organ simulations.

A mathematical model of pacemaker activity recorded from mouse small intestine

The pacemaker activity of interstitial cells of Cajal (ICCs) has been known to initiate the propagation of slow waves along the whole gastrointestinal tract through spontaneous and repetitive generation of action potentials. We studied the mechanism of the pacemaker activity of ICCs in the mouse small intestine and tested it using a mathematical model. The model includes ion channels, exchanger, pumps and intracellular machinery for Ca2+ regulation. The model also incorporates inositol 1,4,5-triphosphate (IP3) production and IP3-mediated Ca2+ release activities. Most of the parameters were obtained from the literature and were modified to fit the experimental results of ICCs from mouse small intestine. We were then able to compose a mathematical model that simulates the pacemaker activity of ICCs. The model generates pacemaker potentials regularly and repetitively as long as the simulation continues. The frequency was set at 20 min(-1) and the duration at 50% repolarization was 639 ms. The resting and overshoot potentials were -78 and +1.2 mV, respectively. The reconstructed pacemaker potentials closely matched those obtained from animal experiments. The model supports the idea that cyclic changes in [Ca2+]i and [IP3] play key roles in the generation of ICC pacemaker activity in the mouse small intestine.