|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Themed Issue Papers |
1 National University of Singapore, Singapore2 Bioengineering Institute4 Department of Engineering Science, The University of Auckland, New Zealand3 Department of Physiology and Cell Biology, University of Nevada School of Medicine, Reno, NV 89557, USA
| Abstract |
|---|
|
|
|---|
(Received 30 October 2005;
accepted after revision 6 January 2006; first published online 11 January 2006)
Corresponding author M. L. Buist: Division of Bioengineering, National University of Singapore 117576. Email: biebml{at}nus.edu.sg
| Introduction |
|---|
|
|
|---|
In the normal situation, slow waves originate from dominant pacemakers along the greater curvature in the corpus and spread, aborally, through the antrum to the pyloric sphincter. The pacemaker region in the corpus is dominant in the stomach because it generates slow waves at the highest frequency. ICC-MY and pacemaker activity are present in a continuous network from corpus to pylorus, but the pacemaker activity in more distal ICC-MY is slower than in the corpus, allowing the corpus to dominate. In humans, the rate of slow waves in the corpus pacemaker is approximately 3 cycles min1, and so this is the dominant slow wave frequency in the entire stomach. Disruptions in the network of ICC-MY or breakdown in the frequency gradient along the ICC-MY network from corpus to pylorus can lead to either physical or functional uncoupling of pacemaker activity and block the proximal to distal spread of gastric peristalsis. This can interrupt normal processing of the gastric contents and delay gastric emptying. Arrhythmias have been associated with a variety of gastric motility disorders, including gastroparesis (Smith et al. 2003), gastric myoelectrical dysrhythmia (Qian et al. 2003), atrophy and hypertrophy (Bortoff & Sillin, 1986) and diabetic gastropathy (Koch, 2001).
The electrical activity of the stomach can be recorded non-invasively with cutaneous electrodes, giving rise to what is termed the electrogastrogram (EGG). This was first recorded by Alvarez (1922) from a little old woman whose abdominal wall was so thin that her gastric peristalsis was easily visible (Smallwood & Brown, 1983). It was independently discovered again by Davis et al. (1957), but it was not until 1975 that the gastric origin of the EGG was conclusively demonstrated (Brown et al. 1975). In the last 25 years much has been learnt about the electrical activity of the stomach, and research into the relationship between gastric electrical activity and the EGG has recently increased substantially.
The cutaneous EGG provides an indirect representation of the electrical activity occurring within the abdomen at a small number of recording electrodes. Unlike the standard 12-lead electrocardiogram (ECG) used in electrocardiology, to date no standard has been universally adopted for the number and placement of the EGG electrodes. One system assumes that the lesser curvature of the stomach begins at the xiphoid process and ends at the point where the mid-clavicular line meets the costal margin (Mintchev et al. 1993). Using this approach, the most proximal electrode was placed 5 cm to the left of the xiphoid process on the costal margin, and a row of four further electrodes (3 cm apart) were placed linearly between the first electrode and the junction of the mid-clavicular line with the right costal margin. An alternative set-up using two electrodes sites one electrode between the umbilicus and the xyphoid process and the second electrode on the left side of the abdomen, one-third of the distance from the ventral mid-line to the left axial mid-line, 1 cm below the lowest rib (Patterson et al. 2001). In yet another approach, Chen et al. (1999) used four recording electrodes. The electrodes were centred on a main electrode located 2 cm above the mid-point between the xiphoid process and the umbilicus. Two more electrodes were located on an upper 45 deg angle, with an additional electrode located 4 cm to the right of the central electrode. Arguably the most commonly used set-up in a clinical setting is that described by Koch & Stern (2004) and Jonderko et al. (2005). With this system, three electrodes are used to record a single-channel bipolar EGG. A reference electrode was placed on the right side of the abdomen at least 3 cm below the right costal margin. One recording electrode was placed 23 cm below the rib cage on the left of the abdomen and the other on the mid-line equidistant from the xiphoid process and the umbilicus, thus providing a linear arrangement aligned on an axial plane. Owing to its clinical adoption, this electrode set-up will be adopted for our simulations.
Before one can use the EGG to infer the electrical health of the stomach, it is important to understand what is recorded at an EGG electrode, so a brief explanation of this is given here. At the level of a single cell, charged ions move selectively across the outer cell membrane in response to a stimulus (smooth muscle cells) or during pacemaking activity (ICC), changing the potential difference across the cell membrane and generating an active cellular response. At the tissue level, when a cell is activated but its neighbour is not, a potential difference is set up between the adjacent cells. The result of this is a local non-zero current density whose magnitude is determined by the size of the potential difference and the resistivity of the pathway between the cells. To a good approximation, away from this local current density fluctuation the torso acts as a resistive network, allowing the current density profile to be reflected instantaneously as a potential field on the skin surface. Within the torso, the principle of superposition applies, so the electrical potential recorded by an electrode on the skin surface is in fact the sum of all the local cellular current density changes within the stomach. It should be noted that this is not a direct measurement; currents originating from different locations will follow different paths and flow through different tissues on the way to the surface electrode. The contribution of each cell to the final signal is therefore biased by both the distance to the recording site and the resistivity of the current pathway.
In addition to the desired gastric electrical signals, the cutaneous electrodes will record information arising from other sources of electrical activity (e.g. the small intestine, colon and heart, and activity due to respiration and motion). In an attempt to isolate the electrical activity of gastric origin, the frequency components in the recording that are known not to correspond to gastric activity are removed. Removing frequencies below 1 cycle min1 (0.016 Hz) and above 15 cycles min1 (0.25 Hz) has been recommended, thereby creating a bandpass filter (Koch & Stern, 2004). In practice, the signal recorded in the time domain is transformed to the frequency domain through the Fast Fourier Transform (FFT), the unwanted frequency components are removed and the filtered signal is reconstructed via the inverse FFT.
Typically EGG analysis revolves around an investigation of the frequency dynamics of the recordings. In a normal subject, dominant activity at 0.05 Hz (3 cycles min1) will be observed. With a breakdown in the ICC network, as can occur for instance in diabetes (see Ordog et al. 2000), regions of ICC-MY within the gastric network can become uncoupled from the dominant pacemaker, and the synchronized spread of electrical activity from corpus to pylorus is disrupted. In such a situation, ICC-MY within the antral region may become local pacemakers, and electrical slow waves may be generated in ectopic sites. An increase in the intrinsic frequency of antral pacemakers can lead to functional uncoupling and ectopic pacemaking. This can result in collisions between slow waves propagating from ectopic sites and the normal pacemaker site, disrupting gastric peristalsis and delaying gastric emptying (gastroparesis). A clinical assessment of this condition can lead to various therapeutic techniques to improve gastric emptying and relieve symptoms of gastroparesis. Unfortunately, it may be extremely difficult to evaluate dysrhythmias with a single or even multiple cutaneous sites of recording.
Few modelling studies have been performed with a focus on generating gastric slow waves using anatomically based models. Most mathematical models used to represent the slow wave can be broadly divided into two categories: those that model the stomach as coupled relaxation oscillators and those that attempt to model the underlying physiology. In 1968, Nelsen and Becker suggested that a chain of relaxation oscillators could simulate the electromechanical activity of the small intestine (Nelsen & Becker, 1968). During the early 1970s, Sarna and coworkers further developed this idea, using an array of bi-directionally coupled oscillators to simulate different aspects of gastrointestinal (GI) activity (Sarna et al. 1971, 1972). Although these models are capable of recreating the general behaviour of a GI slow wave, the parameters prescribing this behaviour cannot easily be related to the underlying electrical activity occurring at the cellular and tissue level within the walls of the GI tract. Attempts have recently been made to develop models of GI electrical activity at the cellular level with a stronger biophysical base, e.g. Miftakhov et al. (1999); Aliev et al. (2000). Both of these studies are primarily focused on the small intestine, and at present good cellular models of gastric electrical activity are somewhat lacking.
The geometrical models used in simulations of GI electrical activity within the abdomen have to date been of an idealized nature. These volume conductor models do not include realistic anatomical information and usually model the abdomen as a homogeneous volume conductor (e.g. Bradshaw et al. 2003) and the stomach with a simplified conoid or ellipsoid geometry (Mirizzi et al. 1986; Mintchev & Bowes, 1998; Irimia & Bradshaw, 2005). In what follows, we describe an anatomically and biophysically based model of the human stomach. This model is used to simulate normal gastric slow waves, and these simulations are shown to be in qualitative agreement with recordings of real slow wave activity. This model is then used to address the issue of functional uncoupling and to examine the issue of whether or not such an event can be detected non-invasively with a commonly used EGG technique.
| Methods |
|---|
|
|
|---|
|
|
|
| (1) |
|
| (2) |
The subscripts i and e represent properties of the intracellular and extracellular domains, respectively. The
terms are tissue conductivities (which in general will be tensors), the
terms are potentials, Vm is the transmembrane potential (the potential difference across the cell membrane), Am is the surface-to-volume ratio of the continuum cell membrane, and Cm is the membrane capacitance. It should be noted that after Aliev et al. (2000), the intracellular current flow between the smooth muscle and ICC layers was modelled using a linear potential difference, as opposed to the diffusive coupling specified by the bidomain equations and employed in the remainder of the stomach. The contributions of the local ionic currents from single cells interact with the continuum through the Iion term in the eqn (1). This allows complex cellular dynamics to be incorporated without a consequential increase in the complexity of the tissue level model.
The large derivative continuous elements that describe the stomach geometry as illustrated in Fig. 1B are insufficient for capturing the gastric slow wave. Therefore each of these volume elements is divided into a large number of smaller hexahedra within the local normalized space of each geometric element. The result is a high-resolution structured hexahedral mesh over which the bidomain equations can be successfully solved using the finite element method.
At each vertex in the high-resolution stomach mesh, a description of cellular electrical activity is placed. A modified FitzHughNagumo (FHN) model (FitzHugh, 1961; Nagumo et al. 1962; Aliev et al. 2000) was used to model the electrical activity of the muscle cells and ICCs, since this currently appears to be the most advanced model that explicitly differentiates between the different cell types (ICCs and smooth muscle). The cellular equations are based on a normalized transmembrane potential, u, that varies from 0 to 1, and a single recovery variable, v. These equations are:
|
| (3) |
|
| (4) |
controls the excitability of the system,
is the recovery rate constant and ß is used to shift the cellular equilibrium from an excitatory to an oscillatory state. The values of the parameters described in the original paper of Aliev et al. (2000) have been modified to match published serosal recordings of gastric slow wave (Bauer et al. 1985), and a full description of these parameter values may be found in Buist et al. (2004).
The above equations are sufficient to simulate a gastric slow wave. However, to calculate the corresponding potentials on the torso surface that arise from this slow wave, two further steps were employed. First, the local current density changes prescribed by the (continuum) cellular electrical activity were expressed as equivalent dipole sources. These dipole sources were then used to compute the electrical potential field on the surface of the torso volume conductor. The equivalent dipole sources (J) are computed from the transmembrane potential gradient as shown in eqn (5).
|
| (5) |
i and
e are the intra- and extracellular conductivities and Vm is the transmembrane potential. Each dipole source has a centre and orientation that varies over time to describe the location, direction and strength of the local electrical activity. These dipole sources are then used to compute the electrical fields within the torso surface by solving the Poisson equation:
|
| (6) |
A typical solution then consists of three sets of calculations that must be performed at each time step. First the cellular models are updated throughout the active stomach. Next the transmembrane and extracellular potentials are calculated from the cellular activity and known diffusive tissue properties. The equivalent gastric dipolar sources are then created and the passive torso volume conductor problem is then solved to calculate the electric field throughout the torso. Using this equation set, the resulting electrical activity at electrodes on the body surface can be determined from the (continuum) cellular level activity. The numerical accuracy of this modelling framework has been demonstrated previously in a number of tests (Pullan et al. 2003).
| Results |
|---|
|
|
|---|
The membrane potential at rest and at the peak of the slow wave was set in accordance with the experimental data from Horiguchi et al. (2001) and Bauer et al. (1985). These studies demonstrated that the resting membrane potential is lower in the regions with dense ICC networks. In the smooth muscle layers, the resting membrane potential was set to 60 mV with a peak potential of 38 mV, and in the ICC layers the resting potential was set to 73 mV with a peak potential of 32 mV, thus providing a degree of electrotonic coupling even under resting conditions.
The stomach volume and torso surface geometries derived from the visible human data set (shown in Fig. 1A) were used for these simulations. In order to solve the bidomain equations accurately over the electrically active stomach, a high-resolution structured finite element mesh was embedded within the volume elements of the stomach wall at an average spatial resolution of less than 1 mm. For both the normal and abnormal situations, 300 s of gastric electrical activity was simulated using a time step of 50 ms, giving a total of 6000 time steps. At every time step, one dipole source was calculated within each of the 320 elements used to describe the stomach geometry, meaning that each source represented the net electrical activity from an element. The resulting 320 dipole sources were placed within the boundary element torso surface model (comprising 300 derivative continuous surface elements), and eqn (6) was solved to generate the body surface potential distribution. From these data, monopolar EGG traces were extracted, providing a simulated sampling frequency of 20 Hz. An example of the results obtained from these simulations is shown in Fig. 3, where the results in the upper panel arise from normal slow wave activity and the results in the lower panel show equivalent results when an ectopic antral pacemaker is present.
|
The potential field was recorded at each solution time point from the three locations indicated on the torso surface in Fig. 3. From these data, a single bipolar trace was constructed for the normal and abnormal situations, as is done in the experimental EGG recording procedure. The resulting bipolar EGG traces are displayed in Fig. 4A and B.
|
| Discussion |
|---|
|
|
|---|
In this study we have described our computational model of the gastric slow wave and the resulting abdominal EGG recordings. Both normal and abnormal slow wave activity have been simulated, and bipolar EGG traces have been generated to mimic what is recorded in the clinical situation. An analysis of the frequency dynamics of these EGG traces, as is standard practice with such recordings, has revealed that this metric was not able to distinguish between normal gastric slow wave activity and functional gastric uncoupling caused by a secondary antral pacemaker. Although this computational evidence is presented here for only a single pacing abnormality, we believe that in fact a wide range of conditions exhibiting functional uncoupling would be undetectable using the EGG technique selected for this study. For example, if the primary pacemaker in the proximal corpus was damaged, an antral pacemaker could produce retrograde slow waves that transverse the length of the active stomach. This would produce a single dominant frequency and return a normal diagnosis if the pacing was within the range of frequencies considered clinically normal (between 2.5 and 3.75 cycles min1; Koch & Stern, 2004). The simulations performed here demonstrate that in fact it is not even necessary for the primary pacing site to be disrupted in order to get EGG traces with a single dominant frequency within the normal range but functional electrical uncoupling.
In the cardiac arena, the sensitivity of ECG recordings to geometric factors has been well documented, along with a relative insensitivity to the composition of the passive torso volume (Bradley et al. 2000). At present these remain open questions in relation to the EGG. We are currently moving to address these issues and have begun the task of constructing patient-specific stomach and torso models from medical images for this purpose.
From a clinical perspective, it is very important to determine whether slow waves in the corpus propagate in a normal manner and drive the pacemaker activity of the distal stomach. This requires high-fidelity recording of electrical activity from more than a single site, placement of electrodes in a serial manner from corpus to distal antrum, and a means to analyse records for aberrant propagation. It may be that these requirements exceed the capacity of cutaneous recordings of gastric electrical activity. Focal changes in intrinsic pacemaker frequency of the order that can lead to functional uncoupling would be extremely difficult to discern with present EGG techniques, and the EGG does not appear to be reliable enough to accomplish evaluation of directional propagation. Gastric serosal surface recording, which would provide the most reliable means of evaluation, is not practical because of its invasive nature. We would suggest that multipoint luminal recording, after placement of electrodes via endoscope, might provide the most suitable means of reliable clinical recording of gastric electrical activity. This technique is now being tried in some clinical gastroenterology units.
| References |
|---|
|
|
|---|
Alvarez WC (1922). The electrogastrogram and what it shows. J Am Med Assoc 78, 11161119.
Bauer AJ, Publicover NG & Sanders KM (1985). Origin and spread of slow waves in canine gastric antral circular muscle. Am J Physiol 249, G800G806.[Medline]
Bortoff A & Sillin LF (1986). Changes in intercellular electrical coupling of smooth muscle accompanying atrophy and hypertrophy. Am J Physiol 250, C292C298.[Medline]
Bradley CP, Pullan AJ & Hunter PJ (2000). Effects of material properties and geometry on electrocardiographic forward simulations. Ann Biomed Eng 28, 721741.[Medline]
Bradshaw LA, Ladipo JK, Staton DJ, Wikswo JP Jr & Richards WO (1999). The human vector magnetogastrogram and magnetoenterogram. IEEE Trans Biomed Eng 46, 959970.[CrossRef][Medline]
Bradshaw LA, Myers A, Wikswo JP & Richards WO (2003). A spatio-temporal dipole simulation of gastrointestinal magnetic fields. IEEE Trans Biomed Eng 50, 836847.[CrossRef][Medline]
Brown BH, Smallwood RH, Duthie HL & Stoddard CJ (1975). Intestinal smooth muscle electrical potentials recorded from surface electrodes. Med Biol Eng 13, 97103.[Medline]
Buist ML, Cheng LK, Yassi R, Bradshaw LA, Richards WO & Pullan AJ (2004). An anatomical model of the gastrointestinal system for producing bioelectric and biomagnetic fields. Physiol Meas 25, 849861.[Medline]
Chen JD, Zou X, Lin X, Ouyang S & Liang J (1999). Detection of gastric slow wave propagation from the cutaneous electrogastrogram. Am J Physiol 277, G424G430.[Medline]
Cousins HM, Edwards FR, Hickey H, Hill CE & Hirst GD (2003). Electrical coupling between the myenteric interstitial cells of Cajal and adjacent muscle layers in the guinea-pig gastric antrum. J Physiol 550, 829844.
Davis RC, Garafolo L & Gault FP (1957). An exploration of abdominal potential. J Comp Physiol Psychol 52, 519523.
Dickens EJ, Hirst GD & Tomita T (1999). Identification of rhythmically active cells in guinea-pig stomach. J Physiol 514, 515531.
Familoni BO, Kingma YJ & Bowes KL (1987). Study of transcutaneous and intraluminal measurement of gastric electrical activity in humans. Med Biol Eng Comput 25, 397402.[CrossRef][Medline]
FitzHugh RA (1961). Impulses and physiological states in theoretical models of nerve membrane. Biophys. J. 1, 445466.
Hirst GDS (2001). An additional role for ICC in the control of gastrointestinal motility?J Physiol 537, 1.
Horiguchi K, Semple GS, Sanders KM & Ward SM (2001). Distribution of pacemaker function through the tunica muscularis of the canine gastric antrum. J Physiol 537, 237250.
Irimia A & Bradshaw LA (2005). Ellipsoidal electrogastrographic forward modelling. Phys Med Biol 50, 44294444.[Medline]
Jonderko K, Kasicka-Jonderko A & Blonska-Fajfrowska B (2005). Does body posture affect the parameters of a cutaneous electrogastrogram?J Smooth Muscle Res 41, 133140.[Medline]
Koch KL (2001). Electrogastrography: physiological basis and clinical application in diabetic gastropathy. Diabetes Tech Therap 3, 5162.
Koch KL & Stern RM (2004). Handbook of Electrogastrography. Oxford University Press, New York.
Lines GT, Buist ML, Grottum P, Pullan AJ & Tveito A (2003). Mathematical models and numerical methods for the forward problem in cardiac electrophysiology. Comput Vis Sci 5, 215239.
Miftakhov RN, Abdusheva GR & Christensen J (1999). Numerical simulation of motility patterns of the small bowel. 1. Formulation of a mathematical model. J Theor Biol 197, 89112.[Medline]
Mintchev MP & Bowes KL (1998). Computer simulation of the effect of changing abdominal thickness on the electrogastrogram. Med Eng Phys 20, 177181.[CrossRef][Medline]
Mintchev MP, Bowes KL & Kingma YJ (1993). Accuracy of cutaneous recordings of gastric electrical activity. Gastroenterology 104, 12731280.[Medline]
Mirizzi N, Stella R & Scafoglieri U (1986). Model to simulate the gastric electrical control and response activity on the stomach wall and on the abdominal surface. Med Biol Eng Comput 24, 157163.[Medline]
Nagumo J, Animoto S & Yoshizawa S (1962). An active pulse transmission line simulating nerve axon. Proc Inst Radio Engineers 50, 20612070.
Nelsen TS & Becker JC (1968). Simulation of the electrical and mechanical gradient of the small intestine. Am J Physiol 214, 749757.
Ordog T, Baldo M, Danko R & Sanders KM (2002). Plasticity of electrical pacemaking by interstitial cells of Cajal and gastric dysrhythmias in W/W mutant mice. Gastroenterology 123, 20282040.[CrossRef][Medline]
Ordog T, Takayama I, Cheung WK, Ward SM & Sanders KM (2000). Remodeling of networks of interstitial cells of Cajal in a murine model of diabetic gastroparesis. Diabetes 49, 17311739.[Abstract]
Ordog T, Ward SM & Sanders KM (1999). Interstitial cells of cajal generate electrical slow waves in the murine stomach. J Physiol 518, 257269.
Ozaki H, Stevens RJ, Blondfield DP, Publicover NG & Sanders KM (1991). Simultaneous measurement of membrane potential, cytosolic Ca2+, and tension in intact smooth muscles. Am J Physiol 260, C917C925.[Medline]
Patterson M, Rintala TR, Lloyd TD, Abernethy L, Houghton TD & Williams J (2001). Validation of electrode placement in neonatal electrogastrography. Dig Dis Sci 46, 22452249.[Medline]
Plonsey R & Barr RC (1987). Mathematical modelling of electrical activity of the heart. J Electrocardiol 20, 219226.[Medline]
Pullan AJ, Cheng LK, Yassi R & Buist ML (2003). Modelling gastrointestinal bioelectric activity. Prog Biophys Mol Biol 85, 523550.
Qian LW, Pasricha PJ & Chen JDZ (2003). Origins and patterns of spontaneous and drug-induced canine gastric myoelectrical dysrhythmia. Digest Dis Sci 48, 508515.[Medline]
Sarna SK, Daniel EE & Kingma YJ (1971). Simulation of slow-wave electrical activity of small intestine. Am J Physiol 221, 166175.
Sarna SK, Daniel EE & Kingma YJ (1972). Premature control potentials in the dog stomach and in the gastric computer model. Am J Physiol 222, 15181523.
Smallwood RH & Brown BH (1983). Non-invasive assessment of gastric activity. In Non-invasive Physiological Measurements, vol. 2, ed. Rolfe P, pp. 179217. Academic Press, London.
Smith DS, Williams CS & Ferris CD (2003). Diagnosis and treatment of chronic gastroparesis and chronic intestinal pseudo-obstruction. Gastroenterol Clin North Am 32, 619658.[Medline]
Spitzer V, Ackerman MJ, Scherzinger AL & Whitlock D (1996). The visible human male: a technical report. J Am Med Inform Assoc 3, 118130.
Szurszewski JH (1987). Electrical basis for gastrointestinal motility. In Physiology of the Gastrointestinal Tract, ed. Johnson LR, pp. 383422. Raven, New York.
This article has been cited by other articles:
![]() |
N. P. Smith, E. J. Crampin, S. A. Niederer, J. B. Bassingthwaighte, and D. A. Beard Computational biology of cardiac myocytes: proposed standards for the physiome J. Exp. Biol., May 1, 2007; 210(9): 1576 - 1583. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Hunter Modelling of biological systems themed issue. Exp Physiol, March 1, 2006; 91(2): 283 - 284. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |